While mean-field approximations, such as the nuclear shell model, provide a good description of many bulk nuclear properties, they fail to capture the important effects of nucleon–nucleon ...correlations such as the short-distance and high-momentum components of the nuclear many-body wave function1. Here, we study these components using the effective pair-based generalized contact formalism2,3 and ab initio quantum Monte Carlo calculations of nuclei from deuteron to 40Ca (refs. 4–6). We observe a universal factorization of the many-body nuclear wave function at short distance into a strongly interacting pair and a weakly interacting residual system. The residual system distribution is consistent with that of an uncorrelated system, showing that short-distance correlation effects are predominantly embedded in two-body correlations. Spin- and isospin-dependent ‘nuclear contact terms’ are extracted in both coordinate and momentum space for different realistic nuclear potentials. The contact coefficient ratio between two different nuclei shows very little dependence on the nuclear interaction model. These findings thus allow extending the application of mean-field approximations to short-range correlated pair formation by showing that the relative abundance of short-range pairs in the nucleus is a long-range (that is, mean field) quantity that is insensitive to the short-distance nature of the nuclear force.Effects of nucleon–nucleon correlations are studied with the generalized contact formalism and ab initio quantum Monte Carlo calculations. For nuclei from deuteron to 40Ca, the many-body nuclear wave function is shown to factorize at short distances.
In recent years local chiral interactions have been derived and implemented in quantum Monte Carlo methods in order to test to what extent the chiral effective field theory framework impacts our ...knowledge of few- and many-body systems. In this Letter, we present Green's function Monte Carlo calculations of light nuclei based on the family of local two-body interactions presented by our group in a previous paper in conjunction with chiral three-body interactions fitted to bound- and scattering-state observables in the three-nucleon sector. These interactions include Δ intermediate states in their two-pion-exchange components. We obtain predictions for the energy levels and level ordering of nuclei in the mass range A=4-12, accurate to ≤2% of the binding energy, in very satisfactory agreement with experimental data.
The muon capture reaction
μ
−
+
d
→
n
+
n
+
ν
μ
in the doublet hyperfine state is studied using nuclear potentials and consistent currents derived in the chiral effective field theory, which are ...local and expressed in coordinate space (the so-called Norfolk models). Only the largest contribution due to the
1
S
0
nn
scattering state is considered. Particular attention is given to the estimate of theoretical uncertainty, for which four sources have been identified: 1) the model dependence, 2) the chiral-order convergence for the weak nuclear current, 3) the uncertainty in the single-nucleon axial form factor, and 4) the numerical technique adopted to solve the bound and scattering
A
= 2 systems. This last source of uncertainty has turned out to be essentially negligible. For the
1
S
0
doublet muon capture rate
Γ
D
(
S
0
1
)
, we obtain
Γ
D
(
S
0
1
)
=
255.8
(
0.6
)
(
4.4
)
(
2.9
)
s
−1
, where the three errors come from the first three sources of uncertainty. The value for
Γ
D
(
S
0
1
)
obtained within this local chiral framework is compared with previous calculations and found in very good agreement.
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 v18 interaction and two of its simplified versions, and four of the new Norfolk Δ -full chiral effective field theory potentials. Primarily because of the advancements in the auxiliary-field diffusion Monte Carlo, we observe good agreement among the three many-body techniques up to nuclear saturation density—the maximum difference in the energy per particle is within 1.5 MeV for all the potentials we consider. At higher densities, the divergences become more important, and are mainly driven by the Fermi hypernetted chain/single-operator calculations. We also study the connection between nucleon-nucleon scattering data and the energy per particle of pure neutron matter. Our results suggest that fitting to higher-energy nucleon-nucleon scattering helps reduce the spread of energies among the models.
In this work, we report quantum Monte Carlo calculations of weak transitions in A≤10 nuclei, based on the Norfolk two- and three-nucleon chiral interactions, and associated one- and two-body axial ...currents. Furthermore, 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 8Li, 8B, and 8He β decay. 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.
In this study, 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 v18 and the Argonne v'6 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 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.