Symmetry simplifications of space types in configuration interaction induced by orbital degeneracy

Symmetry simplifications are introduced in configuration interaction (CI ) by reducing the number of symmetry-allowed space types if there is degeneracy in some of the molecular orbitals by constructing the unique space types. A new symmetry group which we call the configuration symmetry group is defined and is shown to be expressible as a generalized wreath product group. Generating functions are derived for enumerating the equivalence classes of space types. A double coset method is expounded which constructs the representatives of all equivalence classes of space types using the cycle index of generalized wreath product and the double cosets of label subgroup with generalized wreath product in the symmetric group S_n, if n is twice the number of occupied and virtual orbitals. Method is illustrated with CI using the localized orbitals of polyenes, CI in benzene, and atomic CI for several reference states.     

Study of localized molecular orbitals using group theory methods and its approach to the many-electron correlation problem. IV. The symmetry-adaptation of many-center integrals and hamiltonian matrix elements in MCSCF calculations

Group theoretic methods are presented for the transformations of integrals and the evaluation of matrix elements encountered in multiconfigurational self-consistent field (MCSCF) and configuration interaction (CI) calculations. The method has the advantages of needing only to deal with a symmetry unique set of atomic orbitals (AO) integrals and transformation from unique atomic integrals to unique molecular integrals rather than with all of them. Hamiltonian matrix element is expressed by a linear combination of product terms of many-center unique integrals and geometric factors. The group symmetry localized orbitals as atomic and molecular orbitals are a key feature of this algorithm. The method provides an alternative to traditional method that requires a table of coupling coefficients for products of the irreducible representations of the molecular point group. Geometric factors effectively eliminate these coupling coefficients. The saving of time and space in integral computations and transformations is analyzed. © 1994 by John Wiley & Sons, Inc.     

Direct configuration interaction with a reference state composed of many reference configurations

The configuration interaction method where a single reference state is composed of a linear combination of reference configurations is analyzed in detail. In this method single and double replacements are constructed by applying annihilation and creation operators on the reference state. The analysis is based on the recently derived factorization of the direct CI coupling coefficients into internal and external parts. Using the internal coupling coefficients the integrals are transformed to new entities which are used in the diagonalization step. This two-step procedure differs significantly from the usual straightforward one-step direct CI procedure. A number of operations analysis shows that calculations with the present method should be feasible even with a large number of reference configurations in the reference state. Based on first-order perturbation theory the accuracy of the method is predicted to be close to the accuracy obtained with the usual CI method with many reference configurations.     

A method to fast determine the coupling coefficients in CI calculation

A new algorithm for evaluating the coupling coefficients and the addresses of molecular integrals in configuration interaction (CI) calculations is presented, which leads to an improved CI calculation program CGUGA. The validity and efficiency of the new code are compared with other programs, such as MELD and GAUSSIAN-94.     

Beryllium dimer,a critical test case of MBPT and CI methods

The ground-state potential curve for the beryllium dimer is calculated as a critical test case for methods based on many-body perturbation theory (MBPT ) and configuration interaction (CI ). In particular, the recently proposed double excitation (DE ) MBPT method is compared to the standard SCF-CI method including single and double excitations from a single reference determinant. The SCF-CI method is shown to give surprisingly accurate results compared to more complete CI calculations including a larger configuration space, whereas the DE-MBPT method breaks down more or less completely, particularly for larger basis sets. The results thus demonstrate the importance of including the renormalization terms in this case. Finally, Davidson's correction and related methods lead to an even more severe breakdown than the DE-MBPT method.     

Many-Electron,multicenter integrals in superposition of correlated configurations method. I. one- and two-electron integrals

The calculations by means of the superposition of correlated configurations method (Hylleraas-CI ), that is, the combination of configuration interaction with the Hylleraas-type correlation factors, needs the effective evaluation of some nontrivial integrals. This series of papers gives the formulas for all types of integrals needed for molecular calculations when Gaussian lobe functions are used as a basis set. The formulas for two-electron integrals are given in the present paper. The preliminary results for two-electron systems are presented.     

A Hilbert space multi-reference coupled-cluster study of the H4 model system

Summary Employing the Hilbert space ansatz, a fully quadratic coupled-cluster method with a multidimensional reference space is applied to a DZP basis study of the model system, H₄. The reference space is described by two to four configurations at the level of single and double excitations, and single and double excitation operators are included in the expansions for the cluster and wave operator through quadratic terms. The performance of quadratic MRCCSD is investigated for the ground and three excited states of the H₄ system consisting of two stretched hydrogen molecules in a trapezoidal configuration where the degree of quasidegeneracy is varied from a nondegenerate situation to a completely degenerate one. Compared to full CI, in the highly degenerate region, the MRCCSD works quite well. In less degenerate regions, the accuracy is less satisfactory.     

The structure of linear SiCC: An ab initio SCF CI study including vibrational effects

Large-scale SCF CI calculations have been performed for the ground ¹Σ⁺ state of linear SiCC. The calculation includes up to quadruple excitations in the valence space plus all single and double excitations from the valence localized orbitals of the single HF configuration. Vibrational wavefunctions have been derived from the CI potential surface. Vibrational frequencies including anharmonicity corrections are calculated together with the zero-point vibrational correction to the rotational constant. The large dipole moment, 4.62 D, should make this molecule suitable for radioastronomic searches.     

Ab Initio CI Determination of the Exchange Coupling Constant of Doubly-Bridged Nickel(II) Dimers

Ab initio DDCI2 (difference-dedicated configuration interaction) calculations are performed on the exchange coupling constant of the doubly-bridged Ni(II) complexes [Ni(en)(2)Cl](2)(2+) and [Ni(terpy)(N(3))](2)(2+), which are modeled by substituting the external ligands with ammonia groups. The variational CI space is selected on the grounds of the effective Hamiltonian theory and includes all the second-order contributions to the difference between the lowest quintet, triplet, and singlet states. Both complexes are found to be ferromagnetic, with coupling constants of 1.8 and 21.1 cm(-1), in good agreement with the experiment. A transformation of the molecular orbitals is also proposed for large systems, enabling the molecular orbital set to be significantly truncated-as well as the file of two-electron integrals and the DDCI2 space-with no loss of efficiency.     

Local ab initio methods for calculating optical band gaps in periodic systems. I. Periodic density fitted local configuration interaction singles method for polymers

We present a density fitted local configuration interaction singles (CIS) method for calculating optical band gaps in 1D-periodic systems. The method is based on the Davidson diagonalization procedure, carried out in the reciprocal space. The one-electron part of the matrix-vector products is also evaluated in the reciprocal space, where the diagonality of the Fock matrix can be exploited. The contraction of the CIS vectors with the two electron integrals is performed in the direct space in the basis of localized occupied (Wannier) and virtual (projected atomic) orbitals. The direct space approach allows to utilize the sparsity of the integrals due to the local representation and locality of the exciton. The density fitting approximation employed for the two electron integrals reduces the nominal scaling with unit cell size to O(N(4)). Test calculations on a series of prototypical systems demonstrate that the method in its present stage can be used to calculate the excitonic band gaps of polymers with up to a few dozens of atoms in the cell. The computational cost depends on the locality of the exciton, but even relatively delocalized excitons occurring in the polybiphenyl in the parallel orientation, can be routinely treated with this method.     

SPOCK.CI: a multireference spin-orbit configuration interaction method for large molecules

We present SPOCK.CI, a selecting direct multireference spin-orbit configuration interaction (MRSOCI) program based on configuration state functions. It constitutes an extension of the spin-free density functional theory/multireference configuration interaction (DFT/MRCI) code by Grimme and Waletzke [J. Chem. Phys. 111, 5645 (1999)] and includes spin-orbit interaction on the same footing with electron correlation. Key features of SPOCK.CI are a fast determination of coupling coefficients between configuration state functions, the use of a nonempirical effective one-electron spin-orbit atomic mean-field Hamiltonian, the application of a resolution-of-the-identity approximation to computationally expensive spin-free four-index integrals, and the use of an efficient multiroot Davidson diagonalization scheme for the complex Hamiltonian matrix. SPOCK.CI can be run either in ab initio mode or as semiempirical procedure combined with density functional theory (DFT/MRSOCI). The application of these techniques and approximations makes it possible to compute spin-dependent properties of large molecules in ground and electronically excited states efficiently and with high confidence. Second-order properties such as phosphorescence rates are known to converge very slowly when evaluated perturbationally by sum-over-state approaches. We have investigated the performance of SPOCK.CI on these properties in three case studies on 4H-pyran-4-thione, dithiosuccinimide, and free-base porphin. In particular, we have studied the dependence of the computed phosphorescence lifetimes on various technical parameters of the MRSOCI wave function such as the size of the configuration space, selection of single excitations, diagonalization thresholds, etc. The results are compared to the outcome of extensive quasidegenerate perturbation theory (QDPT) calculations as well as experiment. In all three cases, the MRSOCI approach is found to be superior to the QDPT expansion and yields results in very good agreement with experimental findings. For molecules up to the size of free-base porphin, MRSOCI calculations can easily be run on a single-processor personal computer. Total CPU times for the evaluation of the electronic excitation spectrum and the phosphorescence lifetime of this molecule are below 40 h.     

Ab-initio calculations on large molecules and solids by desirable computational procedures

Desirable computational procedures developed here recently for ab-initio calculations on large molecules are outlined. Effective core model potentials (MODPOT) permit calculations of valence electrons only explicitly, yet accurately; a charge-conserving integral prescreening evaluation to decide whether a block of integrals will be larger than a preset threshold and thus be calculated explicitly is effective for spatially extended systems; an efficient MERGE technique to save and reuse common invariant skeletal integrals is useful for geometry variations and for adding basis fcuntions, substituent groups and molecules; and an effective configuration interaction (CI) Hamiltonian into which are folded the effects of the occupied molecular orbitals from which no excitations are allowed is useful for molecular decompositions and intermolecular reactions. These techniques have been extended for CI calculations on breaking a chemical bond in a molecule in a crystal or solid; atom-class/atomic-class potential functions and dispersion calculations have been added. In a new program, POLY-CRYST, all the integral strategies for large molecules are meshed.     

Coupled-cluster method tailored by configuration interaction

A method is presented which combines coupled cluster (CC) and configuration interaction (CI) to describe accurately potential-energy surfaces (PESs). We use the cluster amplitudes extracted from the complete active space CI calculation to manipulate nondynamic correlation to tailor a single reference CC theory (TCC). The dynamic correlation is then incorporated through the framework of the CC method. We illustrate the method by describing the PESs for HF, H2O, and N2 molecules which involve single, double, and triple bond-breaking processes. To the dissociation limit, this approach yields far more accurate PESs than those obtained from the conventional CC method and the additional computational cost is negligible compared with the CC calculation steps. We anticipate that TCC offers an effective and generally applicable approach for many problems.     

Configuration selection as a route towards efficient vibrational configuration interaction calculations

A configuration selective vibrational configuration interaction (CI) approach is presented that efficiently reduces the variational space and thus leads to significant speedups in comparison to standard vibrational CI implementations. Deviations with respect to reference calculations are well below the accuracy of the underlying electronic structure calculations for the potential and hence are essentially negligible. Parallel implementations of the presented configuration selective vibrational CI approaches lead to further significant time savings. Benchmark calculations based on potential energy surfaces of coupled-cluster quality are presented for the fundamental modes of cis- and trans-difluoroethylene. The size-consistency error within the vibrational configuration interaction calculations of the difluoroethylene dimer has been studied in dependence on the excitation level.     

Efficient methods for configuration interaction calculations

Advanced techniques are developed to provide efficient economic treatment of the large scale eigenvalue problem posed when configuration interaction is carried out on SCF basis sets of moderate size. When the characteristic properties of the hamiltonian matrix are examined in light of the type of solution required, partitioning of the configuration space is shown to result in an expansion of the problem about a limited core of states, where the small but cumulative interactions of vast regions of the remaining space are reduced to the form of an effective potential. With proper selection of the core, the evaluation of this potential can be readily and accurately truncated to a level involving minimum expenditure in time and effort. In particular only diagonal elements and a strip of the full CI matrix are required to achieve an accuracy of 1 – 5 kcal/mole with complete treatment for configuration spaces of order tens of thousands. In addition, a close look at current theory on the generation of matrix elements between spin symmetry adapted configurations leads to simplified expressions where the matrix elements are derived in the form of a weighted sum of molecular integrals in which the weighting coefficients represent the integrated value of the wavefunctions over spin coordinates. For typical cases of low multiplicity and limited numbers of open shells the list of unique parameters needed to generate all weights are shown to be readily stored as a program library. Actual times for matrix element generation are believed to be an order of magnitude faster than current techniques. Practical demonstration of the accuracy and efficiency of the method is provided by calculations on formaldehyde, water, and ethylene.     

Spectral density distribution of an N-electron Hamiltonian in a finite-dimensional and spin-adapted model space

Expressions for moments of spectral density distributions of a many-electron Hamiltonian defined in a finite-dimensional, antisymmetric, and spin-adapted model space (as, e.g., a full configuration interaction space) are derived. The moments are expressed in terms of combinations of two-electron integrals corresponding to a symmetric (a two-electron singlet) and antisymmetric (two-electron triplet) two-electron wave functions. A diagrammatic approach based on Hugenholtz-type diagrams and leading to a simple and universal classification scheme of the terms appearing in the expression for a specific moment is proposed. © 1995 John Wiley & Sons, Inc.     

On the number of spin functions in the first order interaction space

A proof is given that in a configuration interaction method the first-order interaction space contains at most only twice as many spin functions as the zeroth-order space. This allows for a dramatic reduction of the size of CI expansion. For most of the high-spin systems only two spin functions for each configuration are needed.     

Die horizontale Korrelation in π-Elektronensystemen und ihre Beschreibung durch Elektronenpaarfunktionen

Different types of pair functions (geminal products and their linear combinations) are tested with respect to their ability to describe the “horizontal correlation” of the π-electrons of butadiene. The validity of the π-electron approximation is not discussed and “full configuration interaction” within the limited LCAO basis is used as the standard to which the model calculations are referred. An APSG-function (APSG = antisymmetrized product of strongly orthogonal geminals) built up from equivalent (localized) geminals, which contains only one variational parameter is able to account for about 90% of the “horizontal correlation energy”. Both APSG and APIG functions constructed from delocalized geminals, are much less favorable. Criteria of the goodness of an approximate wave function are a) the energy b) comparison of its one- and two-particle density matrices with those obtained from “full CI”. The good results with the localized APSG function are related to the fact that electron correlation between electrons of opposite spin is (in this molecule) essential only within either of the “double bonds” of the “canonical structure”. The pertinent results are quite insensitive to different parametrization of the integrals.     

Computing a new family of shape descriptors for protein structures

The large-scale 3D structure of a protein can be represented by the polygonal curve through the carbon alpha atoms of the protein backbone. We introduce an algorithm for computing the average number of times that a given configuration of crossings on such polygonal curves is seen, the average being taken over all directions in space. Hereby, we introduce a new family of global geometric measures of protein structures, which we compare with the so-called generalized Gauss integrals.     

Intracule functional models. V. Recurrence relations for two-electron integrals in position and momentum space

The approach used by Ahlrichs [Phys. Chem. Chem. Phys., 2006, 8, 3072] to derive the Obara-Saika recurrence relation (RR) for two-electron integrals over Gaussian basis functions, is used to derive an 18-term RR for six-dimensional integrals in phase space and 8-term RRs for three-dimensional integrals in position or momentum space. The 18-term RR reduces to a 5-term RR in the special cases of Dot and Posmom intracule integrals in Fourier space. We use these RRs to show explicitly how to construct Position, Momentum, Omega, Dot and Posmom intracule integrals recursively.