Pith. sign in

REVIEW 2 major objections 4 minor 81 references

High-Energy Reaction Dynamics of O$_3$

T0 review · 2 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash

Pith's one-line read The reef near dissociation does not control the temperature dependence of the O+O2 exchange rate.

desk verdict New O+O2 PES and QCT rates are solid; the reef conclusion is plausible but underdetermined by an uncontrolled PES-to-PES comparison. read the letter →

arxiv 2506.06088 v1 pith:GNMCQOLE submitted 2025-06-06 physics.chem-ph

classification physics.chem-ph
keywords O3ozoneO+O2exchangereactionpotentialenergysurfacereefstructurequasi-classicaltrajectoriestemperature-dependentratecoefficientsneuralnetworkstate-to-distributionmodel
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper sets out to test a long-standing attribution: that the 'reef'—a submerged barrier in the O+O2 entrance channel near the dissociation threshold—is what makes the computed O+O2 atom-exchange rate rise with temperature, against experiment. The authors construct a new global potential energy surface from high-level multireference quantum-chemistry energies and run millions of quasi-classical trajectories on it and on an earlier surface of comparable quality. Both surfaces contain reefs, yet both reproduce the measured negative temperature dependence of the exchange rate within error bars. The paper concludes that the reef is not the feature controlling that temperature dependence. The same simulations reproduce the temperature trend for dissociation but underestimate absolute dissociation rates by about two orders of magnitude, and they support a neural-network model that predicts product state distributions for use in hypersonic flow simulations.

What carries the argument

The load-bearing object is the 'reef': a small, angle-dependent submerged barrier in the minimum-energy path for O approaching O2, located within about 10 kcal/mol of the dissociation limit and visible in one-dimensional cuts of the new surface. The argument is carried by comparing quasi-classical trajectory thermal rates on two reef-bearing surfaces—the new reproducing-kernel surface and an earlier permutationally invariant polynomial surface built at a comparable electronic-structure level—with earlier reef-free and reef-bearing surfaces from the literature. Because the two reef-bearing surfaces independently reproduce the measured negative temperature dependence, the reef loses its status as the controlling feature. The new surface itself is built with a reproducing kernel Hilbert space interpolation of about 6300 energies calculated at the MRCI+Q/aug-cc-pVTZ level, with the three O+O2 channels mixed through exponential switching functions.

What would settle it

Flatten only the reef on the new surface, leaving every other part of the potential unchanged, and rerun the quasi-classical exchange rates: if the negative temperature dependence flips sign, the reef is controlling after all. A complementary check is to recompute the near-dissociation region at aug-cc-pVQZ or complete-basis-set level to see whether the reef survives with the same height and angle dependence.

Watch

Extended reading notes

Core claim

The paper's central claim is that the 'reef'—an angle-dependent, submerged barrier in the O2–O entrance channel within about 10 kcal/mol of dissociation—does not control the sign of the temperature dependence of the O+O2 atom-exchange rate. Quasi-classical trajectories on two independent global surfaces, one built here by reproducing-kernel interpolation of high-level multireference configuration-interaction energies and one an earlier permutationally invariant polynomial surface, both contain reefs, and both give the measured negative temperature dependence of the exchange rate within experimental error. Earlier surfaces with a reef had been associated with a positive temperature dependence, and a reef-free surface with the correct negative one; the present comparison breaks that link. The paper also reports that the same trajectories reproduce the experimental temperature trend for O2 dissociation while underestimating absolute dissociation rates by about two orders of magnitude, and that a neural-network state-to-distribution model trained on a subset of initial conditions predicts final translational, vibrational, and rotational distributions well for both surfaces.

Load-bearing premise

The conclusion depends on the near-dissociation region of the new surface being accurate enough, but the surface binds O3 about 4 kcal/mol more weakly than the accepted experimental value, so the reef's height and shape may be distorted by the level of electronic-structure theory.

Editorial extensions

If this is right

  • If the reef is not rate-controlling, improving agreement with measured exchange rates should focus on other features of the O+O2 surface, such as the long-range anisotropy and the barrier region closer to the minimum.
  • The new surface completes a family of surfaces for the NOO, NNO, NNN, and OOO reactive systems built at the same level of quantum-chemical theory, so the four can be used together in reaction-network simulations of air.
  • The neural-network state-to-distribution models reproduce full product translational, vibrational, and rotational distributions from initial conditions not used in training, which means thermally relevant final-state information can be obtained without running new trajectories for every initial state.
  • Dissociation rates computed on both surfaces keep the measured temperature trend but are low by about two orders of magnitude; using a larger electronic degeneracy factor consistent with high-temperature conditions brings them to within roughly one order of magnitude.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The authors do not say this, but the result implies that earlier 'reef removed versus reef present' comparisons did not isolate the reef: those surfaces also differed in other regions, so the attributed change in rate sign may have come from those other differences.
  • Because the new surface binds O3 about 4 kcal/mol more weakly than the accepted experimental dissociation energy and the reef sits within 10 kcal/mol of dissociation, whether a given surface has a reef may partly be a basis-set artifact; a systematic aug-cc-pVQZ or complete-basis-set series would show whether the reef is a genuine feature of the potential.
  • A natural extension is to test the state-to-distribution network in a coarse-grained hypersonic flow code and compare bulk relaxation and dissociation behavior against direct trajectory results; if the networks hold up off-grid, they offer a practical route around the enormous number of trajectories needed for exhaustive state-to-state tables.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 4 minor

Summary. The manuscript reports a new global RKHS potential energy surface for the ground 1A' state of O3 based on MRCI+Q/aug-cc-pVTZ energies, and uses quasi-classical trajectory simulations to compute O(3P)+O2 exchange and dissociation rate coefficients. Comparisons are made with an existing PIP surface and with earlier SSB and DLLJG surfaces, and a neural-network state-to-distribution model is trained on QCT data for both new surfaces. The central scientific claim is that the experimentally observed negative temperature dependence of the exchange rate is reproduced on both RKHS and PIP surfaces even though both contain 'reef' features near dissociation, so the reef is not responsible for the positive temperature dependence previously associated with such features.

Significance. The exchange-rate comparison is a strong, falsifiable test: no parameter is fitted to the experimental rate data, and the RKHS and PIP results bracket the measured k(T). The paper also extends the authors' MRCI+Q/AVTZ family of reactive PESs to O+O2 and provides NN-based state-to-distribution emulators that are useful for coarse-grained hypersonic modeling. If the reef conclusion is correct, it would resolve a long-standing debate about the kinetic role of the near-dissociation reef in O3. However, the reef conclusion is not yet controlled quantitatively: the reef features are characterized only by qualitative 1D cuts, the RKHS surface underbinds O3 by about 4.4 kcal/mol relative to the accepted D_e, and the dissociation rates are two orders of magnitude below experiment. These gaps leave the central claim underdetermined and limit the present significance.

major comments (2)
  1. [Discussion and Conclusions; Fig. 3] The central claim that the reef is not responsible for the positive T-dependence of k_exch is underdetermined by the evidence presented. The argument is that two reef-bearing PESs (RKHS and PIP) give negative T-dependence, while the reef-free DLLJG also gives negative T-dependence. But the reef features are characterized only by 1D cuts at three fixed Jacobi angles, and no barrier heights, widths, or angular extents are reported. Fig. 3 in fact shows that the RKHS reef is absent at θ=105° while the text states that both PESs 'feature' reefs. The RKHS surface also underbinds O3 by 4.4 kcal/mol (D_e=21.7 vs 26.1 kcal/mol, Section 'Construction and Validation'), which distorts the near-dissociation region where the reef sits. An uncontrolled comparison of three qualitatively different electronic-structure levels cannot exclude the threshold scenario that only sufficiently pronounced reefs, such as that in SSB, cause positive T-dependence. I recommend quantifying the reef barriers along the minimum energy path for RKHS, PIP, and SSB, and preferably testing a single PES with and without the reef (for example by morphing) to make the claim load-bearing.
  2. [Dissociation Reaction; Eq. (2)] The abstract's statement that accounting for increased electronic degeneracy reduces the dissociation-rate discrepancy to one order of magnitude is not supported by the paper. In the Dissociation Reaction section, the manuscript proposes g_e=16/27 based on earlier suggestions, but no QCT rate with this value is computed, and the authors themselves use the word 'conceivable.' The reported rates are two orders of magnitude below experiment for both RKHS and PIP surfaces (Fig. 5 and Tables S1–S2), even though PIP has the correct O2 D_e of 120.6 kcal/mol. The persistent factor-of-about-100 shortfall means the high-energy dynamics on these surfaces are not quantitatively validated; since the reef discussion concerns exactly the near-dissociation topography, this quantitative failure weakens the confidence with which reef inferences can be drawn. The authors should either recompute k_diss with the revised degeneracy factor or explicitly label the one-order-of-magnitude reduction as a hypothesis.
minor comments (4)
  1. [Fig. 3 and Results, fourth paragraph] The sentence 'the green solid line do not display this feature' contains a subject-verb agreement error and leaves ambiguous whether the absence occurs only at θ=105° for the RKHS-PES; please rephrase and state explicitly which curves contain a submerged barrier.
  2. [QCT Simulations, Eq. (2)] The quantity b_max appears in the rate expression but is not defined in the text; please specify how b_max is determined and whether it is temperature- or energy-dependent.
  3. [Final State Distribution and NN-based Model] The NN-STD evaluation reports only the best, average, and worst cases, with worst-case R2 values of 0.69 and 0.77; please provide summary statistics for the full test set (mean, median, and range of R2) and state explicitly which initial-condition space is covered by the 34 excluded samples used in training.
  4. [Table 1] The ΔE column for MIN1 is empty; please add a footnote stating that the global minimum is set to zero and report D_e for each surface in a separate column so the reader can directly compare binding energies.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: QCT rates are predicted from ab initio PESs and compared with experiment; no target rate enters the PES construction, and the reef conclusion rests on dynamics rather than on a definitional identity.

full rationale

The paper's central claim is that the presence of a 'reef' near dissociation is not responsible for the positive temperature dependence of the O + O2 exchange rate. That claim is obtained by running quasi-classical trajectories on two independently constructed ab initio PESs (a new RKHS PES and the earlier PIP-PES) and comparing the resulting k_exch(T) with experiments. The RKHS PES is fitted only to MRCI+Q/A VTZ electronic energies, and the PIP-PES was published independently; no parameter of either PES is fitted to the experimental exchange or dissociation rate data. The QCT rate, Eq. (2), is a standard statistical expression that converts reactive trajectory counts into rates, and the degeneracy factors are taken from earlier published work, not tuned to reproduce the target rates. The conclusion about the reef is an empirical inference from the observation that two reef-bearing PESs reproduce the negative T-dependence, not a quantity that is equal to its inputs by construction. The NN state-to-distribution model is trained on QCT final-state distributions and validated on held-out QCT conditions, so it is an interpolation benchmark rather than a disguised restatement of the rates. The paper's acknowledged limitations, such as the RKHS PES underestimating De by about 4 kcal/mol and the dissociation rates being two orders of magnitude low, are accuracy concerns for the reef interpretation, not circularity. Likewise the comparison across PESs of different electronic-structure levels may underdetermine the reef conclusion, but underdetermination is not the same as deriving a result from its own definition. No load-bearing self-citation is used: references to the group's own RKHS toolkit, QCT sampling protocol, and prior STD model are methodological citations with independent technical content, and the central rate comparison is externally falsifiable against measured data. Therefore no circular step meeting the quoted-evidence standard can be identified.

Assumptions & free parameters 2 free parameters · 4 assumptions · 0 invented entities

The central claim rests on the accuracy of the PES in the near-dissociation region, the validity of QCT sampling, and literature degeneracy factors. The two free parameters listed are the only adjustable numbers that directly affect the computed rates, and neither is fitted to the experimental rate data.

free parameters (2)
  • RKHS switching parameters rho_j = 1.40 a0
    Determined by a grid-based search over the mixing dataset to minimize representation error in the channel-overlap regions (Eq. 1, Methods).
  • Electronic degeneracy factor for dissociation g_e = 1/27 (adopted from Ref. 31; authors suggest 16/27)
    Scales the absolute dissociation rates; the paper uses the literature value, which it later argues is too small by an order of magnitude.
assumptions (4)
  • domain assumption MRCI+Q with the A VTZ basis set provides a sufficiently accurate global PES for O3 dynamics despite a ~4 kcal/mol underestimate of the dissociation energy.
    The PES is built from these energies; the paper validates against rates but acknowledges the known De error in the Section 'Construction and Validation'.
  • domain assumption Quasi-classical trajectory sampling with Boltzmann sampling of J, b, Etrans, v, j gives reliable thermal rate coefficients for O+O2 at 100-500 K (exchange) and 1000-20000 K (dissociation).
    Standard QCT assumption; the paper does not provide a quantum benchmark on the same PES at these temperatures, only literature wavepacket results for other PESs.
  • domain assumption The electronic degeneracy factors ge(T)=3[5+3exp(-277.6/T)+exp(-325.9/T)] for exchange and ge=1/27 for dissociation are the correct statistical weights for the ground-state dynamics.
    Taken from Refs. 13,31,52,53; the dissociation factor is later questioned by the authors themselves.
  • ad hoc to paper The 'reef' feature identified in 1D cuts along R at fixed angles corresponds to the same feature discussed in prior literature, and its presence or absence can be compared across PESs with different electronic structure levels.
    The reef is characterized visually from fixed-angle cuts (Figure 3), not from a full minimum-energy path; different PESs differ in many other respects.

how reviews work

0 comments
Cite this review

Pith. "Pith review of High-Energy Reaction Dynamics of O$_3$." pith.science (2026). https://pith.science/paper/GNMCQOLE

@misc{pith2026250606088,
  author       = {Pith},
  title        = {Pith review of: High-Energy Reaction Dynamics of O$_3$},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/GNMCQOLE}},
  note         = {Machine review of arXiv:2506.06088}
}
abstract

The high-temperature atom exchange and dissociation reaction dynamics of the O($^3$P) + O$_2(^3\Sigma_g^{-} )$ system are investigated based on a new reproducing kernel-based representation of high-level multi-reference configuration interaction energies. Quasi-classical trajectory (QCT) simulations find the experimentally measured negative tempe-rature-dependence of the rate for the exchange reaction and describe the experiments within error bars. Similarly, QCT simulations for a recent potential energy surface (PES) at a comparable level of quantum chemical theory reproduce the negative $T-$dependence. Interestingly, both PESs feature a ``reef" structure near dissociation which has been implicated to be responsible for a positive $T-$dependence of the rate inconsistent with experiments. For the dissociation reaction the $T-$dependence correctly captures that known from experiments but underestimates the absolute rates by two orders of magnitude. Accounting for an increased number of accessible electronic states reduces this to one order of magnitude. A neural network-based state-to-distribution model is constructed for both PESs and shows good performance in predicting final translational, vibrational, and rotational product state distributions. Such models are valuable for future and more coarse-grained simulations of reactive hypersonic gas flow.

Figures

Figures reproduced from arXiv: 2506.06088 by the authors.

Figure 1
Figure 1. Validation of the RKHS-represented reactive, 3-dimensional PES for O+O [PITH_FULL_IMAGE:figures/full_fig_p011_1.png] view at source ↗
Figure 2
Figure 2. Contours for the O3 RKHS-PES and PIP-PES.27 Energies are in kcal/mol with the global minimum (MIN1) energy as the zero of energy. Panels A/B and C/D for the RKHS-PES and PIP-PES, respectively. In panels A and C the two bond lengths forming the angle α(OB-OA-OC) are identical. In panels B and D the bond length rBC = 2.283 a0 while scanning the other bond length. Red circles, triangles and crosses designate MIN1 (glob… view at source ↗
Figure 3
Figure 3. Comparison of the R−dependence V (R) of the RKHS-PES (solid lines), the PIP￾PES27 (dashed lines), PES220 (blue circles), and PES3 26 (magenta circles). PES2 features the “reef” whereas PES3 was designed to have the “reef” removed. The Jacobi angles θ for the RKHS and PIP-PES are 105◦ (green), 117◦ (red), and 128◦ (orange), respectively, see Figure 1A. Depending on the angle considered, the RKHS-PES features a “reef”… view at source ↗
Figures from the paper (6 more)
Figure 4
Figure 4. Figure 4: Comparison of k exch(T) for the O2 + O → O2 + O exchange reaction from the present QCT simulations using the RKHS (orange circles) and PIP27 (blue circles) PESs and earlier Jacobi coordinate-based quantum wave packet method simulations using the DLLJG16,17 (violet) and…
Figure 5
Figure 5. Figure 5: Rates k diss(T) for the dissociation reaction from QCT simulations. Simulation results using the RKHS- and PIP-PESs are the red and blue symbols including error bars from bootstrapping. The measurements are in black 69,70 together with a linear regression (dashed black…
Figure 6
Figure 6. Figure 6: Product QCT distribution for the 1A′ RKHS-PES from 3 different initial condi￾tions: [v = 2, j = 193, Etrans = 7.0 eV] (blue), [v = 34, j = 1, Etrans = 4.5 eV] (red) and [v = 38, j = 29, Etrans = 1.0 eV] (green). The final E ′ trans, E ′ int, v ′ , and j ′ are plotted a…
Figure 7
Figure 7. Figure 7: NN-STD Model for final state distributions from QCT simulations using the [PITH_FULL_IMAGE:figures/full_fig_p024_7.png]
Figure 8
Figure 8. Figure 8: NN-STD Model for final state distributions from QCT simulations using the PIP [PITH_FULL_IMAGE:figures/full_fig_p025_8.png]
Figure 9
Figure 9. Figure 9: The cumulative probability distributions [PITH_FULL_IMAGE:figures/full_fig_p027_9.png]

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

81 extracted references · 77 canonical work pages

  1. [1]

    Nonequilibrium Hypersonic Aerothermodynamics; Wiley-Interscience: New York, 1990

    Park, C. Nonequilibrium Hypersonic Aerothermodynamics; Wiley-Interscience: New York, 1990

  2. [2]

    Boyd, I. D. Computation of hypersonic flows using the direct simulation Monte Carlo method. J. Space. Rockets 2015, 52, 38--53

  3. [3]

    D.; Schwartzentruber, T

    Boyd, I. D.; Schwartzentruber, T. E. Nonequilibrium Gas Dynamics and Molecular Simulation; Cambridge University Press, 2017; Vol. 42

  4. [4]

    J.; Meuwly, M

    Koner, D.; Bemish, R. J.; Meuwly, M. Dynamics on multiple potential energy surfaces: Quantitative studies of elementary processes relevant to hypersonics. J. Phys.\ Chem. A 2020, 124, 6255--6269

  5. [5]

    Survey of flight and numerical data of hypersonic rarefied flows encountered in earth orbit and atmospheric reentry

    Schouler, M.; Pr \'e vereaud, Y.; Mieussens, L. Survey of flight and numerical data of hypersonic rarefied flows encountered in earth orbit and atmospheric reentry. Prog. Aerosp. Sci. 2020, 118, 100638

  6. [6]

    T.; Boe, K.; Bemish, R

    Koner, D.; Unke, O. T.; Boe, K.; Bemish, R. J.; Meuwly, M. Exhaustive state-to-state cross sections for reactive molecular collisions from importance sampling simulation and a neural network representation. J. Chem.\ Phys. 2019, 150, 211101

  7. [7]

    L.; Nandi, A.; Bowman, J

    Houston, P. L.; Nandi, A.; Bowman, J. M. A machine learning approach for prediction of rate constants. J. Phys.\ Chem. Lett. 2019, 10, 5250--5258

  8. [8]

    J.; Meuwly, M

    Arnold, J.; Koner, D.; K\"aser, S.; Singh, N.; Bemish, R. J.; Meuwly, M. Machine Learning for Observables: Reactant to Product State Distributions for Atom--Diatom Collisions. J. Phys.\ Chem. A 2020, 124, 7177--7190

Show all 81 references
  1. [9]

    C.; Koner, D.; Singh, N.; Bemish, R

    Arnold, J.; San Vicente Veliz, J. C.; Koner, D.; Singh, N.; Bemish, R. J.; Meuwly, M. Machine Learning Product State Distributions from Initial Reactant States for a Reactive Atom-Diatom Collision System. J. Chem.\ Phys. 2022, 156, 034301

  2. [10]

    Progress towards machine learning reaction rate constants

    Komp, E.; Janulaitis, N.; Valleau, S. Progress towards machine learning reaction rate constants. Phys.\ Chem.\ Chem.\ Phys. 2022, 24, 2692--2705

  3. [11]

    Inelastic N _2 + H _2 collisions and quantum-classical rate coefficients: large datasets and machine learning predictions

    Hong, Q.; Storchi, L.; Bartolomei, M.; Pirani, F.; Sun, Q.; Coletti, C. Inelastic N _2 + H _2 collisions and quantum-classical rate coefficients: large datasets and machine learning predictions. Eur. Phys. J. D 2023, 77, 128

  4. [12]

    State-to-state dynamics and machine learning predictions of inelastic and reactive O( ^3 P) + CO( ^1 ^+ ) collisions relevant to hypersonic flows

    Huang, X.; Cheng, X. State-to-state dynamics and machine learning predictions of inelastic and reactive O( ^3 P) + CO( ^1 ^+ ) collisions relevant to hypersonic flows. J. Chem.\ Phys. 2024, 160

  5. [13]

    Y.; Siebert, R.; Schinke, R.; Halberstadt, N

    Fleurat-Lessard, P.; Grebenshchikov, S. Y.; Siebert, R.; Schinke, R.; Halberstadt, N. Theoretical investigation of the temperature dependence of the O + O _2 exchange reaction. J. Chem.\ Phys. 2003, 118, 610--621

  6. [14]

    K.; Walker, R

    Babikov, D.; Kendrick, B. K.; Walker, R. B.; T Pack, R.; Fleurat-Lesard, P.; Schinke, R. Metastable states of ozone calculated on an accurate potential energy surface. J. Chem.\ Phys. 2003, 118, 6298--6308

  7. [15]

    Global permutationally invariant potential energy surface for ozone forming reaction

    Ayouz, M.; Babikov, D. Global permutationally invariant potential energy surface for ozone forming reaction. J. Chem.\ Phys. 2013, 138

  8. [16]

    Communication: An accurate global potential energy surface for the ground electronic state of ozone

    Dawes, R.; Lolur, P.; Li, A.; Jiang, B.; Guo, H. Communication: An accurate global potential energy surface for the ground electronic state of ozone. J. Chem.\ Phys. 2013, 139, 201103

  9. [17]

    Li, Y.; Sun, Z.; Jiang, B.; Xie, D.; Dawes, R.; Guo, H. Communication: Rigorous quantum dynamics of O + O _2 exchange reactions on an ab initio potential energy surface substantiate the negative temperature dependence of rate coefficients. J. Chem.\ Phys. 2014, 141, 081102

  10. [18]

    reef structure

    Tyuterev, V. G.; Kochanov, R.; Campargue, A.; Kassi, S.; Mondelain, D.; Barbe, A.; Starikova, E.; De Backer, M.; Szalay, P.; Tashkun, S. Does the “reef structure” at the ozone transition state towards the dissociation exist? New insight from calculations and ultrasensitive spe...

  11. [19]

    R.; Pack, R

    Hernandez-Lamoneda, R.; Salazar, M. R.; Pack, R. Does ozone have a barrier to dissociation and recombination? Chem.\ Phys.\ Lett. 2002, 355, 478--482

  12. [20]

    The transition-state region of the O( ^3 P) + O _2 ( ^3 _ g ^- ) potential energy surface

    Schinke, R.; Fleurat-Lessard, P. The transition-state region of the O( ^3 P) + O _2 ( ^3 _ g ^- ) potential energy surface. J. Chem.\ Phys. 2004, 121, 5789--5793

  13. [21]

    G.; M\"uller, T.; Tyuterev, V

    Holka, F.; Szalay, P. G.; M\"uller, T.; Tyuterev, V. G. Toward an improved ground state potential energy surface of ozone. J. Phys.\ Chem. A 2010, 114, 9927--9935

  14. [22]

    Spectroscopy of ozone at the dissociation threshold: Quantum calculations of bound and resonance states on a new global potential energy surface

    Siebert, R.; Schinke, R.; Bittererov \'a , M. Spectroscopy of ozone at the dissociation threshold: Quantum calculations of bound and resonance states on a new global potential energy surface. Phys.\ Chem.\ Chem.\ Phys. 2001, 3, 1795--1798

  15. [23]

    The vibrational energies of ozone up to the dissociation threshold: Dynamics calculations on an accurate potential energy surface

    Siebert, R.; Fleurat-Lessard, P.; Schinke, R.; Bittererov \'a , M.; Farantos, S. The vibrational energies of ozone up to the dissociation threshold: Dynamics calculations on an accurate potential energy surface. J. Chem. Phys. 2002, 116, 9749--9767

  16. [24]

    Potential energy surfaces of ozone in the ground state

    Ju-Xiang, S.; Zheng-He, Z.; Duo-Hui, H.; Jun, W.; Xin-Lu, C.; Xiang-Dong, Y. Potential energy surfaces of ozone in the ground state. Chin. Phys. 2007, 16, 2650

  17. [25]

    K.; Bhandarkar, U

    Mankodi, T. K.; Bhandarkar, U. V.; Puranik, B. P. Dissociation cross sections for N _2 + N 3N and O _2 + O 3O using the QCT method. J. Chem.\ Phys. 2017, 146

  18. [26]

    Communication: Highly accurate ozone formation potential and implications for kinetics

    Dawes, R.; Lolur, P.; Ma, J.; Guo, H. Communication: Highly accurate ozone formation potential and implications for kinetics . J. Chem.\ Phys. 2011, 135, 081102

  19. [27]

    Varga, Z.; Paukku, Y.; Truhlar, D. G. Potential energy surfaces for O + O _2 collisions . J. Chem.\ Phys. 2017, 147, 154312

  20. [28]

    B.; Guo, H.; Truhlar, D

    Shu, Y.; Akher, F. B.; Guo, H.; Truhlar, D. G. Parametrically Managed Activation Functions for Improved Global Potential Energy Surfaces for Six Coupled 5A ' States and Fourteen Coupled 3A ' States of O + O _2 . J. Phys.\ Chem. A 2024, 128, 1207--1217

  21. [29]

    Xie, D.; Guo, H.; Peterson, K. A. Accurate ab initio near-equilibrium potential energy and dipole moment functions of the ground electronic state of ozone. J. Chem.\ Phys. 2000, 112, 8378--8386

  22. [30]

    G.; Kochanov, R

    Tyuterev, V. G.; Kochanov, R. V.; Tashkun, S. A.; Holka, F.; Szalay, P. G. New analytical model for the ozone electronic ground state potential surface and accurate ab initio vibrational predictions at high energy range. J. Chem.\ Phys. 2013, 139

  23. [31]

    A.; Boyd, I

    Andrienko, D. A.; Boyd, I. D. Rovibrational energy transfer and dissociation in O _2 --O collisions. J. Chem.\ Phys. 2016, 144

  24. [32]

    Non equilibrium vibrational population and dissociation rates of oxygen in electrical discharges: the role of atoms and of the recombination process

    Cacciatore, M.; Capitelli, M.; Dilonardo, M. Non equilibrium vibrational population and dissociation rates of oxygen in electrical discharges: the role of atoms and of the recombination process. Beitr \"a ge aus der Plasmaphysik 1978, 18, 279--299

  25. [33]

    L.; Mar, K

    Van Wyngarden, A. L.; Mar, K. A.; Boering, K. A.; Lin, J. J.; Lee, Y. T.; Lin, S.-Y.; Guo, H.; Lendvay, G. Nonstatistical behavior of reactive scattering in the ^ 18 O + ^ 32 O _2 isotope exchange reaction. J. Am.\ Chem.\ Soc. 2007, 129, 2866--2870

  26. [34]

    A.; Zhang, J.; Minton, T

    Lahankar, S. A.; Zhang, J.; Minton, T. K.; Guo, H.; Lendvay, G. Dynamics of the O-atom exchange reaction ^ 16 O( ^3 P) + ^ 18 O ^ 18 O( ^3 _ g ^ - ) ^ 16 O ^ 18 O( ^3 _ g ^ - ) + ^ 18 O( ^3 P) at hyperthermal energies. J. Phys.\ Chem. A 2016, 120, 5348--5359

  27. [35]

    Y.; Schinke, R.; Janssen, C.; Krankowsky, D

    Fleurat-Lessard, P.; Grebenshchikov, S. Y.; Schinke, R.; Janssen, C.; Krankowsky, D. Isotope dependence of the O+ O _2 exchange reaction: Experiment and theory. J. Chem.\ Phys. 2003, 119, 4700--4712

  28. [36]

    Isotope dependence of the O + O _2 exchange reaction: Experiment and theory

    Fleurat-Lessard, P.; Grebenshchikov, S. Y.; Schinke, R.; Janssen, C.; Krankowsky, D. Erratum:“Isotope dependence of the O + O _2 exchange reaction: Experiment and theory”[J. Chem. Phys. 119, 4700 (2003)]. J. Chem.\ Phys. 2004, 120, 4993--4993

  29. [37]

    Mechanism change in the dynamics of the O ' + O _2 O ' O + O atom exchange reaction at high collision energies

    Lendvay, G. Mechanism change in the dynamics of the O ' + O _2 O ' O + O atom exchange reaction at high collision energies. J. Phys.\ Chem. A 2019, 123, 10230--10239

  30. [38]

    T.; Meuwly, M

    Unke, O. T.; Meuwly, M. Toolkit for the Construction of Reproducing Kernel-based Representations of Data: Application to Multidimensional Potential Energy Surfaces. J. Chem. Inf. Model. 2017, 57, 1923--1931

  31. [39]

    R.; Davidson, E

    Langhoff, S. R.; Davidson, E. R. Configuration interaction calculations on the nitrogen molecule. Int.\ J.\ Quant.\ Chem. 1974, 8, 61--72

  32. [40]

    Werner, H.; Knowles, P. J. An efficient internally contracted multiconfiguration–reference configuration interaction method. J. Chem.\ Phys. 1988, 89, 5803--5814

  33. [41]

    Gaussian basis sets for use in correlated molecular calculations

    Dunning, J., Thom H. Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen. J. Chem.\ Phys. 1989, 90, 1007--1023

  34. [42]

    Lin, W.; Varga, Z.; Song, G.; Paukku, Y.; Truhlar, D. G. Global triplet potential energy surfaces for the N _2 (X ^1 ) + O( ^3 P) → NO(X ^2 ) + N( ^4 S) reaction. J. Chem.\ Phys. 2016, 144, 024309

  35. [43]

    J.; Meuwly, M

    Koner, D.; Bemish, R. J.; Meuwly, M. The C( ^3 P) + NO(X ^2 ) O( ^3 P) + CN(X ^2 ^+ ), N( ^2 D)/N( ^4 S) + CO(X ^1 ^+ ) reaction: Rates, branching ratios, and final states from 15 K to 20 000 K . J. Chem.\ Phys. 2018 , 149 , 094305

  36. [44]

    C.; Koner, D.; Schwilk, M.; Bemish, R

    San Vicente Veliz, J. C.; Koner, D.; Schwilk, M.; Bemish, R. J.; Meuwly, M. The N( ^4 S)+ O _2 (X ^3 ) O( ^3 P)+ NO(X ^2 ) reaction: thermal and vibrational relaxation rates for the 2A ' , 4A ' and 2A '' states. Phys. Chem. Chem. Phys. 2020, 22, 3927--3939

  37. [45]

    C.; Koner, D.; Bieske, E

    Patra, S.; San Vicente Veliz, J. C.; Koner, D.; Bieske, E. J.; Meuwly, M. Photodissociation dynamics of N _3^+ . J. Chem.\ Phys. 2022, 156

  38. [46]

    Veliz, J. C. S. V.; Koner, D.; Schwilk, M.; Bemish, R. J.; Meuwly, M. The C( ^3 P)+ O _2 ( ^3 _ g ^- ) CO _2 CO( ^1 ^+ ) + O( ^1 D)/O( ^3 P) Reaction: Thermal and Vibrational Relaxation Rates from 15 K to 20000 K. Phys.\ Chem.\ Chem.\ Phys. 2021, 23, 11251--11263

  39. [47]

    C.; Koner, D.; Schwilk, M.; Bemish, R

    San Vicente Veliz, J. C.; Koner, D.; Schwilk, M.; Bemish, R. J.; Meuwly, M. The N( ^4 S) + O _2 (X ^ 3 ^-_g ) O( ^3 P) + NO(X ^ 2 ) Reaction: Thermal and Vibrational Relaxation Rates for the ^2 A ' , ^4 A ' and ^2 A '' States. Phys.\ Chem.\ Chem.\ Phys. 2020, 22, 3927--3939

  40. [48]

    C.; ; Bemish, R

    Koner, D.; San Vicente Veliz, J. C.; ; Bemish, R. J.; Meuwly, M. Accurate Reproducing Kernel-based Potential Energy Surfaces for the Triplet Ground States of N _2 O and Dynamics for the N + NO O + N _2 Reaction. Phys.\ Chem.\ Chem.\ Phys. 2020, 22, 18488--18498

  41. [49]

    u tz , M.; Celani, P.; Gy\

    Werner, H.-J.; Knowles, P. J.; Knizia, G.; Manby, F. R.; Sch\" u tz , M.; Celani, P.; Gy\"orffy, W.; Kats, D.; Korona, T.; Lindh, R. et al. MOLPRO, version 2020, a package of ab initio programs. 2020

  42. [50]

    almost size consistent

    Rintelman, J. M.; Adamovic, I.; Varganov, S.; Gordon, M. S. Multireference second-order perturbation theory: How size consistent is “almost size consistent”. J. Chem.\ Phys. 2005, 122, 044105

  43. [51]

    Explicitly correlated multireference configuration interaction: MRCI-F12

    Shiozaki, T.; Knizia, G.; Werner, H.-J. Explicitly correlated multireference configuration interaction: MRCI-F12 . J. Chem.\ Phys. 2011, 134, 034113

  44. [52]

    Gross, A.; Billing, G. D. Isotope effects on the rate constants for the processes O _2 + O O + O _2 and O _2 + O + Ar O _3 + Ar on a modified ground-state potential energy surface for ozone. Chem. Phys. 1997, 217, 1--18

  45. [53]

    A.; Boyd, I

    Andrienko, D. A.; Boyd, I. D. High fidelity modeling of thermal relaxation and dissociation of oxygen. Phys. Fluids 2015, 27, 116101

  46. [54]

    C.; Koner, D.; Singh, N.; Bemish, R

    Arnold, J.; San Vicente Veliz, J. C.; Koner, D.; Singh, N.; Bemish, R. J.; Meuwly, M. Machine learning product state distributions from initial reactant states for a reactive atom-diatom collision system . J. Chem.\ Phys. 2022, 156, 034301

  47. [55]

    A.; Bottou, L.; Orr, G

    LeCun, Y. A.; Bottou, L.; Orr, G. B.; M \"u ller, K.-R. Neural networks: Tricks of the trade; Springer, 2012; pp 9--48

  48. [56]

    Understanding the difficulty of training deep feedforward neural networks

    Glorot, X.; Bengio, Y. Understanding the difficulty of training deep feedforward neural networks. Proc. 13th Int. Conf. Artif. Intell. Stat. 2010; pp 249--256

  49. [57]

    Adam: A method for stochastic optimization

    Kingma, D.; Ba, J. Adam: A method for stochastic optimization. arXiv preprint arXiv:1412.6980 2014,

  50. [58]

    Abadi, M.; Barham, P.; Chen, J.; Chen, Z.; Davis, A.; Dean, J.; Devin, M.; Ghemawat, S.; Irving, G.; Isard, M. et al. Tensorflow: A system for large-scale machine learning. 12th USENIX symposium on operating systems Design and Implementation ( OSDI 16). 2016; pp 265--283

  51. [59]

    Schinke, R.; McBane, G. C. Photodissociation of ozone in the Hartley band: Potential energy surfaces, nonadiabatic couplings, and singlet/triplet branching ratio . J. Chem.\ Phys. 2010, 132, 044305

  52. [60]

    A climbing image nudged elastic band method for finding saddle points and minimum energy paths

    Henkelman, G.; Uberuaga, B.; Jonsson, H. A climbing image nudged elastic band method for finding saddle points and minimum energy paths. J. Chem.\ Phys. 2000, 113, 9901--9904

  53. [61]

    L.; Groves, M

    Kolsbjerg, E. L.; Groves, M. N.; Hammer, B. An automated nudged elastic band method . J. Chem.\ Phys. 2016 , 145 , 094107

  54. [62]

    Coriolis interaction and anharmonic potential function of ozone from the microwave spectra in the excited vibrational states

    Tanaka, T.; Morino, Y. Coriolis interaction and anharmonic potential function of ozone from the microwave spectra in the excited vibrational states. J.\ Mol.\ Struct. 1970, 33, 538--551

  55. [63]

    B.; Ignatov, S

    Gadzhiev, O. B.; Ignatov, S. K.; Kulikov, M. Y.; Feigin, A. M.; Razuvaev, A. G.; Sennikov, P. G.; Schrems, O. Structure, Energy, and Vibrational Frequencies of Oxygen Allotropes On ( n 6 ) in the Covalently Bound and van der Waals Forms: Ab Initio Study at the CCSD(T) Level. J...

  56. [64]

    Electronic structure and bonding of ozone

    Kalemos, A.; Mavridis, A. Electronic structure and bonding of ozone. J. Chem.\ Phys. 2008, 129, 054312

  57. [65]

    Theoretical Prediction on the Thermal Stability of Cyclic Ozone and Strong Oxygen Tunneling

    Chen, J.-L.; Hu, W.-P. Theoretical Prediction on the Thermal Stability of Cyclic Ozone and Strong Oxygen Tunneling. J. Am.\ Chem.\ Soc. 2011, 133, 16045--16053

  58. [66]

    W.; Davies, C

    Chase, M. W.; Davies, C. A.; Downey, J. R.; Frurip, D. J.; McDonald, R. A.; Syverud, A. N. JANAF Thermochemical Tables - 3rd edition .1. Al-Co . J. Phys.\ Chem. Ref. Data 1985 , 14 , 1--926

  59. [67]

    S.; Dachsel, H.; Harrison, R

    M \"u ller, T.; Xantheas, S. S.; Dachsel, H.; Harrison, R. J.; Nieplocha, J.; Shepard, R.; Kedziora, G. S.; Lischka, H. A systematic ab initio investigation of the open and ring structures of ozone. Chem.\ Phys.\ Lett. 1998, 293, 72--80

  60. [68]

    R.; Larsen, N

    Wiegell, M. R.; Larsen, N. W.; Pedersen, T.; Egsgaard, H. The temperature dependence of the exchange reaction between oxygen atoms and dioxygen molecules studied by means of isotopes and spectroscopy. Int. J. Chem. Kinet. 1997, 29, 745--753

  61. [69]

    Byron, S. R. Measurement of the Rate of Dissociation of Oxygen. J. Chem. Phys. 1959, 30, 1380--1392

  62. [70]

    Shatalov, O. P. Molecular dissociation of oxygen in the absence of vibrational equilibrium. Combust. Explos. Shock Waves 1973, 9, 610--613

  63. [71]

    Accurate ab initio potential energy curve of O _2

    Bytautas, L.; Ruedenberg, K. Accurate ab initio potential energy curve of O _2 . I. Nonrelativistic full configuration interaction valence correlation by the correlation energy extrapolation by intrinsic scaling method. J. Chem.\ Phys. 2010, 132, 074109

  64. [72]

    E.; Morton, M

    Ruscic, B.; Pinzon, R. E.; Morton, M. L.; von Laszevski, G.; Bittner, S. J.; Nijsure, S. G.; Amin, K. A.; Minkoff, M.; Wagner, A. F. Introduction to active thermochemical tables: Several “key” enthalpies of formation revisited. J. Phys.\ Chem. A 2004, 108, 9979--9997

  65. [73]

    A comprehensive re-analysis of the O _2 (B ^3 _ u ^- - X ^3 _ g ^- ) Schumann-Runge band system

    Creek, D.; Nicholls, R. A comprehensive re-analysis of the O _2 (B ^3 _ u ^- - X ^3 _ g ^- ) Schumann-Runge band system. Proc. R. Soc. London, Ser. A 1975, 341, 517--536

  66. [74]

    Nikitin, E. E. Theory of Elementary Atomic and Molecular Processes in Gases; Oxford University Press: Oxford, 1974

  67. [75]

    B.; Panesi, M

    Munafò, A.; Alberti, A.; Pantano, C.; Freund, J. B.; Panesi, M. A computational model for nanosecond pulse laser-plasma interactions. J. Comput. Phys. 2020, 406, 109190

  68. [76]

    C.; Meuwly, M

    Wang, J.; San Vicente Veliz, J. C.; Meuwly, M. High-Energy Reaction Dynamics of N _3 . J. Phys.\ Chem. A 2024, 128, 8322--8332

  69. [77]

    Meuwly, M.; Hutson, J. M. Morphing ab initio potentials: A systematic study of Ne--HF. J. Chem.\ Phys. 1999, 110, 8338--8347

  70. [78]

    P.; Vazquez-Salazar, L

    Horn, K. P.; Vazquez-Salazar, L. I.; Koch, C. P.; Meuwly, M. Improving potential energy surfaces using measured Feshbach resonance states. Sci. Adv. 2024, 10, eadi6462

  71. [79]

    R.; Glogover, S

    Chang, B.-Y.; Kung, C.-Y.; Kittrell, C.; Hsiao, C.-W.; Johnson, B. R.; Glogover, S. G.; Kinsey, J. L. High-accuracy measurement of vibrational Raman bands of ozone at 266 and 270 nm excitations. J. Chem.\ Phys. 1994, 101, 1914--1922

  72. [80]

    Transfer-learned potential energy surfaces: Toward microsecond-scale molecular dynamics simulations in the gas phase at CCSD (T) quality

    K \"a ser, S.; Meuwly, M. Transfer-learned potential energy surfaces: Toward microsecond-scale molecular dynamics simulations in the gas phase at CCSD (T) quality. J. Chem.\ Phys. 2023, 158, 214301

  73. [81]

    u;m 0[ c& Ϊ0_OX ]AB)dJL[|#Bm eOWvM >ȶzY 3 u ݽv V y[7tB =q `n M s J M!#2 / /x Jԝ:vW5

    Hickson, K. M.; Veliz, J. C. S. V.; Koner, D.; Meuwly, M. Low-temperature kinetics for the N + NO reaction: experiment guides the way. Phys.\ Chem.\ Chem.\ Phys. 2023, 25, 13854--13863 mcitethebibliography arxiv.tex0000664000000000000000000016523615020567762011460 0ustar rootr...

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.