Non-adiabatic quantum reactive scattering in hyperspherical coordinates
A new electronically non-adiabatic quantum reactive scattering methodology is presented based on a time-independent coupled channel formalism and the adiabatically adjusting principal axis hyperspherical coordinates of Pack and Parker [J. Chem. Phys. 87, 3888 (1987)]. The methodology computes the full state-to-state scattering matrix for A + B _{2}(v, j) ↔ AB(v', j') + B and A + AB(v, j) → A + AB(v', j') reactions that involve two coupled electronic states which exhibit a conical intersection. The methodology accurately treats all six degrees of freedom relative to the center-of-mass which includes non-zero total angular momentum J and identical particle exchange symmetry. The new methodology is applied to the ultracold hydrogen exchange reaction for which large geometric phase effects have been recently reported [B. K. Kendrick et al., Phys. Rev. Lett. 115, 153201 (2015)]. Rate coefficients for the H/D + HD(v = 4, j = 0) → H/D + HD(v', j') reactions are reported for collision energies between 1 μK and 100 K (total energy ≈1.9 eV). A new diabatic potential energy matrix is developed based on the Boothroyd, Keogh, Martin, and Peterson (BKMP2) and double many body expansion plus single-polynomial (DSP) adiabatic potential energy surfaces for the ground and firstmore »
- Publication Date:
- Report Number(s):
- LA-UR-17-30356
Journal ID: ISSN 0021-9606
- Grant/Contract Number:
- AC52-06NA25396
- Type:
- Accepted Manuscript
- Journal Name:
- Journal of Chemical Physics
- Additional Journal Information:
- Journal Volume: 148; Journal Issue: 4; Journal ID: ISSN 0021-9606
- Publisher:
- American Institute of Physics (AIP)
- Research Org:
- Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
- Sponsoring Org:
- USDOE Laboratory Directed Research and Development (LDRD) Program; USDOE National Nuclear Security Administration (NNSA)
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; Inorganic and Physical Chemistry; non-adiabatic dynamics; ultracold chemistry; quantum reactive scattering
- OSTI Identifier:
- 1431067
- Alternate Identifier(s):
- OSTI ID: 1418724
Kendrick, Brian K. Non-adiabatic quantum reactive scattering in hyperspherical coordinates. United States: N. p.,
Web. doi:10.1063/1.5014989.
Kendrick, Brian K. Non-adiabatic quantum reactive scattering in hyperspherical coordinates. United States. doi:10.1063/1.5014989.
Kendrick, Brian K. 2018.
"Non-adiabatic quantum reactive scattering in hyperspherical coordinates". United States.
doi:10.1063/1.5014989. https://www.osti.gov/servlets/purl/1431067.
@article{osti_1431067,
title = {Non-adiabatic quantum reactive scattering in hyperspherical coordinates},
author = {Kendrick, Brian K.},
abstractNote = {A new electronically non-adiabatic quantum reactive scattering methodology is presented based on a time-independent coupled channel formalism and the adiabatically adjusting principal axis hyperspherical coordinates of Pack and Parker [J. Chem. Phys. 87, 3888 (1987)]. The methodology computes the full state-to-state scattering matrix for A + B2(v, j) ↔ AB(v', j') + B and A + AB(v, j) → A + AB(v', j') reactions that involve two coupled electronic states which exhibit a conical intersection. The methodology accurately treats all six degrees of freedom relative to the center-of-mass which includes non-zero total angular momentum J and identical particle exchange symmetry. The new methodology is applied to the ultracold hydrogen exchange reaction for which large geometric phase effects have been recently reported [B. K. Kendrick et al., Phys. Rev. Lett. 115, 153201 (2015)]. Rate coefficients for the H/D + HD(v = 4, j = 0) → H/D + HD(v', j') reactions are reported for collision energies between 1 μK and 100 K (total energy ≈1.9 eV). A new diabatic potential energy matrix is developed based on the Boothroyd, Keogh, Martin, and Peterson (BKMP2) and double many body expansion plus single-polynomial (DSP) adiabatic potential energy surfaces for the ground and first excited electronic states of H3, respectively. The rate coefficients computed using the new non-adiabatic methodology and diabatic potential matrix reproduce the recently reported rates that include the geometric phase and are computed using a single adiabatic ground electronic state potential energy surface (BKMP2). The dramatic enhancement and suppression of the ultracold rates due to the geometric phase are confirmed as well as its effects on several shape resonances near 1 K. In conclusion, the results reported here represent the first fully non-adiabatic quantum reactive scattering calculation for an ultracold reaction and validate the importance of the geometric phase on the Wigner threshold behavior.},
doi = {10.1063/1.5014989},
journal = {Journal of Chemical Physics},
number = 4,
volume = 148,
place = {United States},
year = {2018},
month = {1}
}