Towards gravitational wave asteroseismology

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 these 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 relations prove to be surprisingly robust, and we show how they can be used to extract the details of the star from observed modes. The results indicate that, should the various pulsation modes be detected by the new generation of gravitational wave detectors that come online in a few years, the mass and the radius of neutron stars can be deduced with errors no larger than a few per cent.     

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.     

Torsional oscillations of magnetized relativistic stars

Strong magnetic fields in relativistic stars can be a cause of crust fracturing, resulting in the excitation of global torsional oscillations. Such oscillations could become observable in gravitational waves or in high-energy radiation, thus becoming a tool for probing the equation of state of relativistic stars. As the eigenfrequency of torsional oscillation modes is affected by the presence of a strong magnetic field, we study torsional modes in magnetized relativistic stars. We derive the linearized perturbation equations that govern torsional oscillations coupled to the oscillations of a magnetic field, when variations in the metric are neglected (Cowling approximation). The oscillations are described by a single two-dimensional wave equation, which can be solved as a boundary-value problem to obtain eigenfrequencies. We find that, in the non-magnetized case, typical oscillation periods of the fundamental     torsional modes can be nearly a factor of 2 larger for relativistic stars than previously computed in the Newtonian limit. For magnetized stars, we show that the influence of the magnetic field is highly dependent on the assumed magnetic field configuration, and simple estimates obtained previously in the literature cannot be used for identifying normal modes observationally.     

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.     

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.     

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.     

Alfvén polar oscillations of relativistic stars

We study polar Alfvén oscillations of relativistic stars endowed with a strong global poloidal dipole magnetic field. Here, we focus only on the axisymmetric oscillations which are studied by numerically evolving the two-dimensional perturbation equations. Our study shows that the spectrum of the polar Alfvén oscillations is discrete in contrast to the spectrum of axial Alfvén oscillations which is continuous. We also show that the typical fluid modes, such as the f and p modes, are not significantly affected by the presence of the strong magnetic field.     

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.     

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.     

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.     

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.     

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.     

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.     

Radiation pressure supported stars in Einstein gravity: eternally collapsing objects

Even when we consider Newtonian stars, that is, stars with surface gravitational redshift   z ≪ 1  , it is well known that, theoretically, it is possible to have stars supported against self-gravity almost entirely by radiation pressure. However, such Newtonian stars must necessarily be supermassive. We point out that this requirement for excessively large M in the Newtonian case is a consequence of the occurrence of low   z ≪ 1  . However, if we remove such restrictions, and allow for the possible occurrence of a highly general relativistic regime,   z ≫ 1  , we show that it is possible to have radiation pressure supported stars (RPSSs) at an arbitrary value of M . Since RPSSs necessarily radiate at the Eddington limit, in Einstein gravity, they are never in strict hydrodynamical equilibrium. Further, it is believed that sufficiently massive or dense objects undergo continued gravitational collapse to the black hole (BH) stage characterized by   z =∞  . Thus, late stages of BH formation, by definition, would have   z ≫ 1  , and hence would be examples of quasi-stable general relativistic RPSSs. It is shown that the observed duration of such Eddington limited radiation pressure dominated states is   t ≈ 5 × 10⁸ (1 + z ) yr  . Thus,   t →∞  as BH formation  ( z →∞)  takes place. Consequently, such radiation pressure dominated extreme general relativistic stars become eternally collapsing objects (ECOs) and the BH state is preceded by such an ECO phase. This result is also supported by our previous finding that trapped surfaces are not formed in gravitational collapse and the value of the integration constant in the vacuum Schwarzschild solution is zero. Hence the supposed observed BHs are actually ECOs.     

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.     

Monopole gravitational waves from relativistic fireballs driving gamma-ray bursts

Einstein's general relativity predicts that pressure, in general stresses, plays a similar role to energy density,  ε=ρ c ²  (with ρ being the corresponding mass density), in generating gravity. The source of gravitational field, the active gravitational mass density, sometimes referred to as Whittaker's mass density, is  ρ_grav=ρ+ 3 p / c ²  , where p is pressure in the case of an ideal fluid. Whittaker's mass is not conserved, hence its changes can propagate as monopole gravitational waves. Such waves can be generated only by astrophysical sources with varying gravitational mass. Here we show that relativistic fireballs, considered in modelling gamma-ray burst phenomena, are likely to radiate monopole gravitational waves from high-pressure plasma with varying Whittaker's mass. Also, ejection of a significant amount of initial mass-energy of the progenitor contributes to the monopole gravitational radiation. We identify monopole waves with   h ¹¹+ h ²²  waves of Eddington's classification which propagate (in the z -direction) together with the energy carried by massless fields. We show that the monopole waves satisfy Einstein's equations, with a common stress-energy tensor for massless fields. The polarization mode of monopole waves is  Φ₂₂  , i.e. these are perpendicular waves which induce changes of the radius of a circle of test particles only (breathing mode). The astrophysical importance of monopole gravitational waves is discussed.     

General relativistic self-similar solutions in cosmology

We present general relativistic solutions for self-similar spherical perturbations in an expanding cosmological background of cold pressure-less gas. We focus on solutions having shock discontinuities propagating in the surrounding cold gas. The pressure, p , and energy density, μ, in the shock-heated matter are assumed to obey   p = w μ  , where w is a positive constant. Consistent solutions are found for shocks propagating from the symmetry centre of a region of a positive density excess over the background. In these solutions, shocks exist outside the radius marking the event horizon of the black hole which would be present in a shock-less collapse. For large jumps in the energy density at the shock, a black hole is avoided altogether and the solutions are regular at the centre. The shock-heated gas does not contain any sonic points, provided the motion of the cold gas ahead of the shock deviates significantly from the Hubble flow. For shocks propagating in the uniform background, sonic points always appear for small jumps in the energy density. We also discuss self-similar solutions without shocks in fluids with   w < −1/3  .     

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)     

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.     

Possible Alternatives to the Supermassive Black Hole at the Galactic Center

Now there are two basic observational techniques to investigate a gravitational potential at the Galactic Center, namely, (a) monitoring the orbits of bright stars near the Galactic Center to reconstruct a gravitational potential; (b) measuring the size and shape of shadows around black hole giving an alternative possibility to evaluate black hole parameters in mm-band with VLBI-technique. At the moment, one can use a small relativistic correction approach for stellar orbit analysis (however, in the future the approximation will not be precise enough due to enormous progress of observational facilities) while for smallest structure analysis in VLBI observations one really needs a strong gravitational field approximation. We discuss results of observations, their conventional interpretations, tensions between observations and models and possible hints for a new physics from the observational data and tensions between observations and interpretations. We discuss an opportunity to use a Schwarzschild metric for data interpretation or we have to use more exotic models such as Reissner–Nordstrom or Schwarzschild–de-Sitter metrics for better fits.