Fact-checked by Grok 2 weeks ago

Stationary phase approximation

The stationary phase approximation is a mathematical technique in asymptotic analysis for evaluating oscillatory integrals of the form \int f(x) e^{i \omega \phi(x)} \, dx as the large parameter \omega \to \infty, where the dominant contributions to the integral arise from the neighborhoods of stationary points x_0 at which the derivative of the phase function vanishes, \phi'(x_0) = 0.^[1] This method exploits the rapid oscillations of the integrand away from these points, leading to destructive interference and negligible contributions elsewhere, thereby simplifying the integral to a local Gaussian-like approximation around each stationary point.^[2] Originating in the 19th century, the method was first sketched by William Thomson (later Lord Kelvin) in 1887 and further developed by George Gabriel Stokes, who applied it to problems in wave theory and diffraction.^[3] For a non-degenerate stationary point where \phi''(x_0) \neq 0, the leading-order asymptotic formula is \sqrt{\frac{2\pi}{\omega |\phi''(x_0)|}} f(x_0) e^{i \omega \phi(x_0) + i \frac{\pi}{4} \operatorname{sgn}(\phi''(x_0))}, providing an explicit approximation that captures the amplitude and phase of the integral.^[1] Higher-order terms can be derived by expanding the phase and amplitude functions in Taylor series, yielding a full asymptotic series for improved accuracy.^[4] The approximation extends naturally to multiple dimensions, where contributions localize to critical points satisfying \nabla \phi(x_0) = 0 with non-vanishing Hessian determinant, and the leading term involves the determinant of the Hessian matrix in the denominator.^[4] It finds widespread applications in physics and engineering, including wave propagation in optics and acoustics, quantum mechanical semiclassical approximations, and the asymptotic evaluation of special functions such as Bessel and Airy functions.^[2] For instance, the large-argument behavior of the Bessel function J_0(x) is given by \sqrt{\frac{2}{\pi x}} \cos\left(x - \frac{\pi}{4}\right), directly obtained via this method.^[2] Related techniques, such as the method of steepest descent for complex phases, build upon these principles for broader oscillatory and exponential integrals.^[1]

Introduction

Definition and motivation

The stationary phase approximation is a fundamental technique in asymptotic analysis designed to evaluate integrals with highly oscillatory integrands in the limit as a parameter becomes large. It targets integrals consisting of a slowly varying amplitude function multiplied by a rapidly oscillating phase factor, where the oscillations intensify with the parameter's magnitude. The core idea is that the integral's value is predominantly determined by contributions from neighborhoods around stationary points of the phase function, defined as locations where the phase's derivative vanishes.^[1]^[5] This method is motivated by the practical difficulties in computing such oscillatory integrals directly, as the swift phase variations cause substantial cancellation between positive and negative components of the integrand, rendering numerical integration inefficient and prone to errors. By isolating the stationary points, the approximation efficiently captures the leading-order behavior, transforming a potentially intractable problem into a manageable one focused on local expansions near these critical regions. This utility arises in contexts like wave propagation and spectral analysis, where large parameters naturally emerge and demand asymptotic simplification.^[6]^[1] In the framework of asymptotic methods, the stationary phase approximation complements tools for non-oscillatory cases by addressing the unique challenges of complex exponentials. It leverages destructive interference in regions distant from stationary points, where rapid oscillations ensure near-complete cancellation and contribute negligibly to the integral. Consequently, the approximation's error exhibits exponential decay outside these dominant zones, providing robust estimates whose accuracy improves markedly with increasing parameter size.^[5]^[6]

Historical development

The stationary phase approximation traces its origins to 19th-century investigations into wave phenomena in physics. George Gabriel Stokes laid early groundwork in the 1840s and 1850s through his studies of oscillatory integrals in optics and fluid dynamics, particularly in analyzing interference patterns and the asymptotic behavior of wave functions where phase variations are minimal at certain points.^[7]^[8] His work on Airy's integral in the context of rainbow formation implicitly utilized principles akin to stationary phase to identify dominant contributions from non-oscillating regions.^[7]^[8] Lord Kelvin (William Thomson) provided a more explicit formulation in 1887, introducing the principle of stationary phase in his analysis of ship waves. In his paper "On Ship Waves," Kelvin applied the method to approximate oscillatory integrals describing the far-field pattern of water waves generated by a moving vessel, emphasizing that the main contributions arise where the phase function is stationary, i.e., its derivative vanishes.^[3] This heuristic approach, motivated by physical intuition, marked a key advancement in evaluating highly oscillatory integrals asymptotically as the oscillation frequency increases.^[9] Advancements in the late 19th century included contributions from mathematicians refining asymptotic techniques for oscillatory integrals. Henri Poincaré developed foundational ideas for asymptotic expansions in 1886, providing tools to justify divergent series and uniform approximations relevant to stationary phase applications in differential equations.^[10] By the early 20th century, the method gained broader mathematical rigor; G. H. Hardy’s 1949 treatise Divergent Series offered a systematic treatment of asymptotic methods, including stationary phase, framing it within the theory of divergent series and their summation.^[11] A significant milestone came with the inclusion of the stationary phase approximation in Richard Courant and David Hilbert’s Methods of Mathematical Physics (1953, Volume II), where it was presented as an essential technique for asymptotic solutions to partial differential equations in wave propagation.^[12] In the 1970s, Lars Hörmander integrated the method into microlocal analysis, rigorously extending it to multidimensional oscillatory integrals and the propagation of singularities in PDE solutions via complex stationary phase expansions.^[13] Initially a heuristic tool for physical problems, the approximation was fully justified mathematically in the 20th century through complex analysis, transforming it into a cornerstone of asymptotic theory.

Core Concepts

Oscillatory integrals

The stationary phase approximation is primarily applied to oscillatory integrals of the form

I(\lambda) = \int_{\Omega} g(x) \exp(i \lambda f(x)) \, dx,

where \lambda is a large positive real parameter, g: \Omega \to \mathbb{C} is a smooth amplitude function, f: \Omega \to \mathbb{R} is a smooth real-valued phase function, and \Omega is a domain such as \mathbb{R}^n or a compact manifold.^[4]^[1] This standard form captures the essential structure of integrals arising in asymptotic analysis, wave propagation, and Fourier transforms, where the exponential term introduces rapid oscillations as \lambda increases.^[14] For the approximation to be valid, several assumptions must hold: both g and f are infinitely differentiable (C^\infty) on \Omega, with g often having compact support or sufficient decay to ensure integrability.^[4]^[1] Additionally, [\lambda](/page/Lambda) \to +\infty, and the integral converges, which is typically guaranteed by the decay of g at infinity or appropriate boundary conditions on compact domains.^[14] These conditions ensure that the asymptotic behavior can be rigorously analyzed without divergences.^[4] As \lambda becomes large, the integrand oscillates rapidly except near points where the phase f has minimal variation, known as stationary points; away from these, the oscillations lead to near-cancellation of contributions, causing the integral to localize primarily around such points.^[1]^[14] The stationary phase method quantifies this localization, providing an asymptotic expansion that captures the dominant terms.^[4] The large parameter \lambda scales the frequency of oscillation in the exponential, distinguishing this framework from Laplace's method, which involves real exponentials \exp(-\lambda f(x)) for non-oscillatory decay toward minima of the phase.^[1]^[14] This oscillatory nature shifts the focus to phase stationary points rather than exponential suppression.^[4]

Stationary points

In the context of oscillatory integrals, stationary points of the phase function play a central role in determining the asymptotic behavior as the oscillation parameter becomes large. These points are defined as the critical points x_0 where the gradient vanishes, \nabla f(x_0) = 0; in the one-dimensional case, this simplifies to f'(x_0) = 0.^[15] A key assumption for the leading-order stationary phase approximation is the non-degeneracy of these critical points, meaning the Hessian matrix H = \operatorname{Hess} f(x_0) is invertible, or equivalently, \det H \neq 0. This condition ensures that the phase function exhibits quadratic behavior in a neighborhood of x_0, allowing for a local Gaussian-like approximation without additional complications.^[15] At stationary points, the phase varies slowly compared to regions away from them, leading to coherent addition of contributions from nearby integrand values rather than destructive cancellation due to rapid oscillations. The signature of the Hessian, defined as the difference between the number of positive and negative eigenvalues, governs the phase shift in the asymptotic contribution, typically appearing as a factor of e^{i \pi \operatorname{sgn}(H)/4}.^[16]^[15] Degenerate stationary points, where higher-order derivatives vanish beyond the second (e.g., inflection points), necessitate uniform approximations such as those involving Airy functions to capture the asymptotic behavior accurately; however, the standard leading-order analysis assumes non-degeneracy. In integrals over manifolds, these stationary points are further classified by their Morse index—the number of negative eigenvalues of the Hessian—which determines topological properties relevant to the overall asymptotic expansion.^[17]^[18]

Leading-Order Approximation

General formula

The leading-order stationary phase approximation for the multidimensional oscillatory integral \int_D g(\mathbf{x}) \exp(i \lambda f(\mathbf{x})) \, d\mathbf{x}, where \mathbf{x} \in \mathbb{R}^n, D is a domain, g is a smooth amplitude function, f is a smooth real-valued phase function, and \lambda \to \infty, is given by summing the contributions from all non-degenerate interior stationary points \mathbf{x}_0 where \nabla f(\mathbf{x}_0) = 0 and the Hessian matrix \operatorname{Hess} f(\mathbf{x}_0) is invertible. The asymptotic formula is

\int_D g(\mathbf{x}) \exp(i \lambda f(\mathbf{x})) \, d\mathbf{x} \sim \sum_{\mathbf{x}_0 : \nabla f(\mathbf{x}_0)=0} g(\mathbf{x}_0) \exp\left(i \lambda f(\mathbf{x}_0) + i \frac{\pi}{4} \operatorname{sgn}(\operatorname{Hess} f(\mathbf{x}_0))\right) \left(\frac{2\pi}{\lambda}\right)^{n/2} \left|\det \operatorname{Hess} f(\mathbf{x}_0)\right|^{-1/2},

with the remainder term of order o(\lambda^{-n/2}). Here, \operatorname{sgn}(\operatorname{Hess} f(\mathbf{x}_0)) denotes the signature (number of positive eigenvalues minus number of negative eigenvalues) of the Hessian matrix at \mathbf{x}_0, and n is the dimension of the space. Each term in the sum arises from the Gaussian-like localization of the integrand near the stationary point \mathbf{x}_0, where the rapid oscillations away from \mathbf{x}_0 lead to destructive interference and negligible contribution. The prefactor (2\pi / \lambda)^{n/2} reflects the scaling of the effective volume of integration near each \mathbf{x}_0, which shrinks as \lambda^{-n/2} in n dimensions. The factor |\det \operatorname{Hess} f(\mathbf{x}_0)|^{-1/2} accounts for the curvature of the phase at \mathbf{x}_0, determining the width of the contributing region and thus the volume scaling in the local coordinates. This approximation assumes that the stationary points are non-degenerate (i.e., \det \operatorname{Hess} f(\mathbf{x}_0) \neq 0) and isolated, and that the amplitude g(\mathbf{x}) varies slowly compared to the phase oscillations near each \mathbf{x}_0. If the domain D has boundaries, additional contributions may arise from stationary points on the boundary, but the primary focus here is on interior points.

One-dimensional case

In the one-dimensional case, the leading-order stationary phase approximation simplifies to an explicit formula involving the second derivative at the stationary point. Consider the oscillatory integral I(\lambda) = \int_a^b g(x) \exp(i \lambda f(x)) \, dx, where f and g are smooth functions, g has compact support or appropriate decay, and \lambda \to \infty. The main contributions arise from isolated interior stationary points x_0 where f'(x_0) = 0 and f''(x_0) \neq 0. The asymptotic approximation is

I(\lambda) \sim \sum_{x_0 : f'(x_0)=0} g(x_0) \exp\left( i \lambda f(x_0) + i \frac{\pi}{4} \sgn(f''(x_0)) \right) \sqrt{ \frac{2\pi}{\lambda |f''(x_0)|} },

which serves as the one-dimensional limit of the general multidimensional formula. Near a stationary point x_0, the phase function f(x) is locally approximated by its quadratic Taylor expansion f(x) \approx f(x_0) + \frac{1}{2} f''(x_0) (x - x_0)^2, transforming the integral into a shifted Gaussian or Fresnel-type integral over the deviation variable u = x - x_0. Substituting this approximation yields I(\lambda) \approx g(x_0) \exp(i \lambda f(x_0)) \int_{-\infty}^{\infty} \exp\left( i \frac{\lambda f''(x_0)}{2} u^2 \right) \, du, whose evaluation provides the square-root scaling and phase factor. The relative error in this leading-order approximation is o(\lambda^{-1/2}), reflecting the width of the contributing region around each stationary point, which scales as \lambda^{-1/2}. This holds under the non-degeneracy assumption f''(x_0) \neq 0, ensuring the quadratic term dominates the local behavior. The phase shift \exp\left( i \frac{\pi}{4} \sgn(f''(x_0)) \right) arises from the orientation of the quadratic curvature and is rigorously justified by deforming the integration contour in the complex plane to align with the direction of steepest descent, where the integral reduces to a standard Gaussian form.

Derivation

Reduction steps

The reduction of a multidimensional oscillatory integral to a sum of one-dimensional contributions in the stationary phase approximation proceeds through a series of geometric and analytic steps centered on the non-degenerate stationary points of the phase function. These steps transform the original integral over \mathbb{R}^n or a manifold into localized one-dimensional integrals along the principal directions of oscillation, with controlled errors from regions away from the stationary points. The first step involves identifying all non-degenerate stationary points x_0 in the integration domain, where the gradient of the phase function \phi(x) vanishes, \nabla \phi(x_0) = 0, and the Hessian matrix H = \nabla^2 \phi(x_0) is invertible, ensuring \det H \neq 0. This condition guarantees that the phase has a quadratic behavior locally near each x_0, allowing the dominant contributions to the asymptotic expansion to arise solely from these isolated points.^[19] Next, a local change of coordinates is performed around each stationary point to simplify the phase. A linear transformation rotates the coordinate system so that the axes align with the eigenvectors of the Hessian, diagonalizing the quadratic form of the phase expansion: \phi(x) \approx \phi(x_0) + \frac{1}{2} \sum_{i=1}^n \lambda_i y_i^2, where \lambda_i are the eigenvalues of H. This separability exploits the spectral decomposition of the Hessian, reducing cross terms in the phase. With the phase diagonalized, Fubini's theorem is applied to separate the multidimensional integral into a product of one-dimensional oscillatory integrals along each eigenvector direction. The amplitude function is approximated as constant near x_0, effectively reducing the integral to a product of one-dimensional oscillatory integrals.^[1] Finally, the integration domain is localized to small neighborhoods around each x_0 using smooth cutoff functions that are unity near the point and decay rapidly outside. The error from the complementary regions is controlled via integration by parts, leveraging the non-stationary phase principle, which yields higher-order decay in the asymptotic parameter. This localization ensures the approximation captures all leading contributions without significant truncation error.^[19] For integrals over manifolds, the reduction employs normal coordinates adapted to the Riemannian structure at the stationary point, where the Morse lemma facilitates the transformation of the quadratic form into a canonical sum of squares (or differences thereof), preserving the separability while accounting for the manifold's geometry. This application of the Morse lemma provides a rigorous local flattening, enabling the same separation as in the Euclidean case.

Asymptotic methods

In asymptotic analysis of oscillatory integrals following the preparatory localization to neighborhoods of stationary points, the non-stationary regions are handled via repeated integration by parts.^[4] For an integral of the form \int \exp(i \lambda f(x)) a(x) \, dx where f'(x) \neq 0 and |f'(x)| \geq c > 0, integration by parts applied to \exp(i \lambda f(x)) / (i \lambda f'(x)) yields a boundary term and a remainder, with each iteration producing a decay of order O(1/\lambda).^[1] Multiple applications allow for arbitrary polynomial decay O(\lambda^{-N}) in these regions, provided the amplitude a(x) and phase derivatives remain bounded.^[4] For stationary regions in one dimension, where f'(t_0) = 0 but f''(t_0) \neq 0, the contour of integration is deformed in the complex plane to follow the path of steepest descent, transforming the oscillatory phase into a decaying exponential.^[1] This deformation facilitates evaluation of the leading contribution as a Gaussian integral, yielding the asymptotic term

\sqrt{\frac{2\pi}{\lambda |f''(t_0)|}} \exp\left(i \lambda f(t_0) + i \frac{\pi}{4} \operatorname{sgn}(f''(t_0))\right).

^[4] The multidimensional analog employs a change to Laplace-type coordinates that diagonalize the quadratic form given by the Hessian matrix \nabla^2 f(x_0) at the non-degenerate stationary point x_0, where \det(\nabla^2 f(x_0)) \neq 0.^[4] This reduction separates the integral into a product of one-dimensional Gaussian integrals along the principal directions, with the leading asymptotic determined by the determinant of the Hessian and an overall phase factor involving its signature.^[1] The van der Corput lemma complements these techniques by providing uniform frequency-based bounds for oscillatory integrals lacking stationary points, such as |\int_a^b \exp(i \lambda \phi(x)) \, dx| \leq C \lambda^{-1/k} when |\phi^{(k)}(x)| \geq 1 for k \geq 1. Modern microlocal refinements developed in the 1980s, notably by Hörmander, incorporate phase-space localization to handle degenerate or variable-order stationary points, yielding precise asymptotic expansions and wavefront set characterizations for the resulting distributions. Error control in these methods ensures that the remainder after localization to stationary regions decays as o(\lambda^{-n/2}), where n is the dimension, with the non-stationary contributions and higher-order terms confined to smaller orders.^[4]

Advanced Topics

Higher-order terms

The stationary phase approximation admits a full asymptotic expansion for the oscillatory integral I(\lambda) = \int_{\mathbb{R}^n} f(\mathbf{x}) e^{i \lambda \phi(\mathbf{x})} \, d^n \mathbf{x} as \lambda \to \infty, given by

I(\lambda) \sim \sum_{k=0}^\infty \lambda^{-(n/2 + k)} a_k,

where the leading coefficient a_0 is determined by the value of f and the Hessian determinant of \phi at a non-degenerate stationary point \mathbf{x}_0 with \nabla \phi(\mathbf{x}_0) = \mathbf{0} and \det H(\mathbf{x}_0) \neq 0, H being the Hessian matrix of \phi. The higher-order coefficients a_k for k \geq 1 incorporate contributions from higher derivatives of \phi and f; specifically, a_1 involves third-order derivatives of \phi or the gradient of f evaluated at \mathbf{x}_0. This series structure arises from Taylor expanding \phi and f around \mathbf{x}_0 and performing a change of variables to quadratic normal form, followed by term-by-term asymptotic evaluation.^[20]^[1] In the one-dimensional case (n=1), the coefficients a_k are computed via recursive integration by parts, where the phase is expanded as \phi(t) = \phi(t_0) + \frac{1}{2} \phi''(t_0) (t - t_0)^2 + \sum_{j=3}^\infty \frac{\phi^{(j)}(t_0)}{j!} (t - t_0)^j, and each integration by parts isolates powers of \lambda^{-1/2} while handling the remainder. Alternatively, matched asymptotic expansions align inner and outer solutions near the stationary point for systematic coefficient extraction. In multiple dimensions, computation is more complex; for well-behaved amplitudes f, techniques such as the cumulant expansion generate the series by treating deviations from the quadratic phase as perturbations, analogous to how Feynman diagrams organize higher-order cumulants in perturbative quantum field theory, where the leading saddle-point contribution corresponds to tree-level diagrams and subleading terms to loops encoding quantum fluctuations. This diagrammatic analogy clarifies the structure of higher cumulants in the expansion, as detailed in treatments of wave phenomena.^[20]^[21] Truncating the series at order m yields an approximation with remainder O(\lambda^{-(n/2 + m + \varepsilon)}) for some \varepsilon > 0, assuming isolated non-degenerate stationary points and sufficient smoothness of f and \phi. The error bound is uniform in the integration domain provided there is no coalescence of stationary points, where higher-order degeneracies would require rescaling and alter the expansion powers.^[1]^[20]

Multidimensional extensions

In the multidimensional setting, the stationary phase approximation extends beyond non-degenerate critical points, where the Hessian matrix of the phase function has full rank, to handle degenerate cases in which the Hessian has rank less than the dimension n. For such degenerate Hessians, normal form theory is employed to classify the local singularity of the phase function near the critical point and derive tailored asymptotic expansions that capture the reduced contribution from the oscillatory integral. This approach transforms the phase into a canonical form via coordinate changes, allowing for precise leading-order terms even when the standard Gaussian approximation fails.^[22] Degenerate cases often arise when stationary points coalesce, particularly in applications like wave diffraction where the geometry leads to caustics. In these scenarios, uniform asymptotic approximations bridge the gap between separated and merged points; for instance, the Airy function provides a transitional approximation for the fold caustic involving two coalescing stationary points, while the Pearcey integral offers a uniform expansion for the cusp catastrophe with three coalescing points, ensuring accuracy across the parameter space where the number of contributing critical points varies. These special functions arise from rescaling the phase near the degeneracy and evaluating the resulting canonical integrals, which dominate the asymptotics in higher dimensions.^[23]^[24] The method further generalizes to oscillatory integrals over manifolds, such as \int_M g \exp(i \lambda f) \, d\mu, where M is a smooth manifold equipped with measure d\mu. Local charts are used to parametrize neighborhoods of critical points on M, reducing the integral to a Euclidean form locally, while global phase consistency is maintained by incorporating the Maslov index, which tracks topological phase shifts (multiples of \pi) encountered when crossing caustics or Lagrangian submanifolds. This extension is crucial for semiclassical approximations on curved spaces, ensuring the overall phase aligns with the geometry of the problem.^[25] Recent developments since 2020 have integrated these multidimensional extensions with microlocal analysis to study oscillatory integral operators associated with partial differential equations, such as the Schrödinger equation, where degenerate phases model microlocal singularities in the solution's wavefront set. These advances refine operator bounds and propagation estimates, addressing challenges in dispersive PDEs by combining normal forms with pseudodifferential techniques for improved regularity results.^[26] For practical computation of multidimensional oscillatory integrals, numerical strategies focus on identifying stationary points through gradient-based optimization or Newton methods, followed by localized quadrature schemes—such as Gaussian or adaptive rules—applied in the vicinity of each critical point to approximate the dominant contributions, while non-stationary regions are downweighted due to rapid oscillations. This hybrid approach balances accuracy and efficiency, especially for high-dimensional phases encountered in scientific simulations.^[27]

Applications and Examples

Illustrative examples

One of the simplest illustrative examples of the stationary phase approximation is the Gaussian integral, which serves as a benchmark for the leading-order formula in one dimension. Consider the integral

I(\lambda) = \int_{-\infty}^{\infty} \exp\left(i \lambda \frac{x^2}{2}\right) \, dx,

where \lambda > 0 is a large real parameter. The phase function is \phi(x) = x^2 / 2, with amplitude a(x) = 1. The stationary point occurs where \phi'(x) = x = 0, and the second derivative is \phi''(0) = 1 > 0. Applying the leading-order stationary phase approximation yields

I(\lambda) \approx \sqrt{\frac{2\pi}{\lambda}} \exp\left(i \frac{\pi}{4}\right),

which matches the exact value of the Fresnel integral. This result demonstrates how the approximation captures the dominant contribution near the stationary point, with the phase shift \pi/4 arising from the positive curvature of the phase.^[2] A concrete one-dimensional application arises in evaluating oscillatory integrals like that for the Bessel function J_0(x), which has an integral representation as a Fourier-type transform:

J_0(x) = \frac{1}{2\pi} \int_0^{2\pi} \exp\left(i x \cos \theta\right) \, d\theta = \frac{1}{\pi} \int_0^{\pi} \cos(x \cos \theta) \, d\theta.

For large x > 0, the phase \phi(\theta) = \cos \theta has stationary points at \theta = 0 and \theta = \pi, where \phi'( \theta ) = -\sin \theta = 0 and \phi''(0) = -1 < 0, \phi''(\pi) = 1 > 0. Near these points, expand \cos \theta \approx 1 - \theta^2 / 2 around \theta = 0, leading to a Gaussian-like integral \int \exp(i x (1 - \theta^2 / 2)) \, d\theta \approx \exp(i x) \sqrt{2\pi / x} \exp(-i \pi/4). Accounting for contributions from both endpoints and symmetry, the leading-order approximation is

J_0(x) \approx \sqrt{\frac{2}{\pi x}} \cos\left(x - \frac{\pi}{4}\right).

This step-by-step computation highlights the endpoint stationary phase (valid here due to the finite interval), where the oscillation is slowest at the boundaries, yielding the characteristic $1/\sqrt{x} decay and phase shift.^[2] In wave propagation, the stationary phase method elucidates the asymptotic form of dispersive wave packets. Consider the integral representation of a wave field:

\eta(x, t) = \int_{-\infty}^{\infty} A(k) \exp\left(i (k x - \omega(k) t)\right) \, dk,

where \omega(k) is the dispersion relation, A(k) is the amplitude spectrum, and x, t > 0 are large. The phase is \phi(k) = k x - \omega(k) t, with stationary point k_s satisfying \phi'(k_s) = x - t \omega'(k_s) = 0, or x/t = \omega'(k_s) = c_g(k_s), the group velocity. Taylor expanding around k_s, \phi(k) \approx \phi(k_s) + (1/2) \phi''(k_s) (k - k_s)^2, with \phi''(k_s) = -t \omega''(k_s). The leading-order approximation is then

\eta(x, t) \approx A(k_s) \exp\left(i (k_s x - \omega(k_s) t)\right) \sqrt{\frac{2\pi}{t |\omega''(k_s)|}} \exp\left(i \frac{\pi}{4} \operatorname{sgn}(\omega''(k_s))\right).

This shows how the wave packet propagates at the group velocity c_g(k_s), with amplitude decaying as $1/\sqrt{t} due to dispersion.^[28]

Practical applications

The stationary phase approximation finds extensive application in optics, particularly for analyzing diffraction patterns. In the geometrical theory of diffraction developed by Joseph B. Keller, the method is employed to extend geometrical optics by incorporating diffracted rays that account for edge and vertex diffraction effects, providing asymptotic approximations for high-frequency wave propagation.^[29] This approach has been foundational in predicting radar cross-sections and antenna patterns, where it simplifies complex integral representations of scattered fields.^[30] In quantum mechanics, the stationary phase approximation underpins semiclassical methods such as the WKB (Wentzel-Kramers-Brillouin) approximation, which estimates solutions to the Schrödinger equation for slowly varying potentials. By applying the method to path integrals or turning-point analyses, it yields approximate energy levels and tunneling probabilities in systems like atomic potentials or molecular vibrations.^[31] This variant has enabled practical computations of bound-state spectra in quantum systems where exact solutions are intractable.^[32] Seismology utilizes the stationary phase approximation for modeling wave propagation in heterogeneous Earth models, as detailed in the comprehensive treatment by Keiiti Aki and Paul G. Richards. The method approximates travel times and waveforms by evaluating contributions from stationary points in phase integrals, facilitating the inversion of seismic data for subsurface velocity structures.^[33] In large-scale simulations, it enhances efficiency in ray tracing and synthetic seismogram generation for earthquake hazard assessment.^[34] In engineering, particularly signal processing and Fourier optics, the approximation accelerates computations of oscillatory integrals in diffraction-limited systems. For instance, it provides analytic estimates for Fresnel integrals in image reconstruction, reducing numerical overhead in optical design software.^[35] Recent extensions support time-frequency decompositions for chirp signals in radar and communications.^[36] Beyond physics and engineering, the method contributes to number theory through asymptotic evaluations of the Riemann zeta function, notably in the Riemann-Siegel formula, which approximates values on the critical line via saddle-point integrations.^[37] This has implications for zero-distribution estimates and prime number theorems. Emerging uses include applying the method to evaluations of quantum rate constants in molecular dynamics.^[38]

References

[1]
[PDF] Chapter 3 Asymptotic Expansion of Integrals - UC Davis Math
The method of stationary phase provides an asymptotic expansion of integrals with a rapidly oscillating integrand. ecause of cancelation, the behavior of such ...
[2]
[PDF] Method of stationary phase, Physics 2400
May 21, 2017 · The method of stationary phase gives the leading asymptotic behavior of generalized Fourier integrals having stationary points, ψ0 = 0.
[3]
Kelvin's method of stationary phase? - ScienceDirect.com
The method of stationary phase is a standard technique for estimating the value of an oscillatory integral when a certain parameter in the integrand becomes ...
[4]
[PDF] LECTURE NOTES 8 FOR 247B 1. Oscillatory integrals A basic ...
stationary phases, as long as the phase is quadratic at the stationary point: Lemma 2.8 (Asymptotic expansion for non-degenerate phases). Let a be a bump.
[5]
DLMF: §2.3 Integrals of a Real Variable ‣ Areas ‣ Chapter 2 ...
Mar 2, 2012 · For an extension with more general t -powers see Bleistein and Handelsman (1975, §4.1) . Other types of singular behavior in the integrand can ...
[6]
[PDF] ASYMPTOTIC METHODS - DAMTP
The method of stationary phase can be generalised to the case where φ has higher-order stationary points at some c ∈ (a, b). In that case, one needs to evaluate ...
[7]
[PDF] Stokes's Fundamental Contributions to Fluid Dynamics
Jul 12, 2019 · In this paper, Stokes begins to consider friction in fluid flow, discussing no-slip boundary conditions and their consequences. The issue of ...
[8]
[PDF] Stokes' phenomenon ; smoothing a victorian discontinuity - Numdam
He also, almost in passing, invented what we now call the method of stationary phase. The integral (13) has two stationary points, at s = ± (— z)112. When z ...
[9]
Lord Kelvin on fluid mechanics
May 29, 2012 · Here, so briefly demonstrated, is Kelvin's method of stationary phase, of great use in asymptotics. He does not trouble to give a rigorous ...
[10]
[PDF] Riemann Zeta Function - Gauge
Riemann's 1859 Zeta paper defines the Zeta function and uses its properties to approximate the count of prime numbers up to a number t , and the density of ...
[11]
https://academic.oup.com/book/54989
[12]
(PDF) Courant, Hilbert - Methods of Mathematical Physics Vol. II
The present volumes represent Richard Courant's second and final revision of 1953. It will be of interest to mathematicians and physicists. ... Read more ...
[13]
Lars Hormander | What's new - Terry Tao
Nov 30, 2012 · ... stationary phase computation, shows (again heuristically ... microlocal analysis and which was extensively developed in Hörmander's work.
[14]
[PDF] Course notes - OSU Math
1.3 Oscillatory integrals and the stationary phase method. In this setting, an integral of a function against a rapidly oscillating expo- nential becomes ...<|control11|><|separator|>
[15]
[PDF] lecture 5 — 09/30/2020 the method of stationary phase
Sep 30, 2020 · In the case the phase function is complex, there is a similar way to evaluate the asymptotic of the integral: the method of steepest descent.
[16]
Method of Stationary Phase - Richard Fitzpatrick
The method of stationary phase estimates integrals by finding points where the phase has a vanishing derivative, evaluating the integral near these points, and ...Missing: coherence | Show results with:coherence
[17]
Stationary phase, method of the - Encyclopedia of Mathematics
Mar 4, 2022 · The method of the stationary phase has an operator variant: λ=A, where A is the infinitesimal operator of the strongly-continuous group {eitA} ...Missing: historical development
[18]
[PDF] arXiv:1701.01183v1 [math.DG] 5 Jan 2017
Jan 5, 2017 · In §4 we will discuss the Morse-Bott Lemma and stationary phase approximation in the setting of supermanifolds. In §5, we state and prove the ...
[19]
Multidimensional Stationary Phase an Alternative Derivation
The method of multidimensional stationary phase is derived via a technique which makes strong use of integration by parts.
[20]
Multidimensional stationary phase approximation - ScienceDirect.com
These results are then used in a rigorous derivation of the asymptotic expansion of the oscillatory integral I(λ) = ∫Dg(x) eiλf(x) dx (x \te; r n), where a ...Missing: formula | Show results with:formula
[21]
[PDF] ASYMPTOTIC METHODS - DAMTP
The method of stationary phase can be generalised to the case where φ has higher-order stationary points at some c ∈ (a, b). In that case, one needs to ...
[22]
Mathematical Methods for Wave Phenomena - ScienceDirect.com
Mathematical Methods for Wave Phenomena focuses on the methods of applied mathematics, including equations, wave fronts, boundary value problems, and ...Missing: higher- | Show results with:higher-
[23]
An invariant method of stationary phase - MathOverflow
Nov 27, 2012 · A complete expansion is available in Hormander ALPDO first volume, Chapter 7, section devoted to the complex stationary phase method. To sum-up ...Is there something wrong with Hörmander's theorem on stationary ...Applications of microlocal analysis? - MathOverflowMore results from mathoverflow.net
[24]
Convergent and asymptotic expansions of the Pearcey integral
The simplest integral with only one free parameter, that corresponding to the fold catastrophe, involves two coalescing stationary points: the well-known ...
[25]
[PDF] The Pearcey integral in the highly oscillatory region - arXiv
Nov 17, 2015 · responds to the fold catastrophe, involves two coalescing stationary points: the well-known integral representation of the Airy function.
[26]
The Maslov integral representation of slowly varying dispersive ...
The Maslov integral provides an asymptotically valid solution, on the wavelength scale, to the initial value problem under conditions in which exact solutions ...
[27]
[PDF] H\"ormander oscillatory integral operators: a revisit - arXiv
May 6, 2025 · Variable coefficient Wolff-type inequalities and sharp local smoothing estimates for wave equations on manifolds. Analysis & PDE, 2020, 13(2): ...
[28]
[PDF] Approximate calculation of multidimensional first passage times - arXiv
Aug 29, 2025 · A final level of approximation is to perform stationary phase approximations for each of the two remaining integrals in either Eq.(31) or 34 ...
[29]
[PDF] Lecture 2: Gravity waves and the method of stationary phase
ks, a stationary point at which the derivative of the phase function S(k, x, t) ≡ kx − ω(k)t with respect to k vanishes. We approximate S about k = ks by ...<|control11|><|separator|>
[30]
Geometrical Theory of Diffraction* - Optica Publishing Group
The geometrical theory of diffraction is an extension of geometrical optics which accounts for diffraction. It introduces diffracted rays in addition to the ...
[31]
[PDF] Geometrical theory of diffraction - Indian Academy of Sciences
He using the stationary phase approximation further simplified the contour integral to contributions from a finite number of points on the boundary. Using a ...
[32]
[PDF] The WKB Method† 1. Introduction - University of California, Berkeley
The WKB method is important both as a practical means of approximating solutions to the. Schrödinger equation and as a conceptual framework for ...
[33]
[PDF] WKB quantization - ChaosBook.org
So the WKB approximation can interpreted either as a short wavelength/high frequency approximation to a wave-mechanical problem, or as the semiclassical, ħ 1 ...
[34]
(PDF) Aki&Richards_Quantitative_Seismology - Academia.edu
... stationary phase method to approximate the waveform. Similar to the steepest descents method, integration by the stationary phase method also exploits ...
[35]
[PDF] Quantitative Seismology Aki And Richards
The authors delve into Fourier transforms, spectral analysis, and filtering techniques that allow researchers to isolate meaningful signals from noisy data.
[36]
Stationary phase approximations in Fresnel-zone magnitude-only ...
In this paper we use the method of stationary phase[11]–[13] to obtain an analytic approximation to the first iteration in the case of reconstructing an image ...
[37]
[PDF] The Stationary Phase Approximation, Time-Frequency ... - arXiv
Aug 29, 2012 · Abstract. The principle of stationary phase (PSP) is re-examined in the context of linear time-frequency (TF) decomposition.<|control11|><|separator|>
[38]
[PDF] Riemann's Saddle-point Method and the Riemann-Siegel Formula
Riemann's way to calculate his zeta function on the critical line was based on an application of his saddle-point technique for approximating integrals that ...<|separator|>
[39]
[PDF] 1 Deep learning for ultrasound beamforming - arXiv
Sep 23, 2021 · The achievable resolution, contrast, and overall fidelity of B-mode imaging depends on the array aperture and geom- etry, element sensitivity ...<|separator|>
[40]
[PDF] Neural networks can be FLOP-efficient integrators of 1D oscillatory ...
Numerical integral methods crafted for highly oscillatory integrands have been developed. These are, for example, based on the stationary phase approximations ...
[41]
Stationary phase evaluations of quantum rate constants
We compute the quantum rate constant based on two extended stationary phase approximations to the imaginary-time formulation of the quantum rate theory.Missing: multidimensional | Show results with:multidimensional