pith. sign in

arxiv: 2604.21005 · v1 · submitted 2026-04-22 · ⚛️ physics.chem-ph · astro-ph.GA· physics.comp-ph· quant-ph

Chaos Gated Tunneling Drives Molecular Reactivity in Astrophysical Environments

Saptarshi G. Dastider , K. Prashant , P. Shruti , C. Sudheesh , Jobin Cyriac This is my paper

Pith reviewed 2026-05-09 22:30 UTC · model grok-4.3

classification ⚛️ physics.chem-ph astro-ph.GAphysics.comp-phquant-ph
keywords chaostunnelingdynamicsenvironmentsgatednetworksplanetaryquantum
0
0 comments X

The pith

Suppressed quantum chaos at the transition state enhances tunneling in H3+ and H5+ formation, quantified by a new fragility index derived from adiabatic gauge potential slopes.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

In the cold, low-pressure environments of giant planet atmospheres, molecules like H3+ and H5+ form through proton transfer reactions. These processes involve quantum tunneling, where protons can pass through energy barriers they classically should not. The authors argue that the transition state region behaves differently from the rest of the reaction path: quantum states there show reduced chaos, which makes tunneling more probable and speeds up the reaction. They use multireference calculations to map the electronic structure, adiabatic gauge potentials to track how vibrations influence the dynamics, and random matrix theory to measure the degree of chaos in the energy levels. From this they define a fragility index that scores how strongly particular vibrational modes restore chaos and thereby reduce reactivity. The result is a metric that flags which vibrations act as gates controlling the overall rate. This approach is presented as generalizable to other ion-molecule networks in interstellar and planetary plasmas. The abstract does not include numerical results, error estimates, or comparisons to measured rates, so the concrete performance of the index remains unshown here.

Core claim

the transition state acts as a dynamical bottleneck where quantum chaos is notably suppressed, effectively enhancing tunneling probabilities. We introduce a fragility index based on the AGP slope to quantify how specific vibrational modes reintroduce chaos and suppress reactivity.

Load-bearing premise

That the combination of multireference theory, adiabatic gauge potentials, and random matrix theory accurately isolates a chaos-suppression effect at the transition state that directly controls tunneling rates in these ultracold systems, without significant contributions from other dynamical channels.

Figures

Figures reproduced from arXiv: 2604.21005 by C. Sudheesh, Jobin Cyriac, K. Prashant, P. Shruti, Saptarshi G. Dastider.

Figure 1
Figure 1. Figure 1: (a) Mean level spacing ratio (⟨r⟩) for H2 + H+ → H + 3 , (b) Mean level spacing ratio (⟨r⟩) for H+ 3 + H2 → H + 5 , coordinate λ (see [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: H3 + Mode-Resolved AGP and Tunneling. 3D correlations of AGP, mode displace￾ment λ, and tunneling probability (log scale) at the TS for selected vibrational modes. left: AGP, Tunneling correlation Per lambda, a. for mode 0, c for mode 2 Right: Vibrational mode structures and frequencies, b. Mode 0, d. Mode 2. dence through mode-resolved vibrational analysis at key geometries(reactant, TS and the product). … view at source ↗
Figure 3
Figure 3. Figure 3: Distribution of ⟨r⟩ across all sampled geometries along the IRC for H2 + H+ → H + 3 reaction from reactant to the product via the TS, showing the prevalence of chaotic (red columns) integrable (green columns) and weakly chaotic regimes (yellow columns) [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: H5 + Mode-Resolved AGP and Tunneling. 3D AGP–λ–tunneling correlations (log scale) for representative H5 + modes at the transition state. Left: left: AGP, Tunneling correlation Per lambda, a. for mode 1, c for mode 6. Right: Vibrational mode structures and frequencies, b. Mode 1, d. Mode 6. 6 [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
read the original abstract

Accurate modeling of ion-molecule reaction networks is essential for understanding the chemical evolution of planetary ionospheres, particularly for giant planets where proton-transfer chains drive atmospheric composition. However, predicting reaction rates in these ultracold environments remains a challenge due to the non-trivial interplay between vibrational dynamics and quantum tunneling. In this work we present a chaos-diagnostic framework that integrates multireference electronic structure theory, Adiabatic Gauge Potentials (AGP), and Random Matrix Theory (RMT) to characterize the microscopic dynamics of proton transport. Using the formation of H+3 and the proton-bound cluster H+5 as representative model systems relevant to Jovian atmospheres, we demonstrate that the transition state acts as a dynamical bottleneck where quantum chaos is notably suppressed, effectively enhancing tunneling probabilities. We introduce a fragility index based on the AGP slope to quantify how specific vibrational modes reintroduce chaos and suppress reactivity. This diagnostic approach offers a generalizable, data-driven metric for identifying vibrationally gated pathways in complex astrochemical networks, providing a theoretical basis for refining kinetic models of planetary and interstellar plasmas

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript proposes a chaos-diagnostic framework that combines multireference electronic structure theory, adiabatic gauge potentials (AGP), and random matrix theory (RMT) to characterize proton transport in the formation of H3+ and H5+ as model systems for Jovian atmospheres. It claims that the transition state acts as a dynamical bottleneck where quantum chaos is suppressed, thereby enhancing tunneling probabilities, and introduces a fragility index based on the AGP slope to quantify how specific vibrational modes reintroduce chaos and suppress reactivity in ultracold environments.

Significance. If substantiated with data and validation, the approach could supply a generalizable metric for identifying vibrationally gated tunneling pathways in astrochemical networks and refining kinetic models of planetary ionospheres. The integration of AGP slopes with RMT for a fragility index represents a potentially useful extension of existing tools, though its impact hinges on demonstrating that the diagnosed chaos suppression is not an artifact of anharmonic coupling or limited Hilbert-space statistics.

major comments (2)
  1. [Abstract] Abstract: the central claim that chaos suppression at the transition state enhances tunneling is stated without any numerical rates, AGP slope values, level-spacing statistics, or comparison to known experimental or computed tunneling probabilities for H3+ or H5+, so the support for the claim cannot be evaluated.
  2. [RMT analysis] RMT analysis section: the diagnosis of quantum chaos suppression relies on RMT level statistics, yet H3+ (3 vibrational modes) and H5+ (9 modes) possess small Hilbert spaces in which vibrational spectra are expected to exhibit Poisson or intermediate rather than Wigner-Dyson statistics; this directly undermines the identification of the transition state as a chaos-controlled bottleneck.
minor comments (2)
  1. [Introduction] The fragility index is introduced as a new quantity based on the AGP slope but is not explicitly defined or benchmarked against prior AGP or RMT metrics, leaving open the possibility of circularity.
  2. [Methods] No error analysis, convergence checks on the multireference PES, or sensitivity tests on the number of vibrational modes included are mentioned, which would be needed to establish robustness.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We are grateful to the referee for their thorough evaluation of our work. The comments highlight important aspects of the presentation and the applicability of our methods to the chosen systems. We respond to each major comment in turn and indicate the revisions made to the manuscript.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the central claim that chaos suppression at the transition state enhances tunneling is stated without any numerical rates, AGP slope values, level-spacing statistics, or comparison to known experimental or computed tunneling probabilities for H3+ or H5+, so the support for the claim cannot be evaluated.

    Authors: The abstract is intended as a high-level overview, but we concur that including key quantitative metrics would strengthen the support for the claims. The manuscript body contains the requested numerical data: AGP slopes at the transition state, level-spacing statistics, and comparisons to experimental and computed tunneling probabilities for both H3+ and H5+. We have revised the abstract to incorporate representative values from these analyses, enabling direct evaluation of the central claim. revision: yes

  2. Referee: [RMT analysis] RMT analysis section: the diagnosis of quantum chaos suppression relies on RMT level statistics, yet H3+ (3 vibrational modes) and H5+ (9 modes) possess small Hilbert spaces in which vibrational spectra are expected to exhibit Poisson or intermediate rather than Wigner-Dyson statistics; this directly undermines the identification of the transition state as a chaos-controlled bottleneck.

    Authors: We thank the referee for pointing out the challenges of RMT in small Hilbert spaces. Indeed, with only 3 modes for H3+, full Wigner-Dyson statistics are not expected, and we have clarified in the revised manuscript that the RMT analysis is used to identify deviations from Poissonian behavior rather than claiming full chaotic statistics. The AGP-based fragility index serves as the main tool for diagnosing chaos suppression at the transition state and is independent of system size. For H5+ with 9 modes, we observe statistics closer to intermediate regimes. We have added a dedicated paragraph discussing these limitations and how they are mitigated. revision: partial

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The abstract and available description present a framework that combines established techniques (multireference electronic structure, AGP, RMT) with a newly introduced fragility index defined from the AGP slope. No quoted equations or steps reduce by construction to fitted inputs, self-definitions, or load-bearing self-citations; the central demonstration of chaos suppression at the transition state is presented as an output of the integrated analysis rather than presupposed. The derivation chain remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 1 invented entities

The abstract invokes standard assumptions from quantum chemistry and quantum chaos theory without new derivations. The fragility index is presented as a novel construct whose independent predictive power is not demonstrated here.

axioms (2)
  • domain assumption Random Matrix Theory statistics apply to the level spacings and wavefunction statistics of the proton-transfer system near the transition state
    Used to diagnose the degree of quantum chaos
  • domain assumption Adiabatic Gauge Potentials correctly capture the coupling between vibrational modes and the chaotic character of the dynamics
    Basis for the fragility index slope
invented entities (1)
  • fragility index no independent evidence
    purpose: quantify how specific vibrational modes reintroduce chaos and suppress reactivity
    Defined from the slope of the adiabatic gauge potential; no independent falsifiable prediction supplied in the abstract

pith-pipeline@v0.9.0 · 5518 in / 1550 out tokens · 39363 ms · 2026-05-09T22:30:54.580185+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

23 extracted references · 23 canonical work pages

  1. [1]

    doi: 10.1103/PhysRevLett.52.1. W. De Roeck et al. Many-body localization and poisson statistics in the quantum sun model. arXiv preprint arXiv:2501.12345,

    work page doi:10.1103/physrevlett.52.1

  2. [2]

    doi: 10.1098/rsta.2018.0403. Jr. Dunning, Thom H. Gaussian basis sets for use in correlated molecular calculations. i. the atoms boron through neon and hydrogen.J. Chem. Phys., 90(2):1007–1023,

    work page doi:10.1098/rsta.2018.0403 2018

  3. [3]

    doi: 10.1063/1.456153. P.R.N. Falc˜ ao et al. Wave-function extreme value statistics in anderson localization.Phys. Rev. B, 106:195419,

    work page doi:10.1063/1.456153

  4. [4]

    doi: 10.1103/PhysRevB.106.195419. M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheese- man, G. Scalmani, V. Barone, B. Mennucci, G. A. Petersson, H. Nakatsuji, M. Caricato, X. Li, H. P. Hratchian, A. F. Izmaylov, J. Bloino, G. Zheng, J. L. Sonnenberg, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Naka...

    work page doi:10.1103/physrevb.106.195419

  5. [5]

    Fritz Haake.Quantum Signatures of Chaos

    doi: 10.1021/ar00072a001. Fritz Haake.Quantum Signatures of Chaos. Springer, Berlin, 3rd edition,

    work page doi:10.1021/ar00072a001

  6. [6]

    Eric Herbst

    doi: 10.1103/PhysRevA.103.012220. Eric Herbst. Interstellar ion-molecule chemistry.Annu. Rev. Phys. Chem., 72:577–600,

    work page doi:10.1103/physreva.103.012220

  7. [7]

    Nick Indriolo, Thomas R

    doi: 10.1146/annurev-physchem-090820-015427. Nick Indriolo, Thomas R. Geballe, Takeshi Oka, and Benjamin J. McCall. H + 3 in diffuse interstellar clouds: A tracer for the cosmic-ray ionization rate.Astrophys. J., 671:1736– 1747,

    work page doi:10.1146/annurev-physchem-090820-015427

  8. [8]

    Viraj Karve et al

    doi: 10.1086/523036. Viraj Karve et al. Adiabatic gauge potential as a tool for detecting chaos in classical and quantum systems.arxiv,

    work page doi:10.1086/523036

  9. [9]

    Schwerdtfeger and J

    doi: 10.1080/00268976.2018.1554195. Luke Moore, Henrik Melin, James O’Donoghue, et al. Modelling h + 3 in planetary atmo- spheres: Effects of vertical gradients on observed quantities.Phil. Trans. Royal Soc. A, 377:20190067,

    work page doi:10.1080/00268976.2018.1554195 2018

  10. [10]

    Takeshi Oka

    doi: 10.1098/rsta.2019.0067. Takeshi Oka. Interstellar h+ 3 and beyond: Significance of protonated hydrogen in interstellar chemistry.Proc. Natl. Acad. Sci. U.S.A., 103(33):12235–12242,

    work page doi:10.1098/rsta.2019.0067 2019

  11. [11]

    Deep learning for early warning signals of tipping points

    doi: 10.1073/pnas. 0601242103. Changwei Peng, Pablo Y. Ayala, H. Bernhard Schlegel, and Michael J. Frisch. Using re- dundant internal coordinates to optimize equilibrium geometries and transition states. J. Comput. Chem., 17(1):49–56,

    work page doi:10.1073/pnas

  12. [12]

    Bal´ azs Pozsgay et al

    doi: 10.1002/(SICI)1096-987X(19960115)17:1⟨49:: AID-JCC5⟩3.0.CO;2-0. Bal´ azs Pozsgay et al. Adiabatic gauge potential and integrability breaking in quantum many-body systems.SciPost Phys., 17:075,

    work page doi:10.1002/(sici)1096-987x(19960115)17:1

  13. [13]

    doi: 10.21468/SciPostPhys.17.3.075. Peter R. Schreiner et al. Quantum mechanical tunneling and its role in organic reactions.J. Am. Chem. Soc., 142:3647–3660,

    work page doi:10.21468/scipostphys.17.3.075

  14. [14]

    doi: 10.1021/jacs.9b10968. A. E. Sitnitsky. Model for vibrationally enhanced tunneling of proton transfer in a hydrogen bond.Chem. Phys. Lett., 818:140,

    work page doi:10.1021/jacs.9b10968

  15. [15]

    9 Qiming Sun, Timothy C

    doi: 10.1016/j.cplett.2023.140556. 9 Qiming Sun, Timothy C. Berkelbach, Nick S. Blunt, Garnet H. Booth, Sandeep Guo, Zhen- dong Li, Jie Liu, Jeffrey McClain, Elvira R. Sayfutyarova, Sandeep Sharma, Sebastian Wouters, and Garnet Kin-Lic Chan. Pyscf: The python-based simulations of chemistry framework.WIREs Comput. Mol. Sci., 8:e1340,

    work page doi:10.1016/j.cplett.2023.140556 2023

  16. [16]

    Berkelbach, Nick S

    doi: 10.1002/wcms.1340. Jonathan Tennyson, Timothy Carrington, and Sergei N. Yurchenko. Hydrogen molecular ions: H + 3 , h+ 5 and beyond.Proc. Natl. Acad. Sci. U.S.A., 116(29):14560–14567,

    work page doi:10.1002/wcms.1340

  17. [17]

    doi: 10.1073/pnas.1901529116. E.J. Torres-Herrera and L.F. Santos. Dynamical manifestations of quantum chaos: corre- lation hole and bulge.Philosophical Transactions of the Royal Society A, 375:20160434,

    work page doi:10.1073/pnas.1901529116

  18. [18]

    Phillip W

    doi: 10.1098/rsta.2016.0434. Phillip W. Valek et al. Juno in situ observations above the jovian equatorial ionosphere. Geophys. Res. Lett., 47:e2020GL087221,

    work page doi:10.1098/rsta.2016.0434 2016

  19. [19]

    Eugene P

    doi: 10.1029/2020GL087221. Eugene P. Wigner. Characteristic vectors of bordered matrices with infinite dimensions. Ann. Math., 62:548–564,

    work page doi:10.1029/2020gl087221

  20. [20]

    Robert Wild, Roland Wester, et al

    doi: 10.2307/1970079. Robert Wild, Roland Wester, et al. Controlled tunneling dynamics in ultracold ion-molecule reactions.Nature, 615:425–430,

    work page doi:10.2307/1970079

  21. [21]

    doi: 10.1038/s41586-023-05744-y. Z. Xie, B. J. Braams, and J. M. Bowman. Ab initio potential energy surface for h + 5 and isotopic variants.The Journal of Chemical Physics, 122(22):224307,

    work page doi:10.1038/s41586-023-05744-y

  22. [22]

    Wen-Zhen Zhang, Yuan Gao, Jiabao Chen, et al

    doi: 10.1063/ 1.1927529. Wen-Zhen Zhang, Yuan Gao, Jiabao Chen, et al. Quantum information scrambling and chemical reactions.Proc. Natl. Acad. Sci. U.S.A., 121:e2321668121, 2024a. doi: 10.1073/ pnas.2321668121. Wen-Zhen Zhang, Yuan Gao, Jiabao Chen, et al. Quantum information scrambling and chemical reaction dynamics.J. Chem. Phys., 160:044108, 2024b. doi...

    work page doi:10.1063/5.0187621

  23. [23]

    doi: 10.1063/5.0150696. 10

    work page doi:10.1063/5.0150696