B AOInts

B.3 AOInts: Calculating ERIs with Q-Chem

The area of molecular integrals with respect to Gaussian basis functions has recently been reviewed300 and the user is referred to this review for deeper discussions and further references to the general area. The purpose of this short account is to present the basic approach, and in particular, the implementation of ERI algorithms and aspects of interest to the user in the AOInts package which underlies the Q-Chem program.

We begin by observing that all of the integrals encountered in an ab initio calculation, of which overlap, kinetic energy, multipole moment, internuclear repulsion, nuclear-electron attraction and inter electron repulsion are the best known, can be written in the general form

(𝐚𝐛|𝐜𝐝)=ϕ𝐚(𝐫1)ϕ𝐛(𝐫1)θ(r12)ϕ𝐜(𝐫2)ϕ𝐝(𝐫2)𝑑𝐫1𝑑𝐫2 (B.1)

where the basis functions are contracted Gaussians (CGTF)

ϕ𝐚(𝐫)=(x-Ax)ax(y-Ay)ay(z-Az)azi=1KaDaie-αi|𝐫-𝐀|2 (B.2)

and the operator θ is a two-electron operator. Of the two-electron operators (Coulomb, CASE, anti-Coulomb and delta-function) used in the Q-Chem program, the most significant is the Coulomb, which leads us to the ERIs.

An ERI is the classical Coulomb interaction, θ(x)=1/x in Eq. (B.1), between two charge distributions referred to as bras (𝐚𝐛| and kets |𝐜𝐝).