research42 min studyLesson 748 in Physics

Derive a kinetic approximation from two-time dynamics

Separate exact memory evolution from Wigner and quasiparticle approximations, track a first gradient sign, and identify which assumptions make a Boltzmann equation meaningful.

On this page

What you’ll understand

  • Explain what two-time Kadanoff–Baym evolution retains beyond a kinetic equation.
  • Use a specified Wigner convention to compute a first gradient correction.
  • Identify memory, quasiparticle, and conservation assumptions behind relaxation-time models.

A field's correlation function depends on two times, retaining memory of initial conditions and past interactions. A kinetic distribution depending on one time and momentum is a reduced description of that information. The Kadanoff–Baym framework organizes the two-time equations; Wigner transforms, gradient expansions, and quasiparticle assumptions then connect them to transport. These reductions are distinct approximations and should not be treated as interchangeable names for an exact Boltzmann equation.

Retain the initial-value and memory structure

The contour Dyson equation yields coupled equations for statistical correlations and spectral response. Their self-energy integrals run over earlier times and encode memory, scattering, and off-shell structure. Initial correlations supply boundary information. If the self-energy were exact, the resulting equations would be exact too; practical calculations require a truncation or closure. A local thermal occupation should not be imposed on both functions before examining whether their evolution actually supports an equilibrium fluctuation relation.

A two-particle-irreducible effective-action truncation can organize self-consistent self-energies and preserve important conservation properties. That is useful, but conserving an energy or charge does not guarantee all spectral features or long-time relaxation are quantitatively correct. Leading mean-field closures can retain artificial conserved occupations and miss scattering. Numerical memory cutoffs create another approximation. State the initial state, retained self-energy terms, renormalization prescription, and time range before interpreting an apparent late-time plateau as a prediction of the full theory.

X=x+y2,s=x−y,AW(X,p)=∫ds e−ipsA(X+s/2,X−s/2)X=\frac{x+y}{2},\qquad s=x-y,\qquad A_W(X,p)=\int ds\,e^{-ips}A(X+s/2,X-s/2)
This displayed Wigner transform is for one spatial coordinate and its momentum, using the module's negative spatial Fourier sign. X and s have inverse-energy units, while p has energy units. A full spacetime transform uses the separately declared +ωt−k·x convention.↔ Scroll wide equations horizontally.

Expand convolutions only when center variation is slow

In Wigner variables, a convolution becomes a star product. For the one-dimensional spatial convention above, its first correction is plus i over two times the Poisson bracket ∂XA ∂pB−∂pA ∂XB\partial_XA\,\partial_pB-\partial_pA\,\partial_XB. Reversing the Fourier sign would reverse that correction, so a remembered formula must be matched to the transform. The expansion is controlled only when center-coordinate variation is slow relative to the relevant microscopic correlation scales. Near an abrupt quench, that condition can fail even if the interaction coupling is weak.

A gradient expansion alone does not put particles on shell. A quasiparticle description additionally assumes a sufficiently identifiable spectral branch and a justified reduction of its energy dependence. Suppressing memory to obtain a local collision term needs another scale separation. Strong broadening, coherent production, or nonlocal formation processes can obstruct those reductions. A useful kinetic equation should therefore be derived with its collision channels and validity domain, rather than introduced by fitting an instantaneous temperature to an evolving two-time function.

A⋆B=AB+i2{A,B}X,p+O(∂X2),{A,B}X,p=∂XA ∂pB−∂pA ∂XBA\star B=AB+\frac{i}{2}\{A,B\}_{X,p}+O(\partial_X^2),\qquad \{A,B\}_{X,p}=\partial_XA\,\partial_pB-\partial_pA\,\partial_XB
This is the first-gradient expansion of a spatial convolution with the displayed e^(−ips) Wigner convention. Higher terms contain additional paired derivatives and are small only with a justified gradient hierarchy. A complex star product is an operator-composition quantity, not directly a probability distribution.↔ Scroll wide equations horizontally.

Check a reduced collision model

Under suitable on-shell and slow-variation assumptions, a kinetic equation has a drift term and a collision operator. Bose enhancement and Pauli blocking enter that operator with factors 1+f1+f and 1−f1-f, respectively. Conservation laws require the appropriate momentum integrals of collisions to vanish for conserved quantities. A simple homogeneous relaxation ansatz f˙=−(f−feq)/τR\dot f=-(f-f_{\mathrm{eq}})/\tau_R can illustrate a decay scale, but a fixed target distribution need not conserve the evolving system's energy or particle number unless it has been chosen self-consistently.

Comparing a microscopic correlation time with a macroscopic relaxation time is one useful check, not a full proof of kinetic validity. The spectrum, occupations, gradients, and interaction expansion must also be controlled. In weakly coupled gauge theories, repeated soft interactions and collinear formation can require effective collision structures beyond a naive instantaneous two-to-two term. A numerical agreement over a short interval should be tested against a two-time calculation or other constraints when possible. The goal is to identify what information the kinetic approximation discards and whether that information matters to the requested observable.

A relaxation ansatz and one scale-separation check

A homogeneous illustrative occupation obeys f˙=−(f−0.2)/τR\dot f=-(f-0.2)/\tau_R with f(0)=0.8f(0)=0.8 and τR=2 GeV−1\tau_R=2\,\mathrm{GeV^{-1}}. Find f at t=2log⁡2 GeV−1t=2\log2\,\mathrm{GeV^{-1}}. If a correlation time is 0.2 GeV−10.2\,\mathrm{GeV^{-1}}, evaluate its ratio to the relaxation time.

  1. Solve the linear initial-value problem: f(t)=0.2+(0.8−0.2)e−t/τR=0.2+0.6e−t/2f(t)=0.2+(0.8-0.2)e^{-t/\tau_R}=0.2+0.6e^{-t/2}. It gives the specified initial occupation and tends to the chosen target at late time.

  2. At t=2log⁡2t=2\log2, the exponential is one half, so f=0.2+0.3=0.5f=0.2+0.3=0.5. The time argument is dimensionless because both t and the relaxation time use inverse-GeV units.

  3. The correlation-to-relaxation ratio is 0.2/2=0.10.2/2=0.1. This suggests a potentially useful memory hierarchy, but it does not check spectral narrowness, spatial gradients, or whether the chosen relaxation ansatz conserves energy.

The result

The ansatz predicts occupation one half and a time-scale ratio of one tenth. It is an explicitly reduced model, not a solution of the full Kadanoff–Baym field equations.

Explore the model

Interactive lab

Interactive lab · markov

Inspect a deterministic approach to a stationary law

Reset and set transition p = 0.2 and q = 0.3. Starting from state 0, predict probability of state 1 equal to 0.2 after one Step and 0.3 after two Steps; the stationary value is 0.4. The lab computes exact finite-step probabilities, not sampled trajectories. It has no quantum memory kernel, energy spectrum, Wigner transform, or field self-energy. Reset and set both controls before any fresh comparison.

Inspect a deterministic approach to a stationary lawExact finite-step probabilities for a two-state chain starting in state 0, not sampled trajectories. Playback advances two transitions per second. When p = q = 0 every distribution is stationary; when p = q = 1 the chain alternates forever. A single zero transition creates an absorbing state. Current readings: P(state 1) at step 0: 0%; Stationary P(state 1): 60%; Chain behavior: Monotone convergence. All parameters are adjustable with the labeled controls below.Two-state Markov distribution after 0 steps: state zero probability 1, state one probability 0.Initial distribution (1, 0) · step 0p = 0.3q = 0.201100%0%1−p1−qP(0)P(1)Monotone convergence
P(state 0)P(state 1)Stationary distribution
P(state 1) at step 0
0%
Stationary P(state 1)
60%
Chain behavior
Monotone convergence
t = 0 s

Exact finite-step probabilities for a two-state chain starting in state 0, not sampled trajectories. Playback advances two transitions per second. When p = q = 0 every distribution is stationary; when p = q = 1 the chain alternates forever. A single zero transition creates an absorbing state. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forA simple relaxation law can be inspected exactly while its discarded microscopic information stays explicit.

Check your understanding

Which additional step is required beyond a Wigner transform to justify a local quasiparticle Boltzmann description?

Try it yourself

For the displayed spatial star product take A=p2/(2m)A=p^2/(2m) with m=4 GeVm=4\,\mathrm{GeV}. At one phase-space point let p=1 GeVp=1\,\mathrm{GeV}, B=2B=2, and ∂XB=0.4 GeV\partial_XB=0.4\,\mathrm{GeV}. Compute AB and the first-gradient correction, given that A has no X dependence.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Only −(∂pA)(∂XB) survives in the Poisson bracket. Retain the explicit i/2 factor and its sign.

Compare with the solution

At this point A=1/8 GeVA=1/8\,\mathrm{GeV}, so AB=1/4 GeVAB=1/4\,\mathrm{GeV}. Its momentum derivative is p/m=1/4p/m=1/4, giving {A,B}=−(1/4)(0.4)=−0.1 GeV\{A,B\}=-(1/4)(0.4)=-0.1\,\mathrm{GeV}. Therefore A⋆B≃(0.25−0.05i) GeVA\star B\simeq(0.25-0.05i)\,\mathrm{GeV} at first gradient order. The correction's magnitude is twenty percent of the leading product, so higher derivatives deserve examination rather than an automatic claim of excellent convergence. This one-dimensional composition example does not compute a field occupation or a full collision term.

Take it with you

  • Two-time equations retain memory and spectral information absent from a one-time occupation.
  • A Wigner transform does not itself justify gradient, on-shell, or Markovian reductions.
  • Check gradient signs, scale hierarchies, and conservation when using a kinetic closure.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.