TY - JOUR
T1 - Computer program ATOM-MOL-nonBO for performing calculations of ground and excited states of atoms and molecules without assuming the Born-Oppenheimer approximation using all-particle complex explicitly correlated Gaussian functions
AU - Bubin, Sergiy
AU - Adamowicz, Ludwik
PY - 2020/5/29
Y1 - 2020/5/29
N2 - In this work, we describe a computer program called ATOM-MOL-nonBO for performing bound state calculations of small atoms and molecules without assuming the Born-Oppenheimer approximation. All particles forming the systems, electrons and nuclei, are treated on equal footing. The wave functions of the bound states are expanded in terms of all-particle one-center complex explicitly correlated Gaussian functions multiplied by Cartesian angular factors. As these Gaussian functions are eigenfunctions of the operator representing the square of the total angular momentum of the system, the problem separates and calculations of states corresponding to different values of the total rotational quantum number can be solved independently from each other. Due to thorough variational optimization of the Gaussian exponential parameters, the method allows us to generate very accurate wave functions. The optimization is aided by analytically calculated energy gradient determined with respect to the parameters. Three examples of calculations performed for diatomic and triatomic molecules are shown as an illustration of calculations that can be performed with this program. Finally, we discuss the limitations, applicability range, and bottlenecks of the program.
AB - In this work, we describe a computer program called ATOM-MOL-nonBO for performing bound state calculations of small atoms and molecules without assuming the Born-Oppenheimer approximation. All particles forming the systems, electrons and nuclei, are treated on equal footing. The wave functions of the bound states are expanded in terms of all-particle one-center complex explicitly correlated Gaussian functions multiplied by Cartesian angular factors. As these Gaussian functions are eigenfunctions of the operator representing the square of the total angular momentum of the system, the problem separates and calculations of states corresponding to different values of the total rotational quantum number can be solved independently from each other. Due to thorough variational optimization of the Gaussian exponential parameters, the method allows us to generate very accurate wave functions. The optimization is aided by analytically calculated energy gradient determined with respect to the parameters. Three examples of calculations performed for diatomic and triatomic molecules are shown as an illustration of calculations that can be performed with this program. Finally, we discuss the limitations, applicability range, and bottlenecks of the program.
UR - http://www.scopus.com/inward/record.url?scp=85085909514&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85085909514&partnerID=8YFLogxK
U2 - 10.1063/1.5144268
DO - 10.1063/1.5144268
M3 - Article
C2 - 32486658
AN - SCOPUS:85085909514
VL - 152
SP - 204102
JO - Journal of Chemical Physics
JF - Journal of Chemical Physics
SN - 0021-9606
IS - 20
ER -