Communication: Multipole approximations of distant pair energies in local correlation methods with pair natural orbitals

Cluster perturbation theory IX: Perturbation series for the coupled cluster singles and doubles ground state energy

In this paper, we develop and analyze a number of perturbation series that target the coupled cluster singles and doubles (CCSD) ground state energy. We show how classical Møller-Plesset perturbation theory series can be restructured to target the CCSD energy based on a reference CCS calculation and how the corresponding cluster perturbation series differs from the classical Møller-Plesset perturbation series. Subsequently, we reformulate these series using the coupled cluster Lagrangian framework to obtain series, where fourth and fifth order energies are determined only using parameters through second order. To test the methods, we perform a series of test calculations on molecular photoswitches of both total energies and reaction energies. We find that the fifth order reaction energies are of CCSD quality and that they are of comparable accuracy to state-of-the-art approximations to the CCSD energy based on local pair natural orbitals. The advantage of the present approach over local correlation methods is the absence of user defined threshold parameters for neglecting or approximating contributions to the correlation energy. Fixed threshold parameters lead to discontinuous energy surfaces, although this effect is often small enough to be ignored, but the present approach has a differentiable energy that will facilitate derivation and implementation of gradients and higher derivatives. A further advantage is that the calculation of the perturbation correction is non-iterative and can, therefore, be calculated in parallel, leading to a short time-to-solution.

Cluster-in-Molecule Local Correlation Method for Dispersion Interactions in Large Systems and Periodic Systems

ConspectusThe noncovalent interactions, including dispersion interactions, control the structures and stabilities of complex chemical systems, including host-guest complexes and the adsorption process of molecules on the solid surfaces. The density functional theory (DFT) with empirical dispersion correction is now the working horse in many areas of applications. Post-Hartree-Fock (post-HF) methods have been well recognized to provide more accurate descriptions in a systematic way. However, traditional post-HF methods are mainly limited to small- or medium-sized systems, and their applications to periodic condensed phase systems are still very limited due to their expensive computational costs.To extend post-HF calculations to large molecules, the cluster-in-molecule (CIM) local correlation approach has been established, allowing highly accurate electron correlation calculations that are routinely available for very large systems. In the CIM approach, the electron correlation energy of a large molecule could be obtained from electron correlation calculations on a series of clusters, each of which contains a subset of occupied and virtual localized molecular orbitals. The CIM method could be massively and efficiently parallelized on general computer clusters. The CIM method has been implemented at various electron correlation levels, including second-order Mo̷ller-Plesset perturbation theory (MP2), coupled cluster singles and doubles (CCSD), CCSD with perturbative triples correction [CCSD(T)], etc. The CIM-MP2 energy gradient algorithm was developed and applied to the geometry optimizations of large systems. The CIM method has also been extended to condensed-phase systems under periodic boundary conditions (PBC-CIM). For periodic systems, the correlation energy per unit cell could be evaluated with correlation energy contributions from a series of clusters that are built with localized Wannier functions.CIM-based electron correlation calculations have been employed to investigate a number of chemical problems in which the dispersion interaction is important. CIM-based post-HF methods including CIM domain-based local pair natural orbital (DLPNO) CCSD(T) are applied to compute the relative or binding energies of biological systems or supramolecular complexes, the reaction barrier in a relatively complex chemical reaction. The CIM-MP2 method is used to obtain the optimized geometry of large systems. CIM-based post-HF calculations have also been used to compute the cohesive energies of molecular crystals and adsorption energies of molecules on the solid surfaces. The CIM and its PBC variant are expected to become a powerful theoretical tool for accurate calculations of the energies and structures for a broad range of large systems and condensed-phase systems with significant dispersion interactions.

Local Second-Order Møller-Plesset Theory with a Single Threshold Using Orthogonal Virtual Orbitals: Theory, Implementation, and Assessment

It has long been clear that electron correlation methods exhibit unphysical compute scalings with molecular size, which has motivated the development of local correlation methods to discard effectively zero contributions in a controlled way to yield an approximate correlation energy. The ideal local correlation method should have a single numerical threshold that controls the dropping of terms with the ability to have that threshold set small enough so that the correlation energy is reproduced to enough significant figures such that the result is chemically identical. This work reports such a method for the second-order Møller-Plesset (MP2) theory. The theory, implementation, and testing of this local MP2 theory are reported. Thresholds ranging from 10-5 to 10-8 and basis sets ranging from split valence plus polarization through to quadruple-ζ are assessed for local MP2 calculations on a range of molecules, including linear chains and molecules with two- and three-dimensional character. The implementation is shared memory parallel via OpenMP and yields roughly 50% parallel efficiency with 16 cores for a large job. Considerable efforts were made to minimize memory demands, which increased as thresholds were tightened. A variety of relative energy calculations are presented as a function of threshold to provide some guidance to users on how to obtain adequate precision at a low compute cost. It is particularly clear that derivative properties require tighter thresholds in order to achieve an adequate precision.

Accurate Calculation of Isomerization and Conformational Energies of Larger Molecules Using Explicitly Correlated Local Coupled Cluster Methods in Molpro and ORCA

An overview of the approximations in the explicitly correlated local coupled cluster methods PNO-LCCSD(T)-F12 in Molpro and DLPNO-CCSD(T)F12 in ORCA is given. Options to select the domains of projected atomic orbitals (PAOs), pair natural orbitals (PNOs), and triples natural orbitals (TNOs) in both programs are described and compared in detail. The two programs are applied to compute isomerization and conformational energies of the ISOL24 and ACONFL test sets, where the former is part of the GMTKN55 benchmark suite. Thorough studies of basis set effects are presented for selected systems. These revealed large intramolecular basis set superposition effects that make it practically impossible to reliably determine the complete basis set (CBS) limits without including explicitly correlated terms. The latter strongly reduce the basis set dependence and at the same time also errors caused by the local domain approximations. On the basis of these studies, the PNO-LCCSD(T)-F12 method is applied to determine new reference energies for the above-mentioned benchmark sets. We are confident that our results should agree within a few tenths of a kcal mol-1 with the (unknown) CCSD(T)/CBS values, which therefore allowed us to define computational settings for accurate explicitly correlated local coupled cluster methods with moderate computational effort. With these protocols, especially PNO-LCCSD(T)-F12b/AVTZ', reliable reference values for comprehensive benchmark sets can be generated efficiently. This can significantly advance the development and evaluation of the performance of approximate electronic structure methods, especially improved density functional approximations or machine learning approaches.

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.

DLPNO-MP2 second derivatives for the computation of polarizabilities and NMR shieldings

We present a derivation and efficient implementation of the formally complete analytic second derivatives for the domain-based local pair natural orbital second order Møller-Plesset perturbation theory (MP2) method, applicable to electric or magnetic field-response properties but not yet to harmonic frequencies. We also discuss the occurrence and avoidance of numerical instability issues related to singular linear equation systems and near linear dependences in the projected atomic orbital domains. A series of benchmark calculations on medium-sized systems is performed to assess the effect of the local approximation on calculated nuclear magnetic resonance shieldings and the static dipole polarizabilities. Relative deviations from the resolution of the identity-based MP2 (RI-MP2) reference for both properties are below 0.5% with the default truncation thresholds. For large systems, our implementation achieves quadratic effective scaling, is more efficient than RI-MP2 starting at 280 correlated electrons, and is never more than 5-20 times slower than the equivalent Hartree-Fock property calculation. The largest calculation performed here was on the vancomycin molecule with 176 atoms, 542 correlated electrons, and 4700 basis functions and took 3.3 days on 12 central processing unit cores.

Cluster-in-Molecule Local Correlation Method with an Accurate Distant Pair Correction for Large Systems

The cluster-in-molecule (CIM) local correlation approach with an accurate distant pair correlation energy correction is presented. For large systems, the inclusion of distant pair correlation energies is essential for the accurate prediction of absolute correlation energies and relative energies. Here, we propose a simple and efficient scheme for evaluating the distant pair correlation energy correction for the CIM approaches. The corrections can be readily extracted from electron correlation calculations of clusters with almost no additional effort. Benchmark calculations show that the improved CIM approach can recover more than 99.94% of the correlation energy calculated by the parent method. By combining the CIM approach with the domain-based local pair natural orbital (DLPNO) local correlation approach, we have provided accurate binding energies at the CIM-DLPNO-CCSD(T) level for a test set consisting of eight weakly bound complexes ranging in size from 200 to 1027 atoms. With these results as the reference data, the accuracy and applicability of other electron correlation methods and a few density functional methods for large systems have been assessed.

Scalable Electron Correlation Methods. 8. Explicitly Correlated Open-Shell Coupled-Cluster with Pair Natural Orbitals PNO-RCCSD(T)-F12 and PNO-UCCSD(T)-F12

We present explicitly correlated open-shell pair natural orbital local coupled-cluster methods, PNO-RCCSD(T)-F12 and PNO-UCCSD(T)-F12. The methods are extensions of our previously reported PNO-R/UCCSD methods (J. Chem. Theory Comput., 2020, 16, 3135-3151, https://pubs.acs.org/doi/10.1021/acs.jctc.0c00192) with additions of explicit correlation and perturbative triples corrections. The explicit correlation treatment follows the spin-orbital CCSD-F12b theory using Ansatz 3*A, which is found to yield comparable or better basis set convergence than the more rigorous Ansatz 3C in computed ionization potentials and reaction energies using double- to quaduple-ζ basis sets. The perturbative triples correction is adapted from the spin-orbital (T) theory to use triples natural orbitals (TNOs). To address the coupling due to off-diagonal Fock matrix elements, the local triples amplitudes are iteratively solved using small domains of TNOs, and a semicanonical (T0) domain correction with larger domains is applied to reduce the domain errors. The performance of the methods is demonstrated through benchmark calculations on ionization potentials, radical stabilization energies, reaction energies of fragmentations and rearrangements in radical cations, and spin-state energy differences of iron complexes. For a few test sets where canonical calculations are feasible, PNO-RCCSD(T)-F12 results agree with the canonical ones to within 0.4 kcal mol^-1, and this maximum error is reduced to below 0.2 kcal mol^-1 when large local domains are used. For larger systems, results using different thresholds for the local approximations are compared to demonstrate that 1 kcal mol^-1 level of accuracy can be achieved using our default settings. For a couple of difficult cases, it is demonstrated that the errors from individual approximations are only a fraction of 1 kcal mol^-1, and the overall accuracy of the method does not rely on error compensations. In contrast to canonical calculations, the use of spin-orbitals does not lead to a significant increase of computational time and memory usage in the most expensive steps of PNO-R/UCCSD(T)-F12 calculations. The only exception is the iterative solution of the (T) amplitudes, which can be avoided without significant errors by using a perturbative treatment of the off-diagonal coupling, known as (T1) approximation. For most systems, even the semicanonical approximation (T0) leads only to small errors in relative energies. Our program is well parallelized and capable of computing accurate correlation energies for molecules with 100-200 atoms using augmented triple-ζ basis sets in less than a day of elapsed time on a small computer cluster.

Mechanistic Studies of Bioorthogonal ATP Analogues for Assessment of Histidine Kinase Autophosphorylation

Phosphorylation is an essential protein modification and is most commonly associated with hydroxyl-containing amino acids via an adenosine triphosphate (ATP) substrate. The last decades have brought greater appreciation to the roles that phosphorylation of myriad amino acids plays in biological signaling, metabolism, and gene transcription. Histidine phosphorylation occurs in both eukaryotes and prokaryotes but has been shown to dominate signaling networks in the latter due to its role in microbial two-component systems. Methods to investigate histidine phosphorylation have lagged behind those to study serine, threonine, and tyrosine modifications due to its inherent instability and the historical view that this protein modification was rare. An important strategy to overcome the reactivity of phosphohistidine is the development of substrate-based probes with altered chemical properties that improve modification longevity but that do not suffer from poor recognition or transfer by the protein. Here, we present combined experimental and computational studies to better understand the molecular requirements for efficient histidine phosphorylation by comparison of the native kinase substrate, ATP, and alkylated ATP derivatives. While recognition of the substrates by the histidine kinases is an important parameter for the formation of phosphohistidine derivatives, reaction sterics also affect the outcome. In addition, we found that stability of the resulting phosphohistidine moieties correlates with the stability of their hydrolysis products, specifically with their free energy in solution. Interestingly, alkylation dramatically affects the stability of the phosphohistidine derivatives at very acidic pH values. These results provide critical mechanistic insights into histidine phosphorylation and will facilitate the design of future probes to study enzymatic histidine phosphorylation.

Scalable Electron Correlation Methods. 7. Local Open-Shell Coupled-Cluster Methods Using Pair Natural Orbitals: PNO-RCCSD and PNO-UCCSD

We present well-parallelized local implementations of high-spin open-shell coupled cluster methods with single and double excitations (CCSD) using pair natural orbitals (PNOs). The methods are based on the spin-orbital coupled cluster theory using restricted open-shell Hartree-Fock (ROHF) reference functions. Two variants, namely, PNO-UCCSD and PNO-RCCSD are implemented and compared. In PNO-UCCSD, the coupled cluster amplitudes are spin-unrestricted, while in PNO-RCCSD the linear terms are spin-adapted by a spin-projection approach as described in J. Chem. Phys. 1993, 99, 5219-5227. Near linear scaling of the computational cost with the number of correlated electrons is achieved by applying domain and pair approximations. The PNOs are spin-independent and obtained using a semicanonical spin-restricted MP2 approximation with large domains of projected atomic orbitals (PAOs). The pair approximations of our previously described closed-shell PNO-LCCSD method are carefully revised so that they are compatible to the UCCSD theory, and PNO-UCCSD or PNO-RCCSD calculations for closed-shell molecules yield exactly the same results as corresponding spin-free closed-shell PNO-LCCSD calculations. The convergence of the results with respect to the thresholds and options that control the domain and pair approximations is demonstrated. It is found that large domains are required for the single excitations in open-shell calculations in order to obtain converged results. In general, the errors of relative energies caused by the local approximations can be reduced to below 1 kcal mol^-1, even for difficult cases. Presently, PNO-RCCSD and PNO-UCCSD calculations for molecules with 100-200 atoms and augmented triple-ζ basis sets can be carried out in a few hours of elapsed time using ∼100 CPU cores. In addition, the program is also capable of performing distinguishable cluster (PNO-RDCSD and PNO-UDCSC) calculations. The present work is a critical step in developing fully local open-shell PNO-RCCSD(T)-F12 methods.

The Molpro quantum chemistry package

Molpro is a general purpose quantum chemistry software package with a long development history. It was originally focused on accurate wavefunction calculations for small molecules but now has many additional distinctive capabilities that include, inter alia, local correlation approximations combined with explicit correlation, highly efficient implementations of single-reference correlation methods, robust and efficient multireference methods for large molecules, projection embedding, and anharmonic vibrational spectra. In addition to conventional input-file specification of calculations, Molpro calculations can now be specified and analyzed via a new graphical user interface and through a Python framework.

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.

Comprehensive Benchmark Results for the Domain Based Local Pair Natural Orbital Coupled Cluster Method (DLPNO-CCSD(T)) for Closed- and Open-Shell Systems

In this study we examine the accuracy of domain-based local pair natural orbital coupled cluster theory with single, double, and perturbative triple excitations (DLPNO-CCSD(T)) on a large benchmark data set. To this end, we use the recently published GMTKN55 superset of molecules that contains 1505 relative energies and 2462 single-point calculations. To our knowledge this is the most comprehensive benchmark evaluation of any highly correlated wave function based ab initio method to date. In the first part of the study, canonical CCSD(T) reference calculations were carried out on the entire test set in order to guarantee that the reference data are of uniform quality. Second, DLPNO-CCSD(T) calculations were carried out under identical conditions. The main finding is that with the exception of two data sets, all data sets have a MAD of 0.4 kcal/mol or less and the majority of sets have a MAD of less than 0.2 kcal/mol. For open shells, the accuracy of the DLPNO calculations was significantly improved through an iterative version of the triples correction.

Scalable Electron Correlation Methods. 6. Local Spin-Restricted Open-Shell Second-Order Møller-Plesset Perturbation Theory Using Pair Natural Orbitals: PNO-RMP2

We present a (near) linear scaling implementation of high-spin open-shell Møller-Plesset perturbation theory using pair natural orbitals (PNO-RMP2). The theory is based on a new variant of open-shell MP2 which is fully spin-adapted and uses a single set of spin-free amplitudes, as in closed-shell MP2. This method, denoted SROMP2, is invariant to unitary orbital transformations within the closed, open, and virtual orbital subspaces. Accordingly, only a single set of PNOs per spatial orbital pair is needed, and the efficiency is similar to closed-shell calculations. The PNOs are obtained using a semicanonical approximation with large domains of projected atomic orbitals (PAOs). Linear scaling is achieved provided that the open-shell orbitals are local, and distant pairs are treated by multipole approximations. The method is efficiently parallelized. The convergence of ionization and reaction energies as a function of the PAO and PNO domain sizes is demonstrated and found to be very similar as for closed-shell calculations. The suitability of the PNOs for explicitly correlated PNO-RCCSD-F12 calculations is also tested. So far, this method is only simulated using a conventional program with appropriate projections to the PAO and PNO subspaces. It is demonstrated for radical stabilization energies as well as ionization potentials that the errors caused by the local domain approximations with our default thresholds are negligible.

Crystal Field Splittings in Lanthanide Complexes: Inclusion of Correlation Effects beyond Second Order Perturbation Theory

State-averaged complete active space self-consistent field (CASSCF) calculations and a subsequent spin-orbit calculation mixing the CASSCF wave functions (CASSCF/state-interaction with spin-orbit coupling) is the conventional approach used for ab initio calculations of crystal-field splittings and magnetic properties of lanthanide complexes. However, this approach neglects dynamical correlation. Complete active space second-order perturbation theory (CASPT2) can be used to account for dynamical correlation but suffers from the well-known problems of multireference perturbation theory, e.g., intruder state problems. Variational multireference configuration interaction (MRCI) calculations do not show these problems but are usually not feasible due to the large size of real lanthanide complexes. Here, we present a quasi-local projected internally contracted MRCI approach which makes MRCI calculations of lanthanide complexes feasible and allows assessing the influence of dynamical correlation beyond second-order perturbation theory. We apply the method to two well-studied molecules, namely, [Er{N(SiMe₃)₂}₃] and {C(NH₂)₃}₅[Er(CO₃)₄]·11H₂O.