We study the impact of non-standard momentum distributions of cosmic neutrinos on the anisotropy spectrum of the cosmic microwave background and the matter power spectrum of the large scale structure. We show that the neutrino distribution has almost no unique observable imprint, as it is almost entirely degenerate with the effective number of neutrino flavours, $N_{\mathrm{eff}}$, and the neutrino mass, $m_{\nu}$. Performing a Markov chain Monte Carlo analysis with current cosmological data, we demonstrate that the neutrino mass bound heavily depends on the assumed momentum distribution of relic neutrinos. The message of this work is simple and has to our knowledge not been pointed out clearly before: Cosmology allows that neutrinos have larger masses if their average momentum is larger than that of a perfectly thermal distribution. Here we provide an example in which the mass limits are relaxed by a factor of two.

Constraints on inflationary models typically assume only the standard models of cosmology and particle physics. By extending the neutrino sector to include a new interaction with a light scalar mediator ($m_{\phi}\sim$MeV), it is possible to relax these constraints, in particular via opening up regions of the parameter space of the spectral index $n_s$. These new interactions can be probed at IceCube via interactions of astrophysical neutrinos with the Cosmic Neutrino Background for nearly all of the relevant parameter space.

We investigate how a lepton asymmetry impacts the cosmic trajectory in the QCD phase diagram. We study the evolution of chemical potentials during the QCD epoch of the early Universe using susceptibilities from lattice QCD to interpolate between an ideal quark gas and an ideal hadron resonance gas. The lepton asymmetry affects the evolution of all chemical potentials. The standard cosmic trajectory is obtained assuming tiny lepton and baryon asymmetries. For larger lepton asymmetry, the charge chemical potential exceeds the baryon chemical potential before pion annihilation.

We consider the impact of neutrino self-interactions described by an effective four-fermion coupling on cosmological observations. Implementing the exact Boltzmann hierarchy for interacting neutrinos first derived in [arxiv:1409.1577] into the Boltzmann solver CLASS, we perform a detailed numerical analysis of the effects of the interaction on the cosmic microwave background (CMB) anisotropies, and compare our results with known approximations in the literature. While we find good agreement between our exact approach and the relaxation time approximation used in some recent studies, the popular $\left( c_{\text{eff}}^2,c_{\text{vis}}^2 \right)$-parameterisation fails to reproduce the correct scale dependence of the CMB temperature power spectrum. We then proceed to derive constraints on the effective coupling constant $G_{\text{eff}}$ using currently available cosmological data via an MCMC analysis. Interestingly, our results reveal a bimodal posterior distribution, where one mode represents the standard $\Lambda$CDM limit with $G_{\rm eff} \lesssim 10^8 \, G_{\rm F}$, and the other a scenario in which neutrinos self-interact with an effective coupling constant $G_{\rm eff} \simeq 3 \times 10^9 \, G_{\rm F}$.

A cosmic lepton asymmetry $\eta_{\text{l}}=(n_{\text{l}}-n_{\bar{\text{l}}})/n_{\gamma}$ affects the primordial helium abundance and the expansion rate of the early Universe. Both of these effects have an impact on the anisotropies of the cosmic microwave background (CMB). We derive constraints on the neutrino chemical potentials from the Planck 2015 data, assuming equal lepton flavour asymmetries and negligible neutrino masses. We find $\xi=-0.002 ^{+0.114}_{-0.111}$ (95\% CL) for the chemical potentials, which corresponds to $ -0.085 \leq \eta_{\text{l}} \leq 0.084$. Our constraints on the lepton asymmetry are significantly stronger than previous constraints from CMB data analysis and we argue that they are more robust than those from primordial light element abundances. The resulting constraints on the primordial helium and deuterium abundances are consistent with those from direct measurements.

Cosmic reionization and dark matter decay can impact observations of the cosmic microwave sky in a similar way. A simultaneous study of both effects is required to constrain unstable dark matter from cosmic microwave background observations. We compare two reionization models with and without dark matter decay. We find that a reionization model that fits also data from quasars and star forming galaxies results in tighter constraints on the reionization optical depth $\tau_{\text{reio}}$, but weaker constraints on the spectral index $n_{\text{s}}$ than the conventional parametrization. We use the Planck 2015 data to constrain the effective decay rate of dark matter to $\Gamma_{\rm eff} < 2.9 \times 10^{-25}/$s at $95$\% C.L. This limit is robust and model independent. It holds for any type of decaying dark matter and it depends only weakly on the chosen parametrization of astrophysical reionization. For light dark matter particles that decay exclusively into electromagnetic components this implies a limit of $\Gamma < 5.3 \times 10^{-26}/$s at $95$\% C.L. Specifying the decay channels, we apply our result to the case of keV-mass sterile neutrinos as dark matter candidates and obtain constraints on their mixing angle and mass, which are comparable to the ones from the diffuse X-ray background.

Starting from the collisional Boltzmann equation, we derive for the first time and from first principles the Boltzmann hierarchy for neutrinos including interactions with a scalar particle. Such interactions appear, for example, in majoron-like models of neutrino mass generation. We study two limits of the scalar mass: (i) An extremely massive scalar whose only role is to mediate an effective 4-fermion neutrino-neutrino interaction, and (ii) a massless scalar that can be produced in abundance and thus demands its own Boltzmann hierarchy. In contrast to, e.g., the first-order Boltzmann hierarchy for Thomson-scattering photons, our interacting neutrino/scalar Boltzmann hierarchies contain additional momentum-dependent collision terms arising from a non-negligible energy transfer in the neutrino-neutrino and neutrino-scalar interactions. This necessitates that we track each momentum mode of the phase space distributions individually, even if the particles were massless. Comparing our hierarchy with the commonly used $(c_{\rm eff}^2,c_{\rm vis}^2)$-parameterisation, we find no formal correspondence between the two approaches, which raises the question of whether the latter parameterisation even has an interpretation in terms of particle scattering. Lastly, although we have invoked majoron-like models as a motivation for our study, our treatment is in fact generally applicable to all scenarios in which the neutrino and/or other ultrarelativistic fermions interact with scalar particles.