Disentangling the Impact of Quasiparticles and Two-Level Systems on the Statistics of Superconducting Qubit Lifetime

1. 연구 배경 및 문제 제기 (Problem)

• 배경: 초전도 큐비트는 오류 정정이 가능한 양자 컴퓨터 구축을 위한 유망한 후보입니다. 지난 수십 년간 재료 과학 및 나노 가공 기술의 발전으로 큐비트 수명 ($T_1$) 이 나노초에서 밀리초 단위로 크게 향상되었습니다.
• 문제: 그러나 큐비트의 수명 ($T_1$) 및 결맞음 시간 ($T_2$) 이 시간과 주파수 영역에서 변동 (fluctuations) 을 보인다는 것이 관찰되었습니다. 특히 $T_1$이 긴 큐비트일수록 변동 폭이 큰 경향이 있습니다. 이러한 변동은 양자 오류 정정과 확장성에 심각한 장애물이 됩니다.
• 기존 이론의 한계: $T_1$ 변동의 주요 원인으로 표면 유전체 내의 **2-레벨 시스템 (TLS)**과 **준입자 (Quasiparticles, QP)**가 지목되어 왔으나, 두 요인이 동일한 장치에서 어떻게 상호작용하며 변동에 기여하는지를 정량적으로 분리 (disentangle) 하고 이해하는 연구는 부족했습니다.

2. 연구 방법론 (Methodology)

• 실험 장치: 페르미 국립 가속기 연구소 (Fermilab) 의 SQMS 센터에서 제작된 3 개의 고정 주파수 트랜스몬 (transmon) 큐비트 (A, B, C) 를 사용했습니다.
  • 기하학적 차이: 큐비트 A 는 작은 패드 ($120 \times 510 \mu m^2$, 간격 $20 \mu m$), 큐비트 B 와 C 는 큰 패드 ($150 \times 720 \mu m^2$, 간격 $150 \mu m$) 를 가집니다.
  • 재료 차이: A 와 B 는 Nb 박막을 사용했고, C 는 Nb 위에 Ta(탄탈륨) 박막을 증착하여 표면 유전체를 차폐했습니다.
• 측정 조건: 7 mK 에서 153 mK 까지의 온도 범위에서 각 온도당 약 72 시간 동안 연속적인 $T_1$ 측정을 수행하여 데이터의 통계적 특성을 분석했습니다.
• 분석 기법:
  • 큐비트 탈분극율 ($\Gamma_1 = 1/T_1$) 의 시간적 변동 데이터를 푸리에 변환하여 **스펙트럼 밀도 (Spectral Density)**를 추출했습니다.
  • $\Gamma_1$의 변동을 TLS 와 QP 의 기여도로 분리하기 위해, 저주파수에서의 $1/f$ 노이즈 (TLS-TLF 상호작용) 와 고주파수/고온에서의 화이트 노이즈 (QP 수 변동) 를 모델링하여 피팅했습니다.
  • QP 확산 모델 (Diffusion Model) 을 사용하여 비평형 및 평형 QP 가 접합부에 도달하는 밀도와 변동성을 정량화했습니다.

3. 주요 기여 및 결과 (Key Contributions & Results)

가. 변동 원인의 분리 및 주파수 의존성

• 저온 (7 mK): $\Gamma_1$ 변동은 $1/f$ 노이즈 특성을 보이며, 이는 소수의 TLS 가 열적으로 요동치는 TLF(2-레벨 플럭추에이터) 와 상호작용하여 발생함을 확인했습니다.
• 고온 (153 mK): 변동은 화이트 노이즈 (White Noise) 특성을 보이며, 이는 QP 수의 요동 (생성 및 재결합) 이 지배적인 원인이 됨을 확인했습니다.
• 전환: 온도가 상승함에 따라 지배적인 노이즈 원인이 TLS 에서 QP 로 전환되는 과정을 명확히 관측했습니다.

나. 기하학적 크기와 QP 변동의 상관관계

• 작은 패드 (큐비트 A) 의 취약성: 작은 기하학적 크기를 가진 큐비트 A 는 큰 패드 (B, C) 에 비해 QP 및 TLS 변동에 훨씬 더 민감하게 반응했습니다. 특히 QP 로 인한 $\Gamma_1$ 분산 ($\sigma^2_{QP}$) 이 작은 패드에서 최소 10 배 이상 크게 관측되었습니다.
• 효율적 부피 (Effective Volume) 개념 도입:
  • 비평형 QP 밀도 ($x^0_{QP}$) 와 변동 분산 사이의 관계를 설명하기 위해 유효 부피 ($V_{eff}$) 개념을 도입했습니다.
  • 확산 모델에 따르면, 작은 패드에서는 QP 가 접합부까지 이동하는 거리가 짧아 비평형 QP 밀도가 높게 유지되며, 이는 $V_{eff} \propto (1/x^0_{QP})^3$ 관계로 설명됩니다.
  • 실험 데이터와 확산 모델 시뮬레이션은 작은 패드 (A) 와 큰 패드 (B, C) 간의 QP 밀도 비율 (약 2.5 배) 과 변동성 차이를 정량적으로 일치시켰습니다.

다. 재료 공학적 개선 (Ta 캡슐화)

• 큐비트 C(Nb 위에 Ta 증착) 는 표면 유전체 손실 ($\delta_{TaOx}$) 이 Nb 산화물 ($\delta_{NbOx}$) 보다 약 10 배 낮아, TLS 로 인한 평균 탈분극율 ($\Gamma_{TLS}$) 이 크게 개선되었습니다.
• 그러나 QP 변동에 대해서는 기하학적 크기 (패드 면적) 가 재료 차폐보다 더 중요한 영향을 미치는 것으로 나타났습니다.

라. 이론적 모델의 검증

• QP 수의 포아송 분포 (Poisson distribution) 를 가정하여 유도한 분산 공식이 실험적으로 추출된 $\Gamma_1$ 분산 데이터와 잘 일치함을 확인했습니다.
• TLS 에 의한 분산은 온도에 따라 약한 의존성을 보였으며, 이는 TLS 폭 (linewidth) 증가와 TLF 스위칭 빈도 증가라는 두 가지 상반된 효과가 경쟁하기 때문임을 시뮬레이션을 통해 설명했습니다.

4. 의의 및 결론 (Significance)

• 설계 최적화 가이드: 초전도 큐비트의 수명 변동 (fluctuations) 을 줄이기 위해서는 단순히 평균 수명 ($T_1$) 을 높이는 것뿐만 아니라, **기하학적 구조 (패드 크기)**와 **재료 (표면 유전체)**를 동시에 최적화해야 함을 시사합니다.
• 변동성 제어 전략:
  • TLS 변동은 큐비트나 TLS 를 구동 (driving) 하는 방식으로 안정화할 수 있는 가능성이 있으나, QP 변동은 아직 명확한 저감 프로토콜이 부족합니다.
  • 작은 기하학적 구조가 QP 변동에 취약하다는 발견은 향후 대규모 양자 프로세서 설계 시 패드 크기와 QP 트랩 (trap) 배치 전략에 중요한 통찰을 제공합니다.
• 정량적 모델링: "유효 QP 부피" 개념을 통해 서로 다른 기하학적 구조를 가진 큐비트 간의 수명 변동 차이를 정량적으로 설명할 수 있는 이론적 틀을 마련했습니다.

이 연구는 초전도 큐비트의 결맞음 한계를 이해하고, 오류 정정이 가능한 양자 컴퓨팅을 위한 안정적인 큐비트 설계에 필수적인 기초 데이터를 제공한다는 점에서 중요한 의미를 가집니다.