GaN의 육방정계·입방정계에서 스태킹 폴트의 에너지와 전자구조 DFT 전산모델링

초록

본 연구는 DFT를 이용해 육방정계(wz)와 입방정계(zb) GaN의 스태킹 폴트(I1, I2, I3, E 및 intrinsic, extrinsic, twin) 구조와 형성 에너지를 정량화하고, 각 결함이 밴드 구조에 미치는 영향을 분석한다. wz‑GaN에서는 I1이 가장 안정적인 결함이며, zb‑GaN에서는 intrinsic 결함이 우세하다. 전자밀도 분석 결과, 결함면 근처의 원자가 밴드와 전도밴드가 서로 반대쪽 Ga‑N 극성에 국한되며, 이는 type‑II 밴드 정렬을 초래한다. 이러한 결과는 GaN 기반 광·전자소자의 결함 제어와 설계에 중요한 통찰을 제공한다.

상세 분석

이 논문은 대규모 병렬 DFT 코드인 Conquest을 활용해 실험적으로 관측되는 다양한 스태킹 폴트의 원자 수준 구조와 열역학적 안정성을 정밀하게 계산하였다. GGA‑PBE 교환‑상관 함수와 Hamman의 정규화 퍼시듀얼 포텐셜을 사용했으며, 실공간 격자 간격을 0.1 Bohr(≈400 Ha)로 설정해 총 에너지 수렴을 1 µHa 이하로 확보하였다. wurtzite 구조는 ABAB…의 6‑층 주기성을 갖는 (0001)면을 기준으로 I1(AB/CB), I2(AB/CA), I3(AB/C/BA) 및 extrinsic(E) 결함을 모델링했으며, zincblende 구조는 ABCABC…의 (111)면을 기준으로 intrinsic(ABC/BC), extrinsic(ABC/B) 및 twin(ABC/BAC) 결함을 구축하였다. 각 결함은 주기적 경계조건 하에 두 개 이상의 결함을 포함하도록 설계해 상호작용을 최소화하고, 층 수(N)를 늘려 형성 에너지 수렴성을 검증하였다. 결과적으로 wz‑GaN의 I1 형성 에너지는 17 mJ m⁻²로 가장 낮아 실험에서 90 % 이상 관찰되는 사실과 일치한다. I2, I3, E는 차례로 28, 36, 55 mJ m⁻²로 점차 불안정해진다. 반면 zb‑GaN에서는 모든 결함이 음의 형성 에너지(≈‑33 mJ m⁻²)를 보여, wz 삽입이 zb 매트릭스에 자연스럽게 발생함을 의미한다. 전자구조 분석에서는 밴드‑해석 전하밀도(ρₙ(r))를 시각화해 결함면 양쪽에서 전자와 정공이 서로 다른 극성(Ga‑polar vs N‑polar) 쪽에 국한되는 현상을 확인했다. wz‑GaN에서는 가전자대 최상단이 Ga‑polar 쪽에, 전도대 최하단이 N‑polar 쪽에 집중되는 반면, zb‑GaN에서는 그 반대가 나타난다. 이는 결함면 근처 전기적 포텐셜 차이가 반대 부호를 갖기 때문이며, 두 물질 모두 type‑II 밴드 정렬을 보인다. 또한, 밴드 오프셋이 주변 원자 배치와 전하 재분포에 민감하게 변함을 보고했으며, 이는 결함이 국소적으로 밴드갭을 감소시켜 양자우물(QW) 효과를 유도할 수 있음을 시사한다. 이러한 원자‑레벨 인사이트는 GaN 기반 LED·레이저에서 QCSE 억제, 결함‑유도 양자구조 설계, 그리고 비극성 zb‑GaN의 효율 향상 전략에 직접적인 활용 가능성을 제공한다.