에피택셜 성장 전이와 새로운 저이동성 영역: 1+1 차원 KMC 연구

초록

본 논문은 입자 플럭스와 온도를 조절한 1+1 차원 에피택셜 성장 모델을 제시하고, 광범위한 Kinetic Monte Carlo 시뮬레이션을 통해 BD, BDLR, KR, LBL 등 네 가지 전형적인 성장 레짐을 재현한다. 특히 저이동성 한계에서 기존 지표로는 구분되지 않던 새로운 레짐을 발견하고, β 지수, 부착 결합밀도 및 그 변동을 이용해 전이점을 정량화한다. 또한 열적 거칠기 전이와 동적 거칠기 전이 사이의 온도 관계를 제시한다.

상세 분석

이 연구는 에피택셜 박막 성장의 비평형 특성을 최소한의 파라미터(플럭스 F, 온도 비율 φ = J_pp/k_BT, 그리고 확산 장벽 비율 R_E = E_m/J_pp)로 기술한다. 입자 이동은 초기·최종 배위수 차이에 비례하는 활성 장벽 E_i = E_m − (n_f − n_0)J_pp 로 정의하고, Arrhenius 식 r_i = ν exp(−E_i/k_BT) 로 전이율을 계산한다. 이때 모든 과정에 동일한 시도주파수 ν를 가정함으로써 모델의 복잡성을 크게 낮추었다.

시뮬레이션은 10^3 사이트, 10^3 샘플 평균을 사용했으며, 시간 증분 Δt = −ln ξ/∑_i r_i 로 구현하였다. 표면 거칠기 W(t)와 성장 지수 β (W ∝ t^β) 를 주요 관측량으로 삼아 플럭스와 φ 의 함수로 레짐을 구분한다.

1. BD (Ballistic Deposition): 높은 플럭스·낮은 φ 에서 입자들이 거의 재배열 없이 바로 착착해 그림자 효과가 강해 거칠기가 급격히 증가한다. β≈0.5 수준이며, 부착 결합밀도 ρ_b 가 플럭스에 무관하게 일정하게 유지된다.

2. BDLR (Ballistic Deposition with Local Relaxation): 플럭스가 감소하거나 φ 가 증가하면 입자들이 인접한 고배위 자리로 이동하면서 부분적인 재구성이 일어난다. β는 여전히 높은 편이지만, ρ_b 가 감소하기 시작한다. 특히 ρ_b 변동(ϕ_b) 의 피크가 나타나는 지점이 BD→BDLR 전이점으로 해석된다. 피크가 시스템 크기에 따라 지속되는 점은 전이가 열역학적 한계에서 불연속적(1차 전이)임을 시사한다.

3. KR (Kinetic Rough): 플럭스가 더욱 낮아져 입자 확산 시간이 충분히 확보되면, 기존 층 위에 새로운 섬이 nucleation 되어 마운드가 형성된다. 이때 β가 다시 상승하고, 진공(빈자리) 밀도 ρ_v 가 감소한다. KR 레짐은 전형적인 비평형 거칠기 전이이며, β 값이 φ 에 따라 변동한다.

4. LBL (Layer‑by‑Layer): 플럭스가 매우 낮고 φ 가 충분히 커서 입자들이 완전히 평탄화된 후 다음 층이 시작된다. β→0에 수렴하고, 거칠기 진동이 뚜렷하게 나타난다.

5. Thermal Roughening Transition: φ 가 임계값 φ_c≈2 이하로 감소하면, 열적 평형 상태에서 거칠기가 로그적으로 증가한다. 이는 초기 평탄 표면을 일정 온도에서 장시간 유지했을 때 관찰되며, β가 0으로 수렴하지 않는다. φ_c 를 기준으로 열적 거칠기 전이와 동적 전이를 구분한다.

또한 저이동성 한계(R_E ≫ 1, φ → 0)에서 기존의 거칠기 지수나 부착 결합밀도로는 레짐을 구분할 수 없으며, 시뮬레이션 결과는 매우 높은 빈자리 비율과 거의 고정된 구조를 보인다. 이는 새로운 “저이동성 고정” 레짐으로 명명되었다.

마지막으로, 논문은 φ_c 와 β = 0 전이점 사이의 온도 관계를 제시한다. φ = J_pp/k_BT 로 정의되므로, φ_c 가 낮을수록 높은 온도에서 열적 거칠기가 발생하고, 반대로 높은 φ 에서는 동적 전이가 우세한다. 이는 실험적 플럭스·온도 제어를 통해 원하는 성장 모드를 선택할 수 있는 가이드라인을 제공한다.

전반적으로, 최소 파라미터 기반 KMC 모델이 복잡한 다층 성장 현상을 재현하고, 거칠기 지수, 부착 결합밀도, 그 변동 등 다양한 통계량을 이용해 전이점을 정량화함으로써 에피택셜 성장 제어에 실용적인 인사이트를 제공한다.