氢原子在Pt及Pt系双金属催化表面吸附的密度泛函理论研究

采用密度泛函理论(dFT)考察了Pt(100)、(110)、(111)三种表面氢原子的吸附行为, 计算了覆盖度为0.25 ML时氢原子在Pt 三种表面和M-Pt(111)双金属(M=Al, Fe, Co, Ni, Cu, Pd)上的最稳定吸附位、表面能以及吸附前后金属表面原子层间弛豫情况. 分析了氢原子在不同双金属表面吸附前后的局域态密度变化以及双金属表面d 带中心偏离费米能级的程度并与氢吸附能进行了关联. 计算结果表明, 在Pt(100), Pt(110)和Pt(111)表面, 氢原子的稳定吸附位分别为桥位、短桥位和fcc 穴位. 三种表面中以Pt(111)的表面能最低, 结构最稳定. 氢原子在不同M-Pt(111)双金属表面上的最稳定吸附位均为fcc 穴位, 其中在Ni-Pt 双金属表面的吸附能最低, Co-Pt 次之. 表明氢原子在Ni-Pt 和Co-Pt 双金属表面的吸附最稳定. 通过对氢原子在M-Pt(111)双金属表面吸附前后的局域态密度变化的分析, 验证了氢原子吸附能计算结果的准确性. 掺杂金属Ni、Co、Fe 的3d-Pt(111)双金属表面在吸附氢原子后发生弛豫, 第一层和第二层金属原子均不同程度地向外膨胀. 此外, 3d金属的掺入使得其对应的M-Pt(111)双金属表面d带中心与Pt 相比更靠近费米能级, 吸附氢原子能力增强, 表明3d-Pt系双金属表面有可能比Pt具有更好的脱氢活性.     

氢原子在Ni(111)次表面和Ni(211)、(533)台阶面上的吸附与振动

应用原子与表面簇合物相互作用的五参数Morse势(5-MP)方法对氢原子在Ni(111)表面和次表面以及Ni(211), (533)台阶面进行了系统研究, 得到了氢原子在上述各面的吸附位、吸附几何、结合能和本征振动频率. 计算结果表明, 在Ni(111)面上, 氢原子优先吸附在三重位, 随着覆盖度的增加会吸附在次表面八面体位和四面体位. Ni(211), (533)的最优先吸附位都是四重位, 当氢原子的覆盖度增大时占据(111)平台的三重吸附位. 靠近台阶面的吸附位受台阶和平台高度的影响很大. 此外, 我们计算了氢原子在各表面的不同吸附位的扩散势垒, 获得氢原子在各表面的最低能量扩散通道.     

氢原子在Pd(997)台阶面上的吸附和扩散

用对势方法研究了氢原子在Pd（997）台阶面上的吸附和扩散．计算结果表明，氢原子在台阶下形成最稳定的吸附态，并分析了其新的振动模式．理论分析表明，台阶对下台面上扩散的氢原子形成捕获势，对上台面扩散的氢原子形成反射势，很好地支持了实验结果．     

金电极析氢反应过程中吸附氢原子的表面扩散及复合反应动力学研究

用特别设计的以平面金电极为工作电极,以Nafion膜为固体电解质的实验装置,对金电极析氢反应过程中吸附氢原子的表面扩散及吸附氢原子复合为氢分子的反应进行了定量研究,发现吸附氢原子在金表面上的扩散及复合反应引起了电流增加.对这一随时间变化的电流增加数据进行计算机模拟和数值拟合,得到金电极析氢反应中吸附氢原子的表面扩散系数、复合反应常数及其它动力学参数,并对数据进行了讨论.     

氢在Mg_2Ni(100)面的吸附及扩散

运用第一性原理研究氢在清洁和掺杂Al的Mg2Ni(100)面的吸附及扩散.在清洁Mg2Ni(100)表面,氢原子可稳定地吸附于Mg-Ni桥位和Mg-Mg桥位,吸附能为1.19-1.52eV.在掺杂Al的Mg2Ni(100)表面,氢原子可稳定吸附于Al-Ni、Mg-Ni、Mg-Al桥位,吸附能为0.10-0.29eV.氢在掺杂Al的Mg2Ni(100)表面的吸附能低于其在清洁表面的吸附能,说明掺杂Al后氢原子与表面的相互作用减弱.过渡态计算结果表明,氢原子由清洁的Mg2Ni(100)面及掺杂Al的Mg2Ni(100)面扩散至次表层的势垒分别为0.59及-0.04eV,掺杂Al后氢原子的扩散势垒降低,说明氢原子更易由掺杂Al表面扩散至次表层.Al原子替代Mg2Ni(100)面的Mg原子减弱氢原子与表面的相互作用,降低氢原子由表层扩散至次表层的势垒,这可能是Mg2Ni合金掺杂Al可改善其吸氢动力学性能的主要原因之一.     

氢原子在平坦金属锂(100)面上表面扩散行为的ab initio研究

本文用Li7(4,3)-H和Li9(5,4)-H小原子簇模拟氢原子在平坦金属锂(100)面吸附体系, 取小基组作了各吸附位吸附势能曲线及相应分子轨道能级图、吸附和表面扩散势能面的ab initio研究。结果表明, 氢原子优先吸附在配位数较高的吸附位上, 并倾向于由高配位数吸附位向低配位数吸附位迁移, 表面扩散各向异性, 扩散跳跃步长与锂单晶晶格原子间距数量级相同。从吸附和表面扩散势能确定了最低能量表面扩散途径, 分析了原子在平坦金属表面上迁移的微观过程。     

原子H在Cu(100)(111)(110)上的吸附扩散研究

采用5-MP势方法,对原子氢在金属Cu的3个低指数面上的吸附特性,如吸附几何、吸附能、振动频率等以及吸附扩散势能面结构进行了比较系统的研究,计算结果显示低温低覆盖条件下,氢原子在Cu(110)表面上只存在赝式三重位和长桥位吸附态,没有短桥位吸附态,并且获得了实验和理论的支持.     

氢原子在金属锂表面台阶附近吸附和表面扩散行为的ab initio研究

取Li7H和Li9H两个原子簇模拟氢原子与含台阶的金属锂表面的相互作用, 以小基组用ab initip方法计算了体系的吸附和表面扩散势能面(或势能曲线)。结果表明: (1)对Li7H体系, 台阶面附近沿垂直边棱方向存在三种不同的桥位吸附位, 最稳定的吸附位在上台面接近台阶边棱处, 台阶面显著地改变了表面扩散活化能, 台阶边棱处有一个较高的势垒。于是, 迁移原子将会在台阶边棱处受到反射, 并可被捕获于台阶面上及其附近。由势能面确定了最低能量表面扩散途径。(2)对Li9H体系, 在Li7H原子簇基础上增加次表面层两个锂原子后, 表面扩散活化能略有减小, 氢原子在上台面的桥位吸附更趋稳定, 各吸附位相对稳定性及势垒内何位置几无改变, 这些结果显示了台阶面对氢原子的化学吸附和表面扩散发生扰动, 台阶边棱对表面扩散起着重要作用。     

氢原子在金属锂表面台阶附近吸附和表面扩散行为的ab initio研究

取Li_7H和Li_9H两个原子簇模拟氢原子与含台阶的金属锂表面的相互作用,以小基组用abinitio方法计算了体系的吸附和表面扩散势能面(或势能曲线),结果表明:(1)对Li_7H体系,台阶面附近沿垂直于边棱方向存在三种不同的桥位吸附位,最稳定的吸附位在上台面接近台阶边棱处,台阶面显著地改变了表面扩散活化能,台阶边棱处有一个较高的势垒,于是,迁移原子将会在台阶边棱处受到反射,并可被捕获于台阶面上及其附近,由势能面确定了最低能量表面扩散途径,(2)对Li_9H体系,在Li_7H原子簇基础上增加次表面层两个锂原子后,表面扩散活化能略有减小,氢原子在上台面的桥位吸附更趋稳定,各吸附位相对稳定性及势垒几何位置几无改变,这些结果显示了台阶面对氢原子的化学吸附和表面扩散发生扰动,台阶边棱对表面扩散起着重要作用。     

氢原子在钯低指数表面上的吸附和扩散

采用氢－钯相互作用的五参数Ｍｏｒｓｅ势，用对势方法研究了氢原子在Ｐｄ（１００）、Ｐｄ（１１１）和Ｐｄ（１１０）低指数平坦表面上的吸附和扩散，得到氢原子在三个表面上的吸附位、吸附几何、结合能和本征振动等数据，计算结果和实验结果符合得很好。在此基础上，系统地研究了三个系统的吸附扩散势能面结构。     

分形表面上CO氧化反应的蒙特卡罗模拟

众所周知,人们在理想完好的单晶表面上的动力学研究已取得了很大的成功[1],然而对许多实际催化剂来说,传统催化理论和实验之间不符的情况仍然俯拾皆是．现在人们对造成这种情况的原因已经有了更加深入的认识：大量的实验事实表明,在表面催化反应体系中,不但催化剂表面具有复杂的分形结构,而且催化剂表面上的活性中心分布也具有复杂的分形特性[2]．“分形”是指那些具有分数维数的几何对象．这些对象往往是不规则的,不能用通常的欧氏几何来描述[3]．将分形引入催化科学中最早的是Pfeifer和Avnir等[4]人,至今已有近二十年的历史,现…     

Diffusion of particles over dynamically disordered lattice

Using kinetic Monte Carlo simulations we have investigated the diffusion of particles over a lattice with moving atoms in the framework of a simple lattice-gas model of a reconstructive surface. The particle migration over the static and dynamic lattices differs considerably. The dynamic lattice reconstruction changes substantially the particle diffusion coefficient. The Arrhenius dependencies are quantitatively different. An interesting peculiarity of the particle migration is the formation of defectons--local lattice deformations around the adsorbed particles. In certain ranges of the system parameters (jump rates of the substrate atoms and adsorbed particles) the adsorbed particles cause local displacements of the nearest substrate atoms, resulting in deeper adsorption sites and lower jump rates. Such particle self-trapping produces a characteristic minimum on the Arrhenius dependencies.     

O2和CO表面催化反应的活性位分布模型

对于O_2和CO表面催化反应，建立了一个新的不可逆Monte Carlo模拟模型。在 二维格子中，引进了表面活性位和非活性位的要领。模型假设，一定浓度的活性位 随机分布在非活性位上，形成了活性位分布的二维格子模型反应器，并在ZGB模型 的基础上，考虑了氧原子和CO分子的表面扩散，特别是引进了吸附粒子的定向表面 扩散。其中，活性位和活性位最近邻是表面吸附物质反应的活性中心，而非活性位 的作用是通过表面扩散传质。当活性位浓度C_a = 1且考虑扩散时，模型还原为增 加了扩散的ZGB模型。当活性位浓度C_a = 1且只考虑氧的扩散时，模拟结果表明， 扩散几率达到某一数值（0.3）时，二级相变点完全消失。当活性位浓度C_a逐渐减 小时，单位活性位产生的CO_2的速率不断增大，这表明活性位的利用率提高了。     

Kinetic Monte Carlo simulations of the interaction of oxygen with Pt(111)

The adsorption, dissociation, diffusion, and desorption of oxygen interacting with the Pt(111) surface have been studied using kinetic Monte Carlo simulations. This study has been motivated by uncertainties in the theoretical and the experimental derivations of O(2)Pt(111) reaction barriers. The simulations reproduce all known experimental data within basically one set of parameters, thus yielding microscopic insights into the elementary reaction steps occurring in the interaction of oxygen with Pt(111) and providing reliable estimates for adsorption energies and diffusion and desorption barriers. In particular, we confirm that the distance of oxygen atoms directly after dissociation is caused by ballistic hot atom motion rather than by diffusive motion. We address the equilibrium structure of oxygen atoms at high coverages. At low temperatures, chains of oxygen pairs are formed. We show that this mechanism can be explained by a lowered dissociation in the vicinity of already adsorbed atoms. Finally we discuss the role of the lateral interaction between the oxygen atoms in the oxygen desorption process.     

Diffusion Monte Carlo simulations on uracil-water using an anisotropic atom-atom potential model

We have developed an anisotropic atom-atom intermolecular potential model for the interaction of uracil with water. The potential consists of a distributed multipole analysis (DMA) model for the electrostatic energy, and a 6-exp potential to represent the repulsion-dispersion term. The repulsion-dispersion potential parameters are adjusted to yield good agreement with accurate ab initio data on the minima and transition states of the uracil-water complex. We have used this potential in diffusion Monte Carlo simulations of uracil-water, uracil-(water)2 and uracil-(water)3. The uracil-water simulations show that the theoretically based potential gives a qualitatively different picture of uracil hydration than that provided by a standard isotropic atom-atom point charge model, which is shown to underestimate the delocalized motion of the water hydrogen atoms. Plots of the vibrational probability density of the hydrogen atoms show the delocalized motion of the water hydrogen atoms that are not involved in hydrogen bonding.     

Self-diffusion of methane in single-walled carbon nanotubes at sub- and supercritical conditions

The diffusivities of methane in single-walled carbon nanotubes (SWNTs) are investigated at various temperatures and pressures using classical molecular dynamics (MD) simulations complemented with grand canonical Monte Carlo (GCMC) simulations. The carbon atoms at the nanotubes are structured according to the (m, m) armchair arrangement and the interactions between each methane molecule and all atoms of the confining surface are explicitly considered. It is found that the parallel self-diffusion coefficient of methane in an infinitely long, defect-free SWNT decreases dramatically as the temperature falls, especially at subcritical temperatures and high loading of gas molecules when the adsorbed gas forms a solidlike structure. With the increase in pressure, the diffusion coefficient first declines rapidly and then exhibits a nonmonotonic behavior due to the layering transitions of the adsorbed gas molecules as seen in the equilibrium density profiles. At a subcritical temperature, the diffusion of methane in a fully loaded SWNT follows a solidlike behavior, and the value of the diffusion coefficient varies drastically with the nanotube diameter. At a supercritical temperature, however, the diffusion coefficient at high pressure reaches a plateau, with the limiting value essentially independent of the nanotube size. For SWNTs with the radius larger than approximately 2 nm, capillary condensation occurs when the temperature is sufficiently low, following the layer-by-layer adsorption of gas molecules on the nanotube surface. For SWNTs with a diameter less than about 2 nm, no condensation is observed because the system becomes essentially one-dimensional.     

Hydrogen transfer reaction on the surface of an oxide catalyst

RuO2(110) exposes two kinds of active surface species (acidic and basic centers) that govern the interaction of the gas phase in contact with the catalyst's surface. Here we will elucidate the cooperative interplay of these two active surface sites for a simple model reaction, namely the water formation over RuO2 catalysts when supplying hydrogen and oxygen from the gas phase. The bridging O atoms harvest the hydrogen from the gas phase, while the on-top O atoms pick up those adsorbed hydrogen atoms from the bridging O atoms to form water. This mechanism of hydrogen transfer is mediated by a strong hydrogen bond. Hydrogen transfer is expected to play a vital role for the whole class of catalyzed hydrogenation and dehydrogenation reactions of hydrocarbons over RuO2.     

Kinetic study of the "surface explosion" phenomenon in the NO+CO reaction on Pt(100) through dynamic Monte Carlo simulation

The extremely narrow production peak of N2 and CO2 which occurs in the reaction of NO+CO on Pt(100), a phenomenon known as "surface explosion," is studied using a dynamic Monte Carlo method on a square lattice at low pressure under isothermal conditions. This analysis incorporates recent experimental evidence obtained for the same reaction on a Rh(111) surface, which has shown that N2 production occurs either from the classical N+N recombination step or by the formation and successive decay of a (N-NO)* intermediary species. Moreover, the NO dissociation rate is inhibited by coadsorbed NO and CO molecules and is enhanced both by the presence of empty sites and adsorbed N atoms as nearest neighbors. These effects are taken into account in this study, along with the experimental adsorption, desorption, and diffusion rates of the reactants. The "explosive" phenomenon is analyzed through the evolution over time of an adsorbed NO+CO monolayer at a fixed temperature of 400 K. Furthermore, as the diffusion processes of the adsorbates are included, cellular structures are observed. Our simulations show quantitative agreement in the position of maxima with those obtained through experiments using isothermal desorption mass spectroscopy.     

Kinetic Monte Carlo studies of hydrogen abstraction from graphite

We present Monte Carlo simulations on Eley-Rideal abstraction reactions of atomic hydrogen chemisorbed on graphite. The results are obtained via a hybrid approach where energy barriers derived from density functional theory calculations are used as input to Monte Carlo simulations. By comparing with experimental data, we discriminate between contributions from different Eley-Rideal mechanisms. A combination of two different mechanisms yields good quantitative and qualitative agreement between the experimentally derived and the simulated Eley-Rideal abstraction cross sections and surface configurations. These two mechanisms include a direct Eley-Rideal reaction with fast diffusing H atoms and a dimer mediated Eley-Rideal mechanism with increased cross section at low coverage. Such a dimer mediated Eley-Rideal mechanism has not previously been proposed and serves as an alternative explanation to the steering behavior often given as the cause of the coverage dependence observed in Eley-Rideal reaction cross sections.     

Surface diffusion of metal atoms on polymer substrates during physical vapour deposition

 The surface diffusion of physical-vapour-deposited metal atoms on thermoplastic polymer substrates was investigated. In accordance with the hypothesis of the “classical” atomistic diffusion model, diffusion coefficients are derived from a Monte Carlo simulation. Because the “classical” atomistic diffusion model neglects the desorption of the metal atoms, the absolute diffusion data obtained in our investigations should only be considered as rough estimates. It is more the intention of our work to present relative values in order to correlate the metal surface diffusion on polymer substrates with their physical states (morphologies and surface dynamics). As expected, the diffusivity of metal atoms is strongly influenced by the chemical affinity (“reactivity”) between the metal atoms and the polymer substrate. Furthermore, the diffusivity strongly depends on the physical state of the polymer substrate. On polymer surfaces above the glass-transition temperature the surface diffusivity of metal atoms is 1 order of magnitude higher than the diffusivity below the glass-transition temperature. Received: 9 April 1999/Accepted in revised form: 21 October 1999