Michael Fowler, UVa
This is an application of the Saddlepoint method of evaluating integrals discussed earlier.
The WKB semiclassical approximate solution to Schrödinger’s equation,
is reliable in regions where the wavelength (for oscillating solutions) or the decay length (for exponential solutions) changes only slightly over a distance of one wavelength or decay length respectively. For a particle trapped in a (one-dimensional) potential well, classically the particle would bounce back and forth between the two turning points where its kinetic energy vanishes. In the quantum case, these are precisely the points where the wavelength becomes infinite, so the WKB solution fails.
However, for a reasonably smooth potential it may be an adequate approximation to treat a turning point region as one where the potential is increasing linearly with distance over a sufficient range that beyond this point the WKB approximation can be used in both directions.
The solution of Schrödinger’s equation for a linearly increasing or decreasing potential is well known, it is the Airy function, the solution of the differential equation
plotted here at the left-hand turning point :
The strategy is to evaluate this function for large x, both positive and negative, so that we can join together the two WKB solutions, valid in the far regions, in a quantitative fashion.
Following Mathews and Walker (page 116) the differential equation is most simply solved by taking its Fourier transform. If
This is an exact result, but a nontrivial integral! Fortunately, we are only interested in its value for large values of |x|, and this is precisely where saddlepoint methods become accurate: the contour is distorted to be as low-lying as possible, then the only significant contributions to the integral (thanks to the exponential variation of the amplitude) are from those places where we must go over a saddlepoint to get from one valley to the next. To see where the saddlepoints are, and how they relate to an integral along the real axis, we plot below contour maps of the absolute value of the term in the exponential.
Large positive x: In the map below, we take x = 10 (rather small), so the saddlepoints are at . If the path of integration is moved down from the real axis into the (dark) valleys (so the integrand becomes exponentially smaller) the contour from -∞ will come up from the bottom left to the left saddlepoint, over the saddle into the top center valley, then back over the second saddlepoint into the bottom right valley and on to +∞.
Lighter color means higher ground.
Writing the integral as
(dropping irrelevant overall constants) then
(Mathews and Walker take the “large” parameter x out of f, we’ve left it in—this doesn’t affect the final result.) Near the positive saddlepoint
dropping higher-order terms. Since f '' is pure imaginary at the saddlepoint, the appropriate path for a real exponent in the Gaussian integral is at π/4 to the x-axis. So in the path integral
where ds is a real parameter measuring incremental path length, and the sign in the exponent is positive for the saddlepoint on the left. The contributions from the two saddlepoints give the asymptotic (large positive x) solutions as:
Large negative x: In this case, the saddlepoint geography is quite different, although the distant geography is the same, being dominated by the w 3 term.
Just as before, the real axis integration path can be moved down from the real axis into the valleys at bottom left and bottom right. (The extra contributions from linking up the new path with the real axis at infinity are zero.) It is clear from the map above that to get from the valley on the left to the one on the right means just going over the saddlepoint on the negative imaginary axis. Note from the darker shading that the other saddlepoint is at higher elevation.
The integration through the saddlepoint is parallel to the real axis, and gives
This, then, is the decaying wavefunction solution of the Airy equation that we are looking for, and it is clear that it goes smoothly from this exponential to the cosine form
as x is taken along the real axis from large negative to large positive values.