Compact Stars In General Relativity Research Articles

Exact solutions for compact stars in general relativity

Complete, exact, analytic solutions to the spherically symmetric Einstein equations have been constructed. Beginning with the Tolman VII solution we first “complete” the solution by computing an expression for the pressure and subsequently the equation of state. The Tolman VII solution is then extended to include non-zero mass densities at the fluid-vacuum interface. The method developed to construct these solutions is extended to include transverse pressures and charge. Furthermore the quadratic density profile of the Tolman VII type solutions is modified to include quartic density profiles that lead to expressions in terms of Jacobi elliptic integrals. It is shown that, depending on the outer boundary conditions placed on the mass density, the equations of state lead to either a set of polytropes or behaviour that mimics strange quark matter.

Gravitational wave asteroseismology limits from low density nuclear matter and perturbative QCD

We investigate the fundamental mode of non-radial oscillations of non-rotating compact stars in general relativity using a set of equations of state (EOS) connecting state-of-the-art calculations at low and high densities. Specifically, a low density model based on the chiral effective field theory (EFT) and high density results based on perturbative Quantum Chromodynamics (QCD) are matched through different interpolating polytropes fulfilling thermodynamic stability and subluminality of the speed of sound, together with the additional requirement that the equations of state support a two solar mass star. We employ three representative models (EOS I, II and III) presented in ref. [1] such that EOS I gives the minimum stellar radius, EOS II the maximum stellar mass, and EOS III the maximum stellar radius. Using this family of equations of state, we find that the frequency and the damping time of the f-mode are constrained within narrow quite model-independent windows. We also analyze some proposed empirical relations that describe the f-mode properties in terms of the average density and the compactness of the neutron star. We discuss the stringency of these constrains and the possible role of physical effects that cannot be encoded in a mere interpolation between low and high density EOSs.

Embedded class solutions compatible for physical compact stars in general relativity

We have explored a family of new solutions satisfying Einstein’s field equations and Karmarkar condition. We have assumed an anisotropic stress-tensor with no net electric charge. Interestingly, the new solutions yield zero values of all the physical quantities for all even integer $n > 0$ . However, for all $n >0$ ( $n \neq $ even numbers) they yield physically possible solutions. We have tuned the solution for neutron star Vela X-1 so that the solutions matches the observed mass and radius. For the same star we have extensively discussed the behavior of the solutions. The solutions yield a stiffer equation of state for larger values of n since the adiabatic index increases and speed of sound approaches the speed of light. It is also found that the solution is physically possible for Vela X-1 if $1.8 \leq n < 7$ (with $n\neq 2,4,6$ ). All the solutions for $n \geq 7$ violates the causality condition and all the solutions with $0 < n < 1.8$ lead to complex values of transverse sound speed $ v_{t}$ . The range of well-behaved n depends on the mass and radius of compact stars.

Modeling of charged anisotropic compact stars in general relativity

A charged compact star models have been determined for anisotropic fluid distribution. We have solved the Einstein's- Maxwell field equations to construct the charged compact star models by using radial pressure, metric function $e^{\lambda}$ and electric charge function. The generic charged anisotropic solution is verified by exploring different physical conditions like, causality condition, mass-radius relation and stability of the solution (via. adiabatic index, TOV equations and Herrera cracking concept). It is observed that the present charged anisotropic compact star is compatible with the star PSR 1937+21. However we also presented the EOS $\rho=f(p)$ for present charged compact star model.

Anisotropic stellar models admitting conformal motion

We address the problem of finding static and spherically symmetric anisotropic compact stars in general relativity that admit conformal motions. The study is framed in the language of f(R) gravity theory in order to expose opportunity for further study in the more general theory. Exact solutions of compact stars are found under the assumption that spherically symmetric spacetimes which admits conformal motion with matter distribution anisotropic in nature. In this work, two cases have been studied for the existence of such solutions: first, we consider the model given by f(R) = R and then f(R) = aR+b. Finally, specific characteristics and physical properties have been explored by analytically along with graphical representations for conformally symmetric compact stars in f(R) gravity.

Astrophysical applications of the post-Tolman-Oppenheimer-Volkoff formalism

The bulk properties of spherically symmetric stars in general relativity can be obtained by integrating the Tolman-Oppenheimer-Volkoff (TOV) equations. In previous work we developed a "post-TOV" formalism - inspired by parametrized post-Newtonian theory - which allows us to classify in a parametrized, phenomenological form all possible perturbative deviations from the structure of compact stars in general relativity that may be induced by modified gravity at second post-Newtonian order. In this paper we extend the formalism to deal with the stellar exterior, and we compute several potential astrophysical observables within the post-TOV formalism: the surface redshift $z_s$, the apparent radius $R_{\rm app}$, the Eddington luminosity at infinity $L_{\rm E}^\infty$ and the orbital frequencies. We show that, at leading order, all of these quantities depend on just two post-TOV parameters $\mu_1$ and $\chi$, and we discuss the possibility to measure (or set upper bounds on) these parameters.

Post-Tolman-Oppenheimer-Volkoff formalism for relativistic stars

Besides their astrophysical interest, compact stars also provide an arena for understanding the properties of theories of gravity that differ from Einstein's general relativity. Numerous studies have shown that different modified theories of gravity can modify the bulk properties (such as mass and radius) of neutron stars for given assumptions on the microphysics. What is not usually stressed though is the strong degeneracy in the predictions of these theories for the stellar mass and radius. Motivated by this observation, in this paper we take an alternative route and construct a stellar structure formalism which, without adhering to any particular theory of gravity, describes in a simple parametrized form the departure from compact stars in general relativity. This "post-Tolman-Oppenheimer-Volkoff (TOV)" formalism for spherical static stars is inspired by the well-known parametrized post-Newtonian theory, extended to second post-Newtonian order by adding suitable correction terms to the fully relativistic TOV equations. We show how neutron star properties are modified within our formalism, paying special attention to the effect of each correction term. We also show that the formalism is equivalent to general relativity with an "effective" (gravity-modified) equation of state.

I-Love-Q anisotropically: Universal relations for compact stars with scalar pressure anisotropy

Certain physical quantities that characterize neutron stars and quark stars (e.g. their mass, spin angular momentum and quadrupole moment) are interrelated in a way that is approximately insensitive to their internal structure. Such approximately universal relations are useful to break degeneracies in data analysis for future radio, X-ray and gravitational wave observations. Although the pressure inside compact stars is most likely nearly isotropic, certain scenarios have been put forth that suggest otherwise, for example due to phase transitions. We here investigate whether pressure anisotropy affects the approximate universal relations and whether it prevents their use in future observations. We achieve this by numerically constructing slowly-rotating and tidally-deformed, anisotropic, compact stars in General Relativity to third order in spin. We find that anisotropy affects the universal relations only weakly; the relations become less universal by a factor of 1.5-3 relative to the isotropic case, but remain approximately universal to 10%. We succeed in explaining this increase in variability as an increase in the eccentricity variation of isodensity contours, which provides further support for the emergent approximate symmetry explanation of universality. Anisotropy does not affect the universal relations to a sufficient level to prevent their use in gravitational wave astrophysics or in experimental relativity. We provide an explicit example of the latter in dynamical Chern-Simons gravity. The increase in variability of the universal relations due to pressure anisotropy could affect their use in future X-ray observations. Given expected observational uncertainties, however, the relations remain sufficiently universal for use in such observations if the anisotropic modifications to the moment of inertia and the quadrupole moment are less than 10% of their isotropic values.

Approximate analytic expressions for circular orbits around rapidly rotating compact stars

We calculate stationary configurations of rapidly rotating compact stars in general relativity, to study the properties of circular orbits of test particles in the equatorial plane. We search for simple, but precise, analytical formulae for the orbital frequency, specific angular momentum and binding energy of a test particle, valid for any equation of state and for any rotation frequency of the rigidly rotating compact star, up to the mass-shedding limit. Numerical calculations are performed using precise 2-D codes based on multi-domain spectral methods. Models of rigidly rotating neutron stars and the space-time outside them are calculated for several equations of state of dense matter. Calculations are also performed for quark stars consisting of self-bound quark matter. At the mass-shedding limit, the rotational frequency converges to a Schwarzschildian orbital frequency at the equator. We show that orbital frequency for any orbit outside equator is also approximated by a Schwarzschildian formula. Using a simple approximation for the frame-dragging term, we obtain approximate expressions for the specific angular momentum and specific energy on the corotating circular orbits in the equatorial plane of neutron star, which are valid down to the stellar equator. The formulae recover reference numerical values with typically 1% of accuracy for neutron stars with M > 0.5 M_sun. They are less precise for quark stars consisting of self-bound quark matter.

Nonadiabatic oscillations of compact stars in general relativity

We have developed a formalism to study nonadiabatic, nonradial oscillations of nonrotating compact stars in the frequency domain, including the effects of thermal diffusion in the framework of general relativistic perturbation theory. When a general equation of state depending on temperature is used, the perturbations of the fluid result in heat flux which is coupled with the space-time geometry through the Einstein field equations. Our results show that the frequency of the first pressure (p) and gravity (g) oscillation modes is significantly affected by thermal diffusion, while that of the fundamental (f) mode is basically unaltered due to the global nature of that oscillation. The damping time of the oscillations is generally much smaller than in the adiabatic case (more than 2 orders of magnitude for the p- and g-modes) reflecting the effect of thermal dissipation. Both the isothermal and adiabatic limits are recovered in our treatment and we study in more detail the intermediate regime. Our formalism finds its natural astrophysical application in the study of the oscillation properties of newly born neutron stars, neutron stars with a deconfined quark core phase, or strange stars which are all promising sources of gravitational waves with frequencies in the band of the first generation and advanced ground-based interferometric detectors.

The gyromagnetic ratio of rapidly rotating compact stars in general relativity

We numerically calculate equilibrium configurations of uniformly rotating and charged neutron stars, in the case of insulating material and neglecting the electromagnetic forces acting on the equilibrium of the fluid. This allows us to study the behaviour of the gyromagnetic ratio for those objects when varying rotation rate and equation of state for the matter. Under the assumption of low charge and incompressible fluid, we find that the gyromagnetic ratio is directly proportional to the compaction parameter M/R of the star, and very little dependent on its angular velocity. Nevertheless, it seems impossible to have g = 2 for these models with low charge-to-mass ratio, where matter consists of a perfect fluid and the collapse limit is never reached.

New numerical scheme to compute three-dimensional configurations of quasiequilibrium compact stars in general relativity: Application to synchronously rotating binary star systems

We have developed a new numerical scheme to obtain quasiequilibrium structures of nonaxisymmetric compact stars such as binary neutron star systems as well as the spacetime around those systems in general relativity. Concerning quasiequilibrium states of binary systems in general relativity, several investigations have been already carried out by assuming conformal flatness of the spatial part of the metric. However, the validity of the conformally flat treatment has not been fully analyzed except for axisymmetric configurations. Therefore it is desirable to solve quasiequilibrium states by developing totally different methods from the conformally flat scheme. In this paper we present a new numerical scheme to solve directly the Einstein equations for 3D configurations without assuming conformal flatness, although we make use of the simplified metric for the spacetime. This new formulation is the extension of the scheme which has been successfully applied for structures of axisymmetric rotating compact stars in general relativity. It is based on the integral representation of the Einstein equations by taking the boundary conditions at infinity into account. We have checked our numerical scheme by computing equilibrium sequences of binary polytropic star systems in Newtonian gravity and those of axisymmetric polytropic stars in general relativity. We have applied this numerical code to binary star systems in general relativity and have succeeded in obtaining several equilibrium sequences of synchronously rotating binary polytropes with the polytropic indices N = 0.0, 0.5 and 1.0.