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

New Journal of Physics

Quantum dynamics of tunneling dominated reactions at low temperatures

Jisha Hazra, N. Balakrishnan

Published May 1, 2015• 17 citations
Publisher

Abstract

Abstract We report a quantum dynamics study of the Li + HF → LiF + H reaction at low temperatures of interest to cooling and trapping experiments. Contributions from non-zero partial waves are analyzed and results show narrow resonances in the energy dependence of the cross section that survive partial wave summation. The computations are performed using the ABC code and a simple modification of the ABC code that enables separate energy cutoffs for the reactant and product rovibrational energy levels is found to dramatically reduce the basis set size and computational expense. Results obtained using two ab initio electronic potential energy surfaces for the LiHF system show strong sensitivity to the choice of the potential. In particular, small differences in the barrier heights of the two potential surfaces are found to dramatically influence the reaction cross sections at low energies. Comparison with recent measurements of the reaction cross section (Bobbenkamp et al 2011 J. Chem. Phys. 135 204306 ) shows similar energy dependence in the threshold regime and an overall good agreement with experimental data compared to previous theoretical results. Also, usefulness of a recently introduced method for ultracold reactions that employ the quantum close-coupling method at short-range and the multichannel quantum defect theory at long-range, is demonstrated in accurately evaluating product state-resolved cross sections for D + H2 and H + D2 reactions.

DOI
10.1088/1367-2630/17/5/055027

Related papers

Quantum reactive scattering calculations for the cold and ultracold Li + LiNa → Li2 + Na reaction

Brian K. Kendrick

A first-principles based quantum dynamics study of the Li + LiNa(v = 0, j = 0) → Li2(v′, j′) + Na reaction is reported for collision energies spanning the ultracold (1 nK) to cold (1 K) regimes. A full-dimensional ab initio potential energy surface for the ground electronic state of Li2Na is utilized that includes an accurate treatment of the long-range interactions. The Li + LiNa reaction is barrierless and exoergic and exhibits a deep attractive potential well that supports complex formation. Thus, significant reactivity occurs even for collision temperatures approaching absolute zero. The reactive scattering calculations are based on a numerically exact time-independent quantum dynamics methodology in hyperspherical coordinates. Total and rotationally resolved rate coefficients are reported at 56 collision energies and include all contributing partial waves. Several shape resonances are observed in many of the rotationally resolved rate coefficients and a small resonance feature is also reported in the total rate coefficient near 50 mK. Of particular interest, the angular distributions or differential cross sections are reported as a function of both the collision energy and scattering angle. Unique quantum fingerprints (bumps, channels, and ripples) are observed in the angular distributions for each product rotational state due to quantum interference and shape resonance contributions. The Li + LiNa reaction is under active experimental investigation so that these intriguing features could be verified experimentally when sufficient product state resolution becomes feasible for collision energies below 1 K.

The Journal of Chemical Physics•Mar 22, 2021•15 citations

Adiabatic channel capture theory applied to cold atom–molecule reactions: Li + CaH $\to $ LiH + Ca at 1K

Timur V. Tscherbul, Alexei A Buchachenko

We use quantum and classical adiabatic capture theories to study the chemical reaction Li + CaH LiH + Ca. Using a recently developed ab initio potential energy surface, which provides an accurate representation of long-range interactions in the entrance reaction channel, we calculate the adiabatic channel potentials by diagonalizing the Li-CaH Hamiltonian as a function of the atom-molecule separation. The resulting adiabatic channel potentials are used to calculate both the classical and quantum capture probabilities as a function of collision energy, as well as the temperature dependencies of the partial and total reaction rates. The calculated reaction rate agrees well with the measured value at 1 K (V Singh et al 2012 Phys. Rev. Lett. 108 203201), suggesting that the title reaction proceeds without an activation barrier. The calculated classical adiabatic capture rate agrees well with the quantum result in the multiple-partial-wave regime of relevance to the experiment. Significant differences are found only in the ultracold limit ( mK), demonstrating that adiabatic capture theories can predict the reaction rates with nearly quantitative accuracy in the multiple-partial-wave regime.

arXiv:1412.4575•New Journal of Physics•Mar 20, 2015•16 citations

Low-energy limit for tunnelling subject to an Eckart potential barrier

Samuel H. Lipoff, Dudley R. Herschbach

For two-body s-wave collisions subject to tunnelling though an Eckart potential barrier, cross-sections and rate coefficients in the low-energy regime can be evaluated in analytic form. These provide criteria for approach to the Wigner limit and a generic plot applicable to ultracold collisions (<1 mK). The Eckart barrier shape is found to fit fairly well accurate potential curves available for F(2P) + F(2P) and He(1S) + He*(1S). Also discussed are tunnelling-dominated ultracold chemical reactions, exemplified by F + H2. The very slow approach of the reactants allows averaging over rovibrational motions, so the effective adiabatic (‘dressed’) potential surface differs substantially from the Born–Oppenheimer (‘bare’) surface obtained from electronic structure calculations. The dressed barrier is more than two-fold lower and also thinner. For it the Eckart model gives an ultracold rate coefficient for forming HF in its v′ = 2, j′ = 0 rovibrational state (6.5 × 10−14 cm3/s) that is about 200-fold higher than for the bare barrier and of the same magnitude found from a full-scale 3-D quantum scattering calculation.

Molecular Physics•Apr 10, 2010•13 citations

Non-adiabatic quantum dynamics of the ultracold Li+LiNa→ Li2+Na chemical reaction

Brian K. Kendrick, Mingzhe Li, H. Li +3 more

Synopsis We report non-adiabatic dynamics of the Li+LiNa→Li2+Na chemical reaction at cold and ultracold temperatures employing accurate ab initio electronic potential energy surfaces in a quantum dynamics formulation employing a diabatic representation. Results are compared against those from a single adiabatic ground state potential energy surface and a universal model based on the long-range interaction potential. We discuss signatures of non-universal behavior in the total rate coefficients as well as strong non-adiabatic effects in the state-to-state rotationally resolved rate coefficients.

Journal of Physics Conference Series•Jan 1, 2020•2 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

Interactions and dynamics in Li+Li2 ultracold collisions

Marko T. Cvitaš, Pavel Soldán, Jeremy M. Hutson +2 more

A potential energy surface for the lowest quartet electronic state ((4)A(')) of lithium trimer is developed and used to study spin-polarized Li+Li(2) collisions at ultralow kinetic energies. The potential energy surface allows barrierless atom exchange reactions. Elastic and inelastic cross sections are calculated for collisions involving a variety of rovibrational states of Li(2). Inelastic collisions are responsible for trap loss in molecule production experiments. Isotope effects and the sensitivity of the results to details of the potential energy surface are investigated. It is found that for vibrationally excited states, the cross sections are only quite weakly dependent on details of the potential energy surface.

arXiv:physics/0703136•The Journal of Chemical Physics•Aug 20, 2007•71 citations
Abstract We report a quantum dynamics study of the Li + HF → LiF + H reaction at low temperatures of interest to cooling and trapping experiments. Contributions from non-zero partial waves are analyzed and results show narrow resonances in the energy dependence of the cross section that survive partial wave summation. The computations are performed using the ABC code and a simple modification of the ABC code that enables separate energy cutoffs for the reactant and product rovibrational energy levels is found to dramatically reduce the basis set size and computational expense. Results obtained using two ab initio electronic potential energy surfaces for the LiHF system show strong sensitivity to the choice of the potential. In particular, small differences in the barrier heights of the two potential surfaces are found to dramatically influence the reaction cross sections at low energies. Comparison with recent measurements of the reaction cross section (Bobbenkamp et al 2011 J. Chem. Phys. 135 204306 ) shows similar energy dependence in the threshold regime and an overall good agreement with experimental data compared to previous theoretical results. Also, usefulness of a recently introduced method for ultracold reactions that employ the quantum close-coupling method at short-range and the multichannel quantum defect theory at long-range, is demonstrated in accurately evaluating product state-resolved cross sections for D + H2 and H + D2 reactions.
A first-principles based quantum dynamics study of the Li + LiNa(v = 0, j = 0) → Li2(v′, j′) + Na reaction is reported for collision energies spanning the ultracold (1 nK) to cold (1 K) regimes. A full-dimensional ab initio potential energy surface for the ground electronic state of Li2Na is utilized that includes an accurate treatment of the long-range interactions. The Li + LiNa reaction is barrierless and exoergic and exhibits a deep attractive potential well that supports complex formation. Thus, significant reactivity occurs even for collision temperatures approaching absolute zero. The reactive scattering calculations are based on a numerically exact time-independent quantum dynamics methodology in hyperspherical coordinates. Total and rotationally resolved rate coefficients are reported at 56 collision energies and include all contributing partial waves. Several shape resonances are observed in many of the rotationally resolved rate coefficients and a small resonance feature is also reported in the total rate coefficient near 50 mK. Of particular interest, the angular distributions or differential cross sections are reported as a function of both the collision energy and scattering angle. Unique quantum fingerprints (bumps, channels, and ripples) are observed in the angular distributions for each product rotational state due to quantum interference and shape resonance contributions. The Li + LiNa reaction is under active experimental investigation so that these intriguing features could be verified experimentally when sufficient product state resolution becomes feasible for collision energies below 1 K.
We use quantum and classical adiabatic capture theories to study the chemical reaction Li + CaH LiH + Ca. Using a recently developed ab initio potential energy surface, which provides an accurate representation of long-range interactions in the entrance reaction channel, we calculate the adiabatic channel potentials by diagonalizing the Li-CaH Hamiltonian as a function of the atom-molecule separation. The resulting adiabatic channel potentials are used to calculate both the classical and quantum capture probabilities as a function of collision energy, as well as the temperature dependencies of the partial and total reaction rates. The calculated reaction rate agrees well with the measured value at 1 K (V Singh et al 2012 Phys. Rev. Lett. 108 203201), suggesting that the title reaction proceeds without an activation barrier. The calculated classical adiabatic capture rate agrees well with the quantum result in the multiple-partial-wave regime of relevance to the experiment. Significant differences are found only in the ultracold limit ( mK), demonstrating that adiabatic capture theories can predict the reaction rates with nearly quantitative accuracy in the multiple-partial-wave regime.
For two-body s-wave collisions subject to tunnelling though an Eckart potential barrier, cross-sections and rate coefficients in the low-energy regime can be evaluated in analytic form. These provide criteria for approach to the Wigner limit and a generic plot applicable to ultracold collisions (<1 mK). The Eckart barrier shape is found to fit fairly well accurate potential curves available for F(2P) + F(2P) and He(1S) + He*(1S). Also discussed are tunnelling-dominated ultracold chemical reactions, exemplified by F + H2. The very slow approach of the reactants allows averaging over rovibrational motions, so the effective adiabatic (‘dressed’) potential surface differs substantially from the Born–Oppenheimer (‘bare’) surface obtained from electronic structure calculations. The dressed barrier is more than two-fold lower and also thinner. For it the Eckart model gives an ultracold rate coefficient for forming HF in its v′ = 2, j′ = 0 rovibrational state (6.5 × 10−14 cm3/s) that is about 200-fold higher than for the bare barrier and of the same magnitude found from a full-scale 3-D quantum scattering calculation.
Synopsis We report non-adiabatic dynamics of the Li+LiNa→Li2+Na chemical reaction at cold and ultracold temperatures employing accurate ab initio electronic potential energy surfaces in a quantum dynamics formulation employing a diabatic representation. Results are compared against those from a single adiabatic ground state potential energy surface and a universal model based on the long-range interaction potential. We discuss signatures of non-universal behavior in the total rate coefficients as well as strong non-adiabatic effects in the state-to-state rotationally resolved rate coefficients.
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.
A potential energy surface for the lowest quartet electronic state ((4)A(')) of lithium trimer is developed and used to study spin-polarized Li+Li(2) collisions at ultralow kinetic energies. The potential energy surface allows barrierless atom exchange reactions. Elastic and inelastic cross sections are calculated for collisions involving a variety of rovibrational states of Li(2). Inelastic collisions are responsible for trap loss in molecule production experiments. Isotope effects and the sensitivity of the results to details of the potential energy surface are investigated. It is found that for vibrationally excited states, the cross sections are only quite weakly dependent on details of the potential energy surface.