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 manybody 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 manybody 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.

We show that nuclear pairing Hamiltonian exhibits supersymmetry in the strong-coupling limit. The underlying supersymmetric quantum mechanical structure explains the degeneracies between the energies of the N and Nmax − N + 1 pair eigenstates. The supersymmetry transformations connecting these states are given.

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.

We apply the algebraic Bethe ansatz technique to the nuclear pairing problem with orbit dependent coupling constants and degenerate single particle energy levels. We find the exact energies and eigenstates. We show that for a given shell, there are degeneracies between the states corresponding to less and more than half full shell. We also provide a technique to solve the equations of Bethe ansatz.