Weakly Interactive Massive Particles (WIMPs) are the most widely studied candidate particles forming the cold dark matter (CDM) whose existence can be inferred from a wealth of astrophysical and cosmological observations. In the framework of the minimal cosmological model detailed measurements on the cosmic microwave background by the PLANCK collaboration fix the scaled CDM relic density to Ωch 2 = 0.1193 ± 0.0014, with an error of less than 1.5%. In order to fully exploit this observational precision, theoretical calculations should have a comparable or smaller error. In this paper we use recent lattice QCD calculations to improve the description of the thermal plasma. This affects the predicted relic density of "thermal WIMPs", which once were in chemical equilibrium with Standard Model particles. For WIMP masses between 3 and 15 GeV, where QCD effects are most important, our predictions differ from earlier results by up to 9% (12%) for pure S−wave (P −wave) annihilation. We use these results to compute the thermally averaged WIMP annihilation cross section that reproduces the correct CDM relic density, for WIMP masses between 0.1 GeV and 10 TeV.

We investigate dark matter (DM) production in an early matter dominated era where a heavy long-lived particle decays to radiation and DM. In addition to DM annihilation into and thermal DM production from radiation, we include direct DM production from the decay of the long-lived particle. In contrast to earlier treatments the temperature dependence of the number of degrees of freedom g * in the Standard Model (SM) plasma is treated carefully. Besides the well-known cases of thermal hot and cold DM, additional regions of parameter space with the approximately correct DM relic density appear. In some of these regions the temperature dependence of g * can change the final DM density by several hundred percent. Furthermore, we analyze the effect of allowing nonvanishing initial abundances for radiation and DM. We find an upper bound on the mass of the long-lived particle if the DM annihilation cross section is below that corresponding to thermal WIMP (Weakly Interactive Massive Particle) DM in standard cosmology.

Assuming that inflation is followed by a phase where the energy density of the Universe is dominated by a component with a general equation of state, we evaluate the spectrum of primordial gravitational waves induced in the post-inflationary Universe. We show that if the energy density of the Universe is dominated by a component φ before Big Bang nucleosynthesis, its equation of state could be constrained by gravitational wave experiments depending on the ratio of energy densities of φ and radiation, and also the temperature at the end of the φ dominated era. Also, we discuss the impact of scale dependence of tensor modes on the primordial gravitational wave spectrum during the φ-domination. These models are motivated by beyond Standard Model physics and scenarios for non-thermal production of dark matter in the early Universe. We also constrain the parameter space of the tensor spectral index and the tensor-to-scalar ratio, using the experimental limits from gravitational wave experiments.

We study the induced primordial gravitational waves (GW) coming from the effect of scalar perturbation on the tensor perturbation at the second order of cosmological perturbation theory. We use the evolution of the standard model degrees of freedom with respect to temperature in the early Universe to compute the induced gravitational waves bakcground. Our result shows that the spectrum of the induced GW is affected differently by the standard model degrees of freedom than the GW coming from first order tensor perturbation. This phenomenon is due to the presence of scalar perturbations as a source for tensor perturbations and it is effective around the quark gluon deconfinement and electroweak transition. In case of considering a scalar spectral index larger than one at small scales or a non-Gaussian curvature power spectrum this effect can be observed by gravitational wave observatories.

Thermal axion production in the early universe goes through several mass thresholds, and the resulting rate may change dramatically across them. Focusing on the KSVZ and DFSZ frameworks for the invisible QCD axion, we perform a systematic analysis of thermal production across thresholds and provide smooth results for the rate. The QCD phase transition is an obstacle for both classes of models. For the hadronic KSVZ axion, we also deal with production at temperatures around the mass of the heavy-colored fermion charged under the Peccei-Quinn symmetry. Within the DFSZ framework, standard model fermions are charged under this symmetry, and additional thresholds are the heavy Higgs bosons masses and the electroweak phase transition. We investigate the cosmological implications with a specific focus on axion dark radiation quantified by an effective number of neutrino species and explore the discovery reach of future CMB-S4 surveys.