AIDT 온디바이스 STEM 광학 핀셋(Optical Tweezers) 방사압 및 광자 운동량 전달 서보 엔진

1. 개요: 빛으로 나노·마이크로 입자를 포획하고 피코뉴턴(\(pN\))의 힘을 계측하다

아서 아슈킨(Arthur Ashkin, 2018 노벨 물리학상)이 개발한 광학 핀셋(Optical Tweezers)은 고개구수(High-NA) 대물렌즈로 강하게 집속된 단일 레이저 빔을 이용하여 유전체 입자(비드, 박테리아, 바이러스, DNA 등)를 비침습적으로 포집하고 3차원 위치를 나노미터 정밀도로 제어하는 기술입니다.

빛이 입자를 통과하거나 반사될 때 발생하는 광자의 운동량 변화(\(\\Delta \\mathbf{p}\))는 뉴턴 제3법칙에 의해 입자에 반작용 방사압(Radiation Pressure)을 가합니다. 이 힘은 빔의 진행 방향으로 밀어내는 산란력(Scattering Force)과 광 강도가 가장 높은 초점 중심으로 끌어당기는 구배력(Gradient Force)으로 분해되며, 구배력이 산란력을 압도할 때 3차원 안정 포집(Stable Optical Trapping)이 달성됩니다.

2. 레일리 및 기하광학 정권에서의 광학력 수식 모델

입자 반경 \(a\)가 레이저 파장 \(\\lambda\)보다 훨씬 작은 레일리 정권(Rayleigh Regime, \(a \\ll \\lambda\))에서 유전체 입자는 점 유도 쌍극자(Induced Dipole)로 취급됩니다. 입자의 분극률(Polarizability) \(\\alpha\)는 클라우지우스-모소티 관계식에 의해 다음과 같습니다:

\\[ \\alpha = 4 \\pi \\epsilon_m \\epsilon_0 a^3 \\left( \\frac{m^2 - 1}{m^2 + 2} \\right) \\]

여기서 \(m = n_p / n_m\)은 매질 대비 입자의 상대 굴절률입니다. 초점 중심으로 끌어당기는 광학 구배력 \(\\mathbf{F}_{\\text{grad}}\)와 빔 전파 축으로 밀어내는 산란력 \(\\mathbf{F}_{\\text{scat}}\)는 각각 다음과 같습니다:

\\[ \\mathbf{F}_{\\text{grad}} = \\frac{1}{2} \\alpha \\nabla \\langle E^2 \\rangle = \\frac{2 \\pi n_m a^3}{c} \\left( \\frac{m^2 - 1}{m^2 + 2} \\right) \\nabla I(\\mathbf{r}) \\]

\\[ \\mathbf{F}_{\\text{scat}} = \\frac{n_m \\sigma_{\\text{scat}}}{c} I(\\mathbf{r}) \\hat{\\mathbf{z}}, \\quad \\sigma_{\\text{scat}} = \\frac{8 \\pi}{3} k^4 a^6 \\left( \\frac{m^2 - 1}{m^2 + 2} \\right)^2 \\]

축 방향(\(z\)축) 안정 포집 조건은 초점 직전에서 역방향 구배력의 최대값이 순방향 산란력보다 커야 함을 의미합니다:

\\[ \\left| \\frac{F_{z, \\text{grad}}}{F_{z, \\text{scat}}} \\right| > 1 \\implies \\text{NA} > 1.2 \\quad (\\text{고개구수 침액 렌즈 필수}) \\]

3. 광학 핀셋 위치 피드백 및 강성 교정 아키텍처 다이어그램

graph TD A[Nd:YAG 1064nm 레이저 빔] --> B[AOD 음향광학 편향기 / 공간광변조기 SLM] B --> C[고개구수 대물렌즈 NA > 1.3 초점 집속] C --> D[마이크로 비드 3차원 광학 포집] D --> E[포워드 산란광 수집 콘덴서 렌즈] E --> F[4분할 포토다이오드 QPD 간섭 신호 계측] F --> G[후방 초점면 간섭계 BFP 전압 신호 Delta V_x, Delta V_y] G --> H[10,000Hz 0-Alloc C# 나노미터 위치 & 피코뉴턴 힘 서보] H --> B

4. 고성능 0-Alloc C# 광학 핀셋 서보 엔진 구현

QPD 신호로부터 변위와 힘을 계산하고, AOD 피드백을 통해 트랩 강성(\(\\kappa\))과 위치를 실시간 안정화하는 C# 구조체입니다:

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.OpticalTweezers
{
    public readonly record struct TrapParameters(
        double LaserPower,        // Watts (예: 50 mW = 0.05)
        double Wavelength,        // m (1064 nm = 1064e-9)
        double ParticleRadius,    // m (예: 500 nm = 500e-9)
        double MediumIndex,       // n_m = 1.33 (물)
        double ParticleIndex,     // n_p = 1.59 (폴리스티렌 비드)
        double ViscosityEta       // Pa*s (0.001 물)
    );

    public ref struct TrapServoState
    {
        public double TrapStiffnessX; // N/m (pN/nm)
        public double TrapStiffnessZ;
        public double DisplacementX;  // m
        public double DisplacementZ;
        public double RestoringForceX;// N
        public double RestoringForceZ;
        public double EquipartitionThermalNoise; // m RMS
        public bool IsTrapped;
    }

    public static class OpticalTweezersServo
    {
        private const double SpeedOfLight = 299792458.0;
        private const double BoltzmannK = 1.380649e-23; // J/K
        private const double TemperatureK = 298.15;     // 25 C

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public static void StepServo(
            in TrapParameters p,
            double measuredQpdVoltX, // V
            double measuredQpdVoltZ,
            double calibFactorX,     // m/V
            double calibFactorZ,
            ref TrapServoState state)
        {
            // 1. 입자 변위 역산
            state.DisplacementX = measuredQpdVoltX * calibFactorX;
            state.DisplacementZ = measuredQpdVoltZ * calibFactorZ;

            // 2. 상대 굴절률 m 및 분극률 극성 인자
            double m = p.ParticleIndex / p.MediumIndex;
            double m2 = m * m;
            double polarFactor = (m2 - 1.0) / (m2 + 2.0);

            // 3. 트랩 강성 근사 (파워 선형 비례 모델)
            // kappa_x ~ 4 * pi * n_m * a^3 / c * polarFactor * (I_0 / w_0^2)
            double stiffnessScale = (p.LaserPower * polarFactor * Math.Pow(p.ParticleRadius / 500e-9, 3)) * 1e-4;
            state.TrapStiffnessX = Math.Max(1e-6, stiffnessScale);       // 예: 50 pN/um = 5e-5 N/m
            state.TrapStiffnessZ = Math.Max(1e-7, stiffnessScale * 0.25); // 축 방향 강성은 횡방향의 ~25%

            // 4. 복원력 계산: F = -kappa * x
            state.RestoringForceX = -state.TrapStiffnessX * state.DisplacementX;
            state.RestoringForceZ = -state.TrapStiffnessZ * state.DisplacementZ;

            // 5. 열 요동(브라운 운동) 크기:  = k_B * T / kappa
            state.EquipartitionThermalNoise = Math.Sqrt((BoltzmannK * TemperatureK) / state.TrapStiffnessX);

            // 6. 포집 안정성 판정: 복원 위치가 빔 웨이스트 반경 이내인가
            double rDisp = Math.Sqrt(state.DisplacementX * state.DisplacementX + state.DisplacementZ * state.DisplacementZ);
            state.IsTrapped = rDisp < (p.ParticleRadius * 3.0);
        }
    }
}

5. 교육 연수 및 방위 바이오 감시 적용 지침

방통대 직무연수 및 국방 생화학무기(CBRN) 나노 분석

물리 교과 연계: 광자의 운동량(\(p = h/\\lambda\)), 광선 굴절에 의한 모멘텀 역작용, 하모닉 오실레이터 복원력과 열 요동의 등분배 정리를 가상 실험실 교구로 구성.
국방 바이오 방어 체계: 탄저균, 두창 등 고위험성 에어로졸 병원균의 단일 포자 수준 비접촉 포집 및 형광/라만 연계 초신속 동정.