Fast design and scaling of multi-qubit gates in large-scale trapped-ion quantum computers

기술 요약: 대규모 트랩 이온 양자 컴퓨터를 위한 다중 큐비트 게이트의 신속한 설계 및 스케일링

문제 정의
트랩 이온 양자 컴퓨터는 장거리 연결성과 긴 결맞음 시간(coherence times)을 제공하여 양자 컴퓨팅의 선도적인 플랫폼으로 자리 잡고 있다. 그러나 이러한 시스템을 수백 개의 이온으로 확장하는 것은 중대한 개념적, 기술적 과제를 안겨준다. 이온들은 집단적인 운동 모드(phonon bus)를 통해 전방향(all-to-all) 연결성을 가지지만, 대규모 결정(crystal)을 위한 고충실도(high-fidelity)의 프로그래밍 가능한 다중 큐비트 얽힘 게이트를 설계하는 것은 계산적으로 매우 어렵다.

핵심적인 어려움은 게이트 설계 문제에 있다. 즉, 모든 큐비트 쌍 사이의 목표 얽힘 위상($N(N-1)/2$개의 이차 제약 조건)을 동시에 만족하면서, 구동 전력을 최소화하고 운동 모드가 바닥 상태로 돌아오도록 하는 일련의 제어 드라이브를 찾는 것이다. 이는 이온의 수 $N$에 따라 $N^2$으로 스케일링되는 이차 제약 최적화 문제(quadratically constrained optimization problem)로 정식화되며, 이는 NP-난해(NP-hard) 문제이다. 기존의 접근 방식은 문제를 작은 이온 부분 집합으로 제한하거나, 대규모 시스템의 복잡성과 제약 조건에 부딪히기 쉬운 경사 하강법 기반의 최소화 알고리즘에 의존한다.

방법론: 대규모 신속(Large-Scale Fast, LSF) 방법
저자들은 다중 큐비트 게이트 설계의 NP-난해 특성을 완화하여, 수백 개의 이온으로 구성된 결정에 대해서도 다항 시간(polynomial-time) 내에 견고하고 프로그래밍 가능한 게이트를 설계할 수 있게 하는 "Large-Scale Fast"(LSF)라는 방법론적 프레임워크를 도입한다.

1. 스펙트럼 영역 정식화: 게이트 설계를 큐비트 전이 주파수 주변에서 디튜닝된 푸리에 성분의 합으로 나타내는 스펙트럼 영역에서 정식화한다. 문제는 드라이브 톤의 진폭을 나타내는 벡터 $R$을 찾는 문제로 환원된다. 이때 $R$은 선형 제약 조건(운동 변위가 0이 됨)을 만족하면서 $\|R\|$을 최소화해야 하며, 동시에 이차 제약 조건(목표 얽힘 위상 $\phi_{nn'}$)을 충족해야 한다.
2. 제로 위상 솔루션(Zero-Phase Solution, ZPS) 시드: LSF 방법은 "제로 위상 솔루션(ZPS)"이라 불리는 시드 솔루션에 의존한다. ZPS는 모든 목표 얽힘 위상이 0인($\phi^{(0)}_s = 0$) 동차 문제(homogeneous problem)의 비자명한 해이다. 이러한 솔루션들은 특정 목표 게이트를 구하기 전에 미리 축적된다.
3. 변환 및 최적화: 알고리즘은 두 가지 주요 단계로 진행된다.
   • 변환(Conversion): ZPS를 큰 인자 $\lambda$로 스케일링하고 작은 벡터 $D$를 섭동(perturbation)시켜 목표 이차 제약 조건을 만족시킨다. 스케일링된 ZPS 주변에서 제약 조건을 선형화함으로써, 저자들은 목표 위상을 만족하는 근사 해 $R_{conv}$로 변환하기 위한 선형 시스템을 도출한다.
   • 노름 감소(Norm Reduction): 솔루션 벡터의 노름(norm)을 최소화하면서 제약 조건을 엄격하게 준수하는 반복 프로세스를 수행한다. 이는 솔루션 매니폴드를 따라 이동하며 전력을 줄이는 "노름 감소" 단계와, 목표 위상으로부터의 편차를 수정하는 "오차 감소" 단계를 교대로 수행하는 과정을 포함한다. 오차가 지정된 허용 오차(통상 $10^{-4}$) 내에 머물도록 단계 크기는 카르다노(Cardano) 및 페라리(Ferrari)의 근 찾기 방법을 사용하여 분석적으로 계산된다.
4. 집계(Aggregation): 다양한 얽힘 맵(entanglement maps)에 대해 효율적인 솔루션을 찾을 수 있도록 여러 ZPS를 집계하여 검색 공간의 다양성을 확보한다.

주요 기여 및 결과

• 다항 시간 스케일링: LSF 방법은 NP-난해인 설계 문제를 ZPS 생성(특정 대칭성에 대해 효율적으로 수행 가능)과 효율적인 다항 시간 최적화 작업으로 변환한다. 벤치마크 테스트 결과, LSF는 특히 $N > 20$인 시스템에서 표준 신뢰 영역 제한(trust-region constrained) 및 공액 경사(conjugate-gradient) 최적화 도구보다 약 10배 더 빠른 해결 시간을 달성했다.
• 구동 전력 추정: 저자들은 게이트를 구동하는 데 필요한 총 라비 주파수($\Omega_{nuc}$)에 대한 휴리스틱 추정치를 도출한다. 이 추정치는 목표 얽힘 행렬의 핵 노름(nuclear norm) $\|\phi\|_{nuc}$에 의존하며 이온 수 $N$에 따라 선형적으로 스케일링된다. 놀랍게도 시뮬레이션 결과, 요구되는 구동 전력은 특정 모드 구조나 주파수 간격에 크게 무관하며, 평균 모드 밀도 시간 척도 $T_\nu$로 정규화했을 때 하나의 보편적인 곡선으로 수렴한다.
• 최소 게이트 시간 ($T_{min}$): 본 연구는 게이트 지속 시간에 대한 근본적인 하한선을 식별한다. 임계 시간 $T_{min} \approx 2.9 \cdot T_\nu$ (여기서 $T_\nu$는 평균 모드 간격의 역수에 비례함) 미만에서는 고충실도 솔루션을 찾을 수 없다. 결정적으로 $T_{min}$은 이온 수 $N$에 따라 선형적으로 스케일링된다. 이는 동시 얽힘 연산의 수가 $N^2$으로 증가하더라도, 게이트 지속 시간은 $N$에 비례하여 증가함을 의미하며, 이는 순차적 게이트 구현 방식에 비해 스케일링 측면에서 이점을 제공한다.
• 오차 스케일링 분석: 논문은 세 가지 주요 오차 원인이 게이트 충실도에 미치는 영향을 정량화한다.
  • 운동 모드 주파수 드리프트: 오차는 주파수 드리프트와 결합 수($N^2$)에 따라 이차적으로 스케일링된다.
  • 구동 진폭 변동: 오차는 노이즈의 분산 및 결합 수와 함께 스케일링된다.
  • 운동 가열(Motional heating): 가열 유발 오차의 스케일링은 노이즈의 공간적 상관관계에 따라 달라진다. 상관관계가 없는 노이즈의 경우 충실도는 $N$의 제곱에 따라 감소하며, 강한 상관관계가 있는 노이의 경우 선형적으로 감소한다.
    LSF 프레임워크는 이러한 오차들에 대해 게이트를 견고하게 만들기 위해 선형 제약 조건을 포함할 수 있게 하여, 드리프트와 가열에 대한 민감도를 크게 낮춘다.
• 표면 코드 예시: 이 방법은 49개 이온 결정에서 단일 펄스로 33개의 이온을 얽히게 하는 표면 코드 스테빌라이저 게이트를 구현하는 데 사용되었다. 솔루션은 효율적인 전력 배분을 보여주며, "노드(node)" 이온에는 높은 구동 진폭을, "에지(edge)" 이온에는 낮은 진폭을 할당하고, 연결되지 않은 쌍 사이의 원치 않는 교차 얽힘을 효과적으로 억제한다.

의의 및 주장
본 논문은 다중 큐비트 게이트 설계라는 개념적 장벽을 극선함으로써, 수백 개의 큐비트로 트랩 이온 양자 컴퓨터를 확장하기 위한 경로를 제시한다고 주장한다. 빠르고 견고한 알고리즘(LSF)을 제공함으로써, 저자들은 이전에는 NP-난해성 때문에 까다로운 작업으로 여겨졌던 수십 개의 이온과 복잡한 다중 큐비트 게이트를 포함하는 회로의 오프라인 컴파일을 가능하게 한다.

이 연구는 게이트 시간 $T_{min}$이 선형적으로 스케일링되는 임계값보다 길다면, 대규모 결정에서도 고충실도 다중 큐비트 게이트가 물리적으로 가능하다는 것을 입증한다. 또한, 도출된 구동 전력, 게이트 시간 및 오차 원인 간의 관계는 실험 설계자들에게 실질적인 설계 지침과 허용 오차를 제공한다. 저자들은 자신의 방법이 램-디케(Lamb-Dicke) 근사 내에서 지배적인 오차 메커니즘을 다루지만, 고진폭 솔루션에서의 큐비트 캐리어와의 상호작용 및 매그너스 전개(Magnus expansion)의 고차 항에 대해서는 추가적인 검증이 필요하다고 언급한다. 이 프레임워크는 양자 오류 수정 및 최적화 응용 분야를 위한 트랩 이온 아키텍처의 장거리 연결성 잠재력을 실현하기 위한 기초적인 도구로 제시된다.