Tutorials › Graph Theory (Undergraduate) › Critical Scaling Near Mean Degree One

Random Graphs · Tutorial 768 of 799

Critical Scaling Near Mean Degree One

Learn how the critical-window parameter controls exploration drift and why component sizes in this regime are on the \(n^{2/3}\) scale.

Advanced 12 min read

What You'll Learn

  • Define the critical window for \(G(n,p)\) and distinguish \(np\) from the expected degree.
  • Translate between the critical-window parameter and an additive perturbation of \(1/n\).
  • Derive the exploration-drift balance that points to the \(n^{2/3}\) scale.
  • State the standard critical random graph scaling limit and its assumptions.
  • Compare the critical regime with fixed subcritical and fixed supercritical probabilities.

Why Mean Degree One Has a Window

In a random graph, a small change in edge probability can change the typical sizes of connected components. The value \(p=1/n\) is the boundary between two familiar behaviors: below it, components tend to be small; above it, a component containing a positive fraction of the vertices appears. The transition is not confined to the single value \(1/n\). There is a window of probabilities around it where the largest components have an intermediate scale.

Throughout, \(G(n,p)\) is the random simple graph on \(n\) labelled vertices in which each possible edge is included independently with probability \(p\). A component is a maximal connected subgraph, and \(L_1(G)\) denotes the number of vertices in a largest component. The expected degree of any vertex is \((n-1)p\), while \(np\) is a convenient approximation used to locate the transition.

Definition: For a fixed real number \(\lambda\), the critical window is the sequence of random graph models \(G(n,p_n)\) with \(p_n=(1+\lambda n^{-1/3})/n\). Here \(\lambda\) is fixed as \(n\) grows. For every fixed \(\lambda\), this is a valid probability for all sufficiently large \(n\).

The window has width \(n^{-4/3}\) in \(p\): rewriting the definition gives \(p_n=1/n+\lambda n^{-4/3}\). Equivalently, \(np_n-1=\lambda n^{-1/3}\). Thus the mean-degree excess is of order \(n^{-1/3}\), not a fixed positive or negative constant.

Worked Example: The Centre of the Window

Take \(n=1000\) and \(\lambda=0\). Since \(1000^{-1/3}=0.1\), the formula gives \(p=(1+0\cdot0.1)/1000=1/1000=0.001\). Then \(np=1000\cdot0.001=1\), whereas the expected degree is \((n-1)p=999\cdot0.001=0.999\).

The scale \(n^{2/3}\) at this value of \(n\) is \(1000^{2/3}=(1000^{1/3})^2=10^2=100\). The calculation does not say that the largest component has exactly 100 vertices; it identifies the scale on which its size fluctuates.

Check: Substituting \(p=0.001\) back into the model gives \(np=1\) and \((n-1)p=0.999\), as computed. The distinction between these two quantities is a finite-\(n\) correction.

How Exploration Produces the \(n^{2/3}\) Scale

Explore a component by maintaining a queue of discovered but not yet explored vertices. When a vertex is processed, reveal its edges to vertices that have not yet been discovered. If there are \(U_k\) undiscovered vertices at step \(k\), the number \(X_k\) newly discovered at that step has conditional distribution \(\operatorname{Bin}(U_k,p_n)\): the relevant edges have not previously been revealed. The queue changes by \(X_k-1\), because one vertex is removed for processing and \(X_k\) are added.

Suppose \(k\) vertices have been processed and the queue currently contains \(A_k\) vertices. Then \(U_k=n-k-A_k\). During a component exploration, conditional on the history so far, the expected queue change at the next step is therefore

$$ \mathbb{E}[X_k-1\mid\text{history}] =p_n(n-k-A_k)-1 =\lambda n^{-1/3}-p_n k-p_n A_k. $$

This identity explains the competing effects in the critical window. The term \(\lambda n^{-1/3}\) is the initial excess above the branching threshold. As vertices are explored, the supply of undiscovered vertices decreases; the term \(-p_n k\) records that depletion. The active-queue term also contributes, but it is smaller in the leading scale calculation.

Scale calculation: Over \(t\) exploration steps, the initial excess contributes approximately \(\lambda t n^{-1/3}\), while the accumulated depletion contributes approximately \(-p_n\sum_{k=0}^{t-1}k\), or \(-t^2/(2n)\) since \(p_n\) is close to \(1/n\). At \(t\) of order \(n^{2/3}\), both terms are of order \(n^{1/3}\). Fluctuations from roughly \(t\) near-critical discoveries have size of order \(\sqrt{t}=n^{1/3}\) as well. This balance points to exploration segments, and hence large components, on the \(n^{2/3}\) scale.

This is a scaling argument, not a proof of a component-size limit. In particular, the queue is random, components begin at different times, and the exploration must account for all vertices. The rigorous limiting result requires substantially more work.

Worked Example: A Positive Window Parameter

Take \(n=1000\) and \(\lambda=2\). First, \(n^{-1/3}=0.1\), so \(1+\lambda n^{-1/3}=1+2(0.1)=1.2\). Hence \(p=1.2/1000=0.0012\).

The approximate mean degree is \(np=1000(0.0012)=1.2\). The actual expected degree is \((n-1)p=999(0.0012)=1.1988\). The excess \(np-1=1.2-1=0.2\) agrees with \(\lambda n^{-1/3}=2(0.1)=0.2\). Also \(n^{2/3}=100\), as in the first example.

Check: The alternative form gives \(p=1/n+\lambda n^{-4/3}\). Here \(1/n=0.001\) and \(n^{-4/3}=n^{-1}n^{-1/3}=0.001(0.1)=0.0001\), so \(p=0.001+2(0.0001)=0.0012\), matching the direct calculation.

The Critical Random Graph Scaling Result

The scaling argument is made precise by a standard theorem in the theory of random graphs. To state its limiting object, let \(B(t)\), \(t\geq0\), be standard Brownian motion, and define \(W_\lambda(t)=B(t)+\lambda t-t^2/2\). Reflect this path at its running minimum: \(R_\lambda(t)=W_\lambda(t)-\min_{0\leq s\leq t}W_\lambda(s)\). An excursion is a maximal interval on which \(R_\lambda(t)>0\); its length is the interval's duration.

Theorem (Aldous's critical random graph scaling limit): Fix \(\lambda\in\mathbb{R}\), and form \(G(n,p_n)\) with \(p_n=(1+\lambda n^{-1/3})/n\). If the component sizes are arranged in decreasing order, then the sequence of sizes divided by \(n^{2/3}\) converges in distribution to the decreasing sequence of excursion lengths of \(R_\lambda\). In particular, \(L_1(G(n,p_n))/n^{2/3}\) converges in distribution to the largest excursion length, which is positive and finite almost surely. The proof of this deep scaling theorem is beyond this tutorial.

The theorem makes the scale statement precise: for fixed \(\lambda\), the largest component is of order \(n^{2/3}\) in the probabilistic sense, with a nontrivial random limiting size after normalization. It does not assert that every component has that size. Many components are much smaller, while the largest few are represented by the largest excursions in the limit.

The parameter \(\lambda\) affects the limiting process through its linear drift. Positive \(\lambda\) favors longer explorations, negative \(\lambda\) favors shorter ones, and \(\lambda=0\) is the centre of the window. All fixed values of \(\lambda\), however, remain in the same \(n^{-1/3}\) mean-degree range and have the same \(n^{2/3}\) component scale.

Worked Example: A Negative Window Parameter

Take \(n=1000\) and \(\lambda=-2\). Since \(n^{-1/3}=0.1\), \(1+\lambda n^{-1/3}=1+(-2)(0.1)=0.8\), and therefore \(p=0.8/1000=0.0008\). Then \(np=1000(0.0008)=0.8\), and the expected degree is \((n-1)p=999(0.0008)=0.7992\). The mean-degree offset based on \(np\) is \(np-1=0.8-1=-0.2\), which equals \(\lambda n^{-1/3}=(-2)(0.1)=-0.2\).

Although \(np<1\), this is not the fixed subcritical regime \(p=c/n\) with a constant \(c<1\): here the factor \(np=0.8\) depends on \(n\) through the fixed-window parameterization, and is within order \(n^{-1/3}\) of one. The critical scaling theorem still places the largest component on the \(n^{2/3}\) scale.

Check: The additive form gives \(p=0.001+(-2)(0.0001)=0.0008\), using \(n^{-4/3}=0.001(0.1)=0.0001\). This agrees with the original formula.

How the Three Regimes Differ

The critical window is best understood by comparing it with probabilities whose mean degree stays a fixed distance from one.

RegimeTypical probability scaleLargest-component behavior
Subcritical\(p=c/n\), fixed \(0<c<1\)Largest components have logarithmic size with high probability.
Critical window\(p=(1+\lambda n^{-1/3})/n\), fixed \(\lambda\)Largest components have scale \(n^{2/3}\).
Supercritical\(p=c/n\), fixed \(c>1\)A largest component has order \(n\); it contains a positive fraction of the vertices with high probability.

Here “with high probability” means with probability tending to one as \(n\) tends to infinity. The subcritical and supercritical descriptions are the fixed-\(c\) regimes; they should not be applied unchanged when \(c\) approaches one as \(n\) grows. The critical window resolves that approach: its width in \(p\) is \(n^{-4/3}\), and its largest components are much larger than logarithmic but smaller than linear in \(n\).

Common Pitfall: Confusing \(np\) with Expected Degree

A common error is to treat \(np\) as the exact expected degree. For \(G(n,p)\), a vertex has \(n-1\) possible neighbors, so its expected degree is \((n-1)p\), not \(np\). For example, with \(n=1000\) and \(p=0.0012\), \(np=1000(0.0012)=1.2\), but the expected degree is \(999(0.0012)=1.1988\). The difference is small here, but it matters when stating exact finite-\(n\) quantities. The critical-window parameterization is defined through \(np=1+\lambda n^{-1/3}\); the actual expected degree is \((1-1/n)(1+\lambda n^{-1/3})\).

A second error is to interpret the scale \(n^{2/3}\) as a deterministic formula for the largest component. The theorem gives a random limiting size after normalization. At finite \(n\), the largest component fluctuates, and the value of \(\lambda\) changes its limiting distribution. The scale describes the magnitude of those fluctuations, not an exact component size.

Check Your Understanding

Use the parameterization and scale argument above to answer the following.

  1. For fixed \(\lambda=3\), express \(p_n\) as \(1/n\) plus an additive term involving \(n\).
  2. For \(n=1000\) and \(\lambda=-1\), compute \(p_n\), \(np_n\), and the expected degree \((n-1)p_n\).
  3. In the exploration calculation, which two leading contributions balance when \(t\) is of order \(n^{2/3}\)?
  4. Why does the critical scaling theorem not imply that every component has order \(n^{2/3}\)?
  5. Contrast the largest-component scale in the critical window with the fixed-\(c\) subcritical and supercritical regimes.