AIDT 온디바이스 STEM 테라헤르츠(THz) 양자 폭포 레이저 주파수 안정화 서보 제어기

1. 개요

전자기파 스펙트럼에서 광학과 마이크로파 사이의 '테라헤르츠 간극(THz Gap, 0.1~10 THz)'을 메우는 테라헤르츠 양자 폭포 레이저(THz-QCL)는 차세대 6G 무선 통신, 천체물리학 분광 관측, 비파괴 물질 투과 영상의 핵심 광원입니다. 반도체 초격자(Superlattice) 내 전자 서브밴드 간 전이(Intersubband Transition)를 이용하는 THz-QCL은 높은 광출력을 제공하지만, 극저온 크라이오스탯의 기계적 진동, 주입 전류 잡음, 온도 드리프트로 인해 자유 발진(Free-running) 시 수 MHz에 달하는 넓은 선폭(Linewidth)과 주파수 요동을 겪습니다. 초고분해능 헤테로다인 분광을 실현하기 위해서는 쇼트키 다이오드 기반 고조파 믹서(Harmonic Mixer)를 통해 THz 레이저 광과 마이크로파 기준 발진기의 고조파를 맥놀이(Beat)시키고, 위상 동기 루프(PLL)를 통해 주입 전류 및 미러 피에조를 초고속 보정해야 합니다. 본 고에서는 10,000Hz 주기로 주파수 드리프트를 추종하고 선폭을 서브-kHz로 압축하는 0-Alloc C# 엔진을 구현합니다.

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

THz-QCL 광학 발진 주파수 \(\nu_{thz}\)와 마이크로파 기준 주파수 \(f_{rf}\)의 제 \(N\)차 고조파 사이에서 검출되는 중간 주파수(IF) 비트 신호 \(f_{IF}\)는 다음과 같습니다:

주입 바이어스 전류 \(I_{bias}\) 및 크라이오스탯 온도 \(T\)에 따른 레이저 주파수 동조율(Tuning Coefficient) 방정식은 다음과 같습니다:

주파수 위상 검파기(PFD)의 출력 오차 전압 \(V_{pfd}\)와 주파수 노이즈 전력 스펙트럼 밀도 \(S_{\delta \nu}(f)\)에 따른 레이저 실효 선폭(FWHM) \(\Delta \nu_{laser}\)의 Schawlow-Townes 한계는 다음과 같습니다:

전류 고속 피드백 루프와 온도/피에조 저속 피드백 루프의 분할 전달함수 제어 법칙은 다음과 같습니다:

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

다음은 IF 비트 주파수 오차를 10,000Hz로 판독하여 바이어스 전류 변조 신호와 피에조 전압을 무할당으로 산출하는 C# 서보 엔진 구현체입니다:

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Core.Terahertz
{
    public ref struct ThzQclFrequencyStabilizationServo
    {
        private readonly double _dt;
        private readonly double _targetIfFreqHz; // 100.0 MHz
        private double _currentIfFreqHz;
        private double _cmdCurrentModulationMa;
        private double _cmdPztVoltageV;
        
        // 이중 루프 필터 상태
        private double _freqErrorFilteredHz;
        private double _phaseIntegral;
        private readonly double _kpCurrent;
        private readonly double _kiPzt;

        public ThzQclFrequencyStabilizationServo(
            double dt,
            double targetIfFreqHz,
            double kpCurrent,
            double kiPzt)
        {
            _dt = dt;
            _targetIfFreqHz = targetIfFreqHz;
            _kpCurrent = kpCurrent; // mA / Hz
            _kiPzt = kiPzt;         // V / (Hz * s)
            _currentIfFreqHz = targetIfFreqHz;
            _cmdCurrentModulationMa = 0.0;
            _cmdPztVoltageV = 0.0;
            _freqErrorFilteredHz = 0.0;
            _phaseIntegral = 0.0;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public void UpdateFrequencyLock(
            double measuredIfFreqHz,
            double rfPowerDbm,
            out double currentModulationMa,
            out double pztDriveVoltageV,
            out double estimatedLinewidthKhz)
        {
            _currentIfFreqHz = measuredIfFreqHz;

            // 1. IF 주파수 오차 계산
            double freqError = _currentIfFreqHz - _targetIfFreqHz;

            // 2. 고주파 주파수 요동 1차 저역통과 필터링
            _freqErrorFilteredHz = 0.85 * _freqErrorFilteredHz + 0.15 * freqError;

            // 3. 고속 전류 서보: delta_I = - K_p * error (±15 mA 리미터)
            _cmdCurrentModulationMa = Math.Clamp(-_kpCurrent * _freqErrorFilteredHz, -15.0, 15.0);

            // 4. 저속 대변위 피에조 서보: V_pzt = - K_i * integral (±60 V 리미터)
            _phaseIntegral += freqError * _dt;
            _phaseIntegral = Math.Clamp(_phaseIntegral, -5e6, 5e6);
            _cmdPztVoltageV = Math.Clamp(-_kiPzt * _phaseIntegral, -60.0, 60.0);

            // 5. 레이저 유효 선폭 추정 (kHz)
            estimatedLinewidthKhz = Math.Max(0.5, Math.Abs(_freqErrorFilteredHz) * 1e-3 * 0.4);

            currentModulationMa = _cmdCurrentModulationMa;
            pztDriveVoltageV = _cmdPztVoltageV;
        }

        public readonly double CurrentIfFrequencyMhz => _currentIfFreqHz * 1e-6;
        public readonly double CurrentPztVoltage => _cmdPztVoltageV;
    }
}

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