모트-거니 방정식의 유도

공간전하제한전류(SCLC)는 반도체 소자 물리학에서 가장 중요한 전하 수송 메커니즘 중 하나입니다. 이 현상은 전극에서 주입되는 전하 캐리어의 밀도가 반도체 내부에 이미 존재하는 평형 상태의 전하 캐리어 밀도보다 훨씬 커질 때 발생합니다. 이러한 조건 하에서는 전류가 물질 자체의 고유 전도도에 의해 제한되지 않고, 오히려 주입된 전하에 의해 생성된 전기장에 의해 제한됩니다.

공간전하제한전류(SCLC) 측정법은 유기 반도체, OLED, 페로브스카이트 소자 및 무질서한 물질에서 전하 캐리어 이동도를 추정하고 트랩 상태를 연구하는 데 널리 활용된다. 실험적으로는 양공 전용 또는 전자 전용 구조를 이용해 SCLC를 측정하는데, 이때 양쪽 접촉부가 동일한 종류의 캐리어만을 선택적으로 주입한다. 이러한 방식은 전자-양공 재결합을 억제함으로써 단일 캐리어 종류의 이동만이 지배적으로 일어나도록 한다.

트랩 상태가 없는 공간전하제한전류(SCLC) 현상을 기술하는 가장 널리 알려진 해석적 모델은 모트-거니 방정식으로, 이 방정식에 따르면 전류밀도는 인가된 전압의 제곱에 비례하여 변화한다.

$$ J = \frac{9}{8}\varepsilon\mu\frac{V^2}{L^3} $$

여기서 \(J\)은 전류밀도이며, \(\varepsilon\)는 재료의 유전율, \(\mu\)은 전하 캐리어 이동도, \(V\)는 인가된 전압, 그리고 \(L\)는 소자의 두께입니다.

비록 모트-거니 방정식 자체는 간결하지만 이를 도출하는 데 필요한 물리적 가정들이 매우 중요하다. 해당 유도 과정에서는 트랩 상태가 없는 전하 이동, 단일 종류의 전하 캐리어, 오옴성 접촉, 균일한 전하 캐리어 이동도, 미미한 드리프트-확산 현상 및 재결합이 거의 없다는 가정을 전제로 한다. 이러한 가정들을 정확히 이해하는 것은 실험적으로 측정된 공간전하제한전류 데이터를 해석하는 데 필수적이다.

공간전하제한전류의 물리적 기원

일반적인 오옴 전도 상태에서는 반도체 내부의 전하 캐리어 밀도가 대체로 일정하며 전류 밀도는 전압에 따라 선형적으로 변화합니다. 이 영역에서는 소자 내부의 전기장이 거의 균일하고 주입된 전하 캐리어 밀도 역시 비교적 낮은 수준입니다.

하지만 접촉부가 저전도성 반도체 내로 전하 캐리어를 효율적으로 주입할 경우, 주입된 전하가 소자 내부에 축적되기 시작한다. 이러한 축적된 전하는 포아송 방정식을 통해 내부 전기장의 분포를 변화시킨다. 그 결과 전기장은 위치에 따라 달라지게 되며, 주입된 전하 자체의 정전기적 상호작용에 의해 전류가 제한된다.

이러한 현상은 ??에서 설명되어 있다. 오옴 전도와는 달리 공간전하제한전류SCLC 영역에서는 내부 전기장이 균일하지 않다. 대신, 추가로 주입된 전하가 정전기적 전위를 지속적으로 변화시키기 때문에 전기장은 소자 내에서 거리에 따라 증가하게 된다.

기초 방정식들

Mott-Gurney 방정식의 유도 과정은 포아송 방정식과 드리프트 전류 방정식이라는 두 가지 기본 방정식에서 시작된다.

포아송 방정식은 전기장의 기울기와 전하 밀도 간의 관계를 나타냅니다:

$$ \frac{dE}{dx} = \frac{\rho}{\varepsilon} $$

여기서 \(E\)은 전기장이며 \(\rho\)는 전하 밀도이고 \(\varepsilon\)은 반도체의 유전율입니다.

밀도가 \(n\)인 단일 종류의 전하 캐리어를 가정할 때, 전하 밀도는 다음과 같이 표현된다:

$$ \rho = qn $$

여기서 \(q\)은 기본 전하량을 나타냅니다.

드리프트 전류 밀도는 다음과 같이 표현됩니다:

$$ J = q\mu nE $$

여기서 \(\mu\)은 전하 캐리어 이동도를 나타냅니다.

Mott-Gurney 유도 과정에서는 확산 전류가 무시된다. 전기장이 충분히 강하여 드리프트-확산 현상에서 드리프트 전류가 확산 전류보다 우세하게 되는 경우 이 근사가 타당하다.

완전한 드리프트-확산 방정식은 다음과 같습니다:

$$ J = q\mu nE + qD\frac{dn}{dx} $$

하지만 공간전하제한전류(SCLC) 조건 하에서는 확산 항이 드리프트 항에 비해 작다고 가정된다.

$$ qD\frac{dn}{dx} \ll q\mu nE $$

모트-거니 방정식의 유도

먼저 포아송 방정식과 드리프트 전류 방정식을 결합합니다. 그 후 치환을 수행합니다.

$$ n = \frac{\varepsilon}{q}\frac{dE}{dx} $$

드리프트 전류 식에 대입하면:

$$ J = q\mu nE $$

결과:

$$ J = q\mu \left( \frac{\varepsilon}{q}\frac{dE}{dx} \right) E $$

이는 다음과 같이 단순화된다:

$$ J = \varepsilon\mu E\frac{dE}{dx} $$

재배열하면:

$$ E\,dE = \frac{J}{\varepsilon\mu}\,dx $$

소자 전체에 걸쳐 적분하면:

$$ \int_0^{E} E\,dE = \frac{J}{\varepsilon\mu} \int_0^x dx $$

결과:

$$ \frac{E^2}{2} = \frac{Jx}{\varepsilon\mu} $$

따라서:

$$ E(x) = \sqrt{ \frac{2Jx}{\varepsilon\mu} } $$

소자 양단의 전압은 전기장을 적분함으로써 구할 수 있습니다:

$$ V = \int_0^L E(x)\,dx $$

전기장 표현식을 대입하면:

$$ V = \int_0^L \sqrt{ \frac{2Jx}{\varepsilon\mu} } \,dx $$

상수들을 적분 기호 밖으로 꺼내기:

$$ V = \sqrt{ \frac{2J}{\varepsilon\mu} } \int_0^L x^{1/2}\,dx $$

적분 수행 방법:

$$ V = \sqrt{ \frac{2J}{\varepsilon\mu} } \cdot \frac{2}{3}L^{3/2} $$

양변을 제곱하면:

$$ V^2 = \frac{8JL^3}{9\varepsilon\mu} $$

마지막으로 전류밀도에 대해 식을 재배열하면 모트-거니 방정식이 도출된다:

$$ J = \frac{9}{8} \varepsilon\mu \frac{V^2}{L^3} $$

유도 과정에 기반하는 물리적 가정들

수학적으로는 간단해 보이지만, 이 유도 과정은 여러 중요한 물리적 가정에 의존합니다.

트랩이 없는 전하 수송

이 유도 과정에서는 모든 주입된 전하 캐리어가 이동성을 유지하며 전기적 수송에 직접 기여한다고 가정한다. 그러나 실제 무질서 반도체에서는 전하 캐리어들이 종종 국소화된 트랩 상태에 갇히게 된다. 이러한 트랩 현상은 전하 캐리어 분포를 변화시키고 전류의 전압 의존성에도 영향을 미친다.

단일 캐리어 전송

모트-거니 방정식은 오직 한 종류의 전하 캐리어만이 전하 수송에 크게 기여한다고 가정한다. 전자와 정공이 모두 존재하는 양극성 소자의 경우, 재결합 및 전하 중성 조건의 영향으로 인해 단순한 공간전하제한전류(SCLC) 모델의 유도 과정이 타당하지 않게 된다.

오옴 접촉

이 유도 과정에서는 전하 캐리어들이 반도체 내부로 자유롭게 주입될 수 있다고 가정한다. 만약 접촉면에 쇼트키 장벽이나 주입 장벽이 존재할 경우, 전류는 소자 내부의 공간전하가 아닌 접촉 주입에 의해 제한된다.

균일한 전하 캐리어 이동도

이 유도 과정에서는 전하 캐리어 이동도가 전기장 및 캐리어 밀도와 무관하다고 가정한다. 하지만 다수의 유기 반도체에서는 이동도가 전기장, 물질의 무질서도 및 트랩 상태의 충만도에 크게 의존하는 경향이 있다.

확산 전류는 무시할 수 있음

본 유도 과정에서는 확산 전류는 무시된다. 이러한 근사법은 전기장이 약하거나 평형 상태에 가까울 때 그 정확도가 떨어진다.

Mott-Gurney 방정식이 적용되지 않는 경우

실제 반도체 소자들은 이상적인 모트-거니 거동에서 자주 벗어나는 경향이 있다. 가장 흔한 원인 중 하나는 트랩 상태의 존재이다. 트랩에 의해 제한되는 공간전하제한전류의 경우, 전류는 대개 거듭제곱 법칙을 따른다:

$$ J \propto V^m $$

여기서 \(m > 2\)입니다.

전기장에 의존하는 전하 캐리어 이동도 역시 전류-전압 관계에 큰 영향을 미칠 수 있다. 다수의 무질서한 유기 반도체에서는 전기장의 세기에 따라 이동도가 지수함수적으로 증가하는 풀-프렌켈 현상이 나타난다.

극도로 높은 전압 조건에서는 직렬 저항 및 자가발열 효과도 중요한 영향을 미칠 수 있다. 또한 전자와 정공이 동시에 존재하는 경우 재결합 현상으로 인해 전기장 분포가 변화하여 단일 캐리어 가정이 성립하지 않게 된다.

이러한 이유로 실제 반도체 소자의 물리적 거동을 정확히 모델링하기 위해 드리프트-확산 시뮬레이션이 자주 필요합니다.

드리프트-확산 시뮬레이션과의 관계

모트-거니 방정식은 공간전하제한전송 현상에 대한 해석적 통찰을 제공한다는 점에서 매우 유용하다. 그러나 실제 반도체 소자들에는 트랩 분포, 전기장에 의존하는 전하 캐리어 이동도, 비오옴성 접촉면, 재결합 과정 및 공간적으로 변화하는 전기장 등이 존재하며, 이러한 요소들은 단순한 해석적 모델로는 포착될 수 없다.

OghmaNano와 같은 드리프트-확산 시뮬레이션 도구들은 연립된 포아송 방정식 및 연속성 방정식을 수치적으로 해석함으로써 트랩 상태, 이동도 불균일성, 복잡한 다층 구조 등과 같은 비이상적인 효과들을 명시적으로 반영할 수 있게 한다. 이를 통해 무질서한 구조, 트랩 상태, 이동도 구배가 존재하는 실제 공간전하제한전류 소자들을 연구하는 것이 가능해진다.