실리콘/SiGe 양자 우물: 6밴드 루팅거-콘 모델을 이용한 정공 거동 시뮬레이션 튜토리얼
1. 서론
우리의 실리콘 Δ-밴드 튜토리얼에서는 Si/SiGe 양자 우물 내 전자들의 거동을 살펴보았습니다. 본 튜토리얼에서는 반대로 실리콘의 전자 구조 중 정공대에 존재하는 홀들의 거동을 다룹니다. 층 구조는 유사하지만 그 기저가 되는 물리 현상은 상당히 다릅니다. 개별 실리콘 전도대 근처에 있는 전자는 이방성 유효질량을 이용해 기술할 수 있으나, 정공대의 최고점에는 서로 상호작용하는 여러 상태들이 존재합니다. 따라서 단일 홀 분산 곡선만으로는 중요한 정보를 놓칠 수밖에 없습니다.
본 예제에서는 OghmaNano의 6밴드 루팅거-콘 모델(LK6)을 활용하여 두 개의 20nm 두께 SiGe 층 사이에 위치한 10nm 두께 실리콘 양자 우물 내에서 구속된 홀 서브밴드를 계산합니다. 이후 밴드 에너지 계산을 넘어 음향 포논 산란 현상, 홀의 분포 상태 및 이에 따른 면내 전하 캐리어 이동도까지 평가합니다. 사용자는 구조를 처음부터 새로 구성하지 않고도 제공된 기본 예제를 그대로 활용할 수 있습니다.
2. 6밴드 예제 열기
파일 → 새 시뮬레이션을 선택합니다. 예제 브라우저에서 GaAs demos 항목을 더블클릭하세요(??). 이 폴더에는 실리콘 예제를 포함한 반도체 양자 우물 템플릿들도 포함되어 있습니다. 다음 창에서 Strained Si/SiGe Valence-Band Quantum Well (LK6) 항목을 더블클릭하세요(??). 반드시 정공대 예제를 선택해야 합니다. 인접한 Delta-Valley 템플릿은 실리콘의 전도 전자를 계산하는 용도이기 때문입니다.
3. Si/SiGe 구조 및 재료 파라미터 확인
메인 창에는 예시 구조물의 3차원 표현이 나타나 있습니다(??). 전하 캐리어들은 z 방향으로는 구속 효과를 받아 움직임이 제한되지만 x-y 평면에서는 자유롭게 이동할 수 있습니다. ‘레이어 편집기’를 클릭하면 각 층의 실제 두께와 중앙 층이 양자 우물로 지정된 내용을 확인할 수 있습니다.
| 층 | 두께 | 역할 |
|---|---|---|
| 실리콘-게르마늄 | 20 나노미터 | 주변층 |
| 실리콘(양자 우물) | 10 나노미터 | 구속된 정공대 상태 |
| 실리콘-게르마늄 | 20 나노미터 | 주변층 |
??에 제시된 광학-물질 열은 전자 물질 정의 내용을 확인하는 대체 수단이 아닙니다. 밴드 가장자리 오프셋, 합금 조성 및 응력 설정 모두 실제 구속 효과에 영향을 미치므로 초기 시뮬레이션 시 예제에서 제공된 설정값을 그대로 사용하시기 바랍니다.
전기적 파라미터 편집기를 열고 Si (QW)를 선택하십시오(??). 스크린샷에는 루팅거 파라미터 gamma1와 그에 대한 문헌 참고 자료가 표시되어 있습니다. 단일 밴드 전자 모델과는 달리, 정공대 계산에서는 서로 결합된 중정공, 경정공 및 스플릿-오프 상태를 기술하기 위해 세 개의 루팅거 파라미터가 사용됩니다. 또한 스핀-궤도 분열 값이 필요하며, 응력이 적용되는 경우에는 적절한 변형 전위 파라미터도 필요합니다.
4. 왜 정공대 모델에는 여섯 개의 성분이 필요한가요?
실리콘의 정공대 최대값 근처에는 홀 상태로서 잘 알려진 두 가지 계열이 존재하는데, 바로 무거운 홀(HH)과 가벼운 홀(LH)입니다. 스핀-궤도 상호작용에 의해 더 낮은 에너지 위치에 나타나는 분리된 계열(SO)도 존재합니다. 각 계열마다 스핀과 관련된 두 개의 기저 상태가 존재하므로 루팅거-콘 6-밴드 모델에서는 총 6개의 구성 요소가 생성됩니다. 이들은 서로 결합된 파동함수의 기저 성분일 뿐, 반드시 독립적으로 움직이는 6개의 포물선형 밴드를 의미하는 것은 아닙니다.
양자 우물의 파동함수는 6개 성분으로 구성된 열벡터로서 개념적으로 표현할 수 있다.
\[ \Psi(z)=\begin{pmatrix} F_{\mathrm{HH}+}(z) & F_{\mathrm{LH}+}(z) & F_{\mathrm{SO}+}(z) & F_{\mathrm{HH}-}(z) & F_{\mathrm{LH}-}(z) & F_{\mathrm{SO}-}(z) \end{pmatrix}^{\!T}. \]
여기서 각 레이블은 기저가 되는 각운동량의 주요 특성을 나타냅니다. 평면 내 운동량이 변함에 따라 각 성분들은 일반적으로 혼합되므로, 처음에는 중량홀형으로 기술된 고유상태라도 더 큰 파동벡터에서는 반드시 순수한 중량홀형 상태로 남아있지 않을 수 있습니다. 컴퓨터가 해결하는 문제의 형태는 다음과 같습니다
\[ \left[H_{\mathrm{LK}}(k_x,k_y,-i\partial_z) +H_{\mathrm{BP}}(\varepsilon)+V_v(z)I_6\right]\Psi =E(k_x,k_y)\Psi. \]
루팅거-콘 항은 정공대 상태들 간의 운동 에너지 및 상호작용을 기술하며, Bir–Pikus 항은 응력에 따른 이러한 상태들의 반응을 설명한다. 또한 \(V_v(z)\)은 정공대의 구속 효과 프로파일을 나타낸다. 이 모델에서는 HH, LH, SO 상태들에 각각 독립적인 질량값을 부여하는 대신 기저 성분들 간의 결합을 그대로 유지한다. 본 튜토리얼에서는 사용자가 직접 이 행렬을 구성할 필요가 없으며, 상세한 해밀토니안 정보는 the LK6 theory page를 참조하면 된다.
비교를 위해, 단순화된 벌크 근사법만으로도 세 개의 루팅거 파라미터가 왜 중요한지를 알 수 있습니다. 입방 결정의 주축 방향에서 정공과 경공의 에너지 곡률은 각각 해당 파라미터들의 서로 다른 조합과 관련되며, 이는 개념적으로 다음과 같습니다.
\[ \frac{m_0}{m_{\mathrm{HH}}^{*}}\simeq\gamma_1-2\gamma_2, \qquad \frac{m_0}{m_{\mathrm{LH}}^{*}}\simeq\gamma_1+2\gamma_2. \]
이러한 관계식들은 설명용으로 유용한 극한 상황을 나타내는 것일 뿐, 6개 밴드로 구성된 양자 우물 시뮬레이션을 두 개의 별개 포물선 모델로 대체하라는 지시가 아닙니다. 구속 효과를 받는 구조에서는 우물의 폭, 스핀-궤도 결합, 응력 및 면내 운동량이 실제 고유값과 그 상태 특성의 변화를 결정합니다.
5. 계산 수행하기
파일 리본으로 돌아가서 ‘시뮬레이션 실행’ 버튼을 누르세요. OghmaNano는 먼저 여러 평면 내 파동벡터에 대해 구속 효과를 받는 루팅거-콘 모델의 정공대 고유상태들을 계산합니다. 터미널 탭(??)에서 계산 진행 상황을 확인할 수 있습니다. 표시된 target 에너지 값들은 내부 고유값 해법기의 목표값이며, 출력된 수렴 메시지들은 요청된 고유상태들이 수치적으로 산출되었음을 나타냅니다.
솔버가 계산을 완료할 때까지 기다리십시오. 산란 과정은 단순히 서브밴드 에너지만을 구하는 것이 아니라, 점유된 초기 전자 상태와 가능한 최종 전자 상태 간의 전이를 평가해야 하므로 추가적인 시간이 소요됩니다.
6. 계산된 정공대 분산 특성 확인하기
출력 탭을 열어보세요. 생성된 파일들에는 밴드 관련 그래프와 더불어 산란 및 이동도 관련 폴더들이 포함되어 있습니다(??). 정공대 결과 파일을 열면 ??를 확인할 수 있습니다.
\(k_{\parallel}=0\) 지점에서는 구속 효과에 의해 각 서브밴드의 시작 에너지가 이산적으로 결정됩니다. \(k_{\parallel}\) 지점으로 이동하면, 표시된 전자 에너지 체계 하에서 에너지 분포 곡선들이 아래로 휘어지며 곡률과 서브밴드 간 간격도 다소 달라집니다. 서로 가까이 위치한 선들은 시간역전 대칭성을 갖는 계산에서 나타나는 스핀 관련(크레이머스) 쌍들과 일치하며, 이들을 별개의 홀 종류로 간주해서는 안 됩니다.
이 그래프를 이해하는 데 유용한 방법은 정공이란 정공대 내에서 전자가 없는 상태임을 기억하는 것입니다. 따라서 정공 에너지는 선택된 정공대 기준점을 기준으로 표현될 수 있습니다.
\[ E_h(\mathbf{k}_{\parallel})=E_{v,\mathrm{ref}}-E_v(\mathbf{k}_{\parallel}). \]
따라서 최대값에서부터 아래쪽으로 휘어지는 정공대 분지는 정공의 여기 에너지가 증가함을 의미한다. 각 분지의 실제 정공 특성인 중량정공/경량정공/분리된정공 성분은 오직 색상이나 곡률만으로 결정되는 것이 아니라 여섯 개 성분으로 구성된 고유벡터에 의해 결정된다.
이 결과를 실리콘 Δ-밸리 튜토리얼에 제시된 전도대 분산과 비교해보는 것이 유용합니다. 해당 튜토리얼에서는 서로 다른 밸리 방향들을 이방성 전자 질량을 이용해 다루었으나, 본 모델에서는 여러 정공 상태들이 다중밴드 정공 해밀토니안 내에서 직접적으로 결합합니다.
7. 음향 포논 산란과 홀의 수명
밴드 구조를 통해 어떤 정공 상태들이 존재하는지는 알 수 있지만, 이만으로는 정공이 얼마나 빨리 운동량을 잃게 되는지를 예측할 수는 없다. 본 예제에서는 OghmaNano가 계산된 구속 상태들을 활용하여 음향 포논 산란 현상을 평가한다. 음향 포논은 결정 내에 미세한 응력장을 생성함으로써 정공대에 교란을 일으키며, 이로 인해 정공이 초기 상태 \(i\)에서 최종 상태 \(f\)로 이동할 가능성이 생긴다.
si_lk6_acoustic_rate_*.csv, si_lk6_acoustic_tau_*.csv 및 전송 관련 데이터셋들이 포함되어 있습니다.Output 탭에서 산란 관련 폴더(??)를 열어보세요. 번호가 매겨진 파일들은 계산된 고유상태들을 나타냅니다. 예를 들어, si_lk6_acoustic_tau_02.csv에 첨부된 스크린샷은 음향 포논 산란에 의한 수명값을 평면 내 파동벡터에 대해 그래프로 나타낸 것입니다.
물리적 기초는 페르미의 황금법칙입니다.
\[ W_{i\rightarrow f}=\frac{2\pi}{\hbar} \left|\langle\Psi_f|H_{\mathrm{e-ph}}|\Psi_i\rangle\right|^2 \delta(E_f-E_i\pm\hbar\omega_{\mathbf q}). \]
해당 행렬 요소는 격자 진동이 두 개의 홀 고유상태 간에 얼마나 강하게 결합하는지를 나타내며, 디랙 델타 함수는 포논 흡수 또는 방출 과정에서의 에너지 보존을 표현한다. 6-밴드 정공대 모델에서는 이러한 결합이 다성분 홀 파동함수에 의존하므로 허용되는 전이 과정은 HH/LH/SO 상태 간의 혼합 및 이용 가능한 최종 에너지 수준에 의해 영향을 받는다. 본 출력 결과는 이 예제에서 사용된 음향 포논 근사법에 관한 내용이다.
si_lk6_acoustic_tau_02.csv에서 얻은 음향 포논 산란에 의한 수명값 예시입니다. 이 값이 파동벡터에 따라 변하는 현상은 가능한 최종 상태들과 그들 간의 결합 정도를 반영합니다.허용되는 전이들을 모두 합산하면 총 사건 발생률 \(\Gamma_i\)이 산출되며 이에 상응하는 수명값도 얻어진다.
\[ \Gamma_i=\sum_f W_{i\rightarrow f},\qquad \tau_i=\frac{1}{\Gamma_i}\quad(\Gamma_i\gt 0). \]
??에 표시된 수명값은 일정하지 않습니다. 이 값은 파동벡터가 결합된 정공대 분산곡선의 여러 영역을 지날 때마다 변화합니다. 첫 번째 샘플링 지점에서 나타나는 0값은 물리적으로 산란 수명이 정확히 \(k_{\parallel}=0\)에서 0이 된다는 증거가 아니라, 경계 조건 또는 출력 형식상의 표기이거나 검토가 필요한 값으로 간주되어야 합니다.
수송 현상에서는 또 하나의 중요한 구분점이 존재한다. 모든 충돌이 전류를 동일하게 완화시키는 것은 아니다. 소각도 산란의 경우 미시적 상태는 변화할 수 있으나 전체 운동 방향에 미치는 영향은 상대적으로 미미하다. 따라서 수송 가중화된 산란율은 다음과 같이 개략적으로 표현된다.
\[ \Gamma_{\mathrm{tr},i}=\sum_f W_{i\rightarrow f} \left(1-\cos\theta_{if}\right),\qquad \tau_{\mathrm{tr},i}=\frac{1}{\Gamma_{\mathrm{tr},i}}. \]
각도 관련 식은 운동량 완화 원리를 나타내며, 이방성 혼합 밴드의 경우 전체 수송 계산 시 적절한 속도값과 가중치를 반드시 사용해야 한다. 이러한 이유로 출력 결과에서는 음향 포논 산란률 및 수명값이 수송 관련 물리량들과 별도로 제시되며, 모든 용도에 단일 수치만을 적용하지 않는다.
8. 홀 밀도와 페르미 준위로부터 전하 캐리어 이동도 계산하기
측정 가능한 수송 특성을 얻기 위해서는 산란에 관한 정보와 각 서브밴드 내에서 홀이 얼마나 많이 존재하는지를 결합해야 합니다. 따라서 본 튜토리얼에서는 단위가 m−2인 홀 면밀도에 대한 스윕 실험을 포함하고 있습니다 \(p_{2D}\).
si_lk6_hole_fermi_level.csv에서 측정된 홀 페르미 준위의 변화.
si_lk6_hole_mobility.csv을 통해 계산된 홀의 전하 캐리어 이동도.페르미 준위는 어떤 전자 상태가 채워지는지를 결정합니다. 정공대 내의 특정 에너지값 \(E_v\)에 대해 해당 홀 상태의 평형 상태 점유율은 다음과 같이 표현됩니다.
\[ f_h(E_v)=1-f_e(E_v) =\frac{1}{1+\exp\!\left[(E_F-E_v)/(k_B T)\right]}. \]
홀 밀도를 증가시키면 화학적 퍼텐셜이 정공대 상태들에 대해 이동하게 되는데, 이는 왼쪽 그래프에서 보여주는 현상입니다. 실제 수치적 에너지 값은 모델에서 채택한 기준점에 따라 달라지지만, 에너지 스캔의 방향과 그에 따른 홀 분포가 가장 중요한 특징입니다.
OghmaNano는 채워진 상태들의 속도와 운송 완화 시간에 가중치를 부여하여 면내 홀 전도도를 산출한다 \(\sigma_{2D}\). 이에 따른 홀 이동도는
\[ \mu_h=\frac{\sigma_{2D}}{q\,p_{2D}}, \qquad \sigma_{2D}\sim q^2\sum_n\int \frac{d^2 k_{\parallel}}{(2\pi)^2} \,v_{n,x}^{2}\,\tau_{\mathrm{tr},n} \left(-\frac{\partial f_e}{\partial E}\right). \]
여기서 \(v_{n,x}\)은 서브밴드 \(n\)에서의 면내 그룹 속도를 나타내며, 적분식은 서로 다른 파동벡터를 갖는 상태들이 전류에 각기 다르게 기여함을 의미합니다. 전하 캐리어 이동도 그래프(??)를 보면 표시된 밀도 범위 내에서 이동도가 약간 감소하는 경향을 보이는데, 그 값은 대략 0.165에서 0.151 m2 V−1 s−1 사이입니다. 이 결과는 설정된 예제와 그에 포함된 산란 물리 모델에 따른 결과일 뿐이며, 실리콘 내 홀의 이동도에 관한 보편적 법칙이나 모든 실험적 손실 메커니즘을 반영한 예측값은 아닙니다.
9. 직접 시도해보기
💡 권장 실험 방법: 기본 예제 파일을 복사한 뒤 한 번에 하나의 변수만 변경해 보십시오.
- Si 양자 우물의 폭을 10 nm에서 변경한 뒤 계산을 다시 수행하십시오. 구속된 정공대 에너지들을 비교해보면, 우물 폭이 좁아질수록 양자 구속 효과가 커지는 경향이 있으나 HH/LH/SO 상태들의 상호작용에 따른 반응은 단일 질량 모델의 경우보다 훨씬 복잡합니다.
- 홀 모델과 전자 모델을 비교해보세요: 유사한 층 구조를 사용하여 연관된 실리콘 Δ-밴드 예제를 실행하십시오. 정공 해석기에서는 서로 결합된 6개 성분의 홀 상태가 기술되는 반면, 전도대 계산에서는 실리콘의 다양한 밴드 방향들이 구별된다는 점을 확인할 수 있습니다.
si_lk6_acoustic_tau_00.csv과 같은 다른 음향 포논 산란 수명 파일을 열어서si_lk6_acoustic_tau_02.csv와 비교해보십시오. 서로 다른 서브밴드들은 동일한 산란 경로를 갖지 않을 수도 있습니다.- 밀도 변화를 살펴보십시오: 단일의 유효 산란 수명값이 두 곡선 모두를 예측한다고 가정하지 말고, 페르미 준위와 전하 캐리어 이동도의 변화를 비교해보십시오.
이로써 6개 밴드 모델을 다루는 양자 우물 관련 입문 시리즈가 마무리됩니다. 수식적 전개와 재료 파라미터의 의미에 대해서는 6개 밴드 루팅거-콘 이론 장 및 버-피쿠스 응력 장에서 상세히 설명하고 있습니다.