홈 예제 스크린샷 사용자 매뉴얼 Bluesky 로고 YouTube ☰
OghmaNano 유기/페로브스카이트 태양전지, OFET 및 OLED 시뮬레이션 다운로드 Quick Start guide

산란, 완화 시간 및 전하 캐리어 이동도

학습 목표. 이 페이지를 읽은 후 다음을 수행할 수 있어야 합니다:

사전 요구사항: 파동함수 페이지와 Bir–Pikus 페이지가 필요합니다.

예상 소요 시간: 30분.

양자 우물 내에서 전하 캐리어 이동도는 산란 현상에 의해 제한된다. OghmaNano는 현재 변형 퍼텐셜 모델을 통해 양성자 산란을 구현하며, 이때 페르미의 황금법칙을 활용하여 산란률을 계산한다. 이후 얻어진 완화 시간 값을 드리프트-확산 이론의 완화시간 근사법에 따라 전하 캐리어 이동도로 변환한다. 본 문서는 해당 산란 기구에 대해 설명하며, 구현된 기능과 그렇지 않은 기능을 명확히 기술한다.

1. 페르미의 황금법칙

초기 상태 \(|i\rangle\)에서 섭동 \(\delta H\) 하에 최종 상태 \(|f\rangle\)로의 전이율은 다음과 같다.

\[ \Gamma_i = \frac{2\pi}{\hbar}\sum_f \big|\langle f|\delta H|i\rangle\big|^2\, \delta\!\big(E_f - E_i\big), \]

에너지 보존을 만족시키는 탄성 산란 과정을 기술하는 디랙 델타 함수이다. 소자 온도 하에서 발생하는 음향 포논 산란은 거의 탄성적이며, 등분배 법칙이 성립하는 조건에서는 포논의 점유 수가 \(k_B T\)이라는 인수로 작용한다. 여기서 섭동항은 음향 포논의 변형 전위장이 에너지 밴드에 미치는 영향을 나타낸다.

2. 음향 변형 전위

음향 포논은 격자를 국부적으로 변형시키며, 이러한 변형은 변형 퍼텐셜을 통해 밴드 가장자리를 이동시킨다 \(D\). OghmaNano는 파수벡터가 \(\mathbf{q}\)이고 편광이 \(\hat{\mathbf{e}}\)인 포논에 의한 변형장을 계산한 뒤 이를 밴드와 결합시키는데, 특히 정적 변형에 사용되는 동일한 Bir–Pikus 변형 퍼텐셜을 통해 정공 밴드와도 결합한다. 등분배 원리가 적용되는 조건 하에서 포논 분포에 대한 제곱 행렬요소의 평균값을 구하면 앞서 언급된 계수가 산출된다.

\[ \text{prefactor} = \frac{D^2\,k_B T}{\rho\,v_s^2\,\hbar}, \]

여기서 \(\rho\)은 질량 밀도이며 \(v_s\)는 음속입니다. 산란률은 온도에 따라 선형적으로 증가하는데, 이는 음향 변형전위 산란의 특징적인 현상입니다.

3. 최종 상태밀도와 형상 인자

두 가지 추가 요소가 산란률 계산을 완성한다. 각 최종 서브밴드에 대해 파동벡터를 찾음으로써 에너지 델타 함수가 평가된다. \(k_f\) 그 시점에서 \(E_f(k_f)=E_i\)그리고 그 위치의 2차원 상태밀도로 가중치를 부여하는데, 이 값은 다음에 비례한다 \(k_f/|dE/dk|\)그리고 초기와 최종 파동함수의 겹침도 — 형태 인자 \(G\) 이는 두 상태가 양자 우물 내에서 얼마나 동일한 영역을 공유하는지를 나타내며, 상태들 간의 공간적 중첩이 적을 경우 서브밴드 간 산란이 억제된다. 특정 상태에 대한 산란률은 다음과 같다. \(i\) 따라서

\[ \Gamma_i = \text{prefactor}\,\sum_f \frac{k_f\,G_{if}}{|dE/dk|_f}, \qquad \tau_i = \frac{1}{\Gamma_i}. \]

완화 시간 \(\tau_i\)은 특정 상태 \(i\)에서 벗어나는 총 산란률의 역수에 해당합니다.

4. 완화 시에서 전하 캐리어 이동도로

전하 캐리어 이동도는 완화 시간 근사법 하에서 선형화된 볼츠만 방정식으로부터 도출된다. 각 상태는 그룹 속도의 제곱, 완화 시간 및 상태 점유도의 미분값에 비례하여 전도도에 기여한다. 여기서 그룹 속도는

\[ v(k) = \frac{1}{\hbar}\,\frac{dE}{dk}, \]

또한 페르미 준위 근처의 상태들만이 전하 수송에 기여하도록 제한하는 점유도 미분값이 \(-\partial f/\partial E = f(1-f)/k_B T\)입니다. 서브밴드들에 대해 적분을 수행하면

\[ \sigma \propto q^2 \sum_n \int \frac{dk}{(2\pi)}\; k\,v(k)^2\,\tau(k)\,\big[-\partial f/\partial E\big], \]

전하 캐리어 이동도는 단위 전하량 및 2차원 전하 밀도당 전도도를 의미한다.

\[ \mu = \frac{\sigma}{q\,p_{2D}}, \]

동일한 분산 관계에 기반하여 계산된 면밀도 \(p_{2D}=\sum_n\int (k\,dk/2\pi)\,f\)를 이용하면, 단일 질량값으로 추정한 값이 아닌 실제 비포물선적 서브밴드 구조를 반영하는 유효 전하 캐리어 이동도가 산출된다.

5. 가정사항 및 한계점

구현된 모델은 의도적으로 특정 조건에만 적용되도록 설계되었습니다. 이 모델은 탄성 영역 및 등분배 조건 하에서 변형 퍼텐셜을 이용해 음향 포논 산란을 다루며, 완화 시간 근사법을 통해 전하 캐리어 이동도를 산출합니다. 극성 광학(Fröhlich) 산란, 이온화 불순물 산란, 계면 거칠기에 의한 산란, 합금 불규칙성에 기인한 산란, 그리고 등분배 조건을 넘어서는 비탄성 포논 방출/흡수 현상은 이 모델에 포함되지 않습니다. 이러한 산란 기구가 지배적인 경우—즉, 상온에서의 극성 물질이나 불순물 농도가 높거나 계면이 거친 양자 우물의 경우—계산된 이동도 값은 전체적인 전하 수송 예측값이 아니라 음향 포논 산란에 의해 제한된 이동도 성분으로 해석되어야 합니다.

구현 관련 참고사항: 산란 계산 루틴은 정적 변형 퍼텐셜 해밀토니안에서 사용되는 것과 동일한 변형 퍼텐셜을 통해 음향 포논에 의한 변형을 에너지 밴드에 반영합니다. 이후 상태밀도 및 형상 인자들을 활용하여 산란율을 산출하고 각 상태별로 완화 시간 \(\tau=1/\Gamma\)을 저장합니다. 전하 캐리어 이동도 계산 시에는 군속도를 구하기 위해 수치적으로 미분된 분산 관계식을 이용합니다. 현재 모델에서는 음향 변형-퍼텐셜 산란만을 구현하고 있을 뿐이며, 모든 포논 및 캐리어 산란 기구를 포괄한다고 보장하지 않습니다.

주요 내용