从头算 insights into 等离子体与强场效应对银纳米壳上H2解离的贡献

《Physical Chemistry Chemical Physics》：Ab initio insights into plasmonic and strong-field contributions to H2 dissociation on silver nanoshells

　　

研究人员通过将飞秒激光脉冲（10¹³–10¹⁵ W cm^?2）应用于等离子体催化（plasmonic catalysis）的建模，虽然在时域密度泛函理论（time-dependent density functional theory, TDDFT）的时间尺度限制下是合理的，但这可能导致一种与低强度连续光实验（约1 W cm^?2）中等离子体激元激发完全无关的解离机制。在该研究中，研究人员考察了在不同场强、频率和持续时间下，H₂在大型八面体Ag纳米壳上的解离，并探讨了识别当前TDDFT模拟可及的最优建模条件的可能性。研究表明，使用这种大型纳米壳（由Ag₂₃₁团簇的外层组成），即使在外加场强高达(2–8) × 10¹³ W cm^?2时，仍有可能将等离子体激元（plasmon）的作用与强场效应分离开来。具体而言，尽管在这些强度下强场效应始终存在，但研究人员发现，在最低施加强度2 × 10¹³ W cm^?2时，激发的等离子体主导了解离过程。此外，在强场效应占主导地位的最高强度下，等离子体仍有助于加速分子的解离。总的来说，该模拟为弥合等离子体催化中TDDFT建模与实验之间的强度差距铺平了道路。

等离子体催化（plasmonic catalysis）是利用金属纳米颗粒中的局域表面等离子体共振（localized surface plasmon resonance, LSPR）在光的作用下加速和增强化学反应速率的前沿研究领域。实验上已有多项研究表明，当用与纳米颗粒等离子体共振频率相同的光照射时，金属纳米颗粒表面分子解离的速率会增加。对于银纳米颗粒，在共振光照射下观测到化学反应有4至8倍的加速效应。然而，由于等离子体激发和衰减的时间极短（飞秒量级），实验上无法直接分辨，因此仅通过实验无法完全理解等离子体诱导分子解离的精确机制。此外，该机制可能具有系统特异性，仍是当前高度争议的课题。因此，从头计算对于深入理解金属纳米颗粒中等离子体激发与吸附分子解离之间的相互作用至关重要。

时域密度泛函理论（TDDFT）结合埃伦费斯特动力学（Ehrenfest dynamics, ED）是描述纳米尺度飞秒时间尺度上电子和原子激发的先进方法学。实时TDDFT能够准确描述等离子体激发和衰减，而TDDFT-ED可用于描述等离子体衰减后的原子运动。尽管ED作为平均场方法无法描述与不同电子量子态相关的不同路径，但它是能够处理复杂光-物质相互作用以及等离子体衰减后电子和核动力学的少数方法之一，因此仍可为等离子体催化机制提供有价值的洞察。

现有大量TDDFT-ED研究报道等离子体辅助分子解离，但由于计算成本高，这些研究通常对小系统施加高强度激光脉冲（10¹³–10¹⁵ W cm^?2）。这些条件虽加速了解离过程，但会引入强场效应，导致与实验条件（约1 W cm^?2）存在显著差异，使从头计算结果难以外推到实验条件。研究人员此前发现，在小纳米颗粒中，强场效应会抑制等离子体激发的作用并改变分子解离机制。基于此，该研究旨在评估在降低非线性效应的同时，利用当前TDDFT-ED模拟可及的外场强度来研究反应缙变时间尺度的可能性，并深入分析解离机制以识别等离子体及其与强场非线性效应的相互作用。

研究人员选择增大纳米颗粒尺寸来降低强场效应。研究表明，解离阈值强度随纳米颗粒尺寸增大而降低。具体而言，研究人员系统分析了H₂分子吸附在八面体银纳米壳顶点的等离子体效应，该纳米壳由Ag₂₃₁团簇外层构成（记为Ag₂₃₁^L），包含146个原子，内对角线2.45 nm。选择八面体是因其顶点处的场增强效应强于其他形状，可在较低外场强度下增强对分子的作用。通过比较不同脉冲 duration、频率（共振与非共振）和峰值强度，研究人员详细分析了等离子体效应与强场效应的关联。

研究结果表明，在最低施加强度2 × 10¹³ W cm^?2的脉冲-1（pulse-1）作用下，共振频率（?ω_p = 2.48 eV）激发的感应偶极矩幅值显著大于非共振频率（8 eV），且诱导密度分布显示出等离子体激发的集体电子位移特征。然而，偶极矩在脉冲结束后迅速衰减至接近零，表明系统仍处于非线性响应区域。H-H键长演化显示，共振脉冲产生最大伸长，但分子最终脱附而未解离。Mulliken布局分析揭示共振脉冲下分子和纳米壳间存在大振幅电子数振荡，表征等离子体激发特征，同时纳米壳失去超过20个电子，显示强场效应导致的电离。非共振脉冲下分子振荡幅度小、纳米壳电离微弱，证实等离子体可被识别。

延长脉冲 duration的脉冲-2（t₀ = 35 fs, σ = 10 fs）和脉冲-3（t₀ = 55 fs, σ = 15 fs）实现了H₂解离，表明足够长的脉冲 duration是实现解离的关键。脉冲-2和脉冲-3使H₂分别电离约0.6和0.8个电子，纳米壳分别失去35和45个电子。为进一步区分等离子体与强场效应，研究人员比较了总能量相同但峰值强度不同的脉冲-2（I_max = 2 × 10¹³ W cm^?2）和脉冲-4（I_max = 8 × 10¹³ W cm^?2, t₀ = 9 fs, σ = 2.5 fs）。脉冲-2仅在共振时导致解离，是最明确展示等离子体效应的案例；脉冲-4在共振和非共振条件下均导致解离，但共振时解离更快，显示等离子体即使在强场主导时仍加速反应。

诱导电场分析直接证明不同脉冲在系统中的不同扰动。脉冲-2在共振频率下产生高强度感应电场（1.1–1.3 V ?^?1），在分子处局域化；非共振时电场弱且主要局域于纳米壳面。脉冲-4的感应电场更强，非共振时顶点处电场瞬时值（0.9–1.0 V ?^?1）与脉冲-2共振时相当，导致分子解离。这表明外场提供的能量并非唯一决定因素，能量利用效率在共振时更高：等离子体场增强将能量集中于分子附近，实现高效电子发射和键断裂，而非共振激发的能量利用效率较低。

研究所用关键技术方法包括：基于CP2K软件包的密度泛函理论（DFT）几何优化，采用PBE泛函、DZVP基组和GTH赝势；实时TDDFT（RT-TDDFT）计算吸收光谱，通过δ-kick微扰结合离散傅里叶变换获得；TDDFT结合埃伦费斯特分子动力学（TDDFT-ED）模拟电子和原子对外场的响应，使用高斯包络脉冲、ETRS实时传播方案，并添加ghost原子层处理电子发射；数据分析使用NumPy、Matplotlib、VESTA和Gnuplot工具。

研究结论指出，通过增大纳米颗粒尺寸，确实可以区分此前在小纳米颗粒中被强场效应完全抑制的等离子体效应。然而，在所用纳米颗粒尺寸和场强下，非线性效应仍然存在。H₂解离需要分子失去约一个电子，且在短而极强脉冲的非共振条件下也可能发生。但通过与非共振脉冲的系统比较，等离子体效应在所有情况下均可被明确识别：共振激发的等离子体在分子上诱导更大的电荷振荡，产生更强且更局域于分子的感应电场。根据外场强度，等离子体激发可以是引发或加速化学反应的决定性因素。尽管将结果外推到线性区实验条件仍有疑问，但TDDFT-ED已能在此强度和尺寸下解释等离子体诱导光化学过程，为未来研究提供了实用指导：系统比较共振与非共振激发是评估等离子体机制主导性的关键诊断策略；增大纳米颗粒尺寸、延长中等强度脉冲 duration比增大峰值强度更有助于保持等离子体驱动动力学并避免直接多光子电离。这些趋势预期适用于其他材料和几何构型，前提是存在明确的等离子体模式且吸附物位于局域场增强位点附近。