[논문 리뷰] Modeling and Simulation of Strained Heteroepitaxial Growth

Florian Much|arXiv (Cornell University)|2004. 10. 18.

nanoparticles nucleation surface interactions참고 문헌 69인용 수 7

한 줄 요약

이 박사학위논문은 기계적 스트레인을 가진 이종외향성 성장 현상을 모의하기 위해 래티스 자유형 운동역학 몽테카를로(KMC) 시뮬레이션 방법을 개발한다. 주로 격자 불일치 효과에 초점을 맞추며, 분자 정적 계산을 통해 확산 장벽을 산정하고 이를 KMC에 적용함으로써, 스트레인 완화는 불완전성 형성, Stranski-Krastanov 섬 형성, 표면 합금화를 통해 이루어지며, 임계 두께는 격자 불일치에 대해 거듭제곱 법칙으로 스케일링됨을 규명하였다. 이는 ZnSe/GaAs 및 InAs/GaAs 체계에서 실험 데이터와 뛰어난 일치를 보였다.

ABSTRACT

We introduce an off-lattice model with continuous particle distances and pair-potential interactions which allows for the efficient simulation of strained heteroepitaxial growth by means of kinetic Monte Carlo (KMC) simulations. We discuss in detail the application of our model to the simulation of three important strain relaxation mechanisms: First, the formation of misfit dislocations in strained adsorbate films is investigated. We find a power-law misfit dependence of the critical layer thickness, in agreement with various experimental studies. Second, we study the self-assembled island formation on a thin adsorbate wetting-layer (Stranski-Krastanov growth). We show that an increased diffusion barrier on the substrate is a possible kinetic reason for the formation of a stable wetting-layer. The formation of islands can be ascribed to a partial relaxation of adsorbate material which affects the diffusion barriers for adsorbate particles. The dependence of island size and density on misfit, temperature and particle flux is in good qualitative agreement with metal-organic vapor phase epitaxy experiments. Finally, we consider surface alloying as strain relaxation mechanism in multi-component systems. Equilibrium simulations show that the competition between binding and strain energy yields regular stripe patterns. Under non-equilibrium conditions we find pattern formation as well as the experimentally observed island ramification. The comparison with a lattice gas model shows that the stripe formation can solely be due to kinetic effects but the misfit is essential for the ramification and stabilization of the surface structures.

연구 동기 및 목표

• 재료 간 격자 불일치를 고려한 기계적 스트레인을 가진 이종외향성 성장에 대해 물리적으로 정확한 시뮬레이션 프레임워크를 개발하기 위해.
• 특히 큰 격자 불일치가 존재할 경우, 확산 메커니즘(특히 교환 확산)의 역할을 연구하기 위해.
• 불완전성 형성, Stranski-Krastanov 섬 형성, 다성분 체계에서의 표면 합금화를 포함한 주요 스트레인 완화 경로를 모의 및 분석하기 위해.
• 불완전성 핵형성 및 섬 형성에 대한 임계 두께와 격자 불일치, 온도, 기질 공급률 간의 정량적 관계를 수립하기 위해.
• 특히 투과성 ZnSe/GaAs 및 자가조직화된 InAs 양자점에 대한 실험 데이터와 시뮬레이션 모델을 검증하기 위해.

제안 방법

• 입자 간 거리와 격자 불일치에 따라 확산 장벽(이동 및 교환)을 계산하기 위해 분자 정적 방법을 활용하기 위해.
• 쌍체 포텐셜을 사용하여 입자 간 상호작용을 모의하고 표면 확산 사건의 활성화 에너지를 계산하기 위해.
• 계산된 활성화 에너기를 기반으로 개별 입자 사건을 추적하는 래티스 자유형 운동역학 몽테카를로(KMC) 알고리즘을 구현하기 위해.
• 성장 중에 결정 격자의 탄성 완화를 포함하기 위해, 응력 에너지를 최소화하도록 입자 위치를 동적으로 조정하기 위해.
• KMC 시뮬레이션에서 이외향성 성장의 큰 시간 스케일을 처리하기 위해 효율적인 데이터 구조 및 이벤트 선택 전략을 적용하기 위해.
• 표면 합금화를 연구하기 위해 다성분 체계를 시뮬레이션하며, 패턴 형성에서 화학 결합과 격자 스트레인 간의 경쟁을 고려하기 위해.

실험 결과

연구 질문

• RQ1 격자 불일치가 기계적 스트레인을 받는 이종외향성 체계에서 이동 확산 대비 교환 확산 메커니즘의 상대적 중요도에 어떤 영향을 미치는가?
• RQ2 불완전성 핵형성에 대한 임계 필름 두께와 격자 불일치 사이의 스케일링 관계는 무엇이며, 이는 거듭제곱 법칙을 따르는가?
• RQ3 불완전성, 섬, 표면 합금화와 같은 스트레인 완화 메커니즘은 온도, 기질 공급률, 격자 불일치에 따라 어떻게 달라지는가?
• RQ4 ZnSe/GaAs에서 관측된 투과성 성장과 점진적 스트레인 완화는 시뮬레이션을 통해 정량적으로 얼마나 잘 재현될 수 있는가?
• RQ5 주기적인 표면 패턴(예: 줄무늬, 가지치기 형태)의 형성은 운동역학적 효과만으로 설명될 수 있는가, 아니면 이들의 안정화에 격자 불일치가 필수적인가?

주요 결과

• 격자 불일치가 증가할수록 교환 확산 장벽이 크게 감소하여 고스트레인 조건에서 교환이 주요 확산 메커니즘이 된다.
• 불완전성 핵형성에 대한 임계 필름 두께는 격자 불일치에 대해 거듭제곱 법칙으로 스케일링되며, 실험 관측과 일치한다.
• ZnSe/GaAs에서의 투과성 성장과 점진적 스트레인 완화는 수직 격자 상수 변화의 현장 측정 결과와 우수한 정성적 일치를 보였다.
• 투과성 성장 영역은 격자 불일치에 대해 약간의 거듭제곱 법칙 스케일링을 보이며, 이는 모델의 예측 능력을 뒷받침한다.
• Stranski-Krastanov 섬 형성은 섬 표면에서의 느린 확산과 표면 완화로 인한 내향적 기질 공급률에 의해 주도되며, 섬의 크기와 밀도는 온도, 기질 공급률, 격자 불일치에 따라 달라진다.
• 표면 합금화로 인해 주기적인 줄무늬 패턴과 가지치기 형태의 구조가 형성되며, 이는 화학적 농도 기울기가 없더라도 격자 불일치가 가지치기의 안정화 및 촉진에 기여한다.