太赫兹振动光谱与分子模拟揭示α-氢醌包合物中H2包合的光谱指纹

《ChemComm》：Characterization of hydrogen storage in hydroquinone clathrate using terahertz spectroscopyOpen Access

　　

氢醌（hydroquinone, HQ）包合物（clathrates）是固态储氢（solid-state hydrogen storage）的有前景材料，然而对H₂包合（enclathration）的明确光谱表征仍难以实现。在此，研究人员表明，低频（太赫兹）振动光谱（low-frequency (terahertz) vibrational spectroscopy）结合分子模拟（molecular simulations），可揭示α-氢醌（α-hydroquinone, α-HQ）多晶型（polymorph）中H₂包合物形成的清晰光谱指纹（spectroscopic fingerprints），并展示了一条出乎预期的氢掺入（hydrogen incorporation）途径。

**研究背景与问题**
氢作为燃料需要更高效的气体储存方法。包合物（clathrates）是一类由主体晶格包含分子尺度空隙并捕获客体分子形成的多孔材料，被视为H₂储存候选。常见的水合物（hydrates）在自然永久冻土中存在，但室温不稳定。氢醌（hydroquinone, HQ）包合物因室温稳定性受到关注。HQ有四种多晶型（polymorph），其中α和β可形成包合物；α在室温稳定，β在高压并引入客体分子时形成。β-HQ每3个HQ分子有1个孔，α-HQ为1:18。此前β-HQ-H₂形成困难：HQ从溶液升华时倾向形成空β结构而非吸收H₂；H₂与CO或CO₂二元混合物中，HQ优先包合CO或CO₂而不包合H₂；β-HQ-H₂需高于70 MPa并持续多天。α-HQ可在远低于β包合所需压力下将H₂纳入间隙空间，小分子和惰性气体在升华实验中也倾向进入α间隙而不引发β转变。但此前Raman、核磁共振（nuclear magnetic resonance, NMR）和X射线衍射（X-ray diffraction, XRD）难以确认α结构中孔内H₂的包合：X射线散射不易定位质子，常规中红外（mid-infrared, mid-IR）主要探测数百cm^?1以上的分子内动力学，对软主客体动力学信息有限。弱主客体作用在太赫兹（THz）区（低于约100 cm^?1）表现为耦合主客体动力学，可给出框架多晶型和客体身份的光谱指纹，甚至客体取向。因此，表征H₂包合后HQ的THz谱是阐明相关结构的有前景途径。研究人员据此用低频THz振动光谱结合分子模拟，报道α-HQ中H₂包合物形成的光谱指纹，并指出一条出乎预期的氢掺入途径。该工作首次确认负载H₂的α-HQ的THz谱并指认振动动力学，证实相对低压下α-HQ包合物中H₂的存在与动力学，也展示神经网络振动模式计算与传统密度泛函理论（density functional theory, DFT）结合可获得DFT单独难以计算的结果。

**关键技术方法**
实验制备直径5 mm、高1 mm的研磨α-HQ压片，置于同尺寸高密度聚乙烯（high-density polyethylene, HDPE）压片上，以增加THz脉冲有效路径并减少反射影响；腔室充H₂至7.6 MPa，室温保持24 h，每10 min进行太赫兹时域光谱（terahertz time-domain spectroscopy, THz-TDS）测量，每次7272次扫描平均；吸收谱由样品与HDPE组合的THz功率谱除以单独HDPE谱得到。模拟采用MACE-MP-0机器学习原子间势（machine learning interatomic potentials, MLIPs）在无对称约束下生成并优化每孔2、4、6个H₂的α-HQ结构；对每孔2个H₂且收敛为P3?空间群（space group）的结构，用MLIP获得振动本征矢和本征值，用固体DFT计算Born有效电荷（Born effective charges），结合得到红外强度（infrared intensities）；空载结构采用相同流程。无样本队列。

**研究结果**
**实验THz-TDS吸收光谱**：在7.6 MPa H₂原位加压24 h中，α-HQ的1.18 THz声子模频率基本不变；1.67 THz模立即增强、红移并出现不对称线形。先前工作显示1.67 THz峰随β包合物形成而降低，因此此处增强提示α形式改变而非向β转变；THz谱对α→β会剧烈重构，而实测无剧烈重构，说明主体框架保持α，H₂客体引入仅造成较小光谱变化。
**双Gaussian定量分析**：早期1.67 THz特征较对称，可用单Gaussian拟合（Gaussian fit）；后期不对称性需双组分拟合。固定中心频率后，以两个Gaussian振幅作为1.67 THz原特征和约1.64 THz新特征强度指标。两振幅随时间变化表明低频第二特征与包合过程相关。
**混合MLIP-DFT模拟**：常规固体DFT因H₂加载降低对称性、每孔可容纳多个H₂且局部构型多而计算昂贵。研究人员用MLIP生成优化H₂-loaded α-HQ结构，并以MLIP提供晶格动力学、DFT Born有效电荷提供电动力学响应的混合流程计算红外强度。模拟复现实验最重要的定性特征：H₂掺入不造成α-HQ谱均匀位移；约1.2 THz特征基本不变，而约2.1 THz特征加载后移至约1.8 THz。计算绝对频率和位移幅度与实验不同，但复现了响应选择性，即约1.18 THz模基本不变、约1.67 THz特征出现低频组分，提示两模与客体相互作用方式不同。
**振动模式指认**：约1.2 THz模主要由空框架中HQ分子受阻旋转运动（hindered rotational motion）主导，加载框架中同频运动且H₂参与可忽略，决定该模频率的分子间力未被H₂显著扰动；该特征在H₂加载中持续及模拟中空载/加载频率基本不变支持α-HQ框架保留。较高频模在空α-HQ约2.1 THz、加载后约1.8 THz，涉及相同框架运动；H₂掺入后框架振动与客体运动耦合，频率位移来自有效力常数和/或约化质量（reduced mass）改变：弱主客体作用改变框架有效恢复力，客体参与改变简正模（normal mode）有效质量。因此红移是既有α-HQ晶格振动改变而非全新框架模，这为实验1.67 THz特征演化和约1.64 THz低频组分提供原子尺度解释。
**表面贡献评估**：研究人员考虑光谱演化是否主要来自H₂与α-HQ颗粒外表面相互作用而非进入本征孔。对10 μm直径球形颗粒，最外50 ?（5 nm）仅约占总颗粒体积0.30%，按均匀密度也占相同质量分数；若光谱响应仅来自该近表面区，仅很小部分HQ分子贡献体相THz响应。更重要的是实验响应高度模式特异而非高压H₂引入后的宽谱扰动；孔内加载H₂的计算复现该选择性，即客体不敏感模几乎不变、本征矢直接与H₂耦合的模红移。虽不能排除表面结合H₂贡献，但表面材料分数小、体相晶格振动选择性响应及孔加载模拟复现该选择性，使主要表面来源不太可能。

**讨论总结与研究结论翻译**
讨论指出，本研究首次确认负载H₂的α-HQ的THz谱并指认振动动力学。通过比较THz-TDS实验数据与混合MLIP-DFT模拟，研究人员发现实验所见1.67 THz峰移对应H₂气体加载入α-HQ后出现的耦合振动模。THz-TDS适合研究该问题，因为这些振动模对主客体弱相互作用的小变化高度敏感，比X射线散射等测量提供更多信息。这些发现证实相对低压下α-HQ包合物中储存H₂的存在与动力学；也表明将神经网络振动模式计算与传统DFT模拟结合，可获得DFT单独在计算上难以实现的结果。研究结论翻译：本研究首次确认了负载H₂的α-HQ的太赫兹光谱，并伴随振动动力学的指认。通过将THz-TDS实验数据与混合MLIP-DFT模拟进行比较，研究人员发现实验中观察到的1.67 THz峰移对应于H₂气体加载到α-HQ后出现的一个耦合振动模式。THz-TDS是研究该问题的优秀工具，因为这些振动模式对主客体分子之间弱相互作用的微小变化非常敏感，比X射线散射等其他类型测量提供更多信息。这些发现证实了在相对较低压力下储存在α-HQ包合物中的H₂气体的存在和动力学。最后，研究人员的工作还证明了将振动模式的神经网络计算与传统DFT模拟相结合，以获得DFT单独无法在计算上实现的结果的效用。