Nuclear scattering phase shifts with factorized geometric-time RODEO on a quantum processor

기술 요약: 양자 프로세서 상의 인수축적 기하학적-시간 RODEO를 이용한 핵 산란 위상 변화(Phase Shifts)

문제 정의
양자 컴퓨팅으로 계산된 스펙트럼으로부터 핵 산란 위상 변화를 재구성하는 것은 매우 까다로운 과제이다. 양자 고유값 솔버(Quantum Eigensolver)는 가두어진 계(trap) 내의 이산적인 에너지 준위를 추정할 수 있지만, 이러한 이산 상태를 연속체 관측량(자유 공간 산란)과 연결하기 위해서는 고전적인 후처리, 구체적으로 수정 유효 범위 전개(Modified Effective Range Expansion, MERE)가 필요하다. 이 맥락에서 양자 회로의 성능은 단순히 개별 에너지 준위의 제곱평균제곱근(RMS) 오차만으로 판단될 수 없다. 스펙트럼 오차는 이후의 에너지 보간 및 제로 가둠 외삽(zero-confinement extrapolation) 과정에서 증폭되거나, 억제되거나, 혹은 재구성될 수 있기 때문이다. 따라서 자원 절감형 양자 회로는 단순한 스펙트럼 충실도가 아니라, 최종 산란 재구성에 필요한 정보를 보존하는 능력을 기준으로 벤치마킹되어야 한다.

방법론
저자들은 IBM Aachen 양자 프로세서의 데이터를 사용하여 탄성 s-파 중성자-양성자(np) 위상 변화를 재구성한다. 본 연구는 두 가지 RODEO(Repetitive Operation for Discrete Energy Optimization) 알고리즘 구현을 비교한다:

1. 동적 직접 R10 (Dynamic Direct R10): 고정된 가우시안 분위수 시간 단계를 사용하는 10개 사이클의 기준 구현이다. 각 쿼리 에너지마다 제어된 진화(controlled evolution)가 별도로 재합성되고 컴파일된다.
2. 인수 분해된 기하학적-시간 RODEO (FG–R6): 기하학적으로 샘플링된 시간 단계를 사용하는 6개 사이클의 압축 구현이다. 이 방법은 정확한 쿼리-위상 분리(exact query-phase separation)를 채택하며, 해밀토니안 진화는 참조 에너지($E_{ref}$)에 고정되어 한 번만 컴파일된다. 쿼리 의존적 위상($\Delta E = E_q - E_{ref}$)은 컴파일 후 보조 큐비트(ancilla qubit)에 결합된다.

두 방법 모두 IBM Aachen 백엔드 상에서 3-큐비트 레지스터(2개의 시스템 큐비트와 1개의 재사용 가능한 보조 큐비트)를 사용한다. 물리적 모델은 도식적인 s-파 np 사각 우물(square-well) 상호작용이다. 클래식하게 구성된 4차원 Okubo–Lee–Suzuki (OLS) 유효 해밀토니안이 시스템을 표현하는 데 사용된다. 실험에서는 다섯 가지 다른 트랩 강도($\omega = 3.6$ ~ $4.0$ MeV)에서 네 개의 양의 에너지 준위를 측정하여 총 20개의 스펙트럼 입력을 얻는다. 이 입력값들은 자유 공간으로 외삽하여 위상 변화를 계산하기 위한 클래식 MERE 분석에 투입된다.

주요 기여

• 인수 분해를 통한 회로 압축: 본 논문은 보조 큐비트 재사용과 인수 분해된 시간 진화 전략을 결합한 FG–R6 회로를 소개한다. 쿼리 위상을 해밀토니안 진화로부터 분리함으로써, 쿼리 에너지마다 진화 연산자를 다시 컴파일할 필요 없이 회로의 제어 진화 사이클을 10개에서 6개로 줄인다.
• 하드웨어 벤치마킹: 본 연구는 단일 준위 테스트를 넘어, 자유 공간 재구성에 필요한 전체 20개 스펙트럼 입력 세트를 다룸으로써 핵 산란에 대한 완전한 하드웨어 기반 벤치마크를 제공한다.
• 자원 대 정확도 트레이드오프 분석: 회로 깊이와 2-큐비트 게이트 수를 줄이는 것이 중간 단계의 에너지 RMS 오차를 증가시킬 수 있음에도 불구하고, 최종 산란 관측량의 정확도를 저해하지 않을 수 있음을 입증한다.

결과

• 자원 절감: 상대적으로, FG–R6 구현은 R10 구현 대비 중앙 컴파일 회로 깊이를 33.1%(459에서 307으로) 감소시키고, 2-큐비트 게이트 수를 37.2%(113에서 71로) 감소시킨다.
• 에너지 불확실성: FG–R6는 R10(0.610 keV)보다 더 큰 중앙 피팅 진폭과 참조 에너지에 대한 더 큰 중앙 RMS 중심 오프셋($\Delta E_c$, 1.648 keV vs. 0.793 keV)을 가짐에도 불구하고, 더 작은 중앙 유한 샷(finite-shot) 에너지 불확실성(0.534 keV)을 달성한다.
• 산란 위상 변화 정확도: 전체 0.1–30.0 MeV 그리드에서, 재구성된 위상 변화의 정확도(정확한 에너지 MERE 참조값 대비 최대 편차)는 R10의 0.842°에서 FG–R6의 0.163°로 크게 감소한다. 해석적 사각 우물 솔루션과의 최대 편차 또한 0.981°에서 0.156°로 떨어진다.
• 에너지 의존적 성능: 이러한 개선은 주로 저에너지 외삽 영역(10 MeV 미만)에 의해 주도된다. 10.0–30.0 MeV 공통 보간 서브셋에서는 Direct R10 구현이 FG–R6보다 약간 더 작은 잔차(residual)를 나타낸다.
• 피팅 형태 민감도: 연구는 피팅 형태에 대한 민감도가 여전히 상당함을 언급한다. 에너지 및 트랩 피팅에 대한 서로 다른 다항식 차수 간의 결과 편차가 크다는 점(예: 전체 그리드에서 R10의 경우 20.12°, FG–R6의 경우 6.75°의 반범위)은 외삽 모델의 선택이 최종 결과에 중요한 영향을 미침을 시사한다.

의의 및 주장
본 논문은 이 축소된 공간 벤치마크가 RODEO 회로 압축을 핵 연속체 관측량과 성공적으로 연결했음을 주장한다. 주요 의의는 회로의 성능이 스펙트럼 RMS 오차만이 아니라, 산란 관측량의 재구성과 그 특정 에너지 범위를 통해 평가되어야 함을 입증한 데 있다. FG–R6 구현은 테스트된 특정 스펙트럼과 윈도우에 대해 검증된 정확도-자원 타협안으로 제시되는 것이지, 보편적인 최소 사이클 솔루션이 아니다.

저자들은 이 결과가 양자 우위나 다체 확장성을 입증하는 것은 아니라고 명시적으로 밝힌다. 하드웨어 데이터는 위상 컴파일과 시간 샘플링의 효과를 개별적으로 결정하지 않으며, 독립적인 실행 재현성 또한 확립하지 않는다. 본 연구는 더 큰 모델 공간, 실제적인 상호작용, 그리고 공명 고차 부분파 산란(예: 중성자-알파 산란) 및 장거리 쿨롱 상호작용을 포함하는 계와 같은 더 복잡한 반응 역학으로 이 방법들을 확장하기 위한 출발점 역할을 한다.