Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2023 Jan 10;19(1):174-189.
doi: 10.1021/acs.jctc.2c01031. Epub 2022 Dec 28.

Basis Set Limit CCSD(T) Energies for Extended Molecules via a Reduced-Cost Explicitly Correlated Approach

Affiliations

Basis Set Limit CCSD(T) Energies for Extended Molecules via a Reduced-Cost Explicitly Correlated Approach

Mihály Kállay et al. J Chem Theory Comput. .

Abstract

Several approximations are introduced and tested to reduce the computational expenses of the explicitly correlated coupled-cluster singles and doubles with perturbative triples [CCSD(T)] method for both closed and open-shell species. First, the well-established frozen natural orbital (FNO) technique is adapted to explicitly correlated CC approaches. Second, our natural auxiliary function (NAF) scheme is employed to reduce the size of the auxiliary basis required for the density fitting approximation regularly used in explicitly correlated calculations. Third, a new approach, termed the natural auxiliary basis (NAB) approximation, is proposed to decrease the size of the auxiliary basis needed for the expansion of the explicitly correlated geminals. The performance of the above approximations and that of the combined FNO-NAF-NAB approach are tested for atomization and reaction energies. Our results show that overall speedups of 7-, 5-, and 3-times can be achieved with double-, triple-, and quadruple-ζ basis sets, respectively, without any loss in accuracy. The new method can provide, e.g., reaction energies and barrier heights well within chemical accuracy for molecules with more than 40 atoms within a few days using a few dozen processor cores, and calculations with 50+ atoms are still feasible. These routinely affordable computations considerably extend the reach of explicitly correlated CCSD(T).

PubMed Disclaimer

Conflict of interest statement

The authors declare no competing financial interest.

Figures

Figure 1
Figure 1
Palladium catalyzed C–H bond activation leading to the 2-acetaminobenzophenon product of 31 atoms.
Figure 2
Figure 2
Transition state (TS1) and an intermediate (OO) of 40 atoms formed in a model organocatalytic Michael-addition reaction.
Figure 3
Figure 3
Transition state of a halocyclization reaction containing 53 atoms.,

Similar articles

Cited by

References

    1. Purvis G. D. III; Bartlett R. J. A full coupled-cluster singles and doubles model: The inclusion of disconnected triples. J. Chem. Phys. 1982, 76, 1910.10.1063/1.443164. - DOI
    1. Raghavachari K.; Trucks G. W.; Pople J. A.; Head-Gordon M. A fifth-order perturbation comparison of electron correlation theories. Chem. Phys. Lett. 1989, 157, 479.10.1016/S0009-2614(89)87395-6. - DOI
    1. Adamowicz L.; Bartlett R. J. Optimized virtual orbital space for high-level correlated calculations. J. Chem. Phys. 1987, 86, 6314.10.1063/1.452468. - DOI
    1. Adamowicz L.; Bartlett R. J.; Sadlej J. Optimized virtual orbital space for high-level correlated calculations. II. Electric properties. J. Chem. Phys. 1988, 88, 5749.10.1063/1.454721. - DOI
    1. Sosa C.; Geertsen J.; Trucks G. W.; Bartlett R. J.; Franz J. A. Selection of the reduced virtual space for correlated calculations. An application to the energy and dipole moment of H2O. Chem. Phys. Lett. 1989, 159, 148.10.1016/0009-2614(89)87399-3. - DOI