Matrix elements of explicitly correlated Gaussian basis functions with arbitrary angular momentum

Explicitly correlated Gaussians for high-precision variational calculations of S^{e}, P^{e}, and D^{e} states of quantum systems: An efficient algorithm

In this work we consider an efficient algorithm for variational calculations of quantum few-particle systems in S, P, and D states of the even parity using all-particle explicitly correlated Gaussian (ECG) basis sets. We primarily focus on the description of states where the dominant configuration contains either two particles in p states or a single particle in a d state (all other particles are in s states). The basis functions we consider are products of spherically symmetric ECGs and bipolar harmonics. We introduced a scheme for deriving expressions for matrix elements of the overlap, kinetic and potential energy, as well as their derivatives with respect to the nonlinear parameters of the Gaussians. This allowed us to improve the efficiency of numerical calculations of the matrix elements (which is the most critical part of any code that uses ECGs) by one to two orders of magnitude compared to previous implementations. We provide a complete set of formulas for all basic matrix elements using the formalism of matrix differential calculus and discuss some technical details relevant to their efficient implementation. Lastly, we report a few example calculations of the ^{2}D^{e} states of the Li atom, ^{4}P^{e} and ^{2}D^{e} states of the B atom, and ^{3}P^{e} states of the C atom.

Deformed explicitly correlated Gaussians

Deformed explicitly correlated Gaussian (DECG) basis functions are introduced, and their matrix elements are calculated. All matrix elements can be calculated analytically in a closed form, except the Coulomb one, which has to be approximated by a Gaussian expansion. The DECG basis functions can be used to solve problems with nonspherical potentials. One example of such potential is the dipole self-interaction term in the Pauli-Fierz Hamiltonian. Examples are presented showing the accuracy and necessity of deformed Gaussian basis functions to accurately solve light-matter coupled systems in cavity QED.