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.     

General relativistic spectra of accretion discs around rapidly rotating neutron stars: effect of light bending

We present computed spectra, as seen by a distant observer, from the accretion disc around a rapidly rotating neutron star. Our calculations are carried out in a fully general relativistic framework, with an exact treatment of rotation. We take into account the Doppler shift, gravitational redshift and light-bending effects in order to compute the observed spectrum. We find that light bending significantly modifies the high-energy part of the spectrum. Computed spectra for slowly rotating neutron stars are also presented. These results would be important for modelling the observed X-ray spectra of low-mass X-ray binaries containing fast-spinning neutron stars.     

Oscillations of rapidly rotating superfluid stars

Using time evolutions of the relevant linearized equations, we study non-axisymmetric oscillations of rapidly rotating and superfluid neutron stars. We consider perturbations of Newtonian axisymmetric background configurations and account for the presence of superfluid components via the standard two-fluid model. Within the Cowling approximation, we are able to carry out evolutions for uniformly rotating stars up to the mass-shedding limit. This leads to the first detailed analysis of superfluid neutron star oscillations in the fast rotation regime, where the star is significantly deformed by the centrifugal force. For simplicity, we focus on background models where the two fluids (superfluid neutrons and protons) corotate, are in β-equilibrium and co-exist throughout the volume of the star. We construct sequences of rotating stars for two analytical model equations of state. These models represent relatively simple generalizations of single fluid, polytropic stars. We study the effects of entrainment, rotation and symmetry energy on non-radial oscillations of these models. Our results show that entrainment and symmetry energy can have a significant effect on the rotational splitting of non-axisymmetric modes. In particular, the symmetry energy modifies the inertial mode frequencies considerably in the regime of fast rotation.     

On the r-mode spectrum of relativistic stars in the low-frequency approximation

The axial modes for non-barotropic relativistic rotating neutron stars with uniform angular velocity are studied, using the slow-rotation formalism together with the low-frequency approximation, first investigated by Kojima. The time-independent form of the equations leads to a singular eigenvalue problem, which admits a continuous spectrum. We show that for     , it is nevertheless also possible to find discrete mode solutions (the r modes). However, under certain conditions related to the equation of state and the compactness of the stellar model, the eigenfrequency lies inside the continuous band and the associated velocity perturbation is divergent; hence these solutions have to be discarded as being unphysical. We corroborate our results by explicitly integrating the time-dependent equations. For stellar models admitting a physical r-mode solution, it can indeed be excited by arbitrary initial data. For models admitting only an unphysical mode solution, the evolutions do not show any tendency to oscillate with the respective frequency. For higher values of l it seems that in certain cases there are no mode solutions at all.     

Equilibrium configurations of strongly magnetized neutron stars with realistic equations of state

We investigate equilibrium sequences of magnetized rotating stars with four kinds of realistic equations of state (EOSs) of SLy, FPS, Shen and LS, employing the Tomimura–Eriguchi scheme to construct the equilibrium configurations. We study the basic physical properties of the sequences in the framework of Newtonian gravity. In addition, we take a new step by taking into account a general relativistic effect to the magnetized rotating configurations. With these computations, we find that the properties of the Newtonian magnetized stars, e.g. structure of magnetic field, highly depends on the EOSs. The toroidal magnetic fields concentrate rather near the surface for Shen and LS EOSs than those for SLy and FPS EOSs. The poloidal fields are also affected by the toroidal configurations. Paying attention to the stiffness of the EOSs, we analyse this tendency in detail. In the general relativistic stars, we find that the difference due to the EOSs becomes small because all the employed EOSs become sufficiently stiff for the large maximum density, typically greater than  10¹⁵ g cm⁻³  . The maximum baryon mass of the magnetized stars with axis ratio   q ∼ 0.7  increases about up to 20 per cent for that of spherical stars. We furthermore compute equilibrium sequences at finite temperature, which should serve as an initial condition for the hydrodynamic study of newly born magnetars. Our results suggest that we may obtain information about the EOSs from the observation of the masses of magnetars.     

A General Relativistic Calculation of Boundary Layer and Disc Luminosity for Accreting Non-magnetic Neutron Stars in Rapid Rotation

We calculate the disc and boundary layer luminosities for accreting rapidly rotating neutron stars with low magnetic fields in a fully general relativistic manner. Rotation increases the disc luminosity and decreases the boundary layer luminosity. A rapid rotation of the neutron star substantially modifies these quantities as compared with the static limit. For a neutron star rotating close to the centrifugal mass shed limit, the total luminosity has contribution only from the extended disc. For such maximal rotation rates, we find that well before the maximum stable gravitational mass configuration is reached, there exists a limiting central density, for which particles in the innermost stable orbit will be more tightly bound than those at the surface of the neutron star. We also calculate the angular velocity profiles of particles in Keplerian orbits around the rapidly rotating neutron star. The results are illustrated for a representative set of equation of state models of neutron star matter.     

The magnetosphere of oscillating neutron stars in general relativity

Just as a rotating magnetized neutron star has material pulled away from its surface to populate a magnetosphere, a similar process can occur as a result of neutron-star pulsations rather than rotation. This is of interest in connection with the overall study of neutron star oscillation modes but with a particular focus on the situation for magnetars. Following a previous Newtonian analysis of the production of a force-free magnetosphere in this way Timokhin et al., we present here a corresponding general-relativistic analysis. We give a derivation of the general relativistic Maxwell equations for small-amplitude arbitrary oscillations of a non-rotating neutron star with a generic magnetic field and show that these can be solved analytically under the assumption of low current density in the magnetosphere. We apply our formalism to toroidal oscillations of a neutron star with a dipole magnetic field and find that the low current density approximation is valid for at least half of the oscillation modes, similarly to the Newtonian case. Using an improved formula for the determination of the last closed field line, we calculate the energy losses resulting from toroidal stellar oscillations for all of the modes for which the size of the polar cap is small. We find that general relativistic effects lead to shrinking of the size of the polar cap and an increase in the energy density of the outflowing plasma. These effects act in opposite directions but the net result is that the energy loss from the neutron star is significantly smaller than suggested by the Newtonian treatment.     

Magnetic and spin evolution of neutron stars in close binaries

The evolution of neutron stars in close binary systems with a low-mass companion is considered, assuming the magnetic field to be confined within the solid crust. We adopt the standard scenario for the evolution in a close binary system, in which the neutron star passes through four evolutionary phases ('isolated pulsar'–'propeller'– accretion from the wind of a companion – accretion resulting from Roche-lobe overflow). Calculations have been performed for a great variety of parameters characterizing the properties of both the neutron star and the low-mass companion. We find that neutron stars with more or less standard magnetic field and spin period that are processed in low-mass binaries can evolve to low-field rapidly rotating pulsars. Even if the main-sequence life of a companion is as long as 10¹⁰ yr, the neutron star can maintain a relatively strong magnetic field to the end of the accretion phase. The model that is considered can account well for the origin of millisecond pulsars.     

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.     

On the r-mode spectrum of relativistic stars: the inclusion of the radiation reaction

We consider both mode calculations and time-evolutions of axial r modes for relativistic uniformly rotating non-barotropic neutron stars, using the slow-rotation formalism, in which rotational corrections are considered up to linear order in the angular velocity Ω. We study various stellar models, such as uniform density models, polytropic models with different polytropic indices n , and some models based on realistic equations of state. For weakly relativistic uniform density models and polytropes with small values of n , we can recover the growth times predicted from Newtonian theory when standard multipole formulae for the gravitational radiation are used. However, for more compact models, we find that relativistic linear perturbation theory predicts a weakening of the instability compared to the Newtonian results. When turning to polytropic equations of state, we find that for certain ranges of the polytropic index n , the r mode disappears, and instead of a growth, the time-evolutions show a rapid decay of the amplitude. This is clearly at variance with the Newtonian predictions. It is, however, fully consistent with our previous results obtained in the low-frequency approximation.     

Crustal oscillations of slowly rotating relativistic stars

We study low-amplitude crustal oscillations of slowly rotating relativistic stars consisting of a central fluid core and an outer thin solid crust. We estimate the effect of rotation on the torsional toroidal modes and on the interfacial and shear spheroidal modes. The results compared against the Newtonian ones for wide range of neutron star models and equations of state.     

Matching the frequency spectrum of pre-main sequence stars by means of standard and rotating models

We applied the aton evolutionary code to the computation of detailed grids of standard (non-rotating) and rotating pre-main sequence (PMS) models and computed their adiabatic oscillation spectra, with the aim of exploring the seismic properties of young stars. As, until now, only a few frequencies have been determined for ∼40 PMS stars, the way of approaching the interpretation of the oscillations is not unique. We adopt a method similar to the matching mode method by Guenther and Brown making use, when necessary, also of our rotating evolutionary code to compute the models for PMS stars. The method is described by a preliminary application to the frequency spectrum of two PMS stars (85 and 278) in the young open cluster NGC 6530. For the Star 85, we confirm with self-consistent rotating models, previous interpretation of the data, attributing three close frequencies to the mode   n = 4, l = 1  and   m = 0  , +1 and −1. For the Star 278, we find a different fit for the frequencies, corresponding to a model within the original error box of the star, and dispute the possibility that this star has a T _eff much cooler that the red boundary of the radial instability strip.     

On the dynamics of superfluid neutron star cores

We discuss the nature of the various modes of pulsation of superfluid neutron stars using comparatively simple Newtonian models and the Cowling approximation. The matter in these stars is described in terms of a two-fluid model, where one fluid is the neutron superfluid, which is believed to exist in the core and inner crust of mature neutron stars, and the other fluid represents a conglomerate of all other constituents (crust nuclei, protons, electrons, etc.). In our model, we incorporate the non-dissipative interaction known as the entrainment effect, whereby the momentum of one constituent (e.g. the neutrons) carries along part of the mass of the other constituent. We show that there is no independent set of pulsating g-modes in a non-rotating superfluid neutron star core, even though the linearized superfluid equations contain a well-defined (and real-valued) analogue to the so-called Brunt–Väisälä frequency. Instead, what we find are two sets of spheroidal perturbations whose nature is predominately acoustic. In addition, an analysis of the zero-frequency subspace (i.e. the space of time-independent perturbations) reveals two sets of degenerate spheroidal perturbations, which we interpret to be the missing g-modes, and two sets of toroidal perturbations. We anticipate that the degeneracy of all these zero-frequency modes will be broken by the Coriolis force in the case of rotating stars. To illustrate this we consider the toroidal pulsation modes of a slowly rotating superfluid star. This analysis shows that the superfluid equations support a new class of r-modes, in addition to those familiar from, for example, geophysical fluid dynamics. Finally, the role of the entrainment effect on the superfluid mode frequencies is shown explicitly via solutions to dispersion relations that follow from a 'local' analysis of the linearized superfluid equations.     

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.     

Young Close-By Neutron Stars: The Gould Belt Vs. The Galactic Disc

We present new population synthesis calculations of close young neutron stars. In comparison with our previous investigation we use a different neutron star mass spectrum and different initial spatial and velocity distributions. The results confirm that most of ROSAT dim radioquiet isolated neutron stars had their origin in the Gould Belt. We predict that about several tens of young neutron stars can be identified in ROSAT All Sky Survey data at low galactic latitudes. Some of these sources also can have counterparts among EGRET unidentified sources.     

Magnetic fields in axisymmetric neutron stars

We derive general equations for axisymmetric Newtonian magnetohydrodynamics and use these as the basis of a code for calculating equilibrium configurations of rotating magnetized neutron stars in a stationary state. We investigate the field configurations that result from our formalism, which include purely poloidal, purely toroidal and mixed fields. For the mixed-field formalism, the toroidal component appears to be bounded at less than 7 per cent. We calculate distortions induced both by magnetic fields and by rotation. From our non-linear work, we are able to look at the realm of validity of perturbative work: we find for our results that perturbative-regime formulae for magnetic distortions agree to within 10 per cent of the non-linear results if the ellipticity is less than 0.15 or the average field strength is less than 10¹⁷ G. We also consider how magnetized equilibrium structures vary for different polytropic indices.     

General relativistic electromagnetic fields of a slowly rotating magnetized neutron star – II. Solution of the induction equations

We have solved numerically the general relativistic induction equations in the interior background space–time of a slowly rotating magnetized neutron star. The analytic form of these equations was discussed recently (Paper I), where corrections due to both the space–time curvature and the dragging of reference frames were shown to be present. Through a number of calculations we have investigated the evolution of the magnetic field with different rates of stellar rotation, different inclination angles between the magnetic moment and the rotation axis, as well as different values of the electrical conductivity. All of these calculations have been performed for a constant-temperature relativistic polytropic star and make use of a consistent solution of the initial-value problem which avoids the use of artificial analytic functions. Our results show that there exist general relativistic effects introduced by the rotation of the space–time which tend to decrease the decay rate of the magnetic field. The rotation-induced corrections are however generally hidden by the high electrical conductivity of the neutron star matter, and when realistic values for the electrical conductivity are considered, these corrections become negligible even for the fastest known pulsar.     

Simulations of the axisymmetric magnetospheres of neutron stars

In this paper we present the results of time-dependent simulations of the dipolar axisymmetric magnetospheres of neutron stars carried out within the frameworks of both relativistic magnetohydrodynamics (MHD) and resistive force-free electrodynamics. The results of force-free simulations reveal the inability of our numerical method to accommodate the equatorial current sheets of pulsar magnetospheres, and raise a question mark about the robustness of this approach. On the other hand, the MHD approach allows us to make significant progress. We start with a non-rotating magnetically dominated dipolar magnetosphere and follow its evolution as the stellar rotation is switched on. We find that the time-dependent solution gradually approaches a steady state that is very close to the stationary solution of the pulsar equation found in 1999 by Contopoulos, Kazanas & Fendt. This result suggests that other stationary solutions that have the Y-point located well inside the light cylinder are unstable. The role of particle inertia and pressure on the structure and dynamics of MHD magnetospheres is studied in detail, as well as the potential implications of dissipative processes in the equatorial current sheet. We argue that pulsars may have differentially rotating magnetospheres which develop noticeable structural oscillations, and that this may help to explain the nature of the subpulse phenomena.     

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.     

The physics of strong magnetic fields in neutron stars

In this paper we present a new result, namely that the primal magnetic field of the collapsed core during a supernova explosion will, as a result of the conservation of magnetic flux, receive a massive boost to more than 90 times its original value by the Pauli paramagnetization of the highly degenerate relativistic electron gas just after the formation of the neutron star. Thus, the observed super-strong magnetic field of neutron stars may originate from the induced Pauli paramagnetization of the highly degenerate relativistic electron gas in the interior of the neutron star. We therefore have an apparently natural explanation for the surface magnetic field of a neutron star.