Remarkable similarities in distributions of dynamical observables in chaotic systems
Pith reviewed 2026-05-22 15:39 UTC · model grok-4.3
The pith
Dynamical observables in chaotic systems share the same large-deviation rate function when their defining functions differ by a derived one
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We find that certain dynamical observables exhibit a remarkable statistical similarity: even when constructed with distinct functions g1(x) and g2(x), different observables are described by the same rate function. We provide a physical interpretation for this striking similarity by showing that g1(x)−g2(x) belongs to a class of functions that we call derived. Furthermore, we show that if g(x) itself is derived, then the distribution of A becomes independent of N in the large-N limit, and is generally non-Gaussian although mirror-symmetric. We demonstrate that the position observable for certain open maps and the finite-time Lyapunov exponent for the logistic map are of this derived form.
What carries the argument
The class of derived functions, whose defining property forces the large-deviation rate function of an observable to be unchanged when g is replaced by g plus a derived function.
If this is right
- Observables whose g functions differ by a derived function obey the identical large-deviation principle.
- When g is itself derived the probability distribution of the observable becomes independent of trajectory length N for large N.
- Such distributions remain mirror-symmetric yet are typically non-Gaussian.
- The derived property accounts for known statistical similarities between position observables in open maps and finite-time Lyapunov exponents in the logistic map.
Where Pith is reading between the lines
- The derived-function classification may place many additional observables into equivalence classes sharing the same rate function.
- The same structural argument could be checked in other chaotic maps or continuous-time flows to identify further derived relations.
- The unification supplies a route to simplify large-deviation calculations by reducing them to the simplest representative of each derived class.
Load-bearing premise
That the difference between the functions defining two observables belongs to the class of derived functions, without which the rate functions would generally differ.
What would settle it
Numerical extraction of the empirical rate functions for two observables whose g functions differ by a non-derived quantity; identical rates despite a non-derived difference would falsify the necessity of the derived class.
Figures
read the original abstract
The study of chaotic systems, where rare events play a pivotal role, is essential for understanding complex dynamics due to their sensitivity to initial conditions. Recently, tools from large deviation theory, typically applied in the context of stochastic processes, have been used in the study of chaotic systems. Here, we study dynamical observables, $A = \sum_{n=1}^N g(\textbf{x}_n)$, defined along a chaotic trajectory $\{\textbf{x}_1, \textbf{x}_2, \ldots, \textbf{x}_N\}$. For most choices of $g(\textbf{x})$, $A$ satisfies a central limit theorem: At large sequence size $N \gg 1$, typical fluctuations of $A$ follow a Gaussian distribution with a variance that scales linearly with $N$. Large deviations of $A$ are usually described by the large deviation principle, that is, $P(A) \sim e^{- N I(A/N)}$, where $I(a)$ is the rate function. We find that certain dynamical observables exhibit a remarkable statistical similarity: even when constructed with distinct functions $g_1(\textbf{x})$ and $g_2(\textbf{x})$, different observables are described by the same rate function. We provide a physical interpretation for this striking similarity by showing that $g_1(\textbf{x})-g_2(\textbf{x})$ belongs to a class of functions that we call ``derived''. Furthermore, we show that if $g(\textbf{x})$ itself is ``derived'', then the distribution of $A$ becomes independent of $N$ in the large-$N$ limit, and is generally non-Gaussian (although it is mirror-symmetric). We demonstrate that the position observable for certain open maps, used to model random walks and the finite-time Lyapunov exponent (FTLE) for the logistic map are of this derived form, thus providing a simple explanation for some existing results.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies dynamical observables A = sum_{n=1}^N g(x_n) along trajectories of chaotic maps. It claims that distinct observables built from g1 and g2 share the identical large-deviation rate function I(a) whenever g1(x) − g2(x) belongs to a newly introduced class of 'derived' functions. When g itself is derived, the distribution of A becomes N-independent for large N and is non-Gaussian yet mirror-symmetric. Concrete examples are the position observable on certain open maps (modeling random walks) and the finite-time Lyapunov exponent on the logistic map; these are asserted to be derived, thereby explaining prior observations of non-Gaussian or N-independent statistics.
Significance. If the definition and closure properties of the derived class are made rigorous, the result supplies a unifying physical mechanism for the observed coincidence of rate functions across different observables and accounts for the N-independent non-Gaussian limits reported in the literature for open maps and the logistic map. The work thereby strengthens the link between large-deviation theory and deterministic chaos without introducing free parameters or ad-hoc fitting.
major comments (2)
- [§3] §3 (definition of derived functions): the manuscript introduces the class of 'derived' functions to establish both the equality of rate functions I(a) for g1 and g2 and the N-independent limit when g is derived, yet supplies neither an explicit functional equation (e.g., g = f ∘ T − f or a cocycle condition) nor a proof that the class is closed under addition and composition with the observables needed for the large-deviation principle. Without these, the central equivalence does not follow from standard LDT results for deterministic systems.
- [§4.2] §4.2 (logistic-map FTLE example): membership of the finite-time Lyapunov exponent in the derived class is asserted rather than verified by direct substitution into the map's transfer operator or by exhibiting the requisite f such that g = f ∘ T − f. Because this example is used to illustrate the N-independent non-Gaussian claim, an explicit check is required to confirm that the asserted rate-function identity is not circular.
minor comments (2)
- [Abstract] The abstract and introduction use the phrase 'derived' before its formal definition appears; a forward reference or brief parenthetical characterization would improve readability.
- [§2] Notation for the observable A and the scaled variable a = A/N is introduced inconsistently across sections; a single consolidated definition would reduce ambiguity.
Simulated Author's Rebuttal
We thank the referee for the thorough review and valuable suggestions. We have revised the manuscript to address the concerns regarding the rigor of the derived functions definition and the verification of the example. Below we respond point by point.
read point-by-point responses
-
Referee: [§3] §3 (definition of derived functions): the manuscript introduces the class of 'derived' functions to establish both the equality of rate functions I(a) for g1 and g2 and the N-independent limit when g is derived, yet supplies neither an explicit functional equation (e.g., g = f ∘ T − f or a cocycle condition) nor a proof that the class is closed under addition and composition with the observables needed for the large-deviation principle. Without these, the central equivalence does not follow from standard LDT results for deterministic systems.
Authors: We appreciate the referee's observation that the definition of 'derived' functions requires more rigor. In the revised manuscript, we now provide an explicit functional equation: a function g is derived if there exists a function f such that g(x) = f(T(x)) - f(x), where T is the dynamical map. We also prove that this class is closed under addition and under composition with observables in the context of the large deviation principle. This allows the equality of rate functions to follow directly from standard results in large deviation theory for deterministic dynamical systems. revision: yes
-
Referee: [§4.2] §4.2 (logistic-map FTLE example): membership of the finite-time Lyapunov exponent in the derived class is asserted rather than verified by direct substitution into the map's transfer operator or by exhibiting the requisite f such that g = f ∘ T − f. Because this example is used to illustrate the N-independent non-Gaussian claim, an explicit check is required to confirm that the asserted rate-function identity is not circular.
Authors: We agree that an explicit verification for the finite-time Lyapunov exponent (FTLE) is essential to avoid any appearance of circularity. In the revised version, we explicitly construct the function f for the logistic map such that the FTLE observable g satisfies g = f ∘ T - f. We verify this by direct substitution and confirm the resulting N-independent distribution. This strengthens the illustration of the non-Gaussian, mirror-symmetric limit. revision: yes
Circularity Check
No significant circularity; new 'derived' class defined and applied independently
full rationale
The paper introduces the 'derived' function class explicitly to interpret why rate functions coincide for observables differing by such a function, and separately shows that when g itself is derived the large-N distribution becomes N-independent and non-Gaussian. No equation or step in the abstract or described derivation reduces the claimed equivalence to a self-definition, a fitted parameter renamed as prediction, or a load-bearing self-citation. The examples (position observable on open maps, FTLE on logistic map) are asserted to belong to the class after the definition is given, rather than the membership being verified only by assuming the rate-function identity. The derivation chain therefore remains self-contained against external large-deviation results and does not collapse to its inputs by construction.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The large deviation principle applies to the dynamical observable A in chaotic systems.
invented entities (1)
-
derived functions
no independent evidence
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We provide a physical interpretation for this striking similarity by showing that g1(x)−g2(x) belongs to a class of functions that we call “derived”. Furthermore, we show that if g(x) itself is “derived”, then the distribution of A becomes independent of N in the large-N limit
-
IndisputableMonolith/Foundation/AbsoluteFloorClosure.leanabsolute_floor_iff_bare_distinguishability unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
If g∗ satisfies the third property (for all periodic orbits), then, under certain regularity conditions ... it follows from Livshits’s theorem that g∗ is a derived function.
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
If g∗(x) is derived, then its ensemble-average value on the IM vanishes. This may be seen by a direct calcu- lation, ⟨g∗(x)⟩S = Z dxg∗(x)ps(x) = Z dxh(x)ps(x) − Z dxh (f (x)) ps(x) . (29) By changing the integration variable in the second inte- gral on the right-hand side to f (x), it is straightforward to show that the result coincides with the first int...
-
[2]
This may be seen through the Green-Kubo formula (7)
The variance of A∗ grows sublinearly with N. This may be seen through the Green-Kubo formula (7). In- deed, for a derived observable function, one has ⟨g∗(x0)g∗(xδ)⟩S = ⟨h(x0)h(xδ)⟩S − ⟨h(x0)h(xδ+1)⟩S − ⟨ h(x1)h(xδ)⟩S + ⟨h(x1)h(xδ+1)⟩S = 2 cδ − cδ−1 − cδ+1 (30) where cδ = ⟨h(x0)h(xδ)⟩S , (31) and as a result, the sum in (7) vanishes: ∞X δ=−∞ ⟨g∗ (x0) g∗ (...
-
[3]
The average of g∗(x) over any periodic orbit van- ishes. Indeed, if x1, . . . ,xp is a periodic orbit of period p (so f(xp) = x1) then pX i=1 g∗ (xi) = pX i=1 h (xi) − pX i=1 h (f (xi)) = 0 . (33) In the converse direction, if g∗ satisfies the third prop- erty (for all periodic orbits), then, under certain regular- ity conditions on the functions g∗ and f...
-
[4]
which is consistent with the vanishing result of the Green-Kubo formula, Eq. (32). In Sec. III, for the 1d tent and logistic maps, we have used the first and second powers of x as observable func- tions as seen in Eq. (16). Following from those examples, let us now analyze the distribution of the observable gen- erated by the derived function g∗(x) which ...
-
[5]
of ˜f. One can similarly show that I(1/2) = ln φ, which corresponds to an initial condition that is exponentially close to the other fixed point, 2 /(3 + √ 5). The integer part of the RW’s position, A∗, may be described by the observable functiongpos(x) = gright(x)− gleft(x), which equals +1 (-1) when the particle hops to the right (left), as shown in pan...
-
[6]
Quasi-analytic continuation of the power series The SCGF λ(k) can be obtained through a pertur- bative solution of the tilted eigenvalue problem in the small-k limit [26]. This is followed by a quasi-analytic continuation that pushes the result beyond the conver- gence radius for the original expansion, which we intro- 13 duce here. The idea is to write t...
-
[7]
Power method for calculating λ(k) Since the analytical approach requires us to calculate the series expansion to very high terms, we also present an alternative direct semi-analytic route to the largest eigenvalue of the tilted operator Lk, Eq. (12). We will make use of the standard power method for any value of k. Both approaches yield identical results ...
work page 2000
-
[8]
L. De Cruz, S. Schubert, J. Demaeyer, V. Lucarini, and S. Vannitsem, Exploring the Lyapunov instability proper- ties of high-dimensional atmospheric and climate models, Nonlin. Processes Geophys. 25, 387 (2018)
work page 2018
-
[9]
I. Klioutchnikov, M. Sigova, and N. Beizerov, Chaos The- ory in Finance, Procedia Computer Science 119, 368 (2017)
work page 2017
-
[10]
F. Roy, G. Biroli, G. Bunin, and C. Cammarota, Nu- merical implementation of dynamical mean field the- ory for disordered systems: Application to the Lotka– Volterra model of ecosystems, J. Phys. A: Math. Theor. 52, 484001 (2019)
work page 2019
-
[11]
F. Roy, M. Barbier, G. Biroli, and G. Bunin, Complex interactions can create persistent fluctuations in high- diversity ecosystems, PLoS Comput Biol 16, e1007827 (2020)
work page 2020
-
[12]
M. T. Pearce, A. Agarwala, and D. S. Fisher, Stabiliza- tion of extensive fine-scale diversity by ecologically driven spatiotemporal chaos, Proc. Natl. Acad. Sci. U.S.A. 117, 14572 (2020)
work page 2020
-
[13]
O. Postavaru, S. Anton, and A. Toma, COVID-19 pan- demic and chaos theory, Mathematics and Computers in Simulation 181, 138 (2021)
work page 2021
-
[14]
T. Geisel and J. Nierwetberg, Onset of Diffusion and Uni- versal Scaling in Chaotic Systems, Phys. Rev. Lett. 48, 7 (1982)
work page 1982
- [15]
-
[17]
J.-P. Eckmann and I. Procaccia, Fluctuations of dynam- ical scaling indices in nonlinear systems, Phys. Rev. A 34, 659 (1986)
work page 1986
-
[18]
Y. Sato and R. Klages, Anomalous Diffusion in Ran- dom Dynamical Systems, Phys. Rev. Lett. 122, 174101 (2019)
work page 2019
- [19]
-
[20]
J. Wouters and F. Bouchet, Rare event computation in deterministic chaotic systems using genealogical particle analysis, J. Phys. A: Math. Theor. 49, 374002 (2016)
work page 2016
- [21]
-
[22]
F. Ragone and F. Bouchet, Rare Event Algorithm Study of Extreme Warm Summers and Heatwaves Over Eu- rope, Geophysical Research Letters 48, e2020GL091197 (2021)
work page 2021
-
[23]
V. M. G´ alfi, V. Lucarini, F. Ragone, and J. Wouters, Applications of large deviation theory in geophysical fluid dynamics and climate science, Riv. Nuovo Cim. 44, 291 (2021)
work page 2021
-
[24]
V. M. Galfi and V. Lucarini, Fingerprinting Heatwaves and Cold Spells and Assessing Their Response to Climate Change Using Large Deviation Theory, Phys. Rev. Lett. 127, 058701 (2021)
work page 2021
-
[25]
B.-W. Shen, R. A. P. Sr, X. Zeng, J.-J. Baik, S. Faghih- Naini, J. Cui, and R. Atlas, Coexistence of Chaos and Order within a Generalized Lorenz Model,
-
[26]
M. D’Errico, F. Pons, P. Yiou, S. Tao, C. Nardini, F. Lunkeit, and D. Faranda, Present and future synop- tic circulation patterns associated with cold and snowy spells over Italy, Earth Syst. Dynam. 13, 961 (2022)
work page 2022
-
[27]
Y. Ashkenazy and N. R. Smith, Data-driven analysis of annual rain distributions, Phys. Rev. Research 6, 023187 (2024)
work page 2024
-
[28]
P. Grassberger, Information flow and maximum entropy measures for 1-D maps, Physica D: Nonlinear Phenomena 14, 365 (1985)
work page 1985
-
[29]
A. Prasad and R. Ramaswamy, Characteristic distribu- tions of finite-time Lyapunov exponents, Phys. Rev. E 60, 2761 (1999)
work page 1999
-
[30]
Anteneodo, Statistics of finite-time Lyapunov expo- nents in the Ulam map, Phys
C. Anteneodo, Statistics of finite-time Lyapunov expo- nents in the Ulam map, Phys. Rev. E 69, 016207 (2004)
work page 2004
- [31]
-
[32]
Lai, Extreme fluctuations of finite-time lyapunov exponents in chaotic systems, Int
Y.-C. Lai, Extreme fluctuations of finite-time lyapunov exponents in chaotic systems, Int. J. Mod. Phys. B 21, 3941 (2007)
work page 2007
-
[33]
N. R. Smith, Large deviations in chaotic systems: Ex- act results and dynamical phase transition, Phys. Rev. E 106, L042202 (2022)
work page 2022
-
[34]
R. Guti´ errez, A. Canella-Ortiz, and C. P´ erez-Espigares, Finding the Effective Dynamics to Make Rare Events Typical in Chaotic Maps, Phys. Rev. Lett. 131, 227201 (2023)
work page 2023
-
[35]
C. Monthus, Large deviations and conditioning for chaotic non-invertible deterministic maps: Analysis via the forward deterministic dynamics and the backward stochastic dynamics, J. Stat. Mech. 2024, 013208 (2024)
work page 2024
-
[36]
Lippolis, Thermodynamics of chaotic relaxation pro- cesses, Phys
D. Lippolis, Thermodynamics of chaotic relaxation pro- cesses, Phys. Rev. E 110, 024215 (2024)
work page 2024
-
[37]
C. Monthus, Explicit dynamical properties of the Pelikan random map in the chaotic region and at the intermit- tent critical point towards the non-chaotic region, J. Stat. Mech. 2025, 013212 (2025)
work page 2025
- [38]
-
[39]
G. Haller, Distinguished material surfaces and coherent structures in three-dimensional fluid flows, Physica D: Nonlinear Phenomena 149, 248 (2001)
work page 2001
-
[40]
J. Bec, L. Biferale, G. Boffetta, M. Cencini, S. Musacchio, and F. Toschi, Lyapunov exponents of heavy particles in turbulence, Physics of Fluids 18, 091702 (2006)
work page 2006
-
[41]
L. Biferale, C. Meneveau, and R. Verzicco, Deformation statistics of sub-Kolmogorov-scale ellipsoidal neutrally buoyant drops in isotropic turbulence, J. Fluid Mech. 754, 184 (2014)
work page 2014
- [42]
-
[43]
C. De Fazio, J. P. Garrahan, and K. Klobas, Exact results on the dynamics of the stochastic Floquet-East model, J. Phys. A: Math. Theor. 57, 505002 (2024)
work page 2024
-
[44]
J. P. Garrahan, Private communication (2025)
work page 2025
-
[45]
Touchette, The large deviation approach to statistical mechanics, Physics Reports 478, 1 (2009)
H. Touchette, The large deviation approach to statistical mechanics, Physics Reports 478, 1 (2009)
work page 2009
-
[46]
D. Nickelsen and H. Touchette, Anomalous Scaling of Dy- namical Large Deviations, Phys. Rev. Lett. 121, 090602 (2018)
work page 2018
-
[47]
T. Laffargue, K.-D. N. T. Lam, J. Kurchan, and J. Tailleur, Large deviations of Lyapunov exponents, J. Phys. A: Math. Theor. 46, 254002 (2013)
work page 2013
-
[48]
Wikipedia contributors, Tent map, Wikipedia —The Free Encyclopedia
-
[49]
J. C. Leit˜ ao, J. M. V. P. Lopes, and E. G. Altmann, Efficiency of Monte Carlo sampling in chaotic systems, Phys. Rev. E 90, 052916 (2014)
work page 2014
-
[50]
J. C. Leit˜ ao, J. M. V. Parente Lopes, and E. G. Altmann, Importance sampling of rare events in chaotic systems, Eur. Phys. J. B 90, 181 (2017)
work page 2017
- [51]
-
[52]
A. N. Livshits, Homology properties of Y-systems, Math- ematical Notes of the Academy of Sciences of the USSR 10, 758 (1971)
work page 1971
-
[53]
A. N. Livshits, Cohomology of dynamical systems, Math. USSR Izv. 6, 1278 (1972)
work page 1972
-
[54]
Weisstein, Eric W., Complete Elliptic Integral of the First Kind, https://mathworld.wolfram.com/ CompleteEllipticIntegraloftheFirstKind.html (2025)
work page 2025
-
[55]
R. Klages, From deterministic chaos to anomalous diffusion, in Reviews of Nonlinear Dynamics and Complexity, Vol. 3, edited by H. G. Schuster (Wiley – VCH, Weinheim,
-
[56]
H. Fujisaka and S. Grossmann, Chaos-induced diffusion in nonlinear discrete dynamics, Z. Physik B - Condensed Matter 48, 261 (1982)
work page 1982
-
[57]
R. Klages and J. R. Dorfman, Simple Maps with Fractal Diffusion Coefficients, Phys. Rev. Lett. 74, 387 (1995)
work page 1995
-
[58]
R. Klages and J. R. Dorfman, Simple deterministic dy- namical systems with fractal diffusion coefficients, Phys. Rev. E 59, 5361 (1999)
work page 1999
-
[59]
C. P´ erez-Espigares and P. I. Hurtado, Sampling rare events across dynamical phase transitions, Chaos: An Interdisciplinary Journal of Nonlinear Science 29, 10.1063/1.5091669 (2019). 17
- [60]
-
[61]
C. Monthus, Anomalous dynamical large deviations of local empirical densities and activities in the pure and in the random kinetically constrained East model, Eur. Phys. J. B 95, 10.1140/epjb/s10051-022-00281-5 (2022)
-
[62]
C. Monthus, A supersymmetric quantum perspective on the explicit large deviations for reversible Markov jump processes, with applications to pure and random spin chains, J. Stat. Mech. , 073203 (2024)
work page 2024
-
[63]
P. L. Johnson and C. Meneveau, Large-deviation joint statistics of the finite-time Lyapunov spectrum in isotropic turbulence, Physics of Fluids27, 085110 (2015)
work page 2015
- [64]
-
[65]
K. Yoshida, H. Yoshino, A. Shudo, and D. Lippolis, Eigenfunctions of the Perron–Frobenius operator and the finite-time Lyapunov exponents in uniformly hyper- bolic area-preserving maps, J. Phys. A: Math. Theor. 54, 285701 (2021)
work page 2021
-
[66]
P. Grassberger, R. Badii, and A. Politi, Scaling laws for invariant measures on hyperbolic and nonhyperbolic atractors, J Stat Phys 51, 135 (1988)
work page 1988
-
[67]
W. Just, H. Kantz, C. R¨ odenbeck, and M. Helm, Stochas- tic modelling: Replacing fast degrees of freedom by noise, J. Phys. A: Math. Gen. 34, 3199 (2001)
work page 2001
-
[68]
D. Kelly and I. Melbourne, Deterministic homogeniza- tion for fast–slow systems with chaotic noise, Journal of Functional Analysis 272, 4063 (2017)
work page 2017
-
[69]
V. Lucarini and M. Chekroun, Theoretical tools for understanding the climate crisis from Hasselmann’s program and beyond, Nat Rev Phys 5, 744 (2023), arXiv:2303.12009 [physics]
-
[70]
W. R. Inc., Mathematica, Version 13.3, champaign, IL, 2023
work page 2023
-
[71]
C. Sch¨ utte, P. Koltai, and S. Klus, On the numerical ap- proximation of the Perron-Frobenius and Koopman op- erator, JCD 3, 1 (2016)
work page 2016
-
[72]
H. Touchette, Legendre–fenchel transforms in a nutshell, Lecture notes (2014), started July 11, 2005; last compiled October 16, 2014
work page 2014
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.