전달행렬법(TMM)을 이용한 유기 태양전지의 광학 시뮬레이션

1. 서론

이 튜토리얼에서는 전달행렬법(TMM)을 활용하여 박막형 유기태양전지(OPV)의 광학적 특성을 시뮬레이션하게 됩니다. 해당 시뮬레이션은 장치 내 광장 분포, 흡수된 광자 상태밀도, 전하 캐리어 생성 분포, 그리고 반사 및 투과 스펙트럼을 계산합니다. OghmaNano 내 전달행렬법 솔버의 구현에 관한 보다 상세한 정보는 여기에서 확인하실 수 있으며, 이 방법에 대한 전체 유도 과정은 여기에 기술되어 있습니다.

광학적 시뮬레이션은 태양전지 성능을 초기 단계에서 파악하는 데 있어 가장 신뢰할 수 있는 방법 중 하나입니다. 이때 필요한 재료 파라미터는 상대적으로 적기 때문입니다. 각 층에 대한 굴절률 \(n\) 및 흡수 계수 \(\alpha\)만이 요구되며, 이들 값은 문헌에서 널리 찾아볼 수 있고 대체로 정확하게 측정이 가능합니다. 이는 전기적 드리프트-확산 모델과 대조적인 특징인데, 후자의 경우 전하 캐리어 이동도, 재결합률, 트랩 상태밀도 및 계면 특성과 같은 파라미터들이 작동 조건에 따라 크게 달라질 수 있으며 고유한 값을 결정하기가 어려운 경우가 많습니다.

유기 태양전지 소자는 활성층의 두께가 보통 수백 나노미터에 불과하기 때문에 광학적 간섭 현상에 매우 민감하다. 각 물질 계면에서의 반사 현상으로 인해 다층 구조 내부에 정상파 패턴이 형성되며, 이로 인해 광학적 전기장 및 광학적 생성 분포는 파장, 층 두께, 그리고 광학 상수에 따라 크게 변동한다. 따라서 본 광학 모델을 통해 소자 내부에서 광자가 어디에 흡수되고 전하 캐리어가 어디에서 생성되는지를 직접적으로 파악할 수 있다.

전달행렬법(TMM)은 각 층 내에서 맥스웰 방정식을 해석적으로 풀고, 층간 경계면에서 전진 및 후진으로 전파되는 전자기파를 연속적으로 연결한다. 이 방법은 가정 사항이 비교적 적고 시간 영역에서의 직접적인 전파 계산을 필요로 하지 않기 때문에 계산 효율성이 높으며, 유한차분 시간영역법(FDTD)과 같은 전파 해석 기법들에 비해 훨씬 빠른 연산 속도를 보인다. 유한차분 시간영역법(FDTD)전달행렬법(TMM) 시뮬레이션은 전기적 수송 및 재결합 모델을 도입하기에 앞서 박막 태양전지 최적화의 초기 단계로 널리 활용된다.

이 광학 모델은 달성 가능한 단락전류밀도의 상한값도 제시해준다 \(J_{SC}\). 다층 구조가 빛을 효율적으로 흡수하지 못한다면, 아무리 이상적인 전기적 수송 특성을 갖춘다 해도 고성능 소자를 만들어낼 수 없다. 광학적 최적화와 관련된 예시들은 광학 필터, 무반사 코팅, 파브리-페로 공진기, 페로브스카이트 태양전지, 그리고 OLED 광학 방출 튜토리얼에서 확인할 수 있다.

OghmaNano에서는 전달행렬법(TMM) 솔버를 통해 산출된 광학적 생성 분포를 전기적 드리프트-확산 모델에 직접 연동할 수 있습니다. 이를 통해 광학적 흡수 현상, 전하 캐리어 생성 및 이동, 재결합 과정, 그리고 전하 추출 현상까지도 자기일관적으로 시뮬레이션할 수 있습니다.

2. 시뮬레이션 생성하기

??에 표시된 새 시뮬레이션 창을 열어주십시오. ‘유기 태양전지’ 카테고리를 선택한 뒤 ??에 나타나는 PM6:Y6_E6_0hrs 예제를 선택하시기 바랍니다.

시뮬레이션 결과를 디스크 상의 적절한 디렉터리에 저장하십시오. 시뮬레이션이 로드되면 ??에 나타난 메인 인터페이스가 표시됩니다.

이 시뮬레이션에는 투명 전극, 전하 수송층, PM6:Y6 활성층 및 금속 후면 전극으로 구성된 완전한 다층 유기태양전지 구조가 포함되어 있습니다. 소자의 두께가 광학적 파장과 유사하기 때문에 간섭 현상이 내부 광장 분포에 크게 영향을 미칩니다.

3. 소자 구조 검토하기

왼쪽 툴바에서 레이어 편집기를 열면 ??에 표시된 다층 구조가 나타납니다.

이 소자에는 인듐 주석 산화물(ITO) 투명 전극, ZnO 전자 수송층, PM6:Y6 활성층, MoO\(_x\) 정공 수송층, 그리고 은으로 제작된 후면 전극이 포함되어 있습니다.

각 층에는 실험적으로 측정된 광학 상수들이 포함되어 있으며 \(n\) 및 \(k\), 이를 통해 파장에 따른 흡수 현상과 간섭 효과를 정확하게 시뮬레이션할 수 있습니다.

박막 유기 태양전지에서는 정렬 간섭 현상으로 인해 광학장이 활성층 내부에서 최대값에 도달하는 경우가 많습니다. 반면, 소멸 간섭 현상으로 인해 특정 영역에서는 광학장이 억제되기도 합니다. 이러한 간섭 효과들은 전하 캐리어 생성 분포가 광학적 전기장 분포와 직접적으로 연관되어 있기 때문에 소자의 최종 전압-전류 특성에 큰 영향을 미칩니다.

4. 전달행렬법 솔버 열기

??에 표시된 광학 탭을 열어주십시오. 이 리본 메뉴에는 전달행렬법(TMM) 제어 기능, 광학 메쉬 설정, 광학 검출기 및 광원 편집기가 포함되어 있습니다.

전달행렬법(TMM) 버튼을 클릭하면 ??에 나타나는 전달행렬법 시뮬레이션 편집창이 열립니다. 광학 시뮬레이션을 실행하려면 ‘광학 시뮬레이션 실행’ 버튼을 누르십시오. 실행 중에는 ??에 표시된 터미널 출력 내용이 나타납니다. 솔버는 다층 유기태양전지 구조에서의 광자밀도 분포, 흡수된 광자밀도, 생성률 분포, 반사 스펙트럼 및 투과 스펙트럼을 계산합니다.

전달행렬법(TMM)은 각 층 내에서 광학적 문제를 해석적으로 풀어내므로, 여러 파장과 다수의 흡수층을 포함하는 소자의 경우에도 계산 속도가 매우 빠릅니다.

5. 광자 밀도와 흡수된 광자

시뮬레이션이 완료되면 ??에 나타난 광자밀도 분포가 표시됩니다. 가로축은 소자 내부의 위치를 나타내고 세로축은 파장을 나타냅니다. 이 그림은 다층 구조 내에서 전파되는 광자의 공간적 분포를 보여주며, 물질 계면에서의 반사로 인해 발생하는 간섭 현상도 포함합니다.

‘흡수된 광자 분포’ 탭을 선택하면 ??에 나타난 흡수된 광자 분포도가 생성됩니다. 장치 내부의 전체 광장 세기를 보여주는 광자밀도 그래프와는 달리, 이 그림은 물질에 의해 실제로 흡수되어 전자-정공 쌍으로 변환된 광자들을 나타냅니다.

국소적인 흡수된 광자 밀도는 \[ G_{\mathrm{opt}}(x,\lambda)=\alpha(x,\lambda)\,\phi(x,\lambda) \]로 표현되며, 여기서 \(\alpha\)는 파장에 따른 광학적 흡수 계수이고 \(\phi\)는 국소적인 광자 밀도입니다. 대부분의 광학적 흡수는 전하 생성이 주로 일어나는 PM6:Y6 활성층 내부에서 발생합니다.

짧은 파장에서는 광학적 흡수 계수가 비교적 크기 때문에 광자가 전면 계면 근처에서 흡수됩니다. 반면 긴 파장에서는 광학적 흡수 길이가 증가하여 광자가 흡수되기 전에 활성층 내부로 더 깊이 침투할 수 있습니다. 또한 박막 공진기 구조 내에서 전파 방향과 후방으로 전파되는 전자기파 간의 간섭 현상으로 인해 뚜렷한 정상파 패턴도 관찰됩니다.

대부분의 경우 전체 스펙트럼 분포를 직접 살펴보는 것보다는 파장에 대해 적분된 총 생성 분포를 검토하는 것이 더 유용합니다. 이는 ??에 제시된 생성률 탭에서 확인할 수 있으며, 여기서는 장치 내부의 위치에 따른 부피당 전자-정공 쌍 생성률이 그래프로 나타나 있습니다.

생성률 그래프 위에 겹쳐진 다색 영역들은 각 소자 층의 전자 에너지 준위를 나타냅니다. 이 에너지 준위들은 재료 데이터베이스에서 가져온 것으로, 주로 소자 구조에 대한 시각적 참고 자료로 활용됩니다. 이들은 광학적 전달행렬법 시뮬레이션 자체에는 직접적인 영향을 미치지 않습니다. 다만 PM6:Y6 활성층의 경우, 전하 이동도, 재결합 및 전하 추출 과정을 시뮬레이션할 때 전기적 드리프트-확산 모델에 의해 해당 에너지 준위들이 사용됩니다.

6. 반사 및 투과 스펙트럼

반사된 광출력은 ??에 나타나 있으며, 투과된 광출력은 ??에 나타나 있습니다. 반사 손실은 활성층으로 결합되는 광출력의 양을 감소시키는 반면, 투과된 빛은 흡수되지 않은 채 장치에서 빠져나가는 광자들을 의미합니다.

소자 내부의 광출력은 반사, 흡수, 투과로 분배되는데 이는 \[ R(\lambda)+A(\lambda)+T(\lambda)=1 \]에 의해 기술된다. 여기서 \(R\), \(A\), \(T\)는 각각 반사, 흡수, 투과된 광출력의 비율을 나타낸다. 이 관계식은 에너지 보존법칙에서 직접 도출되며 전달행렬법 기반 시뮬레이션의 일관성 검증에 유용하다.

파장이 길어질수록 은으로 제작된 후면 전극이 강력한 광반사체 역할을 하므로 투과율은 0에 가까워집니다. 그 결과 대부분의 장파장 광자는 다층 구조 내에서 흡수되거나 전면 표면 쪽으로 다시 반사되어 활성층 내에서의 유효 광경로 길이가 증가하게 됩니다.

실제 유기태양전지 설계에서는 광학적 간섭 현상, 층의 두께 혹은 반사 손실에 있어 비교적 미미한 변화만으로도 단락전류밀도에 유의미한 변화가 발생할 수 있습니다. 이에 따라 전달행렬법(TMM) 기반 시뮬레이션이 활발히 활용되며, 이를 통해 공진기 효과, 광학적 스페이서 및 전극 두께를 최적화하여 활성층 내 광흡수율을 극대화하고 있습니다.

8. 최대 가능 광전류 계산하기

이 값들은 광학 모델에 의해 예측된 최대 가능한 광전류를 나타냅니다. 이때 모든 생성된 전자-정공 쌍이 성공적으로 수집된다는 가정을 전제로 합니다. 오제 재결합, 캐리어 포획 현상, 수송 과정에서의 손실 및 전하 추출 효율 저하와 같은 요소들은 이 순수 광학적 계산에는 포함되지 않았습니다.

최대 단락전류밀도는 파장 및 위치에 걸쳐 광학적 생성률을 적분함으로써 구해지며, 여기서 \(q\)는 기본 전하량이고 \(G(x,\lambda)\)는 소자 내부의 파장 의존적 부피 생성률입니다 \[ J_{SC,\mathrm{max}} = q \int \int G(x,\lambda)\, d\lambda\, dx \].

최적화가 잘 이루어진 소자의 경우 측정된 단락전류밀도 \(J_{SC}\)는 종종 이러한 광학적 상한값과 비교적 유사한 수준을 보입니다. 따라서 광학적 시뮬레이션은 달성 가능한 소자 성능에 대한 매우 유용한 예측값을 제공합니다. 만약 광학 모델에서 최대 광전류값이 낮게 예측된다면, 전기적 수송 모델을 아무리 최적화해도 소자는 충분한 빛을 흡수하지 못하므로 고효율 태양전지를 제작하는 것은 불가능합니다.

9. 조사 스펙트럼 변경하기

??에 표시된 광학 설정 창으로 돌아가서 등대 모양 아이콘으로 표시된 광원 버튼을 클릭합니다. 그러면 ??에 나타나는 광원 편집기가 열립니다.

기본적으로 시뮬레이션에서는 표준 AM1.5G 태양광 스펙트럼이 사용됩니다. 해당 스펙트럼은 장파장 영역에서 잘려 나가는 형태임을 확인할 수 있습니다. 대부분의 유기태양전지 시뮬레이션에서는 활성층이 적외선 영역에서 상대적으로 약하게 흡수하므로, 대부분의 광전류 생성이 가시광선 스펙트럼 범위에서 일어나기 때문에 이러한 설정만으로도 충분합니다. 파장 범위를 제한함으로써 계산 비용을 줄이고 시뮬레이션 속도를 향상시킬 수 있습니다.

실제로는 대부분의 실험적 측정이 실제 AM1.5G 태양광 대신 실험실용 태양광 시뮬레이터를 이용해 수행됩니다. 크세논 램프가 특히 널리 사용되는데, 이는 그 방사 스펙트럼이 태양 스펙트럼의 가시광 영역을 대략적으로 재현하기 때문입니다. 실험에 사용되는 조사 스펙트럼을 정확히 일치시키는 것은 매우 중요합니다. 왜냐하면 스펙트럼 강도의 변화가 광학적 생성 분포와 측정된 소자 성능에 상당한 영향을 미칠 수 있기 때문입니다.

스펙트럼 드롭다운 메뉴를 이용해 조사 광원을 크세논으로 변경한 뒤 Rebuild 버튼을 누르세요. 그러면 ??에 나타나는 스펙트럼이 생성되는데, 광학 파장 메시지가 아직 크세논 방출 스펙트럼 전체를 포착하기에 충분히 길지 않아 장파장 영역에서 스펙트럼이 다소 잘리는 현상이 발생합니다.

파장 범위를 확장하려면 광학 설정 화면으로 돌아가 ??에 표시된 광학 메쉬 편집기를 열어주십시오. 이후 파장 종료값을 약 1500 nm로 설정한 뒤 Rebuild 버튼을 누르면 됩니다. 이제 광원 편집기로 돌아가면 ??에 나타난 전체 크세논 방사 스펙트럼을 확인할 수 있습니다.

조사 스펙트럼을 변경한 뒤 전달행렬법(TMM) 시뮬레이션을 다시 실행합니다. 이후 생성률 탭으로 돌아가면 ??에 나타난 바와 같은 분포가 도출됩니다. 다층 구조 내에서 전하 캐리어가 생성되는 공간적 분포는 기존의 AM1.5G 조건과 상당히 다르는데, 이는 크세논 램프가 다른 스펙트럼 강도 분포를 생성함에 따라 다층 구조 내에서의 광학적 간섭 조건 역시 달라지기 때문입니다.

이러한 효과들은 활성층 두께가 광파장과 비슷한 수준인 박막 유기 태양전지에서 특히 중요합니다. 따라서 조사 스펙트럼, 공진기 내 간섭 현상 또는 흡수 길이의 미세한 변화만으로도 단락전류밀도와 전체 소자 효율에 측정 가능한 변화가 발생할 수 있습니다.

9. 광학적 생성 현상을 전기적 모델에 연동시키기

전달행렬법(TMM) 솔버를 통해 계산된 흡수된 광자 분포는 드리프트-확산 솔버 내에서 전기적 생성 항으로 변환된다.

이 변환 과정은 ??에 표시된 광자 효율 파라미터에 의해 제어됩니다.

물리적으로 이 매개변수는 흡수된 광자가 전기 전류 생성에 기여하는 자유 전자-정공 쌍을 만들어낼 확률을 나타냅니다.

유기 태양전지의 경우, 흡수된 광자는 처음에는 자유 전하 캐리어가 아닌 엑시톤을 형성하므로 이 값은 1보다 현저히 낮을 수 있습니다. 따라서 엑시톤 해리 과정, 동시 재결합 현상 및 전하 전달 손실 등으로 인해 최종 광전류가 크게 감소할 수 있습니다.

10. 광학적 스냅샷 분석

??에 표시된 출력 탭 내에서 optical_snapshots 디렉터리를 더블클릭하십시오. 그러면 광학 스냅샷 뷰어가 열리며, 이를 통해 다층 유기태양전지 구조 전반에 걸쳐 파장별로 광학적 물리량들을 검토할 수 있습니다.

해당 광학 스냅샷들에는 광학 메쉬 내 각 파장에 대해 광자 밀도, 흡수 계수, 전방으로 전파되는 전기장 \(E^+\), 후방으로 전파되는 전기장 \(E^-\), 흡수된 광자 밀도 및 생성률의 공간 분포가 포함되어 있습니다.

??에 나타난 광자밀도 분포는 광학적 간섭 현상이 소자 내부에서 빛의 공간적 분포에 어떤 영향을 미치는지를 보여줍니다. ??에 제시된 흡수 분포는 다층 구조의 파장 의존적 흡수 계수 \(\alpha\)와 일치합니다.

?? 및 ??에 나타난 전기장 분포는 전달행렬법 솔버를 통해 계산된 전방 및 후방으로 전파되는 광파에 대응한다. 이러한 전기장들 간의 간섭 작용으로 인해 브래그 격자와 유사한 정상파 패턴이 형성되며, 이는 박막 유기 태양전지 내부에서 광학적 생성 분포가 극도로 불균일하게 나타나는 원인이 된다.

11. 요약

이 튜토리얼에서는 전달행렬법(TMM)을 활용하여 PM6:Y6 유기 태양전지 내부의 광학적 흡수 현상, 간섭 효과 및 전하 캐리어 생성 과정을 모델링했습니다. 또한 광자 밀도 분포, 흡수된 광자 분포, 생성률 분포, 반사 및 투과 스펙트럼, 그리고 광학 모델에 의해 예측된 최대 광전류 값 등을 분석했습니다.

또한 조사 광스펙트럼을 변경할 경우 다층 구조 내부의 광장 분포 및 광학적 생성 분포가 어떻게 변화하는지도 살펴보았습니다. 이러한 효과들은 박막형 유기 태양전지에서 특히 중요합니다. 왜냐하면 광학적 간섭 현상 및 공진기 구조의 미미한 변화만으로도 장치 성능에 측정 가능한 수준의 변화가 발생할 수 있기 때문입니다.

전달행렬법(TMM)은 다층 구조를 갖는 태양전지의 모델링을 위해 계산적으로 효율적이면서도 물리적으로 신뢰성 높은 프레임워크를 제공한다. 이 방법은 비교적 적은 수의 재료 파라미터만을 필요로 하면서도 광학적 흡수 현상 및 단락전류밀도의 한계에 대한 상세한 정보를 제공하므로, 보다 복잡한 전기적 수송 및 재결합 모델을 도입하기에 앞서 유기태양전지 최적화의 첫 단계로 널리 활용된다.

OghmaNano에서는 전달행렬법(TMM) 솔버를 통해 계산된 광학적 생성 분포를 드리프트-확산 모델에 직접 연동할 수 있습니다. 이를 통해 광학적 흡수, 전하 캐리어 생성, 이동, 재결합 및 추출 과정들을 자기일관적으로 시뮬레이션할 수 있습니다.