Generalized length-gauge time-dependent Liouville equation

技术摘要：广义长度规范时间相关刘维尔方程

问题陈述
本文针对长度规范（length-gauge, LG）时间相关刘维尔方程在理论表述中长期存在且被忽视的差异进行了研究。该方程是模拟固体超快动力学的核心工具。该方程通过对易子 $[\hat{H}, \hat{\rho}]$ 来控制密度算符 $\hat{\rho}$ 的时间演化，其中长度规范下的相互作用哈密顿量为 $\hat{H}_I = e\mathbf{r} \cdot \mathbf{E}(t)$。

核心冲突在于如何计算位置算符 $\mathbf{r}$ 与密度算符之间的对易子矩阵元，即 $\langle nk | [\mathbf{r}, \hat{\rho}] | mk' \rangle$。文献中存在两种无法调和的版本：

1. Kohn 和 Luttinger (KL) (1954)： 他们的表述产生了两个关于晶格动量 $k$ 和 $k'$ 的偏导数，且没有强制要求 $k=k'$ 的狄拉克 $\delta$ 函数。
2. Aversa 和 Sipe (AS) (1995)： 他们被广泛采用的表述仅产生一个偏导数，并包含一个狄拉克 $\delta$ 函数 $\delta(k-k')$，从而有效地将密度矩阵限制为在晶格动量上是对角化的（$k$-diagonal）。

尽管 KL 的表述在历史上更早，且在物理上与位置算符缺乏平移不变性的特性相一致，但 AS 的表述数十年来一直主导着文献，且未受到审视。作者认为，AS 方法错误地假设了位置算符在 $k$ 空间是对角的，忽略了由于 $\mathbf{r}$ 打破平移对称性而产生的本质的 $k$-非对角元素。

方法论
作者开发了一种广义理论形式，在不假设位置算符在 $k$ 空间是对角的情况下，推导出正确的长度规范刘维尔方程。

• 推导： 从海森堡运动方程和晶格动量表征（CMR）下的密度算符定义出发，作者利用布洛赫波函数的正交归一化条件。通过对正交条件 $\langle \psi_n(k) | \psi_m(q) \rangle = \delta(q-k)\delta_{nm}$ 关于 $q$ 和 $k$ 求偏导，作者推导出了位置算符的矩阵元 $\langle nk | \mathbf{r} | mq \rangle$。
• 核心见解： 推导表明，位置算符矩阵元包含一个与 $\nabla_q \delta(q-k)$ 成正比的项（导致导数）以及一个与 $\delta(q-k)$ 成正比的项（贝里联络/Berry connection）。至关重要的是，在评估对易子 $[\mathbf{r}, \hat{\rho}]$ 时，位置算符分别出现在 $\hat{\rho}$ 的左侧和右侧。这要求必须对初始动量 $k_1$ 和最终动量 $k_2$ 同时进行求导，从而产生一个项 $(\nabla_{k_1} + \nabla_{k_2})\rho$。
• 极限情况分析： 作者证明，只有在相对动量 $Q = (k_1 - k_2)/2$ 被设为零（$Q=0$）时，才能恢复 AS 方程。
• 数值验证： 为了测试 $k$-对角（AS）与全-$k$（广义）形式在物理后果上的差异，作者在三个模型系统中进行了数值模拟：
  1. 二位点模型，用于演示电荷密度的非周期性。
  2. 具有周期性边界条件的一维金属模型（三能带）。
  3. 具有周期性边界条件的半导体模型（三能带）。
  4. Su-Schrieffer-Heeger (SSH) 链，使用开边界条件以测试在量子受限系统中的适用性。

主要贡献

1. 解决分歧： 本文确立了广义长度规范刘维尔方程必须包含两个偏导数 ($\nabla_{k_1} + \nabla_{k_2}$) 且不含狄拉克 $\delta$ 函数，这与 Kohn 和 Luttinger 的原始理论保持一致。
2. 识别误差来源： 作者指出，AS 形式的错误源于假设位置算符在 $k$ 空间是对角的。这一假设迫使密度矩阵保持在 $k$ 上是对角的，忽略了由电场引起的不同晶格动量之间的耦合。
3. 广义形式： 作者提供了一个完整的推导过程，其中密度矩阵 $\rho(nk_1; mk_2)$ 被视为一个全 $k$ 矩阵，允许即使在均匀电场下也能动态生成 $k$-非对角元素。

结果

• 定性差异： 在金属和半导体模型上的数值测试表明，$k$-对角形式无法捕捉到不同 $k$ 点之间的耦合。在广义形式中，电场会产生非零的密度矩阵非对角元素（即 $\rho(k_i, k_j)$，其中 $i \neq j$），即使它们初始为零。
• 定量差异： 最显著的发现是物理观测量的误差量级。
  • 能量吸收： 与全-$k$ 形式相比，$k$-对角形式低估了系统吸收的总能量，其差距约为 7 倍。
  • 极化： 同样地，使用 $k$-对角方法计算出的诱导极化在定性上也有所不同（显示为双峰而非单峰），并且在定量上也小了约 7 倍。
• 鲁棒性： 这些结果在金属和半导体中均成立，且该形式被证明适用于 $k$ 不是良好量子数的系统（如 SSH 模型，具有开边界条件），进一步验证了该方法的普适性。

意义
本文声称，基于 Aversa 和 Sipe 形式的标准 $k$-对角刘维尔方程，在超快模拟中会系统性地低估能量吸收和极化。作者断言，为了准确模拟金属、半导体和绝缘体的超快动力学，必须使用他们提供的这种恢复了全 $k$-非对角结构的广义形式。因此，本文暗示，以往仅基于 $k$-对角方程的研究结果可能需要重新审视。这项工作为长度规范模拟建立了一个修正后的范式，为更准确地研究实际材料中的超快现象开启了大门。