Finding the Prefactor in the Simple Harmonic Oscillator Propagator
Michael Fowler, UVa
Recall that the free particle propagator has the form
From a classical mechanical evaluation of the action, we can show the simple harmonic oscillator propagator has the form
where is a so far unknown function of time. For sufficiently small though, much less than the period of the oscillator, the potential will have no significant effect, so the propagator must tend to the free particle propagator. That is, for tending to zero, But we also know that at time all simple harmonic wavefunctions return to their values, so a particle localized at will be localized at again after one period. In the exponent we do indeed see this cyclic behavior. In the prefactor, we must replace by that is,
Note this expression is also consistent with the free particle propagator in the limit going to zero, in other words, the vanishing of the simple harmonic oscillator potential.
Approximating Integrals by Stationary Phase Techniques
Actually, this result can be derived from the integral over the fluctuations about the classical path. The argument is closely analogous to that for the free particle, and the following equation is a straightforward generalization of that case (discussed in the previous lecture):
Just as for the free particle, the middle term in that last line is zero, as it has to be since it’s the linear term in the path deviation from the stationary point. As before, the path (the deviation from the classical least action path) is zero at the two endpoints, so we can integrate by parts to get
identically zero from the classical equation of motion (after all, this is how it’s derived from least action!)
The hard part is evaluating the integral over paths in the third term. To see how to do this,
it’s worth briefly reviewing a stationary phase integral in ordinary space, that is,
where we assume that the real function has an absolute minimum at the point and take variables Then to leading order near the stationary point, Taking only this term, the integral becomes
The matrix is positive definite and symmetric, it has real orthogonal eigenvectors, and we can choose our coordinate axes in the directions of those eigenvectors. In this case, is a diagonal matrix, and its diagonal elements are just its eigenvalues
Recall that is invariant under an orthogonal transformation, and so is unaffected by transforming back from the eigenvector basis to the original coordinates.
Applying Stationary Phase to the Integral Over Paths
We’re now ready to generalize this result from an integral over an dimensional space of ’s to an integral over the infinite dimensional space of paths with boundary conditions and Note that for the particular case at hand, the Simple Harmonic Oscillator, the leading order quadratic term is in fact the whole story, and so gives the exact result. This is not, of course, true for general quantum systems.
As a preliminary step, we integrate by parts to write
The path integral then becomes:
Comparing this with the finite-dimensional version, we see that the component vector is replaced by the continuous function and the matrix by a differential operator acting on the space of path functions equal to zero at and
The operator is well-defined and Hermitian, with eigenstates and eigenvalues:
Any reasonable path satisfying the boundary conditions can be written as a sum over these eigenstates with Fourier coefficients
The functional integral over all paths now becomes an infinite product of ordinary integrals over the variables This is precisely analogous to the finite integral above, except that there is no maximum value The operator is diagonal with respect to this base formed of its eigenstates, and therefore can be represented as an infinite diagonal matrix with the diagonal elements equal to the eigenvaluesand, just as for the finite matrix case, the gaussian integrals can be carried out in succession, each one giving a term
When we take the infinite product of all these terms, we find in the denominator the square root of the product of all the eigenvaluesand that product is the determinant of the operator. That is to say,
In the last step, we used the infinite-product representation for the sine function:
The rather disconcerting factor which does not depend on and hence has nothing to do with the dynamics, can be thrown in with similar time-dependent and constant factors involved in the measure of the integral over paths, and their overall contribution can be nailed down by the simple observation that for we must recover the known free-particle propagator. So, we can be quite careless about overall multiplying factors!
Therefore, matching to the free-particle propagator in the limit we find for the simple harmonic oscillator.
I used an excellent book, Solitons and Instantons, by R. Rajaraman, for this section. The original work was all by Feynman.