Identifying a First-Order Phase Transition in Neutron-Star Mergers through Gravitational Waves

Impact of Bulk Viscosity on the Postmerger Gravitational-Wave Signal from Merging Neutron Stars

In the violent postmerger of binary neutron-star mergers strong oscillations are present that impact the emitted gravitational-wave (GW) signal. The frequencies, temperatures, and densities involved in these oscillations allow for violations of the chemical equilibrium promoted by weak interactions, thus leading to a nonzero bulk viscosity that can impact dynamics and GW signals. We present the first simulations of binary neutron-star mergers employing the self-consistent and second-order formulation of the equations of relativistic hydrodynamics for dissipative fluids proposed by Müller, Israel, and Stewart. With the spirit of obtaining a first assessment of the impact of bulk viscosity on the structure and radiative efficiency of the merger remnant we adopt a simplified but realistic approach for the viscosity, which we assume to be determined by direct and modified Urca reactions and hence to vary within the stars. At the same time, to compensate for the lack of a precise knowledge about the strength of bulk viscosity, we explore the possible behaviors by considering three different scenarios of low, medium, and high bulk viscosity. In this way, we find that large values of the bulk viscosities damp the collision-and-bounce oscillations that characterize the dynamics of the stellar cores right after the merger. At the same time, large viscosities tend to preserve the m=2 deformations in the remnant, thus leading to a comparatively more efficient GW emission and to changes in the postmerger spectrum that can be up to 100 Hz in the case of the most extreme configurations. Overall, our self-consistent results indicate that bulk viscosity increases the energy radiated in GWs soon after the merger by ≲2% in the (realistic) scenario of small viscosity, and by ≲30% in the (unrealistic) scenario of large viscosity.

Achieving the Fundamental Quantum Limit of Linear Waveform Estimation

Sensing a classical signal using a linear quantum device is a pervasive application of quantum-enhanced measurement. The fundamental precision limits of linear waveform estimation, however, are not fully understood. In certain cases, there is an unexplained gap between the known waveform-estimation quantum Cramér-Rao bound and the optimal sensitivity from quadrature measurement of the outgoing mode from the device. We resolve this gap by establishing the fundamental precision limit, the waveform-estimation Holevo Cramér-Rao bound, and how to achieve it using a nonstationary measurement. We apply our results to detuned gravitational-wave interferometry to accelerate the search for postmerger remnants from binary neutron-star mergers. If we have an unequal weighting between estimating the signal's power and phase, then we propose how to further improve the signal-to-noise ratio by a factor of sqrt[2] using this nonstationary measurement.

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.

Degeneracy in the Inference of Phase Transitions in the Neutron Star Equation of State from Gravitational Wave Data

Gravitational wave (GW) detections of binary neutron star inspirals will be crucial for constraining the dense matter equation of state (EOS). We demonstrate a new degeneracy in the mapping from tidal deformability data to the EOS, which occurs for models with strong phase transitions. We find that there exists a new family of EOS with phase transitions that set in at different densities and that predict neutron star radii that differ by up to ∼500  m but that produce nearly identical tidal deformabilities for all neutron star masses. Next-generation GW detectors and advances in nuclear theory may be needed to resolve this degeneracy.

Gravitational Wave Signal for Quark Matter with Realistic Phase Transition

The cores of neutron stars (NSs) near the maximum mass can realize a transitional change to quark matter (QM). Gravitational waves from binary NS mergers are expected to convey information about the equation of state (EOS) sensitive to the QM transition. Here, we present the first results of gravitational wave simulation with the realistic EOS that is consistent with ab initio approaches of χEFT and pQCD and is assumed to go through smooth crossover. We compare them to results obtained with another EOS with a first-order hadron-quark phase transition. Our results suggest that early collapse to a black hole in the post-merger stage after NS merger robustly signifies softening of the EOS associated with the QM onset in the crossover scenario. The nature of the hadron-quark phase transition can be further constrained by the condition that electromagnetic counterparts should be energized by the material left outside the remnant black hole.

Merger and Postmerger of Binary Neutron Stars with a Quark-Hadron Crossover Equation of State

Fully general-relativistic binary-neutron-star (BNS) merger simulations with quark-hadron crossover (QHC) equations of state (EOS) are studied for the first time. In contrast to EOS with purely hadronic matter or with a first-order quark-hadron phase transition (1PT), in the transition region QHC EOS show a peak in sound speed and thus a stiffening. We study the effects of such stiffening in the merger and postmerger gravitational (GW) signals. Through simulations in the binary-mass range 2.5<M/M_{⊙}<2.75, characteristic differences due to different EOS appear in the frequency of the main peak of the postmerger GW spectrum (f_{2}), extracted through Bayesian inference. In particular, we found that (i) for lower-mass binaries, since the maximum baryon number density (n_{max}) after the merger stays below 3-4 times the nuclear-matter density (n_{0}), the characteristic stiffening of the QHC models in that density range results in a lower f_{2} than that computed for the underlying hadronic EOS and thus also than that for EOS with a 1PT; (ii) for higher-mass binaries, where n_{max} may exceed 4-5n_{0} depending on the EOS model, whether f_{2} in QHC models is higher or lower than that in the underlying hadronic model depends on the height of the sound-speed peak. Comparing the values of f_{2} for different EOS and BNS masses gives important clues on how to discriminate different types of quark dynamics in the high-density end of EOS and is relevant to future kilohertz GW observations with third-generation GW detectors.

Gravitational Waves from Accretion-Induced Descalarization in Massive Scalar-Tensor Theory

Many classes of extended scalar-tensor theories predict that dynamical instabilities can take place at high energies, leading to the formation of scalarized neutron stars. Depending on the theory parameters, stars in a scalarized state can form a solution-space branch that shares a lot of similarities with the so-called mass twins in general relativity appearing for equations of state containing first-order phase transitions. Members of this scalarized branch have a lower maximum mass and central energy density compared to Einstein ones. In such cases, a scalarized star could potentially overaccrete beyond the critical mass limit, thus triggering a gravitational phase transition where the star sheds its scalar hair and migrates over to its nonscalarized counterpart. Such an event resembles, but is distinct from, a nuclear or thermodynamic phase transition. We dynamically track a gravitational transition by first constructing hydrostatic, scalarized equilibria for realistic equations of state, and then allowing additional material to fall onto the stellar surface. The resulting bursts of monopolar radiation are dispersively stretched to form a quasicontinuous signal that persists for decades, carrying strains of order ≳10^{-22}  (kpc/L)^{3/2} Hz^{-1/2} at frequencies of ≲300  Hz, detectable with the existing interferometer network out to distances of L≲10  kpc, and out to a few hundred kpc with the inclusion of the Einstein Telescope. Electromagnetic signatures of such events, involving gamma-ray and neutrino bursts, are also considered.

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).

QED Mesons, the QED Neutron, and the Dark Matter

Schwinger’s boson solution for massless fermions in QED in 1+1D has been applied and generalized to quarks interacting in QED and QCD interactions, leading to stable and confined open-string QED and QCD boson excitations of the quark-QCD-QED system in 1+1D. Just as the open-string QCD excitations in 1+1D can be the idealization of QCD mesons with a flux tube in 3+1D, so the open-string QED excitations in 1+1D may likewise be the idealization of QED mesons with masses in the tens of MeV region, corresponding possibly to the anomalous X17 and E38 particles observed recently. A further search for bound states of quarks interacting in the QED interaction alone leads to the examination on the stability of the QED neutron, consisting of two d quarks and one u quark. Theoretically, the QED neutron has been found to be stable and estimated to have a mass of 44.5 MeV, whereas the analogous QED proton is unstable, leading to a long-lived QED neutron that may be a good candidate for the dark matter.

Triple nuclear collisions – a new method to explore the matter properties under new extreme conditions

We suggest to explore an entirely new method to experimentally and theoretically study the phase diagram of strongly interacting matter based on the triple nuclear collisions (TNC).We simulated the TNC using the UrQMD 3.4 model at the beam center of- mass collision energies √SNN = 200 GeV and √SNN = 2.76 TeV. It is found that in the most central and simultaneous TNC the initial baryonic charge density is about 3 times higher than the one achieved in the usual binary nuclear collisions at the same energies. As a consequence, the production of protons and Λ-hyperons is increased by a factor of 2 and 1.5, respectively. Using the MIT Bag model equation we study the evolution of the central cell in TNC and demonstrate that for the top RHIC energy of collision the baryonic chemical potential is 2-2.5 times larger than the one achieved in the binary nuclear collision at the same time of reaction. Based on these estimates, we show that TNC offers an entirely new possibility to study the QCD phase diagram at very high baryonic charge densities.

Locating the special point of hybrid neutron stars

The special point is a feature unique to models of hybrid neutron stars. It represents a location on their mass–radius sequences that is insensitive to the phase transition density. We consider hybrid neutron stars with a core of deconfined quark matter that obeys a constant–sound–speed (CSS) equation of state model and provide a fit formula for the coordinates of the special point as functions of the squared sound speed (cs2) and pressure scale (A) parameters. Using the special point mass as a proxy for the maximum mass of the hybrid stars we derive limits for the CSS model parameters based on the recent NICER constraint on mass and radius of pulsar PSR J0740+6620, 0.36 < Cs min2 < 0.43 and 80 < A[MeV/fm3] < 160. The upper limit for the maximum mass of hybrid stars depends on the upper limit for cs2 so that choosing cs,max2 = 0.6 results in Mmax < 2.7 M⊙, within the mass range of GW190814.

NICER view on holographic QCD

The holographic models for dense QCD matter work surprisingly well. A general implication seems that the deconfinement phase transition dictates the maximum mass of neutron stars. The nuclear matter phase turns out to be rather stiff which, if continuously merged with nuclear matter models based on effective field theories, leads to the conclusion that neutron stars do not have quark matter cores in the light of all current astrophysical data. We comment that as the perturbative QCD results are in stark contrast with strong coupling results, any future simulations of neutron star mergers incorporating corrections beyond ideal fluid should proceed cautiously. For this purpose, we provide a model which treats nuclear and quark matter phases in a unified framework at strong coupling.

Gravitational wave signatures of phase transition from hadronic to quark matter in isolated neutron stars and binaries

The fundamental constituent of matter at high temperature and density has intrigued physicists for quite some time. Recent results from heavy-ion colliders have enriched the Quantum Chromodynamics phase diagram at high temperatures and low baryon density. However, the phase at low temperatures and finite (mostly intermediate) baryon density remain unexplored. Theoretical Quantum Chromodynamics calculation predicts phase transition from hadronic matter to quark matter at such densities. Presently, the best laboratories available to probe such densities lie at the core of neutron stars. Recent results of how such phase transition signatures can be probed using gravitational waves both in isolated neutron stars and neutron star in binaries. The isolated neutron star would probe the very low-temperature regime, whereas neutron stars in binaries would probe finite baryon density in the intermediate temperature regime. We would also discuss whether the gravitational wave signature of such phase transition is unique and the detector specification needed to detect such signals.

Exploring the Phase Diagram of V-QCD with Neutron Star Merger Simulations

Determining the phase structure of Quantum Chromodynamics (QCD) and its Equation of State (EOS) at densities and temperatures realized inside neutron stars and their mergers is a long-standing open problem. The holographic V-QCD framework provides a model for the EOS of dense and hot QCD, which describes the deconfinement phase transition between a dense baryonic and a quark matter phase. We use this model in fully general relativistic hydrodynamic (GRHD) simulations to study the formation of quark matter and the emitted gravitational wave signal of binary systems that are similar to the first ever observed neutron star merger event GW170817.

Neutron Star Mergers and the Quark Matter Equation of State

As neutron stars merge they can approach very high nuclear density. Here, we summarized recent results for the evolution and gravitational wave emission from binary-neutron star mergers using a a variety of nuclear equations of state with and without a crossover transition to quark matter. We discuss how the late time gravitational wave emission from binary neutron star mergers may possibly reveal the existence of a crossover transition to quark matter.

Phase Conversions in Neutron Stars: Implications for Stellar Stability and Gravitational Wave Astrophysics

We review the properties of hybrid stars with a quark matter core and a hadronic mantle, focusing on the role of key micro-physical properties such as the quark/hadron surface and curvature tensions and the conversion speed at the interface between both phases. We summarize the results of works that have determined the surface and curvature tensions from microscopic calculations. If these quantities are large enough, mixed phases are energetically suppressed and the quark core would be separated from the hadronic mantle by a sharp interface. If the conversion speed at the interface is slow, a new class of dynamically stable hybrid objects is possible. Densities tens of times larger than the nuclear saturation density can be attained at the center of these objects. We discuss possible formation mechanisms for the new class of hybrid stars and smoking guns for their observational identification.

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.

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.

Transport in Strongly Coupled Quark Matter

Motivated by the possible presence of deconfined quark matter in neutron stars and their mergers and the important role of transport phenomena in these systems, we perform the first-ever systematic study of different viscosities and conductivities of dense quark matter using the gauge/gravity duality. Using the V-QCD model, we arrive at results that are in qualitative disagreement with the predictions of perturbation theory, which highlights the differing transport properties of the system at weak and strong coupling and calls for caution in the use of the perturbative results in neutron star applications.

Neutron star merger remnants

Binary neutron star mergers observations are a unique way to constrain fundamental physics and astrophysics at the extreme. The interpretation of gravitational-wave events and their electromagnetic counterparts crucially relies on general-relativistic models of the merger remnants. Quantitative models can be obtained only by means of numerical relativity simulations in 3 + 1 dimensions including detailed input physics for the nuclear matter, electromagnetic and weak interactions. This review summarizes the current understanding of merger remnants focusing on some of the aspects that are relevant for multimessenger observations.

Equation of State Constraints from the Threshold Binary Mass for Prompt Collapse of Neutron Star Mergers

Using hydrodynamical simulations for a large set of high-density matter equations of state (EOSs), we systematically determine the threshold mass M_{thres} for prompt black-hole formation in equal-mass and asymmetric neutron star (NS) mergers. We devise the so far most direct, general, and accurate method to determine the unknown maximum mass of nonrotating NSs from merger observations revealing M_{thres}. Considering hybrid EOSs with hadron-quark phase transition, we identify a new, observable signature of quark matter in NS mergers. Furthermore, our findings have direct applications in gravitational wave searches, kilonova interpretations, and multimessenger constraints on NS properties.

Gravitational-Wave Signature of a First-Order Quantum Chromodynamics Phase Transition in Core-Collapse Supernovae

A first-order quantum chromodynamics (QCD) phase transition (PT) may take place in the protocompact star (PCS) produced by a core-collapse supernova (CCSN). In this work, we study the consequences of such a PT in a nonrotating CCSN with axisymmetric hydrodynamic simulations. We find that the PT leads to the collapse of the PCS and results in a loud burst of gravitational waves (GWs). The amplitude of this GW burst is ∼30 times larger than the postbounce GW signal normally found for nonrotating CCSN. It shows a broad peak at high frequencies (∼2500-4000  Hz) in the spectrum, has a duration of ≲5  ms, and carries ∼3 orders of magnitude more energy than the other episodes. Also, the peak frequency of the PCS oscillation increases dramatically after the PT-induced collapse. In addition to a second neutrino burst, the GW signal, if detected by the ground-based GW detectors, is decisive evidence of the first-order QCD PT inside CCSNe and provides key information about the structure and dynamics of the PCS.

Was GW170817 a Canonical Neutron Star Merger? Bayesian Analysis with a Third Family of Compact Stars

We investigate the possibility that GW170817 was not the merger of two conventional neutron stars (NS), but involved at least one if not two hybrid stars with a quark matter core that might even belong to a third family of compact stars. To this end, we develop a Bayesian analysis method for selecting the most probable equation of state (EoS) under a set of constraints from compact star physics, which now also include the tidal deformability from GW170817 and the first result for the mass and radius determination for PSR J0030+0451 by the NICER Collaboration. We apply this method for the first time to a two-parameter family of hybrid EoS based on the DD2 model with nucleonic excluded volume for hadronic matter and the color superconducting generalized nlNJL model for quark matter. The model has a variable onset density for deconfinement and can mimic the effects of pasta phases with the possibility of producing a third family of hybrid stars in the mass-radius diagram. The main findings of this study are that: (1) the presence of multiple configurations for a given mass (twins or even triples) corresponds to a set of disconnected lines in the Λ 1 – Λ 2 diagram of tidal deformabilities for binary mergers, so that merger events from the same mass range may result in a probability landscape with different peak positions; (2) the Bayesian analysis with the above observational constraints favors an early onset of the deconfinement transition, at masses of M onset ≤ 0.8 M ⊙ with an M–R relationship that in the range of observed neutron star masses is almost indistinguishable from that of a soft hadronic Akmal, Pandharipande, and Ravenhall (APR) EoS; (3) a few, yet fictitious measurements of the NICER experiment two times more accurate than the present value and a different mass and radius that would change the posterior likelihood so that hybrid EoS with a phase transition onset in the range M onset = 1.1–1.6 M ⊙ would be favored.

Postmerger Gravitational-Wave Signatures of Phase Transitions in Binary Mergers

With the first detection of gravitational waves from a binary system of neutron stars GW170817, a new window was opened to study the properties of matter at and above nuclear-saturation density. Reaching densities a few times that of nuclear matter and temperatures up to 100 MeV, such mergers also represent potential sites for a phase transition (PT) from confined hadronic matter to deconfined quark matter. While the lack of a postmerger signal in GW170817 has prevented us from assessing experimentally this scenario, two theoretical studies have explored the postmerger gravitational-wave signatures of PTs in mergers of a binary system of neutron stars. We here extend and complete the picture by presenting a novel signature of the occurrence of a PT. More specifically, using fully general-relativistic hydrodynamic simulations and employing a suitably constructed equation of state that includes a PT, we present the occurrence of a "delayed PT," i.e., a PT that develops only some time after the merger and produces a metastable object with a quark-matter core, i.e., a hypermassive hybrid star. Because in this scenario, the postmerger signal exhibits two distinct fundamental gravitational-wave frequencies-before and after the PT-the associated signature promises to be the strongest and cleanest among those considered so far, and one of the best signatures of the production of quark matter in the present Universe.

Impact of the Nuclear Equation of State on the Stability of Hybrid Neutron Stars

We construct a set of equations of state (EoS) of dense and hot matter with a 1st order phase transition from a hadronic system to a deconfined quark matter state. In this two-phase approach, hadrons are described using the relativistic mean field theory with different parametrisations and the deconfined quark phase is modeled using vBag, a bag–type model extended to include vector interactions as well as a simultaneous onset of chiral symmetry restoration and deconfinement. This feature results in a non–trivial connection between the hadron and quark EoS, modifying the quark phase beyond its onset density. We find that this unique property has an impact on the predicted hybrid (quark core) neutron star mass–radius relations.

Detecting the Hadron-Quark Phase Transition with Gravitational Waves

The long-awaited detection of a gravitational wave from the merger of a binary neutron star in August 2017 (GW170817) marks the beginning of the new field of multi-messenger gravitational wave astronomy. By exploiting the extracted tidal deformations of the two neutron stars from the late inspiral phase of GW170817, it is now possible to constrain several global properties of the equation of state of neutron star matter. However, the most interesting part of the high density and temperature regime of the equation of state is solely imprinted in the post-merger gravitational wave emission from the remnant hypermassive/supramassive neutron star. This regime was not observed in GW170817, but will possibly be detected in forthcoming events within the current observing run of the LIGO/VIRGO collaboration. Numerous numerical-relativity simulations of merging neutron star binaries have been performed during the last decades, and the emitted gravitational wave profiles and the interior structure of the generated remnants have been analysed in detail. The consequences of a potential appearance of a hadron-quark phase transition in the interior region of the produced hypermassive neutron star and the evolution of its underlying matter in the phase diagram of quantum cromo dynamics will be in the focus of this article. It will be shown that the different density/temperature regions of the equation of state can be severely constrained by a measurement of the spectral properties of the emitted post-merger gravitational wave signal from a future binary compact star merger event.

Neutron-Star-Merger Equation of State

In this work, we discuss the dense matter equation of state (EOS) for the extreme range of conditions encountered in neutron stars and their mergers. The calculation of the properties of such an EOS involves modeling different degrees of freedom (such as nuclei, nucleons, hyperons, and quarks), taking into account different symmetries, and including finite density and temperature effects in a thermodynamically consistent manner. We begin by addressing subnuclear matter consisting of nucleons and a small admixture of light nuclei in the context of the excluded volume approach. We then turn our attention to supranuclear homogeneous matter as described by the Chiral Mean Field (CMF) formalism. Finally, we present results from realistic neutron-star-merger simulations performed using the CMF model that predict signatures for deconfinement to quark matter in gravitational wave signals.

Questions Related to the Equation of State of High-Density Matter

Astronomical data about neutron stars can be combined with laboratory nuclear data to give us a strong base from which to infer the equation of state of cold catalyzed matter beyond nuclear density. However, the nuclear and astrophysical communities are largely distinct; each has their own methods, which means that there is often imperfect communication between the communities regarding caveats about claimed measurements and constraints. Here we present a brief summary from one astronomer’s perspective of relevant observations of neutron stars, with warnings as appropriate, followed by a set of questions that are intended to help enhance the dialog between nuclear physicists and astrophysicists.

Signatures of Quark-Hadron Phase Transitions in General-Relativistic Neutron-Star Mergers

Merging binaries of neutron-stars are not only strong sources of gravitational waves, but also have the potential of revealing states of matter at densities and temperatures not accessible in laboratories. A crucial and long-standing question in this context is whether quarks are deconfined as a result of the dramatic increase in density and temperature following the merger. We present the first fully general-relativistic simulations of merging neutron-stars including quarks at finite temperatures that can be switched off consistently in the equation of state. Within our approach, we can determine clearly what signatures a quark-hadron phase transition would leave in the gravitational-wave signal. We show that if after the merger the conditions are met for a phase transition to take place at several times nuclear saturation density, they would lead to a postmerger signal considerably different from the one expected from the inspiral, that can only probe the hadronic part of the equations of state, and to an anticipated collapse of the merged object. We also show that the phase transition leads to a very hot and dense quark core that, when it collapses to a black hole, produces a ringdown signal different from the hadronic one. Finally, in analogy with what is done in heavy-ion collisions, we use the evolution of the temperature and density in the merger remnant to illustrate the properties of the phase transition in a QCD phase diagram.