No Arabic abstract
We study the dynamical evolution of a large amplitude r-mode by numerical simulations. R-modes in neutron stars are unstable growing modes, driven by gravitational radiation reaction. In these simulations, r-modes of amplitude unity or above are destroyed by a catastrophic decay: A large amplitude r-mode gradually leaks energy into other fluid modes, which in turn act nonlinearly with the r-mode, leading to the onset of the rapid decay. As a result the r-mode suddenly breaks down into a differentially rotating configuration. The catastrophic decay does not appear to be related to shock waves at the stars surface. The limit it imposes on the r-mode amplitude is significantly smaller than that suggested by previous fully nonlinear numerical simulations.
Previously we found that large amplitude $r$-modes could decay catastrophically due to nonlinear hydrodynamic effects. In this paper we found the particular coupling mechanism responsible for this catastrophic decay, and identified the fluid modes involved. We find that for a neutron star described by a polytropic equation of state with polytropic index $Gamma=2$, the coupling strength of the particular three-mode interaction causing the decay is strong enough that the usual picture of the $r$-mode instability with a flow pattern dominated by that of an $r$-mode can only be valid for the dimensionless $r$-mode amplitude less than $O(10^{-2})$.
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.
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.
Superfluid hydrodynamics affects the spin-evolution of mature neutron stars, and may be key to explaining timing irregularities such as pulsar glitches. However, most models for this phenomenon exclude the global instability required to trigger the event. In this paper we discuss a mechanism that may fill this gap. We establish that small scale inertial r-modes become unstable in a superfluid neutron star that exhibits a rotational lag, expected to build up due to vortex pinning as the star spins down. Somewhat counterintuitively, this instability arises due to the (under normal circumstances dissipative) vortex-mediated mutual friction. We explore the nature of the superfluid instability for a simple incompressible model, allowing for entrainment coupling between the two fluid components. Our results recover a previously discussed dynamical instability in systems where the two components are strongly coupled. In addition, we demonstrate for the first time that the system is secularly unstable (with a growth time that scales with the mutual friction) throughout much of parameter space. Interestingly, large scale r-modes are also affected by this new aspect of the instability. We analyse the damping effect of shear viscosity, which should be particularly efficient at small scales, arguing that it will not be sufficient to completely suppress the instability in astrophysical systems.
We compute the internal modes of a non-spinning neutron star and its tidal metric perturbation in general relativity, and determine the effect of relativistic corrections to the modes on mode coupling and the criterion for instability. Claims have been made that a new hydrodynamic instability can occur in a neutron star in a binary neutron star system triggered by the nonlinear coupling of the companions tidal field to pairs of p-modes and g-modes in it as the binary inspirals toward merger. This PG instability may be significant since it can influence the binarys inspiral phase by extracting orbital energy, thereby potentially causing large deviations in their gravitational waveforms from those predicted by theoretical models that do not account for it. This can result in incorrect parameter estimation, at best, or mergers going undetected, at worst, owing to the use of deficient waveform models. On the other hand, better modeling of this instability and its effect on binary orbits can unravel a new phenomenon and shed light on stellar instabilities, via gravitational wave observations. So far, all mode-tide coupling instability studies have been formulated in Newtonian perturbation theory. Neutron stars are compact objects, so relativistic corrections might be important. We present and test a new code to calculate the relativistic eigenmodes of nonrotating relativistic stars. We use these relativistic tide and neutron star eigenmodes to compute the mode-tide coupling strength (MTCS) for a few selected equations of state. The MTCS thus calculated can be at most tens of percent different from its purely Newtonian value, but we confirm the dependencies on orbital separation and equation of state found by Newtonian calculations. For some equations of state, the MTCS is very sensitive to the neutron star crust region, demonstrating the importance of treating this region accurately.