The spin-projected Hartree-Fock method: Direct optimization schemes and stability conditions

The necessary theory to construct different direct orbital optimization schemes in the framework of the spin-projected Hartree-Fock (spin-PHF ) method for systems with an even number of electrons and for the proper characterization of the pertinent spin-PHF solutions is developed using the formalism of second quantization. The treatment is based on a unitary transformation of a reference DODS (different orbitals for different spins) Slater determinant. The energy expectation value corresponding to the transformed DODS Slater determinant is expanded in a Taylor series about the reference DODS Slater determinant, including second-order terms. On the basis of this expansion we propose several direct iterative approaches to optimize the spin-PHF energy expectation value: for example, implementations of steepest descent and conjugated gradient techniques, and Newton–Raphson, Fletcher, and limited configuration interaction approaches. The analysis of the second-order term in the Taylor series expansion of the spin-PHF energy expectation value about a spin-PHF solution is used to define stability conditions for spin-PHF solutions. Explicit expressions are reported for all matrix elements between different spin-projected DODS Slater determinants appearing in the treatment.     

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.     

An analysis of VB structures in MO wave functions I

Summary A new method for the resolution of an antisymmetrized product of molecular orbitals into VB structures is proposed. Here VB structures are projected from a single Slater determinant associated with the ground state using the first-order density matrix. The present theory is applied to the ground state of some conjugated hydrocarbons, and special attention is paid to the covalent type VB structures.     

Computation of determinant expansion coefficients within the graphically contracted function method

Most electronic structure methods express the wavefunction as an expansion of N‐electron basis functions that are chosen to be either Slater determinants or configuration state functions. Although the expansion coefficient of a single determinant may be readily computed from configuration state function coefficients for small wavefunction expansions, traditional algorithms are impractical for systems with a large number of electrons and spatial orbitals. In this work, we describe an efficient algorithm for the evaluation of a single determinant expansion coefficient for wavefunctions expanded as a linear combination of graphically contracted functions. Each graphically contracted function has significant multiconfigurational character and depends on a relatively small number of variational parameters called arc factors. Because the graphically contracted function approach expresses the configuration state function coefficients as products of arc factors, a determinant expansion coefficient may be computed recursively more efficiently than with traditional configuration interaction methods. Although the cost of computing determinant coefficients scales exponentially with the number of spatial orbitals for traditional methods, the algorithm presented here exploits two levels of recursion and scales polynomially with system size. Hence, as demonstrated through applications to systems with hundreds of electrons and orbitals, it may readily be applied to very large systems. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Evaluation of expansion coefficients for translation of Slater‐type orbitals using complete orthonormal sets of exponential‐type functions

By the use of exponential‐type functions (ETFs) the simpler formulas for the expansion of Slater‐type orbitals (STOs) in terms of STOs at a displaced center are derived. The expansion coefficients for translation of STOs are presented by a linear combination of overlap integrals. The final results are of a simple structure and are, therefore, especially useful for machine computations of arbitrary multielectron multicenter molecular integrals over STOs that arise in the Hartree–Fock–Roothaan approximation and also in the Hylleraas correlated wave function method for the determination of arbitrary multielectron properties of atoms and molecules. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 126–129, 2001     

Evolutionary algorithm based configuration interaction approach

A new evolutionary algorithm for stochastic configuration interaction (CI) method designed as an affordable approximation to full configuration interaction (FCI) has been described here. The key components of the algorithm are initiation, propagation, and termination steps taking inspiration from the genetic algorithm. The propagation step is performed with cloning (retention of a Slater determinant without change), mutation (single excitation/de‐excitation), and crossover (exchange of α and β strings between two Slater determinants) and termination is selection of few Slater determinants based on certain fitness function (measure of importance of a determinant in the CI space) and rejection of the rest. We find that the absolute value of the CI coefficients is a suitable fitness function when combined with a fixed selection scheme. We have tested its accuracy in 1D Hubbard problem and ground state potential energy surface (PES) has also been constructed for symmetric bond breaking of water molecule, where the errors are found to be around 10 mE_h with low non‐parallelity error, when retaining only a small fraction of the total number of Slater determinants in the final population. This shows that this method has the ability to capture both static and dynamic correlation. Performance and convergence properties of the algorithm are also tested for N₂ triple bond breaking problem. The algorithm opens up a promising way for stochastic sampling of the important determinants in the full Hilbert space.     

On the theory of electron correlation in atoms and molecules. II. General cluster expansion theory and the general correlated wave functions method

An exact cluster expansion of many electron wave functions is derived, beginning with a finite linear combination of Slater determinants rather than the more usual single determinant. This general cluster expansion is found to apply both in the case where all possible Slater determinants from a finite set of spin orbitals are included in the linear combination, and in the case where the number of determinants is restricted. The special properties of that finite linear combination of determinants closest to the exact wave function in the least squares sense are studied. These properties lead to the derivation of a general correlated wave functions method, illustrating again the close relationship between methods of this type and cluster expansion theory. Additional approximations, necessary for practical calculations, are set out.     

First-order properties and the Hellmann–Feynman theorem in the case of a limited CI wave function

Two distinct approaches to the calculation of first-order properties with a limited CI wave function are discussed. One is based on the Hellmann–Feynman theorem and the other on the direct evaluation of the total energy derivative at zero perturbation. Corrections to the Hellmann–Feynman expectation value are given for the CI wave function consisting of a single determinant reference state and all single and double replacements of this. These corrections are the extended Brillouin matrix elements and involve interactions between the zeroth-order wave function and triply substituted configurations. The usefulness of these matrix elements for the generation of MC SCF orbitals and for the calculation of cluster corrections to the wave function is briefly discussed. The formulas for the Brillouin matrix elements expressed in terms of one- and two-electron integrals have been automatically generated using the syntax of the algebraic program SCHOONSCHIP.     

Density matrix in the open shell theory

All the second-order density matrix spin components for the spin-projected Slater determinant are presented as expansions in direct products of powers of unprojected spin- and residual electron density matrices. Spin components of the second-order transition density matrices between spin-projected Slater determinants built of orthogonal orbitals are also obtained.     

A sparse algorithm for the evaluation of the local energy in quantum Monte Carlo

A new algorithm is presented for the sparse representation and evaluation of Slater determinants in the quantum Monte Carlo (QMC) method. The approach, combined with the use of localized orbitals in a Slater-type orbital basis set, significantly extends the size molecule that can be treated with the QMC method. Application of the algorithm to systems containing up to 390 electrons confirms that the cost of evaluating the Slater determinant scales linearly with system size.     

Adjustment of Born‐Oppenheimer electronic wave functions to simplify close coupling calculations

Technical problems connected with use of the Born‐Oppenheimer clamped‐nuclei approximation to generate electronic wave functions, potential energy surfaces (PES), and associated properties are discussed. A computational procedure for adjusting the phases of the wave functions, as well as their order when potential crossings occur, is presented which is based on the calculation of overlaps between sets of molecular orbitals and configuration interaction eigenfunctions obtained at neighboring nuclear conformations. This approach has significant advantages for theoretical treatments describing atomic collisions and photo‐dissociation processes by means of ab initio PES, electronic transition moments, and nonadiabatic radial and rotational coupling matrix elements. It ensures that the electronic wave functions are continuous over the entire range of nuclear conformations considered, thereby greatly simplifying the process of obtaining the above quantities from the results of single‐point Born‐Oppenheimer calculations. The overlap results are also used to define a diabatic transformation of the wave functions obtained for conical intersections that greatly simplifies the computation of off‐diagonal matrix elements by eliminating the need for complex phase factors. © 2013 Wiley Periodicals, Inc.     

The generalized Slater–Condon rules

By relating the blocking structure of the relevant matrix of overlap-integrals to its cofactors, the Slater–Condon rules for the evaluation of an element of a matrix representation of an electronic Hamiltonian in a Slater determinant basis are generalized to the case where not all orbitals are orthogonal. This yields a set of 33 rules, which allows for an efficient implementation of the valence bond theory.     

On the calculation of arbitrary multielectron molecular integrals over Slater‐type orbitals using recurrence relations for overlap integrals I. Single‐center expansion method

The multicenter charge‐density expansion coefficients [I. I. Guseinov, J Mol Struct (Theochem) 417 , 117 (1997)] appearing in the molecular integrals with an arbitrary multielectron operator were calculated for extremely large quantum numbers of Slater‐type orbitals (STOs). As an example, using computer programs written for these coefficients, with the help of single‐center expansion method, some of two‐electron two‐center Coulomb and four‐center exchange electron repulsion integrals of Hartree–Fock–Roothaan (HFR) equations for molecules were also calculated. Accuracy of the results is quite high for the quantum numbers, screening constants, and location of STOs. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 146–152, 2000     

Unified treatment of electric multipole transition radial matrix elements using radial parts of STO and complete orthonormal sets of ψα-ETO

A new method is proposed for evaluating electric multipole transition radial (EMTR) matrix elements of the hydrogen-like atoms using the expansion formulae suggested by one of the authors for the radial parts of ψ^α-exponential type orbitals (ψ^α-ETO) in terms of radial parts of Slater type orbitals (STO), where α is the frictional quantum number (?∞<α≤2). A comparative study is carried out between the proposed analytical computations and other numerical simulations. The effectiveness and accuracy of the obtained analytical expressions are demonstrated by calculation of concrete cases.     

PCILON. perturbative configuration interaction using localized orbitals and numerical integration. I. Numerical integration techniques for the calculation of Hamiltonian matrix elements between localized orbitals

Modern techniques for multidimensional numerical integration, Korobov's and Sobol's formulas namely, are used for the direct computation of matrix elements between the localized molecular orbitals needed for a configuration interaction calculation by a perturbation method. A minimal orbital basis of Slater functions is used for formaldehyde and ethylene taken as example. The resulting precision is satisfactory.     

Projected BCS Tamm-Dancoff method for molecular electronic structures

A new Tamm–Dancoff method for the ground and excited states of molecular electronic systems is developed. The method begins with a number-projected BCS (PBCS ) wave function and is generated by excitations of particle pairs from the degenerate geminals in the PBCS wave function. A direct optimization of the PBCS wave function is accomplished with successive Bogoliubov transformations so that one-pair excitation terms in the Tamm-Dancoff expansion of the ground state vanish (the generalized Brillouin theorem). The spin-symmetry adapted first- and second-order Tamm–Dancoff bases and matrix elements are calculated by means of the CI expansion of the PBCS wave function with natural orbitals that diagonalize the BCS geminal matrix. Numerical calculations are presented for the H₄ system with D_2h and D_4h conformations and for methylene. The PBCS wave function is not a very good approximation for the ground state, accounting for only about half of the correlation energy. The second-order Tamm–Dancoff correction improves the result as much as the double excitation CI . The Tamm–Dancoff terms consisting of two triplet pairs coupled to a singlet, and those relaxing the constraint imposed on the pairwise excitations in the PBCS wave function are important.     

块定域波函数方法及其应用

提出了块定域波函数方法以定量分析分子内的电子定域现象或分子间的电荷传递效应。对于一个假想的严格定域的分子,我们通过将全部的电子和基轨道配分成几个子空间来构造其相应的波函数。其中每一个分子轨道只对某一个子空间展开,各子空间内的分子轨道相互正交,但不同子空间内的分子轨道间是非正交的。Hartree-Fock波函数和块定域波函数之间的能量之差即为分子内的电子定域能或分子间的电荷传递能。我们应用块定域波函数方法讨论了丁二烯分子中的旋转势垒。     

Direct MRCI method for the calculation of relativistic many-electron wavefunctions. I. General formalism

The formalism of a quasi- or full-relativistic multireference CI method has been developed and implemented. The scheme is appropriate for the calculation of molecular systems in which the relativistic effects are of the same order of magnitude as the correlation contributions. In this contribution some important symmetry aspects of a relativistic many-electron wave function are discussed and the consequences for the CI matrix structure are shown. An efficient CI strategy in the form of a direct CI is presented, which avoids the construction of the whole CI matrix. Based on a determinantal expansion of molecular spinor products, the individual one- and two-electron molecular integrals are processed, and the molecular symmetry is easily accounted for by a proper linear combination of Slater determinants in the CI starting vector. For an efficient CI organization some powerful techniques of the graphical unitary group approach have been transferred to the relativistic case.