페로브스카이트 태양전지에서의 CELIV 측정: Part a
1. 개요
이 튜토리얼에서는 먼저 선형 증가 전압에 의한 전하 추출 (CELIV) 기법의 기본 개념을 소개하고, 이 기법이 박막 반도체에서 캐리어 이동도와 밀도를 조사하는 데 어떻게 사용되는지 설명합니다. 이론적 배경을 정리한 뒤, 실제 단계로 빠르게 넘어가 OghmaNano 내에서 CELIV 시뮬레이션을 설정하고 실행하는 방법을 보여줍니다. 이론과 실습 시뮬레이션의 이러한 조합은 방법에 대한 명확한 이해와 이를 자신의 소자에 적용할 수 있는 도구를 모두 제공합니다.
2. CELIV 배경
선형 증가 전압에 의한 전하 추출(CELIV) 기법은 유기 및 페로브스카이트 태양전지에서 전하 캐리어 이동도와 밀도를 연구하는 데 사용되는 실험 방법입니다. 일반적인 구성에서는 태양전지가 선형으로 증가하는 역바이어스를 인가하는 신호 발생기에 연결되고, 그에 따른 전류 과도 응답은 오실로스코프에 기록됩니다. 선택적으로 연속 조명 또는 짧은 레이저 펄스를 추가하여 전하 캐리어를 생성할 수 있습니다. 이 구성의 개략도는 ??에 나와 있습니다.
전압 램프가 인가되면 측정 전류는 두 가지 뚜렷한 기여로 구성됩니다: 소자의 기하학적 정전용량에서 오는 직사각형 기준선과, 반도체 내 이동 가능한 전하 캐리어 추출에 의한 추가 피크입니다. 기하학적 정전용량 항은 항상 존재하는 반면, 캐리어 관련 피크는 이동도와 캐리어 밀도에 대한 직접적인 정보를 제공합니다. 이러한 전압 및 전류 과도 응답의 예는 ??에 나와 있습니다.
캐리어 이동도는 추출 전류가 최대값에 도달하는 시간을 측정함으로써 CELIV 과도 응답으로부터 직접 결정할 수 있습니다. 균일한 캐리어 분포와 무시할 수 있는 재결합이라는 가장 단순한 경우에서, 이동도는 CELIV 방정식으로 주어집니다: \[ \mu = \frac{2 d^{2}}{3 A t_{\text{max}}^{2}} \] 여기서 \(d\)는 active 레이어 두께, \(A\)는 전압 램프 속도 (\(A = \mathrm{d}V/\mathrm{d}t\)), 그리고 \(t_{\text{max}}\)는 전류 피크가 발생하는 시간입니다. 이 관계는 실험적으로 관측된 과도 응답을 관심 대상의 기본 수송 파라미터와 연결합니다.
3. CELIV의 한계
CELIV 이동도 방정식의 유도는 몇 가지 단순화된 가정에 의존합니다. 첫째, 한 종류의 전하 캐리어만 수송을 지배한다고 가정하므로 전압 램프가 인가될 때 오직 한 종의 캐리어만 이동합니다. 둘째, 캐리어는 균일하게 분포되어 있으며, ??의 왼쪽에 표시된 이상화된 경우와 유사하게 소자 밖으로 부드럽게 미끄러져 나간다고 가정합니다. 셋째, 추출 중 재결합 과정은 무시할 수 있다고 간주하므로 유의미한 캐리어 손실이 발생하지 않습니다. 실제로 이러한 조건은 거의 충족되지 않습니다. ??의 오른쪽에 나타난 것처럼, 전하 캐리어는 종종 얕은, 중간, 깊은 트랩을 점유하며, 이들의 지연된 방출은 넓어지거나 왜곡된 추출 과도 응답을 만들어 냅니다. 많은 실험 기법과 마찬가지로, 따라서 CELIV를 사용하여 추출한 이동도는 진정한 미시적 값이라기보다 근사값으로 간주해야 하며, 실제로 겉보기 이동도는 과도 응답 자체 동안 변화할 수 있습니다(doi:10.1063/1.4818267). 이 한계는 특히 무질서 유기 반도체에서 중요한데, 에너지 무질서와 트래핑이 수송에 강하게 영향을 주기 때문입니다. 반대로 페로브스카이트 재료는 더 높은 이동도와 더 낮은 트랩 밀도를 보이는 경우가 많아 표준 CELIV 해석이 더 강건하고 해석하기 쉬운 편입니다.
💡 이해 확인: 간단한 예제에 CELIV 이론을 적용해 보십시오.
두께가 d = 200 nm인 소자,
전압 램프 속도가 A = 2 × 106 V/s,
전류 피크가 tmax = 5 µs에서 관측되었다고 가정해 봅시다.
CELIV 이동도 방정식을 사용하여 캐리어 이동도를 추정하십시오.
답 보기
CELIV 이동도는 다음과 같이 주어집니다
\[
\mu = \frac{2 d^{2}}{3 A t_{\text{max}}^{2}}
\]
값 대입:
d = 200 × 10-9 m,
A = 2 × 106 V/s,
tmax = 5 × 10-6 s.
\[ \mu = \frac{2 (200 × 10^{-9})^{2}}{3 (2 × 10^{6}) (5 × 10^{-6})^{2}} \approx 1.1 × 10^{-4} \; \text{cm}^2/\text{Vs} \]
이 간단한 예는 CELIV 과도 응답의 피크 시간이 어떻게 이동도 추정값으로 변환될 수 있는지를 보여줍니다.
4: OghmaNano에서 CELIV 시뮬레이션 실행
New simulation을 클릭합니다. 그러면 사용 가능한 소자 유형의 라이브러리가 열리며, 이는 ??에 나와 있습니다. Perovskite cells를 더블클릭하면(빨간색으로 강조됨) 페로브스카이트 예제 폴더가 열립니다. MAPbI₃ device, Perovskite solar cell, 그리고 전용 CELIV 데모를 포함한 미리 설정된 시뮬레이션 목록이 ??에 표시된 것처럼 나타납니다. 이 튜토리얼에서는 Perovskite solar cell – CELIV example을 선택합니다. 메시지가 표시되면 쓰기 권한이 있는 폴더에 시뮬레이션을 저장합니다.
Perovskite solar cell – CELIV example을 선택하면 메인 시뮬레이션 창이 열립니다 (참조: ??). 계산을 시작하려면 Run simulation (파란색 재생 아이콘)을 클릭하거나 F9를 누르십시오. 그러면 OghmaNano는 시간 의존 drift–diffusion 방정식을 풀고 CELIV 과도 응답을 생성합니다.
jv.csv (JV 곡선 데이터), time_j.csv (시간에 대한 추출 전류),
time_v.csv (시간에 대한 인가 전압), 그리고 optical_output (필드 분포)가 포함됩니다.
내장 플로팅 도구에서 해당 결과를 보려면 파일을 더블클릭하십시오.
Step 5: CELIV 결과 보기
Output 탭을 열고
(??)
time_v.csv를 더블클릭합니다. 그러면 인가된 전압 프로그램 —
즉, CELIV 램프 — 가
??에 표시된 것처럼 플롯됩니다.
다음으로 time_j.csv를 더블클릭하여 추출 전류 과도 응답을 표시합니다.
이는 ??와 같습니다.
이 과도 응답은 익숙한 CELIV 형태를 가지지만 **반전된** 것처럼 보입니다:
OghmaNano에서는 부호 규약상 소자에서 나가는 전류가 음수입니다.
??의 램프가 인가되면
전하가 소자 밖으로 끌려 나오므로 소자는 음의 전류를 전달합니다(추출 피크).
반대 부호를 갖는 이후 피크는 (
??에서 양수)
바이어스가 되돌아가고 과도 응답이 초기 상태로 복귀한 후 전하가 다시 소자 안으로 급격히 들어올 때 발생합니다.
??에 있는 플롯은 실험에서 일반적으로 사용하는 표현과 맞추기 위해 시뮬레이션 전류에 –1을 곱한 뒤의 CELIV 과도 응답을 보여줍니다. 이 형태에서 기준선은 소자의 기하학적 정전용량에 해당하며, 날카로운 양의 피크는 인가된 전압 램프 동안 이동 가능한 전하 캐리어가 추출되는 것을 나타냅니다. 이 피크의 크기와 발생 시간은 캐리어 이동도와 밀도를 결정하는 데 사용되는 핵심 물리량입니다. 과도 응답을 이러한 표준 방향으로 제시하면 OghmaNano 시뮬레이션을 발표된 CELIV 측정 결과와 직접 비교하기가 더 쉬워집니다. 이동도를 추출하려면 전류가 최대값에 도달하는 시간(\(t_\text{max}\))을 식별하고, 이 값을 CELIV 이동도 방정식 \(\mu = \tfrac{2d^2}{3At_\text{max}^2}\)에 대입하십시오. 여기서 \(d\)는 active 레이어 두께이고 \(A\)는 인가된 전압 램프 속도입니다.
div style="clear: both;">🧪 예제 계산:
??에서 CELIV 피크 시간을 읽고 이동도를 추정하십시오.
active 레이어 두께가 d = 600 nm, 전압 램프 속도가
A = 5 × 105 V/s, 측정된 피크 시간이
tmax = 3 µs라고 가정합니다. 어떤 이동도가 얻어집니까?
답 보기
CELIV 관계식 \[ \mu = \frac{2 d^{2}}{3 A t_{\text{max}}^{2}} \] 를 사용합니다. 여기서 \(d = 600\times10^{-9}\,\mathrm{m}\), \(A = 5\times10^{5}\,\mathrm{V\,s^{-1}}\), \(t_{\text{max}} = 3\times10^{-6}\,\mathrm{s}\)입니다.
\[ \mu = \frac{2(600\times10^{-9})^{2}}{3\,(5\times10^{5})\,(3\times10^{-6})^{2}} \approx 5.3\times10^{-8}\ \mathrm{m^{2}\,V^{-1}\,s^{-1}} \] 이를 \(\mathrm{cm^{2}\,V^{-1}\,s^{-1}}\)로 변환하면(\(1\ \mathrm{m^{2}} = 10^{4}\ \mathrm{cm^{2}}\)): \[ \mu \approx 5.3\times10^{-4}\ \mathrm{cm^{2}\,V^{-1}\,s^{-1}}. \]
이 값은 무질서 유기 반도체에 대해 보고되는 전형적인 범위에 있습니다. (만약 \({t_\text{max}}\)가 램프 속도나 조명에 따라 변한다면 \(\mu\)도 그에 따라 이동하는 것을 보게 될 것입니다.)
👉 다음 단계: 다음 CELIV (C-LIV) 시뮬레이션 섹션 — Part B — 로 이동하여 다음을 수행합니다: 전압 램프 속도를 조정하고, 광 세기를 실험하며, 전자/정공 이동도 변화가 CELIV 과도 응답을 어떻게 재형성하는지 탐구합니다.