Home » Posts tagged 'Lyapunov Theory'
Tag Archives: Lyapunov Theory
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:κ=xinf∫0∞σexcess(ϕt(x))dtδ(x)
where σ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 X be a state space, ϕt(x) the flow of a dynamical system, and A⊆X an attractor set. Let δ(x)=d(x,A) be the distance from x 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) is a scalar function on X satisfying:
- C(x)≥0 for all x
- C(x)=0 if and only if x∈A
- C(ϕt(x))∈L1([0,∞)) for all x in the basin
Definition 2 (Cumulative Persistence Cost): For a finite horizon T>0:DT(x)=∫0TC(ϕt(x))dt
For trajectories that converge to the attractor:D∞(x)=∫0∞C(ϕt(x))dt
3. Existence and Lyapunov Equivalence
Theorem 1 (Existence of the Persistence Functional): Assume C(x)≥0, C=0 only on A, and C(ϕt(x))∈L1([0,∞)) for all x in the basin. Assume f is locally Lipschitz, the flow is continuously differentiable in the initial condition, and C is continuous and locally bounded. Then:
- D∞(x)=∫0∞C(ϕt(x))dt exists and is finite.
- D∞ is continuous.
- D∞ satisfies the transport equation:
∇D∞(x)⋅f(x)=−C(x)
Proof: The integral exists and is finite by the L1 assumption. Continuity follows from the dominated convergence theorem under the stated regularity assumptions. To derive the transport equation, compute:D(ϕh(x))=∫h∞C(ϕt(x))dt=D(x)−∫0hC(ϕt(x))dt
Then:hD(ϕh(x))−D(x)=−h1∫0hC(ϕt(x))dt→−C(x)
as h→0. By the chain rule:∇D(x)⋅f(x)=−C(x)□
Corollary (Equivalence to Lyapunov Theory): Any Lyapunov function V(x) (with V≥0, V=0 on the attractor, and V˙≤0) yields a persistence cost C(x)=−V˙(x). Conversely, any persistence cost C(x) satisfying ∇D⋅f=−C defines a Lyapunov function D(x).
Proof: If V is a Lyapunov function, then V˙=∇V⋅f≤0. Define C=−V˙. Then C≥0, C=0 on the attractor, and DT=∫C=V(x)−V(ϕT(x)). Conversely, if ∇D⋅f=−C, then D˙=−C≤0, so D 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 C. For a detailed treatment of Lipschitz continuity of D∞ 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:dtdSsystem=σ−Φ
where σ≥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:dtdSsystem=0⟹σ=Φ
4.2 Excess Entropy Production
Define the steady-state entropy production rate σss as the rate when the system is at its attractor.
Define the excess entropy production rate:σexcess(x)=σ(x)−σss(x)
Assumption (Excess Entropy Decay): For all trajectories in the basin, there exist constants C<∞ and μ>0 such that:σexcess(ϕt(x))≤Ce−μtσexcess(x)
for all t≥0. This ensures 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 D∞ and to maintain consistency with the prior papers in this series. Generalization to L1 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>0:DT(x)=∫0Tσexcess(ϕt(x))dt
For trajectories that converge to the attractor:D∞(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):κ=x∈B∖AinfD∞(x)δ(x)
where δ(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), where saddle 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:
- B≥0, with equality iff the basin has no barrier (i.e., the boundary coincides with the attractor).
- For gradient systems x˙=−∇V(x), B=V(saddle)−V(A) (the classical energy barrier).
- B is invariant under smooth coordinate changes (coordinate invariance).
- B depends on the chosen persistence cost functional C; different costs yield different barriers.
Proof: (1) follows from non-negativity of D∞. (2) follows from the transport equation ∇D⋅f=−C and the identity f=−∇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), where Ω(x) is the number of microstates. For an isolated system, σss=0 (equilibrium), so σexcess=σ=S˙.κ=xinfS(A)−S(x)δ(x)
Example: A gas returning to equilibrium after compression. The entropy generated is ΔS=nRlog(Vf/Vi).
5.2 Biological Systems: Metabolic Excess Entropy
For a biological system, S(x) is the metabolic entropy. The baseline σss is the resting metabolic rate (homeostasis). The excess is:σexcess=metabolic rate−resting metabolic rateκ=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(y∣x)+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 F 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 is the baseline neural dissipation rate (resting brain activity). The excess is:σexcess=F˙−F˙ssκ=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)
where S˙i(t) is the total entropy production rate of individual i, and S˙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 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˙i≥0 (which follows from the second law), σ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, σisocial≥0; in steady-state, σisocial→0. This is an empirical claim, not a theorem.
The baseline σss is the steady-state social entropy production rate (well-coordinated society). The excess is:σexcess=σsocial−σssκ=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; a turbulent society has σexcess>0; a chronically turbulent society may have settled into a new attractor with a higher σ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 | Excess σexcess | Recovery Rate κ |
|---|---|---|---|---|
| Physical | Thermodynamic entropy | 0 (equilibrium) | S˙ | infΔSδ |
| Biological | Metabolic entropy | Resting metabolic rate | Metabolic rate — resting | inf∫σexcessdtδ |
| Cognitive | Free energy | Baseline neural dissipation | F˙−F˙ss | inf∫σexcessdtδ |
| Social | Social entropy production | Steady-state social dissipation | σsocial−σss | inf∫σexcessdtδ |
6.2 The Universal Structure
Every domain follows the same mathematical structure:
| Component | Expression |
|---|---|
| Excess entropy production | σexcess(x)=σ(x)−σss |
| Cumulative cost | D∞(x)=∫0∞σexcess(ϕt(x))dt |
| Recovery rate | κ=infxδ(x)/D∞(x) |
| Basin depth | B=D∞(saddle) |
| Transport equation | ∇D⋅f=−σ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 |
| Biological | Homeostasis | σ=σss>0 (resting metabolism) |
| Cognitive | Settled Belief | σ=σss>0 (baseline neural dissipation) |
| Social | Coordinated Order | σ=σss>0 (baseline institutional friction) |
Interpretation:
- For equilibrium systems (gases, isolated systems), the attractor is the state where entropy generation reaches zero — the system has nowhere lower to go.
- 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:κ∝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/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)? | Hard |
| Q3: Social second law | Does σ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:κ=xinf∫0∞σexcess(ϕt(x))dtδ(x)
where σ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:
- Uniqueness. Entropy production is not proved to be the unique persistence cost. Many positive functionals C(x) satisfy ∇D⋅f=−C. The identification of entropy production as the canonical cost is a physically motivated hypothesis, not a mathematical theorem.
- Scope. The framework does not imply that all domains obey thermodynamics literally. The cognitive and social realizations are proposed hypotheses requiring empirical validation.
- Decay assumption. Exponential decay of σexcess is a sufficient assumption to ensure finiteness of D∞, not a necessary one. Generalization to L1 integrable decays (e.g., algebraic) is a priority for future work.
- Basin depth. Basin depth B=D∞(saddle) is defined in terms of the persistence cost functional. Its relationship to classical energy barriers is established only for gradient systems.
- Empirical validation. The predictions of the framework — particularly the inverse relationship between κ and D∞ — remain to be tested empirically across domains.
- 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))dt, where δ(x)=d(x,A).
We define:κ=x∈B∖AinfD∞(x)δ(x)
We prove that for linear systems x˙=−Ax with A symmetric positive definite, this definition recovers the slowest eigenvalue λmin(A) — the conventional notion of corrective permeability. We establish a sharp universal persistence bound 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 D∞ is derived. A finite-horizon estimator κ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/τ, where τ is a measured recovery time constant.
This paper derives κ from a single foundational object: the cumulative deviation functional 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 X be a metric space with distance function ∥⋅∥. Let ϕt(x) be the flow of a dynamical system starting from state x∈X at time t=0. Let A⊆X be an attractor set (a compact, invariant set to which trajectories converge). Let B be the basin of attraction of A.
Define the distance from a point to the attractor:δ(x)=d(x,A)=a∈Ainf∥x−a∥
Definition 1 (Cumulative Deviation Functional): For a finite horizon T>0, define:DT(x)=∫0Tδ(ϕt(x))dt
For T→∞, define:D∞(x)=∫0∞δ(ϕt(x))dt
Proposition 1 (Finiteness of D∞D∞): Assume there exist constants C<∞ and μ>0 such that:δ(ϕt(x))≤Ce−μtδ(x)
for all x∈B. Then D∞(x)<∞ for every x∈B.
Proof:D∞(x)=∫0∞δ(ϕt(x))dt≤∫0∞Ce−μtδ(x)dt=μCδ(x)<∞□
Properties (from Galida, 2026a):
| Property | Statement |
|---|---|
| Non-negativity | DT(x)≥0 |
| Monotonicity | DT2(x)≥DT1(x) for T2≥T1 |
| Additivity | DT+S(x)=DT(x)+DS(ϕT(x)) |
| Instantaneous growth | dTdDT(x)=δ(ϕT(x)) |
| Occupation measure | DT(x)=∫δ(y)dμT(y), where μT is the occupation measure |
3. Derivation of Corrective Permeability (κ)
3.1 Variational Definition
Definition 2 (Corrective Permeability):κ=x∈B∖AinfD∞(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 if D∞(x) diverges or if the ratio δ(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.
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) for all t and all α>0, and the distance function satisfies δ(αx)=αδ(x). Then:D∞(αx)δ(αx)=D∞(x)δ(x)
Proof:D∞(αx)=∫0∞δ(ϕt(αx))dt=∫0∞δ(αϕt(x))dt=α∫0∞δ(ϕt(x))dt=αD∞(x)
Corollary: For linear systems, the infimum over all x=0 reduces to an infimum over the unit sphere:κ=∥x∥=1infD∞(x)δ(x)
3.3 Sharp Universal Persistence Bound
Theorem 2 (Sharp Universal Persistence Bound): For any x∈B∖A:D∞(x)≤κδ(x)
Moreover, the constant 1/κ is optimal: it is the smallest constant such that this inequality holds for all x in the basin.
Proof: By definition of κ as the infimum of δ(x)/D∞(x), we have δ(x)/D∞(x)≥κ for all x. Rearranging gives:D∞(x)≤κδ(x)
Optimality follows from Theorem 3: for the slow eigenvector v1, D∞(v1)=δ(v1)/κ, so no smaller constant can work.□
3.4 Consistency with Linear Systems
Consider a linear system x˙=−Ax, with A symmetric positive definite. Let its eigenvalues be 0<λ1≤λ2≤⋯≤λn, with corresponding orthonormal eigenvectors v1,v2,…,vn.
The flow is ϕt(x)=e−Atx. The attractor is A={0}, and the distance to the attractor is δ(x)=∥x∥.
Theorem 3 (Linear Consistency): For x˙=−Ax with A symmetric positive definite,x=0infD∞(x)∥x∥=λmin(A)
Proof:
Since A is symmetric positive definite, e−At is symmetric positive definite with eigenvalues e−λit. Hence its operator norm is ∥e−At∥=e−λ1t. For any x=0:D∞(x)=∫0∞∥e−Atx∥dt≤∫0∞∥x∥e−λ1tdt=λ1∥x∥
Therefore:D∞(x)∥x∥≥λ1
To show equality is achieved, take x=v1 (the eigenvector corresponding to λ1). Then:∥e−Atv1∥=∥v1∥e−λ1t
and:D∞(v1)=∫0∞∥v1∥e−λ1tdt=λ1∥v1∥
Thus:D∞(v1)∥v1∥=λ1
Hence:x=0infD∞(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 f is C1, the flow ϕt is C1, and D∞ is continuously differentiable on B∖A. Then:∇D∞(x)⋅f(x)=−δ(x)
Proof: From the definition:D∞(ϕs(x))=D∞(x)−Ds(x)
Differentiating with respect to s at s=0:dsdD∞(ϕs(x))s=0=−δ(x)
By the chain rule:∇D∞(x)⋅f(x)=−δ(x)□
Interpretation: This is a first-order transport equation, 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 κ=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˙=Ax (where A 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)
for all x∈B, with constants C<∞ and μ>0. Then:κ≥Cμ
Proof: From the stability bound:D∞(x)=∫0∞δ(ϕt(x))dt≤∫0∞Ce−μtδ(x)dt=μCδ(x)
Therefore:D∞(x)δ(x)≥Cμ
Taking the infimum over x:κ=xinfD∞(x)δ(x)≥Cμ□
Interpretation: The variational constant κ is bounded below by the exponential stability constant μ/C.
4. Connections to Existing Theory
4.1 Koopman Operator
The Koopman operator Kt acts on observables as:(Ktf)(x)=f(ϕt(x))
For linear systems x˙=−Ax, the Koopman eigenvalues are e−λit. The dominant nontrivial eigenvalue (largest less than 1) is e−λ1t, corresponding to the slowest decay rate.
For finite-dimensional linear systems, ρ=e−λmint, and therefore:−t1logρ=λ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 has poles at s=−λi. The pole closest to the imaginary axis is s=−λ1.
Since Theorem 3 identifies κ=λmin, and the resolvent poles are si=−λi, we obtain:κ=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=x∈KinfDT(x)δ(x)
where K⊂B is compact and K∩A=∅.
Proposition 2 (Finite-Horizon Estimation): Assume:
- The flow ϕt(x) is jointly continuous in (t,x).
- δ(x) is continuous.
- The exponential stability bound δ(ϕt(x))≤Ce−μtδ(x) holds uniformly for all x∈K, with μ>0.
Then the variational constant κ (from Definition 2) satisfies κ≥μ/C by Theorem 5, and:κT→κas T→∞
with error:∣κT−κ∣=O(e−μT)
Proof: For any x∈K, the tail bound gives:∣D∞(x)−DT(x)∣=∫T∞δ(ϕt(x))dt≤μCe−μTδ(x)
Since δ(x) is bounded on the compact set K, let M=supx∈Kδ(x)<∞. Then:∣D∞(x)−DT(x)∣≤μCMe−μT
The right-hand side is independent of x and tends to zero as T→∞. Hence DT→D∞ uniformly on K.
Moreover, since K is compact and K∩A=∅, continuity of δ gives infx∈Kδ(x)>0. Since DT(x) is continuous (by assumptions 1–2) and monotonically non-decreasing in T (from §2), for any fixed finite T0>0, D∞(x)≥DT0(x), and DT0 is continuous and strictly positive on K. A continuous, strictly positive function on a compact set has a positive infimum:m=x∈KinfDT0(x)>0
Thus:x∈KinfD∞(x)≥m>0
Uniform convergence of DT to D∞ on K therefore implies uniform convergence of δ(x)/DT(x) to δ(x)/D∞(x). Consequently, the infima converge.□
6. Open Questions
| Question | Status | Difficulty |
|---|---|---|
| Q1: Nonlinear systems | Does infD∞δ equal the local Lyapunov exponent? | Hard |
| Q2: Local vs. global consistency | Does limx→AD∞(x)δ(x)=κ hold for general nonlinear systems? | Hard |
| Q3: Non-normal systems | Does the infimum equal the slowest eigenvalue for non-normal A? | 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). The variational definition:κ=xinfD∞(x)δ(x)
is shown to recover the slowest eigenvalue for linear systems, consistent with the conventional empirical definition κ=1/τ. A sharp universal persistence bound D∞(x)≤δ(x)/κ is established. A comparison theorem links κ to classical exponential stability constants. A Hamilton-Jacobi-type transport equation for D∞ is derived. Connections to Koopman theory and resolvent theory are established for finite-dimensional linear systems. A finite-horizon estimator κ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.

