An efficient and near linear scaling pair natural orbital based local coupled cluster method
- PMID: 23343267
- DOI: 10.1063/1.4773581
An efficient and near linear scaling pair natural orbital based local coupled cluster method
Abstract
In previous publications, it was shown that an efficient local coupled cluster method with single- and double excitations can be based on the concept of pair natural orbitals (PNOs) [F. Neese, A. Hansen, and D. G. Liakos, J. Chem. Phys. 131, 064103 (2009)]. The resulting local pair natural orbital-coupled-cluster single double (LPNO-CCSD) method has since been proven to be highly reliable and efficient. For large molecules, the number of amplitudes to be determined is reduced by a factor of 10(5)-10(6) relative to a canonical CCSD calculation on the same system with the same basis set. In the original method, the PNOs were expanded in the set of canonical virtual orbitals and single excitations were not truncated. This led to a number of fifth order scaling steps that eventually rendered the method computationally expensive for large molecules (e.g., >100 atoms). In the present work, these limitations are overcome by a complete redesign of the LPNO-CCSD method. The new method is based on the combination of the concepts of PNOs and projected atomic orbitals (PAOs). Thus, each PNO is expanded in a set of PAOs that in turn belong to a given electron pair specific domain. In this way, it is possible to fully exploit locality while maintaining the extremely high compactness of the original LPNO-CCSD wavefunction. No terms are dropped from the CCSD equations and domains are chosen conservatively. The correlation energy loss due to the domains remains below <0.05%, which implies typically 15-20 but occasionally up to 30 atoms per domain on average. The new method has been given the acronym DLPNO-CCSD ("domain based LPNO-CCSD"). The method is nearly linear scaling with respect to system size. The original LPNO-CCSD method had three adjustable truncation thresholds that were chosen conservatively and do not need to be changed for actual applications. In the present treatment, no additional truncation parameters have been introduced. Any additional truncation is performed on the basis of the three original thresholds. There are no real-space cutoffs. Single excitations are truncated using singles-specific natural orbitals. Pairs are prescreened according to a multipole expansion of a pair correlation energy estimate based on local orbital specific virtual orbitals (LOSVs). Like its LPNO-CCSD predecessor, the method is completely of black box character and does not require any user adjustments. It is shown here that DLPNO-CCSD is as accurate as LPNO-CCSD while leading to computational savings exceeding one order of magnitude for larger systems. The largest calculations reported here featured >8800 basis functions and >450 atoms. In all larger test calculations done so far, the LPNO-CCSD step took less time than the preceding Hartree-Fock calculation, provided no approximations have been introduced in the latter. Thus, based on the present development reliable CCSD calculations on large molecules with unprecedented efficiency and accuracy are realized.
Similar articles
-
Efficient and accurate approximations to the local coupled cluster singles doubles method using a truncated pair natural orbital basis.J Chem Phys. 2009 Aug 14;131(6):064103. doi: 10.1063/1.3173827. J Chem Phys. 2009. PMID: 19691374
-
A new near-linear scaling, efficient and accurate, open-shell domain-based local pair natural orbital coupled cluster singles and doubles theory.J Chem Phys. 2017 Apr 28;146(16):164105. doi: 10.1063/1.4981521. J Chem Phys. 2017. PMID: 28456208
-
Natural triple excitations in local coupled cluster calculations with pair natural orbitals.J Chem Phys. 2013 Oct 7;139(13):134101. doi: 10.1063/1.4821834. J Chem Phys. 2013. PMID: 24116546
-
New coupled cluster approaches based on the unrestricted Hartree-Fock reference for treating molecules with multireference character.Phys Chem Chem Phys. 2011 May 21;13(19):8795-804. doi: 10.1039/c0cp01103g. Epub 2011 Mar 29. Phys Chem Chem Phys. 2011. PMID: 21448471 Review.
-
Calculations on noncovalent interactions and databases of benchmark interaction energies.Acc Chem Res. 2012 Apr 17;45(4):663-72. doi: 10.1021/ar200255p. Epub 2012 Jan 6. Acc Chem Res. 2012. PMID: 22225511 Review.
Cited by
-
Exploring the potential of natural orbital functionals.Chem Sci. 2024 Oct 10. doi: 10.1039/d4sc05810k. Online ahead of print. Chem Sci. 2024. PMID: 39421199 Free PMC article. Review.
-
On the Potential Energy Surface of the Pyrene Dimer.Int J Mol Sci. 2024 Oct 6;25(19):10762. doi: 10.3390/ijms251910762. Int J Mol Sci. 2024. PMID: 39409090 Free PMC article.
-
Insights into the Structural Modification of Selenium-Doped Derivatives with Narrowband Emissions: A Theory Study.Molecules. 2024 Sep 27;29(19):4589. doi: 10.3390/molecules29194589. Molecules. 2024. PMID: 39407519 Free PMC article.
-
Bidentate Substrate Binding Mode in Oxalate Decarboxylase.Molecules. 2024 Sep 17;29(18):4414. doi: 10.3390/molecules29184414. Molecules. 2024. PMID: 39339409 Free PMC article.
-
Enzymatic synthesis of azide by a promiscuous N-nitrosylase.Nat Chem. 2024 Sep 27. doi: 10.1038/s41557-024-01646-2. Online ahead of print. Nat Chem. 2024. PMID: 39333393
LinkOut - more resources
Full Text Sources
Other Literature Sources