Search

Neutron matter is an important many-body system that provides valuable constraints for the equation of state (EOS) of neutron stars. Neutron-matter calculations employing chiral effective field theory (EFT) interactions have been extensively used for this purpose. Among the various many-body methods, quantum Monte Carlo (QMC) methods stand out due to their nonperturbative nature and the achievable precision. However, QMC methods require local interactions as input, which leads to the appearance of stronger regulator artifacts as compared to non-local interactions. To circumvent this, we employ large-cutoff interactions derived within chiral EFT (400 MeV $\leq \Lambda_c \leq$ 700 MeV) for studies of pure neutron matter. These interactions have been adjusted to nucleon-nucleon scattering phase shifts, the triton binding energy, as well as the triton beta-decay half life. We find that regulator artifacts significantly decrease with increasing cutoff, leading to a significant reduction of uncertainties in the neutron-matter EOS. We discuss implications for the symmetry energy and demonstrate how our new calculations lead to a reduction in the theoretical uncertainty of predicted neutron-star radii by up to 30% for low-mass stars., Comment: 7 pages, 3 figures

We calculate the magnetic moments of light nuclei ($A < 20$) using the auxiliary field diffusion Monte Carlo method and local two- and three-nucleon forces with electromagnetic currents from chiral effective field theory. For all nuclei under consideration, we also calculate the ground-state energies and charge radii. We generally find a good agreement with experimental values for all of these observables. For the electromagnetic currents, we explore the impact of employing two different power countings, and study theoretical uncertainties stemming from the truncation of the chiral expansion order-by-order for select nuclei within these two approaches. We find that it is crucial to employ consistent power countings for interactions and currents to achieve a systematic order-by-order convergence., Comment: 7 pages, 4 figures, supplemental material

This white paper was submitted to the 2022 Fundamental Symmetries, Neutrons, and Neutrinos (FSNN) Town Hall Meeting in preparation for the next NSAC Long Range Plan. We advocate to support current and future theoretical and experimental searches for physics beyond the Standard Model using nuclear $\beta$ decay., Comment: 12 pages, 2 figures, White paper submitted to the 2022 Fundamental Symmetries, Neutrons, and Neutrinos (FSNN) Town Hall Meeting at the University of North Carolina Chapel Hill

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., Comment: 14 pages, 7 figures; v2 corresponds to published version

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., Comment: 10 pages, 6 figures; Version accepted for publication, includes Supplemental Material

Thirty years ago, Steven Weinberg published his seminal paper on "Nuclear Forces from chiral Lagrangians" which has revolutionized the field of theoretical nuclear physics. Nowadays, interactions derived from chiral effective field theory are routinely used to describe nuclear systems ranging from atomic nuclei to the dense matter explored in the core of neutron stars with theoretical uncertainty estimates. In our contribution to the special issue "Celebrating 30 years of Steven Weinberg's paper Nuclear Forces from Chiral Lagrangians", we focus on the impact that chiral effective field theory interactions have played in advancing microscopic studies of atomic nuclei and the nuclear-matter equation of state using quantum Monte Carlo methods., Comment: 9 pages, 3 figures, invited contribution to the special issue of Few-Body Systems "Celebrating 30 years of Steven Weinberg's paper on nuclear forces from chiral Lagrangians"

We briefly review the concept of scaling and how it occurs in quasielastic electron and neutrino scattering from nuclei, and then the particular approach to scaling in the short-time approximation. We show that, while two-nucleon currents do significantly enhance the transverse electromagnetic response, they do not spoil scaling, but in fact enhance it. We provide scaling results obtained in the short-time approximation that verify this claim. The enhanced scaling is not ``accidental'' -- as claimed in Ref.~\cite{Benhar:2020jye} -- but rather reflects the dominant role played by pion exchange interactions and currents in the quasielastic regime., Comment: 2 figures, 4 pages

We report Quantum Monte Carlo calculations of weak transitions in $A\leq 10$ nuclei, based on the Norfolk two- and three-nucleon chiral interactions, and associated one- and two-body axial currents. We find that the contribution from two-body currents is at the $2$ - $3\%$ level, with the exception of matrix elements entering the rates of $^8$Li, $^8$B, and $^8$He beta decays. These matrix elements are suppressed in impulse approximation based on the (leading order) Gamow Teller transition operator alone; two-body currents provide a $20$ - $30\%$ correction, which is, however, insufficient to bring theory in agreement with experimental data. For the other transitions, the agreement with the data is satisfactory, and the results exhibit a negligible to mild model dependence when different combinations of Norfolk interactions are utilized to construct the nuclear wave functions. We report a complete study of two-body weak transition densities which reveals the expected universal behavior of two-body currents at short distances throughout the range of $A\,$=$\,3$ to $A\,$=$\,10$ systems considered here., Comment: 14 pages, 9 figures; Corrected experimental values for $^8$Li beta decay in Table V to be consistent with Table VI and the text; Corrected an error in row 2 column 4 of Figure 6, conclusions unchanged; Corrected errors in column 3 rows 3 and 4 of Figure 7, results unchanged; Corrected a typo in Equation 19, results in tables unchanged; Corrected $^8$He ground state isospin from $T=1$ to $T=2$

We carry out an ab initio calculation of the neutrino flux-folded inclusive cross sections, measured on $^{12}$C by the MiniBooNE and T2K collaborations in the charged-current quasielastic (CCQE) regime. The calculation is based on realistic two- and three-nucleon interactions, and on a realistic nuclear electroweak current with one-and two-nucleon terms that are constructed consistently with these interactions and reproduce low-energy electroweak transitions. Numerically exact quantum Monte Carlo methods are utilized to compute the nuclear weak response functions, by fully retaining many-body correlations in the initial and final states and interference effects between one- and two-body current contributions. We employ a nucleon axial form factor of the dipole form with $\Lambda_A = 1.0$ or $1.15$ GeV, the latter more in line with a very recent lattice QCD determination. The calculated cross sections are found to be in good agreement with the neutrino data of MiniBooNE and T2K, and antineutrino MiniBooNE data, yielding a consistent picture of nuclei and their electroweak properties across a wide regime of energy and momenta., Comment: 13 pages, 7 figures

We report a quantum Monte Carlo calculation of the equation of state of symmetric nuclear matter using local interactions derived from chiral effective field theory up to next-to-next-to-leading order fit to few-body observables only. The empirical saturation density and energy are well reproduced within statistical and systematic uncertainties. We have also derived the symmetry energy as a function of the density, finding good agreement with available experimentally derived constraints at saturation and twice saturation density. We find that the corresponding pressure is also in excellent agreement with recent constraints extracted from gravitational waves of the neutron-star merger GW170817 by the LIGO-Virgo detection., Comment: 7 pages, 3 figures, 3 tables

We reply to a Comment on our Letter [Phys. Rev. Lett. 118, 232501 (2017), arXiv:1612.01502] by A. Deltuva and R. Lazauskas [Phys. Rev. Lett 123, 069201 (2019), arXiv:1904.00925]., Comment: 2 pages, published version

High-energy scattering processes, such as deep inelastic scattering (DIS) and quasielastic (QE) scattering provide a wealth of information about the structure of atomic nuclei. The remarkable discovery of the empirical linear relationship between the slope of the European Muon Collaboration (EMC) effect in DIS and the short-range-correlation (SRC) scaling factors $a_2$ in QE kinematics is naturally explained in terms of scale separation in effective field theory. This explanation has powerful consequences, allowing us to calculate and predict SRC scaling factors from ab initio low-energy nuclear theory. We present ab initio calculations of SRC scaling factors for a nucleus $A$ relative to the deuteron $a_2(A/d)$ and relative to $^3\rm He$ $a_2(A/^3\rm He)$ in light and medium-mass nuclei. Our framework further predicts that the EMC effect and SRC scaling factors have minimal or negligible isovector corrections., Comment: 25 pages, 6 figures, 2 tables

Recent observations of neutron-star properties, in particular the recent detection of gravitational waves emitted from binary neutron stars, GW 170817, open the way to put strong constraints on nuclear interactions. In this paper, we review the state of the art in calculating the equation of state of strongly interacting matter from first principle calculations starting from microscopic interactions among nucleons. We then review selected properties of neutron stars that can be directly compared with present and future observations., Comment: 28 pages, 9 figures, Special Issue Article for Journal of Physics G

In recent years, the combination of precise quantum Monte Carlo (QMC) methods with realistic nuclear interactions and consistent electroweak currents, in particular those constructed within effective field theories (EFTs), has lead to new insights in light and medium-mass nuclei, neutron matter, and electroweak reactions. This compelling new body of work has been made possible both by advances in QMC methods for nuclear physics, which push the bounds of applicability to heavier nuclei and to asymmetric nuclear matter and by the development of local chiral EFT interactions up to next-to-next-to-leading order and minimally nonlocal interactions including $\Delta$ degrees of freedom. In this review, we discuss these recent developments and give an overview of the exciting results for nuclei, neutron matter and neutron stars, and electroweak reactions., Comment: 27 pages, 8 figures; Invited review prepared for the Annual Review of Nuclear and Particle Science, Volume 69 (2019)

We present quantum Monte Carlo calculations of the single- and two-nucleon momentum distributions in selected nuclei for $A\le16$. We employ local chiral interactions at next-to-next-to-leading order. We find good agreement at low momentum with the single-nucleon momentum distributions derived for phenomenological potentials. The same agreement is found for the integrated two-nucleon momentum distributions at low relative momentum $q$ and low center-of-mass momentum $Q$. We provide results for the two-nucleon momentum distributions as a function of both $q$ and $Q$. The large ratio of $pn$ to $pp$ pairs around $q=2\,\rm fm^{-1}$ for back-to-back $(Q=0)$ pairs is confirmed up to $^{16}$O, and results are compatible with those extracted from available experimental data., Comment: 10 pages, 14 figures

Quantum Monte Carlo methods have recently been employed to study properties of nuclei and infinite matter using local chiral effective field theory interactions. In this work, we present a detailed description of the auxiliary field diffusion Monte Carlo algorithm for nuclei in combination with local chiral two- and three-nucleon interactions up to next-to-next-to-leading order. We show results for the binding energy, charge radius, charge form factor, and Coulomb sum rule in nuclei with $3\le A\le16$. Particular attention is devoted to the effect of different operator structures in the three-body force for different cutoffs. The outcomes suggest that local chiral interactions fit to few-body observables give a very good description of the ground-state properties of nuclei up to $^{16}$O, with the exception of one fit for the softer cutoff which predicts overbinding in larger nuclei., Comment: 23 pages, 10 figures

Quasielastic neutrino scattering is an important aspect of the experimental program to study fundamental neutrino properties including neutrino masses, mixing angles, the mass hierarchy and CP-violating phase. Proper interpretation of the experiments requires reliable theoretical calculations of neutrino-nucleus scattering. In this paper we present calculations of response functions and cross sections by neutral-current scattering of neutrinos off $^{12}$C. These calculations are based on realistic treatments of nuclear interactions and currents, the latter including the axial-, vector-, and vector-axial interference terms crucial for determining the difference between neutrino and anti-neutrino scattering and the CP-violating phase. We find that the strength and energy-dependence of two-nucleon processes induced by correlation effects and interaction currents are crucial in providing the most accurate description of neutrino-nucleus scattering in the quasielastic regime., Comment: 5 pages, 2 figures

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., Comment: 6 pages, 4 figures

Ab initio calculations of the Gamow-Teller (GT) matrix elements in the $\beta$ decays of $^6$He and $^{10}$C and electron captures in $^7$Be are carried out using both variational and Green's function Monte Carlo wave functions obtained from the Argonne $v_{18}$ two-nucleon and Illinois-7 three-nucleon interactions, and axial many-body currents derived from either meson-exchange phenomenology or chiral effective field theory. The agreement with experimental data is excellent for the electron captures in $^7$Be, while theory overestimates the $^6$He and $^{10}$C data by $\sim 2\%$ and $\sim 10\%$, respectively. We show that for these systems correlations in the nuclear wave functions are crucial to explain the data, while many-body currents increase by $\sim 2$--$3\%$ the one-body GT contributions. These findings suggest that the longstanding $g_A$-problem, i.e., the systematic overprediction ($\sim 20 \%$ in $A\le 18$ nuclei) of GT matrix elements in shell-model calculations, may be resolved, at least partially, by correlation effects., Comment: 6 pages, 2 figures

The properties of cold Bose gases at unitarity have been extensively investigated in the last few years both theoretically and experimentally. In this paper we use a family of interactions tuned to two-body unitarity and very weak three-body binding to demonstrate the universal properties of both clusters and matter. We determine the universal properties of finite clusters up to 60 particles and, for the first time, explicitly demonstrate the saturation of energy and density with particle number and compare with bulk properties. At saturation in the bulk we determine the energy, density, two- and three-body contacts and the condensate fraction. We find that uniform matter is more bound than three-body clusters by nearly two orders of magnitude, the two-body contact is very large in absolute terms, and yet the condensate fraction is also very large, greater than 90%. Equilibrium properties of these systems may be experimentally accessible through rapid quenching of weakly-interacting boson superfluids., Comment: 5 pages, 4 figures

Local chiral effective field theory interactions have recently been developed and used in the context of quantum Monte Carlo few- and many-body methods for nuclear physics. In this work, we go over detailed features of local chiral nucleon-nucleon interactions and examine their effect on properties of the deuteron, paying special attention to the perturbativeness of the expansion. We then turn to three-nucleon interactions, focusing on operator ambiguities and their interplay with regulator effects. We then discuss the nuclear Green's function Monte Carlo method, going over both wave-function correlations and approximations for the two- and three-body propagators. Following this, we present a range of results on light nuclei: Binding energies and distribution functions are contrasted and compared, starting from several different microscopic interactions., Comment: 21 pages, 14 figures, published version, Editor's Suggestion

We present quantum Monte Carlo calculations of few-neutron systems confined in external potentials based on local chiral interactions at next-to-next-to-leading order in chiral effective field theory. The energy and radial densities for these systems are calculated in different external Woods-Saxon potentials. We assume that their extrapolation to zero external-potential depth provides a quantitative estimate of three- and four-neutron resonances. The validity of this assumption is demonstrated by benchmarking with an exact diagonalization in the two-body case. We find that the extrapolated trineutron resonance, as well as the energy for shallow well depths, is lower than the tetraneutron resonance energy. This suggests that a three-neutron resonance exists below a four-neutron resonance in nature and is potentially measurable. To confirm that the relative ordering of three- and four-neutron resonances is not an artifact of the external confinement, we test that the odd-even staggering in the helium isotopic chain is reproduced within this approach. Finally, we discuss similarities between our results and ultracold Fermi gases., Comment: 6 pages, 5 figures, version compatible with published letter

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 measured versus calculated longitudinal response. This is further corroborated by a re-analysis of the Coulomb sum rule, in which the contributions from the low-lying $J^\pi\,$=$\, 2^+$, $0^+_2$ (Hoyle), and $4^+$ states in $^{12}$C are accounted for explicitly in evaluating the total inelastic strength., Comment: 5 pages, 3 figures

Multi-neutron systems are crucial to understanding the physics of neutron-rich nuclei and neutron stars. Neutron drops, neutrons confined in an external field, are investigated systematically in both non-relativistic and relativistic density functional theories and with ab initio calculations. We demonstrate a strong linear correlation, which is universal in the realm of mean-field models, between the rms radii of neutron drops and the neutron skin thickness of Pb-208 and Ca-48; i.e., the difference between the neutron and proton rms radii of a nucleus. Due to its high quality, this correlation can be used to deduce the radii of neutron drops from the measured neutron skin thickness in a model-independent way, and the radii obtained for neutron drops can provide a useful constraint for realistic three neutron forces. This correlation, together with high- precision measurements of the neutron skin thicknesses of Pb-208 and Ca-48, will have an enduring impact on the understanding of multi-neutron interactions, neutron-rich nuclei, neutron stars, etc., Comment: 6 pages, 4 figures

Ab initio calculations provide direct access to the properties of pure neutron systems that are challenging to study experimentally. In addition to their importance for fundamental physics, their properties are required as input for effective field theories of the strong interaction. In this work, we perform auxiliary-field diffusion Monte Carlo calculations of the ground and first excited state of two neutrons in a finite box, considering a simple contact potential as well as chiral effective field theory interactions. We compare the results against exact diagonalizations and present a detailed analysis of the finite-volume effects, whose understanding is crucial for determining observables from the calculated energies. Using the L\"uscher formula, we extract the low-energy S-wave scattering parameters from ground- and excited-state energies for different box sizes., Comment: 11 pages, 8 figures

We present quantum Monte Carlo calculations of light nuclei, neutron-$\alpha$ scattering, and neutron matter using local two- and three-nucleon ($3N$) interactions derived from chiral effective field theory up to next-to-next-to-leading order (N$^2$LO). The two undetermined $3N$ low-energy couplings are fit to the $^4$He binding energy and, for the first time, to the spin-orbit splitting in the neutron-$\alpha$ $P$-wave phase shifts. Furthermore, we investigate different choices of local $3N$-operator structures and find that chiral interactions at N$^2$LO are able to simultaneously reproduce the properties of $A=3,4,5$ systems and of neutron matter, in contrast to commonly used phenomenological $3N$ interactions., Comment: 5 pages, 3 figures, 1 table - updated version: small wording changes, one reference change

Precise measurement of neutrino oscillations, and hence the determination of their masses demands a quantitative understanding of neutrino-nucleus interactions. To this aim, two-body meson-exchange currents have to be accounted for along within realistic models of nuclear dynamics. We summarize our progresses towards the construction of a consistent framework, based on quantum Monte Carlo methods and on the spectral function approach, that can be exploited to accurately describe neutrino interactions with atomic nuclei over the broad kinematical region covered by neutrino experiments., Comment: 8 pages, 4 figure, Proceedings of the 21st International Conference on Few-Body Problems in Physics, Chicago, Illinois, USA

Chiral effective field theory (EFT) enables a systematic description of low-energy hadronic interactions with controlled theoretical uncertainties. For strongly interacting systems, quantum Monte Carlo (QMC) methods provide some of the most accurate solutions, but they require as input local potentials. We have recently constructed local chiral nucleon-nucleon (NN) interactions up to next-to-next-to-leading order (N$^2$LO). Chiral EFT naturally predicts consistent many-body forces. In this paper, we consider the leading chiral three-nucleon (3N) interactions in local form. These are included in auxiliary field diffusion Monte Carlo (AFDMC) simulations. We present results for the equation of state of neutron matter and for the energies and radii of neutron drops. In particular, we study the regulator dependence at the Hartree-Fock level and in AFDMC and find that present local regulators lead to less repulsion from 3N forces compared to the usual nonlocal regulators., Comment: 10 pages, 8 figures, 1 table, published version