Hans-Joachim Werner

An efficient internally contracted multiconfiguration–reference configuration interaction method

A new internally contracted direct multiconfiguration–reference configuration interaction (MRCI) method is described which allows the use of much larger reference spaces than any previous MRCI method. The configurations with two electrons in the external orbital space are generated by applying pair excitation operators to the reference wave function as a whole, while the singly external and internal configurations are standard uncontracted spin eigenfunctions. A new efficient and simple method for the calculation of the coupling coefficients is used, which is well suited for vector machines, and allows the recalculation of all coupling coefficients each time they are needed. The vector H⋅c is computed partly in a nonorthogonal configuration basis. In order to test the accuracy of the internally contracted wave functions, benchmark calculations have been performed for F−, H2O, NH2, CH2, CH3, OH, NO, N2, and O2 at various geometries. The deviations of the energies obtained with internally contracted and unc...

A second order multiconfiguration SCF procedure with optimum convergence

An MCSCF procedure is described which is based on the direct minimization of an approximate energy expression which is periodic and correct to second order in the changes in the orthonormal orbitals. Within this approximation, the CI coefficients are fully optimized, thereby accounting for the coupling between orbital rotations and CI coefficients to higher order than in previous treatments. Additional transformations among the internal orbitals and their associated one‐ and two‐electron integrals are performed which amounts to treating the rotations among internal orbitals to higher than second order. These extra steps are cheap compared to the four index transformation performed in each iteration, but lead to a remarkable enhancement of convergence and overall efficiency. In all calculations attempted to date, convergence has been achieved in at most three iterations. The energy has been observed to converge better than quadratically from the first iteration even when the initial Hessian matrix has many...

An efficient method for the evaluation of coupling coefficients in configuration interaction calculations

Abstract A new method for evaluating one-particle coupling coefficients in a general configuration interaction calculation is presented. Through repeated application and use of resolutions of the identity, two-, three- and four-body coupling coefficients and density matrices may be built in a simple and efficient way. The method is therefore of use in both multiconfiguration SCF (MC SCF) and multireference configuration interaction (MRCI) calculations. Examples show that the approach is efficient for both these applications.

An efficient second-order MC SCF method for long configuration expansions

A new second-order optimisation procedure for general MC SCF wavefunctions is described. The method shows greatly improved convergence as compared to previous methods. Using a determinant-based direct CI procedure which avoids the construction of a formula tape, very long complete active space (CAS SCF) wavefunctions can be handled. Energy averages of several states can also be optimised. Sample calculations for CH2, FeO, and the vinoxy radical CH2CHO with up to 178916 configurations are presented.

Coupled cluster theory for high spin, open shell reference wave functions

The coupled cluster method restricted to single and double excitations (CCSD) is considered for the case of a spin restricted Hartree–Fock open shell reference determinant. A spin–orbital based formulation, in which the cluster operator spans exactly the minimal first order interacting space, is presented, and computationally optimal working equations are given. In the limit of a large number of closed shell orbitals, the cost is identical to that of an optimum treatment of an equivalent closed shell problem, which is obtained as a special case of the formulation presented. The theory is applied to the calculation of a number of diatomic potential energy functions and compared with spin‐unrestricted theory.

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 comparison of the efficiency and accuracy of the quadratic configuration interaction (QCISD), coupled cluster (CCSD), and Brueckner coupled cluster (BCCD) methods

Abstract The coupled-cluster method restricted to single and double excitations from a closed-shell reference function (CCSD) and the corresponding quadratic configuration interaction method (QCISD) are formulated in terms of quantities which can be computed directly from the two-electron integrals in AO basis. A simple yet effective method to accelerate convergence in Brueckner coupled-cluster (BCCD) calculations is also described. Using this procedure BCCD calculations require no more effort than CCSD calculations. In order to compare the accuracy of all three methods potential energy functions and spectroscopic constants have been calculated for N 2 , CO, F 2 , and HF using large basis sets.

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.

Systematically convergent basis sets for explicitly correlated wavefunctions: The atoms H, He, B–Ne, and Al–Ar

Correlation consistent basis sets have been optimized for use with explicitly correlated F12 methods. The new sets, denoted cc-pVnZ-F12 (n=D,T,Q), are similar in size and construction to the standard aug-cc-pVnZ and aug-cc-pV(n+d)Z basis sets, but the new sets are shown in the present work to yield much improved convergence toward the complete basis set limit in MP2-F12/3C calculations on several small molecules involving elements of both the first and second row. For molecules containing only first row atoms, the smallest cc-pVDZ-F12 basis set consistently recovers nearly 99% of the MP2 valence correlation energy when combined with the MP2-F12/3C method. The convergence with basis set for molecules containing second row atoms is slower, but the new DZ basis set still recovers 97%-99% of the frozen core MP2 correlation energy. The accuracy of the new basis sets for relative energetics is demonstrated in benchmark calculations on a set of 15 chemical reactions.