380 likes | 698 Views
GaAs(110) 表面の形成機構と緩和過程の動的モンテカルロシミュレーション. 2006 年 2 月 17 日 鳥取大学応用数理工学科 逢坂研究室 B02T8009B 小田泰丈. 研究背景. 近年半導体デバイスの微細化・高集積化を求めるに当たって、トップダウン型であった微細化も限界に近づいており、ボトムアップ型を用いて半導体デバイスを作成することが求められている。. GaAs(110). 結晶成長が困難とされてきた GaAs(110) 面は近年の技術革新により成長メカニズムが明らかになりつつある。
E N D
GaAs(110)表面の形成機構と緩和過程の動的モンテカルロシミュレーションGaAs(110)表面の形成機構と緩和過程の動的モンテカルロシミュレーション 2006年2月17日 鳥取大学応用数理工学科 逢坂研究室 B02T8009B小田泰丈
研究背景 • 近年半導体デバイスの微細化・高集積化を求めるに当たって、トップダウン型であった微細化も限界に近づいており、ボトムアップ型を用いて半導体デバイスを作成することが求められている。
GaAs(110) • 結晶成長が困難とされてきたGaAs(110)面は近年の技術革新により成長メカニズムが明らかになりつつある。 • 東京大学物性研の秋山英文氏らの実験グループによってGaAs(110)面上に他の面では見られない、3つの奇妙なパターンを発見された[1]。 [1]Surface step-edge kinetics control for an atomically smooth quantum well Mashahiro Yosita,Hidefumi Akiyama,Loren Pfeiffer,Ken W.West
3種類の奇妙なパターン Arrowhead 矢の先のような形。Fishと逆向きで原子2層分窪んでいる。 Fish 魚の形をしている。原子1層分だけ窪んでいる。 Boat 船の形をしている島。原子1層分だけ表面に突き出している。
量子ドットの形成メカニズム • 実験では格子不整合が重要だといわれていた。それに対してシミュレーションによって表面上での原子の動き方も重要であるということが報告されている。
研究目的 • 大規模系の計算が可能である動的モンテカルロ法を用いてナノ構造形成におけるバリアエネルギーの効果を調べる。 • 本研究では窪みのナノ構造形成に着目する。
理論的アプローチ • 表面のバリアエネルギーについては、第一原理計算によっていくつかの場面に対して明らかにされている[2][3]。 [2] 高橋秀昭 鳥取大学工学部応用数理工学科卒業論文 2003 [3] 小池規如 鳥取大学工学部応用数理工学科卒業論文 2004
第一原理計算の結果 それぞれのステップ付近の安定位置から,隣接する安定位置へ移動する際に超えるべき壁の高さは次のようになる. 平坦面での数値は次のように求まっている.
動的モンテカルロ法 • 動的モンテカルロ法とは、系で起こりうるすべてのプロセスに対して確率速度を与え、乱数を用いてプロセスを動かし、シミュレーションを進めていく手法である。
原子の移動確率速度 E =n1 ×ED1+n2×ED2
モデル化 • ナノ構造形成とバリアエネルギーの関係に着目しているので、GaとAsの原子を一成分としてシミュレーションを行った。 • 一成分で近似するにあたり、各局面で一番強く効いていると思われるパラメータを用いた。
一成分化 バリアエネルギーの効果が知りたかったので以下のような規格化を行った。 枠で囲まれた2つの原子を1つの原子として扱う Aサイトを基準として Bサイトを第一近接 Cサイトを第二近接として扱っている
[001]方向のステップでの処理 • バリアエネルギーが高すぎるため、ステップの上段には移動できないとしている。 • バリアエネルギーが高すぎるため、ステップから離れることができないとしている。 • ステップに沿う方向には移動できる。
[001]方向のステップでの処理 • バリアエネルギーが高すぎるため、ステップの上段には移動できないとしている。 • バリアエネルギーが高くなっているが、低い確率でステップから離れることができるとしている。 • ステップに沿う方向には移動できる。
計算条件 • 平坦面での移動の異方性は第一原理計算によって求められたGa原子の値を用いる。 • [001][001]方向のステップのステップに沿う方向の移動の異方性は、第一原理計算によって求められた値を用いる。 • 残りのステップ近傍での動きは窪みを大きく壊さない程度に設定している。 • 初期表面として正方形の窪みを与える。
計算条件 • 第一原理計算で求まっているパラメータでFishができるかを確認する。 • 第一原理計算で得られていない、バリアエネルギー、結合エネルギーの値をFishができる条件になるよう変えていく。
計算モデル1第一原理計算で得られた移動の異方性の効果についての検証計算モデル1第一原理計算で得られた移動の異方性の効果についての検証 計算条件 初期表面:20×20の窪み 温度:873K 第一近接の結合エネルギー:0.45eV 第二近接の結合エネルギー:0.17eV
結果 Fish構造が形成されることが確認できた
計算モデル2温度の効果 計算条件 初期表面:20×20の窪み 温度:923K 第一近接の結合エネルギー:0.45eV 第二近接の結合エネルギー:0.17eV
考察 • 第一原理計算で求められたパラメータでFishを再現することができた。 • 第一原理計算で得られていないパラメータも決定することができた。 • 温度の違いによってFish構造は崩壊していないことが確認できた。
計算モデル3平坦面での異方性の効果 • 第一原理計算で得られた平坦面でのバリアエネルギーの値を変更して、バリアエネルギーの効果について検証した。
平坦面での異方性(1) 873K異方性逆 923K異方性逆
平坦面での異方性(2) 873K異方性無し 923K異方性なし
考察 • 異方性を用いなかったものと、異方性を逆にしたもの、共にFishを再現することができなかった。 • Fish形成のために平坦面の異方性は重要な役割を果たしていることがわかった。
計算モデル4 • [001]方向のステップ近傍の動きの違いとFish形成の関係について調べた。 赤色の部分のステップから平坦面へ出ていくためのバリアエネルギーの値を変えて計算を行った。
[001]方向のステップから離れるためのバリアエネルギーが高すぎる例[001]方向のステップから離れるためのバリアエネルギーが高すぎる例
溝になる理由 バリアが強すぎるため赤線の部分が外れず、黄色の部分のみが崩壊してしまい最終的に溝になってしまう。
まとめ • 第一原理計算で求められたパラメータでFishを再現することができた。 • 第一原理計算で得られていないパラメータも決定することができた。 • 温度の違いによってFish構造は崩壊していないことが確認できた。 • Fish形成のためには平坦面での異方性、ステップ近傍での異方性共に重要である
今後の展望 • FishだけでなくBoatやArrowheadも再現できるようにする
GaAs(110)の模式図 GaAs(110)面の模式図. 2種類の原子が交互に,平坦に,ジグザグに並ぶ.
これまでの理論研究の成果(1) • 同研究室の卒業生である高橋隆亮によって平坦面でのGa,As原子のバリアエネルギーについて求められている[2]。 • 同研究室の卒業生である小池規如によって[00(-1)][001]方向のステップ近傍でのGa,As原子バリアエネルギーについても求められている[3]。