Search

Using time evolutions of the relevant linearised 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) co-rotate, are in beta-equilibrium and coexist throughout the volume of the star. We construct sequences of rotating stars for two analytical model equations of state. These models represent relatively simple generalisations 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., Comment: 18 pages, 7 figures, 5 tables, accepted for publication in MNRAS

Motivated by the possibility of radiation driven instabilities in rotating magnetic stars, we study the stability properties of general linear perturbations of a stationary and axisymmetric, infinitely conducting perfect fluid configuration threaded by a magnetic field and surrounded by vacuum. We develop a Lagrangian perturbation framework which enables us to formulate a strict stability criterion based on the notion of a canonical energy (a functional of the fluid displacement $\xi $ and its first time derivative). For any given choice of $\{\xi,\partial_t \xi \} $, the sign of the canonical energy determines whether the configuration is stable or not at the linear level. Our analysis provides the first complete description of the stability problem for a magnetic star, allowing for both rotation and the presence of a magnetic field in the exterior vacuum region. A key feature of the Lagrangian formulation is the existence of so-called `trivial' fluid displacements, which do not represent true physical perturbations. In order for the stability criterion to make rigorous sense one has to isolate these trivials and consider only the physical `canonical' displacements. We discuss this problem and formulate a condition which must be satisfied by all canonical displacements. Having obtained a well-defined stability criterion we provide examples which indicate that the magnetic field has a stabilising effect on radiation driven instabilities., Comment: 16 pages, accepted for publication in MNRAS

Recent observations of quasi-periodic oscillations in the aftermath of giant flares in soft gamma-ray repeaters suggest a close coupling between the seismic motion of the crust after a major quake and the modes of oscillations in a magnetar. In this paper we consider the purely elastic modes of oscillation in the crust of a neutron star in the relativistic Cowling approximation (disregarding any magnetic field). We determine the axial crust modes for a large set of stellar models, using a state-of-the-art crust equation of state and a wide range of core masses and radii. We also devise useful approximate formulae for the mode-frequencies. We show that the relative crust thickness is well described by a function of the compactness of the star and a parameter describing the compressibility of the crust only. Considering the observational data for SGR 1900+14 and SGR 1806-20, we demonstrate how our results can be used to constrain the mass and radius of an oscillating neutron star., Comment: 14 pages, 1 table, 4 figures. v2: Typos corrected, references added and a bug in the code was fixed which slightly modifies the figures but leaves the conclusions unchanged

We investigate the damping of neutron star r-modes due to the presence of a viscous boundary (Ekman) layer at the interface between the crust and the core. Our study is motivated by the possibility that the gravitational-wave driven instability of the inertial r-modes may become active in rapidly spinning neutron stars, eg. in low-mass X-ray binaries, and the fact that a viscous Ekman layer at the core-crust interface provides an efficient damping mechanism for these oscillations. We review various approaches to the problem and carry out an analytic calculation of the effects due to the Ekman layer for a rigid crust. Our analytic estimates support previous numerical results, and provide further insight into the intricacies of the problem. We add to previous work by discussing the effect that compressibility and composition stratification have on the boundary layer damping. We show that, while stratification is unimportant for the r-mode problem, composition suppresses the damping rate by about a factor of two (depending on the detailed equation of state)., Comment: 16 pages, 2 eps figures. Minor corrections to match published version

We use a simple toy-model to discuss global MHD modes of a neutron star, taking into account the magnetic coupling between the elastic crust and the fluid core. Our results suggest that the notion of pure torsional crust modes is not useful for the coupled system, all modes excite Alfven waves in the core. However, we also show that the modes that are most likely to be excited by a fractured crust, eg. during a magnetar flare, are such that the crust and the core oscillate in concert. For our simple model, the frequencies of these modes are similar to the ``pure crustal'' frequencies. In addition, our model provides a natural explanation for the presence of lower frequency ($< 30$~Hz) QPOs seen in the December 2004 giant flare of SGR 1806-20., Comment: 5 pages, 2 eps figures, minor corrections to match published version

We develop a flux-conservative formalism for a Newtonian multi-fluid system, including dissipation and entrainment (i.e. allowing the momentum of one fluid to be a linear combination of the velocities of all fluids). Maximum use is made of mass, energy, and linear and angular momentum conservation to specify the equations of motion. Also used extensively are insights gleaned from a convective variational action principle, key being the distinction between each velocity and its canonically conjugate momentum. Dissipation is incorporated to second order in the ``thermodynamic forces'' via the approach pioneered by Onsager. An immediate goal of the investigation is to understand better the number, and form, of independent dissipation terms required for a consistent set of equations of motion in the multi-fluid context. A significant, but seemingly innocuous detail, is that one must be careful to isolate ``forces'' that can be written as total gradients, otherwise errors can be made in relating the net internal force to the net externally applied force. Our long-range aim is to provide a formalism that can be used to model dynamical multi-fluid systems both perturbatively and via fully nonlinear 3D numerical evolutions. To elucidate the formalism we consider the standard model for a heat-conducting, superfluid neutron star, which is believed to be dominated by superfluid neutrons, superconducting protons, and a highly degenerate, ultra-relativistic gas of normal fluid electrons. We determine that in this case there are, in principle, 19 dissipation coefficients in the final set of equations., Comment: 21 pages, uses Revtex 4

We discuss the generalisation of the so-called traditional approximation, well known in geophysics, to general relativity. We show that the approximation is applicable for rotating relativistic stars provided that one focuses on relatively thin radial shells. This means that the framework can be used to study waves in neutron star oceans. We demonstrate that, once the effects of the relativistic frame-dragging are accounted for, the angular problem reduces to Laplace's tidal equation. We derive the dispersion relation for various classes of waves in a neutron star ocean and show that the combined effects of the frame-dragging and the gravitational redshift typically lower the frequency of a mode by about 20%., Comment: 16 pages, 1 eps figure

An understanding of the dynamics of differentially rotating systems is key to many areas of astrophysics. We investigate the oscillations of a simple system exhibiting differential rotation, and discuss issues concerning the role of corotation points and the emergence of dynamical instabilities. This problem is of particular relevance to the emission of gravitational waves from oscillating neutron stars, which are expected to possess significant differential rotation immediately after birth or binary merger., Comment: 14 pages, 6 figures, minor error corrected, MNRAS in press

The loss of angular momentum due to unstable r-modes in hot young neutron stars has been proposed as a mechanism for achieving the spin rates inferred for young pulsars. One factor that could have a significant effect on the action of the r-mode instability is fallback of supernova remnant material. The associated accretion torque could potentially counteract any gravitational-wave induced spin-down, and accretion heating could affect the viscous damping rates and hence the instability. We discuss the effects of various external agents on the r-mode instability scenario within a simple model of supernova fallback on to a hot young magnetized neutron star. We find that the outcome depends strongly on the strength of the star's magnetic field. Our model is capable of generating spin rates for young neutron stars that accord well with initial spin rates inferred from pulsar observations. The combined action of r-mode instability and fallback appears to cause the spin rates of neutron stars born with very different spin rates to converge, on a timescale of about a year. The results suggest that stars with magnetic fields $\le10^{13}$ G could emit a detectable gravitational wave signal for perhaps several years after the supernova event. Stars with higher fields (magnetars) are unlikely to emit a detectable gravitational wave signal via the r-mode instability. The model also suggests that the r-mode instability could be extremely effective in preventing young neutron stars from going dynamically unstable to the bar-mode., Comment: 10 pages,5 figures, minor changes. Accepted by MNRAS

The non-linear behaviour of oscillation modes in compact stars is a topic of considerable current interest. Accurate numerical studies of such phenomena are likely to require powerful new approaches to both fluid and spacetime computations. We propose that a key ingredient of such methods will be the non-linear evolution of deviations from the background stationary equilibrium star. We investigate the feasibility of this approach by applying it to non-linear radial oscillations of a neutron star, and explore numerically various non-linear features of this problem, for a large range of amplitudes. Quadratic and higher order mode coupling and non-linear transfer of energy is demonstrated and analysed in detail., Comment: 5 pages,5 eps figures, latex (mn2e style)

We study the r-modes and rotational ``hybrid'' modes of relativistic stars. As in Newtonian gravity, the spectrum of low-frequency rotational modes is highly sensitive to the stellar equation of state. If the star and its perturbations obey the same one-parameter equation of state (as with isentropic stars), there exist {\it no pure r-modes at all} - no modes whose limit, for a star with zero angular velocity, is an axial-parity perturbation. Rotating stars of this kind similarly have no pure g-modes, no modes whose spherical limit is a perturbation with polar parity and vanishing perturbed pressure and density. We compute the post-Newtonian corrections to the $l=m$ r-modes of isentropic and non-isentropic uniform density stars., Comment: 36 pages revtex v3.1, 2 eps figures

We develop a framework for constructing initial data sets for perturbations about spherically symmetric matter distributions. This framework facilitates setting initial data representing astrophysical sources of gravitational radiation involving relativistic stars. The procedure is based on Lichnerowicz-York's conformal approach to solve the constraints in Einstein's equations. The correspondence of these initial data sets in terms of the standard gauge perturbation variables in the Regge-Wheeler perturbation variables is established, and examples of initial data sets of merging neutron stars under the close-limit approximation are presented., Comment: 14 pages, 2 figures

We present a case study for the relevance of the r-mode instability for accreting compact stars. Our estimates are based on approximations that facilitate back-of-the-envelope calculations. We discuss two different cases: 1) For recycled millisecond pulsars we argue that the r-mode instability may be active at rotation periods longer than the Kepler period (which provides the dynamical limit on rotation) as long as the core temperature is larger than about $2\times10^5$ K. Our estimates suggest that the instability may have played a role in the evolution of the fastest spinning pulsars, and that it may be presently active in the recently discovered 2.49 ms X-ray pulsar SAX J1808.4-3658 as well as the rapidly spinning neutron stars observed in low mass X-ray binaries. This provides a new explanation for the remarkably similar rotation periods inferred from kHz quasi-periodic oscillations in the LMXBs. The possibility that the rotation of recycled pulsars may be gravitational-radiation limited is interesting because the gravitational waves from a neutron star rotating at the instability limit may well be detectable with the new generation of interferometric detectors. 2) We also consider white dwarfs, and find that the r-mode instability may possibly be active in short period white dwarfs. Our order of magnitude estimates (for a white dwarf of $M=M_\odot$ and $R=0.01 R_\odot$ composed of $C^{12}$) show that the instability could be operating for rotational periods shorter than $P\approx 27-33$ s. This number is in interesting agreement with the observed periods ($>28$ s) of the rapidly spinning DQ Herculis stars. However, we find that the instability grows too slowly to affect the rotation of these stars significantly ., Comment: 11 pages, 2 figures, to appear in Ap.J

A newly discovered instability in rotating neutron stars, driven by gravitational radiation reaction acting on the stars' $r$-modes, is shown here to set an upper limit on the spin rate of young neutron stars. We calculate the timescales for growth of linear perturbations due to gravitational radiation reaction, and for dissipation by shear and bulk viscosity, working to second order in a slow-rotation expansion within a Newtonian polytropic stellar model. The results are very temperature-sensitive: in hot neutron stars ($T>10^9$ K), the lowest-order $r$-modes are unstable, while in colder stars they are damped by viscosity. These calculations have a number of interesting astrophysical implications. First, the $r$-mode instability will spin down a newly born neutron star to a period close to the initial period inferred for the Crab pulsar, probably between 10 and 20 ms. Second, as an initially rapidly rotating star star spins down, an energy equivalent to roughly 1% of a solar mass is radiated as gravitational waves, which makes the process an interesting source for detectable gravitational waves. Third, the $r$-mode instability rules out the scenario whereby millisecond pulsars are formed by accretion-induced collapse of a white dwarf: the new star would be hot enough to spin down to much slower rates., Comment: 1 ps figure

We present new results for pulsating neutron stars. We have calculated the eigenfrequencies of the modes that one would expect to be the most important gravitational-wave sources: the fundamental fluid f-mode, the first pressure p-mode and the first gravitational-wave w-mode, for twelve realistic equations of state. From this numerical data we have inferred a set of ``empirical relations'' between the mode-frequencies and the parameters of the star (the radius R and the mass M). Some of these relation prove to be surprisingly robust, and we show how they can be used to extract the details of the star (radius, mass, eos) from observed modes with errors no larger than a few percent., Comment: 11 pages, 8 figures, mn.sty

The first numerical study of axial (toroidal) pulsation modes of a slowly rotating relativistic star is presented. The calculation includes terms of first order in $\epsilon \equiv \Omega \sqrt{R^3/M}<<1$ ($R$ is the radius, $M$ is the mass and $\Omega$ is the rotation frequency of the star), and accounts for effects due to the coriolis force. Effects due to the centrifugal flattening of the star enter at order $\epsilon^2$ and are not included in the analysis. It is shown that increased rotation tends to decrease the damping times for prograde modes, while retrograde become longer lived. Specifically, we show that rotation affects the axial gravitational-wave $w$-modes in this way. We also present the first relativistic calculation of the so-called $r$-modes (analogous to Rossby waves in the Earth's oceans). These have frequencies of the same order of magnitude as the rotation frequency of the star. The presented results indicate that the $r$-modes are unstable due to the emission of gravitational radiation for \underline{all} rotating perfect fluid stars. This is interesting since the previously considered gravitational-wave instability associated with (for example) the $f$-mode of the star sets in at a critical rotation rate. Because they are unstable also for the slowest rotating stars the $r$-modes may well be of considerable astrophysical importance., Comment: 14 pages, latex

An extensive survey of gravitational-wave modes for uniform density stars is presented. The study covers stars ranging in compactness from $R/M=100$ to the limit of stability in general relativity: $R/M = 9/4$. We establish that polar and axial gravitational-wave modes exist for all these stellar models. Moreover, there are two distinct families of both axial and polar modes. We discuss the physics of these modes and argue that one family is primarily associated with the interior of the star, while the second family is mainly associated with the stellar surface. We also show that the problem of axial perturbations has all the essential features of the polar problem as far as gravitational waves are concerned. This means that the axial problem is much more important than has previously been assumed. We also find some surprising features, such as avoided crossings between the polar gravitational-wave modes and the Kelvin f-mode as the star becomes very compact. This seems to suggest that the f-mode should be considered on equal footing with the polar w-modes for ultracompact stars. All modes may have the main character of trapped modes inside the curvature potential barrier for $R/M < 3$., Comment: 24 pages, AASMS.STY, 5 postscript figures, Astroph. J., in press

We present results for r-modes of relativistic nonbarotropic stars. We show that the main differential equation, which is formally singular at lowest order in the slow-rotation expansion, can be regularized if one considers the initial value problem rather than the normal mode problem. However, a more physically motivated way to regularize the problem is to include higher order terms. This allows us to develop a practical approach for solving the problem and we provide results that support earlier conclusions obtained for uniform density stars. In particular, we show that there will exist a single r-mode for each permissible combination of 1 and m. We discuss these results and provide some caveats regarding their usefulness for estimates of gravitational-radiation reaction timescales. The close connection between the seemingly singular relativistic r-mode problem and issues arising because of the presence of co-rotation points in differentially rotating stars is also clarified.