We study the dynamics of active-sterile neutrino oscillations in the early universe using full momentum-dependent quantum-kinetic equations. These equations are too complicated to allow for an analytical treatment, and numerical solution is greatly complicated due to very pronounced and narrow structures in the momentum variable introduced by resonances. Here we introduce a novel dynamical discretization of the momentum variable which overcomes this problem. As a result we can follow the evolution of neutrino ensemble accurately well into the stable growing phase. Our results confirm the existence of a "chaotic region" of mixing parameters, for which the final sign of the asymmetry, and hence the SBBN prediction of He(4)-abundance cannot be accurately determined.

We consider active-sterile neutrino oscillations in the early universe in an inhomogeneous isocurvature background. We show that very small initial baryonic seed-inhomogeneities can trigger a growth of very large amplitude spatial fluctuations in lepton asymmetry. Domains of varying asymmetry are observed to persist for a long time despite dissipative effects. The space dependent asymmetry profiles give rise to MSW-resonances within the domain boundaries, enhancing dramatically the equilibration of the sterile neutrino species. According to our one-dimensional toy-model, the effect is so strong that almost the entire parameter space where exponential growth of asymmetry occurs would be ruled out by nucleosynthesis.

We consider the indeterminacy in the sign of the neutrino asymmetry generated by active-sterile neutrino oscillations in the early universe. The dynamics of asymmetry growth is discussed in detail and the indeterminacy in the final sign of the asymmetry is shown to be a real physical phenomenon. Recently published contradicting results are carefully considered and the underlying assumptions leading to the disagreement are resolved.

We consider numerically the growth of neutrino asymmetry in active-sterile neutrino oscillations in the early universe. It is shown that the final sign of the asymmetry can be highly sensitive to small variations of the oscillation parameters. We find regions which are completely or partially chaotic, but also regions where the sign remains very robust. The consequences for atmospheric neutrino oscillations and primordial nucleosynthesis are then discussed. In the completely chaotic region the predicted He(4)-abundance has an inherent arbitrariness \Delta Y \simeq 0.01.