任意经典烧绿石自旋冰中的拓扑相变及其热力学命运

《Physical Review X》：Topological Phase Transitions and Their Thermodynamic Fate in Arbitrary- Classical Pyrochlore Spin Ice

　　我们建立了一个自包含的理论框架，对具有任意自旋 ?? 的经典烧绿石磁体的拓扑相和临界现象进行分类，其中自旋受竞争性交换作用和单离子各向异性支配。我们首先在无磁单极子（冰规则）极限下工作。在小 ?? 区域，单离子项倾向于低自旋振幅，对偶性揭示了一个二分性：整数自旋表现出连续的 3D ???? 退禁闭转变，而半整数自旋则保持在 U(1) 库仑液体中，没有任何转变。在大 ?? 区域，局域自旋振幅达到最大值（|????| = ??），宏观通量被量子化为 2?? 的整数倍。通过将缺陷结构映射到拓扑环气，我们证明了物理冰规则与模 ?2?? 通量守恒之间的相容性在顶点层面成立当且仅当 ?? ≤ 3/2。对于 ?? = 3/2，这将该系统映射到三态波茨模型，其对称性允许的三次不变项驱动了一级相变。对于 ?? ≥ 2，不存在到离散钟模型的映射。利用配分函数的分解，我们表明对离散 ?2?? 结构敏感的唯二过程携带指数小的统计权重，并对应于 3D ???? 不动点处的无关算符，因此退禁闭转变属于 3D ???? 类型。最后，引入热磁单极子后，我们表明它们起到对称性破缺的有效磁场的作用，切断缺陷弦。因此，连续相变被抹圆为交叉行为，而 ?? = 3/2 的一级相变被预测在有限温度下存续，并终止于临界端点。对高达 ?? = 7/2 的经典蒙特卡洛模拟验证了在可达到温度下的预测交叉行为；这些温度高于估算的临界端点，因此 ?? = 3/2 一级相变的预测存续性仍有待直接检验。

I. 引言

对烧绿石晶格上几何阻挫磁性的研究一直是发现奇异物质态的沃土，最著名的是经典和量子自旋液体 [1–4]。典型例子是经典自旋冰，在稀土氧化物如 和 中实现 [5–8]。在这些材料中，作用于大总角动量（ 对 ， 对 ）的强晶体电场隔离出一个低能二重态，使得磁矩可以被视为沿局域 四面体轴取向的有效 伊辛赝自旋 [9]。阻挫交换相互作用迫使自旋进入二进二出的冰规则 [10]，将低能物理映射为一个涨落的、无散度的有效磁场。这种涌现电磁学表征了 U(1) 库仑液相，其特点是代数自旋关联和针点奇异性 [11–15]。

最近的理论和数值进展已超越刚性的 范式，探索了 及更高自旋系统，其中占主导的易轴交换与单离子各向异性竞争 [16–19]。在整篇论文中， 表示位点 处沿局域 轴（连接共享该位点的两个四面体中心的轴）的自旋分量，交换作用的横向分量设为零，模型完全在经典层面处理。对于 情形，占主导的易轴交换与惩罚大 的单离子各向异性 ( ) 之间的竞争产生了丰富的相图 [16–18]。与 模型不同， 自旋可以进入非磁性 态。在无磁单极子（冰规则）极限下，调节各向异性驱动系统经历三个不同的宏观相：平庸顺磁体、U(1) 库仑相和一个 -受限库仑相（一种库仑液体，其宏观极化通量限制为偶数整数，因此携带单个通量单位的弦被禁闭，而成对的弦保持自由），它们由 3D 和 3D 伊辛普适类中的连续相变分隔。最近，Pandey 等人研究了 情形，发现了一个分隔两种不同库仑液体的一级 通量禁闭相变 [19]。这些发展联系到对烧绿石系统中分数化和磁碎裂的更广泛理解 [20–22]。同一类模型也出现在更广泛的系统中，包括高压下的水冰 [23] 和人工纳米磁体阵列 [24]。

这些发现提出了基本问题。宏观拓扑结构和临界现象如何推广到任意自旋 ？为什么 模型表现出连续相变而 表现出连续相变，并且能否从微观映射而不是对称性考虑推断出一级特征 [19]？对于 ，相变和普适类的最终命运是什么？最后，在有限温度下，当严格冰规则不可避免地被热激发的磁单极子违反时，这些拓扑相变的热力学命运是什么？

在本文中，我们通过为自旋 烧绿石磁体建立一个统一的理论框架来回答这些问题，首先处理无磁单极子极限，然后处理热磁单极子。在无磁单极子极限下，我们首先证明相图仅包含上述已识别的相：中间自旋振幅的均匀背景被单离子能量的凸性排除，而“部分退禁闭”情景（其中更高拓扑电荷的弦会在基本弦之前凝聚）通过竞争缺陷弦的自由能比较被排除（第三节）。

然后我们确定每个相变的性质。在小 区域，配分函数在双对相位场中的精确重写表明，所有整数 模型在共同的、与 无关的耦合处经历 3D 相变，而所有半整数 模型约化为普通的 自旋冰，因此根本没有小 相变（第四节）。在大 区域，我们将缺陷弦映射到环气（第五节）。对于 ，环气与三态波茨模型的高温展开一致，这从微观上解释了参考文献 [19] 中数值观察到的一级相变，并定量地定位了它。对于 ，我们表明没有离散钟模型能重现环气；相反，环气是连续 3D 模型的环气，加上统计权重指数小的修正，且这些修正对应的算符在 不动点处无关，因此相变是连续的，属于 3D 类型。由此产生的无磁单极子极限下的分类总结在表 I 和表 II 以及图 1 的示意相图中。

表 I. 无磁单极子极限下整数自旋模型的相图 ( )。宏观极化通量 作为区分三个相的拓扑不变量。从平庸顺磁体出来的相变对于所有整数 都属于 3D 普适类。从 -受限相出来的退禁闭相变对于 是 3D 伊辛，对于 是 3D ，后者由第五节 C 的分解和 -聚变级联的无关性支持。 临界逸度为 ，其中 是金刚石晶格上的 3D 伊辛临界耦合 [18]。对于 ， 且 （见第五节 C）。在有限温度下 ( )，热磁单极子将所有这些相变抹圆为连续交叉行为（见第六节）。值 是第四节的前导解析估计；对于 的测量值为 0.3609(1) [16]。

自旋 平庸相 (相变) U(1) 库仑相 (退禁闭) -受限相
1 (3D ) (3D 伊辛)
(3D ) (3D )

表 II. 无磁单极子极限下半整数自旋模型的相图 ( )。非磁性 态的缺失禁止了平庸顺磁相，将相图减少为两个相（除了 ，它只有一个相）。对于 ，单离子各向异性是常数，系统对所有 都处于 U(1) 库仑相。 退禁闭相变是一级的，由 缺陷弦聚变驱动，映射到三态波茨模型；引用的 是闭式 Bethe-Peierls 估计 ，它高估了参考文献 [19] 的蒙特卡洛值 (即 ) 约 。对于 ，聚变级联的指数抑制恢复了 3D 临界性。在有限温度下 ( )，连续相变变成交叉行为，而 一级相变被预测持续到临界端点（见第六节）。

自旋 U(1) 库仑相 (退禁闭) -受限相
None [始终为 U(1) 相]
(一级；三态波茨) (3D )

图 1. 平面内的示意相图，分别针对 (a) 整数 ，(b) ，和 (c) 半整数 ；每个面板中标记了各个区域， 和 以 为单位表示。虚线表示交叉脊线（在 处无热力学奇异性）；(b) 中的实线表示一级共存线，终止于临界端点（空心圆）。脊线和共存线从原点沿由无磁单极子临界逸度 确定的射线 （表 I 和 II）发出；在每条射线上 是常数，锐利相变仅在无磁单极子极限 下以固定 的方式实现（第二节 C）。对于整数 (a)，U(1) 库仑区域是有限张角楔形，由通过 和 的射线界定。在 平面中的等价表示，其中无磁单极子相变表现为 轴上的独立点，在补充材料的第 I 节中给出 [25]。

最后，通过将热磁单极子纳入配分函数，我们表明它们以对称性破缺场的形式进入对偶描述，允许缺陷弦终止。这将在任何 下将所有连续相变抹圆为交叉行为，并预测 一级线是唯一能在有限温度下存续的边界，终止于临界端点（第六节）。我们用高达 的经典蒙特卡洛模拟验证了交叉行为（第七节）；这些模拟可达到的温度高于估算的临界端点，因此 一级线本身的存续性仍是未来数值研究的开放目标。超越自旋冰的具体背景，我们的分析以定量方式说明了一个微观守恒律为实值（因此为 U(1)）的晶格模型如何在低能下模仿离散规范结构，以及这种模仿何时必须失败。

为清晰起见，我们明确说明本文所依赖的每项近期工作建立了什么，以及此处的新内容是什么。参考文献 [16] 通过直接在冰规则流形中进行的虫算法模拟绘制了 相图，定位了两次相变 [ 和 ]，并通过映射到 配分函数结合有限尺寸标度确定了它们的普适类；参考文献 [18]（由本作者）通过双对 和环气表示解析地推导了相同的 结果，并量化了热磁单极子在这些对偶表示中作为显式对称性破缺场的影响；参考文献 [17] 在强场下研究了一个相关的二维模型；参考文献 [19] 通过相同的虫算法方法发现了 的一级 相变，并通过基于对称性的论证进行了合理化。

本文在单一框架内包含所有 ，我们希望同样明确其新颖性所在和不在何处。变换本身并不是新的：下文使用的对偶性、环气展开和平均场评估是晶格统计力学的标准工具，我们在出现处注明了其来源。新的是一旦为任意 建立这些工具后它们所揭示的结构。

在小 区域，任意 的双对链权重揭示了一个在 处不可见的宇称二分性：前导 耦合，即相变位置，对每个整数 都相同，而对于半整数 ，双对权重是反周期的——半整数自旋的对偶标志——模型约化为普通自旋冰，根本没有相变。在大 区域，组织性结果是第五节 A 的冰规则相容性判据，这是一个基本的几何陈述，据我们所知，文献中没有对应物：离散 描述在顶点层面与物理冰规则一致仅对 成立，这用单一论点解释了为什么伊辛和波茨映射在 和 处成功，以及为什么每个这样的映射对 都必须失败。

在这一分界线的 侧，该判据将参考文献 [19] 的基于对称性的理由升级为微观推导——环气被识别为三态波茨模型的高温展开，因此一级特征先验地由对称性允许的三次不变项得出，而不是从数值推断——然后 Bethe-Peierls 评估以闭式定位相变，与数值值相差 3% 以内。在另一方面，对于此前不存在任何数值或解析结果的情形，我们证明没有任何离散时钟模型能够复现环气模型，通过配分函数的分解约束残余的离散修正，并由此获得整个族的三维普适类和闭合形式的临界线。第三部分的完备性论证和第七部分的蒙特卡洛模拟同样是首次在此报告。第六部分在所述一般性下的有限温度分析也是如此：尽管自旋冰中的热磁单极子传统上通过连续静电类比来处理[8,14,29]，我们却在对偶变量中精确表示磁单极子气体，其中电荷为的磁单极子是局域算符，因此该气体充当显式的U(1)对称性破缺场——这一对应关系在文献[18]中针对引入，并在此处对每个建立，最终归结为时零场自生成的临界终点，据我们所知，这是对受挫磁体的首个此类预测（偶极自旋冰的已知临界终点，即场驱动的磁单极子液-气转变[8]，需要有限磁场）。

本文组织如下。第二部分定义模型并推导精确配分函数。第三部分建立宏观磁通量子化和相图的完备性。小对偶映射及整数与半整数二分性在第四部分给出，而大环气映射和退禁闭转变的分类在第五部分发展。热磁单极子的有限温度效应在第六部分分析。第七部分给出有限温度蒙特卡洛模拟，第八部分给出实验展望，第九部分以未来方向作结。

二、模型表述与精确配分函数

在着手任意自旋的推广理论之前，我们严格定义微观哈密顿量并引入配分函数的精确表示，该表示构成后续所有对偶映射的基础。

A. 晶格几何与有效哈密顿量

烧绿石晶格是由共角四面体组成的网络，其格点位于底层二分金刚石晶格的键中点处[7,30]。金刚石晶格包含个顶点（等价地，烧绿石晶格的个四面体）和条键。我们在每个格点上放置离散自旋变量，其大小为任意正整数或半整数。每条键上的局域量子化轴从子晶格顶点指向子晶格顶点。

我们研究该晶格上的推广Blume-Capel哈密顿量[31,32]，它将反铁磁最近邻交换与单离子各向异性耦合起来：

(1)

其中外层求和遍及每个金刚石晶格顶点，标记其四条关联键（配位数）。由于两项均关于二次，哈密顿量在同时符号翻转（全局时间反演对称性）下不变。此外，由于每个顶点处的交换能是平方散度，其值不依赖于键取向所采用的子晶格符号约定。考虑每条键的子晶格取向，离散晶格散度定义为和（子晶格上的符号翻转源于每条键上的局域量子化轴从指向），这可解释为磁单极子电荷[33]。

B. 精确配分函数

在温度下，正则配分函数为

(2)

我们用两个无量纲逸度参数化玻尔兹曼权重：控制单离子代价，控制磁单极子代价。物理上，（惩罚大自旋振幅（易面倾向），而（有利于最大振幅（易轴倾向）。在每个顶点引入独立的磁单极子电荷变量，并通过克罗内克δ函数施加约束，我们将配分函数重写为

(3)

该式干净地将拓扑权重与局域各向异性权重分离，对任意和均成立。此表示是下文所有解析结果的出发点。

C. 无磁单极子极限与宏观极化磁通

当时，磁单极子逸度指数小，因此任何含的构型都承受禁止性的玻尔兹曼惩罚。在式(3)中仅保留扇区即施加每个顶点处的冰规则。这个无磁单极子极限（）是本文大部分内容的设定：第三至第五部分专门在处工作，而的影响在第六部分分析并在第七部分数值探测。我们要强调，固定时的是对两个耦合的联合极限，而非固定耦合下的零温极限：由于，固定且意味着固定，因此随一起趋于零。变量因此描述了物理上自然的层级，其中约为量级，冰规则流形内的残余热涨落由而非本身控制。等价地，对于具有固定和的材料，冷却使单调地从1偏离，因此由组织的各个区域在冷却到时被遍历——仍远低于，其中指数小。这正是文献[16,19]获得相图时所采用的极限。我们强调，第三至第五部分分析的尖锐转变严格存在于流形上：在任何处磁单极子逸度非零，并且如第六部分所示，连续转变随后被圆化为交叉行为，唯一的例外是处的第一性线。

冰规则流形进一步分解为以宏观极化磁通标记的拓扑不同扇区。在周期边界条件下，沿立方轴的分量为

(4)

其中表示与任何恒定平面相交且贯穿晶格的键集合。符号因子将每个局域自旋投影到全局方向上。由于无散条件迫使所有进入一个薄片的磁通必须离开，因此独立于的选择，从而构成守恒的拓扑量子数[20,34]。

D. 基础情形回顾（和）

对于基础情形，单个自旋只能取值。单离子各向异性项为常数：。由于这均匀移动全局能量，参数在物理上变得无关紧要。热力学仅由拓扑冰规则决定。宏观磁通可自由取任何整数值，（因为任何宏观平面与偶数条键相交，求和偶数个半整数自旋得到整数）。宏观态是典型的U(1)库仑液相[11,14,35]。

极限引入非磁性态。对于（），非磁性态在能量上有利，使系统坍缩为平凡顺磁体，且。随着增大，缺陷弦通过三维连续相变凝聚，产生U(1)库仑相。

对于（），能量惩罚迫使自旋最大化振幅，在所有键上取态。因为二分金刚石晶格上的任何宏观平面与偶数条键相交，求和偶数个态将宏观磁通严格量子化为偶数，稳定了禁闭库仑相。在无磁单极子极限下，冰规则禁止孤立缺陷键，因此基本激发是键的闭合环。由于，这些环缺乏内部取向。由此产生的无向环气映射到三维伊辛模型的高温展开，将退禁闭转变置于三维伊辛普适类[18]，与独立的对称性论证和文献[16]的有限尺寸标度分析一致；环气估计[18]与文献[16]的数值转变点吻合至0.3%。

三、宏观磁通量子化与相图的完备性

当推广到任意时，增大的局域态空间提出了更复杂的相结构——如中间均匀背景或部分退禁闭级联——是否可能出现的问题。在本节中，我们始终在无磁单极子极限下工作。我们首先刻画两个极限区域的宏观拓扑，然后论证中间相被排除。

A. 大极限下的宏观磁通量子化

在大极限下，真空排他性地存在于极性的处。为表征宏观拓扑，我们计算跨宏观截面平面的全局极化磁通：

(5)

其中。

宏观平面完全贯穿二分金刚石晶格。与该平面相交的键总数严格为偶整数（例如，在含个晶胞的周期性晶格中，垂直于立方轴的平面为）。设为具有正投影极性（）的穿越平面的键数，为具有负投影极性（）的键数。我们有几何恒等式。投影极性的代数和为

(6)

因为是偶整数，量显然是偶整数，记为（）。将此拓扑约束代回宏观磁通方程得：

(7)

这个简单的宇称论证表明，大极限下的宏观磁通量子化为的倍数，对任意定义了一个禁闭库仑相。同样的事实曾在文献[16,18]中针对、在文献[19]中针对被注意到；我们在此记录其一般形式。

B. 小极限下的宏观磁通量子化

在相反的极限（），单离子各向异性在每个键上倾向于最小的可达自旋振幅。所得基态流形的性质定性取决于是整数还是半整数。

对于整数，单离子能量的唯一最小值为。在该区域，每条键被非磁性态占据，给出恒为零。热激发是稀薄的（逸度）并在冰规则下形成小的闭合环，因此宏观磁通保持在。系统是平凡顺磁体，具有唯一基态和局域（闭合环）激发。

对于半整数，最小振幅为，基态流形由所有满足冰规则的构型组成。这恰好是自旋冰，与原始值无关。代入大结果立即得到。任何半整数自旋的小相因此是完全退禁闭的U(1)库仑液体，拓扑等价于典型的自旋冰[11,14,35]。

综上所述，两个极限区域表现出以下宏观磁通扇区：对于整数，可达扇区从（小）变为（大）；对于半整数，从（小）变为（大）。在两种情况下，热力学可达磁通扇区的集合在两种极限之间不同，保证随变化时至少存在一个相边界：对于整数，中的非零扇区必须被激活（退禁闭），而对于半整数，中的扇区必须被抑制（禁闭）。剩余的问题是中间相能否介入这两个极端之间。

C. 禁止中间背景相

让我们评估由主要中间自旋振幅（，其中）组成的宏观均匀背景能否作为热力学稳定相出现。

在严格无磁单极子极限（）下，不同均匀背景之间的热力学竞争由局域单离子各向异性能函数决定。对于小区域（），这是严格凸的抛物线，在最小可达振幅处唯一最小化（整数为0，半整数为）。对于大区域（），的凹性将最小值置于边界处。

由于任何中间自旋态（）相对于最小或最大振幅态都承受广泛的能量惩罚，中间均匀和混合背景在热力学上均被禁止。系统必须在最小振幅和最大振幅背景之间直接跃迁。

D. 弦自由能层级与部分退禁闭的缺失

即使排除中间背景，人们仍可能设想在大区域内存在部分退禁闭转变的级联，其中拓扑荷递增的缺陷弦在不同的逸度下凝聚。此处，荷为的缺陷弦指一条连通弦，其上偏离其真空值个单位，满足；该电荷的等价磁通表述在第五部分A节给出。我们现在证明这一情景同样被排除。

在禁闭相（）中，离开真空的退禁闭转变由拓扑缺陷弦的热增殖驱动。携带相对整数拓扑荷的缺陷弦产生弦张力：

(8)

该张力在上形成倒抛物线。绝对最小张力严格在边界处实现：基本单位电荷（ ）及其对偶（ ）。 (9)对于 ，任何中间电荷缺陷（ ）都会带来严格更大的弦张力（ ）。（对于 ，唯一可用的电荷是基本弦 及其对偶 ，它们共享相同的最小张力，因此共同凝聚；不存在中间缺陷。）上述张力比较纯粹是能量层面的，而真正控制哪些弦增殖的量是单位长度自由能，其中还包含弦的构型熵。因此，我们以自由能的形式来表述论证。一条电荷为 长度为 的弦，其统计权重为 ，其中 收集了单位长度的熵因子：金刚石晶格上自回避行走的嵌入数，由联结常数 控制，以及局域顶点权重和内简并因子，这些因子在 无关的、量级为1的常数范围内有上下界。在泡令独立顶点近似下，每个这样的因子都是四配位金刚石晶格顶点处冰规则状态数的比值乘以一个至多为2的链简并度，因此其有界性与 和 均无关（第V A 2节对基本弦进行了显式计算）。关键的是，所有弦种类都位于同一晶格上，因此嵌入熵对所有 是共同的，而 中依赖于种类的部分是有界的。一条电荷为 的弦只有在 时才能增殖——即其单位长度自由能可以变为非正值。由于张力指数 在边界 和 处严格最小，而 仅在有界 范围内变化，增殖阈值最先被达到——即当 从受限相深处降低时，在最大的 值处——由基本弦达到。要让中间电荷弦先凝聚，其单位长度熵必须超过基本弦的熵，超出量随张力差 增长，而有界局域因子无法提供这一点。我们强调这一论证的逻辑地位：张力层次是精确的，而熵界则采用了与第V节中使用的相同且已验证的独立顶点（泡令）处理，该处理相对于文献[16,19]的蒙特卡洛数据在几个百分点水平上得到了验证。在这种受控处理下，部分退禁闭被排除；第VII节的模拟在两个交叉尺度之间没有显示出中间合作反常，这证实了上述结论（对于 存在的宽次级肩部是一个与尺寸无关的单离子肖特基特征，如第VII节所确立的）。我们还注意到，电荷2弦分裂为两条电荷1弦的熵过程有利于这一结论，因为它将高电荷弦转化为基本弦。一旦基本弦作为环绕环面上的缠绕弦增殖，每条这样的弦都会将横截面上的宏观通量平移 ，因此 ——一步到位地完全退禁闭进入U(1)库仑液体。更一般地说，电荷 缺陷弦的凝聚会将可及的通量扇区从 变为 ，因为每条缠绕弦将通量平移 。对于 ， 恢复 。由于基本弦首先达到其增殖阈值并且已经实现完全退禁闭，具有部分通量量子化的中间分数化相（例如 ，这需要 弦在 之前凝聚）不会出现。将此结果与两个极限中的宏观通量量子化（第III A和III B节）以及中间背景的禁止（第III C节）相结合，我们得出结论：对于整数 ，全局相图恰好限于三个相[平凡顺磁体、U(1)库仑液体、 受限相]；对于半整数 ，限于两个相[U(1)库仑液体、 受限相]。IV. 小 区域：精确对偶映射与自旋宇称二象性当 （ ），大自旋幅度在能量上受到惩罚，热力学性质定性依赖于 的宇称。我们通过对无磁单极配分函数的精确对偶变换来揭示这种整数与半整数之间的二象性[36,37]。在本节中，我们继续在 下工作。下面的变换是José等人[38]以及Dasgupta和Halperin[39]在分析 模型时所用的整数取值、无散度链流与紧致相位场之间的标准对偶（另见文献[40]中的教科书处理），此处逆向运行：我们从物理流 出发，重构相位场。同样的到 模型的映射在文献[16,18]中用于 ，在文献[19]中用于 ；下面的推导将此映射推广到任意 ，同时跟踪精确的链权重。A. 精确对偶表示从方程(3)出发，取精确极限 ，配分函数限制在最小阻挫流形上，其中 恒成立： (10)每个克罗内克δ通过拉格朗日乘子相位 来解析： (11)将其代入 并交换求和与积分的顺序，得到 (12)其中 。离散分部求和将指数中的顶点和转换为有向链 上的和： (13)其中我们定义了沿链的相位差为 。此变换解除了全局约束，允许对 的求和分解为独立局域求和的乘积： (14)其中我们定义了精确的局域链权重函数： (15)这个精确对偶连续表示的宏观统计力学因此由有效作用量 支配。B. 整数自旋：3D 普适性的保护对于整数自旋（ ），能量的绝对最小值唯一地定义了平凡非磁性真空态 ，其权重为 。基本热激发是 的基本缺陷，其权重被抑制为 。更高自旋态，如 ，带有被惩罚的权重 。对于整数 ，精确链权重 展开为 (16)将 与 一起展开，我们得到有效作用量： (17)主导项 与 无关（它适用于所有整数 ；见下文），并将小 区域映射到经典3D 模型[38]，其有效耦合为 。临界点 由 确定，其中 是金刚石晶格上的3D 临界耦合。我们对金刚石晶格上经典 模型的蒙特卡洛模拟（铁磁和反铁磁模型在这个二分晶格上共享相同的 ）得到 [41]，对应于 。求解 至 给出 ，与Pandey和Damle[16]对 得到的数值结果 吻合良好。所有非 谐波（ ）在3D 不动点处都是强无关算符[38,42]，因此小 区域对于所有整数 都属于3D 普适类。两点评论阐明了这里的逻辑。首先，3D 普适性最终是对偶作用量的精确U(1)不变性的结果，该不变性编码了实值无散度约束；谐波展开的作用是识别跃迁所连接的两个相（一个 顺磁体，对应于平凡相；一个 “铁磁体”，对应于库仑液体），并定量定位跃迁。其次，通过3D 普适性，我们指的是支配 处跃迁的临界不动点。对偶 模型的有序相不仅仅是辅助构造：其自旋刚度映射到库仑液体的涌现规范场的刚度，因此对偶描述在 两侧都具有物理信息性，尽管对主导谐波的截断仅对小 具有定量可控性。主导 耦合 对于所有整数 完全独立于 ；只有无关的高阶谐波（ ）对可用的自旋态敏感。因此临界点 对所有整数自旋是普适的。对于 ， 态的缺失使得 精确成立，因此高阶谐波系数在次领头阶上与方程(17)不同[例如， 获得额外的 修正]。这些差异仅限于无关扇区，不影响普适的 临界行为。C. 半整数自旋：无相变本小节的实际结论可以简单表述为：当 时，半整数 模型简化为普通的 自旋冰，因为最小二重态 是唯一存活的自由度。因此小 相是U(1)库仑液体，并且在小 处不可能发生热力学相变。下面分析的目的在于使这一表述定量化——特别是表明 中的领头修正不会产生任何相关（禁闭）微扰，而仅仅是重整化库仑液体的刚度。对于半整数自旋（ ），非磁性整数态 被严格禁止。每条链上可用的最低能态是物理二重态 ，两者都携带相同的玻尔兹曼权重 。更高自旋热态（例如 ）携带权重 ，被因子 大幅抑制。对于 ，评估精确链权重函数得到解析展开式： (18)提取因子 ，展开对数，并应用恒等式 其中 ，我们得到有效作用量（忽略常数）： (19)物理调控参数 作为整体全局归一化常数从前导作用量项中分解出来。虽然这移动了全局能量尺度，但不会消除系统的拓扑简并。前导的、与 无关的有效作用量 与经典 自旋冰的规范理论表述相同[11,13]。该作用量描述了两进两出构型的简并流形，该流形承载U(1)库仑液相，其特征是代数偶极关联和无禁闭整数通量（ ）。次领头 修正是全局U(1)对称的 转子耦合 ——与标准3D 模型作用量相同的形式。在对偶框架中，这一项严格保持连续U(1)对称性，仅重整化相位刚度[对应于涌现U(1)规范场的光子刚度]，而不会引入任何对称破缺的相关算符（例如磁单极势 ）来产生禁闭能隙。更广泛地说，3D U(1)库仑相对任意小的局域微扰本质上具有鲁棒性[11,14]， 微扰不会产生任何可能禁闭规范场的相关算符。值得注意的是，在纯 极限中， 光子刚度恒为零，因为 不包含高阶谐波。正是更高自旋态（ ）的混合——仅对 可用——通过虚涨落产生相位刚度。因此，对偶映射确认并量化了本小节开头所述简单图像：对于任何半整数 ，小 区域是 库仑液体，具有弱重整化刚度，且不发生热力学相变。一个技术性说明：任何半整数 的链权重 是反周期的， ——半整数自旋的对偶表现，让人联想到自旋 波函数在 旋转下获得的（ ）相位。对偶被积函数在 下仍保持单值，因为每次这样的平移恰好翻转 条入射链权重和 的符号；金刚石晶格的偶数配位数使得半角表示保持一致，配分函数严格为正实数。或者，从一开始就可以避免所有半角，将 视为整数值守恒流并将对偶相位场耦合到它；由此产生的链权重只涉及整数谐波 并导致相同的结论。我们保留了半角约定，以保持与物理自旋变量的显式联系。V. 大 区域：图形映射与退禁闭跃迁的性质在大 极限（ ）下，系统唯一地处于高度受限的真空态，其中所有自旋取最大绝对值 ，背景极性为 。缺陷弦的特征是自旋幅度减小 ，形成受局域冰规则约束的图。在本节中，我们再次在 下工作，推导配分函数的图形（环气）表示，并利用它们来确定退禁闭跃迁的性质。结果是一个尖锐的二象性：对于 环气与离散自旋模型的高温展开一致（ 时为伊辛模型， 时为三态波茨模型），而对于 它与连续3D 模型的高温展开一致。A. 缺陷图、冰规则相容性和构型熵 1. 冰规则相容性判据为了分析宏观缺陷统计，我们将物理自旋变量映射为离散有效通量。对于每条链 ，我们定义平移后的非负变量 。在这种表示中，最大幅度真空态对应于 （ ）。因此，任何承载基本缺陷（ ）的链都携带非零 通量 。在每个顶点处， 。冰规则（ ）严格地将两个子晶格上的第一项和固定为零，从而给出精确约束 。将此方程取模 ，可得到一个较弱的、纯粹局域的守恒定律： (20)由于 （模 ），冰规则蕴含方程（20），但逆命题不一定成立：模条件允许满足方程（20）但违反 的具有 的解。给定一个在每个顶点处满足方程（20）的缺陷链（即那些 的链）构型，自然产生一个问题：是否总能选择剩余真空链（ ）的取值，使得在每个顶点处严格的冰规则约束 成立？这样的选择称为冰规则相容真空赋值。答案揭示了一种拓扑二分法。冰规则相容性判据。考虑金刚石晶格上的任意一个 守恒缺陷构型及其任意一个顶点 。模守恒律单独保证——对于每个这样的构型和每个这样的顶点—— 的非缺陷链可以被赋予真空值 ，使得冰规则 在 处成立，当且仅当 。对于 ，只要 的少于四条链带缺陷，赋值仍然存在，并且只能在完全激发顶点处失败。证明。设 表示一个顶点处的缺陷键数（ ），并记 为它们的部分和。剩余的 个真空键各携带 ；若其中 个取值为 ，则冰规则 变为 (21) 该方程必须与 一起求解。 (i) （至少存在一条真空键，因此 可调）。逐情况讨论： ：平凡地 且 。 ： 永远不是 的倍数，所以 被 守恒本身所禁止。 ： ； 的唯一倍数是 本身，给出 ，这是可满足的。 ： 且 ，均在[0, 1]内，同样可满足。因此，每个满足 的顶点对所有 都容许冰规则相容真空赋值。 (ii) （所有键均被激发；无可调真空键）。由于 ，我们有 。没有真空键可用时，方程（21）没有自由变量，严格的冰规则要求 。模守恒律 （ ）当且仅当 是运动学允许区间 内 的唯一倍数时保证这一点（表III）。对于 ， 是唯一倍数，因此 守恒保证 。对于 ，区间内出现额外的 倍数，满足 守恒但给出 ——我们称之为冰规则不相容的情形。（我们在此刻意避免使用“单极子”一词：第六节中物理的、热激发的单极子在当前 设定下并不存在，而不相容性纯粹是模守恒律的算术性质。） 结合(i)和(ii)完成证明。 ? 表III. 完全激发顶点（ ）的冰规则相容性。区间 列出 的允许范围。冰规则要求 ；对于 ，区间内出现额外的 倍数，因此模条件不再保证冰规则（冰规则不相容）。自旋 区间 区间内 的倍数 相容？ 1 [4, 4] 4 ? [4, 8] 6 ? 2 [4, 12] 4,8,12 ? ? 该判据在顶点层面建立了鲜明的二分法。对于 ，模守恒律保证逐顶点与冰规则的相容性。必须记录一个注意事项：顶点级相容性本身并不保证全局一致的真空赋值，因为真空链在顶点之间共享。对于 ，此机制不出现（不存在三叉结点，也没有真空值被强制确定）。然而，对于 ，一个立方结点固定其唯一真空链的值（当三条缺陷电流全部指向内时 ，当它们全部指向外时 ），因此两个同向结点的强制真空链在一个公共通过顶点相遇的稀有排列不容许一致赋值，并被排除在物理系综之外。这种排除作为近程的、 对称的硬核约束作用于相邻结点之间；它首次出现在结点密度的二阶，其作用在V B节中说明。最小排除构型以及精确有限晶格验证记录在补充材料[25]的第五节中。对于 ，不相容性已在单顶点层面出现。作为具体例证，考虑 [图2(b)和2(d)]：一个具有 的顶点给出 ，而(3, 3, 3, 3)给出 。两者都满足 守恒但违反冰规则 ，产生电荷 。因此，对于 ，任何仅基于 守恒的环模型都会包含自旋冰模型所禁止的构型。与V C节中量化的缺陷逸度不匹配一起，这阻碍了到离散 自旋或规范模型的映射[37,43,44]。图2. 单个 金刚石晶格顶点处的冰规则相容性；箭头表示定义 的双部 链方向。(a) ：三个基本缺陷（ ）在单个顶点处湮灭， 处剩余的真空链完成冰规则 。(b) ：四个基本缺陷给出 ——模条件满足但冰规则被违反。(c),(d) 对应于(a)和(b)的自旋构型 ，分别具有 和 。四条键被画成平面十字形以便清晰；在金刚石晶格上它们指向四面体的角。 2. 构型熵单个缺陷链相对于真空背景的玻尔兹曼权重改变一个因子 。为了评估完整的环气体配分函数，我们还需要考虑背景自旋的宏观构型熵。我们采用保利的独立顶点近似[10]。在真空中，每个顶点在 个无约束二元赋值中具有 个有效冰规则构型，得到保利分数 。我们现在按照文献[18]计算比率 ，其中 表示在长度为 的缺陷环 存在时满足冰规则的背景自旋构型数， 是无缺陷真空中的相应计数。对于 ，缺陷保留非零振幅 [ ]，因此每条缺陷链携带一个内部箭头（ 的符号），使环成为有向的。对于具有固定取向的这种有向环，每个环顶点处的两条缺陷链是完全确定的——它们的 值由环的通量电荷固定——只留下两条背景链作为自由变量。每条背景链携带 ，冰规则 约束该对：如果缺陷链贡献 给顶点和，则两条背景链必须和为 ，这由一条链携带0而另一条携带 唯一满足。恰好存在两种这样的赋值（任一排序），因此在 个背景构型中，恰好两个满足冰规则，给出顶点分数 。这对所有 都一致成立。此外，每条缺陷链用一个固定的单态替换一个二态真空链（ ），每条缺陷链贡献一个链简并因子 。结合两个因子，每个环链的保利熵减少为 (22)与每条缺陷链的裸玻尔兹曼权重 一起，有效基本弦逸度为（详细推导见补充材料[25]的第二、三节） (23)每个连续环所承载的内部自由度 取决于 。对于 ，缺陷是 ，产生无向环（ ）。对于 ，如上所述，非零缺陷振幅保留位移的符号，产生手性有向环，其承载 个不同物理态。B. 到三态波茨模型的映射 对于特定情况 ，我们的冰规则相容性判据保证几何 晶格约束和 通量守恒相匹配。除了V A节记录的短程结点-结点排斥外，每个 守恒缺陷图在物理冰规则流形内都是可实现的，下面的环展开是一个离散三态自旋模型的展开。基本缺陷 自然地作为共轭弦，携带通量 和 （模3）。严格的 守恒在几何上允许度为2的通过、度为4的交叉，以及唯一地，度为3的结点，其中三个基本弦汇聚 [ （模3）]。这种几何结构容纳了恰好三个拓扑弦在单个顶点处的同时局域湮灭 [图2(a)和2(c)]。顶点权重由每个顶点度数处的冰规则态计数确定。在无缺陷真空（度数 ）中，每个顶点在 种可能性中具有 个有效构型，给出局域态和 。在度为2（通过）顶点处，两条缺陷链由环方向固定，剩余两条真空链必须满足冰规则，给出 。在度为3（立方结点）顶点处， 守恒律已经迫使三条缺陷电流共享一个取向（ ；即全部向内或全部向外）；冰规则另外固定唯一剩余真空链的值，这使得有效构型数相对于度为2的基线减半，给出 。在度为4（交叉）顶点处，所有四条键都被缺陷占据（没有剩余真空链）；有向冰规则然后唯一地固定缺陷构型，给出 。因此，相对于度为2的权重归一化的顶点逸度为 （立方结点处的惩罚）和 [交叉处的增强；分母为 是因为一个交叉取代两个独立的度为2顶点]。利用标准拓扑图恒等式 ，完整的无单极子配分函数形式为 (24)其中 是真空扇区的配分函数，其中每条链携带受冰规则约束的最大振幅 。每个真空构型贡献一个玻尔兹曼权重 ，保利近似（V A 2节）给出冰规则简并度 ，因此 。求和遍及金刚石晶格上的所有图 ，其边携带有向 通量，满足每个顶点 处的守恒律 （模3）[方程(20)]，排除了V A节中识别的短程结点-结点构型； 是总边数， 和 分别是度为3和度为4的顶点数， 是闭合缺陷股数，每股携带两个手性取向（在度为4交叉处，四条入射缺陷链配对成两条贯通股，因此在顶点处交叉的两个环计为两股）。该图系综同构于金刚石晶格上经典三态波茨模型的高温图形展开[36,45]（详细推导见补充材料[25]的第三和第五节）， (25)在 上的标准傅里叶分解将每个键权重展开为 其中 ， ， 。对自旋变量求和，每个顶点处的 通量守恒恰好产生与上述相同的图系综，给出 (26)与自旋冰展开的同一性由对应关系 确立；局域几何权重 和 仅充当裸耦合的短距离重整化，不改变普适类。我们阐明该映射的状态。两个图系综之间的对应关系由 守恒律和顶点级相容性判据固定，直到V A节记录的显式短程排斥：两个同向结点的强制真空链碰撞的构型出现在波茨系综中但不在自旋冰系综中。相比之下，顶点权重在保利独立顶点近似内计算：背景冰流形的相关性在独立顶点水平之外诱发缺陷弦之间的额外相互作用。特别地，一个立方结点，将其单条背景链强制为 ，充当背景高度场的局域微扰（有效电荷），因此结点经历由背景库仑相介导的熵致偶极型相互作用。这些相互作用在有限缺陷密度下是短程的或被屏蔽的——包括在最短距离上，上述严格的结点-结点排斥充当硬核——并且关键的是，它们尊重环气体的 对称性；因此它们重整化有效波茨描述的耦合，但不影响下面朗道论证所依据的对称性分类。它们的数量小性由同一近似在其他地方的准确性事后支持：它再现文献[16]的 转变点 至0.3%以内[18]，以及文献[19]的 转变点至 以内（下文）。最后，对于 ，波茨和时钟描述等价。由于 上的克罗内克δ满足 ，三态波茨模型等价于三态时钟模型——恒等式 确保所有非相同自旋对获得相等权重，使时钟相互作用等同于波茨克罗内克δ。此等价性对 成立，但对 失败，因为不同角间隔产生不同权重（参见V C节）。在连续朗道-金兹堡-威尔逊有效作用量中， 缺陷弦融合产生一个对称允许的三次不变量（ ）。在朗道平均场理论中，三次不变量驱动一级相变，并且在三维情形下这一预期得到证实：事实上，大量蒙特卡洛研究已经确立三维三态Potts模型经历一级相变[45,46]。因此，上述图映射为Pandey等人对[19]所数值观察到的一级相变提供了理论解释。我们注意到文献[19]本身通过一种不同的、基于对称性的路径为一级特征提供了理由：通过映射到配分函数，并观察到大的相展现出长程序但不展现，该转变在那里被等同于三态自旋系统的转变。这里给出的图解推导从一个互补的微观出发点得出相同结论，并且还给出了转变点的定量估计，如下所述。

一级转变的位置可以通过Bethe-Peierls（空腔）近似[47]解析估计。其思想是用具有相同配位数的Bethe晶格（Cayley树）替代金刚石晶格，在该晶格上Potts模型的配分函数可以通过单变量空腔递推精确求解。在这棵树上，每个位点有“子代”；空腔消息——边界自旋匹配参考态的概率与任一非参考态（非参考态由Potts模型的对称性等价）概率之比（ ）——满足自洽方程（27），其中是Potts耦合，且。

无序解总是存在；对于Potts模型，第二个解（有序相）不连续地出现——这是由三次不变量驱动的一级转变的标志。令两个解的Bethe自由能相等（见补充材料[25]第八节）即可闭式确定共存点：，即，其中是通过与自旋冰逸度相关联的Potts高温参数。转换回自旋冰逸度给出闭式估计（28）（注意文献[19]用表示与我们的约定对应的量）。这仅仅高估了蒙特卡洛值[19]一个因子，这是Bethe近似在晶格上预期的精度：在一级相变时关联长度保持有限，所以Bethe晶格上缺失环——其唯一结构缺陷——只有微小的定量影响（完整推导包括自由能评估见补充材料[25]第八节）。

C. ：三维 普适性

1. 钟映射的失效与连续U(1)基线

我们首先回顾一个已在第四节使用的事实：冰规则是一个实值约束——晶格散度精确消失，而非仅模整数——因此微观模型对每个都具有精确的U(1)结构。这正是此处及文献[16,18,19]中对偶表示的基础。因此，本小节所讨论的问题并非U(1)对称性如何涌现，而是相反：大的缺陷系综是否将此微观U(1)结构约化为离散结构，正如它对(Ising)和(Potts)所做的那样。我们证明对于，它不约化：没有离散钟模型能重现环气，且转变由微观U(1)结构主导；即，它属于三维类型。

对于，情形相对于发生定性变化。如冰规则相容性判据（第五节A）所示，高斯定律不再与物理冰规则一致（冰规则不相容性）。此外，缺陷张力的层级结构造成键逸度失配。自然的比较模型是钟模型[38]（而非Potts模型），因为缺陷荷自然对应于钟模型的电流谐波，这些谐波区分不同的角分离——而Potts模型平等处理所有非相同状态，缺乏这一结构。在标准钟模型中，高电流谐波产生多项式逸度衰减（高温下），而在自旋冰模型中，每个缺陷能级带有独立的化学势，其逸度（ ）指数衰减。

一个关键推论是三次交点的算术禁止。对于（ ），三条基本弦满足（模3），允许映射到驱动一级Potts转变的对称性破缺不变量（第五节B）。然而，对于（ ），（ ），因此仅由基本弦构成的三次交点被守恒律本身严格禁止。这种算术障碍内在地抑制了本应驱动一级转变的不变量。

第二个纯粹几何的障碍出现在基本弦的湮灭中。在钟模型中，模守恒律（ ）允许在时多达条基本电流在单个顶点合并；特别地，对于（ ），所有四条基本电流局域湮灭，无中间键，罚值为1（表IV）。即使对于，仅条携带奇电流的中间键的短级联就足够（见补充材料[25]第六节）。

表IV. 在金刚石晶格上湮灭条基本弦所需的中间键数和总级联罚值（ ）。每行显示连接融合顶点的中间键数及其逸度乘积（“罚值”）；罚值为1表示无需中间键。对于两个模型一致且无需级联。对于它们发散：钟模型至多需要条多项式抑制逸度的键（补充材料[25]第四节），而自旋冰模型需要条指数抑制逸度的键；详见正文。

钟 自旋冰 自旋 键数 罚值 键数 罚值
1 0 1 0 1
0 1 0 1
2 0 1 1 1
2
3 1 3
2
4

相比之下，在自旋冰模型中，精确冰规则比模条件更具限制性：它甚至禁止同时合并对于[图2(b)和2(d)]，更一般地，禁止任何顶点吸收超过两条基本弦（当一个顶点有三或四条基本弦时，对于）。因此，基本弦两两融合成复合股线，这些股线通过包含至少条电荷中间键的融合树进一步融合；线性链，其中每次吸收一条基本弦（；图3），是最简单的此类树。

图3. 对于的最简单（线性链）融合通道：在每个融合顶点，一个入射的基本缺陷（ ）与累积的中间缺陷合并，使有效电荷增加一。箭头表示沿其定义的二分子链取向；空心和实心圆分别表示-和-子晶格顶点。显示的链标签沿链在之间交替，因为顶点在两个子晶格之间交替。显示的是半奇数情形；对于整数，相同结构适用且链取向均匀。为清晰起见，链画成直线；在金刚石晶格上，骨架沿之字形路径。

这些组合障碍——冰规则不相容性、指数键逸度失配、以及延伸的级联几何——意味着，与不同，大环气不是任何离散钟模型的高温展开。能够区分离散守恒律与完整U(1)守恒律的唯一过程是融合级联本身，我们现在量化其统计权重。

定量地，每条中间键引入逸度，因此融合通道的总统计代价——桥罚值——是其中间逸度的乘积（见补充材料[25]第二节的自足推导）。对于线性链，（29）。其他融合树带有略不同的指数[对于，最廉价的通道将基本弦两两融合成三条电荷2键，每条逸度，代价为，而非链的——这一最小值由以(4,1,1)终止的成对加倍通道共享；见图4]，但以下界与通道无关。每次成对融合恰好产生一条中间键，因此从条基本弦开始并在吸收条股线的顶点终止的通道含有条融合输出；其中至多一条——带有的共轭基本弦，只能在成对终止中被消耗——落在计数范围之外。因此每条通道至少含有条被计数的中间键，且由于张力指数在此范围边缘取极小值，每条至少花费。因此任何通道的罚值受限于（30）。

图4. 对于（ ）的最廉价湮灭通道：六条基本弦（ ）在三个顶点处成对融合，链值为(1,1,4,6)，三条电荷2的中间键（显示的链值为）在一个公共顶点处湮灭，链值为(4,4,4,0)，其真空链吸收总电荷。箭头表示沿其定义的二分子链取向；实线和虚线分别表示缺陷链和真空链；空心和实心圆分别表示-和-子晶格顶点。通过整体翻转得到等价构型。

对于，；对于，；且抑制随迅速增强，罚值在固定下至少以速度衰减。我们还注意到，指数精确在时消失——任何湮灭通道中中间键的最小数目为——恢复在处无抑制的环对湮灭和无抑制的三次交点。U(1)结构的指数保护因此精确在处启动，将第五节A的 vs 二分性编码在单一公式中。

相比之下，在钟模型中，较弱的模守恒（ ）（而非精确冰规则）允许定性更短的级联。对于（即），所有条基本电流在单个顶点（ ）处湮灭，因此无需中间键且罚值为1。对于，需要级联，但比自旋冰模型中的更短。例如，在（ ）处，三条基本电流（每条）在第一个顶点合并，产生电流的中间键（使用四条链中的三条）；在第二个顶点，此中间键遇到另外三条基本弦，且（模6）以仅一条中间键（携带，逸度为）终止级联，而自旋冰模型中有条中间键。一般地，钟模型级联仅需条多项式抑制逸度的中间键（补充材料[25]第四节）。两个级联代价的比较见表IV。

由于级联被自旋刚性指数抑制，在领头阶，活跃宏观图系综被剥离了度3分叉交点，并坍缩为形成闭合有向环和度4环交叉的基本弦。此受限系综在拓扑上同构于连续经典三维模型[42]的高温图解展开，仅在短距离顶点权重上有差异。

2. 配分函数的分解与三维普适性

为了将此连续U(1)基线转化为对于三维普适类的受控解析论证，我们构造真空归一化配分函数的分解，其中是限制在真空扇区的配分函数，该扇区中每条链携带满足冰规则的最大幅度（即无缺陷的-束缚相）。因此比值衡量所有缺陷构型相对于无缺陷背景的总统计权重。我们将分解为连续U(1)扇区和指数抑制的离散修正。

为了将U(1)守恒扇区与级联修正分离，我们注意到在交叉顶点（度4），自旋冰相对于两条独立穿过弦带有常数Pauling权重。因此我们定义装饰模型，记为，它采用冰模型的键逸度用于电流，并在每个顶点精确U(1)电流守恒（与冰规则相同的高斯定律），同时额外将每个交叉顶点增强一个因子。由于是独立于和的常数，交叉装饰作为短程四次复合算符，在重整化群（RG）下流向零。由于采用与冰模型相同的键逸度，自旋冰环气的每个非级联图在两个配分函数中具有相同权重，权重在此处及本节其余部分都在独立顶点近似内评估。还有两个剩余差异需要追踪。首先，未被装饰模型捕捉的自旋冰构型是上面分析的融合级联，其总贡献由桥惩罚控制[式(29)]，并且在 上呈指数小。其次，装饰 系综包含一些没有自旋冰对应物的局域电流构型[例如，一个具有电流 的顶点，它守恒U(1)电荷但与冰规则不相容，冰规则允许任何顶点处至多两条正偏差链接]；排除它们是对顶点权重的U(1)对称短程修正，与交叉装饰处于同等地位，并在下面的RG论证中如此处理。通过这种记账方式（见补充材料[25]，第七节），配分函数分解为 (31)该分解的物理内容如下。装饰 模型 捕捉自旋冰环气的U(1)守恒扇区，并共享三维 普适类（如下所示）。级联贡献 是唯一对 守恒律离散性敏感的项，其裸权重被晶格几何指数抑制。结合RG论证，该分解将 置于三维 普适类中，具体如下。首先，在装饰 模型中引入的交叉装饰算符 在每个顶点处保持连续U(1)电流守恒，因此对应于一个U(1)对称的无关微扰。在三维 临界不动点处，主导修正标度指数为 [48]，因此该算符的标度维数为 。同样，冰模型键逸度 与贝塞尔比值 之间的差异是涉及高电流算符（ ）的U(1)对称短程微扰，它们在 不动点处更是强无关的。由于无关性，RG流保证装饰 模型流向三维 普适类。其次，级联修正 是区分离散 守恒律与连续U(1)守恒律的唯一机制；在粗粒化 描述中，它由主导离散各向异性 表示。现在两个逻辑上独立的事实结合。（i）RG无关性：在三维 不动点处，各向异性 对所有 具有标度维数 （ 来自蒙特卡洛RG研究[49,50]，且 随 单调增加；此处 ）。因此它是一个（危险的）无关算符[38,49]：对于 不动点吸引域内的任何裸耦合——特别是对于任何足够弱的裸耦合——它不能改变临界普适类。（足够强的裸各向异性可以逃出该吸引域：硬 钟模型，即 处的无限各向异性极限，分解为两个退耦的伊辛模型，并以伊辛而非 指数有序；这正是以下要点重要的原因）。（ii）裸小性：分解表明该算符的裸权重不仅小而且是指数小[ ，式(29)和(30)]，将系统置于吸引域的深处。它还控制交叉标度：超越该长度尺度后渐近 行为变得可见的尺度很短，因此即使在中等规模模拟中也不应观察到离散性的残留。综合（i）和（ii），将 烧绿石自旋冰的解禁闭相变置于三维 普适类中，且渐近临界行为提前出现。 3. 通过O(2)环气映射定量估计 由于离散各向异性被高度抑制且在 处无关，系统在临界点附近渐近趋于连续O(2)对称性。因此，我们可以通过将有向稀释环气映射到金刚石晶格上经典三维 模型的高温图解展开，定量估计解禁闭相变点[42]。相变发生在基本弦逸度 达到临界逸度 时，其中 是金刚石晶格上三维 临界耦合。令 ，我们得到一般标度公式： (32)根据我们的蒙特卡洛测定[41] [即 ]，我们得到 和 ，给出 。当 时，临界逸度渐近于 （即 ），与经典连续自旋极限一致：当 时离散自旋变量 趋近连续变量， 通量量子化变得无限精细，离散 禁闭相消失；本模型中任何弱易轴各向异性，即 略低于 ，足以稳定连续U(1)库仑液体。由上述分解确立的 融合级联的几何抑制确保该公式所基于的连续O(2)近似对所有 保持准确。无单极极限下获得的相、通量量子化和临界现象的完整分类收集在表I和表II中，并在图1中示意性说明。VI. 有限温度效应：磁单极、弦切割与交叉上述各节的分析依赖于严格实施局域冰规则 。这种无单极极限在 取固定 的极限下实现，即在固定 下实现（第二节C）。在任何实际有限温度 下，有限交换 允许磁单极的热激发[8,29]，其特征为一个小但严格非零的逸度 。在本节中，我们回到式(3)中生成的精确配分函数。一般而言，三维拓扑序对热涨落是脆弱的[33,51,52]。基于这一原理，我们具体地表明，单极的热激发显式地破缺了散度自由约束所编码的微观U(1)不变性——在小 区域的对偶描述中作为对称破缺场进入——并在大 区域作为缺陷弦的拓扑切割作用[18,23]，破坏宏观拓扑不变量并将相变圆滑化为对所有 的交叉（除 外）。A. 小 区域：对称破缺正弦-戈尔登场的出现本小节的构造是标准的[36,38]，我们为自洽性简要地保持推导；热单极在对偶 语言中作为有效外场作用的结论在文献[18]中针对 得出。然而，我们注意到这种对偶观点本身是最近才有的：自旋冰中的热单极传统上通过连续静电类比处理，在德拜水平插入屏蔽[8,14,29]。相比之下，在下面的精确对偶表示中，电荷 的单极作为局域顶点算符 出现，因此完整的单极气体——对所有 同时——作为显式U(1)对称破缺位点权重进入配分函数，屏蔽物理由此产生而非被假设。为正确捕捉小 极限下热单极的效应（第四节），我们通过恢复局域散度求和来放宽配分函数中的约束。对偶相位场表示中的配分函数变为 (33)遵循先前引入的精确离散分部积分，精确自旋扇区退耦并如前映射到局域链接权重 。新增加的是对每个位点局域离散单极电荷 的无约束独立求和，我们将其定义为局域位点作用量 ： (34)因为每个金刚石晶格顶点恰好连接四条链接，局域散度 是四个自旋变量的代数和。即使对于半整数自旋，四个半整数的和也是整数，所以对所有 有 。在低温区域 中，该求和由基本单极激发 主导。截断至领头阶， (35)取对数，我们发现热单极在每个晶格顶点产生局域正弦-戈尔登势。对于整数自旋（ ），有效作用量计算为 (36)生成的项 作为沿涌现三维 模型的 方向施加的强度为 的均匀外磁场。该场显式破缺连续U(1)对称性，是一个相关微扰，使戈德斯通模带隙并消除热力学奇异性。因此，三维 相变在任意有限温度 下圆滑化为平滑交叉[18,23]。对于半整数自旋（ ），其中未微扰系统处于U(1)库仑液相，有效作用量获得相同的项。这证实了标准的德拜-休克尔屏蔽机制[29,53]：单极等离子体引入一个屏蔽长度，指数地截断代数偶极关联，将库仑液体平滑转化为平庸顺磁体而无相变。B. 大 区域：推广有向环气与拓扑弦切割现在我们对称地将大 有向环气展开推广到有限温度。热激活局域单极（ ）的几何后果是缺陷弦不再被拓扑力迫成闭合连续环。有向弦现在可以终止，产生对应于热单极位置的开放端点 [54,55]（图5）。图5. 金刚石晶格上的缺陷弦拓扑。实线表示携带 或 的缺陷链接（与亚晶格约定交替）；虚线表示真空链接（ 或 ）；箭头表示定义 的二部 链接取向；空心和实心圆分别表示 和 亚晶格顶点。(a) 在 处，缺陷弦在每一顶点处满足 的闭合有向环。(b) 在 处，热单极（ ）作为弦端点，将闭合环切割成开放段。六边形环在二维中示意性绘制；在真实晶格上它是一个椅形六元环。环气配分函数现在对所有金刚石晶格上的有向图 展开，包括闭合环和由单极端点终止的开放弦段。对于 ，端点顶点权重计算如下。在单极端点处，四条链接中的一条携带基本缺陷（ ，固定），而其余三条链接处于真空扇区（ ）。单极条件 （即 ）约束三条真空链接中有多少条携带 与 。在无缺陷真空中，每个顶点有 个有效冰规则构型（恰好两条链接在 和两条在 的构型）。在这6个构型中，恰好 个与一条链接被提升为缺陷同时保持 相容（见补充材料[25]第二节）。因此端点顶点权重为 。通过用于结和交叉权重的相同归一化约定（ ），端点修正因子为 。包括单极逸度 ，每个端点贡献一个因子 ，给出推广的稀释环气配分函数： (37)其中 是边总数， 是单极端点数（满足 的弦终止位点数）， 是股数——闭合环和开放弦——每条在有向金刚石晶格上携带两种手性取向。该图解展开精确映射到存在均匀外磁场 的经典O(2)（即三维 ）自旋模型的高温展开。在标准O(2)高温展开中，外场 为每个弦端点赋予权重 。令 立即得到有效场强： (38)该结果推广到 （ ，无向环映射到三维伊辛模型）。与 处的有向环不同， 缺陷 是无方向的： 在 下是其自身的共轭，所以一条缺陷链接可以终止于 或 单极。这使每个端点处的有效真空构型从 翻倍到 ，因此端点权重从 到 （详细的 推导见文献[18]）。产生的端点因子 （相对于有向情况 ），精确补偿了伊辛约定中的减半磁耦合 [ ]，给出相同的普适有效场： 。因此，无论 如何，热单极作为施加在涌现规范变量上的均匀外磁场作用。该场切割缺陷弦，允许禁闭拓扑环在体相中终止。第三节建立的宏观通量量子化定律因此被破坏，使所有拓扑扇区绝热连通[18,56]。C. 拓扑相变的阶依赖命运该弦切割场的热力学后果取决于无单极相变的阶。一方面，对于连续相变（ 和 ），任何与序参量共轭的有限对称破缺场破坏二阶临界点[33,57]；因此连续解禁闭相变对任意 被圆滑化为平滑交叉。另一方面，对于一级相变（ ），情况在本质上有所不同：由于该相变涉及被自由能势垒分隔的两个宏观极小值之间的不连续跳跃（由相关的三次不变量驱动），因此预期在弱对称性破缺场下共存线仍会持续存在。因此， 所特有的一级相变预计将在有限温度下作为真正的热力学相边界存续[图1(b)]。这种情况类似于液-气相变：在零外场（ ）下，自由能具有两个被势垒分隔的简并极小值（有序相和无序相）。弱对称性破缺场（ ）会使自由能景观倾斜，但无法消除势垒，因此共存线应持续存在。随着磁单极子逸度 增加，自由能势垒缩小，并在有限的临界磁单极子密度 处消失，此时两个极小值合并为一个。在这一点上，一级相变线终止于一个临界终点，类似于液-气临界点，并预期属于三维伊辛普适类。

为估计 ，我们将第V.B节的Bethe-Peierls分析扩展，纳入磁单极子诱导的对称性破缺场 ，该场以 的形式进入[参见方程(37)]。外场将腔递归修正为 ，其中前因子 使腔消息偏向场偏好的状态。临界终点——即自由能势垒恰好消失的点——由迭代映射 上的三个同时条件确定：(i) （腔消息的自洽性），(ii) （有序和无序不动点已合并，即边缘稳定性），以及(iii) （ 的拐点消失，标志着中间势垒的消失）。这三个关于三个未知量的方程可以闭式求解，得到 和 ，对应于 或 （完整的代数推导见补充材料[25]第VIII节）。

的较小值反映了 处一级相变的微弱性：由于 是Potts模型在三维空间中表现出三级相变的最小整数，自由能势垒本质上很浅，容易被弱扰动破坏。事实上，Bethe-Peierls近似系统性地高估了势垒高度，因为它忽略了进一步侵蚀浅自由能景观的长波涨落。简单立方晶格上三维三态Potts模型在外加有序场下的蒙特卡洛模拟[58,59]已确定一级共存线终止于三维伊辛普适类中的临界终点，对应的临界场为 ——比Bethe-Peierls估计 小一个多数量级。在与当前烧绿石问题相关的金刚石晶格上，由于配位数较低（ 对 ）放大了热涨落，真实临界场预计会更小。

该临界终点提供了独特的实验特征：在候选的 烧绿石材料中，比热和磁化率应在孤立温度处表现出伊辛型临界标度行为，为底层弦融合物理提供了直接的热力学探针。值得注意的是，与常规 自旋冰中需要精细调节外磁场的液-气型临界终点[8]不同，该临界终点严格出现在零外场下。对称性破缺场由热磁单极子等离子体内在产生，使 烧绿石自旋冰成为罕见的系统示例，该系统仅通过热涨落自发产生其自身的拓扑临界终点。

VII. 有限温度下的蒙特卡洛模拟

为在有限温度下检验前几节的解析预测，我们对哈密顿量方程(1)在有限烧绿石晶格上进行经典蒙特卡洛模拟，晶格包含 个自旋（其中 为金刚石晶格顶点数）并采用周期性边界条件。本节中我们固定 作为能量单位。

A. 模拟方法

标准蒙特卡洛扫描包含 次单自旋热浴更新。在 次热化扫描后，在 次扫描上累积测量值。

主要热力学观测量是每个自旋的比热，
(39)
其峰值追踪交叉标度。为明确区分集体拓扑现象与非合作局域单离子物理，我们系统计算每个自旋幅度的热占据分数，
(40)
即携带幅度 的键所占比例，对所有 条键和蒙特卡洛构型求平均。

B. 整数自旋（ ，2，3）

我们在低温（ ）下计算 时 为2和3的比热随 的变化，系统尺寸取 ，4和8；图6(a)和6(b)分别显示 和 的结果。在所有三种情况下， 呈现显著峰值，其位置与解析预测的交叉标度 和 合理一致，其中 和 分别为第IV节和第V节中估计的无磁单极子临界逸度。值得注意的是，峰值高度不随 增大而发散，证实了无磁单极子相变在有限温度下被圆滑化为平滑交叉，如第VI节磁单极子屏蔽分析所预测。对于 ，结果与文献[18]的详细数值研究完全一致，该研究表明比热峰值在热力学极限下趋于有限值，并且 退禁闭转变（映射到高压水冰中冰VII到冰X的转变）是由Debye-Hückel磁单极子屏蔽驱动的连续交叉。

图6. 整数自旋在 时的蒙特卡洛结果， ，4，8；交换耦合固定为 ， 和 以 为单位。(a),(b) 和 时每个自旋的比热 ；插图放大了 峰区域。垂直虚线标出预测的交叉标度 （紫色）和 （暗黄色）以及肖特基位置 （青色）。(c),(d) 和 时的热占据分数 （ ）。(e) 跨越 ，2和3的比热（ ）。

对于 ，除合作交叉峰外，比热在 附近还表现出宽泛的次级肩部，该肩部在 时不存在（因为 只包含两个幅度， 和1，所以没有中间退布居步骤）。如第III节完备性论证所确立，该特征并不对应中间热力学相。相反，它是由单离子各向异性驱动的非合作肖特基反常。两个最高幅度态 和 之间的单离子能级分裂为 。对于 ，最大幅度态 是该二能级系统的激发（较高能量）态。将其近似为独立的二能级系统，其激发态 随 增大而热退布居，肖特基比热在 处达到峰值，得到估计峰位 。该解析估计准确捕捉了蒙特卡洛数据中观察到的次级鼓包的位置和 标度：随着 增大，鼓包向 移近[图6(e)]。此外，显式计算的微观自旋占据数 [图6(c)和6(d)]表明该肩部恰好与从最大幅度态 （ ）到较低自旋幅度（ ）的连续热退布居重合。该肩部缺乏有限尺寸标度，结合占据数数据，清晰地解耦了局域单离子熵释放与集体拓扑现象。

C. 半整数自旋（ ）

对于半整数自旋，第IV节的解析框架预测在小 区域无相变，因为半整数键权重直接映射到 库仑液体。因此唯一的非平凡拓扑交叉是 大区域退禁闭交叉，其无磁单极子标度由 设定，对于 为 （我们的闭式Bethe-Peierls估计，第V.B节；文献[19]测得的转变点为 ），对于 为 [方程(32)]。

在 （ ）和 （ ）下的蒙特卡洛比热证实了该图像[图7(a)和7(b)]：与拓扑交叉相关的尖锐主峰出现在 附近，其高度随 增大而饱和； 对应的峰在 处示于图7(f)（ ）。在大正 （小 ）极限下，单离子各向异性 将每个自旋限制在其最小双态 上，模型约化为有效 伊辛自旋冰，其交换耦合被 重新标度，得到有效温度 。图7(a)和7(b)中的水平虚线显示了独立计算的该有效 模型的比热，蒙特卡洛数据在大 下对所有半整数自旋收敛到这些值，为有效哈密顿量映射提供了直接数值确认。此外，与整数自旋情况类似，对所有 在 侧都清晰可见更宽的次级鼓包。正如分数占据数 [图7(c)和7(d)]所确认的，这些占据数追踪热退布居级联，这是由 驱动的相同非合作肖特基反常。根据我们的标度公式 ，该肖特基峰逐渐向 移动，并随着 增大与拓扑交叉峰合并[图7(f)]。

图7. 半整数自旋的蒙特卡洛结果；交换耦合固定为 ， 和 以 为单位。(a) 在 和(b) 在 时每个自旋的比热 ， ，4，8；插图放大了峰区域。垂直虚线标出预测的交叉标度 （紫色；对于 ， 取自表II）和肖特基位置 （暗黄色）；水平虚线（青色）标出有效 伊辛自旋冰在 时的比热。(c),(d) 和 时的热占据分数 （ ）。(e) 在 和0.4时的比热（ ）。(f) 在 时跨越 ， 和 的比较（ ）。

对于 ，无磁单极子极限（ ）预测了由对称性允许的 Potts三次不变量驱动的一级相变。然而，在我们模拟可及的有限温度下，比热峰没有随 表现出热力学发散迹象，能量Binder累积量也没有显示一级相变特征的负凹陷。这意味着一级共存线要么在模拟温度以下终止于临界终点，要么需要指数更大的系统尺寸来分辨弱一级性质。事实上，在 和 下比热的比较[图7(e)]表明，虽然峰在冷却时变尖，但仍保持交叉行为。定量而言，在 时，磁单极子逸度 比第VI节Bethe-Peierls分析估计的临界终点值 大一个数量级（见补充材料[25]第VIII节）——并且考虑到Bethe-Peierls近似对 的高估相对于三维三态Potts模型精确蒙特卡洛结果[58]超过10倍，很可能比真实晶格值大两个数量级。因此，模拟明确地将系统置于交叉区域内，证实热磁单极子已完全洗掉了宏观相共存。我们强调相应的局限性：由于所有可及温度都高于估计的临界终点，本模拟测试的是交叉区域，但未直接验证预测的一级相变线在较低温度下的存续。

受此局限性推动，我们另外在较小各向异性下对物理（ ）模型进行了平行回火模拟，此时预测的转变温度低于估计的临界终点；所探究的问题不是一级相变的存在性——这在 处已在文献[19]中确立——而是其在 处预测的存续。这些运行显示出与弱一级相变一致的迹象——比热峰在预测位置随系统尺寸变尖并增长、能量Binder累积量的负最小值以及初始双峰能量分布——并且随着 减小（此时转变处磁单极子逸度更小），这些迹象变得更清晰。然而，能够明确鉴定的有限尺寸标志在可及尺寸下尚未分辨；决定性测试需要更精细的采样，留待未来工作。

我们以评论直接模拟无磁单极子极限的前景来结束本节。对于 自旋冰，在 处采样冰规则流形通常使用非局域环更新[9]，文献[16,19]正是使用此类蠕虫算法来描绘 相图——包括与两个转变的三维 和三维伊辛普适类一致的有限尺寸标度分析——并确立 处的一级相变。相同的策略可推广到任意 ：沿闭合回路将 移动 的蠕虫型更新（尊重 ）在冰规则流形内是遍历的。对于 ，其中我们的理论在方程(32)的耦合 处预测三维 临界性，并对整数 在 处预测与 无关的小 转变，此类模拟尚未进行；它们将为这些预测提供直接检验，留待单独研究。在预测的临界终点（，其中）附近的共存线的有限温度命运，同样可以通过补充了磁单极插入的环采样来探测。VIII. 实验展望在得出结论之前，我们讨论模型方程(1)及我们预测的实验背景。A. 高自旋烧绿石磁体方程(1)中的单离子项并非理论上的理想化：每个烧绿石磁体的磁性位点都具有点对称性，其三角轴沿局域方向，因此对于任何具有的离子，主导的对称性允许的单离子各向异性恰好是沿局域轴的。与交换耦合相当的各向异性在两侧都有记录：对钆烧绿石的电子顺磁共振测量发现，其易面单离子项（）约为[60]，而对烧绿石，通过从头算量子化学计算得到了易轴项（），并已应用于磁激发共振非弹性X射线散射分析中[61,62]。高自旋烧绿石磁体现在已覆盖整个自旋值范围——最系统性地体现在氟化物家族中，该家族在同一晶格[66]上实现了有效（[63]）、（[64]）、（，）和（）[65]。我们强调，这些化合物目前均未处于本文所研究的转变点上：氟化物家族的Ni和Mn成员近乎海森堡型（对于，极化中子散射将单离子各向异性限制为交换耦合的一小部分[67]），铁成员的单离子各向异性尚未通过实验确定，而具有主导易轴交换的材料通常深藏于某一相内部。相反，本模型为这一日益增长的材料图景提供了组织相图——典型的自旋冰和占据极端易轴极限，其中晶体场分离出最大的双重态——并确定了哪些可通过化学取代、应变或压力调控的单离子参数会驱动材料跨越交叉脊线——对于，则跨越本文所分类的一阶线。B. 高压冰在磁性之外，方程(1)的成员描述了高压水冰的氢键网络，其中代表氢键上的两个偏心质子位置，代表键中心位置，代表压力调谐的居中能量，而离子缺陷（）扮演磁单极的角色。在这个语言中，质子无序的冰VII对应于大时的受限库仑相，键对称化的冰X对应于小时的平庸顺磁体，而第六节推导的磁单极诱导的圆化解释了为何冰VII–X转变被观察为宽泛的、依赖于探针的交叉而非锐利转变；一篇配套论文详细发展了这一应用[23]。值得注意的是，竞争在那里直接由质子坐标实现，无需微调，而同位素效应（相对于）使交叉压力移动约10 GPa[68]，使模型参数在实验上可扫描。C. 工程化平台Blume-Capel（）自由度最近已在具有光刻可调各向异性的人工纳米磁体阵列中实现[24]，可编程平台——实现人工自旋冰的超导磁通量子退火器[69]和提出的烧绿石晶格上的三维里德堡原子阵列[70]——提供了原则上可以独立控制两个逸度和的环境。表I和II中的定量相边界以及预测的零场临界终点为这类平台提供了具体目标。D. 可测量的特征本工作的三个预测可直接检验。首先，对于实现，一阶线应严格在零外场处终止于临界终点，在此处比热和磁化率在孤立温度处表现出3D伊辛临界标度（第六节）。其次，比热的交叉脊线沿直线延伸，具有表I和II中依赖的斜率（对于和，相应的实验温度尺度已在文献[16,19]中强调）。第三，在处的肖特基异常为在候选高自旋烧绿石磁体中分离单离子特征与集体交叉现象提供了实用方案（第七节）。IX. 结论与展望在本工作中，我们发展了一个自洽的理论框架，对任意自旋的经典烧绿石磁体的宏观拓扑相、涌现规范理论和临界现象进行了分类。该框架结合了热力学能量凸性、微观晶格几何（）、精确图形映射和重整化群分析。在无磁单极极限下，我们确立了相图的完备性：中间自旋幅度的均匀背景被单离子能量的凸性排除，而部分退禁闭情景通过竞争缺陷弦的自由能比较被排除。在小极限下出现自旋奇偶二分性：半整数自旋约化为无散度库仑流体，而整数自旋映射到三维模型。在大极限下，我们从几何上证明了宏观极化通量被量子化为的倍数，确立了涌现的受限库仑相。该极限下的分类总结在表I和II中。退禁闭转变的性质揭示了进一步的二分性。虽然缺陷映射到三维伊辛类，但所有缺陷都产生手性、有向环。冰规则相容性准则表明，高斯定律与顶点层面的物理冰规则相容当且仅当。利用这一几何性质，我们将缺陷气体映射到三态波茨模型，其对称性允许的三次不变量驱动一阶转变。对于，这种离散钟映射失效：模守恒律不再等价于物理冰规则。系统反而被迫进入具有指数抑制中间缺陷的分级融合过程。利用配分函数的分解，我们证明宏观图集在领头阶约化为连续U(1)三维模型。残余融合级联在三维不动点处作为危险无关算符起作用；其指数小的裸权重使系统深陷该不动点的吸引域内，其无关性固定了普适类并保证渐近临界行为的早期出现，将所有退禁闭转变置于三维普适类中。最后，通过放宽严格冰规则，我们证明热激发磁单极在涌现规范理论中充当显式对称破缺场。该场切断拓扑缺陷弦并使所有拓扑扇区绝热连通。其后果依赖于序：连续三维伊辛（）和三维（）转变在任意时被圆化为交叉，而一阶转变被预测在有限温度下存活，由宏观自由能垒保护，并终止于临界终点。对整数（）和半整数（）自旋的经典蒙特卡洛模拟在可达温度下证实了这一交叉图像：比热峰随系统尺寸增大而饱和至有限值，交叉尺度与解析预测的逸度定量一致。然而，我们强调本工作数值部分的一个重要局限：可达温度高于估计的共存线临界终点（在时为，而Bethe-Peierls估计为），因此所预测的一阶转变的存活——以及更一般地，无磁单极转变的普适类——并未被当前模拟直接验证。直接检验需要使用第七节所讨论的环算法对流形进行采样，这是未来工作的一个重要方向。在较小的下的初步并行回火模拟已进入估计终点以下，显示出与预测存活一致的迹象（第七节），但尚未最终确证。自旋幅度占据进一步识别出处的非合作肖特基异常，清晰地将局域单离子物理与集体拓扑现象分离。这种清晰的解耦由解析预测的肖特基峰位置的标度所控制，为实验学家在候选高自旋烧绿石磁体中从集体拓扑现象的特征中分离宽比热异常提供了实用诊断。将此框架扩展到量子区是未来工作的自然方向。量子涨落、涌现离散规范场与所特有的几何性质之间的相互作用可能产生新的分数量子自旋液体和高自旋拓扑相[2,30,71–74]。更广泛地说，本文所识别的机制——微观实值约束使模型保持U(1)对称，而晶格几何指数抑制了会挑出离散子群的过程——可能适用于其他受挫格点模型，在这些模型中离散规范对称性与连续规范对称性竞争，包括非二分格点上的二聚体模型和更高秩张量规范理论。