Neutron Stars In General Relativity Research Articles

Ridges in rotating neutron-star properties due to first order phase transitions

We identify combinations of observables for rotating neutron stars that can one day bear on the question of whether there can be first order phase transitions in the neutron matter therein. We employ the Hartle–Thorne theory for stationary, rotating neutron stars at conventional angular velocities (in the conventional pulsar and millisecond pulsar ranges) and extract three-dimensional sections of the ellipticity or the dynamical angular momentum as function of the star’s mass and angular velocity. An eventual first order phase transition in the equation of state (EoS) leaves a clear ridge (nonanalyticity) in these observables, akin to the sudden kink in popular mass–radius diagrams for static stars. Finally, we observe that static neutron stars in General Relativity (GR) will fail to be compact enough for the light ring’s position at r=3M to be outside the star, except for the most extreme equations of state. The outer light ring of a rotating star might however be formed unless the EoS softens too much, and its eventual detection can then be used to constrain the EoS (or the gravity theory).

Magnetically confined mountains on accreting neutron stars in general relativity

ABSTRACT The general relativistic formulation of the problem of magnetically confined mountains on neutron stars is presented, and the resulting equations are solved numerically, generalizing previous Newtonian calculations. The hydromagnetic structure of the accreted matter and the subsequent magnetic burial of the star’s magnetic dipole moment are computed. Overall, it is observed that relativistic corrections reduce the hydromagnetic deformation associated with the mountain. The magnetic field lines are curved more gently than in previous calculations, and the screening of the dipole moment is reduced. Quantitatively, it is found that the dimensionless dipole moment (md) depends on the accreted mass (Ma) as md = −3.2 × 103Ma/M⊙ + 1.0, implying approximately three times less screening compared to the Newtonian theory. Additionally, the characteristic scale height of the mountain, governing the gradients of quantities like pressure, density, and magnetic field strength, reduces by approximately 40 per cent for an isothermal equation of state.

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.

Quasi-periodic oscillations in precursor flares via seismic aftershocks from resonant shattering

Context. GRB 211211A was accompanied by a kilonova, implying a merger origin for the event. A precursor flare, modulated by quasi-periodic oscillations at ≈22 Hz, was also seen ≈1 s prior to the main emission. Aims. It is suggested here that the precursor resulted from the resonant shattering of one star’s crust prior to coalescence. Seismic aftershocks and low-frequency torsional modes may emanate from the rupture site in this case, explaining the underlying oscillations. This interpretation is directly analogous to proposals for giant flares from soft gamma repeaters, where oscillations at similar frequencies have been observed, involving starquakes followed by crustal vibrations. Methods. Torsional mode properties were computed for sequences of slowly rotating, magnetised neutron stars in general relativity. Results. The ≈22 Hz modulations in the precursor matches that of nodeless, ℓ = 2 torsional modes for a wide variety of stellar parameters. Global or discrete Alfvén modes are also viable candidates.

Neutron stars in massive scalar-Gauss-Bonnet gravity: Spherical structure and time-independent perturbations

The class of scalar-tensor theories with the scalar field coupling to the Gauss-Bonnet invariant has drawn great interest since solutions of spontaneous scalarization were found for black holes in these theories. We contribute to the existing literature a detailed study of the spontaneously scalarized neutron stars (NSs) in a typical theory where the coupling function of the scalar field takes the quadratic form and the scalar field is massive. The investigation here includes the spherical solutions of the NSs as well as their perturbative properties, namely the tidal deformability and the moment of inertia, treated in a unified and extendable way under the framework of spherical decomposition. We find that while the mass, the radius, and the moment of inertia of the spontaneously scalarized NSs show very moderate deviations from those of the NSs in general relativity (GR), the tidal deformability exhibits significant differences between the solutions in GR and the solutions of spontaneous scalarization for certain values of the parameters in the scalar-Gauss-Bonnet theory. As a result, the celebrated universal relation between the moment of inertia and the tidal deformability of neutron stars breaks down. With the mass and the tidal deformability of NSs attainable in the gravitational waves from binary NS mergers, the radius measurable using the X-ray satellites, and the moment of inertia accessible via the high-precision pulsar timing techniques, future multi-messenger observations can be contrasted with the theoretical results and provide us necessary information for building up theories beyond GR.

Quasi-universality of the magnetic deformation of neutron stars in general relativity and beyond

Neutron stars are known to host extremely powerful magnetic fields. Among its effects, one of the consequences of harbouring such fields is the deformation of the neutron star structure, leading, together with rotation, to the emission of continuous gravitational waves. On the one hand, the details of their internal magnetic fields are mostly unknown. Likewise, their internal structure, encoded by the equation of state, is highly uncertain. Here, we present a study of axisymmetric models of isolated magnetised neutron stars for various realistic equations of state considered viable by observations and nuclear physics constraints. We show that it is possible to find simple relations between the magnetic deformation of a neutron star, its Komar mass, and its circumferential radius in the case of purely poloidal and purely toroidal magnetic configurations that satisfy the criterion for equilibrium in the Bernoulli formalism. Such relations are quasi-universal, meaning that they are mostly independent from the equation of state of the neutron star. Thanks to their formulation in terms of potentially observable quantities, as we discuss, our results could help to constrain the magnetic properties of the neutron star interior and to better assess the detectability of continuous gravitational waves by isolated neutron stars, without knowing their equation of state. Our results are derived both in general relativity and in scalar-tensor theories (one of the most promising extensions of general relativity), in this case by also considering the scalar charge. We show that even in this case, general relations that account for deviations from general relativity still hold, which could potentially be used to set constraints on the gravitational theory.

A new class of analytical solutions describing anisotropic neutron stars in general relativity

A new class of solutions describing analytical solutions for compact stellar structures has been developed within the tenets of General Relativity. Considering the inherent anisotropy in compact stars, a stable and causal model for realistic anisotropic neutron stars was obtained using the general theory of relativity. Assuming a physically acceptable nonsingular form of one metric potential and radial pressure containing the curvature parameter [Formula: see text], the constant [Formula: see text] and the radius [Formula: see text], analytical solutions to Einstein’s field equations for anisotropic matter distribution were obtained. Taking the value of [Formula: see text] as −0.44, it was found that the proposed model obeys all necessary physical conditions, and it is potentially stable and realistic. The model also exhibits a linear equation of state, which can be applied to describe compact stars.

Neutron star structure in Hořava-Lifshitz gravity

We present interesting aspects of neutron stars (NSs) from the standpoint of a modified theory of gravity called Ho\v{r}ava-Lifshitz (HL) gravity. A deviation from general relativity (GR) in HL gravity can change typical features of the NS structure. In this study, we investigate the NS structure by deriving the Tolman-Oppenheimer-Volkoff equation in HL gravity. We find that a NS in HL gravity with a larger radius and heavier mass than a NS in GR remains stable without collapsing into a black hole.

Axisymmetric equilibrium models for magnetised neutron stars in scalar-tensor theories

Among the possible extensions of general relativity that have been put forward to address some long-standing issues in our understanding of the Universe, scalar-tensor theories have received a lot of attention for their simplicity. Interestingly, some of these predict a potentially observable non-linear phenomenon, known as spontaneous scalarisation, in the presence of highly compact matter distributions, as in the case of neutron stars. Neutron stars are ideal laboratories for investigating the properties of matter under extreme conditions and, in particular, they are known to harbour the strongest magnetic fields in the Universe. Here, for the first time, we present a detailed study of magnetised neutron stars in scalar-tensor theories. First, we showed that the formalism developed for the study of magnetised neutron stars in general relativity, based on the “extended conformally flat condition”, can easily be extended in the presence of a non-minimally coupled scalar field, retaining many of its numerical advantages. We then carried out a study of the parameter space considering the two extreme geometries of purely toroidal and purely poloidal magnetic fields, varying both the strength of the magnetic field and the intensity of scalarisation. We compared our results with magnetised general-relativistic solutions and un-magnetised scalarised solutions, showing how the mutual interplay between magnetic and scalar fields affect the magnetic and the scalarisation properties of neutron stars. In particular, we focus our discussion on magnetic deformability, maximum mass, and range of scalarisation.

Neutron stars in general relativity and scalar-tensor theory of gravity

The masses and radii of neutron stars are discussed in general relativity and scalar-tensor theory of gravity and the differences are compared with the current uncertainties stemming from the nuclear equation of state in the relativistic mean-field framework. It is shown that astrophysical and gravitational wave observations of radii of neutron stars with masses M lesssim 1.4 M_{odot } constrain only the nuclear equation of state, and in particular the density dependence of the nuclear symmetry energy. Future observations of massive neutron stars may constrain the coupling parameters of the scalar-tensor theory provided that a general consensus on the dense nuclear matter equation of state is reached.

Revisiting the maximum mass of differentially rotating neutron stars in general relativity with realistic equations of state

We study the solution space of general relativistic, axisymmetric, equilibria of differentially rotating neutron stars with realistic, nuclear equations of state. We find that different types of stars, which were identified by earlier works for polytropic equations of state, arise for realistic equations of state, too. Scanning the solution space for the sample of realistic equations of state we treat, we find lower limits on the maximum rest masses supported by cold, differentially rotating stars for each type of stars. We often discover equilibrium configurations that can support more than 2 times the mass of a static star. We call these equilibria "\"ubermassive", and in our survey we find \"ubermassive stars that can support up to 2.5 times the maximum rest mass that can be supported by a cold, non-rotating star with the same equation of state. This is nearly two times larger than what previous studies employing realistic equations of state had found, and which did not uncover \"ubermassive neutron stars. Moreover, we find that the increase in the maximum rest mass with respect to the non-spinning stellar counterpart is larger for moderately stiff equations of state. These results may have implications for the lifetime and the gravitational wave and electromagnetic counterparts of hypermassive neutron stars formed following binary neutron star mergers.

Compactness of neutron stars and Tolman VII solutions in scalar-tensor gravity

We systematically examine the compactness of neutron stars as Tolman VII solutions in scalar-tensor theory of gravity. As a result, when the coupling constant is confined to values provided by astronomical observations we show that the maximum compactness of neutron stars in general relativity is higher than that in scalar-tensor gravity. In addition, we show that although ultra-compact stars, with radius smaller than the Regge-Wheeler potential peak, can exist in general relativity (e.g., Tolman VII solution), their scalarized counterparts cannot {be constructed} even in the limiting case of uniform density stars.

Relativistic MHD modeling of magnetized neutron stars, pulsar winds, and their nebulae

Neutron stars are among the most fascinating astrophysical sources, being characterized by strong gravity, densities about the nuclear one or even above, and huge magnetic fields. Their observational signatures can be extremely diverse across the electromagnetic spectrum, ranging from the periodic and low-frequency signals of radio pulsars, up to the abrupt high-energy gamma-ray flares of magnetars, where energies of are released in a few seconds. Fast-rotating and highly magnetized neutron stars are expected to launch powerful relativistic winds, whose interaction with the supernova remnants gives rise to the non-thermal emission of pulsar wind nebulae, which are known cosmic accelerators of electrons and positrons up to PeV energies. In the extreme cases of proto-magnetars (magnetic fields of G and millisecond periods), a similar mechanism is likely to provide a viable engine for the still mysterious gamma-ray bursts. The key ingredient in all these spectacular manifestations of neutron stars is the presence of strong magnetic fields in their constituent plasma. Here we will present recent updates of a couple of state-of-the-art numerical investigations by the high-energy astrophysics group in Arcetri: a comprehensive modeling of the steady-state axisymmetric structure of rotating magnetized neutron stars in general relativity, and dynamical 3D MHD simulations of relativistic pulsar winds and their associated nebulae.

Gravitational wave content and stability of uniformly, rotating, triaxial neutron stars in general relativity.

Targets for ground-based gravitational wave interferometers include continuous, quasiperiodic sources of gravitational radiation, such as isolated, spinning neutron stars. In this work, we perform evolution simulations of uniformly rotating, triaxially deformed stars, the compressible analogs in general relativity of incompressible, Newtonian Jacobi ellipsoids. We investigate their stability and gravitational wave emission. We employ five models, both normal and supramassive, and track their evolution with different grid setups and resolutions, as well as with two different evolution codes. We find that all models are dynamically stable and produce a strain that is approximately one-tenth the average value of a merging binary system. We track their secular evolution and find that all our stars evolve toward axisymmetry, maintaining their uniform rotation, rotational kinetic energy, and angular momentum profiles while losing their triaxiality.

Numerical models for stationary superfluid neutron stars in general relativity with realistic equations of state

We present a numerical model for uniformly rotating superfluid neutron stars, for the first time with realistic microphysics including entrainment, in a fully general relativistic framework. We compute stationary and axisymmetric configurations of neutron stars composed of two fluids, namely superfluid neutrons and charged particles (protons and electrons), rotating with different rates around a common axis. Both fluids are coupled by entrainment, a non-dissipative interaction which in case of a non-vanishing relative velocity between the fluids, causes the fluid momenta being not aligned with the respective fluid velocities. We extend the formalism by Comer and Joynt (2003) in order to calculate the equation of state (EoS) and entrainment parameters for an arbitrary relative velocity. The resulting entrainment matrix fulfills all necessary sum rules and in the limit of small relative velocity our results agree with Fermi liquid theory ones, derived to lowest order in the velocity. This formalism is applied to two new nuclear equations of state which are implemented in the numerical model. We are able to obtain precise equilibrium configurations. Resulting density profiles and moments of inertia are discussed employing both EoSs, showing the impact of entrainment and the dependence on the EoS.

Slowly rotating anisotropic neutron stars in general relativity and scalar–tensor theory

Some models (such as the Skyrme model, a low-energy effective field theory for quantum chromodynamics) suggest that the high-density matter prevailing in neutron star (NS) interiors may be significantly anisotropic. Anisotropy is known to affect the bulk properties of nonrotating NSs in general relativity (GR). In this paper we study the effects of anisotropy on slowly rotating stars in GR. We also consider one of the most popular extensions of Einstein’s theory, namely scalar–tensor theories allowing for spontaneous scalarization (a phase transition similar to spontaneous magnetization in ferromagnetic materials). Anisotropy affects the moment of inertia of NSs (a quantity that could potentially be measured in binary pulsar systems) in both theories. We find that the effects of scalarization increase (decrease) when the tangential pressure is bigger (smaller) than the radial pressure, and we present a simple criterion to determine the onset of scalarization by linearizing the scalar-field equation. Our calculations suggest that binary pulsar observations may constrain the degree of anisotropy or even, more optimistically, provide evidence for anisotropy in NS cores.

Beyond 2nd order in the simulations of binary neutron stars in general relativity

Beyond 2nd order in the simulations of binary neutron stars in general relativity

Effectively universal behavior of rotating neutron stars in general relativity makes them even simpler than their Newtonian counterparts.

Recently, it was shown that slowly rotating neutron stars exhibit an interesting correlation between their moment of inertia I, their quadrupole moment Q, and their tidal deformation Love number λ (the I-Love-Q relations), independently of the equation of state of the compact object. In the present Letter a similar, more general, universality is shown to hold true for all rotating neutron stars within general relativity; the first four multipole moments of the neutron star are related in a way independent of the nuclear matter equation of state we assume. By exploiting this relation, we can describe quite accurately the geometry around a neutron star with fewer parameters, even if we don't know precisely the equation of state. Furthermore, this universal behavior displayed by neutron stars could promote them to a more promising class of candidates (next to black holes) for testing theories of gravity.

Axisymmetric equilibrium models for magnetized neutron stars in General Relativity under the Conformally Flat Condition

Extremely magnetized neutron stars with magnetic fields as strong as $\sim 10^{15-16}$ G, or magnetars, have received considerable attention in the last decade due to their identification as a plausible source for Soft Gamma Repeaters and Anomalous X-ray Pulsars. Moreover, this class of compact objects has been proposed as a possible engine capable of powering both Long and Short Gamma-Ray Bursts, if the rotation period in their formation stage is short enough (~1 ms). Such strong fields are expected to induce substantial deformations of the star and thus to produce the emission of gravitational waves. Here we investigate, by means of numerical modeling, axisymmetric static equilibria of polytropic and strongly magnetized stars in full general relativity, within the ideal magneto-hydrodynamic regime. The eXtended Conformally Flat Condition (XCFC) for the metric is assumed, allowing us to employ the techniques introduced for the X-ECHO code [Bucciantini & Del Zanna, 2011, Astron. Astrophys. 528, A101], proven to be accurate, efficient, and stable. The updated XNS code for magnetized neutron star equilibria is made publicly available for the community (see www.arcetri.astro.it/science/ahead/XNS). Several sequences of models are here retrieved, from the purely toroidal (resolving a controversy in the literature) or poloidal cases, to the so-called twisted torus mixed configurations, expected to be dynamically stable, which are solved for the first time in the non-perturbative regime.

PULSE PROFILES FROM THERMALLY EMITTING NEUTRON STARS

The problem of computing the pulse profiles from thermally emitting spots on the surface of a neutron star in general relativity is reconsidered. We show that it is possible to extend Beloborodov (2002) approach to include (multiple) spots of finite size in different positions on the star surface. Results for the pulse profiles are expressed by comparatively simple analytical formulas which involve only elementary functions.