A progress report on the status of the COLUMBUS MRCI program system

The COLUMBUS program system is a collection of Fortran programs for performing general multireference single- and double-excitation configuration interaction (MRSDCI) wave function optimization based on the graphical unitary group approach. The program system also includes integral generation, SCF and MCSCF orbital optimization, integral transformation, and wave function analysis programs. The original program system was written in 1980 to 1981. Since that time, it has evolved into one of the most popular MRSDCI program systems used in the computational chemistry community. The discussion of this evolution will include the exploitation of efficient matrix-matrix and matrix-vector product computational kernels, the use of generally contracted symmetry-adapted orbital basis sets, general Hamiltonian diagonalization procedures, energy-based internal walk selection, flexible DRT specification, improved coupling-coefficient evaluation methods, coupled-pair functional and multireference CPF capabilities, and density matrix construction. The numerous versions of the program system that are maintained at different sites and on different computers are now in the process of being merged. The source code for this combined version will be made available to the computational chemistry community. The source code for a specific computer may be generated from the source code for another computer by a single pass through a simple filter utility that is included with the program system. The directly supported computers will initially include various models of VAX, Cray, FPS, IBM, CDC, and ETA machines with the addition of other machines shortly thereafter. The ongoing developments of the COLUMBUS system that are discussed include a new method for computing analytic energy gradients for MRSDCI wave functions. This effective-density-matrix based method avoids the “coupled perturbed MCSCF” solutions for each coordinate direction, avoids the transformation of any derivative-integral quantities from the AO to the MO basis, avoids the transformation of the coupling coefficients from the MO to the AO basis, allows a subset of the MCSCF doubly occupied orbitals to be frozen in the CI wave function, and allows the MRSDCI wave function to be generated from general reference CSFs that are not necessarily related to the MCSCF expansion CSFs.

[1]  Brian T. Smith.  A review and analysis of Fortran 8x , 1988, SGNM.

[2]  Robert J. Gdanitz,et al.  The averaged coupled-pair functional (ACPF): A size-extensive modification of MR CI(SD) , 1988 .

[3]  Distribution of mathematical software via electronic mail , 1987, CACM.

[4]  Peter R. Taylor,et al.  General contraction of Gaussian basis sets. I. Atomic natural orbitals for first‐ and second‐row atoms , 1987 .

[5]  F. B. Brown,et al.  Configuration selection and extrapolation in multireference configuration interaction calculations: The (H2)2 van der waals complex as a benchmark example , 1987 .

[6]  Ron Shepard,et al.  Geometrical energy derivative evaluation with MRCI wave functions , 1987 .

[7]  K. P. Lawley,et al.  Ab initio methods in quantum chemistry , 1987 .

[8]  Poul Jørgensen,et al.  Geometrical derivatives of energy surfaces and molecular properties , 1986 .

[9]  R. Amos,et al.  On the efficient evaluation of analytic energy gradients , 1985 .

[10]  Hans Lischka,et al.  Implementation of an electronic structure program system on the CYBER 205 , 1985 .

[11]  Claus Ehrhardt,et al.  The coupled pair functional (CPF). A size consistent modification of the CI(SD) based on an energy functional , 1985 .

[12]  Henry F. Schaefer,et al.  On the evaluation of analytic energy derivatives for correlated wave functions , 1984 .

[13]  I. Shavitt.  The Treatment of Electron Correlation: Where do We Go from Here? , 1984 .

[14]  Larry McMurchie,et al.  CALCULATION OF INTEGRALS OVER AB INITIO PSEUDOPOTENTIALS , 1981 .

[15]  Charles L. Lawson,et al.  Basic Linear Algebra Subprograms for Fortran Usage , 1979, TOMS.

[16]  Isaiah Shavitt.  The utilization of abelian point group symmetry in the graphical unitary group approach to the calculation of correlated electronic wavefunctions , 1979 .

[17]  M. Inokuti,et al.  Radiological and Environmental Research Division annual report: Fundamental Molecular Physics and Chemistry, October 1977-September 1978. [Summary of research activities at Argonne National Laboratory] , 1978 .

[18]  Isaiah Shavitt,et al.  The Method of Configuration Interaction , 1977 .

[19]  Russell M. Pitzer,et al.  An SCF method for hole states , 1976 .

[20]  Michel Dupuis,et al.  Evaluation of molecular integrals over Gaussian basis functions , 1976 .

[21]  R. Raffenetti,et al.  General contraction of Gaussian atomic orbitals: Core, valence, polarization, and diffuse basis sets; Molecular integral evaluation , 1973 .

[22]  R. Pitzer.  Contribution of atomic orbital integrals to symmetry orbital integrals , 1973 .