# Path Integrals in Quantum Mechanics

*Michael Fowler, UVa*

### Huygen’s Picture of Wave Propagation

If a point source of light is switched on, the wavefront is an expanding sphere centered at the source. Huygens suggested that this could be understood if at any instant in time each point on the wavefront was regarded as a source of secondary wavelets, and the new wavefront a moment later was to be regarded as built up from the sum of these wavelets. For a light shining continuously, this process just keeps repeating.

What use is this idea? For one thing, it explains refraction—the change in direction of a wavefront on entering a different medium, such as a ray of light going from air into glass.

If the light moves more slowly in the glass, velocity $v$ instead of $c,$ with $vc,$ then Huygen’s picture explains Snell’s Law,
that the ratio of the sines of the angles to the normal of incident and
transmitted beams is constant, and in fact is the ratio $c/v.$ This is evident from the diagram below: in
the time the wavelet centered at ** A** has propagated to

**, that from**

*C***has reached**

*B***, the ratio of lengths**

*D***/**

*AC***being $c/v.$ But the angles in Snell’s Law are in fact the angles**

*BD***,**

*ABC***, and those right-angled triangles have a common hypotenuse**

*BCD***, from which the Law follows.**

*BC*### Fermat’s Principle of Least Time

We will now temporarily forget about the wave nature of
light, and consider a narrow ray or beam of light shining from point *A *to
point *B*, where we suppose *A* to be in air, *B* in glass. Fermat showed that the path of such a beam is
given by the Principle of Least Time: a ray of light going from *A* to *B*
by any other path would take longer. How can we see that? It’s obvious that any
deviation from a straight line path in air or in the glass is going to add to
the time taken, but what about moving slightly the point at which the beam
enters the glass?

Where the air meets the glass, the two rays, separated by a
small distance ** CD** =

*d*along that interface, will look parallel:

(Feynman gives a nice illustration: a lifeguard on a beach spots a swimmer in trouble some distance away, in a diagonal direction. He can run three times faster than he can swim. What is the quickest path to the swimmer?)

Moving the point of entry up a small distance $d,$ the light has to travel an extra $d\mathrm{sin}{\theta}_{1}$ in air, but a distance less by $d\mathrm{sin}{\theta}_{2}$ in the glass, giving an extra travel time $\Delta t=d\mathrm{sin}{\theta}_{1}/c-d\mathrm{sin}{\theta}_{2}/v$.
For the *classical* path, Snell’s
Law gives $\mathrm{sin}{\theta}_{1}/\mathrm{sin}{\theta}_{2}=n=c/v$,
so $\Delta t=0$ to first order. But if we look at a series of
possible paths, each a small distance $d$ away from the next at the point of crossing
from air into glass, $\Delta t$ becomes of order $d/c$ away from the classical path.

Suppose now we imagine that the light actually travels along
all these paths with about equal amplitude. What
will be the total contribution of all the paths at *B*? Since the times along the paths are
different, the signals along the different paths will arrive at *B* with
different phases, and to get the total wave amplitude we must add a series of
unit 2*D* vectors, one from each
path. (Representing the amplitude and
phase of the wave by a complex number for convenience$\u2014$for a real wave,
we can take the real part at the end.)

When we map out these unit 2*D* vectors, we find that in the neighborhood of the classical path,
the phase varies little, but as we go away from it the phase spirals more and
more rapidly, so those paths interfere amongst themselves destructively. To formulate this a little more precisely,
let us assume that some close by path has a phase difference $\phi $ from the least time path, and goes from air to
glass a distance $x$ away from the least time path: then for these
close by paths, $\phi =a{x}^{2},$ where $a$ depends on the geometric arrangement and the
wavelength. From this, the sum over the
close by paths is an integral of the form $\int {e}^{ia{x}^{2}}dx$. (We are assuming the wavelength of light is
far less than the size of the equipment.)
This is a standard integral, its value is $\sqrt{\pi /ia},$ all its weight is concentrated in a central
area of width $1/\sqrt{a},$ exactly as for the real function ${e}^{-a{x}^{2}}.$

This is the explanation of Fermat’s Principle$\u2014$only near the
path of least time do paths stay approximately in phase with each other and add
constructively. So this classical path rule has an underlying wave-phase
explanation. In fact, the central role
of phase in this analysis is sometimes emphasized by saying the light beam
follows *the path of stationary phase*.

Of course, we’re not summing over *all* paths here$\u2014$we assume that
the path in air from the source to the point of entry into the glass is a
straight line, clearly the subpath of stationary phase.

### Classical Mechanics: The Principle of Least Action

Confining our attention for the moment to the mechanics of a single nonrelativistic particle in a potential, with Lagrangian $L=T-V,$ the action $S$ is defined by

$S={\displaystyle \underset{{t}_{1}}{\overset{{t}_{2}}{\int}}L(x,\dot{x})dt}.$

*A* at ${t}_{1}$ to *B* at ${t}_{2}$ travels along the path that minimizes the
action. This is called the *Principle
of Least Action*: for example, the parabolic path followed by a ball thrown
through the air minimizes the integral along the path of the action $T-V$ where $T$ is the ball’s kinetic energy, $V$ its gravitational potential energy (neglecting
air resistance, of course). Note here
that the initial and final times are fixed, so since we’ll be summing over
paths with different lengths, necessarily the particles speed will be different
along the different paths. In other words, it will have different energies
along the different paths.

With the advent of quantum mechanics, and the realization
that any particle, including a thrown ball, has wave like properties, the
rather mysterious Principle of Least Action looks a lot like Fermat’s Principle
of Least Time. Recall that Fermat’s Principle
works because the total phase along a path is the integrated time elapsed along
the path, and for a path where that integral is stationary for small path
variations, neighboring paths add constructively, and no other sets of paths
do. If the Principle of Least Action has
a similar explanation, then the wave amplitude for a particle going along a
path from *A* to *B* must have a phase equal to some constant times
the action along that path. If this is the case, then the observed path
followed will be just that of least action, or, more generally, of *stationary* action, for only near that
path will the amplitudes add constructively, just as in Fermat’s analysis of
light rays.

### Going from Classical Mechanics to Quantum Mechanics

Of course, if we write a phase factor for a path ${e}^{icS}$ where $S$ is the action for the path and $c$ is some constant, $c$ must necessarily have the dimensions of inverse action. Fortunately, there is a natural candidate for the constant $c.$ The wave nature of matter arises from quantum mechanics, and the fundamental constant of quantum mechanics, Planck’s constant, is in fact a unit of action. (Recall action has the same dimensions as $Et,$ and therefore the same as $px,$ manifestly the same as angular momentum.) It turns out that the appropriate path phase factor is ${e}^{iS/\hslash}$

That the phase factor is ${e}^{iS/\hslash}$, rather than ${e}^{iS/h}$, say, can be established by considering the double slit experiment for electrons (Peskin page 277).

This is analogous to the light waves going from a source in air to a point in glass, except now we have vacuum throughout (electrons don’t get far in glass), and we close down all but two of the paths.

Suppose electrons from the top slit, Path I, go a distance $D$ to the detector, those from the bottom slit, Path II, go $D+d,$ with $d\ll D.$ Then if the electrons have wavelength $\lambda $ we know the phase difference at the detector is $2\pi d/\lambda .$ To see this from our formula for summing over paths, on Path I the action $S=Et={\scriptscriptstyle \frac{1}{2}}m{v}_{1}^{2}t,$ and ${v}_{1}=D/t,$ so

${S}_{1}={\scriptscriptstyle \frac{1}{2}}m{D}^{2}/t.$

For Path II, we must take ${v}_{2}=\left(D+d\right)/t.$ Keeping only terms of leading order in $d/D,$ the action difference between the two paths

$${S}_{2}-{S}_{1}=mDd/t$$

so the phase difference

$$\frac{{S}_{2}-{S}_{1}}{\hslash}=\frac{mvd}{\hslash}=\frac{2\pi pd}{h}=\frac{2\pi d}{\lambda}.$$

This is the known correct result, and this fixes the
constant multiplying the action/*h* in the expression for the path phase.

In quantum mechanics, such as the motion of an electron in
an atom, we know that the particle does not follow a well-defined path, in
contrast to classical mechanics. Where
does the crossover to a well-defined path take place? Taking the simplest possible case of a free
particle (no potential) of mass $m$ moving at speed $v,$ the action along a straight line path taking
time $t$ from *A* to *B* is ${\scriptscriptstyle \frac{1}{2}}m{v}^{2}t.$ If this
action is of order Planck’s constant $h,$ then the phase factor will not oscillate
violently on moving to different paths, and a range of paths will
contribute. In other words, quantum
rather than classical behavior dominates when ${\scriptscriptstyle \frac{1}{2}}m{v}^{2}t$ is of order $h.$ But $vt$ is the path length $L,$ and $mv/h$ is the wavelength $\lambda ,$ so we conclude that we must use quantum
mechanics when the wavelength $h/p$ is significant compared with the path length. Interference sets in when the *difference*
in path actions is of order $h,$ so in the atomic regime many paths must be
included.

Feynman (in Feynman and Hibbs) gives a nice picture to help
think about summing over paths. He begins with the double slit experiment for
an electron. We suppose the electron is emitted
from some source *A* on the left, and we look for it at a point *B*
on a screen to the right. In the middle
is a thin opaque barrier with the familiar two slits. Evidently, to find the amplitude for the
electron to reach *B *we sum over two paths. Now suppose we add *another* two-slit
barrier. We have to sum over four paths. Now add another. Next, replace the two slits in each barrier by
several slits. We must sum over a
multitude of paths! Finally, increase
the number of barriers to some large number *N*, and at the same time
increase the number of slits to the point that there are no barriers left. We are left with a sum over all possible paths
through space from *A* to *B*, multiplying each path by the
appropriate action phase factor. This is
reminiscent of the original wave propagation picture of Huygens: if one
pictures it at successive time intervals of picoseconds, say, from each point
on the wavefront waves go out 3 mm in all directions, then in the next time
interval each of those sprouts more waves in all directions. One could write this as a sum over all zigzag
paths with random 3 mm steps.

In fact, the sum over paths is even more daunting than
Feynman’s picture suggests. All the
paths going through these many slitted barriers are progressing in a forward
direction, from *A* towards *B*.
Actually, if we’re summing over *all*
paths, we should be including the possibility of paths zigzagging backwards and
forwards as well, eventually arriving at *B*. We shall soon see how to deal systematically
with all possible paths.

### Review: Standard Definition of the Free Electron Propagator

As a warm up exercise, consider an electron confined to one dimension, with no potential present, moving from ${x}^{\prime}$ at time 0 to $x$ at time $T.$ We’ll follow Feynman in using $T$ for the final time, so we can keep $t$ for the continuous (albeit sometimes discretized) time variable over the interval 0 to $T.$

(As explained previously, when we write that the electron is initially at ${x}^{\prime},$ we mean its wave function is a normalizable state, such as a very narrow Gaussian, centered at ${x}^{\prime}.$ The propagator then represents the probability amplitude, that is, the wave function, at point $x$ after the given time $T.$ ) The propagator is given by

$|\psi \left(x,t=T\right)\rangle =U\left(T\right)|\psi \left(x,t=0\right)\rangle ,$

or, in Schrödinger wave function notation,

$\psi \left(x,T\right)={\displaystyle \int U\left(x,T;{x}^{\prime},0\right)\psi \left({x}^{\prime},0\right)\text{\hspace{0.17em}}}d{x}^{\prime}.$

It is clear that for this to make sense, as $T\to 0,\text{\hspace{0.17em}}\text{\hspace{1em}}U\left(x,T;{x}^{\prime},0\right)\to \delta \left(x-{x}^{\prime}\right).$

In the lecture on propagators, we found

$\begin{array}{c}\langle x|U\left(T,0\right)|{x}^{\prime}\rangle \text{\hspace{0.17em}}={\displaystyle \underset{-\infty}{\overset{\infty}{\int}}{e}^{-i\hslash {k}^{2}T/2m}\frac{dk}{2\pi}\langle x|k\rangle \langle k|{x}^{\prime}\rangle}\\ ={\displaystyle \underset{-\infty}{\overset{\infty}{\int}}{e}^{-i\hslash {k}^{2}T/2m}\frac{dk}{2\pi}{e}^{-ik\left(x-{x}^{\prime}\right)}}\\ =\sqrt{\frac{m}{2\pi \hslash iT}}{e}^{im{\left(x-{x}^{\prime}\right)}^{2}/2\hslash T}.\end{array}$

### Summing over Paths

Let us formulate the sum over paths for this simplest one-dimensional case, the free electron, more precisely. Each path is a continuous function of time $x\left(t\right)$ in the time interval $0\le t\le T,$ with boundary conditions $x\left(0\right)={x}^{\prime},\text{\hspace{0.17em}}x\left(T\right)=x.$ Each path contributes a term ${e}^{iS/\hslash},$ where

$S\left[x\left(t\right)\right]={\displaystyle \underset{0}{\overset{T}{\int}}L\left(x\left(t\right),\dot{x}\left(t\right)\right)dt}={\displaystyle \underset{0}{\overset{T}{\int}}{\scriptscriptstyle \frac{1}{2}}m{\dot{x}}^{2}\left(t\right)dt}$

(for the free electron case) evaluated along that path.

The integral over all paths is written:

$\langle x|U\left(T,0\right)|{x}^{\prime}\rangle ={\displaystyle \int D\left[x\left(t\right)\right]}\text{\hspace{0.17em}}{e}^{iS\left[x\left(t\right)\right]/\hslash}$

This rather formal statement begs the question of how, exactly, we perform the sum over paths: what is the appropriate measure in the space of paths?

A natural approach is to measure the paths in terms of their deviation from the classical path, since we know that path dominates in the classical limit. The classical path for the free electron is just the straight line from ${x}^{\prime}$ to $x,$ traversed at constant velocity, since there are no forces acting on the electron.

We write

$x\left(t\right)={x}_{cl}\left(t\right)+y\left(t\right)$

where

${x}_{cl}\left(0\right)={x}^{\prime},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{x}_{cl}\left(T\right)=x$

and therefore

$y\left(0\right)=0,\text{\hspace{1em}}y\left(T\right)=0.$

Then

$\begin{array}{l}\langle x|U\left(T,0\right)|{x}^{\prime}\rangle ={\displaystyle \int D\left[y\left(t\right)\right]}\text{\hspace{0.17em}}{e}^{iS\left[{x}_{cl}\left(t\right)+y\left(t\right)\right]/\hslash},\text{\hspace{1em}}\\ S\left[{x}_{cl}\left(t\right)+y\left(t\right)\right]={\displaystyle \underset{0}{\overset{T}{\int}}{\scriptscriptstyle \frac{1}{2}}m{\left({\dot{x}}_{cl}\left(t\right)+\dot{y}\left(t\right)\right)}^{2}dt}\\ =S\left[{x}_{cl}\left(t\right)\right]+{\displaystyle \underset{0}{\overset{T}{\int}}m{\dot{x}}_{cl}\left(t\right)\dot{y}\left(t\right)dt}+{\displaystyle \underset{0}{\overset{T}{\int}}{\scriptscriptstyle \frac{1}{2}}m{\dot{y}}^{2}\left(t\right)dt}.\end{array}$

The middle term on the bottom line is zero, as it has to be since it is a linear term in the deviation from the minimum path. To see this explicitly, one can integrate by parts: the end terms are zero, from the boundary condition on $y,$ and the other term is the acceleration of the particle along the classical path, which is zero.

Therefore

$$\langle x|U\left(T,0\right)|{x}^{\prime}\rangle ={e}^{iS\left[{x}_{cl}\left(t\right)\right]/\hslash}{\displaystyle \int D\left[y\left(t\right)\right]}\text{\hspace{0.17em}}{e}^{iS\left[y\left(t\right)\right]/\hslash}$$

The $y\text{-}$ paths, being the deviation from the classical
path from ${x}^{\prime}$ to $x,$ necessarily begin and end at the $y\text{-}$ origin, since *all* paths summed over go from ${x}^{\prime}$ to $x.$

The classical path, motion from ${x}^{\prime}$ to $x$ at a constant speed $v=\left({x}^{\prime}-x\right)/T,$ has action $Et,$ with $E$ the classical energy ${\scriptscriptstyle \frac{1}{2}}m{v}^{2},$ so

$U(x,T;{x}^{\prime},0)=A\left(T\right){e}^{im{(x-{x}^{\prime})}^{2}/2\hslash T}.$

This gives the correct exponential term. The prefactor $A,$ representing the sum over the deviation paths $y\left(t\right),$ cannot depend on $x$ or ${x}^{\prime},$ and is fixed by the requirement that as $t$ goes to zero, $U$ must approach a $\delta \text{-}$ function, giving the prefactor found previously.

### Proving that the Sum-Over-Paths Definition of the Propagator is Equivalent to the Sum-Over-Eigenfunctions Definition

The first step is to construct a practical method of summing
over paths. Let us begin with a particle
in one dimension going from ${x}^{\prime}$ at time 0 to $x$ at time $T.$ The
paths can be enumerated in a crude way, reminiscent of Riemann integration: divide
the time interval 0 to $T$ into *N* equal intervals each of duration
$\epsilon ,$ so ${t}_{0}=0,\text{\hspace{0.17em}}{t}_{1}={t}_{0}+\epsilon ,\text{\hspace{0.17em}}{t}_{2}={t}_{0}+2\epsilon ,\dots ,\text{\hspace{0.17em}}{t}_{N}=T.$

Next, define a particular path from $x$ to ${x}^{\prime}$ by specifying the position of the particle at each of the intermediate times, that is to say, it is at ${x}_{1}$ at time ${t}_{1},\text{\hspace{0.17em}}{x}_{2}$ at time ${t}_{2}$ and so on. Then, simplify the path by putting in straight line bits connecting ${x}_{0}$ to ${x}_{1},\text{\hspace{0.17em}}{x}_{1}$ to ${x}_{2},$ etc. The justification is that in the limit of $\epsilon $ going to zero, taken at the end, this becomes a true representation of the path.

The next step is to sum over all possible paths with a
factor ${e}^{iS/\hslash}$ for each one. The sum is accomplished by integrating over
all possible values of the intermediate positions ${x}_{1},{x}_{2},\dots ,{x}_{N-1},$ and
then taking *N* to infinity.

The action on the zigzag path is

$S={\displaystyle \underset{0}{\overset{T}{\int}}dt({\scriptscriptstyle \frac{1}{2}}m{\dot{x}}^{2}-V(x))\to {\displaystyle \sum _{i}\left[\frac{m{({x}_{i+1}-{x}_{i})}^{2}}{2\epsilon}-\epsilon \text{\hspace{0.17em}}V\left(\frac{{x}_{i+1}+{x}_{i}}{2}\right)\right]}}$

We define the “integral over paths” written $\int D[x(t)]$ by

$\underset{N\to \infty}{\underset{\epsilon \to 0}{\mathrm{lim}}}\frac{1}{B(\epsilon )}{\displaystyle \underset{-\infty}{\overset{\infty}{\int}}{\displaystyle \int \mathrm{...}{\displaystyle \int \frac{d{x}_{1}}{B(\epsilon )}\text{\hspace{0.17em}}}}}\mathrm{...}\text{\hspace{0.17em}}\frac{d{x}_{N-1}}{B(\epsilon )}$

where we haven’t yet figured out what the overall weighting factor $B(\epsilon )$ is going to be. (It is standard convention to have that extra $B(\epsilon )$ outside.)

To summarize: the propagator $U(x,T;{x}^{\prime},0)$ is the contribution to the wave function at $x$ at time

$t=T$ from that at ${x}^{\prime}$ at the earlier time $t=0.$

Consequently, $U(x,T;{x}^{\prime},0)$ regarded as a function of $x,T$ is, in fact, nothing but the Schrödinger wave
function $\psi \left(x,T\right),$ and therefore must *satisfy Schrödinger’s
equation*

$i\hslash \frac{\partial}{\partial T}U\left(x,T;{x}^{\prime},0\right)=\left(-\frac{{\hslash}^{2}}{2m}\frac{{\partial}^{2}}{\partial {x}^{2}}+V\left(x\right)\right)U\left(x,T;{x}^{\prime},0\right).$

We shall now show that defining $U\left(x,T;{x}^{\prime},0\right)$ as a sum over paths, it does in fact satisfy Schrödinger’s equation, and furthermore goes to a $\delta \text{-}$ function as time goes to zero.

$\begin{array}{c}U\left(x,T;{x}^{\prime},0\right)={\displaystyle \int D[x(t)]{e}^{iS[x(t)]/\hslash}}\\ =\underset{N\to \infty}{\underset{\epsilon \to 0}{\mathrm{lim}}}\frac{1}{B(\epsilon )}{\displaystyle \underset{-\infty}{\overset{\infty}{\int}}{\displaystyle \int \mathrm{...}{\displaystyle \int \frac{d{x}_{1}}{B(\epsilon )}\text{\hspace{0.17em}}}}}\mathrm{...}\text{\hspace{0.17em}}\frac{d{x}_{N-1}}{B(\epsilon )}{e}^{iS({x}_{1},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{...}\text{\hspace{0.17em}}\text{\hspace{0.17em}},{x}_{N-1})/\hslash}.\end{array}$

We shall establish this equivalence by proving that it satisfies the same differential equation. It clearly has the same initial value$\u2014$as ${t}^{\prime}$ and $t$ coincide, it goes to $\delta \left(x-{x}^{\prime}\right)$ in both representations.

To differentiate $U\left(x,T;{x}^{\prime},0\right)$ with respect to $t,$ we isolate the integral over the last path variable, ${x}_{N-1}:$

$U\left(x,T;{x}^{\prime},0\right)={\displaystyle \int \frac{d{x}_{N-1}}{B(\epsilon )}{e}^{\left[\frac{im{(x-{x}_{N-1})}^{2}}{2\hslash \epsilon}-\frac{i}{\hslash}\epsilon V\left(\frac{x+{x}_{N-1}}{2}\right)\right]}}U\left({x}_{N-1},T-\epsilon ;{x}^{\prime},0\right)$

Now in the limit $\epsilon $ going to zero, almost all the contribution to
this integral must come from close to the point of stationary phase, that is, ${x}_{N-1}=x.$ In that limit, we can take $U({x}_{N-1},t-\epsilon ;{x}^{\prime},{t}^{\prime})$ to be a slowly varying function of ${x}_{N-1},$ and replace it by the leading terms in a *x*, so

$\begin{array}{l}U\left(x,T;{x}^{\prime},0\right)=\\ {\displaystyle \int \frac{d{x}_{N-1}}{B(\epsilon )}{e}^{\frac{im{(x-{x}_{N-1})}^{2}}{2\hslash \epsilon}}\left(1-\frac{i}{\hslash}\epsilon V\left(\frac{x+{x}_{N-1}}{2}\right)\right)}\left(U(x,T-\epsilon )+({x}_{N-1}-x)\frac{\partial U}{\partial x}+\frac{{({x}_{N-1}-x)}^{2}}{2}\frac{{\partial}^{2}U}{\partial {x}^{2}}\right)\end{array}$

The ${x}_{N-1}$ dependence in the potential $V$ can be neglected in leading order$\u2014$that leaves standard Gaussian integrals, and

$U\left(x,T;{x}^{\prime},0\right)=\frac{1}{B(\epsilon )}\sqrt{\frac{2\pi \hslash \epsilon}{-im}}\left(1-\frac{i\epsilon}{\hslash}V(x)+\frac{i\epsilon \hslash}{2m}\frac{{\partial}^{2}}{\partial {x}^{2}}\right)U\left(x,T-\epsilon ;{x}^{\prime},0\right).$

Taking the limit of $\epsilon $ going to zero fixes our unknown normalizing factor,

$B(\epsilon )=\sqrt{\frac{2\pi \hslash \epsilon}{-im}}$

giving

$i\hslash \frac{\partial}{\partial T}U\left(x,T;{x}^{\prime},0\right)=\left(-\frac{{\hslash}^{2}}{2m}\frac{{\partial}^{2}}{\partial {x}^{2}}+V(x)\right)U\left(x,T;{x}^{\prime},0\right)$,

thus establishing that the propagator derived from the sum over paths obeys Schrödinger’s equation, and consequently gives the same physics as the conventional approach.

### Explicit Evaluation of the Path Integral for the Free Particle Case

The required correspondence to the Schrödinger equation result fixes the unknown normalizing factor, as we’ve just established. This means we are now in a position to evaluate the sum over paths explicitly, at least in the free particle case, and confirm the somewhat hand-waving result given above.

The sum over paths is

$U\left(x,T;{x}^{\prime},0\right)={\displaystyle \int D[x(t)]{e}^{iS[x(t)]/\hslash}}=\underset{N\to \infty}{\underset{\epsilon \to 0}{\mathrm{lim}}}\frac{1}{B(\epsilon )}{\displaystyle \underset{-\infty}{\overset{\infty}{\int}}{\displaystyle \int \mathrm{...}{\displaystyle \int \frac{d{x}_{1}}{B(\epsilon )}\text{\hspace{0.17em}}}}}\mathrm{...}\text{\hspace{0.17em}}\frac{d{x}_{N-1}}{B(\epsilon )}{e}^{i{\displaystyle \sum _{i}\frac{m{({x}_{i+1}-{x}_{i})}^{2}}{2\hslash \epsilon}}}.$

Let us consider the sum for small but finite $\epsilon .$ In particular, we’ll divide up the interval
first into halves, then quarters, and so on, into 2* ^{n}* small intervals.
The reason for this choice will become clear.

Now, we’ll integrate over half the paths: those for $i$ odd, leaving the even ${x}_{i}$ values fixed for the moment. The integrals are of the form

$$\begin{array}{c}{\displaystyle \underset{-\infty}{\overset{\infty}{\int}}dy{e}^{\left(ia/2\right)\left[{\left(x-y\right)}^{2}+{\left(y-z\right)}^{2}\right]}}={e}^{\left(ia/2\right)\left({x}^{2}+{z}^{2}\right)}{\displaystyle \underset{-\infty}{\overset{\infty}{\int}}dy{e}^{ia{y}^{2}-iay\left(x+z\right)}}\\ ={e}^{\left(ia/2\right)\left({x}^{2}+{z}^{2}\right)}\sqrt{\frac{\pi}{-ia}}{e}^{\left(-ia/4\right){\left(x+z\right)}^{2}}=\sqrt{\frac{\pi}{-ia}}{e}^{\left(ia/4\right){\left(x-z\right)}^{2}}\end{array}$$

using the standard result $\underset{-\infty}{\overset{\infty}{\int}}dx{e}^{-a{x}^{2}+bx}=\sqrt{\frac{\pi}{a}}{e}^{{b}^{2}/4a}}.$

Now put in the value $a=m/\hslash \epsilon :$ the factor $\sqrt{\frac{\pi}{-ia}}=\sqrt{\frac{\pi \hslash \epsilon}{-im}}$ cancels the normalization factor $B(\epsilon )=\sqrt{\frac{2\pi \hslash \epsilon}{-im}}$ except for the factor of 2 inside the square root. But we need that factor of 2, because we’re left with an integral$\u2014$over the remaining even numbered paths$\u2014$exactly like the one before except that the time interval has doubled, both in the normalization factor and in the exponent, $\epsilon \to 2\epsilon .$

So we’re back where we started. We can now repeat the process, halving the number of paths again, then again, until finally we have the same expression but with only the fixed endpoints appearing.