
2
distributions over a wide range of nuclei from 2H up to
12C, using the Norfolk NN and 3N(NV2+3) forces [25–
29]. These results feature new calculations of the pair
density as a function of both the pair separation and
pair center-of-mass, and calculations of the two-body
momentum distribution coming from short- and long-
range pairs differentiated by a pair separation bound-
ary. The full set of calculations is accessible in graphical
and tabular forms online at www.phy.anl.gov/theory/
research/QMCresults.html.
The paper is structured as follows: a brief review of
Norfolk interactions is given in Sec. II. In Sec. III we
present results for the one- and two-body densities cal-
culated for 3H, 3,4,8He, 6,7Li, 9Be, 10B, and 12C. The pair
density as a function of both the pair separation and pair
center-of-mass is presented for 4He and 12C. In Sec. IV
the results for the one- and two-body momentum distri-
bution are provided for 3H, 3,4,8He, 6,7Li, 9Be, 10B, and
12C. Results for momentum distributions as functions of
the relative momentum and center-of-mass momentum
without and with pair separation boundary are displayed
for 4He and 12C. Additional results are available online.
II. NORFOLK MANY-BODY INTERACTIONS
The Norfolk interactions are obtained from a chiral
effective field theory (χEFT) that uses pions, nucleons
and ∆’s as fundamental degrees of freedom, and con-
sists of long-range parts mediated by one- and two-
pion exchange, and contact terms specified by unknown
low-energy constants (LECs). The LECs entering the
NN contact interactions are constrained to reproduce
NN scattering data from the most recent and up-to-
date database collected by the Granada group [30–32].
The contact terms are regularized via a Gaussian cut-
off function with RSas the Gaussian parameter [25–
27]. The divergences at high-value of momentum trans-
fer in the pion-range operators are removed via a spe-
cial radial function characterized by the cutoff RL[25–
27]. There are two classes of NV2 potentials. Class I
(II) has been fitted to data up to 125 MeV (200 MeV).
For each class, two combinations of short- and long-range
regulators have been used, namely (RS,RL)=(0.8, 1.2)
fm (models NV2-Ia and NV2-IIa) and (RS,RL)=(0.7,
1.0) fm (models NV2-Ib and NV2-IIb). Class I (II) fits
about 2700 (3700) data points with a χ2/datum <
∼1.1
(<
∼1.4) [25, 26]. The short-range component of the 3N
interactions is parametrized in terms of two LECs, cD
and cE. In the first generation of Norfolk potentials
(NV2+3-Ia/b and NV2+3-IIa/b), these LECs have been
determined by simultaneously reproducing the experi-
mental trinucleon ground-state energies and nd doublet
scattering length [33]. Within the χEFT framework, cD
is related to the LEC entering the axial two-body con-
tact current [34–36]. This allows one to adopt a dif-
ferent strategy to constrain cDand cE. In particular,
in Ref. [27] they have been constrained to reproduce the
trinucleon binding energies and the empirical value of the
Gamow-Teller matrix element in tritium βdecay. Nor-
folk models that use this fitting procedure are designated
with a ‘*’ namely, NV2+3-Ia*/b* and NV2+3-IIa*/b*.
These interactions have been recently employed in the
VMC and Green’s function Monte Carlo (GFMC) ap-
proaches [37, 38] to calculate energies [33], charge radii
and electromagnetic form factors [38], beta-decay transi-
tions [27, 39, 40], neutrinoless double beta-decay [41, 42]
of light nuclei, beta decay spectra [43], muon-capture
rates [44] and with the auxiliary field diffusion Monte
Carlo (AFDMC) [38] to study the equation of state of
pure neutron matter [45, 46].
III. DENSITY DISTRIBUTIONS
The one- and two-body densities are evaluated as sim-
ple δ-function expectation values given by
ρN(r) = 1
4πr2Ψ
X
i
PNiδ(r− |ri−Rcm|)
Ψ,(1)
ρNN (r) = 1
4πr2Ψ
X
i<j
PNiPNjδ(r− |ri−rj|)
Ψ,(2)
where PNirepresents the projector operator onto protons
(+) or neutrons (−) defined as PNi= (1 ±τzi)/2, riis
the position of nucleon iand Rcm is the coordinate of
the center of mass.
A detailed survey of one- and two-body densities have
been calculated for a variety of nuclei in the range
A= 2 −12 using variational Monte Carlo wave func-
tions developed for the AV18+UX and the Norfolk local
chiral interactions. The corresponding tables and fig-
ures are available online at www.phy.anl.gov/theory/
research/density/, for the one-nucleon densities, and
at www.phy.anl.gov/theory/research/density2/, for
the two-nucleon densities.
A. One-body density results
In Fig. 1 we present the neutron and proton densities
calculated for 3H, 3,4,8He, 6,7Li, 9Be, 10B, and 12C us-
ing the AV18+UX and the NV2+3-Ia, NV2+3-Ia*, and
NV2+3-IIb* local chiral interactions. Additional den-
sities for 2H, 6He, 8,9Li, 8,10,12Be, 11B and 10,11C may
be found in the online tables, as well as results for the
NV2+3-Ib* and NV2+3-IIa* interactions. We also give
neutron and proton rms radii there.
The VMC wave functions are treated as states of
unique isospin T. Thus for N=Znuclei, proton and
neutron densities are the same and only proton densi-
ties are given in the online tables. However, the wave
functions for nuclei with T > 0 can be different for dif-
ferent isospin projections Tz, so mirror nuclei are not
isospin symmetric. This allows the proton-rich nuclei to