Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:hep-ph/0106195 · Physical review. D. Particles, fields, gravitation, and cosmology/Physical review. D. Particles and fields

Schwinger-Dyson approach to nonequilibrium classical field theory

Krastan B. Blagoev, Fred Cooper, John F. Dawson, Bogdan Mihaila

Published Nov 26, 2001• 67 citations
arXivHTMLPDF

Abstract

In this paper we discuss a Schwinger-Dyson (SD) approach for determining the time evolution of the unequal time correlation functions of a nonequilibrium classical field theory, where the classical system is described by an initial density matrix at time $t=0.$ We focus on $\ensuremath{\lambda}{\ensuremath{\varphi}}^{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 $\ensuremath{\varphi}$ 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.

DOI
10.1103/physrevd.64.125003

Related papers

Quantum dynamics of phase transitions in broken symmetry λφ4 field theory

Fred Cooper, John F. Dawson, Bogdan Mihaila

We perform a detailed numerical investigation of the dynamics of a single component ``explicitly broken symmetry'' $\ensuremath{\lambda}{\ensuremath{\varphi}}^{4}$ field theory in 1+1 dimensions, using a Schwinger-Dyson equation truncation scheme based on ignoring vertex corrections. In an earlier paper, we called this the bare vertex approximation (BVA). We assume here that the initial state is described by a Gaussian density matrix peaked around some nonzero value of $〈\ensuremath{\varphi}(0)〉,$ and characterized by a single particle Bose-Einstein distribution function at a given temperature. We compute the evolution of the system using three different approximations: the Hartree approximation, the BVA, and a related two-particle irreducible (2PI) $1/N$ expansion, as a function of coupling strength and initial temperature. In the Hartree approximation, the static phase diagram shows that there is a first order phase transition for this system. As we change the initial starting temperature of the system, we find that the BVA relaxes to a new final temperature and exhibits behavior consistent with a second order phase transition. We find that the average fields equilibrate for arbitrary initial conditions in the BVA, unlike the behavior exhibited by the Hartree approximation, and we illustrate how $〈\ensuremath{\varphi}(t)〉$ and $〈\ensuremath{\chi}(t)〉$ depend on the initial temperature and on the coupling constant. The Fourier transform of the two-point functions at late times can be fitted by a Bose-Einstein distribution function whose temperature is independent of momentum. We interpret this as evidence for thermalization.

arXiv:hep-ph/0209051•Physical review. D. Particles, fields, gravitation, and cosmology/Physical review. D. Particles and fields•Mar 28, 2003•68 citations

Classical limit of time-dependent quantum field theory — a Schwinger–Dyson approach

Fred Cooper, Avinash Khare, Harvey A. Rose

We rewrite the Martin-Siggia-Rose (MSR) formalism for the statistical dynamics of classical fields in a covariant second order form appropriate for the statistical dynamics of relativistic field theory. This second order formalism is related to a rotation of Schwinger's closed time path (CTP) formalism for quantum dynamics, with the main difference being that certain vertices are absent in the classical theory. These vertices are higher order in an $\hbar$ expansion. The structure of the second order formulation of the Schwinger Dyson (S-D) equations is identical to that of the rotated CTP formalism apart from initial conditions on the Green's functions and the absence of these vertices. We then discuss self-consistent truncation schemes based on keeping certain graphs in the two-particle irreducible effective action made up of bare vertices and exact Green's functions.

arXiv:hep-ph/0106113•Physics Letters B•Aug 1, 2001•21 citations

Exact and truncated dynamics in nonequilibrium field theory

Gert Aarts A, Gian Franco Bonini, CHRISTOF WETTERICH

Nonperturbative dynamics of quantum fields out of equilibrium is often described by the time evolution of a hierarchy of correlation functions, using approximation methods such as Hartree, large ${N}_{f},$ and $n\mathrm{PI}$-effective action techniques. These truncation schemes can be implemented equally well in a classical statistical system, where results can be tested by comparison with the complete nonlinear evolution obtained by numerical methods. For a $(1+1)$-dimensional scalar field we find that the early-time behavior is reproduced qualitatively by the Hartree dynamics. The inclusion of direct scattering improves this to the quantitative level. We show that the emergence of nonthermal temperature profiles at intermediate times can be understood in terms of the fixed points of the evolution equations in the Hartree approximation. The form of the profile depends explicitly on the initial ensemble. While the truncated evolution equations do not seem to be able to get away from the fixed point, the full nonlinear evolution shows thermalization with a (surprisingly) slow relaxation.

arXiv:hep-ph/0007357•Physical review. D. Particles, fields, gravitation, and cosmology/Physical review. D. Particles and fields•Dec 21, 2000•77 citations

Nonequilibrium evolution of correlation functions: A canonical approach

Supratim Sengupta, F. C. Khanna, Sang Pyo Kim

We study nonequilibrium evolution in a self-interacting quantum field theory invariant under space translation only by using a canonical approach based on the recently developed Liouville--von Neumann formalism. The method is first used to obtain the correlation functions both in and beyond the Hartree approximation, for the quantum mechanical analogue of the ${\ensuremath{\varphi}}^{4}$ model. The technique involves representing the Hamiltonian in a Fock basis of annihilation and creation operators. By separating it into a solvable Gaussian part involving quadratic terms and a perturbation of quartic terms, it is possible to find the improved vacuum state to any desired order. The correlation functions for the field theory are then investigated in the Hartree approximation, and those beyond the Hartree approximation are obtained by finding the improved vacuum state corrected up to $O({\ensuremath{\lambda}}^{2}).$ These correlation functions take into account next-to-leading and next-to-next-to-leading order effects in the coupling constant. We also use the Heisenberg formalism to obtain the time evolution equations for the equal-time, connected correlation functions beyond the leading order. These equations are derived by including the connected four-point functions in the hierarchy. The resulting coupled set of equations form part of an infinite hierarchy of coupled equations relating the various connected n-point functions. The connection with other approaches based on the path integral formalism is established, and the physical implications of the set of equations are discussed with particular emphasis on thermalization.

arXiv:hep-ph/0301071•Physical review. D. Particles, fields, gravitation, and cosmology/Physical review. D. Particles and fields•Nov 26, 2003•7 citations

On Thermalization in Classical Scalar Field Theory

Gert Aarts A, Gian Franco Bonini, CHRISTOF WETTERICH

Thermalization of classical fields is investigated in a \\phi^4 scalar field theory in 1+1 dimensions, discretized on a lattice. We numerically integrate the classical equations of motion using initial conditions sampled from various nonequilibrium probability distributions. Time-dependent expectation values of observables constructed from the canonical momentum are compared with thermal ones. It is found that a closed system, evolving from one initial condition, thermalizes to high precision in the thermodynamic limit, in a time-averaged sense. For ensembles consisting of many members with the same energy, we find that expectation values become stationary - and equal to the thermal values - in the limit of infinitely many members. Initial ensembles with a nonzero (noncanonical) spread in the energy density or other conserved quantities evolve to noncanonical stationary ensembles. In the case of a narrow spread, asymptotic values of primary observables are only mildly affected. In contrast, fluctuations and connected correlation functions will differ substantially from the canonical values. This raises doubts on the use of a straightforward expansion in terms of 1PI-vertex functions to study thermalization.

•Jan 1, 2000•34 citations

Approach to thermalization in the classical φ4 theory in 1+1 dimensions: Energy cascades and universal scaling

D. Boyanovsky, C. Destri, Hector J. de Vega

We study the dynamics of thermalization and the approach to equilibrium in the classical ${\ensuremath{\varphi}}^{4}$ theory in $1+1$ spacetime dimensions. At thermal equilibrium we exploit the equivalence between the classical canonical averages and transfer matrix quantum traces of the anharmonic oscillator to obtain exact results for the temperature dependence of several observables, which provide a set of criteria for thermalization. In this context, comparing to the exact results we find that the Hartree approximation is remarkably accurate in equilibrium. The nonequilibrium dynamics is studied by numerically solving the equations of motion in light-cone coordinates for a broad range of initial conditions and energy densities. The long time evolution is described by several distinct stages, all characterized by a cascade of energy towards the ultraviolet. After an initial transient stage, the spatiotemporal gradient terms become larger than the nonlinear term, and there emerges a stage of universal cascade. This cascade starts at a time scale ${t}_{0}$ independent of the initial conditions (except for very low energy density). During this stage the power spectra feature universal scaling behavior and the front of the cascade $\overline{k}(t)$ moves to the ultraviolet as a power law $\overline{k}(t)\ensuremath{\sim}{t}^{\ensuremath{\alpha}}$ with $\ensuremath{\alpha}\ensuremath{\lesssim}0.25$ an exponent weakly dependent on the energy density alone. The wake behind the cascade is described as a state of Local Thermodynamic Equilibrium (LTE) with all correlations being determined by the equilibrium functional form with an effective time dependent temperature ${T}_{\mathrm{eff}}(t),$ which slowly decreases with time as $\ensuremath{\sim}{t}^{\ensuremath{-}\ensuremath{\alpha}}.$ Two well separated time scales emerge: while ${T}_{\mathrm{eff}}(t)$ varies slowly, the wave vectors in the wake with $k<\overline{k}(t)$ attain LTE on much shorter time scales. This universal scaling stage ends when the front of the ultraviolet cascade reaches the cutoff at a time scale ${t}_{1}\ensuremath{\sim}{a}^{\ensuremath{-}1/\ensuremath{\alpha}}.$ Virialization starts to set much earlier than LTE. We find that strict thermalization is achieved only for an infinite time scale.

arXiv:hep-ph/0306124•Physical review. D. Particles, fields, gravitation, and cosmology/Physical review. D. Particles and fields•Feb 9, 2004•56 citations
In this paper we discuss a Schwinger-Dyson (SD) approach for determining the time evolution of the unequal time correlation functions of a nonequilibrium classical field theory, where the classical system is described by an initial density matrix at time t=0.t=0.t=0. We focus on \ensuremathλ\ensuremathφ4\ensuremath{\lambda}{\ensuremath{\varphi}}^{4}\ensuremathλ\ensuremathφ4 field theory in 1+11+11+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/N1/N1/N expansions. The second approximation is a related approximation made to the SD equations of the original formulation in terms of \ensuremathφ\ensuremath{\varphi}\ensuremathφ 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.
We perform a detailed numerical investigation of the dynamics of a single component ``explicitly broken symmetry'' $\ensuremath{\lambda}{\ensuremath{\varphi}}^{4}$ field theory in 1+1 dimensions, using a Schwinger-Dyson equation truncation scheme based on ignoring vertex corrections. In an earlier paper, we called this the bare vertex approximation (BVA). We assume here that the initial state is described by a Gaussian density matrix peaked around some nonzero value of $〈\ensuremath{\varphi}(0)〉,$ and characterized by a single particle Bose-Einstein distribution function at a given temperature. We compute the evolution of the system using three different approximations: the Hartree approximation, the BVA, and a related two-particle irreducible (2PI) 1/N1/N1/N expansion, as a function of coupling strength and initial temperature. In the Hartree approximation, the static phase diagram shows that there is a first order phase transition for this system. As we change the initial starting temperature of the system, we find that the BVA relaxes to a new final temperature and exhibits behavior consistent with a second order phase transition. We find that the average fields equilibrate for arbitrary initial conditions in the BVA, unlike the behavior exhibited by the Hartree approximation, and we illustrate how $〈\ensuremath{\varphi}(t)〉$ and $〈\ensuremath{\chi}(t)〉$ depend on the initial temperature and on the coupling constant. The Fourier transform of the two-point functions at late times can be fitted by a Bose-Einstein distribution function whose temperature is independent of momentum. We interpret this as evidence for thermalization.
We rewrite the Martin-Siggia-Rose (MSR) formalism for the statistical dynamics of classical fields in a covariant second order form appropriate for the statistical dynamics of relativistic field theory. This second order formalism is related to a rotation of Schwinger's closed time path (CTP) formalism for quantum dynamics, with the main difference being that certain vertices are absent in the classical theory. These vertices are higher order in an ℏ\hbarℏ expansion. The structure of the second order formulation of the Schwinger Dyson (S-D) equations is identical to that of the rotated CTP formalism apart from initial conditions on the Green's functions and the absence of these vertices. We then discuss self-consistent truncation schemes based on keeping certain graphs in the two-particle irreducible effective action made up of bare vertices and exact Green's functions.
Nonperturbative dynamics of quantum fields out of equilibrium is often described by the time evolution of a hierarchy of correlation functions, using approximation methods such as Hartree, large Nf,{N}_{f},Nf​, and nPIn\mathrm{PI}nPI-effective action techniques. These truncation schemes can be implemented equally well in a classical statistical system, where results can be tested by comparison with the complete nonlinear evolution obtained by numerical methods. For a (1+1)(1+1)(1+1)-dimensional scalar field we find that the early-time behavior is reproduced qualitatively by the Hartree dynamics. The inclusion of direct scattering improves this to the quantitative level. We show that the emergence of nonthermal temperature profiles at intermediate times can be understood in terms of the fixed points of the evolution equations in the Hartree approximation. The form of the profile depends explicitly on the initial ensemble. While the truncated evolution equations do not seem to be able to get away from the fixed point, the full nonlinear evolution shows thermalization with a (surprisingly) slow relaxation.
We study nonequilibrium evolution in a self-interacting quantum field theory invariant under space translation only by using a canonical approach based on the recently developed Liouville--von Neumann formalism. The method is first used to obtain the correlation functions both in and beyond the Hartree approximation, for the quantum mechanical analogue of the ${\ensuremath{\varphi}}^{4}$ model. The technique involves representing the Hamiltonian in a Fock basis of annihilation and creation operators. By separating it into a solvable Gaussian part involving quadratic terms and a perturbation of quartic terms, it is possible to find the improved vacuum state to any desired order. The correlation functions for the field theory are then investigated in the Hartree approximation, and those beyond the Hartree approximation are obtained by finding the improved vacuum state corrected up to $O({\ensuremath{\lambda}}^{2}).$ These correlation functions take into account next-to-leading and next-to-next-to-leading order effects in the coupling constant. We also use the Heisenberg formalism to obtain the time evolution equations for the equal-time, connected correlation functions beyond the leading order. These equations are derived by including the connected four-point functions in the hierarchy. The resulting coupled set of equations form part of an infinite hierarchy of coupled equations relating the various connected n-point functions. The connection with other approaches based on the path integral formalism is established, and the physical implications of the set of equations are discussed with particular emphasis on thermalization.
Thermalization of classical fields is investigated in a \\phi^4 scalar field theory in 1+1 dimensions, discretized on a lattice. We numerically integrate the classical equations of motion using initial conditions sampled from various nonequilibrium probability distributions. Time-dependent expectation values of observables constructed from the canonical momentum are compared with thermal ones. It is found that a closed system, evolving from one initial condition, thermalizes to high precision in the thermodynamic limit, in a time-averaged sense. For ensembles consisting of many members with the same energy, we find that expectation values become stationary - and equal to the thermal values - in the limit of infinitely many members. Initial ensembles with a nonzero (noncanonical) spread in the energy density or other conserved quantities evolve to noncanonical stationary ensembles. In the case of a narrow spread, asymptotic values of primary observables are only mildly affected. In contrast, fluctuations and connected correlation functions will differ substantially from the canonical values. This raises doubts on the use of a straightforward expansion in terms of 1PI-vertex functions to study thermalization.
We study the dynamics of thermalization and the approach to equilibrium in the classical ${\ensuremath{\varphi}}^{4}$ theory in 1+11+11+1 spacetime dimensions. At thermal equilibrium we exploit the equivalence between the classical canonical averages and transfer matrix quantum traces of the anharmonic oscillator to obtain exact results for the temperature dependence of several observables, which provide a set of criteria for thermalization. In this context, comparing to the exact results we find that the Hartree approximation is remarkably accurate in equilibrium. The nonequilibrium dynamics is studied by numerically solving the equations of motion in light-cone coordinates for a broad range of initial conditions and energy densities. The long time evolution is described by several distinct stages, all characterized by a cascade of energy towards the ultraviolet. After an initial transient stage, the spatiotemporal gradient terms become larger than the nonlinear term, and there emerges a stage of universal cascade. This cascade starts at a time scale t0{t}_{0}t0​ independent of the initial conditions (except for very low energy density). During this stage the power spectra feature universal scaling behavior and the front of the cascade k‾(t)\overline{k}(t)k(t) moves to the ultraviolet as a power law $\overline{k}(t)\ensuremath{\sim}{t}^{\ensuremath{\alpha}}$ with $\ensuremath{\alpha}\ensuremath{\lesssim}0.25$ an exponent weakly dependent on the energy density alone. The wake behind the cascade is described as a state of Local Thermodynamic Equilibrium (LTE) with all correlations being determined by the equilibrium functional form with an effective time dependent temperature Teff(t),{T}_{\mathrm{eff}}(t),Teff​(t), which slowly decreases with time as $\ensuremath{\sim}{t}^{\ensuremath{-}\ensuremath{\alpha}}.$ Two well separated time scales emerge: while Teff(t){T}_{\mathrm{eff}}(t)Teff​(t) varies slowly, the wave vectors in the wake with k<k‾(t)k<\overline{k}(t)k<k(t) attain LTE on much shorter time scales. This universal scaling stage ends when the front of the ultraviolet cascade reaches the cutoff at a time scale ${t}_{1}\ensuremath{\sim}{a}^{\ensuremath{-}1/\ensuremath{\alpha}}.$ Virialization starts to set much earlier than LTE. We find that strict thermalization is achieved only for an infinite time scale.