2025-11-21T12:28:15.020679

Interplay of magnetic and thermodynamic responses in the kagome-triangular system

Jia, Zhang, Duan et al.
Inspired by the recent experimental progress in pyrochlore derivative \ce{RE3Sb3A2O14 (A=Mg, Zn)}, we investigate the Hubbard model on the kagome lattice with an additional hopping $t'/t$, which enables continuous interpolation between the kagome and triangular lattices by using determinant quantum Monte Carlo simulations. We analyze the evolution of magnetic correlations and thermodynamic responses across different values of $t'/t$ and on-site interaction $U$. It is found that increasing $t'/t$ suppresses short-range antiferromagnetic correlations, while the next-nearest-neighbor correlations exhibit a sign change near $t'/t \approx 0.3 \text{--} 0.4$. Within this regime, the specific heat shows a pronounced low-temperature peak, indicating an emergent spin-related energy scale. Increasing $U$ enhances magnetic correlations and shifts the associated $t'/t$ crossover points to larger values. We also discuss the sign problem to clarify which parameter region of our numerical simulations is accessible and reliable. Our results uncover the competition between frustration and correlations and the interplay of magnetic and thermodynamic responses in the kagome lattice, providing insights into correlated states in frustrated materials.
academic

Interplay of magnetic and thermodynamic responses in the kagome-triangular system

基本信息

  • 论文ID: 2510.13530
  • 标题: Interplay of magnetic and thermodynamic responses in the kagome-triangular system
  • 作者: Zixuan Jia, Lufeng Zhang, Qingzhuo Duan, Zenghui Fan, Jingyao Wang, Bing Huang, Tianxing Ma
  • 分类: cond-mat.str-el (强关联电子系统)
  • 发表时间: 2024年10月15日 (arXiv预印本)
  • 论文链接: https://arxiv.org/abs/2510.13530

摘要

本研究受到烧绿石衍生物RE₃Sb₃A₂O₁₄ (A=Mg, Zn)实验进展的启发,通过确定性量子蒙特卡罗(DQMC)模拟,研究了带有额外跳跃参数t'/t的kagome晶格上的Hubbard模型,该参数能够在kagome和三角晶格之间实现连续插值。研究分析了不同t'/t值和在位相互作用U下磁关联和热力学响应的演化。发现增加t'/t会抑制短程反铁磁关联,而次近邻关联在t'/t ≈ 0.3-0.4附近出现符号改变。在此区域内,比热显示出明显的低温峰,表明出现了与自旋相关的能量尺度。增加U会增强磁关联并将相关的t'/t交叉点移向更大的值。

研究背景与动机

问题背景

  1. Kagome晶格的独特性质:kagome晶格具有三个关键特征:高能边缘的平带、费米能级附近的狄拉克锥,以及由鞍点产生的多个van Hove奇点。这些特征导致了非平凡的电子拓扑、关联多体态和奇异磁性相。
  2. 实验动机:烧绿石衍生物RE₃Sb₃A₂O₁₄材料中观察到的RE/A位点混合无序不仅破坏了理想的晶格几何,还产生了有趣的物理现象,如Dy₃Sb₃Zn₂O₁₄中观察到的低温比热异常。
  3. 理论挑战:尽管对kagome晶格进行了广泛的理论研究,其基态性质、全局相图和竞争相互作用机制仍不完全清楚。

研究重要性

  • 理解几何阻挫和电子关联之间的竞争机制
  • 为阻挫材料中的关联态提供理论见解
  • 连接理论模型与实验观测的桥梁

核心贡献

  1. 提出了连续插值模型:通过引入可调跳跃参数t'/t,实现了kagome和三角晶格之间的连续过渡
  2. 发现了磁-热关联:揭示了磁关联和热力学响应之间的内在联系
  3. 识别了临界区域:在t'/t ≈ 0.3-0.4区域发现了次近邻关联的符号改变和比热低温峰的出现
  4. 系统分析了符号问题:为DQMC模拟的可靠性提供了参数依赖的指导
  5. 提供了实验解释:为理解实验中观察到的低温异常提供了理论框架

方法详解

模型定义

研究采用的Hubbard模型定义在二维kagome晶格上,并在d位点添加额外格点:

H=H1+H2+H3H = H_1 + H_2 + H_3

其中:

  • H1H_1:kagome晶格上的最近邻跳跃项,跳跃强度为t
  • H2H_2:kagome格点(a,b,c)与额外格点d之间的跳跃项,跳跃强度为t'
  • H3H_3:Hubbard相互作用项和化学势项

数值方法

采用确定性量子蒙特卡罗(DQMC)方法:

  1. 路径积分表示:将配分函数表示为虚时间离散化的路径积分
  2. Hubbard-Stratonovich变换:引入辅助场解耦相互作用项
  3. 蒙特卡罗采样:使用Metropolis算法更新辅助场
  4. 统计分析:8000次预热扫描 + 10000-30000次测量扫描

观测量计算

  1. 自旋-自旋关联函数cαγ(r)=13Sr0αSr0+rγc^{\alpha\gamma}(r) = \frac{1}{3}\langle S^{\alpha}_{r_0} \cdot S^{\gamma}_{r_0+r}\rangle
  2. 局域磁矩m2=(nr,nr,)2\langle m^2\rangle = \langle(n_{r,\uparrow} - n_{r,\downarrow})^2\rangle
  3. 比热c(T)=1NdHdTc(T) = \frac{1}{N}\frac{d\langle H\rangle}{dT}

实验设置

参数范围

  • 晶格尺寸:主要使用L = 6 (N = 4×L² = 144个格点),部分结果验证L = 8
  • 温度范围:T = t/10 到 T = t/4
  • 相互作用强度:U = 2.0t 到 U = 3.5t
  • 跳跃比值:t'/t = 0.0 到 1.0

技术细节

  • 虚时间步长∆τ选择足够小以确保Trotter误差小于统计误差
  • 使用10个bins进行粗粒化平均和误差估计
  • 通过平均符号⟨sign⟩监控符号问题的严重程度

实验结果

磁关联的演化

最近邻关联

  • 所有检验的U值下,最近邻b-c关联cbc(r=1)c_{bc}(r=1)保持负值,表明短程反铁磁关联持续存在
  • cbc(r=1)c_{bc}(r=1)的幅度随t'/t增加而减小,表明增大t'/t抑制短程反铁磁关联
  • 在固定t'/t下,更强的相互作用U增强了关联

次近邻关联

  • 次近邻c-c关联ccc(r=2)c_{cc}(r=2)在t'/t增加时出现符号改变
  • 交叉点随相互作用强度增强移向更大的t'/t值
  • U = 3.0t时,符号改变发生在t'/t ≈ 0.3-0.4

热力学性质

比热行为

  1. t'/t = 0时:显示与之前QMC研究一致的宽高温峰
  2. t'/t = 0.3时:出现软肩膀结构
  3. t'/t = 0.4时:低温峰变得清晰可见
  4. t'/t增大时:低温峰持续存在并逐渐移向高温

相互作用强度的影响

  • U = 2.0t:在t'/t = 0.3处已出现可辨识的低温特征
  • U增大:低温特征被抑制
  • U = 3.5t:低温特征消失

有限尺寸效应

通过L = 6和L = 8的对比确认:

  • 关联函数对系统尺寸无显著依赖
  • 符号改变的t'/t值在尺寸增大时基本不变
  • L = 6的计算结果收敛可靠

符号问题分析

  • 更大的U和更低的T加剧符号问题
  • U = 2.0t时可达到T = t/10的低温
  • U = 4.0t时在t'/t ≥ 0.6区域出现严重符号问题
  • 研究限制在⟨sign⟩ ≳ 0.1的可靠区域

相关工作

Kagome晶格的理论研究

  1. Heisenberg模型:DMRG、张量网络、变分蒙特卡罗等方法的研究结果存在分歧
  2. 扩展Heisenberg模型:J₁-J₂-J₃相互作用揭示了丰富的竞争相图
  3. Hubbard模型:捕捉电子巡游性和相互作用的竞争

热力学性质研究

  • Heisenberg模型中kagome和三角晶格都显示单峰结构
  • 动力学平均场理论预测kagome晶格出现小的次级峰
  • 连续插值研究揭示了双峰结构的可能性

实验材料

烧绿石衍生物RE₃Sb₃A₂O₁₄展现了丰富的磁性行为,包括长程磁序、标量手征磁序、短程自旋关联等。

结论与讨论

主要结论

  1. 磁-热关联机制:次近邻关联的符号改变与比热低温峰的出现在相同的t'/t区域,表明存在内在联系
  2. 能量尺度分离:低温峰对应自旋相关能量尺度,高温峰对应电荷激发
  3. 阻挫-关联竞争:t'/t有效调节几何阻挫和电子关联之间的竞争

实验意义

研究结果可能为解释Dy₃Sb₃Zn₂O₁₄中观察到的低温比热异常提供理论框架。位点混合无序将理想kagome晶格转化为类似稀释二维三角几何的中间结构。

局限性

  1. 温度限制:由于符号问题,模拟尚未达到实验可及的极低温度
  2. 相互作用强度:限制在U ≤ 3.5t的中等耦合区域
  3. 系统尺寸:主要结果基于相对较小的系统尺寸

未来方向

  1. 发展新的数值方法以克服符号问题,研究更强耦合区域
  2. 扩展到更大系统尺寸和更低温度
  3. 研究其他阻挫几何中的类似现象

深度评价

优点

  1. 方法创新:引入可调参数t'/t实现连续插值的思路新颖且物理图像清晰
  2. 系统研究:全面分析了磁关联、热力学性质和符号问题
  3. 物理洞察:揭示了磁性和热力学响应之间的深层联系
  4. 实验相关性:为理解实际材料中的异常现象提供了理论支持

不足

  1. 模型简化:实际材料中的复杂性(如无序、长程相互作用)未完全考虑
  2. 温度限制:符号问题限制了对极低温行为的探索
  3. 机制解释:对磁-热关联的微观机制需要更深入的理论分析

影响力

该研究为理解阻挫量子磁体中的关联效应提供了重要见解,对于设计和理解新型量子材料具有指导意义。所发展的方法和发现的现象可能在其他阻挫系统中具有普遍性。

适用场景

  1. 阻挫磁性材料的理论建模
  2. 量子蒙特卡罗方法的应用指导
  3. 实验数据的理论解释框架

参考文献

论文引用了43篇重要文献,涵盖了kagome晶格的理论研究、实验进展、数值方法等多个方面,为研究提供了坚实的基础。