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.

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.

Nuclear Charge Radii of Silicon Isotopes.

The nuclear charge radius of ^{32}Si was determined using collinear laser spectroscopy. The experimental result was confronted with ab initio nuclear lattice effective field theory, valence-space in-medium similarity renormalization group, and mean field calculations, highlighting important achievements and challenges of modern many-body methods. The charge radius of ^{32}Si completes the radii of the mirror pair ^{32}Ar-^{32}Si, whose difference was correlated to the slope L of the symmetry energy in the nuclear equation of state. Our result suggests L≤60 MeV, which agrees with complementary observables.

Ab Initio Calculation of the Alpha-Particle Monopole Transition Form Factor.

We present a parameter-free ab initio calculation of the α-particle monopole transition form factor in the framework of nuclear lattice effective field theory. We use a minimal nuclear interaction that was previously used to reproduce the ground state properties of light nuclei, medium-mass nuclei, and neutron matter simultaneously with no more than a few percent error in the energies and charge radii. The results for the monopole transition form factor are in good agreement with recent precision data from Mainz.

New insights into the oscillations of the nucleon electromagnetic form factors.

The electromagnetic form factors of the proton and the neutron in the timelike region are investigated. Electron-positron annihilation into antinucleon-nucleon (N¯N) pairs is treated in distorted wave Born approximation, including the final-state interaction in the N¯N system. The latter is obtained by a Lippmann-Schwinger equation for N¯N potentials derived within SU(3) chiral effective field theory. By fitting to the phase shifts and (differential) cross section data, a high quality description is achieved. With these amplitudes, the oscillations of the electromagnetic form factors of the proton and the neutron are studied. It is found that each of them can be described by two fractional oscillators. One is characterized as "overdamped" and dominates near the threshold, while the other is "underdamped" and plays an important role in the high-energy region. These two oscillators are essential to understand the distributions of polarized electric charges induced by hard photons for the nucleons.

Emergent geometry and duality in the carbon nucleus.

The carbon atom provides the backbone for the complex organic chemistry composing the building blocks of life. The physics of the carbon nucleus in its predominant isotope, 12C, is similarly full of multifaceted complexity. Here we provide a model-independent density map of the geometry of the nuclear states of 12C using the ab initio framework of nuclear lattice effective field theory. We find that the well-known but enigmatic Hoyle state is composed of a "bent-arm" or obtuse triangular arrangement of alpha clusters. We identify all of the low-lying nuclear states of 12C as having an intrinsic shape composed of three alpha clusters forming either an equilateral triangle or an obtuse triangle. The states with the equilateral triangle formation also have a dual description in terms of particle-hole excitations in the mean-field picture.

Revealing the nature of hidden charm pentaquarks with machine learning.

We study the nature of the hidden charm pentaquarks, i.e., the Pc4312,Pc4440 and Pc(4457), with a neural network approach in pionless effective field theory. In this framework, the normal χ2 fitting approach cannot distinguish the quantum numbers of the Pc(4440) and Pc(4457). In contrast to that, the neural network-based approach can discriminate them, which still cannot be seen as a proof of the spin of the states since pion exchange is not considered in the approach. In addition, we also illustrate the role of each experimental data bin of the invariant J/ψp mass distribution on the underlying physics in both neural network and fitting methods. Their similarities and differences demonstrate that neural network methods can use data information more effectively and directly. This study provides more insights about how the neural network-based approach predicts the nature of exotic states from the mass spectrum.

Floating Block Method for Quantum Monte Carlo Simulations.

Quantum Monte Carlo simulations are powerful and versatile tools for the quantum many-body problem. In addition to the usual calculations of energies and eigenstate observables, quantum Monte Carlo simulations can in principle be used to build fast and accurate many-body emulators using eigenvector continuation or design time-dependent Hamiltonians for adiabatic quantum computing. These new applications require something that is missing from the published literature, an efficient quantum Monte Carlo scheme for computing the inner product of ground state eigenvectors corresponding to different Hamiltonians. In this work, we introduce an algorithm called the floating block method, which solves the problem by performing Euclidean time evolution with two different Hamiltonians and interleaving the corresponding time blocks. We use the floating block method and nuclear lattice simulations to build eigenvector continuation emulators for energies of ^{4}He, ^{8}Be, ^{12}C, and ^{16}O nuclei over a range of local and nonlocal interaction couplings. From the emulator data, we identify the quantum phase transition line from a Bose gas of alpha particles to a nuclear liquid.

Perturbative Quantum Monte Carlo Method for Nuclear Physics.

While first order perturbation theory is routinely used in quantum Monte Carlo (QMC) calculations, higher-order terms present significant numerical challenges. We present a new approach for computing perturbative corrections in projection QMC calculations. We demonstrate the method by computing nuclear ground state energies up to second order for a realistic chiral interaction. We calculate the binding energies of several light nuclei up to ^{16}O by expanding the Hamiltonian around the Wigner SU(4) limit and find good agreement with data. In contrast to the natural ordering of the perturbative series, we find remarkably large second-order energy corrections. This occurs because the perturbing interactions break the symmetries of the unperturbed Hamiltonian. Our method is free from the sign problem and can be applied to QMC calculations for many-body systems in nuclear physics, condensed matter physics, ultracold atoms, and quantum chemistry.

New Insights into the Nucleon's Electromagnetic Structure.

We present a combined analysis of the electromagnetic form factors of the nucleon in the space- and timelike regions using dispersion theory. Our framework provides a consistent description of the experimental data over the full range of momentum transfer, in line with the strictures from analyticity and unitarity. The statistical uncertainties of the extracted form factors are estimated using the bootstrap method, while systematic errors are determined from variations of the spectral functions. We also perform a high-precision extraction of the nucleon radii and find good agreement with previous analyses of spacelike data alone. For the proton charge radius, we find r_{E}^{p}=0.840_{-0.002}^{+0.003} _{-0.002}^{+0.002} fm, where the first error is statistical and the second one is systematic. The Zemach radius and third moment are in agreement with Lamb shift measurements and hyperfine splittings. The combined dataset of space- and timelike data disfavors a zero crossing of µ_{p}G_{E}^{p}/G_{M}^{p} in the spacelike region. Finally, we discuss the status and perspectives of modulus and phase of the form factors in the timelike region in the context of future experiments, as well as the onset of perturbative QCD.

Hidden Spin-Isospin Exchange Symmetry.

The strong interactions among nucleons have an approximate spin-isospin exchange symmetry that arises from the properties of quantum chromodynamics in the limit of many colors, N_{c}. However this large-N_{c} symmetry is well hidden and reveals itself only when averaging over intrinsic spin orientations. Furthermore, the symmetry is obscured unless the momentum resolution scale is close to an optimal scale that we call Λ_{large-N_{c}}. We show that the large-N_{c} derivation requires a momentum resolution scale of Λ_{large-N_{c}}∼500 MeV. We derive a set of spin-isospin exchange sum rules and discuss implications for the spectrum of ^{30}P and applications to nuclear forces, nuclear structure calculations, and three-nucleon interactions.

Where Is the Lightest Charmed Scalar Meson?

The lightest charmed scalar meson is known as the D_{0}^{*}(2300), which is one of the earliest new hadron resonances observed at modern B factories. We show here that the parameters assigned to the lightest scalar D meson are in conflict with the precise LHCb data of the decay B^{-}→D^{+}π^{-}π^{-}. On the contrary, these data can be well described by an unitarized chiral amplitude containing a much lighter charmed scalar meson, the D_{0}^{*}(2100). We also extract the low-energy S-wave Dπ phase of the decay B^{-}→D^{+}π^{-}π^{-} from the data in a model-independent way, and show that its difference from the Dπ scattering phase shift can be traced back to an intermediate ρ^{-} exchange. Our work highlights that an analysis of data consistent with chiral symmetry, unitarity, and analyticity is mandatory in order to extract the properties of the ground-state scalar mesons in the singly heavy sector correctly, in analogy to the light scalar mesons f_{0}(500) and K_{0}^{*}(700).

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.

Interpretation of the LHCb P_{c} States as Hadronic Molecules and Hints of a Narrow P_{c}(4380).

Three hidden-charm pentaquark P_{c} states, P_{c}(4312), P_{c}(4440), and P_{c}(4457) were revealed in the Λ_{b}^{0}→J/ψpK^{-} process measured by LHCb using both run I and run II data. Their nature is under lively discussion, and their quantum numbers have not been determined. We analyze the J/ψp invariant mass distributions under the assumption that the crossed-channel effects provide a smooth background. For the first time, such an analysis is performed employing a coupled-channel formalism with the scattering potential involving both one-pion exchange as well as short-range operators constrained by heavy quark spin symmetry. We find that the data can be well described in the hadronic molecular picture, which predicts seven Σ_{c}^{(*)}D[over ¯]^{(*)} molecular states in two spin multiplets, such that the P_{c}(4312) is mainly a Σ_{c}D[over ¯] bound state with J^{P}=1/2^{-}, while P_{c}(4440) and P_{c}(4457) are Σ_{c}D[over ¯]^{*} bound states with quantum numbers 3/2^{-} and 1/2^{-}, respectively. We also show that there is evidence for a narrow Σ_{c}^{*}D[over ¯] bound state in the data which we call P_{c}(4380), different from the broad one reported by LHCb in 2015. With this state included, all predicted Σ_{c}D[over ¯], Σ_{c}^{*}D[over ¯], and Σ_{c}D[over ¯]^{*} hadronic molecules are seen in the data, while the missing three Σ_{c}^{*}D[over ¯]^{*} states are expected to be found in future runs of the LHC or in photoproduction experiments.

Toward a First-Principles Calculation of Electroweak Box Diagrams.

We derive a Feynman-Hellmann theorem relating the second-order nucleon energy shift resulting from the introduction of periodic source terms of electromagnetic and isovector axial currents to the parity-odd nucleon structure function F_{3}^{N}. It is a crucial ingredient in the theoretical study of the γW and γZ box diagrams that are known to suffer from large hadronic uncertainties. We demonstrate that for a given Q^{2} one only needs to compute a small number of energy shifts in order to obtain the required inputs for the box diagrams. Future lattice calculations based on this approach may shed new light on various topics in precision physics including the refined determination of the Cabibbo-Kobayashi-Maskawa matrix elements and the weak mixing angle.

Ab initio Calculations of the Isotopic Dependence of Nuclear Clustering.

Nuclear clustering describes the appearance of structures resembling smaller nuclei such as alpha particles (^{4}He nuclei) within the interior of a larger nucleus. In this Letter, we present lattice Monte Carlo calculations based on chiral effective field theory for the ground states of helium, beryllium, carbon, and oxygen isotopes. By computing model-independent measures that probe three- and four-nucleon correlations at short distances, we determine the shape of the alpha clusters and the entanglement of nucleons comprising each alpha cluster with the outside medium. We also introduce a new computational approach called the pinhole algorithm, which solves a long-standing deficiency of auxiliary-field Monte Carlo simulations in computing density correlations relative to the center of mass. We use the pinhole algorithm to determine the proton and neutron density distributions and the geometry of cluster correlations in ^{12}C, ^{14}C, and ^{16}C. The structural similarities among the carbon isotopes suggest that ^{14}C and ^{16}C have excitations analogous to the well-known Hoyle state resonance in ^{12}C.

Nuclear Binding Near a Quantum Phase Transition.

How do protons and neutrons bind to form nuclei? This is the central question of ab initio nuclear structure theory. While the answer may seem as simple as the fact that nuclear forces are attractive, the full story is more complex and interesting. In this work we present numerical evidence from ab initio lattice simulations showing that nature is near a quantum phase transition, a zero-temperature transition driven by quantum fluctuations. Using lattice effective field theory, we perform Monte Carlo simulations for systems with up to twenty nucleons. For even and equal numbers of protons and neutrons, we discover a first-order transition at zero temperature from a Bose-condensed gas of alpha particles (^{4}He nuclei) to a nuclear liquid. Whether one has an alpha-particle gas or nuclear liquid is determined by the strength of the alpha-alpha interactions, and we show that the alpha-alpha interactions depend on the strength and locality of the nucleon-nucleon interactions. This insight should be useful in improving calculations of nuclear structure and important astrophysical reactions involving alpha capture on nuclei. Our findings also provide a tool to probe the structure of alpha cluster states such as the Hoyle state responsible for the production of carbon in red giant stars and point to a connection between nuclear states and the universal physics of bosons at large scattering length.

Erratum: Spectrum of Three-Body Bound States in a Finite Volume [Phys. Rev. Lett. 114, 091602 (2015)].

This corrects the article DOI: 10.1103/PhysRevLett.114.091602.