Bootstrapping ground state properties of classical frustrated magnets

本論文は、非二次のハミルトニアンと非ブラベー格子を高い効率で扱うことで、従来の解析手法の限界を克服し、並進対称性を持つ古典的フラストレーション磁石の基底状態エネルギー密度および相関関数に対する収束する両側境界を生成するためにラッセル階層を適応させた厳密な半正定値計画法を導入する。

要約

想像してください。巨大で無限のスタジアムに大勢の人々を配置する、完璧な方法を見つけようとしている状況を。各人には特定のルールがあります。彼らは、自分自身のパーソナルスペースの中心から正確に 1 メートルの距離に立たなければなりません（単位長さのスピンのようなもの）。しかし、彼らには矛盾する欲求もあります。隣人の方を向きたい人もいれば、背を向けたい人もいます。これは「フラストレーション（欲求不満）」を伴う系です。なぜなら、すべての人の欲求を同時に満たすことは不可能だからです。

目標は、大勢の人々をできるだけ静か（低エネルギー）にする配置を見つけることです。これは物理学における古典的な問題ですが、人々があまりにも多く、矛盾するルールが多すぎるため、数学が複雑になり、「行き止まり」に満ちてしまうため、解決は極めて困難です。

ここで、ニサールガ・ポールとギル・レファエルという著者たちが、ブートストラッピングと呼ばれる新しい手法を用いてこの問題を解決した方法をご紹介します。

問題：行き止まりの多い迷路

この問題を従来の方法で解くことは、霧のかかった巨大な山脈で最も低い地点を見つけようとするようなものです。丘を下り始めても、小さな谷（局所最小値）に迷い込み、それが底だと思い込んでしまうかもしれませんが、実際にはその近くにずっと深い谷が存在する可能性があります。

• 従来の方法（ロッティガー・ティッサ法）： これは、山を非常に高い位置からぼんやりと眺めるようなものです。単純な山には良い推測を与えますが、地形が奇妙だったりルールが複雑だったりする場合、その推測はしばしば誤りでした。
• シミュレーション手法（モンテカルロ法）： これは、ロボットを山の中を歩かせて探すようなものです。しかし、フラストレーションを伴う系では、ロボットは混乱し、くるくると回り続け、真の底を見つけることができません。

解決策：「影」を用いた手法（ブートストラッピング）

著者たちは、一人ひとりの正確な配置を見つけること（それは不可能です）の代わりに、大勢の人々が投げる影を見ることにしました。

人々がどこに立っているかはわからないが、ゲームのルールはわかっていると想像してください。

1. 正性： 「二人が特定の状態で立っている可能性はどれくらいか？」と問えば、その答えは負にはなり得ません。
2. 規格化： 全員が存在しなければならないため（総確率は 1 です）。
3. 幾何学： 人々は球面上に立っています（伸び縮みはできません）。

著者たちは、数学的な「篩」あるいは一連のフィルターを作成しました。まず、基本的なルールのみをチェックする非常に緩いフィルターから始め、その後、人々の間のより深い関係性をチェックする、より複雑なフィルターを次々と追加していきました。

• 比喩： 隠された物体の形を、その影を見て推測しようとする状況を想像してください。
  • レベル 1： 円のように見える影が見えます。物体はボール、皿、または硬貨の可能性があります。
  • レベル 2： 2 番目の光源を追加します。これで影は両方の角度に一致しなければなりません。物体はボールか皿に絞り込まれます。
  • レベル 3： 3 番目の光源を追加します。これで影は 3 つの角度に一致しなければなりません。物体は間違いなくボールです。

この論文において、「影」とは相関関数（あるスピンが他のスピンとどのように関連しているか）です。「光源」とは、**半正定値計画法（SDP）**と呼ばれる数学的制約です。

実際の実行方法

著者たちは、これらのフィルターの階層を構築しました。

1. 設定： 無限のスタジアムの小さな断片（いくつかの列の座席）を定義しました。
2. 制約： その断片内において、確率と幾何学のルールに従うよう数学に強制しました。
3. 結果： コンピュータは「凸最適化」問題を解きます。これは行き止まりのない数学の問題であり、その特定のフィルターのルール内では常に最良の答えを見つけます。

彼らが断片を大きくし、より複雑なフィルター（階層のより高いレベル）を追加するにつれて、「影」は次第に鮮明になりました。

• 下限： この手法は、大勢の人々がどれだけ静かになり得るかの保証された「床」を提供します。「エネルギーは X より低くなることはあり得ない」と言います。
• 上限： 彼らはまた、標準的なシミュレーションを用いて特定の配置を見つけ、そのエネルギーを計算することで、「天井」を与えました。「エネルギーは Y より高くなることはあり得ない」と。

結果の魔法

多くの場合、「床」と「天井」はほぼ完璧に一致しました。

• 精度： 彼らは基底状態の正確なエネルギーを、驚異的な精度（場合によっては小数点以下 8 桁の精度）で見つけ出しました。
• 推測なし： 他の手法とは異なり、これは出発点を推測することに依存しません。答えが微小な範囲内にあることを厳密に証明します。
• 速度： 数学は複雑ですが、コンピュータは設定ごとにわずか数秒でこれらの問題を解くことができました。
• 大勢の可視化： 「影」を見つけると、それを逆算して、実際の配置（スピンテクスチャ）がどのように見えるかを把握できました。それは他の手法からの最良の推測と完全に一致しました。

なぜこれが重要なのか

この手法は、すべてがぼやけている世界に対する超精密な定規のようなものです。

• スタジアムの任意の形状（単純な格子だけでなく）で機能します。
• 任意の種類のルール（複雑な非線形なものさえも）で機能します。
• 小さなコンピュータシミュレーションだけでなく、無限の極限（理論的に完璧）でも機能します。

著者たちは、「影」（相関）を見てルール（階層）を厳しくすることで、以前は確実性を持って解くことが難しすぎると考えられていた問題を解決できることを示しました。彼らは単に答えを推測したのではなく、答えが存在しなければならない範囲を数学的に証明しました。

技術概要

技術的概要：古典的フラストレーション磁性体の基底状態特性のブートストラップ

問題提起
古典的フラストレーション磁性系の基底状態スピン配置とエネルギー密度を決定することは、困難な最適化問題である。イジングスピンに対しては、この問題は NP 困難であることが知られている。転移不変なハミルトニアンを持つハイゼンベルグスピンに対しては、標準的な解析的手法として、局所スピン正規化制約を緩和することでスピラル和のアナスタス上でのエネルギー最小化を行うルッテガー・ティザ（LT）法がある。しかし、LT 法は、特に非ブラベー格子や非二次ハミルトニアンに対しては有効性が限られている。古典的モンテカルロ法などの数値的代替手法は、多くの密接に競合する局所最小値を有するエネルギー地形のため、フラストレーション系において深刻な平衡化の課題に直面することが多い。さらに、変分法やモンテカルロ法は通常有限サイズ系で動作し、真の熱力学的極限における基底状態への近接性に関する厳密な保証を欠いている。

手法
著者らは、無限格子上の転移不変な古典的スピンモデルの基底状態エネルギー密度および相関関数に対する厳密な両側境界を生成する、半正定値計画法（SDP）に基づくブートストラップ法を導入する。このアプローチは、多項式最適化からのラッセル階層をフラストレーション磁性の文脈に適応させたものである。

1. 多項式最適化としての定式化: 基底状態を見つける問題は、スピンが単位ノルムを持つという制約（$|S_i|^2 = 1$）を伴う多項式最適化問題に帰着される。
2. モーメント階層: スピン配置を直接最適化する（非凸問題）のではなく、転移不変な相関関数（モーメント）の空間上で最適化する。この手法は、実空間パッチ $P$ と正整数 $r$（階層レベル）によってラベル付けされた有限サイズの凸最適化の階層を定義する。
3. 変数と制約:
   • 変数: 最適化変数は、パッチ $P$ 内の次数 $r$ までのスピン単項式 $m_I$ の期待値 $y[I] = \langle m_I \rangle$ である。
   • 正性（モーメント行列）: 単項式のペアでインデックス付けされたモーメント行列 $L$ を構成する。$L$ が半正定値（$L \succeq 0$）であるという条件は、モーメントが有効な確率分布に対応することを保証する。
   • 正規化（局在化行列）: 単位長さの制約は、$\sum_a (S_r^a)^2 - 1 = 0$ という恒等式から導かれるモーメント間の線形関係を課す局在化行列 $C_r$ によって強制される。
4. 収束: 著者らは、パッチサイズ $P$ と階層レベル $r$ が無限大に成長するにつれて、エネルギー密度の下限が単調に真の熱力学的極限における基底状態エネルギー密度に収束することを証明する。この証明には、測度論（コルモゴロフの拡張定理）と実代数幾何学の道具が用いられる。
5. 上限とスピンテクスチャ: 両側境界を得るために、この手法は、大きな周期的クラスター上の反復的な一サイト最小化によって得られる変分上限（$\epsilon_{ub}$）と組み合わされる。さらに、明示的なスピンテクスチャは、二点相関行列の固有値分解を行い、3 つの主要な固有ベクトルへ射影した後に正規化することによって、SDP 解から抽出できる。

主要な貢献

• 厳密な両側境界: この手法は、熱力学的極限において、有限サイズ外挿を必要とせずに、基底状態エネルギー密度および任意の多項式観測量に対する厳密な下限と上限を提供する。
• ルッテガー・ティザ法の一般化: このアプローチは、特にブラベー格子における最下位階層レベルにおいて LT 法を特殊ケースとして包含するが、非ブラベー格子や非二次（例えば、四乗）ハミルトニアンへの適用性を拡張する。
• 効率性: この手法は計算効率が良く、標準的な SDP ソルバー（例：MOSEK）を用いて、パラメータ点あたり通常数秒で実行される。
• テクスチャ抽出: この手法は、高度にフラストレーションされた領域であっても、変分結果とよく一致する明示的な候補基底状態スピンテクスチャの抽出を可能にする。

結果
この手法は、2 つの2次元フラストレーションスピンモデルに適用された。

1. 中心付き正方形格子（二次ハイゼンベルグモデル）: 不等価なサイトと非一様な配位を有し、LT 法が失敗するモデル。SDP 下限と変分上限はパラメータ範囲全体でほぼ一致し、LT 法に対する SDP の優位性と、変分結果を証明する能力を実証した。
2. 三角形格子（二項 - 二項モデル）: 中間結合において広範な基底状態縮退で知られる、高度にフラストレーションされた四乗モデル。SDP は、相図全体にわたり、サイトあたり $10^{-5}$ から $10^{-8}$ の狭いエネルギー境界（括弧）を提供した。高度にフラストレーションされた領域では、本質的な難しさを反映して境界は広くなったが、SDP は変分結果を成功裏に証明した。抽出されたスピンテクスチャは、秩序相において変分配置と一致し、フラストレーション相において縮退多様体から代表者を選択した。

意義と主張
本論文は、このブートストラップアプローチが、既存のヒューリスティックなアプローチやモンテカルロシミュレーションを補完する、古典的フラストレーション磁性における基底状態問題のための実用的で汎用的なツールであると主張する。その主な意義は、変分法やモンテカルロ法には欠けている、熱力学的極限における厳密な保証を提供する点にある。著者らは、この手法が非二次相互作用や非ブラベー格子に関する先行する解析的手法（LT 法等）の限界に対処することを強調している。

著者らは、この作業を手法の導入と、特定の2次元モデルへの適用として謙虚に位置づけている。彼らは、スピンガラスを研究するための転移不変性の緩和、大スピン量子補正のための半古典的ブートストラップの開発、および古典的動的プロセスへのアプローチの適用などの将来の方向性を特定しているが、これらが現在の作業で実現されたとは主張していない。