Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

The Journal of Physical Chemistry A

Global Full-Dimensional Potential Energy Surface for the Reaction 23Na87Rb + 23Na87Rb → 23Na2 + 87Rb2 and the Formation Rate and Lifetime of the 23Na287Rb2 Collision Complex

Yilang Liu, Jing Huang, Dongzheng Yang, Daiqian Xie, Hua Guo

Published Nov 23, 2022• 9 citations
Publisher

Abstract

A full-dimensional global potential energy surface (PES) for the reaction 23 Na 87 Rb + 23 Na 87 Rb → 23 Na 2 + 87 Rb 2 was constructed based on high-level ab initio calculations. The short-range part was expressed as a permutation invariant polynomial-neural network (PIP-NN) fit of 22 003 ab initio points calculated using a coupled cluster method with the one-electron basis 5 s 5 p 5 d 2 f plus effective core potentials and core polarizability potentials, while the long-range part was represented in an asymptotically correct form based on multipole expansion. The formation rate of the 23 Na 2 87 Rb 2 complex calculated using a quantum statistical method is in good agreement with experiment, while the estimated 19.20 μs lifetime of the complex from Rice–Ramsperger–Kassel–Marcus (RRKM) theory is significantly shorter than the measured millisecond decay rate, signaling either the inadequacy of RRKM theory or a yet unresolved loss mechanism.

DOI
10.1021/acs.jpca.2c06438

Related papers

Full-Dimensional Global Potential Energy Surface for the KRb + KRb → K2Rb2* → K2 + Rb2 Reaction with Accurate Long-Range Interactions and Quantum Statistical Calculation of the Product State Distribution under Ultracold Conditions

Jing Huang, Dongzheng Yang, Junxiang Zuo +3 more

A full-dimensional global potential energy surface (PES) for the KRb + KRb → K 2 Rb 2 * → K 2 + Rb 2 reaction is reported based on high-level ab initio calculations. The short-range part of the PES is fit with the permutationally invariant polynomial-neural network method, while the long-range parts of the PES in both the reactant and product asymptotes are represented by an asymptotically correct form. The long- and short-range parts are connected with intermediate-range parts to make them smooth. Within a statistical quantum model, this PES reproduces both the measured loss rates of ultracold KRb molecules and the K 2 and Rb 2 product state distributions, underscoring the important role of tunneling in ultracold chemistry. The PES also correctly predicts the lifetime of the K 2 Rb 2 * intermediate complex within the Rice–Ramsperger–Kassel–Marcus limit. It thus provides a reliable platform for future dynamical studies of the prototypical reaction.

The Journal of Physical Chemistry A•Jul 12, 2021•9 citations

Globally Accurate Gaussian Process Potential Energy Surface and Quantum Dynamics Studies on the Li(2S) + Na2 → LiNa + Na Reaction at Low Collision Energies

Zijiang Yang, Hanghang Chen, Bayaer Buren +1 more

The LiNa2 reactive system has recently received great attention in the experimental study of ultracold chemical reactions, but the corresponding theoretical calculations have not been carried out. Here, we report the first globally accurate ground-state LiNa2 potential energy surface (PES) using a Gaussian process model based on only 1776 actively selected high-level ab initio training points. The constructed PES had high precision and strong generalization capability. On the new PES, the quantum dynamics calculations on the Li(2S) + Na2(v = 0, j = 0) → LiNa + Na reaction were carried out in the 0.001–0.01 eV collision energy range using an improved time-dependent wave packet method. The calculated results indicate that this reaction is dominated by a complex-forming mechanism at low collision energies. The presented dynamics data provide guidance for experimental research, and the newly constructed PES could be further used for ultracold reaction dynamics calculations on this reactive system.

Molecules•Mar 24, 2023•7 citations

Quantum Dynamics Studies of the Li + Na2 (V = 0, j = 0) → Na + NaLi Reaction on a New Neural Network Potential Energy Surface

Bayaer Buren, Jiapeng Zhang, Yongqing Li

The ultracold reaction offers a unique opportunity to elucidate the intricate microscopic mechanism of chemical reactions, and the Na 2 Li system serves as a pivotal reaction system in the investigation of ultracold reactions. In this work, a high-precision potential energy surface (PES) of the Na 2 Li system is constructed based on high-level ab initio energy points and the neural network (NN) method, and a proper asymptotic functional form is adopted for the long-range interaction, which is suitable for the study of cold or ultracold collisions. Based on the new NN PES, the dynamics of the Li + Na 2 ( v = 0, j = 0) → Na + NaLi reaction are studied in the collision energy range of 10 –7 to 80 cm –1 . In the high collision energy range of 8 to 80 cm –1, the dynamics of the reaction is studied using the time-dependent wave packet method and the statistical quantum mechanical (SQM) method. Comparing the results of the two methods, it is found that the SQM method provides a rough description of the product ro-vibrational state distribution but overestimates the integral cross-section values. With the decrease of collision energy, the reaction differential cross section gradually changes from forward–backward symmetric scattering to predominant forward scattering. In the low collision energy range from 10 –7 to 8 cm –1, the SQM method is used to study the reaction dynamics, and the rate constant in the Wigner threshold region is estimated to be 2.87 × 10 –10 cm 3 /s.

The Journal of Physical Chemistry A•Jun 18, 2024•5 citations

Six-dimensional potential energy surface for NaK–NaK collisions: Gaussian process representation with correct asymptotic form

Arthur Christianen, Tijs Karman, Rodrigo A. Vargas–Hernández +2 more

Constructing accurate global potential energy surfaces (PESs) describing chemically reactive molecule-molecule collisions of alkali metal dimers presents a major challenge. To be suitable for quantum scattering calculations, such PESs must represent accurately three- and four-body interactions, describe conical intersections, and have a proper asymptotic form at the long range. Here, we demonstrate that such global potentials can be obtained by Gaussian Process (GP) regression merged with the analytic asymptotic expansions at the long range. We propose an efficient sampling technique, which allows us to construct an accurate global PES accounting for different chemical arrangements with <2500 ab initio calculations. We apply this method to (NaK)2 and obtain the first global PES for a system of four alkali metal atoms. The resulting surface exhibits a complex landscape including a pair and a quartet of symmetrically equivalent local minima and a seam of conical intersections. The dissociation energy found from our ab initio calculations is 4534 cm−1. This result is reproduced by the GP models with an error of less than 3%. The GP models of the PES allow us to analyze the features of the global PES, representative of general alkali metal four-atom interactions. Understanding these interactions is of key importance in the field of ultracold chemistry.

The Journal of Chemical Physics•Feb 13, 2019•51 citations

Universality and chaoticity in ultracold K+KRb chemical reactions

James F. E. Croft, Constantinos Makrides, Ming Li +4 more

Abstract A fundamental question in the study of chemical reactions is how reactions proceed at a collision energy close to absolute zero. This question is no longer hypothetical: quantum degenerate gases of atoms and molecules can now be created at temperatures lower than a few tens of nanokelvin. Here we consider the benchmark ultracold reaction between, the most-celebrated ultracold molecule, KRb and K. We map out an accurate ab initio ground-state potential energy surface of the K 2 Rb complex in full dimensionality and report numerically-exact quantum-mechanical reaction dynamics. The distribution of rotationally resolved rates is shown to be Poissonian. An analysis of the hyperspherical adiabatic potential curves explains this statistical character revealing a chaotic distribution for the short-range collision complex that plays a key role in governing the reaction outcome.

Nature Communications•Jul 19, 2017•81 citations

Quasiclassical method for calculating the density of states of ultracold collision complexes

Arthur Christianen, Tijs Karman, Gerrit C. Groenenboom

We derive a quasiclassical expression for the density of states (DOS) of an arbitrary, ultracold, $N$-atom collision complex, for a general potential energy surface (PES). We establish the accuracy of our quasiclassical method by comparing to exact quantum results for the ${\mathrm{K}}_{2}\text{\ensuremath{-}}\mathrm{Rb}$ and NaK-NaK systems, with isotropic model PESs. Next, we calculate the DOS for an accurate NaK-NaK PES to be 0.124 $\ensuremath{\mu}{\mathrm{K}}^{\ensuremath{-}1}$, with an associated Rice-Ramsperger-Kassel-Marcus sticking time of 6.0 $\ensuremath{\mu}\mathrm{s}$. We extrapolate the DOS and sticking times to all other polar bialkali-bialkali collision complexes by scaling with atomic masses, equilibrium bond lengths, dissociation energies, and dispersion coefficients. The sticking times calculated here are two to three orders of magnitude shorter than those reported by Mayle et al. [Phys. Rev. A 85, 062712 (2012)]. We estimate dispersion coefficients and collision rates between molecules and complexes. We find that the sticking-amplified three-body loss mechanism is not likely the cause of the losses observed in the experiments.

arXiv:1905.06691•Physical Review A•Sep 17, 2019•70 citations
A full-dimensional global potential energy surface (PES) for the reaction 23 Na 87 Rb + 23 Na 87 Rb → 23 Na 2 + 87 Rb 2 was constructed based on high-level ab initio calculations. The short-range part was expressed as a permutation invariant polynomial-neural network (PIP-NN) fit of 22 003 ab initio points calculated using a coupled cluster method with the one-electron basis 5 s 5 p 5 d 2 f plus effective core potentials and core polarizability potentials, while the long-range part was represented in an asymptotically correct form based on multipole expansion. The formation rate of the 23 Na 2 87 Rb 2 complex calculated using a quantum statistical method is in good agreement with experiment, while the estimated 19.20 μs lifetime of the complex from Rice–Ramsperger–Kassel–Marcus (RRKM) theory is significantly shorter than the measured millisecond decay rate, signaling either the inadequacy of RRKM theory or a yet unresolved loss mechanism.
A full-dimensional global potential energy surface (PES) for the KRb + KRb → K 2 Rb 2 * → K 2 + Rb 2 reaction is reported based on high-level ab initio calculations. The short-range part of the PES is fit with the permutationally invariant polynomial-neural network method, while the long-range parts of the PES in both the reactant and product asymptotes are represented by an asymptotically correct form. The long- and short-range parts are connected with intermediate-range parts to make them smooth. Within a statistical quantum model, this PES reproduces both the measured loss rates of ultracold KRb molecules and the K 2 and Rb 2 product state distributions, underscoring the important role of tunneling in ultracold chemistry. The PES also correctly predicts the lifetime of the K 2 Rb 2 * intermediate complex within the Rice–Ramsperger–Kassel–Marcus limit. It thus provides a reliable platform for future dynamical studies of the prototypical reaction.
The LiNa2 reactive system has recently received great attention in the experimental study of ultracold chemical reactions, but the corresponding theoretical calculations have not been carried out. Here, we report the first globally accurate ground-state LiNa2 potential energy surface (PES) using a Gaussian process model based on only 1776 actively selected high-level ab initio training points. The constructed PES had high precision and strong generalization capability. On the new PES, the quantum dynamics calculations on the Li(2S) + Na2(v = 0, j = 0) → LiNa + Na reaction were carried out in the 0.001–0.01 eV collision energy range using an improved time-dependent wave packet method. The calculated results indicate that this reaction is dominated by a complex-forming mechanism at low collision energies. The presented dynamics data provide guidance for experimental research, and the newly constructed PES could be further used for ultracold reaction dynamics calculations on this reactive system.
The ultracold reaction offers a unique opportunity to elucidate the intricate microscopic mechanism of chemical reactions, and the Na 2 Li system serves as a pivotal reaction system in the investigation of ultracold reactions. In this work, a high-precision potential energy surface (PES) of the Na 2 Li system is constructed based on high-level ab initio energy points and the neural network (NN) method, and a proper asymptotic functional form is adopted for the long-range interaction, which is suitable for the study of cold or ultracold collisions. Based on the new NN PES, the dynamics of the Li + Na 2 ( v = 0, j = 0) → Na + NaLi reaction are studied in the collision energy range of 10 –7 to 80 cm –1 . In the high collision energy range of 8 to 80 cm –1, the dynamics of the reaction is studied using the time-dependent wave packet method and the statistical quantum mechanical (SQM) method. Comparing the results of the two methods, it is found that the SQM method provides a rough description of the product ro-vibrational state distribution but overestimates the integral cross-section values. With the decrease of collision energy, the reaction differential cross section gradually changes from forward–backward symmetric scattering to predominant forward scattering. In the low collision energy range from 10 –7 to 8 cm –1, the SQM method is used to study the reaction dynamics, and the rate constant in the Wigner threshold region is estimated to be 2.87 × 10 –10 cm 3 /s.
Constructing accurate global potential energy surfaces (PESs) describing chemically reactive molecule-molecule collisions of alkali metal dimers presents a major challenge. To be suitable for quantum scattering calculations, such PESs must represent accurately three- and four-body interactions, describe conical intersections, and have a proper asymptotic form at the long range. Here, we demonstrate that such global potentials can be obtained by Gaussian Process (GP) regression merged with the analytic asymptotic expansions at the long range. We propose an efficient sampling technique, which allows us to construct an accurate global PES accounting for different chemical arrangements with <2500 ab initio calculations. We apply this method to (NaK)2 and obtain the first global PES for a system of four alkali metal atoms. The resulting surface exhibits a complex landscape including a pair and a quartet of symmetrically equivalent local minima and a seam of conical intersections. The dissociation energy found from our ab initio calculations is 4534 cm−1. This result is reproduced by the GP models with an error of less than 3%. The GP models of the PES allow us to analyze the features of the global PES, representative of general alkali metal four-atom interactions. Understanding these interactions is of key importance in the field of ultracold chemistry.
Abstract A fundamental question in the study of chemical reactions is how reactions proceed at a collision energy close to absolute zero. This question is no longer hypothetical: quantum degenerate gases of atoms and molecules can now be created at temperatures lower than a few tens of nanokelvin. Here we consider the benchmark ultracold reaction between, the most-celebrated ultracold molecule, KRb and K. We map out an accurate ab initio ground-state potential energy surface of the K 2 Rb complex in full dimensionality and report numerically-exact quantum-mechanical reaction dynamics. The distribution of rotationally resolved rates is shown to be Poissonian. An analysis of the hyperspherical adiabatic potential curves explains this statistical character revealing a chaotic distribution for the short-range collision complex that plays a key role in governing the reaction outcome.
We derive a quasiclassical expression for the density of states (DOS) of an arbitrary, ultracold, NNN-atom collision complex, for a general potential energy surface (PES). We establish the accuracy of our quasiclassical method by comparing to exact quantum results for the ${\mathrm{K}}_{2}\text{\ensuremath{-}}\mathrm{Rb}$ and NaK-NaK systems, with isotropic model PESs. Next, we calculate the DOS for an accurate NaK-NaK PES to be 0.124 $\ensuremath{\mu}{\mathrm{K}}^{\ensuremath{-}1}$, with an associated Rice-Ramsperger-Kassel-Marcus sticking time of 6.0 $\ensuremath{\mu}\mathrm{s}$. We extrapolate the DOS and sticking times to all other polar bialkali-bialkali collision complexes by scaling with atomic masses, equilibrium bond lengths, dissociation energies, and dispersion coefficients. The sticking times calculated here are two to three orders of magnitude shorter than those reported by Mayle et al. [Phys. Rev. A 85, 062712 (2012)]. We estimate dispersion coefficients and collision rates between molecules and complexes. We find that the sticking-amplified three-body loss mechanism is not likely the cause of the losses observed in the experiments.