Natural Linear-Scaled Coupled-Cluster Theory with Local Transferable Triple Excitations: Applications to Peptides

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.

Third-Order Many-Body Expansion of OSV-MP2 Wave Function for Low-Order Scaling Analytical Gradient Computation

We present a many-body expansion (MBE) formulation and implementation for efficient computation of analytical energy gradients from the orbital-specific-virtual second-order Møllet-Plesset perturbation theory (OSV-MP2) based on our earlier work (Zhou et al. J. Chem. Theory Comput. 2020, 16, 196-210). The third-order MBE(3) expansion of OSV-MP2 amplitudes and density matrices was developed to adopt the orbital-specific clustering and long-range termination schemes, which avoids term-by-term differentiations of the MBE energy bodies. We achieve better efficiency by exploiting the algorithmic sparsity that allows us to prune out insignificant fitting integrals and OSV relaxations. With these approximations, the present implementation is benchmarked on a range of molecules that show an economic scaling in the linear and quadratic regimes for computing MBE(3)-OSV-MP2 amplitude and gradient equations, respectively, and yields normal accuracy comparable to the original OSV-MP2 results. The MPI-3-based parallelism through shared memory one-sided communication is further developed for improving parallel scalability and memory accessibility by sorting the MBE(3) orbital clusters into independent tasks that are distributed on multiple processes across many nodes, supporting both global and local data locations in which selected MBE(3)-OSV-MP2 intermediates of different sizes are distinguished and accordingly placed. The accuracy and efficiency level of our MBE(3)-OSV-MP2 analytical gradient implementation is finally illustrated in two applications: we show that the subtle coordination structure differences of mechanically interlocked Cu-catenane complexes can be distinguished when tuning ligand lengths; and the porphycene molecular dynamics reveals the emergence of the vibrational signature arising from softened N-H stretching associated with hydrogen transfer, using an MP2 level of electron correlation and classical nuclei for the first time.

iOI: An Iterative Orbital Interaction Approach for Solving the Self-Consistent Field Problem

An iterative orbital interaction (iOI) approach is proposed to solve, in a bottom-up fashion, the self-consistent field problem in quantum chemistry. While it belongs grossly to the family of fragment-based quantum chemical methods, iOI is distinctive in that (1) it divides and conquers not only the energy but also the wave function and that (2) the subsystem sizes are automatically determined by successively merging neighboring small subsystems until they are just enough for converging the wave function to a given accuracy. Orthonormal occupied and virtual localized molecular orbitals are obtained in a natural manner, which can be used for all post-SCF purposes.

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.

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.

Perturbation Improved Natural Linear-Scaled Coupled-Cluster Method and Its Application to Conformational Analysis

The fragment-based coupled-cluster (CC) theory utilizing the transferable functional groups through natural localized molecular orbital (NLMO), that is, the natural linear-scaled coupled-cluster (NLSCC) has been further developed to take the extra-fragment interactions into account. The correction to the interaction energies sacrificed during the fragmentation process for the previous NLSCC method is computed by a computationally efficient perturbation theory that maintains the original scaling. The new linear-scaled coupled-cluster for the singles and doubles (CCSD) method is applied to the analysis of relative energies of delicate conformational problems of polypeptides. By adding a perturbation correction, results accurate to less than a kcal/mol are obtained for the alanine tetramer.

Assessment of DFT for endohedral complexes' dipole moment: PNO-LCCSD-F12 as a reference method

We present a systematic evaluation of the performance of a wide range of exchange-correlation functionals and related dispersion correction schemes for the computation of dipole moments of endohedral complexes, formed through the encapsulation of an AB molecule (AB = LiF, HCl) inside carbon nanotubes (CNTs) of different diameter. The consistency and accuracy of (i) generalized gradient approximation, (ii) meta GGA, (iii) global hybrid, and (iv) range-separated hybrid density functionals are assessed. In total, 37 density functionals are tested. The results obtained using the highly accurate pair natural orbitals based explicitly correlated local coupled cluster singles doubles (PNO-LCCSD-F12) method of Werner and co-workers [Schwilk et al., J. Chem. Theory Comput., 2017, 13, 3650; Ma et al., J. Chem. Theory Comput., 2017, 13, 4871] with the aug-cc-pVTZ basis set serve as a reference. The static electric dipole moment is computed via the finite field response or, when possible, as the expectation value of the dipole operator. Among others, it is shown that functionals belonging to the class of range-separated hybrids, provide results closest to the coupled cluster reference data. In particular, the ωB97X as well as the M11 functional may be considered as a promising choice for computing electric properties of noncovalent endohedral complexes. On the other hand, the worst performance was found for the functionals which do not include the Hartree-Fock exchange. The analysis of both the coupled cluster and the DFT results indicates a strong coupling of dispersion and polarization that may also explain why lower level DFT methods, as well as Hartree-Fock and MP2, cannot yield dipole moments beyond a qualitative agreement with the higher order reference data. Interestingly, the much smaller and less systematically constructed basis sets of Pople of moderate size provide results of accuracy at least comparable with the extended Dunning's aug-cc-pVTZ basis set.

Communication: An improved linear scaling perturbative triples correction for the domain based local pair-natural orbital based singles and doubles coupled cluster method [DLPNO-CCSD(T)]

In this communication, an improved perturbative triples correction (T) algorithm for domain based local pair-natural orbital singles and doubles coupled cluster (DLPNO-CCSD) theory is reported. In our previous implementation, the semi-canonical approximation was used and linear scaling was achieved for both the DLPNO-CCSD and (T) parts of the calculation. In this work, we refer to this previous method as DLPNO-CCSD(T₀) to emphasize the semi-canonical approximation. It is well-established that the DLPNO-CCSD method can predict very accurate absolute and relative energies with respect to the parent canonical CCSD method. However, the (T₀) approximation may introduce significant errors in absolute energies as the triples correction grows up in magnitude. In the majority of cases, the relative energies from (T₀) are as accurate as the canonical (T) results of themselves. Unfortunately, in rare cases and in particular for small gap systems, the (T₀) approximation breaks down and relative energies show large deviations from the parent canonical CCSD(T) results. To address this problem, an iterative (T) algorithm based on the previous DLPNO-CCSD(T₀) algorithm has been implemented [abbreviated here as DLPNO-CCSD(T)]. Using triples natural orbitals to represent the virtual spaces for triples amplitudes, storage bottlenecks are avoided. Various carefully designed approximations ease the computational burden such that overall, the increase in the DLPNO-(T) calculation time over DLPNO-(T₀) only amounts to a factor of about two (depending on the basis set). Benchmark calculations for the GMTKN30 database show that compared to DLPNO-CCSD(T₀), the errors in absolute energies are greatly reduced and relative energies are moderately improved. The particularly problematic case of cumulene chains of increasing lengths is also successfully addressed by DLPNO-CCSD(T).

Combining Accuracy and Efficiency: An Incremental Focal-Point Method Based on Pair Natural Orbitals

In this work, we present a new pair natural orbitals (PNO)-based incremental scheme to calculate CCSD(T) and CCSD(T0) reaction, interaction, and binding energies. We perform an extensive analysis, which shows small incremental errors similar to previous non-PNO calculations. Furthermore, slight PNO errors are obtained by using T_PNO = T_TNO with appropriate values of 10^-7 to 10^-8 for reactions and 10^-8 for interaction or binding energies. The combination with the efficient MP2 focal-point approach yields chemical accuracy relative to the complete basis-set (CBS) limit. In this method, small basis sets (cc-pVDZ, def2-TZVP) for the CCSD(T) part are sufficient in case of reactions or interactions, while some larger ones (e.g., (aug)-cc-pVTZ) are necessary for molecular clusters. For these larger basis sets, we show the very high efficiency of our scheme. We obtain not only tremendous decreases of the wall times (i.e., factors >10²) due to the parallelization of the increment calculations as well as of the total times due to the application of PNOs (i.e., compared to the normal incremental scheme) but also smaller total times with respect to the standard PNO method. That way, our new method features a perfect applicability by combining an excellent accuracy with a very high efficiency as well as the accessibility to larger systems due to the separation of the full computation into several small increments.

Molecular Dipole Moments within the Incremental Scheme Using the Domain-Specific Basis-Set Approach

We present the first implementation of the fully automated incremental scheme for CCSD unrelaxed dipole moments using the domain-specific basis-set approach. Truncation parameters are varied, and the accuracy of the method is statistically analyzed for a test set of 20 molecules. The local approximations introduce small errors at second order and negligible ones at third order. For a third-order incremental CCSD expansion with a CC2 error correction, a cc-pVDZ/SV domain-specific basis set (tmain = 3.5 Bohr), and the truncation parameter f = 30 Bohr, we obtain a mean error of 0.00 mau (-0.20 mau) and a standard deviation of 1.95 mau (2.17 mau) for the total dipole moments (Cartesian components of the dipole vectors). By analyzing incremental CCSD energies, we demonstrate that the MP2 and CC2 error correction schemes are an exclusive correction for the domain-specific basis-set error. Our implementation of the incremental scheme provides fully automated computations of highly accurate dipole moments at reduced computational cost and is fully parallelized in terms of the calculation of the increments. Therefore, one can utilize the incremental scheme, on the same hardware, to extend the basis set in comparison to standard CCSD and thus obtain a better total accuracy.

Third-Order Incremental Dual-Basis Set Zero-Buffer Approach for Large High-Spin Open-Shell Systems

The third-order incremental dual-basis set zero-buffer approach (inc3-db-B0) is an efficient, accurate, and black-box quantum chemical method for obtaining correlation energies of large systems, and it has been successfully applied to many real chemical problems. In this work, we extend this approach to high-spin open-shell systems. In the open-shell approach, we will first decompose the occupied orbitals of a system into several domains by a K-means clustering algorithm. The essential part is that we preserve the active (singly occupied) orbitals in all the calculations of the domain correlation energies. The duplicated contributions of the active orbitals to the correlation energy are subtracted from the incremental expansion. All techniques of truncating the virtual space such as the B0 approximation can be applied. This open-shell inc3-db-B0 approach is combined with the CCSD and CCSD(T) methods and applied to the computations of a singlet-triplet gap and an electron detachment process. Our approach exhibits an accuracy better than 0.6 kcal/mol or 0.3 eV compared with the standard implementation, while it saves a large amount of the computational time and can be efficiently parallelized.

Incremental evaluation of coupled cluster dipole polarizabilities

In this work we present the first implementation of the incremental scheme for coupled cluster linear-response frequency-dependent dipole polarizabilities.

On the relationship between bond-length alternation and many-electron self-interaction error

Predicting accurate bond-length alternations (BLAs) in long conjugated molecular chains has been a major challenge for electronic-structure theory for many decades. While Hartree-Fock (HF) overestimates BLA significantly, second-order perturbation theory and commonly used density functional theory (DFT) approaches typically underestimate it. Here, we discuss how this failure is related to the many-electron self-interaction error (MSIE), which is inherent to both HF and DFT approaches. We use tuned long-range corrected hybrids to minimize the MSIE for a series of polyenes. The key result is that the minimization of the MSIE alone does not yield accurate BLAs. On the other hand, if the range-separation parameter is tuned to yield accurate BLAs, we obtain a significant MSIE that grows with chain length. Our findings demonstrate that reducing the MSIE is one but not the only important aspect necessary to obtain accurate BLAs from density functional theory.

Implementation of the incremental scheme for one-electron first-order properties in coupled-cluster theory

A fully automated parallelized implementation of the incremental scheme for coupled-cluster singles-and-doubles (CCSD) energies has been extended to treat molecular (unrelaxed) first-order one-electron properties such as the electric dipole and quadrupole moments. The convergence and accuracy of the incremental approach for the dipole and quadrupole moments have been studied for a variety of chemically interesting systems. It is found that the electric dipole moment can be obtained to within 5% and 0.5% accuracy with respect to the exact CCSD value at the third and fourth orders of the expansion, respectively. Furthermore, we find that the incremental expansion of the quadrupole moment converges to the exact result with increasing order of the expansion: the convergence of nonaromatic compounds is fast with errors less than 16 mau and less than 1 mau at third and fourth orders, respectively (1 mau=10(-3)ea(0)(2)); the aromatic compounds converge slowly with maximum absolute deviations of 174 and 72 mau at third and fourth orders, respectively.

Fully Automated Implementation of the Incremental Scheme for Correlation Energies

The applicability and performance of a fully automated implementation of the incremental scheme for the evaluation of correlation contributions in wavefunction-based quantum chemical calculations on large molecules is briefly reviewed and a typical example from organic chemistry is discussed in more detail. The accuracy of relative energies is analyzed for various water hexamers at the CCSD(T)/aug-cc-pVDZ level. Furthermore the accuracy of the potential energy surface of octane is analyzed at the CCSD level of theory. Finally, CCSD(T)/aug-cc-pVTZ calculations with 1656 basis functions in C1 symmetry are reported for (H2O)18, where a full calculation is infeasible within reasonable time.

Accurate benchmark calculations on the gas-phase basicities of small molecules

Accurate benchmark calculations of gas-phase basicities of small molecules are presented and compared with available experimental results. The optimized geometries and thermochemical analyses were obtained from MP2/aug-cc-pVTZ calculations. Two different ab initio electron-correlated methods MP2 and CCSD(T) were employed for subsequent gas-phase basicity calculations, and the single-point energies were extrapolated to the complete basis set (CBS) limit. The overall accuracy for different ab initio methods is compared, and the accuracy in descending order is CCSD(T)_CBS > CCSD(T)/aug-cc-pVDZ > (MP2/aug-cc-pVQZ approximately MP2_CBS) > HF/aug-cc-pVQZ. The best root-mean-squared-error obtained was 1.0 kcal mol(-1) at the CCSD(T)_CBS//MP2/aug-cc-pVTZ level for a test set of 41 molecules. Clearly, accurate calculations for the electron correlation energy are important for the theoretical prediction of molecular gas-phase basicities. However, conformational effects were also found to be relevant in several instances when more complicated molecules were examined.