《The Journal of Physical Chemistry C》：When Do Band Gap Calculations Agree with Experiments in Monolayer-Protected Cu14 and Au20 Atomically Precise Nanoclusters? A (TD)-DFT Comparison of HOMO–LUMO, Fundamental, Optical, and Electrochemical Energy Gaps
编辑推荐：
本文聚焦原子精确金属纳米团簇（APNCs）的能隙计算与实验测量的关联难题。为解决不同计算方法（DFT/TD-DFT）预测的HOMO–LUMO、基本、光学、电化学能隙与实验值不一致的问题，研究人员以[Cu14H10(MBN)3(PPh3)8]+和Au20(TBBT)16为模型，系统评估了Hartree–Fock交换、长程分离杂化、基组及溶剂模型对能隙的影响。结果表明，纯PBE泛函结合恰当溶剂模型可较好预测能隙，而HOMO–LUMO能隙对方法参数敏感，不宜直接用作实验能隙的近似。该工作为APNCs的计算研究提供了关键方法论指导。
在材料科学的微观世界里，原子精确金属纳米团簇（APNCs）宛如一颗颗璀璨的明星。它们由特定数量的金属原子构成核心，外层包裹着有机配体，拥有精确的分子式和原子结构。这些特性使其在电催化、光催化、光伏、生物成像和传感等领域展现出巨大应用潜力。然而，要精准“驾驭”这些纳米明星，首先需要透彻理解其最基础的电子性质之一——能隙。能隙如同材料的“电子性格”，决定了它如何吸收光、如何进行电子转移，进而决定了其功能。在实验中，科学家们通过紫外/可见光谱测量光学能隙，通过微分脉冲伏安法等电化学技术测量电化学能隙。在计算中，密度泛函理论（DFT）及其含时版本（TD-DFT）是常用的工具，可以预测最高占据分子轨道与最低未占分子轨道之间的能隙（HOMO–LUMO gap）、基本能隙、垂直激发能等。但一个令人困扰的“能隙迷思”长期存在：为什么不同计算方法得到的“能隙”数值差异巨大？HOMO–LUMO能隙、基本能隙、光学能隙和电化学能隙，这些听起来相似的概念之间究竟有何关联与区别？更重要的是，计算得到的哪个“能隙”才能真正与实验测量值对上号？这不仅是理论计算的精度问题，更直接影响到我们能否通过计算可靠地预测新材料性能，加速材料研发进程。
为了解开这个迷思，并为原子精确纳米团簇的计算研究提供清晰、可靠的方法论指南，一组研究人员在《The Journal of Physical Chemistry C》上发表了一项系统性研究。他们选择了两个具有代表性的“模特”：铜氢化物纳米团簇[Cu14H10(MBN)3(PPh3)8]+（简称Cu14）和硫醇保护的金纳米团簇Au20(TBBT)16（简称Au20）。这两个团簇虽然都小，但“性格”迥异：Cu14的第一激发态具有金属到配体的电荷转移（MLCT）特性，而Au20的激发态则完全局域在金属核心上。研究团队对这两个体系进行了全面的实验测量（紫外/可见光谱和微分脉冲伏安法）和大规模的计算扫描。在计算方面，他们系统考察了从纯泛函到杂化泛函（HF交换比例从0%到50%）、长程分离杂化泛函（调节长程参数ω）、不同关联泛函、多种基组，以及平衡与非平衡溶剂化模型（PCM）对四种计算能隙（HOMO–LUMO能隙、基本能隙、光学能隙、电化学能隙）的影响。其核心目标是回答：在原子精确纳米团簇中，计算能隙何时能与实验能隙达成一致？
为了开展这项研究，作者们主要运用了以下几项关键技术方法：首先，利用X射线单晶衍射确定的晶体结构构建初始计算模型，并对庞大的配体（如PPh3、TBBT）进行了合理的截断简化以降低计算成本（例如替换为PH3、SCH3）。其次，核心的计算工作全部在Gaussian 16软件中完成，广泛使用了密度泛函理论（DFT）和含时密度泛函理论（TD-DFT）。再者，为了评估不同理论方法的性能，他们系统性地扫描了PBE泛函中Hartree–Fock（HF）交换的占比（0%至50%），并测试了LC-ωPBE等长程校正泛函。同时，采用ΔSCF（Delta自洽场）方法计算垂直电离能和垂直电子亲和能以得到基本能隙，并进行了完整的氧化还原电势热力学计算以获得电化学能隙。最后，为了贴近实验环境，所有溶液相计算均使用了二氯甲烷的极化连续介质模型（PCM），并区分了非平衡PCM（用于垂直过程如光谱）和平衡PCM（用于平衡过程如电化学）的不同处理方式。
计算结果与讨论
Cu14纳米团簇的能隙比较
研究首先在纯PBE泛函水平上，比较了不同方法计算的能隙与实验电化学能隙。对于截断模型[Cu14H10(MBN)3(PH3)8]+，计算得到的电化学能隙（Eoox– Eored）为2.29 eV，与实验值2.073eV吻合较好。HOMO–LUMO能隙和TD-DFT计算的垂直激发能（VEE）分别为1.84 eV和1.85 eV，彼此接近但略低于实验电化学能隙。而使用非平衡溶剂化计算的基本能隙（VIE-VEA）高达3.33 eV，显著高于其他值。这表明，简单地使用HOMO–LUMO能隙或气相基本能隙来估算溶液中的电化学性质会造成很大偏差。然而，当使用平衡溶剂化模型计算(V)IE-(V)EA时，能隙降至2.54 eV，更接近电化学能隙，说明溶剂弛豫是关联计算与实验电化学能隙的关键因素。
HF交换对四种计算能隙的影响
通过系统改变PBE泛函中HF交换的比例（0%至50%），研究揭示了不同能隙对方法参数的敏感性存在巨大差异。HOMO–LUMO能隙的变化最为剧烈，从0% HF时的1.84 eV飙升至50% HF时的5.74 eV。基本能隙（VIE-VEA）、光学能隙（VEE）和电化学能隙（Eoox– Eored）虽然也随HF交换增加而增大，但其变化幅度远小于HOMO–LUMO能隙。这意味着，未经校准的HOMO–LUMO能隙极易受到所用泛函的强烈影响，不适合直接用于估算实验能隙。此外，随着HF交换比例的改变，不同能隙之间的能量高低顺序也会发生定性变化，进一步说明了随意选择计算方法和参数可能导致错误的物理图像。
溶剂化效应与轨道特性分析
溶剂化模型的选择对计算结果至关重要。对于涉及电荷变化的过程（如电离、电子附着），使用平衡溶剂化模型相比非平衡模型或气相计算，能使基本能隙降低近2.2 eV。这表明在关联电化学实验时，必须考虑溶剂的平衡弛豫。相比之下，TD-DFT计算的光学能隙对溶剂化模型的细节（平衡vs非平衡）不敏感。这种差异的根源在于Cu14激发态的电荷转移特性。前沿轨道分析和自然跃迁轨道（NTO）分析表明，Cu14的HOMO和LUMO以及其第一激发态的NTO空穴-电子对均显示出金属到配体的电荷转移特征。随着HF交换比例的增加，这种电荷转移特性会减弱，导致激发态的偶极矩变化，进而改变了溶剂化能对激发能的影响。这使得在溶液中，(V)IE-(V)EA与(V)EE两条曲线随着HF交换增加不再平行，而在气相中是基本平行的。
Au20纳米团簇的验证
为了检验结论的普适性，研究在Au20体系上重复了分析。尽管Au20的第一激发态是局域在金属核心的（无电荷转移），但其能隙对HF交换的敏感性模式与Cu14惊人地相似：HOMO–LUMO能隙依旧极度敏感，而VEE和(V)IE-(V)EA的变化则平缓得多。在Au20中，由于激发态无电荷转移，VEE始终低于(V)IE-(V)EA，且两者曲线基本平行，这与理论预期一致（光学能隙应低于基本能隙）。在0% HF交换（纯PBE泛函）下，计算得到的光学能隙为2.02 eV，与实验测得的2.12 eV光学能隙吻合良好。这解释了为何在金纳米团簇研究领域，纯PBE泛函被广泛使用并能取得与实验一致的结果。
其他计算参数的影响
研究还考察了关联泛函、基组、长程分离参数等因素的影响。发现相比HF交换比例，纯泛函之间的差异（如LDA、GGA、meta-GGA）对HOMO–LUMO能隙的影响要小得多（约0.6 eV以内）。基组大小的影响也相对较小（约0.13 eV以内）。然而，长程分离杂化泛函中的范围分离参数ω对能隙的影响极为显著，增大ω会导致能隙急剧升高。例如，使用默认ω=0.4 bohr–1的LC-ωPBE泛函，会严重高估实验能隙。
研究结论与意义
本研究通过对两个特性迥异的原子精确纳米团簇（Cu14和Au20）进行系统的计算与实验对照，得出了一系列对领域具有重要指导意义的结论。首先，文章明确了“能隙”一词在不同语境下的具体指代，强调HOMO–LUMO能隙、基本能隙、光学能隙和电化学能隙是彼此关联但不等价的概念。对于封闭壳层/弱关联的铜或金基APNCs，使用纯PBE泛函（0% HF交换）并结合恰当的溶剂化模型（对于光学过程用非平衡PCM，对于氧化还原过程用平衡PCM），能够较好地预测其实验光学和电化学能隙。这为为何该领域大量研究成功使用纯PBE提供了理论依据。
其次，研究强烈警示，HOMO–LUMO能隙对计算方法的参数（特别是HF交换比例和长程校正参数）异常敏感，在未经验证的情况下，绝不应将其直接当作实验能隙的可靠近似。相比之下，通过ΔSCF计算的基本能隙、通过TD-DFT计算的光学能隙以及通过热力学计算的电化学能隙，由于包含了电子弛豫的更多信息，对方法参数的依赖性较低，是更可靠的比较对象。其中，使用平衡溶剂化模型计算的(V)IE-(V)EA可以作为计算量更大的完整电化学能隙（Eoox– Eored）的一个有效近似。
最后，该工作为未来APNCs的计算研究提供了实用的“方法论清单”。在选择泛函时，应优先基于实验数据或高精度波函数方法进行基准测试。在缺乏基准的情况下，建议系统绘制关键能隙（HOMO–LUMO、VEE、(V)IE-(V)EA）随HF交换比例或ω变化的曲线，以评估其敏感性和一致性。溶剂化模型必须根据所模拟的过程（垂直或平衡）正确选择。本研究揭示，对于APNCs这类介于分子与体材料之间的特殊体系，在气相中通过满足Koopmans定理来“优化”泛函参数（如调大ω或HF比例）可能会得到与溶液相实验严重偏离的结果，因此需要格外谨慎。
总之，这项研究如同为纳米团簇计算领域绘制了一幅详尽的“导航图”，它不仅解答了“何时能算准能隙”的核心问题，更通过清晰的范例和系统的分析，为研究人员选择合适的计算工具、避免常见陷阱、从而更可靠地建立原子精确纳米团簇的结构-功能关联奠定了坚实基础。
