The anisotropic square lattice Ising model with quenched surface disorder

技術要約：無秩序な表面不純性を伴う異方性正方格子イジング模型

問題提起
本研究では、異方的な結合を持つ、長さ $L$、周囲 $M$ の円筒上の正方格子イジング模型を調査する。この系は、一方の側に開端境界条件（BC）を備え、反対側の側には、クエンチされたランダムな表面不純性を有している。不純性は、無限に強い結合を介して追加の層 $L+1$ に配置された境界スピン $\epsilon_m = \pm 1$ によってモデル化されている。本研究の目的は、任意の固定境界条件に対して、すべての可能な不純物構成について平均化を行い、分配関数、自由エネルギー、および臨界カシミール寄与を厳密に計算することである。主な目標は、二次元における表面不理の関連性に関する予測を検証するために、熱力学的極限（$L, M \to \infty$）におけるランダム表面自由エネルギーおよび臨界カシミール振幅を決定することである。

手法
著者らは、以前の開端境界および周期境界に関する手法を、任意の固定境界条件を持つ系へと拡張した、転送行列法に基づく厳密な列挙法を用いている。

• 転送行列定式化: 分配関数は、$2M \times 2M$ サイズのブロック転送行列 $T = U_\updownarrow U_\leftrightarrow$ を用いて導出される。境界の不純物は、境界結合の符号を表す対角行列 $\kappa$ を通じて組み込まれる。無限の境界結合強度の極限を取ることで、分配関数は、不純物構成を含む行列和の行列式を用いて $\det(Q + \kappa)$ と表される。ここで、$Q$ は転送行列の固有値から導かれる歪循環行列である。
• 厳密な平均化: モンテカルロ・サンプリングの代わりに、著者らはすべての $2^M$ 通りの境界不純物構成に対する厳密な平均化を行っている。これは、ウッドベリー行列恒等式およびウッドベリー・ツリー・アルゴリズムを用いることで達成される。これらのアルゴリズムにより、逐次的またはツリーベースの行列式の計算が可能となり、計算複雑度は $O(2^M)$ となり、総当たり的な行列式計算に要する $O(M^3 2^M)$ よりも大幅に高速化されている。
• 外挿: 熱力学的極限に到達するために、有限の $M$（最大 $M=54$）に対する数値的に厳密な結果を、Aitkenの $\Delta^2$ 法、Levin-U変換、一般化されたBulirsch-Stoer法、および高次の補正項を打ち消すように設計された一般化差分演算子（$\Psi$ 法）を含む、いくつかの系列加速技術を用いて外挿している。
• スケーリング解析: 研究は、異方的な臨界点およびハミルトニアン極限（$z_c \to 1$）に焦点を当てている。スケーリング挙動は、減少したアスペクト比 $\rho = L/(M r_\xi)$ （ここで $r_\xi$ はバルク相関長さの振幅比）に関して解析される。

主要な貢献

1. 不規則な境界に対する厳密解: 本論文は、任意のクエンチされた表面不純性を伴うイジング円筒の分配関数および自由エネルギーを計算するための、厳密な解析的枠組みを提供し、モンテカルロ・シミュレーションに固有の統計誤差を回避している。
2. 効率的なアルゴリズム: ウッドベリー・ツリー・アルゴリズムの実装により、最大 $M=54$ のスピン（$2^{54}$ 通りの不純物構成をカバー）を持つ系の自由エネルギー計算が可能となり、これは従来のシミュレーションベースのアプローチよりも数桁高速である。
3. 熱力学的状態密度 (DOTS): 著者らは、不純物アンサンブルにおける自由エネルギーの分布を特徴付け、残留過剰自由エネルギーが熱力学的極限において対数正規分布に従うことを見出した。
4. 表面自由エネルギーの計算: ランダム表面自由エネルギー密度 $\bar{f}_s^{(r)}$ は、等方的結合および最大異方性結合の両方に対して、高い精度（小数点以下6桁）で計算されている。

結果

• 表面自由エネルギー: 計算されたランダム表面自由エネルギー密度 $\bar{f}_s^{(r)}$ は非ゼロであり、結合の異方性に依存する。等方的結合（$z_c^{iso}$）では $\bar{f}_s^{(r)} \approx 0.002662$ であり、ハミルトニアン極限（$z_c \to 1$）では $\bar{f}_s^{(r)} \approx 0.26536$ である。これらの結果は、不純物アンサンブルがすべての構成について平均化されているか、あるいは境界磁化がゼロ（$m_B=0$）に制限されているかに依存しない。
• カシミール振幅: ランダム不純物の場合とスタッガード（ディリクレ）の場合の差を表す過剰カシミール振幅 $\Delta_C^{(ex)}$ は、$\rho \to \infty$ において数値的精度内でゼロと互換性があることが判明した。具体的には、等方的結合に対しては $\lim_{\rho \to \infty} \Delta_C^{(ex)}(\rho) \approx (2.8 \pm 5) \times 10^{-4}$、ハミルトニアン極限に対しては $\approx (1.9 \pm 3) \times 10^{-5}$ である。
• 普遍性クラス: ランダム・アンサンブルおよびゼロ磁化アンサンブルのカシミール力スケーリング関数 $\vartheta$ は、開端・開端 $(o,o)$ 普遍性クラスと区別できないことが示された。逆に、非ゼロの固定磁化（例：$m_B = 3/4$）を持つアンサンブルは、開端プラス $(o,+)$ 普遍性クラスに一致する。
• 対数補正: 不純性は、カシミール振幅が消失するという意味において（$d=2$ におけるハリスの判定基準と一致しており）無関連であるように見えるが、データは、特にアスペクト比 $\rho$ への依存性において、スケーリングに対する対数補正の存在を示唆している。

意義
本論文は、不純物を伴うイジング模型の共形場理論の予測に対し、厳密な検証を提供すると主張している。ランダム表面自由エネルギーが明確に定義されており、かつ臨界カシミール振幅が消失すること（これは2次元における不純物が無関連であることを示唆する）を示すことで、本研究は、開端境界を持つイジング模型には、開端 $(o)$ クラスとプラス $(+)$ クラスの2つの表面普遍性クラスしか存在するという仮説を支持している。著者らは、厳密なアプローチを用いることで、モンテカルロ・シミュレーションにおける有限サイズ効果や統計的サンプリングの限界を回避していることを強調している。対数補正の観察は、不純物が無関連であるとしても、それは純粋な系とは異なる微妙な形でスケーリング挙動を修飾することを示唆している。