Bootstrapping ground state properties of classical frustrated magnets

想象一下，你正在试图找到一种完美的方式，将一大群人安排在一个巨大且无限的体育场中。每个人都有一个特定的规则：他们必须站在自己个人空间中心恰好一米远的地方（就像单位长度的自旋）。然而，他们也有相互冲突的意愿：有些人希望面向邻居，而另一些人则希望背对邻居。这是一个“受挫”系统，因为你无法同时满足所有人的意愿。

目标是找到一种能让人群尽可能平静（能量最低）的排列方式。这是物理学中的一个经典问题，但极其难以解决，因为人数众多且规则冲突，导致数学变得混乱且充满“死胡同”。

以下是作者 Nisarga Paul 和 Gil Refael 使用一种名为**自举法（bootstrapping）**的新方法解决该问题的过程。

问题：一个充满死胡同的迷宫

将传统的解决方法想象成试图在巨大且雾气弥漫的山脉中找到最低点。你可能会开始沿着山坡向下走，但很容易被困在一个小山谷（局部极小值）中，以为那就是底部，而实际上附近有一个更深得多的山谷。

旧方法（Luttinger-Tisza）： 这就像从非常遥远且模糊的距离观察山脉。对于简单的山脉，它能给出一个不错的猜测，但如果地形怪异或规则复杂，这个猜测往往是错误的。
模拟方法（蒙特卡洛）： 这就像派遣一个机器人在山脉中行走。但在受挫系统中，机器人会感到困惑，原地打转，永远找不到真正的底部。

解决方案：“影子”方法（自举法）

作者们决定不去寻找每个人的确切位置（这实际上是不可能的），而是观察人群投射出的影子。

想象一下，你不知道人们站在哪里，但你知道游戏的规则：

正定性： 如果你问“两个人以某种特定方式站立的概率是多少？”，答案不能是负数。
归一化： 每个人都必须存在（总概率为 1）。
几何约束： 人们站在一个球面上（他们不能拉伸或收缩）。

作者们创建了一个数学“筛子”或一系列过滤器。他们从一个非常宽松的过滤器开始，仅检查基本规则。然后，他们添加越来越多复杂的过滤器，以检查人们之间更深层的关系。

类比： 想象试图通过观察一个隐藏物体的影子来猜测其形状。
- 第一层： 你看到一个看起来像圆形的影子。该物体可能是一个球、一个盘子或一枚硬币。
- 第二层： 你添加了第二个光源。现在影子必须同时匹配这两个角度。该物体现在被缩小为球或盘子。
- 第三层： 你添加了第三个光源。现在影子必须匹配三个角度。该物体肯定是一个球。

在这篇论文中，“影子”是关联函数（一个自旋如何与另一个自旋相关）。“光源”是称为**半定规划（SDP）**的数学约束。

实际运作方式

作者们构建了这些过滤器的层级结构：

设置： 他们定义了无限体育场的一小块区域（几排座位）。
约束： 他们强制数学在该区域内遵守概率和几何规则。
结果： 计算机求解一个“凸优化”问题。这是一种没有死胡同的数学问题；它总是能在该特定过滤器的规则范围内找到最佳可能答案。

随着他们扩大区域并添加更复杂的过滤器（层级的更高层级），“影子”变得越来越清晰。

下界： 该方法给出了人群能达到的平静程度的保证“底线”。它表示：“能量不可能低于 X。”
上界： 他们还使用标准模拟找到一种特定排列并计算其能量，从而给出一个“上限”。“能量不可能高于 Y。”

结果的神奇之处

在许多情况下，“底线”和“上限”几乎完美地重合。

精度： 他们以极高的精度找到了基态的确切能量（在某些情况下精确到小数点后 8 位）。
无需猜测： 与其他方法不同，这种方法不依赖于猜测起始点。它提供了严格的证明，表明答案位于一个极小的范围内。
速度： 尽管数学很复杂，但计算机可以在每个设置下仅需几秒钟就能解决这些问题。
可视化人群： 一旦他们找到了“影子”，就可以逆向工程出实际的人群排列（自旋纹理）是什么样子的。它与其他方法得出的最佳猜测完全吻合。

为什么这很重要

这种方法就像在一个一切都很模糊的世界中拥有一把超精准的尺子。

它适用于任何形状的体育场（不仅仅是简单的网格）。
它适用于任何类型的规则（甚至是复杂的非线性规则）。
它在无限极限下有效（理论上完美），而不仅仅是在小型计算机模拟中。

作者们表明，通过观察“影子”（关联）并收紧规则（层级），他们可以解决一个以前被认为太难以确定性解决的问题。他们不仅仅是猜测答案；他们从数学上证明了答案必须存在的范围。

技术摘要：经典阻挫磁体基态性质的自举法

问题陈述
确定经典阻挫磁系统的基态自旋构型与能量密度是一个极其艰巨的优化问题。对于伊辛自旋，该问题已知是 NP 难的。对于具有平移不变哈密顿量的海森堡自旋，标准的解析方法是 Luttinger–Tisza (LT) 方法，该方法通过松弛局部自旋归一化约束，在螺旋和（sums-of-spirals）假设上最小化能量。然而，LT 方法的适用范围有限，特别是在非布拉维晶格和非二次哈密顿量的情况下。数值替代方案（如经典蒙特卡洛）在阻挫系统中常因能量景观中存在众多紧密竞争的局部极小值而面临严重的平衡化挑战。此外，变分法和蒙特卡洛方法通常在有限尺寸系统上运行，缺乏关于其接近真实热力学极限基态的严格保证。

方法论
作者引入了一种基于半定规划（SDP）的自举方法，用于为无限晶格上平移不变的经典自旋模型产生基态能量密度和相关函数的严格双边界限。该方法将多项式优化中的 Lasserre 层级结构适配到阻挫磁学的语境中。

多项式优化表述：寻找基态的问题被简化为一个多项式优化问题，其约束条件为自旋的单位范数条件（ $|S_i|^2 = 1$ ）。
矩层级：该方法不直接对自旋构型进行优化（这是一个非凸问题），而是对平移不变的相关函数（矩）空间进行优化。该方法定义了一个由实空间补丁 $P$ 和正整数 $r$ （层级水平）标记的有限尺寸凸优化层级。
变量与约束：
- 变量：优化变量是补丁 $P$ 内次数高达 $r$ 的自旋单项式 $m_I$ 的期望值 $y[I] = \langle m_I \rangle$ 。
- 正定性（矩矩阵）：构建一个由单项式对索引的矩矩阵 $L$ 。 $L$ 为半正定（ $L \succeq 0$ ）的条件确保了这些矩对应于一个有效的概率分布。
- 归一化（局域化矩阵）：单位长度约束通过局域化矩阵 $C_r$ 强制执行，这些矩阵施加了由恒等式 $\sum_a (S_r^a)^2 - 1 = 0$ 导出的矩之间的线性关系。
收敛性：作者证明，随着补丁大小 $P$ 和层级水平 $r$ 趋于无穷大，能量密度的下界单调收敛至真实的热力学极限基态能量密度。该证明利用了测度论（Kolmogorov 扩张定理）和实代数几何的工具。
上界与自旋织构：为了获得双边界限，该方法与通过在大周期团簇上进行迭代单点最小化获得的变分上界（ $\epsilon_{ub}$ ）相结合。此外，可以通过对两点相关矩阵进行特征值分解，将其投影到前三个主特征向量上，随后进行归一化，从而从 SDP 解中提取显式的自旋织构。

主要贡献

严格的双边界限：该方法直接在热力学极限下提供基态能量密度和任意多项式可观测量的严格上下界，消除了对有限尺寸外推的需求。
Luttinger–Tisza 方法的推广：该方法将 LT 方法作为特例（具体而言是在布拉维晶格的最低层级水平上）纳入其中，但将其适用性扩展到了非布拉维晶格和非二次（例如四次）哈密顿量。
效率：该方法计算高效，使用标准 SDP 求解器（如 MOSEK）时，每个参数点的典型运行时间为数秒。
织构提取：该方法允许提取显式的候选基态自旋织构，即使在高度阻挫的机制下，这些织构也能与变分结果很好地吻合。

结果
该方法被应用于两个二维阻挫自旋模型：

中心方格晶格（二次海森堡模型）：这是一个具有不等价位点和非均匀配位数的模型，LT 方法在此失效。SDP 下界和变分上界在整个参数范围内几乎重合，证明了 SDP 优于 LT 方法，并具备认证变分结果的能力。
三角晶格（双线性 - 双二次模型）：这是一个高度阻挫的四次模型，已知在中间耦合下具有广泛的基态简并。SDP 在整个相图中提供了紧密的能量界限（每个位点的括号范围为 $10^{-5}$ 到 $10^{-8}$ ）。在高度阻挫区域，界限较宽，反映了固有的困难，但 SDP 成功认证了变分结果。提取的自旋织构在有序相中与变分构型匹配，并在阻挫相中从简并流形中选择了代表。

意义与主张
本文主张，这种自举方法可作为经典阻挫磁学中基态问题的实用、通用工具，是对现有启发式方法和蒙特卡洛模拟的补充。其主要意义在于在热力学极限下提供严格保证，这是变分法和蒙特卡洛方法所缺乏的特征。作者强调，该方法解决了先前解析方法（如 LT 方法）在处理非二次相互作用和非布拉维晶格方面的局限性。

作者谦逊地将这项工作定位为该方法及其在特定二维模型中应用的介绍。他们指出了未来的方向，例如放松平移不变性以研究自旋玻璃、开发用于大自旋量子修正的半经典自举法，以及将该方法应用于经典动力学过程，但并未声称这些已在当前工作中实现。