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.

Studies of non-equilibrium dynamics of cosmological phase transitions may involve a scalar field interacting weakly with the energy-momentum tensor of a thermal plasma. At late times, when the scalar field is approaching equilibrium, it experiences both damping and thermal fluctuations. We show that thermal fluctuations induce a shear viscosity and a gravitational wave production rate, and propose that including this tunable contribution may help in calibrating the measurement of the gravitational wave production rate in hydrodynamic simulations. Furthermore it may enrich their physical scope, permitting in particular for a study of the instability of growing bubbles.

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

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.

We determine the abundance of the lightest (dark matter) sterile neutrinos created in the Early Universe due to active-sterile neutrino transitions from the thermal plasma. Our starting point is the field-theoretic formula for the sterile neutrino production rate, derived in our previous work [JHEP 06(2006)053], which allows to systematically incorporate all relevant effects, and also to analyse various hadronic uncertainties. Our numerical results differ moderately from previous computations in the literature, and lead to an absolute upper bound on the mixing angles of the dark matter sterile neutrino. Comparing this bound with existing astrophysical X-ray constraints, we find that the Dodelson-Widrow scenario, which proposes sterile neutrinos generated by active-sterile neutrino transitions to be the sole source of dark matter, is only possible for sterile neutrino masses lighter than 3.5 keV (6 keV if all hadronic uncertainties are pushed in one direction and the most stringent X-ray bounds are relaxed by a factor of two). This upper bound may conflict with a lower bound from structure formation, but a definitive conclusion necessitates numerical simulations with the non-equilibrium momentum distribution function that we derive. If other production mechanisms are also operative, no upper bound on the sterile neutrino mass can be established.

It has been pointed out in ref.[1] that in the nuMSM (Standard Model extended by three right-handed neutrinos with masses smaller than the electroweak scale), there is a corner in the parameter space where CP-violating resonant oscillations among the two heaviest right-handed neutrinos continue to operate below the freeze-out temperature of sphaleron transitions, leading to a lepton asymmetry which is considerably larger than the baryon asymmetry. Consequently, the lightest right-handed (``sterile'') neutrinos, which may serve as dark matter, are generated through an efficient resonant mechanism proposed by Shi and Fuller [2]. We re-compute the dark matter relic density and non-equilibrium momentum distribution function in this situation with quantum field theoretic methods and, confronting the results with existing astrophysical data, derive bounds on the properties of the lightest right-handed neutrinos. Our spectra can be used as an input for structure formation simulations in warm dark matter cosmologies, for a Lyman-alpha analysis of the dark matter distribution on small scales, and for studying the properties of haloes of dwarf spheroidal galaxies.

Sterile neutrinos with masses in the keV range are considered to be a viable candidate for warm dark matter. The rate of their production through active-sterile neutrino transitions peaks, however, at temperatures of the order of the QCD scale, which makes it difficult to estimate their relic abundance quantitatively, even if the mass of the sterile neutrino and its mixing angle were known. We derive here a relation, valid to all orders in the strong coupling constant, which expresses the production rate in terms of the spectral function associated with active neutrinos. The latter can in turn be expressed as a certain convolution of the spectral functions related to various mesonic current-current correlation functions, which are being actively studied in other physics contexts. In the naive weak coupling limit, the appropriate Boltzmann equations can be derived from our general formulae.

Relativistic scalar field theories with a conserved global charge Q possess often (meta)stable spherically symmetric soliton solutions, called Q-balls. We elaborate on the perfect formal analogy which exists between Q-balls, and spherically symmetric solitons in certain non-relativistic atomic Bose-Einstein condensates, for which the dominant interatomic interaction can be tuned attractive. In a harmonic trap, present in existing experiments, the Q-ball solution is modified in an essential way. If the trap is significantly prolongated in one direction, however, then genuine solitons do appear, and actual experimental data can be obtained for some of the Q-ball properties studied numerically in the relativistic cosmological context, such as their formation and collisions. We also suggest conditions under which the same cosmologically relevant analogies could be extended to the fully three-dimensional case.

Scalar electrodynamics can be used to investigate the formation of cosmic strings in the early universe. We present the results of lattice Monte Carlo simulations of an effective three-dimensional U(1)+Higgs theory that describes the equilibrium properties of finite-temperature scalar electrodynamics near the transition. A gauge-invariant criterion for the existence of a vortex is used in measuring the properties of the vortex network in the equilibrium state both in the Coulomb and in the Higgs phase of the system. The naive definition of the vortex density becomes meaningless in the continuum limit and special care is needed in extracting physical quantities. Numerical evidence for a physical discontinuity in the vortex density is given.

We study numerically the real-time history of the cosmological electroweak phase transition, as it may take place in the Standard Model or in MSSM for m_H < m_W according to recent lattice results. We follow the nucleated bubbles from the initial stages of acceleration and rapid growth, through collisions with compression waves resulting in slowing down and reheating to T_c, until the final stages of slow growth and evaporation. We find that collisions with compression waves may make the bubble walls oscillate in the radial direction, and that reheating to T_c takes generically place.

We study spherically symmetric bubble growth and droplet decay in first order cosmological phase transitions, using a numerical code including both the complete hydrodynamics of the problem and a phenomenological model for the microscopic entropy producing mechanism at the phase transition surface. The small-scale effects of finite wall width and surface tension are thus consistently incorporated. We verify the existence of the different hydrodynamical growth modes proposed recently and investigate the problem of a decaying quark droplet in the QCD phase transition. We find that the decaying droplet leaves behind no rarefaction wave, so that any baryon number inhomogeneity generated previously should survive the decay.

The classification of the hydrodynamical growth mechanisms for the spherical bubbles of the low-temperature phase in cosmological phase transitions is completed by showing that the bubbles can grow as supersonic deflagrations. Such deflagrations consist of a Jouguet deflagration, followed by a rarefaction wave. Depending on the amount of supercooling, the maximal velocity of supersonic deflagrations varies between the sound and the light velocities. The solutions faster than supersonic deflagrations are weak detonations.

We examine the first-order cosmological QCD phase transition for a large class of parameter values, previously considered unlikely. We find that the hadron bubbles can nucleate at very large distance scales, they can grow as detonations as well as deflagrations, and that the phase transition may be completed without reheating to the critical temperature. For a subset of the parameter values studied, the inhomogeneities generated at the QCD phase transition might have a noticeable effect on nucleosynthesis.

We study how bubbles grow after the initial nucleation event in generic first-order cosmological phase transitions characterised by the values of latent heat, interface tension and correlation length, and driven by a scalar order parameter $\phi$. Equations coupling $\phi$ and the fluid variables $v$ and $T$ and depending on a dissipative constant $\Gamma$ are derived and solved numerically in the 1+1 dimensional case starting from a slightly deformed critical bubble configuration. Parameters corresponding to QCD and electroweak phase transitions are chosen and the whole history of the bubble with formation of combustion and shock fronts is computed as a function of $\Gamma$. Both deflagrations and detonations can appear depending on the values of the parameters. Reheating due to collisions of bubbles is also computed.