GaAs/AlGaAs 양자 우물: 8-밴드 아연석형 결정 구조 k·p 모델 튜토리얼

1. 서론

이 튜토리얼에서는 OghmaNano의 8-밴드 아연블렌드 구조 k·p 모델을 활용하여 GaAs/AlGaAs 양자 우물의 전자적 및 광학적 특성을 계산합니다. 먼저 미리 준비된 예제를 살펴보고 층 구조를 검토한 뒤 시뮬레이션을 실행하여 얻어진 전도대 및 정공대 서브밴드를 분석합니다. 마지막으로 해당 고유상태들로부터 계산된 광학 스펙트럼을 검토하게 됩니다.

이 포물선형 양자 우물 튜토리얼에서는 비교적 단순한 유효질량 밴드 모델을 사용하여 전자와 정공을 다룬다. 여기서 해밀토니안은 전도대 기저와 무거운 정공(HH), 가벼운 정공(LH), 스핀-궤도 결합에 의한 분리된 정공(SO) 밴드 기저로 구성되며 각각에 대해 두 개의 스핀 성분이 존재하므로 총 8개의 기저 상태가 형성된다. 이러한 상태들 간의 상호작용을 통해 본 모델은 독립적인 포물선형 밴드 모델로는 설명할 수 없는 현상들, 특히 밸런스 밴드 혼합 및 면내 파동벡터에 따른 밴드 곡률 변화까지도 기술할 수 있다.

2. 8-밴드 예제 열기

OghmaNano에서 파일 리본 메뉴의 ‘새 시뮬레이션’을 클릭합니다. 그 다음 GaAs 데모 폴더(??)를 더블클릭합니다. 양자 우물 예제들 중에서 ‘GaAs/AlGaAs QW 8 band k.p’(??)를 더블클릭하면 미리 설정된 시뮬레이션 창이 열립니다. 8×8 크기의 해밀토니안을 직접 입력할 필요는 없습니다.

3. 양자 우물 구조 및 재료 파라미터

메인 창에는 3개 층으로 구성된 AlGaAs/GaAs/AlGaAs 구조가 표시된다(??). 왼쪽에 있는 버튼을 이용해 레이어 편집기를 열 수 있다. 제공된 예제에서는 10 nm 두께의 GaAs 양자 우물을 40 nm 두께의 AlGaAs 장벽 두 개가 둘러싸고 있다.

GaAs의 밴드 가장자리들은 주변 AlGaAs의 밴드 갭 내에 위치하여 층 성장 방향으로 전하 캐리어의 구속을 가능하게 한다. 레이어 편집기에서 중앙에 위치한 GaAs 층은 양자 우물로서 특별히 표시되어 있다(??). 본 예제에서는 광학 재료 설정값 및 층 유형들이 이미 지정된 상태이다.

8-띠 모델 해석기에는 스칼라형 포물선띠 모델에 비해 더 많은 재료 관련 정보가 필요합니다. 즉, 띠 가장자리 위치와 유효질량 외에도 띠 간 결합을 기술하는 매개변수, 스핀-궤도 분리 에너지 및 필요 시 응력 값도 해밀토니안 계산에 사용됩니다. Lua 재료 파일 내의 GaAs 관련 함수들은 전기적 매개변수 설정 창 및 재료용 코드 편집기(??)를 통해 확인할 수 있습니다. 스크린샷에는 정공대 변형 퍼텐셜 값인 b 및 d가 나타나 있으며, 초기 시뮬레이션 수행 시 이 값들을 수정할 필요는 없습니다.

8-밴드 계산에서 어떤 변화가 발생하는가?

포물선 모델은 각 밴드에 별도의 곡률 값을 할당한다. 반면 k·p 모델은 대각선상이 아닌 요소들이 기저 상태들을 결합시키도록 행렬을 구성한다. 개략적으로, 구속 상태 문제는 다음과 같다.

\[ H_{8\times8}(k_x,k_y,-i\partial_z)\,\Psi_n(z) =E_n(k_x,k_y)\,\Psi_n(z). \]

여기서 \(z\)은 구속 방향을 나타내며 \(k_x\)와 \(k_y\)는 층에 평행한 운동을 기술합니다. 포락파 함수 \(\Psi_n\)는 8개의 성분으로 구성되는데 이는 서로 다른 기저 상태들의 기여도를 나타냅니다. 본 튜토리얼에서는 전체 행렬 대수보다는 계산된 에너지값에 초점을 맞출 것입니다. 보다 상세한 설명은 8-밴드 아연석 이론 장에서 확인할 수 있습니다.

4. 시뮬레이션 실행

파일 리본으로 돌아가서 파란색 ‘시뮬레이션 실행’ 버튼을 클릭하십시오. 그러면 8-밴드 해밀토니안이 구성된 뒤 여러 평면 내 파동벡터에 대해 해석이 수행됩니다. 터미널 탭(??)에서 계산 과정을 확인할 수 있습니다.

예시 터미널 화면에서는 희소 행렬이 생성되는 과정과 함께 고유상태 에너지값 및 block converged after 12 iterations과 같은 메시지들이 표시됩니다. 전도대 및 정공대 상태들은 연속적인 파동벡터 값에 대해 각각 계산됩니다. 계산이 완료되어 버튼이 다시 ‘시뮬레이션 실행’ 상태로 돌아올 때까지 기다린 후 결과값들을 확인하시기 바랍니다.

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

해당 솔버는 각 평면 내 파동벡터에 대해 구속된 고유상태들의 집합을 계산한다 \(k_{\parallel}\). 따라서 그림에서는 층들에 걸친 공간적 밴드 경계 형태가 아닌 서브밴드 분산 현상이 나타난다. Output 탭을 열어 bands 폴더를 확인하거나 ??에 표시된 요약 밴드 파일들을 열어보기 바란다.

5.1 전도대 서브밴드

개방전압 그래프(??)에는 위쪽으로 휘어진 두 개의 곡선이 나타나 있습니다. \(k_{\parallel}=0\) 지점에서 두 곡선 간의 수직 간격은 서로 다른 구속 고유상태에 기인합니다. \(k_{\parallel}\) 값이 커짐에 따라 에너지가 증가하는 현상은 양자 우물 내에서 캐리어가 구속되지 않은 채로 평면 방향으로 이동함을 의미합니다.

전도대 서브밴드의 최소점 근처에서는 분산 관계를 다음과 같이 근사할 수 있다

\[ E_n(k_{\parallel})\simeq E_n(0) +\frac{\hbar^2 k_{\parallel}^{2}}{2m_{n,\parallel}^{*}}. \]

다중밴드 모델에서는 이를 국소적 근사로 간주합니다. 즉, 도표에 나타난 전체 파동벡터 범위에 걸쳐 곡률 및 그에 따른 유효질량이 일정하게 유지될 필요는 없습니다. 전도대와 정공대 간의 상호작용 역시 이러한 특성에 영향을 미칩니다.

5.2 정공 서브밴드: 밸런스 밴드 혼합과 곡률

밸런스 밴드 그래프가 더 흥미롭습니다(??). 여기에 표시된 두 서브밴드는 각기 다른 에너지에서 시작하여 아래로 휘어지며 \(k_{\parallel}\) 값에서 서로 만나는 듯 보입니다. 아연블렌드 구조의 밸런스 밴드에서는 무거운 정공 상태와 가벼운 정공 상태가 서로 완전히 독립적인 포물선 형태로 나타나지 않습니다. 즉, 유한한 파수에서 해밀토니안의 항들이 이들의 특성을 서로 혼합시킬 수 있습니다.

이를 표현하는 유용한 방법은 행렬 요소에 의해 결합된 두 개의 인접한 고유상태를 고려하는 것입니다 \(V\).

\[ H_{\mathrm{illustrative}}= \begin{pmatrix}E_{\mathrm{HH}}(k)&V(k)\\ V^{*}(k)&E_{\mathrm{LH}}(k)\end{pmatrix}. \]

만약 두 상태가 서로 결합할 수 있다면, 만나는 지점 근처에서 각각의 에너지가 서로 밀려나게 되며 HH/LH 특성도 서로 교체될 수 있습니다. 대칭성에 의해 관련 결합이 사라지는 경우에도 실제적인 에너지 교차 현상이 발생할 수 있습니다. 단지 현재 스크린샷만으로는 이러한 경우들을 구별할 수 없습니다. 상태의 특성을 확인하려면 대응하는 고유벡터를 살펴보고 필요하다면 실제 교차 지점 근처에서 더 많은 계산을 수행해야 합니다. 위에 제시된 2×2 행렬은 단지 예시일 뿐이며, 실제 8-밴드 해밀토니안 전체를 나타내지는 않습니다.

5.3 전도대 및 정공대 상태를 함께 확인하기

결합된 그래프(??)에서는 전자 상태와 정공대 상태가 동일한 에너지 축 상에 표시된다. 이 예시에서 정공대 상태들은 에너지 기준점 근처에 위치하며, 구속된 전도대 상태들은 약 1.5 eV 이상의 에너지 영역에 분포한다. 특정 두 상태 간의 에너지 차이는 수직 방향의 광학적 전이가 발생할 때의 광자 에너지가 되는데, 이는 해당 광학 행렬 요소와 캐리어 점유도에 의해 결정된다.

광자 운동량이 비교적 작기 때문에, 광학적 전이 시 평면 내 파수벡터가 보존된다는 것이 유용한 1차 근사치입니다.

\[ \hbar\omega=E_{c,n}(k_{\parallel})-E_{v,m}(k_{\parallel}), \qquad \lambda=\frac{hc}{\hbar\omega}. \]

따라서 서로 다른 전자 및 정공 상태 쌍들이 각기 다른 파장에서 기여하게 된다. 다중밴드 파동함수들은 오직 에너지 값만으로는 파악할 수 없는 각 전이의 세기와 편광 특성까지도 결정한다.

6. 양자 우물의 광학적 특성

8-띠 모델 계산을 통해 광학 모델에 필요한 전자 에너지값 및 혼합 상태의 파동함수가 산출된다. 본 예제에서는 자극 방출, 흡수, 순 이득 및 자발 방출에 대한 각각의 스펙트럼이 생성된다. 각 그래프에 나타난 여러 곡선들은 제공된 시뮬레이션에서 사용된 전하 캐리어 상태의 순서를 나타내는 것이지, 8개의 개별 해밀토니안 밴드를 의미하는 것은 아니다.

6.1 자극 방출과 흡수

자극 방출은 입사 광자가 전자를 유도하여 전도대의 채워진 상태에서 정공대의 빈 상태로 천이시키는 현상입니다. 흡수는 그 반대 과정으로, 정공대의 전자가 광자를 흡수하여 이용 가능한 전도대 상태로 이동하는 현상입니다. 이들 상태의 점유율은 개념적으로 다음과 같습니다.

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

자극 방출 곡선(??)과 흡수 곡선(??)은 기저에 존재하는 띠 간 전이 에너지가 동일하지만, 캐리어 분포의 변화에 따라 그 강도는 달라진다. 본 결과에서 대부분의 스펙트럼 구조는 790~850 nm 범위에 나타나며, 정확한 형태는 설정된 서브밴드, 행렬 요소, 캐리어 밀도 및 선폭 확장에 따라 결정된다.

6.2 순 이득: 자극 방출이 흡수보다 우세하는 경우

순 이득 그래프(??)는 자극 방출과 흡수라는 서로 상반되는 효과들을 결합한 결과입니다. 각각의 점유 확률을 서로 빼주면 매우 간단한 결과가 도출됩니다:

\[ f_c(1-f_v)-f_v(1-f_c)=f_c-f_v. \]

따라서 해당 전이 과정은 참여하는 전도대 상태의 전하 캐리어 밀도가 대응하는 정공대 상태의 밀도보다 높을 경우 순 이득을 발생시킨다. 음수 값의 순 이득은 흡수가 우세한 조건을 나타내며 양수 값은 자극 방출이 흡수를 초과하는 상황에 해당한다. 여기에 제시된 순 이득 계수는 양자 우물에서 얻어진 광학적 결과치이므로 이를 단순히 레이저 공진기 전체의 모드 이득으로 간주해서는 안 된다. 모드 이득은 구속 효과 및 기타 광학적 손실 요인에도 의존하기 때문이다.

6.3 자발 방출

자발 방출 현상은 전이를 유발하기 위해 외부 광자가 필요하지 않다. 이는 전자가 점유하고 있던 상태에서 이용 가능한 정공 상태로 방사적으로 재결합할 때 발생한다. 이때 해당 점유 인자는 다시 \(f_c(1-f_v)\)이 되지만, 자발 방출과 자극 방출은 각각 서로 다른 스펙트럼 선행 인자와 물리적 의미를 지닌다. ??에 제시된 스펙트럼들은 로그 형태의 수직 축을 사용하므로 수직 간격이 동일하다고 해서 방출되는 광자 수가 절대적으로 동일하게 증가하는 것은 아니다.

전이 강도, 캐리어 점유율 및 선폭 확장에 관한 보다 상세한 설명은 광학적 전이 및 이득 이론 페이지를 참조하시기 바랍니다.

7. 출력 파일 찾기

시뮬레이션이 완료된 후 메인 창에서 Output 항목을 선택하십시오(??). 표시되는 결과 중에는 밴드 구조를 나타내는 Ec.csv, Ev.csv, Ec_Ev.csv가 있으며, 이와 함께 gain.csv, absorption.csv, net_gain.csv, spontaneous.csv도 확인할 수 있습니다. bands 및 gain 폴더에는 추가적인 해석 결과들이 저장되어 있습니다. CSV 파일을 더블클릭하면 OghmaNano 내에서 해당 그래프가 열리며, 수치 데이터를 다른 프로그램에서 분석하기 위해 내보낼 수도 있습니다.

8. 직접 시도해보기: 포물선형 모델과 비교해보십시오

두 번째 실험으로서, 40 nm 두께의 AlGaAs 장벽층은 그대로 유지한 채 GaAs 양자 우물의 두께만 10 nm에서 8 nm로 변경합니다. 이후 모델을 다시 실행하여 구속된 고유상태들의 에너지 값과 광학적 전이가 일어나는 파장 범위를 비교합니다.