AIDT 온디바이스 STEM 광학 속삭임의 회랑 모드(WGM) 마이크로 공진기 열-광학 쌍안정성 서보 엔진

1. 개요: 속삭임의 회랑 모드(WGM)와 극초고감도 광자 포집

세인트 폴 대성당의 돔 벽면을 따라 작은 속삭임이 반대편까지 또렷하게 전달되는 음향 현상에서 착안된 속삭임의 회랑 모드(Whispering Gallery Mode, WGM) 광학 공진기는 실리카 마이크로토로이드(Microtoroid)나 구형 미소구체 내부에서 빛이 연속 전반사(TIR)를 거치며 직경 수십 마이크로미터의 미소 체적 내에 광자를 수십만 번 순환 가둡니다. 품위 계수(Quality Factor) \(Q > 10^8\)에 달하는 극한의 광 에너지 집속으로 인해 미약한 입력 레이저에서도 광자 흡수에 의한 줄 발열과 열-광학 효과(Thermo-optic Effect, \(dn/dT > 0\))가 공진 피크를 비대칭으로 찌그러뜨리는 열-광학 쌍안정성(Thermal Bistability)이 발생합니다.

2. 결합 모드 이론(CMT) 및 열역학 공진 편이 수식 모델

테이퍼 광섬유와 마이크로 공진기 사이의 결합 광파 진폭 \(a(t)\)는 시간 종속 결합 모드 방정식(Coupled Mode Theory)을 따릅니다:

\frac{da}{dt} = \left( i(\omega_L - \omega_0 - \Delta \omega_{th}) - \frac{\gamma_0 + \gamma_{ex}}{2} \right) a + i \sqrt{\gamma_{ex}} s_{in}

여기서 \(\omega_L\)은 입사 레이저 주파수, \(\omega_0\)는 고유 공진 주파수, \(\gamma_0\)는 진성 손실률, \(\gamma_{ex}\)는 테이퍼 결합 손실률입니다. 공진기 내부의 순환 광자수 \(|a|^2\)에 비례하여 온도 상승 \(\Delta T\)가 유도되며, 공진 주파수는 장파장(적색 편이) 쪽으로 밀려납니다:

\frac{d\Delta T}{dt} = -\frac{\Delta T}{\tau_{th}} + \frac{\gamma_0 |a|^2}{C_{th}}, \quad \Delta \omega_{th} = -\omega_0 \left( \frac{1}{n} \frac{dn}{dT} + \alpha_L \right) \Delta T

투과율 \(T_{opt} = |s_{out}|^2 / |s_{in}|^2\)은 전형적인 비대칭 삼각형 로렌츠 왜곡 곡선을 형성하며, 레이저 주파수를 적색에서 청색으로 주파수 스캔할 때 점프 불연속이 발생합니다.

3. 0-Alloc C# WGM 공진 주파수 능동 잠금 서보 구현

열 편이를 예측하고 레이저 피에조 파장을 공진점 최저점에 10,000Hz 주기로 고정(Lock-in)하는 ref struct 아키텍처입니다:

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.Nanophotonics
{
    public readonly ref struct WgmResonatorState
    {
        public readonly float OpticalTransmission;
        public readonly float CavityTempDeltaK;
        public readonly float ThermalShiftGhz;
        public readonly bool IsLocked;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public WgmResonatorState(float trans, float dtK, float shiftGhz, bool locked)
        {
            OpticalTransmission = trans;
            CavityTempDeltaK = dtK;
            ThermalShiftGhz = shiftGhz;
            IsLocked = locked;
        }
    }

    public ref struct WgmMicroresonatorServo
    {
        private const float Gamma0 = 1.2e8f;  // Intrinsic damping (rad/s), Q ~ 1.5e8
        private const float GammaEx = 1.0e8f; // External coupling (near critical coupling)
        private const float TauThSec = 0.0012f; // 1.2 ms thermal relaxation
        private const float ThermoOpticCoeff = 1.2e-5f; // dn/dT of silica (1/K)

        private float _deltaTK;
        private float _laserDetuningGhz;

        public WgmMicroresonatorServo(float initialDetuningGhz = -0.5f)
        {
            _deltaTK = 0.0f;
            _laserDetuningGhz = initialDetuningGhz;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public WgmResonatorState UpdateServo(float laserPowerMw, float targetDetuningGhz, float dt = 0.0001f)
        {
            _laserDetuningGhz += (targetDetuningGhz - _laserDetuningGhz) * (dt * 100.0f);

            // Resonant intra-cavity energy: |a|^2
            float totalGamma = Gamma0 + GammaEx;
            float thermalShiftRadS = (_deltaTK * ThermoOpticCoeff * 1.93e14f); // 1550nm omega0
            float effectiveDetuningRadS = (_laserDetuningGhz * 1.0e9f * 2.0f * MathF.PI) - thermalShiftRadS;

            float denom = (effectiveDetuningRadS * effectiveDetuningRadS) + (totalGamma * totalGamma * 0.25f);
            float intraCavityEnergy = (GammaEx * laserPowerMw * 1e-3f) / denom;

            // Thermal equation: dDeltaT/dt = -DeltaT/tau + P_abs / C_th
            float heatGen = Gamma0 * intraCavityEnergy * 8.5e3f; // conversion factor to K/s
            _deltaTK += (- (_deltaTK / TauThSec) + heatGen) * dt;

            // Optical transmission |s_out|^2 / |s_in|^2
            float trans = 1.0f - ((Gamma0 * GammaEx) / denom);
            trans = MathF.Clamp(trans, 0.01f, 1.0f);

            float shiftGhz = thermalShiftRadS / (2.0f * MathF.PI * 1e9f);
            bool locked = trans < 0.15f; // Resonance dip locked

            return new WgmResonatorState(trans, _deltaTK, shiftGhz, locked);
        }
    }
}

4. 실시간 WGM 열-광학 쌍안정성 및 로렌츠 피크 왜곡 시뮬레이터

입사 레이저 출력과 디튜닝을 조절하며 대칭형 로렌츠 딥이 열효과에 의해 적색 편이된 비대칭 톱니형(Sawtooth) 쌍안정 곡선으로 변화하는 과정을 관측하십시오.

입사 레이저 출력 (Laser Power): 2.5 mW 레이저 파장 디튜닝 (Laser Detuning Δf): +0.20 GHz

5. 교육 연계 및 첨단 양자 나노포토닉스 실무적 의의

2022 개정 교육과정 물리학 II: 전반사와 광섬유 전파, 간섭과 공명(Fabry-Perot 및 원형 링 공진기), 빛의 흡수와 열역학적 굴절률 변화.
화학 II: 단일 바이러스/단백질 나노입자 부착에 따른 미세 모드 분할(Mode Splitting) 분자 감지 메커니즘.
차세대 광 집적회로 / 국방 전자전: 온칩 초저위상잡음 마이크로콤(Microcomb) 광학 주파수 빗, 위성 레이저 자이로스코프.