ABSTRACT
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.
