화석연료를 대체하여 태양광 에너지만으로 물을 직접 분해해 청정 그린 수소를 생산하는 인공광합성 기술의 핵심은 광전기화학(Photoelectrochemical, PEC) 소자입니다. 가시광 흡수 밴드갭($E_g \approx 2.4\,\text{eV}$)을 가진 비스무트 바나데이트($\text{BiVO}_4$) 광양극(Photoanode)과 전해질 계면에서 발생하는 쇼트키 접합 공간 전하 영역(Depletion Width)의 밴드 굽힘(Band Bending)과 전하 분리 효율을 1,000Hz 주기로 0-Alloc C# 수치해석하는 고성능 솔라-투-하이드로겐(STH) 수전해 시뮬레이터를 공개합니다.

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

  • 화학 (산화-환원 반응, 광화학 반응 및 물 분해 열역학): 광양극 물 산화 반쪽 반응 $2\text{H}_2\text{O} + 4h^+ \rightarrow \text{O}_2 + 4\text{H}^+$ ($E^0 = 1.23\,\text{V vs RHE}$), 전자-정공 쌍의 표면 재결합 억제 및 보조 촉매(Co-catalyst) 반응 속도론.
  • 물리학 (광전효과, 반도체 에너지 밴드 구조 및 밴드 굽힘): 페르미 준위 정렬에 따른 반도체-액체 접합(SCLJ) 공핍층 폭 $W = \sqrt{\frac{2\varepsilon_s (V - V_{fb})}{q N_D}}$, 내장 전계(Built-in field)에 의한 전자-정공 표류 분리.
  • 수학 (가르트너 광전류 방정식 및 STH 변환 효율 수치 적분): 단색광 흡수 계수 $\alpha(\lambda)$와 정공 확산 길이 $L_p$를 포함하는 가르트너(Gärtner) 식 $J_{ph} = q \Phi_0 \left( 1 - \frac{e^{-\alpha W}}{1 + \alpha L_p} \right)$ 및 태양 스펙트럼(AM1.5G) 파장 적분.
  • 정보 (무할당 실시간 1,000Hz STH 수소 생산 수율 및 최대 효율 전위 추적 코어): 1ms 제어 루프 내에서 GC 오버헤드 없이 인가 바이어스 전압 대비 수소 발생량과 STH 변환 효율을 0-Alloc으로 최적화하는 C# 온디바이스 엔진.

2. PEC 광양극 밴드 굽힘 및 광전류 수식 모델

인가 전위 $V$, 평탄대 전위(Flat-band potential) $V_{fb}$, 도너 도핑 농도 $N_D$에 대해 공간 전하 공핍층 폭 $W$와 가르트너 광전류 밀도 $J_{ph}$는 다음과 같이 유도됩니다:

$$W = \sqrt{\frac{2 \varepsilon_0 \varepsilon_r (V - V_{fb})}{q N_D}}$$ $$J_{ph} = q I_0 \left( 1 - \frac{\exp(-\alpha W)}{1 + \alpha L_p} \right) \cdot \eta_{sep} \cdot \eta_{inj}$$

여기서 $\eta_{sep}$는 전하 분리 효율, $\eta_{inj}$는 표면 정공 주입 효율이며, 전체 솔라-투-하이드로겐(STH) 효율은 패러데이 효율 $\eta_F \approx 1.0$에 대해 다음과 같이 정의됩니다:

$$\text{STH} = \frac{J_{ph} \times 1.23\,\text{V} \times \eta_F}{P_{sun}} \times 100\%$$

3. 1,000Hz 0-Alloc C# 인공광합성 광수전해 코어

가비지 컬렉터 중단을 배제하고 ref struct 구조체를 통해 1ms 타임스텝마다 광전류 밀도와 시간당 수소 생성량을 초고속 연산합니다. 

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.ArtificialPhotosynthesis;

/// 
/// PEC 인공광합성 셀 상태 구조체 (100% Stack Allocation)
/// 
public readonly ref struct PecCellState
{
    public readonly double BiasVoltageVolt;        // 인가 전위 (V vs RHE)
    public readonly double DepletionWidthNm;       // 공핍층 폭 W (nm)
    public readonly double PhotocurrentDensityMa;  // 광전류 밀도 (mA/cm^2)
    public readonly double SthEfficiencyPercent;   // STH 수소 전환 효율 (%)
    public readonly double HydrogenProductionMmolH;// 시간당 수소 생산량 (mmol/h)
    public readonly double RecombinationRate;      // 표면 전하 재결합 손실률

    [MethodImpl(MethodImplOptions.AggressiveInlining)]
    public PecCellState(
        double biasVoltageVolt,
        double depletionWidthNm,
        double photocurrentDensityMa,
        double sthEfficiencyPercent,
        double hydrogenProductionMmolH,
        double recombinationRate)
    {
        BiasVoltageVolt = biasVoltageVolt;
        DepletionWidthNm = depletionWidthNm;
        PhotocurrentDensityMa = photocurrentDensityMa;
        SthEfficiencyPercent = sthEfficiencyPercent;
        HydrogenProductionMmolH = hydrogenProductionMmolH;
        RecombinationRate = recombinationRate;
    }
}

/// 
/// 1,000Hz 0-Alloc C# BiVO4 광양극 PEC 수전해 수치 시뮬레이터
/// 
public static class PecHydrogenEngine
{
    private const double ElectronCharge = 1.602176634e-19;
    private const double Epsilon0 = 8.8541878128e-12;
    private const double FaradayConst = 96485.33;
    private const double RelativePermittivity = 68.0;   // BiVO4 유전율
    private const double DopingDensity = 1.2e25;        // 도핑 농도 (m^-3)
    private const double FlatBandPotential = 0.22;      // V vs RHE
    private const double AbsorptionCoeff = 3.5e6;       // 450nm 흡수계수 (m^-1)
    private const double HoleDiffLength = 70e-9;        // 정공 확산 거리 (70 nm)
    private const double IncidentPhotonFlux = 2.8e21;   // AM1.5G 유효 광자속 (photons/m^2·s)
    private const double PhotoelectrodeAreaCm2 = 25.0;  // 전극 면적 (25 cm^2)

    /// 
    /// 1ms 타임스텝 인가 전위에 따른 광전기화학 수소 생산 연산 (0-Alloc)
    /// 
    [MethodImpl(MethodImplOptions.AggressiveInlining)]
    public static PecCellState StepSimulation(
        double biasVoltVsRhe)
    {
        // 1. 공핍층 폭 연산: W = sqrt(2 * eps * (V - Vfb) / (q * Nd))
        double deltaV = Math.Max(biasVoltVsRhe - FlatBandPotential, 0.0);
        double eps = Epsilon0 * RelativePermittivity;
        double wMeters = Math.Sqrt(2.0 * eps * deltaV / (ElectronCharge * DopingDensity));
        double wNm = wMeters * 1e9;

        // 2. 가르트너 식 기반 전하 수집 효율
        double expTerm = Math.Exp(-AbsorptionCoeff * wMeters);
        double gartnerFactor = 1.0 - (expTerm / (1.0 + AbsorptionCoeff * HoleDiffLength));
        double chargeCollection = Math.Clamp(gartnerFactor, 0.0, 1.0);

        // 3. 표면 산화 반응 과전압에 따른 정공 주입 효율 (버틀러-볼머형)
        double overpotential = Math.Max(biasVoltVsRhe - 0.6, 0.0);
        double injectionEff = 1.0 / (1.0 + Math.Exp(-overpotential * 6.5));

        // 4. 광전류 밀도 (mA/cm^2)
        double maxCurrentAmpsM2 = ElectronCharge * IncidentPhotonFlux;
        double currentAmpsM2 = maxCurrentAmpsM2 * chargeCollection * injectionEff;
        double currentMaCm2 = currentAmpsM2 * 0.1; // A/m^2 -> mA/cm^2

        // 5. STH 효율 (%): STH = (J_ph * 1.23V / 100mW/cm^2) * 100
        double sthPercent = (currentMaCm2 * 1.23 / 100.0) * 100.0;

        // 6. 패러데이 법칙 수소 생산 속도 (mmol/h): n_dot = (I / 2F) * 3600 * 1000
        double totalCurrentAmps = (currentMaCm2 * 1e-3) * PhotoelectrodeAreaCm2;
        double h2MmolPerSec = (totalCurrentAmps / (2.0 * FaradayConst)) * 1000.0;
        double h2MmolPerHour = h2MmolPerSec * 3600.0;

        double recombLoss = 1.0 - (chargeCollection * injectionEff);

        return new PecCellState(
            biasVoltVsRhe,
            wNm,
            currentMaCm2,
            sthPercent,
            h2MmolPerHour,
            recombLoss);
    }
}

4. 교육 및 산업적 시사점

2022 개정 교육과정 화학의 산화-환원 및 전기화학 셀 단원과 물리학의 광전효과, 에너지 밴드 갭 이론이 결합된 궁극의 친환경 수소 기술을 깊이 있게 다룹니다. 반도체 내부의 전계 형성과 전하 수집 효율이 태양광 수소 생산 수율(STH)로 직접 연결되는 물리적 과정을 1,000Hz 0-Alloc C# 코드로 실시간 탐구함으로써, 탄소중립 수소 경제를 선도할 융합적 공학 사고력을 배양합니다.