Keplerian Limit Research Articles

Rapidly rotating neutron stars in f(R,T)=R+2 lambda T gravity

In this work, we study the influence of f(R, T) gravity on rapidly rotating neutron stars. First we discuss the main aspects of this modified theory of gravity where the gravitational Lagrangian is an arbitrary function of the Ricci scalar R and of the trace of the energy–momentum tensor T. Then we present the basic equations for neutron stars including the equations of state used in the present work to describe the hadronic matter. Some physical quantities of interest are calculated such as mass–radius relations, moments of inertia, angular momentum, and compactness. By considering four different rotation regimes, we obtain results that indicate substantial modifications in the physical properties of neutron stars in f(R, T) gravity when compared to those in the context of general relativity. In particular, the mass–radius relation for sequences of stars indicates that f(R, T) gravity increases the mass and the equatorial radius of the neutron stars for stars rotating with an angular velocity smaller than Kepler limit.

Revised formalism for slowly rotating superfluid neutron stars in general relativity

We discuss slowly-rotating, general relativistic, superfluid neutron stars in the Hartle-Thorne formulation. The composition of the stars is described by a simple two-fluid model which accounts for superfluid neutrons and all other constituents. We apply a perturbed matching framework to derive a new formalism for slowly-rotating superfluid neutron stars, valid up to second-order perturbation theory, building on the original formulation reported by Andersson and Comer in 2001. The present study constitutes an extension of previous work in the single-fluid case where it was shown that the Hartle-Thorne formalism needs to be amended since it does not provide the correct results when the energy density does not vanish at the surface of the star. We discuss in detail the corrections that need to be applied to the original two-fluid formalism in order to account for non vanishing energy densities at the boundary. In the process, we also find a correction needed in the computation of the deformation of the stellar surface in the original two-fluid model in all cases (irrespective of the value of the energy density at the surface). The discrepancies found between the two formalisms are illustrated by building numerical stellar models, focusing on the comparison in the calculation of the stellar mass, the deformation of the star, and in the Kepler limit of rotation. In particular, using a toy-model equation of state for which the energy density does not vanish at the boundary of the star we demonstrate that the corrections to the formalism we find impact the structure of slowly-rotating superfluid neutron stars in a significant way.

Fast Rotating Neutron Stars: Oscillations and Instabilities

In this review article, we present the main results from our most recent research concerning the oscillations of fast rotating neutron stars. We derive a set of time evolution equations for the investigation of non-axisymmetric oscillations of rapidly rotating compact objects in full general relativity, taking into account the contribution of a dynamic spacetime. Using our code, which features high accuracy at comparably low computational expense, we are able to extract the frequencies of non-axisymmetric modes of compact objects with rotation rates up to the Kepler limit. We propose various universal relations combining bulk properties of isolated neutron stars as well as of binary systems before and after merger; these relations are independent of the true equation of state and may serve as a valuable tool for gravitational wave asteroseismology. We also present an introductory example using a Bayesian analysis.

Synthetic Spectra of Rotating Stars

Abstract Many early-type stars have oblate surfaces, spatial temperature variations, and spectral line broadening that indicate large rotational velocities. Rotation ought to have a significant effect on the full spectra of such stars. To infer structural and life history parameters from their spectra, one must integrate specific intensity over the two-dimensional surfaces of corresponding stellar models. Toward this end, we offer Paint the Atmospheres of Rotating Stars—an integration scheme based on models that incorporate solid body rotation, Roche mass distribution, and collinearity of gravity and energy flux. The scheme features a closed-form expression for the azimuthal integral, a high-order numerical approximation of the longitudinal integral, and a precise calculation of surface effective temperature at rotation rates up to 99.9% of the Keplerian limit. Extensions of the scheme include synthetic color–magnitude diagrams and planetary transit curves.

Dynamics of fast rotating neutron stars: An approach in the Hilbert gauge

We describe a set of time evolution equations and its numerical implementation for the investigation of non-axisymmetric oscillations of rapidly rotating compact objects in full General Relativity, taking into account the contribution of a dynamic spacetime. We derive the perturbation equations for the spacetime in the Hilbert gauge, while the hydrodynamical evolution is based on perturbations of the energy-momentum tensor. In our numerical implementation, we use Kreiss-Oliger dissipation in order to achieve a stable time evolution. Our code features high accuracy at comparably low computational expense and we are able to extract the frequencies of non-axisymmetric modes of compact objects with rotation rates up to the Kepler limit.

Effect of spin on the inspiral of binary neutron stars

We perform long-term simulations of spinning binary neutron stars, with our highest dimensionless spin being $\chi \sim 0.32$. To assess the importance of spin during the inspiral we vary the spin, and also use two equations of state, one that consists of plain nuclear matter and produces compact stars (SLy), and a hybrid one that contains both nuclear and quark matter and leads to larger stars (ALF2). Using high resolution that has grid spacing $\Delta x\sim 98$ m on the finest refinement level, we find that the effects of spin in the phase evolution of a binary system can be larger than the one that comes from tidal forces. Our calculations demonstrate explicitly that although tidal effects are dominant for small spins ($\lesssim 0.1$), this is no longer true when the spins are larger, but still much smaller than the Keplerian limit.

On gravitational echoes from ultracompact exotic stars

At the dawn of a golden age for gravitational wave astronomy, we must leave no stone unturned in our quest for new phenomena beyond our current understanding of General Relativity (GR), particle physics and nuclear physics. In this paper we discuss gravitational echoes from ultracompact stars. We restrict our analysis to exact solutions of Einstein field equations in GR that are supported by an isotropic fluid with a physically motivated equations of state (EoS), and in particular we impose the constraint of causality. Our main conclusion is that ultracompact objects belonging to this category are not able to generate gravitational echoes like those that characterize the relaxation phase of a putative black hole mimicker. Nevertheless, we identify a class of physical exotic objects that are compact enough to accommodate the presence of an external unstable light ring, thus opening the possibility of trapping gravitational radiation and affecting the ringdown phase of a merger event. Most importantly, we show that once rotation is included these stars—contrary to what usually expected for ultracompact objects—are not plagued by any ergoregion instability. We extend our analysis for arbitrary values of angular velocity up to the Keplerian limit, and we comment about potential signals relevant for gravitational wave interferometers.

Oscillation modes of rapidly rotating neutron stars in scalar-tensor theories of gravity

We perform the first study of the oscillation frequencies of rapidly rotating neutron stars in alternative theories of gravity, focusing mainly on the fundamental $f$-modes. We concentrated on a particular class of alternative theories - the (massive) scalar-tensor theories. The generalization to rapid rotation is important because on one hand the rapid rotation can magnify the deviations from general relativity compared to the static case and on the other hand some of the most efficient emitters of gravitational radiation, such as the binary neutron star merger remnants, are supposed to be rotating close to their Kepler (mass-shedding) limits shortly after their formation. We have constructed several sequences of models starting from the nonrotating case and reaching up to the Kepler limit, with different values of the scalar-tensor theory coupling constant and the scalar field mass. The results show that the deviations from pure Einstein's theory can be significant especially in the case of nonzero scalar field mass. An important property of the oscillation modes of rapidly rotating stars is that they can become secularly unstable due to the emission of gravitational radiation, that is so-called Chandrasekhar-Friedman-Schutz instability. Such unstable modes are efficient emitters of gravitational radiation. Our studies show that the inclusion of nonzero scalar field would decrease the threshold value of the normalized angular momentum where this instability stars to operate, but the growth time of the instability seems to be increased compared to pure general relativity.

Rotating hybrid stars with the Dyson-Schwinger quark model

We study rapidly rotating hybrid stars with the Dyson-Schwinger model for quark matter and the Brueckner-Hartree-Fock many-body theory with realistic two-body and three-body forces for nuclear matter. We determine the maximum gravitational mass, equatorial radius, and rotation frequency of stable stellar configurations by considering the constraints of the Keplerian limit and the secular axisymmetric instability, and compare with observational data. We also discuss the rotational evolution for constant baryonic mass, and find a spinup phenomenon for supramassive stars before they collapse to black holes.

The stability of tidal equilibrium for hierarchical star–planet–moon systems

Motivated by the current search for exomoons, this paper considers the stability of tidal equilibrium for hierarchical three-body systems containing a star, a planet, and a moon. In this treatment, the energy and angular momentum budgets include contributions from the planetary orbit, lunar orbit, stellar spin, planetary spin, and lunar spin. The goal is to determine the optimized energy state of the system subject to the constraint of constant angular momentum. Due to the lack of a closed form solution for the full three-body problem, however, we must use use an approximate description of the orbits. We first consider the Keplerian limit and find that the critical energy states are saddle points, rather than minima, so that these hierarchical systems have no stable tidal equilibrium states. We then generalize the calculation so that the lunar orbit is described by a time-averaged version of the circular restricted three-body problem. In this latter case, the critical energy state is a shallow minimum, so that a tidal equilibrium state exists. In both cases, however, the lunar orbit for the critical point lies outside the boundary (roughly half the Hill radius) where (previous) numerical simulations indicate dynamical instability. These results suggest that star-planet-moon systems have no viable long-term stable states analogous to those found for two-body systems.

Rapidly rotating neutron stars in dilatonic Einstein-Gauss-Bonnet theory

We construct sequences of rapidly rotating neutron stars in dilatonic Einstein-Gauss-Bonnet theory, employing two equations of state for the nuclear matter. We analyze the dependence of the physical properties of these neutron stars on the Gauss-Bonnet coupling strength. For a given equation of state we determine the physically relevant domain of rapidly rotating neutron stars, which is delimited by the set of neutron stars rotating at the Kepler limit, the set of neutron stars along the secular instability line, and the set of static neutron stars. As compared to Einstein gravity, the presence of the Gauss-Bonnet term decreases this domain, leading to lower values for the maximum mass as well as to smaller central densities. The quadrupole moment is decreased by the Gauss-Bonnet term for rapidly rotating neutron stars, while it is increased for slowly rotating neutron stars. The universal relation between the quadrupole moment and the moment of inertia found in General Relativity appears to extend to dilatonic Einstein-Gauss-Bonnet theory with very little dependence on the coupling strength of the Gauss-Bonnet term. The neutron stars carry a small dilaton charge.

THE MOST LUMINOUS SUPERNOVA ASASSN-15LH: SIGNATURE OF A NEWBORN RAPIDLY ROTATING STRANGE QUARK STAR

ABSTRACT In this paper we show that the most luminous supernova discovered very recently, ASASSN-15lh, could have been powered by a newborn ultra-strongly magnetized pulsar, which initially rotates near the Kepler limit. We find that if this pulsar is a neutron star, its rotational energy could be quickly lost as a result of gravitational-radiation-driven r-mode instability; if it is a strange quark star (SQS), however, this instability is highly suppressed due to a large bulk viscosity associated with the nonleptonic weak interaction among quarks and thus most of its rotational energy could be extracted to drive ASASSN-15lh. Therefore, we conclude that such an ultra-energetic supernova provides a possible signature for the birth of an SQS.

Slowly rotating superfluid neutron stars with isospin dependent entrainment in a two-fluid model

We investigate the slowly rotating general relativistic superfluid neutron stars including the entrainment effect in a two-fluid model, where one fluid represents the superfluid neutrons and the other is the charge-neutral fluid called the proton fluid, made of protons and electrons. The equation of state and the entrainment effect between the superfluid neutrons and the proton fluid are computed using a relativistic mean field (RMF) model where baryon-baryon interaction is mediated by the exchange of $\sigma$, $\omega$, and $\rho$ mesons and scalar self interactions are also included. The equations governing rotating neutron stars in the slow rotation approximation are second order in rotational velocities of neutron and proton fluids. We explore the effects of the isospin dependent entrainment and the relative rotation between two fluids on the global properties of rotating superfluid neutron stars such as mass, shape, and the mass shedding (Kepler) limit within the RMF model with different parameter sets. It is observed that for the global properties of rotating superfluid neutron stars in particular, the Kepler limit is modified compared with the case that does not include the contribution of $\rho$ mesons in the entrainment effect.

GRAVITATIONAL WAVE EMISSION AND SPIN-DOWN OF YOUNG PULSARS

The rotation frequencies of young pulsars are systematically below their theoretical Kepler limit. R-modes have been suggested as a possible explanation for this observation. With the help of semi-analytic expressions that make it possible to assess the uncertainties of the r-mode scenario due to the impact of uncertainties in underlying microphysics, we perform a quantitative analysis of the spin-down and the emitted gravitational waves of young pulsars. We find that the frequency to which r-modes spin down a young neutron star is surprisingly insensitive both to the microscopic details and the saturation amplitude. Comparing our result to astrophysical data, we show that for a range of sufficiently large saturation amplitudes r-modes provide a viable spindown scenario and that all observed young pulsars are very likely already outside the r-mode instability region. Therefore the most promising sources for gravitational wave detection are unobserved neutron stars associated with recent supernovae, and we find that advanced LIGO should be able to see several of them. We find the remarkable result that the gravitational wave strain amplitude is completely independent of both the r-mode saturation amplitude and the microphysics, and depends on the saturation mechanism only within some tens of per cent. However, the gravitational wave frequency depends on the amplitude and we provide the required expected timing parameter ranges to look for promising sources in future searches.

Phase Diagram and Structural Diversity of the Densest Binary Sphere Packings

The densest binary sphere packings have historically been very difficult to determine. The only rigorously known packings in the α-x plane of sphere radius ratio α and relative concentration x are at the Kepler limit α=1, where packings are monodisperse. Utilizing an implementation of the Torquato-Jiao sphere-packing algorithm [S. Torquato and Y. Jiao, Phys. Rev. E 82, 061302 (2010)], we present the most comprehensive determination to date of the phase diagram in (α,x) for the densest binary sphere packings. Unexpectedly, we find many distinct new densest packings.

Alternative derivation of the relativistic contribution to perihelic precession

An alternative derivation of the first-order relativistic contribution to perihelic precession is presented. Orbital motion in the Schwarzschild geometry is considered in the Keplerian limit, and the orbit equation is derived for approximate elliptical motion. The method of solution makes use of coordinate transformations and the correspondence principle rather than the standard perturbative approach. The form of the resulting orbit equation is similar to that derived from Newtonian mechanics and includes first-order corrections to Kepler’s orbits due to general relativity. The associated relativistic contribution to perihelic precession agrees with established first-order results. The reduced radius for the circular orbit is in agreement to first-order with that calculated from the Schwarzschild effective potential. The method of solution is understandable by undergraduate students.

Relativistic numerical models for stationary superfluid neutron stars

We have developed a theoretical model and a numerical code for stationary rotating superfluid neutron stars in full general relativity. The underlying two-fluid model is based on Carter's covariant multifluid hydrodynamic formalism. The two fluids, representing the superfluid neutrons on one hand, and the protons and electrons on the other, are restricted to uniform rotation around a common axis, but are allowed to have different rotation rates. We have performed extensive tests of the numerical code, including quantitative comparisons to previous approximative results for these models. The results presented here are the first ``exact'' calculations of such models in the sense that no approximations (other than that inherent in a discretized numerical treatment) are used. Using this code we reconfirm the existence of prolate-oblate shaped configurations. We studied the dependency of the Kepler rotation limit and of the mass-density relation on the relative rotation rate. We further demonstrate how one can simulate a (albeit fluid) neutron-star crust by letting one fluid extend further outwards than the other, which results in interesting cases where the Kepler limit is actually determined by the outermost but slower fluid.

Statistical mechanics and phase diagrams of rotating self-gravitating fermions

We compute statistical equilibrium states of rotating self-gravitating systems enclosed within a box by maximizing the Fermi-Dirac entropy at fixed mass, energy and angular momentum. We increase the rotation up to the Keplerian limit and describe the flattening of the configuration until mass shedding occurs. At the maximum rotation, the system develops a cusp at the equator. We draw the equilibrium phase diagram of the rotating self-gravitating Fermi gas and discuss the structure of the caloric curve as a function of degeneracy parameter and angular velocity. We argue that systems described by the Fermi-Dirac distribution in phase space do not bifurcate to non-axisymmetric structures, in continuity with the case of polytropes with index n>0.808 (the Fermi gas at T=0 corresponds to n=3/2). This contrasts with the study of Votyakov et al. (2002) who consider a Fermi-Dirac distribution in configuration space and find ``double star'' structures (their model at T=0 corresponds to n=0). We also discuss the influence of rotation on the onset of the gravothermal catastrophe for globular clusters. On general grounds, we complete previous investigations concerning the nature of phase transitions in self-gravitating systems. We emphasize the inequivalence of statistical ensembles regarding the formation of binaries (or low-mass condensates) in the microcanonical ensemble and Dirac peaks (or massive condensates) in the canonical ensemble. We also describe an hysteretic cycle between the gaseous phase and the condensed phase that are connected by a ``collapse'' or an ``explosion''. This notion of hysteresis in self-gravitating systems is new.

Slowly rotating superfluid Newtonian neutron star model with entrainment

We develop a formalism that can be used to model slowly rotating superfluid Newtonian neutron stars. A simple two-fluid model is used to describe the matter, where one fluid consists of the superfluid neutrons that are believed to exist in the inner crust and core of mature neutron stars, while the other fluid is a charge neutral conglomerate of the remaining constituents (crust nuclei, core superconducting protons, electrons, etc). We include the entrainment effect, which is a non-dissipative interaction between the two fluids whereby a momentum induced in one of the fluids will cause part of the mass of the other fluid to be carried along. The equations that describe rotational equilibria (i.e. axisymmetric and stationary configurations) are approximated using the slow-rotation approximation; an expansion in terms of the rotation rates of the two fluids where only terms up to second-order are kept. Our formalism allows the neutrons to rotate at a rate different from that of the charged constituents. For a particular equation of state that is quadratic in the two mass-densities and relative velocities of the fluids, we find an analytic solution to the slow-rotation equations. This result provides an elegant generalisation to the two-fluid problem of the standard expressions for the one-fluid polytrope E / ρ 2 . The model equation of state includes entrainment and is general enough to allow for realistic values for, say, mass and radius of the star. It also includes a mixed term in the mass densities that can be related to symmetry energy terms that appear in more realistic equations of state. We use the analytic solution to explore how relative rotation between the two fluids, the symmetry energy term, and entrainment affect the neutron star's local distribution of particles, as well as global quantities as the Kepler limit, ellipticity, and moments of inertia.

Quadrupole Moments of Rotating Neutron Stars

Numerical models of rotating neutron stars are constructed for four equations of state using the computer code RNS written by Stergioulas. For five selected values of the star's gravitational mass (in the interval between 1.0 and 1.8 solar masses) and for each equation of state, the star's angular momentum is varied from J=0 to the Keplerian limit J=J_{max}. For each neutron-star configuration we compute Q, the quadrupole moment of the mass distribution. We show that for given values of M and J, |Q| increases with the stiffness of the equation of state. For fixed mass and equation of state, the dependence on J is well reproduced with a simple quadratic fit, Q \simeq - aJ^2/M c^2, where c is the speed of light, and a is a parameter of order unity depending on the mass and the equation of state.