Does Pressure Increase or Decrease Active Gravitational Mass Density?

It is known that, for a static fluid sphere, the GeneralRelativistic (GR) Effective Mass Energy Density (EMD) appears to be (rho + 3 p), where rho is the bare mass density, p is the isotropic pressure, from a purely localized view point. But since there is no truly local definition of ``gravitational field'', such a notion could actually be misleading. On the other hand, by using the Tolman mass formula, we point out that, from a global perspective, the Active Gravitational Mass Energy Density (AGMD) is sqrt{g_{00}} (rho + 3 p) and which is obviously smaller than (rho + 3p) because g_{00} < 1. Then we show that the AGMD eventually is (rho - 3p), i.e., exactly opposite to what is generally believed. We further identify the AGMD to be proportional to the Ricci Scalar. By using this fundamental and intersting property, we obtain the GR virial theorem in terms of appropriate ``proper energies''.Comment: The originally accepted journal version, subsequent modifications remove

General Relativistic Stars : Polytropic Equations of State

In this paper, the gravitational field equations for static spherically symmetric perfect fluid models with a polytropic equation of state, $p=k\rho^{1+1/n}$, are recast into two complementary 3-dimensional {\it regular} systems of ordinary differential equations on compact state spaces. The systems are analyzed numerically and qualitatively, using the theory of dynamical systems. Certain key solutions are shown to form building blocks which, to a large extent, determine the remaining solution structure. In one formulation, there exists a monotone function that forces the general relativistic solutions towards a part of the boundary of the state space that corresponds to the low pressure limit. The solutions on this boundary describe Newtonian models and thus the relationship to the Newtonian solution space is clearly displayed. It is numerically demonstrated that general relativistic models have finite radii when the polytropic index $n$ satisfies $0\leq n \lesssim 3.339$ and infinite radii when $n\geq 5$. When $3.339\lesssim n<5$, there exists a 1-parameter set of models with finite radii and a finite number, depending on $n$, with infinite radii.Comment: 31 pages, 10 figure

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.Comment: 6 pages, 6 figures, A&A in pres

The Stability of an Isentropic Model for a Gaseous Relativistic Star

We show that the isentropic subclass of Buchdahl's exact solution for a gaseous relativistic star is stable and gravitationally bound for all values of the compactness ratio $u [\equiv (M/R)$, where $M$ is the total mass and $R$ is the radius of the configuration in geometrized units] in the range, $0 < u \leq 0.20$, corresponding to the {\em regular} behaviour of the solution. This result is in agreement with the expectation and opposite to the earlier claim found in the literature.Comment: 9 pages (including 1 table); accepted for publication in GR

Gaussian integration with rescaling of abscissas and weights

An algorithm for integration of polynomial functions with variable weight is considered. It provides extension of the Gaussian integration, with appropriate scaling of the abscissas and weights. Method is a good alternative to usually adopted interval splitting.Comment: 14 pages, 5 figure

Generalized Fermi-Dirac Functions and Derivatives: Properties and Evaluation

The generalized Fermi-Dirac functions and their derivatives are important in evaluating the thermodynamic quantities of partially degenerate electrons in hot dense stellar plasmas. New recursion relations of the generalized Fermi-Dirac functions have been found. An effective numerical method to evaluate the derivatives of the generalized Fermi-Dirac functions up to third order with respect to both degeneracy and temperature is then proposed, following Aparicio. A Fortran program based on this method, together with a sample test case, is provided. Accuracy and domain of reliability of some other, popularly used analytic approximations of the generalized Fermi-Dirac functions for extreme conditions are investigated and compared with our results.Comment: accepted for publication in Comp. Phys. Com

Truncated post-Newtonian neutron star model

As a preliminary step towards simulating binary neutron star coalescing problem, we test a post-Newtonian approach by constructing a single neutron star model. We expand the Tolman-Oppenheimer-Volkov equation of hydrostatic equilibrium by the power of $c^{-2}$, where $c$ is the speed of light, and truncate at the various order. We solve the system using the polytropic equation of state with index $\Gamma=5/3, 2$ and 3, and show how this approximation converges together with mass-radius relations. Next, we solve the Hamiltonian constraint equation with these density profiles as trial functions, and examine the differences in the final metric. We conclude the second `post-Newtonian' approximation is close enough to describe general relativistic single star. The result of this report will be useful for further binary studies. (Note to readers) This paper was accepted for publication in Physical Review D. [access code dsj637]. However, since I was strongly suggested that the contents of this paper should be included as a section in our group's future paper, I gave up the publication.Comment: 5 pages, RevTeX, 3 eps figs, epsf.sty, accepted for publication in PRD (Brief Report), but will not appea

Excitation of the odd-parity quasi-normal modes of compact objects

The gravitational radiation generated by a particle in a close unbounded orbit around a neutron star is computed as a means to study the importance of the $w$ modes of the neutron star. For simplicity, attention is restricted to odd parity (``axial'') modes which do not couple to the neutron star's fluid modes. We find that for realistic neutron star models, particles in unbounded orbits only weakly excite the $w$ modes; we conjecture that this is also the case for astrophysically interesting sources of neutron star perturbations. We also find that for cases in which there is significant excitation of quadrupole $w$ modes, there is comparable excitation of higher multipole modes.Comment: 18 pages, 21 figures, submitted to Phys. Rev.

Post-Newtonian Models of Binary Neutron Stars

Using an energy variational method, we calculate quasi-equilibrium configurations of binary neutron stars modeled as compressible triaxial ellipsoids obeying a polytropic equation of state. Our energy functional includes terms both for the internal hydrodynamics of the stars and for the external orbital motion. We add the leading post-Newtonian (PN) corrections to the internal and gravitational energies of the stars, and adopt hybrid orbital terms which are fully relativistic in the test-mass limit and always accurate to PN order. The total energy functional is varied to find quasi-equilibrium sequences for both corotating and irrotational binaries in circular orbits. We examine how the orbital frequency at the innermost stable circular orbit depends on the polytropic index n and the compactness parameter GM/Rc^2. We find that, for a given GM/Rc^2, the innermost stable circular orbit along an irrotational sequence is about 17% larger than the innermost secularly stable circular orbit along the corotating sequence when n=0.5, and 20% larger when n=1. We also examine the dependence of the maximum neutron star mass on the orbital frequency and find that, if PN tidal effects can be neglected, the maximum equilibrium mass increases as the orbital separation decreases.Comment: 53 pages, LaTex, 9 figures as 10 postscript files, accepted by Phys. Rev. D, replaced version contains updated reference

Dynamics of a thin shell in the Reissner-Nordstrom metric

We describe the dynamics of a thin spherically symmetric gravitating shell in the Reissner-Nordstrom metric of the electrically charged black hole. The energy-momentum tensor of electrically neutral shell is modelled by the perfect fluid with a polytropic equation of state. The motion of a shell is described fully analytically in the particular case of the dust equation of state. We construct the Carter-Penrose diagrams for the global geometry of the eternal black hole, which illustrate all possible types of solutions for moving shell. It is shown that for some specific range of initial parameters there are possible the stable oscillating motion of the shell transferring it consecutively in infinite series of internal universes. We demonstrate also that this oscillating type of motion is possible for an arbitrary polytropic equation of state on the shell.Comment: 17 pages, 7 figure