MS吸附能计算

一、项目背景与科学问题

分子筛（沸石）因其规整的微孔结构与可调的表面化学，被广泛用于气体分离、吸附储运与催化转化。MFI 型分子筛具有三维交叉孔道体系（直孔道与正弦孔道），对甲烷、二氧化碳等小分子表现出选择性吸附能力。理解客体分子在分子筛不同孔道位点的吸附能差异，是设计高效吸附分离材料的关键。

本项目的科学问题：CH₄ 在 MFI 分子筛不同吸附位点（直孔道、正弦孔道、孔道交叉口）的吸附能分别为多少？CH₄ 与 CO₂、N₂、H₂O 等分子在 MFI 中的吸附强度如何排序？力场方法能否给出合理的吸附能并解释实验上的吸附选择性？

二、计算方法与理论背景

2.1 分子力学与力场

与第一性原理不同，Materials Studio 的 Sorption 模块采用分子力学与蒙特卡罗方法，用经验力场（如 COMPASS、CVFF）描述分子筛骨架与客体分子的相互作用。力场将能量分解为键伸缩、键角弯曲、二面角扭转与范德华/静电非键项。对甲烷这类非极性分子，范德华（Lennard-Jones）项占主导，因此吸附能计算速度快、可处理大体系与多个吸附位点。

2.2 吸附能计算

吸附能 E_ads = E(分子筛+客体) − E(分子筛) − E(客体)，单位为 kJ/mol 或 eV（1 eV ≈ 96.5 kJ/mol）。对周期孔道体系，通常采用巨正则蒙特卡罗（GCMC）或固定负载构型优化的方式，对吸附位点逐一进行能量评价。吸附能绝对值越大，吸附越强。

2.3 MFI 孔道结构

MFI 骨架由十元环窗口构成直孔道（约 5.3×5.6 Å）与正弦孔道（约 5.1×5.5 Å），两者交叉处形成较大的笼状空间。不同位点的空间限制与骨架原子排布不同，导致客体分子的吸附能存在位点依赖。

三、计算设置与模型构建

• 软件：Materials Studio，Sorption 模块，COMPASS 力场。
• 结构：MFI 分子筛单胞（正交晶系），骨架 Si 与 O 原子按实验结构固定，电荷用 QEq 方法分配。
• 计算：固定负载（1 分子/胞）下对 CH₄ 在直孔道、正弦孔道、交叉口的吸附构型进行能量最小化，统计吸附能；并比较 CH₄、CO₂、N₂、H₂O 四种客体分子在交叉口位点的吸附能。
• 细节：范德华截断 12.5 Å，Ewald 求和处理静电，能量收敛阈值 1×10⁻⁵ kcal/mol。

四、结果分析与案例图解读

图 6 的左面板（a）给出 CH₄ 在 MFI 三个吸附位点的吸附能柱状图。孔道交叉口吸附能最强，E_ads = −0.44 eV，因为交叉口空间最大，甲烷可与更多骨架氧原子形成范德华接触；直孔道次之，为 −0.38 eV；正弦孔道最弱，为 −0.31 eV，其弯曲几何限制了甲烷的取向自由度。位点间约 0.13 eV 的差异表明，孔道拓扑对吸附有显著影响，交叉口是 CH₄ 的优先吸附位点。

图 6 的右面板（b）给出四种客体分子在 MFI 中吸附能的对比：H₂O 最强（−0.72 eV），因为水与骨架氧形成氢键且分子小、易进入孔道；CO₂ 次之（−0.44 eV），其四极矩与骨架产生较强静电与色散作用；CH₄ 为 −0.38 eV，主要靠范德华力；N₂ 最弱（−0.21 eV）。该排序与实验上 MFI 对 CO₂/CH₄ 分离系数大于 1、对 N₂ 吸附量低的观测一致，说明力场计算能定性乃至半定量重现分子筛的吸附选择性，为气体分离工艺设计提供了直接参考。

两面板结合，左面板展示了吸附能对孔道位点的依赖，右面板展示了吸附能对客体分子的选择性，共同说明 MS 力场吸附能计算在分子筛吸附研究中的实用价值。

五、结论与展望

本项目用 Materials Studio 的 Sorption/COMPASS 力场计算了 CH₄ 在 MFI 分子筛三个位点的吸附能（交叉口 −0.44 eV、直孔道 −0.38 eV、正弦孔道 −0.31 eV），并比较了四种气体的吸附强度。案例图以位点吸附能柱状图与客体对比柱状图双面板呈现。

后续可拓展：(1) 用 GCMC 模拟吸附等温线与吸附热，与实验等温线对比；(2) 考察 Si/Al 比与阳离子交换对吸附选择性的影响；(3) 对同族分子筛（CHA、FAU、LTA）进行高通量吸附能筛选。