---
title: 'Isomorph: Invariance in Physics, Logic & Supply Chains'
url: https://www.emergentmind.com/topics/isomorph
type: topic
---

# Isomorph: Invariance in Physics, Logic & Supply Chains

In contemporary research, **isomorph** denotes several distinct but technically related ideas. In condensed-matter physics, it names a curve in a thermodynamic phase diagram along which structure and dynamics are invariant when expressed in reduced units, a central concept of isomorph theory for Roskilde-simple systems [1406.2216]. In logic and theoretical computer science, it appears in formal treatments of isomorphism, including forcing-based notions of “probably isomorphic” structures, set-theoretic semantics for dependent type theory, and mechanized data-type transformations in ACL2 [2507.01518] [1407.7274] [2009.13771]. More recently, **ISOMORPH** has also been adopted as the name of a supply-chain digital twin for simulation, dataset generation, and forecasting benchmarks [2605.12768].

## 1. Isomorphs in condensed-matter theory

In isomorph theory, an isomorph is a curve in the phase diagram \((\rho,T)\) consisting of state points with constant excess entropy \(S_{\rm ex}\equiv S-S_{\rm ideal}\) [1307.5237]. Two state points \((\rho_1,T_1)\) and \((\rho_2,T_2)\) are defined as isomorphic if for physically relevant pairs of microconfigurations with coordinates that match when scaled by density,  
\[
\rho_1^{1/3} \mathbf{R}_1 = \rho_2^{1/3} \mathbf{R}_2,
\]
the Boltzmann factors are proportional,
\[
\exp(-U(\mathbf{R}_1)/k_B T_1)=C_{1,2}\exp(-U(\mathbf{R}_2)/k_B T_2).
\]
A Roskilde-simple system is defined by the property that the order of the potential energies of configurations at one density is maintained when these are scaled uniformly to a different density [1406.2216].

The underlying physical statement is hidden scale invariance. A convenient representation is
\[
U(\mathbf{R}) \approx h(\rho)\,\Phi(\rho^{1/3}\mathbf{R}) + g(\rho),
\]
where \(h(\rho)\) sets the density-dependent energy scale and \(g(\rho)\) is an additive term that does not affect forces, structure, or dynamics [1810.07255]. This implies that, along curves satisfying \(h(\rho)/T=\mathrm{const}\), the canonical configurational distributions are the same in reduced units. Standard reduced units are
\[
\tilde r = \rho^{1/3} r,\qquad
\tilde t = t\,\rho^{1/3}\sqrt{\frac{k_B T}{m}},\qquad
\tilde f = \frac{f}{\rho^{1/3}\sqrt{k_B T/m}},
\]
with corresponding reduced transport quantities such as \(D^*=D\rho^{1/3}/\sqrt{k_B T/m}\) and \(\eta^*=\eta/(\rho^{2/3}\sqrt{m k_B T})\) [1307.5237].

A practical characterization uses equilibrium virial–potential-energy correlations. The Pearson correlation coefficient is
\[
R=\frac{\langle \Delta W \Delta U\rangle}{\sqrt{\langle (\Delta W)^2\rangle \langle (\Delta U)^2\rangle}},
\]
and the density-scaling exponent is
\[
\gamma=\frac{\langle \Delta W \Delta U\rangle}{\langle (\Delta U)^2\rangle}
=\left(\frac{\partial \ln T}{\partial \ln \rho}\right)_{S_{\rm ex}}.
\]
Roskilde-simple systems typically satisfy \(R>0.9\), and along an isomorph one has \(d\ln T=\gamma\, d\ln \rho\) [1004.5142]. A recurrent misconception is that isomorphs require exact inverse-power-law interactions. The Buckingham-liquid study showed instead that strong correlations and isomorphs do not depend critically on the mathematical form of the repulsion being an inverse power law [1106.2973].

## 2. Invariants, response functions, and experimental tests

A central strength of isomorph theory is that it yields explicit invariants for measurable response functions. For dielectric spectroscopy of isotropic molecular liquids, one can combine linear response, the fluctuation-dissipation theorem, and the isomorph condition \(\Gamma\equiv \rho^\gamma/T=\mathrm{const}\) to obtain four equivalent isomorph-invariant quantities:  
\[
\langle (V P(\tilde t)-V P(0))^2\rangle,\qquad
\chi_e(\tilde t)\,T/\rho,\qquad
\chi_e(\tilde t)\,\rho^{\gamma-1},\qquad
F g,
\]
where \(\chi_e\) is the electric susceptibility and \(F g\) is the Kirkwood–Frölich factor product [1405.6098].

The explicit experimental prediction tested for the van der Waals liquid 5-phenyl-4-ether (5PPE) was
\[
-C_2''(f)=-C_1''(f)\left(\frac{\rho_1}{\rho_2}\right)^{\gamma-1},
\]
using the experimentally justified approximation \(C_{\rm empty}(2)\approx C_{\rm empty}(1)\) [1405.6098]. The liquid was studied at isochronal states in the temperature range \(266\!-\!333\) K and pressure range \(0.1\!-\!300\) MPa, with relaxation times around \(10^{-3}\) s and \(10^{-4}\) s. The dielectric setup achieved amplitude reproducibility \(\pm 0.1\%\), the empty capacitance was measured as \(51\) pF at ambient pressure, and 42 pairs of isochronal states were identified from 155 measured state points. Using \(\gamma=5.5\) from prior work and densities from a Tait equation of state, the predicted and measured loss spectra showed good match, with a relative deviation at the loss peak within \(-3\%\) to \(+19\%\) for the \(\chi_e \rho^{\gamma-1}\) prediction [1405.6098].

Dynamic mechanical analysis yields an analogous invariance statement for viscoelastic response. In simulations of the Kob–Andersen binary Lennard-Jones system, the standard reduced modulus is
\[
\tilde G = \frac{G}{\rho k_B T},
\]
with reduced frequency
\[
\tilde\omega=\frac{\omega}{\rho^{1/3}\sqrt{k_B T/m}}.
\]
Using SLLOD with time-dependent strain rates to impose sinusoidal shear, reduced loss-modulus curves \(\tilde G''\) were found to collapse along isomorphic state points when plotted against the unscaled temperatures, provided the isomorphic temperatures were chosen so that reduced forces remain invariant [2204.06962]. The study considered densities \(1.3\), \(1.35\), \(1.50\), and \(2.00\), a reference angular frequency \(10^{-4}\), and strain amplitude \(0.02\). A force-based construction using temperature-matched configurations remained applicable at \(\rho=2.0\), whereas a quenched-configuration variant broke down for the largest density rescalings because of force-vector decorrelation [2204.06962].

## 3. Scope across liquids, polymers, crystals, plasmas, and confinement

The theory was first developed for strongly correlating liquids and generalized Lennard-Jones systems, for which isomorphs were shown to be curves along which structure, dynamics, and some thermodynamic properties are invariant in reduced units [1004.5142]. For generalized Lennard-Jones pair potentials with exponents \(m\) and \(n\), the density-scaling function has the form
\[
h(\rho)=C_m \rho^{m/3}+C_n \rho^{n/3},
\]
and in the standard \(12\!-\!6\) case all isomorphs can be scaled onto a single curve in the \(W\!-\!U\) plane [1004.5142]. This already implied that melting and freezing lines are closely tied to isomorphs on the liquid side.

Subsequent work broadened the class of model systems. A binary Buckingham liquid, despite its exponential repulsion, was shown to be strongly correlating and to possess isomorphs; reduced radial distribution functions and both incoherent and coherent intermediate scattering functions were approximately invariant, and the viscous-state dynamics were closely mimicked by a purely repulsive inverse-power-law reference system with exponent \(n\approx 14.71\) [1106.2973]. Flexible Lennard-Jones chains likewise exhibited isomorphs: segmental and chain-center-of-mass incoherent intermediate scattering functions, end-to-end vector autocorrelation functions, most Rouse-mode correlators, and mean-square displacements collapsed across density changes up to about \(25\%\) in reduced units, while jumps between isomorphic state points produced instantaneous equilibration without slow relaxation [1307.5237].

The concept also extends to crystals. Simulations of face-centered-cubic Lennard-Jones crystals showed that reduced radial distribution functions, velocity autocorrelation functions, phonon dynamics, and even slow vacancy-jump dynamics are approximately invariant along isomorphs [1406.1911]. Other crystalline systems with isomorphs included the Wahnström binary Lennard-Jones crystal with the \({\rm MgZn_2}\) Laves structure, monatomic FCC Buckingham crystals, a purely repulsive finite-separation model, and an ortho-terphenyl molecular crystal. In contrast, a NaCl crystal model and SPC/E hexagonal ice did not exhibit isomorph invariances, supporting the broader conjecture that crystalline solids with isomorphs include most or all formed by atoms or molecules interacting via metallic or van der Waals forces, whereas covalently- or hydrogen-bonded crystals are not expected to have isomorphs [1406.1911].

High-pressure metallic crystals provided a further extension. Molecular-dynamics simulations using effective medium theory for Au, Ni, Cu, Pd, Ag, and Pt found strong hidden scale invariance at condensed-state densities: reduced-unit radial distribution functions collapsed almost perfectly along isomorphs, reduced velocity autocorrelation functions and vibrational density-of-states proxies showed good collapse, and jumps between isomorphic points led to instantaneous equilibration [1810.07255]. A notable difference from simple Lennard-Jones behavior is that \(\gamma\) varies substantially with density. For the Au crystal isomorph, as density increased from \(19.3\) to \(38.6\) g/cm\(^3\) and pressure from \(10\) to \(710\) GPa, \(\gamma\) decreased from \(6.45\) to \(1.92\), while \(R\) increased from \(0.985\) to \(0.998\) [1810.07255].

Soft-matter and plasma systems also fall within the scope of isomorph theory when virial–potential-energy correlations are strong. For the Yukawa fluid, \(R>0.99\) at all simulated state points, and isomorphs identified by both direct isomorph check and an analytical construction displayed invariance of the reduced radial distribution function, static structure factor, mean-square displacement, and incoherent intermediate scattering function [1505.06706]. The analytically derived form
\[
\Gamma(\kappa)=\Gamma_0 \,\frac{2 e^{\Lambda \kappa}}{(\Lambda \kappa)^2 + 2\Lambda \kappa + 2}
\]
reproduces the known melting-line shape, which the theory interprets as an isomorph [1505.06706].

Confinement does not automatically destroy isomorphs. In slit-pore simulations with crystalline walls, both the single-component Lennard-Jones liquid and the Kob–Andersen binary Lennard-Jones mixture retained good isomorph behavior for pore widths \(H\gtrsim 4\): reduced density profiles parallel and perpendicular to the walls, reduced mean-square displacements, and higher-order structures from topological cluster classification were nearly invariant along confined isomorphs [2102.06949]. The breakdown at \(H\approx 2\), where \(R<0.90\), is a useful counterexample to the assumption that reduced-unit invariance is automatic in strongly inhomogeneous environments [2102.06949].

## 4. Nonequilibrium, shear, and topological extensions

Isomorph theory is not restricted to equilibrium. For Couette shear flows generated by the SLLOD equations of motion, the reduced equations are identical along an isomorph provided the reduced strain rate
\[
\tilde{\dot\gamma}= t_0 \dot\gamma,\qquad
t_0=\rho^{-1/3}\sqrt{m/k_B T},
\]
is held fixed [1212.4480]. Under this condition, simulations of both the single-component Lennard-Jones liquid and the Kob–Andersen binary Lennard-Jones mixture showed collapse of the reduced radial distribution function, the transverse self-intermediate scattering function, and the reduced viscosity
\[
\tilde \eta=\frac{\eta}{\rho^{2/3}T^{1/2}}
\]
as a function of reduced strain rate, in both linear and shear-thinning regimes [1212.4480].

A more general nonequilibrium reformulation introduces the **systemic temperature** \(T_{\rm s}\), defined as the temperature of the equilibrium state point with average potential energy equal to \(U(\mathbf{R})\) [2008.02590]. Systemic isomorphs are lines of constant excess entropy in the phase diagram defined by density and systemic temperature, and the reduced dynamics is invariant along a systemic isomorph if there is a constant ratio between the systemic and the bath temperature. In thermal equilibrium, \(T_{\rm s}=T_{\rm b}\) and the original formalism is recovered [2008.02590]. This framework rationalizes earlier observations of isomorph invariance in nonlinear steady-state shear flows, zero-temperature plastic flows, and glass-state isomorphs.

A different extension uses topological information rather than thermodynamic excess entropy as the scaling variable. For soft-sphere fluids with repulsive \(n\!-\!6\) Mie potentials, the Shannon entropy of the Voronoi-cell-topology distribution,
\[
S_V = -\sum_i p_i \ln p_i,
\]
was shown to provide a scaling law for reduced transport properties comparable to conventional excess-entropy scaling [1901.02772]. Across \(n=8,12,16,20,24\), the Voronoi excess entropy and thermodynamic excess entropy were almost linearly related, with \(R^2\) between \(0.994\) and \(0.996\). The work further suggested that the Frenkel line is a topological isomorphic line, marked by \(\hat{\mathcal H}\approx 12\) and \(\Pi_{\rm solid}\approx 0.116\), where the functional form of the scaling relation changes qualitatively [1901.02772].

At the same time, higher-order structure can expose the limits of approximate isomorph invariance. In the Kob–Andersen Lennard-Jones glassformer and its mapped inverse-power-law reference system, two-point structure and dynamics were nearly identical, but topological cluster classification showed that bicapped square antiprisms, the locally favored 11A structures, had populations up to \(80\%\) higher in the Lennard-Jones system and lifetimes up to \(40\%\) higher than in the inverse-power-law reference system [1307.5516]. The structural relaxation times were almost identical, while the four-point dynamical susceptibility was marginally higher in the inverse-power-law system. This indicates that higher-order structural observables need not be as tightly constrained by isomorph theory as two-point reduced-unit observables [1307.5516].

## 5. Isomorphism in logic, type theory, and program transformation

Outside condensed matter, the term appears in formal semantics and theorem proving. In a set-theoretic formulation of dependent type theory, types are divided into small and large types—sets and proper classes respectively—and each proper class, such as “group” or “topological space,” has an associated notion of isomorphism [1407.7274]. Isomorphism is handled by defining a groupoid structure on the space of all definable values. The values are simultaneously objects and morphisms—“morphoids”—which supports sound inference rules for deriving isomorphisms and for substitution of isomorphics [1407.7274].

A forcing-theoretic generalization was introduced under the name **probably isomorphic**. Two structures \(M,N\) in the same language are probably isomorphic if they, or in the metric case their completions, are isomorphic after forcing with the Lebesgue measure algebra [2507.01518]. For discrete structures, or extremal models of a non-degenerate simplicial theory, the paper proved the equivalence
\[
M \text{ and } N \text{ are probably isomorphic}
\quad\Longleftrightarrow\quad
L^1([0,1],M)\cong L^1([0,1],N),
\]
thereby linking forcing-based isomorphism to randomization structures in continuous logic [2507.01518].

In mechanized program derivation, isomorphism is treated operationally. In ACL2, “types” are represented by predicates `old` and `new`, and an isomorphism is a pair of total ACL2 functions `iso` and `osi` that are inverse bijections when restricted to those domains [2009.13771]. The APT tools implement this through `defiso`, `isodata`, and `propagate-iso`. Once versions of the interface functions of a data type have been derived on the isomorphic representation, higher-level functions can be generated by substitution, and the tools automatically produce proofs of equivalence [2009.13771]. The paper gives examples ranging from refinement of finite sets to duplicate-free ordered lists or bit vectors to record extensions that cache derived fields.

## 6. ISOMORPH as a supply-chain digital twin

In a distinct and acronymic usage, **ISOMORPH** names a public, open-source digital twin of a multi-echelon supply-chain logistics network [2605.12768]. The simulator advances a directed routing graph in discrete time: demand arrives at the destination, is served from stock or recorded as backlog, and triggers replenishment through the network. The state vector
\[
\xi_t=(OH_t,B_t,Out_t,IT_t,\tilde\lambda_t)
\]
tracks on-hand inventory, destination backlog, outstanding orders, in-transit shipments, and a smoothed demand estimate, closing the dynamics as a Markov chain on a hybrid state space [2605.12768].

The released benchmark includes catalogue sizes \(C=50\) and \(C=200\), a horizon \(T=52{,}560\), six one-at-a-time scenario sweeps producing 30 additional rollouts, and 20 Latin-hypercube perturbations over demand-side parameters [2605.12768]. The demand process is
\[
y_{i,t}\sim \mathrm{Poisson}(\lambda_{i,t}),
\]
with \(\lambda_{i,t}\) composed of yearly and weekly seasonality, clipped AR(1) drift, per-item bursts, and shared macro-shocks. Replenishment follows \((s,S)\) rules, routing uses Dijkstra weights \(\tau_e/(K_eV_e)\), and dispatch to the destination is governed by
\[
\mathrm{target}_{i,t}=\max\{0, B^{d^\star,i}_t + m \tilde\lambda^i_t - IT^i_t - OH^{d^\star,i}_t\},
\]
with pipeline multiplier \(m=7\) [2605.12768].

The system was designed to encode three pathwise conservation laws, including per-node mass conservation
\[
OH^{n,i}_{t+1}=OH^{n,i}_t + R^{n,i}_t - D^{n,i}_t
\]
and global internal-network mass conservation
\[
I_i(t+1)=I_i(t)+A^{\rm src}_{i,t}-S_{i,t},
\]
which serve as verification tools for simulator extensions [2605.12768]. The released data reproduces the bullwhip effect at empirically consistent magnitudes: tier-level monthly means on \(C=50\) fall in \([1.07,1.57]\), within the \([0.97,1.90]\) interval spanning the \(50\)th–\(90\)th percentiles reported in the literature [2605.12768].

The benchmark was also used for zero-shot forecasting with Chronos-T5, Moirai-1.1-R, TimesFM-2.0, and Lag-Llama. Using context length \(L=512\), horizons \(h\in\{1,7,14,30\}\), and MASE relative to a seasonal-naive baseline, TimesFM achieved \(0.742, 0.786, 0.843,\) and \(0.978\) on the \(C=50\) baseline, while Lag-Llama yielded \(1.027, 1.292, 1.396,\) and \(1.601\) [2605.12768]. The same digital twin and Latin-hypercube perturbations were used to generate forward uncertainty-quantification bands, illustrating a role for foundation models as fast surrogates for forward UQ under parameter uncertainty [2605.12768].

## 7. Conceptual unity and points of divergence

Across these domains, the common theme is not a single formalism but a recurring principle: a complicated object or process is organized by an invariance under a structured transformation. In condensed matter, the relevant transformation is uniform scaling in density and temperature, and the main content of the theory is reduced-unit invariance along configurational adiabats [1406.2216]. In logic and computer science, the transformation is relabeling, transport, or representation change under a bijection, with correctness expressed by commuting diagrams or equivalence theorems [1407.7274] [2009.13771]. In the supply-chain digital twin, the name functions as an acronym rather than a direct statement of mathematical invariance, although the simulator itself is organized around an explicit state-space and conservation structure [2605.12768].

Several boundary conditions recur. In physics, isomorph invariance is strongest for Roskilde-simple liquids, many van der Waals systems, and many metallic crystals, but may fail for hydrogen-bonded, ionic, network-forming, or ultra-confined systems [1405.6098] [1406.1911] [2102.06949]. Power-law density scaling is only an approximation to the more general \(h(\rho)/T\) scaling and breaks down across sufficiently large density ranges [1307.5237]. Higher-order structural observables may vary even when two-point reduced-unit observables collapse well [1307.5516]. In formal settings, by contrast, isomorphism is exact once the relevant domains and transport maps are specified, but the scope of admissible constructions is determined by the logic or theorem-proving environment [2507.01518] [2009.13771].

The term **isomorph** therefore names a family of ideas rather than a single object. Its most developed physical usage is the one codified by isomorph theory: curves of constant excess entropy along which reduced structure and dynamics are invariant. Its broader technical use retains the same structural intuition—equivalence under a transformation that preserves the relations of interest—even when the transformation is a forcing extension, a data-type refinement, or a configurable logistics simulator.

Source: https://www.emergentmind.com/topics/isomorph