Gerald Knizia

Molpro: a general-purpose quantum chemistry program package

Molpro (available at http://www.molpro.net) is a general‐purpose quantum chemical program. The original focus was on high‐accuracy wave function calculations for small molecules, but using local approximations combined with explicit correlation treatments, highly accurate coupled‐cluster calculations are now possible for molecules with up to approximately 100 atoms. Recently, multireference correlation treatments were also made applicable to larger molecules. Furthermore, an efficient implementation of density functional theory is available.

A simple and efficient CCSD(T)-F12 approximation

A new explicitly correlated CCSD(T)-F12 approximation is presented and tested for 23 molecules and 15 chemical reactions. The F12 correction strongly improves the basis set convergence of correlation and reaction energies. Errors of the Hartree-Fock contributions are effectively removed by including MP2 single excitations into the auxiliary basis set. Using aug-cc-pVTZ basis sets the CCSD(T)-F12 calculations are more accurate and two orders of magnitude faster than standard CCSD(T)/aug-cc-pV5Z calculations.

Simplified CCSD(T)-F12 methods: Theory and benchmarks

The simple and efficient CCSD(T)-F12x approximations (x = a,b) we proposed in a recent communication [T. B. Adler, G. Knizia, and H.-J. Werner, J. Chem. Phys. 127, 221106 (2007)] are explained in more detail and extended to open-shell systems. Extensive benchmark calculations are presented, which demonstrate great improvements in basis set convergence for a wide variety of applications. These include reaction energies of both open- and closed-shell reactions, atomization energies, electron affinities, ionization potentials, equilibrium geometries, and harmonic vibrational frequencies. For all these quantities, results better than the AV5Z quality are obtained already with AVTZ basis sets, and usually AVDZ treatments reach at least the conventional AVQZ quality. For larger molecules, the additional cost for these improvements is only a few percent of the time for a standard CCSD(T) calculation. For the first time ever, total reaction energies with chemical accuracy are obtained using valence-double-zeta basis sets.

Explicitly correlated RMP2 for high-spin open-shell reference states.

We present an explicitly correlated version of the high-spin open-shell RMP2 method. The theory is derived in a unitarily invariant form, which is suitable for the insertion of local approximations. It is demonstrated that the rapid basis set convergence of closed-shell MP2-F12 is also achieved in RMP2-F12, and similar Ansatze and approximations can be employed. All integrals are computed using efficient density fitting approximations, and many-electron integrals are avoided using resolution of the identity approximations. The performance of the method is demonstrated by benchmark calculations on a large set of ionization potentials, electron affinities and atomization energies. Using triple-zeta basis sets RMP2-F12 yields results that are closer to the basis set limit than standard RMP2 with augmented quintuple-zeta basis sets for all properties. Different variants of perturbative corrections for the open-shell Hartree-Fock treatment are described and tested.

Extrapolating MP2 and CCSD explicitly correlated correlation energies to the complete basis set limit with first and second row correlation consistent basis sets

Accurate extrapolation to the complete basis set (CBS) limit of valence correlation energies calculated with explicitly correlated MP2-F12 and CCSD(T)-F12b methods have been investigated using a Schwenke-style approach for molecules containing both first and second row atoms. Extrapolation coefficients that are optimal for molecular systems containing first row elements differ from those optimized for second row analogs, hence values optimized for a combined set of first and second row systems are also presented. The new coefficients are shown to produce excellent results in both Schwenke-style and equivalent power-law-based two-point CBS extrapolations, with the MP2-F12/cc-pV(D,T)Z-F12 extrapolations producing an average error of just 0.17 mE(h) with a maximum error of 0.49 for a collection of 23 small molecules. The use of larger basis sets, i.e., cc-pV(T,Q)Z-F12 and aug-cc-pV(Q,5)Z, in extrapolations of the MP2-F12 correlation energy leads to average errors that are smaller than the degree of confidence in the reference data (approximately 0.1 mE(h)). The latter were obtained through use of very large basis sets in MP2-F12 calculations on small molecules containing both first and second row elements. CBS limits obtained from optimized coefficients for conventional MP2 are only comparable to the accuracy of the MP2-F12/cc-pV(D,T)Z-F12 extrapolation when the aug-cc-pV(5+d)Z and aug-cc-pV(6+d)Z basis sets are used. The CCSD(T)-F12b correlation energy is extrapolated as two distinct parts: CCSD-F12b and (T). While the CCSD-F12b extrapolations with smaller basis sets are statistically less accurate than those of the MP2-F12 correlation energies, this is presumably due to the slower basis set convergence of the CCSD-F12b method compared to MP2-F12. The use of larger basis sets in the CCSD-F12b extrapolations produces correlation energies with accuracies exceeding the confidence in the reference data (also obtained in large basis set F12 calculations). It is demonstrated that the use of the 3C(D) Ansatz is preferred for MP2-F12 CBS extrapolations. Optimal values of the geminal Slater exponent are presented for the diagonal, fixed amplitude Ansatz in MP2-F12 calculations, and these are also recommended for CCSD-F12b calculations.

Intrinsic Atomic Orbitals: An Unbiased Bridge between Quantum Theory and Chemical Concepts.

Modern quantum chemistry can make quantitative predictions on an immense array of chemical systems. However, the interpretation of those predictions is often complicated by the complex wave function expansions used. Here we show that an exceptionally simple algebraic construction allows for defining atomic core and valence orbitals, polarized by the molecular environment, which can exactly represent self-consistent field wave functions. This construction provides an unbiased and direct connection between quantum chemistry and empirical chemical concepts, and can be used, for example, to calculate the nature of bonding in molecules, in chemical terms, from first principles. In particular, we find consistency with electronegativities (χ), C 1s core-level shifts, resonance substituent parameters (σR), Lewis structures, and oxidation states of transition-metal complexes.

Explicitly correlated multireference configuration interaction: MRCI-F12

An internally contracted multireference configuration interaction is developed which employs wave functions that explicitly depend on the electron-electron distance (MRCI-F12). This MRCI-F12 method has the same applicability as the MRCI method, while having much improved basis-set convergence with little extra computational cost. The F12b approximation is used to arrive at a computationally efficient implementation. The MRCI-F12 method is applied to the singlet-triplet separation of methylene, the dissociation energy of ozone, properties of diatomic molecules, and the reaction barrier and exothermicity of the F + H(2) reaction. These examples demonstrate that already with basis sets of moderate size the method provides near complete basis set MRCI accuracy, and hence quantitative agreement with the experimental data. As a side product, we have also implemented the explicitly correlated multireference averaged coupled pair functional method (MRACPF-F12).

Accurate calculation of vibrational frequencies using explicitly correlated coupled-cluster theory.

The recently proposed explicitly correlated CCSD(T)-F12x (x = a,b) approximations [T. B. Adler, G. Knizia, and H.-J. Werner, J. Chem. Phys. 127, 221106 (2007)] are applied to compute equilibrium structures and harmonic as well as anharmonic vibrational frequencies for H(2)O, HCN, CO(2), CH(2)O, H(2)O(2), C(2)H(2), CH(2)NH, C(2)H(2)O, and the trans-isomer of 1,2-C(2)H(2)F(2). Using aug-cc-pVTZ basis sets, the CCSD(T)-F12a equilibrium geometries and harmonic vibrational frequencies are in very close agreement with CCSD(T)/aug-cc-pV5Z values. The anharmonic frequencies are evaluated using vibrational self-consistent field and vibrational configuration interaction methods based on automatically generated potential energy surfaces. The mean absolute deviation of the CCSD(T)-F12a/aug-cc-pVTZ anharmonic frequencies from experimental values amounts to only 4.0 cm(-1).

Explicitly correlated coupled cluster methods with pair-specific geminals

Explicitly correlated MP2-F12 and CCSD(T)-F12 methods with orbital-pair-specific Slater-type geminals are proposed. The fixed amplitude ansatz of Ten-no is used, and different exponents of the Slater geminal functions can be chosen for core–core, core–valence, and valence–valence pairs. This takes care of the different sizes of the correlation hole and leads to improved results when inner-shell orbitals are correlated. The complications and the extra computational cost as compared to corresponding calculations with a single geminal are minor. The improved accuracy of the method is demonstrated for spectroscopic properties of Br2, As2, Ga2, Cu2, GaCl, CuCl, and CuBr, where the d-orbitals are treated as core.

A new internally contracted multi-reference configuration interaction method

We present a new internally contracted multi-reference configuration interaction (MRCI) method which, at the same time, efficiently handles large active orbital spaces, long configuration expansions, and many closed-shell orbitals in the reference function. This is achieved by treating the closed-shell orbitals explicitly, so that all required coupling coefficients and density matrices only depend on active orbital labels. As a result, closed-shell orbitals are handled as efficiently as in a closed-shell single-reference program, and this opens up the possibility to perform high-accuracy MRCI calculations for much larger molecules than before. The enormously complex equations are derived using a new domain-specific computer algebra system and semi-automatically implemented using a newly developed integrated tensor framework. The accuracy and efficiency of the MRCI method is demonstrated with applications to dioxygen-copper complexes with different ligands, some of which involve more than 30 atoms, and to spin-state splittings of ferrocene.