AIDT 온디바이스 STEM X선 자유 전자 레이저(XFEL) 펨토초 펄스 위상 정렬 피에조 서보 엔진

1. 개요

X선 자유전자레이저(X-ray Free-Electron Laser, XFEL)는 원자 단위의 공간 분해능(옹스트롬, \(\text{Å}\))과 분자 진동의 시간 분해능(펨토초, \(10^{-15}\text{s}\))을 동시에 제공하는 궁극의 꿈의 현미경입니다. 화학 결합의 파괴와 재결합, 단백질 분자의 실시간 구조 동역학을 촬영하는 펌프-프로브(Pump-Probe) 실험에서는 광학 레이저 펄스와 XFEL X선 펄스 사이의 상대적 도착 시간 지터(Arrival Time Jitter)가 수 펨토초 이내로 엄격히 유지되어야 합니다. 그러나 수 킬로미터에 달하는 가속관의 온도 드리프트, RF 위상 요동, 미세 진동으로 인해 수십~수백 펨토초의 타이밍 오차가 불가피하게 발생합니다. 본 고에서는 전광 샘플링(Electro-Optic Sampling, EOS) 결정을 통과한 레이저의 복굴절 위상차를 계측하고, 광경로 딜레이 라인의 고속 압전 액추에이터를 10,000Hz 주기로 능동 피드백 제어하는 0-Alloc C# 서보 엔진을 구현합니다.

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

전자빔이 언듈레이터를 통과하며 포켈스 효과(Pockels Effect)를 갖는 비선형 광학 결정(ZnTe 또는 GaP)을 스쳐 지나갈 때, 전자 다발의 쿨롱 전기장 \(E_{coulomb}(t)\)에 의해 유도되는 광학 위상차 \(\Gamma(t)\)는 다음과 같습니다:

프로브 광학 펄스가 편광 빔 스플리터를 거쳐 균형 광검출기(Balanced Photodetector, BPD)로 수신될 때 검출되는 차동 전압 신호 \(\Delta V_{BPD}(t)\)와 도착 시간 오차 \(\tau_{arr}\)의 관계식은 다음과 같습니다:

광경로 피에조 딜레이 스테이지의 전압 구동 역학 모델은 고유진동수 \(\omega_n\), 감쇠비 \(\zeta\), 피에조 변위 상수 \(d_{piezo}\)를 포함한 2차 시스템 미분방정식으로 지배됩니다:

광경로 길이 보정량 \(\Delta L = 2 \Delta x_{piezo}\)에 따른 시간 지연 보정량 \(\Delta \tau_{comp} = 2 \Delta x_{piezo} / c\)과 폐루프 오차 최소화 목적함수는 다음과 같습니다:

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

다음은 전광 샘플링 BPD 신호를 10,000Hz로 판독하여 펨토초 타이밍 지터를 나노미터 피에조 변위로 즉각 상쇄하는 무할당 C# 엔진 구현체입니다:

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Core.Accelerators
{
    public ref struct XfelFemtosecondPulsePhaseServo
    {
        private readonly double _dt;
        private readonly double _speedOfLight; // m/s
        private readonly double _piezoGainMPerV;
        private readonly double _eosSensitivityVPerFs;
        
        // 칼만 필터 및 지터 추정기 상태
        private double _estimatedJitterFs;
        private double _jitterVelocityFsPerS;
        private double _piezoVoltageCmd;
        private double _piezoPositionM;
        
        // 제어기 이득
        private readonly double _kp;
        private readonly double _ki;
        private double _integralErr;

        public XfelFemtosecondPulsePhaseServo(
            double dt,
            double piezoGainMPerV,
            double eosSensitivityVPerFs,
            double kp,
            double ki)
        {
            _dt = dt;
            _speedOfLight = 299792458.0;
            _piezoGainMPerV = piezoGainMPerV;
            _eosSensitivityVPerFs = eosSensitivityVPerFs;
            _kp = kp;
            _ki = ki;
            _estimatedJitterFs = 0.0;
            _jitterVelocityFsPerS = 0.0;
            _piezoVoltageCmd = 0.0;
            _piezoPositionM = 0.0;
            _integralErr = 0.0;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public double Update(
            double bpdDifferentialVoltageV,
            double rfRefPhaseRad,
            out double residualJitterFs,
            out double piezoDisplacementNm)
        {
            // 1. 전광 샘플링 전압으로부터 순간 도착 시간 오차(fs) 계산
            double measuredArrivalErrorFs = bpdDifferentialVoltageV / Math.Max(1e-6, _eosSensitivityVPerFs);

            // 2. 2상태 알파-베타 추정 필터링 (가속관 진동 및 열 드리프트 분리)
            double alpha = 0.25;
            double beta = 0.05;
            double predJitter = _estimatedJitterFs + _jitterVelocityFsPerS * _dt;
            double innovation = measuredArrivalErrorFs - predJitter;

            _estimatedJitterFs = predJitter + alpha * innovation;
            _jitterVelocityFsPerS += (beta / _dt) * innovation;

            residualJitterFs = innovation;

            // 3. 딜레이 보정 피에조 필요 변위 계산: Delta_x = (c * Delta_tau) / 2
            double requiredDelaySec = _estimatedJitterFs * 1e-15;
            double targetDisplacementM = (requiredDelaySec * _speedOfLight) / 2.0;

            // 4. PI 서보 제어기 연산
            double errorM = targetDisplacementM - _piezoPositionM;
            _integralErr += errorM * _dt;
            _integralErr = Math.Clamp(_integralErr, -1e-6, 1e-6); // 안티 와인드업

            double deltaV = _kp * errorM + _ki * _integralErr;
            _piezoVoltageCmd += deltaV;
            _piezoVoltageCmd = Math.Clamp(_piezoVoltageCmd, -150.0, 150.0); // ±150V 드라이버

            // 피에조 액추에이터 동역학 모사 (1차 저역통과)
            double actualTargetM = _piezoVoltageCmd * _piezoGainMPerV;
            _piezoPositionM += (_dt / (_dt + 1e-4)) * (actualTargetM - _piezoPositionM);

            piezoDisplacementNm = _piezoPositionM * 1e9;
            return _piezoVoltageCmd;
        }

        public readonly double EstimatedJitterFemtoseconds => _estimatedJitterFs;
        public readonly double CurrentPiezoPositionNm => _piezoPositionM * 1e9;
    }
}

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