Shortcuts for Adiabatic and Variational Algorithms in Molecular Simulation

문제 정의

현재의 분자 시뮬레이션을 위한 양자 알고리즘은 근미래적 하드웨어의 제한된 결맞음 시간(coherence times)과 게이트 오류로 인해 상당한 제약에 직면해 있다. 양자 위상 추정(QPE)은 결함 허용(fault-tolerant) 이점을 향한 경로를 제시하지만, 깊은 회로 요구 사항으로 인해 현재의 장치에서는 여전히 비실용적이다. 결과적으로, 단열 양자 컴퓨팅(AQC)과 변분 양자 고유치 솔버(VQE)가 선호되는 대안으로 꼽힌다. 그러나 AQC 구현은 유한 시간 진화 과정 중 디아바틱 전이(diabatic transitions)를 피하기 위해 종종 깊은 회로를 필요로 하는 반면, VQE의 성능은 안사츠(ansatz)의 선택에 결정적으로 의존한다. 즉, 노이즈 축적과 최적화의 어려움을 최소화하기 위해 압축적이면서도 표현력이 풍부한 회로가 필요하다. 저자들은 이러한 노이즈 제한 환경에서 정확도와 구현 비용 사이의 균형을 맞추는 과제를 다룬다.

방법론

본 논문은 디지털화된 단열 방식과 변분 방식을 모두 향상시키기 위해 단열로의 지름길(Shortcuts to Adiabaticity, STA), 구체적으로는 반단열(Counterdiabatic, CD) 구동에 기반한 통합 전략을 제안한다.

1. 반단열 형식론 (Counterdiabatic Formalism):

   • 저자들은 근사적인 CD 해밀토니안, $\hat{H}_{CD}(t) = \dot{\lambda}(t)\hat{A}_\lambda$를 활용한다. 여기서 $\hat{A}_\lambda$는 단열 게이지 포텐셜(Adiabatic Gauge Potential, AGP)이다.
   • 정확한 AGP는 다루기 까다로우므로, 저자들은 **절단된 중첩 교환자 전개(truncated nested-commutator expansion)**를 사용하여 근사적인 AGP를 구축한다:
     $$\hat{A}^{(l)}_\lambda = i \sum_{k=1}^l \alpha_k [\hat{H}_\lambda, [\hat{H}_\lambda, \dots [\hat{H}_\lambda, \partial_\lambda \hat{H}_\lambda]]]$$
     여기서 $l$은 절단 차수이다. 저자들은 표현력과 비용 사이의 실질적인 타협점으로 $l=1$에 집중한다.
   • 계수 $\vec{\alpha}$는 작용 범함수(action functional)를 최소화함으로써 변분적으로 결정된다.

2. 디지털화된 단열 양자 컴퓨팅 (DAQC):

   • 상태 준비를 가속화하기 위해 CD 해밀토니안을 분자 해밀토니안에 추가한다.
   • 진화 과정은 1차 트로터 분해(first-order Trotter decomposition)를 사용하여 디지털화된다. 저자들은 CD 보정으로 인한 해밀토니안 항의 증가와 목표 정확도에 도달하기 위해 필요한 총 진화 시간의 감소 사이의 트레이드오프를 분석한다.

3. 변분 양자 고유치 솔버 (VQE) - AGA 및 AGAR:

   • 저자들은 CD 구조를 **단열 게이지 안사츠(Adiabatic Gauge Ansatz, AGA)**라는 변분 안사츠로 변환한다.
   • AGA에서, 근사적 AGP(중첩 교환자 전개로부터 생성됨)에 의해 생성된 파울리 문자열(Pauli-string) 성분들은 독립적인 변분 생성자로 승격된다:
     $$\hat{U}^{(l)}_{AGA}(\vec{\theta}) = \prod_{j=1}^m \exp(-i\theta_j \hat{A}_j)$$
   • 회로 깊이와 게이트 수를 줄이기 위해, 저자들은 가중치(비항등 연산자의 개수)가 2 이하인 파울리 생성자만을 유지하는 축소된 변형인 AGAR를 도입한다.

주요 기여

• CD 보조 DAQC: CD 보정이 상태 준비를 가속화하여, 구현된 연산자의 총 개수가 증가함에도 불구하고 빠른 진화 영역(fast-evolution regime)에서 목표 정확도에 도달할 수 있음을 입증하였다.
• AGA 및 AGAR 구축: AGP 연산자 구조로부터 직접 유도된 물리적으로 동기화된 변분 안사츠(AGA)와 하드웨어 효율성을 위해 설계된 축소된 버전(AGAR)을 제안하였다.
• 하드웨어 벤치마킹: 이 안사츠들을 IBM basquecountry(156-큐비트 Heron r2) 초전도 프로세서 상에서 구현하였으며, 탈분극 노이즈(depolarizing noise) 하에서의 시뮬레이션도 병행하였다.

결과

• 단열 시뮬레이션: LiH 및 BeH$_2$에 대한 벤치마크 결과, CD 보조 진화가 광범위한 트로터 스텝 크기에 걸쳐 에너지 수렴을 개선함을 보여준다. 빠른 진화 영역에서, 특정 정확도에 도달하기 위해 필요한 순 연산자 수의 실질적인 감소가 달성되었다.
• 변분 시뮬레이션 (무노이즈):
  • LiH의 경우, 1차 AGA(AGA(1))는 상당한 결합 거리 범위에서 화학적 정확도(오차 < 1 kcal/mol)를 달성하며, 일반적으로 UCCSD, 3-UpCCGSD, ADAPT-Givens보다 우수한 성능을 보인다.
  • BeH$_2$의 경우, AGA(1)와 AGAR(1) 모두 평형 구조에서 화학적 정확도에 도달한다. 그러나 늘어난(stretched) 구조에서는 모든 안사츠의 성능이 저하되며, AGAR가 AGA보다 약간 더 큰 오차를 보이는데, 이는 늘어난 구조의 난이도가 증가함에 따른 결과와 일치한다.
  • 자원 효율성: AGAR는 특히 늘어난 구조의 BeH$_2$에 대해 UCCSD 및 전체 AGA와 비교하여 변분 파라미터와 파울리 생성자의 수를 크게 줄인다.
• 하드웨어 및 노이즈 성능:
  • IBM 프로세서 실험 결과, 오류 완화 기술(resilience level 1 및 2)을 사용할 때 AGAR가 테스트된 안사츠들(UCCSD 포함) 중 가장 작은 에너지 불일치를 나타냈다.
  • 탈분극 노이즈를 적용한 시뮬레이션은 축소된 AGAR 안사츠가 다양한 노이즈 강도($p \sim 10^{-3} - 10^{-2}$) 범위에서 UCCSD보다 더 작은 에너지 오차를 유지함을 확인시켜 주었다. 이는 감소된 회로 복잡성 덕분에 게이트 오류의 축적이 적기 때문이다.

의의 및 주장

본 논문은 AGP에서 영감을 받은 연산자 구조가 어떻게 역학적(단열적) 접근법과 변분적 접근법 사이의 가교 역할을 할 수 있는지를 보여준다고 주장한다. 주요 의의는 정확도와 구현 비용 사이의 조절 가능한 트레이드오프를 제공한다는 점에 있다:

• AQC 프레임워크에서, CD 구동은 디아바틱 전이를 억제하는 메커니즘을 제공하여, 빠른 영역에서 더 적은 총 자원으로 더 빠른 진화를 가능하게 한다.
• VQE 프레임워크에서, AGA와 AGAR는 회로 설계를 가이드하기 위해 역학적 정보를 활용하는 물리적으로 동기화된 압축 안사츠를 제공한다.
• 결정적으로, 결과는 현재의 노이즈 제한 하드웨어에서 안사츠의 복잡성을 줄이는 것(AGAR를 통해)이 UCCSD와 같은 더 크고 표현력이 풍부한 회로를 사용하는 것보다 정확도 측면에서 더 유익할 수 있음을 시사한다. 이는 오류 축적을 제한하기 때문이다.

저자들은 현재의 결과가 소규모 분자(LiH, BeH$_2$)에 국한되어 다소 미미한 수준이지만, AGP 프레임워크가 미래의 양자 프로세서를 위한 효율적인 알고리즘을 설계하는 유망한 경로를 제공한다고 결론짓는다. 그들은 회로의 압축성을 더욱 최적화하기 위해 AGP 유도 공간과 적응형 선택 전략을 결합하는 후속 연구를 제안한다.