Recent theoretical work indicates that the neutrino radiation in core-collapse supernovae may be susceptible to flavor instabilities that set in far behind the shock, grow extremely rapidly, and have the potential to profoundly affect supernova dynamics and composition. Here we analyze the nonlinear collective oscillations that are prefigured by these instabilities. We demonstrate that a zero-crossing in $n_{\nu_e} - n_{\bar{\nu}_e}$ as a function of propagation angle is not sufficient to generate instability. Our analysis accounts for this fact and allows us to formulate complementary criteria. Using Fornax simulation data, we show that fast collective oscillations qualitatively depend on how forward-peaked the neutrino angular distributions are.

An extensive, growing body of work has been penned on cosmologies that include one or more sterile neutrinos. Early entries in the literature formulated a Boltzmann-like equation describing sterile-neutrino production in a way that bypasses the numerical tracking of high-frequency complex phases, and meticulous quantum-kinetic analyses shortly thereafter put the formula on firmer ground. A new and more direct derivation is given here, showing that the equation follows almost immediately from a quantum relaxation-time approximation and an expansion in the mixing angle. Besides reproducing the desired result, the relaxation ansatz captures to a high degree of accuracy the interlaced dynamics of oscillations, decoherence, and plasma repopulation. Successes and limitations of the semiclassical equation are illustrated numerically and are shown to reflect the accuracy of the approximations employed in the derivation. The inclusion of interactions among the sterile neutrinos is also briefly addressed.

For active-sterile mixing to be responsible for the full relic abundance of dark matter additional new physics is needed beyond the keV-scale sterile neutrino itself. The extra ingredient we consider here is the presence of self-interactions among the sterile neutrinos. We examine whether active-to-sterile conversion is amplified enough in this scenario that the observed abundance of dark matter can be obtained with a subconstraint mixing angle. This turns out never to be the case in the region we explore: either self-interactions have too small an impact and cannot escape bounds on the mass and mixing angle, or they have too great an impact and cause dark matter to be overproduced. The sharp transition from marginal to excessive effectiveness occurs because a resonance criterion is met in the effective in-medium mixing angle. Once the system goes resonant the game is as good as over, as nonlinearity in the Boltzmann equation leads to runaway production of sterile neutrinos, beginning at a plasma temperature of a few hundred MeV and typically ending at a few tens of MeV. The scenario is therefore ruled out largely by its own dynamics. In this study we focus exclusively on mediators heavier than $\sim 1$ GeV; future work will extend the analysis to lighter mediators, allowing for contact to be made with the kinds of scenarios motivated by issues of small-scale structure.

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.

The thermodynamic arrow-of-time problem is thought to be resolved by the observation that our universe initially was---and still is---far from equilibrium. The psychological arrow-of-time problem is often attributed the same resolution, but the connection has not been thoroughly established. I argue that a compelling explanation of the psychological arrow requires an understanding of the physical conditions necessary for life to emerge from a prebiotic environment. A simple calculation illustrates how life-sustaining energy fluxes from the Sun and the Earth's interior bias the development of life in the direction of the arrow inherited from cosmic initial conditions.

We identify in the flavor transformation of astrophysical neutrinos a new class of phenomena, a common outcome of which is the suppression of flavor conversion. Appealing to the equivalence between a bipolar neutrino system and a gyroscopic pendulum, we find that these phenomena have rather striking interpretations in the mechanical picture: in one instance, the gyroscopic pendulum initially precesses in one direction, then comes to a halt and begins to precess in the opposite direction---a counterintuitive behavior that we analogize to the motion of a toy known as a rattleback. We analyze these behaviors in the early universe, wherein a chance connection to sterile neutrino dark matter emerges, and we briefly suggest how they might manifest in compact-object environments.

We present our determination of the nuclear supermassive black hole (SMBH) mass function for spiral galaxies in the Local Universe, established from a volume-limited sample consisting of a statistically complete collection of the brightest spiral galaxies in the Southern Hemisphere. Our SMBH mass function agrees well at the high-mass end with previous values given in the literature. At the low-mass end, inconsistencies exist in previous works that still need to be resolved, but our work is more in line with expectations based on modeling of SMBH evolution. This low-mass end of the spectrum is critical to our understanding of the mass function and evolution of SMBHs since the epoch of maximum quasar activity. A luminosity distance $\leq$ 25.4 $Mpc$ and an absolute B-band magnitude $\leq$ -19.12 define the sample. These limits define a sample of 140 spiral galaxies, with 128 measurable pitch angles to establish the pitch angle distribution for this sample. This pitch angle distribution function may be useful in the study of the morphology of late-type galaxies. We then use an established relationship between the pitch angle and the mass of the central SMBH in a host galaxy in order to estimate the mass of the 128 respective SMBHs in this sample. This result effectively gives us the distribution of mass for SMBHs residing in spiral galaxies over a lookback time $\leq$ 82.1 $h_{67.77}^{-1}$ $Myr$ and contained within a comoving volume of 3.37 $\times$ $10^4$ $h_{67.77}^{-3}$ $Mpc^3$. We estimate the density of SMBHs residing in spiral galaxies in the Local Universe is $5.54_{-2.73}^{+6.55}$ $\times$ $10^4$ $h_{67.77}^3$ $M_{\odot}$ $Mpc^{-3}$. Thus, our derived cosmological SMBH mass density for spiral galaxies is $\Omega_{BH} = 4.35_{-2.15}^{+5.14}$ $\times$ $10^{-7}$ $h_{67.77}$.

We assess the utility of an optimization-based data assimilation (D.A.) technique for treating the problem of nonlinear neutrino flavor transformation in core collapse supernovae. D.A. uses measurements obtained from a physical system to estimate the state variable evolution and parameter values of the associated model. Formulated as an optimization procedure, D.A. can offer an integration-blind approach to predicting model evolution, which offers an advantage for models that thwart solution via traditional numerical integration techniques. Further, D.A. performs most optimally for models whose equations of motion are nonlinearly coupled. In this exploratory work, we consider a simple steady-state model with two mono-energetic neutrino beams coherently interacting with each other and a background medium. As this model can be solved via numerical integration, we have an independent consistency check for D.A. solutions. We find that the procedure can capture key features of flavor evolution over the entire trajectory, even given measurements of neutrino flavor only at the endpoint, and with an assumed known initial flavor distribution. Further, the procedure permits an examination of the sensitivity of flavor evolution to estimates of unknown model parameters, locates degeneracies in parameter space, and can identify the specific measurements required to break those degeneracies.

Neutrinos propagating in dense astrophysical environments sustain nonlinear refractive effects due to neutrino-neutrino forward scattering. We study geometric phases in neutrino oscillations that arise out of cyclic evolution of the potential generated by these forward-scattering processes. We perform several calculations, exact and perturbative, that illustrate the robustness of such phases, and of geometric effects more broadly, in the flavor evolution of neutrinos. The scenarios we consider are highly idealized in order to make them analytically tractable, but they suggest the possible presence of complicated geometric effects in realistic astrophysical settings. We also point out that in the limit of extremely high neutrino densities, the nonlinear potential in three flavors naturally gives rise to non-Abelian geometric phases. This paper is intended to be accessible to neutrino experts and non-specialists alike.

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.