비정질 실리콘(a-Si:H) 태양전지(1D) — 박막 p/i/n(결함 제한)

1. 소개

수소화 비정질 실리콘(a-Si:H) 태양전지는 대표적인 박막 광전지 기술입니다. 결정질 실리콘과 달리 a-Si:H는 구조적으로 무질서하며 높은 밀도의 결함 상태와 테일 상태를 포함합니다. 그 결과, 소자 거동은 일반적으로 결함에 의해 재결합 제한되며(Shockley–Read–Hall형 과정) 낮은 이동도로 인해 수송 또한 제한됩니다. 이러한 특성 때문에 a-Si:H는 재결합과 캐리어 밀도가 어떻게 V_oc를 결정하는지, 그리고 박막 광학이 어떻게 J_sc를 형성하는지를 학습하기에 매우 적합한 시스템입니다.

이 튜토리얼은 OghmaNano를 사용하여 1D에서 a-Si:H 태양전지를 시뮬레이션하기 위한 실용적이고 물리 기반의 워크플로를 제공합니다. 이 모델은 깊이 의존적인 광 생성과 SRH형 결함 재결합이 지배하는 물리적 재결합을 포함하여, Poisson 정전기학과 결합된 드리프트–확산 캐리어 수송을 해석합니다. 목표는 표준 출력물—JV 곡선, Voc, Suns–Voc 거동—을 내부 변수(밴드, 준-페르미 준위, 전하, 그리고 유효 재결합 시간)와 연결하는 것입니다.

표준적인 p/i/n 구조를 갖는 1차원 박막 a-Si:H 태양전지를 구성하고 시뮬레이션할 것입니다 (?? 참조). 소자는 다음과 같이 정의됩니다: 구조(전면 → 후면): p a-Si:H / i a-Si:H / n a-Si:H / 금속 접촉. 이 기준 소자를 사용하여 조명하의 JV 곡선을 생성하고, 핵심 성능 지표를 추출하며, 캐리어 밀도에 대한 의존성과 함께 재결합 제한 전압 손실을 식별하기 위해 Voc 및 Suns–Voc 스윕을 수행합니다.

2. 새 시뮬레이션 만들기

시작하려면 OghmaNano 메인 창에서 새 시뮬레이션을 생성합니다. 툴바에서 New simulation 버튼을 클릭합니다. 그러면 시뮬레이션 유형 선택 대화상자가 열립니다 (?? 참조).

시뮬레이션 유형 대화상자에서 Si demos를 더블클릭한 다음, Amorphous silicon (a-Si:H)을 더블클릭합니다 (?? 참조). OghmaNano는 미리 정의된 a-Si:H 박막 태양전지 시뮬레이션을 로드합니다.

메인 창은 (?? 참조) 소자 구조의 3차원 보기를 제공합니다. 마우스를 사용하여 장면을 회전, 이동, 확대하여 기하 구조를 확인할 수 있습니다. 현재 튜토리얼은 1차원 전기 모델을 사용하지만, 3D 보기는 박막 구조와 접촉을 시각화하는 편리한 방법을 제공합니다.

층 표를 열려면 Layer editor 탭을 클릭합니다 (?? 참조). 여기에서 p형 층, 고유 흡수층(i층), n형 층, 후면 접촉을 포함한 수직 소자 구조와, 각 층의 두께 및 할당된 재료를 확인할 수 있습니다.

3. 도핑 프로파일 살펴보기

도핑 프로파일은 조명과 바이어스가 인가되기 전 소자의 정전기 구조를 설정합니다. 박막 비정질 실리콘(a-Si:H) 태양전지에서 도핑은 도핑된 흡수층 내부에 넓은 공핍 영역을 만들기보다는, 주로 캐리어 선택적 접촉을 형성하는 데 사용됩니다. 중심적인 동작 원리는 광생성된 캐리어를 분리하고 수송하는 고유층 전반에 걸친 내장 전기장입니다.

OghmaNano에서 도핑을 보려면 메인 창의 Electrical 리본을 열고 Doping / Ions를 클릭합니다 (?? 참조). 그러면 프로파일 편집기가 열리며 (?? 참조), 여기에는 깊이에 따른 도너 및 억셉터 전하 밀도가 플롯되고 각 층에 할당된 수치 값이 나열됩니다.

이 소자는 표준적인 p+/i/n+ a-Si:H 구조를 사용합니다. 얇고 고농도로 도핑된 p+ 층이 조명되는 전면 표면에 배치되며 (억셉터 밀도는 \(10^{24}\,\mathrm{m^{-3}}\) 정도), 정공 선택적 접촉을 제공하고 전면 밴드 굽힘을 설정합니다. 그 다음에는 훨씬 더 두꺼운 고유 a-Si:H 흡수층이 오며, 이 층은 매우 약하게만 도핑되어 있고(여기서는 \(\sim 10^{20}\,\mathrm{m^{-3}}\)) 광 생성의 대부분을 포함합니다. 소자의 후면에는 얇고 고농도로 도핑된 n+ 층이 있으며 (도너 밀도 역시 \(\sim 10^{24}\,\mathrm{m^{-3}}\)), 금속에 인접한 전자 선택적 후면 접촉을 형성합니다.

프로파일 편집기에서 인식해야 할 핵심 특징은 강한 내장 전기장이 도핑된 영역이 아니라 고유층을 가로질러 형성된다는 점입니다. 광생성된 캐리어는 주로 고유 a-Si:H에서 생성되며 이 내부 전기장에 의해 분리됩니다. 전자는 n+ 후면 접촉 쪽으로 구동되고, 정공은 p+ 전면 접촉 쪽으로 구동됩니다. 도핑된 층은 좋은 전기적 선택성과 접촉 전도도를 유지하면서도 결함 보조 재결합을 최소화하기 위해 얇게 유지됩니다.

4. 전기적 파라미터: a-Si:H에서의 수송, 정전기학, 재결합

메인 창에서 전기적 파라미터 편집기를 엽니다: Device structure → Electrical parameters. 파라미터는 p+, i, n+ 탭을 사용하여 영역별로 정의되며, 이는 a-Si:H 적층의 전면 접촉층, 고유 흡수층, 후면 접촉층에 해당합니다. 이 편집기의 설정은 드리프트–확산 수송 계수, 에너지 정렬, 그리고 JV 및 Suns–Voc 스윕 동안 사용되는 재결합 메커니즘을 결정합니다.

비정질 실리콘에서의 수송은 무질서와 강한 산란으로 인해 낮은 이동도로 특징지어집니다. 따라서 이 데모는 결정질 실리콘보다 훨씬 낮은 이동도를 사용합니다. 고유층에서 (?? 참조), 전자 이동도는 μ_n = 7×10⁻⁴ m²V⁻¹s⁻¹로, 정공 이동도는 μ_p = 7×10⁻⁶ m²V⁻¹s⁻¹로 설정되어 있습니다. 접촉 인접 층은 서로 다른(일반적으로 더 낮은) 이동도를 사용하며 (?? 및 ?? 참조), 이는 도핑/선택적 층이 고이동도 흡수층이라기보다 주로 얇은 추출 영역으로 작용한다는 사실을 반영합니다. 드리프트–확산에서 이동도가 감소하면 캐리어 체류 시간이 증가하고 개방회로 근처에서 재결합이 더 중요해집니다.

적층 전체에 사용되는 정전기 기준값은 모든 탭에서 확인할 수 있습니다: 전자 친화도 χ = 4.05 eV, 밴드갭 E_g = 1.75 eV, 그리고 상대 유전율 ε_r = 11입니다. 더 넓은 밴드갭은 a-Si:H의 광학 밴드갭을 반영하며 달성 가능한 광전압의 상한을 설정하고, 전자 친화도와 유전율은 밴드 정렬과 전계 침투를 정의합니다. 3절의 p+/i/n+ 도핑 프로파일과 함께 ??, 이러한 파라미터들은 고유층 전반에 걸친 강한 내장 전기장을 형성하며, 이는 박막 p–i–n 소자에서 주요 캐리어 분리 메커니즘입니다.

a-Si:H에서 재결합은 결함 매개 과정이 지배하므로 중심 손실 모델은 국소화된 트랩 상태를 통한 Shockley–Read–Hall (SRH) 재결합입니다. 이는 모든 영역에서 Equilibrium SRH traps 블록을 통해 활성화됩니다 (??, ??, 그리고 ??의 SRH 파라미터 필드 참조). 이 정식화에서 트랩 밀도 N_t, 포획 단면적 σ_n 및 σ_p, 그리고 트랩 에너지 준위 E_t는 캐리어 포획 속도를 설정하며 따라서 유효 수명을 결정합니다: \[ \tau_n = \frac{1}{\sigma_n v_{\mathrm{th}} N_t}, \qquad \tau_p = \frac{1}{\sigma_p v_{\mathrm{th}} N_t}, \] 여기서 \(v_{\mathrm{th}}\)는 열속도입니다. 이에 대응하는 SRH 재결합 속도는 표준적인 형태를 가집니다 \[ R_{\mathrm{SRH}} = \frac{np - n_i^2}{\tau_p (n+n_1) + \tau_n (p+p_1)}, \] 여기서 \(n_1\)과 \(p_1\)은 트랩 에너지 준위 \(E_t\)에 의해 결정됩니다. 대부분의 광생성이 고유 흡수층에서 일어나므로, i 층에서의 SRH 재결합은 이 튜토리얼에서 준-페르미 준위 분리와 따라서 V_oc를 제한하는 주요 메커니즘입니다.

Auger 재결합 파라미터는 p+ 및 n+ 탭에도 존재하며 (?? 및 ?? 참조), 계수는 C_n = C_p = 1×10⁻⁴³ m⁶s⁻¹입니다. Auger 재결합은 다음 형태의 고밀도 3입자 손실 과정을 나타냅니다 \[ R_{\mathrm{Auger}} = C_n\,n^2 p + C_p\,p^2 n, \] 이는 고농도로 도핑된 접촉 영역이나 Suns–Voc 시뮬레이션에서 매우 강한 조명 조건과 같이 캐리어 밀도가 클 때 중요해집니다. 비정질 실리콘의 표준 1-sun 동작에서는 Auger 재결합은 일반적으로 2차적인 효과이며, 전압 손실은 결함 매개 Shockley–Read–Hall(SRH) 재결합이 지배합니다. 그럼에도 Auger 항을 포함하면 물리적으로 일관된 고주입 채널을 제공하고, 소자가 매우 높은 캐리어 밀도로 구동될 때 (예를 들어 수백~수천 suns에서) 비물리적인 캐리어 축적을 방지합니다. Auger 재결합과 그 구현에 대한 자세한 논의는 Auger recombination theory에 제시되어 있습니다.

종합하면, 낮은 캐리어 이동도, 상대적으로 넓은 유효 밴드갭, 그리고 높은 밀도의 국소화된 결함 상태로 인해 a-Si:H는 분명히 재결합 제한 동작 영역에 놓이게 됩니다. 그 결과 JV 곡선과 Suns–Voc 응답은 모두 주로 캐리어 밀도 증가에 따라 유효 SRH 수명이 어떻게 붕괴하는가에 의해 지배됩니다. 따라서 이 튜토리얼은 전압 거동을 결정질 실리콘식의 고이동도 및 수송 제한 동작 가정보다, 준-페르미 준위 분리와 결함 제어 재결합의 관점에서 직접 해석합니다. SRH 재결합 물리와 그것이 반도체 소자 모델링에서 수행하는 역할에 대한 자세한 설명은 이론 섹션 Shockley–Read–Hall recombination theory에 제공되어 있습니다.

명확성과 수치적 투명성을 위해, 이 튜토리얼은 결함 재결합을 고정된 포획 단면적을 갖는 단일 유효 SRH 트랩 준위로 모델링합니다. 이는 현재 튜토리얼을 위한 의도적인 모델링 선택입니다. 즉, a-Si:H 태양전지에 관련된 지배적 재결합 물리를 포착하면서도 파라미터 공간을 간결하게 유지합니다. OghmaNano는 또한 트랩 분포가 전하를 능동적으로 저장하고 바이어스 및 조명 하에서 동적으로 진화하는 더 발전된 처리도 지원합니다. 이러한 비평형 처리는 시간 영역 거동이 고려될 때, 소자 충전이 중요할 때, 또는 과도 응답이나 히스테리시스가 있는 실험 데이터와의 정량적 일치를 요구할 때 필수적입니다. 이러한 확장과 그 물리적 동기는 Why traps are required in disordered semiconductors 및 Non-equilibrium Shockley–Read–Hall recombination에 설명되어 있습니다.

5. 광 생성 프로파일

이 데모는 드리프트–확산 해석기를 위한 깊이 의존 소스 항을 생성하기 위해 1차원 광학 계산을 사용합니다. 실제로 광학 모델은 입사 AM1.5 스펙트럼이 박막 적층 내부로 어떻게 전파되는지, 파장과 깊이의 함수로 얼마나 흡수되는지 계산하고, 그 흡수된 전력을 전자–정공 쌍 생성 속도로 변환합니다. 박막 a-Si:H에서는 흡수층의 두께가 일반적으로 수백 나노미터에 불과하므로, 광학 프로파일은 적층 설계와 밀접하게 결합됩니다. 유용한 흡수의 상당 부분은 조명되는 쪽 근처에서 발생하며, 후면 반사/광 포획 전략은 J_sc에 강하게 영향을 줍니다.

광학 해를 보려면(및 다시 계산하려면) 메인 창의 Optical 리본으로 가서 Transfer matrix를 클릭합니다. 그러면 광학 시뮬레이션 편집기가 열립니다. 파란색 Run optical simulation 버튼(재생 아이콘)을 눌러 광학장을 계산하고 플롯을 업데이트합니다 (?? 참조).

??의 플롯은 파장–깊이 맵입니다. 가로축은 소자 내부의 깊이(y-position으로 표시)이고, 세로축은 파장입니다. 조명된 표면 근처의 밝은 색은 높은 광자 개체수를 나타내고, 깊이에 따라 희미해지는 것은 적층 내에서의 흡수를 나타냅니다.

a-Si:H의 경우 가시광에서의 흡수가 강하므로, 유용한 광자의 큰 부분이 전면으로부터 비교적 짧은 거리 내에서 흡수됩니다. 장파장 한계는 a-Si:H 유효 밴드갭에 의해 제어됩니다. 밴드갭 이하의 광자는 캐리어를 효율적으로 생성하지 못하므로, 광자 맵에는 그 이후 흡수와 생성이 강하게 감소하는 스펙트럼 경계가 나타납니다. 이것이 박막 a-Si:H 소자가 일반적으로 절대적인 J_sc를 더 높은 밴드갭과 향상된 스펙트럼 선택성과 맞바꾸는 이유 중 하나입니다.

광학 편집기는 “광자 밀도”가 “생성 속도”와 같지 않기 때문에 여러 탭을 제공합니다. Photon distribution 보기는 광학장이 어디에 존재하는지 알려줍니다. 전기 모델에 중요한 것은 Photon distribution absorbed 및 Generation rate 탭입니다. 흡수된 광자 맵은 깊이 의존적인 전자–정공 쌍 생성 속도로 변환되며, 이것이 드리프트–확산 방정식에 주입되는 소스 항입니다. 이 a-Si:H 전지의 경우 결과적인 생성 프로파일은 전면에 강하게 치우쳐 있고 얇은 흡수층 내부에 제한되며, 바로 이 영역에서 개방회로 근처의 재결합과 수송이 V_oc에 결정적이 됩니다.

6. 시뮬레이션 실행, JV 곡선, 파라미터 추출

소자 구조, 전기적 파라미터, 광 생성이 정의되면 메인 창에서 직접 시뮬레이션을 실행할 수 있습니다. 해석기를 시작하려면 툴바의 Run simulation 버튼을 클릭합니다. 실행 중에는 해석기 진행 상황과 수렴 정보가 터미널 창에 기록됩니다 (?? 참조).

각 바이어스 지점마다 먼저 상부 접촉에 인가된 전압이 나열되고, 그 뒤에 결과 전류 밀도가 표시됩니다. 조명 하에서는 전류가 처음에 음수(전력 생성)입니다. 인가 전압이 증가함에 따라 전류 크기는 감소하고, 개방회로 전압에서 0을 지나게 됩니다. 이 지점을 넘으면 전류는 양수가 되며, 이는 순방향 바이어스된 다이오드 동작에 해당합니다. 터미널 출력은 또한 각 바이어스 단계에 대한 해석기 잔차와 수렴 시간도 보고합니다. 작은 잔차는 결합된 드리프트–확산 및 Poisson 방정식이 정확하게 풀리고 있음을 나타냅니다.

전기적 특성을 확인하려면 Output 탭을 열고 jv.csv를 더블클릭합니다. 그러면 전류 밀도–전압(JV) 곡선이 표시됩니다 (?? 참조). JV 곡선은 소자 거동에 대한 주요 진단 도구입니다. 0 바이어스에서 단락 전류를 지나야 하고, 개방회로 전압에서 전류 0을 교차해야 합니다.

siminfo.dat를 더블클릭하면 시뮬레이션 정보 창이 열립니다 (?? 참조). 여기에는 충전율, 전력 변환 효율, 최대 전력점, V_oc, J_sc를 포함한 추출된 성능 지표가 보고됩니다. 개방회로에서의 자유 캐리어 밀도와 같은 추가 진단량도 a-Si:H에서의 결함 제한 거동을 해석하는 데 유용합니다.

실용적인 규칙으로서, 수치 지표를 해석하기 전에 항상 JV 곡선을 먼저 확인해야 합니다. JV 곡선이 단락 및 개방회로를 깔끔하게 지나지 않거나, 전류의 부호나 형상이 예상과 다르면, siminfo.dat의 유도 물리량 또한 신뢰할 수 없습니다. JV 곡선을 시각적으로 확인하는 것은 구성 또는 모델링 문제를 진단하는 가장 빠른 방법입니다.

7. Suns–Voc 분석: a-Si:H에서의 재결합 제한 전압

Suns–Voc 측정은 개방회로 전압이 전류가 엄격히 0인 조건에서 조명 강도의 함수로 어떻게 변화하는지를 조사합니다. 단자 전류가 0으로 제한되므로, Suns–Voc는 수송 및 접촉 저항 효과로부터 재결합 물리를 분리합니다. 따라서 결과 곡선은 태양전지에서 전압 손실 메커니즘을 진단하는 가장 직접적인 방법 중 하나입니다.

OghmaNano에서 Suns–Voc는 전용 시뮬레이션 모드로 구현되어 있습니다. 인가 전압을 스윕하는 대신, 해석기는 조명 강도를 스윕하고 각 강도에서 순 전류가 0이 될 때까지 단자 전압을 조정합니다. 이렇게 하면 각 조명 수준에서 개방회로 동작점을 직접 계산할 수 있습니다.

Suns–Voc를 활성화하려면 메인 창의 Simulation type 리본을 열고 (?? 참조) Suns–Voc를 선택합니다. 그런 다음 Run simulation을 클릭합니다. 해석기는 결과를 자동으로 디스크에 기록하며, suns_voc.csv 파일에는 조명 강도의 함수로서의 V_oc가 포함됩니다.

낮은 조명에서 V_oc의 증가는 고유 흡수층 내의 과잉 캐리어 밀도 축적에 의해 지배됩니다. 비축퇴 한계에서 준-페르미 준위 분리는 \[ qV_{\mathrm{oc}} = E_{Fn} - E_{Fp} = kT \ln\!\left(\frac{np}{n_i^2}\right), \] 이므로 V_oc는 캐리어 밀도에 대해서만 로그적으로 증가합니다.

비정질 실리콘에서는 이러한 낮은 강도 영역이 결함 점유 효과에 의해 종종 수정됩니다. 매우 낮은 생성률에서는 광생성된 캐리어의 상당 부분이 자유 캐리어 개체수에 기여하기 전에 깊은 상태 및 테일 상태에 포획됩니다. 이는 ??에서 관찰되듯이 Suns–Voc 곡선에서 약하거나 다소 평탄한 초기 응답을 만들 수 있습니다. 이러한 거동은 수치적 인공산물이 아니라 트랩 충전의 물리적 결과입니다.

고주입 영역을 탐구하려면 조명 범위를 확장합니다. Editors 리본에서 Suns–Voc editor를 열고 stop intensity를 1000 suns로 증가시킵니다. 그런 다음 시뮬레이션을 다시 실행합니다.

시뮬레이션을 다시 실행한 후 suns_voc.csv를 열면 전체 조명 범위의 거동이 드러납니다 (?? 참조). 높은 조명에서 V_oc는 계속 증가하지만, 기울기는 점차 감소하며 결국 포화에 접근합니다.

전압 포화의 기원은 Suns–Voc 곡선을 저장 전하 및 유효 수명과 함께 볼 때 분명해집니다. 조명이 증가함에 따라 총 과잉 캐리어 밀도는 증가하지만 (?? 참조), 동시에 유효 재결합 시간은 감소합니다 (?? 참조).

a-Si:H에서 재결합은 결함 매개 SRH형 과정이 지배합니다. 캐리어 밀도가 증가함에 따라 트랩 보조 재결합이 가속되고 유효 수명 \[ \tau_{\mathrm{eff}} = \frac{\Delta n}{R} \] 은 감소합니다. 이 지점을 넘어서면 추가 광생성은 준-페르미 준위 분리를 증가시키기보다는 재결합을 주로 증가시킵니다.

생성과 재결합 사이의 이러한 균형은 V_oc에 대한 본질적인 상한을 부여합니다. 최대 전압은 전자 및 정공 준-페르미 준위가 각각의 밴드 에지에 동시에 접근하는 것을 재결합이 막기 때문에 유효 a-Si:H 밴드갭보다 낮게 유지됩니다. 따라서 E_g/q − V_oc의 차이는 이 소자에서의 근본적인 재결합 유도 전압 손실을 나타냅니다.