AIDT 온디바이스 STEM 그래핀 나노리본(GNR) 양자 점 접촉(QPC) 전도도 양자화 서보 제어기

1. 개요

2차원 탄소 동소체인 그래핀(Graphene)을 원자 몇 개 너비로 정밀 재단한 그래핀 나노리본(GNR)은 양자 가둠(Quantum Confinement) 효과로 인해 밴드갭이 열리며 뛰어난 전계효과 트랜지스터 및 양자 수송 소자로 주목받고 있습니다. 특히 나노스케일의 협착부(Constriction)를 형성하는 양자 점 접촉(Quantum Point Contact, QPC) 구조에서는 페르미 파장 수준의 채널 폭에 의해 통과 가능한 1차원 서브밴드(Subband) 모드의 개수가 정수배로 제한되며, 미분 전도도가 양자 콘덕턴스 단위 \(G_0 = 2e^2/h \approx 77.48 \mu\text{S}\)의 정수배로 양자화되는 계단형 플래토(Plateau)를 나타냅니다. 그러나 전하 불순물에 의한 쿨롱 산란, 열적 요동 및 접촉 저항 변동은 전도도 플래토를 왜곡시킵니다. 본 고에서는 Landauer 수송 모델을 기반으로 분할 스플릿 게이트 전압을 정밀 제어하여 양자화 전도도 플래토를 10,000Hz 주기로 능동 고정(Lock-in)하는 0-Alloc C# 서보 엔진을 구현합니다.

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

탄도 전자 수송(Ballistic Transport) 영역에서 1차원 양자 채널을 통과하는 미분 전도도 \(G\)는 Landauer-Büttiker 공식으로 기술됩니다:

여기서 \(\mathcal{T}_n(E_F)\)는 제 \(n\)번째 1차원 서브밴드의 투과 확률입니다. 새들 포인트(Saddle-Point) 정전기 포텐셜 장벽 모델에서 투과 계수는 Büttiker 공식에 의해 다음과 같이 주어집니다:

게이트 전압 \(V_g\)에 따른 미분 전도도의 변화율 \(dG/dV_g\)는 각 플래토(Plateau)의 중앙에서 0에 도달합니다:

양자 점 접촉 서보의 제어 목적은 측정된 전류 \(I_{sd}\) 및 전압 \(V_{sd}\)로부터 계산된 \(G\)를 목표 정수 플래토 \(N^* \cdot (2e^2/h)\)에 고정하고 쿨롱 전하 드리프트를 보상하는 것입니다:

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

다음은 전류 계측 신호로부터 전도도 양자화 계단을 추적하고 분할 게이트 전압을 10,000Hz 주기로 보정하는 무할당 C# 엔진 구현체입니다:

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Core.Quantum
{
    public ref struct GrapheneQpcQuantizationServo
    {
        private readonly double _dt;
        private readonly double _g0; // 2e^2 / h = 77.480917 umho
        private double _gateVoltageV;
        private int _targetModeN;
        
        // 저역통과 필터 및 미분기
        private double _filteredConductance;
        private double _prevConductance;
        private double _conductanceSlope;
        
        // 제어기 이득
        private readonly double _kp;
        private readonly double _ki;
        private double _integralError;

        public GrapheneQpcQuantizationServo(
            double dt,
            double initialGateV,
            int targetModeN,
            double kp,
            double ki)
        {
            _dt = dt;
            _g0 = 7.748091729e-5; // Siemens (2e^2/h)
            _gateVoltageV = initialGateV;
            _targetModeN = targetModeN;
            _kp = kp;
            _ki = ki;
            _filteredConductance = targetModeN * _g0;
            _prevConductance = _filteredConductance;
            _conductanceSlope = 0.0;
            _integralError = 0.0;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public double Update(
            double drainCurrentA,
            double sourceDrainVoltageV,
            out double measuredConductanceSiemens,
            out double plateauErrorRatio)
        {
            // 1. 순시 미분 전도도 계산
            double rawG = drainCurrentA / Math.Max(1e-6, sourceDrainVoltageV);

            // 2. 1차 저역통과 필터링
            _filteredConductance = 0.90 * _filteredConductance + 0.10 * rawG;
            measuredConductanceSiemens = _filteredConductance;

            // 3. 전도도 슬로프(미분치) 연산
            _conductanceSlope = (_filteredConductance - _prevConductance) / _dt;
            _prevConductance = _filteredConductance;

            // 4. 목표 양자화 플래토 (N * G0) 오차 계산
            double targetG = _targetModeN * _g0;
            double errorG = _filteredConductance - targetG;
            plateauErrorRatio = errorG / _g0;

            // 5. 게이트 전압 보정 연산 (PI 서보)
            _integralError += errorG * _dt;
            _integralError = Math.Clamp(_integralError, -1e-5, 1e-5);

            double deltaGateV = -(_kp * errorG + _ki * _integralError);
            _gateVoltageV += deltaGateV;
            _gateVoltageV = Math.Clamp(_gateVoltageV, -5.0, 5.0); // ±5V 게이트 바이어스 범위

            return _gateVoltageV;
        }

        public readonly double CurrentGateVoltage => _gateVoltageV;
        public readonly double QuantizedPlateauIndex => _filteredConductance / _g0;
    }
}

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