New第八章计算机分子模拟应用实例.ppt

计计 算算 机机 分子模拟分子模拟中国石油大学计算机分子模拟应用计算机分子模拟应用-咪唑啉类缓蚀剂缓咪唑啉类缓蚀剂缓蚀机理的研究蚀机理的研究 研究背景研究背景研究方法研究方法研究对象研究对象研究内容研究内容结结 论论码头上浮船及锚链的腐蚀3水中打捞的摩托车的锈蚀状况4齿轮片的腐蚀表面5铜器表面生成一器表面生成一层薄薄的薄薄的铜绿 铜绿的主要成分是的主要成分是CuCu2 2(OH)(OH)2 2COCO3 3 酸雨腐蚀的文物雕像表面 8四川省宜宾市城区的南门大桥四川省宜宾市城区的南门大桥钢铁钢铁生锈生锈造成大桥造成大桥断裂断裂！铁锈斑斑斑斑某居民楼水管腐蚀状况11管道内壁腐蚀12下水道局部腐蚀状况13管道外壁局部腐蚀状况14研究背景研究背景腐蚀是困扰油气工业发展的一个极为突出的问题添加缓蚀剂是一种有效的防腐手段新型缓蚀剂的设计开发迫切需要理论指导缓蚀机理的实验研究工作的局限性分子模拟技术为深入研究缓蚀机理创造了条件研究对象研究对象CO2腐蚀是石油工业中一种常见的腐蚀类型腐蚀是石油工业中一种常见的腐蚀类型CO2是油气伴生气是油气伴生气注注CO2强化开采工艺强化开采工艺咪唑啉类缓蚀剂咪唑啉类缓蚀剂研究对象研究对象研究方法研究方法量子化学计算量子化学计算分子力学分子力学分子动力学模拟分子动力学模拟开展缓蚀机理的开展缓蚀机理的多尺度模拟研究多尺度模拟研究 以期通过系统的理论研究，丰富和完善缓蚀剂作用理论，并为新型有机缓蚀剂的设计与开发提供理论指导。

咪唑啉类缓蚀剂咪唑啉类缓蚀剂研究内容研究内容1咪唑啉分子反应活性的密度泛函理论研究 2咪唑啉缓蚀剂自组装膜的成膜机理研究3液相条件下咪唑啉在金属表面吸附的MD研究4腐蚀介质在缓蚀剂膜中扩散行为的MD研究1咪唑啉分子反应活性的密度泛函理论研究 研究目标：研究目标：1 1、明确咪唑啉分子的活性、明确咪唑啉分子的活性区域和活性位点分布；区域和活性位点分布；2 2、分析烷基链长对分子反、分析烷基链长对分子反应活性的影响应活性的影响反应活性反应活性全局反应活性全局反应活性局部反应活性局部反应活性全全局局反反应应活活性性化学势化学势()总能量（总能量（E）、电子数（）、电子数（N） 软度软度(S)硬度硬度() 体系总能量体系总能量E E对对N N的一阶导数的一阶导数体系总能量体系总能量E E对对N N的二阶导数的二阶导数亲电指数亲电指数()衡量一个分子的亲电子能力大小衡量一个分子的亲电子能力大小 分子的硬度越小，软度越大、化学势的绝对值越大、分子的硬度越小，软度越大、化学势的绝对值越大、亲电指数越大，表明分子的反应活性越强、越容易得到亲电指数越大，表明分子的反应活性越强、越容易得到（失去）电子失去）电子。

表表2 AH缓蚀剂分子的缓蚀剂分子的化学势化学势()、硬度硬度()、软度软度(S)和亲电指数和亲电指数()Molecule/eV/eVS/eV-1/eVA-1.813.26 0.310.50 B-1.85 3.26 0.310.52 C-1.82 3.26 0.310.51 D-1.85 3.26 0.310.52 E-1.85 3.26 0.310.53 F-1.84 3.26 0.310.52 G-1.803.260.310.50H-1.803.260.310.501.2 1.2 局部反应活性局部反应活性( ( 前线轨道分布前线轨道分布 ) )最高占有轨道最高占有轨道 (HOMO) 最低未占轨道最低未占轨道 (LUMO) MoleculeA4N0.051N0.127N0.164N0.198C0.307N0.18B4N0.041N0.127N0.154N0.198C0.337N0.18C4N0.051N0.127N0.164N0.198C0.307N0.18D4N0.041N0.127N0.154N0.198C0.337N0.18E4N0.051N0.127N0.124N0.198C0.357N0.18F4N0.051N0.127N0.154N0.198C0.327N0.18G4N0.041N0.127N0.164N0.198C0.337N0.18H4N0.051N0.127N0.134N0.198C0.357N0.18表表3 AH缓蚀剂分子的缓蚀剂分子的Fukui指数指数1.3 1.3 反应活性位点反应活性位点 AH分子的全电子密度 Fukui指数可描述分子中各原子的活性强弱定义为电子密度定义为电子密度对电子数对电子数N N的一阶偏导的一阶偏导分别为亲核攻击指数和亲电攻击指数，分别为亲核攻击指数和亲电攻击指数，表示分子中原子表示分子中原子i i得得/ /给电子能力的强弱。

给电子能力的强弱8种咪唑啉分子的反应活性区域均分布在分子的咪唑环及其极性官能团上；其中3个N原子为亲电反应中心，咪唑环上成双键结构的C、N两原子为亲核反应中心分子碳链长度对分子的整体反应活性和活性区域分布基本不产生影响 2咪唑啉缓蚀剂自组装膜的成膜机理研究研究目标：研究目标： 1 1、缓蚀剂膜自身的稳定性以及膜与金属表面的结合强度；、缓蚀剂膜自身的稳定性以及膜与金属表面的结合强度； 2 2、分子烷基链长对自组装膜的影响规律分子烷基链长对自组装膜的影响规律2.1 2.1 吸附表面的选取吸附表面的选取 碳钢表面碳钢表面CO2腐蚀产物膜的结构示意图腐蚀产物膜的结构示意图 碳钢在含碳钢在含CO2环境中腐蚀产环境中腐蚀产物主要是物主要是FeCO3，沉积在金属表，沉积在金属表面形成面形成疏松多孔疏松多孔的的FeCO3膜 模拟过程中选择模拟过程中选择Fe和和FeCO3作为吸附表面作为吸附表面FeCO3Fe2.2 2.2 缓蚀剂分子与金属表面的相互作用缓蚀剂分子与金属表面的相互作用 MoleculeEadsorption/kJmol-1/d/nmFeFeCO3FeFeCO3FeFeCO3A510.63821.197.63 6.46 0.297 0.336B588.41 858.496.71 5.54 0.296 0.350C661.42 894.505.18 1.65 0.299 0.348D714.02 932.234.27 1.80 0.291 0.346E811.03 957.484.98 1.97 0.293 0.326F866.96 981.204.10 1.48 0.296 0.327G968.91 1025.734.06 1.45 0.306 0.331H1021.56 1058.304.00 1.370.289 0.328H2O27.38 267.4133.1727.74 0.295 0.308表4 AH分子在Fe和FeCO3表面的吸附能、取向角( )和距离(d)的统计平均值金属表面金属表面 dFeFeCO3 MoleculeEadsorption/kJmol-1/d/nmFeFeCO3FeFeCO3FeFeCO3A510.63821.197.63 6.46 0.297 0.336B588.41 858.496.71 5.54 0.296 0.350C661.42 894.505.18 1.65 0.299 0.348D714.02 932.234.27 1.80 0.291 0.346E811.03 957.484.98 1.97 0.293 0.326F866.96 981.204.10 1.48 0.296 0.327G968.91 1025.734.06 1.45 0.306 0.331H1021.56 1058.304.00 1.370.289 0.328H2O27.38 267.4133.1727.74 0.295 0.308表4 AH分子在Fe和FeCO3表面的吸附能、取向角( )和距离(d)的统计平均值FeFeCO3 2.3 2.3 缓蚀剂在金属表面的自组装单层膜缓蚀剂在金属表面的自组装单层膜 俯视图俯视图 侧视图侧视图 FeFeCO3 MoleculeECohesive/kJmol-1/ D/nmFeFeCO3FeFeCO3FeFeCO3A23.4228.2131.4736.600.6350.599B28.4530.6739.5143.140.5860.561C40.5243.9656.0961.860.5260.453D50.2655.3057.9764.740.4830.439E65.5169.0960.7864.790.4560.432F78.77105.5359.3467.220.4530.417G92.11119.4460.2867.770.4550.414H112.92133.9564.3068.360.4540.412表表5 AH分子在分子在Fe和和FeCO3表面单层膜的内聚能表面单层膜的内聚能(ECohesive),、吸附角、吸附角() 和链间距和链间距(D)MoleculeECohesive/kJmol-1/ D/nmFeFeCO3FeFeCO3FeFeCO3A23.4228.2131.4736.600.6350.599B28.4530.6739.5143.140.5860.561C40.5243.9656.0961.860.5260.453D50.2655.3057.9764.740.4830.439E65.5169.0960.7864.790.4560.432F78.77105.5359.3467.220.4530.417G92.11119.4460.2867.770.4550.414H112.92133.9564.3068.360.4540.412表表5 AH分子在分子在Fe和和FeCO3表面单层膜的内聚能表面单层膜的内聚能(ECohesive),、吸附角、吸附角() 和链间距和链间距(D)1 1、吸附发生时，分子极性的头部优先吸附于金属表面；、吸附发生时，分子极性的头部优先吸附于金属表面；而非极性的烷基长链则背离表面形成一层密排结构的疏而非极性的烷基长链则背离表面形成一层密排结构的疏水膜。

水膜2 2、随分子烷基链长的增加，缓蚀剂自组装膜自身的稳定、随分子烷基链长的增加，缓蚀剂自组装膜自身的稳定性以及膜与金属基体的结合强度逐步增强性以及膜与金属基体的结合强度逐步增强3 3、同种缓蚀剂在、同种缓蚀剂在FeCOFeCO3 3表面比在表面比在FeFe表面形成表面形成 的自组装膜的自组装膜更加稳定更加稳定 真空条件下的简单模型尚不能真空条件下的简单模型尚不能完全准确完全准确地描述地描述缓蚀剂在金属表面的自组装成膜行为缓蚀剂在金属表面的自组装成膜行为理论评价结论为：理论评价结论为：ABCDEFGH 2.4 2.4 缓蚀性能的理论评价结果与实验结果的对比缓蚀性能的理论评价结果与实验结果的对比 原因：没有考虑溶剂对缓蚀剂吸附行为的影响文献实验结果为：文献实验结果为：ABCDFGH 3液相条件下咪唑啉在金属表面吸附的MD研究研究目标：溶剂对缓蚀剂吸附、成膜行为的影响规律研究目标：溶剂对缓蚀剂吸附、成膜行为的影响规律3.1 3.1 计算模型计算模型 (a)(b)(c)(d)液相条件下缓蚀剂分子在金属表面的初始吸附构型 3.2 3.2 液相条件下单分子在金属表面的吸附液相条件下单分子在金属表面的吸附 水溶液中水溶液中A(7)、C(11)、E(15)和和H(21)分子在分子在Fe 和和FeCO3表面的平衡吸附构型表面的平衡吸附构型 FeFeCO3氨基上的氨基上的N原子易被极性水分子极化成为原子易被极性水分子极化成为(NH3)(NH3)+ +水溶液中水溶液中C(11)分子在分子在Fe 和和FeCO3表面的平衡吸附构型表面的平衡吸附构型 表表6 AH分子在金属表面分子在金属表面吸附能吸附能、溶剂化能和形变能溶剂化能和形变能的统计平均值的统计平均值MoleculeEadsorption /kJmol-1Esolvent /kJmol-1Edeform /kJmol-1FeFeCO3FeFeCO3FeFeCO3A235.79439.69143.0313。