A numerical study of the r-mode instability of rapidly rotating nascent neutron stars

The first results of numerical analysis of classical r-modes of rapidly rotating compressible stellar models are reported. The full set of linear perturbation equations of rotating stars in Newtonian gravity is solved numerically without the slow rotation approximation. A critical curve of gravitational wave emission induced instability, which restricts the rotational frequencies of hot young neutron stars, is obtained. Taking the standard cooling mechanisms of neutron stars into account, we also show the 'evolutionary curves' along which neutron stars are supposed to evolve as cooling and spinning down proceed. Rotational frequencies of 1.4-M_⊙ stars suffering from this instability decrease to around 100 Hz when the standard cooling mechanism of neutron stars is employed. This result confirms the results of other authors, who adopted the slow rotation approximation.     

Low T/|W| dynamical instability in differentially rotating stars: diagnosis with canonical angular momentum

We study the nature of non-axisymmetric dynamical instabilities in differentially rotating stars with both linear eigenmode analysis and hydrodynamic simulations in Newtonian gravity. We especially investigate the following three types of instability; the one-armed spiral instability, the low   T /| W |  bar instability, and the high   T /| W |  bar instability, where T is the rotational kinetic energy and W is the gravitational potential energy. The nature of the dynamical instabilities is clarified by using a canonical angular momentum as a diagnostic. We find that the one-armed spiral and the low   T /| W |  bar instabilities occur around the corotation radius, and they grow through the inflow of canonical angular momentum around the corotation radius. The result is a clear contrast to that of a classical dynamical bar instability in high   T /| W |  . We also discuss the feature of gravitational waves generated from these three types of instability.     

Models of rapidly rotating neutron stars: remnants of accretion-induced collapse

Equilibrium models of differentially rotating nascent neutron stars are constructed, which represent the result of the accretion-induced collapse of rapidly rotating white dwarfs. The models are built in a two-step procedure: (1) a rapidly rotating pre-collapse white dwarf model is constructed; (2) a stationary axisymmetric neutron star having the same total mass and angular momentum distribution as the white dwarf is constructed. The resulting collapsed objects consist of a high-density central core of size roughly 20 km, surrounded by a massive accretion torus extending over 1000 km from the rotation axis. The ratio of the rotational kinetic energy to the gravitational potential energy of these neutron stars ranges from 0.13 to 0.26, suggesting that some of these objects may have a non-axisymmetric dynamical instability that could emit a significant amount of gravitational radiation.     

Bulk viscosity of mixed nucleon–hyperon–quark matter in neutron stars

We calculate the coefficient of bulk viscosity by considering the non-leptonic weak interactions in the cores of hybrid stars with both hyperons and quarks. We first determine the dependence of the production rate of neutrons on the reaction rate of quarks in the non-leptonic processes, that is,  Γ _n = K _s Γ _s +Γ_Λ+ 2Γ_Σ−  . The conversion rate,   K _s   , in our scenario is a complicated function of baryon number density. We also consider the medium effect of quark matter on bulk viscosity. Using these results, we estimate the limiting rotation of the hybrid stars, which may suppress the r-mode instability more effectively. Hybrid stars should be the candidates for the extremely rapid rotators.     

Pseudo-Newtonian models for the equilibrium structures of rotating relativistic stars

We obtain equilibrium solutions for rotating compact stars, including special relativistic effects. The gravity is assumed to be Newtonian, but we use the active mass density, which takes into account all energies such as the motion of the fluid, internal energy and pressure energy in addition to the rest-mass energy, in computing the gravitational potential using Poisson's equation. Such a treatment could be applicable to neutron stars with relativistic motions or a relativistic equation of state. We applied Hachisu's self-consistent field (SCF) method to find spheroidal as well as toroidal sequences of equilibrium solutions. Our solutions show better agreement with general relativistic solutions than the Newtonian relativistic hydrodynamic approach, which does not take into account the active mass. Physical quantities such as the peak density and equatorial radii in our solutions agree with the general relativistic ones to within 5 per cent. Therefore our approach can be used as a simple alternative to the fully relativistic one when a large number of model calculations is necessary, as it requires much fewer computational resources.     

Crust–core coupling and r-mode damping in neutron stars: a toy model

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 & Ushomirsky (BU) under the assumption of a perfectly rigid crust, sets the maximum spin frequency for neutron stars spun up by accretion in low-mass X-ray binaries (LMXBs). 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 having a crust with a constant shear modulus. 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 with the rigid-crust limit. We therefore revise down, by as much as a factor of 10²–10³ , the damping rate computed by BU, significantly reducing the maximal possible spin frequency of neutron stars with solid crusts. 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 multivalued for some temperatures; this is related to avoided crossings between the r-mode and higher-order torsional modes in the crust. The critical frequencies are comparable to the observed spins of neutron stars in LMXBs and millisecond pulsars.     

Axisymmetric modes of rotating relativistic stars in the Cowling approximation

Axisymmetric pulsations of rotating neutron stars can be excited in several scenarios, such as core collapse, crust- and core-quakes or binary mergers, and could become detectable in either gravitational waves or high-energy radiation. Here, we present a comprehensive study of all low-order axisymmetric modes of uniformly and rapidly rotating relativistic stars. Initial stationary configurations are appropriately perturbed and are numerically evolved using an axisymmetric, non-linear relativistic hydrodynamics code, assuming time-independence of the gravitational field (Cowling approximation). The simulations are performed using a high-resolution shock-capturing finite-difference scheme accurate enough to maintain the initial rotation law for a large number of rotational periods, even for stars at the mass-shedding limit. Through Fourier transforms of the time evolution of selected fluid variables, we compute the frequencies of quasi-radial and non-radial modes with spherical harmonic indices l =0 , 1, 2 and 3, for a sequence of rotating stars from the non-rotating limit to the mass-shedding limit. The frequencies of the axisymmetric modes are affected significantly by rotation only when the rotation rate exceeds about 50 per cent of the maximum allowed. As expected, at large rotation rates, apparent mode crossings between different modes appear. In addition to the above modes, several axisymmetric inertial modes are also excited in our numerical evolutions.     

Relativistic simulations of the phase-transition-induced collapse of neutron stars

An increase in the central density of a neutron star may trigger a phase transition from hadronic matter to deconfined quark matter in the core, causing it to collapse to a more compact hybrid star configuration. We present a study of this, building on previous work by Lin et al.. We follow them in considering a supersonic phase transition and using a simplified equation of state, but our calculations are general relativistic (using 2D simulations in the conformally flat approximation) as compared with their 3D Newtonian treatment. We also improved the treatment of the initial phase transformation, avoiding the introduction of artificial convection. As before, we find that the emitted gravitational wave spectrum is dominated by the fundamental quasi-radial and quadrupolar pulsation modes but the strain amplitudes are much smaller than suggested previously, which is disappointing for the detection prospects. However, we see significantly smaller damping and observe a non-linear mode resonance which substantially enhances the emission in some cases. We explain the damping mechanisms operating, giving a different view from the previous work. Finally, we discuss the detectability of the gravitational waves, showing that the signal-to-noise ratio for current or second generation interferometers could be high enough to detect such events in our Galaxy, although third generation detectors would be needed to observe them out to the Virgo cluster, which would be necessary for having a reasonable event rate.     

On the solution space of differentially rotating neutron stars in general relativity

A highly accurate, multidomain spectral code is used in order to construct sequences of general relativistic, differentially rotating neutron stars in axisymmetry and stationarity. For bodies with a spheroidal topology and a homogeneous or an   N = 1  polytropic equation of state, we investigate the solution space corresponding to broad ranges of degree of differential rotation and stellar densities. In particular, starting from static and spherical configurations, we analyse the changes of the corresponding surface shapes as the rate of rotation is increased. For a sufficiently weak degree of differential rotation, the sequences terminate at a mass-shedding limit, while for moderate and strong rates of differential rotation they exhibit a continuous parametric transition to a regime of toroidal fluid bodies. In this article, we concentrate on the appearance of this transition, analyse in detail its occurrence and show its relevance for the calculation of astrophysical sequences. Moreover, we find that the solution space contains various types of spheroidal configurations, which were not considered in previous work, mainly due to numerical limitations.