Perturbation theory of the electron correlation cusp based on a partitioning of the electron–electron interaction into long- and short-range components

Resolutions of the Coulomb operator. VI. Computation of auxiliary integrals

We discuss the efficient computation of the auxiliary integrals that arise when resolutions of two-electron operators (specifically, the Coulomb operator [T. Limpanuparb, A. T. B. Gilbert, and P. M. W. Gill, J. Chem. Theory Comput. 7, 830 (2011)] and the long-range Ewald operator [T. Limpanuparb and P. M. W. Gill, J. Chem. Theory Comput. 7, 2353 (2011)]) are employed in quantum chemical calculations. We derive a recurrence relation that facilitates the generation of auxiliary integrals for Gaussian basis functions of arbitrary angular momentum and propose a near-optimal algorithm for its use.

Unified treatment of chemical and van der Waals forces via symmetry-adapted perturbation expansion

We propose a symmetry-adapted perturbation theory (SAPT) expansion of the intermolecular interaction energy which in a finite order provides the correct values of the constants determining the asymptotics of the interaction energy (the van der Waals constants) and is convergent when the energy of the interacting system is submerged in the continuum of Pauli-forbidden states-the situation common when at least one of the monomers has more than two electrons. These desirable features are achieved by splitting the intermolecular electron-nucleus attraction terms of the Hamiltonian into regular (long-range) and singular (short-range) parts. In the perturbation theory development, the regular part is treated as in the conventional polarization theory, which guarantees the correct asymptotics of the interaction energy, while the singular part is weakened sufficiently by an application of permutational symmetry projectors so that a convergent perturbation series is obtained. The convergence is demonstrated numerically, for both the chemical and van der Waals minima, by performing high-order calculations of the interaction energy of the ground-state lithium and hydrogen atoms-the simplest system for which the physical ground state is submerged in the Pauli-forbidden continuum. The obtained expansion enables a systematic extension of SAPT calculations beyond second order with respect to the intermolecular interaction operator.

Convergence Behavior of Symmetry-Adapted Perturbation Expansions for Excited States. A Model Study of Interactions Involving a Triplet Helium Atom

The convergence behavior of symmetry-adapted perturbation theory (SAPT) expansions is investigated for an interacting system involving an excited, open-shell monomer. By performing large-order numerical calculations for the interaction of the lowest, 1s2s, triplet state of helium with the ground state of the hydrogen atom we show that the conventional polarization and symmetrized Rayleigh-Schrödinger expansions diverge in this case. This divergence is attributed to the continuum of intruder states appearing when the hydrogen electron is falling on the helium 1s orbital and the 2s electron is ejected from the interacting system. One of the dimer states resulting from the interaction becomes then a resonance, which presents a hard case to treat by a perturbation theory. We show that the SAPT expansions employing the strong symmetry-forcing procedure, such as the Eisenschitz-London- Hirschfelder-van der Avoird or the Amos-Musher theories, can cope with this situation and lead to convergent series when the permutational symmetry of the bound, quartet state is forced. However, these theories suffer from a wrong asymptotic behavior of the second- and higher-order energies when the interatomic distance R grows to infinity, which makes them unsuitable for practical applications. We show that by a suitable regularization of the Coulomb potential and by treating differently the regular, long-range and the singular, short-range parts of the interatomic electron-nucleus attraction terms in the Hamiltonian one obtains a perturbation expansion which has the correct asymptotic behavior in each order and which converges fast for a wide range of interatomic distances.

Correlation energy extrapolation by intrinsic scaling. I. Method and application to the neon atom

Remarkably accurate scaling relations are shown to exist between the correlation energy contributions from various excitation levels of the configuration interaction approach, considered as functions of the size of the correlating orbital space. These relationships are used to develop a method for extrapolating a sequence of smaller configuration interaction calculations to the full configuration-interaction energy. Calculations of the neon atom ground state with the Dunning's quadruple zeta basis set demonstrate the ability of the method to obtain benchmark quality results.