氢活化与吸附机制在铬修饰B??N??纳米笼中的研究：结合DFT、分子动力学与CI-NEB计算

《Materials Today Communications》：Hydrogen Activation and Adsorption Mechanisms on Cr-Modified B 12N 12 Nanocages: A Combined DFT, Molecular Dynamics and CI-NEB Study

　　

分子氢（H?）作为一种高能量密度且燃烧无污染的清洁能源载体，其高效储存与安全利用是实现大规模应用的关键挑战。理解过渡金属修饰纳米材料上氢活化与吸附的调控机制，对于氢相关技术的理性设计具有基础性意义。本研究采用B3LYP-D3/6-31G(d,p)水平的密度泛函理论（density functional theory, DFT）计算，系统探究了H?与本征及铬（chromium, Cr）修饰的B??N??纳米笼在掺杂（doped）、修饰（decorated）及包封（encapsulated）构型中的相互作用。通过前线分子轨道、态密度（density of states, DOS）、电荷转移、电子定域函数（electron localization function, ELF）、量子理论原子（quantum theory of atoms in molecules, QTAIM）及约化密度梯度（reduced density gradient, RDG）分析，系统考察了结构、电子及热力学性质。铬功能化显著改变了纳米笼的电子结构，将本征B??N??的HOMO-LUMO能隙从6.88 eV降至修饰体系的2.69–4.95 eV，并增强了其与氢的相互作用。在所有构型中，Cr@b66表现出最强的H?吸附能力，吸附能为?1.63 eV，电子调制率达85%，电荷转移量为?0.379 |e|，表明通过强的金属-氢耦合实现了显著的氢活化。多分子吸附、分子动力学模拟及初步的爬坡弹性带（climbing-image nudged elastic band, CI-NEB）计算进一步揭示了复杂的吸附-脱附行为及氢相互作用伴随的显著结构重组和电子重构。上述发现为铬功能化调控B??N??纳米笼氢吸附提供了机理层面的深入认识，有助于增进对氮化硼基纳米结构中氢活化的基础理解。

一、研究背景与问题

清洁可持续能源的开发推动了以分子氢（H?）为代表的化石燃料替代方案研究。H?凭借高能量密度和燃烧过程零污染物排放的优势受到广泛关注，然而其高效、安全的储存仍是制约大规模应用的关键瓶颈。氮化硼纳米笼（B??N??）凭借高热稳定性、低介电常数、高热导率、抗氧化性及电绝缘特性，在储氢领域展现出应用潜力。实验上，Oku等人成功合成了本征及过渡金属修饰的B??N??和B??N??纳米笼，并探索了其储氢性能。然而，本征B??N??与H?分子间的弱相互作用严重限制了其储氢效率。过渡金属功能化被证实可通过调控B??N??纳米笼的电子结构和几何构型来增强H?吸附。已有DFT研究表明，Ni、Li、Na、K等金属修饰可显著提高纳米笼的H?吸附容量。尽管取得了一系列进展，但在吸附强度、可逆性与稳定性之间取得最优平衡仍是重大挑战。以往研究涉及多种过渡金属和碱金属修饰，但据研究人员所知，铬（Cr）功能化B??N??纳米笼的H?吸附与活化行为尚未被系统研究，不同金属构型的影响、动态稳定性及多分子气体吸附行为亦鲜有报道。Cr因其部分填充的3d电子构型（3d?4s1）、可及的自旋态和多氧化态，能够提供促进电荷再分布和分子活化的多功能电子环境，前期研究已证实Cr修饰B??N??对SO?吸附具有显著的电子响应（电子敏感度79%，Eads = ?0.96 eV），这为探究Cr中心能否活化H?提供了合理依据。基于此，本研究在DFT水平上对Cr修饰的B??N??纳米笼（掺杂、修饰和包封构型）进行理论分析，评估金属位置对H?吸附效率和复合物稳定性的影响，旨在为理解Cr修饰B??N??纳米笼上H?的相互作用机制、电子调制行为和吸附特性提供基础见解。

二、研究内容与主要结论

研究人员对B??N??纳米笼及五种Cr修饰构型——内包封（Cr@B??N??）、B或N替代掺杂（CrB??N??和B??N??Cr）以及外修饰（Cr@b64和Cr@b66）——进行了全面理论计算。结果表明，铬功能化显著调控了纳米笼的电子结构，缩小了HOMO-LUMO能隙（从本征的6.88 eV降至2.69–4.95 eV），增强了体系的反应活性和氢亲和力。在所有构型中，Cr@b66表现出最强的H?吸附能力（Eads = ?1.63 eV）、最高的电子调制率（ΔEgap = 85.1%）和显著的电荷转移（Q_H? = ?0.379 |e|），并诱导H?发生解离吸附，形成Cr-H和B-H部分共价相互作用。分子动力学模拟证实了吸附复合物在近环境温度下的动态稳定性，而CI-NEB计算揭示H?脱附伴随显著的结构和电子重组过程（能垒2.48 eV），表明Cr@b66兼具高效氢活化中心的特征。多分子吸附研究表明，Cr活性位点可稳定吸附约2个H?分子，随后饱和，对应重力储氢容量约1.14 wt%。上述发现阐明了Cr功能化在调控B??N??纳米笼氢吸附中的关键作用，为氮化硼基纳米结构中氢活化的基础理解和氢相关纳米材料的理性设计提供了重要的理论依据。

三、关键研究方法

本研究采用多种理论计算方法。以B3LYP-D3/6-31G(d,p)水平的DFT计算为基础，使用ORCA 5.0软件包，对所有体系进行几何优化和性质计算，并采用Grimme的D3色散校正描述范德华相互作用。研究中使用MultiWfn软件包进行态密度（DOS）分析和拓扑分析（ELF和LOL）。电荷转移和电子结构性质通过Mulliken电荷分析、前线分子轨道分析获得。为进一步揭示相互作用的本质，开展了QTAIM（量子理论原子）分析、ELF（电子定域函数）分析和RDG（约化密度梯度）非共价相互作用分析。分子动力学模拟（MD）采用xTB包在GFN-1哈密顿量下进行，用于评估Cr@b66及Cr@b66–H?体系在298.15 K下的动态稳定性。此外，采用爬坡弹性带（CI-NEB）方法计算H?脱附反应的动力学能垒。研究还通过BSSE校正评估了吸附能，通过Gibbs自由能变化评估了热力学稳定性。样本来源方面，研究对象为本征B??N??纳米笼及其铬修饰衍生物，未涉及实验样本。

四、主要研究结果

结构分析与稳定性：几何优化表明B??N??纳米笼具有高对称性和零偶极矩，B-N键长（b64 = 1.484 ?、b66 = 1.437 ?）与文献吻合。Cr修饰体系的自旋多重度研究表明，Cr@B??N??以单重态最稳定，掺杂体系以三重态最稳定，修饰体系以五重态最稳定。Cr修饰降低了化学硬度（η）并增加了亲电性（ω），表明反应活性增强。所有体系均具有负的内聚能（Ecoh），证实了能量的稳定性，其中CrB??N??最稳定（Ecoh = 7.428 eV）。

H?在本征及Cr修饰B??N??上的吸附行为：几何分析表明，H?在掺杂和修饰纳米笼中吸附于Cr原子上，在包封纳米笼中吸附于B原子上，而在本征B??N??中吸附于b64键上。Cr@b66–H?中H?发生解离吸附，Cr-H和B-H键同时形成。Cr修饰显著缩短了笼-H?距离，Cr@b64–H?具有最短距离（1.443 ?），本征B??N??–H?距离最大（3.016 ?）。H-H键伸长和ν(H-H)振动频率软化证实了H?的渐进活化，Cr@b66–H?中ν(H-H)降至2022.03 cm?1。Gibbs自由能分析表明Cr@b66–H?具有最有利的热力学稳定性（ΔGads = ?1.39 eV）。

态密度和前线轨道分析：DOS分析显示Cr掺入强烈改变纳米笼电子结构，促进HOMO去稳定化和LUMO稳定化，前线轨道主要由Cr贡献，表明存在强烈的金属-笼杂化。H?吸附后，特别是在Cr@b66–H?中，H? σ轨道与Cr 3d轨道之间出现显著重叠，产生新的杂化态，增强电荷离域。Cr@b66的β自旋态Egap变化达85.1%，表明其对H?具有极高的电子敏感度。吸附能计算进一步表明Cr@b66具有最强的相互作用（?1.63 eV），属于强化学吸附。

QTAIM、ELF和NCI分析：QTAIM分析揭示Cr@b66–H?的Cr-H键具有最高电子密度ρ(r) = 0.0946 a.u.，?G(r)/V(r)比值0.6794，总能量密度H(r) = ?0.0310，表明部分共价相互作用。ELF分析显示Cr修饰增强了金属中心附近的电子定域化，Cr@b66–H?中解离氢物种上具有高电子密度，支持强部分共价气体-笼相互作用。RDG分析进一步证实了Cr-H强吸引相互作用（蓝色等值面）和范德华弱相互作用（绿色等值面）在不同体系中的分布。

分子动力学分析：在298.15 K下对Cr@b66和Cr@b66–H?进行了1000 ps的NVT-MD模拟。能量曲线仅表现出小幅度热振荡，无突变或持续漂移。Cr@b66的RMSD为0.3–0.4 ?，表明高动态稳定性。Cr@b66–H?在初始结构重排后稳定在活性位点，RMSD为0.7–1.1 ?，H-H距离波动约0.60 ?，证实了H?的持续活化状态。

脱附时间分析：基于吸附能估算的脱附时间表明，Cr@b66–H?在450 K（黄光频率下）和550 K（紫外频率下）可实现脱附。针对强相互作用体系的CI-NEB计算则揭示出2.48 eV的动力学能垒，对应于H-H重组和Cr-H键断裂过程，远高于简化模型预测值，表明脱附过程涉及显著的结构重排。

H?浓度研究：多分子逐步吸附研究表明，第一个H?分子呈强化学吸附（Eads = ?1.63 eV），第二个处于中等吸附区间（Eads = ?0.33 eV），第三、四个为弱物理吸附，第五个变热力学不稳定。Cr活性位点早期饱和限制了有效储氢容量。基于两个有利吸附H?分子，Cr@b66的重力储氢容量约为1.14 wt%，低于美国能源部（DOE）设定的5.5 wt%实用储氢最低目标。

五、讨论与结论

讨论部分指出，虽然Cr@b66表现出最强的H?吸附能力和显著的活化能力，但其?1.63 eV的吸附能超出了环境条件下可逆储氢的理想范围（约?0.2至?0.6 eV/H?）。Cr@b66本质上更接近高效的分子活化位点，而非最优的可逆储氢材料。第二个H?分子的中等吸附能力暗示有限的储氢可逆性。Cr@b66的主要局限在于尽管对第一个H?表现出强活化能力，但整体吸附容量偏低。优化活性金属中心的数量与分布，或调控纳米笼的拓扑结构，可能成为控制氮化硼纳米结构氢吸附行为的有效策略。该研究深化了对过渡金属功能化纳米笼中氢相互作用的基础理解，为下一代氢相关纳米材料的理性设计提供了理论指导。