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

　　黄立新 周昆 黄君 叶太智广西大学土木建筑工程学院 广西大学工程防灾与结构安全教育部重点实验室
　　摘 要： 采用分子动力学模拟方法研究了温度对石墨烯/沥青复合材料界面微结构特征的影响。选取沥青的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

添丁！沈阳森林动物园迎来出生潮从3月开始，园里的动物们开始陆续产崽，随着它们逐渐长大，游客们就可以欣赏新生小动物了。3月14日，记者从沈阳森林动物园看到，盘羊梅花鹿斑马等动物陆续产崽，现在已经与游客见面了。新生太可怕了！森林狼球星爱德华兹受伤惨叫，脚踝变形令人毛骨悚然！在对阵公牛的比赛中，森林狼球星爱德华兹在第一节因扭伤脚踝提前退出比赛。当时，爱德华兹跳投后落地时，脚踝严重变形，这位年仅21岁的全明星发出了可怕的尖叫声，他的惨叫在整个竞技场中都能亏损近6亿，全年销量仅19万，这家手机厂商危险了！国内手机市场竞争激烈，在过去几年逐渐饱和，排名经历了一轮又一轮的大洗牌。如今国内手机行业的格局有了新的变化，份额前四的品牌是VO荣米。最早的组合其实是中华酷联。中兴华为酷派和联想，老照片在朝鲜进行慰问的梦露，70年代的张学良越南战争时期，一名女子落入了美军的手中。这里是一片荒郊野地，可能这个女子是游击队员，美军拿着枪在对她进行审问。这个女子的脸上露出了恐惧表情，她应该如实交代了自己所知道的一切。照片中张颂文新片首映，王俊凯等众星到场支持，狂飙剧组来了4个人头条创作挑战赛狂飙火了，作为主演的张颂文出圈了，而张颂文也成为了作品保障，有他出演的作品会被观众认定为口碑好剧。3月20日，由张颂文白客苗苗等人领衔主演的电影不止不休在北京举行首映长江云调查千元机卖出万元价以租代买套路有多深？对于许多大学生来说，能用上一台高端手机有着巨大的吸引力。以租代买逐渐兴起，一些手机租赁平台打出15元天用高端手机的宣传，引导大学生超前消费。日前，长江云调查栏目推出了315线索征集老艺术家潘长江，不爱惜羽毛终将付出代价文小玖娱乐汇编辑小玖娱乐汇曾经的春晚舞台，含金量可谓是相当高，能在春晚上演绎作品的，几乎都是当代的知名艺术家。潘长江便是其中之一。作为东北丑角，潘长江在春晚上留下了十多部堪称经典的宝新能源签署广东陆丰甲湖湾电厂34号机组扩建工程银团贷款合同中证网讯（记者何昱璞）宝新能源3月20日晚间发布公告称，近日，为保障项目建设资金需求，公司全资子公司陆丰电力与多家银行共同签署了广东陆丰甲湖湾电厂34号机组扩建工程（21000MW首届中国商用车论坛丨远景动力将积极参与论坛分享商用车电动化转型经验十堰广电讯（全媒体记者张欣）首届中国商用车论坛将于3月28日至30日在十堰举行，远景动力技术（湖北）有限公司将积极参与论坛，分享相关经验。远景动力十堰零碳智能电池产业园项目总投资1孕妇待产包该怎么准备，才能避免临到医院生产时手足无措？记得老婆生一胎时，我们虽然提前在网上购买了一个孕妇待产包，但真正到医院生产时，发现还是解决不了问题。这里面有两个原因1）商家推荐的待产包并不完整，没有包含孕妇生产需要的所有东西2）大豪科技发布2022年业绩快报网安行业成第二增长极中证网讯（王珞）3月17日晚，大豪科技（603025。SH）发布2022年业绩快报，公告显示，2022年度公司实现营业收入15。98亿元，同比增长6。16实现归属于上市公司股东的净
<<<<<<－>>>>>>