扩展至混合物的巨等容绝热(Grand-Isochoric Adiabatic, μV L)蒙特卡洛模拟：氩-氪(Ar-Kr)二元混合物在体相及MCM-41介孔中的吸附

《The Journal of Physical Chemistry C》：Rapid Exploration of Mixture Adsorption via Adiabatic Sampling

　　

摘要：研究人员将先前发展的单组分巨等容绝热(grand-isochoric adiabatic, μVL)蒙特卡洛(Monte Carlo, MC)系综推广至多组分体系，推导了二元氩(Ar)–氪(Kr)混合物在体相及MCM-41圆柱形介孔中吸附的MC接受准则，并与传统等温巨正则(grand-canonical, μVT)模拟结果进行对比验证。通过ab initio(+三体Axilrod–Teller–Muto)势能模型描述Ar–Kr相互作用，研究人员发现两种系综所得温度、组成、势能及粒子数均吻合；靠近共存条件时巨正则模拟因需越过成核自由能垒而收敛显著变慢，而巨等容绝热模拟因允许瞬时升温可快速填充，收敛速度快1–2个数量级且不受过饱和度影响；该系综亦可有效捕捉吸附等otherm及选择性，并为多孔材料高通量筛选提供高效途径。

论文解读：扩展至混合物的巨等容绝热(μVL)蒙特卡洛模拟研究Ar–Kr混合物吸附

一、研究背景与意义

传统等温巨正则(grand-canonical, μVT)系综模拟靠近气液共存或介孔吸附台阶处时，系统须克服较大自由能势垒方能从低密度相转变为高密度相，导致收敛严重迟滞。巨等容绝热(grand-isochoric adiabatic, 亦称μVL)系综固定化学势μ_i、体积V及Hill能量L（而非温度T），粒子数N_i与动能可涨落，此前仅应用于单组分系统。本文将其推广至双组分Ar–Kr混合物，推导相应MC移动接受概率，验证其热力学一致性，并证明其在体相及限域吸附中可规避势垒诱导的收敛延迟，具重要方法学价值。论文发表于《The Journal of Physical Chemistry C》。

二、主要关键技术方法

研究人员采用Deiters等开发的ab initio两体(2B)势能加Axilrod–Teller–Muto三体(3B)项描述Ar–Ar、Kr–Kr及Ar–Kr相互作用，无组合规则；MCM-41圆柱孔用Steele型固体–流体(solid–fluid, sf)积分公式描述；通过微正则(microcanonical)相空间体积Ω(N₁,N₂,V,E)的Dirichlet积分导出多组分理想气体因子及密度态ω(μ₁,μ₂,V,L)；据此写出平移、插入(type i)、删除(type i)三类MC移动的接受比(acceptance ratio)，其中K=L+μ₁N₁+μ₂N₂?U为可用动能；先运行μVT模拟获取?L?=?U?+(3/2)(?N_Ar?+?N_Kr?)k_BT?μ_Ar?N_Ar??μ_Kr?N_Kr?确定对应L再开展μVL模拟，反向由μVL模拟算?T?；体相立方盒边长18 ?，MCM-41孔半径R=17 ?、长L_z=40 ?(部分长孔L_z=100 ?)，周期边界，截断加长程校正，各态平均基于2×10⁷MC步。

三、研究结果

Extension to Mixtures：通过对双组分微正则相空间动量积分应用Dirichlet公式，获得混合物的相空间体积Ω(N₁,N₂,V,E)及密度态ω(μ₁,μ₂,V,L)，单组分μVL系综的理想气体前置因子及动能项K=L+μN?U分别被双组分形式(b₁V)^N₁(b₂V)^N₂/(N₁!N₂!h^3(N₁+N₂))Γ[3(N₁+N₂)/2+1]与K=L+μ₁N₁+μ₂N₂?U取代。

Acceptance Rules for Grand-Isochoric Adiabatic Simulations of Mixtures：推导出三类MC移动接受准则——粒子平移(不改变N₁,N₂，仅含K_n^{3(N₁+N₂)/2?1}/K_o^{3(N₁+N₂)/2?1}之比)、组分i插入(含b_iVΓ(3(N₁+N₂)/2)/((N_i+1)Γ(3(N₁+N₂+1)/2))及动能幂次比)、组分i删除(含N_iΓ(3(N₁+N₂)/2)/(b_iVΓ(3(N₁+N₂?1)/2))及动能幂次比)。

Grand-Canonical and Grand-Isochoric Adiabatic Simulations of Mixtures：在等同(μ_Ar,μ_Kr,V)下运行μVT与μVL模拟，通过Eq 22与Eq 23互算?L?与?T?。对等摩尔液态Ar–Kr混合物(P=40–125 bar, T=143.15 K)，两系综所得摩尔Hill能L?、内能、温度、Ar摩尔分数x_Ar、势能U及总粒子数N一致，验证μVL系综正确性。

Extracting Hill Energy from Available Data：由实验或模拟给出的P,T,S?,V?按L?=T S??P V?计算摩尔Hill能，或由μVT模拟的Ξ(μ_Ar,μ_Kr,V,T)、?N_i?、内部能E算出同样吻合，单组分Ar与混合物端点均与文献一致。

Results and Discussion部分图示分析：纯Ar与Ar–Kr混合物的气液共存线在T–ρ与T–L?平面形状相似但分支翻转；混合物的P–x_Kr呈封闭环，P–L?呈开放两支(左低L?液相、右高L?气相)，端点分别对应纯Kr与纯Ar的L?值，证明ab initio+3B势能可正确描述Hill能变化。

Convergence Comparison (Bulk & Adsorption)：体相充填接近共存时(过饱和度s=P/P*从7降至6.4)，μVT模拟特征步数n_τ急剧增大(成核势垒效应)，μVL模拟n_τ基本恒定且较μVT快约10–100倍；原因是μVL允许初始温度大幅跃升(图4右下)，促进粒子随机插入被接受，随后冷却至平衡T。MCM-41内等摩尔Ar–Kr吸附(T=143.15 K, P=1–10 bar)两系综所得吸附等温线、各组分吸附数、势能一致；μVL较μVT收敛快1–2个数量级，且在吸附台阶区(P≈6 bar)仍快速达平衡。更长孔(L_z=100 ?)更低T(120 K)出现N与L双稳态(metastability)，μVT可滞留于亚稳低L值，μVL因直接指定平衡L可跳出亚稳直接达真平衡，再次佐证效率优势。

四、讨论与结论翻译

研究人员成功将巨等容绝热(grand-isochoric adiabatic, μVL)蒙特卡洛模拟方法推广至双组分混合物，严格导出了平移、插入与删除移动的Metropolis接受准则。通过与等温巨正则(grand-canonical, μVT)模拟对体相Ar–Kr液体及MCM-41限域吸附的详尽比对，证实两系综所得热力学性质（温度、组成、势能、吸附量、Hill能）高度一致，验证了该方法可靠性。关键优势在于：当热力学条件趋近共存或吸附台阶（低过饱和度）时，μVT模拟因需跨越自由能垒致收敛严重迟滞，而μVL模拟允许系统瞬时升温以高效采样高能组态，规避势垒，收敛速率基本不受过饱和度影响且较μVT快一至两个数量级。该μVL系综特别适用于存在亚稳态或需长时间 equilibration 的气–液共存及多孔材料内混合物吸附过程之模拟，可为大规模候选多孔吸附/分离材料的高通量计算筛选提供显著加速手段。