The multi-messenger astrophysics of compact objects presents a vast range of environments where neutrino flavor transformation may occur and may be important for nucleosynthesis, dynamics, and a detected neutrino signal. Development of efficient techniques for surveying flavor evolution solution spaces in these diverse environments, which augment and complement existing sophisticated computational tools, could leverage progress in this field. To this end we continue our exploration of statistical data assimilation (SDA) to identify solutions to a small-scale model of neutrino flavor transformation. SDA is a machine learning (ML) formula wherein a dynamical model is assumed to generate any measured quantities. Specifically, we use an optimization formulation of SDA wherein a cost function is extremized via the variational method. Regions of state space in which the extremization identifies the global minimum of the cost function will correspond to parameter regimes in which a model solution can exist. Our example study seeks to infer the flavor transformation histories of two mono-energetic neutrino beams coherently interacting with each other and with a matter background. We require that the solution be consistent with measured neutrino flavor fluxes at the point of detection, and with constraints placed upon the flavor content at various locations along their trajectories, such as the point of emission, and the locations of the Mikheyev-Smirnov-Wolfenstein (MSW) resonances. We show how the procedure efficiently identifies solution regimes and rules out regimes where solutions are infeasible. Overall, results intimate the promise of this \lq\lq variational annealing\rq\rq\ methodology to efficiently probe an array of fundamental questions that traditional numerical simulation codes render difficult to access.

We investigate the importance of going beyond the mean-field approximation in the dynamics of collective neutrino oscillations. To expand our understanding of the coherent neutrino oscillation problem, we apply concepts from many-body physics and quantum information theory. Specifically, we use measures of nontrivial correlations (otherwise known as "entanglement") between the constituent neutrinos of the many-body system, such as the entanglement entropy and the Bloch vector of the reduced density matrix. The relevance of going beyond the mean field is demonstrated by comparisons between the evolution of the neutrino state in the many-body picture vs the mean-field limit, for different initial conditions.

We demonstrate a method to systematically obtain eigenvalues and eigenstates of a many-body Hamiltonian describing collective neutrino oscillations. The method is derived from the Richardson-Gaudin framework, which involves casting the eigenproblem as a set of coupled nonlinear "Bethe Ansatz equations", the solutions of which can then be used to parametrize the eigenvalues and eigenvectors. The specific approach outlined in this paper consists of defining auxiliary variables that are related to the Bethe-Ansatz parameters, thereby transforming the Bethe-Ansatz equations into a different set of equations that are numerically better behaved and more tractable. We show that it is possible to express not only the eigenvalues, but also the eigenstates, directly in terms of these auxiliary variables without involving the Bethe Ansatz parameters themselves. In this paper, we limit ourselves to a two-flavor, single-angle neutrino system.

We consider systems where dynamical variables are the generators of the SU(2) group. A subset of these Hamiltonians is exactly solvable using the Bethe ansatz techniques. We show that Bethe ansatz equations are equivalent to polynomial relationships between the operator invariants, or equivalently, between eigenvalues of those invariants.

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.

We examine the flavor evolution of neutrinos emitted from the disk-like remnant (hereafter called \lq\lq neutrino disk\rq\rq) of a binary neutron star (BNS) merger. We specifically follow the neutrinos emitted from the center of the disk, along the polar axis perpendicular to the equatorial plane. We carried out two-flavor simulations using a variety of different possible initial neutrino luminosities and energy spectra, and for comparison, three-flavor simulations in specific cases. In all simulations, the normal neutrino mass hierarchy was used. The flavor evolution was found to be highly dependent on the initial neutrino luminosities and energy spectra; in particular, we found two broad classes of results depending on the sign of the initial net electron neutrino lepton number (i.e., the number of neutrinos minus the number of antineutrinos). In the antineutrino dominated case, we found that the Matter-Neutrino Resonance (MNR) effect dominates, consistent with previous results, whereas in the neutrino dominated case, a bipolar spectral swap develops. The neutrino dominated conditions required for this latter result have been realized, e.g, in a BNS merger simulation that employs the \lq\lq DD2\rq\rq\ equation of state for neutron star matter[Phys. Rev. D 93, 044019 (2016)]. For this case, in addition to the swap at low energies, a collective Mikheyev-Smirnov-Wolfenstein (MSW) mechanism generates a high-energy electron neutrino tail. The enhanced population of high-energy electron neutrinos in this scenario could have implications for the prospects of $r$-process nucleosynthesis in the material ejected outside the plane of the neutrino disk.

We present neutrino bulb model simulations of Majorana neutrino coherent spin transformation (i.e., neutrino-antineutrino transformation), coupled to neutrino flavor evolution, for conditions corresponding to the neutronization burst epoch of an Oxygen-Neon-Magnesium (O-Ne-Mg) core collapse supernova. Significant neutrino spin transformation in, for example, the neutronization burst, could alter the fluences of neutrinos and antineutrinos in a way which is potentially detectable for a Galactic core collapse supernova. Our calculations for the first time incorporate geometric dilution in the spin evolution of the neutrinos and combine two-flavor and three-flavor evolution with spin mixing physics. We find that significant spin transformations can occur, but only with a large neutrino luminosity and an electron fraction ($Y_e$) profile which facilitates adiabatic conditions for the spin-channel resonance. Using our adopted parameters of neutrino energy spectra, luminosity, density and $Y_e$ profiles, our calculations require an unrealistically large neutrino rest mass to sustain the spin transformation. It is an open question whether examining different density profiles or incorporating other sources of nonlinear feedback, such as $Y_e$ feedback, could mitigate this need. We find that spin transformations are not sensitive to the flavor structure of neutrinos, i.e., the spin transformations occur regardless of whether we simulate two- or three-flavor transformations. In the two-flavor case, spin transformations were insensitive to the choice of solar or atmospheric mass-squared splitting as well as the choice of the Majorana phase. Importantly, our three-flavor simulations, as well as our two-flavor simulations done with the atmospheric mass-squared splitting, show that the inclusion of spin degrees of freedom can significantly and qualitatively alter neutrino flavor evolution.

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

We present a model where sterile neutrinos with rest masses in the range ~ keV to ~ MeV can be the dark matter and be consistent with all laboratory, cosmological, large-scale structure, as well as x-ray constraints. These sterile neutrinos are assumed to freeze out of thermal and chemical equilibrium with matter and radiation in the very early Universe, prior to an epoch of prodigious entropy generation ("dilution") from out-of-equilibrium decay of heavy particles. In this work, we consider heavy, entropy-producing particles in the ~ TeV to ~ EeV rest-mass range, possibly associated with new physics at high-energy scales. The process of dilution can give the sterile neutrinos the appropriate relic densities, but it also alters their energy spectra so that they could act like cold dark matter, despite relatively low rest masses as compared to conventional dark matter candidates. Moreover, since the model does not rely on active-sterile mixing for producing the relic density, the mixing angles can be small enough to evade current x-ray or lifetime constraints. Nevertheless, we discuss how future x-ray observations, future lepton number constraints, and future observations and sophisticated simulations of large-scale structure could, in conjunction, provide evidence for this model and/or constrain and probe its parameters.

We show that a particular class of postrecombination phase transitions in the vacuum can lead to localized overdense regions on relatively small scales, roughly 10^6 to 10^10 M_sun, potentially interesting for the origin of large black hole seeds and for dwarf galaxy evolution. Our study suggests that this mechanism could operate over a range of conditions which are consistent with current cosmological and laboratory bounds. One byproduct of phase transition bubble-wall decay may be extra radiation energy density. This could provide an avenue for constraint, but it could also help reconcile the discordant values of the present Hubble parameter (H_0) and sigma_8 obtained by cosmic microwave background (CMB) fits and direct observational estimates. We also suggest ways in which future probes, including CMB considerations (e.g., early dark energy limits), 21-cm observations, and gravitational radiation limits, could provide more stringent constraints on this mechanism and the sub-eV scale beyond-standard-model physics, perhaps in the neutrino sector, on which it could be based. Late phase transitions associated with sterile neutrino mass and mixing may provide a way to reconcile cosmological limits and laboratory data, should a future disagreement arise.