Second-order correction to perfect pairing: An inexpensive electronic structure method for the treatment of strong electron-electron correlations

Low-Scaling, Efficient and Memory Optimized Computation of Nuclear Magnetic Resonance Shieldings within the Random Phase Approximation Using Cholesky-Decomposed Densities and an Attenuated Coulomb Metric

An efficient method for the computation of nuclear magnetic resonance (NMR) shielding tensors within the random phase approximation (RPA) is presented based on our recently introduced resolution-of-the-identity (RI) atomic orbital RPA NMR method [Drontschenko, V. J. Chem. Theory Comput. 2023, 19, 7542-7554] utilizing Cholesky decomposed density type matrices and employing an attenuated Coulomb RI metric. The introduced sparsity is efficiently exploited using sparse matrix algebra. This allows for an efficient and low-scaling computation of RPA NMR shielding tensors. Furthermore, we introduce a batching method for the computation of memory demanding intermediates that accounts for their sparsity. This extends the applicability of our method to even larger systems that would have been out of reach before, such as, e.g., a DNA strand with 260 atoms and 3408 atomic orbital basis functions.

Geminal-based electronic structure methods in quantum chemistry. Toward a geminal model chemistry

In this review, we discuss the recent progress in developing geminal-based theories for challenging problems in quantum chemistry. Specifically, we focus on the antisymmetrized geminal power, generalized valence bond, antisymmetrized product of strongly orthogonal geminals, singlet-type orthogonal geminals, the antisymmetric product of 1-reference orbital geminal, also known as the pair coupled cluster doubles ansatz, and geminals constructed from Richardson-Gaudin states. Furthermore, we review various corrections to account for the missing dynamical correlation effects in geminal models and possible extensions to target electronically excited states and open-shell species. Finally, we discuss some numerical examples and present-day challenges for geminal-based models, including a quantitative and qualitative analysis of wave functions, and software availability.

Cluster many-body expansion: A many-body expansion of the electron correlation energy about a cluster mean field reference

The many-body expansion (MBE) is an efficient tool that has a long history of use for calculating interaction energies, binding energies, lattice energies, and so on. In the past, applications of MBE to correlation energy have been unfeasible for large systems, but recent improvements to computing resources have sparked renewed interest in capturing the correlation energy using the generalized nth order Bethe-Goldstone equation. In this work, we extend this approach, originally proposed for a Slater determinant, to a tensor product state (TPS) based wavefunction. By partitioning the active space into smaller orbital clusters, our approach starts from a cluster mean field reference TPS configuration and includes the correlation contribution of the excited TPSs using the MBE. This method, named cluster MBE (cMBE), improves the convergence of MBE at lower orders compared to directly doing a block-based MBE from a RHF reference. We present numerical results for strongly correlated systems, such as the one- and two-dimensional Hubbard models and the chromium dimer. The performance of the cMBE method is also tested by partitioning the extended π space of several large π-conjugated systems, including a graphene nano-sheet with a very large active space of 114 electrons in 114 orbitals, which would require 10⁶⁶ determinants for the exact FCI solution.

Ring coupled cluster doubles at the multireference level

A ring approximation within an internally contracted multireference (MR) Coupled Cluster (CC) framework is worked out and tested. Derivation of equations utilizes MR based, generalized normal ordering and the corresponding generalized Wick-theorem (MR-GWT). Contractions among cluster operators are avoided by adopting a normal ordered exponential ansatz. The original version of the MR ring CC doubles (MR-rCCD) equations [Á. Szabados and Á. Margócsy, Mol. Phys. 115, 2731 (2017)] is rectified in two aspects. On the one hand, over-completeness of double excitations is treated by relying on the concept of frames. On the other hand, restriction on the maximal cumulant rank is lifted from two to four. This is found essential for obtaining reliable correlation corrections to the energy. The MR function underlying the approach is provided by the Generalized Valence Bond (GVB) model. The pair structure of the reference ensures a fragment structure of GVB cumulants. This represents a benefit when evaluating cumulant contractions appearing as a consequence of MR-GWT. In particular, cumulant involving terms remain less expensive than their traditional, pair-contracted counterpart, facilitating an O(N⁶) eventual scaling of the proposed MR-rCCD method. Pilot applications are presented for covalent bond breaking, deprotonation energies, and torsional potentials.

Locating multiple self-consistent field solutions: an approach inspired by metadynamics

We propose a method to locate the solutions to the self-consistent field (SCF) equations, using an approach based upon metadynamics. We define a distance function between density matrices. Within an SCF calculation, when a solution is found, a biasing potential based on distance from the solution is added to the energy to avoid reconvergence to the same solution. Multiple solutions can therefore be relatively easily found. Using this method we locate all known solutions and one unknown solution of the H4 model. The set of restricted Hartree-Fock (RHF) solutions for the nitrogen molecule is located, and a broken-symmetry solution lower in energy than the symmetric RHF solution is found corresponding to dissociation into doublet fragments.

An approximate diatomics in molecules formulation of generalized valence bond theory

The slow computational speed of the generalized valence bond perfect pairing method (GVB-PP) has been an impediment to its routine use. We have addressed this problem by employing a diatomics in molecules Hamiltonian derived from a second quantization perturbation approach. This results in all three- and four-centered two-electron integrals being dropped from the traditional GVB-PP calculation. For moderate sized molecules, as for example C20 computed with a double zeta + polarization basis, there is on average a fifty-fold decrease in computational times. In this article, we present the theory behind our approach and analyze the accuracy and speed of this approximate GVB-PP method for several cases where density functional methods have produced ambivalent results.

Geminal model chemistry. IV. Variational and size consistent pure spin states

We present a computationally inexpensive method that yields ground state wave functions of pure spin symmetry. The method is variational and rigorously size consistent, free from adjustable parameters, and has a favorable scaling with system size. It is based on the recently introduced partially spin restricted geminal wave functions with limited spin contamination. Computations of a bond breaking, a transition metal compound, and a symmetric hydrogen cluster confirm the properties of this method.

Geminal model chemistry III: Partial spin restriction

The authors define an ab initio electronic structure model that uses partial spin restriction. It is an intermediate case between the so-called spin-restricted and spin-unrestricted formulations, which are popular in electronic structure methodology. Partial spin restriction arises naturally when the wave function is represented as an antisymmetrized product of two-electron functions, as it is done in generalized valence bond and antisymmetrized product of strongly orthogonal geminal theories. The authors show that the new model is size consistent, and it improves the description of transition metal compounds.

Resolution of Identity Density Functional Theory Augmented with an Empirical Dispersion Term (RI-DFT-D): A Promising Tool for Studying Isolated Small Peptides

Resolution of identity standard density functional theory augmented with a damped empirical dispersion term (RI-DFT-D) calculations have been carried out on a set of lowest energy minima of tryptophyl-glycine (Trp-Gly) and tryptophyl-glycyl-glycine (Trp-Gly-Gly) peptides. RI-DFT-D (TPSS/TZVP) results are in excellent agreement with benchmark data based on the CCSD(T) method. Experimental spectra could be assigned according to the calculated IR frequencies. Central processing unit (CPU) time requirements are only slightly higher than those needed for the DFT calculations. Consequently, RI-DFT-D theory seems to be a promising methodology for studying oligopeptides with accuracy comparable to ab initio quantum chemical calculations.

Multireference correlation in long molecules with the quadratic scaling density matrix renormalization group

We have devised a local ab initio density matrix renormalization group algorithm to describe multireference correlations in large systems. For long molecules that are extended in one of their spatial dimensions, we can obtain an exact characterization of correlation, in the given basis, with a cost that scales only quadratically with the size of the system. The reduced scaling is achieved solely through integral screening and without the construction of correlation domains. We demonstrate the scaling, convergence, and robustness of the algorithm in polyenes and hydrogen chains. We converge to exact correlation energies (in the sense of full configuration interaction, with 1-10 microE(h) precision) in all cases and correlate up to 100 electrons in 100 active orbitals. We further use our algorithm to obtain exact energies for the metal-insulator transition in hydrogen chains and compare and contrast our results with those from conventional quantum chemical methods.

Hybrid correlation models based on active-space partitioning: Seeking accurate O(N5) ab initio methods for bond breaking

Moller-Plesset second-order (MP2) perturbation theory remains the least expensive standard ab initio method that includes electron correlation, scaling as O(N5) with the number of molecular orbitals N. Unfortunately, when restricted Hartree-Fock orbitals are employed, the potential energy curves calculated with this method are of little use at large interatomic separations because of the divergent behavior of MP2 in these regions. In our previous study [J. Chem. Phys. 122, 234110 (2005)] we combined the MP2 method with the singles and doubles coupled cluster (CCSD) method to produce a hybrid method that retains the computational scaling of MP2 and improves dramatically the shape of the MP2 curves. In this work we expand the hybrid methodology to several other schemes. We investigate a new, improved MP2-CCSD method as well as a few other O(N5) methods related to the Epstein-Nesbet pair correlation theory. Nonparallelity errors across the dissociation curve as well as several spectroscopic constants are computed for BH, HF, H2O, CH+, CH4, and Li2 molecules with the 6-31G* basis set and compared with the corresponding full configuration interaction results. We show that among the O(N5) methods considered, our new hybrid MP2-CCSD method is the most accurate and significantly outperforms MP2 not only at large interatomic separations, but also near equilibrium geometries.

On the Nature of Unrestricted Orbitals in Variational Active Space Wave Functions

Active space coupled cluster methods exhibit unusual, nonsmooth spin symmetry-breaking behavior where the unrestricted minimum lies higher in energy at short bond distances and crosses below the restricted solution at longer distances. The restricted solution is also observed to be a stable minimum slightly beyond the symmetry-breaking point. This behavior arises due to differences in the optimal active spaces defining the restricted and unrestricted wave functions and results in unrestricted wave functions that are not strictly size consistent. We suggest a new, size-consistent model that allows the orbitals to break spin symmetry only within the active space.

Valence bond theory for chemical dynamics

This essay provides a perspective on several issues in valence bond theory: the physical significance of semilocal bonding orbitals, the capability of valence bond concepts to explain systems with multireferences character, the use of valence bond theory to provide analytical representations of potential energy surfaces for chemical dynamics by the method of semiempirical valence bond potential energy surfaces (an early example of specific reaction parameters), by multiconfiguration molecular mechanics, by the combined valence bond-molecular mechanics method, and by the use of valence bond states as coupled diabatic states for describing electronically nonadiabatic processes (photochemistry). The essay includes both ab initio and semiempirical approaches.