The ellipsoidal gaussian basis in molecular orbital theory

New analytic integral formulas are presented for the potential energy integrals over ellipsoidal Gaussian basis functions [ exp (-x² - y² - z²)] that enter into solving the conventional expansion self-consistent field equations. Near minimal atomic orbital bases combined from large nuclear-centered primitive Gaussian sets are used in test calculations on the HF and CO molecules. The ellipsoidal exponential parameters for the valence atomic orbitals are fully optimized using a single scale factor for each atomic orbital and nuclear coordinate. The results are compared with those obtained using an unoptimized nuclear centered double-zeta spherical Gaussian basis.     

On the mathematical formulation and interpretation of LACO MO theory

The mathematical basis of LCAO MO theory is studied, both within the Hartree-Fock approximation and in more exact formulations. The basic LCAO expansion for molecular orbitals ¦> in terms of atomic orbitals ¦x> is conveniently written ¦> = ¦x> S ^–1 B where S is the overlap matrix for atomic orbitals and B is the matrix of atomic orbital-molecular orbital overlaps. It is suggested that matrices P and Q, defined by P=B B and Q=BnB where n is the matrix of molecular orbital occupation numbers, are appropriate to the interpretation of molecular calculations in terms of atomic orbital components, electronic populations and the degree of bonding. Implications for Hartree-Fock calculations are investigated.     

A novel interpretation of reduced density matrix and cumulant for electronic structure theories

We propose a novel interpretation of the reduced density matrix (RDM) and its cumulant that combines linear and exponential parametrizations of the wavefunction. Any n-particle RDM can be written as a weighted average of "configuration interaction" amplitudes. The corresponding n-particle cumulant is represented in terms of two types of contributions: "connected" (statistical averages of substitution amplitudes) and "disconnected" (cross-correlations of substitution amplitudes). A diagonal element of n-RDM represents the average occupation number of the orbital n-tuple. The diagonal elements of 2- and 3-cumulants take particularly elegant forms in the natural spin-orbital basis: they represent the covariances (correlated fluctuations) of the occupation numbers of the orbital pair and triples, respectively. Thus, the diagonal elements of the cumulants quantify the correlation between the orbital occupation numbers. Our interpretation is used to examine the weak to strong correlation transition in the "two electrons in two orbitals" problem.     

Equations de Hartree–Fock dependant du temps

The time-dependent Hartree–Fock (TDHF ) equations are derived up to the second order when the system is perturbed by a monochromatic plane wave. The solutions of the equations are subjected to the orthonormalization conditions satisfied by the orbitals. In the equations, these conditions are expressed by the appearing of coefficients λ playing the part of Lagrangian multipliers. Relations between the coefficients λ are established. These relations are equivalent to the above-mentioned orthonormalization conditions. This equivalence enables us to substitute for the solution of an integrodifferential equation system subject to constraint conditions, that of a free system. The TDHF equations obtained determine the first- and second-order orbital perturbations, which no doubt verify the orthonormalization conditions. These orbitals can be used in the calculation, up to second order, of different nonlinear optical effects.     

Natural spin orbital analysis of diatomic molecular wave functions in terms of generalized diatomic orbitals I. Outline of the method. Results for the ground state of H2

The method of linear combinations of generalized diatomic orbitals (LCGDO) is combined with the method of configuration interaction (CI). CI wave functions obtained in this way are finally submitted to a natural spin orbital analysis; the resulting natural spin orbitals are expansions in terms of generalized diatomic orbitals.For the ground state of H₂, a one-determinantal-approach with a single completely optimized one-electron basis function nearly reproduces the Hartree-Fock-result. The two-determinantal approach with two optimized basis functions of type _g and _u nearly gives the optimized double configuration SCF result.     

Perturbation theory for the density matrix in the MO LCAO method

The problem of calculating the density matrix of perturbed electron shells in complex molecules has been solved in the MO LCAO [molecular orbital linear combination atomic orbitals] approximation. Considering the electron interaction in correcting the density matrix in each order of the perturbation theory yields one matrix equation uniquely defining the correction. This equation is formulated as an ordinary, linear nonhomogeneous system with respect to the individual matrix elements; this provided the possibility of establishing the convergence region of the self-consistent method. The solution of this system is accomplished by the method of steepest descent which, as has been shown, is always convergent in this case. It was established that just as in the conventional perturbation theory, the (2i + 1)th correction for the energy in terms of the Hartree-Fock method may be expressed by a correction of the density matrix not exceeding the i-th order. Detailed expressions are presented up to the eighth energy correction.     

The Bond Order of C2 from a Strictly N‐Representable Natural Orbital Energy Functional Perspective

The bond order of the ground electronic state of the carbon dimer has been analyzed in the light of natural orbital functional theory calculations carried out with an approximate, albeit strictly N‐representable, energy functional. Three distinct solutions have been found from the Euler equations of the minimization of the energy functional with respect to the natural orbitals and their occupation numbers, which expand upon increasing values of the internuclear coordinate. In the close vicinity of the minimum energy region, two of the solutions compete around a discontinuity point. The former, corresponding to the absolute minimum energy, features two valence natural orbitals of each of the following symmetries, σ, σ*, π and π*, and has three bonding interactions and one antibonding interaction, which is very suggestive of a bond order large than two but smaller than three. The latter, features one σ–σ* linked pair of natural orbitals and three degenerate pseudo‐bonding like orbitals, paired each with one triply degenerate pseudo‐antibonding orbital, which points to a bond order larger than three. When correlation effects, other than Hartree–Fock for example, between the paired natural orbitals are accounted for, this second solution vanishes yielding a smooth continuous dissociation curve. Comparison of the vibrational energies and electron ionization energies, calculated on this curve, with their corresponding experimental marks, lend further support to a bond order for C ₂ intermediate between acetylene and ethylene.     

An analysis of the complex molecular orbitals method

Numerical results for the ground state of the HN ₂ ⁺ and HCO⁺ molecular ions at their near equilibrium geometry, obtained by the complex molecular orbitals (CMO) method in the extended basis set, are reported. The CMO wavefunction of the HN ₂ ⁺ ion is compared with the CI wavefunction obtained in the same basis set. This reveals the nature of approximations inherent in the CMO method. A peculiar feature of the occupation numbers of the CMO natural orbitals is also explained.Alexander von Humboldt Fellow. On leave from the Institute Rudjer Bokovi, Zagreb, Croatia, Yugoslavia.     

The use of population localised orbitals to interpret molecular orbital calculations

An efficient and general method is derived to calculate population localised molecular orbitals (LMO's) from a given SCF eigenvector matrix, by reduction to an eigenvalue problem. Applications to both localised molecules (NH₃ and C₂H₂) and delocalised ones (B₂H₆, C₆H₆ and butadiene) are discussed in some detail. It is shown that unequal occupation of atomic energy levels leads to non-orthogonal LMO's. The consequences of non-orthogonal atomic hybrid orbitals are discussed, formulas for their overlap in terms of atomic occupation numbers are derived and it is shown that the occupation numbers are connected to LMO atomic orbital coefficients by various sum rules.     

The structure of the second-order reduced density matrix in density matrix functional theory and its construction from formal criteria

Some formal requirements for the second-order reduced density matrix are discussed in the context of density matrix functional theory. They serve as a basis for the ad hoc construction of the second-order reduced density matrix in terms of the first-order reduced density matrix and lead to implicit functionals where the occupation numbers of the natural orbitals are obtained as diagonal elements of an idempotent matrix the elements of which represent the variational parameters to be optimized. The numerical results obtained from a first realization of such an implicit density matrix functional give excellent agreement with the results of full configuration interaction calculations for four-electron systems like LiH and Be. Results for H2O taken as an example for a somewhat larger molecule are numerically less satisfactory but still give reasonable occupation numbers of the natural orbitals and indicate the capability of density matrix functional theory to cope with static electron correlation.     

Molecular states of atoms. I. Orbital energy parameters

Atomic valence state energies are analyzed to obtain values of orbital energy parameters that may be used in semiempirical molecular orbital calculations. Difficulty in defining the interaction between orbitals with non-integer electron populations is systematically avoided by distinguishing between a valence state and a molecular state of an atom, only the latter state having non-integer spin paired orbital occupancy. Application of the virial theorem to the molecular state enables a value for the orbital kinetic energy to be obtained from the valence state orbital energy parameters once an arbitrary configuration is defined as reference. The orbitals then are eigenfunctions of the atomic Fock operator for that reference molecular state and, with their energy parameters, may be employed as a fixed basis set for molecular orbital calculations.     

The perturbation theory of the extended Hartree–Fock approximation for two-electron atoms

The first-order 1/Z perturbation theory of the extended Hartree–Fock approximation for two-electron atoms is described. A number of unexpected features emerge: (a) it is proved that the orbitals must be expanded in powers of Z^?1/2, rather than in Z^?1 as expected; (b) it is shown that the restricted Hartree–Fock and correlation parts of the orbitals can be uncoupled to first order, so that second-order energies are additive; (c) the equation describing the first-order correlation orbital has an infinite number of solutions of all angular symmetries in general, rather than only one of a single symmetry as expected; (d) the first-order correlation equation is a homogeneous linear eigenvalue-type equation with a non-local potential. It involves a parameter μ and an eigenvalue ω(μ) which may be interpreted as the probability amplitude and energy of a virtual correlation state. The second-order correlation energy is 2μ²ω. Numerical solutions for the first-order correlation orbitals, obtained variationally, are presented. The approximate second-order correlation energy is nearly 90% of the exact value. The first-order 1/Z perturbation theory of the natural-orbital expansion is described, and the coupled first-order integro-differential perturbation equations are obtained. The close relationship between the first-order extended Hartree–Fock correlation orbitals and the first-order natural correlation orbitals is discussed. A comparison of the numerical results with those of Kutzelnigg confirms the similarity.     

Comparison of Magnetic Properties between Spin Singlet and Triplet Li<Subscript>3</Subscript>Al<Stack><Subscript>4</Subscript><Superscript>−</Superscript></Stack> Clusters

With a second-order Møller–Plesset perturbation theory and Hartree–Fock nuclear magnetic resonance calculations, we investigated the magnetic properties of spin singlet and triplet Li₃Al ₄ ^? clusters. The obtained gauge-independent atomic orbital magnetic shielding tensors confirm the paramagnetism of singlet Li₃Al ₄ ^? and diamagnetism of the triplet. The planar rings composed of four aluminum atoms make the magnetic properties of Li₃Al ₄ ^? clusters versatile. The localized molecular orbital, low symmetry of geometric conformation and narrow gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital are found to correlate with the paramagnetism of singlet Li₃Al ₄ ^? . The origin of the paramagnetism is explained. In triplet Li₃Al ₄ ^? , the two outmost orbitals are degenerate, causing a conversion from the paramagnetism to diamagnetism.     

Application of the many-body perturbation theory to normal saturated hydrocarbons in the localized representation

The second-order energy corrections are calculated for some normal saturated hydrocarbons by using the many body-perturbation theory (MBPT) based on localized orbitals. The correlation energies are expressed as the sum of contributions from virtual orbital pairs. We have found that these contributions are transferable and have interesting structural features: the trans-coplanar effects are relatively large. Partitioning the correlation energies according to the order of neighbourhood we have found that the zero order effects are the largest but the first and second neighbour contributions are also important.Dedicated to Professor J. Koutecký on the occasion of his 65th birthday     

A variational perturbation treatment of the hydrogen molecular ion using the green function method on the united atom

We calculated the second-order perturbed energy and the first-order perturbed wave function of the hydrogen molecular ion by the Green function method based on a united atom with effective nuclear charge, Z. Since the perturbation is given by (Z/r) — (1/r_A) — (1/r_B), its magnitude depends on the parameter Z. We obtain an analytical expression for the second-order perturbed energy as a function of Z and the internuclear distance R. We discuss a variety of methods to determine the best values for Z and to derive the potential energy curve and the force constant for H; these methods vary in the accuracy of their results.     

A local second-order Møller-Plesset method with localized orbitals: a parallelized efficient electron correlation method

Using orthogonal localized occupied orbitals we have developed and implemented a parallelized local second-order M?ller-Plesset (MP2) method based on the idea developed by Head-Gordon and co-workers. A subset of nonorthogonal correlation functions (the orbital domain) was assigned to each of the localized occupied orbitals using a distance criterion and excitations from localized occupied orbitals that were arranged into subsets. The correlation energy was estimated using a partial diagonalization and an iterative efficient method for solving large-scale linear equations. Some illustrative calculations are provided for molecules with up to 1484 Cartesian basis sets. The orbital domain sizes were found to be independent of the molecular size, and the present local MP2 method covered about 98%-99% of the correlation energy of the conventional canonical MP2 method.     

Calculation of Two-Center Nuclear Attraction Integrals over Integer and Noninteger n-Slater Type Orbitals in Nonlined-Up Coordinate Systems

Two-center nuclear attraction integrals over Slater type orbitals with integer and noninteger principal quantum numbers in nonlined up coordinate systems have been calculated by means of formulas in our previous work (T. Özdoan and M. Orbay, Int. J. Quant. Chem. 87 (2002) 15). The computer results for integer case are in best agreement with the prior literature. On the other hand, the results for noninteger case are not compared with the literature due to the scarcity of the literature, but also compared with the limit of integer case and good agreements are obtained. The proposed algorithm for the calculation of two-center nuclear attraction integrals over Slater type orbitals with noninteger principal quantum numbers in nonlined-up coordinate systems permits to avoid the interpolation procedure used to overcome the difficulty introduced by the presence of noninteger principal quantum numbers. Finally, numerical aspects of the presented formulae are analyzed under wide range of quantum numbers, orbital exponents and internuclear distances.     

SCF perturbation theory for unimolecular reactions

A CNDO/2 SCF perturbation theory is presented for interpreting the form of CNDO/2 potential energy surfaces of unimolecular reactions. The analysis is performed by calculating the energy change E arising from a distortion of the molecular geometry along the reaction coordinate. E is decomposed into different perturbational contributions which are appropriate for an interpretation of the perturbation energy E. Moreover, E is resolved into energy parts arising from a single occupied orbital and contributions due to pairwise orbital interactions. In this way one evaluates numerically how the form of the occupied and unoccupied orbitals determines the magnitude of E. If the distortion occurs along a definite symmetry coordinate, group-theoretical arguments can be applied to discuss the magnitude of characteristic components of the perturbation energy. The SCF perturbation theory is used to analyze the isomerization of ethylene, cis-2-butene and cis-2-butenenitrile.This work was partially supported by Nato-Grant No. 1072     

Molecules under external electric field: On the changes in the electronic structure and validity limits of the theoretical predictions

A detailed analysis of the electronic structure of the ground and first excited spin state of three diatomic molecules ( and ) under static applied electric field is performed at CCSD(T), DFT, MRCI and MRCI(Q) levels of theory. Our findings have revealed that by boosting the applied field one induces changes in the occupation numbers of molecular orbitals, giving rise to changes in the equilibrium geometry and in the HOMO–LUMO energy gap. Specifically, singlet to triplet spin transition can be induced by increasing the applied electric field beyond a critical value. Accordingly, affecting the accuracy of the widely used expression of energy expanded in Taylor series with respect to the applied electric field. © 2018 Wiley Periodicals, Inc.