4.7 Article

Optimization of the Linear-Scaling Local Natural Orbital CCSD(T) Method: Improved Algorithm and Benchmark Applications

Journal

JOURNAL OF CHEMICAL THEORY AND COMPUTATION
Volume 14, Issue 8, Pages 4193-4215

Publisher

AMER CHEMICAL SOC
DOI: 10.1021/acs.jctc.8b00442

Keywords

-

Funding

  1. National Research, Development, and Innovation Office (NKFIH) [KKP126451]
  2. BME-Biotechnology FIKP grant of EMMI (BME FIKP-BIO)
  3. New National Excellence Program of the Ministry of Human Capacities [UNKP-17-4-BME-55]
  4. BME HPC Cluster [TAMOP-4.2.2.B-10/1-2010-0009]

Ask authors/readers for more resources

An optimized implementation of the local natural orbital (LNO) coupled-cluster (CC) with single-, double-, and perturbative triple excitations [LNO-CCSD(T)] method is presented. The integral-direct, in-core, highly efficient domain construction technique of our local second-order Moller-Plesset (LMP2) scheme is extended to the CC level. The resulting scheme, which is also suitable for general-order LNO-CC calculations, inherits the beneficial properties of the LMP2 approach, such as the asymptotically linear-scaling operation count, the asymptotically constant data storage requirement, and the completely independent domain calculations. In addition to integrating our recent redundancy-free LMP2 and Laplace-transformed (T) algorithms with the LNO-CCSD(T) code, the memory demand, the domain and LNO construction, the auxiliary basis compression, and the previously rate-determining two external integral transformation have been significantly improved. The accuracy of all of the approximations is carefully examined on medium-sized to large systems to determine reasonably tight default truncation thresholds. Our benchmark calculations, performed on molecules of up to 63 atoms, show that the optimized method with the default settings provides average correlation and reaction energy errors of less than 0.07% and 0.34 kcal/mol, respectively, compared to the canonical CCSD(T) reference. The efficiency of the present LNO-CCSD(T) implementation is demonstrated on realistic, three-dimensional examples. Using the new code, an LNO-CCSD(T) correlation energy calculation with a triple-sigma basis set is feasible on a single processor for a protein molecule including 2380 atoms and more than 44000 atomic orbitals.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.7
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

No Data Available
No Data Available