Relativistic Quantum Simulation under Periodic and Dirichlet Boundary Conditions: A First-Quantised Framework for Near-Term Devices

기술 요약: 주기적 및 디리클레 경계 조건 하에서의 상대론적 양자 시뮬레이션

문제 정의
양자 시스템을 시뮬레이션하는 것은 시스템 크기에 따라 복잡성이 기하급수적으로 증가하기 때문에 계산 집약적이며, 물리적 현상을 정확하게 설명하기 위해 상대론적 효과를 포함해야 할 경우 이러한 어려움은 더욱 가중됩니다. 섭동 이론은 비상대론적 프레임워크를 넘어 상대론적 효과를 연구하기 위한 표준적인 접근법이지만, 이를 근미래 양자 장치에 구현하는 것은 여전히 어렵습니다. 기존 방법들은 대개 양자 화학에서 표준적인 제2 양자화에 의존하지만, 이는 비선형성을 가진 시스템을 직관적으로 처리하는 데 부족할 수 있습니다. 따라서 서로 다른 경계 조건 하에 있는 시스템의 상대론적 바닥 상태 에너지를 근미래 하드웨어에서 추정할 수 있는 실용적인 제1 양자화 프레임워크가 필요합니다.

방법론
저자들은 시스템 큐비트로 표현되는 유한한 격자 상에 파동 함수를 이산화하는 상대론적 양자 시뮬레이션을 위한 제1 양자화 프레임워크를 제안합니다. 핵심 방법론은 다음 구성 요소들을 포함합니다:

1. 이산화 및 이동 연산자 (Discretization and Translation Operators):
   연속적인 1차원 공간 파동 함수는 $L$개의 큐비트를 사용하여 $2^L$개의 점을 가진 이산 격자로 매핑되며, 이때 위치 $x$는 도메인 $0 \le x < 1$로 스케일링됩니다. 제곱 운동량 연산자 $\hat{p}^2$은 양자 이동 연산자에 기반한 유한 차분법을 사용하여 근사됩니다. 이동 연산자 $\hat{A}_x$(양자 가산기)는 위치 상태 $|j\rangle$를 $|j+\delta x\rangle$로 이동시키며, 그 켤레인 $\hat{A}_x^\dagger$(양자 감산기)는 반대 방향으로 이동시킵니다. 제곱 운동량은 $\hat{p}^2 = -(mc L_m)^2 (\hat{A}^{(1)} - 2\hat{1})$로 표현됩니다. 여기서 $L_m = \lambda_m / \delta x$는 축소된 컴프턴 파장과 관련된 무차원 매개변수입니다.

2. 운동 에너지의 섭동 전개 (Perturbative Expansion of Kinetic Energy):
   상대론적 운동 에너지 해밀토니안 $\hat{H}_{re}^K$는 상대론적 에너지-운동량 관계의 제곱근 형태로부터 유도됩니다. 이는 $\langle \hat{p}^2 \rangle \ll m^2 c^2$인 영역에서 유효한 고차 운동량 항($\hat{p}^2, \hat{p}^4, \hat{p}^6, \dots$)의 섭동 급수로 전개됩니다. 전체 상대론적 운동 에너지는 클라인-고든 방정식에서 유도된 계수들을 포함하여 이 급수를 절단함으로써 근사됩니다.

3. 경계 조건 처리 (Boundary Condition Handling):
   본 프레임워크는 두 가지 뚜렷한 경계 조건을 다룹니다:

• 주기적 경계 조건 (Periodic Boundary Conditions, PBC): 이동 연산자는 격자를 감싸도록 하여 유니타리성을 보존합니다 ($\hat{A}_x |2^L-1\rangle = |0\rangle$). 짝수 차수 운동량의 기댓값은 $\langle \hat{A}^{(l)} \rangle$의 조합으로 표현됩니다.
• 디리클레 경계 조건 (Dirichlet Boundary Conditions, DBC): 여기서는 파동 함수가 감싸지지 않습니다 ($\psi(0) \neq \psi(1)$). 저자들은 경계를 처리하기 위해 보정 연산자 $\hat{E}_0 = |2^L-1\rangle\langle 0| + |0\rangle\langle 2^L-1|$를 도입합니다. DBC를 위한 기댓값은 PBC 항에 격자 경계에서의 확률 진폭에 의존하는 관련 연산자 $\hat{E}_0, \hat{E}_1, \hat{E}_2$를 포함하는 추가 항들을 더하여 공식화됩니다.

4. 양자 시뮬레이션 알고리즘:
   저자들은 필요한 기댓값을 추정하기 위한 양자 회로를 설계합니다:

• 운동 에너지: 제어된 이동 게이트(controlled-translation gates)가 제어 큐비트에 반복적으로 적용됩니다. 제어 큐비트를 Z-기저에서 측정하면 $\langle (\hat{A}_x^\dagger)^l \rangle$의 실수부에 해당하는 값을 얻으며, 이는 필요한 운동량 항에 대응합니다.
• 퍼텐셜 에너지: 퍼텐셜을 나타내는 밀도 행렬 $\hat{\rho}_V$를 사용하여 시스템 상태의 확률을 통해 $\langle \hat{V} \rangle$를 추정합니다.
• 변분 최법 (Variational Optimization): 전체 에너지는 안사츠(ansatz) 상태 $|\psi(\{\bar{\kappa}\}, x)\rangle$에 대한 변분 매개변수 세트 $\{\bar{\kappa}\}$에 대해 운동 에너지와 퍼텐셜 에너지 기댓값의 합을 최소화함으로써 추정됩니다.

주요 기여

• 제1 양자화 레시피: 본 논문은 제1 양자화를 사용한 상대론적 시뮬레이션을 위한 구체적인 레시피를 제공하며, 이는 널리 사용되는 제2 양자화 방법의 대안을 제시합니다. 이는 특히 비선형성을 가진 시스템에 유익합니다.
• 통합된 경계 프레임워크: 주기적 및 디리클레 경계 조건을 모두 처리하는 통합된 접근 방식을 제시합니다. 결정적으로, PBC에서 얻은 데이터를 특정 경계 보정 항을 더함으로써 DBC 에너지를 추정하는 데 재사용할 수 있음을 보여주며, 이를 통해 계산 오버헤드를 줄입니다.
• 근미래 장치 호환성: 본 방법은 표준 양자 회로 프리미티브(제어된 가산기, 하다마르 게이트 및 측정)와 변분 최적화를 활용하므로 근미래 양자 장치에 구현하기 적합합니다.
• 명시적 회로 설계: 저자들은 이동 연산자의 기댓값과 섭동 전개에 필요한 경계 보정 항을 계산하기 위한 상세한 회로 설계를 제공합니다.

결과 및 주장
본 논문은 근미래 양자 장치에서 상대론적 효과를 시뮬레이션하기 위한 "실용적인 경로"를 제공한다고 주장합니다. 이론적 결과는 다음을 보여줍니다:

• 특정 차수까지의 상대론적 운동 에너지 항은 이동 연산자($\langle \hat{A}^{(l)} \rangle$)와 경계 연산자($\langle \hat{E}_d \rangle$)의 선형 결합으로 분해될 수 있습니다.
• 이러한 전개를 위한 계수들은 물리적 매개변수($m, c, \hbar, \delta x$)를 기반으로 사전 계산될 수 있습니다.
• 본 방법은 안사츠 상태를 변분 최적화하여 비상대론적 및 상대론적 바닥 상태 에너지를 모두 추정할 수 있게 합니다.

의의
저자들은 이 연구를 미래의 양자 물리학 및 화학 발전을 지원하기 위한 기초적인 단계로 설정합니다. 근미래 하드웨어에서 상대론적 효과를 시뮬레이션할 수 있게 함으로써, 이 프레임워크는 현재의 양자 시뮬레이션 능력의 공백을 메웁니다. 본 논문은 격자 해상도를 높이기 위해 더 많은 큐비트를 늘릴 수 있지만, 섭동 전개가 유효하게 유지되도록 무차원 매개변수 $L_m$이 충분히 작게 유지되도록 주의해야 함을 강조합니다. 이 작업은 완전한 디락 방정식을 해결하거나 초상대론적 극한을 다루는 것을 목표로 하지 않으며, 특히 비상대론적 프레임워크에 대한 상대론적 보정이 유효한 영역에 초점을 맞추고 있습니다. 향후 연구로는 더 높은 차수의 라플라시안 연산자와 더 정확한 안사츠 구축 방법을 탐구하는 것이 제안되었습니다.