AIDT 온디바이스 STEM 음향유체학(Acoustofluidics) 표면탄성파(SAW) 세포 분리 서보 엔진

1. 개요: 초음파 정재파의 음압 노드(Node)를 이용해 세포를 손상 없이 크기별로 분리하다

형광 또는 자성 입자를 부착하는 전통적인 세포 분리 기술(FACS, MACS)은 형광체 결합 과정에서 세포 변성을 유발하고 비용과 시간이 많이 소요됩니다. 음향유체학(Acoustofluidics)은 미세유체 채널에 초음파 음향장(Acoustic field)을 인가하여 세포의 크기, 밀도, 압축률(Compressibility) 차이만을 기반으로 비침습적, 무표지(Label-free) 방식으로 고순도 연속 분리를 수행합니다.

압전 기판(예: \(128^\circ\) 회전 \(Y\)-절단 \(\text{LiNbO}_3\)) 표면에 인터디지털 트랜스듀서(Interdigital Transducer, IDT)를 마주보게 배치하고 RF 고주파를 인가하면 정재 표면탄성파(Standing Surface Acoustic Wave, SSAW)가 발생하여 유체 채널 내에 주기적인 음압 절(Pressure Node)과 배(Antinode)를 형성합니다. 본 백서에서는 고로프(Gor'kov) 음향 방사력 포텐셜을 정밀 모델링하고 RF 주파수와 위상을 실시간 제어하는 10,000Hz 0-Alloc C# 서보를 구현합니다.

2. 고로프 포텐셜 및 음향 콘트라스트 계수 수식 모델

파장보다 작은 반경 \(a\)의 단일 입자가 1차 음향장(음압 \(p_1\), 속도장 \(\mathbf{v}_1\)) 내에서 받는 음향 방사력(Acoustic Radiation Force, \(\mathbf{F}_{\text{rad}}\))은 고로프 포텐셜(Gor'kov Potential, \(U\))의 음의 그래디언트로 주어집니다:

\\[ \mathbf{F}_{\text{rad}} = -\nabla U \\]

\\[ U = V_0 \left[ f_1 \frac{\langle p_1^2 \rangle}{2 \rho_m c_m^2} - f_2 \frac{3 \rho_m \langle v_1^2 \rangle}{4} \right] \\]

여기서 \(V_0 = \frac{4}{3}\pi a^3\)은 입자 부피, \(\rho_m, c_m\)은 유체 매질의 밀도와 음속입니다. 모노폴 및 다이폴 산란 인자 \(f_1, f_2\)는 다음과 같습니다:

\\[ f_1 = 1 - \frac{\kappa_p}{\kappa_m}, \quad f_2 = \frac{2(\rho_p - \rho_m)}{2\rho_p + \rho_m} \\]

여기서 \(\kappa = \frac{1}{\rho c^2}\)는 등온 압축률입니다. 1차원 평면 정재파 \(p_1(y) = p_0 \cos(k y)\) 하에서 횡방향 음향 방사력은 다음과 같이 전개됩니다:

\\[ F_{\text{rad}}(y) = 4 \pi \Phi(\beta, \rho) k a^3 E_{\text{ac}} \sin(2 k y) \\]

여기서 \(E_{\text{ac}} = \frac{p_0^2}{4 \rho_m c_m^2}\)는 음향 에너지 밀도이며, \(\Phi\)는 음향 콘트라스트 계수(Acoustic Contrast Factor)입니다:

\\[ \Phi = \frac{5 \rho_p - 2 \rho_m}{2 \rho_p + \rho_m} - \frac{\kappa_p}{\kappa_m} \\]

대부분의 생체 세포(백혈구, 암세포)는 물 속에서 \(\Phi > 0\)이므로 음압 절(Node, \(p_1 = 0\))로 이동하며, 이동 속도는 입자 반경의 제곱(\(a^2\))에 비례하므로 크기가 큰 순환종양세포(CTC, \(\sim 20\,\mu\text{m}\))가 작은 적혈구(\(\sim 6\,\mu\text{m}\))보다 훨씬 빠르게 중심 유로로 포집 분리됩니다.

3. 실시간 음향유체 세포 분리 서보 아키텍처 다이어그램

graph TD A[혈액 및 완충 희석액 주입 펌프 Q_in] --> B[PDMS 마이크로 채널 입구 3-inlet 시스 유동] C[양방향 IDT RF 신호 발생기 f_RF ~ 20 MHz] --> D[LiNbO3 기판 정재 표면탄성파 SSAW 생성] B & D --> E[채널 내 음향 방사력 F_rad ~ a^3 * Phi * sin 2ky 작용] E --> F[CTC 대형 암세포 음압 절 중심 정렬 / 적혈구 통과] F --> G[출구 3-outlet 분기 및 고속 라인스캔 카메라 관측] G --> H[CTC 순도 및 회수율 계측 신호] H --> I[10,000Hz 0-Alloc C# IDT RF 파워 & 위상 가변 서보] I --> C

4. 고성능 0-Alloc C# 음향 영동 분리 서보 구현

유량과 세포 크기 분포에 맞추어 RF 인가 파워 및 분리 효율을 실시간 제어하는 0-Alloc C# 클래스입니다:

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.Acoustofluidics
{
    public readonly record struct AcoustofluidicParameters(
        double ChannelWidthM,      // m (예: 200 um = 200e-6)
        double RfFrequencyHz,      // Hz (예: 20 MHz = 20e6)
        double MediumDensityRhoM,  // kg/m^3 (물 1000)
        double MediumSpeedOfSound, // m/s (물 1490)
        double MediumViscosityEta  // Pa*s (0.001)
    );

    public ref struct AcousticSeparationState
    {
        public double AcousticEnergyDensityJm3;
        public double ContrastFactorPhi;
        public double RadiationForceNewton;
        public double TerminalVelocityMps;
        public double LateralDeflectionUm;
        public bool IsSeparated;
    }

    public static class AcousticTweezersServo
    {
        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public static void CalculateCellTrajectory(
            in AcoustofluidicParameters p,
            double particleRadiusM,    // m (예: 10 um 암세포)
            double particleDensityRho, // kg/m^3 (1080)
            double particleCompressKp, // Pa^-1 (3.8e-10)
            double acousticPressurePa, // Pa (예: 1.0 MPa = 1e6)
            double transitTimeSeconds, // 체류 시간
            ref AcousticSeparationState state)
        {
            // 1. 매질 압축률 kappa_m = 1 / (rho_m * c_m^2)
            double kappaM = 1.0 / (p.MediumDensityRhoM * p.MediumSpeedOfSound * p.MediumSpeedOfSound);

            // 2. 음향 콘트라스트 계수 Phi = (5*rho_p - 2*rho_m) / (2*rho_p + rho_m) - (kappa_p / kappa_m)
            double densTerm = (5.0 * particleDensityRho - 2.0 * p.MediumDensityRhoM) / (2.0 * particleDensityRho + p.MediumDensityRhoM);
            double compTerm = particleCompressKp / kappaM;
            state.ContrastFactorPhi = densTerm - compTerm;

            // 3. 음향 에너지 밀도: E_ac = p0^2 / (4 * rho_m * c_m^2)
            state.AcousticEnergyDensityJm3 = (acousticPressurePa * acousticPressurePa) / (4.0 * p.MediumDensityRhoM * Math.Pow(p.MediumSpeedOfSound, 2));

            // 4. 최대 음향 방사력: F_rad_max = 4 * pi * Phi * k * a^3 * E_ac
            double k = (2.0 * Math.PI * p.RfFrequencyHz) / p.MediumSpeedOfSound;
            double a3 = Math.Pow(particleRadiusM, 3);
            state.RadiationForceNewton = 4.0 * Math.PI * state.ContrastFactorPhi * k * a3 * state.AcousticEnergyDensityJm3;

            // 5. 종단 속도 (스토크스 항력 평형: 6 * pi * eta * a * v = F_rad)
            double dragFactor = 6.0 * Math.PI * p.MediumViscosityEta * particleRadiusM;
            state.TerminalVelocityMps = Math.Abs(state.RadiationForceNewton) / Math.Max(1e-12, dragFactor);

            // 6. 횡방향 변위
            state.LateralDeflectionUm = state.TerminalVelocityMps * transitTimeSeconds * 1e6;

            // 분리 판정: 변위가 채널 폭의 20% 이상 이동
            state.IsSeparated = state.LateralDeflectionUm > (p.ChannelWidthM * 0.2 * 1e6);
        }
    }
}

5. 교육 연수 및 조기 암 진단/화생방 감시 국방 적용 지침

방통대 직무연수 및 국방 야전 혈액 생화학 분석

물리/생물 교과 융합: 압전 효과와 표면탄성파 전파, 음압의 절과 배, 고로프 음향 방사력과 스토크스 점성 유체 항력의 평형을 컴퓨터 그래픽 교구로 시연.
국방 군진의학 및 야전 독소 감별: 전투 현장에서 채취한 전혈로부터 원심분리기 없이 미세유체 칩 상에서 백혈구 및 병원균(탄저, 페스트)을 무표지 초고속 정제 분리.