Joint Scaling Law in Complex Systems

• Joint scaling laws are defined as principled relationships that couple observables via power-law, lognormal, or universal forms, linking their scaling exponents.
• They are applied across disciplines, unifying insights in linguistics, genomics, neural scaling, and city-size hierarchies for improved predictive theory.
• Methodologies like data collapse and closed-form optimality validate these joint laws, revealing constraints that refine theory and model precision.

A joint scaling law expresses a principled relationship wherein two or more observable quantities in complex systems exhibit interdependent scaling behavior—typically via power-law, lognormal, or universal functional forms—whose exponents or parameters are not independent but are mathematically or mechanistically linked. Joint scaling laws unify phenomena that otherwise appear as separate or marginal scaling patterns, providing deeper insight into universality, constraint, and predictive theory across linguistics, genomics, neural scaling, statistical physics, and large-scale deep learning.

1. Conceptual Foundations and Formal Definition

A joint scaling law models the multivariate distribution or coupled evolution of observables, typically via an explicit relationship between their conditional (or marginal) distributions, leading to non-independent constraints on scaling exponents or collapse of data under universal forms. The hallmark is that knowledge of one scaling exponent (or law) constrains or determines the others through analytically tractable links. This framework generalizes classical scaling—e.g., $P(k) \propto k^{-q}$—to multidimensional or functionally coupled variables, such as (city) rank and class size, word frequency and vocabulary size, or model and data size.

A prototypical formulation is:

• Determine the joint probability (or moment) structure of variables $X$, $Y$ such that:

$$f_{X|Y}(x|y) = y^{-\alpha}\psi_X\left(\frac{x}{y^{\alpha}}\right), \qquad f_{Y|X}(y|x) = x^{-\beta}\psi_Y\left(\frac{y}{x^{\beta}}\right)$$

leading to bivariate functional equations and explicit linking of exponents or functional forms (Aoyama et al., 2010).

The joint scaling law is also formalized in multivariable risk or performance curves for complex models:

$L(N_1,N_2,\ldots) \sim \mathcal{F}(N_1,N_2,\ldots)$

where the scaling of loss $L$ as a function of several system parameters (e.g., model size, data size, architectural factors, data ratio) is not separable, but governed by a higher-level constraint (Zhao et al., 28 Sep 2025, Zhang et al., 10 Feb 2026, He et al., 2024).

2. Canonical Instances Across Disciplines

Joint scaling laws have been established across sciences; key exemplars with explicit mathematical relationships include:

• Quantitative Linguistics: The frequency distribution $P(n;L)$ of word types as a function of frequency $n$ and text length $L$ obeys a scaling ansatz

$P(n;L) = \frac{1}{L\,V_L}g\left(\frac{n}{L}\right)$

where $X$0 (vocabulary size) is itself a function of $X$1, and both Zipf's law and Heaps' law emerge as special cases linked through the behavior of the scaling function $X$2 (Font-Clos et al., 2013).

• Multivariate Production/ Econophysics: The bivariate distribution of firm sales $X$3 and labor $X$4 in Japanese firms satisfies two joint scaling laws, yielding a unique lognormal joint PDF determined by the exponents $X$5, $X$6 of the conditional scaling, encoding a micro-macro equilibrium (Aoyama et al., 2010).
• Statistical Linguistics (Length–Frequency): The joint distribution of word length $X$7 and frequency $X$8 is found to be governed by a bivariate scaling form

$X$9

and the coupled exponents $Y$0, $Y$1 dictate the emergent Zipf exponent for the marginal $Y$2 with $Y$3 (Corral et al., 2019).

• City-Size Hierarchies: Zipf's law for city ranks connects to the hierarchical scaling law via two geometric/exponential intermediate relations, producing the equivalence

$Y$4

with $Y$5 cities per class $Y$6, class-average size $Y$7, and Zipf exponent $Y$8; thus, fractal dimension and distribution exponent are reciprocally tied (Chen, 2011).

• Genomic Evolution: The scaling of gene families and functional categories is coupled via the correlated duplication (recipe) model, predicting that the exponent $Y$9 of the family size distribution within functional category $$f_{X|Y}(x|y) = y^{-\alpha}\psi_X\left(\frac{x}{y^{\alpha}}\right), \qquad f_{Y|X}(y|x) = x^{-\beta}\psi_Y\left(\frac{y}{x^{\beta}}\right)$$0 satisfies

$$f_{X|Y}(x|y) = y^{-\alpha}\psi_X\left(\frac{x}{y^{\alpha}}\right), \qquad f_{Y|X}(y|x) = x^{-\beta}\psi_Y\left(\frac{y}{x^{\beta}}\right)$$1

where $$f_{X|Y}(x|y) = y^{-\alpha}\psi_X\left(\frac{x}{y^{\alpha}}\right), \qquad f_{Y|X}(y|x) = x^{-\beta}\psi_Y\left(\frac{y}{x^{\beta}}\right)$$2 is the scaling exponent for the category's size with total genome size—superlinear growth enforces flatter distributions (Grilli et al., 2011).

• MoE and Multilingual Neural Scaling: Joint loss $$f_{X|Y}(x|y) = y^{-\alpha}\psi_X\left(\frac{x}{y^{\alpha}}\right), \qquad f_{Y|X}(y|x) = x^{-\beta}\psi_Y\left(\frac{y}{x^{\beta}}\right)$$3 in LLMs or mixture-of-expert (MoE) architectures is expressed as a function of multiple coupled variables (e.g., model size, data size, number of experts, sampling ratios) with closed-form optimal configurations, universal exponents, and cross-factor dependencies (Zhao et al., 28 Sep 2025, He et al., 2024).

3. Mathematical Structure and Analytic Linkages

The core distinguishing feature is the presence of coupled or functional equations that relate the scaling of marginal and conditional distributions, often reducing the degrees of freedom relative to naive or independent power-law fits.

Examples:

• Linguistics (Zipf–Heaps Joint Law):

$$f_{X|Y}(x|y) = y^{-\alpha}\psi_X\left(\frac{x}{y^{\alpha}}\right), \qquad f_{Y|X}(y|x) = x^{-\beta}\psi_Y\left(\frac{y}{x^{\beta}}\right)$$4

with $$f_{X|Y}(x|y) = y^{-\alpha}\psi_X\left(\frac{x}{y^{\alpha}}\right), \qquad f_{Y|X}(y|x) = x^{-\beta}\psi_Y\left(\frac{y}{x^{\beta}}\right)$$5 scaling linearly in $$f_{X|Y}(x|y) = y^{-\alpha}\psi_X\left(\frac{x}{y^{\alpha}}\right), \qquad f_{Y|X}(y|x) = x^{-\beta}\psi_Y\left(\frac{y}{x^{\beta}}\right)$$6, so that the frequency and vocabulary growth exponents, and the transition between logarithmic and power-law growth, are analytically constrained by the exponents of $$f_{X|Y}(x|y) = y^{-\alpha}\psi_X\left(\frac{x}{y^{\alpha}}\right), \qquad f_{Y|X}(y|x) = x^{-\beta}\psi_Y\left(\frac{y}{x^{\beta}}\right)$$7 (Font-Clos et al., 2013).

• Bivariate Scaling and Productivity:

The double scaling law for $$f_{X|Y}(x|y) = y^{-\alpha}\psi_X\left(\frac{x}{y^{\alpha}}\right), \qquad f_{Y|X}(y|x) = x^{-\beta}\psi_Y\left(\frac{y}{x^{\beta}}\right)$$8 and $$f_{X|Y}(x|y) = y^{-\alpha}\psi_X\left(\frac{x}{y^{\alpha}}\right), \qquad f_{Y|X}(y|x) = x^{-\beta}\psi_Y\left(\frac{y}{x^{\beta}}\right)$$9 yields a bivariate lognormal with variance and correlation matrix determined by the scaling exponents, and the marginal distribution for productivity $L(N_1,N_2,\ldots) \sim \mathcal{F}(N_1,N_2,\ldots)$0 has scaling collapse determined by these indices (Aoyama et al., 2010).

• MoE Joint Law:

$L(N_1,N_2,\ldots) \sim \mathcal{F}(N_1,N_2,\ldots)$1

with optimal $L(N_1,N_2,\ldots) \sim \mathcal{F}(N_1,N_2,\ldots)$2, $L(N_1,N_2,\ldots) \sim \mathcal{F}(N_1,N_2,\ldots)$3, and $L(N_1,N_2,\ldots) \sim \mathcal{F}(N_1,N_2,\ldots)$4 determined analytically via minimization, showing nontrivial coupling across variables (Zhao et al., 28 Sep 2025).

• Multilingual Scaling:

$L(N_1,N_2,\ldots) \sim \mathcal{F}(N_1,N_2,\ldots)$5

where the loss in language family $L(N_1,N_2,\ldots) \sim \mathcal{F}(N_1,N_2,\ldots)$6 depends only on the family’s data fraction $L(N_1,N_2,\ldots) \sim \mathcal{F}(N_1,N_2,\ldots)$7 and not on joint mixture, validated via controlled sweep experiments (He et al., 2024).

4. Methodologies, Empirical Confirmations, and Universality

Joint scaling laws are empirically validated by data collapse under rescaling transformations, quantitative fits of observable exponents, and universality under domain shifts and mechanism variations:

• Data collapse: For vocabulary growth, plotting $L(N_1,N_2,\ldots) \sim \mathcal{F}(N_1,N_2,\ldots)$8 vs $L(N_1,N_2,\ldots) \sim \mathcal{F}(N_1,N_2,\ldots)$9 demonstrates single-master-curve behavior (Font-Clos et al., 2013).
• Closed-form optimality: For MoE systems, optimal configurations computed from joint laws robustly predict settings used in large-scale deployments (DeepSeek, Qwen, Kimi models) (Zhao et al., 28 Sep 2025).
• Universality: Transform invariance, mixture extensions, and transferability of exponents and functional forms are demonstrated across linearized NTK, finite-width, feature-learning, and neural scaling settings (Bi et al., 25 Sep 2025, He et al., 2024).
• Statistical physics analogy: The number of scaling indices is small, with macroscopic (aggregate) statistics emergent from the joint law, analogous to temperature and density defining equilibrium states (Aoyama et al., 2010).

5. Implications, Constraints, and Limitations

Joint scaling laws expose structural constraints in multivariate data: the measurement or modeling of one axis automatically specifies others, reducing ambiguity and increasing predictive power. This has implications for:

• Resource allocation: In neural scaling, joint exponents dictate optimal tradeoffs between compute, model size, and data (Ngo et al., 10 Oct 2025).
• Theory building: The existence of closed-form joint laws connects empirically observed scaling regimes to underlying stochastic or mechanistic models, as seen in genomics (coupled duplication/innovation), linguistics (unified statistical structure), and econophysics (micro–macro bridges) (Grilli et al., 2011, Font-Clos et al., 2013, Aoyama et al., 2010).
• Limitations: Functional forms and analytic links may break under regime changes, unmodeled coupling, nonpolynomial spectral tails, or when variables lose independence or the system departs from assumptions (e.g., very small data fractions, boundary artifacts, mis-tokenization, phase transitions in models) (Zhao et al., 28 Sep 2025, Bi et al., 25 Sep 2025, Chen, 2011).

6. Cross-Domain Extensions and Generalizations

Joint scaling law frameworks extend naturally to:

• Critical phenomena: Analogies are drawn with finite-size scaling and universality in statistical physics, where joint laws dictate phase behavior and finite-size corrections (Corral et al., 2019).
• Reinforcement of universality: Mixture-of-experts, multi-modal, and configuration-to-performance mapping in deep learning robustly fit into the joint-scaling paradigm, enabling predictive modeling of performance across multiple axes under large-scale heterogeneity (Zhang et al., 10 Feb 2026).
• Future directions: Open questions remain in understanding artifact regime transitions (largest outliers, foothill anomalies), interaction with quantization or radical architecture shifts, and the limiting behavior for highly non-polynomial or ultrahigh-dimensional data (Chen, 2011, Bi et al., 25 Sep 2025).

7. Representative Joint Scaling Laws Across Fields

Domain	Observables/Parameters	Joint Scaling Formulation	Key Constraint or Analytic Link
Linguistics	Word freq $L$0, length $L$1	$L$2	$L$3 and $L$4 determine both Heaps' and Zipf's law
Econophysics	Sales $L$5, labor $L$6	$L$7 and $L$8	Exponents $L$9 uniquely determine joint and marginals
Multilingual LMs	$P(n;L)$0, $P(n;L)$1, $P(n;L)$2	$P(n;L)$3 as explicit function in all variables	Loss per family depends only on $P(n;L)$4, exponents universal
City-size hierarchies	Rank $P(n;L)$5, class $P(n;L)$6	$P(n;L)$7 vs. $P(n;L)$8	Exponents $P(n;L)$9, $n$0 reciprocally tied via construction
Genomics	Family size $n$1, category $n$2	$n$3, $n$4	$n$5 links evolutionary and functional exponents

These joint scaling laws demonstrate the universal underlying regularities that emerge from complex system interactions and foster trans-disciplinary advances in theory, modeling, and predictive analytics.