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
26、H)）、水分子中的氧原子（O(H2O)）沿颗粒径向的数密度分布。从图 4 可以看出，当水分子含量为 100 且温度为 260K 时，5 种原子在不同半径位置处均有分布，而且温度的升高并未对原子的分布和相对位置产生明显影响。当水分子含量为 500 且温度为260K 时，丁二酸中 4 种原子的密度分布在颗粒表面（r1.8nm）出现波峰，在 r 1nm 的区域接近0，然而水分子的径向分布函数则表明颗粒中形成中心水滴。随着温度的升高，丁二酸径向数密度在颗粒表面的波峰高度逐渐变低，在颗粒内部逐渐增大。当水分子含量由 500 增加到 1000 时，组分径向分布函数的差异主要体现为在相同温度条件下小半径区域
27、丁二酸径向数密度值的减小。这表明颗粒中丁二酸分布于表面区域的趋势随着温度的降低和水分含量的增多而增强。此外，通过比较羧基氢原子和烷基氢原子沿颗粒径向的相对位置可以看出，对于不同工况下的颗粒表面，丁二酸分子的羧基（亲水基）朝向颗粒内部，烷基（疏水基）朝向气相一侧。2.3 溶剂化分子数NsolSA-W为了定量解释不同工况下丁二酸水颗粒平衡后的结构差异，计算了丁二酸分子周围水分子的径向分布函数 g(r)和溶剂化分子数（，SA 和W 分别表示丁二酸分子和水分子）。原子（或分子）的径向分布函数定义为一个原子（或分子）周围另一种原子（或分子）的区域数密度与平均密度的比值。计算过程中丁二酸用羧基碳（C(CO
28、OH)）代表，水分子用其中的氧原子（O(H2O)）代表。溶剂化分子数用于描述溶剂对溶解物分子的影响。在160260280300320020406080100wSA/%表面区中间区核心区T/K(a)Nw=100 H2O(b)Nw=500 H2O(c)Nw=1 000 H2O160260280300320020406080100wSA/%T/K表面区中间区核心区160260280300320020406080100wSA/%T/K表面区中部区核心区图 3 不同温度下颗粒中丁二酸在各区域（核心区、中间区和表面区）的质量分数Fig.3 Mass fraction of succinic acid in
29、 each region (core,intermediate ring and surface regions)at differenttemperatures第4期布龙祥，等：丁二酸气溶胶吸湿生长后结构和界面的分子动力学模拟349NsolSA-W本文中表示丁二酸分子周围水分子平均数目。表达式为NsolSA-W=12NrsolCO=2wrsol0r2g(r)dr（3）NrsolCO式中：g(r)为羧基碳原子（C（COOH）周围氧原子（O（H2O）分子的径向密度函数；为羧基碳原 子 与 周 围 氧 原 子 的 平 均 数 目；为 羧 基 碳（C（COOH）的数密度；rsol为溶剂化壳半径。图
30、5 给出了包含 500 个水分子颗粒中羧基碳（C（COOH）周围氧原子（O（H2O）的径向分布函数 g(r)。NsolSA-W从图 5 可以看出，不同温度条件下 g(r)在 r0.5nm 时均包含 2 个较为明显的波峰，代表羧基碳（C(COOH)）周围氧原子（O(H2O)）的 2 种相对平衡位置。随着温度的升高，g(r)中波峰的高度增加，表明丁二酸周围水分子数增多。式（3）中rsol为 g(r)中第二个波谷出现的位置，对于不同的工况，rsol均在 0.4550.495nm 之间。图 6 为不同工况下丁二酸分子的溶剂化分子数。从图 6 中可以看出，随着温度的升高，丁二酸的溶剂化分子数增大，这表明
31、丁二酸被水分子溶解的程度随着温度的升高而增大，这导致了图 3和图 4 中颗粒中组分分布随温度的变化，即温度从 160K升高到 260K 时，丁二酸颗粒表面的水分00.51.0r/nm100 H2O1.52.0120r/nm500 H2O1230r/nm1 000 H2O10203001020(r)/nm3(r)/nm3(r)/nm3300102030010203040T=320 KO(C=O)H(HO)C(CH)H(CH)O(H2O)T=300 KT=280 KT=260 KT=320 KT=280 KT=260 KT=320 KT=300 KT=280 KT=260 KT=300 K图 4
32、不同温度和水含量条件下丁二酸和水沿颗粒径向的数密度分布Fig.4 Number density distribution of succinic acid and water along the radial direction of the particles at different temperaturesand water contents0.20.40.60.81.0050100150200250300350400g(r)r/nm 160 K 260 K 280 K 300 K 320 K0.455 nm0.495 nm图 5 500 个 水 分 子 颗 粒 中 羧 基 碳（C(CO
33、OH)）周 围 氧 原 子（O(H2O)）的径向分布函数Fig.5 Radial distribution function of oxygen atoms(O(H2O)around carboxyl carbon (C(COOH)in particlescontaining 500 water molecules1602602803003200.51.01.52.02.5T/KNsolSA-W100 H2O500 H2O1 000 H2O340图 6 不同温度及水分子含量下，丁二酸分子的溶剂化分子数Fig.6 Solvation numbers of succinic acid molecu
34、le atdifferent temperatures and water molecule contents350上海理工大学学报2023年第45卷子渗入内部，而高温条件（T260K）下，液滴表面丁二酸随着温度升高向颗粒中心区域移动。NsolSA-W在温度较低时（T=160K），与颗粒中水分子含量没有明显的相关关系，主要是因为此时水分子以团簇形式包裹在丁二酸颗粒表面，只有表面丁二酸分子和水分子接触。在温度较高（T260K）条件下，水分子含量增加会导致N(r0)CO增大，这是由于水分含量的增加会导致水滴表面积增大，呈现酸包水结构的颗粒中和丁二酸接触的水分子增多。3 结论采用分子动力学方法模拟研
35、究了丁二酸气溶胶颗粒吸湿生长后的平衡状态，通过分析模拟快照、组分沿颗粒径向的数密度分布，揭示了温度和水含量对颗粒结构和界面特性的综合影响机制，基于溶剂化分子数对其中的机理进行了分析，主要结论为：随着温度的升高，包含 100 个水分子的二元颗粒结构由水包酸的分层结构转变成混合结构，当水含量增长为 500 和 1000 时，二元颗粒先后呈现水包酸的分层结构、混合结构和酸包水的分层结构；丁二酸在二元颗粒核心区域的质量分数随着温度的升高先减少后增多，且随着水含量的增加而减少，而在中间环状区域和表面区域则呈现相反的变化趋势；较高温度条件下（T=260320K），二元颗粒表面丁二酸的亲水基朝向颗粒内部，疏
36、水基朝向气相一侧；颗粒中丁二酸溶剂化分子数与温度（T260K）和水分含量均为正相关关系。参考文献：ZHANGC,MAN,FANFX,etal.Hygroscopicgrowthofaerosolparticlesconsistedofoxalicacidanditsinternalmixturewithammoniumsulfatefortherelativehumidityrangingfrom80%to99.5%J.AtmosphericEnvironment,2021,252:118318.1CHENXL,MAR,ZHONGWQ,etal.Numericalstudyoftheeffec
37、tsoftemperatureandhumidityonthetransportand deposition of hygroscopic aerosols in a G3-G6airwayJ.InternationalJournalofHeatandMassTransfer,2019,138:545552.2PSCHLU.Atmosphericaerosols:composition,transformation,climateandhealtheffectsJ.AngewandteChemieInternationalEdition,2005,44(46):75207540.3郭阳,凡凤仙
38、,张超,等.氨法脱硫系统排放细颗粒物的4异质核化特性 J.动力工程学报,2022,42(1):4955.钟佳利,王炜罡,彭超,等.大气气溶胶吸湿性及其对环境的影响 J.化学进展,2022,34(4):801814.5PETTERSSS,PETTERSMD.Surfactanteffectoncloudcondensation nuclei for two-component internally mixedaerosolsJ.JournalofGeophysicalResearch:Atmospheres,2016,121(4):18781895.6MCNEILLVF,SAREENN,SCHW
39、IERAN.Surface-active organics in atmospheric aerosolsJ.Topics inCurrentChemistry,2014,339:201259.7BOUZIDI H,ZUEND A,ONDREK J,et al.Hygroscopic behavior of inorganic organic aerosolsystems including ammonium sulfate,dicarboxylic acids,and oligomerJ.Atmospheric Environment,2020,229:117481.8CHOCZYNSKIJ
40、M,KOHLIRK,SHELDONCS,etal.Adual-dropletapproachformeasuringthehygroscopicityofaqueousaerosolJ.AtmosphericMeasurementTechniques,2021,14(7):50015013.9WANGY,VOLIOTISA,HUDW,etal.Ontheevolutionof sub-and super-saturated water uptake of secondaryorganic aerosol in chamber experiments from mixedprecursorsJ.
41、AtmosphericChemistryandPhysics,2022,22(6):41494166.10SORJAMAAR,SVENNINGSSONB,RAATIKAINENT,et al.The role of surfactants in Khler theoryreconsideredJ.Atmospheric Chemistry and Physics,2004,4(8):21072117.11PRISLENL,RAATIKAINENT,LAAKSONENA,etal.Surfactants in cloud droplet activation:mixed organic-inor
42、ganicparticlesJ.AtmosphericChemistryandPhysics,2010,10(12):56635683.12ZHANGC,BULX,FANFX,etal.Surfactanteffectonthehygroscopicityofaerosolparticlesatrelativehumidityrangingfrom80%to99.5%:Internallymixedadipicacid-ammoniumsulfateparticlesJ.AtmosphericEnvironment,2021,266:118725.13MAXF,CHAKRABORTYP,HEN
43、ZBJ,etal.Moleculardynamic simulation of dicarboxylic acid coated aqueousaerosol:structure and processing of water vaporJ.Physical Chemistry Chemical Physics,2011,13(20):93749384.14VARDANEGA D,PICAUD S.Water and formic acidaggregates:amoleculardynamicsstudyJ.TheJournalofChemicalPhysics,2014,141(10):1
44、04701.15RADOLA B,PICAUD S,VARDANEGA D,et al.MoleculardynamicssimulationsoftheinteractionbetweenwatermoleculesandaggregatesofaceticorpropionicacidmoleculesJ.TheJournalofPhysicalChemistryB,2015,119(51):1566215674.16 下转第 363 页 第4期布龙祥，等：丁二酸气溶胶吸湿生长后结构和界面的分子动力学模拟351BAEK S,KIM J H,JEONG S,et al.Development
45、 ofhighly effective cryogenic printed circuit heat exchanger(PCHE)with low axial conductionJ.Cryogenics,2012,52(7/9):366374.17张海燕,郭江峰,淮秀兰,等.PCHE 内轴向导热对局部换热 性 能 的 影 响 研 究 J.化 工 学 报,2019,70(12):45904598.18王绩德,冯岩,韩东江.超临界二氧化碳动力循环在钠冷快堆中的应用综述 J.核科学与工程,2019,39(2):289297.19SHIHY,LIMJ,WANGWQ,etal.Heattransfe
46、randfrictionofmoltensaltandsupercriticalCO2flowinginanairfoil channel of a printed circuit heat exchangerJ.International Journal of Heat and Mass Transfer,2020,150:119006.20何雅玲,王文奇,邱羽,等.熔盐在复杂换热结构内的对流换热特性实验研究及进展 J.科学通报,2019,64(28/29):30073019.21KATOY,NITAWAKIT,MUTOY.Mediumtemperature22carbondioxidega
47、sturbinereactorJ.NuclearEngineeringandDesign,2004,230(1/3):195207.高毅超,夏文凯,龙颖,等.管径和转折角对 Z 型PCHE 换热及压降影响的研究 J.热能动力工程,2019,34(2):94100.23WANG B,SHEN J F,CHENG J D,et al.Numericalimprovementusingflowandheattransfercalculationsofthe zigzag geometry for carbon dioxide PCHEsJ.Energies,2022,15(8):2831.24顾维藻
48、,神家锐,马重芳,等.强化传热 M.北京:科学出版社,1990:10.25史以奇,李凌.螺纹管管内流动与传热的数值模拟 J.上海理工大学学报,2016,38(2):133136,141.26KIM I H,NO H C,LEE J I,et al.Thermal hydraulicperformanceanalysisoftheprintedcircuitheatexchangerusing a helium test facility and CFD simulationsJ.NuclearEngineeringandDesign,2009,239(11):23992408.27（编辑：石瑛
49、）上接第 351 页 ZHANGC,WANGYS,WANGHJ.Interactionbetweenwaterandaceticacid-sodiumhalideaerosol:AmoleculardynamicsstudyJ.PowderTechnology,2017,314:919.17ZHANGC,WANGYS,WANGHJ,etal.Microscopicmechanism of the interaction between water and formicacid-sodium chloride aerosolJ.Powder Technology,2019,344:4657.18
50、KARADIMAKS,MAVRANTZASVG,PANDISSN.Insights into the morphology of multicomponent organicandinorganicaerosolsfrommoleculardynamicssimulationsJ.AtmosphericChemistryandPhysics,2019,19(8):55715587.19JORGENSENWL,MAXWELLDS,TIRADO-RIVESJ.DevelopmentandtestingoftheOPLSall-atomforcefieldon conformational ener
