We consider the many-body system of neutrinos interacting with each other through neutral current weak force. Emerging many-body effects in such a system could play important roles in some astrophysical sites such as the core collapse supernovae. In the literature this many-body system is usually treated within the mean field approximation which is an effective one-body description based on omitting entangled neutrino states. In this paper, we consider the original many-body system in an effective two flavor mixing scenario under the single angle approximation and present a solution without using the mean field approximation. Our solution is formulated around a special class of many-body eigenstates which do not undergo any level crossings as the neutrino self interaction rate decreases while the neutrinos radiate from the supernova. In particular, an initial state which consists of electron neutrinos and antineutrinos of an orthogonal flavor can be entirely decomposed in terms of those eigenstates. Assuming that the conditions are perfectly adiabatic so that the evolution of these eigenstates follow their variation with the interaction rate, we show that this initial state develops a spectral split at exactly the same energy predicted by the mean field formulation.

We study the effects of collective neutrino oscillations on $\nu p$ process nucleosynthesis in proton-rich neutrino-driven winds by including both the multi-angle $3\times3$ flavor mixing and the nucleosynthesis network calculation. The number flux of energetic electron antineutrinos is raised by collective neutrino oscillations in a $1$D supernova model for $40 M_{\odot}$ progenitor. When the gas temperature decreases down to $\sim2-3\times10^{9}$ K, the increased flux of electron antineutrinos promotes $\nu p$ process more actively, resulting in the enhancement of $p$-nuclei. In the early phase of neutrino-driven wind, blowing at $0.6$ s after core bounce, oscillation effects are prominent in inverted mass hierarchy and $p$-nuclei are synthesized up to $^{106}\mathrm{Cd}$ and $^{108}\mathrm{Cd}$. On the other hand, in the later wind trajectory at $1.1$ s after core bounce, abundances of $p$-nuclei are increased remarkably by $\sim10-10^{4}$ times in normal mass hierarchy and even reaching heavier $p$-nuclei such as $^{124}\mathrm{Xe}$, $^{126}\mathrm{Xe}$ and $^{130}\mathrm{Ba}$. The averaged overproduction factor of $p$-nuclei is dominated by the later wind trajectories. Our results demonstrate that collective neutrino oscillations can strongly influence $\nu p$ process, which indicates that they should be included in the network calculations in order to obtain precise abundances of $p$-nuclei. The conclusions of this paper depend on the difference of initial neutrino parameters between electron and non-electron antineutrino flavors which is large in our case. Further systematic studies on input neutrino physics and wind trajectories are necessary to draw a robust conclusion. However, this finding would help understand the origin of solar-system isotopic abundances of $p$-nuclei such as $^{92,94}\mathrm{Mo}$ and $^{96,98}\mathrm{Ru}$.

We show that the spectral split of a neutrino ensemble which initially consists of electron type neutrinos, is analogous to the BCS-BEC crossover already observed in ultra cold atomic gas experiments. Such a neutrino ensemble mimics the deleptonization burst of a core collapse supernova. Although these two phenomena belong to very different domains of physics, the propagation of neutrinos from highly interacting inner regions of the supernova to the vacuum is reminiscent of the evolution of Cooper pairs between weak and strong interaction regimes during the crossover. The Hamiltonians and the corresponding many-body states undergo very similar transformations if one replaces the pair quasispin of the latter with the neutrino isospin of the former.

In a core collapse supernova, collective oscillations of neutrinos emitted by the proto-neutron star significantly modifies the partition of energy between different flavors in a way which is potentially important for the nucleosynthesis of heavy elements at the outer layers. As the outcome of collective oscillations depends nonlinearly on the chosen initial conditions proto-neutron star, their systematical study across the proto-neutron star model space may provide answers to some key questions.

We consider collective oscillations of neutrinos, which are emergent nonlinear flavor evolution phenomena instigated by neutrino-neutrino interactions in astrophysical environments with sufficiently high neutrino densities. We investigate the symmetries of the problem in the full three flavor mixing scheme and in the exact many-body formulation by including the effects of CP violation and neutrino magnetic moment. We show that, similar to the two flavor scheme, several dynamical symmetries exist for three flavors in the single-angle approximation if the net electron background in the environment and the effects of the neutrino magnetic moment are negligible. Moreover, we show that these dynamical symmetries are present even when the CP symmetry is violated in neutrino oscillations. We explicitly write down the constants of motion through which these dynamical symmetries manifest themselves in terms of the generators of the SU(3) flavor transformations. We also show that the effects due to the CP-violating Dirac phase factor out of the many-body evolution operator and evolve independently of nonlinear flavor transformations if neutrino electromagnetic interactions are ignored. In the presence of a strong magnetic field, CP-violating effects can still be considered independently provided that an effective definition for neutrino magnetic moment is used.

The effective many-body Hamiltonian which describes vacuum oscillations and self interactions of neutrinos in a two flavor mixing scheme under the single angle approximation has the same dynamical symmetries as the well known BCS pairing Hamiltonian. These dynamical symmetries manifest themselves in terms of a set of constants of motion and can be useful in formulating the collective oscillation modes in an intuitive way. In particular, we show that a neutrino spectral split can be simply viewed as an avoided level crossing between the eigenstates of a mean field Hamiltonian which includes a Lagrange multiplier in order to fix the value of an exact many-body constant of motion. We show that the same dynamical symmetries also exist in the three neutrino mixing scheme by explicitly writing down the corresponding constants of motion.

Using the primordial helium abundance, an upper limit to the magnetic moments for Dirac neutrinos had been provided by imposing restrictions on the number of the additional helicity states. Considering non-thermal photons produced in the decay of the heavy sterile mass eigenstates due to the neutrino magnetic moment, we explore the constraints imposed by the observed abundances of all the light elements produced during the Big Bang nucleosynthesis.

Following a recent suggestion of axion cooling of photons between the nucleosynthesis and recombination epochs in the Early Universe, we investigate a hybrid model with both axions and relic supersymmetric particles. In this model we demonstrate that the 7Li abundance can be consistent with observations without destroying the important concordance of deuterium abundance.

We consider the flavor evolution of a dense neutrino gas by taking into account both vacuum oscillations and self interactions of neutrinos. We examine the system from a many-body perspective as well as from the point of view of an effective one-body description formulated in terms of the neutrino polarization vectors. We show that, in the single angle approximation, both the many-body picture and the effective one-particle picture possess several constants of motion. We write down these constants of motion explicitly in terms of the neutrino isospin operators for the many-body case and in terms of the polarization vectors for the effective one-body case. The existence of these constants of motion is a direct consequence of the fact that the collective neutrino oscillation Hamiltonian belongs to the class of Gaudin Hamiltonians. This class of Hamiltonians also includes the (reduced) BCS pairing Hamiltonian describing superconductivity. We point out the similarity between the collective neutrino oscillation Hamiltonian and the BCS pairing Hamiltonian. The constants of motion manifest the exact solvability of the system. Borrowing the well established techniques of calculating the exact BCS spectrum, we present exact eigenstates and eigenvalues of both the many-body and the effective one-particle Hamiltonians describing the collective neutrino oscillations. For the effective one-body case, we show that spectral splits of neutrinos can be understood in terms of the adiabatic evolution of some quasi-particle degrees of freedom from a high density region where they coincide with flavor eigenstates to the vacuum where they coincide with mass eigenstates. We write down the most general consistency equations which should be satisfied by the effective one-body eigenstates and show that they reduce to the spectral split consistency equations for the appropriate initial conditions.

Tabulated rates for astrophysical photodisintegration reactions make use of Boltzmann statistics for the photons involved as well as the interacting nuclei. Here we derive analytic corrections for the Planck-spectrum quantum statistics of the photon energy distribution. These corrections can be deduced directly from the detailed-balance condition without the assumption of equilibrium as long as the photons are represented by a Planck spectrum. Moreover we show that these corrections affect not only the photodisintegration rates but also modify the conditions of nuclear statistical equilibrium as represented in the Saha equation. We deduce new analytic corrections to the classical Maxwell-Boltzmann statistics which can easily be added to the reverse reaction rates of existing reaction network tabulations. We show that the effects of quantum statistics, though generally quite small, always tend to speed up photodisintegration rates and are largest for nuclei and environments for which Q/kT ~ 1. As an illustration, we examine possible effects of these corrections on the r-process, the rp-process, explosive silicon burning, the $\gamma$-process and big bang nucleosynthesis. We find that in most cases one is quite justified in neglecting these corrections. The correction is largest for reactions near the drip line for an r-process with very high neutron density, or an rp-process at high-temperature.

An algebraic approach to the neutrino propagation in dense media is presented. The Hamiltonian describing a gas of neutrinos interacting with each other and with background fermions is written in terms of the appropriate SU(N) operators, where N is the number of neutrino flavors. The evolution of the resulting many-body problem is formulated as a coherent-state path integral. Some commonly used approximations are shown to represent the saddle-point solution of the path integral for the full many-body system.