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.     

Secular instability of g-modes in rotating neutron stars

Gravitational radiation tends to drive gravity modes in rotating neutron stars to become unstable. For an inviscid star, the instability sets in when the rotation frequency is about 0.7 times the corresponding mode frequency of the non-rotating star. Neutron stars with spin frequencies ≳100 Hz are susceptible to this instability, with a growth time of the order of years. However, it is likely that viscous dissipation suppresses the instability except for a narrow range of temperatures around 10⁹ K. We also show that the viscosity-driven instability of g-modes is absent.     

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.     

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.     

Do superfluid instabilities prevent neutron star precession?

We discuss short wavelength (inertial wave) instabilities present in the standard two-fluid neutron star model when there is sufficient relative flow along the superfluid neutron vortex array. We demonstrate that these instabilities may be triggered in precessing neutron stars, since the angular velocity vectors of the neutron and proton fluids are misaligned during precession. Our results suggest that the standard (Eulerian) slow precession that results for weak drag between the vortices and the charged fluid (protons and electrons) is not seriously affected by the instability. In contrast, the fast precession, which results when vortices are strongly coupled to the charged component, is generally unstable. The presence of this instability renders the standard (solid body) rotation model for free precession inconsistent and makes unsafe conclusions that have recently been drawn regarding neutron star interiors based on observations of precession in radio pulsars.     

Almost analytic models of ultramagnetized neutron star envelopes

Recent ROSAT measurements show that the X-ray emission from isolated neutron stars is modulated at the stellar rotation period. To interpret these measurements, one needs precise calculations of the heat transfer through the thin insulating envelopes of neutron stars. We present nearly analytic models of the thermal structure of the envelopes of ultramagnetized neutron stars. Specifically, we examine the limit in which only the ground Landau level is filled. We use the models to estimate the amplitude of modulation expected from non-uniformities in the surface temperatures of strongly magnetized neutron stars. In addition, we estimate cooling rates for stars with fields B  ∼ 10¹⁵ − 10¹⁶ G, which are relevant to models that invoke 'magnetars' to account for soft γ-ray emission from some repeating sources.     

Non-linear hydrodynamical evolution of rotating relativistic stars: numerical methods and code tests

We present numerical hydrodynamical evolutions of rapidly rotating relativistic stars, using an axisymmetric, non-linear relativistic hydrodynamics code. We use four different high-resolution shock-capturing (HRSC) finite-difference schemes (based on approximate Riemann solvers) and compare their accuracy in preserving uniformly rotating stationary initial configurations in long-term evolutions. Among these four schemes, we find that the third-order piecewise parabolic method scheme is superior in maintaining the initial rotation law in long-term evolutions, especially near the surface of the star. It is further shown that HRSC schemes are suitable for the evolution of perturbed neutron stars and for the accurate identification (via Fourier transforms) of normal modes of oscillation. This is demonstrated for radial and quadrupolar pulsations in the non-rotating limit, where we find good agreement with frequencies obtained with a linear perturbation code. The code can be used for studying small-amplitude or non-linear pulsations of differentially rotating neutron stars, while our present results serve as testbed computations for three-dimensional general-relativistic evolution codes.     

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.     

Protodiscs around hot magnetic rotator stars

We develop equations and obtain solutions for the structure and evolution of a protodisc region that is initially formed with no radial motion and super-Keplerian rotation speed when wind material from a hot rotating star is channelled towards its equatorial plane by a dipole-type magnetic field. Its temperature is around 10⁷ K because of shock heating and the inflow of wind material causes its equatorial density to increase with time. The centrifugal force and thermal pressure increase relative to the magnetic force and material escapes at its outer edge. The protodisc region of a uniformly rotating star has almost uniform rotation and will shrink radially unless some instability intervenes. In a star with angular velocity increasing along its surface towards the equator, the angular velocity of the protodisc region decreases radially outwards and magnetorotational instability (MRI) can occur within a few hours or days. Viscosity resulting from MRI will readjust the angular velocity distribution of the protodisc material and may assist in the formation of a quasi-steady disc. Thus, the centrifugal breakout found in numerical simulations for uniformly rotating stars does not imply that quasi-steady discs with slow outflow cannot form around magnetic rotator stars with solar-type differential rotation.     

The effect of r-mode instability on the evolution of isolated strange stars

We studied the evolution of isolated strange stars (SSs) synthetically, considering the influence of r -mode instability. Our results show that the cooling of SSs with non-ultrastrong magnetic fields is delayed by heating due to r -mode damping for millions of years, while the spin-down of the stars is dominated by gravitational radiation (GR). Especially for the SSs in a possible existing colour–flavour locked (CFL) phase, the effect of r -mode instability on the evolution of stars becomes extremely important because the viscosity, neutrino emissivity and specific heat involving pairing quarks are blocked. It leads to the cooling of these colour superconducting stars being very slow and the stars can remain at high temperature for millions of years, which differs completely from previous understanding. In this case, an SS in CFL phase can be located at the bottom of its r -mode instability window for a long time, but does not spin-down to a very low frequency for hours.     

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.     

Excitation and visibility of slow modes in rotating B-type stars

We use the traditional approximation to describe oscillations with frequencies comparable to the angular rotation rate. Validity of this approximation in application to main-sequence B stars is discussed. Numerical results regarding mode stability and visibility are presented for a model of the Be star HD 163868. For this object, Walker et al. detected a record number of mode frequencies using data from the small space telescope MOST . Our interpretation of these data differs from that of Walker et al. In particular, we interpret peaks in the lowest frequency range as retrograde g modes. We find instability in a large number of modes that remain undetectable because of unfavourable aspect and/or effect of cancellation. There is no clear preference to excitation of prograde modes.     

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.     

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.     

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.     

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.     

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.     

Differential rotation in F stars

Differential rotation can be detected in single line profiles of stars rotating more rapidly than about v sin i = 10km s-1 with the Fourier transform technique. This allows to search for differential rotation in large samples to look for correlations between differential rotation and other stellar parameters. I analyze the fraction of differentially rotating stars as a function of color, rotation, and activity in a large sample of F-type stars. Color and rotation exhibit a correlation with differential rotation in the sense that more stars are rotating differentially in the cooler, less rapidly rotating stars. Effects of rotation and color, however, cannot be disentangled in the underlying sample. No trend with activity is found. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A search for magnetic fields in cool sdB stars

Hot cluster horizontal branch (HB) stars and field subdwarf B (sdB) stars are core helium burning stars that exhibit abundance anomalies that are believed to be due to atomic diffusion. Diffusion can be effective in these stars because they are slowly rotating. In particular, the slow rotation of the hot HB stars (T_eff > 11000 K), which show abundance anomalies, contrasts with the fast rotation of the cool HB stars, where the observed abundances are consistent with those of red giants belonging to the same cluster. The reason why sdB stars and hot HB stars are rotating slowly is unknown. In order to assess the possible role of magnetic fields on abundances and rotation, we investigated the occurrence of such fields in sdB stars with T_eff < 30 000 K, whose temperatures overlap with those of the hot HB stars. We conclude that large‐scale organised magnetic fields of kG order are not generally present in these stars but at the achieved accuracy, the possibility that they have fields of a few hundred Gauss remains open. We report the marginal detection of such a field in SB 290; further observations are needed to confirm it (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Rotation and activity in the solar-type stars of NGC 2547

We present high-resolution spectroscopy of a sample of 24 solar-type stars in the young (15–40 Myr), open cluster NGC 2547. We use our spectra to confirm cluster membership in 23 of these stars, to determine projected equatorial velocities and chromospheric activity, and to search for the presence of accretion discs. We find examples of both fast ( v _e sin  i >50 km s⁻¹) and slow ( v _e sin  i <10 km s⁻¹) rotators, but no evidence for active accretion in any of the sample. The distribution of projected rotation velocities is indistinguishable from the slightly older IC 2391 and IC 2602 clusters, implying similar initial angular momentum distributions and circumstellar disc lifetimes. The presence of very slow rotators indicates either that long (10–40 Myr) disc lifetimes or internal differential rotation are needed, or that NGC 2547 (and IC 2391/2602) were born with more slowly rotating stars than are presently seen in even younger clusters and associations. The solar-type stars in NGC 2547 follow a similar rotation–activity relationship to that seen in older clusters. X-ray activity increases until a saturation level is reached for v _e sin  i >15–20 km s⁻¹. We are unable to explain why this saturation level, of log( L _x L _bol)≃−3.3, is a factor of 2 lower than in other clusters, but rule out anomalously slow rotation rates or uncertainties in X-ray flux calculations.