Laplace's method and steepest descent
Extract a Gaussian contribution near a nondegenerate minimum and track the analytic requirements behind a complex contour deformation.
On this page
The question for today
What you’ll understand
- Derive the leading Gaussian scale at an interior minimum.
- State nondegeneracy and localization assumptions for a Laplace expansion.
- Justify a simple steepest-descent contour shift without crossing singularities.
A large positive parameter in a decaying exponential concentrates an integral near the lowest value of its phase. Laplace's method replaces the phase near an isolated nondegenerate minimum by a quadratic and rescales the narrow contribution to a Gaussian. Steepest descent extends the geometry to complex contours, where the path is deformed so that the real part of the exponent decays away from a saddle. Both methods depend on global information: an overlooked endpoint, singularity, or competing saddle can invalidate a locally correct calculation.
A quadratic minimum selects an inverse-square-root width
For a clear real-variable theorem, take a smooth compactly supported amplitude on the real line and a real smooth phase near its support. Assume the phase has a unique global minimum on that support at an interior point where the amplitude is nonzero, its first derivative vanishes, and its second derivative is strictly positive. Compactness and uniqueness separate the phase from its minimum outside any small neighborhood. Thus the distant contribution is exponentially suppressed, while the nearby quadratic expansion determines the leading term.
Set the distance from the minimum equal to a new variable divided by the square root of the large parameter. The quadratic phase becomes order one and the integration measure contributes the inverse square root. Smooth expansion of the amplitude and higher phase terms then produces corrections. Under the stated smooth compact-support assumptions, the leading Gaussian term has an absolute remainder of order the parameter to minus three-halves times the exponential minimum factor. Boundary minima or degenerate critical points require different powers and constants, so check the geometry before using the formula.
A Gaussian polynomial provides an exact check
For a quadratic phase and a polynomial amplitude, some of the asymptotic calculation can be exact. The integral of a Gaussian is the square root of pi divided by the parameter. Differentiating that identity with respect to a positive parameter, justified by Gaussian domination on compact positive parameter ranges, gives its second-moment integral. Adding those two contributions evaluates an amplitude equal to one plus the squared coordinate. This is an improper-integral example with explicit tail control rather than an unproved extension of the compact-support theorem.
The correction is smaller than the leading term by the reciprocal of twice the large parameter. Relative error must specify its denominator: the omitted correction divided by the leading approximation differs from the omitted correction divided by the exact positive integral. Exact examples are valuable for checking these distinctions and normalizing constants. They do not establish that a leading-term approximation has the same finite-parameter error for every amplitude or phase with a similar-looking minimum.
Complex saddles require permitted contour geometry
For a holomorphic amplitude and phase, a steepest-descent contour through a nondegenerate saddle keeps the phase's imaginary part constant locally and makes its real part increase away from the saddle for an integrand with negative exponential phase. Cauchy's theorem permits a deformation only through a region where the integrand is holomorphic, with endpoint and infinite-tail contributions controlled. Crossing a pole adds residue information; crossing a branch cut changes the continuation problem. Locating a saddle alone does not establish that its steepest path belongs to the original integral.
The orientation of the local contour fixes the phase of the Gaussian square root. Several accessible saddles may contribute, and their relative importance can change as parameters move. A coalescing saddle or a saddle approaching an endpoint needs a uniform canonical approximation rather than separate nondegenerate formulas with exploding coefficients. The contour-shift practice below is deliberately entire and quadratic, so the connecting vertical edges can be bounded explicitly. That gives a complete deformation argument without concealing singularities or unverified tail behavior.
Work it through
Quantify a leading Laplace approximation
Evaluate , compare it with the leading Gaussian term, and calculate the relative error using the exact integral as denominator.
The phase is , with minimum at zero and second derivative two. The leading term is .
The exact second-moment correction multiplies this by , giving .
The omitted amount is the leading term divided by fifty. Dividing instead by the exact integral gives relative error .
The leading term underestimates this exact integral by about 1.96%. The two-term expression is exact for this chosen amplitude and quadratic phase.
The living notebook
Explore the model
Interactive lab · function
A quadratic phase controls concentration width
Reset and set both . Predict phase , so gives exponent size one there. Reset and set both ; the same phase threshold occurs at . The instrument plots the phase , not the exponential weight or its integral. Computing from the height shows why quadrupling the large parameter halves the dominant width.
- y at chosen x
- 0.608
- Tangent slope · 2ax
- 1.35
- Signed area · ax³/3
- 0.182
This lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.
Look forNear a nondegenerate minimum, the relevant region narrows like the inverse square root of the large parameter.
Pause & reason
Check your understanding
Your turn
Try it yourself
For real and , evaluate by shifting to the horizontal line through its saddle. Explain why the connecting edges vanish.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
The integrand is entire. On a vertical edge with real part ±R and bounded imaginary part between zero and b, its modulus has a factor e to minus lambda times R squared.
Compare with the solution
Shift the contour to , where the exponent becomes . On each connecting vertical edge, , so its bounded-length integral vanishes as for fixed . Cauchy's theorem gives the exact result . No poles or branch cuts are crossed.
Keep the idea
Take it with you
- A nondegenerate minimum gives a Gaussian scale and its normalization.
- Error order and relative-error denominators must be stated precisely.
- Steepest descent requires an allowed analytic contour deformation and controlled tails.
Read further
Original sources and deeper explanations to follow at your own pace.
Optional progress, saved in this browser when available. Revisit any lesson at any time.
Lesson not marked complete.