A new functional form of the explicitly correlated Gaussian-type functions (later called Gaussians or ECGs) for performing non-Born-Oppenheimer (BO) calculations of molecular systems with an arbitrary number of nuclei is presented. In these functions, the exponential part explicitly depends on all interparticle distances and the preexponential part depends only on the distances between the nuclei.
The new Gaussians are called sin/cos-Gaussians and their preexponential part is a product of sin and/or cos factors. The effectiveness of the new Gaussians in describing non-BO pure vibrational states is investigated by comparing them with r(m)-Gaussians containing preexponential multipliers in the form of non-negative powers of internuclear distances (the internuclear distance in the diatomic case).
The testing is performed for a diatomic system with the nuclei interacting through a Morse potential. It shows that the new sin/cos-Gaussian basis set is capable of providing equally accurate results as obtained with the r(m)-Gaussians.
However, especially for lower vibrational states, more sin/cos-Gaussians are needed to reach a similar accuracy level as obtained with the r(m)-Gaussians. Implementation of the sin/cos-Gaussians in non-BO calculations of diatomic and, in particular, of triatomic systems, which will follow, will provide further assessment of the efficiency of the new functions.