Seismic signatures of strange stars with crust

We study acoustic oscillations (eigenfrequencies, velocity distributions, damping times) of normal crusts of strange stars. These oscillations are very specific because of huge density jump at the interface between the normal crust and the strange matter core. The oscillation problem is shown to be self-similar. For a low (but non-zero) multipolarity l , the fundamental mode (without radial nodes) has a frequency of ∼300 Hz and mostly horizontal oscillation velocity; other pressure modes have frequencies ≳20 kHz and almost radial oscillation velocities. The latter modes are similar to radial oscillations (having approximately the same frequencies and radial velocity profiles). The oscillation spectrum of strange stars with crust differs from the spectrum of neutron stars. If detected, acoustic oscillations would allow one to discriminate between strange stars with crust and neutron stars and constrain the mass and radius of the star.     

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.     

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.     

The magnetohydrodynamics model of twin kilohertz quasi-periodic oscillations in low-mass X-ray binaries

We suggest an explanation for the twin kilohertz quasi-periodic oscillations (kHz QPOs) in low-mass X-ray binaries (LMXBs) based on magnetohydrodynamics (MHD) oscillation modes in neutron star magnetospheres. Including the effect of the neutron star spin, we derive several MHD wave modes by solving the dispersion equations, and propose that the coupling of the two resonant MHD modes may lead to the twin kHz QPOs. This model naturally relates the upper, lower kHz QPO frequencies with the spin frequencies of the neutron stars, and can well account for the measured data of six LMXBs.     

Ekman layer damping of r modes revisited

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, for example, 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 2 (depending on the detailed equation of state).     

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.     

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.     

Natural MHD oscillations of a neutron star

Natural, low-frequency, hydromagnetic oscillations of an isolated, nonrotating neutron star, which are localized in the peripheral crust, the structure of which is determined by the electron-nuclear plasma (the Ae phase), are studied. The plasma medium of the outer crust is treated as a homogeneous, infinitely conducting, incompressible continuum, the motions of which are determined by the equations of magnetohydrodynamics. In the approximation of a constant magnetic field inside the crust (the magnetic field outside the star is assumed to have a dipole structure), the spectrum of normal poloidal and toroidal hydromagnetic oscillations, due to presumed residual fluctuations of flow and their associated fluctuations in magnetic field strength, is calculated. Numerical estimates given for the periods of MHD oscillations fall in the range of periods of radio pulsar emission, indicating a close connection between the residual hydromagnetic oscillations and the electromagnetic activity of neutron stars. Translated from Astrofizika, Vol. 40, No. 1, pp. 77–86, January–March, 1997.     

Impact of neutron star oscillations on the accelerating electric field in the polar cap of pulsar

Pulsar “standard model”, that considers a pulsar as a rotating magnetized conducting sphere surrounded by plasma, is generalized to the case of oscillating star. We developed an algorithm for calculation of the Goldreich-Julian charge density for this case. We consider distortion of the accelerating zone in the polar cap of pulsar by neutron star oscillations. It is shown that for oscillation modes with high harmonic numbers (l,m) changes in the Goldreich-Julian charge density caused by pulsations of neutron star could lead to significant altering of an accelerating electric field in the polar cap of pulsar. In the moderately optimistic scenario, that assumes excitation of the neutron star oscillations by glitches, it could be possible to detect altering of the pulsar radioemission due to modulation of the accelerating field. This work was partially supported by RFBR grant 04-02-16720, and by the grants N.Sh.-5218.2006.2 and RNP-2.1.1.5940.     

High frequency oscillations during magnetar flares

The recent discovery of high frequency oscillations during giant flares from the Soft Gamma Repeaters SGR 1806-20 and SGR 1900+14 may be the first direct detection of vibrations in a neutron star crust. If this interpretation is correct it offers a novel means of testing the neutron star equation of state, crustal breaking strain, and magnetic field configuration. We review the observational data on the magnetar oscillations, including new timing analysis of the SGR 1806-20 giant flare using data from the Ramaty High Energy Solar Spectroscopic Imager and the Rossi X-ray Timing Explorer. We discuss the implications for the study of neutron star structure and crust thickness, and outline areas for future investigation.      

Freely precessing neutron stars: model and observations

We present a model of a freely precessing neutron star, which is then compared against pulsar observations. The aim is to draw conclusions regarding the structure of the star, and to test theoretical ideas of crust–core coupling and superfluidity. We argue that, on theoretical grounds, it is likely that the core neutron superfluid does not participate in the free precession of the crust. We apply our model to the handful of proposed observations of free precession that have appeared in the literature. Assuming crust-only precession, we find that all but one of the observations are consistent with there being no pinned crustal superfluid at all; the maximum amount of pinned superfluid consistent with the observations is about 10⁻¹⁰ of the total stellar moment of inertia. However, the observations do not rule out the possibility that the crust and neutron superfluid core precess as a single unit. In this case the maximum amount of pinned superfluid consistent with the observations is about 10⁻⁸ of the total stellar moment of inertia. Both of these values are many orders of magnitude less than the 10⁻² value predicted by many theories of pulsar glitches. We conclude that superfluid pinning, at least as it affects free precession, needs to be reconsidered.     

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.     

Wave beams in the magnetoelastic crust of pulsars

The propagation of axially symmetric magnetoelastic waves near the equatorial plane of the crust of a neutron star embedded in a transverse magnetic field is examined. The crust is treated as a solid-state plasma and waves are excited in it in the form of a transverse magnetic field applied to the inner boundary of the star’s crust. The time dependent equation is solved in a linear approximation assuming that the perturbing magnetic field is small compared to the unperturbed field. A simple, exact solution in the form of linear gaussian beams is obtained without additional conditions being imposed on the dissipation, dispersion, and narrowness of the beam, provided only that the velocity c_n of these waves depends weakly on position. This last condition is satisfied for the plasma in the crust of a neutron star. As it propagates to the star’s surface, the radius of the beam remains constant. The electric currents generated by the wave beam on the star’s surface are also calculated. __________ Translated from Astrofizika, Vol. 50, No. 4, pp. 547–556 (November 2007).     

Spin-down of millisecond pulsars owing to gravitational wave emission

The gravitational radiation from millisecond pulsars owing to glitches in their angular velocity is examined. It is assumed that the energy transferred from interior superfluid regions to the crust of a neutron star is converted into gravitational wave energy by damping oscillations of the matter in the star. The gravitational wave intensity and amplitude are calculated for fourteen millisecond pulsars. Gravitational radiation can explain the observed spin-down of millisecond pulsars and an estimate is given for the magnetic field at which the proposed mechanism predominates in the spin-down of these pulsars. __________ Translated from Astrofizika, Vol. 51, No. 3, pp. 479–486 (August 2008).     

Oscillations in the neutron star crust

We investigate the spectrum of torsional modes in the neutron star crust and discuss what conclusions may be drawn about the global properties of the star from observations of such modes.      

Maximum mass of a hot neutron star with a quark core

We have considered a hot neutron star with a quark core,a mixed phase of quark-hadron matter,and a hadronic matter crust and have determined the equation of state of the hadronic phase and the quark phase.We have then found the equation of state of the mixed phase under the Gibbs conditions.Finally,we have computed the structure of a hot neutron star with a quark core and compared our results with those of the neutron star without a quark core.For the quark matter calculations,we have used the MIT bag model...     

Quasi-periodic oscillations in low-mass X-ray binaries and constraints on the equation of state of neutron star matter

Recently discovered quasi-periodic oscillations in the X-ray brightness of low-mass X-ray binaries are used to derive constraints on the mass of the neutron star component and the equation of state of neutron star matter. The observations are compared with models of rapidly rotating neutron stars which are calculated by means of an exact numerical method in full relativity. For the equations of state we select a broad collection of models representing different assumptions about the many-body structure and the complexity of the composition of superdense matter. The mass constraints differ from their values in the approximate treatment by ∼10 per cent. Under the assumption that the maximum frequency of the quasi-periodic oscillations originates from the innermost stable orbit, the mass of the neutron star is in the range M ∼1.92–2.25 M_⊙. The quasi-periodic oscillation in the Atoll-source 4U 1820−30 in particular is only consistent with equations of state that are rather stiff at high densities, which is explainable, so far, only with pure nucleonic/leptonic composition. This interpretation contradicts the hypothesis that the protoneutron star formed in SN 1987A collapsed to a black hole, since this would demand a maximum neutron star mass below 1.6 M_⊙. The recently suggested identification of quasi-periodic oscillations with frequencies of about 10 Hz with the Lense–Thirring precession of the accretion disc is found to be inconsistent with the models studied in this work, unless it is assumed that the first overtone of the precession is observed.     

The inverse problem for pulsating neutron stars: a 'fingerprint analysis' for the supranuclear equation of state

We study the problem of detecting, and inferring astrophysical information from, gravitational waves from a pulsating neutron star. We show that the fluid f and p modes, as well as the gravitational-wave w modes, may be detectable from sources in our own Galaxy, and investigate how accurately the frequencies and damping rates of these modes can be inferred from a noisy gravitational-wave data stream. Based on the conclusions of this discussion we propose a strategy for revealing the supranuclear equation of state using the neutron star fingerprints: the observed frequencies of an f and a p mode. We also discuss how well the source can be located in the sky using observations with several detectors.     

Constraints on the magnetic field geometry of magnetars

We study the effect of the magnetic field geometry on the oscillation spectra of strongly magnetized stars. We construct a configuration of magnetic field where a toroidal component is added to the standard poloidal one. We consider a star with a type I superconductor core so that both components of the magnetic field are expelled from the core and confined in the crust. Our results show that the toroidal contribution does not influence significantly the torsional oscillations of the crust. On the contrary, the confinement of the magnetic field in the crust drastically affects the torsional oscillation spectrum. A comparison with estimations for the magnetic field strength, from observations, excludes the possibility that magnetars will have a magnetic field solely confined in the crust, that is, our results suggest that the magnetic field in whatever geometry has to permeate the whole star.     

Vortex lattice oscillations in rotating neutron stars with quark “CFL” cores

Collective elastic oscillations of a lattice of nonabelian quark semisuperfluid vortex filaments in the superfluid core of a rotating neutron star are examined. It is shown that in the incompressible fluid approximation, transverse long wavelength oscillations (Tkachenko oscillations) owing to shear deformation of the vortex lattice propagate in a plane perpendicular to the axis of rotation. The periods of these oscillations are consistent with rotational variations on the order of 100-1000 days observed in the pulsars PSR B0531+21 and PSR B1828-11. Translated from Astrofizika, Vol. 52, No. 1, pp. 165–169 (February 2009).