2022 개정 교육과정 연계 온디바이스 테라헤르츠물리·차세대스핀트로닉스 STEM 시리즈 #70

반강자성체(AFM) THz 스핀토크 발진기 서보 엔진

서브피코초 닐(Néel) 벡터 세차운동 및 0-Alloc C# 0.5THz 스핀 궤도 토크(SOT) 위상 잠금(PLL) 제어 루프

ZeliDesk 양자자성·테라헤르츠통신팀 초고속스핀소자·6G테라헤르츠 분과

2026년 9월 12일

1. 반강자성 닐 벡터 동역학 및 확장된 LLG 방정식

두 서브래티스 자기 모멘트 $\mathbf{m}_1, \mathbf{m}_2$를 갖는 반강자성체의 정규화된 닐 벡터 $\mathbf{n} = (\mathbf{m}_1 - \mathbf{m}_2)/2$와 총 자화 벡터 $\mathbf{m} = (\mathbf{m}_1 + \mathbf{m}_2)/2$에 대한 확장된 란다우-립시츠-길버트(LLG) 방정식은 다음과 같습니다:

교환 자기장($H_{ex} \sim 1000\text{ T}$)과 자기 이방성장($H_{an} \sim 1\text{ T}$)의 결합으로 유효 공진 각주파수는 테라헤르츠 영역으로 도약합니다:

중금속(Pt/W) 층에서 주입되는 스핀 홀 각도 $\theta_{SH}$ 매개 스핀 전류 $J_s$에 의해 유도되는 스핀 궤도 토크 $\mathbf{T}_{SOT}$는 닐 벡터를 연속 세차운동시켜 테라헤르츠 교류 전압을 출력합니다:

2. 0-Alloc C# 스택 기반 AFM THz 발진기 PLL 서보

테라헤르츠 다운컨버전 IF 신호로부터 순간 주파수와 위상 오차를 산출하여 주입 전류를 고속 변조하는 ref struct 엔진입니다.

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.Spintronics
{
    public readonly ref struct AfmStoTelemetry
    {
        public readonly double OutputFrequencyThz;
        public readonly double OutputPowerMicroWatts;
        public readonly double PhaseJitterRmsRad;
        public readonly double DeviceTemperatureKelvin;
        public readonly double BiasCurrentDensityMaCm2;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public AfmStoTelemetry(double freqThz, double pUf, double jitterRad, double tempK, double jBias)
        {
            OutputFrequencyThz = freqThz;
            OutputPowerMicroWatts = pUf;
            PhaseJitterRmsRad = jitterRad;
            DeviceTemperatureKelvin = tempK;
            BiasCurrentDensityMaCm2 = jBias;
        }
    }

    public readonly ref struct StoPllControlResult
    {
        public readonly double FrequencyErrorGhz;
        public readonly double BiasCurrentDeltaMaCm2;
        public readonly double PhaseNoiseReductionDb;
        public readonly bool IsFrequencyLocked;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public StoPllControlResult(double freqErrGhz, double jDelta, double snrImpDb, bool isLocked)
        {
            FrequencyErrorGhz = freqErrGhz;
            BiasCurrentDeltaMaCm2 = jDelta;
            PhaseNoiseReductionDb = snrImpDb;
            IsFrequencyLocked = isLocked;
        }
    }

    public ref struct AntiferromagneticPllServo
    {
        private readonly double _targetFrequencyThz;
        private readonly double _kpCurrent;
        private readonly double _kiCurrent;
        private double _integratedFreqError;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public AntiferromagneticPllServo(double targetThz, double kp, double ki)
        {
            _targetFrequencyThz = targetThz;
            _kpCurrent = kp;
            _kiCurrent = ki;
            _integratedFreqError = 0.0;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public StoPllControlResult Step(in AfmStoTelemetry telemetry, double dt)
        {
            // 1. 목표 주파수(예: 0.500 THz = 500 GHz) 대비 오차 산출 (GHz 단위)
            double freqErrorGhz = (telemetry.OutputFrequencyThz - _targetFrequencyThz) * 1000.0;

            // 2. 오차 적분 및 PI 바이어스 전류 튜닝 (주입 전류 증가 -> 유효 감쇠 보상 및 주파수 블루시프트)
            _integratedFreqError += freqErrorGhz * dt;
            double currentDelta = -(_kpCurrent * freqErrorGhz + _kiCurrent * _integratedFreqError);

            // 3. 온도 변화에 따른 스핀 홀 전도도 열 드리프트 보정 (T = 300K 기준)
            double tempFactor = 1.0 + 0.0015 * (telemetry.DeviceTemperatureKelvin - 300.0);
            currentDelta *= tempFactor;

            // 4. 위상 잠금(PLL Lock) 상태 판별: 주파수 오차 0.5 GHz 미만 & 위상 지터 0.1 rad 미만
            bool isLocked = Math.Abs(freqErrorGhz) < 0.5 && telemetry.PhaseJitterRmsRad < 0.10;

            // 5. 위상 노이즈 억제율 추정
            double noiseReductionDb = isLocked ? 18.5 - 20.0 * telemetry.PhaseJitterRmsRad : 3.0;

            return new StoPllControlResult(freqErrorGhz, currentDelta, noiseReductionDb, isLocked);
        }
    }
}

3. 교육적 가치 및 중고교 실험 시뮬레이션 적용