
We compute the slope and curvature, at vanishing fourmomentum transfer
squared, of the leading order hadron vacuum polarization function, using
lattice QCD. Calculations are performed with 2+1+1 flavors of staggered
fermions directly at the physical values of the quark masses and in volumes of
linear extent larger than 6fm. The continuum limit is carried out using six
different lattice spacings. All connected and disconnected contributions are
calculated, up to and including those of the charm.

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 postinflation 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.

We study the finite temperature transition in QCD with two flavors of
dynamical fermions at a pseudoscalar pion mass of about 350 MeV. We use
lattices with temporal extent of $N_t$=8, 10 and 12. For the first time in the
literature a continuum limit is carried out for several observables with
dynamical overlap fermions. These findings are compared with results obtained
within the staggered fermion formalism at the same pion masses and extrapolated
to the continuum limit. The presented results correspond to fixed topology and
its effect is studied in the staggered case. Nice agreement is found between
the overlap and staggered results.

The existence and stability of atoms rely on the fact that neutrons are more
massive than protons. The measured mass difference is only 0.14\% of the
average of the two masses. A slightly smaller or larger value would have led to
a dramatically different universe. Here, we show that this difference results
from the competition between electromagnetic and mass isospin breaking effects.
We performed lattice quantumchromodynamics and quantumelectrodynamics
computations with four nondegenerate Wilson fermion flavors and computed the
neutronproton masssplitting with an accuracy of $300$ kiloelectron volts,
which is greater than $0$ by $5$ standard deviations. We also determine the
splittings in the $\Sigma$, $\Xi$, $D$ and $\Xi_{cc}$ isospin multiplets,
exceeding in some cases the precision of experimental measurements.

While electromagnetic and updown quark mass difference effects on octet
baryon masses are very small, they have important consequences. The stability
of the hydrogen atom against beta decay is a prominent example. Here we include
these effects by adding them to valence quarks in a lattice QCD calculation
based on $N_f=2+1$ simulations with 5 lattice spacings down to 0.054 fm,
lattice sizes up to 6 fm and average updown quark masses all the way down to
their physical value. This allows us to gain control over all systematic
errors, except for the one associated with neglecting electromagnetism in the
sea. We compute the octet baryon isomultiplet mass splittings, as well as the
individual contributions from electromagnetism and the updown quark mass
difference. Our results for the total splittings are in good agreement with
experiment.

We determine the equation of state of QCD for nonzero chemical potentials via
a Taylor expansion of the pressure. The results are obtained for N_f=2+1
flavors of quarks with physical masses, on various lattice spacings. We present
results for the pressure, interaction measure, energy density, entropy density,
and the speed of sound for small chemical potentials. At low temperatures we
compare our results with the Hadron Resonance Gas model. We also express our
observables along trajectories of constant entropy over particle number. A
simple parameterization is given (the Matlab/Octave script parameterization.m,
submitted to the arXiv along with the paper), which can be used to reconstruct
the observables as functions of T and mu, or as functions of T and S/N.

We present the equation of state (pressure, trace anomaly, energy density and
entropy density) of the SU(3) gauge theory from lattice field theory in an
unprecedented precision and temperature range. We control both finite size and
cutoff effects. The studied temperature window (0.7...1000 T_c) stretches from
the glueball dominated system into the perturbative regime, which allows us to
discuss the range of validity of these approaches. We also determine the
preferred renormalization scale of the Hard Thermal Loop scheme and we fit the
unknown g^6 order perturbative coefficient at extreme high temperatures T>100
T_c. We furthermore quantify the nonperturbative contribution to the trace
anomaly using a simple functional form. Our high precision data allows one to
have a complete theoretical description of the equation of state from T=0 all
the way to the phase transition, through the transition region into the
perturbative regime up to the StefanBoltzmann limit. We will discuss this
description, too.

We present the equation of state (pressure, trace anomaly, energy density and
entropy density) of the SU(3) gauge theory from lattice field theory in an
unprecedented precision and temperature range. We control both finite size and
cutoff effects. The studied temperature window ($0.7... 1000 T_c$) stretches
from the glueball dominated system into the perturbative regime, which allows
us to discuss the range of validity of these approaches. From the critical
couplings on fine lattices we get $T_c/\Lambdamsbar=1.26(7)$ and use this ratio
to express the perturbative free energy in $T_c$ units. We also determine the
preferred renormalization scale of the Hard Thermal Loop scheme and we fit the
unknown $g^6$ order perturbative coefficient at extreme high temperatures
$T>100T_c$. We furthermore quantify the nonperturbative contribution to the
trace anomaly using two simple functional forms.

We extend our previous study [Phys. Lett. B643 (2006) 46] of the crossover
temperatures (T_c) of QCD. We improve our zero temperature analysis by using
physical quark masses and finer lattices. In addition to the kaon decay
constant used for scale setting we determine four quantities (masses of the
\Omega baryon, K^*(892) and \phi(1020) mesons and the pion decay constant)
which are found to agree with experiment. This implies that independently of
which of these quantities is used to set the overall scale the same results
are obtained within a few percent. At finite temperature we use finer lattices
down to a <= 0.1 fm (N_t=12 and N_t=16 at one point). Our new results confirm
completely our previous findings. We compare the results with those of the
'hotQCD' collaboration.

We discuss the possible extension of the bosonic classical field theory
simulations to include fermions. This problem has been addressed in terms of
the inhomogeneous mean field approximation by Aarts and Smit. By performing a
stochastic integration of an equivalent set of equations we can extend the
original 1+1 dimensional calculations so that they become feasible in higher
dimensions. We test the scheme in 2 + 1 dimensions and discuss some classical
applications with fermions for the first time, such as the decay of oscillons.

We review recent developments and open questions for the description of
nonequilibrium quantum fields, continuing hepph/0302210 and hepph/0410330.

We investigate lattice simulations of scalar and nonabelian gauge fields in
Minkowski spacetime. For SU(2) gaugetheory expectation values of link
variables in 3+1 dimensions are constructed by a stochastic process in an
additional (5th) ``Langevintime''. A sufficiently small Langevin step size and
the use of a tilted realtime contour leads to converging results in general.
All fixed point solutions are shown to fulfil the infinite hierarchy of
DysonSchwinger identities, however, they are not unique without further
constraints. For the nonabelian gauge theory the thermal equilibrium fixed
point is only approached at intermediate Langevintimes. It becomes more stable
if the complex time path is deformed towards Euclidean spacetime. We analyze
this behavior further using the realtime evolution of a quantum anharmonic
oscillator, which is alternatively solved by diagonalizing its Hamiltonian.
Without further optimization stochastic quantization can give accurate
descriptions if the realtime extend of the lattice is small on the scale of
the inverse temperature.

Calculations of nonequilibrium processes become increasingly feasable in
quantum field theory from first principles. There has been important progress
in our analytical understanding based on 2PI generating functionals. In
addition, for the first time direct lattice simulations based on stochastic
quantization techniques have been achieved. The quantitative descriptions of
characteristic farfromequilibrium time scales and thermal equilibration in
quantum field theory point out new phenomena such as prethermalization. They
determine the range of validity of standard transport or semiclassical
approaches, on which most of our ideas about nonequilibrium dynamics were based
so far. These are crucial ingredients to understand important topical phenomena
in highenergy physics related to collision experiments of heavy nuclei, early
universe cosmology and complex manybody systems.

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 relaxationtime approximations, where
isotropization and thermal equilibration rates agree. For a weakcoupling
$\phi^4$model the relaxationtime 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 closetoequilibrium value. Our results are obtained from the
twoparticle irreducible effective action, which includes offshell and memory
effects and does not involve a gradient expansion. This allows us to determine
the range of validity of a description to lowestorder in gradients, which is
typically employed in kinetic equations.

Hightemperature resummed perturbation theory is plagued by poor convergence
properties. The problem appears for theories with bosonic field content such as
QCD, QED or scalar theories. We calculate the pressure as well as other
thermodynamic quantities at high temperature for a scalar onecomponent field
theory, solving a threeloop 2PI effective action numerically without further
approximations. We present a detailed comparison with the twoloop
approximation. One observes a strongly improved convergence behavior as
compared to perturbative approaches. The renormalization employed in this work
extends previous prescriptions, and is sufficient to determine all counterterms
required for the theory in the symmetric as well as the spontaneously broken
phase.

The dominance of the direct production of elementary Goldstone waves is
demonstrated in tachyonic preheating by numerically determining the evolution
of the dispersion relation, the equation of state and the kinetic power spectra
for the angular degree of freedom of the complex matter field. The importance
of the domain structure in the order parameter distribution for the
quantitative understanding of the excitation mechanism is emphasized. Evidence
is presented for the very early decoupling of the lowmomentum Goldstone modes.

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 nontrivial order, which includes scattering as well as memory
and offshell effects. The dynamics is solved numerically without further
approximation, for different farfromequilibrium initial conditions. The
latetime 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
FermiDirac and BoseEinstein distributions from the nonequilibrium dynamics.

We study the onset of Goldstone phenomenon in a hybrid inflation scenario.
The physically motivated range of parameters is analyzed in order to meet the
cosmological constraints. Classical equations of motion are solved and the
evolution through the spontaneous symmetry breaking is followed. We emphasize
the role of topological defects that partially maintain the disordered phase
well after the waterfall. We study the emergence of the Goldstone excitations
and their role in the onset of the radiation dominated universe.

The squared mass of a complex scalar field is turned dynamically into
negative by its O(2)invariant coupling to a real field slowly rolling down in
a quadratic potential. The emergence of gapless excitations is studied in real
time simulations after spinodal instability occurs. Careful tests demonstrate
that the Goldstone modes appear almost instantly after the symmetry breaking is
over, much before thermal equilibrium is established.

An increase in the size of coherent domains in the one component $\Phi^4$
field theory under the influence of a uniformly changing external magnetic
field near the critical endpoint $T_{\Phi}=T_c, h_{\Phi}=0$ was proposed
recently as an estimate also for the variation of the chiral correlation length
of QCD near its respective hypothetical end point in the $T_{QCD}\mu_{QCD}$
plane. The present detailed numerical investigation of the effective model
suggests that passing by the critical QCD end point with realistic rate of
temperature change will trigger large amplitude oscillations in the temporal
variation of the chiral correlation length. A simple mechanism for producing
this phenomenon is suggested.

Real time thermalization and relaxation phenomena are studied in the low
energy density phase of the 2+1 dimensional classical O(2) symmetric scalar
theory by solving numerically its dynamics. The nearequilibrium decay rate of
onshell waves and the power law governing the large time asymptotics of the
offshell relaxation agree with the analytic results based on linear response
theory. The realisation of the MerminWagner theorem is also studied in the
final equilibrium ensemble.

Thermalisation of configurations with initial white noise power spectrum is
studied in numerical simulations of a classical onecomponent $\Phi^4$ theory
in 2+1 dimensions, coupled to a small amplitude homogenous external field. The
study is performed for energy densities corresponding to the broken symmetry
phase of the system in equilibrium. The effective equation of the order
parameter motion is reconstructed from its trajectory which starts from an
initial value near the metastable point and ends in the stable ground state.
This phenomenological theory quantitatively accounts for the decay of the false
vacuum. The large amplitude transition of the order parameter between the two
minima displays characteristics reflecting dynamical aspects of the Maxwell
construction.

We have calculated the damping rate $\Gamma ({\bf k})$ for classical
onshell Goldstone modes of the O(2) symmetric scalar fields propagating in a
thermal medium of the broken symmetry phase taking into account the effect of
the explicit symmetry breaking. The result of the oneloop analysis can be
expanded around $\Gamma (0)$, which depends nonanalytically on the parameter
of the explicit symmetry breaking, h. $\Gamma (0)$ vanishes when $h\to 0$,
demonstrating in this way the absence of the restoring force, when the
equilibrium direction of the symmetry breaking is modulated homogeneously.