Round 3! Fight… with path integrals! XD

## Introduction: generic aspects

Consider a particle moving in one dimension (the extension to ND is trivial), the hamiltonian being of the usual form:

The fundamental question and problem in the path integral (PI) approach of quantum mechanics (QM) is:

If the particle is at a position q (fluctuating, of course) at a time , what is the probability amplitude that it will be at some other position, say , at a time ?

Any other “amplitude” can be written in terms of the previous and fundamental amplitude above. It is easy to get a formal expression for this amplitude in terms of the usual formulation of QM simply solving the Schrödinger equation. Furthermore, you can make a guess using the evolution operator formalism as well! Dirac did it himself, but we had to wait for Feynman to “clean” and understand the formalism. Let me introduce the eigenstates of the position operator . It forms what physicists and mathematicians call “complete, orthonormal set”. Without entering into the technical (and axiomatic) issues, we have:

Then, the initial state is

Letting the state evolve in time and projecting on the state , we have the amplitude

when there is the possibility of an ambiguity, operators will be denoted with a “hat”. Otherwise, we will drop that symbol. We also fix the rationalized Planck constant to the unit, i.e., we write .

The object is called the **propagator (sometimes Green function) **from the initial spacetime point to the final point . The propagator is clearly independent of the origin of time, for time independent hamiltonians:

We will derive an expression for this amplitude in the form of a sum over all possible path between the initial and final points. For this purpose, let us separate the time evolution in the above amplitude into two smaller time evolutions, with the aid of the exponential function:

Thus, the amplitude becomes

Inserting a factor “1” in the form of the previously seen integrals for an orthonormal set, we get the next sum over the position eigenstates:

This formula is just an expression of the quantum mechanical rule for combining amplitudes: if a process can occur a number of ways, the amplitudes for each of these will add. And since amplitude is a complex number (!), therefore interference (constructive and destructive) will arise at some points in general. In particular, any particle, in its propagation from to , must be somewhere at an intermediate time . However, this “somewhere” is a bit “fuzzy”, since it is spread along every possible QM path, not only the classical trajectories are possible. Of course, the weights of paths are not equal in general, but, in principle, every path does contribute! Labeling an intermediate position with , we can compute the amplitude for the propagation via the point , and it is the product of the two propagators as it is shown above. In general, we have always integrated over all possible intermediate (generally non-observable and virtual) positions. This seems a bit of mathematical trickery but it is all true. Quantum particles “smell” every possible path. We can repeat the division of the time interval ; let me divide it up into a large number, say N, of time intervals of tiny duration . Then, we can write for the propagator the next formula:

We can again insert some complete sets of states between every exponential function, and it yields

We can rewrite this expression as follows:

or, equivalently,

with

and

and where we have defined and . Note that these initial and final positions are NOT integrated over, but only the intermediate states. This expression says that the amplitude IS the integral of the amplitude of all N-leg paths. Apart from some additional mathematical tools and details concerning the limit when , this is clearly going to become a sum over every possible path of the amplitude for each path. Namely, we have

where at formal level, we identify

and

A final and important result is the so-called Lie-Kato-Trotter formula. For any pair of matrices/operators A,B, under very general conditions and assumptions, we have the identity:

Let us make use of the Lie-Kato-Trotter formula to derive the so-called Dirac formula. For the hamiltonian

with (mass equal to one) and

where

The amplitude reads

We insert the set

and thus

Carrying out the operations, we get

or

And now, using

we finally obtain

namely,

i.e., the modern language version of the Dirac formula in terms of the classical action reads as follows

Now we are ready to writhe the propagator and the PI in the configuration space. Let me write

or

We apply now the above Dirac formula:

and where

The argument inside of the exponential function is a discrete approximation of the action of a path passing through the points . Moreover, in the continuum limit, we obtain the formal definition of a Feynman PI as follows:

and where and are the initial and the final points in the amplitude. This final result is the configuration space path integral. It should be viewed as a notation of an infinite dimensional integral, or infinite differential form, of the more precise discrete expression as . Furthermore, this formula says that any particle, in going from one position to another, takes ALL possible path between these two points. It can be shown that, while all paths contribute, the classical path is the dominant one and it minimizes (or maximizes) the path integral. A key observation for these facts is to consider how different and close paths interfere to each other.

We can perform all this formalism in the phase space as well! Take

as the propagator for one subinterval. We can expand the exponential function for small

It yields

The first term is a delta function, and we write

The second term can be worked out as well inserting a complete set of momentum eigenstates between H and . This gives:

or equivalently

where we used that . We must view the momentum in the hamiltonian as the momentum operator operating to the right, while the potential (operator valued!) should operate to the left. Our last expression is asymmetric between and , but the origin of this is our choice of putting the factor “1” to the right of the hamiltonian H. Had we write it to the left instead to the right of the hamiltonian, we would obtain the potential instead . If we introduce the shift

and thus , the expression would be symmetric in this variable. The exact choice does NOT matter in the continuum limit of a path integral! Now, we calculate

and

There, we have N such factors in the amplitude. We can combine them, and write

and thus

and where we have neglected a multiplicative factor

which will tend to one in the continuum limit . Then, the propagator finally becomes

**Remark:** there is one momentum integral for every interval (N in total), while there is one position integral for every intermediate position (N-1 in total). In the continuum limit, , this last expression approximates an integral over functions , and we write it as

where

is the classical action in phase space! This phase integral is over all the classical functions , with such as and .

Now we are ready for two simple and frequently found examples of path integrals, very important indeed in practical applications!

## Examples of path integrals

The first easy example is the free particle propagator and its path integral. The propagator will be related to the hamiltonian

We can compute the propagator with simple QM:

thus

so

The integral is a gaussian, so we can integrate it inmediately, to obtain

and where the classical action for a given time reads

Of course, this approach is the normal QM procedure. We can get the same result using our brand new path integral approach!!! The configuration space PI reads:

or expanding the last sum

where the initial and final points are and . The integrals are all gaussian and can be evaluated exactly. The only issue is that they are “coupled” and it complicates the integration. You can do it as an interesting exercise. The final result is

and taking into account that the total time is we write and check the previously obtained result

with as the classical action as before. If we restore the Planck constant, the free particle propagator reads

and the propagator separates into two factors, one of which is the “classical phase” and the remaining one is a weird term depending on the square root of the mass and the time T.

The second example is the calculation of the path integral and propagator for the harmonic oscillator potential. We have:

The harmonic oscillator action reads

The paths are integrated out with and . To perform the integral, suppose we do know the solution of the classical problem, . That is:

so

and we change variables in the path integral to

Integrating out over all deviations from the classical path is equivalent to integrating out all the possible paths.

Remark: The jacobian of the above transformation is 1.

Remark(II): Since the classical solution obeys the correct boundary conditions, the paths y(t) over which we integrate go from

to

The action can be written in terms of the path as a power series in :

The terms linear in vanish by construction. Then, we write

By direct substitution, we get

The path integral separates into two pieces:

1st. The term of the action with the classical action.

2nd. The path integral over the deviations from the classical path. This second factor is independent of the initial and final points!

We can calculate the path integral in position space, as we did above for the harmonic oscillator or, alternatively, we can compute it in the Fourier space, writing

and integrating over the coefficients , the result will be

Remark: For tiny frequencies, the harmonic oscillator propagator becomes the free particle propagator.

Remark(II): The classical action can be evaluated in a straightforward fashion, to get

A final comment follows: the path integral for ANY quadratic action can be evaluated EXACTLY, essentially performing a path integral with the aid of gaussian integrals!

## Statistical Mechanics with Path Integrals

The central object in Statistical Mechanics is the partition function Z. Indeed, there is a beautiful correspondence between Statistical Mechanics and the path integral approach to QM by realizing the link between the partition function and the path integral (or the propagator). The classical partition function reads

and where is the reciprocal temperature (in natural units) and are the energies of the jth levels of the system, . We can write,

(1)

But, the propagator definition was

Suppose we consider T to be a COMPLEX parameter (it sounds a bit weird since time is real, not complex) and consider it to be a pure imaginary number. In particular, we write a “Wick” rotation of time as the temperature!

where t is the absolute temperature (be aware, here t is NOT time but temperature and, though, they are closely related) and it is a REAL number. Therefore,

Inserting and integrating over the position variable q, we obtain

And thus, we have the formal identity and link between the path integral (the propagator) and the partition function of statistical mechanics:

Put it into (equivalent) words:

**1st. The propagator evaluated at negative imaginary time (a Wick rotated time) IS the partition function. **

2nd. Temperature is a Wick rotation of time or time is a Wick rotated temperature.

3rd. The propagator in diagonal position representation IS the partition function.

Fascinating! We can easily work out an elementary example with the harmonic oscillator above. Recall the PI for it:

Write the Wick rotation and (diagonal representation of the partition function/propagator):

This propagator gives a partition function

and then

But this is precisely the classical (quantum) result

and where we reinserted the Planck constant in the last step! Path integrals are disguised partition functions and vice versa!

We can rewrite the partition function in terms of a path integral. In ordinary (real) time,

and where the integral is taken over all paths from to . In the diagonal representation, and is the Wick rotated time/temperature, so

and the integral in now along the negative imaginary time/temperature axis! Let us define a real variable for this integration, say . Here, is called the imaginary time, since when the time is IMAGINARY, then is REAL. Then, the integral over the IMAGINARY time is along the real axis : when runs from 0 to , then runs from 0 to . We can write q as a function of the variable : . It defines the so-called imaginary time formalism. Then,

Note the similarity of the above last equation to the time dependent Schrödinger equation and the classical Hamiltonian Dynamics! Then, the propagator becomes

The integral is taken over all function such that

This final result is a Wick rotated “imaginary time” or “euclidean” path integral, defined by associating to any path an amplitude (statistical weight) , where is the so-called euclidean action, obtained from the usual (or Minkovskian, spacetime) action by changing the sign of the potential energy term.

In summary, we have got two main formulae:

1st. The Minkovski or real time propagator, in which we have an oscillating integral of the classical action

2nd. The Euclidean or imaginary time propagator, in which we have an exponential integral of the classical action

The integral here is over all functions such as

We can use one of another formula, but usually euclidean time are preferred in most of the applications due to convergence features. We could perform all the calculations in a Euclidean theory and analytically continue the results to the real time when computing the physical quantities as well!

From the partition function, the statistical partition function, it is possible to extract the ground state energy. In fact, from the definition of Z

and

We can see that the contribution of every state decreases exponentially with , the temperature (reciprocal of). However, the contribution of the energy of the ground state decreases less slowly than any other state contribution. Thus, in the limit of large , i.e. low temperature, the ground state contribution will dominate. We arrive at

We can do the same trickery in an alternative way. Let us look at the euclidean time propagator from to :

We can insert a complete set of eigenstates of H:

or

and where are the wave functions of H. Again, the ground state dominates at low temperatures (), and we get

Namely, in the low temperature regime the difference between the terms and goes to zero. And this result is important…

## Density matrix

The density matrix is an important object. It has a well known expression in quantum statistical mechanics:

Using the previous results, we have

Indeed, every observable can be related to the partition function, the path integral, and the density matrix with these definitions! The expectation value for some operator is given by

Using the path integral, we can write

## Correlators in QM and QFT

The position operator in the usual Heisenberg representation is defined in terms of the Schrödinger operator picture by means of the identity

The time eigenstates read

and its relation with the time-independent eigenstates is

We can express the path integral with these states:

We can define the time ordering of operators as follows. Let us have two time dependent operators A(t), B(t). The time ordered operator (TOO) is:

We also define the n-order Green function, or n-point correlator function as the following expectation value in QM:

An important case is the so-called 2-point (correlation) function. The two point function has a very special importance in QM. Let us now calculate the 2-point function via the path integral formalism:

Firstly, we will calculate

and then, we will take the limits and . Suppose that , then

Every matrix element is itself a path integral:

This last expression have 3 different pieces:

1st. A path integral form the initial position to an arbitrary position .

2nd. A path integral from to a second arbitrary position .

3rd. A path integral from an arbitrary position to the final position at .

We are integrating over all paths from q to q’, subject to the restriction that the paths pass through the intermediate points at and . We then integrate over the two arbitrary positions, so that in fact we are integrating over ALL paths: we can combine these three path integrals plus the integrations over and into one single path integral. The factos in the above integral can be incorporated into this simple integral by including a factor into the path integral. Therefore, we read

Mimicking this calculation, in the case of :

That is, the same expression is valid in both cases! The path integral does the time ordering automatically!!!! Taking now the limit , we get

For the n-point correlation functions, we obtain similar results:

## The generating functional Z(J)

In order to define a generating function Z(J), we have to build up some preliminary definitions. The general amplitude

The vacuum amplitude or 0-0 vacuum persistence amplitude is:

with

where the integration is calculated over all the closed paths such that

We will see how this amplitude is related to the ground state of survival (or persistence) amplitude. Firstly, we add a new “source” to the action of an arbitrary system. Thus, we modify the action with an extra term/piece:

and where is an arbitrary function such as

The generating functional Z(J) is defined now to be the quantity

with as before. The functional derivatives of the generating functional Z(J) are very interesting objects:

evaluated at provides the one-point correlation function! The proof is simple:

or equivalently

Iterating this process, we obtain a path integral with several q’s in the numerator, i.e., we get the n-point (correlation) function!

Indeed, using the functional derivatives and expanding the generating functional Z(J) in Taylor-like series with these derivatives in a series on J, we get the formal series:

where

and so on! You can understand now why Z(J) is called the generating functional: its expansion provides the n-point correlation functions!

We can compute the generating functional of the harmonic oscillator as a working example. Suppose that initially the action of this system reads

Then, the corresponding denominator, N(0), is the propagator K we computed for the harmonic oscillator before! For the numerator in the generating functional, we write

By definition, we also get

Now, we use the same mathematical magic we did for the unforced harmonic oscillator, defining

and we integrate over all the possible paths and the classical solution satisfies

and then

The linear terms in y vanish since satisfies the equation of motion. Therefore, the path integral reduces to

And now, the path integral over y is a constant, independent from J, that we will call C. We write

where

We can compute this last integral with the aid of the equation of motion and integrating the first term by parts:

and it yields

We can also determine the classical path and define the Green function with this new formalism. The classical equation of motion is

The classical path can be written in terms of the Green functions defined by the formal identities

and

We can solve these equations for the Green function by going into the momentum space representation by Fourier transformations. The result is

However, there are poles on the axis of integration. The Green function is ambiguous unless we give it a pole prescription, i.e., a boundary condition (namely, we have to avoid the pole singularity). We require that G go to zero as . The correct pole prescription then turns out to be

Note that in the case of equal times, the integral is FINITE and the Green function simplifies to

We can now write

Dividing by the denominator merely cancels the constant C factor, and thus

and the final result will be

and where G is the Green function defined above!

Now we can study some general features or properties. Let me take any functional K(J):

We can compute the first functional derivative for any time

The second derivative can be also calculated further

Indeed, we can generalize this result to an arbitrary functional F in the following way

This expression can be proved correct in two steps:

1st. Put inside the PI, since it does not depend on the path.

2nd. Expand F in a functional Taylor series, and every functional derivative with respect to J provides a q term in front of the exponential term.

Remark: As an example, take

We get in this case

The previous expression is the starting point of perturbation theory in QM! Of course, we are making perturbation theory with the path integral formalism, not in the usual way but it looks pretty similar and it is done in a completely covariant way. That is the hidden power of path integrals. Let us expand the exponential in the right-handed side above and write the expression as follows

or

or

This is the analogue form of the so-called Dyson series in the PI formalism! For small values of the parameters , we get an approximation to the PI truncating the series at any desired order. At first order in we write

If we plug the condition J=0 at the end of our calculation, it yields

A very important popular (VIP) example is provided by the quartic interaction potential. It provides a functional term

The potential term is very important in the Standard Model and Quantum Mechanics. We obtain

and where

Have you enjoyed this post? See you in my last introductory path integral TSOR post…

P.S.: A special post, the 150th, is coming soon! Stay tuned!

First, nice blog!

Second, I’ve watched a conference held by Prof. Unruh a few months ago, where he claimed the experimental proof of his effect in the vibration field on the surface of a flowing liquid, i.e. an “acoustic black hole”. He argued that his and the Hawking effect are in fact one and the same, and that there’s no real “quantum” origin behind the black hole thermal radiation. He also dismissed the whole “virtual particle pair creation at the event horizon” as something that he considers unphysical nonsense. I don’t remember the details though as GR and quantum gravity are not my field, I guess his papers probably explain it better. What do you think about this? In a sense, considering this to be a natural effect that takes place in all fields without any intervention from quantum effects, and that it arises in acoustic fields as well as in extreme conditions like black holes would be a reassuringly elegant solution to the issues at hand. In another, it seems weird that something so trivial as waves on a liquid surface can arise from the same reasons that cause something as crazy as black holes to irradiate (and potentially evaporate), when we know that at the event horizon there might as well be effects we don’t understand at play, given how we can’t reconcile GR and quantum mechanics.

BTW, the paper by Unruh describing the experiment can be found here:

http://arxiv.org/abs/1008.1911

Well, as I have told you in facebook. I am skeptical about what Unruh is defending now. Who knows? Who nose? Unless experimental evidence of Unruh effect and other quantum gravitational experiments involving the vacuum appears, we have only theories and conjectures…:) Best wishes!