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 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.
- The characteristic function of is , defined for every random variable, for all real .
- always exists (no integrability needed), is uniformly continuous, and satisfies , .
- Uniqueness: everywhere implies and have the same distribution.
- Independence becomes multiplication: if , then .
- Lévy's Continuity Theorem: pointwise convergence of characteristic functions to a function continuous at 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 in a Taylor series around and taking the -th power reduces the CLT to an elementary limit of the form , 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 near get combined with the Sum Theorem to prove the Central Limit Theorem step by step?
Let be a real-valued random variable on a probability space , with distribution (law) on . The characteristic function of is the function defined by
where , and the integral is understood as
For a discrete random variable with p.m.f. , ; for a continuous random variable with density , . More generally, for a random vector , the joint characteristic function is
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 , indexed by frequency . The value measures how much the distribution of resonates with the oscillation of frequency : it is a weighted average, over all outcomes, of the phase . Because always, no integrability condition on is needed — the expectation 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 and are independent, the density of is the convolution , a genuinely complicated integral operation. But — 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.
Basic properties (established before the main theorems)
For any random variable , the characteristic function satisfies: 1. . 2. for all . 3. is uniformly continuous on . 4. (the complex conjugate). 5. For constants , . 6. If has an -th moment, is times differentiable with for .
(1) .
(2) By the triangle inequality for expectations (Jensen's inequality applied to the convex function , or directly for complex ),
since pointwise.
(3) Fix . Then
using and the triangle inequality. The bound does not depend on . The integrand is bounded by and tends to pointwise as (for every fixed outcome , by continuity of ). By the Dominated Convergence Theorem (dominating function the constant , integrable since ),
Since the bound is uniform in , this shows uniform continuity: for every there is such that implies .
(4) , using that conjugation commutes with expectation (linearity applied to real and imaginary parts).
(5) .
(6) Suppose . We show is times differentiable by induction, differentiating under the expectation sign. For , define ; we claim . The base case is the definition. For the inductive step, suppose for some . Consider the difference quotient
Using the elementary inequality for all real (which follows from ), we get
which is integrable since and implies for (by Lyapunov/monotonicity of moments on a probability space). As , pointwise. By Dominated Convergence,
completing the induction. Setting : .
∎The central theorems
We now state and prove the three pillars: Uniqueness, behavior under sums, and Lévy's Continuity Theorem.
Let and be random variables with characteristic functions and . If
then and have the same distribution, i.e. for every Borel set .
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 real numbers that are continuity points of (i.e. ),
Fix and define
Substituting gives
The integrand, as a function of , is bounded: since ... more precisely we bound each piece using (shown in the Lemma's part 6 computation), so the integrand is bounded uniformly in by , and is a jointly measurable function of on the finite-measure space (finite in , probability measure in ). Fubini's Theorem applies, so we may swap the order of integration:
Define, for real , ... Let us instead directly compute the bracketed inner integral. Write
It is cleaner to use the classical Dirichlet integral representation: define
and use the known limit (Dirichlet integral):
with the integral equal to when ; moreover is bounded uniformly in by a universal constant (a standard fact about the Dirichlet integral). Since has imaginary part , an odd function of integrating to over , and real part , an even function, one computes directly
Applying this with and (noting ) gives
As , using the Dirichlet limit, this tends to
where for , for or , and at or . The whole quantity is uniformly bounded in by the universal Dirichlet bound . By the Bounded (Dominated) Convergence Theorem applied to the probability measure ,
If are continuity points, , so this equals . This proves the inversion formula.
Step 2 (Uniqueness follows). Suppose for all . By Step 1, for all continuity points of both and (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),
since . Thus and agree on all intervals with common continuity points; by right-continuity of the CDFs and density of continuity points, and agree on all intervals , hence (since these generate the Borel -algebra and form a -system) on all of by the standard - uniqueness argument for measures agreeing on a generating -system with total mass . Hence and have the same distribution.
∎Let be independent random variables (not necessarily identically distributed) with characteristic functions . Let . Then for all ,
The technique is direct expansion via independence and induction on .
Base case . For independent ,
Write and similarly for ; the product is a (complex) linear combination of the four real products , , , . Since , any bounded measurable function of is independent of any bounded measurable function of (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:
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,
(This factorization can also be seen directly: for independent complex random variables , since this identity for real and imaginary parts is exactly the four real identities above.)
Inductive step. Suppose the claim holds for independent random variables. Given independent , note that is a function of and is independent of (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 ,
By the inductive hypothesis, , so
By induction the identity holds for all (and trivially for ).
∎If are i.i.d. with common characteristic function , then
Immediate from the theorem with .
∎Let be a sequence of random variables with characteristic functions , and let be a random variable with characteristic function .
(Forward direction) If (convergence in distribution), then for every , and the convergence is uniform on every bounded interval.
(Converse direction) Suppose for every , for some function , and suppose is continuous at . Then is the characteristic function of some random variable , and .
Forward direction. For fixed , the map and are bounded and continuous functions of . By the definition of convergence in distribution (equivalently, the Portmanteau characterization: iff for every bounded continuous ),
Hence for every fixed . (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 . We first show the family is tight: for every there exists such that for all . This uses the elementary inequality, valid for any random variable and any ,
which is proved by writing the right side as (interpreting the value at as by continuity), and observing that whenever (using ), so
Applying this to : since pointwise and , by Dominated Convergence,
Because is continuous at with , the right side tends to as (the integrand as , so the average over vanishes as ). Choosing small makes , and then choosing large controls the finite- discrepancy; combined with the finitely many remaining (each individually tight, being a single random variable), this yields tightness of the whole sequence. We take this construction as establishing: for 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 and a random variable such that . By the forward direction already proved, for every . But by hypothesis along the full sequence, hence along the subsequence too. By uniqueness of pointwise limits, for all , so is indeed a characteristic function, namely of . Step 3: The full sequence converges. We must show along the entire sequence, not just the subsequence. Suppose not: then there is a bounded continuous and with along some further subsequence . But 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 ; by the forward direction, , so by the Uniqueness Theorem, . This contradicts persisting along a further subsequence converging to in distribution. Hence no such subsequence exists, and along the full sequence.
∎Worked examples
Example 1: A sum of independent, non-identically-distributed exponential and Bernoulli-scaled variables. Let be (density ) and let be independent, taking value with probability and value with probability . We compute . First, , valid since ensures convergence. So
Second, directly from the definition for discrete variables. By the Sum Theorem (independence of ),
Sanity check at : , , so , as required by Lemma property (1).
Sanity check via moments: (mean of Exponential()), , so . Differentiating at should give . Using the product rule,
At : first term . Second term . Sum: . ✓ 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 be i.i.d. random variables (so each has mean ), and let . We identify the distribution of 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 variable is
By the i.i.d. Sum Corollary,
This is exactly the characteristic function of a random variable (by the same computation with parameter in place of ). By the Uniqueness Theorem, since characteristic functions determine distributions, we conclude
a fact that would otherwise require an inductive convolution argument on the p.m.f.s (summing over all ways ), here obtained in one line. Now consider a triangular array: for each , let with as (e.g. , an asymmetric, vanishing-probability regime — the classical Poisson approximation setting). We verify via the Continuity Theorem. The characteristic function of is
Write ; since , for fixed we have as . Using the standard limit whenever (a fact from complex analysis, provable by taking logarithms and using for small ),
The limit function is continuous everywhere (in particular at , where ), and is exactly the characteristic function of computed above (with replaced by ). By Lévy's Continuity Theorem, . 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.
The converse direction of Lévy's Continuity Theorem requires that the pointwise limit be continuous at . 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 for . Its characteristic function is
Fix . As , . So
This limit exists for every — pointwise convergence genuinely holds. But is discontinuous at ( while ). Correspondingly, does not converge in distribution to any random variable: the mass of spreads out over an interval of growing length , so for any fixed , — all the probability mass escapes to infinity (the family is not tight). There is no limiting random variable, defective or otherwise, in the usual sense; the "limit" is not the characteristic function of any distribution on (a genuine characteristic function is always continuous, by the Lemma above — but is discontinuous at , hence cannot be one). This shows precisely why the theorem's hypothesis " continuous at " 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 to establish.
Scope, edge cases & common traps
| Situation | What holds | What can fail |
|---|---|---|
| has no finite mean/variance (e.g. Cauchy) | still exists, uniformly continuous, everywhere-defined | Moments do not exist for ; e.g. Cauchy has , not differentiable at |
| dependent | is not generally | The Sum Theorem's multiplicativity strictly requires independence |
| Convergence in distribution | pointwise, uniformly on compacts | Converse needs continuous at — see counterexample |
| Two distinct distributions | Can share some but not all moments | Uniqueness Theorem requires equality of at every , not just matching moments (moment problem can be indeterminate) |
| MGF exists | by analytic continuation | MGF may not exist at all (heavy tails); always does |
The mathematician behind it
Origins · 1 minPaul Lévy and how this was discovered
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.
Why use characteristic functions instead of moment generating functions?
Moment generating functions can fail to exist for any — for example, heavy-tailed distributions like the Cauchy or certain Pareto distributions have for all . Characteristic functions replace the real exponential with the complex exponential , which has modulus exactly , so is always a well-defined, bounded complex number for every random variable and every real . 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 for all does not always imply and 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 , 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 and are independent random variables, the density (or p.m.f.) of 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: . 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 characterizes the joint distribution of , 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 useful, beyond just "niceness"?
Uniform continuity guarantees that cannot oscillate arbitrarily fast anywhere on , 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 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 with mean and variance , forms , and shows its characteristic function converges pointwise to , the characteristic function of the standard normal. This uses exactly the machinery built here: the Sum Theorem to write , a second-order Taylor expansion of near (using the moment property of the Lemma), and then Lévy's Continuity Theorem to conclude convergence in distribution to .
Can two different characteristic functions ever agree on an interval but not everywhere?
Yes — agreement on a proper subset of does not force full agreement, and does not imply equal distributions; the Uniqueness Theorem strictly requires equality of and for every real . This is a common trap: checking equality "for small " 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).