Random Matrix Theory

the moment method and non-crossing pairings

The moment method is the first and most hands-on proof that the eigenvalue histogram of a Wigner matrix is a semicircle. Its idea is borrowed straight from Volume I's method of moments for ordinary convergence in distribution: a measure on the line (with light tails) is determined by its sequence of moments, so to show a random measure converges to a target, you show every moment converges to the target's moment. For random matrices the moments of the eigenvalue distribution are computable without ever finding an eigenvalue — through the trace.

The key identity is that the k-th moment of the empirical spectral distribution is a normalised trace: the integral of x^k against mu_{H_N} equals (1/N) sum over i of lambda_i^k = (1/N) trace(H_N^k). Now expand the trace of the k-th power as a sum over closed sequences of indices i_1 -> i_2 -> ... -> i_k -> i_1: trace(H_N^k) = sum over such closed walks of the product H_{i_1 i_2} H_{i_2 i_3} ... H_{i_k i_1}. Taking the expectation, every walk contributes E of a product of entries. Because the entries are independent with mean zero, a walk contributes nothing unless every edge it uses is traversed an even number of times — each entry must be paired with itself to survive. Counting which pairings survive the N -> infinity scaling is then a pure combinatorics problem: when k = 2m is even, the dominant contributions come from walks that are encoded by non-crossing pair partitions of the 2m steps, walks that trace out a tree (a planar structure). The number of such non-crossing pairings of 2m objects is exactly the Catalan number C_m, and the odd moments vanish — which is precisely the moment sequence of the semicircle.

Non-crossing pairings are the heart of the matter. Draw the 2m steps as points on a line and connect paired steps by arcs above the line: a pairing is non-crossing if no two arcs intersect. Non-crossing pairings of 2m points number C_m = (1/(m+1)) binomial(2m, m), the Catalan numbers, and they are in bijection with planar trees, balanced parenthesisations, and Dyck paths. The crossing pairings, by contrast, force the closed walk to revisit fewer distinct indices and so come with an extra power of 1/N that kills them in the limit. This combinatorial dichotomy — non-crossing survives, crossing dies — is the same structure that reappears in free probability, where freeness is defined by exactly the vanishing of crossing (alternating) moments. A caveat: the moment method needs all moments of the entries to be finite (or careful truncation), and it proves convergence of the bulk density but is poorly suited to the edge, where the largest eigenvalue requires moments of growing order k ~ N^(2/3) and far more delicate counting.

Compute E[(1/N) trace(H_N^4)] by hand. A closed walk of length 4 on the index set is i->j->k->l->i. For its expectation to survive, every entry must appear an even number of times, forcing the walk to use exactly two distinct edges each traversed twice. The surviving walks correspond to the two non-crossing pairings of 4 steps (matching 1-2,3-4 or 1-4,2-3), each contributing 1; the single crossing pairing 1-3,2-4 forces a longer-range constraint and is suppressed by 1/N. Total: 2 = C_2, the fourth semicircle moment.

Surviving walks = non-crossing pairings; their count is the Catalan number, the semicircle moment.

The method works for the bulk because non-crossing pairings dominate; crossings carry an extra 1/N. It is the wrong tool for the edge, where the largest eigenvalue needs trace(H_N^k) with k growing like N^(2/3) and the simple counting breaks down.

Also called
combinatorial methodmethod of moments for random matricestrace method動差法跡法