• We compute the slope and curvature, at vanishing four-momentum 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 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.
  • 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 quantum-chromodynamics and quantum-electrodynamics computations with four nondegenerate Wilson fermion flavors and computed the neutron-proton mass-splitting with an accuracy of $300$ kilo-electron 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 up-down 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 up-down 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 up-down 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 cut-off 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 Stefan-Boltzmann 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 cut-off 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 cross-over 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 hep-ph/0302210 and hep-ph/0410330.
  • We investigate lattice simulations of scalar and nonabelian gauge fields in Minkowski space-time. For SU(2) gauge-theory expectation values of link variables in 3+1 dimensions are constructed by a stochastic process in an additional (5th) ``Langevin-time''. A sufficiently small Langevin step size and the use of a tilted real-time contour leads to converging results in general. All fixed point solutions are shown to fulfil the infinite hierarchy of Dyson-Schwinger identities, however, they are not unique without further constraints. For the nonabelian gauge theory the thermal equilibrium fixed point is only approached at intermediate Langevin-times. It becomes more stable if the complex time path is deformed towards Euclidean space-time. We analyze this behavior further using the real-time 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 real-time 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 far-from-equilibrium 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 semi-classical approaches, on which most of our ideas about nonequilibrium dynamics were based so far. These are crucial ingredients to understand important topical phenomena in high-energy physics related to collision experiments of heavy nuclei, early universe cosmology and complex many-body 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 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.
  • High-temperature 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 one-component field theory, solving a three-loop 2PI effective action numerically without further approximations. We present a detailed comparison with the two-loop 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 low-momentum 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 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.
  • 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 end-point $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 near-equilibrium decay rate of on-shell waves and the power law governing the large time asymptotics of the off-shell relaxation agree with the analytic results based on linear response theory. The realisation of the Mermin-Wagner 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 one-component $\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 on-shell 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 one-loop analysis can be expanded around $\Gamma (0)$, which depends non-analytically 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.