AIDT 온디바이스 STEM 테라헤르츠 시분할 분광(THz-TDS) 광전도 안테나 서보 엔진

1. 개요

테라헤르츠(0.1 ~ 10 THz) 주파수 대역은 전자파와 광파의 경계에 위치하여 물질의 분자 간 진동, 거대 고분자 뒤틀림 모드, 반도체 자유 캐리어 동역학의 핑거프린트를 직접 검출할 수 있는 차세대 분광 영역입니다. 테라헤르츠 시분할 분광학(Terahertz Time-Domain Spectroscopy, THz-TDS)은 저온 성장 갈륨비소(LT-GaAs) 광전도 안테나(Photoconductive Antenna, PCA)에 펨토초 레이저를 조사하여 단일 사이클 피코초 THz 전계 펄스를 발생시키고, 전기광학 샘플링(ZnTe 결정)을 통해 진폭과 위상을 동시 계측합니다. 이를 통해 크라머스-크로니히(Kramers-Kronig) 관계식에 의존하지 않고 시료의 복소 굴절률(\(\tilde{n} = n + i \kappa\))을 절대적으로 산출합니다. 그러나 레이저 타이밍 지터와 안테나 열잡음으로 인한 시간축 표류는 위상 오차를 초래합니다. 본 고에서는 광학 딜레이 라인 10,000Hz 위상 지터 보정 C# 서보 엔진을 구현합니다.

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

광전도 안테나에서 방출되는 THz 전계 \(E_{THz}(t)\)는 안테나 과도 광전류 \(J(t)\)의 시간 미분에 비례합니다:

E_{THz}(t) \propto \frac{d J(t)}{dt} = e \frac{d}{dt} \left[ n_c(t) v_d(t) \right]

여기서 \(n_c(t)\)는 광유도 캐리어 밀도, \(v_d(t)\)는 바이어스 전계 하의 전자 표류 속도입니다. 두께 \(d\)인 시료를 투과한 THz 신호와 기준 신호의 푸리에 변환 복소 전달 함수 \(H(\omega)\)는 다음과 같습니다:

H(\omega) = \frac{E_{sample}(\omega)}{E_{ref}(\omega)} = \frac{4 \tilde{n}}{(\tilde{n} + 1)^2} \exp\left( -i \frac{(\tilde{n} - 1)\omega d}{c} \right)

펄스 피크 도착 시간 \(t_{peak}\)의 지터를 10,000Hz 보이스 코일 정밀 광 지연선 피드백 루프로 서브 펨토초 단위에서 잠금(Lock-in)하여 THz 위상 분해능을 보장합니다.

3. 실시간 피드백 제어 아키텍처

+-----------------------------------------------------------------------------+
|    AIDT THz-TDS Time-Domain Spectroscopy 10,000Hz Real-Time Servo Engine    |
+-----------------------------------------------------------------------------+
|  [ZnTe Balanced Photodiode Electro-Optic Detector & Fast Lock-In Amplifier] |
|                                       |                                     |
|                                       v                                     |
|     [THz Pulse Peak Time t_{peak} & Complex Permittivity (\epsilon_r, \epsilon_i)]|
|                                       |                                     |
|                                       v                                     |
|    [Zero-Alloc Discrete Optical Delay Line Sub-Femtosecond Tracking Servo]  |
|                                       |                                     |
|                                       v                                     |
|   [Sub-Nanometer Voice Coil Stage Actuator & PCA High-Voltage Bias Driver]  |
+-----------------------------------------------------------------------------+

4. Zero-Allocation C# 서보 엔진 구현

THz 펄스 피크 위상 드리프트를 실시간 보정하여 복소 광학 파라미터를 추출하는 C# 클래스입니다.

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Photonics.ThzTds
{
    public readonly struct ThzTdsConfig
    {
        public readonly double TargetPeakTimePs;
        public readonly double SampleThicknessUm;
        public readonly double OpticalC;
        public readonly double Kp;
        public readonly double Ki;

        public ThzTdsConfig(double targetPeakTimePs, double sampleThicknessUm,
                            double opticalC, double kp, double ki)
        {
            TargetPeakTimePs = targetPeakTimePs;
            SampleThicknessUm = sampleThicknessUm;
            OpticalC = opticalC;
            Kp = kp;
            Ki = ki;
        }
    }

    public ref struct ThzTdsServo
    {
        private readonly ThzTdsConfig _config;
        private double _timeJitterIntegral;
        private double _controlledDelayUm;

        public ThzTdsServo(in ThzTdsConfig config)
        {
            _config = config;
            _timeJitterIntegral = 0.0;
            _controlledDelayUm = 0.0;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public (double delayCorrectionUm, double realRefractiveIndex) Update(double measuredPeakTimePs, double dtSec)
        {
            double jitterErrorPs = measuredPeakTimePs - _config.TargetPeakTimePs;
            _timeJitterIntegral += jitterErrorPs * dtSec;
            _timeJitterIntegral = Math.Clamp(_timeJitterIntegral, -5.0, 5.0);

            // 시간 지터 상쇄 광 지연선 거리 보정 (\Delta L = c * \Delta t / 2)
            double deltaUm = (_config.Kp * jitterErrorPs) + (_config.Ki * _timeJitterIntegral);
            _controlledDelayUm = Math.Clamp(_controlledDelayUm + deltaUm, -500.0, 500.0);

            // 실시간 굴절률 n 추정치 계산
            double realN = 1.0 + (measuredPeakTimePs * 1e-12 * 3e8) / (_config.SampleThicknessUm * 1e-6);

            return (_controlledDelayUm, realN);
        }
    }
}

5. 학술적 의의 및 교육적 활용 가치

본 설계는 2022 개정 물리 및 화학 교과 내 '전자기파 스펙트럼과 파동', '물질의 유전 분극과 복소 굴절률', '분자 진동과 분광학'을 실시간 데이터로 탐구하는 첨단 STEM 교구입니다. 빛의 파동성과 전자기학적 유전 분극 메커니즘을 온디바이스 상에서 정밀 수치 제어로 시뮬레이션할 수 있습니다.