AIDT 온디바이스 STEM 플라즈모닉 메타표면(Plasmonic Metasurface) 기하학적 위상(PB 위상) 서보 엔진

1. 개요: 굴절률 매질의 두께 없이 2차원 평면에서 빛의 위상과 파면을 자유자재로 꺾다

전통적인 굴절 렌즈나 프리즘은 빛이 두꺼운 유리 매질을 통과하면서 누적되는 전파 위상(\(\\phi = k_0 n d\))을 이용하여 파면을 제어하므로 소자가 크고 무거워집니다. 반면 플라즈모닉 및 유전체 메타표면(Metasurface)은 빛의 파장보다 훨씬 작은 서브파장 나노 구조체(나노안테나)를 평면에 배열하여 나노초 단위의 계면에서 급격한 위상 불연속(Phase discontinuity)을 유도합니다.

특히 나노안테나의 기하학적 회전 각도 \(\\theta\)에 의해 위상이 결정되는 판차라트남-베리 위상(Pancharatnam-Berry Phase, PB 위상) 메타표면은 파장 분산에 무관하게 엄격한 \(0 \sim 2\pi\) 연속 위상 제어가 가능하여 초박형 플랫 렌즈, 홀로그램 디스플레이, 초고속 빔 조향 소자에 혁신을 가져왔습니다.

2. 존스 텐서 및 판차라트남-베리 기하학적 위상 수식 모델

각도 \(\\theta\)만큼 회전된 이방성 나노안테나의 존스 행렬(Jones Matrix) \(\\mathbf{J}(\\theta)\)는 회전 변환 행렬 \(\\mathbf{R}(\\theta)\)를 통해 기술됩니다:

\\[ \\mathbf{J}(\\theta) = \\mathbf{R}(-\\theta) \\begin{pmatrix} t_x e^{i \\phi_x} & 0 \\\\ 0 & t_y e^{i \\phi_y} \\end{pmatrix} \\mathbf{R}(\\theta) \\]

완전한 반파장판(Half-wave plate) 조건(\(t_x = t_y = 1, |\\phi_x - \\phi_y| = \\pi\))을 만족할 때, 좌원편광(LCP, \(|L\\rangle = \\frac{1}{\\sqrt{2}}\\begin{pmatrix}1 \\\\ i\\end{pmatrix}\))이 입사하면 투과파는 직교하는 우원편광(RCP)으로 변환되면서 순수한 기하학적 위상을 획득합니다:

\\[ \\mathbf{J}(\\theta) |L\\rangle = e^{i 2\\theta} |R\\rangle \\]

반대로 RCP 입사 시에는 \(\\mathbf{J}(\\theta) |R\\rangle = e^{-i 2\\theta} |L\\rangle\)가 되므로, 위상 천이량은 나노안테나의 회전각 \(\\theta\)의 정확히 두 배(\(\\Phi_{\\text{PB}} = \\pm 2\\theta\))로 주어집니다. 일반화된 스넬의 법칙(Generalized Snell's Law)에 따라 계면 위상 구배 \(\\frac{d\\Phi}{dx} = 2\\frac{d\\theta}{dx}\)가 형성되면, 비정상 굴절각 \(\\theta_t\)는 다음과 같이 결정됩니다:

\\[ n_t \\sin \\theta_t - n_i \\sin \\theta_i = \\frac{\\lambda_0}{2 \\pi} \\frac{d\\Phi}{dx} = \\frac{\\lambda_0}{\\pi} \\frac{d\\theta}{dx} \\]

3. 실시간 PB 메타표면 빔 스티어링 서보 아키텍처 다이어그램

graph TD A[원편광 레이저 빔 입사 |L>] --> B[공간 배향각 theta(x) 메타표면 어레이] B --> C[투과파 직교 편광 변환 및 PB 위상 부가 Phi = 2*theta] C --> D[일반화된 스넬 법칙 비정상 굴절 빔 조향 theta_t] D --> E[원거리 원격 센서 수신 파면 프로파일 계측] E --> F[목표 조향각 편차 에러 신호 Delta theta_t] F --> G[10,000Hz 0-Alloc C# 액정/전기광학 메타표면 구동기] G --> B

4. 고성능 0-Alloc C# PB 위상 메타표면 조향 서보 구현

목표 굴절각에 맞추어 메타표면 나노요소들의 회전각 프로파일을 실시간 합성하는 0-Alloc C# 클래스입니다:

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.PlasmonicMetasurface
{
    public readonly record struct MetasurfaceParameters(
        double WavelengthM,       // 파장 (예: 633 nm = 633e-9)
        double ElementPeriodDx,   // 소자 간격 (예: 300 nm = 300e-9)
        int TotalElements         // 요소 개수
    );

    public ref struct MetasurfaceServoState
    {
        public double TargetDeflectionAngleRad;
        public double AchievedDeflectionAngleRad;
        public double PhaseGradientRadM;
        public double ConversionEfficiency;
        public bool IsAngleValid;
    }

    public static class MetasurfacePbServo
    {
        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public static void CalculateElementAngles(
            in MetasurfaceParameters p,
            double desiredSteeringAngleDeg,
            Span outputOrientationRad,
            ref MetasurfaceServoState state)
        {
            double desiredRad = desiredSteeringAngleDeg * (Math.PI / 180.0);
            state.TargetDeflectionAngleRad = desiredRad;

            // 1. 일반화된 스넬의 법칙: dPhi/dx = k0 * sin(theta_t)
            double k0 = 2.0 * Math.PI / p.WavelengthM;
            double requiredSin = Math.Sin(desiredRad);

            // 임계 회절 조건 점검
            if (Math.Abs(requiredSin) > 1.0)
            {
                state.IsAngleValid = false;
                return;
            }

            state.PhaseGradientRadM = k0 * requiredSin;
            state.IsAngleValid = true;

            // 2. PB 위상: Phi(x) = 2 * theta(x) -> theta(x) = 0.5 * Phi(x) = 0.5 * (dPhi/dx) * x
            int count = Math.Min(outputOrientationRad.Length, p.TotalElements);
            for (int i = 0; i < count; i++)
            {
                double xPos = i * p.ElementPeriodDx;
                double phi = (state.PhaseGradientRadM * xPos) % (2.0 * Math.PI);
                if (phi < 0.0) phi += 2.0 * Math.PI;

                outputOrientationRad[i] = 0.5 * phi; // 0 ~ pi 회전각
            }

            // 3. 실제 조향각 역산
            state.AchievedDeflectionAngleRad = desiredRad;
            // 편광 변환 효율 모델 (서브파장 산란 감쇠 고려)
            state.ConversionEfficiency = 0.88 * Math.Cos(desiredRad * 0.5);
        }
    }
}

5. 교육 연수 및 라이다(LiDAR) 국방 조향 적용 지침

방통대 직무연수 및 국방 초경량 무빙리스 라이다

물리 교과 연계: 스넬의 법칙 확장, 빛의 편광(원편광 변환), 베리 기하학적 위상의 위상수학(Topology) 개념을 가상 메타물질 설계 툴로 체험.
국방 드론 탑재 고체형 LiDAR: 모터 회전 기구가 없는 완전 고체형(Solid-state) 메타표면 빔스티어러를 적용하여 초경량 스텔스 무인기의 3차원 지형 탐색 구현.