Multimessenger constraints on the neutron-star equation of state and the Hubble constant

Time-Delay Cosmography: Measuring the Hubble Constant and Other Cosmological Parameters with Strong Gravitational Lensing

Multiply lensed images of a same source experience a relative time delay in the arrival of photons due to the path length difference and the different gravitational potentials the photons travel through. This effect can be used to measure absolute distances and the Hubble constant ( H 0 ) and is known as time-delay cosmography. The method is independent of the local distance ladder and early-universe physics and provides a precise and competitive measurement of H 0 . With upcoming observatories, time-delay cosmography can provide a 1% precision measurement of H 0 and can decisively shed light on the current reported 'Hubble tension'. This manuscript details the general methodology developed over the past decades in time-delay cosmography, discusses recent advances and results, and, foremost, provides a foundation and outlook for the next decade in providing accurate and ever more precise measurements with increased sample size and improved observational techniques.

Low-latency gravitational wave alert products and their performance at the time of the fourth LIGO-Virgo-KAGRA observing run

Multimessenger searches for binary neutron star (BNS) and neutron star-black hole (NSBH) mergers are currently one of the most exciting areas of astronomy. The search for joint electromagnetic and neutrino counterparts to gravitational wave (GW)s has resumed with ALIGO's, AdVirgo's and KAGRA's fourth observing run (O4). To support this effort, public semiautomated data products are sent in near real-time and include localization and source properties to guide complementary observations. In preparation for O4, we have conducted a study using a simulated population of compact binaries and a mock data challenge (MDC) in the form of a real-time replay to optimize and profile the software infrastructure and scientific deliverables. End-toend performance was tested, including data ingestion, running online search pipelines, performing annotations, and issuing alerts to the astrophysics community. We present an overview of the low-latency infrastructure and the performance of the data products that are now being released during O4 based on the MDC. We report the expected median latency for the preliminary alert of full bandwidth searches (29.5 s) and show consistency and accuracy of released data products using the MDC. We report the expected median latency for triggers from early warning searches (-3.1 s), which are new in O4 and target neutron star mergers during inspiral phase. This paper provides a performance overview for LIGO-Virgo-KAGRA (LVK) low-latency alert infrastructure and data products using theMDCand serves as a useful reference for the interpretation of O4 detections.

An updated nuclear-physics and multi-messenger astrophysics framework for binary neutron star mergers

The multi-messenger detection of the gravitational-wave signal GW170817, the corresponding kilonova AT2017gfo and the short gamma-ray burst GRB170817A, as well as the observed afterglow has delivered a scientific breakthrough. For an accurate interpretation of all these different messengers, one requires robust theoretical models that describe the emitted gravitational-wave, the electromagnetic emission, and dense matter reliably. In addition, one needs efficient and accurate computational tools to ensure a correct cross-correlation between the models and the observational data. For this purpose, we have developed the Nuclear-physics and Multi-Messenger Astrophysics framework NMMA. The code allows incorporation of nuclear-physics constraints at low densities as well as X-ray and radio observations of isolated neutron stars. In previous works, the NMMA code has allowed us to constrain the equation of state of supranuclear dense matter, to measure the Hubble constant, and to compare dense-matter physics probed in neutron-star mergers and in heavy-ion collisions, and to classify electromagnetic observations and perform model selection. Here, we show an extension of the NMMA code as a first attempt of analyzing the gravitational-wave signal, the kilonova, and the gamma-ray burst afterglow simultaneously. Incorporating all available information, we estimate the radius of a 1.4M⊙ neutron star to be [Formula: see text] km.

Strongly interacting matter exhibits deconfined behavior in massive neutron stars

Neutron-star cores contain matter at the highest densities in our Universe. This highly compressed matter may undergo a phase transition where nuclear matter melts into deconfined quark matter, liberating its constituent quarks and gluons. Quark matter exhibits an approximate conformal symmetry, predicting a specific form for its equation of state (EoS), but it is currently unknown whether the transition takes place inside at least some physical neutron stars. Here, we quantify this likelihood by combining information from astrophysical observations and theoretical calculations. Using Bayesian inference, we demonstrate that in the cores of maximally massive stars, the EoS is consistent with quark matter. We do this by establishing approximate conformal symmetry restoration with high credence at the highest densities probed and demonstrating that the number of active degrees of freedom is consistent with deconfined matter. The remaining likelihood is observed to correspond to EoSs exhibiting phase-transition-like behavior, treated as arbitrarily rapid crossovers in our framework.

Constraints on the Hubble constant from Supernova Refsdal's reappearance

The gravitationally lensed Supernova Refsdal appeared in multiple images, produced through gravitational lensing by a massive foreground galaxy cluster. After the supernova appeared in 2014, lens models of the galaxy cluster predicted an additional image of the supernova would appear in 2015, which was subsequently observed. We use the time delays between the images to perform a blinded measurement of the expansion rate of the Universe, quantified by the Hubble constant (H₀). Using eight cluster lens models, we infer [Formula: see text], where Mpc is the megaparsec. Using the two models most consistent with the observations, we find [Formula: see text]. The observations are best reproduced by models that assign dark-matter halos to individual galaxies and the overall cluster.

Plausible presence of new state in neutron stars with masses above 0.98MTOV

We investigate the neutron star (NS) equation of state (EOS) by incorporating multi-messenger data of GW170817, PSR J0030 + 0451, PSR J0740 + 6620, and state-of-the-art theoretical progresses, including the information from chiral effective field theory (χEFT) and perturbative quantum chromodynamics (pQCD) calculation. Taking advantage of the various structures sampling by a single-layer feed-forward neural network model embedded in the Bayesian nonparametric inference, the structure of NS matter's sound speed c_s is explored in a model-agnostic way. It is found that a peak structure is common in the c_s² posterior, locating at 2.4-4.8ρ_sat (nuclear saturation density) and c_s² exceeds c²/3 at 90% credibility. The non-monotonic behavior suggests evidence of the state deviating from the hadronic matter inside the very massive NSs. Assuming the new/exotic state is featured as it is softer than typical hadronic models or even with hyperons, we find that a sizable (⩾10^-3M_⊙) exotic core, likely made of quark matter, is plausible for the NS with a gravitational mass above about 0.98M_TOV, where M_TOV represents the maximum gravitational mass of a non-rotating cold NS. The inferred M_TOV=2.18_-0.13^+0.27M_⊙ (90% credibility) is well consistent with the value of 2.17_-0.12^+0.15M_⊙ estimated independently with GW170817/GRB 170817A/AT2017gfo assuming a temporary supramassive NS remnant formed after the merger. PSR J0740 + 6620, the most massive NS detected so far, may host a sizable exotic core with a probability of ≈0.36.

Nonideal effect of free electrons on ionization equilibrium and radiative property in dense plasmas

The thermodynamic as well as optical properties of strongly coupled plasmas depend crucially on the average degree of ionization and the ionic state composition, which, however, cannot be determined by using the normal Saha equation usually used for the ideal plasmas. Hence, an adequate treatment of the ionization balance and the charge state distribution of strongly coupled plasmas is still a challenge for theory due to the interactions between the electrons and ions and among the electrons themselves. Based on a local density temperature-dependent ion-sphere model, the Saha equation approach is extended to the regime of strongly coupled plasmas by taking into account the free-electron-ion interaction, the free-free-electron interaction, the nonuniform free-electron space distribution, and the free-electron quantum partial degeneracy. All the quantities, including the bound orbitals with ionization potential depression, free-electron distribution, and bound and free-electron partition function contributions, are calculated self-consistently in the theoretical formalism. This study shows that the ionization equilibrium is evidently modified by considering the above nonideal characteristics of the free electrons. Our theoretical formalism is validated by the explanation of a recent experimental measurement of the opacity of dense hydrocarbon.

Nuclear Equation of State for Arbitrary Proton Fraction and Temperature Based on Chiral Effective Field Theory and a Gaussian Process Emulator

We calculate the equation of state of asymmetric nuclear matter at finite temperature based on chiral effective field theory interactions to next-to-next-to-next-to-leading order. Our results assess the theoretical uncertainties from the many-body calculation and the chiral expansion. Using a Gaussian process emulator for the free energy, we derive the thermodynamic properties of matter through consistent derivatives and use the Gaussian process to access arbitrary proton fraction and temperature. This enables a first nonparametric calculation of the equation of state in beta equilibrium, and of the speed of sound and the symmetry energy at finite temperature. Moreover, our results show that the thermal part of the pressure decreases with increasing densities.

Computational challenges for multimodal astrophysics

In the coming decades, we will face major computational challenges, when the improved sensitivity of third-generation gravitational wave detectors will be such that they will be able to detect a high number (of the order of 7 × 104 per year) of multi-messenger events from binary neutron star mergers, similar to GW 170817. In this Perspective, we discuss the application of multimodal artificial intelligence techniques for multi-messenger astrophysics, fusing the information from different signal emissions.

Incorporating a Radiative Hydrodynamics Scheme in the Numerical-Relativity Code BAM

To study binary neutron star systems and to interpret observational data such as gravitational-wave and kilonova signals, one needs an accurate description of the processes that take place during the final stages of the coalescence, for example, through numerical-relativity simulations. In this work, we present an updated version of the numerical-relativity code BAM in order to incorporate nuclear-theory-based equations of state and a simple description of neutrino interactions through a neutrino leakage scheme. Different test simulations, for stars undergoing a neutrino-induced gravitational collapse and for binary neutron stars systems, validate our new implementation. For the binary neutron stars systems, we show that we can evolve stably and accurately distinct microphysical models employing the different equations of state: SFHo, DD2, and the hyperonic BHBΛϕ. Overall, our test simulations have good agreement with those reported in the literature.

Constraining neutron-star matter with microscopic and macroscopic collisions

Interpreting high-energy, astrophysical phenomena, such as supernova explosions or neutron-star collisions, requires a robust understanding of matter at supranuclear densities. However, our knowledge about dense matter explored in the cores of neutron stars remains limited. Fortunately, dense matter is not probed only in astrophysical observations, but also in terrestrial heavy-ion collision experiments. Here we use Bayesian inference to combine data from astrophysical multi-messenger observations of neutron stars1-9 and from heavy-ion collisions of gold nuclei at relativistic energies10,11 with microscopic nuclear theory calculations12-17 to improve our understanding of dense matter. We find that the inclusion of heavy-ion collision data indicates an increase in the pressure in dense matter relative to previous analyses, shifting neutron-star radii towards larger values, consistent with recent observations by the Neutron Star Interior Composition Explorer mission5-8,18. Our findings show that constraints from heavy-ion collision experiments show a remarkable consistency with multi-messenger observations and provide complementary information on nuclear matter at intermediate densities. This work combines nuclear theory, nuclear experiment and astrophysical observations, and shows how joint analyses can shed light on the properties of neutron-rich supranuclear matter over the density range probed in neutron stars.

Multi-Messenger Constraints on the Hubble Constant Through Combination of Gravitational Waves, Gamma-Ray Bursts and Kilonovae from Neutron Star Mergers

The simultaneous detection of gravitational waves and light from the binary neutron star merger GW170817 led to independent measurements of distance and redshift, providing a direct estimate of the Hubble constant H0 that does not rely on a cosmic distance ladder, nor assumes a specific cosmological model. By using gravitational waves as “standard sirens”, this approach holds promise to arbitrate the existing tension between the H0 value inferred from the cosmic microwave background and those obtained from local measurements. However, the known degeneracy in the gravitational-wave analysis between distance and inclination of the source led to a H0 value from GW170817 that was not precise enough to resolve the existing tension. In this review, we summarize recent works exploiting the viewing-angle dependence of the electromagnetic signal, namely the associated short gamma-ray burst and kilonova, to constrain the system inclination and improve on H0. We outline the key ingredients of the different methods, summarize the results obtained in the aftermath of GW170817 and discuss the possible systematics introduced by each of these methods.

How Perturbative QCD Constrains the Equation of State at Neutron-Star Densities

We demonstrate in a general and analytic way how high-density information about the equation of state (EOS) of strongly interacting matter obtained using perturbative quantum chromodynamics constrains the same EOS at densities reachable in physical neutron stars. Our approach is based on utilizing the full information of the thermodynamic potentials at the high-density limit together with thermodynamic stability and causality. This requires considering the pressure as a function of chemical potential p(μ) instead of the commonly used pressure as a function of energy density p(ε). The results can be used to propagate the perturbative quantum chromodynamics calculations reliable around 40n_{s} to lower densities in the most conservative way possible. We constrain the EOS starting from only a few times the nuclear saturation density n≳2.2n_{s}, and at n=5n_{s} we exclude at least 65% of otherwise allowed area in the ε-p plane. This provides information complementary to astrophysical observations that should be taken into account in any complete statistical inference study of the EOS. These purely theoretical results are independent of astrophysical neutron-star input, and hence, they can also be used to test theories of modified gravity and beyond the standard model physics in neutron stars.

Constraints on the Maximum Densities of Neutron Stars from Postmerger Gravitational Waves with Third-Generation Observations

Using data from 289 numerical relativity simulations of binary neutron star mergers, we identify, for the first time, a robust quasiuniversal relation connecting the postmerger peak gravitational-wave frequency and the value of the density at the center of the maximum mass nonrotating neutron star. This relation offers a new possibility for precision equation-of-state constraints with next-generation ground-based gravitational-wave interferometers. Mock Einstein Telescope observations of fiducial events indicate that Bayesian inferences can constrain the maximum density to ∼15% (90% credibility level) for a single signal at the minimum sensitivity threshold for a detection. If the postmerger signal is included in a full-spectrum (inspiral-merger-postmerger) analysis of such a signal, the pressure-density function can be tightly constrained up to the maximum density, and the maximum neutron star mass can be measured with an accuracy better than 12% (90% credibility level).

Ranking Love Numbers for the Neutron Star Equation of State: The Need for Third-Generation Detectors

Gravitational-wave measurements of the tidal deformability in neutron-star binary coalescences can be used to infer the still unknown equation of state (EOS) of dense matter above the nuclear saturation density. By employing a Bayesian-ranking test, we quantify the ability of current and future gravitational-wave observations to discriminate among families of nuclear-physics based EOS which differ in particle content and ab initio microscopic calculations. While the constraining power of GW170817 is limited, we show that even twenty coalescences detected by LIGO-Virgo at design sensitivity are not enough to discriminate between EOS with similar softness but distinct microphysics. However, just a single detection with a third-generation detector such as the Einstein Telescope or Cosmic Explorer will rule out several families of EOS with very strong statistical significance and can discriminate among models which feature similar softness, hence, constraining the properties of nuclear matter to unprecedented levels.

Anisotropic Compact Stars in D → 4 Limit of Gauss–Bonnet Gravity

In the frame of Gauss–Bonnet gravity and in the limit of D→4, based on the fact that spherically symmetric solution derived using any of regularization schemes will be the same form as the original theory, we derive a new interior spherically symmetric solution assuming specific forms of the metric potentials that have two constants. Using the junction condition we determine these two constants. By using the data of the star EXO 1785-248, whose mass is M=1.3±0.2M⊙ and radius l=8.849±0.4 km, we calculate the numerical values of these constants, in terms of the dimensionful coupling parameter of the Gauss–Bonnet term, and eventually, we get real values for these constants. In this regard, we show that the components of the energy–momentum tensor have a finite value at the center of the star as well as a smaller value to the surface of the star. Moreover, we show that the equations of the state behave in a non-linear way due to the impact of the Gauss–Bonnet term. Using the Tolman–Oppenheimer–Volkoff equation, the adiabatic index, and stability in the static state we show that the model under consideration is always stable. Finally, the solution of this study is matched with observational data of other pulsars showing satisfactory results.

Neutron Stars with Baryon Number Violation, Probing Dark Sectors

The neutron lifetime anomaly has been used to motivate the introduction of new physics with hidden-sector particles coupled to baryon number, and on which neutron stars provide powerful constraints. Although the neutron lifetime anomaly may eventually prove to be of mundane origin, we use it as motivation for a broader review of the ways that baryon number violation, be it real or apparent, and dark sectors can intertwine and how neutron star observables, both present and future, can constrain them.

Effect of the Nuclear Equation of State and Relativistic Turbulence on Core-Collapse Supernovae

The nuclear Equation of State (EoS) is an important component in the evolution and subsequent explosion of core collapse supernovae. We make a survey of various equations of state that can be found in the literature and analyze their effect on the explosion. To simulate the supernovae, we use the general relativistic spherically-symmetric code GR1D, modified to take into account the effects of three-dimensional turbulence through a new mixing length theory approach (STIR). We show that the viability of the explosion is quite EoS dependent and that the strength of explosions correlate best with the central entropy density right after bounce and the onset of turbulent mixing in the proto-neutron star.

Translating Neutron Star Observations to Nuclear Symmetry Energy via Deep Neural Networks

One of the most significant challenges involved in efforts to understand the equation of state of dense neutron-rich matter is the uncertain density dependence of the nuclear symmetry energy. In particular, the nuclear symmetry energy is still rather poorly constrained, especially at high densities. On the other hand, detailed knowledge of the equation of state is critical for our understanding of many important phenomena in the nuclear terrestrial laboratories and the cosmos. Because of its broad impact, pinning down the density dependence of the nuclear symmetry energy has been a long-standing goal of both nuclear physics and astrophysics. Recent observations of neutron stars, in both electromagnetic and gravitational-wave spectra, have already constrained significantly the nuclear symmetry energy at high densities. The next generation of telescopes and gravitational-wave observatories will provide an unprecedented wealth of detailed observations of neutron stars, which will improve further our knowledge of the density dependence of nuclear symmetry energy, and the underlying equation of state of dense neutron-rich matter. Training deep neural networks to learn a computationally efficient representation of the mapping between astrophysical observables of neutron stars, such as masses, radii, and tidal deformabilities, and the nuclear symmetry energy allows its density dependence to be determined reliably and accurately. In this work, we use a deep learning approach to determine the nuclear symmetry energy as a function of density directly from observational neutron star data. We show, for the first time, that artificial neural networks can precisely reconstruct the nuclear symmetry energy from a set of available neutron star observables, such as masses and radii as measured by, e.g., the NICER mission, or masses and tidal deformabilities as measured by the LIGO/VIRGO/KAGRA gravitational-wave detectors. These results demonstrate the potential of artificial neural networks to reconstruct the symmetry energy and the equation of state directly from neutron star observational data, and emphasize the importance of the deep learning approach in the era of multi-messenger astrophysics.

Early deconfinement of asymptotically conformal color-superconducting quark matter in neutron stars

We present a relativistic density functional approach to color superconducting quark matter that mimics quark confinement by a fast growth of the quasiparticle selfenergy in the confining region. The approach is shown to be equivalent to a chiral model of quark matter with medium dependent couplings. While the (pseudo)scalar sector of the model is fitted to the vacuum phenomenology of quantum chromodynamics, the strength of interaction in the vector and diquark channels is varied in order to provide the best agreement with the observational constraints on the mass-radius relation and tidal deformability of neutron stars modelled with our approach. In order to recover the conformal behavior of quark matter at asymptotically high densities we introduce a medium dependence of the vector and diquark couplings motivated by the non-perturbative gluon exchange. Our analysis signals that the onset of deconfinement to color superconducting quark matter is likely to occur in neutron stars with masses below 1.0 M⊙.

Special point "trains" in the M-R diagram of hybrid stars

We present a systematic investigation of the possible locations for the special point (SP), a unique feature of hybrid neutron stars in the massradius diagram. The study is performed within the two-phase approach where the high-density (quark matter) phase is described by the covariant nonlocal Nambu–Jona-Lasinio (nlNJL) model equation of state (EOS) which is shown to be equivalent to a constant-sound-speed (CSS) EOS. For the nuclear matter phase around saturation density different relativistic density functional EOSs are used: DD2p00, its excluded-volume modification DD2p40 and the hypernuclear EOS DD2Y-T. In the present contribution we apply the Maxwell construction scheme for the deconfinement transition and demonstrate that a simultaneous variation of the vector and diquark coupling constants results in the occurrence of SP "trains" which are invariant against changing the nuclear matter EOS. We propose that the SP train corresponding to a variation of the diquark coupling at constant vector coupling is special since it serves as a lower bound for the line of maximum masses and accessible radii of massive hybrid stars.

The Stochastic Gravitational Wave Background from Magnetars

Magnetars have already been a potential candidate as gravitational wave sources that could be detected by current and future terrestrial as well as ground-based gravitational wave detectors. In this article, we focus on the gravitational wave emission from the distorted rotating neutron stars. The deformation is assumed to be symmetric around an axis that is perpendicular to the rotation axis. The form is applied in the context of a neutron star whose magnetic field has been deformed on its own. By introducing the effects from all magnetars in the Universe, based on various proposed magnetic field configurations, such as poloidal and toroidal, the stochastic gravitational wave background can be generated. We choose to figure out exactly how the observations of the stochastic gravitational wave background should be used to understand much more about physics correlated with the magnetar behavior, based on the restriction on the ellipticity of the magnetar.

A Neutron Star Is Born

A neutron star was first detected as a pulsar in 1967. It is one of the most mysterious compact objects in the universe, with a radius of the order of 10 km and masses that can reach two solar masses. In fact, neutron stars are star remnants, a kind of stellar zombie (they die, but do not disappear). In the last decades, astronomical observations yielded various contraints for neutron star masses, and finally, in 2017, a gravitational wave was detected (GW170817). Its source was identified as the merger of two neutron stars coming from NGC 4993, a galaxy 140 million light years away from us. The very same event was detected in γ-ray, X-ray, UV, IR, radio frequency and even in the optical region of the electromagnetic spectrum, starting the new era of multi-messenger astronomy. To understand and describe neutron stars, an appropriate equation of state that satisfies bulk nuclear matter properties is necessary. GW170817 detection contributed with extra constraints to determine it. On the other hand, magnetars are the same sort of compact object, but bearing much stronger magnetic fields that can reach up to 1015 G on the surface as compared with the usual 1012 G present in ordinary pulsars. While the description of ordinary pulsars is not completely established, describing magnetars poses extra challenges. In this paper, I give an overview on the history of neutron stars and on the development of nuclear models and show how the description of the tiny world of the nuclear physics can help the understanding of the cosmos, especially of the neutron stars.

Nuclear Physics and Astrophysics Constraints on the High Density Matter Equation of State

(1) This review has been written in memory of Steven Moszkowski who unexpectedly passed away in December 2020. It has been inspired by our many years of discussions. Steven’s enthusiasm, drive and determination to understand atomic nuclei in simple terms of basic laws of physics was infectious. He sought the fundamental origin of nuclear forces in free space, and their saturation and modification in nuclear medium. His untimely departure left our job unfinished but his legacy lives on. (2) Focusing on the nuclear force acting in nuclear matter of astrophysical interest and its equation of state (EoS), we take several typical snapshots of evolution of the theory of nuclear forces. We start from original ideas in the 1930s moving through to its overwhelming diversity today. The development is supported by modern observational and terrestrial data and their inference in the multimessenger era, as well as by novel mathematical techniques and computer power. (3) We find that, despite the admirable effort both in theory and measurement, we are facing multiple models dependent on a large number of variable correlated parameters which cannot be constrained by data, which are not yet accurate, nor sensitive enough, to identify the theory closest to reality. The role of microphysics in the theories is severely limited or neglected, mostly deemed to be too difficult to tackle. (4) Taking the EoS of high-density matter as an example, we propose to develop models, based, as much as currently possible, on the microphysics of the nuclear force, with a minimal set of parameters, chosen under clear physical guidance. Still somewhat phenomenological, such models could pave the way to realistic predictions, not tracing the measurement, but leading it.

Progress in Constraining Nuclear Symmetry Energy Using Neutron Star Observables Since GW170817

The density dependence of nuclear symmetry energy is among the most uncertain parts of the Equation of State (EOS) of dense neutron-rich nuclear matter. It is currently poorly known especially at suprasaturation densities partially because of our poor knowledge about isovector nuclear interactions at short distances. Because of its broad impacts on many interesting issues, pinning down the density dependence of nuclear symmetry energy has been a longstanding and shared goal of both astrophysics and nuclear physics. New observational data of neutron stars including their masses, radii, and tidal deformations since GW170817 have helped improve our knowledge about nuclear symmetry energy, especially at high densities. Based on various model analyses of these new data by many people in the nuclear astrophysics community, while our brief review might be incomplete and biased unintentionally, we learned in particular the following: (1) The slope parameter L of nuclear symmetry energy at saturation density ρ0 of nuclear matter from 24 new analyses of neutron star observables was about L≈57.7±19 MeV at a 68% confidence level, consistent with its fiducial value from surveys of over 50 earlier analyses of both terrestrial and astrophysical data within error bars. (2) The curvature Ksym of nuclear symmetry energy at ρ0 from 16 new analyses of neutron star observables was about Ksym≈−107±88 MeV at a 68% confidence level, in very good agreement with the systematics of earlier analyses. (3) The magnitude of nuclear symmetry energy at 2ρ0, i.e., Esym(2ρ0)≈51±13 MeV at a 68% confidence level, was extracted from nine new analyses of neutron star observables, consistent with the results from earlier analyses of heavy-ion reactions and the latest predictions of the state-of-the-art nuclear many-body theories. (4) While the available data from canonical neutron stars did not provide tight constraints on nuclear symmetry energy at densities above about 2ρ0, the lower radius boundary R2.01=12.2 km from NICER’s very recent observation of PSR J0740+6620 of mass 2.08±0.07M⊙ and radius R=12.2–16.3 km at a 68% confidence level set a tight lower limit for nuclear symmetry energy at densities above 2ρ0. (5) Bayesian inferences of nuclear symmetry energy using models encapsulating a first-order hadron–quark phase transition from observables of canonical neutron stars indicated that the phase transition shifted appreciably both L and Ksym to higher values, but with larger uncertainties compared to analyses assuming no such phase transition. (6) The high-density behavior of nuclear symmetry energy significantly affected the minimum frequency necessary to rotationally support GW190814’s secondary component of mass (2.50–2.67) M⊙ as the fastest and most massive pulsar discovered so far. Overall, thanks to the hard work of many people in the astrophysics and nuclear physics community, new data of neutron star observations since the discovery of GW170817 have significantly enriched our knowledge about the symmetry energy of dense neutron-rich nuclear matter.

Combining Electromagnetic and Gravitational-Wave Constraints on Neutron-Star Masses and Radii

We perform a joint Bayesian inference of neutron-star mass and radius constraints based on GW170817, observations of quiescent low-mass x-ray binaries (QLMXBs), photospheric radius expansion x-ray bursting sources, and x-ray timing observations of J0030+0451. With this dataset, the form of the prior distribution still has an impact on the posterior mass-radius curves and equation of state (EOS), but this impact is smaller than recently obtained when considering QLMXBs alone. We analyze the consistency of the electromagnetic data by including an "intrinsic scattering" contribution to the uncertainties, and find only a slight broadening of the posteriors. This suggests that the gravitational-wave and electromagnetic observations of neutron-star structure are providing a consistent picture of the neutron-star mass-radius curve and the EOS.

Multimessenger constraints on the neutron-star equation of state and the Hubble constant

Observations of neutron-star mergers with distinct messengers, including gravitational waves and electromagnetic signals, can be used to study the behavior of matter denser than an atomic nucleus and to measure the expansion rate of the Universe as quantified by the Hubble constant. We performed a joint analysis of the gravitational-wave event GW170817 with its electromagnetic counterparts AT2017gfo and GRB170817A, and the gravitational-wave event GW190425, both originating from neutron-star mergers. We combined these with previous measurements of pulsars using x-ray and radio observations, and nuclear-theory computations using chiral effective field theory, to constrain the neutron-star equation of state. We found that the radius of a 1.4-solar mass neutron star is [Formula: see text] km at 90% confidence and the Hubble constant is [Formula: see text] at 1σ uncertainty.