6. 結論

6.1 蒸着法の優位性

蒸着法は凝集法(標準急冷)に対して、中間温度域で3-6%高い密度、約38%高い$T_g$(1.8 vs 1.3)、より低いPE、より秩序化されたRDFを示した。蒸着法は「逐次的表面最適化」メカニズムにより、バルク加熱なしにアニーリング処理と同等の高密度ガラスを形成する。

6.2 最適プロセス条件

蒸着法の最適条件: - 蒸着速度: $V_{dep} = 3.0$(熱速度の約3-5倍) - 基板温度: $T_{dep} = 0.5$($T_g$の約0.5-0.7倍) - 基板相互作用: $\varepsilon_{wall} = 1.5\text{-}2.5$

凝集法の改善策(効果順): 1. アニーリング($T_g$直上、T=1.5): +12%密度向上 2. 高圧処理(P=5.0): HDA形成、+16% 3. 冷却速度低減(100kステップ): +4%

6.3 分子設計指針

  • π-πスタッキング: $\varepsilon_{CC} \geq 3.0$が熱的安定性に必須。共役系の拡張が有効。
  • 主鎖剛性: バイモーダル依存性を示し、非常に柔軟(フォールディング充填)または非常に剛直(配向充填)の極端な設計が有利。中間的な剛性は最も不利。

6.4 今後の課題

  1. より大規模な系(1000分子以上)での検証と有限サイズ効果の評価
  2. 全原子モデルとの比較による粗視化モデルの定量的妥当性検証
  3. 蒸着膜の異方的構造(面内 vs 面外配向)の詳細解析
  4. 電荷輸送特性(移動度、伝導度)との相関解析
  5. 多成分系(ホスト-ゲスト系)への拡張

付録

A. シミュレーション環境

  • ソフトウェア: LAMMPS(Python LAMMPS インターフェース使用)
  • 単位系: LJ還元単位($\varepsilon$, $\sigma$, mass = 1.0基準)
  • 時間刻み: dt = 0.001
  • 乱数シード: 42

B. 全パラメトリック研究の結果画像一覧

各パラメトリック研究で以下のプロットが生成されている:

  • density_vs_T.png — 全パラメータ値の密度-温度曲線オーバーレイ
  • pe_vs_T.png — ポテンシャルエネルギー-温度曲線
  • tg_analysis.png — $d(\text{density})/dT$による$T_g$推定

主要比較で追加生成: - snapshot_comparison.png — 最終構造のスナップショット比較 - rdf_comparison.png — 動径分布関数比較 - specific_volume_vs_T.png — 比体積の温度応答

C. 用語集

用語 説明
PVD Physical Vapor Deposition(物理蒸着法)
CG Coarse-Grained(粗視化)
LJ Lennard-Jones(レナード・ジョーンズ)
FENE Finitely Extensible Nonlinear Elastic
NPT 定温定圧アンサンブル
NVT 定温定積アンサンブル
NVE ミクロカノニカルアンサンブル
$T_g$ ガラス転移温度
PE ポテンシャルエネルギー
RDF 動径分布関数 (Radial Distribution Function)
HDA High-Density Amorphous(高密度アモルファス)
OPLS Optimized Potentials for Liquid Simulations

本研究は、ローカル環境で動作するLAMMPS分子動力学シミュレーションを用いた計算研究であり、再現可能な研究設計を特徴とする。