Schwinger-Dyson approach to non-equilibrium classical field theory

In this paper we discuss a Schwinger-Dyson [SD] approach for determining the time evolution of the unequal time correlation functions of a non-equilibrium classical field theory, where the classical system is described by an initial density matrix at time $t=0$. We focus on $\lambda \phi^4$ field theory in 1+1 space time dimensions where we can perform exact numerical simulations by sampling an ensemble of initial conditions specified by the initial density matrix. We discuss two approaches. The first, the bare vertex approximation [BVA], is based on ignoring vertex corrections to the SD equations in the auxiliary field formalism relevant for 1/N expansions. The second approximation is a related approximation made to the SD equations of the original formulation in terms of $\phi$ alone. We compare these SD approximations as well as a Hartree approximation with exact numerical simulations. We find that both approximations based on the SD equations yield good agreement with exact numerical simulations and cure the late time oscillation problem of the Hartree approximation. We also discuss the relationship between the quantum and classical SD equations.Comment: 36 pages, 5 figure

Turbulent Thermalization

We study, analytically and with lattice simulations, the decay of coherent field oscillations and the subsequent thermalization of the resulting stochastic classical wave-field. The problem of reheating of the Universe after inflation constitutes our prime motivation and application of the results. We identify three different stages of these processes. During the initial stage of ``parametric resonance'', only a small fraction of the initial inflaton energy is transferred to fluctuations in the physically relevant case of sufficiently large couplings. A major fraction is transfered in the prompt regime of driven turbulence. The subsequent long stage of thermalization classifies as free turbulence. During the turbulent stages, the evolution of particle distribution functions is self-similar. We show that wave kinetic theory successfully describes the late stages of our lattice calculation. Our analytical results are general and give estimates of reheating time and temperature in terms of coupling constants and initial inflaton amplitude.Comment: 27 pages, 13 figure

The matter power spectrum in redshift space using effective field theory

The use of Eulerian 'standard perturbation theory' to describe mass assembly in the early universe has traditionally been limited to modes with k <= 0.1 h/Mpc at z=0. At larger k the SPT power spectrum deviates from measurements made using N-body simulations. Recently, there has been progress in extending the reach of perturbation theory to larger k using ideas borrowed from effective field theory. We revisit the computation of the redshift-space matter power spectrum within this framework, including for the first time for the full one-loop time dependence. We use a resummation scheme proposed by Vlah et al. to account for damping of the baryonic acoustic oscillations due to large-scale random motions and show that this has a significant effect on the multipole power spectra. We renormalize by comparison to a suite of custom N-body simulations matching the MultiDark MDR1 cosmology. At z=0 and for scales k <~ 0.4 h/Mpc we find that the EFT furnishes a description of the real-space power spectrum up to ~ 2%, for the ell=0 mode up to ~ 5% and for the ell = 2, 4 modes up to ~ 25%. We argue that, in the MDR1 cosmology, positivity of the ell = 0 mode gives a firm upper limit of k ~ 0.74 h/Mpc for the validity of the one-loop EFT prediction in redshift space using only the lowest-order counterterm. We show that replacing the one-loop growth factors by their Einstein-de Sitter counterparts is a good approximation for the ell = 0 mode, but can induce deviations as large as 2% for the ell = 2, 4 modes. An accompanying software bundle, distributed under open source licenses, includes Mathematica notebooks describing the calculation, together with parallel pipelines capable of computing both the necessary one-loop SPT integrals and the effective field theory counterterms