1. 研究概要
该研究由来自德累斯顿工业大学（Technische Universität Dresden）化学与食品化学学院的Hongde Yu和Thomas Heine共同完成。通讯作者为Thomas Heine，其同时隶属于亥姆霍兹德累斯顿罗森多夫研究中心（Helmholtz Zentrum Dresden-Rossendorf）的资源生态研究所和韩国延世大学（Yonsei University）化学系。这项原创性研究成果于2023年8月23日发表在《美国化学会志》（Journal of the American Chemical Society）第145卷，第19303-19311页。
2. 研究背景与目标
本研究属于分子自旋电子学领域，聚焦于无金属磁性材料的前沿探索。长久以来，实现纯粹由有机分子构成、不依赖过渡金属元素的磁性材料一直是该领域的重大挑战。三角烯（triangulene）作为一种具有本征自旋极化的非凯库勒（non-Kekulé）多环芳烃，其最小的纳米石墨烯片段拥有三重态基态，被视为构建扩展磁性结构的理想单体。然而，一个长期存在的难题是，已合成的三角烯二聚体（dimers）仅表现出弱的反铁磁（antiferromagnetic, AFM）耦合，其磁交换耦合强度（J）远低于室温应用所需的朗道尔极限（Landauer limit, 18 meV）。例如，直接相连的二聚体J值仅为-14 meV，苯基桥连的二聚体更是只有-2 meV。
为突破这一瓶颈，研究者设定了清晰的目标：系统性地探索并增强三角烯二聚体中的磁耦合作用，并尝试实现铁磁（ferromagnetic, FM）耦合。研究旨在通过第一性原理计算，不仅寻找超越朗道尔极限的强反铁磁体系，还要挑战传统的奥夫钦尼科夫规则（Ovchinnikov rule），实现对磁耦合类型和强度的有效调控，为无金属有机自旋电子学器件的设计提供分子层面的理论指导和实用策略。
3. 研究流程与方法
本研究的工作流程由四个核心步骤构成：体系构建、结构优化、磁性与电子结构计算、以及数据关联分析与机理探索。
步骤一：模型体系构建 研究团队选取了六种具有代表性的三角烯类似物作为单体构建模块。这包括原型三角烯（tri），以及氮掺杂（tri(n)）、硼掺杂（tri(b)）、三芳基甲基（tam）、三氧三角烯（tot）和酚烯基（ply）衍生物。这些单体在对称性和电子结构上具有多样性，例如tri具有D3h对称性，而tri(n)因姜-泰勒畸变（Jahn-Teller distortion）变为C2v对称性。利用这些单体，他们构建了24种二聚体，其连接方式包括直接相连，以及通过乙炔基（-C≡C-）、丁二炔基（-C≡C-C≡C-）和苯基（phenyl）作为桥连基团。
步骤二：结构优化 所有二聚体分子均首先作为自由分子，使用Gaussian 16软件包在PBE0/def2-TZVP理论水平下进行几何结构优化。优化结果显示，由于氢原子的空间位阻效应，直接相连的二聚体均呈现扭曲构象，而非平面结构，例如tri-tri二聚体的二面角为37°。
步骤三：多态能量与磁耦合计算 为评估磁耦合强度，研究者基于海森堡-狄拉克-范弗莱克（Heisenberg-Dirac-Van Vleck, HDVV）哈密顿量（Hamiltonian），采用破缺对称密度泛函理论（broken-symmetry DFT）计算了二聚体在低自旋态（即开壳层单重态，open-shell singlet, OSS）和高自旋态下的能量差，以此量化磁耦合常数J。其中，J<0表示反铁磁相互作用，J>0表示铁磁相互作用。同时，他们还计算了关键的闭壳层单重态（closed-shell singlet, CSS）能量，用于提取哈伯德模型（Hubbard model）中的现场库仑排斥能（U）。
步骤四：电子结构分析与机理探索 通过对分子前线轨道（尤其是单占据分子轨道SOMO）的分析，研究者量化了单体片段间的电子耦合作用（即跃迁积分t），该值与最高占据分子轨道（HOMO）与最低非占据分子轨道（LUMO）的能隙（Eg）成正比，满足关系式Eg = 2|t|。依据古迪纳夫-金森规则（Goodenough–Kanamori rule），磁耦合J被分解为铁磁性的势交换（potential exchange, K）和反铁磁性的动能交换（kinetic exchange, -4t²/U）两部分，从而阐明了各体系的磁性起源。为进一步调控磁性，他们通过改变二面角（φ）模拟了从扭曲到平面的构型变化，并使用CASSCF/NEVPT2高精度多参考方法对关键体系进行验证。此外，他们还通过引入甲基（-CH₃）作为空间位阻控制单元（SCU）来锁定特定的二面角。
4. 主要研究结果与分析
磁耦合的广泛分布与强耦合体系发现 计算结果显示，24种二聚体的磁耦合J值跨越了从-0.07 eV至-144 meV的巨大范围，且全部为OSS基态，表明基态下均为反铁磁相互作用。其中，tam-tam二聚体表现出惊人的强反铁磁耦合，J值高达-144 meV，是原型tri-tri二聚体（-9.7 meV）的近15倍。此外，tot-tot（-60 meV）、tot-cc-tot（-54 meV）和ply-ply（-20 meV）等多个体系的J值也显著超越了室温下能量耗散最低阈值的朗道尔极限（18 meV）。这一发现证明，通过合理的分子设计，获得具有强健反铁磁耦合、可在室温工作的无金属磁性材料是完全可行的。
结构与磁性关联及内在机理 对电子结构的分析揭示了这些现象背后的深刻物理图像。所有二聚体的SOMO本质上都是其单体片段前线轨道的叠加。tam-tam拥有D3h高对称性的SOMO，其单体间的重叠积分（0.448）和电子耦合跃迁积分t（0.56 eV）远大于C2v对称性的tri-tri（重叠积分0.025，t=0.14 eV）。这种增强的电子耦合直接导致了动能交换项（-4t²/U）的急剧增加，从而极大地强化了反铁磁相互作用。研究者进一步建立了一个明确的正相关关系：在能隙大于0.5 eV的体系中，HOMO-LUMO能隙越大，电子耦合越强，反铁磁耦合J值就越大。这为设计强磁性分子提供了一个清晰的判据：寻找具有大能隙和高对称性的非凯库勒多自由基。
平面化策略对磁耦合的显著增强 研究发现，通过从自由分子的扭曲构型转变为平面构型，可以进一步增强反铁磁耦合。例如，tri-tri的J值从-9.7 meV增至-20 meV，这与先前表面实验观测到的-14 meV吻合。tam-tam的J值更是在平面构型下达到惊人的-198 meV。这种增强归因于平面化后更大的SOMO重叠和电子耦合。
氮掺杂二聚体中反铁磁到铁磁的颠覆性转变 最具突破性的发现来自tri(n)-tri(n)二聚体。随着二面角的变化，该体系的磁耦合符号发生了翻转。当二面角在约30°至70°之间时，表现为反铁磁耦合（J<0），而在小于30°（包括平面结构）或大于70°时，则转变为铁磁耦合（J>0）。在平面构型下，其J值为+1.2 meV，而更高精度的CASSCF/NEVPT2计算给出了+23.6 meV的显著铁磁耦合值。这一现象是由于势交换（K）和动能交换（-4t²/U）之间的精妙平衡被打破所致。这一发现打破了奥夫钦尼科夫规则的限制，代表了首个无法被该规则预测的三角烯二聚体铁磁耦合实例，为实现无金属铁磁性开辟了全新道路。通过引入甲基作为空间控制单元，他们成功将二面角锁定在80°，也确实观测到了J=0.3 meV的铁磁相互作用。
5. 研究结论与价值
本研究通过系统性的理论计算，成功揭示了24种三角烯二聚体中磁耦合的可调性，实现了高达-144 meV的强反铁磁耦合。研究不仅在分子层面阐明了能隙、电子耦合与磁相互作用的定量关系，更重要的是，通过巧妙的分子扭曲策略，首次在氮掺杂三角烯二聚体中实现了从反铁磁到铁磁的调控，获得了无法被传统理论预测的铁磁耦合。这一工作为增强和调控无金属磁体中的磁相互作用提供了切实可行的分子设计策略，对推动有机自旋电子学、量子计算和数据存储等领域的发展具有重要的科学价值和深远的应用前景。
6. 研究亮点
本研究的亮点可归纳为以下三点： * 性能突破：发现了J值高达-144 meV的tam-tam二聚体，远超室温应用阈值，为实际应用奠定了基础。 * 机理创新：建立了能隙、电子耦合与反铁磁相互之间的正相关性，提供了一个简洁而强大的分子设计指导原则。 * 概念革新：在tri(n)-tri(n)中首次实现了对奥夫钦尼科夫规则的超越，通过调节二面角实现了反铁磁到铁磁的转变，创造了首例铁磁三角烯二聚体，颠覆了传统认知。