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.     

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.     

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.     

Localization of molecular orbitals on fragments

A non‐iterative algorithm for the localization of molecular orbitals (MOs) from complete active space self consistent field (CASSCF) and for single‐determinantal wave functions on predefined moieties is given. The localized fragment orbitals can be used to analyze chemical reactions between fragments and also the binding of fragments in the product molecule with a fragments‐in‐molecules approach by using a valence bond expansion of the CASSCF wave function. The algorithm is an example of the orthogonal Procrustes problem, which is a matrix optimization problem using the singular value decomposition. It is based on the similarity of the set of MOs for the moieties to the localized MOs of the molecule; the similarity is expressed by overlap matrices between the original fragment MOs and the localized MOs. For CASSCF wave functions, localization is done independently in the space of occupied orbitals and active orbitals, whereas, the space of virtual orbitals is mostly uninteresting. Localization of Hartree–Fock or Kohn–Sham density functional theory orbitals is not straightforward; rather, it needs careful consideration, because in this case some virtual orbitals are needed but the space of virtual orbitals depends on the basis sets used and causes considerable problems due to the diffuse character of most virtual orbitals. © 2012 Wiley Periodicals, Inc.     

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.     

虚轨道定域化

通过引入虚轨道定域化函数，扩展了集居数的定域化方法，在从头算STO－3G水平上用护展后的方法计算无机小分子、烷烃、卤代烷、醇、胺及共轭烃等，结果表明，定域占据轨道与定域虚轨道的系数和能级分布相似的规律。     

Application of many-body perturbation theory in the localized representation for the all-trans conjugated polyenes

Ab initio self-consistent field (SCF ) calculations are performed with the standard 6-31G* basis set for all-trans conjugated oligomers C_2n+2H_2n+4. The canonical occupied and virtual molecular orbitals (MO s) are separately localized by unitary transformations. Due to the localization, the perturbation operator is partitioned into two-particle and into single-particle terms; the MBPT is, therefore, a double-perturbation expansion in this case. By using the localized representation of the MBPT , the correlation energy contributions can be partitioned into local and nonlocal effects. It can be shown that the local effects are very important and well transferable, which makes possible the calculation of the correlation energy of larger molecules if the localized molecular orbitals (occupied and virtual) of smaller related molecules are known. © 1994 John Wiley & Sons, Inc.     

Direct selected multireference configuration interaction calculations for large systems using localized orbitals

A selected multireference configuration interaction (CI) method and the corresponding code are presented. It is based on a procedure of localization that permits to obtain well localized occupied and virtual orbitals. Due to the local character of the electron correlation, using local orbitals allows one to neglect long range interactions. In a first step, three topological matrices are constructed, which determine whether two orbitals must be considered as interacting or not. Two of them concern the truncation of the determinant basis, one for occupied/virtual, the second one for dispersive interactions. The third one concerns the truncation of the list of two electron integrals. This approach permits a fine analysis of each kind of approximation and induces a huge reduction of the CI size and of the computational time. The procedure is tested on linear polyene aldehyde chains, dissociation potential energy curve, and reaction energy of a pesticide-Ca(2+) complex and finally on transition energies of a large iron system presenting a light-induced excited spin-state trapping effect.     

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.     

Local Hartree–Fock orbitals using a three‐level optimization strategy for the energy

Using the three‐level energy optimization procedure combined with a refined version of the least‐change strategy for the orbitals—where an explicit localization is performed at the valence basis level—it is shown how to more efficiently determine a set of local Hartree–Fock orbitals. Further, a core–valence separation of the least‐change occupied orbital space is introduced. Numerical results comparing valence basis localized orbitals and canonical molecular orbitals as starting guesses for the full basis localization are presented. The results show that the localization of the occupied orbitals may be performed at a small computational cost if valence basis localized orbitals are used as a starting guess. For the unoccupied space, about half the number of iterations are required if valence localized orbitals are used as a starting guess compared to a canonical set of unoccupied Hartree–Fock orbitals. Different local minima may be obtained when different starting guesses are used. However, the different minima all correspond to orbitals with approximately the same locality. © 2013 Wiley Periodicals, Inc.     

Full spectral decomposition of ring currents

Within the ipsocentric method for calculation of molecular magnetic response, projection of perturbed orbitals onto the virtual orbital space allows partition of induced current density into contributions from individual virtual excitations between occupied and unoccupied orbitals, enabling detailed assignment of the origin of currents in, e.g., benzene, cyclooctatetraene, borazine, coronene, and corannulene. Whereas delocalized currents in benzene and planar cyclooctatetraene are described by transitions within the valence space, localized currents in the borazine pi system involve excitations outside the valence space.     

Localization scheme for relativistic spinors

A new method to determine localized complex-valued one-electron functions in the occupied space is presented. The approach allows the calculation of localized orbitals regardless of their structure and of the entries in the spinor coefficient matrix, i.e., one-, two-, and four-component Kramers-restricted or unrestricted one-electron functions with real or complex expansion coefficients. The method is applicable to localization schemes that maximize (or minimize) a functional of the occupied spinors and that use a localization operator for which a matrix representation is available. The approach relies on the approximate joint diagonalization (AJD) of several Hermitian (symmetric) matrices which is utilized in electronic signal processing. The use of AJD in this approach has the advantage that it allows a reformulation of the localization criterion on an iterative 2 × 2 pair rotating basis in an analytical closed form which has not yet been described in the literature for multi-component (complex-valued) spinors. For the one-component case, the approach delivers the same Foster-Boys or Pipek-Mezey localized orbitals that one obtains from standard quantum chemical software, whereas in the multi-component case complex-valued spinors satisfying the selected localization criterion are obtained. These localized spinors allow the formulation of local correlation methods in a multi-component relativistic framework, which was not yet available. As an example, several heavy and super-heavy element systems are calculated using a Kramers-restricted self-consistent field and relativistic two-component pseudopotentials in order to investigate the effect of spin-orbit coupling on localization.     

An extension of the pairing theorem

An extended pairing scheme is presented which ensures the fulfillment of pairing conditions not only between the sets of occupied orbitals for spin α and for spin β, but also between their orthogonal complements, i.e., the sets of virtual orbitals for spin α and spin β, as well as between occupied orbitals for spin α and virtual orbitals for spin β and between virtual orbitals for spin α and occupied orbitals for spin β. It is shown that the extended pairing properties are suggested by some aspects of the construction of alternant molecular orbitals. The algorithm for singular value decomposition of rectangular matrices is proposed for use in practical implementations of the (extended) pairing scheme.     

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.     

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.     

A localized basis that allows fast and accurate second-order Moller-Plesset calculations

We present a method for computing a basis of localized orthonormal orbitals (both occupied and virtual), in whose representation the Fock matrix is extremely diagonal dominant. The existence of these orbitals is shown empirically to be sufficient for achieving highly accurate second-order Moller-Plesset (MP2) energies, calculated according to Kapuy's method. This method (which we abbreviate KMP2) involves a different partitioning of the n-electron Hamiltonian and scales at most quadratically, with potential for linearity, in the number of electrons. As such, we believe the KMP2 algorithm presented here could be the basis of a viable approach to local-correlation calculations.     

Fast noniterative orbital localization for large molecules

We use Cholesky decomposition of the density matrix in atomic orbital basis to define a new set of occupied molecular orbital coefficients. Analysis of the resulting orbitals ("Cholesky molecular orbitals") demonstrates their localized character inherited from the sparsity of the density matrix. Comparison with the results of traditional iterative localization schemes shows minor differences with respect to a number of suitable measures of locality, particularly the scaling with system size of orbital pair domains used in local correlation methods. The Cholesky procedure for generating orthonormal localized orbitals is noniterative and may be made linear scaling. Although our present implementation scales cubically, the algorithm is significantly faster than any of the conventional localization schemes. In addition, since this approach does not require starting orbitals, it will be useful in local correlation treatments on top of diagonalization-free Hartree-Fock optimization algorithms.     

Large-scale RPA calculations of chiroptical properties of organic molecules: Program RPAC

The computational considerations involved in calculating ordinary and rotatory intensities and electronic excitation energies in the random phase approximation (RPA ) are examined. We employ a localized orbital formulation in order to analyze the results in terms of local and charge-transfer excitations. Occupied orbitals are localized by the Foster–Boys procedure. The virtual space is transformed into a localized “valence” set that maximizes dipole strengths with the occupied counterparts, and a delocalized remainder. The two-electron integral transformation is performed with an efficient algorithm, based on Diercksen's, that generates only the particle–hole-type integrals required in the RPA . The lowest solutions of the RPA equations are obtained iteratively using a modification of the Davidson-Liu simultaneous vector expansion method. This allows the inclusion of the entire set of particle–hole states supported by a basis set of up to 102 orbitals. Calculations at this level give better excitation energies and intensities than SDCI methods, at substantial savings in computational effort. Comparative timings, computed results and analysis in terms of localized orbitals are given for planar and distorted ethylene using extended atomic orbital bases including diffuse functions. The results for planar ethylene are in excellent agreement with experiment.     

Localization maps by orbital partitioning of the electron density

Summary We define a localization measure for one-determinantal wave-functions based on the partitioning of the total electron density to orbital contributions. The set of occupied orbitals is the more localized the fewer terms are necessary to describe the total density. This measure varies from point to point in space which leads to characteristic localization maps for molecules.Supported in part by the grant Hungarian Research Fund OTKA 517/1990