Roussin红盐和Roussin黑盐的电子结构

Roussin黑盐簇阴离子及其"元件化合物"Roussin红盐簇阴离子,是固氮酶活性中心福州模型I(网兜状原子簇模型)的模型物.本文用闭壳层CNDO/2(S,D方案)法计算了它们的电子结构.根据计算所得的Mulliken重叠集居,电荷密度,分子轨道能量和轨道特征等数据,对成键性质进行了分析,得出如下主要结论:两种簇阴离子骨架电子的非定域性都比较强,桥硫原子Sb在由红盐形成黑盐的电子转移过程中起施主作用,两种簇阴离子中都存在M-M键,强度与M-Sb键相近,其主要贡献都来源于金属的s,pz,dz2轨道与硫原子的s,pz轨道之间的σ作用,金属d轨道的π作用对整个骨架的成键贡献很小.     

铍的夹心和半夹心化合物中铍原子的共价

用DV-XαSCC方法和自然键轨道法研究了CpBeH和Cp_2Be的电子结构和成键情况,进而根据作者之一提出的共价定义研究了CpBeX(X=H,Cl,Br,CH_3,C≡CH,……)和Cp_2Be等化合物中铍原子的共价.结果表明,在这些化合物中,铍原子的共价均为6.     

用电荷和构型自洽的EHMO法计算二茂锆卤化物的电子结构

黄贻琛等曾用未经电荷自洽的 EHMO 法计算了 Cp_2ZrX_2(X=F、Cl、I)的电子结构,得到这类络合物基本上是离子型的结论.我们采用改进的电荷和构型自洽的 EHMO 法对该络合物进行研究,认为它们是具有一定离子性的共价型络合物.计算所得偶极矩与实测值相近.     

超价分子SO2和SO3成键性的研究

SO2,SO3和H2SO4这样的分子,由于其中S原子的"超价”(supervalent)而被称为"超价分子”.对于它们的成键特性,除早期价键理论的共振论观点外,一般都采用简单分子轨道理论观点加以解释,认为在SO2和SO3分子中,诸原子在同一平面上,S原子sp2杂化,与O原子形成σ键,S的另一个p轨道与分子平面垂直,其余O原子也各有一个与分子平面垂直的p轨道,这些p轨道线性组合形成离域π轨道.SO2中的离域π键是Π43,SO3中则是Π64.之后对S原子d轨道是否参与成键的研究表明,S的dπ轨道可以与O的pπ轨道形成π键,认为S原子d轨道参与了π键的形成[1,2].本文通过分子轨道理论的自洽场法从头算和价键理论的广义价键法(Generalized Valence Bond,GVB)计算,根据计算结果对SO2和SO3的成键性做进一步的讨论.......     

分子轨道法研究硅酸盐分子结构

Pauling提出的硅酸盐结构规则,原则上只适用于纯粹离子键化合物。硅酸盐中硅酸根与金属原子之间的化学键是离子键,但是硅氧键却没有证据说明它是离子键。 Stewart等从X射线衍射实验和轨道电负性理论论证,α-石英中Si—O键仅有25％的离子性。还有人认为Si—O键的离子性为31～37.5％。我们用CNDO/2法计算得     

三茚基钐(Sm(C9H7)3)的电子结构和化学键

采用INDO方法研究了三茚基钐的电子结构和化学键。结果表明,钐原子的6s、6p和5d轨道与配位体轨道有不同程度的混合;HOMO和LUMO是由钐原子的5d轨道和配位体原子轨道构成的π型分子轨道;4f轨道是高度定域的,参与成键不超过2%。三茚基钐的化学键具有相当程度的共价性,茚基中碳原子上的净电荷分布不均匀。讨论了三茚基钐的四氢呋喃加合物中的Sm-C键。     

稀土夹心化合物的SCF-X_α-SW研究 Ⅰ.Cp_2Sm、Cp_2Yb和Cp_3Sm

对Cp_2Sm、Cp_2Yb和Cp_3Sm进行了非相对论和相对论SCF-X_α-SW计算,用轨道相互作用、分子轨道图形、布居数分析等方法讨论了化学键图象。在Cp_2Ln(Ln=稀土元素)中以Cp为主要成分的轨道能级两种方案结果变化不大。而相对论间接效应的存在,使Ln4f能级明显升高,与Cp_2Ln易被氧化的实验结果一致。二价的Cp_2Ln成键轨道中Ln成分是d>f>p>s,与三价的Cp_3Sm、LnF_3比较,Ln的s、p、d成分变化不大,而Lnf成分明显减少,使Cp_2Ln共价性明显地低于Cp_3Sm和LnF_3。     

哌嗪二酮的气相Hel光电子能谱及其量化计算

给出了哌嗪二酮的气相HeI紫外光电子能谱(UPS), 并进行了化合物分子的HAM/3, MNDO, MINDO/3, INDO, CNDO/2和EHMO等量子化学计算研究. UPS谱低电离能(<11.00 eV)区的四重峰被指认为分子体系中氧-氧, 氮-氮原子孤对轨道间的通过键相互作用导致的分裂峰. 表明HAM/3和MNDO计算法是预指该化合物实验电离能正确次序、轨道对称性类型以及通过键相互作用导致分裂大小的较好方法.     

镱硫属化合物的密度泛函理论研究

用密度泛函理论(DFT)研究镱硫属化合物的电子结构和性质,通过与实验比较考察了现有的几种近似密度泛函公式对镧系元素化合物的适用程度和相对论效应的影响.结果表明,用DFT计算的YbO键长对实验值的偏差约为0.002nm;但得到的键能即使在考虑梯度校正和相对论效应之后,仍比实验值高,在定域密度近似基础上引入交换梯度校正使键能计算值减小,其中PW86x使键能计算值减小稍多些,结果更接近实验值;相关梯度校正使键能计算值升高.相对论效应使键长缩短0.004～0.006nm,键能减小约0.5eV.计算结果的分析表明,Yb的5d轨道和配体的np轨道间形成σ键和π键.在所研究的分子体系中,配体原子从O到Te、Yb原子的5d轨道布后数依次减少,同键能减弱的顺序一致.相对论效应使键能减小的主要原因是在成键过程中发生了Yb的6s电子向5d轨道的转移,而相对论效应使该过程能量增加.偶极矩和电荷分布的计算表明,Yb-L键以共价性为主,相对论效应使共价性成份增加.     

系列杂氮钛三环中的分子内弱相互作用的紫外光电子能谱和理论研究

首次报道了四种杂氮钛三环的紫外光电子能谱,结合量子化学从头算和和电子密度拓扑分析方法对这些化合物的电子结构和分子内存在的弱相互作用进行了分析研究。通过优化STO-6G(d)基组(C:1s=5.66,ζ2s=1.675,ζ2p=1.66;O:ζ1s=7.67,ζ2s=2.25,ζ2p=2.20;N:ζ1s=6.66,ζ2s=1.91,ζ2p=1.89;H:ζ1s=1.24)对所研究的化合物进行了分子轨道计算,结合计算结果对化合物的紫外光电子能谱进行了解析和指认。实验得到的体现分子内弱相互作用σN-Ti的轨道电子的电离能表明了N,Ti间相互作用的弱成键特性,电子密度拓扑分析显示N,Ti原子间存在体现成键作用的键鞍点,键鞍点处电子密度的拓扑性质也进一步支持了N,Ti间弱成键作用的特性。     

用HAM/3法研究胞嘧啶及其甲基衍生物的光电子能谱与紫外光谱

本文用半经验SCF-MO-HAM/3方法计算了胞嘧啶和它的某些甲基衍生物的电离能、激发能和振子强度.指认了这些分子的紫外光电子能谱和紫外电子光谱.讨论了在紫外光电子能谱指认上与CNDO/S的不同之处.分析了胞嘧啶在磷酸三甲酯中可能存在的主要异构体形式.     

Calculation of excited-state properties of some small molecules: Comparative study of the CNDO/S and CNDO/2-vn−1 potential methods

The transition energy and geometry of the lowest excited (nπ*) singlet and triplet states of CO, CS, HNO, H₂CO, HFCO, and F₂CO molecules are calculated by CNDO /S and CNDO /2-V_N?1 potential methods, and the results are compared with those of experimental and ab initio theoretical studies, wherever available. In the calculation of the vertical transition energy, the performance of the CNDO /S method is seen to be generally more satisfactory than that of the CNDO /2-V_N?1 potential method, while the reverse is true for the excited-state geometry. The CNDO /S method as such fails to describe the geometry of the excited state, but a combined version (CNDO /S-2) of CNDO /S and CNDO /2, as well as the CNDO /2-V_N?1 potential method is fairly successful in this regard.     

Valence shell calculations on polyatomic molecules

Dipole moments and charge distributions for twenty molecules of widely different types have been calculated using (a) the CNDO/2 method and (b) a CNDO/2D method in which the orbitals from the CNDO/2 method are deorthogonalized by a Löwdin transformation and are then used to calculate the dipole moments in a rigorous manner. A statistical analysis of the results for the dipole moments calculated by the CNDO/2D method shows that they are in very slightly better agreement with experiment than those from the CNDO/2 method. The net charge distributions from the CNDO/2D method follow more closely the trends of ab initio calculations than do the CNDO/2 net charges.

---

Zusammenfassung Dipolmomente und Ladungsdichten von Molekülen unterschiedlichen Typs wurden mittels des CNDO/2- und CNDO/2D-Verfahrens (d. i. mit delokalisierten Löwdin-Orbitalen als AO's) berechnet. Eine statistische Analyse zeigt, daß die Resultate der zweiten Methode etwas besser als die der ersten den experimentellen Ergebnissen folgen. Das Analoge gilt für die Nettoladungsverteilungen in bezug auf die Trends bei ab initio-Rechnungen.

---

Résumé Les moments dipolaires et les distributions de charge pour vingt molécules de types divers ont été calculés par: a) la méthode CNDO/2; b) une méthode CNDO/2D où les orbitales de CNDO/2 sont déorthogonalisées par une transformation de Löwdin. Une analyse statistique montre que les moments dipolaires calculés par CNDO/2D sont légèrement en meilleur accord avec l'expérience que ceux calculés par CNDO/2. Les distributions de charge de CNDO/2D sont plus ressemblantes à celles de calculs ab-initio que ne le sont les distributions de CNDO/2.

---

---

This work represents part of the Ph.D. Dissertation submitted to the University of Virginia by D. D. S. and was supported by Grants No. 1-F01-GM41986-01 from the National Institutes of Health, Bethesda, Maryland, U.S.A., and No. AF-AFOSR-1184-67 from the Air Force Directorate of Scientific Research.

---

NASA Research Trainee.     

Deorthogonalization of atomic originals in the CNDO approach

The procedure for deorthogonalization (D) of atomic orbitals in the semiempirical CNDO approach is reviewed. For comparative studies, CNDO/2, CNDO/2D, and STO -3G calculations of molecular dipole moments and Mulliken populations are carried out on 35 prototype molecules containing H, C, N, O, and F atoms. The calculated values are assessed on the basis of how well they agree with experimental trends, chemical bonding theories, and ab initio molecular orbital (MO) values. Results of analyses indicate that the CNDO/2D values for dipole moments are in reasonable agreement with experimental values, and those for net atomic charges and electron populations bear greater resemblance to the ab initio (STO -3G and 6-31G**) values than the original CNDO/2 values. These findings, together with those of previous investigators, demonstrate unequivocally the advantages of incorporating deorthogonalization into routine CNDO/2 or INDO calculations as a means to obtain reasonable estimates of charge distributions.     

Cndo/2 conformational calculation for conjugated and non-conjugated molecular systems

In order to check the effectiveness of the well known CNDO/2 method for predicting molecular conformations, twenty-six molecular systems, involving B, C, N, O, F, Si, S and Cl atoms, were studied. If these molecules are classified into three groups, according to simple rules, the CNDO/2 method is found to fail systematically in one of them. In this group the twisted bond is expected to be delocalized. It is concluded that the CNDO/2 method is not useful in predicting conformations of molecules in which the twisted bond is delocalized.     

Force constants in hydrogen-bonded crystals

Hydrogen bond stretching force constants in crystals of imidazole, urea and cyanuric acid are calculated using both a modified CNDO/2 method and a constrained least squares fitting of interatomic pair potentials to the lattice vibrations. Results show that the modified CNDO/2 method gives closer agreement with the experimentally derived force constants than the normal CNDO/2 parameterization.     

Torsional isomerization of biologically active bicyclic molecules

Rotational isomerization of bipyridines C₅H₄N? C₅H₄N was studied by CNDO /2, PPP /CI , and CNDO /CI methods. It is shown that CNDO /2 overestimates the angle of rotation ? between the pyridine rings ca. two times. The angle ? was determined for 2,2′-bipyridine by means of correlation of the theoretical (CNDO /2) and experimental dipole moment. It was also found from the correlation between the theoretical and experimental UV spectra. It is shown that there is an explicit dependence of the results upon the distance between heteroatoms (PPP /CI ). It has been found that the CNDO /CI method correctly predicts the value of the rotational angles and their sequence in bipyridines.     

Treatment of hydrogen bonding within CNDO/2 and MINDO/3: CNDO/2H and MINDO/3H

A formalism has been developed to treat hydrogen-bonded A—H…?B systems within the CNDO /2 and the MINDO /3 methodologies. In this formalism the interactions are divided into three distinct classes; those between (a) two hydrogen-bonded atoms, (b) one hydrogen-bonded and non-hydrogen-bonded atom, and (c) two non-hydrogen-bonded atoms. The last class of interactions is treated solely by the existing CNDO /2 or MINDO /3 method. For A –H…?B systems, the core resonance integrals are individually parametrized depending upon the class of the interaction. Three types of A—H…?B systems have been thus far parametrized. Nine hydrogen-bonded dimers have been studied using the new formalism and the current CNDO /2 and the MINDO /3 methods. MINDO /3 predicts very large interatomic (A –B) distances for the equilibrium geometry, and relatively small stabilization values for the hydrogen-bond energies. CNDO/2 predicts the reverse. The new formalism for both CNDO /2 and MINDO /3 predicts accurate geometries as well as energies for all nine dimers. The new formalisms are called CNDO /2H and MINDO /3H. A general discussion of the nature of hydrogen bonding as exhibited by CNDO /2H and MINDO /3H is presented.     

A comparative study of the structure and bonding of HOO,HOS, HSO,and HSS radicals by CNDO/2 and INDO methods

Equilibrium geometries, force constants, barriers to linearity, charge distributions, dipole moments, and electron spin density of HOO, HOS, HSO, and HSS radicals are calculated by CNDO/2 and INDO methods using respectively the original and some recently introduced scheme of parametrization. Three sets of calculations, namely, CNDO/2(sp), CNDO/2(spd), and INDO, are performed, and the results are compared with the ab initio and experimental values, wherever available. A good agreement is obtained for geometry in the case of CNDO/2 (sp) and INDO calculations. The performance of CNDO/2 (spd) calculations in this regard is quite unreliable. The stretching force constants are considerably overestimated by all the methods, while the bending force constants are in reasonable agreement with the ab initio values. With respect to dipole moments, the CNDO/2 values are in better agreement with the ab initio results than the INDO values. In all the cases, the dipole moment vector directions are in complete disagreement with the ab initio predictions.     

Calculation of chlorine Kα1,2 X-ray emission shifts in molecules

Shifts of the chlorine Kα_1,2 X-ray emission line in gaseous molecules and in crystalline NaCl are correlated with effective atomic charges from CNDO/2 molecular orbital calculations made (i) for the neutral atomic ground states. and (ii) within the equivalent core approximation. For the latter case CNDO/2 parameters are derived for argon. The best agreement with experimental chemical shifts is obtained with the equivalent cores approximation. The results for NaCl depend strongly upon the CNDO/2 treatment of sodium.