Finite-temperature formation of magnetic plateaus and simplex liquid states on the frustrated ruby lattice

信念伝播法を用いて最適化された無限テンソルネットワーク状態を用いた本研究は、ルビー格子上のフラストレートしたスピン1/2ハイゼンベルク反強磁性体が、低温において残留エントロピーを伴う新規な「simplex liquid state（シンプレックス・リキッド状態）」を宿す安定な磁気プラトーを形成しており、この過程は系が冷却されるにつれて相転移を伴わずに連続的に起こることを明らかにしている。

要約

混み合ったダンスフロアを想像してみてください。そこでは誰もがパートナーを探そうとしていますが、ダンスのルールが非常にトリッキーです。これは、「ルビー格子（ruby lattice）」と呼ばれる特別な種類の磁性材料に関する新しい論文の物語です。

科学者たちが発見した内容を、簡単な比喩を用いて解説します：

トリッキーなダンスフロア（幾何学的フラストレーション）

通常の磁石では、原子の微小なスピンは、小さなコンパスの針のように、すべて同じ方向に並ぼうとします（行進する兵士のように）。しかし、この特定の「ルビー格子」構造では、その幾材学的な形状があまりにねじれているため、すべてのスピンが同時に「幸せ（安定）」になることができません。これは、椅子取りゲームのようなものです。椅子は人数よりも多いのですが、椅子の配置のせいで、誰かが誰かにぶつかることなく、全員が快適に座ることが不可能な状態になっています。これは幾何学的フラストレーションと呼ばれます。

通常、これらの材料を冷却すると、フラストレーションを感じ、最終的には秩序ある硬いパターン（結晶のように）へと「パチン」と収まります。しかし、科学者たちは、もし非常にゆっくりと、かつ慎重に冷却したらどうなるのかを知りたいと考えました。

「信念伝播（Belief Propagation）」の魔法

これを解明するために、研究者たちはテンソルネットワーク、特に**信念伝播（BP）**と呼ばれる強力なコンピュータ手法を用いました。

信念伝播を、大勢の群衆の中に広がる「噂」のようなものだと考えてください。群衆の一人ひとりに何をしているのかを聞き回る代わりに、数人の人に聞き、その人たちが隣の人に伝え、さらにその隣の人へと伝えていくのです。最終的に、全員が一人ひとりに直接確認することなく、グループ全体がどのような状態にあるのかを理解できるようになります。研究者たちは、この「噂を広める」数学的手法を用いて、システムが無限に大きい場合でも、異なる温度における磁気スピンの振る舞いをシミュレートしました。

驚きの結果： 「スナップ」ではなく、「液体」

システムを冷却したとき、彼らはスピンが突然、硬い秩序ある結晶へと「スナップ」して固まる（相転移する）ことを予想していました。しかし、実際に発見されたのは、もっと流動的なものでした。

温度が下がっても、スピンは単一のパターンへと凍りつくことはありませんでした。代わりに、彼らは**「シンプレックス・リキッド状態（Simplex Liquid State）」**を形成しました。

• 比喩： パーティー会場にいる人々の集まりを想像してください。全員が完璧な格子状に立っている（結晶）のではなく、彼らは3人組の小さな結束力の強いグループ（シンプレックス）を作ります。これらのグループは一緒に踊っていますが、その「グループの配置」自体は常に変化し続けています。
• 結果： 極低温においても、システムは無秩序な状態を保ちます。これは、踊るグループの「液体」なのです。これらのグループには多くの異なる配置の仕方があるため、システムは多くの「残留エントロピー（無秩序さの尺度）」を保持しています。それは、見るたびに完璧にシャッフルされるトランプのデッキのように、決して特定の順序に落ち着くことがない状態です。

磁気プラトー（平坦部）

研究者たちは、磁場（ダンスフロアに吹き付ける強い風のようなもの）もオンにしました。磁場を強くしていくにつれ、スピンはそれに沿って整列しようとします。

スムーズに回転する代わりに、スピンは特定の「プラトー（平坦部）」で止まってしまいます。

• 比喩： 階段を想像してください。システムをより強く押し進めると、磁化（整列の度合い）は跳ね上がった後、しばらく平坦（プラトー）になり、再び跳ね上がります。
• 彼らは、磁化が最大値の1/3、1/2、または2/3となる安定した「ステップ」を発見しました。
• ひねり： これらの平坦な「ステップ」の上であっても、材料は硬い結晶にはなりませんでした。特定の平均的な整列を持ちながらも、依然としてその「液体」状態に留まっていたのです。

「ラムダ」ピークとスイッチ

階段の中ほど（1/2プラトー付近）に、非常に興味深い瞬間がありました。

• 比喩： ダンスフロアが二つに分かれていると想像してください。片側ではグループがある方向に踊り、もう片側では別の方向に踊っています。特定の温度と風速において、フロア全体が一方のダンススタイルからもう一方のスタイルへと突如として切り替わります。
• この切り替えはスムーズなものではありませんでした。それは、ギリシャ文字のラムダ（λ）のような形をした、鋭い「熱容量（システムが吸収するエネルギーの量）」のスパイク（急上昇）を生み出しました。これは、プラトーの極めて端の部分において、システムが量子ゆらぎによって駆動される大きな変化の瀬戸際にあることを示唆していますました。

大きな教訓

最も重要な発見は、このシステムは決して伝統的な結晶へと「凍結」することはないということです。

極限まで冷やしても、システムは無秩序で液体のような状態を保ち、多くの可能性に満ちています。研究者たちは、「熱容量（システムが温度変化にどのように反応するかを示す尺度）」が滑らかで連続していることを示すことで、これを証明しました。もしシステムが結晶へと凍りついていたならば、相転移を示す鋭くギザギザしたスパイクが現れたはずです。しかし、代わりに、システムは未知の、エキゾチックな状態へと滑らかに流れ込んでいったのです。

要約すると： 研究者たちは、巧妙な「噂を広める」数学的トリックを用いることで、フラストレーションを受けた磁性材料は、冷却されても硬い結晶には凍りつかないことを証明しました。代わりに、それは踊るスピン・グループの「液体」へと姿を変え、温度が絶対零点に近い状態であっても、無秩序で可能性に満ちた状態であり続けるのです。

技術概要

技術要約：ルビー格子における磁性プラトーおよびシンプレックス・リキッド状態の有限温度形成

問題提起
幾何学的フラストレーションを持つ量子系は、低温において伝統的な磁気秩序を回避する量子スピン液体のような、非従来型の相を安定化させることが知られている。しかし、有限温度におけるこれらの状態の形成とその熱力学的シグネチャを理解することは、依然として大きな課題である。具体的には、次近接相互作用を持つルビー格子上のスピン1/2ハイゼンベルク反強磁性体のゼロ温度（$T=0$）基底状態は、以前の研究において分数化された磁気プラトー（$m_z = 0, 1/3, 1/2, 2/3$）を示すことが示されていた。これらのプラトーは、格子対称性を破る離散的な結晶相として解釈されてきた。未解決の課題として、正確な基底状態の縮退、有限温度におけるプラトーの安定性、およびこれらに到達するための過程が従来の相転移を経るのかどうかという点が残されていた。これらの問題を解決するには、熱的揺らぎと量子揺らぎが競合する実験的実現と計算結果との架け橋となるシミュレーションが必要である。

手法
著者らは、平衡密度行列の平方根の非正規化表現 $\rho^{1/2}(\beta) = \exp(-\beta H/2)$ （ここで $\beta = 1/T$）を記述するために、精製された無限テンソルネットワーク状態（iTNS）アンザッツを採用している。本研究は、ルビー格子上のハイゼンベルク・ハミルトニアンの等方的なケース（$J_t = J_h = J_d = J$）に焦点を当てている。

• 格子表現: 三角単位の頂点は三角形（シンプレックス）にグループ化され、ルビー格子を配位数 $z=3$ の粗視化されたハニカム幾何学へと写像している。iTNSは2サイトの単位胞を使用しており、各テンソルは6つの物理インデックス（各三角形内の3つのスピンを表す）と3つの仮想インデックスを持つ。
• アルゴリズム: システムは無限温度（$\rho^{1/2}(0) = I$）から初期化され、プロパゲーターの2次トロッター化を用いた虚時間発展によって冷却される。
• 観測量: 熱力学的量は2つのアプローチを用いて算出される：
  1. 信念伝播法 (BP): ノルムネットワークに対して標準的なBP縮退を適用し、期待値を計算する。
  2. ループ補正BP: 平均場近似的なBPを超えて精度を高めるため、著者らは自由エネルギー密度に対して第1次クラスター補正（ループ補正）を適用している。これにより、格子上の六角形周囲の相関を考慮に入れている。
• スケール: 著者らはGPUハードウェアを利用して、ボンド次元 $\chi \sim 200$ まで到達しており、これにより広い温度範囲にわたる熱力学的観測量の熱力学的極限における正確な測定を可能にしている。

主要な結果

1. ゼロ磁場におけるシンプレックス・リキッド状態:

   • 無限温度から冷却する過程で、システムは結晶秩序への相転移を起こさない。代わりに、非ゼロの残留エントロピー（$s_0 \approx \frac{1}{3}\lambda$、ここで $\lambda$ はハニカム格子上のダイマー被覆の残留エントロピー）を特徴とする、無秩序な「シンプレックス・リキッド状態」を形成する。
   • 熱容量（$C_v$）はすべての温度で有限かつ連続であり、相転移を示す発散するピークは見られない。
   • 簡約密度行列の解析により、リーディング固有値 $\lambda_1 \approx 1/3$ が得られた。これは、2つのシンプレックスが強くペアリングされている状態に対応する。残りの固有値は、弱く相関した状態に対応している。システムは、ペアリングされたスピン・シンプレックスを含む指数関数的に多数の結晶構成 $|\psi_d\rangle$ の完全な無秩序混合物である。
   • このリキッド状態のエネルギー密度（$e/J \approx -0.498$）は、以前の $T=0$ 研究で見出された特定の結晶構成の基底状態エネルギーと一致するが、システムは単一の結晶を選択するのではなく、並進対称性を保持したリキッド状態に留まっている。

2. 磁気プラトーと磁場誘起挙動:

   • 外部磁場の存在下で、システムは磁化密度 $m_z = 0, 1/3, 1/2, 2/3$ において安定な磁気プラトーを示す。
   • 熱力学的シグネチャ: これらのプラトーの形成は、従来の相転移を経ることなく滑らかに進行する。熱容量は連続的であるが、浅いピークが観察される。
   • 残留エントロピー: プラトーは広範な縮退を持つ「シンプレックス・リキッド状態」を宿している。特に $m_z = 1/2$ プラトーは、$m_z = 0$ プラトーの約2倍の大きな残留エントロピーを示す。
   • 安定性: $m_z = 0$ プラトーは熱的揺らぎに対して最も安定している。$m_z = 1/3$ プラトーは安定性が低く、非常に低い温度（$T/J \le 10^{-2}$）でのみ狭い磁場範囲に現れる。$m_z = 2/3$ プラトーは基底状態計算では観察されているが、現在の有限温度シミュレーションで完全に解像するには、より低い温度が必要である。
   • フェーズ・スイッチング: $m_z = 1/2$ と $m_z = 1/3$ のプラトーの界面において、システムは熱的に駆動されたリキッド状態間のスイッチングを示す。これは、熱容量における非対称な「ラムダ」ピークと、磁化率の急激な上昇によって特徴付けられ、おそらくゼロ温度量子相転移付近の臨界揺らぎを反映している。

意義と主張
本論文は、信念伝播法とループ補正によって最適化された無限テンソルネットワーク状態、特にテンソルネットワーク技術が、有限温度におけるフラストレートした量子磁性体を理解するための強力かつスケーラブルな経路を提供することを実証している。

• 手法の検証: 本研究は、BPベースのテンソルネットワーク・アプローチが、エネルギーギャップや残留エントロピーの検出を含め、多くのオーダーの温度範囲にわたって熱力学的特性を正確に捉えることができることを検証している。
• 基底状態の性質: 本研究は、ルビー格子上の磁気プラトーが単純な結晶秩序ではなく、「シンプレックス・リキッド状態」、すなわち非ゼロの残留エントロピーを保持し、格子対称性を破らない指数関数的に多数の結晶構成の無秩序混合物であることを明らかにしている。
• 相転移の不在: 本研究の重要な発見は、システムが有限温度の相転移を経ることなく、これらの安定な磁性相に到達することであり、これは連続的かつ有限な熱容量によって裏付けられている。
• 今後の方向性: 著者らは、本手法を3次元のフラストレートした系（パイロクロア格子など）の研究や、基底状態の縮退をトポロジカルに秩序化した重ね合わせ状態へと破る摂動の研究へと拡張できることを示唆している。