A theoretical model chemistry designed to achieve high accuracy for enthalpies of formation of atoms and small molecules is described. This approach is entirely independent of experimental data and contains no empirical scaling factors, and includes a treatment of electron correlation up to the full coupled-cluster singles, doubles, triples and quadruples approach. Energies are further augmented by anharmonic zero-point vibrational energies, a scalar relativistic correction, first-order spin-orbit coupling, and the diagonal Born-Oppenheimer correction. The accuracy of the approach is assessed by several means. Enthalpies of formation (at 0 K) calculated for a test suite of 31 atoms and molecules via direct calculation of the corresponding elemental formation reactions are within 1 kJ mol(-1) to experiment in all cases. Given the quite different bonding environments in the product and reactant sides of these reactions, the results strongly indicate that even greater accuracy may be expected in reactions that preserve (either exactly or approximately) the number and types of chemical bonds.

We are proposing a new computational thermochemistry protocol denoted W3 theory, as a successor to W1 and W2 theory proposed earlier [Martin and De Oliveira, J. Chem. Phys. 111, 1843 (1999)]. The new method is both more accurate overall (error statistics for total atomization energies approximately cut in half) and more robust (particularly towards systems exhibiting significant nondynamical correlation) than W2 theory. The cardinal improvement rests in an approximate account for post-CCSD(T) correlation effects. Iterative T3 (connected triple excitations) effects exhibit a basis set convergence behavior similar to the T3 contribution overall. They almost universally decrease molecular binding energies. Their inclusion in isolation yields less accurate results than CCSD(T) nearly across the board: It is only when T4 (connected quadruple excitations) effects are included that superior performance is achieved. T4 effects systematically increase molecular binding energies. Their basis set convergence is quite rapid, and even CCSDTQ/cc-pVDZ scaled by an empirical factor of 1.2532 will yield a quite passable quadruples contribution. The effect of still higher-order excitations was gauged for a subset of molecules (notably the eight-valence electron systems): T5 (connected quintuple excitations) contributions reach 0.3 kcal/mol for the pathologically multireference X 1Sigmag+ state of C2 but are quite small for other systems. A variety of avenues for achieving accuracy beyond that of W3 theory were explored, to no significant avail. W3 thus appears to represent a good compromise between accuracy and computational cost for those seeking a robust method for computational thermochemistry in the kJ/mol accuracy range on small systems.

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.

Articles you may be interested inInclusion of selected higher excitations involving active orbitals in the state-specific multireference coupledcluster theory J. Chem. Phys. 133, 234110 (2010); 10.1063/1.3515478Second-and third-order triples and quadruples corrections to coupled-cluster singles and doubles in the ground and excited states Extended benchmark studies of coupled cluster theory through triple excitations The viability of treating higher excitations in coupled-cluster theory is discussed. An algorithm is presented for solving coupled-cluster ͑CC͒ equations which can handle any excitation. Our method combines the formalism of diagrammatic many-body perturbation theory and string-based configuration interaction ͑CI͒. CC equations are explicitly put down in terms of antisymmetrized diagrams and a general method is proposed for the factorization of the corresponding algebraic expressions. Contractions between cluster amplitudes and intermediates are evaluated by a string-based algorithm. In contrast to our previous developments ͓J. Chem. Phys. 113, 1359 ͑2000͔͒ the operation count of this new method scales roughly as the (2nϩ2)nd power of the basis set size where n is the highest excitation in the cluster operator. As a by-product we get a completely new CI formalism which is effective for solving both truncated and full CI problems. Generalization for approximate CC models as well as multireference cases is also discussed.

The recently developed high-accuracy extrapolated ab initio thermochemistry method for theoretical thermochemistry, which is intimately related to other high-precision protocols such as the Weizmann-3 and focal-point approaches, is revisited. Some minor improvements in theoretical rigor are introduced which do not lead to any significant additional computational overhead, but are shown to have a negligible overall effect on the accuracy. In addition, the method is extended to completely treat electron correlation effects up to pentuple excitations. The use of an approximate treatment of quadruple and pentuple excitations is suggested; the former as a pragmatic approximation for standard cases and the latter when extremely high accuracy is required. For a test suite of molecules that have rather precisely known enthalpies of formation {as taken from the active thermochemical tables of Ruscic and co-workers [Lecture Notes in Computer Science, edited by M. Parashar (Springer, Berlin, 2002), Vol. 2536, pp. 25-38; J. Phys. Chem. A 108, 9979 (2004)]}, the largest deviations between theory and experiment are 0.52, -0.70, and 0.51 kJ mol(-1) for the latter three methods, respectively. Some perspective is provided on this level of accuracy, and sources of remaining systematic deficiencies in the approaches are discussed.

A new method is presented for treating the effects of quadruple excitations in coupled-cluster theory. In the approach, quadruple excitation contributions are computed from a formula based on a non-Hermitian perturbation theory analogous to that used previously to justify the usual noniterative triples correction used in the coupled cluster singles and doubles method with a perturbative treatment of the triple excitations (CCSD(T)). The method discussed in this paper plays a parallel role in improving energies obtained with the full coupled-cluster singles, doubles, and triples method (CCSDT) by adding a perturbative treatment of the quadruple excitations (CCSDT(Q)). The method is tested for an extensive set of examples, and is shown to provide total energies that compare favorably with those obtained with the full singles, doubles, triples, and quadruples (CCSDTQ) method.

An improved version of our general-order local coupled-cluster (CC) approach [Z. Rolik and M. Kállay, J. Chem. Phys. 135, 104111 (2011)] and its efficient implementation at the CC singles and doubles with perturbative triples [CCSD(T)] level is presented. The method combines the cluster-in-molecule approach of Li and co-workers [J. Chem. Phys. 131, 114109 (2009)] with frozen natural orbital (NO) techniques. To break down the unfavorable fifth-power scaling of our original approach a two-level domain construction algorithm has been developed. First, an extended domain of localized molecular orbitals (LMOs) is assembled based on the spatial distance of the orbitals. The necessary integrals are evaluated and transformed in these domains invoking the density fitting approximation. In the second step, for each occupied LMO of the extended domain a local subspace of occupied and virtual orbitals is constructed including approximate second-order Mo̸ller-Plesset NOs. The CC equations are solved and the perturbative corrections are calculated in the local subspace for each occupied LMO using a highly-efficient CCSD(T) code, which was optimized for the typical sizes of the local subspaces. The total correlation energy is evaluated as the sum of the individual contributions. The computation time of our approach scales linearly with the system size, while its memory and disk space requirements are independent thereof. Test calculations demonstrate that currently our method is one of the most efficient local CCSD(T) approaches and can be routinely applied to molecules of up to 100 atoms with reasonable basis sets.

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 Møller-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-ζ basis set is feasible on a single processor for a protein molecule including 2380 atoms and more than 44000 atomic orbitals.