Equivalence of residual entropy of hexagonal and cubic ices from tensor network methods

高精度なテンソルネットワーク法を用いてアイス則を明示的に符号化し、転送演算子の正規性を検証することにより、本研究は、六角形アイスと立方体アイスの残留エントロピーが等しいという厳密な数値的証拠を提示する。

要約

あなたは、水分子で作られた巨大な三次元のパズルを想像してみてください。このパズルでは、酸素原子が硬い骨格を形成していますが、水素原子は、酸素原子の間にあるさまざまな場所に座ることができる、小さくて落ち着きのないゲストのような存在です。

これらのゲストには、座り方に関する厳格なルール（「アイスルール」と呼ばれます）があります。すなわち、すべての酸素原子には、ちょうど2つのゲストが近くに座り、さらに2つのゲストが少し離れた場所に座っていなければなりません。温度が非常に低くなり、水が固体の氷に変わったとしても、ゲストたちは一つの完璧な配置に凍りつくことはありません。代わりに、彼らはルールに従いながら、何兆通りもの異なる方法でシャッフルし続けることができます。

この「シャッフル」は、残留エントロピーとして知られる、残された無秩序を生み出します。科学者たちは数十年にわたり、ある特定の問いについて議論してきました。それは、**「無秩序の量は、氷の骨格の形によって異なるのか？」**という問いです。

氷には主に2つの形状があります：

1. 六方晶氷 (Ih): 自然界で最も一般的な形態（雪の結晶など）。
2. 立方晶氷 (Ic): わずかに異なる3D構造を持つ、より珍しい形態。

長年、数学者たちは、六方晶氷は必ず立方晶氷と同等以上の無秩序を持つはずである（$S_h \ge S_c$）ことを証明してきました。しかし、コンピュータ・シミュレーションは、その数値があまりにも近いため、実際には同一である可能性を示唆していました。問題は、これらを検証するためにコンピュータが使用した手法（「モンテカルロ法」と呼ばれます）が、ランダムに推測することであらゆるシャッフルを数えようとするものであり、全体像を明確に見通すことができなかった点にありました。

新しいアプローチ：「テンソルネットワーク」のレンズ

論文の著者たちは、強力な新しい数学的ツールであるテンソルネットワークを使用しました。これは、単に答えを推測するのではなく、可能性の全景を一括してマッピングする高精細なレンズのようなものです。

ゲストをランダムにシャッフルする代わりに、彼らは数学的な「転送マシン」（転送演算子と呼ばれます）を構築しました。このマシンは、氷の一層を取り込み、ルールを適用し、それを次の層へと渡していきます。このマシンから出てくる「最も強い信号」（最大固有値）を見つけることで、彼らは推測することなく、無秩序の正確な量を計算することができました。

大発見： 「鏡」テスト

ここが彼らの発見の巧妙な部分です。二種類の氷が全く同じ無秩序を持つためには、立方晶氷に使用される数学的なマシンが、非常に特殊な挙動、すなわち「正規（normal）」である必要がありました。

簡単に言えば、「正規」のマシンとは、そのステップを実行する順序によって最終的な結果が変わらないもののことです。それは、光を完璧に反射する鏡のようなものです。正面から見ても横から見ても、反射は一貫しています。

著者たちは、立方晶氷のマシンが「正規」であるかどうかを確認するために、高精度なテストを実施しました。その結果、それは99.99% 正規であることが分かりました。完全な鏡ではありませんが（ごくわずかな欠陥があります）、実用上は完璧であるかのように振る舞います。

最終結果

マシンが「正規」に極めて近いため、著者たちは、数値を特定の形に無理やり合わせる（以前の研究者が行っていたテクニック）ことなく、直接計算を行うことができました。

計算を行った結果：

• 六方晶氷の無秩序（$S_h$）は 0.4104251 となりました。
• 立方晶氷の無秩序（$S_c$）は 0.4104248 となりました。

これら二つの数値の差は極めて小さく（約500万分の1）、計算手法による微小な誤差である可能性が高く、物理学的な差異ではありません。

結論

日常的な言葉で言えば、六方晶氷と立方晶氷は、全く同じ量の残された無秩序を持っています。

著者たちは単に推測したのではなく、洗練された数学的な「レンズ」を用いて、二種類の氷を支配するルールがあまりにも似ているため、結果として同じレベルの混沌をもたらすことを証明し、物理学における長年の論争に終止符を打ちました。また、彼らは、この手法が最近科学者によって発見された他の奇妙な形態の氷の研究にも使用できる可能性があることも指摘しましたが、それは将来の研究の課題です。

技術概要

技術要約：テンソルネットワーク法による六方晶氷と立方晶氷における残留エントロピーの等価性

問題提起
本論文は、水の残留エントロピーに関する統計物理学における長年の未解決問題に取り組んでいる。具体的には、六方晶氷の残留エントロピー（$S_h$）が立方晶氷の残留エントロピー（$S_c$）と等しいか否かを調査している。解析的研究では不等式 $S_h \geq S_c$ が確立されている一方で、モンテカルロ法や級数展開を用いた数値的研究では、両者の値が極めて近いことが一貫して示唆されており、厳密な等価性の問題は未解決のまま残されていた。この問題を解決する上での主な課題は、標準的なモンテカルロ法では、残留エントロピーを得るために基底状態の縮退空間を直接サンプリングすることができず、代わりに有限温度における比熱の積分に依存しなければならない点にある。さらに、これまでのテンソルネットワークによる試みは、立方晶氷の転送演算子が非エルミートであるため、人工的な対称性の制約を課さざるを得ず、それが二つの氷相の区別を不明瞭にしていた。

手法
著者らは、高精度なテンソルネットワーク法を用いて、残留エントロピー問題を再定式化している。

1. テンソルネットワーク表現： 「アイスルール」（酸素の近くに水素原子が2つ、遠くに2つある状態）を局所テンソルにエンコードしている。六方晶（$I_h$）および立方晶（$I_c$）氷の3次元格子構造は、全構成数が層間転送演算子（$\hat{T}$）の最大固有値によって決定されるテンソルネットワークへとマッピングされる。
   • 立方晶氷の場合、$\hat{T}_c = \hat{M}$ である。
   • 六方晶氷の場合、$\hat{T}_h = \hat{M}\hat{M}^T$ である。
2. 正規性（Normality）の解析： 著者らは、立方晶氷の転送演算子 $\hat{M}$ の「正規性」（行列 $\hat{M}$ が $\hat{M}\hat{M}^T = \hat{M}^T\hat{M}$ を満たすとき、その行列は正規である）を調査している。著者らは、Variational Uniform Matrix Product State (VUMPS) アルゴリズムを用いて、$\hat{M}\hat{M}^T$ と $\hat{M}^T\hat{M}$ の間のフィデリティ（忠実度）を数値的に計算している。高い正規性は、$\hat{M}\hat{M}^T$ の固有値が $\hat{M}$ の固有値のノルムの二乗に対応することを意味し、理論的に $S_h = S_c$ を支持するものである。
3. 変分最適化： 残留エントロピーを直接計算するために、著者らは最近開発された split Corner Transfer Matrix Renormalization Group (split-CTMRG) アルゴリズムを利用している。重要な点として、$\hat{M}$ は非負であり、かつ高度に正規であるため、従来の研究で使用されていた空間対称性の制約を課すことなく、変分原理（レイリー商）を適用して最大固有ベクトルを見つけることができると彼らは主張している。これにより、$S_h$ と $S_c$ の直接的かつ制約のない比較が可能となる。

主な貢献

• 直接計算： 著者らは、モンテカルロ法で必要とされる有限温度の積分という制限を回避し、テンソルネットワークを用いて基底状態の縮退空間において残留エントロピーを直接計算する枠組みを提供している。
• 正規性の視点： 著者らは、転送演算子の正規性の解析を、理論的かつ数値的なツールとして導入した。彼らは、立方晶氷の転送演算子が「極めて近い」正規性を持っていることを示し、これが人工的な対称性を課さずに変分法を用いることを正当化することを実証した。
• アルゴリズムの適用： 本研究は、非エルミート転送演算子に対して無限投影絡み合いペア状態（iPEPS）を最適化するための split-CTMRG 法を適用しており、高いボンド次元（$D=7$）およびCTM次元（$\chi=150$）を達成している。

結果

• 正規性の尺度： $\hat{M}\hat{M}^T$ と $\hat{M}^T\hat{M}$ の間の数値的なフィデリティは、無限ボンド次元極限へと外挿され、値 $F \approx 0.999903$ を与えた。これは、厳密には完全ではないものの、非常に高い程度の正規性を示している。
• 残留エントロピーの値： split-CTMRG 法を用い、熱力学的極限（$D, \chi \to \infty$）への外挿により、著者らは以下の値を得た。
  • $S_h = 0.4104251(6)$
  • $S_c = 0.4104248(14)$
• 比較： 両者の値の相対的な差は $5 \times 10^{-6}$ のオーダーである。彼らの手法の数値精度範囲内において、結果は $S_h$ と $S_c$ が等しいことを示している。制約のない最適化は、従来の制約付きテンソルネットワーク研究と比較して、六方晶氷に対してわずかに高い値を示しており、これはより大きな変分空間の重要性を裏付けている。

意義と主張
本論文は、$S_h = S_c$ という等価性を支持する強力な数値的証拠を提供すると主張している。本研究の意義は以下の通りである。

1. 不等式の解決： 解析的な下限（$S_h \geq S_c$）を超え、そのギャップが現在の計算精度において事実上ゼロであることを示す数値的証拠を提示した。
2. 手法の進展： 転送演算子の正規性が、人工的な対称性を課すことなく非エルミート系に変分テンソルネットワーク法を適用するための十分条件であることを示した。これは、層間転送行列がエルミートではない3次元統計物理モデルを研究するための新しい視点を提供する。
   3.将来の適用可能性： 著者らは、このアプローチが、酸素原子が異なる規則的な格子を形成しながらも配位数4を維持する、他のより複雑な氷相（20種類以上存在する）の残留エントロピーを調査する上で有望であると示唆している。

著者らは謙虚な姿勢を崩さず、高い正規性が等価性を支持しているものの、厳密な正規性の微小な違反は数学的に等価性を否定するものではないが、彼らが用いたような高精度な数値的アプローチが必要であることを示していると述べている。