Neutron Star Equation Of State Research Articles

Impact of the PSR J0740+6620 radius constraint on the properties of high-density matter

X-ray pulse profile modeling of PSR J0740+6620, the most massive known pulsar, with data from the NICER and XMM-Newton observatories recently led to a measurement of its radius. We investigate this measurement's implications for the neutron star equation of state (EoS), employing a nonparametric EoS model based on Gaussian processes and combining information from other x-ray, radio and gravitational-wave observations of neutron stars. Our analysis mildly disfavors EoSs that support a disconnected hybrid star branch in the mass-radius relation, a proxy for strong phase transitions, with a Bayes factor of $6.9$. For such EoSs, the transition mass from the hadronic to the hybrid branch is constrained to lie outside ($1,2$) $M_{\odot}$. We also find that the conformal sound-speed bound is violated inside neutron star cores, which implies that the core matter is strongly interacting. The squared sound speed reaches a maximum of $0.75^{+0.25}_{-0.24}\, c^2$ at $3.60^{+2.25}_{-1.89}$ times nuclear saturation density at 90% credibility. Since all but the gravitational-wave observations prefer a relatively stiff EoS, PSR J0740+6620's central density is only $3.57^{+1.3}_{-1.3}$ times nuclear saturation, limiting the density range probed by observations of cold, nonrotating neutron stars in $\beta$-equilibrium.

Bayesian Nonparametric Inference of the Neutron Star Equation of State via a Neural Network

Abstract We develop a new nonparametric method to reconstruct the equation of state (EoS) of a neutron star with multimessenger data. As a universal function approximator, the feed-forward neural network (FFNN) with one hidden layer and a sigmoidal activation function can approximately fit any continuous function. Thus, we are able to implement the nonparametric FFNN representation of the EoSs. This new representation is validated by its capability of fitting the theoretical EoSs and recovering the injected parameters. Then, we adopt this nonparametric method to analyze the real data, including the mass–tidal deformability measurement from the binary neutron star merger gravitational-wave event GW170817 and mass–radius measurement of PSR J0030+0451 by NICER. We take the publicly available samples to construct the likelihood and use the nested sampling to obtain the posteriors of the parameters of the FFNN according to the Bayesian theorem, which in turn can be translated to the posteriors of the EoS parameters. Combining all of these data for a canonical 1.4 M ⊙ neutron star, we get a radius R 1.4 = 11.83 − 1.08 + 1.25 km and tidal deformability Λ 1.4 = 323 − 165 + 334 (90% confidence interval). Furthermore, we find that in the high-density region (≥3ρ sat), the 90% lower limits of c s 2 / c 2 (where c s is the sound speed and c is the velocity of light in vacuum) are above 1/3, which means that the so-called conformal limit (i.e., c s 2 / c 2 < 1 / 3 ) is not always valid in the neutron stars.

Holographic bound on area of a compact binary merger remnant

Using concomitantly the generalized second law of black hole thermodynamics and the holographic Bekenstein entropy bound embellished by loop quantum gravity corrections to quantum black hole entropy, we show that the boundary cross sectional area of the postmerger remnant formed from the compact binary merger in gravitational wave detection experiments like GW150914 et seq., by the LIGO-VIRGO collaboration, is bounded from below. This lower bound is more general than the bound obtained from application of Hawking's classical area theorem for black holes, since it does not depend on whether the inspiralling compact binary pair or the postmerger remnant consists of black holes or other exotic compact objects. The derivation of the bound entails an estimate of the entropy of the gravitational waves emitted during the binary merger, which adapts to gravitational waves an extant formalism proposed originally for particle ensembles. The results for the minimal cross sectional area of the merger remnant due to binary compact mergers observed recently by the LIGO-VIRGO collaboration are discussed. While accurate measurement of the mass of the remnant for the binary neutron star merger GW170817 remains a challenge, we provide a proof of principle that for binary neutron star mergers our lower bound on the cross sectional area of the remnant provides an alternative approach to probe the validity of neutron star equations of state, independent of the tidal deformabilities of the components.

Tight multimessenger constraints on the neutron star equation of state from GW170817 and a forward model for kilonova light-curve synthesis

ABSTRACT We present a rapid analytic framework for predicting kilonova light curves following neutron star (NS) mergers, where the main input parameters are binary-based properties measurable by gravitational wave detectors (chirp mass and mass ratio, orbital inclination) and properties dependent on the nuclear equation of state (tidal deformability, maximum NS mass). This enables synthesis of a kilonova sample for any NS source population, or determination of the observing depth needed to detect a live kilonova given gravitational wave source parameters in low latency. We validate this code, implemented in the public mosfit package, by fitting it to GW170817. A Bayes factor analysis overwhelmingly (B > 1010) favours the inclusion of an additional luminosity source in addition to lanthanide-poor dynamical ejecta during the first day. This is well fit by a shock-heated cocoon model, though differences in the ejecta structure, opacity or nuclear heating rate cannot be ruled out as alternatives. The emission thereafter is dominated by a lanthanide-rich viscous wind. We find the mass ratio of the binary is q = 0.92 ± 0.07 (90 per cent credible interval). We place tight constraints on the maximum stable NS mass, MTOV $=2.17^{+0.08}_{-0.11}$ M⊙. For a uniform prior in tidal deformability, the radius of a 1.4-M⊙ NS is R1.4 ∼ 10.7 km. Re-weighting with a prior based on equations of state that support our credible range in MTOV, we derive a final measurement R1.4 $=11.06^{+1.01}_{-0.98}$ km. Applying our code to the second gravitationally detected NS merger, GW190425, we estimate that an associated kilonova would have been fainter (by ∼0.7 mag at 1 d post-merger) and declined faster than GW170817, underlining the importance of tuning follow-up strategies individually for each GW-detected NS merger.

Towards mitigation of apparent tension between nuclear physics and astrophysical observations by improved modeling of neutron star matter

Observations of neutron stars (NSs) by the LIGO-Virgo and NICER collaborations have provided reasonably precise measurements of their various macroscopic properties. In this paper, we employ a Bayesian framework to combine them and place improved joint constraints on the properties of the NS equation of state (EOS). We use a hybrid EOS formulation that employs a parabolic expansion-based nuclear empirical parametrization around the nuclear saturation density augmented by a generic 3-segment piecewise polytrope model at higher densities. Within the 90% credible level this parametrization predicts ${R}_{1.4}={12.57}_{\ensuremath{-}0.92}^{+0.73}\text{ }\text{ }\mathrm{km}$ and ${\mathrm{\ensuremath{\Lambda}}}_{1.4}={550}_{\ensuremath{-}225}^{+223}$ for the radius and dimensionless tidal deformability, respectively, of a $1.4\text{ }\text{ }{M}_{\ensuremath{\bigodot}}$ NS. Finally, we show how the construction of the full NS EOS based solely on the nuclear empirical parameters at saturation density leads to certain tension with the astrophysical data, and how the hybrid approach provides a resolution to it.

GW190814: on the properties of the secondary component of the binary

ABSTRACT We show that the odds of the mass-gap (secondary) object in GW190814 being a neutron star (NS) improve if one allows for a stiff high-density equation of state (EoS) or a large spin. Since its mass is $\in (2.50,2.67) \, \mathrm{M}_{\odot }$, establishing its true nature will make it either the heaviest NS or the lightest black hole (BH), and can have far-reaching implications on NS EoS and compact object formation channels. When limiting oneself to the NS hypothesis, we deduce the secondary’s properties by using a Bayesian framework with a hybrid EoS formulation that employs a parabolic expansion-based nuclear empirical parametrization around the nuclear saturation density augmented by a generic 3-segment piecewise polytrope (PP) model at higher densities and combining a variety of astrophysical observations. For the slow-rotation scenario, GW190814 implies a very stiff EoS and a stringent constraint on the EoS specially in the high-density region. On the other hand assuming the secondary object is a rapidly rotating NS, we constrain its rotational frequency to be $f=1170^{+389}_{-495}$ Hz, within a 90 per cent confidence interval (CI). In this scenario, the secondary object in GW190814 would qualify as the fastest rotating NS ever observed. However, for this scenario to be viable, rotational instabilities would have to be suppressed both during formation and the subsequent evolution until merger, otherwise the secondary of GW190814 is more likely to be a BH.

AT2017gfo: Bayesian inference and model selection of multicomponent kilonovae and constraints on the neutron star equation of state

ABSTRACTThe joint detection of the gravitational wave GW170817, of the short γ-ray burst GRB170817A and of the kilonova AT2017gfo, generated by the the binary neutron star (NS) merger observed on 2017 August 17, is a milestone in multimessenger astronomy and provides new constraints on the NS equation of state. We perform Bayesian inference and model selection on AT2017gfo using semi-analytical, multicomponents models that also account for non-spherical ejecta. Observational data favour anisotropic geometries to spherically symmetric profiles, with a log-Bayes’ factor of ∼104, and favour multicomponent models against single-component ones. The best-fitting model is an anisotropic three-component composed of dynamical ejecta plus neutrino and viscous winds. Using the dynamical ejecta parameters inferred from the best-fitting model and numerical–relativity relations connecting the ejecta properties to the binary properties, we constrain the binary mass ratio to q < 1.54 and the reduced tidal parameter to $120\lt \tilde{\Lambda }\lt 1110$. Finally, we combine the predictions from AT2017gfo with those from GW170817, constraining the radius of a NS of 1.4 M⊙ to 12.2 ± 0.5 km (1σ level). This prediction could be further strengthened by improving kilonova models with numerical-relativity information.

The Role of Core-collapse Physics in the Observability of Black Hole Neutron Star Mergers as Multimessenger Sources

Abstract Recent 1D core-collapse simulations indicate a nonmonotonicity of the explodability of massive stars with respect to their precollapse core masses, which is in contrast to commonly used prescriptions. In this work, we explore the implications of these results on the formation of coalescing black hole (BH)–neutron star (NS) binaries. Furthermore, we investigate the effects of natal kicks and the NS’s radius on the synthesis of such systems and potential electromagnetic counterparts (EMCs) linked to them. Models based on 1D core-collapse simulations result in a BH–NS merger detection rate ( ∼ 2.3 yr−1), 5–10 times larger than the predictions of “standard” prescriptions. This is primarily due to the formation of low-mass BHs via direct collapse, and hence no natal kicks, favored by the 1D simulations. The fraction of observed systems that will produce an EMC, with the supernova engine from 1D simulations, ranges from 2% to 25%, depending on the NS equation of state. Notably, in most merging systems with EMCs, the NS is the first-born compact object, as long as the NS’s radius is ≲ 12 km. Furthermore, models with negligible kicks for low-mass BHs increase the detection rate of GW190426_152155-like events to ∼ 0.6 yr−1, with an associated probability of EMC ≤10% for all supernova engines. Finally, models based on 1D core-collapse simulations predict a ratio of BH–NSs to binary BHs’ merger rate density that is at least twice as high as other prescriptions, but at the same time overpredicting the measured local merger density rate of binary black holes.

Investigation of $$\varXi ^- nn$$ ($$S=-2$$) hypernucleus in low-energy pionless halo effective theory

We study the \(\varXi ^- nn\) (\(S=-2,\,I=3/2,\,J^P={1/2}^+\)) three-body system using low-energy effective field theory (EFT). Due to the acute inadequacy of empirical information in this sector, there exists substantial degree of ambiguity in determining various few-body observables, some of which are expected to yield vital clues to resolving longstanding contentious issues in hypernuclear physics. Moreover, in astrophysical studies, a precise determination of neutron star equation of state (EoS) of putative hyperonic cores relies on essential input from the \(S=-2\) sector. In this obscure current scenario, a pionless EFT analysis provides a systematic model-independent framework for assessing the feasibility of light three-particle-stable bound states, utilizing low-energy universality. Here we take recourse to a simplistic speculation of the three-body system by eliminating the repulsive spin-singlet \(\varXi ^- n\) sub-system, while retaining the predominantly attractive (possibly bound) spin-triplet \(\varXi ^- n\) and the virtual bound spin-singlet nn sub-systems. In particular, a qualitative leading order EFT investigation by introducing a sharp momentum ultraviolet cut-off parameter \(\varLambda _{\mathrm{reg}}\) into the coupled integral equations indicates a discrete scaling behavior akin to a renormalization group limit cycle, thereby suggesting the formal existence of Efimov states in the unitary limit, as \(\varLambda _{\mathrm{reg}}\rightarrow \infty \). Our subsequent non-asymptotic analysis indicates that the three-body binding energy \(B_3\) is sensitively dependent on the cut-off without the inclusion of three-body contact interactions. Furthermore, our analysis reproduces several values of the binding energy \(B_3\sim 3{-}4\) MeV, predicted in context of existing potential models, with the regulator \(\varLambda _{\mathrm{reg}}\) in the range \(\sim 350{-}460\) MeV. Finally, based on these model inputs for \(B_3\), a ballpark estimate of the three-body scattering length in the range \(2.6{-}4.9\) fm, is naively constrained by our EFT analysis. Despite approximations, the resulting Phillips line is expected to yield a robust feature of the halo-bound \(\varXi ^- nn\) system. For pedagogical reasons, using a simple toy model interacting three-bosons system, we highlight in the appendices the typical universal features leading to emergence of RG limit cycle and Efimov states which are amenable to a low-energy EFT formalism.

Constraining neutron star radii in black hole–neutron star mergers from their electromagnetic counterparts

ABSTRACT Mergers of black hole (BH) and neutron star (NS) binaries are of interest since the emission of gravitational waves (GWs) can be followed by an electromagnetic (EM) counterpart, which could power short gamma-ray bursts. Until now, LIGO/Virgo has only observed a candidate BH–NS event, GW190426_152155, which was not followed by any EM counterpart. We discuss how the presence (absence) of a remnant disc, which powers the EM counterpart, can be used along with spin measurements by LIGO/Virgo to derive a lower (upper) limit on the radius of the NS. For the case of GW190426_152155, large measurement errors on the spin and mass ratio prevent from placing an upper limit on the NS radius. Our proposed method works best when the aligned component of the BH spin (with respect to the orbital angular momentum) is the largest, and can be used to complement the information that can be extracted from the GW signal to derive valuable information on the NS equation of state.

Extensive studies of the neutron star equation of state from the deep learning inference with the observational data augmentation

We discuss deep learning inference for the neutron star equation of state (EoS) using the real observational data of the mass and the radius. We make a quantitative comparison between the conventional polynomial regression and the neural network approach for the EoS parametrization. For our deep learning method to incorporate uncertainties in observation, we augment the training data with noise fluctuations corresponding to observational uncertainties. Deduced EoSs can accommodate a weak first-order phase transition, and we make a histogram for likely first-order regions. We also find that our observational data augmentation has a byproduct to tame the overfitting behavior. To check the performance improved by the data augmentation, we set up a toy model as the simplest inference problem to recover a double-peaked function and monitor the validation loss. We conclude that the data augmentation could be a useful technique to evade the overfitting without tuning the neural network architecture such as inserting the dropout.

Toward observing neutron star collapse with gravitational wave detectors

Gravitational waves from binary neutron star inspirals have been detected along with the electromagnetic transients coming from the aftermath of the merger in GW170817. However, much is still unknown about the post-merger dynamics that connects these two sets of observables. This includes if, and when, the post-merger remnant star collapses to a black hole, and what are the necessary conditions to power a short gamma-ray burst, and other observed electromagnetic counterparts. Observing the collapse of the post-merger neutron star would shed led on these questions, constraining models for the short gamma-ray burst engine and the hot neutron star equation of state. In this work, we explore the scope of using gravitational wave detectors to measure the timing of the collapse either indirectly, by establishing the shut-off of the post-merger gravitational emission, or---more challengingly---directly, by detecting the collapse signal. For the indirect approach, we consider a kilohertz high-frequency detector design that utilises a previously studied coupled arm cavity and signal recycling cavity resonance. This design would give a signal-to-noise ratio of 0.5\,-\,8.6 (depending on the variation of waveform parameters) for a collapse gravitational wave signal occurring at 10\,ms post-merger of a binary at 50\,Mpc and with total mass $2.7 M_\odot$. For the direct approach, we propose a narrow-band detector design, utilising the sensitivity around the frequency of the arm cavity free spectral range. The proposed detector achieves a signal-to-noise ratio of 0.3\,-\,1.9, independent of the collapse time. This detector is limited by both the fundamental classical and quantum noise with the arm cavity power chosen as 10\,MW.

Properties of the remnant disk and the dynamical ejecta produced in low-mass black hole-neutron star mergers

We systematically perform numerical-relativity simulations for low-mass black hole-neutron star mergers for the models with seven mass ratios $Q=M_{\rm BH}/M_{\rm NS}$ ranging from 1.5 to 4.4, and three neutron-star equations of state, focusing on properties of matter remaining outside the black hole and ejected dynamically during the merger. We pay particular attention to the dependence on the mass ratio of the binaries. It is found that the rest mass remaining outside the apparent horizon after the merger depends only weakly on the mass ratio for the models with low mass ratios. It is also clarified that the rest mass of the ejecta has a peak at $Q \sim 3$, and decreases steeply as the mass ratio decreases for the low mass-ratio case. We present a novel analysis method for the behavior of matter during the merger, focusing on the matter distribution in the phase space of specific energy and specific angular momentum. Then we model the matter distribution during and after the merger. Using the result of the analysis, we discuss the properties of the ejecta.

Optimized Statistical Approach for Comparing Multi-messenger Neutron Star Data

Abstract The neutron star equation of state is now being constrained from a diverse set of multi-messenger data, including gravitational waves from binary neutron star mergers, X-ray observations of the neutron star radius, and many types of laboratory nuclear experiments. These measurements are often mapped to a common domain for comparison with one another or are used to constrain the predictions of theoretical equations of state. We explore here the statistical biases that can arise when such multi-messenger data are compared or combined across different domains. We find that placing Bayesian priors individually in each domain of measurement can lead to biased constraints. We present a new prescription for defining Bayesian priors consistently across different experiments, which will allow for robust cross-domain comparisons. Using the first two binary neutron star mergers as an example, we show that a uniform prior in the tidal deformability can produce inflated evidence for large radii, while a uniform prior in the radius points toward smaller radii. Finally, using this new prescription, we provide a status update on multi-messenger constraints on the neutron star radius.

Probing the Nuclear Equation of State from the Existence of a ∼2.6 M⊙ Neutron Star: The GW190814 Puzzle

On 14 August 2019, the LIGO/Virgo collaboration observed a compact object with mass ∼2.59−0.09+0.08M⊙, as a component of a system where the main companion was a black hole with mass ∼23M⊙. A scientific debate initiated concerning the identification of the low mass component, as it falls into the neutron star–black hole mass gap. The understanding of the nature of GW190814 event will offer rich information concerning open issues, the speed of sound and the possible phase transition into other degrees of freedom. In the present work, we made an effort to probe the nuclear equation of state along with the GW190814 event. Firstly, we examine possible constraints on the nuclear equation of state inferred from the consideration that the low mass companion is a slow or rapidly rotating neutron star. In this case, the role of the upper bounds on the speed of sound is revealed, in connection with the dense nuclear matter properties. Secondly, we systematically study the tidal deformability of a possible high mass candidate existing as an individual star or as a component one in a binary neutron star system. As the tidal deformability and radius are quantities very sensitive on the neutron star equation of state, they are excellent counters on dense matter properties. We conjecture that similar isolated neutron stars or systems may exist in the universe and their possible future observation will shed light on the maximum neutron star mass problem.

Neutron Star Equation of State in Light of GW190814.

The observation of gravitational waves from an asymmetric binary opens the possibility for heavy neutron stars, but these pose challenges to models of the neutron star equation of state. We construct heavy neutron stars by introducing nontrivial structure in the speed of sound sourced by deconfined QCD matter, which cannot be well recovered by spectral representations. Their moment of inertia, Love number, and quadrupole moment are very small, so a tenfold increase in sensitivity may be needed to test this possibility with gravitational waves, which is feasible with third generation detectors.

Inference of the Neutron Star Equation of State from Cosmological Distances.

Finite-size effects on the gravitational wave signal from a neutron star merger typically manifest at high frequencies where detector sensitivity decreases. Proposed sensitivity improvements can give us access both to stronger signals and to a myriad of weak signals from cosmological distances. The latter will outnumber the former and the relevant part of the signal will be redshifted towards the detector's most sensitive band. We study the redshift dependence of information about neutron star matter and find that single-scale properties, such as the star radius or the postmerger frequency, are better measured from the distant weak sources from z∼1.

Weighing massive neutron star with screening gravity: a look on PSR J0740\xa0+\xa06620 and GW190814 secondary component

Neutron stars (NSs) are excellent natural laboratories to constrain gravity on strong field regime and nuclear matter in extreme conditions. Motivated by the recent discovery of a compact object with 2.59^{+0.08}_{-0.09} M_odot in the binary merger GW190814, if this object was a NS, it serves as a strong constraint on the NS equation of state (EoS), ruling out several soft EoSs favored by GW170817 event. In this work, we revisit the question of the maximum mass of NSs considering a chameleon screening (thin-shell effect) on the NS mass-radius relation, where the microscopic physics inside the NS is given by realistic soft EoSs. We find that from appropriate and reasonable combination of modified gravity, rotation effects and realistic soft EoSs, that it is possible to achieve high masses and explain GW190814 secondary component, and in return also NSs like PSR J0740 + 6620 (the most NS massive confirmed to date). It is shown that gravity can play an important role in estimating maximum mass of NSs, and even with soft EoSs, it is possible to generate very high masses. Therefore, in this competition of hydrostatic equilibrium between gravity and pressure (from EoS choice), some soft EoSs, in principle, cannot be completely ruled out without first taking into account gravitational effects.

A new consistent neutron star equation of state from a generalized Skyrme model

We propose a new equation of state for nuclear matter based on a generalized Skyrme model which is consistent with all current constraints on the observed properties of neutron stars. This generalized model depends only on two free parameters related to the ranges of pressure values at which different submodels are dominant, and which can be adjusted so that mass-radius and deformability constraints from astrophysical and gravitational wave measurements can be met. Our results support the Skyrme model and its generalizations as good candidates for a low energy effective field-theoretic description of nuclear matter even at extreme conditions such as those inside neutron stars.

Speed of sound constraints from tidal deformability of neutron stars

The upper bound of the speed of sound in dense nuclear matter is one of the most interesting but still unsolved problems in Nuclear Physics. Theoretical studies in connection with recent observational data of isolated neutron stars as well as binary neutron stars systems offer an excellent opportunity to shed light on this problem. In the present work, we suggest a method to directly relate the measured tidal deformability (polarizability) of binary neutron stars system (before merger) to the maximum neutron star mass scenario and possible upper bound on the speed of sound. This method is based on the simple but efficient idea that while the upper limit of the effective tidal deformability favors soft equations of state, the recent high measured values of neutron star mass favor stiff ones. In the present work, firstly, using a simple well established model we parametrize the stiffness of the equation of state with the help of the speed of sound. Secondly, in comparison with the recent observations by LIGO/VIRGO collaboration of two events, GW170817 and GW190425, we suggest possible robust constraints. Moreover, we evaluate and postulate, in the framework of the present method, what kind of future measurements could help us to improve the stringent of the constraints on the neutron star equation of state.