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

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).

Size-consistent explicitly correlated triple excitation correction.

A new approach is proposed to reduce the basis set incompleteness error of the triple excitation correction in explicitly correlated coupled-cluster singles and doubles with perturbative triples calculations. Our method is similar to the intuitive triples correction approach of Knizia et al. [J. Chem. Phys. 130, 054104 (2009)] but, in contrast to the latter, is size-consistent. The new approximation is easy to implement, and its overhead is negligible with respect to the conventional (T) correction. The performance of the approach is assessed for atomization, reaction, and interaction energies as well as for bond lengths and harmonic vibrational frequencies. The advantages of its size consistency are also demonstrated.

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.

The MRCC program system: Accurate quantum chemistry from water to proteins.

MRCC is a package of ab initio and density functional quantum chemistry programs for accurate electronic structure calculations. The suite has efficient implementations of both low- and high-level correlation methods, such as second-order Møller-Plesset (MP2), random-phase approximation (RPA), second-order algebraic-diagrammatic construction [ADC(2)], coupled-cluster (CC), configuration interaction (CI), and related techniques. It has a state-of-the-art CC singles and doubles with perturbative triples [CCSD(T)] code, and its specialties, the arbitrary-order iterative and perturbative CC methods developed by automated programming tools, enable achieving convergence with regard to the level of correlation. The package also offers a collection of multi-reference CC and CI approaches. Efficient implementations of density functional theory (DFT) and more advanced combined DFT-wave function approaches are also available. Its other special features, the highly competitive linear-scaling local correlation schemes, allow for MP2, RPA, ADC(2), CCSD(T), and higher-order CC calculations for extended systems. Local correlation calculations can be considerably accelerated by multi-level approximations and DFT-embedding techniques, and an interface to molecular dynamics software is provided for quantum mechanics/molecular mechanics calculations. All components of MRCC support shared-memory parallelism, and multi-node parallelization is also available for various methods. For academic purposes, the package is available free of charge.

Integral-Direct and Parallel Implementation of the CCSD(T) Method: Algorithmic Developments and Large-Scale Applications.

A completely integral-direct, disk I/O, and network traffic economic coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] implementation has been developed relying on the density-fitting approximation. By fully exploiting the permutational symmetry, the presented algorithm is highly operation count and memory-efficient. Our measurements demonstrate excellent strong scaling achieved via hybrid MPI/OpenMP parallelization and a highly competitive, 60-70% utilization of the theoretical peak performance on up to hundreds of cores. The terms whose evaluation time becomes significant only for small- to medium-sized examples have also been extensively optimized. Consequently, high performance is also expected for systems appearing in extensive data sets used, e.g., for density functional or machine learning parametrizations, and in calculations required for certain reduced-cost or local approximations of CCSD(T), such as in our local natural orbital scheme [LNO-CCSD(T)]. The efficiency of this implementation allowed us to perform some of the largest CCSD(T) calculations ever presented for systems of 31-43 atoms and 1037-1569 orbitals using only four to eight many-core CPUs and 1-3 days of wall time. The resulting 13 correlation energies and the 12 corresponding reaction energies and barrier heights are added to our previous benchmark set collecting reference CCSD(T) results of molecules at the applicability limit of current implementations.

Reaction fronts of the autocatalytic hydrogenase reaction.

We have built a model to describe the hydrogenase catalyzed, autocatalytic, reversible hydrogen oxidation reaction where one of the enzyme forms is the autocatalyst. The model not only reproduces the experimentally observed front properties, but also explains the found hydrogen ion dependence. Furthermore, by linear stability analysis, two different front types are found in good agreement with the experiments.

Definitive thermochemistry and kinetics of the interconversions among conformers of n-butane and n-pentane.

The focal-point analysis (FPA) technique is used for the definitive characterization of conformational interconversion parameters, including activation energy barriers, activation free energies, and kinetic rate coefficients at 298 K, of two n-alkanes, n-butane, and n-pentane, yielding the first complete analysis of their interconversion kinetics. The FPA implementation developed in this study is based on geometry optimizations and harmonic frequency computations carried out with density functional theory methods and single-point energy computations up to the CCSD(T) level of electronic structure theory using atom-centered Gaussian basis sets as large as cc-pV5Z. The anharmonic vibrational computations are carried out, at the MP2/6-31G* level of theory. Reflecting the convergence behavior of the Gibbs free-energy terms and the interconversion parameters, well-defined uncertainties, mostly neglected in previous theoretical studies, are provided. Finally, the effect of these uncertainties on the concentrations of the conformers of n-butane and n-pentane is examined via a global Monte-Carlo uncertainty analysis. © 2017 Wiley Periodicals, Inc.

Enthalpy Differences of the n-Pentane Conformers.

The energy and enthalpy differences of alkane conformers in various temperature ranges have been the subject for both experimental and theoretical studies over the last few decades. It was shown previously for the conformers of butane [G. Tasi et al., J. Chem. Theory Comput. 2012, 8, 479-486] that quantum chemical results can compete with spectroscopic techniques and results obtained even from the most carefully performed experiments could be biased due to the improper statistical model utilized to evaluate the raw experimental data. In the current study, on one hand, the experimental values and their uncertainties for the enthalpy differences for pentane conformers are re-evaluated using the appropriate statistical model. On the other hand, a coupled-cluster-based focal-point analysis has been performed to calculate energy and enthalpy differences for the conformers of pentane. The model chemistry defined in this study includes contributions up to the perturbative quadruple excitations augmented with further small correction terms beyond the Born-Oppenheimer and nonrelativistic approximations. Benchmark quality energy and enthalpy differences for the pentane conformers are given at temperatures 0 and 298.15 K as well as for the various temperature ranges used in the gas-phase experimental measurements. Furthermore, a slight positive shift for the experimental enthalpy differences is also predicted due to an additional Raman active band belonging to the gauche-gauche conformer.