Home » Posts tagged 'recovery rate'

Tag Archives: recovery rate

The Attractor Framework in Astrophysics: Persistence, Entropy, and Gravitational Systems; Robert Galida (July 2026) [A]

Abstract

The attractor framework provides a domain-general vocabulary for describing persistence and change across physical, biological, cognitive, and social systems. This paper extends the framework to astrophysical dissipative systems. We distinguish between conservative gravitational dynamics — which define families of stable invariant solutions — and dissipative processes — which select and can stabilize particular configurations within those families.

The central thesis is:

Gravity defines the landscape. Dissipation selects the configuration.

We provide an operational definition of the excess entropy production functional σexcessσexcess​ for gravitational systems, grounding the persistence functional D=σexcessdtD∞​=∫σexcess​dt in physical dissipation rates above steady-state baselines. We show that:

  • Orbital circularization is a dissipative process driven by gravitational radiation and tidal friction
  • Tidal locking is an asymptotically stable state reached through dissipative evolution
  • Planetary systems settle into metastable low-dissipation configurations through dissipative processes in protoplanetary disks
  • Binary inspirals provide a natural setting for the framework’s persistence functional

The framework’s contribution is not a new mechanism of orbital evolution, but a unifying description of persistence across disparate dissipative systems using a common mathematical quantity: the persistence functional.

Keywords: attractor framework, astrophysics, gravitational radiation, tidal locking, orbital circularization, dissipative structures, Hamiltonian dynamics, planetary systems, binary inspirals, excess entropy production


1. Introduction

The attractor framework has been developed to describe persistence and change across physical, biological, cognitive, and social systems. The core claim is that every dissipative system maintains its attractor through continuous reconfiguration, and that reconfiguration generates excess entropy.

This paper extends the framework to astrophysical dissipative systems. The key insight is a distinction that is often blurred in the literature:

ConceptRole
Conservative gravitational dynamicsDefines the landscape of possible configurations (orbits, resonances, stable solutions)
Dissipative processesSelect and can stabilize particular configurations within that landscape

Gravity does not provide attractors in the dynamical systems sense — Hamiltonian systems conserve phase-space volume and do not have attractors. However, when dissipative processes are added, the system evolves toward particular asymptotically stable configurations within the family of invariant solutions. The circular orbit is not a dynamical attractor of pure Newtonian gravity; it is the endpoint of dissipative evolution (tidal friction, gravitational radiation, gas drag).

This distinction is central to the paper. Gravity defines the landscape; dissipation determines which configuration is reached.

What is new: Existing astrophysical theory explains how dissipative mechanisms drive orbital evolution. The attractor framework proposes a common mathematical quantity — the persistence functional — that measures the cumulative irreversible cost of approaching an asymptotically stable configuration. The novelty is therefore not a new mechanism of orbital evolution, but a unifying description of persistence across disparate dissipative systems.


2. Conservative vs. Dissipative Systems

2.1 Hamiltonian Dynamics

A conservative Hamiltonian system preserves phase-space volume (Liouville’s theorem). It does not have attractors in the dynamical systems sense. Orbits are determined by initial conditions and remain on their invariant tori (Arnold, 1989).

PropertyImplication
No phase-space contractionNo attractors
Time-reversibleNo arrow of time
Energy conservedNo dissipation

2.2 Dissipative Dynamics

When dissipative processes are added, the system loses energy and angular momentum. Phase-space volume contracts, and asymptotically stable states can emerge. For foundational treatments of irreversible thermodynamics, see Onsager (1931) and Prigogine (1947).

PropertyImplication
Phase-space contractionAsymptotically stable states appear
Time-irreversibleArrow of time
Energy lostEntropy generated

2.3 The Framework’s Position

The framework treats gravity as defining the landscape of possible configurations. Dissipation determines which of those configurations are actually reached.

Gravity defines the landscape. Dissipation selects the configuration.

This is the core insight of the paper.


3. The Gravitational Persistence Functional

3.1 Excess Entropy Production

Following Galida (2026c), the excess entropy production rate is defined as:σexcess(x)=σ(x)σss(x)σexcess​(x)=σ(x)−σss​(x)

where σ(x)σ(x) is the total entropy production rate and σss(x)σss​(x) is the steady-state baseline rate at the attractor.

For gravitational systems, we propose:σexcess=E˙irrevE˙ssTeffσexcess​=Teff​E˙irrev​−E˙ss​​

where E˙irrevE˙irrev​ is the total irreversible energy loss rate, E˙ssE˙ss​ is the steady-state baseline loss rate at the attractor, and TeffTeff​ is an effective temperature.

This decomposition ensures σexcess→0σexcess​→0 at the attractor, avoiding the divergence problem that would arise from integrating raw dissipation rates over infinite time. Systems that continue to dissipate at a steady baseline (e.g., a circular binary emitting GWs, a tidally locked moon with residual eccentricity-driven heating) contribute only their excess above baseline to the persistence cost.

3.2 Domain-Specific Definitions

ProcessTotal E˙E˙Baseline E˙ssE˙ssσexcessσexcess​
Orbital circularizationLGW(e)LGW​(e)LGW(e=0)LGW​(e=0)[LGW(e)LGW(0)]/Teff[LGW​(e)−LGW​(0)]/Teff​
Tidal lockingPtide(Ω,e)Ptide​(Ω,e)Ptide(Ω=n,e)Ptide​(Ω=n,e)[Ptide(Ω,e)Ptide(n,e)]/Teff[Ptide​(Ω,e)−Ptide​(n,e)]/Teff​
Disk dissipationLdiskLdisk​Ldisk, steadyLdisk, steady​[LdiskLdisk, ss]/Teff[Ldisk​−Ldisk, ss​]/Teff​

3.3 The Persistence Functional

Definition 1 (Gravitational Persistence 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: D(x)D∞​(x) measures the total excess entropy generated during the approach to an asymptotically stable configuration — the cumulative cost of reconfiguration above the steady-state baseline.

Note on gravitational wave entropy: Classical gravitational waves are coherent radiation and do not automatically carry large thermodynamic entropy. The entropy associated with gravitational wave emission arises from coarse-graining the wave’s phase space or from the generalized entropy increase of the sources (e.g., black hole horizons). The proposed definition σexcess=[LGW(e)LGW(0)]/Teffσexcess​=[LGW​(e)−LGW​(0)]/Teff​ isolates the eccentricity-specific excess above the circular-orbit baseline. Constructing an explicit entropy functional for gravitational radiation remains an open problem.


4. Orbital Circularization

4.1 The Phenomenon

Binary systems (stars, black holes, planets) often have elliptical orbits. Over time, these orbits tend to circularize — the eccentricity decreases and the orbit becomes more circular.

This is a dissipative process. The system loses energy and angular momentum through:

  • Gravitational radiation (for compact objects)
  • Tidal friction (for fluid bodies)
  • Gas drag (for protoplanetary disks)

4.2 Framework Interpretation

ComponentRole
The landscapeFamily of Keplerian orbits (all ellipses)
Asymptotically stable stateCircular orbit (endpoint of dissipative evolution)
The dissipationGravitational radiation, tidal friction, gas drag
The costσexcess=[LGW(e)LGW(0)]/Teffσexcess​=[LGW​(e)−LGW​(0)]/Teff​

The framework proposes:κ1DκD∞​1​

where κκ is the circularization rate and D=σexcessdtD∞​=∫σexcess​dt is the cumulative excess entropy production during circularization.

4.3 The Peters & Mathews Formula

The foundational computation of the gravitational-wave power from a Keplerian orbit was given by Peters & Mathews (1963). The secular decay of semi-major axis and eccentricity was derived by Peters (1964):dadt=645G3m1m2(m1+m2)c5a3(1e2)7/2(1+7324e2+3796e4)dtda​=−564​c5a3(1−e2)7/2G3m1​m2​(m1​+m2​)​(1+2473​e2+9637​e4)dedt=30415G3m1m2(m1+m2)c5a4(1e2)5/2e(1+121304e2)dtde​=−15304​c5a4(1−e2)5/2G3m1​m2​(m1​+m2​)​e(1+304121​e2)

Framework Interpretation: The decay of eccentricity e0e→0 is the approach to the asymptotically stable state. The excess entropy production is the eccentricity-dependent component of the gravitational wave luminosity:σexcess=LGW(e)LGW(0)Teffσexcess​=Teff​LGW​(e)−LGW​(0)​

This quantity vanishes as e0e→0, consistent with the ee-proportionality of the de/dtde/dt equation. Orbital eccentricity may serve as an experimentally accessible proxy for the cumulative excess entropy production.


5. Binary Inspirals

5.1 The Phenomenon

Binary systems of compact objects (neutron stars, black holes) lose energy through gravitational radiation. The orbit shrinks and the binary inspirals.

This is one of the most direct applications of the framework. The inspiral is a dissipative process driven by gravitational wave emission. For general relativistic treatments of binary dynamics and the geometry of spacetime, see Carroll (2004), Schutz (2009), Wald (1984), and Misner, Thorne & Wheeler (1973).

5.2 Framework Interpretation

ComponentRole
The landscapeFamily of binary orbits
Asymptotically stable stateQuasi-circular orbit (endpoint of circularization)
The dissipationGravitational radiation
The costσexcess=[LGW(e)LGW(0)]/Teffσexcess​=[LGW​(e)−LGW​(0)]/Teff​

5.3 The Persistence Functional

The persistence functional for a binary inspiral is:D=0σexcess(t)dt=0LGW(e(t))LGW(0)TeffdtD∞​=∫0∞​σexcess​(t)dt=∫0∞​Teff​LGW​(e(t))−LGW​(0)​dt

Note on circularization: For compact-object binaries, eccentricity damps on a much shorter timescale than the inspiral itself. Gravitational radiation circularizes the orbit well before merger, so the system reaches a quasi-circular state as a near-asymptotic limit before the final coalescence.

Hypothesis: The inspiral time ττ is inversely proportional to DD∞​:κ=1τ1Dκ=τ1​∝D∞​1​


6. Tidal Locking

6.1 The Phenomenon

Tidal locking occurs when a body’s rotational period equals its orbital period. The Moon is tidally locked to Earth. Many exoplanets in the habitable zone are expected to be tidally locked.

Tidal locking is a dissipative process. Tidal friction converts rotational energy into heat, gradually slowing the body’s rotation until it matches its orbital period.

6.2 Framework Interpretation

ComponentRole
The landscapeFamily of rotational states
Asymptotically stable stateTidal lock (rotational period = orbital period)
The dissipationTidal friction (heat generation)
The costσexcess=[Ptide(Ω,e)Ptide(Ω=n,e)]/Teffσexcess​=[Ptide​(Ω,e)−Ptide​(Ω=n,e)]/Teff​

Hypothesis: The tidally locked state is a low-dissipation configuration for the system. Once locked, tidal dissipation approaches a minimum. The excess entropy production is the despinning-specific component above whatever baseline eccentricity-driven heating persists after lock.

6.3 The Tidal Locking Timescale

The timescale for tidal locking is commonly given as (see, e.g., Murray & Dermott, 1999):τlock221Qk2mM(aR)61Ωτlock​≈212​k2​QMm​(Ra​)6Ω1​

where:

  • QQ is the tidal dissipation factor
  • k2k2​ is the Love number
  • mm is the mass of the body
  • MM is the mass of the primary
  • aa is the semi-major axis
  • RR is the radius of the body
  • ΩΩ is the rotation rate

(Different derivations use different prefactors depending on the assumed dissipation model; the (a/R)6(a/R)6 scaling is robust.)

Hypothesis: κ=1/τlockκ=1/τlock​. The recovery rate is the inverse of the locking timescale. The cumulative excess entropy production is the total tidal heat dissipated during despinning above the post-lock baseline.


7. Planetary Systems

7.1 Formation and Evolution

Planetary systems form from protoplanetary disks. The disk is a dissipative structure: it loses energy through radiation, viscosity, and accretion.

Over time, the system approaches a stable configuration:

  • Planets on nearly circular orbits
  • Resonances between orbits
  • Stable spin-orbit states

For a comprehensive treatment of solar system dynamics and tidal evolution, see Murray & Dermott (1999).

7.2 Framework Interpretation

ComponentRole
The landscapeFamily of possible planetary configurations
Metastable configurationLow-dissipation planetary system
The dissipationDisk viscosity, radiation, accretion
The costσexcess=[LdiskLdisk, ss]/Tdiskσexcess​=[Ldisk​−Ldisk, ss​]/Tdisk​

Hypothesis: Mature planetary systems approach metastable low-dissipation configurations. The cumulative excess entropy production is the total disk dissipation above the steady-state baseline integrated over the formation epoch.


8. Entropy Generation in Gravitational Systems

8.1 The Subtlety of Gravitational Entropy

Gravitational waves carry energy. Whether they carry entropy is a more subtle question. Classical gravitational waves are coherent radiation; coherent radiation is not obviously high-entropy. Binary mergers ultimately increase the generalized entropy of spacetime, but the bookkeeping is subtle.

Note: Throughout this paper, entropy generation refers to the irreversible processes associated with tidal heating, viscous dissipation, and the generalized entropy increase accompanying gravitational-wave emission. The precise entropy carried by gravitational radiation remains an active topic.

8.2 Operational Definition of σexcessσexcess​

For the purposes of this framework, we propose the following operational definition:σexcess=E˙irrevE˙ssTeffσexcess​=Teff​E˙irrev​−E˙ss​​

where:

  • E˙irrevE˙irrev​ is the total irreversible energy loss rate
  • E˙ssE˙ss​ is the steady-state baseline loss rate at the attractor
  • TeffTeff​ is an effective temperature for the dissipative process

This definition ensures σexcess0σexcess​≥0 and vanishes when the system reaches its attractor. For specific astrophysical contexts:

ContextE˙irrevE˙irrev​E˙ssE˙ssTeffTeff​
Orbital circularizationLGW(e)LGW​(e)LGW(0)LGW​(0)Effective GW temperature
Tidal lockingPtide(Ω,e)Ptide​(Ω,e)Ptide(Ω=n,e)Ptide​(Ω=n,e)Effective body temperature
Disk dissipationLdiskLdisk​Ldisk, ssLdisk, ss​Disk temperature
Black hole mergersLGWLGW​0Hawking temperature of final black hole

Note: This is a working hypothesis. Constructing an explicit entropy functional for relativistic gravitational systems remains an open problem. The effective temperature TeffTeff​ is the primary underdetermined quantity in the framework; its derivation from first principles is a priority for future work.


9. The Boundary

The framework’s boundary is not absolute zero. It is the absence of irreversible processes. At the boundary, the system becomes conservative and no entropy is generated. Hamiltonian systems exist at nonzero temperature; the boundary is dynamical, not thermal.


10. Testable Predictions

10.1 Core Prediction

Prediction: The circularization rate κκ is inversely proportional to the cumulative excess entropy production during circularization.κ1DκD∞​1​

10.2 Specific Predictions

PredictionFalsification
Tidal locking timescale correlates with total tidal heat dissipated above baselineIf no correlation, the prediction is falsified
Circularization rate correlates with total eccentricity-dependent GW energy emittedIf no correlation, the prediction is falsified
Planetary system stability correlates with total disk dissipation above steady stateIf no correlation, the prediction is falsified

11. Open Questions

QuestionStatus
Q1: Gravitational entropyWhat is the entropy of a gravitational system? (Penrose, 1965; Hawking & Ellis, 1973)
Q2: Black hole entropyHow does black hole entropy fit into the framework?
Q3: Entropy of gravitational radiationDoes gravitational radiation carry entropy, and if so, how is it defined? (Zeldovich, 1972)
Q4: Cosmological stabilityDo cosmological models admit asymptotically stable late-time solutions?
Q5: Effective temperature for GWsWhat is the correct TeffTeff​ for gravitational wave entropy production? (Galida, 2026d)
Q6: Coarse-grainingWhat coarse-graining scheme defines the entropy of classical gravitational waves? (Galida, 2026d)

12. Conclusion

The attractor framework extends naturally to astrophysical dissipative systems. The key insight is a distinction that is often blurred:

Gravity defines the landscape. Dissipation selects the configuration.

Conservative gravitational dynamics define families of stable invariant solutions. Dissipative processes — gravitational radiation, tidal friction, gas drag — select and can stabilize particular configurations within those families.

The framework does not claim that gravity provides attractors. It claims that the combination of conservative dynamics and dissipative processes produces asymptotically stable states. This is a more accurate and defensible position.

The contribution is not a new mechanism of orbital evolution, but a unifying description of persistence across disparate dissipative systems using a common mathematical quantity: the persistence functional D=σexcessdtD∞​=∫σexcess​dt, with σexcessσexcess​ operationally defined as the rate of irreversible energy loss above steady-state baseline divided by an effective temperature.


References

Arnold, V. I. (1989). Mathematical Methods of Classical Mechanics. Springer.

Carroll, S. M. (2004). Spacetime and Geometry: An Introduction to General Relativity. Addison-Wesley.

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.

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

Galida, R. (2026d). “Deep Research Questions on the Attractor Framework.” Fantasy Attractor.

Goldreich, P., & Soter, S. (1966). “Q in the Solar System.” Icarus, 5(1-6), 375-389.

Hawking, S. W., & Ellis, G. F. R. (1973). The Large Scale Structure of Space-Time. Cambridge University Press.

Misner, C. W., Thorne, K. S., & Wheeler, J. A. (1973). Gravitation. W. H. Freeman.

Murray, C. D., & Dermott, S. F. (1999). Solar System Dynamics. Cambridge University Press.

Onsager, L. (1931). “Reciprocal Relations in Irreversible Processes.” Physical Review, 37(4), 405-426.

Penrose, R. (1965). “Gravitational Collapse and Space-Time Singularities.” Physical Review Letters, 14(3), 57-59.

Peters, P. C. (1964). “Gravitational Radiation and the Motion of Two Point Masses.” Physical Review, 136(4B), B1224-B1232.

Peters, P. C., & Mathews, J. (1963). “Gravitational Radiation from Point Masses in a Keplerian Orbit.” Physical Review, 131(1), 435-440.

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

Schutz, B. F. (2009). A First Course in General Relativity (2nd ed.). Cambridge University Press.

Wald, R. M. (1984). General Relativity. University of Chicago Press.

Zeldovich, Y. B. (1972). “A Hypothesis Unifying the Structure and the Entropy of the Universe.” Monthly Notices of the Royal Astronomical Society, 160(1), 1P-4P.


Suggested citation: Galida, R. S. (2026). The Attractor Framework in Astrophysics: Persistence, Entropy, and Gravitational Systems (Final Edition). Fantasy Attractor.

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

DomainEntropy FunctionalBaseline σssσssExcess σexcessσexcess​Recovery Rate κκ
PhysicalThermodynamic entropy0 (equilibrium)S˙S˙infδΔSinfΔSδ
BiologicalMetabolic entropyResting metabolic rateMetabolic rate — restinginfδσexcessdtinf∫σexcess​dtδ
CognitiveFree energyBaseline neural dissipationF˙F˙ssF˙−F˙ssinfδσexcessdtinf∫σexcess​dtδ
SocialSocial entropy productionSteady-state social dissipationσsocialσssσsocial−σssinfδσexcessdtinf∫σexcess​dtδ

6.2 The Universal Structure

Every domain follows the same mathematical structure:

ComponentExpression
Excess entropy productionσexcess(x)=σ(x)σssσexcess​(x)=σ(x)−σss
Cumulative costD(x)=0σexcess(ϕt(x))dtD∞​(x)=∫0∞​σexcess​(ϕt​(x))dt
Recovery rateκ=infxδ(x)/D(x)κ=infxδ(x)/D∞​(x)
Basin depthB=D(saddle)B=D∞​(saddle)
Transport equationDf=σ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.

DomainAttractorEntropy Generation at Attractor
PhysicalEquilibriumσ=0σ=0
BiologicalHomeostasisσ=σss>0σ=σss​>0 (resting metabolism)
CognitiveSettled Beliefσ=σss>0σ=σss​>0 (baseline neural dissipation)
SocialCoordinated 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

DomainPredictionFalsification
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

QuestionStatusDifficulty
Q1: Uniqueness of S(x)S(x)Are there multiple valid entropy functionals for a given domain?Hard
Q2: Variational principleIs there a universal variational principle that yields S(x)S(x)?Hard
Q3: Social second lawDoes σsocial0σsocial≥0 always hold during recovery?Very Hard
Q4: Cross-level entropyHow does entropy generation at one level relate to entropy generation at another?Hard
Q5: MeasurementCan we measure excess entropy generation in cognitive and social systems directly?Moderate
Q6: UnificationCan 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):

PropertyStatement
Non-negativityDT(x)0DT​(x)≥0
MonotonicityDT2(x)DT1(x)DT2​​(x)≥DT1​​(x) for T2T1T2​≥T1​
AdditivityDT+S(x)=DT(x)+DS(ϕT(x))DT+S​(x)=DT​(x)+DS​(ϕT​(x))
Instantaneous growthddTDT(x)=δ(ϕT(x))dTdDT​(x)=δ(ϕT​(x))
Occupation measureDT(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

QuestionStatusDifficulty
Q1: Nonlinear systemsDoes infδDinfD∞​δ​ equal the local Lyapunov exponent?Hard
Q2: Local vs. global consistencyDoes limxAδ(x)D(x)=κlimx→A​D∞​(x)δ(x)​=κ hold for general nonlinear systems?Hard
Q3: Non-normal systemsDoes the infimum equal the slowest eigenvalue for non-normal AA?Moderate
Q4: Multiple timescalesDoes the infimum isolate the slowest timescale?Hard
Q5: Stochastic systemsHow does noise affect the finite-horizon estimator?Hard
Q6: Multiple attractorsHow 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.

The Persistence Functional: A Mathematical Measure of Attractor Resilience

Author: Robert Galida
Date: May 2026


Abstract

The attractor framework says that persistence under disturbance is the basic mark of reality.
To turn this idea into a formal science, we introduce the persistence functional P(x)P(x).

P(x)P(x) is a single number that measures:

  • How deep a state is inside an attractor basin.
  • How quickly it returns after a knock.

We define P(x)P(x) for three different kinds of systems:

  1. Deterministic dissipative systems – here PP is linked to Lyapunov exponents and basin stability.
  2. Stochastic systems – here PP is linked to escape time and quasipotential.
  3. Information‑theoretic systems – here PP is linked to negative free energy or mutual information.

The recovery rate P˙/PP˙/P is a universal sign of critical slowing down – a warning that a system is about to tip.

We also discuss limitations: resilience may depend on direction (“anisotropic”), and multiple timescales may need vector or tensor persistence. We list open mathematical problems.
This paper is a roadmap, not a finished theory.


1. Introduction

In the attractor framework, persistence under disturbance is central. But we have not had a single number to say how persistent a state is.
The persistence functional P(x)P(x) aims to fill that gap.

What P(x)P(x) should do:

  • P(x)>0P(x)>0 for states inside an attractor basin.
  • For a conservative attractor (like a free electron), PP is maximal (normalised to 1).
  • For a dissipative attractorPP drops after a disturbance and then recovers.
    The recovery rate P˙/PP˙/P equals:
    • the negative of the largest Lyapunov exponent (for deterministic systems)
    • the inverse return time (for stochastic systems)
    • the rate of information loss (for informational systems)
  • PP falls as the system approaches a bifurcation, giving early warning.

We do not give one universal formula. Instead, we give a family of definitions, each suited to a different type of system, all united by the same purpose – measuring resilience.


2. Deterministic Dissipative Systems

Consider a smooth system x˙=f(x)x˙=f(x) with a stable attractor AA and its basin B(A)B(A).

A natural candidate for P(x)P(x) uses a Lyapunov function V(x)V(x) – a kind of energy that always decreases inside the basin (V˙<0V˙<0).

We define:P(x)=1V(x)VAVmaxVAP(x)=1−Vmax​−VAV(x)−VA​​

This gives P=1P=1 on the attractor and P0P→0 at the basin boundary.
Near the attractor, the recovery rate is related to the largest Lyapunov exponent λ1λ1​:P˙/Pλ1P˙/P≈−λ1​

When the system approaches a tipping point, λ10λ1​→0−, so the recovery rate slows down – this is critical slowing down.

Conclusion: For deterministic systems, PP can be built from a Lyapunov function. The recovery rate equals the negative of the largest Lyapunov exponent.


3. Stochastic Systems

When noise is present, persistence is about how long it takes to escape from the basin.
The mean first passage time τ(x)τ(x) – the average time to leave – is a natural measure.

We define:P(x)=τ(x)τmaxP(x)=τmax​τ(x)​

where τmaxτmax​ is the value at the attractor.

For weak noise, τ(x)τ(x) grows exponentially with the quasipotential U(x)U(x) (Freidlin–Wentzell theory):τ(x)eU(x)/ϵτ(x)∼eU(x)/ϵ

So:P(x)e(UmaxU(x))/ϵP(x)∝e−(Umax​−U(x))/ϵ

The recovery rate is the inverse of the return time. As a tipping point is approached, the return time diverges, and the recovery rate goes to zero. This again gives critical slowing down – rising variance and autocorrelation.

Conclusion: For stochastic systems, PP is proportional to the mean exit time (or the exponential of the quasipotential). This connects persistence to large deviation theory.


4. Information‑Theoretic Systems

For systems where information matters (neural, cognitive, social), we can define persistence using mutual information between past and future.

Let Ipast,futureIpast,future​ be the predictive information. Then:P(t)=I(past;future at time t)orP=esurprisalP(t)=I(past;future at time t)orP=e−surprisal

The decay of P(t)P(t) over time measures memory loss.
Landauer’s principle connects information loss to entropy production:P˙/PS˙kBln2P˙/P≤−kB​ln2S˙​

Alternatively, in the free energy principle (Friston), the negative free energy FF acts like a Lyapunov function. We can set:P=eF/kTorP=FP=eF/kTorP=−F

Then P˙/PP˙/P is the rate of free energy minimisation, which slows near bifurcations.

Conclusion: For information‑theoretic systems, PP can be defined via mutual information decay or negative free energy, linking persistence to entropy production and predictive coding.


5. Unifying Recovery Rate and Critical Slowing Down

Across all types of systems, the recovery rate λrec=P˙/Pλrec​=−P˙/P (just after a small disturbance) is a universal indicator:

  • Deterministic dissipative: λrec=λ1λrec​=∣λ1​∣ (absolute value of the largest Lyapunov exponent)
  • Stochastic: λrecλrec​ = inverse of the return time, related to the quasipotential’s curvature
  • Information‑theoretic: λrecλrec​ = rate of free energy minimisation or information loss

As the system approaches a bifurcation, λrec0λrec​→0. This is critical slowing down.
It shows up as rising lag‑1 autocorrelation and variance (Scheffer et al., 2009).
So PP and its recovery rate give early warnings.


6. Normalisation for Conservative Attractors

For a perfect conservative attractor (e.g., an electron in its ground state, no decay), the persistence functional should be constant and maximal:Pcons=1for all timesPcons​=1for all times

No recovery rate is defined (or it is zero). This anchors the scale.

For emergent approximate conservative systems (like atomic clocks), PP is very close to 1 and decays extremely slowly.


7. Limitations – Scalar Collapse and Anisotropic Resilience

A single scalar P(x)P(x) may not be enough for systems where resilience is anisotropic – that is, recovery speed depends on the direction of the perturbation.
High‑dimensional systems can have multiple timescales (fast and slow modes). A scalar average can miss important structure.

Future work may need:

  • Vector persistence – a list of recovery rates along different directions.
  • Tensor persistence – a metric that captures the full shape of the basin.
  • Persistence manifold – the geometry of the basin in state space.

We accept this limitation. The scalar PP is a useful first approximation for systems with isotropic resilience or for early‑warning applications where a single number is enough. For complex systems, a multidimensional generalisation is an open research problem.


8. Open Mathematical Problems

  1. Derive P(x)P(x) from first principles for a given class of systems (e.g., from a variational principle).
  2. Prove that −P˙/P=∣λ1∣−P˙/P=∣λ1​∣ for a wide class of dissipative systems.
  3. Extend the definition to systems with multiple attractors and chaotic basins (where basin stability is fractal).
  4. Establish a rigorous relationship between PP and the mutual information decay rate for non‑equilibrium processes.
  5. Formulate a universal persistence functional that works across all regimes – or prove it’s impossible.
  6. Test the predictive power of PP in controlled experiments (e.g., ecological microcosms, neural cultures, social media sentiment).
  7. Develop vector/tensor persistence for anisotropic resilience.

9. Conclusion

The persistence functional P(x)P(x) gives a mathematical language for attractor resilience.

We have given operational definitions for three regimes:

  • Deterministic dissipative → Lyapunov / basin stability
  • Stochastic → escape time / quasipotential
  • Information‑theoretic → mutual information / free energy

The recovery rate P˙/PP˙/P unifies critical slowing down across all these domains.

We have explicitly noted limitations (scalar collapse, anisotropy) as open problems.

This paper is a roadmap, not a final theory. The framework now has a quantitative step.


Suggested citation: Galida, R. S. (2026). The Persistence Functional: Towards a Mathematical Measure of Attractor Resilience (Reader‑Friendly Version). Fantasy Attractor.

image_pdfimage_print