温度影响石墨烯沥青复合材料界面微结构特征的分子动力学模拟

　　黄立新 周昆 黄君 叶太智广西大学土木建筑工程学院 广西大学工程防灾与结构安全教育部重点实验室
　　摘 要： 采用分子动力学模拟方法研究了温度对石墨烯/沥青复合材料界面微结构特征的影响。选取沥青的4组分模型，进行分子动力学模拟，得到合理的沥青玻璃化转换温度。然后加入石墨烯，共同搭建石墨烯/沥青复合材料模型并进行分子动力学模拟。在此基础上，基于复合材料切片的相对浓度，定义界面区间的范德华间隙和致密区间。根据复合材料相对浓度曲线，计算得到范德华间隙厚度，然后采用统计学方法，得到沥青基体累计浓度以及一阶导数的曲线，计算得到致密区间厚度。在各种温度条件下，通过分子动力学模拟得到石墨烯/沥青复合材料界面区间范德华间隙和致密区间的厚度。研究结果表明，温度对致密区间厚度影响较小，而对范德华间隙厚度影响较大；并且在玻璃化转换温度之前随着温度升高范德华间隙厚度波动下降，玻璃化转换温度之后随温度升高范德华间隙厚度逐步增大。
　　关键词： 石墨烯改性沥青;分子动力学模拟;范德华间隙;致密区间;温度影响;
　　基金： 国家自然科学基金项目，项目编号11262002;
　　沥青主要是来自原油蒸馏的副产物，是一种高度复杂的材料，包含105～106种不同的分子[1]。近年来，常加入纳米材料对沥青进行改性，以获得性能更优异、更稳定的沥青混合料[2,3]。沥青纳米复合材料在双材料界面是薄弱环节，界面微结构特征对复合材料力学性能、破坏强度都有显著影响[4]。因此界面微结构特征是研究复合材料力学性能的重点。
　　文献[5,6,7]分别采用扫描电镜、X射线、原子力显微镜对沥青复合材料微观界面进行直接观察，以得到复合材料界面破坏机制。其试验费时费力，并且微观上不便于了解沥青与其他成分的作用机制。分子动力学(Molecular Dynamics, 简称MD)模拟是研究这些微观作用机制的有效途径之一。Xu等[8]利用MD方法对沥青的密度、溶解度参数等物理性质进行预测，并通过对黏附功的构成进行分析，得出范德华(van der Waals, 简称vdW)力对沥青结合料的黏结性能起着关键作用的结论。Ramezani等[9]由MD模拟得出，在碳纳米管改性沥青界面中，vdW力作用对黏结性能做出了主要贡献。因此，vdW力作用对复合材料界面力学性能的影响显著。
　　沥青作为一种热敏感性材料，环境温度的变化会使沥青料面层处于非定常和非均匀的状态[10]。研究变温度下沥青的力学性能是当前的热点。Zhou等[11]通过试验测量了碳纳米管改性沥青和石墨烯改性沥青(Graphene Modified Asphalt, 简称GMA)的玻璃化转换温度等热力学性能，后与MD模拟结果对比，计算得到的温度值和试验结果吻合良好；模拟发现随着温度升高，改性沥青的杨氏模量、剪切模量会逐渐减小。Liu等[12]对沥青混合料不同温度下静态动态荷载作用的力学性能进行分析，结果表明温度的升高会使沥青力学性能下降。Ouyang等[13]建立数学模型描述不同温度下沥青的抗压强度和杨氏模量；但该方法由于模型使用温度范围内不能存在玻璃化转换温度而具有明显的局限性。Zhu等[14]研究了二氧化硅对改性沥青的热力学和力学性能方面的影响，结果表明改性沥青的玻璃化转换温度随着二氧化硅掺入量的增加而降低，且力学性能随着填料的增加而逐渐提升；但对于温度对力学性能影响的微观机制，文献并没有进一步说明。
　　综上所述，沥青复合材料的界面微结构特征对研究复合材料力学性能具有重要意义，而温度对界面微结构特征的影响研究甚少。本文选取沥青的4组分模型，将石墨烯作为增强材料加入沥青，采用软件Material Studio 2017(简称MS)进行分子动力学模拟，首先计算得到沥青的玻璃化转换温度并与试验值对比，然后通过切片方法得到石墨烯/沥青复合材料相对浓度曲线，求出vdW间隙厚度，最后采用统计学方法求得致密区间厚度，并研究了在不同温度下vdW间隙和致密区间的变化情况。   1 分子动力学模拟理论
　　MD模拟的基本方法是将模拟物质看作一个集合体，之后通过数值求解该集合体体系经典力学运动方程，并统计体系的结构特征与性质[15]。通常，统计体系有以下几种，包括微正则系综(NVE系综)、正则系综(NVT系综)、等温-等压系综(NPT系综)、等焓-等压系综(NPH系综)和广义系综[16]。本文采用的NVT系综是指粒子数  N  、体积  V  、温度  T  都相同的热力学体系，晶胞内温度恒定，总能量发生变化。NPT系综内部体系压力恒定，但体积会发生改变，模拟中常用于压缩晶胞体积以达到真实材料的密度。
　　在GMA的动力学模拟中，采用凝聚相分子势能优化力场(Condensed-phase Optimized Molecular Potentials for Atomistic Simulation Studies, 简称COMPASS)。该力场用于描述键合相互作用和非键合相互作用，对于聚合物、无机小分子和大多数有机物，已被验证能够非常精确地拟合试验数据[17]。力场可表达为：
　　Ebonded  =  Eb  +  Eθ  +  Et  +  Eo  +  Eb  ,  b  ′+  Eb  ,  θ  +  Eθ  ,  θ  ′+  Eb  ,  t  +  Eθ  ,  t  +  Eθ  ,  t  ,  θ  ′ (1)
　　Enon  -  bonded  =  EvdW  +  Ecoul   (2)
　　Etotal  =  Ebonded  +  Enon  -  bonded   (3)
　　式中：  Eb  为键伸缩势；  Eθ  为键角弯曲势；  Et  为二面角扭曲势；  Eo  为离面弯曲势能；交叉耦合项包括伸缩-伸缩耦合势  Eb  ,  b  ′、伸缩-弯曲耦合势  Eb  ,  θ  、弯曲-弯曲耦合势  Eθ  ,  θ  ′、伸缩-扭转耦合势  Eb  ,  t  、弯曲-扭曲耦合势  Eθ  ,  t  、弯曲-扭曲-弯曲耦合势  Eθ  ,  t  ,  θ  ′;非键合相互作用包括vdW相互作用势  EvdW  、静电相互作用势  Ecoul  。
　　模型搭建使用Amorphous Cell(简称AC),基于MD的Forcite模块进行操作计算。vdW相互作用采用基于原子的求和法，静电相互作用采用Ewald求和法，截断半径取1.55 nm。除GMA模型弛豫外，其余过程均选用Nose恒温器和Berendsen恒压器分别用于维持目标温度和压力。所有的模拟都是在1 fs的时间步长、1个大气压的压力下进行的。   2 模型建立2.1沥青模型建立
　　美国战略公路研究计划(SHRP)根据溶解度将沥青分为了沥青质、极性芳香烃、环烷芳香烃、饱和烃等4种组分[18,19]。Derek等[20]在此基础上提出4组分模型。本文选取该模型是因为其能很好地模拟真实情况下的沥青，而且它们的分子量和分子面积基本相同，便于计算。模型每个单独的分子结构如图1所示，模型采用相关参数见表1。
　　图1 沥青4组分模型分子结构 下载原图
　　表1 沥青4组分模型参数 导出到EXCEL
　　组分
　　分子名称
　　化学式
　　分子量g/mol分子量g/mol
　　分子个数个分子个数个
　　沥青质
　　Asphaltene-phenol
　　C42H54O
　　574.89
　　3
　　Asphaltene-pyrrole
　　C66H81N
　　888.38
　　2
　　Asphaltene-thiophene
　　C51H62S
　　707.11
　　3
　　饱和烃
　　Squalane
　　C30H62
　　422.82
　　4
　　Hophane
　　C35H62
　　482.88
　　4
　　环烷
　　芳香烃
　　PHPN
　　C35H44
　　464.73
　　11
　　DOCHN
　　C30H46
　　406.69
　　13
　　极性
　　芳香烃
　　Quinolinohopane
　　C40H59N
　　553.91
　　4
　　Thioisorenieratane
　　C40H60S
　　572.97
　　4
　　Trimethylbenzeneoxane
　　C29H50
　　414.71
　　5
　　Pyridinohopane
　　C36H57N
　　503.85
　　4
　　Benzobisbenzothiophene
　　C18H10S2
　　290.39
　　15
　　在选取了沥青组分以后，需要建立沥青晶胞，具体步骤如下。
　　(1)在MS中建立好沥青的12个分子，然后分别对每个分子进行几何优化和能量最小化。
　　(2)构建沥青分子的AC模型。为防止分子链相互缠绕，必须使上述分子分布完全随机，并在AC模块以0.1 g/cm3的初始密度进行沥青三维周期模型的构建。建立好的沥青初始晶胞尺寸为8.05 nm×8.05 nm×8.05 nm。然后对该晶胞再进行一次几何优化，迭代次数取10 000次，精度取最高。
　　(3)沥青分子AC模型的退火(Annealing)操作。由于AC模块采用的是蒙特卡洛法，所有分子是完全随机分布的，所以由此搭成的模型能量极高，不利于后续计算，要对其进行退火操作。退火时，温度区间取300 K～800 K,精度取最高，进行50次迭代，并在迭代过程中再次进行分子优化。完成后在这50个构型中选取能量最低的一帧进行后续操作。
　　(4)经过退火操作后，为了让沥青体系在目标温度下从初始状态进入更为平衡的状态，实现体系的预平衡，模型进行500 ps的NVT系综弛豫，之后再进行500 ps的NPT系综弛豫。弛豫过程如图2所示。随着时间增加，沥青三维周期模型晶胞尺寸减小，密度逐步上升至真实密度。在温度  T  为298 K、1个大气压下(  p  =1 atm),弛豫后的沥青模型如图3所示。弛豫后尺寸为3.752 nm×3.752 nm×3.752 nm, 密度稳定在0.99g/cm3,与SHRP试验[19]所测得的1.03 g/cm3相差不大，结果理想。
　　图2 NPT弛豫下的密度变化 下载原图   2.2石墨烯模型的建立
　　在建立石墨烯模型时，取用MS中自带的周期性石墨烯三维模型。由于在后续操作中进行复合材料界面搭建模型时，石墨烯片晶胞尺寸需同沥青模型晶胞尺寸相近，才能够进行搭建，故需要在单个石墨烯的基础上建立超晶胞。石墨烯超晶胞平面尺寸为3.935 nm×3.834 nm, 如图4所示。
　　图3 沥青晶胞模型 下载原图
　　图4 石墨烯超晶胞模型 下载原图   2.3复合材料模型的建立
　　界面搭建时，为了使石墨烯在系综弛豫后仍在晶胞范围内，研究中放置两片石墨烯加以稳定。为防止底层石墨烯对沥青产生影响，导致结果错误，需对底层石墨烯坐标进行固定。搭建好界面模型后，需要对界面搭建模型进行几何优化、NVT系综弛豫，使能量预稳定。最后进行500 ps的NPT系综弛豫，恒温器选择Nose, 恒压器选择Souza-Martins, 将  y  方向压力设置为1个大气压(  p  =1 atm),其余方向全部设置为0。平衡后的晶胞如图5所示。   3 结果与分析3.1沥青玻璃化转换温度
　　沥青作为一种无定形物质，其分散体中存在大量胶体，致使存在一种特殊的温度，即玻璃化转换温度  Tg  。玻璃化温度  Tg  对沥青力学性能的影响很大。低于该温度之时沥青处于玻璃态，其本构关系为线性；而在高于该温度之时沥青处于高弹态，表现出黏弹性体的性质。对于沥青这样一种物理性质随温度变化的材料，需要对其玻璃化转换温度进行分析。
　　图5 石墨烯改性沥青模型 下载原图
　　在MS中编写脚本计算沥青模型的玻璃化转换温度  Tg  ,温度区间取在188.15 K～368.15 K之间，以20 K为温度间隔进行计算，模拟得到10个温度下对应的沥青密度，结果见表2。根据表2数据画出沥青温度～密度散点图，如图6所示。从图6中可以看出，沥青密度随温度的升高而降低，并且温度～密度曲线的斜率有明显变化的区间。玻璃化转换温度可以理解为温度～密度曲线上两条线性回归线相交点的温度。因此首先目测数据点，粗略估计出两条回归线相交位置。然后以该位置作为分界，将数据点分为左右两个部分，即分为188 K～248 K和268 K～368 K两个区间。分别用最小二乘法对两部分数据点进行回归分析，得到两条回归直线方程，即：
　　D  =-9.5×10-5  T  +1.047 (188≤  T  ≤248) (4)
　　D  =-4.7×10-4  T  +1.146 (268≤  T  ≤368) (5)
　　式中：  T  表示温度；  D  表示密度。
　　表2 不同温度下的沥青密度 导出到EXCEL
　　温度/K
　　188.15
　　208.15
　　228.15
　　248.15
　　268.15
　　288.15
　　308.15
　　328.15
　　348.15
　　368.15
　　密度/(g·cm-3)
　　1.028
　　1.028
　　1.027
　　1.022
　　1.018
　　1.011
　　1.002
　　0.993
　　0.984
　　0.970
　　图6 沥青温度～密度散点分布 下载原图
　　联立求解方程可以得到交点为(263.76,1.021 9),故沥青玻璃化转换温度  Tg  为263.76 K。
　　本文得到的沥青玻璃化温度与文献结果的比较列于表3。可以看出测量方法不同，各文献所得  Tg  差别较大。但由于  Tg  属于一个温度范围，故本文所得结果  Tg  =263.76 K是合理的，证明所采用的模型可以进行下一步的模拟计算。
　　表3 沥青玻璃化温度结果比较 导出到EXCEL
　　本文结果
　　文献结果
　　方法
　　Tg/K
　　方法
　　Tg/K
　　分子动力学
　　263.76
　　试验[21]
　　248.55
　　分子动力学[22]
　　298.15
　　试验[23]
　　223～303
　　分子动力学[11]
　　250
　　分子动力学[14]
　　278.66 3.2范德华间隙与致密区间
　　vdW间隙是石墨烯与沥青基体材料之间存在的由vdW力和静电力共同作用产生的真空带。随着作用距离的增加，vdW作用力会产生吸引势，作用于沥青基体，导致作用区域内密度增加。该作用区域即为致密区间。由于所用系统属于微观层面，故采用相对浓度对复合材料界面区间性质进行讨论更为精确。首先将晶胞等分为  N  个盒子，则相对浓度表达式为：
　　cR=cic¯ (0

去退保的人都是什么心理？为什么要退保？您好，我是读心解万事，是一名心理咨询师。对于很多人都把保险保费退了，这个事实是劝阻不了的，如果有也很少，尤其今年受疫情影响，生意不景气，人们收入降低，可能会有很多人退保费。其主要原华北平原地下有没有化石？地表以下多少米是坚硬岩石？肯定有化石存在的。首先说什么是化石。地壳中保存的历史久远的动物或植物的遗体遗物或生物留下的痕迹叫化石。华北平原自从地球形成开始就已经形成了，气候气温适宜，淡水资源充足，而且物产丰富原子弹之父是钱学森还是邓稼先？在基础科学领域，其实有很多科技成果是同宗同源的，但是每一位科学家在刻苦卓越的努力之上，各自侧重的研究方向又是不一样的原子弹之父钱三强过去以来，老一辈科学家钱三强的知名度好像不是很大人像摄影的照片不好看，你觉得是你拍摄的水平问题还是模特儿长得不好看？没镜头感的人，就要靠抓拍，最隹抓拍应该是在被拍者不注意，形态自然，这样的效果较好一些。人像拍摄，在摄影表现中，是一种常见的拍摄题材，对于摄影师来说，没有拍过人像的摄影师几乎没有。在在沈阳拍身份证照片能自拍吗？身份证是证明身份的重要凭证，如果它丢失了或者过期了，你不光跑各大银行要排队，去办理新的身份证也要排队好不容易排到了还没等你酝酿情绪闪光灯一个咔嚓于是就有了好的，下一位谁没被自己的身茄子饼怎么煎好吃又简单？不知道头条友们是不是喜欢吃茄子，美食坐家我是不怎么喜欢吃茄子的，因为常见的红烧茄子鱼香茄子做起来都要用很多的油，所以吃起来总是油腻腻的。可是前几天去朋友家吃了一次茄子饼以后我就彻底怎么样把手机图片转换成文字？这得看手机里的这张图是来至于何处？1来至于在网上下载的，可以在电脑上用360和UC浏览器打开你下载这张图的链接地址，然后这两个浏览器的设置里打开源码模式，很多都可以看到这张图的文字哪些女明星的穿搭简单又有气质？张瑜。孙丽穿搭非常时尚，漂亮，善良，有气质，我觉的还是一个非常好的演员。哪些女明星的穿搭简单又有气质？答，一般的女明星全都艰苦朴素，不人多显贵，明星穿搭什么衣服都有气质，因为她是明怎样在人头攒动的景点拍摄出没有人的空景照片？关于这个问题我给出四个解决方案，相信总有一个适合你。1拍摄时机这个很好理解，就是错开旅游的高峰期，同时最好选择景区一天中刚刚开始营业或者即将结束营业的时间拍摄。比如下面这张照片，我有哪些看着好厉害做法好简单的美食？我是爱烹饪的豆妈maomao，特别爱做饭，四年多来已分享过上千个菜谱，所以这个问题我来回答，绝对让你满意。我们以1日3餐为例，早餐2道，午餐3道，晚餐2道，一共分享7道颜值高，做法自己做面条，如何制作好吃又简单的卤子？夏天天气热，有时候就是馋这样一碗过水的凉面，浇上两勺喷香的汆卤，那叫一个好吃。北京人把各种拌面的料，称为汆，这种汆面，就是打卤面的意思。而这其中，茄丁汆面，是我最喜欢的。小时候，母
<<<<<<－>>>>>>