InGaNBi 합금은 차세대 포토닉 소자를 위한 유망한 소재로, 가시광에서 중적외선까지 밴드갭 튜닝이 가능한 다용도 플랫폼을 제공합니다. 제1원리 밀도범함수 이론을 사용해 인듐(In)과 비스무트(Bi) 농도가 GaN 기판 상 성장 시 격자 상수, 밴드갭 에너지, 스핀-오빗 분할, 변형에 미치는 영향을 체계적으로 조사했습니다. 격자 상수는 In과 Bi 증가에 따라 거의 선형적으로 증가하며, 밴드갭은 Bi 9.375% 이하, In 50% 이하에서 3.273 eV에서 0.651 eV(0.38-1.9 µm)까지 튜닝 가능합니다. 계산된 스핀-오빗 분할은 Bi 3.125%, 6.25%, 9.375%에서 각각 0.220 eV, 0.360 eV, 0.600 eV에 도달합니다. In/Bi 비율을 조절해 GaN 상에서 허용 가능한 변형을 유지하도록 합금을 설계할 수 있어 고효율 LED와 레이저 다이오드 개발의 길을 엽니다.
서론
버츠아이트 InₓGa₁₋ₓN 합금과 InGaN/GaN 양자우물(QW)은 튜닝 가능한 밴드갭과 GaN 기판과의 호환성으로 고효율 발광다이오드(LED)와 레이저 다이오드(LD)용으로 큰 관심을 받아왔습니다[1-10]. 그러나 [0001] 방향 InₓGa₁₋ₓN/GaN QW는 이축 압축 변형으로 인한 강한 내장 전기장으로 전자-정공쌍의 방출 에너지와 진동자 세기가 감소합니다[11]. 또한 기하학적 결함(적층 결함, 스레딩 전위)의 높은 밀도는 비복사 재결합 중심이 됩니다[12], 전자 누설과 오제 재결합은 InGaN LED의 악명 높은 효율 드룹에 기여합니다[13].
적외선 영역에서 비스무트(Bi)와 GaAs 합금화는 밴드갭 에너지(E_g) 감소와 스핀-오빗(SO) 분할 증대 전략으로 등장해 오제 재결합을 완화합니다[14]. Bi는 5족 최대 원소로 III-V 비스마이드 합금의 밴드 구조에 심대한 영향을 미치며, AlNBi[15], GaNBi[16,17], GaSbBi[18,19], InPBi[20,21], InSbBi[22-24] 등에서 입증되었습니다. 이들 삼원계에서 Bi 혼입은 Bi 불순물 상태와 중공/경공 밴드 간 상호작용을 통해 가격대(VB)를 교란시키고, Bi의 큰 원자 크기로 변형 관련 밴드갭 수정을 유도합니다[21]. 최근에는 GaAsNBi[25-27], InGaAsBi[28,29], GaAsPBi[30] 같은 사원소 비스마이드 시스템이 밴드갭 튜닝 범위 확장으로 주목받고 있습니다. Bi를 다른 III 또는 V 원소와 결합해 밴드갭, SO 분할, 변형을 동시에 엔지니어링함으로써 광전자 소자의 기능적 설계 공간을 확장할 수 있습니다[25].
InGaNBi가 가시광과 중적외선 응용을 연결할 잠재력을 고려해, 본 연구는 제1원리 계산으로 In과 Bi 농도가 합금의 구조적, 전자적, 변형 특성을 어떻게 형성하는지 규명합니다. 큰 격자 부정합과 제한된 Bi 용해도로 결정 품질이 저하되는 것을 피하기 위해 In ≤ 50%, Bi ≤ 9.375%로 제한했습니다[32,33]. 논문 구성: 2절 계산 방법, 3절 격자 상수, 밴드갭, SO 분할, 변형 결과 및 논의, 4절 주요 시사점 요약.
연구 방법
모든 계산은 Vienna Ab-Initio Simulation Package(VASP)[33,34]로 구현된 밀도범함수 이론(DFT) 프레임워크 내에서 수행했습니다. 구조 완화에는 투영 증강파(PAW) 방법[35,36]과 Perdew-Burke-Ernzerhof(PBE) 일반화 기울기 근사(GGA)를 교환-상관함수로 사용했습니다. PBE의 체계적인 밴드갭 과소평가를 보정하기 위해 수정된 Becke-Johnson(MBJLDA) 교환 퍼텐셜과 국부 밀도 근사(LDA) 상관을 결합해 사용했습니다[38]. Bi가 강한 스핀-오빗 결합(SOC)을 보이므로 모든 전자 구조 계산에 SOC를 명시적으로 포함했습니다.
원자가 전자 구성은 In: 4d¹⁰5s²5p¹, Ga: 3d¹⁰4s²4p¹, N: 2s²2p³, Bi: 5d¹⁰6s²6p³로 설정했습니다. 원자 힘은 <0.02 eV/Å, 총 에너지는 10⁻⁴ eV까지 수렴시켰습니다. 평면파 컷오프 450 eV, 4×4×4 Monkhorst-Pack k-포인트 메시를 사용해 정확한 총 에너지와 전자적 특성을 보장했습니다.
결과 및 논의
구조적 특성
36가지 조성의 InᵧGa₁₋ᵧN₁₋ₓBiₓ를 모델링하기 위해 4×2×2 슈퍼셀(64 원자)을 구성했습니다(0≤x≤0.09375, 0≤y≤0.5). 각 경우 In과 Bi 원자를 가능한 한 균일하게 분포시켜 클러스터링 효과를 최소화했습니다. 삼원계 InᵧGa₁₋ᵧN과 사원소 InᵧGa₁₋ᵧN₁₋ₓBiₓ 합금의 계산된 격자 상수를 그림 1에 도시했습니다. 순수 GaN의 경우 a=3.211 Å, c=5.235 Å를 얻었으며, 이는 기존 이론값(a=3.155-3.22 Å, c=5.144-5.24 Å)과 실험값(a=3.19 Å, c=5.19 Å)과 잘 일치합니다[39-42]. InᵧGa₁₋ᵧN의 격자 상수는 y에 따라 거의 선형적으로 증가해 y=0.25에서 a=3.304 Å, c=5.365 Å, y=0.5에서 a=3.397 Å, c=5.509 Å에 도달하며 기존 연구와 일치합니다[39,40,43,44]. 사원소 합금에서는 In과 Bi 함량 증가에 따라 a와 c 모두 거의 선형적으로 상승하며, 이는 In(1.42 Å)과 Bi(1.48 Å)의 이온 반경이 Ga(1.22 Å)와 N(0.71 Å)보다 크기 때문입니다[45].
a 삼원소 합금 InᵧGa₁₋ᵧN(0≤y≤0.5)과 b 사원소 합금 InᵧGa₁₋ᵧN₁₋ₓBiₓ(0≤x≤0.09375, 0≤y≤0.5)의 격자 상수. 비교를 위해 그림 1a에 다른 계산 및 실험 데이터[39-44] 추가. 실선은 a, 점선은 c
대표 조성 In₀.₂₅Ga₀.₇₅N₀.₉₃₇₅Bi₀.₀₆₂₅의 결합 길이 분석(그림 2) 결과 Ga-N 결합이 2.009 Å로 단축되는 반면 In-N 결합은 2.195 Å로 연장되고, Ga-Bi와 In-Bi 결합은 각각 2.592 Å, 2.704 Å에 도달했습니다. 이러한 경향은 In과 Bi 원자가 더 큰 공유 반경으로 국소 격자를 확장함을 확인하며, 이러한 격자 왜곡과 전기음성도 차이는 전자적, 광학적 거동에 영향을 미칠 것으로 예상됩니다.
In₀.₂₅GaNBi₀.₀₆₂₅의 결합 길이 히스토그램. 패널의 값은 네 가지 결합 유형의 평균 길이
전자적 특성
밴드갭 에너지는 III-V 합금에 신뢰성이 입증된 MBJLDA+SOC 접근법으로 계산했습니다. 그림 3은 InᵧGa₁₋ᵧN에서 In 함량에 따른 E_g 의존성을 실험 및 다른 이론 벤치마크(HSE06, mBJ, LMTO-CPA-MBJ)와 함께 보여줍니다. 계산된 GaN 밴드갭 3.273 eV는 실험값(3.40-3.50 eV) 및 이전 계산(3.261-3.33 eV)과 잘 일치합니다[39-49]. y가 0에서 0.5로 증가함에 따라 밴드갭이 3.273 eV에서 1.546 eV로 부드럽게 감소하며 이론과 실험과 일치합니다[50,51].
InᵧGa₁₋ᵧN에서 In 조성에 따른 예측 밴드갭 에너지(E_g, 빨간 실선) 및 데이터 피팅(검정 점선). 다른 이론[39,40,46] 및 실험[47-51] 결과도 함께 플롯
사원소 InᵧGa₁₋ᵧN₁₋ₓBiₓ 합금으로 확장하면, 그림 4의 등고선 지도에서 In과 Bi 함량 모두에 따라 밴드갭이 비선형적으로 크게 감소함을 볼 수 있습니다. Bi 9.375% 이하, In 50% 이하에서 밴드갭은 3.273~0.651 eV 범위로, 0.38~1.9 µm 파장에 해당해 가시광 전체와 중적외선 영역을 커버합니다. 중요하게도 Bi 혼입은 E_g를 낮출 뿐만 아니라 스핀-오빗 분할 Δ_SO를 향상시켜 Bi 3.125%에서 0.220 eV, 6.25%에서 0.360 eV, 9.375%에서 0.600 eV에 도달합니다. 이 값들은 In 비율에 크게 의존하지 않아 Bi가 Δ_SO 향상의 주된 동인임을 시사합니다. 더 높은 조성 합금에서 Δ_SO > E_g를 달성하면 오제 재결합을 억제할 수 있어 고효율 LED와 LD의 핵심 병목 현상을 해결할 수 있습니다.
InᵧGa₁₋ᵧN₁₋ₓBiₓ 합금의 밴드갭 등고선 플롯, Bi(x)와 In(y) 조성의 함수
GaN, In₀.₂₅GaN, In₀.₂₅GaN₀.₉₆₈₇₅Bi₀.₀₃₁₂₅의 투영 밴드 구조와 총 상태 밀도(TDOS)(그림 5)는 In과 Bi의 뚜렷한 역할을 보여줍니다. In 혼입은 전도대 최소(CBM)와 가격대 최대(VBM)를 모두 낮은 에너지로 이동시켜 밴드갭을 효과적으로 좁힙니다. 반면 Bi는 페르미 준위 바로 아래에 결함 밴드를 도입해 가격대 상단과 강하게 혼성화하며 VBM을 추가로 높입니다. TDOS는 -1.0~-0.5 eV 사이에 뚜렷한 Bi 유래 피크를 보여 캐리어 동역학을 조절할 수 있는 국소 상태 형성을 확인합니다.
a GaN, b, c In₀.₂₅Ga₀.₇₅N, d, e In₀.₂₅Ga₀.₇₅N₀.₉₆₈₇₅Bi₀.₀₃₁₂₅의 투영 밴드 구조와 해당 TDOS. 검은 점선은 페르미 준위(0 설정). In(파란색)과 Bi(빨간색) 기여 상태를 색으로 강조
GaN 상 InGaNBi의 변형
실제 소자 통합을 위해 GaN 위에 성장한 InᵧGa₁₋ᵧN₁₋ₓBiₓ 층의 변형 상태가 관리 가능한 범위 내에 있어야 합니다. 그림 6은 In과 Bi 비율에 따른 이축 변형을 매핑합니다. 최대 고려 조성(y=0.5, x=0.09375)에서 합금은 8.5% 압축 변형을 경험하는데, 이는 상당하지만 적절한 버퍼 층으로 잠재적으로 수용 가능합니다. 더 중요하게는 In 비율 ≤0.0625, Bi 비율 ≤0.028에서 변형이 1% 미만으로 유지되어, 적당한 합금 조성으로 준 격자 정합 시스템을 구현해 고품질 에피택시 성장이 가능함을 보여줍니다.
다양한 In(0-0.5) 함량에서 Bi 비율에 따른 GaN 기판 상 InGaNBi 합금의 변형. 양수 값은 InGaNBi가 압축 변형 상태임을 의미
결론
제1원리 DFT 계산은 InGaNBi 합금이 가시광에서 중적외선 스펙트럼까지 광전자 응용을 위한 고도로 튜닝 가능한 플랫폼을 제공함을 보여줍니다. 격자 상수는 In과 Bi 함량에 따라 거의 선형적으로 증가하며 In과 Bi의 더 큰 원자 크기를 반영합니다. 밴드갭은 Bi 9.375%, In 50%까지 3.273 eV에서 0.651 eV로 엔지니어링 가능하며, 스핀-오빗 분할은 최대 Bi 비율에서 0.600 eV에 도달해 Δ_SO > E_g 시 오제 재결합을 억제할 수 있습니다. 밴드 구조 분석에서 In은 주로 CBM과 VBM을 이동시키는 반면, Bi는 가격대 상단을 높이는 국소 결함 상태를 도입합니다. 변형 분석 결과 적당한 In(≤6.25%)과 Bi(≤2.8%) 농도에서 합금이 GaN과 거의 격자 정합되어 고품질 에피택시 성장이 가능함을 보였습니다. 이러한 발견은 맞춤형 방출 파장과 감소된 비복사 손실을 가진 고효율 InGaNBi LED 및 레이저 다이오드 설계의 토대를 마련합니다.