Spin evolution of accreting neutron stars: Nonlinear development of the r-mode instability
- Center for Radiophysics and Space Research, Cornell University, Ithaca, New York 14853 (United States)
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 modes in conjunction with the spin and temperature evolution equations. We observe that very quickly the mode amplitudes settle into quasistationary states that change slowly as the temperature and spin of the star evolve. 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. The evolution of the neutron star angular velocity and temperature follow easily calculated trajectories along these sequences of quasistationary states. The outcome depends on whether or not the star will reach thermal equilibrium, where the viscous heating by the three modes is equal to the neutrino cooling (H=C curve). If, when the r-mode becomes unstable, the star spins at a frequency below the maximum of the H=C curve, then it will reach a state of thermal equilibrium. It can then either (1) undergo a cyclic evolution with a small cycle size 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 time scale of {approx_equal}10{sup 6} years. If the star does not reach a state of thermal equilibrium, then a faster thermal runaway (time scale of {approx_equal}100 years) occurs and the r-mode amplitude increases above the second parametric instability threshold. Following this evolution requires more inertial modes to be included. The sources of damping considered are shear viscosity, hyperon bulk viscosity, and viscosity within the core-crust boundary layer. We vary proprieties of the star such as the hyperon superfluid transition temperature T{sub 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. We focus on T{sub c} > or approx. 5x10{sup 9} K where nonlinear effects are important. Wagoner [R. Wagoner, Astrophys. J. 578, L63 (2002).] has shown that a very low r-mode amplitude arises at smaller T{sub c}. For all our bounded evolutions the r-mode amplitude remains small {approx}10{sup -5}. The spin frequency of accreting neutron stars is limited by boundary layer viscosity to {nu}{sub max}{approx_equal}800 Hz[S{sub ns}/(M{sub 1.4}R{sub 6})]{sup 4/11}T{sub 8}{sup -2/11}. Fast rotators are allowed for [S{sub ns}/(M{sub 1.4}R{sub 6})]{sup 4/11}T{sub 8}{sup -2/11}{approx}1 and we find that in this case 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.
- OSTI ID:
- 21032249
- Journal Information:
- Physical Review. D, Particles Fields, Vol. 76, Issue 6; Other Information: DOI: 10.1103/PhysRevD.76.064019; (c) 2007 The American Physical Society; Country of input: International Atomic Energy Agency (IAEA); ISSN 0556-2821
- Country of Publication:
- United States
- Language:
- English
Similar Records
Merger of black hole and neutron star in general relativity: Tidal disruption, torus mass, and gravitational waves
Effects of background gas pressure on the dynamics of a nonneutral electron plasma confined in a Malmberg-Penning trap
Related Subjects
AMPLITUDES
ANGULAR VELOCITY
BINARY STARS
BOUNDARY LAYERS
COSMOLOGY
EMISSION
GRAVITATIONAL WAVES
HEATING
HYPERONS
MASS
NEUTRINOS
NEUTRON STARS
NONLINEAR PROBLEMS
PARAMETRIC INSTABILITIES
SPIN
SUPERFLUIDITY
THERMAL EQUILIBRIUM
TRANSITION TEMPERATURE
VISCOSITY
X RADIATION