Characteristic Functions

How do you prove that averages of independent messy random variables become Gaussian? Direct manipulation of densities via repeated convolution quickly becomes unworkable. Convolving nn arbitrary densities has no closed form. The fix is to transform the problem. Replace each random variable by a complex-valued function that turns convolution into multiplication. This transform, the characteristic function, is the Fourier transform of a probability distribution. It exists for every random variable, discrete or continuous, with or without a density. The characteristic function matters for two reasons. First, it is a lossless encoding of a distribution — two random variables have the same law if and only if they have the same characteristic function. Second, independence translates into the simplest possible algebraic operation: multiplication of characteristic functions. Sums of independent random variables are analytically painful to study directly. They become trivial once transformed. This single fact is the engine behind the proof of the Central Limit Theorem. It is also behind stability properties of families like the Gaussian, Cauchy, Poisson, and Gamma distributions. It underlies the entire subject of weak convergence of probability measures. The big picture is a dictionary. Probabilistic operations on distributions (convolution, scaling, independence, weak limits) correspond to simple analytic operations. These are operations on characteristic functions (multiplication, rescaling of argument, uniform continuity, pointwise convergence). Once this dictionary is built and shown to be faithful — via the Uniqueness Theorem — hard probabilistic questions become questions in complex analysis. There, tools like Taylor expansion, dominated convergence, and pointwise limits are far more tractable. This article builds that dictionary rigorously. We give the precise definition and basic properties. We prove the Uniqueness Theorem (via a Fourier inversion argument). We prove Lévy's Continuity Theorem, connecting pointwise convergence of characteristic functions to convergence in distribution. We also prove the behavior of characteristic functions under sums of independent random variables. These three results are exactly the machinery needed to state and prove the Central Limit Theorem, which we foreshadow at the end.

Key takeaways
  • The characteristic function of XX is φX(t)=E[eitX]\varphi_X(t) = \mathbb{E}[e^{itX}], defined for every random variable, for all real tt.
  • φX\varphi_X always exists (no integrability needed), is uniformly continuous, and satisfies φX(0)=1\varphi_X(0)=1, ∣φX(t)∣≤1|\varphi_X(t)|\le 1.
  • Uniqueness: φX=φY\varphi_X = \varphi_Y everywhere implies XX and YY have the same distribution.
  • Independence becomes multiplication: if X⊥YX\perp Y, then φX+Y(t)=φX(t)φY(t)\varphi_{X+Y}(t) = \varphi_X(t)\varphi_Y(t).
  • Lévy's Continuity Theorem: pointwise convergence of characteristic functions to a function continuous at 00 is equivalent to convergence in distribution.
  • Key idea: characteristic functions turn convolution (hard) into multiplication (easy), giving a faithful, always-defined transform that is the engine of limit theorems.

Where this lives today

Characteristic functions are the analytic backbone of the proof of the Central Limit Theorem: expanding φX(t/n)\varphi_{X}(t/\sqrt n) in a Taylor series around t=0t=0 and taking the nn-th power reduces the CLT to an elementary limit of the form (1+an)n→ea(1+a_n)^n \to e^a, exactly the technique used in the Poisson-limit worked example above. In statistics, characteristic functions underlie the theory of stable distributions (used to model heavy-tailed financial returns), and inversion formulas are used to recover densities numerically from empirical characteristic functions (a technique used in high-frequency finance and signal processing). In physics, the same Fourier-transform structure appears as the partition function in statistical mechanics and as wavefunctions in quantum mechanics, where the characteristic function of position is (up to constants) the momentum-space wavefunction. In cryptography and coding theory, Fourier-analytic methods over finite groups (discrete analogues of characteristic functions) underlie Fourier-based cryptanalysis and coding bounds. In computer science, characteristic-function-based tail bounds (via analytic continuation to moment generating functions) give concentration inequalities used in randomized algorithm analysis.

Ask the tutor: How exactly does the Taylor expansion of φX\varphi_X near t=0t=0 get combined with the Sum Theorem to prove the Central Limit Theorem step by step?

DefinitionCharacteristic Function

Let XX be a real-valued random variable on a probability space (Ω,F,P)(\Omega,\mathcal F,\mathbb P), with distribution (law) μX\mu_X on (R,B(R))(\mathbb R,\mathcal B(\mathbb R)). The characteristic function of XX is the function φX:R→C\varphi_X:\mathbb R\to\mathbb C defined by

φX(t)  =  E[eitX]  =  ∫Reitx dμX(x),t∈R,\varphi_X(t) \;=\; \mathbb E\big[e^{itX}\big] \;=\; \int_{\mathbb R} e^{itx}\,d\mu_X(x), \qquad t\in\mathbb R,

where eitx=cos⁡(tx)+isin⁡(tx)e^{itx} = \cos(tx) + i\sin(tx), and the integral is understood as

φX(t)=E[cos⁡(tX)]+i E[sin⁡(tX)].\varphi_X(t) = \mathbb E[\cos(tX)] + i\,\mathbb E[\sin(tX)].

For a discrete random variable with p.m.f. p(xk)p(x_k), φX(t)=∑keitxkp(xk)\varphi_X(t) = \sum_k e^{itx_k} p(x_k); for a continuous random variable with density fXf_X, φX(t)=∫−∞∞eitxfX(x) dx\varphi_X(t) = \int_{-\infty}^{\infty} e^{itx} f_X(x)\,dx. More generally, for a random vector X=(X1,…,Xn)\mathbf X = (X_1,\dots,X_n), the joint characteristic function is

φX(t1,…,tn)=E[ei(t1X1+⋯+tnXn)],(t1,…,tn)∈Rn.\varphi_{\mathbf X}(t_1,\dots,t_n) = \mathbb E\big[e^{i(t_1X_1+\cdots+t_nX_n)}\big], \qquad (t_1,\dots,t_n)\in\mathbb R^n.

Intuition

Think of the characteristic function as a spectral fingerprint of a random variable. Just as a musical chord can be decomposed into pure sine waves of various frequencies, a probability distribution can be "decomposed" against the family of pure oscillations eitxe^{itx}, indexed by frequency tt. The value φX(t)\varphi_X(t) measures how much the distribution of XX resonates with the oscillation of frequency tt: it is a weighted average, over all outcomes, of the phase eitXe^{itX}. Because ∣eitx∣=1|e^{itx}|=1 always, no integrability condition on XX is needed — the expectation E[cos⁡(tX)]+iE[sin⁡(tX)]\mathbb E[\cos(tX)] + i\mathbb E[\sin(tX)] is a bounded quantity for every random variable, unlike the moment generating function, which can fail to exist (e.g. for heavy-tailed distributions). This universality is the first great advantage of characteristic functions over moment generating functions. The second advantage is algebraic: convolution of densities corresponds to multiplication of characteristic functions. If XX and YY are independent, the density of X+YX+Y is the convolution fX∗fYf_X * f_Y, a genuinely complicated integral operation. But φX+Y(t)=φX(t)φY(t)\varphi_{X+Y}(t) = \varphi_X(t)\varphi_Y(t) — ordinary multiplication of two complex numbers. This is exactly analogous to how logarithms turn multiplication into addition, except here the transform turns convolution into multiplication.

Figure
The characteristic function as a Fourier "fingerprint": real and imaginary parts of φ for an asymmetric two-point distribution
t value Re φ_X(t) = 0.75cos(t)+0.25cos(2t) Im φ_X(t) = 0.75sin(t)−0.25sin(2t) φ(0)=1

Basic properties (established before the main theorems)

LemmaElementary properties of φX\varphi_X

For any random variable XX, the characteristic function φX\varphi_X satisfies: 1. φX(0)=1\varphi_X(0) = 1. 2. ∣φX(t)∣≤1|\varphi_X(t)| \le 1 for all t∈Rt \in \mathbb R. 3. φX\varphi_X is uniformly continuous on R\mathbb R. 4. φX(−t)=φX(t)‾\varphi_X(-t) = \overline{\varphi_X(t)} (the complex conjugate). 5. For constants a,b∈Ra,b\in\mathbb R, φaX+b(t)=eibtφX(at)\varphi_{aX+b}(t) = e^{ibt}\varphi_X(at). 6. If XX has an nn-th moment, φX\varphi_X is nn times differentiable with φX(k)(0)=ikE[Xk]\varphi_X^{(k)}(0) = i^k \mathbb E[X^k] for k≤nk\le n.

Proof

(1) φX(0)=E[ei0⋅X]=E[1]=1\varphi_X(0) = \mathbb E[e^{i0\cdot X}] = \mathbb E[1] = 1.

(2) By the triangle inequality for expectations (Jensen's inequality applied to the convex function ∣⋅∣|\cdot|, or directly ∣E[Z]∣≤E[∣Z∣]|\mathbb E[Z]|\le \mathbb E[|Z|] for complex ZZ),

∣φX(t)∣=∣E[eitX]∣≤E[∣eitX∣]=E[1]=1,|\varphi_X(t)| = \big|\mathbb E[e^{itX}]\big| \le \mathbb E\big[|e^{itX}|\big] = \mathbb E[1] = 1,

since ∣eitX∣=1|e^{itX}| = 1 pointwise.

(3) Fix t,h∈Rt,h\in\mathbb R. Then

∣φX(t+h)−φX(t)∣=∣E[eitX(eihX−1)]∣≤E[∣eihX−1∣],|\varphi_X(t+h) - \varphi_X(t)| = \big|\mathbb E\big[e^{itX}(e^{ihX}-1)\big]\big| \le \mathbb E\big[|e^{ihX}-1|\big],

using ∣eitX∣=1|e^{itX}|=1 and the triangle inequality. The bound E[∣eihX−1∣]\mathbb E[|e^{ihX}-1|] does not depend on tt. The integrand ∣eihX−1∣|e^{ihX}-1| is bounded by 22 and tends to 00 pointwise as h→0h\to 0 (for every fixed outcome X=xX=x, eihx→1e^{ihx}\to 1 by continuity of t↦eitxt\mapsto e^{itx}). By the Dominated Convergence Theorem (dominating function the constant 22, integrable since E[2]=2<∞\mathbb E[2]=2<\infty),

E[∣eihX−1∣]→h→00.\mathbb E\big[|e^{ihX}-1|\big] \xrightarrow[h\to 0]{} 0.

Since the bound is uniform in tt, this shows uniform continuity: for every ε>0\varepsilon>0 there is δ>0\delta>0 such that ∣h∣<δ|h|<\delta implies sup⁡t∣φX(t+h)−φX(t)∣<ε\sup_t |\varphi_X(t+h)-\varphi_X(t)|<\varepsilon.

(4) φX(−t)=E[e−itX]=E[eitX‾]=E[eitX]‾=φX(t)‾\varphi_X(-t) = \mathbb E[e^{-itX}] = \mathbb E[\overline{e^{itX}}] = \overline{\mathbb E[e^{itX}]} = \overline{\varphi_X(t)}, using that conjugation commutes with expectation (linearity applied to real and imaginary parts).

(5) φaX+b(t)=E[eit(aX+b)]=eibt E[ei(at)X]=eibtφX(at)\varphi_{aX+b}(t) = \mathbb E[e^{it(aX+b)}] = e^{ibt}\,\mathbb E[e^{i(at)X}] = e^{ibt}\varphi_X(at).

(6) Suppose E[∣X∣n]<∞\mathbb E[|X|^n]<\infty. We show φX\varphi_X is nn times differentiable by induction, differentiating under the expectation sign. For k≤nk \le n, define gk(t)=E[(iX)keitX]g_k(t) = \mathbb E[(iX)^k e^{itX}]; we claim φX(k)(t)=gk(t)\varphi_X^{(k)}(t) = g_k(t). The base case k=0k=0 is the definition. For the inductive step, suppose φX(k−1)(t)=E[(iX)k−1eitX]\varphi_X^{(k-1)}(t) = \mathbb E[(iX)^{k-1}e^{itX}] for some k≤nk \le n. Consider the difference quotient

φX(k−1)(t+h)−φX(k−1)(t)h=E[(iX)k−1eitX⋅eihX−1h].\frac{\varphi_X^{(k-1)}(t+h) - \varphi_X^{(k-1)}(t)}{h} = \mathbb E\left[(iX)^{k-1} e^{itX}\cdot \frac{e^{ihX}-1}{h}\right].

Using the elementary inequality ∣eiu−1u∣≤1\left|\frac{e^{iu}-1}{u}\right| \le 1 for all real u≠0u\ne 0 (which follows from ∣eiu−1∣=2∣sin⁡(u/2)∣≤∣u∣|e^{iu}-1| = 2|\sin(u/2)| \le |u|), we get

∣(iX)k−1eitX⋅eihX−1h∣≤∣X∣k−1⋅∣X∣=∣X∣k,\left|(iX)^{k-1} e^{itX}\cdot \frac{e^{ihX}-1}{h}\right| \le |X|^{k-1}\cdot |X| = |X|^{k},

which is integrable since k≤nk \le n and E[∣X∣n]<∞\mathbb E[|X|^n]<\infty implies E[∣X∣k]<∞\mathbb E[|X|^k]<\infty for k≤nk\le n (by Lyapunov/monotonicity of moments on a probability space). As h→0h\to 0, eihX−1h→iX\frac{e^{ihX}-1}{h}\to iX pointwise. By Dominated Convergence,

lim⁡h→0φX(k−1)(t+h)−φX(k−1)(t)h=E[(iX)k−1⋅iX eitX]=E[(iX)keitX]=gk(t),\lim_{h\to 0}\frac{\varphi_X^{(k-1)}(t+h)-\varphi_X^{(k-1)}(t)}{h} = \mathbb E[(iX)^{k-1}\cdot iX\, e^{itX}] = \mathbb E[(iX)^k e^{itX}] = g_k(t),

completing the induction. Setting t=0t=0: φX(k)(0)=E[(iX)k]=ikE[Xk]\varphi_X^{(k)}(0) = \mathbb E[(iX)^k] = i^k\mathbb E[X^k]. ■\blacksquare

∎

The central theorems

We now state and prove the three pillars: Uniqueness, behavior under sums, and Lévy's Continuity Theorem.

TheoremUniqueness Theorem

Let XX and YY be random variables with characteristic functions φX\varphi_X and φY\varphi_Y. If

φX(t)=φY(t)for all t∈R,\varphi_X(t) = \varphi_Y(t) \quad \text{for all } t \in \mathbb R,

then XX and YY have the same distribution, i.e. P(X∈B)=P(Y∈B)\mathbb P(X \in B) = \mathbb P(Y\in B) for every Borel set BB.

Proof

We use the Fourier inversion formula, proved as a lemma first, and then deduce uniqueness as a corollary. The technique is direct computation and interchange of integrals (justified by Fubini), followed by a limiting/dominated convergence argument.

Step 1 (Inversion formula for intervals). We claim that for any a<ba<b real numbers that are continuity points of μX\mu_X (i.e. μX({a})=μX({b})=0\mu_X(\{a\})=\mu_X(\{b\})=0),

μX((a,b])=lim⁡T→∞12π∫−TTe−ita−e−itbit φX(t) dt.\mu_X\big((a,b]\big) = \lim_{T\to\infty} \frac{1}{2\pi}\int_{-T}^{T} \frac{e^{-ita}-e^{-itb}}{it}\,\varphi_X(t)\,dt.

Fix T>0T>0 and define

IT:=12π∫−TTe−ita−e−itbitφX(t) dt.I_T := \frac{1}{2\pi}\int_{-T}^{T}\frac{e^{-ita}-e^{-itb}}{it}\varphi_X(t)\,dt.

Substituting φX(t)=∫Reitx dμX(x)\varphi_X(t) = \int_{\mathbb R} e^{itx}\,d\mu_X(x) gives

IT=12π∫−TT∫Reit(x−a)−eit(x−b)it dμX(x) dt.I_T = \frac{1}{2\pi}\int_{-T}^{T}\int_{\mathbb R}\frac{e^{it(x-a)}-e^{it(x-b)}}{it}\,d\mu_X(x)\,dt.

The integrand, as a function of (t,x)(t,x), is bounded: since ∣eiu−eivi∣≤∣u−v∣\left|\frac{e^{iu}-e^{iv}}{i}\right| \le |u-v|... more precisely we bound each piece using ∣eitu−1it∣≤∣u∣\left|\frac{e^{itu}-1}{it}\right|\le |u| (shown in the Lemma's part 6 computation), so the integrand is bounded uniformly in t∈[−T,T]t\in[-T,T] by (b−a)(b-a), and is a jointly measurable function of (t,x)(t,x) on the finite-measure space [−T,T]×R[-T,T]\times \mathbb R (finite in tt, probability measure in xx). Fubini's Theorem applies, so we may swap the order of integration:

IT=∫R[12π∫−TTeit(x−a)−eit(x−b)it dt]dμX(x).I_T = \int_{\mathbb R} \left[\frac{1}{2\pi}\int_{-T}^{T} \frac{e^{it(x-a)}-e^{it(x-b)}}{it}\,dt\right] d\mu_X(x).

Define, for real uu, JT(u):=12π∫−TTeitu−1it dt⋅(−1)J_T(u) := \frac{1}{2\pi}\int_{-T}^{T}\frac{e^{itu}-1}{it}\,dt \cdot (-1)... Let us instead directly compute the bracketed inner integral. Write

eit(x−a)−eit(x−b)it=∫x−bx−a...\frac{e^{it(x-a)}-e^{it(x-b)}}{it} = \int_{x-b}^{x-a}...

It is cleaner to use the classical Dirichlet integral representation: define

DT(u):=12π∫−TTsin⁡(tu)t dt=1π∫0Tsin⁡(tu)t dt⋅sgn⁡(u)-adjusted, D_T(u) := \frac{1}{2\pi}\int_{-T}^{T} \frac{\sin(tu)}{t}\,dt = \frac{1}{\pi}\int_0^T \frac{\sin(tu)}{t}\,dt \cdot \operatorname{sgn}(u)\text{-adjusted, }

and use the known limit (Dirichlet integral):

lim⁡T→∞∫0Tsin⁡(tu)t dt=π2sgn⁡(u),u≠0,\lim_{T\to\infty} \int_0^T \frac{\sin(tu)}{t}\,dt = \frac{\pi}{2}\operatorname{sgn}(u), \qquad u\ne 0,

with the integral equal to 00 when u=0u=0; moreover ∣∫0Tsin⁡(tu)tdt∣\left|\int_0^T \frac{\sin(tu)}{t}dt\right| is bounded uniformly in T,uT,u by a universal constant MM (a standard fact about the Dirichlet integral). Since eitu−1it\frac{e^{itu}-1}{it} has imaginary part −cos⁡(tu)−1t-\frac{\cos(tu)-1}{t}, an odd function of tt integrating to 00 over [−T,T][-T,T], and real part sin⁡(tu)t\frac{\sin(tu)}{t}, an even function, one computes directly

12π∫−TTeitu−1it dt=1π∫0Tsin⁡(tu)t dt.\frac{1}{2\pi}\int_{-T}^{T} \frac{e^{itu}-1}{it}\,dt = \frac{1}{\pi}\int_{0}^{T}\frac{\sin(tu)}{t}\,dt.

Applying this with u=x−au = x-a and u=x−bu=x-b (noting eit(x−a)−eit(x−b)it=eit(x−a)−1it−eit(x−b)−1it\frac{e^{it(x-a)}-e^{it(x-b)}}{it} = \frac{e^{it(x-a)}-1}{it} - \frac{e^{it(x-b)}-1}{it}) gives

12π∫−TTeit(x−a)−eit(x−b)it dt=1π∫0Tsin⁡(t(x−a))t dt−1π∫0Tsin⁡(t(x−b))t dt.\frac{1}{2\pi}\int_{-T}^{T}\frac{e^{it(x-a)}-e^{it(x-b)}}{it}\,dt = \frac{1}{\pi}\int_0^T \frac{\sin(t(x-a))}{t}\,dt - \frac{1}{\pi}\int_0^T \frac{\sin(t(x-b))}{t}\,dt.

As T→∞T\to\infty, using the Dirichlet limit, this tends to

12sgn⁡(x−a)−12sgn⁡(x−b)=:h(x),\frac{1}{2}\operatorname{sgn}(x-a) - \frac{1}{2}\operatorname{sgn}(x-b) =: h(x),

where h(x)=1h(x) = 1 for a<x<ba<x<b, h(x)=0h(x)=0 for x<ax<a or x>bx>b, and h(x)=1/2h(x) = 1/2 at x=ax=a or x=bx=b. The whole quantity is uniformly bounded in (T,x)(T,x) by the universal Dirichlet bound 2M/π2M/\pi. By the Bounded (Dominated) Convergence Theorem applied to the probability measure μX\mu_X,

lim⁡T→∞IT=∫Rh(x) dμX(x)=μX((a,b))+12μX({a})+12μX({b}).\lim_{T\to\infty} I_T = \int_{\mathbb R} h(x)\, d\mu_X(x) = \mu_X((a,b)) + \tfrac12\mu_X(\{a\}) + \tfrac12\mu_X(\{b\}).

If a,ba,b are continuity points, μX({a})=μX({b})=0\mu_X(\{a\}) = \mu_X(\{b\}) = 0, so this equals μX((a,b))=μX((a,b])\mu_X((a,b)) = \mu_X((a,b]). This proves the inversion formula.

Step 2 (Uniqueness follows). Suppose φX(t)=φY(t)\varphi_X(t) = \varphi_Y(t) for all tt. By Step 1, for all continuity points a<ba<b of both μX\mu_X and μY\mu_Y (all but at most countably many reals are continuity points of a given measure, since a measure has at most countably many atoms; the union of two countable atom-sets is still countable, so continuity points of both are co-countable and in particular dense),

μX((a,b])=lim⁡T→∞12π∫−TTe−ita−e−itbitφX(t) dt=lim⁡T→∞12π∫−TTe−ita−e−itbitφY(t) dt=μY((a,b]),\mu_X((a,b]) = \lim_{T\to\infty} \frac{1}{2\pi}\int_{-T}^T \frac{e^{-ita}-e^{-itb}}{it}\varphi_X(t)\,dt = \lim_{T\to\infty}\frac{1}{2\pi}\int_{-T}^{T}\frac{e^{-ita}-e^{-itb}}{it}\varphi_Y(t)\,dt = \mu_Y((a,b]),

since φX≡φY\varphi_X\equiv \varphi_Y. Thus μX\mu_X and μY\mu_Y agree on all intervals (a,b](a,b] with a,ba,b common continuity points; by right-continuity of the CDFs and density of continuity points, μX\mu_X and μY\mu_Y agree on all intervals (a,b](a,b], hence (since these generate the Borel σ\sigma-algebra and form a π\pi-system) μX=μY\mu_X = \mu_Y on all of B(R)\mathcal B(\mathbb R) by the standard π\pi-λ\lambda uniqueness argument for measures agreeing on a generating π\pi-system with total mass 1<∞1<\infty. Hence XX and YY have the same distribution. ■\blacksquare

∎
TheoremCharacteristic Function of a Sum of Independent Random Variables

Let X1,…,XnX_1,\dots,X_n be independent random variables (not necessarily identically distributed) with characteristic functions φX1,…,φXn\varphi_{X_1},\dots,\varphi_{X_n}. Let Sn=X1+⋯+XnS_n = X_1+\cdots+X_n. Then for all t∈Rt\in\mathbb R,

φSn(t)=∏k=1nφXk(t).\varphi_{S_n}(t) = \prod_{k=1}^{n} \varphi_{X_k}(t).
Proof

The technique is direct expansion via independence and induction on nn.

Base case n=2n=2. For independent X1,X2X_1, X_2,

φS2(t)=E[eit(X1+X2)]=E[eitX1eitX2].\varphi_{S_2}(t) = \mathbb E[e^{it(X_1+X_2)}] = \mathbb E[e^{itX_1}e^{itX_2}].

Write eitX1=cos⁡(tX1)+isin⁡(tX1)e^{itX_1} = \cos(tX_1)+i\sin(tX_1) and similarly for X2X_2; the product eitX1eitX2e^{itX_1}e^{itX_2} is a (complex) linear combination of the four real products cos⁡(tX1)cos⁡(tX2)\cos(tX_1)\cos(tX_2), cos⁡(tX1)sin⁡(tX2)\cos(tX_1)\sin(tX_2), sin⁡(tX1)cos⁡(tX2)\sin(tX_1)\cos(tX_2), sin⁡(tX1)sin⁡(tX2)\sin(tX_1)\sin(tX_2). Since X1⊥X2X_1 \perp X_2, any bounded measurable function of X1X_1 is independent of any bounded measurable function of X2X_2 (a standard consequence of the definition of independence of random variables, established when studying Expectation of Functions of Several Random Variables), and expectation of a product of independent random variables factors:

E[cos⁡(tX1)cos⁡(tX2)]=E[cos⁡(tX1)] E[cos⁡(tX2)],\mathbb E[\cos(tX_1)\cos(tX_2)] = \mathbb E[\cos(tX_1)]\,\mathbb E[\cos(tX_2)],

and similarly for the other three terms (each factor is bounded, hence integrable, so the factorization identity for expectations of products of independent bounded random variables applies without extra integrability hypotheses). Recombining,

E[eitX1eitX2]=E[eitX1] E[eitX2]=φX1(t)φX2(t).\mathbb E[e^{itX_1}e^{itX_2}] = \mathbb E[e^{itX_1}]\,\mathbb E[e^{itX_2}] = \varphi_{X_1}(t)\varphi_{X_2}(t).

(This factorization can also be seen directly: E[ZW]=E[Z]E[W]\mathbb E[ZW] = \mathbb E[Z]\mathbb E[W] for independent complex random variables Z=eitX1,W=eitX2Z=e^{itX_1}, W=e^{itX_2}, since this identity for real and imaginary parts is exactly the four real identities above.)

Inductive step. Suppose the claim holds for n−1n-1 independent random variables. Given independent X1,…,XnX_1,\dots,X_n, note that Sn−1=X1+⋯+Xn−1S_{n-1} = X_1+\cdots+X_{n-1} is a function of (X1,…,Xn−1)(X_1,\dots,X_{n-1}) and is independent of XnX_n (a function of independent random variables is independent of the remaining ones — a standard grouping property of independence). By the base case applied to the pair (Sn−1,Xn)(S_{n-1}, X_n),

φSn(t)=φSn−1+Xn(t)=φSn−1(t) φXn(t).\varphi_{S_n}(t) = \varphi_{S_{n-1}+X_n}(t) = \varphi_{S_{n-1}}(t)\,\varphi_{X_n}(t).

By the inductive hypothesis, φSn−1(t)=∏k=1n−1φXk(t)\varphi_{S_{n-1}}(t) = \prod_{k=1}^{n-1}\varphi_{X_k}(t), so

φSn(t)=∏k=1n−1φXk(t)⋅φXn(t)=∏k=1nφXk(t).\varphi_{S_n}(t) = \prod_{k=1}^{n-1}\varphi_{X_k}(t)\cdot \varphi_{X_n}(t) = \prod_{k=1}^{n}\varphi_{X_k}(t).

By induction the identity holds for all n≥2n\ge 2 (and trivially for n=1n=1). ■\blacksquare

∎
Corollaryi.i.d. sums

If X1,…,XnX_1,\dots,X_n are i.i.d. with common characteristic function φ\varphi, then

φSn(t)=(φ(t))n.\varphi_{S_n}(t) = \big(\varphi(t)\big)^n.
Proof

Immediate from the theorem with φX1=⋯=φXn=φ\varphi_{X_1}=\cdots=\varphi_{X_n}=\varphi. ■\blacksquare

∎
TheoremLévy's Continuity Theorem

Let (Xn)n≥1(X_n)_{n\ge1} be a sequence of random variables with characteristic functions φn\varphi_n, and let XX be a random variable with characteristic function φ\varphi.

(Forward direction) If Xn→dXX_n \xrightarrow{d} X (convergence in distribution), then φn(t)→φ(t)\varphi_n(t) \to \varphi(t) for every t∈Rt\in\mathbb R, and the convergence is uniform on every bounded interval.

(Converse direction) Suppose φn(t)→ψ(t)\varphi_n(t) \to \psi(t) for every t∈Rt \in \mathbb R, for some function ψ\psi, and suppose ψ\psi is continuous at t=0t=0. Then ψ\psi is the characteristic function of some random variable XX, and Xn→dXX_n \xrightarrow{d} X.

Proof

Forward direction. For fixed tt, the map x↦cos⁡(tx)x \mapsto \cos(tx) and x↦sin⁡(tx)x\mapsto \sin(tx) are bounded and continuous functions of xx. By the definition of convergence in distribution (equivalently, the Portmanteau characterization: Xn→dXX_n\xrightarrow d X iff E[g(Xn)]→E[g(X)]\mathbb E[g(X_n)] \to \mathbb E[g(X)] for every bounded continuous g:R→Rg:\mathbb R\to\mathbb R),

E[cos⁡(tXn)]→E[cos⁡(tX)],E[sin⁡(tXn)]→E[sin⁡(tX)].\mathbb E[\cos(tX_n)] \to \mathbb E[\cos(tX)], \qquad \mathbb E[\sin(tX_n)]\to \mathbb E[\sin(tX)].

Hence φn(t)=E[cos⁡(tXn)]+iE[sin⁡(tXn)]→E[cos⁡(tX)]+iE[sin⁡(tX)]=φ(t)\varphi_n(t) = \mathbb E[\cos(tX_n)]+i\mathbb E[\sin(tX_n)] \to \mathbb E[\cos(tX)]+i\mathbb E[\sin(tX)] = \varphi(t) for every fixed tt. (Uniformity on bounded intervals follows from a standard strengthening of Portmanteau via tightness, which we do not detail further as it is not needed for the converse direction proved next.)

Converse direction. This is the deeper half; we sketch the standard argument via tightness, citing the two classical facts it rests on. Step 1: Tightness from continuity at 00. We first show the family {Xn}\{X_n\} is tight: for every ε>0\varepsilon>0 there exists M<∞M<\infty such that P(∣Xn∣>M)<ε\mathbb P(|X_n|>M) < \varepsilon for all nn. This uses the elementary inequality, valid for any random variable YY and any u>0u>0,

P(∣Y∣>2/u)≤1u∫−uu(1−φY(t)) dt,\mathbb P(|Y| > 2/u) \le \frac{1}{u}\int_{-u}^{u} \big(1-\varphi_Y(t)\big)\,dt,

which is proved by writing the right side as 1u∫−uuE[1−cos⁡(tY)] dt=E[1u∫−uu(1−cos⁡(tY)) dt]=E[2−2sin⁡(uY)uY]\frac{1}{u}\int_{-u}^u \mathbb E[1-\cos(tY)]\,dt = \mathbb E\left[\frac{1}{u}\int_{-u}^u(1-\cos(tY))\,dt\right] = \mathbb E\left[2 - \frac{2\sin(uY)}{uY}\right] (interpreting the value at Y=0Y=0 as 00 by continuity), and observing that 2−2sin⁡(uy)uy≥2−2u∣y∣≥12 - \frac{2\sin(uy)}{uy} \ge 2 - \frac{2}{u|y|} \ge 1 whenever ∣y∣>2/u|y|>2/u (using ∣sin⁡(uy)∣≤1|\sin(uy)|\le 1), so

E[2−2sin⁡(uY)uY]≥E[(2−2u∣Y∣)1∣Y∣>2/u]≥P(∣Y∣>2/u).\mathbb E\left[2-\frac{2\sin(uY)}{uY}\right] \ge \mathbb E\left[\Big(2-\tfrac{2}{u|Y|}\Big)\mathbf 1_{|Y|>2/u}\right] \ge \mathbb P(|Y|>2/u).

Applying this to Y=XnY=X_n: since φn(t)→ψ(t)\varphi_n(t)\to \psi(t) pointwise and ∣1−φn(t)∣≤2|1-\varphi_n(t)|\le 2, by Dominated Convergence,

1u∫−uu(1−φn(t)) dt→1u∫−uu(1−ψ(t)) dtas n→∞.\frac1u\int_{-u}^u (1-\varphi_n(t))\,dt \to \frac1u \int_{-u}^u (1-\psi(t))\,dt \quad \text{as } n\to\infty.

Because ψ\psi is continuous at 00 with ψ(0)=lim⁡φn(0)=1\psi(0)=\lim \varphi_n(0) = 1, the right side tends to 00 as u→0+u\to 0^+ (the integrand 1−ψ(t)→01-\psi(t)\to 0 as t→0t\to0, so the average over [−u,u][-u,u] vanishes as u→0u\to 0). Choosing uu small makes 1u∫−uu(1−ψ(t))dt<ε/2\frac1u\int_{-u}^u(1-\psi(t))dt < \varepsilon/2, and then choosing nn large controls the finite-nn discrepancy; combined with the finitely many remaining X1,…,Xn0X_1,\dots,X_{n_0} (each individually tight, being a single random variable), this yields tightness of the whole sequence. We take this construction as establishing: sup⁡nP(∣Xn∣>2/u)≤ε\sup_n \mathbb P(|X_n| > 2/u) \le \varepsilon for uu small enough, which is tightness. Step 2: Extracting a weak limit via Helly's Selection Theorem. Tightness plus Helly's Selection Theorem (every sequence of CDFs has a subsequence converging weakly to some possibly-defective distribution function; tightness rules out mass escaping to infinity, ensuring the limit is a genuine CDF) gives a subsequence (Xnk)(X_{n_k}) and a random variable XX such that Xnk→dXX_{n_k} \xrightarrow{d} X. By the forward direction already proved, φnk(t)→φX(t)\varphi_{n_k}(t) \to \varphi_X(t) for every tt. But by hypothesis φnk(t)→ψ(t)\varphi_{n_k}(t)\to \psi(t) along the full sequence, hence along the subsequence too. By uniqueness of pointwise limits, φX(t)=ψ(t)\varphi_X(t) = \psi(t) for all tt, so ψ\psi is indeed a characteristic function, namely of XX. Step 3: The full sequence converges. We must show Xn→dXX_n \xrightarrow{d} X along the entire sequence, not just the subsequence. Suppose not: then there is a bounded continuous gg and ε>0\varepsilon>0 with ∣E[g(Xn)]−E[g(X)]∣≥ε|\mathbb E[g(X_n)] - \mathbb E[g(X)]| \ge \varepsilon along some further subsequence (nj)(n_j). But (Xnj)(X_{n_j}) is itself tight (as a subsequence of a tight sequence), so by Helly's theorem again it has a further subsequence converging weakly to some X′X'; by the forward direction, φX′=ψ=φX\varphi_{X'} = \psi = \varphi_X, so by the Uniqueness Theorem, X′XX' \stackrel d= X. This contradicts ∣E[g(Xnj)]−E[g(X)]∣≥ε|\mathbb E[g(X_{n_j})]-\mathbb E[g(X)]|\ge \varepsilon persisting along a further subsequence converging to XX in distribution. Hence no such subsequence exists, and Xn→dXX_n \xrightarrow{d} X along the full sequence. ■\blacksquare

∎

Worked examples

Example 1: A sum of independent, non-identically-distributed exponential and Bernoulli-scaled variables. Let XX be Exponential(λ=2)\mathrm{Exponential}(\lambda=2) (density fX(x)=2e−2x,x>0f_X(x) = 2e^{-2x}, x>0) and let YY be independent, taking value 33 with probability 0.60.6 and value −1-1 with probability 0.40.4. We compute φX+Y\varphi_{X+Y}. First, φX(t)=∫0∞eitx⋅2e−2x dx=2∫0∞e(it−2)x dx=22−it\varphi_X(t) = \int_0^\infty e^{itx}\cdot 2e^{-2x}\,dx = 2\int_0^\infty e^{(it-2)x}\,dx = \frac{2}{2-it}, valid since Re(it−2)=−2<0\mathrm{Re}(it-2) = -2<0 ensures convergence. So

φX(t)=22−it.\varphi_X(t) = \frac{2}{2-it}.

Second, φY(t)=0.6 e3it+0.4 e−it\varphi_Y(t) = 0.6\,e^{3it} + 0.4\,e^{-it} directly from the definition for discrete variables. By the Sum Theorem (independence of X,YX,Y),

φX+Y(t)=22−it(0.6 e3it+0.4 e−it).\varphi_{X+Y}(t) = \frac{2}{2-it}\Big(0.6\,e^{3it}+0.4\,e^{-it}\Big).

Sanity check at t=0t=0: φX(0)=2/2=1\varphi_X(0) = 2/2 = 1, φY(0)=0.6+0.4=1\varphi_Y(0) = 0.6+0.4=1, so φX+Y(0)=1⋅1=1\varphi_{X+Y}(0) = 1\cdot 1 = 1, as required by Lemma property (1).

Sanity check via moments: E[X]=1/2\mathbb E[X] = 1/2 (mean of Exponential(22)), E[Y]=0.6(3)+0.4(−1)=1.8−0.4=1.4\mathbb E[Y] = 0.6(3)+0.4(-1) = 1.8-0.4=1.4, so E[X+Y]=1.9\mathbb E[X+Y] = 1.9. Differentiating φX+Y\varphi_{X+Y} at t=0t=0 should give i⋅1.9i\cdot 1.9. Using the product rule,

φX+Y′(t)=2i(2−it)2(0.6e3it+0.4e−it)+22−it(1.8i e3it−0.4i e−it).\varphi_{X+Y}'(t) = \frac{2i}{(2-it)^2}\Big(0.6e^{3it}+0.4e^{-it}\Big) + \frac{2}{2-it}\Big(1.8i\,e^{3it} - 0.4i\,e^{-it}\Big).

At t=0t=0: first term =2i4(0.6+0.4)=2i4=0.5i= \frac{2i}{4}(0.6+0.4) = \frac{2i}{4} = 0.5i. Second term =22(1.8i−0.4i)=1⋅1.4i=1.4i= \frac{2}{2}(1.8i - 0.4i) = 1\cdot 1.4i = 1.4i. Sum: 0.5i+1.4i=1.9i=iE[X+Y]0.5i + 1.4i = 1.9i = i\mathbb E[X+Y]. ✓ This confirms property (6) of the Lemma and validates the computation.

Example 2: Using the Uniqueness Theorem to identify a distribution, and the Continuity Theorem to verify a limit. Let X1,…,XnX_1,\dots,X_n be i.i.d. Poisson(λ=3)\mathrm{Poisson}(\lambda=3) random variables (so each has mean 33), and let Sn=X1+⋯+XnS_n = X_1+\cdots+X_n. We identify the distribution of SnS_n via characteristic functions, and then verify a Poisson-limit-type statement using the Continuity Theorem — illustrating both theorems on an asymmetric (non-symmetric-around-zero) example. The characteristic function of a single Poisson(3)\mathrm{Poisson}(3) variable is

φX1(t)=∑k=0∞eitk e−33kk!=e−3∑k=0∞(3eit)kk!=e−3e3eit=exp⁡(3(eit−1)).\varphi_{X_1}(t) = \sum_{k=0}^{\infty} e^{itk}\,\frac{e^{-3}3^k}{k!} = e^{-3}\sum_{k=0}^\infty \frac{(3e^{it})^k}{k!} = e^{-3}e^{3e^{it}} = \exp\big(3(e^{it}-1)\big).

By the i.i.d. Sum Corollary,

φSn(t)=(exp⁡(3(eit−1)))n=exp⁡(3n(eit−1)).\varphi_{S_n}(t) = \Big(\exp(3(e^{it}-1))\Big)^n = \exp\big(3n(e^{it}-1)\big).

This is exactly the characteristic function of a Poisson(3n)\mathrm{Poisson}(3n) random variable (by the same computation with parameter 3n3n in place of 33). By the Uniqueness Theorem, since characteristic functions determine distributions, we conclude

Sn∼Poisson(3n),S_n \sim \mathrm{Poisson}(3n),

a fact that would otherwise require an inductive convolution argument on the p.m.f.s (summing over all ways k1+⋯+kn=kk_1+\cdots+k_n=k), here obtained in one line. Now consider a triangular array: for each nn, let Yn∼Binomial(n,pn)Y_n \sim \mathrm{Binomial}(n, p_n) with npn→6np_n \to 6 as n→∞n\to\infty (e.g. pn=6/np_n = 6/n, an asymmetric, vanishing-probability regime — the classical Poisson approximation setting). We verify Yn→dPoisson(6)Y_n \xrightarrow{d} \mathrm{Poisson}(6) via the Continuity Theorem. The characteristic function of YnY_n is

φYn(t)=(1−pn+pneit)n=(1+pn(eit−1))n.\varphi_{Y_n}(t) = (1-p_n + p_n e^{it})^n = \big(1+p_n(e^{it}-1)\big)^n.

Write pn(eit−1)=npn(eit−1)np_n(e^{it}-1) = \frac{np_n(e^{it}-1)}{n}; since npn→6np_n\to 6, for fixed tt we have npn(eit−1)→6(eit−1)n p_n(e^{it}-1) \to 6(e^{it}-1) as n→∞n\to\infty. Using the standard limit (1+ann)n→ea\left(1+\frac{a_n}{n}\right)^n \to e^{a} whenever an→aa_n\to a (a fact from complex analysis, provable by taking logarithms and using log⁡(1+z)=z+O(z2)\log(1+z) = z + O(z^2) for small ∣z∣|z|),

φYn(t)=(1+npn(eit−1)n)n⟶exp⁡(6(eit−1))for every fixed t.\varphi_{Y_n}(t) = \left(1 + \frac{np_n(e^{it}-1)}{n}\right)^n \longrightarrow \exp\big(6(e^{it}-1)\big) \quad \text{for every fixed } t.

The limit function ψ(t)=exp⁡(6(eit−1))\psi(t) = \exp(6(e^{it}-1)) is continuous everywhere (in particular at t=0t=0, where ψ(0)=1\psi(0)=1), and is exactly the characteristic function of Poisson(6)\mathrm{Poisson}(6) computed above (with 33 replaced by 66). By Lévy's Continuity Theorem, Yn→dPoisson(6)Y_n \xrightarrow{d} \mathrm{Poisson}(6). This is the Poisson limit theorem (Poisson approximation to the Binomial), derived purely from pointwise convergence of characteristic functions — no direct manipulation of binomial coefficients was needed.

CounterexampleWhy continuity of the limit function at 00 is needed in Lévy's Continuity Theorem

The converse direction of Lévy's Continuity Theorem requires that the pointwise limit ψ(t)=lim⁡nφn(t)\psi(t) = \lim_n \varphi_n(t) be continuous at t=0t=0. Without this hypothesis, the conclusion fails: pointwise convergence of characteristic functions can occur even when mass escapes to infinity, so that no probability distribution is being approximated at all. Let Xn∼Uniform(−n,n)X_n \sim \mathrm{Uniform}(-n, n) for n=1,2,3,…n=1,2,3,\dots. Its characteristic function is

φn(t)=12n∫−nneitx dx=sin⁡(nt)nt,t≠0,φn(0)=1.\varphi_n(t) = \frac{1}{2n}\int_{-n}^{n} e^{itx}\,dx = \frac{\sin(nt)}{nt}, \qquad t \ne 0, \qquad \varphi_n(0)=1.

Fix t≠0t\ne 0. As n→∞n\to\infty, ∣sin⁡(nt)nt∣≤1n∣t∣→0\left|\frac{\sin(nt)}{nt}\right| \le \frac{1}{n|t|} \to 0. So

φn(t)→ψ(t):={0,t≠01,t=0.\varphi_n(t) \to \psi(t) := \begin{cases} 0, & t\ne 0\\ 1, & t=0.\end{cases}

This limit ψ\psi exists for every t∈Rt\in\mathbb R — pointwise convergence genuinely holds. But ψ\psi is discontinuous at t=0t=0 (ψ(0)=1\psi(0)=1 while lim⁡t→0ψ(t)=0\lim_{t\to0}\psi(t) = 0). Correspondingly, (Xn)(X_n) does not converge in distribution to any random variable: the mass of Uniform(−n,n)\mathrm{Uniform}(-n,n) spreads out over an interval of growing length 2n2n, so for any fixed MM, P(∣Xn∣≤M)=Mn→0\mathbb P(|X_n|\le M) = \frac{M}{n}\to 0 — all the probability mass escapes to infinity (the family {Xn}\{X_n\} is not tight). There is no limiting random variable, defective or otherwise, in the usual sense; the "limit" ψ\psi is not the characteristic function of any distribution on R\mathbb R (a genuine characteristic function is always continuous, by the Lemma above — but ψ\psi is discontinuous at 00, hence cannot be one). This shows precisely why the theorem's hypothesis "ψ\psi continuous at 00" cannot be dropped: it is exactly the condition that rules out escape of mass to infinity (equivalently, tightness), which the tightness argument in the proof used the continuity at 00 to establish.

Scope, edge cases & common traps

SituationWhat holdsWhat can fail
XX has no finite mean/variance (e.g. Cauchy)φX\varphi_X still exists, uniformly continuous, everywhere-definedMoments φX(k)(0)=ikE[Xk]\varphi_X^{(k)}(0)=i^k\mathbb E[X^k] do not exist for k≥1k\ge1; e.g. Cauchy has φ(t)=e−∣t∣\varphi(t)=e^{-|t|}, not differentiable at 00
X1,…,XnX_1,\dots,X_n dependentφSn(t)\varphi_{S_n}(t) is not generally ∏φXk(t)\prod\varphi_{X_k}(t)The Sum Theorem's multiplicativity strictly requires independence
Convergence in distribution Xn→XX_n\to Xφn(t)→φ(t)\varphi_n(t)\to\varphi(t) pointwise, uniformly on compactsConverse needs ψ\psi continuous at 00 — see counterexample
Two distinct distributionsCan share some but not all momentsUniqueness Theorem requires equality of φ\varphi at every tt, not just matching moments (moment problem can be indeterminate)
MGF existsφX(t)=MX(it)\varphi_X(t) = M_X(it) by analytic continuationMGF may not exist at all (heavy tails); φX\varphi_X always does

The mathematician behind it

Origins · 1 minPaul Lévy and how this was discovered
PL
Paul Lévy
1886–1971 · Paris, École Polytechnique

Paul Lévy was a French mathematician whose work founded much of the modern theory of probability, including The Central Limit Theorem, Lévy processes, stable distributions, and the continuity theorem for characteristic functions that bears his name. Working largely in the 1920s–1930s at the École Polytechnique, Lévy built probability theory into a rigorous branch of analysis by systematically exploiting Fourier-analytic tools — characteristic functions chief among them — to convert delicate limiting arguments about sums of random variables into tractable analytic statements about complex-valued functions.

History
- 1810: Laplace uses proto-Fourier/generating-function methods in early limit theorems for sums of random variables.
- 1922: Lindeberg refines conditions for the Central Limit Theorem using direct (non-characteristic-function) methods.
- 1925: Paul Lévy publishes foundational work establishing the continuity theorem linking convergence of characteristic functions to convergence in distribution.
- 1930s: Lévy and Khinchin develop the Lévy–Khinchin representation for infinitely divisible distributions, built entirely on characteristic function machinery.
- 1949: Cramér's textbook cements characteristic functions as the standard tool for proving central limit theorems in mathematical statistics.
Why use characteristic functions instead of moment generating functions?

Moment generating functions MX(t)=E[etX]M_X(t) = \mathbb E[e^{tX}] can fail to exist for any t≠0t\ne 0 — for example, heavy-tailed distributions like the Cauchy or certain Pareto distributions have E[etX]=∞\mathbb E[e^{tX}]=\infty for all t≠0t\ne0. Characteristic functions replace the real exponential etXe^{tX} with the complex exponential eitXe^{itX}, which has modulus exactly 11, so E[eitX]\mathbb E[e^{itX}] is always a well-defined, bounded complex number for every random variable and every real tt. This universality is what lets characteristic functions serve as the basis for the Uniqueness Theorem and Continuity Theorem for arbitrary distributions, not just light-tailed ones.

Does equality of all moments imply equal distributions?

Not necessarily. Matching moments E[Xk]=E[Yk]\mathbb E[X^k]=\mathbb E[Y^k] for all kk does not always imply XX and YY have the same distribution — this is the classical "moment problem," and there exist distinct distributions (e.g. certain log-normal-type constructions) sharing all moments. However, if the characteristic functions agree at every real tt, the Uniqueness Theorem guarantees the distributions are identical, with no such pathology. Characteristic functions are strictly more informative than the moment sequence because they encode the full distribution via an integral transform, not merely a countable list of numbers.

What does it mean that "independence becomes multiplication"?

If XX and YY are independent random variables, the density (or p.m.f.) of X+YX+Y is the convolution of the individual densities — a complicated integral (or sum) over all splittings of the total. Characteristic functions convert this operation into ordinary multiplication: φX+Y(t)=φX(t)φY(t)\varphi_{X+Y}(t) = \varphi_X(t)\varphi_Y(t). This is analogous to how logarithms convert multiplication into addition; here, the Fourier transform converts convolution into multiplication. It is the single most useful property of characteristic functions for studying sums of independent random variables, including in the proof of the Central Limit Theorem.

Is the characteristic function unique to real-valued random variables, or does it generalize?

The same construction generalizes immediately to random vectors: the joint characteristic function φX(t1,…,tn)=E[ei(t1X1+⋯+tnXn)]\varphi_{\mathbf X}(t_1,\dots,t_n) = \mathbb E[e^{i(t_1X_1+\cdots+t_nX_n)}] characterizes the joint distribution of (X1,…,Xn)(X_1,\dots,X_n), with an analogous uniqueness theorem. It also generalizes to random variables valued in more general groups (Fourier analysis on locally compact abelian groups), which underlies applications in cryptography and coding theory. The real-line case treated here is the foundational instance from which these generalizations are built.

Why is uniform continuity of φX\varphi_X useful, beyond just "niceness"?

Uniform continuity guarantees that φX\varphi_X cannot oscillate arbitrarily fast anywhere on R\mathbb R, which is essential for the inversion formula's limiting arguments and for controlling error terms in the tightness argument used in Lévy's Continuity Theorem. It also rules out certain pathological "limit functions": since every genuine characteristic function is uniformly continuous, a pointwise limit ψ\psi that fails to be continuous (as in the escaping-mass counterexample) is immediately disqualified from being a characteristic function at all.

How does this connect to the Central Limit Theorem?

The standard proof of the Central Limit Theorem takes i.i.d. random variables X1,…,XnX_1,\dots,X_n with mean 00 and variance σ2\sigma^2, forms Sn/(σn)S_n/(\sigma\sqrt n), and shows its characteristic function converges pointwise to e−t2/2e^{-t^2/2}, the characteristic function of the standard normal. This uses exactly the machinery built here: the Sum Theorem to write φSn/(σn)(t)=(φX1(t/(σn)))n\varphi_{S_n/(\sigma\sqrt n)}(t) = \big(\varphi_{X_1}(t/(\sigma\sqrt n))\big)^n, a second-order Taylor expansion of φX1\varphi_{X_1} near 00 (using the moment property of the Lemma), and then Lévy's Continuity Theorem to conclude convergence in distribution to N(0,1)N(0,1).

Can two different characteristic functions ever agree on an interval but not everywhere?

Yes — agreement on a proper subset of R\mathbb R does not force full agreement, and does not imply equal distributions; the Uniqueness Theorem strictly requires equality of φX(t)\varphi_X(t) and φY(t)\varphi_Y(t) for every real tt. This is a common trap: checking equality "for small tt" or on a bounded range is insufficient, since characteristic functions are generally not determined by their values on a sub-interval (they need not be real-analytic globally, though they often are for specific families).

Before this

Where this leads