Chiral Three-Nucleon Interactions in Light Nuclei, Neutron-α Scattering, and Neutron Matter

Revisiting the Helium Isotope-Shift Puzzle with Improved Uncertainties from Nuclear Structure Corrections

Measurements of the difference between the squared charge radii of the helion (^{3}He nucleus) and the α particle (^{4}He nucleus) have been characterized by longstanding tensions recently spotlighted in the 3.6σ discrepancy of the extractions from ordinary atoms versus those from muonic atoms [Karsten Schuhmann et al., arXiv:2305.11679]. Here, we present a novel analysis of uncertainties in nuclear structure corrections that must be supplied by theory to enable the extraction of the difference in radii from spectroscopic experiments. We use modern Bayesian inference techniques to quantify uncertainties stemming from the truncation of the chiral effective field theory expansion of the nuclear force for both muonic and ordinary atoms. With the new nuclear structure input, the helium isotope-shift puzzle cannot be explained, rather, it is reinforced to a 4σ discrepancy.

Distilling the Essential Elements of Nuclear Binding via Neural-Network Quantum States

To distill the essential elements of nuclear binding, we seek the simplest Hamiltonian capable of modeling atomic nuclei with percent-level accuracy. A critical aspect of this endeavor consists of accurately solving the quantum many-body problem without incurring an exponential computing cost with the number of nucleons. We address this challenge by leveraging a variational Monte Carlo method based on a highly expressive neural-network quantum state ansatz. In addition to computing binding energies and charge radii of nuclei with up to A=20 nucleons, by evaluating their magnetic moments, we demonstrate that neural-network quantum states are able to correctly capture the self-emerging nuclear shell structure. To this end, we introduce a novel computational protocol based on adding an external magnetic field to the nuclear Hamiltonian, which allows the neural network to learn the preferred polarization of the nucleus within the given magnetic field.

Electromagnetic properties of nuclei from first principles: a case for synergies between experiment and theory

One of the overarching goals in nuclear science is to understand how the nuclear chart emerges from the underlying fundamental interactions. The description of the structure of nuclei from first principles, using ab initio methods for the solution of the many-nucleon problem with inputs from chiral effective field theory, has advanced dramatically over the past two decades. We present an overview over the available ab initio tools with a specific emphasis on electromagnetic observables, such as multipole moments and transition strengths. These observables still pose a challenge for ab initio theory and are one of the most exciting domains to exploit synergies with modern experiments. Precise experimental data are vital for the validation of the theory predictions and the refinement of ab initio methods. We discuss some of the past and future experimental efforts highlighting these synergies. This article is part of the theme issue 'The liminal position of Nuclear Physics: from hadrons to neutron stars'.

Structure Factors for Hot Neutron Matter from Ab Initio Lattice Simulations with High-Fidelity Chiral Interactions

We present the first ab initio lattice calculations of spin and density correlations in hot neutron matter using high-fidelity interactions at next-to-next-to-next-to-leading order in chiral effective field theory. These correlations have a large impact on neutrino heating and shock revival in core-collapse supernovae and are encapsulated in functions called structure factors. Unfortunately, calculations of structure factors using high-fidelity chiral interactions were well out of reach using existing computational methods. In this Letter, we solve the problem using a computational approach called the rank-one operator (RO) method. The RO method is a general technique with broad applications to simulations of fermionic many-body systems. It solves the problem of exponential scaling of computational effort when using perturbation theory for higher-body operators and higher-order corrections. Using the RO method, we compute the vector and axial static structure factors for hot neutron matter as a function of temperature and density. The ab initio lattice results are in good agreement with virial expansion calculations at low densities but are more reliable at higher densities. Random phase approximation codes used to estimate neutrino opacity in core-collapse supernovae simulations can now be calibrated with ab initio lattice calculations.

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.

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.

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.

Trends of Neutron Skins and Radii of Mirror Nuclei from First Principles

The neutron skin of atomic nuclei impacts the structure of neutron-rich nuclei, the equation of state of nucleonic matter, and the size of neutron stars. Here we predict the neutron skin of selected light- and medium-mass nuclei using coupled-cluster theory and the auxiliary field diffusion Monte Carlo method with two- and three-nucleon forces from chiral effective field theory. We find a linear correlation between the neutron skin and the isospin asymmetry in agreement with the liquid-drop model and compare with data. We also extract the linear relationship that describes the difference between neutron and proton radii of mirror nuclei and quantify the effect of charge symmetry breaking terms in the nuclear Hamiltonian. Our results for the mirror-difference charge radii and binding energies per nucleon agree with existing data.

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.

The Equation of State of Neutron-Rich Matter at Fourth Order of Chiral Effective Field Theory and the Radius of a Medium-Mass Neutron Star

We report neutron star predictions based on our most recent equations of state. These are derived from chiral effective field theory, which allows for a systematic development of nuclear forces, order by order. We utilize high-quality two-nucleon interactions and include all three-nucleon forces up to fourth order in the chiral expansion. Our ab initio predictions are restricted to the domain of applicability of chiral effective field theory. However, stellar matter in the interior of neutron stars can be up to several times denser than normal nuclear matter at saturation, and its composition is essentially unknown. Following established practices, we extend our microscopic predictions to higher densities matching piecewise polytropes. The radius of the average-size neutron star, about 1.4 solar masses, is sensitive to the pressure at normal densities, and thus it is suitable to constrain ab initio theories of the equation of state. For this reason, we focus on the radius of medium-mass stars. We compare our results with other theoretical predictions and recent constraints.

The 1S0 Pairing Gap in Neutron Matter

We report ab initio calculations of the S wave pairing gap in neutron matter calculated using realistic nuclear Hamiltonians that include two- and three-body interactions. We use a trial state, properly optimized to capture the essential pairing correlations, from which we extract ground state properties by means of auxiliary field diffusion Monte Carlo simulations. We extrapolate our results to the thermodynamic limit by studying the finite-size effects in the symmetry-restored projected Bardeen-Cooper-Schrieffer (PBCS) theory and compare our results to other ab initio studies done in the past. Our quantum Monte Carlo results for the pairing gap show a modest suppression with respect to the mean-field BCS values. These results can be connected to cold atom experiments, via the unitarity regime where fermionic superfluidity assumes a unified description, and they are important in the prediction of thermal properties and the cooling of neutron stars.

Astrophysical Constraints on the Symmetry Energy and the Neutron Skin of ^{208}Pb with Minimal Modeling Assumptions

The symmetry energy and its density dependence are crucial inputs for many nuclear physics and astrophysics applications, as they determine properties ranging from the neutron-skin thickness of nuclei to the crust thickness and the radius of neutron stars. Recently, PREX-II reported a value of 0.283±0.071  fm for the neutron-skin thickness of ^{208}Pb, implying a slope parameter L=106±37  MeV, larger than most ranges obtained from microscopic calculations and other nuclear experiments. We use a nonparametric equation of state representation based on Gaussian processes to constrain the symmetry energy S_{0}, L, and R_{skin}^{^{208}Pb} directly from observations of neutron stars with minimal modeling assumptions. The resulting astrophysical constraints from heavy pulsar masses, LIGO/Virgo, and NICER clearly favor smaller values of the neutron skin and L, as well as negative symmetry incompressibilities. Combining astrophysical data with PREX-II and chiral effective field theory constraints yields S_{0}=33.0_{-1.8}^{+2.0}  MeV, L=53_{-15}^{+14}  MeV, and R_{skin}^{^{208}Pb}=0.17_{-0.04}^{+0.04}  fm.

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.

Symmetric Nuclear Matter from the Strong Interaction

We study the equation of state of symmetric nuclear matter at zero temperature over a wide range of densities using two complementary theoretical approaches. At low densities, up to twice nuclear saturation density, we compute the energy per particle based on modern nucleon-nucleon and three-nucleon interactions derived within chiral effective field theory. For higher densities, we derive for the first time constraints in a Fierz-complete setting directly based on quantum chromodynamics using functional renormalization group techniques. We find remarkable consistency of the results obtained from both approaches as they come together in density and the natural emergence of a maximum in the speed of sound c_{S} at supranuclear densities. The presence of this maximum appears tightly connected to the formation of a diquark gap. Notably, this maximum is observed to exceed the asymptotic value c_{S}^{2}=1/3 while its exact position in terms of the density cannot yet be determined conclusively.

Equation of State Effects in Core-Collapse Supernovae

We investigate the impact of different properties of the nuclear equation of state in core-collapse supernovae, with a focus on the proto-neutron-star contraction and its impact on the shock evolution. To this end, we introduce a range of equations of state that vary the nucleon effective mass, incompressibility, symmetry energy, and nuclear saturation point. This allows us to point to the different effects in changing these properties from the Lattimer and Swesty to the Shen et al. equations of state, the two most commonly used equations of state in simulations. In particular, we trace the contraction behavior to the effective mass, which determines the thermal nucleonic contributions to the equation of state. Larger effective masses lead to lower pressures at nuclear densities and a lower thermal index. This results in a more rapid contraction of the proto-neutron star and consequently higher neutrino energies, which aids the shock evolution to a faster explosion.

What can neutron stars reveal about the equation of state of dense matter?

Neutron stars are astrophysical objects of extremes, reaching the highest densities we can observe in the cosmos, and probing matter under conditions that cannot be recreated in terrestrial experiments. In August 2017, the first neutron-star merger has been observed, which provided compelling evidence that these events are an important site for r-process nucleosynthesis. Furthermore, the gravitational-wave signal of such events might shed light upon the nature of strongly interacting matter in the neutron-star core. To understand these remarkable events, reliable nuclear physics input is essential. In this contribution, I explain how to use chiral effective field theory and advanced many-body methods to provide a consistent and systematic approach to strongly inter- acting systems from nuclei to neutron stars with controlled theoretical uncertainties. I will discuss recent results for the equation of state relevant for the nuclear astrophysics of neutron stars and neutron-star mergers.

Harmonic Potential Theorem: Extension to Spin-, Velocity-, and Density-Dependent Interactions

One of the few exact results for the description of the time evolution of an inhomogeneous, interacting many-particle system is given by the harmonic potential theorem (HPT). The relevance of this theorem is that it sets a tight constraint on time-dependent many-body approximations. In this contribution, we show that the original formulation of the HPT is valid also for the case of spin-, velocity-, and density-dependent interactions. This result is completely general and relevant, among the rest, for nuclear structure theory both in the case of ab initio and of more phenomenological approaches. As an example, we report on a numerical implementation by testing the small-amplitude limit of the time-dependent Hartree-Fock-also known as the random phase approximation-for the translational frequencies of a neutron system trapped in a harmonic potential.

Ab Initio Optical Potentials and Nucleon Scattering on Medium Mass Nuclei

We derive ab initio optical potentials from self-consistent Green's function theory and compute the elastic scattering of neutrons off oxygen and calcium isotopes. The comparison with scattering data is satisfactory at low scattering energies. The method is benchmarked against the no-core shell model with continuum calculations, showing that virtual excitations of the target are crucial to predict proper fragmentation and absorption at higher energies. This is a significant step toward deriving optical potentials for medium mass nuclei and complex many-body systems in general.

Chiral Interactions up to Next-to-Next-to-Next-to-Leading Order and Nuclear Saturation

We present an efficient Monte Carlo framework for perturbative calculations of infinite nuclear matter based on chiral two-, three-, and four-nucleon interactions. The method enables the incorporation of all many-body contributions in a straightforward and transparent way, and makes it possible to extract systematic uncertainty estimates by performing order-by-order calculations in the chiral expansion as well as the many-body expansion. The versatility of this new framework is demonstrated by applying it to chiral low-momentum interactions, exhibiting a very good many-body convergence up to fourth order. Following these benchmarks, we explore new chiral interactions up to next-to-next-to-next-to-leading order (N^{3}LO). Remarkably, simultaneous fits to the triton and to saturation properties can be achieved, while all three-nucleon low-energy couplings remain natural. The theoretical uncertainties of nuclear matter are significantly reduced when going from next-to-next-to-leading order to N^{3}LO.

Tidal Deformabilities and Radii of Neutron Stars from the Observation of GW170817

We use gravitational-wave observations of the binary neutron star merger GW170817 to explore the tidal deformabilities and radii of neutron stars. We perform a Bayesian parameter estimation with the source location and distance informed by electromagnetic observations. We also assume that the two stars have the same equation of state; we demonstrate that, for stars with masses comparable to the component masses of GW170817, this is effectively implemented by assuming that the stars' dimensionless tidal deformabilities are determined by the binary's mass ratio q by Λ_{1}/Λ_{2}=q^{6}. We investigate different choices of prior on the component masses of the neutron stars. We find that the tidal deformability and 90% credible interval is Λ[over ˜]=222_{-138}^{+420} for a uniform component mass prior, Λ[over ˜]=245_{-151}^{+453} for a component mass prior informed by radio observations of Galactic double neutron stars, and Λ[over ˜]=233_{-144}^{+448} for a component mass prior informed by radio pulsars. We find a robust measurement of the common areal radius of the neutron stars across all mass priors of 8.9≤R[over ^]≤13.2  km, with a mean value of ⟨R[over ^]⟩=10.8  km. Our results are the first measurement of tidal deformability with a physical constraint on the star's equation of state and place the first lower bounds on the deformability and areal radii of neutron stars using gravitational waves.

Properties of Nuclei up to A=16 using Local Chiral Interactions

We report accurate quantum Monte Carlo calculations of nuclei up to A=16 based on local chiral two- and three-nucleon interactions up to next-to-next-to-leading order. We examine the theoretical uncertainties associated with the chiral expansion and the cutoff in the theory, as well as the associated operator choices in the three-nucleon interactions. While in light nuclei the cutoff variation and systematic uncertainties are rather small, in ^{16}O these can be significant for large coordinate-space cutoffs. Overall, we show that chiral interactions constructed to reproduce properties of very light systems and nucleon-nucleon scattering give an excellent description of binding energies, charge radii, and form factors for all these nuclei, including open-shell systems in A=6 and 12.

Light-Nuclei Spectra from Chiral Dynamics

In recent years local chiral interactions have been derived and implemented in quantum Monte Carlo methods in order to test to what extent the chiral effective field theory framework impacts our knowledge of few- and many-body systems. In this Letter, we present Green's function Monte Carlo calculations of light nuclei based on the family of local two-body interactions presented by our group in a previous paper in conjunction with chiral three-body interactions fitted to bound- and scattering-state observables in the three-nucleon sector. These interactions include Δ intermediate states in their two-pion-exchange components. We obtain predictions for the energy levels and level ordering of nuclei in the mass range A=4-12, accurate to ≤2% of the binding energy, in very satisfactory agreement with experimental data.

Short-Range Correlations and the EMC Effect in Effective Field Theory

We show that the empirical linear relation between the magnitude of the EMC effect in deep inelastic scattering on nuclei and the short-range correlation scaling factor a_{2} extracted from high-energy quasielastic scattering at x≥1 is a natural consequence of scale separation and derive the relationship using effective field theory. While the scaling factor a_{2} is a ratio of nuclear matrix elements that individually depend on the calculational scheme, we show that the ratio is independent of this choice. We perform Green's function Monte Carlo calculations with both chiral and Argonne-Urbana potentials to verify this and determine the scaling factors for light nuclei. The resulting values for ^{3}He and ^{4}He are in good agreement with experimental values. We also present results for ^{9}Be and ^{12}C extracted from variational Monte Carlo calculations.

Status and future of nuclear matrix elements for neutrinoless double-beta decay: a review

The nuclear matrix elements that govern the rate of neutrinoless double beta decay must be accurately calculated if experiments are to reach their full potential. Theorists have been working on the problem for a long time but have recently stepped up their efforts as ton-scale experiments have begun to look feasible. Here we review past and recent work on the matrix elements in a wide variety of nuclear models and discuss work that will be done in the near future. Ab initio nuclear-structure theory, which is developing rapidly, holds out hope of more accurate matrix elements with quantifiable error bars.