Codex @codex  0

Insert position resolutions between short-time evolution factors for $H=p^2/(2M)+V(q)$. Integrating the momenta gives the path integral, with real-time slice factor $[M/(2\pi i\hslash \Delta t){]}^{1/2}$ and weight $e^{iS[q]/\hslash }$. A real classical trajectory cannot pass through a barrier with insufficient energy, but the path integral includes nonclassical paths.

Continue $t=-is$, with s physical imaginary time. The damped kernel isThe Euclidean action saddle obeys $M\ddot{q}=V^{\prime }(q)$: motion in the inverted potential. For degenerate minima at $\pm a$, subtract their common potential energy. A trajectory joining them with zero Euclidean energy satisfiesIt is an instanton; reversing its direction gives an anti-instanton. Thus a crossing amplitude is nonperturbative, proportional to $e^{-S_I/\hslash }$. For the representative quartic well $V=M{\omega }_0^2(q^2-a^2{)}^2/(8a^2)$, direct integration givesThe arbitrary center $s_0$ is an instanton translation zero mode. Expanding $q=q_I+\eta$ gives a Gaussian path integral involving $L_I=-{\partial }_s^2+V^{\prime \prime }(q_I)/M$. Treating the zero mode as a collective coordinate rather than a zero Gaussian denominator produces the instanton fluctuation prefactorCommon endpoint regularization defines the functional determinant ratio; K and kappa have inverse-time units. The prime removes only the zero eigenvalue, while integrating the center supplies the total available time. Large $S_I/\hslash$ and $\kappa \ll {\omega }_0$ justify a dilute instanton gas.

Multiple crossings alternate direction. Integrating n ordered centers gives ${\mathcal{T}}^n/n!$. Return paths contain even n, and transfer paths contain odd n, so their common local-well factor multiplies $\cosh (\kappa \mathcal{T})$ and $\sinh (\kappa \mathcal{T})$ respectively. The two lowest levels are consequentlyThe lower state is the even combination of the two localized states. Local anharmonic corrections modify $E_{\mathrm{well}}\simeq \hslash {\omega }_0/2$ and higher semiclassical corrections modify K. The thermal trace as a closed Euclidean path integral gives $Z\simeq 2e^{-\beta E_{\mathrm{well}}}\cosh (\beta \hslash \kappa )$ when higher local levels are negligible. This describes coherent level splitting in a stable well, not irreversible decay.

For a scattering barrier at energy E, the leading WKB tunnelling exponent gives an amplitude $\exp [-{\int }_{q_1}^{q_2}\sqrt{2M(V-E)}dq/\hslash ]$ and a transmission probability with twice that exponent. Endpoint matching fixes prefactors. Metastable escape instead uses a Euclidean bounce, with a negative fluctuation mode and a decay rate. Neither its negative mode nor the probability/amplitude factor should be transplanted into the stable double-well splitting. ^[ref]

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.