Skip to content

Asymptotic Scales, Remainders, Uniformity, and Optimal Truncation

An expansion is asymptotic only after its limiting process, scale, region, and remainder have been specified. For every fixed number of retained terms, the remainder must be smaller than the last retained scale in the stated limit. Uniformity strengthens this statement by requiring one estimate to work over a whole parameter set. Optimal truncation is stronger again: it needs information about how the remainder bound depends on the truncation order. The symbol \sim by itself supplies none of that growing-order control.

This page develops those distinctions and then proves them in a zero-dimensional Euclidean ϕ4\phi^4 integral. The example gives a large-parameter expansion, an exact finite-order remainder bound, factorially growing coefficients, and a justified least-term prescription for that bound. It is a controlled analogue of a perturbative field integral, not a theorem about continuum QFT.

Helpful background. Limits, Completeness, and Modes of Convergence supplies the pointwise and uniform convergence language used in remainder estimates.

Start with a dimensionless parameter ϵ0+\epsilon\to0^+; other limits can be reduced to this form, for example by setting ϵ=1/Λ\epsilon=1/\Lambda when Λ+\Lambda\to+\infty. A complete asymptotic statement must say:

  • how the parameter approaches its limit;
  • which auxiliary variables are fixed and which are allowed to vary;
  • the real interval or complex sector on which the estimate holds;
  • the norm in which the remainder is measured;
  • and, for complex powers or logarithms, which branch is used.

For functions gg and hh, with hh nonzero sufficiently near the limit,

g(ϵ)=o ⁣(h(ϵ))g(ϵ)h(ϵ)0,g(ϵ)=O ⁣(h(ϵ))g(ϵ)h(ϵ) is bounded.\begin{aligned} g(\epsilon)=o\!\left(h(\epsilon)\right) &\quad\Longleftrightarrow\quad \frac{g(\epsilon)}{h(\epsilon)}\longrightarrow0,\\ g(\epsilon)=O\!\left(h(\epsilon)\right) &\quad\Longleftrightarrow\quad \left|\frac{g(\epsilon)}{h(\epsilon)}\right| \text{ is bounded.} \end{aligned}

Little-oo is a limiting statement. Big-OO asserts the existence of a bound, but the hidden constant and the neighborhood on which it holds still matter. Neither notation is a numerical error estimate until those quantities are known.

An ordered family {φn}n0\{\varphi_n\}_{n\geq0} is an asymptotic scale in the declared limit and region if

φn+1(ϵ)=o ⁣(φn(ϵ))(n=0,1,2,).\varphi_{n+1}(\epsilon) =o\!\left(\varphi_n(\epsilon)\right) \qquad(n=0,1,2,\ldots).

The standard power scale is 1,ϵ,ϵ2,1,\epsilon,\epsilon^2,\ldots as ϵ0\epsilon\to0. At a logarithmic threshold, or in an exponential problem, a different scale may be more natural. The ordering is part of the claim; a list of functions is not an asymptotic scale independently of a limit and region.

Poincaré expansions are fixed-order statements

Section titled “Poincaré expansions are fixed-order statements”

Define the NN-term partial sum and its remainder by

SN(ϵ)=n=0N1anφn(ϵ),RN(ϵ)=f(ϵ)SN(ϵ).\begin{aligned} S_N(\epsilon) &=\sum_{n=0}^{N-1}a_n\varphi_n(\epsilon),\\ R_N(\epsilon) &=f(\epsilon)-S_N(\epsilon). \end{aligned}

Then

f(ϵ)n=0anφn(ϵ)f(\epsilon)\sim \sum_{n=0}^{\infty}a_n\varphi_n(\epsilon)

means that, for every fixed N1N\geq1,

RN(ϵ)=o ⁣(φN1(ϵ))(ϵ0+).R_N(\epsilon) =o\!\left(\varphi_{N-1}(\epsilon)\right) \qquad(\epsilon\to0^+).

The quantifiers are decisive: choose NN, hold it fixed, and only then take the limit. Because the expansion is asserted to every order,

RN(ϵ)=aNφN(ϵ)+o ⁣(φN(ϵ)),R_N(\epsilon) =a_N\varphi_N(\epsilon) +o\!\left(\varphi_N(\epsilon)\right),

so RN=O(φN)R_N=O(\varphi_N) at each fixed order. For a finite expansion known only through φN1\varphi_{N-1}, however, the definition gives only RN=o(φN1)R_N=o(\varphi_{N-1}) unless a next-scale estimate is proved separately.

For a fixed nonvanishing scale and a fixed approach region, the coefficients are unique. They can be recovered recursively:

aN=limϵ0+f(ϵ)n=0N1anφn(ϵ)φN(ϵ).a_N= \lim_{\epsilon\to0^+} \frac{ f(\epsilon)-\displaystyle\sum_{n=0}^{N-1} a_n\varphi_n(\epsilon)} {\varphi_N(\epsilon)}.

The function represented by the coefficients is not unique. If A>0A>0, then

eA/ϵ=o(ϵN)for every fixed N.e^{-A/\epsilon}=o(\epsilon^N) \qquad\text{for every fixed }N.

Thus f(ϵ)f(\epsilon) and f(ϵ)+CeA/ϵf(\epsilon)+C e^{-A/\epsilon} have the same asymptotic power expansion as ϵ0+\epsilon\to0^+. The extra term is beyond all algebraic orders. This observation alone does not identify a saddle, a new physical sector, or a summation prescription.

The fixed-order definition, coefficient uniqueness, and this beyond-all-orders nonuniqueness are developed in Hunter 2004, Chapter 2, pp. 19–27, PDF.

For complex ϵ\epsilon, even this elementary statement needs a sector. On the closed subsector

argϵπ2δ,0<δ<π2,|\arg\epsilon|\leq\frac{\pi}{2}-\delta, \qquad 0<\delta<\frac{\pi}{2},

one has

eA/ϵexp ⁣(Asinδϵ),\left|e^{-A/\epsilon}\right| \leq \exp\!\left( -\frac{A\sin\delta}{|\epsilon|} \right),

uniformly. On the boundary argϵ=±π/2\arg\epsilon=\pm\pi/2, its modulus is 11, so the exponentially small estimate fails. Sector boundaries are therefore mathematical data, not decorative qualifiers.

Convergence and asymptoticity take different limits:

  • convergence studies SN(ϵ)S_N(\epsilon) as NN\to\infty with ϵ\epsilon fixed;
  • asymptoticity studies SN(ϵ)S_N(\epsilon) as ϵ0\epsilon\to0 with NN fixed.

An asymptotic series may converge, but it need not. If it diverges for every nonzero ϵ\epsilon, its partial sums can still approximate ff increasingly well for several orders before eventually getting worse. Conversely, a convergent series centered at one point need not be useful asymptotically in a different limit.

Accordingly,

f(ϵ)n=0anϵnf(\epsilon)\sim\sum_{n=0}^{\infty}a_n\epsilon^n

does not assert equality to an infinite sum. The right-hand side is a formal record of a family of remainder statements unless convergence or some separate summation procedure has been established.

Suppose f(ϵ,y)f(\epsilon,y) depends on an auxiliary parameter yYy\in Y. With a scale independent of yy, the expansion is uniform on YY if, for every fixed N1N\geq1,

supyYRN(ϵ,y)φN1(ϵ)0.\sup_{y\in Y} \left| \frac{R_N(\epsilon,y)} {\varphi_{N-1}(\epsilon)} \right| \longrightarrow0.

A common stronger result is a bound

supyYRN(ϵ,y)CN,YφN(ϵ),\sup_{y\in Y}|R_N(\epsilon,y)| \leq C_{N,Y}|\varphi_N(\epsilon)|,

where the constant is independent of yy. It may still depend strongly on NN, on YY, or on the distance from a sector boundary.

The elementary function

F(ϵ,y)=11+ϵy,ϵ>0,y0,F(\epsilon,y)=\frac{1}{1+\epsilon y}, \qquad \epsilon>0,\quad y\geq0,

exposes the distinction. The finite geometric identity gives

F(ϵ,y)=n=0N1(ϵy)n+(ϵy)N1+ϵy.F(\epsilon,y) =\sum_{n=0}^{N-1}(-\epsilon y)^n +\frac{(-\epsilon y)^N}{1+\epsilon y}.

On every fixed interval 0yM0\leq y\leq M,

sup0yMRN(ϵ,y)(ϵM)N,\sup_{0\leq y\leq M}|R_N(\epsilon,y)| \leq(\epsilon M)^N,

so the power expansion is uniform there. It is not uniform on [0,)[0,\infty). Already at leading order,

supy0F(ϵ,y)1=1\sup_{y\geq0} \left|F(\epsilon,y)-1\right|=1

for every ϵ>0\epsilon>0. The distinguished scaling y=c/ϵy=c/\epsilon makes the reason visible:

F(ϵ,c/ϵ)=11+c,F(\epsilon,c/\epsilon)=\frac{1}{1+c},

which is not 1+o(1)1+o(1). Near such a transition, the original ordering of terms has ceased to be uniform and a rescaled or matched description is needed.

For the corresponding need for uniform integral approximations near transition points, see Temme 1995, pp. 395–399.

The same caution applies in the complex plane. An expansion may be uniform on every closed subsector while its constants diverge as the subsector approaches a boundary. Powers such as zαz^\alpha and terms containing logz\log z also require a declared branch. Uniform control of the function itself does not automatically justify termwise differentiation or integration over an unbounded domain; those operations need derivative estimates, domination, or an appropriate norm bound.

Poincaré asymptoticity controls each fixed NN separately. It does not say how the constants in

RN(ϵ)CNφN(ϵ)|R_N(\epsilon)|\leq C_N|\varphi_N(\epsilon)|

grow with NN. Consequently, it says nothing by itself about a choice N=N(ϵ)N=N(\epsilon)\to\infty. A least-term rule becomes justified only when a theorem or an explicit calculation relates the remainder to the terms over the relevant growing range of NN.

This distinction between a formal fixed-order series and a justified least-term prescription is illustrated in Mariño 2026, § 2.1, pp. 2–4, PDF.

It helps to separate three objects:

  1. the actual optimal index, which minimizes RN|R_N| but is usually unknown;
  2. the index that minimizes a proved upper bound BNB_N;
  3. the least term of the formal series, which is only a proxy unless linked to the remainder.

For example, suppose one has proved, uniformly in the required region and for the needed range of NN,

RN(ϵ)BN(ϵ)=KΓ(N+β)(Aϵ)N,K>0,A>0.|R_N(\epsilon)| \leq B_N(\epsilon) =K\,\Gamma(N+\beta)(A\epsilon)^N, \qquad K>0,\quad A>0.

Here KK, AA, and the real number β\beta are fixed, and N+β>0N+\beta>0 throughout the range under consideration.

The ratio of consecutive bounds is

BN+1BN=Aϵ(N+β).\frac{B_{N+1}}{B_N} =A\epsilon(N+\beta).

The smallest bound therefore lies near

N(ϵ)1Aϵβ.N_*(\epsilon)\simeq \frac{1}{A\epsilon}-\beta.

Stirling’s formula then gives

BN=O ⁣((Aϵ)1/2βe1/(Aϵ)).B_{N_*} =O\!\left( (A\epsilon)^{1/2-\beta} e^{-1/(A\epsilon)} \right).

The exponentially small scale follows from the assumed NN-dependent bound, not from the definition of an asymptotic expansion. Without such control, even decreasing terms can give a misleading picture of the remainder; NIST DLMF 2026, § 2.11(i) gives explicit warnings and counterexamples.

A controlled zero-dimensional scalar benchmark

Section titled “A controlled zero-dimensional scalar benchmark”

Consider the normalized Euclidean integral

Z(Λ)=Λ2πRexp ⁣[Λ(ϕ22+ϕ424)]dϕ,Λ+.Z(\Lambda) = \sqrt{\frac{\Lambda}{2\pi}} \int_{\mathbb R} \exp\!\left[ -\Lambda\left( \frac{\phi^2}{2}+\frac{\phi^4}{24} \right) \right]\mathrm d\phi, \qquad \Lambda\to+\infty.

Here Λ>0\Lambda>0 is dimensionless. After x=Λϕx=\sqrt{\Lambda}\,\phi,

Z(Λ)=12πRex2/2ex4/(24Λ)dx.Z(\Lambda) = \frac{1}{\sqrt{2\pi}} \int_{\mathbb R} e^{-x^2/2} e^{-x^4/(24\Lambda)} \mathrm dx.

Equivalently, if XX is a standard normal random variable,

Z(Λ)=E ⁣[eX4/(24Λ)].Z(\Lambda) =\mathbb E\!\left[ e^{-X^4/(24\Lambda)} \right].

Expanding the second exponential through N1N-1 and evaluating the Gaussian moments gives

Z(Λ)=n=0N1(1)ncnΛn+RN(Λ),cn=E[X4n]24nn!=(4n)!96nn!(2n)!.\begin{aligned} Z(\Lambda) &= \sum_{n=0}^{N-1} (-1)^n c_n\Lambda^{-n} +R_N(\Lambda),\\ c_n &= \frac{\mathbb E[X^{4n}]}{24^n n!} = \frac{(4n)!} {96^n n!(2n)!}. \end{aligned}

Thus

Z(Λ)118Λ+35384Λ23853072Λ3+.Z(\Lambda) \sim 1-\frac{1}{8\Lambda} +\frac{35}{384\Lambda^2} -\frac{385}{3072\Lambda^3} +\cdots.

This is not merely a formal integration. For u0u\geq0, Taylor’s integral remainder is

eu=n=0N1(u)nn!+(u)N(N1)!01(1s)N1esuds.e^{-u} = \sum_{n=0}^{N-1}\frac{(-u)^n}{n!} + \frac{(-u)^N}{(N-1)!} \int_0^1 (1-s)^{N-1}e^{-su}\,\mathrm ds.

Substitute u=X4/(24Λ)u=X^4/(24\Lambda) and take the Gaussian expectation. The integrand in the remainder has a fixed sign and esu1e^{-su}\leq1, so for every integer N1N\geq1 and every Λ>0\Lambda>0,

0(1)NRN(Λ)cNΛN.0\leq (-1)^N R_N(\Lambda) \leq c_N\Lambda^{-N}.

This exact estimate proves the fixed-order asymptotic expansion and, in this example, bounds the error by the first omitted term.

Now let Tn=cnΛnT_n=c_n\Lambda^{-n} be the magnitude of the nnth term. Directly from the factorial formula,

Tn+1Tn=(4n+1)(4n+3)24(n+1)Λ2n3Λ.\frac{T_{n+1}}{T_n} = \frac{(4n+1)(4n+3)} {24(n+1)\Lambda} \sim \frac{2n}{3\Lambda}.

For every fixed Λ\Lambda, the terms eventually grow, so the formal series diverges. The first-omitted-term bound is smallest near

N3Λ2.N_*\sim\frac{3\Lambda}{2}.

Stirling’s formula makes the corresponding scale explicit:

cn1π2(23)nΓ(n).c_n \sim \frac{1}{\pi\sqrt2} \left(\frac{2}{3}\right)^n \Gamma(n).

Choosing the nearest integer to 3Λ/23\Lambda/2 therefore gives

RN(Λ)TN23πΛe3Λ/2.|R_{N_*}(\Lambda)| \leq T_{N_*} \sim \sqrt{\frac{2}{3\pi\Lambda}}\, e^{-3\Lambda/2}.

The final \sim describes the upper bound TNT_{N_*} as Λ+\Lambda\to+\infty; it does not assert that the actual remainder is asymptotic to that bound. This is precisely the extra evidence that the generic least-term slogan lacks.

The ordinary integral is often called zero-dimensional ϕ4\phi^4 theory. More concretely, E[X4n]=(4n1)!!\mathbb E[X^{4n}]=(4n-1)!! counts Gaussian Wick pairings, while 1/[n!(4!)n]1/[n!(4!)^n] is the factor obtained by expanding nn indistinguishable quartic vertices. The coefficients therefore reproduce the combinatorics of perturbative vacuum diagrams, so the example is genuinely QFT-facing; compare Zinn-Justin 2021, § 7.4, p. 132. But it has no spacetime modes, ultraviolet limit, renormalization scheme, or physical observable beyond the finite integral. The interpretation of control parameters, saddle sectors, loop counting, and failure conditions belongs to Saddles, Control Parameters, and Loop Counting.

A trustworthy asymptotic approximation should make the following information recoverable:

  1. Limit: Which dimensionless quantity tends to zero or infinity, and along what path?
  2. Scale: Which ordered functions define “successively smaller”?
  3. Region: Which interval, parameter set, or closed complex subsector is covered?
  4. Uniformity: Which variables are fixed, and which vary under a supremum or norm?
  5. Remainder: Is the claim little-oo, big-OO, or an explicit inequality?
  6. Order dependence: Are constants known when NN grows, or only for each fixed NN?
  7. Branches and boundaries: Which branch is used, and where can the estimate fail?

If an “optimal” truncation is quoted, also state whether it minimizes the true error, a proved bound, or only the displayed term magnitudes.

Reading \sim as an infinite equality. A Poincaré series records fixed-order remainder behavior. It need not converge or define a unique function.

Letting the truncation order drift silently. A proof for every fixed NN does not cover N=N(ϵ)N=N(\epsilon). Growing-order truncation requires bounds uniform in the relevant range of NN.

Dropping the uniformity set. Pointwise validity at every fixed auxiliary parameter does not imply one estimate works over an unbounded or transition region. Test distinguished scalings in which the auxiliary parameter depends on the small parameter.

Treating hidden constants as universal. A constant in O(φN)O(\varphi_N) may depend on NN, the parameter set, a branch, or the distance from a sector boundary.

Using the first omitted term without a theorem. The zero-dimensional example has a sign-definite integral remainder, so the rule is proved there. A generic asymptotic series does not inherit that bound.

Overinterpreting the benchmark. A finite-dimensional scalar integral illustrates perturbative combinatorics and large-order control; it does not establish continuum-QFT existence, renormalized error bounds, or physical nonperturbative sectors.

State the asymptotic-scale condition and the remainder condition for an NN-term Poincaré expansion. Which quantity is held fixed in the defining limit?

Solution

The scale satisfies φn+1=o(φn)\varphi_{n+1}=o(\varphi_n). With

SN=n=0N1anφn,RN=fSN,S_N=\sum_{n=0}^{N-1}a_n\varphi_n, \qquad R_N=f-S_N,

the expansion means RN=o(φN1)R_N=o(\varphi_{N-1}) for every fixed N1N\geq1. The truncation order is fixed while the asymptotic parameter approaches its limit; no NN\to\infty assertion follows.

2. Find information invisible to all powers

Section titled “2. Find information invisible to all powers”

For A>0A>0, prove that eA/ϵ=o(ϵN)e^{-A/\epsilon}=o(\epsilon^N) as ϵ0+\epsilon\to0^+ for every fixed NN. What changes if ϵ\epsilon approaches zero in the complex plane?

Solution

Set t=A/ϵt=A/\epsilon. Then

ϵNeA/ϵ=ANtNet0(t+).\epsilon^{-N}e^{-A/\epsilon} =A^{-N}t^Ne^{-t}\longrightarrow0 \qquad(t\to+\infty).

For complex ϵ\epsilon, decay along an approach requires Re(A/ϵ)+\operatorname{Re}(A/\epsilon)\to+\infty. This holds uniformly on any closed subsector argϵπ/2δ|\arg\epsilon|\leq\pi/2-\delta, but it fails on the imaginary-axis boundary. Hence the power coefficients can be unique even though the represented function is not.

For

F(ϵ,y)=11+ϵy,F(\epsilon,y)=\frac{1}{1+\epsilon y},

derive the exact NN-term geometric remainder. Decide whether the expansion is uniform on [0,M][0,M] and on [0,)[0,\infty).

Solution

The finite identity is

F(ϵ,y)=n=0N1(ϵy)n+(ϵy)N1+ϵy.F(\epsilon,y) =\sum_{n=0}^{N-1}(-\epsilon y)^n +\frac{(-\epsilon y)^N}{1+\epsilon y}.

For fixed MM, the remainder is bounded by (ϵM)N(\epsilon M)^N, so the expansion is uniform on [0,M][0,M]. On [0,)[0,\infty), the leading error has supremum 11 for every ϵ>0\epsilon>0, so uniformity fails. The scaling y=c/ϵy=c/\epsilon keeps ϵy\epsilon y of order one and exposes the transition.

4. Transfer the method to the scalar benchmark

Section titled “4. Transfer the method to the scalar benchmark”

Rescale the zero-dimensional integral, compute c1c_1 and c2c_2, derive Tn+1/TnT_{n+1}/T_n, and locate the least first-omitted-term bound.

Solution

The substitution x=Λϕx=\sqrt{\Lambda}\phi turns the integral into a standard-normal expectation. From

cn=(4n)!96nn!(2n)!,c_n=\frac{(4n)!}{96^n n!(2n)!},

one obtains

c1=18,c2=35384.c_1=\frac18, \qquad c_2=\frac{35}{384}.

For Tn=cnΛnT_n=c_n\Lambda^{-n},

Tn+1Tn=(4n+1)(4n+3)24(n+1)Λ2n3Λ.\frac{T_{n+1}}{T_n} = \frac{(4n+1)(4n+3)} {24(n+1)\Lambda} \sim\frac{2n}{3\Lambda}.

The bound is therefore least near N3Λ/2N_*\sim3\Lambda/2. Taylor’s sign-definite integral remainder is what converts this term calculation into a rigorous error bound. The conclusion applies to this regulated zero-dimensional model, not automatically to continuum QFT.

This page supplied a criterion and an error discipline, not a mechanism for deriving coefficients from a general integral. Laplace Method and Steepest Descent develops that mechanism, including saddle selection and contour geometry. Stationary Phase, Coalescing Saddles, and Stokes Phenomena treats oscillatory integrals and changing critical-point structure, while WKB, Eikonal Expansions, and Turning Points carries the same error discipline into differential equations. Large-Order Growth and the Borel Transform develops the later summation questions that a divergent series raises.