金属氧化物SPC簇模型方法——嵌入簇点电量的确定

探讨和比较了确定金属氧化物晶体点电荷量的几种自洽条件的优劣 ,并使用abinitioHF方法对在不同点电荷量嵌入环境下MgO簇电子性质进行了研究 .计算结果表明 :点电荷量的大小会明显影响嵌入簇的电子性质 ;即使MgO被通常认为是纯的离子型氧化物 ,点电荷量取为表观化合价± 2 ,仍过高地估计了Madelung势 ,嵌入簇与点电荷之间应满足一定的自洽条件 ;使用点电荷为嵌入环境时 ,即使达到电荷自洽、偶极矩自洽或势自洽 ,簇模型计算结果仍与固体性质存在一定的偏差 ;将点电荷进行球谐展开 ,赋予其连续变化的电荷密度分布后 ,得到了一个对固体整体性质描述较好的SPC簇模型 .     

O2在具有氧和镁缺陷MgO(001)表面的吸附

在密度泛函理论的框架下，采用嵌入点电荷簇模型研究了O2在具有氧缺陷和镁缺陷MgO(001)表面上的吸附.用电荷自洽的方法确定了点电荷的值.计算结果表明，O2倾向吸附在具有氧缺陷的MgO(001)表面上.通过和我们近期研究过的O2在低配位的边、角上吸附结果相比较，发现具有氧缺陷的MgO(001)表面更加有利于O2的吸附和解离. Mülliken电荷分析表明，电荷由底物向吸附的O2反键轨道上转移是导致O2键强削弱的主要原因.势能曲线表明，O2在具有氧缺陷的MgO(001)表面上发生解离所需要克服的能垒比在角阳离子端发生解离所需克服的能垒有大幅度降低.     

簇模型ZnO(100)面上吸附CO_2的密度泛函理论研究

采用电荷自洽方法, 以嵌入原子簇Zn4O4为模型, 使用量子化学的密度泛函理论, 研究了二氧化碳在六方ZnO非极化的(1010)面的可能吸附态。计算表明, CO2垂直底物表面吸附, 氧原子只能与Zn原子配位, 并且吸附能为很弱的1.8 kJ/mol;吸附质分子平行于底物表面时, 得到了5种平衡吸附构型, 其中采用CZn配位和η2O, O二齿配位时, 吸附很弱, 经BSSE校正后的吸附能在8.8~6.6 kJ/mol。 采用η2C, O方式分别与O和Zn配位时, 吸附能为31.1 kJ/mol; C原子与表面O配位时计算得到了唯一的一个化学吸附态, 吸附能为139.6 kJ/mol, 与实验结果一致。     

簇模型ZnO(10(1-)0)面上吸附CO2的密度泛函理论研究

采用电荷自洽方法，以嵌入原子簇Zn4O4为模型，使用量子化学的密度泛函理论，研究了二氧化碳在六方ZnO非极化的(101^-0)面的可能吸附态。计算表明，CO2垂直底物表面吸附，氧原子只能与Zn原子配位，并且吸附能为很弱的1．8kJ／mol；吸附质分子平行于底物表面时，得到了5种平衡吸附构型，其中采用C-Zn配位和η^2-O,O二齿配位时，吸附很弱，经BSSE校正后的吸附能在8．8～6．6kJ／mol。采用η^2-C,O方式分别与O和Zn配位时，吸附能为31．1kJ／mol；C原子与表面O配位时计算得到了唯一的一个化学吸附态，吸附能为139．6kJ/mol,与实验结果一致。     

金属氧化物表面化学吸附和反应的量子化学簇模型方法研究

综述了本研究小组利用量子化学簇模型方法研究金属氧化物表面化学吸附和反应的工作.提出了选簇的三个原则,即电中性原则、化学配比原则和配位原则.发现在符合前两个原则的基础上,一个具有最饱和配位、或最少悬空键的簇往往是一个用于化学吸附研究的好的簇模型.与此同时,探讨了如何恰当地考虑大块固体本底的长程影响,提出了用球电荷模拟簇模型的环境、环境与簇体进行电荷自洽的SPC簇模型方法.利用该模型研究了一系列具有催化背景的重要体系,包括H2/ZnO、O/MgO、NO/MgO、N2O/MgO、N2O/Li/MgO、CO/MgO、CO/NiO等.     

CO在M55（M=Cu,Ag,Au）团簇上吸附的密度泛函研究

在全电子相对论BVP86/DNP水平下对CO在Au55,Ag55和Cu55团簇上的吸附进行了比较研究,并考察了电荷对吸附的影响.计算结果表明,CO在Au55团簇上吸附能最大,其次为Cu55团簇,最弱的为Ag55团簇.团簇电荷对C—O键活化和CO与团簇表面原子成键影响较小.金团簇的电荷对吸附能影响较大,而银和铜团簇的电荷对吸附能影响较小.CO吸附到团簇上导致团簇上电子向CO转移.C—O键活化强度与吸附位置密切相关,其中孔位吸附导致C—O键活化程度最大,最弱的为顶位吸附.CO在金团簇上吸附具有较好选择性,而在银和铜团簇上吸附无选择性.     

MoO2储锂性能及循环容量反常特性的第一性原理研究

基于第一性原理,对MoO2作为电极材料的储锂性能进行了计算,并探讨了其储锂容量在一定循环次数内呈上升的反常现象微观机理.计算了MoO2材料中Li的单键能,态密度(DOS)及其嵌锂电压,结果表明MoO2中Li的吸附能较大,储锂结构稳定.嵌锂结构呈金属性,嵌锂电压变化规律与文献实验结果一致.针对循环容量反常特性,计算了Mo的空位形成能,LiMoO2的差分电荷密度以及电荷布居情况,计算结果表明Li的嵌入能为O提供电荷,减弱了Mo—O键间的相互作用,另一方面嵌入的Li能减弱Mo空位形成后的电荷极化作用,从而大大降低Mo空位的形成能.形成的Mo空位能为Li的嵌入提供了新的吸附位点,提高了嵌锂的容量.计算结果与实验符合得很好,能为电极材料储锂性能的改善提供一定的理论指导.     

Pt原子在y-Al203(001)表面的吸附及迁移

采用密度泛函理论(DFT)中广义梯度近似(GGA)方法,对Pt原子与y-Al2O3(001)面的相互作用及迁移性能进行了研究.分析了各种可能吸附位及吸附构型的松弛和变形现象,吸附能和迁移能垒的计算结果表明:Pt团簇能够稳定吸附在该表面.Pt原子在表面O位的吸附能明显较高,这主要是由Pt向基底O原子转移了电子所致.电荷布居分析表明,Pt原子显电正性,Pt和Al原子之间存在排斥作用,导致与Al原子产生较弱相互作用.计算的平均吸附能大小依赖于Pt团簇的大小和形状,总体趋势是随着Pt原子数增多,吸附能降低.Pt原子在y-Al2O3(001)表面迁移过程所需克服的迁移能垒最高值为0.51 eV.随着吸附的Pt原子数增多,更倾向于形成Pt团簇.因此,Pt原子在y-Al2O,(001)表面的吸附演变不可能形成光滑、均匀平铺的吸附构型,而在一定条件下容易出现团聚.     

第一性原理模拟氢气分子在氧化硅团簇上的吸附

基于密度泛函理论,采用广义梯度近似(GGA)分析了H2分子吸附在氧化硅团簇上的几何结构、电子性质以及吸附能.结果发现:H2分子与Si3O4团簇相互作用时,H2分子被分解,游离的H原子优先吸附在末端Si原子上,表明Si3O4团簇体系对氢气的存储主要依赖于末端存在悬挂键的Si原子,接着H2分子才以分子的形式以较小吸附能吸附在Si3O4H4团簇上.氢气分子主要引起与其邻近的原子电荷的重新分布.该团簇体系的红外、拉曼光谱图有效地鉴定了H2分子的吸附状态,为理论上确定团簇的稳定结构和实验上对观测结果的分析提供有力的途径.     

NiO(001)表面吸附CO的从头算研究

用量子化学B3LYP方法,以外加点电荷来封闭边界效应的簇为模型,计算了CO在NiO(001)面上不同吸附位置的吸附情况,并计算了振动频率。结果表明:1) CO的最佳吸附方式为C端垂直吸附在Ni位;2)吸附后CO间的振动频率蓝移13 cm-1;3)在O空缺、边和角等位置的吸附不如在完整表面的吸附稳定,这些均与实验结果一致。吸附后CO把主要起反键作用的C2s电子给予簇表面,使得吸附后CO键级加大,导致吸附后振动频率蓝移。并比较了Gaussian98 和Crystal98的计算结果,两者的结果能较好地符合。     

DFT Study of Metal Atoms Adsorbed at Low-coordinated Sites of MgO （001） Surface

1 INTRODUCTION The interfaces between metals and oxide play a vital role in many industrial applications: hetero- geneous catalysis, microelectronics, thermal barriers, corrosion protection, metal processing and so on[1]. In catalysis, the choice of metal and oxide support is critical in order to obtain a desired reactivity and selectivity[2]. This is due in part to the inherent reac- tivity of the two components. Also the size and shape of the metal particle, which depend on the choice…     

NO在MgO(001)缺陷表面吸附解离的量子化学研究

在密度泛函理论框架下,采用嵌入点电荷簇模型研究了NO在MgO(001)完整和缺陷表面上的吸附。研究结果表明:具有氧缺陷结构表面的催化活性较高,有利于NO键的削弱;当另一个NO分子进攻已吸附的NO分子时,NO键将进一步削弱,直致断裂,并伴有N2O产生,这与UPS和MIES实验观察到的现象一致。Mulliken布居分析指出,底物电子向NO转移,并填充到NO的*反键轨道上,从而导致NO键的削弱,并形成NO-。这也是可能导致形成NO-的原因。研究还表明,具有镁缺陷的MgO(001)表面对NO的解离没有催化活性。     

Conversion of N2O to N2 on MgO （001） Surface with Vacancy： A DFT Study

The adsorption and decomposition of NzO at regular and defect sites of MgO (001) surface have been studied using cluster models embedded in a large array of point charges (PCs) by DFT/B3LYP method. The results indicate that the MgO (001)surface with oxygen vacancies exhibits high catalytic reactivity toward N2O adsorptive-decomposition. It is different from the regular MgO surface or the surface with magnesium vacancies.Much elongation of O—N bond of N2O after adsorption at oxy-gen vacancy site with O end down shows that O—N bond has been broken with concurrent production of N2, leaving a regu-lar site instead of the original oxygen vacancy site (F center ).The MgO (001) surface with magnesium vacancies hardly ex-hibits catalytic reactivity. It can be concluded that N2O dissoci-ation likely occurs at oxygen vacancy sites of MgO (001) sur-face, which is consistent with the generally accepted viewpoint in the experiments. The potential energy surface (PES) reflects that the dissociation process of N2O does not virtually need to surmount a given energy barrier.     

FTIR spectroscopy combined with quantum chemical calculations to investigate adsorbed nitrate on aluminium oxide surfaces in the presence and absence of co-adsorbed water

Surface reactions of nitrogen oxides with aluminium oxide particles result in the formation of adsorbed nitrate. Specifically, when alpha-Al(2)O(3) and gamma-Al(2)O(3) particles are exposed to gas-phase NO(2) and HNO(3) adsorbed nitrate forms on the surface. In this study, Fourier transform infrared (FTIR) spectroscopy is combined with quantum chemical calculations to further our understanding of the adsorbed nitrate product on aluminium oxide particle surfaces in the presence and absence of co-adsorbed water at 296 K. FTIR spectra of adsorbed nitrate on alpha-Al(2)O(3) and gamma-Al(2)O(3) particles are interpreted using calculated vibrational frequencies of nitrate coordinated to binuclear Al oxide cluster models. Comparison of the calculated and experimental vibrational frequencies of adsorbed nitrate establishes different modes of coordination (monodentate, bidentate and bridging) of the nitrate ion to the surface in the absence of adsorbed water. In the presence of co-adsorbed water, the nitrate ion becomes fully solvated, as shown by a comparison of the experimental nitrate infrared spectra as a function of relative humidity with the calculated nitrate vibrational frequencies for binuclear Al cluster compounds which contain both coordinated nitrate ions and water molecules. These calculations also suggest that adsorbed water can displace nitrate from direct coordination to the surface, leading to an outer-sphere nitrate adsorption complex as well as an inner-sphere complex. Furthermore, the relative humidity dependence of the spectra suggest that water does not evenly wet the surface even at high relative humidity, as there are open or bare surface sites where nitrate ions are not solvated. Besides adsorbed mondendate, bidendate, bridging and solvated nitrate, the presence of ion bound nitrate ion, partially solvated nitrate, molecular nitric acid, hydronium ion and H(3)O(+):NO(3)(-) ion pairs on the oxide surface are also discussed.     

Adsorption and diffusion of Mg, O, and O2 on the MgO(001) flat surface

A thorough investigation of the adsorption and diffusion of Mg, O, and O(2) on MgO(001) terraces is performed by first-principles calculations. The single Mg adatom weakly binds to surface oxygens, diffuses, and evaporates easily at room temperatures. Atomic O strongly binds to surface oxygens, forming peroxide groups. The diffusion of the O adatom is strongly influenced by the spin polarization, since energy barriers are significantly different for the singlet and triplet states. The crossing of the two Born-Oppenheimer surfaces corresponding to the distinct spin states is also analyzed. Although the O(2) molecule does not stick to the perfect surface, it chemisorbs on surface nonstoichiometric point defects such as O vacancies or Mg adatoms, forming in the latter case new chemical species on the surface. We show that the oxidation rate limiting factor in an O(2) atmosphere is the concentration of point defects (O vacancies and Mg adatoms) in the growing surface. The simulated O core-level shifts for the various adsorption configurations enable a meaningful comparison with the measured values, suggesting the presence of peroxide ions on growing surfaces. Finally, the computed energy barriers are used to estimate the Mg and O surface lifetimes and diffusion lengths, and some implications for the homoepitaxial growth of MgO are discussed.     

Theoretical study of MgO(001) surfaces: Pure,doped with Fe,Ca, and Al,and with and without adsorbed water

Ab initio calculations of large cluster models have been performed in order to study water adsorption at the five‐fold coordinated adsorption site on pure Mg(001) and MgO(001) surfaces doped with Fe, Ca, and Al. The geometric parameters of the adsorbed water molecule have been optimized preparatory to analysis of binding energies, charge transfer, preferential sites of interaction, and bonding distances. We have used Mulliken population analysis methods in order to analyze charge distributions and the direction of charge transfer. We have also investigated energy gaps, HOMO energies, and SCF orbital energies as well as the acid‐base properties of our cluster model. Numerical results are compared, where possible, with experiment and interpreted in the framework of various analytical models. © 2001 John Wiley & Sons, Inc. Int J Quant Chem, 2001     

EXAFS and DFT studies of microscopic structure with different density upon Zn(II) adsorption on anatase TiO2

Microscopic structures of Zn(II) adsorbed on anatase TiO₂ surface with different densities were studied using extended X-ray absorption fine structure (EXAFS) spectroscopy and density functional theory (DFT) calculation. Quantitative analysis of the EXAFS spectra showed that microscopic structures of Zn(II) were fourfold coordinated complexes, and different microscopic structures were present of the solid surface. Three modes of corner–corner/sharing-corner/sharing-edge adsorptions on anatase (101) face cluster were calculated by the DFT method. The results from DFT method were consistent with the EXAFS fittings. The optimized Zn–O average distance of the Zn–O tetrahedron was determined as about 2.00 Å. The Zn–Ti distance was 3.69 Å for the corner–corner adsorption, 3.35 Å for the sharing-corner adsorption, and 3.02 Å for the sharing-edge adsorption. According to the adsorption energies calculated by the DFT method, the microscopic structure of corner–corner adsorption was less stable than those of the other adsorption modes. With the increasing adsorption density, the corner–corner adsorption mode would be enhanced more intensively than the other adsorption modes.     

Method of local increments for the calculation of adsorption energies of atoms and small molecules on solid surfaces. 2. CO/MgO(001)

The method of local increments is used in connection with an embedded cluster approach and wave function based quantum chemical ab initio methods to describe the adsorption of a single CO molecule on the MgO(001) surface. The first step in this approach is a conventional Hartree-Fock calculation. The occupied orbitals are then localized by means of the Foster-Boys localization procedure, and the full system is decomposed into several "subunits" that consist of the orbitals localized at the CO molecule and at the Mg and O atoms of the MgO cluster. The correlation energy is expanded into a series of local n-body increments that are evaluated separately and independently. In this way, big savings in computer time can be achieved because (a) the treatment of a large system is replaced with a series of much faster calculations for small subsystems and (b) the big basis sets necessary for describing dispersion effects are only needed for the atoms in the respective subsystem while all other atoms can be treated by medium size Hartree-Fock type basis sets. The coupled electron pair approach, CEPA, an approximate coupled cluster method, is used to calculate the correlation energies of the various subsystems. For the vertical adsorption of CO on top a Mg atom of the MgO(001) surface with the C atom toward Mg, the individual one- and two-body increments are calculated as functions of the CO-MgO separation and a full potential energy curve is constructed from them. A very shallow minimum with an adsorption energy of 0.016 eV at a Mg-C distance of 3.04 ? is found at the Hartree-Fock level, while inclusion of correlation (dispersion) effects shortens the Mg-C distance to 2.59 ? and yields a much larger adsorption energy of 0.124 eV. This is in very good agreement with the best experimental value of 0.14 eV. The basis set superposition error, BSSE, was fully corrected for by the counterpoise method and the bonding mechanism was analyzed at the Hartree-Fock level by means of the constrained space orbital variation, CSOV, analysis.