Competing ferromagnetic and antiferromagnetic phases on the frustrated Ising honeycomb lattice

本研究利用簇平均场方法，研究了蜂窝晶格上受挫的 $J_1$-$J_2$-$J_3$ 伊辛模型，旨在揭示增加第二近邻铁磁耦合如何使系统的三临界点向强受挫极限移动，并最终导致铁磁相、反铁磁相与顺磁相共存的双临界点。

要点

想象一个巨大的、六边形蜂窝状晶格，就像一个巨大的蜂巢，每个交点都持有一个微小的磁铁（“自旋”），它们可以指向“上”或“下”。这就是这篇论文所讲述故事的舞台。

科学家们正在玩一场关于这些磁铁的拔河游戏，受三组不同的规则（相互作用）支配，这些规则在不同的方向上拉扯它们：

1. 最好的朋友 ($J_1$)： 这些是紧挨着的磁铁。它们是铁磁性的，这意味着它们非常希望手牵手，指向相同的方向（全部向上或全部向下）。
2. 表亲 ($J_2$)： 这些是隔开一步远的磁铁。在本研究中，它们也是铁磁性的，所以它们也希望每个人都能达成共识，指向相同的方向。
3. 对手 ($J_3$)： 这些是隔开三步远的磁铁。它们是反铁磁性的，意味着它们是固执的对手。它们希望自己的邻居指向相反的方向。

伟大的冲突

本文关注一个特定的、棘手的场景：当“最好的朋友”和“表亲”很强，想要大家保持统一（铁磁性），而“对手”也很强，想要创造一种黑白相间的棋盘格图案（反铁磁性）时。

当这两种欲望发生冲突时，系统会产生挫折感（frustrated）。这就像一群人试图决定去哪家餐厅吃饭，一半人想吃披萨，另一半人想吃寿司，但喜欢披萨的人同时也和喜欢寿司的人是好朋友。这群人无法简单地做出选择；他们必须找到一个折中方案，否则决策过程就会变得一团糟。

科学家的工具箱：“团簇”法

为了弄清楚当加热这个系统（增加热能）时会发生什么，研究人员使用了名为**团簇平均场（Cluster Mean-Field）**的方法。

想象一下尝试预测整个体育场的气氛。

• 旧方法（单点法）： 你观察一个人，并根据他来猜测整个观众席。这通常过于简单，容易忽略混乱。
• 本文的方法（团簇法）： 他们一次观察一小组（6人或18人）的人。他们精确计算这些小团体如何相互作用，然后利用平均值来推测其余的部分。这能更清晰地描绘出人群中发生的“挫折感”。

他们的发现

1. “由无序产生有序”的惊喜
通常，我们认为热量（无序）会破坏有序。如果你加热一块磁铁，它就会停止指向一个方向，变得随机。
然而，在“对手”力量最强的点附近，科学家们发现了一个奇特的现象，叫做**“由无序产生有序”（Order-by-Disorder）**。

• 类比： 想象一个房间里的人，大家对于站成一个圆圈还是一个正方形都感到同样满意。这是一个平局。但是，如果你开始摇晃房间（增加热量），站成正方形的人可能会发现，比起圆圈，他们在晃动时更容易活动而不至于互相碰撞。突然间，“正方形”变成了首选位置，这并不是因为它在静止时最舒适，而是因为它在混乱发生时最具灵活性。
• 结果： 在这个模型中，增加热量实际上稳定了铁磁相（全部方向相同）而非反铁磁相。热量帮助系统“选择了”一个赢家。

2. 形态变化的转变
科学家们绘制了随着改变“对手”磁铁的强度 ($J_3$) 和温度，系统如何从有序（磁铁排列整齐）变为无序（随机）的图谱。

• 平滑转变（二阶相变）： 有时，这种变化就像冰融化一样。它发生得非常渐进。随着温度升高，磁铁逐渐失去它们的排列。
• 突然跳跃（一阶相变）： 有时，这种变化就像水沸腾一样。它是一个突然的、剧烈的切换。系统会从一个状态猛然跳到另一个状态。
• 三临界点（Tricritical Point）： 这是“金发姑娘”式的完美平衡点，即转变从平滑变为突然的转折点。这是游戏规则发生改变的临界点。

3. “表亲” ($J_2$) 的影响
研究人员发现，如果让“表亲”磁铁 ($J_2$) 变得更强，系统就会变得不那么混乱。

• 类比： 如果“对手” ($J_3$) 非常强，系统就会处于一场混乱的拔河中，导致那些突然的、剧烈的跳跃（一阶相变）。但如果加强“表亲” ($J_2$)，他们就能帮助“最好的朋友” ($J_1$) 守住阵线。
• 结果： 随着“表亲”变得更强，那些混乱的、突然的跳跃消失了。系统开始再次进行平滑的变化。最终，如果“表亲”足够强，系统将只进行平滑变化，而“三临界点”会消失，取而代之的是一个双临界点（Bicritical Point），在那里三种不同的相和平静地汇聚在一起。

核心结论

这篇论文是一张复杂的磁性景观地图。它展示了当存在竞争力量（朋友想要统一 vs 对手想要对立）时，系统会表现出各种奇特行为：

• 它可能会从一个状态突然跳跃到另一个状态。
• 它可以用热量来挑选赢家（由无序产生有序）。
• 它拥有特殊的“汇合点”（三临界点和双临界点），在这些点上物理规则会发生转变。

这项研究证实，虽然“最好的朋友” ($J_1$) 和“表亲” ($J_2$) 总是会导致向随机性的平滑转变，但“对手” ($J_3$) 的引入创造了一个丰富且复杂的领域，充满了突然的跳跃和特殊的临界点，尤其是在“表亲”的力量还不足以完全压制冲突的时候。

技术摘要

技术摘要：受挫伊辛蜂窝晶格上竞争的铁磁与反铁磁相

问题陈述
本研究调查了受挫 $J_1$–$J_2$–$J_3$ 伊辛模型在蜂窝晶格上的相变。该系统由铁磁性第一近邻（$J_1 > 0$）和第二近邻（$J_2 > 0$）耦合，以及反铁磁第三近邻相互作用（$J_3 < 0$）定义。尽管之前的理论研究广泛探讨了 $J_2$ 较弱或为负，以及 $J_3$ 为铁磁性的情形，但在 $g_2 \equiv J_2/J_1 > 1/2$ 的参数范围在很大程度上被忽视了。在此机制下，$J_1$ 和 $J_2$ 的铁磁性质与由反铁磁 $J_3$ 引入的受挫性相互竞争，创造了一个复杂的景观，其中铁磁（FE）相和反铁磁（AF）相在稳定性方面展开竞争。作者旨在表征热相变，特别是序-无序边界的性质以及是否存在多重临界点。

方法论
作者采用簇平均场（CMF）近似来分析该模型。这种方法通过精确处理簇内相互作用来纳入短程相关性，同时通过平均场解耦来近似簇间相互作用，从而改进了标准的单点平均场理论。

• 簇大小： 计算使用两种不同的簇几何结构进行：6位簇（$n_s=6$）和18位簇（$n_s=18$）。6位簇允许通过自由能导数的解析处理来识别临界点，而18位簇则通过包含更多拓扑非等效位点，为结果的稳健性提供更严格的检查。
• 分析工具： 通过分析以下内容来确定相变的性质：
  • 相之间的自由能差（$\Delta f$）。
  • 自由能对序参数（$\eta_2$ 和 $\eta_4$）的导数，以区分二阶相变、一阶相变和三临界点。
  • 序参数（铁磁性的磁化强度和反铁磁性的交错磁化强度）以及每个自旋的熵。
  • 通过基态能量比较来构建零温相图。

主要贡献与结果
研究绘制了在不同固定 $g_2$ 值（$0.6, 0.75, 1.0$）下，温度（$T$）相对于第三近邻耦合（$g_3$）平面的相图。主要发现如下：

1. 相变性质：

   • FE–PM 转变： 铁磁相与顺磁（PM）相之间的转变在所研究的整个相互作用范围内始终被发现是二阶的。这一结果与近期的蒙特卡洛模拟及先前的 CMF 研究一致，表明当 $J_1$ 为铁磁性时，该连续特征是该模型的稳健特征。
   • AF–PM 转变： 反铁磁相与顺磁相之间的边界表现出丰富的行为。对于 $g_2 = 0.6$ 和 $g_2 = 0.75$，AF–PM 转变可以是一阶或二阶的，并由一个三临界点分隔。
     • 在6位近似中，一阶转变区域存在于 $-1.23 < g_3 < -1.04$ 之间（在 $g_2=0.6$ 时）。
     • 18位簇证实了存在一阶转变，但预测发生一阶转变的 $g_3$ 范围比6位簇更宽，这表明该现象并非近似规模带来的伪影。

2. 增加 $J_2$ ($g_2$) 的影响：

   • 增加铁磁性第二近邻耦合（$g_2$）会系统地缩小发生一阶 AF–PM 转变的 $g_3$ 范围。
   • 在足够高的 $g_2$ 下（具体为 $g_2 = 1.0$），一阶 AF–PM 转变被完全抑制，仅留下二阶转变。因此，三临界点向 $g_3 \approx -1$ 移动，并最终转化为一个双临界点，在该点处 FE、AF 和 PM 相通过二阶线相遇。

3. “通过无序实现有序”（Order-by-Disorder）与连续相变：

   • 在强受挫极限（$g_3 \approx -1$）附近，系统表现出通过无序实现有序的状态选择。热涨落由于熵的选择（具有更高的激发态密度）使铁磁相优于反铁磁相。
   • 在强竞争机制下（靠近 $g_3 = -1$），系统随着温度升高经历两次连续相变：首先是从 AF 到 FE 的一阶转变，随后是从 FE 到 PM 的二阶转变。

4. 临界终点：

   • 二阶 FE–PM 边界终止于一个临界终点，该点与一阶 AF–PM 转变线相遇（对于较低的 $g_2$）。

意义与主张
本文声称提供了对受挫蜂窝伊辛模型在先前未被充分探索的参数范围（$g_2 > 0.5$）内更深入的理解。作者断言，通过两个簇规模验证的 CMF 结果表明：

• 纳入反铁磁第三近邻耦合不仅对于稳定 AF 相至终关重要，而且对于驱动 AF–PM 转变中三临界性的出现也至关重要，而这一特征在仅具有铁磁耦合的简单 $J_1$–$J_2$ 模型中是不存在的。
• 铁磁与反铁磁相互作用之间的竞争导致了复杂的现象，包括通过无序实现有序机制以及连续相变，这些现象与描述具有强伊辛各向异性和竞争交换相互作用的真实范德华磁性材料（如 FePS$_3$, VBr$_3$ 和 CrI$_3$）相关。
• 二阶 FE–PM 转变的稳健性以及增加 $J_2$ 对一阶 AF–PM 转变的抑制作用，被确定为该受挫系统的关键特征。

作者总结道，虽然其 CMF 方法捕捉到了本质的受挫效应，但仍需使用补充技术（蒙特卡洛模拟、张量网络等）进行进一步研究，以充分解析这些系统的临界行为和潜在的场诱导现象。