실리콘 PN 접합 다이오드 (1D) — Drift–Diffusion (암 I–V, SRH 재결합)
1. 소개
실리콘 PN 접합 다이오드는 전형적인 반도체 소자입니다. 이는 개별 소자로 명시적으로 사용되기도 하고, 전력 정류기부터 로직 및 아날로그 회로에 이르기까지 거의 모든 집적 전자 시스템 전반에 암묵적으로 존재합니다. 대표적인 응용 맥락은 ??에 나타나 있으며, 여기서 PN 접합은 독립형 소자로 사용되기보다 집적 회로 내부에 내장되어 있습니다.
이 튜토리얼에서 시뮬레이션하는 소자는 단순하지만, 기본적인 구성 블록으로 이해해야 합니다. 동일한 접합 물리는 다이오드 접속 소자, 트랜지스터 접합, 그리고 집적 실리콘 기술 내부의 절연 구조를 지배합니다. 여기서 사용하는 층상 도핑 구조는 ??에 개략적으로 나타나 있습니다.
이 튜토리얼에서는 OghmaNano의 결합된 drift–diffusion + Poisson 솔버를 사용하여 1차원 실리콘 PN 접합 다이오드를 시뮬레이션합니다. 이상적인 Shockley 방정식에만 의존하는 대신, 이 접근법은 내장 전기장, 공핍 영역, 그리고 캐리어 밀도 및 전류의 공간 분포를 해석합니다.
특히 Shockley–Read–Hall (SRH) 재결합을 포함한 재결합 물리의 명시적 처리는 도핑과 수명의 변화를 턴온 거동과 이상성의 변화에 직접 연결할 수 있게 합니다. 암 I–V 곡선을 생성하고, 바이어스 하에서 밴드 에지와 준페르미 준위를 확인한 다음, 다이오드 응답에서 재결합 제한 영역을 식별하기 위해 수명 스윕을 수행할 것입니다.
2. 새 시뮬레이션 만들기
시작하려면 OghmaNano의 메인 창에서 새 시뮬레이션을 생성합니다. 툴바에서 New simulation 버튼을 클릭합니다. 그러면 시뮬레이션 유형 선택 대화상자가 열립니다 ( ?? 참조).
시뮬레이션 유형 대화상자에서 Si demos를 더블 클릭한 다음, 실리콘 접합/다이오드 예제를 선택합니다 ( ?? 참조). OghmaNano는 미리 정의된 실리콘 접합 구조를 불러오며, 이를 PN 다이오드로 취급할 것입니다.
불러온 소자 구조는 메인 시뮬레이션 창에 표시됩니다 ( ?? 참조). 이 튜토리얼에서 푸는 전기적 문제는 1차원이지만, 3D 뷰는 수직 층 구조와 캐리어 수송 및 재결합에 참여하는 영역을 명확하게 시각화합니다.
다이오드는 수직으로 적층된 실리콘 층들의 연속으로 구현되며, 고농도 도핑된 p+ 영역, 더 낮게 도핑된 p 영역, 낮게 도핑된 n 영역, 그리고 고농도 도핑된 n+ 영역으로 구성됩니다. 이 구조는 Layer editor에 명시적으로 나열되어 있으며 ( ?? 참조), 각 층에는 두께, 재료, 전기적 역할이 지정되어 있습니다.
중심의 p 및 n 층이 활성 PN 접합을 형성합니다. 평형 상태에서는 이 계면을 가로질러 공핍 영역이 형성되며, 이는 캐리어 분리와 수송을 제어하는 내장 전기장을 발생시킵니다. 얇고 고농도 도핑된 p+ 및 n+ 층은 저항이 낮은 접촉 영역으로 작용하여, 인가된 바이어스가 접촉이 아니라 주로 접합 전체에 걸리도록 합니다.
다음 절들에서 이 구조는 1차원 소자로 취급됩니다: 모든 변화는 성장 방향을 따라 해석되며, 측면 변화는 무시됩니다. 이러한 단순화에도 불구하고, 이 모델은 실제 전자 소자에 내장된 실리콘 PN 접합 다이오드의 암 I–V 거동을 지배하는 본질적인 정전기, 캐리어 수송, 그리고 재결합 물리를 포착합니다.
3. 도핑 프로파일 확인
도핑 프로파일은 실리콘 PN 접합을 정의하며, 따라서 다이오드의 기본 정전기적 특성을 결정합니다. 이는 접합 위치, 내장 전위, 공핍 폭, 그리고 평형 상태 및 바이어스 하에서 형성되는 내부 전기장을 결정합니다.
도핑 구성을 확인하려면 Electrical 리본에서 Doping / Ions 편집기를 엽니다 ( ?? 참조). 편집기는 깊이의 함수로서 이온화된 도너와 억셉터의 공간 분포를 표시합니다 ( ?? 참조).
이 튜토리얼에서 다이오드는 전통적인 p+/p/n/n+ 실리콘 도핑 프로파일을 사용하여 구성됩니다. 중심의 p 및 n 영역은 중간 정도로 도핑되어 있으며 활성 PN 접합을 형성하고, 여기서 공핍 영역과 내장 전기장이 발달합니다.
얇고 고농도 도핑된 p+ 및 n+ 층은 저항이 낮은 접촉 영역으로 작용합니다. 이들의 역할은 캐리어의 전기적 주입 및 추출을 원활하게 제공하는 것이며, 인가된 전압의 대부분이 접촉이 아니라 접합 자체에 걸리도록 하는 것입니다.
이 튜토리얼의 목적상 핵심 점검 사항은 소자에 하나의 주된 억셉터 도핑 영역과 하나의 주된 도너 도핑 영역이 존재하고, 그 사이에 뚜렷한 전이가 있는지 여부입니다. 도핑 밀도의 정확한 수치 값은 주로 공핍 폭과 내장 전기장의 세기에 영향을 미치며, 이는 이후 절에서 다이오드의 암 I–V 특성을 통해 간접적으로 살펴볼 것입니다.
4. 전기적 매개변수와 재결합 메커니즘 확인
전기적 재료 매개변수는 영역별로 정의되며 다이오드 전반의 캐리어 수송, 재결합, 그리고 정전기적 특성을 제어합니다. 메인 창에서 Device structure → Electrical parameters를 통해 전기적 매개변수 편집기를 엽니다. 소자 구조의 각 층은 고유한 매개변수 탭을 가집니다. 이 튜토리얼에서는 네 개 영역 모두에 동일한 실리콘 재료 모델을 사용하지만, 이를 다르게 해석합니다: p+ 및 n+는 저항이 낮은 접촉 영역으로 작용하고, p 및 n은 활성 접합을 형성합니다.
그림 ??– ??는 각 영역 (p+, p, n, n+)에 대한 전기적 매개변수 편집기를 보여줍니다. 이 데모에서는 전자 및 정공 이동도가 각각 약 0.135 및 0.048 m2V−1s−1로 설정되어 있습니다 (예를 들어, ?? 참조). 이러한 값은 결정질 실리콘에 전형적입니다: 비정질 실리콘이나 다수의 용액 공정 반도체와 같은 무질서하거나 결함 제한적인 재료보다 훨씬 높지만, GaAs와 같은 고이동도 III–V 재료보다는 낮습니다. 그 결과, 이 튜토리얼에서의 다이오드 거동은 주로 벌크 수송 제한이 아니라 접합 정전기와 재결합에 의해 제어됩니다.
여기서 사용된 유효 상태 밀도도 편집기에서 확인할 수 있습니다 (예: ??): 유효 자유전자 상태 밀도는 약 2.8×1025 m−3, 유효 자유정공 상태 밀도는 1.04×1025 m−3로 설정되어 있습니다. 이러한 매개변수는 주어진 밴드 구조에 대한 캐리어 통계와 평형 캐리어 농도를 정하며, \(n\)과 \(p\)를 통해 재결합 및 주입 거동에 암묵적으로 작용합니다.
Shockley–Read–Hall (SRH) 재결합
SRH 재결합은 밴드갭 내 전자 상태를 통한 결함 매개 재결합을 포착합니다. OghmaNano에서는 편집기에 표시된 평형 SRH 트랩 매개변수로 이를 제어합니다 ( ?? 참조): 트랩 에너지 \(E_t\) (중간 밴드갭 기준), 트랩 밀도 \(N_t\), 그리고 전자 및 정공 포획 단면적 \(\sigma_n\) 및 \(\sigma_p\)입니다. 여기에 표시된 설정에서는 트랩 에너지가 중간 밴드갭에 가깝고 (\(E_t \approx 0\)), 트랩 밀도는 \(N_t \approx 10^{21}\,\mathrm{m^{-3}}\), 포획 단면적은 \(\sigma_n \approx \sigma_p \approx 10^{-21}\,\mathrm{m^2}\)입니다.
이러한 미시적 매개변수는 다음과 같이 SRH 수명을 정의합니다
\[ \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}}\)는 열적 캐리어 속도입니다. 따라서 트랩 밀도나 포획 단면적을 증가시키면 캐리어 수명이 감소하고 재결합이 강화됩니다.
drift–diffusion 형태에서, 결과적인 SRH 재결합 속도는 다음과 같습니다
\[ R_{\mathrm{SRH}} = \frac{np - n_i^2} {\tau_p (n + n_1) + \tau_n (p + p_1)} . \]이후 절에서 “수명 스윕”을 수행할 때, 실제로는 트랩 보조 재결합 세기를 정의하는 \(N_t\), \(\sigma_n\), 그리고 \(\sigma_p\)를 수정하게 됩니다. 가장 큰 영향은 공핍 영역 내부 및 그 근처의 재결합이 주입된 캐리어 개체군을 제어할 때 나타나며, 이는 재결합 제한 순방향 바이어스 거동으로 이어집니다.
Auger 재결합
Auger 재결합은 높은 캐리어 밀도에서의 손실 메커니즘이며, 따라서 고농도 도핑된 p+ 및 n+ 접촉 영역에서 가장 중요합니다. 매개변수 편집기에는 Auger 계수 \(C_n\) 및 \(C_p\)가 지정되어 있습니다 ( ?? 및 ?? 참조). 이 실리콘 데모에서 계수의 크기는 \(C_n \approx 2.8\times10^{-31}\,\mathrm{m^6\,s^{-1}}\) 및 \(C_p \approx 9.9\times10^{-32}\,\mathrm{m^6\,s^{-1}}\) 정도이며, 이는 결정질 실리콘에 전형적인 값입니다.
Auger 재결합 속도는 다음과 같이 주어집니다
\[ R_{\mathrm{Auger}} = C_n n^2 p + C_p p^2 n . \]실제로 Auger 재결합은 고주입 조건에서 캐리어 축적을 제한하고, 접촉 영역이 캐리어에 대한 효율적인 싱크로 동작하도록 보장하지만, 중간 정도로 도핑된 접합 자체의 재결합 물리를 지배하지는 않습니다.
정전기와 밴드 매개변수
마지막으로, 실리콘을 정의하는 데 사용되는 밴드 구조 및 정전기적 매개변수는 각 영역 탭에서 확인할 수 있습니다 (예: ??): 전자 친화도는 약 \(\chi \approx 4.05\,\mathrm{eV}\), 밴드갭은 \(E_g \approx 1.12\,\mathrm{eV}\), 상대 유전율은 \(\varepsilon_r \approx 11.7\)로 설정되어 있습니다.
5. 시뮬레이션 실행, 암 I–V 곡선, 그리고 매개변수 추출
소자 구조, 도핑 프로파일, 전기적 매개변수가 정의되면 다이오드 시뮬레이션을 메인 창에서 직접 실행할 수 있습니다. 솔버를 시작하려면 Run simulation을 클릭합니다. 실행 중에는 각 바이어스 점에 대한 수렴 정보가 터미널에 기록되므로, 솔버의 안정성과 진행 상황을 모니터링할 수 있습니다 ( ?? 참조).
jv.csv입니다.
다이오드 특성을 확인하려면 Output 탭을 열고
jv.csv를 더블 클릭합니다
( ?? 참조).
올바르게 설정된 실리콘 다이오드의 경우 I–V 곡선은 매끄럽고 단조적이어야 합니다.
역방향 바이어스에서는 전류가 작고 전압 의존성이 약하며,
이는 재결합 제한 포화를 반영합니다.
순방향 바이어스에서는 인가 전압과 함께 전류가 급격히 증가하며,
이는 PN 접합을 가로지르는 캐리어 주입에 해당합니다.
순방향 바이어스 영역의 형상에는 유용한 물리적 정보가 포함되어 있습니다. 반로그 플롯에서 지수 영역의 기울기를 사용하면 이상성 인자를 추출할 수 있으며, 이는 전류가 확산 제한 수송(\(n \approx 1\))에 지배되는지 혹은 재결합 제한 과정 (\(n \approx 2\))에 지배되는지를 나타냅니다. 이 영역의 외삽 절편은 역포화 전류의 추정값을 제공하며, 이는 4절에서 논의한 SRH 및 Auger 재결합 매개변수와 직접 연결됩니다.
실용적인 규칙으로, 파생량을 해석하기 전에 항상 I–V 곡선을 확인하십시오. 불연속성, 예상치 못한 부호 규약, 또는 비물리적인 전류 도약은 대개 경계 조건, 바이어스 스텝, 재결합 설정, 또는 솔버 수렴 문제를 나타냅니다. 이와 같은 단순한 실리콘 PN 다이오드의 경우 암 I–V 곡선은 물리적으로 직관적이고 해석하기 쉬워야 합니다.
6. 시뮬레이션 스냅샷 확인: 밴드, 재결합, 그리고 전류 흐름
I–V 스윕 동안 OghmaNano는 drift–diffusion 방정식의 내부 해를 각 바이어스 점에서 snapshots 디렉터리에 저장합니다. 이 파일들은 솔버가 다이오드 내부에서 무엇을 예측하는지 보여줍니다: 밴드 굽힘, 준페르미 준위 분리, 재결합 활동, 그리고 전류 수송입니다. 이러한 양을 확인하는 것은 특정 I–V 특성이 왜 나타나는지 이해하는 데 필수적입니다.
이 절에서는 세 개의 대표적인 바이어스 점을 살펴봅니다: 평형에 가까운 역방향 바이어스 (−0.1 V), 턴온 근처의 중간 순방향 바이어스 (≈0.45 V), 그리고 높은 순방향 바이어스 (0.8 V)입니다. 이 스냅샷들은 함께 평형 상태에서 주입 제한 수송을 거쳐 고주입 동작으로의 전이를 보여줍니다.
6.1 밴드 에지와 준페르미 준위
밴드 다이어그램을 재현하려면 스냅샷 뷰어를 열고
Ec.csv, Ev.csv, Fn.csv, 그리고 Fp.csv 파일을 추가합니다.
이들은 각각 전도대 에지, 가전자대 에지,
전자 준페르미 준위, 그리고 정공 준페르미 준위에 해당합니다.
−0.1 V에서 (그림 ??), 다이오드는 평형에 가깝습니다. 밴드 굽힘은 도핑 프로파일이 부과하는 내장 전위를 반영하며, 준페르미 준위는 거의 평탄하고 서로 일치하여, 순 전류 흐름이 무시할 정도로 작음을 나타냅니다. 공핍 영역은 접합에서 강한 밴드 곡률을 갖는 영역으로 명확하게 보입니다. ≈0.45 V에서 (그림 ??), 순방향 바이어스가 접합 장벽을 낮춥니다. 전자 및 정공 준페르미 준위는 공핍 영역을 가로질러 분리되며, 이는 캐리어 주입의 내부적 신호입니다. 이러한 준페르미 준위 분리는 I–V 곡선에서 관찰되는 전류의 지수적 증가에 직접적으로 책임이 있습니다. 0.8 V에서 (그림 ??), 접합은 깊은 순방향 바이어스 상태에 있습니다. 장벽은 크게 억제되고, 준페르미 준위는 넓게 분리되며, 소자는 구조의 상당 부분에서 캐리어 밀도가 큰 고주입 영역에서 동작합니다.
6.2 Shockley–Read–Hall 재결합
재결합을 확인하려면 다이오드 내부의 공간적으로 분해된
Shockley–Read–Hall 재결합 속도를 보여주는 R_ss_srh.csv를 플롯합니다.
아래의 세 플롯은 밴드 다이어그램 분석에서 사용한 것과 동일한 바이어스 점,
즉 −0.1 V, ≈0.45 V, 그리고 0.8 V에 해당합니다.
여기서 주목할 핵심은 재결합 속도의 절대 크기가 아니라,
인가 바이어스에 따라 그 공간적 국소화가 어떻게 변하는지입니다.
−0.1 V에서 (그림 ??), 다이오드는 평형에 가깝습니다. 전자는 n형 쪽에서 우세하고 정공은 p형 쪽에서 우세하므로, 두 캐리어 유형이 동시에 존재하는 접합 주변의 좁은 영역에서만 유의미한 재결합이 일어날 수 있습니다. 그 결과, SRH 재결합 속도는 소자의 중심부에 강하게 국소화되며, 공핍 영역과 일치합니다. ≈0.45 V에서 (그림 ??), 순방향 바이어스가 접합을 가로질러 캐리어를 주입하고 전자 및 정공 밀도의 국소 곱을 증가시킵니다. 재결합 피크의 크기는 크게 증가하지만, 여전히 소자의 중앙 영역에 공간적으로 제한됩니다. 이는 이 바이어스 범위에서 SRH 재결합이 여전히 주로 공핍 영역 내부 및 그 근처의 캐리어 중첩에 의해 제어되는 접합 중심 과정임을 나타냅니다. 0.8 V에서 (그림 ??), 거동은 질적으로 변합니다. 캐리어 주입이 충분히 강해져 전자와 정공이 다이오드의 상당 부분 전반에 걸쳐 높은 농도로 존재하게 됩니다. SRH 재결합 속도는 더 이상 접합에 제한되지 않고, 소자의 넓은 영역으로 퍼집니다. 이러한 공간적 확장은 재결합이 더 이상 좁은 중심 영역에만 제한되지 않는 고주입 조건의 시작을 나타냅니다.
날카롭게 국소화된 재결합 피크에서 공간적으로 확장된 재결합 프로파일로의 진행은 다이오드가 평형 상태에서 접합 제한 동작을 거쳐 구조 전반에서 재결합이 발생하는 영역으로 전이하는 방식을 명확한 내부 그림으로 제공합니다. 이 변화는 밴드 다이어그램에서 보이는 변화와 대응하며 순방향 I–V 특성의 기울기 변화를 뒷받침합니다.
6.3 전자 및 정공 전류 밀도
마지막으로 Jn.csv와 Jp.csv를 플롯하여 캐리어 전류를 확인합니다.
이 파일들은 각각 공간적으로 분해된 전자 및 정공 전류 밀도를 보여줍니다.
이 플롯들은 다양한 바이어스 조건에서 전하가 다이오드를 통해 어떻게 수송되는지,
그리고 소자가 평형 상태에서 정상 상태 순방향 전도로 어떻게 전이하는지를
직접적으로 보여줍니다.
−0.1 V에서 (그림 ??), 다이오드는 평형에 가깝고 실제 물리적 전류는 극히 작습니다. 전자 및 정공 플럭스는 소자 전체에서 거의 균형을 이루므로, 순전류는 거의 동일한 두 양의 차이로부터 발생합니다. 이 영역에서는 수치 문제가 본질적으로 악조건이며, 전류 프로파일에 작은 진동이나 겉보기 잡음이 나타나는 것이 예상됩니다. 이러한 특징은 수치적 기원이며 실제 캐리어 수송에 해당하지 않습니다. ≈0.45 V에서 (그림 ??), 순방향 바이어스는 접합을 가로지르는 캐리어 주입을 유도합니다. 전자 전류는 n측에서 우세하고 정공 전류는 p측에서 우세하지만, 두 전류 모두 소자를 통해 연속적이며, 이는 정상 상태 전하 보존을 반영합니다. 전류 밀도는 평형 근처 경우와 비교하여 급격히 증가하지만, 접합 근처에서는 여전히 공간 구조를 가지며, 이는 재결합에 의해 제어되는 주입 제한 수송과 일치합니다. 0.8 V에서 (그림 ??), 다이오드는 깊은 순방향 바이어스에서 동작합니다. 캐리어 밀도는 구조 전반에 걸쳐 높고, 전자 및 정공 전류는 모두 커지고 매끄러우며 준중성 영역 전반에서 거의 균일해집니다. 이 영역에서 소자는 강하게 전도하는 요소처럼 동작하며, 전류는 장벽 주입보다 주로 수송과 재결합에 의해 제한됩니다.
종합하면, 이러한 전류 밀도 플롯은 다이오드 동작에 대한 일관된 내부 그림을 제공합니다: 평형 상태에서 전자 및 정공 플럭스의 거의 완전한 상쇄로부터, 주입 제한 순방향 전도를 거쳐, 큰 순방향 바이어스에서의 고전류 정상 상태 수송까지입니다.
6.4 고농도 도핑 영역에서의 Auger 재결합
위치의 함수로서 국소 Auger 재결합 속도를 보고하는
R_auger.csv를 플롯하여 Auger 재결합을 확인할 수 있습니다.
전자와 정공이 비슷한 밀도로 공존하는 곳에서 가장 강한 SRH 재결합과 달리,
Auger 재결합은 캐리어 밀도에 강하게 비례하므로
고농도 도핑된 영역에서 우세해집니다.
매우 낮은 바이어스에서 (그림 ??), Auger 재결합은 이미 p+ 및 n+ 접촉층에서 가시적입니다. 이는 평형 근처에서도 이러한 영역이 의도적으로 고농도 도핑되어 있어, 국소적으로 삼입자(Auger) 과정이 우세할 만큼 캐리어 농도가 높기 때문입니다.
다이오드가 순방향 바이어스로 구동되면 (그림 ??), Auger 재결합 속도는 크기가 급격히 증가하지만 여전히 공간적으로 접촉 영역에 제한됩니다. 이러한 국소화는 Auger 재결합의 강한 밀도 의존성을 반영합니다: 캐리어가 접합을 가로질러 주입되더라도, 가장 높은 캐리어 밀도는 여전히 축퇴된 접촉층에 존재합니다.
높은 순방향 바이어스에서 (그림 ??), Auger 재결합은 접촉에서 극도로 커지며, 크기 면에서 SRH 재결합을 압도합니다. 이러한 거동은 예상 가능하며 물리적으로 타당합니다. 여기서 Auger 재결합의 역할은 접합 전류를 직접 제한하는 것이 아니라, 고주입 조건에서 캐리어 밀도가 무한정 증가하는 상황을 막아 비물리적인 캐리어 축적을 방지하는 것입니다.
중요한 점은, Auger 재결합 속도가 수치적으로 SRH 속도보다 훨씬 크더라도, 다이오드의 이상성이나 턴온 거동을 지배하지는 않는다는 것입니다. 이러한 특성은 여전히 주로 공핍 영역 내부 및 그 근처의 재결합에 의해 제어됩니다. 대신 Auger 재결합은 접촉층에서의 고밀도 안정화 메커니즘으로 작용하여, 큰 전류에서 현실적인 거동을 보장합니다.