BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
R. B. Wiringa appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
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 $β$ decay.
We present the equation of state of infinite neutron matter as obtained from highly-realistic Hamiltonians that include nucleon-nucleon and three-nucleon coordinate-space potentials. We benchmark three independent many-body methods: Brueckner-Bethe-Goldstone (BBG), Fermi hypernetted chain/single-operator chain (FHNC/SOC), and auxiliary-field diffusion Monte Carlo (AFDMC). We find them to provide similar equations of state when the Argonne $v_{18}$ and the Argonne $v_{6}^\prime$ nucleon-nucleon potentials are used in combination with the Urbana IX three-body force. Only at densities larger than about 1.5 the nuclear saturation density ($ρ_0 = 0.16\,\rm{fm}^{-3}$) the FHNC/SOC energies are appreciably lower than the other two approaches. The AFDMC calculations carried out with all of the Norfolk potentials fitted to reproduce the experimental trinucleon ground-state energies and $nd$ doublet scattering length yield unphysically bound neutron matter, associated with the formation of neutron droplets. Including tritium $β$-decay in the fitting procedure, as in the second family of Norfolk potentials, mitigates but does not completely resolve this problem. An excellent agreement between the BBG and AFDMC results is found for the subset of Norfolk interactions that do not make neutron-matter collapse, while the FHNC/SOC equations of state are moderately softer.
We perform first-principle calculations of electron-nucleus scattering on $^3$He and $^3$H using the Green's function Monte Carlo method and two approaches based on the factorization of the final hadronic state: the spectral-function formalism and the short-time approximation. These three methods are benchmarked among each other and compared to the experimental data for the longitudinal and transverse electromagnetic response functions of $^3$He, and the inclusive cross sections of both $^3$He and $^3$H. Since these three approaches are based on the same description of nuclear dynamics of the initial target state, comparing their results enables a precise quantification of the uncertainties inherent to factorization schemes. At sufficiently large values of the momentum transfer, we find an excellent agreement of the Green's function Monte Carlo calculation with experimental data and with both the spectral-function formalism and the short-time approximation. We also analyze the relevance of relativistic effects, whose inclusion becomes crucial to explain data at high momentum and energy transfer.
Ab-initio Quantum Monte Carlo (QMC) calculations of nuclei from deuterium to 40Ca, obtained using four different phenomenological and local chiral nuclear potentials, are analyzed using the Generalized Contact Formalism (GCF). We extract spin- and isospin-dependent "nuclear contact terms" for each interaction in both coordinate and momentum space. The extracted contact terms, that count the number of short-range correlated (SRC) pairs with different quantum numbers, are dependent on the nuclear interaction model used in the QMC calculation. However, the ratios of contact terms for a nucleus A to deuterium (for spin-1 pn pairs) or to 4He (for all NN pairs) are independent of the nuclear interaction model and are the same for both short-distance and high-momentum pairs. This implies that the relative abundance of short-range pairs in the nucleus is a long-range (mean-field) quantity that is insensitive to the short-distance nature of the nuclear force. Measurements of exclusive (e,e'NN) pair breakup processes are instead more sensitive to short-range dynamics
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.
We report benchmark calculations of the energy per particle of pure neutron matter as a function of the baryon density using three independent many-body methods: Brueckner-Bethe-Goldstone, Fermi hypernetted chain/single-operator chain, and auxiliary-field diffusion Monte Carlo. Significant technical improvements are implemented in the latter two methods. The calculations are made for two distinct families of realistic coordinate-space nucleon-nucleon potentials fit to scattering data, including the standard Argonne $v_{18}$ interaction and two of its simplified versions, and four of the new Norfolk $Δ$-full chiral effective field theory potentials. The results up to twice nuclear matter saturation density show some divergence among the methods, but improved agreement compared to earlier work. We find that the potentials fit to higher-energy nucleon-nucleon scattering data exhibit a much smaller spread of energies.
We study knockout reactions with proton probes within a theoretical framework where {\it ab initio} Quantum Monte Carlo wave functions are combined with the Faddeev/Alt-Grassberger-Sandhas few-body reaction formalism. New Quantum Monte Carlo wave functions are used to describe $^{12}$C, yielding, for the first time, results consistent with the experimental point rms radii, electron scattering data and (p,2p) total cross sections data. Our results for $\mathrm{A}\leq 12$ and $(N-Z) \leq 3$ nuclei show that the theoretical ratios between the (i) {\it ab initio} and Mean Field Approximation theoretical cross sections, $\mathcal{R}_σ$, (ii) corresponding ratios between the spectroscopic factors, $\mathcal{R}_Σ$, summed over states below particle emission, depend moderately on the nucleon separation energy S$_{\rm N}$. These ratios are determined by a delicate interplay between the radii of the parent and the residual nuclei and the nucleon separation energy, and were found to be always smaller for the knockout of the more correlated deficient species nucleon. In the case of the symmetric $^{12}$C nucleus, the theoretical ratios still appear to indicate that protons are more correlated than neutrons.
Improvements in theoretical modeling of Short Range structures and phenomena, and comparisons with data, will require sustained collaboration between nuclear theorists and neutrino experimentalists. The extensive history of studying this area of nuclear physics in electron- and hadron-scattering experiments, coupled with the transformative capabilities of LArTPCs to identify neutrinos, will provide a ripe opportunity for new discoveries that will further our understanding of the nucleus.
Quantum Monte Carlo methods have proved very valuable to study the structure and reactions of light nuclei and nucleonic matter starting from realistic nuclear interactions and currents. These ab-initio calculations reproduce many low-lying states, moments and transitions in light nuclei, and simultaneously predict many properties of light nuclei and neutron matter over a rather wide range of energy and momenta. We review the nuclear interactions and currents, and describe the continuum Quantum Monte Carlo methods used in nuclear physics. These methods are similar to those used in condensed matter and electronic structure but naturally include spin-isospin, tensor, spin-orbit, and three-body interactions. We present a variety of results including the low-lying spectra of light nuclei, nuclear form factors, and transition matrix elements. We also describe low-energy scattering techniques, studies of the electroweak response of nuclei relevant in electron and neutrino scattering, and the properties of dense nucleonic matter as found in neutron stars. A coherent picture of nuclear structure and dynamics emerges based upon rather simple but realistic interactions and currents.
We report quantum Monte Carlo calculations of electromagnetic transitions in 8Be. The realistic Argonne v18 two-nucleon and Illinois-7 three-nucleon potentials are used to generate the ground state and nine excited states, with energies that are in excellent agreement with experiment. A dozen M1 and eight E2 transition matrix elements between these states are then evaluated. The E2 matrix elements are computed only in impulse approximation, with those transitions from broad resonant states requiring special treatment. The M1 matrix elements include two-body meson-exchange currents derived from chiral effective field theory, which typically contribute 20--30% of the total expectation value. Many of the transitions are between isospin-mixed states; the calculations are performed for isospin-pure states and then combined with empirical mixing coefficients to compare to experiment. Alternate mixings are also explored. In general, we find that transitions between states that have the same dominant spatial symmetry are in reasonable agreement with experiment, but those transitions between different spatial symmetries are often underpredicted.
We report Green's function Monte Carlo calculations of isospin-mixing (IM) matrix elements for the 2+, 1+, and 3+ T=0,1 pairs of states at 16--19 MeV excitation in 8Be. The realistic Argonne v18 (AV18) two-nucleon and Illinois-7 three-nucleon potentials are used to generate the nuclear wave functions. Contributions from the full electromagnetic interaction and strong class III charge-symmetry-breaking (CSB) components of the AV18 potential are evaluated. We also examine two theoretically more complete CSB potentials based on rho-omega mixing, tuned to give the same neutron-neutron scattering length as AV18. The contribution of these different CSB potentials to the 3H-3He, 7Li-7Be, and 8Li-8B isovector energy differences is evaluated and reasonable agreement with experiment is obtained. Finally, for the 8Be IM calculation we add the small class IV CSB terms coming from one-photon, one-pion, and one-rho exchange, as well as rho-omega mixing. The expectation values of the three CSB models vary by up to 20% in the isovector energy differences, but only by 10% or less in the IM matrix element. The total matrix element gives 85--90% of the experimental IM value of -145 keV for the 2+ doublet, with about two thirds coming from the Coulomb interaction. We also report the IM matrix element to the first 2+ state at 3 MeV excitation, which is the final state for various tests of the Standard Model for beta-decay.
An earlier measurement on the 4$^+$ to 2$^+$ radiative transition in $^8$Be provided the first electromagnetic signature of its dumbbell-like shape. However, the large uncertainty in the measured cross section does not allow a stringent test of nuclear structure models. The present paper reports a more elaborate and precise measurement for this transition, via the radiative capture in the $^4$He+$^4$He reaction, improving the accuracy by about a factor of three. The {\it ab initio} calculations of the radiative transition strength with improved three-nucleon forces are also presented. The experimental results are compared with the predictions of the alpha cluster model and {\it ab initio} calculations.
We consider the role of the internal kinetic energy of bound systems of matter in tests of the Einstein equivalence principle. Using the gravitational sector of the standard model extension, we show that stringent limits on equivalence principle violations in antimatter can be indirectly obtained from tests using bound systems of normal matter. We estimate the bound kinetic energy of nucleons in a range of light atomic species using Green's function Monte Carlo calculations, and for heavier species using a Woods-Saxon model. We survey the sensitivities of existing and planned experimental tests of the equivalence principle, and report new constraints at the level of between a few parts in $10^{6}$ and parts in $10^{8}$ on violations of the equivalence principle for matter and antimatter.
We present recent Green's function Monte Carlo calculations of magnetic moments and M1 transitions in $A \leq 9$ nuclei, which include corrections arising from two-body meson-exchange electromagnetic currents. Two-body effects provide significant corrections to the calculated observables, bringing them in excellent agreement with the experimental data. In particular, we find that two-body corrections are especially large in the A=9, $T=3/2$ systems, in which they account for up to $\sim 20%$ ($\sim 40 %$) of the total predicted value for the $^9$Li ($^9$C) magnetic moment.
We report variational Monte Carlo calculations of single-nucleon momentum distributions for A <= 12 nuclei and nucleon-pair and nucleon-cluster momentum distributions for A <= 8. The wave functions have been generated for a Hamiltonian containing the Argonne v18 two-nucleon and Urbana X three-nucleon potentials. The single-nucleon and nucleon-pair momentum distributions exhibit universal features attributable to the one-pion-exchange tensor interaction. The single-nucleon distributions are broken down into proton and neutron components and spin-up and spin-down components where appropriate. The nucleon-pair momentum distributions are given separately for pp and pn pairs. The nucleon-cluster momentum distributions include dp in 3He, tp and dd in 4He, alpha-d in 6Li, alpha-t in 7Li, and alpha-alpha in 8Be. Detailed tables are provided on-line for download.
Quantum Monte Carlo calculations of electromagnetic moments and transitions are reported for A <= 9 nuclei. The realistic Argonne v18 two-nucleon and Illinois-7 three-nucleon potentials are used to generate the nuclear wave functions. Contributions of two-body meson-exchange current (MEC) operators are included for magnetic moments and M1 transitions. The MEC operators have been derived in both a standard nuclear physics approach and a chiral effective field theory formulation with pions and nucleons including up to one-loop corrections. The two-body MEC contributions provide significant corrections and lead to very good agreement with experiment. Their effect is particularly pronounced in the A=9, T=3/2 systems, in which they provide up to 20% (40 %) of the total predicted value for the 9Li (9C) magnetic moment.
We review the calculation of the equation of state of pure neutron matter using quantum Monte Carlo (QMC) methods. QMC algorithms permit the study of many-body nuclear systems using realistic two- and three-body forces in a nonperturbative framework. We present the results for the equation of state of neutron matter, and focus on the role of three-neutron forces at supranuclear density. We discuss the correlation between the symmetry energy, the neutron star radius and the symmetry energy. We also combine QMC and theoretical models of the three-nucleon interactions, and recent neutron star observations to constrain the value of the symmetry energy and its density dependence.
The lifetime of the J=2+ state in 10C was measured using the Doppler Shift Attenuation Method following the inverse kinematics p(10B,n)10C reaction at 95 MeV. The 2+ state, at 3354 keV, has tau = 219\pm(7)stat \pm(10)sys fs corresponding to a B(E2) # of 8.8(3) e2fm4. This measurement,combined with that recently determined for 10Be (9.2(3) e2fm4), provides a unique challenge to abinitio calculations, testing the structure of these states, including the isospin symmetry of the wave functions. Quantum Monte Carlo calculations using realistic two- and three-nucleon Hamiltonians that reproduce the 10Be B(E2) value generally predict a larger 10C B(E2) probability but with considerable sensitivity to the admixture of different spatial symmetry components in the wave functions, and to the three-nucleon potential used.
We study the elastic scattering of polarized protons from He isotopes. The central and spin-orbit parts of the optical potential are derived using the Glauber theory that can naturally take account of the breakup effect of the He isotopes. Both the differential cross section and the vector analyzing power for $p+^{4,6,8}$He scattering at 71 MeV are in reasonable agreement with experiment. Scattering observables at 300 MeV are predicted. The Pauli blocking effect is examined at 71 MeV.
We present Quantum Monte Carlo calculations of magnetic moments and M1 transitions in $A\le 9$ nuclei which take into account contributions of two-body electromagnetic currents. The Hamiltonian utilized to generate the nuclear wave functions includes the realistic Argonne-{\it v}$_{18}$ two-nucleon and the Illinois-7 three-nucleon interactions. The nuclear two-body electromagnetic currents are derived from a pionful chiral effective field theory including up to one-loop corrections. These currents involve unknown Low Energy Constants which have been fixed so as to reproduce a number of experimental data for the two- and three-nucleon systems, such as $np$ phase shifts and deuteron, triton, and $^3$He magnetic moments. This preliminary study shows that two-body contributions provide significant corrections which are crucial to bring the theory in agreement with the experimental data in both magnetic moments and M1 transitions.
We present calculations of asymptotic normalization coefficients (ANCs) for one-nucleon removals from nuclear states of mass numbers 3 to 9. Our ANCs were computed from variational Monte Carlo solutions to the many-body Schroedinger equation with the combined Argonne v18 two-nucleon and Urbana IX three-nucleon potentials. Instead of computing explicit overlap integrals, we applied a Green's function method that is insensitive to the difficulties of constructing and Monte Carlo sampling the long-range tails of the variational wave functions. This method also allows computation of the ANC at the physical separation energy, even when it differs from the separation energy for the Hamiltonian. We compare our results, which for most nuclei are the first ab initio calculations of ANCs, with existing experimental and theoretical results and discuss further possible applications of the technique.
We use simple analytic considerations and a Monte Carlo calculation of nucleons in a box to argue that the use of Jastrow functions as short-range correlators in the commonly employed two-body-cluster approximation causes significant errors in the matrix elements for double-beta decay. The Jastrow approach appears to agree with others, however, if many-body clusters are included. A careful treatment of the charge-changing analog of the nuclear pair density shows, in addition, that differences between Unitary Correlator Operator Method and Brueckner methods for treating short-range correlations in double-beta are less significant than suggested by previous work.
We present Green's function Monte Carlo calculations of spectroscopic overlaps for $A \leq 7$ nuclei. The realistic Argonne v18 two-nucleon and Illinois-7 three-nucleon interactions are used to generate the nuclear states. The overlap matrix elements are extrapolated from mixed estimates between variational Monte Carlo and Green's function Monte Carlo wave functions. The overlap functions are used to obtain spectroscopic factors and asymptotic normalization coefficients, and they can serve as an input for low-energy reaction calculations.
In order to test {\it ab-initio} calculations of light nuclei, we have remeasured lifetimes in $^{10}$Be using the Doppler Shift Attenuation Method (DSAM) following the $^{7}$Li($^7$Li,$α$)$^{10}$Be reaction at 8 and 10 MeV. The new experiments significantly reduce systematic uncertainties in the DSAM technique. The J$^π$ = $2^+_1$ state at 3.37 MeV has $τ$ = 205${\pm}$(5)$_{stat}\pm$(7)$_{sys}$ fs corresponding to a $B$($E2\downarrow$) of 9.2(3) $e^2$fm$^4$ in broad agreement with many calculations. The J$^π$ = $2^+_2$ state at 5.96 MeV was found to have a $B$($E2\downarrow$) of 0.11(2) $e^2$fm$^4$ and provides a more discriminating test of nuclear models. New Green's Function Monte Carlo (GFMC) calculations for these states and transitions with a number of Hamiltonians are also reported and compared to experiment.
Green's function Monte Carlo calculations of magnetic moments and M1 transitions including two-body meson-exchange current (MEC) contributions are reported for A<=7 nuclei. The realistic Argonne v18 two-nucleon and Illinois-2 three-nucleon potentials are used to generate the nuclear wave functions. The two-body meson-exchange operators are constructed to satisfy the continuity equation with the Argonne v18 potential. The MEC contributions increase the A=3,7 isovector magnetic moments by 16% and the A=6,7 M1 transition rates by 17--34%, bringing them into very good agreement with the experimental data.
We review recent progress in the theory of neutron matter with particular emphasis on its superfluid properties. Results of quantum Monte Carlo calculations of simple and realistic models of uniform superfluid neutron gas are discussed along with those of neutrons interacting in a potential well chosen to approximate neutron-rich oxygen isotopes. The properties of dilute superfluid Fermi gases that may be produced in atom traps, and their relations with neutron matter, are illustrated. The density dependence of the effective interaction between neutrons, used to describe neutron-rich systems with the mean field approximations, is also discussed.
In the past, several efficient methods have been developed to solve the Schroedinger equation for four-nucleon bound states accurately. These are the Faddeev-Yakubovsky, the coupled-rearrangement-channel Gaussian-basis variational, the stochastic variational, the hyperspherical variational, the Green's function Monte Carlo, the no-core shell model and the effective interaction hyperspherical harmonic methods. In this article we compare the energy eigenvalue results and some wave function properties using the realistic AV8' NN interaction. The results of all schemes agree very well showing the high accuracy of our present ability to calculate the four-nucleon bound state.
The energies of $^{3}H$, $^{3}He$, and $^{4}He$ ground states, the ${\frac{3}{2}}^{-}$ and ${\frac{1}{2}}^{-}$ scattering states of $^{5}He$, the ground states of $^{6}He$, $^{6}Li$, and $^{6}Be$ and the $3^{+}$ and $0^{+}$ excited states of $^{6}Li$ have been accurately calculated with the Green's function Monte Carlo method using realistic models of two- and three-nucleon interactions. The splitting of the $A=3$ isospin $T=\frac{1}{2}$ and $A=6$ isospin $T=1$, $J^π = 0^{+}$ multiplets is also studied. The observed energies and radii are generally well reproduced, however, some definite differences between theory and experiment can be identified.