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.