We report a quantum Monte Carlo calculation of the equation of state of symmetric nuclear matter using local interactions derived from chiral effective field theory up to next-to-next-to-leading order fit to few-body observables only. The empirical saturation density and energy are well reproduced within statistical and systematic uncertainties. We have also derived the symmetry energy as a function of the density, finding good agreement with available experimentally derived constraints at saturation and twice saturation density. We find that the corresponding pressure is also in excellent agreement with recent constraints extracted from gravitational waves of the neutron-star merger GW170817 by the LIGO-Virgo detection.

Neutrinos emitted deep within a supernova explosion experience a self-induced index of refraction. In the stationary, one-dimensional (1D) supernova "bulb model", this self-induced refraction can lead to a collective flavor transformation which is coherent among different neutrino momentum modes. Such collective oscillations can produce partial swaps of the energy spectra of different neutrino species as the neutrinos stream away from the proto-neutron star. However, it has been demonstrated that the spatial symmetries (such as the spherical symmetry in the bulb model) can be broken spontaneously by collective neutrino oscillations in multi-dimensional models. Using a stationary, 2D neutrino ring model we demonstrate that there exist two limiting scenarios where collective oscillations may occur. In one limit, the collective flavor transformation begins at a radius with relatively high neutrino densities and develops small-scale flavor structures. The loss of the spatial correlation in the neutrino flavor field results in similar (average) energy spectra for the anti-neutrinos of almost all energies and the neutrinos of relatively high energies. In the other limit, the flavor transformation starts at a radius where the neutrino densities are smaller (e.g., due to the suppression of the high matter density near the proto-neutron star). Although the spatial symmetry is broken initially, it is restored as the neutrino densities decrease, and the neutrinos of different flavors partially swap their energy spectra as in the 1D bulb model. This finding may have interesting ramifications in other aspects of supernova physics.

We use the quantum Monte Carlo (QMC) techniques to calculate the static structure function $S(q)$ of a one-component ion lattice and use it to calculate the thermal conductivity $\kappa$ of high-density solid matter expected in the neutron star crust. By making detailed comparisons with the results for the thermal conductivity obtained using standard techniques based on the one-phonon approximation (OPA) valid at low temperature, and the multi-phonon harmonic approximation expected to be valid over a wide range of temperatures, we asses the temperature regime where $S(q)$ from QMC can be used directly to calculate $\kappa$. We also compare the QMC results to those obtained using classical Monte Carlo to quantitatively asses the magnitude of the quantum corrections. We find that quantum effects became relevant for the calculation of $\kappa$ at temperature $T \lesssim 0.3 ~\Omega_\mathrm{P}$, where $\Omega_\mathrm{P}$ is the ion plasma frequency. At $T \simeq 0.1 ~\Omega_\mathrm{P}$ the quantum effects suppress $\kappa$ by about $30\%$. The comparison with the results of the OPA indicates that dynamical information beyond the static structure is needed when $T \lesssim 0.1~ \Omega_\mathrm{P}$. These quantitative comparisons help to establish QMC as a viable technique to calculate $\kappa$ at moderate temperatures in the range $T=0.1-1~\Omega_\mathrm{P}$ of relevance to the study of accreting neutron stars. This finding is especially important because QMC is the only viable technique so far for calculating $\kappa$ in multi-component systems at low-temperatures.

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.

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.

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.

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.

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.

We study the spin response of cold dense neutron matter in the limit of zero momentum transfer, and show that the frequency dependence of the long-wavelength spin response is well constrained by sum-rules and the asymptotic behavior of the two-particle response at high frequency. The sum-rules are calculated using Auxiliary Field Diffusion Monte Carlo technique and the high frequency two-particle response is calculated for several nucleon-nucleon potentials. At nuclear saturation density, the sum-rules suggest that the strength of the spin response peaks at $\omega \simeq$ 40--60 MeV, decays rapidly for $\omega \geq $100 MeV, and has a sizable strength below 40 MeV. This strength at relatively low energy may lead to enhanced neutrino production rates in dense neutron-rich matter at temperatures of relevance to core-collapse supernova.

We present a new formulation of Lagrangian perturbation theory which allows accurate predictions of the real- and redshift-space correlation functions of the mass field and dark matter halos. Our formulation involves a non-perturbative resummation of Lagrangian perturbation theory and indeed can be viewed as a partial resummation of the formalism of Matsubara (2008a,b) in which we keep exponentiated all of the terms which tend to a constant at large separation. One of the key features of our method is that we naturally recover the Zel'dovich approximation as the lowest order of our expansion for the matter correlation function. We compare our results against a suite of N-body simulations and obtain good agreement for the correlation functions in real-space and for the monopole correlation function in redshift space. The agreement becomes worse for higher multipole moments of the redshift-space, halo correlation function. Our formalism naturally includes non-linear bias and explains the strong bias-dependence of the multipole moments of the redshift-space correlation function seen in N-body simulations.

We calculate the equation of state of neutron matter with realistic two- and three-nucleon interactions using quantum Monte Carlo techniques, and illustrate that the short-range three-neutron interaction determines the correlation between neutron matter energy at nuclear saturation density and higher densities relevant to neutron stars. Our model also makes an experimentally testable prediction for the correlation between the nuclear symmetry energy and its density dependence -- determined solely by the strength of the short-range terms in the three neutron force. The same force provides a significant constraint on the maximum mass and radius of neutron stars.

We argue that the small fraction of neutrinos that undergo direction-changing scattering outside of the neutrinosphere could have significant influence on neutrino flavor transformation in core-collapse supernova environments. We show that the standard treatment for collective neutrino flavor transformation is adequate at late times, but could be inadequate in the crucial shock revival/explosion epoch of core-collapse supernovae, where the potentials that govern neutrino flavor evolution are affected by the scattered neutrinos. Taking account of this effect, and the way it couples to entropy and composition, will require a new paradigm in supernova modeling.

We show that the bump in the electron number density profile at the base of the hydrogen envelope in O-Ne-Mg core-collapse supernovae causes an interesting interplay between neutrino-electron and neutrino-neutrino forward scattering effects in the flavor evolution of low-energy nu_e in the neutronization burst. The bump allows a significant fraction of the low-energy nu_e to survive by rendering their flavor evolution nonadiabatic. Increasing the luminosity of the neutronization burst shifts the bump-affected nu_e to lower energy with reduced survival probability. Similarly, lowering the luminosity shifts the bump-affected neutrinos to higher energies. While these low energy neutrinos lie near the edge of detectability, the population of bump-affected neutrinos has direct influence on the spectral swap formation in the neutrino signal at higher energies.

We investigate the effect of matter density fluctuations on supernova collective neutrino flavor oscillations. In particular, we use full multi-angle, 3-flavor, self-consistent simulations of the evolution of the neutrino flavor field in the envelope of an O-Ne-Mg core collapse supernova at shock break-out (neutrino neutronization burst) to study the effect of the matter density "bump" left by the He-burning shell. We find a seemingly counterintuitive increase in the overall electron neutrino survival probability created by this matter density feature. We discuss this behavior in terms of the interplay between the matter density profile and neutrino collective effects. While our results give new insights into this interplay, they also suggest an immediate consequence for supernova neutrino burst detection: it will be difficult to use a burst signal to extract information on fossil burning shells or other fluctuations of this scale in the matter density profile. Consistent with previous studies, our results also show that the interplay of neutrino self-coupling and matter fluctuation could cause a significant increase in the electron neutrino survival probability at very low energy

We report results of the first 3-by-3 "multi-angle" simulation of the evolution of neutrino flavor in the core collapse supernova environment. In particular, we follow neutrino flavor transformation in the neutronization neutrino burst of an O-Ne-Mg core collapse event. Though in a qualitative sense our results are consistent with those obtained in 3-by-3 single-angle simulations, at least in terms of neutrino mass hierarchy dependence, performing multi-angle calculations is found to reduce the adiabaticity of flavor evolution in the normal neutrino mass hierarchy, resulting in lower swap energies. Differences between single-angle and multi-angle results are largest for the normal neutrino mass hierarchy. Our simulations also show that current uncertainties in the measured mass-squared and mixing angle parameters translate into uncertainties in neutrino swap energies. Our results show that at low theta-13 it may be difficult to resolve the neutrino mass hierarchy using the O-Ne-Mg neutronization neutrino burst.

Connecting cosmological simulations to real-world observational programs is often complicated by a mismatch in geometry: while surveys often cover highly irregular cosmological volumes, simulations are customarily performed in a periodic cube. We describe a technique to remap this cube into elongated box-like shapes that are more useful for many applications. The remappings are one-to-one, volume-preserving, keep local structures intact, and involve minimal computational overhead.

We present a set of ultra-large particle-mesh simulations of the LyA forest targeted at understanding the imprint of baryon acoustic oscillations (BAO) in the inter-galactic medium. We use 9 dark matter only simulations which can, for the first time, simultaneously resolve the Jeans scale of the intergalactic gas while covering the large volumes required to adequately sample the acoustic feature. Mock absorption spectra are generated using the fluctuating Gunn-Peterson approximation which have approximately correct flux probability density functions (PDFs) and small-scale power spectra. On larger scales there is clear evidence in the redshift space correlation function for an acoustic feature, which matches a linear theory template with constant bias. These spectra, which we make publicly available, can be used to test pipelines, plan future experiments and model various physical effects. As an illustration we discuss the basic properties of the acoustic signal in the forest, the scaling of errors with noise and source number density, modified statistics to treat mean flux evolution and misestimation, and non-gravitational sources such as fluctuations in the photo-ionizing background and temperature fluctuations due to HeII reionization.

Recently a number of analytic prescriptions for computing the non-linear matter power spectrum have appeared in the literature. These typically involve resummation or closure prescriptions which do not have a rigorous error control, thus they must be compared with numerical simulations to assess their range of validity. We present a direct side-by-side comparison of several of these analytic approaches, using a suite of high-resolution N-body simulations as a reference, and discuss some general trends. All of the analytic results correctly predict the behavior of the power spectrum at the onset of non-linearity, and improve upon a pure linear theory description at very large scales. All of these theories fail at sufficiently small scales. At low redshift the dynamic range in scale where perturbation theory is both relevant and reliable can be quite small. We also compute for the first time the 2-loop contribution to standard perturbation theory for CDM models, finding improved agreement with simulations at large redshift. At low redshifts however the 2-loop term is larger than the 1-loop term on quasi-linear scales, indicating a breakdown of the perturbation expansion. Finally, we comment on possible implications of our results for future studies.

We discuss a new kind of astrophysical transport problem: the coherent evolution of neutrino flavor in core collapse supernovae. Solution of this problem requires a numerical approach which can simulate accurately the quantum mechanical coupling of intersecting neutrino trajectories and the associated nonlinearity which characterizes neutrino flavor conversion. We describe here the two codes developed to attack this problem. We also describe the surprising phenomena revealed by these numerical calculations. Chief among these is that the nonlinearities in the problem can engineer neutrino flavor transformation which is dramatically different than in standard Mikheyev-Smirnov-Wolfenstein treatments. This happens even though the neutrino mass-squared differences are measured to be small, and even when neutrino self-coupling is sub-dominant. Our numerical work has revealed potential signatures which, if detected in the neutrino burst from a Galactic core collapse event, could reveal heretofore unmeasurable properties of the neutrinos, such as the mass hierarchy and vacuum mixing angle theta_13.

We present results of 3-neutrino flavor evolution simulations for the neutronization burst from an O-Ne-Mg core-collapse supernova. We find that nonlinear neutrino self-coupling engineers a single spectral feature of stepwise conversion in the inverted neutrino mass hierarchy case and in the normal mass hierarchy case, a superposition of two such features corresponding to the vacuum neutrino mass-squared differences associated with solar and atmospheric neutrino oscillations. These neutrino spectral features offer a unique potential probe of the conditions in the supernova environment and may allow us to distinguish between O-Ne-Mg and Fe core-collapse supernovae.

We examine a phenomenon recently predicted by numerical simulations of supernova neutrino flavor evolution: the swapping of supernova $\nu_e$ and $\nu_{\mu,\tau}$ energy spectra below (above) energy $\EC$ for the normal (inverted) neutrino mass hierarchy. We present the results of large-scale numerical calculations which show that in the normal neutrino mass hierarchy case, $\EC$ decreases as the assumed $\nu_e\rightleftharpoons\nu_{\mu,\tau}$ effective $2\times 2$ vacuum mixing angle ($\simeq \theta_{1 3}$) is decreased. However, these calculations also indicate that $\EC$ is essentially independent of the vacuum mixing angle in the inverted neutrino mass hierarchy case. With a good neutrino signal from a future Galactic supernova, the above results could be used to determine the neutrino mass hierarchy even if $\theta_{13}$ is too small to be detected in terrestrial neutrino oscillation experiments.

We study the flavor evolution of a dense gas initially consisting of pure mono-energetic $\nu_e$ and $\bar\nu_e$. Using adiabatic invariants and the special symmetry in such a system we are able to calculate the flavor evolution of the neutrino gas for the cases with slowly decreasing neutrino number densities. These calculations give new insights into the results of recent large-scale numerical simulations of neutrino flavor transformation in supernovae. For example, our calculations reveal the existence of what we term the ``collective precession mode''. Our analyses suggest that neutrinos which travel on intersecting trajectories subject to destructive quantum interference nevertheless can be in this mode. This mode can result in sharp transitions in the final energy-dependent neutrino survival probabilities across all trajectories, a feature seen in the numerical simulations. Moreover, this transition is qualitatively different for the normal and inverted neutrino mass hierarchies. Exploiting this difference, the neutrino signals from a future galactic supernova can potentially be used to determine the actual neutrino mass hierarchy.

We calculate coherent neutrino and antineutrino flavor transformation in the supernova environment, for the first time including a self-consistent treatment of forward scattering-induced coupling and entanglement of intersecting neutrino/antineutrino trajectories. For the atmospheric neutrino mass-squared difference we find that in the normal (inverted) mass hierarchy the more tangentially-propagating (radially-propagating) neutrinos and antineutrinos can initiate collective, simultaneous medium-enhanced flavor conversion of these particles across broad ranges of energy and propagation direction. Accompanying alterations in neutrino/antineutrino energy spectra and/or fluxes could affect supernova nucleosynthesis and the expected neutrino signal.

We present results of large-scale numerical simulations of the evolution of neutrino and antineutrino flavors in the region above the late-time post-supernova-explosion proto-neutron star. Our calculations are the first to allow explicit flavor evolution histories on different neutrino trajectories and to self-consistently couple flavor development on these trajectories through forward scattering-induced quantum entanglement. Employing the atmospheric-scale neutrino mass-squared difference and values of theta_13 allowed by current bounds, we find transformation of neutrino and antineutrino flavors over broad ranges of energy and luminosity in roughly the ``bi-polar'' collective mode. We find that this large-scale flavor conversion, largely driven by the flavor off-diagonal neutrino-neutrino forward scattering potential, sets in much closer to the proto-neutron star than simple estimates based on flavor-diagonal potentials and Mikeheyev-Smirnov-Wolfenstein evolution would indicate. In turn, this suggests that models of r-process nucleosynthesis sited in the neutrino-driven wind could be affected substantially by active-active neutrino flavor mixing, even with the small measured neutrino mass-squared differences.