Skip to content

Particle Creation in Time-Dependent Backgrounds

Time-dependent backgrounds create an in/out particle interpretation by mixing asymptotic positive and negative frequencies. The Bogoliubov coefficient βk\beta_{\mathbf k} gives out occupation in the in-state when both asymptotic bases exist and are normalized. It does not by itself give a local energy density during the evolution, and a different out complex structure defines a different number observable.

Required background. Bogoliubov Transformations and Unitary Implementability supplies the representation criterion; Mode Bases, In/Out States, and Number Operators supplies asymptotic normalization.

Helpful background. In–Out versus In–In Expectation Values separates amplitudes from real-time expectation values; Adiabatic States, WKB Order, and Regularity supplies controlled approximate bases.

For a homogeneous scalar background, each canonical mode can often be reduced to

vk+ωk(η)2vk=0,vkvkvkvk=i.v_{\mathbf k}'' +\omega_{\mathbf k}(\eta)^2v_{\mathbf k}=0, \qquad v_{\mathbf k}v_{\mathbf k}^{*'} -v_{\mathbf k}'v_{\mathbf k}^*=i.

Assume ωkωkin/out>0\omega_{\mathbf k}\to\omega_{\mathbf k}^{\mathrm{in/out}}>0 in stationary asymptotic regions. A normalized in-mode has future form

vkinαkeiωkoutη+βke+iωkoutη2ωkout.v_{\mathbf k}^{\mathrm{in}} \longrightarrow \frac{ \alpha_{\mathbf k}e^{-i\omega_{\mathbf k}^{\mathrm{out}}\eta} +\beta_{\mathbf k}e^{+i\omega_{\mathbf k}^{\mathrm{out}}\eta} }{ \sqrt{2\omega_{\mathbf k}^{\mathrm{out}}} }.

Wronskian conservation yields

αk2βk2=1,nkout=βk2\lvert\alpha_{\mathbf k}\rvert^2 -\lvert\beta_{\mathbf k}\rvert^2=1, \qquad n_{\mathbf k}^{\mathrm{out}} = \lvert\beta_{\mathbf k}\rvert^2

in the in-vacuum. The background supplies the energy. In a prescribed geometry or external field, backreaction is neglected only while the produced stress is small relative to the source that drives ωk\omega_{\mathbf k}.

The late-time excess out energy for a diagonal free scalar state is

Δρout=dd1k(2π)d1ωkoutβk2,\Delta\rho_{\mathrm{out}} = \int\frac{\mathrm d^{d-1}k}{(2\pi)^{d-1}}\, \omega_{\mathbf k}^{\mathrm{out}} \lvert\beta_{\mathbf k}\rvert^2,

after the out-vacuum contribution is treated in the chosen renormalization scheme. During the quench, a local Tμνren\langle T_{\mu\nu}\rangle_{\mathrm{ren}} can also contain phase-sensitive interference terms and local vacuum-polarization terms. Number and local energy are therefore checked separately.

First application: an exactly smooth frequency quench

Section titled “First application: an exactly smooth frequency quench”

Take

ωk2(t)=(ωkin)2+(ωkout)22+(ωkout)2(ωkin)22tanhtτ.\omega_{\mathbf k}^2(t) = \frac{ (\omega_{\mathbf k}^{\mathrm{in}})^2 +(\omega_{\mathbf k}^{\mathrm{out}})^2 }{2} +\frac{ (\omega_{\mathbf k}^{\mathrm{out}})^2 -(\omega_{\mathbf k}^{\mathrm{in}})^2 }{2} \tanh\frac{t}{\tau}.

Hypergeometric connection formulas give

βk2=sinh2 ⁣[πτ2(ωkoutωkin)]sinh(πτωkin)sinh(πτωkout).\lvert\beta_{\mathbf k}\rvert^2 = \frac{ \sinh^2\!\left[ \frac{\pi\tau}{2} \left( \omega_{\mathbf k}^{\mathrm{out}} -\omega_{\mathbf k}^{\mathrm{in}} \right) \right] }{ \sinh(\pi\tau\omega_{\mathbf k}^{\mathrm{in}}) \sinh(\pi\tau\omega_{\mathbf k}^{\mathrm{out}}) }.

Das, Galante, and Myers give the exact modes, coefficients, and distinction between smooth-fast and instantaneous quenches Das, Galante, and Myers 2015, §2.

Two limits check the result. For τ0\tau\to0 at fixed mode,

βk2(ωkoutωkin)24ωkinωkout.\lvert\beta_{\mathbf k}\rvert^2 \longrightarrow \frac{ \left( \omega_{\mathbf k}^{\mathrm{out}} -\omega_{\mathbf k}^{\mathrm{in}} \right)^2 }{ 4\omega_{\mathbf k}^{\mathrm{in}} \omega_{\mathbf k}^{\mathrm{out}} }.

For a slow quench with both asymptotic frequencies nonzero, production is exponentially suppressed. At large momentum the finite-τ\tau result also falls exponentially, while an ideal instantaneous quench can have a much harder ultraviolet tail. The order of the τ0\tau\to0 and cutoff limits is physical.

The number and energy densities are

nout=dd1k(2π)d1βk2,Δρout=dd1k(2π)d1ωkoutβk2.n_{\mathrm{out}} = \int\frac{\mathrm d^{d-1}k}{(2\pi)^{d-1}} \lvert\beta_{\mathbf k}\rvert^2, \qquad \Delta\rho_{\mathrm{out}} = \int\frac{\mathrm d^{d-1}k}{(2\pi)^{d-1}} \omega_{\mathbf k}^{\mathrm{out}} \lvert\beta_{\mathbf k}\rvert^2.

Convergence of the first does not imply convergence of the second because the energy integral has an extra power of frequency.

Keep the evolved state fixed and define a second out basis

u~kout=Akukout+Bkukout,Ak2Bk2=1.\widetilde u_{\mathbf k}^{\mathrm{out}} = A_{\mathbf k}u_{\mathbf k}^{\mathrm{out}} +B_{\mathbf k}u_{-\mathbf k}^{\mathrm{out}*}, \qquad \lvert A_{\mathbf k}\rvert^2-\lvert B_{\mathbf k}\rvert^2=1.

The new production coefficient is a composition of the physical evolution with (A,B)(A,B), so β~k2\lvert\widetilde\beta_{\mathbf k}\rvert^2 generally differs. If the original out basis was selected by a future Killing generator and the rotated one is not, the latter number is still a mathematically valid basis occupation but has lost the asymptotic energy interpretation.

Local stress-energy, computed from the fixed two-point function with one fixed renormalization prescription, does not change under this relabeling. This is the decisive comparison: basis-dependent number changes, whereas the same local observable does not.

The construction map places Bogoliubov evolution after asymptotic basis selection and before interpretation as number or energy.

A smooth time-dependent frequency connects normalized in and out modes, producing Bogoliubov occupation that is then checked against out energy

In/out production is controlled by asymptotic normalization, canonical identities, ultraviolet convergence, and background energy balance; the map is schematic and not to scale.

The failure map identifies an unselected instantaneous or rotated basis as a failure witness, not as new invariant production.

A production-number claim is downgraded when asymptotic stationarity, basis selection, ultraviolet convergence, or source energy accounting fails

Only the number and energy associated with the declared asymptotic basis and controlled quench are licensed; the map is schematic and not to scale.

Use Domain and failure conditions. This page’s checks are Wronskian normalization, two asymptotic regimes, α2β2=1\lvert\alpha\rvert^2-\lvert\beta\rvert^2=1, wavepacket or density normalization, ultraviolet convergence, local-energy translation, and negligible backreaction.

Can a quench have finite number density but divergent excess energy in an idealized sudden limit?

Solution

Yes. The energy integral weights each produced mode by ωkout\omega_{\mathbf k}^{\mathrm{out}}. A high-momentum tail can therefore be integrable for nk\int n_{\mathbf k} but not for ωknk\int\omega_{\mathbf k}n_{\mathbf k}. A finite switching time or physical cutoff must remain when the instantaneous idealization exceeds its ultraviolet domain.

Adiabaticity, Stokes Phenomena, and Production Rates analyzes the complex-time origin of β\beta. Chapter 2 owns representation selection, Chapter 6 collapse-induced production, and Chapter 14 cosmological deployment.

  • N. D. Birrell and P. C. W. Davies, Quantum Fields in Curved Space, Cambridge University Press (1982), DOI, Chapter 3.
  • Sumit R. Das, Damián A. Galante, and Robert C. Myers, “Smooth and Fast versus Instantaneous Quenches in Quantum Field Theory,” Journal of High Energy Physics 2015 (2015), article 73, DOI, arXiv:1505.05224.
  • Leonard Parker, “Particle Creation in Expanding Universes,” Physical Review Letters 21 (1968), 562–564, DOI.