GaN/AlGaN 양자우물: 8-밴드 육방정계 k·p 모델 튜토리얼
1. 서론
이 튜토리얼에서는 OghmaNano의 8-밴드 육방정계 k·p 모델을 활용하여 두 개의 AlGaN 장벽층 사이에 위치한 GaN 양자우물의 전자상태 및 광학적 특성을 계산합니다. 해당 예제에서는 두께 10 nm인 GaN 양자우물과 두께 20 nm인 AlGaN 층 두 개가 사용되었습니다. 먼저 제공된 예제를 실행한 뒤 전도대 및 밸런스 밴드의 부대역 분산 특성과 이에 따른 광학 스펙트럼을 검토할 것입니다.
육방정계와 스펀지형 결정 구조 간의 차이는 이 맥락에서 매우 중요합니다. 두 구조 모두 반도체 결정 구조이지만, 육방정계는 육방형 c축을 갖습니다. 이러한 결정 대칭성의 차이는 밸런스 밴드의 분리 정도, 전하 캐리어 상태들 간의 결합 강도, 그리고 서로 다른 방향으로 편광된 빛의 작용 강도에 영향을 미칩니다. 따라서 스펀지형 구조를 기반으로 한 8-밴드 계산 결과를 단순히 갈륨나이트라이드 재료 파라미터에 적용할 수는 없습니다.
2. 육방정계 8-밴드 예제 열기
파일 리본 메뉴에서 ‘새 시뮬레이션’을 클릭합니다. 나타나는 창에서 ‘GaAs demos’ 항목(??)을 더블클릭하세요. 카테고리 명칭과는 달리 이 항목에는 질화물 및 실리콘 양자우물 예제들도 포함되어 있습니다. 그 다음 창에서 ‘GaN Wurtzite 8 band k.p’ 항목(??)을 더블클릭하면 미리 설정된 재료 층 구조가 열리며 해당 데모가 선택됩니다.
3. GaN/AlGaN 양자우물 검토하기
메인 창에는 3개 층으로 구성된 구조가 표시된다(??). 실제 양자우물 계산에서는 이 층들의 두께 방향으로 전하 캐리어가 가두어지게 되는데, 전자 엔벨로프 함수는 각 층을 따라 변화하고 분산 그래프는 층들에 평행한 방향의 운동을 나타낸다.
왼쪽 툴바에서 ‘레이어 편집기’를 클릭하십시오. ??에 나타난 바와 같이, 예제에서는 다음과 같은 기하구조를 사용합니다:
| 층 | 두께 | 유형 |
|---|---|---|
| AlGaN | 20 나노미터 | 활성층 / 장벽층 |
| GaN (양자우물) | 10 나노미터 | 양자우물 |
| AlGaN | 20 나노미터 | 활성층/장벽층 |
첫 번째 시뮬레이션 실행 시 제공된 재료 설정값을 그대로 사용하십시오. 스크린샷에서는 층의 두께만 나타나며 알루미늄 몰 분율, 응력 구성 또는 분극 파라미터는 표시되지 않습니다. 따라서 특정 실험용 웨이퍼를 모사하기 위해 해당 예제를 활용하기 전에 재료 설정값을 반드시 확인해야 합니다.
4. 왜 8-밴드 육방정계 해밀토니안을 사용해야 하는가?
단순한 포물선 모델에서는 각 전하 캐리어에 대해 유효질량을 할당하고, 그 에너지가 면내 파동벡터의 제곱에 비례하여 변화한다고 가정합니다. 이 모델은 밴드 극값 근처에서는 유용하지만, 파동벡터가 증가할 때 여러 밸런스 밴드 상태들이 서로 어떻게 혼합되는지는 설명할 수 없습니다. 육방정계 구조를 갖는 물질에서는 육방정계 결정 구조로 인해 c축 방향의 운동 및 빛의 편광 현상이 수직 평면 내에서의 현상과 구별되기도 합니다.
8-밴드 모델에서는 전도대와 유사한 궤도 하나와 두 개의 스핀 성분을 갖는 밸런스 밴드와 유사한 궤도 세 개가 유지된다. 개념적으로 고유상태는 다음과 같은 형태를 갖는다
\[ \Psi_{n,\mathbf{k}_{\parallel}}(\mathbf{r})= e^{i\mathbf{k}_{\parallel}\cdot\mathbf{r}_{\parallel}} \sum_{\alpha=1}^{8}F_{n\alpha}(z)u_{\alpha}(\mathbf{r}), \]
여기서 함수 \(F_{n\alpha}(z)\)은 양자우물 내에서의 전하 캐리어의 구속 상태를 기술하며, \(u_\alpha\)는 8개의 블로흐 기저 상태를 나타냅니다. 따라서 특정 구속 상태는 하나 이상의 밸런스 밴드 궤도로부터 기여를 받을 수 있습니다. 이러한 이유로 밴드 가장자리에서 벗어난 영역에서는 각 색상으로 표시된 곡선을 단일한 무거운정공 또는 가벼운정공 라벨로만 설명하는 것은 오해의 소지가 있습니다.
육방정계 구조에서 밸런스 밴드의 배열 순서는 결정장 효과와 스핀-궤도 상호작용의 영향을 받습니다. 간략한 표기법을 사용할 경우, 해당 계산 과정은 다음과 같이 수행됩니다.
\[ H_{8\times8}\!\left(k_x,k_y,-i\frac{\partial}{\partial z}\right) \mathbf{F}_n(z)=E_n(\mathbf{k}_{\parallel})\,\mathbf{F}_n(z). \]
면내 파동벡터는 명시적으로 유지되는 반면 성장 방향의 파동벡터는 미분 연산자로 대체된다는 점에 유의하시기 바랍니다. 여러 \(k_{\parallel}\) 값에 대해 이 고유값 문제를 반복 수행하면 아래와 같은 부대역 분산 그래프가 생성됩니다. 기저 함수, 결정장 항 및 재료 계수에 관한 보다 상세한 논의는 육방정계 k·p 이론 장에서 확인하실 수 있습니다.
5. 시뮬레이션 실행하기
파일 리본으로 돌아가서 ‘시뮬레이션 실행’ 버튼을 누르십시오. 터미널 탭을 열어 진행 상황을 확인할 수 있습니다. ??에서 터미널은 WZ KP8 솔버를 식별하고 희소 해밀토니안 행렬의 조립 과정을 보고하며 수렴 정보를 출력하고 샘플링된 파동벡터별로 계산을 진행합니다. 이 솔버는 각 파동벡터 지점에서 전자 및 정공 고유상태를 개별적으로 계산합니다.
여러 성분으로 구성된 엔벨로프 함수 문제이므로, 각 밴드에 대해 별도의 스칼라 포물선형 양자우물을 해석하는 것보다 훨씬 복잡한 과정입니다. 결과 파일을 확인하기 전에 고유값 해석기와 그 이후의 광학 계산이 모두 완료될 때까지 기다려주시기 바랍니다.
5.1 광학 전이 강도 확인하기
시뮬레이션이 진행된 후 터미널에는 광학 행렬 요소의 세기 값이 출력됩니다(??). 일반적인 c면 육방정계 구조에서는 TE가 c축에 수직인 광학 전기장을 의미하며, TM은 c축에 평행한 성분을 나타냅니다. 이러한 편광들은 특정 전자 전이에 대해 동일한 강도로 결합하지 않을 수도 있습니다.
광학적 강도를 이해하는 데 유용한 방법은 전이 행렬 요소의 제곱값을 활용하는 것입니다.
\[ M_{cv}^{(\mathbf{e})}(\mathbf{k}_{\parallel})= \left\langle\Psi_{c,\mathbf{k}_{\parallel}}\middle| \mathbf{e}\!\cdot\!\mathbf{p}\middle| \Psi_{v,\mathbf{k}_{\parallel}}\right\rangle, \qquad S_{cv}^{(\mathbf{e})}\propto\left|M_{cv}^{(\mathbf{e})}\right|^2. \]
여기서 \(\mathbf{e}\)은 빛의 편광 상태를 나타내며 \(\mathbf{p}\)는 운동량 연산자입니다. 특정 상태의 광학 행렬 요소 값이 작을 경우, 해당 상태가 전자 밴드 구조 내에 존재하더라도 특정 광학 편광에 거의 기여하지 않을 수 있습니다. 따라서 최종 결과에서 나타나는 TE/TM 값의 비율은 에너지 준위 위치 정보 외에도 추가적인 물리적 정보를 제공합니다.
6. 계산된 전자 밴드 구조의 이해
각 밴드 그래프에서 가로축은 양자우물 내 위치가 아니라 면내 파동벡터\(k_{\parallel}\)입니다. 세로축은 국소화된 부대역의 에너지를 나타냅니다. \(k_{\parallel}=0\) 지점에서 각 곡선은 개별적인 국소화 에너지 값 중 하나에서 시작됩니다. 오른쪽으로 이동함에 따라 전하 캐리어가 층들에 평행하게 이동할 때 동일한 고유상태가 어떻게 분산되는지를 확인할 수 있습니다.
6.1 전도대 고유상태
개방전압 출력 결과를 열어 ??에서 밴드 분산 현상을 살펴보십시오. GaN 밴드 가장자리 근처의 여러 에너지 지점에서 시작되는 여러 개의 국소화된 준위들이 관찰됩니다. 파동벡터가 작은 영역에서는 전도대 부대역에 대해 다음과 같은 유용한 근사식이 성립합니다
\[ E_{c,n}(k_{\parallel})\simeq E_{c,n}(0) +\frac{\hbar^2 k_{\parallel}^2}{2m_{e,\parallel}^{*}}. \]
\(E_{c,n}(0)\) 항에는 양자우물 내에서의 전하 캐리어 구속 에너지가 포함되어 있으며, \(m_{e,\parallel}^{*}\) 항은 평면 내 운동에 대한 전하 캐리어의 반응을 나타냅니다. 이 관계식은 대략적으로 상승하는 곡선 형태를 설명하는 데 기여하지만, 실제 8-밴드 모델 계산에서는 고정된 포물선 형태를 가정하는 대신 밴드 간 상호작용을 그대로 고려합니다.
6.2 밸런스 밴드 고유상태
밸런스 밴드 그래프(??)는 보다 복잡한 형태를 보인다. 상위 대역은 처음에는 \(k_{\parallel}=0\)로부터 상승하다가 이후 하강하는 경향을 보이며, 여러 하위 대역들은 서로 다른 곡률을 나타낸다. 이러한 형상들은 결정장 분리, 스핀-궤도 결합, 양자구속 효과 및 밴드 혼합이 상호작용하는 육방정계 밸런스 밴드 모델에서 실현 가능하다.
이 그래프를 통해 모든 색상으로 표시된 곡선들이 파동벡터 범위 전체에 걸쳐 각각 고유하게 A-, B- 또는 C-대역에 해당한다고 해석해서는 안 됩니다. 대역을 명확히 식별하려면 기저 고유벡터 성분들의 기여도를 검토해야 합니다. 마찬가지로, 이 모델을 정량적으로 활용할 경우 비정상적으로 중심에서 벗어난 극값들도 선택된 재료 및 응력 파라미터와 비교해볼 필요가 있습니다.
7. 광학적 전이 현상, 유도방출 및 흡수
전자 밴드 구조는 가능한 광학적 전이의 에너지를 결정하지만, 전이 강도는 전자 및 정공의 파동함수, 편광 상태, 그리고 캐리어 분포에도 의존한다. 동일한 평면상 파동벡터에서 수직 전이에 해당하는 광자 에너지는
\[ \hbar\omega=E_{c,n}(\mathbf{k}_{\parallel}) -E_{v,m}(\mathbf{k}_{\parallel}). \]
페르미의 황금법칙은 터미널에 출력된 광학 행렬 요소들이 왜 중요한지를 설명합니다. 개략적으로 표현하자면, 전이 발생률은 다음과 같습니다.
\[ W_{v\rightarrow c}^{(\mathbf{e})}= \frac{2\pi}{\hbar}\left|H_{cv}^{\prime(\mathbf{e})}\right|^2 \delta\!\left(E_c-E_v-\hbar\omega\right). \]
섭동 행렬 요소 \(H'_{cv}\)는 빛의 편광 상태와 관련된 양자 상태들에 의존한다. 디랙 델타 함수는 에너지 보존 법칙을 나타내며, 실제 스펙트럼을 구성할 때 허용된 전이 과정들의 기여도는 샘플링된 양자 상태들에 대해 합산된 후 폭이 넓어진다. 특정 파장에서 흡수 현상이 주로 일어나는지 아니면 자극 방출 현상이 주로 일어나는지는 전하 캐리어 점유도에 의해 결정된다. 편리한 계산용 순수 인자로는 \(f_c-f_v\)가 사용되는데, 여기서 \(f_c\)와 \(f_v\)는 각각 전도대 및 밸런스 밴드 상태에 존재하는 전자의 점유도를 의미한다.
GaN 양자우물은 전자 밴드갭이 넓기 때문에 계산된 스펙트럼은 근자외선 영역에 나타나며, 제공된 예시에서는 주로 330~370 nm 범위에 분포한다. 파장과 광자 에너지의 관계는 \(E_{\gamma}=hc/\lambda\)에 의해 결정된다.
7.1 자극방사와 순이득
자극방사 플롯(??)에는 다양한 전하 캐리어 조건에 따른 여러 스펙트럼이 나타나 있습니다. 순이득 계산(??)은 방사 과정과 흡수 과정 간의 상호작용을 고려합니다. 양의 값을 갖는 물질 이득은 해당 전자 모델 내에서의 증폭 현상을 의미하지만, 이만으로는 레이저 공진기가 임계값에 도달했음을 입증할 수 없습니다. 구속 인자와 공진기 손실 역시 반드시 고려되어야 합니다.
7.2 흡수
흡수 결과는 ??에 나타나 있습니다. 흡수 스펙트럼의 경계 및 특성들은 점유된 밸런스 밴드 상태와 가용한 전도대 상태 간의 에너지 차이에 의해 결정되며, 이때 광학 행렬 요소가 가중치로 작용합니다. 육방정계 양자우물 내에서는 각 상태의 궤도적 특성에 따라 TE 및 TM 편광에 대한 결합 강도가 달라지므로 이에 따른 스펙트럼 역시 서로 차이를 보일 수 있습니다.
해당 스펙트럼에 나타나는 특징들은 설정된 모델 및 스펙트럼 폭화 현상의 결과입니다. 이를 특정 엑시톤 피크와 자동적으로 연관시켜서는 안 됩니다. 단일 입자 기반의 k·p 간대 계산만으로는 전자-정공 간 쿨롱 상호작용에 의해 결합된 엑시톤을 완벽하게 기술할 수 없기 때문입니다.
8. 자발적 방사
자극방출과는 달리, 자발방출은 전이를 유발하는 입사 광자가 없이도 발생한다. 이 과정은 전도대에 존재하는 전자와 밸런스 밴드에 빈 상태인 정공의 존재 여부에 의존한다. 따라서 전자 점유 인자는 다음과 같이 비례한다.
\[ R_{\mathrm{sp}}(\omega)\ \propto\ \sum_{c,v,\mathbf{k}_{\parallel}} \left|M_{cv}(\mathbf{k}_{\parallel})\right|^2 f_c(1-f_v)\,L\!\left(E_c-E_v-\hbar\omega\right), \]
여기서 \(L\)은 선택된 스펙트럼 확산 정도를 나타내며, 비례 관계식에서는 광자 관련 계수들이 생략되어 있다. 계산 결과는 ??에 도시되어 있다. 로그 스케일을 사용함으로써 여러 수십 배에 달하는 범위에 걸친 스펙트럼들을 비교할 수 있다. 스펙트럼의 단파장 영역에서 나타나는 광학적 특성들과 여러 조건 하에서 전체 강도의 변화 양상을 주목해보기 바란다.
9. 직접 시도해보기
💡 유용한 후속 계산 두 가지:
- GaN 양자우물의 폭을 변경해보세요. 중앙에 위치한 10 nm 두께의 층을 7 nm로 줄이고 장벽층의 두께는 그대로 유지해보시기 바랍니다. 그 후 고유상태를 재계산하여 가장 낮은 에너지 준위를 갖는 전도대 상태들과 광학적 전이 에너지들을 비교해보세요. 강한 양자구속 효과는 일반적으로 구속된 에너지 준위들 간의 간격을 넓혀주지만, 정공대의 반응 양상은 여러 밴드 간 상호작용에 의해 복잡하게 나타납니다.
- TE와 TM 전이 강도를 비교해보십시오. 터미널에서 제공되는 광학 강도 값을 활용하여 특정 편광 방향을 선호하는 전이들을 식별하십시오. 모든 전이가 스펙트럼에 동일하게 기여한다고 가정하지 말고, 이와 같은 편광 선호도가 정공대 상태의 특성과 어떻게 연관되는지 살펴보십시오.
✅ 무엇을 주목해야 하는가?
좁은 폭의 양자우물은 전자들의 국한 에너지를 상승시키고 정공 에너지 준위에도 변화를 일으켜 가능한 간밴드 전이 에너지값을 변동시킨다. TE/TM 비율은 육방정계 구조에서의 파동함수가 지닌 서로 다른 궤도 성분들을 반영하므로 상태 쌍마다 현저하게 달라질 수 있다. 정량적인 재료 분석을 위해서는 합금 조성, 응력 및 편극장이 미치는 영향도 함께 검토해야 한다.