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

Second quantization-based symmetry-adapted perturbation theory: Generalizing exchange beyond single electron pair approximation

This paper presents a general second-quantized form of a permutation operator interchanging n pairs of electrons between interacting subsystems in the framework of the symmetry-adapted perturbation theory (SAPT). We detail the procedure for constructing this operator through the consecutive multiplication of single-pair permutation operators. This generalized form of the permutation operator has enabled the derivation of universal formulas for S2n approximations of the exchange energies in the first and second order of the interaction operator. We present expressions for corrections of S4 approximations and assess its efficacy on a selection of systems anticipated to exhibit a slowly converging overlap expansion. Additionally, we outline a method to sum the overlap expansion series to infinity in second-quantization, up to the second order in V. This new approach offers an alternative to the existing formalism based on density-matrix formulations. When combined with a symbolic algebra program for automated derivations, it paves the way for advancements in SAPT theory, particularly for intricate wavefunction theories.

Symmetry-Adapted Perturbation Theory Based on Multiconfigurational Wave Function Description of Monomers

We present a formulation of the multiconfigurational (MC) wave function symmetry-adapted perturbation theory (SAPT). The method is applicable to noncovalent interactions between monomers which require a multiconfigurational description, in particular when the interacting system is strongly correlated or in an electronically excited state. SAPT(MC) is based on one- and two-particle reduced density matrices of the monomers and assumes the single-exchange approximation for the exchange energy contributions. Second-order terms are expressed through response properties from extended random phase approximation (ERPA). The dispersion components of SAPT(MC) have been introduced in our previous works [HapkaM.J. Chem. Theory Comput.2019, 15, 1016−102730525591; HapkaM.J. Chem. Theory Comput.2019, 15, 6712–672331670950]. SAPT(MC) is applied either with generalized valence bond perfect pairing (GVB) or with complete active space self-consistent field (CASSCF) treatment of the monomers. We discuss two model multireference systems: the H₂ ··· H₂ dimer in out-of-equilibrium geometries and interaction between the argon atom and excited state of ethylene. Using the C₂H₄* ··· Ar complex as an example, we examine second-order terms arising from negative transitions in the linear response function of an excited monomer. We demonstrate that the negative-transition terms must be accounted for to ensure qualitative prediction of induction and dispersion energies and develop a procedure allowing for their computation. Factors limiting the accuracy of SAPT(MC) are discussed in comparison with other second-order SAPT schemes on a data set of small single-reference dimers.

Second-Order Dispersion Energy Based on Multireference Description of Monomers

We propose a method for calculating a second-order dispersion energy for weakly interacting multireference systems in arbitrary electronic states. It is based on response properties obtained from extended random phase approximation equations. The introduced formalism is general and requires only one- and two-particle reduced density matrices of monomers. We combine the new method with either generalized valence bond perfect pairing (GVB) or complete active space (CAS) self-consistent field description of the interacting systems. In addition to a general scheme, three approximations, leading to significant reduction of the computational cost, are developed by exploiting Dyall partitioning of the monomer Hamiltonians. For model multireference systems (H₂···H₂ and Be···Be) the method is accurate, unlike its single-reference-based counterpart. Neither GVB nor CAS description of single-reference monomers improves the dispersion energy with respect to the Hartree-Fock-based results.

First-order symmetry-adapted perturbation theory for multiplet splittings

We present a symmetry-adapted perturbation theory (SAPT) for the interaction of two high-spin open-shell molecules (described by their restricted open-shell Hartree-Fock determinants) resulting in low-spin states of the complex. The previously available SAPT formalisms, except for some system-specific studies for few-electron complexes, were restricted to the high-spin state of the interacting system. Thus, the new approach provides, for the first time, a SAPT-based estimate of the splittings between different spin states of the complex. We have derived and implemented the lowest-order SAPT term responsible for these splittings, that is, the first-order exchange energy. We show that within the so-called S² approximation commonly used in SAPT (neglecting effects that vanish as fourth or higher powers of intermolecular overlap integrals), the first-order exchange energies for all multiplets are linear combinations of two matrix elements: a diagonal exchange term that determines the spin-averaged effect and a spin-flip term responsible for the splittings between the states. The numerical factors in this linear combination are determined solely by the Clebsch-Gordan coefficients: accordingly, the S² approximation implies a Heisenberg Hamiltonian picture with a single coupling strength parameter determining all the splittings. The new approach is cast into both molecular-orbital and atomic-orbital expressions: the latter enable an efficient density-fitted implementation. We test the newly developed formalism on several open-shell complexes ranging from diatomic systems (Li⋯H, Mn⋯Mn, …) to the phenalenyl dimer.

Symmetry-adapted perturbation theory utilizing density functional description of monomers for high-spin open-shell complexes

We present an implementation of symmetry-adapted perturbation theory (SAPT) to interactions of high-spin open-shell monomers forming high-spin dimers. The monomer spin-orbitals used in the expressions for the electrostatic and exchange contributions to the interaction energy are obtained from density functional theory using a spin-restricted formulation of the open-shell Kohn-Sham (ROKS) method. The dispersion and induction energies are expressed through the density-density response functions predicted by the time-dependent ROKS theory. The method was applied to several systems: NH...He, CN...Ne, H2O...HO2, and NH...NH. It provides accuracy comparable to that of the best previously available methods such as the open-shell coupled-cluster method with single, double, and noniterative triple excitations, RCCSD(T), with a significantly reduced computational cost.

Third-order interactions in symmetry-adapted perturbation theory

We present an extension of many-body symmetry-adapted perturbation theory (SAPT) by including all third-order polarization and exchange contributions obtained with the neglect of intramonomer correlation effects. The third-order polarization energy, which naturally decomposes into the induction, dispersion, and mixed, induction-dispersion components, is significantly quenched at short range by electron exchange effects. We propose a decomposition of the total third-order exchange energy into the exchange-induction, exchange-dispersion, and exchange-induction-dispersion contributions which provide the quenching for the corresponding individual polarization contributions. All components of the third-order energy have been expressed in terms of molecular integrals and orbital energies. The obtained formulas, valid for both dimer- and monomer-centered basis sets, have been implemented within the general closed-shell many-electron SAPT program. Test calculations for several small dimers have been performed and their results are presented. For dispersion-bound dimers, the inclusion of the third-order effects eliminates the need for a hybrid SAPT approach, involving supermolecular Hartree-Fock calculations. For dimers consisting of strongly polar monomers, the hybrid approach remains more accurate. It is shown that, due to the extent of the quenching, the third-order polarization effects should be included only together with their exchange counterparts. Furthermore, the latter have to be calculated exactly, rather than estimated by scaling the second-order values.

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.