4. 結果

4.1 蒸着法と凝集法の直接比較(100分子系)

蒸着膜(左)と凝集膜(右)の最終構造スナップショット

4.1.1 密度の温度依存性

密度-温度曲線:蒸着法 vs 凝集法

表7: 密度-温度データ(100分子系)

$T^*$ 蒸着法 密度 ($\sigma^{-3}$) 凝集法 密度 ($\sigma^{-3}$) 差異 (%)
0.3 1.006 1.011 -0.5
0.5 1.027 0.990 +3.7
0.7 1.004 0.983 +2.1
0.9 0.994 0.959 +3.6
1.1 0.988 0.942 +4.9
1.3 0.968 0.914 +5.9
1.5 0.951 0.895 +6.3
1.7 0.925 0.875 +5.7
2.0 0.884 0.840 +5.2

蒸着法で形成された薄膜は、$T^ = 0.3$を除くほぼ全温度域で凝集法に対して高い密度を示した。特に中間温度域($T^ = 0.5 \sim 1.5$)で3-6%の密度優位性が顕著であり、$T^* = 1.5$で最大差6.3%を記録した。

$T^* = 0.3$での密度がほぼ同等(蒸着: 1.006, 凝集: 1.011)であることは興味深い。これは最低温度での密な充填には上限があり、両手法とも類似の局所構造に到達していることを示唆する。しかし、高温側では蒸着膜が顕著に密度を維持しており、より安定な構造であることが明らかである。

4.1.2 ポテンシャルエネルギーの温度依存性

ポテンシャルエネルギー-温度曲線

表8: 主要温度でのPE比較

項目 蒸着法 凝集法
PE ($T^*$=0.3) 10.15 10.64
PE ($T^*$=2.0) 13.12 13.74

蒸着法は全温度範囲でより低いPEを示した。$T^* = 0.3$での差(10.15 vs 10.64)は、蒸着膜がよりエネルギー的に有利なパッキング状態にあることを示す。

4.1.3 ガラス転移温度

ガラス転移温度の推定:密度微分と比熱

比体積の温度依存性

密度-温度曲線の傾き変化から推定した$T_g$は:

  • 蒸着法: $T_g$ ~ 1.8
  • 凝集法: $T_g$ ~ 1.3

蒸着法は凝集法に対して$T_g$が約38%向上した。これは蒸着膜が高温まで秩序構造を維持できる、より安定なガラス状態にあることを意味する。

4.1.4 動径分布関数(RDF)

動径分布関数(RDF)比較:蒸着法 vs 凝集法

蒸着膜のRDFは$r \sim 1.0\text{-}1.15\sigma$に鋭いピーク($g(r) \sim 32$)を示し、高度に秩序化されたπ-πスタッキング構造の形成を示唆した。一方、凝集膜のRDFはより広いピーク領域($g(r) \sim 3\text{-}5$)を示し、構造がより無秩序であることが確認された。

4.1.5 小括

100分子系の直接比較により、蒸着法は凝集法に対して(1)中間温度域で3-6%高い密度、(2)全温度域で低いPE、(3)約38%高い$T_g$、(4)より秩序化されたRDFパターンを示すことが定量的に確認された。これは実験で報告されている超安定ガラスの特徴と定性的に一致する。


4.2 蒸着法パラメトリック研究

4.2.1 蒸着速度($V_{dep}$)の影響

表9: 蒸着速度と薄膜特性

$V_{dep}$ 密度 ($T^*$=0.3) 密度 ($T^*$=2.0) PE ($T^*$=0.3) $T_g$推定値
0.5 0.969 0.830 11.44 1.5
1.0 0.686 0.679 10.41 0.9
3.0 1.018 0.890 10.28 1.7
5.0 0.983 0.810 11.29 1.7
10.0 0.884 0.883 10.92 1.5

蒸着速度は薄膜品質に劇的な影響を与えた。

  • $V_{dep} = 0.5$(過度に低速): 分子の運動エネルギーが不十分で基板到達が不安定。密度は0.969と中程度。
  • $V_{dep} = 1.0$(熱速度近傍): 最も劣る結果。分子が基板から散逸し、極めて低密度(0.686)の膜を形成。
  • $V_{dep} = 3.0$(最適値): 強い指向性運動により効率的に基板に到達。最高密度(1.018)かつ最低PE(10.28)を達成。
  • $V_{dep} = 5.0$: 良好だが$V=3.0$にやや劣る。高い運動エネルギーによる衝撃で表面構造にわずかな損傷。
  • $V_{dep} = 10.0$(過度に高速): 衝撃エネルギーが形成中の構造を破壊。密度低下(0.884)。

最適蒸着速度は熱速度の約3-5倍であり、これは十分な基板到達性と表面構造の保存のバランスに対応する。

4.2.2 基板温度($T_{dep}$)の影響

表10: 基板温度と薄膜特性

$T_{dep}$ 密度 ($T^*$=0.3) 密度 ($T^*$=2.0) PE ($T^*$=0.3) $T_g$推定値
0.3 0.874 0.847 11.13 1.7
0.5 1.026 0.904 10.20 1.7
0.7 1.018 0.890 10.28 1.7
0.9 1.005 0.912 10.68 1.5
1.1 0.967 0.883 11.00 0.7

基板温度は表面拡散の活性化を通じて薄膜品質を支配する。

  • $T_{dep} = 0.3$(過度に低温): 表面移動度が極端に低く、分子が到達位置で凍結。密度は0.874と低い。
  • $T_{dep} = 0.5$(最適値): 適度な表面移動度により最適パッキングを達成。最高密度(1.026)かつ最低PE(10.20)。
  • $T_{dep} = 0.7$: 0.5に匹敵する良好な結果。
  • $T_{dep} = 0.9$: 熱揺らぎが増大し、わずかに秩序低下。
  • $T_{dep} = 1.1$(過度に高温): $T_g$が0.7に急落。高温による過度な分子運動が構造安定性を破壊。

最適基板温度は$T_g$の0.5-0.7倍の範囲にあり、これは実験で報告されている0.85$T_g$付近の最適条件と定性的に一致する。

4.2.3 基板相互作用($\varepsilon_{wall}$)の影響

表11: 基板相互作用強度と薄膜特性

$\varepsilon_{wall}$ 密度 ($T^*$=0.3) 密度 ($T^*$=2.0) PE ($T^*$=0.3) $T_g$推定値
0.5 0.776 0.833 11.00 1.7
1.5 1.030 0.902 10.11 1.5
2.5 1.018 0.890 10.28 1.7
4.0 1.011 0.878 10.49 1.5

基板相互作用は分子の接着と配向に影響する。最適範囲は$\varepsilon_{wall} = 1.5 \sim 2.5$であり、「吸着安定性」と「表面移動度」のトレードオフの最適点に対応する。

4.2.4 成膜過程の解析

表12: 成膜段階と膜構造

分子数 ($N_{mol}$) 総原子数 $z_{max}$ ($\sigma$) $z_{mean}$ ($\sigma$) 表面被覆率 (%)
5 35 2.0 1.26 15.3
10 70 2.0 1.23 32.6
20 140 14.0 2.49 36.1
50 350 19.1 4.62 43.1

成膜過程は明確に4つのステージに分かれる:

  1. 核形成(N=5): 基板上に孤立した単分子層アイランドが形成。被覆率15%。
  2. アイランド成長(N=10): 複数のアイランドが成長・接近。被覆率33%。z方向にはまだ単層。
  3. 合体(N=20): アイランドの合体が始まり、多層膜が出現。被覆率は36%で頭打ち傾向。
  4. 連続膜形成(N=50): 準連続膜が形成。平均膜厚$z_{mean} \sim 4.6\sigma$。

4.3 凝集法パラメトリック研究

4.3.1 冷却速度の影響

表13: 冷却速度と薄膜特性

急冷ステップ数 冷却速度 (dT/step) 密度 ($T^*$=0.3) 密度 ($T^*$=2.0) $T_g$推定値
5,000 $1.2 \times 10^{-4}$ 0.915 0.883 0.7
10,000 $6.0 \times 10^{-5}$ 0.917 0.879 0.7
20,000 $3.0 \times 10^{-5}$ 0.913 0.862 1.7
50,000 $1.2 \times 10^{-5}$ 0.919 0.873 0.7
100,000 $6.0 \times 10^{-6}$ 0.954 0.863 1.7

最遅冷却(100,000ステップ)は急速冷却(5,000ステップ)に対して約4%の密度向上を達成した。ただし、100,000ステップでも蒸着法の最適条件には及ばない。

4.3.2 圧縮圧力の影響

表14: 圧縮圧力と薄膜特性

$P_{final}$ 密度 ($T^*$=0.3) 密度 ($T^*$=2.0) $T_g$推定値
0.1 0.998 0.827 1.7
0.5 0.989 0.865 0.7
1.0 0.913 0.862 1.7
2.0 1.017 0.874 1.3
5.0 1.059 0.966 1.3

P = 5.0では全手法中最高の低温密度(1.059)を達成。HDA(High-Density Amorphous)状態に相当する圧力緻密化ガラスを形成した。

4.3.3 C-C相互作用強度(π-πスタッキング)の影響

表15: C-C相互作用強度と薄膜特性

$\varepsilon_{CC}$ 密度 ($T^*$=0.3) PE ($T^*$=0.3) 密度 ($T^*$=2.0) 密度降下率 (%)
1.0 0.983 13.4 0.722 26.5
2.0 0.913 10.7 0.862 5.6
3.0 0.997 8.4 0.860 13.7
4.0 0.960 5.7 0.882 8.1

$\varepsilon_{CC} \geq 3.0$が有機半導体の熱的安定性に必要であり、π-πスタッキングの強化が膜安定化の最も直接的な材料戦略である。

4.3.4 主鎖剛性($K_{CCC}$)の影響

表16: 主鎖剛性と薄膜特性

$K_{CCC}$ 柔軟性 密度 ($T^*$=0.3) 密度 ($T^*$=2.0) $T_g$推定値
10 非常に柔軟 1.082 0.874 1.7
30 やや柔軟 1.054 0.914 1.5
50 デフォルト 0.913 0.862 1.7
100 非常に剛直 1.067 0.927 1.7

主鎖剛性は予想外のバイモーダル(二峰性)依存性を示した。中間剛性を避け、「非常に柔軟(折り畳み充填)」か「非常に剛直(配向充填)」の極端な設計が有利。

4.3.5 アニーリング効果

表17: アニーリング条件と薄膜特性

条件 密度 ($T^*$=0.3) PE ($T^*$=0.3) $T_g$推定値
急冷(10kステップ) 0.917 10.58 0.7
急冷(100kステップ) 0.919 10.62 0.7
アニール($T^*$=1.5) 1.037 10.31 1.5
アニール($T^*$=2.0) 1.025 10.42 1.7

アニーリングは凝集法における最も効果的な改善策(急冷比+12%密度向上)。


4.4 手法間の総合比較

表18: 全手法・全条件の横断的比較

手法 条件 密度 ($T^*$=0.3) PE ($T^*$=0.3) $T_g$推定値
蒸着 最適条件 ($T_{dep}$=0.5) 1.026 10.20 1.7
蒸着 最適条件 ($\varepsilon_{wall}$=1.5) 1.030 10.11 1.5
凝集 急速急冷 (5kステップ) 0.915 0.7
凝集 低速急冷 (100kステップ) 0.954 1.7
凝集 アニール (T=1.5) 1.037 10.31 1.5
凝集 高圧 (P=5.0) 1.059 1.3
凝集 柔軟分子 (K=10) 1.082 1.7
凝集 剛直分子 (K=100) 1.067 1.7
凝集 強スタッキング ($\varepsilon$=4.0) 0.960 5.7