232B Lecture Notes

Path Integral I
Hitoshi Murayama
1

Feynmans Path Integral Formulation

Feynmans formulation of quantum mechanics using the so-called path integral is arguably the most elegant. It can be stated in a single line:
hxf , tf |xi , ti i =

Dx(t)eiS[x(t)]/h .

(1)

The meaning of this equation is the following. If you want to know the
quantum mechanical amplitude for a point particle at a position xi at time
ti to reach a position xf at time tf , you integrate over all possible paths
connecting the points with a weight factor given by the classical action for
each path. Hence the name path integral. This is it. Note that the position
kets form a complete set of basis, and knowing this amplitude for all x is
enough information to tell you everything about the system. The expression
is generalized for more dimensions and more particles in a straightforward
manner.
As we will see, this formulation is completely equivalent to the usual
formulation of quantum mechanics. On the other hand, there are many
reasons why this expression is just beautiful.
First, the classical equation of motion comes out in a very simple way.
If you take the limit h
0, the weight factor eiS/h oscillates very rapidly.
Therefore, we expect that the main contribution to the path integral comes
from paths that make the action stationary. This is nothing but the derivation of EulerLagrange equation from the classical action. Therefore, the
classical trajectory dominates the path integral in the small h
limit.
Second, we dont know what path the particle has chosen, even when we
know what the initial and final positions are. This is a natural generalization
of the two-slit experiment. Even if we know where the particle originates from
and where it hit on the screen, we dont know which slit the particle came
through. The path integral is an infinite-slit experiment. Because you cant
specificy where the particle goes through, you sum them up.

Third, we gain intuition on what quantum fluctuation does. Around the

classical trajectory, a quantum particle explores the vicinity. The trajectory can deviate from the classical trajectory if the difference in the action is
roughly within h
. When a classical particle is confined in a potential well, a
quantum particle can go on an excursion and see that there is a world outside
the potential barrier. Then it can decide to tunnel through. If a classical
particle is sitting at the top of a hill, it doesnt fall; but a quantum particle
realizes that the potential energy can go down with a little excursion, and
decides to fall.
Fourth, whenever we have an integral expression for a quantity, it is often
easier to come up with an approximation method to work it out, compared
to staring at a differential equation. Good examples are the perturbative
expansion and the steepest descent method. One can also think of useful
change of variables to simplify the problem. In fact, some techniques in
quantum physics couldnt be thought of without the intuition from the path
integral.
Fifth, the path integral can also be used to calculate partition functions
in statistical mechanics.
Sixth, the connection between the conservation laws and unitarity transformations become much clearer with the path integral. We will talk about
it when we discuss symmetries.
There are many more but I stop here.
Unfortunately, there is also a downside with the path integral. The actual calculation of a path integral is somewhat technical and awkward. One
might even wonder if an integral over paths is mathematically well-defined.
Mathematicians figured that it can actually be, but nonetheless it makes some
people nervous. Reading off energy eigenvalues is also less transparent than
with the Schrodinger equation.
Below, we first derive the path integral from the conventional quantum
mechanics. Then we show that the path integral can derive the conventional
Schrodinger equation back. After that we look at some examples and actual
calculations.
By the way, the original paper by Feynman on the path integral Rev.
Mod. Phys. 20, 367-387 (1948) is quite readable for you, and I recommend
it. Another highly recommended read is Feynmans Nobel Lecture. You
will see that Feynman invented the path integral with the hope of replacing
quantum field theory with particle quantum mechanics; he failed. But the
path integral survived and did mighty good in the way he didnt imagine.
2

The derivation of path integral in this lecture note is based on J.J. Sakurai, Modern Quantum Mechanics.

Physical Intuition

Take the two-slit experiment. Each time an electron hits the screen, there is
no way to tell which slit the electron has gone through. After repeating the
same experiment many many times, a fringe pattern gradually appears on the
screen, proving that there is an interference between two waves, one from one
slit, the other from the other. We conclude that we need to sum amplitudes
of these two waves that correspond to different paths of the electron. Now
imagine that you make more slits. There are now more paths, each of which
contributing an amplitude. As you increase the number of slits, eventually
the entire obstruction disappears. Yet it is clear that there are many paths
that contribute to the final amplitude of the electron propagating to the
screen.
As we generalize this thought experiment further, we are led to conclude
that the amplitude of a particle moving from a point xi to another point
xf consists of many components each of which corresponds to a particular
path that connects these two points. One such path is a classical trajectory.
However, there are infinitely many other paths that are not possible classically, yet contribute to the quantum mechanical amplitude. This argument
leads to the notion of a path integral , where you sum over all possible paths
connecting the initial and final points to obtain the amplitude.
The question then is how you weight individual paths. One point is clear:
the weight factor must be chosen such that the classical path is singled out
in the limitR h
0. The correct choice turns out to be eiS[x(t)]/h , where
tf
S[x(t)] = ti dtL(x(t), x(t))

is the classical action for the path x(t) that

satisfies the boundary condition x(ti ) = xi , x(tf ) = xf . In the limit h
0,
the phase factor oscillates so rapidly that nearly all the paths would cancel
each other out in the final amplitude. However, there is a path that makes the
action stationary, whose contribution is not canceled. This particular path
is nothing but the classical trajectory. This way, we see that the classical
trajectory dominates the path integral in the h
0 limit.
As we increase h
, the path becomes fuzzy. The classical trajectory
still dominates, but there are other paths close to it whose action is within
S ' h
and contribute significantly to the amplitude. The particle does an
3

excursion around the classical trajectory.

Propagator

The quantity
K(xf , tf ; xi , ti ) = hxf , tf |xi , ti i

(2)

is called a propagator . It knows everything about how a wave function propagates in time, because
(xf , tf ) = hxf , tf |i =
=

hxf , tf |xi , ti idxi hxi , ti |i

K(xf , tf ; xi , ti )(xi , ti )dxi .

(3)

In other words, it is the Greens function for the Schrodinger equation.

The propagator can also be written using energy eigenvalues and eigenstates (if the Hamiltonian does not depend on time),
K(xf , tf ; xi , ti ) = hxf |eiH(tf ti )/h |xi i =
=

iEn (tf ti )/
h

hxf |nieiEn (tf ti )/h hn|xi i

n (xf )n (xi ).

(4)

In particular, Fourier analyzing the propagator tells you all energy eigenvalues, and each Fourier coefficients the wave functions of each energy eigenstates.
The propagator is a nice package that contains all dynamical information
about a quantum system.

Derivation of the Path Integral

The basic point is that the propagator for a short interval is given by the
classical Lagrangian
2 )/
h

hx1 , t + t|x0 , ti = cei(L(t)t+O(t)

(5)

Note that the expression here is in the Heisenberg picture. The base kets |xi and the
operator x depend on time, while the state ket |i doesnt.

where c is a normalization constant. This can be shown easily for a simple

Hamiltonian
p2
H=
+ V (x).
(6)
2m
The quantity we want is
iHt/
h

hx1 , t + t|x0 , ti = hx1 |e

|x0 i =

dphx1 |pihp|eiHt/h |x0 i.

(7)

Because we are interested in the phase factor only at O(t), the last factor
can be estimated as
hp|eiHt/h |x0 i
= hp|1 iHt/
h + (t)2 |x0 i
!
ipx0 /
h
i
i p2
2 e

t V (x0 )t + O(t)
= 1
h
2m
h

2
h
!
2
1
p
i
2
=
px0 +
t + V (x0 )t + O(t) .
exp
h

2m
2h

(8)

Then the p-integral is a Fresnel integral

hx1 , t + t|x0 , ti
Z
dp ipx1 /h i(px0 + p2 t+V (x0 )t+O(t)2 )/h
2m
=
e
e
2
h
!
r
i m (x1 x0 )2
m
2
exp
V (x0 )t + O(t) .
=
2i
ht
h
2
t

(9)

The quantity in the parentheses is nothing but the classical Lagrangian

times t by identifying x 2 = (x1 x0 )2 /(t)2 .
Once Eq. (5) is shown, we use many time slices to obtain the propagator
for a finite time interval. Using the completeness relation many many times,
hxf , tf |xi , ti i =

hxf , tf |xN 1 , tN 1 idxN 1 hxN 1 , tN 1 |xN 2 , tN 2 idxN 2

dx2 hx2 , t2 |x1 , t1 idx1 hx1 , t1 |xi , ti i.

(10)

The expression is asymmetrical between x1 and x0 , but the difference between V (x0 )t
2
and V (x1 )t is approximately V 0 t(x0 x1 ) ' V 0 x(t)

and hence of higher oder.

The time interval for each factor is t = (tf ti )/N . By taking the limit
N , t is small enough that we can use the formula Eq. (5), and we
find
hxf , tf |xi , ti i =

Z NY
1

dxi ei

PN 1
i=0

L(ti )t/
h

(11)

i=1

up to a normalization. (Here, t0 = ti .) In the limit N , the integral

over positions at each time slice can be said to be an integral over all possible
paths. The exponent becomes a time-integral of the Lagrangian, namely the
action for each path.
This completes the derivation of the path integral in quantum mechanics.
As clear from the derivation, the overall normalization of the path integral
is a tricky business.
A useful point to notice is that even matrix elements of operators can be
written in terms of path integrals. For example,
hxf , tf |x(t0 )|xi , ti i =
=
=

dx(t0 )

Z
tf >t>ti

Z
tf >t>t0

dx(t0 )hxf , tf |x(t0 ), t0 ix(t0 )hx(t0 ), t0 |xi , ti i

Dx(t)eiS[x(t)]/h x(t0 )

Dx(t)eiS[x(t)]/h

t0 >t>ti

Dx(t)eiS[x(t)]/h x(t0 ),

(12)

At the last step, we used the fact that an integral over all paths from xi
to x(t0 ), all paths from x(t0 ) to xf , further integrated over the intermediate
position x(t0 ) is the same as the integral over all paths from xi to xf . The
last expression is literally an expectation value of the position in the form of
an integral. If we have multiple insertions, by following the same steps,
Z

hxf , tf |x(t2 )x(t1 )|xi , ti i =

tf >t>ti

Dx(t)eiS[x(t)]/h x(t2 )x(t1 ).

(13)

Here we assumed that t2 > t1 to be consistent with successive insertion of

positions in the correct order. Therefore, expectation values in the path
integral corresponds to matrix elements of operators with correct ordering
in time. Such a product of operators is called timed-ordered T x(t2 )x(t1 )
defined by x(t2 )x(t1 ) as long as t2 > t1 , while by x(t1 )x(t2 ) if t1 > t2 .
In general,
hxf , tf |T O1 O2 On |xi , ti i =

Z
tf >t>ti

Dx(t)O1 O2 On eiS[x(t)]/h .

(14)

Another useful point is that EulerLagrange equation is obtained by the

change of variable x(t) x(t) + x(t) with xi = x(ti ) and xf = x(tf ) held
fixed. A change of variable of course does not change the result of the integral,
and we find
Z
Z
Dx(t)eiS[x(t)+x(t)]/h = Dx(t)eiS[x(t)]/h ,
(15)
and hence
Z

iS[x(t)+x(t)]/
h

Dx(t)e

iS[x(t)]/
h

Dx(t)e

i
Dx(t)eiS[x(t)]/h S = 0. (16)
h

Recall (see Note on Classical Mechanics I)

S = S[x(t) + (t)] S[x(t)] =

d L
L

x(t)dt.
x dt x

(17)

Therefore,
Z

Dx(t)e

iS[x(t)]/
h

d L
L

x(t)dt = 0.
x dt x

iZ
h

(18)

Because x(t) is an arbitrary change of variable, the expression must be zero

at all t independently,
Z

Dx(t)e

iS[x(t)]/
h

d L
L

x dt x

i
h

= 0.

(19)

Therefore, the EulerLagrange equation must hold as an expectation value,

nothing but the Ehrenfests theorem.

Schr
odinger Equation from Path Integral

Here we would like to see that the path integral contains all information we
need. In particular, we rederive Schrodinger equation from the path integral.
Let us first see that the momentum is given by a derivative. Starting
from the path integral
hxf , tf |xi , ti i =

Dx(t)eiS[x(t)]/h ,

(20)

we shift the trajectory x(t) by a small amount x(t)+x(t) with the boundary
condition that xi is held fixed (x(ti ) = 0) while xf is varied (x(tf ) 6= 0).
Under this variation, the propagator changes by
hxf + x(tf ), tf |xi , ti i hxf , tf |xi , ti i =

hxf , tf |xi , ti ix(tf ).

xf

(21)

On the other hand, the path integral changes by

Z

Dx(t)eiS[x(t)+x(t)]/h

Dx(t)eiS[x(t)]/h =

Dx(t)eiS[x(t)]/h

iS
.
h

(22)

Recall that the action changes by (see Note on Classical Mechanics II)
S = S[x(t) + x(t)] S[x(t)]
!
Z tf
L
L
dt
=
x +
x
x
x
ti
x
!
L
L  f Z tf
d L
dt
=
x +

x.
x xi
x dt x
ti

(23)

The last terms vanishes because of the equation of motion (which holds as
an expectation value, as we saw in the previous section), and we are left with
S =

L
x(tf ) = p(tf )x(tf ).
x

(24)

By putting them together, and dropping x(tf ), we find

Z

i
hxf , tf |xi , ti i = Dx(t)eiS[x(t)]/h p(tf ).
xf
h

(25)

This is precisely how the momentum operator is represented in the position

space.
Now the Schrodinger equation can be derived by taking a variation with
respect to tf . Again recall (see Classical Mechanics II)
S
= H(tf )
tf

(26)

after using the equation of motion. Therefore,

Z

i
hxf , tf |xi , ti i = Dx(t)eiS[x(t)]/h H(tf ).
tf
h

(27)

If

p2
H=
+ V (x),
(28)
2m
the momentum can be rewritten using Eq. (25), and we recover the Schrodinger
equation,

ih hxf , tf |xi , ti i =
t
2m

h

i xf

!2

+ V (xf ) hxf , tf |xi , ti i.

(29)

In other words, the path integral contains the same information as the
conventional formulation of the quantum mechanics.

6
6.1

Examples
Free Particle

Here we calculate the path integral for a free particle in one-dimension. We

need to calculate
hxf , tf |xi , ti i =

Dx(t)ei

m 2
x dt/
h
2

(30)

over all paths with the boundary condition x(ti ) = xi , x(tf ) = xf .

6.1.1

The sloppy way

First, we do it somewhat sloppily, but we can obtain dependences on important parameters nonetheless. We will come back to much more careful
calculation with close attention to the overall normalization later on.
The classical path is
xc (t) = xi +

xf xi
(t ti ).
tf ti

(31)

We can write x(t) = xc (t) + x(t), where the left-over piece (a.k.a. quantum
fluctuation) must vanish at the initial and the final time. Therefore, we can
expand x(t) in Fourier series
x(t) =

an sin

n=1

n
(t ti ).
tf ti

(32)

The integral over all paths can then be viewed as integrals over all an ,
Z

Dx(t) = c

Z Y

dan .

(33)

n=1

The overall normalization factor c that depends on m and tf ti will be fixed

later. The action, on the other hand, can be calculated. The starting point
is

xf xi X
n
n
x = x c + x =
+
an cos
(t ti ).
(34)
tf ti
tf ti
n=1 tf ti
Because different modes are orthogonal upon t-integral, the action is
S=

1 (n)2 2
m (xf xi )2 m X
+
a .
2 tf ti
2 n=1 2 tf ti n

(35)

The first term is nothing but the classical action. Then the path integral
reduces to an infinite collection of Fresnel integrals,
c

Z Y

"

im
dan exp
h
2
n=1

1 (n)2 2
(xf xi )2 X
+
an
tf ti
n=1 2 tf ti

!#

(36)

We now obtain
hxf , tf |xi , ti i = c

Y
n=1

i m 1 (n)2

h 2 2 tf ti

!1/2

i m (xf xi )2
exp
.
h
2 tf ti
"

(37)

Therefore, the result is simply

i m (xf xi )2
hxf , tf |xi , ti i = c (tf ti ) exp
.
h
2 tf ti
"

(38)

The normalization constant c can depend only on the time interval tf ti ,

and is determined by the requirement that
Z

dxhxf , tf |x, tihx, t|xi , ti i = hxf , tf |xi , ti i.

(39)

And hence
c0 (tf t)c0 (t ti

v
u
u 2i
h(tf
)t

t)(t ti )
= c0 (tf ti ).
m(tf ti )

(40)

This argument does not eliminate the possible factor ei0 (tf ti ) , which would correspond to a zero point in the energy h
0 . Mahiko Suzuki pointed out this problem to
me. Here we take the point of view that the normalization is fixed by comparison to the
conventional method.

10

We therefore find

m
,
(41)
2i
ht
recovering Eq. (2.5.16) of Sakurai precisely.
In order to obtain this result directly from the path integral, we should
have chosen the normalization of the measure to be
0

c (t) =

Dx(t) =
6.1.2

Y
m
2ih(tf ti ) n=1

m
n
dan .
2i
h(tf ti ) 2

(42)

The careful way with discretized time

Now we do it much more carefully with close attention to the overall normalization. The
path integral for the time interval tf ti is divided up into N time slices each with
t = (tf ti )/N ,
1
 m N/2 Z NY
K=
dxn eiS/h ,
(43)
2iht
n=1
where
S=

N 1
1 X (xn+1 xn )2
m
.
2 n=1
t

(44)

The point is that this is nothing but a big Gaussian (to be more precise, Fresnel) integral,
because the action is quadratic in the integration variables x1 , , xN 1 . To make this
point clear, we rewrite the action into the following form:

xN
0

0 1m
m
1m 2

2
(x + x0 )
(xN 1 , xN 2 , xN 3 , , x2 , x1 ) . +

S=
.. 2 t
2 t N
t

0
x0

2 1 0 0
0
xN 1
1 2 1 0

xN 2
0 1 2 0
xN 3
0

(xN 1 , xN 2 , xN 3 , , x2 , x1 ) .

..
..
..
..
..
..
..

.
.
.
.
.
.

0
0 2 1
x2
0
0
0 1 2
x1
(45)
Comapring to the identity Eq. (111),
Z Y
N

dxn e 2 x

AxxT y

= (2)N/2 (detA)1/2 e+ 2 y

n=1

11

A1 y

(46)

we identify

i 2m

h t

1
1/2
0
..
.

1/2
1
1/2
..
.

0
1/2
1
..
.

..
.

0
0
0
..
.

0
0

0
0

0
0

1
1/2

0
0
0
..
.

1/2
1

(47)

This is nothing but the matrix KN 1 in Eq. (111) with a = 1/2 up to a factor of hi 2m
t .
One problem is that, for a = 1/2, + = = 1/2 and hence Eq. (116) is singular.
Fortunately, it can be rewritten as
detKN 1 =

N
N
+
1
N 2
N 2
N 1
= N
+ +
+ + + +
+
=N
+
+

 N 1
1
. (48)
2

Therefore,

detA =

2m
i
ht

N 1
N

1
2N 1

=N

 m N 1
.
iht

(49)

Therefore, the prefactor of the path integral is



m N/2
(2)(N 1)/2 (detA)1/2
2iht
 m N/2
 m (N 1)/2
=
(2)1/2 N 1/2
i
ht
iht
1/2



1/2
m
m
=
.
=
2i
hN t
2ih(tf ti )

(50)

The exponent is given by

1  m 2

i 1 im 2
(x + x20 ) +
2
h
ht N
2

xN
0
0
..
.

1
(xN , 0, 0, , 0, x0 )A

0
x0

(51)

We need only (1, 1), (N 1, N 1), (1, N 1), and (N 1, 1) components of A1 .

(A1 )1,1 = (A1 )N 1,N 1 =

iht N 1
iht detKN 2
=
.
2m detKN 1
m
N

On the other hand,

(A1 )1,N 1 = (A1 )N 1,1

12

(52)

i
ht
1

(1)N 2 det
2m detKN 1

1/2
0
..
.

1
1/2
..
.

1/2
1
..
.

..
.

0
0
..
.

1/2

N 1

iht 2
1
2m N 2N 2
i
ht
=
.
Nm
We hence find the exponent
=

(53)

xN
0
0
..
.

i1m 2

1
2
(xN , 0, 0, , 0, x0 )A
(xN + x0 ) +

h 2 t

2 t

0
x0

2 

i1m 2
1 im
iht N 1 2
iht
=
(xN + x20 ) +
(xN + x20 ) + 2
xN x0
h 2 t

2 ht
m
N
Nm
1
im
i1m 1 2
(x + x20 )
2xN x0
=
h 2 t N N

2 hN t
i m (xN x0 )2
=
h 2

N t
i m (xf xi )2
=
.
(54)
h 2 tf ti

1  m 2

Putting the prefactor Eq. (50) and the exponent Eq. (81) together, we find the propagator

1/2
2
i m (xf xi )
m
h
2
tf ti
K=
.
(55)
e
2ih(tf ti )
We recovered Eq. (2.5.16) of Sakurai without any handwaving this time!

6.2

Harmonic Oscillator

Now that we have fixed the normalization of the path integral, we calculate
the path integral for a harmonic oscillator in one-dimension. We need to
calculate
Z
R m 2 1 2 2
hxf , tf |xi , ti i = Dx(t)ei ( 2 x 2 m x )dt/h ,
(56)
over all paths with the boundary condition x(ti ) = xi , x(tf ) = xf . The
classical path is
xc (t) = xi

sin (tf t)
sin (t ti )
+ xf
.
sin (tf ti )
sin (tf ti )
13

(57)

The action along the classical path is


(x2 + x2f ) cos (tf ti ) 2xi xf
1
1
2 2
Sc =
x m x dt = m i
.
2
2
2
sin (tf ti )
ti
(58)
The quantum fluctuation around the classical path contributes as the
O(
h) correction to the amplitude, relative to the leading piece eiSc /h . We
expand the quantum fluctuation in Fourier series as in the case of free particle.
Noting that the action is stationary with respect to the variation of x(t)
around xc (t), there is no linear piece in an . Because different modes are
orthogonal upon t-integral, the action is
Z tf 
m

m
S = Sc +
n=1 2

1 2
(n)2
2 (tf ti )
a .
tf ti
2 n
!

(59)

Therefore the path integral is an infinite number of Fresnel integrals over an

using the measure in Eq. (42)
hxf , tf |xi , ti i = e

n=1

= eiSc /h

iSc /
h

m
2i
h(tf ti )
"

n Z
im
m

dan exp
2i
h(tf ti ) 2
h
2

Y
m
(tf ti )
1
2i
h(tf ti ) n=1
n

(n)2
1 2
a
2 (tf ti )
tf ti
2 n
!

!2 1/2

(60)

Now we resort to the following infinite product representation of the sine

function
!

Y
x2
sin x
1 2 =
.
(61)
n
x
n=1
We find
hxf , tf |xi , ti i = eiSc /h
= e

iSc /
h

v
u
u
t

(tf ti )
sin (tf ti )

m
2i
h(tf ti )

m
2i
h sin (tf ti )

(62)

This agrees with the result from the more conventional method as given in
Sakurai Eq. (2.5.18).
One can obtain exactly the same result using the discretized time as
shown in the Appendix.
14

Partition Function

The path integral is useful also in statistical mechanics to calculate partition

functions. Starting from a conventional definition of a partition function
Z=

eEn ,

(63)

where = 1/(kB T ), we rewrite it in the following way using completeness

relations in both energy eigenstates and position eigenstates.
Z=

hn|eH |ni =

dx

dx

hn|xihx|eH |ni

hx|e

|nihn|xi =

dxhx|eH |xi.

(64)

The operator eH is the same as eiHt/h except the analytic continuation

t i = ih. Therefore, the partition function can be written as a path
integral for all closed paths, i.e., paths with the same beginning and end
points, over a time interval i
h. For a single particle in the potential
V (x), it is then

Z=

1 I h m
d
Dx( ) exp
h
0
2

!2

+ V (x) .

(65)

The position satisfies the periodic boundary condition x(

h) = x(0).
For example, in the case of a harmonic oscillator, we can use the result
in the previous section
hxf , tf |xi , ti i
s

(x2i + x2f ) cos (tf ti ) 2xi xf

m
i1
exp
m
,
2i
h sin (tf ti )
h
2
sin (tf ti )
(66)
"

and take xf = xi = x, tf ti = ih
s

hx|e

|xi =

"

m
11
2 cosh
h 2
exp mx2
.
2
h sinh
h
h
2
sinh h

15

(67)

Then integrate over x,

Z =

dxhx|eH |xi

m
2
sinh h
h
2h sinh h
m 2 cosh
h 2

1
4 sinh h/2
1
=
2 sinh h/2
eh/2
=
1 eh
=

e(n+ 2 )h .

(68)

n=0

Of course, we could have redone the path integral to obtain the same result.
The only difference is that the classical path you expand the action around
is now given by sinh rather than sin.
Let us apply the path-integral representation of the partition function to
a simple harmonic oscillator. The expression is

Z=

1 I h m
d
Dx( ) exp
h
0
2

!2

m
+ 2 x2 .
2

(69)

We write all possible paths x( ) in Fourier series as

x( ) = a0 +

an cos

n=1

X
2n
2n
+
bn sin
,
h

n=1

(70)

satisfying the periodic boundary condition x(

h) = x(0). Then the exponent
of the path integral becomes

1 I h m
d
h
0
2

!2

m
+ 2 x2
2

m 2 1 X
2n

=
a0 +
2
2 n=1
h

16

!2

(a2n

b2n ) .

(71)

The integration over all possible paths can be done by integrating over the
Fourier coefficients an and bn . Therefore the partition function is
Z=c

m 2 2 1 X
2n

da0
dan dbn exp
a0 +
2
2 n=1
h

n=1

!2

+ 2 (a2n + b2n ) .
(72)

This is an infinite collection of Gaussian integrals and becomes

2 Y
4 2n
Z=c
2
m n=1 m
h

!2

+ 2

(73)

Here, c is a normalization constant which can depend on because of normalization of Fourier modes, but not on Lagrangian parameters such as m or
. Now we use the infinite product representation of the hyperbolic function

Y
n=1

x2
1+ 2
n

sinh x
.
x

(74)

The partition function Eq. (82) can be rewritten as

h
2
h

2 Y
Z = c
1+
2
2
m n=1 mn
2n

2
sinh h
/2

h/2
e
= c0
,
1 eh

!2 1

= c0

(75)

where c0 is an overall constant, which does not depend on , is not important

when evaluating various thermally averages quantities. If you drop c0 , this
is exactly the partition function for the harmonic oscillator, including the
zero-point energy.

Schr
odinger Field

We now apply the path integral to field theory. As a starter, we compute

the partition function of the non-relativistic identical bosons. This calculation shows non-trivially that the quantized Schrodinger field theory indeed
17

contains multi-body states automatically and is also useful from a practical

point of view. It is also meant to be a demonstration that path integral in
field theory makes sense.
The path integral in the particle quantum mechanics is given by the
summation over all possible paths in the configuration space xi (t). In the field
theory, (~x) is the canonical coordinate, and it can follow various paths
(~x, t). Therefore, the path integral in field theory is a summation over all
possible field configurations (~x, t). This discussion defines the path integral
Z

D(~x, t)D (~x, t)eiS/h ,

(76)

where the action is that of the field theory, such as

S=

d~xdt i
h +
2m

(77)

in the free case.

To calculate the partition function, we go to the imaginary time t = i ,
Z=

1 Z h Z
h2
d d~x h
+

D(~x, )D (~x, ) exp

h
0
2m
"

!#

(78)
We impose the periodic boundary condition (~x, t + h
) = (~x) for =
1/kT .
Because of the periodic boundary condition in and also in space due to
the box normalization, we can expand in Fourier series both in the imaginary
time as well as in space,
1 X
zp~,n ei~p~x/h e2in /h
3/2
L p~,n

(79)

1 X i~p~x/h 2in /h
z e
e
.
L3/2 p~,n p~,n

(80)

(~x, ) =
and simiarly for
(~x, ) =

Then the exponent in the partition function Eq. (97) is

X
1 Z h Z
h2
p~2
d d~x h
+
=
2in +
zp~,n zp~,n .

h
0
2m
2m
p
~,n
(81)
!

18

Therefore, the path integral is simply a product of many many Gaussian

integrals
Z=

YZ

2 /2m)z z
~,n
p
~,n p

dzp~,n dzp~,n e(2in+~p

p
~,n

Y
p
~,n

.
2in + ~p2 /2m

(82)

Now we use the infinite product representation of the hyperbolic functions

Y
n=1

Y
n=1

x2
1+ 2
n

sinh x
x

x2
1+
(2n 1)2

= cosh

x
.
2

(83)

We use the first identity here. The partition function Eq. (82) can be rewritten as
Z =

Y
p
~

Y
p
~

2
2
~p /2m n=1 (2n) + (~p2 /2m)2

(~p2 /2m/2)2
1
+
~p2 /2m n=1 (2n)2
n2

Y
p
~

= c

Y
p
~

2
n=1 (2n)

!1

~p2 /4m

~p2 /2m sinh ~p2 /4m

2
sinh ~p2 /4m
2

= c

Y
p
~

e~p /4m
.
1 e~p2 /2m

(84)

where c is an (infinite) overall constant, which is not important when evaluating various thermally averages quantities. This is a grand partition function
which contains summation over all possible number of particles but with no
chemical potential. One small difference from the conventional calculation is
that the path integral automatically includes the zero-point energy p~2 /4m,
which results in the exponential factor in the numerator. The reason why
there is this zero-point energy is because the Fock space is a collection of
(infinite number of) harmonic oscillators and each harmonic oscillator has
the zero-point energy. But the zero-point energy does not lead to any physical consequences (while the zero-point fluctuation does) and we can always
19

shift the energy of the system by an infinite constant p~ p~2 /4m to remove
the zero-point energy in the expression.
The inclusionR of the chemical potential is obvious. Because the number
operator is N = d~x (~x)(~x) and the grand canonical emsenble is summed
with a factor of e(EN ) , the path integral Eq. (97) is modified to
P

Z =

D(~x, )D (~x, )

h2
1 Z h Z
exp
d d~x h
+

h
0
2m
"

!#

(85)

The result of the path integral is then

2

Z=c

Y
p
~

e(~p /2m)/2
.
1 e(~p2 /2m)

(86)

Note that Z = e with = pV = F N for the grand partition

function. Namely,
2

= kT

X
p
~

"

X
p
~

1
2

e(~p /2m)/2
ln
1 e(~p2 /2m)
p~2
2
+ kT ln(1 e(~p /2m) ) .
2m
!

(87)

The thermodynamic quantities are related by

F = U T S,
= F N,

(88)
(89)

and have the following differentials

dU = T dS pdV + dN,
dF = SdT pdV + dN,
d = SdT pdV N d.

(90)
(91)
(92)

Therefore, the number is given by

2

X
X 1

1
e(~p /2m)
1

=
+
=
+
.
N =
2
2
T,V
2 1 e(~p /2m)
2 e(~p /2m) 1
p
~
p
~
(93)


"

20

Apart from the zero-point piece 1/2, this is the familiar expression for the
occupation number for Bose-Einstein distribution.
The thermal average energy U is given by the Legendre transform U =
+ T S + N with S = /T |V, . We obtain
U=

X
p
~

1
p~2 1
+ (~p2 /2m)
.
2m 2 e
1


(94)

This must be a famililar expression to you, again except the zero-point energy
term that can be dropped without changing physical content.

Non-relativistic Fermions

In the case of fermions, their fields are supposed to satisfy canonical anticommutation relation. In the classical limit, they become anti-commuting
numbers which are called Grassmann odd. Clearly, a Grassmann-odd
classical Fermi field is not an observable as no measurement would yield a
Grassmann-odd number as a result. However, this issue is not completely
academic because the path integral requires the classical integration variable:
we need to integrate over Grassmann-odd Fermi field (~x, t).
What is the definition of an integral over a Grassmann-odd number? I
do not go into details, but rather give you a consistent definition for our use:
the integration is the same as the differentiation! To make it explicit what
I mean, let me take one Grassmann-odd number (no ~x, t dependence).
Since it anti-commutes with itself, {, } = 0, we find 2 = 0. Namely, any
higher powers than one give vanishing results. Therefore, any function of
can always be expanded as
f () = f (0) + f 0 (0).

(95)

The definition of the Grassmann integral is

Z

df () = f 0 (0).

(96)

One could have tried to take f (0) instead as the definition of the integral
rather than f 0 (0); however, it would violate the shift invariance of the integration variable 0 = +. We would like toRretain this property
to be analR

ogous to the bosonic integral which satisfies dxf (x) = dxf (x + a).
21

To see this, let us imagine

defining the integral to pick up the piece
f (0) inR
R
0
stead of f (0). Then dfR() = f (0) by definition.
However, df ( + )
R
for a constant would be df ( + ) = d(f (0) + ( + )f 0 (0)). Since
this (wrong) definition would pick all independent pieces, the result would
be f (0) + f 0 (0), which is different from the result without the shift. On the
other hand, picking up the piece proportional to yields f 0 (0) in both cases.
The partition function is given by the same expression as the bosonic case
Z=

h2
1 Z h Z
d d~x h
+

D(~x, )D (~x, ) exp

h
0
2m
"

!#

(97)
However, there is another peculiarity with Fermi field: we need to take the
anti-periodic boundary condition (
h) = (0) instead of the periodic
boundary condition. This can be shown by carefully working out the trace
TreH in terms of fermionic path integral (see the appendix). We do not
go into the detail of the derivation here, but proceed with the anti-periodic
boundary condition. Fourier expansion is then modified
(~x, ) =

1 X
zp~,n ei~p~x/h ei(2n1) /h .
L3/2 p~,n

(98)

The partition function is

Z =

YZ

2 /2m))z z
~,n
p
~,n p

dzp~,n dzp~,n e(i(2n1)+(~p

p
~,n

"

Y
p
~,n

p~2

i(2n 1) +
2m

!#

(99)

where we used the rule for the Grassmann integrals after Taylor expanding
the exponential up to order zp~,n zp~,n for each p~, n. This time we use the
formula
!

Y
x2
x
1+
=
cosh
,
(100)
(2n 1)2
2
n=1
and find
Z=c

Y
p
~

cosh
2


Y
p~2
2
2
=c
e(~p /2m) + 1 e(~p /2m)/2 .
2m
p
~
!

22

(101)

The energy expectation value is hence




X p
~2

1
1

hEi =
ln Z  =
+ (~p2 /2m)
,

2m
2 e
+1
p
~

(102)

which must be familiar to you except the zero-point energy contribution. It

is interesting to note that the zero-point energy is negative for fermions.

Path Integrals in Phase Space

As we saw above, the normalization of the path integral is a somewhat tricky issue. The
source of p
the problem was when we did p integral in Eq. (9), that produced a singular
prefactor m/2i
ht. On the other hand, if we do not do the p integral, and stop at
Z
dp ipx1 /h i(px0 + p2 t+V (x0 )t+O(t)2 )/h
2m
e
e
hx1 , t + t|x0 , ti =
2h
Z
dp i(p x1 x0 p2 V (x0 ))t
t
2m
e
,
(103)
=
2h
we do not obtain any singular prefactor except 1/2h for every momentum integral. The
path integral can then be written as
Z
hxf , tf |xi , ti i = Dx(t)Dp(t)eiS[x(t),p(t)]/h ,
(104)
where the action is given in the phase space
Z tf
S[x(t), p(t)] =
(px H(x, p)) dt.

(105)

ti

The expression (104) is an integral over all paths in the phase space (x, p). It is understood
that there is one more p integral than the x integral (albeit both infinite) because of the
number of times the completeness relations are inserted.
Using the path integral in the phase space, we can work out the precise normalization
of path integrals up to numerical constants that depend on t, but not on m, or any other
parameters in the Lagrangian. For instance, for the harmonic oscillator, we expand
s

2 X
n
xn sin
(t ti ),
(106)
x(t) = xc (t) +
tf ti n=1
tf ti
s
s

1
2 X
n
p(t) = pc (t) +
p0 +
pn cos
(t ti ).
(107)
tf ti
tf ti n=1
tf ti

If you have both bosons and fermions with the same mass, the zero-point enegies
precisely cancel. This is what happens with supersymmetry.

23

The action is then

p2
1X
S = Sc 0
(pn , xn )
2m 2 n=1

1
m
n
tf ti

n
tf ti
2

!

pn
xn


,

(108)

where Sc is the action for the classical solution Eq. (58). Using the integration volume
Z
Z

dp0 Y dxn dpn

,
(109)
Dx(t)Dp(t) = c(tf ti )
2h n=1 2h
we find the path integral to be

1/2
2mh
1
2h
i



1/2

1/2
Y 1
1
2mh
2h

. (110)
2h
i
im 2 (n/(tf ti ))2
n=1

hxf , tf |xi , ti i = eiSc /h c(tf ti )

Up to an overall constant c(tf ti ) that depends only on tf ti , both m and dependences

are obtained correctly.

A Few Useful Identities

To carry out the integral, a few identities would be useful. For N 1-dimensional column
vectors x, y, and a symmetric matrix A,
Z Y
N

dxn e 2 x

AxxT y

= (2)N/2 (detA)1/2 e+ 2 y

A1 y

(111)

n=1

It is easy to prove this equation first without the linear term y by diagonalizatin the matrix
A = OT DO, where O is a rotation
Q matrix and D = diag(1 , , N ) has N eigenvalues.
Under a rotation, the measure
dxn is invariant because the Jacobian is detO = 1.
Therefore, after rotation of the integration variables, the integral is
N Z
N r
N
Y
Y
Y
2
n x2n /2
dxn e
=
= (2)N/2 (
(112)
n )1/2

n
n=1
n=1
n=1
The product of eigenvalues is nothing but the determinant of the matrix detA = detOT detDdetO =
Q
N
n=1 n , and hence
Z Y
N

dxn e 2 x

Ax

= (2)N/2 (detA)1/2 .

(113)

n=1

In the presence of the linear term, all we need to do is to complete the square in the
exponent,
1
1
1
xT Ax xT y = (xT yT A1 )A(x A1 y) + yT A1 y,
2
2
2

24

(114)

and shift x to eliminate A1 y, and diagonalize

For a N N matrix of the form

1 a 0
a 1 a

0 a 1

KN = .
..
..
..
.
.

0
0
0
0
0
0
one can show
detKN =

+1
+1
N
N
+

,
+

A. This proves the identity Eq. (111).

..
.

0
0
0
..
.

0
0
0
..
.

a
1

1
a

(115)

p
1
(1 1 4a2 ).
2

(116)

This is done by focusing on the top two-by-two block. The determinant that picks
the (1, 1) element 1 comes with the determinant of remaining (N 1) (N 1) block,
namely detKN 1 . If you pick the (1, 2) element a, it is necessary to also pick the (2, 1)
element which is also a and the rest comes from the remaining (N 2) (N 2) block,
namely detKN 2 . Therefore, we find a recursion relation
detKN = detKN 1 a2 detKN 2 .
We rewrite this recursion relation as

 
detKN
1
=
detKN 1
1

a2
0



(117)

detKN 1
detKN 2


.

(118)

The initial conditions are detK1 = 1, detK2 = 1 a2 . Using the recursion relation
backwards, it is useful to define detK0 = 1. We diagonalize the matrix

 



1
1 a2
+
+ 0
1
=
.
(119)
1
0
1
1
0
1 +
+
Hence,


detKN
detKN 1

1




+
0

+
1

+
1



1
+

1
+

N
+
N 1
+

1
+

+1
+1
N
N
+

N
N
+

1
N

25

N 1 
N 1
+
0


1
1
0
1
N

+


1
+

1
1 + +

detK1
detK0


.

(120)

Harmonic Oscillator with Discretized Time

We perform path integral for a harmonic oscillator again using discretized time. This way,
we can pay closer attention to the normalization by going through the same steps as for a
free particle using discretized time slices. The action is
S=



N
1
N 1
X
1
1
1 2 1 2
1 X (xn+1 xn )2
m
m 2
x2n m 2
x0 + xN .
2 n=1
t
2
2
2
2
n=1

(121)

The last term Ris there to correctly acccount for all N 1 time slices for the potential
energy term V dt, but to retain the symmetry between the initial and final positions,
we took their average. Once again, this is nothing but a big Fresnel (complex Gaussian)
integral:

xN
0

0 1m
m
1m 2

2
(x + x0 )
(xN 1 , xN 2 , xN 3 , , x2 , x1 ) . +

S=
.. 2 t
2 t N
t

0
x0

xN 1
2 1 0 0
0

1 2 1 0
0
xN 2

0 1 2 0

0 xN 3

(xN 1 , xN 2 , xN 3 , , x2 , x1 ) .

.
.
.
..
..
..
..
..
..

..

.
.
.

0
0
0 2 1 x2
0
0
0 1 2
x1

xN 1
xN 2

xN 3
1
1

m 2 (xN 1 , xN 2 , xN 3 , , x2 , x1 )
t m 2 (x2N + x20 )t
..

2
4
.

x2
x1
(122)
Comapring to the identity Eq. (111),
Z Y
N

dxn e 2 x

AxxT y

= (2)N/2 (detA)1/2 e+ 2 y

n=1

26

A1 y

(123)

we identify

i
i
2m

A = m 2 tIN 1

h t

1
1/2
0
..
.

1/2
1
1/2
..
.

0
1/2
1
..
.

..
.

0
0
0
..
.

0
0

0
0

0
0

1
1/2

0
0
0
..
.

1/2
1

(124)

This is nothing but the matrix KN 1 in Eq. (111) up to an overall factor of

i 2m
i
m
+ m 2 t =
(2 2 (t)2 )
h t

h
iht

with
a=

1
2

2m
t

2m/t
1
=
.
2
2 2 (t)2
m t

(125)

(126)

Using Eq. (116),

1
=
2

s
1

4
1
(2 2 (t)2 )2

!
(127)

Therefore,
detA =

 m N 1
N N

.
(2 2 (t)2 )N 1 +
i
ht
+

(128)

Therefore, the prefactor of the path integral is



m N/2
(2)(N 1)/2 (detA)1/2
2i
ht
r
 N
1/2


+ N
m

2
2 (N 1)/2
=
2 (t)
.
2i
ht
+

(129)

Now we use the identity

lim

1+

x N
= ex .
N

We find in the limit of t = (tf ti )/N and N ,

s
!
1
4
1
=
1 1
= (1 it) + O(N 3 ),
2
(2 2 (t)2 )2
2

(130)

(131)

and hence
N
N
+
+

1 ei(tf ti ) ei(tf ti )
+ O(N 1 )
2N
it
N sin (tf ti )
+ O(N 1 ).
2N 1 (tf ti )

27

(132)

In the same limit,

2 2 (t)2

(N 1)/2


(N 1)/2
1
2(N 1)/2 1 2 (t)2
2

2(N 1)/2 + O(N 1 ).

(133)

Therefore, the prefactor Eq. (129) is

r

1/2 r
N sin (tf ti )
m
m
=
2(N 1)/2
.
N
1
2i
ht
2
(tf ti )
2ih sin (tf ti )

(134)

The exponent is given by

i 1 im 2
i1
1
(x + x20 )
m 2 (x2N + x20 ) +
2
h
ht N
h4

im
t
h

2

xN
0
0
..
.

(xN , 0, 0, , 0, x0 )A1

0
x0

(135)

We need only (1, 1), (N 1, N 1), (1, N 1), and (N 1, 1) components of A1 .

(A1 )1,1 = (A1 )N 1,N 1
=

N 1
N 1
detKN 2
iht
i
ht
(2 2 (t)2 )1
(2 2 (t)2 )1 + N
=
m
detKN 1
m
+ N

(136)
On the other hand,
(A1 )1,N 1 = (A1 )N 1,1

a 1
0 a
..
..
.
.
0
0

i
ht
1

(2 2 (t)2 )1
(1)N 2 det
m
detKN 1

+ N 2
i
ht
(2 2 (t)2 )1 N
a
m
+ N

i
ht 1 2N 1 (tf ti )
1
iht (tf ti )
2
=
.
m
N sin (tf ti ) 2N 2
N m sin (tf ti )

a 0
1 0
..
..
..
.
.
.
0 a

(137)

We hence find the exponent

i
h

1  m 2

1 im
1
m 2 (x2N + x20 ) +
2
ht 4
2

28

xN
0
0
..
.

1
(xN , 0, 0, , 0, x0 )A

0
x0

i
h



2
1 im
1
1 im
iht
m 2 (x2N + x20 ) +

2
ht 4
2 ht
m

(2

N 1
(t)2 )1 + N
+

N 1

2 (tf ti )
(x2N + x20 ) +
xN x0
N
N
sin (tf ti )

(138)
Taking the limit N , the second term is obviously regular, and the last term in the
parentheses is also:

2
i
ht 2 (tf ti )
i
m
1 im
xN x0 =
xN x0 .
(139)
2
ht
m N sin (tf ti )
h sin (tf ti )
The other terms are singular and we need to treat them carefully.
N 1
N 1
1 im
i1m 2
(xN + x20 )
(2 2 (t)2 )1 + N
(x2N + x20 )
h 2 t

2h
t
+ N

!
N 1
N 1

im

=
1 (2 2 (t)2 )1 + N
(x2N + x20 ).
2ht
+ N

(140)

We need O(N 1 ) term in the parentheses. The factor (2 2 (t)2 )1 = 21 + O(N 2 )

is easy. Using Eq. (131), the next factor is
1
1
N
N

+
N
N
+

2N (1 + it)N 1 (1 it)N 1 + O(N 2 )

(1 + it)N (1 it)N + O(N 2 )

2N 1

(141)

Therefore, dropping O(N 2 ) corrections consistently,

1 (2 2 (t)2 )1

1
N 1
N

+
N
N
+

(1 + it)N (1 it)N (1 + it)N 1 + (1 it)N 1

(1 + it)N (1 it)N
(1 + it)N 1 + (1 it)N 1
= it
(1 + it)N (1 it)N
(tf ti ) 2 cos (tf ti )
1 (tf ti ) cos (tf ti )
= i
=
.
N
2i sin (tf ti )
N
sin (tf ti )

(142)

Collecting all the terms, the exponent is

i
Sc

=
=

i
m
im 1 (tf ti ) cos (tf ti ) 2
xN x0 +
(xN + x20 )
sin (tf ti )
h
2ht N
sin (tf ti )

i m (x2N + x20 ) cos (tf ti ) 2xN x0

h 2

sin (tf ti )

which is the same as the classical action Eq. (58).

29

(143)

Anti-periodic Boundary Condition

Here is the derivation why fermion fields must satisfy the anti-periodic boundary condition.
It turns out it is a somewhat long story.
Let us start with a Lagrangian
L = ih h .

(144)

(t) depends only on time but not on space. The classical EulerLagrange equation is
ih = 0,

(145)

(t) = (0)eit .

(146)

with an obvious solution

To see what this Lagrangian describes quantum mechanically, note that the can be regarded as the canonical coordinate, while its canonically conjugate momentum is L/ =
Therefore, we find the canonical anti-commutation relation
i.
{, } = 1.

(147)

This defines a two-dimensional Hilbert space spanned by

|1i = |0i,

|0i,

(148)

where the state |0i is defined by the requirement |0i = 0. The Hamiltonian of the system
is


1

.
(149)
H=h

2
(The zero-point energy was added to be consistent with the result from the path integral.)
In other words, this is a two-state system with energies 21 h.
R
Now we need to find the decomposition of unity analogous to |xidxhx| = 1 to be
used to derive the path integral for . Note that the conventional path integral was
derived using the position eigenstates x|qi = q|qi. Likewise, we need the eigenstates of
the operator . It can be obtained in analogy to coherent states
|i = |i,

(150)

where the coherent state |i is defined by

|i = |0i |1i.

(151)

Note that is a Grassman number, and anti-commutes with the annihilation operator
even though it is not an operator. That guarantees that
|i = |0i |1i = |1i = |0i = |i,

(152)

where we used 2 = {, }/2 = 0 at the very last step. This coherent state, however, is
not normalized, because

h|i = 1 + = e ,
(153)

30

where the last expression is more convenient for later purposes. Now it is easy to see that
the following decomposition of unity holds:
Z

1 = d d|ie h|.
(154)
Starting from the r.h.s., we can check that
Z
Z

d d|ie
h| =
d d(|0i |1i)e (h0| h1|)
Z
=
d d ( |0ih0| + |1ih1|)
= |0ih0| + |1ih1|.

(155)

Here and
R below, we stick
R to the convention that the integrals are done from right to left,
so that dd = d = 1.
We need another quantity, h 0 |eiHt/h |i, analogous to hx0 |eiHt/h |xi for small t for
the ordinary path integral as in Eq. (9). This can be calculated easily
h 0 |eiHt/h |i = h 0 |ei

|i = h 0 |ei

|i = ei

t 0

(156)

for small t.
Now we can keep inserting the decomosition of unity between the initial and the final
states.
h(t)|eiHt/h |(0)i
Z

iHt/
h
= lim
hN |eiHt/h |N 1 ieN 1 N 1 dN
|N 2 i
1 dN 1 hN 1 |e
N

1 1
eN 2 N 2 dN
d1 d1 h1 |eiHt/h |0 i
2 dN 2 e
Z

iN
1 N 2 eN 1 N 2
eiN N 1 eN N 1 eN 1 N 1 dN
lim
1 dN 1 e

1 1
eN 2 N 2 dN
d1 d1 ei1 0 e1 0
2 dN 2 e

 Z
Z
i

dt( ih h ) .
d (t)d(t) exp
h

(157)

Note that (t) = N

and (0) = 0 are not integrated over. This defines the path integral.
Now comes the question of the trace. The trace of an operator O is defined by

TrO = h0|O|0i + h1|O|1i.

(158)

The point is to rewrite it in terms of an integral. It is easy to show that

Z
Z

d dh|O|ie =
d d(h0| h1|)O(|0i |1i)e
Z
=
d d( h0|O|0i + h1|O|1i)
= h0|O|0i h1|O|1i.

31

(159)

This is off by a sign for the second term for the trace we wanted to calculate. We can
avoid this if we had taken not h| but rather h|:
Z

d dh|O|ie = h0|O|0i + h1|O|1i.

(160)
This is the origin of the anti-periodic boundary condition.
We finally calculate hN |eH |0 i in terms of a path integral, and to obtain the trace
from this result, we do the following integral
#
"
Z
Z
1 h

d ( (ih) + h )
(161)
Z = d (t)d(t) exp
h 0
with the understanding that (t), (t) follow the anti-periodic boundary conditions.

32