Thermodynamic geometry of three coupled Ising spins: exact curvature, frustration and bond anisotropy

技術要約：3つの結合したイジングスピンの熱力学的幾何学

問題提起
熱力学的曲率、特にフィッシャー・ラーオ（ルピナー）曲率 $R$ は、相互作用、相転移、および臨界スケーリングを特定するための標準的な診断ツールである。しかし、その解釈は通常、熱力学的極限（$N \to \infty$）においてのみ検証される。小さなスピンクラスターは、分子トリマーや少数の量子ビットシミュレータとして物理的に重要であるが、熱力学的極限を持たないため、標準的な曲率の解釈（臨界点での発散や相関長へのスケーリングなど）を適用できない。課題は、相互作用のトポロジーを連続的に変化させることが可能な有限系の正確な幾何学的特性を決定し、数値近似の曖昧さなしに、フラストレーション、縮退、および結合異方性がどのように曲率に現れるかを理解することである。

手法
著者らは、以下のハミルトニアンに従う3つのイジングスピン（$S_i = \pm 1$）の系を分析する：
$$H = -J (S_1S_2 + u S_2S_3 + a S_1S_3) - h (S_1 + S_2 + S_3)$$
ここで、$J$ は結合スケール、$h$ は縦磁場、$u, a$ は結合強度を定義する無次元比である。このモデルは、開いた鎖（$a=0$）から対称な三角形（$a=u=1$）までを補間する。

本研究では以下のステップを採用している：

1. 厳密な分配関数： すべての構成に対して、無次元変数 $x = \beta J$ および $y = \beta h$ を用いて、分配関数 $Z$ を解析的に導出する。
2. 計量の構築： フィッシャー・ラーオ計量 $g_{ij}$ を、自然パラメータ $(x, y)$ に関する $\ln Z$ のヘッセ行列として構築する。この計量は、結合変数 $B$ と磁化 $M$ の共分散行列に対応する。
3. 曲率の計算： 2次元ヘッセ行列の行列式公式を用いて、スカラー曲率 $R$ を計算する。式は、$B$ と $M$ の2次および3次の同時累積モーメントのみに依存するように簡約される。
4. 極限の分析： 対称な三角形（$a=u=1$）を分析して閉形式の解を得た後、弱い結合異方性（$a = 1+\delta$）によって対称性を摂動させ、基底状態の縮退がどのように解消されるかを研究する。最後に、開いた鎖から三角形への幾何学的進化を追跡する。

主な貢献と結果

• 対称な三角形の閉形式解：
  対称な三角形について、曲率は以下のコンパクトな閉形式で導出される：
  $$R_\triangle(x, y) = -8 \frac{e^{4x} \sinh(4x) \cosh^2 y (\cosh 2y - 2)}{(3e^{4x} + 2 \cosh 2y - 1)^2}$$
  $R$ の符号は、2つのゼロ線によって決定される：$x=0$（非相互作用極限）および $y=y_c$（ここで $\cosh 2y_c = 2$、すなわち $y_c \approx 0.658$）。$y=y_c$ の線は水平であり、温度や結合強度に依存しない。

• ゼロ線の解釈：

  • $x=0$： 非相互作用のスピンに対して期待される通り、曲率は消失する。
  • $y=y_c$： この線は単一スピン感受率の変曲点に対応する。曲率は、単一スピン磁化分布がプラティクルティック（平坦）からレプトクルティック（尖鋭）へと遷移する場所で正確に符号を変える。
  • スケーリング挙動： $y_c$ 付近において、曲率は線形に振る舞う（$R \propto |y-y_c|^1$）。著者らは、この指数 $\gamma=1$ が、臨界スケーリングからではなく、有限系の解析性に由来する厳密なものであることを強調している。

• フラストレートした反強磁性領域：
  対称な反強磁性三角形（$J<0, a=u=1$）では、システムは $h=0$ において6重の基底状態縮退を示す。曲率はゼロ温度まで負の有限値を保ち、「プラトー」を形成する。

  • 低温公式： 一般的な公式により、このプラトーと、縮退した基底多様体と第1励起状態との間の磁化成分の変化との関係が示される。
  • 異方性への感度： 弱い結合異方性（$\delta \neq 0$）を導入すると、縮退が解消される。曲率のプラトーは、異方性に依存する対数的なクロスオーバー温度 $T^*$ 以下で破壊される：$k_B T^* \propto 2|J| / \ln(1/|\delta|)$。これは、$R$ が正確な基底状態の縮退に対して指数関数的に敏感であることを示している。

• 強磁性領域と異方性：

  • 普遍性： 強磁性結合（$J>0$）の場合、ゼロ温度における曲率は普遍的であり、単一スピン感受率の変曲点のみによって決定される極限に接近する。
  • 励起レベルの交差： 異方性を持つクラスターでは、強磁性セクターにおいて追加のゼロ線が現れる。これらの線は、磁場誘起による励起エネルギー準位（例：1スピン反転状態と2スピン反転状態）の交差に対応し、$h/J$ の固定された比に固定されている。

• 開いた鎖 vs 三角形：
  開いた鎖（$a=0$）はフラストレーションを起こさない。その曲率は、反強磁性領域の低温において指数関数的に小さい。鎖から三角形への遷移は、第3の結合（$a>0$）を導入することで、ゼロ線の構造がいかにして温度依存の傾斜を持つ線から、低温度における水平線へと変換され、フラストレーションの発生を知らせるかを明らかにしている。

意義と主張
本論文は、3スピンクラスターが、有限系における熱力学的曲率を解釈するための厳密なベンチマークとなることを主張している。主な結論は以下の通りである：

1. 解析性と臨界性： 有限系において、曲率の符号変化は解析関数の単純なゼロであり、指数は正確に1である。ゼロ付近で見られるべき冪乗則を、臨界スケーリングの証拠と誤認してはならない。
2. 縮退の検出： 不安定な（フラストレートした）領域における有限の負の曲率プラトーは、正確な基底状態の縮退と、第1励起における磁化成分の変化が組み合わさった幾何学的な指紋である。この特徴は対称性の破れに対して極めて敏感であり、微小な異方性であっても、エネルギー分裂自体の温度よりもはるかに高い温度でプラトーを破壊する。
3. 多様体への依存性： 非相互作用系において曲率が消失するという記述は、幾何学的多様体が非相互作用のパラメータから構築されている場合にのみ成立する。相互作用する方向（例えば $J$ を変化させる方向）を含む多様体では、共役な統計の3次相関により、 $J=0$ においても曲率が非ゼロになり得る。
4. 実験的関連性： これらの正確な結果は、独立した分子トリマーの集合体や捕捉イオン量子シミュレータに対する、テスト可能な厳密なシグネチャを提供する。具体的には、温度に依存しないゼロ線 $\cosh 2y_c = 2$ と、異方性によって制御されるフラストレーション・プラトーの崩壊が、磁化および感受率のデータを通じて測定可能なシグネチャとして提案されている。

著者らは、これらの正確な結果が、小さなスピン系における相互作用、フラストレーション、および縮退を診断するために熱力学的曲率を用いるための、制御された参照を提供すると結論づけている。