BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
M. Piarulli 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.
Searches for neutrinoless-double beta decay rates are crucial in addressing questions within fundamental symmetries and neutrino physics. The rates of these decays depend not only on unknown parameters associated with neutrinos, but also on nuclear properties. In order to reliably extract information about the neutrino, one needs an accurate treatment of the complex many-body dynamics of the nucleus. Neutrinoless-double beta decays take place at momentum transfers on the order of 100 MeV/$c$ and require both nuclear electroweak vector and axial current matrix elements. Muon capture, a process in the same momentum transfer regime, has readily available experimental data to validate these currents. In this work, we present results of {\it ab initio} calculations of partial muon capture rates for $^3$He and $^6$Li nuclei using variational and Green's Function Monte Carlo computational methods. We estimate the impact of the three-nucleon interactions, the cutoffs used to regularize two-nucleon ($2N$) interactions, and the energy range of $2N$ scattering data used to fit these interactions.
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.
Experimental form factors of the hydrogen and helium isotopes, extracted from an up-to-date global analysis of cross sections and polarization observables measured in elastic electron scattering from these systems, are compared to predictions obtained in three different theoretical approaches: the first is based on realistic interactions and currents, including relativistic corrections (labeled as the conventional approach); the second relies on a chiral effective field theory description of the strong and electromagnetic interactions in nuclei (labeled $χ$EFT); the third utilizes a fully relativistic treatment of nuclear dynamics as implemented in the covariant spectator theory (labeled CST). For momentum transfers below $Q \lesssim 5$ fm$^{-1}$ there is satisfactory agreement between experimental data and theoretical results in all three approaches. However, at $Q \gtrsim 5$ fm$^{-1}$, particularly in the case of the deuteron, a relativistic treatment of the dynamics, as is done in the CST, is necessary. The experimental data on the deuteron $A$ structure function extend to $Q \simeq 12$ fm$^{-1}$, and the close agreement between these data and the CST results suggests that, even in this extreme kinematical regime, there is no evidence for new effects coming from quark and gluon degrees of freedom at short distances.
We construct a coordinate-space chiral potential, including $Δ$-isobar intermediate states in its two-pion-exchange component. The contact interactions entering at next-to-leading and next-to-next-to-next-to-leading orders ($Q^2$ and $Q^4$, respectively, $Q$ denoting generically the low momentum scale) are rearranged by Fierz transformations to yield terms at most quadratic in the relative momentum operator of the two nucleons. The low-energy constants multiplying these contact interactions are fitted to the 2013 Granada database, consisting of 2309 $pp$ and 2982 $np$ data (including, respectively, 148 and 218 normalizations) in the laboratory-energy range 0--300 MeV. For the total 5291 $pp$ and $np$ data in this range, we obtain a $χ^2$/datum of roughly 1.3 for a set of three models characterized by long- and short-range cutoffs, $R_{\rm L}$ and $R_{\rm S}$ respectively, ranging from $(R_{\rm L},R_{\rm S})=(1.2,0.8)$ fm down to $(0.8,0.6)$ fm. The long-range (short-range) cutoff regularizes the one- and two-pion exchange (contact) part of the potential.
The muon capture reactions 2H(μ^-,ν_μ)nn and 3He(μ^-,ν_μ)3H are studied with conventional or chiral realistic potentials and consistent weak currents. The initial and final A=2 and 3 nuclear wave functions are obtained from the Argonne v18 or chiral N3LO two-nucleon potential, in combination with, respectively, the Urbana IX or chiral N2LO three-nucleon potential in the case of A=3. The weak current consists of polar- and axial-vector components. The former are related to the isovector piece of the electromagnetic current via the conserved-vector-current hypothesis. These and the axial currents are derived either in a meson-exchange or in a chiral effective field theory (chiEFT) framework. There is one parameter (either the N-to-Δaxial coupling constant in the meson-exchange model, or the strength of a contact term in the chiEFT model) which is fixed by reproducing the Gamow-Teller matrix element in tritium beta-decay. The model dependence relative to the adopted interactions and currents (and cutoff sensitivity in the chiEFT currents) is weak, resulting in total rates of 392.0 +/- 2.3 Hz for A=2, and 1484 +/- 13 Hz for A=3, where the spread accounts for this model dependence.
This work investigates the muon capture reactions 2H(μ^-,ν_μ)nn and 3He(μ^-,ν_μ)3H and the contribution to their total capture rates arising from the axial two-body currents obtained imposing the partially-conserved-axial-current (PCAC) hypothesis. The initial and final A=2 and 3 nuclear wave functions are obtained from the Argonne v_{18} two-nucleon potential, in combination with the Urbana IX three-nucleon potential in the case of A=3. The weak current consists of vector and axial components derived in chiral effective field theory. The low-energy constant entering the vector (axial) component is determined by reproducting the isovector combination of the trinucleon magnetic moment (Gamow-Teller matrix element of tritium beta-decay). The total capture rates are 393.1(8) s^{-1} for A=2 and 1488(9) s^{-1} for A=3, where the uncertainties arise from the adopted fitting procedure.