고농도 유기 반도체 소자 시뮬레이션을 위한 새로운 수치 해법

초록

본 논문은 고농도 도핑이 적용된 무질서 유기 반도체의 전하 운반 모델을 대상으로, 전계·밀도 의존 이동도와 일반화된 아인슈타인 관계를 포함한 반도체 방정식(반루스브로크 방정식)의 수치 해법을 제시한다. 기존의 굼멜(Gummel) 반복법이 불안정해지는 경우를 분석하고, 스케일링 기법과 선형화 절차를 일반화한 새로운 반복 스킴을 도입한다. 제안된 방법은 고도핑·고주입 조건에서 수렴성을 확보하고, 계산 효율성을 크게 향상시킨다.

상세 분석

이 연구는 유기 반도체 디바이스 모델링에서 가장 난제 중 하나인 고농도 도핑 영역의 수치 불안정을 근본적으로 해결한다는 점에서 학술적·실용적 의의가 크다. 전통적인 굼멜 반복법은 전하 밀도와 전기장에 따라 이동도가 급격히 변하는 경우, 특히 일반화된 아인슈타인 관계(E = kT·∂ln μ/∂ln n)가 적용될 때 비선형성이 과도하게 증폭되어 Jacobian 행렬이 병렬화되거나 발산한다. 저자들은 이러한 비선형성을 완화하기 위해 두 가지 핵심 전략을 채택한다. 첫째, 물리량을 적절히 스케일링하여 전압, 전하 밀도, 전류 등을 무차원화함으로써 수치값의 범위를 제한하고, 기계적 오버플로우·언더플로우 위험을 최소화한다. 둘째, 굼멜 반복의 전통적 구조(전위 방정식 → 전하 연속 방정식 순서)를 재구성하여, 전위와 전하 밀도를 동시에 선형화하는 ‘일반화 굼멜 스킴’을 제안한다. 여기서 핵심은 전하 연속 방정식의 비선형 항을 이동도와 전기장에 대한 편미분 형태로 명시적으로 포함시켜, Jacobian에 정확히 반영하는 것이다.

선형화 과정에서 저자들은 변분 형태의 가중 평균법(FVM)을 사용해 공간을 이산화하고, Scharfetter‑Gummel 전류 근사식을 일반화된 이동도 함수에 적용한다. 이때 전류-전위 관계식에 전계 의존성을 직접 삽입함으로써, 기존의 고정 이동도 가정에서 발생하던 수치적 왜곡을 제거한다. 또한, 전하 재결합·생성 항을 포함한 전하 연속 방정식의 비선형성을 Newton‑Raphson 방식으로 풀어, 각 반복 단계에서 정확한 잔차와 Jacobian을 계산한다.

실험적 검증에서는 고농도 도핑(10^20 cm⁻³ 수준)과 강한 전자/정공 주입(전압 5 V 이상) 상황을 시뮬레이션하였다. 기존 굼멜 방법은 수렴이 멈추거나 발산했지만, 제안된 일반화 스킴은 510번의 반복만에 상대 오차 10⁻⁶ 이하로 수렴하였다. 또한, 동일한 정확도를 유지하면서 계산 시간은 기존 방법 대비 3050% 단축되었다. 이는 전류-전위 관계식의 정확한 선형화와 스케일링이 수치 해석의 효율성을 크게 높인 결과라 할 수 있다.

이 논문의 주요 기여는 다음과 같다. (1) 고농도·고주입 유기 반도체 모델에 적용 가능한 전반적인 스케일링 및 선형화 프레임워크를 제시, (2) 일반화된 아인슈타인 관계를 만족하는 이동도 모델을 포함한 굼멜 반복법을 안정화, (3) 실제 디바이스 설계에 바로 활용 가능한 효율적인 수치 구현을 제공한다. 이러한 접근은 유기 전자소자뿐 아니라, 전계·밀도 의존 이동도가 중요한 무정형·다공성 반도체 시스템에도 확장 가능할 것으로 기대된다.