AIDT 온디바이스 STEM 광전자 방출 분광학(ARPES) 2D 페르미 표면 매핑 서보 엔진

1. 개요

양자 물질, 위상 절연체 및 고온 초전도체의 전자 구조를 규명하는 가장 직접적이고 강력한 실험 기법은 각도 분해 광전자 분광학(Angle-Resolved Photoemission Spectroscopy, ARPES)입니다. 단색 극자외선(EUV) 또는 연X선 광자를 고체 표면에 조사하여 방출되는 광전자의 운동 에너지(\(E_{kin}\))와 방출 각도(\(\theta, \phi\))를 측정함으로써, 결정 내부 전자의 결합 에너지와 결정 운동량(\(\mathbf{k}_{\|}\))을 직접 매핑합니다. 그러나 고분해능 페르미 표면 토폴로지를 획득하기 위해서는 10K 이하의 극저온 6축 고니오미터가 각도 오차 0.01도 미만으로 정밀 회전해야 하며, 지구 자기장 및 미세 표면 대전(Charging)에 의한 전자 궤적 굴절을 실시간 정전기 편향판(Deflector)으로 동적 보상해야 합니다. 본 고에서는 2차원 검출기 스펙트럼 강도 행렬로부터 페르미 준위 교차점(\(\mathbf{k}_F\))을 추출하고 편향 전압을 10,000Hz 주기로 제어하는 0-Alloc C# 엔진을 구현합니다.

2. 수학적 모델링 및 지배 방정식

아인슈타인의 광전효과 및 표면 평행 운동량 보존 법칙에 따른 전자 결합 에너지 \(E_B\)와 결정 평행 운동량 \(k_{\|}\)는 다음과 같습니다:

반구형 정전 에너지 분석기(Hemispherical Deflector Analyzer) 내부의 전자 궤적 반경 \(R_0\) 및 통과 에너지(Pass Energy) \(E_p\)에 따른 미분 분해능 \(\Delta E\)와 각도 분해능 \(\Delta \theta\)는 다음과 같습니다:

페르미-디락 분포 \(f(E, T) = [1 + \exp((E - E_F)/k_B T)]^{-1}\)와 단일 전자 스펙트럼 함수 \(A(\mathbf{k}, \omega)\)의 곱으로 모델링되는 광전자 강도 분포식은 다음과 같습니다:

전자 빔 궤적 편향 각도 \(\theta_{def}\)를 조절하는 정전 편향판 전압 \(V_{def}\)와 잔류 표면 전하 포텐셜 \(\Phi_{surf}\) 보상 방정식은 다음과 같습니다:

3. 실시간 C# 서보 제어기 구현

다음은 2D ARPES 에너지-각도 프레임 데이터를 10,000Hz 주기로 분석하여 페르미 표면 횡단 좌표를 추종하고 편향 전압을 보정하는 0-Alloc C# 엔진 구현체입니다:

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Core.QuantumMaterials
{
    public ref struct ArpesFermiSurfaceServo
    {
        private readonly double _dt;
        private readonly double _photonEnergyEv;
        private readonly double _workFunctionEv;
        private readonly double _passEnergyEv;
        private double _targetBindingEnergyEv;
        
        // 정전 편향 상태
        private double _currentDeflectorVoltageV;
        private double _surfaceChargeOffsetV;
        private double _goniometerThetaDeg;
        private double _measuredKineticEnergyEv;
        
        // PI 제어기 이득
        private readonly double _kp;
        private readonly double _ki;
        private double _integralErr;

        public ArpesFermiSurfaceServo(
            double dt,
            double photonEnergyEv,
            double workFunctionEv,
            double passEnergyEv,
            double kp,
            double ki)
        {
            _dt = dt;
            _photonEnergyEv = photonEnergyEv; // 예: 21.2 eV (He-I)
            _workFunctionEv = workFunctionEv; // 약 4.3 eV
            _passEnergyEv = passEnergyEv;     // 약 10.0 eV
            _targetBindingEnergyEv = 0.0;     // Fermi level E_F
            _kp = kp;
            _ki = ki;
            _currentDeflectorVoltageV = 0.0;
            _surfaceChargeOffsetV = 0.0;
            _goniometerThetaDeg = 0.0;
            _measuredKineticEnergyEv = photonEnergyEv - workFunctionEv;
            _integralErr = 0.0;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public double UpdateFermiTracking(
            double detectedPeakKineticEnergyEv,
            double measuredEmissionAngleDeg,
            double targetWavevectorKxInvA,
            out double currentKxInvA,
            out double fermiEnergyDeviationEv)
        {
            _measuredKineticEnergyEv = detectedPeakKineticEnergyEv;

            // 1. 순시 페르미 준위 편차 계산 (E_kin_peak - (h*nu - Phi))
            double theoreticalFermiKinetic = _photonEnergyEv - _workFunctionEv;
            fermiEnergyDeviationEv = _measuredKineticEnergyEv - theoreticalFermiKinetic;

            // 표면 대전에 의한 포텐셜 드리프트 적분
            _surfaceChargeOffsetV += 0.01 * fermiEnergyDeviationEv;

            // 2. 평행 운동량 k_x 역산 (Inv Angstrom)
            // k_x = 0.512 * sqrt(E_kin) * sin(theta_rad)
            double thetaRad = (measuredEmissionAngleDeg * Math.PI) / 180.0;
            double sqrtE = Math.Sqrt(Math.Max(0.1, _measuredKineticEnergyEv));
            currentKxInvA = 0.5123167 * sqrtE * Math.Sin(thetaRad);

            // 3. 목표 파수 k_x 오차 도출
            double kxError = targetWavevectorKxInvA - currentKxInvA;

            // 4. 고속 정전 편향판 전압 PI 서보 연산
            _integralErr += kxError * _dt;
            _integralErr = Math.Clamp(_integralErr, -0.5, 0.5);

            double deltaV = _kp * kxError + _ki * _integralErr;
            _currentDeflectorVoltageV += deltaV;
            _currentDeflectorVoltageV = Math.Clamp(_currentDeflectorVoltageV, -150.0, 150.0);

            return _currentDeflectorVoltageV;
        }

        public readonly double CurrentDeflectorVoltage => _currentDeflectorVoltageV;
        public readonly double EstimatedSurfacePotential => _surfaceChargeOffsetV;
    }
}

4. 교육 및 국방 현장 적용 방안