AIDT 온디바이스 STEM 광자기 효과(Faraday & Voigt) 테라헤르츠 편광 회전 서보 엔진

1. 개요: 자기장에 의해 빛의 편광면이 비가역적으로 회전하는 광자기 상호작용

일반적인 광학 매질에서 반사된 빛은 동일한 경로를 되돌아오며 편광 회전이 상쇄(가역적)되지만, 매질에 외부 정자기장 \(\mathbf{B}\)가 인가되면 시간 역전 대칭성(Time-reversal symmetry)이 깨지면서 광자기 효과(Magneto-optical Effect)가 발생합니다. 자기장 방향과 나란하게 전파하는 빛의 선편광면이 회전하는 현상을 패러데이 효과(Faraday Effect)라 하며, 이는 좌원편광(LCP)과 우원편광(RCP)의 굴절률(\(n_+, n_-\))이 달라지는 원편광 복굴절(Circular Birefringence)에 기인합니다.

패러데이 효과의 비가역적 특성은 레이저 발진기에서 후방 반사광이 이득 매질을 파괴하지 못하도록 막아주는 광 아이솔레이터(Optical Isolator)의 기본 작동 원리이며, 고감도 자력계 및 테라헤르츠 위상 변조기에도 널리 활용됩니다. 본 백서에서는 테르븀 갈륨 가넷(TGG) 결정 및 페라이트 매질의 베르데 상수를 모델링하고 45도 회전각을 일정하게 유지하는 10,000Hz 0-Alloc C# 전자기 서보를 구현합니다.

2. 유전율 텐서 및 패러데이 회전각 수식 모델

\(z\)축 방향 정자기장 \(B_z\)가 인가된 자이로트로픽(Gyrotropic) 매질의 복소 유전율 텐서 \(\boldsymbol{\epsilon}\)는 비대각 성분(Off-diagonal components)을 갖습니다:

\\[ \\boldsymbol{\epsilon} = \\begin{pmatrix} \\epsilon_{xx} & -i\\epsilon_{xy} & 0 \\\\ i\\epsilon_{xy} & \\epsilon_{xx} & 0 \\\\ 0 & 0 & \\epsilon_{zz} \\end{pmatrix} \\]

이때 좌우 원편광의 유효 굴절률 \(n_\\pm\)은 다음과 같이 분열됩니다:

\\[ n_\\pm = \\sqrt{\\epsilon_{xx} \\pm \\epsilon_{xy}} \\approx n_0 \\pm \\frac{\\epsilon_{xy}}{2 n_0} \\]

길이 \(L\)의 매질을 통과한 선편광의 패러데이 회전각 \(\\theta_F\)는 두 원편광 모드의 위상차로부터 도출되며, 베르데 상수(Verdet Constant, \(V\))와 자기장 \(B\)의 곱으로 선형 근사됩니다:

\\[ \\theta_F = \\frac{\\omega L}{2 c} (n_+ - n_-) = V(\\lambda, T) \\cdot B \\cdot L \\]

비가역성(Non-reciprocity)으로 인해 왕복 경로(\(2L\))를 진행할 경우 회전각은 상쇄되지 않고 \(2\\theta_F\)로 누적됩니다. 한편 자기장에 수직인 방향으로 전파하는 빛의 선형 복굴절 현상인 보이트 효과(Voigt Effect / Cotton-Mouton Effect)의 위상 지연 \(\\delta_V\)는 자기장의 제곱에 비례합니다:

\\[ \\delta_V = C_M \\cdot B^2 \\cdot L \\]

3. 실시간 광 아이솔레이터 45도 회전각 서보 아키텍처 다이어그램

graph TD A[입사 선편광 레이저 theta = 0 deg] --> B[TGG 광자기 결정 소자 L] C[영구자석 & 보정 전자석 코일 B_z] --> B B --> D[패러데이 편광 회전 theta_F = V * B * L] D --> E[편광 빔 분할기 PBS 및 밸런스드 포토다이오드] E --> F[회전각 편차 오차 신호 Delta theta = theta_F - 45 deg] F --> G[10,000Hz 0-Alloc C# 코일 전류 제어기] G --> C D --> H[45도 편광자 투과 / 후방 반사광 완벽 90도 차단]

4. 고성능 0-Alloc C# 패러데이 회전각 서보 구현

온도 변화에 따른 베르데 상수 변동을 코일 보정 자기장으로 실시간 상쇄하여 45도 정밀도를 유지하는 0-Alloc C# 클래스입니다:

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.MagnetoOptics
{
    public readonly record struct FaradayParameters(
        double CrystalLengthL,     // m (예: 20 mm = 0.02)
        double BaseVerdetRadTm,    // rad / (T * m) (TGG 1064nm 기준 ~ -40 rad/(T*m))
        double VerdetTempCoeff,    // 1/K 베르데 상수 온도 계수 (예: -0.003 / K)
        double TargetAngleRad      // 목표 편광 회전각 (45도 = pi / 4)
    );

    public ref struct FaradayServoState
    {
        public double CurrentBField;
        public double MeasuredAngleRad;
        public double AngleErrorDeg;
        public double CorrectionCurrentAmp;
        public double OpticalIsolationDb;
        public bool IsLocked;
    }

    public static class FaradayRotatorServo
    {
        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public static void StepServo(
            in FaradayParameters p,
            double temperatureK,
            double nominalBField,  // 기본 영구자석 자기장 (Tesla)
            double dt,
            ref FaradayServoState state)
        {
            // 1. 온도에 따른 베르데 상수 보정: V(T) = V_0 * (1 + beta * (T - 295))
            double deltaT = temperatureK - 295.15;
            double actualVerdet = p.BaseVerdetRadTm * (1.0 + p.VerdetTempCoeff * deltaT);

            // 2. 현재 자기장 하에서의 회전각: theta = V * B * L
            double effectiveB = nominalBField + (state.CorrectionCurrentAmp * 0.05); // 코일 상수 0.05 T/A
            state.CurrentBField = effectiveB;
            state.MeasuredAngleRad = Math.Abs(actualVerdet * effectiveB * p.CrystalLengthL);

            // 3. 목표 45도와의 각도 오차
            double errorRad = p.TargetAngleRad - state.MeasuredAngleRad;
            state.AngleErrorDeg = errorRad * (180.0 / Math.PI);

            // 4. 비례-적분(PI) 코일 전류 보정
            double reqDeltaB = errorRad / (actualVerdet * p.CrystalLengthL);
            state.CorrectionCurrentAmp += reqDeltaB * 20.0 * dt; // 수렴 제어

            // 5. 광 격리도(Optical Isolation) 계산: -10 * log10(sin^2(Delta theta))
            double sinErr = Math.Sin(errorRad);
            double leakage = Math.Max(1e-6, sinErr * sinErr);
            state.OpticalIsolationDb = -10.0 * Math.Log10(leakage);

            state.IsLocked = Math.Abs(state.AngleErrorDeg) < 0.05; // 0.05도 이내 제어
        }
    }
}

5. 교육 연수 및 레이저 무기/전자전 적용 지침

방통대 직무연수 및 고출력 레이저 방호

물리/정보 교과 지도: 전자기파의 편광 벡터, 대각화 유전율 텐서, 베르데 상수 온도 보정 및 비가역 광소자 원리를 컴퓨터 그래픽스 실험으로 이해.
국방 고출력 레이저 요격기(HEL): 수십 킬로와트급 광섬유 레이저 증폭기 간 상호 발진 간섭 차단 및 테라헤르츠 고속 전자기 펄스(EMP) 편광 방어막 구축.