[论文解读] Introduction to Hubbard Model and Exact Diagonalization

S. A. Jafari|ArXiv.org|Jul 30, 2008

Magnetism in coordination complexes参考文献 5被引用 15

一句话总结

本文介紹了 Hubbard 模型和精确对角化（ED）作为研究强关联电子体系的数值精确方法。它展示了如何通过在对称性自适应基底下对哈密顿量进行对角化，使 ED 能够在小簇上实现精确求解，揭示了自旋单重态基态等关键物理现象，并在大 U 极限下展现出 t-J 模型的涌现，其中自旋交换相互作用 J ≈ 4t²/U。

ABSTRACT

Hubbard model is an important model in theory of strongly correlated electron systems. In this contribution we introduce this model along with numerically exact method of diagonalization of the model.

研究动机与目标

• 为 Hubbard 模型及其在强关联电子物理中的作用提供教学性介绍。
• 将精确对角化（ED）作为在小簇上求解 Hubbard 模型的数值精确方法进行展示。
• 演示如何构建对称性自适应基底，以简化哈密顿量矩阵的表示。
• 通过 ED 结果说明自旋物理以及在大 U 极限下 t-J 模型的涌现。
• 通过一个简单的 Fortran 程序引导读者实现 ED，以供教育目的。

提出的方法

• 在第二量子化形式下形式化 Hubbard 哈密顿量，包含跃迁项（−t）和原子内排斥项（U），使用晶格上的产生和湮灭算符。
• 对非相互作用部分（H₀）应用傅里叶变换，以对角化并获得余弦能带色散关系 εₖ = −2t cos k。
• 利用群表示理论，通过投影算符 Γₚ(R)*Rϕ 构造对称性自适应的基态，确保具有确定的空间宇称和自旋量子数。
• 在对称性自适应基底下表示完整哈密顿量，将其分块对角化为具有确定 S_z 和空间对称性的子空间。
• 对每个子块求解本征值问题，以获得精确的能级和波函数，特别是在 S_z = 0 子空间中。
• 通过分析自旋单重态与三重态之间的能级分裂，在大 U 极限下推导出有效 t-J 模型，得到 J ≈ 4t²/U。

实验结果

研究问题

• RQ1 如何通过数值对角化精确求解小簇中的 Hubbard 模型？
• RQ2 对称性自适应在简化哈密顿量矩阵和提高计算效率方面起到什么作用？
• RQ3 两体 Hubbard 模型的基态性质是什么？其依赖于 U 和 t 的行为如何？
• RQ4 Hubbard 模型在大 U 极限下的行为如何映射到 t-J 模型？有效交换相互作用 J 是什么？
• RQ5 在两体模型中，自旋单重态与三重态之间的能量分裂是多少？其随 t 和 U 的变化规律如何？

主要发现

• 两体 Hubbard 模型的基态为自旋单重态，能量为 Eₛ⁻ = U/2 − √(U²/4 + 4t²)，低于三重态能量 Eₜ = 0。
• 在大 U 极限下，自旋单重态相对于三重态的能量分裂约为 −J = −4t²/U，表明存在反铁磁关联。
• 对称性自适应基底将哈密顿量简化为分块对角形式，其中 S_z = 0 子空间包含一个 2×2 的单重态块和一个 2×2 的三重态块，从而实现精确对角化。
• 三重态（S_z = 0）处于零能级，而双占态 |ψ₅⟩ 的能量为 U，证实了在大 U 区域电荷涨落被抑制。
• 在大 U 极限下，单重态波函数主要由配置 |ψ₃⟩ = (|ϕ₃⟩ + |ϕ₄⟩)/√2 主导，该配置无双占，确认了电荷涨落的抑制。
• 在大 U 极限下，有效 t-J 模型涌现，其中低能物理由自旋交换主导，J ≈ 4t²/U，且希尔伯特空间被限制在无双占的态中。