AIDT 온디바이스 STEM 스퀴즈드 광(Squeezed Light) 간섭계 표준 양자 한계(SQL) 극복 서보 엔진

1. 개요

차세대 거대 레이저 간섭계 중력파 검출기(LIGO, Virgo, KAGRA 및 Cosmic Explorer)는 블랙홀 병합과 중성자별 충돌에 의해 시공간에 발생하는 \(10^{-21}\) 미만의 극미세 변형률(Strain \(h\))을 포착합니다. 그러나 간섭계 감도는 하이젠베르크 불확정성 원리(\(\Delta X_1 \Delta X_2 \ge 1/4\))에 기인한 광학 양자 잡음의 벽에 가로막혀 있습니다. 고주파수 영역에서는 광자수 요동인 샷 노이즈(Shot Noise)가 지배하고, 저주파수 영역에서는 광자가 거울을 때리는 광 복사압 잡음(Quantum Radiation Pressure Noise)이 미러를 요동치게 만듭니다. 이를 동시에 돌파하기 위해 광학 파라메트릭 발진기(OPO)에서 생성된 직교위상 압축 진공(Squeezed Vacuum) 상태를 300m 급 디튠 필터 공동(Detuned Filter Cavity)에 통과시켜 주파수에 따라 압축 위상각(\(\theta_{sq}(\Omega)\))을 90도 연속 회전시키는 주파수 의존 스퀴징(Frequency-Dependent Squeezing) 기술이 필수적입니다. 본 고에서는 필터 공동의 디튜닝 공진 주파수를 10,000Hz 주기로 추종하며 양자 역작용(Quantum Back-action)을 상쇄하는 0-Alloc C# 엔진을 구현합니다.

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

간섭계 암 캐비티(Arm Cavity) 내의 등가 변위 잡음 전력 스펙트럼 밀도 \(S_h(\Omega)\)와 표준 양자 한계(SQL) \(S_h^{SQL}(\Omega)\)는 다음과 같습니다:

여기서 광-기계 결합 파라미터인 Kimble 인자 \(\mathcal{K}(\Omega)\)는 간섭계 레이저 파워 \(P_0\), 미러 질량 \(m\), 캐비티 반치폭 \(\gamma_{arm}\)에 의해 결정됩니다:

이상적인 주파수 의존 압축 위상각 \(\theta_{opt}(\Omega)\)는 저주파수(방사압 영역)에서의 위상 직교위상 압축과 고주파수(샷 노이즈 영역)에서의 진폭 직교위상 압축을 매끄럽게 연결합니다:

길이 \(L_{fc}\), 디튜닝 \(\Delta \omega_{fc}\), 선폭 \(\gamma_{fc}\)를 갖는 광학 필터 공동에 반사된 압축광의 주파수 의존 위상 회전각 \(\Phi_{fc}(\Omega)\)는 다음과 같이 달성됩니다:

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

다음은 필터 공동 공진 오차 신호(PDH 복조)를 10,000Hz로 연산하여 피에조 미러를 제어하고 최적 압축각을 고정하는 0-Alloc C# 엔진 구현체입니다:

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Core.QuantumOptics
{
    public ref struct SqueezedLightInterferometerServo
    {
        private readonly double _dt;
        private readonly double _filterCavityLinewidthHz;
        private readonly double _targetDetuningHz;
        private double _currentDetuningHz;
        private double _piezoVoltageCmd;
        private double _squeezingLevelDb;
        
        // PDH 오차 적분기
        private double _pdhErrorFiltered;
        private double _integralError;
        private readonly double _kp;
        private readonly double _ki;

        public SqueezedLightInterferometerServo(
            double dt,
            double filterCavityLinewidthHz,
            double targetDetuningHz,
            double squeezingLevelDb,
            double kp,
            double ki)
        {
            _dt = dt;
            _filterCavityLinewidthHz = filterCavityLinewidthHz; // 약 50 Hz
            _targetDetuningHz = targetDetuningHz;               // 약 35 Hz
            _squeezingLevelDb = squeezingLevelDb;               // 10 dB 스퀴징
            _kp = kp;
            _ki = ki;
            _currentDetuningHz = targetDetuningHz;
            _piezoVoltageCmd = 0.0;
            _pdhErrorFiltered = 0.0;
            _integralError = 0.0;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public double UpdateCavityLock(
            double pdhDemodulatedErrorVolts,
            double opticalPowerInCavityWatts,
            out double rotationAngleDegAt100Hz,
            out double quantumNoiseReductionDb)
        {
            // 1. PDH 록인 오차 신호 1차 저역통과 필터링
            _pdhErrorFiltered = 0.85 * _pdhErrorFiltered + 0.15 * pdhDemodulatedErrorVolts;

            // 2. PI 서보 제어를 통한 필터 공동 피에조 미러 전압 보정
            _integralError += _pdhErrorFiltered * _dt;
            _integralError = Math.Clamp(_integralError, -10.0, 10.0);

            double deltaV = _kp * _pdhErrorFiltered + _ki * _integralError;
            _piezoVoltageCmd += deltaV;
            _piezoVoltageCmd = Math.Clamp(_piezoVoltageCmd, -75.0, 75.0); // ±75V 피에조 드라이버

            // 3. 피에조 변위에 따른 순시 디튜닝 주파수 변화 계산
            double piezoFreqShiftPerVolt = 1.25; // Hz/V
            _currentDetuningHz = _targetDetuningHz + (_piezoVoltageCmd * piezoFreqShiftPerVolt);

            // 4. 100Hz(중력파 관측 핵심 대역)에서의 압축 위상 회전각 계산
            double omega = 100.0; // Hz
            double gamma = _filterCavityLinewidthHz;
            double delta = _currentDetuningHz;

            double phi1 = Math.Atan2(omega - delta, gamma);
            double phi2 = Math.Atan2(omega + delta, gamma);
            double totalRotationRad = Math.PI - (phi1 + phi2);
            rotationAngleDegAt100Hz = (totalRotationRad * 180.0) / Math.PI;

            // 5. 양자 잡음 저감도 dB 계산 (위상 오차 페널티 포함)
            double phaseError = Math.Abs(totalRotationRad - (Math.PI / 4.0)); // 45도 최적각 대비 편차
            double degradationFactor = Math.Cos(phaseError);
            quantumNoiseReductionDb = _squeezingLevelDb * Math.Max(0.1, degradationFactor);

            return _piezoVoltageCmd;
        }

        public readonly double CurrentDetuningHz => _currentDetuningHz;
        public readonly double PiezoDriveVoltage => _piezoVoltageCmd;
    }
}

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