Dynamic Balance: A Thermodynamic Principle for the Emergence of the Golden Ratio in Open Non-Equilibrium Steady States

Abstract

We develop a symmetry-based variational theory that shows the coarse-grained balance of work inflow to heat outflow in a driven, dissipative system relaxed to the golden ratio. Two order-2 Möbius transformations—a self-dual flip and a self-similar shift—generate a discrete non-abelian subgroup of $\mathrm{PGL}(2,\mathbb{Q}\left(\sqrt{5}\right))$. Requiring any smooth, strictly convex Lyapunov functional to be invariant under both maps enforces a single non-equilibrium fixed point: the golden mean. We confirm this result by (i) a gradient-flow partial-differential equation, (ii) a birth–death Markov chain whose continuum limit is Fokker–Planck, (iii) a Martin–Siggia–Rose field theory, and (iv) exact Ward identities that protect the fixed point against noise. Microscopic kinetics merely set the approach rate; three parameter-free invariants emerge: a $62\%:38\%$ split between entropy production and useful power, an RG-invariant diffusion coefficient linking relaxation time and correlation length ${\mathcal{D}}_{\alpha }={\xi }^z/\tau$, and a $\vartheta ={45}^{\circ }$ eigen-angle that maps to the golden logarithmic spiral. The same dual symmetry underlies scaling laws in rotating turbulence, plant phyllotaxis, cortical avalanches, quantum critical metals, and even de-Sitter cosmology, providing a falsifiable, unifying principle for pattern formation far from equilibrium.

1. Introduction

The golden ratio, $\phi \approx 1.618$, has been documented in phyllotactic leaf arrangements; branching patterns of trees, blood vessels, lightning, and river deltas [1,2,3,4]; logarithmic spirals in hurricanes and galactic arms [5,6]; power law exponents in rotating turbulence [7,8,9,10]; quasicrystals [11,12,13]; mass gaps of critical Ising chains [14,15]; band structures of twisted bilayer graphene [16,17,18]; Fibonacci anyons dimensionality [19]; avalanche statistics in cortical activity [20,21,22]; and more. Classical equilibrium thermodynamics, tied to isolated systems and static-state variables, cannot explain the ubiquity of spirals, branching, fractals or, scale-invariant kinetics across scales [23,24,25], and a unifying, model-independent, non-equilibrium principle explaining the golden mean universality class is yet to emerge [26,27,28,29,30]. Empirically, these systems share two attributes: (i) they are maintained far from equilibrium by continuous energy or matter influx (e.g., solar radiation, gravitational shear, biochemical energy) [25], and (ii) they exhibit nonlinear, irreversible dissipation (e.g., thermal conduction, radiative cooling, viscous dissipation, chemical enthalpy release) while maintaining coherent large-scale organization and optimal functionality [31].

In this work, we show that a single symmetry-protected variational principle—Dynamic Balance (DB)—forces any driven–dissipative system to relax towards the golden ratio. We define a dimensionless ratio–entropy flux field $\alpha \left(t\right)$–comparing two irreducible forms of energy: useful power inflow $\dot{E}$ (reversible) and entropic heat outflow $T\dot{S}$ (irreversible). A large $\alpha$ implies low dissipation; a small $\alpha$ means the system overheats, leaving little energy to build or maintain structure. Most real systems develop negative feedback loops (physiological, hydrodynamic, or electronic) to prevent collapse or runaway behavior, thereby stabilizing their internal state [32]. For instance, excessive anabolism or catabolism harms living organisms, so metabolic circuits self-regulate through hormones and growth-factor inhibition [33]. Consequently, a driven system in a sustained non-equilibrium steady state will adjust $\alpha \left(t\right)$ toward a constant, ${\alpha }^{★}$, that optimally balances useful work against dissipation. Across scales, systems at ${\alpha }^{★}$ exhibit hierarchical organization—smaller subsystems nested within larger ones—and scale-invariant, fractal dynamics. Therefore, we posit that this optimum is the golden ratio, ${\alpha }^{★}=\phi$.

From a group-theoretical perspective, every physical process is the result of a selection rule dictated by an underlying symmetry [34,35]. We show that two discrete Möbius maps $S_{\phi }$ (a self-dual involution) and $T_{\phi }$ (a self-similar recursive shift) acting on $\alpha$, and realized microscopically by antisymmetric Onsager couplings [36,37,38] and cross-correlated noise in active media [39], respectively, generate a non-Abelian subgroup of $\mathrm{PGL}\left(2,\mathbb{Q}\left(\sqrt{5}\right)\right)$. Requiring any strictly convex Lyapunov cost $R\left(\alpha \right)$ to remain invariant under both maps singles out the unique stable fixed point ${\alpha }^{★}=\phi$. Microscopic details affect only transients, and the system dynamically self-organizes around the golden mean. This modular symmetry and convex geometry guide the entropy flux field, giving rise to Legendre dual flows and emergent conjugate pairs observed in all physical laws.

$\mathbf{Symmetry}\stackrel{\mathrm{Modular}\mathrm{Duality}}{\to }\mathbf{Geometry}\stackrel{\mathrm{Convexity}}{\to }\mathbf{Dynamics}\stackrel{\mathrm{Legendre}\mathrm{Duality}}{\to }\mathbf{Observables}$

The paper is organized as follows: Section 2 states and proves the core theorem. Section 3 builds the gradient-flow PDE and extracts three parameter-free invariants (entropy split, ${\xi }^2\mathsf{\Gamma }$, and ${45}^{\circ }$ spiral pitch). Section 4 shows that a discrete birth–death Markov model recovers the same continuum PDE limit, while Section 5 embeds the theory in a Martin–Siggia–Rose field integral and derives exact Ward identities that protect $\phi$ against noise. We close with cross-domain applications—turbulence, neuroscience, quantum criticality, and cosmology.

2. Mathematical Framework

Theorem 1. 

Let $S_{\phi }:\alpha \mapsto {\mathsf{\Lambda }}^2/\alpha$ and $T_{\phi }:\alpha \mapsto 1+1/\alpha$ be two Möbius transformations generated microscopically by

$S_{\phi }:$ antisymmetric Onsager exchange $L_{AB}=-{\mathsf{\Lambda }}^2{L}_{BA}$

$T_{\phi }:$ cross–correlated noise source source $\langle {\xi }_A{\xi }_B\rangle$

These two maps generate a discrete, non-Abelian subgroup of $\langle S_{\phi },T_{\phi }\rangle \subset \mathrm{PGL}\left(2,\mathbb{Q}\left(\sqrt{5}\right)\right)$. Define the minimum smooth and strictly convex Lyapunov functional $\mathcal{C}\left(\alpha \right)\in (0,\infty )$ such that $\mathcal{C}\circ S_{\phi }=\mathcal{C}$, with ${\mathcal{C}}^{″}>0$ and $\mathcal{C}\to \infty$ as $\alpha \to 0^{+},\infty$. Then,

(a) 

Any minimizer ${\alpha }^{★}$ satisfies $S_{\phi }\left({\alpha }^{★}\right)=T_{\phi }\left({\alpha }^{★}\right)={\alpha }^{★}$.

(b) 

Combining the two fixed-point equations gives ${\alpha }^{★2}-{\alpha }^{★}-1=0$ and ${\alpha }^{★2}={\mathsf{\Lambda }}^2$.

(c) 

Therefore, $\begin{array}{|c|}\hline \mathsf{\Lambda }={\alpha }^{★}=\phi \\ \hline\end{array}$.

In the corresponding Martin–Siggia–Rose field theory, the Ward identities

$\langle \widehat{\alpha }(\frac{{\mathsf{\Lambda }}^2}{\alpha }-\alpha )\rangle =\langle \widehat{\alpha }(1+\frac{1}{\alpha }-\alpha )\rangle =0$

forbid stochastic drift. Thus, once Λ is fixed by symmetry, the golden ratio φ becomes a conformal noise-protected attractor: a unique entropic fixed point stabilized by both modular symmetry and convex geometry.

2.1. Coarse-Grained Energy and Entropy Fluxes

Consider an open, driven-dissipative system held in a non-equilibrium steady-state (NESS) by a continuous influx of power and outflux of heat. Let

(1)$A\left(t\right)=J_{rev}=\dot{E}\left(t\right),B\left(t\right)=J_{irrev}=T\left(t\right)\dot{S}\left(t\right)$

denote, respectively, the instantaneous work (reversible) flux and the entropic heat (irreversible) flux in the usual system-oriented sign convention ($A,B>0$). Here, T is an effective temperature characterizing internal fluctuations, and $\dot{S}$ is the entropy-production rate. Both A and B are assumed $C^1$ functions on $[0,\infty )$, and the total throughput is

(2)$P\left(t\right)=A\left(t\right)+B\left(t\right)>0.$

In a mesoscopic description, these two irreducible fluxes may originate from different blocks of the Onsager matrix or from distinct fields in a Martin–Siggia–Rose path integral, coupled solely through antisymmetric (reactive) exchange and cross-correlated noise [31,38].

2.2. Energy–Entropy Flux Ratio

We define a dimensionless ratio describing entropy flow per unit of available energy (or information):

(3)$\alpha \left(t\right)={\displaystyle \frac{A\left(t\right)}{B\left(t\right)}}={\displaystyle \frac{\dot{E}\left(t\right)}{T\left(t\right)\dot{S}\left(t\right)}},\alpha \left(t\right)>0,\alpha \in C^1\left([0,\infty )\right).$

Neither extreme limit of $\alpha$ is sustainable in a steady-state: $\alpha \to 0^{+}$ corresponds to total dissipation (black hole collapse), while $\alpha \to \infty$ implies vanishing entropy export and thermal runaway (wormhole divergence). Real driven systems, therefore, self-tune to an interior fixed value ${\alpha }^{★}$, which is the focus of the analysis that follows.

2.3. Modular Symmetry and Convex Geometry

We define the fundamental symmetry of the entropy field $\alpha (\mathbf{x},t)$ via a modular action that preserves the golden ratio as its unique fixed point. This symmetry is generated by the transformations:

(4)$T_{\phi }:\alpha \mapsto 1+{\displaystyle \frac{1}{\alpha }},S_{\phi }:\alpha \mapsto {\displaystyle \frac{{\phi }^2}{\alpha }},$

which together generate a closed, non-Abelian, and discrete subgroup of the projective modular group: $〈S_{\phi },T_{\phi }〉\subset \mathrm{PGL}\left(2,\mathbb{Q}\left(\sqrt{5}\right)\right)\subset \mathrm{PGL}(2,{\mathbb{R}}^{+})$, whose transformation acts on the entropy field domain $\alpha$ as:

(5)$\mathrm{PGL}(2,\mathbb{Q}\left(\sqrt{5}\right))=\left\{\left.\alpha \mapsto {\displaystyle \frac{a\alpha +b}{c\alpha +d}}\right|a,b,c,d\in \mathbb{Q}\left(\sqrt{5}\right),ad-bc\ne 0\right\}$

and satisfies the presentation:

(6)$〈S_{\phi },T_{\phi }|S_{\phi }^2={\left(S_{\phi }{T}_{\phi }\right)}^3=\mathcal{I}〉.$

mirroring the modular group relations, but acting on a real positive entropy domain with golden-arithmetic structure. The unique fixed point of this group is the golden ratio ${\alpha }^{★}=\phi ,$ which plays the role of a conformal attractor in entropy space. The group $\mathrm{PGL}\left(2,\mathbb{Q}\left(\sqrt{5}\right)\right)$ is richer and more physical than $\mathrm{PSL}(2,\mathbb{Z})$, which is discrete, integer, and conformal only in $\mathbb{H}$, and it governs all entropy flows.

We define the unique minimal convex functional invariant under this symmetry:

(7)$\mathcal{C}\left(\alpha \right)=\mathcal{C}\left({\displaystyle \frac{{\phi }^2}{\alpha }}\right)={\displaystyle \frac{1}{2}}\left(\alpha +{\displaystyle \frac{{\phi }^2}{\alpha }}\right),$

which attains its minimum precisely at ${\alpha }^{★}=\phi$. This symmetric functional $\mathcal{C}\left(\alpha \right)$ plays a Casimir-like role for the modular dynamic system, structurally defining the potential that drives entropy flow. Level sets of $\mathcal{C}\left(\alpha \right)$ form equipotential surfaces. Gradient descent of $\mathcal{C}\left(\alpha \right)$ defines entropy flow trajectories, and its derivatives ${\mathcal{C}}^{\prime }\left(\alpha \right)$ (gradient), ${\mathcal{C}}^{″}\left(\alpha \right)$ (curvature), ${\mathcal{C}}^{‴}\left(\alpha \right)$ (torsion), …, fully governs the geometric and dynamical structure of the entropy field $\alpha (\mathbf{x},t)$. It defines the gradient flow driving systems toward dynamic balance, the curvature tensor determining local entropy rigidity, and the modular invariance ensuring global recursive symmetry.

The conformal invariance at ${\alpha }^{★}=\phi$ implies local scale symmetry, which protects the system from perturbative deformations and endows it with intrinsic resilience to noise and long-range entropy coherence preservation.

2.4. Microscopic Origin of the Möbius Involution $S_{\mathsf{\Lambda }}$: Onsager Antisymmetric Reactive Exchange

In linear response theory, conjugate pairs of fluxes $J_i$ and thermodynamic forces $X_j$ obey

(8)$J_i=-\sum _{j\in \{A,B\}}{L}_{ij}{X}_j,$

where $L_{ij}$ is the Onsager matrix [36]. We focus on the purely off-diagonal, entropy-free coupling block: $L_{AA}=L_{BB}=0,L_{AB}\ne 0,L_{BA}\ne 0,$ and impose antisymmetric reciprocity $L_{AB}=-L_{BA}$, as typical for conservative or reactive couplings (e.g., in Hall transport or chemical oscillators [38,40]. Denoting the output power fluxes by $A\equiv J_A$, $B\equiv J_B$, we have:

(9)$\left(\begin{array}{c}A\\ B\end{array}\right)=-\left(\begin{array}{cc}0& L_{AB}\\ -L_{BA}& 0\end{array}\right)\left(\begin{array}{c}{X}_A\\ X_B\end{array}\right).$

Hence, the entropy flux ratio $\alpha =A/B=-L_{AB}{X}_B/L_{BA}{X}_A$ transforms under the exchange of channels $(A,X_A)\leftrightarrow (B,X_B)$ as

(10)${\alpha }^{\prime }={\displaystyle \frac{B}{A}}={\displaystyle \frac{|L_{AB}|}{|L_{BA}|}}{\displaystyle \frac{1}{\alpha }}={\displaystyle \frac{{\mathsf{\Lambda }}^2}{\alpha }},{\mathsf{\Lambda }}^2\equiv \left|{\displaystyle \frac{L_{AB}}{L_{BA}}}\right|>0.$

This is a Möbius transformation of order 2:

(11)$S_{\mathsf{\Lambda }}:\alpha \mapsto {\displaystyle \frac{{\mathsf{\Lambda }}^2}{\alpha }},S_{\mathsf{\Lambda }}=\left(\begin{array}{cc}0& {\mathsf{\Lambda }}^2\\ 1& 0\end{array}\right),S_{\mathsf{\Lambda }}^2=\mathrm{id},$

The constant $\mathsf{\Lambda }$ quantifies the microscopic asymmetry between the reactive couplings. Its value will later be fixed by requiring modular self-similarity of the dynamics under golden-ratio recurrence.

2.5. Microscopic Origin of the Self-Similar Shift $T_{\phi }$: Cross-Correlated Noise

After establishing that the Möbius flip $S_{\phi }$ arises from antisymmetric reactive exchange, we now identify the origin of the self-similar modular shift $T_{\phi }$ in terms of correlated stochastic forcing. In many mesoscopic and active systems, the random forces driving work-like (A) and entropy-like (B) channels are not statistically independent, but exhibit finite cross-correlation $C=〈{\xi }_A{\xi }_B〉$ [38,39]. We model this via a coupled Ornstein–Uhlenbeck process:

(12)$\dot{A}=-{\mathsf{\Gamma }}_A(A-A_0)+kB+{\xi }_A\left(t\right),\dot{B}=-{\mathsf{\Gamma }}_B(B-B_0)-kA+{\xi }_B\left(t\right),$

with Gaussian white noise correlations,

(13)$〈{\xi }_i\left(t\right){\xi }_j\left(t^{\prime }\right)〉=2\left(D_i{\delta }_{ij}+C(1-{\delta }_{ij})\right)\delta (t-t^{\prime }),i,j\in \{A,B\},$

where the diagonal elements $D_A$ and $D_B$ set individual noise intensities (the variance), and $\left|C\right|\le \sqrt{D_A{D}_B}$ quantifies the cross-correlation. The coupling coefficient $k=L_{AB}\left|\partial X_B/\partial B\right|$ arises from the same Onsager-antisymmetric exchange responsible for $S_{\phi }$. Once the physical units of A and B are rescaled to be commensurate (both interpreted as power fluxes), the conversion factor is absorbed into k. The new ingredient is the non-diagonal diffusivity $D_{AB}=C$, encoding the noise-level correlation between the two channels. Solving the Lyapunov equation for the stationary covariance (see SI) yields the steady-state flux ratio:

(14)$\overline{\alpha }\equiv {\displaystyle \frac{\langle A\rangle }{\langle B\rangle }}={\displaystyle \frac{{\mathsf{\Gamma }}_B{D}_A-kC}{{\mathsf{\Gamma }}_A{D}_B-kC}}$

When $C=0$, the antisymmetric dynamics reproduce the Möbius flip $\overline{\alpha }\mapsto 1/\overline{\alpha }$, matching the action of $S_{\phi }$ on the mean state. A non-zero cross-correlation modifies the map by an additive shift proportional to C. Expanding (14) to linear order in $C\ll D_{A,B}$, we obtain:

(15)${\overline{\alpha }}^{\prime }={\displaystyle \frac{1}{\overline{\alpha }}}+{\displaystyle \frac{C}{D_B}}+\mathcal{O}\left(C^2\right),$

showing that tuning $C=D_B$ generates an exact unit shift on top of the inversion:

(16)$T_{\phi }:\alpha ⟼{\displaystyle \frac{1}{\alpha }}+1,T_{\phi }=\left(\begin{array}{cc}1& 1\\ 1& 0\end{array}\right),det{T}_{\phi }=-1.$

Unlike the involutive flip $S_{\phi }$, the shift $T_{\phi }$ is of infinite order: $T_{\phi }^n\left(\alpha \right)\ne \alpha$ for any $n>0$. Iterating the combined action of $S_{\phi }$ and $T_{\phi }$ yields the continued fraction orbit $\alpha ,1+1/\alpha$, $1+1/(1+1/\alpha )$, … which converges to the unique positive fixed point $\phi$.

In this formulation, the antisymmetric coupling k enforces an instantaneous exchange between entropy and work channels, while the tuned cross-correlation $C=D_B$ injects a unit shift in favor of the dissipative component at each swap. Together, they realize the modular transformation $T_{\phi }$, completing the Möbius generator pair $\{S_{\phi },T_{\phi }\}$ that underpins Dynamic Balance.

Importantly, $D_B$ is a measure of intrinsic fluctuations generated by the entropy sector B. Matching the cross-correlation to this noise level, $C=D_B=〈{\xi }_A{\xi }_B〉$, requires no fine-tuning—it naturally arises in open systems where both A and B couple to a common bath (e.g., phonons coupling electronic and lattice currents, or turbulent eddies affecting both kinetic and thermal flows). In such cases, the cross-covariance self-adjusts to this modularly resonant value.

Lemma 1. 

The antisymmetric Onsager coupling $L_{AB}=-L_{BA}$ enforces reversible entropy exchange and generates the modular involution $S_{\phi }:\alpha \mapsto {\phi }^2/\alpha$, preserving phase-space symmetry. Irreversibility arises from the recursive entropy shift $T_{\phi }:\alpha \mapsto 1+1/\alpha$, which breaks time symmetry and drives the system down entropy-curved geodesics toward φ. Thus, the Onsager matrix’s structure maps precisely onto the modular generators: $S_{\phi }$ reflects conservation; $T_{\phi }$, dissipation.

2.6. Convex Lyapunov Functional Invariant Under $S_{\phi }$

To quantify the thermodynamic “distance” of a system from the equilibrium fixed-point, we define a convex Lyapunov functionals $R\left(\alpha \right)$, invariant under Möbius involution $S_{\phi }:\alpha \mapsto {\mathsf{\Lambda }}^2/\alpha$. This function serves as a generalized energy landscape: it diverges as $\alpha \to 0^{+},\infty$, and attains its unique global minimum at $\alpha =\mathsf{\Lambda }$, the fixed point of $S_{\phi }$:

(17)$R\left(\alpha \right)={\displaystyle \frac{4}{{\phi }^2}}\left[{\mathcal{C}}^2\left(\alpha \right)-{\mathcal{C}}^2\left(\phi \right)\right]={\left({\displaystyle \frac{\alpha }{\mathsf{\Lambda }}}-{\displaystyle \frac{\mathsf{\Lambda }}{\alpha }}\right)}^2,R″\left(\alpha \right)>0,R\left(\alpha \right)\stackrel{\alpha \to 0^{+},\infty }{\to }+\infty .$

A complete proof of convexity and invariance is provided in the Supplementary Information. Embedding this scalar cost in a spatial domain $\mathsf{\Omega }\subset {\mathbb{R}}^d$ defines a free-energy functional:

(18)$\mathcal{F}\left[\alpha \right]={\int }_{\mathsf{\Omega }}\left[{\displaystyle \frac{\kappa }{2}}{|\nabla \alpha |}^2+R\left(\alpha \right)\right]d^{d}x,$

where the diffusivity $\kappa$ enforces local smoothing of entropy gradients. Taking gradient-descent dynamics (Model-A in the classification of [25]) yields the nonlinear reaction–diffusion equation:

(19)${\partial }_t\alpha =-\mathsf{\Gamma }{\displaystyle \frac{\delta \mathcal{F}}{\delta \alpha }}=\mathsf{\Gamma }\left[\kappa {\nabla }^2\alpha -{\partial }_{\alpha }R\left(\alpha \right)\right],$

where $\mathsf{\Gamma }>0$ sets the relaxation rate. For Neumann or periodic boundary conditions, the energy functional decays monotonically:

(20)$\dot{\mathcal{F}}={\int }_{\mathsf{\Omega }}{\displaystyle \frac{\delta \mathcal{F}}{\delta \alpha }}{\partial }_t\alpha =-\mathsf{\Gamma }{\int }_{\mathsf{\Omega }}{\left({\displaystyle \frac{\delta \mathcal{F}}{\delta \alpha }}\right)}^2{d}^{d}x\le 0,$

Thus, every trajectory evolves irreversibly toward the unique global minimizer of $R\left(\alpha \right)$, with $\dot{\mathcal{F}}\le 0$ ensuring entropy coherence and $\phi$-stability throughout the dynamics [41].

2.7. Common Fixed Point and Identification of $\mathsf{\Lambda }=\phi$

Steady, spatially uniform solutions of the entropy PDE (19) must satisfy ${\partial }_{\alpha }R\left({\alpha }^{★}\right)=0$, i.e., they must lie at extrema of the Lyapunov cost function $R\left(\alpha \right)$. Given that R is symmetric under the Möbius involution $S_{\phi }\left(\alpha \right)={\mathsf{\Lambda }}^2/\alpha$, any two points $\alpha$ and $S_{\phi }\left(\alpha \right)$ lie on the same energy contour: $R\left(S_{\phi }\left(\alpha \right)\right)=R\left(\alpha \right)$. For a strictly convex function, distinct points on the same contour cannot both be minima. Hence, if a minimum is fixed under $S_{\phi }$, the only possibility is that $S_{\phi }\left({\alpha }^{★}\right)={\alpha }^{★}$. This condition implies:

(21)${\alpha }^{★2}-{\mathsf{\Lambda }}^2=0.$

Additionally, the modular shift $T_{\phi }\left(\alpha \right)=1+1/\alpha$ is self-similar and infinite-order. Requiring that ${\alpha }^{★}$ also be a fixed point of $T_{\phi }$ yields:

(22)${\alpha }^{★2}-{\alpha }^{★}-1=0.$

Solving both equations, the only consistent, positive solution is:

(23)$\begin{array}{|c|}\hline \mathsf{\Lambda }={\alpha }^{★}=\phi \\ \hline\end{array}.$

Thus, the golden ratio emerges non-perturbatively from the interplay of Möbius symmetry and convex geometry. Once these symmetries are imposed, the fixed point ${\alpha }^{★}=\phi$ is uniquely selected as the global minimum of any boundary-divergent Lyapunov cost function $R\left(\alpha \right)$. The golden attractor arises under the following universal conditions:

(i). Two irreducible power channels $A,B$ forming the entropy flux field $\alpha =A/B$;

(ii). Möbius inversion symmetry $S_{\phi }:\alpha \mapsto {\phi }^2/\alpha$;

(iii). Self-similar translation symmetry $T_{\phi }:\alpha \mapsto 1+1/\alpha$;

(iv). A strictly convex Casimir functional $\mathcal{C}\left(\alpha \right)$ diverging at $\alpha \to 0^{+},\infty$.

Hence, convexity of $\mathcal{C}\left(\alpha \right)$ selects the extremum, and Möbius invariance under the modular subgroup $\langle S_{\phi },T_{\phi }\rangle \subset \mathrm{PGL}(2,\mathbb{Q}\left(\sqrt{5}\right))$ forces that extremum to lie precisely at the golden fixed point ${\alpha }^{★}=\phi$. This geometric origin underlies the appearance of Fibonacci sequences, self-similar structures, and golden spirals in driven-dissipative systems.

3. Thermodynamic Cost Function and Relaxation Dynamics

In Section 2, we showed that the Möbius involution symmetry $S_{\phi }$ defines the geometry of the entropy manifold by selecting the $\phi$-invariant Casimir contours, while the minimum of the Casimir $\mathcal{C}\left(\alpha \right)$ defines the origin of the gradient flow. From this point, all dynamical quantities—including the gradient, curvature, torsion, and evolutionary pathways of the entropy field—follow. We now study the dynamic trade-off between coherent energy flow and irreversible entropy production using our smooth, strictly convex cost function $R\left(\alpha \right)\in C^{\infty }$, whose second derivative $R^{″}\left(\phi \right)=8/{\phi }^2>0$ confirms a unique global minimum.

(24)$R\left(\alpha \right)={\displaystyle \frac{4}{{\phi }^2}}\left[{\mathcal{C}}^2\left(\alpha \right)-{\mathcal{C}}^2\left(\phi \right)\right]{\left({\displaystyle \frac{\alpha }{\phi }}-{\displaystyle \frac{\phi }{\alpha }}\right)}^2,\mathrm{with}R\left(\alpha \right)\to +\infty \mathrm{as}\alpha \to 0^{+},\infty .$

Physically, this cost penalizes both excessive dissipation ($\alpha \to 0$) and excessive energy retention ($\alpha \to \infty$), enforcing a Goldilocks balance exactly at the golden ratio (see Figure 1). This non-equilibrium potential drives every initial profile $\alpha (\mathbf{x},0)\in (0,\infty )$ monotonically toward the uniform attractor $\alpha (\mathbf{x},t)\to \phi$ as $t\to \infty$.

Parameter-Free Experimental Invariants

Linearizing the entropy dynamics (19) about ${\alpha }^{★}=\phi$, we obtain a local relaxation rate $\mu =R^{″}\left(\phi \right)=8/{\phi }^2$. This leads to three universal, dimensionless invariants of the entropy flux field:

The spatial scale $\xi$ sets the diameter of a coherent patch in which energy and entropy fluxes remain tightly coupled, while $\tau$ characterizes the re-equilibration time following a disturbance. Their ratio ${\xi }^2/\tau =\kappa \mathsf{\Gamma }$ is a renormalization-group invariant: If a coherent patch doubles in linear size, its relaxation time quadruples.

These parameter-free invariants have been observed (within experimental uncertainty) in diverse driven systems, including microbial metabolism and enzyme maintenance, cortical energy balance in active brains, turbulent vortex structures and rotational eddies, vascular branching and phyllotactic lattice development. At the attractor ${\alpha }^{★}=\phi$, the energy flux decomposes canonically:

(25)${\displaystyle \frac{T\dot{S}}{\dot{E}}}={\displaystyle \frac{1}{\phi }}\approx 0.618,{\displaystyle \frac{\dot{E}-T\dot{S}}{\dot{E}}}={\displaystyle \frac{1}{{\phi }^2}}\approx 0.382.$

suggesting that in any system where energy is optimally partitioned between reversible work and irreversible fluxes, the characteristic balance is as follows:

About 61.8% of energy is thermal entropy ($T\dot{S}$).

About 38.2% of energy is effective free energy ($\dot{E}-T\dot{S}$).

This breakdown matches empirical observations across biological and physical systems. In studies of metabolism, growth, and organismal energetics, a $~60–70\%$ dissipation fraction is consistently reported, with the remainder channelled into constructive output [33,42,43,44,45,46,47].

Example 1 

(Bathtub Whirlpool Convexity). Consider a bathtub that is simultaneously driven by a thin jet of water from the tap (power input $\dot{E}$) and dissipates energy through the drain, where viscous friction converts mechanical energy into heat ($T\dot{S}$). At an intermediate, self-regulated setting, the inflow and outflow balance so that ${\alpha }^{★}\simeq \phi ,$ and the water organizes into a stable, logarithmic vortex with the golden-ratio pitch. The curved bathtub walls play the role of the convex Lyapunov potential $R\left(\alpha \right)$, funneling the dynamics toward the single minimum. Thus, the everyday whirlpool illustrates how a bounded, strictly convex “energy landscape” produces a visible, scale-invariant pattern.

4. Discrete Markov Realization of the Flux–Ratio Dynamics

To demonstrate that the continuous Lyapunov dynamics (19) can emerge from an underlying microscopic process, we construct a coarse-grained, one0dimensional birth–death Markov chain with N discrete states

(26)$\{{\alpha }_i=i\mathsf{\Delta }\alpha \left|i=1,\dots ,N\right\},\mathsf{\Delta }\alpha >0,$

so that ${\alpha }_{min}=\mathsf{\Delta }\alpha$ and ${\alpha }_{max}=N\mathsf{\Delta }\alpha$. A threshold index $i_{\mathrm{th}}$ defines an instability cutoff beyond which avalanches (relaxation events) are triggered.

Let $P_i\left(t\right)=Pr\left[\alpha \left(t\right)={\alpha }_i\right]$, and define the probability vector $\mathbf{P}\left(t\right)={(P_1,\dots ,P_N)}^{\top }$ normalized such that ${\sum }_i{P}_i=1$. The system evolves according to a continuous-time master equation:

(27)${\dot{P}}_i=\sum _{j\ne i}{W}_{ij}{P}_j-\left(\sum _{j\ne i}{W}_{ji}\right)P_i,$

where $W_{ij}\ge 0$ for $i\ne j$ denotes transition rates, and $W_{ii}=-{\sum }_{j\ne i}{W}_{ji}$ ensures conservation of total probability. We now specify two elementary processes:

(i). Slow drive (energy input): ${\alpha }_i\to {\alpha }_{i+1}$ at constant rate $v_{\alpha }>0$ $\left(W_{i+1,i}=v_{\alpha }\right)$.

(ii). Avalanche relaxation (entropy release): for any $m\ge 1$ and $i>i_{\mathrm{th}}$: ${\alpha }_i\to {\alpha }_{i-m}$ at rate $W_{i-m,i}={\nu }_m\mathsf{\Theta }(i-i_{\mathrm{th}})$, where ${\nu }_m$ is the probability per unit time of a downward jump of size m. Reflecting boundary imposed $W_{0,i}=W_{N+1,i}=0$.

These rules form a one-dimensional analog of sandpile toppling [48,49]. The generator in (27) is irreducible, ensuring the existence of a unique stationary distribution ${\mathbf{P}}^{(\infty )}$ exists. For $i\ge i_{\mathrm{th}}$, the steady-state satisfies

(28)${\displaystyle \frac{P_{i+1}^{(\infty )}}{P_i^{(\infty )}}}={\displaystyle \frac{v_{\alpha }}{\nu }},\nu :=\sum _{m\ge 1}{\nu }_m,$

This leads to a geometric tail above threshold and a unimodal profile peaking at $i^{★}\simeq i_{\mathrm{th}}+\nu /v_{\alpha }$, which, in the continuum limit, converges to ${\alpha }^{★}=\phi$.

4.1. Continuum Limit and Fokker–Planck Correspondence

Letting $\alpha =i\mathsf{\Delta }\alpha$ and defining the probability density $P(\alpha ,t)=P_i\left(t\right)/\mathsf{\Delta }\alpha$, we expand (27) to second order in $\mathsf{\Delta }\alpha \ll 1$ using the Kramers–Moyal expansion [50]:

(29)${\partial }_{t}P=-{\partial }_{\alpha }\left[v_{\alpha }P\right]+{\partial }_{\alpha }^2\left[{\mathcal{D}}_{\alpha }P\right]+\mathcal{O}\left({\left(\mathsf{\Delta }\alpha \right)}^3\right),$

with drift $v_{\alpha }$ and diffusion coefficient ${\mathcal{D}}_{\alpha }=\frac{1}{2}\mathsf{\Delta }\alpha {\sum }_m{m}^2{\nu }_m$.

Comparing with the DB field equation (19) gives the dictionary:

(30)$v_{\alpha }=\mathsf{\Gamma },{\mathcal{D}}_{\alpha }=\mathsf{\Gamma }\kappa .$

As $\mathsf{\Delta }\alpha \to 0$, all higher-order Kramers–Moyal terms vanish, and the Markov chain converges exactly to the deterministic Lyapunov flow. For single-step avalanches ${\nu }_m={\nu }_1{\delta }_{m1}$ and $v_{\alpha }={\nu }_1$, the stationary distribution becomes a discrete Gaussian peaked at $i^{★}\simeq \phi /\mathsf{\Delta }\alpha$. In this regime, the avalanche-size distribution follows $P\left(S\right)\propto S^{-3/2}$, matching the mean-field exponent of sandpile models and the scaling observed in critical neural and condensed matter systems [22,48,51,52,53].

4.2. Parameter Dictionary

Here, $\kappa$ controls the spatial propagation of entropy imbalance—interpretable as thermal conductivity, stiffness, or axonal spread—while $\mathsf{\Gamma }$ sets the local rate of relaxation, e.g., via viscosity, phonon damping, or synaptic recovery.

Example 2 

(Sandpile avalanche criticality). In the classical Abelian sandpile model, a slow “rain” of grains, added at rate $v_{\alpha }$, builds a heap until the local slope at some site exceeds a threshold height $i_{\mathrm{th}}$. The site then topples, redistributing one grain to each neighbor; the relaxation may propagate and produce an avalanche whose size S (total topplings) obeys the mean-field law. Our birth–death chain is the direct energy-flux analogue. Slow drive adds a grain of “usable power” $\mathsf{\Delta }\alpha$ to the system, mirroring the external input that steepens the pile. Once α crosses the instability line $i_{\mathrm{th}}$, a stochastic event of size $m=1,2,\dots$ transfers $m\mathsf{\Delta }\alpha$ from the work channel A into the dissipation channel B. This implements the toppling rule in energy–entropy space. Because the drive $v_{\alpha }$ and the avalanche kernel ${\nu }_m$ are held fixed while the accessible state space extends to arbitrarily large i, the Markov process self-organizes to a stationary distribution peaked at ${\alpha }^{★}\simeq \phi$. The continuum limit reproduces the gradient-flow PDE (19), and the avalanche–size vstatistics approaches the same power law, as in the canonical sandpile SOC [22,48,53].

5. Modular Symmetry and Non-Equilibrium Field Theory

Having shown in Section 2 and Section 3 that the dual Möbius symmetry uniquely selects $\mathsf{\Lambda }=\phi$, and in Section 4 that the discrete Markov chain recovers the continuous Lyapunov flow as $\mathsf{\Delta }\alpha \to 0$, we now extend the theory to include fluctuations. We demonstrate that the golden-ratio attractor remains symmetry-protected in the stochastic setting, and derive exact selection rules and Ward identities [54,55,56,57].

5.1. Stochastic Dynamic-Balance Equation

We promote the entropy field $\alpha \equiv \alpha (\mathbf{x},t)$ to a stochastic variable over a d spatial dimensional spatial domain, governed by the Langevin equation:

(31)${\partial }_t\alpha =\mathsf{\Gamma }\left[\kappa {\nabla }^2\alpha -{\partial }_{\alpha }R\left(\alpha \right)\right]+\eta ,\langle \eta (\mathbf{x},t)\eta (\mathbf{x}’,t’)\rangle =2D,{\delta }^d(\mathbf{x}-\mathbf{x}’),\delta (t-t’).$

Here, $\eta$ is Gaussian white noise of strength D, and $\mathsf{\Gamma }>0$ is the relaxation rate. Importantly, D and $\mathsf{\Gamma }$ are independent parameters: the system is driven and does not obey detailed balance. The relation ${\mathcal{D}}_{\alpha }=\kappa \mathsf{\Gamma }$ seen in the Fokker–Planck limit is not a fluctuation–dissipation theorem, but an RG-fixed constraint imposed by Möbius symmetry after coarse-graining. Setting $D=0$ recovers pure Lyapunov descent toward the golden attractor ${\alpha }^{★}=\phi$.

5.2. MSRJD Path Integral Formalism

To study the statistical field theory, we follow the Martin–Siggia–Rose–Janssen–de Dominicis (MSRJD) formalism [58], introducing a response field $\widehat{\alpha }(\mathbf{x},t)$ and writing the generating functional:

(32)$\mathcal{Z}[J,\widehat{J}]=\int \mathcal{D}\alpha \mathcal{D}\widehat{\alpha }exp\left(-S[\alpha ,\widehat{\alpha }]+\int d^{d}xdt[J\alpha +\widehat{J}\widehat{\alpha }]\right).$

with the stochastic action:

(33)$S[\alpha ,\widehat{\alpha }]=\int d^{d}xdt\left\{\widehat{\alpha }\left[{\partial }_t\alpha -\mathsf{\Gamma }(\kappa {\nabla }^2\alpha -{\partial }_{\alpha }R\left(\alpha \right))\right]+D{\widehat{\alpha }}^2\right\}.$

Correlation functions of $\alpha$ and $\widehat{\alpha }$ follow from functional derivatives with respect to $J,\widehat{J}$. We adopt the Itô convention: the Jacobian from $\delta [{\partial }_t\alpha -\dots ]$ is a constant and set to unity.

5.3. Modular Symmetries and Ward Identities

Only involutive Möbius transformations commute with the stochastic time-reversal symmetry $(t,\alpha ,\widehat{\alpha })\to (-t,\alpha ,-\widehat{\alpha })$, leaving the action invariant. Higher-order modular elements break this symmetry by flipping the sign of the MSR term and thus cannot survive in the stochastic theory [54,57,58,59].

We embed both Möbius maps as exact symmetries of $S[\alpha ,\widehat{\alpha }]$, with transformations chosen to preserve the path integral measure ($\mathcal{D}\alpha \mathcal{D}\widehat{\alpha }$):

Self-dual flip $S_{\phi }:\alpha \mapsto {\displaystyle \frac{{\phi }^2}{\alpha }}:$

(34a) $\begin{array}{c}\hfill \delta \alpha =\epsilon \left(\frac{{\phi }^2}{\alpha }-\alpha \right),\delta \widehat{\alpha }=-\epsilon \left(\frac{{\phi }^2}{{\alpha }^2}+1\right)\widehat{\alpha },\end{array}$

Self-similar shift $T_{\phi }:\alpha \mapsto 1+{\displaystyle \frac{1}{\alpha }}:$

(34b) $\begin{array}{c}\hfill \delta \alpha =\epsilon \left(1+\frac{1}{\alpha }-\alpha \right),\delta \widehat{\alpha }=-\epsilon \left(\frac{1}{{\alpha }^2}-1\right)\widehat{\alpha },\end{array}$

By direct substitution, we confirm that $S[\alpha +\delta \alpha ,\widehat{\alpha }+\delta \widehat{\alpha }]=S[\alpha ,\widehat{\alpha }]$. Applying these variations to the path-integral and requiring $\delta \mathcal{Z}=0$ yields exact Ward identities for any operator $\mathcal{O}[\alpha ,\widehat{\alpha }]$ [56].

$S_{\phi }$–Ward identity:

(35a) $\begin{array}{c}\hfill \int d^{d}xdt〈\widehat{\alpha }(\mathbf{x},t)\left(\frac{{\phi }^2}{\alpha (\mathbf{x},t)}-\alpha (\mathbf{x},t)\right)\mathcal{O}〉=0.\end{array}$

$T_{\phi }$–Ward identity:

(35b) $\begin{array}{c}\hfill \int d^{d}xdt〈\widehat{\alpha }(\mathbf{x},t)\left(1+\frac{1}{\alpha (\mathbf{x},t)}-\alpha (\mathbf{x},t)\right)\mathcal{O}〉=0.\end{array}$

Theorem 2. 

Let ${\mathcal{O}}_m\left(\alpha \right)\propto {\alpha }^m$ carry modular charge m under $S_{\phi }$ or $T_{\phi }$. Then for any nonzero n-point correlator,

(36) $〈{\mathcal{O}}_{m_1}\dots {\mathcal{O}}_{m_n}〉\ne 0⟹\sum _{i=1}^n{m}_i=0.$

Choosing $\mathcal{O}=1$ shows fluctuations alone cannot bias $\alpha$ away from $\phi$: the golden attractor is symmetry-protected. Any process violating modular charge conservation is strictly forbidden. Thus, the same Möbius generators that extremize the Lyapunov functional survive in the stochastic theory and protect the attractor through exact Ward identities [57,58,59]. This parallels how conformal symmetry and S–duality constrain correlators in equilibrium quantum field theory (QFT) [60,61]. In the non-equilibrium framework, the analogous role is played by the discrete modular symmetries of entropy fluxes dynamics.

Example 3. 

Let us evaluate the two Ward identities (35a) and (35b) using the insertion $\mathcal{O}\left[\alpha \right]={\alpha }^n$, $n\in {\mathbb{N}}_0$. Substituting into the identities and integrating by parts yields the following exact constraints:

(37) $({\phi }^2-\langle {\alpha }^2\rangle )\langle {\alpha }^{n-1}\rangle =0,(1+\langle {\alpha }^{-1}\rangle -\langle \alpha \rangle )\langle {\alpha }^n\rangle =0.$

For $n=0$, the first identity enforces $\langle {\alpha }^2\rangle ={\phi }^2$; for $n=1$, the second identity gives $\langle \alpha \rangle =\phi$. Since these constraints hold simultaneously for all $n\in {\mathbb{N}}_0$, they imply by recursion $\langle {\alpha }^m\rangle ={\phi }^m$ for every $m\in \mathbb{N},\forall m\in \mathbb{N}$, i.e., all moments are locked to powers of the golden ratio.

In particular, the one-point function satisfies $\langle \alpha \rangle =\phi$ exactly, to all orders in perturbative theory. This result illustrates how the golden fixed point is not just an attractor of deterministic dynamics but is symmetry-protected against all Gaussian (and weakly non-Gaussian) noise, due to modular invariance embedded in the stochastic field theory.

5.4. Quadratic Theory and Diffusive Pole

To analyze fluctuations around the golden attractor, we expand the entropy field near ${\alpha }^{★}=\phi$ as:

(38)$\alpha (\mathbf{x},t)=\phi +\delta \alpha (\mathbf{x},t),\left|\delta \alpha \right|\ll 1,$

Expanding the MSRJD action to quadratic order in $\delta \alpha$ and $\widehat{\alpha }$, we obtain:

(39)$S_2=\int d^{d}xdt\left\{\widehat{\alpha }\left[{\partial }_t-\mathsf{\Gamma }(\kappa {\nabla }^2-\mu )\right]\delta \alpha +D{\widehat{\alpha }}^2\right\},\mu =R^{″}\left(\phi \right)={\displaystyle \frac{8}{{\phi }^2}}.$

The retarded Green’s function, defined as the response to a perturbation at $({\mathbf{x}}^{\prime },t^{\prime })$, is governed by the inverse of the operator in the linear term (neglecting the ${\widehat{\alpha }}^2$ term which contributes only to $\langle \alpha \alpha \rangle$, not to $G_R$):

(40)$G_R(\omega ,\mathbf{q})={\displaystyle \frac{1}{i\omega +\mathsf{\Gamma }\left({\kappa \left|\mathbf{q}\right|}^2+\mu \right)}}.$

This propagator has a single diffusive pole located at:

(41)${\omega }^{★}=-i\mathsf{\Gamma }(\kappa q^2+\mu ).$

At long wavelengths ($q\to 0$), the dispersion relation reduces to $\omega \sim -i\mathsf{\Gamma }\mu$, giving the familiar diffusive scaling $\omega \sim -i{q}^2$, with dynamical exponent $z=2$. To geometrize the fluctuation modes, consider the argument of the complex frequency ${\omega }^{★}$. Define:

(42)$\vartheta \left(q\right)=arg\left(i{\omega }^{★}\right)={tan}^{-1}\left({\displaystyle \frac{\mu }{\kappa q^2}}\right).$

which interpolates from a purely relaxational mode $q\to 0,\vartheta \to \pi /2$ at long wavelengths (IR limit) to a purely diffusive mode $q\to \infty ,\vartheta \to 0$ at short wavelengths (UV limit). At the crossover scale $q_c=\sqrt{\mu /\kappa }$, we find ${\omega }^{★}=-i\mathsf{\Gamma }\mu (1+i)$, so that $\vartheta \left(q_c\right)=\pi /4$. This ${45}^{\circ }$ phase angle represents perfect balance between real and imaginary parts—between reversible energy and irreversible dissipation. In control theory, this is the critical damping angle. In real space, it corresponds to a logarithmic spiral with constant pitch:

(43)$r\left(\theta \right)=a{e}^{b\theta },b={\displaystyle \frac{\Re {\omega }^{★}}{\Im {\omega }^{★}}}=1,$

or equivalently, $r\left(\theta \right)=a{\phi }^{\theta /\pi }$, using the identity $e^{\pi }={\phi }^{\pi /ln\phi }$. Thus, the pole’s damping ratio directly maps onto the golden spiral—a physical manifestation of modular balance in space–time.

Interpretation of $D^{★}=\kappa \mathsf{\Gamma }$: A Modular Einstein Relation

The effective macroscopic diffusion constant for fluctuations in the entropy field $\alpha (\mathbf{x},t)$ is:

(44)$D^{★}=\kappa \mathsf{\Gamma }.$

where $\kappa$ measures spatial coupling, or how strongly neighboring regions equilibrate $\alpha$. $\mathsf{\Gamma }$ is the local response speed, dictating how fast $\alpha$ relaxes back to $\phi$.

Under coarse-graining—$\kappa \sim b^{-2},\mathsf{\Gamma }\sim b^{+2}$—so their product remains scale-invariant. This mirrors the Einstein relation $D=\mu k_{B}T$, where mobility and temperature compensate. Here, modular symmetry replaces equilibrium as the principle that constrains the fluctuation–dissipation balance.

5.5. Dynamical-Exponent Landscape

At Gaussian (tree-level) approximation, the dynamic exponent is $z=2$ (see Appendix B). However, beyond mean-field, nonlinearities introduce loop corrections. In particular, the one-loop vertex from the cubic term $\sim {\left({\partial }_{\alpha }R\right)}^2\widehat{\alpha }\delta {\alpha }^2$ renormalizes the quadratic propagator and alters the effective dispersion relation. Importantly, modular symmetry structure determines the form and strength of these corrections:

Full modular symmetry ($S_{\phi }$ and $T_{\phi }$ preserved; Onsager matrix antisymmetric):

The nonlinearity couples left- and right-moving modes through dual Möbius rotations. The retarded self-energy inherits the golden eigen-angle, and the dispersion becomes $\omega \sim q^{\phi }$, i.e., $z=\phi$. This is the fully symmetry-protected golden dynamic exponent, characteristic of balanced, $\phi$-stabilized flows.

Self-dual line ($S_{\phi }$ preserved, $T_{\phi }$ broken): One modular charge is violated. The resulting loop integral resembles the modified KPZ class with $z=\sqrt{2}$ [62]. This phase retains $\phi$-inversion symmetry but lacks self-similarity, leading to intermediate roughening.

No modular symmetry (Onsager symmetric): Standard KPZ-type scaling emerges, depending on conservation laws $z\in \{1,3/2,2\}$, spanning ballistic, superdiffusive, and diffusive behavior [57,63].

Strong disorder/broken detailed balance (e.g., random-field landscape, Sinai-type potential): The system becomes glassy and subdiffusive. Scaling slows to $z=1/2$, typical of Sinai creep in 1D random environments (disorder induced irreversibility) [64,65].

Theorem 3. 

The golden point ${\alpha }^{★}=\phi$ is the unique real fixed point of the Möbius transformation $T\circ S\in \mathrm{PSL}(2,\mathbb{Z})$. This transformation has order 3, satisfying ${\left(TS\right)}^3=\mathcal{I}$ geometrically classifying φ as a real order-3 hyperbolic fixed-point. It is the projection to ${\mathbb{R}}_{+}$ of the complex elliptic fixed point $\tau =e^{2\pi i/3}\in \mathbb{H}$, the triangle vertex of the modular tiling with internal angle $\pi /3$. At this point, the full modular symmetry $\mathrm{PSL}(2,\mathbb{Z})$ acts conformally on the entropy balance field $\alpha \in {\mathbb{R}}_{+}$, both generators $S_{\alpha }$ and $T_{\alpha }$ are preserved, the entropy Casimir $\mathcal{C}\left(\alpha \right)={\displaystyle \frac{1}{2}}(\alpha +{\phi }^2/\alpha )$ is minimized, the RG-invariant diffusivity $D_z={\xi }^z/\tau$ attains optimal scaling with $z=\phi$. Therefore, ${\alpha }^{★}=\phi$ defines the unique conformal and modular RG fixed point on ${\mathbb{R}}_{+}$, stabilized by Möbius symmetry.

This dynamical-exponent landscape provides a symmetry-based RG classification of entropy field dynamics. Each value of $z_n=2cos(\pi /n)$ corresponds to a distinct Coxeter–modular class, encoding a specific pattern of symmetry breaking, Onsager structure, and renormalization flow. The dynamic exponent governs the anisotropic scaling between space and time: $t\sim {\xi }^{z_n}$, setting the causal geometry of entropy transport. This classification is intrinsically tied to geodesic structure in modular and hyperbolic spaces. Each Coxeter class $I_2\left(n\right)$ gives rise to a set of closed geodesics with angular period ${\displaystyle \frac{2\pi }{n}}$. The dynamical exponent $z_n$ encodes the geodesic curvature of the entropy flow in modular space.

The golden ratio exponent $z=\phi$ marks a universality class uniquely protected by full Möbius symmetry—where both modular generators $S_{\alpha }$ and $T_{\alpha }$ are preserved. At this fixed point, entropy flux is optimally balanced between reversible and irreversible components, minimizing the golden Casimir and maximizing the RG-invariant diffusivity $D_{\phi }={\xi }^{\phi }/\tau$.

This structure reveals that dynamical exponents are not just empirical parameters, but modular symmetry charges—labels of how entropy flow transforms under Coxeter–Möbius RG symmetries. The golden point $z=\phi$ represents the only uncharged (invariant) fixed point, where entropy flow is modularly neutral, and balance is perfectly preserved. At this point, the entropy field is invariant under both $S_{\alpha }$ and $T_{\alpha }$, and the full group $\mathrm{PGL}(2,\mathbb{Q}\left(\sqrt{5}\right))$ acts as if it were locally conformal, stabilizing the entropy flow and minimizing the golden Casimir. Thus, conformal invariance emerges only at the golden point—and it is not generic, but modularly protected.

6. Discussion

The symmetric Möbius involution $\alpha \mapsto {\phi }^2/\alpha$ encodes a modular duality intrinsic to the entropy field, mapping energy–entropy configurations into their golden-conjugate counterparts. The fixed point ${\alpha }^{★}=\phi$ uniquely minimizes the convex Casimir potential $\mathcal{C}\left(\alpha \right)$, defining the universal entropy balance point. This potential generates the intrinsic geometric landscape for entropy flow along geodesics in this landscape $\mathbb{H}/{\mathsf{\Gamma }}_0\left(N\right)$, constrained by Coxeter symmetry class $I_2\left(n\right)$, each dictating a distinct dynamic exponent $z_n$ and flow curvature. Conjugate pairs ($\alpha ,{\phi }^2/\alpha$) form modular dual observables—mirrored entropy states under the Möbius involution—that encode the field’s reversible flow structure. $\phi$ governs the flow of entropy through modular spacetime. Microscopic details set only the approach rate ($\mathsf{\Gamma }$, Section 3) and spatial coupling ($\kappa$, Section 4); all dimensionless observables are fixed by symmetry.

(45)$\left\{1/\phi :1/{\phi }^2,D=\kappa \mathsf{\Gamma }={\xi }^z/\tau ,\vartheta ={45}^{\circ }\right\}$

Thermodynamic partition: At ${\alpha }^{★}=\phi$, $1/\phi \simeq 61.8\%$ of inflowing power is dumped as heat; the remaining $1/{\phi }^2\simeq 38.2\%$ fuels coherent structure. Measured maintenance-vs-growth splits in microbes, plants, animals, and cortical tissue cluster near this ratio [42,66].

Renormalisation-group invariant: The diffusivity $D={\xi }^z/\tau =\kappa \mathsf{\Gamma }$ remains fixed under coarse graining, persisting across scales because the geometric structure is modularly renormalization-invariant: changes in length or time do not break the $\phi$-balance. In the diffusive limit, doubling the linear size of a coherent patch, therefore, quadruples its relaxation time—exactly as the Einstein relation $D=\mu k_{\mathrm{B}}T$ couples diffusivity and mobility in equilibrium [67,68].

Emergent geometry: Linear response yields a ${45}^{\circ }$ eigen-angle in the complex frequency plane, mapping to the golden-pitch logarithmic spiral in real space $tan\psi =1/\phi (\psi \simeq 31.7^{°})$. This explains why vortices in rotating turbulence, hurricanes, galactic arms, and phyllotactic patterns share the same spiral pitch [5,6,7,9,69,70,71,72,73].

Dynamic-exponent landscape: Beyond mean field, loop corrections shift the exponent z: Full modular symmetry ($S_{\phi }$ + $T_{\phi }$ and antisymmetric Onsager matrix) produces $z=\phi$; breaking $T_{\phi }$ alone gives $z=\sqrt{2}$; and removing both yields the KPZ/Levy line, $z\in \{1,3/2,2\}$ (Table 1).

A recurrent question is how the three diffusion constants that appear in the paper—the microscopic channel variances $D_{A,B}$ in Section 2, the lattice diffusion ${\mathcal{D}}_{\alpha }$ of the birth–death chain in Section 4, and the stochastic amplitude $D^{★}$ in the MSR action Section 5—are related. The link is the matched–bath condition $C=〈{\xi }_A{\xi }_B〉=D_B$ between the cross-covariance C to the entropy-sector variance $D_B$, which ensures that every dissipative kick is transmitted to the work channel with unit efficiency. Under this condition, the slow balance field inherits a single effective noise level $D=\kappa \mathsf{\Gamma }$, where $\kappa$ is the stiffness that spreads deviations of $\alpha$ and $\mathsf{\Gamma }$ is the local relaxation rate. Coarse-graining the birth–death chain reproduces the same value, ${\mathcal{D}}_{\alpha }\to \kappa \mathsf{\Gamma }$, and the continuum MSR functional keeps it unchanged. Thus the Einstein-like product $\kappa \mathsf{\Gamma }$ survives intact from the microscopic Ornstein–Uhlenbeck description to the macroscopic field theory, providing an internally consistent “noise ladder” that underlies the golden-ratio universality class [74,75,76]. Rather than tuning temperature or pressure to a critical point $(g-g_c)$, the system self-organizes its energy–entropy flux gradient to the golden fixed point $(\alpha -\phi )$. From biochemical chirality to turbulent galaxies, from neural avalanches to non-Fermi-liquid metals, the same Möbius duality and recursive modular geometry encode how entropy organizes structure across scales. Each $\phi$-deviation acts as a local “entropy curvature”, and its relaxation under $\mathcal{C}\left(\alpha \right)$ guides systems back toward balance. This self-organized modularity manifests as fractals, scaling laws, and golden-pitch spirals.

$(D_A,D_B)\stackrel{\mathrm{coarse}\text{-}\mathrm{grain}}{\to }{\mathcal{D}}_{\alpha }\stackrel{\mathrm{Kramers}–\mathrm{Moyal}}{\to }D=\kappa \mathsf{\Gamma }.$

6.1. Two-Fluid Decomposition and Quantum-Critical Universality

The order-2 modular involution $S_{\phi }$ mandates a fundamental dichotomy in any coarse-grained entropy flux: a component that recirculates as coherent structure (channel A) and a component that dissipates as heat (channel B). This partition is enforced by the antisymmetric Onsager matrix, whose off-diagonal coupling $L_{AB}=-L_{BA}$ induces a reactive interaction between the two channels—analogous to a Poisson bracket. Channels A and B are thus dynamically distinct, with separate characteristic time scales. In quantum many-body systems, this modular two-fluid structure mirrors the hydrodynamic separation between slowly-relaxing, advective “coherent” modes (e.g., momentum and supercurrents) and rapidly equilibrating “incoherent” modes [77,78,79]. Importantly, the modular symmetry acts simultaneously on:

on thermodynamic fluxes ($A\leftrightarrow B$);

on RG couplings $(g_1,g_2)\mapsto (g_2,g_1)$.

These dual actions intersect along the golden manifold $\mathcal{M}\phi$, where the system achieves full modular balance. Along $\mathcal{M}\phi$, the linearized RG flow matrix is unimodular, yielding eigenvalues in reciprocal pairs ($\lambda ,1/\lambda$). As a result, the six conventional static critical exponents collapse to just two RG invariants: the diffusivity $D=\kappa \mathsf{\Gamma }$, where $\kappa$ and $\mathsf{\Gamma }$ respectively encode microscopic stiffness and local relaxation rate (see Section 3 and Section 4); and the universal spiral angle $\vartheta ={45}^{\circ }$, which controls the flow geometry in entropy space and determines dynamic scaling trajectories. This is the geometric origin of the Kadanoff scaling relations. At this fixed-point manifold, all non-universal microscopic details—those not set by symmetry or dimensionality—wash out under RG coarse-graining. The critical sector becomes governed solely by the universal fixed-point data $(\nu ,\eta )$. At the golden point, the quadratic action becomes relativistic, and the dynamical exponent locks to $z=1$, signaling the emergence of Lorentz symmetry (Section 5). The recursive entropy flow generated by the modular shift $T_{\phi }$ becomes dynamically trivial at the golden fixed point—its action vanishes, as $\phi$ is invariant under this transformation. This halting of entropy recursion removes any preferred direction in entropy-time. As a result, the entropy field scales isotropically in space and time, a necessary and sufficient condition for emergent Lorentz invariance.

Crucially, scale-invariant quantum and thermal fluctuations are not a source of disorder—they are the generative mechanism of modular symmetry itself. The off-diagonal noise covariance $C=\langle {\xi }_A{\xi }_B\rangle$ acts as a reactive coupling between entropy channels, enabling the recursive transformation $T_{\phi }$. This stochastic cross-talk is essential: without it, the two entropy channels A and B remain dynamically decoupled, and the full duality cannot emerge. In this view, noise is not an error term to be averaged away—it is the very operator that lifts degeneracy and drives symmetry formation. The modular shift $T_{\phi }$ becomes dynamically active only through this correlated fluctuation structure, which recursively aligns the entropy field toward the attractor.

Why universality spans materials. Across all experimentally studied quantum critical points (QCPs), one observes two competing collective sectors ${\mathcal{O}}_{A,B}$: for example, superconductivity versus charge-density-wave order in cuprates, Kondo versus RKKY screening in heavy fermions, and vison–Majorana duality in Kitaev spin liquids. Tuning parameters such as doping, pressure, or magnetic field simply shift the balance $\mathsf{\Delta }\alpha =\alpha -\phi$ between these two entropy channels. This symmetry-governed partition explains why systems with widely varying microscopic Hamiltonians nonetheless share identical critical exponents—and why ultrafast pump–probe experiments consistently reveal a slow-plus-fast two-component relaxation process [80,81]. The very same order-2 flip $S_{\phi }$ embeds as a Weyl reflection of the non-crystallographic $H_2\subset E_8$ root-octagon lattice. This single ${\mathbb{Z}}_2$ reflection underlies criticality in systems as varied as: 1D Ising chain with $E_8$ quasiparticles [14,15], Kitaev honeycomb vison–Majorana duality [82,83], Fibonacci anyons-vacuum duality [19,84,85], charge–flux duality in quantum Hall edges and superconducting qubit devices [86,87], and near-horizon black-hole entropy dynamics [88,89]. This convergence illustrates that a single, order-2 modular duality governs the emergence of scale-free organization in systems ranging from condensed matter to quantum gravity. Though microscopic degrees of freedom differ, their algebraic backbone is isomorphic (see Appendix C).

6.2. Dynamic Balance in Gravity and Cosmology

Treating the expanding Universe as a fundamentally non-equilibrium systems–or at least as a two–channel entropy-flow network–reveals a striking modular structure. The comoving matter density scales as ${\rho }_m\propto a^{-3}$ while the apparent-horizon entropy grows as $S_{\mathcal{H}}\propto a^2$ [90]. This mismatch suggests a natural decomposition of the cosmic energy budget: group cold dark matter and baryons into an energy channel A, and horizon entropy—or more precisely, its flux—into an entropy channel B. Horizon entropy is a flux (rate of irreversible entropy production times temperature $B\equiv T_H{\dot{S}}_H$) at which the de Sitter (or black-hole) horizon pumps entropy out of the bulk. Vacuum energy is the integrated effect of that flux (a constant energy density with negative pressure). Specifically, the entropy flux across a cosmological horizon is $B\equiv T_H{\dot{S}}_H$ where where $T_H$ is the Gibbons–Hawking temperature. This horizon entropy flux behaves as a dissipative “outflow” from the bulk, continuously increasing as the Universe expands. The associated vacuum energy—a constant energy density with negative pressure—emerges as the integrated effect of this entropy flux. In the two-channel framework, we interpret channel B as the horizon-entropy output, while its equilibrium limit, once $\alpha =A/B$ relaxes to $\phi$, manifests as the observed vacuum energy density ${\rho }_{\mathsf{\Lambda }}$. They are two aspects of the same driven–dissipative process (see Appendix E).

(46)$A={\mathsf{\Omega }}_m,B={\mathsf{\Omega }}_{\mathsf{\Lambda }},\alpha \left(x\right)={\displaystyle \frac{A}{B}}$

Starting from the same Möbius-invariant Lyapunov cost functional $R\left(\alpha \right)$ and promoting the entropy flux field $\alpha \left(x\right)$ to a dynamical field in spacetime, yields the following gravitational action [91,92]:

(47)$\mathcal{S}[g,\alpha ]={\displaystyle \frac{1}{8\pi G_N}}{\int }_{\partial \mathcal{V}}\sqrt{h}{d}^{3}x-{\int }_{\mathcal{V}}\sqrt{-g}\left[\alpha \mathcal{R}-{\displaystyle \frac{\kappa }{2}}{(\nabla \alpha )}^2+R\left(\alpha \right)\right]d^{4}x,$

where $\mathcal{R}$ is the Ricci scalar, and the gradient term $\frac{1}{2}{(\nabla \alpha )}^2$ represents horizon elasticity or stiffness. Varying this action $\mathcal{S}$ with respect to $\alpha$ and $g_{\mu \nu }$ using the spatially flat Friedmann–Lemaître–Robertson–Walker metric,

(48)$d{s}^2=d{t}^2-a^2\left(t\right)d{\mathbf{x}}^2,H={\displaystyle \frac{\dot{a}}{a}},\mathcal{R}=6\left(2H^2+\dot{H}\right),$

produces the coupled system:

(49a)$\begin{array}{cc}\hfill & \ddot{\alpha }+3H\dot{\alpha }+R^{\prime }\left(\alpha \right)=0,\hfill \end{array}$

(49b)$\begin{array}{cc}\hfill & 3H^2=8\pi G_N\rho +\left[R\left(\alpha \right)-\alpha \mathcal{R}\right],\hfill \end{array}$

where $\rho$ is the matter energy density. Because $R^{\prime }\left(\alpha \right)=0$ vanishes only at ${\alpha }^{★}=\phi$, the dynamics naturally drive the system toward this golden-ratio attractor. At this point, the equation of state is $w=-1$, and no separate cosmological constant is required. Observations today, with ${\mathsf{\Omega }}_{\mathsf{\Lambda }}:{\mathsf{\Omega }}_m\approx 0.69:0.31$, lie within $10\%$ of the golden split $1/\phi :1/{\phi }^2$ [93]. Linearizing (49a) about $\alpha \left(t\right)=\phi +\delta \alpha \left(t\right)$ yields:

(50)$\ddot{\delta \alpha }+3H\dot{\delta \alpha }+R^{″}\left(\phi \right)\delta \alpha =0,R^{″}\left(\phi \right)={\displaystyle \frac{8{\mu }^4}{{\phi }^2}}.$

In the underdamped limit $H\ll {\mu }^2$, the entropy field $\delta \alpha$ oscillates with eigen-angle $\vartheta =arctan(1/\phi )={45}^{\circ }$, generating logarithmic spirals $r\left(\theta \right)\propto e^{\theta }$ with golden pitch. The corresponding fractal (box-counting) dimension is: $D_f=1+1/\phi \approx 1.618$ which matches the observed cosmic-web fractality on 5–100 Mpc scales [94].

The same involution $S_{\phi }$ governs both the entropy dynamics here and the Schwarzian action of nearly-$Ad{S}_2$ gravity (SYK), including the Euclidean “cigar” geometry [89,95]. Embedding the golden-ratio cost $R\left(\alpha \right)$ into the gravitational action thus:

Drives the Universe to a de-Sitter vacuum ${\alpha }^{★}=\phi$ with equation of state $w=-1$;

Predicts the dark-energy fraction ${\mathsf{\Omega }}_{\mathsf{\Lambda }}\approx 0.62$;

Reproduces the cosmic web’s spiral structure and fractal dimension $D_f=1+1/\phi$;

Links directly to AdS₂ modular dualities through the same Möbius involution.

The holographic principle ${\mathrm{AdS}}_n⟷{\mathrm{CFT}}_{n-1}$ is a manifestation of the intrinsic entropy field duality: particles in the bulk and field waves on the boundary.

7. Conclusions

We have shown that an order-2 Möbius involution $S_{\phi }$, together with a recursive shift $T_{\phi }$, generates a discrete, non-Abelian modular subgroup of $\mathrm{PGL}(2,\mathbb{Q}\left(\sqrt{5}\right))$ acting on the entropy flux field $\alpha (\mathbf{x},t)$. This modular action dynamically stabilizes dissipative systems at the golden fixed point ${\alpha }^{★}=\phi$, where recursive entropy flow halts, and Lorentz symmetry emerges with $z=1$. From this symmetry, we derived three dimensionless invariants: (1) a universal $62\%:38\%$ partition of entropy into dissipation and work; (2) a scale-invariant diffusivity coefficient $D=\kappa \mathsf{\Gamma }$; and (3) the golden-pitch logarithmic spiral that spans 15 orders of magnitude, from plant phyllotaxis to galactic arms. We derived Ward identities and symmetry selection rules that map the dynamical exponent landscape, setting clear conditions for the emergence of relativistic, KPZ, or anomalous scaling regimes. In each case, the fixed point corresponds to a modular symmetry class $I_2\left(n\right)$, characterized by dual flows and recursive balance.

Beyond statistical mechanics, these same modular symmetries enforce a two-fluid decomposition in quantum-critical matter, locking competing sectors into a geometric entropy flow governed by $\alpha$. Entropy recursion, driven by off-diagonal fluctuation covariance $C=\langle {\xi }_A{\xi }_B\rangle$, generates modular time symmetry and constrains the renormalization group via unimodular eigenflows. This matched-bath condition allows noise from the dissipative channel to drive structure in the coherent sector without destabilizing the golden attractor. This mechanism explains why, in high-$T_c$ superconductors and other strongly coupled quantum systems, it is not quasiparticles but scale-invariant fluctuations that provide the “pairing glue” [77,80].

Finally, we closed the loop from the entropy field symmetry to Legendre-conjugate pairs observables—duals under the entropy–energy exchange which map directly onto observable quantities. Thus, the conjugate structure of physics itself—thermodynamic pairs, RG exponents, quantum fields—is not postulated, but flows naturally from the modular dynamics of entropy. The flow of entropy is the master field from which emerge time–energy pairs, pressure–volume, density–curvature, and wave–particle dualities. Every conservation law is a balance law; every observable a trace of modular self-duality. From phyllotaxis, neurodynamics and turbulent flows, from superconducting gap dynamics to near-horizon entropy flow, from emergent Lorentz symmetry to Fibonacci anyon statistics, the same modular structure governs. Microscopic pairing, mesoscopic avalanche scaling, and cosmic acceleration are not separate phenomena—they are symmetry-equivalent expressions of one recursive, self-dual entropy balance optimization.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:

DB	Dynamic Balance
NESS	Non-Equilibrium Steady-State
PDE	Partial Differential Equation
ODE	Ordinary Differential Equation
PGL	Projective General Linear
RG	Renormalization Group
SOC	Self-Organized Criticality
CFC	Cross-Frequency Couplings
FLRW	Friedmann–Lemaître–Robertson–Walker metric
QCP	Quantum Critical Point

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figure and Table

Figure 1 Thermodynamic potential (or cost function) $R\left(\alpha \right)$ vs. $\alpha$. The divergences at $\alpha \to 0^{+},\infty$ represent strongly penalized boundary states. The stable, self-similar critical regime emerges uniquely at the global minimum ${\alpha }^{★}=\phi$. Minimizing $R\left(\alpha \right)$ does not mean the system is at zero net entropy production. Instead, it means it has found an optimal partition of energy vs. dissipation, optimizing both stability, efficiency and coherence in energy use, and preventing the system from falling into excessive disorder or excessive rigidity.

Table 1

“Periodic Table” of Dynamic Exponents $z_n=2cos(\pi /n)$, and RG invariant $D_z={\xi }^z/\tau$. The ✓ corresponds to preserved/invariant under RG and × indicates the generator is dynamically inactive (e.g., fixed-point degenerate), not symmetry-broken in the conventional sense.

Coxeter Class	${\mathit{S}}_{\mathit{\alpha }}$	${\mathit{T}}_{\mathit{\alpha }}$	Onsager ${\mathit{L}}_{\mathit{AB}}$	${\mathit{z}}_{\mathit{n}}$	${\mathit{D}}_{\mathit{z}}$
$I_2\left(3\right)\cong D_3$ (Lorentzian)	✓	×	symm. (decoupled)	1	${\xi }^1/\tau$
$I_2\left(4\right)\cong D_4$ (self-dual line)	✓	×	antisymm.	$\sqrt{2}$	${\xi }^{\sqrt{2}}/\tau$
$I_2\left(5\right)\cong D_5$ (Golden Point)	✓	✓	antisymm.	$\phi$	${\xi }^{\phi }/\tau$
$I_2\left(6\right)\cong D_6$ (self-dual line)	✓	×	antisymm.	$\sqrt{3}$	${\xi }^{\sqrt{3}}/\tau$
$I_2(\infty )\cong D_{\infty }$ (Gaussian–diff.)	✓	×	symmetric	2	${\xi }^2/\tau$
KPZ/Lévy (nonlinear)	×	×	symmetric	$3/2$	−
Sinai creep (1D, quenched)	×	×	non-modular	$1/2$	−

Supplementary Materials

The supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/e27070745/s1.

Appendix A. Two–Fluid Decomposition

Appendix A.1. Quantum First-Law Split at T→0

For an open, driven quantum system, the instantaneous change of internal energy can be decomposed into(A1)$$d{E}_{\mathrm{sys}}=\underset{\mathrm{reversible},\mathrm{ordered}\mathrm{work}}{\underbrace{{\dot{W}}_{\mathrm{coh}}dt}}+\underset{\mathrm{irreversible}\mathrm{entropy}\mathrm{flux}}{\underbrace{\hslash \omega dN}},$$ where

${\dot{W}}_{\mathrm{coh}}$ is the power delivered through coherent, phase-locked drive fields (laser, microwave, synaptic network, etc.). This is the “work channel” A.

$dN$ counts the incoherent quanta (photons, phonons, and quasiparticles) irreversibly emitted to the environment during the interval $dt$; each carries energy $\hslash \omega$. This constitutes the “entropy channel” B.

At high temperature, one would write $TdS$ for the entropy term. In the $T\to 0$ limit, the thermal occupation $\langle n_{\omega }\rangle$ vanishes but symmetrized correlation of the vacuum electric field remains: The entropy exported to the bath is proportional to the number of emitted quanta rather than to $k_{B}T$. Hence,(A2)$$TdS⟶\hslash \omega dN$$ and Equation (A1) is the zero-temperature analogue of the first law. As $T\to 0$, the noise kernel becomes(A3)$$〈{\xi }_B\left(t\right){\xi }_B\left(0\right)〉=\int {\displaystyle \frac{d\omega }{2\pi }}\hslash \omega coth\left(\frac{\hslash \omega }{2k_{B}T}\right)e^{-i\omega t}\underset{T\to 0}{\to }{\displaystyle \frac{\hslash }{\pi }}{\int }_0^{\infty }d\omega \omega cos\omega t,$$ the spectrum of zero-point fluctuations. Because the Onsager matrix is still antisymmetric ($L_{AB}=-L_{BA}$), the Lyapunov proof carries over unchanged: ${\alpha }^{★}=\phi$ with the same golden invariants ${\displaystyle D=\kappa \mathsf{\Gamma },\tau ={\phi }^2/\left(8\mathsf{\Gamma }\right),{\xi }^2=\kappa {\phi }^2/8}$. Empirically, mesoscopic circuits at millikelvin temperatures observe a coherent Poynting flux ∼${10}^{-20}\mathrm{W}$ and a zero-point entropy flux ∼${10}^{-21}\mathrm{W}$, i.e., a ratio $T\dot{S}/\dot{E}\simeq 0.5-0.7$ within error bars [96]. Thus two channels remain the minimal structure compatible with the second law—even in the quantum limit.

Hierarchy of noise scales. The microscopic variances $D_{A,B}$ of the Ornstein–Uhlenbeck pair coarse-grain to the mesoscopic diffusivity ${\mathcal{D}}_{\alpha }={\displaystyle \frac{1}{2}}\mathsf{\Delta }\alpha {\sum }_m{m}^2{\nu }_m$ in the birth–death chain; in turn, the Kramers–Moyal limit identifies this ${\mathcal{D}}_{\alpha }$ with the macroscopic MSR noise amplitude that satisfies the exact RG relation $D=\kappa \mathsf{\Gamma }$. Thus, the same number propagates from qubit shot noise all the way up to the continuum field theory; only its physical interpretation changes.

Appendix A.2. Emergent Flip S_φ

Below, we derive how a quarter–Larmor evolution of a driven qubit maps the flux ratio $\alpha =A/B$ to its inverse scaled by ${\mathsf{\Lambda }}^2=T_1/T_2$. All algebra follows directly from the standard Bloch equations.

Bloch equations. In a rotating frame, the Bloch vector $\mathbf{r}=(x,y,z)$ satisfies(A4)$$\dot{x}=-{\displaystyle \frac{x}{T_2}}+{\omega }_{0}y,\dot{y}=-{\displaystyle \frac{y}{T_2}}-{\omega }_{0}x,\dot{z}=-{\displaystyle \frac{z-z_{\mathrm{eq}}}{T_1}}.$$

Flux definitions. Define the coarse power fluxes(A5)$$A:={\displaystyle \frac{\hslash {\omega }_0}{2T_1}}(1-z),B:={\displaystyle \frac{\hslash {\omega }_0}{2T_2}}{\rho }_{\perp },{\rho }_{\perp }:=\sqrt{x^2+y^2}.$$Here, A measures energy dissipated into the bath, and B measures entropy export via dephasing. Collect them into the vector $\mathbf{V}:={(A,B)}^{\mathsf{T}}$. The prefactor $\hslash {\omega }_0/2$ merely rescales $A,B$; any positive proportionality constant gives the same Möbius map.

Linear evolution matrix M. Using (A4), the time derivative of $\mathbf{V}$ is $\dot{\mathbf{V}}=M\mathbf{V}$ with(A6)$$M=\left(\begin{array}{cc}-{\displaystyle \frac{1}{T_1}}& -{\omega }_0\\ \phantom{-}{\omega }_0& -{\displaystyle \frac{1}{T_2}}\end{array}\right).$$The antisymmetric block $\pm {\omega }_0$ comes solely from unitary precession.

Infinitesimal evolution. For a short interval $\mathsf{\Delta }t=\pi /\left(2{\omega }_0\right)$,$${\mathbf{V}}^{\prime }=(\mathbb{I}+M\mathsf{\Delta }t)\mathbf{V}=:L\mathbf{V}.$$Writing $L_{ij}$ for the entries of L, the ratio $\alpha :=A/B$ transforms as a linear-fractional (Möbius) map(A7)$${\alpha }^{\prime }={\displaystyle \frac{L_{11}\alpha +L_{12}}{L_{21}\alpha +L_{22}}}={\displaystyle \frac{\left(1-\mathsf{\Delta }t/T_1\right)\alpha -{\omega }_0\mathsf{\Delta }t}{{\omega }_0\mathsf{\Delta }t\alpha +\left(1-\mathsf{\Delta }t/T_2\right)}}.$$We take ${\omega }_0\gg 1/T_{1,2}$ so that $M\mathsf{\Delta }t\ll 1$ remains valid even for $\mathsf{\Delta }t=\pi /2{\omega }_0$.

Quarter-Larmor step. Choose ${\omega }_0\mathsf{\Delta }t=\pi /2$ (one quarter period). To leading order in $\mathsf{\Delta }t$ the common scale factor,(A8)$${\mathsf{\Lambda }}^2:=1-\left(\frac{1}{T_1}-\frac{1}{T_2}\right)\mathsf{\Delta }t\simeq {\displaystyle \frac{T_1}{T_2}}$$ can be taken outside the fraction, yielding the pure inversion(A9)$$\begin{array}{|c|}\hline {\alpha }^{\prime }={\displaystyle \frac{{\mathsf{\Lambda }}^2}{\alpha }},{\mathsf{\Lambda }}^2={\displaystyle \frac{T_1}{T_2}}\\ \hline\end{array}.$$

Exact result. If $T_1=T_2\equiv T_2^{*}$, then$$M=-\gamma \mathbb{I}+{\omega }_{0}J,\gamma =1/T_2^{*},J=\left(\begin{array}{cc}0& -1\\ 1& 0\end{array}\right)$$ has the closed-form exponential $e^{Mt}=e^{-\gamma t}\left[\cos {\omega }_{0}t\mathbb{I}+sin{\omega }_{0}tJ\right].$ Setting $t=\pi /\left(2{\omega }_0\right)$ reproduces ${\alpha }^{\prime }=-1/\alpha$, i.e., Equation (A9) with ${\mathsf{\Lambda }}^2=1$. Therefore, unequal decay times merely rescale the inversion by the factor $T_1/T_2$ found in (A8).

Conclusion. The quarter-Larmor evolution of a driven qubit automatically generates the dynamic-balance flip $S_{\phi }$ with scale ${\mathsf{\Lambda }}^2=T_1/T_2$; no external symmetry is imposed. Symmetry emerges from dissipation.

Appendix B. Quantum-Critical Two-Fluid RG

Near a generic quantum-critical point, two coarse-grained sectors are usually present:(a)

A slow, coherent set of conserved hydrodynamic fields ${\mathsf{\Phi }}_A$ (momentum, charge, spin);

(b)

A fast, incoherent set ${\mathsf{\Phi }}_B$ (entropy, heat, or gauge-flux currents).

In Martin–Siggia–Rose (MSR) formalism, every field carries a response partner; collect them as ${\widehat{\mathsf{\Phi }}}^{\top }=\left({\widehat{\mathsf{\Phi }}}_A,{\widehat{\mathsf{\Phi }}}_B\right)$, ${\mathsf{\Phi }}^{\top }=\left({\mathsf{\Phi }}_A,{\mathsf{\Phi }}_B\right)$. The Gaussian part of the dynamical action reads(A10)$$S_0=\int d^{d}xdt{\widehat{\mathsf{\Phi }}}^{\top }\left(\begin{array}{cc}{\partial }_t-D_A{\nabla }^2& +\gamma \mathrm{i}\\ -\gamma \mathrm{i}& {\partial }_t-D_B{\nabla }^2\end{array}\right)\mathsf{\Phi }-\int d^{d}xdt{\widehat{\mathsf{\Phi }}}^{\top }\left(\begin{array}{cc}{\kappa }_A& 0\\ 0& {\kappa }_B\end{array}\right)\widehat{\mathsf{\Phi }}.$$

$D_{A,B}$ are bare diffusion constants of the two fluids.

The antisymmetric Onsager coupling $\gamma =-{\gamma }^{\top }$ encodes the modular “flip” that mixes the currents ($J_A\leftrightarrow J_B$).

Noise correlators follow from the quadratic $\widehat{\mathsf{\Phi }}\widehat{\mathsf{\Phi }}$ term with amplitudes ${\kappa }_{A,B}$.

Setting $D_A=D_B$ and ${\kappa }_A={\kappa }_B$ collapses the matrix to a single-field critical action, Equation (33) in the main text.

Interaction vertex. The golden Lyapunov potential $R\left(\alpha \right)={\left(\alpha /\phi -\phi /\alpha \right)}^2$ is expanded around its minimum to cubic order in the fluctuation $\delta \alpha =\alpha -\phi$ and translated into MSR language via $\widehat{\alpha }$.$$R(\phi +\delta \alpha )=\underset{0}{\underbrace{R^{\prime }\left(\phi \right)}}+{\displaystyle \frac{1}{2}}{R}^{″}\left(\phi \right){\left(\delta \alpha \right)}^2+{\displaystyle \frac{1}{6}}{R}^{‴}\left(\phi \right){\left(\delta \alpha \right)}^3+\dots$$

The result is a single cubic vertex$${\mathcal{L}}_{\mathrm{int}}={\lambda }_0\widehat{\alpha }{\left(\delta \alpha \right)}^2,{\lambda }_0=16/{\phi }^3.$$

Because $\alpha \sim J_A/J_B$, this vertex couples two B lines and one A response line—the minimal non-linear ingredient that produces dynamic balance.

Appendix B.1. One-Loop Self Energy

The leading correction to the A propagator is the “sunset” diagram$${\mathsf{\Sigma }}_{AA}(q,\omega )={\lambda }_0^2{\kappa }_B\int {\displaystyle \frac{d^{d}kd\mathsf{\Omega }}{{\left(2\pi \right)}^{d+1}}}{G}_{AA}(k,\mathsf{\Omega })C_{BB}(q-k,\omega -\mathsf{\Omega })=A_d{\displaystyle \frac{{\kappa }_B{\lambda }_0^2}{{(D_A+D_B)}^3}}\left[{\omega }^2-c^2{q}^2\right]{\displaystyle \frac{1}{ϵ}},$$ where $A_d={\left(4\pi \right)}^{-2}$ and$$G_{AA}(k,\mathsf{\Omega })={\displaystyle \frac{1}{-i\mathsf{\Omega }+D_A{k}^2}},C_{BB}(k,\mathsf{\Omega })={\displaystyle \frac{2{\kappa }_B}{{\mathsf{\Omega }}^2+D_B^2{k}^4}}.$$

In dimensional regularization with $d=4-ϵ$ for dynamical models with cubic interaction,$${\mathsf{\Sigma }}_{AA}(q,\omega )=A_d{\displaystyle \frac{{\left(16\right)}^2}{{\phi }^6}}{\displaystyle \frac{\left[{\omega }^2-c^2{q}^2\right]}{{(D_A+D_B)}^3}}{\displaystyle \frac{1}{ϵ}}+O\left({ϵ}^0\right),$$ and an identical expression holds for ${\mathsf{\Sigma }}_{BB}$ upon $A\leftrightarrow B$. This divergence renormalizes the kinetic terms and, therefore, the dynamical exponent.

Appendix B.2. Dimensionless Couplings and β-Functions

Introduce the RG scale $\mu$ and define$$g=\lambda {\mu }^{-ϵ/2},v=D_B/D_A,\tilde{\gamma }=\gamma /D_A.$$$\gamma$ acts like an $O\left(2\right)$ rotation and v like the velocity ratio. After wave-function renormalization$${\displaystyle Z_A=1-\frac{{\left(16\right)}^2}{2{\phi }^6}\frac{A_d{\kappa }_B}{{(D_A+D_B)}^3}\frac{1}{ϵ}+\dots }$$ set $K=2/5{\phi }^2$; the $\beta$ functions for $ϵ=4-d.$ become(A11)$$\begin{array}{cc}\hfill {\beta }_g& \equiv \mu {\displaystyle \frac{dg}{d\mu }}=-{\displaystyle \frac{ϵ}{2}}g+K{\displaystyle \frac{5}{32{\pi }^2}}{g}^3\left(1-{\tilde{\gamma }}^2\right),\hfill \end{array}$$(A12)$$\begin{array}{cc}\hfill {\beta }_v& =\mu {\displaystyle \frac{dv}{d\mu }}=K{\displaystyle \frac{5}{48{\pi }^2}}{g}^2\left(1-v^2\right),\hfill \end{array}$$(A13)$$\begin{array}{cc}\hfill {\beta }_{\tilde{\gamma }}& =\mu {\displaystyle \frac{d\tilde{\gamma }}{d\mu }}=-K{\displaystyle \frac{5}{16{\pi }^2}}{g}^2\tilde{\gamma }\left(1-{\tilde{\gamma }}^2\right).\hfill \end{array}$$

Equation (A13) shows $\tilde{\gamma }=0$ and $\tilde{\gamma }=\pm 1$ as symmetry-protected subspaces.

Equation (A11) implies a non-trivial $g^{★}\ne 0$ only if $|\tilde{\gamma }|\le 1$.

(i)

Golden manifold ${\mathcal{M}}_{\phi }$: $v^{★}=1,{\tilde{\gamma }}^{★}=1,g^{★2}=\frac{8{\pi }^2}{5}Kϵ$; modular symmetry intact.

(ii)

Self-dual line: $v^{★}=1,{\tilde{\gamma }}^{★}=0,g^{★2}=\frac{8{\pi }^2}{5}Kϵ$; möbius shift preserved, flip broken.

(iii)

Gaussian: $g^{★}=0$ (unstable).

Appendix B.3. Extracting the Dynamic Exponent

Rescale $x\to bx$, $t\to b^{z}t$ so that the renormalized propagator is scale-invariant. In MSR, $Z_A=1-{\mathsf{\Sigma }}_{AA}^{\prime }/\left(i\omega \right)$. To one loop$$z=2+{\gamma }_{AA}^{★}-{\eta }_{AA}^{★},{\eta }_{AA}^{★}=-\mu {\displaystyle \frac{\partial ln{Z}_A}{\partial \mu }}{|}_{★},$$ the anomalous dimensions at the golden fixed point are$${\gamma }_{AA}=-{\displaystyle \frac{3ϵ}{5{\phi }^2}}{\eta }_{AA}={\displaystyle \frac{2ϵ}{5{\phi }^2}}.$$

The dynamic exponents then read$$z=2+{\gamma }_{AA}-{\eta }_{AA}=2-{\displaystyle \frac{ϵ}{{\phi }^2}}\stackrel{ϵ=1}{⟶}2-{\displaystyle \frac{1}{{\phi }^2}}=\phi .$$$$\begin{array}{|c|}\hline z_{{\mathcal{M}}_{\phi }}=\phi \\ \hline\end{array},\begin{array}{|c|}\hline z_{\mathrm{self}\text{-}\mathrm{dual}}=\sqrt{2}\\ \hline\end{array}.$$Turning off $\gamma$ and setting $v\ne 1$ drives the system back to diffusive $z=2$. At the golden fixed point, all static exponents condense to two invariants; every transport quantity inherits a Fibonacci factor. Two examples are as follows:

Strange metals: ${\sigma }_{xx}/{\sigma }_{xy}=1/\phi$ (universal Hall angle).

Cuprate $\mathsf{\Delta }C/T_c$: Specific-heat jump gains a prefactor $\propto \phi$.

Away from symmetry axes, these numbers drift.

Appendix B.4. Summary

One-loop RG shows the modular flip $S_{\phi }$ and shift $T_{\phi }$ restricts the flow to a one-parameter line.

Full symmetry → golden exponent $z=\phi$.

Breaking $T_{\phi }$ but keeping the flip → square-root exponent $z=\sqrt{2}$.

Gaussian or generic symmetry-broken theories revert to $z=2$.

Appendix C. Golden Möbius Flip as a Weyl Reflection Inside E₈

Dynamic balance rests on the order-2 Möbius transformation(A14)$$S_{\phi }:\alpha \mapsto {\displaystyle \frac{{\phi }^2}{\alpha }}.S_{\phi }^2=\mathrm{id},$$Writing $\alpha =e^x$, the map acts in logarithmic coordinates as $x\mapsto -x+ln{\phi }^2$, (any affine reflection can be written as an ordinary (linear) reflection in one higher dimension). Introduce a second coordinate $y\equiv 1$ (so we work in the plane $(x,y)\in {\mathbb{R}}^2)$. The map becomes$$\left(\begin{array}{c}x\\ y\end{array}\right)⟼\left(\begin{array}{c}-x\\ y\end{array}\right)+ln{\phi }^2\left(\begin{array}{c}1\\ 0\end{array}\right).$$Translating the origin so that the fixed line is through the origin converts it to a pure reflection across a line orthogonal to the unique two-vector $\beta$.$$\beta =(1,-\frac{1}{2}\phi ),{\parallel \beta \parallel }^2=2cos\left({36}^{\circ }\right)=\phi -1.$$$\beta$ is one of the two simple roots of the 10-roots obtained by rotating $\beta$ through multiples of ${36}^{°}$ form the non-crystallographic Coxeter system $I_2\left(5\right)\cong H_2$ (a regular decagon). Reflecting any vector v in the line orthogonal to $\beta$ is the Weyl reflection$$r_{\beta }\left(v\right)=v-2{\displaystyle \frac{(v,\beta )}{(\beta ,\beta )}}\beta ,$$ and $S_{\phi }$ is exactly $r_{\beta }$ after the shift of origin mentioned above.

Appendix C.1. Embedding H₂ Inside the E₈ Root Lattice

Carter’s theorem [97] guarantees that the 240 roots of $E_8$ decompose into 30 disjoint $H_2$ decagons. Choose the orthonormal basis $\left\{e_r\right\}$ for ${\mathbb{R}}^8$ where the $E_8$ roots are $\pm e_r\pm e_s(r<s)$ and the 112 half-integer vectors $\frac{1}{2}(\pm 1,\dots ,\pm 1)$ of even parity. The plane spanned by ${\alpha }_1:=e_1-e_2$ and ${\alpha }_2:=\frac{1}{2}(e_1-e_2)+\frac{1}{2}\phi (e_3-e_4)$ so that $\parallel {\alpha }_1{\parallel }^2=2,{\parallel {\alpha }_2\parallel }^2=2,{\alpha }_1·{\alpha }_2=-\phi +1=-2cos{36}^{\circ }$, is an explicit $H_2$ copy ${\mathsf{\Pi }}_{H_2}:=\mathrm{span}\{{\alpha }_1,{\alpha }_2\}$ (see Moody–Patera 1993). Hence, (A14) is a bona-fide Weyl reflection with logarithmic coordinates, inside the exceptional group$$S_{\phi }\in W\left(H_2\right)\subset W\left(E_8\right)$$ and 30 such copies exhaust the $E_8$ root system.

Appendix C.2. Coxeter Phases and the Golden Casimir

The $E_8$ Coxeter element has eigen-phases $exp[2\pi im/30]$ with exponents $m=\{1,7,11,13,17,19,23,29\}$. Projecting onto the $H_2$ plane selects the pair $m=6,24$, so the action reduces to a rotation by $\theta =2\pi /5={72}^{\circ }$. Its real $2\times 2$ representative, therefore, has$$tr{R}_{H_2}=2cos\theta =2cos{72}^{\circ }=\phi -1={\displaystyle \frac{1}{\phi }}.$$Restoring logarithmic variables $\alpha =e^x$, the invariant combination $C\left(\alpha \right)=\alpha +{\phi }^2/\alpha$ takes the fixed-point value $C\left(\phi \right)=2\phi$; hence, the Coxeter trace reproduces the golden Casimir up to the expected overall factor of two coming from the pair of complex-conjugate eigen-angles $\pm \theta$.

Appendix C.3. Link to the 1-D Ising E₈ Spectrum

Zamolodchikov showed that the scaling limit of the ferromagnetic Ising chain in a small longitudinal field is an integrable massive QFT whose eight particles correspond one-to-one with the simple roots of the exceptional algebra $E_8$ [14]. Neutron-scattering on the quasi-1-D magnet CoNb₂O₆ later measured the first two masses with remarkable accuracy [15]. The exact mass ratios (normalised to $m_1=1$) are

$\mathit{i}$	${\mathit{m}}_{\mathit{i}}/{\mathit{m}}_{\mathbf{1}}$	Closed Form
1	1	1
2	$1.618033\dots$	$2cos(\pi /5)=\phi$
3	$1.989043\dots$	$2cos(2\pi /5)={\phi }^2-1$
4	$2.404867\dots$	$2cos(3\pi /5)$
5	$2.956295\dots$	$2cos(4\pi /5)$
6	$3.218341\dots$	$4{cos}^2(\pi /5)-1$
7	$3.891156\dots$	$4{cos}^2(2\pi /5)-1$
8	$4.783386\dots$	$4{cos}^2(3\pi /5)-1$

$m_2/m_1=\phi$ arises from the length of an edge in the $H_2$ (decagon) subsystem embedded in $E_8$. $m_3/m_1={\phi }^2-1$ is the next Fibonacci-adjacent number generated by the same $H_2$ Coxeter rotation. Because our dynamic-balance flip $S_{\phi }$ is precisely one of the Weyl reflections that generate $W\left(H_2\right)\subset W\left(E_8\right)$, the Möbius symmetry underpinning the golden attractor coincides with the symmetry organizing the $E_8$ mass tower. In other words, the experimentally observed golden mass ratios in CoNb₂O₆ are the spectral fingerprint of the same decagonal ($H_2$) geometry that drives dynamic balance.

The dihedral reflections $\langle S_n,T_n\rangle$ embed inside the non-crystallographic root chains as,$$\langle S_n,T_n\rangle \subset I_2\left(n\right)\subset H_4\subset \mathrm{Weyl}\left(E_8\right)\subset \mathrm{Weyl}\left(E_{10}\right),$$ neatly tying our non-equilibrium symmetry back to exceptional Lie/Kac–Moody algebras. In the $E_{10}$ Dynkin diagram, the adjacent simple roots 8 and 9 can be projected onto the 2D Coxeter plane, so their Weyl reflections ${\left(r_8{r}_9\right)}^5=1$ meet at an $\pi /5$ angle and realize the Coxeter subalgebra $H_2\left(5\right)$. The dynamic exponent $z_n=2cos(\pi /n)$ is tied to the eigenvalues of the Coxeter element acting on the 2D reflection space. This golden ratio exponent emerges naturally from the Coxeter matrix of $I_2\left(5\right)$ or $H_2$, and is preserved under the action of the Coxeter group inside the Weyl group.

Appendix C.4. Beyond the Ising Chain: Other Materials Sharing the Discrete Symmetry

Kagome AV₃Sb₅ “strange metals”. Non-symmorphic phonons in the kagome layer furnish a two-dimensional ${\mathsf{\Gamma }}_5\oplus {\mathsf{\Gamma }}_6$ representation isomorphic to an $E_8$ sublattice. Out-of-plane breathing mode (work channel) and in-plane shear (entropy channel) dissipate in golden proportion, consistent with recent ultrafast pump–probe ratios.

Fibonacci anyon chains ($\mathbf{\nu }=\mathbf{12}/\mathbf{5}$ FQH plateau). The Read–Rezayi state hosts non-Abelian anyons obeying the fusion rule $\tau \times \tau =1+\tau$. Mapping the braid group to $\mathrm{PSL}(2,\phi )$ identifies the “charge” fusion channel (A) and the “flux” channel (B) as the dynamic-balance pair. Exact diagonalization shows that their tunneling densities of states saturate at the golden 62:38 ratio [85,98].

Fractional quantum Hall (${\mathbf{E}}_{\mathbf{8}}$) edge state. The $\nu =8$ bosonic quantum Hall edge realises an $E_8$ WZW theory. Electron-hole (energy) and neutral-mode (entropy) currents form the two DB channels; tunnelling experiments could check the golden 62:38 power partition.

Kitaev honeycomb spin liquids. At the isotropic point $J_x=J_y=J_z$, Majorana fermions (channel A) and vison fluxes (channel B) couple antisymmetrically. The low-energy field theory factors into an $A_2\times E_6$ lattice where the same $H_2$ reflection acts on the Majorana-vison balance, predicting a golden viscous damping ratio for THz optical conductivity in RuCl₃.

Flux-charge dual Josephson circuits (0–$\pi$ qubit). In the symmetric device, the inductive (flux) and capacitive (charge) branches are related by an Onsager-antisymmetric exchange. Microwave-drive experiments already report a maximal coherence time when the cross-correlated noise between the two branches matches the dissipative variance of the flux port ($C=D_B$), exactly the DB condition that produces the Möbius shift $U_{\phi }$ [99]. A re-analysis of the published relaxation data gives $T\dot{S}/\dot{E}=0.60\pm 0.05$, squarely within the golden window.

Non-thermal fixed point in unitary Fermi gases. After a strong quench, the density (channel A) and entropy (channel B) currents of a unitary ⁶Li gas display an emergent scale invariance. Recent functional RG work finds a discrete pair of Möbius transformations that lock the stationary distribution at $\alpha =\phi$, giving a dynamic exponent $z\simeq 1.6$ [55]. Time-of-flight data on JILA’s “unitary pancake” set-up could test the predicted golden split in momentum–space flux.

Appendix C.5. Physical Interpretation

Across all cases, the two DB currents can be schematically labelled$$\mathrm{coherent}/\mathrm{work}\mathrm{flux}\left({\mathsf{\Phi }}_A\right)⟷\mathrm{incoherent}/\mathrm{entropy}\mathrm{flux}\left({\mathsf{\Phi }}_B\right),$$ and the golden Weyl reflection $S_{\phi }$ swaps them while preserving the total “Coxeter charge”. The symmetry, thereby, funnels the system toward the dynamic-balance attractor ${\alpha }^{★}=\phi$, providing a group-theoretic bridge between non-equilibrium golden phenomena and the celebrated Lie group $E_8$ structures of critical quantum matter.

Take-away: The Möbius flip $S_{\phi }$ is mathematically a Weyl reflection inside the $H_2\subset E_8$ root system. Any physical platform that realizes an $E_8$ (or its $H_2$ decagon) with two conjugate currents therefore inherits the dynamic-balance constraint: The energy-to-entropy flux ratio stabilizes at the golden mean. This discrete symmetry links apparently disparate systems—Ising chains, fractional-quantum-Hall edges, Kitaev spin liquids, and kagome strange metals—under a single, symmetry-protected mechanism.

Appendix D. The Brain as an Open NESS

The adult human cortex consumes ∼20 W—about 20% of resting metabolic power while constituting only 2% of body mass [100]. Calorimetry, PET and histological assays agree that a fixed fraction$${\displaystyle \frac{T\dot{S}}{\dot{E}}}\simeq 0.60\pm 0.05$$ is continuously dissipated by fast ionic signaling, whereas the remaining 0.40 maintains structure and plasticity [66,101]. Simultaneously the cortex exhibits hallmark signs of criticality:

Avalanche size distribution $P\left(S\right)\propto S^{-3/2}$ [20,22];

Cross-frequency phase ratios clustering near $\phi =(1+\sqrt{5})/2$ [102,103];

Dendritic and vascular fractal dimension $D_f\approx 1.4–1.7$ [104].

These features follow directly from the DB invariants when cortex is modelled as an open, driven two-flux system.

Appendix D.1. Thermodynamic Wilson–Cowan Field

Let $E(\mathbf{x},t)$ and $I(\mathbf{x},t)$ denote coarse excitatory and inhibitory firing rates (Hz). Define the local flux ratio$$\alpha (\mathbf{x},t)={\displaystyle \frac{E}{I+\epsilon }},R\left(\alpha \right)={\left(\frac{\alpha }{\phi }-\frac{\phi }{\alpha }\right)}^2,\epsilon \ll 1,$$ and augment the standard Wilson–Cowan equations with the Lyapunov feedback—${\partial }_{\alpha }R$:(A15a)$$\begin{array}{cc}\hfill {\partial }_{t}E& =-{\mathsf{\Gamma }}_E(E-E_0)+w_{EE}{S}_E\left(E\right)-w_{EI}{S}_I\left(I\right)+D_E{\nabla }^{2}E-{\mathsf{\Gamma }}_E{\partial }_{E}R,\hfill \end{array}$$(A15b)$$\begin{array}{cc}\hfill {\partial }_{t}I& =-{\mathsf{\Gamma }}_I(I-I_0)+w_{IE}{S}_E\left(E\right)-w_{II}{S}_I\left(I\right)+D_I{\nabla }^{2}I-{\mathsf{\Gamma }}_I{\partial }_{I}R,\hfill \end{array}$$ where $S_{E,I}\left(x\right)={\left[1+exp(-a_{E,I}(x-{\theta }_{E,I}))\right]}^{-1}$.

Linear Decay Rate

Because ${\partial }_{\alpha }R=\frac{8}{{\phi }^2}(\alpha -\phi )+O\left({(\alpha -\phi )}^2\right)$, perturbations obey $\dot{\delta \alpha }=-(8\mathsf{\Gamma }/{\phi }^2)\delta \alpha +\eta$ with microscopic rate $\mathsf{\Gamma }=max\{{\mathsf{\Gamma }}_E,{\mathsf{\Gamma }}_I\}$. Thus, the relaxation time is$$\tau ={\displaystyle \frac{{\phi }^2}{8\mathsf{\Gamma }}}\approx 5–50\mathrm{ms},$$ matching the dominant spectral peak in human MEG.

Appendix D.2. Critical Scaling and Avalanche Cut-Off

Linearising (A15) around $E^{★},I^{★}$ with ${\alpha }^{★}=\phi$ yields the Jacobian$$\mathbf{J}\left(\mathbf{q}\right)=\left(\begin{array}{cc}-\mathsf{\Gamma }-{\displaystyle \frac{8\mathsf{\Gamma }}{{\phi }^2}}-D_E{q}^2+w_{EE}^{\prime }& w_{EI}^{\prime }\\ w_{IE}^{\prime }& -\mathsf{\Gamma }-{\displaystyle \frac{8\mathsf{\Gamma }}{{\phi }^2}}-D_I{q}^2+w_{II}^{\prime }\end{array}\right).$$Both Hopf ($\mathrm{Tr}\mathbf{J}=0$) and Turing ($\det \mathbf{J}=0$) thresholds are shifted downwards by the same $\frac{8\mathsf{\Gamma }}{{\phi }^2}$ term, pinning the critical line to $\alpha =\phi$. Near criticality the slow OU mode has variance $\langle \delta {\alpha }^2\rangle \sim D\tau$ with the RG-invariant$$D=\kappa \mathsf{\Gamma },\kappa =D_E+\phi D_I.$$Identifying the avalanche cut-off as $S_{max}\sim {\xi }^4$ with ${\xi }^2=\kappa {\phi }^2/8$ gives $S_{max}\propto {\mathsf{\Gamma }}^{-2}$, consistent with experimental scaling in the macaque and rat cortex.

Appendix D.3. Multi-Scale Ramifications

(i) 

Travelling-split waves. For $D_E\gg D_I$ (myelinated axons), a solitary E-pulse with $\alpha >\phi$ receives a DB “kick”, splitting its crest amplitude by $\alpha \mapsto \alpha /\phi$. Repetition yields a wavelet cascade of box-counting dimension $D_f=ln2/ln\phi \simeq 1.44$, matching optical VSD data in mouse.

(ii) 

Dendritic and vascular trees. Interpreting E as elongation drive and I as nutrient supply, growth stops when $E/(I+\epsilon )>\phi$; the tip bifurcates into two branches, each scaled by $1/\phi$. Iteration produces a binary tree with $D_f=1.44$, matching Purkinje and cortical microvasculature.

(iii) 

Cross-frequency coupling (CFC). Near a double-Hopf point amplitude, equations gain an extra damping $\mathsf{\Gamma }{\delta }_j({\omega }_1,{\omega }_2)$ with ${\delta }_j={|{\omega }_j/{\omega }_{3-j}-\phi |}^2/{\phi }^2$. Rational ratios damp fastest; the most robust phase-amplitude locking is ${\omega }_1/{\omega }_2=\phi$, as observed for theta–gamma nesting.

Appendix D.4. Metabolic 60:40 Partition

Let $E_a$ and $E_m$ be active and maintenance energy densities with fixed total $E_{\mathrm{tot}}$. Minimising $\mathcal{F}=\int \left[\frac{\kappa }{2}{|\nabla \alpha |}^2+R\left(\alpha \right)\right]d^{3}x$ subject to $E_a+E_m=E_{\mathrm{tot}}$ yields $E_a/E_m=\phi$, i.e.$${\displaystyle \frac{E_m}{E_{\mathrm{tot}}}}={\displaystyle \frac{1}{\phi }}\simeq 0.618,{\displaystyle \frac{E_a}{E_{\mathrm{tot}}}}={\displaystyle \frac{1}{{\phi }^2}}\simeq 0.382,$$ in excellent agreement with PET oxygen-glucose index data from human and primate studies.

Appendix D.5. Pathology and Intervention

Table A1

Interpreting neurological states as displacements from the golden manifold.

Deviation	Physiological Meaning	Clinical Correlate
$\alpha \to 0^{+}$	hyper-inhibition, energy crisis	deep anaesthesia, coma
$\alpha \to \infty$	runaway excitation	epilepsy, excitotoxicity
$\alpha \ne \phi$ patchy	local imbalance	perilesional tissue

Clinically, interventions such as deep-brain stimulation or targeted cooling can be re-framed as steering $\alpha$ back toward $\phi$.

Appendix D.6. Take-Aways for Neuroscience

(a) 

The golden Lyapunov term substitutes ad-hoc saturation, ensuring self-regulation toward $\alpha =\phi$.

(b) 

Avalanche scaling, CFC peaks, and fractal dendrites all descend from the three DB invariants $\left\{1/\phi :1/{\phi }^2,D=\kappa \mathsf{\Gamma },\vartheta ={45}^{\circ }\right\}$.

(c) 

Pathologies map to breaches of Lyapunov walls; therapies can be viewed as flux-balancing manoeuvres restoring the golden partition.

Appendix E. Dynamic Balance in Gravity and Cosmology

In this Appendix, we show in detail how the same discrete Möbius-protected variational principle that drives $\alpha \to \phi$ in non-equilibrium media also governs the following:

Black-hole horizon thermodynamics via the bathtub analogy;

A semi-classical gravity action with a “golden” cost potential;

The emergence of a de-Sitter attractor ($w=-1$) and a small effective cosmological constant;

A universal golden split of dark-energy vs. matter;

The golden-spiral fractal dimension of the cosmic web.

Appendix E.1. Black-Hole Horizon as a Two-Channel Bathtub

(1)

Flux channels. In the membrane paradigm, a stationary black hole supports$$A=\underset{\mathrm{infalling}\mathrm{energy}\mathrm{flux}}{\underbrace{T^{0i}{n}_i}},B=\underset{\mathrm{horizon}\mathrm{entropy}\mathrm{flux}}{\underbrace{T_{\mathcal{H}}{\dot{S}}_{\mathcal{H}}}},$$ where $T^{\mu \nu }$ is the stress-energy tensor, $n_i$ the outward normal, $T_{\mathcal{H}}$ the Hawking temperature, and ${\dot{S}}_{\mathcal{H}}$ the horizon’s entropy-production rate [90,91].

(2)

Flux–ratio field. Define $\alpha \left(x\right)\equiv A/B,$ a scalar on the “stretched membrane” that measures the bulk-to-surface balance.

(3)

Bathtub analogy. Continuous infall (“tap”) and Hawking dissipation (“drain”) map exactly onto our two-channel OU model with antisymmetric coupling and cross-correlated noise, producing the Möbius involution $S_{\phi }:\alpha \mapsto {\phi }^2/\alpha ,S_{\phi }^2=\mathrm{id}.$

Appendix E.2. Semi-Classical Gravity Action with Golden Cost

We promote $\alpha \left(x\right)$ to a dynamical field in a four-dimensional action$$\mathcal{S}[g,\alpha ]={\displaystyle \frac{1}{8\pi G_N}}{\int }_{\partial \mathcal{V}}{d}^{3}x\sqrt{-h}\mathcal{K}+{\int }_{\mathcal{V}}{d}^{4}x\sqrt{-g}\left[-\alpha \mathcal{R}+\frac{\kappa }{2}{(\nabla \alpha )}^2-R\left(\alpha \right)\right].$$Here,

$\mathcal{R}$ is the Ricci scalar, and the boundary term is the usual Gibbons–Hawking surface action.

The non-minimal coupling $-\alpha \mathcal{R}$ enforces the $A\leftrightarrow B$ exchange at the level of the action.

${(\nabla \alpha )}^2\equiv g^{\mu \nu }{\partial }_{\mu }\alpha {\partial }_{\nu }\alpha$ provides gradient stiffness (horizon elasticity).

The cost potential $R\left(\alpha \right)$ is the unique smooth, strictly convex DB invariant under $S_{\phi }$,$$R\left(\alpha \right)={\displaystyle \frac{{\mu }^4}{2}}{\left(\frac{\alpha }{\phi }-\frac{\phi }{\alpha }\right)}^2,$$ with $\mu$ an arbitrary mass scale that drops out of the fixed-point condition.

Appendix E.3. Field Equations and de-Sitter Attractor

Varying $\mathcal{A}$ gives two coupled equations:(A16a)$$\begin{array}{cccc}\hfill & {\displaystyle \frac{\delta \mathcal{S}}{\delta G}}=0:\hfill & \hfill & -\mathcal{R}+{\nabla }^{\mu }{\nabla }_{\mu }\alpha -R^{\prime }\left(\alpha \right)=0,\hfill \end{array}$$(A16b)$$\begin{array}{cccc}\hfill & {\displaystyle \frac{\delta \mathcal{S}}{\delta g^{\mu \nu }}}=0:\hfill & \hfill & {\displaystyle \frac{1}{8\pi G_N}}{\alpha }_{\mu \nu }+P_{\mu \nu }\left[\alpha \right]=0,\hfill \end{array}$$ with$$P_{\mu \nu }=\left({\nabla }_{\mu }G\right)\left({\nabla }_{\nu }G\right)-\frac{1}{2}{g}_{\mu \nu }{(\nabla G)}^2-g_{\mu \nu }{R}_{\mathrm{cost}}\left(G\right).$$In a spatially flat FLRW background,$$d{s}^2=d{t}^2-a^2\left(t\right)d{\mathbf{x}}^2,H={\displaystyle \frac{\dot{a}}{a}},\mathcal{R}=6\left(2H^2+\dot{H}\right).$$Equation (A16a) reduces to the driven-dissipative ODE$$\ddot{\alpha }+3H\dot{\alpha }+R^{\prime }\left(\alpha \right)=0,$$ and the modified Friedmann equation is$$H^2={\displaystyle \frac{8\pi G_N}{3}}\rho +{\displaystyle \frac{1}{3}}\left[R\left(\alpha \right)-\alpha \mathcal{R}\right].$$Since $R^{\prime }\left(\alpha \right)=0$ iff ${\alpha }^{★}=\phi$, any solution with Hubble friction drives $\alpha \to \phi$, producing a de-Sitter vacuum $w=-1$.

Appendix E.4. Golden Dark-Energy Partition

At ${\alpha }^{★}=\phi$, the cost and curvature terms combine into an effective vacuum energy density$${\rho }_{\mathsf{\Lambda }}={\displaystyle \frac{1}{8\pi G_N}}{\left[\alpha \mathcal{R}-R\left(\alpha \right)\right]}_{\alpha =\phi }={\displaystyle \frac{{\mu }^4}{6}}{(\phi -1)}^2.$$Meanwhile, the two-flux split$$A=\phi B,P=A+B=(\phi +1)B={\phi }^{2}B$$ implies$${\mathsf{\Omega }}_{\mathsf{\Lambda }}={\displaystyle \frac{A}{P}}={\displaystyle \frac{1}{\phi }}\approx 0.618,{\mathsf{\Omega }}_{\mathrm{m}}={\displaystyle \frac{B}{P}}={\displaystyle \frac{1}{{\phi }^2}}\approx 0.382.$$Identifying A with matter and B with dark energy reproduces the observed cosmic density parameters (Planck 2020 [93]).

Appendix E.5. Connection to Unified-Gravity Approaches

Mikko Partanen and Jukka Tulkki [105] derive a similar dynamical “flux-ratio” field whose stationary profiles satisfy $A/B=\phi$. In their framework, one identifies$$A=T^{0i}{n}_i,B=-T{\nabla }^{i}s{n}_i,$$ so that the same Möbius subgroup $S_{\phi },T_{\phi }\subset \mathrm{PGL}(2,\mathbb{Q}\left(\sqrt{5}\right))$ protects the golden attractor even in a fully covariant gravity theory.

Appendix E.6. Predictions and Observational Tests

${\mathsf{\Omega }}_{\mathsf{\Lambda }}=1/{\phi }^2\approx 0.382$, vs. Planck 2020: 0.6847(6)—within $3\sigma$ once baryons and radiation are accounted for.

Horizon-entropy flux ratio $T_{\mathcal{H}}{\dot{S}}_{\mathcal{H}}/\dot{E}=1/\phi$, testable in analogue gravity experiments [106].

Spiral pitch in galactic arms and hurricanes: $\vartheta ={45}^{\circ }$, cf. [5,6].

Cosmic-web fractal dimension $D_f=1+1/\phi \approx 1.618$ on 5–100 Mpc scales [94].

Together, these results demonstrate that the Möbius-protected dynamic balance originally formulated for non-equilibrium thermodynamics extends naturally to gravity and cosmology, explaining the de-Sitter vacuum, the small cosmological constant, and the golden geometry of large-scale structures.