[논문 리뷰] Quantum algorithms for computing observables of nonlinear partial differential equations

Shi Jin, Nana Liu|arXiv (Cornell University)|2022. 02. 16.

Model Reduction and Neural Networks인용 수 15

한 줄 요약

이 논문은 레벨 세트 방법을 통해 비선형 편미분방정식(PDE), 예를 들어 하미르톤-자코비 및 스칼라 초구형 PDE 등에서 물리적 관측가능량을 계산하기 위한 정확한 양자 알고리즘을 제시한다. 이 방법은 비선형 PDE를 선형 PDE로 매핑함으로써, 초기 자료 수 $M$에 관계없이 양자 우월성을 달성하며, 차원과 오차에 대해 지수적 속도 향상을 제공한다. 또한 $M$에 관계없이 단일 양자 계산으로 앙상블 평균을 계산할 수 있다. 이 방법은 정확하고 시간에 대해 전역적으로 유효하며, HHL을 통한 행렬 역행렬 계산을 활용해 결과 선형 시스템을 효율적으로 양자적으로 해결한다.

ABSTRACT

We construct quantum algorithms to compute physical observables of nonlinear PDEs with M initial data. Based on an exact mapping between nonlinear and linear PDEs using the level set method, these new quantum algorithms for nonlinear Hamilton-Jacobi and scalar hyperbolic PDEs can be performed with a computational cost that is independent of M, for arbitrary nonlinearity. Depending on the details of the initial data, it can also display up to exponential advantage in both the dimension of the PDE and the error in computing its observables. For general nonlinear PDEs, quantum advantage with respect to M is possible in the large M limit.

연구 동기 및 목표

• 비선형성의 정도에 관계없이 임의의 비선형성과 다수의 초기 자료를 가진 비선형 PDE의 물리적 관측가능량을 계산할 수 있는 양자 알고리즘을 개발한다.
• 비선형성 증가에 따라 성능이 급격히 악화되거나 오차를 유발하는 근사 선형화 방법의 한계를 극복한다.
• 초기 자료 수 $M$에 관계없이 계산 비용에 대한 양자 우월성을 확보한다.
• 단일 양자 시뮬레이션을 통해 $M$개의 초기 조건에 대한 앙성 평균을 효율적으로 계산할 수 있도록 한다.
• 자르기나 근사 없이 선형 표현을 통해 비선형 PDE의 정확하고 시간에 대해 전역적인 양자 시뮬레이션을 위한 프레임워크를 구축한다.

제안 방법

• 레벨 세트 형식을 활용하여 $(d+1)$차원 비선형 PDE를 최대 $(2d+1)$차원 선형 운반 PDE로 매핑함으로써 모든 물리적 정보를 정확히 유지한다.
• 초기 자료 수 $M$의 시나리오 문제를 고차원 선형 PDE에서 단일 초기 자료 문제로 축소함으로써, 한 번의 양자 실행으로 앙상블 관측가능량을 계산할 수 있도록 한다.
• 안정적인 유한차분 스킴을 사용하여 선형 PDE를 이산화하며, 안정성 조건 $D\frac{\triangle t}{h} \leq 1$을 확보함으로써 수치적 안정성을 확보한다.
• 이산화된 시스템을 큰 희소 행렬 $\mathcal{K}_{\text{ODE}}$로 표현하고, 이를 HHL 기반의 양자 행렬 역행렬에 적합한 에르미트 행렬 $\mathcal{M}_{\text{ODE}}$로 변환한다.
• HHL 알고리즘을 적용하여 $\mathcal{M}_{\text{ODE}}$를 역행렬화함으로써, 조건수 $\kappa \lesssim DNT$ 및 희소성 $s \sim D$를 가진 선형 시스템의 해를 빠르게 얻을 수 있도록 한다.
• 이산화된 해 $\Phi_{n,\boldsymbol{j}}$를 사용하여 양자 기댓값을 통해 물리적 관측가능량을 추정하며, 오차는 $\epsilon_c \leq C(Dh)^{1/3}$로 제한된다.

실험 결과

연구 질문

• RQ1 근사나 절단 없이도 비선형 PDE에 대해 양자 알고리즘이 계산 우월성을 달성할 수 있는가?
• RQ2 초기 자료 설정 수 $M$에 관계없이 관측가능량을 계산할 수 있는 비용 구조를 확보할 수 있는가?
• RQ3 레벨 세트 방법을 통해 비선형 PDE를 근사나 절단 없이 정확하고 시간에 대해 전역적인 선형 표현으로 변환할 수 있는가?
• RQ4 이 방법을 통해 비선형 PDE 관측가능량의 차원과 오차에 대해 달성 가능한 양자 속도 향상은 어느 정도인가?
• RQ5 결과 선형 시스템에 대해 HHL 기반 행렬 역행렬 계산을 적용할 때, 차원 증가에도 불구하고 양자 우월성이 유지되는가?

주요 결과

• 제안된 양자 알고리즘은 초기 자료 수 $M$에 관계없이 계산 비용이 독립적이며, 단일 양자 실행으로 $M$개의 구성에 대한 효율적인 앙상블 평균 계산을 가능하게 한다.
• 하미르톤-자코비 및 스칼라 초구형 방정식을 포함한 특정 비선형 PDE 클래스에 대해서는, 공간 차원과 오차 내성에 대해 지수적 양자 우월성을 제공한다.
• 레벨 세트 매핑은 비선형 PDE를 최대 $(2d+1)$차원 선형 PDE로 변환하여 근사나 절단 없이 정확한 표현을 보장한다.
• 물리적 관측가능량 추정 오차는 $\epsilon_c \leq C(Dh)^{1/3}$로 제한되며, $h$는 공간 그리드 크기이고 $D$는 공간 차원 수이다.
• 양자 우월성에 핵심적인 역할을 하는 행렬 역행렬 단계는 조건수 $\kappa \lesssim DNT$ 및 희소성 $s \sim D$를 가진 에르미트 행렬 $\mathcal{M}_{\text{ODE}}$를 대상으로 HHL을 적용하여 수행되며, 양자 우월성이 유지된다.
• 이 방법은 시간에 대해 전역적으로 유효하며, 비선형성의 정도에 관계없이 물리적 정확성을 유지하여 카를레만 선형화 및 기타 근사 기법의 한계를 극복한다.