[論文レビュー] Low-frequency vibrational modes of stable glasses

Lijin Wang, Andrea Ninarello+4 more|RePEc: Research Papers in Economics|Apr 23, 2018

Spectroscopy and Quantum Chemical Studies被引用数 5

ひとこと要約

本研究は、非常に冷却された状態から超安定状態まで広い運動的安定性を示すコンピュータ生成ガラスにおける低周波数振動モードを、モンテカルロ平衡化を用いて実験的に関連する状態に到達させることで調査した。その結果、準局在モードがすべての安定性において普遍的な ω⁴ 定律に従うことが明らかになった。また、安定性が高くなるにつれてその密度は急激に減少し、ガラスの安定性と軟らかく局在した振動モードとの直接的な関連が確立された。

ABSTRACT

We numerically study the evolution of the vibrational density of states $D(ω)$ of zero-temperature glasses when their kinetic stability is varied over an extremely broad range, ranging from poorly annealed glasses obtained by instantaneous quenches from above the onset temperature, to ultrastable glasses obtained by quenching systems thermalised below the experimental glass temperature. The low-frequency part of the density of states splits between extended and quasi-localized modes. Extended modes exhibit a boson peak crossing over to Debye behaviour ($D(ω) \sim ω^2$) at low-frequency, with a strong correlation between the two regimes. Quasi-localized modes instead obey $D(ω) \sim ω^4$, irrespective of the glass stability. However, the prefactor of this quartic law becomes smaller in more stable glasses, and the corresponding modes become more localized and sparser. Our work is the first numerical observation of quasi-localized modes in a regime relevant to experiments, and it establishes a direct connection between glass stability and soft vibrational motion in amorphous solids.

研究の動機と目的

• コンピュータガラスにおける広い運動的安定性の範囲に到達することで、シミュレーションによるガラスの性質と実験的観察とのギャップを埋めること。
• 準局在モードの低周波数振動状態密度 D(ω) が、さまざまなガラス安定性においても普遍的であるかどうかを特定すること。
• 運動的安定性が高くなるにつれて、軟らかく局在した振動モードの密度と局在性がどのように変化するかを調査すること。
• ガラスの安定性と低周波数振動モードの性質、特に準局在モードとの直接的な関連を確立すること。
• 実際の実験的ガラスに相当する領域で、以前のシミュレーション結果における準局在モードの妥当性を検証すること。

提案手法

• 密度 ρ = 1.0 で、N = 48,000 から 450,000 個の粒子を含む多分散 Lennard-Jones 型モデルガラスを用いた。
• 母体温度 Tp を Tp = ∞ から Tp = 0.062 まで低下させ、モード結合温度 Tc ≈ 0.108 よりも低い範囲で、スワップモンテカルロシミュレーションを用いて過冷液体を平衡化した。
• 即時クエンチングを FIRE アルゴリズムを用いて行い、零温度における固有構造を生成した。
• ヘシアン行列の対角化により振動モードを計算し、状態密度 D(ω) = (1/(3N−3)) Σ δ(ω − ωl) および参加比 P(ωl) を算出した。
• 参加比 P(ωl) を基準にモードを拡張的か準局在的かに分類し、拡張的モードでは P ≈ 2/3、局在的モードでは P ≈ 1/N となるようにした。
• P(ω) 対 ω のヒストグラムとピーク検出を用いて拡張的モードの寄与 D_ex(ω) を分離・分析し、D_loc(ω) を ω⁴ スケーリングにフィットした。

実験結果

リサーチクエスチョン

• RQ1 運動的安定性が実験的ガラスと同等のガラスにおいて、準局在振動モードの ω⁴ スケーリングが成立するか？
• RQ2 D_loc(ω) = A4ω⁴ における比例定数 A4 は、ガラス安定性が向上するにつれてどのように変化するか？
• RQ3 ガラス安定性が向上するにつれて、軟らかく局在したモードの局在性および空間的分布に体系的な変化があるか？
• RQ4 低周波数振動性質が、平衡化プロトコルや母体温度 Tp に依存する程度はどの程度か？
• RQ5 観察された D(ω) の変化が、Tp で定量化されたガラスの運動的安定性と直接関連しているか？

主な発見

• すべてのガラス安定性において、低周波数振動状態密度 D(ω) は拡張的モードと準局在モードに分離する。
• 拡張的モードは低周波数領域で D_ex(ω) ∼ ω² のデュイの法則的挙動を示し、デュイ領域とボーズピーク領域の間に強い相関がある。
• 準局在モードは、ガラス安定性に依存せず、普遍的に ω⁴ スケーリング則 D_loc(ω) ∼ ω⁴ に従う。
• ω⁴ 定律の比例定数 A4 は母体温度 Tp が低下するにつれて急激に減少し、より安定なガラスでは準局在モード密度が強く抑制されていることが示された。
• 安定性が高くなるにつれてモードはより局在的かつ希釈的になり、参加比の低下によって確認された。
• 本研究は、実験的ガラス安定性に一致する領域における準局在モードの最初の数値的証拠を提供し、以前のシミュレーション結果を物理的に意味のある文脈で検証した。