丁二酸氣溶膠吸濕生長后結(jié)構(gòu)和界面的分子動力學(xué)模擬

布龍祥　張超　凡鳳仙　楊楊　熊威

摘要：基于分子動力學(xué)方法模擬了丁二酸氣溶膠的吸濕生長現(xiàn)象，研究了溫度和水含量對顆粒平衡后結(jié)構(gòu)和界面特性的綜合影響機(jī)制。結(jié)果表明，溫度和水含量會影響丁二酸的溶劑化分子數(shù)，從而導(dǎo)致顆粒吸濕生長后呈現(xiàn)水包酸的分層結(jié)構(gòu)、混合結(jié)構(gòu)和酸包水的分層結(jié)構(gòu)。顆粒表面丁二酸含量隨溫度的升高先增加后減少，且隨水含量的增加而增加。此外，溫度高于 260 K 時，顆粒表面丁二酸分子的親水基和疏水基分別朝向顆粒內(nèi)部和氣相。

關(guān)鍵詞：氣溶膠；表面活性劑；吸濕生長；分子動力學(xué)

中圖分類號：? O 648???????????? 文獻(xiàn)標(biāo)志碼：?? A

Molecular dynamics simulation on the structure characteristics of succinic acid aerosols after hygroscopic growth

BU Longxiang，ZHANG Chao，F(xiàn)AN Fengxian，YANG Yang，XIONG Wei

（School of Energy and Power Engineering， University of Shanghai for Science and Technology， Shanghai 200093， China）

Abstract： The hygroscopic growth of aerosols has significant effects on their environmental behaviors and health effects. However， the knowledge on the microscopic kinetic mechanism for hygroscopic growth of aerosols containing surface active components has been largely lacking. In this study， the hygroscopic? growth? of? succinic? acid? aerosol? particle? was? simulated? using? the? molecular? dynamics method. The key point is to investigate the combined effects of temperature and water content on the structure and interface characteristics of particles at the equilibrium state. Results show that temperature and water content can affect the solvation number of succinic acid， resulting in the occurrence of the layered structure with water wrapping around acid particle， the mixed structure and the layered structure with? acid? wrapping? around? water? droplet.? The? distribution? of? succinic? acid? on? the? surface? region increases first and then decreases， and it increases with the increasing water content. Additionally， when the temperature is greater than 260 K， the hydrophilic and hydrophobic group of succinic acid molecules on the particle surface are toward the particle inside and the gas phase， respectively.

Keywords： aerosols; surfactant; hygroscopic growth; molecular dynamics

吸濕生長會使氣溶膠顆粒的物理化學(xué)性質(zhì)（包括粒徑、相態(tài)和化學(xué)組分等）發(fā)生改變，從而影響其環(huán)境效應(yīng)和健康效應(yīng)[1-5]，主要體現(xiàn)為：通過改變氣溶膠對光的吸收和散射特性影響輻射平衡和大氣能見度；通過改變云凝結(jié)核和冰核的數(shù)濃度分布從而影響降水過程；通過影響氣粒化學(xué)反應(yīng)從而改變大氣和顆粒的化學(xué)成分；通過改變氣溶膠動力學(xué)特性進(jìn)而影響其在人體呼吸系統(tǒng)內(nèi)部的沉積行為。表面活性有機(jī)物是實際大氣氣溶膠顆粒中一類常見且性質(zhì)特殊的物質(zhì)，通常包括有機(jī)酸（主要為二元羧酸、芳香酸、大分子脂肪酸）和脂肪酸鈉鹽等[6-7]，其來源主要有化石燃料燃燒、生物質(zhì)燃料燃燒、植物排放以及大氣中氣?；瘜W(xué)反應(yīng)。表面活性成分由于同時包含親水基和憎水基，趨向于分布在液滴的表面，在吸濕生長過程中其濃度會沿液滴徑向呈現(xiàn)梯度分布（分配動力學(xué)過程），進(jìn)而影響氣溶膠吸濕生長的穩(wěn)態(tài)特性和非穩(wěn)態(tài)過程。

當(dāng)前，對含表面活性成分氣溶膠顆粒吸濕特性的研究包括試驗研究和理論研究。試驗研究主要包括：基于串聯(lián)差分電遷移率分析儀、電動力學(xué)天平等分析相對濕度、粒徑、組分對水蒸氣次飽和條件下顆粒吸濕生長因子的影響機(jī)制[8-9]；基于云凝結(jié)核計數(shù)器對顆?；罨R界過飽和度的影響因素進(jìn)行分析[10]。理論研究方面，對含表面活性成分氣溶膠吸濕特性的描述經(jīng)歷了由簡單到詳細(xì)的過程，由最初在氣粒平衡理論（科勒理論）的基礎(chǔ)上只考慮表面活性成分對表面張力的影響，逐漸發(fā)展到同時考慮表面活性劑分子在顆粒中的分配動力學(xué)過程以及液滴溶液的非理想特性[11-13]。

雖然試驗和理論研究均取得了較大進(jìn)展，但受限于現(xiàn)有實驗儀器在測量微觀信息層面的缺陷以及理論模型中引入的諸多假設(shè)，當(dāng)前對于含表面活性成分氣溶膠顆粒吸濕特性機(jī)理的認(rèn)識仍存在不足。分子動力學(xué)模擬基于牛頓運(yùn)動定律對模擬體系中原子的運(yùn)動狀態(tài)進(jìn)行計算，在不引入其他假設(shè)條件的基礎(chǔ)上對研究對象的動力學(xué)過程進(jìn)行分析，因此，該方法在微觀機(jī)理描述方面有其獨特的優(yōu)勢。近年來，隨著計算機(jī)技術(shù)的發(fā)展，分子動力學(xué)模擬在氣溶膠吸濕生長微觀動力學(xué)機(jī)理研究領(lǐng)域的應(yīng)用日益增多。

Ma 等[14]建立了不同種類二元羧酸包裹水滴的分子動力學(xué)模型，觀察液滴在溫度為300 K 條件下達(dá)到平衡后的結(jié)構(gòu)和界面特性，發(fā)現(xiàn)二元羧酸和水的相分離狀態(tài)與分子的鏈長和水溶性有關(guān)。以大氣中冰核形成為背景， Vardanega 等[15]基于分子動力學(xué)模擬研究了甲酸?水二元顆粒的結(jié)構(gòu)特性，指出顆粒在不同溫度（100～250 K）和水含量（120～640 H2O）條件下會出現(xiàn)水包酸、混合和酸包水的結(jié)構(gòu)。Radola等[16]研究了乙酸和丙酸顆粒與水分子混合后的結(jié)構(gòu)（100～250 K ，120～720 H2O），發(fā)現(xiàn)顆粒出現(xiàn)了水包酸、部分溶解和酸包水的結(jié)構(gòu)，通過和先前研究的對比證明了羧酸顆粒作為凝結(jié)核的能力不僅與分子中氧碳比有關(guān)，還與羧酸分子中羧基數(shù)目以及分子結(jié)構(gòu)有關(guān)。 Zhang 等[17-18]采用分子動力學(xué)模擬方法，研究了鹵化鈉種類和濃度、溫度以及水含量對甲酸和乙酸顆粒吸濕生長后結(jié)構(gòu)和界面特性的影響機(jī)制，結(jié)果發(fā)現(xiàn)，初始混合顆粒被水溶解程度與溫度和濕度正相關(guān)，初始顆粒穩(wěn)定性和顆粒表面組分構(gòu)成與無機(jī)鹽種類和濃度密切相關(guān)。Karadima等[19]基于分子動力學(xué)模擬研究了相對濕度、表面活性劑種類和組分質(zhì)量比對多組分氣溶膠顆粒形態(tài)結(jié)構(gòu)的影響。

二元羧酸是大氣氣溶膠中一類常見的表面活性有機(jī)物，由于其較為簡單的分子構(gòu)型被廣泛應(yīng)用于氣溶膠吸濕特性機(jī)理研究中。然而，已有針對顆粒結(jié)構(gòu)特性的分子動力學(xué)研究中對于溫度和濕度條件綜合影響下二元羧酸顆粒吸濕生長后結(jié)構(gòu)和界面特性的認(rèn)識還很缺乏。本文選取丁二酸（C4H6O4）為代表性二元羧酸，建立了丁二酸顆粒吸濕生長的二元體系分子動力學(xué)模型，模擬研究不同溫度及水含量（對應(yīng)環(huán)境相對濕度）條件下顆粒結(jié)構(gòu)和界面特性的演化機(jī)制。

1 分子動力學(xué)模擬方法

模擬體系中丁二酸采用 OPLS-AA（optimized potentialsforliquidsimulationsinanall-atom representation）[20]力場參數(shù)，水分子選用 SPC/E 模型，該水分子模型的優(yōu)點是計算成本低，適用于水的相變過程研究。在模擬過程中采用 SHAKE 算法將水分子的 O?H 鍵長與 H?O?H 鍵角分別穩(wěn)定在0.1 nm 和109.47°。范德華作用力采用 Lennard- Jones 12-6勢能來計算。所有計算均基于開源程序 Large-scaleAtomic/MolecularMassivelyParallelSimulator（LAMMPS）[21]進(jìn)行。

模擬體系中分子間總勢能

式中：第一項為非鍵結(jié)項（n），包含 Lennard-Jones 12-6勢能和庫倫勢能；θ為鍵角；rij為2個原子的間距；σij為2個原子間勢能等于零時的間距；εij為勢能的最小值；ε0為真空介電常數(shù)； qi ，qj為2個原子的帶電量；在同一個分子中，1-4原子對應(yīng)系數(shù)fij =0.5，否則，系數(shù)fij =1。1-4原子是指通過共價鍵相互結(jié)合并被3個中間原子隔開的原子對。這種方法可以減少1-4原子對非鍵結(jié)作用的影響，這樣做的目的是為了更好地模擬分子動力學(xué)過程，并產(chǎn)生更準(zhǔn)確的計算結(jié)果[20]；后 3項分別為鍵伸縮項（b）、鍵角彎曲項（a）和二面角扭曲項（t）； Kr 和 Kθ為化學(xué)鍵鍵長和鍵角參數(shù)； r0和θ0分別是鍵長和鍵角的平衡值； V 為傅里葉級數(shù)系數(shù)；φ為二面角角度；庫倫勢能和 L- J 勢能的截斷半徑 r0=1.4 nm ，大于截斷半徑時，長程庫侖力則采用 PPPM（particle particleparticlemesh）算法計算。

模擬區(qū)域為立方體單元，x， y， z 方向均采用周期性邊界條件。在模擬開始時，將60個離散的丁二酸分子放置于模擬盒子（邊長為10 nm）內(nèi)，在 NVT 系綜下經(jīng)過2 ns后形成初始丁二酸干顆粒。隨后在干顆粒周圍放入不同數(shù)目的離散水分子（水分子個數(shù) NW=100， 500， 1000）以模擬濕度環(huán)境，此時，丁二酸的質(zhì)量分?jǐn)?shù)為80%，44%，28%。在溫度 T=100 K 的溫度環(huán)境下采用 NVT 系綜進(jìn)行模擬?？梢园l(fā)現(xiàn)在經(jīng)過2 ns后水分子均已經(jīng)完全吸附在有機(jī)物顆粒上。在此過程中，溫度設(shè)置為100 K 是為了讓水蒸氣分子在較短的時間內(nèi)凝結(jié)到顆粒表面。然后將混合顆粒轉(zhuǎn)移到新的模擬盒子中心，并且在盒子中添加氮氣和氧氣分子（比值為3∶1）模擬實際大氣組分，在不同的溫度下（T =160，260，280，300，320 K）采用 NPT 系綜分別對該體系進(jìn)行平衡模擬，壓強(qiáng) P=0.1 MPa。

在模擬中，分別使用 Nosé-Hoover 恒溫器和恒壓器來控制系統(tǒng)溫度和壓力，積分步長為1 fs。模擬過程中通過監(jiān)測丁二酸和水分子之間的能量變化判斷體系是否達(dá)到平衡狀態(tài)。圖1為水分子含量為500時不同溫度條件下模擬體系中丁二酸和水間勢能 E 的監(jiān)測結(jié)果， t 為時間。從圖1中可以看出，不同溫度條件下模擬體系達(dá)到平衡的時間有明顯差異。在勢能穩(wěn)定后，每經(jīng)過1 ps統(tǒng)計一次顆粒的相關(guān)信息，每種工況的統(tǒng)計時長為1 ns。

2 模擬結(jié)果分析

2.1 顆粒結(jié)構(gòu)圖

圖2為不同溫度和水含量條件下顆粒平衡結(jié)構(gòu)圖。從圖中可以看出，在 T=160 K 時，水蒸氣分子在丁二酸顆粒表面形成多個小團(tuán)簇，這種結(jié)構(gòu)被稱為分層結(jié)構(gòu)1（水分子以團(tuán)簇形式包裹酸顆粒）。隨著溫度升高到260 K，水分子會滲入初始的丁二酸顆粒內(nèi)部，但是，水分子的滲入程度隨著水分子含量的不同會存在較大的差異。當(dāng)水分子含量為100時，水分子分散分布在丁二酸顆粒中，整體呈現(xiàn)混合結(jié)構(gòu)。當(dāng)水分子含量為500和1000時，滲入顆粒內(nèi)部的水分子聚合形成中心液滴，丁二酸分布在顆粒表面。需要注意的是，液滴表面的丁二酸分子并非均勻分布，部分丁二酸分子聚集形成多個“島嶼”狀。這種顆粒結(jié)構(gòu)被稱為分層結(jié)構(gòu)2（酸包水）。丁二酸在顆粒表面的非均勻分布狀態(tài)和 Ma等[14]的模擬結(jié)果一致，主要原因在于丁二酸分子中包含2個羧基（親水基），這導(dǎo)致2個羧基中間的烷基（疏水基）和水滴的距離被迫減小，根據(jù)能量最小化原理，水滴表面的丁二酸會聚集在一起以減小總體和水的接觸面積。當(dāng)溫度繼續(xù)升高時（T=280～320 K），3種水含量條件下顆粒的形態(tài)并未進(jìn)一步發(fā)生明顯變化。

2.2 顆粒的組分分布

在上述定性描述的基礎(chǔ)上，對丁二酸?水混合顆粒在達(dá)到平衡后的結(jié)構(gòu)特性開展定量分析。沿徑向?qū)㈩w粒依次劃分成核心區(qū)、中間區(qū)和表面區(qū)，各區(qū)域體積分別占總體積的35%，35%和30%。圖3為不同工況下各區(qū)域中丁二酸質(zhì)量與顆粒中丁二酸總質(zhì)量的比值，即丁二酸質(zhì)量分?jǐn)?shù)wSA。計算過程中對于橫跨2個區(qū)域的丁二酸分子，將其拆分并將原子質(zhì)量納入所在區(qū)域丁二酸質(zhì)量的計算中。在較高溫度條件下，顆粒中液態(tài)水會出現(xiàn)蒸發(fā)現(xiàn)象，為了保證不同溫度條件下顆粒內(nèi)水分子數(shù)目一致，達(dá)到分析溫度對顆粒結(jié)構(gòu)影響的目的，首先統(tǒng)計獲得高溫條件下到達(dá)氣粒平衡狀態(tài)時體系中的水蒸氣分子數(shù) n ，然后將吸附了 NW 個水分子的丁二酸顆粒置于包含 n 個水蒸氣分子的環(huán)境中進(jìn)行模擬計算，根據(jù) K?hler理論[22]，溫度和粒徑恒定時，水蒸氣平衡蒸汽壓恒定，在環(huán)境中加入 n 個水蒸氣分子是為了抑制顆粒中水分子的蒸發(fā)，通過這種方法，可以將不同溫度條件下顆粒中包含的水分子數(shù)目差異控制在小于5的范圍內(nèi)。

如圖3所示，當(dāng)溫度 T 為160 K 時，3種不同水分子含量（100/500/1000 H2O）條件下顆粒中丁二酸主要分布在顆粒核心區(qū)，由于此時顆粒呈現(xiàn)水分子團(tuán)簇包裹丁二酸顆粒的情形，顆粒形狀的不規(guī)則程度較大，因此，不能從本圖分析水分含量對丁二酸在顆粒內(nèi)部分配的影響規(guī)律。隨著溫度的升高，顆粒核心區(qū)的丁二酸分子被滲入的水分子驅(qū)趕到中間區(qū)和表面區(qū)。在溫度為260 K 時，包含3種水分子數(shù)目顆粒的核心區(qū)丁二酸質(zhì)量分?jǐn)?shù)均降低，此后隨著溫度的繼續(xù)升高，顆粒核心區(qū)的丁二酸質(zhì)量分?jǐn)?shù)逐漸增多，中間區(qū)和表面區(qū)丁二酸的質(zhì)量分?jǐn)?shù)逐漸減少。由圖2可以看出，在160 K 時，顆粒呈現(xiàn)“水包酸”結(jié)構(gòu)，水以團(tuán)簇的形式包裹在丁二酸顆粒表面，這是由于低溫狀態(tài)下，水分子能量較低，水并未溶解丁二酸，丁二酸始終處于顆粒核心區(qū)，因此，在核心區(qū)占比較大。當(dāng)溫度升高時，即溫度為260 K 時，水分子會有更高的能量溶解丁二酸分子，且丁二酸分子中存在親水基與疏水基，開始分布于顆粒表面，因此，導(dǎo)致核心區(qū)丁二酸質(zhì)量急劇減小。而隨著溫度升高，丁二酸向核心區(qū)擴(kuò)散，導(dǎo)致丁二酸在核心區(qū)占比逐漸增大，表明溫度有增大丁二酸趨向顆粒內(nèi)部的能力。在溫度較高時（T=260 K），當(dāng)顆粒中水含量為100，500和1000時，顆粒核心區(qū)丁二酸質(zhì)量分?jǐn)?shù)分別為0.60，0.37和0.16，呈現(xiàn)減少趨勢，中間區(qū)丁二酸質(zhì)量分?jǐn)?shù)分別為0.38，0.58和0.77，呈現(xiàn)增長趨勢，表面區(qū)丁二酸質(zhì)量分?jǐn)?shù)分別為0.02，0.04和0.07，呈現(xiàn)增長趨勢。不同溫度條件（T=260～320 K）下水含量對不同區(qū)域內(nèi)丁二酸分布結(jié)果的影響一致。

為了進(jìn)一步對不同工況下顆粒的結(jié)構(gòu)和界面特性進(jìn)行定量研究，對 T=260，280，300和320 K，NW=100，500和1000時顆粒中各原子的徑向數(shù)密度分布函數(shù)進(jìn)行了計算分析。徑向數(shù)密度ρ（r）表示以顆粒的質(zhì)心為中心，選定一個原子為“目標(biāo)原子”，與顆粒質(zhì)心的距離由 r 至r+dr間的分子個數(shù)為dN，則徑向數(shù)密度即為dN與厚度為dr的薄層體積的比值，它的物理意義為每個徑向區(qū)域內(nèi)分子的數(shù)密度。其表達(dá)式為

ρ（r）=（2）

對于低溫工況（ T =160 K），由于顆粒非球形，并未對其進(jìn)行徑向分布函數(shù)的計算。圖4給出了羰基中的氧原子（O（O=C））、羧基中的氫原子（H（COOH））、烷基中碳原子（C（C—H））、烷基中的氫原子（H（C—H））、水分子中的氧原子（O（H2O））沿顆粒徑向的數(shù)密度分布。

從圖4可以看出，當(dāng)水分子含量為100且溫度為260 K 時，5種原子在不同半徑位置處均有分布，而且溫度的升高并未對原子的分布和相對位置產(chǎn)生明顯影響。當(dāng)水分子含量為500且溫度為260 K 時，丁二酸中4種原子的密度分布在顆粒表面（r ≈1.8 nm）出現(xiàn)波峰，在 r <1 nm 的區(qū)域接近0，然而水分子的徑向分布函數(shù)則表明顆粒中形成中心水滴。隨著溫度的升高，丁二酸徑向數(shù)密度在顆粒表面的波峰高度逐漸變低，在顆粒內(nèi)部逐漸增大。當(dāng)水分子含量由500增加到1000時，組分徑向分布函數(shù)的差異主要體現(xiàn)為在相同溫度條件下小半徑區(qū)域丁二酸徑向數(shù)密度值的減小。這表明顆粒中丁二酸分布于表面區(qū)域的趨勢隨著溫度的降低和水分含量的增多而增強(qiáng)。此外，通過比較羧基氫原子和烷基氫原子沿顆粒徑向的相對位置可以看出，對于不同工況下的顆粒表面，丁二酸分子的羧基（親水基）朝向顆粒內(nèi)部，烷基（疏水基）朝向氣相一側(cè)。

2.3 溶劑化分子數(shù)

為了定量解釋不同工況下丁二酸?水顆粒平衡后的結(jié)構(gòu)差異，計算了丁二酸分子周圍水分子的徑向分布函數(shù) g（r）和溶劑化分子數(shù)（NS（s）A（ol）-W ， SA 和W 分別表示丁二酸分子和水分子）。原子（或分子）的徑向分布函數(shù)定義為一個原子（或分子）周圍另一種原子（或分子）的區(qū)域數(shù)密度與平均密度的比值。計算過程中丁二酸用羧基碳（C（COOH））代表，水分子用其中的氧原子（O（H2O））代表。溶劑化分子數(shù)用于描述溶劑對溶解物分子的影響。在本文中表示丁二酸分子周圍水分子平均數(shù)目。

NS（s）A（ol）-W表達(dá)式為

NS（s）A（ol）-W = NCrs—ol O =2πρ l0rsolr2g （r）dr（3）

式中： g（r）為羧基碳原子（C（COOH））周圍氧原子（O（H2O））分子的徑向密度函數(shù)；NCrs—ol O 為羧基碳原子與周圍氧原子的平均數(shù)目；ρ為羧基碳（C（COOH））的數(shù)密度；rsol為溶劑化殼半徑。

圖5給出了包含500個水分子顆粒中羧基碳（C（COOH））周圍氧原子（O（H2O））的徑向分布函數(shù) g（r）。

從圖5可以看出，不同溫度條件下 g（r）在 r <0.5 nm 時均包含2個較為明顯的波峰，代表羧基碳（C（COOH））周圍氧原子（O（H2O））的2種相對平衡位置。隨著溫度的升高， g（r）中波峰的高度增加，表明丁二酸周圍水分子數(shù)增多。式（3）中rsol為 g（r）中第二個波谷出現(xiàn)的位置，對于不同的工況，rsol均在0.455～0.495 nm 之間。圖6為不同工況下丁二酸分子的溶劑化分子數(shù)NS（s）A（ol）-W。

從圖6中可以看出，隨著溫度的升高，丁二酸的溶劑化分子數(shù)增大，這表明丁二酸被水分子溶解的程度隨著溫度的升高而增大，這導(dǎo)致了圖3和圖4中顆粒中組分分布隨溫度的變化，即溫度從160 K 升高到260 K 時，丁二酸顆粒表面的水分子滲入內(nèi)部，而高溫條件（T ≥260 K）下，液滴表面丁二酸隨著溫度升高向顆粒中心區(qū)域移動。

在溫度較低時（T =160 K），NS（s）A（ol）-W 與顆粒中水分子含量沒有明顯的相關(guān)關(guān)系，主要是因為此時水分子以團(tuán)簇形式包裹在丁二酸顆粒表面，只有表面丁二酸分子和水分子接觸。在溫度較高（ T ≥260 K）條件下，水分子含量增加會導(dǎo)致 N（r0）C — O 增大，這是由于水分含量的增加會導(dǎo)致水滴表面積增大，呈現(xiàn)酸包水結(jié)構(gòu)的顆粒中和丁二酸接觸的水分子增多。

3 結(jié)論

采用分子動力學(xué)方法模擬研究了丁二酸氣溶膠顆粒吸濕生長后的平衡狀態(tài)，通過分析模擬快照、組分沿顆粒徑向的數(shù)密度分布，揭示了溫度和水含量對顆粒結(jié)構(gòu)和界面特性的綜合影響機(jī)制，基于溶劑化分子數(shù)對其中的機(jī)理進(jìn)行了分析，主要結(jié)論為：隨著溫度的升高，包含100個水分子的二元顆粒結(jié)構(gòu)由水包酸的分層結(jié)構(gòu)轉(zhuǎn)變成混合結(jié)構(gòu)，當(dāng)水含量增長為500和1000時，二元顆粒先后呈現(xiàn)水包酸的分層結(jié)構(gòu)、混合結(jié)構(gòu)和酸包水的分層結(jié)構(gòu)；丁二酸在二元顆粒核心區(qū)域的質(zhì)量分?jǐn)?shù)隨著溫度的升高先減少后增多，且隨著水含量的增加而減少，而在中間環(huán)狀區(qū)域和表面區(qū)域則呈現(xiàn)相反的變化趨勢；較高溫度條件下（T=260～320 K），二元顆粒表面丁二酸的親水基朝向顆粒內(nèi)部，疏水基朝向氣相一側(cè)；顆粒中丁二酸溶劑化分子數(shù)與溫度（T ≥260 K）和水分含量均為正相關(guān)關(guān)系。

參考文獻(xiàn)：

[1] ZHANG C， MA N， FAN F X， et al. Hygroscopic growth of aerosol particles consisted of oxalic acid and its internal mixture with ammonium sulfate for the relative humidity ranging from 80% to 99.5%[J]. Atmospheric Environment， 2021， 252：118318.

[2] CHEN X L， MA R， ZHONG W Q， et al. Numerical study of the effects of temperature and humidity on the transport anddepositionofhygroscopicaerosolsinaG3-G6 airway[J]. International Journal of Heat and Mass Transfer， 2019， 138：545–552.

[3] P?SCHLU.Atmosphericaerosols：composition， transformation， climate and health effects[J]. AngewandteChemie International Edition， 2005， 44（46）：7520–7540.

[4]郭陽， 凡鳳仙， 張超， 等.氨法脫硫系統(tǒng)排放細(xì)顆粒物的異質(zhì)核化特性[J].動力工程學(xué)報， 2022， 42（1）：49–55.

[5]鐘佳利， 王煒罡， 彭超， 等.大氣氣溶膠吸濕性及其對環(huán)境的影響[J].化學(xué)進(jìn)展， 2022， 34（4）：801–814.

[6] PETTERS S S， PETTERS M D. Surfactant effect on cloud condensationnucleifortwo-componentinternallymixed aerosols[J].JournalofGeophysicalResearch： Atmospheres， 2016， 121（4）：1878–1895.

[7] MCNEILL V F，SAREEN N，SCHWIER A N. Surface-activeorganicsinatmosphericaerosols[J]. Topicsin Current Chemistry， 2014， 339：201–259.

[8] BOUZIDIH，ZUENDA，ONDR??EKJ，etal. Hygroscopicbehaviorofinorganic –organicaerosol systems including ammoniumsulfate， dicarboxylic acids， andoligomer[J]. AtmosphericEnvironment， 2020， 229：117481.

[9] CHOCZYNSKI J M， KOHLI R K， SHELDON C S， et al. A dual-droplet approach for measuring the hygroscopicity ofaqueousaerosol[J].AtmosphericMeasurement Techniques， 2021， 14（7）：5001–5013.

[10] WANG Y， VOLIOTIS A， HU D W， et al. On the evolution ofsub-andsuper-saturatedwateruptakeofsecondary organicaerosolinchamberexperimentsfrommixed precursors[J]. Atmospheric Chemistry and Physics， 2022， 22（6）：4149–4166.

[11] SORJAMAA R， SVENNINGSSON B， RAATIKAINEN T， etal.TheroleofsurfactantsinK?hlertheory reconsidered[J]. AtmosphericChemistryandPhysics， 2004， 4（8）：2107–2117.

[12] PRISLE N L， RAATIKAINEN T， LAAKSONEN A， et al. Surfactantsinclouddropletactivation： mixedorganic- inorganic particles[J]. Atmospheric Chemistry and Physics， 2010， 10（12）：5663–5683.

[13] ZHANG C， BU L X， FAN F X， et al. Surfactant effect on the hygroscopicity of aerosol particles at relative humidity ranging from 80% to 99.5%： Internally mixed adipic acid- ammonium sulfate particles[J]. Atmospheric Environment， 2021， 266：118725.

[14] MA X F， CHAKRABORTY P， HENZ B J， et al. Molecular dynamicsimulationof dicarboxylicacidcoatedaqueous aerosol： structureandprocessingofwatervapor[J]. PhysicalChemistryChemicalPhysics， 2011，13（20）：9374–9384.

[15] VARDANEGAD，PICAUDS. Waterandformicacid aggregates： a molecular dynamics study[J]. The Journal of Chemical Physics， 2014， 141（10）：104701.

[16] RADOLAB，PICAUDS，VARDANEGAD，etal. Molecular dynamics simulations of the interaction between water molecules and aggregates of acetic or propionic acid molecules[J]. The Journal of Physical ChemistryB， 2015， 119（51）：15662–15674.

[17] ZHANG C， WANG Y S， WANG H J. Interaction between water and acetic acid-sodium halide aerosol： A molecular dynamics study[J]. Powder Technology， 2017， 314：9–19.

[18] ZHANG C， WANG Y S， WANG H J， et al. Microscopic mechanism of the interaction between water and formic acid-sodium chloride aerosol[J]. Powder Technology， 2019， 344：46–57.

[19] KARADIMA K S， MAVRANTZAS V G， PANDIS S N. Insights into the morphology of multicomponent organic and inorganic aerosols from molecular dynamics simulations[J]. Atmospheric Chemistry and Physics， 2019， 19（8）：5571–5587.

[20] JORGENSEN W L， MAXWELL D S， TIRADO-RIVES J. Development and testing of the OPLS all-atom force field on conformational energetics and properties of organic liquids[J]. Journal of the American Chemical Society， 1996， 118（45）：11225–11236.

[21] PLIMPTON S. Fast parallel algorithms for short-range molecular dynamics[J]. Journal of Computational Physics， 1995， 117（1）：1–19.

[22] K?HLER H. The nucleus in and the growth of hygroscopic droplets[J]. Transactions of the Faraday Society， 1936， 32（32）：1152–1161.

（編輯：石瑛）