Measuring structural restraint in an amorphous material using the variance of speckle intensity

技術要約：スペックル強度の分散を用いた非晶質材料における構造的拘束の測定

問題提起
ガラス転移の根本的な謎は、過冷却液体から非晶質固体へと至る劇的な機械的変化と、静的な構造（例：X線散乱振幅）に見られる微小かつ連続的な変化との間の乖離にある。特定の「好ましい」局所配位幾何学と構造を結びつけようとする従来のアプローチは、非晶質材料における局所構造の多様性を考慮すると、構造と機械的特性の間の因果関係を確立する上で困難に直面している。著者らは、冷却時に重要な変化が起こるのは構造そのものではなく、構成原子の運動を抑制する「構造の能力」であると考えている。しかし、既存のシミュレーションに基づく構造的拘束の定義（ポテンシャルエネルギーの調和近似に基づくもの）は、実験的にアクセス可能ではない。課題は、個々の反射強度の長期的かつ連続的な追跡を必要とせずに（これは実験的な不安定性を招きやすい）、実験的な散乱データから導出できる構造的拘束の尺度を定義することである。

手法
著者らは、二成分レナード・ジョーンズ混合物（Kob-Andersenモデル）の分子動力学シミュレーションを利用して、構造的拘束の新しい指標である $R(q,T)$ を開発し、テストした。

1. 拘束の定義： この尺度は、コヒーレント放射の弾性散乱から生じるスペックルパターンの強度の分散に基づいている。特定の波数の温度依存性を経時的に追跡する必要がある従来の提案とは異なり、本手法は各制御パラメータ（温度）における散乱を独立して扱う。
   構造的拘束 $R(q,T)$ は、温度 $T$ における散乱強度の角度分散を、参照温度（固体を代表するのに十分な低温度）における分散との比として定義される：
   $$R(q,T) = \frac{\text{Var}[I(q, \phi, T)]}{\text{Var}[I(q, \phi, T_{ref})]}$$
   ここで、強度 $I$ は測定時間 $\tau$ にわたって平均化され、角度分散は固定された散乱ベクトル $q$ の大きさにおける角度 $\phi$ に関して取られる。

2. シミュレーション・パラメータ： シミュレーションは、5000個の粒子を含む系に対し、LAMMPSソフトウェアとNose-Hoover熱浴を用いて実施された。測定時間は $\tau = 2 \times 10^4$（簡約単位）に設定された。研究は構造因子の第1ピーク $q_1 = 7.2$ に焦点を当てているが、他の波数についてもテストが行われた。

3. 比較分析： 提案された構造的拘束 $R(q,T)$ を、原子の拘束性と相関することが期待される2つの確立された物理量と比較した：

   • デバイ・ワラー因子 ($D_0$): 平均二乗変位から算出。
   • せん断弾性率 ($G$): 同一の測定時間 $\tau$ における瞬時せん断応力の分散を含む、Squire, Holt, および Hooverによる揺らぎ公式を用いて算出。

4. 実験的実現可能性のテスト： 著者らは、実用的な実験変数である以下の2点に対する $R(q,T)$ の感度を調査した：

   • 散乱波数 $q$ の選択。
   • 強度データを収集する際の角度的粗視化（画素サイズ）の度合い（ショットノイズや検出器分解能の影響をシミュレート）。

主な結果

• 機械的特性との相関： 構造的拘束 $R(q_1, T)$ は、高温域と低温域の間で大幅な変化を示す。極めて重要なことに、$R(q_1, T)$ の温度依存性は、簡約せん断弾性率 $G(T)/G(T=0.01)$ と驚くほど一致している。デバイ・ワラー因子との相関は正であるが、$0.15 < T < 0.35$ の温度範囲において系統的な偏差が見られる。原子の位置の揺らぎとせん断応力の間の非線形な結合があるにもかかわらず、せん断弾性率との一致が得られたことは驚くべきことである。
• 波数依存性： 構造的拘束が始まる特徴的な温度 $T_c$ （$R=0.5$ となる温度）は、$q$ によって変化する。$T_c$ は $q = q_1$ （構造因子の第1ピーク）において鋭い極大値を示し、拘束の開始がブラッグ波波数における密度揺らぎと最も強く関連していることを示唆している。$q < 5$ では、散乱強度の分散が激減するため、拘束の計算は信頼性が低くなる。
• 粗視化に対する堅牢性： 本研究では、角度的粗視化（より大きな画素サイズ）が増加するにつれて、スペックル強度の生の分散は単調に減少するものの、構造的拘束を定義する比 $R(q,T)$ は比較的安定していることが示された。構造因子の分散が閾値（約0.1）を上回っている限り、測定される拘束はほとんど変化しない。これは、本指標が実際の実験で必要とされる画素平均化に対して堅牢であることを示している。

意義と主張
本論文は、凝縮系における構造的拘束を測定するための、実用的かつ実験的にアクセス可能な手法を提示していると主張している。個々の反射を追跡するのではなく、スペックル強度の分散を利用することで、本手法は試料の配向シフト、光源のコヒーレンス性の変化、および放射損傷に伴う実験上の困難を回避している。

著者らは、$R(q,T)$ が流体と非晶質固体の間の構造転移を定量化するための有効な秩序パラメータとして機能すると主張している。その主要な意義は、弾性散乱のみを用いて、固体的な原子の拘束が始まる温度を特定できる点にある。散乱から導出されたこの拘束とせん断弾性率との間に見られた強い相関は、この手法を用いて、散乱測定から直接、機械的特性が温度とともにどのように変化するかを決定できる可能性を示唆している。また、このアプローチは熱的な制御パラメータに限定されず、密度（コロイドガラス）、印加歪み、あるいは照射や化学変化を受けた材料のような非平衡状態にも拡張可能であると述べている。