2022 개정 교육과정 연계 온디바이스 생체물리화학·극저온나노구조 STEM 시리즈 #69

초저온 전자현미경(Cryo-EM) 급랭 유리화 액체 에탄 서보 엔진

단백질 손상 없는 비정질 얼음막 형성을 위한 초당 $10^6\text{ K/s}$ 냉각 역학 및 0-Alloc C# 펠티어-블로팅 제어 루프

ZeliDesk 구조생물학·저온물리화학팀 초저온전자현미경(Cryo-EM)·단분자동결 분과

2026년 9월 12일

1. 얼음 결정 핵형성 속도론 및 에탄 플런지 급랭 열전도 모델

물 분자가 과냉각(Supercooling)될 때 균일 핵형성 속도 $J_{nuc}(T)$는 고전 핵형성 이론(CNT)에 따라 열역학적 활성화 배리어 $\Delta G^*$와 점성 확산 활성화 에너지 $\Delta G_D$의 함수로 주어집니다:

결정화 임계 온도 구간($-40^\circ\text{C} \sim -137^\circ\text{C}$)을 통과하는 냉각 속도 $\frac{dT}{dt}$가 임계 급랭 속도 $\dot{T}_{crit} \approx 10^6\text{ K/s}$를 초과해야만 비정질 상태($T_g \approx -137^\circ\text{C}$)로 동결됩니다. 두께 $2L$인 액막의 1차원 푸리에 비정상 열전도는 다음과 같습니다:

액체 에탄 온도가 $-183^\circ\text{C}$ 미만으로 내려가면 고체로 응고되어 그리드를 찌그러뜨리고, $-160^\circ\text{C}$를 초과하면 기화로 인해 결정 얼음이 생기므로, 서보 루프는 에탄 포트 온도를 $-175^\circ\text{C} \sim -180^\circ\text{C}$의 좁은 윈도우로 고정합니다.

2. 0-Alloc C# 스택 기반 Cryo-EM 급랭 항온 서보

백금 RTD 센서와 습도 챔버 광학 센서로부터 에탄 액화 상태와 수용액 막 두께를 감지하여 펠티어 전류 및 플런지 타이밍을 지휘하는 ref struct 루프입니다.

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.CryoBiology
{
    public readonly ref struct CryoTelemetry
    {
        public readonly double EthaneCupTempC;
        public readonly double LiquidNitrogenTempC;
        public readonly double ChamberRelativeHumidityPercent;
        public readonly double OpticalGridReflectance;
        public readonly double PlungeVelocityMPerSec;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public CryoTelemetry(double ethaneTempC, double ln2TempC, double rhPercent, double reflectance, double plungeVel)
        {
            EthaneCupTempC = ethaneTempC;
            LiquidNitrogenTempC = ln2TempC;
            ChamberRelativeHumidityPercent = rhPercent;
            OpticalGridReflectance = reflectance;
            PlungeVelocityMPerSec = plungeVel;
        }
    }

    public readonly ref struct CryoControlResult
    {
        public readonly double PeltierCurrentAmps;
        public readonly double EstimatedIceThicknessNm;
        public readonly bool ReadyForPlunge;
        public readonly bool IsVitrificationGuaranteed;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public CryoControlResult(double peltierAmps, double thicknessNm, bool ready, bool guaranteed)
        {
            PeltierCurrentAmps = peltierAmps;
            EstimatedIceThicknessNm = thicknessNm;
            ReadyForPlunge = ready;
            IsVitrificationGuaranteed = guaranteed;
        }
    }

    public ref struct CryoVitrificationServo
    {
        private const double TargetEthaneTempC = -178.0; // 에탄 최적 액체 윈도우 (-178도)
        private readonly double _kpPeltier;
        private readonly double _kiPeltier;
        private double _integratedTempError;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public CryoVitrificationServo(double kp, double ki)
        {
            _kpPeltier = kp;
            _kiPeltier = ki;
            _integratedTempError = 0.0;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public CryoControlResult Step(in CryoTelemetry telemetry, double dt)
        {
            // 1. 에탄 포트 온도 오차 및 펠티어 양방향 전류 계산 (냉각/가열)
            double tempError = telemetry.EthaneCupTempC - TargetEthaneTempC;
            _integratedTempError += tempError * dt;
            double peltierCurrent = -(_kpPeltier * tempError + _kiPeltier * _integratedTempError);
            peltierCurrent = Math.Clamp(peltierCurrent, -3.0, 3.0); // +-3A 제한

            // 2. 광학 반사도 기반 시료 수막 두께 추정 (간섭 박막 두께 모델)
            // 반사도 0.05 ~ 0.20 구간이 최적 20~50nm 무정형 얼음층
            double thicknessNm = Math.Max(10.0, 250.0 * telemetry.OpticalGridReflectance);

            // 3. 플런지 준비 완료 조건 (온도 편차 1도 이내, 챔버 습도 95% 이상으로 증발 억제)
            bool tempStable = Math.Abs(tempError) < 1.0;
            bool humiditySafe = telemetry.ChamberRelativeHumidityPercent >= 95.0;
            bool readyForPlunge = tempStable && humiditySafe && (thicknessNm <= 60.0);

            // 4. 유리화(Vitrification) 보장 여부 (두께 < 50nm & 플런지 속도 > 1.5 m/s)
            bool vitrificationGuaranteed = readyForPlunge && (thicknessNm <= 50.0) && (telemetry.PlungeVelocityMPerSec >= 1.5);

            return new CryoControlResult(peltierCurrent, thicknessNm, readyForPlunge, vitrificationGuaranteed);
        }
    }
}

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