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. © 2011 John Wiley & Sons, Ltd. This article is categorized under: Software > Quantum Chemistry

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 uncontracted MRCI wave functions are mostly smaller than 1 mH and typically 3–5 times smaller than the deviations between the uncontracted MRCI and the full CI. Dipole moments, electric dipole polarizabilities, and electronic dipole transition moments calculated with uncontracted and contracted MRCI wave functions also are found to be in close agreement. The efficiency of the method is demonstrated in large scale calculations for the CN, NH3, CO2, and Cr2 molecules. In these calculations up to 3088 reference configurations and up to 154 orbitals were employed. The biggest calculation is equivalent to an uncontracted MRCI with more than 78 million configurations.

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 negative eigenvalues.

No abstract

No abstract