
We present quantum Monte Carlo calculations of the single and twonucleon
momentum distributions in selected nuclei for $A\le16$. We employ local chiral
interactions at nexttonexttoleading order. We find good agreement at low
momentum with the singlenucleon momentum distributions derived for
phenomenological potentials. The same agreement is found for the integrated
twonucleon momentum distributions at low relative momentum $q$ and low
centerofmass momentum $Q$. We provide results for the twonucleon 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 backtoback $(Q=0)$ pairs is
confirmed up to $^{16}$O, and results are compatible with those extracted from
available experimental data.

We report accurate quantum Monte Carlo calculations of nuclei up to $A=16$
based on local chiral two and threenucleon interactions up to
nexttonexttoleading 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 threenucleon interactions. While in
light nuclei the cutoff variation and systematic uncertainties are rather
small, in $^{16}$O these can be significant for large coordinatespace cutoffs.
Overall, we show that chiral interactions constructed to reproduce properties
of very light systems and nucleonnucleon scattering give an excellent
description of binding energies, charge radii, and form factors for all these
nuclei, including openshell systems in $A=6$ and 12.

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 threenucleon interactions up to nexttonexttoleading 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 threebody force for
different cutoffs. The outcomes suggest that local chiral interactions fit to
fewbody observables give a very good description of the groundstate
properties of nuclei up to $^{16}$O, with the exception of one fit for the
softer cutoff which predicts overbinding in larger nuclei.

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 CPviolating phase. Proper interpretation
of the experiments requires reliable theoretical calculations of
neutrinonucleus scattering. In this paper we present calculations of response
functions and cross sections by neutralcurrent 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
vectoraxial interference terms crucial for determining the difference between
neutrino and antineutrino scattering and the CPviolating phase. We find that
the strength and energydependence of twonucleon processes induced by
correlation effects and interaction currents are crucial in providing the most
accurate description of neutrinonucleus scattering in the quasielastic regime.

We present the first ab initio calculations of neutrinoless double beta decay
matrix elements in $A=6$$12$ nuclei using Variational Monte Carlo wave
functions obtained from the Argonne $v_{18}$ twonucleon potential and
Illinois7 threenucleon interaction. We study both light Majorana neutrino
exchange and potentials arising from a large class of multiTeV mechanisms of
lepton number violation. Our results provide benchmarks to be used in testing
manybody methods that can be extended to the heavy nuclei of experimental
interest. In light nuclei we have also studied the impact of twobody short
range correlations and the use of different forms for the transition operators,
such as those corresponding to different orders in chiral effective theory.

Ab initio calculations of the GamowTeller (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}$ twonucleon and Illinois7
threenucleon interactions, and axial manybody currents derived from either
mesonexchange 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 manybody currents
increase by $\sim 2$$3\%$ the onebody 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
shellmodel calculations, may be resolved, at least partially, by correlation
effects.

We examine the groundstate properties of bosons with shortrange two and
threebody interactions where the twobody interactions are tuned to unitarity
and threebody interactions to weakly bound trimers. We determine the
properties of both finite clusters and uniform matter including the energy,
density, two and threebody contacts and the condensate fraction in the bulk,
and demonstrate that these properties are universal in the limit of shortrange
interactions. We find that uniform matter is more bound than threebody
clusters by nearly two orders of magnitude, yet the condensate fraction in the
ground state remains very large, greater than 90\%. Equilibrium properties of
these systems may be experimentally accessible through rapid quenching of
weaklyinteracting boson superfluids.

Local chiral effective field theory interactions have recently been developed
and used in the context of quantum Monte Carlo few and manybody methods for
nuclear physics. In this work, we go over detailed features of local chiral
nucleonnucleon interactions and examine their effect on properties of the
deuteron, paying special attention to the perturbativeness of the expansion. We
then turn to threenucleon 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 wavefunction correlations and
approximations for the two and threebody 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.

The longitudinal and transverse electromagnetic response functions of
$^{12}$C are computed in a "firstprinciples" Green's function Monte Carlo
calculation, based on realistic two and threenucleon interactions and
associated one and twobody 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 reanalysis of the Coulomb sum rule, in which the contributions from the
lowlying $J^\pi\,$=$\, 2^+$, $0^+_2$ (Hoyle), and $4^+$ states in $^{12}$C are
accounted for explicitly in evaluating the total inelastic strength.

We present quantum Monte Carlo calculations of light nuclei, neutron$\alpha$
scattering, and neutron matter using local two and threenucleon ($3N$)
interactions derived from chiral effective field theory up to
nexttonexttoleading order (N$^2$LO). The two undetermined $3N$ lowenergy
couplings are fit to the $^4$He binding energy and, for the first time, to the
spinorbit 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.

Precise measurement of neutrino oscillations, and hence the determination of
their masses demands a quantitative understanding of neutrinonucleus
interactions. To this aim, twobody mesonexchange 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.

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 abinitio calculations reproduce many
lowlying 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 spinisospin, tensor, spinorbit, and
threebody interactions. We present a variety of results including the
lowlying spectra of light nuclei, nuclear form factors, and transition matrix
elements. We also describe lowenergy 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.

The US neutrino community gathered at the Workshop on the Intermediate
Neutrino Program (WINP) at Brookhaven National Laboratory February 46, 2015 to
explore opportunities in neutrino physics over the next five to ten years.
Scientists from particle, astroparticle and nuclear physics participated in the
workshop. The workshop examined promising opportunities for neutrino physics in
the intermediate term, including possible new small to midscale experiments,
US contributions to large experiments, upgrades to existing experiments, R&D
plans and theory. The workshop was organized into two sets of parallel working
group sessions, divided by physics topics and technology. Physics working
groups covered topics on Sterile Neutrinos, Neutrino Mixing, Neutrino
Interactions, Neutrino Properties and Astrophysical Neutrinos. Technology
sessions were organized into Theory, ShortBaseline Accelerator Neutrinos,
Reactor Neutrinos, Detector R&D and Source, Cyclotron and Meson Decay at Rest
sessions.This report summarizes discussion and conclusions from the workshop.

Neutron matter is an intriguing nuclear system with multiple connections to
other areas of physics. Considerable progress has been made over the last two
decades in exploring the properties of pure neutron fluids. Here we begin by
reviewing work done to explore the behavior of very low density neutron matter,
which forms a strongly paired superfluid and is thus similar to cold Fermi
atoms, though at energy scales differing by many orders of magnitude. We then
increase the density, discussing work that ties the study of neutron matter
with the determination of the properties of neutronrich nuclei and
neutronstar crusts. After this, we review the impact neutron matter at even
higher densities has on the massradius relation of neutron stars, thereby
making contact with astrophysical observations.

Ab initio calculations of the quasielastic electromagnetic and neutralweak
response functions of 4He and 12C are carried out for the first time. They are
based on a realistic approach to nuclear dynamics, in which the strong
interactions are described by two and threenucleon potentials and the
electroweak interactions with external fields include one and twobody terms.
The Green's function Monte Carlo method is used to calculate directly the
Laplace transforms of the response functions, and maximumentropy techniques
are employed to invert the resulting imaginarytime correlation functions with
associated statistical errors. The theoretical results, confirmed by experiment
in the electromagnetic case, show that twobody currents generate excess
transverse strength from threshold to the quasielastic to the dip region and
beyond. These findings challenge the conventional picture of quasielastic
inclusive scattering as being largely dominated by singlenucleon knockout
processes.

We present the first Green's function Monte Carlo calculations of light
nuclei with nuclear interactions derived from chiral effective field theory up
to nexttonexttoleading order. Up to this order, the interactions can be
constructed in a local form and are therefore amenable to quantum Monte Carlo
calculations. We demonstrate a systematic improvement with each order for the
binding energies of $A=3$ and $A=4$ systems. We also carry out the first
fewbody tests to study perturbative expansions of chiral potentials at
different orders, finding that higherorder corrections are more perturbative
for softer interactions. Our results confirm the necessity of a threebody
force for correct reproduction of experimental binding energies and radii, and
pave the way for studying few and manynucleon systems using quantum Monte
Carlo methods with chiral interactions.

Recent results in d+Au and p+Pb collisions at RHIC and the LHC provide
evidence for collective expansion and flow of the created medium. We propose a
control set of experiments to directly compare particle emission patterns from
p+Au, d+Au, and He3+Au or t+Au collisions at the same sqrt(sNN). Using Monte
Carlo Glauber we find that a He3 or triton projectile, with a realistic
wavefunction description, induces a significant intrinsic triangular shape to
the initial medium and that, even with viscous damping, this survives into a
significant third order flow moment v3. By comparing systems with one, two, and
three initial hot spots, one can disentangle the effects from the initial
spatial distribution of the deposited energy and viscous damping. These are key
tools to answering the question of how small a droplet of matter is necessary
to form a quarkgluon plasma described by nearly inviscid hydrodynamics.

We investigate the inhomogeneous unitary Fermi gas and use the
longwavelength properties to predict the energies of small clusters of unitary
fermions trapped in harmonic potentials. The large pairing gap and scale
invariance place severe restrictions on the form of the density functional. We
determine the relevant universal constants needed to constrain the functional
from calculations of the bulk in oscillating external potentials. Comparing
with exact Quantum Monte Carlo calculations, we find that the same functional
correctly predicts the lack of shell closures for small clusters of fermions
trapped in harmonic wells as well as their absolute energies. A rapid
convergence to the bulk limit in three dimensions, where the surface to volume
ratio is quite large, is demonstrated. The resulting functional can be tested
experimentally, and is a key ingredient in predicting possible polarized
superfluid phases and the properties of the unitary Fermi gas in optical
lattices.

We present microscopic calculations of light and medium mass nuclei and the
equation of state of symmetric and asymmetric nuclear matter using different
nucleonnucleon forces, including a new Argonne version that has the same
spin/isospin structure as local chiral forces at nexttonexttoleading order
(N2LO). The calculations are performed using Auxiliary Field Diffusion Monte
Carlo (AFDMC) combined with an improved variational wave function. We show that
the AFDMC method can now be used to successfully calculate the energies of very
light to medium mass nuclei as well as the energy of isospinasymmetric nuclear
matter, demonstrating microscopically the quadratic dependence of the energy on
the symmetry energy.

The Facility for Rare Isotope Beams (FRIB) will be a worldleading laboratory
for the study of nuclear structure, reactions and astrophysics. Experiments
with intense beams of rare isotopes produced at FRIB will guide us toward a
comprehensive description of nuclei, elucidate the origin of the elements in
the cosmos, help provide an understanding of matter in neutron stars, and
establish the scientific foundation for innovative applications of nuclear
science to society. FRIB will be essential for gaining access to key regions of
the nuclear chart, where the measured nuclear properties will challenge
established concepts, and highlight shortcomings and needed modifications to
current theory. Conversely, nuclear theory will play a critical role in
providing the intellectual framework for the science at FRIB, and will provide
invaluable guidance to FRIB's experimental programs. This article overviews the
broad scope of the FRIB theory effort, which reaches beyond the traditional
fields of nuclear structure and reactions, and nuclear astrophysics, to explore
exciting interdisciplinary boundaries with other areas.
\keywords{Nuclear Structure and Reactions. Nuclear
Astrophysics. Fundamental Interactions. High Performance
Computing. Rare Isotopes. Radioactive Beams.

An ab initio calculation of the sum rules of the neutral weak response
functions in $^{12}$C is reported, based on a realistic Hamiltonian, including
two and threenucleon potentials, and on realistic currents, consisting of
one and twobody terms. We find that the sum rules of the response functions
associated with the longitudinal and transverse components of the (spacelike)
neutral current are largest and that a significant portion ($\simeq 30$ %) of
the calculated strength is due to twobody terms. This fact may have
implications for the MiniBooNE and other neutrino quasielastic scattering data
on nuclei.

We report variational Monte Carlo calculations of singlenucleon momentum
distributions for A <= 12 nuclei and nucleonpair and nucleoncluster momentum
distributions for A <= 8. The wave functions have been generated for a
Hamiltonian containing the Argonne v18 twonucleon and Urbana X threenucleon
potentials. The singlenucleon and nucleonpair momentum distributions exhibit
universal features attributable to the onepionexchange tensor interaction.
The singlenucleon distributions are broken down into proton and neutron
components and spinup and spindown components where appropriate. The
nucleonpair momentum distributions are given separately for pp and pn pairs.
The nucleoncluster momentum distributions include dp in 3He, tp and dd in 4He,
alphad in 6Li, alphat in 7Li, and alphaalpha in 8Be. Detailed tables are
provided online for download.

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
manybody nuclear systems using realistic two and threebody forces in a
nonperturbative framework. We present the results for the equation of state of
neutron matter, and focus on the role of threeneutron 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 threenucleon interactions, and recent neutron star observations to
constrain the value of the symmetry energy and its density dependence.

An "ab initio" calculation of the Carbon12 elastic form factor, and sum
rules of longitudinal and transverse response functions measured in inclusive
(e,e') scattering, is reported, based on realistic nuclear potentials and
electromagnetic currents. The longitudinal elastic form factor and sum rule are
found to be in satisfactory agreement with available experimental data. A
direct comparison between theory and experiment is difficult for the transverse
sum rule. However, it is shown that the calculated one has large contributions
from twobody currents, indicating that these mechanisms lead to a significant
enhancement of the quasielastic transverse response. This fact may have
implications for the anomaly observed in recent neutrino quasielastic
chargechanging scattering data off Carbon12.

We give the first selfconsistent calculation of the effect of the scattered
neutrino halo on flavor evolution in supernovae. Our example case is an ONeMg
core collapse supernova neutronization neutrino burst. We find that the
addition of the halo neutrinos produces qualitative and quantitative changes in
the final flavor states of neutrinos. We also find that the halo neutrinos
produce a novel distortion of the neutrino flavor swap. Our results provide
strong motivation for tackling the full multidimensional and
compositiondependent aspects of this problem in the future.