There exist a range of exciting scientific opportunities for Big Bang Nucleosynthesis (BBN) in the coming decade. BBN, a key particle astrophysics "tool" for decades, is poised to take on new capabilities to probe beyond standard model (BSM) physics. This development is being driven by experimental determination of neutrino properties, new nuclear reaction experiments, advancing supercomputing/simulation capabilities, the prospect of high-precision next-generation cosmic microwave background (CMB) observations, and the advent of 30m class telescopes.

We perform a self-consistent calculation of collective neutrino oscillations including the effect of back scattered neutrinos (Halo effect) in the `single-angle' approximation, within a spherically symmetric supernova model. We find that due to the Halo effect the onset of flavor transformations is pushed to smaller radii, by a few kilometers. The celebrated phenomenon of the spectral split is found to be robust under the present inclusion of the Halo effect.

We investigate the stability, with respect to spatial inhomogeneity, of a two-dimensional dense neutrino gas. The system exhibits growth of seed inhomogeneity due to nonlinear coherent neutrino self-interactions. In the absence of incoherent collisional effects, we observe a dependence of this instability growth rate on the neutrino mass spectrum: the normal neutrino mass hierarchy exhibits spatial instability over a larger range of neutrino number density compared to that of the inverted case. We further consider the effect of elastic incoherent collisions of the neutrinos with a static background of heavy, nucleon-like scatterers. At small scales, the growth of flavor instability can be suppressed by collisions. At large length scales we find, perhaps surprisingly, that for inverted neutrino mass hierarchy incoherent collisions fail to suppress flavor instabilities, independent of the coupling strength.

We calculate Boltzmann neutrino energy transport with self-consistently coupled nuclear reactions through the weak-decoupling-nucleosynthesis epoch in an early universe with significant lepton numbers. We find that the presence of lepton asymmetry enhances processes which give rise to nonthermal neutrino spectral distortions. Our results reveal how asymmetries in energy and entropy density uniquely evolve for different transport processes and neutrino flavors. The enhanced distortions in the neutrino spectra alter the expected big bang nucleosynthesis light element abundance yields relative to those in the standard Fermi-Dirac neutrino distribution cases. These yields, sensitive to the shapes of the neutrino energy spectra, are also sensitive to the phasing of the growth of distortions and entropy flow with time/scale factor. We analyze these issues and speculate on new sensitivity limits of deuterium and helium to lepton number.

We investigate neutrino flavor transformation in the early universe in the presence of a lepton asymmetry, focusing on a two-flavor system with 1 - 3 mixing parameters. We identify five distinct regimes that emerge in an approximate treatment neglecting collisions as the initial lepton asymmetry at high temperature is varied from values comparable to current constraints on the lepton number down to values at which the neutrino-neutrino forward-scattering potential is negligible. The characteristic phenomena occurring in these regimes are (1) large synchronized oscillations, (2) minimal flavor transformation, (3) asymmetric (neutrino- or antineutrino-only) MSW, (4) partial MSW, and (5) symmetric MSW. We examine our numerical results in the framework of adiabaticity, and we illustrate how they are modified by collisional damping. Finally, we point out the existence of matter-neutrino resonances in the early universe and show that they suffer from non-adiabaticity.

We calculate the evolution of the early universe through the epochs of weak decoupling, weak freeze-out and big bang nucleosynthesis (BBN) by simultaneously coupling a full strong, electromagnetic, and weak nuclear reaction network with a multi-energy group Boltzmann neutrino energy transport scheme. The modular structure of our code provides the ability to dissect the relative contributions of each process responsible for evolving the dynamics of the early universe in the absence of neutrino flavor oscillations. Such an approach allows a detailed accounting of the evolution of the $\nu_e$, $\bar\nu_e$, $\nu_\mu$, $\bar\nu_\mu$, $\nu_\tau$, $\bar\nu_\tau$ energy distribution functions alongside and self-consistently with the nuclear reactions and entropy/heat generation and flow between the neutrino and photon/electron/positron/baryon plasma components. This calculation reveals nonlinear feedback in the time evolution of neutrino distribution functions and plasma thermodynamic conditions (e.g., electron-positron pair densities), with implications for: the phasing between scale factor and plasma temperature; the neutron-to-proton ratio; light-element abundance histories; and the cosmological parameter \neff. We find that our approach of following the time development of neutrino spectral distortions and concomitant entropy production and extraction from the plasma results in changes in the computed value of the BBN deuterium yield. For example, for particular implementations of quantum corrections in plasma thermodynamics, our calculations show a $0.4\%$ increase in deuterium. These changes are potentially significant in the context of anticipated improvements in observational and nuclear physics uncertainties.

We discuss how small neutrino rest masses can increase the expansion rate near the photon decoupling epoch in the early universe, causing an earlier, higher temperature freeze-out for ionization equilibrium compared to the massless neutrino case. This yields a larger free-electron fraction. A larger ratio of the sound horizon to the photon diffusion length follows, implying a smaller inferred Neff. This neutrino-mass/recombination effect depends strongly on the neutrino rest masses. Though below current sensitivity, this effect could be probed by next-generation cosmic microwave background experiments, giving an observational handle of neutrino mass physics.

We show that a self-consistent and coupled treatment of the weak decoupling, big bang nucleosynthesis, and photon decoupling epochs can be used to provide new insights and constraints on neutrino sector physics from high-precision measurements of light element abundances and cosmic microwave background observables. Implications of beyond-standard-model physics in cosmology, especially within the neutrino sector, are assessed by comparing predictions against five observables: the baryon energy density, helium abundance, deuterium abundance, effective number of neutrinos, and sum of the light neutrino mass eigenstates. We give examples for constraints on dark radiation, neutrino rest mass, lepton numbers, and scenarios for light and heavy sterile neutrinos.

Recent activity in solving the 'lithium problem' in big bang nucleosynthesis has focused on the role that putative resonances may play in resonance-enhanced destruction of 7Li. Particular attention has been paid to the reactions involving the 9B compound nuclear system, d+7Be->9B. These reactions are analyzed via the multichannel, two-body unitary R-matrix method using code (EDA) developed by Hale and collaborators. We employ much of the known elastic and reaction data, in a four-channel treatment. The data include elastic 3He+6Li differential cross sections from 0.7 to 2.0 MeV, integrated reaction cross sections for energies from 0.7 to 5.0 MeV for 6Li(3He,p)8Be* and from 0.4 to 5.0 MeV for the 6Li(3He,d)7Be reaction. Capture data have been added to an earlier analysis with integrated cross section measurements from 0.7 to 0.825 MeV for 6Li(3He,gamma)9B. The resulting resonance parameters are compared with tabulated values, and previously unidentified resonances are noted. Our results show that there are no near d+7Be threshold resonances with widths that are 10's of keV and reduce the likelihood that a resonance-enhanced mass-7 destruction mechanism, as suggested in recently published work, can explain the 7Li problem.

We show that a hybrid (nuclear + quark matter) star can have a mass-radius relationship very similar to that predicted for a star made of purely nucleonic matter. We show this for a generic parameterization of the quark matter equation of state, and also for an MIT bag model, each including a phenomenological correction based on gluonic corrections to the equation of state. We obtain hybrid stars as heavy as 2 M_solar for reasonable values of the bag model parameters. For nuclear matter, we use the equation of state calculated by Akmal, Pandharipande, and Ravenhall using many-body techniques. Both mixed and homogeneous phases of nuclear and quark matter are considered.