A stationary property of the APSG wave function

The method of antisymmetrized product of strongly orthogonal geminals (APSG) emerges by optimizing the one-electron functions used to construct two-electron functions (geminals), the latter being expanded (with coefficients selected variationally) in mutually exclusive subspaces of the former ones. Accordingly, E _APSG is stationary with respect to the expansion coefficients and to unitary transformations of the one-electron orbitals. We show that the APSG energy is also stationary to the unitary transformation of identical geminals. For non-identical geminals, this statement holds only approximately.     

Semiempirical implementation of strictly localized geminals for analysis of molecular electronic structure

Approximate electronic trial wave function taken as the antisymmetrized product of strictly localized geminals (APSLG) is implemented for semiempirical analysis of molecular electronic structure of “organic” compounds and for calculations of their heats of formation. This resulted in an O(N)‐scaling method. Using the MINDO/3 form of the semiempirical Hamiltonian with reparameterized β_AB values in combination with the APSLG form of the wave function yields the computational procedure BF'98. Calculations on the heats formation and the equilibrium geometries for a wide range of molecules show that the APSLG‐MINDO/3 approach is more favorable than its self‐consistent field‐based counterpart. Also, the APSLG formalism allows to interpret molecular electronic wave function in chemically sensible terms. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 752–764, 2001     

Bare Coulomb field and atomic reciprocal form factor

The reciprocal form factor of N‐electron closed shells systems in a bare Coulomb field is shown to be a spherically symmetric, positive, and decreasing function of the radial distance. Nonmonotonicities of the reciprocal form factor appear when studying bare Coulomb field open‐shell systems. Analysis of the weight of the interelectronic repulsion term is carried out for some isoelectronic series as well as neutral atoms with N = 1–103. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

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.     

Applications of n-electron wave functions built up with two kinds of geminals

The ground states of atoms and molecules Li^?, Be, LiH, LiH, and Li₂ have been calculated using the n-electron wave functions built up with two kinds of geminals. As a comparison, the above systems have been calculated with the Hartree–Fock self-consistent field and the multi-configuration self-consistent field method as well. The results show that the wave functions in this work are capable of describing the electron correlations, and both kinds of geminals can be taken as a starting point in building up n-electron ground states.     

Spin adaptation of antisymmetrized geminal product wave functions

A spin-free polynomial representation of antisymmetrized geminal products is presented for several cases. In particular, products of identical geminals, which possess different spin multiplicity, are considered. The cases of singlet geminals, singlet geminals with one or two triplet geminals coupled to the lowest possible spin multiplet, and triplet geminals coupled to an arbitrary multiplet are considered in detail, and explicit polynomial representation is given.     

On the behaviour of the wave function for the 23S state of the two-electron atom in the region where both electron–nuclear distances are small

The expansion of the wave function for the 2³S state of the two-electron atom in the neighbourhood of the singularity at r₁ = r₂ = 0 is considered. The restrictions imposed on the variational functions by this expansion are discussed. For the 2³S state of He, Li⁺, N⁵⁺ the behaviour of the variational function based on the Fock expansion in the neighbourhood of this singularity is investigated. The agreement of the variational coefficients with the theoretical coefficients is satisfactory. The calculated values of E and 〈δ(r₂)〉 for He, Li⁺, N⁵⁺ are given.     

The interaction of chemical bonds. III. Perturbed strictly localized geminals in LMO basis

The formalism of strictly localized geminals (SLGs ) is summarized. It is shown that the SLG wave function serves as an appropriate multiconfigurational reference state that can easily be improved by perturbational, CI - or coupled cluster-type procedures. The possibility of expanding the geminals in the basis set of localized Hartree-Fock molecular orbitals (LMOs ) is discussed. Sample calculations on H₄, CH₄, H₂O, and He…?He systems are reported. © 1994 John Wiley & Sons, Inc.     

Supersymmetric solutions of <Emphasis Type="Italic">PT</Emphasis>-/non-<Emphasis Type="Italic">PT</Emphasis>-symmetric and non-Hermitian screened Coulomb potential via Hamiltonian hierarchy inspired variational method

The supersymmetric solutions of PT -symmetric and Hermitian/non-Hermitian forms of quantum systems are obtained by solving the Schr?dinger equation for the Exponential-Cosine Screened Coulomb potential. The Hamiltonian hierarchy inspired variational method is used to obtain the approximate energy eigenvalues and corresponding wave functions.      

Geminal model chemistry II. Perturbative corrections

We introduce and investigate a chemical model based on perturbative corrections to the product of singlet-type strongly orthogonal geminals wave function. Two specific points are addressed (i) Overall chemical accuracy of such a model with perturbative corrections at a leading order; (ii) Quality of strong orthogonality approximation of geminals in diverse chemical systems. We use the Epstein-Nesbet form of perturbation theory and show that its known shortcomings disappear when it is used with the reference Hamiltonian based on strongly orthogonal geminals. Application of this model to various chemical systems reveals that strongly orthogonal geminals are well suited for chemical models, with dispersion interactions between the geminals being the dominant effect missing in the reference wave functions.     

Quasi-relativistic approximation in the SCF-MO-LCAO method

A quasi-relativistic approach to the MO-LCAO method is formulated taking into account the relativistic effects with an accuracy up to (v/c)² terms, the relativistic part of the electronic interaction in the Hamiltonian being neglected. In the framework of this approximation a set of SCF equations of the Roothaan form is derived; here only the relativistic analogue to the closed shell systems with one-determinant wave functions is considered. In so doing three types of relativistic corrections arise which are quite similar to those of the Pauli equation for one-electron atoms. The new matrix elements appearing due to these corrections can be reduced to some common integrals, which have to be calculated with relativistic radial atomic functions. The method allows a semi-empirical approach to the problem and does not require the Dirac four-component atomic functions (unknown in the most cases), thus making possible approximate quasi-relativistic electronic structure calculations of heavy-atom compounds.     

Equation for the direct determination of the density matrix: Time-dependent density equation and perturbation theory

The density equation proposed previously for the direct determination of the density matrix, i.e. for the wave mechanics without wave, is extended to a time-dependent case. The time-dependent density equation has been shown to be equivalent to the time-dependent Schr?dinger equation so long as the density matrix, included as a self-contained variable, is N-representable. Formally, it is obtainable from the previous time-independent equation by replacing the energy E with . The perturbation theory formulas for the density equation have also been given for both the time-dependent and time-independent cases. Received: 16 June 1998 / Accepted: 2 September 1998 / Published online: 8 February 1999     

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.     

Product of group-function expansions for correlated wave functions

A theorem is proved which demonstrates the relationship between a product of group functions describing the correlated motion of a particular group of electrons in an N-electron system and a wave function obtained from the exact wave function which describes the correlation of the same group of electrons. By considering such products of group functions as elements in a variational wave function, an expansion for correlated wave functions is suggested, which emphasizes the correlated motion of groups of electrons in the whole system.     

A new look at correlations in atomic and molecular systems. I. Application of fermion monte carlo variational method

We are engaged in research directed toward the development of compact and accurate correlation functions for many-electron systems. Our computational tool is the variational method in which the many-electron integrals are calculated by Monte Carlo using the fermion Metropolis sampling algorithm. That is, a many-fermion system is simulated by sampling the square of a correlated antisymmetric wave function. The principal advantage of the method is that interelectronic distance r_ij may be included directly in the wave function without adding significant computational complexity. In addition, other quantities of physical and theoretical interest such as electron correlation functions and representations of Coulomb and Fermi “holes” are very easily obtained. Preliminary results are reported for He, H₂, and Li₂.     

The virial theorem scaling model for estimating the charge sensitivities of hydrogens in molecules

A simple variational model of hydrogens in molecules is presented, using the virial theorem (Fock) scaling of the wave function to account for the orbital relaxation due to a change in the number of electrons. The resulting interpolative formulas for the energy allow realistic predictions of the spin-nonpolarized or spin-polarized hardness (electron repulsion) parameters and other sensitivities of the hydrogen systems in molecules, including the realistic bare nuclei limit (N → 0) data.     

Shell model calculations: An efficient new algorithm

We describe an efficient new algorithm which extends the range of feasible shell model calculations. This algorithm is applicable to single shell and multiple shell configurations, where two or more quantum numbers (e.g., L and S) are required to label the states within each shell. The algorithm proceeds by factoring the shell model Hilbert space into a product of subspaces, one for each angular momentum. N-particle wave functions are built up recursively from N – 1 particle wave functions. Three kinds of N – 1- to N-particle coefficients are required to carry out the construction of N-particle electron (or fermion) states from N – 1 particle states. These are (1) coefficients of fractional parentage (CFP s) within a single shell, (2) outerproduct isoscalar factors (OISF s) within a single angular momentum subspace, and (3) innerproduct isoscalar factors (IISF s) which describe how multishell states within the complementary angular momentum subspaces are combined to form totally antisymmetric wave functions. All three types of N – 1- to N-particle coefficients are generated recursively using a single powerful and efficient matrix diagonalization algorithm. Matrix elements of single particle creation and annihilation operators are expressed in terms of single particle CFP s, OISF s, and IISF s. We also describe an efficient algorithm for computing matrix elements of products of creation and anihilation operators by inserting and summing over complete sets of intermediate states. This is the Feynman-like sum over path overlaps procedure. Timing benchmarks are presented comparing the new Drexel University shell model (DUSM ) code with a state of the art shell model code.     

Approximate ansatz for the expansion of the spherically averaged wave function in terms of interelectronic separation r12 for the Hookean atom,atomic ions,and the H2 molecule

For the two‐electron Hookean atom, it is first emphasized that, for a specific force constant k = 1/4, the ground‐state wave function has a simple dependence on the interelectronic separation r₁₂, namely, (1 + ½r₁₂)exp(??r). For this two‐electron model, therefore, the study of Rassolov and Chipman on the electron–electron cusp conditions on the spherically averaged wave function for the N electron atomic ions can be generalized to all orders in the interelectronic separation r₁₂. This Hookean model has therefore been used to give some justification for an ansatz for the spherically averaged wave function in atomic ions with N electrons for N ≥ 2. Several approximate two‐electron wave functions satisfying the Rassolov and Chipman conditions were tested and found to give excellent results. Another ansatz has been tested numerically on the ground state of two‐electron atomic ions and the H₂ molecule. Finally, for the Hookean atom a partial differential equation that is essentially for the pair correlation density is given in the Appendix . © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 95: 21–29, 2003     

Variational Monte Carlo calculations for some cations and anions of the first‐row atoms using explicitly correlated wave functions

The variational Monte Carlo method is applied to calculate ground‐state energies of some cations and anions of the first‐row atoms. Accurate values providing between 80 and 90% of the correlation energy are obtained. Explicitly correlated wave functions including up to 42 variational parameters are used. The nondynamic correlation due to the 2s ? 2p near degeneracy effect is included by using a multideterminant wave function. The variational free parameters have been fixed by minimizing the energy that has shown to be a more convenient functional than the variance of the local energy, which is the most commonly employed method in variational Monte Carlo calculations. The energies obtained improve previous works using similar wave functions. © 2002 Wiley Periodicals, Inc.; DOI 10.1002/qua.10125