A Unified Formulation for $\langle \hat{S}^2 \rangle$ in Two-Component TDDFT

本文提出了一种用于在二分量含时密度泛函理论（TDDFT）中计算自旋期望值$\langle \hat{S}^2 \rangle$的统一形式，通过推导针对共线参考态的专用方程，证明了激发态自旋污染既源于参考态，也源于激发过程本身。

要点

想象一下，你正在描述电子的“自旋”。在量子世界中，自旋并非仅仅是一个微小陀螺的旋转；它是一种基本属性，决定了电子在磁场中的行为以及它们彼此之间的成键方式。通常，科学家将电子视为要么“自旋向上”，要么“自旋向下”，就像硬币要么是正面要么是反面。这种方法在简单情况下非常有效。

然而，在复杂分子中或处理重原子时，电子可以做一些更棘手的事情：它们可以同时处于正面和反面的混乱混合状态，或者它们的自旋可以指向奇怪的对角线方向。这被称为“非共线”态。为了处理这种情况，科学家使用一种称为**双分量含时密度泛函理论（TDDFT）**的复杂数学框架。可以将此框架想象为一台高科技相机，能够以三维方式捕捉这些混乱的、对角线的自旋，而不仅仅是平面的二维图像。

问题：“自旋混乱”
当科学家使用这台高科技相机观察激发态（被踢到更高能级的电子）时，他们遇到了一个问题。数学有时会“受到污染”。这就像试图计算罐子里红色和蓝色弹珠的数量，但罐子略微透明，你不小心将一些背景光也算作了弹珠。

在量子力学中，我们有一个想要计算的具体数值，称为$\langle \hat{S}^2 \rangle$（总自旋平方的期望值）。这个数值告诉我们“自旋多重度”——本质上，电子是像一对平静的配对（单重态）那样行为，还是像一群喧闹的未配对群体（三重态）那样行为。如果数学受到污染，这个数值就会出错，使得难以确定实际发生的是何种化学反应。

解决方案：统一的配方
本文作者张晓宇撰写了一种新的“配方”（统一公式），用于正确计算这个自旋数值，无论电子自旋多么混乱。

以下是该论文如何利用简单类比进行分解：

1. 蓝图（二次量子化）：
   作者首先使用一种称为“二次量子化”的语言重写自旋规则。想象电子是舞台上的演员。这种方法不是一次性描述整场戏，而是描述每一位演员的入场和退场。通过这样做，作者表明计算自旋（$\hat{S}^2$）的数学公式与计算能量（$\hat{H}$）的数学公式几乎完全相同。这就像意识到蛋糕的配方只是面包配方的一个略微修改版本。

2. 自旋的两个来源：
   论文发现，激发态的总自旋来自两个不同的地方：

   • 基础自旋（$\langle \hat{S}^2 \rangle_0$）： 这是分子在被激发之前所具有的自旋。它是建筑物的“地基”。
   • 激发变化（$\Delta \langle \hat{S}^2 \rangle$）： 这是电子跃迁到更高能级时额外增加或改变的自旋。它是对建筑物进行的“翻新”。
     论文提供了一种分别计算这两部分然后将它们相加以获得真实总值的方法。

3. “卡希达”机器：
   作者使用一种称为“卡希达方程”（Casida equation）的数学机器（这就像化学中寻找激发态的标准计算器）。通常，这台机器计算能量。作者的大技巧是将机器上的“能量”设置替换为“自旋”设置。由于数学如此相似，这台机器现在可以像吐出能量数值一样轻松地吐出自旋数值。

4. 测试配方：
   为了证明该配方有效，作者在三种不同类型的分子上进行了测试：

   • 水（$H_2O$）： 一种标准的稳定分子。
   • 水离子（$H_2O^+$）： 带电的水分子版本。
   • 氢三重态（$H_3$）： 一种棘手的、不稳定的分子，其中自旋变得非常混乱。

   结果显示，对于简单分子，自旋数值非常干净。但对于混乱的$H_3$分子，该方法正确地识别出自旋是“受污染的”（混合的），这对于试图理解这些分子如何反应的化学家来说是一个至关重要的信息。

为何这很重要
在这篇论文之前，如果你想知道复杂、非直接系统中的激发电子的自旋，你可能不得不根据具体情况使用不同且不一致的方法。这篇论文提供了一本单一的、统一的规则手册，适用于所有这些情况。

这就像拥有一个通用翻译器，可以完美地讲每一种方言，而以前你需要为每个村庄使用不同的翻译器。这使得科学家在研究化学反应、光与物质的相互作用或分子在磁场中的行为时能够更加自信，确保他们不会被数学“噪声”所愚弄。

总结
这篇论文为科学家提供了一种新的、可靠的工具，用于测量复杂系统中激发电子的“自旋”。它将测量分解为“基础”部分和“激发”部分，利用巧妙的数学替换来高效计算，并在各种测试分子上证明了其有效性。它并不承诺直接治愈疾病或制造新电池，但它修复了化学家工具箱中的基本工具，确保他们用于导航量子世界的理论地图是准确的。

技术摘要

以下是 Xiaoyu Zhang 所著论文《两分量含时密度泛函理论中$\langle \hat{S}^2 \rangle$的统一公式》的详细技术总结。

1. 问题陈述

线性响应含时密度泛函理论（TDDFT）是计算激发态的标准工具。然而，在处理非共线磁性和相对论效应的两分量形式中，自旋算符$\hat{S}^2$并不一定守恒，从而导致激发态中出现自旋污染。

• 差距：虽然已有方法可计算自旋守恒 TDDFT 或特定自旋翻转情况下的$\langle \hat{S}^2 \rangle$，但在更广泛的两分量 TDDFT 框架内，尚缺乏评估期望值$\langle \hat{S}^2 \rangle$的统一、通用公式。
• 需求：准确计算$\langle \hat{S}^2 \rangle$对于确定自旋多重度、分析自旋污染以及量化非绝热分子动力学中的自旋混合（例如单重态 - 三重态混合）至关重要。现有方法往往缺乏在一个单一理论框架内涵盖共线和非共线参考态的系统推导。

2. 方法论

作者基于**二次量子化运动方程（EOM）**理论，建立了一个严谨的理论框架。

• 二次量子化表示：
  • 总自旋平方算符$\hat{S}^2$以二次量子化形式表达。作者证明其代数结构与电子哈密顿量$\hat{H}$平行。
  • $\hat{S}^2$被分解为涉及自旋算符$\hat{s}_u$以及产生/湮灭算符的单体和二体项。
• 通用公式推导：
  • 参考态：利用 Slater-Condon 规则推导了一般 Slater 行列式（参考态）的期望值$\langle \hat{S}^2 \rangle_0$，该值用两分量基组中的单粒子约化密度矩阵（$D$）和重叠矩阵（$S$）表示。
  • 激发态：通过将标准 TDDFT 线性响应（Casida）方程中的哈密顿量$\hat{H}$替换为算符$\hat{S}^2$，推导出自旋期望值的变化量$\Delta \langle \hat{S}^2 \rangle$。
  • 类 Casida 本征值方程：作者推导出了一个特定的本征值方程：
    $$\mathbf{M}^s \mathbf{X} = \Delta \langle \hat{S}^2 \rangle \mathbf{N} \mathbf{X}$$
    其中$\mathbf{M}^s$是由自旋算符二阶导数构建的“自旋 - 福克”和“自旋 - 核”矩阵，$\mathbf{X}$代表激发向量。
• 共线参考态的特化：
  • 该通用公式被特化用于共线参考态（源自非限制性 Kohn-Sham (UKS) 或限制性开壳层 Kohn-Sham (ROKS)）。
  • 在此极限下，两分量公式分解为自旋守恒和自旋翻转通道。作者推导了这些特定块简化的工作方程，从而能够与先前文献进行直接比较。
• 非共线核：该框架包含了通过多共线方法推导出的非共线交换 - 相关核，确保理论保留开壳层系统的零能激发（这对于维持自旋态之间的简并性至关重要）。

3. 主要贡献

1. 统一理论：本文提供了两分量 TDDFT 中$\langle \hat{S}^2 \rangle$的首个统一推导，适用于非共线参考态和共线参考态（涵盖自旋守恒和自旋翻转激发）。
2. 自旋污染的分解：分析表明，激发态中的总$\langle \hat{S}^2 \rangle$源于两个不同的来源：
   • $\langle \hat{S}^2 \rangle_0$：参考态固有的自旋污染。
   • $\Delta \langle \hat{S}^2 \rangle$：激发过程专门引入的额外自旋污染。
3. 消除歧义：推导出的共线自旋翻转 TDDFT 方程修正并澄清了先前的公式（特别是解决了 Li 等人以及 Ipatov 等人在度量矩阵$\mathbf{N}$的定义和反对易子处理方面的差异）。
4. 处理零能模式：作者对零本征值模式（由对称性破缺引起的非物理模式）提供了严谨的处理，强制这些物理上平凡的模式满足$\Delta \langle \hat{S}^2 \rangle = 0$，以确保数值稳定性。

4. 结果与数值基准

该方法已在一个基于 JAX 的自研量子化学代码（IQC）中实现，并在三个系统上进行了测试：H$_2$O、H$_2$O$^+$和H$_3$。

• 系统：计算使用了各种泛函（HF、SVWN、PBE、B3LYP）和参考形式（GKS、UKS、ROKS）。
• 观察结果：
  • H$_2$O / H$_2$O$^+$：这些系统表现出可忽略的自旋污染，$\langle \hat{S}^2 \rangle$值接近单重态和双重态的理论极限。
  • H$_3$：这个非共线分子在激发态中显示出严重的自旋污染，归因于基组中自旋空间的不完整性。
  • ROKS 异常：在使用 ROKS 参考态进行的 H$_2$O$^+$计算中，出现了负的激发能。作者将此解释为 ROKS 参考态未能严格维持零能激发（自旋翻转跃迁）的结果，而 UKS/GKS 由于缺乏完整的福克矩阵对角化而保留了该特性。
• 一致性：自旋守恒跃迁的结果与既定的ab initio公式（Myneni 和 Casida）相符，验证了通用推导的正确性。

5. 意义

• 光谱学与动力学：准确计算$\langle \hat{S}^2 \rangle$的能力对于解释光谱数据以及非绝热分子动力学至关重要，因为在后者中，自旋混合驱动了电子跳跃概率。
• 相对论兼容性：虽然该公式是在非相对论极限下提出的，但其设计旨在可迁移至相对论两分量方法（例如 X2C），并可能适用于四分量 TDDFT（用于$\langle \hat{J}^2 \rangle$）。
• 理论基础：通过建立统一的二次量子化框架，这项工作使得系统计算各种可观测量（不仅限于自旋）的期望值成为可能，增强了 TDDFT 对复杂磁性和相对论系统的预测能力。