{"id":"f51ddb81-8f3a-49c2-af36-4235f1da5c1f","arxiv_id":"2603.24675","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":4,"one_line_summary":"Trajectory Shannon entropy mirrors Lyapunov chaos trends with energy but decreases with N while λ_max stays nearly constant, better capturing phase-space mixing in finite-N gravity.","lead":"Shannon entropy of orbital trajectories tracks chaos like the largest Lyapunov exponent in Hénon-Heiles and N-body systems, but falls with particle number while Lyapunov stays flat. This gives a practical chaos diagnostic when tangent-space methods are unavailable, as for Solar System minor bodies.","discovery_kind":"new_application","skeptic_critique":{"model":"grok-4.5","headline":"The N-dependence claim rests on coarse-grained Shannon entropy whose decline with N may largely track shrinking graininess rather than a distinct diagnostic of phase-space mixing.","rationale":"The Reader correctly flags coarse-graining and live-potential energy non-conservation as the weakest assumption. That assumption is load-bearing for the paper's distinctive claim: not merely that entropy correlates with λ_max (which Hénon–Heiles supports), but that entropy better tracks global mixing via an N-dependence that λ_max lacks. The supplied text shows λ_max nearly constant in N (Fig. 7 and surrounding discussion) and asserts a monotonic entropy decrease, yet does not demonstrate that the decrease survives N-adaptive binning or frozen-potential controls. If those checks preserve the trend, the claim strengthens; if not, the complementarity is overstated. No internal inconsistency or critical error is evident, so the verdict remains CONDITIONAL rather than REJECT, pending exactly those robustness tests (and public code). Confidence stays moderate because the cacheable source includes methods and key figures but not exhaustive binning tables.","tokens_in":6392,"tokens_out":660,"duration_ms":7476,"concrete_test":"Recompute bSSh for the same N-body ensembles with (i) bin widths scaled to the measured phase-space volume of each trajectory (or to N-dependent graininess), and (ii) a frozen Plummer potential at each N so energy is conserved. If the monotonic decline of entropy with N disappears or reverses under either change while λ_max stays flat, the claimed complementary N-dependence is not robust.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's strongest claim is that trajectory Shannon entropy complements λ_max and may better capture global phase-space mixing because, in live Plummer N-body runs, λ_max is nearly N-independent while bSSh falls monotonically with N (Abstract; Results on N-body; Conclusions). That contrast is load-bearing: without a physically meaningful N-dependence of entropy, the claim reduces to the weaker (and already known) statement that entropy tracks energy-driven chaos similarly to λ_max in Hénon–Heiles. The estimator is a fixed common binning of phase-space samples (Freedman–Diaconis reference width, nearby scans, Nsamp/Nocc ≥ 5; Methods on bSSh). As N rises, the live potential becomes smoother (Fig. 6), so trajectories explore a less grainy effective phase space and occupy fewer bins under that fixed coarse-graining; the reported entropy drop can therefore be an artifact of bin occupancy under decreasing discreteness noise rather than a measure of mixing that λ_max misses. Energy non-conservation of test particles further couples median energy to N, so the N-trend is not cleanly isolated. The mutual-information path-length entropy IN is defined but not used for the N-scaling that drives the claim.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","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.","tokens_in":6731,"tokens_out":1330,"duration_ms":11892,"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":[{"comment":"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","section":null},{"comment":"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.","section":null},{"comment":"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.","section":null}],"minor_comments":[{"comment":"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.","section":null},{"comment":"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.","section":null},{"comment":"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.","section":null},{"comment":"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.","section":null},{"comment":"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.","section":null}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid numerical methods/diagnostics contribution appropriate for A&A if the N-dependence control is added. The main risk is overclaiming that entropy 'better captures mixing' when the reported N-trend may be largely a fixed-bin discreteness effect; that is fixable within the existing experimental setup. Scope fit is good for celestial mechanics / computational dynamics."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The one new result worth knowing is the side-by-side comparison in live Plummer N-body: for the same test-particle setups, λ_max stays roughly flat with N while coarse-grained Shannon entropy falls monotonically, while both quantities track energy the same way (tighter binding → more chaos). In Hénon–Heiles the two diagnostics simply move together across the weak-to-strong chaos transition. That contrast is the paper’s actual contribution; the rest is careful application of known tools.\n\nWhat they do well is straightforward. Ensembles with 1σ scatter, common binning within comparison sets, Freedman–Diaconis reference widths plus nearby scans, and an occupancy cut are all reported. The Hénon–Heiles energy trends look clean and the λ_max–energy relation in live N-body matches earlier noisy-potential and low-N work they cite. The methods are independent: neither diagnostic is fitted to the other. Citations are appropriate and not circular.\n\nThe soft spot is real but not fatal. The N-scaling that drives the “entropy better captures global mixing” claim uses fixed common phase-space bins. As N rises the live potential smooths (their Fig. 6), so trajectories under that fixed coarse-graining occupy fewer bins; the entropy drop can therefore be discreteness graininess rather than a mixing measure that λ_max misses. Energy is not conserved, so they use median energy; that couples N and energy somewhat. Mutual-information path-length entropy is defined but not used for the N-trend. These are acknowledged limitations, not hidden ones, and they do not erase the comparative result.\n\nThis is for people who already compute chaos diagnostics in celestial mechanics or stellar dynamics and want a practical alternative when tangent-space work is expensive or unavailable (minor-body time series, large surveys). It does not reorganize the field. I would send it to peer review; a referee should demand clearer isolation of the N-effect (adaptive binning, frozen vs live controls, or continuum-limit comparison) and public code/data. Worth reading if you work on these diagnostics; not urgent otherwise.","headline":"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.","tokens_in":7343,"tokens_out":590,"would_cite":false,"duration_ms":5472,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Shannon entropy tracks chaos like the Lyapunov exponent, but also sees how N-body mixing changes with particle number.","keywords":["chaos","Lyapunov exponent","Shannon entropy","Hénon-Heiles","N-body dynamics","phase-space mixing","Plummer model","celestial mechanics"],"falsifier":"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.","tokens_in":7281,"feed_emoji":"✨","tokens_out":637,"duration_ms":5659,"temperature":0.7,"pith_summary":"The paper asks whether a coarse-grained Shannon entropy computed from orbital trajectories can diagnose chaos in gravitational systems as usefully as the largest Lyapunov exponent, and whether it captures more of the story when phase space is mixed or the system has finitely many particles. In the Hénon-Heiles potential, the entropy follows the same energy-driven transition from weak to widespread chaos that the Lyapunov exponent does. For test particles in live Plummer N-body models, both measures say tightly bound orbits are more chaotic, but only the entropy falls as particle number rises while the Lyapunov exponent stays roughly flat. The authors conclude that information entropy can complement or replace the Lyapunov exponent when tangent-space dynamics are unavailable or expensive, and that it is a natural tool for densely sampled trajectories such as minor bodies in the Solar System.","feed_headline":"Entropy tracks chaos where Lyapunov alone misses N-body mixing","feed_subtitle":"In live Plummer models the entropy falls with particle number while the largest Lyapunov exponent stays flat.","key_machinery":"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.","core_discovery":"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.","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["Shannon entropy tracks Lyapunov but captures N-body mixing changes","Entropy mirrors Lyapunov on energy yet falls with N as Lyapunov stays flat","Trajectory entropy diagnoses chaos where Lyapunov misses global mixing","Entropy reveals N-dependent phase-space mixing missed by largest Lyapunov","Coarse-grained Shannon entropy tracks chaos onset and drops as N grows"],"cache_read_input_tokens":128,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Shannon entropy tracks Lyapunov but captures N-body mixing changes","Entropy mirrors Lyapunov on energy yet falls with N as Lyapunov stays flat","Trajectory entropy diagnoses chaos where Lyapunov misses global mixing","Entropy reveals N-dependent phase-space mixing missed by largest Lyapunov","Coarse-grained Shannon entropy tracks chaos onset and drops as N grows"]},"model":"grok-4.5","effort":"low","cost_usd":0.005766,"raw_usage":{"total_tokens":1605,"prompt_tokens":870,"num_sources_used":0,"completion_tokens":69,"cost_in_usd_ticks":57660000,"prompt_tokens_details":{"text_tokens":870,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":666,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":870,"tokens_out":69,"duration_ms":5464,"temperature":1.0,"reasoning_tokens":666,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-13T18:41:50.110455+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"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.","supporting_citations":[],"review_version":1}