Affiliation:
1. Department of Physics, Stony Brook University, Stony Brook, NY 11794, USA
2. Institute for Advanced Computational Science, Stony Brook University, Stony Brook, NY 11794, USA
Abstract
We present Python Quasi-classical atom–molecule scattering (PyQCAMS v0.1.0), a new Python package for atom–diatom scattering within the quasi-classical trajectory approach. The input consists of the mass, collision energy, impact parameter, and pair-wise/three-body interactions. As the output, the code provides the vibrational quenching, dissociation, and reactive cross sections along with the rovibrational energy distribution of the reaction products. We benchmark the program for a reaction involving a molecular ion in a high-density ultracold gas, RbBa+ + Rb. Furthermore, we treat H2 + Ca → CaH + H reactions as a prototypical example to illustrate the properties and performance of the software. Finally, we study the parallelization performance of the code by looking into the speedup of the program as a function of the number of CPUs used.
Funder
the United States Air Force Office of Scientific Research
Reference38 articles.
1. Exchange Reactions with Activation Energy. I. Simple Barrier Potential for (H, H2);Karplus;J. Chem. Phys.,1965
2. Truhlar, D.G., and Muckerman, J.T. (1979). Atom-Molecule Collision Theory: A Guide for the Experimentalist, Plenum Press.
3. Ríos, J. (2020). An Introduction to Cold and Ultracold Chemistry: Atoms, Molecules, Ions and Rydbergs, Springer International Publishing.
4. Quasiclassical state to state reaction cross sections for D+H2(v = 0, j = 0)→HD(v′,j′) + H. Formation and characteristics of short-lived collision complexes;Aoiz;J. Chem. Phys.,1992
5. Recent results from quasiclassical trajectory computations of elementary chemical reactions;Aoiz;J. Chem. Soc. Faraday Trans.,1998