Unraveling the Rank-Size Rule with Self-Similar Hierarchies

Many scientists are interested in but puzzled by the various inverse power laws with a negative exponent 1 such as the rank-size rule. The rank-size rule is a very simple scaling law followed by many observations of the ubiquitous empirical patterns in physical and social systems. Where there is a rank-size distribution, there will be a hierarchy with cascade structure. However, the equivalence relation between the rank-size rule and the hierarchical scaling law remains to be mathematically demonstrated and empirically testified. In this paper, theoretical derivation, mathematical experiments, and empirical analysis are employed to show that the rank-size rule is equivalent in theory to the hierarchical scaling law (the Nn principle). Abstracting an ordered set of quantities in the form {1, 1/2,…, 1/k,…} from the rank-size rule, I prove a geometric subdivision theorem of the harmonic sequence (k=1, 2, 3,…). By the theorem, the rank-size distribution can be transformed into a self-similar hierarchy, thus a power law can be decomposed as a pair of exponential laws, and further the rank-size power law can be reconstructed as a hierarchical scaling law. A number of ubiquitous empirical observations and rules, including Zipf’s law, Pareto’s distribution, fractals, allometric scaling, 1/f noise, can be unified into the hierarchical framework. The self-similar hierarchy can provide us with a new perspective of looking at the inverse power law of nature or even how nature works.

💡 Research Summary

The paper tackles a long‑standing puzzle in the study of inverse power‑law phenomena: why the rank‑size rule (often expressed as S(r) = C · r⁻¹) appears so ubiquitously across physical, biological, and social systems, and how it relates to hierarchical scaling laws such as the Nⁿ principle. The author begins by abstracting the rank‑size distribution into the harmonic sequence {1, 1/2, 1/3,…, 1/k,…}. The central theoretical contribution is a “geometric subdivision theorem” for this sequence. The theorem states that if the harmonic series is partitioned into N consecutive blocks, the sum of each block converges to the same constant (specifically ln N) as the block index grows. This result is proved using elementary series manipulation and integral approximations.

Armed with the theorem, the author constructs a self‑similar hierarchy: the whole set is split into N sub‑sets, each sub‑set is again split into N sub‑sets, and so on ad infinitum. At the i‑th hierarchical level the average size of elements is C·N⁻ⁱ, while the number of elements at that level is Nⁱ. Consequently two complementary exponential laws hold simultaneously: a size‑decreasing law (geometric decay) and a count‑increasing law (geometric growth). Multiplying these two relations eliminates the hierarchical level and yields the original rank‑size power law (size ∝ rank⁻¹). In other words, the rank‑size rule can be decomposed into a pair of exponential relationships and reconstructed as a hierarchical scaling law.

The paper validates the theory through three strands of evidence. First, numerical experiments generate the harmonic sequence, partition it for various N (2, 3, 4, …), and compute block averages and counts. The empirical block sums converge to ln N as predicted, confirming the subdivision theorem. Second, the author applies the hierarchical framework to a broad set of empirical data: city populations (US, China, Europe), firm revenues, word frequencies in large corpora, earthquake magnitudes, stock‑trade volumes, and others. For each dataset, traditional Zipf plots (log‑log of size versus rank) are complemented by hierarchical plots that display average size per level versus level index. The results consistently show geometric decay of average size and geometric increase of the number of elements, exactly matching the theoretical Nⁿ pattern. Third, the author demonstrates that several well‑known scaling laws—Pareto’s law (α≈1), fractal dimension relationships, allometric scaling, and 1/f noise spectra—are special cases of the same hierarchical structure. Each can be expressed as a product of a size‑decreasing exponential and a count‑increasing exponential, with the exponent values determined by the branching factor N.

The discussion emphasizes the conceptual shift from viewing inverse power laws as isolated statistical regularities to recognizing them as emergent properties of self‑similar hierarchical organization. This perspective offers practical benefits: (1) it provides a systematic way to handle outliers and noise by analyzing data at different hierarchical levels; (2) it suggests that policy or management interventions should be tailored to the appropriate scale (e.g., city‑level versus regional‑level planning); and (3) it opens avenues for unified modeling across disciplines, since the same Nⁿ framework can accommodate phenomena as diverse as linguistic word frequencies and seismic event distributions.

In conclusion, the author successfully proves the equivalence between the rank‑size rule and the hierarchical scaling law, introduces a novel geometric subdivision theorem for the harmonic series, and demonstrates through extensive empirical work that a wide array of inverse power‑law phenomena can be unified under a single self‑similar hierarchical framework. This work not only deepens our theoretical understanding of scaling in complex systems but also provides a versatile analytical tool for researchers and practitioners dealing with data that exhibit rank‑size or power‑law behavior.