H_%282%29S在Cr%28111%29面上吸附与解离的第一性原理研究.pdf

第41卷第2 期2024年4月J.At.Mol.Phys.,2024,41:021002(6pp)H,S在Cr(111)面上吸附与解离的第一性原理研究原子与分子物理学报JOURNAL OF ATOMIC AND MOLECULAR PHYSICSVol.41 No.2Apr.2024江佐禹，黄本生，何子涛，罗霞1（1.西南石油大学新能源与材料学院，成都6 10 50 0；2.西南石油大学四川省玄武岩纤维复合材料及开发应用工程技术研究中心，成都6 10 50 0）摘要：采用了第一性原理研究了H,S在Cr(111)面的吸附解离过程,利用吸附能、吸附构型和偏态密度图(PDOS)研究了H,S及其解离产物在Cr(111)面上的吸附情况,都偏向倾斜吸附在Cr(111)面同时研究了HS/H和 S/H共吸附情况,得到共吸附物质在Cr(111)面上有明显的相互作用。最后使用线性同步和二次同步变换方法确定了解离反应的过渡态，了解到第一、二步解离的活化能分别为1.6 5eV、0.8 2 e V，H,S分子在Cr(111)面上的解离过程是放热反应,反应能为-2.90 eV.关键词：H,S；Cr(111)；吸附；解离；第一性原理中图分类号：0 6 47.31文献标识码：AFirst-principles study on adsorption anddissociation of H,S on Cr(111)surfaceDOI:10.19855/j.1000-0364.2024.021002JIANG Zuo-Yu,HUANG Ben-Sheng2,HE Zi-Tao,LUO Xia(1.School of New Energy and Materials,Southwest Petroleum University,Chengdu 610500,China;2.Sichuan Engineering TechnologyResearch Center of Basalt Fiber Composites Development and Application,Southwest Petroleum University,Chengdu 610500,China)Abstract:First principles were used to study the adsorptive dissociation process of H,S on the Cr(111)plane,and adsorptions of H,S molecules and their dissociated products on the Cr(111)plane were studied using ad-sorption energy,adsorption configuration,and Partial density of states(PDOS)plots,all of which are biased to-ward inclined adsorption on the Cr(111)plane.At the same time,the co-adsorption of HS/H and S/H wasstudied,and it was found that the co-adsorbed species had obvious interactions on the Cr(111)surface.Final-ly,the transition state of the dissociation reaction was determined using linear synchronization and quadratic syn-chronization transformation methods.It is known that the activation energies of the first and second steps dissocia-tion are 1.65 eV and 0.82 eV respectively,and the dissociation process of H,S molecule on the Cr(111)surfaceis an exothermic reaction,and the reaction energy is-2.90 eV.Key words:H,S;Cr(111);Adsorption;Dissociation;First-principles体广泛存在于石油和天然气开采过程中，在干燥1 引 言环境中不具有腐蚀性，但是一旦溶于水会使水溶随着科技的发展，人们对能源的需要越来越液产生酸性，这种物质不仅会引起管道的氢致开大。为了满足需求，越来越多的含硫油气田被开裂和硫化物应力腐蚀开裂，还会导致大多数的金发出来，同时也带来了原油、天然气开采过程中属催化剂失活2-5因此研究H,S在金属表面的的管道等器材腐蚀问题 H,S作为一种伴生气反应过程对防腐材料的开发具有重要意义：收稿日期：2 0 2 2-0 8-0 3基金项目：四川省科技厅国际合作基金（2 0 2 0 YFH0151）作者简介：江佐禹（1998 一），男，江西南城人，硕士研究生，主要从事气体分子在金属表面吸附的研究。E-mail：10 7 52 7 8 8 6 4 q q.c o m通讯作者：黄本生.E-mail：h b s l x p 16 3.c o m021002-1第41卷为了提高材料的耐蚀性，人们往其中添加Cr、Cu 和Mo等合金元素，Yang等人研究Cu的添加对 Ni-Co-CrM o 合金耐蚀性的影响，其中Cu含量为2%时，此时形成的钝化膜最厚，耐蚀性能最好6，Ajit研究发现，在镍基合金中加人Cr、M o 元素，在盐酸和硫酸溶液中表现出很强的均匀耐腐蚀性，同时其合金的局部耐腐蚀性高于其他的在售合金7。Liu等人通过对比分析添加不同 Cr 含量的 Ni-Cr-Mo涂层在 H,S 下的腐蚀性能，发现随着Cr含量的增加其耐腐蚀性能显著增加，当Cr含量为15wt.%时，其耐腐蚀性能和用量是最优的8 从这些研究可以看出，Cr对材料的耐腐蚀性的提升具有很大的作用然而令人惊的是关于Cr金属表面微观机理研究却非常稀缺，现仅发表了一篇在Cr表面上的研究，Anil等人运用了DFT模拟研究了氢同位素在体心立方Cr和Cr(100)表面上的化学吸附过程9如今，有很多文献研究了H,S 在金属10-19及其合金12 0-2 表面上的吸附解离过程H,S在这些表面上的吸附过程通常被描述成一个跨越低能垒的简单过程，而H,S在Cr表面上的具体吸附机理并无研究本文选用 Cr(111)面研究 H,S 吸附解离过程，为后续的含铬防腐材料的开发提供更深入的了解。2计算方法本文采用Material Studio软件中的CASTEP模块2 3进行计算，所有计算过程都是利用自旋无限制的平面波密度泛函进行采用的是广义梯度近似(GGA）下的 Revised-Perdew-Burke-Ernzerh(RPBE)泛函2 4进行研究体计算和表面计算的布里渊区积分分别使用(444)和（441)的K点网格，截断能设置为40 0 eV.为了保证收敛，能量收敛准则设置为1.0 10-eV/atom，原子作用力收敛准则设置为0.0 3eV/A，应力偏差小于0.05GPa，公差偏移小于0.0 0 1A.在基于优化的Cr晶体结构，建立了4层的Cr(111)表面模型，采用了p(22)的超晶胞结构，吸附物的覆盖率为1/4ML.其中最下面两层原子被固定住，其余层的原子完全弛豫，为了保证周期性平板之间不会发生相互作用，在平板的z方向加上了15A的真空层为了保证计算的可靠性，我们计算了Cr的晶格常数（0.2 8 6 nm）与实验值(0.2 8 8 nm)非常接近，同时 H,S、H S 和 H的键长和键角理论值与实验值差别如表1所示，原子与分子物理学报误差值均小于2.2%.通过线性同步变换（LST）和二次同步变换(QST)方法得到反应的最小能量路径，研究H,S在Cr(111)表面上的吸附解离过程，并且确定过渡态(TS).吸附能公式表示如下：其中Eslal EadobalesEasorbaes/sla分别表示洁净的表面，在气相中的吸附物分子，和具有吸附物分子的表面的能量由定义可知，吸附能为负值则为放热反应数值的绝对值越大，放热能量越多，说明吸附作用越强：对于在Cr(111)面上计算的活化能E和反应能E可以表示如下：E。=E(T S)-E(RE)E=E(PR)-E(RE)其中E(RE)代表反应物在Cr(111)面上初始态的总能量，E(TS)代表过渡态在Cr(111)面上的总能量，20E(PR)代表最终产物在Cr(111)面上的总能量.3结果与讨论我们考虑了几种不同的吸附位置研究H,S及其解离产物在Cr(111)面上的吸附：代表Cr原子顶部的位点Top，代表两个Cr原子之间的桥接位点Bridge.如图1所示，Cr(111)面有3层原子裸露在外面，在Cr(111)面上顶部的位点就有三种，分别为：第一层Cr原子的顶部（T）、第二层Cr原子的顶部（T）和第三层Cr原子的顶部（T,）；同时桥接位点也有三种，分别是：第一层Cr原子与第二层Cr原子之间的桥接位点（B12）、第二层Cr原子与第三层Cr原子之间的桥接位点（B23）和第一层Cr原子与第三层Cr原子之间的桥接位点（Bi 3）本次研究黄色圆球代表S，白色圆球代表H，其他颜色原子均为Cr，仅考虑了垂直吸附.表1气相H,S、H S、H,的几何参数Table 1 The geometrical parameters of gas-phase H,S,HS,H2.H,SCal.Exp.25r(A)1.3561.3280(deg)92.03.1H,S及其解离产物的单一吸附情况对于 H,S 及其解离产物在 Cr(111)表面单一021002-2第2 期(1)(2)HSH2Cal.Exp.251.3601.3460.75491.6一Cal.Exp.250.740一一第41卷(a)江佐禹，等：H,S在Cr(111)面上吸附与解离的第一性原理研究(b)第2 期和2.499 A.而H在 Cr表面的最优吸附位为T,位点（图2（d)），仅与一个Cr原子成键，吸附能和键长分别为-0.45eV和1.7 6 1A.(a)(b)图1（2 2）Cr(111）表面示意图：（a）俯视图，（b）侧视图T代表顶部站点，B代表桥接站点侧视图中的1、2 和3表示Cr(111)表面上的第一层、第二层和第三层原子。Fig.1(2 2)Cr(111)surface schematic diagram:(a)top view,(b)side view.The T stands forthe top site and the B for the bridge site.1,2and 3 on the side view represent the first,sec-ond and third layer atoms on the Cr(111)sur-face.吸附的结构图如图2 所示，相应最优吸附位置、吸附能以及结构信息如表2 所示。研究发现，H,S的稳定吸附位点是T,位点（图2（a），而吸附在其他位点是都发生了不同程度的分解和偏移当H,S 吸附在T,位点，由一开始的垂直吸附变成了倾斜吸附，这意味着顶部的倾斜吸附更加稳定同时H-S键和ZH-S-H发生了不同程度的增大，表明 Cr表面将H,S活化，导致H-S键被激活且强度减弱.对于HS，由于其的不稳定性，并没有在实验上观察到其吸附在金属表面上的吸附结构，根据吸附能来看，在T位点的吸附能最小为-2.81eV,在B23位点的吸附能最大为-3.8 4eV.并且HS 吸附在T、T,和Bi2位点时，都往 Cr(111)表面上的B23位点发生了偏移而当HS吸附在Bi3时，往Cr(111)表面上的T,位点发生偏移.这进一步说明了 B23位点是HS 在Cr(111)表面上最稳定的吸附位（图2(b)）根据H-S键吸附前后对比，结果表明，H-S键更倾向于倾斜吸附于基体表面，H-S键长由大气中1.353.A延长到了1.37 3A，说明Cr表面活化了H-S键导致H-S键增长，H-S键减弱且容易断裂.对于S 和H，S在Cr(111)面上的6 种吸附位置都较为稳定，均呈垂直吸附于基体表面上，但当S处于Cr(111)面Biz位点和Bi3位点时，经过结构优化后分别转移到T,位点和T,位点；其中在T位点的吸附能是最大的（图2（c），为-6.76 eV.吸附时与周围Cr原子形成了4个S-Cr键，键长分别为2.343A，2.49 6 A，2.49 6 A(c)图2 1H,S及其解离产物在Cr(111)面上的最优吸附位:(a)H,S,(b)HS，(c)S,(d)HFig.2 The optimal adsorption sites of H,S and its dis-sociation products on the Cr（1l 11)s u r f a c e:(a)H,S,(b)HS,(c)S,(d)H表2 H,S及其解离产物在Cr(111)面上的最优吸附结构参数以及吸附能Table 2 The optimal adsorption structure parameters and ad-sorption energies of H,S and its dissociation productson Cr(111)surfaceOptimalAdsorptionadsorptionenergy/evpositionH,STHSB23ST2HT2为了进一步理解H,S及其解离产物在Cr(111)面上的单一吸附情况，我们绘制了H,S、HS、S和H在Cr(111)面吸附前后的偏态密度图(PDOS)如图3所示从态密度图中可以看出，吸附前后S元素的s轨道波峰、p轨道波峰和H元素的s轨道波峰均显著往低能级方向移动，说明H,S及其解离产物在Cr(111)表面上均为化学吸附.根据图3(a）、图3(b)和图3（c)中吸附后的态密度图可以看出部分S-p轨道电子与Cr-d轨021002-3(d)dH-s/Adm-Cr(m=s.m)/A0.171.3633.841.373-6.76-0.452.5522.483/2.4852.343/2.496/2.496/2.4991.761第41卷道电子处在同一能级表明H,S、H S和S吸附时S-Cr键的形成主要是S-p轨道电子和Cr-d轨道电子的杂化共轭作用由图3(d)吸附后的态密度图可以看出少量的H-s轨道电子与Cr-d轨道处于同一能级，H吸附时H-Cr键的形成主要是H-s轨道电子和Cr-d轨道电子的杂化共轭作用.50Before adsorption40-30-20-1050-40-30-20-10-0-1650-Before adsorption40-30-2010050-Afteradsorption40-30-20-10-0-16图3H,S及其解离产物吸附前后的偏态密度图：（a）H,S，（b）H S，（c）S，（d）HFig.3 Partial density of states(PDOS)plots of H,S and its dissociation products before and after adsorp-tion:(a)H,S,(b)HS,(c)S,(d)H3.2HS/H 和 S/H 在 Cr(111)面上的共吸附情况图4(a)显示了 HS 和H在 Cr 表面的共吸附构型，HS处于B2s位置，H处于T,位置与HS/H单一吸附的情况不同的是，在共吸附情况下，S-Ni键长和HNi 键长相对于单一吸附情况下均有所增加，而H-S键长从1.37 3A降到了1.37 1A略有降低其共吸附的吸附能为4.2 8 4eV，小于单一吸附情况下HS与H吸附的吸附能总和-4.309 eV.说明H与HS之间在 Cr(111)面上吸附时具有一定的相互作用。图4(b)显示了S和H在Cr表面的共吸附构型.S处于T,位置，H处于S相邻的T,位置，不同于单一吸附情况的是，S只与周围3个Cr原子成键，S-Cr键长（2.2 50 A-2.552A)不同于单原子与分子物理学报通过对比 S、H S和H,S吸附时态密度图中S-p轨道与Cr-d轨道的重叠面积可以发现，说明S吸附体系上S-p轨道与Cr-d轨道的杂化共轭作用是最强的，其次是HS分子，最后是H,S分子，进一步证明了在Cr(111)表面上的吸附强度大小排序:SHS H,S.50S-sS-PCr-dAfteradsorptionS-sS-pCr-d-12-8-12-8第2 期Before adsorptionS-s40S-pCr-d.30-2010050-40-30-20-10-40Energy(ev)(a)40Energy(ev)(c)Afleradsorplion488S-sS-pCr-d1216S-sS-pCr-dS-sS-PpCr-d1216021002-4-1650-Beforeadsorption40-30-201040-30-20-10-16(a)图4Cr(111)面上的共吸附构型：（a）H S/H，（b)S/HFig.4 Co-adsorption configurations on Cr(111)sur-face:(a)HS/H,(b)S/H一吸附情况的S-Cr键长(2.343A-2.499A).H在T,位置也是倾斜吸附状态，且 H-Cr 键的键长-12Afteradsorption-12-8-8-4Energy(ev)(b)-40Energy(ev)(d)(b)044812812H-sCr-dH-sCr-d1616第41卷从1.7 6 1A增加到了1.7 6 7 A.在Cr(111)表面,S与H的吸附能为之和为-7.2 0 5 eV，与共吸附能（-6.931 eV)不同这种现象说明S与H之间在Cr(111)表面上也具有一定的相互作用.4H,S 的解离步骤我们在确定了H,S及其解离产物在Cr表面上的稳定吸附位置之后，通过LST/QST方法确认了H,S在Cr(111)面上解离过程的最小能量路径如图5所示，同时绘制了H,S在Cr(111)面上解离的能量分布图（图6)H,S分子的解离步骤有两步，分别为：H,SHS+H，H S-S+H 两个步骤由图4可知，H,S首先倾斜吸附于第一层Cr原子的顶位（T,），S与一个Cr原子成键，随后发生了第一步解离：第一个H从H,S上脱离，向二层Cr原子的顶部位点（T,）方向移动，此时H还没有与Cr原子成键同时HS也向B23位置偏移。在TS1之后，H与一个Cr 原子成键，呈倾斜吸附.而HS也完全倾斜吸附在B23位点上，其中S与周围三个Cr原子成键，在第二步解离过程开始时，HS的最优稳定位置是处于第二、三层Cr原子之间的桥接位点（B23），随着解离的发生，H往另外一个二层Cr原子的顶位方向（T,）偏移，形成TS2，这时H-S键还没发生断裂随着反应的继续进行，H-S键发生断裂，H移动到了二层Cr原子的顶部位点，与一个Cr原子成键，呈倾斜吸附；S移动到了另一个二层Cr原子的顶部位点，与周围的3个Cr原子成键在第一步解离过程中，反应的活化能为1.6 5eV，反应为放热反应（-1.9 5 eV），说明第一步解离反应在热力学与动力学层面上是可以自发进行的第二步解离反应的活化能较第一步反应的活化能偏低，为0.82eV，反应能为-0.9 5eV，说明第二步解离反应比第一步反应更加容易进行.5结论(1)对于单一吸附在 Cr(111)面的情况，H,S偏稳定吸附在第一层Cr原子的顶部位置，HS偏稳定吸附在B23位置，而S与H都偏稳定吸附在T,位置，这几种物质在Cr(111）面上的吸附能大小顺序为：S HS HH,S.(2)H,S、H S、S 和H能稳定吸附在 Cr(111)面主要是因为 S-p轨道、H-s轨道和Cr=d轨道发生杂化共轭作用：同时通过S-p轨道与Cr江佐禹，等：H,S在Cr(111)面上吸附与解离的第一性原理研究TS1HS(ad)TS2图5H,S的解离路径Fig.5The dissociation pathway of H,STS11.51.65eVH,Sea)0.0H,S(ats)0.17eV1.95cV-1.5-3.0-图6 H,S 在Cr(111)面上解离的能量分布图Fig.6 Energy distribution of dissociation of H,S onCr(111)surface-d轨道的重叠面积进一步证明S的吸附能最大，HS 其次,H,S 最后.(3)在共吸附情况下，HS的吸附情况变化不大，两个共吸附构型中H都变成了倾斜吸附S吸附时与单一吸附情况下不同，倾斜吸附在T,位置，仅与周围3个Cr原子成键两种共吸附情况的吸附能均与单一吸附情况有所不同，说明共吸附时吸附物质之间存在相互作用。(4)H,S 在Cr(111)面上的解离过程是放热的,反应能为-2.90 eV.同时第一步解离的活化能为1.65eV，约为第二步解离活化能(0.8 2 eV)的两倍.说明H,S若想在Cr表面上发生解离过程，首先需要跨越一个相对较大的能垒，才能进一步解离.参考文献：1Wu M,Xie F,Chen X,et al.Reserch progress and021002-5第2 期H,S(ads)HS+HH+STS20.82eVHS+HHS0.95VReactioncoordinateS+H第41卷thinking on corrsion failure of buried oil and gas pipe-linesJ.Oil Gas Sto.Trans.,2022,41:712(inChinese）吴明,谢飞,陈旭,等埋地油气管道腐蚀失效研究进展及思考J油气储运，2 0 2 2，41：7122Wang H B,Zhang C,Lu C S,et al.Understanding ofthe corrosion behavior of carbon steel in the NH,Cl-H,S-CO,coexisting environmentJ.Int.J.Electro-chem.Sci.,2022,17:22035.3Sui P F,Sun J B,Hua Y,et al.Effect of temperatureand pressure on corrosion behavior of X65 carbon steelin water-saturated CO,transport environments mixedwith H,SJ.Int.J.Greenhouse Gas Control.,2022,73:60.4Yang S Z,Li C F,Li H,et al.Stress corrosion behav-ior and cracking mechanism of 2205 duplex stainlesssteel in acidic H,S environmentJ.Rare Met.Mater.Eng.,2018,47;904.5 Dou X M,Veksha A,Chan W P,et al.Poisoningeffects of H,S and HCl on the naphthalene steam refor-ming and water-gas shift activities of Ni and Fe cata-lystsJ.Fuel,2019,241:1008.6Yang B B,Shi C Y,Li Y P,et al.Effect of Cu on thecorrosion resistance and electrochemical response of aNi-Co-Cr-Mo alloy in acidic chloride solution J.J.Mater.Res.,2018,33:3801.7Ajit M.Performance of corrosion-resistant alloys inconcentrated acidsJ.Acta Metall.Sin.,2017,30:306.8Liu CC,Liu Z D,Gao Y,et al.Effect of Cr contenton corrosion resistance of Ni-xCr-Mo laser-clad-ding coatings under H,S-induced high-temperaturecorrosion atmosphereJ.Materials,2022,15:1885.9Anil B,Sakshi B,Sk Musharaf A,et al.Chemisorp-tion,diffusion and permeation of hydrogen isotopes inbcc bulk Cr and Cr(100)surface:First-principlesdft simulationsJ.J.Nucl.Mater.,2021,543:152538.10 Chen S H,Sun S Q,Lian B J,et al.The adsorptionand dissociation of H,S on Cu(100)surface:A DTFstudyJ.Surf.Sci.,2014,620:51.11 Li F,Zhou Z H,He C Z,et al.H,S dissociation ondefective or strained Fe(110)and subsequent forma-tion of iron sulfides:a density functional theory studyJ.Surf.Sci.,2021,709:121835.12 Jiang Z,Qin P,Fang T.Investigation on adsorption原子与分子物理学报and decomposition of H,S on Pd(100)surface:A DFTstudyJ.Surf.Sci.,2015,632:195.13 Zhang T T,Tang Q L,Yao M Y,et al.Quantumchemical DFT study of molecular adsorption of H,S onclean and chemically modified Au(110)surfacesJ.Appl.Surf.Sci.,2021,542:148595.14Dominic R A,Anthony V C,Dan C S.Adsorption anddecomposition of H,S on Pd(111)surface:a first-principles studyJ.Catal.Today,2005,99:315.15Zhang M H,Fu Z Z,Yu Y Z.Adsorption and decom-position of H,S on the Ni(111)and Ni(211)surfaces:A first-principles density functional study JJ.Appl.Surf.Sci.,2019,473:657.16Peng S F,Ho J J.Theoretical study of H,S dissocia-tion and sulfur oxidation on a W(111)surfaceJ.J.Phys.Chem.C,2010,114:19489.17Dominic R A.First-principles studies of H,S adsorp-tion and dissociation on metal surfacesJ.Surf.Sci.,2008,602:2758.18 Jiang D E,Emily A C.First principles study of H,Sadsorption and dissociation on Fe（110）J.Su f.Sci.,2005,583:60.19Jiang D E,Emily A C.Adsorption,diffusion,anddissociation of H,S on Fe(100)from first principlesJ.J.Phys.Chem.B,2004,108:19140.20 Liu N,Wang X Y,Wan Y L.First principle caleula-tions of hydrogen sulfide adsorption and dissociation onpure Pd（111)a n d A u（111),a n d a l l o y Pd/A u(111)and Au/Pd(111)surfaces J.J.Theor.Comput.Chem.,2014,13;14.21Ekin O,Jennifer W.Investigation of H,and H,S ad-sorption on niobium-and copper-doped palladiumsurfacesJ.J.Phys.Chem.B,2010,114:12851.22Erik J A,Abolghasem S.Density functional theorystudy of hydrogen sulfide dissociation on bi-metallicNi-Mo catalystsJ.Surf.Sci.,2006,600:3202.23 Clark S J,Segall M D,Pickard C J,et al.First princi-ples methods using CASTEP J.Z.Kristllgr.,2005,220:567.24Hammer B.Improved adsorption energetics within den-sity-functional theory using revised Perdew-Burke-Ernzerhof functionals J.Phys.Rev.B,1999,59:7413.25Villars P,Calvert L D.Pearsons handbook of crystal-lographic data for intermetallic phase M .2th ed.Mater Park,OH:ASM Int.,1991.第2 期021002-6