Percolation Criticality of Amorphous-Amorphous Transitions in Compressed Glasses

本研究は、大規模な分子動力学シミュレーションとパーコレーション理論を用いることで、圧縮されたシリカガラスにおける低密度から高密度への転移が構造的クラスターの臨界パーコレーションによって駆動されていることを明らかにし、剛性パーコレーション機構を示唆する臨界指数を提示するとともに、結合性および非結合性の両アモルファス材料間における共通の変形原理を強調している。

要約

ガラスのコップや窓ガラスを想像してみてください。それらは硬くて動かないものとして知られています。しかし、もしあなたが原子の大きさにまで縮小してその内部を覗き込んだとしたら、そこには小さな積み木が入り混じった、混沌としたもつれ合ったネットワークが見えるはずです。シリカガラス（窓の材料）では、これらのブロックはシリコンと酸素で作られたピラミッド型（四面体）をしています。

この論文は、このガラスに凄まじい圧力——大気圧の35万倍にも達する圧力——をかけて押しつぶしたときに何が起こるのかをズームアップした、ハイテクな映画のようなものです。科学者たちは、ガラスが溶けたり壊れたりすることなく、どのようにその形を変えるのか、すなわち「アモルファス・アモルファス転移」と呼ばれるプロセスを解明しようとしました。

以下に、彼らの発見の物語を、簡単な比喩を用いて説明します。

1. コンサート会場の群衆（構造）

ガラスを、コンサート会場にいる巨大な群衆だと考えてみてください。

• 通常時の圧力： 全員がゆったりとした、開けた隊列で立っています。シリカガラスの場合、「人々」はシリコン原子であり、彼らは4人の隣人と手をつなぎ、完璧なピラミッド形状（四面体）を形成しています。彼らの間には多くの空隙があり、その構造は「ぐにゃぐにゃ」としていて、圧縮しやすい状態にあります。
• 圧力が増加すると： コンサートホールが縮小し始めたと想像してください。群衆が押しつぶされます。人々はもはや完璧なピラミッド形状を維持できなくなります。彼らは互いにぶつかり合い、何人の隣人と手をつなぐかを変え始めます。ある人は5人と、またある人は6人と手をつなぎ始めます。

2. 「パーコレーション」ゲーム（大きな変化）

科学者たちは、**パーコレーション（浸透）**という概念を用いました。スポンジに水を注ぐ場面を想像してください。

• 低圧時： スポンジには穴がありますが、それらはすべてバラバラです。水を注いでも、水は小さなポケットの中に留まってしまいます。最後まで通り抜けることはできません。ガラスの中では、この「ピラミッド」形状は孤立した島のようなものです。
• 臨界圧力： もっと強く押しつぶしていくと、穴同士がつながり始めます。突然、スポンジの上から下まで続く、巨大で連続的な経路が形成されます。水が流れ通るのです！
• ガラスの中では： 科学者たちは、特定の圧力点において、新しい形状（5角形や6角形のブロックなど）が突然つながり、ガラス全体に広がる巨大で連続的な鎖を形成することを発見しました。これが「パーコレーション転移」です。これは、ガラスが根本的に再編成され、より高密度の状態へと移行する瞬間です。

3. 群衆を見る2つの視点

研究者たちは、2つの異なるカメラアングルから映画を観るように、この群衆を2つの異なる方法で観察しました。

• 「結合」による視点（握手）： 彼らは、誰が直接手をつないでいるか（化学結合）を見ました。そして、ピラミッド形状がその「握手」を変えていく様子を観察しました。
• 「非結合」による視点（パーソナルスペース）： 彼らは握手のことは無視し、たとえ触れていなくても、誰が誰の近くに立っているのかだけを見ました。これは、人々が手をつないではいないけれど、ただ近くに立っている群衆を見ているようなものです。

驚きの事実： 両方のカメラは、全く同じストーリーを示していました。「握手」の視点も「パーソナルスペース」の視点も、ガラスが同じ順序で変容することを示していました。つまり、まず緩やかな形状が連結し、次に高密度の形状が支配的になるのです。このことは、原子が「手をつないでいる」場合（シリカのように）でも、単に「ぶつかり合っている」場合（凍った水や氷のように）でも、ガラスが変化するルールは普遍的であることを示唆しています。

4. 「魔法の数字」とゲームのルール

科学者たちは、この変容が標準的なルールブック（確率に基づいたゲームのようなもの）に従っているのか、それとも独自の特別なルールを持っているのかを知りたがりました。

• 四面体（4角形の形状）： 元のピラミッド形状（4つの手をつなぐ形）が崩壊するとき、それはまさにランダムな確率ゲームのように崩壊しました。これは「標準的」な挙動でした。
• 高次の形状（5つ、6つ、あるいはそれ以上の手）： 新しい、より高密度の形状が形成され、連結していくとき、それらは標準的なルールを破りました。彼らは異なる、より複雑なルールに従っていました。科学者たちはこれを「剛性パーコレーション」と呼んでいます。それは、群衆が単にランダムにつながったのではなく、構造全体を突如として非常に硬く、強固なものにするような方法でつながったことを意味します。

5. まとめ

この論文は、ガラスを押しつぶすと、単に小さくなるだけでなく、内部構造が新しい、より高密度の「状態」へと再編成される劇的な相転移のようなイベントが起こることを結論づけています。

• この転移は、特定の「臨界」圧力で起こります。
• 新しい構造がどのように連結するかは、ランダムな偶然（古い形状の場合）と、より硬く構造化されたルール（高密度の新しい形状の場合）の混合です。
• この振る舞いはシリカガラスとアモルファス氷において似ており、自然界が異なる種類のガラス状物質を圧力下で再編成するために、同様の「設計図」を使用していることを示唆しています。

要約すると、この論文は、ガラスを押しつぶしたときに、その微視的な「骨格」がどのように折れ、シフトし、そして再構築されるのかを正確に描き出しています。そして、緩くて動きやすいガラスから、高密度で硬いガラスへの転移は、特定の予測可能な転換点を経て起こることを明らかにしています。

技術概要

技術要約：圧縮ガラスにおけるアモルファス・アモルファス転移のパーコレーション臨界性

問題提起
圧縮されたガラスにおける構造変化、特にテトラヘドラル（四面体）構造に基づく材料（$a$-SiO$_2$やアモルファス氷（$a$-H$_2$O）など）におけるアモルファス・アモルファス転移を支配するメカニズムは、依然として重要な課題である。結晶相転移のような鋭い転移とは異なり、これらの転移は、局所的な構造モチーフが異なる多形へと漸進的に進化するプロセスを伴う。これまでの研究では、高配位多面体からなる大規模なパーコレーティング・クラスターの出現が確認されてきたが、ガラスにおけるこれらの転移の臨界性は不明なままであった。具体的には、これらの構造転移が標準的な（ランダムな）無相関パーコレーションの普遍性クラスに属するのか、あるいは剛性パーコレーション（rigidity percolation）のような異なる臨界挙動を示すのかが明らかになっていなかった。ab initio（第一原理）的手法に固有の小さなシステムサイズによる制約から、これまでの臨界指数の推定は困難であった。

手法
これらの制限に対処するため、著者らは大規模な古典分子動力学（MD）シミュレーションを用いた。

• 系とポテンシャル： シリカガラスについてはSHIKペアポテンシャルを用いたLAMMPSパッケージを使用し、アモルファス氷についてはTIP4P力場を用いたGromacsを使用。システムサイズは、シリカについては約1,000から$10^6$原子、氷については最大49,000分子とした。
• プロトコル： サンプルを平衡化し、室温（シリカは300 K、氷は124 K）まで冷却した後、準静的に最大35 GPa（シリカ）または15 kbar（氷）まで圧縮した。
• 構造記述子： 2つの補完的なアプローチを用いた。
  1. 結合モデル（Bonded Model）： 共有結合的なSi-O連結性に基づき、SiO$_Z$多面体（$Z$は配位数）とその連結性（角共有、辺共有、面共有）を解析する。
  2. 非結合モデル（Non-bonded Model）： カットオフ半径内における中心原子の配置（シリカではSi-Si、氷ではO-O）に基づき、SiSi$_Z$およびOO$_Z$多面体を定義する。このアプローチは相互作用を非結合力として扱うことで、結合型ガラス（シリカ）と非結合型ガラス（氷）の直接的な比較を可能にする。
• パーコレーション解析： クラスターを特定するために、Union-Findアルゴリズムを実装した専用のPythonコード（Nexus-CAT）を使用した。パーコレーション特性として、相関長（$\xi$）、秩序変数（$P_\infty$）、平均クラスターサイズ（$\langle S \rangle$）、および最大クラスターサイズ（$S_{max}$）を算出した。
• 臨界指数： 臨界点付近での臨界指数（$\nu, \beta, \gamma$）およびフラクタル次元（$D_f, D_f^*$）を決定するために、有限サイズスケーリング法（データ崩壊手法）を適用した。

主な結果

1. パーコレーション転移のシーケンス： 圧縮下において、一連の明確なパーコレーション・イベントが特定された。$a$-SiO$_2$では、低密度（LD）のテトラヘドラルネットワーク（$Z=4$）が脱パーコレーション（depercolation）する一方で、高密度（HD）クラスター（$Z=5, 6$）が出現し、パーコレートする。具体的には以下の通りである：

   • SiO$_5$クラスターは $\approx$12 GPaでパーコレートする。
   • SiO$_6$クラスターは $\approx$23 GPaでパーコレートする。
   • 「スティショバイト様（Stishovite-like）」クラスター（2つの辺共有接続を持つSiO$_6$）は $\approx$30 GPaでパーコレートする。
     同様のシーケンスが$a$-H$_2$O（LD、HD、およびVHDクラスター）でも観察されたが、これらはシリカよりも約20倍低い圧力で発生しており、物質固有の相互作用によってスケールされた普遍的なメカニズムを示唆している。

2. モデルの相補性： 結合モデル（SiO$_Z$）と非結合モデル（SiSi$_Z$/OO$_Z$）のアプローチは、一貫した補完的な知見を提供した。非結合モデルは構造変化を捉えることに成功し、シリカとアモルファス氷を比較するための統一された枠組みを可能にした。

3. 臨界指数と普遍性：

   • 相関長 $\xi$ は、システムサイズ $L$ に対して $\xi \propto L^a$ （$a \approx 1$）としてスケールし、臨界挙動を確認した。
   • テトラヘドラの脱パーコレーション： 低密度テトラヘドラルクラスター（$Z=4$）の脱パーコレーションは、標準的な3Dランダムパーコレーションに非常に近い臨界指数（例：フィッシャー指数 $\tau \approx 2.12-2.14$）を示した。
   • 高配位クラスターのパーコレーション： 対照的に、高配位クラスター（$Z=5, 6$）およびVHD構造のパーコレーションは、標準的なパーコレーションからの逸脱を示した。フラクタル次元（$D_f$）は標準的な値（2.53）よりも系統的に低く、その指数は異なる普遍性クラスを示唆していた。

4. フラクタル次元： 有限クラスター（ラティス・アニマル）のフラクタル次元 $D_f^*$ は、様々な配位数に対して $\approx 2.29$ から $2.46$ の範囲で計算され、一般に標準的なパーコレーションの値よりも低かった。

意義と主張
本論文は、長距離記述子の導入と大規模シミュレーションの使用により、圧縮ガラスにおけるパーコレーション臨界指数の初の堅牢な推定が可能になったと主張している。主要な意義は、アモルファス・アモルファス転移の臨界挙動が一様ではないという発見にある：

• 本来のテトラヘドラルネットワークの崩壊は、標準的なランダムパーコレーションに従う。
• 高い連結性を持つクラスターが出現する中で、高配位かつ剛性を持つ多面体ネットワーク（高密度化に関連するもの）の出現とパーコレーションは、おそらく異なる普遍性クラスに従う。これは、低連結性の媒体の中に高連結性のクラスターが出現する際に生じる、トポロジー的および弾性的不均一性に起因すると考えられる。

本研究は、これらのパーコレーション転移と、圧縮率のピークや塑性化の開始といったガラスの機械的特性との間の関連性を確立している。さらに、パーコレーション駆動型の構造変化を追跡することが、アモルファス状態と結晶多形（例：$\text{coesite IV}$, $\text{coesite V}$, $\text{stishovite}$）との間の類推の基礎を提供することを提案しており、ガラス状態がこれらの結晶構造に関連するエネルギー・メガベイシンに対応している可能性を示唆している。本研究は、結合型および非結合型双方のガラスシステムにおける圧力誘起変換の理解を統一する上で、非結合アプローチの有用性を強調している。