Local orbitals by minimizing powers of the orbital variance

It is demonstrated that a set of local orthonormal Hartree-Fock (HF) molecular orbitals can be obtained for both the occupied and virtual orbital spaces by minimizing powers of the orbital variance using the trust-region algorithm. For a power exponent equal to one, the Boys localization function is obtained. For increasing power exponents, the penalty for delocalized orbitals is increased and smaller maximum orbital spreads are encountered. Calculations on superbenzene, C(60), and a fragment of the titin protein show that for a power exponent equal to one, delocalized outlier orbitals may be encountered. These disappear when the exponent is larger than one. For a small penalty, the occupied orbitals are more local than the virtual ones. When the penalty is increased, the locality of the occupied and virtual orbitals becomes similar. In fact, when increasing the cardinal number for Dunning's correlation consistent basis sets, it is seen that for larger penalties, the virtual orbitals become more local than the occupied ones. We also show that the local virtual HF orbitals are significantly more local than the redundant projected atomic orbitals, which often have been used to span the virtual orbital space in local correlated wave function calculations. Our local molecular orbitals thus appear to be a good candidate for local correlation methods.     

A simple extension of the external magnasco-perico localization procedure to the virtual MO-Space

The external localization procedure of Magnasco and Perico is extended to the unoccupied molecular orbitals of the Fock-operator. The formal correspondence between bonding orbitals and localized antibonding MOs is demonstrated. Localized occupied and virtual one-electron functions are calculated within a semiempirical INDO-Hamiltonian and are analyzed; the externally localized occupied MOs are compared with energy localized orbitals computed by the Edmiston and Ruedenberg procedure. Various applications of the fully localized (occupied and virtual) MO set are discussed.     

A regionally contracted multireference configuration interaction method: general theory and results of an incremental version

This work proposes to take benefit of the localizability of both occupied and virtual inactive molecular orbitals (MOs) in the context of complete active space singles and doubles configuration interaction (CAS-SDCI). The doubly occupied MOs are partitioned into blocks, or regions, corresponding to a subset of adjacent bonds and lone pairs. The localized virtual MOs are attributed to these regions from a spatial criterion. Then a series of limited post-CAS-CI calculations is performed, using the same reference space, one for each block, and then one per pair of blocks. From these independent CI calculations contracted external functions are defined for each block or for each pair of blocks, and for each state. A general multistate formalism is proposed, the CI matrix being expressed in the space defined by the CAS and the contracted functions. Preliminary numerical studies, resting on the evaluation of single-block and two-block contributions to the dynamical correlation energy of each state, are presented. Provided that size-consistency corrections are taken into account the results of the procedure are shown to be in excellent agreement with those of the nonpartitioned post-CAS-CI. The computational benefits of this evidently parallelizable procedure are underlined.     

A perspective on the localizability of Hartree–Fock orbitals

A common perception about molecular systems with a nonlocal electronic structure (as manifested by a nonlocal Hartree–Fock (HF) density matrix), such as conjugated π-systems, is that they can only be described in terms of nonlocal molecular orbitals. This view is mostly founded on chemical intuition, and further, this view is strengthened by traditional approaches for obtaining local occupied and virtual orbital spaces, such as the occupied Pipek–Mezey orbitals, and projected atomic orbitals. In this article, we discuss the limitations for localizability of HF orbitals in terms of restrictions posed by the delocalized character of the underlying density matrix for the molecular system and by the orthogonality constraint on the molecular orbitals. We show that the locality of the orbitals, in terms of nonvanishing charge distributions of orbitals centered far apart, is much more strongly affected by the orthogonality constraint than by the physical requirement that the occupied orbitals must represent the electron density. Thus, the freedom of carrying out unitary transformations among the orbitals provides the flexibility to obtain highly local occupied and virtual molecular orbitals, even for molecular systems with a nonlocal density matrix, provided that a proper localization function is used. As an additional consideration, we clear up the common misconception that projected atomic orbitals in general are more local than localized orthogonal virtual orbitals.     

Analyzing the electric response of molecular conductors using “electron deformation” orbitals and occupied‐virtual electron transfer

The concept of “electron deformation orbitals” (EDOs) is used to investigate the electric response of conducting metals and oligophenyl chains. These orbitals and their eigenvalues are obtained by diagonalization of the deformation density matrix (difference between the density matrices of the perturbed and unperturbed systems) and can be constructed as linear combinations of the unperturbed molecular orbitals within “frozen geometry” conditions. This form of the EDOs allows calculating the part of the electron deformation density associated to an effective electron transfer from occupied to virtual orbitals (valence to conduction band electron transfer in the band model of conductivity). It is found that the “electron deformation” orbitals pair off, displaying the same eigenvalue but opposite sign. Each pair represents an amount of accumulation/depletion of electron charge at different molecular regions. In the oligophenyl systems investigated only one pair contributes effectively to the charge flow between molecular ends, resulting from the promotion of electrons from occupied orbitals to close in energy virtual orbitals of appropriate symmetry and overlapping. Analysis of this pair along explains the differences in conductance of olygophenyl chains based on phenyl units. © 2014 Wiley Periodicals, Inc.     

The divide-expand-consolidate family of coupled cluster methods: numerical illustrations using second order Møller-Plesset perturbation theory

Previously, we have introduced the linear scaling coupled cluster (CC) divide-expand-consolidate (DEC) method, using an occupied space partitioning of the standard correlation energy. In this article, we show that the correlation energy may alternatively be expressed using a virtual space partitioning, and that the Lagrangian correlation energy may be partitioned using elements from both the occupied and virtual partitioning schemes. The partitionings of the correlation energy leads to atomic site and pair interaction energies which are term-wise invariant with respect to an orthogonal transformation among the occupied or the virtual orbitals. Evaluating the atomic site and pair interaction energies using local orbitals leads to a linear scaling algorithm and a distinction between Coulomb hole and dispersion energy contributions to the correlation energy. Further, a detailed error analysis is performed illustrating the error control imposed on all components of the energy by the chosen energy threshold. This error control is ultimately used to show how to reduce the computational cost for evaluating dispersion energy contributions in DEC.     

Localizability of dynamic electron correlation

The convergence of the intrapair correlation energy for a localized internal orbital is investigated as the virtual subspace is enlarged. At variance with previous investigations of this kind, the virtual subspace is represented in atomic orbitals. This allows to define spatial relations between the orbitals involved. Typically, over 98% of the pair correlation energy is recovered by a small local basis set, consisting of the valence orbitals of the atoms with which the electron pair is associated. This opens the possibility of an efficient Cl procedure based on localized pairs.     

Usefulness of modified virtual orbitals in multireference CI procedure illustrated by calculations on lithium clusters

The maximization of the exchange interaction between the canonical Hartree–Fock virtual and occupied orbitals leads to a transformed set of virtual orbitals which are well suited as one-electron functions for CI calculations. The procedure, generally known for a long time is seldom applied, despite its simplicity and very low computational demand. However, it is found to be particularly useful in the case of multireference CI, since an improved energy is obtained with a considerable shortening of the CI expansion. Moreover, in the final CI wave function, several configurations appear with considerable weight, thus allowing an easy choice of additional configurations to be inserted in the definition of a new zero-order wave function. The efficiency of the computational procedure is discussed for the case of a Li₆ cluster of D_3h symmetry and for the NaCO and PdCO complexes. Results are reported for the relative stability of four different geometrical arrangements of the Li₆ cluster.     

Fast localized orthonormal virtual orbitals which depend smoothly on nuclear coordinates

We present here an algorithm for computing stable, well-defined localized orthonormal virtual orbitals which depend smoothly on nuclear coordinates. The algorithm is very fast, limited only by diagonalization of two matrices with dimension the size of the number of virtual orbitals. Furthermore, we require no more than quadratic (in the number of electrons) storage. The basic premise behind our algorithm is that one can decompose any given atomic-orbital (AO) vector space as a minimal basis space (which includes the occupied and valence virtual spaces) and a hard-virtual (HV) space (which includes everything else). The valence virtual space localizes easily with standard methods, while the hard-virtual space is constructed to be atom centered and automatically local. The orbitals presented here may be computed almost as quickly as projecting the AO basis onto the virtual space and are almost as local (according to orbital variance), while our orbitals are orthonormal (rather than redundant and nonorthogonal). We expect this algorithm to find use in local-correlation methods.     

Pipek–Mezey localization of occupied and virtual orbitals

Recent advances in orbital localization algorithms are used to minimize the Pipek–Mezey localization function for both occupied and virtual Hartree–Fock orbitals. Virtual Pipek–Mezey orbitals for large molecular systems have previously not been considered in the literature. For this work, the Pipek–Mezey (PM) localization function is implemented for both the Mulliken and a Löwdin population analysis. The results show that the standard PM localization function (using either Mulliken or Löwdin population analyses) may yield local occupied orbitals, although for some systems the occupied orbitals are only semilocal as compared to state‐of‐the‐art localized occupied orbitals. For the virtual orbitals, a Löwdin population analysis shows improvement in locality compared to a Mulliken population analysis, but for both Mulliken and Löwdin population analyses, the virtual orbitals are seen to be considerably less local compared to state‐of‐the‐art localized orbitals. © 2013 Wiley Periodicals, Inc.     

Virtual orbital transformation prior to configuration interaction with localized orbitals

A virtual orbital transformation is proposed involving pairing of localized occupied orbitals with virtual orbitals. The virtual orbitals are transformed so that each virtual orbital is as “close” as possible to its occupied counterpart, where closeness is the inverse of the particular definition of localization. The appropriate transformation is derived for the special case of Foster–Boys localization, and an illustrative CNDO /2 calculation on HNO is presented. INDO CI results on the series N₂, CO, BF indicate that use of this transformation may reduce the number of energetically significant configurations.     

Valence-only model potential calculations on copper hydride molecule

Different sets of one-electron functions obtained according to the strong-orthogonal geminal theory (GEM) [1], the Generalized Molecular Orbital (GMO) method [2] and the exchange maximization between virtual and occupied orbitals (EVO) [3], are tested as basis for CI calculations. The efficiency of the three procedures is discussed investigating the electronic structure of the CuH molecule using an effective-core potential. The values computed for the bond length, the dissociation energy and the vibrational frequency of the ground electronic state are compared with the experimental ones. The charge distribution is examined to estimate the contribution of the d electrons to the Cu-H bond. Comparisons are made with the results obtained by other theoretical works in which the copper atom is treated as a one valence electron atom.     

Ab initio quantum chemistry calculations on the electronic structure of heavier alkyne congeners: diradical character and reactivity

The electronic structure of the heavier congeners of alkynes has been studied with emphasis on characterizing their extent of diradical character. Four orbitals play a crucial role in determining the electronic structure in planar trans-bent geometries. Two are associated with an out-of-plane pi interaction, pi and pi, and two are associated with in-plane interactions and/or in-plane lone pairs, LP(n-) and LP(n+). The ordering of these orbitals can change depending upon geometry. One extreme, corresponding to the local minimum for Si-Si and Ge-Ge, is a diradicaloid multiple-bonding configuration where LP and pi are nominally occupied. Another extreme, corresponding to a local minimum for Sn-Sn, is a relatively closed-shell single-bond configuration where LP and LP are nominally occupied. This ordering leads to predicted bond shortening upon excitation from singlet to triplet state. For the heavier elements, there appears to be very little energy penalty for large geometric distortions that convert from one ordering to the other on the singlet surface. The implications of these results with respect to experimental observations are discussed.     

Optimal virtual orbitals to relax wave functions built up with transferred extremely localized molecular orbitals

Extremely localized molecular orbitals (ELMOs), namely orbitals strictly localized on molecular fragments, are easily transferable from one molecule to another one. Hence, they provide a natural way to set up the electronic structure of large molecules using a data base of orbitals obtained from model molecules. However, this procedure obviously increases the energy with respect to a traditional MO calculation. To gain accuracy, it is important to introduce a partial electron delocalization. This can be carried out by defining proper optimal virtual orbitals that supply an efficient set for nonorthogonal configurations to be employed in VB-like expansions.     

One electron orbitals intrinsic to the reduced hamiltonian

One electron orbitals are determined from the reduced hamiltonian by a simple one-step diagonalization. These reduced hamiltonian orbitals (RHO's) are uniquely determined and virtual orbitals obtained in this procedure are on a par with filled orbitals. These RHO's appear well suited for CI calculations. Minimum basis set calculations are presented for H₂O and compared with similar SCF studies.     

Computational developments in generalized valence bond calculations

In this article a procedure for generating starting orbitals for generalized valence bond (GVB) calculations is presented. This is achieved by selecting orbitals which correspond to specific bonds or electron pairs. These orbitals can be identified from the localized molecular orbitals, for both occupied and virtual orbitals, which are obtained through a unitary transformation of the Hartree-Fock canonical molecular orbitals using the Boys's localization method. A scheme has also been implemented which achieves optimum convergence of the pairwise orbital optimization. An object-oriented GVB program is developed which automatically generates reliable initial GVB orbitals, leading to proper and fast convergence. © 1996 by John Wiley & Sons, Inc.     

Generalized-molecular-orbital theory: Simple multiconfiguration self-consistent-field method

Even after completing a multiconfiguration self-consistent-field (MCSCF ) calculation, one must often include additional configuration interaction (CI ) to obtain quantitative or semiquantitative results. There is some question of whether the prior MCSCF calculation is worthwhile, if additional CI is needed later. We have developed a new MCSCF computational method, which, because of our assumptions about the nature of the configurations, yields one Fock-like operator for all the “filled” orbitals (high occupation numbers) and a second Fock-like operator for all the “virtual” orbitals (low occupation numbers). Since there are only two matrices to build, our method is considerably faster than other MCSCF approaches. Because of these similarities to standard molecular-orbital (MO ) calculations, we have termed our approach generalized-molecular-orbital (GMO ) theory. However, the “virtual” orbitals, unlike those of standard MO theory, are optimized to correlate the “filled” ones and can he used in a subsequent CI calculation. Results are presented for the correlation energy of H₂O, the spectroscopic constants of N₂, the singlet–triplet energy separations in CH₂, and the nature of the chromium–chromium quadruple bond. Although these results are at a very low level of CI , the GMO approach appears to correct for the gross deficiencies of the single-determinant SCF procedure.     

Local orbitals for the study of the π→π* excitation in polyenes

Localized orbitals have recently been employed in large ab initio calculations, but their use has generally been restricted to ground‐state problems. In this work, we analyze the molecular orbitals of the excited states, optimized with a recently proposed local procedure. This method produces local orbitals of the CAS–SCF type, which permits its application to the study of excited states. In particular, we focus on the π→π* triplet excited state in polyenes, calculated using a 2/2 CAS space which includes two electrons in one π and one π* orbitals. In small polyenes, these two singly occupied active orbitals are delocalized all along the molecule. The extent of the delocalization is analyzed by studying polyenes of increasing size. Different polyenes have been studied, going from C₁₄H₁₆ to the C₇₀H₇₂ polyene. The relation of the π→π* excitation with the cation and anion systems is also discussed. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

DFT and CI Studies of Electronic Structure and Photoionization of Sc,Ti, V,Cr, and Co Tris-β-diketonate Complexes

Orbital structure calculations were performed in the density functional theory (DFT) approximation for neutral complexes of Sc, Ti, V, Cr, and Co tris-β-diketonates; for the first three compounds, the structures of the ground ionic states and ionization energies were calculated in the CI approximation with decomposition on the orbitals of DFT. The sequence of the highest occupied orbitals found by this procedure coincides with the order of bands in the PES spectrum, while in the SCF-HF ab initio method, it does not. After the electron removal, all orbitals are stabilized by about 4.5 eV; for the vanadium complex, the removal of one d electron leads to the greatest stabilization of the remaining occupied orbital, which is essentially a d orbital. In CI calculations, using the DFT orbitals for decomposition does not lead to significantly better agreement with experiment when compared to the single-determinantal approximation and to the CI method with orbitals of the ab initio approximation.Original Russian Text Copyright © 2004 by I. S. Osmushko and V. I. Vovna__________Translated from Zhurnal Strukturnoi Khimii, Vol. 45, No. 5, pp. 783–791, September–October, 2004.