Numerical approach for high precision 3-D relativistic star models

A multi-domain spectral method for computing very high precision 3-D stellar models is presented. The boundary of each domain is chosen in order to coincide with a physical discontinuity (e.g. the star's surface). In addition, a regularization procedure is introduced to deal with the infinite derivatives on the boundary that may appear in the density field when stiff equations of state are used. Consequently all the physical fields are smooth functions on each domain and the spectral method is absolutely free of any Gibbs phenomenon, which yields to a very high precision. The power of this method is demonstrated by direct comparison with analytical solutions such as MacLaurin spheroids and Roche ellipsoids. The relative numerical error reveals to be of the order of $10^{-10}$. This approach has been developed for the study of relativistic inspiralling binaries. It may be applied to a wider class of astrophysical problems such as the study of relativistic rotating stars too.Comment: Minor changes, Phys. Rev. D in pres

MHD of rotating compact stars with spectral methods: description of the algorithm and tests

A flexible spectral code for the study of general relativistic magnetohydrodynamics is presented. Aiming at investigating the physics of slowly rotating magnetized compact stars, this new code makes use of various physically motivated approximations. Among them, the relativistic anelastic approximation is a key ingredient of the current version of the code. In this article, we mainly outline the method, putting emphasis on algorithmic techniques that enable to benefit as much as possible of the non-dissipative character of spectral methods, showing also a potential astrophysical application and providing a few illustrative tests.Comment: 15 pages, 4 figures (new figure added, misprints corrected) Article accepted for publication in a special issue of Classical and Quantum Gravity "New Frontiers in Numerical Relativity

Maximal mass of uniformly rotating homogeneous stars in Einsteinian gravity

Using a multi domain spectral method, we investigate systematically the general-relativistic model for axisymmetric uniformly rotating, homogeneous fluid bodies generalizing the analytically known Maclaurin and Schwarzschild solutions. Apart from the curves associated with these solutions and a further curve of configurations that rotate at the mass shedding limit, two more curves are found to border the corresponding two parameter set of solutions. One of them is a Newtonian lens shaped sequence bifurcating from the Maclaurin spheroid sequence, while the other one corresponds to highly relativistic bodies with an infinite central pressure. The properties of the configuration for which both the gravitational and the baryonic masses, moreover angular velocity, angular momentum as well as polar red shift obtain their maximal values are discussed in detail. In particular, by comparison with the static Schwarzschild solution, we obtain an increase of 34.25% in the gravitational mass. Moreover, we provide exemplarily a discussion of angular velocity and gravitational mass on the entire solution class.Comment: 4 pages, 4 figures, 1 table, submitted to A&A, corrected eq. for W, W' in 3.

Constraint violation in free evolution schemes: comparing BSSNOK with a conformal decomposition of Z4

We compare numerical evolutions performed with the BSSNOK formulation and a conformal decomposition of a Z4-like formulation of General Relativity. The important difference between the two formulations is that the Z4 formulation has a propagating Hamiltonian constraint, whereas BSSNOK has a zero-speed characteristic variable in the constraint subsystem. In spherical symmetry we evolve both puncture and neutron star initial data. We demonstrate that the propagating nature of the Z4 constraints leads to results that compare favorably with BSSNOK evolutions, especially when matter is present in the spacetime. From the point of view of implementation the new system is a simple modification of BSSNOK.Comment: Published in PR

Rotating star initial data for a constrained scheme in numerical relativity

A new numerical code for computing stationary axisymmetric rapidly rotating stars in general relativity is presented. The formulation is based on a fully constrained-evolution scheme for 3+1 numerical relativity using the Dirac gauge and maximal slicing. We use both the polytropic and MIT bag model equations of state to demonstrate that the code can construct rapidly rotating neutron star and strange star models. We compare numerical models obtained by our code and a well-established code, which uses a different gauge condition, and show that the two codes agree to high accuracy.Comment: Minor changes and one figure added. Version accepted for publication in Class. Quant. Gra

Keplerian frequency of uniformly rotating neutron stars and quark stars

We calculate Keplerian (mass shedding) configurations of rigidly rotating neutron stars and quark stars with crusts. We check the validity of empirical formula for Keplerian frequency, f_K, proposed by Lattimer & Prakash, f_K(M)=C (M/M_sun)^1/2 (R/10km)^-3/2, where M is the (gravitational) mass of Keplerian configuration, R is the (circumferential) radius of the non-rotating configuration of the same gravitational mass, and C = 1.04 kHz. Numerical calculations are performed using precise 2-D codes based on the multi-domain spectral methods. We use a representative set of equations of state (EOSs) of neutron stars and quark stars. We show that the empirical formula for f_K(M) holds within a few percent for neutron stars with realistic EOSs, provided 0.5 M_sun < M < 0.9 M_max,stat, where M_max,stat is the maximum allowable mass of non-rotating neutron stars for an EOS, and C=C_NS=1.08 kHz. Similar precision is obtained for quark stars with 0.5 M_sun < M < 0.9 M_max,stat. For maximal crust masses we obtain C_QS = 1.15 kHz, and the value of C_QS is not very sensitive to the crust mass. All our C's are significantly larger than the analytic value from the relativistic Roche model, C_Roche = 1.00 kHz. For 0.5 M_sun < M < 0.9 M_max,stat, the equatorial radius of Keplerian configuration of mass M, R_K(M), is, to a very good approximation, proportional to the radius of the non-rotating star of the same mass, R_K(M) = aR(M), with a_NS \approx a_QS \approx 1.44. The value of a_QS is very weakly dependent on the mass of the crust of the quark star. Both a's are smaller than the analytic value a_Roche = 1.5 from the relativistic Roche model.Comment: 6 pages, 6 color figures, submitted to A&

Constraining crystalline color superconducting quark matter with gravitational-wave data

We estimate the maximum equatorial ellipticity sustainable by compact stars composed of crystalline color-superconducting quark matter. For the theoretically allowed range of the gap parameter $\Delta$, the maximum ellipticity could be as large as $10^{-2}$, which is about 4 orders of magnitude larger than the tightest upper limit obtained by the recent science runs of the LIGO and GEO600 gravitational wave detectors based on the data from 78 radio pulsars. We point out that the current gravitational-wave strain upper limit already has some implications for the gap parameter. In particular, the upper limit for the Crab pulsar implies that $\Delta$ is less than O(20) MeV for a range of quark chemical potential accessible in compact stars, assuming that the pulsar has a mass $1.4 M_{\odot}$, radius 10 km, breaking strain $10^{-3}$, and that it has the maximum quadrupole deformation it can sustain without fracturing.Comment: Minor changes to match the published versio

The bar-mode instability in differentially rotating neutron stars: Simulations in full general relativity

We study the dynamical stability against bar-mode deformation of rapidly spinning neutron stars with differential rotation. We perform fully relativistic 3D simulations of compact stars with $M/R \geq 0.1$, where $M$ is the total gravitational mass and $R$ the equatorial circumferential radius. We adopt an adiabatic equation of state with adiabatic index $\Gamma=2$. As in Newtonian theory, we find that stars above a critical value of $\beta \equiv T/W$ (where $T$ is the rotational kinetic energy and $W$ the gravitational binding energy) are dynamically unstable to bar formation. For our adopted choices of stellar compaction and rotation profile, the critical value of $\beta = \beta_{dGR}$ is $\sim 0.24-0.25$, only slightly smaller than the well-known Newtonian value $\sim 0.27$ for incompressible Maclaurin spheroids. The critical value depends only very weakly on the degree of differential rotation for the moderate range we surveyed. All unstable stars form bars on a dynamical timescale. Models with sufficiently large $\beta$ subsequently form spiral arms and eject mass, driving the remnant to a dynamically stable state. Models with moderately large $\beta \gtrsim \beta_{dGR}$ do not develop spiral arms or eject mass but adjust to form dynamically stable ellipsoidal-like configurations. If the bar-mode instability is triggered in supernovae collapse or binary neutron star mergers, it could be a strong and observable source of gravitational waves. We determine characteristic wave amplitudes and frequencies.Comment: 17 pages, accepted for publication in AP

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 where the collapse limit is never reached.Comment: 11 pages, 6 figures, accepted for publication in Classical and Quantum Gravit