AIDT 온디바이스 STEM 원자간력 현미경(AFM) qPlus 센서 피에조 피드백 서보 제어기

1. 개요

표면 물리 및 분자 나노기술에서 단일 분자의 화학 결합 구조를 시각화하는 비접촉 원자간력 현미경(NC-AFM)은 현대 기초과학의 핵심 도구입니다. 특히 석영 소리굽쇠 기반의 qPlus 센서는 수천 N/m에 달하는 극히 높은 스프링 상수(\(k \approx 1800 \text{ N/m}\))와 높은 기계적 품질 계수(\(Q > 30,000\))를 지녀, 피코미터(\(10^{-12}\text{m}\)) 단위의 미세 진폭으로 탐침 팁을 진동시키면서 팁-시편 간 단일 화학 결합 인력을 비파괴적으로 측정할 수 있습니다. 주파수 변조 AFM(FM-AFM) 방식에서는 원자간력 기울기(Force Gradient)에 의한 공진 주파수 이동량(\(\Delta f\))을 계측하여 팁-시편 간격을 일정하게 유지해야 합니다. 그러나 온도 드리프트와 기계적 크립(Creep)에 의한 나노스케일 충돌(Crash)을 방지하기 위해서는 초고속 Z축 피에조 변위 서보 제어가 필수적입니다. 본 고에서는 Giessibl 섭동 공식을 통해 상호작용력을 실시간 역산하고 10,000Hz 주기로 Z축 피에조 변위를 능동 제어하는 0-Alloc C# 엔진을 구현합니다.

2. 수학적 모델링 및 지배 방정식

qPlus 캔틸레버의 기계적 고유진동수 \(f_0\)와 감쇠 상수 \(\gamma\)를 갖는 조화 진동자 방정식에 팁-시편 상호작용력 \(F_{ts}(z)\)가 가해질 때의 운동 방정식은 다음과 같습니다:

진동 진폭 \(A\)가 상호작용력 감쇠 길이보다 작은 작은 진폭 극한에서 Giessibl 섭동 이론에 따른 공진 주파수 이동량 \(\Delta f\)는 힘의 기울기(Force Gradient) \(\frac{\partial F_{ts}}{\partial z}\)에 직접 비례합니다:

단일 원자 수준의 Lennard-Jones 포텐셜 상호작용력 모델은 인력(van der Waals)과 척력(Pauli 배타 원리)의 합으로 기술됩니다:

Z축 피에조 피드백 서보의 제어 편차 \(e_f(t) = \Delta f(t) - \Delta f_{setpoint}\)에 대한 2차 저역통과 안티앨리어싱 및 PI 피드백 전압 제어 법칙은 다음과 같습니다:

3. 실시간 C# 서보 제어기 구현

다음은 PLL 주파수 시프트 신호로부터 원자간력 구배를 역산하고 10,000Hz로 Z축 피에조 드라이버를 구동하는 0-Alloc C# 엔진 구현체입니다:

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Core.Nanotech
{
    public ref struct AfmQplusPiezoServo
    {
        private readonly double _dt;
        private readonly double _f0; // 고유 진동수 (약 30 kHz)
        private readonly double _springConstantK; // N/m (약 1800)
        private readonly double _piezoSensitivityNmPerV;
        private double _freqSetpointHz;
        
        // Z축 상태 변수
        private double _currentDeltaF;
        private double _filteredDeltaF;
        private double _zPiezoVoltageCmd;
        private double _zPiezoPosNm;
        
        // PID 제어기 상태
        private double _prevError;
        private double _integralError;
        private readonly double _kp;
        private readonly double _ki;
        private readonly double _kd;

        public AfmQplusPiezoServo(
            double dt,
            double f0,
            double springConstantK,
            double piezoSensitivityNmPerV,
            double freqSetpointHz,
            double kp,
            double ki,
            double kd)
        {
            _dt = dt;
            _f0 = f0;
            _springConstantK = springConstantK;
            _piezoSensitivityNmPerV = piezoSensitivityNmPerV;
            _freqSetpointHz = freqSetpointHz; // 음의 값 (인력 영역, 예: -15.0 Hz)
            _kp = kp;
            _ki = ki;
            _kd = kd;
            _currentDeltaF = 0.0;
            _filteredDeltaF = 0.0;
            _zPiezoVoltageCmd = 0.0;
            _zPiezoPosNm = 0.0;
            _prevError = 0.0;
            _integralError = 0.0;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public double UpdateZ(
            double pllMeasuredFreqShiftHz,
            double oscillationAmplitudeNm,
            out double estimatedForceGradientNPerM,
            out double actualZDisplacementNm)
        {
            // 1. PLL 주파수 시프트 1차 저역통과 필터링
            _filteredDeltaF = 0.85 * _filteredDeltaF + 0.15 * pllMeasuredFreqShiftHz;
            _currentDeltaF = _filteredDeltaF;

            // 2. 상호작용력 구배(Force Gradient) 역산: k_ts = -2 * k * (Delta_f / f0)
            estimatedForceGradientNPerM = -2.0 * _springConstantK * (_filteredDeltaF / _f0);

            // 3. 주파수 오차 계산 (설정점 추종)
            double error = _filteredDeltaF - _freqSetpointHz;

            // 4. PID 제어기 연산
            _integralError += error * _dt;
            _integralError = Math.Clamp(_integralError, -20.0, 20.0); // 안티 와인드업

            double derivative = (error - _prevError) / _dt;
            _prevError = error;

            double deltaV = -(_kp * error + _ki * _integralError + _kd * derivative);
            _zPiezoVoltageCmd += deltaV;
            _zPiezoVoltageCmd = Math.Clamp(_zPiezoVoltageCmd, -100.0, 100.0); // ±100V 피에조 앰프

            // 5. 피에조 액추에이터 변위 계산
            _zPiezoPosNm = _zPiezoVoltageCmd * _piezoSensitivityNmPerV;
            actualZDisplacementNm = _zPiezoPosNm;

            return _zPiezoVoltageCmd;
        }

        public readonly double CurrentDeltaF => _currentDeltaF;
        public readonly double ZPositionNm => _zPiezoPosNm;
    }
}

4. 교육 및 국방 현장 적용 방안