《Molecules》：Functionalized Fullerene Nanomaterials: Evaluating Heteroatom Identity for Enhanced Charge-Transfer and Reactivity
Abdullah M. S. Alhuthali,
Khaled S. Amin,
Hanan Elhaes and
Medhat A. Ibrahim
编辑推荐：
本研究致力于解决如何有效调控富勒烯（C60）的电子特性以提升其功能化应用潜力的问题。研究人员通过密度泛函理论（DFT）计算，系统研究了O、S、Se杂原子功能化对C60结构、电子及反应性质的调控机制。结果表明，功能化有效降低了HOMO-LUMO能隙，增加了总偶极矩，并通过QTAIM、NCI等分析揭示了二聚体中的弱相互作用。这项研究为设计基于C60的高性能光电材料和传感器提供了理论基础。
引言：探寻“巴基球”的无限可能
自1985年被发现以来，足球状的C60分子（又称巴基球）就以其完美的球形结构和独特的电子性质，在材料科学、生物医学和能源领域展现出巨大潜力。它像一个由60个碳原子编织成的微型“笼子”，兼具高对称性、高电子亲和力和化学多功能性。然而，纯净的C60虽然稳定，但其完美的对称性也限制了其化学活性和电子特性的进一步开发。科学家们一直在思考：如何“撬动”这个稳定的碳笼，精准地调控其电子特性，从而为更高效的光电器件、传感器或药物递送系统服务？
其中一种有效策略是功能化，即在C60的表面引入其他原子。在元素周期表的“氧族”（又称硫属元素）中，氧（O）、硫（S）、硒（Se）因原子尺寸和电负性（衡量原子吸引电子能力的标度）各异，成为理想的“改造工具”。例如，氧的电负性最强，能强烈极化电荷分布；而硫和硒拥有更大的原子半径和更易极化的电子云，可能有利于电荷转移。尽管已有大量研究关注单一杂原子的修饰，但对于多杂原子组合修饰，特别是利用这些杂原子将两个C60笼连接成二聚体的结构，其详细的电子特性和笼间的结合性质仍然有待深入探索。这项发表于《Molecules》的研究，正是为了系统地解答这些问题：用O、S、Se功能化C60，究竟如何改变其电子结构、反应活性和分子间作用力？这些改变又预示着哪些潜在的应用前景？
关键技术方法
研究人员主要运用了基于密度泛函理论（DFT， Density Functional Theory）的计算化学方法。所有模型结构均在B3LYP/6-31G(d,p)的理论水平下进行几何结构优化和电子性质计算。研究计算了包括总偶极矩（TDM）、前线分子轨道（FMOs, HOMO-LUMO能隙ΔE）、全局反应性描述符（如电离势、电子亲和能、化学势、硬度、软度、亲电性指数）、总态密度（TDOS）与重叠态密度（OPDOS）、分子静电势（MESP）。此外，还采用了基于电子密度的量子理论（QTAIM, Quantum Theory of Atoms in Molecules）和非共价相互作用（NCI, Non-Covalent Interaction）分析，结合约化密度梯度（RDG, Reduced Density Gradient）映射，来精确表征体系中的化学键和非共价作用力的本质。模型分子从富勒烯C60出发，通过在其C-C键上引入杂原子桥（如环氧结构）或作为笼间连接体，构建了单桥、双桥、三桥（包括O、S、Se组合）以及C60-X-C60（X=O, Se）二聚体等共十种衍生物。
研究结果
1. 优化结构
几何优化显示，C-X键长和C-X-C键角随杂原子尺寸增大而规律性变化（O < S < Se），表明杂原子身份显著扭曲了富勒烯骨架。在二聚体C60-O-C60和C60-Se-C60中，优化的结构显示出笼间分离距离较大（约3.5 ?），预示着非共价相互作用。
2. 计算的电子参数
2.1. 总偶极矩和HOMO-LUMO能隙
原始C60具有完美的对称性，总偶极矩（TDM）为0德拜，HOMO-LUMO能隙（ΔE）为2.762 eV。所有功能化都打破了对称性，增加了TDM并降低了ΔE。例如，C60-O-S-Se的TDM最高（2.156 Debye），而C60-Se的ΔE最低（2.532 eV），表明极性和反应性增强。
2.2. 前线分子轨道（FMOs）
FMOs分析表明，大多数体系的最高占据分子轨道（HOMO）和最低未占分子轨道（LUMO）在整个C60骨架上广泛离域。但在某些体系中观察到轨道重分布，如C60-O-S的HOMO集中在硫桥上，而C60-O-C60二聚体中HOMO和LUMO分别定域在不同的笼上，这建立了定向的电荷转移路径。
2.3. 总态密度（TDOS）
TDOS图显示了所有体系的能隙。重叠态密度（OPDOS）分析显示，在导带区存在负的OPDOS值，表明杂原子与C60片段之间存在反键特性。这在传感器应用中可能有重要意义，因为分析物吸附或光激发可能占据这些反键态，从而削弱连接键，导致电导率变化。
2.4. 全局反应性描述符
计算得到的全局反应性描述符表明，与原始C60相比，所有功能化衍生物的电离能（I）降低，电子亲和能（A）和亲电性指数（ω）增加，化学硬度（η）降低而软度（S）增加。这证实了功能化总体上增强了C60的化学反应性。
3. 分子静电势（MESP）
MESP图以颜色编码（红负蓝正）显示了分子表面的静电势分布。功能化，特别是C60-O-Se，在杂原子位点周围诱导了最强的表面电势变化（激活），表明反应性增强。氧的强电负性产生更明显的负电势区域，而硫和硒则产生相对更正或更中性的区域。
4. 拓扑与非共价相互作用分析
4.1. 分子中原子量子理论（QTAIM）
QTAIM分析显示，在单桥和双桥衍生物中，杂原子与碳原子之间存在键临界点（BCPs），其电子密度ρ(r)较高，且拉普拉斯值?2ρ(r)和总能密度H(r)为负，表明主要是共价相互作用。然而，在二聚体C60-O-C60和C60-S-C60中，连接两个笼的“桥”处观察到较低的ρ(r)和正的?2ρ(r)与H(r)值，结合较大的笼间距离，表明存在的是非共价、范德华（van der Waals）类型的弱相互作用，而非稳定的共价桥。
4.2. 非共价相互作用（NCI）与约化密度梯度（RDG）
NCI-RDG分析进一步证实了上述结论。在所有结构中，C60笼周围可见红色等值面，表明存在空间位阻。更重要的是，在二聚体C60-O-C60和C60-Se-C60中，两个笼之间观察到广泛的绿色等值面，明确指示了稳定结构的范德华相互作用。C60-Se-C60显示出比C60-O-C60更显著的绿色区域，表明其具有更强的色散稳定作用。
5. 结合能
结合能计算显示了一个清晰的周期性趋势：单原子功能化的结合强度遵循O > S > Se的顺序。多杂原子体系（如C60-O-S-Se）显示出更高的结合能（-12.87 eV），表明存在协同稳定效应。二聚体结构（C60-X-C60）的结合能与其单原子对应物相当，这与其结构中存在的强范德华相互作用而非共价键的结论相一致。
结论与讨论
本项研究通过系统的DFT计算，深入揭示了氧、硫、硒杂原子功能化对富勒烯C60结构、电子和反应性质的深刻影响。核心结论可归纳为以下几点：
首先，所有类型的杂原子功能化都有效打破了C60的完美对称性，导致总偶极矩（TDM）增加和HOMO-LUMO能隙（ΔE）减小，其中三桥联的C60-O-S-Se具有最高的TDM（2.156 Debye），而C60-Se具有最低的ΔE（2.532 eV）。这种变化与全局反应性描述符（如降低的电离能、增加的亲电性和软度）所揭示的趋势一致，共同指向功能化显著增强了C60的化学活性和极性。
其次，电子结构分析提供了更微观的图景。前线分子轨道（FMO）的重分布，特别是在某些多杂原子和二聚体体系中出现的轨道局域化现象，暗示了定向电荷转移途径的可能性。总态密度（TDOS）和重叠态密度（OPDOS）分析进一步表明，杂原子对导带附近的电子态有重要贡献，并且在导带区存在反键特性，这为开发基于电荷转移机制的光电或传感应用提供了理论依据。分子静电势（MESP）图谱直观地展示了功能化如何在分子表面创造出具有不同静电势的反应“热点”，尤其是C60-O-Se表现出最强的表面活化。
尤为重要的是，通过量子理论原子在分子中（QTAIM）和非共价相互作用（NCI）分析，本研究澄清了二聚体体系中的键合本质。研究发现，初始建模为共价桥联的C60-O-C60和C60-S-C60在几何优化后，两个富勒烯笼发生分离，最终由范德华相互作用而非共价键稳定在一起。这表明，试图通过简单的杂原子“桥”来共价连接两个C60笼在热力学上可能并不稳定，但可以形成由弱相互作用稳定的独特组装体。
最后，结合能分析量化了这些体系的稳定性，单原子相互作用的强度顺序为O > S > Se，而多杂原子结构由于笼-原子和杂原子-杂原子之间的协同作用而表现出额外的稳定化。
综上所述，这项研究通过一系列互补的计算分析方法，建立了一套从几何结构、电子性质到键合本质的完整评估框架。研究结果表明，通过精心选择杂原子（O, S, Se）及其排列方式（单、双、三桥联），可以像“调音”一样，精准地调控C60的电子特性、表面反应性和分子间作用力。其中，C60-O-S-Se和C60-Se等衍生物在增强极性、降低能隙和活化表面方面表现突出。这些理论见解为后续实验上设计和合成高性能的富勒烯基材料奠定了坚实基础，特别是在需要高效电荷分离和转移的有机光伏（OPV, Organic Photovoltaics）、光电器件以及高灵敏度气体传感器等领域，具有明确的应用前景和指导价值。
