We consider the cosmological production of fermionic dark matter during inflation and a post-inflationary radiation dominated era. This fermion \emph{only interacts gravitationally}, has a mass $m$ much smaller than the Hubble scale during inflation (but is otherwise arbitrary) and is in its Bunch-Davies vacuum state during inflation. We focus on superhorizon modes at the end of inflation, and assume instantaneous reheating. We obtain the full energy momentum tensor discussing its renormalization, and show that the contribution from particle production is of the kinetic-fluid form near matter-radiation equality. We find \emph{exactly} the distribution function of produced particles $|B(k)|^2=\frac{1}{2}\Big[1-(1-e^{-\frac{k^2}{2mT_H}})^{1/2}\Big]$ which exhibits an " emergent temperature" $T_H=H_0\sqrt{\Omega_R}\simeq 10^{-36}(\mathrm{eV})$. The abundance of the produced particles $\Omega_{pp}$ is very similar to that of a non-relativistic degree of freedom thermalized at temperature $T_H$, $\Omega_{pp} \propto m \, (m\,T_H)^{3/2}\simeq \big(m/10^8\mathrm{GeV})\big)^{5/2}$ and "cold" equation of state $w(a) \simeq (T_H/m a^2)$, both dominated by superhorizon modes at the end of inflation. We discuss subtle aspects of isocurvature perturbations.

We study the non-adiabatic cosmological production of ultra light dark matter (ULDM) under a minimal set of assumptions: a free ultra light real scalar as a spectator field in its Bunch-Davies vacuum state during inflation and instantaneous reheating into a radiation dominated era. For (ULDM) fields minimally coupled to gravity, non-adiabatic particle production yields a \emph{distribution function} peaked at \emph{low} comoving momentum $\mathcal{N}_k \propto 1/k^3$. The infrared behavior is a remnant of the infrared enhancement of light minimally coupled fields during inflation. We obtain the full energy momentum tensor, show explicity its equivalence with the fluid-kinetic one in the adiabatic regime, and extract the abundance, equation of state and free streaming length (cutoff in the matter power spectrum). Taking the upper bound on the scale of inflation from Planck, the (UDLM) saturates the dark matter abundance for $m \simeq 1.5\,\times 10^{-5}\mathrm{eV}$ with an equation of state parameter $w \simeq 10^{-14}$ and a free streaming length $\lambda_{fs} \simeq 70\,\mathrm{pc}$. Thus this cosmologically produced (ULDM) yields a \emph{cold} dark matter particle. We argue that the abundance from non-adiabatic production yields a \emph{lower bound} on generic (ULDM) and axion-like particles that must be included in any assessment of (ULDM) as a dark matter candidate.

We study the decay of a Higgs-like scalar Yukawa coupled to massless fermions in post-inflationary cosmology, combining a non-perturbative method with an adiabatic expansion. The renormalized survival probability $\mathcal{P}_\Phi(t)$ of a (quasi) particle ``born'' at time $t_b$ and decaying at rest in the comoving frame, $\mathcal{P}_\Phi(t) = \Big[\frac{t}{t_b}\Big]^{-\frac{Y^2}{8\pi^2}}~ e^{ \frac{Y^2}{4\pi^2}\,\big(t/t_b\big)^{1/4} } \,e^{-\Gamma_0\,(t-t_b)}~ \mathcal{P}_\Phi(t_b) $, with $\Gamma_0$ the decay rate at rest in Minkowski space-time. For an ultrarelativistic particle we find $\mathcal{P}_\Phi(t) = e^{-\frac{2}{3}\Gamma_0\,t_{nr}\,(t/t_{nr})^{3/2}}~ \mathcal{P}_\Phi(t_b)$ before it becomes non-relativistic at a time $t_{nr}$ as a consequence of the cosmological redshift. For $t\gg t_{nr}$ we find $\mathcal{P}_\Phi(t) = \Big[\frac{t}{t_{nr}}\Big]^{-\frac{Y^2}{8\pi^2}}~ e^{ \frac{Y^2}{4\pi^2}\,\big(t/t_{nr}\big)^{1/4} }~\Big[\frac{t}{t_{nr}}\Big]^{\Gamma_0 t_{nr}/2} \,e^{-\Gamma_0\,(t-t_{nr})}~ \mathcal{P}_\Phi(t_{nr})$. The extra power is a consequence of the memory on the past history of the decay process. We compare these results to an S-matrix inspired phenomenological Minkowski-like decay law modified by an instantaneous Lorentz factor to account for cosmological redshift. Such phenomenological description \emph{under estimates the lifetime of the particle}. For very long lived, very weakly coupled particles, we obtain an \emph{upper bound} for the survival probability as a function of redshift $z$ valid throughout the expansion history $\mathcal{P}_\Phi(z) \gtrsim e^{-\frac{\Gamma_0}{H_0}\,\Upsilon(z,z_b)}\,\mathcal{P}_\Phi(z_b)$, where $\Upsilon(z,z_b)$ only depends on cosmological parameters and $t_{nr}$.

We study scalar particle decay during the radiation and matter dominated epochs of a standard cosmological model. An adiabatic approximation is introduced that is valid for degrees of freedom with typical wavelengths much smaller than the particle horizon ($\propto$~Hubble radius) at a given time. We implement a non-perturbative method that includes the cosmological expansion and obtain a cosmological Fermi's Golden Rule that enables one to compute the decay law of a parent particle of mass $m_1$, along with the build up of the population of daughter particles of mass $m_2$. The survival probability of the decaying particle is $P(t)=e^{-\widetilde{\Gamma}_k(t)\,t}$ with $\widetilde{\Gamma}_k(t)$ being an \emph{effective momentum and time dependent decay rate}. It features a transition time scale $t_{nr}$ between the relativistic and non-relativistic regimes and for $k \neq 0$ is always smaller than the analogous rate in Minkowski spacetime, as a consequence of (local) time dilation and the cosmological redshift. For $t \ll t_{nr}$ the decay law is a "stretched exponential" $P(t) = e^{-(t/t^*)^{3/2}}$, whereas for the non-relativistic stage with $t \gg t_{nr}$, we find $P(t) = e^{-\Gamma_0 t}\,(t/t_{nr})^{\Gamma_0\,t_{nr}/2}$. The Hubble time scale $\propto 1/H(t)$ introduces an energy uncertainty $\Delta E \sim H(t)$ which relaxes the constraints of kinematic thresholds. This opens new decay channels into heavier particles for $2\pi E_k(t) H(t) \gg 4m^2_2-m^2_1$, with $E_k(t)$ the (local) comoving energy of the decaying particle. As the expansion proceeds this channel closes and the usual two particle thresholds restrict the decay kinematics.