2025-11-19T17:58:14.286609

Charge Susceptibility and Kubo Response in Hatsugai-Kohmoto-related Models

Ma, Zhao, Huang et al.
We study in depth the charge susceptibility for the band Hatsugai-Kohmoto (HK) and orbital (OHK) models. As either of these models describes a Mott insulator, the charge susceptibility takes on the form of a modified density response function with lower and upper Hubbard bands, thereby giving rise to a multi-pole structure. The particle-hole continuum consists of hot spots along the $ω$ vs $q$ axis arising from inter-band transitions. Such transitions, which are strongly suppressed in non-interacting systems, obtain here because of the non-rigidity of the Hubbard bands. This modified density response function gives rise to a plasmon dispersion that is inversely dependent on the momentum, resulting in an additional contribution to the conventional f-sum rule. This extra contribution originates from a long-range diamagnetic contribution to the current. This results in a non-commutativity of the long-wavelength ($q\rightarrow 0$) and thermodynamic ($L\rightarrow\infty$) limits. When the correct limits are taken, we find that the Kubo response computed with either open or periodic boundary conditions yields identical results that are consistent with the continuity equation contrary to recent claims. We also show that the long wavelength pathology of the current noted previously also plagues the Anderson impurity model interpretation of dynamical mean-field theory (DMFT).
academic

Charge Susceptibility and Kubo Response in Hatsugai-Kohmoto-related Models

基本信息

  • 论文ID: 2409.07522
  • 标题: Charge Susceptibility and Kubo Response in Hatsugai-Kohmoto-related Models
  • 作者: Yuhao Ma, Jinchao Zhao, Edwin W. Huang, Dhruv Kush, Barry Bradlyn, Philip W. Phillips
  • 分类: cond-mat.str-el (强关联电子系统)
  • 发表时间: 2024年9月 (arXiv预印本)
  • 论文链接: https://arxiv.org/abs/2409.07522

摘要

本文深入研究了带状Hatsugai-Kohmoto (HK)模型和轨道HK (OHK)模型的电荷磁化率。由于这两个模型都描述了Mott绝缘体，其电荷磁化率呈现为具有下、上Hubbard带的修正密度响应函数形式，从而产生多极结构。粒子-空穴连续谱在ω vs q轴上包含由带间跃迁产生的热点。这些跃迁在非相互作用系统中被强烈抑制，但在此处由于Hubbard带的非刚性而得以实现。这种修正的密度响应函数产生了与动量成反比的等离子体色散，导致对传统f-sum规则的额外贡献。这一额外贡献源于电流的长程反磁性贡献，导致长波长(q→0)和热力学(L→∞)极限的非对易性。

研究背景与动机

问题背景

  1. 强关联电子系统的挑战: 在强关联电子系统(如铜氧化物超导体)中，传统的Fermi液体理论失效，需要新的理论框架来理解其电荷响应性质。
  2. Mott绝缘体的特殊性: Mott绝缘体由于强电子关联效应，其能带结构分裂为上、下Hubbard带，这与传统的能带理论预期大相径庭。
  3. 实验观测的异常现象: 最新的动量分辨电子能量损失谱(MEELS)实验在铜氧化物中观察到：
    • 中等能量尺度上的平坦温度和动量无关连续谱
    • 低能时χ''∝q²，高能时衰减为q²/ω²
    • 与密度无关的等离子体模式

研究动机

  1. 理论空白: 尽管HK模型作为精确可解的Mott物理平台被广泛研究，但其密度-密度响应函数此前从未被深入研究。
  2. 基础理论问题: 需要理解在强关联系统中：
    • 粒子-空穴连续谱如何在非刚性Hubbard带存在下发生变化
    • 导致等离子体频率的动量因子抵消是否仍然成立
    • f-sum规则是否需要修正

核心贡献

  1. 首次完整计算了HK模型的电荷磁化率: 提供了带状HK和轨道HK模型密度响应函数的精确解析表达式。
  2. 发现了异常的等离子体色散关系: 证明了HK模型中等离子体频率ωₚ∼1/q，与传统金属中的行为截然不同。
  3. 揭示了修正的f-sum规则: 发现由于长程反磁性电流贡献，f-sum规则获得额外项，破坏了传统的q²标度律。
  4. 解决了电流算符定义的争议: 通过最小耦合原理严格推导了规范不变的电流算符，证明了连续性方程的成立。
  5. 建立了与DMFT的联系: 证明了单杂质Anderson模型在DMFT框架下也存在类似的长波长病理问题。

方法详解

任务定义

研究HK相关模型中的电荷磁化率χ(q,ω)，其中：

  • 输入: 动量q和频率ω
  • 输出: 密度-密度关联函数
  • 约束: 满足电荷守恒和连续性方程

模型架构

1. 带状HK模型(bHK)

哈密顿量为：

H_bHK = Σ_kσ ξ_k n_kσ + Σ_k U n_k↑ n_k↓

其中ξ_k = ε_k - μ是动能项，U是动量无关的密度-密度排斥相互作用。

2. 轨道HK模型(OHK)

H_OHK = Σ_k,α,α',σ g_αα'(k) c†_kασ c_kα'σ - μΣ_k,α n_kα,σ + Σ_k,α,α' U_αα' n_kα↑ n_kα'↓

3. 密度响应函数计算

利用HK模型的动量因式分解性质：

χ⁰_HK(q,ω) = (1/N) Σ_kσ Σ_i,j=L,U [w^i_k+q w^j_k (f(ξ^j_k) - f(ξ^i_k+q))/(ω + i0⁺ - ξ^i_k+q + ξ^j_k)]

技术创新点

1. 多带结构的精确处理

  • 将HK的双极结构映射为有效的多带系统
  • 区分直接项(L,L)和(U,U)与交叉项(L,U)和(U,L)
  • 交叉项导致q⁰阶贡献，这是Mott系统特有的

2. 长波长极限的细致分析

发现了关键的非对易性：

lim_q→0 lim_L→∞ ≠ lim_L→∞ lim_q→0

这要求在有限尺寸系统中先取q→0极限，再取热力学极限。

3. 规范不变电流算符的构建

通过最小耦合原理：

U_αβ(R₁+R₃-R₂-R₄) → U^A_αβ(R₁+R₃-R₂-R₄) = U_αβ(R₁+R₃-R₂-R₄)e^(-i∫A·dx)

导出了满足连续性方程的电流密度算符。

实验设置

理论计算框架

  1. 近似条件:
    • 长波长高频极限：q·P_i/(m(ω+U)) << 1
    • 低温极限：∂f/∂ξ_k ≈ -δ(ξ_k)
    • 抛物色散：ξ_k = |k|²/2m - μ
  2. RPA处理:
    • HK相互作用精确求解
    • Coulomb相互作用在RPA近似下处理
    • 介电函数：ε_RPA(q,ω) = 1 - V_q χ⁰_HK(q,ω)

数值参数

  • 方格晶格紧束缚色散
  • t = 0.1eV, t_p = 0.03eV, U = 0.14eV
  • 系统尺寸和边界条件的系统性研究

实验结果

主要结果

1. 等离子体色散的异常行为

发现等离子体频率遵循：

ω_p ∼ 1/q

这与传统金属中的ω_p = constant形成鲜明对比。

2. 修正的f-sum规则

传统f-sum规则：

∫ dω ω Im χ(q,ω) = -π n_e q²/m

在HK模型中获得额外贡献：

∫ dω ω Im χ_HK(q,ω) = (4πU/N) Σ_k (⟨n_k↓n_k↑⟩ - ⟨n_k↓⟩⟨n_k↑⟩)

3. 粒子-空穴连续谱的热点结构

  • 在ω = U处出现显著的热点，对应带间跃迁
  • 与非相互作用系统中被抑制的带间跃迁形成对比
  • 热点强度与单占据区域体积成正比

消融实验

通过分离直接项和交叉项的贡献：

  • 直接项：χ_dir ∝ q²/ω²，与自由电子类似
  • 交叉项：χ_cro ∝ 1/ω²，提供主导贡献

边界条件的一致性验证

证明了开放边界条件和周期边界条件下：

  • 电流算符形式一致
  • 连续性方程均成立
  • 长波长极限行为相同

相关工作

HK模型的发展历程

  1. 原始HK模型(Hatsugai & Kohmoto, 1992): 提出了精确可解的强关联模型
  2. 轨道推广(近年来): 发展了多轨道版本，改善了对Hubbard模型的逼近
  3. 拓扑性质研究: 在拓扑Mott绝缘体等方面的应用

密度响应函数研究

  1. DCA数值研究(Dong et al.): 在虚频上的数值计算，需要解析延拓
  2. MEELS实验: 在铜氧化物中的最新观测结果
  3. 传统理论: RPA和Fermi液体理论的局限性

电流算符争议

最近Guerci等人质疑了HK模型中电流算符的定义，本文给出了完整的澄清。

结论与讨论

主要结论

  1. Mott绝缘体的独特响应: HK模型展现了与传统金属截然不同的电荷响应特征
  2. 带间跃迁的重要性: 非刚性Hubbard带使得带间跃迁不再被抑制
  3. 长程效应的关键作用: 长程反磁性电流导致了f-sum规则的修正
  4. 理论一致性: 正确处理极限顺序后，理论框架完全自洽

局限性

  1. RPA近似: 对Coulomb相互作用的处理仍在RPA水平
  2. 模型简化: HK模型虽然精确可解，但相对于真实材料仍有简化
  3. 维度限制: 主要结果在d=∞下严格成立，有限维度需要修正

未来方向

  1. 实验验证: 需要在真实材料中寻找类似的异常等离子体行为
  2. 有限维度效应: 研究维度对结果的影响
  3. 超导态: 扩展到掺杂Mott绝缘体的超导相
  4. 输运性质: 研究相关的电导率和霍尔效应

深度评价

优点

  1. 理论严谨性: 提供了HK模型电荷响应的首个完整理论框架
  2. 数学精确性: 所有计算都是精确的，没有不受控近似
  3. 物理洞察深刻: 揭示了Mott物理中的新现象和机制
  4. 争议澄清: 解决了电流算符定义的重要争议
  5. 写作清晰: 逻辑严密，表述清楚

不足

  1. 实验联系: 缺乏与具体实验的直接对比
  2. 数值验证: 某些解析结果缺乏独立的数值验证
  3. 适用范围: 对于真实材料的适用性需要进一步验证

影响力

  1. 理论贡献: 为强关联电子系统理论提供了新的理论工具
  2. 方法学价值: 处理长程相互作用系统的方法具有普适性
  3. 实验指导: 为相关实验提供了理论预期和解释框架

适用场景

  1. 强关联电子系统: 特别是Mott绝缘体和相关材料
  2. 理论模型研究: 作为研究强关联效应的重要工具模型
  3. 数值方法验证: 为DMFT等数值方法提供基准测试

参考文献

本文引用了48篇重要参考文献，涵盖了HK模型、强关联电子系统、DMFT理论等相关领域的关键工作，为研究提供了坚实的理论基础。

总体评价: 这是一篇在强关联电子系统理论方面具有重要贡献的高质量论文。作者通过精确的理论分析，揭示了Mott绝缘体中电荷响应的新奇性质，不仅解决了重要的理论争议，还为理解强关联材料的输运性质提供了新的理论框架。论文的技术水平高，物理洞察深刻，对该领域的发展具有重要推动作用。