1. 개요
양자 암호 키 분배(QKD) 및 장거리 우주 레이저 양자 통신에서 1550nm 통신 파장의 단일 광자를 98% 이상의 시스템 양자 효율(SDE)과 피코초 단위 지터로 검출하는 핵심 소자는 초전도 나노와이어 단일광자 검출기(Superconducting Nanowire Single-Photon Detector, SNSPD)입니다. 2K 미만의 극저온에서 나노와이어에 임계 전류(\(I_c\)) 직전의 직류 바이어스 전류(\(I_b \approx 0.95 I_c\))를 인가할 때, 단일 광자가 입사하면 국소 저항성 핫스팟(Hotspot)이 형성되며 고속 전압 펄스를 방출합니다. 그러나 높은 광자 계수율 환경이나 열 방출 불균일로 인해 줄 발열(Joule Heating)이 극저온 기판 냉각 능력을 초과하면, 나노와이어가 초전도 상태로 복귀하지 못하고 영구적인 상전도 저항 상태에 고착되는 '열적 래칭(Thermal Latching)' 현상이 발생합니다. 래칭이 발생하면 소자는 광자 검출 기능을 상실하므로 바이어스 전류를 즉각 수십 나노초 동안 차단(Quench)한 뒤 서냉 복귀시켜야 합니다. 본 고에서는 래칭 감지 및 고속 자가회복 전류 램핑을 10,000Hz 주기로 관리하는 0-Alloc C# 엔진을 구현합니다.
2. 수학적 모델링 및 지배 방정식
나노와이어 국소 핫스팟 형성 및 전자-포논 2온도 열전도 편미분방정식은 다음과 같습니다:
\[
C_e(T_e) \frac{\partial T_e}{\partial t} = \nabla \cdot (k_e \nabla T_e) - \frac{C_e}{\tau_{e-ph}} (T_e - T_p) + \rho_n J_b^2 \Theta(r_{hs} - r)
\]
핫스팟 발생 시 50\(\Omega\) 고주파 전송선로로 방출되는 전압 펄스 진폭 수식은 다음과 같습니다:
\[
V_{pulse}(t) = I_b \frac{R_{hs}(t) Z_0}{R_{hs}(t) + Z_0 + R_L}, \quad \tau_{decay} \approx \frac{L_{kinetic}}{Z_0 + R_L}
\]
줄 발열과 기판 냉각 열コンダクタンス \(G_{th}\) 간의 영구 래칭(Latching) 임계 조건식은 다음과 같습니다:
\[
P_{joule} = I_b^2 R_n \ge G_{th} (T_c - T_{bath}) \implies I_{latch} = \sqrt{\frac{G_{th} (T_c - T_{bath})}{R_n}}
\]
래치 감지 시 고속 능동 소광(Active Quench) 후 초전도 회복 지수 램프 전류 서보식은 다음과 같습니다:
\[
I_b(t) = \begin{cases} 0, & t \in [t_{latch}, t_{latch} + \Delta t_{quench}] \\ I_{opt} \left( 1 - \exp\left(-\frac{t - t_{reset}}{\tau_{ramp}}\right) \right), & t > t_{latch} + \Delta t_{quench} \end{cases}
\]
3. 실시간 C# 서보 제어기 구현
다음은 RF 검출 펄스 전압 및 DC 바이어스 모니터링 신호를 10,000Hz 주기로 샘플링하여 나노와이어 래칭을 즉시 차단하고 최적 바이어스 전류로 복귀시키는 0-Alloc C# 엔진 구현체입니다:
SnspdCurrentLatchRecoveryServoEngine.cs
using System;
using System.Runtime.CompilerServices;
namespace ZeliDesk.Stem.Snspd
{
public readonly ref struct SnspdSensorFrame
{
public readonly double MeasuredDcCurrentMicroAmps;
public readonly double TargetBiasCurrentMicroAmps;
public readonly double NanowireVoltageDropMv;
public readonly double CryoStageTemperatureKelvin;
public readonly double PhotonCountRateMcps;
public SnspdSensorFrame(
double measuredDcCurrentMicroAmps,
double targetBiasCurrentMicroAmps,
double nanowireVoltageDropMv,
double cryoStageTemperatureKelvin,
double photonCountRateMcps)
{
MeasuredDcCurrentMicroAmps = measuredDcCurrentMicroAmps;
TargetBiasCurrentMicroAmps = targetBiasCurrentMicroAmps;
NanowireVoltageDropMv = nanowireVoltageDropMv;
CryoStageTemperatureKelvin = cryoStageTemperatureKelvin;
PhotonCountRateMcps = photonCountRateMcps;
}
}
public ref struct SnspdControlCommand
{
public double CommandedBiasCurrentMicroAmps;
public bool TriggerActiveQuench;
public double HotspotResistanceOhms;
public bool IsNanowireSuperconducting;
public bool IsDetectionEfficiencyOptimal;
}
public ref struct SnspdCurrentLatchRecoveryServoEngine
{
private const double LatchVoltageThresholdMv = 15.0; // Steady DC voltage drop indicates normal state latch
private const double QuenchDurationSec = 0.0005; // 500 us quench
private const double TauRampSec = 0.001; // 1 ms recovery ramp
private readonly double _dt;
private double _timeInQuench;
private double _timeInRecovery;
private bool _isLatched;
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public SnspdCurrentLatchRecoveryServoEngine(double dt)
{
_dt = dt;
_timeInQuench = 0.0;
_timeInRecovery = 0.0;
_isLatched = false;
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public void ComputeBiasRegulation(in SnspdSensorFrame frame, out SnspdControlCommand cmd)
{
cmd = default;
double currentUa = Math.Max(0.1, frame.MeasuredDcCurrentMicroAmps);
double rHotspot = (frame.NanowireVoltageDropMv * 1000.0) / currentUa; // mV / uA = kOhm * 1000 -> Ohm
// Detect DC latching
if (frame.NanowireVoltageDropMv > LatchVoltageThresholdMv && !_isLatched)
{
_isLatched = true;
_timeInQuench = 0.0;
_timeInRecovery = 0.0;
}
double cmdBiasUa = 0.0;
bool triggerQuench = false;
if (_isLatched)
{
if (_timeInQuench < QuenchDurationSec)
{
_timeInQuench += _dt;
cmdBiasUa = 0.0; // Cut bias to zero to quench normal state
triggerQuench = true;
}
else
{
_timeInRecovery += _dt;
// Exponential recovery ramp
double rampRatio = 1.0 - Math.Exp(-_timeInRecovery / TauRampSec);
cmdBiasUa = frame.TargetBiasCurrentMicroAmps * rampRatio;
if (rampRatio >= 0.99 && frame.NanowireVoltageDropMv <= 1.0)
{
_isLatched = false; // Successfully recovered to superconducting state
cmdBiasUa = frame.TargetBiasCurrentMicroAmps;
}
}
}
else
{
// Normal superconducting operation
cmdBiasUa = frame.TargetBiasCurrentMicroAmps;
}
cmd.CommandedBiasCurrentMicroAmps = cmdBiasUa;
cmd.TriggerActiveQuench = triggerQuench;
cmd.HotspotResistanceOhms = rHotspot;
cmd.IsNanowireSuperconducting = !_isLatched && frame.NanowireVoltageDropMv <= 2.0;
cmd.IsDetectionEfficiencyOptimal = cmd.IsNanowireSuperconducting && Math.Abs(cmdBiasUa - frame.TargetBiasCurrentMicroAmps) <= 0.5;
}
}
}4. 교육 및 실무 적용 방안
- 한국방송통신대학교 양자통신/초전도물리: 단일 광자 검출 원리, 핫스팟 비선형 열전도 및 래칭 자가회복 회로 해석 교구.
- 한국전자통신연구원(ETRI) 및 국방과학연구소: 양자암호통신(QKD) 위성-지상 링크용 고속 SNSPD 바이어스 안정화 제어기 탑재.