《The Journal of Physical Chemistry Letters》：Real-Time CASSCF Ehrenfest (TD-CASSCF) Modeling of Spintronics in a Tetrathiafulvalene-Based Nitronyl Nitroxide (BTBN) Model: The Competition between Charge Transfer and Spin Permutation
编辑推荐：
利用实时全电子TD-CASSCF/Ehrenfest模拟，研究人员研究了有机材料BTBN（一种四硫富瓦烯-氮氧自由基双自由基）中磁电阻的分子起源。计算表明，电荷输运与自旋动力学存在强耦合。在高自旋态下，局域自旋平行排列，电荷通过相干的波状迁移在π堆叠体系中高效
利用实时全电子TD-CASSCF/Ehrenfest模拟，研究人员研究了有机材料BTBN（一种四硫富瓦烯-氮氧自由基双自由基）中磁电阻的分子起源。计算表明，电荷输运与自旋动力学存在强耦合。在高自旋态下，局域自旋平行排列，电荷通过相干的波状迁移在π堆叠体系中高效传播。相反，在含有反平行自旋排列的低自旋态中，电荷迁移与局域自由基自旋之间的自旋置换过程相互竞争。这些额外的自旋交换途径不传递电荷，但会将布居从电荷输运通道中分流，在所研究的示例中，相对于高自旋态，电荷转移效率降低了约50%。因此，模拟确定了电子迁移与自旋置换之间的竞争是BTBN中磁电阻的基本分子机制，从而解释了磁场诱导的自旋排列为何会增强电荷输运并降低电阻。
**研究背景与问题**
磁电阻（MR）现象是指材料电阻随外加磁场变化的现象，是自旋电子学领域的核心效应之一。BTBN是一种由四硫富瓦烯（TTF）与氮氧自由基（NN）单元共价偶联而成的纯有机双自由基材料，在空穴注入后表现出磁电阻行为。现有对分子材料中磁电阻的机理解释通常依赖预先设定的物理假设，例如交换相互作用、自旋离域等概念，这些方法限制了真实机制的描述。为突破此局限，研究人员采用实时全电子动力学方法，该方法不施加先验约束、不引入自由参数，能够直接提供真实动力学的分子图景。研究目的在于揭示BTBN中磁电阻的微观分子机制，特别是电荷输运与局域自旋动力学之间的竞争关系。
**研究内容与结论**
研究人员利用实时全电子TD-CASSCF/Ehrenfest模拟，构建了包含四个BTBN单元的π堆叠模型体系，该模型源自实验晶体结构，并进行了适当的简化处理（以F替代Br，以H替代长链烷基）。通过求解含时薛定谔方程，分别模拟了高自旋（六重态）与低自旋（四重态）初始条件下的电子与自旋动力学过程。计算结果显示，在高自旋态下，电荷以波状相干方式在TTF单元间高效迁移；而在低自旋态下，电荷迁移过程与NN局域自旋之间的自旋置换过程相互竞争，导致电荷转移效率下降约50%。这一竞争机制被确认为BTBN磁电阻的分子起源，解释了磁场为何通过稳定自旋平行排列态来增强电荷输运。该研究发表在《The Journal of Physical Chemistry Letters》。
**主要技术方法**
研究人员采用实时CASSCF（Complete Active Space Self-Consistent Field）方法，在Ehrenfest动力学框架下传播电子波函数。模拟使用CASSCF(11,8)/6-31G*级别，活性空间包含4个TTF最高占据分子轨道（HOMO）和4个NN单占据分子轨道（SOMO），共11个电子。电子与自旋动力学通过含时波函数在局域分子轨道构建的组态态函数（CSF）上投影的布居数演化进行分析。轨道局域化采用Boys算法，电子耦合矩阵元由CASSCF组态相互作用哈密顿量给出。核位置在10飞秒动力学过程中保持冻结，忽略超快时间尺度上的结构弛豫效应。
**研究结果**
**1. 高自旋态下的电荷迁移**
模拟从六重态初始构型出发，空穴定域于顶部BTBN-1单元的TTF部分。随时间演化，布居以波状方式从BTBN-1迁移至底部BTBN-4单元，表明自旋平行排列时电荷可以通过相干波状迁移实现高效输运。
**2. 低自旋态下的电荷迁移与自旋置换竞争**
在四重态初始构型中，相邻NN单元上存在反平行自旋排列。虽然电荷仍以波状方式迁移，但顶部BTBN-1单元布居的衰减速度相比六重态减慢约50%，表明电荷输运效率下降。原因在于低自旋态下出现了自旋置换过程，即局域NN自旋之间发生翻转而不伴随电荷转移。模拟中观察到处在NN1与NN3单元之间的自旋置换具有较大的电子哈密顿量矩阵元，而其他置换通道的贡献小一个数量级。
**3. 自旋置换的机制分析**
电荷迁移与自旋置换的速率由电子哈密顿量矩阵元K
ij控制。该矩阵元包含单电子轨道耦合项、库仑积分项和双电子交换积分项。电荷迁移主要受单电子项支配，而自旋置换仅涉及双电子交换项。由于自旋置换不产生净电荷转移，其占用了电荷迁移的竞争通道，从而提高了体系的有效电阻。
**讨论与结论**
模拟结果表明，BTBN中磁电阻的分子机制源于电荷迁移与自旋置换之间的竞争。当NN局域自旋平行排列时，TTF迁移电子无需改变自旋方向即可在相邻单元间转移；而当存在反平行自旋排列时，自旋置换成为独立的竞争通道，耗散了部分电荷迁移的布居流，降低了电荷输运效率。应用磁场使自旋趋于平行排列，即可关闭自旋置换通道，恢复高效电荷输运，从而降低材料电阻。这一发现为分子尺度磁电阻现象提供了直接的量子动力学证据，并且该模拟方法不依赖任何预设物理模型，具有普遍的适用性，为设计新型有机自旋电子材料提供了理论指导。
