Aleksi Vuorinen

We perform a model-agnostic Bayesian analysis of the neutron-star-matter equation of state (EoS), using known ab-initio constraints and astrophysical observations to limit its behavior at intermediate densities. Permitting explicit first-order phase transitions allows us to systematically search for twin-star solutions, i.e. the existence of stars degenerate in mass but differing in radius. We find that current observational constraints exclude all but two classes of twin stars. The first is characterized by a first-order transition occurring at a very low density, where the material properties of the system either stay largely intact or move away from the conformal limit. In the second, more interesting class, the discontinuity in the mass-radius curve emerges after a rapid crossover transition at a significantly higher density, with the speed of sound exhibiting two sharp peaks at distinct densities. Since neither class shows clear conformalization upon entering the second branch, the standard twin-star scenario linking the mass-radius discontinuity to deconfinement can be firmly ruled out, while even the remaining solutions -- disfavored by per-mille Bayes factors and in tension with theoretical bounds -- are likely to be excluded in the future.

Knowledge of the pressure of cold and dense quark matter is known to significantly constrain the equation of state of neutron-star matter, a quantity playing a key role in deciphering the stars’ internal structure. In this Letter, we make important progress toward determining the last unknown contribution to the pressure at next-to-next-to-next-to-leading order in the strong coupling constant α s , available through the sum of all four-loop vacuum diagrams of dense quantum chromodynamics. We demonstrate the cancellation of both the covariant gauge parameter and all infrared divergences in the sum, showcasing the effective-field-theory paradigm in action. For the remaining infrared-finite four-loop integrals, we demonstrate the efficacy of the dense loop tree duality method—a novel numerical framework for multiloop calculations in thermal field theory. Together, our results show that completing the next-to-next-to-next-to-leading order pressure of cold quark matter is no longer merely possible but for the first time within reach.

Modern hydrodynamic simulations of core-collapse supernovae and neutron-star mergers require knowledge not only of the equilibrium properties of strongly interacting matter, but also of the system’s response to perturbations, encoded in various transport coefficients. Using perturbative and holographic tools, we derive here an improved weak-coupling and a new strong-coupling result for the most important transport coefficient of unpaired quark matter, its bulk viscosity. These results are combined in a simple analytic pocket formula for the quantity that is rooted in perturbative quantum chromodynamics at high densities but takes into account nonperturbative holographic input at neutron-star densities, where the system is strongly coupled. This expression can be used in the modeling of unpaired quark matter at astrophysically relevant temperatures and densities. Published by the American Physical Society 2024

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.

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.

We study current bounds on strong first-order phase transitions (PTs) along the equation of state (EOS) of dense strongly interacting matter in neutron stars, under the simplifying assumption that on either side of the PT, the EOS can be approximated by a simple polytropic form. We construct a large ensemble of possible EOSs of this form, anchor them to chiral effective field theory calculations at nuclear density and perturbative Quantum Chromodynamics at high densities, and subject them to astrophysical constraints from high-mass pulsars and gravitational-wave observations. Within this setup, we find that a PT permits neutron-star solutions with larger radii, but only if the transition begins below twice nuclear saturation density. We also identify a large parameter space of allowed PTs currently unexplored by numerical-relativity studies. Additionally, we locate a small region of parameter space allowing twin-star solutions, though we find them to only marginally pass the current astrophysical constraints. Finally, we find that sizeable cores of high-density matter beyond the PT may be located in the centers of some stable neutron stars, primarily those with larger masses.