Asymptotic Expansions

⚡ Key Takeaways

Asymptotic expansions represent a function by a (often divergent) series whose partial sums approximate the function with controllable error in a limit (e.g. $x\to\infty$ or $\varepsilon\to 0$ ). Key result: a divergent series $\sum a_n x^{-n}$ can give superb numerical accuracy if truncated optimally, even though it never converges. Where it applies: singular perturbation theory, integral evaluation (Laplace, stationary phase, saddle point), special functions, WKB analysis, statistical mechanics. Who/when: Poincaré gave the rigorous definition (1886), building on Stirling, Euler, and Laplace; Stokes and others developed the analytic machinery.

💡 The Core Insight

A series need not converge to be useful. What matters is not the behavior of the partial sums as $N\to\infty$ at fixed $x$ , but the behavior of the remainder as $x\to x_0$ at fixed $N$ . Asymptotic expansions trade convergence for uniform smallness of the truncation error in a limiting regime, and the best accuracy is obtained by truncating at the smallest term (optimal truncation / superasymptotics).

DefinitionAsymptotic Sequence

A sequence of functions $\{\phi_n(x)\}_{n=0}^\infty$ is an asymptotic sequence (or scale) as $x\to x_0$ if for every $n$ ,

1

\quad \phi_{n+1}(x) = o\big(\phi_n(x)\big), \qquad x\to x_0,

that is $\lim_{x\to x_0}\phi_{n+1}(x)/\phi_n(x) = 0$ . Common scales: $\{x^{-n}\}$ as $x\to\infty$ , $\{(x-x_0)^n\}$ as $x\to x_0$ , $\{\varepsilon^n\}$ as $\varepsilon\to 0$ , and mixed scales like $\{\varepsilon^n\log^m\varepsilon\}$ .

DefinitionAsymptotic Expansion (Poincaré sense)

Let $\{\phi_n\}$ be an asymptotic sequence as $x\to x_0$ . The formal series $\sum_{n=0}^\infty a_n \phi_n(x)$ is an asymptotic expansion of $f(x)$ , written

2

\quad f(x) \sim \sum_{n=0}^\infty a_n \phi_n(x), \qquad x\to x_0,

if for every fixed $N\ge 0$ ,

3

\quad f(x) - \sum_{n=0}^{N} a_n \phi_n(x) = o\big(\phi_N(x)\big), \qquad x\to x_0.

Equivalently the remainder is of the order of the first omitted term:

4

\quad f(x) - \sum_{n=0}^{N-1} a_n \phi_n(x) = O\big(\phi_N(x)\big).

TheoremUniqueness of Coefficients

Given an asymptotic sequence $\{\phi_n\}$ , the coefficients $a_n$ in [2] are uniquely determined by $f$ via

5

\quad a_N = \lim_{x\to x_0} \frac{\,f(x) - \sum_{n=0}^{N-1} a_n\phi_n(x)\,}{\phi_N(x)}.

However, the function is not uniquely determined by its expansion: two distinct functions may share the same asymptotic series.

TheoremOperations on Asymptotic Series

If $f\sim\sum a_n x^{-n}$ and $g\sim\sum b_n x^{-n}$ as $x\to\infty$ , then:

Linearity: $\alpha f + \beta g \sim \sum(\alpha a_n + \beta b_n)x^{-n}$ .
Multiplication: $fg \sim \sum c_n x^{-n}$ with $c_n=\sum_{k=0}^n a_k b_{n-k}$ (Cauchy product).
Integration: $\int_x^\infty\!\big(f(t)-a_0-a_1 t^{-1}\big)\,dt$ is asymptotic to the termwise integral.
Differentiation is NOT generally valid without extra hypotheses (e.g. $f'$ must itself possess an asymptotic expansion).

mathematician · Henri Poincaré | 1854–1912

Poincaré supplied the rigorous definition (1886) in his study of celestial mechanics, where divergent series of perturbation theory demanded justification. He recognized that astronomers had been computing with divergent series for a century and that their utility lay in the order of the remainder, not convergence. His framing — now called the Poincaré definition — separates the analytic from the formal and remains the standard. The same insight powered his qualitative theory of differential equations and his work on the three-body problem.

ProofThe Asymptotic Series of an Integral — Watson/Laplace prototype

Consider the exponential integral and derive its expansion to show how [2]–[4] arise.

Step 1. Define, for $x>0$ ,

6

\quad f(x) = \int_0^\infty \frac{e^{-xt}}{1+t}\,dt.

Step 2. Repeated integration by parts (or expanding $1/(1+t)$ with remainder) gives the geometric-with-remainder identity

7

\quad \frac{1}{1+t} = \sum_{n=0}^{N-1}(-1)^n t^n + \frac{(-1)^N t^N}{1+t}.

Step 3. Insert [7] into [6] and use $\int_0^\infty t^n e^{-xt}\,dt = n!\,x^{-n-1}$ :

8

\quad f(x) = \sum_{n=0}^{N-1} \frac{(-1)^n n!}{x^{n+1}} + R_N(x),\qquad R_N(x)=(-1)^N\!\int_0^\infty \frac{t^N e^{-xt}}{1+t}\,dt.

Step 4. Bound the remainder using $1/(1+t)\le 1$ :

9

\quad |R_N(x)| \le \int_0^\infty t^N e^{-xt}\,dt = \frac{N!}{x^{N+1}}.

Thus $R_N = O(x^{-N-1})$ , which is the order of the first omitted term — exactly condition [4]. Hence

10

\quad f(x) \sim \sum_{n=0}^\infty \frac{(-1)^n n!}{x^{n+1}}, \qquad x\to\infty.

Step 5. Convergence test: the ratio of consecutive terms is $|a_{n+1}/a_n| = (n+1)/x \to\infty$ , so the series diverges for every $x$ . Yet [9] shows the error is rigorously controlled. ∎

∎

TheoremWatson's Lemma

If $g(t)\sim \sum_{n=0}^\infty a_n t^{(n+\lambda-\mu)/\mu}$ as $t\to 0^+$ (with $\lambda,\mu>0$ ) and the Laplace integral converges, then

11

\quad \int_0^\infty e^{-xt} g(t)\,dt \sim \sum_{n=0}^\infty a_n\,\frac{\Gamma\!\big((n+\lambda)/\mu\big)}{x^{(n+\lambda)/\mu}}, \qquad x\to\infty.

The dominant contribution comes entirely from a neighborhood of $t=0$ .

🧠 Intuition

Think of an asymptotic series as a sequence of successively finer rulers calibrated for a particular regime. Each term is much smaller than the last when $x$ is large, so adding terms improves the picture — up to a point. Because the coefficients $a_n = (-1)^n n!$ grow factorially while the powers $x^{-n}$ only shrink geometrically, eventually $|a_n x^{-n}|$ stops decreasing and starts growing. The smallest term occurs near $n\approx x$ , and that minimal term sets the intrinsic accuracy floor: you cannot do better than $\sim e^{-x}$ (the size of the least term) using the series alone. This exponentially small residual is the realm of exponential asymptotics / resurgence, where the "missing" information is hidden across the Stokes lines of the complex plane.

Figure

Optimal truncation: error decreases then explodes; least term sets the accuracy floor

ExampleExponential Integral at x = 10

Given: Approximate $f(10)=\int_0^\infty \frac{e^{-10t}}{1+t}\,dt$ using the divergent series [10]. The true value is $f(10)=0.0915633\ldots$

Steps. Terms $a_n = (-1)^n n!/10^{n+1}$ :

$n=0$ : $+0.100000$
$n=1$ : $-0.010000$ → partial $0.090000$
$n=2$ : $+0.002000$ → $0.092000$
$n=3$ : $-0.000600$ → $0.091400$
$n=4$ : $+0.000240$ → $0.091640$
$n=5$ : $-0.000120$ → $0.091520$
$n=6$ : $+0.0000720$ → $0.091592$
$n=7$ : $-0.0000504$ → $0.091542$
...

Least term near $n\approx 10$ : $|a_{10}|=10!/10^{11}=3.6\times10^{-5}$ . Truncating just before it (around $n=9$ , $|a_9|=3.6\times 10^{-5}$ ) gives the best accuracy.

Boxed answer.

\boxed{f(10)\approx 0.09156,\quad \text{error}\sim 3.6\times10^{-5}\ \text{(set by least term)}}

Adding more terms beyond $n\approx 10$ worsens the result — the divergence asserts itself.

ExampleStirling's Formula via Laplace's Method

Given: Find the leading asymptotics of $\Gamma(x+1)=x!=\int_0^\infty t^x e^{-t}\,dt$ as $x\to\infty$ .

Steps. 1. Write $t^x e^{-t}=e^{x\ln t - t}$ and substitute $t=xs$ : integrand $=e^{x(\ln(xs)-s)}=x^{x+1}e^{x(\ln s - s)}\cdot x^{-... }$ . Cleaner: $\Gamma(x+1)=x^{x+1}\int_0^\infty e^{x(\ln s - s)}\,ds$ . 2. The exponent $\phi(s)=\ln s - s$ has maximum at $s=1$ , $\phi(1)=-1$ , $\phi''(1)=-1$ . 3. Gaussian (Laplace) approximation around $s=1$ :

12

\quad \int_0^\infty e^{x\phi(s)}\,ds \approx e^{-x}\sqrt{\frac{2\pi}{x\,|\phi''(1)|}} = e^{-x}\sqrt{\frac{2\pi}{x}}.

4. Multiply by $x^{x+1}$ :

13

\quad \Gamma(x+1)\sim x^{x+1}e^{-x}\sqrt{\tfrac{2\pi}{x}} = \sqrt{2\pi x}\;x^x e^{-x}.

Numerical check, $x=10$ : Stirling gives $\sqrt{20\pi}\cdot 10^{10}e^{-10}=3598695.6$ ; true $10!=3628800$ . Relative error $0.83\%\approx 1/(12\cdot 10)$ , consistent with the next correction factor $(1+\tfrac{1}{12x}+\cdots)$ .

Boxed answer.

\boxed{n!\sim\sqrt{2\pi n}\,\left(\frac{n}{e}\right)^{n}\Big(1+\frac{1}{12n}+\frac{1}{288n^2}-\cdots\Big)}

⚠️ Common mistake

"Asymptotic series converge." They usually do not. Convergence concerns $N\to\infty$ at fixed $x$ ; asymptoticity concerns $x\to x_0$ at fixed $N$ . These are independent. A series can be convergent but not asymptotic, asymptotic but divergent (the common case), both, or neither. Never add "more terms" expecting indefinite improvement.

⚠️ Common mistake

"You can always differentiate an asymptotic expansion term by term." False without justification. The relation $f\sim\sum a_n\phi_n$ does not imply $f'\sim\sum a_n\phi_n'$ . Integration is safe (smoothing), differentiation is not (roughening). One must verify that $f'$ itself admits an expansion.

CounterexampleSame expansion, different functions

The function $g(x)=e^{-x}$ satisfies $e^{-x}=o(x^{-n})$ for every $n$ as $x\to\infty$ , so $e^{-x}\sim 0$ in the power scale $\{x^{-n}\}$ . Therefore $f(x)$ and $f(x)+e^{-x}$ have the identical asymptotic power series. The expansion cannot detect the exponentially small term — it lives "beyond all orders." This non-uniqueness (Theorem [5]) is precisely the gap that exponential asymptotics and resurgence theory aim to fill via Stokes phenomena.

connection

Asymptotic expansions are the analytic engine of Perturbation Theory (regular and singular), where solutions are expanded in a small parameter $\varepsilon$ . In Method of Matched Asymptotic Expansions one stitches together inner and outer expansions valid in different regions of a boundary layer. The WKB Approximation is an asymptotic expansion of solutions to differential equations in $\hbar$ (or wavelength). The evaluation of Laplace's Method, Method of Steepest Descent (saddle point), and Stationary Phase Approximation all produce asymptotic series for integrals. The deep structure connecting divergent series across Stokes Phenomenon is formalized by Resurgence Theory and Borel summation, linking asymptotics to Analytic Continuation. The factorial growth of coefficients ties directly to Borel Summation.

🔬 Where this appears

⏳ History

- 1730 (Stirling/De Moivre): Stirling's series for

n!

— historically the first widely used divergent asymptotic series.

- 1754 (Euler): Manipulates divergent series ("the series

1-1!+2!-\cdots

") obtaining numerically sensible values; intuits summation methods.

- 1812 (Laplace): Develops the method of approximating integrals by the contribution of their dominant peak — Laplace's method.

- 1847–1857 (Stokes): Discovers the Stokes phenomenon: asymptotic representations change discontinuously across rays in the complex plane; studies the Airy function.

- 1863 (Riemann): Saddle-point / steepest-descent ideas in his work, later named after him and Debye.

- 1886 (Poincaré): Gives the rigorous definition of an asymptotic expansion and the notion of an asymptotic series, establishing the field.

- 1909 (Debye): Systematizes the method of steepest descent for Bessel-function asymptotics.

- 1918 (Watson): Proves Watson's lemma, the workhorse for Laplace-type integrals.

- 1920s–30s (Borel, Hardy): Theory of divergent series and Borel summation; Hardy's Divergent Series (1949).

- 1980s–present (Dingle, Écalle, Berry): Resurgence theory, hyperasymptotics, and smoothing of the Stokes phenomenon (Berry's error-function transition) — recovering exponentially small terms beyond all orders.

open-problem

The exponentially small terms invisible to a power-series expansion ("beyond all orders") carry genuine physical and mathematical content — instantons, tunneling amplitudes, and selection rules in pattern formation. Resurgence theory (Écalle) provides a framework in which the divergence of a perturbative series encodes its own non-perturbative completion, but a fully general, constructive theory of resummation for arbitrary nonlinear PDEs remains incomplete. Open questions include: precise transseries representations for solutions of broad classes of differential equations; rigorous control of multi-instanton sectors in realistic quantum field theories; effective algorithms for hyperasymptotic summation past several Stokes lines; and the interplay between the Stokes phenomenon and singularities in the Borel plane for systems lacking known integral representations.

Formal Statement

Why It Holds / Derivation

Building Intuition

Worked Examples

Common Misconceptions

Connections

Real-World Applications

Historical Development

Subtleties & Limits

See Also