含时密度泛函理论

玻尔百科

核心要点

含时密度泛函理论 (TD-DFT) 将 DFT 扩展至模拟电子对含时场的响应，从而可以计算电子激发能。
由于常见交换相关泛函的局域性，TD-DFT 中的标准近似在处理长程电荷转移激发和具有双重激发特性的态时会系统性地失效。
范围分离杂化 (RSH) 泛函通过修正长程交换相互作用，极大地提高了对电荷转移态的计算精度，提供了一个强大的解决方案。
TD-DFT 是多个领域中的关键工具，能够用于设计 OLED 材料、通过 QM/MM 方法研究视觉等生物过程以及解读 X 射线光谱。

引言

理解光与物质的相互作用是几乎所有现代科学分支的基础，从自然界绚丽的色彩到驱动我们屏幕的先进技术。虽然密度泛函理论 (DFT) 为描述分子和材料在其最低能量状态（即基态）下的电子结构提供了一个卓越的框架，但当面对由光引发的动态过程时，它却无能为力。吸收一个光子会将系统激发到激发态，这需要一个更强大的理论透镜，能够追踪电子随时间的运动。

本文深入探讨了含时密度泛函理论 (TD-DFT)，它是为捕捉这些动态现象而对 DFT 进行的扩展。它解决了基态方法留下的关键空白，提供了一种计算上可行的方法来预测和理解电子激发。在接下来的章节中，您将学习 TD-DFT 的核心原则，探索其优美的原理、意想不到的失败以及使其成为不可或缺的工具的巧妙解决方案。我们将首先揭示驱动该理论的基本原理和机制。随后，我们将探索其变革性的应用和跨学科联系，揭示 TD-DFT 如何被用于设计新材料、解开生物之谜和解读复杂的光谱数据。

原理与机制

在上一章中，我们接触到了这样一个思想：密度泛函理论 (DFT) 不仅可以描述世界最安静的状态，还可以描述它与光剧烈相互作用时的情景。这种扩展就是含时密度泛函理论 (TD-DFT)，它是我们理解颜色、光化学和一系列动态过程的关键。但它是如何工作的呢？幕后转动的齿轮和杠杆是什么？让我们拉开帷幕一探究竟。就像任何精良的机器一样，它的设计遵循一些优美的原则，而它偶尔的故障与其平稳运行时一样具有启发性。

最初的猜测及其不足之处

当一个分子吸收一个光子时，一个电子从一个已占据的轨道被提升到一个先前未占据的轨道。如果我们已经进行过基态 DFT 计算，我们就会有一张包含所有分子轨道及其能量的清晰图表。要激发一个分子，其最低能量最显而易见的猜测就是最高占据分子轨道 (HOMO) 和最低未占分子轨道 (LUMO) 之间的能量差。这看起来如此简单，如此直观。你有一个能量阶梯，HOMO 及以下的所有梯级都被电子占据，而从 LUMO 开始向上的梯级都是空的。电子最小的可能“跳跃”就是从 HOMO 到 LUMO，对吗？

然而，物理学很少如此简单。这个 HOMO-LUMO 能量差是一个不错的初步猜测，一个“零阶”近似，但它通常错得离谱。为什么？因为它忘记了谜题中一个关键的部分：当电子跃迁时，它在 HOMO 中留下了一个带正电的“空穴”。处于 LUMO 的受激电子和它留下的空穴现在彼此相互作用，通常是吸引作用。真正的激发能必须考虑到这种电子-空穴相互作用，而简单的轨道能量差完全忽略了这一点。

想象一下我们正在研究一种新的有机染料分子。一次基态 DFT 计算可能会告诉我们 HOMO-LUMO 能隙是 $3.542$ eV。但是一次考虑了完整电子响应的、恰当的 TD-DFT 计算，可能会揭示出第一个真实的激发能是 $3.851$ eV——相差超过 $0.3$ eV！。这不是一个小的舍入误差；这可能是一个分子被预测为黄色还是绿色的区别。为了做得更好，我们需要一个理论，它不仅仅是观察静态的能量阶梯，而是描述整个系统在响应光子“踢”一下时的行为。

让电子起舞

TD-DFT 的基本见解，由 Runge-Gross 定理正式确立，与其基态理论的对应部分一样深刻。它指出，对于一个给定的初始量子态，随时间演化的电子密度 $n(\mathbf{r}, t)$ 和导致其演化的含时外势 $v_{ext}(\mathbf{r}, t)$ 是唯一相关的。换句话说，电子密度的舞蹈包含了关于系统的所有信息。

这使我们能够再次使用基态 DFT 中的绝妙技巧：Kohn-Sham 系统。我们虚构一个无相互作用的电子系统，通过设计，它能重现与我们真实的、有相互作用的系统完全相同的含时密度 $n(\mathbf{r}, t)$ 。这些虚构的电子在一个精心构造的含时有效势 $v_{KS}[n](\mathbf{r}, t)$ 中运动，该势包括外势、经典的 Hartree 排斥势，以及我们那位至关重要又神秘的朋友——含时交换相关 (XC) 势。通过求解这些无相互作用电子的含时薛定谔方程，我们可以观察真实密度是如何演化的。

这就是 TD-DFT 的核心机器。现在的问题是：我们如何使用这个让电子跳舞的机器来找到分子可以演奏的特定能量“音符”？

两种方式聆听音乐

事实证明，有两种优美且互补的方法可以从我们的含时 Kohn-Sham 系统中提取电子光谱。想象一口钟。你如何找出它的固有振铃频率？

一种方法是用锤子敲击它然后聆听。这就是实时 (RT) TD-DFT的精神。我们从分子的基态开始，然后施加一个非常短而尖锐的电场脉冲——这相当于计算中的锤击。这个“踢”包含很宽的频率范围，原则上可以同时激发所有可能的电子跃迁。然后我们只需让 Kohn-Sham 系统随时间演化，并追踪分子偶极矩 $\boldsymbol{\mu}(t)$ 来回摆动的情况。这个时间信号 $\boldsymbol{\mu}(t)$ 的傅里叶变换揭示了一个谱图，其峰值恰好出现在分子的共振频率上——也就是它的激发能！。

这种实时方法非常强大。一次模拟原则上就可以给出宽能量范围内的整个吸收光谱。其分辨率仅受我们愿意“聆听”摆动偶极矩的时间长短限制。此外，因为它求解了完整的含时方程，RT-TD-DFT 不局限于弱微扰；它是模拟强激光场和其他高度非线性现象中电子行为的首选方法。它甚至可以通过让模拟的电子密度飞走来描述电离过程，即电子被完全从分子中射出的过程。

找到钟的频率的第二种方法更为巧妙。不是敲击它，而是对着它哼唱，并改变你的音高。当你的哼唱与钟的一个固有频率匹配时，它会突然开始强烈振动以响应。这就是共振。这也是线性响应 (LR) TD-DFT的精神。该方法提出了一个数学问题：“对于一个频率为 $\omega$ 的微小振荡微扰，电子密度的响应在哪些频率下会激增至无穷大？”这些频率就是系统的激发能。

这个问题被优雅地重新表述为一个矩阵本征值问题，在量子化学中被称为Casida 方程。该矩阵由基态 Kohn-Sham 轨道能和从 Hartree 势及 XC 势导出的耦合项构成。这个矩阵的本征值是激发能 $\omega$ ，而相应的本征矢告诉我们每次激发的特性（例如，“95% 是一个 HOMO 到 LUMO 的跃迁”）。如果你只需要前几个最低能量的激发，LR-TD-DFT 通常更高效，并且它为每个激发态提供了清晰的定量描述，这对于化学解释来说是无价的。

机器中的幽灵：一个不完美理论的美丽缺陷

没有任何物理理论是完美的，除非是完整的的多体薛定谔方程。我们在 TD-DFT 中所做的近似使其在计算上变得可行，但也留下了系统性的“幽灵”——这些失败不是随机的，而是具有深刻的启发性。事实上，理解 TD-DFT 为何在某些情况下会失败，比它总能成功更能教会我们关于量子力学的知识。最常见的近似是使用半局域 XC 泛函（如 GGA）并假设 XC 势是绝热的——也就是说，它对密度的变化做出瞬时响应，没有对过去的记忆。让我们看看这会带来什么麻烦。

长程关系问题

考虑一个由两部分组成的分子：一个富电子的给体 (D) 和一个缺电子的受体 (A)，它们相距很远，距离为 $R$ 。现在想象一个激发过程，一个电子从给体转移到受体，形成一个 $D^{+} \dots A^{-}$ 态。这个电荷转移 (CT) 激发的真实能量必须考虑三件事：从 D 移走电子所需的能量（其电离势 $I_D$ ），A 捕获电子时释放的能量（其电子亲和能 $A_A$ ），以及至关重要的，D 上新形成的正电荷与 A 上的负电荷之间的库仑吸引力。因此，精确的能量对分离距离有一个非常特定的依赖关系：

E_{\mathrm{CT}}^{\mathrm{exact}}(R) \approx I_{\mathrm{D}} - A_{\mathrm{A}} - \frac{e^2}{4\pi\varepsilon_0 R}

随着 $R$ 减小，能量降低（吸引力增强）。那么，标准的 TD-DFT 会预测什么呢？“绝热”XC 核是“短视的”。它是一个局域或半局域函数，意味着它只关心密度在当前位置的情况，而不关心远处。当电子和它留下的空穴相距很远时，这个短视的核无法看到它们的长程库仑相互作用。它完全忽略了 $-1/R$ 项！。结果是灾难性的失败：TD-DFT 预测的 CT 能量随距离几乎不变，并且通常严重偏低。

这不仅仅是一个学术上的好奇。从摄影到生物学，花青素染料无处不在，它们是这个问题的完美例子。它们的最长波长吸收是一个具有显著 CT 特性的跃迁。当你增加染料中共轭链的长度时，你增加了电子和空穴的平均分离距离。正如所预测的那样，标准的 TD-DFT 计算系统性地低估了激发能，并且随着染料变长，误差越来越严重——这是这一根本缺陷的直接体现。其根本原因在于大多数 XC 泛函固有的自相互作用误差，这使得 XC 势在长程范围内衰减得太快，从而搞乱了轨道能量，并使响应核对长程物理现象“视而不见”。

“一次两个”的难题

标准形式的 LR-TD-DFT 是基于单电子跃迁——即单粒子-单空穴 ( $1p1h$ ) 组态——来构建其激发态的。如果分子的一个真实激发态涉及两个电子同时被激发（一个双重激发，或 $2p2h$ 态）会发生什么？标准的绝热 TD-DFT 对它们完全“视而不见”。罪魁祸首是“绝热”近似，它假设 XC 核是频率无关的。一个频率相关的核会拥有对其他激发的“记忆”，使其能够从单个激发的组合中构建出像双重激发这样的状态。而没有记忆的绝热核却做不到。在实时模拟中，你根本不会在光谱中应该出现双重激发的地方找到任何峰值。

多米诺效应：断裂的路径

这些失败对光化学有可怕的后果。当分子吸收光时，它们通常通过超快的非辐射路径耗散能量。这种快速传输的枢纽是锥形交叉——几何空间中两个电子势能面接触的点，为分子从较高能态转换到较低能态提供了一个通道。这些通道的存在和位置决定了光激发分子的命运。

但是，如果锥形交叉涉及一个具有强电荷转移或双重激发特性的状态呢？TD-DFT 对这些状态已知存在缺陷，因此会得到完全错误的势能面。它可能会错误地定位交叉点，或者在应该有真实交叉的地方预测出一个“避免交叉”，或者完全搞错通道周围的拓扑结构。这可能导致对分子光化学行为的预测完全错误。

物理学家如技师：调校泛函

故事并没有以失败告终。通过理解机器为何会坏，我们可以设计出更好的零件。特别是电荷转移问题，催生了一项绝妙的创新：范围分离杂化 (RSH) 泛函。

这个想法体现了优美的物理直觉。我们知道，半局域 XC 泛函（如 GGAs）在描述电子靠得很近的短程电子关联方面相当不错。我们也知道它们在长程上表现糟糕。相反，来自 Hartree-Fock 理论的精确交换能是非局域的，能正确描述长程相互作用（比如我们的 $-1/R$ 问题），但常常忽略重要的短程关联。那么，为什么不取二者之长呢？

RSH 泛函正是这样做的。它们使用一个平滑的数学函数将电子-电子相互作用分为短程和长程部分。然后对这两部分采用不同的处理方式：

在短程：它们使用标准的半局域泛函，有时会混合少量精确交换。
在长程：它们切换到使用 100% 的非局域、精确的 Hartree-Fock 交换。

这种双管齐下的方法效果奇佳。长程精确交换解决了电荷转移中的两个关键问题。首先，它修正了基态 XC 势，迫使其具有正确的 $-1/r$ 渐近衰减行为，这反过来又修正了由自相互作用误差导致的错误轨道能。其次，它为 TD-DFT 响应核提供了必要的非局域分量，使其能够正确地捕捉到远处电子和空穴之间的 $-1/R$ 吸引相互作用。这证明了理论物理学的力量：对一个根本性缺陷的深刻理解，导出了一个优美的、基于物理动机的解决方案，极大地扩展了我们模型的预测能力。它向我们展示，即使在电子复杂的舞蹈中，简单而优美的原则依然占据主导。

应用与跨学科联系

在上一章中，我们深入探讨了含时密度泛函理论的核心，探索了使我们能够模拟电子对光线闪烁响应的原理和机制。现在我们拥有了一个强大的理论透镜。然而，真正的冒险始于我们将这个透镜转向我们周围的世界。我们能看到什么？我们能学到什么？我们能构建怎样的新世界？这不仅仅是一次计算练习；这是一段探索颜色量子起源、设计未来技术、生命复杂机制以及物质本质的旅程。

用电子作画：世界的色彩

你是否曾想过是什么让胡萝卜呈橙色，或让夏日天空呈蓝色？在最根本的层面上，答案在于电子对光的选择性吸收。分子，就像微小的乐器，只能在特定频率的光下共振。我们的眼睛感知到的是那些未被吸收的频率。几个世纪以来，我们凭经验理解这一点，但有了 TD-DFT，我们可以从第一性原理出发进行预测。

想象一个简单的分子家族：多烯。它们只是一串交替出现单键和双键的碳原子链。链越长， $\pi$ -电子可以漫游的“空间”就越大。这个延展的“跑道”产生了深远的影响。随着链的增长，最高占据分子轨道 (HOMO) 和最低未占分子轨道 (LUMO) 之间的能隙会缩小。由于光的能量与其波长成反比 ( $E=hc/\lambda$ )，较小的能隙意味着分子吸收更长波长的光。TD-DFT 使我们能够以非凡的精度计算这个吸收波长 $\lambda_{\max}$ 。通过对长度递增的多烯进行一系列计算，我们可以追踪 $\lambda_{\max}$ 从不可见的紫外区移动到可见光谱区。当我们延长分子的共轭链时，我们简直可以在电脑屏幕上看到一个分子获得颜色。这个过程不仅仅是一个思想实验；它是材料科学中用于设计新染料和颜料的标准且可靠的计算方案。

理论家的工具箱：一个关于成功、失败与更深层次真理的故事

没有一种科学工具是完美的，而它的局限性往往比它的成功更能给人以启发。TD-DFT 的故事就是对此的一个绝佳例证。它迫使我们更深入地思考电子激发的本质。

考虑电子被激发的两种不同方式。在像甲醛这样的分子中，一个电子可以从氧原子上的一个非键轨道被轻轻推到附近的某个反键 $\pi^*$ 轨道。这是一个局域事件，一次小小的跳跃。标准的 TD-DFT，甚至像 $\Delta$ SCF 方法这样更简单的方法，都能很好地描述这种 $n \to \pi^*$ 跃迁。

但是，当一个电子被转移到更远的距离时会发生什么？想象一个“给体”分子与一个“受体”分子连接在一起。一次激发可以是一个电子离开给体并到达受体，形成一个“电荷转移”(CT) 态。在这里，标准的 TD-DFT，特别是使用我们最初遇到的那些更简单的局域或半局域泛函时，会遭遇惨败。为什么？原因很深刻。这些泛函有一种“近视”。它们产生的交换相关势（包含了所有棘手的多体物理）随距离衰减得太快。就好像电子被拴在一条短绳上；该理论根本无法“看到”远处的电子和它留下的空穴之间应该感受到的强烈的长程库仑吸引力 ( $V(R) \propto 1/R$ )。结果，它极大地低估了这种电荷分离所需的能量。

这一失败的发现并非一次挫败，而是一次突破。它促成了“范围分离”泛函的开发，这些泛函被巧妙地设计成在长程时切换到正确的行为。这些不仅仅是随意的数学修正；它们代表了一种已被明确编码进理论的更深层次的物理理解。

还有另一种电子运动让标准的 TD-DFT 感到棘手：两个电子同步起舞。大多数激发涉及单个电子从占据轨道跳到未占据轨道。但有些态，比如长链多烯中著名的 $2\,{}^{1}A_g$ 态，具有主要的“双重激发”特性。标准的绝热 TD-DFT 建立在单激发框架之上，难以描述这些难以捉摸的状态。奇怪的是，在这个问题上，一些更古老、更简单的模型，如 Pariser–Parr–Pople (PPP) 方法，虽然只关注 $\pi$ -电子，但通过组态相互作用更明确地处理了它们的相互作用，反而能提供更清晰的见解。这教给我们科学中一个至关重要的教训：一个使用不当的复杂工具，其价值不如一个能抓住核心物理的简单工具。

前沿与跨学科联系

手握一个成熟且被充分理解的工具，我们现在可以涉足迷人而复杂的领域，从工程新技术到破译生命和物质的秘密。

光之工程：OLED中的单重态-三重态探戈

你手机或电视中的现代显示器可能使用了有机发光二极管 (OLED)。OLED 设计中的一个关键挑战是，电激发会同时产生单重态（自旋配对， $\uparrow\downarrow$ ）和三重态（自旋平行， $\uparrow\uparrow$ ）激发态，比例通常为 1:3。三重态是“暗态”，通常在不发光的情况下衰变，浪费了 75% 的能量。一个革命性的解决方案在于那些表现出热活化延迟荧光 (TADF) 的分子。这些分子被设计成在其最低单重态 ( $S_1$ ) 和三重态 ( $T_1$ ) 之间具有非常小的能隙 $\Delta E_{S-T}$ 。这个微小的能隙使得暗态的三重态可以利用环境的热能转换回明亮的单重态，从而极大地提高了器件的效率。

TD-DFT 已成为理性设计 TADF 材料不可或缺的工具。化学家现在可以通过计算 $\Delta E_{S-T}$ 来对数百种候选分子进行计算筛选。这是一项高风险的计算，通常涉及电荷转移态，并且需要一个复杂的工作流程：使用范围分离泛函来获得正确的能量，仔细检查三重态的自旋态，有时还将 TD-DFT 与其他方法结合以确保稳健性。这是基础理论直接指导前沿技术的一个绝佳范例。这些计算的准确性对于预测动力学速率（例如对能隙呈指数敏感的系间窜越速率）也至关重要。

解码生物学：视觉的量子力学

我们如何将我们的量子透镜应用于广阔而杂乱的生物世界？一个蛋白质由数千个原子组成。一次完整的量子计算是不可能的。解决方案既优雅又强大：多尺度模拟。量子力学/分子力学 (QM/MM) 方法用高水平的量子方法处理系统中最重要的部分——发生化学反应的“活性位点”，而广阔的蛋白质其余部分及其环境则用更简单的经典物理来处理。

考虑视紫红质，我们眼中负责视觉的蛋白质。当一个光子撞击到藏在蛋白质内部的一个小小的视黄醛分子时，一切就开始了。这个单一的量子事件引发了视黄醛的形状改变（异构化），从而启动了成为视觉的信号级联。为了模拟这个过程，我们可以在视黄醛及其直接的关键邻居（如“调节”其颜色的带电谷氨酸残基）周围画一个小的气泡。这个区域用量子方法处理。对于像异构化这样涉及断键的过程，通常需要更高级的多参考方法，但 TD-DFT 可以在表征态或模拟 QM 区域中不太关键的部分中发挥关键作用。蛋白质的其余部分被建模为一个经典支架，提供正确的结构和静电环境。由 TD-DFT 等方法驱动的 QM/MM，让我们能够逐个原子、逐个电子地观察一个基本生物过程的展开。

探测内层：来自X射线的瞥见

到目前为止我们讨论的光——可见光和紫外光——能量足以激发最外层的价电子。如果我们使用能量高得多的光，比如X射线，会怎样呢？X射线可以从原子核附近的深层内层轨道中敲出电子。这就是X射线吸收谱的世界，一种探测元素组成和化学环境的强大工具。

用 TD-DFT 模拟这些内层激发带来了一个新的挑战：轨道弛豫。当一个内层电子被剥离时，就像从一栋建筑中移走一根柱子。整个电子结构会“颤抖”并收缩，以响应那个高度集中的正电荷空穴。建立在基态轨道之上的标准 TD-DFT 难以捕捉这种剧烈的重排。一个巧妙的解决方案见于像跃迁势DFT (TP-DFT) 这样的方法中，它在一个内层轨道仅被半占据的“跃迁态”中计算激发。这个优雅的技巧将初始态和最终态的弛豫效应的平均值直接构建到计算中，从而得出惊人准确的X射线吸收能。

对于更重的原子，物理的另一层规律也开始显现：爱因斯坦的相对论。对于在重原子核附近高速运动的电子，相对论效应变得显著。其中最重要的是自旋-轨道耦合，它将电子的自旋与其轨道运动联系起来。这种耦合足够强，可以将 $2p$ 内层轨道分裂成两个不同的能级， $2p_{3/2}$ 和 $2p_{1/2}$ 。这导致X射线吸收信号分裂成两个峰，即 $L_3$ 和 $L_2$ 边。为了模拟这一点，TD-DFT 必须与相对论量子力学相结合，从而产生了明确包含自旋-轨道耦合的双分量和四分量形式。这是物理学统一性的惊人展示，一个重金属配合物中电子的舞蹈由量子力学和狭义相对论的联合定律支配。

平坦（且充满激子）的世界：二维材料

我们的旅程在材料科学的前沿结束：二维材料，如石墨烯和过渡金属二硫族化物。在这些原子级薄的半导体中，一个吸收的光子产生一个电子和一个空穴，它们之间相互强烈吸引，形成一个被称为激子的紧密束缚的准粒子。

在这里，我们再次面临 TD-DFT 的长程问题。这些激子中的电子-空穴吸引力是一种长程力，标准 TD-DFT 的简单、短程核再次显得力不从心。虽然可以开发更高级的核，但这正是多体微扰理论以 GW-Bethe-Salpeter 方程 (GW-BSE) 方法的形式大放异彩的领域。GW-BSE 方法明确包含了一个被屏蔽的长程库仑相互作用，为这些强束缚激子提供了更严谨、更可靠的描述。TD-DFT 与 GW-BSE 在这些系统上的比较是一个活跃的研究领域，推动着理论家们开发出更好的泛函，以弥合效率和准确性之间的差距。

从一朵花的颜色到太阳能电池的效率，从视觉行为到能够改变我们世界的新材料的创造，电子响应光的故事就是我们世界的故事。TD-DFT，尽管有其种种近似和局限性，已被证明是讲述这个故事的一个无比宝贵的工具。它证明了一个好的物理思想的力量，一个继续将充满活力、动态的量子世界带入更清晰焦点的透镜。