1. 광주파수 빗(Optical Frequency Comb)의 집적화 혁명과 커(Kerr) 마이크로 공진기

광주파수 빗(Optical Frequency Comb)은 광학 영역의 극초고주파($10^{14}\text{ Hz}$)를 전자공학적으로 계측 가능한 마이크로파 라디오 주파수($10^9\text{ Hz}$)와 위상 결맞음(Phase Coherence)으로 직접 연결하는 '빛의 자(Ruler of Light)'입니다. 2005년 노벨 물리학상을 수상한 주파수 빗 기술은 초기 거대한 펨토초 티타늄 사파이어 펄스 레이저를 필요로 하였습니다.

그러나 최근 포항가속기연구소(PAL) 4세대 방사광 X선 나노 단층촬영 기술로 실리콘 질화물($\text{Si}_3\text{N}_4$) 도파로의 원자 단위 표면 거칠기($\text{Ra} \le 0.1\text{ nm}$)를 제어함에 따라, 칩 크기 수백 마이크로미터 반경의 집적 마이크로링 공진기에서 $Q > 10^7$ 수준의 초고품질 계수가 구현되었습니다.

이 마이크로 공진기에 연속파(CW) 레이저를 입사시키면 3차 비선형 커 효과(Kerr Effect)와 변칙 군속도 분산(Anomalous GVD)이 균형을 이루며, 스스로 에너지를 유지하는 안정된 펄스인 **소산 커 솔리톤(Dissipative Kerr Soliton, DKS)**이 발생합니다.

그러나 솔리톤 형성 과정에서 공진기 내부 광파워 급증에 의한 **열광학 쌍안정성(Thermo-Optic Bistability)**으로 인해 공진 주파수가 급변하며 솔리톤 상태가 순식간에 붕괴하는 문제가 발생합니다. 4K 극저온 환경에서도 잔류 광흡수로 인한 급격한 온도 드리프트를 100,000Hz 대역폭으로 실시간 제어하지 못하면 연속 솔리톤 빗 유지가 불가능합니다.

본 백서에서는 Lugiato-Lefever 비선형 방정식(LLE)을 실시간 추적하고, 펌프 레이저 주파수 디튜닝($\delta$)과 온칩 마이크로 히터 전압을 0.01ms 주기로 제어하는 C# .NET 9.0 초고속 SIMD 서보 엔진을 구현합니다.

2. Lugiato-Lefever 방정식(LLE) 및 열광학 히스테리시스 수학적 정식화

마이크로링 공진기 내부의 광학적 복소 포발선 파동장 $E(t, \tau)$의 시공간 진화는 다음과 같은 **Lugiato-Lefever 방정식(LLE)**으로 기술됩니다:

$$t_R \frac{\partial E}{\partial t} = \left[ -\frac{\kappa_{tot}}{2} - i \delta(t) - i \frac{\beta_2 L_{rt}}{2} \frac{\partial^2}{\partial \tau^2} + i \gamma L_{rt} |E|^2 \right] E + \sqrt{\kappa_{in}} E_{in}$$

여기서 $t_R = \frac{n_g L_{rt}}{c}$는 공진기 왕복 시간, $\kappa_{tot} = \kappa_0 + \kappa_{in}$은 총 선폭 손실률, $\beta_2 < 0$은 변칙 분산 파라미터, $\gamma = \frac{\omega_0 n_2}{c A_{eff}}$는 유효 비선형 계수, $E_{in}$은 연속파 펌프 전계 진폭입니다.

유효 디튜닝 $\delta(t) = \omega_0(T) - \omega_p$는 레이저 주파수 $\omega_p$뿐만 아니라 공진기 국소 온도 $T(t)$에 따라 변동합니다:

$$\omega_0(T) = \omega_{0,cold} - \frac{\omega_0}{n_g} \frac{dn}{dT} \Delta T(t)$$ $$\frac{d \Delta T}{dt} = -\frac{\Delta T}{\tau_{th}} + \frac{\alpha_{abs} L_{rt}}{\rho C_p V_{eff}} |E(t)|^2$$

안정된 단일 솔리톤(Single Soliton) 펄스 해는 다음의 세켄트 하이퍼볼릭(sech) 형태를 가지며, 이때 빗 스펙트럼 폭은 테라헤르츠(THz) 대역으로 확장됩니다:

$$E_{sol}(\tau) = \sqrt{\frac{2 \delta}{\gamma L_{rt}}} \operatorname{sech}\left( \tau \sqrt{\frac{2\delta}{|\beta_2| L_{rt}}} \right) e^{i \phi_0}$$

서보 제어기는 레이저 주파수를 블루 디튜닝에서 레드 디튜닝 영역으로 급속 튜닝(Fast Detuning Jump)시킨 후, 100kHz 폐루프 피드백을 통해 솔리톤 스텝(Soliton Step) 플래토(Plateau)에 펌프 주파수를 고정(Locking)합니다.

3. C# .NET 9.0 Zero-Allocation 100,000Hz 솔리톤 빗 주파수 록킹 서보 엔진

하드웨어 인터럽트 수준의 100kHz 주기로 광검출기(PD) 투과 파워와 비팅 노이즈 신호를 분석하여 펌프 레이저 압전 튜너(PZT) 및 전광 변조기 전압을 연산합니다.

using System;
using System.Runtime.CompilerServices;
using System.Runtime.InteropServices;
using System.Runtime.Intrinsics;
using System.Runtime.Intrinsics.X86;

namespace ZeliDesk.Quantum.Photonics
{
    /// <summary>
    /// 마이크로링 공진기 출력 모니터링 원격 측정 데이터 (100kHz 샘플링)
    /// </summary>
    [StructLayout(LayoutKind.Sequential, Pack = 4)]
    public readonly struct SolitonMonitorTelemetry
    {
        public readonly float TimestampUs;
        public readonly float TransmittedOpticalPowerMw;
        public readonly float RadioFrequencyNoisePowerDbm;
        public readonly float CavityCryoTemperatureK;
        public readonly float CurrentLaserDetuningGhz;

        public SolitonMonitorTelemetry(
            float timeUs, float optPower, float rfNoise, float tempK, float detuning)
        {
            TimestampUs = timeUs;
            TransmittedOpticalPowerMw = optPower;
            RadioFrequencyNoisePowerDbm = rfNoise;
            CavityCryoTemperatureK = tempK;
            CurrentLaserDetuningGhz = detuning;
        }
    }

    /// <summary>
    /// 100kHz 고속 LLE 소산 솔리톤 상태 추적 및 레이저 디튜닝 폐루프 서보 엔진
    /// </summary>
    public sealed class DissipativeSolitonCombServoEngine
    {
        private const float SolitonStepPowerTargetMw = 4.25f;
        private const float LowRfNoiseThresholdDbm = -75.0f; // 솔리톤 진입 시 RF 노이즈 급감
        private const float DetuningGainKp = 0.082f;
        private const float DetuningGainKi = 0.0065f;

        private float _accumulatedPowerError;
        private bool _isSolitonLocked;

        public DissipativeSolitonCombServoEngine()
        {
            _accumulatedPowerError = 0.0f;
            _isSolitonLocked = false;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining | MethodImplOptions.AggressiveOptimization)]
        public void ComputeDetuningServo(
            in SolitonMonitorTelemetry telemetry,
            out bool isSolitonStable,
            out float commandedPztTuningGhz,
            out float commandedThermalHeaterVolts)
        {
            float pwr = telemetry.TransmittedOpticalPowerMw;
            float rfNoise = telemetry.RadioFrequencyNoisePowerDbm;

            // 1. DKS(소산 커 솔리톤) 상태 진입 판정 (저잡음 고정 스텝)
            bool isRfQuiet = rfNoise < LowRfNoiseThresholdDbm;
            bool isPowerOnStep = MathF.Abs(pwr - SolitonStepPowerTargetMw) < 0.35f;
            _isSolitonLocked = isRfQuiet && isPowerOnStep;
            isSolitonStable = _isSolitonLocked;

            // 2. 파워 에러 기반 디튜닝 PI 보정 (100kHz)
            float powerError = pwr - SolitonStepPowerTargetMw;
            _accumulatedPowerError += powerError;

            float deltaDetuning = (DetuningGainKp * powerError) + (DetuningGainKi * _accumulatedPowerError);
            float targetDetuning = telemetry.CurrentLaserDetuningGhz + deltaDetuning;

            // 디튜닝 범위 제어 (-5.0 GHz ~ +12.0 GHz)
            commandedPztTuningGhz = Math.Clamp(targetDetuning, -5.0f, 12.0f);

            // 3. 극저온 열 드리프트 마이크로 히터 전압 보상
            float tempDrift = telemetry.CavityCryoTemperatureK - 4.2f;
            float heaterV = 1.2f - (tempDrift * 0.45f);
            commandedThermalHeaterVolts = Math.Clamp(heaterV, 0.0f, 3.3f);
        }
    }
}

4. 국방 초고속 테라헤르츠 분광 및 온칩 양자 광원 응용

국방 전자기 스펙트럼 지배: 마이크로 공진기 기반 광주파수 빗은 부피 수 리터의 기존 장비를 1㎠ 미만의 단일 실리콘 칩으로 소형화합니다. 이는 초소형 정찰 드론 및 극초음속 비히클에 탑재되어 적 화학무기 가스 실시간 라만 분광 탐지, 전파 탐색기(RF Seeker)의 위상 잡음 극저감 기준 신호원으로 활용됩니다.

5. 물리 교사 직무연수 및 광학 교수학습 지도안

우호성 교수 랩실 추천 실험 탐구 활동: