Implementation of High-Order Multireference Coupled-Cluster Methods on Intel Many Integrated Core Architecture

Accurate Reduced-Cost CCSD(T) Energies: Parallel Implementation, Benchmarks, and Large-Scale Applications

The accurate and systematically improvable frozen natural orbital (FNO) and natural auxiliary function (NAF) cost-reducing approaches are combined with our recent coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] implementations. Both of the closed- and open-shell FNO-CCSD(T) codes benefit from OpenMP parallelism, completely or partially integral-direct density-fitting algorithms, checkpointing, and hand-optimized, memory- and operation count effective implementations exploiting all permutational symmetries. The closed-shell CCSD(T) code requires negligible disk I/O and network bandwidth, is MPI/OpenMP parallel, and exhibits outstanding peak performance utilization of 50-70% up to hundreds of cores. Conservative FNO and NAF truncation thresholds benchmarked for challenging reaction, atomization, and ionization energies of both closed- and open-shell species are shown to maintain 1 kJ/mol accuracy against canonical CCSD(T) for systems of 31-43 atoms even with large basis sets. The cost reduction of up to an order of magnitude achieved extends the reach of FNO-CCSD(T) to systems of 50-75 atoms (up to 2124 atomic orbitals) with triple- and quadruple-ζ basis sets, which is unprecedented without local approximations. Consequently, a considerably larger portion of the chemical compound space can now be covered by the practically "gold standard" quality FNO-CCSD(T) method using affordable resources and about a week of wall time. Large-scale applications are presented for organocatalytic and transition-metal reactions as well as noncovalent interactions. Possible applications for benchmarking local CCSD(T) methods, as well as for the accuracy assessment or parametrization of less complete models, for example, density functional approximations or machine learning potentials, are also outlined.

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook.

Perturbative universal state-selective correction for state-specific multi-reference coupled cluster methods

In this work, we report an extension of our previous development of the universal state-selective (USS) multireference coupled-cluster (MRCC) formalism. It was shown [Brabec et al., J. Chem. Phys. 136, 124102 (2012)] and [Banik et al., J. Chem. Phys. 142, 114106 (2015)] that the USS(2) approach significantly improves the accuracy of Brillouin-Wigner and Mukherjee MRCC formulations, however, the numerical and storage costs associated with calculating highly excited intermediates pose a significant challenge, which can restrict the applicability of the USS(2) method. Therefore, we introduce a perturbative variant of the USS(2) approach (USS(pt)), which substantially reduces numerical overhead of the full USS(2) correction while preserving its accuracy. Since the new USS(pt) implementation calculates the triple and quadruple projections in on-the-fly manner, the memory bottleneck associated with the need of storing expensive recursive intermediates is entirely eliminated. On the example of several benchmark systems, we demonstrate accuracies of USS(pt) and USS(2) approaches and their efficiency in describing quasidegenerate electronic states. It is also shown that the USS(pt) method significantly alleviates problems associated with the lack of invariance of MRCC theories upon the rotation of active orbitals.