《Journal of Chemical Theory and Computation》：Spin Dynamics in Open Quantum Systems: A DLvN-TDDFT Approach
编辑推荐：
研究人员提出了一种在含时密度泛函理论（TDDFT）框架内发展的自旋不平衡驱动Liouville–von Neumann（DLvN）方法，用于模拟开放量子系统中的共线电子自旋输运。在引入并验证该方法后，通过简单分子结中自旋极化输运的基准模拟，研究人员将其应用于外
研究人员提出了一种在含时密度泛函理论（TDDFT）框架内发展的自旋不平衡驱动Liouville–von Neumann（DLvN）方法，用于模拟开放量子系统中的共线电子自旋输运。在引入并验证该方法后，通过简单分子结中自旋极化输运的基准模拟，研究人员将其应用于外部电场下的磁性锯齿形石墨烯纳米带结模型。模拟揭示了丰富的自旋分辨电流动力学特性，凸显了DLvN-TDDFT框架作为探索低维开放量子系统中动态自旋电子学现象的有前途工具。
自旋电子学作为利用电子自旋自由度进行信息处理的前沿领域，其核心目标是在纳米尺度器件中实现对自旋流的精确调控。自Aviram和Ratner提出分子电子整流器概念以来，分子电子学取得了长足发展，为电子器件小型化、功能集成化提供了新路径。其中，分子干涉仪和手性诱导自旋选择性（CISS）效应等实验进展，展示了在分子结中操控自旋自由度的潜力，为实现基于自旋的新型功能器件奠定了基础。然而，实验手段在揭示单分子结输运背后的复杂物理机制方面存在局限，难以满足设计与控制新型动态自旋器件的需求。基于第一性原理的准确电子动力学模拟，为此提供了关键的理论研究平台。
驱动Liouville–von Neumann（DLvN）方法是近十年来发展的一种概念简洁、计算高效的电子动力学模拟方案，特别适用于研究分子结的瞬态响应与外部时域调控效应。此前，该方法已在紧束缚和扩展休克尔等简单电子结构模型下得到验证，并初步扩展至TDDFT框架，实现了从第一性原理出发对真实分子结模型电子输运的动态研究。然而，早期实现仅限于自旋平衡体系，无法描述磁化动力学和自旋输运过程。为突破这一局限，本研究正式引入了共线自旋不平衡DLvN-TDDFT模拟方案的形式体系、实现方法及其实际应用。
该方法的实施包含四个核心步骤。首先，对模型体系进行空间划分，通常分为左电极（L）、扩展分子区（EM）和右电极（R），并定义相应的自旋分辨Kohn–Sham（KS）哈密顿量矩阵和重叠矩阵。其次，为隔离边界条件对分子区动力学的影响，对扩展分子区的原子中心基函数进行块正交化处理，通过非幺正变换使扩展分子区与电极区的基函数正交。接着，进行从实空间（站点）表示到能谱（状态）表示的变换，以便在状态表示下施加开放边界条件。最后，在状态表示中求解自旋不平衡的DLvN运动方程，该方程包含两部分：第一项描述基于Liouville–von Neumann方程的幺正动力学演化；第二项则通过引入驱动项，模拟电极作为隐式热库向系统注入平衡电子并抽取流出电子的开放边界条件，其中驱动速率Γσ表征了电极的弛豫速率。完成边界条件施加后，驱动项需变换回实空间表示，以在实空间进行密度矩阵的时间演化。整个流程通过自洽迭代实现：基于当前密度矩阵构建KS哈密顿量，然后进行上述变换与演化步骤。自旋分辨电流则可通过变换后的密度矩阵块与耦合矩阵块计算获得。
在实现层面，研究团队开发了一个Python驱动程序，调用Gaussian程序包在每个时间步基于瞬时密度矩阵构建KS哈密顿量。计算成本主要来源于KS哈密顿量的构建以及与位点-状态变换相关的线性代数运算，当前实现的算法复杂度为O(N3)。
研究首先在两个简单模型体系上对自旋不平衡DLvN-TDDFT方法进行了验证。第一个体系是由氢原子链电极桥接的氢分子（H–H）结，其电子结构为自旋平衡。模拟显示，在施加对称偏压时，自旋不平衡方案得到的两个自旋通道电流与自旋平衡方案的结果完全一致，验证了方法在自旋平衡极限下的自洽性。第二个体系是采用垂直取向H–He–H线性分子作为桥的模型结，该分子具有简并的反铁磁和铁磁基态，是典型的自旋不平衡桥模型。研究表明，即使施加对称偏压，该体系的电流响应也强烈依赖于初始自旋构型。当从反铁磁构型开始时，两个自旋通道的电流大小相等但流经不同的空间路径：↑自旋电子主要流经具有↑自旋极化的下方氢原子路径，而↓自旋电子则流经具有↓自旋极化的上方氢原子路径。相反，从铁磁构型开始时，体系表现出净自旋电流，因为↑自旋电子可通过中心非极化路径和两个极化侧路径传输，而↓自旋电子仅能通过中心路径传输，导致传输通道数目不平衡。进一步地，研究人员演示了施加自旋相关偏压的情形。对于铁磁构型，当仅对↑自旋通道施加偏压时，该通道达到稳态电流，而未偏压的↓自旋通道则表现出显著的瞬态上游电流，随后衰减至零，这种非平衡现象值得未来深入研究。
在验证方法有效性后，研究人员将其应用于一个更实际的石墨烯纳米带（GNR）结模型。该模型由两个扶手椅型GNR电极桥接一段锯齿形GNR构成，其自旋极化特性可通过外加横向电场进行调控。模拟结果表明，即使施加对称偏压，↓自旋通道在瞬态阶段也表现出显著的上游电流，几乎镜像了↑自旋通道的下游动力学。这归因于锯齿形GNR段在初始反铁磁构型下，其自旋极化向电极延伸，在结入口和出口分别形成了对↓自旋和↑自旋电子的泡利传导势垒，导致了电子反射与电流振荡。最终，在系统的静电响应下，两个自旋通道的电流趋于一致的下游稳态。当施加垂直于锯齿形GNR轴向的面内电场时，两个自旋通道的瞬态电流镜像对称性被破坏，稳态电流发生劈裂。由于电场降低了↑自旋的能隙并增大了↓自旋的能隙，↑自旋电流在瞬态和近稳态均高于↓自旋电流，这直接演示了利用该方法模拟外部扰动对纳米尺度自旋电子器件动态自旋过滤效应的调控能力。
综上所述，本研究提出了基于TDDFT框架的共线自旋不平衡DLvN方法，实现了对开放量子系统中自旋输运的实时空间分辨模拟。通过在简单分子结和实际石墨烯纳米带结上的应用，该方法成功捕捉了依赖于自旋构型、偏压方案及外电场调控的丰富自旋分辨电流动力学行为。这些结果确立了DLvN-TDDFT框架作为探索开放磁性量子体系自旋动力学的稳健而多功能的平台，为分子自旋电子学和量子输运工程领域未来通过外部时域控制进行研究铺平了道路。
