Skip to content

Time slicing is the finite-partition construction behind the path-integral symbol

∫Dx(t) eiS[x]/ℏ.\int \mathcal D x(t)\,e^{iS[x]/\hbar}.

The continuum notation is compact, but the real definition begins with a finite number of time steps, ordinary integrals over intermediate positions, and short-time kernels with normalization factors. Time slicing is where the measure, action discretization, and operator-ordering choices enter.

Fix endpoints

(xi,ti),(xf,tf),(x_i,t_i), \qquad (x_f,t_f),

and define

T=tf−ti>0.T=t_f-t_i>0.

Choose NN equal time steps

ϵ=TN,\epsilon=\frac{T}{N},

with

tj=ti+jϵ,j=0,1,…,N.t_j=t_i+j\epsilon, \qquad j=0,1,\ldots,N.

The endpoint positions are fixed:

x0=xi,xN=xf.x_0=x_i, \qquad x_N=x_f.

The intermediate positions x1,…,xN−1x_1,\ldots,x_{N-1} are variables to be integrated over.

Time-sliced path with fixed endpoints and intermediate positions

A path integral starts from a finite partition tj=ti+jϵt_j=t_i+j\epsilon, fixed endpoints x0=xix_0=x_i and xN=xfx_N=x_f, and integrations over intermediate positions x1,…,xN−1x_1,\ldots,x_{N-1}.

The propagator kernel is

K(xf,tf;xi,ti)=⟨xf∣U(tf,ti)∣xi⟩.K(x_f,t_f;x_i,t_i) = \langle x_f|U(t_f,t_i)|x_i\rangle.

Use the composition law

U(tf,ti)=U(tN,tN−1)⋯U(t1,t0).U(t_f,t_i) = U(t_N,t_{N-1})\cdots U(t_1,t_0).

Between each short-time factor, insert the position resolution of identity:

I=∫−∞∞dxj ∣xj⟩⟨xj∣.I=\int_{-\infty}^{\infty}dx_j\,|x_j\rangle\langle x_j|.

This gives

K(xf,tf;xi,ti)=∫dx1⋯dxN−1×∏j=0N−1Kϵ(xj+1,xj;tj),K(x_f,t_f;x_i,t_i) = \int dx_1\cdots dx_{N-1} \times \prod_{j=0}^{N-1} K_\epsilon(x_{j+1},x_j;t_j),

where

Kϵ(xj+1,xj;tj)=⟨xj+1∣U(tj+1,tj)∣xj⟩.K_\epsilon(x_{j+1},x_j;t_j) = \langle x_{j+1}|U(t_{j+1},t_j)|x_j\rangle.

The product order follows time order: earlier factors act first and appear to the right in the operator product.

For

H=p22m+V(x),H=\frac{p^2}{2m}+V(x),

one uses a short-time approximation to

U(tj+1,tj)=e−iHϵ/ℏ.U(t_{j+1},t_j) = e^{-iH\epsilon/\hbar}.

A first-order kinetic-potential splitting gives

e−iHϵ/ℏ=e−ip2ϵ/(2mℏ)e−iV(x)ϵ/ℏ+O(ϵ2)e^{-iH\epsilon/\hbar} = e^{-ip^2\epsilon/(2m\hbar)} e^{-iV(x)\epsilon/\hbar} +O(\epsilon^2)

in operator norm for sufficiently controlled bounded models, and formally for many standard unbounded wave-mechanics models. In coordinate representation this leads to

Kϵ(xj+1,xj)≈(m2πiℏϵ)1/2×exp⁡{iℏϵ[m2(xj+1−xjϵ)2−V(xj)]}.K_\epsilon(x_{j+1},x_j) \approx \left( \frac{m}{2\pi i\hbar\epsilon} \right)^{1/2} \times \exp\left\{ \frac{i}{\hbar}\epsilon \left[ \frac{m}{2} \left( \frac{x_{j+1}-x_j}{\epsilon} \right)^2 -V(x_j) \right] \right\}.

The kinetic factor is the free-particle short-time kernel. The potential factor is diagonal in position space.

The short-time splitting is a Trotter product formula in disguise. If

Tp=p22m,T_p=\frac{p^2}{2m},

then

e−i(Tp+V)ϵ/ℏ≈e−iTpϵ/ℏe−iVϵ/ℏ.e^{-i(T_p+V)\epsilon/\hbar} \approx e^{-iT_p\epsilon/\hbar} e^{-iV\epsilon/\hbar}.

A symmetric second-order step is

e−iHϵ/ℏ≈e−iVϵ/(2ℏ)e−iTpϵ/ℏe−iVϵ/(2ℏ).e^{-iH\epsilon/\hbar} \approx e^{-iV\epsilon/(2\hbar)} e^{-iT_p\epsilon/\hbar} e^{-iV\epsilon/(2\hbar)}.

Its coordinate kernel is exactly the free kernel times two endpoint phases:

Ksym(x′,x;ϵ)=K0(x′,x;ϵ)exp⁡ ⁣[−iϵ2ℏ(V(x′)+V(x))].K_{\mathrm{sym}}(x',x;\epsilon) =K_0(x',x;\epsilon) \exp\!\left[-\frac{i\epsilon}{2\hbar}\bigl(V(x')+V(x)\bigr)\right].

The potential contribution is the endpoint average, not V((x′+x)/2)V((x'+x)/2). These expressions agree for linear potentials, but not generally. For mixed position–momentum Hamiltonians, a midpoint phase-space prescription has a separate connection to Weyl ordering; it must not be inferred merely from a symmetric kinetic–potential split.

Using the first-order short-time kernel, the product over slices gives

K(xf,tf;xi,ti)=lim⁡N→∞(m2πiℏϵ)N/2∫dx1⋯dxN−1 eiSN/ℏ,K(x_f,t_f;x_i,t_i) = \lim_{N\to\infty} \left( \frac{m}{2\pi i\hbar\epsilon} \right)^{N/2} \int dx_1\cdots dx_{N-1}\, e^{iS_N/\hbar},

where

SN=∑j=0N−1ϵ[m2(xj+1−xjϵ)2−V(xj)].S_N = \sum_{j=0}^{N-1} \epsilon \left[ \frac{m}{2} \left( \frac{x_{j+1}-x_j}{\epsilon} \right)^2 -V(x_j) \right].

The expression SNS_N is the discretized action. The formal continuum statement

SN→∫titfdt [m2x˙2−V(x)]S_N\to \int_{t_i}^{t_f}dt\, \left[ \frac{m}{2}\dot x^2-V(x) \right]

is a mnemonic for this limiting construction, not a license to forget the discretization.

The path-integral measure symbol hides the finite-dimensional normalization:

∫Dx(t)means schematicallylim⁡N→∞(m2πiℏϵ)N/2∏j=1N−1∫dxj.\int\mathcal D x(t) \quad\text{means schematically}\quad \lim_{N\to\infty} \left( \frac{m}{2\pi i\hbar\epsilon} \right)^{N/2} \prod_{j=1}^{N-1} \int dx_j.

This factor has dimensions. For one coordinate, each short-time kernel has dimension 1/length1/\text{length}, and after N−1N-1 intermediate integrations the full kernel still has dimension 1/length1/\text{length}. Dropping the normalization changes the amplitude, not just an irrelevant convention.

For multidimensional systems, constrained systems, curved configuration spaces, fields, or fermions, the measure changes. Propagators in Multiple Dimensions gives the flat-space dd-coordinate formula and explains why coordinate Jacobians, kinetic operators, and short-time factors must be transformed together.

The continuum notation

K(xf,tf;xi,ti)=∫x(ti)=xix(tf)=xfDx(t) eiS[x]/ℏK(x_f,t_f;x_i,t_i) = \int_{x(t_i)=x_i}^{x(t_f)=x_f} \mathcal D x(t)\, e^{iS[x]/\hbar}

is useful because it displays the physical structure: amplitudes are weighted by the action phase and summed over histories. But the notation suppresses:

  • the time partition;
  • normalization factors;
  • endpoint conditions;
  • the choice of where each potential is evaluated;
  • operator-ordering conventions;
  • convergence or regularization prescriptions.

In real time, the weight eiS/ℏe^{iS/\hbar} is oscillatory. It is not a probability measure on paths.

  • Writing only eiSN/ℏe^{iS_N/\hbar} and dropping the product of short-time prefactors.
  • Integrating over the endpoints even though the propagator fixes them.
  • Forgetting that NN short-time kernels give NN prefactors but only N−1N-1 intermediate position integrals.
  • Treating Dx\mathcal D x as if it were an ordinary translation-invariant measure on a space of smooth functions.
  • Changing Fourier-transform conventions without updating the short-time prefactor.
  • Assuming a midpoint rule, left-endpoint rule, and right-endpoint rule always define the same quantum operator.
  • M. Chaichian and A. Demichev, Path Integrals in Physics, Volume I, Institute of Physics Publishing, 2001.
  • R. P. Feynman and A. R. Hibbs, Quantum Mechanics and Path Integrals, McGraw-Hill, 1965.
  • C. Grosche, “An Introduction into the Feynman Path Integral,” lecture notes (1993), arXiv:hep-th/9302097, operator ordering and path integrals in curved coordinates.
  • H. Kleinert, Path Integrals in Quantum Mechanics, Statistics, Polymer Physics, and Financial Markets, 5th ed., World Scientific, 2009.
  • L. S. Schulman, Techniques and Applications of Path Integration, Wiley, 1981.
  • J. Zinn-Justin, Path Integrals in Quantum Mechanics, Oxford University Press, 2005.
  1. For NN time intervals, why are there N−1N-1 intermediate position integrations?
Solution

The positions x0=xix_0=x_i and xN=xfx_N=x_f are fixed by the propagator endpoints. The remaining positions x1,…,xN−1x_1,\ldots,x_{N-1} are not fixed, so each is integrated over. There are N−1N-1 such intermediate variables.

  1. Check the dimension of the one-dimensional free short-time prefactor.
Solution

The combination ℏϵ/m\hbar\epsilon/m has dimensions of length squared:

ℏϵm∼ML2T−1 TM=L2.\frac{\hbar\epsilon}{m} \sim \frac{ML^2T^{-1}\,T}{M} = L^2.

Therefore

(m2πiℏϵ)1/2\left( \frac{m}{2\pi i\hbar\epsilon} \right)^{1/2}

has dimension 1/L1/L, as a coordinate-space kernel should in one dimension.

  1. For V(x)=λx4V(x)=\lambda x^4, compare the potential in a symmetric split from x=−ax=-a to x′=ax'=a with its value at the spatial midpoint.
Solution

The endpoint average is [V(−a)+V(a)]/2=λa4[V(-a)+V(a)]/2=\lambda a^4, whereas the midpoint value is V(0)=0V(0)=0. Thus the two finite-step phases differ by −ϵλa4/ℏ-\epsilon\lambda a^4/\hbar. Symmetry of an operator split does not make these prescriptions identical.

  1. In the first-order discretization above, why does V(xj)V(x_j) appear rather than an average over the whole slice?
Solution

The displayed short-time kernel used a left-endpoint splitting,

e−iHϵ/ℏ≈e−iTpϵ/ℏe−iV(x)ϵ/ℏ.e^{-iH\epsilon/\hbar} \approx e^{-iT_p\epsilon/\hbar} e^{-iV(x)\epsilon/\hbar}.

Since V(x)V(x) acts on the ket ∣xj⟩|x_j\rangle, the potential factor contributes e−iV(xj)ϵ/ℏe^{-iV(x_j)\epsilon/\hbar}. A symmetric Trotter splitting gives the average of the potential values at the two endpoints. This is distinct from evaluating the potential at the midpoint of those positions.