반도체의 드리프트-확산 현상에서의 아인슈타인 관계

1. 오제 관계란 무엇인가?

반도체 내부에서는 에너지 밴드의 기울기로 인한 드리프트 현상과 캐리어 밀도의 기울기로 인한 확산 현상이라는 두 가지 메커니즘을 통해 전하 캐리어가 이동한다. 아인슈타인 관계식은 각각의 전하 캐리어 이동도와 확산계수를 서로 연관시켜 열적 평형 상태에서 두 메커니즘이 정확히 균형을 이루도록 한다. 전자와 정공에 대한 드리프트-확산 방정식은 다음과 같다.

\[ \begin{aligned} \boldsymbol{J}_n &= \underbrace{\color{red}{\mu_n n\nabla E_c}}_{\text{drift}} + \underbrace{\color{blue}{qD_n\nabla n}}_{\text{diffusion}},\\[5pt] \boldsymbol{J}_p &= \underbrace{\color{red}{\mu_p p\nabla E_v}}_{\text{drift}} - \underbrace{\color{blue}{qD_p\nabla p}}_{\text{diffusion}}. \end{aligned} \]

여기서 \(\boldsymbol{J}_n\)과 \(\boldsymbol{J}_p\)는 각각 전자와 정공의 전류 밀도이며, \(q\)은 양의 기본 전하량이다. \(n\)와 \(p\)는 자유 캐리어 밀도, \(\mu_n\)와 \(\mu_p\)은 전하 캐리어 이동도, \(D_n\)과 \(D_p\)는 확산 계수에 해당한다. \(E_c\)과 \(E_v\)은 각각 전도대 및 가전대 가장자리의 에너지값으로 단위는 줄이며, 이들의 기울기는 전기장 항을 대체한다. 만약 대역 가장자리 값이 볼트 단위로 저장된다면(에너지를 \(q\)로 나눈 값), 전자에 대한 드리프트 항은 \(q\mu_n n\nabla (E_c/q)\)이 되며 정공의 경우도 마찬가지이다. 본 식들은 등온 조건 하에서 대역 구조가 국소적으로 일정하다고 가정하므로 온도 구배나 공간적으로 변하는 물질 특성은 별도로 고려해야 한다. 배경 지식은 전기정역학, 전하 캐리어 이동도, 그리고 드리프트-확산 이론의 전체 유도 과정을 참조하기 바란다.

각 방정식에는 두 가지 항이 존재합니다. 빨간색으로 표시된 드리프트 항은 에너지 밴드 경사에 의해 발생하는 전하 캐리어의 이동을 나타내며, 파란색 확산 항은 캐리어 밀도 차이에 의해 발생하는 이동을 설명합니다. 그림 ??에서 전자들은 전도대 에너지 경사를 따라 이동하여 오른쪽에 축적됩니다. 확산 현상에 의해 전자들은 밀도가 낮은 영역으로 다시 이동하게 됩니다. 전자들의 이동 방향이 일반적인 전류 방향과 반대이지만, 각각의 전류 기여분은 서로 반대 방향이어서 평형 상태에서는 상쇄됩니다.

열적 평형 상태에 있는 조명이 없는 소자, 즉 내부 전기장이 존재하는 p-n 접합을 포함한 소자에서는 전자 및 정공에 대한 순전류가 각각 국소적으로 0이 되어야 한다:

\[ \boldsymbol{J}_{n,\mathrm{drift}}+\boldsymbol{J}_{n,\mathrm{diff}}=\boldsymbol{0}, \qquad \boldsymbol{J}_{p,\mathrm{drift}}+\boldsymbol{J}_{p,\mathrm{diff}}=\boldsymbol{0}. \]

만약 아인슈타인 관계식이 캐리어 통계와 일치하지 않는다면 드리프트와 확산 현상이 서로 상쇄되지 않을 수 있다. 이 경우 시뮬레이션 결과로서 외부 전압이 0인 상태에서도 허위의 암전류가 예측되거나, 그 암전류-전압 곡선에 오프셋이 발생할 수 있으며 이는 그림 ??에 예시되어 있다. 이러한 현상은 특히 매우 작은 전류를 시뮬레이션할 때 심각한 문제가 된다. 장치에 광조사가 가해지거나 기타 방식으로 구동되는 경우 외부 바이어스가 0이라 하더라도 평형 상태가 성립하지 않으므로, 여기서 언급된 ‘전류가 0인 상태’는 열적으로 평형을 이룬 암전 상태의 장치를 의미한다.

일반적인 비축퇴 상태(볼츠만 통계 적용 시)에서는 아인슈타인 관계식은 다음과 같이 간단히 표현됩니다.

\[ \boxed{D_n=\mu_n\frac{k_BT}{q},\qquad D_p=\mu_p\frac{k_BT}{q}.} \]

여기서 \(k_B\)은 볼츠만 상수이며 \(T\)는 절대온도입니다. \(k_BT/q\)의 값은 열적전압으로, 300 K에서는 약 25.9 mV에 해당합니다. 따라서 전하 캐리어 이동도가 높을수록 확산 계수 역시 비례적으로 커져야 하며, 일관된 소자 모델에서는 이 둘을 독립적으로 설정할 수 없습니다. 이후 절에서는 이 식을 유도하고 축퇴 상태나 비포물선형 밴드 구조에서 어떤 변화가 발생하는지 설명합니다.

2. 간단한 유도 과정: 고전적인 아인슈타인 관계식

균일하고 등온적이며 비축퇴 상태인 반도체 내의 전자들을 고려해보자. 열적 평형 상태에서는 전자의 준페르미 준위 \(E_F\)가 공간적으로 일정하며, 캐리어 밀도는 볼츠만 분포를 따른다. 볼츠만 분포와 페르미-디랙 분포 간의 차이점에 대해서는 페르미-디랙 통계와 맥스웰-볼츠만 통계에서 설명하고 있다.

\[ n=N_c\exp\!\left(\frac{E_F-E_c}{k_BT}\right). \]

여기서 \(E_c\)은 전도대 가장자리를 나타내며 \(N_c\)는 유효 전도대 상태밀도입니다. \(N_c\)와 \(T\)가 일정할 때 공간 기울기를 취하면 다음과 같은 결과가 얻어집니다.

\[ \nabla n=-\frac{n}{k_BT}\nabla E_c. \]

이 상태밀도 기울기를 전도대 기울기 표기법을 활용해 전자 전류식에 직접 대입하면 다음과 같은 결과가 도출된다.

\[ \boldsymbol{J}_n =\mu_n n\nabla E_c+qD_n\nabla n =\mu_n n\nabla E_c-qD_n\frac{n}{k_BT}\nabla E_c =n\nabla E_c\left(\mu_n-\frac{qD_n}{k_BT}\right). \]

내부 전기장이 존재하는 경우라도 평형 상태에서는 전류가 0이어야 하므로 대괄호 안의 식은 반드시 0이 되어야 한다. 따라서

\[ \boxed{D_n=\mu_n\frac{k_BT}{q}.} \]

정공에 대해서도 동일한 논리를 적용하면 \(D_p=\mu_p k_BT/q\)과 같은 결과가 도출된다. 이는 널리 알려진 고전적인 아인슈타인 관계식이다. 이 관계식은 일반적으로 포물선형 에너지 밴드를 가정하여 유도되지만, 여기서 핵심적인 근사 조건은 볼츠만 통계(비축퇴 조건)이다. 즉 캐리어들의 밀도가 충분히 낮아서 가능한 상태들에 대한 점유 확률이 미미한 수준인 경우에만 성립한다.

3. 일반화된 아인슈타인 관계식

단순한 수식 결과가 항상 적용되는 것은 아니다. 고농도 도핑된 물질이거나 페르미 준위가 에너지 대역에 근접하거나 그 안으로 진입하는 경우, 전하 캐리어들은 볼츠만 근사 대신 페르미-디랙 통계를 따르게 된다. 이때 캐리어 밀도는 페르미 준위의 단순한 지수 함수 형태로 나타나지 않는다. 또한 상태밀도 역시 비포물선형이거나 무질서한 형태를 띨 수 있으며, 밴드 구조 계산을 통해 수치적으로 도출될 수도 있다. 캐리어 밀도 계산에 관해서는 Charge carrier density에서 상세히 설명한다.

전자 밀도가 일정한 온도 조건 하에서 그리고 국부적인 상태밀도가 변하지 않는다고 가정할 때, 에너지 차이 \(\xi_n=E_{Fn}-E_c\)에 따라 결정된다고 하자: \(n=n(\xi_n,T)\). 평형 상태에서는 전자의 준페르미 준위 \(E_{Fn}\)가 일정한 값을 유지한다. 따라서 이 준위에 대한 공간적 기울기를 계산하면 다음과 같은 결과가 얻어진다.

\[ \nabla n =\left.\frac{\partial n}{\partial\xi_n}\right|_T (\nabla E_{Fn}-\nabla E_c) =-\left.\frac{\partial n}{\partial\xi_n}\right|_T\nabla E_c. \]

이 결과를 전류 방정식에 대입하고 평형 상태에서 전류가 0이 되도록 조건을 설정하면 일반화된 아인슈타인 관계식이 도출된다.

\[ \boxed{ \frac{D_n}{\mu_n} =\frac{n}{q\left(\partial n/\partial\xi_n\right)_T} =\frac{1}{q}\left(\frac{\partial\xi_n}{\partial\ln n}\right)_T. } \]

정공의 경우, 정공 밀도가 증가함에 따라 값이 커지는 \(\xi_p=E_v-E_{Fp}\)을 정의한다. 그러면 다음과 같은 등가 관계식이 성립한다.

\[ \boxed{ \frac{D_p}{\mu_p} =\frac{p}{q\left(\partial p/\partial\xi_p\right)_T}. } \]

이 방정식들은 물리적으로 유의미한 의미를 지닙니다. 즉, 화학적 포텐셜이 변할 때 캐리어 밀도가 얼마나 쉽게 변하는지를 이 도함수가 나타내는 것입니다. 준페르미 에너지가 아주 조금만 변해도 캐리어 밀도가 크게 변한다면 이 도함수의 값은 커지게 되며 \(D/\mu\)은 상대적으로 작아집니다. 반면, 상태들이 모두 채워져서 캐리어 밀도의 변화가 적어진다면 이 비율값은 증가하게 됩니다. 따라서 아인슈타인 관계식은 단순히 이동도에 곱해지는 수치적 계수가 아니라 캐리어 통계와 상태밀도의 속성을 반영하는 것입니다.

임의의 전도대 상태밀도 \(g_c(E)\)에 대해 전자 밀도는 다음과 같이 계산할 수 있다.

\[ n=\int_{E_c}^{\infty}g_c(E) \frac{1}{1+\exp[(E-E_{Fn})/(k_BT)]}\,dE. \]

일정한 밴드 구조와 온도 조건 하에서는, 그 도함수를 동일한 상태들로부터 산출할 수 있다.

\[ \left.\frac{\partial n}{\partial E_{Fn}}\right|_T =\frac{1}{k_BT}\int_{E_c}^{\infty}g_c(E)f(E)[1-f(E)]\,dE, \qquad f(E)=\frac{1}{1+\exp[(E-E_{Fn})/(k_BT)]}. \]

고정된 \(E_c\) 조건 하에서 \(\partial\xi_n/\partial E_{Fn}=1\)이 성립하므로 이 도함수를 일반화된 아인슈타인 관계식에 직접 적용할 수 있다. 이 방법은 비포물선형 수치 상태밀도의 경우에도 유효하다. 정공의 경우에는 해당 가전대 상태밀도와 정공 점유도 \(1-f\)를 사용하여 \(E_v-E_{Fp}\)에 대한 도함수를 계산하면 된다.

4. 포물선형 밴드: 고전적 결과의 재도출

3차원 파라볼릭 밴드의 경우, 전자 밀도는 페르미-디랙 적분을 이용해 표현될 수 있습니다:

\[ n=N_c F_{1/2}(\eta_n), \qquad \eta_n=\frac{E_{Fn}-E_c}{k_BT}. \]

우리는 정규화된 정의 \(F_j(\eta)=\Gamma(j+1)^{-1} \int_0^\infty x^j/[1+\exp(x-\eta)]\,dx\)를 사용합니다. \(dF_{1/2}/d\eta=F_{-1/2}\)를 활용하면 일반화된 아인슈타인 관계식은 다음과 같이 표현됩니다.

\[ \boxed{ \frac{D_n}{\mu_n} =\frac{k_BT}{q} \frac{F_{1/2}(\eta_n)}{F_{-1/2}(\eta_n)}. } \]

정공의 경우에도 \(\eta_p=(E_v-E_{Fp})/(k_BT)\)을 이용해 유사한 식이 성립한다. 이 결과는 3차원 파라볼릭 밴드에서 밴드가 축퇴 상태일 때도 적용된다. 두 페르미-디랙 적분값의 비율은 상태 충만도의 영향을 반영한다.

비축퇴 조건 하에서는 \(\eta_n\ll-1\)에서 두 페르미-디랙 적분값이 모두 \(\exp(\eta_n)\)에 수렴하게 된다. 이에 따라

\[ \frac{F_{1/2}(\eta_n)}{F_{-1/2}(\eta_n)}\longrightarrow 1, \qquad \boxed{\frac{D_n}{\mu_n}\longrightarrow\frac{k_BT}{q}.} \]

이를 통해 우리가 익히 알고 있는 수식이 어디서 유래하는지 명확히 알 수 있다. 단순히 포물선형 밴드만으로는 충분하지 않으며, \(D/\mu=k_BT/q\)을 얻기 위해서는 캐리어가 비축퇴 상태여야 한다. 축퇴된 포물선형 밴드의 경우 페르미-디랙 보정항이 필요하지만, 일반적인 비포물선형 밴드에서는 밀도 미분값을 활용한 표현식이 보다 유용한 출발점이 된다.

5. 캐리어 압력: 확산 현상을 이해하는 또 다른 방법

확산 현상이 발생하는 원인을 이해하는 데 또 하나의 유용한 방법이 존재한다. 볼츠만 수송 방정식으로부터 도출된 운동량 균형 방정식에 따르면, 전하 캐리어들은 각기 다른 속도와 운동량을 가지므로 압력을 발생시킨다. 이러한 캐리어 압력의 기울기는 캐리어들이 압력이 높은 영역에서 낮은 영역으로 이동하도록 하는 힘을 만들어내며, 이것이 바로 확산 현상의 미시적 기원이다. 운동량 균형과의 연관성에 대해서는 볼츠만 수송 방정식 유도를 참조하기 바란다.

등방성인 3차원 포물선형 전도대의 경우, 전자 압력을 \(P_n\)과 같이 정의한다.

\[ P_n=\frac{2}{3}\int_{E_c}^{\infty}(E-E_c)g_c(E)f(E)\,dE. \]

여기서 \(g_c(E)\)은 전도대의 상태밀도이며 \(f(E)\)는 입자 점유 확률입니다. \(2/3\) 인수는 3차원에서의 포물선형 분산 관계로부터 도출됩니다. 비축퇴 상태의 고전적 조건 하에서는 이는 잘 알려진 이상기체 방정식으로 전환됩니다.

\[ P_n=nk_BT. \]

균일하고 등온인 물질의 경우, 전자 전류는 밀도 기울기 대신 압력 기울기를 이용해 표현할 수 있다.

\[ \boldsymbol{J}_n=\mu_n n\nabla E_c+\mu_n\nabla P_n =\mu_n n\nabla E_c+qD_n\nabla n. \]

두 확산 항을 비교해보면, 중요한 것은 캐리어 밀도에 따라 압력이 얼마나 빠르게 변하는가이며 단순히 압력을 밀도로 나눈 값이 아닙니다.

\[ \boxed{\frac{D_n}{\mu_n}=\frac{1}{q}\left(\frac{\partial P_n}{\partial n}\right)_T.} \]

따라서 \(P_n=nk_BT\) 볼츠만 한계 하에서 \(\partial P_n/\partial n=k_BT\)그러면 고전적인 아인슈타인 관계식이 즉시 도출된다. 전자 기체가 축퇴 상태가 되면 압력은 더 이상 밀도에 대해 선형적으로 증가하지 않는다. 3차원의 포물선형 밴드 구조에서 위에서 사용한 동일한 정규화된 페르미-디랙 적분식을 적용하면

\[ P_n=N_c k_BT F_{3/2}(\eta_n),\qquad \frac{\partial P_n}{\partial n} =k_BT\frac{F_{1/2}(\eta_n)}{F_{-1/2}(\eta_n)}. \]

이를 통해 4절에서 도출된 일반화된 아인슈타인 관계식이 정확히 재현된다. \(P_n/n=k_BT F_{3/2}(\eta_n)/F_{1/2}(\eta_n)\)과 \(\partial P_n/\partial n\)를 혼동해서는 안 되는데, 이 둘은 오직 고전적 한계 조건에서만 서로 같기 때문이다. 비포물선형 밴드의 경우 캐리어 압력을 계산할 때는 속도 및 밴드 구조에 의존하는 보다 일반적인 처리가 필요하다. 따라서 일관된 드리프트-확산 모델 내에서 아인슈타인 관계식을 산출하는 데에는 3절에서 제시된 캐리어 밀도에 대한 미분값을 활용하는 방법이 보다 편리하다.

6. 수치적 정확도: 허위 평형 전류 발생 방지

드리프트-확산 시뮬레이터에서는 일반화된 아인슈타인 관계식을 작성하는 것은 간단하지만 이를 정확하게 계산하는 것은 어려울 수 있다. 주된 난제는 준페르미 에너지에 대한 캐리어 밀도의 도함수를 구하는 것이다. 일반적으로 사용되는 근사법으로는 유한차분법이 있다:

\[ \frac{\partial n}{\partial E_{Fn}} \approx\frac{n(E_{Fn}+\Delta E)-n(E_{Fn}-\Delta E)}{2\Delta E}. \]

만약 \(\Delta E\) 만약 이 값이 지나치게 크면 절단 오차로 인해 도함수의 값이 부정확해진다. 반대로 지나치게 작으면 거의 같은 값을 갖는 두 캐리어 밀도가 서로 빼어지게 되어 부동소수점 반올림 오차나 도파관 구조 및 상태밀도 계산에서 발생하는 수치적 노이즈가 증폭될 수 있다. 이로 인해 발생한 오차는 이후 계산 과정에서 더욱 증폭된다. \(D/\mu=n/[q(\partial n/\partial E_{Fn})]\)특히 미분값이 작은 경우에 그러하다.

아주 작은 상대적 오차라도 문제가 될 수 있는데, 이는 평형 상태에서는 드리프트 전류와 확산 전류가 각각 상당한 크기를 가지지만 그 합은 정확히 0이어야 하기 때문입니다. 이러한 미세한 불일치로 인해 허위의 잔류 전류가 발생할 수 있습니다. 이러한 현상은 비포물선형 밴드 구조, 상태밀도의 수치적 계산, 낮은 캐리어 밀도 또는 밴드 구조가 급격히 변하는 이종접합 구조를 다루는 시뮬레이션에서 특히 두드러지게 나타납니다.

가능하다면 3절에서 설명한 바와 같이 상태밀도 적분 하에서 점유 함수를 해석적으로 미분하는 것이, 각각 따로 계산된 캐리어 밀도 값들에 유한차분법을 적용하는 것보다 바람직합니다. 캐리어 밀도가 매우 낮은 경우에는 적절할 때 볼츠만 근사식을 활용하는 것이 좋습니다. 또한 수치적 적분 방식, 캐리어 통계 모델, 밴드 가장자리 정의 방식 및 드리프트-확산 과정의 이산화 처리 방식이 서로 일관성을 유지하도록 하는 것도 매우 중요합니다.

마지막으로, 올바른 아인슈타인 관계식만으로는 부적절한 전류를 제거하기에 충분하지 않으며, 공간 이산화 과정 역시 평형 상태를 유지해야 한다. 특히, 플럭스 이산화 방식은 평형 상태에서의 일정한 페르미 준위를 준수해야 하는데, 전형적인 드리프트-확산 시뮬레이션에서는 Scharfetter–Gummel과 같은 평형 보존형 기법이 흔히 활용된다. 실제 시뮬레이터 검증 시에는 암전 상태에서 바이어스가 0인 평형 계산을 수행하여 소자 전체에 걸쳐, 특히 물질 계면을 포함하여 전자 및 정공 전류가 요구되는 수치적 허용오차 범위 내에서 소멸하는지 확인해야 한다. 접합부 및 밴드 오프셋 처리에 관해서는 Interfaces 및 The built-in potential을 참조하기 바란다.

핵심 내용은 간단합니다. 바로 아인슈타인 관계식이 드리프트-확산 현상에 대한 기술적 설명들을 열역학적으로 일관되게 만든다는 점입니다. 널리 알려진 고전적 형태의 아인슈타인 관계식은 사용이 용이하지만, 보다 일반적인 캐리어 통계 모델을 적용하려면 그에 상응하는 일반적이고 수치적으로 정확한 관계식이 필요합니다. 보다 포괄적인 운송 현상 프레임워크는 볼츠만 방정식에서 드리프트-확산 모델으로의 유도에서 상술되어 있습니다.