来源论文: https://arxiv.org/abs/2602.17959v1 生成时间: Feb 23, 2026 05:16

0. 执行摘要

本研究论文提出了一种创新性的方法来计算二维哈伯德模型的光学和霍尔电导率，特别关注于铜氧化物超导体中观察到的非费米液体行为。传统上，在强关联电子系统中研究电磁响应面临两大挑战：一是理论模型中常见的非物理性发散，二是在量子蒙特卡罗模拟中臭名昭著的“符号问题”。作者通过引入一种基于局部矩涨落的有效理论描述，并巧妙地利用哈伯德-斯特拉托诺维奇变换，成功地将系统的电磁响应表达为非相互作用系统在涨落局部矩背景下的系综平均。更值得注意的是，他们证明了在许多物理相关情况下，特别是当局部矩涨落的有效作用采用广泛使用的Millis-Monien-Pines（MMP）唯象磁化率形式时，蒙特卡罗模拟可以避免符号问题。这为在传统方法难以处理的强关联系统中进行可靠的非扰动计算打开了大门。研究结果揭示了铜氧化物超导体的光学和霍尔电导率具有显著的双组分结构，包括低能区的Drude分量和高能区的近红外分量。这些分量被解释为反铁磁有序态中谱特征的残余，并且霍尔响应中的符号变化被认为是费米面重构（特别是电子口袋出现）的前兆。尽管本研究集中于热涨落主导的极限，并为定性理解提供了重要见解，但该框架的通用性预示着未来在考虑量子涨落和直流输运性质方面的进一步扩展，为解决铜氧化物超导体这一复杂难题奠定了坚实基础。

1. 核心科学问题、理论基础、技术难点、方法细节

1.1 核心科学问题

铜氧化物超导体（cuprates）的正常态表现出显著的非费米液体（non-Fermi liquid）行为，这被认为是理解高温超导机制的关键。本文的核心科学问题在于，如何从理论上准确描述和计算这些强关联电子系统的电磁响应，特别是光学电导率 ($\sigma_{xx}(\omega)$) 和霍尔电导率 ($\sigma_{xy}(\omega)$)，并解释其独特的实验现象。具体挑战包括：

  1. 非Drude光学吸收谱： 铜氧化物超导体的光学电导率随频率衰减远慢于传统费米液体金属，表现出低频Drude峰和宽广的近红外峰的双组分结构。这种谱重分布与温度和掺杂水平的变化广泛相关，与传统费米液体金属截然不同。
  2. 奇异的霍尔响应： 正常态下的霍尔数 ($n_H = 1/R_H$) 随温度几乎线性增加，并在低温下（强磁场抑制超导）随掺杂水平在$x^*$处发生戏剧性交叉。这暗示着费米面重构，但实验上缺乏平移对称性破缺的证据。这些现象与传统费米液体的行为形成鲜明对比，引发了对有效载流子密度和超导凝聚体起源的疑问。
  3. 局部矩的涨落作用： 作为掺杂Mott绝缘体的标志，电子谱重分布伴随着涨落局部矩的出现。这些局部矩在RIXS实验中被证实普遍存在，即使系统远离磁有序态。电子的双重性质（巡游准粒子和局部矩）构成了理论描述的严峻挑战。

1.2 理论基础

本研究采用了基于自旋-费米子模型的有效理论描述框架，其核心思想是将强关联哈伯德模型中的电子相互作用，通过引入涨落的局部矩场进行解耦。具体步骤如下：

  1. 哈伯德模型的重新表述： 传统的二维哈伯德模型 Hamiltonian 为：

    $$H = - \sum_{i,j,\sigma} t_{i,j} c_{i,\sigma}^{\dagger} c_{j,\sigma} + U \sum_i n_{i,\uparrow} n_{i,\downarrow} - \mu \sum_{i,\sigma} n_{i,\sigma}$$

    为了引入局部矩场，作者将相互作用项 $U n_{i,\uparrow} n_{i,\downarrow}$ 重新表述为：

    $$U n_{i,\uparrow} n_{i,\downarrow} = \frac{2U}{3} S_i \cdot S_i + \frac{U}{2} (n_{i,\uparrow} + n_{i,\downarrow})$$

    其中 $S_i = \frac{1}{2} c_{i,\sigma}^{\dagger} \boldsymbol{\sigma}_{\sigma,\sigma'} c_{i,\sigma'}$ 是在位点 $i$ 的自旋密度算符。通过将最后一项吸收到化学势的重新定义中，Hamiltonian变为：

    $$H = - \sum_{i,j,\sigma} t_{i,j} c_{i,\sigma}^{\dagger} c_{j,\sigma} - \frac{2U}{3} \sum_i S_i \cdot S_i - \mu \sum_{i,\sigma} n_{i,\sigma}$$

  2. 电磁响应的引入（Peierls 替代）： 为了计算电磁响应，电子与电磁场通过Peierls替代进行耦合。即，跳跃积分 $t_{i,j}$ 被修改为 $t_{i,j} e^{i A_{i,j} \cdot (r_i - r_j)}$，其中 $A_{i,j}$ 是连接位点 $i$ 和 $j$ 的矢量势。Hamiltonian 变为 $H[A]$：

    $$H[A] = - \sum_{i,j,\sigma} t_{i,j}^{A} c_{i,\sigma}^{\dagger} c_{j,\sigma} - \frac{2U}{3} \sum_i S_i \cdot S_i - \mu \sum_{i,\sigma} n_{i,\sigma}$$

    其中 $t_{i,j}^A = t_{i,j} e^{i A_{i,j} \cdot (r_i - r_j)}$。

  3. 哈伯德-斯特拉托诺维奇 (HS) 变换和有效作用量： 为了处理 $S_i \cdot S_i$ 项，引入矢量HS场 $\phi_i$ 来描述局部矩的涨落。通过HS变换，相互作用项被解耦，得到一个含有费米子场和HS场的新Hamiltonian $H[\psi, \psi^{\dagger}, \phi, A]$。然后，通过对费米子场进行高斯积分，得到仅依赖于HS场 $\phi$ 和矢量势 $A$ 的有效作用量 $S_{\text{eff}}[\phi, A]$：

    $$Z[A] = \int D[\phi] e^{-S_{\text{eff}}[\phi,A]}$$

    其中 $S_{\text{eff}}[\phi, A] = \int_0^{\beta} d\tau \left( \frac{2U}{3} \sum_i \phi_i \cdot \phi_i - \text{Tr} \ln[-G^{-1}[\phi]] \right)$。 $G^{-1}[\phi]$ 是在HS场存在且没有EM场的情况下，电子的逆格林函数。

  4. 电导率的计算： 电磁响应内核 ($\kappa^{xx}(\omega)$ 和 $\kappa^{xyy}(\omega)$) 是通过对自由能 $F[A] = -T \ln Z[A]$ 进行关于EM势 $A$ 的泛函导数展开得到的，精确到三阶。具体地，电流 $j^{\alpha}(q)$ 与势 $A^{\alpha}(q)$ 的关系为 $j^{\alpha}(q) = -\frac{\delta F[A]}{\delta A^{\alpha}(-q)}$。光学电导率 $\sigma_{xx}(\omega)$ 和霍尔电导率 $\sigma_{xy}(\omega)$ 则由这些内核通过Wick旋转 ($i\omega_n \to \omega + i0^+$) 得到：

    $$\sigma_{xx}(\omega + i0^+) = \frac{i}{\omega + i0^+} \kappa^{xx}(\omega + i0^+)$$

    $$\sigma_{xy}(\omega + i0^+) = -\frac{1}{i\omega_n q_x} \frac{\partial}{\partial q_x} \kappa^{xyy}(q_x, \omega + i0^+) \big|_{q_x\to 0}$$

    其中内核 $\kappa$ 的表达式涉及格林函数 $G$ 和电磁顶点函数 $v^{(n)}$ 的迹运算，这些迹运算通过求和电子在给定 $\phi$ 场配置下的本征态和本征能量来实现。这是核心计算步骤，通过对HS场进行蒙特卡罗采样，并对每个采样配置下的费米子迹进行平均来完成。

1.3 技术难点与解决方案

  1. 符号问题（Sign Problem）：

    • 难点： 在量子蒙特卡罗模拟中，当 $S_{\text{eff}}[\phi]$ 不是实数时，$e^{-S_{\text{eff}}[\phi]}$ 将不是正定的，导致蒙特卡罗权重出现负值，使得统计采样无法进行。这在强关联系统中是一个普遍且臭名昭著的难题。
    • 解决方案： 本文的重大突破在于识别出几种可以避免符号问题的情况：
      • 高温极限 ($\beta \to 0$)： 在此极限下，HS场 $\phi$ 的时间依赖性可以忽略，从而 $S_{\text{eff}}[\phi]$ 变为实数，模拟无符号问题。论文的数值计算主要关注此极限。
      • 高斯极限（扰动展开）： 在对相互作用强度 $U$ 的扰动展开中，系统自由能的所有项都是实数，因此 $e^{-S_{\text{eff}}[\phi]}$ 在高斯近似下是正定的。
      • Millis-Monien-Pines (MMP) 唯象磁化率： 最关键的是，作者指出，当涨落局部矩的有效作用量采用广泛使用的MMP形式时，即 $S_{\text{eff}}[\phi] = \frac{1}{2} \sum_q \chi^{-1}(q, i\omega_n) \phi_q \phi_{-q}$，其中 $\chi(q, i\omega_n)$ 是唯象磁化率，且在虚频率域是实数，则 $e^{-S_{\text{eff}}[\phi]}$ 始终是正定的。这使得在非扰动水平上对局部矩涨落进行蒙特卡罗模拟成为可能，从而避免了符号问题。即使在有效作用量中包含四次项，此结论也成立。

  2. 非物理性发散和顶点修正：

    • 难点： 之前基于自旋-费米子模型的扰动理论在计算低频电导率时，由于格林函数极点在顶点函数中的合并，常遇到非物理性发散。
    • 解决方案： 本文的框架通过在给定 $\phi$ 场配置下精确积分出费米子，并随后对 $\phi$ 场进行蒙特卡罗平均，自然地以非扰动方式包含了这些效应。因此，其导出的电导率公式在热涨落主导的情况下不会出现这些发散，能够产生任意频率下的数值准确预测。

  3. 数值解析延拓：

    • 难点： 从虚频率数据通过Wick旋转到实频率，通常需要进行数值解析延拓，这是一个病态问题，且易受玻尔兹曼权重在虚时间域内指数抑制的影响。
    • 解决方案： 为了避免解析延拓的复杂性和不确定性，本研究主要聚焦于热涨落主导的极限，即忽略HS场 $\phi$ 的时间依赖性。在这种情况下，电子的格林函数在Matsubara频率上是对角的，使得实频率电导率的计算可以避免解析延拓，从而直接在实频率域进行。

  4. 有限尺寸效应：

    • 难点： 蒙特卡罗模拟通常在有限尺寸晶格上进行，可能引入有限尺寸效应。
    • 解决方案： 论文使用 $L \times L = 400$ 个格点的正方晶格（即 $20 \times 20$）并施加周期性边界条件。通过在计算霍尔电导率时使用小但非零的波矢 $q_x = 2\pi/L$ 来近似导数，以恢复正确的带内光学谱重并处理有限尺寸带来的对称性破缺。然而，这也可能导致一些微小的数值伪影，例如在 $\omega = 0$ 附近出现的负过冲。

1.4 方法细节总结

该方法的整体流程是：

  1. 从二维哈伯德模型出发，通过HS变换将其转化为费米子与经典（或半经典）局部矩场耦合的模型。
  2. 引入电磁场，通过Peierls替代将Hamiltonian扩展到 $H[A]$。
  3. 精确积分出费米子场，得到一个仅依赖于局部矩场 $\phi$ 和电磁势 $A$ 的有效作用量 $S_{\text{eff}}[\phi, A]$。
  4. 通过对 $S_{\text{eff}}[\phi, A]$ 关于 $A$ 的泛函导数展开，得到电磁响应内核 $\kappa^{xx}$ 和 $\kappa^{xyy}$ 的表达式，这些表达式包含了在 $\phi$ 场背景下的费米子迹。
  5. 利用蒙特卡罗方法对 $\phi$ 场配置进行采样，避免符号问题的关键在于选择合适的 $S_{\text{eff}}[\phi]$ 形式（例如MMP唯象磁化率或高温极限）。
  6. 对于每个采样配置，计算费米子格林函数及其迹，然后进行系综平均。
  7. 最后，通过Wick旋转或直接计算（在热涨落极限下）获得实频率下的光学和霍尔电导率。

这一框架提供了一个强大的工具，可以在非扰动水平上研究强关联电子系统的电磁响应，同时解决了长期困扰该领域的符号问题和非物理发散问题。

2. 关键benchmark体系、计算所得数据、性能数据

2.1 Benchmark体系与参数设定

本研究的数值计算是在一个有限尺寸的正方晶格上进行的，格点数为 $L \times L = 400$ (即 $20 \times 20$ 的二维晶格)。为了模拟无限大体系，在 $x$ 和 $y$ 方向均施加了周期性边界条件。

  • 模型： 二维哈伯德模型。
  • 跳跃积分：
    • 最近邻跳跃积分：$t$ (被用作能量单位)。
    • 次近邻跳跃积分：$t' = -0.3t$。
  • 化学势： $\mu = -t$。在非相互作用极限下，这对应于约15%的空穴掺杂水平。作者明确指出，为了在计算中固定掺杂水平，必须根据温度和耦合强度 $\bar{U}$ 调整化学势，但本文未进行这种细节处理，而是旨在进行定性讨论。
  • 温度： $k_B T = 0.1t$。
  • 局域矩有效作用量形式：
    1. 高斯作用量（MMP型磁化率）： $S_{\text{eff}}[\phi] = \frac{1}{2} \sum_q \chi^{-1}(q, i\omega_n) \phi_q \phi_{-q}$，其中 $\chi(q, i\omega_n) = \frac{\chi_0}{1 + (q-Q)^2\xi^2 + |\omega_n|/W_{sf}}$。
      • 在本文的热涨落主导极限下，忽略了时间依赖性，因此 $|\omega_n|/W_{sf}$ 项为零。因此 $\chi(q, i\omega_n) = \frac{\chi_0}{1 + (q-Q)^2\xi^2}$。在格点正则化后， $(q-Q)^2$ 替换为 $4 + 2[\cos(q_x) + \cos(q_y)]$。因此，有效作用量变为 $S_{\text{eff}}[\phi] = \frac{1}{2T\chi_0} \sum_i (4 + \xi^{-2}) \phi_i^2 - \frac{1}{2T\chi_0 \xi^2} \sum_{i,\delta} \phi_i \cdot \phi_{i+\delta}$。
      • 参数： $\bar{U}$（局部矩与电子的整体耦合强度）和 $\xi$（反铁磁关联长度）被视为独立参数。$\bar{U}$ 通过 $U = \bar{U}/\sqrt{T\chi_0}$ 间接影响耦合强度。
    2. 非高斯作用量（非线性 $\sigma$ 模型, NLSM）： $S_{\text{eff}}[\phi] = \eta \sum_{i,\delta} \phi_i \cdot \phi_{i+\delta}$，并伴随约束条件 $\phi_i \cdot \phi_i = 1$。
      • 参数： $\eta$ 是一个常数，用于调节关联长度。
  • 计算细节：
    • 蒙特卡罗采样：生成6400个统计独立的局部矩场 $\phi$ 样本。
    • $\delta$ 函数峰展宽：为了平滑谱线，所有的 $\delta$ 函数峰都被展宽为洛伦兹峰，其宽度为 $0.03t$。
    • 霍尔电导率导数近似：由于在 $q_x \to 0$ 极限下 $q_x$ 导数难以精确计算，文中采用了有限差分近似，即 $\frac{\partial}{\partial q_x} \kappa^{xyy}(q_x, \omega+i0^+) \big|_{q_x \to 0} \approx \frac{1}{q_m} \text{Im}[\kappa^{xyy}(q_x=q_m, \omega+i0^+) - \kappa^{xyy}(0, \omega+i0^+)]$，其中 $q_m = 2\pi/L$ 是布里渊区中最小的非零动量。

2.2 计算所得数据与结果

2.2.1 光学电导率 ($\sigma_{xx}(\omega)$)

  1. 高斯作用量结果 (图1)：

    • 小 $\bar{U}$： 光学电导率由低能Drude峰主导。
    • 大 $\bar{U}$： 光谱由宽广的近红外峰主导，延伸到能带宽度。
    • 中间区域： 出现双组分结构，包括低频Drude分量和高频近红外分量。随着局部矩关联长度 $\xi$ 和耦合强度 $\bar{U}$ 的增加，这种双组分特征变得越来越明显。
    • Drude分量的解释： 归因于热自旋涨落存在下费米能级附近的残余电子态密度。可以理解为在涨落局部矩凝聚成静态反铁磁有序态时，电子能带会分裂为上下SDW（自旋密度波）带，Drude分量对应于带内跃迁。
    • 近红外分量的解释： 对应于SDW分裂带之间的带间跃迁。
    • 与实验比较： 预测的Drude权重相对于强耦合区的实验观测值偏弱。这被归因于模型中忽略了局部矩的量子性质。量子效应预计会恢复费米能级附近的电子态密度，从而增强Drude峰的谱重。

  2. 非高斯作用量（NLSM）结果 (图2)：

    • 与高斯作用量相比，NLSM作用量下的双组分结构（Drude峰和近红外峰）更为明显。这与电子掺杂铜氧化物超导体的实验观测结果更吻合，NLSM描述被认为更适用于这些材料。
    • 无发散性： 本文的计算避免了先前扰动处理中遇到的低频非物理性发散，这对于揭示光学谱的双组分结构至关重要。

2.2.2 霍尔电导率 ($\text{Im}\sigma_{xy}(\omega)$)

  1. 高斯作用量结果 (图3)：

    • 与光学电导率类似，霍尔电导率的虚部 $\text{Im}\sigma_{xy}(\omega)$ 也呈现双组分结构。
    • 小 $\bar{U}$： 低能峰主导，且为正值。
    • 大 $\bar{U}$： 谱由从近红外区开始并延伸至带宽的宽广峰主导。当 $\bar{U}$ 较大时，低能峰和近红外峰之间出现两次符号变化。
    • 符号变化的解释： 被认为是反铁磁长程有序和相关费米面重构（特别是电子口袋出现）的前兆。在本文参数选择下，重构费米面上的空穴口袋通常主导电子口袋，因此低能峰通常为正。这与一些早期工作中电子口袋主导空穴口袋从而导致低能峰可能为负的情况不同。
    • 直流霍尔响应的微妙性： 结果暗示了低频霍尔响应非常微妙，可能敏感地依赖于系统低能物理的细节。

  2. 非高斯作用量（NLSM）结果 (图4)：

    • 基本的趋势与高斯作用量相似，双组分结构在NLSM作用量下表现得更为清晰。

  3. 霍尔电导率实部 ($\text{Re}\sigma_{xy}(\omega)$) (图5)：

    • 同样呈现双组分结构。当 $\bar{U}$ 较大时，观察到三次符号反转，其中前两次发生在低能峰和近红外峰之间，第三次发生在近红外峰之上。

2.2.3 反铁磁有序态中的光学和霍尔电导率 (图6, 7)

为了更好地理解涨落局部矩模拟结果中的双组分结构和符号变化，作者在平均场近似下计算了静态反铁磁有序态（SDW态）下的光学和霍尔电导率。SDW态具有一个SDW能隙 $\Delta_{AF}$，导致电子能带分裂为上下SDW带，并发生费米面重构（空穴口袋和电子口袋）。

  • 光学电导率 ($\sigma_{xx}(\omega)$) (图7a)： 清楚地显示出双组分结构。
    • 低能Drude峰：归因于SDW带内的带内跃迁。
    • 近红外峰：归因于跨SDW能隙的带间跃迁（起始能量约为 $2\Delta_{AF}$）。
    • 结论： 这表明在涨落局部矩模拟中观察到的双组分结构是反铁磁长程有序的先驱效应，可以视为SDW态的“涂抹”版本。
  • 霍尔电导率 ($\text{Im}\sigma_{xy}(\omega)$) (图7b)： 也显示出双组分结构和符号变化。
    • 低能峰：归因于带内跃迁。在本文参数下，空穴口袋主导电子口袋，因此低能峰为正。
    • 近红外峰：归因于带间跃迁。
    • 结论： 低能峰以下的 $\text{Im}\sigma_{xy}(\omega)$ 符号变化被理解为费米面重构（特别是贡献负霍尔响应的电子口袋出现）的前兆。这与电子掺杂铜氧化物超导体中的实验观测一致。

2.3 性能数据

论文中并未提供具体的计算性能数据（例如运行时间、CPU/GPU使用量、内存消耗等）。但是，从其描述中可以推断出一些关于性能的特点：

  • 计算量大： 对于每个蒙特卡罗样本，都需要计算费米子格林函数及其迹。这通常涉及对 $L \times L$ 矩阵（$400 \times 400$）的对角化或求逆，对于每个 $\phi$ 场配置，这都是一个计算密集型操作。
  • 采样数量： 6400个统计独立样本的数量，对于蒙特卡罗模拟来说是一个中等规模的样本量，表明单次模拟可以在合理的时间内完成。
  • 算法选择： 高斯作用量可以通过直接采样多元高斯分布进行，而NLSM作用量则采用热浴算法和过松弛技巧。这些都是标准且经过优化的蒙特卡罗采样方法。
  • 无符号问题优势： 避免符号问题是性能上的巨大优势，因为这使得蒙特卡罗采样能够有效探索配置空间，而不需要复杂的重加权或退火技术，否则这些技术会指数级地增加计算成本。
  • 热涨落极限简化： 聚焦于热涨落主导的极限，即忽略 $\phi$ 场的时间依赖性，极大地简化了格林函数计算，因为它们在Matsubara频率上变得对角化。这避免了昂贵的非对角计算和数值解析延拓。

尽管缺乏具体的性能指标，但该方法的可行性和所获得的数据量表明，在现代计算资源下，这种基于蒙特卡罗的有效理论框架对于研究强关联系统电磁响应是高效且可管理的。

3.1 代码实现细节 (基于论文描述推断)

论文并未直接提供源代码或详细的伪代码，但根据其方法描述，可以推断出以下关键的代码实现细节：

  1. 哈伯德模型初始化：

    • 晶格设置： 创建 $L \times L$ (例如 $20 \times 20$) 的二维正方晶格，设置周期性边界条件。
    • Hamiltonian构建： 根据给定的跳跃积分 $t, t'$ 和化学势 $\mu$ 构建非相互作用的Hamiltonian。在有磁场时，需要通过Peierls替代修改跳跃项。
    • 自旋密度算符 $S_i$： 实现自旋密度算符的矩阵元，通常通过Pauli矩阵和创建/湮灭算符表示。

  2. 哈伯德-斯特拉托诺维奇 (HS) 场 $\phi$ 的蒙特卡罗采样：

    • 高斯作用量采样：
      • 逆磁化率矩阵 $$\chi^{-1}(q, \omega_n=0)$$ 的构建： 根据MMP形式和格点正则化后的 $(q-Q)^2$ 项，构建此矩阵。
      • 多变量高斯分布采样： 由于高斯作用量 $S_{\text{eff}}[\phi]$ 是 $\phi$ 的二次型，可以直接从其对应的多变量高斯分布中采样 $\phi$ 场配置。这通常涉及对协方差矩阵（$\chi$）的Cholesky分解或奇异值分解。
    • 非高斯作用量 (NLSM) 采样：
      • 热浴算法 (Heat Bath Algorithm)： 对于每个位点 $i$，根据给定其邻居 $\phi_{\delta}$ 的场配置，从条件概率分布 $P(\phi_i | \{\phi_{\delta}\})$ 中重新采样 $\phi_i$。由于约束 $\phi_i \cdot \phi_i = 1$，采样的 $\phi_i$ 必须保持单位长度。
      • 过松弛技巧 (Over-relaxation Trick)： 为了提高蒙特卡罗链的混合效率，减少自相关时间，可以结合使用过松弛技巧。这通常涉及将 $\phi_i$ 绕某个轴旋转，以保持能量不变但改变其方向。
      • 关联长度调节： 通过调整参数 $\eta$ 来改变NLSM中的关联长度。
    • 样本管理： 记录和存储多个统计独立的 $\phi$ 场配置样本 (例如6400个)。

  3. 费米子格林函数和迹的计算：

    • 单粒子Hamiltonian构建： 对于每个采样的 $\phi$ 场配置和给定的电磁场 $A$，构建包含电子跳跃项和与 $\phi$ 场耦合项的单粒子Hamiltonian矩阵 $H_{\text{eff}}[\phi, A]$。
    • 对角化： 对 $H_{\text{eff}}[\phi, A]$ 进行对角化，得到本征能量 $E_m$ 和本征波函数 $\psi_{m,\sigma}$。
    • 费米分布函数： 根据温度 $T$ 和化学势 $\mu$ 计算费米分布函数 $f(E_m) = 1/(e^{(E_m - \mu)/k_BT} + 1)$。
    • 迹计算： 利用本征能量和波函数，以及费米分布函数，计算表达式中涉及格林函数 $G$ 和电磁顶点 $v^{(n)}$ 的各种迹。这些迹通常涉及形如 $\sum_{m,m'} \text{Tr}[\psi_m \psi_{m'}^{\dagger} \dots]$ 的求和。例如，Eqs. 63-66给出了迹的通用形式。
    • 电磁顶点 $v^{(n)}$： 根据Peierls替代，计算电磁顶点算符的矩阵元，如Eq. 34-36所示。

  4. 电导率计算和平均：

    • 光学电导率 ($\sigma_{xx}(\omega)$)： 根据Eq. 68，计算每个 $\phi$ 场样本的 $\kappa^{xx}$，然后对所有样本进行系综平均。$\delta$ 函数用洛伦兹函数替换进行展宽。
    • 霍尔电导率 ($\sigma_{xy}(\omega)$)：
      • 计算 $\kappa^{xyy}(q, \omega)$： 根据Eq. 62计算每个 $\phi$ 场样本的 $\kappa^{xyy}(q, \omega)$，然后进行系综平均。
      • 有限差分近似： 使用Eq. 77或78的有限差分近似计算关于 $q_x$ 的导数。
      • 反对称化： 为了恢复霍尔电导率的反对称性，可能需要手动对结果进行反对称化处理。
    • 实时频率转换： 由于专注于热涨落极限，HS场是时间无关的，可以直接在实频率下进行计算，避免了数值解析延拓。

  5. 数据后处理和可视化：

    • 将计算得到的电导率数据保存，并使用如matplotlib等工具进行绘图。

3.2 复现指南

复现这项工作需要对强关联物理、量子场论、蒙特卡罗方法和数值线性代数有深入理解。以下是复现的关键步骤和建议：

  1. 环境设置：

    • 编程语言： C++ 或 Python（结合NumPy/SciPy进行高性能数值计算）。C++通常用于性能关键的蒙特卡罗模拟。
    • 库： 线性代数库 (如Eigen, BLAS/LAPACK, SciPy.linalg)，用于矩阵对角化、求逆和高斯分布采样。
    • 计算集群： 考虑到 $L=20$ 晶格的Hamiltonian矩阵大小 ($400 \times 400$ 或更大，取决于内部自由度)，以及6400个样本的需求，高性能计算资源（多核CPU或GPU）将是必要的。

  2. 核心模块实现：

    • 晶格生成器： 用于构建周期性边界条件下的二维正方晶格。
    • Hamiltonian构造器： 实现根据参数 $t, t', \mu, \bar{U}, \phi$ 和电磁势 $A$ 构建费米子Hamiltonian。
    • 对角化器： 调用BLAS/LAPACK或Eigen等库进行Hamiltonian矩阵的对角化。
    • 蒙特卡罗采样器： 实现高斯场采样或NLSM的热浴/过松弛算法。
    • 迹计算器： 实现Eq. 63-66中各种迹的计算。
    • 电导率计算器： 将上述部件组合起来，计算 $\sigma_{xx}(\omega)$ 和 $\sigma_{xy}(\omega)$。

  3. 参数校准与验证：

    • 非相互作用极限： 首先在 $U=0$ 的非相互作用极限下验证代码，确保费米面、能带结构和电导率与已知结果一致。
    • 高斯极限： 验证高斯作用量采样和计算结果的稳定性。
    • 收敛性测试： 运行蒙特卡罗模拟时，确保体系达到平衡，并生成足够多的统计独立样本，以保证结果的收敛性。
    • 有限尺寸效应： 尝试不同 $L$ 值来评估有限尺寸效应，尽管本文只用了 $L=20$。

  4. 调试：

    • 数值稳定性： 注意浮点精度问题，特别是涉及到接近奇异的矩阵操作时。
    • 边界条件处理： 仔细处理周期性边界条件，避免错误索引。

论文中没有明确提及所使用的具体软件包或任何开源代码库链接。这是学术论文中常见的现象，特别是在方法细节较多的计算物理研究中。作者通常会自行实现核心算法。

然而，基于论文描述的计算任务，可以推断出可能需要或使用的通用数值库：

  • 高性能线性代数库：

    • LAPACK (Linear Algebra PACKage)： Fortran编写，提供了用于解决线性代数问题的例程（如矩阵对角化、特征值分解、求解线性方程组）。
    • BLAS (Basic Linear Algebra Subprograms)： 提供了矩阵-矩阵乘法、向量点积等基本线性代数操作，通常是LAPACK等高级库的底层优化实现。
    • Eigen (C++模板库)： 如果使用C++，Eigen是一个流行的、高性能的线性代数模板库，易于使用且功能强大。
    • NumPy/SciPy (Python库)： 如果使用Python，NumPy提供高效的数组操作，SciPy则提供了包括线性代数在内的许多科学计算功能，可以直接调用底层的BLAS/LAPACK。

  • 蒙特卡罗采样：

    • 对于高斯场，可以直接使用NumPy或SciPy的随机数生成器生成多变量高斯分布的样本。
    • 对于NLSM的热浴/过松弛，需要自定义实现，但可以利用Python的随机数模块。

  • 数据分析与绘图：

    • Python (Matplotlib/Seaborn)： 最常用的科学绘图库。
    • Gnuplot： 经典的命令行绘图工具。

总结： 尽管没有直接的开源repo链接，但复现者可以利用上述标准的数值计算库和编程工具，根据论文中详尽的方法描述自行实现。这是一项复杂的任务，需要投入大量时间和专业知识。

4. 关键引用文献，以及对这项工作局限性的评论

4.1 关键引用文献及其贡献

该论文引用了大量文献，这些文献为本研究的理论框架、数值方法以及对铜氧化物超导体物理的理解提供了重要背景。以下是一些关键文献及其贡献：

  1. [16] A. Millis, H. Monien, and D. Pines, Phenomenological model of nuclear relaxation in the normal state of YBa2Cu3O7, Phys. Rev. B 42, 197(1990).

    • 贡献： 这篇论文提出了著名的Millis-Monien-Pines (MMP) 唯象磁化率形式。本研究的关键突破之一就是，当局部矩涨落的有效作用量采用这种MMP形式时，量子蒙特卡罗模拟可以避免符号问题。这使得在强关联系统中对涨落局部矩进行非扰动研究成为可能，而MMP形式本身也被广泛用于描述铜氧化物超导体中的自旋涨落。

  2. [17] A. Abanov, A. V. Chubukov and J. Schmalian, Quantum-critical theory of the spin-fermion model and its application to cuprates: Normal state analysis, Advances in Physics 52, 117(2003).

    • 贡献： 这是一篇关于自旋-费米子模型的综述性文章，详细讨论了其量子临界理论及其在铜氧化物超导体正常态分析中的应用。自旋-费米子模型是本文有效理论描述的基础，它将电子视为巡游准粒子并与局部矩耦合。这篇综述为理解局部矩在铜氧化物超导体物理中的作用提供了理论背景。

  3. [20] J. Lin and A. Millis, Optical and Hall conductivities of a thermally disordered two-dimensional spin-density wave: Two-particle response in the pseudogap regime of electron-doped high-Tc superconductors, Phys. Rev. B 83, 125108(2011).

    • 贡献： 本文直接引用了Lin和Millis的这项工作，它也研究了电子掺杂铜氧化物超导体的光学和霍尔电导率在热涨落自旋密度波背景下的行为。然而，Lin和Millis的扰动处理在低频极限下出现了非物理性发散。本研究的重要意义在于，其提出的框架能够克服这种发散，并在任意频率下产生数值准确的预测，从而更可靠地揭示电导率的双组分结构。

  4. [21] A. Eberlein, W. Metzner, S. Sachdev, and H. Yamase, Fermi surface reconstruction and drop in the Hall number due to spiral antiferromagnetism in high-Tc cuprates, Phys. Rev. Lett. 117, 187001 (2016).

    • 贡献： 这篇论文通过假设赝能隙相中存在非零螺旋磁序，解释了铜氧化物超导体中霍尔数随掺杂的变化。本文借用了其平均场计算结果来解释其蒙特卡罗模拟中霍尔电导率符号变化的物理起源，将其与费米面重构（尤其是电子口袋的出现）联系起来。这为理解涨落局部矩效应如何导致类似平均场有序态下的谱特征提供了重要联系。

  5. [25] E. Berg, M. A. Metlitski and S. Sachdev, Sign-problem-free quantum Monte Carlo of the onset of antiferromagnetism in metals, Science 338, 1606(2012).

    • 贡献： 这篇Science论文展示了如何通过修改自旋-费米子模型，使其在反铁磁量子相变附近实现无符号问题的量子蒙特卡罗模拟。虽然他们的修改是为了研究量子临界行为而非输运性质，且采用了不同的费米子自由度处理，但它启发了本研究在特定条件下实现无符号问题蒙特卡罗的思路，特别是关于有效作用量形式的选择。

4.2 对这项工作的局限性评论

尽管本研究提出了一个强大的框架并取得了重要进展，但其也存在一些局限性，值得在未来的工作中加以解决：

  1. 量子涨落的忽略：

    • 局限： 论文主要关注热涨落主导的极限，即假设局部矩场 $\phi$ 是时间无关的。作者明确承认，这种简化导致预测的Drude权重相对于实验观测值偏弱，尤其是在强耦合区域。对于低频电磁响应，忽略局部矩的量子性质（即其时间动力学）会变得不准确。
    • 评论： 铜氧化物超导体中，特别是空穴掺杂的体系，自旋涨落是高度动态的。量子涨落预计会部分恢复费米能级附近的电子态密度，从而增强Drude峰的谱重。因此，当前模型在描述低频行为和与实验进行定量比较时存在固有不足。虽然文章提及未来将考虑量子效应，但这是当前研究的一个主要简化。

  2. 直流（DC）极限的复杂性：

    • 局限： 本文选择避免直接计算DC霍尔响应，而是专注于光学频率区域。作者指出DC霍尔响应的理论计算和分析更为复杂，且强烈依赖于系统低能物理的细节（如费米面几何、准粒子寿命的动量依赖性）。此外，在低频极限下，热涨落主导的假设最终会失效。
    • 评论： 尽管光学霍尔响应具有一定的普遍趋势，但DC输运是实验上最容易获得的量，其物理意义也至关重要。回避DC极限虽然简化了计算，但限制了模型与一些最直接实验数据的比较能力。

  3. 唯象磁化率的依赖：

    • 局限： 论文依赖MMP形式的唯象磁化率来确保蒙特卡罗模拟的无符号问题特性。虽然MMP形式在铜氧化物研究中被广泛接受，但它并非从第一性原理在强耦合机制下严格推导出来的。
    • 评论： 采用唯象参数有助于克服计算障碍，但可能限制了模型从更基本层面揭示物理机制的能力。未来的工作可能需要探索如何在不引入唯象参数的情况下，从更微观的理论中导出无符号问题的有效作用量。

  4. 数值近似的引入：

    • 局限： 霍尔电导率的计算中使用了有限差分近似来处理 $q_x \to 0$ 极限下的导数。这种近似在有限尺寸晶格上可能引入微小的数值伪影，例如在 $\omega = 0$ 附近出现的微弱负过冲。
    • 评论： 尽管这些伪影可能很小，但它们表明了数值方法本身可能存在的限制。更精确的导数计算方法（例如，通过对 $q_x$ 的多项式拟合）可能有助于消除这些误差。

  5. 有限尺寸效应：

    • 局限： 模拟在 $L=20$ 的有限晶格上进行。对于某些输运性质，特别是当关联长度 $\xi$ 变长时，这种尺寸可能不足以完全消除有限尺寸效应。
    • 评论： 虽然 $20 \times 20$ 晶格对于许多蒙特卡罗模拟来说是合理的，但对于精确捕捉长程关联或在热力学极限下的行为可能需要更大的系统尺寸，这会显著增加计算成本。

  6. 定性而非定量研究：

    • 局限： 作者明确表示，本研究旨在对热自旋涨落的影响进行定性讨论，而非对实验数据进行定量拟合。
    • 评论： 这种选择虽然合理，但在直接与实验结果进行比较时，其解释力会受到限制。未来工作如果能结合更全面的参数空间探索和量子效应，有望实现更强的定量预测能力。

总的来说，这项工作通过引入无符号问题的蒙特卡罗框架，为研究强关联电子系统的电磁响应提供了一条有前途的途径。尽管存在上述局限，但它为未来的研究奠定了坚实的基础，特别是在将量子涨落和更复杂的低能物理纳入考量方面。

5. 其他你认为必要的补充

5.1 铜氧化物超导体背景下的非费米液体行为

铜氧化物超导体以其在高临界温度下表现出的奇异金属（strange metal）行为而闻名，其正常态物理远非传统的费米液体理论所能描述。典型的非费米液体特征包括：

  • 线性电阻率： 在宽广的温度范围内，电阻率随温度呈线性关系，而非费米液体预期的 $T^2$ 关系。
  • 非Drude光学电导率： 如论文所指出的，光学吸收谱不服从简单的Drude模型，而是展现出低能Drude峰和高能近红外峰的双组分结构，谱重分布覆盖整个能带宽度。
  • 奇异的霍尔响应： 霍尔数 $n_H$ 随温度的复杂依赖性，以及在掺杂相图中观察到的费米面重构迹象，但缺乏对称性破缺的证据。

这些现象强烈指向强关联效应在这些材料中的主导作用。掺杂的Mott绝缘体是描述铜氧化物超导体的一个核心概念，其中电子-电子相互作用非常强，导致电子局域化并形成局部磁矩。当体系被掺杂（无论是空穴还是电子）时，这些局部矩会与新引入的巡游载流子相互作用，形成一个复杂的系统，具有巡游准粒子和局部矩的双重性质。本研究正是建立在这种理解之上，试图通过涨落局部矩的有效理论来解释这些非费米液体行为。

5.2 局部矩涨落与自旋-费米子模型

局部矩涨落是强关联Mott绝缘体衍生体系中的一个核心物理图像。在铜氧化物中，即使在宏观磁有序消失的区域，短程反铁磁关联仍然存在，并表现为动态的自旋涨落（或称准磁子）。RIXS（共振非弹性X射线散射）实验已证实，准磁子激发在铜氧化物相图中普遍存在。这些涨落的局部矩被认为是理解铜氧化物物理的关键要素，不仅影响输运性质，还可能在超导配对机制中发挥作用。

自旋-费米子模型是处理这种双重性质的有效工具。它将电子视为在费米面附近巡游的准粒子，并与具有特定动量和频率的集体自旋涨落模式耦合。该模型在解释d波超导配对和正常态异常行为方面取得了成功。然而，传统自旋-费米子模型的一个主要挑战是顶点修正问题，这使得计算输运性质变得复杂且易出现非物理发散。本文的贡献在于，它提供了一个在非扰动水平上处理局部矩涨落的框架，有效地解决了这些问题，特别是通过其无符号蒙特卡罗模拟方法。

5.3 无符号问题量子蒙特卡罗的突破

量子蒙特卡罗（QMC）是研究强关联系统最强大的非扰动工具之一。然而，著名的“符号问题”长期以来一直是其主要障碍，使得无法直接模拟费米子系统在许多物理相关条件下的基态和有限温度性质。本研究通过识别出在特定条件下（特别是采用MMP型唯象磁化率或在高温极限下）可以避免符号问题，为QMC在强关联输运问题中的应用开辟了新途径。这一突破意义重大：

  • 可靠的非扰动结果： 摆脱了扰动理论的限制，可以获得更可靠的结果，尤其是在强耦合区域。
  • 处理顶点修正： 自动包含了复杂的顶点修正效应，而无需手动处理。
  • 更广泛的应用前景： 这种方法不仅限于本研究中的电磁响应，还可以扩展到其他集体涨落（如电荷密度波、配对涨落）和更广泛的强关联电子系统。

5.4 与实验观测的联系和未来方向

本研究的定性结果与铜氧化物超导体的实验观测具有良好的吻合：

  • 双组分电导率： 预测的光学和霍尔电导率的双组分结构与实验中普遍观察到的Drude峰和近红外峰相符。
  • 霍尔响应的符号变化： 霍尔电导率中的符号变化被解释为费米面重构的先兆，特别是电子口袋的出现，这与电子掺杂铜氧化物超导体中的一些实验观察一致。

尽管当前工作主要侧重于热涨落主导的极限，但该框架的通用性为未来的研究提供了明确的方向：

  1. 纳入量子涨落： 这是最重要的未来工作，特别是对于描述低频电磁响应和恢复更强的Drude权重。通过在蒙特卡罗模拟中包含 $\phi$ 场的时间依赖性，可以更完整地捕捉局部矩的动力学。
  2. 直流输运性质： 结合量子涨落，研究DC极限下的纵向和霍尔电导率，将模型与最直接的实验测量进行比较。
  3. 其他集体涨落： 将该框架推广到电荷密度波和超导配对涨落等其他集体模式，从而更全面地理解铜氧化物超导体的复杂相图。
  4. 单粒子性质： 将电子与Grassmannian外部场耦合，计算系统的单粒子性质（如谱函数），以研究赝能隙的起源和费米弧等现象。

这项研究不仅为理解铜氧化物超导体的非费米液体行为提供了一个新的强大工具，也为解决更广泛的强关联物理问题开辟了新的途径。通过逐步克服现有模型的局限性，有望最终揭示高温超导的深层机制。