无机盐对C12E5-HPAM复合体系界面性质影响的分子模拟研究论文

1、本科毕业设计（论文）题 目：无机盐对c12e5-hpam复合体系界面性质影响的分子模拟研究学生姓名：学 号：专业班级：材料物理 指导教师:无机盐对c12e5-hpam复合体系界面性质影响的分子模拟研究摘要木论文主要以部分水解的聚丙烯酰胺（hpam）和五甘醇单十二烷基瞇（c12e5）为研究对象，采用分子动力学模拟的方法构建c12e5-hpam复配体系 模型，研究hpam浓度和无机盐种类（cacb、mgc12、nacl）对ci2e5-hpam 复配体系界面性质的影响。研究表明：在c12e5-hpam复配体系中，随着hpam 浓度增大，g2e5和hpam周围水分子数目增多，体系界面水层厚度增加；hp

2、am 与na+之间的静电吸引作用增强，界面双电层的厚度增大。然而，无机盐种类对 g2e5头基周围水分子的分布没有明显影响，但无机盐离子能够穿透hpam竣基 周围第一水化层，且无机盐影响hpam竣基与水分子间相互作用的强弱顺序依 次为 mg2+> na+> ca2关键词：无机盐；表面活性剂聚合物复配体系；界面性质；相互作用effects of inorganic salts on the interfacial properties of c12e5-hpam mixtures by molecular dynamics simulationabstractthis paper ma

3、inly centers on the partially hydrolyzed polyacrylamide (hpam) and pentaethylene glycol monododecyl (cues'by constructing the c12e5-hpam mixtures,we use the molecular dynamics simulation to study the effects of the concentration of hpam and inorganic salts (cacb mgcb nacl) on the interfacial pro

4、perties of mixtures. the results show that the concentration of hpam has effect on the water distribution around the carboxyl group of hpam and the head group of cues. with the increase of hpam concentration, the number of water molecules around c12e5 and hpam increases, and the thickness of interfa

5、cial water layer also increases. in the c12e5-hpam mixtures, inorganic salts has little effect on the water distribution around head groups of ces. however, salt ions can penetrate the hydration layer of hpam carboxyl group, and the order of the influence of inorganic salts on the interaction betwee

6、n hpam carboxyl and water molecules is as follows : mg2+> na+> ca2+.keywords: inorganic salt; surfactant-polymer mixtures; interfacial properties;interaction目录第1章引言11.1表面活性剂简介11.1.1表面活性剂的定义及结构特征11.1.2表面活性剂的分类21.2表面活性剂与聚合物的相互作用31.2.1影响因素31.2.2相互作用机理5第2章理论计算方法或模拟软件介绍62.1分子力学方法62.2分子动力学模拟62.3mate

7、rial studio 软件简介7第3章hpam浓度对ci2e5-hpam复配体系界面性质的影响83.1引言83.2模型构建方法与过程93.3结果讨论分析103.3.1 hpam浓度对体系界面层的影响103.3.2 hpam浓度对体系水分子分布的影响123.3.3 hpam浓度对体系界面双电层的影响153.4本章小结17第4章无机盐对c12e5-hpam复配体系界面性质的影响184模型构建方法与过程184.2结果与讨论184.2.1无机盐种类对体系界面层的影响184.2.2无机盐种类对c12e5的影响204.2.3无机盐种类对hpam的影响224.2.4无机盐种类对组分扩散行为的影响254.3

8、木章小结29第5章结论30致谢31参考文献32第1章引言三次油cei开釆技术，亦可称为“强化采油提高采收率技术”，主要利用一系 列的生物，物理和化学等方法开采油田中不动油的技术。当今已形成的三次开 采技术主要包含化学驱、热驱、气驱和微牛物驱以。中国作为能源大国，人民对 能源的需求不断增长，提高原油釆收率势在必行。我国在复合驱、蒸汽吞吐和聚 合物驱的应用上远超其他国家，其中聚合物驱主要应用在大庆油田，二元复合驱 和高温高盐聚合物驱已在胜利油实现大规模应用。聚合物驱能够有效地改善 流度比，表面活性剂驱在改变残余油和表面张力方面更具优势，但是由于聚合物 驱和表面活性剂驱相对于气驱和热驱的成木更高，并

9、没有得到大规模地推广。 因此聚合物-表面活性剂复合驱具有更高的可行性，一方面聚合物具有粘弹性， 能够增加驱替液的粘度，控制流度，提高波及范阖；另一方面表面活性剂能够降 低界面张力，改善微观驱油效率，二者的协同作用能够极大地改善采油效率并且 降低化学剂成本。由于实际实验条件的限制，对于聚合物-表面活性剂驱的研 究十分有限，因此木论文采用md方法研究ci2e5.hpam复配体系的微观性质， 探究体系中各组分的作用机理以及无机盐和不同浓度hpam对c12e5.hpam复配 体系的影响。1.1表面活性剂简介1.1.1表面活性剂的定义及结构特征表面活性剂是一种能够以极低的浓度将溶剂的表面张力显著降低的物

10、质，其 单体的基本组成是具有极性的亲水基（与水等极性物质相亲）和不具有极性的亲 油基（与油等非极性物质相亲），这样的分子组成决定它具有双亲性，但是又不 会完全亲水或亲油，其亲油亲水性质的平衡影响它的微观性质。一般而言，亲 油基为长绘链，按照怪链中碳的个数，可以将其大致分为水溶性和油溶性。若碳 少于12个，则表现为水溶性，亲水基能够把整个分子拉入水中。若碳大于14 个，则表现为油溶性。作为一种化学驱油剂，它不仅具有应用范围广，提高油藏 采收率幅度大，发展潜力大的优点，而ii能够改变油层的润湿性，显著降低油水 界面张力。由于结构中同时具有亲水基和亲油基，这样的组成使得它能够在油水界而处 定向排列，

11、形成独特的吸附层，进而降低溶剂的界面张力。它溶于水后，亲水 基受到水的吸引甚至会将炷链亲油基拉入水溶液中，而亲油基由于被水排斥，有 逃离水溶液的趋势。为了维持体系的稳定状态，表面活性剂就会吸附在溶液的界 面，使得亲油基向气相延仲，亲水基向水相延仲，形成定向单分子层吸附，导致界 而张力下降,表现出表而活性。这种特殊的分子组成使其易于富集在表而，但 是能够富集在表面的含量仍存在一个最大的限度，即临界胶束浓度（cmc） ,01, 大于cmc,它就会在溶液中形成胶束，不能继续溶于水中，此时再增加它的浓度 也不会改变溶液的表面性质，不能够再继续降低表面张力。1.1.2表面活性剂的分类通常根据它溶于水后，

12、生成离子与否及其极性基电离特性，将它划分为离子 和非离子两大基本类型。（1）阴离子表面活性剂它是用途最广泛，使用量最大的一类表面活性剂,约占总产量的40%以上】, 溶于水后会电离出带负电的阴离子活性基团，根据基团的特性不同，又可将其分 为几种不同类型：具有乳化和洗涤作用的r-coona,即竣酸盐型；具有良好洗涤 性和发泡性的r-oso3na,其可以在硬水屮稳定存在，使溶液显屮性和弱碱性； 根据不同的烷基链长可表现出不同的活性，这一类型是r-so3na,用途更加广 泛，可作发泡剂，乳化剂，洗涤剂等；另外一种是r-oposna,具有抗静电和增 稠作用i。在三次采油中磺酸盐类比较常用，主要包含石油磺

13、酸盐，烷基磺酸 盐和烷基芳基磺酸盐等w（2）阳离子表面活性剂这一类型的显著特点是其溶于水电离后，亲水基带正电，它一-般不会用作洗 涤剂，主要因为它易与基质上的油污发生作用，不具有洗涤功效。除此之外，它 的用途广泛，可以抗静电，用作护发素和煽油膏等。目前常见的阳离子表面活性 剂仍是含氮化合物i，按照氮的位置可将其划分为三种类型，分别是杂环型， 胺盐和季钱盐w（3）两性离子表面活性剂结构中具有彼此不可被电离的阴、阳离子电荷中心，其性质与阴、阳离子表 而活性剂均不同。普遍意义上，它可以按照阴阳离子电荷中心氮原子分成两 类，分别是含弱碱性氮和含强碱性氮的表面活性剂"i。它具有许多优良的性质,

14、 如乳化性，低毒性，低刺激性，抗静电性和抗菌抑霉等。(4) 非离子表面活性剂它在水中不发生电离，亲水基以-0h或醞键为主，醯基中的氧及0h通过与 水中的氢形成氢键而产生亲水性9201。一般-oh或瞇键要达到一定的数量才能 够发挥亲水基的亲水性，一个-oh或是醯键的亲水性太弱。它能够在溶液屮稳 定存在，耐强电解质、强酸、强磁 并能够与其他物质混合使用，具有相容性好 的优点，相较于其他类型，它具有更高的乳化能力，以及更好的润湿性和除污洗 涤性能】。1.2表面活性剂与聚合物的相互作用在二元复合驱中，表面活性剂与聚合物存在相互作用，能够改变聚合物的构 象，使其能够更加容易吸附在固态物质的表面，改变界面

15、行为；也可以改变表面 活性剂的物理化学性质昭，即使在低于cmc的情况下，也能够增加溶解度；二 者的相互作用会导致整个体系的粘度，表面张力，界面吸附行为等性质发生改变 23o通常情况下，它们的相互作用可能是由氢键，疏水作用，偶极作用和库伦力 产生的23,受表面活性剂或聚合物的分子结构，电荷等因素的影响。一般能够 与表面活性剂发生相互作用的聚合物主要是水溶性的，包括聚电解质，水溶性非 离子均聚物和疏水缔合的非离子共聚物】,二者可能会产生协同效应，也会产生 非协同效应，这主要取决于两种物质的类型以及其他相关因素。1.2.1影响因素影响因素大致分为两类，外部因素和内部因素，其中内部因素起决定作用， 包

16、括两种物质本身的性质，而外部因素主要有温度，无机电解质，溶液的酸碱度 和氢键等。(1)温度通常，温度越高，相互作用越弱，越不容易发牛缔合，离子型表面活性剂的 cmc就越大。一般来说，温度对表面活性剂的cmc具有复杂的影响卩叫在很大 的温度范围内，ch1c是随着温度的升高先减小后增大3】，在某一温度时达到最低 cmco这主要是因为温度对亲水基和亲油基的影响存在差异。温度越高，胶束越 容易形成，不利于亲水基的水合，而高温会破坏亲油基周围水层的形成，不利于 胶束的形成，因此cmc取决于温度对这两种基团作用的综合影响。低温对亲水 基的水合产牛较大影响；高温主要对亲油基的疏水性产生较大影响。此外，当达

17、到浊点时，中性聚合物的疏水性提高，更容易与表面活性剂发生相互作用绚。(2) 无机盐体系屮加入的无机盐能够起促进作用，增加二者发生相互作用的概率0】。 无机盐能够促进离子表面活性剂和中性聚合物的结合，然而聚电解质与带相反电 荷的表面活性剂的作用会被无机盐屏蔽，相互作用减弱12叫本论文中研究的g2e5溶于水后，瞇键中的氧与水中的氢通过氢键作用实现 亲水效果，随着温度升高，氧键逐渐断裂，溶液开始岀现浑浊，即g2e5达到浊 点温度，开始由完全溶解转化为部分溶解，出现浑浊现象。浊点常用來表征非离 子表面活性剂的亲水性，亲水性越好，浊点越高。无机盐会对浊点产生影响，不 利于氢键形成而造成“脱水”现象，降低

18、溶解度。无机盐对hpam的柔顺性和粘度等性质具有显著的影响。柔顺性是指高分 子链能够改变其构象的能力，无机盐与竣基通过静电吸引产生相互作用，能够在 竣基周阖形成无机盐离子层，屏蔽hpam ±的竣基间的相互排斥，致使分子链 收缩，溶液的黏度也会降低。(3) 表面活性剂的影响它本身的结构和性质会影响其与聚合物的相互作用，主要包括以下儿个方 面：链长。通常情况下，分子链越长，cmc就越小，能够在较低浓度时发生作 用。链结构。链的结构不同，二者相互作用也存在差异。比如，sds的硫酸碳 氢链中加入氧乙基后，sds与pvp的相互作用相对减弱种类。如非离子 型与中性聚合物几乎不发牛或只存在很弱的相

19、互作用卩气极性基团所带电荷 的正负及电荷量悶。浓度。未达到cmc时，它以单分子在体系中存在，与聚 合物几乎不发生相互作用。浓度增加至cac时，它开始吸附在聚合物上并发生相 互作用，当达到cmc时，它在溶液屮逐渐形成胶束，与聚合物的作用饱和，此 后再增加它的浓度也不会改变作用的强弱。(4)聚合物的影响分子量。当分子量达到某一确定值时，其与表而活性剂的作用达到最强， 大于或小于这一确定值，都会减弱相互作用。浓度。浓度对cac几乎没有影响, 但却与作用达到饱和吋的浓度成正比关系w疏水性，电性和柔顺性都会产 牛影响，且强弱顺序是疏水性电性柔顺性圖01.2.2相互作用机理二元驱中所用聚合物和表而活性剂是

20、具有一定的性质的，并不是所有类型都 可以发生作用，而作用强弱也取决于两者类型以及所带电荷，分子结构，所带基 团种类等因素卩2】。高分子量的水溶性化合物，由于分子量较大，在溶液体系中 通常会达到胶体尺寸大小，能够跨越很大长度与溶剂分子发牛相互作用，提高其 在溶液中与其他分子的缔合能力。它在溶液体系中占据空间尺寸较大且不交叉， 因此增加分子量和浓度可以提高它们的缠结倾向1舛。二者作用主要是由三种作用力引起的，分别是静电力，分子间作用力和疏水 作用力，其中疏水力和静电力占主导【绚，实际上是受静电力、分子间作用力和 疏水作用力三者的复杂影响，能够影响这几种力的因素都能够改变作用力的大小 l22jo聚电

21、解质与带相反电荷的表面活性剂主要受静电力的驱动，两种物质的离 子相互影响，而每个吸附在聚电解质上的离子都能够成为离子簇的核心必。聚 电解质的类型和所带电荷决定相互作用的强弱，电荷越多，簇化及发生的作用就 越强。聚电解质与带同种电荷的表面活性剂存在强烈的静电排斥作用，它们几乎 不会发牛缔合。由于非离子表面活性剂分子间不存在静电排斥作用，在溶液中很容易形成胶 束，cmc值一般要低于离子型的cmc值，同时它的分子结构中又含有较多的氧乙 基基团，体积较大会屏蔽疏水相互作用，因此非离子表面活性剂与聚合物几乎不 会发生或仅存在很弱的作用由于疏水作用力的驱动，表面活性剂会优先以 胶束聚集体的形式结合到聚合物

22、链上，形成类似于“聚电解质”复合物i绚，形 态类似于一条聚合物骨架上带有胶束的珍珠项链a】，即串珠模型。第2章 理论计算方法或模拟软件介绍分子模拟是基于量子力学理论的一种计算机模拟方法，通过对某-种特定体 系或是过程建立一种数学模型以简化表达，计算或是预测出分子体系的行为或性 质。分子动力学模拟已经发展成为一种成熟的技术，可以用来有效地理解大分子 的结构与功能的关系。本章中主要介绍分子动力学以及分子动力学在模拟 c12e5-hpam复配体系中的应用。2.1分子力学方法分子力学方法（molecularmechanics）以经典力学为基础来模拟分子体系。 所有原子的分子力学模拟方法都有如下的特点：

23、每一个原子都被模拟成一个微 粒；每一个原子都被分配了一个半径范围（特别是范德瓦尔斯瓦半径范围，即化 合键结构中范德瓦尔斯引力和斥力处于平衡的区间）；极性和恒定净电荷（通常 是来自于量子计算或是实验）；键与键的相互作用，也相当于弹簧上的作用，它 的平衡距离等于实验或理论计算出来的键长。化学体系屮势函数或力场可以抽彖为（21）的函数：e二ecovalcnt+enoncovalcnt（ 2-1 ）式中：e为总化学势，ecovalent为共价分子化学势，enoncovaienl为非共价分子化学势。通过（21）可以计算出整个分子体系的势函数，其相当于体系中存在的每 一组成部分各自的势能的总和。其中共价分

24、子（或原子）和非共价分子（或原子）(2-2)(2-3)对于整个势能的贡献可以通过（2-2）和（2-3）的表达式给出:ecovalent' -eb()nd+eanglc+e(jihcclnilenoncovaienl-eclcctrosiaiicevandcrwaals式屮：eb°nd为键与键的化学势，eangc为键角化学势，edihcdral为二面角化学势。eelectrostatic为静电化学势， evanderwaais 为范德瓦尔斯化学势具体的势函数或是力场的函数形式却取决于使用的某种具体的模拟方法。2.2分子动力学模拟原子和分子的运动轨迹是通过数值求解牛顿运动方程来确

25、定的，其中粒子的 相互所用力以及它们的势能是利用原子间势能或分子力学力场来计算的。基于x 射线晶体学或nmr光谱学的实验限制，分子动力学模拟经常用于改进蛋口质和 其他大分子的三维结构。在物理学领域，分子动力学模拟也用來观察直接用肉眼 难以观察到的原子尺度的动力学现象。然而设计出分子动力学模拟方法必须要考 虑计算机的计算能力，也要考虑模拟的范围大小（粒子的数目），时间步长和总 的持续时间，以保证能够在合理的时间段里完成模拟计算。在分子动力学模拟中, 要注意选取可以处理稳定系统的合适系综，不同的系综处理分子动力学模拟问题 的应用范围不同。2.3 material studio 软件简介materi

26、al studio是用来模拟和构建材料模型的一种软件，它是由biovia （之 前的accelrys）研发与公布的，这个软件被各大高校，研究屮心和高科技公司所 采用，广泛应用于计算和模拟各种各样先进材料，比如聚合物，碳纳米，晶体， 金属，陶瓷等等。material studio的基木工作流程是：1.materials visualizer用来 创建或导入材料的图像模型。2量子力学，半经验理论或是经典模拟方法决定精 确的材料结构。3.分析或预测各种各样的所需要的性质。本论文采用material studio 6.0软件，首先利用visualizer模块构建一系列 c12e5-hpam复配体系模型

27、，然后通过discover模块对c12e5-hpam复配体系模 型进行优化和动力学计算，最后通过forcite模块分析模型，得到与c12e5-hpam 复配体系界面性质相关的数据，实现对c12e5-hpam复配体系界面性质的研究。第3章hpam浓度对c12e5-hpam复配体系界面性质的影响3.1引言c12e5是一种由亲水极性醯键r-o-r和疏水非极性烷基组成的非离子表而活 性剂,它的结构屮含有足够数量的亲水基(-r-o-r)使其在溶液体系屮具有良好的 亲水性。由于g2e5结构中仅在分子的端基存在一个oh基团，瞇键的数量弥补 了0h亲水性小的不足，因此c12e5在水中和某些有机溶剂中亲水性好，

28、稳定性 高，并且不易受强酸强碱影响。hpam是一种结构中除了含有极性基团酰胺基 (-co-nh2),还含有竣基(cocr)的阴离子型聚合物，相对于pam, hpam 能够更好地改善水的黏度，提高采收率。本文主要是从微观角度探究c|2e5-hpam复配体系的界面性质，采用md方 法探究不同hpam浓度和不同种类无机盐对c12e5-hpam体系界面性质的影响, 首先要构建c12e5和hpam分子模型，如图31 (a) (b)所示。(b)图3-1(a) c12e5分子构型图(b) hpam分子构型图fig. 3-1 (a) model of pentaethylene glycol monodode

29、cy (c12e5)(b) model of partially hydrolyzed polyacrylamide(hpam)3.2模型构建方法与过程木文釆用material studio 6.0软件构建ci2e5.hpam复配体系模型，设置周期 性边界条件，将模型建立为镜面对称的g2e5单层膜结构，即在c12e5_hpam体 系模型中，形成g2e5一层-水层-g2e5一层的三层结构，而hpam链分别处于水层 与g2e5一层界面处，且水层屮分散着一定数量的na4"以保持电屮性。首先构建分子模型，利用materials visualizer子模块，分别构建c12e5分子, hpam分

30、子和水分子模型，并利用discover子模块中的setup对建立的分子 分配力场，以及【minimized对分子进行初步优化，使分子能量最小化，消除 在构建体系模型时可能造成的分子间发生重叠。建立盒子，为保证©止5能够完 全处于盒子z屮，要建立合适尺寸的g2e5单分子层和水盒子，其中g2e5分子 层含有36个c12e5分子。利用build layer构建三层体系结构模型，即c12e5.层- 水层-c12e5-层，同时各有30 a的真空层位于c12e5-层的两侧。如图32所示，其 中红色是水分子，绿色是hpam链，蓝色是nas白灰色是g2e5烷基链，红白 色是c12e5头基。图3-2

31、c12e5-hpam复配体系的初始构型fig. 3-2 the initial model of c12e5-hpam mixtures本章主要研究不同hpam浓度对c12e5-hpam复配体系的影响，各体系中 分别含有一条hpam链，两条hpam链和三条hpam链。建立体系模型并进行 参数设置以及动力学模拟后，可以得到含有不同浓度hpam的c12e5-hpam复 配体系的平衡构型图，如图33所示。(a) 一条 hpam 链(b) 两条hpam链 (c)三条hpam链图3-3不同浓度hpam体系中ci2e5-hpam复配模型的平衡构型图fig. 3-3 balanced models of c

32、12e5-hpam mixtures in systems with different concentrationsof hpam3.3结果讨论分析 331 hpam浓度对体系界面层的影响为了分析动力学模拟后体系中各组分的分布，分别做各体系中g2e5分子, hpam分子，g2e5头基，c12e5烷基链，水分子和na*沿z轴的密度分布曲线, 如图3-4所示。140160r-uowj、a=suuqh2o head chainci2e5hpamo0 8 6 4 2 081 00000s2uoob>.esuocl20looz/a140160(a) cpes-lhpam(b) cpes-2hpa

33、m6()nz/a(c) ci2e5-3hpam图34不同浓度hpam体系中各组分沿z轴的密度分布曲线fig. 3-4 the density curves of different components in systems with different concentrations ofhpam图3-4是不同hpam浓度体系中各组分沿z轴的密度分布曲线。由图3-4 可知g2e5在体系中的分布发生明显变化，c2e5头基几乎完全进入水层中，c2e5 烷基链大部分位于空气屮，在图屮表示为g2e5头基的曲线几乎全部位于水分子 曲线范围内。由图中g2e5头基和hpam的密度分布曲线可知，二者的密度分布

34、 曲线范围大部分重合，说明g2e5头基与hpam存在较强的相互作用。各个体系 中水的密度均接近真实密度，并在体系中占据很宽的界面分布，存在着明显的界 面水层，在密度分布曲线中表示为水密度逐渐减小的部分。形成界面水层的原因 在于体系屮水分子间是通过氢键结合的，在逐渐达到界面处水分子间的氢键作用 减弱而使得密度减小。为了能够详细地了解ci2e5-hpam复配体系的界面层，采用"10%-90%厚度 原则”划分界面厚度，即水相密度的10%到90%的范围，如图34中红色虚线范 围所示，计算不同浓度hpam体系中厚度如表31所示。由表可得，随着hpam 浓度增大，体系的界面水层厚度从11.224

35、 a增大到12.219 a,即界面水层厚度随聚 合物浓度增加而增大。表3-1不同hpam浓度体系的界面水层厚度table 3-1 tickness of interfacial hydration layer in systems with different concentration of hpam体系c12e5-ihpamc12e5-2hpamc12e5-3hpam界面水层厚度/a11.22411.23312.2193.3.2 hpam浓度对体系水分子分布的影响图3-5 g2e5头基屮eo5基团与水屮氧原子的rdf曲线fig. 3-5 rdf curves between eo5grou

36、p of c12e5 headgroup and oxygen atom of water molecules表3-2 c12e5头基中eo5周围第一水化层内水分子配位数table 3-2 number of h2o in the first hydration layer of eo5group of c12e5 headgroup体系ci2e5-ihpamc12e5-2hpamc|2e5-3hpam水配位数0.5110.5460.869径向分布函数（rdf）可以描述指定分子周围其他粒子的空间分布情况。图 3-5是不同浓度hpam体系屮c12e5头基屮eo5与水中氧的rdf曲线，表3-2 是

37、e0基团周围水配位数。由图35可知，当体系中加入一条hpam链和两条 hpam链时，rdf曲线上的峰值和eo基团周围水配位数变化不明显。当体系中 加入三条hpam链时，rdf曲线上第一个峰的峰强明显增大，eo5基团周圉水 配位数明显增多，即g2e5头基周围结合水分子数目随hpam浓度的增大而增加。 这可能是由于随hpam浓度增大，具有更强亲水能力的hpam和g2e5间的相 互作用增强，导致c12e5头基和水分子间的相互作用增强，使得c12e5头基屮eo5 基团周围聚集更多的水分子。54u 32100246810r/a图3-6 hpam竣基与水中氧原子的rdf曲线fig. 3-5 rdf cur

38、ves between carboxyl of hpam and oxygen atom of water molecules表3-3 hpam竣基周围第一水化层内水分子配位数table 3-2 number of water in the first hydration layer of carboxyl of hpam体系c12e5-ihpamci2e5-2hpamc|2e5-3hpam水配位数0.7180.9130.968图3-6是hpam竣基与水中的氧原子的rdf曲线，表33是竣基周围水分子配位数。由此可知，随着hpam浓度增大，hpam周围水分子的数目从0.718 增大到0.968,

39、即hpam周围结合水分子数目随hpam浓度的增大而增加。这主 要是因为随着hpam浓度增大，带负电荷的竣酸根离子数冃增多，hpam与水 分子之间的相互作用增强，导致hpam周围聚集更多的水分子。为了进一步探究hpam浓度对hpam竣基周围水分子分布的影响，可以做 hpam竣基周围水分子的均方位移曲线。均方位移曲线（msd）可以用来研究 体系中不同组分的扩散行为，进而说明影响其扩散行为的因素以及所受作用力大 小。扩散系数是描述粒子扩散能力的物理量，其与粒子的msd曲线斜率存在如 （3-1）对应：” nd = limsk -耳（0）（3-1）6n w dt =其中微分项是msd曲线的斜率，n是参与

40、扩散的粒子个数，扩散系数是msd曲线斜率的1/6。图3-7 hpam竣基周围第一水化层中水分子的msd曲线fig. 3-7 msd curves of water in the first hydration layer of carboxyl of hpam表3-4 hpam竣基周围第一水化层中水分子扩散系数table 3-4 diffusion coefficients of water in the first hydration layer of carboxyl of hpam体系c12e5-ihpamc12e5-2hpamc|2e5-3hpam扩散系数(10'4cm2/s)

41、0.3610.3390.338图3-7是hpam竣基周圉水分子的msd曲线，表34是hpam竣基周围水 分子的扩散系数。由此可知，体系中hpam浓度越大，竣基周围水分子的扩散 系数逐渐减小（0.361、0.339、0.338）o这是因为随着体系中hpam竣基数目增多, hpam能够结合更多的水分子，同时对竣基周围水分子的束缚作用加强，使水的 扩散受到限制。由图3-5和图3-6可知，导致不同浓度hpam体系中界面水层厚度不同的原 因在于，随着体系中hpam浓度增大，g2e5与水分子相互作用增强，c12e5头 基结合的水分子数量增多；同时hpam周围结合水分子数目也会增加，两者同 时作用导致界面水

42、层厚度增加。3.3.3 hpam浓度对体系界面双电层的影响hpam在水溶液中会形成双电层，双电层主要包括两个部分，一部分是在 hpam表面上的离子固定吸附层，即stem区；另一部分是离子扩散分布的diffuse 区。为了进一步了解hpam与nf的相互作用，分析不同浓度hpam体系中 hpam竣基与na*的rdf曲线，积分求得hpam竣基周围na*配位数，如表3-5, 同时得到双电层的stern层和diffuse层厚度，如表36。图3-8 hpam竣基与na'的rdf曲线fig. 3-8 rdf curves between carboxyl of hpam and sodium ion

43、表3-5 hpam竣基周围na*配位数table 3-5 number of sodium ion of carboxyl of hpam体系ci2e5-1hpamci2e5-2hpamc|2e5-3hpam配位数7.069.6012.89表36不同hpam浓度体系双电层厚度table 3-6 tickness of electric double layer in systems with different concentration of hpam体系c12e5-ihpamc12e5-2hpamc|2e5-3hpamstern 层/a1.021.0414diffuse 层/a1.302.

44、382.42图3-8是不同浓度hpam体系中hpam竣基与na*的rdf曲线，表3-5是 根据rdf曲线得到的hpam竣基周围na*配位数，表3-6是根据rdf曲线得到 的双电层厚度。由图38可知，在不同浓度hpam体系中竣基与nf的rdf曲 线上，第一峰位均位于2.04a左右，竣基与nf发生强的静电吸引作用，但是竣 基周围nf配位数不同，而且不同体系的界面双电层的厚度也不同，如表35和 表3-6所示。由此可知，随着体系中hpam浓度增大，hpam竣基周r1 na+配位 数增大，界面双电层的厚度增大。这主要是因为，随着hpam浓度增大，体系屮 竣基数目增多，与nf的静电相互作用增强，使双电层的

45、厚度增加。由于体系中 nf与hpam竣基存在强烈的静电吸引作用，这必然会对体系中nf的扩散行为 产生影响，如图39是hpam竣基周围na*的msd曲线，并再根据msd曲线 得到不同浓度hpam体系中n£的扩散系数，如表3-7所示。图3-9 hpam竣基周围na4的msd曲线fig. 3-9 msd curves of sodium ion of carboxyl of hpam表3-7 hpam竣基周围na*扩散系数table 3-7 diffusion coefficients of sodium ion of carboxyl of hpam体系c12e5-ihpamc12e5-

46、2hpamc|2e5-3hpam扩散系数(10'5cm2/s)0.1780.1150.105图3-9是不同浓度hpam体系中hpam竣基周围na*的msd曲线，表3-7 是根据图3-9得到的na*扩散系数。由表37可知，随着体系中hpam浓度增大, hpam竣基周围nf扩散系数逐渐减小，这主要是因为hpam竣基与nf存在强 的静电吸引作用。体系中hpam浓度越大，与nf的静电相互作用越强，对na* 扩散行为的束缚能力越强，导致hpam竣基周围nf扩散系数越小。3.4本章小结本章主要采用分子动力学模拟的方法研究hpam浓度对c12e5-hpam复配 体系界面性质的影响，通过分析复配体系平

47、衡构型、界面水层结构、界面双电层 结构等，可以得到以下结论：(1) 在c12e5-hpam复配体系中，hpam浓度越大，界面水层厚度越大。(2) 随着体系中hpam浓度增大，c12e5头基和hpam竣基周围结合水分 子数增多，且hpam周围结合水的扩散系数减小。(3) 随着hpam浓度增大，hpam与nf的静电相互作用增强，界面双电 层厚度增大，且na+的扩散系数逐渐减小。第4章 无机盐对c12e5-hpam复配体系界面性质的影响4.1模型构建方法与过程木章主要探究无机盐种类(cacb、mgcl2> nacl)对ci2e5-hpam复配体系 界面性质的影响。在复配体系中添加不同种类的无机

48、盐，分别是36个cacb, 36个mgcb和72个nacl,其中hpam链均为一条，其他分子的种类及数量均 与第三章屮相同。体系经过动力学模拟后，可以得到含有不同种类无机盐的 ci2e5-hpam复配体系的平衡构型图，其中绿色代表c&2+,蓝色代表mg2+,紫色 代表nf,如图41所示。(a)加入c产(b)加入mg*(c)加入nf图41加入ca2+ ,mg2+ ,na+的ci2e5-hpam复配体系的平衡构型图fig. 41 balanced models of c12e5-hpam mixtures in systems with calcium ion,magnesian iona

49、nd sodium ion4.2结果与讨论4.2.1无机盐种类对体系界面层的影响(a)未加入无机盐0.4°%1.00.80.6().2h2o-head 一一chainc12e5-hpam100120z/a140160(c)加入 mg?*图4-2未加入无机盐和加入ca2+、mg1.0density/g.cm p0005i-bsdo av : th2o head chain c12e5 hpamna+80100120140160z/a(d)加入nf、nf体系中不同组分的密度分布曲线(b) ijq a ca2+fig. 4-2 the density curves of different

50、 components in systems with no ion and calciumion,magnesian ion,sodium ion图42是c12e5-hpam体系中各组分沿z轴的密度分布曲线，其中各组分分 别为水分子，c12e5头基，c12e5烷基链，c12e5分子，hpam分子和添加的无羽l 盐离子。由图41和图42可知，体系经过动力学平衡后，g2e5的分布相对于 初始构型发牛变化，界面无序性增加，g2e5头基伸向水相，g2e5非极性烷基链 由于具有疏水作用，大部分伸向空气中，只有一小部分能够伸向水相。为了更加详细地了解无机盐种类对c12e5-hpam复配体系界面层的影响，

51、分 别计算含有c0、mg2+ . n"体系的界面水层厚度，如表41所示。由此可知, 与未加入无机盐的体系相比(11.22 a),加入ca2 mg2+、na*体系的界面水层 厚度明显减小(9.77 a、8.53 a、9.31 a), >1界面水层的减小程度由大到小依次 为 mg2+> na+> ca2%表4-1未加入无机盐和加入ca2+,mg2+,na+体系的界面水层厚度table 41 tickness of interfacial hydration layer in systems with no ion and calciumion,magnesian ion

52、,sodium ion.体系加入ca2+加入mg2+加入na*未加入无机盐界面水层厚度/a9.778.539.3111.224.2.2无机盐种类对g2e5的影响图43 c|2e5头基屮eo5基团与ca2+、mg2 na4的rdf曲线fig 4-3 rdf curves between e05 group of c12e5 and calcium ion,magnesian ion,sodium ion.图43是未加入无机盐和加入ca2+ , mg2+, na*的体系中，g2e5头基中eo5 基团与ca2+ , mg2+, na*的rdf曲线，其屮eo5基团是与c12e5屮端基拜基距 离最近的氧

53、乙基(-occ-)基团。在未加入无机盐和加入na*的体系中，第一个 峰值出现在2.53a处，第二峰值出现在4.73a处。而加入mg2 c产的体系中, 第一个峰值分别出现在4.29a和4.53入处。这是因为加入na+,体系中nf的浓 度增大，与c12e5头基的离子-偶极作用增强，所以峰位较未添加无机盐体系没 有改变，只是峰值略有増加。而体系屮加入mg2ca2+后，使g2e5头基与体系 中的nf之间的相互作用减弱，所以第一个峰值消失。根据第二峰值大小可知， 与g2e5头基发生相互作用的强弱顺序是mg2+> ca2%原因在于，与c/+相比， mg?+具有更小的离子半径和更小的质量，因此在加入c

54、a2+ , mg?+的体系中，mg2+ 在体系中更加容易发生移动，与g2e5头基之间的离子-偶极作用更强。c12e5头基屮的eo (-occ-)基团具有极性，体系屮的水会在c12e5头基周 围形成水化层，并将rdf曲线上第一个峰称为第-水化层，并将第一个峰的起始位置到第一个峰谷位置的距离跨度称为第一水化层的厚度。图4-4是体系中加 入ca2+ . mg2+> na+后，g2e5头基与水分子中氧的rdf曲线。r/a图44 c12e5头基与水屮氧原子的rdf曲线fig. 44 rdf curves between g2e5 head group and oxygen atom of wate

55、r molecules从图44可以看出，无论体系中是否加入c/s mg2 na+, g2e5头基与水 屮的氧原子的rdf曲线几乎完全一样，说明不同种类的无机盐对c12e5头基与 水的相互作用影响不大，这主要是由于g2e5对盐不敏感，受盐影响较小。为了 进一步说明g2e5与水的相互作用，将各体系中g2e5中的五个e0基团分别与 水中的氧原子作rdf曲线，如图45所示，其中e05, e04, e03, e02, e01 分别代表与g2e5端基疑基最近邻，次紧邻，依次距离渐远的氧乙基基团。051()15202530r/a(b)加入mg"(a) in a ca2+(3)加入na+(4)未加入

56、无机盐图4-5各eo基团与水中氧原子的rdf曲线fig 45 rdf curves between every eo group of c12e5 and oxygen atom of water molecules图45是未加入无机盐和加入ca2+ . mg2 na+的复配体系中，各e0基团 与水中的氧原子的rdf曲线。从图45可以发现，e0基团与c12e5端基-0h越 远，与水的相互作用越小，其中最近邻的e05基团与水的相互作用最强，rdf 曲线上出现较为明显的两个峰，而与端基oh距离最远的eo1基团与水相互作用最弱。这主要是因为g2e5烷基尾链的疏水性导致g2e5有脱离水的倾向，会 减弱

57、各e0基团与水的相互作用。eo1基团与疏水尾链距离最近，e05基团距离 最远，所以e01与水的相互作用最弱，而e05与水的相互作用最强。4.2.3无机盐种类对hpam的影响图4-6 hpam竣基与ca2+、m+、na*的rdf曲线fig. 4-6 rdf curves between carboxyl of hpam and calcium ion,magnesian ion,sodium ion.图4-6是加入ca2+ >mg2+>na+的复配体系屮，hpam竣基屮氧原子与ca2+、 mg*、na*的rdf曲线。由图46可知，在加入ca2+ > mg2 na*的体系中， rdf曲线上存在两个明显的峰，竣基周围ca2+ > mg2 nf的第一峰位分别为 2.07入、1.69