We calculate and discuss the implications of neutrino self-interactions for the physics of weak decoupling and big bang nucleosynthesis (BBN) in the early universe. In such neutrino-sector extensions of the standard model, neutrinos may not free-stream, yet can stay thermally coupled to one another. Nevertheless, the neutrinos exchange energy and entropy with the photon, electron-positron, and baryon component of the early universe only through the ordinary weak interaction. We examine the effects of neutrino self-interaction for the primordial helium and deuterium abundances and $N_{\rm eff}$, a measure of relativistic energy density at photon decoupling. These quantities are determined in, or may be influenced by, the physics in the weak decoupling epoch. Self interacting neutrinos have been invoked to address a number of anomalies, including as a possible means of ameliorating tension in the Hubble parameter. Our calculations show that surprisingly subtle changes in BBN accompany some of these neutrino self-interaction schemes. Such minute signals require high-precision measurements, making deuterium the best abundance for BBN constraints in the models explored here.

Electron-positron annihilation largely occurs in local thermal and chemical equilibrium after the neutrinos fall out of thermal equilibrium and during the Big Bang Nucleosynthesis (BBN) epoch. The effects of this process are evident in BBN yields as well as the relativistic degrees of freedom. We self-consistently calculate the collision integral for electron-positron creation and annihilation using the Klein-Nishina amplitude and appropriate statistical factors for Fermi-blocking and Bose-enhancement. Our calculations suggest that this annihilation freezes out when the photon-electron-positron-baryon plasma temperature is approximately 16 keV, after which its rate drops below the Hubble rate. In the temperature regime near 16 keV, we break the assumption of chemical equilibrium between the electrons, positrons, and photons to independently calculate the evolution of the chemical potentials of the electrons and positrons while computing the associated collision integrals at every time step. We find that the electron and positron chemical potentials deviate from the case with chemical equilibrium. While our results do not affect the interpretation of precision cosmological measurements in elucidating the standard cosmological model, these out of equilibrium effects may be important for testing physics beyond the standard model.

We provide an overview of the science case, instrument configuration and project plan for the next-generation ground-based cosmic microwave background experiment CMB-S4, for consideration by the 2020 Decadal Survey.

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 present the science case, reference design, and project plan for the Stage-4 ground-based cosmic microwave background experiment CMB-S4.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

Opening up a new window of millimeter-wave observations that span frequency bands in the range of 30 to 500 GHz, survey half the sky, and are both an order of magnitude deeper (about 0.5 uK-arcmin) and of higher-resolution (about 10 arcseconds) than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. In particular, such a survey would allow for major advances in measuring the distribution of dark matter and gas on small-scales, and yield needed insight on 1.) dark matter particle properties, 2.) the evolution of gas and galaxies, 3.) new light particle species, 4.) the epoch of inflation, and 5.) the census of bodies orbiting in the outer Solar System.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

Motivated by the possible existence of other universes, this paper considers the evolution of massive stars with different values for the fundamental constants. We focus on variations in the triple alpha resonance energy and study its effects on the resulting abundances of $^{12}$C, $^{16}$O, and larger nuclei. In our universe, the $0^{+}$ energy level of carbon supports a resonant nuclear reaction that dominates carbon synthesis in stellar cores and accounts for the observed cosmic abundances. Here we define $\Delta{E}_R$ to be the change in this resonant energy level, and show how different values affect the cosmic abundances of the intermediate alpha elements. Using the state of the art computational package $MESA$, we carry out stellar evolution calculations for massive stars in the range $M_\ast$ = $15-40M_\odot$, and for a wide range of resonance energies. We also include both solar and low metallicity initial conditions. For negative $\Delta{E}_R$ , carbon yields are increased relative to standard stellar models, and such universes remain viable as long as the production of carbon nuclei remains energetically favorable, and stars remain stable, down to $\Delta{E}_R\approx-300$ keV. For positive $\Delta{E}_R$, carbon yields decrease, but significant abundances can be produced for resonance energy increments up to $\Delta{E}_R\approx+500$ keV. Oxygen yields tend to be anti-correlated with those of carbon, and the allowed range in $\Delta{E}_R$ is somewhat smaller. We also present yields for neon, magnesium, and silicon. With updated stellar evolution models and a more comprehensive survey of parameter space, these results indicate that the range of viable universes is larger than suggested by earlier studies.

Motivated by the possibility that the laws of physics could be different in other regions of space-time, we consider nuclear processes in universes where the weak interaction is either stronger or weaker than observed. We focus on the physics of both Big Bang Nucleosynthesis (BBN) and stellar evolution. For sufficiently ineffective weak interactions, neutrons do not decay during BBN, and the baryon-to-photon ratio $\eta$ must be smaller in order for protons to survive without becoming incorporated into larger nuclei. For stronger weak interactions, neutrons decay before the onset of BBN, and the early universe is left with nearly a pure hydrogen composition. We then consider stellar structure and evolution for the different nuclear compositions resulting from BBN, a wide range of weak force strengths, and the full range of stellar masses for a given universe. We delineate the range of this parameter space that supports working stars, along with a determination of the dominant nuclear reactions over the different regimes. Deuterium burning dominates the energy generation in stars when the weak force is sufficiently weak, whereas proton-proton burning into helium-3 dominates for the regime where the weak force is much stronger than in our universe. Although stars in these universes are somewhat different, they have comparable surface temperatures, luminosities, radii, and lifetimes, so that a wide range of such universes remain potentially habitable.

We investigate a class of universes in which the weak interaction is not in operation. We consider how astrophysical processes are altered in the absence of weak forces, including Big Bang Nucleosynthesis (BBN), galaxy formation, molecular cloud assembly, star formation, and stellar evolution. Without weak interactions, neutrons no longer decay, and the universe emerges from its early epochs with a mixture of protons, neutrons, deuterium, and helium. The baryon-to-photon ratio must be smaller than the canonical value in our universe to allow free nucleons to survive the BBN epoch without being incorporated into heavier nuclei. At later times, the free neutrons readily combine with protons to make deuterium in sufficiently dense parts of the interstellar medium, and provide a power source before they are incorporated into stars. Almost all of the neutrons are incorporated into deuterium nuclei before stars are formed. As a result, stellar evolution proceeds primarily through strong interactions, with deuterium first burning into helium, and then helium fusing into carbon. Low-mass deuterium-burning stars can be long-lived, and higher mass stars can synthesize the heavier elements necessary for life. Although somewhat different from our own, such universes remain potentially habitable.

We present calculations showing how electron rest mass influences entropy flow, neutrino decoupling, and Big Bang Nucleosynthesis (BBN) in the early universe. To elucidate this physics and especially the sensitivity of BBN and related epochs to electron mass, we consider a parameter space of rest mass values larger and smaller than the accepted vacuum value. Electromagnetic equilibrium, coupled with the high entropy of the early universe, guarantees that significant numbers of electron-positron pairs are present, and dominate over the number of ionization electrons to temperatures much lower than the vacuum electron rest mass. Scattering between the electrons-positrons and the neutrinos largely controls the flow of entropy from the plasma into the neutrino seas. Moreover, the number density of electron-positron-pair targets can be exponentially sensitive to the effective in-medium electron mass. This entropy flow influences the phasing of scale factor and temperature, the charged current weak-interaction-determined neutron-to-proton ratio, and the spectral distortions in the relic neutrino energy spectra. Our calculations show the sensitivity of the physics of this epoch to three separate effects: finite electron mass, finite-temperature quantum electrodynamic (QED) effects on the plasma equation of state, and Boltzmann neutrino energy transport. The ratio of neutrino to plasma component energy scales manifests in Cosmic Microwave Background (CMB) observables, namely the baryon density and the radiation energy density, along with the primordial helium and deuterium abundances. Our results demonstrate how the treatment of in-medium electron mass (i.e., QED effects) could translate into an important source of uncertainty in extracting neutrino and beyond-standard-model physics limits from future high-precision CMB data.

In both stars and in the early universe, the production of deuterium is the first step on the way to producing heavier nuclei. If the strong force were slightly weaker, deuterium would not be stable, and many authors have noted that nuclesynthesis would be compromised so that helium production could not proceed through standard reaction chains. Motivated by the possibility that other regions of space-time could have different values for the fundamental constants, this paper considers stellar evolution in universes without stable deuterium and argues that such universes can remain habitable. Even in universes with no stellar nucleosynthesis, stars can form and will generate energy through gravitational contraction. We show that such stars can be sufficiently luminous and long-lived to support life. Stars with initial masses that exceed the Chandrasekhar mass cannot be supported by degeneracy pressure and explode at the end of their contraction phase. The resulting explosive nucleosynthesis can provide the universe with some heavy elements. We also explore the possibility that helium can be produced in stellar cores through a triple-nucleon reaction (roughly analogous to the triple-alpha reaction). Next we show that with even trace amounts of heavy elements --- produced through the triple-nucleon process or by explosive nucleosynthesis --- the CNO cycle can operate and allow stars to function. Finally, we consider Big Bang Nucleosynthesis without stable deuterium and find that only trace amounts of helium are produced, with even smaller abundances of other nuclei. With stars evolving through gravitational contraction, explosive nucleosynthesis, the triple-nucleon reaction, and the CNO cycle, universes with no stable deuterium are thus potentially habitable, contrary to many previous claims.

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.

This paper derives an upper limit on the density $\rho_{\scriptstyle\Lambda}$ of dark energy based on the requirement that cosmological structure forms before being frozen out by the eventual acceleration of the universe. By allowing for variations in both the cosmological parameters and the strength of gravity, the resulting constraint is a generalization of previous limits. The specific parameters under consideration include the amplitude $Q$ of the primordial density fluctuations, the Planck mass $M_{\rm pl}$, the baryon-to-photon ratio $\eta$, and the density ratio $\Omega_M/\Omega_b$. In addition to structure formation, we use considerations from stellar structure and Big Bang Nucleosynthesis (BBN) to constrain these quantities. The resulting upper limit on the dimensionless density of dark energy becomes $\rho_{\scriptstyle\Lambda}/M_{\rm pl}^4<10^{-90}$, which is $\sim30$ orders of magnitude larger than the value in our universe $\rho_{\scriptstyle\Lambda}/M_{\rm pl}^4\sim10^{-120}$. This new limit is much less restrictive than previous constraints because additional parameters are allowed to vary. With these generalizations, a much wider range of universes can develop cosmic structure and support observers. To constrain the constituent parameters, new BBN calculations are carried out in the regime where $\eta$ and $G=M_{\rm pl}^{-2}$ are much larger than in our universe. If the BBN epoch were to process all of the protons into heavier elements, no hydrogen would be left behind to make water, and the universe would not be viable. However, our results show that some hydrogen is always left over, even under conditions of extremely large $\eta$ and $G$, so that a wide range of alternate universes are potentially habitable.

Motivated by the possible existence of other universes, with different values for the fundamental constants, this paper considers stellar models in universes where $^8$Be is stable. Many previous authors have noted that stars in our universe would have difficulty producing carbon and other heavy elements in the absence of the well-known $^{12}$C resonance at 7.6 MeV. This resonance is necessary because $^8$Be is unstable in our universe, so that carbon must be produced via the triple alpha reaction to achieve the requisite abundance. Although a moderate change in the energy of the resonance (200 -- 300 keV) will indeed affect carbon production, an even smaller change in the binding energy of beryllium ($\sim100$ keV) would allow $^8$Be to be stable. A stable isotope with $A=8$ would obviate the need for the triple alpha process in general, and the $^{12}$C resonance in particular, for carbon production. This paper explores the possibility that $^8$Be can be stable in other universes. Simple nuclear considerations indicate that bound states can be realized, with binding energy $\sim0.1-1$ MeV, if the fundamental constants vary by a $\sim {\rm few}-10$ percent. In such cases, $^8$Be can be synthesized through helium burning, and $^{12}$C can be produced later through nuclear burning of beryllium. This paper focuses on stellar models that burn helium into beryllium; once the universe in question has a supply of stable beryllium, carbon production can take place during subsequent evolution in the same star or in later stellar generations. Using both a semi-analytic stellar structure model as well as a state-of-the-art stellar evolution code, we find that viable stellar configurations that produce beryllium exist over a wide range of parameter space. Finally, we demonstrate that carbon can be produced during later evolutionary stages.

This book lays out the scientific goals to be addressed by the next-generation ground-based cosmic microwave background experiment, CMB-S4, envisioned to consist of dedicated telescopes at the South Pole, the high Chilean Atacama plateau and possibly a northern hemisphere site, all equipped with new superconducting cameras. CMB-S4 will dramatically advance cosmological studies by crossing critical thresholds in the search for the B-mode polarization signature of primordial gravitational waves, in the determination of the number and masses of the neutrinos, in the search for evidence of new light relics, in constraining the nature of dark energy, and in testing general relativity on large scales.

The weak interaction charged current processes ($\nu_e+n\leftrightarrow p+e^-$, $\bar\nu_e +p\leftrightarrow n+e^+$, $n\leftrightarrow p+e^-+\bar\nu_e$) interconvert neutrons and protons in the early universe and have significant influence on Big Bang Nucleosynthesis (BBN) light-element abundance yields, particulary that for $^{4}{\rm He}$. We demonstrate that the influence of these processes is still significant even when they operate well below temperatures $T\sim0.7\,{\rm MeV}$ usually invoked for "weak freeze-out," and in fact down nearly into the alpha-particle formation epoch ($T \approx 0.1\,{\rm MeV}$). This physics is correctly captured in commonly used BBN codes, though this late-time, low-temperature persistent effect of the isospin-changing weak processes, and the sensitivity of the associated rates to lepton energy distribution functions and blocking factors are not widely appreciated. We quantify this late-time influence by analyzing weak interaction rate dependence on the neutron lifetime, lepton energy distribution functions, entropy, the proton-neutron mass difference, and Hubble expansion rate. The effects we point out here render BBN a keen probe of any beyond-standard-model physics that alters lepton number/energy distributions, even subtly, in epochs of the early universe all the way down to near $T=100\,{\rm keV}$.

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 examine the physics of the early universe when Majorana neutrinos (electron neutrino, muon neutrino, tau neutrino) possess transition magnetic moments. These extra couplings beyond the usual weak interaction couplings alter the way neutrinos decouple from the plasma of electrons/positrons and photons. We calculate how transition magnetic moment couplings modify neutrino decoupling temperatures, and then use a full weak, strong, and electromagnetic reaction network to compute corresponding changes in Big Bang Nucleosynthesis abundance yields. We find that light element abundances and other cosmological parameters are sensitive to magnetic couplings on the order of 10^{-10} Bohr magnetons. Given the recent analysis of sub-MeV Borexino data which constrains Majorana moments to the order of 10^{-11} Bohr magnetons or less, we find that changes in cosmological parameters from magnetic contributions to neutrino decoupling temperatures are below the level of upcoming precision observations.

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.