Branching Processes
A single bacterium divides, or fails to. Its offspring divide or fail, generation after generation. Will the family line survive forever, or is extinction certain? This question is simple to state, but subtle to answer. It is the subject of branching processes. Its consequences reach far beyond biology: the fate of surnames, the spread of an epidemic, the growth of a nuclear chain reaction, and the survival of a mutant gene all obey the same mathematical skeleton. The Galton–Watson process is named for the two Victorians who first studied the extinction of family surnames. It models a population in which each individual independently produces a random number of offspring, according to a fixed distribution. The remarkable fact is this: the entire long-run behavior of this process is governed by a single number. That number is the mean number of offspring per individual. It determines whether the population dies out with probability one, or survives with positive probability. This is a strikingly clean dichotomy. It is proved using an equally clean tool: the probability generating function of the offspring distribution, iterated once per generation. The big picture is that branching processes sit at the crossroads of two ideas the reader already knows. They are Markov chains, since the population size at generation depends on the past only through the population size at generation . But their state space is unbounded. Their transition structure is generated by a sum of i.i.d. random variables. This is exactly where generating functions become indispensable. They turn a probabilistic recursion into an algebraic one (iteration of a fixed function). This article builds the Galton–Watson process from its probabilistic definition. It develops the generating-function machinery that makes it tractable. It proves the central extinction theorem in full: extinction is certain if and only if the mean offspring number is at most 1. This excludes the degenerate case of no variability. Along the way, we compute extinction probabilities exactly in two asymmetric examples. We exhibit the boundary (critical) case where the naive intuition breaks. Finally, we point toward the martingale viewpoint that unifies and extends everything here.
- A Galton–Watson process models a population where each individual independently has a random number of children, i.i.d. across all individuals and generations.
- Population sizes form a Markov chain on with an absorbing state.
- The probability generating function of the offspring law controls everything: the pgf of is the -fold iterate .
- The extinction probability is the smallest fixed point in of .
- Extinction Theorem: if the mean offspring (and is not a.s. equal to ) then ; if then .
- Key idea: iterating a single generating function replaces an intractable probabilistic recursion with a tractable algebraic fixed-point problem.
Where this lives today
Branching processes are a genuine working tool, not just a historical curiosity. In epidemiology, the basic reproduction number of an infectious disease plays exactly the role of : if an outbreak dies out with certainty (in the idealized branching approximation), while allows a positive-probability major epidemic, with the same extinction-probability equation used to estimate the chance a single imported case fizzles out. In nuclear physics, the same mathematics (originally analyzed via Monte Carlo methods by Ulam and von Neumann) governs whether a chain reaction of neutron-induced fissions is sub-critical, critical, or super-critical — the terminology is inherited directly. In genetics, the survival probability of a new mutant allele in a large population is computed via the extinction probability of an associated branching process (Haldane's classical approximation near criticality for slightly advantageous mutants uses exactly this framework). In computer science, branching processes model the size of search trees (e.g., in randomized algorithms and parsing), and the analysis of random recursive trees and random graphs (e.g. component sizes in Erdős–Rényi random graphs near the phase transition) reduces to Galton–Watson extinction/survival dichotomies. The natural next step in the theory — recasting as a martingale — gives sharp almost-sure growth-rate theorems (the Kesten–Stigum theorem) refining everything proved here.
Ask the tutor: Can you walk through why the extinction probability equals when starting from ancestors instead of ?
Let be independent and identically distributed random variables taking values in , with common probability mass function , . Assume (offspring are not deterministically zero) to avoid triviality. Define a stochastic process by and, for ,
where for each is an independent family of i.i.d. copies of , independent across as well, and the empty sum () is . Interpretation: is the population size at generation ; each of the individuals alive at generation independently produces -many offspring, and is the total. The process is called a Galton–Watson process with offspring distribution . Extinction is the event
and the extinction probability is .
Intuition
Picture the process as a tree growing downward. The root is the single ancestor at generation . Each node sprouts a random number of children, drawn independently from the offspring law, and the process continues layer by layer. If at some generation every individual has zero children, the tree stops growing there and the population is extinct forever, because is an absorbing state: once , every subsequent automatically, since summing zero copies of anything gives zero. The key quantity is the mean offspring number . Intuitively, (proved below), so if the expected population shrinks geometrically to , if it stays constant at , and if it grows geometrically. One might guess naively that forces extinction and prevents it — and the first half of that guess is right, but the second needs care: even when the population can still die out, just not with certainty. The precise threshold is captured not by the mean directly, but by a fixed point of the generating function . The picture below shows the graph of for a fixed offspring law, together with the "staircase" (cobweb) construction that generates the iterates , converging to the smallest root of in — precisely the extinction probability.
The central theorem
Let be a Galton–Watson process with offspring pgf for , and suppose (non-trivial offspring) and (the offspring law is not degenerate at , ruling out a trivial process that never changes size). Let
and let . Then: 1. is the smallest root in of the fixed-point equation . 2. If , then (extinction is certain). 3. If , then (extinction is not certain); more precisely is the unique root of in .
Let be the pgf of . Then for all ,
where denotes the -fold functional composition of with itself ().
(Lemma: generating function iteration.) We use the technique of conditioning on and applying independence of offspring. Condition on . Given this event, is a sum of i.i.d. copies of , independent of the history up to time (this is exactly the defining recursion). By the standard product rule for pgfs of sums of independent random variables,
Hence, by the law of total expectation,
This proves . We now show separately that , i.e. that composition in the other order also works, by induction on . For : since , so . Assume for some . Then using the recursion just derived twice,
and separately
by the inductive hypothesis applied inside. Since composition of functions is associative, trivially (both denote applying then , versus then... to avoid confusion we verify directly): write . By the base recursion applied with the whole process shifted, one instead argues more directly as follows. A cleaner induction: define as above and prove (the -fold self-composition) directly by induction on , which immediately gives both composition identities since for iterates of a single fixed function (iterates of one function always commute with themselves). Base case : . Inductive step: assuming , the relation proved above gives
This completes the induction and proves for all , from which both stated identities follow immediately, the second equality because self-iterates of one function commute with one more application of that function on either side.
∎(Theorem, part 1: is the smallest fixed point of .) We use the technique of monotone limits together with continuity of pgfs. Let . Since is impossible (), . The event implies , because is absorbing (if generation is already extinct, generation is too, as the empty sum is ). Hence is an increasing sequence of events in , so is non-decreasing, and
by continuity of probability measures along increasing unions. By the Lemma, , i.e. is obtained by starting at and applying times: , with . Since is a power series with non-negative coefficients summing to , it is continuous on (indeed real-analytic on and continuous up to by Abel's theorem for power series with non-negative coefficients and finite value at ). Taking in and using continuity of ,
So is a fixed point of in . To see is the smallest such fixed point, let be any fixed point, . We show for all by induction. Base case: . Inductive step: if , then since is non-decreasing on (all coefficients , so for ),
By induction for all , hence . So is every fixed point in , i.e. is the smallest one. This proves part 1.
(Theorem, part 2 and 3: the role of .) We use the technique of convexity analysis of on . Since all and , differentiating termwise (valid inside the radius of convergence, which is at least since ) gives
So is convex and non-decreasing on , and by monotone convergence, as (allowing if the offspring distribution has infinite mean). Also , so is always a fixed point of . We must determine whether there is a smaller fixed point in , and the answer is governed by the slope of at . Case . We claim for all , which forces the smallest fixed point in to be itself (if throughout , no point there can be a fixed point, so for finite but the increasing sequence can only converge to the fixed point , by part 1's characterization — more directly, since is a fixed point of by part 1, and no point in is a fixed point in this case, we must have ). To prove the claim, define on . Then and . Since is non-decreasing (as ) and , we get for all , hence on , so is non-increasing on . Since and is non-increasing, for , i.e. for . We need the strict inequality to rule out other fixed points. Suppose for contradiction for some . Since is non-increasing and , must be identically on (a non-increasing function that returns to its starting value must be constant on the interval). Then for all , i.e. identically there. Since is non-decreasing and equals on this whole subinterval while , and is analytic (a power series), forces constant on all of by analytic continuation (a power series agreeing with a constant on an interval has all higher coefficients zero). But ; if this is identically then and for , i.e. , meaning almost surely — excluded by our hypothesis (this hypothesis says , incompatible with ; note also rules out the other degenerate case separately handled). Contradiction. Hence no such exists, proving strictly on , and therefore . Case . Now , so since is continuous and non-decreasing near with limit , there exists such that for , hence there, so is strictly increasing just to the left of . Since , strict increase just before forces , i.e. , for slightly less than . On the other hand, (and by hypothesis, with the boundary case handled below). By the Intermediate Value Theorem applied to the continuous function on , since and for some close to , there exists a point with , i.e. . (If exactly, note already gives a fixed point at itself, and the argument below shows it is not the smallest fixed point strictly less than is still attained since we need the smallest root; we handle uniqueness next which resolves this cleanly regardless.) We now show this fixed point in is unique, using strict convexity. Since , is not the identity, so (as argued above) has at least one term for some or the mean could not exceed the coefficient alone in the required way — more directly: , and if on then is affine, with , giving , contradicting . So for some, hence (since it's a power series with non-negative coefficients) for all in some interval , actually for all once any coefficient for is active, making strictly convex on . A strictly convex function on with can cross the line at most twice on (a strictly convex function minus a linear function, , is itself strictly convex, and a strictly convex function has at most two zeros). We have exhibited one crossing at and one crossing at ; strict convexity of forbids a third, so is the unique fixed point in . Finally we confirm using part 1 (smallest fixed point) together with monotonicity. Since is the smallest fixed point in by part 1, and the only fixed points are and , with , we conclude . This proves parts 2 and 3.
∎Worked examples
Example 1 — Asymmetric geometric offspring, supercritical. Let the offspring distribution be geometric-shifted: for , with . So an individual has children with probability : The pgf is
The mean is
so by the theorem this is a supercritical process: extinction is not certain, and is the unique root in of . Solve , i.e. , i.e.
Factor: , giving roots and . The smallest root is . So starting from a single ancestor, this population goes extinct with probability and survives forever (population size ) with probability , despite each individual expecting children on average. This illustrates that only makes survival possible, not certain — extinction can still be the more likely outcome.
Example 2 — Asymmetric finite support, subcritical, exact extinction probability at every generation. Let take values in with , , (so ; note the support need not be consecutive integers). The pgf is
The mean is . This is subcritical: by the theorem , extinction is certain regardless of any finer detail of the distribution. We can still extract quantitative information before extinction is reached. Compute . Compute . Compute . These values increase monotonically toward , as guaranteed by part 1 of the theorem (), illustrating the cobweb picture concretely: each new value is applied to the previous one, climbing the staircase toward the fixed point at . Also, exponentially, consistent with — but a strictly weaker fact than — almost-sure extinction; the mean going to zero does not by itself prove in general unbounded settings, though here the full theorem gives the stronger conclusion directly.
The theorem assumes , i.e. , precisely to exclude the degenerate offspring law (that is, , all other ). If instead exactly, then identically for all , so , and every is a fixed point of . The theorem's conclusion "" would be false here: since every individual has exactly one child forever, for all , so the population never goes extinct, giving , not . Tracing through the proof shows exactly where it breaks: in the case , we needed to be strictly negative-sloped-enough to force strictly on ; but here everywhere, so for all (since is literally the identity map, never moves from ), giving . The "smallest fixed point" is still correctly identified as by part 1 of the theorem (which needs no non-degeneracy assumption), but part 2's claim that this smallest fixed point equals specifically requires ruling out the totally-degenerate offspring law, which is exactly the hypothesis dropped. This shows the hypothesis is not a technicality: without it, the mean-based dichotomy statement literally fails.
Scope, edge cases & common traps
| Situation | What happens | Note |
|---|---|---|
| No individual ever has zero children | always if , so automatically (theorem still consistent: is itself the fixed point when typically, or check directly) | |
| exactly (critical case), | Extinction still certain, | But for all — the mean does not go to zero, yet extinction is still a.s.; a classic place naive "mean " intuition misleads |
| (infinite mean offspring) | Theorem's conclusion still applies whenever in the extended sense | Proof goes through since we only used , and trivially |
| Multiple generations, non-i.i.d. environment | Not covered here | This is the "branching process in varying/random environment," a genuine generalization |
| ancestors instead of | Extinction probability becomes | Because the initial lines evolve independently, each with individual extinction probability |
The mathematician behind it
Origins · 1 minFrancis Galton and how this was discovered
Francis Galton, a cousin of Charles Darwin, was a prolific Victorian scientist working across meteorology, psychology, and statistics — a founder of regression analysis and correlation, and a controversial pioneer (and namer) of eugenics. In 1873 he posed the "problem of the extinction of families" in an educational journal, worried about the disappearance of aristocratic surnames, and the mathematician Henry William Watson supplied a generating-function solution the following year. Their joint 1875 paper is the origin of the process now bearing both their names, though it contained a subtle error about the extinction probability at that was only fully corrected decades later.
Why is the smallest fixed point of the extinction probability, rather than any other fixed point?
Because the extinction probability is literally defined as a limit of an increasing sequence starting at . The proof shows this sequence is monotone increasing (since extinction by generation implies extinction by generation ) and bounded above by , so it converges to some fixed point of . A separate monotonicity argument (using that is non-decreasing) shows this limit cannot exceed any other fixed point , because propagates inductively through the recursion . Hence the limit is the smallest fixed point available, not an arbitrary one. Larger fixed points (like , which is always a fixed point) correspond to non-extinction scenarios and are simply not reached by this particular increasing sequence starting from .
Does guarantee the population grows to infinity if it doesn't go extinct?
Yes, in a precise sense, though proving it fully uses the martingale machinery beyond this article's scope. On the survival event (probability when ), one can show almost surely: a Galton–Watson process either dies out or its size diverges, it cannot hover at a positive finite value forever, because the chain has no other absorbing or recurrent states besides — every positive integer is transient. The martingale converges almost surely to a non-negative limit , and on the survival event typically (Kesten–Stigum theorem, under a mild extra moment condition), giving the precise exponential growth rate .
What happens exactly at the critical case ?
Extinction is still certain (), by part 2 of the theorem, provided the offspring law is not degenerate at . However the rate of extinction is qualitatively different and much slower than in the subcritical case: instead of decaying geometrically, a classical result (Kolmogorov) shows as , where , i.e. only a polynomial decay. This is the borderline, most delicate regime, and its analysis genuinely requires more refined tools (Tauberian theorems or direct recursive estimates on ) beyond the basic fixed-point argument given here.
Can the offspring distribution have infinite mean, and does the theorem still apply?
Yes. If , we interpret using the monotone limit definition of used in the proof (which never assumed finiteness, only that is non-decreasing with some limit in as ). Since trivially, part 3 of the theorem applies verbatim: , and is still the unique fixed point of in . The convexity argument used to prove uniqueness never required to be finite, only that strictly, so the whole proof survives unchanged.
How is this related to the Markov chain machinery I already know?
A Galton–Watson process is a special Markov chain on the countable state space , but with two features that make the ordinary finite-chain toolkit (eigenvalues of a finite transition matrix, stationary distributions) inapplicable: the state space is infinite, and the transition mechanism is a sum of i.i.d. offspring rather than an arbitrary stochastic matrix. This special branching structure is exactly what allows generating functions to substitute for the transition-matrix algebra: instead of matrix powers , we use functional iterates , which are far easier to analyze because composition of one-dimensional functions is more tractable than matrix multiplication in infinite dimensions.
Why does the article foreshadow martingales — what's the connection?
The quantity satisfies exactly (a direct computation from the branching recursion and ), which is precisely the defining property of a martingale. This reformulation is more powerful than the generating-function approach alone: martingale convergence theorems immediately give almost-sure convergence of to a limiting random variable , refining the crude dichotomy "extinction vs. survival" into a precise a.s. growth rate on the survival event. The generating-function extinction theorem proved here is the natural first result; the martingale viewpoint is the natural next chapter.