전달행렬법(TMM)을 이용한 페로브스카이트 태양전지의 광학 시뮬레이션

전달행렬법(TMM)을 활용하여 다층 구조의 페로브스카이트 태양전지에서 광흡수, 간섭, 반사, 투과 및 전하 캐리어 생성 과정을 모델링하는 방법을 배워보십시오. 본 튜토리얼에서는 박막 광학 시뮬레이션을 통해 광전지 소자의 성능을 분석하고 OghmaNano 환경에서 달성 가능한 최대 광전류 값을 산출하는 절차를 설명합니다.

1. 서론

이 튜토리얼에서는 전달행렬법(TMM)을 이용하여 다층 구조의 페로브스카이트 태양전지의 광학적 특성을 시뮬레이션하게 됩니다. 해당 시뮬레이션을 통해 광장 분포, 흡수된 광자 밀도, 전하 캐리어 생성 분포, 그리고 장치의 반사 및 투과 스펙트럼을 산출할 수 있습니다.

광학 시뮬레이션은 박막 태양전지를 이해하는 데 가장 유용하고 신뢰성 높은 방법 중 하나입니다. 이는 주로 굴절률의 실수부와 허수부인 \(n\) 및 \(k\)에 의존하기 때문인데, 이 값들은 문헌에서 널리 찾아볼 수 있으며 대체로 정확하게 측정이 가능합니다. 반면 전기적 드리프트-확산 시뮬레이션에는 보다 많은 재료 및 수송 관련 파라미터가 필요하며 이들 값을 정확히 결정하는 것은 어려운 경우가 많습니다. 따라서 광학 시뮬레이션은 페로브스카이트 및 유기 태양전지에 대한 물리적 특성을 파악하는 데 훨씬 유용한 정보를 제공합니다.

광생성 분포 전달행렬법(TMM) 각 광학층 내부에서 맥스웰 방정식을 해석적으로 풀고 물질 계면에서 전자기파를 연속적으로 연결함으로써 이 방법은 매우 빠르고 수치적으로 안정적이며, 강한 간섭 효과가 나타나는 다층 구조의 광전지 장치에 특히 적합하다. 수학적 유도 과정에 대한 상세한 설명은 다음에 제시되어 있다. 전달행렬법(TMM)의 유도 과정.

태양광 소자에서 광생성 분포는 해당 태양전지가 생성할 수 있는 광전류의 상한값을 나타낸다. 광학 시뮬레이션 결과 흡수가 미미하거나 생성률이 낮게 예측된다면, 전하 수송 특성과 무관하게 소자는 높은 효율을 달성할 수 없다. OghmaNano에서는 광생성 분포를 전기적 드리프트-확산 해석기와 직접 연동시켜 자기일관적인 전기광학 시뮬레이션을 수행할 수 있다. 관련 예시로는 무반사 코팅, 파브리-페로 공진기, 유기 태양전지, OLED 광방출 특성 등이 있다.

2. 새로운 시뮬레이션 생성하기

새로운 시뮬레이션 창을 열고 ??에 표시된 페로브스카이트 셀 카테고리를 선택합니다. 그 다음 ??에 표시된 페로브스카이트 태양전지(MAPI) 예제를 선택합니다.

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

이 소자는 투명 전극, 전자 수송층, 페로브스카이트 흡수층, 정공 수송층 및 금속 후면 전극을 포함하는 다층 평면형 페로브스카이트 태양전지로 구성된다. 이러한 층들 간의 광학적 간섭 현상은 소자 내부의 전자기장 분포를 크게 변화시키며, 그 결과 국부적인 전하 캐리어 생성률에도 영향을 미친다.

??에 표시된 레이어 편집기를 열려면 Device structure 탭의 왼쪽에 있는 레이어 편집기 버튼을 클릭하십시오. 해당 소자는 투명 전도성 산화물층, TiO\(_2\) 전자 수송층, MAPI 페로브스카이트 흡수층, Spiro 정공 수송층, 그리고 금 금속 접촉층으로 구성됩니다. 각 레이어의 두께는 Thickness 열에서 직접 수정할 수 있습니다.

각 층에 할당된 광학 재료는 재료명 오른쪽에 위치한 점 세 개 버튼을 클릭하여 변경할 수 있습니다. 이때 재료 데이터베이스가 열리며, 실험적으로 측정된 광학 상수값인 \(n\) 및 \(k\)를 각 층에 선택적으로 지정할 수 있습니다. 재료 데이터베이스에 관한 보다 상세한 정보 및 OghmaNano에 새로운 광학 재료를 추가하는 방법은 여기에서 확인할 수 있습니다.

3. 전달행렬 해석기 열기

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

전달행렬법(TMM) 버튼을 클릭하여 광학 시뮬레이션 편집기를 열고, 이후 ‘광학 시뮬레이션 실행’ 버튼(재생 버튼)을 눌러 광학 시뮬레이션을 수행하십시오.

전달행렬 해석기는 다음 값들을 계산한다:

• 소자 내부의 광자 밀도 분포.
• 흡수된 광자 분포
• 생성률 분포
• 반사된 광출력
• 투과된 광출력

전달행렬법(TMM)은 다층 구조를 해석적으로 풀어내므로 시뮬레이션이 대개 거의 즉시 완료된다.

시뮬레이션이 완료되면 광학 시뮬레이션 편집기가 여러 출력 탭과 함께 표시됩니다. 첫 번째 탭인 광자 분포 탭에는 ??에 나타난 바와 같이 소자 내부의 광자 밀도가 표시됩니다.

수평축은 소자 내부의 위치를 나타내고 수직축은 파장을 나타냅니다. 밝게 표시된 영역은 광장 세기가 높은 부분입니다. 그래프에 나타난 수직 경계선들은 굴절률이 변하는 소자의 각 층을 나타냅니다. 전방 및 후방으로 전파되는 전자기파 간의 간섭 현상으로 인해 다층 구조 내부에 정상파 패턴이 형성됩니다.

‘흡수된 광자 분포’ 탭을 선택하면 ??에 나타난 흡수된 광자 분포도가 생성됩니다. 이 그래프는 장치 내부에서 실제로 광자가 흡수되는 위치를 보여주며, 이는 전체 광자 밀도 분포와는 차이가 있습니다.

국소적인 광흡수율은 \[ G_{\mathrm{opt}}(y,\lambda)=\alpha(y,\lambda)\,\phi(y,\lambda) \]로 표현되며, 여기서 \(\alpha\)는 광흡수계수이고 \(\phi\)는 광자유속밀도입니다. 광흡수는 주로 전자-정공 쌍이 생성되는 페로브스카이트 흡수층 내부에서 일어납니다.

파장이 짧은 경우 광흡수계수가 크기 때문에 광장은 전면 계면 근처에서 급격히 감쇠한다. 파장이 길어질수록 광흡수계수는 감소하여 광자가 구조 내부로 더 깊이 침투할 수 있게 된다. 따라서 이 그림들은 파장 및 소자 내 위치에 따른 광흡수 분포를 나타낸다.

많은 경우에 파장 의존성을 단일 공간 분포 형태로 압축하는 것이 유용하다. 이 분포는 ‘생성률’ 탭을 통해 확인할 수 있으며, 여기서는 파장에 걸쳐 적분된 전하 캐리어 생성률이 페로브스카이트 소자 내부의 위치에 따라 나타난다.

4. 반사 및 투과 스펙트럼

반사된 광출력은 ??에 나타나 있으며, 투과된 광출력은 ??에 나타나 있습니다. 반사 및 투과 현상은 비교적 간단한 분광 기술을 이용해 정확하게 측정할 수 있기 때문에 이러한 스펙트럼들은 광학 시뮬레이션을 실험적 측정값과 비교 교정하는 데 가장 유용한 수단 중 하나입니다.

반사 손실은 활성층으로 결합되는 광량을 감소시키는 반면, 투과된 빛은 장치 내에서 흡수되지 않고 완전히 통과한 광자에 해당합니다. 고성능 태양전지에서는 일반적으로 이러한 두 값을 모두 최소화하여 입사하는 광자의 대부분이 활성 영역 내에서 흡수되어 전하 캐리어로 변환되도록 하는 것이 목표입니다.

반사 및 투과 스펙트럼은 다층 구조의 광학적 특성을 직접적으로 파악하는 데에도 활용된다. 스펙트럼에서 나타나는 급격한 진동 현상은 대개 박막 간 간섭 현상이나 광학적 공진 효과, 혹은 층 두께에 따른 공명 현상의 징후이다. 층 두께나 굴절률에 미미한 변화만으로도 이러한 특성들이 크게 변할 수 있으므로, 광학 스펙트럼은 소자 구조를 분석하는 데 매우 민감한 지표가 된다.

전달행렬법(TMM)은 박막 태양전지에서 광흡수 및 광전류 생성을 극대화하기 위해 층 두께, 굴절률 분포 및 무반사 코팅을 최적화하는 데 특히 유용하다.

광학 시뮬레이션에서는 입사 광에너지가 반사, 투과 및 소자 내부에서의 흡수로 분배됩니다. 수치적 오차를 무시할 경우 에너지 보존법칙에 따라 \[ R(\lambda)+T(\lambda)+A(\lambda)=1 \]이 성립해야 하며, 여기서 \(R\)는 파장 \(\lambda\)에서의 반사 광에너지, \(T\)은 투과 광에너지, \(A\)는 흡수된 광에너지입니다. 이 관계식은 전달행렬법 기반 시뮬레이션의 타당성을 검증할 때 매우 유용한데, 계산된 광스펙트럼의 일관성을 직접적으로 확인할 수 있기 때문입니다.

5. 달성 가능한 최대 광전류

광학 시뮬레이션이 완료되면 생성된 파일들이 ??에 표시된 출력 탭에 나타납니다. 출력 디렉터리에는 광자 밀도 맵, 흡수된 광자 분포, 반사 및 투과 스펙트럼, 그리고 파장별 생성률 데이터가 포함되어 있습니다.

sim_info.dat 파일을 더블클릭하면 ??에 나타나는 시뮬레이션 정보 창이 열립니다. 이 파일에는 광흡수만을 기반으로 예측된 최대 광전류 밀도 및 광전류에 관한 정량적 정보가 포함되어 있습니다.

이 값들은 오직 광흡수만을 기준으로 산출된 최대 가능 광전류량에 해당하며, 생성된 모든 전자-정공 쌍이 재결합이나 수송 손실 없이 성공적으로 수집된다고 가정합니다. 실제로는 오제 재결합, 캐리어 트래핑 현상, 전하 캐리어 이동도의 한계 및 불완전한 전하 추출 과정으로 인해 실제 전류값은 이 광학적 상한값보다 낮게 나타납니다.

전달행렬법(TMM) 기반 시뮬레이션은 전기적 드리프트-확산 및 재결합 효과가 장치 모델에 반영되기 전에 달성 가능한 최대 광전류값에 대한 유용한 기본적 추정치를 제공한다.

광학 모델에서 얻을 수 있는 총 광전류는 소자 전체에 걸쳐 전하 캐리어 생성 분포를 적분함으로써 결정된다. 단순화된 형태에서 최대 광전류 밀도는 \[ J_{\mathrm{ph}} = q \int G(x)\,dx \]로 표현되며, 여기서 \(q\)는 기본 전하량이고 \(G(x)\)은 국부적인 부피당 생성률이다. 이 식은 박막형 태양전지 내부에서의 광생성 분포가 광전지 성능 예측에 얼마나 중요한지를 잘 보여준다.

6. 조사 광원 변경하기

??에 표시된 광학 탭으로 돌아가서 광원 메뉴를 클릭하십시오. 그러면 ??에 나타나는 광원 편집기가 열립니다. 광학 조명원 설정에 관한 보다 상세한 정보는 여기에서 확인할 수 있습니다.

기본적으로 시뮬레이션에서는 AM1.5G 태양광 스펙트럼이 사용됩니다. 현재 스펙트럼은 약 900 nm까지만 나타나는 것을 확인할 수 있습니다. 이러한 범위는 대부분의 페로브스카이트 시뮬레이션에 충분한데, 페로브스카이트 흡수층은 장파장 영역에서 흡수 능력이 상대적으로 약하며 대부분의 광전류 생성은 스펙트럼의 가시광선 영역에서 일어나기 때문입니다.

스펙트럼 드롭다운 메뉴를 이용하여 조사원을 AM1.5G에서 크세논으로 변경합니다. 그 결과 ??에 나타나는 스펙트럼은 여전히 절단된 형태인데, 이는 광학 파장 메쉬의 범위가 아직 크세논 방출 스펙트럼 전체를 포착하기에 충분하지 않기 때문입니다.

파장 범위를 확장하려면 광학 설정 화면으로 돌아가 ??에 표시된 광학 메쉬 편집기를 열어야 합니다. 파장 종료값을 약 1500 nm로 설정한 뒤 재구성 버튼을 클릭하십시오. 이후 광원 편집기로 돌아가면 ??에 나타나는 전체 크세논 스펙트럼을 확인할 수 있습니다.

광학 파장 메쉬는 전달행렬법(TMM) 시뮬레이션 시 사용되는 스펙트럼 해상도와 파장 범위를 모두 결정합니다. 부적절하게 좁은 파장 범위를 선택할 경우, 특히 광대역 조명원이나 적외선 영역의 광학적 특성을 모델링할 때 광흡수 및 생성률 분포가 불완전하게 산출될 수 있습니다.

조사 광원을 변경한 후 광학 시뮬레이션을 다시 실행하십시오. 이제 생성률 탭으로 돌아가면 ??에 나타난 바와 같은 광생성 분포가 산출됩니다. 크세논 램프의 스펙트럼 분포는 AM1.5G 태양광과 상당히 다르기 때문에, 소자 내부의 광투과 깊이 및 전하 캐리어 생성률 분포 역시 크게 변화하게 됩니다.

실험실에서 사용되는 광원의 스펙트럼을 시뮬레이션 조건과 정확히 일치시키는 것은 계산 결과와 실제 태양광 소자 성능을 비교하는 데 매우 중요하다. 실제로 대부분의 태양광 시뮬레이터는 이상적인 AM1.5G 스펙트럼과 상당한 차이를 보인다. 이에 따라 OghmaNano에서는 광원 데이터베이스 편집기를 통해 사용자 정의 조명 스펙트럼을 임포트하는 기능도 제공한다.

7. 광생성을 전기적 모델에 연동시키기

전달행렬 해석기에 의해 산출된 광생성 분포는 드리프트-확산 해석기 내에서 전기적 생성 항으로 변환된다.

광생성 분포와 드리프트-확산 전기 모델 간의 변환 과정은 ??에 표시된 광자 효율 파라미터에 의해 제어된다. 이 파라미터는 광학 시뮬레이션 편집기 내의 구성 버튼을 눌러 설정할 수 있다.

드리프트-확산 해석기로 입력되는 전기적 생성률은 대략적으로 \[ G_{\mathrm{dd}}=\eta\,G_{\mathrm{opt}} \] 식으로 표현되며, 여기서 \(G_{\mathrm{opt}}\)는 광학적 전달행렬법에 의해 계산된 광생성 분포이고 \(\eta\)는 광자 효율 계수입니다.

물리적으로 \(\eta\)은 흡수된 광자가 최종적으로 전류 생성에 기여하는 자유 전자-정공 쌍을 만들어낼 확률을 나타냅니다. 이상적인 태양전지 시스템에서는 이 값이 1에 가까워질 수 있으나 실제 소자에서는 쌍생 재결합, 엑시톤의 불완전한 해리, 부차적 흡수 현상, 계면 손실 또는 극초단 시간 내에 일어나는 재결합 과정 등으로 인해 이 값이 감소하는 경우가 많습니다.

이러한 형태의 매개변수들은 광흡수량만을 계산하는 광학 시뮬레이션 결과와 전하 분리 과정의 미시적 물리 현상을 연결하는 데 유용하므로 광전지 소자 모델링에서 널리 활용된다. 따라서 광자-효율 계수는 광학 모델과 전기적 수송 모델 간의 효과적인 연결 고리 역할을 한다.

실제로 \(\eta\)은 실험 데이터와 시뮬레이션 결과를 정합시키는 과정에서 작은 피팅 파라미터로도 활용될 수 있다. 특히 재료의 광학 상수인 \(n\) 및 \(k\), 층 구조 형태, 혹은 계면 손실 등이 정확히 알려져 있지 않은 경우에 그러하다. 다만 이 값이 1에서 크게 벗어날 경우, 이는 효율 파라미터에 모든 물리적 현상이 흡수된 것이 아니라 전기적 모델에 추가적인 물리 메커니즘을 명시적으로 포함시켜야 함을 시사한다.

8. 광학 스냅샷 검토하기

출력 탭으로 돌아가서 optical_snapshots 디렉터리를 더블클릭합니다. 그러면 다층 구조의 소자 내 특정 파장에서 각각의 광학량을 확인할 수 있는 광학 스냅샷 뷰어가 열립니다.

그래프 아래에 위치한 슬라이더를 이용하면 파장을 변경할 수 있으며, 창 하단의 드롭다운 메뉴에서는 다양한 광학량을 선택할 수 있다. 선택 가능한 항목으로는 광자 밀도, 흡수된 광자 밀도, 광흡수계수, 생성률, 그리고 전방 및 후방으로 전파되는 전기장 성분 등이 있다.

??에 제시된 광자밀도 분포는 소자 내부의 국부적인 전자기장 세기를 나타낸다. 다층 구조 내 각 계면에서 발생하는 반사 현상으로 인해 정상파 간섭 패턴이 명확하게 관찰된다.

??에 나타난 광자 흡수 분포는 구조 내부에서 광에너지가 실제로 흡수되는 위치를 나타냅니다. 이렇게 흡수된 광에너지는 ??에 표시된 전하 캐리어 생성 분포로 전환됩니다.

??에 나타난 광흡수계수 분포는 각 물질층별로 파장에 따라 변하는 광흡수계수 \(\alpha\)를 나타낸다. \(\alpha\) 값이 더 큰 층일수록 광자를 더 강하게 흡수하여 광장의 감쇠 속도가 더 빨라진다.

?? 및 ??에 나타난 전진 및 후진 방향으로 전파되는 전기장 성분들은 전달행렬법(TMM) 내부에서 사용되는 두 개의 역방향으로 전파되는 전자기파에 해당합니다. 이들 파동 간의 간섭 작용으로 인해 소자 전체에 걸쳐 정상파 형태의 광학 패턴이 형성됩니다.

이러한 광학적 스냅샷들은 박막형 태양전지 구조 내에서 빛이 전파되고 반사되며 간섭하고 흡수되는 과정을 상세히 보여줍니다. 따라서 광학적 공진기 효과를 이해하고 다층 구조의 태양전지 설계를 최적화하는 데 특히 유용합니다.

9. 요약

이 튜토리얼에서는 전달행렬법(TMM)을 활용하여 OghmaNano 내에서 다층 구조의 페로브스카이트 태양전지의 광학적 특성을 시뮬레이션했습니다. 광자 밀도 분포, 흡수된 광자 분포, 파장에 따른 광생성 분포, 반사 및 투과 스펙트럼, 그리고 광학 모델로부터 도출된 최대 광전류 추정값 등을 분석했습니다.

또한 조사 광스펙트럼 및 광파장 메쉬의 변경이 소자 내부의 전하 캐리어 생성 분포에 어떤 영향을 미치는지도 살펴보았습니다. 이러한 효과들은 실험 측정값과 시뮬레이션 결과를 비교할 때 특히 중요합니다. 왜냐하면 박막형 태양전지 소자들은 광학적 간섭 현상, 조사 조건 및 층 두께에 매우 민감하기 때문입니다.

광학적 전달행렬법 기반 시뮬레이션은 박막형 태양전지의 특성을 이해하고 최적화하는 데 활용될 수 있는 가장 강력한 도구 중 하나입니다. 이 방법은 비교적 적은 양의 재료 파라미터만을 필요로 하면서도 광흡수 및 광전류 한계에 관한 신뢰도 높은 물리적 정보를 제공합니다. OghmaNano에서는 전달행렬 해석기가 전기-광학적 드리프트-확산 시뮬레이션의 광학적 기반을 구성하며, 이를 통해 광생성 분포를 전하 수송 및 재결합 모델과 직접 연동할 수 있습니다.