N掺杂石墨烯负载Cu–Sc纳米团簇用于CO2捕获的计算筛选研究

《Sustainability》：Computational Screening of N-Doped Graphene-Supported Cu-Sc Nanoclusters for CO2 Capture Katherine Liset Ortiz Paternina and Joaquín Hernández Fernández

　　

将二氧化碳(CO2)转化为高附加值化学品或在排放前进行捕集是缓解大气中CO2浓度上升的互补策略。铜基材料因Cu可结合并电子活化CO2及相关中间体，被广泛研究用于CO2转化。本研究通过构象筛选及298.15、350和400 K下的热化学吸附分析，评估了锚定于三氮原子掺杂石墨烯双层(graphene-3N)上的CuxScγ纳米团簇(Cu3Sc、Cu2Sc2和CuSc3)对CO2的活化作用。首先优化并表征Cu3Sc、Cu2Sc2和CuSc3纳米团簇(相对能量、自旋多重度及电子特性)，并通过比较自由几何与部分受限几何验证载体模型(graphene-3N双层)，借助振动分析确认极小值点。随后评估CO2在CuxScγ@graphene-3N上的吸附/活化，计算ΔH和ΔG值。最终基于ΔG(T)值建立Sabatier(Sabatier)区，发现Cu3Sc呈现中等程度放能吸附(298.15、350和400 K下ΔG分别为?76.07、?67.31和?58.92 kJ·mol?1)；Cu2Sc2呈强吸附(?165.02、?156.36和?148.04 kJ·mol?1)；CuSc3则近乎不可逆固定(?293.98、?287.32和?279.09 kJ·mol?1)。综合活化—再生平衡考量，Cu3Sc为最优团簇组成。

《Sustainability》刊载论文解读：N掺杂石墨烯负载Cu–Sc纳米团簇用于CO₂捕获的计算筛选研究

一、研究背景与意义

大气CO₂持续升高是引发全球变暖的关键因素，单纯收集CO₂不足以解决问题，需将其转化为有价值产物，但CO₂热力学稳定性高、本征反应性低限制了转化过程，核心难点在于将线性CO₂(O–C–O≈180°)活化为极化弯曲的吸附态——即削弱C–O键并向反键轨道转移电荷。Cu因其可结合碳物种并促进多路径CO₂转化被广泛研究，但选择性易受活性位点配位环境等影响；纳米结构化与双金属化(Cu–Sc体系，Sc为具路易斯酸特征的亲氧元素)可调控电子排布(配体效应)、局部几何(集团效应)及吸附键强，但吸附过强会导致活性位堵塞无法再生，需遵循Sabatier原则——相互作用强到足以活化，又不至于阻碍脱附再生。此外氮(N)掺杂石墨烯可提供更强锚定位点和促进电荷转移至吸附物的电子环境。现有文献对纳米团簇上CO₂吸附常忽略随温度升高关键的熵贡献及自旋多重度敏感性(尤其能隙小的体系)，可能改变"最优候选"。因此研究人员采用密度泛函理论(DFT，Density Functional Theory)分析锚于graphene-3N的Cu₃Sc、Cu₂Sc₂、CuSc₃纳米团簇，引入熵校正与自旋态考察，依热力学及电子判据确立具位点再生潜力的CO₂活化体系，提出Sabatier准则优选Cu:Sc配比。

二、主要关键技术方法

研究人员采用ORCA v6.1.0软件进行电子结构全从头算，选用r²SCAN-3c方法(含meta-GGA r²SCAN泛函、def2-mTZVP基组、D4色散校正及几何抵消基组重叠误差(gCP，geometrical CounterPoise)校正)统一优化几何、作用能及振动贡献。构建四原子Cu_4?nSc_n(n=1,2,3)双金属纳米团簇多种初始构型并扫描自旋多重度(单重态/三重态/五重态)；以带空位缺陷且边缘氢饱和的有限石墨烯双层为基底，中心空位处三个C原子替换为N得graphene-3N载体模型。通过谐振频率计算验证极小点、获热化学校正，298.15/350/400 K下用标准平动–转动–振动配分函数求焓(H)与吉布斯自由能(G)，玻尔兹曼(Boltzmann)分布求异构体相对布居。CO₂吸附热力学以ΔX_ads=X(载体+团簇+CO₂)?X(载体+团簇)?X(CO₂(g))计算(X∈{H,?TΔS,G})，假设CO₂为1 atm理想气体并加零点能(ZPE，Zero-Point Energy)校正。电子分析中由前线轨道能(E_HOMO、E_LUMO)求能隙ΔE_gap=E_LUMO?E_HOMO，进而得全局反应性描述符：硬度η=ΔE_gap/2、化学势μ=(E_HOMO+E_LUMO)/2、亲电指数ω=μ²/(2η)；对吸附后CO₂做Hirshfeld电荷分析及量子理论原子中原子(QTAIM，Quantum Theory of Atoms In Molecules)拓扑验证(键临界点BCP，Bond Critical Point处电子密度ρ_c、拉普拉斯?²ρ_c、总能量密度H、动能密度G、势能密度V及|V|/G比)。

三、研究结果

3.1. Analysis of the Structure, Stability, and Electronic Characteristics of Cu–Sc Clusters(Cu–Sc团簇的结构、稳定性及电子特性分析)

研究人员对Cu₃Sc、Cu₂Sc₂、CuSc₃各多初始异构体做DFT全优化及自旋态扫描，借振动分析确证极小点。以室温k_BT≈0.025 eV为参照，相对能ΔE<~24.1 kJ·mol^?1(0.25 eV)的异构体可共存。Cu₃Sc全局最小为异构体5(单重态M=1，准平面菱形¹C_2v)，与异构体6近简并(ΔE=0.03 kJ·mol^?1)，另有ΔE≈0.21 kJ·mol^{?1局域极小及~5.7–5.8 kJ·mol?1三重态可热访问；Cu₂Sc₂基态为三重态异构体4(3C_2v，ΔE=0.00)，伴ΔE≈1.06–2.55 kJ·mol?1低能竞争及~6.7–21.5 kJ·mol?1其他异构体；CuSc₃基态为三重态异构体9(ΔE=0.00)，存近简并(ΔE≈0.01–0.05 kJ·mol?1)及~1.6–1.7 kJ·mol?1低能极小。表明各化学计量比最稳团簇为全局极小，但因存低于~24 kJ·mol?1竞争异构体，室温应描述为多低能态(依赖多重度)共存而非单一构象主导。}

3.2. Thermal Population of Isomers by Boltzmann Distribution(基于玻尔兹曼分布的异构体热布居)

依ΔE作Boltzmann分布求各异构体分数。Cu₃Sc单重态组份298.15/350/400 K分别占91.77%/88.86%/86.21%，均分于7个低能异构体(各约12–14%)，三重态随温升贡献增大(8.26%→13.77%)，五重态因ΔE>~58 kJ·mol^?1基本无布居；Cu₂Sc₂三重态为主(77.81%→70.86%)，集中于3个近等能极小(各约23–26%)，单重态恒定(~16–18%)分布于4个紧邻异构体(ΔE≈1.65–2.55 kJ·mol^{?1)，五重态随温升从6.02%增至11.31%(ΔE≈6.7–6.8 kJ·mol?1)；CuSc₃三重态占优(82.15%→79.42%)均分4个近简并极小，单重态占17.83–20.35%(ΔE≈1.60–1.68 kJ·mol?1)，五重态因ΔE≈20.35–20.37 kJ·mol?1布居≤0.21%。结论：平衡布居受ΔE≤6 kJ·mol?1有限组极小控制，高能组(ΔE≥35–45 kJ·mol?1)贡献可忽略。}

3.3. Nitrogen-Modified Graphene Surface(氮改性石墨烯表面)

3.3.1 完全几何优化较部分约束优化省时且无虚假振动模，选为标淮方案。3.3.2 限放松热力学(ΔX_relax=X_free?X_restricted)示未掺杂空位ΔH_relax≈?40.5至?40.7 kJ·mol^?1、ΔG_relax≈?38.9至?38.6 kJ·mol^?1；graphene-3N对应值更小(ΔH_relax≈?5.48至?5.51 kJ·mol^?1，ΔG_relax≈?4.01至?3.55 kJ·mol^{?1)，温变影响微小，ΔG变化主因熵项(?TΔS_relax)。3.3.3 前线轨道分析得空位石墨烯ΔE_gap=0.872 eV，graphene-3N显著扩至1.581 eV；全局描述符示graphene-3N具最高硬度η=0.791 eV、最低亲电指数ω=7.99 eV(电子级最稳、极化/受荷倾向低)，空位石墨烯μ=?4.0915 eV(相对电子可用性更高)。静电势(MEP，Molecular Electrostatic Potential)图显示掺杂诱导缺陷区静电梯度差异。负载团簇后Cu₃Sc@3N-graphene与Cu₂Sc₂@3N-graphene具较小HOMO–LUMO能隙(ΔE=0.219、0.231 eV)与硬度(η=0.109、0.115 eV)，CuSc₃@3N-graphene ΔE=0.518 eV、η=0.259 eV，三者化学势接近(μ=?3.028至?3.041 eV)，亲电指数ω分别为0.502、0.498、1.198 eV。说明前两者电子结构更软易重整。}

3.4. Nanocluster Adsorption and CO₂Adsorption in Graphene-3N-Nanocluster Complex(graphene-3N–纳米团簇复合物中的团簇吸附与CO₂吸附)

选各团簇全局极小锚于graphene-3N评估CO₂吸附热力学。ΔH_ads很负且温漂<1 kJ·mol^{?1(化学作用主导)，熵项?TΔS_ads随T增大(约49→65 kJ·mol?1)使ΔG_ads负性减弱但仍放能。吸附强度随Sc含量增加：Cu₃Sc ΔG_ads(298.15/350/400 K)=?76.07/?67.31/?58.92 kJ·mol?1(中等放能)；Cu₂Sc₂为?165.02/?156.36/?148.04 kJ·mol?1(强吸附)；CuSc₃为?293.98/?287.32/?279.09 kJ·mol?1(近不可逆固定)。依Sabatier解读三体系均处放能区，Cu₃Sc吸附最适中利活化–再生循环，后二者过强增活性位阻塞风险。}

3.5. Activation of CO₂Adsorbed on Cu_xSc_γ@3N-Graphene: Geometric Descriptors, Charge Transfer, and QTAIM Topological Validation(Cu_xSc_γ@3N-石墨烯上吸附CO₂的活化：几何描述符、电荷转移及QTAIM拓扑验证)

气态CO₂为线性(∠O–C–O=180°，d(C–O)≈1.18 ?)；吸附后均失线性(∠O–C–O≈125.55–130.63°)、C–O伸长(d(C–O)≈1.218–1.326 ?(Cu₃Sc、Cu₂Sc₂)及1.270–1.282 ?(CuSc₃))，CO₂片段获净负电荷(ΔQ(CO₂): Cu₃Sc≈?0.449 e，Cu₂Sc₂≈?0.316 e，CuSc₃≈?0.412至?0.425 e)，属活化化学吸附。Cu₃Sc电荷转移最显且C–O不对称伸长示更强极化吸附机制；CuSc₃活化特征略弱、C–O伸长均衡。几何/电荷指标在298–400 K内几乎不变，吸附态稳。QTAIM示Cu₃Sc与Cu₂Sc₂体系存CO₂–金属键径及BCP，?²ρ_c>0、H<0、|V|/G≈1.24–1.25(闭壳层极化/配位相互作用，典型化学吸附)；Cu₃Sc的BCP处ρ_c=2.45×10^?1、H=?5.84×10^?2略高于Cu₂Sc₂(ρ_c=2.00×10^?1、H=?5.10×10^?2)，|V|/G比近同，印证Cu₃Sc上CO₂–团簇相互作用稍强且极化更明显，与Hirshfeld电荷转移一致。

四、讨论与结论翻译

讨论指出本研究并非宣称Cu₃Sc为新成分发现(与Szalay等人气相筛选一致)，而是验证并将此活性基元拓展至负载材料情境——团簇行为现受团簇–载体耦合调控；N掺杂载体非惰性旁观者，提供强锚定、调局域电子分布、稳化不同于气相有利之团簇几何/电子态，故吸附热力学与活化信号应视为Cu_xSc_γ@graphene-N体系的涌现属性。此视角可分析气相筛选未捕捉的两点：(i)N位点锚定如何重构异构体/自旋态层级；(ii)载体如何调节工作条件下吸附物稳定与有效吸附强度的平衡，从而建立基于团簇–石墨烯–N协同作用的(而非仅金属组分)设计框架。

结论翻译：Cu_xSc_γ化学计量比与graphene-3N载体工程共同调控CO₂黏附强度与电子可及性。N掺杂使HOMO–LUMO能隙扩至ΔE=1.58 eV(空位为0.87 eV)，增硬η=0.79 eV、降亲电ω=7.99 eV，稳化载体并使活化集中于纳米团簇。所评体系中Cu₃Sc在催化循环中达最佳折衷：中等放能吸附(298.15、350、400 K下ΔG_ads=?76.07、?67.31、?58.92 kJ·mol^?1)，熵惩罚随T升高削弱自发但不完全抑制；Sc含量升高致过强键合——Cu₂Sc₂与CuSc₃使吸附阱加深(