Home » Posts tagged 'free energy principle'
Tag Archives: free energy principle
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.
Free Will as Attractor Autonomy: A Dynamical Account of Agency
Author: Robert Galida https://fantasyattractor.com/
Date: May 2026
Abstract
Free will is often seen as either a magical mystery (libertarianism) or an illusion (hard determinism).
This paper offers a third view using the attractor framework.
In this framework, your mind is a dissipative, self‑referential attractor of your whole body.
Free will is redefined as attractor autonomy:
- The ability to generate behaviour from your own internal dynamics.
- To keep yourself stable over time.
- To model yourself.
- And to reshape your own attractor landscape over time.
Agency comes in degrees – it is not a simple yes/no.
We give a mathematical formula for an agency index A that combines three factors:
- Attractor dimensionality D (complexity of your brain’s activity)
- Recursive self‑modification R (your ability to change your own habits)
- Self‑reference strength S (how well you have a persistent self‑model)
The paper makes a falsifiable prediction: an inverted‑U relationship between attractor dimensionality and sense of agency – too low or too high reduces agency.
We describe how to test this with EEG, intentional binding tasks, and statistical methods. We also engage with classic compatibilist philosophers (Frankfurt, Dennett) and address Pereboom’s manipulation argument.
We even provide an explicit rule to avoid the “liver problem” (a false positive for self‑reference).
1. Introduction
The attractor framework says that persistence under disturbance is the basic mark of reality.
Minds are dissipative attractors – patterns that need constant energy flow, integrating the whole body.
In this view, free will cannot be a supernatural break from cause and effect. Instead, it must be a dynamical property of certain attractors.
We do not claim to solve the ancient free will debate. We offer a naturalistic, testable redefinition that adds new empirical content to compatibilism.
2. What Free Will Is Not – And What It Is
2.1 Rejecting supernatural libertarianism
Libertarian free will requires an uncaused choice – a break in the chain of cause and effect.
The attractor framework rejects this: there is no evidence for it, and it contradicts physical laws.
2.2 The error of hard determinism
Hard determinism says freedom is an illusion because everything is determined. But it confuses “determined” with “externally coerced”.
A system can be internally determined – by its own attractor – yet still be free. That is the core of compatibilism.
2.3 Free will as attractor autonomy
We define free will (or agency) as the degree to which a system has four properties:
- Dissipative persistence – it stays alive by using energy and exporting waste (measured by energy use and recovery speed).
- Self‑reference – it has an internal subsystem (an “indexical locus”) that models the whole system and is stable.
- Trajectory selection – it can choose among different possible futures (measured by policy entropy H(π)).
- Recursive self‑engineering – it can change its own attractor shape (measured by learning‑to‑learn or metacognitive accuracy).
These four are jointly necessary. If any is missing, agency is at best primitive.
Because they are necessary, we combine them with a multiplicative formula (if any factor is zero, agency is zero).A=(Dmax−DminD−Dmin)α(RmaxR)β(Smax−SminS−Smin)γ
Where:
- D = attractor dimensionality (e.g., from EEG)
- R = recursive modification capacity (e.g., improvement in a meta‑learning task)
- S = self‑reference strength (normalised mutual information)
The constants (Dmin,Dmax, etc.) are set from a reference population.
The exponents α,β,γ are estimated from data (e.g., comparing healthy people with patients).
A threshold Acrit (e.g., the 5th percentile of healthy humans) decides where agency begins.
Agency is graded:
- Rock: A≈0
- Thermostat: A≈0
- Worm: A≈0.1 (some learning, little self‑model)
- Human: A≈0.8
3. The Indexical Locus: Defining the “Self” and Avoiding the “Liver Problem”
The indexical locus L is the part of the system that acts as a persistent self‑model.
To avoid trivial cases (like a liver having high mutual information with the rest of the body), we add three extra conditions:
- Top‑down causal influence – L can change the rest of the body in ways that serve the body’s goals (measured by variance explained beyond bottom‑up effects).
- Informational closure – L’s own dynamics are relatively independent of the rest over short timescales (conditional mutual information > 0).
- Self‑referential loop – L influences the body, and the body influences L back (bidirectional Granger causality).
These criteria rule out livers, pacemakers, and simple homeostats. The indexical locus is a recursive self‑model, not just a predictive subsystem.
4. Active Inference and Policy Entropy
In active inference (Friston), agents try to minimise “free energy” – they pick policies (sequences of actions).
Each policy is a trajectory through the agent’s attractor landscape.
Policy entropy H(π)=−∑p(π)logp(π) measures how many different policies are available.
- Low entropy → rigid, one‑track mind.
- High entropy → flexible, but possibly noisy.
Free will is the ability to access many low‑energy policies. The agent’s choices are not random; they are constrained by the attractor geometry. But if several attractor basins are open, the agent can choose among them – that is what we feel as free choice.
Policy entropy can be measured in behavioural tasks where multiple choices are equally good (e.g., probabilistic reversal learning, two‑armed bandit tasks).
5. The Inverted‑U Prediction and Falsification
5.1 Core prediction
We predict an inverted‑U relationship between attractor dimensionality D and the subjective sense of agency (e.g., from intentional binding experiments).
- Very low D → chaotic, unstable (like schizophrenia) → low agency.
- Very high D → rigid, stuck (like OCD) → low agency.
- In the middle → flexible but stable → high agency.
The agency index A also includes R and S, which we think increase agency across the board. So to test the inverted‑U for D alone, you need to control for R and S (e.g., study people matched on those, or use partial correlation).
5.2 How to measure and test
- Attractor dimensionality DD – use the Grassberger‑Procaccia algorithm on 5‑min resting‑state EEG/MEG.
- Sense of agency – use the intentional binding paradigm: press a key, then a tone sounds; participants estimate the time between action and tone. Stronger binding means higher agency.
- Statistical test – fit a quadratic regression: agency = β0+β1D+β2D2.
If β2<0 and the vertex lies inside the observed range of D, the inverted‑U is supported. Use bootstrap (1000 resamples) to check confidence intervals.
5.3 Falsification condition
The framework is falsified if:
- The quadratic coefficient β2 is not negative (no inverted‑U).
- Or, in a clinical experiment (e.g., increasing D in OCD patients with NMDA drugs), agency does not decrease but keeps increasing.
6. Experimental Proxies – Summary Table
| Construct | Measure | How to record | Expected relation to agency |
|---|---|---|---|
| Attractor dimensionality D | Correlation dimension (Grassberger‑Procaccia) | Resting‑state EEG/MEG (5 min) | Inverted‑U |
| Policy entropy H(π) | Entropy of choice distribution | Probabilistic reversal learning (200 trials) | Inverted‑U |
| Sense of agency | Intentional binding magnitude | Action‑outcome interval compression (50 trials) | Max at intermediate D |
| Recursive self‑modification R | Learning‑to‑learn improvement | Meta‑learning task (pre‑post difference) | Positive (more is better) |
| Self‑reference strength S | Normalised mutual info In(L;S) | Resting‑state fMRI or MEG | Threshold > θ |
7. Hierarchical Constraints and Social Attractors
Free will is nested inside larger attractors – society, culture, laws, economy. Your range of choices is partly set by these.
This is not an objection; it is just the fact that freedom is always constrained autonomy.
We predict that societies with more cultural diversity (higher “cultural entropy”) allow more individual agency, other things being equal. This can be tested by cross‑cultural comparisons of policy entropy in decision tasks.
8. Engagement with Compatibilist Literature
8.1 Standard compatibilists (Frankfurt, Dennett)
- Frankfurt (1971): freedom is about your will aligning with your own desires. Our framework adds that those desires must be encoded in a persistent self‑referential attractor. The recursive self‑engineering component R maps directly to Frankfurt’s “second‑order volitions”.
- Dennett (1984): freedom is about being able to respond to reasons. Our framework adds that this requires a certain basin geometry and recursive plasticity.
8.2 Addressing Pereboom’s manipulation argument
Pereboom argues: if a neuroscientist engineers your brain, you are not free – even if your behaviour comes from internal dynamics.
Our reply: agency requires recursive self‑modification (R>0) at some point in your history.
- A perfectly manipulated agent that never changed its own attractor would have R≈0 and thus A≈0.
- A healthy human who learned and adapted has R>0 and genuine agency.
The origin of the initial attractor does not matter – only the presence of self‑modification over time.
9. Open Questions and Limitations
- Calibrating exponents – α,β,γ and the threshold θ need to be estimated from large‑scale data (e.g., Human Connectome Project) using maximum likelihood.
- The liver problem – our exclusion criteria need empirical validation; we must show that organs like the liver do not satisfy them.
- Inverted‑U for policy entropy – the same shape is predicted but may be hidden by decision noise.
- Moral responsibility – the framework gives a basis for responsibility (if A>Acrit), but it does not settle all normative questions – it only gives a scientific starting point.
10. Conclusion
Free will is not a supernatural escape from physics. It is a dynamical property of certain dissipative, self‑referential attractors:
- The ability to act from your own internal dynamics.
- To keep a stable self‑model over time.
- And to reshape your own attractor landscape.
This account is compatibilist, testable, and graded.
The inverted‑U prediction, with a specified statistical test, gives a clear falsification criterion.
The dance of free will is the dance of a self that persists under perturbation.
Suggested citation: Galida, R. S. (2026). Free Will as Attractor Autonomy: A Dynamical Account of Agency in the Attractor Framework (Reader‑Friendly Version). 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) 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) for three different kinds of systems:
- Deterministic dissipative systems – here P is linked to Lyapunov exponents and basin stability.
- Stochastic systems – here P is linked to escape time and quasipotential.
- Information‑theoretic systems – here P is linked to negative free energy or mutual information.
The recovery rate −P˙/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) aims to fill that gap.
What P(x) should do:
- P(x)>0 for states inside an attractor basin.
- For a conservative attractor (like a free electron), P is maximal (normalised to 1).
- For a dissipative attractor, P drops after a disturbance and then recovers.
The recovery rate −P˙/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)
- P 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) with a stable attractor A and its basin B(A).
A natural candidate for P(x) uses a Lyapunov function V(x) – a kind of energy that always decreases inside the basin (V˙<0).
We define:P(x)=1−Vmax−VAV(x)−VA
This gives P=1 on the attractor and P→0 at the basin boundary.
Near the attractor, the recovery rate is related to the largest Lyapunov exponent λ1:−P˙/P≈−λ1
When the system approaches a tipping point, λ1→0−, so the recovery rate slows down – this is critical slowing down.
Conclusion: For deterministic systems, P 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) – the average time to leave – is a natural measure.
We define:P(x)=τmaxτ(x)
where τmax is the value at the attractor.
For weak noise, τ(x) grows exponentially with the quasipotential U(x) (Freidlin–Wentzell theory):τ(x)∼eU(x)/ϵ
So: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, P 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,future be the predictive information. Then:P(t)=I(past;future at time t)orP=e−surprisal
The decay of P(t) over time measures memory loss.
Landauer’s principle connects information loss to entropy production:P˙/P≤−kBln2S˙
Alternatively, in the free energy principle (Friston), the negative free energy −F acts like a Lyapunov function. We can set:P=e−F/kTorP=−F
Then −P˙/P is the rate of free energy minimisation, which slows near bifurcations.
Conclusion: For information‑theoretic systems, P 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 (just after a small disturbance) is a universal indicator:
- Deterministic dissipative: λrec=∣λ1∣ (absolute value of the largest Lyapunov exponent)
- Stochastic: λrec = inverse of the return time, related to the quasipotential’s curvature
- Information‑theoretic: λrec = rate of free energy minimisation or information loss
As the system approaches a bifurcation, λrec→0. This is critical slowing down.
It shows up as rising lag‑1 autocorrelation and variance (Scheffer et al., 2009).
So P 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 times
No recovery rate is defined (or it is zero). This anchors the scale.
For emergent approximate conservative systems (like atomic clocks), P is very close to 1 and decays extremely slowly.
7. Limitations – Scalar Collapse and Anisotropic Resilience
A single scalar 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 P 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
- Derive P(x)P(x) from first principles for a given class of systems (e.g., from a variational principle).
- Prove that −P˙/P=∣λ1∣−P˙/P=∣λ1∣ for a wide class of dissipative systems.
- Extend the definition to systems with multiple attractors and chaotic basins (where basin stability is fractal).
- Establish a rigorous relationship between PP and the mutual information decay rate for non‑equilibrium processes.
- Formulate a universal persistence functional that works across all regimes – or prove it’s impossible.
- Test the predictive power of PP in controlled experiments (e.g., ecological microcosms, neural cultures, social media sentiment).
- Develop vector/tensor persistence for anisotropic resilience.
9. Conclusion
The persistence functional 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˙/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.

