POET：一种用于解析磁各向异性特定位置电子起源的软件套件

《The Journal of Physical Chemistry A》：POET: A Software Suite for Mapping the Site-Specific Electronic Origins of Magnetic Anisotropy

　　

研究人员提出Palacky光电子工具包（POET），这是一个开放获取的软件套件，可将复杂的模拟数据分解为直观的图形表示，将电子重排映射到原子及轨道贡献上，并阐明成键与磁各向异性之间的相互作用。为展示其效用，研究人员研究了在多种衬底上实现电场可调磁各向异性能（MAE）的两种互补策略。研究表明，在MgO支撑的氮修饰单空位石墨烯上，Pt吸附原子的碘化可产生强烈的、电场可调的平面内各向异性，而OsPt和OsPd异质二聚体则可产生约150 meV的 exceptional 垂直MAE。POET所启发的微观洞见有助于通过化学工程实现对磁各向异性的理性控制，为原子极限下的能效信息存储与处理铺平道路。

本研究旨在开发一种能够解析磁各向异性能（MAE）特定电子起源的计算工具，并应用于设计基于石墨烯的原子尺度磁体。随着全球信息技术的能源需求日益增长，开发能够在原子尺度保留和处理信息的可持续材料成为迫切需求。将磁比特微型化至原子尺度是实现高密度、低能耗信息存储的有力途径，而缺陷石墨烯特别是氮修饰单空位（NSV）石墨烯，为过渡金属（TM）基原子尺度磁体（ASMs）提供了理想宿主，其MAE可通过化学和结构方式在原子精度上进行调控。然而，将这些磁体整合到实际器件中需要与固体衬底相结合，尽管已有大量关于石墨烯在金属和绝缘体表面合成的研究，但衬底尤其是含缺陷石墨烯对这些原子磁体电子和磁学性质的影响尚缺乏深入理解。此外，传统的二阶微扰理论（PT2）虽为解释MAE的电子起源提供了框架，但其固有缺陷在于忽略磁化切换时电子态的重排，且在近简并态占主导时会失效，而完全自洽的非共线（SC NCL）自旋-轨道耦合（SOC）计算结果又过于复杂，难以揭示驱动各向异性的具体轨道机制。因此，研究人员开发了Palacky光电子工具包（POET），以弥合PT2框架的解释清晰度与密度泛函理论（DFT）-SOC结果定量准确性之间的差距，实现对复杂体系中电子结构变化的直观分析。

研究人员采用自旋极化密度泛函理论，利用维也纳从头算模拟软件包（VASP）进行计算。交换-关联效应采用广义梯度近似（GGA）下的Perdew-Burke-Ernzerhof（PBE）泛函描述，芯-价电子相互作用采用投影缀加波（PAW）方法处理。对于标量相对论（SR）和非共线SOC计算，平面波基组的动能截断均设为400 eV。部分占据采用高斯展宽处理，展宽宽度为0.02 eV。电子自洽收敛标准为总能量变化低于10^?6 eV，结构弛豫直至残余力小于10^?2 eV/?。对于卤化系统和含衬底系统，采用DFT-D3修正在内的色散相互作用。布里渊区采样根据各模型系统的横向尺寸进行定制：自由悬挂NSV石墨烯结构弛豫采用6×6×1的Monkhorst-Pack k点网格，支撑系统采用2×2×1网格（矩形MgO(001)表面晶胞为2×4×1）；静态SR计算对所有结构采用6×6×1网格；NCL-SOC计算则需要更密的采样，自由悬挂系统为9×9×1，支撑系统为3×3×1（MgO(001)表面为2×4×1）。MAE通过NCL DFT结合SOC从静态总能量计算中提取，定义为MAE = min(E_x, E_xy, E_y) ? E_z。为确保识别正确的电子和磁学基态，NCL-SOC计算从多个不同的总磁矩初始值出发进行；在标准自洽NCL-SOC弛豫导致非磁性解或弛豫回易磁化方向的情况下，采用约束磁矩方案固定磁矩的方向和/或大小。

研究结果显示，在卤化策略方面，石墨烯框架的直接卤化在能量上是不利的，但锚定在NSV位点的TM原子的卤化则具有强放热特性，且相对于X₂在整个TM@NSV系列中保持能量优势。多卤素吸附（2-3个原子）虽渐趋不利但仍保持放热特征，表明可通过控制卤素暴露程度来调节功能化程度。在衬底效应方面，X-Pt@NSV（X = Cl和I）在Cu(111)、Ni(111)和MgO(001)上的结合强度仅略有降低（最多7%），而在Ir(111)和MgO(111)上则有更强抑制（I-Pt@NSV@Ir(111)高达55%），这一趋势与Pt上的预吸附Bader电荷相关（R=0.77），Pt上正电荷越多，可用于Pt-X键形成的电子密度越少，从而削弱卤素结合。

在TM二聚体方面，TM-TM二聚化通过第二个金属原子（TMT）在预先锚定于NSV位点的TM原子（TMB）上的强放热吸附进行。大多数二聚体采取直立几何构型，仅在MgO(111)上出现倾斜构型，后者因与缺陷边界原子的相互作用而猝灭磁性。OsPd@NSV系统展示了NSV缺陷处的对称性降低如何实现非常规轨道耦合，其高自旋基态（5.00 μ_B）比低自旋构型（1.00 μ_B）低0.46 eV，由自旋向下的σ_s,d杂化（Pd d_z²与Os d_z²及s轨道之间）稳定化。化学取代通过改变价电子数系统性地调控这一成键框架：以Pt取代Pd可减少s-d分离，增加s轨道占据，将自旋向上的π_p,d*能级推至费米能级（E_F）以上，降低反键填充；同时Os d_x²_-y²在自旋向下通道中移至E_F以下，稳定低自旋基态（1.33 μ_B）。

在磁各向异性能方面，自由悬挂NSV石墨烯上的单TM吸附原子通常表现出可忽略的MAE（<1 meV），仅Ir@NSV例外，达到-10 meV且具有平面内易轴。衬底支撑后，大多数TM@NSV系统仍保持有效各向同性；仅Pt@NSV@MgO(111)例外，发展出-25 meV的较大MAE。结构畸变不足以解释此效应，因去除衬底后MAE反转为+8 meV，表明大的负各向异性源于衬底驱动的电子重排而非单纯的结构变形。化学功能化为大各向异性提供了替代途径：Pt@NSV的碘化产生-53 meV的创纪录高MAE，主要源于占据与未占据d_yzd_xz态之间的二阶自旋翻转SOC。外加电场提供了不同的MAE调控机制：与Pt@NSV@MgO(111)最多36%的MAE降低相比，I-Pt@NSV表现出近单调的场依赖性，在-0.5 eV/?时达到-71 meV；而Pt@NSV@MgO(111)的碘化则抑制了MAE，伴随总磁矩的降低。

为实现大的正MAE值，需要改变基本结构单元，从卤化单原子转为锚定在NSV位点的TM二聚体。在自由悬挂系统中，OsPd@NSV具有本研究识别的最大MAE（133 meV），PT2分析将此各向异性归因于E_F附近占据的自旋向上π_p,d*与未占据的自旋向下d_xyd_x²_-y²态之间的自旋翻转耦合。磁力定理（FT）结合POET能带轨迹重建进一步阐明了MAE的轨道起源：磁化垂直于Os-Pd键时，最高占据分子轨道/最低未占据分子轨道（HOMO/LUMO）分别对应π_p,d*和简并的非键d_xy与d_x²_-y²对；沿二聚体键方向排列磁化则解除此简并并重排前线能级，降低总能量。在整个3d-5d二聚体系列中，MAE的大小和符号由前线π_p,d*和d_xyd_x²_-y²轨道的相对能量和占据情况决定。OsPt@NSV和OsPd@NSV因关键耦合通道最大化而 yield >100 meV的MAE。在-1.00至+0.50 eV/?的变化电场下，自由悬挂OsPd@NSV的MAE从62 meV单调增至155 meV；OsPt@NSV的MAE则在相反极性的场下从67单调增至140 meV，且在-0.25至-0.75 eV/?范围内几乎不变，伴随场诱导的低自旋态（~3 μ_B）向高自旋态（~5 μ_B）转变。衬底吸附通常使OsPd@NSV的零场MAE猝灭60 meV以上，仅MgO(111)和MgO(001)表面分别仅降低3和5 meV。OsPt@NSV在MgO(001)上则表现出从120至151 meV的MAE显著增强，源于N原子被Mg²⁺不对称配位导致自旋向上π_p,d*态的简并解除。

在讨论部分，研究人员总结了POET套件的核心价值：通过将磁化依赖的电子重排映射到原子与轨道贡献的直观图形表示，阐明支配磁行为的物理机制。研究得出两个主要结论：一是Os吸附原子在MgO支撑石墨烯的Pd和Pt修饰NSV缺陷处产生exceptionally大的垂直MAE，对应的 predicted 十年阻塞温度约为30-40 K，可与最先进的Ho@MgO单原子磁体相媲美，且远高于直接吸附在石墨烯上的镧系元素；二是单Pt吸附原子的化学功能化提供了更具实验可及性的大各向异性、场可调易轴途径。前者具有极高的 projected 比特密度（440 TBit/in.²）和强磁电响应，后者则因其平面内易轴可实现低耗散自旋输运、超高灵敏度磁力测量和概率自旋逻辑等替代器件范式。