ピットパターンSi(001)基板の表面拡散駆動形態進化：相場モデル研究

概要

ヘテロエピタキシャル島の側方秩序化は、堆積前の基板ピットパターン形成で常套化している。しかし、ピットの形状、配向、深さの制御は自明ではない。アニーリングや後続成長時にメタスタブル幾何構造が進化するため。本研究では、初期深さと形状の関数としてSi(001)上の典型的ピット形態を予測する連続体相場モデルを提示。モデルは表面拡散とシリコンのウルフ形状を再現する強異方性表面エネルギー密度を組み込む。実験で観察されるメタスタブル構成を再現し、初期深さによらずピット平坦化への共通運動論的経路を明らかにした。さらに、異なる平衡ファセットを持つ二材料ヘテロエピタキシーへモデル拡張し、例えばGe堆積時のピット輪郭回転を再現、最近の実験観測と合致。

背景

Ge/SiやInGaAs/GaAsなど格子ミスマッチ半導体ヘテロエピタキシーでは、ストランスキー・クラスタノフ（SK）成長モードで三次元島が形成される。自己組織化で量子ドットが得られるが、ランダム核生成固有のサイズ・形状分布がデバイス性能を制限。過去数十年、秩序構造へ成長を誘導する様々な手法が開発された。中でもピットパターン基板は、島の高い側方秩序化と精密なサイズ制御を両立する最も汎用的なアプローチの一つとして浮上。ピットパターンはナノインプリントリソグラフィ、電子ビームリソグラフィ+反応性イオンエッチング/ウェットエッチング、ナノインデンテーションなどトップダウン法で作製。ナノメートル精度の有序ピットアレイ設計が可能で、最適成長条件下でほぼ完全な島の側方秩序化を実現。ピット形状は全表面エネルギーに影響し、結果として島核生成を制御するため、ピット形態の精密制御が必須。高温では毛管力がピットをエネルギー的に有利な構成へ駆動、完全治癒に至ることも。そのためアニーリングや追加材料堆積でメタスタブルピット形状を安定化。安定化後もヘテロエピタキシャル成長中に更なる進化が起こり得る。本研究では、表面拡散による表面エネルギー低減駆動でピットパターン基板が進化する様相を記述。長さと時間スケールで実験に関連する拡散律速ヘテロエピタキシャル系に適用済みの相場定式化を採用。現実的な異方性表面エネルギーで平衡形状を再現可能、実験データで検証済み。ピットパターンSi(001)表面に焦点を当て、文献で広く研究されている系。論文構成：次節で相場モデル導入、第3節でSi(001)ピット平滑化と平衡への運動論的経路、第4節でGe堆積によるピット回転への応用、最後に知見の意義を議論。

実験方法

相場モデル

連続秩序パラメータφを導入、φ=1（固体）からφ=0（真空）へ滑らかに変化。φの進化は、強異方性表面エネルギーを尊重しながら表面拡散を捉ける退化カーン・ヒラード方程式で制御。全自由エネルギー汎関数：F = ∫Ω γ(n̂)[ε/2|∇φ|² + 1/ε B(φ)] d𝑟 + ∫Ω β/(2ε)[-ε∇²φ + 1/ε B'(φ)]² d𝑛。Ωはシミュレーション領域、εは界面厚、B(φ)=18φ²(1-φ)²は二重井戸ポテンシャル、βはウィルモア正則化パラメータ（強異方性領域で鋭いコーナーを丸める）。化学ポテンシャルμ=δF/δφが∂φ/∂t = D∇·[M(φ)∇μ]で進化を駆動。Dは拡散係数、M(φ)=36/ε φ²(1-φ)²は表面移動度。μの式には異方性表面エネルギー密度γ(n̂)とウィルモア正則化項を含み、ファセット進化の物理的現実的記述を保証。

異方性表面エネルギー

表面エネルギー密度をγ(n̂) = γ₀[1 - Σᵢ αᵢ(n̂·mᵢ)^{wᵢ} Θ(n̂·mᵢ)]でモデル化。ベクトルmᵢは最小エネルギー結晶方位（シリコンでは〈001〉, 〈113〉, 〈110〉, 〈111〉）に対応。各最小値の深さαᵢは文献報告の相対表面エネルギー値γᵢから導出、幅パラメータwᵢがファセット角度範囲を制御。シリコンでγ₀=1、α_{〈001〉}=0.15、w値50（〈113〉のみ100）。ウィルモア正則化パラメータβ=0.005で丸いコーナー分解能と計算可能性を両立。

初期形態とシミュレーション設定

(001)面上の初期ピット幾何構造を、半径L、深さHの円形窪みで定義、周囲平面へガウシアンテール（幅σ=L/2）で滑らかに接続。界面からの符号付き距離d(r)をハイパーボリックタンジェントプロファイルφ(r)=½[1-tanh(3d(r)/ε)]に代入し初期条件生成。典型的領域サイズL=1（任意単位）、界面厚ε=0.2。有限要素ツールボックスAMDiSを用い、半陰的時間積分スキームと界面メッシュリファインメント。側方に周期境界条件、垂直方向に無流束条件。拡散係数D=1とし、シミュレーション時間を任意単位で報告。

結果と考察

Si(001)ピットの平滑化

Fig.2でアスペクト比R=0.25とR=0.5のピット初期進化を示す。当初無特徴プロファイルが自発的にシリコン表面エネルギー密度の極小に対応するファセットを形成。浅ピット（R=0.25）では正方形〈001〉底面を〈113〉エッジと狭い〈110〉三角形が囲む。深ピット（R=0.5）では底部付近に追加〈111〉ファセット出現、〈110〉ファセット拡大。これら予測はピット形状Si(001)表面の実験観測と定量的に一致。Fig.3で深ピット（R=0.5）の長時間進化を追跡。透視図と俯瞰図の系列で、急峻な〈111〉ファセット消失、隣接〈113〉ファセット粗大化、〈110〉ファセット消滅、最終的な平坦化を示す。総表面エネルギー単調減少（Fig.4）、全R≤0.5で共通運動論的経路。数十nmピットの1100–1200°Cでの特性時間スケールは数時間オーダー、文献値のシリコン表面拡散と整合。

Ge堆積による形状変化の模倣

異なる平衡ファセットセットを持つ材料堆積の影響を評価するため、R≈0.1のピットパターンSi(001)基板上へのGe堆積をシミュレーション。t≈5.0の形態（Fig.3）から開始、シリコン表面エネルギー密度をGe/Si(001)系特有の〈105〉方向深い極小を含むものに置換。幅パラメータw_{〈105〉}=500で浅い〈105〉ファセットを捕捉。Fig.5で表面が急速再編成：残存〈113〉ファセット間に〈105〉ファセット核生成、後者犠牲に拡大。最終ピットは〈105〉ファセットのみで囲まれ、俯瞰図で初期配向から45°回転したピット輪郭を示す。この回転はGe堆積時に実験報告されており、純粋に表面エネルギー最小化から生じることを示した。

結論

ピットパターンSi(001)基板の表面拡散駆動進化を捉える連続体相場フレームワークを開発。シリコン表面エネルギーの強異方性を明示的に組み込み、メタスタブルファセット形成、ピット最終平坦化、これら状態へ至る運動論的経路を再現。表面エネルギー記述を〈105〉ファセット含むよう拡張し、Ge堆積時のピット輪郭回転が表面エネルギー低減のみで説明可能ことを実証。モデルは汎用的で、異方性パラメータ再設定で他基板へ容易適用可能、ナノメートル精度でヘテロエピタキシャル島配置を誘導するピット幾何設計の予測ツールを提供。