Path integral formulation

This article is about a formulation of quantum mechanics. For integrals along a path, also known as line or contour integrals, see line integral.

Quantum field theory

Feynman diagram

History of... Background Quantum mechanics Poincaré symmetry Gauge theory • Field theory Spontaneous symmetry breaking Tools Quantization • Renormalization Partition function • Vacuum state Anomaly • Expectation value Wightman axioms • Feynman diagram Effective field theory LSZ reduction formula Wick's theorem Equations Klein-Gordon equation Dirac equation Wheeler-deWitt equation Standard model Quantum electrodynamics Electroweak interaction Quantum chromodynamics Higgs mechanism Faddeev-Popov ghosts Yang-Mills Incomplete theories Quantum gravity Supersymmetry String theory Theory of everything Scientists Feynmann • Schwinger • Tomonaga • Dyson Gross • Politzer •Wilczek Yang • Mills • 't Hooft Klein • Weinberg • Salam Kibble • Higgs Bekenstein· Hawking · others

This box: view • talk • edit

The path integral formulation of quantum mechanics is a description of quantum theory which generalizes the action principle of classical mechanics. It replaces the classical notion of a single, unique history for a system with a sum, or functional integral, over an infinity of possible histories to compute a quantum amplitude.

The path integral formulation was developed in 1948 by Richard Feynman. Some preliminaries were worked out earlier, in the course of his doctoral thesis work with John Archibald Wheeler.

This formulation has proved crucial to the subsequent development of theoretical physics, since it provided the basis for the grand synthesis of the 1970s called the renormalization group which unified quantum field theory with statistical mechanics. If we realize that the Schrödinger equation is essentially a diffusion equation with an imaginary diffusion constant, then the path integral is a method for the enumeration of random walks. For this reason path integrals had also been used in the study of Brownian motion and diffusion before they were introduced in quantum mechanics.

These are just three of the paths that contribute to the quantum amplitude for a particle moving from point A at some time t₀ to point B at some other time t₁.

Formulating quantum mechanics

The path integral method is an alternative formulation of quantum mechanics. The canonical approach, pioneered by Erwin Schrödinger, Werner Heisenberg and Paul Dirac paid great attention to wave-particle duality and the resulting uncertainty principle by replacing Poisson brackets of classical mechanics by commutators between operators in quantum mechanics. The Hilbert space of quantum states and the superposition law of quantum amplitudes follows. The path integral starts from the superposition law, and exploits wave-particle duality to build a generating function for quantum amplitudes.

Abstract formulation

Feynman proposed the following postulates:

1. The probability for any fundamental event is given by the square modulus of a complex amplitude.
2. The amplitude for some event is given by adding together the contributions of all the histories which include that event.
3. The amplitude a certain history contributes is proportional to Failed to parse (Missing texvc executable; please see math/README to configure.): e^{i S/\hbar}

, where Failed to parse (Missing texvc executable; please see math/README to configure.): \hbar

is reduced Planck's constant and S is the action of that history, given by the time integral of the Lagrangian along the corresponding path in the phase space of the system.

In order to find the overall probability amplitude for a given process, then, one adds up, or integrates, the amplitude of postulate 3 over the space of all possible histories of the system in between the initial and final states, including histories that are absurd by classical standards. In calculating the amplitude for a single particle to go from one place to another in a given time, it would be correct to include histories in which the particle describes elaborate curlicues, histories in which the particle shoots off into outer space and flies back again, and so forth. The path integral assigns all of these histories amplitudes of equal magnitude but with varying phase, or argument of the complex number. The contributions that are wildly different from the classical history are suppressed only by the interference of similar, canceling histories (see below).

The mathematical technique of path integrals does not imply that real particles must actually follow the paths so constructed. Mathematical expansions of functions by other functions are a general technique, and as such the functions used are not required to have any physical interpretation at all. They are usually constructed for mathematical convenience, with no necessary analogy to the physical model that they are modeling. Indeed, in the case of the Feynman path integral, the integration is over imaginary time, so the relevance of the paths to the particle's real physical path is open to debate.

Feynman showed that his formulation of quantum mechanics is equivalent to the canonical approach to quantum mechanics. An amplitude computed according to Feynman's principles will also obey the Schrödinger equation for the Hamiltonian corresponding to the given action.

Recovering the action principle

Feynman was extending a classic paper by Paul Dirac on the quantum equivalent of the classical action principle. Dirac knew from earlier semiclassical work that the Maupertuis action, the action-angle variable action, is the classical analog of the phase of the wavefunction in an energy eigenstate. He noted that the transformation between the Maupertuis action and the usual Lagrange action replaces a description of a motion at fixed energy with description of a motion evolving between two fixed times. The phase interpretation of the Maupertuis action implied that the Lagrange action was also a phase, and that it was associated to an infinitesimal path starting and ending at two nearby positions at two nearby times. He noted that the full quantum transition amplitude between two points at two times can be heuristically thought of as a sum of this phase factor over all possible paths, where the phase over each infinitesimal segment is given by the Lagrange action.

Dirac's work was only heuristic because he did not provide a precise prescription to calculate the sum, and he did not show that one could recover the Schrödinger equation or the canonical commutation relations from this rule. This was done by Feynman. Both noted that in the limit of action that is large compared to Planck's constant Failed to parse (Missing texvc executable; please see math/README to configure.): \hbar , the path integral is dominated by solutions which are in the neighbourhood of stationary points of the action, since there the amplitudes of similar histories will tend to constructively interfere with one another. Conversely, for paths that are far from the stationary points of the action, the complex phase of the amplitude calculated according to postulate 3 will vary rapidly for similar paths, and amplitudes will tend to cancel. Therefore the important parts of the integral—the significant possibilities—in the limit of large action simply consist of solutions of the Euler-Lagrange equation, and classical mechanics is correctly recovered.

Classical action principles are puzzling because of their seemingly teleological quality: instead of predicting the future from initial conditions, they give a combination of initial conditions and final conditions to find the path in between, as if the system somehow knows where it's going to end up. The path integral explains why this works. The system doesn't have to know in advance where it's going; the path integral simply calculates the probability amplitude for any given process, and the path goes everywhere. After a long enough time, interference effects guarantee that only the contributions from the stationary points of the action give histories with appreciable probabilities.

Concrete formulation

Feynman's postulates are somewhat ambiguous in that they do not define what an "event" is or the exact proportionality constant in postulate 3. The proportionality problem can be solved by simply normalizing the path integral by dividing the amplitude by the square root of the total probability for something to happen (resulting in that the total probability given by all the normalized amplitudes will be 1, as we would expect). Generally speaking one can simply define the "events" in an operational sense for any given experiment.

The equal magnitude of all amplitudes in the path integral tends to make it difficult to define it such that it converges and is mathematically tractable. For purposes of actual evaluation of quantities using path-integral methods, it is common to give the action an imaginary part in order to damp the wilder contributions to the integral, then take the limit of a real action at the end of the calculation. In quantum field theory this takes the form of Wick rotation.

There is some difficulty in defining a measure over the space of paths. In particular, the measure is concentrated on "fractal-like" distributional paths.

Time-slicing definition

For a particle in a smooth potential, the path integral is approximated by Feynman as the small-step limit over zig-zag paths, which in one dimension is a product of ordinary integrals. For the motion of the particle from position Failed to parse (Missing texvc executable; please see math/README to configure.): x_a

at time Failed to parse (Missing texvc executable; please see math/README to configure.): t_a
to Failed to parse (Missing texvc executable; please see math/README to configure.): x_b
at time Failed to parse (Missing texvc executable; please see math/README to configure.): t_b

, the time sequence Failed to parse (Missing texvc executable; please see math/README to configure.): t_a<t_1<...<t_{n-1}<t_n<t_b

 can be divided up into n little segments Failed to parse (Missing texvc executable; please see math/README to configure.): t_j-t_{j-1},\,\, j=1,..., n\,,
  of fixed duration Failed to parse (Missing texvc executable; please see math/README to configure.): \Delta t\,\,(=\epsilon )=\frac{t_b-t_a}{n+1}

. (The one remaining segment can be neglected, since finally the limit Failed to parse (Missing texvc executable; please see math/README to configure.): n\to\infty

is considered.)

This process is called time slicing.
An approximation for the path integral can be computed as proportional to

Failed to parse (Missing texvc executable; please see math/README to configure.): \int\limits_{x_1=x_a}^{x_1=x_b} \ldots \int\limits_{x_n=x_a}^{x_n=x_b} \ \exp \left(\frac{{\rm i}}{\hbar}\int\limits_{t_a}^{t_b} \mathcal L(x(t),v(t), t)\,\mathrm{d}t\right) \, \mathrm{d}x_1 \cdots \mathrm{d}x_{n}

where Failed to parse (Missing texvc executable; please see math/README to configure.): \mathcal L(x,v,t)

is the Lagrangian of the 1d-system with position variable x(t) and velocity Failed to parse (Missing texvc executable; please see math/README to configure.):  v=\dot x(t)
considered (see below), and Failed to parse (Missing texvc executable; please see math/README to configure.): {\rm d}x_j
corresponds to the position at the j-th time step, if the time integral is approximated by a sum of n terms.

In the limit of Failed to parse (Missing texvc executable; please see math/README to configure.): n

going to infinity, this becomes a functional integral, which - apart from a nonessential factor - is directly the product of the probability amplitudes Failed to parse (Missing texvc executable; please see math/README to configure.): \langle x_a,t_a|x_b, t_b\rangle
- more precisely, since one must work with a continuous spectrum, the respective densities - to find the quantum mechanical particle at Failed to parse (Missing texvc executable; please see math/README to configure.):  t_a
in the initial state Failed to parse (Missing texvc executable; please see math/README to configure.): x_a
and at Failed to parse (Missing texvc executable; please see math/README to configure.): t_b
in the final state Failed to parse (Missing texvc executable; please see math/README to configure.): x_b

.

Actually Failed to parse (Missing texvc executable; please see math/README to configure.): \mathcal L

is the classical Lagrangian of the one-dimensional system considered,  Failed to parse (Missing texvc executable; please see math/README to configure.):  \mathcal L(x,\dot x , t)=p\cdot \dot x -\mathcal H(x,p,t)

, where Failed to parse (Missing texvc executable; please see math/README to configure.): \mathcal H

is the Hamiltonian, with Failed to parse (Missing texvc executable; please see math/README to configure.):  p=\frac {\partial \mathcal L}{\partial \dot x} 

, and the above-mentioned "zigzagging" corresponds to the appearance of the terms Failed to parse (Missing texvc executable; please see math/README to configure.): \exp\left (\frac{{\rm i}}{\hbar}\epsilon\, \,\sum_{j=1}^{n}\mathcal L(\tilde x_{j},\frac{x_j-x_{j-1}}{\epsilon},j)\right )

in the Riemannian sum approximating the time integral, which are finally integrated over Failed to parse (Missing texvc executable; please see math/README to configure.): x_1
to Failed to parse (Missing texvc executable; please see math/README to configure.): x_n
with the integration measure Failed to parse (Missing texvc executable; please see math/README to configure.): {\rm d}x_1\cdot ...\cdot {\rm d}x_n\,.\,\,\tilde x_j
is an arbitrary value of the interval corresponding to j, e.g. its center, Failed to parse (Missing texvc executable; please see math/README to configure.): \frac{x_j+x_{j-1}}{2}\,.
 

Thus in contrast to classical mechanics not only the stationary path contributes but actually all virtual paths between the initial and the final point.

The limit Failed to parse (Missing texvc executable; please see math/README to configure.): n\to\infty

does not, however, exist for the most important quantum-mechanical  systems, the atoms, due to the 

singularity of the Coulomb potential Failed to parse (Missing texvc executable; please see math/README to configure.): e^2/r \,

at the origin. The problem was solved in 1979 

by H. Duru and Hagen Kleinert ^[1][2] by choosing Failed to parse (Missing texvc executable; please see math/README to configure.): \Delta t\,\,(=\epsilon )

proportional to Failed to parse (Missing texvc executable; please see math/README to configure.): r
 and going to new coordinates whose square length is equal to Failed to parse (Missing texvc executable; please see math/README to configure.): r

(Duru-Kleinert transformation).

Free Particle

The path integral representation gives the quantum amplitude to go from point x to point y as an integral over all paths. For a free particle action (Failed to parse (Missing texvc executable; please see math/README to configure.): \scriptstyle m=1,\hbar=1

Failed to parse (Missing texvc executable; please see math/README to configure.): S= \int {\dot{x}^2\over 2} dt \,

the integral can be evaluated explicitly.

To do this, it is conceptually convenient to start without the factor i in the exponential, so that large deviations are suppressed by small numbers, not by cancelling oscillatory contributions.

Failed to parse (Missing texvc executable; please see math/README to configure.): K(x-y;T) = \int_{x(0)=x}^{x(T)=y} e^{-\int_0^T {\dot{x}^2\over 2} dt} Dx \,

Splitting the integral into time slices:

Failed to parse (Missing texvc executable; please see math/README to configure.): K(x,y;T) = \int_{x(0)=x}^{x(T)=y} \Pi_t e^{-{1\over 2} ({x(t+\epsilon) -x(t) \over \epsilon})^2 \epsilon } Dx \,

where the Dx is interpreted as a finite collection of integrations at each integer multiple of Failed to parse (Missing texvc executable; please see math/README to configure.): \epsilon . Each factor in the product is a Gaussian as a function of Failed to parse (Missing texvc executable; please see math/README to configure.): x(t+\epsilon)

centered at x(t) with variance Failed to parse (Missing texvc executable; please see math/README to configure.): \epsilon

. The multiple integrals are a repeated convolution of this Gaussian Failed to parse (Missing texvc executable; please see math/README to configure.): G_\epsilon

with copies of itself at adjacent times.

Failed to parse (Missing texvc executable; please see math/README to configure.): K(x-y;T) = G_\epsilon*G_\epsilon ... *G_\epsilon \,

Where the number of convolutions is Failed to parse (Missing texvc executable; please see math/README to configure.): T/\epsilon . The result is easy to evaluate by taking the fourier transform of both sides, so that the convolutions becomes multiplications.

Failed to parse (Missing texvc executable; please see math/README to configure.): K(p;T) = G_\epsilon(p)^{T \over \epsilon} \,

The Fourier transform of the Gaussian G is another Gaussian of reciprocal variance:

Failed to parse (Missing texvc executable; please see math/README to configure.): G_\epsilon(p) = e^{-\epsilon {p^2\over 2} } \,

and the result is:

Failed to parse (Missing texvc executable; please see math/README to configure.): K(p;T) = e^{-T {p^2 \over 2}} \,

The Fourier transform gives K, and it is a Gaussian again with reciprocal variance:

Failed to parse (Missing texvc executable; please see math/README to configure.): K(x-y;T) \propto e^{ (x-y)^2\over 2T} \,

The proportionality constant is not really determined by the time slicing approach because it depends on an arbitrary choice of normalization constant at each time. If the slicing is nonuniform, the constant can easily be shifted this way and that. Only the ratio of values for different endpoint choices is determined. This only really matters when comparing functions at different times, but it is important when finding differential equations for the state.

The result has a probability interpretation. The sum over all paths of the exponential factor can be seen as the sum over each path of the probability of selecting that path. The probability is the product over each segment of the probability of selecting that segment, so that each segment is probabilistically independently chosen. The fact that the answer is a Gaussian spreading linearly in time is the central limit theorem, which can be interpreted as the first historical evaluation of a statistical path integral.

The probability interpretation gives a natural normalization choice. The path integral should be defined so that:

Failed to parse (Missing texvc executable; please see math/README to configure.): \int K(x-y;T) dy = 1 \,

This condition normalizes the Gaussian, and produces a Kernel which obeys the diffusion equation:

Failed to parse (Missing texvc executable; please see math/README to configure.): {d\over dt} K(x;T) = {\nabla^2 \over 2} K \,

For oscillatory path integrals, ones with an i in the numerator, the time-slicing produces convolved Gaussians, just as before. Now, however, the convolution product is marginally singular since it requires careful limits to evaluate the oscillating integrals. To make the factors well defined, the easiest way is to add a small imaginary part to the time increment Failed to parse (Missing texvc executable; please see math/README to configure.): \epsilon . Then the same convolution argument as before gives the propagation kernel:

Failed to parse (Missing texvc executable; please see math/README to configure.): K(x-y;T) \propto e^{i(x-y)^2 \over 2T} \,

Which, with the same normalization as before (not the sum-squares normalization! this function has a divergent norm), obeys a free Schrödinger equation

Failed to parse (Missing texvc executable; please see math/README to configure.): {d\over dt} K(x;T) = {\rm i} {\nabla^2 \over 2} K \,

This means that any superposition of K's will also obey the same equation, by linearity. Defining

Failed to parse (Missing texvc executable; please see math/README to configure.): \psi_t(y) = \int \psi_0(x) K(x-y;t) dx = \int \psi_0(x) \int_{x(0)=x}^{x(t)=y} e^{iS} Dx \,

then Failed to parse (Missing texvc executable; please see math/README to configure.): \psi_t

obeys the free Schrödinger equation just as K does:

Failed to parse (Missing texvc executable; please see math/README to configure.): {\rm i}{\partial \over \partial t} \psi_t = - {\nabla^2\over 2} \psi_t \,

Schrödinger Equation

The path integral reproduces the Schrödinger equation for the initial and final state even when a potential is present. This is easiest to see by taking a path-integral over infinitsimally separated times.

Failed to parse (Missing texvc executable; please see math/README to configure.): \psi(y;t+\epsilon) = \int_{-\infty}^\infty dx\;\;\psi(x;t)\int_{x(t)=x}^{x(t+\epsilon)=y} e^{{\rm i}\int_t^{t+\epsilon} (\dot{x}^2 - V(x)) dt} Dx(t) \qquad (1)

Since the time separation is infinitesimal and the cancelling oscillations become severe for large values of Failed to parse (Missing texvc executable; please see math/README to configure.): \scriptstyle \dot{x} , the path integral has most weight for y close to x. In this case, to lowest order the potential energy is constant, and only the kinetic energy contribution is nontrivial. The exponential of the action is

Failed to parse (Missing texvc executable; please see math/README to configure.): e^{{\rm i}\epsilon V(x)} e^{{\rm i}\dot{x}^2\epsilon} \,

The first term rotates the phase of Failed to parse (Missing texvc executable; please see math/README to configure.): \psi(x)

locally by an amount proportional to the potential energy. The second term is the free particle propagator, corresponding to i times a diffusion process. To lowest order in Failed to parse (Missing texvc executable; please see math/README to configure.): \epsilon
they are additive; in any case  one has with (1):

Failed to parse (Missing texvc executable; please see math/README to configure.): \psi(y;t+\epsilon) \approx \int \psi(x;t) e^{{\rm i}\epsilon V(x)} e^{{\rm i}(x-y)^2 \over 2\epsilon} dx \,.

As mentioned, the spread in Failed to parse (Missing texvc executable; please see math/README to configure.): \psi

is diffusive from the free particle propagation, with an extra infinitesimal rotation in phase which slowly varies from point to point from the potential:

Failed to parse (Missing texvc executable; please see math/README to configure.): \frac{\partial\psi}{\partial t} = {\rm i}\cdot \{ \frac{1}{2}\nabla^2 - V(x)\}\psi \,

and this is the Schrödinger equation. Note that the normalization of the path integral needs to be fixed in exactly the same way as in the free particle case. An arbitrary continuous potential does not affect the normalization, although singular potentials require careful treatment.

Equations of Motion

Since the states obey the Schrödinger equation, the path integral must reproduce the Heisenberg equations of motion for the averages of x and Failed to parse (Missing texvc executable; please see math/README to configure.): \dot{x}

variables, but it is instructive to see this directly. The direct approach shows that the expectation values calculated from the path integral reproduce the usual ones of quantum mechanics.

Start by considering the path integral with some fixed initial state

Failed to parse (Missing texvc executable; please see math/README to configure.): \int \psi_0(x) \int_{x(0)=x} e^{{\rm i}S(x,\dot{x})} Dx \,

Now note that Failed to parse (Missing texvc executable; please see math/README to configure.): x(t)

at each separate time is a separate integration variable. So it is legitimate to change variables in the integral by shifting: Failed to parse (Missing texvc executable; please see math/README to configure.):  x(t)=u(t)+\epsilon(t)
where Failed to parse (Missing texvc executable; please see math/README to configure.): \epsilon(t)
is a different shift at each time but Failed to parse (Missing texvc executable; please see math/README to configure.): \epsilon(0)=\epsilon(T)=0

, since the endpoints are not integrated:

Failed to parse (Missing texvc executable; please see math/README to configure.): \int \psi_0(x) \int_{u(0)=x} e^{{\rm i}S(u+\epsilon,\dot{u}+\dot{\epsilon})} Du \,

The change in the integral from the shift is, to first infinitesimal order in epsilon:

Failed to parse (Missing texvc executable; please see math/README to configure.): \int \psi_0(x) \int_{u(0)=x} \left( \int {\partial S \over \partial u } \epsilon + { \partial S \over \partial \dot{u} } \dot{\epsilon} dt \right) e^{iS} Du \,

which, integrating by parts in t, gives:

Failed to parse (Missing texvc executable; please see math/README to configure.): - \int \psi_0(x) \int_{u(0)=x} \left( \int ({d\over dt} {\partial S\over \partial \dot{u}} - {\partial S \over \partial u})\epsilon(t) dt \right) e^{iS} Du \,

But this was just a shift of integration variables, which doesn't change the value of the integral for any choice of Failed to parse (Missing texvc executable; please see math/README to configure.): \epsilon(t) . The conclusion is that this first order variation is zero for an arbitrary initial state and at any arbitrary point in time:

Failed to parse (Missing texvc executable; please see math/README to configure.): \langle \psi_0| {\delta S \over \delta x}(t) |\psi_0 \rangle = 0 \,

this is the Heisenberg equations of motion.

If the action contains terms which multiply Failed to parse (Missing texvc executable; please see math/README to configure.): \scriptstyle \dot{x}

and Failed to parse (Missing texvc executable; please see math/README to configure.): x

, at the same moment in time, the manipulations above are only heuristic, because the multiplication rules for these quantities is just as noncommuting in the path integral as it is in the operator formalism.

Canonical Commutation Relations

The formulation of the path integral does not make it clear at first sight that the quantities x and p do not commute. In the path integral, these are just integration variables and they have no obvious ordering. Feynman discovered that the non-commutativity is still there ^[3] .

To see this, consider the simplest path integral, the brownian walk. This is not yet quantum mechanics, so in the path-integral the action is not multiplied by i:

Failed to parse (Missing texvc executable; please see math/README to configure.): S= \int \left( {dx \over dt} \right)^2 dt \,

The quantity x(t) is fluctuating, and the derivative is defined as the limit of a discrete difference.

Failed to parse (Missing texvc executable; please see math/README to configure.): {dx \over dt} = {x(t+\epsilon) - x(t) \over \epsilon} \,

Note that the distance that a random walk moves is proportional to Failed to parse (Missing texvc executable; please see math/README to configure.): \sqrt{t} , so that:

Failed to parse (Missing texvc executable; please see math/README to configure.): x(t+\epsilon) - x(t) \approx \sqrt{\epsilon} \,

This shows that the random walk is not differentiable, since the ratio that defines the derivative diverges with probability one.

The quantity Failed to parse (Missing texvc executable; please see math/README to configure.): \scriptstyle x\dot{x}

is ambiguous, with two possible meanings:

Failed to parse (Missing texvc executable; please see math/README to configure.): [1] = x { dx\over dt} = x(t) {(x(t+\epsilon) - x(t)) \over \epsilon } \,

Failed to parse (Missing texvc executable; please see math/README to configure.): [2] = x {dx \over dt} = x(t+\epsilon) {(x(t+\epsilon) - x(t)) \over \epsilon} \,

In ordinary calculus, the two are only different by an amount which goes to zero as Failed to parse (Missing texvc executable; please see math/README to configure.): \epsilon

goes to zero. But in this case, the difference between the two is not zero:

Failed to parse (Missing texvc executable; please see math/README to configure.): [2] - [1] = {( x(t + \epsilon) - x(t) )^2 \over \epsilon} \approx {\epsilon \over \epsilon} \,

give a name to the value of the difference for any one random walk:

Failed to parse (Missing texvc executable; please see math/README to configure.): {(x(t+\epsilon)- x(t))^2 \over \epsilon} = f(t) \,

and note that Failed to parse (Missing texvc executable; please see math/README to configure.): f(t)

is a rapidly fluctuating statistical quantity, whose average value is 1,  i.e. a normalized "Gaussian process". The fluctuations of such a quantity can be described by a statistical Lagrangian Failed to parse (Missing texvc executable; please see math/README to configure.): \mathcal L =  (f(t)-1)^2 
, and the equations of motion for f derived from extremizing the action S corresponding to Failed to parse (Missing texvc executable; please see math/README to configure.):  \mathcal L
just set it equal to 1. In physics, such a quantity is "equal to 1 as an operator identity". In mathematics, it "weakly converges to 1". In either case, it is 1 in any expectation value, or when averaged over any interval, or for all practical purpose.

Defining the time order to be the operator order:

Failed to parse (Missing texvc executable; please see math/README to configure.): [x, \dot x] = x {dx\over dt} - {dx \over dt} x = 1 \,

This is called the Ito lemma in stochastic calculus, and the (euclideanized) canonical commutation relations in physics.

For a general statistical action, a similar argument shows that

Failed to parse (Missing texvc executable; please see math/README to configure.): [x , {\partial S \over \partial \dot x} ] = 1 \,

And in quantum mechanics, the extra imaginary unit in the action converts this to the canonical commutation relation.

Failed to parse (Missing texvc executable; please see math/README to configure.): [x,p ] ={\rm i} \,

Particle in curved space

For a particle in curved space the kinetic term depends on the position and the above time slicing cannot be applied, this being a manifestation of the notorious operator ordering problem in Schrödinger quantum mechanics. One may, however, solve this problem by transforming the time-sliced flat-space path integral to curved space using a multivalued coordinate transformation (nonholonomic mapping explained here).

The path integral and the partition function

The path integral is just the generalization of the integral above to all quantum mechanical problems—

Failed to parse (Missing texvc executable; please see math/README to configure.): Z = \int Dx\, e^{{\rm i}\mathcal{S}[x]/\hbar}

 where  Failed to parse (Missing texvc executable; please see math/README to configure.): \mathcal{S}[x]=\int_0^T \mathrm{d}t L[x(t)]

is the action of the classical problem in which one investigates the path starting at time t=0 and ending at time t = T, and Dx denotes integration over all paths. In the classical limit, Failed to parse (Missing texvc executable; please see math/README to configure.): \mathcal{S}[x] \gg \hbar , the path of minimum action dominates the integral, because the phase of any path away from this fluctuates rapidly and different contributions cancel.

The connection with statistical mechanics follows. Considering only paths which begin and end in the same configuration, perform the Wick rotation Failed to parse (Missing texvc executable; please see math/README to configure.): t\to{\rm i}t , i.e., make time imaginary, and integrate over all possible beginning/ending configurations. The path integral now resembles the partition function of statistical mechanics defined in a canonical ensemble with temperature Failed to parse (Missing texvc executable; please see math/README to configure.): 1/T\hbar . Strictly speaking, though, this is the partition function for a statistical field theory.

Clearly, such a deep analogy between quantum mechanics and statistical mechanics cannot be dependent on the formulation. In the canonical formulation, one sees that the unitary evolution operator of a state is given by

Failed to parse (Missing texvc executable; please see math/README to configure.): |\alpha;t\rangle=e^{-{\rm i}Ht / \hbar}|\alpha;0\rangle

where the state α is evolved from time t=0. If one makes a Wick rotation here, and finds the amplitude to go from any state, back to the same state in (imaginary) time iT is given by

Failed to parse (Missing texvc executable; please see math/README to configure.): Z={\rm Tr} [e^{-HT / \hbar}]

which is precisely the partition function of statistical mechanics for the same system at temperature quoted earlier. One aspect of this equivalence was also known to Schrödinger who remarked that the equation named after him looked like the diffusion equation after Wick rotation.

Quantum field theory

The path integral formulation was very important for the development of quantum field theory. Both the Schrödinger and Heisenberg approaches to quantum mechanics single out time, and are not in the spirit of relativity. For example, the Heisenberg approach requires that scalar field operators obey the commutation relation

Failed to parse (Missing texvc executable; please see math/README to configure.): [\phi(x),\partial_t \phi(y) ] = {\rm i} \delta(x-y)

for x and y two simulataneous spatial positions, and this is not a relativistically invariant concept. The results of a calculation are covariant at the end of the day, but the symmetry is not apparent in intermediate stages. If naive field theory calculations did not produce infinite answers in the continuum limit, this would not have been such a big problem--- it would just have been a bad choice of coordinates. But the lack of symmetry means that the infinite quantities must be cut off, and the bad coordinates makes it nearly impossible to cut off the theory without spoiling the symmetry. This makes it difficult to extract the physical predictions, which require a careful limiting procedure.

The problem of lost symmetry also appears in classical mechanics, where the Hamiltonian formulation also superficially singles out time. The Lagrangian formulation makes the relativistic invariance apparent. In the same way, the path integral is manifestly relativistic. It reproduces the Schrödinger equation as the evolution equation for the state at each time slice, and it also includes the Heisenberg equations of motion as a differential identity obeyed by averages of the integration variables at nearby times. It includes the canonical commutation relations in a natural way, and it extends them using the relativistic symmetry to operator product rules which are new relations difficult to extract from the old formalism.

Further, different choices of canonical variables lead to very different seeming formulations of the same theory. The transformations between the variables can be very complicated, but the path integral makes them into reasonably straightforward changes of integration variables. For these reasons, the Feynman path integral has made earlier formalisms largely obsolete.

The price of a path integral representation is that the unitarity of a theory is no longer self evident, but it can be proven by changing variables to some canonical representation. The path integral itself also deals with larger mathematical spaces than is usual, which requires more careful mathematics not all of which has been fully worked out. The path integral historically was not immediately accepted, partly because it took many years to incorporate fermions properly. This required physicists to invent an entirely new mathematical object--- the grassman variable--- which also allowed changes of variables to be done naturally, as well as allowing constrained quantization.

The integration variables in the path integral are subtly non-commuting. The value of the product of two field operators at what looks like the same point depends on how the two points are ordered in space and time. This makes some naive identities fail.

The propagator

In relativistic theories, there is both a particle and field representation for every theory. The field representation is a a sum over all field configurations, and the particle representation is a sum over different particle paths.

The nonrelativistic formulation is traditionally given in terms of particle paths, not fields. There, the path integral in the usual variables, with fixed boundary conditions, gives the probability amplitude for a particle to go from point x to point y in time T.

Failed to parse (Missing texvc executable; please see math/README to configure.): K(x,y;T) = \langle x;T|y;0 \rangle = \int_{x(0)=x}^{x(T)=y} e^{i S[x]} Dx \,

This is called the propagator. Superposing different values of the initial position Failed to parse (Missing texvc executable; please see math/README to configure.): x

with an arbitrary initial state Failed to parse (Missing texvc executable; please see math/README to configure.): \psi_0(x)
constructs the final state.

Failed to parse (Missing texvc executable; please see math/README to configure.): \psi_T(y) = \int_{x} \psi_0(x) K(x,y;T) dx = \int^{x(T)=y} \psi_0(x(0)) e^{i S[x]} Dx \,

For a spatially homogenous system, where K(x, y) is a only a function of (x-y), the integral is a convolution, the final state is the initial state convolved with the propagator.

Failed to parse (Missing texvc executable; please see math/README to configure.): \psi_T = \psi_0 * K(;T) \,

For a free particle of mass m, the propagator can be evaluated either explicitly from the path integral or by noting that the Schrödinger equation is a diffusion equation in imaginary time and the solution must be a normalized Gaussian:

Failed to parse (Missing texvc executable; please see math/README to configure.): K(x,y;T) \propto e^{i m(x-y)^2\over 2T}

Taking the Fourier transform in (x-y) produces another Gaussian:

Failed to parse (Missing texvc executable; please see math/README to configure.): K(p;T) = e^{i T p^2\over 2m}

and in p-space the proportionality factor here is constant in time, as will be verified in a moment. The Fourier transform in time, extending K(p;T) to be zero for negative times, gives the Green's Function, or the frequency space propagator:

Failed to parse (Missing texvc executable; please see math/README to configure.): G_F(p,E) = {-i \over E - {\vec{p}^2\over 2m} + i\epsilon}

Which is the reciprocal of the operator which annihilates the wavefunction in the Schrödinger equation, which wouldn't have come out right if the proportionality factor weren't constant in the p-space representation.

The infinitesimal term in the denominator is a small positive number which guarantees that the inverse Fourier transform in E will be nonzero only for future times. For past times, the inverse Fourier transform contour closes toward values of E where there is no singularity. This guarantees that K propagates the particle into the future and is the reason for the subscript on G. The infinitesimal term can be interpreted as an infinitesimal rotation toward imaginary time.

It is also possible to reexpress the nonrelativistic time evolution in terms of propagators which go toward the past, since the Schrödinger equation is time-reversible. The past propagator is the same as the future propagator except for the obvious difference that it vanishes in the future, and in the gaussian Failed to parse (Missing texvc executable; please see math/README to configure.): t

is replaced by Failed to parse (Missing texvc executable; please see math/README to configure.): -t

. In this case, the interpretation is that these are the quantities to convolve the final wavefunction so as to get the initial wavefunction.

Failed to parse (Missing texvc executable; please see math/README to configure.): G_B(p,E) = { - i \over - E - {i\vec{p}^2\over 2m} + i\epsilon}

Given the nearly identical only change is the sign of E and ε. The parameter E in the Green's function can either be the energy if the paths are going toward the future, or the negative of the energy if the paths are going toward the past.

For a nonrelativistic theory, the time as measured along the path of a moving particle and the time as measured by an outside observer are the same. In relativity, this is no longer true. For a relativistic theory the propagator should be defined as the sum over all paths which travel between two points in a fixed proper time, as measured along the path. These paths describe the trajectory of a particle in space and in time.

Failed to parse (Missing texvc executable; please see math/README to configure.): K(x-y,\Tau) = \int_{x(0)=x}^{x(\Tau)=y} e^{i \int_0^\Tau \sqrt{{\dot x}^2} - \alpha d\tau} \,

The integral above is not trivial to interpret, because of the square root. Fortunately, there is a heuristic trick. The sum is over the relativistic arclength of the path of an oscillating quantity, and like the nonrelativistic path integral should be interpreted as slightly rotated into imaginary time. The function K(x-y,\tau) can be evaluated when the sum is over paths in Euclidean space.

Failed to parse (Missing texvc executable; please see math/README to configure.): K(x-y,\Tau) = e^{-\alpha \Tau} \int_{x(0)=x}^{x(\Tau)=y} e^{-L} \,

This describes a sum over all paths of length Failed to parse (Missing texvc executable; please see math/README to configure.): \Tau

of the exponential of minus the length. This can be given a probability interpretation. The sum over all paths is a probability average over a path constructed step by step. The total number of steps is proportional to Failed to parse (Missing texvc executable; please see math/README to configure.): \Tau

, and each step is less likely the longer it is. By the central limit theorem, the result of many independent steps is a Gaussian of variance proportional to Failed to parse (Missing texvc executable; please see math/README to configure.): \Tau .

Failed to parse (Missing texvc executable; please see math/README to configure.): K(x-y,\Tau) = e^{-\alpha \Tau} e^{-(x-y)^2\over \Tau} \,

The usual definition of the relativistic propagator only asks for the amplitude is to travel from x to y, after summing over all the possible proper times it could take.

Failed to parse (Missing texvc executable; please see math/README to configure.): K(x-y) = \int_0^{\infty} K(x-y,\Tau) W(\Tau) d\Tau \,

Where Failed to parse (Missing texvc executable; please see math/README to configure.): W(\Tau)

is a weight factor, the relative importance of paths of different proper time. By the translation symmetry in proper time, this weight can only be an exponential factor, and can be absorbed into the constant Failed to parse (Missing texvc executable; please see math/README to configure.): \alpha

.

Failed to parse (Missing texvc executable; please see math/README to configure.): K(x-y) = \int_0^{\infty} e^{-{(x-y)^2\over\Tau} -\alpha \Tau} d\Tau \,

This is the Schwinger representation. Taking a Fourier transform over the variable Failed to parse (Missing texvc executable; please see math/README to configure.): x-y

can be done for each value of Failed to parse (Missing texvc executable; please see math/README to configure.): \Tau
separately, and because each separate Failed to parse (Missing texvc executable; please see math/README to configure.): \Tau
contribution is a Gaussian, gives whose fourier transform is another Gaussian with reciprocal width. So in p-space, the propagator can be reexpressed simply:

Failed to parse (Missing texvc executable; please see math/README to configure.): K(p) = \int_0^{\infty} e^{-\Tau p^2 - \Tau \alpha} d\Tau = {1\over p^2 + \alpha } \,

Which is the Euclidian propagator for a scalar particle. Rotating Failed to parse (Missing texvc executable; please see math/README to configure.): p_0

to be imaginary gives the usual relativistic propagator, up to a -i and an ambiguity which will be clarified below.

Failed to parse (Missing texvc executable; please see math/README to configure.): K(p) = {i\over p_0^2 - \vec{p}^2 - m^2} \,

This expression can be interpreted in the nonrelativistic limit, where it is convenient to split it by partial fractions:

Failed to parse (Missing texvc executable; please see math/README to configure.): 2 p_0 K(p) = {i \over p_0 - \sqrt{\vec{p}^2 + m^2}} + {i \over p_0 + \sqrt{\vec{p}^2 + m^2}}

For states where one nonrelativistic particle is present, the initial wavefunction has a frequency distribution concentrated near Failed to parse (Missing texvc executable; please see math/README to configure.): p_0 = m . When convolving with the propagator, which in p space just means multiplying by the propagator, the second term is suppressed and the first term is enhanced. For frequencies near Failed to parse (Missing texvc executable; please see math/README to configure.): p_0=m , the dominant first term has the form:

Failed to parse (Missing texvc executable; please see math/README to configure.): 2m K_\mathrm{NR}(p) = {i \over (p_0-m) - {\vec{p}^2\over 2m} }

This is the expression for the Nonrelativistic Green's function of a free Schrödinger particle.

The second term has a nonrelativistic limit also, but this limit is concentrated on frequencies which are negative. The second pole is dominated by contributions from paths where the proper time and the coordinate time are ticking in an opposite sense, which means that the second term is to be interpreted as the antiparticle. The nonrelativistic analysis shows that with this form the antiparticle still has positive energy.

The proper way to express this mathematically is that, adding a small suppression factor in proper time, the limit where Failed to parse (Missing texvc executable; please see math/README to configure.): t \rightarrow -\infty

of the first term must vanish, while the Failed to parse (Missing texvc executable; please see math/README to configure.): t\rightarrow \infty
limit of the second term must vanish. In the fourier transform, this means shifting the pole in Failed to parse (Missing texvc executable; please see math/README to configure.): p_0
slightly, so that the inverse fourier transform will pick up a small decay factor in one of the time directions:

Failed to parse (Missing texvc executable; please see math/README to configure.): K(p) = {i \over p_0 - \sqrt{\vec{p}^2 + m^2} + i\epsilon} + {i \over p_0 - \sqrt{\vec{p}^2+m^2} - i\epsilon}

Without these terms, the pole contribution could not be unambiguously evaluated when taking the inverse Fourier transform of Failed to parse (Missing texvc executable; please see math/README to configure.): p_0 . The terms can be recombined:

Failed to parse (Missing texvc executable; please see math/README to configure.): K(p) = { i \over {p^2 - m^2 + i\epsilon}}

Which when factored, produces opposite sign infinitesimal terms in each factor. This is the mathematically precise form of the relativistic particle propagator, free of any ambiguities. The Failed to parse (Missing texvc executable; please see math/README to configure.): \epsilon

term introduces a small imaginary part to the Failed to parse (Missing texvc executable; please see math/README to configure.): \alpha=m^2

, which in the Minkowski version is a small exponential suppresion of long paths.

So in the relativistic case, the Feynman path-integral representation of the propagator includes paths which go backwards in time, which describe antiparticles. The paths which contribute to the relativistic propagator go forward and backwards in time, and the interpretation of this is that the amplitude for a free particle to travel between two points includes amplitudes for the particle to fluctuate into an antiparticle, travel back in time, then forward again.

Unlike the nonrelativistic case, it is impossible to produce a relativistic theory of local particle propagation without including antiparticles. All local differential operators have inverses which are nonzero outside the lightcone, meaning that it is impossible to keep a particle from travelling faster than light. Such a particle cannot be have a Greens function which is only nonzero in the future in a relativistically invariant theory.

Functionals of fields

However, the path integral formulation is also extremely important in direct application to quantum field theory, in which the "paths" or histories being considered are not the motions of a single particle, but the possible time evolutions of a field over all space. The action is referred to technically as a functional of the field: Failed to parse (Missing texvc executable; please see math/README to configure.): S[\phi] \,

where the field Failed to parse (Missing texvc executable; please see math/README to configure.): \phi (x^\mu) \,
is itself a function of space and time, and the square brackets are a reminder that the action depends on all the field's values everywhere, not just some particular value.  In principle, one integrates Feynman's amplitude over the class of all possible combinations of values that the field could have anywhere in space-time.

Much of the formal study of QFT is devoted to the properties of the resulting functional integral, and much effort (not yet entirely successful) has been made toward making these functional integrals mathematically precise.

Such a functional integral is extremely similar to the partition function in statistical mechanics. Indeed, it is sometimes called a partition function, and the two are essentially mathematically identical except for the factor of Failed to parse (Missing texvc executable; please see math/README to configure.): i

in the exponent in Feynman's postulate 3.  Analytically continuing the integral to an imaginary time variable (called a Wick rotation) makes the functional integral even more like a statistical partition function, and also tames some of the mathematical difficulties of working with these integrals.

Expectation values

In quantum field theory, if the action is given by the functional Failed to parse (Missing texvc executable; please see math/README to configure.): \mathcal{S}

of field configurations (which only depends locally on the fields), then the time ordered vacuum expectation value of polynomially bounded functional F, <F>, is given by

Failed to parse (Missing texvc executable; please see math/README to configure.): \left\langle F\right\rangle=\frac{\int \mathcal{D}\phi F[\phi]e^{i\mathcal{S}[\phi]}}{\int\mathcal{D}\phi e^{i\mathcal{S}[\phi]}}

The symbol Failed to parse (Missing texvc executable; please see math/README to configure.): \int \mathcal{D}\phi

here is a concise way to represent the infinite-dimensional integral over all possible field configurations on all of space-time.  As stated above, we put the unadorned path integral in the denominator to normalize everything properly.

Schwinger-Dyson equations

Since this formulation of quantum mechanics is analogous to classical action principles, one might expect that identities concerning the action in classical mechanics would have quantum counterparts derivable from a functional integral. This is often the case.

In the language of functional analysis, we can write the Euler-Lagrange equations as Failed to parse (Missing texvc executable; please see math/README to configure.): \frac{\delta \mathcal{S}[\phi]}{\delta \phi}=0

(the left-hand side is a functional derivative; the equation means that the action is stationary under small changes in the field configuration).  The quantum analogues of these equations are called the Schwinger-Dyson equations.

If the functional measure Failed to parse (Missing texvc executable; please see math/README to configure.): \mathcal{D}\phi

turns out to be translationally invariant (we'll assume this for the rest of this article, although this does not hold for, let's say nonlinear sigma models) and if we assume that after a Wick rotation 

Failed to parse (Missing texvc executable; please see math/README to configure.): e^{i\mathcal{S}[\phi]},

which now becomes

Failed to parse (Missing texvc executable; please see math/README to configure.): e^{-H[\phi]}\,

for some H, goes to zero faster than any reciprocal of any polynomial for large values of φ, we can integrate by parts (after a Wick rotation, followed by a Wick rotation back) to get the following Schwinger-Dyson equations:

Failed to parse (Missing texvc executable; please see math/README to configure.): \left\langle \frac{\delta F[\phi]}{\delta \phi} \right\rangle = -i \left\langle F[\phi]\frac{\delta \mathcal{S}[\phi]}{\delta\phi} \right\rangle

for any polynomially bounded functional F.

Failed to parse (Missing texvc executable; please see math/README to configure.): \left\langle F_{,i} \right\rangle = -i \left\langle F \mathcal{S}_{,i} \right\rangle

in the deWitt notation.

These equations are the analog of the on shell EL equations.

If J (called the source field) is an element of the dual space of the field configurations (which has at least an affine structure because of the assumption of the translational invariance for the functional measure), then, the generating functional Z of the source fields is defined to be:

Failed to parse (Missing texvc executable; please see math/README to configure.): Z[J]=\int \mathcal{D}\phi e^{i(\mathcal{S}[\phi] + \left\langle J,\phi \right\rangle)}.

Note that

Failed to parse (Missing texvc executable; please see math/README to configure.): \frac{\delta^n Z}{\delta J(x_1) \cdots \delta J(x_n)}[J] = i^n \, Z[J] \, {\left\langle \phi(x_1)\cdots \phi(x_n)\right\rangle}_J

or

Failed to parse (Missing texvc executable; please see math/README to configure.): Z^{,i_1\dots i_n}[J]=i^n Z[J] {\left \langle \phi^{i_1}\cdots \phi^{i_n}\right\rangle}_J

where

Failed to parse (Missing texvc executable; please see math/README to configure.): {\left\langle F \right\rangle}_J=\frac{\int \mathcal{D}\phi F[\phi]e^{i(\mathcal{S}[\phi] + \left\langle J,\phi \right\rangle)}}{\int\mathcal{D}\phi e^{i(\mathcal{S}[\phi] + \left\langle J,\phi \right\rangle)}}.

Basically, if Failed to parse (Missing texvc executable; please see math/README to configure.): \mathcal{D}\phi e^{i\mathcal{S}[\phi]}

is viewed as a functional distribution (this shouldn't be taken too literally as an interpretation of QFT, unlike its Wick rotated statistical mechanics analogue, because we have time ordering complications here!), then Failed to parse (Missing texvc executable; please see math/README to configure.): \left\langle\phi(x_1)\cdots \phi(x_n)\right\rangle
are its moments and Z is its Fourier transform.

If F is a functional of φ, then for an operator K, F[K] is defined to be the operator which substitutes K for φ. For example, if

Failed to parse (Missing texvc executable; please see math/README to configure.): F[\phi]=\frac{\partial^{k_1}}{\partial x_1^{k_1}}\phi(x_1)\cdots \frac{\partial^{k_n}}{\partial x_n^{k_n}}\phi(x_n)

and G is a functional of J, then

Failed to parse (Missing texvc executable; please see math/README to configure.): F\left[-i\frac{\delta}{\delta J}\right] G[J] = (-i)^n \frac{\partial^{k_1}}{\partial x_1^{k_1}}\frac{\delta}{\delta J(x_1)} \cdots \frac{\partial^{k_n}}{\partial x_n^{k_n}}\frac{\delta}{\delta J(x_n)} G[J].

Then, from the properties of the functional integrals, we get the "master" Schwinger-Dyson equation:

Failed to parse (Missing texvc executable; please see math/README to configure.): \frac{\delta \mathcal{S}}{\delta \phi(x)}\left[-i \frac{\delta}{\delta J}\right]Z[J]+J(x)Z[J]=0

or

Failed to parse (Missing texvc executable; please see math/README to configure.): \mathcal{S}_{,i}[-i\partial]Z+J_i Z=0.

If the functional measure is not translationally invariant, it might be possible to express it as the product Failed to parse (Missing texvc executable; please see math/README to configure.): M\left[\phi\right]\,\mathcal{D}\phi

where M is a functional and Failed to parse (Missing texvc executable; please see math/README to configure.): \mathcal{D}\phi
is a translationally invariant measure. This is true, for example, for nonlinear sigma models where the target space is diffeomorphic to Rn. However, if the target manifold is some topologically nontrivial space, the concept of a translation does not even make any sense.

In that case, we would have to replace the Failed to parse (Missing texvc executable; please see math/README to configure.): \mathcal{S}

in this equation by another functional Failed to parse (Missing texvc executable; please see math/README to configure.): \hat{\mathcal{S}}=\mathcal{S}-i\ln(M)

If we expand this equation as a Taylor series about J=0, we get the entire set of Schwinger-Dyson equations.

Functional identity

If we perform a Wick rotation inside the functional integral, professors J. Garcia and Gerard 't Hooft showed using a functional differential equation that:

Failed to parse (Missing texvc executable; please see math/README to configure.): \int D[x]e^{-\mathcal{S}[x]/\hbar}=-A[x]\sum_{n=0}^{\infty}(\hbar)^{n+1}\delta^{n} e^{-J/\hbar}

where S is the Wick-rotated classical action of the particle, J is the classical action with an extra term "x" and delta here is the functional derivative operator

Failed to parse (Missing texvc executable; please see math/README to configure.): A[x]=\exp\left({1/\hbar}\int X(t)\,\mathrm{d}t\right).

Ward-Takahashi identities

See main article Ward-Takahashi identity

Now how about the on shell Noether's theorem for the classical case? Does it have a quantum analog as well? Yes, but with a caveat. The functional measure would have to be invariant under the one parameter group of symmetry transformation as well.

Let's just assume for simplicity here that the symmetry in question is local (not local in the sense of a gauge symmetry, but in the sense that the transformed value of the field at any given point under an infinitesimal transformation would only depend on the field configuration over an arbitrarily small neighborhood of the point in question). Let's also assume that the action is local in the sense that it is the integral over spacetime of a Lagrangian, and that Failed to parse (Missing texvc executable; please see math/README to configure.): Q[\mathcal{L}(x)]=\partial_\mu f^\mu (x)

for some function f where f only depends locally on φ (and possibly the spacetime position).

If we don't assume any special boundary conditions, this would not be a "true" symmetry in the true sense of the term in general unless f=0 or something. Here, Q is a derivation which generates the one parameter group in question. We could have antiderivations as well, such as BRST and supersymmetry.

Let's also assume Failed to parse (Missing texvc executable; please see math/README to configure.): \int \mathcal{D}\phi Q[F][\phi]=0

for any polynomially bounded functional F. This property is called the invariance of the measure. And this does not hold in general. See anomaly (physics) for more details.

Then,

Failed to parse (Missing texvc executable; please see math/README to configure.): \int \mathcal{D}\phi\, Q[F e^{iS}][\phi]=0,

which implies

Failed to parse (Missing texvc executable; please see math/README to configure.): \left\langle Q[F]\right\rangle +i\left\langle F\int_{\partial V} f^\mu ds_\mu\right\rangle=0

where the integral is over the boundary. This is the quantum analog of Noether's theorem.

Now, let's assume even further that Q is a local integral

Failed to parse (Missing texvc executable; please see math/README to configure.): Q=\int d^dx q(x)

where

Failed to parse (Missing texvc executable; please see math/README to configure.): q(x)[\phi(y)] = \delta^{(d)}(X-y)Q[\phi(y)] \,

so that

Failed to parse (Missing texvc executable; please see math/README to configure.): q(x)[S]=\partial_\mu j^\mu (x) \,

where

Failed to parse (Missing texvc executable; please see math/README to configure.): j^{\mu}(x)=f^\mu(x)-\frac{\partial}{\partial (\partial_\mu \phi)}\mathcal{L}(x) Q[\phi] \,

(this is assuming the Lagrangian only depends on φ and its first partial derivatives! More general Lagrangians would require a modification to this definition!). Note that we're NOT insisting that q(x) is the generator of a symmetry (i.e. we are not insisting upon the gauge principle), but just that Q is. And we also assume the even stronger assumption that the functional measure is locally invariant:

Failed to parse (Missing texvc executable; please see math/README to configure.): \int \mathcal{D}\phi\, q(x)[F][\phi]=0.

Then, we would have

Failed to parse (Missing texvc executable; please see math/README to configure.): \left\langle q(x)[F] \right\rangle +i\left\langle F q(x)[S]\right\rangle=\left\langle q(x)[F]\right\rangle +i\left\langle F\partial_\mu j^\mu(x)\right\rangle=0.

Alternatively,

Failed to parse (Missing texvc executable; please see math/README to configure.): q(x)[S][-i \frac{\delta}{\delta J}]Z[J]+J(x)Q[\phi(x)][-i \frac{\delta}{\delta J}]Z[J]=\partial_\mu j^\mu(x)[-i \frac{\delta}{\delta J}]Z[J]+J(x)Q[\phi(x)][-i \frac{\delta}{\delta J}]Z[J]=0.

The above two equations are the Ward-Takahashi identities.

Now for the case where f=0, we can forget about all the boundary conditions and locality assumptions. We'd simply have

Failed to parse (Missing texvc executable; please see math/README to configure.): \left\langle Q[F]\right\rangle =0.

Alternatively,

Failed to parse (Missing texvc executable; please see math/README to configure.): \int d^dx\, J(x)Q[\phi(x)][-i \frac{\delta}{\delta J}]Z[J]=0.

The path integral in quantum-mechanical interpretation

In one philosophical interpretation of quantum mechanics, the "sum over histories" interpretation, the path integral is taken to be fundamental and reality is viewed as a single indistinguishable "class" of paths which all share the same events. For this interpretation, it is crucial to understand what exactly an event is. The sum over histories method gives identical results to canonical quantum mechanics, and Sinha and Sorkin (see the reference below) claim the interpretation explains the Einstein-Podolsky-Rosen paradox without resorting to nonlocality.

Some advocates of interpretations of quantum mechanics emphasizing decoherence have attempted to make more rigorous the notion of extracting a classical-like "coarse-grained" history from the space of all possible histories.

References

See also

• Theoretical and experimental justification for the Schrödinger equation
• Feynman checkerboard

Suggested reading

• Feynman, R. P., and Hibbs, A. R., Quantum Physics and Path Integrals, New York: McGraw-Hill, 1965 [ISBN 0-07-020650-3]. The historical reference, written by the inventor of the path integral formulation himself and one of his students.
• Hagen Kleinert, Path Integrals in Quantum Mechanics, Statistics, Polymer Physics, and Financial Markets, 4th edition, World Scientific (Singapore, 2004); Paperback ISBN 981-238-107-4 (also available online: PDF-files)
• Zinn Justin, Jean ; Path Integrals in Quantum Mechanics, Oxford University Press (2004), [ISBN 0-19-856674-3]. A highly readable introduction to the subject.
• Schulman, Larry S.  ; Techniques & Applications of Path Integration, John Wiley & Sons (New York-1981) [ISBN]. A modern reference on the subject.
• Grosche, Christian & Steiner, Frank ; Handbook of Feynman Path Integrals, Springer Tracts in Modern Physics 145, Springer-Verlag (1998) [ISBN 3-540-57135-3]
• Ryder, Lewis H. ; Quantum Field Theory (Cambridge University Press, 1985), [ISBN 0-521-33859-X] Highly readable textbook; introduction to relativistic Q.F.T. for particle physics.
• Rivers, R.J. ; Path Integrals Methods in Quantum Field Theory, Cambridge University Press (1987) [ISBN 0-521-25979-7]
• Albeverio, S. & Hoegh-Krohn. R. ; Mathematical Theory of Feynman Path Integral, Lecture Notes in Mathematics 523, Springer-Verlag (1976) [ISBN].
• Glimm, James, and Jaffe, Arthur, Quantum Physics: A Functional Integral Point of View, New York: Springer-Verlag, 1981. [ISBN 0-387-90562-6].
• Gerald W. Johnson and Michel L. Lapidus ; The Feynman Integral and Feynman's Operational Calculus, Oxford Mathematical Monographs, Oxford University Press (2002) [ISBN 0-19-851572-3].
• Etingof, Pavel ; Geometry and Quantum Field Theory, M.I.T. OpenCourseWare (2002). This course, designed for mathematicians, is a rigorous introduction to perturbative quantum field theory, using the language of functional integrals.

Papers on-line

• Grosche, Christian  ; An Introduction into the Feynman Path Integral, lecture given at the graduate college Quantenfeldtheorie und deren Anwendung in der Elementarteilchen- und Festkörperphysik, Universität Leipzig, 16-26 November 1992. Full text available at : hep-th/9302097.
• MacKenzie, Richard ; Path Integral Methods and Applications, lectures given at Rencontres du Vietnam: VIth Vietnam School of Physics, Vung Tau, Vietnam, 27 December 1999 - 8 January 2000. Full text available at : quant-ph/0004090.
• DeWitt-Morette, Cécile ; Feynman's path integral - Definition without limiting procedure, Communication in Mathematical Physics 28(1) (1972) pp. 47–67. Full text available at : Euclide Project.
• Sukanya Sinha and Rafael D. Sorkin, "A Sum-over-histories Account of an EPR(B) Experiment", Found. of Phys. Lett. 4:303-335 (1991). Full text available at :"Sinha-Sorkin 1991".
• Cartier, Pierre & DeWitt-Morette, Cécile ; A new perspective on Functional Integration, Journal of Mathematical Physics 36 (1995) pp. 2137-2340. Full text available at : funct-an/9602005.de:Pfadintegral

fr:Intégrale de chemin it:Integrale sui cammini he:אינטגרלי מסלול nl:Padintegraal ja:経路積分 pt:Integração funcional fi:Feynmanin polkuintegraali