《The Journal of Physical Chemistry A》：Hydrogen Adsorption on Transition-Metal-Decorated Graphene: Thermodynamic and AIMD Insights from DFT and Post-Hartree–Fock Methods
编辑推荐：
研究人员采用有限尺寸石墨烯纳米片层(nanoflake)模型结合revPBE-D4几何优化与ωB97M-V单点能计算，并以PNO–CCSD(T)和DLPNO–CCSD(T)对代表性体系进行基准校验，研究了过渡金属(TM)修饰石墨烯上的氢(H2)吸附行为，特别关
研究人员采用有限尺寸石墨烯纳米片层(nanoflake)模型结合revPBE-D4几何优化与ωB97M-V单点能计算，并以PNO–CCSD(T)和DLPNO–CCSD(T)对代表性体系进行基准校验，研究了过渡金属(TM)修饰石墨烯上的氢(H2)吸附行为，特别关注常被周期性模型忽略的自旋态效应与有限尺寸电子结构的作用。电子结构分析揭示石墨烯纳米片层(G较大、G2较小)存在尺寸依赖的开壳层特征——G呈现稳健的多重组态单重态(multiconfigurational singlet)，G2则表现出较弱且泛函依赖的开壳层倾向；CASSCF/NEVPT2计算支持上述图像并验证了Sc–G与Y–G的高自旋指派。所考察金属中Y、Nb和Pd与石墨烯结合较强(?0.67至?1.02 eV)，Sc、Zr和Ni呈中等稳定化。H2吸附具强烈体系依赖性：首分子H2吸附对Y、Zr、Nb、Ni和Pd为放能(?0.56至?1.63 eV)，Sc则需第二分子H2吸附才变得有利；相反Y、Nb和Ni的二次H2摄取通常不利，仅Zr二次步骤仍适度有利(?0.32 eV)。选取跨越不同成键机制——分子吸附(Pd)、解离氢化物形成(Zr)及高自旋构型(Sc)——的三个代表性体系评估方法学，ωB97M-V//revPBE-D4与耦合簇相互作用能在这些情形下吻合良好。对Sc–G–2H2和Pd–G–H2的AIMD模拟揭示不同的有限温度响应：Sc体系中一分子H2脱附且金属发生表面迁移，而Pd体系中H2保持配位伴随H–H键伸长。综上结果表明，TM修饰石墨烯纳米片层上的氢吸附敏感依赖于金属种类、自旋态及有限尺寸电子结构的协同作用。
论文解读：过渡金属修饰石墨烯上氢吸附的热力学与AIMD洞察——基于DFT与后Hartree–Fock方法的研究
该文发表于The Journal of Physical Chemistry A，针对氢能储存中固体吸附材料的设计需求，指出本征石墨烯与H2作用过弱（≈0.05–0.1 eV），需经过渡金属(TM)修饰以增强结合。现有密度泛函理论(DFT)研究多采用周期性模型且泛函与色散校正选择不一，导致吸附能报道差异大，且普遍忽视石墨烯纳米片层的有限尺寸效应(finite-size effect)、边缘诱导开壳层电子结构及TM–石墨烯复合物自旋态对吸附热力学的影响，亦缺乏跨金属的一致性高精度框架比较吸附热力学与动力学。为此，研究人员采用有限尺寸石墨烯纳米片层模型结合高精度DFT与后Hartree–Fock基准计算，系统考察Sc、Y、Zr、Nb、Ni、Pd修饰石墨烯的电子结构、TM–石墨烯结合能及分步H2吸附的自由能与电子相互作用能，辅以从头算分子动力学(AIMD)探讨有限温度行为，明确自旋态与纳米片层尺寸效应对氢吸附的关键影响，并为TM–石墨烯储氢体系提供接近化学精度的能量基准。
主要关键技术方法
研究人员选用ORCA 6.0.1程序包，以92原子(C92，记作G)与84原子(C84，记作G2)扶手/锯齿混合型石墨烯纳米片层为基底，采用revPBE-D4/def2-mTZVPP进行几何优化并经无虚频验证，以ωB97M-V/def2-mTZVPP做单点能校正；选取Pd–G–H2、Zr–G–H2、Sc–G–H2在DLPNO–CCSD(T)与PNO–CCSD(T)/def2-mTZVPP水平作基准比对。所有相关自旋多重态均被探索，基态自旋由总电子能与Gibbs自由能确定；CASSCF/NEVPT2(活性空间CAS(12,12)或CAS(11,12))验证石墨烯纳米片层及Sc–G/Y–G自旋态指派。H2吸附逐步进行几何优化，按公式(2)计算298.15 K与1 atm下逐次吸附的Gibbs自由能ΔG(n)，按公式(3)计算平均电子相互作用能ΔE(n)；热化学校正取准谐转子近似(quasi-RRHO)。NBO(Natural Bond Orbital)分析由TURBOMOLE完成。AIMD模拟于revPBE-D4/def2-mTZVPP水平在NVT系综以Nosé–Hoover链恒温器进行14 ps(Δt=1 fs)，初始速度按300 K与500 K抽样。
3.1. Electronic Structure of Graphene Nanoflakes（石墨烯纳米片层的电子结构）
研究人员计算两种尺寸石墨烯纳米片层(G与G2)不同自旋态的相对能量，发现较大G纳米片层在所有DFT水平均偏好开壳层单重态(OSS)，较三重态与闭壳单重态分别低0.54 eV与1.69 eV(ωB97M-V//revPBE-D4)，SCF稳定性分析确认OSS比闭壳解更稳定；CASSCF(12,12)/NEVPT2给出自旋匹配单重态最低（较三重态低0.19 eV），自然轨道占据数(如1.3366/0.6650)表明明显双自由基(biradicaloid)及中度多组态特征。较小G2的OSS与闭壳单重态能差缩小(ωB97M-V//revPBE-D4下OSS仅低0.10 eV)，开壳层倾向较弱且泛函敏感——revPBE-D4下OSS解坍缩为自旋纯单重态，ωB97M-D4下二者近简并。自旋密度图中G的自旋密度在石墨烯骨架上较离域，G2则更局域化。结论：石墨烯纳米片层电子结构具尺寸依赖性，大尺寸G具稳健开壳层/多组态单重态特征，可调制金属结合与吸附能。
3.2. Metal–Graphene Interaction Energies（金属–石墨烯相互作用能）
研究人员计算各TM在最优自旋态下与G的结合能ΔE，并借CASSCF(11,12)/NEVPT2验证Sc–G(六重态基态，二重态与四重态分别高0.39 eV与0.83 eV)与Y–G(四重态基态，二重态与六重态分别高0.21 eV与0.72 eV)。Y(?1.02 eV, S=3/2)、Nb(?0.67 eV, S=5/2)、Pd(?0.80 eV, 开壳层单重态S=0)结合强；Sc(?0.57 eV, S=5/2)与Zr(?0.35 eV, S=0)中等；Ni(?0.46 eV, S=0*)偏弱；Cr、Fe、Co、Cu、Mo、W结合能正值或近零，吸附不稳定。早期TM(Sc、Zr)常保持高自旋反映石墨烯弱配体场；Pd呈闭壳行为适于催化活化。结论：Y、Nb、Pd、Sc、Zr、Ni可在石墨烯上形成稳定吸附位点，其中Y与Pd结合最强。
3.3. Hydrogen Adsorption on TM–Graphene Complexes（TM–石墨烯复合物上的氢吸附）
研究人员对稳定TM–G依次加H2优化并计算ΔG与ΔE。首步H2吸附对Y(ΔG=?1.09 eV, 解离氢化物形成，Y–H≈1.96–1.97 ?)、Zr(ΔG=?1.63 eV, 完全解离为氢化物Zr–H≈1.86–1.87 ?)、Nb(ΔG=?1.33 eV, 解离氢化物Nb–H≈1.77–1.81 ?)、Ni(ΔG=?0.56 eV, Kubas型部分活化H–H=0.89 ?)和Pd(ΔG=?0.62 eV, 分子吸附H–H=0.86 ?, Pd–H=1.71 ?)为放能；Sc首步ΔG=+0.06 eV不利，但第二步Sc–G–2H2ΔG=?0.29 eV变有利(H–H伸至0.91 ?, Sc–H≈1.91 ?)。二次吸附仅Zr仍有适度放能(ΔG=?0.32 eV, 混合氢化物+弱配位)，其余(Y、Nb、Ni)二次ΔG>0不利，Pd未见稳定二聚H2吸附。Zr第三、四步ΔG微正(0.03、0.02 eV)表明饱和。ωB97M-V//revPBE-D4与DLPNO/PNO–CCSD(T)相互作用能偏差≤0.06 eV(Pd–H2–G、Zr–H2–G、Sc–H2–G)。NBO显示解离型(Zr/Y/Nb)电荷转移至H显著，Kubas型(Pd/Ni)转移小。结论：H2吸附机制因金属而异——Y/Zr/Nb为解离氢化物型，Pd/Ni为分子(Kubas)活化型，Sc需双H2配位才有利；ωB97M-V//revPBE-D4可高精度复现耦合簇基准能。
3.4. Ab Initio Molecular Dynamics Simulations（从头算分子动力学模拟）
研究人员对Sc–G–2H2与Pd–G–H2做300 K与500 K AIMD。Sc体系中一条H2于500 K下脱附(H–Sc距离剧增>400 ?)，另一条保持近Sc配位；Sc中心在石墨烯上有表面迁移迹象。Pd体系中H2保持配位(H–H波动于0.8–1.0 ?, H–Pd≈1.7–1.9 ?)，未完全解离，Pd亦显示表面迁移。势能涨落围绕均值有限振荡，数值轨迹稳定。结论：两体系有限温响应不同——Sc倾向部分H2脱附与金属扩散，Pd维持H2配位伴键伸长，短时AIMD未表现完全解离。
讨论与结论总结（翻译结论部分）
研究人员指出，石墨烯纳米片层(尤其较大G)具稳健开壳层/多组态单重态特征，较小G2开壳层倾向减弱且泛函敏感；CASSCF/NEVPT2确认G具双自由基特征及Sc–G/Y–G高自旋基态指派。Y、Nb、Pd、Ni、Zr和Sc与石墨烯形成有利复合物(?0.35至?1.02 eV)，Cr、Fe、Co、Cu、Mo、W结合弱或不利。首步H2吸附对Y、Zr、Nb、Ni、Pd放能，Sc需二次H2方有利；Y、Nb、Ni二次吸附不利，Zr二次仍适度有利但三、四步近热中性。选取Pd(分子吸附)、Zr(解离氢化物)、Sc(高自旋)体系作基准，ωB97M-V//revPBE-D4与耦合簇相互作用能高度吻合。自旋密度分析显示金属磁性主要定域于TM，但石墨烯尺寸调控自旋离域程度。AIMD揭示Sc–G–2H2中H2部分脱附与金属迁移，Pd–G–H2中H2持续配位伴键伸长。整体表明氢吸附敏感依赖于金属身份、自旋态及石墨烯有限尺寸电子结构；ωB97M-V//revPBE-D4联合纳米片层模型给出与周期性DFT趋势、实验参考值及CCSD(T)相符的相互作用能，是研究TM修饰石墨烯储氢应用的可靠实用方法。
