AIDT 온디바이스 STEM 원자물리 광학 포획 단일 원자 트위저(Optical Tweezer) 음향광학 편향기(AOD) 위치 지터 서보 제어기

1. 개요

중성 원자(Neutral Atom) 기반 양자 컴퓨터 및 양자 시뮬레이터는 고진공 챔버 내부에서 레이저 냉각된 개별 루비듐(\(^{87}\text{Rb}\)) 또는 세슘(\(^{133}\text{Cs}\)) 원자를 수 마이크로미터 간격의 정밀 2차원/3차원 배열로 포획하기 위해 광학 다이폴 트랩(Optical Tweezer)을 사용합니다. 극초단 레이저 빔을 음향광학 편향기(Acousto-Optic Deflector, AOD)에 통과시킨 후 고개구수(High-NA) 현미경 대물렌즈로 집속하면, AOD에 인가되는 고주파(RF) 신호의 주파수에 비례하여 집속 초점의 공간 위치를 임의로 조작할 수 있습니다. 그러나 빔 전달 광학계 미러의 기계적 열팽창 표류 및 서브-마이크론 진동은 트랩 위치의 지터(Position Jitter)를 유발하여 원자의 파라메트릭 가열(Parametric Heating)과 트랩 이탈을 초래하며, 리드버그(Rydberg) 양자 게이트 충실도(Fidelity)를 급격히 저하시킵니다. 본 고에서는 10,000Hz 주기로 위치 검출기(PSD) 및 EMCCD 형광 신호를 처리하여 AOD RF 주파수를 밀리헤르츠 단위로 실시간 보정하는 0-Alloc C# 엔진을 구현합니다.

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

가우시안 빔 초점에서의 광학 다이폴 포텐셜 \(U_{dipole}\) 및 트랩 깊이 방정식은 다음과 같습니다:

U_{dip}(\vec{r}) = -\frac{3\pi c^2}{2\omega_0^3} \left( \frac{\Gamma}{\omega_0 - \omega_L} + \frac{\Gamma}{\omega_0 + \omega_L} \right) I(\vec{r}), \quad I(r, z) = \frac{2 P_{laser}}{\pi w(z)^2} \exp\left(-\frac{2 r^2}{w(z)^2}\right)

음향광학 편향기(AOD)의 브래그 1차 회절각 \(\theta_{diff}\) 및 대물렌즈 초점면 원자 트랩 위치 \(x_{trap}\) 관계식은 다음과 같습니다:

\theta_{diff} = \frac{\lambda_{laser} f_{RF}}{v_{acoustic}} \implies x_{trap} = f_{obj} \frac{\lambda_{laser}}{v_{acoustic}} f_{RF} = K_{AOD} \cdot f_{RF}

미세 위치 편차 \(\Delta x_{jitter}\)를 영점화하기 위한 AOD 구동 RF 주파수 보정 비례식은 다음과 같습니다:

\Delta f_{RF}^{cmd}(t) = -\frac{1}{K_{AOD}} \left( K_p \Delta x(t) + K_i \int \Delta x(t) dt + K_d \frac{d\Delta x}{dt} \right)

트랩 강도 변동에 의한 원자의 파라메트릭 가열률 \(\Gamma_{heat}\) 억제 조건식은 다음과 같습니다:

\Gamma_{heat} = \pi^2 \nu_{trap}^2 S_I(2\nu_{trap}) \le 1.0\text{ phonon/s} \quad \left( S_I \le 10^{-8}\text{ Hz}^{-1} \right)

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

다음은 4분할 포토다이오드(QPD) 빔 위치 오차와 단일 원자 형광 카운트를 10,000Hz 주기로 분석하여 DDS(Direct Digital Synthesizer) RF 주파수 레지스터를 실시간 갱신하는 0-Alloc C# 엔진 구현체입니다:

OpticalTweezerAodServoEngine.cs

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.OpticalTweezer
{
    public readonly ref struct TweezerSensorFrame
    {
        public readonly double BeamPositionXNanometers;
        public readonly double BeamPositionYNanometers;
        public readonly double TargetPositionXNanometers;
        public readonly double TargetPositionYNanometers;
        public readonly double FluorescenceCountRateKhz;
        public readonly double TrapLaserPowerMw;

        public TweezerSensorFrame(
            double beamPositionXNanometers,
            double beamPositionYNanometers,
            double targetPositionXNanometers,
            double targetPositionYNanometers,
            double fluorescenceCountRateKhz,
            double trapLaserPowerMw)
        {
            BeamPositionXNanometers = beamPositionXNanometers;
            BeamPositionYNanometers = beamPositionYNanometers;
            TargetPositionXNanometers = targetPositionXNanometers;
            TargetPositionYNanometers = targetPositionYNanometers;
            FluorescenceCountRateKhz = fluorescenceCountRateKhz;
            TrapLaserPowerMw = trapLaserPowerMw;
        }
    }

    public ref struct TweezerControlCommand
    {
        public double RfFrequencyXShiftKHz;
        public double RfFrequencyYShiftKHz;
        public double CommandedLaserPowerMw;
        public double PositionErrorMagnitudeNm;
        public bool IsAtomTrapped;
        public bool IsPositionSubPixelLocked;
    }

    public ref struct OpticalTweezerAodServoEngine
    {
        private const double KaodNmPerKhz = 12.5; // 12.5 nm / kHz
        private const double BaseLaserPowerMw = 15.0;
        private readonly double _kp;
        private readonly double _ki;
        private readonly double _kd;
        private readonly double _dt;
        private double _integralX;
        private double _integralY;
        private double _prevErrX;
        private double _prevErrY;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public OpticalTweezerAodServoEngine(double kp, double ki, double kd, double dt)
        {
            _kp = kp;
            _ki = ki;
            _kd = kd;
            _dt = dt;
            _integralX = 0.0;
            _integralY = 0.0;
            _prevErrX = 0.0;
            _prevErrY = 0.0;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public void ComputeTweezerPositionServo(in TweezerSensorFrame frame, out TweezerControlCommand cmd)
        {
            cmd = default;
            double errX = frame.TargetPositionXNanometers - frame.BeamPositionXNanometers;
            double errY = frame.TargetPositionYNanometers - frame.BeamPositionYNanometers;

            _integralX += errX * _dt;
            _integralY += errY * _dt;
            _integralX = Math.Clamp(_integralX, -200.0, 200.0);
            _integralY = Math.Clamp(_integralY, -200.0, 200.0);

            double dErrX = (errX - _prevErrX) / _dt;
            double dErrY = (errY - _prevErrY) / _dt;
            _prevErrX = errX;
            _prevErrY = errY;

            // Frequency shift calculation in kHz
            double shiftX = ((_kp * errX) + (_ki * _integralX) + (_kd * dErrX)) / KaodNmPerKhz;
            double shiftY = ((_kp * errY) + (_ki * _integralY) + (_kd * dErrY)) / KaodNmPerKhz;

            shiftX = Math.Clamp(shiftX, -500.0, 500.0);
            shiftY = Math.Clamp(shiftY, -500.0, 500.0);

            double errMag = Math.Sqrt(errX * errX + errY * errY);

            cmd.RfFrequencyXShiftKHz = shiftX;
            cmd.RfFrequencyYShiftKHz = shiftY;
            cmd.CommandedLaserPowerMw = BaseLaserPowerMw;
            cmd.PositionErrorMagnitudeNm = errMag;
            cmd.IsAtomTrapped = frame.FluorescenceCountRateKhz >= 8.0; // Single atom threshold
            cmd.IsPositionSubPixelLocked = errMag <= 10.0; // Under 10nm jitter
        }
    }
}

4. 교육 및 실무 적용 방안

한국방송통신대학교 양자정보/원자물리학: 중성 원자 트랩 원리, AOD 음향광학 회절 및 리드버그 양자 큐비트 배열 제어 교구.
한국표준과학연구원(KRISS) 및 KIST 양자컴퓨팅연구단: 100큐비트급 중성원자 양자컴퓨터 광학 트위저 서보 탑재 실장.