Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Subexponentiality

Statement

Let X1,X2,…X_1,X_2,\dots be iid nonnegative random variables with distribution FF and unbounded right support, so Fˉ(x)=P(X>x)>0\bar F(x)=\mathbb P(X>x)>0 for every finite xx. The distribution is subexponential when, for two independent copies,

P(X1+X2>x)2P(X>x)→1,x→∞.\frac{\mathbb P(X_1+X_2>x)}{2\mathbb P(X>x)} \to 1, \qquad x\to\infty.

Equivalently, for every fixed integer n≥2n\ge2 (and trivially for n=1n=1), with Sn=X1+⋯+XnS_n=X_1+\cdots+X_n and Mn=max⁡1≤i≤nXiM_n=\max_{1\le i\le n}X_i,

P(Sn>x)∼P(Mn>x)∼nFˉ(x).\mathbb P(S_n>x) \sim \mathbb P(M_n>x) \sim n\bar F(x).

For those same fixed integers, equivalently,

P(Sn>x,  Mn≤x)=o(Fˉ(x)).\mathbb P(S_n>x,\; M_n\le x)=o(\bar F(x)).

This is the one-big-jump principle: for a fixed number of summands, a very large sum is usually caused by one summand being very large, not by all summands moving moderately together.

The maximum-tail asymptotic P(Mn>x)∼nFˉ(x)\mathbb P(M_n>x)\sim n\bar F(x) is not the distinctive part of the theorem: it follows for every iid distribution with an unbounded right tail and fixed nn. The subexponential content is that the sum tail has the same first-order asymptotic as the maximum tail.

Regularly varying survival tails with positive exponent, Fˉ∈RV−α\bar F\in RV_{-\alpha} for α>0\alpha>0, are a standard subexponential class. We cite that theorem and use the exact Pareto family as a checkable example.

Here FF is the common CDF, nn is fixed while xx grows, and α>0\alpha>0 and xm>0x_m>0 are the exponent and cutoff in the Pareto example. The notation a(x)∼b(x)a(x)\sim b(x) means a(x)/b(x)→1a(x)/b(x)\to1; o(b(x))o(b(x)) denotes a term negligible after division by b(x)b(x).

One-big-jump intuition

Subexponential does not mean “lighter than exponential,” and it is unrelated to the concentration-theory phrase “sub-exponential random variable” for a light-tailed ψ1\psi_1 condition. In this context it means the opposite: the tail is heavy enough that convolution barely changes the leading-order tail probability.

For two iid summands, a sufficiently large sum is about twice as likely as one summand exceeding the same threshold. For nn fixed summands, the leading factor is nn.

This is the mathematical version of a recurring Incerto warning. In thin-tail settings, a large aggregate deviation is often the result of many small coordinated deviations. In subexponential settings, the aggregate tail is usually dominated by a single large term. That changes how sums, ruin events, insurance losses, and sample moments should be read.

What we can prove directly

The statement collects several equivalent forms and one important sufficient class. We prove the elementary maximum-tail asymptotic and an exact Pareto two-summand check. The full equivalence theory and the theorem that regularly varying tails are subexponential are cited, because their proofs need more convolution-tail machinery than we develop here.

For the general fixed-nn extension and its converse, we cite Foss, Korshunov, and Zachary (2013), §3.3, Corollary 3.20 and Theorem 3.21, p. 52. Their nonnegative case agrees with the definition here. The fixed-size Pareto proof for arbitrary nn is already on the Pareto page.

The theorem that regularly varying tails are subexponential goes back to Chistyakov’s convolution-tail criterion for sums of independent positive random variables 1964. Standard modern references with proofs are Embrechts, Klueppelberg, and Mikosch, Appendix A3 1997, and Foss, Korshunov, and Zachary, second edition, §3.5, Theorem 3.29 and its following paragraph, p. 57 2013 for the regular-variation implication. We do not reprove the full theorem, but the exact Pareto case gives the right calibration.

Universal maximum tail

For any distribution with unbounded right tail, Fˉ(x)→0\bar F(x)\to0 as x→∞x\to\infty. For fixed nn, iid independence gives

P(Mn≤x)=P(X1≤x,…,Xn≤x)=F(x)n=(1−Fˉ(x))n.\mathbb P(M_n\le x) =\mathbb P(X_1\le x,\dots,X_n\le x) =F(x)^n =\left(1-\bar F(x)\right)^n.

Therefore

P(Mn>x)=1−(1−Fˉ(x))n\mathbb P(M_n>x) =1-\left(1-\bar F(x)\right)^n

Set u=Fˉ(x)u=\bar F(x). Since u→0u\to0 and nn is fixed, the binomial expansion gives

1−(1−u)n=nu−(n2)u2+O(u3),1-(1-u)^n =nu-\binom{n}{2}u^2+O(u^3),

and therefore

1−(1−u)nnu=1−n−12u+O(u2)→1.\frac{1-(1-u)^n}{nu} =1-\frac{n-1}{2}u+O(u^2) \to1.

Substituting back,

P(Mn>x)∼nFˉ(x).\mathbb P(M_n>x) \sim n\bar F(x).

This maximum formula only uses iid independence, fixed nn, and Fˉ(x)→0\bar F(x)\to0. For nonnegative variables, Mn>xM_n>x implies Sn=X1+⋯+Xn>xS_n=X_1+\cdots+X_n>x, so it is a lower bound for the sum tail. What is special about subexponential tails is the matching upper bound:

P(Sn>x,  Mn≤x)=o(Fˉ(x)).\mathbb P(S_n>x,\;M_n\le x)=o(\bar F(x)).

Those are the configurations where several observations are moderately large but none exceeds xx. Subexponentiality says this residual event is negligible at the scale of one tail probability.

Pareto upper-bound check

For the Pareto tail Fˉ(x)=(xm/x)α\bar F(x)=(x_m/x)^\alpha with x≥xmx\ge x_m, the missing upper-bound step can be checked directly for n=2n=2. Let

Ax={X1+X2>x,  M2≤x}.A_x=\{X_1+X_2>x,\;M_2\le x\}.

Choose h=xm(x/xm)γh=x_m(x/x_m)^\gamma with 1/2<γ<11/2<\gamma<1. This keeps the units consistent and gives h/x→0h/x\to0 and h→∞h\to\infty.

The residual event AxA_x is the region in [0,x]2[0,x]^2 above X1+X2=xX_1+X_2=x. It is covered by two thin strips and a region where both coordinates exceed hh. The following table replaces the geometric plot.

Cover eventProbability or bound
x−h<X1≤xx-h<X_1\le xFˉ(x−h)−Fˉ(x)\bar F(x-h)-\bar F(x)
x−h<X2≤xx-h<X_2\le xFˉ(x−h)−Fˉ(x)\bar F(x-h)-\bar F(x)
X1>hX_1>h and X2>hX_2>hFˉ(h)2\bar F(h)^2 by independence

Outside the cover region, a point has neither coordinate in (x−h,x](x-h,x] and not both coordinates above hh. Thus at least one coordinate is ≤h\le h, and the other is ≤x−h\le x-h, so X1+X2≤xX_1+X_2\le x. Therefore any point in AxA_x must lie in the cover: either one summand lies in (x−h,x](x-h,x], or both summands exceed hh. Hence

P(Ax)≤2P(x−h<X≤x)+P(X1>h,  X2>h).\mathbb P(A_x) \le 2\mathbb P(x-h<X\le x) +\mathbb P(X_1>h,\;X_2>h).

For the first term,

P(x−h<X≤x)Fˉ(x)=Fˉ(x−h)−Fˉ(x)Fˉ(x)=(xx−h)α−1→0,\frac{\mathbb P(x-h<X\le x)}{\bar F(x)} = \frac{\bar F(x-h)-\bar F(x)}{\bar F(x)} = \left(\frac{x}{x-h}\right)^\alpha-1 \to0,

because h/x→0h/x\to0. For the second term, independence gives

P(X1>h,  X2>h)Fˉ(x)=Fˉ(h)2Fˉ(x)=(x/xm)α(1−2γ)→0,\frac{\mathbb P(X_1>h,\;X_2>h)}{\bar F(x)} = \frac{\bar F(h)^2}{\bar F(x)} = (x/x_m)^{\alpha(1-2\gamma)} \to0,

because γ>1/2\gamma>1/2. Thus P(Ax)=o(Fˉ(x))\mathbb P(A_x)=o(\bar F(x)). Since {S2>x}\{S_2>x\} is the disjoint union of {M2>x}\{M_2>x\} and AxA_x,

P(X1+X2>x)∼2Fˉ(x).\mathbb P(X_1+X_2>x)\sim 2\bar F(x).

The regularly varying theorem cited above generalizes this Pareto calculation: the maximum formula is universal, while the negligible residual event is the heavy-tail property.

A useful comparison table is:

DistributionSubexponential?Regularly varying?
ParetoYesYes
LognormalYesNo
Weibull with shape 0<β<10<\beta<1YesNo
ExponentialNoNo

The Tail Class Catalog gives a broader comparison. The lognormal and stretched Weibull entries are cited examples, not consequences of the Pareto proof.

The exponential is an exact nonexample. If Xi∼Exp⁡(λ)X_i\sim\operatorname{Exp}(\lambda), then the universal maximum formula still holds:

P(M2>x)=2e−λx−e−2λx∼2e−λx.\mathbb P(M_2>x)=2e^{-\lambda x}-e^{-2\lambda x} \sim2e^{-\lambda x}.

But the sum tail is much larger:

P(X1+X2>x)=e−λx(1+λx),\mathbb P(X_1+X_2>x)=e^{-\lambda x}(1+\lambda x),

so

P(X1+X2>x)2P(X>x)=1+λx2→∞,\frac{\mathbb P(X_1+X_2>x)}{2\mathbb P(X>x)} =\frac{1+\lambda x}{2}\to\infty,

not 1.

The fixed-nn assumption is load-bearing. If nn grows with xx, additional large-deviation regimes can appear.

A static bound on the residual tail

The proof supplies a finite-threshold bound. Put t=x/xmt=x/x_m and take α=1.5\alpha=1.5, γ=0.75\gamma=0.75. For sufficiently large tt with x−h≥xmx-h\ge x_m,

0≤P(Ax)Fˉ(x)≤2[(1−t−1/4)−3/2−1]+t−3/4.0\le\frac{\mathbb P(A_x)}{\bar F(x)} \le 2\left[(1-t^{-1/4})^{-3/2}-1\right]+t^{-3/4}.
t=x/xmt=x/x_mUpper bound (rounded upward)
1040.343428
1080.030381
10120.003004

Together with P(M2>x)=2Fˉ(x)−Fˉ(x)2\mathbb P(M_2>x)=2\bar F(x)-\bar F(x)^2, this bounds the sum tail around its leading term. It illustrates how an asymptotic statement can require large thresholds even when the model is exactly Pareto. These are conservative analytic bounds, not measured convolution probabilities.

The fixed-nn, large-threshold limit differs from the growing-sample max-to-sum ratio. A finite-mean Pareto law is subexponential even though Mn/Sn→0M_n/S_n\to0 almost surely as n→∞n\to\infty.

The chapter 5 diagnostic figure evaluates R(x)R(x) by convolution for Pareto, stretched Weibull, and exponential laws. The first two illustrate the cited classification; the exponential curve is checked against the exact formula above. The quadrature has an explicit error tolerance and cannot prove a limiting statement.

Caveats

References

References
  1. Chistyakov, V. P. (1964). A Theorem on Sums of Independent Positive Random Variables and Its Applications to Branching Random Processes. Theory of Probability and Its Applications, 9(4), 640–648. 10.1137/1109088
  2. Embrechts, P., Klüppelberg, C., & Mikosch, T. (1997). Modelling Extremal Events. Springer. 10.1007/978-3-642-33483-2
  3. Foss, S., Korshunov, D., & Zachary, S. (2013). An Introduction to Heavy-Tailed and Subexponential Distributions (2nd ed.). Springer. 10.1007/978-1-4614-7101-1
  4. Taleb, N. N. (2025). Statistical Consequences of Fat Tails (3rd ed.). https://arxiv.org/abs/2001.10488v4
  5. Bingham, N. H., Goldie, C. M., & Teugels, J. L. (1987). Regular Variation. Cambridge University Press. 10.1017/CBO9780511721434
  6. Resnick, S. I. (2007). Heavy-Tail Phenomena: Probabilistic and Statistical Modeling. Springer. 10.1007/978-0-387-45024-7