关于预测轻质气体在超临界二氧化碳中自扩散系数的分子动力学方法研究的学术报告
本文旨在介绍一项近期发表于国际期刊《The Journal of Supercritical Fluids》（2025年）的研究工作，题为“Predicting Self-Diffusion Coefficients of Light Gases in Supercritical Carbon Dioxide: A Molecular Dynamics Method”。该研究由西安交通大学动力工程多相流国家重点实验室的Tongjia Zhang, Bowei Zhang, Jie Zhang, Xiaoran Rong及通讯作者Hui Jin*合作完成。研究聚焦于碳捕集、利用与封存及超临界二氧化碳压裂等关键能源与环境技术中的基础传质问题，运用分子动力学模拟探究了氢气、氧气、一氧化碳和甲烷在超临界二氧化碳溶剂中的扩散行为，并开发了高精度的预测模型。
一、 研究的学术背景与目标 随着全球温室气体排放持续增加，降低大气二氧化碳浓度已成为一项紧迫的全球性挑战。二氧化碳的捕集与有效利用是应对这一挑战的核心策略。其中，二氧化碳甲烷化、碳捕集利用与封存以及超临界二氧化碳压裂等技术展现出巨大潜力。在这些技术过程中，超临界二氧化碳常作为反应介质或工作流体，而氢气、氧气、一氧化碳、甲烷等轻质气体作为溶质广泛参与其中。理解这些气体在超临界二氧化碳中的扩散机制，对于优化工艺设计、提高反应效率与封存安全性至关重要。
然而，目前针对高压条件下轻质气体在超临界二氧化碳中扩散行为的研究十分有限，尤其缺乏专门针对超临界二氧化碳溶剂体系的自扩散系数数据与预测方程。实验测量高压条件下的扩散系数往往成本高昂且技术难度大。因此，分子动力学模拟作为一种有效的补充手段，可以突破实验条件的限制，在原子尺度上深入研究扩散过程。本研究的主要目标正是：利用分子动力学模拟方法，在宽温度范围和高压力条件下，计算几种轻质气体在超临界二氧化碳（无限稀释条件）中的自扩散系数；系统分析温度、密度、粘度、分子质量及溶质-溶剂相互作用能等因素对扩散行为的影响；并基于模拟数据，构建能够准确预测自扩散系数的经验方程，以填补该领域的数据与模型空白。
二、 详细的研究流程与方法 本研究遵循了一套严谨、系统化的分子动力学模拟与数据分析流程，主要包含以下几个核心步骤：
1. 力场验证与模拟方法可靠性确认 这是确保后续所有模拟结果准确性的基石。研究首先在超临界区域（以373 K为代表温度）比较了四种常用的二氧化碳力场模型：EPM2、TraPPE、TraPPE-flex和Cygan模型。通过构建包含1000个纯二氧化碳分子的体系，计算自扩散系数并与文献实验数据对比。结果表明，EPM2模型在所有模型中预测最为准确，因此被选为本研究中二氧化碳溶剂的力场。
随后，研究进一步验证所选EPM2力场及模拟方法的普适性。首先，在更宽的温度和压力范围（包括接近二氧化碳临界点的敏感区域）模拟纯二氧化碳的自扩散系数，并与多个来源的实验数据进行对比。为了避免在近临界区因压力波动导致的密度偏差，研究采用了NVT（恒温恒容）系综而非NPT（恒温恒压）系综进行模拟，密度值从美国国家标准与技术研究院数据库获取。对比结果显示，EPM2力场的计算结果与实验数据高度吻合，验证了其在典型超临界条件及近临界区的可靠性。
此外，研究还构建了一个二元体系（5000个CO₂分子和51个H₂O分子），模拟了低压力下无限稀释水在二氧化碳中的扩散，并将模拟结果与文献中的实验数据进行比较。结果再次显示良好的一致性，从而全面证实了所选力场和模拟方法在不同压力条件下的准确性和可靠性。
2. 超临界条件下自扩散系数的计算与体系表征 在验证了基础方法后，研究进入主体模拟阶段。为模拟接近工程实际的无限稀释条件，研究构建了溶剂为5000个CO₂分子、溶质气体摩尔分数为1%的二元混合体系。模拟在宽温度范围（313 K至713 K，间隔50 K）和代表性高压条件（75.994 bar）下进行。在此阶段，研究首先对体系的稳定性进行了评估，通过观察势能随时间的变化确认体系已达到平衡。接着，通过径向分布函数和配位数分析表征了分子的空间分布。结果显示，由于无限稀释条件，所有气体的配位数都非常低，且气体分子在体系中均匀分布，无聚集现象，证实了模拟结果的物理合理性。最后，通过计算分子的均方位移，并应用爱因斯坦关系式以及Yeh-Hummer有限尺寸校正公式，精确计算了各气体溶质的自扩散系数。
3. 扩散行为影响因素的系统分析 获得自扩散系数数据后，研究系统地分析了各种因素对其的影响。 * 温度效应：在恒定压力下，所有四种气体的自扩散系数均随温度升高而单调增加。其中，氢气由于分子质量最小，其自扩散系数最大且对温度变化的敏感性最高。 * 密度与粘度效应：研究发现，自扩散系数的对数与密度的对数之间存在强烈的线性关系，表明自扩散系数随体系密度增加而系统性下降。同时，研究引入了一个组合参数D_self * η / T，发现该参数与温度呈近似线性关系，这揭示了粘度对扩散行为的显著影响，表明温度效应与粘性阻力之间存在竞争关系。 * 分子质量与相互作用能：自扩散系数总体上随溶质分子质量的增加而减小，符合经典扩散理论。溶质-溶剂相互作用能随温度升高而减弱，这与自扩散系数的增加趋势一致。相互作用能的强弱排序与自扩散系数的大小排序相符，表明了内在逻辑的一致性。 * 阿伦尼乌斯关系：研究发现，自扩散系数的对数与温度的倒数呈线性关系，符合阿伦尼乌斯方程。通过拟合得到的自扩散活化能与之前分析的相互作用能排序一致，进一步支持了模拟结果的内在一致性。
4. 经验预测方程的构建与验证 基于以上分子动力学模拟数据，研究团队致力于开发实用且准确的预测模型。 * 快速-安吉尔型幂律方程：首先，研究重新拟合了基于Speedy-Angell关系的幂律方程。该方程仅以温度为自变量，形式简洁。通过对四种气体的模拟数据进行拟合，得到了各自的最优参数。该模型预测的平均相对误差为4.52%，在追求计算简便性的场合具有实用价值。 * 新构建的多参数经验方程：为了更清晰地揭示物理机制并提高预测能力，研究借鉴先前在超临界水体系的工作，构建了一个新的经验方程。该方程明确引入了温度、密度和粘度三个关键物理量，形式为D_self = A0 * T^a / (ρ^b * η^c)。其中，A0是依赖于具体气体的拟合前因子，反映了分子的本征特性及其与溶剂的相互作用；指数a, b, c则反映了各物理量影响的相对强弱。多元线性回归拟合结果显示，温度指数最大，表明自扩散系数对温度变化最为敏感。该方程的整体预测平均相对误差降至4.44%，精度略有提升。
为了检验新方程的普适性，研究进一步在更高压力（101.325 bar和126.656 bar）和扩展温度范围（413 K–713 K）下进行了额外的分子动力学模拟。将新方程的预测结果与这些新的模拟数据对比，平均相对误差分别为4.44%和4.67%，验证了模型在不同压力条件下的良好预测能力。
最后，研究将新建的两个方程与六个常用的经验或基于模拟的关联式进行了对比。结果表明，无论是实验领域常用的Chapman–Enskog、Lusis–Ratcliff、Hayduk–Laudie方程，还是分子动力学领域针对其他体系（如CO₂-H₂O）开发的Moultos、Lu、Zeebe方程，在预测本研究关注的轻质气体在超临界CO₂中的扩散时，均表现出较大偏差。而本研究提出的两个方程，尤其是新构建的多参数方程，展现了显著优越的预测准确性。
三、 主要研究结果 1. 扩散行为的定量规律：研究成功获得了H₂、O₂、CO、CH₄在313–713 K温度范围和75.994–126.656 bar压力范围内于超临界CO₂中的自扩散系数数据库。数据明确揭示了自扩散系数随温度升高而增加、随密度增加而降低的基本规律，并符合阿伦尼乌斯关系。 2. 关键影响因素识别：研究定量分析了温度、密度、粘度、分子质量和溶质-溶剂相互作用能对自扩散系数的综合影响，明确了温度是主导因素，而密度和粘度是重要的制约因素。 3. 高精度预测模型：研究提出了两个新的经验预测方程。重新拟合的幂律方程（平均误差4.52%）适用于以温度为主要变量的快速估算。新构建的包含温度、密度、粘度的方程（平均误差~4.44%-4.67%）具有更高的精度和物理清晰度，其预测性能显著优于文献中现有的多种经典及现代关联式。 4. 方法学验证：研究通过系统性的力场比对和与大量实验数据的交叉验证，确立了采用EPM2力场结合NVT系综的分子动力学模拟方法在研究超临界CO₂体系中气体扩散问题上的可靠性与准确性。
四、 研究结论与价值 本研究的核心结论是：通过经过严格验证的分子动力学模拟方法，首次系统获得了多种轻质气体在宽温压范围内超临界二氧化碳中的自扩散系数数据，并成功建立了两个高精度的经验预测方程。这些方程为相关工程过程（如CCUS、超临界CO₂压裂、CO₂甲烷化）中的传质计算与工艺设计提供了重要的理论工具和基础数据。
其科学价值在于深化了对超临界流体中溶质扩散微观机制的理解，特别是明确了温度、密度、粘度等多物理量的耦合作用规律。所构建的预测模型兼具物理明确性和实用准确性，弥补了该特定体系预测工具的不足。应用价值体现在，这些模型可以用于预测实际工艺条件下气体的扩散速率，从而优化反应器设计、提高封存效率或评估过程动力学，对推动碳中和相关技术的发展具有积极意义。
五、 研究亮点 1. 研究对象的针对性：聚焦于CCUS和超临界CO₂压裂等前沿技术中实际涉及的轻质气体（H₂, O₂, CO, CH₄），填补了该领域基础物性数据与模型的空白。 2. 模拟方法的严谨性：进行了详尽的力场筛选和方法学验证，包括在近临界区采用NVT系综以保持密度稳定的策略，确保了模拟结果的可靠性。 3. 影响因素的全面性：不仅关注温度，还系统分析了密度、粘度、分子特性等多因素的协同影响，对扩散机制提供了更全面的解释。 4. 模型构建的创新性与实用性：提出的新经验方程巧妙地将气体特性（通过A0表征）与溶剂的环境物理量（T, ρ, η）解耦，形式简洁而物理意义明确，且在宽泛条件下保持了优异的预测精度，实用性强。 5. 广泛的对比验证：不仅内部验证了模型在不同压力下的性能，还与多个领域的经典及现代预测公式进行了横向对比，凸显了本研究所开发模型的优越性。
六、 其他值得关注的内容 研究也坦诚指出了当前工作的一个局限性：由于缺乏对应条件下轻质气体在超临界CO₂中扩散的直接实验测量数据，因此所提出的模拟结果和预测方程尚未得到实验的直接验证。然而，作者强调，他们在低压力、高压力以及近临界点等条件下与现有CO₂和H₂O扩散实验数据进行的广泛、一致性良好的对比，已为模拟方法和结果的可靠性提供了强有力的支撑。作者建议，未来结合相应的实验测量将进一步增强所提模型的公信力和普适性。这一实事求是的态度体现了研究的科学性。