AIDT 온디바이스 STEM 냉각 원자 광학 쌍극자 트랩(ODT) 증발 냉각 서보 엔진

1. 개요

레이저 냉각(MOT)의 도플러 및 편광 구배 한계(수 \(\mu\text{K}\))를 넘어 나노켈빈(nK) 극저온에 도달하여 보즈-아인슈타인 응축(BEC) 또는 페르미 축퇴(Fermi Degeneracy)를 달성하기 위한 최종 관문은 강하게 집속된 1064nm 적외선 레이저를 이용한 광학 쌍극자 트랩(Optical Dipole Trap, ODT) 내에서의 강제 증발 냉각(Forced Evaporative Cooling)입니다. 레이저 광출력을 지수함수적으로 감쇠시켜 트랩 깊이(\(U_0\))를 낮춤으로써 가장 에너지가 높은 원자들을 선택적으로 배출시키고, 잔류 원자들 간의 탄성 충돌(Elastic Collisions)을 통해 원자 기체의 온도를 급속 냉각시킵니다. 그러나 레이저 강도 잡음(Intensity Noise)에 의한 원자 가열과 3체 재결합 손실(Three-Body Inelastic Loss)은 위상 공간 밀도(\(\text{PSD} = n \lambda_{dB}^3\)) 상승을 저해합니다. 본 고에서는 10,000Hz 주기로 증발 경로 램프다운을 추종하는 C# 서보 엔진을 구현합니다.

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

원자 분극도 \(\alpha(\omega)\)와 레이저 빔 허리 \(w_0\) 하에서 광학 쌍극자 트랩 깊이 \(U_0\) 및 트랩 진동 주파수 \(\omega_r, \omega_z\)는 레이저 파워 \(P(t)\)에 비례합니다:

U_0(t) = \frac{2 \alpha P(t)}{\pi c \varepsilon_0 w_0^2}, \quad \omega_r = \sqrt{\frac{4 U_0}{m w_0^2}}, \quad \omega_z = \sqrt{\frac{2 U_0}{m z_R^2}}

증발 냉각 효율 파라미터 \(\gamma_{evap} = -\frac{d(\ln \text{PSD})}{d(\ln N)}\)를 극대화하는 최적 트랩 파워 램프 곡선은 오하라(O'Hara) 스케일링 법칙을 따릅니다:

P(t) = P_0 \left( 1 + \frac{t}{\tau_{evap}} \right)^{-\beta_{evap}}, \quad \eta = \frac{U_0(t)}{k_B T(t)} \approx 8 \sim 10

광세기 요동에 의한 파라메트릭 가열률 \(\Gamma_{heat} = \pi^2 \nu_{trap}^2 S_P(2\nu_{trap})\)을 억제하기 위해, 10,000Hz 초고속 음향광학 변조기(AOM) 피드백 루프로 광학 파워 노이즈 스펙트럼 밀도 \(S_P\)를 -120 dBc/Hz 이하로 고정합니다.

3. 실시간 피드백 제어 아키텍처

+-----------------------------------------------------------------------------+
|    AIDT Optical Dipole Trap Evaporative Cooling 10,000Hz Real-Time Servo    |
+-----------------------------------------------------------------------------+
|  [InGaAs Fast Photodiode Laser Power Monitor (100kHz) & Absorption Imaging] |
|                                       |                                     |
|                                       v                                     |
|     [Phase Space Density PSD & Parametric Heating Rate S_P(2\nu) Solver]    |
|                                       |                                     |
|                                       v                                     |
|     [Zero-Alloc Discrete Logarithmic Power Rampdown & Noise Damping Servo]  |
|                                       |                                     |
|                                       v                                     |
|   [Sub-Microsecond RF Driver] -> [1064nm High-Power Fiber Laser AOM]        |
+-----------------------------------------------------------------------------+

4. Zero-Allocation C# 서보 엔진 구현

ODT 레이저 램프다운 프로파일을 10,000Hz 주기로 추종하는 C# 엔진 클래스입니다.

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.UltracoldAtoms.DipoleTrap
{
    public readonly struct OdtEvapConfig
    {
        public readonly double InitialPowerW;
        public readonly double FinalPowerW;
        public readonly double TimeConstantSec;
        public readonly double DecayExponent;
        public readonly double Kp;
        public readonly double Ki;

        public OdtEvapConfig(double initialPowerW, double finalPowerW, double timeConstantSec,
                             double decayExponent, double kp, double ki)
        {
            InitialPowerW = initialPowerW;
            FinalPowerW = finalPowerW;
            TimeConstantSec = timeConstantSec;
            DecayExponent = decayExponent;
            Kp = kp;
            Ki = ki;
        }
    }

    public ref struct OdtEvaporativeServo
    {
        private readonly OdtEvapConfig _config;
        private double _powerErrorIntegral;
        private double _elapsedTimeSec;

        public OdtEvaporativeServo(in OdtEvapConfig config)
        {
            _config = config;
            _powerErrorIntegral = 0.0;
            _elapsedTimeSec = 0.0;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public (double targetPowerW, double aomDriveVoltageV) Update(double measuredPowerW, double dtSec)
        {
            _elapsedTimeSec += dtSec;

            // 이상적 증발 냉각 감쇠 곡선 연산
            double baseRamp = Math.Pow(1.0 + (_elapsedTimeSec / _config.TimeConstantSec), -_config.DecayExponent);
            double targetPower = Math.Max(_config.FinalPowerW, _config.InitialPowerW * baseRamp);

            double error = targetPower - measuredPowerW;
            _powerErrorIntegral += error * dtSec;
            _powerErrorIntegral = Math.Clamp(_powerErrorIntegral, -5.0, 5.0);

            // AOM 제어 전압 (0~10V)
            double deltaV = (_config.Kp * error) + (_config.Ki * _powerErrorIntegral);
            double aomVoltage = Math.Clamp((targetPower / _config.InitialPowerW) * 10.0 + deltaV, 0.05, 10.0);

            return (targetPower, aomVoltage);
        }
    }
}

5. 학술적 의의 및 교육적 활용 가치

본 설계는 2022 개정 물리 교과 내 '맥스웰-볼츠만 분포와 온도', '원자간 탄성 충돌과 열평형', '양자 통계와 보즈-아인슈타인 응축'을 아우르는 최정상급 현대 물리학 실험 모델입니다. 레이저 빛으로 형성된 그릇(포텐셜)의 턱을 서서히 낮추어 물질을 절대 0도에 근접시키는 정교한 양자 냉각 기술을 직접 경험할 수 있습니다.