非球形不純物を用いた非晶質固体における延性-脆性遷移の有限無秩序臨界点

背景と学術的系譜

起源と学術的系譜

非晶質固体の機械的破壊の研究は、産業および日常生活における広範な応用から、極めて重要な分野である。歴史的に、外部変形を受けた結晶性固体の機械的挙動は欠陥の観点から包括的に解明されてきたが、非晶質固体は長距離構造秩序を欠くため、その機械的応答を正確に描写することは困難である。このため、非晶質固体の降伏メカニズムに関する包括的な理解なしに、数十年にわたる厳密な研究が行われてきた。

重要な研究領域の一つは、材料が徐々に変形する（延性）状態から破壊的に破壊する（脆性）状態へと移行する延性-脆性遷移である。特に非熱的準静的ひずみ（AQS）条件下での最近の研究は、この遷移が材料固有の無秩序強度によって制御されており、延性と脆性の降伏を分離する臨界点が存在することを示唆している。しかしながら、このような臨界点の存在は学術界では論争の的となっており、熱力学極限においては大規模系では常に降伏は脆性であると主張する者もいる。

従来の解析手法の根本的な限界は、この遷移を研究するために固有の無秩序を決定論的かつ実験的にアクセス可能な方法で制御できないことにある。例えば、分子シミュレーションで用いられる熱アニーリングのような手法は、有限サイズ効果に悩まされ、非ブラウン運動するコロイド系には適用できない。粒子固定化も、提案されている手法であるが、分子ガラスへの実装は依然として困難である。さらに、微量合金化（少量の不純物の添加）は材料強度を高めるための一般的な工業的実践であるが、延性-脆性遷移に対するその効果の背後にある微視的メカニズムは依然として不明瞭である。無秩序を調整し、この臨界遷移をプローブするための、堅牢で実験的に実行可能で微視的に理解された方法論の欠如が、著者らを動機づけた痛点である。彼らは、これらの限界に対処するために非球形不純物を用いた新しいプロトコルを提供することを目指している。

直感的領域用語

• 非晶質固体 (Amorphous Solids): ランダムに積み重ねられたレゴブロックの山を想像してください。壁のようにきれいに積み上げられたものではありません。固体ですが、完全に構築されたレゴの家（結晶性固体）とは異なり、繰り返しパターンや組織化された構造はありません。
• 延性-脆性遷移 (Ductile-to-Brittle Transition): 柔らかい粘土と乾燥したクラッカーを考えてみてください。粘土（延性）は、壊れる前にかなりの量伸びて成形できますが、クラッカー（脆性）は、ほとんど変形せずに突然パキッと割れます。この遷移は、材料がその破壊モードにおいて、粘土のような性質からクラッカーのような性質へと変化することを表します。
• 非熱的準静的ひずみ (Athermal Quasistatic Straining, AQS): 温度変化が影響しない部屋で、信じられないほど穏やかで一定の手で、非常にゆっくりとゼリーの塊を押す様子を想像してください。「非熱的」とは熱が関与しないことを意味し、「準静的」とはプロセスが非常に遅いため、系が常にほぼ平衡状態にあり、変形に適応できることを意味します。
• 回転自由度 (Rotational Degrees of Freedom, rDoF): 大理石の入ったボウルに、小さな細長いビーズを追加することを想像してください。大理石は互いに滑るだけですが、細長いビーズは転がったり回転したりすることもでき、移動や再配置のための追加の方法を提供します。これらの追加の「回転」運動が、それらの回転自由度です。
• せん断帯 (Shear Band): トランプのデッキを横方向に押すことを考えてみてください。すべてのカードがスムーズに滑るのではなく、数層が突然非常に速く互いに滑り、ほとんどの変形が発生する局所的な「断層線」または帯を形成し、突然の壊滅的な破壊につながります。

記法表

記法	説明	タイプ
$\gamma$	せん断ひずみ、変形の尺度。	変数
$\sigma_{xy}$	せん断応力、表面に平行に作用する単位面積あたりの力。	変数
$\chi_{dis}$	非連結感受率、応力変動の尺度であり、脆性を示す。	変数
$\Theta$	構造秩序パラメータ、非晶質固体の局所構造秩序または安定性を定量化する。	変数
$S_r$	回転緩和関数、棒状不純物の回転移動度を特徴づける。	変数
$D_{min}^2$	非アフィン変位、巨視的変形では説明されない粒子再配置の尺度。	変数
$c_s$	球形不純物の数分率。	パラメータ
$c_d$	二量体不純物の数分率。	パラメータ
$c_r$	棒状不純物の数分率。	パラメータ
$L_r$	棒状不純物の長さ。	パラメータ
$\sigma_{AA}$	親粒子タイプAの直径。	パラメータ
$\sigma_{BB}$	親粒子タイプBの直径。	パラメータ
$\sigma_s$	球形不純物の直径。	パラメータ
$\sigma_b$	二量体/棒状不純物を形成するビーズの直径。	パラメータ
$T$	温度。	パラメータ
$N$	系中の粒子総数（系サイズ）。	パラメータ
$\Delta \sigma_{max}$	最大塑性降下、延性-脆性遷移の秩序パラメータ。	変数
$\chi_d$	$\Delta \sigma_{max}$変動の感受率、臨界性を示す。	変数
$c_r^*$	臨界棒状不純物分率、延性-脆性遷移が発生する不純物濃度。	パラメータ

問題定義と制約

中核問題の定式化とジレンマ

本論文が取り組む中核問題は、非晶質固体の延性-脆性遷移を理解し、正確に特徴づけることである。

入力/現状: 長距離構造秩序を欠く非晶質固体は、外部変形下で延性降伏（連続的な応力応答を伴う段階的な材料流動）と脆性降伏（不連続な応力降下を伴う突然のせん断帯形成による壊滅的な破壊）という2つの異なる破壊モードを示す。この降伏の性質は、材料固有の無秩序強度によって制御されており、延性挙動と脆性挙動の間には有限の臨界無秩序強度で遷移が存在すると考えられている。微量合金化（少量の不純物の添加）は、材料の機械的特性を向上させる一般的な工学的実践であるが、非晶質固体におけるこれらの改善を支配する微視的メカニズムは十分に理解されていない。

出力/目標状態: 本論文は、延性降伏と脆性降伏の境界を示す有限無秩序臨界点を正確に定義し、特徴づけることを目指す。これは、様々な回転自由度（rDoF）を持つ非球形不純物を導入することによって、固有の無秩序を系統的に制御することによって達成される。最終的な目標は、非晶質固体、特にコロイド系のようなソフトガラスの機械的特性を調整し、望ましい延性または脆性応答を誘発するための、新規で実験的にアクセス可能なプロトコルを確立することである。

欠落しているリンクまたは数学的ギャップ:
1. 論争のある臨界点: 以前の研究では有限無秩序臨界点を示唆していたにもかかわらず、その存在と性質、特に熱力学極限におけるものは、文献において依然として論争の的となっている。熱アニーリングや弾塑性モデルのような以前の手法は、不確かな結果をもたらしたり、有限サイズ効果に悩まされたりした。
2. 微量合金化の微視的メカニズム: 微量合金化、特に非球形不純物によるものが、非晶質固体の降伏遷移と機械的安定性にどのように影響するかという正確な微視的メカニズムは、ほとんど知られていない。非球形不純物によって導入される追加の回転自由度の役割は、重要な未踏の側面である。

痛みを伴うトレードオフまたはジレンマ:
中心的なジレンマは、不純物の特性、回転自由度、および機械的応答の間の複雑な関係にある。不純物の添加は一般的に機械的強度を向上させることができるが、脆性または延性への顕著なシフトを達成することは容易ではない。球形不純物は、降伏ひずみの改善は最小限である。非球形不純物、特に棒状のものは、延性を大幅に向上させ、耐荷重能力を高めることができる回転自由度を導入する。しかし、本論文は直感に反するトレードオフを明らかにしている：この回転自由度を減少させること（不純物の縦横比を増加させるか、その回転を人工的に凍結することによって）は、より脆性な応答と超安定性様の機械的挙動につながる。したがって、研究者は、rDoFを介して延性を向上させる不純物と、rDoFを制限することによって脆性を誘発する不純物の間で、臨界遷移の正確な制御を困難にするジレンマに陥っている。

制約と破壊モード

著者らは、この問題を解決することを非常に困難にする、いくつかの厳しい現実的な壁に遭遇した。

• 物理的/構造的制約:
  • 長距離秩序の不在: 欠陥を正確に描写できる結晶性固体とは異なり、非晶質固体は長距離構造秩序を欠いており、欠陥の定義と追跡が本質的に困難であり、それによって機械的応答の理解が複雑になる。
  • 標準アニーリング手法の不適用性: 無秩序を調整するための従来の熱アニーリングや蒸着などの方法は、すべての系、特に非ブラウン運動するコロイドガラスには効果的または適用可能ではない。これは、これらの系における延性-脆性遷移の実験的アクセス可能性を制限する。
• 計算上の制約:
  • 大きな不純物の遅いダイナミクス: 非常に大きな球形不純物または大径二量体不純物を含む系は、極めて遅いダイナミクスを示す。これにより、現実的なシミュレーション時間スケール内でそのような系を平衡化することが計算上困難になり、研究可能な不純物サイズの範囲が制限される。
  • 棒状長さに対するシミュレーション限界: シミュレーション可能な棒状長さ ($L_r = 2.5\sigma_{AA}$) には実用的な上限がある。これにより、自然に回転拡散率がゼロになる非常に長い棒状の完全な回転停止を直接観測することができなくなる。
  • 長い棒状系のアニーリングの困難さ: より長い棒状の系はダイナミクスが遅いため、標準的な準備プロトコルはアニーリングが不十分な状態をもたらす。このシミュレーションの困難さは、そのような系で観測される脆性の過小評価につながる可能性がある。
  • 自由回転自由度に対する計算コスト: 完全な回転自由度を許容しながら、非常に大きな縦横比の不純物を含む系をシミュレートすることは計算集約的である。これを回避するために、著者らは非常に長い、回転不可能な棒状の効果を模倣するために、rDoFを人工的に凍結することに頼った。
• データ駆動型/分析的制約:
  • 不明瞭な微視的メカニズム: 微量合金化が降伏ひずみを改善できるという実験的証拠にもかかわらず、この改善の背後にある正確な微視的メカニズムは不明瞭であり、詳細な計算分析によって明らかにされる必要がある。
  • パラメータの限定された範囲: 本研究は主に縦横比と回転自由度に焦点を当てている。粒子形状（例：楕円形対棒状）や境界相互作用の性質（例：棒状を形成する球の重なりの程度）などの他の重要なパラメータは認識されているが、複雑さのため将来の研究に残されている。
  • 非熱的準静的仮定: 現在の解析は、非熱的準静的ひずみ（AQS）条件下で実行される。著者らは、結果が熱効果や有限ひずみ率の包含によって変化する可能性があることを認めているが、これらは本論文では十分に探求されていない。

なぜこのアプローチなのか

選択の必然性

著者らが、特に回転自由度が凍結された非球形、棒状不純物を用いた分子動力学（MD）シミュレーションを採用するという決定は、単なる好みではなく、既存手法の限界と問題の特定の要件によって駆動される戦略的な必要性であった。中核的な課題は、非晶質固体の延性-脆性遷移を調査し、実験的にアクセス可能で固有無秩序に対する正確な制御を提供するプロトコルを使用して、有限無秩序臨界点の存在を確立することであった。

従来の「最先端」（SOTA）手法、例えば熱アニーリングは、いくつかの理由で不十分であることが証明された。熱アニーリングは無秩序強度を変化させることができるが、以前の研究[15]は「有限サイズ効果のために部分的な支持しか提供しなかった」し、大規模弾塑性モデル[22]からの結論は「決定的ではなかった」。より決定的なのは、熱アニーリングは「非ブラウン運動するコロイド系には適用できない」ことであり、これらはガラス遷移を研究するための重要な実験プラットフォームである。同様に、粒子固定化[28]は無秩序を制御する別のルートを提供し、コロイドガラス実験[29]で達成できるが、「分子ガラスでは依然として困難」である。著者らは明確に、「特に実験的にアクセス可能なものを含む、この遷移を調査するための新しいプロトコルを特定することは、この分野に大きな影響を与える可能性がある」と述べている。

非球形不純物、特に棒状のものの導入は、このギャップに直接対処する。それらの縦横比と回転自由度を操作することにより、著者らは「分子ガラスとコロイドガラス実験の両方で容易にアクセス可能」な固有無秩序強度を調整する新しい方法を見出した。さらに、「回転拡散率がゼロになる非常に長い棒状」をシミュレートする計算上の困難さは、応力解放中のrDoFの「手動凍結」を必要とした。このシミュレーション技術により、利用可能なシミュレーション時間スケールでは扱いにくいそのような長い棒状の効果を模倣することができ、この特定のアプローチは望ましい領域を探索するための唯一実行可能な経路となった。

比較優位性

制御された回転自由度を持つ棒状の非球形不純物を組み込むという選択された方法は、球形不純物や熱アニーリングを含む以前のゴールドスタンダードに対して、圧倒的な質的優位性を示している。

球形不純物と非球形（二量体）不純物を比較すると、本論文は、球形包有物による機械的特性の「最小限の」改善を強調している。例えば、球形不純物 $c_s = 0.1$ を持つ 3DKA 系における降伏ひずみ ($\gamma_y$) は、純粋な系での $\gamma_y = 0.09$ から $\gamma_y = 0.107$ へとわずかに増加したにすぎない（図2(a,b)）。対照的に、$c_d = 0.1$ の非球形二量体不純物は、降伏ひずみを $\gamma_y = 0.127$ に達させ、「ほぼ100%の改善」をもたらした（図3(a,b)）。これは実質的な質的違いである。

構造的な利点は、非球形粒子によって導入される「追加の回転自由度（rDoF）」にある。これらのrDoFは、「内部応力を散逸するための追加の局所経路」を提供し、系が「追加の応力を解放し、より高い負荷を維持し、非局所的なせん断帯形成による破壊前に、より多くの応力を蓄積する」ことを可能にする。このメカニズムにより、系は球形不純物よりも大幅に高い負荷に耐え、降伏遷移をより効果的に遅延させることができる。球形不純物はこれらの回転自由度を欠いている。

さらに、シミュレーションでこれらのrDoFを凍結し、非常に長い棒状の効果を模倣する能力は、「超安定性様の機械的応答」と「極めて脆い」破壊（図5(e,f)）につながる。これは強力な構造制御メカニズムを示している。rDoFの存在は延性を向上させ、それらの制限または不在は脆性および超安定性を促進する。構造秩序パラメータ $\Theta$ は、不純物分率と棒状長さの増加とともに一貫して減少する（図2(d)、図4(e)）ため、構造安定性が向上していることを示しており、これは制御が不十分な方法よりも直接的な質的利点である。この方法は、不純物の縦横比と回転自由度を調整するだけで、材料の機械的応答を延性から高度に脆性まで、正確に調整可能な制御を提供する。

制約との整合性

選択された方法は、問題定義に概説された制約と完全に整合しており、厳しい要件と解決策のユニークな特性との間の明確な「結婚」を示している。

主な制約の一つは、「熱アニーリングが適用できない非ブラウン運動コロイド系を含む、より広範な実験系でアクセス可能な」延性-脆性遷移を調査するためのプロトコルが必要であるというものであった。非球形不純物、特に棒状のものの使用は、この制約に直接対処する。著者らは明確に、彼らのアプローチが「分子ガラスとコロイドガラス実験の両方で容易にアクセス可能な新しい延性-脆性遷移」を提供すると述べている。これにより、ソフトガラスを扱う実験家にとって、従来の熱アニーリング法が実行不可能な場合に、非常に適切で実用的な方法となる。

もう一つの暗黙の制約は、「固有無秩序強度」を正確に制御して臨界点をプローブするための方法が必要であるというものであった。非球形不純物の数分率と縦横比、およびそれらの回転自由度を制御することにより、著者らは系の構造安定性と機械的応答に対する微細な制御を得る。構造秩序パラメータ $\Theta$ は、この制御された無秩序の定量的な尺度として機能し、観測された機械的挙動との明確な相関を示す。

最後に、非常に長い棒状のシミュレーションという計算上の制約、「非常に遅いダイナミクス」と「ゼロに近い回転拡散率」を示すものは、応力解放プロセス中にrDoFを人工的に凍結することによって克服された。この巧妙な回避策により、著者らは「非常に長い棒状の効果を模倣」し、それ以外では利用可能なシミュレーション時間スケールでは不可能であった結果として生じる超安定性、脆性挙動を研究することができた。これは、計算上の障壁にもかかわらず、現実的な物理現象を研究する必要性との実用的な整合性を示している。

代替案の却下

本論文は、非晶質固体の延性-脆性遷移を調査するという特定の問題に対して、実験的にアクセス可能な方法で、特に多様な系に対して、それらの限界に基づいて、いくつかの代替アプローチを暗黙的かつ明示的に却下している。

1. 熱アニーリング: 著者らは、熱アニーリングが無秩序強度を変化させることができる一方で、「有限サイズ効果のために部分的な支持しか提供せず」、その結論は「決定的ではない」[15, 22]と指摘している。さらに重要なのは、「非ブラウン運動コロイド系には適用できない」ことであり、これらは実験的検証の主要な対象である。分子ガラスであっても、アニーリングを通じて「現実的な系」を達成することは困難であり、「利用可能な最低冷却速度でさえ、アニーリングが不十分な状態をもたらす」。スワップモンテカルロ[34, 35]のような手法は超安定ガラスを生成できるが、「実験には容易に適用できない」。これにより、熱アニーリングは普遍的なアプローチとしては不適切である。

2. 粒子固定化: 粒子固定化は最近の研究[28]でコロイドガラス実験[29]における無秩序制御に使用されているが、著者らはそれが「分子ガラスでは依然として困難」であると述べている。有望な技術ではあるが、提案された微量合金化戦略ほど、異なる材料タイプにわたる広範な適用性を提供しない。

3. 球形不純物: 本論文は、その非球形不純物アプローチを球形不純物と比較している。球形不純物は「より高い降伏ひずみと脆性の増加」をもたらすが、著者らは「改善は最小限」であると結論付けている（例：3DKAでは $\gamma_y$ が0.09から0.107に増加）。これは、延性-脆性遷移を堅牢に制御するための十分に効果的な方法として、球形不純物を却下する明確でデータに基づいた理由である。球形粒子における回転自由度の欠如は、非球形粒子と比較して、応力散逸と安定性向上能力を制限する。

本論文は、生成モデリングではなく材料特性の物理シミュレーションに焦点を当てているため、GANや拡散モデルのような機械学習アプローチについては議論していない。代替案の却下は、調査対象の現象に対して、十分な制御、実験的アクセス可能性、または望ましい機械的特性向上の大きさを提供できないことに根ざしている。

FIG. 6. Ultra-stability with frozen rDoF: Kinetic as- pect. (a) The potential energy per particle e is plotted for a system with 10% dimers subjected to heating-cooling cycles at the rate of dT/dt = 10−4 (same as the rate with which the sample was prepared). The solid lines are for a system with frozen rotational DoFs; the observed hysteresis during the first cycle indicates the ultra-stable character. The sec- ond heating cycle does not show any hysteresis. The dashed lines are for the system evolved with rotational DoFs, and the absence of hysteresis is seen as expected. (b) The specific heat CV = de/dT is calculated from data in (a) by numerical differentiation and conveys the same

数学的・論理的メカニズム

マスター方程式

本研究における粒子の挙動、ひいては非晶質固体の機械的特性を支配する基本的な相互作用は、レンナード・ジョーンズ（L-J）ポテンシャルによって記述される。このポテンシャルエネルギー関数は、系内の任意の2つの粒子がどのように相互作用するかを決定し、実行された分子動力学シミュレーションの基盤を形成する。これは、系の緩和を固有状態に導く際に最小化される目的関数である。

$$V_{\alpha\beta}(r) = 4\epsilon_{\alpha\beta} \left[ \left(\frac{\sigma_{\alpha\beta}}{r}\right)^{12} - \left(\frac{\sigma_{\alpha\beta}}{r}\right)^6 \right]$$

項ごとの解剖

この方程式を詳細に分析し、各成分の役割を理解しよう。

• $V_{\alpha\beta}(r)$:
  1) 数学的定義: この項は、距離 $r$ で分離された $\alpha$ 型と $\beta$ 型の2つの粒子間の相互作用ポテンシャルエネルギーを表す。
  2) 物理的/論理的役割: 粒子ペア間の相互作用に蓄えられるエネルギーを定量化する。このエネルギーはそれらの間の力を決定し、それがそれらの動き、配置、および材料全体の構造的および機械的応答を決定する。
  3) なぜ減算か: 括弧内の2つの項は、それぞれ異なる物理現象、すなわち反発と引力を表す。減算は、これらが反対の力であることを意味する。最初の項（反発）は粒子を押し離し、2番目の項（引力）は粒子を引き寄せる。

• $4\epsilon_{\alpha\beta}$:
  1) 数学的定義: 相互作用エネルギーのスケーリング係数であり、$\epsilon_{\alpha\beta}$ は $\alpha$ 型と $\beta$ 型の粒子間の相互作用のポテンシャル井戸の深さである。
  2) 物理的/論理的役割: この係数は相互作用の強度または大きさを設定する。より大きな $\epsilon_{\alpha\beta}$ は、より強い引力と反発力を意味し、より堅牢または緊密に結合した系につながる。それは相互作用のエネルギー尺度を定義する。
  3) なぜ乗算か: 全てのポテンシャル関数の全体的な強度を比例してスケーリングすることを保証するために、全体的な乗数として機能する。

• $\sigma_{\alpha\beta}$:
  1) 数学的定義: $\alpha$ 型と $\beta$ 型の粒子間のポテンシャルエネルギーがゼロとなる距離。それは実質的に粒子の「サイズ」または衝突直径を表す。
  2) 物理的/論理的役割: このパラメータは相互作用の長さスケールを定義する。それは粒子の実効的な「ハードコア」サイズを決定し、それらがどれだけ密にパッキングできるか、および引力が有意な範囲を決定する。それは材料の密度と局所構造を確立するために重要である。
  3) なぜ $r$ で割るのか: それは、ポテンシャルを相対的な分離の観点から表現できるように、粒子距離 $r$ を特徴的な粒子サイズ $\sigma_{\alpha\beta}$ で正規化し、括弧内の項を無次元化する。

• $r$:
  1) 数学的定義: 相互作用する2つの粒子の中心間の瞬間的な距離。
  2) 物理的/論理的役割: これは基本的な空間変数である。ポテンシャルエネルギー、したがって粒子間の力は、この分離に直接依存する。
  3) なぜ逆べき乗則 ($r^{-12}$ および $r^{-6}$) なのか: これらの逆べき乗則は、原子間および分子間力の典型的な性質をモデル化するために選択される。$r^{-12}$ 項は電子雲の重なりによる非常に短距離で強い反発を表し、$r^{-6}$ 項はより長距離で弱い引力（例：ファンデルワールス力）を表す。

• $\left(\frac{\sigma_{\alpha\beta}}{r}\right)^{12}$:
  1) 数学的定義: レンナード・ジョーンズポテンシャルの反発成分であり、$r^{-12}$ に比例する。
  2) 物理的/論理的役割: この項は、粒子が互いに重なったり、同じ物理空間を占めたりするのを防ぐ、強く短距離の反発をモデル化する。それは粒子の「ハードコア」性質を保証する。
  3) なぜ指数12なのか: この急峻なべき乗則は、パウリの排他原理を模倣し、粒子が非常に近くなったときの反発の急激な増加を効果的に近似する経験的な選択である。

• $\left(\frac{\sigma_{\alpha\beta}}{r}\right)^{6}$:
  1) 数学的定義: レンナード・ジョーンズポテンシャルの引力成分であり、$r^{-6}$ に比例する。
  2) 物理的/論理的役割: この項は、材料を結合させる原因となる、より長距離の引力（例：ファンデルワールス力）をモデル化する。
  3) なぜ指数6なのか: このより緩やかなべき乗則は、多くの原子および分子系で一般的な誘起双極子-双極子相互作用の特性である。

ステップバイステップの流れ

レンナード・ジョーンズポテンシャルは、非晶質固体の延性-脆性遷移を探求するために使用される分子動力学シミュレーションを駆動するエンジンである。抽象的なデータポイント（粒子の位置）がシミュレーションを流れる様子を以下に示す。

1. 系設定: 各粒子は、そのタイプ（$\alpha$ または $\beta$）と初期空間座標 $(X_i, Y_i, Z_i)$ によって定義される粒子の集まりを想像してください。これらのタイプは、その相互作用の特定の $\epsilon_{\alpha\beta}$ および $\sigma_{\alpha\beta}$ パラメータを決定します。
2. 高温平衡化: 系は高温に「加熱」される。各粒子の位置は、全ポテンシャルエネルギー（近傍粒子との全ての $V_{\alpha\beta}(r)$ 相互作用の合計）の勾配から導出される、それに作用する正味の力に基づいて連続的に更新される。粒子は激しく動き、エネルギーランドスケープを探索し、安定した液体状態が達成されるまで続く。
3. 冷却と急冷: 次に、粒子の運動エネルギーを徐々に減少させることによって、系は「冷却」される。温度が低下すると、$V_{\alpha\beta}(r)$ からの引力がより支配的になり、粒子は無秩序な非晶質固体状態に落ち着く。最後に、系はポテンシャルエネルギーを最小化することによって固有状態に「急冷」される。これは、粒子が全 $V_{\alpha\beta}(r)$ の合計が局所的最小値になる位置に移動され、全ての運動エネルギーが除去され、機械的に安定な構成に達することを意味する。
4. 非熱的準静的ひずみ（AQS）サイクル:
   • アフィン変換: 小さな増分せん断ひずみ $\delta\gamma$ が系全体に適用される。粒子 $(X_i, Y_i, Z_i)$ について、そのx座標はアフィン変換されて $X_i \leftarrow X_i + \delta\gamma Y_i$ となる。このステップは粒子の位置を近傍粒子に対して変位させ、一時的に系を機械的平衡から外す。
   • エネルギー最小化（緩和）: アフィンステップの後、系はもはやエネルギー最小値にない。次に、粒子は緩和できるようにされる。これには、共役勾配法を使用して、全ポテンシャルエネルギーを最小化するために位置を反復的に調整することが含まれる。力（$V_{\alpha\beta}(r)$ の勾配）は、各粒子を新しい局所的エネルギー最小値に導く。この緩和プロセスは、材料が応力下でどのように変形し、再配置するかをシミュレートする。
5. 反復変形: ステップ4aと4bが多くの小さなひずみ増分に対して繰り返される。系の応答（例：応力）は、各緩和ステップ後に記録される。この反復プロセスにより、研究者は材料の応力-ひずみ曲線を描き、塑性イベントと降伏挙動を観察できる。

最適化ダイナミクス

この文脈での「最適化」とは、主にエネルギー最小化を通じて、系が機械的に安定な構成を見つけるプロセスを指す。

• エネルギーランドスケープ: 系全体のポテンシャルエネルギー $U = \sum_{\text{all pairs } \alpha\beta} V_{\alpha\beta}(r)$ は、複雑で高次元のエネルギーランドスケープを定義する。このランドスケープは、多数の局所的最小値（安定な構成または「固有状態」を表す）とエネルギー障壁（鞍点）によって隔てられている。
• 勾配を駆動力として: 「学習」または「更新」メカニズムは、各粒子に作用する力、すなわちそれらの位置に対する全ポテンシャルエネルギーの負の勾配（$\mathbf{F}_i = -\nabla_i U$）に依存する。これらの力は、緩和中の粒子移動の方向と大きさを決定する。
• 共役勾配法: この反復アルゴリズムは、エネルギーランドスケープをナビゲートし、局所的最小値を見つけるために使用される。それは以下のように動作する。
  1. 勾配計算: 現在の粒子位置で、力（勾配）が計算される。
  2. 探索方向の決定: 以前の方向と「共役」な探索方向が選択され、最小値への効率的な降下を目指す。これはジグザグ運動を避け、収束を速めるのに役立つ。
  3. 線形探索: その特定の方向でエネルギー最小値を見つけるために、粒子はその方向に沿って移動される。
  4. 位置の更新: 粒子位置が更新され、プロセスが繰り返される。
• 状態更新と収束: 系は、全ての粒子に対する正味の力が所定の許容値以下になるまで、状態（粒子位置）を反復的に更新する。この時点で、系は局所的エネルギー最小値に収束し、機械的平衡を示す。この反復プロセスにより、系はアフィンひずみによって摂動された後に「緩和」し、安定な構成を見つけることができる。
• 塑性ときず: ひずみが蓄積するにつれて、系は最終的に、それが占めている局所的最小値が不安定になる点に達する。緩和するためには、エネルギー障壁を乗り越えて、異なる局所的最小値に遷移する必要がある。これらの遷移は「塑性イベント」であり、突然の不可逆的な粒子再配置に対応する。これらのイベントの集合的な挙動が、材料の降伏応答を形成する。延性降伏では、塑性イベントは多数かつ空間的に分布しており、段階的な応力応答につながる。脆性降伏では、塑性イベントはせん断帯に局在化し、壊滅的な不連続な応力降下につながる。L-Jポテンシャルによって決定されるエネルギーランドスケープの形状は、これらの再配置の障壁と経路を制御し、それによって材料の機械的挙動を制御する。不純物、特に非球形不純物の導入は、このランドスケープを変更し、これらの塑性イベントの容易さと性質を変化させる。

結果、限界、結論

実験設計とベースライン

著者らは、主に2次元および3次元のコブ・アンダーセン（KA）モデルガラス形成体を用いたレンナード・ジョーンズ（L-J）粒子の分子動力学シミュレーションを使用して、実験を細心の注意を払って設計した。中心的なアイデアは、非熱的準静的ひずみ（AQS）条件下での降伏遷移に及ぼす影響を観察するために、バイナリガラスマトリックスに様々な種類の不純物を導入することであった。

実験設定は、高温での平衡化による静穏状態の準備、それに続く冷却と固有状態への急冷を含んでいた。その後、増分的に変形が適用された。提案されたメカニズムがテストされた「犠牲者」またはベースラインモデルには以下が含まれる。
1. 純粋系: 不純物なしのガラス系 ($c_s = 0.00$, $c_d = 0.00$, $c_r = 0.00$)。これらは延性挙動の基本的な参照として機能した。
2. 球形不純物を含む系: より大きな直径 ($\sigma_s = 2.0 \sigma_{AA}$) を持つ第3の粒子タイプが追加され、$N_T = 100000$ 個の粒子からなる系で、その数分率 ($c_s \in [0, 0.1]$) が変化した。
3. 非球形二量体不純物を含む系: 2つの重なり合った球形ビーズ（直径 $\sigma_b = 2.0 \sigma_{AA}$、縦横比 $L_d : \sigma_b = 2.3 : 2$）で構成された棒状粒子が、数分率 ($c_d$) を変化させて導入された。これらは回転自由度（rDoF）を持つように設計された。
4. より長い棒状不純物を含む系: 複数のビーズ（直径 $\sigma_r = \sigma_{AA}$）でできた、より非球形な棒状粒子が、固定分率 ($c_r = 0.1$) で、その長さ ($L_r$) を変化させて導入され、回転自由度を制御した。
5. 人工的に凍結されたrDoFを持つ系: 回転の効果を分離するために、棒状不純物は、アフィン変換を許容しながらも、エネルギー最小化ステップ中に人工的に回転を妨げられた。これは重要な制御実験であった。
6. アニーリングが不十分なサンプル: 延性-脆性遷移の研究のために、サンプルは比較的高い温度 ($T=0.6$) で準備され、急速に冷却 ($T=10^{-1}$) され、アニーリングが不十分で固有に延性な状態を保証した。

機械的応答は、いくつかの主要な観測量を使用して厳密に特徴づけられた。
* 応力-ひずみ曲線 ($\sigma_{xy}$ vs. $\gamma$): 降伏ひずみ、せん断弾性率、および応力オーバーシュートまたは不連続降下の存在を直接観察するため。
* 非連結感受率 ($\chi_{dis}(\gamma)$): $\chi_{dis}(\gamma) = N_T (\langle \sigma_{xy}^2 \rangle - \langle \sigma_{xy} \rangle^2)$ と定義され、そのピークの大きさとしばしば性は、降伏遷移の脆性を示す。より大きく、より狭いピークは、より脆性な破壊を示唆する。
* 構造秩序パラメータ ($\Theta$): 静穏状態の構造安定性と無秩序を定量化するための局所的な尺度（式(4)）。低い $\Theta$ はより良い安定性を示す。
* 回転非相関関数 ($S_r(\gamma)$): 機械的負荷下での棒状不純物の回転緩和の程度を定量化するため。より速い減衰は、より独立した移動度を示唆する。
* 非アフィン変位 ($D_{min}^2$): 塑性イベントの空間分布とせん断帯の形成を視覚化および定量化するため。
* 粒子あたりのポテンシャルエネルギー ($e$) および比熱 ($C_v = de/dT$): 加熱-冷却サイクルで使用され、ガラスの運動安定性と超安定性を評価するため。

最後に、臨界点に関する数学的言明を断固として証明するために、著者らは有限サイズスケーリング（FSS）解析を採用した。彼らは、回転凍結棒状の系について、粒子総数 ($N_T = [25000, 50000, 100000, 200000]$) を変化させ、$\chi_{dis}$ ピーク高さ、$\langle \Delta \sigma_{max} \rangle$、およびその変動 ($\chi_d$) のスケーリングを系サイズで解析した。これにより、熱力学極限 ($N \to \infty$) への外挿と、ランダム場イジングモデル（RFIM）の普遍性クラスとの比較が可能になった。

証拠が証明すること

本論文で提示された証拠は、非晶質固体の機械的特性と延性-脆性遷移の調整における非球形不純物と回転自由度の役割について、説得力のある物語を提供している。

1. 球形不純物は脆性と安定性を向上させる:

   • 応力-ひずみ曲線（図2a、b）は、球形不純物を添加すると降伏ひずみとせん断弾性率が著しく増加することを示している。例えば、3DKAモデルはより顕著な応力オーバーシュートを示す。
   • $\chi_{dis}$ プロット（図2c）は、球形不純物分率 ($c_s$) の増加に伴い、ピークがより高いひずみにシフトし、著しくシャープになることを示しており、これはより脆性的な破壊への遷移の決定的な証拠である。
   • 構造秩序パラメータ $\Theta$（図2d）は、$c_s$ の増加とともに減少しており、これらの不純物がガラスマトリックスに構造安定性と秩序の向上をもたらすことを示している。これは、冷却速度を単に低減することによって達成されるよりも実質的な効果である。

2. rDoFを持つ非球形二量体は優れた機械的強化を提供する:

   • わずかに非球形な二量体不純物（図3a、b）の導入は、球形不純物と比較して、降伏ひずみとせん断弾性率のさらなる顕著な増加をもたらす。3DKAモデルの場合、$c_d = 0.1$ の二量体分率により、球形不純物で見られた18%の増加のほぼ2倍である、降伏ひずみが40%増加する。
   • 応力寄与の微視的解析（図3d）は、棒状-球形相互作用がより高いひずみで応力を維持し続けることを否定できない証拠を提供しており、球形-球形相互作用が早期に飽和するのとは対照的である。これは、棒状を含む領域が構造的に安定しており、より高い負荷に耐えられることを示唆している。
   • この優れた強化は、非球形粒子の「追加の回転自由度（rDoF）」に直接起因する。これらのrDoFは、「内部応力を散逸するための追加の局所経路」を提供し、系が「追加の応力を維持し、より高い負荷を維持し、壊滅的な破壊を引き起こす非局所的なせん断帯形成前に、より多くの応力を蓄積する」ことを可能にする。

3. rDoFの制限（より長い棒状）は脆性につながる:

   • 棒状不純物の縦横比が増加するにつれて（つまり、より長い棒状）、回転自由度は減少し、これは回転非相関関数 $S_r(\gamma)$ のより速い減衰（図4b）によって証明される。
   • このrDoFの減少は、降伏点を低い値に戻し、機械的応答はますます脆性になる（図4a）。
   • $\chi_{dis}$ ピーク（図4d）は、棒状長さの増加とともに、よりシャープで大きくなり、脆性相の出現を確認する。
   • 決定的なことに、非アフィン変位マップ ($D_{min}^2$)（図4f、g、h）は明確な遷移を示す。短い棒状は塑性イベントが空間的に広がることを可能にし（延性挙動）、一方、長い棒状は局所化された、系全体に広がるせん断帯の形成につながる（脆性破壊）。これは、中心的なメカニズム（rDoF）が破壊モードに直接影響することの強力な証拠である。

4. 人工的に凍結されたrDoFは超安定性と極端な脆性を誘発する:

   • 最も説得力のある証拠は、不純物のrDoFを人工的に凍結することから得られる。この操作は機械的応答を劇的に変化させ、不連続な応力降下を伴う極めて脆性な破壊につながる（図5a、c、e）。自由rDoFで観測された顕著な降伏ひずみ強化は完全に失われ、回転の重要な役割を再確認する。
   • $\chi_{dis}$ ピーク（図5b、d、f）は例外的にシャープで高くなり、高度に脆性な応答を示唆する。
   • 運動安定性解析（図6a、b）は、凍結rDoFを持つ系では加熱-冷却サイクル中にポテンシャルエネルギーに明確なヒステリシスを示し、これは超安定ガラスの特徴である。自由rDoFを持つ系は、そのようなヒステリシスを示さない。これは、rDoFを凍結することが機械的および運動的超安定性の両方につながることを証明している。

5. 有限無秩序臨界点が延性-脆性遷移を支配する:

   • アニーリングが不十分なサンプルで回転凍結棒状の分率 ($c_r$) を変化させることにより、著者らは延性から脆性降伏への系統的な遷移を実証する（図7a、b）。
   • 有限サイズスケーリング解析は、熱力学極限における有限無秩序臨界点の決定的な、否定できない証拠を提供する。
   • ドーピングされた系では $\chi_{dis}$ ピーク振幅は、純粋な系では小さく非応答性のままであるのに対し、系サイズとともにべき乗則 $N^{1.1 \pm 0.03}$ として増加する（図8b、c）。
   • 最大塑性降下 $\langle \Delta \sigma_{max} \rangle$ は、延性系 ($c_r=0$) では $N^{-0.4 \pm 0.02}$ としてゼロになるが、脆性系 ($c_r=0.1$) では有限で $N^{0.2 \pm 0.04}$ として増加する（図8d、e）。
   • $\langle \Delta \sigma_{max} \rangle$ の変動、すなわち $\chi_d$ は、系サイズとともに増加し、 $N^{0.36 \pm 0.027}$ として発散する非単調ピークを示し、その半値幅（FWHM）は $N^{-0.46 \pm 0.036}$ としてゼロになる（図8f、g）。
   • これらのべき乗則を用いたデータコラプス（図8h）は、遷移の臨界性を確認する。
   • スケーリング関係 $\chi_{dis}^{peak} \sim (\chi_{con}^{peak})^2$（図9a）は、この臨界性がランダム場イジングモデル（RFIM）普遍性クラスに属することを示唆している。
   • 熱力学極限 ($N \to \infty$) への $c_r^*(N)$ の外挿は、$c_r^* \sim 0.045 \pm 0.0003$ の有限臨界棒状分率を生成する（図9b）。この結果は、熱力学極限において延性-脆性遷移が存在しないという考えを否定し、真の無秩序駆動臨界遷移の堅牢な証拠を提供する。

限界と将来の方向性

この論文は、見事かつ徹底的な分析を提示しているが、あらゆる科学的努力と同様に、それは特定の制約の中で機能し、数多くの将来の探求の道を開く。

一つの顕著な限界は、回転拘束された不純物の高濃度における微視的メカニズムの不完全な理解である。著者らは、特に高不純物濃度において、せん断誘起ネマチック秩序が機械的応答に及ぼす影響は不明瞭なままであることを認めている。著者らは、初期配向が等方的で、回転拡散率が低い、またはrDoFが凍結された系では、これらの効果は最小限であると述べているが、この側面をより深く掘り下げることで、より微妙な相互作用が明らかになる可能性がある。

もう一つの制約は、非熱的準静的ひずみ（AQS）条件に主に焦点を当てていることである。著者ら自身が指摘するように、熱効果と有限ひずみ率は、特に回転凍結棒状の場合、延性から脆性への臨界遷移を変化させる可能性がある。彼らは、微量合金化された系が有限ひずみ率と温度でも機械的強度が増強されると予想しているが、これらの変数を組み込んだ包括的な研究は、完全な理解のために不可欠である。

さらに、本研究は、粒子形状（楕円形対棒状）と境界相互作用の性質（例：棒状を形成する球の重なりの程度）が機械的応答を制御する重要なパラメータであるが、これらの効果は将来の研究に残されていることを強調している。これは、現在のモデルにおける単純化を示唆しており、拡張される可能性がある。計算上の限界により、回転拡散率がゼロに近い非常に長い棒状のシミュレーションも不可能であり、代わりに「凍結rDoF」アプローチが必要となった。効果的ではあるが、そのような系を直接シミュレートすることは、さらなる洞察を提供する可能性がある。

将来に向けて、本論文の発見はエキサイティングな議論のトピックと研究方向性を提供する。

1. 不純物設計の完全なパラメータ空間の探求: 将来の研究では、様々な不純物形状（例：楕円体、分岐ポリマー）や、ホストマトリックスとの異なる相互作用ポテンシャルの影響を系統的に調査できる。これにより、不純物の特性が機械的特性と延性-脆性遷移にどのように影響するかについてのより包括的な理解が得られるだろう。
2. 熱効果とひずみ率の統合: 現在のAQS解析を有限温度とひずみ率に拡張することが不可欠である。これにより、より現実的な条件下での延性-脆性遷移の相図をマッピングでき、新しい臨界点やクロスオーバーが明らかになる可能性がある。熱アニーリング法（スワップモンテカルロなど）と棒状不純物ドーピングを組み合わせることは、固有無秩序、不純物誘起無秩序、および動的効果の相互作用に関する貴重な洞察を提供する可能性がある。
3. コロイド系および分子系における実験的検証: 本論文は、特にコロイド系におけるこの新しい延性-脆性遷移の実験的アクセス可能性を強調している。将来の研究は、非球形粒子でドープされたコロイドガラスにおける塑性イベント、せん断帯形成の観察、および機械的特性の測定に高度なイメージング技術を使用することに焦点を当てるべきである。これは、シミュレーションと現実世界の材料との間のギャップを埋めるだろう。従来の熱アニーリングが困難なソフトガラスの場合、この微量合金化戦略はゲームチェンジャーとなる可能性がある。
4. 超安定性とせん断帯の微視的メカニズム: 回転拘束された不純物の高濃度における遷移を駆動する微視的メカニズムのさらなる調査が求められる。これには、せん断誘起ネマチック秩序が機械的応答にどのように影響するか、および長い棒状から課される静的相関長が、過冷却領域における超安定ガラスの点-設定（PTS）相関長と比較してどのように見えるかを理解することが含まれる。これには、高度な理論モデリングと高解像度シミュレーションが必要になる可能性がある。
5. 材料設計への応用: 本研究の結果は、非晶質固体の機械的特性を調整するためのシンプルで制御された微量合金化戦略を提供する。将来の研究では、特定の産業用途（頑丈な複合材料から自己修復材料まで）のために、延性または脆性を調整した新しい材料を設計するために、これらの原則をどのように適用できるかを検討できる。凍結rDoFによる超安定性と脆性破壊を誘発する能力は、高い機械的安定性を必要とする用途に特に適している可能性がある。
6. 普遍性クラスの探求: さらに詳細な有限サイズスケーリング解析は、より多くの系サイズと異なる縦横比で実行される可能性があり、RFIM普遍性クラスとの関連を固め、異なる条件下で他の普遍性クラスが出現するかどうかを探求できる。これは、非平衡無秩序系における相転移の基本的な理解を深めるだろう。

FIG. 1. Components of the system. (a, b) Particle type ‘A’ and ‘B’ of the parent KA system. (c) Spherical impurity with twice the diameter of the particle ‘A’. It has a variable number fraction of cs in the system. (d) Slightly aspherical dimers composed of larger particles with σb = 2.0, added at number fraction cd in the system. (e) Rods with larger aspect ratio, made by attaching A-type particles, studied across dif- ferent aspect ratios at a fixed number fraction of cr = 0.1 FIG. 4. Mechanical properties with aspherical impurities: Effect of the length of rod impurity. (a) Stress-strain curves for 3dKA systems doped with 10% rods (Fig. 1(e)) of different lengths, Lr. The yield point shifts back as Lr increases, while the response becomes increasingly brittle due to restricted rotations. (b) Rotational de-correlation function for rods of different lengths with mechanical loading; the red points indicate the equal net D2 min . (c) Single ensemble stress-strain plots for systems with Lr = 2.5 (Lr = 1.3) show abrupt-brittle (continuous-ductile) yielding. (d) Corresponding susceptibility plots with increasing magnitude and sharpness indicate the sample’s emerging brittle behavior. (e) The structural order parameter decreases with increasing Lr, indicating greater structural stability. (f) D2 min, averaged in spatial-strips perpendicular to the shear band (at dx = 0) for systems with rods of varying Lr. Strain values are selected to yield a similar area under the curve, ensuring equal plasticity. Larger spread for smaller Lr and large displacement peak for larger Lr advocates the emerging brittle behavior. (g, h) D2 min maps obtained at equal net displacements for Lr = 1.3 and Lr = 2.5, respectively, illustrating the contrast between ductile and brittle responses FIG. 8. Impurity driven ductile-to-brittle transition, FSS study. (a) System size dependence of stress-strain curves for the 3dKA system prepared at a high temperature of T = 0.6. The pure system, shown in lighter colors, exhibits a ductile macroscopic response with no system size dependence, while systems with 10% rod impurities (aspect ratio 1.9 : 1) show progressively brittle behavior as system size increases. (b) Susceptibility plots for both the pure and doped systems at various system sizes. (c) System-size dependence of peak height for different impurity concentrations. For cr ≤0.02, the χp dis saturates, while it follows the indicated power law for cr ≥0.05. (d) System size dependence of ⟨σmax⟩with increasing cr. The response shifts from vanishing to not-vanishing in the thermodynamic limit on increasing cr. (e) The increase in ⟨σmax⟩with increasing cr becomes sharper with increasing system size. (f) Ensemble-level fluctuations of ⟨σmax⟩reveal a nonmonotonic peak structure, marking the transition point. As system size increases, the peak narrows, and its height increases, highlighting the critical nature of the transition. (g) The peak height diverges as χp d ∼N 0.36±0.027, while the FWHM vanishes as FWHM ∼N −0.46±0.036. (h) The obtained collapse of data points in panel (f) using the determined power laws