This article is about a formulation of quantum mechanics. For integrals along a path, also known as line or contour integrals, see
line integral.
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 trajectory for a system with a sum, or functional integral, over an infinity of possible trajectories to compute a quantum amplitude.
The basic idea of the path integral formulation can be traced back to Norbert Wiener, who introduced the Wiener integral for solving problems in diffusion and Brownian motion.^{[1]} This idea was extended to the use of the Lagrangian in quantum mechanics by P. A. M. Dirac in his 1933 paper.^{[2]} The complete method 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. The original motivation stemmed from the desire to obtain a quantummechanical formulation for the Wheeler–Feynman absorber theory using a Lagrangian (rather than a Hamiltonian) as a starting point.
This formulation has proven crucial to the subsequent development of theoretical physics, because it is manifestly symmetric between time and space. Unlike previous methods, the pathintegral allows a physicist to easily change coordinates between very different canonical descriptions of the same quantum system.
The path integral also relates quantum and stochastic processes, and this provided the basis for the grand synthesis of the 1970s which unified quantum field theory with the statistical field theory of a fluctuating field near a secondorder phase transition. The Schrödinger equation is a diffusion equation with an imaginary diffusion constant, and the path integral is an analytic continuation of a method for summing up all possible random walks. For this reason path integrals were used in the study of Brownian motion and diffusion a while before they were introduced in quantum mechanics.^{[3]}
Quantum action principle
In quantum mechanics, as in classical mechanics, the Hamiltonian is the generator of timetranslations. This means that the state at a slightly later time differs from the state at the current time by the result of acting with the Hamiltonian operator (multiplied by the negative imaginary unit, −i). For states with a definite energy, this is a statement of the De Broglie relation between frequency and energy, and the general relation is consistent with that plus the superposition principle.
But the Hamiltonian in classical mechanics is derived from a Lagrangian, which is a more fundamental quantity considering special relativity. The Hamiltonian tells you how to march forward in time, but the time is different in different reference frames. So the Hamiltonian is different in different frames, and this type of symmetry is not apparent in the original formulation of quantum mechanics.
The Hamiltonian is a function of the position and momentum at one time, and it tells you the position and momentum a little later. The Lagrangian is a function of the position now and the position a little later (or, equivalently for infinitesimal time separations, it is a function of the position and velocity). The relation between the two is by a Legendre transform, and the condition that determines the classical equations of motion (the Euler–Lagrange equations) is that the action is a minimum.
In quantum mechanics, the Legendre transform is hard to interpret, because the motion is not over a definite trajectory. So what does the Legendre transform mean? In classical mechanics, with discretization in time,
 $\backslash epsilon\; H\; =\; p\; (q(t+\backslash epsilon)\; \; q(t))\; \; \backslash epsilon\; L\; \backslash ,$
and
 $p\; =\; \{\backslash partial\; L\; \backslash over\; \backslash partial\; \backslash dot\{q\}\; \}\; \backslash ,$
where the partial derivative with respect to q holds q(t + ε) fixed. The inverse Legendre transform is:
 $\backslash epsilon\; L\; =\; p\; \backslash epsilon\; \backslash dot\{q\}\; \; \backslash epsilon\; H\; \backslash ,$
where
 $\backslash dot\; q\; =\; \{\backslash partial\; H\; \backslash over\; \backslash partial\; p\}\; \backslash ,$
and the partial derivative now is with respect to p at fixed q.
In quantum mechanics, the state is a superposition of different states with different values of q, or different values of p, and the quantities p and q can be interpreted as noncommuting operators. The operator p is only definite on states that are indefinite with respect to q. So consider two states separated in time and act with the operator corresponding to the Lagrangian:
 $e^\{i(\; p\; (q(t+\backslash epsilon)\; \; q(t))\; \; \backslash epsilon\; H(p,q)\; )\}\backslash ,$
If the multiplications implicit in this formula are reinterpreted as matrix multiplications, what does this mean?
It can be given a meaning as follows: The first factor is
 $e^\{ip\; q(t)\}\; \backslash ,$
If this is interpreted as doing a matrix multiplication, the sum over all states integrates over all q(t), and so it takes the Fourier transform in q(t), to change basis to p(t). That is the action on the Hilbert space – change basis to p at time t.
Next comes:
 $e^\{i\backslash epsilon\; H(p,q)\}\; \backslash ,$
or evolve an infinitesimal time into the future.
Finally, the last factor in this interpretation is
 $e^\{i\; p\; q(t+\backslash epsilon)\}\; \backslash ,$
which means change basis back to q at a later time.
This is not very different from just ordinary time evolution: the H factor contains all the dynamical information – it pushes the state forward in time. The first part and the last part are just doing Fourier transforms to change to a pure q basis from an intermediate p basis.
Another way of saying this is that since the Hamiltonian is naturally a function of p and q, exponentiating this quantity and changing basis from p to q at each step allows the matrix element of H to be expressed as a simple function along each path. This function is the quantum analog of the classical action. This observation is due to Paul Dirac.
"...we see that the integrand in (11) must be of the form e^{iF/h} where F is a function of q_{T},q_{1},q_{2} ... q_{m},q_{t}, which remains finite as h tends to zero. Let us now picture one of the intermediate qs, say q_{k}, as varying continuously while the other ones are fixed. Owing to the smallness of h, we shall then in general have F/h varying extremely rapidly. This means that e^{iF/h} will vary periodically with a very high frequency about the value zero, as a result of which its integral will be practically zero. The only important part in the domain of integration of q_{k} is thus that for which a comparatively large variation in q_{k} produces only a very small variation in F. This part is the neighbourhood of a point for which F is stationary with respect to small variations in q_{k}. We can apply this argument to each of the variables of integration ....and obtain the result that the only important part in the domain of integration is that for which F is stationary for small variations in all intermediate qs. ...We see that F has for its classical analogue Template:Intmath L dt, which is just the action function which classical mechanics requires to be stationary for small variations in all the intermediate qs. This shows the way in which equation (11) goes over into classical results when h becomes extremely small."
Dirac (1932) op. cit., p. 69
Dirac further noted that one could square the timeevolution operator in the S representation
 $e^\{i\backslash epsilon\; S\}\; \backslash ,$
and this gives the time evolution operator between time t and time t + 2ε. While in the H representation the quantity that is being summed over the intermediate states is an obscure matrix element, in the S representation it is reinterpreted as a quantity associated to the path. In the limit that one takes a large power of this operator, one reconstructs the full quantum evolution between two states, the early one with a fixed value of q(0) and the later one with a fixed value of q(t). The result is a sum over paths with a phase which is the quantum action. Crucially, Dirac identified in this paper the deep quantum mechanical reason for the principle of least action controlling the classical limit (see quote box).
Feynman's interpretation
Dirac's work did not provide a precise prescription to calculate the sum over paths, 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.^{[4]}
Feynman showed that Dirac's quantum action was, for most cases of interest, simply equal to the classical action, appropriately discretized. This means that the classical action is the phase acquired by quantum evolution between two fixed endpoints. He proposed to recover all of quantum mechanics from the following postulates:
 The probability for an event is given by the squared length of a complex number called the "probability amplitude".
 The probability amplitude is given by adding together the contributions of all paths in configuration space.
 The contribution of a path is proportional to $e^\{i\; S/\backslash hbar\}$. while S is the action given by the time integral of the Lagrangian along the path.
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 paths of the system in between the initial and final states, including those 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 is correct to include paths in which the particle describes elaborate curlicues, curves in which the particle shoots off into outer space and flies back again, and so forth. The path integral assigns to all these amplitudes equal magnitude but varying phase, or argument of the complex number. Contributions from paths wildly different from the classical trajectory may be suppressed by interference (see below).
Feynman showed that this formulation of quantum mechanics is equivalent to the canonical approach to quantum mechanics when the Hamiltonian is quadratic in the momentum. An amplitude computed according to Feynman's principles will also obey the Schrödinger equation for the Hamiltonian corresponding to the given action.
Concrete formulation
Feynman's postulates can be interpreted as follows:
Timeslicing definition
For a particle in a smooth potential, the path integral is approximated by zigzag paths, which in one dimension is a product of ordinary integrals. For the motion of the particle from position x_{a} at time t_{a} to x_{b} at time t_{b}, the time sequence
 $t\_a=t\_0...\{n1\}=t\_\{n+1\}\; math>$
can be divided up into n + 1 little segments t_{j} − t_{j − 1}, where j = 1,...,n + 1, of fixed duration
 $\backslash epsilon\; =\; \backslash Delta\; t=\backslash tfrac\{t\_bt\_a\}\{n+1\}\backslash ,.$
This process is called timeslicing.
An approximation for the path integral can be computed as proportional to
 $\backslash int\backslash limits\_\{\backslash infty\}^\{+\backslash infty\}\backslash ,\backslash ldots\; \backslash int\backslash limits\_\{\backslash infty\}^\{+\backslash infty\}\backslash ,$
\ \exp \left(\fracTemplate:\rm i{\hbar}\int\limits_{t_a}^{t_b} L(x(t),v(t), t)\,\mathrm{d}t\right)dx_0 \, \ldots \, dx_n
where $L(x,v,t)$ is the Lagrangian of the 1d system with position variable x(t) and velocity v = ẋ(t) considered (see below), and dx_{j} corresponds to the position at the jth time step, if the time integral is approximated by a sum of n terms.^{[note 1]}
In the limit n → ∞, this becomes a functional integral, which – apart from a nonessential factor – is directly the product of the probability amplitudes $\backslash langle\; x\_a,t\_ax\_b,\; t\_b\backslash rangle$ – more precisely, since one must work with a continuous spectrum, the respective densities – to find the quantum mechanical particle at t_{a} in the initial state x_{a} and at t_{b} in the final state x_{b}.
Actually $L$ is the classical Lagrangian of the onedimensional system considered, also
 $L(x,\backslash dot\; x\; ,\; t)=p\backslash cdot\; \backslash dot\; x\; \; H(x,p,t)\backslash ,,$
where $H$ is the Hamiltonian,
 $p=\backslash frac\; \{\backslash partial\; L\}\{\backslash partial\; \backslash dot\; x\}$, and the abovementioned "zigzagging" corresponds to the appearance of the terms:
 $\backslash exp\backslash left\; (\backslash frac$Template:\rm i{\hbar}\epsilon\, \,\sum_{j=1}^{n} L \left (\tilde x_{j},\frac{x_jx_{j1}}{\epsilon},j \right )\right )
In the Riemannian sum approximating the time integral, which are finally integrated over x_{1} to x_{n} with the integration measure dx_{1}...dx_{n}x̃_{j} is an arbitrary value of the interval corresponding to j, e.g. its center, (x_{j} + x_{}j − 1)/2.
Thus, in contrast to classical mechanics, not only does the stationary path contribute, but actually all virtual paths between the initial and the final point also contribute.
File:Path integral example.webm
Feynman's timesliced approximation does not, however, exist for the most important quantummechanical path integrals of atoms, due to the singularity of the Coulomb potential e^{2}/r at the origin. Only after replacing the time t by another pathdependent pseudotime parameter
 $s=\backslash int\; \backslash frac\{dt\}\{r(t)\}$
the singularity is removed and a timesliced approximation exists, that is exactly integrable, since it can be made harmonic by a simple coordinate transformation, as discovered in 1979 by İsmail Hakkı Duru and Hagen Kleinert.^{[5]}^{[6]} The combination of a pathdependent time transformation and a coordinate transformation is an important tool to solve many path integrals and is called generically the 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 (m = 1, ħ = 1):
 $$
S= \int {\dot{x}^2\over 2} dt
\,
the integral can be evaluated explicitly.
To do this, it is convenient to start without the factor i in the exponential, so that large deviations are suppressed by small numbers, not by cancelling oscillatory contributions.
 $$
K(xy;T) = \int_{x(0)=x}^{x(T)=y} \exp\left\{\int_0^T {\dot{x}^2\over 2} dt\right\} Dx
\,
Splitting the integral into time slices:
 $$
K(x,y;T) = \int_{x(0)=x}^{x(T)=y} \Pi_t \exp\left\{{1\over 2} \left({x(t+\epsilon) x(t) \over \epsilon}\right)^2 \epsilon \right\} Dx
\,
where the Dx is interpreted as a finite collection of integrations at each integer multiple of ε. Each factor in the product is a Gaussian as a function of x(t + ε) centered at x(t) with variance ε. The multiple integrals are a repeated convolution of this Gaussian G_{ε} with copies of itself at adjacent times.
 $$
K(xy;T) = G_\epsilon*G_\epsilon ... *G_\epsilon
\,
Where the number of convolutions is T/ε. The result is easy to evaluate by taking the fourier transform of both sides, so that the convolutions become multiplications.
 $$
\tilde{K}(p;T) = \tilde{G}_\epsilon(p)^{T/\epsilon}
\,
The Fourier transform of the Gaussian G is another Gaussian of reciprocal variance:
 $$
\tilde{G}_\epsilon(p) = e^{\epsilon {p^2/2} }
\,
and the result is:
 $$
\tilde{K}(p;T) = e^{T {p^2/2}}
\,
The Fourier transform gives K, and it is a Gaussian again with reciprocal variance:
 $$
K(xy;T) \propto e^{ {(xy)^2/(2T)}}
\,
The proportionality constant is not really determined by the time slicing approach, only the ratio of values for different endpoint choices is determined. The proportionality constant should be chosen to ensure that between each two timeslices the timeevolution is quantummechanically unitary, but a more illuminating way to fix the normalization is to consider the path integral as a description of a stochastic process.
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:
 $$
\int K(xy;T) dy = 1
\,
This condition normalizes the Gaussian, and produces a Kernel which obeys the diffusion equation:
 $$
{d\over dt} K(x;T) = {\nabla^2 \over 2} K
\,
For oscillatory path integrals, ones with an i in the numerator, the timeslicing 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 $\backslash epsilon$. This is closely related to Wick rotation. Then the same convolution argument as before gives the propagation kernel:
 $$
K(xy;T) \propto e^{i(xy)^2 / (2T)}
\,
Which, with the same normalization as before (not the sumsquares normalization – this function has a divergent norm), obeys a free Schrödinger equation
 $$
{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
 $$
\psi_t(y) = \int \psi_0(x) K(xy;t) dx = \int \psi_0(x) \int_{x(0)=x}^{x(t)=y} e^{iS} Dx
\,
then ψ_{t} obeys the free Schrödinger equation just as K does:
 $$
{\rm i}{\partial \over \partial t} \psi_t =  {\nabla^2\over 2} \psi_t
\,
The 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 pathintegral over infinitesimally separated times.
 $$
\psi(y;t+\epsilon) = \int_{\infty}^\infty \;\;\psi(x;t)\int_{x(t)=x}^{x(t+\epsilon)=y} e^  {\partial S \over \partial u}\right)\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 ε(t). The conclusion is that this first order variation is zero for an arbitrary initial state and at any arbitrary point in time:
 $$
\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 ẋ and 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.
Stationary phase approximation
If the variation in the action exceeds ħ by many orders of magnitude, we typically have destructive phase interference other than in the vicinity of those trajectories satisfying the Euler–Lagrange equation, which is now reinterpreted as the condition for constructive phase interference.
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 noncommutativity is still there.^{[7]}
To see this, consider the simplest path integral, the brownian walk. This is not yet quantum mechanics, so in the pathintegral the action is not multiplied by i:
 $$
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.
 $$
{dx \over dt} = {x(t+\epsilon)  x(t) \over \epsilon}
\,
Note that the distance that a random walk moves is proportional to √t, so that:
 $$
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 x ẋ is ambiguous, with two possible meanings:
 $[1]\; =\; x\; \{\; dx\backslash over\; dt\}\; =\; x(t)\; \{(x(t+\backslash epsilon)\; \; x(t))\; \backslash over\; \backslash epsilon\; \}\; \backslash ,$
 $[2]\; =\; x\; \{dx\; \backslash over\; dt\}\; =\; x(t+\backslash epsilon)\; \{(x(t+\backslash epsilon)\; \; x(t))\; \backslash over\; \backslash epsilon\}\; \backslash ,$
In elementary calculus, the two are only different by an amount which goes to zero as ε goes to zero. But in this case, the difference between the two is not zero:
 $$
[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:
 $$
{(x(t+\epsilon) x(t))^2 \over \epsilon} = f(t)
\,
and note that 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
 $\backslash mathcal\; L\; =\; (f(t)1)^2\; \backslash ,,$
and the equations of motion for f derived from extremizing the action S corresponding to $\backslash mathcal\; L$ just set it equal to 1. In physics, such a quantity is "equal to 1 as an identity operator". 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:
 $$
[x, \dot x] = x {dx\over dt}  {dx \over dt} x = 1
\,
This is called the Itō lemma in stochastic calculus, and the (euclideanized) canonical commutation relations in physics.
For a general statistical action, a similar argument shows that
 $$
\left[x , {\partial S \over \partial \dot x} \right] = 1
\,
and in quantum mechanics, the extra imaginary unit in the action converts this to the canonical commutation relation,
 $$
[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 timesliced flatspace 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—
 $Z\; =\; \backslash int\; Dx\backslash ,\; e^$
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.
Measure theoretic factors
Sometimes (e.g. a particle moving in curved space) we also have measuretheoretic factors in the functional integral.
 $\backslash int\; \backslash mu[x]\; e^\{iS[x]\}\; \backslash mathcal\{D\}x$
This factor is needed to restore unitarity.
For instance, if
 $S=\backslash int\; \backslash left[\; \backslash frac\{m\}\{2\}g\_\{ij\}\backslash dot\{x\}^i\backslash dot\{x\}^j\; \; V(x)\; \backslash right]\; dt$,
then it means that each spatial slice is multiplied by the measure √g. This measure can't be expressed as a functional multiplying the $\backslash mathcal\{D\}x$ measure because they belong to entirely different classes.
Quantum field theory
Quantum field theory 

Feynman diagram 
History 








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
 $$
[\phi(x),\partial_t \phi(y) ] = {\rm i} \delta^3(xy) \,
for x and y two simultaneous 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 make 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, the Heisenberg equations of motion, and the canonical commutation relations and shows that they are compatible with relativity. It extends the Heisenberg type operator algebra to operator product rules which are new relations difficult to see in 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 selfevident, 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 Grassmann 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 noncommuting. 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 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.
 $$
K(x,y;T) = \langle y;Tx;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 $x$ with an arbitrary initial state $\backslash psi\_0(x)$ constructs the final state.
 $$
\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 only a function of (x − y), the integral is a convolution, the final state is the initial state convolved with the propagator.
 $$
\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:
 $$
K(x,y;T) \propto e^{i m(xy)^2\over 2T}
Taking the Fourier transform in (x − y) produces another Gaussian:
 $$
K(p;T) = e^{i T p^2\over 2m}
and in pspace 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:
 $$
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 pspace 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 timereversible. 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 t is replaced by (−t). In this case, the interpretation is that these are the quantities to convolve the final wavefunction so as to get the initial wavefunction.
 $$
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.
 $$
K(xy,\Tau) = \int_{x(0)=x}^{x(\Tau)=y} e^{i \int_0^\Tau \sqrt + {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 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 p_{0} = m, the dominant first term has the form:
 $$
2m K_\mathrm{NR}(p) = {i \over (p_0m)  {\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 t → −∞ of the first term must vanish, while the t → +∞ limit of the second term must vanish. In the fourier transform, this means shifting the pole in p_{0} slightly, so that the inverse fourier transform will pick up a small decay factor in one of the time directions:
 $$
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 p_{0}. The terms can be recombined:
 $$
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 ε term introduces a small imaginary part to the α = m^{2}, which in the Minkowski version is a small exponential suppression of long paths.
So in the relativistic case, the Feynman pathintegral 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 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: S[ϕ] where the field ϕ(x^{μ}) 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 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 $\backslash 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
 $\backslash left\backslash langle\; F\backslash right\backslash rangle=\backslash frac\{\backslash int\; \backslash mathcal\{D\}\backslash phi\; F[\backslash phi]e^\{i\backslash mathcal\{S\}[\backslash phi]\}\}\{\backslash int\backslash mathcal\{D\}\backslash phi\; e^\{i\backslash mathcal\{S\}[\backslash phi]\}\}$
The symbol $\backslash int\; \backslash mathcal\{D\}\backslash phi$ here is a concise way to represent the infinitedimensional 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.
As a probability
Strictly speaking the only question that can be asked in physics is: "What fraction of states satisfying condition A also satisfy condition B?" The answer to this is a number between 0 and 1 which can be interpreted as a probability which is written as P(BA). In terms of path integration, since $P(BA)\; =\; \backslash frac\{P(A\; \backslash cap\; B)\}\{P(A)\}$ this means:
 $P(BA)\; =\; \backslash frac\{\backslash sum\_\{F\backslash subset\; A\; \backslash cap\; B\}\backslash left\; \backslash int\; \backslash mathcal\{D\}\backslash phi\; O\_\{in\}[\backslash phi]e^\{i\backslash mathcal\{S\}[\backslash phi]\}\; F[\backslash phi]\backslash right^2\}\{\backslash sum\_\{F\backslash subset\; A\}\; \backslash left\backslash int\backslash mathcal\{D\}\backslash phi\; O\_\{in\}[\backslash phi]\; e^\{i\backslash mathcal\{S\}[\backslash phi]\}\; F[\backslash phi]\backslash right^2\}$
where the functional O_{in}[ϕ] is the superposition of all incoming states that could lead to the states we are interested in. In particular this could be a state corresponding to the state of the Universe just after the big bang although for actual calculation this can be simplified using heuristic methods. Since this expression is a quotient of path integrals it is naturally normalised.
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
 $\backslash frac\{\backslash delta\; \backslash mathcal\{S\}[\backslash phi]\}\{\backslash delta\; \backslash phi\}=0$
(the lefthand 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 $\backslash mathcal\{D\}\backslash 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
 $e^\{i\backslash mathcal\{S\}[\backslash phi]\},$
which now becomes
 $e^\{H[\backslash phi]\}\backslash ,$
for some H, goes to zero faster than a 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 for the expectation:
 $\backslash left\backslash langle\; \backslash frac\{\backslash delta\; F[\backslash phi]\}\{\backslash delta\; \backslash phi\}\; \backslash right\backslash rangle\; =\; i\; \backslash left\backslash langle\; F[\backslash phi]\backslash frac\{\backslash delta\; \backslash mathcal\{S\}[\backslash phi]\}\{\backslash delta\backslash phi\}\; \backslash right\backslash rangle$
for any polynomially bounded functional F.
 $\backslash left\backslash langle\; F\_\{,i\}\; \backslash right\backslash rangle\; =\; i\; \backslash left\backslash langle\; F\; \backslash mathcal\{S\}\_\{,i\}\; \backslash right\backslash 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:
 $Z[J]=\backslash int\; \backslash mathcal\{D\}\backslash phi\; e^\{i(\backslash mathcal\{S\}[\backslash phi]\; +\; \backslash left\backslash langle\; J,\backslash phi\; \backslash right\backslash rangle)\}.$
Note that
 $\backslash frac\{\backslash delta^n\; Z\}\{\backslash delta\; J(x\_1)\; \backslash cdots\; \backslash delta\; J(x\_n)\}[J]\; =\; i^n\; \backslash ,\; Z[J]\; \backslash ,\; \{\backslash left\backslash langle\; \backslash phi(x\_1)\backslash cdots\; \backslash phi(x\_n)\backslash right\backslash rangle\}\_J$
or
 $Z^\{,i\_1\backslash dots\; i\_n\}[J]=i^n\; Z[J]\; \{\backslash left\; \backslash langle\; \backslash phi^\{i\_1\}\backslash cdots\; \backslash phi^\{i\_n\}\backslash right\backslash rangle\}\_J$
where
 $\{\backslash left\backslash langle\; F\; \backslash right\backslash rangle\}\_J=\backslash frac\{\backslash int\; \backslash mathcal\{D\}\backslash phi\; F[\backslash phi]e^\{i(\backslash mathcal\{S\}[\backslash phi]\; +\; \backslash left\backslash langle\; J,\backslash phi\; \backslash right\backslash rangle)\}\}\{\backslash int\backslash mathcal\{D\}\backslash phi\; e^\{i(\backslash mathcal\{S\}[\backslash phi]\; +\; \backslash left\backslash langle\; J,\backslash phi\; \backslash right\backslash rangle)\}\}.$
Basically, if $\backslash mathcal\{D\}\backslash phi\; e^\{i\backslash mathcal\{S\}[\backslash 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 $\backslash left\backslash langle\backslash phi(x\_1)\backslash cdots\; \backslash phi(x\_n)\backslash right\backslash 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
 $F[\backslash phi]=\backslash frac\{\backslash partial^\{k\_1\}\}\{\backslash partial\; x\_1^\{k\_1\}\}\backslash phi(x\_1)\backslash cdots\; \backslash frac\{\backslash partial^\{k\_n\}\}\{\backslash partial\; x\_n^\{k\_n\}\}\backslash phi(x\_n)$
and G is a functional of J, then
 $F\backslash left[i\backslash frac\{\backslash delta\}\{\backslash delta\; J\}\backslash right]\; G[J]\; =\; (i)^n\; \backslash frac\{\backslash partial^\{k\_1\}\}\{\backslash partial\; x\_1^\{k\_1\}\}\backslash frac\{\backslash delta\}\{\backslash delta\; J(x\_1)\}\; \backslash cdots\; \backslash frac\{\backslash partial^\{k\_n\}\}\{\backslash partial\; x\_n^\{k\_n\}\}\backslash frac\{\backslash delta\}\{\backslash delta\; J(x\_n)\}\; G[J].$
Then, from the properties of the functional integrals
 $\{\backslash left\; \backslash langle\; \backslash frac\{\backslash delta\; \backslash mathcal\{S\}\}\{\backslash delta\; \backslash phi(x)\}\backslash left[\backslash phi\; \backslash right]+J(x)\backslash right\backslash rangle\}\_J=0$
we get the "master" Schwinger–Dyson equation:
 $\backslash frac\{\backslash delta\; \backslash mathcal\{S\}\}\{\backslash delta\; \backslash phi(x)\}\backslash left[i\; \backslash frac\{\backslash delta\}\{\backslash delta\; J\}\backslash right]Z[J]+J(x)Z[J]=0$
or
 $\backslash mathcal\{S\}\_\{,i\}[i\backslash partial]Z+J\_i\; Z=0.$
If the functional measure is not translationally invariant, it might be possible to express it as the product $M\backslash left[\backslash phi\backslash right]\backslash ,\backslash mathcal\{D\}\backslash phi$ where M is a functional and $\backslash mathcal\{D\}\backslash phi$ is a translationally invariant measure. This is true, for example, for nonlinear sigma models where the target space is diffeomorphic to R^{n}. 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 $\backslash mathcal\{S\}$ in this equation by another functional $\backslash hat\{\backslash mathcal\{S\}\}=\backslash mathcal\{S\}i\backslash ln(M)$
If we expand this equation as a Taylor series about J = 0, we get the entire set of Schwinger–Dyson equations.
Localization
The path integrals are usually thought of as being the sum of all paths through an infinite space–time. However, in Local quantum field theory we would restrict everything to lie within a finite causally complete region, for example inside a double lightcone. This gives a more mathematically precise and physically rigorous definition of quantum field theory.
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
$\backslash int\; D[x]e^\{\backslash mathcal\{S\}[x]/\backslash hbar\}=A[x]\backslash sum\_\{n=0\}^\{\backslash infty\}(\backslash hbar)^\{n+1\}\backslash delta^\{n\}\; e^\{J/\backslash hbar\}\; \backslash text\{,\}$
where S is the Wickrotated classical action of the particle, J is the classical action with an extra term "x", delta (here) is the functional derivative operator and
 $A[x]=\backslash exp\backslash left(\backslash frac\{1\}\{\backslash hbar\}\backslash int\; X(t)\backslash ,\backslash mathrm\{d\}t\backslash right)\; \backslash text\{.\}$
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 $Q[\backslash mathcal\{L\}(x)]=\backslash partial\_\backslash mu\; f^\backslash 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 $\backslash int\; \backslash mathcal\{D\}\backslash phi\; Q[F][\backslash 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,
 $\backslash int\; \backslash mathcal\{D\}\backslash phi\backslash ,\; Q\backslash left[F\; e^\{iS\}\backslash right][\backslash phi]=0,$
which implies
 $\backslash left\backslash langle\; Q[F]\backslash right\backslash rangle\; +i\backslash left\backslash langle\; F\backslash int\_\{\backslash partial\; V\}\; f^\backslash mu\; ds\_\backslash mu\backslash right\backslash 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
 $Q=\backslash int\; d^dx\; q(x)$
where
 $q(x)[\backslash phi(y)]\; =\; \backslash delta^\{(d)\}(Xy)Q[\backslash phi(y)]\; \backslash ,$
so that
 $q(x)[S]=\backslash partial\_\backslash mu\; j^\backslash mu\; (x)\; \backslash ,$
where
 $j^\{\backslash mu\}(x)=f^\backslash mu(x)\backslash frac\{\backslash partial\}\{\backslash partial\; (\backslash partial\_\backslash mu\; \backslash phi)\}\backslash mathcal\{L\}(x)\; Q[\backslash phi]\; \backslash ,$
(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:
 $\backslash int\; \backslash mathcal\{D\}\backslash phi\backslash ,\; q(x)[F][\backslash phi]=0.$
Then, we would have
 $\backslash left\backslash langle\; q(x)[F]\; \backslash right\backslash rangle\; +i\backslash left\backslash langle\; F\; q(x)[S]\backslash right\backslash rangle=\backslash left\backslash langle\; q(x)[F]\backslash right\backslash rangle\; +i\backslash left\backslash langle\; F\backslash partial\_\backslash mu\; j^\backslash mu(x)\backslash right\backslash rangle=0.$
Alternatively,
 $q(x)[S]\backslash left[i\; \backslash frac\{\backslash delta\}\{\backslash delta\; J\}\backslash right]Z[J]+J(x)Q[\backslash phi(x)]\backslash left[i\; \backslash frac\{\backslash delta\}\{\backslash delta\; J\}\backslash right]Z[J]=\backslash partial\_\backslash mu\; j^\backslash mu(x)\backslash left[i\; \backslash frac\{\backslash delta\}\{\backslash delta\; J\}\backslash right]Z[J]+J(x)Q[\backslash phi(x)]\backslash left[i\; \backslash frac\{\backslash delta\}\{\backslash delta\; J\}\backslash right]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
 $\backslash left\backslash langle\; Q[F]\backslash right\backslash rangle\; =0.$
Alternatively,
 $\backslash int\; d^dx\backslash ,\; J(x)Q[\backslash phi(x)]\backslash left[i\; \backslash frac\{\backslash delta\}\{\backslash delta\; J\}\backslash right]Z[J]=0.$
The need for regulators and renormalization
Path integrals as they are defined here require the introduction of regulators. Changing the scale of the regulator leads to the renormalization group. In fact, renormalization is the major obstruction to making path integrals welldefined.
The path integral in quantummechanical 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. (Note that the Copenhagen/pragmatism interpretation claims there is no paradox—only a sloppy materialism motivated question on the part of EPR—Joseph Wienberg a lecture. On the other hand, the fact that the EPR thought experiment (and its result) does represent the results of a QM experiment says that (despite the path dependence of parallelness/antiparallelness in curved space) all contributions of paths close to black holes cancel in the action for an EPR style experiment here on earth.)
Some advocates of interpretations of quantum mechanics emphasizing decoherence have attempted to make more rigorous the notion of extracting a classicallike "coarsegrained" history from the space of all possible histories.
Quantum Gravity
Whereas in Quantum Theory the path integral formulation is fully equivalent to other formulations, it may be that it can be extended to quantum gravity, which would make it different from the Hilbert space model. Feynman had some success in this direction and his work has been extended by Hawking and others.^{[9]} The most developed approach is that of spinfoams.
See also
References
Notes
Suggested reading
 The historical reference, written by the inventor of the path integral formulation himself and one of his students.


 A highly readable introduction to the subject.
 A modern reference on the subject.


 Highly readable textbook; introduction to relativistic QFT for particle physics.




 This course, designed for mathematicians, is a rigorous introduction to perturbative quantum field theory, using the language of functional integrals.
 A great introduction to Path Integrals (Chapter 1) and QFT in general.
 Grosche, Christian (1992). "An Introduction into the Feynman Path Integral". hepth/9302097.
 MacKenzie, Richard (2000). "Path Integral Methods and Applications". quantph/0004090.



External links
 Path integral on Scholarpedia
 Path Integrals in Quantum Theories: A Pedagogic 1st Step
This article was sourced from Creative Commons AttributionShareAlike License; additional terms may apply. World Heritage Encyclopedia content is assembled from numerous content providers, Open Access Publishing, and in compliance with The Fair Access to Science and Technology Research Act (FASTR), Wikimedia Foundation, Inc., Public Library of Science, The Encyclopedia of Life, Open Book Publishers (OBP), PubMed, U.S. National Library of Medicine, National Center for Biotechnology Information, U.S. National Library of Medicine, National Institutes of Health (NIH), U.S. Department of Health & Human Services, and USA.gov, which sources content from all federal, state, local, tribal, and territorial government publication portals (.gov, .mil, .edu). Funding for USA.gov and content contributors is made possible from the U.S. Congress, EGovernment Act of 2002.
Crowd sourced content that is contributed to World Heritage Encyclopedia is peer reviewed and edited by our editorial staff to ensure quality scholarly research articles.
By using this site, you agree to the Terms of Use and Privacy Policy. World Heritage Encyclopedia™ is a registered trademark of the World Public Library Association, a nonprofit organization.