In the standard model extended with a seesaw mass matrix, we study the production of sterile neutrinos from the decay of vector bosons at $T\simeq M_{W,Z}$. We derive a general quantum kinetic equation for the production of sterile neutrinos and their effective mixing angles valid in a wide range of temperature, to all orders in interactions of the standard model, and to leading order mixing angle $\theta \ll 1$. Production rates and effective mixing angles depend sensitively on helicity. Positive helicity states interact more weakly with the medium and their effective mixing angle is not modified significantly whereas the mixing angle for negative helicity is strongly suppressed by the medium. If $M_s \lesssim 8.35\,\mathrm{MeV}$, there are fewer states with negative helicity produced than those with positive helicity. There is an MSW resonance in the absence of lepton asymmetry, but is screened by the damping rate, without production enhancement. Negative helicity states freeze-out at $T^-_f\simeq 5\,\mathrm{GeV}$ and positive helicity states freeze-out at $T^+_f \simeq 8\,\mathrm{GeV}$, both distributions are far from thermal. Negative helicity states feature a broader momentum distribution than that for those with positive helicity. Sterile neutrinos produced via vector boson decay do not satisfy the abundance, lifetime and cosmological constraints to be the sole dark matter component in the universe but might solve the $^{7}Li$ problem, albeit at the very edge of the possible parameter space. A heavy sterile neutrino with a mass of a few MeV could decay into light sterile neutrinos, of a few keV in mass, that contribute to warm dark matter. We argue that heavy sterile neutrinos with lifetime $\leq 1/H_0$ reach local thermodynamic equilibrium.

Sterile neutrinos are SU(2) singlets that mix with active neutrinos via a mass matrix, its diagonalization leads to mass eigenstates that couple via standard model vertices. We study the cosmological production of heavy neutrinos via standard model charged and neutral current vertices under a minimal set of assumptions: i) the mass basis contains a hierarchy of heavy neutrinos, ii) these have very small mixing angles with the active (flavor) neutrinos, iii) standard model particles, including light (active-like) neutrinos are in thermal equilibrium. The same weak interaction processes that produce active-like neutrinos also produce the heavier species. We introduce the kinetic equations that describe their production, freeze out and decay and discuss the various processes that lead to their production in a wide range of temperatures assessing their feasibility as dark matter candidates. We identify processes in which finite temperature collective excitations lead to the production of the heavy species. As a specific example, we consider the production of heavy neutrinos from pion decay shortly after the QCD crossover including finite temperature corrections to the pion form factors and mass. We consider the different decay channels that allow for the production of heavy neutrinos showing that their frozen distribution functions exhibit effects from "kinematic entanglement" and argue for their viability as mixed dark matter candidates. We discuss abundance, phase space density and stability constraints and argue that heavy neutrinos with lifetime $\tau> 1/H_0$ freeze out of local thermal equilibrium, and \emph{conjecture} that those with lifetimes $\tau \ll 1/H_0$ may undergo cascade decay into lighter DM candidates and/or inject non-LTE neutrinos into the cosmic neutrino background. We provide a comparison with non-resonant production via active-sterile mixing.

We study the production of sterile neutrinos in the early universe from $\pi \rightarrow l \nu_s$ shortly after the QCD phase transition in the absence of a lepton asymmetry while including finite temperature corrections to the $\pi$ mass and decay constant $f_{\pi}$. Sterile neutrinos with masses $\lesssim 1 MeV$ produced via this mechanism freeze-out at $T_f \simeq 10 MeV$ with a distribution function that is highly non-thermal and features a sharp enhancement at low momentum thereby making this species \emph{cold} even for very light masses. Dark matter abundance constraints from the CMB and phase space density constraints from the most dark matter dominated dwarf spheroidal galaxies provide upper and lower bounds respectively on combinations of mass and mixing angles. For $\pi \rightarrow \mu \nu_s$, the bounds lead to a narrow region of compatibility with the latest results from the $3.55 \mathrm{KeV}$ line. The non-thermal distribution function leads to free-streaming lengths (today) in the range of $\sim \mbox{few kpc}$ consistent with the observation of cores in dwarf galaxies. For sterile neutrinos with mass $\lesssim 1 eV$ that are produced by this reaction, the most recent accelerator and astrophysical bounds on $U_{ls}$ combined with the non-thermal distribution function suggests a substantial contribution from these sterile neutrinos to $N_{eff}$.

We study corrections to the power spectra of curvature and tensor perturbations and the tensor-to-scalar ratio in single field slow roll inflation due to initial conditions imprinted by a fast-roll stage prior to slow roll. For a wide range of initial inflaton kinetic energy, this stage lasts only a few e-folds and merges smoothly with slow-roll leading to non-Bunch-Davies initial conditions for modes that exit the Hubble radius during slow roll. We describe a program that yields the dynamics in the fast-roll stage while matching to the slow roll stage independent of the inflationary potentials. Corrections to the power spectra are encoded in a transfer function $\mathcal{T}_\alpha(k)$. We obtain $\mathcal{T}_\alpha(k)$ to leading order in a Born approximation valid for modes of observational relevance today. A fit yields $\mathcal{T}_\alpha(k) =1+ A_{\alpha}k^{-p}\cos[2\pi \omega k/H_{sr}+\varphi_\alpha]$, with $1.5 \lesssim p \lesssim 2$, $\omega \simeq 1$ and $H_{sr}$ the Hubble scale during slow roll inflation, where curvature and tensor perturbations feature the same $p,\omega$ for a wide range of initial conditions. These corrections lead to both a suppression of the quadrupole and oscillatory features in both $P_R(k)$ and $r(k_0)$ with a period of the order of the Hubble scale during slow roll inflation. The results are independent of the specific inflationary potentials, depending solely on the ratio of kinetic to potential energy $\kappa$ and the slow roll parameters to leading order. For a wide range of $\kappa$ and the values of $\epsilon_V;\eta_V$ corresponding to the upper bounds from Planck, we find that the low quadrupole is consistent with the results from Planck, and the oscillations in $r(k_0)$ could be observable if the modes corresponding to the quadrupole and the pivot scale crossed the Hubble radius a few e-folds after the onset of slow roll.

If the large scale anomalies in the temperature power spectrum of the cosmic microwave background are of primordial origin, they may herald modifications to the slow roll inflationary paradigm on the largest scales. We study the possibility that the origin of the large scale power suppression is a modification of initial conditions during slow roll as a result of a pre-slow roll phase during which the inflaton evolves rapidly. This stage is manifest in a potential in the equations for the Gaussian fluctuations during slow roll and modify the power spectra of scalar perturbations via an initial condition transfer function $\mathcal{T}(k)$. We provide a general analytical study of its large and small scale properties and analyze the impact of these initial conditions on the infrared aspects of typical test scalar fields. The infrared behavior of massless minimally coupled test scalar field theories leads to the dynamical generation of mass and anomalous dimensions, both depend non-analytically on $\mathcal{T}(0)$. During inflation all quanta decay into many quanta even of the same field because of the lack of kinematic thresholds. The decay leads to a quantum entangled state of sub and superhorizon quanta with correlations across the horizon. We find the modifications of the decay width and the entanglement entropy from the initial conditions. In all cases, initial conditions from a ``fast-roll'' stage that lead to a suppression in the scalar power spectrum at large scales also result in a suppression of the dynamically generated masses, anomalous dimensions and decay widths.

In inflationary cosmology all particle states decay as a consequence of the lack of kinematic thresholds. The decay of an initial single particle state yields an \emph{entangled quantum state of the product particles}. We generalize and extend a manifestly unitary field theoretical method to obtain the time evolution of the quantum state. We consider the decay of a light scalar field with mass $M\ll H$ with a cubic coupling in de Sitter space-time. Radiative corrections feature an infrared enhancement manifest as poles in $\Delta=M^2/3H^2$ and we obtain the quantum state in an expansion in $\Delta$. To leading order in $\Delta$ the pure state density matrix describing the decay of a particle with sub-horizon wavevector is dominated by the emission of superhorizon quanta, describing \emph{entanglement between superhorizon and subhorizon fluctuations and correlations across the horizon}. Tracing over the superhorizon degrees of freedom yields a mixed state density matrix from which we obtain the entanglement entropy. Asymptotically this entropy grows with the \emph{physical} volume as a consequence of more modes of the decay products crossing the Hubble radius. A generalization to localized wave packets is provided. The cascade decay of single particle states into many particle states is discussed. We conjecture on \emph{possible} impact of these results on non-gaussianity and on the ``low multipole anomalies'' of the CMB.