We present a full result for the equation of state (EoS) in 2+1+1 (up/down, strange and charm quarks are present) flavour lattice QCD. We extend this analysis and give the equation of state in 2+1+1+1 flavour QCD. In order to describe the evolution of the universe from temperatures several hundreds of GeV to several tens of MeV we also include the known effects of the electroweak theory and give the effective degree of freedoms. As another application of lattice QCD we calculate the topological susceptibility (chi) up to the few GeV temperature region. These two results, EoS and chi, can be used to predict the dark matter axion's mass in the post-inflation scenario and/or give the relationship between the axion's mass and the universal axionic angle, which acts as a initial condition of our universe.

Nonperturbative approximation schemes based on two-particle irreducible (2PI) effective actions provide an important means for our current understanding of (non-)equilibrium quantum field theory. A remarkable property is their renormalizability, since these approximations involve selective summations to infinite perturbative orders. In this paper we show how to renormalize all $n$-point functions of the theory, which are given by derivatives of the 2PI-resummed effective action $\Gamma[\phi]$ for scalar fields $\phi$. This provides a complete description in terms of the generating functional for renormalized proper vertices, which extends previous prescriptions in the literature on the renormalization for 2PI effective actions. The importance of the 2PI-resummed generating functional for proper vertices stems from the fact that the latter respect all symmetry properties of the theory and, in particular, Goldstone's theorem in the phase with spontaneous symmetry breaking. This is important in view of the application of these techniques to gauge theories, where Ward identities play a crucial role.

Isotropization occurs on time scales much shorter than the thermal equilibration time. This is a crucial ingredient for the understanding of collision experiments of heavy nuclei or other nonequilibrium phenomena in complex many body systems. We discuss in detail the limitations of estimates based on standard ``linear'' or relaxation-time approximations, where isotropization and thermal equilibration rates agree. For a weak-coupling $\phi^4$-model the relaxation-time approximation underestimates the thermal equilibration time by orders of magnitude, in contrast to the isotropization time. The characteristic nonequilibrium isotropization rate can be enhanced as compared to the close-to-equilibrium value. Our results are obtained from the two-particle irreducible effective action, which includes off-shell and memory effects and does not involve a gradient expansion. This allows us to determine the range of validity of a description to lowest-order in gradients, which is typically employed in kinetic equations.

We solve the nonequilibrium dynamics of a 3+1 dimensional theory with Dirac fermions coupled to scalars via a chirally invariant Yukawa interaction. The results are obtained from a systematic coupling expansion of the 2PI effective action to lowest non-trivial order, which includes scattering as well as memory and off-shell effects. The dynamics is solved numerically without further approximation, for different far-from-equilibrium initial conditions. The late-time behavior is demonstrated to be insensitive to the details of the initial conditions and to be uniquely determined by the initial energy density. Moreover, we show that at late time the system is very well characterized by a thermal ensemble. In particular, we are able to observe the emergence of Fermi--Dirac and Bose--Einstein distributions from the nonequilibrium dynamics.