Excess Entropy Production as a Candidate Universal Cost of Persistence: A Thermodynamic Foundation for the Attractor Framework; Robert Galida (July 2026) [F]

Abstract

Every dissipative system maintains its attractor through continuous reconfiguration. Reconfiguration requires work; work generates entropy. The recovery rate κκ — corrective permeability — is the rate at which a system reconfigures to return to its attractor after perturbation. This paper proposes that κκ is a measure of excess entropy generation rate.

We develop an abstract persistence cost framework and prove its equivalence to Lyapunov theory. We then identify entropy production as a physical realization of this cost, deriving:κ=infxδ(x)0σexcess(ϕt(x))dtκ=xinf​∫0∞​σexcess​(ϕt​(x))dtδ(x)​

where σexcess=σσssσexcess​=σσss​ is the excess entropy production rate above the system’s steady-state baseline. For physical systems, the baseline is zero (equilibrium); for biological, cognitive, and social systems, the baseline is the steady-state dissipation rate of the healthy, well-coordinated attractor.

This unifies physical, biological, cognitive, and social systems. The framework is grounded in the second law of thermodynamics and non-equilibrium steady-state thermodynamics, not analogy. Empirical predictions are provided for each domain.

Keywords: entropy generation, excess entropy production, corrective permeability, attractor framework, dissipative structures, reconfiguration, Lyapunov theory, free energy principle, allostatic load


1. Introduction

The attractor framework defines persistence as the ability of a system to maintain its attractor under perturbation. Historically, persistence has been measured kinematically — as distance traveled or time spent away from equilibrium. This paper proposes that the true cost of persistence is thermodynamic: it is the excess entropy generated during reconfiguration and recovery.

Every dissipative system maintains its attractor through continuous reconfiguration. A bacterium reconfigures its metabolism to maintain homeostasis. A brain reconfigures its synaptic connections to maintain predictive models. A society reconfigures its institutions to maintain order. Reconfiguration requires work; work generates entropy. The second law of thermodynamics applies at every level of organization.

We develop an abstract persistence cost framework first, establishing its equivalence to Lyapunov theory. We then identify entropy production as a physical realization of this cost, deriving the relationship between corrective permeability and excess entropy generation.

The framework unifies physical, biological, cognitive, and social systems. It is grounded in the second law of thermodynamics and non-equilibrium steady-state thermodynamics, not analogy.


2. The Persistence Cost Functional

Let XX be a state space, ϕt(x)ϕt​(x) the flow of a dynamical system, and AXA⊆X an attractor set. Let δ(x)=d(x,A)δ(x)=d(x,A) be the distance from xx to the attractor. For a treatment of state-space constraints in viability theory, see Aubin (1991).

Definition 1 (Persistence Cost Functional): A persistence cost functional C(x)C(x) is a scalar function on XX satisfying:

  1. C(x)0C(x)≥0 for all xx
  2. C(x)=0C(x)=0 if and only if xAx∈A
  3. C(ϕt(x))L1([0,))C(ϕt​(x))∈L1([0,∞)) for all xx in the basin

Definition 2 (Cumulative Persistence Cost): For a finite horizon T>0T>0:DT(x)=0TC(ϕt(x))dtDT​(x)=∫0TC(ϕt​(x))dt

For trajectories that converge to the attractor:D(x)=0C(ϕt(x))dtD∞​(x)=∫0∞​C(ϕt​(x))dt


3. Existence and Lyapunov Equivalence

Theorem 1 (Existence of the Persistence Functional): Assume C(x)0C(x)≥0, C=0C=0 only on AA, and C(ϕt(x))L1([0,))C(ϕt​(x))∈L1([0,∞)) for all xx in the basin. Assume ff is locally Lipschitz, the flow is continuously differentiable in the initial condition, and CC is continuous and locally bounded. Then:

  1. D(x)=0C(ϕt(x))dtD∞​(x)=∫0∞​C(ϕt​(x))dt exists and is finite.
  2. DD∞​ is continuous.
  3. DD∞​ satisfies the transport equation:

D(x)f(x)=C(x)D∞​(x)⋅f(x)=−C(x)

Proof: The integral exists and is finite by the L1L1 assumption. Continuity follows from the dominated convergence theorem under the stated regularity assumptions. To derive the transport equation, compute:D(ϕh(x))=hC(ϕt(x))dt=D(x)0hC(ϕt(x))dtD(ϕh​(x))=∫h∞​C(ϕt​(x))dt=D(x)−∫0hC(ϕt​(x))dt

Then:D(ϕh(x))D(x)h=1h0hC(ϕt(x))dtC(x)hD(ϕh​(x))−D(x)​=−h1​∫0hC(ϕt​(x))dt→−C(x)

as h0h→0. By the chain rule:D(x)f(x)=C(x)D(x)⋅f(x)=−C(x)

Corollary (Equivalence to Lyapunov Theory): Any Lyapunov function V(x)V(x) (with V0V≥0, V=0V=0 on the attractor, and V˙0V˙≤0) yields a persistence cost C(x)=V˙(x)C(x)=−V˙(x). Conversely, any persistence cost C(x)C(x) satisfying Df=CDf=−C defines a Lyapunov function D(x)D(x).

Proof: If VV is a Lyapunov function, then V˙=Vf0V˙=∇Vf≤0. Define C=V˙C=−V˙. Then C0C≥0, C=0C=0 on the attractor, and DT=C=V(x)V(ϕT(x))DT​=∫C=V(x)−V(ϕT​(x)). Conversely, if Df=CDf=−C, then D˙=C0D˙=−C≤0, so DD is a Lyapunov function.

Interpretation: The persistence cost framework is mathematically equivalent to classical Lyapunov stability theory. For the connection to contraction analysis, see Lohmiller & Slotine (1998). For control Lyapunov functions, see Freeman & Kokotovic (1996). Entropy production is one physically meaningful realization of the cost function CC. For a detailed treatment of Lipschitz continuity of DD∞​ under a Lipschitz-flow hypothesis, see Galida (2026a), Proposition 4.


4. Entropy Production as Persistence Cost

4.1 Entropy Balance

For an open system, the entropy balance equation is:dSsystemdt=σΦdtdSsystem​​=σ−Φ

where σ0σ≥0 is the entropy production rate (always non-negative by the second law) and ΦΦ is the entropy export rate to the environment. For foundational treatments of stochastic thermodynamics and entropy production, see Seifert (2012) and Sekimoto (2010).

For a system in a steady state:dSsystemdt=0    σ=ΦdtdSsystem​​=0⟹σ

4.2 Excess Entropy Production

Define the steady-state entropy production rate σssσss​ as the rate when the system is at its attractor.

Define the excess entropy production rate:σexcess(x)=σ(x)σss(x)σexcess​(x)=σ(x)−σss​(x)

Assumption (Excess Entropy Decay): For all trajectories in the basin, there exist constants C<C<∞ and μ>0μ>0 such that:σexcess(ϕt(x))Ceμtσexcess(x)σexcess​(ϕt​(x))≤Ceμtσexcess​(x)

for all t0t≥0. This ensures D(x)<D∞​(x)<∞ and is the standard hypothesis under which the persistence functional and its associated bounds are well-defined, consistent with Galida (2026a, 2026b). The decay rate μμ may be domain-specific and is empirically measurable.

Note on generalization: The exponential decay assumption is adopted here to ensure finiteness of DD∞​ and to maintain consistency with the prior papers in this series. Generalization to L1L1 integrable decays (e.g., algebraic) is a priority for future work.

4.3 The Entropy Persistence Functional

Definition 3 (Cumulative Excess Entropy Functional): For a finite horizon T>0T>0:DT(x)=0Tσexcess(ϕt(x))dtDT​(x)=∫0Tσexcess​(ϕt​(x))dt

For trajectories that converge to the attractor:D(x)=0σexcess(ϕt(x))dtD∞​(x)=∫0∞​σexcess​(ϕt​(x))dt

Interpretation: The persistence functional is the total excess entropy generated during reconfiguration and recovery.

4.4 Corrective Permeability

Definition 4 (Corrective Permeability):κ=infxBAδ(x)D(x)κ=x∈B∖Ainf​D∞​(x)δ(x)​

where δ(x)=d(x,A)δ(x)=d(x,A) is the distance to the attractor.

Interpretation: κκ is the minimum excess entropy cost per unit distance. It measures the efficiency of reconfiguration: a system that returns with minimal excess entropy generation has high κκ; a system that generates excess entropy has low κκ.


4.5 Basin Depth

Proposition 1 (Properties of Basin Depth): Define B=D(saddle)B=D∞​(saddle), where saddlesaddle is the lowest point on the basin boundary (the separatrix between attractors). For the connection to large-deviation theory and escape rates, see Freidlin & Wentzell (2012). Then:

  1. B0B≥0, with equality iff the basin has no barrier (i.e., the boundary coincides with the attractor).
  2. For gradient systems x˙=V(x)x˙=−∇V(x), B=V(saddle)V(A)B=V(saddle)−V(A) (the classical energy barrier).
  3. BB is invariant under smooth coordinate changes (coordinate invariance).
  4. BB depends on the chosen persistence cost functional CC; different costs yield different barriers.

Proof: (1) follows from non-negativity of DD∞​. (2) follows from the transport equation Df=CDf=−C and the identity f=Vf=−∇V. (3) follows from the invariance of the integral under diffeomorphisms. (4) is self-evident.


5. Domain-Specific Realizations

5.1 Physical Systems: Thermodynamic Excess Entropy

For a thermodynamic system, S(x)=kBlogΩ(x)S(x)=kB​logΩ(x), where Ω(x)Ω(x) is the number of microstates. For an isolated system, σss=0σss​=0 (equilibrium), so σexcess=σ=S˙σexcess​=σ=S˙.κ=infxδ(x)S(A)S(x)κ=xinf​S(A)−S(x)δ(x)​

Example: A gas returning to equilibrium after compression. The entropy generated is ΔS=nRlog(Vf/Vi)ΔS=nRlog(Vf​/Vi​).

5.2 Biological Systems: Metabolic Excess Entropy

For a biological system, S(x)S(x) is the metabolic entropy. The baseline σssσss​ is the resting metabolic rate (homeostasis). The excess is:σexcess=metabolic rateresting metabolic rateσexcess​=metabolic rate−resting metabolic rateκ=infxδ(x)0σexcess(ϕt(x))dtκ=xinf​∫0∞​σexcess​(ϕt​(x))dtδ(x)​

Example: A cell returning to homeostasis after a nutrient shock. The excess entropy generated is the metabolic cost of restoring homeostasis above baseline. For the dissipative-structures framework underlying biological self-organization, see Nicolis & Prigogine (1989).

5.3 Cognitive Systems: Free Energy Dissipation

For a cognitive system, variational free energy F=logp(yx)+DKL[q()p(x)]F=−logp(yx)+DKL​[q(⋅)∥p(⋅∣x)] is adopted here as one candidate persistence functional. We do not claim variational free energy is uniquely correct; it is adopted as the most developed existing candidate persistence functional for cognitive systems. Other candidates (Bayesian surprise, expected free energy, predictive information) are possible; this paper focuses on FF due to its established role in the free-energy principle (Friston, 2010). For the thermodynamics of information and its connection to free-energy minimization, see Parrondo, Horowitz & Sagawa (2015) and Sagawa & Ueda (2008).

The baseline σssσss​ is the baseline neural dissipation rate (resting brain activity). The excess is:σexcess=F˙F˙ssσexcess​=F˙−F˙ssκ=infxδ(x)0σexcess(ϕt(x))dtκ=xinf​∫0∞​σexcess​(ϕt​(x))dtδ(x)​

Example: A cognitive system updating its beliefs after a prediction error. The excess entropy generated is the free energy dissipated during belief updating above baseline.

5.4 Social Systems: Coordination Excess Entropy

For a social system, define the aggregate social entropy production rate as:σsocial(t)=i(S˙i(t)S˙irest)σsocial(t)=i∑​(S˙i​(t)−S˙irest​)

where S˙i(t)S˙i​(t) is the total entropy production rate of individual ii, and S˙irestS˙irest​ is the individual’s baseline entropy production rate in a resting, minimally socially constrained state. This is measured via physiological proxies such as basal metabolic rate, resting allostatic load, or cortisol baseline (McEwen, 1998; Sterling & Eyer, 1988).

Interpretation: σsocialσsocial measures the excess dissipation attributable to social constraints: the additional entropy generated by coordination, communication, conflict, norm enforcement, and institutional friction.

Non-Negativity: Unlike total entropy production S˙i0S˙i​≥0 (which follows from the second law), σisocialσisocial​ is not guaranteed to be non-negative. Division of labor, infrastructure, and specialization may reduce an individual’s metabolic burden relative to a solitary baseline. The hypothesis is that during recovery from social disruption, σisocial0σisocial​≥0; in steady-state, σisocial0σisocial​→0. This is an empirical claim, not a theorem.

The baseline σssσss​ is the steady-state social entropy production rate (well-coordinated society). The excess is:σexcess=σsocialσssσexcess​=σsocial−σssκ=infxδ(x)0σexcess(ϕt(x))dtκ=xinf​∫0∞​σexcess​(ϕt​(x))dtδ(x)​

Example: A society recovering from a shock (economic crisis, political upheaval). The excess entropy generated is the coordination cost of restructuring above baseline. A harmonious society has σexcess=0σexcess​=0; a turbulent society has σexcess>0σexcess​>0; a chronically turbulent society may have settled into a new attractor with a higher σssσss​. This illustrates the framework’s central distinction: the attractor is the state of minimum entropy generation for that class of system.


6. The Unified Framework

6.1 Summary Table

Domain Entropy Functional Baseline σssσss Excess σexcessσexcess​ Recovery Rate κκ
Physical Thermodynamic entropy 0 (equilibrium) S˙S˙ infδΔSinfΔSδ
Biological Metabolic entropy Resting metabolic rate Metabolic rate — resting infδσexcessdtinf∫σexcess​dtδ
Cognitive Free energy Baseline neural dissipation F˙F˙ssF˙−F˙ss infδσexcessdtinf∫σexcess​dtδ
Social Social entropy production Steady-state social dissipation σsocialσssσsocial−σss infδσexcessdtinf∫σexcess​dtδ

6.2 The Universal Structure

Every domain follows the same mathematical structure:

Component Expression
Excess entropy production σexcess(x)=σ(x)σssσexcess​(x)=σ(x)−σss
Cumulative cost D(x)=0σexcess(ϕt(x))dtD∞​(x)=∫0∞​σexcess​(ϕt​(x))dt
Recovery rate κ=infxδ(x)/D(x)κ=infxδ(x)/D∞​(x)
Basin depth B=D(saddle)B=D∞​(saddle)
Transport equation Df=σexcessDf=−σexcess​

6.3 The Low-Energy Attractor Benchmark (Proposed Hypothesis)

We propose the following benchmark as an additional hypothesis: the attractor is the state of minimum entropy generation for that class of system.

Domain Attractor Entropy Generation at Attractor
Physical Equilibrium σ=0σ=0
Biological Homeostasis σ=σss>0σ=σss​>0 (resting metabolism)
Cognitive Settled Belief σ=σss>0σ=σss​>0 (baseline neural dissipation)
Social Coordinated Order σ=σss>0σ=σss​>0 (baseline institutional friction)

Interpretation:

  1. For equilibrium systems (gases, isolated systems), the attractor is the state where entropy generation reaches zero — the system has nowhere lower to go.
  2. For dissipative systems (cells, brains, societies), the attractor is the state where entropy generation reaches its lowest non-zero steady-state value — the minimum entropy generation the system can sustain while maintaining its functional organization.

Important caveats:

  • This is a proposed benchmark, not a derived theorem.
  • For cognitive systems in particular, minimizing entropy production rate (a thermodynamic quantity) and minimizing free energy/surprise (the actual claim in the free-energy principle) are distinct minimization principles. The framework does not establish a bridge between them; this is an open question.
  • The benchmark is an empirical hypothesis that requires domain-specific validation.

In all cases, the attractor is the lowest entropy-generating state that system can have while remaining itself.


7. Testable Predictions

7.1 Core Prediction

Prediction: The recovery rate κκ is inversely proportional to the excess entropy generated during reconfiguration:κ1DκD∞​1​

Falsification: If a system returns to its attractor with high excess entropy generation but high recovery rate, the prediction is falsified.

7.2 Secondary Prediction

Prediction: Systems that maintain their attractor with minimal excess entropy generation are more “efficient.” Systems that generate excess entropy are “inefficient” or “stressed.”

Falsification: If an inefficient system has lower excess entropy generation than an efficient system, the prediction is falsified.

7.3 Domain-Specific Predictions

Domain Prediction Falsification
Physical κκ correlates with thermal efficiency κκ high but efficiency low
Biological κκ correlates with metabolic efficiency κκ high but metabolic cost high
Cognitive κκ correlates with learning efficiency κκ high but learning cost high
Social κκ correlates with institutional efficiency κκ high but coordination cost high

8. Experimental Design

8.1 Physical Systems

  • System: Gas in a piston
  • Perturbation: Compression
  • Measurement: Excess entropy generation (heat measurement) and recovery time
  • Test: Correlation between κκ and 1/D1/D∞​

8.2 Biological Systems

  • System: Cell culture
  • Perturbation: Nutrient shock
  • Measurement: Metabolic rate above resting (oxygen consumption) and recovery time
  • Test: Correlation between κκ and metabolic cost

8.3 Cognitive Systems

  • System: Human participants in a learning task
  • Perturbation: Prediction error
  • Measurement: Free energy dissipation above baseline (EEG complexity, pupil dilation) and belief updating rate
  • Test: Correlation between κκ and free energy dissipation

8.4 Social Systems

  • System: Institutional response to shocks
  • Perturbation: Economic or political crisis
  • Measurement: Social entropy production above baseline (allostatic load, cortisol, institutional friction) and recovery time
  • Test: Correlation between κκ and social entropy production

9. Open Questions

Question Status Difficulty
Q1: Uniqueness of S(x)S(x) Are there multiple valid entropy functionals for a given domain? Hard
Q2: Variational principle Is there a universal variational principle that yields S(x)S(x)? Hard
Q3: Social second law Does σsocial0σsocial≥0 always hold during recovery? Very Hard
Q4: Cross-level entropy How does entropy generation at one level relate to entropy generation at another? Hard
Q5: Measurement Can we measure excess entropy generation in cognitive and social systems directly? Moderate
Q6: Unification Can all domain-specific entropy functionals be derived from a single universal functional? Very Hard

10. Conclusion

Every dissipative system maintains its attractor through continuous reconfiguration. Reconfiguration requires work; work generates excess entropy. The recovery rate κκ — corrective permeability — is the rate at which a system reconfigures to return to its attractor after perturbation. We have proposed that κκ is a measure of excess entropy generation rate.

We developed an abstract persistence cost framework and proved its equivalence to Lyapunov theory. We then identified entropy production as a physical realization of this cost, deriving:κ=infxδ(x)0σexcess(ϕt(x))dtκ=xinf​∫0∞​σexcess​(ϕt​(x))dtδ(x)​

where σexcess=σσssσexcess​=σσss​ is the excess entropy production rate above the system’s steady-state baseline — thermodynamic entropy for physical systems, metabolic entropy for biological systems, free energy dissipation for cognitive systems, and social entropy production for social systems.

We proposed a unified benchmark: the attractor is the state of minimum entropy generation for that class of system — zero for equilibrium systems, non-zero steady-state for dissipative systems. This provides a unified criterion for identifying attractors across domains: an attractor is a state from which the system cannot reduce its entropy generation further without losing its defining structure or function.

This unifies physical, biological, cognitive, and social systems. In each domain, persistence requires reconfiguration; reconfiguration generates excess entropy; κκ measures the entropy cost of that reconfiguration. The framework is grounded in the second law of thermodynamics and non-equilibrium steady-state thermodynamics, not analogy.

Social Application: The framework provides a thermodynamic interpretation of social dynamics: harmony is a low-entropy attractor state; turbulence is a high-entropy state generated by excess dissipation during reconfiguration. The recovery rate κκ measures how efficiently a society transitions from turbulence back to harmony — that is, how quickly it reduces its excess entropy production to zero.


11. Limitations

This paper establishes an abstract persistence cost framework with a proposed thermodynamic realization. Several limitations should be explicitly acknowledged:

  1. Uniqueness. Entropy production is not proved to be the unique persistence cost. Many positive functionals C(x)C(x) satisfy Df=CDf=−C. The identification of entropy production as the canonical cost is a physically motivated hypothesis, not a mathematical theorem.
  2. Scope. The framework does not imply that all domains obey thermodynamics literally. The cognitive and social realizations are proposed hypotheses requiring empirical validation.
  3. Decay assumption. Exponential decay of σexcessσexcess​ is a sufficient assumption to ensure finiteness of DD∞​, not a necessary one. Generalization to L1L1 integrable decays (e.g., algebraic) is a priority for future work.
  4. Basin depth. Basin depth B=D(saddle)B=D∞​(saddle) is defined in terms of the persistence cost functional. Its relationship to classical energy barriers is established only for gradient systems.
  5. Empirical validation. The predictions of the framework — particularly the inverse relationship between κκ and DD∞​ — remain to be tested empirically across domains.
  6. Low-energy attractor benchmark. The benchmark proposed in §6.3 is a hypothesis, not a derived theorem. For cognitive systems, it risks conflating thermodynamic entropy production with free-energy minimization — distinct principles whose relationship remains open.

References

Aubin, J. P. (1991). Viability Theory. Birkhäuser.

Boltzmann, L. (1877). “Über die Beziehung zwischen dem zweiten Hauptsatz der mechanischen Wärmetheorie und der Wahrscheinlichkeitsrechnung.” Wiener Berichte, 76, 373-435.

Clausius, R. (1865). “Über verschiedene für die Anwendung bequeme Formen der Hauptgleichungen der mechanischen Wärmetheorie.” Annalen der Physik, 125(7), 353-400.

Freeman, R. A., & Kokotovic, P. V. (1996). Robust Nonlinear Control Design: State-Space and Lyapunov Techniques. Birkhäuser.

Freidlin, M. I., & Wentzell, A. D. (2012). Random Perturbations of Dynamical Systems (3rd ed.). Springer.

Friston, K. (2010). “The free-energy principle: a unified brain theory?” Nature Reviews Neuroscience, 11(2), 127-138.

Galida, R. (2026a). “The Persistence Functional: A Candidate Formal Foundation for the Attractor Framework.” Fantasy Attractor.

Galida, R. (2026b). “Deriving Corrective Permeability from the Cumulative Deviation Functional.” Fantasy Attractor.

Jaynes, E. T. (1957). “Information Theory and Statistical Mechanics.” Physical Review, 106(4), 620-630.

Khalil, H. K. (2002). Nonlinear Systems (3rd ed.). Prentice Hall.

Kondepudi, D., & Prigogine, I. (1998). Modern Thermodynamics: From Heat Engines to Dissipative Structures. Wiley.

Lohmiller, W., & Slotine, J. J. E. (1998). “On contraction analysis for non-linear systems.” Automatica, 34(6), 683-696.

Lyapunov, A. M. (1892). The General Problem of the Stability of Motion.

McEwen, B. S. (1998). “Stress, Adaptation, and Disease: Allostasis and Allostatic Load.” Annals of the New York Academy of Sciences, 840(1), 33-44.

Nicolis, G., & Prigogine, I. (1989). Exploring Complexity: An Introduction. W. H. Freeman.

Parrondo, J. M. R., Horowitz, J. M., & Sagawa, T. (2015). “Thermodynamics of information.” Nature Physics, 11(2), 131-139.

Prigogine, I. (1947). Étude Thermodynamique des Phénomènes Irréversibles. Dunod.

Prigogine, I., & Nicolis, G. (1977). Self-Organization in Non-Equilibrium Systems. Wiley.

Sagawa, T., & Ueda, M. (2008). “Second law of thermodynamics with discrete quantum feedback control.” Physical Review Letters, 100(8), 080403.

Seifert, U. (2012). “Stochastic thermodynamics, fluctuation theorems and molecular machines.” Reports on Progress in Physics, 75(12), 126001.

Sekimoto, K. (2010). Stochastic Energetics. Springer.

Shannon, C. E. (1948). “A Mathematical Theory of Communication.” Bell System Technical Journal, 27(3), 379-423.

Sterling, P., & Eyer, J. (1988). “Allostasis: A New Paradigm to Explain Arousal Pathology.” In Handbook of Life Stress, Cognition and Health, 629-649.


Suggested citation: Galida, R. S. (2026). Excess Entropy Production as a Candidate Universal Cost of Persistence: A Thermodynamic Foundation for the Attractor Framework. Fantasy Attractor.




Deriving Corrective Permeability from the Cumulative Deviation Functional; Robert Galida (June 2026) [F]

Abstract

The attractor framework defines κκ (corrective permeability) as the rate at which a system returns to its attractor after perturbation. Historically, κκ has been treated as an empirical parameter — fitted to data rather than derived from first principles. This paper derives κκ from the framework’s foundational object: the cumulative deviation functional DT(x)=0Tδ(ϕt(x))dtDT​(x)=∫0Tδ(ϕt​(x))dt, where δ(x)=d(x,A)δ(x)=d(x,A).

We define:κ=infxBAδ(x)D(x)κ=x∈B∖Ainf​D∞​(x)δ(x)​

We prove that for linear systems x˙=Axx˙=−Ax with AA symmetric positive definite, this definition recovers the slowest eigenvalue λmin(A)λmin​(A) — the conventional notion of corrective permeability. We establish a sharp universal persistence bound D(x)δ(x)/κD∞​(x)≤δ(x)/κ, show homogeneity and scale invariance of the variational ratio, and demonstrate consistency with Koopman spectral theory and resolvent poles for finite-dimensional linear systems. A comparison theorem links κκ to classical exponential stability constants. A Hamilton-Jacobi-type transport equation for DD∞​ is derived. A finite-horizon estimator κT=infxδ(x)DT(x)κT​=infxDT​(x)δ(x)​ is provided with exponential convergence under explicit assumptions.

The derivation is rigorous for linear systems and testable. Open questions for nonlinear, multiscale, and stochastic systems are identified.

Keywords: corrective permeability, cumulative deviation functional, attractor framework, Koopman operator, trajectory functional


1. Introduction

The attractor framework has been applied across physics, biology, cognition, and social systems. Its central variable — corrective permeability κκ — measures the rate at which a system returns to its attractor after perturbation. Historically, κκ has been defined empirically as κ=1/τκ=1/τ, where ττ is a measured recovery time constant.

This paper derives κκ from a single foundational object: the cumulative deviation functional DT(x)DT​(x). Within the present framework, κκ is defined variationally rather than introduced as an empirical fitting parameter. We show that κκ is a consequence of the trajectory geometry — specifically, the ratio of initial distance to total cumulative deviation.

The derivation is rigorous for linear systems, connects to established theory (Koopman operators, resolvent poles), and provides a finite-horizon estimator for empirical use. Open questions for nonlinear and stochastic systems are identified.


2. The Cumulative Deviation Functional

Let XX be a metric space with distance function ∥⋅∥. Let ϕt(x)ϕt​(x) be the flow of a dynamical system starting from state xXx∈X at time t=0t=0. Let AXA⊆X be an attractor set (a compact, invariant set to which trajectories converge). Let BB be the basin of attraction of AA.

Define the distance from a point to the attractor:δ(x)=d(x,A)=infaAxaδ(x)=d(x,A)=a∈Ainf​∥xa

Definition 1 (Cumulative Deviation Functional): For a finite horizon T>0T>0, define:DT(x)=0Tδ(ϕt(x))dtDT​(x)=∫0Tδ(ϕt​(x))dt

For TT→∞, define:D(x)=0δ(ϕt(x))dtD∞​(x)=∫0∞​δ(ϕt​(x))dt

Proposition 1 (Finiteness of D∞D∞​): Assume there exist constants C<C<∞ and μ>0μ>0 such that:δ(ϕt(x))Ceμtδ(x)δ(ϕt​(x))≤Ceμtδ(x)

for all xBx∈B. Then D(x)<D∞​(x)<∞ for every xBx∈B.

Proof:D(x)=0δ(ϕt(x))dt0Ceμtδ(x)dt=Cμδ(x)<D∞​(x)=∫0∞​δ(ϕt​(x))dt≤∫0∞​Ceμtδ(x)dt=μCδ(x)<∞

Properties (from Galida, 2026a):

Property Statement
Non-negativity DT(x)0DT​(x)≥0
Monotonicity DT2(x)DT1(x)DT2​​(x)≥DT1​​(x) for T2T1T2​≥T1​
Additivity DT+S(x)=DT(x)+DS(ϕT(x))DT+S​(x)=DT​(x)+DS​(ϕT​(x))
Instantaneous growth ddTDT(x)=δ(ϕT(x))dTdDT​(x)=δ(ϕT​(x))
Occupation measure DT(x)=δ(y)dμT(y)DT​(x)=∫δ(y)dμT​(y), where μTμT​ is the occupation measure

3. Derivation of Corrective Permeability (κκ)

3.1 Variational Definition

Definition 2 (Corrective Permeability):κ=infxBAδ(x)D(x)κ=x∈B∖Ainf​D∞​(x)δ(x)​

Interpretation: κκ is the effective recovery rate — the smallest ratio of initial distance to total cumulative deviation. It serves as a global measure of the slowest recovery mode in the basin.

Remark on κκ: The definition allows κ=0κ=0 if D(x)D∞​(x) diverges or if the ratio δ(x)/D(x)δ(x)/D∞​(x) can be made arbitrarily small. Throughout the remainder of this paper, we assume hypotheses (such as the exponential stability in Proposition 1) that guarantee κ>0κ>0.

Remark on attainment: The infimum in the definition of κκ need not be attained; minimizing sequences may exist without a minimizing state. For linear systems, the infimum is attained on the slow eigenspace.


3.2 Homogeneity and Scale Invariance

Theorem 1 (Homogeneity and Scale Invariance): Suppose the flow satisfies ϕt(αx)=αϕt(x)ϕt​(αx)=αϕt​(x) for all tt and all α>0α>0, and the distance function satisfies δ(αx)=αδ(x)δ(αx)=αδ(x). Then:δ(αx)D(αx)=δ(x)D(x)D∞​(αx)δ(αx)​=D∞​(x)δ(x)​

Proof:D(αx)=0δ(ϕt(αx))dt=0δ(αϕt(x))dt=α0δ(ϕt(x))dt=αD(x)D∞​(αx)=∫0∞​δ(ϕt​(αx))dt=∫0∞​δ(αϕt​(x))dt=α∫0∞​δ(ϕt​(x))dt=αD∞​(x)

Corollary: For linear systems, the infimum over all x0x=0 reduces to an infimum over the unit sphere:κ=infx=1δ(x)D(x)κ=∥x∥=1inf​D∞​(x)δ(x)​


3.3 Sharp Universal Persistence Bound

Theorem 2 (Sharp Universal Persistence Bound): For any xBAx∈B∖A:D(x)δ(x)κD∞​(x)≤κδ(x)​

Moreover, the constant 1/κ1/κ is optimal: it is the smallest constant such that this inequality holds for all xx in the basin.

Proof: By definition of κκ as the infimum of δ(x)/D(x)δ(x)/D∞​(x), we have δ(x)/D(x)κδ(x)/D∞​(x)≥κ for all xx. Rearranging gives:D(x)δ(x)κD∞​(x)≤κδ(x)​

Optimality follows from Theorem 3: for the slow eigenvector v1v1​, D(v1)=δ(v1)/κD∞​(v1​)=δ(v1​)/κ, so no smaller constant can work.


3.4 Consistency with Linear Systems

Consider a linear system x˙=Axx˙=−Ax, with AA symmetric positive definite. Let its eigenvalues be 0<λ1λ2λn0<λ1​≤λ2​≤⋯≤λn​, with corresponding orthonormal eigenvectors v1,v2,,vnv1​,v2​,…,vn​.

The flow is ϕt(x)=eAtxϕt​(x)=eAtx. The attractor is A={0}A={0}, and the distance to the attractor is δ(x)=xδ(x)=∥x∥.

Theorem 3 (Linear Consistency): For x˙=Axx˙=−Ax with AA symmetric positive definite,infx0xD(x)=λmin(A)x=0inf​D∞​(x)∥x∥​=λmin​(A)

Proof:

Since AA is symmetric positive definite, eAteAt is symmetric positive definite with eigenvalues eλiteλit. Hence its operator norm is eAt=eλ1teAt∥=eλ1​t. For any x0x=0:D(x)=0eAtxdt0xeλ1tdt=xλ1D∞​(x)=∫0∞​∥eAtxdt≤∫0∞​∥xeλ1​tdt=λ1​∥x∥​

Therefore:xD(x)λ1D∞​(x)∥x∥​≥λ1​

To show equality is achieved, take x=v1x=v1​ (the eigenvector corresponding to λ1λ1​). Then:eAtv1=v1eλ1teAtv1​∥=∥v1​∥eλ1​t

and:D(v1)=0v1eλ1tdt=v1λ1D∞​(v1​)=∫0∞​∥v1​∥eλ1​tdt=λ1​∥v1​∥​

Thus:v1D(v1)=λ1D∞​(v1​)∥v1​∥​=λ1​

Hence:infx0xD(x)=λ1x=0inf​D∞​(x)∥x∥​=λ1​

Corollary: For linear systems, the variational definition of κκ recovers the slowest eigenvalue — the conventional notion of corrective permeability.


3.5 Transport Equation

Theorem 4 (Transport Equation): Assume the vector field ff is C1C1, the flow ϕtϕt​ is C1C1, and DD∞​ is continuously differentiable on BAB∖A. Then:D(x)f(x)=δ(x)D∞​(x)⋅f(x)=−δ(x)

Proof: From the definition:D(ϕs(x))=D(x)Ds(x)D∞​(ϕs​(x))=D∞​(x)−Ds​(x)

Differentiating with respect to ss at s=0s=0:ddsD(ϕs(x))s=0=δ(x)dsdD∞​(ϕs​(x))​s=0​=−δ(x)

By the chain rule:D(x)f(x)=δ(x)D∞​(x)⋅f(x)=−δ(x)

Interpretation: This is a first-order transport equation, fD=δf⋅∇D=−δ, which belongs to the broader Hamilton-Jacobi family but lacks a Hamiltonian in the usual sense. It may serve as a foundation for numerical computation and further theoretical development.


3.6 Local vs. Global Interpretation

The variational definition κ=infxδ(x)D(x)κ=infxD∞​(x)δ(x)​ is global — it is the slowest recovery rate over the entire basin. This is not necessarily the same as the local recovery rate near the attractor (the slowest eigenvalue of the linearization). For linear systems, they coincide. For nonlinear systems, they may differ if transient excursions produce slower effective recovery than the local linearization predicts.

This distinction is important: κκ is a global invariant of the basin, not merely a local property of the attractor. The relationship between the global κκ and the local Lyapunov exponent is an open question (see §6).


3.7 Non-Symmetric Linear Systems

For a general linear system x˙=Axx˙=Ax (where AA is stable, i.e., all eigenvalues have negative real parts), the same principle holds in the diagonalizable case. The slowest mode corresponds to the eigenvalue with the largest real part (closest to zero).

Conjecture: An analogous result holds for non-normal linear systems under additional assumptions on the semigroup, such as a uniformly exponentially stable semigroup satisfying suitable norm bounds. This remains an open question.


3.8 Comparison with Exponential Stability

Theorem 5 (Comparison with Exponential Stability): Suppose the system satisfies the exponential stability bound:δ(ϕt(x))Ceμtδ(x)δ(ϕt​(x))≤Ceμtδ(x)

for all xBx∈B, with constants C<C<∞ and μ>0μ>0. Then:κμCκCμ

Proof: From the stability bound:D(x)=0δ(ϕt(x))dt0Ceμtδ(x)dt=Cμδ(x)D∞​(x)=∫0∞​δ(ϕt​(x))dt≤∫0∞​Ceμtδ(x)dt=μCδ(x)

Therefore:δ(x)D(x)μCD∞​(x)δ(x)​≥Cμ

Taking the infimum over xx:κ=infxδ(x)D(x)μCκ=xinf​D∞​(x)δ(x)​≥Cμ

Interpretation: The variational constant κκ is bounded below by the exponential stability constant μ/Cμ/C.


4. Connections to Existing Theory

4.1 Koopman Operator

The Koopman operator KtKt acts on observables as:(Ktf)(x)=f(ϕt(x))(Ktf)(x)=f(ϕt​(x))

For linear systems x˙=Axx˙=−Ax, the Koopman eigenvalues are eλiteλit. The dominant nontrivial eigenvalue (largest less than 1) is eλ1teλ1​t, corresponding to the slowest decay rate.

For finite-dimensional linear systems, ρ=eλmintρ=eλmin​t, and therefore:1tlogρ=λmin=κt1​logρ=λmin​=κ

Thus, under the hypotheses of Theorem 3, the variational constant equals the exponential decay rate associated with the dominant Koopman eigenvalue.


4.2 Resolvent Poles

For finite-dimensional stable linear systems, the resolvent (sI+A)1(sI+A)−1 has poles at s=λis=−λi​. The pole closest to the imaginary axis is s=λ1s=−λ1​.

Since Theorem 3 identifies κ=λminκ=λmin​, and the resolvent poles are si=λisi​=−λi​, we obtain:κ=mini(si)κ=imin​∣ℜ(si​)∣

for finite-dimensional linear systems.


5. Finite-Horizon Estimation

In practice, we can only measure finite trajectories. Define the finite-horizon estimator:κT=infxKδ(x)DT(x)κT​=x∈Kinf​DT​(x)δ(x)​

where KBK⊂B is compact and KA=K∩A=∅.

Proposition 2 (Finite-Horizon Estimation): Assume:

  1. The flow ϕt(x)ϕt​(x) is jointly continuous in (t,x)(t,x).
  2. δ(x)δ(x) is continuous.
  3. The exponential stability bound δ(ϕt(x))Ceμtδ(x)δ(ϕt​(x))≤Ceμtδ(x) holds uniformly for all xKx∈K, with μ>0μ>0.

Then the variational constant κκ (from Definition 2) satisfies κμ/Cκμ/C by Theorem 5, and:κTκas TκT​→κas T→∞

with error:κTκ=O(eμT)κT​−κ∣=O(eμT)

Proof: For any xKx∈K, the tail bound gives:D(x)DT(x)=Tδ(ϕt(x))dtCeμTδ(x)μD∞​(x)−DT​(x)∣=∫T∞​δ(ϕt​(x))dtμCeμTδ(x)​

Since δ(x)δ(x) is bounded on the compact set KK, let M=supxKδ(x)<M=supx∈K​δ(x)<∞. Then:D(x)DT(x)CMeμTμD∞​(x)−DT​(x)∣≤μCMeμT

The right-hand side is independent of xx and tends to zero as TT→∞. Hence DTDDT​→D∞​ uniformly on KK.

Moreover, since KK is compact and KA=K∩A=∅, continuity of δδ gives infxKδ(x)>0infx∈K​δ(x)>0. Since DT(x)DT​(x) is continuous (by assumptions 1–2) and monotonically non-decreasing in TT (from §2), for any fixed finite T0>0T0​>0, D(x)DT0(x)D∞​(x)≥DT0​​(x), and DT0DT0​​ is continuous and strictly positive on KK. A continuous, strictly positive function on a compact set has a positive infimum:m=infxKDT0(x)>0m=x∈Kinf​DT0​​(x)>0

Thus:infxKD(x)m>0x∈Kinf​D∞​(x)≥m>0

Uniform convergence of DTDT​ to DD∞​ on KK therefore implies uniform convergence of δ(x)/DT(x)δ(x)/DT​(x) to δ(x)/D(x)δ(x)/D∞​(x). Consequently, the infima converge.


6. Open Questions

Question Status Difficulty
Q1: Nonlinear systems Does infδDinfD∞​δ​ equal the local Lyapunov exponent? Hard
Q2: Local vs. global consistency Does limxAδ(x)D(x)=κlimx→A​D∞​(x)δ(x)​=κ hold for general nonlinear systems? Hard
Q3: Non-normal systems Does the infimum equal the slowest eigenvalue for non-normal AA? Moderate
Q4: Multiple timescales Does the infimum isolate the slowest timescale? Hard
Q5: Stochastic systems How does noise affect the finite-horizon estimator? Hard
Q6: Multiple attractors How does κκ behave in basins with multiple attractors? Moderate

7. Conclusion

This paper derives corrective permeability κκ from the cumulative deviation functional DT(x)DT​(x). The variational definition:κ=infxδ(x)D(x)κ=xinf​D∞​(x)δ(x)​

is shown to recover the slowest eigenvalue for linear systems, consistent with the conventional empirical definition κ=1/τκ=1/τ. A sharp universal persistence bound D(x)δ(x)/κD∞​(x)≤δ(x)/κ is established. A comparison theorem links κκ to classical exponential stability constants. A Hamilton-Jacobi-type transport equation for DD∞​ is derived. Connections to Koopman theory and resolvent theory are established for finite-dimensional linear systems. A finite-horizon estimator κTκT​ is provided with exponential convergence under explicit assumptions.

Key contribution: Within the present framework, κκ is defined variationally rather than introduced as an empirical fitting parameter — at least for the class of systems analyzed here.

Next steps: Extend the derivation to nonlinear systems (Q1–Q2), non-normal systems (Q3), multiple timescales (Q4), and stochastic dynamics (Q5).


References

Crandall, M. G., Ishii, H., & Lions, P. L. (1992). “User’s Guide to Viscosity Solutions of Second Order Partial Differential Equations.” Bulletin of the American Mathematical Society, 27(1), 1-67.

Evans, L. C. (2010). Partial Differential Equations. American Mathematical Society.

Galida, R. (2026a). “The Persistence Functional: A Candidate Formal Foundation for the Attractor Framework.” Fantasy Attractor.

Hale, J. K. (1988). Asymptotic Behavior of Dissipative Systems. American Mathematical Society.

Hirsch, M. W., Smale, S., & Devaney, R. L. (2004). Differential Equations, Dynamical Systems, and an Introduction to Chaos (2nd ed.). Elsevier Academic Press.

Khalil, H. K. (2002). Nonlinear Systems (3rd ed.). Prentice Hall.

Koopman, B. O. (1931). “Hamiltonian Systems and Transformations in Hilbert Space.” Proceedings of the National Academy of Sciences, 17(5), 315-318.

Lyapunov, A. M. (1892). The General Problem of the Stability of Motion. (English translation: 1992, Taylor & Francis).

Mezić, I. (2005). “Spectral Properties of Dynamical Systems, Model Reduction and Decompositions.” Nonlinear Dynamics, 41(1-3), 309-325.

Pazy, A. (1983). Semigroups of Linear Operators and Applications to Partial Differential Equations. Springer.

Vidyasagar, M. (1993). Nonlinear Systems Analysis (2nd ed.). Prentice Hall.


Suggested citation: Galida, R. S. (2026). Deriving Corrective Permeability from the Cumulative Deviation Functional. Fantasy Attractor.