Pith. sign in

REVIEW 3 major objections 5 minor 13 references

Non-adiabatically driven quantum interference effects in the ultracold K + KRb $\longrightarrow$ Rb + K$_{2}$ chemical reaction

T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Non-adiabatic coupling to an excited electronic state brings the ultracold K + KRb reaction rate into line with experiment.

desk verdict A solid first non-adiabatic calculation for K+KRb that improves agreement with experiment, but the hyperfine-free comparison makes the headline claim less clean than it looks. read the letter →

arxiv 2506.01831 v1 pith:A5PAY6DZ submitted 2025-06-02 physics.chem-ph physics.atm-clusquant-ph

H. da Silva Jr. , B. K. Kendrick , H. Li , S. Kotochigova , N. Balakrishnan This is my paper · ORCID
classification physics.chem-phphysics.atm-clusquant-ph PACS 34.50.-s34.50.Lf
keywords ultracoldchemistrynon-adiabaticquantumdynamicsconicalintersectiongeometricphaseK+KRbreactionratecoefficientsinterferencecoupled-channelscattering
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

The paper reports the first non-adiabatic quantum scattering study of the ultracold K + KRb → Rb + K2 reaction, the first atom-diatom reaction measured below 1 μK. Previous single-surface calculations that ignored the first excited electronic state produced a rate coefficient about 35% below the measured value. By coupling the ground and first excited electronic states through a 2×2 diabatic treatment, the authors obtain a zero-temperature total rate of roughly 1.32×$10^{-10}$ cm3/s, close to the lower edge of the measured 1.7±0.3×$10^{-10}$ cm3/s. They argue that short-range dynamics near a conical intersection introduces quantum interference that changes both the overall rate and product state-to-state distributions, and that such non-adiabatic effects are confined to collision energies below about a millikelvin.

What carries the argument

The central object is a 2×2 diabatic representation of the time-independent Schrödinger equation, written in hyperspherical coordinates, in which the ground and first excited electronic states become the diagonal potentials and are connected by an off-diagonal coupling. The diabatic potentials are built by transforming ab initio adiabatic surfaces using an adiabatic-to-diabatic transformation angle obtained by integrating the non-adiabatic coupling term. In the region where the excited-state channels are closed but strongly coupled, this 2×2 system is propagated, then transformed back to a single-surface form at a large hyperradius. Physically, the scattering amplitude splits into a direct path and a path encircling the conical intersection; the non-adiabatic treatment changes the sign of their interference cross term, converting constructive interference into destructive interference and vice versa.

What would settle it

A hyperfine-resolved scattering calculation or measurement using the same initial hyperfine states that changes the zero-temperature rate by roughly 40% or more would show that the agreement is not caused by the non-adiabatic coupling. Alternatively, a product-state-resolved measurement of the K2 rotational distribution below 1 μK, if it matched the single-surface prediction rather than the two-state interference pattern, would rule out the claimed mechanism.

Watch

Extended reading notes

Core claim

The central claim is that the long-standing theory-experiment gap for the ultracold K + KRb reaction is mostly accounted for by non-adiabatic coupling to the first excited electronic state, which is energetically accessible even in the zero-temperature limit. In a conventional Born-Oppenheimer no-geometric-phase calculation the total rate at 250 nK is about 0.9×$10^{-10}$ cm3/s; including the excited state in a two-state coupled-channel calculation raises it to about 1.32×$10^{-10}$ cm3/s, within the experimental error bars of the 2010 measurement. The authors further claim that the coupling flips the sign of the quantum interference term between the direct reactive path and the path that loops around the conical intersection, so the effect is genuinely quantum mechanical and grows as the temperature approaches zero.

Load-bearing premise

The calculation omits hyperfine structure in the entrance channel, so it assumes the measured rate, taken with K and KRb in specific hyperfine states, can be compared to a scattering calculation without any hyperfine resolution; if hyperfine effects shift the rate by tens of percent, the improved agreement could be coincidental.

Editorial extensions

If this is right

  • At zero temperature the 2×2 total rate converges to about 1.32×10^-10 cm3/s, versus about 0.9×10^-10 cm3/s in the no-geometric-phase calculation, placing theory at the lower edge of the experimental 1.7±0.3×10^-10 cm3/s.
  • Above roughly 1 mK the two-state and single-surface calculations agree, so non-adiabatic interference is an ultracold-limit phenomenon tied to the Wigner threshold regime.
  • The products remain dominated by highly rotationally excited K2 in the v'=0 and 1 manifolds, meaning the main non-adiabatic effect is on the magnitude and interference pattern of the rate rather than on the gross energy disposal.
  • Short-range coupling to the excited electronic state, not only long-range forces, must be included to reproduce ultracold reactive rates, with consequences for other near-universal-loss bialkali systems.

Reading between the lines

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

  • The authors' own caveat implies the natural next test is a hyperfine-resolved version of this calculation; until that exists, the agreement with experiment should be read as conditional.
  • If the mechanism is general, other ultracold alkali atom + alkali dimer reactions with accessible excited states may show similar non-adiabatic corrections, and the present treatment offers a template for computing them.
  • A product-state-resolved ultracold experiment on K + KRb could distinguish the two-state interference pattern from the single-surface one, giving a direct observable signature of the geometric-phase sign flip.
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

3 major / 5 minor

Summary. The manuscript reports the first non-adiabatic quantum reactive scattering study of the ultracold K + KRb(X 1Σ+, v=0, j=0) -> Rb + K2(X 1Σg+, v', j') reaction. The authors construct new diabatic potential energy surfaces for the lowest two electronic states of K2Rb from MRCI/ECP calculations, propagate a 2x2 diabatic coupled-channel system in hyperspherical coordinates, and compare the resulting total rate coefficients with the NGP (no-geometric-phase) calculation of Croft et al. and with the experimental value of Ospelkaus et al. They obtain a zero-temperature 2x2 total rate of about 1.32 x 10^-10 cm3/s, closer to the measured 1.7 ± 0.3 x 10^-10 cm3/s than the previous NGP value (0.9–1.1 x 10^-10 cm3/s), and attribute the improvement to non-adiabatic coupling and quantum interference near a conical intersection. The authors explicitly acknowledge that the calculation does not include hyperfine structure in the entrance channel.

Significance. If the calculation is correct, it would be the first non-adiabatic treatment of this benchmark ultracold atom–diatom reaction and would indicate that short-range non-adiabatic effects can substantially close a long-standing theory–experiment gap. The work has notable strengths: no parameter is fitted to the experimental rate, the rate emerges from ab initio PESs and coupled-channel scattering, and the authors reproduce the earlier NGP calculation to within about 20%. However, the central quantitative claim is not yet fully established: the calculation omits hyperfine structure, the authors themselves cite experiments showing hyperfine-dependent effects that 'cannot be neglected', and the claimed non-adiabatic correction is comparable in size to the experimental uncertainty. In addition, the specific claim that non-adiabatic interference changes state-to-state product distributions is not directly demonstrated by any NGP-versus-2x2 state-to-state comparison.

major comments (3)
  1. [Penultimate paragraph (hyperfine caveat)] The central comparison with the Ospelkaus et al. measurement is made against a scattering calculation that omits all hyperfine structure, as the authors state: 'it does not yet account for the hyperfine structure within the entrance channel.' At ultracold temperatures the total loss rate is controlled by the short-range quantum boundary condition, and the cited K + NaK and Rb + KRb experiments show hyperfine-dependent collisional outcomes of order unity. The non-adiabatic correction claimed here is an increase from 0.9 x 10^-10 to 1.32 x 10^-10 cm3/s, i.e. about 0.4 x 10^-10 cm3/s, which is only slightly larger than the experimental error bar of ±0.3 x 10^-10 cm3/s. Without a quantitative estimate or at least a bound on the hyperfine shift for the specific K + KRb initial state used in the experiment, the improved agreement cannot be uniquely attributed to non-adiabatic coupling; the paper's main conclusion is therefore not yet established.
  2. [Abstract and Fig. 4] The manuscript claims that 'short-range dynamics mediated by coupling with the excited electronic state introduces quantum interference effects that influence both the state-to-state rate coefficients and the overall reaction rates.' However, no NGP-versus-2x2 comparison of state-to-state rate coefficients is presented. Figure 3 shows only the 2x2 rotationally resolved rates at 211 nK, while Figure 4 compares only total rates. To support the state-to-state interference claim, the authors should provide state-to-state NGP and 2x2 rate coefficients (or an equivalent decomposition of the scattering amplitude into direct, loop, and interference terms) and show specifically where the interference changes sign or redistributes population.
  3. ['As a first step' paragraph and Fig. 4] The NGP baseline computed in this work is 0.9 x 10^-10 cm3/s at 250 nK, whereas the original Croft et al. value is 1.1 x 10^-10 cm3/s; the authors attribute this 18% discrepancy to differences in the adiabatic PES V1. This baseline uncertainty is roughly half the size of the claimed non-adiabatic correction (about 0.42 x 10^-10 cm3/s). The quantitative statement of 'better agreement with experiment' would be considerably strengthened by an error estimate for the present PESs or by repeating the 2x2 calculation with the original Croft et al. V1 to separate PES sensitivity from genuine non-adiabatic effects.
minor comments (5)
  1. [Page 6, electronic-structure paragraph] The text reads 'the NATC matrix was evaluated'; this should be 'NACT matrix' (non-adiabatic coupling term).
  2. [Title and abstract] The notation 'KRb−' appears to be a typographical artifact; the reaction is K + KRb, not K + KRb−. The minus sign is also inconsistent with the body of the paper.
  3. [Page 6, basis-set description] 'spd fdiffuse functions' should read 'spdf diffuse functions'.
  4. [Reference 45] The author name 'Fraschio' appears to be misspelled; please verify against the original publication.
  5. [Page 7, discussion of Fig. 2] The sentence 'presents every hundred V1 effective coupled potential curves ... compared to the first and the hundredth corresponding V2 curves' is unclear; please specify exactly which curves are shown and whether every 100th curve of V1 is plotted.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the non-adiabatic rate coefficients are computed from ab initio PESs and coupled-channel scattering, with no parameter fitted to the experimental datum.

full rationale

The paper's central claim is that including the first excited electronic state through a 2x2 diabatic coupled-channel calculation raises the zero-temperature K + KRb rate coefficient to about 1.32 x 10^-10 cm3/s, thereby improving agreement with the measured 1.7 +/- 0.3 x 10^-10 cm3/s. This is a direct dynamical result: the diabatic potentials are constructed from freshly computed adiabatic PESs and ab initio NACT/ADT angles (Eq. 1 and the subsequent definitions of V11, V22, V12), and the scattering matrix is obtained by numerically propagating the Schrodinger equation. No parameter is adjusted to reproduce the Ospelkaus measurement, and the comparison to experiment is used only as an external benchmark. Self-citations to Kendrick's APH3D non-adiabatic scattering formalism and to the earlier Li + LiNa study are method inheritance, not circularity: those methods have been validated on independent systems and the present calculation is externally falsifiable, for instance by hyperfine-resolved or temperature-resolved measurements. The authors' explicit admission that the calculation omits hyperfine structure is a correctness and uncertainty concern, not a circularity, because the theoretical prediction would have identical content even absent the experimental number. No equation is defined in terms of its own predicted outcome, and no fitted parameter is renamed as a prediction. Therefore the paper derives its rates from its stated inputs without any circular step.

Assumptions & free parameters 0 free parameters · 6 assumptions · 0 invented entities

The calculation does not fit any parameter to the experimental rate; the numerical result is an output of ab initio surfaces and scattering. The main unverified inputs are the two-state truncation, the accuracy of the new MRCI/ECP surfaces and NACT-derived diabatic coupling, and the neglect of hyperfine structure.

assumptions (6)
  • domain assumption The lowest two adiabatic states of K2Rb capture all relevant non-adiabatic dynamics; higher electronic states are negligible.
    The 2x2 diabatic Hamiltonian in Eq. 1 contains only two surfaces, and no test with a third state is reported. The authors state the excited state is accessible but do not quantify contributions from any further states.
  • domain assumption The adiabatic-to-diabatic transformation angles beta from numerical NACT integration are accurate and lead to well-defined diabatic potentials.
    The diabatic surfaces V11, V22, and V12 are constructed from V1, V2, and beta; no test of path independence or gauge consistency of the NACT integration is shown.
  • domain assumption The new MRCI/ECP/CPP potential surfaces are accurate enough to determine the conical intersection region and the short-range coupling.
    The PESs are new, and the only indirect check is the NGP rate of 0.9 x 10^-10 versus Croft's 1.1 x 10^-10 cm3/s. No direct comparison of V1 with the prior Croft PES is given.
  • domain assumption Neglecting hyperfine structure does not change the rate coefficient enough to invalidate the comparison with the Ospelkaus experiment.
    The authors state hyperfine structure is not included and cite experiments showing hyperfine-dependent outcomes 'cannot be neglected' in related systems. This is the most fragile assumption for the comparison to experiment.
  • domain assumption The hyperspherical coupled-channel scattering computation is numerically converged at the chosen basis sizes and propagation radii.
    The authors report convergence tests for rho_max (300 vs 325 a.u.) and rho_match (38.1, 49.9, 60.3 a.u.) but not for the channel basis size or the PES grid resolution.
  • standard math The coupled-channel scattering formalism in hyperspherical coordinates is a correct representation of the two-state Schrodinger equation.
    This is the established formalism of reference 26 and is not in dispute; the paper relies on it without re-derivation.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Non-adiabatically driven quantum interference effects in the ultracold K + KRb $\longrightarrow$ Rb + K$_{2}$ chemical reaction." pith.science (2026). https://pith.science/paper/A5PAY6DZ

@misc{pith2026250601831,
  author       = {Pith},
  title        = {Pith review of: Non-adiabatically driven quantum interference effects in the ultracold K + KRb $\longrightarrow$ Rb + K$_2$ chemical reaction},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/A5PAY6DZ}},
  note         = {Machine review of arXiv:2506.01831}
}
abstract

The K + KRb $\longrightarrow$ Rb + K$_{2}$ chemical reaction is the first ultracold atom-diatom chemical reaction for which experimental results have been reported for temperatures below 1 $\mu$K more than a decade ago. The reaction occurs through coupling with an excited electronic state that is accessible even in the ultracold limit. A previous quantum dynamics study, excluding non-adiabatic effects, has reported a rate coefficient that is about 35\% below the experimental value. Here, we report the first non-adiabatic quantum dynamics study of this reaction and obtain rate coefficients in better agreement with experiments. Our results show that short-range dynamics mediated by coupling with the excited electronic state introduces quantum interference effects that influence both the state-to-state rate coefficients and the overall reaction rates.

Figures

Figures reproduced from arXiv: 2506.01831 by the authors.

Figure 1
Figure 1. Global PESs V1 (bottom) and V2 (top), in cm−1 , for the K + KRb(X1Σ +) arrangement as functions of the Jacobi coordinates (R, θ) at a fixed KRb internuclear distance of r = 7.7 atomic units. The same energy scale and color coding are used on both surfaces. Isolines varying every 500 cm−1 from -6500 to 2000 cm−1 are displayed. Finally, the global diabatic PESs and the off-diagonal coupling surface are constructed by … view at source ↗
Figure 2
Figure 2. A subset of the lowest APH effective coupled potential curves (cm [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. State-to-state rate constants (cm3/s) for the K + KRb(X1Σ +, υ = 0, j = 0) −→ Rb + K2 [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: Total rate constants (cm3/s) as functions of the collision energy (K). Red curves are for the 2×2 results computed at the displayed values of ρdat, ρmatch, and ρmax, all in a.u. The green stripe tags the experimental value of Ospelkaus et al.18 within its error bar and…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

13 extracted references · 13 canonical work pages

  1. [1]

    (1) Krems, R. V. Cold controlled chemistry.Phys. Chem. Chem. Phys. 2008, 10, 4079–4092. (2) Carr, L. D.; Ye, J. Focus on Cold and Ultracold Molecules.New J. Phys. 2009, 11, 055009. (3) Kendrick, B. K.; Hazra, J.; Balakrishnan, N. The geometric phase controls ultracold chem- istry. Nature Communications 2015, 6,

    work page 2008

  2. [6]

    Chemical Physics 2018, 515, 387–399

    → H+H2(v′, j′). Chemical Physics 2018, 515, 387–399. (15) Kendrick, B. K. Nonadiabatic Ultracold Quantum Reactive Scattering of Hydrogen with Vibrationally Excited HD(v = 5 − 9). The Journal of Physical Chemistry A 2019, 123, 9919–9933. (16) Hazra, J.; Kendrick, B. K.; Balakrishnan, N. Importance of Geometric Phase Effects in Ultracold Chemistry.The Journ...

    work page 2018

  3. [8]

    K.; Balakrishnan, N

    (36) da Silva, H.; Kendrick, B. K.; Balakrishnan, N. On the use of stereodynamical effects to control cold chemical reactions: The H + D2 → D + HD case study.The Journal of Chemical Physics 2022, 156, 044305. (37) da Silva, H.; Yao, Q.; Morita, M.; Kendrick, B. K.; Guo, H.; Balakrishnan, N. The Li + CaF → Ca + LiF chemical reaction under cold conditions.P...

    work page 2022

  4. [73]

    K.; Hazra, J.; Balakrishnan, N

    (10) Kendrick, B. K.; Hazra, J.; Balakrishnan, N. Geometric Phase Appears in the Ultracold Hydrogen Exchange Reaction.Phys. Rev. Lett. 2015, 115, 153201. (11) Hazra, J.; Kendrick, B. K.; Balakrishnan, N. Geometric phase effects in ultracold hydrogen exchange reaction.Journal of Physics B: Atomic, Molecular and Optical Physics 2016, 49, 194004. (12) Kendri...

    work page 2015

  5. [113]

    W.; Will, S

    (55) Park, J. W.; Will, S. A.; Zwierlein, M. W. Ultracold Dipolar Gas of Fermionic 23Na40K Molecules in Their Absolute Ground State.Phys. Rev. Lett. 2015, 114, 205302. (56) Voges, K. K.; Gersema, P.; Meyer zum Alten Borgloh, M.; Schulze, T. A.; Hartmann, T.; Zenesini, A.; Ospelkaus, S. Ultracold Gas of Bosonic23Na39K Ground-State Molecules.Phys. Rev. Lett...

    work page 2015

  6. [142]

    On the reliability of semi-empirical pseudopotentials: simulation of Hartree-Fock and Dirac-Fock results.Journal of Physics B: Atomic and Molecular Physics 1983, 16, L323

    (44) Fuentealba, P.; Stoll, H.; von Szentpaly, L.; Schwerdtfeger, P.; Preuss, H. On the reliability of semi-empirical pseudopotentials: simulation of Hartree-Fock and Dirac-Fock results.Journal of Physics B: Atomic and Molecular Physics 1983, 16, L323. 18 (45) Fuentealba, P.; Szentpály, L.; Stoll, H.; Fraschio, F. X.; Preuss, H. Pseudopotential calcu- lat...

    work page 1983

  7. [145]

    L.; Rey, A

    (5) Bohn, J. L.; Rey, A. M.; Ye, J. Cold molecules: Progress in quantum engineering of chemistry and quantum matter.Science 2017, 357,

    work page 2017

  8. [853]

    X.; Zhu, L.; Luke, J.; Houwman, J

    (19) Liu, Y. X.; Zhu, L.; Luke, J.; Houwman, J. J. A.; Babin, M. C.; Hu, M. G.; Ni, K. K. Quantum interference in atom-exchange reactions.Science 2024, 384,

    work page 2024

Show all 13 references
  1. [1002]

    D.; Lin, Y.-W.; Gheorghe, A

    (6) Hu, M.-G.; Liu, Y.; Grimes, D. D.; Lin, Y.-W.; Gheorghe, A. H.; Vexiau, R.; Bouloufa- Maafa, N.; Dulieu, O.; Rosenband, T.; Ni, K.-K. Direct observation of bimolecular reactions of ultracold KRb molecules.Science 2019, 366, 1111–1115. (7) Liu, Y.; Hu, M.-G.; Nichols, M. A....

    2019

  2. [1117]

    (20) Croft, J. F. E.; Balakrishnan, N.; Kendrick, B. K. Long-lived complexes and signatures of chaos in ultracoldK2+Rb collisions.Phys. Rev. A 2017, 96, 062707. (21) Yarkony, D. R. Diabolical conical intersections.Rev. Mod. Phys. 1996, 68, 985–1013. (22) Domcke, W.; Yarkony, D...

    2017

  3. [3104]

    K.; Meyer zum Alten Borgloh, M.; Koch, L.; Hartmann, T.; Zen- esini, A.; Ospelkaus, S.; Lin, J.; He, J.; Wang, D

    (60) Gersema, P.; Voges, K. K.; Meyer zum Alten Borgloh, M.; Koch, L.; Hartmann, T.; Zen- esini, A.; Ospelkaus, S.; Lin, J.; He, J.; Wang, D. Probing Photoinduced Two-Body Loss of Ultracold Nonreactive Bosonic23Na87Rb and 23Na39K Molecules. Phys.Rev. Lett. 2021, 127, 163401. (...

    2021

  4. [3888]

    (40) Delves, L. M. Tertiary and general-order collisions.Nuclear Physics 1958, 9, 391–399. (41) Delves, L. M. Tertiary and general-order collisions (II).Nuclear Physics 1960, 20, 275–308. (42) Johnson, B. The multichannel log-derivative method for scattering calculations.Journ...

    1958

  5. [7918]

    Perspective: Ultracold molecules and the dawn of cold controlled chemistry

    14 (4) Balakrishnan, N. Perspective: Ultracold molecules and the dawn of cold controlled chemistry. The Journal of Chemical Physics 2016,

    2016

Pith tools

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