Transferability of some properties of localized molecular orbitals

The localized molecular orbitals of some related ten- and eighteen-electron systems have been studied. The transferability of the kinetic, self-interaction, Coulomb and exchange interaction energies on localized orbitals have been shown. The standard deviation of the kinetic and of the interaction energies (including exchange) are less than 2.5% except for lone pair orbitals of the oxygen atoms where the standard deviation is close to 4%.     

On localized orbitals in infinite,periodic systems

It is shown that in two- or three-dimensional periodic systems with an odd number of electrons per unit cell no unitary transformation of the canonical Hartree - Fock orbitals can produce equivalent, localized orbitals for all electrons. However, in one-dimensional systems such orbitals can be found. The physical implications are discussed.     

An efficient linear scaling procedure for constructing localized orbitals of large molecules based on the one-particle density matrix

We have developed a linear-scaling algorithm for obtaining the Boys localized molecular orbitals from the one-particle density matrix. The algorithm is made up of two steps: the Cholesky decomposition of the density matrix to obtain Cholesky molecular orbitals and the subsequent Boys localization process. Linear-scaling algorithms have been proposed to achieve linear-scaling calculations of these two steps, based on the sparse matrix technique and the locality of the Cholesky molecular orbitals. The present algorithm has been applied to compute the Boys localized orbitals in a number of systems including α-helix peptides, water clusters, and protein molecules. Illustrative calculations demonstrate that the computational time of obtaining Boys localized orbitals with the present algorithm is asymptotically linear with increasing the system size.     

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.     

一种计算量小的构造合理最紧缩定域轨道的方法

提出一种计算量小的构造合理最紧缩定域化轨道的方法.非正交定域轨道(NOLMO)没有“正交化尾巴”,比正交定域轨道(OLMO)更局域化、更紧缩、更具可移植性,从而更适合用于化学问题的理论研究.但若在用变分法确定最大限度局域化NOLMO时单纯取消正交条件而不附加其它的限制,会得到趋于线性相关的不合理结果.提出用强制NOLMO与投影自然键轨道(PNBO)的重心重合代替正交化条件确定合理的最大限度局域化NOLMO(合理最紧缩定域轨道)的方法.对一系列不同类型分子的计算结果表明,用该方法可以得到空间分布合理及线性独立的NOLMO,延伸度与文献已有的最佳结果接近,而计算量大幅度降低.由于构造投影自然键轨道的计算量随分子中原子数目的增加只是线性增长,因此,该方法可以用来构造较大体系的合理最紧缩定域轨道     

Superexchange in magnetic insulators: an interpretation of the metal-metal charge transfer energy in the Anderson theory

The superexchange interactions in four three-center model systems A-L-B, for A and B being paramagnetic centers and L a diamagnetic bridging ligand, are analyzed by valence bond configuration interaction models in combination with fourth-order perturbation theory. We analyze the four distinct cases where a bridging ligand orbital simultaneously interacts with half-filled orbitals localized on A and B (case i), a half-filled orbital localized on A and an empty orbital localized on B (case ii), a full orbital localized on A and a half-filled orbital localized on B (case iii), and finally a full orbital localized on A and an empty orbital localized on B (case iv). In all four cases we compare our new results using localized orbitals with the equivalent results obtained using the Anderson ansatz of delocalized (magnetic) orbitals. The effective metal-to-metal electron transfer energy Ueff in the old formalism with delocalized orbitals is expressed in terms of the metal-to-metal electron transfer energy U and the ligand-to-metal electron transfer energy delta using localized orbitals. We find that the old formalism containing only Ueff is in general not sufficient. For cases i and ii we show that Ueff can be regarded as an effective U strongly reduced with respect to the free ion as a result of hybridization effects, whereas the same reduction of U for the cases iii and iv is not possible. The relevance and applicability of our theoretical results is demonstrated on magnetochemical data from the literature.     

LocalSCF method for semiempirical quantum-chemical calculation of ultralarge biomolecules

A linear-scaling semiempirical method, LocalSCF, has been proposed for the quantum-chemical calculations of ultralarge molecular systems by treating the large-scale molecular task as a variational problem. The method resolves the self-consistent field task through the finite atomic expansion of weakly nonorthogonal localized molecular orbitals. The inverse overlap matrix arising from the nonorthogonality of the localized orbitals is approximated by preserving the first-order perturbation term and applying the second-order correction by means of a penalty function. This allows for the separation of the orbital expansion procedure from the self-consistent field optimization of linear coefficients, thereby maintaining the localized molecular orbital size unchanged during the refinement of linear coefficients. Orbital normalization is preserved analytically by the variation of virtual degrees of freedom, which are orthogonal to the initial orbitals. Optimization of linear coefficients of localized orbitals is performed by a gradient procedure. The computer program running on a commodity personal computer was applied to the GroEL-GroES chaperonin complex containing 119,273 atoms.     

Localization of the filled and virtual orbitals in the nucleotide bases

For the four nucleotide bases cytosine, uracil, adenine and guanine both Boys (B) and Edminston-Ruedenberg (ER) localization procedures of the ab initio canonical orbitals have been performed. The results obtained for both σ-π separation and by treating all electrons together show a very good localization for all electrons (one-center lone-pairs and two-center localized orbitals even for π-electrons) and a rather good localization for the virtuals applying both B and ER criteria. The results of the two methods are essentially identical. These results suggest that the application of localized orbitals will open new possibilities for the calculation of correlation in extended systems.     

Localized orbitals,electric field gradients and diamagnetic shielding in hydrogen bonded systems

Ab-initio localized molecular orbitals are used to evaluate the diamagnetic shieldings and the electric field gradients in hydrogen bonded systems.     

Distortion of interacting atoms and ions

The distortion of atoms in HeNe and ions in LiF and NaCl is studied by use of non-orthogonal, least distorted, localized molecular orbitals in the restricted Hartree-Fock approximation. We find that the least distortion criterion is effective in producing localized orbitals. In addition we find that point-by-point the localized orbitals differ little from the Hartree-Fock orbitals of the corresponding atoms or ions, although the differences are energetically important. We infer from our calculations that the effective potential which distorts atomic orbitals into localized molecular orbitals must be quite weak.     

An efficient method for constructing nonorthogonal localized molecular orbitals

A new method for constructing nonorthogonal localized molecular orbitals (NOLMOs) is presented. The set of highly localized NOLMOs is obtained by minimization of the spread functional starting from an initial set of canonical orthogonal molecular orbitals. To enhance the stability and efficiency, the centroids of the NOLMOs are constrained to be those of the corresponding orthogonal localized molecular orbitals (OLMOs), which are obtained with the Boys criterion in advance. In particular, these centroid constraints make the optimization for each NOLMO independent of the others, which is an attractive feature for application to large systems. The minimization with the constraints incorporated through the multiplier-penalty function method is stable and efficient in convergence. While exhibiting the classical bonding pattern in chemistry and sharing a spatial distribution similar to that of the corresponding OLMOs, the obtained NOLMOs are more compact than the corresponding OLMOs with about 10%-28% reduction in the value of the spread functional and devoid of the troublesome "orthogonalization tails."     

Application of the many-body perturbation theory by using localized orbitals

Diagrammatic formulation of the many-body perturbation theory is investigated when both the occupied orbitals and the virtual ones are localized, i.e., they are unitary transforms of the canonical Hartree–Fock orbitals. All diagrams representing ground state correlation energy can be generated through fifth order. For cyclic polyenes C₆H₆ and C₁₀H₁₀ as model systems, the energy corrections are calculated in the Pariser–Parr–Pople approximation for a wide range of the coupling constant β^?1, through fourth order including some fifth order terms. The results are compared to those obtained by other methods: perturbation theory by using canonical orbitals and full CI. The effect of neglecting contributions from orbitals localized into neighboring sites is also studied.     

A study of the molecular orbital localization into an extended Hartree—Fock approach: Application to the BeH2 ground state

In this article, we present a study of the localization and properties of the molecular orbitals (MOs) of polyatomic systems by using a comprehensive version of the G1 model. In this version, the wave function is written as a DODS product of univocally determined spin orbitals (MOs), “projected” on the singlet ground state. A procedure for determining the MOs is given and applied to the BeH₂ ground state. Equivalent split shell and localized MOs are found. The Be orbitals are seen to exhibit sp hybridization and the localized valence MOs are found to produce − 13.7 kcal/mol localization energy. Multistructural calculations are carried out and show that the present approach is able to describe localized and well-oriented bonds whenever the molecule under study presents only a single well-defined nonresonant chemical structure. © 1996 John Wiley & Sons, Inc.     

An ab initio method for approximation of the frozen molecular fragment

An ab initio method for calculation on many-electron molecular systems with the approximation of the inactive part of a molecule by frozen molecular fragment is presented. In the following method the SCF calculations are performed in two series. First the molecular orbitals resulting from the first SCF calculation (modest basis set) are localized. In the second SCF run, the basis set is extended for the active part of the molecule, while molecular orbitals of the inactive part, selected from the localized set, are kept frozen. The results are in good agreement with the extended basis set calculation.     

Determination of extremely localized molecular orbitals in the framework of density functional theory

Extremely localized molecular orbitals are rigorously localized on only a preselected set of atoms and do not have any tails outside the localization region. The importance of these orbitals lies in their ability to be transferred from one molecule to another one. A new algorithm to determine extremely localized molecular orbitals in the framework of the density functional theory method is presented. This could also be a valuable tool in the quantum mechanics/molecular mechanics methodology where localized molecular orbitals are used to describe covalent bonds across the frontier region. The present approach is used to build up the electron density of thymopentin, a polypeptide constituted by five residues, starting from extremely localized molecular orbitals determined on a set of model molecules. The results obtained confirm good transferability properties for these orbitals.Proceedings of the 11th International Congress of Quantum Chemistry satellite meeting in honor of Jean-Louis Rivail     

Extremely localized molecular orbitals (ELMO): a non-orthogonal Hartree-Fock method

A new optimization method for extremely localized molecular orbitals (ELMO) is derived in a non-orthogonal formalism. The method is based on a quasi Newton-Raphson algorithm in which an approximate diagonal-blocked Hessian matrix is calculated through the Fock matrix. The Hessian matrix inverse is updated at each iteration by a variable metric updating procedure to account for the intrinsically small coupling between the orbitals. The updated orbitals are obtained with approximately n ² operations. No n ³ processes such as matrix diagonalization, matrix multiplication or orbital orthogonalization are employed. The use of localized orbitals allows for the creation of high-quality initial “guess” orbitals from optimized molecular orbitals of small systems and thus reduces the number of iterations to converge. The delocalization effects are included by a Jacobi correction (JC) which allows the accurate calculation of the total energy with a limited number of operations. This extension, referred to as ELMO(JC), is a variational method that reproduces the Hartree-Fock (HF) energy with an error of less than 2 kcal/mol for a reduced total cost compared to standard HF methods. The small number of variables, even for a very large system, and the limited number of operations potentially makes ELMO a method of choice to study large systems. Received: 30 December 1996 / Accepted: 5 June 1997     

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.     

Potential magnetic properties of nanotubes (n, 0) with Klein and Fujita edges

Analytical solutions for localized states of zigzag-type nanotube (NT) fragments with various combinations of Klein and Fujita borders are considered using the Hückel approach. It is shown that the equations for determining molecular orbitals (MOs) in systems with two Klein edges are similar to equations for systems with two Fujita edges. An analytical formula for the energies of all ?? MOs is obtained for systems that have a Klein edge on one side and a Fujita edge on the other. It is established that these systems have n orbitals with energy ?? that are localized on the Fujita and Klein edges in dependence on the MO symmetry. The degeneracy of edge orbitals indicates that there is a tendency toward single occupancy of them and to the appearance of spin (magnetic) properties. In addition, the energies of the states of different multiplicity for NT fragments (8, 0) are calculated using the CASSCF approach. It is shown that the ground state has a multiplicity of 9, as was also indicated by estimates obtained using the density functional method (B3LYP). It is concluded that zigzag-type NTs with asymmetric edges have a tendency to exhibit spin properties. It is noted that the construction of nanoscale magnetic materials based on them is very promising.