说明:本文华算科技主要介绍高自旋与低自旋的电子占据、配位场竞争、DFT 能量比较及温度和光照下的自旋态变化。

高自旋与低自旋描述同一金属中心、同一 d 电子数下的不同总自旋组态。高自旋保留较多未成对电子,低自旋增加电子配对。单个电子的自旋量子数始终是 1/2;改变的是多个电子的占据与自旋耦合,金属元素名称或形式价态本身并没有唯一指定自旋态。
在理想八面体配位中,五个 d 轨道分为能量较低的 t2g 和较高的 eg 两组,间隔记为 Δoct。电子进入较高轨道要付出晶场能量,在同一轨道内配对也有能量代价。晶场分裂与成对能、交换作用的竞争决定占据:弱场通常有利于高自旋,强场通常有利于低自旋。

以八面体 Fe2+ 的 d6 为例,高自旋占据为 t2g4eg2,有四个未成对电子,S=2,属于五重态;低自旋为 t2g6eg0,电子全部配对,S=0,属于单重态。多重度为 2S+1。首图上部给出两种占据,下部的封装与显微图描述这些分子的器件环境。
“低自旋”并不要求所有电子配对。八面体 d5 的高自旋 S=5/2,低自旋仍有一个未成对电子,S=1/2;d7 则可由 S=3/2 变为 S=1/2。低自旋也可以具有局域磁矩。八面体 d4—d7 最常讨论这类竞争;四面体、平方平面或明显畸变的配位应按各自轨道分裂重新排列电子。
这些 dn 占据是离子模型中的记号。金属与配体形成共价键后,金属 d 轨道的投影电子数可以偏离整数,自旋密度也可延伸到配体。判断真实材料中的高低自旋,应检查轨道占据、局域自旋密度和配位构型,而不能把一个原子投影数值直接当成整数电子构型。

反键轨道占据怎样影响键长?
八面体中的 eg 轨道沿金属—配体方向分布,常与配体 σ 轨道形成明显的反键组分。Fe2+ 从低自旋转为高自旋时,eg 相关反键态增加占据,金属—配体键常随之伸长。结构弛豫又会改变轨道重叠和晶场分裂,电子占据与键长相互影响。
三联吡啶铁配合物的光激发实验中,Fe—N 距离从约 2.0 Å 增至约 2.2 Å,配体也发生位移和转动。键伸长与配体整体重排具有不同坐标,图中的四组曲线分别记录它们的演化。这里的约 0.2 Å 属于该配合物的结构变化,并非所有高低自旋转变都具有相同键长差。

配体的 σ 给电子和 π 相互作用会改变两组 d 轨道的相对能量。π 接受配体可降低金属 t2g 相关态的能量,从而增大分裂;π 给电子配体可能产生相反趋势。金属价态、键长和配位畸变共同参与晶场变化,单按配体名称给所有材料指定高低自旋,会遗漏晶体或表面上的实际构型。
限域和电场怎样改变能量差?
分子进入碳纳米管后,朝向、分子—管壁接触及轨道杂化都会改变两种自旋态的能量。该 Fe 配合物模型中,低自旋态与配体、管壁的轨道混合较强,限域提高了高自旋相对低自旋的能量。自旋态排序具有具体环境条件,孤立分子的结果与封装体系应分别比较。

同一模型的电场计算还显示,沿不同方向施加电场,EHS−ELS 可向不同方向变化。电场方向与分子取向共同影响两态的能量响应,原因包括电荷分布、偶极矩和极化响应的差异。图内场强采用原文单位;将其解释为器件工作电压,需要电极几何与实际电势分布。

分别弛豫后,比较的是什么能量?
在相同原子组成、电荷、泛函和数值精度下,分别优化高低自旋结构,常用的绝热电子能差为 ΔEad=EHS(RHS)−ELS(RLS)。正值表示低自旋的电子能较低。若把两态都放在同一结构 R 上比较,得到的是垂直能量差,其中尚未包含另一自旋态的结构弛豫。
多金属模型中还需指定每个中心的自旋。三唑桥联铁链的七核模型,只有中间五个 Fe 位点切换高低自旋,末端两个位点始终固定为高自旋。全低自旋、全高自旋和高低自旋混合的五位点组态,具有不同键长与弹性应变;比较结果应对应这些具体组态。

初始磁矩只参与自洽迭代的起点选择。优化结束后,应核对收敛态的局域磁矩与占据矩阵,确认两次计算确实保留了两种电子组态。基态搜索还应比较配位模型允许的其他自旋组态。金属—配体杂化可能降低原子投影磁矩;反铁磁排列则可使晶胞总磁矩为零,而各金属中心仍保持高自旋。局域组态与位点间磁序要分别记录。
泛函差异怎样影响自旋态排序?
杂化泛函中的精确交换比例改变金属 d 电子的交换贡献,DFT+U 修正其局域占据,两者都会影响高低自旋的相对电子能。不同方法可能给出不同的基态排序。固定同一套几何比较方法差异,可以辨认电子能贡献;重新优化后还会叠加结构变化。某个 U 值或泛函的结论应由相应材料的结构和谱学数据检验。

铁链模型中,改变绝热电子能参数后,中间五个铁位点从全低自旋变为全高自旋的 ΔG 曲线整体移动,与零线相交的温度也改变。电子能误差会传递到预测转变温度。这张图比较的是能量参数敏感程度,参数按单个可切换铁中心计量,分别取 10、15、20、25 kJ mol−1,并非四种泛函的直接对照。

热平衡与光激发各自比较什么?
高低自旋转变的热平衡由焓差与温度加权的熵差共同决定,按 ΔG=ΔH−TΔS 比较两态自由能。高自旋通常具有较多自旋简并态;在孤立、自旋简并未分裂的 Fe2+ 模型中,自旋熵差为 R ln 5。较长的键还常带来较低振动频率和较大的振动熵。热膨胀、压力下的体积项及晶格弹性相互作用,也会改变自由能。
光照则可以从低自旋基态激发到金属向配体电荷转移态,再经系间窜越形成高自旋激发态。三联吡啶铁的占据示意与电子衍射装置描述这一过程:光脉冲产生激发,电子脉冲记录结构变化。图中的高自旋布居属于非平衡激发态,寿命受电子弛豫、分子振动和周围晶格影响。

三联吡啶铁分子的 Fe—N 键伸长会改变分子在晶胞内的几何匹配,邻近分子对配体进一步重排形成约束。链状材料中的相邻中心发生键长变化,也会使高低自旋混合组态产生弹性应变。升温和降温时出现的热滞回涉及协同作用与亚稳态;孤立位点的零温电子能差,尚未包含晶体长程相互作用、自由能势垒及扫描过程。
磁响应与热响应怎样区分?
共线计算常报告占据差对应的自旋磁矩;在自旋贡献占主导的孤立离子近似中,饱和磁矩约为 2SμB,顺磁有效磁矩则约为 2√[S(S+1)]μB。对 S=2,两者分别约为 4 和 4.90 μB。DFT 磁矩与变温磁化率采用不同数值定义;共价杂化、轨道贡献和磁交换会进一步改变测量值。

该三唑铁链在 10 K min−1 扫描时,升温约 387 K 出现吸热峰,降温约 343 K 出现放热峰,三次循环保留了热滞回。磁化率记录自旋响应,DSC 记录转变热量,两种测量分别检验自旋组态与热力学变化。该材料的转变焓和熵分析还使用了穆斯堡尔谱确定的位点转换比例。
