AIDT 온디바이스 STEM 상대론적 레이저 후류장 가속(LWFA) 플라즈마 기포 서보 엔진

1. 개요: 수 킬로미터 거대 가속기를 테이블탑(Table-top) 수 센티미터 칩 크기로 압축하다

전통적인 고주파(RF) 공동 선형 입자 가속기(예: SLAC)는 금속 표면의 방전 파괴(Dielectric Breakdown) 한계로 인해 가속 전기장이 수십 \(\text{MV/m}\)로 제한되며, GeV~TeV 고에너지 전자를 얻기 위해 수 킬로미터의 거대한 시설이 필요합니다. 이미 전리된 매질인 플라즈마는 물질 파괴 한계가 존재하지 않습니다.

레이저 후류장 가속(Laser Wakefield Acceleration, LWFA)은 테라와트(TW)에서 페타와트(PW)급 초단 레이저 펄스를 초음속 가스 젯에 집속시켜 발생시킵니다. 강력한 폰더로모티브 힘(Ponderomotive Force)이 플라즈마 전자를 레이저 경로 밖으로 밀어내고 무거운 이온만 남겨놓아 레이저 뒤편에 완벽한 구형 이온 공동인 블로아웃 기포(Bubble / Blowout Regime)를 형성합니다. 이 기포 후면에 형성되는 정전기 가속 전기장은 \(100\,\text{GV/m}\)에 달하여 불과 수 센티미터 내에서 전자를 GeV 에너지로 가속합니다.

2. 폰더로모티브 포텐셜 및 블로아웃 가속 전기장 수식 모델

레이저의 정규화된 벡터 포텐셜 \(a_0\)는 전기장 피크 \(E_0\)와 레이저 주파수 \(\omega_0\)에 의해 다음과 같이 정의됩니다:

\\[ a_0 = \frac{e E_0}{m_e \omega_0 c} \approx 0.85 \times 10^{-9} \lambda_0 [\mu\text{m}] \sqrt{I_0 [\text{W/cm}^2]} \\]

\(a_0 > 2\)인 초상대론적(Ultra-relativistic) 강도 하에서 전자에 작용하는 상대론적 폰더로모티브 힘 \(\mathbf{F}_p\)는 다음과 같습니다:

\\[ \mathbf{F}_p = -m_e c^2 \nabla \sqrt{1 + \frac{\langle a^2 \rangle}{2}} \\]

플라즈마 파수 \(k_p = \frac{\omega_p}{c} = \sqrt{\frac{n_e e^2}{\epsilon_0 m_e c^2}}\)에 대해 플라즈마 파장 \(\lambda_p = \frac{2\pi}{k_p}\) 수준의 크기를 갖는 기포(Bubble) 내부에서 전자가 느끼는 최대 가속 전기장 \(E_{\text{acc}}\)과 에너지 이득 한계(디페이징 길이 \(L_{\text{deph}}\) 기준)는 다음과 같습니다:

\\[ E_{\text{acc}} \approx E_0 \sqrt{a_0} = \frac{m_e c \omega_p}{e} \sqrt{a_0} \approx 96 [\text{GV/m}] \sqrt{a_0} \sqrt{\frac{n_e [10^{18}\text{cm}^{-3}]}{1}} \\]

\\[ \Delta E_{\text{max}} \approx e E_{\text{acc}} L_{\text{deph}} \approx \frac{2}{3} m_e c^2 a_0 \left( \frac{\omega_0}{\omega_p} \right)^2 \\]

가속된 전자가 레이저 후류 위상보다 빨라져 감속 영역으로 진입하는 디페이징(Dephasing) 현상을 방지하기 위해 가스 노즐의 밀도 테이퍼링(Density Tapering)을 정밀 제어해야 합니다.

3. 실시간 LWFA 전자 빔 가속 제어 서보 아키텍처 다이어그램

graph TD A[100TW 티타늄 사파이어 펨토초 레이저] --> B[축외 포물면 거울 OAP 집속 a_0 > 3] C[초음속 He/N2 가스 노즐 밸브] --> D[플라즈마 밀도 n_e ~ 2x10^18 cm^-3 형성] B & D --> E[블로아웃 기포 형성 및 이온화 주입 Ionization Injection] E --> F[수 센티미터 가속 거리 내 100 GV/m 초강력 가속] F --> G[쌍극자 자석 스펙트로미터 & Lanex 형광판 계측] G --> H[전자 빔 에너지 E_beam & 에너지 퍼짐 Delta E / E 계측] H --> I[10,000Hz 0-Alloc C# 노즐 배압 및 레이저 초점 서보] I --> C

4. 고성능 0-Alloc C# 레이저 후류장 가속 서보 구현

플라즈마 밀도와 레이저 강도에 따른 가속 전기장과 전자 에너지 이득을 실시간 산출하는 0-Alloc C# 클래스입니다:

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.AcceleratorPhysics
{
    public readonly record struct LwfaParameters(
        double LaserWavelengthUm, // um (0.8 um)
        double LaserPulseEnergyJ, // Joules (예: 3.0 J)
        double PulseDurationFs,   // fs (예: 30 fs)
        double BeamWaistW0Um,     // um (초점 반경 예: 15 um)
        double AccelerationLengthMm // mm (예: 10 mm = 1 cm)
    );

    public ref struct LwfaServoState
    {
        public double NormalizedA0;
        public double PlasmaDensityCm3;
        public double PeakAccelFieldGvM;
        public double DephasingLengthMm;
        public double FinalElectronEnergyGev;
        public double EnergySpreadPercent;
        public bool IsGeVReachable;
    }

    public static class LaserWakefieldServo
    {
        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public static void CalculateAcceleration(
            in LwfaParameters p,
            double backingPressureBar, // 가스 노즐 배압 (bar)
            double nozzleSlitWidthMm,
            ref LwfaServoState state)
        {
            // 1. 배압에 따른 플라즈마 전자 밀도 추정 (1 bar ~ 1.5e18 cm^-3)
            state.PlasmaDensityCm3 = Math.Max(1e17, backingPressureBar * 1.5e18);

            // 2. 레이저 피크 강도 및 a_0 계산
            double peakPowerWatts = p.LaserPulseEnergyJ / (p.PulseDurationFs * 1e-15);
            double waistCm = p.BeamWaistW0Um * 1e-4;
            double areaCm2 = 0.5 * Math.PI * waistCm * waistCm;
            double intensityWcm2 = peakPowerWatts / areaCm2;

            // a_0 = 0.85e-9 * lambda_um * sqrt(I)
            state.NormalizedA0 = 0.85e-9 * p.LaserWavelengthUm * Math.Sqrt(intensityWcm2);

            // 3. 플라즈마 파수 및 파장
            // omega_p = sqrt(n_e * e^2 / (eps_0 * m_e))
            // lambda_p [um] ~ 33.4 / sqrt(n_e [10^18 cm^-3])
            double ne18 = state.PlasmaDensityCm3 / 1e18;
            double lambdaPUm = 33.4 / Math.Sqrt(ne18);

            // 4. 가속 전기장 E_acc ~ 96 [GV/m] * sqrt(a_0 * n_e18)
            state.PeakAccelFieldGvM = 96.0 * Math.Sqrt(state.NormalizedA0 * ne18);

            // 5. 디페이징 길이 L_deph ~ (2/3) * (omega_0 / omega_p)^2 * lambda_p * sqrt(a_0)
            double freqRatio = lambdaPUm / p.LaserWavelengthUm; // omega_0 / omega_p
            double lDephUm = (2.0 / 3.0) * Math.Pow(freqRatio, 2) * lambdaPUm * Math.Sqrt(state.NormalizedA0);
            state.DephasingLengthMm = lDephUm * 0.001;

            // 6. 유효 가속 거리 및 최종 전자 에너지: E = e * E_acc * min(L_acc, L_deph)
            double effectiveLMm = Math.Min(p.AccelerationLengthMm, state.DephasingLengthMm);
            double energyGainGev = (state.PeakAccelFieldGvM) * (effectiveLMm * 0.001); // GV/m * m = GV = GeV
            state.FinalElectronEnergyGev = energyGainGev;

            // 에너지 퍼짐 추정
            state.EnergySpreadPercent = Math.Clamp(1.5 / Math.Max(0.5, state.FinalElectronEnergyGev), 0.5, 8.0);
            state.IsGeVReachable = state.FinalElectronEnergyGev >= 1.0;
        }
    }
}

5. 교육 연수 및 탁상형 방사광/핵사찰 국방 적용 지침

방통대 직무연수 및 국방 초소형 X-선 자유전자레이저(XFEL)

물리 교과 연계: 특수 상대성 이론의 질량 증가, 전자기파 폰더로모티브 힘, 플라즈마 집단 진동과 파동 타기(Surf-riding) 가속 메커니즘을 3D 그래픽스로 교육.
국방 항만 컨테이너 투시 및 방사광 검사: 수 킬로미터 가속기 없이 트럭 탑재형 페타와트 레이저를 이용하여 은닉 핵물질 및 중금속 차폐 화물을 투과하는 고휘도 펨토초 감마선 빔 발생.