腐蚀介质在缓蚀剂膜中扩散行为的分子动力学模拟

1、NovemberArticleActaPhys.鄄Chim.Sin.,2010,26(11)：3041-3046物理化学学报(WuliHuaxueXuebao)腐蚀介质在缓蚀剂膜中扩散行为的分子动力学模拟乔贵民任振甲张军胡松青燕友果257061)提阳(中国石油大学物理科学与技术学院,山东东营摘要：采用分子动力学模拟方法,从缓蚀剂膜阻碍腐蚀介质粒子(H2O、H3O+和HCO3-)向金属表面扩散的角度,研究了4种1R12十一烷基咪唑啉缓蚀剂(R1:羧甲基(A),羟乙基(B),氨乙基(C),氢(D)抑制碳钢CO2腐蚀的粒子与膜的相互作用缓蚀机理,并对其缓蚀性能进行了理论评价.腐蚀介质粒子在不同缓蚀剂

2、膜中的扩散系数、能以及膜的自扩散性能的计算结果表明:4种缓蚀剂均可形成稳定的缓蚀剂膜,能有效阻碍腐蚀介质粒子向金属表面的扩散,达到抑制或延缓腐蚀的目的;随亲水支链(R1)极性的增加,缓蚀剂膜对腐蚀介质粒子扩散行为的-抑制能力逐渐增强;同种缓蚀剂膜对正负离子H3O+和HCO3比对中性的H2O分子具有更强的扩散抑制能力.综合计算及分析结果,4种缓蚀剂缓蚀性能的理论评价结果为ABCD,与文献实验结果吻合.关键词：分子动力学模拟;中图分类号：O641缓蚀剂膜;腐蚀介质;扩散系数MolecularDynamicsSimulationofCorrosiveMediumDiffusioninCorrosio

3、nInhibitorMembrane(CollegeofPhysicsScienceandTechnology,ChinaUniversityofPetroleum,DongyingQIAOGuiMinRENZhenJiaZHANGJunHUSongQing257061,ShandongProvince,P.R.China)YANYouGuoTIYangAbstract：Thecorrosioninhibitionmechanismoffour1R12undecylimidazolineinhibitors(R1:CH2COOH(A),CH2CH2OH(B),CH2CH2NH2(C),H(D)

4、forcarbonsteelagainstcarbondioxidecorrosionwasinvestigatedbymolecular-dynamicssimulation,fromtheaspectofcorrosivemediumparticle(H2O,H3O+,andHCO3)diffusiontothemetalsurfacehinderedbythecorrosioninhibitormembrane.Thecorrosioninhibitionperformanceofthesecompoundswasalsoevaluatedbythetheoreticalmethod.T

5、hediffusioncoefficientsinvariouscorrosioninhibitormembranes,theinteractionenergiesbetweenparticlesandmembranes,andtheselfdiffusionperformanceofthemembranesindicatedthatthefourimidazolineinhibitorscouldformstablemembranes,whichcouldeffectivelylimitedthediffusionofcorrosivemediumparticlestothemetalsur

6、face,inordertoinhibitordelaycorrosion.Withanincreaseinthepolarityofthehydrophilicchain(R1),theabilityofthemembranetohinderparticlediffusionenhanced.Themembranewas-betteratlimitingthediffusionofchargedparticles(H3O+andHCO3)thanthatofaneutralparticle(H2O).Basedontheaboveanalysis,wefoundthattheoretical

7、lythecorrosioninhibitionperformanceofthefourimidazolineinhibitorsdecreasedasfollows:ABCD,whichagreeswithpreviousexperimentalresults.KeyWords：Moleculardynamicssimulation;DiffusioncoefficientCorrosioninhibitormembrane;Corrosivemedium;Received:May6,2010;Revised:July26,2010;PublishedonWeb:September9,201

8、0.Correspondingauthor.Email:dynamic_zh;Tel:+865468391592.TheprojectwassupportedbytheChinaNationalPetroleumCorporationInnovationFoundation(07E1021,2008D500602)andNaturalScienceFoundationofShandongProvince,China(Y2006B35).中国石油中青年创新基金(07E1021,2008D500602)及山东省自然科学基金(Y2006B35)资助项目EditorialofficeofActaPhy

9、sicoChimicaSinica3042ActaPhys.鄄Chim.Sin.,2010Vol.26腐蚀是工业生产过程中一个极为突出的问题,造成了巨大的经济损失和安全隐患.在众多的防腐蚀措施之中,添加缓蚀剂是广泛采用的方法之一.常用的缓蚀剂主要为吸附型缓蚀剂,它可稳定吸附在金属表面并在金属表面形成致密的保护膜,阻碍腐蚀介质向金属表面迁移扩散,以达到抑制或延缓腐电蚀的目的1.目前研究人员2-6主要通过腐蚀电流、化学阻抗及缓蚀效率等宏观参量的变化来评价缓蚀剂膜性能的优劣,由于实验技术本身的限制,还不能揭示缓蚀剂膜阻碍腐蚀介质扩散的微观机理,缓蚀机理还有待于进一步完善.目前,分子动力学模拟(MD)

10、方法已发展成为一种能从分子水平对复杂体系进行研究的有效手段,可跟踪复杂体系随时间的动态演化过程7-9,描述微小时间尺度内粒子的动态扩散行为.国内外学者10-16曾采用MD方法对聚合物膜、反渗透膜等体系中小分子的扩散动力学过程进行研究,为深入分析小分子扩散行为的微观机制提供了依据.但针对腐蚀介质在缓蚀剂膜中扩散行为的MD研究却鲜有报道.鉴于此,本文以4种含有不同亲水支链的1R12十一烷基咪唑啉缓蚀剂为研究对象(见表117),采用分子动力学模拟方法,研究CO2腐蚀环境中存在的3种腐蚀介质粒子在不同缓蚀剂膜中的扩散行为;粒子与缓蚀剂通过计算腐蚀介质粒子的扩散系数、膜的相互作用能和缓蚀剂膜的自扩散系数

11、等来分析缓蚀剂膜抑制腐蚀介质粒子扩散的微观机理,并在此基础上对4种咪唑啉缓蚀剂抑制CO2腐蚀的缓蚀性能进行理论评价.的分子模型,通过discover模块下的minimizer工具采用smartminimizer方法对其结构进行优化,然后利用amorphouscell模块分别构建包含80种缓蚀剂分子的无定形组织结构,采用等温等压系综(NPT)进行分子动力学模拟,平衡后统计体系密度的平均值,并把该值作为计算体系中缓蚀剂膜的密度.模拟体系由三层结构组成:第一层为Fe(001)面,其尺寸为3.153nm×3.153nm×1.576nm;第二层为包含80个缓蚀剂分子和1个腐蚀介质粒子

12、的无定形组织结构;第三层是厚度为2nm的真空层;然后“buildlayer”通过visualizer模块下的命令将三层结构组合在一起即可获取本文模拟所需的计算模型,如图1所示.在室温条件下,金属表面原子的热振动可以忽表面体系中的所有Fe略,因此在计算过程中“冻结”原子20.腐蚀介质粒子在缓蚀剂膜中的扩散行为的模拟通过discover模块的正则系综(NVT)21来实现,模拟温度为298K,温度和压力分别采用Andersen22和Berendsen23方法控制,各分子起始速度由MaxwellBoltzmann分布随机产生,在周期性边界条件和时间平均等效于系综平均等假设基础上,运用velocityv

13、erlet算法24求解牛顿运动方程.范德华和库仑相互作用采用chargegroup方法25计算.截断半径选取1.5nm(splinewidth:0.10nm,bufferwidth:0.05nm),截断距离之外的分子间相互作用按平均密度近似方法进行校正.模拟时间为2000ps,步长为1.0fs,每1000步记录一次体系的轨迹信息(帧),共输出2000帧.通过温度和能量判据来判断体系是否已达到平衡,图2为H2O分子在D缓蚀剂膜中扩散时体系的能量和温度随时间演化曲线,可以看出500ps后体系温度和能量波动平缓,能量偏差仅为2.35%,这表明此后体系已达到平衡,分析其它模拟体系也可得出相同的结论.因

14、此各体系均取500-1900ps内的1计算方法本文所有模拟工作均由Accelrys公司开发的MaterialsStudio软件包18完成,力场为COMPASS力场19,它能够准确给出孤立体系和凝聚态的结构与性质.首先利用visualizer模块构建A、B、C和D表1四种咪唑啉类缓蚀剂的分子结构和缓蚀效率实验值Table1Conformationsandexperimentalinhibitionefficiencies(IE)forfourimidazolineinhibitorsaCH2COOHCH2CH2OHHa1carboxymethyl2undecylimidazoline1(2hyd

15、roxyethyl)2undecyl-imidazoline1H2undecylimidazolineABCDb93.7%91.4%88.3%83.2%CH2CH2NH21(2aminoethyl)2undecylimidazoline2103bTheinhibitionefficiency(IE)istakenfromRef.17,whichwasinvestigatedthroughlinearpolarisationresistancemethod.No.11乔贵民等：腐蚀介质在缓蚀剂膜中扩散行为的分子动力学模拟3043图13种腐蚀介质粒子在缓蚀剂膜D中扩散的模型Fig.1Modelso

16、fthreecorrosiveparticlesdiffusionincorrosioninhibitormembraneD(a)H2O;(b)H3O+;(c)HCO3-(a)(b)图2H2O在缓蚀剂膜D中扩散时温度(a)和能量(b)随时间的演化曲线Fig.2Temperature(a)andtotalenergy(b)fluctuationcurvesofH2OdiffusionincorrosioninhibitormembraneD数据(最后100ps数据存在统计误差故舍去)计算腐粒子与膜的相互作用能及蚀介质粒子的扩散系数、缓蚀剂分子的自扩散系数等参量.2.1腐蚀介质粒子在缓蚀剂膜中的扩

17、散系数扩散系数可定量描述腐蚀介质粒子在缓蚀剂膜中扩散能力的强弱.腐蚀介质粒子的扩散系数越小,说明粒子的扩散迁移能力越弱,缓蚀剂膜的缓蚀性能越好;反之,缓蚀性能则越差.腐蚀介质粒子的扩散系数(D)可通过爱因斯坦关系式26求得,公式如下:D=limRi(t)-Ri(t0)2(1)ti其中,Ri(t0)和Ri(t)分别为第i个分子或粒子在初始时刻和t时刻的位置,n为分子或粒子的总数.-腐蚀介质粒子(H2O、H3O+和HCO3)在纯水体系2结果与讨论界面型缓蚀剂缓蚀性能的优劣,关键在于缓蚀剂膜能否有效抑制腐蚀介质向金属表面扩散.腐蚀介质在缓蚀剂膜中的扩散行为是一个极其复杂的过程,它不仅取决于粒子自身的

18、结构、电荷分布和迁移能力等,同时也与腐蚀介质和膜之间的相互作用以及缓蚀剂膜的自扩散性能等密切相关.n3044ActaPhys.鄄Chim.Sin.,2010Vol.26表2腐蚀介质粒子在纯水体系及缓蚀剂膜中的扩散系数Table2DiffusioncoefficientsofcorrosiveparticlesinaqueousphaseandcorrosioninhibitormembranesMembrane92·s-123+-3A0.00630.00340.0024B0.02530.00430.0036C0.03070.00490.0040D0.06180.00610.0058和

19、缓蚀剂膜中的扩散系数列于表2.由表2可见,纯水体系中H2O分子自扩散系数的计算值为2.3658×10-9m2·s-1,与实验室测得扩散系数的数值2.30×10-9-2.50×10-9m2·s-1相吻合27-29,说明本文所采用的计算方法和模型是合理和可靠的.与纯水体系相比,3种腐蚀介质粒子在缓蚀剂膜中的扩散系数均大幅度下降,说明4种咪唑啉缓蚀剂所形成的缓蚀剂膜均能有效抑制腐蚀介质向金属表面扩散,都可以起到缓蚀的效果.对于同一缓蚀剂膜,分析3种腐蚀介质粒子在其中的扩散系数可以发现,H3O+和HCO-3的扩散系数均比H2O的扩散系数小,这说明缓蚀剂

20、膜对带电粒子扩散的抑制能力更强,可有效提高腐蚀原电池的电阻,减缓腐蚀反应速率.从表2中还可以看出,3种腐蚀介质在缓蚀剂膜A-D中的扩散系数逐渐变大,说明缓蚀剂膜抑制粒子扩散的能力逐渐降低,4种咪唑啉缓蚀剂的缓蚀效率应满足:ABCD,这与实验室测得缓蚀效率的变化规律相吻合17.其原因可能与亲水支链R1中官能团COOH,OH,NH2和H的极性有关.羧基(COOH)是一个强极性基团,极性比羟基(OH)和氨基(NH2)都强,羧基中的羟基和羰基都可以形成氢键,提高缓蚀剂膜自身的稳定性,从而形成对腐蚀介质粒子紧密地“包裹”,因此缓蚀剂膜A对腐蚀介质粒子扩散的抑制能力应最强;羟基(OH)和氨基(NH2)也可

21、形成分子间氢键,但是羟基的极性比氨基强,其分子间氢键也比较强30,因此缓蚀剂膜B的缓蚀性能应比缓蚀剂膜C好,而缓蚀剂膜D的性能最差.2.2缓蚀剂膜与腐蚀介质粒子之间的相互作用腐蚀介质在缓蚀剂膜中运移,将不可避免地与缓蚀剂膜发生相互作用,两者之间相互作用强弱将直接影响腐蚀介质在膜中的扩散能力.为此,文中计算了缓蚀剂膜与腐蚀介质粒子之间的相互作用图3缓蚀剂膜与腐蚀介质粒子相互作用能的变化曲线Fig.3inhibitorCurvesmembranesofinteractionandenergycorrosivebetweenparticlescorrosion能(EparticlefilmE),公式

22、如下:particlefilm=(Efilm+Eparticle)-Efilm+particle(2)其中,EEfilm是不包含腐蚀介质粒子的膜体系能量,particle是孤立腐蚀介质粒子的能量,Efilm+particle是包含腐蚀介质粒子的膜体系能量.图3为缓蚀剂膜与腐蚀介质粒子相互作用能的变化曲线.可见,同种腐蚀介质粒子与4种缓蚀剂膜的相互作用能数值均满足:ABCD,这说明亲水支链官能团的极性越强,缓蚀剂膜对腐蚀介质粒子的“包裹”越紧密,越能有效抑制腐蚀介质向金属表面扩散,缓蚀性能也越强,因此4种咪唑啉缓蚀剂的缓蚀性能应为ABCD.此外,从图3可以看出,同种缓蚀剂膜与H3O+和HCO-3

23、的相互作用能远大于与H2O分子的作用能,这可能与两者之间相互作用的方式有关.从分子动力学模拟角度来说,缓蚀剂膜与腐蚀介质粒子之间的相互作用包括库仑作用和范德华作用,而粒子本身的性质将显著影响这两种作用方式的强弱.表3给出了缓蚀剂膜与腐蚀介质粒子的库仑相互作用能(ECoulomb)和范德华相互作用能(EvdW)数值.可知,3种腐蚀介质粒子与缓蚀剂膜之间的库仑相互作用均占主导地位,占两者相互作用能总量的80%以上;对于H3O+和HCO-3来说,粒子均带一个单位的电荷,与缓蚀剂膜之间的库仑相互作用明显高于H2O.因此,与中性粒子(H2O)相比,缓蚀剂膜对带电粒子(H3O+和HCO-3)向金属表面扩散

24、的抑制能力更强,有效阻碍带电粒子在金属表面聚集,减缓腐蚀反应的速率,从而起到缓蚀效果.2.3缓蚀剂膜的自扩散性能缓蚀剂膜的自扩散行为对膜内孔洞结构的稳定No.11乔贵民等：腐蚀介质在缓蚀剂膜中扩散行为的分子动力学模拟3045表3腐蚀介质粒子与缓蚀剂膜的库仑相互作用能(ECoulomb)和范德华相互作用能(EvdW)Table3Coulombinteractionenergies(ECoulomb)andvanderWaalsinteractionenergies(EvdW)betweencorrosiveparticlesandcorrosioninhibitormembranesMembra

25、neBCD2Coulomb-148.8237.1735.10vdW-1-5.82-2.31-5.61825.94786.15775.833+Coulomb-1vdW-1-15.64-24.39-28.72889.43817.07744.98-3Coulomb-1vdW-1-16.7812.6421.76图44种缓蚀剂膜扩散的均方位移(MSD)曲线和自扩散系数Fig.4Meansquardisplacement(MSD)curvesandself鄄diffusioncoefficientsoffourcorrosioninhibitormembranes性具有重要影响,可在一定程度上反映缓蚀剂膜

26、抑制腐蚀介质扩散能力的强弱31.如果缓蚀剂分子的自扩散能力越强,缓蚀剂膜的瞬时孔洞变化越剧烈,所包裹的腐蚀介质粒子的运动也会随之加剧,缓蚀性能则越差.缓蚀剂膜的自扩散性能可通过均方位移(MSD)和自扩散系数曲线来衡量.MSD曲线26定义如下:MSD=Ri(t)-Ri(t0)2(3)根据MSD曲线利用公式(1)可计算出缓蚀剂分子的自扩散系数,缓蚀剂膜的MSD曲线越陡,自扩散系数越大,自扩散能力就越强.4种缓蚀剂膜的MSD曲线和自扩散系数的变化规律见图4.从图4(a)中可以看出,4种缓蚀剂膜的MSD曲线均随时间增加呈上升趋势,这说明由于缓蚀剂分子自身的热运动,膜内的孔洞是不断变化的,可能导致腐蚀介

27、质粒子在膜内的迁移.此外,缓蚀剂膜A的MSD曲线增加幅度最小,说明缓蚀剂膜A内的瞬时孔洞变化最小,膜内孔洞的连通性最差,膜结构最稳定,最不利于腐蚀介质粒子在膜内的扩散迁移,而B次之,再次之是C和D.从图4(b)可以看出,缓蚀剂膜A、B、C和D的自扩散系数分别为4.9×10-12、6.2×10-12、7.7×10-12和15.2×10-12m2·s-1.4种缓蚀剂膜的自扩散系数大小为ABCD,数值逐渐变大,表明膜内缓蚀剂分子的运动性能逐渐增强,膜结构稳定性逐渐降低.这主要是由于亲水支链官能团的极性造成的,官能团极性越强,缓蚀剂分子之间交织作用越剧

28、烈,形成缓蚀剂膜的致密性越好,因此膜内缓蚀剂分子彼此受到的牵制作用越明显,运动迁移能力越弱,缓蚀性能应越好.从缓蚀剂膜的稳定性可以推断出4种缓蚀剂的缓蚀效率顺序为:ABCD.3结论采用分子动力学模拟的方法,研究了4种1R12十一烷基咪唑啉缓蚀剂膜(R1:羧甲基(A),羟乙基(B),氨乙基(C),氢(D)阻碍腐蚀介质粒子(H2O、-H3O+和HCO3)向金属表面扩散的机理,并对其抑制碳钢CO2腐蚀的缓蚀性能进行了理论评价.扩散系数的计算结果显示,缓蚀剂膜均可有效阻碍腐蚀介质向金属表面扩散,从而达到缓蚀效果,且腐蚀介质粒子在不同膜内的扩散系数顺序为ABCD.相互作用能数据表明缓蚀剂膜与腐蚀介质粒子

29、之间存在强烈的库仑作用,同种缓蚀剂膜对带电粒子(H3O+和HCO-3)扩散的抑制能力明显强于对中性粒子(H2O).3046ActaPhys.Chim.Sin.,2010Vol.26缓蚀剂膜的自扩散系数说明,4种缓蚀剂均可形成稳定的膜结构,且稳定性大小为ABCD.综合上述计算结果,4种缓蚀剂缓蚀性能的理论评价结果为ABCD,这与实验结果完全吻合.References1Zhang,T.S.Corrosioninhibitor.Beijing:ChemicalIndustryPress,2008:131-143张天胜.缓蚀剂.北京:化学工业出版社,2008:131-1432Jovancicevic,

30、V.;Ramachandran,S.;Prince,P.Corrosion,1998,55:4493Zhang,G.A.;Chen,C.F.;Lu,M.X.;Chai,C.W.;Wu,Y.S.Mater.Chem.Phys.,2007,105:3314Liu,X.Y.;Chen,S.H.;Tian,F.;Ma,H.Y.;Shen,L.X.;Zhai,H.Y.Surf.InterfaceAnal.,2007,39:3175Zhang,Z.;Chen,S.H.;Li,Y.H.;Li,S.H.;Wang,L.CorrosionSci.,2009,51:291Zhang,J.;Du,M.;Yu,H.H

31、.;Wang,N.ActaPhys.Chim.Sin.,2009,25:525张静,杜敏,于会华,王宁.物理化学学报,2009,25:5257Zhang,J.;Zhao,W.M.;Guo,W.Y.;Wang,Y.;Li,Z.P.ActaPhys.Chim.Sin.,2008,24:1239张军,赵卫民,郭文跃,王勇,李中谱.物理化学学报,2008,24:12398Domiínguez,H.Langmuir,2009,25:90069Kornherr,A.;Nauer,G.E.;Sokol,A.A.;French,S.A.;Catlow,C.R.A.;Zifferer,G.Langmu

32、ir,2006,22:803610Huang,Y.;Liu,Q.L.JournalofXiamenUniversity:NatureScience,2006,45:664黄宇,刘庆林.厦门大学学报:自然科学版,2006,45:66411Zhou,J.H.;Zhu,R.X.;Zhou,J.M.Polymer,2006,47:520612Pan,F.S.;Peng,F.B.;Jiang,Z.Y.Chem.Eng.Sci.,2007,62:70313Yang,J.Z.;Liu,Q.L.;Wang,H.T.J.MembraneSci.,2007,291:114Liu,Q.Z.;Yang,D.F.;Hu

33、,Y.D.Chem.J.Chin.Univ.,2009,30:568刘清芝,杨登峰,胡仰栋.高等学校化学学报,2009,30:56815Tao,C.G.;Feng,H.J.;Zhou,J.;Lü,L.H.;Lu,X.H.ActaPhys.Chim.Sin.,2009,25:1373陶长贵,冯海军,周健,吕玲红,陆小华.物理化学学报,2009,25:137316Pan,F.S.;Ma,J.;Cui,L.;Jiang,Z.Y.Chem.Eng.Sci.,2009,64:519217Liu,X.;Zheng,Y.G.Corros.Eng.Sci.Technol.,2008,43:8718MaterialsStudio.Revision4.2W.SanDiego,USA:AccelrysInc.,200519S