水稳定金属有机框架对小型有毒分子吸附性能的计算研究：基于DFT+U方法的电子结构与吸附机制解析

《Physical Chemistry Chemical Physics》：Computational investigation of small toxic molecule adsorption on water-stable metal–organic frameworksOpen Access

　　

高效捕获NH3、H2S、NO2和SO2等有毒气体对环境和人类健康至关重要。金属有机框架（metal–organic frameworks，MOFs）是一类有前景的吸附剂，但其结构、电子性质与吸附性能之间的关联尚未明确，限制了其理性设计。此外，MOFs在潮湿条件下的应用要求设计水稳定的MOFs用于去除目标有毒分子。在本工作中，研究人员采用含Hubbard U校正的密度泛函理论（density-functional theory with Hubbard U correction，DFT+U）系统研究了多种水稳定MOFs（即UiO-66、UiO-66-NH2、UiO-67、MIL-53（Al， Cr， Fe）和MFM-300（Al， Cr， Fe， V））对选定有毒气体的吸附机制。引入Hubbard U校正显著提高了带隙和结合能等性质的预测精度，使计算结果与实验数据更加吻合。吸附主要由吸附质与水稳定MOFs羟基之间的弱物理吸附和氢键作用主导，但NO2、NH3和SO2在UiO-66-NH2上的吸附除外，此时NH2基团中的氢原子是最活泼的结合位点。对UiO-66进行氨基官能团修饰可略微增强其对NH3、NO2和SO2的亲和力，表明官能团在调控相互作用强度和选择性方面具有重要作用。量子原子分子理论（quantum theory of atoms-in-molecules，QTAIM）和非共价相互作用（non-covalent interaction，NCI）分析证实了水稳定MOFs与目标分子之间存在非共价氢键相互作用，且不涉及与金属中心的直接成键。然而，水分子可能有效竞争吸附位点，从而在潮湿条件下限制所研究水稳定MOFs的选择性。本研究强调了Hubbard U校正在改善过渡金属MOFs电子性质和吸附性能预测方面的关键作用，同时也凸显了通过靶向配体修饰来增强有毒分子在水稳定MOFs上弱结合的必要性，同时保持其结构稳定性。

论文解读

一、研究背景与问题

NO₂、SO₂、NH₃和H₂S等有毒气体的排放对生态环境和人体健康构成严重威胁，涉及臭氧层损耗、全球变暖、酸雨及生态破坏等多重危害。活性炭作为传统吸附剂存在物理吸附作用弱、需化学改性且效率随时间下降等固有缺陷。金属有机框架（MOFs）凭借其高比表面积、孔隙率和可调孔结构，在过去十年中成为有毒化学品吸附、储存、分离和传感领域的研究热点。然而，多数高性能MOFs在潮湿环境中易发生晶体相变或结构坍塌，严重制约其实际工业应用。因此，开发兼具水稳定性和高效吸附性能的MOFs材料成为该领域的关键挑战。

尽管已有大量关于水稳定MOFs吸附小型有毒分子的研究，但MOFs化学与结构特性与其吸附性能之间缺乏明确关联，金属节点和配体对吸附强度的影响尚未被系统计算研究。更为关键的是，准确描述过渡金属d电子一直是理论研究的难点——标准密度泛函理论（DFT）中的近似交换关联泛函对含d电子体系描述不佳，导致带隙、结合能等性质预测严重偏离实验值。此外，在潮湿环境下，水分子与有毒分子对吸附位点的竞争吸附问题尚未得到定量评估。基于上述背景，研究人员开展本研究，旨在通过DFT+U方法系统评估多种水稳定MOFs对NH₃、NO₂、SO₂和H₂S的吸附性能，并阐明Hubbard U校正对电子性质和吸附能预测的影响。

二、主要技术方法

研究人员采用自旋极化的DFT计算，使用广义梯度近似（GGA）下的Perdew–Burke–Ernzerhof（PBE）交换关联泛函，并结合D3色散校正。计算采用GBRV高通量超软赝势和布里渊区Gamma点采样，平面波（PW）和电荷密度截断能分别为100 Ry和1000 Ry。对于过渡金属d电子，研究人员通过线性响应方法自洽迭代计算Hubbard U参数：Zr在UiO系列中约为2.8 eV，V在MFM-300(V)中为4.5 eV，Cr在MIL-53(Cr)和MFM-300(Cr)中分别为6.0和6.1 eV，Fe在MIL-53(Fe)和MFM-300(Fe)中分别为6.1和6.2 eV。Al³⁺因无部分填充d轨道而无需U校正。所有计算采用原胞以降低计算成本，原子位置优化至原子力低于0.002 eV ?^?1，使用Quantum ESPRESSO开源软件包执行，并采用Critic2程序进行QTAIM和NCI分析。

三、研究结果

水稳定MOFs的结构与电子性质

UiO-66、UiO-66-NH₂和UiO-67为立方晶系Fmm空间群，MIL-53(M)为正交晶系Imcm空间群，MFM-300(M)为四方晶系I4₁22空间群。研究人员构建了包含114、126和174原子的菱面体原胞分别对应UiO-66、UiO-66-NH₂和UiO-67，以及包含38和72原子的三斜原胞对应MIL-53(M)和MFM-300(M)。DFT+U优化的晶格参数与实验值高度吻合，平均绝对误差（MAE）仅为0.19 ?（1.6%）。引入Hubbard U校正后，带隙预测值增大，并修正了标准DFT对MIL-53(Fe)和MFM-300(V)预测的虚假金属性行为，DFT+U带隙与实验值的MAE为0.44 eV（17%）。MFM-300(Cr)的带隙偏差归因于电子带隙与光学带隙的差异及可能的激子效应。

吸附/结合能

所有吸附质均与MOFs羟基的氢原子相互作用，唯UiO-66-NH₂例外——NO₂、SO₂和NH₃优先与NH₂基团的氢原子结合，而H₂S仍倾向与羟基作用。结合能整体处于中等偏低水平（5.1–13.1 kcal mol^?1），表明物理吸附为主导机制。UiO-66和UiO-67呈现几乎一致的结合能趋势，说明配体长度对吸附亲和力影响甚微。MFM-300(Fe)和MFM-300(V)对多数吸附质表现出最强结合能力，其中NH₃的结合能最高达13.1 kcal mol^?1。氨基官能团化使UiO-66-NH₂对NH₃、NO₂和SO₂的吸附增强，但对H₂S的吸附反而弱于UiO-66。计算得到的NH₃在MFM-300系列上的结合能与实验吸附热数据（7.2–15.5 kcal mol^?1）吻合良好。与MOF-74和HKUST-1等材料的高结合能（30–69 kcal mol^?1）相比，水稳定MOFs虽亲和力适中，但结构稳定性优势显著。

竞争吸附分析

将目标分子结合能与水分子结合能对比发现：NH₃在UiO-66-NH₂、MIL-53(Al)和MFM-300(V)上的结合能仅比水高1.2–2.2 kcal mol^?1；H₂S和SO₂的结合能普遍低于水，仅SO₂在UiO-66-NH₂和MIL-53(Cr)上略高于水；NO₂在UiO-66-NH₂、MIL-53(Al， Cr， Fe)和MFM-300(Cr)上略强于水。UiO-66-NH₂是唯一能在水存在下仍有效吸附NH₃、NO₂和SO₂的材料。总体而言，潮湿环境下水分子对吸附位点的有效竞争将降低这些水稳定MOFs的选择性和吸附效率。

Hubbard U校正的影响

不含U校正的DFT计算普遍高估结合能。U校正使NO₂在MFM-300(V)、MFM-300(Fe)和MFM-300(Cr)上的结合能分别降低12.6、7.0和5.7 kcal mol^?1，SO₂在MFM-300(V)和MFM-300(Cr)上分别降低4.0和2.2 kcal mol^?1，NH₃在MFM-300(Fe)、MFM-300(V)和MIL-53(Cr)上分别降低2.7、2.6和2.5 kcal mol^?1。DFT+U预测的NH₃在MFM-300(V)和MFM-300(Fe)上的结合能（13.1 kcal mol^?1）与实验吸附热（8.4–10.8和8.4–9.6 kcal mol^?1）更吻合，而DFT值（15.7和15.8 kcal mol^?1）明显高估。DFT的过估源于GGA泛函的自相互作用误差对d电子体系静态关联效应的夸大，U校正通过引入在位库仑排斥项改善了局域d电子的描述。

电荷分布分析

DFT+U计算的金属离子电荷（如Zr约2.7，Cr约2.0，Fe约2.0，V约2.1）比DFT结果更接近理论氧化态（Zr⁴⁺、M³⁺），证实U校正提供了更合理的金属中心电子环境描述。吸附前后金属中心电荷几乎不变，排除了吸附质与金属直接成键的可能。羟基氢原子在吸附后仅发生微小电荷变化，与弱氢键相互作用一致。DFT预测的羟基氢电荷变化可达1.0e，而DFT+U预测的变化极小，后者与实验观测的弱结合行为相符。电子密度差分析表明，U校正增强了金属d电子的局域化，减少了吸附过程中的电子转移，从而降低结合能。

四、讨论与结论

QTAIM和NCI分析进一步证实水稳定MOFs与有毒分子之间为非共价氢键相互作用，不涉及金属中心的直接配位。研究结果表明，水稳定MOFs对目标有毒分子的吸附以弱物理吸附为主，金属节点和配体类型对吸附强度的影响有限，但氨基官能团化可适度增强对NH₃、NO₂和SO₂的亲和力。在潮湿环境中，水分子对吸附位点的竞争构成实际应用的重大挑战，其中UiO-66-NH₂表现出相对最优的耐湿选择性。

本研究的核心结论可归纳为：第一，Hubbard U校正对过渡金属MOFs的电子结构和吸附性能预测至关重要，能有效修正标准DFT对带隙和结合能的高估偏差，使计算结果更贴近实验数据；第二，水稳定MOFs对有毒分子的吸附机制以弱氢键物理吸附为主，缺乏与金属中心的直接强相互作用；第三，配体官能团修饰（如氨基化）是调控吸附选择性的有效策略，但需在增强吸附亲和力与保持水稳定性之间寻求平衡。该研究为理性设计兼具水稳定性和高效有毒气体捕获能力的MOFs材料提供了理论基础，并强调了DFT+U方法在含过渡金属MOFs计算研究中的必要性。