---
title: 'Density Scaling Laws: Fundamentals & Applications'
url: https://www.emergentmind.com/topics/density-scaling-laws
type: topic
---

# Density Scaling Laws: Fundamentals & Applications

Density scaling laws are relations in which observables collapse when expressed as functions of density, inverse density, or density-weighted fields rather than as independent functions of multiple control parameters. In the cited literature, “density” denotes specific volume inverse in glass-forming liquids, packing fraction in dense particulate systems, mass density in fluids and self-gravitating matter, population density in settlement systems, and dataset density in machine learning corpora. The common mathematical theme is reduction to a master curve or a piecewise power law, but the physical content ranges from hidden scale invariance and isomorph theory to jamming constraints, inverse cascades, hydrodynamic self-similarity, and breakpoint structure in heterogeneous spatial systems [1012.1612] [1705.06856] [2509.14258].

## 1. Canonical forms of density scaling

A basic form is the single power law
\[
D_i(\rho)=a\,\rho^b,
\]
where \(D_i(\rho)\) is an indicator density, \(a\) is a pre-exponential constant, and \(b\) is the scaling exponent. In the rural–urban setting, \(b<1\) is sublinear, \(b=1\) linear, and \(b>1\) superlinear. The same literature also uses a segmented law with a breakpoint \(\rho_c\),
\[
D_i(\rho)=
\begin{cases}
a_1\,\rho^{b_1}, & \rho<\rho_c,\\
a_2\,\rho^{b_2}, & \rho\ge \rho_c,
\end{cases}
\]
which distinguishes rural and urban regimes [2509.14258].

In viscous liquids and polymers, the canonical variable is not \(\rho\) itself but the combined density–temperature scaling variable
\[
X=T\,V^\gamma,
\]
with \(V\) the specific volume. Relaxation time, viscosity, and diffusion are then written as
\[
\tau(T,P)=F(TV^\gamma),\qquad \eta(T,P)=G(TV^\gamma),\qquad D(T,P)=H(TV^\gamma).
\]
For strongly correlating liquids, a more general form replaces the power law by an isomorph function,
\[
\tilde\tau=f\!\bigl(g(\rho/\rho_*)/T\bigr),
\]
so that density enters through a nontrivial function \(g(\rho)\) rather than a constant exponent \(\gamma\) [1012.1612] [1112.1602].

In disordered granular packing, density scaling is local and geometric. The average distance of the first \(n\) neighbors is
\[
R_n=\frac{1}{n}\sum_{i=1}^n r_i,
\]
and the packing fraction obeys
\[
\phi \propto R_n^{-\alpha_n},\qquad
D_R(n)=-\frac{\partial \ln \phi}{\partial \ln R_n}.
\]
Here the exponent is shell-dependent and reaches a minimum at the first-shell boundary [1705.06856].

In compressible MHD turbulence and compressible turbulence more broadly, density appears as a weighting factor. For solar-wind Elsässer fields,
\[
W^\pm(x)=\rho(x)^{1/3} Z^\pm(x),
\]
and the third-order law is written in terms of \(W^\pm\) rather than \(Z^\pm\). In compressible turbulence spectra, the density-weighted velocity \(v=\rho^{1/3}u\) plays an analogous role [1003.0533] [1502.02815].

Stationary driven fluids admit yet another master-curve representation. For uniaxial compressible Navier–Stokes–Fourier flows, the density field can be written as
\[
\rho(x)=p^\alpha F\!\bigl(p^\beta x,\omega(x)/p^\gamma\bigr),
\]
and in the separable case the exponents are \(\alpha=0\), \(\beta=-1\), \(\gamma=1\), with \(\omega(x)=\tfrac12[v(x)+J/\sigma]^2\) the kinetic field [2512.13205].

## 2. Liquids, supercooled states, and glass-transition scaling

Density scaling in glass-forming liquids and polymers is the empirical observation that structural relaxation times, viscosities, and diffusion constants measured over a range of temperatures and pressures superpose when plotted against \(T V^\gamma\). This has been demonstrated for roughly a hundred materials. The theoretical rationale is the inverse power-law approximation
\[
U(r)=\varepsilon (\sigma/r)^n,
\]
for which the scaling exponent is
\[
\gamma=\frac{n}{3}.
\]
However, the strict derivation applies to reduced quantities, not raw experimental observables. The reduced dynamical variables are
\[
\tau^*=V^{-1/3}\Bigl(\frac{k_B T}{m}\Bigr)^{1/2}\tau,\qquad
\eta^*=V^{2/3}\Bigl(\frac{m}{k_B T}\Bigr)^{1/2}\eta,\qquad
D^*=V^{-1/3}\Bigl(\frac{k_B T}{m}\Bigr)^{1/2}D.
\]
Only the exponent obtained from reduced quantities, \(\gamma^*\), can be sensibly related to the intermolecular potential. In the deeply supercooled regime the difference between reduced and unreduced scaling is small, but above the melting point it can be substantial; in dodecane, for example, unreduced viscosity gives \(\gamma\approx 6.5\), whereas reduced viscosity gives \(\gamma^*\approx 5.2\) [1012.1612].

Power-law density scaling is not exact over arbitrarily wide density ranges. In the Kob–Andersen binary Lennard–Jones mixture, a single exponent collapses any two isochores well but fails to collapse the three isochores \(\rho=1.2\), \(1.6\), and \(2.0\) simultaneously. Isomorph theory replaces
\[
\tilde\tau=F(\rho^\gamma/T)
\]
by
\[
\tilde\tau=f\!\bigl(g(\rho/\rho_*)/T\bigr).
\]
For a standard 12–6 Lennard–Jones potential,
\[
g(\tilde\rho)=c\,\tilde\rho^4+(1-c)\tilde\rho^2,
\qquad \tilde\rho=\rho/\rho_*.
\]
This generalization preserves single-variable collapse while allowing the local effective exponent \(\gamma(\rho)=d\ln g/d\ln \rho\) to vary with density [1112.1602].

Volumetric data can be incorporated into the same framework. In Lennard–Jones-based model systems, the same exponent that scales dynamics also scales \(P\)–\(v\) data through an equation of state derived from the short-range effective potential \(U_{\rm eff}(r)\simeq A r^{-m_{\rm IPL}}+\) weak background, with \(m_{\rm IPL}/3\equiv \gamma_{\rm EOS}\). In the KABLJ model, \(\gamma_{\rm dyn}\simeq 4.84\pm0.13\) from raw \(\tau\), \(\gamma_{\rm dyn,reduced}\simeq 4.98\pm0.11\), and \(\gamma_{\rm EOS}=4.83\pm0.16\); in the Lewis–Wahnström OTP model, \(\gamma_{\rm dyn}\simeq 8.0\) and \(\gamma_{\rm EOS}=8.20\pm0.23\) [1112.4001].

Generalized equations of state make the density dependence explicit through a one-parameter function \(h(\rho)\). For KABLJ,
\[
h(\rho)=c\rho^4+(1-c)\rho^2,\qquad
\gamma(\rho)=\frac{d\ln h}{d\ln \rho}
=\frac{4c\rho^4+2(1-c)\rho^2}{c\rho^4+(1-c)\rho^2}.
\]
Fitting volumetric data gives \(c=1.74\pm0.03\) from the configurational EOS and \(c=1.70\pm0.03\) from the total-pressure EOS, and the same \(h(\rho)\) collapses all eleven isotherms of reduced structural relaxation times when plotted against \(h(\rho)/T\) [1311.3910].

Packing fraction enters analytically in generalized mode-coupling theory through the reduced control parameter
\[
\epsilon=\frac{\phi-\phi^c}{\phi^c}.
\]
Near the ideal glass transition, GMCT yields two divergent time scales,
\[
\tau_\beta=t_0(\sigma|\epsilon|)^{-1/(2a)},\qquad
\tau_\alpha=t_0 B^{-1/b}(\sigma|\epsilon|)^{-(1/(2a)+1/(2b))},
\]
with
\[
\gamma=\frac{1}{2a}+\frac{1}{2b}.
\]
For the Percus–Yevick hard-sphere system under MF closures, the critical packing fractions are \(\phi^c\approx0.5159\), \(0.5319\), and \(0.5442\) for \(N=2\), \(3\), and \(4\), respectively, while the exponents shift systematically as higher-order correlations are included [2008.10550].

## 3. Granular packing, jamming, and non-cubic density–distance scaling

In disordered granular packing of spherical particles, X-ray tomography identifies a non-cubic scaling law between packing fraction and characteristic nearest-neighbor distance,
\[
\phi\propto r^{-\alpha},\qquad \alpha\approx 2.5.
\]
More precisely, for the average distance \(R_n\) of the first \(n\) neighbors,
\[
\phi\propto R_n^{-\alpha_n},
\]
and the minimum exponent occurs at \(n=13\), corresponding to the first-shell boundary. Over packings with global packing fraction \(\phi\in[0.57,0.64]\),
\[
\phi\propto R_{13}^{-D_R(13)},\qquad D_R(13)=2.54\pm0.03.
\]
The first-shell boundary is located at \(\langle r_{13}\rangle\approx 1.37\,d\) [1705.06856].

The experimental system consists of \(200\,\mu{\rm m}\) glass beads with \(\sim 3\%\) polydispersity, packed in a cylindrical cell by tapping, hopper deposition, and flow-pulse protocols. Synchrotron X-ray CT at \(5.5\,\mu{\rm m}\) voxel resolution yields particle centroids and radii with \(\sim0.3\%\) diameter precision, and the analysis is restricted to beads at least \(3\) diameters from walls, amounting to \(\sim1.7\times10^4\) beads per scan. Packing fraction is defined as
\[
\phi = \frac{N\,V_{\rm bead}}{V_{\rm sample}},\qquad V_{\rm bead}=\frac{\pi}{6}d^3.
\]
These definitions matter because the observed exponent is extracted from a strictly local neighbor statistic rather than from an assumed continuum field [1705.06856].

The non-cubic exponent is traced to the internal structure of the first neighbor shell. Within \(n=1\ldots13\), the six closest neighbors change distance very little with \(\phi\), corresponding to local exponents \(D_i>3\), whereas neighbors \(n=7\ldots13\) shrink strongly with \(\phi\), corresponding to \(D_i<3\). The measured \(D_R(13)\approx2.5\) is the weighted average of these two groups. Conditional on the quasi-contact number \(z\), defined as neighbors whose surface separation is \(<0.01d\), the local law is
\[
\phi\big|_z \propto R_{13}^{-d_z},
\]
with \(d_z\) increasing from \(\sim2.5\) at \(z\approx6\) toward \(3\) for smaller \(z\). The decomposition
\[
d=\langle d_z\rangle_z+\text{(small cross-terms)}
\]
shows that more than \(80\%\) of the deviation from \(3\) arises from the \(z\)-dependence alone [1705.06856].

The physical interpretation is therefore jamming-related. The exponent \(\sim2.5\) emerges close to the isostatic, marginally jammed state \(z\approx6\), while as \(\phi\) decreases and \(z\) drops below \(6\), the local exponent rises toward \(3\) and the global exponent rises from \(\sim2.6\) to \(\sim2.9\). Metallic glasses report \(\alpha\approx2.3\)–\(2.5\) under density changes, and these values had been attributed to “fractal” medium-range order. The granular results show that no true fractal is required: the non-cubic law is fully explained by local non-affine displacements of contact versus non-contact neighbors and their isostatic count. Friction, bead stiffness, thermal effects, and departures from the hard-sphere jamming point shift \(z\) and therefore modify \(\alpha\) [1705.06856].

## 4. Density weighting in turbulence and driven fluids

In fast polar solar wind, the incompressible MHD analogue of the Kolmogorov–Yaglom law is
\[
Y^\pm(\ell)=\left\langle |\Delta Z^\pm(\ell)|^2\,\Delta Z_\parallel^\mp(\ell)\right\rangle
= -\frac{4}{3}\epsilon^\pm |\ell|,
\]
where \(Z^\pm=v\pm b/\sqrt{4\pi\rho}\). To incorporate compressible fluctuations, a phenomenological density-weighted form introduces
\[
W^\pm(x)=\rho(x)^{1/3}Z^\pm(x),
\]
and the corresponding third-order law becomes
\[
W^\pm(\ell)=\frac{\left\langle |\Delta W^\pm(\ell)|^2\,\Delta W_\parallel^\mp(\ell)\right\rangle}{\langle \rho\rangle}
=-\frac{4}{3}\epsilon^\pm |\ell|.
\]
The Ulysses analysis uses 8-minute averages in the first half of 1996, 11-day sliding windows, and the Taylor hypothesis \(\ell\approx-\langle v_R\rangle \tau\). The compressible law is observed in \(\simeq50\%\) of the time for \(W^+\) and \(\simeq25\%\) for \(W^-\), whereas the incompressible law appears in only \(\simeq33\%\) of the time. Despite density fluctuations of only \(\delta\rho/\langle\rho\rangle\lesssim10\%\), the compressible cascade rate is an order of magnitude larger than the incompressible one and is comparable to the heating rate required to explain the non-adiabatic solar wind; the compressible derivation, however, remains purely phenomenological and lacks an exact proof in compressible MHD [1003.0533].

Dimensional analysis gives analogous density-weighted scaling in compressible turbulence. The ordinary kinetic-energy spectrum takes the form
\[
E(k,\epsilon,\nu,\rho,M)\simeq \epsilon^{2/3}k^{-5/3} f[(k\eta)^{4/3},M],
\]
while the density-weighted spectrum satisfies
\[
E_\rho(k,\epsilon,\nu,\rho,M)\simeq \rho^{1/3}\epsilon^{2/3}k^{-5/3} F[(k\eta)^{4/3},M],
\]
with \(\eta=(\nu^3/\epsilon)^{1/4}\). Using Barenblatt incomplete similarity,
\[
E_\rho(k,\epsilon,M)=C_2(M)\,\epsilon^{2/3+\beta_\rho/\ln M}\,
k^{-5/3+(4\beta_\rho)/(3\ln M)}.
\]
As \(M\to0\), \(M\to1\), and \(M\to\infty\), the effective slope tends to \(-5/3\). The density-weighted velocity \(v=\rho^{1/3}u\) is introduced because it more robustly follows the \(k^{-5/3}\) law across a wide range of Mach numbers [1502.02815].

Stationary compressible Navier–Stokes–Fourier flows admit a different form of density scaling. In uniaxial steady states, pressure \(P\), shear stress \(\sigma\), and heat current \(J\) are constants, and elimination of temperature yields
\[
{\cal G}_P(\rho(x))=\omega(x)+\xi,\qquad
\omega(x)=\frac12\Bigl(v(x)+\frac{J}{\sigma}\Bigr)^2.
\]
Inverting \({\cal G}_P\) gives
\[
\rho(x)={\cal R}_P(\omega(x)+\xi),
\]
while spatial scaling of the kinetic field yields
\[
\omega(x)={\cal W}_{P,\xi}(\pm \sigma x+\zeta).
\]
For density–temperature separable constitutive laws, this simplifies to
\[
\rho(x)=p^0\,F\!\Bigl(p^{-1}x,\omega(x)/p\Bigr),
\]
with \(\alpha=0\), \(\beta=-1\), \(\gamma=1\). Large-scale molecular dynamics simulations of 2D hard disks and 3D Lennard–Jones fluids show excellent data collapse in the bulk after discarding boundary layers of order one cell wide [2512.13205].

Hydrodynamic self-similarity also governs ATI-formed plasma channels for laser wakefield accelerators. After the ionization pulse, the over-pressured plasma drives a cylindrical shock, and the late-time density profile depends only on
\[
r^*=\frac{r}{r_{\rm ion}^{1/2}},\qquad
\rho^*=\frac{\rho}{n_g}.
\]
For parabolic channels matched to a Gaussian driver, the on-axis density scales linearly with initial gas density,
\[
n_0\propto n_g,
\]
with fitted relations \(n_0=1.00\,n_g\) for \(N_2\) and \(n_0=0.86\,n_g\) for He. The matched spot radius obeys
\[
r_m\propto r_{\rm ion}^{1/2} n_g^{-1/4},
\]
with fit constants \(C=4.91\) for nitrogen and \(C=3.10\) for helium. Simulation profiles at \(t=0.5\), \(1.5\), and \(4.0\,{\rm ns}\) collapse onto a single curve over \(n_g=10^{17}\)–\(10^{19}\,{\rm cm}^{-3}\) and \(r_{\rm ion}=20\)–\(50\,\mu{\rm m}\), with deviations from ideal exponents \(\lesssim10\%\) [2508.11238].

## 5. Gravitational, network, and settlement-scale density laws

A cascade-based theory of dark matter proposes an inverse kinetic-energy cascade with constant rate
\[
\varepsilon_u \simeq -4.6\times10^{-7}\,{\rm m}^2\,{\rm s}^{-3}.
\]
Dimensional analysis then gives the two-thirds law
\[
v_r^2\propto (\varepsilon_u r)^{2/3},
\]
which, combined with the virial theorem, yields the five-thirds law for enclosed mass and the four-thirds law for mean density,
\[
m_r\sim \varepsilon_u^{2/3}G^{-1}r^{5/3},\qquad
\rho_r\sim \varepsilon_u^{2/3}G^{-1}r^{-4/3}.
\]
At the scale radius \(r_s\),
\[
\rho_s\propto \varepsilon_u^{2/3}G^{-1}r_s^{-4/3}.
\]
For fully virialized haloes with vanishing radial flow, the asymptotic inner slope is
\[
\gamma=\frac{d\ln \rho}{d\ln r}=-\frac{4}{3}.
\]
Rotation-curve fits from SPARC, DMS, and SOFUE yield a best-fit slope \(\simeq-1.36\pm0.05\), but the same framework states explicitly that nonzero radial flow or non-steady accretion can steepen or flatten \(\gamma\) relative to \(-4/3\) [2209.03313].

Finite-density scaling laws also govern condensation in zero-range processes on scale-free networks. With \(N=\rho L\) particles and hopping rate \(u(n_i)=n_i^\delta\), condensation occurs for
\[
\delta<\delta_c\equiv \frac{1}{\gamma-1}.
\]
In the condensed regime,
\[
k_c\sim \frac{L^{\delta_c-\delta}}{\rho^\delta},
\]
and the average occupation takes the unified form
\[
m_k\sim \rho^\delta\,{\cal G}(k\rho^\delta),
\]
where \({\cal G}(y)\sim y\) for \(y\ll y_c\) and \({\cal G}(y)\sim y^{1/\delta}\) for \(y\gg y_c\). Relaxation is hierarchical, with
\[
\frac{I_t}{I_\infty}\sim \frac{t}{\rho^{1-\delta}},\qquad t\ll t_R,
\]
for the inverse participation ratio. Monte Carlo simulations on Barabási–Albert networks with \(\gamma=3\), \(\delta_c=1/2\), and \(\delta=0.2\) validate both the steady-state collapse and the transient scaling [1704.06750].

Population density introduces a segmented scaling paradigm in rural–urban systems. For Middle Layer Super Output Areas in England and Wales, \(n=7{,}080\), 117 indicators are converted to per-hectare densities and fitted in log-space by OLS. Segmented models are compared to single power laws using Davies’ test, AIC, and BIC, and 92 of the 117 indicators exhibit a significant breakpoint at
\[
\rho_c \approx 33\pm 5 \text{ persons per hectare}.
\]
Typical exponents have median rural \(b_1\approx0.9\) and urban \(b_2\approx1.2\) for accelerating phenomena, while some mortality indicators show urban exponents \(b_2\approx0.8\)–\(0.9\). Crime, property transactions, road accidents, and mortality display distinct rural-to-urban transitions, and finer MSOA resolution reveals segmented behaviors not visible in coarser units. For dementia and ischaemic heart disease stratified by older age groups, \(b_1\) remains near unity whereas \(b_2<1\), which is interpreted as an urban protective effect [2509.14258].

## 6. Universality, invariance, and breakdown

In biological continuum mechanics, density scaling appears as an invariance statement rather than as a fitted exponent. For a coupled system of incompressible Navier–Stokes fluid dynamics, nonlinear elasticity, and rigid-body mechanics, the scaling group is
\[
L\mapsto \lambda L,\qquad
T\mapsto \alpha T,\qquad
\rho\mapsto \eta \rho,
\]
with
\[
\nu\mapsto \frac{\lambda^2}{\alpha}\nu,\qquad
\mathbf g\mapsto \frac{\lambda}{\alpha^2}\mathbf g,\qquad
E,\mu,\kappa\mapsto \frac{\eta\lambda^2}{\alpha^2}(E,\mu,\kappa).
\]
Demanding invariance of the full PDE system forces all densities—fluid, rigid, stiff tissue, and soft tissue—to scale with the same factor \(\eta\). In a fixed terrestrial gravitational field, \(\lambda/\alpha^2=1\), so \(\alpha=\sqrt{\lambda}\), and the natural choice is \(\eta=1\). The resulting terrestrial density law is therefore size-invariant: \(\rho\propto L^0 T^0\) [2502.11398].

In large language models, “data density” denotes redundancy in embedding space rather than a physical density. Cluster density is defined by
\[
\rho_i=\frac{N_i\,\Gamma(n/2+1)}{\pi^{n/2} r_i^n},
\qquad
r_i=\frac{1}{N_i}\sum_{x\in C_i}\|x-c_i\|,
\]
and the overall dataset density is
\[
\rho=\frac{N\,\Gamma(n/2+1)}{\pi^{n/2}R^n},
\]
with a weighted centroid distance \(R\). High \(\rho\) indicates tightly packed, less diverse data. Classical scaling laws,
\[
L(C)=\lambda_C C^{-\alpha_C},\qquad
L(N,D)=E+\frac{\lambda_N}{N^{\alpha_N}}+\frac{\lambda_D}{D^{\alpha_D}},
\]
are extended by a sub-optimal law,
\[
L(N,D)=E+\frac{\lambda_N R_N}{N^{\alpha_N}}+\frac{\lambda_D R_D}{D^{\alpha_D}},
\]
with \(R_D\) and \(R_N\) logistic in the over-training ratio \({\rm OTR}=D/N\). Over more than 400 runs spanning 20 M to 7.03 B parameters, high-density regimes \((\rho\gtrsim0.56)\) and high OTR \((>50)\) show sub-scaling, the loss–compute exponent \(\alpha_C\) falls from \(\sim0.07\) at OTR \(=5\) to \(\sim0.052\) by OTR \(=50\), and the sub-optimal law reduces MAPE by 40–90% in dense-data regimes [2507.10613].

Across these literatures, observed exponents are frequently regime-dependent rather than universal. The granular \(\alpha\approx2.5\) law is tied to the isostatic contact number and varies with coordination, friction, and related parameters rather than indicating a true fractal structure [1705.06856]. Power-law density scaling in supercooled liquids is accurate over modest density ranges but fails over larger density variations, where \(g(\rho)\) replaces \(\rho^\gamma\) [1112.1602]. The compressible Yaglom law in the solar wind is phenomenological rather than exact [1003.0533]. The dark-matter inner slope \(-4/3\) requires vanishing radial flow and is modified by accretion history [2209.03313]. Rural–urban exponents depend on spatial granularity and demographic stratification, with a consistent breakpoint emerging only after fine-grained analysis [2509.14258]. This suggests that density scaling laws are best understood as symmetry-based, cascade-based, or empirically stabilized effective descriptions whose validity is set by coordinate choice, constitutive assumptions, and the dynamical regime under study.

Source: https://www.emergentmind.com/topics/density-scaling-laws