[論文レビュー] Asymmetry in crystal facet dynamics of homoepitaxy by a continuum model

Jian‐Guo Liu, Jianfeng Lu+2 more|arXiv (Cornell University)|Mar 12, 2017

nanoparticles nucleation surface interactions被引用数 3

ひとこと要約

本稿では、非線形表面拡散力学に基づく指数的移動度を有する退化型放物型PDEを用いた連続体モデルを提案し、ヘテロエpitaxial結晶表面の非対称な側面ダイナミクスを記述する。このモデルは、原子論的破断結合理論とミクロスケールのステップ拡散理論に基づいて導出され、解析的に、凸型（上側）の側面がサイズにかかわらず急速に拡大する一方、凹型（下側）の側面は臨界長を超過する場合にのみ進化することを示しており、非線形表面拡散力学に起因する凸型と凹型側面進化の根本的な非対称性を明らかにする。

ABSTRACT

In the absence of external material deposition, crystal surfaces usually relax to become flat by decreasing their free energy. We study an asymmetry in the relaxation of macroscopic plateaus, facets, of a periodic surface corrugation in 1+1 dimensions via a continuum model below the roughening transition temperature. The model invokes a highly degenerate parabolic partial differential equation (PDE) for surface diffusion, which is related to the weighted-$H^{-1}$ (nonlinear) gradient flow of a convex, singular surface free energy in homoepitaxy. The PDE is motivated both by an atomistic broken-bond model and a mesoscale model for steps. By constructing an explicit solution to the PDE, we demonstrate the lack of symmetry in the evolution of top and bottom facets in periodic surface profiles. Our explicit, analytical solution is compared to numerical simulations of the PDE via a regularized surface free energy.

研究の動機と目的

• ヘテロエpitaxyにおける結晶表面緩和の過程で生じる側面進化の非対称性の起源を理解すること。
• ステップの化学ポテンシャルに指数的依存する吸着原子フラックスの非線形的・指数的依存性を捉える連続体モデルを構築すること。
• 同じ物理的条件下でも凸型（上側）と凹型（下側）の側面が異なるように進化することを解析的に示すこと。
• キネティックモンテカルロシミュレーションやSiグレーティングの実験的アニールで観察された直感に反する非対称性の理論的説明を提供すること。
• 今後、2+1次元への拡張および弾性相互作用の組み込みを含めたモデルの発展の基盤を築くこと。

提案手法

• 凸型で特異な表面自由エネルギーを持つ非線形的・重み付き $H^{-1}$ 勾配流れとして特異拡散PDEを定式化する。
• 部分勾配形式を用いて、1+1次元周期的表面粗さにおける高さプロファイルの明示的解析解を構築する。
• アレニウス則とギブズ＝トムソン則に基づいて導出された指数的移動度関数を組み込み、強いステップ剛性と非線形吸着原子フラックスを反映する。
• 特異性を扱うために正則化された表面自由エネルギーを用いた数値シミュレーションと解析解を検証する。
• 連続体極限におけるバーバー＝キャブレラ＝フラスク（BCF）モデルと原子論的破断結合モデルに依拠して、PDEの物理的根拠を裏付ける。
• 周期的表面プロファイルにフレームワークを適用し、上側と下側の側面のダイナミクスを分離して分析する。

実験結果

リサーチクエスチョン

• RQ1 対称的な初期条件のもとで、周期的結晶表面粗さにおける上側と下側の側面が異なる速度で進化するのはなぜか？
• RQ2 ステップの化学ポテンシャルに依存する吸着原子フラックスの指数的依存性が、側面ダイナミクスの非対称性をどのように引き起こすのか？
• RQ3 キネティックモンテカルロシミュレーションで観察された非対称性は、非線形表面拡散を有する連続体モデルで説明可能か？
• RQ4 下側側面の運動を開始する臨界サイズ閾値は何か？ そして、それはPDE構造からどのように生じるのか？
• RQ5 化学ポテンシャル関係式の線形化を用いない場合、標準的なフィック型拡散モデルと比較して進化にどのような影響を与えるか？

主な発見

• PDEにおける指数的移動度のおかげで、上側の側面は初期サイズにかかわらず急速かつ連続的に拡大する。
• 下側の側面は、初期長さが臨界閾値を超過する場合にのみ進化する。この閾値はPDEの非線形構造から自然に生じる。
• 非対称性は直接的に指数的移動度項に起因し、フラックスがステップの化学ポテンシャルの指数関数に結合していることに起因する。
• 解析解は、特に直感に反する非対称性がキネティックモンテカルロシミュレーションで観察された定性的な挙動を正しく捉えていることを確認した。
• モデルはBCF理論の連続体極限および原子論的破断結合モデルと整合しており、物理的根拠が裏付けられている。
• 正則化された自由エネルギーを用いた数値シミュレーションは、解析予測を確認し、凸型と凹型領域における顕著な異なるダイナミクスを示した。