Stabilizer Thermodynamic Systems

• Stabilizer thermodynamic systems are defined by enforcing stability through the curvature of energy potentials, concave entropy, and geometric criteria at equilibrium.
• They integrate dynamic stabilization via nonnegative entropy production and Lyapunov functionals to ensure robust hydrodynamic and control responses.
• These systems extend to quantum many-body settings and memory architectures, highlighting diverse applications across theoretical and applied thermodynamics.

Searching arXiv for the supplied topic and cited papers to ground the synthesis. “Stabilizer thermodynamic systems” is not a single uniformly defined doctrine. In the literature represented here, the phrase is used for several closely related constructions in which stability is enforced by thermodynamic structure itself: by the curvature of equilibrium potentials and their Hessians; by concave entropy and nonnegative entropy production in nonequilibrium continuum theories; by Lyapunov or availability functionals for open systems; by geometric criteria for invertibility of thermodynamic maps; and, in quantum many-body settings, by stabilizer-based observables whose scaling and phase structure are thermodynamic-like. Taken together, these uses identify a common theme: stability is not treated as an external add-on, but as a consequence of the admissible thermodynamic geometry, constitutive laws, or stabilizer structure (Lima et al., 2019, Somogyfoki et al., 2024, Schaft, 2021).

1. Equilibrium stabilization by curvature of thermodynamic potentials

In the equilibrium setting, the most explicit meaning of a stabilizer is given by the curvature of thermodynamic potentials near equilibrium. For a sufficiently smooth potential $f(x)$, Taylor expansion about an equilibrium state $x_0$ gives a quadratic form controlled by the Hessian. For internal energy per mole $u(s,v)$, with $s$ entropy and $v$ volume, the paper writes

$$\tilde{u}(s,v)\equiv 2![u(s,v)-u(s_0,v_0)]\cong a(\Delta s)^2+b(\Delta v)^2+c(\Delta s)(\Delta v) = \Delta x^{T}M\Delta x,$$

where $M$ is the symmetric Hessian-like matrix of second derivatives. Diagonalization yields

$$\tilde{u}(s,v)=\lambda_1 \Delta s'^2+\lambda_2 \Delta v'^2,$$

so the eigenvalues $\lambda_1,\lambda_2$ are the curvatures along principal directions. In this framework, the stabilizer of the thermodynamic system is precisely this local curvature data: second derivatives, Hessian definiteness, and eigenvalue signs (Lima et al., 2019).

For internal energy, equilibrium under the relevant constraints is a minimum, so the quadratic form must be positive semidefinite. In two variables this gives the standard definiteness conditions

$a>0,\qquad b>0,\qquad ab-c^2>0,$

or, equivalently,

$x_0$0

Entropy $x_0$1 obeys the dual statement: equilibrium is a maximum, so the Hessian is negative semidefinite, and $x_0$2 is concave in its extensive variables. Helmholtz free energy $x_0$3 and enthalpy $x_0$4 acquire mixed curvature and are saddle surfaces, while Gibbs free energy $x_0$5 is concave in its intensive variables (Lima et al., 2019).

These sign conditions are not merely geometric. They generate the usual thermodynamic stability inequalities. The same analysis yields

$x_0$6

and, through the standard relations

$x_0$7

also $x_0$8 and $x_0$9. A common misconception is that heat capacities and compressibilities are primary axioms of stability; in this formulation they are consequences of local curvature conditions on the relevant potential (Lima et al., 2019).

2. Thermodynamic laws as dynamical stabilizers

A second major usage treats thermodynamic structure itself as the stabilizer of dynamics. In extended heat conduction, the state space is enlarged from the equilibrium variable $u(s,v)$0 to $u(s,v)$1, where $u(s,v)$2 is the heat flux and $u(s,v)$3 a second-order tensor interpreted as flux of heat flux. The entropy density is chosen as

$u(s,v)$4

with $u(s,v)$5 and $u(s,v)$6 positive definite. The Hessian of $u(s,v)$7 is then negative semidefinite, so concavity of entropy directly encodes thermodynamic stability, while nonnegative entropy production constrains admissible transport coefficients. In one dimension, linear constitutive closure produces Onsager-type relations with coefficients $u(s,v)$8 and $u(s,v)$9, and the second law imposes positivity conditions such as

$s$0

together with the analogous conditions for $s$1. Linearization about homogeneous equilibrium leads to a cubic dispersion relation, and Routh–Hurwitz analysis shows that the entropy structure and entropy production almost entirely guarantee linear stability; only the extra coupling inequality

$s$2

must be added in the most general case (Somogyfoki et al., 2024).

The same principle appears at the level of general hydrodynamics. For hydrodynamic variables $s$3 and thermodynamic sources $s$4, linear perturbations obey

$s$5

The first law yields a Hermitian positive-definite susceptibility matrix $s$6, while the second law constrains the dissipative part of the constitutive matrix $s$7 to have positive symmetric part. The paper’s conclusion is that these two inputs are sufficient to place all hydrodynamic poles in the lower half-plane, so linear dynamical stability is not an independent assumption but a consequence of the thermodynamic laws. This statement is extended to systems with spontaneously or softly broken symmetries and in the presence of magnetic fields (Goutéraux et al., 2024).

Taken together, these results sharpen the equilibrium picture. Curvature stabilizes static states, while concave entropy plus nonnegative entropy production stabilize linear dynamics. A plausible implication is that equilibrium and hydrodynamic stability are different levels of the same structural constraint rather than separate criteria.

3. Lyapunov, port-thermodynamic, and control formulations

A third line of work develops explicit stabilizing functionals and control architectures from thermodynamic quantities themselves. For heat conduction in a rigid body,

$s$8

the usual “energy method” employs the artificial quadratic functional $s$9, but the thermodynamic construction replaces it by a physically motivated Lyapunov functional. For an equilibrium rest state $v$0, one obtains

$v$1

and for an inhomogeneous nonequilibrium steady state $v$2,

$v$3

Its time derivative satisfies

$v$4

so the functional is a strict Lyapunov stabilizer for the open-system steady state (Bulíček et al., 2017).

Port-thermodynamic systems generalize this logic geometrically. In that framework, thermodynamic state properties are encoded by a Liouville submanifold $v$5, and dynamics by a homogeneous Hamiltonian $v$6. Power-conserving interconnection via power ports produces another port-thermodynamic system, and the same holds for rate-of-entropy-increasing interconnection via entropy flow ports. Stability is then studied by Lyapunov functions built from generating functions for the submanifold and from additional conserved quantities, with canonical point transformations on the symplectized thermodynamic phase space serving as the shaping mechanism (Schaft, 2021).

A control-theoretic realization appears for disturbed continuous stirred tank reactors. In the sidSPHS formulation,

$v$7

the opposite entropy function is used as the modeling Hamiltonian, while the availability function is obtained by a generalized canonical transformation and then used for stabilization. Under a stochastic passivity condition, the feedback law

$v$8

asymptotically stabilizes the equilibrium in probability (Lu et al., 2017).

The same stabilizing role of thermal variables appears in abstract thermo-elasticity. For the Cattaneo-law thermo-elastic system, the thermal part acts as damping, with energy identity

$v$9

Exponential stability is characterized by an observability inequality, and exponential or polynomial stability for the Cattaneo system implies polynomial stability for the corresponding Fourier-law thermoelastic system (Hassi et al., 2013).

4. Geometric, driven, and ensemble-based stability limits

A distinct literature studies stability as a property of thermodynamic maps, driven free energies, or admissible ensembles. For classical discrete systems, the forward map from potential energy surface $$\tilde{u}(s,v)\equiv 2![u(s,v)-u(s_0,v_0)]\cong a(\Delta s)^2+b(\Delta v)^2+c(\Delta s)(\Delta v) = \Delta x^{T}M\Delta x,$$0 to equilibrium structure $$\tilde{u}(s,v)\equiv 2![u(s,v)-u(s_0,v_0)]\cong a(\Delta s)^2+b(\Delta v)^2+c(\Delta s)(\Delta v) = \Delta x^{T}M\Delta x,$$1 is the thermodynamic average $$\tilde{u}(s,v)\equiv 2![u(s,v)-u(s_0,v_0)]\cong a(\Delta s)^2+b(\Delta v)^2+c(\Delta s)(\Delta v) = \Delta x^{T}M\Delta x,$$2. Bidirectional stability asks whether the inverse problem $$\tilde{u}(s,v)\equiv 2![u(s,v)-u(s_0,v_0)]\cong a(\Delta s)^2+b(\Delta v)^2+c(\Delta s)(\Delta v) = \Delta x^{T}M\Delta x,$$3 is also well-conditioned. The paper introduces the anharmonicity vector field $$\tilde{u}(s,v)\equiv 2![u(s,v)-u(s_0,v_0)]\cong a(\Delta s)^2+b(\Delta v)^2+c(\Delta s)(\Delta v) = \Delta x^{T}M\Delta x,$$4, derived entirely from non-interacting configurational geometry, and shows that the local breaking of bidirectional stability satisfies

$$\tilde{u}(s,v)\equiv 2![u(s,v)-u(s_0,v_0)]\cong a(\Delta s)^2+b(\Delta v)^2+c(\Delta s)(\Delta v) = \Delta x^{T}M\Delta x,$$5

for $$\tilde{u}(s,v)\equiv 2![u(s,v)-u(s_0,v_0)]\cong a(\Delta s)^2+b(\Delta v)^2+c(\Delta s)(\Delta v) = \Delta x^{T}M\Delta x,$$6, with a sum over principal minors in higher dimensions. In this usage, a stabilizer thermodynamic system is one for which $$\tilde{u}(s,v)\equiv 2![u(s,v)-u(s_0,v_0)]\cong a(\Delta s)^2+b(\Delta v)^2+c(\Delta s)(\Delta v) = \Delta x^{T}M\Delta x,$$7, its divergence, and the relevant Jacobian terms remain small, so the thermodynamic map is nearly linear and robustly invertible (Yuge et al., 2018).

Driven reactive mixtures provide another stabilization mechanism. For a generalized Cahn–Hilliard–Allen–Cahn system, the relevant nonequilibrium potential is a driven Gibbs free energy $$\tilde{u}(s,v)\equiv 2![u(s,v)-u(s_0,v_0)]\cong a(\Delta s)^2+b(\Delta v)^2+c(\Delta s)(\Delta v) = \Delta x^{T}M\Delta x,$$8 whose second variation satisfies

$$\tilde{u}(s,v)\equiv 2![u(s,v)-u(s_0,v_0)]\cong a(\Delta s)^2+b(\Delta v)^2+c(\Delta s)(\Delta v) = \Delta x^{T}M\Delta x,$$9

Linear stability is governed by the competition of the chemical diffusion tensor $M$0 and the autocatalytic rate tensor $M$1, while electrochemical control enters through solo-autocatalysis and differential reaction resistance. In the one-species case, the growth rate takes the form

$M$2

and the theory predicts that spinodal decomposition can be suppressed above a critical current $M$3. This is the paper’s “electrochemical freezing” of phase separation, with the converse “electrochemical melting” under reversed conditions (Bazant, 2017).

Stability can also fail because the choice of control variables is itself destabilizing. For a binary ideal gas in the $M$4-ensemble, equilibrium requires

$M$5

When $M$6, a finite equilibrium exists. At the threshold $M$7, $M$8, and for $M$9 no physical finite solution exists. The paper therefore gives a counterexample to the widely accepted claim that any equilibrium state of a two-component system is determined by specifying four thermodynamic variables including at least one extensive variable. The controversy is not about Gibbs–Duhem counting itself, but about the stronger interpretation that any such choice necessarily leads to a stable equilibrium (Norizoe et al., 2018).

5. Quantum stabilizer observables and hidden thermodynamic structure

In several quantum many-body problems, “stabilizer thermodynamic systems” refers not to classical Hessians or Lyapunov functions, but to stabilizer-derived observables that scale extensively and diagnose phase structure invisible to ordinary thermodynamics. In the coupled Sachdev–Ye–Kitaev model, the stabilizer Rényi entropy

$$\tilde{u}(s,v)=\lambda_1 \Delta s'^2+\lambda_2 \Delta v'^2,$$0

is treated as a thermodynamic-like probe of quantum magic. A large-$$\tilde{u}(s,v)=\lambda_1 \Delta s'^2+\lambda_2 \Delta v'^2,$$1 path-integral and saddle-point analysis shows that $$\tilde{u}(s,v)=\lambda_1 \Delta s'^2+\lambda_2 \Delta v'^2,$$2 is extensive and exhibits first-order transitions as temperature varies. One of these, at $$\tilde{u}(s,v)=\lambda_1 \Delta s'^2+\lambda_2 \Delta v'^2,$$3 for $$\tilde{u}(s,v)=\lambda_1 \Delta s'^2+\lambda_2 \Delta v'^2,$$4, is an intrinsic SRE transition not visible in conventional thermodynamic quantities; it corresponds to a change in replica connectivity geometry rather than in the ordinary thermal saddle (Zhang et al., 22 Sep 2025).

For the open critical transverse-field Ising chain, the finite-temperature stabilizer Rényi entropy at Rényi index $$\tilde{u}(s,v)=\lambda_1 \Delta s'^2+\lambda_2 \Delta v'^2,$$5 is controlled by the Pauli-spectrum moment

$$\tilde{u}(s,v)=\lambda_1 \Delta s'^2+\lambda_2 \Delta v'^2,$$6

At $$\tilde{u}(s,v)=\lambda_1 \Delta s'^2+\lambda_2 \Delta v'^2,$$7, the resulting exponentially large sum over absolute values of all square minors of the correlation matrix is reduced exactly to a single Pfaffian. In the crossover window $$\tilde{u}(s,v)=\lambda_1 \Delta s'^2+\lambda_2 \Delta v'^2,$$8, the stabilizer entropy factorizes into a saturated extensive contribution and a universal finite-size scaling function, and that scaling function is a level-eight eta quotient rather than the ordinary free-boundary Majorana thermal factor. The paper’s interpretation is that finite-temperature stabilizer entropy reveals hidden defect-like conformal boundary data invisible to ordinary thermodynamic probes (Khasseh et al., 7 Jun 2026).

A related metrological usage appears in stabilizer codes under monitoring. There, nonlocal observables constructed from dual Ising spins have quantum Fisher information

$$\tilde{u}(s,v)=\lambda_1 \Delta s'^2+\lambda_2 \Delta v'^2,$$9

and the monitored dynamics exhibit transitions from an extensive regime, where long-range string order prevails, to an intensive regime driven by competing single-site measurements. In this sense, stabilizer structure supports a thermodynamic-like macroscopic resource whose scaling changes across measurement-induced phases (Lira-Solanilla et al., 16 Apr 2026).

The limitations of stabilizer thermodynamic behavior are also sharply characterized. For stabilizer zero-energy eigenstates of two-body Hamiltonians at infinite temperature, the paper proves a no-go theorem: such states cannot satisfy $\lambda_1,\lambda_2$0-body microscopic thermal equilibrium for any $\lambda_1,\lambda_2$1. The bound is tight, since explicit two-body nonintegrable Hamiltonians are constructed whose stabilizer eigenstates reproduce thermal expectation values for all two-body and all three-body observables. The structural origin is that any parent Hamiltonian must be built from factorizations of stabilizers into few-body Pauli strings, which bounds the achievable MITE order (Hokkyo, 22 Jan 2026).

6. Memory stabilizer structures and finite-size thermodynamics

The oldest usage in the supplied literature concerns stabilizer Hamiltonians for classical memory. A concatenated classical triple modular redundancy code is realized as an Ising-like stabilizer system, and a sequence of CNOT gates maps the interacting Hamiltonian into uncoupled single spins and, in some cases, uncoupled pairs. This gives an exact partition function, closed-form relative magnetization, and susceptibility for several “memory stabilizer structures” and for the canonical stabilizer Hamiltonian. For the pure Ising tree, the transformed Hamiltonian is

$\lambda_1,\lambda_2$2

and the relative magnetization obeys

$\lambda_1,\lambda_2$3

The analytical solutions show that all of the memory structures have no finite critical temperature in the thermodynamic limit,

$\lambda_1,\lambda_2$4

even though finite systems exhibit spontaneous magnetization below an apparent finite critical temperature. The mismatch arises because the infinite system is a poor approximation even for astronomically large finite systems. An additional, and somewhat counterintuitive, conclusion is that Hamiltonians with two-body interactions have a higher apparent critical temperature than the many-body canonical stabilizer Hamiltonian (Tomita et al., 2010).

This result clarifies a broader misconception relevant to stabilizer thermodynamic systems. Finite-size thermodynamic signatures—susceptibility peaks, apparent ordering temperatures, robust relative magnetization—need not imply a genuine finite-$\lambda_1,\lambda_2$5 ordered phase in the strict thermodynamic limit. In the stabilizer-memory setting, the practical finite-size system can be strongly ordered while the exact infinite-size critical temperature remains zero (Tomita et al., 2010).

Across these literatures, the phrase “stabilizer thermodynamic systems” therefore denotes a family of structurally constrained systems in which stability is encoded by curvature, entropy, dissipation, interconnection, geometric invertibility, or stabilizer observables. The common content is rigorous rather than terminological: stability is read off from Hessians, susceptibilities, entropy production matrices, Lyapunov or availability functionals, interconnection laws, or stabilizer-based scaling functions, and instability appears when those structures lose definiteness, invertibility, or dissipativity.