# 从磁性视角看分子芳香性

《The Journal of Physical Chemistry Letters》：Molecular Aromaticity from a Magnetic Perspective

　　

芳香性概念最初被用于描述类似苯的化合物。后来，芳香性被扩展用于描述具有离域共轭键的一般环状分子，这类分子具有增强的热力学稳定性、较小的键长交替、特定的反应活性以及独特的磁性。在本展望中，研究人员聚焦于芳香性的磁性判据。通过计算磁感应电流密度（MICD）磁化率来评估分子芳香性。芳香性程度通过计算磁感应环电流（MIRC）的强度来确定，而芳香性质则由MIRC涡旋的循环方向获得。通过计算不同尺寸的反式环烯烃环的（反）芳香性程度，研究了芳香性和反芳香性的尺寸极限。研究人员展示了如何利用由MICD磁化率计算得到的磁屏蔽密度来理解核磁共振（NMR）屏蔽张量的空间起源。研究人员还考虑了分子结构的拓扑性质及其对MICD的影响。此外，还讨论了MIRC的轨道贡献、群理论方面以及反芳香性环的芳香键合。

# 磁性视角下的分子芳香性研究：《物理化学快报》展望文章解读

## 一、研究背景与科学问题

芳香性是化学中最基础也是最富争议的概念之一。该概念最初用于描述与苯性质相似的化合物，后经扩展用于描述具有离域共轭键的环状分子体系。根据国际纯粹与应用化学联合会（IUPAC）的定义，芳香性的评估主要包括四个判据：能量判据、反应判据、结构判据和磁性判据。然而，这四种判据往往不能给出完全一致的芳香性判断，导致芳香性概念在应用中存在诸多模糊之处。Hoffmann在十年前就曾撰文指出，芳香性概念的过度使用虽然威胁到其学术严谨性，但这一概念本身仍将存续。

从能量判据来看，尽管芳香稳定化能是评估芳香性的最直接指标，但其计算无论在实验还是理论层面都面临困难，尤其对于较大分子体系，找到合适的非芳香参考化合物本身就十分困难。反芳香性稳定化能的测定则更为复杂，因为反芳香性环和芳香性环都比非芳香性环更稳定。反应判据和动力学稳定性与芳香性质的关系通常较为复杂。结构判据虽可应用于大分子和固态材料，但对于键长交替的解读存在局限——反芳香性分子在某些情况下也可能表现出较小的键长交替。重要的是，对于无机团簇来说，结构判据难以有效评估其芳香性。

在众多判据中，核磁共振（NMR）波谱是最为实用的实验手段，其理论基础正是磁性判据：当外部磁场垂直于芳香性分子环平面时，共轭键中的离域电子会维持一个磁感应环电流（MIRC），环外的质子被去屏蔽，环内的质子被屏蔽。然而，MIRC对其他原子NMR化学位移的影响较小且难以测定。正是基于这一背景，本研究从磁性判据出发，系统探讨了芳香性的本质及其量化方法。

## 二、研究内容与主要方法

研究人员的核心工作集中在磁感应电流密度（MICD）磁化率的计算与应用上。他们采用的GIMIC方法建立在包含规范因子的原子轨道（GIAO）基础上，通过单电子密度矩阵和磁微扰密度矩阵来计算MICD磁化率张量，这一方法可适用于从Hartree-Fock（HF）、密度泛函理论（DFT）、二阶M?ller-Plesset微扰理论（MP2）到耦合簇理论（CCSD、CCSD(T)、CCSDT）以及完全活性空间自洽场（CASSCF）等多种电子结构理论水平。样本体系涵盖苯、反式环烯烃系列、circum[n]coronene系列、多并苯及其衍生物、含过渡金属的ScC5H5、Co2O2、堆叠的norcorrole三聚体等。

## 三、主要研究结果

### MIRC强度与芳香性程度的量化

研究人员通过GIMIC方法计算了多种分子的MIRC强度。苯的MIRC强度约为12 nA/T，与理论水平无关；而反芳香性的四氧杂异酞菁在DFT/B3LYP水平下其顺磁环电流强度为-62.5 nA/T，在MP2水平下为-27.4 nA/T。基于经验，当MIRC强度超过3 nA/T的阈值时，可判定为由顺磁环电流主导的芳香性分子。

### 芳香环的尺寸极限

研究人员计算了不同尺寸的反式环烯烃环的MIRC强度。结果表明，在B3LYP水平下，芳香性环的MIRC强度在C30H30处达到最大值48 nA/T，此后随环尺寸增大而下降。C90H90和反芳香性的C92H92的MIRC强度绝对值均小于3 nA/T。在CAM-B3LYP水平下，C50H50和C48H48的MIRC强度绝对值已达到该阈值以下。这说明芳香性存在尺寸极限，且该极限依赖于所使用的泛函。

### 拓扑效应与特殊芳香性

研究人员通过群论分析发现，芳香性本质上是一种拓扑性质。球形芳香性可用2(n+1)2的幻数描述，典型的球形芳香性分子C601??在分子框架内外均呈现顺磁环电流。柱状芳香性适用于具有D∞h对称性的分子。对于M?bius扭曲环，其芳香性规则取决于环绕数（Lk）的奇偶性。此外，Craig型芳香性概念在当代计算化学视角下存在争议，因为原子轨道相位并非可观测量。

### 轨道贡献与应变电流

通过将MIRC分解为σ和π轨道贡献，研究人员发现应变环状分子如1,6-亚甲基环十碳五烯的σ轨道贡献（18.4 nA/T）远大于π轨道贡献（5.6 nA/T），表明应变是其主要MIRC的来源。环丙烯正离子（C3H3?）和环丁二烯（C4H4）的MIRC均由σ和π轨道共同维持。

### 面对面堆叠引起的芳香性质转变

反芳香性分子环可形成具有短键距（约3 ?）的芳香性堆叠二聚体。对于四氧杂异酞菁二聚体，由于前线轨道不可约表示的变化，HOMO-LUMO跃迁从转动允许变为电偶极允许，导致二聚体呈现顺磁环电流，实现了从反芳香性到芳香性的转变。norcorrole三聚体的计算表明，中心环承担顺磁环电流（21.3 nA/T），两侧反芳香性环承担顺磁环电流（-20.3 nA/T），且垂直于堆叠方向也存在净顺磁环电流。

## 四、讨论与结论

研究人员指出，与其他芳香性判据相比（能量判据虽被认为最佳但需参考化合物，结构判据对于三维结构不适用，反应判据的应用范围有限），磁性判据的主要优势在于NMR波谱的强有力实验支持。他们建议MICD研究应使用多个DFT泛函来检验MIRC强度对HF交换百分比的敏感性，并以MP2水平作为可行时的参考。对于含重元素的分子，自旋-轨道耦合效应对MICD的贡献不容忽视，且目前GIMIC尚不能考虑此类效应。他们总结道，MICD磁化率是分子磁性质的基本响应函数，其计算为广泛磁性质的计算研究开辟了新途径。