We compute the production rate of the energy density carried by gravitational waves emitted by a Standard Model plasma in thermal equilibrium, consistently to leading order in coupling constants for momenta $k\sim \pi T$. Summing up the contributions from the full history of the universe, the highest temperature of the radiation epoch can be constrained by the so-called $N_{\rm eff}$ parameter. The current theoretical uncertainty $\Delta N_{\rm eff} \le 10^{-3}$ corresponds to $T_{\rm max} \le 2\times 10^{17}$ GeV. In the course of the computation, we show how a subpart of the production rate can be determined with the help of standard packages, even if subsequently an IR subtraction and thermal resummation need to be implemented.

The putative recent indication of an unidentified 3.55 keV X-ray line in certain astrophysical sources is taken as a motivation for an improved theoretical computation of the cosmological abundance of 7.1 keV sterile neutrinos. If the line is interpreted as resulting from the decay of Warm Dark Matter, the mass and mixing angle of the sterile neutrino are known. Our computation then permits for a determination of the lepton asymmetry that is needed for producing the correct abundance via the Shi-Fuller mechanism, as well as for an estimate of the non-equilibrium spectrum of the sterile neutrinos. The latter plays a role in structure formation simulations. Results are presented for different flavour structures of the neutrino Yukawa couplings and for different types of pre-existing lepton asymmetries, accounting properly for the charge neutrality of the plasma and incorporating approximately hadronic contributions.

Because of physical processes ranging from microscopic particle collisions to macroscopic hydrodynamic fluctuations, any plasma in thermal equilibrium emits gravitational waves. For the largest wavelengths the emission rate is proportional to the shear viscosity of the plasma. In the Standard Model at T > 160 GeV, the shear viscosity is dominated by the most weakly interacting particles, right-handed leptons, and is relatively large. We estimate the order of magnitude of the corresponding spectrum of gravitational waves. Even though at small frequencies (corresponding to the sub-Hz range relevant for planned observatories such as eLISA) this background is tiny compared with that from non-equilibrium sources, the total energy carried by the high-frequency part of the spectrum is non-negligible if the production continues for a long time. We suggest that this may constrain (weakly) the highest temperature of the radiation epoch. Observing the high-frequency part directly sets a very ambitious goal for future generations of GHz-range detectors.