[论文解读] Numerical Studies of the 2D Hubbard Model

D. J. Scalapino|arXiv (Cornell University)|Oct 25, 2006

Physics of Superconductivity and Magnetism被引用 8

一句话总结

本文综述了数值证据，表明二维 Hubbard 模型能够捕捉铜氧化物超导体的关键物理现象，包括反铁磁性、d波配对、条纹态以及赝能隙。通过使用行列式量子蒙特卡洛、动态簇近似和 DMRG 方法，研究识别出有效粒子-粒子相互作用 Γ⁠pp 是配对的主要媒介，揭示了由 U、t′ 和掺杂浓度驱动的各相之间微妙的竞争关系。

ABSTRACT

Numerical studies of the two-dimensional Hubbard model have shown that it exhibits the basic phenomena seen in the cuprate materials. At half-filling one finds an antiferromagnetic Mott-Hubbard groundstate. When it is doped, a pseudogap appears and at low temperature d-wave pairing and striped states are seen. In addition, there is a delicate balance between these various phases. Here we review evidence for this and then discuss what numerical studies tell us about the structure of the interaction which is responsible for pairing in this model.

研究动机与目标

• 确定二维 Hubbard 模型是否在真实参数范围内重现了铜氧化物超导体的关键物理行为，包括反铁磁性、d 波配对、条纹态和赝能隙现象。
• 研究 Hubbard 模型中配对机制所依赖的有效相互作用的本质，重点关注不可约粒子-粒子顶点 Γ⁠pp。
• 通过分析边界条件、晶格尺寸和数值偏差对相竞争的影响，解决不同方法（如 DMRG 与 QMC）之间结果不一致的问题。
• 建立一个框架，将数值结果与实验可观测量（如频率依赖的能隙和重正化函数）关联，以识别高温超导体中的配对机制。
• 探讨是否存在一种共同机制（如自旋涨落或量子临界性）解释铜氧化物、重费米子材料和锕系元素中的非传统超导性。

提出的方法

• 采用行列式量子蒙特卡洛方法计算单粒子和双粒子格林函数，并提取不可约粒子-粒子顶点 Γ⁠pp。
• 应用动态簇近似（DCA）研究具有受控边界条件和自能反馈的有限尺寸簇，以获取低能物理性质。
• 在 n 通道梯状结构上使用密度矩阵重正化群（DMRG）研究有限掺杂下的基态，并分析在不同 U/t 和 t′ 条件下条纹态的形成。
• 分析 Γ⁠pp 的动量和 Matsubara 频率依赖性，以判断配对是否由自旋涨落（S=1）、电荷涨落（S=0）或更复杂的机制介导。
• 对比不同方法的结果，评估数值偏差以及纵横比和开放边界条件对相稳定性的影响。
• 利用提取出的 Γ⁠pp 推断有效配对相互作用的结构，重点关注其短程特性及其在 Tc 以上的主导作用。

实验结果

研究问题

• RQ1 在真实参数范围内，二维 Hubbard 模型是否表现出反铁磁性、d 波超导配对、条纹序和赝能隙现象？
• RQ2 Hubbard 模型中有效配对相互作用 Γ⁠pp 的结构是什么？其主要由自旋、电荷还是其他涨落介导？
• RQ3 为何不同数值方法（如 DMRG 与 QMC）在相似模型参数下对基态（条纹态 vs. d 波配对）的结论存在冲突？
• RQ4 边界条件和有限尺寸效应对数值模拟中竞争相（如条纹态和 d 波配对）的稳定性有何影响？
• RQ5 Hubbard 模型中的有效相互作用 Γ⁠pp 是否可与实验可观测量（如频率依赖的能隙 Δ(ω) 和重正化函数 Z(ω)）关联？

主要发现

• 二维 Hubbard 模型再现了关键的铜氧化物物理现象：在半满时表现为反铁磁性，在欠掺杂区域表现出赝能隙行为，以及在掺杂后出现 d 波配对和条纹序。
• 数值研究表明，反铁磁相、条纹相和 d 波超导相之间存在微妙的平衡，U、t′ 或掺杂浓度的微小变化即可导致系统在各相之间转换。
• 在掺杂的 8 根梯状 t-J 模型中，DMRG 计算得到条纹态基态，而变分法和格林函数蒙特卡洛在 t-J 晶格上的结果则更倾向于 d 波超导序，仅表现出微弱的条纹特征。
• 对于 Hubbard 梯状结构，当 U/t > 3 时条纹态出现；然而，簇 DMRG 和簇动态平均场研究显示，d 波配对与反铁磁性在相同参数区域内相互竞争。
• 行列式量子蒙特卡洛和簇计算均证实了在欠掺杂区域存在赝能隙行为，与铜氧化物中的实验观测一致。
• 有效配对相互作用 Γ⁠pp 具有短程特性，并在 Tc 以上占主导地位，使其可在较小簇上进行研究；其动量和频率依赖性揭示了以 S=1（磁性）交换机制为主导的配对机制。