프뢰리히 산란: 반도체 내 극성 광학 포논과의 상호작용

프뢰리히 산란은 GaAs, InP, GaN, ZnO 및 금속-할라이드 페로브스카이트와 같은 극성 반도체 내에서 전자와 정공이 종방향 광학(LO) 포논과 상호작용하는 현상을 설명한다. 이러한 상호작용의 물리적 기원은 ??에 도시되어 있다. 극성 결정 구조에서는 인접한 원자들이 서로 반대 부호의 부분 전하를 지닌다. 결정 격자가 광학 모드로 진동할 때 양이온과 음이온이 서로 상대적으로 움직이면서 진동하는 전기 분극이 발생한다. 종방향 광학 포논의 경우 이 현상은 거시적인 전기장을 생성하게 되며, 이 전기장은 결정 내부를 이동하는 전하 캐리어들과 상호작용하게 된다.

이 전기장을 통과하는 전자는 초기 상태의 파동벡터 \(\mathbf{k}\)에서 최종 상태 \(\mathbf{k}'\)로 산란될 수 있으며, 이때 에너지가 \(\hbar\omega_{LO}\)인 LO 포논을 흡수하거나 방출하게 된다. 이러한 전이 과정들은 ??에 제시된 전도대 다이어그램에서 설명된다. 자유 전하 캐리어와 결정격자의 분극장 사이의 정전기적 상호작용을 전하-분극 상호작용이라 부른다. 이 현상은 극성 반도체 내에서 특히 상온 이상에서 전하 캐리어 이동도를 제한하는 중요한 메커니즘이다. 이는 고순도 III-V족 반도체에서 전자의 산란 현상을 지배하며, 고에너지 상태의 캐리어가 결정격자로 에너지를 전달하는 효율적인 경로이기도 하다.

프뢰리히 산란 현상을 이해하는 것은 전하 캐리어 이동도 계산, 광여기된 캐리어의 냉각 현상 모델링, 그리고 에너지 균형 기반 수송 모델에서 사용되는 에너지 완화 시간 결정에 있어 매우 중요하다. 이 상호작용의 에너지 의존성은 ??에서 확인할 수 있으며, 여기서는 GaAs 내 전자에 대한 계산된 흡수율, 방출율 및 총 산란율이 나타나 있다. 이후 절들에서는 LO 포논의 흡수 및 방출 현상을 보다 상세히 살펴본 뒤, 프뢰리히 행렬 요소와 산란율 계산에 활용되는 3차원 \(k\) 공간 적분식에 대해 설명한다.

프뢰리히 상호작용이란 무엇인가?

극성 반도체 내에서 종방향 광학(LO) 진동이 발생할 때 양이온과 음이온의 상대적 변위로 인해 진동하는 전기장이 생성되며, 이는 ??에 도시된 바와 같이 \(\mathbf{E}_{LO}\)을 형성한다. 전하 캐리어는 이 전기장과 전기적 퍼텐셜 \(\phi_{LO}\)을 통해 상호작용하며, 전자의 경우 상호작용 에너지는 \(-e\phi_{LO}\)이 된다. 이 상호작용은 쿨롱 힘에 기인하므로 장거리에 걸쳐 작용하며, 특히 작은 포논 파벡터 영역에서 가장 강하게 나타난다 \(q\).

프뢰리히 상호작용의 강도는 결정체가 전기장에 얼마나 강하게 반응하는지에 따라 결정된다. 이러한 반응에는 두 가지 요소가 기여하는데, 매우 빠르게 반응하는 전자들과 그에 비해 훨씬 느리게 반응하는 이온들이 있다. 이로 인해 두 가지 유전 상수가 생성되는데, 전자 반응만을 기술하는 고주파 유전 상수 \(\varepsilon_\infty\)와 이온 운동의 영향까지 포함하는 정적 유전 상수 \(\varepsilon_s\)가 그것이다.

극성 반도체 내에서 LO 포논에 의해 이온들이 변위하면 전기장이 생성되며 이 전기장은 자유 전하 캐리어들과 상호작용하게 됩니다. 이러한 상호작용의 강도는 전자적 반응만을 나타내는 고주파 유전 상수 \(\varepsilon_\infty\)와 이온 운동의 영향까지 포함한 정적 유전 상수 \(\varepsilon_s\) 간의 차이에 의해 결정됩니다.

\[ \frac{1}{\varepsilon_\infty} - \frac{1}{\varepsilon_s}. \]

이 값이 클수록 극성 결합의 강도가 더욱 커진다. 실리콘과 같은 비극성 반도체에서는 광학적 격자 진동에 따른 장거리 전기장이 존재하지 않는다. 따라서 이 모델에서는 유전 반응에 대한 이온 기여도가 없어 \(\varepsilon_s=\varepsilon_\infty\)이 발생하며 프뢰리히 전하-분극 상호작용도 소멸하게 된다. 이 유전 항은 아래에서 도출되는 산란 행렬 요소에 직접적으로 포함된다.

결합 강도를 단일 무차원 수인 프뢰리히 결합 상수로 표현하는 것이 종종 편리합니다.

\[ \alpha = \frac{e^2}{4\pi\varepsilon_0 \hbar}\sqrt{\frac{m^*}{2\hbar\omega_{LO}}}\left(\frac{1}{\varepsilon_\infty}-\frac{1}{\varepsilon_s}\right), \]

여기서 \(m^*\)은 전하 캐리어의 유효질량입니다. GaAs 내 전자의 경우 \(\alpha\approx 0.07\) 값은 약간의 결합만을 나타내며; GaN, ZnO 및 할라이드 페로스카이트와 같은 이온성이 더 강한 물질에서는 \(\alpha\) 값이 0.4에서 2 사이입니다. \(\alpha\) 값이 작을 때는 이 상호작용을 섭동 이론으로 다룰 수 있으며, 이를 통해 아래에 기술된 산란률을 얻게 됩니다. 반면 \(\alpha\) 값이 1에 가깝거나 그 이상일 경우 캐리어는 상당한 양의 분극 구름을 함께 끌고 다니게 되어 폴라론으로 간주하는 것이 적절하며, 이 경우 재정규화된 질량 및 에너지를 갖게 됩니다. 그러므로 이 페이지에 제시된 섭동 이론 기반 산란률은 이러한 경우에는 신중하게 사용되어야 합니다.

밴드 구조에서의 흡수 및 방출 현상

프뢰리히 산란은 전자가 LO 포논과 에너지 및 운동량을 교환하는 비탄성적 전이 현상입니다. ??에 나타난 바와 같이, 초기에 에너지가 \(E(\mathbf{k})\)인 상태 \(\mathbf{k}\)에 있는 전자는 다음과 같은 두 가지 과정 중 하나를 겪을 수 있습니다:

포논 흡수: 전자는 격자로부터 에너지 \(\hbar\omega_{LO}\)를 얻어 더 높은 에너지 상태인 \(\mathbf{k}'=\mathbf{k}+\mathbf{q}\)로 이동하며, 이때 \(E(\mathbf{k}')=E(\mathbf{k})+\hbar\omega_{LO}\)가 된다.

포논 방출: 전자는 에너지를 \(\hbar\omega_{LO}\) 결정격자로 전달하여 더 낮은 에너지 상태인 \(\mathbf{k}'=\mathbf{k}-\mathbf{q}\)로 이동하며, 이때 \(E(\mathbf{k}')=E(\mathbf{k})-\hbar\omega_{LO}\)가 된다.

두 과정 모두 에너지와 결정격자 운동량이 보존된다. 영역 중심 근처에서 LO 포논은 거의 분산이 없으므로 그 에너지 \(\hbar\omega_{LO}\)는 일정한 값으로 간주될 수 있으며, GaAs에서는 약 36 meV, GaN에서는 약 92 meV이다. 이는 포논 방출 과정에 중요한 제약을 부여하는데, 전자가 LO 포논을 방출하려면 반드시 밴드 최저점보다 최소한 \(\hbar\omega_{LO}\)만큼의 운동 에너지를 보유해야만 한다. 이 임계값 이하의 전자는 오직 포논 흡수만을 겪을 수 있다. 임계값을 초과하면 방출 과정이 가능해져 고에너지 전자들이 \(\hbar\omega_{LO}\) 크기만큼의 에너지를 격자에 방출할 수 있게 된다. 이러한 현상은 GaAs 내 전자의 흡수율, 방출율 및 총 산란율을 나타낸 ??에서 명확히 확인할 수 있다.

각 과정의 발생 확률은 격자 내에 존재하는 LO 포논의 수에도 의존합니다. 열평형 상태에서는 평균 포논 점유수가 보스-아인슈타인 분포에 의해 기술됩니다.

\[ N_{LO} = \frac{1}{\exp\!\left(\hbar\omega_{LO}/k_B T_L\right)-1}, \]

여기서 \(T_L\)은 격자 온도를 나타냅니다. 흡수 확률은 \(N_{LO}\)에 비례하는데, 이는 전자가 포논을 흡수하려면 해당 포논이 이미 존재해야 하기 때문입니다. 방출 현상은 \(N_{LO}+1\)에 비례합니다. 추가적으로 \(+1\)는 자발적 포논 방출을 의미하는데, 이는 전자가 충분한 에너지를 지니고 있는 한 영온도에서도 발생할 수 있습니다.

300K 온도에서의 GaAs의 경우 \(N_{LO}\approx 0.33\)에 따라, 방출에 의한 포논 점유 인자는 흡수 시의 값보다 약 4배 높게 나타납니다. 실제 산란률은 ??에서 보듯이 이용 가능한 최종 전자 상태에도 의존합니다. 격자 온도가 상승하면 포논의 수가 증가하여 전하-분극 상호작용이 강화되며, 이는 극성 반도체에서 일반적으로 관찰되는 전하 캐리어 이동도의 감소 현상에 기여합니다.

프뢰리히 결합 행렬 요소

이제 양자역학적 관점에서 이러한 상호작용을 고찰해보겠습니다. 전하 캐리어와 종방향 광학 포논 간의 결합은 프뢰리히 행렬 요소로 기술되며, 극성 반도체의 경우 이는 다음과 같이 표현됩니다.

\[ M(q) = -\frac{i}{q} \left[ \frac{e^2\hbar\omega_{LO}}{2V\varepsilon_0} \left( \frac{1}{\varepsilon_\infty} - \frac{1}{\varepsilon_s} \right) \right]^{1/2}. \]

여기서 \(V\)은 결정 정규화 부피이며 \(q\)는 포논의 파동벡터입니다. 전자 상태를 평면파로 가정할 경우, 이 상호작용은 결정 운동량 보존 조건에 따라 초기 상태 \(|\mathbf{k}\rangle\)과 최종 상태 \(|\mathbf{k}'\rangle\)를 연결합니다. 포논 흡수 및 방출 과정에 따라 최종 파동벡터는 각각 \(\mathbf{k}'=\mathbf{k}\pm\mathbf{q}\)가 됩니다.

포논 분포 인자를 고려하고, 초기 및 최종 블로흐 함수의 주기적 성분들 간의 중첩이 대략 1에 가깝다고 가정할 경우, 전이 행렬 요소의 제곱값은 다음과 같이 표현된다.

\[ |M_{\mathrm{abs}}(q)|^2 = |M(q)|^2 N_{LO}, \]

\[ |M_{\mathrm{em}}(q)|^2 = |M(q)|^2 (N_{LO}+1). \]

흡수 항은 존재하는 LO 포논의 수에 의존하는 반면, 방사 항에는 자발적 방사에 의한 추가적인 기여분이 포함된다. 이러한 행렬 요소들은 페르미의 황금법칙에서 전자 상태 간 전이 확률을 계산하는 데 활용된다.

페르미의 황금률에서 3차원 k-공간 적분까지

프뢰리히 결합 상수 행렬 요소를 도출했으므로, 이제 전자가 초기 상태 \(|\mathbf{k}\rangle\)에서 최종 상태 \(|\mathbf{k}'\rangle\)로 산란하는 비율을 계산할 수 있다. ??에 나타난 바와 같이, 이러한 산란은 LO 포논의 흡수 또는 방출 과정을 통해 일어나며, 그 결과 전자는 에너지 \(\hbar\omega_{LO}\)를 얻거나 잃게 된다.

페르미의 황금률을 적용하면, 모든 허용되는 포논 파벡터에 대해 합산하고 두 가지 산란 과정을 모두 고려함으로써 총 산란률을 구할 수 있다 \(\mathbf{q}\).

\[ \begin{aligned} \frac{1}{\tau(\mathbf{k})} = \frac{2\pi}{\hbar} \sum_{\mathbf{q}} \Big[ &|M_{\mathrm{abs}}(q)|^2 \delta\big( E(\mathbf{k}+\mathbf{q})-E(\mathbf{k})-\hbar\omega_{LO} \big) \\ +{}& |M_{\mathrm{em}}(q)|^2 \delta\big( E(\mathbf{k}-\mathbf{q})-E(\mathbf{k})+\hbar\omega_{LO} \big) \Big]. \end{aligned} \]

여기서 \(\delta\)은 에너지 보존을 만족시키는 디랙 델타 함수입니다. 흡수 및 방출 과정에 관여하는 행렬 요소인 \(|M_{\mathrm{abs}}(q)|^2=|M(q)|^2N_{LO}\)와 \(|M_{\mathrm{em}}(q)|^2=|M(q)|^2(N_{LO}+1)\)는 이전 절에서 소개되었습니다. 최종 전자 상태들이 비점유 상태라고 가정하므로 파울리 차단 효과는 무시할 수 있습니다.

충분히 큰 크기의 결정체의 경우 허용되는 포논의 파벡터들이 서로 매우 근접해 있으므로, 이 합산 과정은 브릴루앙 영역에 대한 적분으로 대체될 수 있다.

\[ \sum_{\mathbf{q}} \longrightarrow \frac{V}{(2\pi)^3}\int d^3q. \]

\(V\) 인수는 제곱된 프뢰리히 행렬 요소에서 나타나는 \(1/V\) 의존성을 상쇄시켜, 결정체 정규화 부피와 무관한 산란률을 산출한다. \(M_{\mathrm{abs}}(q)\)에 대한 식을 대입하면 흡수률이 다음과 같이 표현된다.

\[ \begin{aligned} \frac{1}{\tau_{\mathrm{abs}}(\mathbf{k})} ={}& \frac{2\pi}{\hbar} \frac{V}{(2\pi)^3} \int d^3q\; \frac{e^2\hbar\omega_{LO}}{2V\varepsilon_0} \left( \frac{1}{\varepsilon_\infty} - \frac{1}{\varepsilon_s} \right) \frac{N_{LO}}{q^2} \\[4pt] &\times \delta\big( E(\mathbf{k}+\mathbf{q})-E(\mathbf{k})-\hbar\omega_{LO} \big). \end{aligned} \]

이 적분값을 구하기 위해 초기 전자 파동벡터 방향을 극축으로 하는 구면좌표계를 사용한다 \(\mathbf{k}\). 방위각에 대해 적분하면

\[ d^3q = 2\pi q^2\,dq\,d(\cos\theta), \]

여기서 \(\theta\)은 전자와 포논의 파동벡터 사이의 각도를 나타내며, \(\mathbf{k}\)와 \(\mathbf{q}\)도 마찬가지입니다. 주목할 점은 부피 요소에서 기인하는 \(q^2\) 항이 프뢰리히 행렬 요소의 \(1/q^2\) 의존성을 상쇄한다는 것입니다. 이로 인해 남은 적분 과정이 상당히 단순화됩니다.

포물선형이며 등방성인 전도대를 가정할 때, \(E(\mathbf{k})=\hbar^2k^2/2m^*\) 흡수 델타 함수의 인수를 전개하면 다음과 같다.

\[ \begin{aligned} &E(\mathbf{k}+\mathbf{q})-E(\mathbf{k})-\hbar\omega_{LO} \\[4pt] &\qquad = \frac{\hbar^2}{2m^*} \left( 2kq\cos\theta+q^2 \right) -\hbar\omega_{LO}. \end{aligned} \]

\(\cos\theta\)에 대한 적분을 수행하면 델타 함수가 제거되며, 해가 \(-1\leq\cos\theta\leq1\)을 만족할 경우 \(m^*/(\hbar^2kq)\)라는 인수가 추가된다. 이 조건은 에너지 보존 및 결정 격자 운동량 보존에 따라 허용되는 포논 파벡터의 범위를 제한한다.

포논 흡수의 경우 허용되는 범위는 다음과 같다

\[ \begin{aligned} q_{\min} &= k'-k,\\ q_{\max} &= k'+k,\\ k' &= \sqrt{k^2+\frac{2m^*\omega_{LO}}{\hbar}}. \end{aligned} \]

포논 방출에는 \(E\geq\hbar\omega_{LO}\)이 필요하며, 이에 따른 한계값은 다음과 같다.

\[ \begin{aligned} q_{\min} &= k-k',\\ q_{\max} &= k+k',\\ k' &= \sqrt{k^2-\frac{2m^*\omega_{LO}}{\hbar}}. \end{aligned} \]

남은 방사형 적분 과정은 간단합니다:

\[ \int_{q_{\min}}^{q_{\max}} \frac{dq}{q} = \ln\left(\frac{q_{\max}}{q_{\min}}\right). \]

각 항들을 모아서 \(\hbar k=\sqrt{2m^*E}\)을 활용하면, 여기서 \(E\)는 밴드 최솟값을 기준으로 측정된 전하 캐리어의 운동에너지이므로 프뢰리히 산란률에 대한 표준적인 폐형식 표현식들을 얻게 됩니다. 흡수 과정의 경우

\[ \boxed{ \begin{aligned} \frac{1}{\tau_{\mathrm{abs}}(E)} ={}& C(E)N_{LO} \\ &\times \ln\left| \frac{ \sqrt{E+\hbar\omega_{LO}}+\sqrt{E} }{ \sqrt{E+\hbar\omega_{LO}}-\sqrt{E} } \right|, \end{aligned} } \]

반면에 방출의 경우에는

\[ \boxed{ \begin{aligned} \frac{1}{\tau_{\mathrm{em}}(E)} ={}& C(E)(N_{LO}+1) \\ &\times \ln\left| \frac{ \sqrt{E}+\sqrt{E-\hbar\omega_{LO}} }{ \sqrt{E}-\sqrt{E-\hbar\omega_{LO}} } \right|, \qquad E\geq\hbar\omega_{LO}. \end{aligned} } \]

LO 포논 에너지 임계값보다 낮은 경우 방출률은 0이 됩니다. 위 식들에서 공통적으로 나타나는 전구항은 다음과 같습니다.

\[ C(E)= \frac{e^2\omega_{LO}}{4\pi\varepsilon_0\hbar} \left( \frac{1}{\varepsilon_\infty} - \frac{1}{\varepsilon_s} \right) \sqrt{\frac{m^*}{2E}}. \]

전체 프뢰리히 방출 산란률은 두 가지 기여분의 합으로 간단히 표현됩니다.

\[ \boxed{ \frac{1}{\tau(E)} = \frac{1}{\tau_{\mathrm{abs}}(E)} + \frac{1}{\tau_{\mathrm{em}}(E)}. } \]

이 공통 선행인수는 무차원 프뢰리히 결합 상수 \(\alpha\)를 이용해 표현할 수 있으며, 그 결과는 \(C(E)=\alpha\omega_{LO}\sqrt{\hbar\omega_{LO}/E}\)와 같다. 두 가지 한계 상황을 주목할 필요가 있다. \(E\rightarrow0\)일 때 흡수율은 유한한 값인 \(2\alpha\omega_{LO}N_{LO}\)에 수렴한다. 한편 방출 임계값인 \(E=\hbar\omega_{LO}\)에서는 방출 기여도가 0에서 시작하여 캐리어 에너지가 증가함에 따라 급격히 상승한다.

이러한 현상은 ??에서 직접 확인할 수 있으며, 이 그림은 300 K 온도에서 GaAs 내 전자들의 흡수율, 방출율 및 총 산란율을 나타낸다. 고에너지 전자의 경우 산란율은 일반적으로 \(10^{12}\)–\(10^{13}\,\mathrm{s^{-1}}\) 수준에 달하며, LO 포논 방출에 소요되는 시간은 약 0.1–0.2 ps에 이른다. 이러한 빠른 방출 과정은 고에너지 전하 캐리어가 격자에 에너지를 전달하는 중요한 메커니즘으로 작용한다.

전하 캐리어 산란부터 이동도까지

전하 캐리어의 산란 현상은 반도체 내 전하 캐리어 이동도를 결정하는 데 중요한 역할을 한다. 전기장이 인가되면 전자와 정공은 그 힘에 의해 가속된다. 그러나 이러한 운동은 포논, 불순물, 결함 및 기타 캐리어들과의 상호작용에 의해 지속적으로 방해받는다. 이러한 산란 과정들은 캐리어가 무한히 가속되는 것을 막으며 그들의 평균 드리프트 속도를 결정하게 된다.

가장 단순한 수송 모델에서는 전하 캐리어 이동도와 산란 현상 간의 관계가 드루드 식에 의해 기술된다.

\[ \boxed{ \mu = \frac{e\tau_{\mathrm{tr}}}{m^*} } \]

여기서 \(e\)은 기본 전하량이며 \(m^*\)는 전하 캐리어의 유효 질량, \(\tau_{\mathrm{tr}}\)은 운동량 완화 시간입니다. 물리적으로 \(\tau_{\mathrm{tr}}\)는 외부 전기장에 의해 생성된 캐리어의 방향성 운동이 산란 현상에 의해 얼마나 빠르게 소실되는지를 나타냅니다. 완화 시간이 길수록 전하 캐리어 이동도는 증가하는 반면, 운동량 완화 현상이 강해질수록 이동도는 감소하게 됩니다.

산란률과 운동량 완화 과정의 관계

이전 절에서 계산된 총 산란률과 산란 현상에 의해 전하 캐리어의 운동량이 감소하는 속도 사이에는 중요한 차이가 존재한다. 총 외부 산란률인 \(1/\tau(\mathbf{k})\)은 전하 캐리어의 운동 방향이 얼마나 변하는지와 무관하게 초기 전자 상태에서 발생하는 모든 전이 현상을 모두 계산한다.

초기 파동벡터가 \(\mathbf{k}\)인 전자가 있다고 가정하자. 산란 사건 이후 그 파동벡터는 \(\mathbf{k}'\)로 변하며 원래 방향과 \(\theta_s\) 각도를 이룬다. 등방성 포물선형 밴드 구조에서 탄성 산란이 일어날 경우 전자의 속도는 변하지 않으나 최종 속도 중 원래 방향과 평행한 성분만이 남게 된다. 이 성분은

\[ v_{\parallel}' = v\cos\theta_s. \]

따라서 산란 과정에서 원래의 방향성 속도 중 소실되는 비율은 다음과 같다.

\[ \frac{v-v_{\parallel}'}{v} = \boxed{1-\cos\theta_s}. \]

이 간단한 기하학적 논증은 서로 다른 각도로 발생하는 산란 현상이 운동량 완화에 동일하게 기여하지 않는 이유를 설명한다. 예를 들어:

• 정면 산란(\(\theta_s=0^\circ\)): 가중치 값은 0이며, 이는 전자가 원래 진행 방향을 그대로 유지하기 때문이다.
• \(90^\circ\) 방향으로의 산란: 가중치는 1이며, 이는 전자가 원래의 운동 방향을 더 이상 유지하지 않기 때문입니다.
• 역산산란(\(\theta_s=180^\circ\)): 가중치는 2이며, 이는 전자가 운동 방향을 완전히 반전시키기 때문입니다.

이러한 각도 의존성을 활용하여 근사적인 운송 산란률을 도출할 수 있다. 등방성 밴드 구조의 경우 일반적으로 사용되는 각도 가중치 근사법은 다음과 같다.

\[ \begin{aligned} \frac{1}{\tau_{\mathrm{tr}}(\mathbf{k})} \simeq \frac{2\pi}{\hbar} \sum_{\mathbf{q}} (1-\cos\theta_s) \Big[ &|M_{\mathrm{abs}}(q)|^2 \delta(E'-E-\hbar\omega_{LO}) \\ +{}& |M_{\mathrm{em}}(q)|^2 \delta(E'-E+\hbar\omega_{LO}) \Big]. \end{aligned} \]

여기서 \(\theta_s\)은 초기 및 최종 전하 캐리어 속도 벡터 사이의 각도입니다. 등방성 포물선형 밴드 구조에서는 이러한 속도 벡터들이 각각 대응하는 파동벡터와 평행합니다. 행렬 요소 \(M_{\mathrm{abs}}\)와 \(M_{\mathrm{em}}\)는 이전 절에서 도출된 값들입니다.

비탄성 산란 현상에서는 캐리어의 속도 역시 변화하므로 앞서 언급된 단순 각도 가중치 모델은 근사치에 불과합니다. 등방성 밴드 구조를 갖는 경우 보다 일반적인 단일사건 투영 인자는 \(1-(v'/v)\cos\theta_s\)로 표현됩니다. 비탄성 수송 현상을 엄밀하게 다루기 위해서는 초기 및 최종 캐리어 상태 간의 상호작용을 고려한 선형화된 볼츠만 수송 방정식을 해석해야 합니다.

이러한 차이는 특히 프뢰리히 산란의 경우에 중요하게 작용한다. 해당 상호작용은 작은 파동벡터를 갖는 포논을 선호하므로 \(q\) 많은 산란 과정에서 캐리어의 진행 방향 변화가 비교적 미미하게 나타난다. 따라서 전체 산란률은 캐리어의 방향성 운동이 감쇠되는 속도와 상당히 다를 수 있다.

밴드 구조를 이용한 전하 캐리어 이동도 계산

드루드 모델 식은 유용한 출발점이 되지만, 실제 반도체에서는 단일의 유효질량이나 일정한 완화시간이 존재하지 않을 수 있습니다. 완화율과 전하 캐리어 속도 모두 전자 밴드 구조 전반에 걸쳐 변화할 수 있습니다.

이를 고려하기 위해 전하 캐리어 이동도는 수송에 관여하는 전자 상태들의 기여도를 평균하여 계산할 수 있다. 완화시간 근사법에 따르면 전자 이동도 텐서는

\[ \boxed{ \mu_{ij} = \frac{e}{n} \sum_b \int \frac{d^3k}{(2\pi)^3} v_{b,i}(\mathbf{k}) v_{b,j}(\mathbf{k}) \tau_{\mathrm{tr},b}(\mathbf{k}) \left( -\frac{\partial f}{\partial E} \right)_{E_b(\mathbf{k})} } \]

여기서 \(n\)은 전자 밀도이며 \(f\)는 페르미-디랙 분포입니다. 합산은 관련된 전자 대역들과 명시적으로 표현된 스핀 또는 밸리 상태들에 대해 이루어집니다. 각 상태의 군속도는 대역 분산 관계식으로부터 구해집니다.

\[ \mathbf{v}_b(\mathbf{k}) = \frac{1}{\hbar} \nabla_{\mathbf{k}}E_b(\mathbf{k}). \]

페르미-디랙 분포의 도함수는 전기적 응답에 기여하는 전자 상태들을 결정한다. 이후 적분 과정을 통해 해당 상태들의 속도와 운송 완화 시간을 결합함으로써 전체적인 전하 캐리어 이동도를 산출하게 된다.

이 방법은 특히 \(k\cdot p\) 이론을 이용해 계산된 유한차분 시간영역법(FDTD) 기반의 전자 밴드 구조에 유용한데, 이는 에너지에 따른 전자 속도를 일정한 유효질량을 가정하는 대신 계산된 분산 관계식으로부터 직접 도출할 수 있기 때문입니다. 또한 이 방법을 통해 다양한 산란 메커니즘들이 전하 캐리어 이동도에 미치는 영향을 각각의 전달완화시간 기여도를 통해 반영할 수 있습니다. 밴드 구조 기반의 이동도 계산 예시로는 Si/SiGe Δ-밸리 전자 이동도 튜토리얼 및 Si/SiGe 6-밴드 루팅거-콘 정공 이동도 튜토리얼이 있습니다. 해당 실리콘 관련 예제들에서는 프뢰리히 산란 대신 비극성 산란 메커니즘들이 사용되었으나, 전자 밴드 구조와 산란률로부터 전하 캐리어 이동도를 도출하는 과정은 동일하게 적용됩니다.