엑시톤 시뮬레이션 튜토리얼

1. 소개

이 튜토리얼에서는 유기 광전(OPV) 소자에서 엑시톤 동역학을 시뮬레이션하는 방법을 살펴봅니다. 유기 반도체에서 광자가 흡수되면, 이는 엑시톤으로 알려진 결합된 전자–정공 쌍을 생성합니다. OPV 소자에서는 이러한 엑시톤이 도너–억셉터 계면으로 드리프트하거나 확산해야 하며, 그곳에서 자유 전하 캐리어로 해리될 수 있습니다. 이 과정은 엑시톤 해리 방정식을 사용하여 모델링할 수 있으며, 이는 엑시톤 수송과 이들이 결국 전자와 정공으로 전환되는 과정을 모두 설명합니다. OghmaNano에서는 이러한 엑시톤 동역학이 전체 광물리적 그림을 포착하기 위해 표준 drift–diffusion 방정식 위에 도입됩니다. 엑시톤 모델링이 언제 필요한지에 대한 더 자세한 논의는 (그리고 언제 생략할 수 있는지는) 엑시톤과 geminate 재결합에 제공되어 있습니다.

2. 시작하기

파일 리본의 New simulation 탭에서 클릭하여 New simulation 창을 엽니다 (그림 1a 참조). 이 창에서 미리 구성된 예제의 여러 범주를 찾을 수 있습니다. Exciton simulations 폴더를 더블클릭하여 사용 가능한 엑시톤 관련 예제의 하위 메뉴를 엽니다 (그림 1b 참조). 이 튜토리얼에서는 Exciton device 템플릿을 사용할 것이며, 이는 OPV 구조에서 엑시톤 생성, 수송, 해리를 시뮬레이션하기 위한 단순한 시작점을 제공합니다.

3. 시뮬레이션 실행 및 엑시톤 출력 확인

소자가 열리면 메인 창에서 적층 구조를 볼 수 있습니다. 시뮬레이션을 실행하기 전에 Electrical 탭으로 가서 Exciton solver 버튼이 눌려 있는지 확인하십시오 — 일반적으로 기본값으로 그렇지만, 확인하는 것이 좋습니다 (??). 이렇게 하면 모델에서 엑시톤 동역학이 활성화됩니다. 이 작업이 완료되면 File 메뉴로 돌아가 Run simulation 버튼을 눌러 모델을 실행합니다 (??).

시뮬레이션이 끝나면 Output 탭에서 결과를 확인할 수 있습니다 (??). 여기에는 엑시톤 관련 출력을 포함하여 디스크에 기록된 모든 파일이 나열됩니다. jv.csv를 더블클릭하면 JV 곡선을 플롯할 수 있습니다 (??). 이 플롯은 일반적인 소자 시뮬레이션에서 얻는 표준 JV 곡선처럼 보이지만, 이 경우에는 Exciton solver가 활성화되어 있으므로 엑시톤 동역학이 기본 물리에 포함되어 있습니다.

출력 디렉터리 내의 exciton_output 폴더를 더블클릭하면 엑시톤 솔버의 상세 결과에 접근할 수 있습니다 (?? 참조). 이 디렉터리에는 위치의 함수로서의 상수와 엑시톤 동역학으로부터 도출된 계산량을 포함하여 모든 솔버 출력이 들어 있습니다. 예를 들어, exciton.csv를 더블클릭하면 소자 두께 전반의 엑시톤 분포 플롯이 생성됩니다 (?? 참조). 마찬가지로 Gn.csv는 위치의 함수로서 전자 생성률을, Gp.csv는 정공 생성률을 보여줍니다 (?? 참조). 단순한 1D 소자의 경우 이 두 생성률 파일은 사실상 동일합니다.

4. 엑시톤 수송 방정식과 매개변수

소자 내 엑시톤 분포는 다음 엑시톤 수송 방정식에 의해 지배됩니다:

\[ \frac{\partial X}{\partial t} = \nabla \!\cdot \!\big(D\,\nabla X\big) + G_{\mathrm{optical}} - k_{\mathrm{dis}}\,X - k_{\mathrm{FRET}}\,X - k_{\mathrm{PL}}\,X - \alpha\,X^{2} \]

여기서 \(X(\mathbf{r},t)\)는 엑시톤 밀도 (m\(^{-3}\)), \(D\)는 엑시톤 확산 계수 (m\(^2\)s\(^{-1}\))입니다; \(G_{\mathrm{optical}}\)은 광학 모델이 제공하는 국소 엑시톤 생성률(흡수된 광자에 비례)이고; \(k_{\mathrm{dis}}\)는 자유 전하로의 해리 속도, \(k_{\mathrm{FRET}}\)는 Förster 공명 에너지 전달 속도이며; \(k_{\mathrm{PL}}\)는 복사 감쇠 속도, \(\alpha\)는 엑시톤–엑시톤 소멸 계수 (m\(^3\)s\(^{-1}\))입니다. 이 모델이 활성화되면 전자 drift–diffusion 방정식은 \(G = k_{\mathrm{dis}}\,X\)로 정의된 생성 항을 사용하며, 설정된 경우 계면 영역으로 제한될 수 있습니다.

이 모든 매개변수는 메인 메뉴에서 접근 가능한 Electrical parameter editor에서 보고 수정할 수 있습니다 (?? 참조). 엑시톤 관련 필드는 Excitons 제목 아래에 그룹화되어 있습니다. 이는 아래 표에 요약되어 있습니다.

5. 솔버가 시뮬레이션 과정에 맞물리는 방식

엑시톤 솔버가 꺼져 있을 때 (??의 상단), transfer-matrix 광학은 캐리어 생성률을 계산하여 바로 drift–diffusion 솔버에 전달합니다. 엑시톤 솔버가 켜져 있을 때 (??의 하단), 광학은 대신 엑시톤 생성률을 엑시톤 솔버에 제공합니다. 그런 다음 엑시톤 솔버는 Electrical parameter editor에 설정된 매개변수를 사용하여 그 집단을 확산, 전달(FRET), 해리, 복사 감쇠, 소멸시키며 진화시키고, 최종 캐리어 생성률을 drift–diffusion 솔버에 출력합니다. 즉, 엑시톤 솔버는 광학과 전기 수송 사이에 자신을 삽입하여 광학 설정이나 drift–diffusion 방정식을 변경하지 않고도 엑시톤 물리를 모델링할 수 있게 합니다.

엑시톤이 활성화된 시뮬레이션의 연산 순서는 ??에 설명되어 있습니다. 먼저 광학 솔버가 실행되어 소자 전반에 걸쳐 슬라이스별로 광자 흡수 프로파일을 계산합니다. 다음으로 엑시톤 솔버가 실행되며, 생성된 엑시톤은 솔버가 수렴할 때까지 전달되거나 해리되며, 보통 수십 단계 내에 수렴합니다. 마지막으로 엑시톤 생성률은 transfer matrix 모델이 계산한 대로 소자 깊이와 파장의 함수로서 조사됩니다. 이러한 출력들은 함께 광학 모델을 엑시톤 솔버와, 궁극적으로 drift–diffusion 방정식과 연결하는 파이프라인을 형성합니다.