# Method of stationary phase

Many useful results in mathematical physics are the result of asymptotic approximations.
That is, the are formulae and results that are obtained by deriving an asymptotic series and keeping the leading order term as the approximation, and using the next order
term as an error estimate.^{[1]} ^{[2]}^{[3]}

On such result is called **The method of stationary phase,** which applies to integrals that resemble the Fourier transform, but have
a more general phase function. The method of stationary phase applies to one dimensional Fourier-like integrals as well as to multi-dimensional
Fourier integrals.

Many wave phenomena are the result of the preferential constructive interference of the wavefield. These phenomena are examples of result that would be described by extrema of the phase function of a Fourier representation.

## The method of stationary phase in 1-dimension

We consider integrals of the form

.

Here, is a *large parameter,* which may be frequency or wave number in problems of wave propagation.
The function is called the *amplitude* and the real-valued function is called the *phase.* As the name of the
technique implies our points of interest are places where the phase function is slowly varying.

### Simple critical point

The endpoints of integration, places where the derivatives of fail to be continuous, and places where the
derivatives of vanish are called *critical points*.

A *simple critical point* is a point , where but . Such a simple critical point is also called a *stationary point* because
this is the place where the phase function has a minimum or a maximum and is thus, stationary.

Through a number of coordinate transformations, the Fourier-like integral can be repetitively integrated by parts to yield an asymptotic series with large parameter . The leading order term of that series is called the *stationary phase formula* and the second term or third term is the estimate of
the asymptotic error of the approximation.

We consider three possibilities, that the stationary point is on the lower endpoint of integration , the upper endpoint of integration , or an interior stationary point , where .

#### Stationary phase formula -- stationary point at the lower endpoint of integration

In practice the stationary phase formula for a stationary point at the lower limit of integration is the leading order term in inverse powers of

,

as .

#### Stationary phase formula -- upper endpoint of integration

In practice the stationary phase formula for a stationary point at the upper limit of integration is the leading order term in inverse powers of

,

as .

#### Stationary phase formula -- interior stationary point , where

The stationary phase expansion at an interior stationary point is found by combining the upper and lower limit stationary phase series expansions, which causes the terms of order to cancel yielding the familiar result

,

as .

When applying these formulas, it is important to remember that this is not a simple identity, but rather a process of mathematical analysis. The procedure for applying the method of stationary phase is to

- Identify the large parameter .
- Identify the phase function .
- Find the stationary point(s) (there may be more than one) by setting .

## What about the endpoint contribution?

When we learned to compute definite integrals in undergraduate calculus, we were evaluating the integral
at the limits of integration. We call these limits of integration the *endpoints*. Thus, we are accustomed
to evaluating the *endpoint contribution*. What is asymptotic order of the endpoint contribution compared with the stationarity contribution?

We can answer that question by performing repetitive integration by parts formally to yield the following series representation. Here, we assume that there are no singularities of the amplitude or the phase of our Fourier-like integral

and apply integration by parts repetitively so as to bring down factors of the large parameter in the denominators of the resulting terms. To do this we multiply and divide the integrant by . The first application of integration by parts (integrating the exponential) yields

.

Applying this operation repetitively yields

as .

This formal result assumes that all of the parts are sufficiently differentiable, and there are no divisions by zero. The first term of the summation is as , whereas a stationary point is of an asymptotically more slowly decaying contribution of .

Endpoint contributions occur wherever there are discontinuities in the data. These constitute such phenomena
as ringing, or *diffraction smiles* seen in processed seismic data.

## Higher order stationary points in 1D

Similar formulas may be derived for higher order stationary points. In general, such a stationary point would be represented as critical points in the amplitude and in the phase .

### Higher order stationary point at the lower endpoint of integration at .

In general, we may consider the amplitude factor to represented by the more general power series representations

for

and the phase is represented as

for .

Here and are the first non vanishing coefficients of the power series representations of the functions.

The resulting stationary phase formula is

as .

### Higher order stationary point at the upper endpoint of integration at .

In general, we may consider the amplitude factor to represented by the more general power series representations

for

and the phase is represented as

for .

The resulting stationary phase formula is

as .

Here, the Gamma function is

## Stationary phase analysis

A common misunderstanding of the method of stationary phase is that this merely a formula lookup. A better approach is to consider that applying the method of stationary phase is a method of analysis. To perform this analysis, given a Fourier-like integral, the following steps must be applied

- identify the large parameter
- identify the phase function, it's derivative, and next highest order non-vanishing derivative at the zeros of the first derivative of the phase. If this is the second derivative of the phase, then this is a simple critical point, also known as a
*stationary point*. - find the stationary point(s)
- apply the appropriate stationary phase formula.

### Example, asymptotic form of the Bessel function for large

We consider the integral representation of a Bessel function

- is the large parameter. Here, both and
- , , and
- the critical points are simple interior stationary points at
- we apply the stationary phase formula for simple interior stationary points to each stationary point and combine the results

.

Simplifying, yields the common asymptotic form of the Bessel function for large

.

Thus, at the end of this analysis, we have the asymptotic approximation as well a an order estimate of the error. In this case, the product might further be reduced into a large approximation, or a large approximation, where either the or the , respectively, become part of the phase.

# Multidimensional stationary phase

For problems involving Fourier-like integrals in higher dimensions

.

as . Here, is a subdomain in dimensions, and is an dimensional hypervolume (or hypersurface element).

### The Multidimensional stationary phase formula

For an interior stationary point, the asymptotic representation of the multidimensional Fourier-like integral, for a simple stationary point, which is to say a point where but for which the determinant of the Hessian matrix . Here . As with other discussions, the large parameter

as . Here is the *signature* of
the matrix , which is the number of positive eigenvalues of minus the
number of its negative eigenvalues.