2025-11-13T11:31:10.165583

Localized fermions on the triangular lattice with Ising-like interactions

Regeciová, Kapcia
The model of localized fermions on the triangular lattice is analyzed in means of the Monte Carlo simulations in the grand canonical ensemble. The Hamiltonian of the system has a form of the extended Hubbard model (at the atomic limit) with nearest-neighbor Ising-like magnetic $J$ interactions and onsite Coulomb $U$ interactions. The model is investigated for both signs of $J$, arbitrary $U$ interaction and arbitrary chemical potential $μ$ (or, equivalently, arbitrary particle concentration $n$). Based on the specific heat capacity and sublattice magnetization analyses, the phase diagrams of the model are determined. For ferromagnetic case ($J<0$), the transition from the ordered phase (which is a standard ferromagnet and can be stable up to $k_{B}T/|J| \approx 0.61$) is found to be second-order (for sufficiently large temperatures $k_{B}T/|J| \gtrsim 0.2$) or first-order (for $-1<U/|J|<-0.65$ at the half-filling, i.e., $n=1$). In the case of $J>0$, the ordered phase occurs in a range of $-1/2<U/|J|<0$ (for $n=1$), while for larger $U$ the state with short-range order is also found (also for $n \neq 1$). The ordered phase is characterized by an antiferromagnetic arrangement of magnetic moments in two sublattices forming the hexagonal lattice. The transition from this ordered phase, which is found also for $μ\neq 0$ ($n \neq 1$) and $U/|J|>-1/2$ is always second-order for any model parameters. The ordered phase for $J>0$ can be stable up to $k_{B}T/|J| \approx 0.06$.
academic

Localized fermions on the triangular lattice with Ising-like interactions

基本信息

  • 论文ID: 2412.10808
  • 标题: Localized fermions on the triangular lattice with Ising-like interactions
  • 作者: Lubomíra Regeciová (Slovak Academy of Sciences), Konrad Jerzy Kapcia (Adam Mickiewicz University)
  • 分类: cond-mat.stat-mech cond-mat.other cond-mat.quant-gas cond-mat.str-el physics.comp-ph
  • 发表时间: 2025年10月17日
  • 论文链接: https://arxiv.org/abs/2412.10808

摘要

本研究通过蒙特卡罗模拟在正则系综中分析了三角晶格上局域化费米子模型。系统的哈密顿量具有扩展Hubbard模型（原子极限）的形式，包含最近邻伊辛型磁相互作用J和在位库仑相互作用U。研究考察了J的两种符号、任意U相互作用和任意化学势μ（或等价地，任意粒子浓度n）。基于比热容和子晶格磁化强度分析，确定了模型的相图。对于铁磁情况（J<0），从有序相的转变可以是二级相变（在足够高温度kBT/|J|≳0.2时）或一级相变（在半填充时-1<U/|J|<-0.65）。对于J>0的情况，有序相出现在-1/2<U/|J|<0范围内（n=1时），而对于更大的U值，还发现了短程有序态。

研究背景与动机

  1. 核心问题: 研究几何受挫系统中局域化费米子的磁性相变行为，特别是三角晶格上的扩展Hubbard模型在原子极限下的相图结构。
  2. 重要性:
    • 三角晶格具有固有的几何受挫特性，导致反铁磁相互作用无法同时满足所有相互作用
    • 该模型是经典伊辛模型向费米子系统的非平凡推广
    • 对理解强关联电子系统和磁性材料具有重要意义
  3. 现有方法局限性:
    • 超立方晶格的研究较为充分，但三角晶格的精确解析结果缺乏
    • 平均场近似高估了相变温度
    • 缺乏对几何受挫效应的详细数值研究
  4. 研究动机:
    • 填补三角晶格上该模型的数值研究空白
    • 探索几何受挫对相变行为的影响
    • 为实验材料（如Li2MnTeO6、KCeO2等）提供理论基础

核心贡献

  1. 首次系统性蒙特卡罗研究: 对三角晶格上带伊辛型磁相互作用的扩展Hubbard模型进行了全面的数值模拟
  2. 完整相图构建: 确定了铁磁（J<0）和反铁磁（J>0）情况下的温度-化学势-相互作用强度三维相图
  3. 相变性质分析: 精确确定了一级和二级相变的边界，发现反铁磁情况下所有相变均为二级
  4. 受挫效应揭示: 发现几何受挫显著影响相图结构，反铁磁情况下有序相稳定温度远低于铁磁情况
  5. 短程有序态识别: 在反铁磁情况下发现并表征了短程有序态的存在区域

方法详解

任务定义

研究三角晶格上局域化费米子系统的相变行为，输入为模型参数（J, U, μ, T），输出为系统的相图和相变性质。

模型架构

哈密顿量: H^=U∑in^i↑n^i↓+2Jz∑⟨i,j⟩s^is^j−(μ+U2)∑in^i\hat{H} = U\sum_i \hat{n}_{i\uparrow}\hat{n}_{i\downarrow} + \frac{2J}{z}\sum_{\langle i,j\rangle} \hat{s}_i\hat{s}_j - \left(\mu + \frac{U}{2}\right)\sum_i \hat{n}_i

其中：

  • s^i=(n^i↑−n^i↓)/2\hat{s}_i = (\hat{n}_{i\uparrow} - \hat{n}_{i\downarrow})/2：局域自旋算符
  • n^i=n^i↑+n^i↓\hat{n}_i = \hat{n}_{i\uparrow} + \hat{n}_{i\downarrow}：粒子数算符
  • z=6z = 6：三角晶格配位数

格点状态: 每个格点可处于四种状态：

  • ∣0⟩|0\rangle：空位
  • ∣↑⟩|\uparrow\rangle：单占据（自旋向上）
  • ∣↓⟩|\downarrow\rangle：单占据（自旋向下）
  • ∣↑↓⟩|\uparrow\downarrow\rangle：双占据

技术创新点

  1. 正则系综模拟: 采用粒子数涨落的正则系综，更准确描述化学势控制的系统
  2. 有限尺寸分析: 系统研究了L = 12×12到36×36的有限尺寸效应
  3. 相变判识方法:
    • 使用Challa等人的能量分布函数P(E)方法区分一级和二级相变
    • 结合比热容和子晶格磁化强度的标度行为分析
  4. 子晶格分解: 将三角晶格分解为三个等价子晶格，精确分析磁有序结构

实验设置

数据集

  • 系统尺寸: 三角形团簇，L = 12×12, 24×24, 36×36格点
  • 参数范围:
    • 温度：kBT/|J| = 0.005 - 1.0
    • 相互作用：U/|J| = -1.0 - 5.0
    • 化学势：μ/|J| = -0.6 - 0（利用粒子-空穴对称性）

评价指标

  1. 比热容: C=⟨ε2⟩−⟨ε⟩2/T2C = \langle\varepsilon^2\rangle - \langle\varepsilon\rangle^2/T^2
  2. 子晶格磁化强度: mα=(6/L)∑i∈α⟨s^i⟩m_\alpha = (6/L)\sum_{i\in\alpha}\langle\hat{s}_i\rangle
  3. 粒子浓度: n=⟨n^i⟩n = \langle\hat{n}_i\rangle
  4. 能量分布函数: P(E)P(E)用于相变性质判断

实现细节

  • 算法: Metropolis蒙特卡罗算法
  • 平衡化步数: 10610^6步
  • 统计步数: 10610^6步
  • 温度扫描: 从高温向低温逐步冷却
  • 边界条件: 周期性边界条件

实验结果

主要结果

铁磁情况 (J < 0):

  • 最高相变温度：kBT/|J| ≈ 0.605（U → ∞时，接近经典伊辛模型的0.607）
  • 三临界点位置：kBT/|J| ≈ 0.18-0.20，与U和μ几乎无关
  • 相变性质：低温时为一级相变，高温时为二级相变
  • 有序相：标准铁磁相，所有子晶格磁化强度相等

反铁磁情况 (J > 0):

  • 最高相变温度：kBT/|J| ≈ 0.06（远低于铁磁情况）
  • 有序相存在范围：-1/2 < U/|J| < 0（半填充时）
  • 相变性质：所有相变均为二级
  • 有序相结构：两个子晶格形成反铁磁序，第三个子晶格磁化强度为零

消融实验

通过改变系统尺寸验证了相变的热力学极限行为：

  • 一级相变：能量分布出现双峰结构
  • 二级相变：比热容出现λ型奇异性，峰值随系统尺寸增长
  • 短程有序：比热容出现不随系统尺寸标度的宽峰

案例分析

典型相变行为:

  1. U = -0.7, μ = -0.10：一级F-NO相变，伴随磁化强度和粒子数的不连续跳跃
  2. U = -0.6, μ = -0.12：二级F-NO相变，连续的序参量变化
  3. U = -0.2, μ = -0.15：二级AF-NO相变，两个子晶格的反铁磁有序

实验发现

  1. 几何受挫的强烈影响: 反铁磁情况下相变温度比铁磁情况低一个数量级
  2. 相变性质的不对称性: 铁磁情况既有一级又有二级相变，反铁磁情况仅有二级相变
  3. 短程有序的普遍存在: 在反铁磁情况下广泛观察到短程有序态

相关工作

  1. 经典伊辛模型: 三角晶格上的精确解（Houtappel, 1950）为本研究提供了基准
  2. 扩展Hubbard模型: 超立方晶格上的变分和数值研究为方法学提供了基础
  3. 受挫磁性系统: 三角晶格反铁磁体的理论和实验研究为结果解释提供了背景
  4. 冷原子实验: 光晶格中超冷原子气体的实现为模型的实验验证提供了可能

结论与讨论

主要结论

  1. 三角晶格的几何受挫显著影响局域化费米子系统的相变行为
  2. 铁磁和反铁磁相互作用导致完全不同的相图结构和相变性质
  3. 反铁磁情况下的有序相具有独特的两子晶格反铁磁结构
  4. 所有反铁磁相变均为连续相变，而铁磁情况下存在一级-二级相变的转换

局限性

  1. 系统尺寸限制: 受计算资源限制，最大系统仅为36×36
  2. 原子极限近似: 忽略了粒子跳跃项，可能影响实际材料的定量预测
  3. 经典自旋近似: 在强量子涨落区域可能不够准确

未来方向

  1. 包含更长程相互作用的研究
  2. 量子蒙特卡罗方法的应用
  3. 与实际材料的定量比较
  4. 冷原子系统中的实验验证

深度评价

优点

  1. 方法学严谨: 采用了成熟的蒙特卡罗方法和相变判识技术
  2. 结果全面: 构建了完整的三维相图，涵盖了多个参数区域
  3. 物理洞察深刻: 揭示了几何受挫对相变行为的本质影响
  4. 实验相关性强: 与多种实际材料和冷原子系统相关

不足

  1. 理论分析不足: 缺乏对相变机制的深入理论解释
  2. 动力学性质缺失: 未涉及系统的动力学和输运性质
  3. 量子效应忽略: 原子极限下忽略了重要的量子涨落效应

影响力

  1. 学术价值: 为受挫磁性系统理论提供了重要的数值基准
  2. 实用价值: 为相关材料的理解和设计提供了指导
  3. 可复现性: 方法描述详细，结果具有良好的可重现性

适用场景

  1. 三角晶格磁性材料的理论建模
  2. 冷原子光晶格实验的理论指导
  3. 强关联电子系统的相变研究
  4. 受挫量子多体系统的基准测试

参考文献

本论文引用了71篇重要文献，涵盖了伊辛模型理论、Hubbard模型研究、蒙特卡罗方法、受挫磁性系统以及相关实验材料等多个方面，为研究提供了坚实的理论基础和实验背景。