环碳[6]中的高自旋态诱导芳香性反转

《ACS Omega》：High-Spin-State-Induced Aromaticity Reversal in Cyclo[6]carbon

　　

环碳[6]（Cyclo[6]carbon，C6）是最小的偶成员环碳，为考察高度应变的sp杂化碳环如何响应电子态变化提供了严格平台。研究人员采用密度泛函理论（density functional theory，DFT）、含时密度泛函理论（time-dependent density functional theory，TD-DFT）、从头算分子动力学（ab initio molecular dynamics，AIMD）以及CASSCF(12,12)占据数分析，考察了C6的结构、电子、磁响应、激发态及稳定性/反应性行为。优化所得S0态结构严格平面，呈现显著键角交替而C–C键长近乎均匀，表明其成键无法用简单定域交替键图像描述。实空间描述符IRI-π与ELF-π勾勒出S0中连续双π离域骨架，并为分析态依赖磁响应提供电子参照。基于ACID、ICSSzz与NICSzz的磁响应分析进一步显示，所选高自旋五重态Q1由S0的抗磁（diatropic）响应逆转为顺磁（paratropic）响应。CASSCF(12,12)占据数分析支持S0 C6主要为闭壳层描述，无显著开壳层多参考特征，从而支持基于DFT描述符的定性判读。态密度与代表性激发态分析进一步表明电子态占据变化使电子密度在πin与πout流形内重新分布，为双π骨架内的态依赖变化提供额外电子结构背景。AIMD模拟指示S0环骨架在低温下动力学稳定而高温下柔性增强，ESP与ALIE分析确认C–C骨架为主要反应区域。上述结果揭示了高自旋Q1态中显著的芳香性反转，并凸显C6双π骨架的自旋态依赖磁响应。

研究背景方面，环碳仅由sp杂化碳原子排成环，几何极简却呈现非经典成键与离域，其价电子可形成两个相互正交的π电子系统：一个近似位于分子平面内（π_in），一个垂直于分子平面（π_out），由此引出双芳香性概念。自环碳[18]直接成像以来，不同环尺寸偶成员环碳的合成与理论阐释受到关注，表面合成研究进一步显示环尺寸微调会显著改变结构与磁芳香性特征。近期Sun等人报道了线型与环型C₆的表面合成与表征，使C₆从假想模型变为实验可达极限成员。已有研究表明环尺寸强烈调制电子离域与芳香响应，且C₆作为最小偶成员环碳须在极端紧凑骨架中承受角应变同时保持电子相干，因此单一几何指标不足以指派芳香性，需结构、实空间电子拓扑与磁响应多方证据一致。在此背景下，考察极小sp碳环是否仍能经双正交通道维持全局π离域，以及电子态占据变化能否逆转芳香磁响应，成为核心动机；近期环碳[16]中经振动效应发现的芳香性反转进一步激励了对最小成员C₆自旋态依赖磁响应切换的聚焦。该研究发表于《ACS Omega》，意义在于以最小环碳为极限体系揭示自旋态对双π骨架磁芳香性的选择性逆转机制。

主要关键技术方法方面，研究人员以ωB97X-D/def2-TZVP级别进行气相DFT几何优化与频率计算确认S₀为真极小，并以同级别TD-DFT计算吸收谱与激发态；在S₀优化几何上做CASSCF(12,12)单点计算（12电子落入12个前沿轨道对应双π流形）作多参考诊断；用Multiwfn与AICD 2.0做IRI-π、ELF-π、LOL-π实空间分析以及ACID、ICSS_zz、NICS_zz磁响应分析；用CP2K/Quickstep在NVT系综（CSVR thermostat）下于4.7、10、20、50、100 K做5 ps AIMD（PBE+D3(BJ)，DZVP-MOLOPT-SR-GTH-q4基组，GTH-PBE赝势）评估有限温度结构稳定性；以ESP与ALIE分析反应区域。样本为气相中性C₆，无外部样本队列。

研究结果部分，2.1节分子结构表明优化S₀严格平面，内键角约90°与150°交替，所有C–C键长约1.32 ?近乎一致，键角交替与键长均匀解耦，说明应变经方向性几何变形而非键长分化表达，电子离域维持相对均匀成键框架。2.2节基态双π电子离域中，IRI-π填色图与等值面显示π相关作用区沿环连通分布，未被键角交替中断；补充LOL-π与ELF-π支持两正交π流形共存且π_out共轭略强，构成S₀电子参照。2.3节态依赖磁响应与芳香性反转中，ACID显示S₀具双正交π回路抗磁环电流，高自旋五重态Q₁（多重度5）逆转为顺磁环电流；ICSS_zz与NICS_zz扫描一致显示S₀抗磁屏蔽分布、Q₁图案反转；三重态T₁–T₃保留双通道贡献但未现Q₁级逆转，表明磁响应改变具态选择性。2.4节态依赖磁响应电子结构背景中，TDOS/PDOS显示σ居低能、前沿由π_in/π_out主导，HOMO–LUMO隙8.01 eV（?10.61与?2.60 eV）；CASSCF(12,12)占据数居1.901–1.966（占轨）与0.041–0.143（虚轨），近理想闭壳层，无显著开壳层多参考特征。2.5节双π骨架内激发态电子重排中，吸收谱50–450 nm有多离散峰，98.44 nm主峰来自简并S₅₁/S₅₂；TDM与CDD显示代表激发S₃、S₅₁、S₉₈偏π_in，S₂₁偏π_out，均为π→π*局域伴部分电荷转移，电子重排发生于预存双π骨架内而非新成键框架。2.6节热稳定性与反应性中，AIMD显示4.7 K下RMSD小、环骨架紧凑，升温至100 K波动增强；ESP正面区表面积大于负面区，ALIE低值沿C–C键分布，确认C–C骨架为主电子活性与亲电作用区。

讨论与结论翻译部分，研究人员在结论中指出：C₆ S₀平面、键角显著交替而键长近均匀，不能简化为定域交替键模型；IRI-π、ELF-π等支持两正交π流形共存，构成磁响应电子框架。磁响应分析确认Q₁为芳香性反转关键态：S₀双π骨架抗磁环电流，高自旋五重态Q₁环电流反向呈顺磁响应；T₁–T₃未现同等逆转，反转具态选择性且在Q₁最显著。CASSCF(12,12)占据数支持S₀主闭壳层，DOS/PDOS与激发态分析显示电子重排主要在π_in/π_out内发生，与Q₁反转源于预存双π骨架内电子态占据改变而非新成键骨架一致。AIMD示低温保环、高温波动增大；ESP与ALIE示C–C骨架为主电子活性区。整体表明C₆高自旋Q₁相对抗磁基态发生显著芳香性反转，为该极小碳环自旋态依赖磁芳香性提供态分辨图像。