No Arabic abstract
R-modes in neutron stars with crusts are damped by viscous friction at the crust-core boundary. The magnitude of this damping, evaluated by Bildsten and Ushomirsky (BU) under the assumption of a perfectly rigid crust, sets the maximum spin frequency for a neutron star spun up by accretion in a Low-Mass X-ray binary (LMXB). In this paper we explore the mechanical coupling between the core r-modes and the elastic crust, using a toy model of a constant density neutron star with a constant shear modulus crust. We find that, at spin frequencies in excess of ~50 Hz, the r-modes strongly penetrate the crust. This reduces the relative motion (slippage) between the crust and the core compared to the rigid crust limit. We therefore revise down, by as much as a factor of 10^2-10^3, the damping rate computed by BU, significantly reducing the maximal possible spin frequency of neutron star with a solid crust. The dependence of the crust-core slippage on the spin frequency is complicated, and is very sensitive to the physical thickness of the crust. If the crust is sufficiently thick, the curve of the critical spin frequency for the onset of the r-mode instability becomes multi-valued for some temperatures; this is related to the avoided crossings between the r-mode and the higher-order torsional modes in the crust. The critical frequencies are comparable to the observed spins of neutron stars in LMXBs and millisecond pulsars.
Recent work has raised the exciting possibility that r-modes (Rossby waves) in rotating neutron star cores might be strong gravitational wave sources. We estimate the effect of a solid crust on their viscous damping rate and show that the dissipation rate in the viscous boundary layer between the oscillating fluid and the nearly static crust is >10^5 times higher than that from the shear throughout the interior. This increases the minimum frequency for the onset of the gravitational r-mode instability to at least 500 Hz when the core temperature is less than 10^10 K. It eliminates the conflict of the r-mode instability with the accretion-driven spin-up scenario for millisecond radio pulsars and makes it unlikely that the r-mode instability is active in accreting neutron stars. For newborn neutron stars, the formation of a solid crust shortly after birth affects their gravitational wave spin-down and hence detectability by ground-based interferometric gravitational wave detectors.
It is commonly believed that the magnetic field threading a neutron star provides the ultimate mechanism (on top of fluid viscosity) for enforcing long-term corotation between the slowly spun down solid crust and the liquid core. We show that this argument fails for axisymmetric magnetic fields with closed field lines in the core, the commonly used `twisted torus field being the most prominent example. The failure of such magnetic fields to enforce global crust-core corotation leads to the development of a persistent spin lag between the core region occupied by the closed field lines and the rest of the crust and core. We discuss the repercussions of this spin lag for the evolution of the magnetic field, suggesting that, in order for a neutron star to settle to a stable state of crust-core corotation, the bulk of the toroidal field component should be deposited into the crust soon after the neutron stars birth.
We model the nonlinear saturation of the r-mode instability via three-mode couplings and the effects of the instability on the spin evolution of young neutron stars. We include one mode triplet consisting of the r-mode and two near resonant inertial modes that couple to it. We find that the spectrum of evolutions is more diverse than previously thought. The evolution of the star is dynamic and initially dominated by fast neutrino cooling. Nonlinear effects become important when the r-mode amplitude grows above its first parametric instability threshold. The balance between neutrino cooling and viscous heating plays an important role in the evolution. Depending on the initial r-mode amplitude, and on the strength of the viscosity and of the cooling this balance can occur at different temperatures. If thermal equilibrium occurs on the r-mode stability curve, where gravitational driving equals viscous damping, the evolution may be adequately described by a one-mode model. Otherwise, nonlinear effects are important and lead to various more complicated scenarios. Once thermal balance occurs, the star spins-down oscillating between thermal equilibrium states until the instability is no longer active. For lower viscosity we observe runaway behavior in which the r-mode amplitude passes several parametric instability thresholds. In this case more modes need to be included to model the evolution accurately. In the most optimistic case, we find that gravitational radiation from the r-mode instability in a very young, fast spinning neutron star within about 1 Mpc of Earth may be detectable by advanced LIGO for years, and perhaps decades, after formation. Details regarding the amplitude and duration of the emission depend on the internal dissipation of the modes of the star, which would be probed by such detections.
It has recently been suggested that the tidal deformation of a neutron star excites daughter p- and g-modes to large amplitudes via a quasi-static instability. This would remove energy from the tidal bulge, resulting in dissipation and possibly affecting the phase evolution of inspiralling binary neutron stars and hence the extraction of binary parameters from gravitational wave observations. This instability appears to arise because of a large three-mode interaction among the tidal mode and high-order p- and g-modes of similar radial wavenumber. We show that additional four-mode interactions enter into the analysis at the same order as the three-mode terms previously considered. We compute these four-mode couplings by finding a volume-preserving coordinate transformation that relates the energy of a tidally deformed star to that of a radially perturbed spherical star. Using this method, we relate the four-mode coupling to three-mode couplings and show that there is a near-exact cancellation between the destabilizing effect of the three-mode interactions and the stabilizing effect of the four-mode interaction. We then show that the equilibrium tide is stable against the quasi-static decay into daughter p- and g-modes to leading order. The leading deviation from the quasi-static approximation due to orbital motion of the binary is considered; while it may slightly spoil the near-cancellation, any resulting instability timescale is at least of order the gravitational-wave inspiral time. We conclude that the p-/g-mode coupling does not lead to a quasi-static instability, and does not impact the phase evolution of gravitational waves from binary neutron stars.
The nonlinear saturation of the r-mode instability and its effects on the spin evolution of Low Mass X-ray Binaries (LMXBs) are modeled using the triplet of modes at the lowest parametric instability threshold. We solve numerically the coupled equations for the three mode amplitudes in conjunction with the spin and temperature evolution equations. We observe that very quickly the mode amplitudes settle into quasi-stationary states. Once these states are reached, the mode amplitudes can be found algebraically and the system of equations is reduced from eight to two equations: spin and temperature evolution. Eventually, the system may reach thermal equilibrium and either (1) undergo a cyclic evolution with a frequency change of at most 10%, (2) evolve toward a full equilibrium state in which the accretion torque balances the gravitational radiation emission, or (3) enter a thermogravitational runaway on a very long timescale of about $10^6$ years. Alternatively, a faster thermal runaway (timescale of about 100 years) may occur. The sources of damping considered are shear viscosity, hyperon bulk viscosity and boundary layer viscosity. We vary proprieties of the star such as the hyperon superfluid transition temperature T_c, the fraction of the star that is above the threshold for direct URCA reactions, and slippage factor, and map the different scenarios we obtain to ranges of these parameters. For all our bound evolutions the r-mode amplitude remains small $sim 10^{-5}$. The spin frequency is limited by boundary layer viscosity to $ u_{max} sim 800 Hz [S_{ns}/(M_{1.4} R_6)]^{4/11} T_8^{-2/11}$. We find that for $ u > 700$ Hz the r-mode instability would be active for about 1 in 1000 LMXBs and that only the gravitational waves from LMXBs in the local group of galaxies could be detected by advanced LIGO interferometers.