Charge form factor and sum rules of electromagnetic response functions in 12C

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.

Variational Monte Carlo Calculations of A≤4 Nuclei with an Artificial Neural-Network Correlator Ansatz

The complexity of many-body quantum wave functions is a central aspect of several fields of physics and chemistry where nonperturbative interactions are prominent. Artificial neural networks (ANNs) have proven to be a flexible tool to approximate quantum many-body states in condensed matter and chemistry problems. In this work we introduce a neural-network quantum state ansatz to model the ground-state wave function of light nuclei, and approximately solve the nuclear many-body Schrödinger equation. Using efficient stochastic sampling and optimization schemes, our approach extends pioneering applications of ANNs in the field, which present exponentially scaling algorithmic complexity. We compute the binding energies and point-nucleon densities of A≤4 nuclei as emerging from a leading-order pionless effective field theory Hamiltonian. We successfully benchmark the ANN wave function against more conventional parametrizations based on two- and three-body Jastrow functions, and virtually exact Green's function Monte Carlo results.

Nuclear Charge Radii of ^{10,11}B

The first laser spectroscopic determination of the change in the nuclear charge radius for a five-electron system is reported. This is achieved by combining high-accuracy ab initio mass-shift calculations and a high-accuracy measurement of the isotope shift in the 2s^{2}2p  ^{2}P_{1/2}→2s^{2}3s  ^{2}S_{1/2} ground state transition in boron atoms. Accuracy is increased by orders of magnitude for the stable isotopes ^{10,11}B and the results are used to extract their difference in the mean-square charge radius ⟨r_{c}^{2}⟩^{11}-⟨r_{c}^{2}⟩^{10}=-0.49(12)  fm^{2}. The result is qualitatively explained by a possible cluster structure of the boron nuclei and quantitatively used to benchmark new ab initio nuclear structure calculations using the no-core shell model and Green's function Monte Carlo approaches. These results are the foundation for a laser spectroscopic determination of the charge radius of the proton-halo candidate ^{8}B.

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.

Electromagnetic Response of ^{12}C: A First-Principles Calculation

The longitudinal and transverse electromagnetic response functions of ^{12}C are computed in a "first-principles" Green's function Monte Carlo calculation, based on realistic two- and three-nucleon interactions and associated one- and two-body currents. We find excellent agreement between theory and experiment and, in particular, no evidence for the quenching of the measured versus calculated longitudinal response. This is further corroborated by a reanalysis of the Coulomb sum rule, in which the contributions from the low-lying J^{π}=2^{+}, 0_{2}^{+} (Hoyle), and 4^{+} states in ^{12}C are accounted for explicitly in evaluating the total inelastic strength.

Unified Description of Electron-Nucleus Scattering within the Spectral Function Formalism

The formalism based on factorization and nuclear spectral functions has been generalized to treat transition matrix elements involving two-nucleon currents, whose contribution to the nuclear electromagnetic response in the transverse channel is known to be significant. We report the results of calculations of the inclusive electron-carbon cross section, showing that the inclusion of processes involving two-nucleon currents appreciably improves the agreement between theory and data in the dip region, between the quasielastic and Δ-production peaks. The relation to approaches based on the independent particle of the nucleus and the implications for the analysis of flux-integrated neutrino-nucleus cross sections are discussed.

Relativistic and Nuclear Medium Effects on the Coulomb Sum Rule

In light of the forthcoming high precision quasielastic electron scattering data from Jefferson Lab, it is timely for the various approaches to nuclear structure to make robust predictions for the associated response functions. With this in mind, we focus here on the longitudinal response function and the corresponding Coulomb sum rule for isospin-symmetric nuclear matter at various baryon densities. Using a quantum field-theoretic quark-level approach which preserves the symmetries of quantum chromodynamics, as well as exhibiting dynamical chiral symmetry breaking and quark confinement, we find a dramatic quenching of the Coulomb sum rule for momentum transfers |q|≳0.5  GeV. The main driver of this effect lies in changes to the proton Dirac form factor induced by the nuclear medium. Such a dramatic quenching of the Coulomb sum rule was not seen in a recent quantum Monte Carlo calculation for carbon, suggesting that the Jefferson Lab data may well shed new light on the explicit role of QCD in nuclei.

Quantum Monte Carlo calculations of light nuclei using chiral potentials

We present the first Green's function Monte Carlo calculations of light nuclei with nuclear interactions derived from chiral effective field theory up to next-to-next-to-leading order. Up to this order, the interactions can be constructed in a local form and are therefore amenable to quantum Monte Carlo calculations. We demonstrate a systematic improvement with each order for the binding energies of A=3 and A=4 systems. We also carry out the first few-body tests to study perturbative expansions of chiral potentials at different orders, finding that higher-order corrections are more perturbative for softer interactions. Our results confirm the necessity of a three-body force for correct reproduction of experimental binding energies and radii, and pave the way for studying few- and many-nucleon systems using quantum Monte Carlo methods with chiral interactions.

Neutral weak current two-body contributions in inclusive scattering from 12C

An ab initio calculation of the sum rules of the neutral weak response functions in 12C is reported, based on a realistic Hamiltonian, including two- and three-nucleon potentials, and on realistic currents, consisting of one- and two-body terms. We find that the sum rules of the response functions associated with the longitudinal and transverse components of the (spacelike) neutral current are largest and that a significant portion (≃30%) of the calculated strength is due to two-body terms. This fact may have implications for the MiniBooNE and other neutrino quasielastic scattering data on nuclei.