[논문 리뷰] Efficient quantum algorithm for nonlinear reaction-diffusion equations and energy estimation

Dong An, Di Fang외 4명|arXiv (Cornell University)|2022. 05. 02.

Quantum Computing Algorithms and Architecture인용 수 11

한 줄 요약

이 논문은 비선형 반응-확산 방정식을 해결하기 위한 효율적인 양자 알고리즘을 제시한다. 개선된 Carleman 선형화 기법을 활용하여 ℓ∞ 노름을 사용해 더 온건한 수렴 조건 R_D < 1를 정의함으로써, 공간 그리드 크기와 무관하게 수렴 조건을 확보한다. 이 방법은 고차원 PDE를 해결하는 데 양자 속도 향상을 가능하게 하며, 해 상태의 도함수 정보로부터 운동 에너지를 추정하기 위한 후처리를 허용한다.

ABSTRACT

Nonlinear differential equations exhibit rich phenomena in many fields but are notoriously challenging to solve. Recently, Liu et al. [1] demonstrated the first efficient quantum algorithm for dissipative quadratic differential equations under the condition $R &lt; 1$, where $R$ measures the ratio of nonlinearity to dissipation using the $\ell_2$ norm. Here we develop an efficient quantum algorithm based on [1] for reaction-diffusion equations, a class of nonlinear partial differential equations (PDEs). To achieve this, we improve upon the Carleman linearization approach introduced in [1] to obtain a faster convergence rate under the condition $R_D &lt; 1$, where $R_D$ measures the ratio of nonlinearity to dissipation using the $\ell_{\infty}$ norm. Since $R_D$ is independent of the number of spatial grid points $n$ while $R$ increases with $n$, the criterion $R_D&lt;1$ is significantly milder than $R&lt;1$ for high-dimensional systems and can stay convergent under grid refinement for approximating PDEs. As applications of our quantum algorithm we consider the Fisher-KPP and Allen-Cahn equations, which have interpretations in classical physics. In particular, we show how to estimate the mean square kinetic energy in the solution by postprocessing the quantum state that encodes it to extract derivative information.

연구 동기 및 목표

• 물리학, 생물학, 데이터 과학 등 광범위한 응용을 갖는 PDE의 한 종류인 비선형 반응-확산 방정식을 해결하기 위한 효율적인 양자 알고리즘을 개발하기 위해.
• 이전의 양자 알고리즘이 ℓ2 노름에 기반한 엄격한 R < 1 조건을 요구하여 그리드 세분화에 따라 악화되는 한계를 극복하기 위해.
• Carleman 선형화 접근법을 향상시켜 ℓ∞-노름 기반 수렴 조건 R_D < 1를 도입함으로써, 공간 그리드 세분화에 대해 안정성을 확보하기 위해.
• 양자 해 상태의 후처리를 통해 평균 제곱 운동 에너지와 같은 고전적 관측량을 계산할 수 있도록 하기 위해.
• 암호화된 도함수 준비 절차를 통해 양자 해 상태에서 도함수 정보를 활용하여 평균 제곱 운동 에너지를 추정하는 양자 알고리즘의 가능성을 탐색하기 위해.

제안 방법

• 비선형 반응-확산 방정식을 무한차원 선형 시스템으로 매핑하는 Carleman 선형화 프레임워크를 적응하여, 원래 PDE의 구조를 유지하기 위해.
• 해의 ℓ∞ 노름에 기반한 새로운 수렴 조건 R_D < 1을 도입하여, 이는 공간 그리드 세분화에 대해 유한하게 유지되며, ℓ2 기반의 R < 1 조건과는 달리 그렇지 않다.
• 양자 위상 추정과 앰플리티 증폭을 사용하여 복잡도가 poly(log n, T, 1/ε)인 선형화된 시스템의 해를 코딩하는 양자 상태를 준비하기 위해.
• 암호화된 도함수 준비 절차를 통해 도함수 정보를 포함하는 양자 상태에서 평균 제곱 운동 에너지와 같은 고전적 정보를 추출하기 위해 앰플리티 추정을 적용하기 위해.
• 최대 원리와 해의 정규성에 기반하여 ℓ∞ 의미에서 해의 유한성과 안정성을 확보함으로써, 새로운 수렴 조건을 가능하게 하기 위해.
• 해의 ℓ2 노름에 대한 의존도를 최소화함으로써 그리드 세분화에 대해 효율성을 유지하는 양자 알고리즘 설계하기 위해.

(a) $f(u)=0.4u-u^{2}$ and $u_{0}(x)=0.1-0.1\cos(4\pi x)$ .

실험 결과

연구 질문

• RQ1 PDE에서 공간 그리드 크기와 무관하게 수렴 조건을 개선할 수 있는가?
• RQ2 비선형성 대 소산 비율을 측정하는 데 ℓ2 노름 대신 ℓ∞ 노름을 사용하면, 반응-확산 방정식에 대한 더 강건하고 확장 가능한 양자 알고리즘이 될 수 있는가?
• RQ3 기하급수적 오버헤드 없이, 양자 해 상태에서 고전적 관측량(예: 운동 에너지)을 효율적으로 추정할 수 있는가?
• RQ4 암호화된 도함수 준비와 앰플리티 추정을 사용하여, 비선형 PDE의 고전적 출력에 대해 초다항 속도 향상을 달성할 수 있는가?
• RQ5 알고리즘의 성능은 그리드 세분화에 따라 어떻게 변화하는가? 그리고 공간 이산화에 대해 완전히 민감하지 않게 만들 수 있는가?

주요 결과

• 알고리즘은 공간 그리드 포인트 수 n과 무관하게 유지되는 ℓ∞ 노름 기반 수렴 조건 R_D < 1를 달성하여, ℓ2 기반 R < 1 조건보다 훨씬 온건한 조건을 제공한다.
• 복잡도는 O(T²q/(gε) · poly(log T, log n, log 1/ε))로 표현되며, 비선형 ODE에 대한 이전 양자 알고리즘 대비 시간 복잡도에서 제곱 개선을 이룬다.
• 도함수 정보를 코딩하는 양자 상태를 준비하고 앰플리티 추정을 적용함으로써 평균 제곱 운동 에너지를 효율적으로 추정할 수 있다.
• ℓ∞-노름 기반 조건 덕분에 그리드 세분화에 대해 안정성과 수렴성을 유지하여, 이전 접근법의 핵심 한계를 해결한다.
• Fisher-KPP 및 Allen-Cahn과 같은 물리적으로 의미 있는 방정식에 적용 가능하며, 인구 역학 및 상 분리 현상을 모델링한다.
• 고차원 비선형 PDE를 해결하는 데서 양자 우위를 달성할 수 있는 길을 제시하며, 후처리를 통해 고전적 출력을 효율적으로 추출할 수 있다는 전제 하에.

(b) $f(u)=0.16u-u^{3}$ and $u_{0}(x)=0.1\sin(4\pi x)$ .