We present quantum Monte Carlo calculations of the single- and two-nucleon momentum distributions in selected nuclei for $A\le16$. We employ local chiral interactions at next-to-next-to-leading order. We find good agreement at low momentum with the single-nucleon momentum distributions derived for phenomenological potentials. The same agreement is found for the integrated two-nucleon momentum distributions at low relative momentum $q$ and low center-of-mass momentum $Q$. We provide results for the two-nucleon 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 back-to-back $(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 three-nucleon interactions up to next-to-next-to-leading 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 three-nucleon interactions. While in light nuclei the cutoff variation and systematic uncertainties are rather small, in $^{16}$O these can be significant for large coordinate-space cutoffs. Overall, we show that chiral interactions constructed to reproduce properties of very light systems and nucleon-nucleon scattering give an excellent description of binding energies, charge radii, and form factors for all these nuclei, including open-shell 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 three-nucleon interactions up to next-to-next-to-leading 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 three-body force for different cutoffs. The outcomes suggest that local chiral interactions fit to few-body observables give a very good description of the ground-state 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 CP-violating phase. Proper interpretation of the experiments requires reliable theoretical calculations of neutrino-nucleus scattering. In this paper we present calculations of response functions and cross sections by neutral-current 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 vector-axial interference terms crucial for determining the difference between neutrino and anti-neutrino scattering and the CP-violating phase. We find that the strength and energy-dependence of two-nucleon processes induced by correlation effects and interaction currents are crucial in providing the most accurate description of neutrino-nucleus 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}$ two-nucleon potential and Illinois-7 three-nucleon interaction. We study both light Majorana neutrino exchange and potentials arising from a large class of multi-TeV mechanisms of lepton number violation. Our results provide benchmarks to be used in testing many-body methods that can be extended to the heavy nuclei of experimental interest. In light nuclei we have also studied the impact of two-body 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 Gamow-Teller (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}$ two-nucleon and Illinois-7 three-nucleon interactions, and axial many-body currents derived from either meson-exchange 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 many-body currents increase by $\sim 2$--$3\%$ the one-body 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 shell-model calculations, may be resolved, at least partially, by correlation effects.

We examine the ground-state properties of bosons with short-range two- and three-body interactions where the two-body interactions are tuned to unitarity and three-body interactions to weakly bound trimers. We determine the properties of both finite clusters and uniform matter including the energy, density, two- and three-body contacts and the condensate fraction in the bulk, and demonstrate that these properties are universal in the limit of short-range interactions. We find that uniform matter is more bound than three-body 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 weakly-interacting boson superfluids.

Local chiral effective field theory interactions have recently been developed and used in the context of quantum Monte Carlo few- and many-body methods for nuclear physics. In this work, we go over detailed features of local chiral nucleon-nucleon interactions and examine their effect on properties of the deuteron, paying special attention to the perturbativeness of the expansion. We then turn to three-nucleon 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 wave-function correlations and approximations for the two- and three-body 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 "first-principles" Green's function Monte Carlo calculation, based on realistic two- and three-nucleon interactions and associated one- and two-body 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 re-analysis of the Coulomb sum rule, in which the contributions from the low-lying $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 three-nucleon ($3N$) interactions derived from chiral effective field theory up to next-to-next-to-leading order (N$^2$LO). The two undetermined $3N$ low-energy couplings are fit to the $^4$He binding energy and, for the first time, to the spin-orbit 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 neutrino-nucleus interactions. To this aim, two-body meson-exchange 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 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.

The US neutrino community gathered at the Workshop on the Intermediate Neutrino Program (WINP) at Brookhaven National Laboratory February 4-6, 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 mid-scale 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, Short-Baseline 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 neutron-rich nuclei and neutron-star crusts. After this, we review the impact neutron matter at even higher densities has on the mass-radius relation of neutron stars, thereby making contact with astrophysical observations.

Ab initio calculations of the quasi-elastic electromagnetic and neutral-weak 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 three-nucleon potentials and the electroweak interactions with external fields include one- and two-body terms. The Green's function Monte Carlo method is used to calculate directly the Laplace transforms of the response functions, and maximum-entropy techniques are employed to invert the resulting imaginary-time correlation functions with associated statistical errors. The theoretical results, confirmed by experiment in the electromagnetic case, show that two-body currents generate excess transverse strength from threshold to the quasi-elastic to the dip region and beyond. These findings challenge the conventional picture of quasi-elastic inclusive scattering as being largely dominated by single-nucleon 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 next-to-next-to-leading 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 few-body tests to study perturbative expansions of chiral potentials at different orders, finding that higher-order corrections are more perturbative for softer interactions. Our results confirm the necessity of a three-body force for correct reproduction of experimental binding energies and radii, and pave the way for studying few- and many-nucleon 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 quark-gluon plasma described by nearly inviscid hydrodynamics.

We investigate the inhomogeneous unitary Fermi gas and use the long-wavelength 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 nucleon-nucleon forces, including a new Argonne version that has the same spin/isospin structure as local chiral forces at next-to-next-to-leading 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 isospin-asymmetric nuclear matter, demonstrating microscopically the quadratic dependence of the energy on the symmetry energy.

The Facility for Rare Isotope Beams (FRIB) will be a world-leading 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 three-nucleon potentials, and on realistic currents, consisting of one- and two-body terms. We find that the sum rules of the response functions associated with the longitudinal and transverse components of the (space-like) neutral current are largest and that a significant portion ($\simeq 30$ %) of the calculated strength is due to two-body terms. This fact may have implications for the MiniBooNE and other neutrino quasi-elastic scattering data on nuclei.

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.

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.

An "ab initio" calculation of the Carbon-12 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 two-body currents, indicating that these mechanisms lead to a significant enhancement of the quasi-elastic transverse response. This fact may have implications for the anomaly observed in recent neutrino quasi-elastic charge-changing scattering data off Carbon-12.

We give the first self-consistent calculation of the effect of the scattered neutrino halo on flavor evolution in supernovae. Our example case is an O-Ne-Mg 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 composition-dependent aspects of this problem in the future.