氢原子在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(100), Ni(111)和Ni(110)面上吸附扩散势能面的结构

本文构造了氢-镍相互作用的5参数Morse势, 用经典的对势方法研究氢原子在Ni(100), Ni(111)和Ni(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原子簇基础上增加次表面层两个锂原子后,表面扩散活化能略有减小,氢原子在上台面的桥位吸附更趋稳定,各吸附位相对稳定性及势垒几何位置几无改变,这些结果显示了台阶面对氢原子的化学吸附和表面扩散发生扰动,台阶边棱对表面扩散起着重要作用。     

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

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

二维TiC层表面H2的化学吸附与物理吸附

第一性原理计算研究发现由于二维TiC单原子层具有高的比表面积与大量的暴露在表面的Ti原子,其是一种非常有潜力的储氢材料.计算结果显示H2可以在二维TiC单原子层表面进行物理吸附与化学吸附.其中化学吸附能为每个氢分子0.36 eV,物理吸附能是每个氢分子0.09 eV.覆盖度为1和1/4层(ML)时,H2分子在二维TiC单原子层表面的离解势垒分别为1.12和0.33 eV.因此,除了物理吸附与化学吸附,TiC表面还存在H单原子吸附.最大的H2储存率可以达到7.69％(质量分数).其中,离解的H原子、化学吸附的H2、物理吸附的H2的储存率分别为1.54％、3.07％、3.07％.符合Kubas吸附特征的储存率为3.07％.化学吸附能随覆盖度的变化非常小,这有利于H2分子的吸附与释放.     

氢原子在Be(0001)表面吸附的密度泛函理论研究

采用第一性原理的密度泛函理论研究单个氢原子和多个氢原子在Be(0001)表面吸附性质.给出了氢吸附Be(0001)薄膜表面的原子结构、吸附能、饱和度、功函数、偶极修正等特性参数.同时也讨论了相关吸附性质与氢原子覆盖度(0.06-1.33ML)的关系.计算结果表明:氢原子的吸附位置与覆盖度之间有强烈的依赖关系,覆盖度低于0.67ML时,氢原子能量上易于占据fcc或hcp的中空位置;覆盖度为0.78ML时,中空位与桥位为氢原子的最佳吸附位;覆盖度在0.89到1.00ML时,桥位是氢原子吸附能量最有利的位置;以上覆盖度中Be(0001)表面最外层铍原子的结构均没有发生明显变化.当覆盖度为1.11-1.33ML,高覆盖度下Be(0001)表面的最外层铍原子部分发生膨胀,近邻氢原子渗入到铍表面次层,氢原子易于占据在hcp和桥位.吸附结构中的氢原子比氢分子中的原子稳定.当覆盖度大1.33ML时,计算结果没有发现相对于氢分子更稳定的吸氢结构.同时从分析偶极修正和氢原子吸附垂直高度随覆盖度的变化关系判断氢覆盖度为1.33ML时,在Be(0001)表面吸附达到饱和.     

CO在α-U(001)表面的吸附

采用密度泛函理论中的广义梯度近似,计算了CO在α-U(001)表面的吸附、解离和扩散.结果表明:CO分子以CU3OU2构型化学吸附在α-U(001)表面,吸附能为1.78-1.99eV;吸附后表层U原子向上迁移,伴随着褶皱的产生;CO分子与表面U原子的相互作用主要是U原子的电子向CO分子最低空轨道2π*转移,以及CO2π*/5σ/1π-U6d轨道间杂化而生成新的化学键;CO解离吸附较分子吸附在能量上更为有利,h1(C)+h2(O)和h1(C)+h1(O)(h:空位)解离态吸附能分别为2.71和3.08eV;近邻三重穴位之间C、O原子的扩散能垒分别为0.57和0.14eV,预示O原子较C原子更易在U(001)表面扩散迁移.     

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

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

H_2分子在Ni(110)表面吸附的位能面研究

本文采用改进的LEPS方法计算了H_2分子在Ni(110)表面的吸附位能面,结果表明H_2分子在表面吸附过程形成的吸附形态与分子接近表面的方式有关。分子垂直接近于表面时,仅有分子形态的吸附发生;而分子平行接近于表面时,可以形成分子吸附态和原子吸附态,表面吸附态的转变需要克服一个活化能为0.9～1.3eV的势垒。     

密度泛函理论研究十二烷硫醇在Au(111)面上的吸附

采用第一性原理方法研究了十二烷硫醇(C₁₂H₂₅SH)分子在Au(111)面上未解离和解离吸附的结构、能量和吸附性质,在此基础上分析判断长链硫醇分子在Au(111)面吸附时S―H键的解离, 以及分子链长度对吸附结构和能量的影响. 计算了S原子在不同位置以不同方式吸附的系列构型, 结果表明在S―H键解离前和解离后,均存在两种可能的表面结构, 直立吸附构型和平铺吸附构型; 未解离的C₁₂H₂₅SH分子倾向于吸附在top位, 吸附能为0.35-0.38 eV; H原子解离后C₁₂H₂₅S基团倾向于吸附在bri-fcc位, 吸附能量为2.01-2.09 eV. 比较分析未解离吸附和解离吸附, 发现C₁₂H₂₅SH分子未解离吸附相较于解离吸附要稳定, 未解离吸附属于弱化学吸附.局域电子态密度和差分电荷密度分析进一步验证了S―H解离后S原子与表面之间成键的数目增加, 而且键合更强. 同时我们发现长链硫醇的吸附能量较短链硫醇的吸附能量略大, S原子与表面Au原子之间的距离略小.     

NiO表面磁性对分子吸附的影响的第一性原理研究

过渡金属氧化物广泛应用在当今能源与环境相关的催化领域,理解其表面化学性质以及结构-反应活性之间的关系对于先进催化材料的进一步发展以至理性设计至关重要.3d后过渡系金属(Mn,Fe,Co,Ni)的氧化物以其中金属离子独特的自旋状态和由此产生的铁磁/反铁磁性为典型特征.研究过渡金属氧化物的自旋状态以及磁性对表面化学的影响将使我们更加完整了解这些材料的表面化学.以NiO为代表的后过渡系金属岩盐结构一元氧化物具有反铁磁性,被经常作为反铁磁研究的模型体系.尽管在低温(低于其Neel温度)下NiO体相的完整晶体具有确定的反铁磁序,但是一系列最新研究表明,在条件变化时NiO表面的Ni离子可以产生不同的磁序.以此为背景,本工作以NiO为模型体系,采用DFT+U的第一性原理方法研究了NiO表面磁序对表面的小分子吸附活性的影响,包括表面吸附活性对各磁性相的表面取向以及吸附物种磁性的依赖关系.我们考察了NiO的5种反铁磁相和一种铁磁相,两个晶面NiO(001)和NiO(011),顺磁性分子NO和非顺磁性分子CO.我们发现表面能受磁性的影响较轻微,NiO(001)面上从49到54 meV/?2,NiO(011)面上从162到172 meV/?2.在NiO(001)面上,CO与NO都倾向于在Ni离子的顶位吸附.对于不同的体相磁序与表面取向,CO吸附能的变化范围为-0.33~-0.37 eV,NO吸附能的变化范围为-0.42~-0.46 eV.在NiO(011)表面,两种分子都倾向于吸附在由两个Ni离子构成的桥位.我们发现相对于NiO不同磁性相的体相长程磁序,吸附位点处构成桥位的两个Ni离子的局部磁矩相对取向对于分子的吸附具有更加显著的影响.计算得到NO在局部磁矩相对取向反平行(↑↓)吸附位点处的吸附能为-0.99~-1.05 eV,在局部磁矩相对取向平行(↑↑)吸附位点处吸附会增强,吸附能为-1.21~-1.30 eV.对于CO,尽管计算的吸附能在(↑↓)吸附位点(-0.73~-0.75 eV)与在(↑↑)吸附位点(-0.71~-0.72 eV)非常接近,两种吸附位点处的CO吸附时分子轨道杂化方式以及吸附后CO的局域电子态密度却具有明显不同的特征.本工作突出揭示了分子在过渡金属氧化物表面的多重吸附位点上吸附时吸附位点的局域磁矩相对取向对吸附性能的影响.     

Adsorption and dissociation of H2O2 on Pt and Pt-alloy clusters and surfaces

The adsorption of H(2)O(2) on Pt and Pt-M alloys, where M is Cr, Co, or Ni, is investigated using density functional theory. Binding energies calculated with a hybrid DFT functional (B3PW91) are in the range of -0.71 to -0.88 eV for H(2)O(2) adsorbed with one of the oxygen atoms on top Pt positions of Pt(3), Pt(2)M, and PtM(2), and enhanced values in the range of -0.81 to -1.09 eV are found on top Ni and Co sites of the Pt(2)M clusters. Adsorption on top sites of Pt(10) yields a weaker binding of -0.48 eV, whereas on periodic Pt(111) and Pt(3)Co(111) surfaces, H(2)O(2) generally dissociates into two OH radicals. On the other hand, attempts to attach H(2)O(2) on bridge sites cause spontaneous dissociation of H(2)O(2) into two adsorbed OH radicals, suggesting that stable adsorptions on bridge sites are not possible for any of the clusters or extended surfaces that are being studied. We also found that the water-H(2)O(2) interaction reduces the strength of the adsorption of H(2)O(2) on these clusters and surfaces.     

Hydrogen adsorption on nickel (100) single-crystal face. A Monte Carlo study of the equilibrium and kinetics

We propose a model of the dissociative adsorption of hydrogen on nickel single-crystal face. In this model, we treat the Ni(100) surface as a strongly correlated energetically heterogeneous surface, because the density functional theory (DFT) studies indicate that hydrogen atoms may adsorb either on hollow sites (energetically more favorable, binding energy 48 kJ/mol H) or bridge sites with the binding energy less by 11 kJ/mol H. The essential assumption of the proposed model is that the dissociation of the hydrogen molecule is possible only over the topmost Ni atom, and the resulting H atoms may adsorb either on two free hollow sites (but the adjacent bridge sites must be free) or two bridge sites (the adjacent hollow sites must be free). If the above condition is not fulfilled, then the dissociation and adsorption are impossible. The second assumption is that the rate (probability) of the associative desorption is limited by the rate of diffusion of H atoms on the surface. This is because the two H atoms desorb, giving an H2 molecule, only when they meet on two adjacent hollow-bridge sites. Our model recovers very well the behavior of the experimental equilibrium adsorption isotherms as well as kinetic isotherms. As a result, we stated that hydrogen atoms are not completely free on the surface, but they cannot also be considered localized at room and elevated temperatures. Additionally, while analyzing the kinetic adsorption isotherms, we stated that the rate-limiting step during the dissociative adsorption of H2 is the disintegration of the activated complex and the subsequent adsorption of hydrogen atoms.     

Adsorption and dissociation of the methanethiol (CH3SH) molecule on the Fe(100) surface

These contributions explore interaction modes between the methanethoil (CH₃SH) molecule and the Fe(100) surface via implementing accurate density functional theory (DFT) calculations with the inclusion of van der Waals corrections. We consider three adsorption sites over the Fe(100) surface, namely, top(T), bridge (B), and hollow (H) sites as potential catalytic active sites for the molecular and dissociative adsorption of the CH₃SH molecule. The molecular adsorption structures are found to occupy either B or T sites with former sites holding higher stability by 0.17 eV. The inclusion of van der Waals corrections refound to slightly alter adsorption energies. For instance, adsorption energies increased by ~ 0.18 and ~ 0.21 eV for B and T structure, respectively, in reference to values obtained by the plain generalized gradient approximation (GGA) functional. A stability ordering of the dissociation products was found to follow the sequence (CH₄, S) > (CH₃, S, H) > (─SCH_3, H) > (─CH₃, SH). The differential charge density distributions were examined to underpin prominent electronic contributing factors. Direct fission of C─S bond in the CH₃SH molecule attains exothermic values in the range 2.0 to 2.1 eV. The most energetically favorable sites for the surface-mediated fission of the thiol's S─H bond correspond to the structure where the ─SCH₃ and H are both situated on hollow sites with an adsorption energy of −2.43 eV. Overall, we found that inclusion of van der Waals functional to change the binding energies more noticeably in case of dissociative adsorption structures. The results presented herein should be instrumental in efforts that aim to design stand-alone Fe desulfurization catalysts.     

A Density Functional Theory Study on the Effect of Zero‐Point Energy Corrections on the Methanation Profile on Fe(100)

The thermodynamics and kinetics of the surface hydrogenation of adsorbed atomic carbon to methane, following the reaction sequence C+4 H?CH+3 H?CH₂+2 H?CH₃+H?CH₄, are studied on Fe(100) by means of density functional theory. An assessment is made on whether the adsorption energies and overall energy profile are affected when zero‐point energy (ZPE) corrections are included. The C, CH and CH₂ species are most stable at the fourfold hollow site, while CH₃ prefers the twofold bridge site. Atomic hydrogen is adsorbed at both the twofold bridge and fourfold hollow sites. Methane is physisorbed on the surface and shows neither orientation nor site preference. It is easily desorbed to the gas phase once formed. The incorporation of ZPE corrections has a very slight, if any, effect on the adsorption energies and does not alter the trends with regards to the most stable adsorption sites. The successive addition of hydrogen to atomic carbon is endothermic up to the addition of the third hydrogen atom resulting in the methyl species, but exothermic in the final hydrogenation step, which leads to methane. The overall methanation reaction is endothermic when starting from atomic carbon and hydrogen on the surface. Zero‐point energy corrections are rarely provided in the literature. Since they are derived from C? H bonds with characteristic vibrations on the order of 2500–3000 cm^?1, the equivalent ZPE of 1/2 hν is on the order of 0.2–0.3 eV and its effect on adsorption energy can in principle be significant. Particularly in reactions between CH_x and H, the ZPE correction is expected to be significant, as additional C? H bonds are formed. In this instance, the methanation reaction energy of +0.77 eV increased to +1.45 eV with the inclusion of ZPE corrections, that is, less favourable. Therefore, it is crucial to include ZPE corrections when reporting reactions involving hydrogen‐containing species.     

Theoretical study of the geometries and dissociation energies of molecular water on neutral aluminum clusters Al(n) (n = 2-25)

Geometries and dissociation energies of water molecules on Al(n) (n = 2-25) clusters were investigated using density functional theory with all electron relativistic spin-polarized calculations under the generalized gradient approximation. An extensive structure search was performed to identify the low-energy conformations of Al(n)H(2)O complexes for each size. Optimal adsorption sites were assigned for low-energy isomers of the clusters. Size and site specific dependences were studied for the Al(n)H(2)O complexes in stabilities, geometries, adsorption energies, dissociation energies, Al-O bond lengths, and other characteristic quantities. The stabilities and geometries revealed that H atom in H(2)O is not inclined to bond with Al atoms. The most stable Al(n)H(2)O configurations for each size tend to correspond to the most stable bare Al(n) cluster except of Al(6) and Al(24) clusters. The HO bond lengths increase generally 0.01 ? with respect to the isolated H(2)O in all of the adsorption complexes. The dissociation energy of an isolated H(2)O into HO and H was 5.39 eV, which decreased about two-thirds to the energy range of 0.83-2.12 eV with the help of Al(n) clusters. In spite of the fluctuations, the dissociation energies of Al(n)H(2)O complexes rise with the size increasing as a whole. In addition, we also found that the bare Al(n) clusters with high vertical ionization potentials usually have high dissociation energies of H(2)O in the corresponding adsorption models. The energetically preferred spin-multiplicity of all the odd-n Al(n)H(2)O complexes is doublet, and it is singlet for all the even-n complexes with exception of Al(2)H(2)O which is triplet.     

The Surface Chemistry of Water on Fe(100): A Density Functional Theory Study

The formation of water by hydrogenation of atomic oxygen is studied using density functional theory. Atomic oxygen preferentially adsorbs at the four‐fold hollow site, the hydroxyl group prefers the bridge site in a tilted configuration, and water is most stable when adsorbed at the top site with the two O? H bonds parallel to the Fe surface. Water formation by the hydrogenation of oxygen is a highly activated process on the Fe(100) surface, with similar activation energies, in the order of 1.1 eV, for the first and second hydrogen additions. A more favourable route for the addition of the second hydrogen atom involves the disproportionation of hydroxyl groups to form water and adsorbed oxygen. Dissociation of the OH is also likely since the activation energy is similar to that for disproportionation of 0.65 eV. Furthermore, the results show that the dissociation of water on Fe(100) is a non‐activated process: 0.16 eV for the zero‐coverage limit and 0.03 eV when surface oxygen is present. Herein, adsorption energies, structures and vibrational frequencies are presented for several adsorption states at 0.25 ML coverage, as well as the potential energy surface for water formation on Fe(100).     

Pd掺杂对ZnO(1120)面上水解离的影响

采用密度泛函理论计算研究了清洁的以及Pd掺杂的ZnO(1120)面上水分子的吸附和解离.结果表明,在清洁ZnO(1120)上,水分子倾向于分子吸附,解离吸附较为困难.在Pd掺杂的ZnO上,水分子仍倾向吸附在Zn原子上,且吸附能与其在清洁ZnO表面的相当.然而,Pd的掺杂可增强水解离产物OH和H的吸附,从而显著提高了水的解离活性,相应的水解离能垒为0.36eV,放热0.21eV.     

A DFT study of H2O dissociation on metal‐precovered Fe (100) surface

The H₂O adsorption and dissociation on the Fe (100) surface with different precovered metals are studied by density functional theory. On both kinds of metal‐precovered surface, H₂O molecules prefer adsorb on hollow sites than bridge and top sites. The impurity energy difference is proportional to the adsorption energy, but the adsorbates are not sensitive to the adsorption orientation and height relative to the surface. The Hirshfeld charge analysis shows that water molecules act as an electron donor while the surface Fe atoms act as an electron acceptor. The rotation and dissociation of H₂O molecule occur on the Co‐ and Mn‐precovered surfaces. Some H₂O molecules are dissociated into OH and H groups. The energy barriers are about 0.5 to 1.0 eV, whose are consistence with the experimental data. H₂O molecules can be dissociated more easily at the top site on Co‐precovered surface 1 than that at bridge site on Mn‐precovered surface 2 because of the lower reaction barrier. The dispersion correction effects on the energies and adsorption configurations on Co‐precovered surface 1 were calculated by OBS + PW91. The dispersion contributions can improve a bit of the bond energy of adsorbates and weaken the hydrogen bond effect between adsorption molecules a little.