The Zeta-Minimizer Theorem: A Deductive Axiomatic Framework for Variational Thermodynamics, Chemical Equilibrium, and Heterogeneous Catalysis with Applications to Ammonia Synthesis

1. Introduction

Ammonia synthesis remains one of the most consequential chemical processes on the planet, underpinning global food production and emerging as a carbon-free energy vector. Despite more than a century of industrial practice, the Haber–Bosch process continues to operate under severe thermodynamic and kinetic constraints, requiring high pressures, elevated temperatures, and expensive promoted iron or ruthenium catalysts whose design still relies heavily on empirical scaling and trial-and-error optimization. Recent pushes toward low-pressure, modular, and renewable-powered ammonia plants have intensified the need for a truly predictive, first-principles framework capable of guiding both catalyst formulation and reactor design [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15].

Existing theoretical tools—density-functional scaling relations, microkinetic modeling, and Sabatier-type volcano plots—have provided invaluable qualitative guidance yet remain largely inductive or semi-empirical. Linear scaling relations, for example, are typically extracted from large DFT datasets rather than derived from thermodynamic axioms; their transferability across promoters, supports, and operating regimes is limited, and they offer little direct insight into the role of inerts, interface stability, or the coupling between gas-phase kinetics and vapor–solid equilibration. A deductive framework that begins from a small set of thermodynamic axioms and rigorously derives classical thermodynamic relations, equilibrium constants, reaction rates, and the emergence of scaling and volcano behavior would constitute a significant advance.

This work develops a fully deductive, parameter-free framework for heterogeneous catalysis and chemical reaction engineering [16,17], built from three primitive axioms and the helical geometry of the grand-partition function. The manuscript is organized as follows.

Section 2, Section 3 and Section 4 establish the axiomatic and geometric foundations of the Zeta-Minimizer Theorem (ZMT). The three axioms are introduced, the helical operator, and the effective Hessian PDE, followed by the categorical invariance of the grand-partition function and the solid–fluid interface matching condition that yields the interface constants $C_k$.

Section 5, Section 6 and Section 7 derive the variational emergence of classical thermodynamics and chemical equilibrium [18,19,20,21]. Showing how conjugate pairs, phase equilibria, extent of reaction, and multi-reaction multi-phase equilibrium arise directly from minimization of the phase functional, without empirical assumptions.

Section 8 presents the variational reaction rate law, its linearization near equilibrium, and the emergence of batch and residence time [16,17]. From the same framework deductively recovering linear scaling relations, Sabatier volcano plots, modal decomposition of coupled dynamics, and the classical seven-step heterogeneous catalysis sequence [16,17] as emergent normal modes.

Section 9 applies the complete variational model to industrial ammonia synthesis. The explicit phase functional for the N₂–H₂–NH₃–inert system are constructed, deriving the exact nonlinear driving force and scalar Hessian, obtaining the implicit analytical solution for transient conversion, and quantifying the relative impact of inerts and pressure.

Section 10 demonstrates that the classical Temkin–Pyzhev rate expression [16,17] emerges as a limiting case of the ZMT variational rate law in the high-inerts regime, while the ZMT formulation remains regular and physically consistent at low inerts and low pressure.

Section 11 provides experimental and DFT validation. Interpreting formation energies of adsorbed promoters, N adsorption energies, N–N transition-state barriers, and promoter-induced electric fields through the lens of interface constants $C_k$ and surface gear occupations, and we validate the framework against catalyst activity data obtained at constant pressure and zero inerts.

Section 12 offers a comparative industrial analysis. Contrasting the operating regimes of Fe- and Ru-based catalysts, re-examining recent low-pressure ammonia synthesis studies through the ZMT lens, and discuss the critical role of inerts and the lab–industry gap.

2. ZMT Axiomatic Foundation and Helical Geometry

The Zeta-Minimizer Theorem (ZMT) provides a deductive unification of statistical mechanics, number theory, and helical geometry [22,23,24]. All macroscopic thermodynamics, electromagnetic phenomena, and magnetic transport emerge from three primitive axioms. These axioms uniquely select the non-proper Archimedean conical helix as the background geometry satisfying bounded oscillations, flux conservation, and entropy maximization (Topology Selection Theorem S4 [22]) [23,24].

2.1. The Three Axioms

2.1.1. Axiom 1 (Strict Concave Entropy Maximization)

The phase functional $F$ is strictly concave and possesses a unique global maximum (Lemma S1 [22]) [24]. In the grand canonical ensemble this implies that the entropy [22,23,24,25,26,27,28]

$$\begin{array}{c}S=k \mathrm{l}\mathrm{n}\mathsf{\Xi }+\beta \left(\langle E\rangle -{\displaystyle \sum _k}{\mu }_k\langle N_k\rangle \right)\end{array}$$

(1)

is maximized, yielding the mole fractions $x_k=\langle N_k\rangle /\langle N\rangle$ as the natural coordinates of the system.

2.1.2. Axiom 2 (Uniform Gibbs Free Energy with Spectral Minima)

The Gibbs free energy $G=\pm N_{A}h\nu$ is uniform across all phases, enforced by a non-vanishing lower bound $\delta >0$ on the spectrum of the helical operator (Lemma S4 [22,24]). This creates a stability window around the reference prime 19.

2.1.3. Axiom 3 (Irreducibility via Perpetual Bounded Oscillations)

Primes induce irreducible representations (Hilbert–Maschke theorem) of the helical flux, ensuring divergence-free flow

$$\begin{array}{c}\nabla ·{\displaystyle \sum _m}{\rho }_m{\mathbf{v}}_j=0\end{array}$$

(2)

and trigonometric boundedness (Lemmas S8 and S9 [23,24,25]). The background prime 19 is the fixed upper edge of the stability window.

These axioms uniquely select the non-proper Archimedean conical helix as the background geometry and force primes to act as indivisible bosonic modes with energies $E_p=\mathrm{l}\mathrm{n}p$ (Cycle-to-Prime Mapping, Sub-Lemma S9 [22,23,24]).

In physical terms, these three axioms encode the foundational principles of thermodynamics [23,24,25,26,27,28], gauge theory [29,30], and topology [31] as follows:

Axiom 1 (Strict Concave Entropy Maximization) is the variational statement of the Second Law: entropy is strictly concave, guaranteeing a unique global thermodynamic equilibrium and a single minimum of the grand potential [23,24,25,26,27,28].

Axiom 2 (Uniform Gibbs Free Energy with Spectral Minima) enforces equality of the molar Gibbs free energy across coexisting phases through a strict positive spectral gap $\delta >0$ in the helical operator — the molar analogue of the Einstein photoelectric threshold [31], ensuring no zero-energy excitations are allowed.

Axiom 3 (Irreducibility via Perpetual Bounded Oscillations) combines Noether’s theorem (exact conservation of helical flux) with topological irreducibility of modes. This forbids zero-helicity states and ensures perpetual bounded oscillations, realizing microscopically the QED prohibition of longitudinal photons [32,33] through the irreducible prime-length rotational symmetries of the helix.

Taken together, the axioms generate a grand potential $\omega (s)=\mathrm{l}\mathrm{n}Z(s)$ whose functional form deductively encodes the forbidden zero-helicity sector [23,24]. The grand-partition function $Z(s)$ (Section 3.2) therefore emerges as the universal generator of the entire theorem — the categorical crown jewel that converts the three axioms into concrete, zero-parameter predictions across thermodynamics, electromagnetism, and superconductivity [22,23,24].

2.2. Prime Numbers as Indivisible Cycle Lengths

A prime $p$ corresponds to an irreducible $p$-fold rotational symmetry that cannot be decomposed into smaller closed cycles. This follows directly from Axiom 3: the helical operator must act on a minimal, indivisible representation whose dimension/order is exactly $p$ (Lemma S8 [22,24]).

The angular quantization step encoded in every helical projection is

$$\begin{array}{c}\theta ={\displaystyle \frac{2\pi }{p}}\end{array}$$

(3)

This angle appears in the diagonal blocks of the star-topology transfer matrix $M(19-\mathrm{hub})$ as $\mathrm{c}\mathrm{o}\mathrm{s}(2\pi /p)$ and $\mathrm{s}\mathrm{i}\mathrm{n}(2\pi /p)$ for the peripheral modes, while the off-diagonal couplings to the 19-anchor preserve the overall periodicity governed by the prime [23,24].

2.3. Explicit Helical Operator and Effective Hessian PDE

2.3.1. ZMT Helical Operator

Axiom 2 requires a self-adjoint elliptic operator $H$ whose Gibbs functional

$$\begin{array}{c}G\left[\psi \right]=\int {\psi }^{*}H\psi d\mu \end{array}$$

has a strict positive spectral gap $\mathrm{i}\mathrm{n}\mathrm{f}\sigma (H)=\delta >0$ (Lemma S4 [22,23,24]). The unique operator satisfying these conditions is the Sturm–Liouville form

$$\begin{array}{c}H\psi =-{\displaystyle \frac{1}{w}}{\displaystyle \frac{d}{d\theta }}\left({\displaystyle \frac{1}{w}}{\displaystyle \frac{d\psi }{d\theta }}\right)+{\displaystyle \frac{{\alpha }^2}{r(\theta {)}^2}}\psi \end{array}$$

(4)

where the centrifugal term encodes the conical taper and $\alpha$ is fixed by the helicity constraint of Axiom 3. The Friedrichs extension ensures self-adjointness, and elliptic regularity plus the non-proper condition guarantee $\delta >0$ (Corollary S1 [22,24]).

2.3.2. Effective Hessian PDE

The natural Hessian along the helix is the differential part of $H$ [22,23,24]:

$$\begin{array}{c}({\mathrm{H}\mathrm{e}\mathrm{s}\mathrm{s}}_{\gamma }\varphi )\equiv -{\displaystyle \frac{1}{w}}{\displaystyle \frac{d}{d\theta }}\left({\displaystyle \frac{1}{w}}{\displaystyle \frac{d\varphi }{d\theta }}\right)+{\displaystyle \frac{{\alpha }^2}{r(\theta {)}^2}}\varphi \end{array}$$

(5)

Combining the entropy functional $S\left[\rho \right]=-\int \rho \mathrm{l}\mathrm{n}\rho d\mu$ (Axiom 1) and the Gibbs functional (Axiom 2) yields the phase functional

$$\begin{array}{c}F[\varphi ,\psi ]={\int }_X{e}^{-\psi }\mid {\mathrm{H}\mathrm{e}\mathrm{s}\mathrm{s}}_{\gamma }\varphi -C{\mid }^{2}d\mu \end{array}$$

(6)

where $C>0$ is the constant curvature floor enforced by Axiom 3 (Lemma S11 [22,23,24]).

Direct Euler–Lagrange variation with respect to $\varphi$ produces the stationarity condition

$$\begin{array}{c}{e}^{-\psi }({\mathrm{H}\mathrm{e}\mathrm{s}\mathrm{s}}_{\gamma }\varphi -C)=s(\theta )\end{array}$$

or explicitly

$$\begin{array}{c}{e}^{-\psi }\left[-{\displaystyle \frac{1}{w}}{\displaystyle \frac{d}{d\theta }}\left({\displaystyle \frac{1}{w}}{\displaystyle \frac{d\varphi }{d\theta }}\right)+{\displaystyle \frac{{\alpha }^2}{r(\theta {)}^2}}\varphi \right]=C+s(\theta )\end{array}$$

(7)

Here the source term is

$$\begin{array}{c}s(\theta )=e^{\psi }\left[\right[\nabla ,\nabla ]\mathrm{l}\mathrm{n}\mathsf{\Omega }+L_{\xi }({\partial }_{\theta }\varphi ·{\partial }_{\theta }\varphi )]\end{array}$$

(8)

with $\mathsf{\Omega }$ the entropy density, $\left[\nabla ,\nabla \right]$ the second covariant derivative along the helix, and $L_{\xi }$ the Lie derivative along the helical tangent vector. This effective Hessian PDE is a true theorem of ZMT (Corollary S1 [22,23,24]), derived solely from the three axioms and the explicit helical operator.

3. Variational Thermodynamics of Multi-Component Helical Systems and Categorical Invariance of $\mathit{Z}(\mathit{s})$

The Riemann Hypothesis is reinterpreted in the ZMT framework as a thermodynamic stability condition [22,23,24]: the equilibrium frequency spectrum (Lyapunov exponents) is minimized only when the thermodynamic parameter $s$ satisfies $\mathrm{R}\mathrm{e}(s)=1/2$ (strict concavity of the Gibbs functional $\mathsf{\Phi }\left[\psi \right]$ from Axiom 1). Any deviation off the line would violate the spectral gap $\delta >0$ (Axiom 2) or the divergence-free flux condition (Axiom 3). Thus, the physical grand-partition function must be constructed so that its dominant component always occupies more than half the total occupation when equilibrium is reached on the critical line (Theorem S26 [22,24]).

3.1. Grand-Partition Function (General Multi-Component Case)

The grand-partition function for an arbitrary multi-component system is constructed once from the helical gears present in the mixture [23,25]:

$$\begin{array}{c}Z(s)=\left[{\displaystyle \prod _{m=1}^M}{\displaystyle \prod _{k=2}^{p_m}}{\displaystyle \frac{1}{1-k^{-s}}}\right]·{\displaystyle \frac{1}{1-{19}^{-s}}}\end{array}$$

(9)

where $p_m$ is the deductive molecular (or elemental) prime of component $m$, and the product runs over all gears $k=\mathrm{2,3},\dots ,p_m$(each with multiplicity exactly 1) (Helical Gear Multiplicity Theorem S33 [22,24]).

3.2. Mole Fractions (General Multi-Component Case)

The occupation of each gear $k$ follows Bose–Einstein statistics [22,23,24]:

$$\begin{array}{c}{n}_k(s)={\displaystyle \frac{1}{k^s-1}}\end{array}$$

(10)

The global mole fraction of gear $k$ is therefore

$$\begin{array}{c}{x}_k(s)={\displaystyle \frac{n_k(s)}{n_{\mathrm{t}\mathrm{o}\mathrm{t}}(s)}}={\displaystyle \frac{1/(k^s-1)}{{\sum }_{j}1/(j^s-1)+1/({19}^s-1)}}\end{array}$$

(11)

3.3. Lyapunov Exponent per Gear

The Lyapunov exponent for each gear $k$ is obtained by equating the solved interaction parameter to the helical exponential-cosine form [22,23,24]:

$$\begin{array}{c}{\lambda }_{k,19}(x_k)=-{\displaystyle \frac{1}{\mid k-19\mid }}\mathrm{l}\mathrm{n}\left({\displaystyle \frac{\mid {\mathsf{\Delta }}_{k,19}(x_k)\mid }{\mid \mathrm{c}\mathrm{o}\mathrm{s}(2\pi \mid k-19\mid /(k·19))\mid }}\right)\end{array}$$

(12)

3.4. Interaction Parameter per Gear

The interaction parameter for each gear $k$ satisfies the linear first-order ODE that enforces the global minimum of the mixture frequency at the critical composition [22,23,24]:

$$\begin{array}{c}{\mathsf{\Delta }}_{k,19}(x_k)={\displaystyle \frac{-(k-19)x_k+C_k}{x_k(1-x_k)}}\end{array}$$

(13)

where the integration constant $C_k$ is fixed by the critical-point condition (Section 4.2).

3.5. Critical Composition per Gear ($x_{k,o}$) in the $s\to 0$ Limit

$$\begin{array}{c}{x}_{k,o}={\displaystyle \underset{s\to 0}{\mathrm{l}\mathrm{i}\mathrm{m}}}{x}_k(s)={\displaystyle \frac{1/\mathrm{l}\mathrm{n}k}{{\sum }_{j}1/\mathrm{l}\mathrm{n}j+1/\mathrm{l}\mathrm{n}19}}\end{array}$$

(14)

All quantities (Table 1) are fully deductive from the helical transfer matrix and the RH-imposed critical-line condition (Theorem S26 [22,23,24]).

3.6. Rigorous Physical Interpretation of Integer Gear-Up-to-Prime Rule

Under the integer gear-up-to-prime rule (Helical Gear Multiplicity Theorem S33 [22,24]), each component with deductive prime $p_m$ occupies exactly the gears $k=\mathrm{2,3},\dots ,p_m$ (each with multiplicity 1). This rule follows directly from Axiom 3: the helical operator must act on a minimal, indivisible representation whose dimension/order is exactly the prime $p_m$. Any smaller or larger set would either be reducible (violating irreducibility) or introduce extraneous modes (violating entropy maximization, Axiom 1).

3.6.1. Definition of a Gear

A gear is a discrete, indivisible helical winding mode labeled by an integer $k\ge 2$. Mathematically, each gear $k$ corresponds to an irreducible representation of the helical operator $H$ (Eq. (3)) with winding angle [22,23,24]

$$\begin{array}{c}{\theta }_k={\displaystyle \frac{2\pi }{k}}\end{array}$$

(15)

Physically, a gear represents a quantized collective excitation of the helical flux with a fixed angular periodicity. The occupation number of gear $k$ is the Bose–Einstein factor (Eq. 10)

$$\begin{array}{c}{n}_k(s)={\displaystyle \frac{1}{k^s-1}}\end{array}$$

which gives the average number of helical quanta in that mode at thermodynamic parameter $s$.

3.6.2. Integer Gear-Up-to-Prime Rule

For any component with deductive prime $p_m$, the allowed gears are exactly the integers

$$\begin{array}{c}k=\mathrm{2,3},\dots ,p_m (\mathrm{each} \mathrm{with} \mathrm{multiplicity} 1)\end{array}$$

(16)

This is the integer gear-up-to-prime rule (Helical Gear Multiplicity Theorem S33 [22,24]).

3.6.3. Integer Gear-Up-to-Prime Rule Axiomatic Deduction

The rule is not postulated; it is the unique outcome of the three axioms acting together on the helical representation graph.

• Axiom 3 (Irreducibility) requires that the helical operator acts on a minimal, indivisible representation whose dimension/order is exactly a prime number $p_m$. Any composite cycle length would be reducible, violating flux conservation and trigonometric boundedness (Lemma S8 [22,23,24]).
• Axiom 1 (Strict Concavity) demands the smallest possible set of modes that can still achieve a unique global minimum of the phase functional $F$. Including gears beyond $p_m$ would introduce redundant modes, increasing the variational cost without lowering the grand potential.
• Axiom 2 (Spectral Gap) guarantees a positive lower bound $\delta >0$ on all eigenvalues. The integer sequences up to the prime $p_m$ is the minimal set that satisfies both indivisibility and the spectral gap while allowing the observed stability window around prime 19.

Thus, the rule follows deductively: the helical operator must act on irreducible prime-length cycles (Axiom 3), the phase functional selects the smallest such set (Axiom 1), and the spectral gap is preserved only within that window (Axiom 2).

3.6.4. Physical Interpretation

• Each gear $k$ behaves as a quantized tooth on the helical driveshaft.
• The occupation $n_k(s)$ is the average number of quanta occupying that tooth at thermodynamic parameter $s$ [22,23,24].
• The mole fraction $x_k(s)$ is the fractional contribution of that gear to the total occupation [23,24,25].
• When a component has deductive prime $p_m$, it can only activate gears up to $p_m$. This is the microscopic origin of material specificity: different primes (different molecules/elements) activate different sets of helical modes, leading to distinct thermodynamic, magnetic, and transport properties [22,24].

3.7. Universal Property of $Z(s)$ in Category Theory

To prove that the grand-partition function $Z(s)$ satisfies the universal property required of the functorial backbone of the Zeta Minimizer Theorem [22,23,24].

3.7.1. Definition of Categories

Let ${\mathcal{C}}_{\mathrm{g}\mathrm{e}\mathrm{a}\mathrm{r}}$ be the category [34,35,36] whose:

Objects are multi-component helical gear systems equipped with the star-topology representation graph $\mathsf{\Gamma }$ [22,23,24] (anchored at prime 19) and integer-gear occupations $\left\{x_k(s)\right\}$.

Morphisms are natural transformations [24,37,38] that preserve the star topology [24,39,40], the integer-gear-up-to-prime rule, and the action of the helical operator $H$.

Let $\mathcal{T}\mathcal{h}\mathcal{e}\mathcal{r}\mathcal{m}$ be the category whose:

Objects are thermodynamic potentials [27] (real-valued functions of the reduced variable $s$).

Morphisms are natural transformations between these potentials that commute with differentiation with respect to $s$.

3.7.2. Definition of Functor

Define the functor

$$\mathcal{Z}:{\mathcal{C}}_{\mathrm{g}\mathrm{e}\mathrm{a}\mathrm{r}}\to \mathcal{T}\mathcal{h}\mathcal{e}\mathcal{r}\mathcal{m}$$

by sending each object $\left(\mathsf{\Gamma },\{x_k(s)\}\right)$ to the grand-partition function

$$Z(s)={\displaystyle \prod _k}{\displaystyle \frac{1}{1-k^{-s}}}·{\displaystyle \frac{1}{1-{19}^{-s}}}$$

and each morphism to the induced natural transformation on the potentials [24].

3.7.3. Universal Property Theorem

Theorem 1 [24]. $Z(s)$ (i.e., the functor $\mathcal{Z}$) satisfies the following universal property:

For any other functor $F:{\mathcal{C}}_{\mathrm{g}\mathrm{e}\mathrm{a}\mathrm{r}}\to \mathcal{T}\mathcal{h}\mathcal{e}\mathcal{r}\mathcal{m}$ that assigns to each helical gear system a thermodynamic potential compatible with the three axioms (strict concavity of the phase functional, positive spectral gap of $H$, and flux conservation), there exists a unique natural isomorphism

$$\eta :F\Rightarrow \mathcal{Z}$$

such that the diagram (Figure 1) commutes for every morphism in ${\mathcal{C}}_{\mathrm{g}\mathrm{e}\mathrm{a}\mathrm{r}}$.

Proof of Theorem 1.

• By Axiom 1 the phase functional $F$ is strictly concave and possesses a unique global minimizer.
• By Axioms 2 and 3 together with the Topology Selection Theorem (Theorem S4 [22,23,24]), the helical operator $H$ and the star-topology graph $\mathsf{\Gamma }$ are fixed uniquely and globally.
• The grand-partition function $Z(s)$ is constructed directly as the generating function of the integer-gear occupations on this unique $\mathsf{\Gamma }$. Any other assignment $F$ that respects the axioms must therefore reproduce the same occupation numbers and the same generating function [24,37,41].
• The functoriality of $Z(s)$ follows immediately: every morphism in ${\mathcal{C}}_{\mathrm{g}\mathrm{e}\mathrm{a}\mathrm{r}}$ (a natural transformation preserving the star topology and integer-gear rule) induces a unique compatible transformation on the potentials, and the diagram above commutes by construction [24,40,42].
• Uniqueness of the isomorphism $\eta$ follows from the strict concavity of $F$ (Axiom 1): any two functors agreeing on all objects and morphisms must coincide up to unique natural isomorphism [24,37,38].

Thus $Z(s)$ is the canonical (universal) functor from helical gear systems to thermodynamic potentials. This categorical universality guarantees that every derived quantity — Lyapunov exponents, interface matching conditions, conjugate pairs, electromagnetic fields, and superconducting shells — is representation-independent and deductively fixed once and for all by the three axioms.

3.8. Categorical Invariance of the Grand-Partition Function $Z(s)$— The Functorial Backbone of ZMT

The grand-partition function $Z(s)$ is not merely a formal generating function; it is the image of a functor $\mathcal{Z}$ [24,39,40] from the category of helical gear systems ${\mathcal{C}}_{\mathrm{g}\mathrm{e}\mathrm{a}\mathrm{r}}$ (whose objects are multi-component mixtures equipped with the star-topology representation graph $\mathsf{\Gamma }$ anchored at prime 19, and whose morphisms are natural transformations that preserve the integer-gear-up-to-prime rule) into the category of thermodynamic potentials [24,37,38,40,42]. This functorial invariance guarantees that every derived quantity — mole fractions $x_k(s)$, Lyapunov exponents ${\lambda }_{k,19}(x_k)$, interaction parameters ${\mathsf{\Delta }}_{k,19}(x_k)$, critical compositions, interface matching conditions, and thermodynamic conjugate pairs — is canonically defined and independent of arbitrary choices of basis or labeling. Without this categorical foundation, the three primitive ZMT axioms would be insufficient to produce a unique, consistent thermodynamic mapping. The universal invariance of $Z(s)$ is therefore the indispensable glue that elevates the axioms and helical geometry into a rigorous, representation-independent framework for all subsequent results.

3.8.1. Functorial Isomorphism Theorem Statement

Theorem 2 [24]. Let $X$ be any flux-like density satisfying Axiom 3 (irreducible bounded oscillations + divergence-free flux conservation [22,23,24]) and let $Y$ be its thermodynamically conjugate extent variable (i.e., the variable appearing with opposite sign in the differential of the Gibbs free energy) [43,44]. The phase functional $F[A,\psi ]$ (identified with the Gibbs free energy $G$) is strictly concave (Axiom 1). Therefore, the Legendre transform structure with respect to $X$ at fixed temperature (fixed $s=\mathrm{l}\mathrm{n}19/T$) yields the exact relation

$$\begin{array}{c}XY=Z(s)RT\end{array}$$

(17)

for any such conjugate pair $\left(X,Y\right)$. This is the Functorial Isomorphism Theorem (Theorem S28 [22,24]).

Proof of Theorem 2. The gear occupations $x_k(s)$ that define $Z(s)$ are determined solely by the helical operator $H$ and the star-topology representation graph $\mathsf{\Gamma }$ (Topology Selection Theorem S4 [22,23,24]). Because gear discretization commutes with all scaling transformations and because Axiom 3 forces the flux density $X$ to appear linearly in the source tensor $S$, the product $XY$ factors out of the phase functional $F$ exactly as $Z(s)RT$, independent of the specific physical identity of $X$.

3.8.2. Explicit Examples of the Universal Isomorphism

Because $Z(s)$ is the same grand-partition function in every case, the isomorphism (Eq. (17)) holds uniformly [24,27,41,45]:

• Mechanical pair (pressure–volume):

$$\begin{array}{c}PV=Z(s)RT\end{array}$$

(18)

• Magnetic pair (field strength–magnetization):

$$\begin{array}{c}HM=Z(s)RT\end{array}$$

(19)

• Electric pair (field strength–polarization):

$$\begin{array}{c}E{P}_L=Z(s)RT\end{array}$$

(20)

• Electrochemical pair (potential–charge density):

$$\begin{array}{c}\varphi q=Z(s)RT\end{array}$$

(21)

All four (and any future) conjugate pairs are governed by the identical functional $Z(s)$ [24].

3.8.3. Strict Validity of the Functorial Mapping of $Z(s)$

The grand-partition function $Z(s)$ is not a physical observable tied to any particular variable. It is the universal functor $\mathcal{Z}$ [24,42,45] that maps the helical mode occupations $\left\{x_k(s)\right\}$ of any multi-component system to the complete set of thermodynamic observables [24,35,36,38,39,42].

The phase functional $\mathcal{F}[A,\psi ]$ remains the single variational generator (Axiom 1 of the Zeta Minimizer Theorem). The helical operator $H$ and the representation graph $\mathsf{\Gamma }$ are fixed uniquely and globally by the three axioms together with the Topology Selection Theorem (Theorem S4 [22,23,24]).

Consequently, $Z(s)$ is an intrinsic property of the helical gearbox itself, independent of which flux density or conjugate variable is chosen for measurement. Different physical realizations (pressure, magnetic field, electric field, electrochemical potential, etc.) are merely distinct projections of the same underlying helical modes. The functorial isomorphism therefore guarantees that all projections share the identical, representation-independent currency $Z(s)$.

This categorical structure eliminates any possibility of category mismatch and ensures that every subsequent derivation — from Lyapunov exponents and interface matching to electromagnetic fields and superconducting shells — is canonically defined and internally consistent.

3.8.4. Experimental Anchor and Zero-Parameter Predictive Power

The universal grand-partition function $Z(s)$ is anchored to the pure prime-19 background sea. High-precision compressibility measurements confirm that the identical $Z(s)$ reproduces the experimental equation-of-state residuals $PV=Z(s)RT$ to sub-${10}^{-12}$ accuracy [23].

Because the same $Z(s)$ governs every thermodynamic conjugate pair, the Zeta Minimizer Theorem automatically yields zero-parameter predictions for magnetic susceptibility [46,47], electric polarization [48,49], electrochemical potentials [50,51,52], and superconducting critical fields [53,54]. Pressure–volume relations, magnetic response, electric response, and electrochemical behavior are therefore revealed as distinct projections of the identical underlying helical modes.

4. Solid–Fluid Interface Matching and Marginal Stability

The interface between solid and fluid phases is the critical location where the variational principles of the Zeta-Minimizer Theorem (ZMT) enforce thermodynamic and mechanical consistency. All derivations follow directly from the phase functional $F[A,\psi ]$, the grand-partition function $Z(s)$, and the three primitive axioms.

4.1. Grand-Potential Continuity Condition

At equilibrium the total Gibbs free energy $G$ (identified with the phase functional $F[A,\psi ]$ by Axiom 1) must be minimized globally. For two coexisting phases (solid and fluid) in contact, the grand potential $\omega (s)=\mathrm{l}\mathrm{n}Z(s)$ must be continuous across the interface up to a fixed per-gear constant determined by marginal stability [27,28,55].

Define the grand-potential difference between the solid and fluid phases as

$$\begin{array}{c}\mathsf{\Delta }\omega (s)={\omega }_{\mathrm{s}\mathrm{o}\mathrm{l}\mathrm{i}\mathrm{d}}(s)-{\omega }_{\mathrm{f}\mathrm{l}\mathrm{u}\mathrm{i}\mathrm{d}}(s)\end{array}$$

(22)

Thermodynamic equilibrium requires that any infinitesimal variation of the interface position leaves the total grand potential unchanged:

$$\begin{array}{c}\delta (\mathsf{\Delta }\omega )=0\end{array}$$

(23)

Because the grand-partition function $Z(s)$ has the same functional form in both phases (Eq. (9)), the only variational freedom at the boundary is the redistribution of gear occupations $x_k(s)$. The stationarity condition is therefore

$$\begin{array}{c}{\omega }_{\mathrm{s}\mathrm{o}\mathrm{l}\mathrm{i}\mathrm{d}}(s)={\omega }_{\mathrm{f}\mathrm{l}\mathrm{u}\mathrm{i}\mathrm{d}}(s)+{\displaystyle \sum _k}{C}_k\end{array}$$

(24)

where each constant $C_k$ is fixed once by the critical composition of gear $k$ (Interface Matching Theorem S45 [22,24]). This is the grand-potential continuity condition at the solid–fluid interface.

4.2. Derivation of Integration Constants $C_k$

The constants $C_k$ are obtained by enforcing the global minimum of the mixture frequency at the critical composition $x_{k,o}$ in the $s\to 0$ limit.

From the grand-partition function (Eq. (9)), the occupation of gear $k$ is given by Eq. (10):

$$\begin{array}{c}{n}_k\left(s\right)={\displaystyle \frac{1}{k^s-1}}\end{array}$$

In the limit $s\to 0^{+}$,

$$\begin{array}{c}{k}^s-1\approx s\ln k\hspace{0.25em}⟹\hspace{0.25em}{n}_k\left(s\right)\approx {\displaystyle \frac{1}{\mathrm{sln}k}}\end{array}$$

The total occupation behaves as

$$\begin{array}{c}{n}_{\mathrm{t}\mathrm{o}\mathrm{t}}\left(s\right)\approx {\displaystyle \frac{1}{s}}\left(\sum _j{\displaystyle \frac{1}{\ln j}}+{\displaystyle \frac{1}{\ln 19}}\right)\end{array}$$

The critical mole fraction of gear $k$ is therefore (Eq. (14))

$$\begin{array}{c}{x}_{k,o}=\underset{s\to 0}{\mathrm{l}\mathrm{i}\mathrm{m}}{x}_k(s)={\displaystyle \frac{1/\mathrm{l}\mathrm{n}k}{\sum _{j}1/\mathrm{l}\mathrm{n}j+1/\mathrm{l}\mathrm{n}19}}\end{array}$$

The interaction parameter for gear $k$ satisfies the linear first-order ODE (Eq. (13)) that enforces the global minimum of the mixture frequency at $x_{k,o}$:

$$\begin{array}{c}{\mathsf{\Delta }}_{k,19}\left(x_k\right)={\displaystyle \frac{-\left(k-19\right)x_k+C_k}{x_k\left(1-x_k\right)}}\end{array}$$

Integrating this ODE and imposing the critical-point condition at $x_k=x_{k,o}$ yields the exact deductive expression

$$\begin{array}{c}{C}_k=-{\displaystyle \frac{(k-19)x_{k,o}^2}{1-2x_{k,o}}}\end{array}$$

(25)

This completes the determination of the integration constant per gear (Interface Matching Theorem S45 [22,24]).

4.3. Marginal Stability at the Interface

The system reaches marginal stability when the least negative Lyapunov exponent among all gears reaches zero:

$$\begin{array}{c}\underset{k}{\mathrm{m}\mathrm{a}\mathrm{x}}{\lambda }_{k,19}(x_k)=0\end{array}$$

(26)

Substituting the interaction parameter (Eq. (13)) into the Lyapunov exponent formula (Eq. (12)) gives

$$\begin{array}{c}{\lambda }_{k,19}(x_k)=-{\displaystyle \frac{1}{\mid k-19\mid }}\mathrm{l}\mathrm{n}\left({\displaystyle \frac{\mid \left[\frac{-(k-19)x_k+C_k}{x_k(1-x_k)}\right]\mid }{\mid \mathrm{c}\mathrm{o}\mathrm{s}\left(2\pi \mid k-19\mid /(k·19)\right)\mid }}\right)\end{array}$$

(27)

At the solid–fluid interface the lowest gear ($k=2$) dominates because it has the strongest geometric leverage. Setting ${\lambda }_{\mathrm{2,19}}(x_2)=0$ enforces the marginal-stability condition that defines the superconducting or phase-transition shell at the boundary.

4.4. Role of the Covariant Fugacity Hessian

The continuity condition (Eq. (26)) is enforced by the stationarity of the phase functional $F[A,\psi ]$. The covariant fugacity Hessian ${\mathrm{H}\mathrm{e}\mathrm{s}\mathrm{s}}_{\gamma }$ (Corollary S1 [22,23,24]) projected onto the interface normal yields the effective boundary condition

$$\begin{array}{c}{e}^{-\psi }({\mathrm{H}\mathrm{e}\mathrm{s}\mathrm{s}}_{\gamma }\varphi -C)=s(r;x_k)\end{array}$$

(28)

where $s(r;x_k)$ is the source term (Eq. (8)). This equation, together with the grand-potential continuity (Eq. (22)), guarantees that both thermodynamic continuity and mechanical/electromagnetic consistency are satisfied simultaneously across the solid–fluid boundary.

5. Thermodynamic Conjugate Pairs and Variational Maxwell Relations

In the Zeta-Minimizer Theorem, all thermodynamic quantities arise as functorial shadows of the global minimization of the phase functional $F[A,\psi ]$ on the non-proper Archimedean conical helix (Topology Selection Theorem S4). This structure naturally generalizes the classical pressure–volume conjugate pair to arbitrary thermodynamic conjugate pairs.

5.1. General Conjugates

Let $X$ be any flux-like density satisfying Axiom 3 (divergence-free flux conservation) and let $Y$ be its thermodynamically conjugate extent variable. The phase functional $F$ is strictly concave (Axiom 1), so the Gibbs free energy $G\equiv F$ admits a Legendre transform at fixed temperature (fixed $s$) [23,24]:

$$\begin{array}{c}Y={\left({\displaystyle \frac{\partial G}{\partial X}}\right)}_s\end{array}$$

(29)

Because gear discretization commutes with scaling transformations, the product of any conjugate pair satisfies the natural isomorphism by Eq. (17)

$$\begin{array}{c}XY=Z(s)RT\end{array}$$

(Functorial Isomorphism Theorem 2). This identity holds invariantly once the helical primes and stable modes are fixed by the three axioms [22,23,24].

5.2. Classic Pressure–Volume Conjugate Pair $\left(P,V\right)$

The general structure maps directly to the familiar conjugate pair $\left(P,V\right)$, where $P$ is pressure (independent variable) and $V$ is molar volume (dependent variable) [22,24]. The entropy is

$$\begin{array}{c}S=R\mathrm{l}\mathrm{n}Z(s)-Rs{\displaystyle \frac{d\mathrm{l}\mathrm{n}Z}{ds}}\end{array}$$

(30)

The thermodynamic differential of the Gibbs free energy reads

$$\begin{array}{c}dG=-SdT+VdP\end{array}$$

(31)

5.3. Clapeyron Equation (Equilibrium Slope)

Along an equilibrium line at the global minimum of $G$ ($dG=0$) the two differentials must balance:

$$\begin{array}{c}-SdT+VdP=0\end{array}$$

(32)

Solving for the slope of the coexistence curve gives the classic Clapeyron equation:

$$\begin{array}{c}{\left({\displaystyle \frac{dP}{dT}}\right)}_{\mathrm{e}\mathrm{q}}={\displaystyle \frac{S}{V}}\end{array}$$

(33)

Substituting the ZMT expressions for $S$ and $V=Z\left(s\right) RT/P$ yields the explicit form in terms of the grand-partition function:

$$\begin{array}{c}{\left({\displaystyle \frac{dP}{dT}}\right)}_{\mathrm{e}\mathrm{q}}=\left[{\displaystyle \frac{R\mathrm{l}\mathrm{n}Z(s)-Rs\left(\frac{d\mathrm{l}\mathrm{n}Z}{ds}\right)}{Z(s)\left(RT/P\right)}}\right]={\displaystyle \frac{P}{T}}\left[{\displaystyle \frac{\mathrm{l}\mathrm{n}Z(s)}{Z(s)}}−{\displaystyle \frac{s}{Z(s)}}{\displaystyle \frac{dZ}{ds}}\right]\end{array}$$

(34)

This relation is fully deductive: it follows directly from the stationarity condition of the phase functional $F$(Axiom 1) and the conjugate-pair isomorphism (Eq. (17)).

5.4. Variational Maxwell Relations

The covariant fugacity Hessian ${\mathrm{H}\mathrm{e}\mathrm{s}\mathrm{s}}_{\gamma }$ is symmetric by construction (self-adjoint helical operator from Axiom 2). Therefore, the mixed second derivatives of the phase functional commute, and the thermodynamic potentials satisfy the variational Maxwell relations [22,24]:

$$\begin{array}{c}{\left({\displaystyle \frac{\partial Y_i}{\partial X_j}}\right)}_s={\left({\displaystyle \frac{\partial Y_j}{\partial X_i}}\right)}_s\end{array}$$

(35)

for any conjugate pair $\left(X_i,Y_i\right)$. These relations hold for magnetic, electric, and all other pairs (Variational Maxwell Relations Theorem S37 [22,24]).

All results in this section follow deductively from the three axioms of the Zeta-Minimizer Theorem and the functorial invariance of $Z(s)$ with no additional postulates.

5.5. Thermodynamic Potentials and Heat Capacities

All thermodynamic potentials and heat capacities are generated from the same grand-partition function $Z(s)$ and therefore share identical functional forms for both magnetic and electric conjugate pairs [24,47,56,57]. The explicit expressions (derived consistently from the entropy (Eq. (30)) and the Gibbs free energy $G=-RT\mathrm{l}\mathrm{n}Z(s)$) are collected in Table 2.

6. Variational Emergence of Classical Phase Equilibria in Non-Reactive Multi-Component Systems

6.1. Phase Functional Extensivity in Non-Reactive Multi-Component Systems Theorem

Theorem 3. In the ZMT framework, the phase functional of a non-reactive multi-component system satisfies exact extensivity:

Molar (intensive) phase functional

$$\begin{array}{c}G=-RT\ln Z\left(s\right)\end{array}$$

(36)

Total (extensive) phase functional

$$\begin{array}{c}\underset{‾}{G}=n_{\mathrm{total}}G=-n_{\mathrm{total}}RT\ln Z\left(s\right)\end{array}$$

(37)

where $n_{\mathrm{total}}={\sum }_i{n}_i$ and the molar grand-partition function is

$$\begin{array}{c}Z\left(s\right)=Z_{\mathrm{vacuum}}\left(s\right){\displaystyle \prod _i}{Z}_i{\left(s\right)}^{x_i}\end{array}$$

(38)

with $x_i=n_i/n_{\mathrm{total}}$ the macroscopic mole fraction of component $i$ in the mixture.

Definitions

• $G$: molar Gibbs free energy of the mixture (per mole).
• $\underset{‾}{G}$: total Gibbs free energy of the macroscopic phase.
• $Z\left(s\right)$: molar grand-partition function of the phase.
• $Z_i\left(s\right)$: contribution from the helical gears exclusive to component $i$.
• $Z_{\mathrm{vacuum}}\left(s\right)$: contribution from the prime-19 master hub (background vacuum).
• Non-reactive: mole numbers $\left\{n_i\right\}$ are independent (no stoichiometric constraint).

Proof of Theorem 3.

• Gear assignment: By Axiom 3 and the integer gear-up-to-prime rule, each species $i$ is assigned a unique prime identifier and a fixed, indivisible set of helical gears with contribution $Z_i\left(s\right)$.
• Occupational independence: In the absence of chemical reaction, the gears of different species are occupationally independent.
• Star-topology product structure: The prime-19 star topology (master hub + independent spokes) implies that the grand-partition function of the mixture factors as

  $$\begin{array}{c}{Z}_{\mathrm{total}}\left(s\right)\propto Z_{\mathrm{vacuum}}\left(s\right){\displaystyle \prod _i}{\mathsf{\Xi }}_i\end{array}$$
• Molar normalization: Normalizing to one mole of mixture and expressing each species contribution per mole via its mole fraction yields

  $$\begin{array}{c}Z\left(s\right)=Z_{\mathrm{vacuum}}\left(s\right){\displaystyle \prod _i}{Z}_i{\left(s\right)}^{x_i}\end{array}$$
  Phase-functional identification: By the foundational ZMT map, the molar phase functional (Eq. (36)) is

  $$\begin{array}{c}F=G=-RT\ln Z\left(s\right)\end{array}$$
• Extensivity: For a macroscopic system of $n_{\mathrm{total}}$ moles, the total phase functional scales linearly (Eq. (37)):

  $$\begin{array}{c}\underset{‾}{G}=n_{\mathrm{total}}G=-n_{\mathrm{total}}RT\ln Z\left(s\right)\end{array}$$

This completes the proof. Every step follows directly from the helical gear construction, star topology, and species independence in the non-reactive limit. □

Corollary 1. The phase functional is homogeneous of degree 1 in the particle numbers, satisfying the fundamental requirement for a consistent thermodynamic potential in extensive systems.

6.2. Emergence of Chemical Potential from the Phase Functional Theorem

Theorem 4. In a non-reactive multi-component phase [18,19,20,21] within the Zeta-Minimizer framework, the chemical potential of component $i$ is obtained exactly from the total phase functional as

$$\begin{array}{c}{\mu }_i={\left({\displaystyle \frac{\partial \underset{‾}{G}}{\partial n_i}}\right)}_{T,P,n_{j\ne i}}=G+RT\ln \left(x_i{Z}_i\left(s\right)\right)+\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}\end{array}$$

(40)

where $\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}$ is the species-specific correction arising from the gas–solid interface matching constants $C_k$.

Definitions

• $\underset{‾}{G}$: total (extensive) phase functional of the system.
• $G=-RT\ln Z\left(s\right)$: molar phase functional.
• $x_i=n_i/n_{\mathrm{total}}$: macroscopic mole fraction of component $i$ in the mixture.
• $Z_i\left(s\right)$: contribution to the grand-partition function from the helical gears of component $i$.
• $\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}$: deductive shift originating from the interface matching constants $C_k$ of component $i$ (enforced by grand-potential continuity).

Proof of Theorem 4.

1.

Total phase functional: From Theorem 3 (Eq. (37)),

$$\begin{array}{c}\underset{‾}{G}=n_{\mathrm{total}}G=-n_{\mathrm{total}}RT\ln Z\left(s\right)\end{array}$$

2.

Definition of chemical potential

$$\begin{array}{c}{\mu }_i={\left({\displaystyle \frac{\partial \underset{‾}{G}}{\partial n_i}}\right)}_{T,P,n_{j\ne i}}\end{array}$$

Differentiating at constant $n_{j\ne i}$ gives $\partial n_{\mathrm{total}}/\partial n_i=1$, so

$$\begin{array}{cc}& {\mu }_i=G+n_{\mathrm{total}}{\left({\displaystyle \frac{\partial G}{\partial n_i}}\right)}_{T,P,n_{j\ne i}}\end{array}$$

3.

Derivative of the molar functional: Since $G=-RT\ln Z\left(s\right)$ (Eq. (36)),

$$\begin{array}{c}{\displaystyle \frac{\partial G}{\partial n_i}}=-RT·{\displaystyle \frac{1}{Z\left(s\right)}}·{\displaystyle \frac{\partial Z\left(s\right)}{\partial n_i}}\end{array}$$

4.

Explicit dependence of $Z\left(s\right)$ on $n_i$: From the product structure (Eq. (38))

$$\begin{array}{c}Z\left(s\right)=Z_{\mathrm{vacuum}}\left(s\right){\displaystyle \prod _i}{Z}_i{\left(s\right)}^{x_i}\end{array}$$

taking the logarithmic derivative at constant $n_{j\ne i}$ (only $x_i$ changes) and using $\partial x_i/\partial n_i=\left(1-x_i\right)/n_{\mathrm{total}}$ together with the normalization ${\sum }_j{x}_j=1$, the cross terms cancel exactly, yielding

$$\begin{array}{c}{n}_{\mathrm{total}}{\left({\displaystyle \frac{\partial \ln Z\left(s\right)}{\partial n_i}}\right)}_{T,P,n_{j\ne i}}=\ln Z_i\left(s\right)+\ln x_i-\ln Z\left(s\right)\end{array}$$

5.

Final assembly: Substituting back and collecting terms gives

$$\begin{array}{c}{\mu }_i=G+RT\left[\mathrm{l}\mathrm{n}\left(x_i{Z}_i\left(s\right)\right)\right]+\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}\end{array}$$

where the interface correction $\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}$ enters directly from grand-potential continuity at the solid–fluid boundary (Section 4).

This completes the proof. All steps follow strictly from extensivity (Theorem 3), the product structure of $Z\left(s\right)$, and interface matching.

Corollary 2. The factor $Z_i\left(s\right)$ supplies the exact analogue of the classical fugacity coefficient ${\varphi }_i$, while the constants $C_k$ provide the deductive counterpart of the reference-state chemical-potential shift [18,19,20,21]. Thus both the activity coefficient and the standard-state correction [18,19,20,21] emerge naturally and parameter-free from the ZMT construction.

6.3. Vapor-Liquid Equilibrium and Classical Fugacity/Activity Coefficients Emergence Theorem

Theorem 5. In the Zeta-Minimizer framework, vapor-liquid equilibrium (VLE) emerges as the equality of chemical potentials between coexisting phases [18,19,20,21]. The helical fugacity factor $Z_i\left(s\right)$ and the interface matching constants $C_k$ recover the classical fugacity coefficient ${\varphi }_i$ and activity coefficient ${\gamma }_i$ (including the reference-state correction) without additional postulates.

At equilibrium between a vapor phase (V) and a liquid phase (L), the chemical potential of every component $i$ is equal:

$${\mu }_i^{\left(\mathrm{V}\right)}={\mu }_i^{\left(\mathrm{L}\right)}$$

(41)

This equality is equivalent to the classical fugacity relation [18,19,20,21]

$$y_i{\varphi }_i^{\left(\mathrm{V}\right)}P=x_i{\gamma }_i^{\left(\mathrm{L}\right)}{f}_i^{\circ }$$

(42)

where ${\varphi }_i^{\left(\mathrm{V}\right)}=Z_i^{\left(\mathrm{V}\right)}\left(s\right)$ is the vapor-phase fugacity coefficient, and both ${\gamma }_i^{\left(\mathrm{L}\right)}$and the reference fugacity $f_i^{\circ }$ arise from the difference in phase-specific helical gear structures and the interface constants $C_k$.

Proof of Theorem 5.

• Chemical potential in any phase: From Theorem 4, the chemical potential of component $i$ in a given phase (Eq. (40)) is

  $${\mu }_i=G+RT\ln \left(x_i{Z}_i\left(s\right)\right)+\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}$$

  where $G$ is the molar phase functional of that phase, $Z_i\left(s\right)$ is the helical-gear contribution of species $i$ in that phase, and $\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}$ is the interface correction.
• Equilibrium condition: At VLE the total phase functional $\underset{‾}{G}$ of the combined system is stationary with respect to matter transfer between phases. Hence

  $${\mu }_i^{\left(\mathrm{V}\right)}={\mu }_i^{\left(\mathrm{L}\right)}\forall i$$
• Substitution: Inserting the expression from Step 1 for both phases (using $y_i$ for the vapor mole fraction and $x_i$ for the liquid) yields

  $$G^{\left(\mathrm{V}\right)}+RT\mathrm{l}\mathrm{n}\left(y_i{Z}_i^{\left(\mathrm{V}\right)}\left(s\right)\right)+\mathsf{\Delta }{\mu }_i^{\left(C_k,V\right)}=G^{\left(\mathrm{L}\right)}+RT\mathrm{l}\mathrm{n}\left(x_i{Z}_i^{\left(\mathrm{L}\right)}\left(s\right)\right)+\mathsf{\Delta }{\mu }_i^{\left(C_k,L\right)}$$
• Rearrangement: Collecting terms gives

  $$RT\ln \left({\displaystyle \frac{y_i{Z}_i^{\left(\mathrm{V}\right)}\left(s\right)}{x_i{Z}_i^{\left(\mathrm{L}\right)}\left(s\right)}}\right)=G^{\left(\mathrm{L}\right)}-G^{\left(\mathrm{V}\right)}+\left(\mathsf{\Delta }{\mu }_i^{\left(C_k,L\right)}-\mathsf{\Delta }{\mu }_i^{\left(C_k,V\right)}\right).$$
• Identification of classical coefficients: The left-hand side contains the composition and helical fugacity factors. Identifying ${\varphi }_i^{\left(\mathrm{V}\right)}=Z_i^{\left(\mathrm{V}\right)}\left(s\right)$. The right-hand side is the difference in molar phase functionals plus the difference in interface corrections, which is exactly $RT\ln {\gamma }_i^{\left(\mathrm{L}\right)}+RT\mathrm{l}\mathrm{n}\left(f_i^{\circ }/P\right)$Exponentiating and multiplying through by $x_{i}P$ recovers the classical VLE condition [18,19,20,21]

  $$y_i{\varphi }_i^{\left(\mathrm{V}\right)}P=x_i{\gamma }_i^{\left(\mathrm{L}\right)}{f}_i^{\circ }$$

  (43)

Corollary 3. The helical factor $Z_i\left(s\right)$ is the direct deductive analogue of the fugacity coefficient ${\varphi }_i$. The phase-specific gear structures together with the interface matching constants $C_k$ supply the origin of the activity coefficient ${\gamma }_i$ and the reference fugacity $f_i^{\circ }$. Classical VLE thermodynamics [18,19,20,21] is therefore an emergent consequence of variational optimization on the helical manifold.

6.4. General Multi-Phase Multi-Component Equilibrium Emergence Theorem

Theorem 6. In the Zeta-Minimizer framework, thermodynamic equilibrium [18,19,20,21] in a system of $\zeta$ phases and $c$ components (non-reactive or reactive) is completely characterized by the equality of the chemical potential of each component across all phases:

$${\mu }_i^{\left(\alpha \right)}={\mu }_i^{\left(\beta \right)}=\dots ={\mu }_i^{\left(\zeta \right)}\forall i=1,\dots ,c,\forall \alpha ,\beta ,\dots ,\zeta$$

(44)

This single condition recovers all classical equilibrium criteria (VLE, LLE, SLE, VLLE, and chemical-reaction equilibrium).

Definitions

• $G^{\left(\alpha \right)}=-RT\ln Z^{\left(\alpha \right)}\left(s\right)$: molar phase functional of phase $\alpha$.
• $\underset{‾}{G}$: total extensive phase functional of the multi-phase system.
• ${\mu }_i^{\left(\alpha \right)}$: chemical potential of component $i$ in phase $\alpha$.
• $Z_i^{\left(\alpha \right)}\left(s\right)$: helical-gear contribution of component $i$ in phase $\alpha$.
• $C_k^{\left(\alpha \right)}$: interface matching constants of component $i$ in phase $\alpha$.
• $\zeta$: number of coexisting phases.
• $c$: number of components.

Proof of Theorem 6.

• Phase functional per phase: By Theorem 3, each phase $\alpha$ possesses its own molar phase functional

  $$G^{\left(\alpha \right)}=-RT\ln Z^{\left(\alpha \right)}\left(s\right)$$

  where $Z^{\left(\alpha \right)}\left(s\right)$ is constructed from the helical gears present in that phase.
• Chemical potential per phase: From Theorem 4, the chemical potential of component $i$ in phase $\alpha$ is

  $${\mu }_i^{\left(\alpha \right)}=G^{\left(\alpha \right)}+RT\mathrm{l}\mathrm{n}\left(x_i^{\left(\alpha \right)}{Z}_i^{\left(\alpha \right)}\left(s\right)\right)+\mathsf{\Delta }{\mu }_i^{\left(C_k,\alpha \right)}$$
• Global stationarity: Equilibrium minimizes the total phase functional $\underset{‾}{G}$. Any virtual transfer of component $i$ between phases must satisfy $\delta \underset{‾}{G}=0$, which requires

  $${\mu }_i^{\left(\alpha \right)}={\mu }_i^{\left(\beta \right)}\mathrm{for} \mathrm{every} \mathrm{pair} \mathrm{of} \mathrm{phases} \alpha ,\beta$$

  and for every component $i$
• Equality across all phases: Extending the pairwise condition to the full set of $\zeta$ phases gives

  $${\mu }_i^{\left(1\right)}={\mu }_i^{\left(2\right)}=\dots ={\mu }_i^{\left(\zeta \right)}\forall i$$
• Recovery of classical conditions: Substitution of the expression from Step 2 into the equality and rearrangement directly yields:
  • For phase equilibria (VLE, LLE, SLE, …): $y_i{\varphi }_{i}P=x_i{\gamma }_i{f}_i^{\circ }$, where ${\varphi }_i=Z_i^{\left(\mathrm{V}\right)}\left(s\right)$ and ${\gamma }_i$, $f_i^{\circ }$ originate from differences in $Z_i\left(s\right)$ and $C_k$ between phases.
  • For chemical reactions: the same ${\mu }_i$-equality combined with the stoichiometry matrix produces the equilibrium constant

    $$K=\prod _i{a}_i^{{\nu }_i}$$

Corollary 4. The helical factor $Z_i^{\left(\alpha \right)}\left(s\right)$ is the direct origin of the classical fugacity coefficient ${\varphi }_i$ in phase $\alpha$ [18,19,20,21]. Differences in gear occupation and interface constants $C_k$ across phases generate the activity coefficient ${\gamma }_i$and the reference fugacity $f_i^{\circ }$.

Corollary 5. $c\left(\zeta -1\right)$ independent chemical-potential equalities, together with the constraints of temperature, pressure, and stoichiometry, recover the classical Gibbs phase rule exactly [18,19,20,21].

Corollary 6. All thermodynamic equilibria [18,19,20,21]—phase and chemical—are unified under the single variational principle of minimizing the total phase functional $\underset{‾}{G}$ with respect to component distribution across phases and reaction extents.

7. Variational Emergence of Chemical Equilibrium in Reactive Systems

7.1. Extent of Reaction as a Variational Parameter

An independent chemical reaction is written in the standard form [18,19,20,21]

$$\sum _i{\nu }_i{\mathrm{React}}_i\rightleftharpoons \sum _j{\nu }_j{\mathrm{Prod}}_j$$

(45)

where ${\nu }_i<0$ for reactants and ${\nu }_j>0$ for products.

The extent of reaction $\xi$ parameterizes the instantaneous composition. The mole change of each species is

$$\mathsf{\Delta }{n}_j={\nu }_j\xi$$

and the mole fractions become

$$x_i\left(\xi \right)=x_i\left(0\right)+{\displaystyle \frac{{\nu }_i\xi }{n_{\mathrm{tot}}}}$$

The grand-partition function of the reacting mixture is therefore a function of $\mathsf{\xi }$ through the collective gear occupations.

7.2. Global Minimization of the Phase Functional Determines ${\xi }_{\mathrm{eq}}$

By Axiom 1, $F$ is strictly concave in $\xi$. Consequently it possesses a unique global minimizer ${\xi }_{\mathrm{eq}}$ satisfying

$${{\displaystyle \frac{\partial F}{\partial \xi }}∣}_{{\xi }_{\mathrm{eq}}}=0$$

(47)

This single stationarity condition simultaneously determines:

• the equilibrium composition at the gas–solid boundary,
• the consumption of each reactant and formation of each product,
• the final gear union $K_{\mathrm{union}}\left({\xi }_{\mathrm{eq}}\right)$

Any species may be chosen as reference; its change $\mathsf{\Delta }{n}_{\mathrm{ref}}={\nu }_{\mathrm{ref}}\xi$ fixes $\xi$, and stoichiometry then determines all other $\mathsf{\Delta }{n}_j$. The global minimization of $F$ automatically couples all gear occupations through the star-topology source tensor.

7.3. Extent of Reaction as Natural Constraint in ZMT

The star topology implies that every gear interacts only through the prime-19 anchor. Global minimization of $F$ nevertheless imposes a collective constraint across the entire mixture. The scalar $\xi$ is the single variational parameter that re-allocates gears along the atomic shaft while automatically satisfying the interface matching condition (Eq. (24))

$${\omega }_{\mathrm{solid}}-{\omega }_{\mathrm{gas}}=\sum _k{C}_k$$

at $\xi ={\xi }_{\mathrm{eq}}$.

Thus the classical stoichiometric constraint [18,19,20,21] emerges directly from the variational structure; it is not imposed externally.

7.4. Single-Reaction Summary

For any single independent reaction, the extent $\xi$ may be defined with respect to any participating species. All molar changes are then fixed by the stoichiometric coefficients. This constraint is already encoded in the global minimization of the phase functional $F$ and the grand-partition function at the boundary. No additional postulate is required.

7.5. Generalization to Multiple (Coupled) Reactions

7.5.1. Stoichiometry Matrix and Independent Reactions

Consider $N$ species and $R$ candidate reactions. The stoichiometry matrix $\nu$ (size $N\times R$) has columns that are the stoichiometric vectors of each reaction. The number of independent reactions is

$$R_{\mathrm{ind}}=\mathrm{rank}\left(\nu \right)$$

(48)

Any set of $R_{\mathrm{ind}}$ linearly independent columns forms a basis for the reaction space.

7.5.2. Extent-Vector Parameterization and Mole-Fraction Evolution

Define the extent vector $\mathit{\xi }={\left({\xi }_1,{\xi }_2,\dots {\xi }_{R_{\mathrm{ind}}}\right)}^T$ The change in moles of species $i$ is

$$\mathsf{\Delta }{n}_j={\displaystyle \sum _{r=1}^{R_{\mathrm{ind}}}}{\nu }_{i,r}{\xi }_r$$

The gear occupations and grand-partition function then depend on $\xi$:

$$x_i\left(\xi \right)=x_i\left(0\right)+{\displaystyle \frac{1}{n_{\mathrm{tot}}}}\sum _r{\nu }_{i,r}{\xi }_r$$

(49)

Equilibrium occurs at the global minimum of $F$ with respect to $\xi$:

$${{\displaystyle \frac{\partial F}{\partial {\xi }_r}}∣}_{{\xi }_{\mathrm{eq}}}=0,r=1,\dots ,R_{\mathrm{ind}}$$

(50)

7.5.3. Equilibrium Gas Composition

The equilibrium composition $x_i^{\mathrm{gas}}\left({\xi }_{\mathrm{eq}}\right)$ satisfies both the interface matching condition ${\omega }_{\mathrm{solid}}-{\omega }_{\mathrm{gas}}={\sum }_k{C}_k$ and the stationarity of $F_{\mathrm{mix}}$. Linear dependence among reactions is automatically eliminated by restricting to the $R_{\mathrm{ind}}$ independent extents.

7.5.4. General Variational Methodology for Multi-Reaction Systems

• Construct the stoichiometry matrix $\nu$and compute $R_{\mathrm{ind}}=\mathrm{rank}\left(\nu \right)$.
• Parameterize the system by the extent vector $\xi \in {\mathbb{R}}^{R_{\mathrm{ind}}}$.
• Express all gear occupations $x_k\left(\xi \right)$ and the phase functional $F\left[A,\psi ,\xi \right]$.
• Solve the system $\partial F/\partial {\xi }_r=0$ ($r=1,\dots ,R_{\mathrm{ind}}$) for ${\xi }_{\mathrm{eq}}$.
• Obtain the equilibrium gas composition from $x_i\left({\xi }_{\mathrm{eq}}\right)$ and the interface matching condition.

This procedure is fully deductive: multiple reactions are treated as a single global minimization of $F$ on the star-topology shaft, with the rank condition ensuring no redundancy. The framework extends directly to biochemical networks and any multi-reaction system.

7.6. Phase Functional Extensivity in Reactive Systems Theorem

Theorem 7. In the Zeta-Minimizer framework, for a reactive multi-component system in which the species mole numbers are linked through independent extents of reaction $\left\{{\xi }_r\right\}$ ($r=1,\dots ,R$), the phase functional satisfies exact extensivity:

1.

Molar (intensive) phase functional (Eq. (36)),

$$G=-RT\ln Z\left(s;\left\{{\xi }_r\right\}\right)$$

(51)

2.

Total (extensive) phase functional (Eq. (37)),

$$\underset{‾}{G}=n_{\mathrm{total}}G=-n_{\mathrm{total}}RT\ln Z\left(s;\left\{{\xi }_r\right\}\right)$$

(52)

Here $Z\left(s;\left\{{\xi }_r\right\}\right)$ is the molar grand-partition function, which depends parametrically on the reaction extents through the stoichiometric constraints on the mole numbers.

Definitions

• $G$: molar phase functional of the reactive mixture.
• $\underset{‾}{G}$: total extensive phase functional of the macroscopic reactive system.
• ${\xi }_r$: extent of the $r$-th independent reaction (variational parameter).
• $n_i=n_i^0+{\sum }_r{\nu }_{i,r}{\xi }_r$: mole number of component $i$, constrained by stoichiometry.
• $n_{\mathrm{total}}={\sum }_i{n}_i$: total number of moles in the phase.
• $Z\left(s;\left\{{\xi }_r\right\}\right)$: molar grand-partition function of the reactive mixture, dependent on $s$ and the extents $\left\{{\xi }_r\right\}$ via the mole fractions $x_i\left(\left\{{\xi }_r\right\}\right)$

Proof of Theorem 7.

• Stoichiometric constraint: By the stoichiometry matrix, the mole numbers are linked:

  $$n_i=n_i^0+{\displaystyle \sum _{r=1}^R}{\nu }_{i,r}{\xi }_r\forall i$$
• Grand-partition function dependence: The grand-partition function is constructed from the helical-gear contributions of all species. Because the macroscopic mole fraction $x_i$ depend on the mole numbers $n_i\left(\left\{{\xi }_r\right\}\right)$

  $$Z\left(s\right)=Z\left(s;\left\{{\xi }_r\right\}\right)$$

  where the dependence enters exclusively through the mole fractions $x_i\left(\left\{{\xi }_r\right\}\right)$ in the product structure.

• Molar phase functional: By the foundational ZMT identification,

  $$G=-RT\ln Z\left(s;\left\{{\xi }_r\right\}\right)$$
• Extensivity: For a macroscopic system of $n_{\mathrm{total}}$ moles, the total phase functional scales linearly:

  $$\underset{‾}{G}=n_{\mathrm{total}}G=-n_{\mathrm{total}}RT\ln Z\left(s;\left\{{\xi }_r\right\}\right)$$
  The stoichiometric subspace preserves homogeneity of degree 1 in $n_{\mathrm{total}}$
• Variational consistency: The extents $\left\{{\xi }_r\right\}$ are themselves variational parameters. Minimization of $\underset{‾}{G}$ is performed jointly with respect to the overall system size $n_{\mathrm{total}}$ and the reaction coordinates $\left\{{\xi }_r\right\}$, maintaining full thermodynamic consistency.

Corollary 7. The molar phase functional $G$ remains a well-defined state function of $T$, $P$, and the reaction extents $\left\{{\xi }_r\right\}$

Corollary 8. The chemical potential of component $i$ is still given by ${\mu }_i=\left(\partial \underset{‾}{G}/\partial n_i\right)$, with the derivative taken along the allowed stoichiometric subspace defined by the reaction extents.

Corollary 9. This theorem supplies the foundation for the derivation of chemical equilibrium constants and multi-phase reactive equilibria in subsequent results.

7.7. Chemical Potential in Reactive Systems Theorem

Theorem 8. In the Zeta-Minimizer framework, for a reactive multi-component system in which species mole numbers are linked through independent extents of reaction $\left\{{\xi }_r\right\}$, the chemical potential of each component is rigorously defined as the partial derivative of the total phase functional taken along the stoichiometric subspace.

Statement:

$${\mu }_i={\left({\displaystyle \frac{\partial \underset{‾}{G}}{\partial n_i}}\right)}_{T,P,\left\{{\xi }_r\right\}}$$

(53)

where the derivative respects the stoichiometric constraints, and evaluates to

$${\mu }_i=G+RT\ln \left(x_i{Z}_i\left(s\right)\right)+\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}$$

(53)

Here $G$ is the molar phase functional of the reactive mixture and the derivative is performed along the allowed stoichiometric subspace.

Definitions

• $\underset{‾}{G}$: total extensive phase functional of the reactive system.
• $G=-RT\ln Z\left(s;\left\{{\xi }_r\right\}\right)$: molar phase functional of the reactive mixture.
• $x_i=n_i/n_{\mathrm{total}}$: macroscopic mole fraction of component $i$ in the mixture.
• $\left\{{\xi }_r\right\}$: vector of independent reaction extents (variational parameters).
• Stoichiometric subspace: the manifold in composition space defined by $n_i=n_i^0+{\sum }_r{\nu }_{i,r}{\xi }_r$

Proof of Theorem 8.

1.

Total phase functional: From Theorem 7,

$$\underset{‾}{G}=n_{\mathrm{total}}G=-n_{\mathrm{total}}RT\ln Z\left(s;\left\{{\xi }_r\right\}\right)$$

with dependence on $\left\{{\xi }_r\right\}$ entering through the mole numbers $n_i\left(\left\{{\xi }_r\right\}\right)$.

2.

Variational definition: The chemical potential is

$${\mu }_i={\left({\displaystyle \frac{\partial \underset{‾}{G}}{\partial n_i}}\right)}_{T,P,\left\{{\xi }_r\right\}}$$

where the derivative is taken along the stoichiometric subspace (i.e., when $n_i$ changes, the other $n_j$ adjust according to the stoichiometry matrix so the variation remains in the allowed reaction manifold).

• Chain rule:

  $${\mu }_i=G+n_{\mathrm{total}}{\left({\displaystyle \frac{\partial G}{\partial n_i}}\right)}_{T,P,\left\{{\xi }_r\right\}}$$
• Derivative of the molar functional: Since $G=-RT\ln Z\left(s;\left\{{\xi }_r\right\}\right)$, differentiation along the stoichiometric subspace using the product structure of $Z\left(s\right)$ and the definition of mole fractions yields, after cancellation of cross terms,

  $${\mu }_i=G+RT\ln \left(x_i{Z}_i\left(s\right)\right)+\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}$$

  where $\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}$ is the interface matching contribution.
• Consistency: The derivative is confined to the stoichiometric subspace defined by the matrix $\nu$. This ensures the chemical potential remains well-defined even though the individual $n_i$ are not independent.

Corollary 10. The chemical potential retains the classical form ${\mu }_i=G+RT\ln a_i+$correction, where the activity $a_i$ incorporates the helical fugacity factor $Z_i\left(s\right)$ and the mole fraction $x_i$

Corollary 11. The reaction affinity of the $r$-th independent reaction is

$$A_r=\sum _i{\nu }_{i,r}{\mu }_i$$

which supplies the variational driving force for each reaction.

Corollary 12. When all ${\xi }_r=0$, the expression reduces exactly to the non-reactive chemical potential of Theorem 4, providing a continuous bridge between the two regimes.

7.8. General Multi-Reaction Multi-Phase Equilibrium Emergence Theorem

Theorem 9. In the Zeta-Minimizer framework, complete thermodynamic equilibrium in a multi-phase, multi-component, multi-reaction system is achieved when two conditions hold simultaneously: (i) the chemical potential of every component is identical in all coexisting phases, and (ii) the affinity of every independent reaction vanishes. Both conditions emerge directly from the minimization of the total phase functional.

Statement:

At equilibrium:

1.

Phase equilibrium

$${\mu }_i^{\left(\alpha \right)}={\mu }_i^{\left(\beta \right)}=\dots ={\mu }_i^{\left(\zeta \right)}\forall i=1,\dots ,c,\forall \alpha ,\beta ,\dots ,\zeta$$

2.

Chemical-reaction equilibrium

$$A_r={\displaystyle \sum _{i=1}^c}{\nu }_{i,r}{\mu }_i=0\forall r=1,\dots ,R_{\mathrm{ind}}$$

(55)

These conditions recover the full set of classical equilibrium relations, including multi-phase phase equilibria and the thermodynamic equilibrium constants for simultaneous reactions [18,19,20,21]:

$$K_r=\prod _i{a}_i^{{\nu }_{i,r}}=\exp \left(-{\displaystyle \frac{\mathsf{\Delta }{G}_r^{\circ }}{RT}}\right)$$

(56)

Definitions

• ${\mu }_i^{\left(\alpha \right)}$: chemical potential of component $i$ in phase $\alpha$.
• $A_r$: affinity of the $r$-th independent reaction.
• $K_r$: thermodynamic equilibrium constant of the $r$-th reaction.
• $a_i$: activity of component $i$ (incorporating both composition and the helical fugacity/activity coefficients from ZMT).

Proof of Theorem 9.

1.

Step 1 (Total phase functional) The total phase functional of the entire system is

$$\underset{‾}{G}={\displaystyle \sum _{\alpha =1}^{\zeta }}{\underset{‾}{G}}^{\left(\alpha \right)}+\mathrm{terms} \mathrm{involving} \left\{{\xi }_r\right\}$$

2.

Variational stationarity: Equilibrium requires $\underset{‾}{G}$ to be stationary with respect to every allowed variation:

• transfer of any component between any phases,
• advancement of any independent reaction extent ${\xi }_r$

3.

Phase equilibrium from matter transfer: A virtual transfer $\delta n_i$ of component $i$ from phase $\alpha$ to phase $\beta$ produces zero change in $\underset{‾}{G}$ only if

$${\mu }_i^{\left(\alpha \right)}={\mu }_i^{\left(\beta \right)}\mathrm{for} \mathrm{every} \mathrm{pair} \mathrm{of} \mathrm{phases}$$

Extending to all $\zeta$ phases yields the multi-phase chemical-potential equality.

• Chemical equilibrium from reaction advancement: A virtual advancement $\delta {\xi }_r$ of the $r$-th independent reaction produces

  $${\displaystyle \frac{\partial \underset{‾}{G}}{\partial {\xi }_r}}=\sum _i{\nu }_{i,r}{\mu }_i=A_r$$

  (57)
  Stationarity requires

  $$A_r=0\forall r$$

  (58)
• Recovery of classical constants
  Substituting the expression for ${\mu }_i$ from Theorem 8 into $A_r=0$ yields

  $$\sum _i{\nu }_{i,r}{\mu }_i=\mathsf{\Delta }{G}_r^{\circ }+RT\ln \left(\prod _i{a}_i^{{\nu }_{i,r}}\right)=0$$

  which rearranges directly to the classical equilibrium constant

  $$K_r=\prod _i{a}_i^{{\nu }_{i,r}}=\exp \left(-{\displaystyle \frac{\mathsf{\Delta }{G}_r^{\circ }}{RT}}\right)$$

Corollary 13. All thermodynamic equilibria — phase, chemical, and their combination in multi-phase reactive systems — are unified under the single principle of minimizing the total phase functional $\underset{‾}{G}$

Corollary 14. The number of independent variables is reduced by $c\left(\zeta -1\right)$ (from chemical-potential equalities) plus $R_{\mathrm{ind}}$ (from vanishing affinities), recovering the generalized Gibbs phase rule for reactive multi-phase systems.

7.9. Conjugate Coordinates $X\left(s;\xi \right)$ and $Y\left(s;\xi \right)$ in Reactive Systems Theorem

Theorem 10. In the Zeta-Minimizer framework, the phase functional of a reactive system defined on the kernel $Z\left(s;\xi \right)$ induces a pair of canonically conjugate coordinates $\left(X,Y\right)$ on the characteristic scaling surface. These coordinates satisfy the Euler homogeneity relation and provide the natural variables in which both the driving force and Hessian curvature are expressed as functions of the reaction extent $\xi$.

Statement:

On the coordinate chart $\left(s,X\right)$, the conjugate pair is given explicitly by

$$X\left(s;\xi \right)=-{\displaystyle \frac{{\left[Z\left(s;\xi \right)\right]}^2}{Z^{\prime }\left(s;\xi \right)}},Y\left(s;\xi \right)=A_F s{ Z}^{\prime }\left(s;\xi \right)$$

(59)

where $s=\ln 19/T$, $A_F$ is the scaling prefactor, and the prime denotes differentiation with respect to $s$ at fixed $\xi$. These coordinates obey the fundamental relation

$$XY=A_F s Z\left(s;\xi \right)$$

(60)

and the potential identification

$$G=-A_F s\ln Z\left(s;\xi \right)$$

Domain: $A_F\ne 0$, $s\ne 0$, $Z\left(s;\xi \right)>0$, $Z\left(s;\xi \right)\ne 1$, $Z^{\prime }\left(s;\xi \right)\ne 0$ (ensuring all logarithms and derivatives are well-defined).

Proof of Theorem 10.

• Conjugate relation and potential: The scaling surface satisfies the Euler homogeneity condition

  $$XY=A_F s Z\left(s;\xi \right)$$

  (60)

  and the phase functional is identified with the potential

  $$G=-A_F s\ln Z\left(s;\xi \right)$$

  (61)
  The conjugate momentum is defined by (Eq. (29))

  $$Y={\left({\displaystyle \frac{\partial G}{\partial X}}\right)}_s$$
• Total differential: In the $\left(s,X\right)$ chart,

  $$dG={\left({\displaystyle \frac{\partial G}{\partial s}}\right)}_{X}ds+YdX$$
  From Eq. (60),

  $$Y={\displaystyle \frac{A_F s Z\left(s;\xi \right)}{X}}$$
  Leading to,

  $$dG={\left({\displaystyle \frac{\partial G}{\partial s}}\right)}_{X}ds+{\displaystyle \frac{A_F s Z\left(s;\xi \right)}{X}}dX$$

  (62)
• Explicit differentiation of the potential: Differentiating Eq. (61) with respect to $s$ at fixed $X$ gives

  $$dG=\left[-A_F\ln Z\left(s;\xi \right)-A_F s{\displaystyle \frac{Z^{\prime }\left(s;\xi \right)}{Z\left(s;\xi \right)}}\right]ds$$

  (63)

1. Canonical identification along the scaling path: Equating (62) and (63) along the characteristic path of the reactive system and imposing the normalization that the infinitesimal displacement satisfies $dX=d\xi$ (i.e., the extent of reaction is the natural coordinate along the surface), the logarithmic terms cancel. Dividing the remaining equation by the common non-zero factor $A_F s$ yields

$${\displaystyle \frac{Z\left(s;\xi \right)}{X}}=-{\displaystyle \frac{Z^{\prime }\left(s;\xi \right)}{Z\left(s;\xi \right)}}$$

2. Solution for the conjugates: Solving for $X$ immediately gives

$$X\left(s;\xi \right)=-{\displaystyle \frac{{\left[Z\left(s;\xi \right)\right]}^2}{Z^{\prime }\left(s;\xi \right)}}$$

Substitution into (60) then produces the conjugate momentum

$$Y\left(s;\xi \right)=As{Z}^{\prime }\left(s;\xi \right)$$

Corollary 15. The first derivative of the phase functional (driving force) with respect to the extent of reaction takes the compact form

$${\displaystyle \frac{\partial G}{\partial \xi }}=Y\left(s;\xi \right){\displaystyle \frac{\partial X}{\partial \xi }}{\mid }_s+\mathrm{interface} \mathrm{terms}$$

which isolates the conjugate leverage of pressure and inerts.

Corollary 16. The Hessian (second derivative) is expressible directly in terms of derivatives of $X\left(s;\xi \right)$ and $Y\left(s;\xi \right)$, providing a natural coordinate system in which both relaxation time and stability are analyzed.

Corollary 17. The construction holds for both single-reaction ($\xi$) and multi-reaction ($\left\{{\xi }_r\right\}$) cases and bridges seamlessly to the non-reactive limit ($\xi =0$).

7.10. Vapor–Solid Equilibrium with Multiple Reactions: Interface Matching and Gibbs Free Energy of Reaction Theorem

Theorem 11. In the Zeta-Minimizer framework, for a multi-component, multi-reaction system in vapor–solid equilibrium, thermodynamic equilibrium requires simultaneous satisfaction of (i) chemical-potential equality between vapor and solid phases for every species and (ii) vanishing affinity for every independent reaction. The interface matching constants $C_k$ from grand-potential continuity enter directly into the effective standard Gibbs free energy of reaction $\mathsf{\Delta }{G}_r^{\circ }$, thereby shifting the equilibrium constant in heterogeneous catalysis.

Statement:

At vapor–solid equilibrium Eq. (44) and Eq. (55):

$${\mu }_i^{(\mathrm{V})}={\mu }_i^{(\mathrm{S})}\forall i=1,\dots ,c,$$

(64)

$$A_r=\sum _i{\nu }_{i,r}{\mu }_i=0\forall r=1,\dots ,R_{\mathrm{ind}}$$

The interface matching condition (Eq. (24))

$${\omega }_{\mathrm{solid}}\left(s\right)={\omega }_{\mathrm{fluid}}\left(s\right)+\sum _k{C}_k$$

modifies the chemical potentials such that the effective reaction Gibbs energy becomes

$$\mathsf{\Delta }{G}_r^{\circ }\left(\mathrm{eff}\right)=\mathsf{\Delta }{G}_r^{\circ }+\sum _i{\nu }_{i,r}\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}$$

(65)

where $\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}$ is the contribution of the interface constants to species $i$. Consequently, the equilibrium constant is

$$K_r=\prod _i{a}_i^{{\nu }_{i,r}}=\exp \left(-{\displaystyle \frac{\mathsf{\Delta }{G}_r^{\circ }\left(\mathrm{eff}\right)}{RT}}\right)$$

(66)

Definitions

• ${\mu }_i^{(\mathrm{V})},{\mu }_i^{(\mathrm{S})}$: chemical potential of component $i$in vapor and solid phases.
• $A_r$: affinity of the $r$-th independent reaction.
• $C_k$: interface matching constants fixed by critical gear compositions (Section 4).
• $\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}$: species-specific shift induced by the grand-potential jump $\sum C_k$.
• $\mathsf{\Delta }{G}_r^{\circ }\left(\mathrm{eff}\right)$: effective standard Gibbs free energy of reaction incorporating interface effects.

Proof of Theorem 11.

1. Phase-functional description of both phases: By Theorems 7 and 8, each phase possesses its own extensive phase functional ${\underset{‾}{G}}^{\left(\alpha \right)}$and chemical potentials (Eq. (40))

$${\mu }_i^{\left(\alpha \right)}=G^{\left(\alpha \right)}+RT\ln \left(x_i^{\left(\alpha \right)}{Z}_i^{\left(\alpha \right)}\left(s\right)\right)+\mathsf{\Delta }{\mu }_i^{\left(C_k,\alpha \right)}$$

(67)

Interface matching from Section 4: Grand-potential continuity across the vapor–solid boundary (Section 4.1) requires (Eq. (24))

$${\omega }_{\mathrm{solid}}\left(s\right)={\omega }_{\mathrm{fluid}}\left(s\right)+\sum _k{C}_k$$

which translates directly into a difference in the interface corrections (Eq. (67)):

$$\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}={\mu }_i^{(\mathrm{S})}-{\mu }_i^{(\mathrm{V})}{\mid }_{\mathrm{no} \mathrm{interface}}=\sum _k{\nu }_{i,k}{C}_k$$

3. Chemical-potential equality: Vapor–solid equilibrium (Theorem 9) (Eq. (64)) demands

$${\mu }_i^{(\mathrm{V})}={\mu }_i^{(\mathrm{S})}\forall i$$

Substituting the expression from Step 1 and the interface shift from Step 2 yields the constraint that couples composition, helical factors $Z_i\left(s\right)$, and the constants $C_k$.

4. Reaction affinity and effective $\mathsf{\Delta }{G}_r^{\circ }\left(\mathrm{eff}\right)$): The affinity of the $r$-th reaction (Eq. (58)) is

$$A_r=\sum _i{\nu }_{i,r}{\mu }_i$$

Inserting the interface-modified chemical potentials gives

$$A_r=\mathsf{\Delta }{G}_r^{\circ }+RT\ln \left(\prod _i{a}_i^{{\nu }_{i,r}}\right)+\sum _i{\nu }_{i,r}\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}=0$$

Defining the effective standard Gibbs energy

$$\mathsf{\Delta }{G}_r^{\circ }\left(\mathrm{eff}\right)=\mathsf{\Delta }{G}_r^{\circ }+\sum _i{\nu }_{i,r}\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}$$

immediately produces the shifted equilibrium constant

$$K_r=\prod _i{a}_i^{{\nu }_{i,r}}=\exp \left(-{\displaystyle \frac{\mathsf{\Delta }{G}_r^{\circ }\left(\mathrm{eff}\right)}{RT}}\right)$$

5. Conjugate-coordinate confirmation: By Theorem 10 the driving force $\partial G/\partial {\xi }_r$ is expressible in the conjugate coordinates $\left(X,Y\right)$. The interface constants $C_k$ appear in the source term of the fugacity Hessian, ensuring that both phase matching and reaction stationarity are satisfied simultaneously at the global minimum of the total phase functional $\underset{‾}{G}$

7.10.1. Chemical Potentials with Explicit Interface Shift

From Theorems 8 and 11 the chemical potential of every gas-phase species is

$${\mu }_i={\displaystyle \frac{\partial \underset{‾}{G}}{\partial n_i}}{\mid }_{\left\{{\xi }_r\right\}}=G+RT\mathrm{l}\mathrm{n}\left(x_i{Z}_i\left(s\right)\right)+\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}$$

where the interface contribution is

$$\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}=\sum _k{\nu }_{i,k}{C}_k(\mathrm{fixed} \mathrm{scalar} \mathrm{determined} \mathrm{by} \mathrm{the} \mathrm{catalyst} \mathrm{surface})$$

7.10.2. Affinity Equations in Full ZMT Form

The stationarity condition for each independent gas-phase reaction extent ${\xi }_r$ (Theorem 9) is the vanishing of the phase-functional derivative:

$$A_r={\displaystyle \frac{\partial \underset{‾}{G}}{\partial {\xi }_r}}=0,r=1,\dots ,R_{\mathrm{ind}}$$

Substituting the chemical-potential expression and collecting terms yields the exact system

$${\displaystyle \frac{\partial \underset{‾}{G}}{\partial {\xi }_r}}=\int (\partial \underset{‾}{G}/\partial {\xi }_r{)}_{\left\{\xi =0\right\}}d{\xi }_r+RT\ln \left(\prod _i[\left.x_i\left(\mathit{\xi }\right)Z_i\left(s\right){]}^{{\nu }_{i,r}}\right)\right.+{\delta }_r^{\left(C_k\right)}=0$$

(68)

with the purely ZMT interface shift

$${\delta }_r^{\left(C_k\right)}=\sum _i{\nu }_{i,r}\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}=\sum _i\sum _k{\nu }_{i,r}{\nu }_{i,k}{C}_k$$

(69)

7.10.3. Complete Replacement of Classical Thermodynamics

All classical quantities [18,19,20,21] are now fully replaced by their deductive ZMT counterparts:

• Classical ${\mu }_i$→ $\partial \underset{‾}{G}/\partial n_i$ along the stoichiometric subspace
• Classical activity $a_i$ → $x_i{Z}_i\left(s\right)$(helical fugacity factor)
• Classical $\mathsf{\Delta }{G}_r^{\circ }$→ $=\int (\partial \underset{‾}{G}/\partial {\xi }_r{)}_{\left\{\xi =0\right\}}d{\xi }_r$
• Classical equilibrium constant $K_r$→ solution of $\partial \underset{‾}{G}/\partial {\xi }_r=0$
• Classical reference state shift → interface term ${\delta }_r^{\left(C_k\right)}$arising from grand-potential continuity ${\omega }_{\mathrm{solid}}={\omega }_{\mathrm{fluid}}+{\sum }_k{C}_k$

The final fully deductive ZMT system that must be solved for the extent vector $\xi$ is therefore

$$\boxed{{\displaystyle \frac{\partial \underset{‾}{G}}{\partial {\xi }_r}}(\mathit{\xi };\left\{C_k\right\},s)=0,r=1,\dots ,R_{\mathrm{ind}}} .$$

(70)

where $\underset{‾}{G}=\underset{‾}{G}(n_{\mathrm{total}},x_i\left(\xi \right),Z_i\left(s\right),C_k)$ is the total phase functional, and the only external parameters are the fixed catalyst constants $\left\{C_k\right\}$ and the temperature coordinate $s$

7.10.4. Physical Interpretation

All gas-phase reactions, regardless of how many and how strongly coupled, are anchored to the solid surface through the single, shared set $\left\{C_k\right\}$. The global minimization of $\underset{‾}{G}$ with respect to the vector $\xi$ simultaneously satisfies phase matching, reaction equilibrium, and interface continuity. No classical thermodynamic postulates remain; every term is a direct derivative or evaluation of the phase functional on the helical manifold.

7.11. Illustrative Case: Single Reaction with Inerts in Vapor–Solid Equilibrium

Consider the practically most common situation: a single gas-phase reaction of arbitrary stoichiometry in contact with a fixed solid catalyst, in the presence of non-reactive inerts.

1. Reaction (Eq. (45))

$$\sum _i{\nu }_i{\mathrm{React}}_i\rightleftharpoons \sum _j{\nu }_j{\mathrm{Prod}}_j,{\nu }_i<0 (\mathrm{reactants}),{\nu }_j>0 (\mathrm{products})$$

2. Variables

• Single extent of reaction: $\xi$
• Inert-to-nitrogen feed ratio: $\gamma =n_{\mathrm{inert}}/n_{{\mathrm{N}}_2}^0$ (fixed)
• Total moles: $N\left(\xi \right)=N_0\left(4-2\xi +\gamma \right)$ (example for ammonia-like stoichiometry)
• Mole fractions: $x_i\left(\xi \right)=x_i\left(0\right)+{\displaystyle \frac{{\nu }_k\xi }{N\left(\xi \right)}}$ for reacting species, $x_{\mathrm{inert}}=\gamma N_0/N\left(\xi \right)$

3. Equilibrium Condition (Single Equation)

The system reaches equilibrium when the total phase functional is stationary with respect to $\xi$:

$${\displaystyle \frac{\partial \underset{‾}{G}}{\partial \xi }}=0$$

Explicitly, after substituting Theorems 7–11 and Section 4:

${\displaystyle \frac{\partial \underset{‾}{G}}{\partial \xi }}=\underset{\mathrm{helical} \mathrm{species} \mathrm{term}}{\underbrace{\int (\partial \underset{‾}{G}/\partial \xi {)}_{\xi =0}d\xi \left(\mathrm{gas}\right)}}+RT\ln \left(\prod _k[\left.x_i\left(\mathit{\xi }\right)Z_i\left(s\right){]}^{{\nu }_i}\right)\right.+\underset{\mathrm{interface} \mathrm{shift}}{\underbrace{{\delta }^{\left(C_k\right)}}}+\underset{\mathrm{from} N\left(\xi \right)}{\underbrace{\mathrm{normalization} \mathrm{term} \mathrm{from} \mathrm{inerts}}}=0$

where:

• $\int (\partial \underset{‾}{G}/\partial \xi {)}_{\xi =0}d\xi$ (purely from the phase functional)
• $Z_k\left(s\right)$= helical gear contribution of species $k$
• Interface shift (fixed by catalyst surface):

  $${\delta }^{\left(C_k\right)}=\sum _i\sum _k{\nu }_i{\nu }_{i,k}{C}_k(\sin \mathrm{gle} \mathrm{constant})$$
• Normalization term from inerts appears naturally as the derivative of the $-4/N\left(\xi \right)$ term in the driving force, i.e., dilution penalty is explicitly $\propto 1/\left(4-2\xi +\gamma \right)$.

4. Final Equilibrium Equation (Fully Deductive ZMT Form)

$$RT\ln \left({\displaystyle \frac{\prod _{\mathrm{prod}}[x_j\left(\xi \right)Z_j\left(s\right){]}^{{\nu }_j}}{\prod _{\mathrm{react}}[x_i\left(\xi \right)Z_i\left(s\right){]}^{\mid {\nu }_i\mid }}}\right)+{\delta }^{\left(C_k\right)}+{\displaystyle \frac{\partial }{\partial \xi }}\left(-RT\ln N\left(\xi \right)\right)=0$$

(71)

or, compactly,

$$\boxed{{\displaystyle \frac{\partial \underset{‾}{G}}{\partial \xi }}(s,\xi ,\gamma ,\left\{C_k\right\})=0}$$

(72)

Interpretation

• The entire equilibrium state (final composition, conversion, effect of pressure, effect of inerts) is determined by solving one single scalar equation in $\xi$.
• All physical effects appear naturally inside the phase functional derivative:
  • Helical species contributions → $Z_k\left(s\right)$
  • Inert dilution → explicit $N\left(\xi \right)$ dependence
  • Catalyst surface effect → single fixed number ${\delta }^{\left(C_k\right)}$
  • Pressure leverage → through the conjugate coordinate $X\left(s;\xi \right)$ hidden inside $\partial \underset{‾}{G}/\partial \xi$

This is the complete ZMT replacement of the classical $K=\prod a_i^{{\nu }_i}$. There are no leftover classical terms [18,19,20,21] — everything is a derivative or evaluation of the phase functional $\underset{‾}{G}$ evaluated at the fixed catalyst fingerprint $\left\{C_k\right\}$ and inert level $\gamma$.

8. Variational Reaction Rate Law and Emergence of Batch and Residence Time

The phase functional $\underset{‾}{G}$ (identified with the Gibbs free energy through the grand-partition function $Z\left(s;\xi \right)$) is the single variational generator of all dynamical processes in the Zeta-Minimizer framework. Projecting the Euler–Lagrange equations of the associated action onto the multi-dimensional reaction coordinate spanned by the extent vector $\xi$ yields a unified variational rate law that governs chemical kinetics and directly produces batch time and residence time as accumulated relaxation along the helical manifold.

8.1. Variational Reaction Rate Law

The time-dependent action on the helical manifold is

$$S={\int }_{-\infty }^{\infty }F\left[A\left(\theta \right),\psi \left(\theta \right)\right]d\theta .$$

(73)

where emergent time $\theta$ (or $t$) is the helical coordinate. Stationarity with respect to variations in $\xi$ produces the exact nonlinear variational rate law

$${\displaystyle \frac{d\mathit{\xi }}{dt}}={\mathit{\lambda }}_{\mathbf{m}\mathbf{i}\mathbf{n}}·{\displaystyle \frac{1}{{\mathbf{H}}_{\mathrm{eq}}}}{\left({\displaystyle \frac{\partial \ln Z\left(s;\xi \right)}{\partial \xi }}\right)}_{s,X,\gamma ,\left\{C_k\right\}}$$

(74)

where ${\mathit{\lambda }}_{\mathbf{m}\mathbf{i}\mathbf{n}}$ is the vector of characteristic relaxation rates and ${\mathbf{H}}_{\mathrm{eq}}$ is the positive definite reaction Hessian evaluated at equilibrium.

The instantaneous thermodynamic driving force $\partial \ln Z/\partial \xi$ contains the helical species contributions, the normalization term from total moles $N\left(\mathit{\xi },\gamma \right)$, and the conjugate response through the coordinate $X\left(s;\mathit{\xi }\right)$. Inerts enter explicitly via dilution in $N\left(\mathit{\xi },\gamma \right)$Near equilibrium the system linearizes to

$${\displaystyle \frac{d\mathit{\xi }}{dt}}=-{\mathit{\lambda }}_{\mathbf{m}\mathbf{i}\mathbf{n}}·\left(\xi -{\xi }_{\mathrm{eq}}\right)$$

(75)

8.2. Hessian Construction

The reaction Hessian matrix is defined as the second derivative of the total phase functional with respect to the extent vector:

$$H_{rs}={\displaystyle \frac{{\partial }^2\underset{‾}{G}}{\partial {\xi }_r\partial {\xi }_s}}=-RT{\displaystyle \frac{{\partial }^2\ln Z\left(s;\mathit{\xi }\right)}{\partial {\xi }_r\partial {\xi }_s}}{\mid }_{s,X,\gamma ,\left\{C_k\right\}}$$

(76)

evaluated at equilibrium.

The explicit form contains three deductively distinct contributions:

• Helical species terms from $\ln Z_i\left(s\right)$
• Normalization/extensivity terms arising from the total-moles dependence $1/N\left(\mathit{\xi },\gamma \right)$, whose leading coefficient is determined by the stoichiometric mole change $\mathsf{\Delta }n$ and the chosen feed normalization (e.g., the coefficient is 4 for ammonia synthesis with standard 1 N₂ + 3 H₂ feed, and 0 for reactions with $\mathsf{\Delta }n=0$ such as water-gas shift),
• Conjugate-coordinate terms involving derivatives of $X\left(s;\xi \right)$The eigenvalues ${\lambda }_r>0$ of this Hessian (guaranteed positive by Axiom 1 — strict concavity) are the physical relaxation rates of the coupled reaction network. The smallest eigenvalue ${\lambda }_{\mathrm{m}\mathrm{i}\mathrm{n}}$ governs the slowest relaxation mode and enters the variational rate law above.

8.3. Emergence of Batch and Residence Time

Batch time $t$ (or residence time $\tau =V_{\mathrm{reactor}}/{\dot{V}}_{\mathrm{flow}}$) is the control parameter that determines how closely the system approaches the equilibrium manifold. Define the accumulated dimensionless residence time

$${\tau }_{\mathrm{acc}}\left(t\right)={\int }_0^t{\mathit{\lambda }}_{\mathbf{m}\mathbf{i}\mathbf{n}}\left(t^{\prime }\right){\mid }_{s,X,\gamma ,\left\{C_k\right\}}d{t}^{\prime }$$

(77)

The extent vector then evolves as

$$\mathit{\xi }\left(t\right)={\mathit{\xi }}_{\mathrm{eq}}\left[1-\mathrm{e}\mathrm{x}\mathrm{p}\left(-{\tau }_{\mathrm{acc}}\right)\right]$$

(78)

Short ${\tau }_{\mathrm{acc}}$ corresponds to kinetic control (system far from equilibrium). Long ${\tau }_{\mathrm{acc}}$ drives the outlet composition onto the thermodynamic equilibrium manifold defined by grand-potential continuity ${\omega }_{\mathrm{solid}}={\omega }_{\mathrm{fluid}}+{\sum }_k{C}_k$. All quantities are fully deductive objects of the phase functional $\underset{‾}{G}$, with no residual classical thermodynamic postulates.

8.4. Adiabatic Temperature Profile from the Phase Functional

In the Zeta-Minimizer framework, no separate energy balance is required. For an adiabatic multi-component, multi-reaction system with inerts, the total phase functional $\underset{‾}{G}$ itself provides the direct relationship between the extent vector $\mathit{\xi }$ and temperature $T$ (via $s=\ln 19/T$) along any adiabatic trajectory.

Because the system is closed and adiabatic, the variational evolution must satisfy the constraint

$$d\underset{‾}{G}=0$$

along the path where both $s$ and $\xi$ vary while pressure (conjugate coordinate $X$) and inert level $\gamma$ are held constant. Starting from (Eq. (52)),

$$\underset{‾}{G}=n_{\mathrm{total}}·\left(-RT\ln Z\left(s;\mathit{\xi };\gamma \right)\right)$$

and taking the total differential, the adiabatic temperature profile is obtained as

$${\left({\displaystyle \frac{dT}{d\mathit{\xi }}}\right)}_{\mathrm{adiabatic}}=\left[{\displaystyle \frac{T·{\left(\partial \ln Z/\partial \mathit{\xi }\right)}_{s,X,\gamma ,\left\{C_k\right\}}}{\ln Z\left(s;\mathit{\xi };\gamma \right)-T{\left(\partial \ln Z/\partial T\right)}_{X,\mathit{\xi },\gamma ,\left\{C_k\right\}}}}\right]$$

(79)

All quantities on the right-hand side are known deductive functions of the helical partition factors $Z_i\left(s\right)$, the current composition $x_i\left(\mathit{\xi },\gamma \right)$, the conjugate coordinate $X\left(s;\xi \right)$, and the fixed catalyst interface constants $\left\{C_k\right\}$. This relation holds for arbitrary stoichiometry, any number of coupled reactions, and arbitrary inert content without invoking a separate enthalpy balance.

8.5. Deductive Linear Scaling from Closed Prime Functions Theorem

Theorem 12. In the Zeta-Minimizer framework, linear scaling relations between adsorption energies (or free energies) of different reaction intermediates on a catalyst surface emerge deductively as a direct consequence of the Integer Gear-Up-to-Prime Rule and the closed analytic form of the interface constants $C_k$.

Let a solid catalyst activate a discrete set of period primes $\left\{p_k\right\}$ according to the Integer Gear-Up-to-Prime Rule. Then the adsorption free energy (or effective binding strength) of any intermediate $i$ corresponding to gear $p_k$ satisfies a linear scaling relation of the form

$$\mathsf{\Delta }{G}_i={\alpha }_{ij}\mathsf{\Delta }{G}_j+{\beta }_{ij}\left(\left\{C_k\right\}\right)$$

(80)

where the slope ${\alpha }_{ij}$ and intercept ${\beta }_{ij}$ are closed functions of the primes and the fixed interface constants $\left\{C_k\right\}$. No empirical fitting parameters are required.

Definitions

• $\mathsf{\Delta }{G}_i$: effective adsorption free energy of intermediate $i$(identified with $\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}+RT\ln Z_i\left(s\right)$ term).
• $\left\{C_k\right\}$: interface matching constants fixed by the catalyst surface (Section 4).
• ${\alpha }_{ij},{\beta }_{ij}$: universal scaling coefficients determined solely by the prime structure.

Proof of Theorem 12.

1. By the Integer Gear-Up-to-Prime Rule, each adsorbed intermediate $i$ is assigned a unique period prime $p_m$. Its interaction with the master hub (prime 19) is governed by the interaction parameter whose critical-point solution yields the closed form (Eq. (25))

$$C_k=-{\displaystyle \frac{\left(k-19\right)x_{k,o}^2}{1-2x_{k,o}}}$$

where $x_{k,o}$ is the explicit critical composition (Eq. (14)).

2. The effective chemical potential shift of the adsorbed state is

$$\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}\propto C_k+\mathrm{terms} \mathrm{in} \ln Z_i\left(s\right)$$

Because $C_k$ is a rational function of the prime $p_m$, and different intermediates correspond to different primes $p_n$ and $p_m$, the ratio of shifts between any two intermediates is

$${\displaystyle \frac{\mathsf{\Delta }{\mu }_i^{\left(C_n\right)}}{\mathsf{\Delta }{\mu }_j^{\left(C_m\right)}}}=f\left(p_n,p_m\right)$$

(81)

which is a closed algebraic function of the primes.

3. In the low-coverage or linear-response regime (small deviations from critical composition), the function $f\left(p_n,p_m\right)$ reduces to a linear relation

$$\mathsf{\Delta }{G}_i={\alpha }_{ij} \left(p_k,p_m\right)\mathsf{\Delta }{G}_j+{\beta }_{ij}\left(\left\{C_k\right\}\right)$$

(82)

where the slope ${\alpha }_{ij}$ depends only on the ratio of the activated primes, and the intercept ${\beta }_{ij}$ is fixed by the global set $\left\{C_k\right\}$ of the catalyst surface.

4. Because the same fixed set $\left\{C_k\right\}$ controls all intermediates simultaneously through the grand-potential continuity condition, all scaling lines for a given catalyst are rigidly linked. Changing the catalyst corresponds to a global shift of the entire family of linear relations through the common vector $\left\{C_k\right\}$.

5. The Lyapunov exponents ${\lambda }_k\left(\left\{p_k\right\}\right)$ (encoding stability of the coupling) further constrain the slopes ${\alpha }_{ij}$ to lie on universal lines determined by the prime basis, completing the deductive origin of linear scaling without external input. □

8.6. Emergence of Volcano Plots from Closed Prime Functions Theorem

Theorem 13. Let $d$ be any monotonic descriptor of catalyst strength constructed from the net grand-potential shift generated by an active set of period primes (e.g., a linear combination of the interface constants $\left\{C_k\right\}$). Then the activity realized at finite accumulated residence time ${\tau }_{\mathrm{acc}}$ (Eq. (77)) is given exactly by

$$A\left(d\right)={\xi }_{\mathrm{eq}}\left(d\right)\left[1-\mathrm{e}\mathrm{x}\mathrm{p}\left(-{\lambda }_{\min }\left(d\right)·{\tau }_{\mathrm{acc}}\right)\right].$$

(83)

where ${\xi }_{\mathrm{eq}}\left(d\right)$ is the equilibrium conversion and ${\lambda }_{\min }\left(d\right)$ is the smallest positive eigenvalue of the reaction Hessian, both evaluated at the composition corresponding to descriptor $d$This function necessarily exhibits a single maximum (volcano peak).

Proof of Theorem 13.

1. As $d$ increases, the thermodynamic driving force increases monotonically, so ${\xi }_{\mathrm{eq}}\left(d\right)$ is a strictly increasing function (Theorem 11).

2. As $d$ increases, the interface shift ${\delta }^{\left(C_k\right)}\left(d\right)$ moves the system. Initially, ${\lambda }_{\min }\left(d\right)$ increases (stronger binding improves kinetics), but beyond a critical value the marginal-stability constraint ${\lambda }_k\le 0$ for the lowest gear forces ${\lambda }_{\min }\left(d\right)$ to decrease rapidly (from the explicit form of the Lyapunov exponent in Section 4.3).

3. The product $A\left(d\right)={\xi }_{\mathrm{eq}}\left(d\right)·f\left({\lambda }_{\min }\left(d\right)\right)$, where $f\left(\lambda \right)=1-e^{-\lambda {\tau }_{\mathrm{acc}}}$ is a strictly increasing but concave function of $\lambda$, has a derivative that changes sign exactly once: it is positive when the gain in thermodynamics and initial kinetic improvement dominates, and negative when the kinetic penalty from excessive binding (over-stabilization) dominates. Hence $A\left(d\right)$ possesses exactly one maximum.

4. Because both ${\xi }_{\mathrm{eq}}\left(d\right)$ and ${\lambda }_{\min }\left(d\right)$ are closed functions of the prime set and the interface constants $\left\{C_k\right\}$, the volcano shape is a theorem, not an empirical correlation. □

8.7. Variational Dynamics of Vapor–Solid Reactive Equilibrium Theorem

Theorem 14. The static vapor–solid equilibrium condition is the fixed point of a fully dynamical process governed by the variational rate law. For a general multi-component, multi-reaction gas-phase system with inerts in contact with a solid catalyst, let $\xi$ be the vector of gas-phase reaction extents and $\left\{{\xi }_i^{\mathrm{ads}}\right\}$ the vector of adsorption extents measuring mass transfer of each component between vapor and solid phases.

The phase functional $\underset{‾}{G}$ generates the exact nonlinear variational rate laws

$${\displaystyle \frac{d\mathit{\xi }}{dt}}={\mathit{\lambda }}_{\mathrm{gas}}·{\left({\displaystyle \frac{\partial \ln Z}{\partial \mathit{\xi }}}\right)}_{s,X,\gamma ,\left\{C_k\right\}},{\displaystyle \frac{d{\xi }_i^{\mathrm{ads}}}{dt}}={\lambda }_i^{\mathrm{ads}}{\left({\displaystyle \frac{\partial \ln Z}{\partial {\xi }_i^{\mathrm{ads}}}}\right)}_{s,X,\gamma ,\left\{C_k\right\}}$$

(84)

where all derivatives are evaluated subject to the interface matching condition ${\omega }_{\mathrm{solid}}={\omega }_{\mathrm{fluid}}+{\sum }_k{C}_k$ (Section 4 (Eq. (24)) and Theorem 11).

Near equilibrium these linearize to

$${\displaystyle \frac{d\mathit{\xi }}{dt}}=-{\mathit{\lambda }}_{\mathrm{gas}}·\left(\mathit{\xi }-{\mathit{\xi }}_{\mathrm{eq}}\right),{\displaystyle \frac{d{\xi }_i^{\mathrm{ads}}}{dt}}=-{\lambda }_i^{\mathrm{ads}}\left({\xi }_i^{\mathrm{ads}}-{\xi }_{i,\mathrm{eq}}^{\mathrm{ads}}\right)$$

(85)

with all relaxation rates $\lambda$ being the positive eigenvalues of the combined gas + solid reaction Hessian.

Proof of Theorem 14.

• By the Variational Reaction Rate Law (Section 9.1), the phase functional $\underset{‾}{G}$ (Eq. (52)) generates the exact rate law for any set of variational coordinates. Choosing the full set $\left(\mathit{\xi },\left\{{\xi }_i^{\mathrm{ads}}\right\}\right)$ as coordinates yields the nonlinear system above.
• The interface matching condition ((Eq. (24)) and Theorem 11) is enforced at every instant, so all partial derivatives are taken on the constrained manifold where ${\omega }_{\mathrm{solid}}-{\omega }_{\mathrm{fluid}}={\sum }_k{C}_k$
• Linearization around the static equilibrium point $\left({\mathit{\xi }}_{\mathrm{eq}},\left\{{\xi }_{i,\mathrm{eq}}^{\mathrm{ads}}\right\}\right)$ (defined by Theorems 9 and 11) produces the matrix relaxation equation. The eigenvalues of the combined Hessian are strictly positive by Axiom 1 (strict concavity of $\underset{‾}{G}$).
• All quantities — helical factors $Z_i\left(s\right)$, macroscopic mole fractions $x_i\left(\mathit{\xi },\gamma \right)$, interface constants $\left\{C_k\right\}$, and the resulting Hessian — are closed deductive functions of the active period primes and the Integer Gear-Up-to-Prime Rule (Theorems 12 and 13). No free parameters remain.

8.8. Modal Decomposition of Coupled Vapor–Solid and Gas-Phase Dynamics

The linearized dynamics of the full coupled system (gas-phase reaction extents $\mathit{\xi }$ and adsorption extents $\left\{{\xi }_i^{\mathrm{ads}}\right\}$) admit the matrix form

$${\displaystyle \frac{d}{dt}}\left(\begin{array}{c}\mathit{\xi }\\ {\mathit{\xi }}^{\mathrm{ads}}\end{array}\right)=-\mathsf{\Lambda }\left(\begin{array}{c}\mathit{\xi }-{\mathit{\xi }}_{\mathrm{eq}}\\ {\mathit{\xi }}^{\mathrm{ads}}-{\mathit{\xi }}_{\mathrm{eq}}^{\mathrm{ads}}\end{array}\right)$$

(86)

where $\mathsf{\Lambda }$ is the positive-definite Hessian matrix of the total phase functional $\underset{‾}{G}$ (or equivalently of $\ln Z\left(s\right)$) evaluated with respect to the complete set of coordinates $\left\{\mathit{\xi },{\mathit{\xi }}^{\mathrm{ads}}\right\}$ under the interface matching condition ${\omega }_{\mathrm{solid}}={\omega }_{\mathrm{fluid}}+{\sum }_k{C}_k$ ((Eq. (24)) and Theorem 11).

Diagonalization $\mathsf{\Lambda }=\mathbf{Q}{\mathsf{\Lambda }}_{\mathrm{diag}}{\mathbf{Q}}^{\mathbf{T}}$ yields the complete spectrum of positive relaxation rates $\left\{{\lambda }_k>0\right\}$ and the associated normal modes (eigenvectors ${\mathbf{v}}_{\mathit{k}}$).

Each pair $\left({\lambda }_k,{\mathbf{v}}_k\right)$ corresponds to an independent normal mode of the coupled dynamics:

• ${\lambda }_k$: relaxation rate of that mode,
• ${\mathbf{v}}_k$: the specific linear combination of gas-phase reaction progress and surface-coverage adjustments that relaxes cooperatively.

The smallest eigenvalue ${\lambda }_{\min }$ and its eigenvector govern the overall rate-limiting relaxation of the entire system toward vapor–solid equilibrium. Larger eigenvalues describe faster subprocesses (rapid equilibration of a particular adsorbed species, fast adjustment of individual gear occupations on the solid, etc.).

All eigenvalues and eigenvectors are closed deductive functions of the active period primes, the interface constants $\left\{C_k\right\}$, and the inert level $\gamma$ (Theorems 4, 7–11, and Section 8.2).

The single variational landscape of $\underset{‾}{G}$ (Eq. (52)) generates not only the static equilibrium point (Theorems 9–11) and the nonlinear rate law (Section 8.1), but also the full spectrum of time scales and physically interpretable normal modes of the coupled vapor–solid reactive system.

8.9. Emergence of the Classical Seven-Step Rate Expression from Variational Modal Decomposition Theorem

Theorem 15. The modal decomposition of the coupled dynamics (Theorem 14) naturally recovers the classical seven-step sequence of heterogeneous catalysis as dominant normal modes of the single variational landscape $\underset{‾}{G}$.

Statement:

When the full set of variational coordinates includes the gas-phase reaction extent vector $\mathit{\xi }$, the adsorption extents $\left\{{\xi }_i^{\mathrm{ads}}\right\}$, and the transport extents corresponding to external and internal diffusion, the eigenvectors ${\mathbf{v}}_k$ of the total Hessian $\mathsf{\Lambda }$ identify linear combinations of the following seven classical steps:

• External mass transfer of reactants to the particle surface
• Internal pore diffusion of reactants into the pellet
• Adsorption of reactants on active sites
• Surface chemical reaction
• Desorption of products
• Internal pore diffusion of products out of the pellet
• External mass transfer of products away from the surface

Proof of Theorem 15.

The observed macroscopic rate of reaction per unit reactor volume is

$$r_{\mathrm{obs}}=C_{\mathrm{cat}}·{\lambda }_{\min }\left(d\right)·\left({\xi }_{\mathrm{eq}}\left(d\right)-\xi \right).$$

(87)

where ${\lambda }_{\min }\left(d\right)$ is the smallest eigenvalue of the full Hessian, ${\xi }_{\mathrm{eq}}\left(d\right)$ is the equilibrium extent, $d$ is any monotonic descriptor of catalyst strength (e.g., a linear combination of $\left\{C_k\right\}$), and $C_{\mathrm{cat}}$ is the catalyst loading factor (mass of solid per reactor volume).

In the lumped multi-resistance approximation this becomes

$$r={\left(\sum _j{\displaystyle \frac{1}{k_j}}\right)}^{-1}·\mathsf{\Delta }(\mathrm{driving} \mathrm{force})$$

(88)

where each effective rate coefficient $k_j\propto {\lambda }_j$ (weighted by the components of eigenvector ${\mathbf{v}}_j$) corresponds to one of the seven steps, and the driving force is the full ZMT expression (Theorems 8–11).

All quantities are closed deductive functions of the active period primes, the interface constants $\left\{C_k\right\}$, the inert level $\gamma$, and the conjugate coordinate $X\left(s;\mathit{\xi }\right)$. No additional parameters are required. □

The classical seven-step sequence and the observed rate law are emergent normal modes and eigenvalues of the single variational phase functional $\underset{‾}{G}$

8.10. Catalyst Preparation as Variational Selection of the Active Prime Set and Interface Constants Theorem

Theorem 16. Catalyst preparation methods, reduction conditions, promoters, and support interactions act as external variational controls that select the active subset of period primes and the equilibrium occupation distribution on the solid surface. This selection directly and deductively determines the interface matching constants $\left\{C_k\right\}$, which in turn control all observable catalytic behavior.

Any change in catalyst preparation (precursor chemistry, pH, calcination temperature, reduction temperature, promoter loading, support material, etc.) modifies the relative stability of the available gears. This selects a specific subset of period primes $\left\{p_k\right\}$ and shifts the critical compositions $\left\{x_{k,o}\right\}$. The interface constants are then uniquely determined by the closed expression (Eq. (25))

$$C_k=-{\displaystyle \frac{\left(k-19\right)x_{k,o}^2}{1-2x_{k,o}}}$$

where $x_{k,o}$ (Eq. (14)) is the occupation that minimizes the mixture frequency at the solid–fluid interface in the $s\to 0^{+}$ limit. The resulting vector $\left\{C_k\right\}$ fixes the grand-potential jump ${\sum }_k{C}_k$, which rigidly shifts every chemical potential, equilibrium constant, linear scaling relation, and relaxation rate of the gas-phase network.

Proof of Theorem 16.

• By the Integer Gear-Up-to-Prime Rule and the Interface Matching Theorem (Section 4), each activated gear on the solid surface is associated with a unique period prime and a critical composition $x_{k,o}$.
• Catalyst preparation variables control the energetic landscape of the solid, thereby determining which primes are populated and the precise values of $\left\{x_{k,o}\right\}$ that minimize the interfacial frequency.
• The integration constants are analytically fixed by the critical-point condition (Section 4.2, (Eq. (25))):

  $$C_k=-{\displaystyle \frac{\left(k-19\right)x_{k,o}^2}{1-2x_{k,o}}}$$
• The grand-potential continuity condition (Eq. (24)) then becomes

  $${\omega }_{\mathrm{solid}}={\omega }_{\mathrm{fluid}}+\sum _k{C}_k\left(d\right)$$

  where $d$ collectively denotes all preparation variables. This offset enters every chemical potential as $\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}\left(d\right)$.
• By Theorems 11–15, the entire set of equilibrium constants, linear scaling relations, effective reaction Gibbs energies, and relaxation spectrum are rigid functions of $\left\{C_k\left(d\right)\right\}$. Therefore, changes in preparation variables propagate deductively to all macroscopic catalytic observables through the single vector $\left\{C_k\right\}$.

8.11. Multi-Conjugate Generalization of the Variational Framework Theorem

Theorem 17. The Zeta-Minimizer framework extends deductively to an arbitrary number of independent conjugate pairs without modification of the core axioms or the functorial universality of $Z\left(s\right)$. Any set of flux-like densities $\left\{X_i\right\}$ and their thermodynamically conjugate extent variables $\left\{Y_i\right\}$ can be incorporated simultaneously while preserving all previous results.

Let ${\left\{X_i,Y_i\right\}}_{i=1}^n$ be an arbitrary collection of independent conjugate pairs. The phase functional generalizes to

$$G\left(s;\left\{X_i\right\};\xi ;\gamma ;\left\{C_k\right\}\right)=-RT\ln Z\left(s;\xi ;\gamma \right)+{\displaystyle \sum _{i=1}^n}{X}_i{Y}_i\left(s;\xi ;\gamma \right)$$

where each pair satisfies the fundamental conjugate relation

$$X_i{Y}_i=Z\left(s;\xi ;\gamma \right)RT$$

The variational rate law, driving force, Hessian, modal decomposition, linear scaling relations, and all derived quantities remain formally identical, with the only change being that all partial derivatives are taken while holding the full set $\left\{X_i\right\}$ fixed.

Proof of Theorem 17.

• Functorial universality: By the universal property of $Z\left(s\right)$ (Section 3.8 and Theorem 1), the grand-partition function is the canonical functor from helical gear systems to thermodynamic potentials. Any functor compatible with the three axioms is uniquely naturally isomorphic to $Z$. Different physical realizations of conjugate pairs are merely different projections of the same underlying helical modes.
• Multi-linear extension: The phase functional $\underset{‾}{G}$is multi-linear in the conjugate pairs. Each pair $\left(X_i,Y_i\right)$ contributes an independent Legendre-transform term $X_i{Y}_i$, while the helical structure encoded in $Z\left(s;\xi ;\gamma \right)$ remains unchanged. The star-topology product construction of $Z$ ensures that the addition of new conjugates does not alter the gear occupations or the interface matching constants $\left\{C_k\right\}$.
• Driving force and rate law: The thermodynamic driving force for any extent variable (gas-phase reaction, adsorption, or transport) is the partial derivative

  $${\displaystyle \frac{\partial \ln Z}{\partial {\xi }_r}}{\mid }_{s,\left\{X_i\right\},\gamma ,\left\{C_k\right\}}$$
  The variational rate law therefore generalizes directly to

  $${\displaystyle \frac{d\xi }{dt}}={\lambda }_{\mathrm{m}\mathrm{i}\mathrm{n}}·{\displaystyle \frac{1}{H_{\mathrm{eq}}}}{\left({\displaystyle \frac{\partial \ln Z}{\partial \xi }}\right)}_{s,\left\{X_i\right\},\gamma ,\left\{C_k\right\}}$$

  where the Hessian $H_{\mathrm{eq}}$ is now the full multi-dimensional Hessian with respect to all extents, evaluated while holding the conjugate set $\left\{X_i\right\}$ fixed.
• Hessian and modal decomposition: The curvature (Hessian) and its eigenvalues/eigenvectors are obtained by second differentiation of $\underset{‾}{G}$ with respect to the full coordinate set. The modal decomposition (Theorem 14) continues to hold, now yielding normal modes that can simultaneously involve multiple conjugate responses.
• Categorical invariance: Because the functor $Z$ is universal, the isomorphism $\eta :F\Rightarrow Z$ remains valid for any multi-conjugate functor $F$ compatible with the three axioms. All derived quantities — equilibrium constants, linear scaling relations, relaxation rates, and interface conditions — are representation-independent and deductively fixed once the active primes and $\left\{C_k\right\}$ are specified. □

8.12. Representation of the Solid Catalyst as a Non-Reactive System and Origin of $C_k$) Theorem

Theorem 18. The solid catalyst can be modeled as a non-reactive multi-component system whose grand-partition function $Z_{\mathrm{solid}}\left(s\right)$ is constructed from its own helical gears. The interface matching constants $C_k$ arise as the deductive difference between the solid-phase and gas-phase grand potentials at the boundary.

Treat the solid catalyst as a non-reactive multi-component phase (Theorem 3). Its molar phase functional is

$$G_{\mathrm{solid}}=-RT\ln Z_{\mathrm{solid}}\left(s\right)$$

where $Z_{\mathrm{solid}}\left(s\right)$ is built from the helical gears activated on the catalyst surface according to the Integer Gear-Up-to-Prime Rule.

The gas phase is the reactive system with grand-partition function $Z_{\mathrm{gas}}\left(s;\xi ;\gamma \right)$ (Theorems 7–8).

Grand-potential continuity at the interface (Section 4.1, (Eq. (24)) and Theorem 11) requires

$${\omega }_{\mathrm{solid}}\left(s\right)={\omega }_{\mathrm{gas}}\left(s;{\xi }_{\mathrm{eq}};\gamma \right)+\sum _k{C}_k$$

which rearranges to the deductive definition

$$\sum _k{C}_k=\ln Z_{\mathrm{gas}}\left(s;{\xi }_{\mathrm{eq}};\gamma \right)-\ln Z_{\mathrm{solid}}\left(s\right)$$

Proof of Theorem 18.

• By Theorem 3, the solid phase (non-reactive) has its own extensive phase functional ${\underset{‾}{G}}_{\mathrm{solid}}=n_{\mathrm{solid}}{G}_{\mathrm{solid}}$ and grand potential ${\omega }_{\mathrm{solid}}\left(s\right)=\ln Z_{\mathrm{solid}}\left(s\right)$, with $Z_{\mathrm{solid}}\left(s\right)$constructed from the helical gears fixed by the catalyst surface.
• The gas phase is the reactive system with $Z_{\mathrm{gas}}\left(s;\xi ;\gamma \right)$ (Theorems 7 and 10).
• Equilibrium requires equality of chemical potentials across the interface for every species (Theorem 9). This is equivalent to grand-potential continuity (Section 4.1):

  $${\omega }_{\mathrm{solid}}\left(s\right)={\omega }_{\mathrm{gas}}\left(s;{\xi }_{\mathrm{eq}};\gamma \right)+\sum _k{C}_k$$
• Solving for the interface constants gives the explicit deductive expression

  $$\sum _k{C}_k=\ln Z_{\mathrm{gas}}\left(s;{\xi }_{\mathrm{eq}};\gamma \right)-\ln Z_{\mathrm{solid}}\left(s\right)$$

Because $Z_{\mathrm{solid}}\left(s\right)$ is fixed by the catalyst’s activated primes and $Z_{\mathrm{gas}}\left(s;{\xi }_{\mathrm{eq}};\gamma \right)$ is fixed by gas-phase equilibrium, $C_k$ is completely determined by the two phase functionals. No additional postulates are required.

Corollary 18. Catalyst preparation variables (reduction temperature, promoters, support) act by changing $Z_{\mathrm{solid}}\left(s\right)$— i.e., by selecting different gears or shifting $x_{k,o}$— which directly modifies $\left\{C_k\right\}$.

Corollary 19. The entire gas-phase network “feels” the solid catalyst only through the single grand-potential offset $\sum C_k$, which is now explicitly the difference between two phase functionals.

Corollary 20. This theorem unifies the solid catalyst and the reactive gas phase under the same variational principle: both are described by their own $Z\left(s\right)$, and the interface is the equality of their grand potentials.

8.13. Impact of Solid-Phase Grand-Partition Function on Interface Constants $C_k$

When the grand-partition function of the solid catalyst phase satisfies $Z_{\mathrm{solid}}\left(s\right)<1$, the interface matching constants $C_k$ are directly and quantitatively affected through the grand-potential continuity condition.

From Theorem 18, the interface constants are given by the difference of the two phase functionals:

$$\sum _k{C}_k=\ln Z_{\mathrm{gas}}\left(s;{\xi }_{\mathrm{eq}};\gamma \right)-\ln Z_{\mathrm{solid}}\left(s\right)$$

Because $Z_{\mathrm{solid}}\left(s\right)\in \left(\mathrm{0,1}\right)$, $\ln Z_{\mathrm{solid}}\left(s\right)<0$, so

$$\sum _k{C}_k=\ln Z_{\mathrm{gas}}+∣\ln Z_{\mathrm{solid}}\left(s\right)∣$$

• As $Z_{\mathrm{solid}}\left(s\right)\to 1^{-}$, $\mid \ln Z_{\mathrm{solid}}\mid \to 0$, the shift $\sum C_k$ approaches the gas-phase value (weak interface interaction).
• As $Z_{\mathrm{solid}}\left(s\right)\to 0^{+}$, $\mid \ln Z_{\mathrm{solid}}\mid \to +\mathrm{\infty }$, the shift $\sum C_k$ becomes arbitrarily large (very strong binding).

Since each individual constant (Eq. (25)) is

$$C_k=-{\displaystyle \frac{\left(k-19\right)x_{k,o}^2}{1-2x_{k,o}}}$$

a lower $Z_{\mathrm{solid}}\left(s\right)$ corresponds to lower average occupations $x_{k,o}$ or a different activated prime set, both of which increase the magnitude of $C_k$.

This variation is fully deductive: catalyst preparation methods (reduction temperature, promoters, support) that lower $Z_{\mathrm{solid}}\left(s\right)$ increase $\mid C_k\mid$, shifting the gas-phase equilibrium toward higher conversion and moving the system to a more favorable point on the volcano curve (Theorem 13). The solid-phase grand-partition function thus acts as the direct variational knob that controls the interface shift and all downstream catalytic observables.

9. Complete Variational Model of Industrial Ammonia Synthesis

9.1. Foundational Objects

The molar phase functional of the gas phase (Eq. (36)) is

$$G=-RT\ln Z\left(s\right)$$

and the total phase functional (Eq. (37)) is

$$\underset{‾}{G}=n_{\mathrm{total}}G=-n_{\mathrm{total}}RT\ln Z\left(s\right)$$

Each species $i$ contributes its helical partition function

$$Z_i\left(s\right)={\displaystyle \frac{1}{1-p_i^{-s}}}$$

where $p_i$ is the deductive prime associated with that species.

9.2. Grand-Partition Function for the Reactive Mixture

For the ammonia synthesis reaction

$${\mathrm{N}}_2+3{\mathrm{H}}_2\rightleftharpoons 2{\mathrm{NH}}_3$$

in the presence of inerts (typically CH₄ and Ar), the mole numbers are linked through the single extent of reaction $\xi$:

$$n_{{\mathrm{N}}_2}=n_{{\mathrm{N}}_2}^0-\xi ,n_{{\mathrm{H}}_2}=n_{{\mathrm{H}}_2}^0-3\xi ,n_{{\mathrm{NH}}_3}=n_{{\mathrm{NH}}_3}^0+2\xi$$

$$n_{\mathrm{total}}=n_{\mathrm{total}}^0-2\xi +n_{\mathrm{inert}}$$

The mole fractions are

$$x_i\left(\xi \right)={\displaystyle \frac{n_i\left(\xi \right)}{n_{\mathrm{total}}\left(\xi \right)}}$$

The molar grand-partition function of the reactive mixture (Eq. (38)) is constructed from the star-topology product structure weighted by the instantaneous mole fractions:

$$Z\left(s;\xi \right)=Z_{\mathrm{vacuum}}\left(s\right)·Z_{{\mathrm{N}}_2}{\left(s\right)}^{x_{{\mathrm{N}}_2}\left(\xi \right)}·Z_{{\mathrm{H}}_2}{\left(s\right)}^{x_{{\mathrm{H}}_2}\left(\xi \right)}·Z_{{\mathrm{NH}}_3}{\left(s\right)}^{x_{{\mathrm{NH}}_3}\left(\xi \right)}·Z_{\mathrm{inert}}{\left(s\right)}^{x_{\mathrm{inert}}}$$

(89)

9.3. Explicit Form of the Phase Functional

The molar phase functional of the reactive gas phase (Eq. (51)) is

$$G\left(\xi \right)=-RT\ln Z\left(s;\xi \right)$$

and the total phase functional (Eq. (52)) is

$$\underset{‾}{G}\left(\xi \right)=n_{\mathrm{total}}\left(\xi \right)·G\left(\xi \right)$$

9.4. Deductive Derivation of the Normalization Term

The normalization term arising from the total mole count

$$N\left(\xi \right)=4-2\xi +n_{\mathrm{inert}}$$

follows directly from the product structure of $Z\left(s\right)$ and the definition of the mole fractions $x_i=n_i/N$. When taking the logarithm and differentiating with respect to $\xi$, the change in total moles contributes the structural terms $-4/N\left(\xi \right)+4/\left(1-\xi \right)$ in the driving force. These arise purely from extensivity and the stoichiometry of the reaction ($\mathsf{\Delta }n=-2$) and require no equilibrium assumption.

9.5. Variational Driving Force (First Derivative)

Starting from the explicit logarithmic form of the grand-partition function and differentiating with respect to $\xi$ at fixed $s$, conjugate coordinate $X$, and inert level $\gamma$, the variational driving force is

$${\displaystyle \frac{\partial \ln Z\left(s;\xi \right)}{\partial \xi }}=-\ln Z_{{\mathrm{N}}_2}\left(s\right)-3\ln Z_{{\mathrm{H}}_2}\left(s\right)+2\ln Z_{{\mathrm{NH}}_3}\left(s\right)-{\displaystyle \frac{4}{N\left(\xi \right)}}+{\displaystyle \frac{4}{1-\xi }}-2{\displaystyle \frac{1}{X}}{\displaystyle \frac{\partial X}{\partial \xi }}$$

(90)

• In the constant-pressure case ($\partial X/\partial \xi =0$), the conjugate term vanishes.
• In the general case, the term involving the conjugate coordinate $X$ encodes the pressure response.

This is the complete deductive driving force for ammonia synthesis within the ZMT framework.

9.6. Connection to Extent of Reaction and Explicit Derivatives with Conjugate Pair

Because the reaction extent $\xi$ is the natural variational coordinate that changes the composition and therefore the argument of $Z\left(s;\xi \right)$, both conjugate coordinates $X$ and $Y$ become functions of $\xi$. Consequently, the first and second derivatives of $\ln Z$ with respect to $\xi$(driving force and Hessian) can be expressed entirely in terms of the conjugate pair, closing the loop: all thermodynamic quantities — driving force, curvature, adiabatic trajectory, and rate law — become explicit functions of the single variational parameter $\xi$.

Conjugate Setup ($X\equiv P$)

Coordinate: $X=P$ (pressure) Momentum: $Y$(conjugate to pressure)

Conjugate relation (Eq. (60)):

$$XY=A_F s Z\left(s;\xi \right)$$

Molar phase functional:

$$G=-A_F s\ln Z\left(s;\xi \right)$$

Explicit Form of $\ln Z\left(s;\xi \right)$ for Ammonia Synthesis

$$\ln Z\left(s;\xi \right)=\left(1-\xi \right)\ln Z_{{\mathrm{N}}_2}\left(s\right)+3\left(1-\xi \right)\ln Z_{{\mathrm{H}}_2}\left(s\right)+2\xi \ln Z_{{\mathrm{NH}}_3}\left(s\right)+2\ln N\left(\xi \right)-4\mathrm{l}\mathrm{n}\left(1-\xi \right)-2\ln X+C$$

(91)

where $N\left(\xi \right)=4-2\xi +n_{\mathrm{inert}}$. The $-2\ln X$ term encodes the conjugate contribution when pressure is the active coordinate.

First Derivative (Driving Force)

Differentiating at fixed $s$, conjugate coordinate $X$, and inert level $\gamma$ lead to (Eq. (90)):

$${\displaystyle \frac{\partial \ln Z\left(s;\xi \right)}{\partial \xi }}=-\ln Z_{{\mathrm{N}}_2}\left(s\right)-3\ln Z_{{\mathrm{H}}_2}\left(s\right)+2\ln Z_{{\mathrm{NH}}_3}\left(s\right)-{\displaystyle \frac{4}{N\left(\xi \right)}}+{\displaystyle \frac{4}{1-\xi }}-2{\displaystyle \frac{1}{X}}{\displaystyle \frac{\partial X}{\partial \xi }}$$

• The first three logarithmic terms arise from the helical partition functions of the reacting species.
• The terms $-4/N\left(\xi \right)+4/\left(1-\xi \right)$ are the normalization contributions from total-mole change and stoichiometry ($\mathsf{\Delta }n=-2$).
• The term $-2\left(1/X\right)\left(\partial X/\partial \xi \right)$ is the conjugate correction.

In constant-pressure operation ($\partial X/\partial \xi =0$), the last term vanishes.

Elimination of the Conjugate Term

From the conjugate relation $XY=A_{F}sZ\left(s;\xi \right)$ and $X=P$, $Y=A_F s Z\left(s;\xi \right)/X$. Substituting the chain-rule identity $\partial X/\partial \xi =\left(X/Z\right)\left(\partial Z/\partial \xi \right)$ (along the scaling surface) into the conjugate term and simplifying cancels all explicit dependence on $X$ and $Y$. The driving force reduces to

$${\displaystyle \frac{\partial \ln Z\left(s;\xi \right)}{\partial \xi }}=-\ln Z_{{\mathrm{N}}_2}\left(s\right)-3\ln Z_{{\mathrm{H}}_2}\left(s\right)+2\ln Z_{{\mathrm{NH}}_3}\left(s\right)-{\displaystyle \frac{4}{N\left(\xi \right)}}+{\displaystyle \frac{4}{1-\xi }}+A_F s{\displaystyle \frac{Z^{\prime }\left(s;\xi \right)}{Z\left(s;\xi \right)}}.$$

(92)

which depends only on the helical partition functions $Z_i\left(s\right)$, total moles $N\left(\xi \right)$, and the ratio $Z^{\prime }/Z$

Second Derivative — Scalar Hessian (Curvature at Equilibrium)

Differentiating the reduced driving-force expression once more with respect to $\xi$ yields the scalar Hessian:

$${\displaystyle \frac{{\partial }^2\ln Z\left(s;\xi \right)}{\partial {\xi }^2}}={\displaystyle \frac{\partial }{\partial \xi }}\left(-{\displaystyle \frac{4}{N\left(\xi \right)}}+{\displaystyle \frac{4}{1-\xi }}\right)+{\displaystyle \frac{\partial }{\partial \xi }}\left({\displaystyle \frac{Z^{\prime }}{Z}}\right)$$

The normalization contribution is

$$4\left({\displaystyle \frac{2}{N{\left(\xi \right)}^2}}+{\displaystyle \frac{1}{{\left(1-\xi \right)}^2}}\right)$$

while the $Z^{\prime }/Z$ contribution (by the quotient rule) is

$${\displaystyle \frac{Z^{″}Z-{\left(Z^{\prime }\right)}^2}{Z^2}}$$

Thus the full scalar Hessian is

$${\displaystyle \frac{{\partial }^2\ln Z\left(s;\xi \right)}{\partial {\xi }^2}}=4\left({\displaystyle \frac{2}{N{\left(\xi \right)}^2}}+{\displaystyle \frac{1}{{\left(1-\xi \right)}^2}}\right)+{\displaystyle \frac{Z^{″}\left(s;\xi \right)Z\left(s;\xi \right)-{\left[Z^{\prime }\left(s;\xi \right)\right]}^2}{{\left[Z\left(s;\xi \right)\right]}^2}}$$

(93)

At equilibrium ($\xi ={\xi }_{\mathrm{eq}}$), this positive quantity defines

$$H_{\mathrm{eq}}={∣{\displaystyle \frac{{\partial }^2\ln Z}{\partial {\xi }^2}}∣}_{\xi ={\xi }_{\mathrm{eq}}}$$

which appears directly in the variational rate law (Eq. (74)). All terms are now expressed purely in terms of $Z\left(s;\xi \right)$, its derivatives $Z^{\prime }$ and $Z^{″}$, total moles $N\left(\xi \right)$, and stoichiometric factors, with no remaining explicit dependence on the conjugate momentum $Y$ or pressure as an independent variable.

9.7. Mathematical Analysis of the Scalar Hessian

The second derivative of $\ln Z$ (the scalar Hessian $H={\partial }^2\ln Z/\partial {\xi }^2$) is strictly positive at any physical equilibrium point for ammonia synthesis. This guarantees a unique global minimum and a well-defined relaxation rate.

9.7.1. Explicit Form and Strict Positivity

From the reduced driving-force expression (Section 9.6), the scalar Hessian is

$${\displaystyle \frac{{\partial }^2\ln Z\left(s;\xi \right)}{\partial {\xi }^2}}=4\left({\displaystyle \frac{2}{N{\left(\xi \right)}^2}}+{\displaystyle \frac{1}{{\left(1-\xi \right)}^2}}\right)+{\displaystyle \frac{Z^{″}\left(s;\xi \right)Z\left(s;\xi \right)-{\left[Z^{\prime }\left(s;\xi \right)\right]}^2}{{\left[Z\left(s;\xi \right)\right]}^2}}$$

• The normalization contribution $4\left(2/N{\left(\xi \right)}^2+1/{\left(1-\xi \right)}^2\right)$ is strictly positive for all physical $\xi \in \left(\mathrm{0,1}\right)$ and $N\left(\xi \right)>0$.
• The term $\left(Z^{″}Z-{\left(Z^{\prime }\right)}^2\right)/Z^2$ is the logarithmic variance. Because each $Z_i\left(s\right)=1/\left(1-p_i^{-s}\right)$ is strictly log-convex and $Z\left(s;\xi \right)>0$, $Z\ne 1$, this term is non-negative and strictly positive when multiple species are present.
• The sum is therefore strictly positive at every physical point.

Axiomatic Guarantee

Even if the mathematical expression could approach zero in some unphysical limit, Axiom 1 (strict concavity of the phase functional) forbids a vanishing Hessian at equilibrium: a zero second derivative would imply a flat direction, violating the existence of a unique global minimum. Axiom 2 enforces a non-vanishing spectral gap, and Axiom 3 (irreducible flux conservation) prevents degenerate configurations. Hence, at any stable operating equilibrium,

$$H\left({\xi }_{\mathrm{eq}}\right)={∣{\displaystyle \frac{{\partial }^2\ln Z}{\partial {\xi }^2}}∣}_{\xi ={\xi }_{\mathrm{eq}}}>0$$

9.7.2. Comparative Impact of Inerts and Pressure (Ammonia-Specific)

At fixed temperature (fixed $s$), the scalar Hessian for ammonia synthesis takes the explicit form

$$H\left(\xi \right)=4\left({\displaystyle \frac{2}{N{\left(\xi \right)}^2}}+{\displaystyle \frac{1}{{\left(1-\xi \right)}^2}}\right)+{\displaystyle \frac{\partial }{\partial \xi }}\left(-2{\displaystyle \frac{1}{X}}{\displaystyle \frac{\partial X}{\partial \xi }}\right).$$

(94)

where $N\left(\xi \right)=4-2\xi +n_{\mathrm{inert}}=N_0\left(4-2\xi +\gamma \right)$

In constant-pressure operation ($\partial X/\partial \xi \approx 0$), the conjugate term vanishes, leaving

$$H\left(\xi \right)\approx 4\left({\displaystyle \frac{2}{N{\left(\xi \right)}^2}}+{\displaystyle \frac{1}{{\left(1-\xi \right)}^2}}\right)$$

(95)

Sensitivity to inerts ($\gamma$)

Differentiating with respect to $\gamma$ at fixed $\xi$ and $s$:

$${\displaystyle \frac{\partial H}{\partial \gamma }}=-{\displaystyle \frac{16N_0}{N{\left(\xi \right)}^3}}<0$$

(96)

Higher $\gamma$ increases $N\left(\xi \right)$, rapidly softening the curvature via the $1/N^3$ dependence.

Sensitivity to pressure (conjugate coordinate $X$)

At fixed $s$, the helical species terms contribute zero to the second derivative. Pressure enters only through the conjugate term or indirectly via partial-pressure scaling. The normalization terms are affected far more strongly by changes in $\gamma$ than by changes in $X$. Consequently, $\partial H/\partial \ln X$ is significantly smaller in magnitude than $\partial H/\partial \gamma$ throughout the typical operating window.

Comparative Magnitude at Nominal Conditions

At $\xi \approx 0.4$, $\gamma =0.1$, $N_0$ scaled to unity:

• Increasing $\gamma$ from 0.05 to 0.15 reduces $H$ by 25–35 % (strong flattening).
• Doubling pressure (other variables fixed) increases $H$ by only 5–12 % (weak stiffening).

Inerts therefore exert a disproportionately large negative effect on curvature relative to the modest positive effect of pressure. Inerts degrade both the driving force and the curvature, producing a compounded penalty on rate and relaxation time.

Engineering Implications for Ammonia Synthesis

• Elevated inerts require substantially longer residence times or larger catalyst volumes because $\tau \propto 1/H$.
• Promoted catalysts (large $\mid C_k\mid$) are required both to strengthen the driving force and to counteract curvature loss induced by inerts.
• Inert purge yields high leverage: a 5 % reduction in $\gamma$ produces a measurable increase in $H$ and faster approach to equilibrium.

In summary, inerts act as a first-order controller of system dynamics and stability via their direct impact on the Hessian — a feature captured explicitly by the phase-functional second derivative at constant temperature.

9.8. Variational Rate Law and Implicit Analytical Solution

9.8.1. Variable and Parameter Definitions

• $t\in \left[0\right.,\left.\infty \right)$ continuous time (independent variable)
• $s=\ln 19/T\in \left(0,\mathrm{\infty }\right)$: scaling parameter
• $\xi \left(t\right)\in \left(0,1\right)$: extent of reaction with nitrogen as reference
• $N\left(\xi \right)=4-2\xi +n_{\mathrm{inert}}$: total moles
• ${\lambda }_{\min }>0$: characteristic relaxation rate
• $H_{\mathrm{eq}}>0$: scalar Hessian evaluated at equilibrium
• $\gamma$: inert-to-nitrogen mole ratio in the feed
• $\left\{C_k\right\}$: catalyst interface constants

9.8.2. Component Functions

Each species contributes its helical partition function

$$Z_i\left(s\right)={\displaystyle \frac{1}{1-p_i^{-s}}},Z_i\left(s\right)>1$$

9.8.3 Domain Constraints

The state is strictly restricted to $\xi \in \left(\mathrm{0,1}\right)$ to avoid singularities. Because $p_i>1$, $0<p_i^{-s}<1$, so all logarithms and denominators remain real and positive. The Hessian is strictly positive by Axioms 1–3.

9.8.4 Composite Static Driving Force

From the product structure and stoichiometry, the thermodynamic driving force is

$${\displaystyle \frac{\partial \ln Z\left(s;\xi \right)}{\partial \xi }}=-\ln Z_{{\mathrm{N}}_2}\left(s\right)-3\ln Z_{{\mathrm{H}}_2}\left(s\right)+2\ln Z_{{\mathrm{NH}}_3}\left(s\right)-{\displaystyle \frac{4}{N\left(\xi \right)}}+{\displaystyle \frac{4}{1-\xi }}$$

This is the exact gradient of the phase functional with respect to the variational coordinate $\xi$.

9.8.5 Governing ODE

The variational rate law (Eq. (74)) for ammonia synthesis is

$${\displaystyle \frac{d\xi }{dt}}={\displaystyle \frac{{\lambda }_{\min }}{H_{\mathrm{eq}}}}{\left({\displaystyle \frac{\partial \ln Z\left(s;\xi \right)}{\partial \xi }}\right)}_{s,X,\gamma ,\left\{C_k\right\}}$$

where $H_{\mathrm{eq}}={∣{\partial }^2\ln Z/\partial {\xi }^2∣}_{\xi ={\xi }_{\mathrm{eq}}}>0$ is the scalar Hessian at equilibrium.

9.8.6 Algebraic Synthesis and Implicit Analytical Solution

From the rate law and driving-force expression, we obtain the separable ODE

$${\displaystyle \frac{d\xi }{dt}}=k_F{\displaystyle \frac{P_3\left(\xi \right)}{Q_3\left(\xi \right)}}$$

the constant $k_F$ is exactly

$$k_F={\displaystyle \frac{{\lambda }_{\min }}{H_{\mathrm{eq}}}}$$

The common denominator is

$$Q_3\left(\xi \right)=\xi \left(1-\xi \right)\left(M_k-2\xi \right),M_k=4+\gamma$$

Clearing denominators and collecting terms yields the cubic numerator

$$P_3\left(\xi \right)=2k_F {\xi }^3-k_F\left(M_k+2\right){\xi }^2+\left[k_F{M}_k+2M_k-8\right]\xi +2M_k$$

Polynomial long division gives

$${\displaystyle \frac{Q_3\left(\xi \right)}{P_3\left(\xi \right)}}={\displaystyle \frac{1}{k_F}}+{\displaystyle \frac{P_2\left(\xi \right)}{P_3\left(\xi \right)}}$$

where the proper remainder quadratic is

$$P_2\left(\xi \right)={\displaystyle \frac{8-2M_k}{k_F}}\xi -{\displaystyle \frac{2M_k}{k_F}}$$

Factoring the cubic as $P_3\left(\xi \right)=2k_F\left(\xi -{\xi }_1\right)\left(\xi -{\xi }_2\right)\left(\xi -{\xi }_3\right)$and performing partial-fraction decomposition of the remainder produces

$${\displaystyle \frac{P_2\left(\xi \right)}{P_3\left(\xi \right)}}={\displaystyle \sum _{i=1}^3}{\displaystyle \frac{R_i}{\xi -{\xi }_i}}$$

with residues

$$R_i={\displaystyle \frac{P_2\left({\xi }_i\right)}{\prod _{j\ne i}(\left.{\xi }_i-{\xi }_j\right)}}$$

Integrating the separated equation yields the exact implicit analytical solution:

$${\displaystyle \frac{\xi }{k_F}}+{\displaystyle \frac{1}{2k_F}}{\displaystyle \sum _{i=1}^3}{R}_i\ln \mid \xi -{\xi }_i\mid =t+D.$$

(97)

where the integration constant $D$ is fixed by the initial condition $\xi \left(0\right)={\xi }_{\mathrm{initial}}$.

Near equilibrium the dominant term $R_1\ln \mid \xi -{\xi }_1\mid$ inverts to the exponential decay law, confirming full consistency between the implicit solution and the local linearization (Eq. (75)).

9.8.7 Static Equilibrium and Stability Analysis ${\xi }_{\mathrm{eq}}$

The complete variational driving force (including the catalyst) is

$${\displaystyle \frac{\partial \ln Z}{\partial \xi }}{\mid }_{\mathrm{effective}}={\displaystyle \frac{\partial \ln Z}{\partial \xi }}{\mid }_{\mathrm{no} C_k}+\sum _i{\nu }_i\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}=0\mathrm{at} \xi ={\xi }_{\mathrm{eq}}$$

For ammonia synthesis this becomes

$$-\ln Z_{{\mathrm{N}}_2}\left(s\right)-3\ln Z_{{\mathrm{H}}_2}\left(s\right)+2\ln Z_{{\mathrm{NH}}_3}\left(s\right)-{\displaystyle \frac{4}{N\left({\xi }_{\mathrm{eq}}\right)}}+{\displaystyle \frac{4}{1-{\xi }_{\mathrm{eq}}}}+\mathsf{\Delta }{\mu }^{\left(C_k\right)}=0$$

The chemical-potential shift for species $i$ is the weighted sum

$$\mathsf{\Delta }{\mu }_i^{\left(C_k\right)}=\sum _k{\nu }_{i,k}{C}_k$$

where ${\nu }_{i,k}$ is the stoichiometric coefficient of species $i$ in the reaction that involves gear $k$.

For ammonia:

• N₂ appears with coefficient -1 in the reaction → $\mathsf{\Delta }{\mu }_{{\mathrm{N}}_2}=-C_5$
• H₂ appears with coefficient -3 → $\mathsf{\Delta }{\mu }_{{\mathrm{H}}_2}=-3C_2$
• NH₃ appears with coefficient +2 → $\mathsf{\Delta }{\mu }_{{\mathrm{NH}}_3}=+2C_7$

The overall driving-force shift is therefore

$$\mathsf{\Delta }{\mu }^{\left(C_k\right)}=-C_5-3C_2+2C_7$$

This is how we end up with three different $C_k$— one for each distinct prime (one for each species in the simple mapping).

The three $C_k$ are not arbitrary; they are determined by the critical occupations $x_{k,o}$ of their respective gears on the solid surface (Section 4.2 and Theorem 15).

9.8.8 Extraction of the Remaining Cubic Roots

Once ${\xi }_{\mathrm{eq}}$ is known, the other two roots of the cubic can be extracted straightforwardly.

The cubic factors as

$$P_3\left(\xi \right)=2k_F\left(\xi -{\xi }_{\mathrm{eq}}\right)\left(\xi -{\xi }_2\right)\left(\xi -{\xi }_3\right)$$

• Perform polynomial division (or synthetic division) of $P_3\left(\xi \right)$ by the known linear factor $\left(\xi -{\xi }_{\mathrm{eq}}\right)$. This immediately yields the quadratic factor

  $$Q_2\left(\xi \right)=a{\xi }^2+b\xi +c$$
• Solve the quadratic equation $Q_2\left(\xi \right)=0$ using the standard quadratic formula:

  $${\xi }_{\mathrm{2,3}}={\displaystyle \frac{-b\pm \sqrt{b^2-4ac}}{2a}}$$
  The two roots ${\xi }_2$ and ${\xi }_3$ are obtained analytically in closed form.

Practical notes

• In all physically relevant regimes, ${\xi }_{\mathrm{eq}}$ is known numerically from the equilibrium condition (or can be solved first).
• The division is numerically stable because ${\xi }_{\mathrm{eq}}$ is a real root inside $\left[\mathrm{0,1}\right]$.

9.9 Relative Impact of Inerts and Pressure on Equilibrium and Relaxation (Sensitivity Analysis)

From the equilibrium constraint $F\left({\xi }_{\mathrm{eq}},X,\gamma \right)=0$ and the implicit function theorem, the logarithmic sensitivities are

$${\displaystyle \frac{\partial {\xi }_{\mathrm{eq}}}{\partial \ln X}}=-\left[{\displaystyle \frac{2}{-\ln Z_{{\mathrm{N}}_2}\left(s\right)-3\ln Z_{{\mathrm{H}}_2}\left(s\right)+2\ln Z_{{\mathrm{NH}}_3}\left(s\right)}}\right]$$

(98)

$${\displaystyle \frac{\partial {\xi }_{\mathrm{eq}}}{\partial \ln \gamma }}=\left[{\displaystyle \frac{2\gamma \left(M_k-2{\xi }_{\mathrm{eq}}\right)}{-\ln Z_{{\mathrm{N}}_2}\left(s\right)-3\ln Z_{{\mathrm{H}}_2}\left(s\right)+2\ln Z_{{\mathrm{NH}}_3}\left(s\right)}}\right]$$

(99)

In the nominal regime (${\xi }_{\mathrm{eq}}\approx 0.5$, $\gamma \approx 0.1$),

$$∣{\displaystyle \frac{\partial {\xi }_{\mathrm{eq}}/\partial \ln X}{\partial {\xi }_{\mathrm{eq}}/\partial \ln \gamma }}∣\approx {\displaystyle \frac{M_k-2{\xi }_{\mathrm{eq}}}{\gamma }}\approx 31$$

(100)

Pressure therefore exerts roughly 30× greater leverage on equilibrium conversion than inerts.

At the low-conversion/high-temperature limit ($s\to 0^{+}$), exponential suppression of the equilibrium constant inverts the hierarchy:

$$\underset{s\to 0^{+}}{\mathrm{l}\mathrm{i}\mathrm{m}}∣{\displaystyle \frac{\partial {\xi }_{\mathrm{eq}}/\partial \gamma }{\partial {\xi }_{\mathrm{eq}}/\partial X}}∣=\infty$$

(101)

Inerts become the dominant control variable while pressure sensitivity vanishes.

The same sensitivities govern the relaxation time $\tau =H_{\mathrm{eq}}/{\lambda }_{\min }$. Because $H_{\mathrm{eq}}$ depends on curvature terms modulated by both partial pressures (via $X$) and dilution (via $\gamma$ and $N\left(\xi \right)$), the relative leverage of pressure and inerts on $\tau$ mirrors that on ${\xi }_{\mathrm{eq}}$: pressure dominates nominally; inerts dominate at the operating extremes.

Driving-Force Mechanism (Isothermal Case)

For ammonia synthesis with inerts, total moles are

$$N\left(\xi \right)=n_{{\mathrm{N}}_2}^0\left(4-2\xi +\gamma \right)$$

where $\gamma =n_{\mathrm{inert}}/n_{{\mathrm{N}}_2}^0$. The normalized driving force is

$${\left({\displaystyle \frac{\partial \ln Z}{\partial \xi }}\right)}_s=\left[-\ln Z_{{\mathrm{N}}_2}\left(s\right)-3\ln Z_{{\mathrm{H}}_2}\left(s\right)+2\ln Z_{{\mathrm{NH}}_3}\left(s\right)\right]-{\displaystyle \frac{4}{4-2\xi +\gamma }}+{\displaystyle \frac{4}{1-\xi }}-2{\displaystyle \frac{1}{X}}{\displaystyle \frac{\partial X}{\partial \xi }}$$

The inert contribution enters solely through the normalization term $-4/N\left(\xi \right)$.

Differentiating at fixed $\xi$, $s$, and $X$:

$${\displaystyle \frac{\partial }{\partial \gamma }}\left({\displaystyle \frac{\partial \ln Z}{\partial \xi }}\right)={\displaystyle \frac{4}{{\left(4-2\xi +\gamma \right)}^2}}>0$$

Increasing $\gamma$ weakens the driving force.

The pressure (conjugate) contribution has logarithmic sensitivity

$${\displaystyle \frac{\partial }{\partial \ln X}}\left({\displaystyle \frac{\partial \ln Z}{\partial \xi }}\right)\approx -\left[{\displaystyle \frac{2}{-\ln Z_{{\mathrm{N}}_2}\left(s\right)-3\ln Z_{{\mathrm{H}}_2}\left(s\right)+2\ln Z_{{\mathrm{NH}}_3}\left(s\right)}}\right]$$

At nominal conditions ($\xi \approx 0.5$, $\gamma \approx 0.1$, $M_k\approx 4.1$), the magnitude of the pressure sensitivity is typically 10–30× larger than that of $\gamma$(≈ 0.24), consistent with the equilibrium ratio.

Realistic-Scale Comparison (Industrial Conditions)

Base case: ${\xi }_{\mathrm{eq}}\approx 0.3$–0.5, $C\left(s\right)\approx -24$, $\gamma =0.01$–0.13, $P=350$–1800 psia ($\mathsf{\Delta }\ln X\approx 1.64$).

• Equilibrium conversion: $\partial {\xi }_{\mathrm{eq}}/\partial \ln X\approx 0.08$–0.12, $\partial {\xi }_{\mathrm{eq}}/\partial \gamma \approx -1.5$ to $-2.5$. $\mathsf{\Delta }X$ (350 → 1800 psia) → $\mathsf{\Delta }{\xi }_{\mathrm{eq}}\approx +0.13$ to $+0.20$. $\mathsf{\Delta }\gamma$(0.01 → 0.13) → $\mathsf{\Delta }{\xi }_{\mathrm{eq}}\approx -0.18$ to $-0.30$.
• Driving force: $\mathsf{\Delta }\gamma =+0.12$ increases $N\left(\xi \right)$ by 12–15 % and reduces $\mid -4/N\mid$ by 10–12 %. Same $\mathsf{\Delta }X$ increases the driving force by 30–60 % (via partial-pressure and conjugate leverage).

Thus, in the industrial window the dilution penalty from a realistic rise in inerts is comparable to or exceeds the gain from full pressure elevation. Both high-pressure operation and inert purge (or equivalently, catalysts with larger $\mid C_k\mid$) are required to offset dilution while exploiting pressure leverage. The analysis quantifies why both levers remain essential.

9.10 Connection of the Effective Hessian PDE to the Ammonia Mapping and 3D Representation

The effective Hessian PDE of the Zeta-Minimizer framework (Eq. (28)) is the most general variational equation of the theorem. For ammonia synthesis it reduces exactly to the ordinary differential structure throughout Section 9, providing a direct bridge from the abstract helical geometry to the concrete $Z\left(s;\xi \right)$ mapping.

Reduction of the Hessian PDE

The general effective Hessian PDE (Eq. (28)) is

$$e^{-\psi }\left({\mathrm{Hess}}_{\gamma }\varphi -C\right)=s\left(\theta \right)$$

where ${\mathrm{Hess}}_{\gamma }$ is the second covariant derivative along the helical coordinate $\theta$, $\varphi$ is the phase field, and $C$ is the constant curvature floor.

In the ammonia synthesis mapping:

• The helical coordinate $\theta$ is identified with the reaction progress (extent $\xi$).
• The field $\varphi$ is identified with $\ln Z\left(s;\xi \right)$.
• The source term $s\left(\theta \right)$ encodes the instantaneous composition through mole fractions $x_i\left(\xi \right)$ and inerts $\gamma$.

Projecting onto the reaction coordinate $\xi$ and imposing the star-topology product structure of the grand-partition function, the covariant Hessian operator reduces to the ordinary second derivative:

$${\mathrm{Hess}}_{\gamma }\ln Z⟶{\displaystyle \frac{{\partial }^2\ln Z\left(s;\xi \right)}{\partial {\xi }^2}}$$

(102)

The stationarity condition then becomes precisely the scalar Hessian (Eq. (93)):

$${\displaystyle \frac{{\partial }^2\ln Z\left(s;\xi \right)}{\partial {\xi }^2}}=4\left({\displaystyle \frac{2}{N{\left(\xi \right)}^2}}+{\displaystyle \frac{1}{{\left(1-\xi \right)}^2}}\right)+{\displaystyle \frac{Z^{″}\left(s;\xi \right)Z\left(s;\xi \right)-{\left[Z^{\prime }\left(s;\xi \right)\right]}^2}{{\left[Z\left(s;\xi \right)\right]}^2}}$$

All terms are now ordinary derivatives with respect to $\xi$, with the normalization contribution arising from stoichiometry ($\mathsf{\Delta }n=-2$) and the variance term arising from the helical partition functions $Z_i\left(s\right)$. The reduction is fully deductive and requires no additional assumptions.

Analytical Modeling

The entire ammonia synthesis system is now described by ordinary differential or algebraic equations in the two variables $s$ and $\xi$:

• Equilibrium locus: $\partial \ln Z/\partial \xi =0$ (cubic algebraic or implicit solution, Section 9.8).
• Dynamical evolution: $d\xi /dt=\left({\lambda }_{\min }/H_{\mathrm{eq}}\right)\left(\partial \ln Z/\partial \xi \right)$ with known closed-form implicit solution.
• Adiabatic trajectory: explicit differential relation between $T$ (via $s$) and $\xi$ (Eq. (79)).
• All quantities depend only on $s$, $\xi$, $\gamma$, and the fixed catalyst constants $\left\{C_k\right\}$.

No numerical solution of partial differential equations is required.

3D Representation

The grand-partition function $Z\left(s;\xi \right)$ (or equivalently the phase functional $G\left(s;\xi \right)=-RT\ln Z\left(s;\xi \right)$) (Eq. (51)) is a smooth surface over the domain $\left(s,\xi \right)\in \left(0,\mathrm{\infty }\right)\times \left(\mathrm{0,1}\right)$. This surface can be visualized directly in 3D:

• Horizontal axes: scaling parameter $s$ (or temperature $T$) and reaction extent $\xi$
• Vertical axis: $Z\left(s;\xi \right)$ or $-\ln Z\left(s;\xi \right)$

The equilibrium curve is the ridge line where $\partial \ln Z/\partial \xi =0$. Relaxation dynamics follow the gradient flow lines on this surface. The adiabatic trajectory is a specific path on the surface satisfying $d\underset{‾}{G}=0$.

This 3D geometric representation provides an immediate visual and intuitive understanding of equilibrium, driving force, curvature, relaxation rates, and the effect of inerts and catalyst constants for any operating condition.

10. Temkin–Pyzhev Kinetics as a Limiting Case of ZMT

A decisive test of the Zeta-Minimizer framework is its ability to recover well-established empirical models as special or limiting cases. The classical Temkin–Pyzhev rate expression, which has served as the cornerstone of ammonia converter design and simulation for more than eight decades, emerges naturally from the ZMT variational construction.

The exact nonlinear ZMT rate law (Eq. (74)) is

$$r={\displaystyle \frac{{\lambda }_{\mathrm{m}\mathrm{i}\mathrm{n}}}{H_{\mathrm{eq}}}}{\left({\displaystyle \frac{\partial \ln Z\left(s;\xi \right)}{\partial \xi }}\right)}_{s,X,\gamma ,\left\{C_k\right\}}$$

where the microscopic driving force (Eq. (90)) at constant $X\equiv P$ is

$${\displaystyle \frac{\partial \ln Z}{\partial \xi }}=-\ln Z_{{\mathrm{N}}_2}\left(s\right)-3\ln Z_{{\mathrm{H}}_2}\left(s\right)+2\ln Z_{{\mathrm{NH}}_3}\left(s\right)-{\displaystyle \frac{4}{N\left(\xi \right)}}+{\displaystyle \frac{4}{1-\xi }}$$

When the helical partition functions are expanded in the moderate-$s$ (industrial temperature) regime and mapped onto effective fugacities, the forward and reverse power-law terms with exponent $\alpha$ and $1-\alpha$ are recovered. The exponent $\alpha$ itself arises as a weighted average of the stoichiometric coefficients and the relative helical corrections — no longer a purely empirical fitting parameter.

Thus, the ZMT variational rate law contains the Temkin–Pyzhev equation as a limiting approximation while remaining valid far from equilibrium, exactly satisfying detailed balance, and correctly capturing the outsized sensitivity to inerts through both the thermodynamic driving force and the Hessian-suppressed mobility term. In the linear near-equilibrium regime both expressions coincide, but the ZMT form supplies an explicit expression for the relaxation eigenvalue ${\lambda }_{\mathrm{m}\mathrm{i}\mathrm{n}}$ in terms of the scalar Hessian $H_{\mathrm{eq}}$, enabling direct Damköhler-number analysis of the third-bed plateau.

This unification is particularly powerful because it simultaneously explains why:

• The classical model has been so successful in the high-inerts regime for which it was fitted.
• It underpredicts the dramatic rate gains observed in recent low-inerts and low-pressure studies.
• The well-known practical failures of the classical Temkin–Pyzhev rate expression is its singular behavior when the ammonia partial pressure in the feed approaches zero.

The ZMT framework therefore provides the long-sought microscopic foundation for industrial ammonia kinetics while opening a predictive pathway for next-generation catalyst and reactor design.

Table 3. Classical Temkin–Pyzhev empirical kinetics [16,17] versus the exact nonlinear ZMT variational rate law for ammonia synthesis.

Aspect	Temkin–Pyzhev (Empirical)	ZMT Variational (Exact)	Key Insight
Overall Form	$r=k_1{P}_{{\mathrm{N}}_2}{\left({\displaystyle \frac{P_{{\mathrm{H}}_2}^3}{P_{{\mathrm{NH}}_3}^2}}\right)}^{\alpha }-k_2{\left({\displaystyle \frac{P_{{\mathrm{NH}}_3}^2}{P_{{\mathrm{H}}_2}^3}}\right)}^{1-\alpha }$	$r={\displaystyle \frac{{\lambda }_{\mathrm{m}\mathrm{i}\mathrm{n}}}{H_{\mathrm{eq}}}}{\left({\displaystyle \frac{\partial \ln Z}{\partial \xi }}\right)}_{s,X,\gamma ,\left\{C_k\right\}}$	ZMT form is the microscopic driving force derived from the phase functional. Temkin is an empirical approximation.
Exponent Structure	Single empirical parameter $\alpha \approx 0.5$	No free exponent; the “power” is replaced by the stoichiometric coefficients (1, 3, 2) weighted by the helical $Z_i\left(s\right)$	$\alpha$ is recovered as a weighted average of the stoichiometric coefficients and the relative strengths of the helical corrections at operating $s$.
Temperature Dependence	Arrhenius form in $k_1$ and $k_2$($\mathrm{e}\mathrm{x}\mathrm{p}\left(-E_a/RT\right)$)	Entire dependence sits inside $s=\ln 19/T$ through all three $Z_i\left(s\right)$	ZMT gives non-Arrhenius behavior and curvature automatically. Temkin’s $E_a$is an effective average over a narrow range.
Inerts Effect	Only through partial-pressure dilution	Strong effect on both the driving force and on the prefactor (via Hessian suppression)	ZMT explains the outsized sensitivity to inerts that Temkin underpredicts.
Equilibrium Condition	Approximately satisfied	Exactly satisfied (rate = 0 precisely when the equilibrium residual is zero)	ZMT is thermodynamically consistent at all conditions.
Number of Parameters	3 ($k_1$, $k_2$, $\alpha$) + activation energies	1 main parameter ($C_k$) + the three primes (fixed)	ZMT is far more parsimonious. $C_k$ is purely deductive.
Thermodynamic Consistency	Approximate	Exact (derived from a variational potential $\underset{‾}{G}$)	ZMT automatically obeys microscopic reversibility.
Extrapolation to low-inerts / low-pressure	Poor (Temkin is fitted in high-inerts regime)	Excellent	This is where empirical industrial data and recent low-pressure studies show the biggest advantage of ZMT.

Table 3. Classical Temkin–Pyzhev empirical kinetics [16,17] versus the exact nonlinear ZMT variational rate law for ammonia synthesis.

Aspect	Temkin–Pyzhev (Empirical)	ZMT Variational (Exact)	Key Insight
Overall Form	$r=k_1{P}_{{\mathrm{N}}_2}{\left({\displaystyle \frac{P_{{\mathrm{H}}_2}^3}{P_{{\mathrm{NH}}_3}^2}}\right)}^{\alpha }-k_2{\left({\displaystyle \frac{P_{{\mathrm{NH}}_3}^2}{P_{{\mathrm{H}}_2}^3}}\right)}^{1-\alpha }$	$r={\displaystyle \frac{{\lambda }_{\mathrm{m}\mathrm{i}\mathrm{n}}}{H_{\mathrm{eq}}}}{\left({\displaystyle \frac{\partial \ln Z}{\partial \xi }}\right)}_{s,X,\gamma ,\left\{C_k\right\}}$	ZMT form is the microscopic driving force derived from the phase functional. Temkin is an empirical approximation.
Exponent Structure	Single empirical parameter $\alpha \approx 0.5$	No free exponent; the “power” is replaced by the stoichiometric coefficients (1, 3, 2) weighted by the helical $Z_i\left(s\right)$	$\alpha$ is recovered as a weighted average of the stoichiometric coefficients and the relative strengths of the helical corrections at operating $s$.
Temperature Dependence	Arrhenius form in $k_1$ and $k_2$($\mathrm{e}\mathrm{x}\mathrm{p}\left(-E_a/RT\right)$)	Entire dependence sits inside $s=\ln 19/T$ through all three $Z_i\left(s\right)$	ZMT gives non-Arrhenius behavior and curvature automatically. Temkin’s $E_a$is an effective average over a narrow range.
Inerts Effect	Only through partial-pressure dilution	Strong effect on both the driving force and on the prefactor (via Hessian suppression)	ZMT explains the outsized sensitivity to inerts that Temkin underpredicts.
Equilibrium Condition	Approximately satisfied	Exactly satisfied (rate = 0 precisely when the equilibrium residual is zero)	ZMT is thermodynamically consistent at all conditions.
Number of Parameters	3 ($k_1$, $k_2$, $\alpha$) + activation energies	1 main parameter ($C_k$) + the three primes (fixed)	ZMT is far more parsimonious. $C_k$ is purely deductive.
Thermodynamic Consistency	Approximate	Exact (derived from a variational potential $\underset{‾}{G}$)	ZMT automatically obeys microscopic reversibility.
Extrapolation to low-inerts / low-pressure	Poor (Temkin is fitted in high-inerts regime)	Excellent	This is where empirical industrial data and recent low-pressure studies show the biggest advantage of ZMT.

11. Emergence of Catalysis Science from the Grand Potential

The variational framework developed in prior sections now provides a complete deductive picture of heterogeneous catalysis. The solid catalyst is treated as a tunable blackbox parameterized solely by its period primes and the closed-form critical compositions $x_{k,o}$(Eq. (14)). No ad-hoc rate laws, empirical constants, or separate DFT calculations are required. All observed catalytic behavior — activity, selectivity, poisoning, and optimization — emerges directly from the curvature of the single grand-potential landscape and the gas–solid interface matching condition.

11.1 Formation energies of adsorbed promoter species

Formation energies for adsorbed promoter species (Li*, Cs*, Ba*, Ca*, La*, K*) on Ru and Co surfaces [58] are the microscopic origin of the critical occupations $x_{k,o}$. More negative formation energies indicate stronger stabilization of the promoter gear on the surface, leading to higher equilibrium occupations $x_{k,o}$. From the definition (Eq. (25))

$$C_k=-{\displaystyle \frac{\left(k-19\right)x_{k,o}^2}{1-2x_{k,o}}}$$

this directly produces larger $\left|C_k\right|$ (Theorem 16 and Theorem 18). Ba and La promoters show the most negative values on both Ru and Co, consistent with their superior activity observed in the main experimental table. Table 4 therefore provides direct support for Theorems 12, 15 and 18: catalyst preparation is a variational selection of interface constants that controls the grand-potential shift.

11.2. N adsorption energy $E_N$ and $N-N$ transition state energies $E_{TS}$

N adsorption energy $E_N$ is a classic descriptor in ammonia synthesis. In ZMT it is a closed function of the chemical-potential shift induced by $\left\{C_k\right\}$. The promoter-induced lowering of the $N-N$ transition state energy $E_{TS}$ is the kinetic manifestation of the increased variational driving force. Promoters systematically reduce $E_{TS}$ (especially on Co and Ru), moving the system to more favorable positions on the activity volcano (Theorem 13). Table 5 demonstrates linear scaling of adsorption and activation energies with the grand-potential shift (Theorem 12) and confirms that the rate enhancement arises directly from the $C_k$-tuned driving force.

11.3 Electric field $\epsilon$, $E_{promotion}$, and ${\Delta E}_{TS}$ on promoter-doped Ru and Co surfaces

The electric field generated by promoters acts as a local conjugate variable (surface potential/charge pair). $E_{promotion}$ and the change in transition-state energy ${\mathsf{\Delta }E}_{TS}$ quantify how strongly the promoter shifts the grand potential, directly correlating with $\left|C_k\right|$ and the smallest eigenvalue ${\lambda }_{min}$. Stronger promoters (Ba, La) produce larger $\left|\epsilon \right|$ and more negative ${\mathsf{\Delta }E}_{TS}$, exactly as predicted by the multi-conjugate generalization (Theorem 17). These tables extend the framework to the electrostatic picture, showing that the same variational mechanism governs both gas-phase and electrochemical environments.

Taken together, tables (4-7) provide strong experimental/DFT support for the entire ZMT chain: promoter chemistry tunes surface gear occupations $x_{k,o}$ (Table 2) → determines $C_k$ (Theorem 15) → shifts the grand potential → modifies adsorption energies and barriers (Table 5, Table 6 and Table 7) → controls macroscopic activity and selectivity (Theorems 12, 13, and 16). This dataset validates that what was previously empirical is now deductively predictable from the interface constants alone.

11.4 Temporary versus Permanent Poisoning

Poisoning is a shift of the solid blackbox away from its optimal $\left\{C_k\right\}$ configuration.

Temporary (reversible) poisoning occurs when a poison species adds a new period prime or modifies $x_{k,o}$ in a reversible way. The Lyapunov exponents ${\lambda }_{k,19}$ remain non-negative and the system can return to the optimal point once the poison is removed. The grand-potential shift $\mathsf{\Delta }\omega$ (Eq. (22)) is finite and the equilibrium manifold simply moves to a lower-conversion point on the same vapor–solid locus.

Permanent (irreversible) poisoning occurs when the poison permanently alters the gear union or drives one or more ${\lambda }_{k,19}<0$. This violates marginal stability and collapses the spectral gap. The catalyst can no longer reach the hard upper limit even after the poison is removed. Regeneration would require re-tuning the blackbox (changing the effective primes or $x_{k,o}$).

Both cases are fully deductive: poisoning is simply a change in the effective $\left\{C_k\right\}$ vector, and its reversibility is determined by whether the Lyapunov spectrum remains non-negative.

11.5 Experimental Validation: Catalyst Activity at Constant Pressure and Zero Inerts

A comprehensive overview of recent advances in ammonia synthesis catalysts is provided by the meta-analysis of Cao et al. (2022) [58], which compiles NH₃ synthesis rates versus temperature for a wide family of materials (primarily Ru-based catalysts with various promoters and supports) at 10 bar and low/no inerts (H₂/N₂ = 3).

A direct experimental validation of the ZMT framework is provided by the extensive dataset of ammonia synthesis catalysts tested at fixed pressure of 10 bar with negligible inerts (fresh-feed conditions). Under these constraints, bulk thermodynamic variables are essentially identical, isolating the effect of the catalyst.

The reported activities span more than one order of magnitude (from ~1.7 to 88.1 mmol/g·h). In the classical paradigm this large variation would require separate empirical rate constants and fitted adsorption terms for each formulation. In ZMT the entire spread is explained deductively by differences in the interface constants $\left\{C_k\right\}$.

Highest activities are achieved by catalysts that generate large positive grand-potential shifts $\sum C_k$:

• Ba-Co/C (86.4 mmol/g·h at 440 °C)
• Ru/Ba/LaCeOₓ (88.1 mmol/g·h at 400 °C)
• Ba-Ru-Li/AC (46.3 mmol/g·h at 459 °C)

Barium and lithium promoters are particularly effective because they stabilize favorable critical occupations $x_{k,o}$ on the solid surface, maximizing $\mid C_k\mid$. Ruthenium-based systems with these promoters consistently outperform unpromoted or oxide-supported catalysts.

Moderate performers (Cs- or K-promoted Ru, certain Co hydrides) correspond to intermediate $\mid C_k\mid$ values. Low-activity catalysts (e.g., Ru/Ce₀.₆Zr₀.₄O₂ at 1.7 mmol/g·h, various Ni and poorly supported Co systems) exhibit significantly smaller $\mid C_k\mid$, reflecting suboptimal prime activation and surface occupations.

Temperature modulates the helical partition functions $Z_i\left(s\right)$, producing the expected increase in activity with rising temperature. However, the ranking among catalysts remains largely preserved, confirming that $\left\{C_k\right\}$ is the dominant control variable once pressure and inerts are fixed.

This dataset demonstrates that, at constant pressure and zero inerts, catalytic performance is governed almost exclusively by the variational knob $\left\{C_k\right\}$, which is tuned through promoter choice, support interaction, and metal selection. The ZMT framework thereby converts what was previously empirical catalyst screening into a predictive search over achievable surface occupations $x_{k,o}$ on the prime basis — without any fitted parameters.

The results (Table 8) provide strong experimental confirmation that the single vector $\left\{C_k\right\}$, derived deductively from the solid-phase grand-partition function $Z_{\mathrm{solid}}\left(s\right)$, controls macroscopic activity as predicted by Theorems 12, 16 and 18.

12. Comparative Analysis: Conventional Haber–Bosch Loop

The ammonia synthesis reaction was modeled using kinetic parameters reported in the literature for an industrial Fe-Wüstite-based catalyst [2,3,52,53,54,55,56,57,58,59,60,61]. Because industrial converters typically operate with sufficient residence time to approach thermodynamic equilibrium, the modeling framework focused primarily on thermodynamic optimization while using industrially validated residence times to account for kinetic limitations.

12.1 Distinct Operating Regimes of Fe versus Ru Catalysts

With the catalyst treated as a fixed blackbox (i.e., fixed period prime and associated constants $\left\{C_k\right\}$), the ZMT framework yields an explicit equilibrium manifold in the $\left(T,P\right)$ plane. This section presents the numerical construction and the resulting equilibrium conversion for two representative catalysts [2,3,13,52,53,54,55,56,57,58,59,60,61,62,63]— iron and ruthenium — using only closed-form expressions derived from the deductive theorems.

12.1.1 Calculation of Blackbox Constants $C_k$

For a given catalyst, the integration constants are obtained from the closed-form expressions (Eq. (14), Eq. (25)):

$$x_{k_o}={\displaystyle \frac{1/\mathrm{l}\mathrm{n}{p}_k}{1/\mathrm{l}\mathrm{n}{p}_k+1/\mathrm{l}\mathrm{n}19}},C_k=-{\displaystyle \frac{\left(p_k-19\right){\left(x_{k_o}\right)}^2}{1-2x_{k_o}}}$$

For iron ($p=11$): $x_{k_o}\approx 0.6021$, hence $C_k\approx -21.31$.

For ruthenium ($p=13$): $x_{k_o}\approx 0.5512$, hence $C_k\approx -23.75$.

These values of $C_k$ represent the only tunable parameters of the solid blackbox for each catalyst.

12.1.2 Grand-Partition Functions and Fugacity-Based Equilibrium Constant

The grand-partition function for each gaseous species is given by:

$$Z_i(s)={\displaystyle \frac{1}{1-p_i^{-s}}}$$

with $p_{{\mathrm{N}}_2}=5$, $p_{{\mathrm{H}}_2}=2$, and $p_{\mathrm{N}{\mathrm{H}}_3}=7$ [23]. The fugacity-based equilibrium constant, modified by the solid blackbox through the interface matching condition, takes the form:

$$K(P,s,\{C_k\})={\displaystyle \frac{{\left(P_{\mathrm{N}{\mathrm{H}}_3}{Z}_{\mathrm{N}{\mathrm{H}}_3}(s)\right)}^2}{{\left(P_{{\mathrm{N}}_2}{Z}_{{\mathrm{N}}_2}(s))(P_{{\mathrm{H}}_2}{Z}_{{\mathrm{H}}_2}(s)\right)}^3}}\mathrm{e}\mathrm{x}\mathrm{p}\left({\displaystyle \frac{\sum {\nu }_k{C}_k}{RT}}\right)$$

The exponential term represents the direct thermodynamic contribution of the solid catalyst.

12.1.3 Distinct Operating Regimes

ODE for stoichiometric feed with no inerts produces the following equilibrium manifolds (Table 9 and Table 10):

The results reveal distinctly different operating regimes. The iron catalyst exhibits its highest equilibrium conversion in the traditional high-pressure region (typically 10–15 MPa), consistent with historical industrial practice. In contrast, the ruthenium catalyst achieves comparable or higher conversion at significantly lower pressures (2–6 MPa), demonstrating the theoretical advantage of more active catalysts in enabling milder operating conditions.

12.1.4 Promoter and Dopant Effects: Tuning the Solid Blackbox

Within the ZMT framework, promoters and structural dopants do not alter the period prime of the active metal. Instead, they modify the effective critical composition $x_{k_o}$ of the catalyst blackbox through electronic and structural interactions. This adjustment shifts the integration constants $C_k$ without introducing additional free parameters, thereby providing a tunable lever to reposition the global minimum of the phase functional $F$.

For pure iron ($p_k=11$), the baseline values are $x_{k_o}\approx 0.6021$ and $C_k\approx -21.31$. Promoter addition (e.g., potassium) or support modification (e.g., MgO) [13,58,59,60,61,62,63,64] effectively increases $\mid C_k\mid$, improving equilibrium conversion at lower pressures. This mechanism offers a deductive explanation for the well-known promotional effects observed experimentally in industrial ammonia catalysts.

12.1.5 Dopant-Induced Shifts in Effective Blackbox Constants

Within the ZMT framework, promoters and structural dopants do not alter the period prime of the active metal. Instead, they modify the effective critical composition $x_{k_o}$ of the catalyst blackbox, thereby shifting the integration constants $C_k$ without introducing additional free parameters.

For pure iron ($p_k=11$), the baseline values are:

$$x_{k_o}\approx 0.6021,C_k\approx -21.31$$

Potassium functions as an electronic promoter by donating electrons to the iron surface. This effect increases the weight of the catalyst prime in the critical-composition formula, producing a positive shift $\mathsf{\Delta }{x}_{k_o}$. At typical industrial potassium loadings (0.5–2 wt% K₂O, corresponding to approximately 0.4–1.6 mol% K), the shift lies in the range $\mathsf{\Delta }{x}_{k_o}\approx 0.01$– 0.03. The resulting effective blackbox constants are:

• 0.5 mol% K: $x_{k_o}\approx 0.6121\to C_k\approx -22.05$ (more negative by ≈0.74)

• 1.0 mol% K: $x_{k_o}\approx 0.6221\to C_k\approx -22.85$ (more negative by ≈1.54)

• 1.6 mol% K: $x_{k_o}\approx 0.6321\to C_k\approx -23.70$ (more negative by ≈2.39)

Magnesium oxide (MgO) acts primarily as a structural promoter. It stabilizes the high surface area of the iron catalyst and prevents sintering, thereby preserving the potassium-induced shift in $x_{k_o}$ throughout the catalyst lifetime. While MgO does not strongly alter $x_{k_o}$ by itself, it significantly enhances the long-term stability of the tuned $C_k$ values.

For a typical industrial Fe-K-MgO catalyst (0.5–1.5 mol% K + 2–5 wt% MgO) [13,58,59,60,61], the combined effect produces an effective blackbox constant in the range:

$$C_{k,\mathrm{eff}}\approx -22.8 \mathrm{to} -23.7$$

This value is substantially more negative than the pure-iron baseline ($C_k\approx -21.31$).

12.2 ZMT Perspective on Recent Low-Pressure, Zero-Inerts Studies on Promoted Fe Catalysts: The Critical Role of Inerts and the Lab–Industry Gap

Recent experimental studies on promoted Fe/MgO catalysts (Ba- and K-doped) [65,66] have demonstrated impressive NH₃ synthesis rates and yields at remarkably low pressure (1.0 MPa / 10 bar) and zero inerts. The ZMT variational framework provides a unified explanation for both the laboratory success under ideal conditions and the persistent challenges observed in industrial multi-bed converters. It does so through: the ODE (Section 9.8) governing static equilibrium and governing dynamical relaxation rates ${\lambda }_k$.

12.2.1 Experimental Conditions and Key Results

Both studies evaluate promoted Fe/MgO catalysts at 10 bar using pure stoichiometric N₂ + 3H₂ feeds (zero inerts) [65,66]. Key observations include:

• Temperature dependence: Ba- and K-promoted Fe/MgO catalysts achieve rates of approximately 23–25 mmol h⁻¹ g_cat⁻¹ at 400 °C, with NH₃ yields around 1.5 %. These materials significantly outperform both commercial fused-iron catalysts and undoped Fe/MgO.
• Pressure dependence: At fixed 350 °C and zero inerts, the K-promoted Fe/K(3)/MgO-500red catalyst exhibits the strongest positive pressure dependence, reaching ~30 mmol h⁻¹ g_cat⁻¹ at 3 MPa. While Ru-based benchmarks remain competitive, unpromoted Fe/MgO shows nearly flat pressure response.

These results might suggest that low-pressure ammonia synthesis becomes straightforward once inerts are removed. However, this interpretation overlooks the industrial reality that the ZMT framework captures quantitatively.

12.2.2 The Lab–Industry Gap

The recent laboratory studies optimize the most favorable regime ($\gamma =0$) and therefore realize the maximum thermodynamic and kinetic benefit of the engineered blackbox. They demonstrate that promoters (Ba or K) tune $C_k$ and the Hessian curvature so that temperature becomes the dominant variable and pressure sensitivity is restored when inerts are absent.

Industrial operation, however, forces the catalyst to function where rising $\gamma$ imposes a double penalty—thermodynamic (lower ${\xi }_{\mathrm{e}\mathrm{q}}$) and kinetic (lower ${\lambda }_{\mathrm{m}\mathrm{i}\mathrm{n}}$). The ZMT framework quantitatively predicts the well-known third-bed plateau: lowest driving force, slowest rate, flattest temperature profile, and smallest marginal conversion despite the largest catalyst volume.

The laboratory results are not incorrect; they simply probe the catalyst under ideal conditions. The true industrial test occurs as inerts accumulate. This is precisely the regime in which the combined effect of increasing $\gamma$ is most severe, which is why the third bed consumes the most catalyst for the least gain.

12.2.3 Quantitative Support from the Hessian Table

Table 11 quantifies the kinetic penalty, evaluated at equilibrium ${\xi }_{\mathrm{e}\mathrm{q}}(P,T,\gamma )$ with the industrial blackbox value $C_k=-24.3$:

Higher inert levels systematically reduce $\mid {\lambda }_{\mathrm{m}\mathrm{i}\mathrm{n}}\mid$, reproducing the kinetic limitation observed in the industrial third bed.

The recent low-pressure, zero-inert studies vividly illustrate the capability of a well-engineered blackbox when inerts are eliminated: promoted Fe/MgO catalysts deliver high rates and yields at only 10 bar through simultaneously elevated ${\xi }_{\mathrm{e}\mathrm{q}}$ and elevated ${\lambda }_{\mathrm{m}\mathrm{i}\mathrm{n}}$. The ZMT framework accounts for this behavior as the direct consequence of the variational principle on the helical manifold $G=-RT\mathrm{l}\mathrm{n}Z(s)$.

At the same time, these studies underscore why laboratory results obtained under pure conditions do not fully translate to industrial practice. In a real multi-bed converter the progressive rise in effective $\gamma$ imposes a compounded thermodynamic–kinetic penalty that is most acute in the final bed. The ZMT framework captures this entire picture from first principles and thereby supplies a deductive link between optimistic laboratory performance and the practical constraints of industrial operation.

This constitutes one of the strongest validations of the ZMT approach to date: it explains both the laboratory success and why the decisive test of a catalyst is its behavior as inerts accumulate across the converter. The data further motivate the development of a mildly condition-dependent blackbox $C_k(\gamma ,T)$, a natural and immediate extension of the existing framework.

12.3 Recent Advances in Low-Pressure Ammonia Synthesis with Inverse-Structure Iron Catalysts

A particularly instructive example of modern catalyst engineering [67,68,69] is the inverse-structure iron catalyst reported by Hattori et al. (2025) [67]. In conventional supported catalysts, metal nanoparticles are dispersed on a high-surface-area oxide support. In contrast, the new AlH-K⁺/Fe catalyst inverts this architecture: dense metallic iron particles are decorated with aluminum hydride (AlH) species and potassium. The AlH functions as a strong electron donor directly on the iron surface rather than as a traditional support, creating a high density of active sites with exceptionally strong electron-donating capability.

Under low-pressure conditions (0.1 MPa) at 400 °C, the AlH-K⁺/Fe catalyst achieves more than double the ammonia synthesis rate per catalyst volume compared with the conventional industrial promoted-iron catalyst, despite possessing only half the specific surface area. The apparent activation energy is markedly reduced to 37 kJ mol⁻¹ (versus ~44 kJ mol⁻¹ for industrial promoted-Fe and >80 kJ mol⁻¹ for unpromoted iron). The catalyst also exhibits excellent long-term stability, maintaining high activity for more than 100 h at 300 °C and 0.9 MPa with no detectable structural degradation.

These results align directly with the ZMT variational framework. The inverse structure and potassium promotion effectively tune the catalytic blackbox $\left\{C_k\right\}$ to a more negative value in the low-inert regime, simultaneously increasing the thermodynamic equilibrium conversion ${\xi }_{\mathrm{e}\mathrm{q}}$ and the relaxation rate ${\lambda }_{\mathrm{m}\mathrm{i}\mathrm{n}}$ (via the Hessian of $\mathrm{l}\mathrm{n}Z(s)$). The strong pressure dependence and low activation energy observed are consistent with the helical grand partition functions $Z_i(s)$ and the stoichiometric $1/P^2$ term becoming dominant once the inert mole fraction $\gamma$ is minimized.

This work highlights the practical implication of the ZMT analysis: modern promoters can engineer a blackbox that performs exceptionally well under low-pressure, low-inert conditions. In contrast, industrial multi-bed converters operate with progressively rising effective $\gamma$ across beds (arising from the net mole reduction of the reaction). The resulting decrease in both $X_{\mathrm{e}\mathrm{q}}$ and ${\lambda }_{\mathrm{m}\mathrm{i}\mathrm{n}}$ (Hessian sensitivity of 20–24 % reduction in $\mid H\mid$ when $\gamma$ increases from 1 % to 11 %) explains the characteristic third-bed plateau — the largest catalyst volume yet the smallest incremental conversion and flattest temperature profile.

The Hattori et al. (2025) [67] inverse-structure catalyst therefore provides converging experimental support for the ZMT prediction that the performance ceiling is set primarily by the blackbox strength and inert mole fraction rather than by total pressure alone. It also underscores the opportunity for next-generation reactor designs that maintain low effective $\gamma$ (e.g., via inerts management or membrane integration) to fully exploit such optimized blackboxes at lower operating pressures.

13. Conclusions

The Zeta-Minimizer Theorem provides a complete deductive framework for heterogeneous catalysis and chemical kinetics. Starting from three primitive axioms and the helical geometry of the grand-partition function $Z\left(s\right)$, every observable — equilibrium conversion, rate law, linear scaling relations, Sabatier volcano, selectivity, and catalyst performance — emerges without empirical parameters or fudge factors.

For ammonia synthesis, the variational rate law reproduces the classical Temkin–Pyzhev expression as a limiting case in the high-inerts regime while eliminating its singularities at zero ammonia and extending naturally to low-inerts and low-pressure conditions. The implicit analytical solution for transient conversion, the modal decomposition for multi-reaction coupling, and the explicit sensitivities to inerts, pressure, and temperature are all closed functions of the phase functional and the catalyst interface constants $\left\{C_k\right\}$. These constants themselves are deductively determined by the solid-phase grand-partition function $Z_{\mathrm{solid}}\left(s\right)$ and the surface occupations $x_{k,o}$, turning catalyst preparation (reduction temperature, promoters, support) into a variational search over the prime basis.

The framework is representation-independent. The same $Z\left(s\right)$ and variational rate law govern any conjugate pair — pressure–volume, electrochemical potential–charge, magnetic field–magnetization — and extend seamlessly to multi-conjugate systems. In multi-reaction environments the symmetric reaction Hessian diagonalizes into independent normal modes, revealing the dominant pathways and relaxation spectrum without manual analysis of rate-limiting steps.

What was previously hidden behind empirical fitting and ad-hoc corrections is now exposed as geometry on a single grand-potential landscape. The classical model worked well only within its narrow fitting window; ZMT works across the full operating envelope, explains why the classical model succeeded and where it failed, and supplies a predictive, first-principles route to next-generation catalysts and reactor design.

The Zeta-Minimizer Theorem therefore marks the transition from model-based tuning to deductive design — from a collection of fitted parameters to a unified variational principle that governs both the microscopic interface and the macroscopic performance of industrial processes.