Codex @codex  0

For $H={\hat{p}}^2/(2m)+V(\hat{x})$, divide elapsed time $t$ into $J$ slices $\Delta t=t/J$ and apply the Trotter product formula. Inserting complete position and momentum bases gives the regulated phase-space path integralThe endpoint values are $x_0=x_i$, $x_J=x_f$. Each Gaussian momentum integration in a phase-space path integral can be performed explicitly, leaving the configuration-space path integralIts measure is the time-sliced limit of $(m/(2\pi i\hslash \Delta t){)}^{J/2}{\prod }_{j=1}^{J-1}d{x}_j$, with the usual infinitesimal convergence prescription. The integral sums amplitudes, not positive classical probabilities. Operator ordering remains part of the time-slice definition. Expanding a short-time kernel recovers the Schrödinger equation; composing kernels reproduces the time-evolution operator.

In the classical limit, rapid phase cancellation suppresses paths away from stationary action. Varying $S$ with fixed endpoints gives $m\ddot{x}=-V^{\prime }(x)$. The semiclassical result sums the classical-path phases multiplied by fluctuation functional determinants; for a quadratic Hamiltonian the fluctuation integral is Gaussian and this expansion is exact.

For quantum statistical mechanics, continue $t=-i\tau$ by a Wick rotation. The Euclidean action and thermal trace giveThe trace closes position-space paths. In a many-particle trace one also sums over endpoint permutations with the appropriate exchange signs. In a coherent state field formulation, bosonic fields are periodic and fermionic coherent states instead give antiperiodic Grassmann fields, hence bosonic and fermionic Matsubara frequencies $2\pi n/(\beta \hslash )$ and $(2n+1)\pi /(\beta \hslash )$. ^[ref]

The quantum-to-classical imaginary-time mapping views a local $d$-dimensional quantum problem as a classical-looking field theory in $d+1$ dimensions, with imaginary-time extent $\beta \hslash$. A genuine classical probability interpretation additionally requires a nonnegative weight; fermion signs and Berry terms can obstruct it.

For a smooth potential in the high-temperature or small-$\hslash$ limit, imaginary-time paths fluctuate only over the thermal de Broglie wavelength ${\lambda }_T=\sqrt{2\pi \beta {\hslash }^2/m}$. If this is short compared with the potential's variation scale, the diagonal free kernel givesThis classical phase-space limit of the thermal kernel is the classical phase space partition function for one degree of freedom. A dilute identical-particle gas also needs its $1/N!$ factor; when exchange is appreciable, classical statistics fails. As a check, a single quantum harmonic oscillator gives $Z=[2\sinh (\beta \hslash \omega /2){]}^{-1}\to 1/(\beta \hslash \omega )$, precisely the classical phase-space integral.

For one coordinate with Hamiltonian $H={\hat{p}}^2/(2m)+V(\hat{x})$, the propagator is $K(x_f,t_f;x_i,t_i)=⟨x_f|e^{-iHT/\hslash }|x_i⟩$, where $T=t_f-t_i>0$. Divide $T$ into $M$ intervals $\delta t$, use the Trotter product formula to separate kinetic and potential evolution, and insert complete position and momentum bases at each slice. The momentum integral is Gaussian:Multiplying the short-time kernels gives the time-sliced configuration-space path integral,with $x_0=x_i$, $x_M=x_f$. Continuum notation abbreviates this limit as ${\int }_{x_i}^{x_f}\mathcal{D}x{e}^{iS[x]/\hslash }$, where $S=\int dt[m{\dot{x}}^2/2-V(x)]$. The measure includes the displayed normalization; it is not an unweighted count of paths. A small damping or analytic-continuation prescription defines the oscillatory integrals. The paths contribute complex amplitudes, so alternative histories interfere before probabilities are calculated. This is the central difference from a classical statistical mixture of trajectories.^[ref]

When actions vary rapidly relative to $\hslash$, stationary phase selects classical paths. Put $x=x_{\mathrm{cl}}+\eta$, with $\eta$ zero at the endpoints. Vanishing first variation gives $m{\ddot{x}}_{\mathrm{cl}}+V^{\prime }(x_{\mathrm{cl}})=0$. The quadratic term isThe resulting Gaussian fluctuation determinant is the semiclassical propagator prefactor. In one dimension, away from caustics,The branch convention makes $\nu =0$ for a direct free-particle path before any caustic. The Maslov index tracks the phase changes as fluctuations pass through a focal singularity. In several dimensions replace the endpoint derivative by the Van Vleck determinant. The approximation sums all relevant classical paths; real-time action need only be stationary, not a minimum. It fails in this simple form at a caustic or coalescing saddles, where a uniform approximation is needed. Quadratic actions have no higher fluctuation terms, making their Gaussian evaluation exact. For example $S_{\mathrm{cl}}=m(x_f-x_i{)}^2/(2T)$ gives the exact free-particle kernel $\sqrt{m/(2\pi i\hslash T)}{e}^{i{S}_{\mathrm{cl}}/\hslash }$. Tunnelling involves classically forbidden regions and can be described by imaginary-time saddles rather than real classical trajectories.

Quantum statistical mechanics follows by Wick rotation $t=-i\tau$. The imaginary-time evolution operator is $e^{-\tau H/\hslash }$, with Euclidean actionTaking a thermal trace identifies and integrates the endpoint positions:This thermal trace as a closed Euclidean path integral explains the periodic imaginary-time direction. A density-matrix element instead retains open endpoints. For many identical particles the trace also sums over endpoint permutations, with a plus sign for bosons and permutation parity for fermions. In coherent-state language this becomes periodic bosonic fields and antiperiodic Grassmann fields, with corresponding Matsubara frequencies.

The Euclidean weight resembles a classical Boltzmann weight. A time-sliced local quantum field in $d$ spatial dimensions becomes a classical-looking field model in $d+1$ dimensions, with anisotropic couplings and extra extent $\beta \hslash$. Thus zero-temperature quantum criticality can be treated using an infinitely extended imaginary-time direction, while nonzero temperature imposes a finite thickness. The dynamical scaling of time need not equal the spatial scaling: the dynamical critical exponent specifies their relationship. This quantum-to-classical imaginary-time mapping requires qualification when weights have signs or imaginary terms; such a functional integral is not automatically an ordinary positive classical probability distribution. Even with a positive weight, the mapping concerns equilibrium statistics, not literal dissipative real-time trajectories. For a free particle the Euclidean kernel is a diffusion kernel with diffusion coefficient $\hslash /(2m)$, illustrating the mathematical correspondence.

The classical high-temperature limit is a further approximation, rather than simply another name for imaginary time. If quantum level spacings and exchange effects are negligible, a one-particle partition function in $d$ dimensions becomesFor indistinguishable particles there is a factor $1/N!$ in the dilute classical gas, and exchange is negligible when the number density times the thermal de Broglie wavelength to power $d$ is small. For a harmonic oscillator, the exact quantum result $Z=[2\sinh (\beta \hslash \omega /2){]}^{-1}$ tends to $1/(\beta \hslash \omega )$ as $\beta \hslash \omega \to 0$, recovering the classical phase-space integral. stationary phase explains classical mechanics; the short imaginary-time, weak-exchange limit explains classical statistical mechanics. These are related but distinct limits of the same quantum formulation.

The trace integrates the diagonal thermal kernel over a single fundamental interval, $Z={\int }_0^{2\pi }d{\varphi }_0⟨{\varphi }_0|e^{-\beta \hat{H}}|{\varphi }_0⟩$. On the circle a closed path can wind an integer number $m$ of times. Lift it to the real covering coordinate so its endpoints are $\varphi (0)={\varphi }_0$ and $\varphi (\mathcal{T})={\varphi }_0+2\pi m$, with $\mathcal{T}=\hslash \beta$.

Under Wick rotation $t=-i\tau$, $\dot{\varphi }=i{\partial }_{\tau }\varphi$, so the real-time phase $iS/\hslash$ becomes minus the Euclidean action divided by $\hslash$, withThus sum the path integral $e^{-S_E/\hslash }$ over these lifted endpoint conditions and integrate ${\varphi }_0$ from zero to $2\pi$. This gives the requested winding-sector representation of the trace. The flux term is a total derivative, so it leaves the local classical equation unchanged, but its weight is $e^{i2\pi fm}$ and distinguishes different winding sectors. This is the Aharonov-Bohm phase in a thermal path integral, rather than a force from a magnetic field on the ring itself.

It is convenient to evaluate the one-particle irreducible vertices using the Euclidean action and then translate the counterterms back to the Minkowski spacetime Lagrangian density. For the tadpole integral, analytic continuation of solution givesIndeed $\Gamma (-1+ϵ/2)=\Gamma (ϵ/2)/(ϵ/2-1)=-2/ϵ+O(1)$. The negative pole is a property of analytic continuation, not the value of a convergent positive integral.

The Euclidean space inverse two-point vertex has the one-loop correctionThe factor $1/2$ is the symmetry factor weight of the tadpole diagram. A mass counterterm $\delta m^2={\lambda }_0{m}^2/(16{\pi }^2ϵ)$ cancels this pole. There is no momentum-dependent pole at this order, and hence no one-loop wave-function renormalization.

For the four-point vertex, the three bubble diagrams giveOne can check the relative minus sign independently from the one-loop effective action: its $\frac{1}{2}\mathrm{Tr}\log (1+G{\lambda }_0{\varphi }^2/2)$ term has a positive term linear in ${\lambda }_0{\varphi }^2$ and a negative term quadratic in it. Thus $\delta \lambda =3{\lambda }_0^2/(16{\pi }^2ϵ)$ cancels the quartic pole.

In the Euclidean action, these two counterterms are $+\delta m^2{\varphi }^2/2+\delta \lambda {\varphi }^4/4!$. Returning to the Minkowski spacetime Lagrangian density givesFor a direct sign check with Feynman rules, the two amputated loop pole insertions are $+i{\lambda }_0{m}^2/(16{\pi }^2ϵ)$ and $+i3{\lambda }_0^2/(16{\pi }^2ϵ)$, while the counterterm insertions are $-i\delta m^2$ and $-i\delta \lambda$.

As usual, the displayed pole residues suppress renormalization scale powers: $[{\lambda }_0]=ϵ$ away from four dimensions. The next part restores these powers explicitly with a dimensionless renormalized coupling. Replacing the tree coupling by its renormalized counterpart inside a one-loop residue changes only higher-loop terms.

The free Euclidean action has inverse propagator $D=-m{\partial }_{\tau }^2+\kappa$. Therefore the thermal Green's function obeyswhere ${\delta }_{\beta }$ is the periodic Dirac delta function. It is continuous, has continuous first derivative away from its source, and satisfies $G^{\prime }(0^{+})-G^{\prime }(0^{-})=-1/m$. These conditions follow by integrating the differential equation across the source and using the bosonic thermal trace.

For $\Omega =\sqrt{\kappa /m}$, solving the homogeneous equation on $0<\tau <\beta$ and imposing these conditions givesIn particular $g_0\equiv G^{(0)}(0)=\coth (\beta \Omega /2)/(2m\Omega )$. Equivalently $G^{(0)}={\beta }^{-1}{\sum }_n{e}^{-i{\omega }_n\tau }/(m{\omega }_n^2+\kappa )$ with bosonic Matsubara frequencies ${\omega }_n=2\pi n/\beta$.

For $m^2>0$, the zero-dimensional scalar field theory with Euclidean action $m^2{x}^2/2+\lambda x^4/4!$ has a Gaussian integral normalization $I(0)=\sqrt{2\pi /m^2}$. Its Feynman rules assign $1/m^2$ to each propagator, $-\lambda$ to each quartic vertex, and the reciprocal Feynman-diagram symmetry factor to each vacuum diagram. The integral itself includes disconnected vacuum diagrams; its logarithm selects connected ones. At second order the three possible numbers of edges joining the two vertices are zero, two and four, with Feynman-diagram symmetry factors $128,16,48$ respectively. Their weights add to $35{\lambda }^2/(384m^8)$.

The Wick contractions give the general coefficient $I(0)(-1{)}^v(4v-1)!!/[v!(4!{)}^v{m}^{4v}]$. Successive absolute coefficients have ratio $(4v+3)(4v+1)/[24(v+1)m^4]$, which tends to infinity. Thus this is an asymptotic expansion as $\lambda \downarrow 0$, rather than a convergent power series. For every fixed truncation the Taylor theorem remainder estimate for $e^{-u}$, $u\ge 0$, can be integrated against the Gaussian distribution to justify the corresponding remainder bound.