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

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

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."     

The Fock Matrix Analysis for Atomic Orbitals in Molecular Orbitals I. A New Look on the Covalent Bond in H2?

The character of the molecular orbitals can be better accounted for in terms of molecular adapted atomic orbitals and the Fock matrix expanded in these atomic orbital sets. A clean‐cut and unique criterion for the diradicals and the covalent bonds can be given for the molecular orbitals in both restricted and unrestricted Hartree‐Fock wavefunctions. Instead of the picture that overlap charge migrates into the bonding region, the new analysis displays another picture that the charge densities for the electrons with α and β spins give rise to two opposite spin density shifts. If the α one shifts from atom A toward atom B then it is vice versa for the β one. The spin density shifts proceed until the bonding molecular orbitals form.     

Analytic projection from plane‐wave and PAW wavefunctions and application to chemical‐bonding analysis in solids

Quantum‐chemical computations of solids benefit enormously from numerically efficient plane‐wave (PW) basis sets, and together with the projector augmented‐wave (PAW) method, the latter have risen to one of the predominant standards in computational solid‐state sciences. Despite their advantages, plane waves lack local information, which makes the interpretation of local densities‐of‐states (DOS) difficult and precludes the direct use of atom‐resolved chemical bonding indicators such as the crystal orbital overlap population (COOP) and the crystal orbital Hamilton population (COHP) techniques. Recently, a number of methods have been proposed to overcome this fundamental issue, built around the concept of basis‐set projection onto a local auxiliary basis. In this work, we propose a novel computational technique toward this goal by transferring the PW/PAW wavefunctions to a properly chosen local basis using analytically derived expressions. In particular, we describe a general approach to project both PW and PAW eigenstates onto given custom orbitals, which we then exemplify at the hand of contracted multiple‐ζ Slater‐type orbitals. The validity of the method presented here is illustrated by applications to chemical textbook examples—diamond, gallium arsenide, the transition‐metal titanium—as well as nanoscale allotropes of carbon: a nanotube and the fullerene. Remarkably, the analytical approach not only recovers the total and projected electronic DOS with a high degree of confidence, but it also yields a realistic chemical‐bonding picture in the framework of the projected COHP method. © 2013 Wiley Periodicals, Inc.     

Intuitive local picture for spin polarization of chemical bonds

The spin polarization of chemical bonds near radical centers is investigated by an analytical spin unrestricted Hartree–Fock model with numerical examples. The centroid analysis of localized molecular orbitals is also introduced to obtain an intuitive local picture for the spin polarization. The alternation of spin alignments in molecules are discussed with orbital symmetries and Hund's rule through chemical bonds. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

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.     

磷的d轨道在H~3PO分子中的作用

用分子片轨道在分子环境中发生极化的概念研究d轨道在H~3PO分子中的作用。H~3PO分子被分为两个分子片---H~3P和O.在RHF/6-31G^*水平上计算出分子环境中的极化了的分子片轨道(FOM)。再剔除d函数为主的FOM,用剩余的FOM为基进行构型优化,得到与RHF/6-31G^*相近的结果。这一结果说明磷原子的d函数在H~3PO分子中仅仅起一个极化函数的作用,而不是起价轨道作用。     

Orbital deletion procedure and its applications

The orbital deletion procedure is introduced, which is suited to quantitatively investigating the electronic delocalization effiect in earboeations and boranes. While the routine, ab initio molecular orbital methods can generate wavefunetions for real systems where all electrons are delocalized, the present orbital deletion procedure can generate wavefunctions for hypothetical reference molecules where electronic delocalization effect is deactivated. The latter wavefunetion normlly corresponds In the most stable resonance structure in terms of the resonance theory. By comparing and analyzing the delocalized and the localized wavefunetions, one can obtain a quantitative and instinct pieture to show how electronic deloealizalion inside a molecule affects the molecular structure, energy as well as other physical properties. Two examples are detailedly discussed. The first is related to the hypercoujugation of alkyl groups in carbocations and a comparison of the order of stability of carbocations is made, T     

Role of s-p orbital mixing in the bonding and properties of second-period diatomic molecules

Qualitative molecular orbital theory is widely used as a conceptual tool to understand chemical bonding. Symmetry-allowed orbital mixing between atomic or fragment orbitals of different energies can greatly complicate such qualitative interpretations of chemical bonding. We use high-level Amsterdam Density Functional calculations to examine the issue of whether orbital mixing for some familiar second-row homonuclear and heteronuclear diatomic molecules results in net bonding or antibonding character for a given molecular orbital. Our results support the use of slopes of molecular orbital energy versus bond distance plots (designated radial orbital-energy slope: ROS) as the most useful criterion for making this determination. Calculated atomic charges and frontier orbital properties of these molecules allow their acid-base chemistry, including their reactivities as ligands in coordination chemistry, to be better understood within the context of the Klopman interpretation of hard and soft acid-base theory. Such an approach can be extended to any molecular species.     

Intramolecular interactions of L‐phenylalanine: Valence ionization spectra and orbital momentum distributions of its fragment molecules

Intramolecular interactions between fragments of L ‐phenylalanine, i.e., phenyl and alaninyl, have been investigated using dual space analysis (DSA) quantum mechanically. Valence space photoelectron spectra (PES), orbital energy topology and correlation diagram, as well as orbital momentum distributions (MDs) of L ‐phenylalanine, benzene and L ‐alanine are studied using density functional theory methods. While fully resolved experimental PES of L ‐phenylalanine is not yet available, our simulated PES reproduces major features of the experimental measurement. For benzene, the simulated orbital MDs for 1e_1g and 1a_2u orbitals also agree well with those measured using electron momentum spectra. Our theoretical models are then applied to reveal intramolecular interactions of the species on an orbital base, using DSA. Valence orbitals of L ‐phenylalanine can be essentially deduced into contributions from its fragments such as phenyl and alaninyl as well as their interactions. The fragment orbitals inherit properties of their parent species in energy and shape (ie., MDs). Phenylalanine orbitals show strong bonding in the energy range of 14‐20 eV, rather than outside of this region. This study presents a competent orbital based fragments‐in‐molecules picture in the valence space, which supports the fragment molecular orbital picture and building block principle in valence space. The optimized structures of the molecules are represented using the recently developed interactive 3D‐PDF technique. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

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.     

Binary (e,2e) spectroscopic study of carbonyl sulfide and the momentum-space chemistry of CO2, CS2 and OCS

The valence-shell binding energy spectra (8–44 eV) and molecular orbital momentum distributions of OCS have been studied by non-coplanar symmetric binary (e,2e) spectroscopy. Existing theoretical binding energy spectra calculated using the many-body 2ph-TDA Green's function (GF) method and using the symmetry-adapted cluster (SAC) on method are compared with the experiment. Intense many-body structure in the measured and calculated binding energy spectra indicates the general breakdown of the independent particle ionization picture. Experimental momentum distributions are compared with those calculated using ab initio SCF wavefunctions of minimal basis set quality and of near Hartree—Fock quality. Excellent agreement between the experimental momentum distributions and those calculated by the near Hartree—Fock wavefunction is obtained for the three innermost valence orbitals: 8σ, 7σ and 6σ. The correct order of the close lying outer-valence 2π and 9σ orbitals is unambiguously identified from the shapes of the measured momentum distributions. Momentum and position contour density maps computed from theoretical wavefunctions of near Hartree—Fock quality are used to interpret the shapes and atomic characters of the observed momentum distributions. The momentum densities of the outermost-valence antibonding π orbitals and of the outermost-valence bonding σ orbitals of the linear triatomic group: CO₂, CS₂ and OCS are compared respectively with each other. The associated chemical trends are discussed within the existing framework of momentum-space chemical principles.     

The Bonding Nature of the Canonical Molecular Orbitals in Simple Diatomic Molecules

The bonding nature of the canonical molecular orbitals 2σg, 2σu and 3σg in the molecules N₂,O₂, F₂ and the related analogous molecular orbitals in the molecules P₂ and CO, is analysed using Weinhold's natural bond orbital set. When the canonical molecular orbitals can be well localized into natural bond orbitals, the covalent bond can be completely attributed to the bonding type natural bond orbitals. The decomposition of canonical molecular orbitals into the natural bond orbital basis then gives the weighted bond order as the component of the bonding portion in the canonical molecular orbital. The weighted bond order results match the photoelectron spectroscopy assignment quite satisfactorily.     

用分子片轨道逐步剔除法研究电子离域

估计分子中不同类型的电子离域作用（如p-π → d-π和p-π → σ~*）的相 对强弱对理解分子中化学键的本质有关键的作用。剔除某一分子片轨道（d-π或σ ~*）后分子体系能量的改变可用来计算电子离域到该轨道的离域能。由于轨道之间 的相互作用,使离域能的计算与轨道剔除的次序有关。为克服这种不确定性,可以 逐步轮流增加某一对特定轨道（d-π和σ~*）的库仑积分,以使这对轨道在分子波 函数中的比重逐步减少,即将这对轨道轮流逐步剔除。这样可将轨道间的相互影响 减小以至消除,从而得到各轨道的精确的离域能。以H_3PO中P-O键为例说明了轨道 逐步剔除方法的应用。     

Maximum bond order hybrid orbitals

Summary Based on the simplified calculation scheme of the maximum bond order principle and the basic idea of the maximum overlap symmetry orbital method, a simple procedure is suggested for constructing systematically the bonding hybrid orbitals, called maximum bond order hybrid orbitals, for a given molecule from the first-order density matrix obtained from a molecular orbital calculation. As an example, the proposed procedure is performed for some typical small molecules by use of the density matrix obtained from CNDO/2 calculation. It is shown that the bonding hybrid orbitals constructed by using the procedure are extremely close to those by using the natural hybrid orbital procedure and in good agreement with chemical intuition, and that the proposed procedure can be performed more easily than the natural hybrid orbital procedure and can given simultaneously the values of the maximum bond order for all bonds in molecules.The project was supported by National Natural Science Foundation of China and also supported partly by Foundation of Hubei Education Commission     

On the suitability of strictly localized orbitals for hybrid QM/MM calculations

In the QM/MM method we have developed (LSCF/MM), the QM and the MM parts are held together by means of strictly localized bonding orbitals (SLBOs). Generally these SLBOs are derived from localized bond orbitals (LBOs) that undergo tails deletion, resulting in a nonpredictable change of their properties. An alternative set of SLBOs is provided by the extremely localized molecular orbitals (ELMOs) approach, where the orbitals are rigorously localized on some prefixed atoms without tails on the other atoms of the molecule. A comparative study of SLBOs arising from various localization schemes and ELMOs is presented to test the reliability and the transferability of these functions within the Local Self-Consistent Field (LSCF) framework. Two types of chemical bonds were considered: C--C and C--O single bonds. The localized functions are obtained on the ethane and the methanol molecules, and are tested on beta-alanine and diethyl ether molecules. Moreover, the various protonation forms of beta-alanine have been investigated to illustrate how well the polarity variation of the chemical bond can be handled throughout a chemical process. At last, rotation energy profiles around C--C and C--O bonds are reproduced for butane and fluoromethanol. Energetic, geometric, as well as electronic factors all indicate that ELMO functions are much more transferable from one molecule to another, leading to results closer to the usual SCF reference than any other calculations involving any other localized orbitals. When the shape of the orbital is the most important factor then ELMO functions will perform as well as any other localized orbital.