REVIEW 3 major objections 5 minor 1 cited by
Beyond the Largest Lyapunov Exponent: Entropy-Based Diagnostics of Chaos in Henon-Heiles and N-Body Dynamics
T0 review · 3 major / 5 minor · reviewed 2026-07-13 · grok-4.5
Pith's one-line read Shannon entropy tracks chaos like the Lyapunov exponent, but also sees how N-body mixing changes with particle number.
desk verdict Solid comparative numerics: entropy tracks energy-driven chaos like λ_max in Hénon–Heiles, but the load-bearing N-dependence claim is soft because fixed-bin Shannon entropy can just track smoother potentials. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
Coarse-grained Shannon entropy of a trajectory, estimated from the occupation probabilities of phase-space bins (Freedman-Diaconis reference width, common binning across ensembles, minimum occupancy cut), which measures how broadly an orbit populates accessible phase space rather than the local exponential divergence rate.
What would settle it
Recompute the same ensembles with a substantially different stable binning (or a binning-independent entropy estimator) and check whether the N-dependence of the Shannon entropy still falls while the Lyapunov exponent stays flat; if the entropy trend disappears or reverses under that change, the central complementarity claim fails.
Extended reading notes
Core claim
Trajectory-based coarse-grained Shannon entropy diagnoses chaos in gravitational systems in a way that tracks the largest Lyapunov exponent with energy, yet additionally registers how global phase-space mixing changes with particle number: in Hénon-Heiles the two quantities mirror each other across the onset of widespread chaos, while in live Plummer models both rise for more tightly bound orbits but only the entropy decreases monotonically as N increases.
Load-bearing premise
That the chosen binning and occupancy rules produce an entropy whose trends with energy and particle number reflect real phase-space mixing, not artifacts of how the bins are drawn or of the fact that test-particle energy is not conserved in live N-body runs.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript compares the largest Lyapunov exponent λ_max with trajectory-based Shannon entropy (and a path-length mutual-information entropy I_N) as diagnostics of chaos in the Hénon–Heiles potential and for test particles in live Plummer N-body realizations. In Hénon–Heiles, coarse-grained Shannon entropy tracks the energy-driven transition from weak to widespread chaos and correlates with λ_max. In live N-body runs both diagnostics indicate stronger chaos for more tightly bound orbits, but their N-dependence differs: λ_max is nearly constant over the explored particle numbers while Shannon entropy decreases monotonically with N. The authors conclude that information entropy can complement λ_max, may better capture global phase-space mixing when the leading Lyapunov exponent is uninformative, and is useful when tangent-space dynamics is unavailable.
Significance. If the N-dependence contrast is physical rather than an artifact of fixed coarse-graining, the paper offers a practical, trajectory-only diagnostic that is especially relevant for Solar System minor bodies and other densely sampled orbits where variational equations are costly or unavailable. The dual-system design (smooth mixed phase space plus live finite-N gravity), ensemble averages with scatter, and explicit nearby-binning checks are strengths. The work sits in a mature literature on N-body chaos and Lyapunov scaling; its incremental contribution is the side-by-side comparison of entropy and λ_max under controlled energy and N variation, not a new theorem.
major comments (3)
- The central claim that Shannon entropy 'may better capture changes in global phase-space mixing' rests on the N-body result that λ_max is nearly N-independent while bS_Sh falls monotonically with N (Abstract; Results around Figs. 6–8; Conclusions). The estimator uses a common fixed binning (Freedman–Diaconis reference width, nearby scans, N_samp/N_occ ≥ 5). As N increases the live potential smooths (Fig. 6), so under fixed bins the trajectory occupies fewer cells simply because discreteness noise shrinks. That can produce an entropy drop without any change in mixing that λ_max misses. The manuscript needs a control that isolates this: e.g. (i) recompute bS_Sh in the continuum/smooth Plummer limit with the same binning, (ii) scale bin widths with the local force fluctuation or with N, or (iii) report occupied-bin counts and phase-space volume explored versus N. Without such a control the
- In live N-body runs test-particle energy is not conserved (explicitly noted in the N-body Results). The paper therefore plots λ_max against median energy (Fig. 8) and reports that both diagnostics strengthen for more tightly bound orbits. Because the energy distribution itself can depend on N (and on how long the particle is integrated), the claimed pure N-dependence of bS_Sh is not cleanly separated from energy drift. A clearer isolation is needed: either restrict the comparison to trajectories whose median (or time-averaged) energy lies in a narrow common window across N, or show bS_Sh versus N at fixed median energy bins with the same occupancy cuts. Otherwise the monotonic decline of entropy with N may partly track a shift in the energy sampling rather than a distinct mixing diagnostic.
- The path-length mutual-information entropy I_N is carefully defined (Eqs. 19–26) and shown versus λ_max for Hénon–Heiles (Fig. 3), yet the N-scaling that drives the abstract and conclusions is reported only for the coarse-grained Shannon entropy bS_Sh. Given that I_N is constructed to be sensitive to dynamical divergence of nearby trajectories (and is therefore closer in spirit to a Lyapunov diagnostic), the manuscript should either (a) show I_N versus N for the same Plummer test-particle ensembles, or (b) state explicitly why I_N is omitted from the N-body comparison and why bS_Sh alone is the appropriate complement. Leaving I_N unused for the key claim weakens the argument that entropy-based measures as a class outperform λ_max for global mixing.
minor comments (5)
- Figure numbering in the supplied text jumps (Fig. 3 then Fig. 6); ensure all intermediate figures are present and consistently referenced in the final PDF.
- Notation for the coarse-grained Shannon estimator appears as bS_Sh / \hat{S}_Sh; pick one symbol and use it uniformly in text, equations, and figure captions.
- The abstract and conclusions emphasize suitability for Solar System minor bodies; a short quantitative remark on typical sampling cadence or phase-space dimension for which the occupancy cut remains feasible would make that claim more concrete.
- Clarify early whether the N-body experiments are pure test particles in live Plummer realizations (as the text indicates) or include self-gravity of the tracer; the distinction matters for energy non-conservation.
- A few references in the introduction (e.g. on few-body chaos and galactic potentials) are dense; a brief sentence separating Hamiltonian chaos diagnostics from collisional relaxation timescales would help non-specialist readers.
Circularity Check
No significant circularity: λ_max and Shannon entropy are independently computed from the same trajectories; reported correlations are post-hoc, not definitional or fitted predictions.
full rationale
The paper’s load-bearing comparison is between two separately defined numerical diagnostics. The largest Lyapunov exponent is obtained from tangent-space evolution (Eq. 1–3 and the standard renormalization algorithm). The coarse-grained Shannon entropy bSSh is a histogram estimator on stored phase-space samples (Eq. 18), with binning fixed by Freedman–Diaconis plus stability scans and an occupancy cut; the mutual-information path-length entropy IN is defined from arc-length densities of a reference and a nearby trajectory (Eqs. 19–26). Neither entropy is defined in terms of λ_max, nor is λ_max fitted from entropy. Figure 3’s dashed lines are explicitly “best-fitting relation[s] between the entropy and λ_max” after both quantities have been measured—descriptive correlations, not inputs that force a “prediction.” The N-body claim that λ_max is nearly N-independent while Shannon entropy falls with N is likewise an empirical contrast from independent measurements on the same live Plummer realizations, not a quantity recovered by construction from a fitted parameter. Self-citations (e.g. Di Cintio & Casetti 2019; Sartorello et al. 2025) supply background that λ_max depends weakly on N; they do not define or force the new entropy–N trend. Methodological concerns about whether the entropy drop with N is partly a binning/graininess artifact are correctness risks, not circular reductions. The derivation chain is therefore self-contained against its own inputs.
Assumptions & free parameters
free parameters (4)
- phase-space bin widths (Freedman-Diaconis reference + nearby scan)
- minimum average occupancy Nsamp/Nocc ≥ 5
- initial separation d0 for nearby trajectories (I_N and Lyapunov)
- explored N range and energy window for Plummer test particles
assumptions (4)
- domain assumption Largest Lyapunov exponent defined by the standard tangent-space limit (Eq. 1) correctly quantifies local exponential divergence for the systems studied.
- domain assumption Coarse-grained Shannon entropy of occupied phase-space bins (Eq. 18) is a meaningful measure of orbital complexity / phase-space exploration.
- domain assumption Path-length mutual-information entropy I_N of Núñez et al. (1996) is a valid finite-time instability diagnostic distinct from λ_max.
- domain assumption Live N-body Plummer realizations with test particles adequately probe finite-N chaos despite non-conserved test-particle energy.
Cite this review
Pith. "Pith review of Beyond the Largest Lyapunov Exponent: Entropy-Based Diagnostics of Chaos in Henon-Heiles and N-Body Dynamics." pith.science (2026). https://pith.science/paper/PQ2ODOV3
@misc{pith2026260324675,
author = {Pith},
title = {Pith review of: Beyond the Largest Lyapunov Exponent: Entropy-Based Diagnostics of Chaos in Henon-Heiles and N-Body Dynamics},
year = {2026},
howpublished = {\url{https://pith.science/paper/PQ2ODOV3}},
note = {Machine review of arXiv:2603.24675}
}
read the original abstract
The largest Lyapunov exponent is widely used to diagnose chaos in gravitational dynamics, but in mixed phase spaces and finite-N systems it does not always provide a complete description of orbital complexity and phase-space transport. Entropy-based diagnostics may offer a complementary perspective. We investigate whether trajectory-based information entropy can provide a useful diagnostic of chaos in gravitational systems and how it relates to the largest Lyapunov exponent as a function of orbital energy and of the number of degrees of freedom. We computed the largest Lyapunov exponent and a coarse-grained Shannon entropy for ensembles of trajectories in the Henon-Heiles potential and for test-particle orbits in live N-body realizations of a Plummer model. We then compared the dependence of both quantities on orbital energy and, for the N-body case, on particle number. In the Henon-Heiles system, the Shannon entropy follows the transition from weak to widespread chaos and exhibits an energy dependence that closely mirrors that of the largest Lyapunov exponent. For test-particle orbits in live N-body potentials, both diagnostics indicate stronger chaos for more tightly bound trajectories. However, their dependence on N differs: the largest Lyapunov exponent remains nearly constant over the explored range of particle numbers, whereas the Shannon entropy decreases monotonically as N increases. These results show that the information entropy can complement the largest Lyapunov exponent and may better capture changes in global phase-space mixing, especially in systems where the leading Lyapunov exponent alone is not sufficiently informative. It therefore provides a promising alternative for diagnosing chaos when tangent-space dynamics is unavailable or computationally expensive, and it is naturally suited to systems with densely sampled trajectories, such as minor bodies in the Solar System.
Forward citations
Cited by 1 Pith paper
-
On an empirical method to build near-equilibrium axisymmetric and triaxial galaxy models
An improved adiabatic-squeezing method with an Ornstein-Uhlenbeck noise term lets modelers prescribe the final axial ratios of near-equilibrium triaxial N-body galaxy models.
Reviewed July 13, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.