分子动力学扩散模拟：纳米通道中分子传输的量化分析路径

分子动力学扩散模拟在膜分离和纳米流体领域是核心手段，但受限空间中扩散行为偏离自由扩散模型、统计采样充分性难以保证、力场对扩散系数影响显著，如何输出可靠扩散数据指导膜材料设计，是计算纳米科学领域的关键技术挑战。

项目背景是一个碳纳米管(CNT)膜的气体分离性能评估任务。客户制备了垂直排列的CNT膜，实验测得H₂/CH₄选择性约15，需要用分子动力学扩散模拟量化解释CNT孔径对气体扩散速率和选择性的影响。分子动力学扩散模拟在纳米受限空间中是研究传输行为的直接手段，但受限扩散与自由扩散的模型差异是核心技术难点。

模型构建是第一步。选了(8,8)碳纳米管（直径10.85Å）和(12,12)碳纳米管（直径16.27Å），分别代表小孔和大孔两种情况。管长5nm，两端各接一个气体储库。储库用碳壁围成，体积约2×10⁻²⁰ m³，初始填充50个H₂和50个CH₄分子（混合比1:1）。分子动力学扩散模拟中储库设计的关键是保证浓度梯度稳定——管两端储库的浓度差驱动气体通过CNT，但模拟中浓度会随气体传输而变化，需要足够大的储库使浓度波动<5%。

力场选择和气体模型设置是技术核心。CNT用AIREBO力场（碳纳米管专用势函数），H₂用 Buch势参数（量子修正的Lennard-Jones，考虑零点能），CH₄用TraPPE势（联合原子模型，CH₄视为一个粒子）。交叉交互参数用Lorentz-Berthelot混合规则：σ_H₂-CNT = (σ_H₂+σ_CNT)/2 = 3.13Å，ε_H₂-CNT = sqrt(ε_H₂·ε_CNT) = 0.28 meV。分子动力学扩散模拟中交叉交互参数对受限扩散的影响极大——气体-管壁交互太强会吸附在壁面不传输，太弱则自由穿行无选择性。

扩散系数提取用三种方法交叉验证。方法一是MSD（均方根位移）：D = MSD(t)/(2nt)。在CNT中气体沿轴向传输，n=1（一维扩散）。方法二是均方根速度关联（Green-Kubo）：D = ∫⟨v(0)·v(t)⟩dt。方法三是粒子计数法：从穿过CNT的粒子数N随时间的变化率推算传输通量J，再由J=-D·∇c推算D。分子动力学扩散模拟中三种方法的一致性是数据可信度的核心指标。

(8,8) CNT中的扩散系数计算结果：H₂ D=1.8×10⁻⁵ cm²/s，CH₄ D=1.2×10⁻⁸ cm²/s。H₂/CH₄扩散选择性=1500。但实验选择性只有15——差距巨大的原因在于选择性不是由扩散系数比决定的，而是由通量比决定。分子动力学扩散模拟中通量=扩散×浓度梯度×面积。CH₄虽然扩散极慢但在CNT中的吸附密度远高于H₂——CH₄在管壁上的吸附使管内CH₄浓度比H₂高8倍。通量选择性=扩散比1500/浓度比8≈188，仍然偏高。

偏差排查揭示了实验和模拟的条件差异。模拟是纯CNT理想通道，实际CNT膜有管间空隙（inter-tube voids），这些空隙的直径约3-5Å，允许H₂快速通过但阻挡CH₄。实验测量的选择性是CNT通道和空隙通道的混合贡献。分子动力学扩散模拟中如果只计算CNT内的传输，忽略了空隙贡献，选择性会系统性偏高。补建了管间空隙模型（5Å缝隙），H₂在空隙中D=2.5×10⁻⁵ cm²/s（比CNT内更快），CH₄ D=4×10⁻¹² cm²/s（几乎不传输）。综合CNT和空隙的通量后，H₂/CH₄选择性=18，与实验值15偏差仅20%。分子动力学扩散模拟的这个修正说明模型完整性对结果的准确性至关重要——漏掉次要通道会导致主要指标偏差一个量级。

(12,12) CNT的结果对比：H₂ D=3.5×10⁻⁵ cm²/s（比小管快2倍，因为受限效应减弱），CH₄ D=8×10⁻⁶ cm²/s（比小管快近700倍，因为大管内CH₄不再被完全阻挡）。H₂/CH₄扩散选择性降到4.4——大管失去选择性。分子动力学扩散模拟的数据对客户的设计决策有直接指导意义：如果目标是H₂/CH₄高选择性分离，必须用(8,8)或更小的CNT；如果目标是高H₂通量（牺牲选择性），(12,12)更优。

温度效应分析也做了。25°C和80°C对比：(8,8)管内H₂扩散系数从1.8×10⁻⁵提升到2.3×10⁻⁵ cm²/s（增幅28%），CH₄从1.2×10⁻⁸提升到3.5×10⁻⁸ cm²/s（增幅近3倍）。CH₄扩散的温度敏感性远高于H₂——因为CH₄的受限扩散受吸附控制，升温减弱吸附释放更多CH₄。分子动力学扩散模拟中温度效应的分析为客户评估膜在不同操作温度下的性能变化提供了量化依据。

计算资源方面，单次50ns MD在GPU上跑约6小时，4种管径×2种温度=8个工况共48小时。分子动力学扩散模拟的计算量在纳米通道体系中属于中等规模，但多工况对比是必要的——单一工况的数据无法建立性能-结构关系。

https://www.keyanxueshu.com/