Universality in quasi-normal modes of neutron stars

We study universality in gravitational waves emitted from non-rotating neutron stars characterized by different equations of state (EOSs). We find that the quasi-normal mode frequencies of such waves, including the w -modes and the f -mode, display similar universal scaling behaviours that hold for most EOSs. Such behaviours are shown to stem from the mathematical structure of the axial and the polar gravitational wave equations, and the fact that the mass distribution function can be approximated by a cubic–quintic polynomial in the radius. As a benchmark for other realistic neutron stars, a simple model of neutron stars is adopted here to reproduce the pulsation frequencies and the generic scaling behaviours mentioned above with good accuracy.     

Oscillations of rapidly rotating stratified neutron stars

We use time evolutions of the linear perturbation equations to study the oscillations of rapidly rotating neutrons stars. Our models account for the buoyancy due to composition gradients and we study, for the first time, the nature of the resultant g modes in a fast spinning star. We provide detailed comparisons of non-stratified and stratified models. This leads to an improved understanding of the relationship between the inertial modes of a non-stratified star and the g modes of a stratified system. In particular, we demonstrate that each g mode becomes rotation dominated, i.e. approaches a particular inertial mode, as the rotation rate of the star is increased. We also discuss issues relating to the gravitational wave driven instability of the various classes of oscillation modes.     

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.     

Influence of the r-mode instability on hypercritically accreting neutron stars

We have investigated the influence of the r-mode instability on hypercritically accreting neutron stars in close binary systems during their common envelope phases, based on the scenario proposed by Brown et al. On the one hand, neutron stars are heated by the accreted matter at the stellar surface, but on the other hand they are also cooled down by the neutrino radiation. At the same time, the accreted matter transports its angular momentum and mass to the star. We have studied the evolution of the stellar mass, temperature and rotational frequency.
The gravitational-wave-driven instability of the r-mode oscillation strongly suppresses spinning up of the star, the final rotational frequency of which is well below the mass-shedding limit, in fact typically as low as 10 per cent of that of the mass-shedding state. On a very short time-scale the rotational frequency tends to approach a certain constant value and saturates there, as long as the amount of accreted mass does not exceed a certain limit to collapse to a black hole. This implies that a similar mechanism of gravitational radiation to that in the so-called 'Wagoner star' may work in this process. The star is spun up by accretion until the angular momentum loss by gravitational radiation balances the accretion torque. The time-integrated dimensionless strain of the radiated gravitational wave may be large enough to be detectable by gravitational wave detectors such as LIGO II.     

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.     

Masses and luminosities of O‐ and B‐type stars and red supergiants

Massive stars are of interest as progenitors of supernovae, i.e. neutron stars and black holes, which can be sources of gravitational waves. Recent population synthesis models can predict neutron star and gravitational wave observations but deal with a fixed supernova rate or an assumed initial mass function for the population of massive stars. Here we investigate those massive stars, which are supernova progenitors, i.e. with O‐ and early B‐type stars, and also all supergiants within 3 kpc. We restrict our sample to those massive stars detected both in 2MASS and observed by Hipparcos, i.e. only those stars with parallax and precise photometry. To determine the luminosities we calculated the extinctions from published multi‐colour photometry, spectral types, luminosity class, all corrected for multiplicity and recently revised Hipparcos distances. We use luminosities and temperatures to estimate the masses and ages of these stars using different models from different authors. Having estimated the luminosities of all our stars within 3 kpc, in particular for all O‐ and early B‐type stars, we have determined the median and mean luminosities for all spectral types for luminosity classes I, III, and V. Our luminosity values for supergiants deviate from earlier results: Previous work generally overestimates distances and luminosities compared to our data, this is likely due to Hipparcos parallaxes (generally more accurate and larger than previous ground‐based data) and the fact that many massive stars have recently been resolved into multiples of lower masses and luminosities. From luminosities and effective temperatures we derived masses and ages using mass tracks and isochrones from different authors. From masses and ages we estimated lifetimes and derived a lower limit for the supernova rate of ≈20 events/Myr averaged over the next 10 Myr within 600 pc from the sun. These data are then used to search for areas in the sky with higher likelihood for a supernova or gravitational wave event (like OB associations) (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Biases for neutron star mass, radius and distance measurements from Eddington-limited X-ray bursts

Eddington-limited X-ray bursts from neutron stars can be used in conjunction with other spectroscopic observations to measure neutron star masses, radii and distances. In order to quantify some of the uncertainties in the determination of the Eddington limit, we analysed a large sample of photospheric radius-expansion thermonuclear bursts observed with the Rossi X-ray Timing Explorer . We identified the instant at which the expanded photosphere 'touches down' back on to the surface of the neutron star and compared the corresponding touchdown flux to the peak flux of each burst. We found that for the majority of sources, the ratio of these fluxes is smaller than ≃1.6, which is the maximum value expected from the changing gravitational redshift during the radius expansion episodes (for a  2 M_⊙  neutron star). The only sources for which this ratio is larger than ≃1.6 are high-inclination sources that include dippers and Cyg X-2. We discuss two possible geometric interpretations of this effect and show that the inferred masses and radii of neutron stars are not affected by this bias. On the other hand, systematic uncertainties as large as ∼50 per cent may be introduced to the distance determination.     

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 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.     

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.     

Self-similar pressure oscillations in neutron star envelopes as probes of neutron star structure

We study eigenmodes of acoustic oscillations of high multipolarity l ∼ 100–1000 and high frequency (∼100 kHz), localized in neutron star envelopes. We show that the oscillation problem is self-similar. Once the oscillation spectrum is calculated for a given equation of state (EOS) in the envelope and given stellar mass M and radius R , it can be rescaled to a star with any M and R (but the same EOS in the envelope). For l ≳ 300, the modes can be subdivided into the outer and inner ones. The outer modes are mainly localized in the outer envelope. The inner modes are mostly localized near the neutron drip point, being associated with the softening of the EOS after the neutron drip. We calculate oscillation spectra for the EOSs of cold-catalyzed and accreted matter and show that the spectra of the inner modes are essentially different. A detection and identification of high-frequency pressure modes would allow one to infer M and R and determine also the EOS in the envelope (accreted or ground state) providing a new, potentially powerful method to explore the main parameters and internal structure of neutron stars.     

Comparison of the large-scale clustering in the APM and the EDSGC galaxy surveys

We present the results of a numerical study of the fluid f, p and the gravitational w modes for increasingly relativistic non-rotating polytropes. The results for f and w modes are in good agreement with previous data for uniform density stars, which supports an understanding of the nature of the gravitational wave modes based on the uniform density data. We show that the p modes can become extremely long-lived for some relativistic stars. This effect is attributed to the change in the perturbed density distribution as the star becomes more compact.     

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.     

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.     

Simulating a stochastic background of gravitational waves from neutron star formation at cosmological distances

We develop a temporal simulation of the potentially detectable gravitational wave background from neutron star formation at cosmological distances. By using a recent model for the evolving star formation rate, we investigate the statistical distribution of gravitational wave amplitudes due to supernovae that result in neutron star formation in the Einstein–de Sitter cosmology. We find that the gravitational wave amplitude distribution in our frame is highly skewed, with skewness related to the distribution of sources, and that the potentially detectable gravitational wave strain is dominated by sources at a redshift of     . Time traces of the simulation, using selected waveforms, are presented graphically and are also made available as web-based audio files. The method developed can readily be extended to different cosmologies, as well as to incorporate other waveforms and source types. This type of simulation will be useful in testing and optimizing detection strategies for gravitational wave backgrounds due to various types of individually undetectable astrophysical sources.     

Newtonian hydrodynamics of the coalescence of black holes with neutron stars – II. Tidally locked binaries with a soft equation of state

We present a numerical study of the hydrodynamics in the final stages of inspiral of a black hole–neutron star binary, when the binary separation becomes comparable to the stellar radius. We use a Newtonian three-dimensional smooth particle hydrodynamics (SPH) code, and model the neutron star with a soft (adiabatic index Γ=5/3) polytropic equation of state, and the black hole as a Newtonian point mass that accretes matter via an absorbing boundary at the Schwarzschild radius. Our initial conditions correspond to tidally locked binaries in equilibrium, and we have explored configurations with different values of the mass ratio q M _NS M _BH, ranging from q =1 to 0.1. The dynamical evolution is followed for approximately 23 ms, and in every case studied here we find that the neutron star is tidally disrupted on a dynamical time-scale, forming a dense torus around the black hole that contains a few tenths of a solar mass. A nearly baryon-free axis is present in the system throughout the coalescence, and only modest beaming of a fireball that could give rise to a gamma-ray burst would be sufficient to avoid excessive baryon contamination. We find that some mass (of the order of 10⁻³–10⁻² M_⊙) may be dynamically ejected from the system, and could thus contribute substantially to the amount of observed r-process material in the galaxy. We calculate the gravitational radiation waveforms and luminosity emitted during the coalescence in the quadrupole approximation.     

The spin evolution of nascent neutron stars

The loss of angular momentum owing 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 magnetic field of the star. 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 time-scale of approximately 1 year. The results suggest that stars with magnetic fields ≤10¹³ 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.     

Estimating the effect of intermediate-range forces on the mass of rapidly-rotating neutron stars

Intermediate-range gravitational forces have been predicted by certain grand unified theories. If such forces exist, they would naturally affect the structure of neutron stars. Here, a simple rotating neutron star model is constructed which, under fairly mild assumptions, can be integrated exactly for the pressure. According to this model, the effect on neutron star masses by intermediate range forces is negligible, except when the range approaches the radius of the star and the coupling constant is close to the usual gravitation constant. In addition, extremely short range forces can be shown to have negligible effect, even when the coupling constant is many orders of magnitude greater thanG. Thus, there appears to be little hope of using neutron star mass measurements to test such grand unified theories.     

用中子星限制暗物质粒子散射截面

中子星可以通过重子物质和暗物质的相互作用吸积暗物质,且在一定条件下, 中子星吸积的暗物质粒子可以引发自引力塌缩形成小型黑洞, 生成的黑洞可能会进一步吞噬中子星.依据文献已有模型, 基于以上物理过程给出了在暗物质粒子不同质量下对暗物质粒子--中子的散射截面的限制.使用弱相互作用大质量粒子(Weakly Interacting Massive Particle, WIMP)模型, 并考虑暗物质粒子是玻色子的情形, 讨论了暗物质粒子有无自相互作用以及有无湮灭等条件下对限制暗物质参数的影响.既考虑了已发现的两个中子星系统来给出对暗物质参数空间的限制,也考虑了两个可能存在的年老中子星来预测未来观测可能对暗物质参数空间的限制.对于考虑玻色--爱因斯坦凝聚(Bose-Einstein Condensate, BEC)的玻色子暗物质, 在无自相互作用或有弱自相互作用, 无湮灭或有很小湮灭截面的条件下,中子星给出的间接观测对暗物质粒子-中子散射截面的限制的强度比XENON1T直接探测实验来得更强.未来, 如果在银心附近能观测到年老中子星, 其观测结果可以提升模型给出的对暗物质粒子--中子散射截面的限制, 从而帮助人们进一步理解暗物质.     

Chromothermal oscillations and collapse of strange stars to black holes: astrophysical implications

We study the effects of temperature on strange stars. It is found that the maximum mass of the star decreases with the increase of temperature, as at high temperatures the equations of state become softer. Moreover, if the temperature of a strange star increases, keeping its baryon number fixed, its gravitational mass increases and its radius decreases. This leads to a limiting temperature, where it turns into a black hole. These features are the result of a combined effect of the change of gluon mass and the quark distribution with temperature. We report on a new type of radial oscillation of strange stars, driven by what we call 'chromothermal' instability. We also discuss the relevance of our findings in the astrophysics of core collapse supernovae and gamma-ray bursts.