GaInNAs 희석 질화물 양자우물: 10-밴드 징크블렌드 k·p 튜토리얼

1. 서론

이 튜토리얼에서는 OghmaNano의 10-밴드 징크블렌드 k·p 모델을 활용하여 희석 질화물로 구성된 GaInNAs 양자우물을 시뮬레이션합니다. 기존 예제 파일을 열어 층 구조를 확인한 뒤 전자상태 계산을 수행하고 예측된 광학 스펙트럼을 분석하게 됩니다. 해당 예제에서는 3개 층만으로 구성된 독립형 양자우물이 아니라, 보다 큰 GaAs/AlGaAs 구조 내에 위치한 두께 7 nm의 GaInNAs 양자우물이 포함되어 있습니다.

10-밴드 모델이 왜 유용한가? 일반적인 8-밴드 k·p 이론은 스핀까지 포함하여 전도대와 중정공대, 경정공대, 스플릿-오프 가전대를 기술한다. III-V 화합물 반도체에 소량의 질소를 첨가하면 질소와 관련된 공명 준위가 생성되며 이 준위는 기저체의 전도대와 특히 강하게 상호작용한다. 10-밴드 모델은 이러한 질소 관련 준위의 두 가지 스핀 성분까지 추가로 기술한다. 이러한 상호작용은 기본 전이 에너지를 현저히 낮추고 전도대의 분산곡선 형태도 변화시킬 수 있다. 이러한 특성 덕분에 희석 질화물 반도체는 장파장 광학 소자 개발에 매우 유용하다.

2. 10-밴드 예제 열기

파일 → 새 시뮬레이션을 선택합니다. 새 시뮬레이션 창에서 GaAs 데모 항목(??)을 더블클릭하세요. 그 다음 창에는 여러 개의 양자우물 시뮬레이션 예제들이 나타납니다. AlGaAs 희석 질화물 10-밴드 k·p 항목(??)을 선택하면 사전 설정된 10-밴드 시뮬레이션이 로드되며, 여기에는 적절한 층 구조 정의와 재료 파라미터 파일들이 포함됩니다.

3. 소자 구조 및 재료 파라미터

메인 창에는 다층 구조가 표시됩니다(??). 이전 튜토리얼들에서 사용된 단순한 GaAs/AlGaAs 양자우물과는 달리, 이 예제에서는 희석 질화물 양자우물 주위에 클래딩층, 농도구배층 및 장벽층이 추가로 존재합니다. 구조 형상을 확인하려면 왼쪽 툴바에서 ‘레이어 편집기’를 클릭하십시오.

레이어 편집기(??)에는 p-GaAs부터 시작해 AlGaAs, 활성층을 거쳐 n-GaAs에 이르기까지 총 8개의 레이어 항목이 표시됩니다. 특히 중앙에 위치한 GaInNAs(QW) 레이어의 두께는 7.0 nm이며 레이어 유형은 양자우물로 설정되어 있습니다. 해당 예시 구조에서는 활성층 양측에 각각 0.16 µm 두께의 AlGaAs 레이어와 0.02 µm 두께의 GaAs 장벽/SCL 레이어도 포함되어 있습니다. 제공된 템플릿을 활용할 경우 스택 구조를 직접 재구성할 필요가 없습니다.

3.1 재료 정의 사항 확인

전기적 파라미터 메뉴를 클릭한 뒤 GaInNAs (QW) 탭을 선택하여 코드 편집기(??)를 엽니다. 스크린샷에서는 물질 파일에 유효 가전대 상태밀도 및 전하 캐리어 이동도와 같은 물리량들이 포함되어 있는 것을 확인할 수 있습니다. 예를 들어, 낮은 전기장 하에서 측정된 전자 이동도 값은 0.05 m² V⁻¹ s⁻¹로 표시되며 이는 500 cm² V⁻¹ s⁻¹에 해당합니다. 이 값들은 도식화된 양자우물의 분산곡선으로부터 추론된 수치가 아니라 물질 모델에 입력되는 실제 물성값입니다.

희석 질화물 시뮬레이션 수행 시 조성비, 밴드 가장자리 에너지 및 질소 결합 파라미터들이 서로 일관성을 유지하는 것이 매우 중요합니다. 첫 번째 실행에서는 템플릿에 지정된 값들을 그대로 사용하시기 바랍니다. 페이지 말미에서 유용한 조성비 실험에 대해 설명하도록 하겠습니다.

4. 질소를 도입하면 어떤 변화가 발생하는가?

모델을 실행하기에 앞서, 왜 해밀토니안이 8개가 아닌 10개의 대역을 갖는지 이해해두는 것이 중요합니다. 전자가 움직이는 평면 내 파동벡터가 증가할 때 기저 반도체의 전도대 가장자리 에너지가 상승한다고 가정해보겠습니다. 질소에 의해 생성된 공명 상태 역시 고유한 에너지값을 지닙니다. 만약 이 두 상태 간에 상호작용이 없다면 각각의 에너지값은 서로 가까워지거나 교차할 수 있습니다. 하지만 희석 질화물에서는 질소가 기저 반도체의 전도대 상태들과 결합하여 최종적인 에너지 준위가 변화하게 됩니다.

이 물리 현상을 설명하는 유용한 예시가 바로 2개의 상태로 구성된 밴드-안티크로싱(BAC) 행렬입니다:

\[ H_{\mathrm{BAC}}(k)= \begin{pmatrix} E_C(k) & V_N\\ V_N & E_N \end{pmatrix}, \qquad V_N=\beta_N\sqrt{x_N}. \]

여기서 \(E_C(k)\)은 기질의 전도대 분산곡선이며, \(E_N\)는 질소에 의해 생성된 에너지 준위의 값, \(x_N\)은 질소 함량, 그리고 \(\beta_N\)는 두 준위 간의 결합 강도를 조절합니다. 이 작은 행렬을 대각화하면 두 개의 분산곡선이 도출됩니다.

\[ E_{\pm}(k)=\frac{E_C(k)+E_N}{2} \pm\frac{1}{2}\sqrt{\big[E_C(k)-E_N\big]^2+4V_N^2}. \]

제곱근 항은 \(V_N\neq0\)일 때 결합된 에너지 분기들이 서로 교차하는 것을 방지한다. 특히 하위에 위치한 전도대 유사 분기의 에너지가 원래의 기질 전도대 에너지보다 낮아질 수 있으며, 이로 인해 기본 밴드갭이 감소하게 된다. 해당 상호작용은 \(E_C(k)\)와 \(E_N\) 사이의 거리에 의존하므로, 에너지 분산 특성 역시 강하게 비포물선적인 형태를 띠게 된다.

이 2×2 행렬 예시는 핵심 개념을 설명하기 위한 것이며, 본 튜토리얼에서 다루는 완전한 양자우물 해밀토니안은 아닙니다. 10-밴드 k·p 모델은 스핀까지 고려한 징크블렌드 구조의 전도대 및 가전대 기저를 그대로 유지하면서 질소와 관련된 구성 요소들도 추가합니다. 이후 성장 방향에 따른 양자 구속 효과도 포함됩니다. 완전한 이론적 설명은 10-밴드 징크블렌드 관련 장을 참조하십시오.

5. 시뮬레이션 실행하기

파일 리본으로 돌아가서 ‘시뮬레이션 실행’을 클릭합니다. 터미널 탭을 선택하여 계산 과정을 확인합니다(??). 솔버는 희소 해밀토니안을 구성한 뒤, 제한된 전도대 및 가전대의 고유상태들을 산출하며, 일련의 평면 내 파벡터 값에 대해 계산을 반복합니다(\(k_{\parallel}\)).

스크린샷에서 로그 기록을 통해 KP10 솔버가 식별되며, 계산된 에너지 준위가 출력되고 일정 횟수의 반복 계산 후 수렴이 보고됩니다. 서로 근접하게 배치된 에너지 쌍들에는 스핀 관련 축퇴 현상이 포함될 수 있으므로 이를 질소에 의한 밴드 반교차 현상으로 나타나는 별개의 분지로 자동 해석해서는 안 됩니다. 출력 파일을 확인하기 전에 시뮬레이션이 완료될 때까지 기다려주시기 바랍니다.

6. 계산된 밴드 구조 이해하기

생성된 그래프들은 장치 내 깊이에 따른 에너지 변화가 아니라, 평면 내 파수벡터에 따른 에너지 변화를 나타낸다. 양자우물의 퍼텐셜은 층 성장 방향으로 전하 캐리어들을 가두어 이산적인 에너지 상태들을 형성하며, 층들에 평행한 운동에 의해 각각의 가두어진 상태에 분산곡선이 생성된다 \(E_n(k_{\parallel})\). 단순한 포물선형 모델에서는 이러한 분산곡선이 이차함수 형태를 띠지만, 10개의 밴드가 상호작용하는 해밀토니안 모델에서는 포물선 형태에서 벗어날 수 있다.

6.1 전도대 상태와 감소된 밴드갭

개방전압 그래프(??)에는 서로 다른 에너지 값에서 시작하는 두 개의 계산된 전도대 분지가 나타나 있습니다. 가장 낮은 에너지를 갖는 전도대 상태는 선택된 시뮬레이션 에너지 기준점에서 약 0.89 eV의 값을 가집니다. 이 수치는 벌크 상태의 GaInNAs 밴드갭과 혼동해서는 안 되며, 이는 특정한 제한 구조에서의 고유에너지값입니다. 또한 두 분지 곡선의 곡률이 동일하지 않다는 점에도 주목해야 합니다. 이러한 특성은 8-밴드 계산 결과와 희석 질화물 모델을 비교할 때 검토해야 할 중요한 요소 중 하나입니다.

위에 제시된 단순화된 BAC 방정식들은 질소가 함유된 물질에서 전도대 분산곡선이 어떻게 변화할 수 있는지를 설명하는 데 도움이 됩니다. 다만 이 방정식들만으로는 두 색상으로 표시된 곡선 각각에 질소가 미치는 영향을 입증할 수는 없습니다. 이를 확인하려면 대응하는 고유벡터들을 검토하거나 \(x_N\) 값을 변화시켜 전하 캐리어 상태들이 어떻게 변하는지 관찰해야 합니다.

6.2 가전대 상태

개별 상태 분포 그래프(??)에는 비포물선 형태를 띠는 여러 분기들이 나타나 있습니다. 유한한 평면 파수 범위 내에서 중량정공과 경량정공의 특성이 혼합될 수 있으며, 양자 구속 에너지, 스핀-궤도 상호작용 및 기저 상태 간의 결합 역시 이러한 분포 곡선에 영향을 미칩니다. 곡선이 위쪽으로 꺾이거나 곡률이 변한다고 해서 반드시 수치적 오류가 발생한 것은 아니며, 반대로 매끄럽게 그려진 곡선만으로도 재료 파라미터가 올바르다는 증거가 되지 않습니다.

각 색상들은 단순히 고정된 HH 또는 LH 라벨을 전체 파동벡터 범위에 걸쳐 부여하기보다는, 실제로 도시된 고유상태 분산곡선으로 간주하는 것이 바람직합니다.

6.3 두 대역을 함께 살펴보기

결합된 분산곡선 그래프(??)를 열면 전자 상태와 정공대 상태가 동일한 에너지 척도 상에 표시됩니다. 기본적인 밴드 간 전이 현상은 각 상태의 절대적 에너지값이 아닌, 두 상태 간의 에너지 차이와 관련이 있습니다.

\[ \hbar\omega_{nm}(k_{\parallel}) =E_{c,n}(k_{\parallel})-E_{v,m}(k_{\parallel}). \]

이 식은 평면 내 파동벡터 방향에서 일어나는 직접적이고 대략 수직 방향의 광학적 전이를 나타냅니다. 즉, 광자의 운동량은 전자 결정의 운동량 규모에 비해 매우 작습니다. 다만 실제 관측되는 스펙트럼은 일반적으로 하나의 날카로운 선 형태가 아닙니다. 여러 허용된 전이들이 기여하며, 각각의 전이에 대한 강도와 전하 캐리어 점유도는 서로 다르게 나타납니다.

7. 전자 상태로부터 광학 스펙트럼까지

계산된 상태들이 확보되면 OghmaNano는 광학 스펙트럼을 구성할 수 있다. 여기서 중요한 물리적 질문은 전자와 정공의 에너지 차이가 적절한지 여부뿐만 아니라, 두 상태가 전자기장과 얼마나 강하게 결합하는지 그리고 해당 초기 상태와 최종 상태가 적절히 채워져 있는지 여부이다.

7.1 왜 일부 전이 현상이 다른 전이 현상보다 더 크게 기여하는가?

페르미의 황금법칙은 전자 전이율을 이해하는 데 유용한 출발점이 됩니다. 약한 광학적 섭동에 의해 초기 전자 상태 \(i\)에서 최종 상태 \(f\)로의 전이에 대한 이 법칙의 개략적 형태는 다음과 같습니다:

\[ W_{i\rightarrow f}=\frac{2\pi}{\hbar} \left|\langle f|H_{\mathrm{int}}|i\rangle\right|^2 \delta(E_f-E_i-\hbar\omega). \]

행렬 요소 \(\langle f|H_{\mathrm{int}}|i\rangle\)는 각 상태와 빛의 상호작용을 기술합니다. 이 요소는 상태의 궤도 특성, 포락함수 및 빛의 편광 등에 의해 결정됩니다. 델타 함수는 에너지 보존 법칙을 만족시키는 역할을 합니다. 수치적 스펙트럼에서는 이상적으로 급격한 전이 현상이 적절한 선폭 확장 함수를 통해 표현됩니다. 허용된 모든 하위 대역 및 평면 내 파벡터에 대해 합산을 수행한 결과로 단일 전이 에너지가 아닌 전체 스펙트럼이 도출됩니다.

전하 캐리어의 점유 상태에 따라 흡수와 방출이 구분된다. \(f_c\)이 전도대 상태의 점유 수를, \(f_v\)가 가전대 상태의 점유 수를 나타낼 때 각각의 특성 인수는 다음과 같다:

\[ \begin{aligned} \text{absorption:}\quad &(1-f_c)f_v,\\ \text{stimulated emission:}\quad &f_c(1-f_v),\\ \text{net stimulated contribution:}\quad &f_c-f_v. \end{aligned} \]

마지막 식은 상향 및 하향 광학 과정들 간의 차이를 나타냅니다. 순이득에 긍정적인 기여가 발생하려면 충분한 인구 역전 현상이 필요합니다. 이러한 요인들이 다음 그래프들의 해석을 설명해줍니다. 보다 상세한 내용은 광학 전이 이론을 참조하십시오.

7.2 자극방출과 흡수

두 개의 패널(?? 및 ??)은 동일한 양자우물의 밴드 구조에서 발생하는 상호 경쟁적인 광학적 과정들을 보여준다. 여러 곡선들은 예제에서 사용된 조건들의 순서를 나타내며, 전하 캐리어 분포에 따라 어떤 전이 과정이 가능한지가 결정되므로 각 곡선의 상대적 크기도 달라진다.

시뮬레이션 결과 스펙트럼은 근적외선 영역까지 나타나며, 대략 1.1~1.5 µm 범위에 걸쳐 분포합니다. 참고로 광자 에너지와 진공 파장은 다음 관계식으로 연결됩니다:

\[ E_{\mathrm{photon}}=\frac{hc}{\lambda} \simeq \frac{1240\ \mathrm{eV\,nm}}{\lambda\,[\mathrm{nm}]}. \]

따라서 1300 nm 파장은 약 0.954 eV에 해당합니다. 이 값은 도파관 효과, 기여하는 전이 현상들의 범위 및 선택된 에너지 기준점을 고려할 때, 결합된 밴드 분포도상의 에너지 차이와 광학 스펙트럼의 척도 간의 일관성을 확인하는 데 유용합니다.

7.3 순이득: 방출 현상이 흡수 현상을 압도하는 경우

증폭 매질에서 주요 관심 대상은 자극방출과 흡수 간의 균형입니다. 해당 그래프에서 사용된 부호 규칙에 따라 이 관계는 다음과 같이 개략적으로 표현될 수 있습니다:

\[ g_{\mathrm{net}}(\lambda) =g_{\mathrm{stim}}(\lambda)-\alpha(\lambda). \]

계산된 순이득 그래프(??)에는 음수값과 양수값을 갖는 곡선들이 모두 포함되어 있다. 흡수 현상이 우세할 경우 \(g_{\mathrm{net}}\) 값은 음수가 된다. 전자와 정공의 분포가 변화함에 따라 자극방출 현상이 충분히 강해져 스펙트럼의 일부 영역에서 양의 값을 갖는 순이득이 나타날 수 있다. 여러 개의 국소화된 상태 간 전이 현상들이 복합적으로 작용하기 때문에 그래프상의 전이 영역은 매끄러운 단일 피크 형태가 아닌 복잡한 구조를 보인다.

중요한 점은, 순이득이 양수 값을 나타낸다고 해서 반드시 레이저가 임계값 이상에서 작동할 것이라고 예측할 수는 없다는 점입니다. 공진기 손실, 광학적 밀폐성, 거울의 반사율 및 전하 캐리어 주입량 역시 반드시 고려되어야 합니다.

8. 시뮬레이션 결과 확인하기

계산이 완료되면 Output 탭(??)을 열어보세요. 스크린샷에는 주요 출력 파일들이 나타나는데, 여기에는 Ec.csv (전도대 분산곡선), Ev.csv (가전대 분산곡선), Ec_Ev.csv (결합된 밴드 구조), gain.csv, absorption.csv, net_gain.csv, spontaneous.csv이 포함됩니다. 그 외 결과 데이터들은 bands 및 gain 디렉터리 내에 정리되어 있습니다. CSV 파일을 더블클릭하면 그래프가 생성되며, 저장된 수치 데이터를 활용해 추가적인 분석도 수행할 수 있습니다.

이 튜토리얼에서는 자극방출, 흡수 및 순이득에 대해 설명하지만, 해당 예제에서는 자발방출에 관한 출력 결과도 생성합니다. 자발방출이란 외부 광자가 존재하지 않아도 발생하는 방사성 재결합 현상이며, 이에 대한 특성 점유 인자는 \(f_c(1-f_v)\)입니다. 자발방출의 절대 스펙트럼 선행계수는 자극방출의 경우와 서로 다릅니다.

9. 직접 시도해보기: 질소와의 상호작용을 연구해보십시오

BAC 시각화 결과에 따르면 질소 결합 강도를 증가시키면 결합된 에너지 분지 간의 분리 정도가 커지고 하위 전도대 에너지도 변화하는 것으로 예측됩니다. 그러나 전체 합금 계산에서는 조성 변화가 다른 물질 파라미터에도 영향을 미칠 수 있습니다. 따라서 시뮬레이션에서 활용 가능한 실제 물질 함수를 기반으로 결과 변화를 해석해야 하며, 단순한 2상태 공식이 전체 합금 모델이라고 가정해서는 안 됩니다.

두 번째 실험에서는 재료 파라미터를 고정한 상태에서 양자우물의 폭을 7 nm에서 6 nm로 변경합니다. 양자 구속 효과가 증가하면 일반적으로 가장 낮은 에너지 준위를 갖는 전자의 에너지가 상승하며 허용되는 간대 전이 에너지들도 변화합니다. 이러한 에너지 변화의 크기는 각 상태에 단일한 유효질량 값을 할당하는 방식이 아니라 10개 대역을 모두 고려한 계산에 의해 결정됩니다.