Domain-Based Local Pair Natural Orbital Version of Mukherjee’s State-Specific Coupled Cluster Method

Frozen Natural Orbitals for the State-Averaged Driven Similarity Renormalization Group

We present a reduced-cost implementation of the state-averaged driven similarity renormalization group (SA-DSRG) based on the frozen natural orbital (FNO) approach. The natural orbitals (NOs) are obtained by diagonalizing the one-body reduced density matrix from SA-DSRG second-order perturbation theory (SA-DSRG-PT2). We consider three criteria to truncate the virtual NOs for the subsequent electron correlation treatment beyond SA-DSRG-PT2. An additive second-order correction is applied to the SA-DSRG Hamiltonian to reintroduce correlation effects from the discarded orbitals. The FNO SA-DSRG method is benchmarked on 35 small organic molecules in the QUEST database. When keeping 98-99% of the cumulative occupation numbers, the mean absolute error in the vertical transition energies due to FNO is less than 0.01 eV. Using the same FNO threshold, we observe a speedup of 9 times compared to the conventional SA-DSRG implementation for nickel carbonyl with a quadruple-ζ basis set. The FNO approach enables nonperturbative SA-DSRG computations on chloroiron corrole [FeCl(C19H11N4)] with more than 1000 basis functions, surpassing the current limit of a conventional implementation.

The Adsorption of Small Molecules on the Copper Paddle-Wheel: Influence of the Multi-Reference Ground State

We report a theoretical study of the adsorption of a set of small molecules (C₂H₂, CO, CO₂, O₂, H₂O, CH₃OH, C₂H₅OH) on the metal centers of the “copper paddle-wheel”—a key structural motif of many MOFs. A systematic comparison between DFT of different rungs, single-reference post-HF methods (MP2, SOS–MP2, MP3, DLPNO–CCSD(T)), and multi-reference approaches (CASSCF, DCD–CAS(2), NEVPT2) is performed in order to find a methodology that correctly describes the complicated electronic structure of paddle-wheel structure together with a reasonable description of non-covalent interactions. Apart from comparison with literature data (experimental values wherever possible), benchmark calculations with DLPNO–MR–CCSD were also performed. Despite tested methods show qualitative agreement in the majority of cases, we showed and discussed reasons for quantitative differences as well as more fundamental problems of specific cases.

An efficient Fock space multi-reference coupled cluster method based on natural orbitals: Theory, implementation, and benchmark

We present a natural orbital-based implementation of the intermediate Hamiltonian Fock space coupled-cluster method for the (1, 1) sector of Fock space. The use of natural orbitals significantly reduces the computational cost and can automatically choose an appropriate set of active orbitals. The new method retains the charge transfer separability of the original intermediate Hamiltonian Fock space coupled-cluster method and gives excellent performance for valence, Rydberg, and charge-transfer excited states. It offers significant computational advantages over the popular equation of motion coupled cluster method for excited states dominated by single excitations.

Electron Attachment to Cytosine: The Role of Water

We present an EOM-CCSD-based quantum mechanical/molecular mechanical (QM/MM) study on the electron attachment process to solvated cytosine. The electron attachment in the bulk solvated cytosine occurs through a doorway mechanism, where the initial electron is localized on water. The electron is subsequently transferred to cytosine by the mixing of electronic and nuclear degrees of freedom, which occurs on an ultrafast time scale. The bulk water environment stabilizes the cytosine-bound anion by an extensive hydrogen-bond network and drastically enhances the electron transfer rate from that observed in the gas phase. Microhydration studies cannot reproduce the effect of the bulk water environment on the electron attachment process, and one needs to include a large number of water molecules in the calculation to obtain converged results. The predicted adiabatic electron affinity and electron transfer rate obtained from our QM/MM calculations are consistent with the available experimental results.

Electronic and spin structures of CaMn4Ox clusters in the S0 state of the oxygen evolving complex of photosystem II. Domain-based local pair natural orbital (DLPNO) coupled-cluster (CC) calculations using optimized geometries and natural orbitals (UNO) by hybrid density functional theory (HDFT) calculations

Domain-based local pair natural orbital (DLPNO) coupled cluster single and double (CCSD) with triple perturbation (T) correction methods were performed to elucidate the relative stabilities of ten different intermediate structures of the CaMn4Ox cluster in the S0 state of the oxygen evolving complex (OEC) of photosystem II (PSII). Full geometry optimizations of all the S0 intermediates were performed by the UB3LYP-D3/Def2-TZVP methods, providing the assumed geometrical structures and starting natural orbitals (UNO) for DLPNO-CCSD(T)/Def2TZVP calculations. The effective exchange integrals (J) for the spin Hamiltonian models for the ten intermediates were obtained by the UB3LYP/Def2-TZVP calculations followed by the general spin projections. DLPNO-CCSD(T) calculations followed by the CBS extrapolation procedure elucidated that the (II, III, IV, IV) and (III, III, III, IV) valence states in the CaMn4O5 cluster of the OEC of the PS II were nearly degenerated in energy in the S0 state, indicating an important role of dynamical electron correlation effects for the valence and spin fluctuations in strongly correlated electron systems (SCESs) consisting of 3d transition metals.

Explicitly correlated coupled cluster method for accurate treatment of open-shell molecules with hundreds of atoms

We present a near-linear scaling formulation of the explicitly correlated coupled-cluster singles and doubles with the perturbative triples method [CCSD(T)_F12¯] for high-spin states of open-shell species. The approach is based on the conventional open-shell CCSD formalism [M. Saitow et al., J. Chem. Phys. 146, 164105 (2017)] utilizing the domain local pair-natural orbitals (DLPNO) framework. The use of spin-independent set of pair-natural orbitals ensures exact agreement with the closed-shell formalism reported previously, with only marginally impact on the cost (e.g., the open-shell formalism is only 1.5 times slower than the closed-shell counterpart for the C₁₆₀H₃₂₂ n-alkane, with the measured size complexity of ≈1.2). Evaluation of coupled-cluster energies near the complete-basis-set (CBS) limit for open-shell systems with more than 550 atoms and 5000 basis functions is feasible on a single multi-core computer in less than 3 days. The aug-cc-pVTZ DLPNO-CCSD(T)_F12¯ contribution to the heat of formation for the 50 largest molecules among the 348 core combustion species benchmark set [J. Klippenstein et al., J. Phys. Chem. A 121, 6580-6602 (2017)] had root-mean-square deviation (RMSD) from the extrapolated CBS CCSD(T) reference values of 0.3 kcal/mol. For a more challenging set of 50 reactions involving small closed- and open-shell molecules [G. Knizia et al., J. Chem. Phys. 130, 054104 (2009)], the aug-cc-pVQ(+d)Z DLPNO-CCSD(T)_F12¯ yielded a RMSD of ∼0.4 kcal/mol with respect to the CBS CCSD(T) estimate.

Near-Linear Scaling in DMRG-Based Tailored Coupled Clusters: An Implementation of DLPNO-TCCSD and DLPNO-TCCSD(T)

We present a new implementation of density matrix renormalization group based tailored coupled clusters method (TCCSD), which employs the domain-based local pair natural orbital approach (DLPNO). Compared to the previous local pair natural orbital (LPNO) version of the method, the new implementation is more accurate, offers more favorable scaling, and provides more consistent behavior across the variety of systems. On top of the singles and doubles, we include the perturbative triples correction (T), which is able to retrieve even more dynamic correlation. The methods were tested on three systems: tetramethyleneethane, oxo-Mn(Salen), and iron(II)-porphyrin model. The first two were revisited to assess the performance with respect to LPNO-TCCSD. For oxo-Mn(Salen), we retrieved between 99.8 and 99.9% of the total canonical correlation energy which is an improvement of 0.2% over the LPNO version in less than 63% of the total LPNO runtime. Similar results were obtained for iron(II)-porphyrin. When the perturbative triples correction was employed, irrespective of the active space size or system, the obtained energy differences between two spin states were within the chemical accuracy of 1 kcal/mol using the default DLPNO settings.

A multireference coupled-electron pair approximation combined with complete-active space perturbation theory in local pair-natural orbital framework

The Complete-Active Space Second-order Perturbation Theory (CASPT2) has been one of the most widely-used methods for reliably calculating electronic structures of multireference systems. Because of its lowest level treatment of dynamic correlation, it has a high computational feasibility; however, its accuracy in some cases falls short of needs. Here, as a simple yet higher-order alternative, we introduce a hybrid theory of the CASPT2 and a multireference variant of the Coupled-Electron Pair Approximation (CEPA), which is a class of high level correlation theory. A central feature of our theory (CEPT2) is to use the two underlying theories for describing different divisions of correlation components based on the full internal contraction framework. The external components, which usually give a major contribution to the dynamic correlation, are intensively described using the CEPA Ansatz, while the rests are treated at the CASPT2 level. Furthermore, to drastically reduce the computational demands, we have incorporated the pair-natural orbital (PNO) method into our multireference implementations. This development, thus, requires highly complex derivations and coding, while it has been largely facilitated with an automatic expression and code generation technique. To highlight the accuracy of the CEPT2 approach and to assess the errors caused by the PNO truncation, benchmark calculations are shown on small- to medium-size molecules, illustrating the high accuracy of the present CEPT2 model. By tightening the truncation thresholds, the PNO-CEPT2 energy converges toward the canonical counterpart and is more accurate than that of PNO-CASPT2 as long as the same truncation thresholds are used.

Linear scaling perturbative triples correction approximations for open-shell domain-based local pair natural orbital coupled cluster singles and doubles theory [DLPNO-CCSD(T0/T)]

The coupled cluster method with single-, double-, and perturbative triple excitations [CCSD(T)] is considered to be one of the most reliable quantum chemistry theories. However, the steep scaling of CCSD(T) has limited its application to small or medium-sized systems for a long time. In our previous work, the linear scaling domain based local pair natural orbital CCSD variant (DLPNO-CCSD) has been developed for closed-shell and open-shell. However, it is known from extensive benchmark studies that triple-excitation contributions are important to reach chemical accuracy. In the present work, two linear scaling (T) approximations for open-shell DLPNO-CCSD are implemented and compared: (a) an algorithm based on the semicanonical approximation, in which off-diagonal Fock matrix elements in the occupied space are neglected [referred to as DLPNO-(T₀)]; and (b) an improved algorithm in which the triples amplitudes are computed iteratively [referred to as DLPNO-(T)]. This work is based on the previous open-shell DLPNO-CCSD algorithm [M. Saitow et al., J. Chem. Phys. 146, 164105 (2017)] as well as the iterative (T) correction for closed-shell systems [Y. Guo et al., J. Chem. Phys. 148, 011101 (2018)]. Our results show that the new open-shell perturbative corrections, DLPNO-(T₀/T), can predict accurate absolute and relative correlation energies relative to the canonical reference calculations with the same basis set. The absolute energies from DLPNO-(T) are significantly more accurate than those of DLPNO-(T₀). The additional computational effort of DLPNO-(T) relative to DLPNO-(T₀) is a factor of 4 on average. We report calculations on systems with more than 4000 basis functions.

Toward the efficient local tailored coupled cluster approximation and the peculiar case of oxo-Mn(Salen)

We introduce a new implementation of the coupled cluster method with single and double excitations tailored by the matrix product state wave functions (DMRG-TCCSD), which employs the local pair natural orbital (LPNO) approach. By exploiting locality in the coupled cluster stage of the calculation, we were able to remove some of the limitations that hindered the application of the canonical version of the method to larger systems and/or with larger basis sets. We assessed the accuracy of the approximation using two systems: tetramethyleneethane (TME) and oxo-Mn(Salen). Using the default cut-off parameters, we were able to recover over 99.7% and 99.8% of the canonical correlation energy for the triplet and singlet state of TME, respectively. In the case of oxo-Mn(Salen), we found that the amount of retrieved canonical correlation energy depends on the size of the complete active space (CAS)-we retrieved over 99.6% for the larger 27 orbital CAS and over 99.8% for the smaller 22 orbital CAS. The use of LPNO-TCCSD allowed us to perform these calculations up to quadruple-ζ basis set, amounting to 1178 basis functions. Moreover, we examined dependence of the ground state of oxo-Mn(Salen) on the CAS composition. We found that the inclusion of 4d_xy orbital plays an important role in stabilizing the singlet state at the DMRG-CASSCF level via double-shell effect. However, by including dynamic correlation, the ground state was found to be triplet regardless of the size of the basis set or the composition of CAS, which is in agreement with previous findings by canonical DMRG-TCCSD in smaller basis.

Data-Driven Acceleration of the Coupled-Cluster Singles and Doubles Iterative Solver

Solving the coupled-cluster (CC) equations is a cost-prohibitive process that exhibits poor scaling with system size. These equations are solved by determining the set of amplitudes (t) that minimize the system energy with respect to the coupled-cluster equations at the selected level of truncation. Here, a novel approach to predict the converged coupled-cluster singles and doubles (CCSD) amplitudes, thus the coupled-cluster wave function, is explored by using machine learning and electronic structure properties inherent to the MP2 level. Features are collected from quantum chemical data, such as orbital energies, one-electron Hamiltonian, Coulomb, and exchange terms. The data-driven CCSD (DDCCSD) is not an alchemical method because the actual iterative coupled-cluster equations are solved. However, accurate energetics can also be obtained by bypassing solving the CC equations entirely. Our preliminary data show that it is possible to achieve remarkable speedups in solving the CCSD equations, especially when the correct physics are encoded and used for training of machine learning models.

Quantum information-based analysis of electron-deficient bonds

Recently, the correlation theory of the chemical bond was developed, which applies concepts of quantum information theory for the characterization of chemical bonds, based on the multiorbital correlations within the molecule. Here, for the first time, we extend the use of this mathematical toolbox for the description of electron-deficient bonds. We start by verifying the theory on the textbook example of a molecule with three-center two-electron bonds, namely, diborane(6). We then show that the correlation theory of the chemical bond is able to properly describe the bonding situation in more exotic molecules which have been synthesized and characterized only recently, in particular, the diborane molecule with four hydrogen atoms [diborane(4)] and a neutral zerovalent s-block beryllium complex, whose surprising stability was attributed to a strong three-center two-electron π bond stretching across the C-Be-C core. Our approach is of high importance especially in the light of a constant chase after novel compounds with extraordinary properties where the bonding is expected to be unusual.

Polynomial-product states: A symmetry-projection-based factorization of the full coupled cluster wavefunction in terms of polynomials of double excitations

Our goal is to remedy the failure of symmetry-adapted coupled-cluster theory in the presence of strong correlation. Previous work along these lines has taken us from a diagram-level analysis of the coupled-cluster equations to an understanding of the collective modes which can occur in various channels of the coupled-cluster equations to the exploration of non-exponential wavefunctions in efforts to combine coupled-cluster theory with symmetry projection. In this manuscript, we extend these efforts by introducing a new, polynomial product wavefunction ansatz that incorporates information from symmetry projection into standard coupled-cluster theory in a way that attempts to mitigate the effects of the lack of size extensivity and size consistency characteristic of symmetry-projected methods. We describe the new approach in detail within the context of our previous efforts, explore some illustrative calculations, and consider one route for reducing the computational cost of the new method.

Perturbative triples correction to domain-based local pair natural orbital variants of Mukherjee's state specific coupled cluster method

In this article we report an implementation of the perturbative triples correction to Mukherjee's state-specific multireference coupled cluster method based on the domain-based pair natural orbital approach (DLPNO-MkCC). We tested the performance of DLPNO-MkCCSD(T) in calculations involving tetramethyleneethane and isomers of naphthynes. These tests show that more than 97% of triples energy was recovered with respect to the canonical MkCCSD(T) method, which together with the DLPNO-MkCCSD part accounts for about 99.70-99.85% of the total correlation energy. The applicability of the method was demonstrated on calculations of singlet-triplet gaps for several large systems: triangulene, dynemicin A, and a beryllium complex.