Device variability of Josephson junctions induced by interface roughness

기술 요약: 계면 거칠기에 의해 유도된 조셉슨 접합의 소자 가변성

문제 정의
결함 허용 양자 컴퓨팅(예: 쇼어 알고리즘)을 위해 수천 개의 큐비트가 요구되는 규모로 초전도 양자 프로세서가 확장됨에 따라, 소자 간 가변성이 중요한 병목 현상으로 부상했습니다. 초전도 큐비트, 특히 트랜스몬(transmon)은 터널링 영역에서 작동하는 Al/AlOx/Al 조셉슨 접합에 의존합니다. 이 영역에서 임계 전류 $I_c$—그리고 결과적으로 조셉슨 에너지 $E_J$—는 미시적인 장벽 세부 사항에 지수적으로 의존합니다. 국부적인 두께 변동이나 오염과 같은 Al/AlOx 계면 기하학의 미세한 변화조차 큐비트 전이 주파수와 게이트 충실도(fidelity)에 상당한 변화를 유발할 수 있습니다. 미세 구조의 재현성을 높이기 위한 공정 최적화 및 현미경 관찰 등 광범위한 실험적 노력이 진행되어 왔으나, 미세 구조적 무질서(특히 계면 거칠기)와 접합 수준의 에너지 통계 사이를 연결하는 정량적인 이론적 프레임워크는 부족한 실정이었습니다.

방법론
저자들은 원자 규모의 계면 무질서와 거시적인 접합 수송 특성을 연결하는 정량적 모델링 프레임워크를 개발했습니다. 이 접근 방식은 세 가지 주요 구성 요소를 통합합니다:

1. 계면 거칠기 모델링: Al/AlOx 계면을 제곱 평균 제곱근(RMS) 거칠기 진폭 $\sigma$와 횡방향 상관 길이 $\xi$로 특징지어지는 가우시안 랜덤 필드 $h(x,y)$로 모델링합니다. 이 매개변수들은 기존 문헌(예: 결정립 크기 및 두께 히스토그램에 대한 TEM 연구)과 새로운 분자 역학(MD) 시뮬레이션을 통해 도출됩니다. MD 시뮬레이션은 CHGNet 머신 러닝 포텐셜을 사용하여 이상적인 조건 하에서의 Al 표면 산화를 모델링하여 $\sigma$를 추정합니다.
2. 양자 수송 근사: 실제 접합 크기($\sim 100 \text{ nm} \times 100 \text{ nm}$)에 대해 보고리우보-드 브라위너(Bogoliubov–de Gennes, BdG) 방정식을 푸는 것이 계산적으로 불가능하므로, 저자들은 두 가지 핵심 근사를 사용합니다:
   • 암베가카르-바로토프(Ambegaokar–Baratoff, AB) 관계식: 임계 전류 $I_c$를 정상 저항 $R_N$과 연결하여 계산합니다. 이는 전체 에너지 스펙트럼에 대한 적분을 피하고 페르미 에너지에서의 투과 계수에 집중함으로써 계산 복잡성을 줄여줍니다.
   • 국부 두께 근사: 페르미 파장($\lambda_F \approx 3.6$ Å)과 장벽 감쇄 길이($\lambda_D \approx 1.9$ Å)가 거칠기의 상관 길이($\xi \sim 10$ nm)보다 훨씬 작기 때문에, 접합은 독립적인 터널링 채널들의 앙상블로 취급됩니다. 총 컨덕턴스는 국부적인 장벽 두께가 거칠기 프로파일에 따라 변하는 국부 컨덕턴스 밀도의 적분으로 근사됩니다.
3. 수치 시뮬레이션: 모델은 가우시안 랜덤 필드 모델로 생성된 다양한 터널 장벽 프로파일을 가진 5,000개의 조셉슨 접합 앙상블에 적용되었습니다.

주요 결과
수치 시뮬레이션은 $E_J$의 통계적 거동에 대한 몇 가지 정량적 통찰을 제공합니다:

• 로그 정규 분포(Log-Normal Distribution): 앙상블 전체에서 조셉슨 에너지 $E_J$의 분포는 가우시안이 아니라 로그 정규 분포를 따릅니다. 이는 $E_J$가 정규 분포를 따르는 장벽 두께에 지수적으로 의존하기 때문입니다. 결과적으로 확률 밀도 함수는 더 큰 값 쪽으로 치우친 형태를 보입니다.
• RMS 거칠기($\sigma$)의 영향: RMS 거칠기 진폭 $\sigma$가 증가하면 평균값 $\langle E_J \rangle$과 분산 $\text{Var}(E_J)$이 모두 크게 증가합니다. 평균값이 증가하는 것은 터널링의 지수적 민감도 때문이며, 드물게 발생하는 비정상적으로 얇은 영역들이 전체 컨덕턴스를 지배하게 됩니다.
• 상관 길이($\xi$)의 영향: 횡방향 상관 길이 $\xi$가 증가하면 평균 $\langle E_J \rangle$은 약간 감소하지만, 분산 $\text{Var}(E_J)$은 현저히 증가합니다. 분산이 증가하는 이유는 더 큰 $\xi$가 통계적으로 독립적인 터널링 채널의 수를 줄여(자기 평균화 효과 감소) 샘플 간의 변동성을 높이기 때문입니다. 평균이 약간 감소하는 것은 유한 크기 효과 때문인데, $\xi$가 접합 치수에 가까워질수록 터널링을 지배하는 극도로 얇은 장벽 영역이 실현될 확률이 감소합니다.
• 큐비트 주파수 가변성: 도출된 $E_J$ 통계와 전형적인 충전 에너지 $E_C$를 사용하여, 모델은 테스트된 특정 매개변수 하에서 평균 6.145 GHz 주변의 큐비트 전이 주파수 $\omega_{01}/2\pi$의 표준 편차가 약 0.364 GHz가 될 것이라고 예측합니다.

의의 및 주장
본 논문은 표면 거칠기가 어떻게 조셉슨 에너지의 가변성을 유도하는지에 대해 "정량적이고 직관적인 그림"을 제공한다고 주장합니다. 미세 구조 매개로 $\sigma$와 $\xi$를 직접 연결함으로써, 이 연구는 대규모 양자 프로세서의 수율 및 균일성 문제를 이해하기 위한 프레임워크를 제공합니다. 저자들은 자신들의 결과가 "접합 설계에 직접적인 관련성"이 있음을 강조하며, 이는 계면 거칠기 매개변수를 제어하는 것이 소자 가변성을 최소화하는 데 매우 중요하다는 것을 시사합니다.

저자들은 본 연구의 범위에 대해 신중한 입장을 유지합니다. 그들은 모델이 접합 내의 두 계면에 대해 $\sigma$와 $\xi$가 독립적이고 동일하다고 가정하며, 이 가정은 다양한 성장 조건 하에서 실험적 검증이 필요하다고 언급합니다. 또한, 향후 연구 방향으로 이 모델을 제작된 접합의 실험적 특성 분석과 연결하고, 접합과 회로 소자의 공동 설계를 가능하게 하기 위해 자신들의 초전도 솔버를 마이크로웨이브 엔지니어링 소프트웨어와 통합하는 것을 제시합니다. 이 논문은 새로운 실험 프로토콜을 제안하거나 가변성 문제를 완전히 해결한다고 주장하는 것이 아니라, 계면 거칠기의 영향을 정량화하기 위한 이론적 기준을 확립하는 데 목적이 있습니다.