Constraining neutron-star matter with microscopic and macroscopic collisions

Nuclear structure opportunities with GeV radioactive beams at FAIR

The Facility for Antiproton and Ion Research (FAIR) is in its final construction stage next to the campus of the Gesellschaft für Schwerionenforschung Helmholtzzentrum for heavy-ion research in Darmstadt, Germany. Once it starts its operation, it will be the main nuclear physics research facility in many basic sciences and their applications in Europe for the coming decades. Owing to the ability of the new fragment separator, Super-FRagment Separator, to produce high-intensity radioactive ion beams in the energy range up to about 2 GeV/nucleon, these can be used in various nuclear reactions. This opens a unique opportunity for various nuclear structure studies across a range of fields and scales: from low-energy physics via the investigation of multi-neutron systems and halos to high-density nuclear matter and the equation of state, following heavy-ion collisions, fission and study of short-range correlations in nuclei and hypernuclei. The newly developed reactions with relativistic radioactive beams (R3B) set up at FAIR would be the most suitable and versatile for such studies. An overview of highlighted physics cases foreseen at R3B is given, along with possible future opportunities, at FAIR. This article is part of the theme issue 'The liminal position of Nuclear Physics: from hadrons to neutron stars'.

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.

Signatures of Prethermalization in a Quenched Cavity-Mediated Long-Range Interacting Fermi Gas

The coupling of ultracold quantum gases to an optical cavity provides an ideal system for studying the novel long-range interacting nonequilibrium dynamics. Here we report an experimental observation of the out-of-equilibrium dynamics of a degenerate Fermi gas in the cavity after quenching the pump strength over a superradiant quantum phase transition. The relaxation dynamics exhibits impressively different stages of a delay, violent relaxation, long-lifetime prethermalization, and slowly final thermalization due to the photon-mediated long-range interaction with dissipation. Importantly, we reveal that the lifetime of the system stayed on the prethermalization exhibits the superlinear scaling of the atom number. Furthermore, we show that the backaction of the superradiant cavity field on the gas causes the exchange of atoms between the normal and superradiant state in the early evolution and then induces the prethermalization. This work opens an avenue to explore complex nonequilibrium dynamics of the dissipatively long-range interacting Fermi gases.

QCD Equation of State of Dense Nuclear Matter from a Bayesian Analysis of Heavy-Ion Collision Data

Bayesian methods are used to constrain the density dependence of the QCD equation of state (EOS) for dense nuclear matter using the data of mean transverse kinetic energy and elliptic flow of protons from heavy ion collisions (HICs), in the beam energy range sqrt[s_{NN}]=2-10  GeV. The analysis yields tight constraints on the density dependent EOS up to 4 times the nuclear saturation density. The extracted EOS yields good agreement with other observables measured in HIC experiments and constraints from astrophysical observations both of which were not used in the inference. The sensitivity of inference to the choice of observables is also discussed.

Gravitational Waves from a Core g Mode in Supernovae as Probes of the High-Density Equation of State

Using relativistic supernova simulations of massive progenitor stars with a quark-hadron equation of state (EOS) and a purely hadronic EOS, we identify a distinctive feature in the gravitational-wave signal that originates from a buoyancy-driven mode (g mode) below the proto-neutron star convection zone. The mode frequency lies in the range 200≲f≲800  Hz and decreases with time. As the mode lives in the core of the proto-neutron star, its frequency and power are highly sensitive to the EOS, in particular the sound speed around twice saturation density.

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.

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.

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.