丁二酸气溶胶吸湿生长后结构和界面的分子动力学模拟.pdf
《丁二酸气溶胶吸湿生长后结构和界面的分子动力学模拟.pdf》由会员分享,可在线阅读,更多相关《丁二酸气溶胶吸湿生长后结构和界面的分子动力学模拟.pdf(8页珍藏版)》请在文库网上搜索。
1、文章编号:1007 6735(2023)04 0345 07DOI:10.13255/ki.jusst.20220604003丁二酸气溶胶吸湿生长后结构和界面的分子动力学模拟布龙祥,张超,凡凤仙,杨杨,熊威(上海理工大学能源与动力工程学院,上海200093)摘要:基于分子动力学方法模拟了丁二酸气溶胶的吸湿生长现象,研究了温度和水含量对颗粒平衡后结构和界面特性的综合影响机制。结果表明,温度和水含量会影响丁二酸的溶剂化分子数,从而导致颗粒吸湿生长后呈现水包酸的分层结构、混合结构和酸包水的分层结构。颗粒表面丁二酸含量随温度的升高先增加后减少,且随水含量的增加而增加。此外,温度高于 260K 时,颗粒
2、表面丁二酸分子的亲水基和疏水基分别朝向颗粒内部和气相。关键词:气溶胶;表面活性剂;吸湿生长;分子动力学中图分类号:O648文献标志码:AMolecular dynamics simulation on the structure characteristics ofsuccinic acid aerosols after hygroscopic growthBU Longxiang,ZHANG Chao,FAN Fengxian,YANG Yang,XIONG Wei(School of Energy and Power Engineering,University of Shanghai f
3、or Science and Technology,Shanghai 200093,China)Abstract:Thehygroscopicgrowthofaerosolshassignificanteffectsontheirenvironmentalbehaviorsandhealtheffects.However,theknowledgeonthemicroscopickineticmechanismforhygroscopicgrowthofaerosolscontainingsurfaceactivecomponentshasbeenlargelylacking.Inthisstu
4、dy,thehygroscopic growth of succinic acid aerosol particle was simulated using the molecular dynamicsmethod.Thekeypointistoinvestigatethecombinedeffectsoftemperatureandwatercontentonthestructureandinterfacecharacteristicsofparticlesattheequilibriumstate.Resultsshowthattemperatureandwatercontentcanaf
5、fectthesolvationnumberofsuccinicacid,resultingintheoccurrenceofthelayeredstructurewithwaterwrappingaroundacidparticle,themixedstructureandthelayeredstructurewith acid wrapping around water droplet.The distribution of succinic acid on the surface regionincreasesfirstandthendecreases,anditincreaseswit
6、htheincreasingwatercontent.Additionally,whenthetemperatureisgreaterthan260K,thehydrophilicandhydrophobicgroupofsuccinicacidmoleculesontheparticlesurfacearetowardtheparticleinsideandthegasphase,respectively.Keywords:aerosols;surfactant;hygroscopic growth;molecular dynamics上 海 理 工 大 学 学 报第45卷第4期J.Univ
7、ersityofShanghaiforScienceandTechnologyVol.45No.42023收稿日期:20220604基金项目:国家自然科学基金资助项目(52106207)第一作者:布龙祥(1997),男,硕士研究生研究方向:气溶胶分子动力学模拟E-mail:通信作者:张超(1989),男,讲师研究方向:多相流动与传热E-mail:吸湿生长会使气溶胶颗粒的物理化学性质(包括粒径、相态和化学组分等)发生改变,从而影响其环境效应和健康效应1-5,主要体现为:通过改变气溶胶对光的吸收和散射特性影响辐射平衡和大气能见度;通过改变云凝结核和冰核的数浓度分布从而影响降水过程;通过影响气粒化学
8、反应从而改变大气和颗粒的化学成分;通过改变气溶胶动力学特性进而影响其在人体呼吸系统内部的沉积行为。表面活性有机物是实际大气气溶胶颗粒中一类常见且性质特殊的物质,通常包括有机酸(主要为二元羧酸、芳香酸、大分子脂肪酸)和脂肪酸钠盐等6-7,其来源主要有化石燃料燃烧、生物质燃料燃烧、植物排放以及大气中气粒化学反应。表面活性成分由于同时包含亲水基和憎水基,趋向于分布在液滴的表面,在吸湿生长过程中其浓度会沿液滴径向呈现梯度分布(分配动力学过程),进而影响气溶胶吸湿生长的稳态特性和非稳态过程。当前,对含表面活性成分气溶胶颗粒吸湿特性的研究包括试验研究和理论研究。试验研究主要包括:基于串联差分电迁移率分析仪
9、、电动力学天平等分析相对湿度、粒径、组分对水蒸气次饱和条件下颗粒吸湿生长因子的影响机制8-9;基于云凝结核计数器对颗粒活化临界过饱和度的影响因素进行分析10。理论研究方面,对含表面活性成分气溶胶吸湿特性的描述经历了由简单到详细的过程,由最初在气粒平衡理论(科勒理论)的基础上只考虑表面活性成分对表面张力的影响,逐渐发展到同时考虑表面活性剂分子在颗粒中的分配动力学过程以及液滴溶液的非理想特性11-13。虽然试验和理论研究均取得了较大进展,但受限于现有实验仪器在测量微观信息层面的缺陷以及理论模型中引入的诸多假设,当前对于含表面活性成分气溶胶颗粒吸湿特性机理的认识仍存在不足。分子动力学模拟基于牛顿运动
10、定律对模拟体系中原子的运动状态进行计算,在不引入其他假设条件的基础上对研究对象的动力学过程进行分析,因此,该方法在微观机理描述方面有其独特的优势。近年来,随着计算机技术的发展,分子动力学模拟在气溶胶吸湿生长微观动力学机理研究领域的应用日益增多。Ma 等14建立了不同种类二元羧酸包裹水滴的分子动力学模型,观察液滴在温度为 300K 条件下达到平衡后的结构和界面特性,发现二元羧酸和水的相分离状态与分子的链长和水溶性有关。以大气中冰核形成为背景,Vardanega 等15基于分子动力学模拟研究了甲酸水二元颗粒的结构特性,指出颗粒在不同温度(100250K)和水含量(120640H2O)条件下会出现水
11、包酸、混合和酸包水的结构。Radola 等16研究了乙酸和丙酸颗粒与水分子混合后的结构(100250K,120720H2O),发现颗粒出现了水包酸、部分溶解和酸包水的结构,通过和先前研究的对比证明了羧酸颗粒作为凝结核的能力不仅与分子中氧碳比有关,还与羧酸分子中羧基数目以及分子结构有关。Zhang 等17-18采用分子动力学模拟方法,研究了卤化钠种类和浓度、温度以及水含量对甲酸和乙酸颗粒吸湿生长后结构和界面特性的影响机制,结果发现,初始混合颗粒被水溶解程度与温度和湿度正相关,初始颗粒稳定性和颗粒表面组分构成与无机盐种类和浓度密切相关。Karadima 等19基于分子动力学模拟研究了相对湿度、表面
12、活性剂种类和组分质量比对多组分气溶胶颗粒形态结构的影响。二元羧酸是大气气溶胶中一类常见的表面活性有机物,由于其较为简单的分子构型被广泛应用于气溶胶吸湿特性机理研究中。然而,已有针对颗粒结构特性的分子动力学研究中对于温度和湿度条件综合影响下二元羧酸颗粒吸湿生长后结构和界面特性的认识还很缺乏。本文选取丁二酸(C4H6O4)为代表性二元羧酸,建立了丁二酸颗粒吸湿生长的二元体系分子动力学模型,模拟研究不同温度及水含量(对应环境相对湿度)条件下颗粒结构和界面特性的演化机制。1 分子动力学模拟方法模拟体系中丁二酸采用 OPLS-AA(optimizedpotentials for liquid simul
13、ations in an all-atomrepresentation)20力场参数,水分子选用 SPC/E 模型,该水分子模型的优点是计算成本低,适用于水的相变过程研究。在模拟过程中采用 SHAKE 算法将水分子的 OH 键长与 HOH 键角分别稳定在0.1nm 和109.47。范德华作用力采用Lennard-Jones12-6 势能来计算。所有计算均基于开源程序Large-scaleAtomic/MolecularMassivelyParallel346上海理工大学学报2023年第45卷Simulator(LAMMPS)21进行。模拟体系中分子间总势能Utot=n4ij12ijr12ij6
14、ijr6ij+140qiqjrijfij+bKr(rijr0)2+aK(0)2+tV121+cos+V221cos(2)+V321+cos(3)+V421cos(4)(1)式中:第一项为非键结项(n),包含 Lennard-Jones12-6 势能和库伦势能;为键角;rij为 2 个原子的间距;ij为 2 个原子间势能等于零时的间距;ij为势能的最小值;0为真空介电常数;qi,qj为 2 个原子的带电量;在同一个分子中,1-4 原子对应系数 fij=0.5,否则,系数 fij=1。1-4 原子是指通过共价键相互结合并被 3 个中间原子隔开的原子对。这种方法可以减少 1-4 原子对非键结作用的影
15、响,这样做的目的是为了更好地模拟分子动力学过程,并产生更准确的计算结果20;后3 项分别为键伸缩项(b)、键角弯曲项(a)和二面角扭曲项(t);Kr和 K为化学键键长和键角参数;r0和 0分别是键长和键角的平衡值;V 为傅里叶级数系数;为二面角角度;库伦势能和 L-J 势能的截断半径 r0=1.4nm,大于截断半径时,长程库仑力则采用 PPPM(particleparticleparticlemesh)算法计算。模拟区域为立方体单元,x,y,z 方向均采用周期性边界条件。在模拟开始时,将 60 个离散的丁二酸分子放置于模拟盒子(边长为 10nm)内,在NVT 系综下经过 2ns 后形成初始丁二
16、酸干颗粒。随后在干颗粒周围放入不同数目的离散水分子(水分子个数 NW=100,500,1000)以模拟湿度环境,此时,丁二酸的质量分数为 80%,44%,28%。在温度 T=100K 的温度环境下采用 NVT 系综进行模拟。可以发现在经过 2ns 后水分子均已经完全吸附在有机物颗粒上。在此过程中,温度设置为 100K是为了让水蒸气分子在较短的时间内凝结到颗粒表面。然后将混合颗粒转移到新的模拟盒子中心,并且在盒子中添加氮气和氧气分子(比值为31)模拟实际大气组分,在不同的温度下(T=160,260,280,300,320K)采用 NPT 系综分别对该体系进行平衡模拟,压强 P=0.1MPa。在模
17、拟中,分别使用 Nos-Hoover恒温器和恒压器来控制系统温度和压力,积分步长为 1fs。模拟过程中通过监测丁二酸和水分子之间的能量变化判断体系是否达到平衡状态。图 1 为水分子含量为 500 时不同温度条件下模拟体系中丁二酸和水间势能 E 的监测结果,t 为时间。从图 1 中可以看出,不同温度条件下模拟体系达到平衡的时间有明显差异。在势能稳定后,每经过 1ps 统计一次颗粒的相关信息,每种工况的统计时长为 1ns。02468102 4002 0001 6001 200800400200E/(kcalmol1)t/ns160 K 260 K 280 K 300 K320 K图 1 模拟体系中
18、丁二酸和水分子间的能量(NW=500)Fig.1 Interaction energy between succinic acid and watermolecules in the system(NW=500)2 模拟结果分析 2.1 颗粒结构图图 2 为不同温度和水含量条件下颗粒平衡结构图。从图中可以看出,在 T=160K 时,水蒸气分子在丁二酸颗粒表面形成多个小团簇,这种结构被称为分层结构 1(水分子以团簇形式包裹酸颗粒)。随着温度升高到 260K,水分子会渗入初始的丁二酸颗粒内部,但是,水分子的渗入程度随着水分子含量的不同会存在较大的差异。当水分子含量为 100 时,水分子分散分布在丁
19、二酸颗粒中,整体呈现混合结构。当水分子含量为500 和1000时,渗入颗粒内部的水分子聚合形成中心液滴,丁二酸分布在颗粒表面。需要注意的是,液滴表面的丁二酸分子并非均匀分布,部分丁二酸分子聚集形成多个“岛屿”状。这种颗粒结构被称为分层结构 2(酸包水)。丁二酸在颗粒表面的非均匀分布状态和 Ma 等14的模拟结果一致,主要原因在于丁二酸分子中包含 2 个羧基(亲水基),这导致 2 个羧基中间的烷基(疏水基)和水滴的距离被第4期布龙祥,等:丁二酸气溶胶吸湿生长后结构和界面的分子动力学模拟347迫减小,根据能量最小化原理,水滴表面的丁二酸会聚集在一起以减小总体和水的接触面积。当温度继续升高时(T=2
20、80320K),3 种水含量条件下颗粒的形态并未进一步发生明显变化。2.2 颗粒的组分分布在上述定性描述的基础上,对丁二酸水混合颗粒在达到平衡后的结构特性开展定量分析。沿径向将颗粒依次划分成核心区、中间区和表面区,各区域体积分别占总体积的 35%,35%和30%。图 3 为不同工况下各区域中丁二酸质量与颗粒中丁二酸总质量的比值,即丁二酸质量分数wSA。计算过程中对于横跨 2 个区域的丁二酸分子,将其拆分并将原子质量纳入所在区域丁二酸质量的计算中。在较高温度条件下,颗粒中液态水会出现蒸发现象,为了保证不同温度条件下颗粒内水分子数目一致,达到分析温度对颗粒结构影响的目的,首先统计获得高温条件下到达
21、气粒平衡状态时体系中的水蒸气分子数 n,然后将吸附了 NW个水分子的丁二酸颗粒置于包含 n 个水蒸气分子的环境中进行模拟计算,根据 Khler 理论22,温度和粒径恒定时,水蒸气平衡蒸汽压恒定,在环境中加入 n 个水蒸气分子是为了抑制颗粒中水分子的蒸发,通过这种方法,可以将不同温度条件下颗粒中包含的水分子数目差异控制在小于5 的范围内。如图 3 所示,当温度 T 为 160K 时,3 种不同水分子含量(100/500/1000H2O)条件下颗粒中丁二酸主要分布在颗粒核心区,由于此时颗粒呈现水分子团簇包裹丁二酸颗粒的情形,颗粒形状的不规则程度较大,因此,不能从本图分析水分含量对丁二酸在颗粒内部分
22、配的影响规律。随着温度的升高,颗粒核心区的丁二酸分子被渗入的水分子驱赶到中间区和表面区。在温度为 260K 时,包含 3 种水分子数目颗粒的核心区丁二酸质量分数均降低,此后随着温度的继续升高,颗粒核心区的丁二酸质量分数逐渐增多,中间区和表面区丁二酸的质量分数逐渐减少。由图 2 可以看出,在 160K 时,颗粒呈现“水包酸”结构,水以团簇的形式包裹在丁二酸颗粒表面,这是由于低温状态下,水分子能量较低,水并未溶解丁二酸,丁二酸始终处于颗粒核心区,因此,在核心区占比较大。当温度升高时,即温度为 260K 时,水分子会有更高的能量溶解丁二酸分子,且丁二酸分子中存在亲水基与疏水基,开始分布于颗粒表面,因
23、此,导致核心区丁二酸质量急剧减小。而随着温度升高,丁二酸向核心区扩散,导致丁二酸在核心区占比逐渐增大,表明温度有增大丁二酸趋向颗粒内部的能力。在温度较高时(T=260K),当颗粒中水含量为 100,500和 1000 时,颗粒核心区丁二酸质量分数分别为 0.60,0.37 和 0.16,呈现减少趋势,中间区丁二酸质量分数分别为 0.38,0.58 和 0.77,呈现增长趋势,表面区丁二酸质量160 K260 K280 K300 K320 K100 H2O500 H2O1 000 H2O图 2 不同温度及湿度下丁二酸和水混合颗粒的模拟快照(蓝色表示水分子,绿色表示丁二酸分子)Fig.2 Simu
24、lated snapshots of particles mixed with succinic acid and water at different temperatures and humidity(blue indicateswater molecules,green indicates succinic acid molecules)348上海理工大学学报2023年第45卷分数分别为 0.02,0.04 和 0.07,呈现增长趋势。不同温度条件(T=260320K)下水含量对不同区域内丁二酸分布结果的影响一致。为了进一步对不同工况下颗粒的结构和界面特性进行定量研究,对 T=260,2
25、80,300 和 320K,NW=100,500 和 1000 时颗粒中各原子的径向数密度分布函数进行了计算分析。径向数密度(r)表示以颗粒的质心为中心,选定一个原子为“目标原子”,与颗粒质心的距离由 r 至 r+dr 间的分子个数为 dN,则径向数密度即为 dN 与厚度为 dr 的薄层体积的比值,它的物理意义为每个径向区域内分子的数密度。其表达式为(r)=3dN4(r+dr)3r3(2)对于低温工况(T=160K),由于颗粒非球形,并未对其进行径向分布函数的计算。图 4 给出了羰基中的氧原子(O(O=C))、羧基中的氢原子(H(COOH))、烷基中碳原子(C(CH))、烷基中的氢原子(H(C
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 丁二酸 气溶胶 吸湿 生长 结构 界面 分子 动力学 模拟