Electronic correlation studies. II. Self-correlated field method. Application to ground state and first 1P, 3P excited states of two-electron atomic systems

An antisymmetric pair function can be built upon two kinds of monoelectronic functions, the former ones being correlated local functions and the second ones nonlocal functions taking external effects into account. This function, brought into the generalised product function procedure by means of the density matrix formulation, makes possible the study of correlation within N-electronic systems. The results of a first application of this method to the fundamental and to ¹P and ³P excited states of two-electron systems are given.     

Electron correlation study for two-electron atoms by a simple correlated wave function

The difference in the electron correlation between H^- and other two-electron atoms is clarified by the introduction of the r₁₂ term in the wave function. By using the expansion of r₁₂, a certain modification of the usual electron correlation factor 1 + Cr₁₂ is introduced and its effectiveness is examined. Calculations are carried out for the ground state and the three lowest excited states (2³S, 2³P and 2¹P). The peculiar electron correlation in the ground state of H^? is shown by looking at the Coulomb hole for closed- and open-shelf models in comparison with those for other two-electron atoms.     

Analytical gradients for Singer's multicenter n‐electron explicitly correlated Gaussians

Analytical gradients for Singer's basis of n‐electron multicenter explicitly correlated Gaussian functions are derived and implemented to variationally optimize the energy and wave function of molecular systems within the Born–Oppenheimer approximation. Wave functions are optimized with respect to (½n(n+1)+3n) nonlinear variational parameters and one linear coefficient per term in the basis set. Preliminary results for the ground states of H₃⁺ and H₃ suggest that the method can be more flexible and can achieve lower energies than previously reported calculations. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 82: 151–159, 2001     

Half-Projected and Projected Hartree-Fock Calculations for Singlet Ground States. i. four-Electron Atomic Systems

The half-projected Hartree-Fock function (HPHF ) for singlet states is defined as a linear combination of two Slater determinants, which contains only spin eigenfunctions with even quantum number. Using a self-consistent procedure based on the generalized Brillouin's theorem, the RHF , HPHF and PHF functions are deduced for the ground states of the Li^?, Be, B⁺, and C²⁺ systems, in a limited basis set. It is found that the HPHF function yields better energy values than the RHF function, very close to that of the PHF one. The HPHF scheme seems thus to be useful as a substitute for the PHF model, specially in the case of large electronic systems in which the latter method becomes unmanageable.     

Numerical investigation of a lower bound to the Thomas—Fermi kinetic energy and the validity of the Lieb conjecture for H−, He,Li+, B3+, O6+, Ne8+, and Mg12+

The purpose of the present work was twofold: First, a lower bound of Gálvez and Dehesa to the Thomas-Fermi kinetic energy was investigated for the two-electron systems H⁻, He, Li⁺, B³⁺, O⁶⁺, Ne⁸⁺, and Mg¹²⁺. Second, a conjecture of Lieb, relating the kinetic energy to the Thomas-Fermi kinetic energy, was examined. For both investigations, the analytical approximations of Benesch to the radial electron densities of the two-electron systems were used. These approximations are based on the explicitly correlated (Hylleraas-type) 20 variational parameter wave functions of Hart and Herzberg. © 1996 John Wiley & Sons, Inc.     

Calculations of autoionization states of He and H−

The positions and lifetimes of several ¹S and ^1,3P autoionizing states of He and H^? are obtained by two methods involving standard techniques of electronic structure calculation which can be extended to more complicated systems. The first method involves an approximate evaluation of Miller's “golden rule” formula; the second is an application of the recently developed complex coordinate method.     

Physical nature of the chemical bond II. Valence atomic orbital and energy partitioning studies of linear nitriles

The valence atomic orbitals (VAO 's) of several linear nitriles are determined using non-empirical SCF –LCAO –MO wave functions expanded in a minimal (CN^?, HCN, FCN, C₂N₂), double-zeta (CN^?, HCN), or double-zeta + polarization (HCN) basis of Slater atomic orbitals (AO 's). The molecular energy of each system (except the double-zeta + polarization HCN system) is partitioned according to the procedure of Ruedenberg to obtain numerical values of nitrile C and N atomic and C?N bond components of the energy. In addition, the nitrile results are compared with minimal AO basis results obtained previously by other authors for homonuclear diatomics, diatomic hydrides and H₂O. The numerical data are used to test the internal self-consistency of the various definitions entering the partitioning method, i.e. whether or not analogous quantities assume similar values in chemically similar situations. The analysis of nitrile SCF –MO wave functions in terms of the set of VAO 's characteristic of the system under consideration is shown to be a promising approach to the problem of extracting useful information from the wave functions. In general, numerical results for the nitrile systems studied are fairly consistent with the concepts on which the partitioning method is based: promotion, quasi-classical interaction, sharing penetration, sharing interference and charge transfer. However, the VAO expansions for several energy components need to be investigated further and possibly revised.     

Periodic orbits in atomic hydrogen exposed to circularly polarised laser light

Exact circular periodic orbits for a hydrogen atom in a strong circularly polarised laser field are derived. A stability analysis shows one orbit to be stable or unstable depending on the laser parameters, whereas the other orbit is always unstable. These orbits are expected to manifest themselves as resonances in laser assistede ^? ?H ⁺ scattering. The binding mechanism for these resonances is provided by the trapping of the quantum mechanical wave function onto the neighbourhood of a classical periodic orbit. This mechanism is analogous to that for a Wannier ridge resonance in a doubly excited two-electron atom.     

Factored wave function for bound S-type states of two-electron atomic systems

A simple and accurate variational wave function in which the dependence in the interelectronic distance is factored is proposed to describe S-type states of two-electron atomic systems. We introduce a parameterization which generalizes the previous ones used in this same framework and which allows us to obtain in a simple way the wave function of both symmetric and antisymmetric excited states. We performed a systematic analysis of some exact properties such as the virial theorem and the cusp conditions and a study of both the one- and two-body densities. Finally, a comparison among the different correlation functions for these states was performed for helium. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 68: 405–413, 1998     

Ground state for two‐electron and electron‐muon three‐body atomic systems

In this article, the angular correlated configuration interaction method previously introduced by some of the authors is extended to three‐body atomic systems with general masses. A recently proposed angularly correlated basis set is used to construct ground state wave functions which: (i) satisfy exactly Kato cusp conditions at the two‐body coalescence points; (ii) have only linear coefficients; and (iii) show a fast convergency rate for the energy. The efficiency of the construction is illustrated by the study of the negatively charged hydrogen‐like systems (^∞H^?, T^?, D^?, ¹H^?, and Mu^?), neutral helium‐like systems (e^?e^{? ∞}He⁺²,e^?e^{? 4}He⁺², e^?e^{? 3}He⁺², e^?μ^{? ∞}He⁺², e^?μ ^?4He⁺², and e^?μ^{? 3}He⁺²), and positively charged lithium‐like systems (e^?e^{? ∞}Li⁺³, e^?e^{? 7}Li⁺³, e^?e^{? 6}Li⁺³, e^?μ^{? ∞}Li⁺³, e^?μ^{? 7}Li⁺³, and e^?μ^{? 6}Li⁺³). The ground state energies and other mean values are compared with those given in the literature, when available. Wave functions with a moderate number of (20 maximum) linear coefficients are given explicitly; they are sufficiently simple and accurate to be used in practical calculations of atomic collision in which multidimensional integrations are involved. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Analytical potential from ab initio calculations for the Fe+-H2O and Feo-H2O systems

Analytical intermolecular potentials for the Fe⁺?H₂O and Fe^o?H₂O systems have been determined from ab initio calculations. Interaction energies for a lot of points along the two potential energy surfaces were calculated using Huzinga's MINI ?2 basis set. The results obtained were fitted to an analytical function containing 11 adjustable parameters that we have already used with success for the Fe²⁺?H₂O system. The goodness of the generated intermolecular potentials is discussed.     

A Green's function approach to the nonrelativistic radial wave equation of hydrogen atom

We report an efficient and consistent method for the solution of Schrödinger wave equation using the Green's function technique and its successful application to the nonrelativistic radial wave equation of hydrogen atom. For the radial wave equation, the Green's function is worked out analytically by means of Laplace transform method and the wave function is proposed under the boundary conditions. Computationally the product of potential term, the proposed wave function and the Green's function are integrated iteratively to get the nonrelativistic radial wave function. The resultant wave after each iteration is normalized and plotted against the standard nonrelativistic radial wave function. The solution converges to the standard wave with the increasing number of iterations. Results are verified for the first 15 states of hydrogen atom. The method adopted here can be extended to many‐body problem and hope that it can enhance our knowledge about complex systems.     

A three-dimensional finite element approach towards molecular SCF computations

Fully three-dimensional, SCF ground-state computations for the Hartree equation are carried out by a finite element approach that completely avoids forming or storing the Fock matrix. A combination of strategies is used to reduce storage and computational requirements by orders of magnitude over the traditional finite element approach, which makes three-dimensional molecular orbital computations feasible. Results using the three-dimensional formulation and computer program are shown for one-electron systems: He⁺ and H₂⁺, and for two-electron systems: He and H₂. The best results are within about 30–100 micro-Hartrees of the exact values of the total energies for the ground states of these systems, indicating that our three-dimensional approach has been correctly implemented in the computer code.     

A semiempirical method based on geminal functions

An attempt has been made to develop a semiempirical method which considers only the n- and π-electrons, with the eigenfunctions expressed as an antisymmetrized product of two-electron functions or geminals. These geminals are expressed as a linear combination of products of Hückel-type MO's and the matrix elements are evaluated assuming the strong orthogonality condition to hold among the geminals, with an average effective Hamiltonian where the interaction between paired electrons is explicitly included. A first application of the method to N₂, HCN, C₂H₂, and C₂H^? was carried out without the introduction of parameters adjustable to best fit. A parametric approximation was then used for studying the electronic structure of diazabarrelene, an unknown molecule containing three pairs of π-electrons and two pairs of n-electrons. It is concluded that the explicit introduction of (π,π) and (n-π)-electron interaction describes the ground state more realistically than the simplest Hückel method. The separation of the energy levels was also affected, but the calculated transitions compared rather poorly with the spectroscopic observations. Diazabarrelene is expected to be not less stable than barrelene and attempts to its synthesis should be considered worthy of being made.     

Note on the space part of anti-symmetric wave functions in the many-electron problem

It is pointed out that if a many-electron antisymmetric wave function is expanded as a sum of spin-product functions, each multiplied by a function of coordinates, the resulting functions of coordinates have many of the same useful features found with the symmetric and antisymmetric functions representing singlet and triplet states in a two-electron system. For finding the energy, or any function of coordinates only, in the approximation in which spin-orbit interaction is neglected, one such function of coordinates can be used, the spins being disregarded. Simple procedures allow one to find matrix components of such operators as S ² and L . S from the functions of coordinates. These procedures are much easier to visualize than the use of projection operators, the permutation group, or other methods in current use. The general procedures are illustrated by application to the three-electron problem of the lithium atom, as treated by Lunell, Kaldor, and Harris, and their application to the contact hyperfine structure is pointed out.     

On the SCF calculation of excited states: Singlet states in the two-electron problem

The problem of determining SCF wave functions for excited electronic states is examined for singlet states of two-electron systems using a Lowdin natural orbital transformation of the full CI wave function. This analysis facilitates the comparison of various SCF methods with one another. The distribution of the full CI states among the natural orbital MCSCF states is obtained for the S states of helium using a modest Gaussian basis set. For SCF methods that are not equivalent to the full CI wave functions, it is shown that the Hartree-Fock plus all single excitation wave functions are equivalent to that of Hartree-Fock plus one single excitation. It is further shown that these wave functions are equivalent to the perfect pair or TCSCF wave functions in which the CI expansion coefficients are restricted to have opposite signs. The case of the natural orbital MCSCF wave function for two orbitals is examined in greater detail. It is shown that the first excited state must always be found on the lower natural orbital MCSCF CI root, thus precluding the use of the Hylleras-Undeim-MacDonald (HUM) theorem in locating this state. It is finally demonstrated that the solution obtained by applying the HUM theorem (minimizing the upper MCSCF CI root with respect to orbital mixing parameters) is an artifact of the MCSCF method and does not correspond to any of the full CI states.     

Simple extended STO basis sets. Helium

The two-parameter function, φ = (C₁ + C₂r^n?1) exp (?ζr), (n = 2–5), has been used as a basis function to determine the independent particle model energy of two-electron atomic systems in their ground state. The best energy is found for n = 3 (He—B³⁺) and for n = 4 (H^?). Our energy values are significantly close to Hartree-Fock results.     

Some properties of the bivariational Hartree–Fock scheme for complex symmetric many-particle operators

The bivariational Hartree–Fock scheme for a general many-body operator T is discussed with particular reference to the complex symmetric case: T^? = T*. It shown that, even in the case when the complex symmetric operator T is real and hence also self-adjoint, the complex symmetric Hartree–Fock scheme does not reduce to the conventional real form, unless one introduces the constraint that the N-dimensional space spanned by the Hartree–Fock functions ? should be stable under complex conjugation, so that ?* = ?α. If one omits this constraint, one gets a complex symmetric formulation of the Hartree–Fock scheme for a real N-electron Hamiltonian having the properties H = H* = H^?, in which the effective Hamiltonian H_eff (1) may have complex eigenvalues ε_k. By using the method of complex scaling, it is indicated that these complex eigenvalues—at least for certain systems—may be related to the existence of so-called physical resonance states, and a simple example is given. Full details will be given elsewhere.