초고색순도(BT.2020 색영역 95% 이상)와 저비용 용액 공정으로 OLED를 잇는 차세대 디스플레이의 강력한 후보로 떠오른 페로브스카이트 발광다이오드(Perovskite LED, PeLED)는 얇은 발광층 내에서 전자와 정공의 거동에 의해 효율이 지배됩니다. 결정 결함에 의한 트랩 보조 비복사 재결합(SRH, $A$), 쌍분자 복사 재결합(Radiative, $B$), 그리고 고전류 영역에서의 오제 재결합(Auger, $C$) 간의 경쟁을 나타내는 ABC 재결합 속도론 모델과 외부양자효율(EQE) 드룹(Efficiency Droop) 현상을 1,000Hz 주기로 0-Alloc C# 수치해석하는 광전자 반도체 시뮬레이터를 공개합니다.

1. 2022 개정 교육과정 교과 연계 분석

  • 물리학 (반도체 에너지 밴드갭, 전도대-가전자대 전이 및 엑시톤 발광): 직접 천이형 밴드갭에서의 전자-정공 재결합 광자 방출 $h\nu = E_g$, 쿨롱 인력에 의한 엑시톤 속박 에너지 $E_b$와 실온 열에너지($k_B T \approx 26\,\text{meV}$) 비교.
  • 화학/신소재 (유-무기 할라이드 페로브스카이트 결정 구조 및 결함 제어): 3차원 $CH_3NH_3PbBr_3$와 준2D Ruddlesden-Popper 상($n=1, 2, 3$) 나노플레이틀릿 에너지 깔때기(Energy Funneling) 현상 및 트랩 밀도 저감.
  • 수학 (비선형 캐리어 밀도 3차 미분방정식 및 효율 피크 해석): 캐리어 생성율 $G$와 소멸율 $R(n) = An + Bn^2 + Cn^3$의 정상 상태 연립해석, 내부 양자효율 $\text{IQE} = \frac{Bn^2}{An + Bn^2 + Cn^3}$의 극대값 $n_{opt} = \sqrt{A/C}$ 증명.
  • 정보 (무할당 실시간 1,000Hz 주입 전류 밀도별 EQE 및 휘도 연산 코어): 디스플레이 구동 전압-전류 곡선을 1ms마다 샘플링하여 가비지 컬렉터 없이 광자 방출 플럭스와 루미넌스를 0-Alloc 계산하는 C# 온디바이스 엔진.

2. PeLED ABC 재결합 속도론 및 EQE 수식 모델

활성층 내 주입 전류 밀도 $J$, 막 두께 $d$, 기본 전하량 $e$에 대해 캐리어 연속 방정식의 정상 상태는 다음과 같이 주어집니다:

$$\frac{J}{e \cdot d} = A n + B n^2 + C n^3$$

여기서 $A$는 쇼클리-리드-홀(SRH) 단분자 트랩 계수 ($\text{s}^{-1}$), $B$는 복사 재결합 쌍분자 계수 ($\text{cm}^3/\text{s}$), $C$는 3체 오제(Auger) 재결합 계수 ($\text{cm}^6/\text{s}$)입니다. 내부 양자효율(IQE)은 순수 복사 재결합의 분율입니다:

$$\text{IQE}(n) = \frac{B n^2}{A n + B n^2 + C n^3}$$

최대 IQE는 $n_{opt} = \sqrt{\frac{A}{C}}$에서 발생하며, 외부 양자효율(EQE)은 광추출 효율 $\eta_{out}$과의 곱으로 결정됩니다:

$$\text{EQE} = \eta_{out} \cdot \text{IQE} = \eta_{out} \frac{B}{\frac{A}{n} + B + C n}$$

3. 1,000Hz 0-Alloc C# PeLED 재결합 및 EQE 해석 코어

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Education.StemOptoelectronics.PerovskiteLed
{
    public readonly record struct PeledConfig(
        float ActiveLayerThicknessNm,// 활성층 두께 (nm, 예: 45.0 nm)
        float TrapCoeffA,            // A 계수 (1/s, 예: 1.5e6)
        float RadiativeCoeffB,       // B 계수 (cm^3/s, 예: 2.2e-10)
        float AugerCoeffC,           // C 계수 (cm^6/s, 예: 3.5e-28)
        float OutcouplingEfficiency  // 광추출 효율 (예: 0.22)
    );

    public ref struct PeledTelemetry
    {
        public float CurrentDensityMAperCm2;
        public float CarrierDensityCm3;
        public float InternalQuantumEfficiencyPercent;
        public float ExternalQuantumEfficiencyPercent;
        public float RadiativeFluxPhotonsPerS;
        public float LuminanceNit;
    }

    public static class PeledKineticsEngineCore
    {
        private const float ElementaryCharge = 1.60217663e-19f;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public static void CalculateKinetics(
            in PeledConfig config,
            float currentDensityMAperCm2,
            ref PeledTelemetry telemetry)
        {
            telemetry.CurrentDensityMAperCm2 = currentDensityMAperCm2;

            // 1. 주입 생성율 G (cm^-3 * s^-1) (0-Alloc)
            // J in A/cm^2 = mA/cm^2 * 1e-3
            float jAcm2 = currentDensityMAperCm2 * 1e-3f;
            float dCm = config.ActiveLayerThicknessNm * 1e-7f;
            float gGen = jAcm2 / (ElementaryCharge * dCm);

            // 2. 캐리어 밀도 n 수치 근사 (정상상태 G = An + Bn^2 + Cn^3)
            // 저전류: G ~ An -> n0 = G/A, 고전류: G ~ Bn^2 -> n0 = sqrt(G/B)
            float nEst = MathF.Sqrt(gGen / config.RadiativeCoeffB);
            
            // 뉴턴-랩슨 2회 반복으로 0-Alloc 고정밀 수렴
            for (int i = 0; i < 2; i++)
            {
                float n = nEst;
                float f = (config.TrapCoeffA * n) + (config.RadiativeCoeffB * n * n) + (config.AugerCoeffC * n * n * n) - gGen;
                float fPrime = config.TrapCoeffA + (2.0f * config.RadiativeCoeffB * n) + (3.0f * config.AugerCoeffC * n * n);
                nEst = MathF.Max(1e12f, n - (f / MathF.Max(1e-10f, fPrime)));
            }
            telemetry.CarrierDensityCm3 = nEst;

            // 3. 재결합 속도 성분 계산
            float rTrap = config.TrapCoeffA * nEst;
            float rRad = config.RadiativeCoeffB * nEst * nEst;
            float rAuger = config.AugerCoeffC * nEst * nEst * nEst;
            float rTotal = rTrap + rRad + rAuger;

            // 4. 내부 및 외부 양자효율 (IQE, EQE)
            float iqe = (rRad / MathF.Max(1e-15f, rTotal)) * 100.0f;
            telemetry.InternalQuantumEfficiencyPercent = iqe;

            float eqe = iqe * config.OutcouplingEfficiency;
            telemetry.ExternalQuantumEfficiencyPercent = eqe;

            // 5. 광자 방출 플럭스 및 휘도 (nit = cd/m^2) 근사
            float photonRate = rRad * dCm; // photons / (cm^2 * s)
            telemetry.RadiativeFluxPhotonsPerS = photonRate;

            // 525nm 녹색 페로브스카이트 시감도 계수 환산
            float luminance = eqe * currentDensityMAperCm2 * 18.5f;
            telemetry.LuminanceNit = MathF.Max(0.0f, luminance);
        }
    }
}

4. 교육 및 연구 현장 활용 제언

본 시뮬레이터는 한국나노기술원(KANC) 광소자 연구단 및 방통대 우호성 교수 연구실의 차세대 디스플레이 반도체 교원 연수에서 페로브스카이트 소재의 결함 밀도와 전류 주입에 따른 효율 드룹(Droop) 메커니즘을 규명하는 실험 교구로 활용될 수 있습니다. 교사와 학생들은 트랩 계수 $A$와 오제 계수 $C$에 따른 피크 EQE 전류 밀도 위치 변화를 1,000Hz 실시간으로 시뮬레이션하며 차세대 발광소자의 양자 효율 극대화 설계를 완벽히 학습할 수 있습니다.