Light-Nuclei Spectra from Chiral Dynamics

Benchmarking Λ NN three-body forces and first predictions for A = 3 - 5 hypernuclei

Explicit expressions for the leading chiral hyperon-nucleon-nucleon three-body forces have been derived by Petschauer et al (Phys Rev C93:014001, 2016). An important prerequisite for including these three-body forces in few- and many-body calculations is the accuracy and efficiency of their partial-wave decomposition. A careful benchmark of the Λ NN potential matrix elements, computed using two robust and efficient partial-wave decomposition methods, is presented. In addition, results of a first quantitative assessment for the contributions of Λ NN forces to the separation energies in A = 3 - 5 hypernuclei are reported.

Quantum Monte Carlo Calculations of Magnetic Form Factors in Light Nuclei

We present quantum Monte Carlo calculations of magnetic form factors in A=6-10 nuclei, based on Norfolk two- and three-nucleon interactions, and associated one- and two-body electromagnetic currents. Agreement with the available experimental data for ^{6}Li, ^{7}Li, ^{9}Be, and ^{10}B up to values of momentum transfer q∼3  fm^{-1} is achieved when two-nucleon currents are accounted for. We present a set of predictions for the magnetic form factors of ^{7}Be, ^{8}Li, ^{9}Li, and ^{9}C. In these systems, two-body currents account for ∼40%-60% of the total magnetic strength. Measurements in any of these radioactive systems would provide valuable insights on the nuclear magnetic structure emerging from the underlying many-nucleon dynamics. A particularly interesting case is that of ^{7}Be, as it would enable investigations of the magnetic structure of mirror nuclei.

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.

Wavefunction matching for solving quantum many-body problems

Ab initio calculations have an essential role in our fundamental understanding of quantum many-body systems across many subfields, from strongly correlated fermions1-3 to quantum chemistry4-6 and from atomic and molecular systems7-9 to nuclear physics10-14. One of the primary challenges is to perform accurate calculations for systems where the interactions may be complicated and difficult for the chosen computational method to handle. Here we address the problem by introducing an approach called wavefunction matching. Wavefunction matching transforms the interaction between particles so that the wavefunctions up to some finite range match that of an easily computable interaction. This allows for calculations of systems that would otherwise be impossible owing to problems such as Monte Carlo sign cancellations. We apply the method to lattice Monte Carlo simulations15,16 of light nuclei, medium-mass nuclei, neutron matter and nuclear matter. We use high-fidelity chiral effective field theory interactions17,18 and find good agreement with empirical data. These results are accompanied by insights on the nuclear interactions that may help to resolve long-standing challenges in accurately reproducing nuclear binding energies, charge radii and nuclear-matter saturation in ab initio calculations19,20.

Convergence of Eigenvector Continuation

Eigenvector continuation is a computational method that finds the extremal eigenvalues and eigenvectors of a Hamiltonian matrix with one or more control parameters. It does this by projection onto a subspace of eigenvectors corresponding to selected training values of the control parameters. The method has proven to be very efficient and accurate for interpolating and extrapolating eigenvectors. However, almost nothing is known about how the method converges, and its rapid convergence properties have remained mysterious. In this Letter, we present the first study of the convergence of eigenvector continuation. In order to perform the mathematical analysis, we introduce a new variant of eigenvector continuation that we call vector continuation. We first prove that eigenvector continuation and vector continuation have identical convergence properties and then analyze the convergence of vector continuation. Our analysis shows that, in general, eigenvector continuation converges more rapidly than perturbation theory. The faster convergence is achieved by eliminating a phenomenon that we call differential folding, the interference between nonorthogonal vectors appearing at different orders in perturbation theory. From our analysis we can predict how eigenvector continuation converges both inside and outside the radius of convergence of perturbation theory. While eigenvector continuation is a nonperturbative method, we show that its rate of convergence can be deduced from power series expansions of the eigenvectors. Our results also yield new insights into the nature of divergences in perturbation theory.

Energy-Dependent π^{+}π^{+}π^{+} Scattering Amplitude from QCD

Focusing on three-pion states with maximal isospin (π^{+}π^{+}π^{+}), we present the first nonperturbative determination of an energy-dependent three-hadron scattering amplitude from first-principles QCD. The calculation combines finite-volume three-hadron energies, extracted using numerical lattice QCD, with a relativistic finite-volume formalism, required to interpret the results. To fully implement the latter, we also solve integral equations that relate an intermediate three-body K matrix to the physical three-hadron scattering amplitude. The resulting amplitude shows rich analytic structure and a complicated dependence on the two-pion invariant masses, represented here via Dalitz-like plots of the scattering rate.

Ab Initio Nuclear Thermodynamics

We propose a new Monte Carlo method called the pinhole trace algorithm for ab initio calculations of the thermodynamics of nuclear systems. For typical simulations of interest, the computational speedup relative to conventional grand-canonical ensemble calculations can be as large as a factor of one thousand. Using a leading-order effective interaction that reproduces the properties of many atomic nuclei and neutron matter to a few percent accuracy, we determine the location of the critical point and the liquid-vapor coexistence line for symmetric nuclear matter with equal numbers of protons and neutrons. We also present the first ab initio study of the density and temperature dependence of nuclear clustering.

Clustering of Four-Component Unitary Fermions

Ab initio nuclear physics tackles the problem of strongly interacting four-component fermions. The same setting could foreseeably be probed experimentally in ultracold atomic systems, where two- and three-component experiments have led to major breakthroughs in recent years. Both due to the problem's inherent interest and as a pathway to nuclear physics, in this Letter we study four-component fermions at unitarity via the use of quantum Monte Carlo methods. We explore novel forms of the trial wave function and find one which leads to a ground state of the eight-particle system whose energy is almost equal to that of two four-particle systems. We investigate the clustering properties involved and also extrapolate to the zero-range limit. In addition to being experimentally testable, our results impact the prospects of developing nuclear physics as a perturbation around the unitary limit.

Nonperturbative Extraction of the Effective Mass in Neutron Matter

We carry out nonperturbative calculations of the single-particle excitation spectrum in strongly interacting neutron matter. These are microscopic quantum Monte Carlo computations of many-neutron energies at different densities as well as several distinct excited states. As input, we employ both phenomenological and chiral two- and three-nucleon interactions. We use the single-particle spectrum to extract the effective mass in neutron matter. With a view to systematizing the error involved in this extraction, we carefully assess the impact of finite-size effects on the quasiparticle dispersion relation. We find an effective-mass ratio that drops from 1 as the density is increased. We conclude by connecting our results with the physics of ultracold gases as well as with energy-density functional theories of nuclei and neutron-star matter.

Structure of the Lightest Tin Isotopes

We link the structure of nuclei around ^{100}Sn, the heaviest doubly magic nucleus with equal neutron and proton numbers (N=Z=50), to nucleon-nucleon (NN) and three-nucleon (NNN) forces constrained by data of few-nucleon systems. Our results indicate that ^{100}Sn is doubly magic, and we predict its quadrupole collectivity. We present precise computations of ^{101}Sn based on three-particle-two-hole excitations of ^{100}Sn, and we find that one interaction accurately reproduces the small splitting between the lowest J^{π}=7/2^{+} and 5/2^{+} states.

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.