[논문 리뷰] Quantum simulation of Maxwell's equations via Schrödingersation

Shi Jin, Nana Liu외 1명|arXiv (Cornell University)|2023. 08. 16.

Quantum Computing Algorithms and ArchitectureComputer Science인용 수 3

한 줄 요약

이 논문은 슈뢰딩거화 방법을 사용하여 맥스웰 방정식을 양자 알고리즘으로 시뮬레이션하는 방법을 제시한다. 이 방법은 비유니터리 전자기 동역학을 확장된 공간에서 왜곡된 위상 변환을 통해 유니터리 진동으로 변환한다. 이 접근법은 물리적 경계 조건과 인터페이스 조건을 정확하게 처리할 수 있으며, 이는 이산(예: 예이 및 업윈드 스킴) 및 연속 변수 양자 프레임워크에서 고정밀도를 달성한다. 수치 결과는 QLA 및 정확한 해석적 해와 일치한다.

ABSTRACT

We present quantum algorithms for electromagnetic fields governed by Maxwell's equations. The algorithms are based on the Schrödingersation approach, which transforms any linear PDEs and ODEs with non-unitary dynamics into a system evolving under unitary dynamics, via a warped phase transformation that maps the equation into one higher dimension. In this paper, our quantum algorithms are based on either a direct approximation of Maxwell's equations combined with Yee's algorithm, or a matrix representation in terms of Riemann-Silberstein vectors combined with a spectral approach and an upwind scheme. We implement these algorithms with physical boundary conditions, including perfect conductor and impedance boundaries. We also solve Maxwell's equations for a linear inhomogeneous medium, specifically the interface problem. Several numerical experiments are performed to demonstrate the validity of this approach. In addition, instead of qubits, the quantum algorithms can also be formulated in the continuous variable quantum framework, which allows the quantum simulation of Maxwell's equations in analog quantum simulation.

연구 동기 및 목표

• 표준 양자 시뮬레이션에서 비유니터리 동역학으로 인해 어려운 물리적 경계 및 인터페이스 조건을 포함하는 맥스웰 방정식에 대한 양자 알고리즘 개발.
• 이전에 미분방정식에 적용된 바 있는 슈뢰딩거화 방법을 맥스웰 방정식에 의해 지배되는 전자기장 시스템으로 확장.
• 양자 시뮬레이션 프레임워크 내에서 발산이 없는 자속밀도 및 총 전자기 에너지와 같은 핵심 물리적 성질의 유지 검증.
• 근접한 양자 장치에서의 접근성을 높이기 위해 이산(큐비트 기반) 및 연속 변수 양자 계산 플랫폼 모두에서의 양자 시뮬레이션 가능화.
• 완전 도체, 임피던스 및 비균일 매질을 포함한 물리적 인터페이스가 존재하는 시스템에 대해 해석적 및 기준 수치 해와의 비교를 통한 수치 실험을 통한 방법 검증.

제안 방법

• 슈뢰딩거화 방법은 보조 공간 차원을 도입하고 왜곡된 위상 변환을 통해 맥스웰 방정식을 유니터리 진동 문제로 변환함으로써 비유니터리 동역학을 유니터리 동역학으로 전환한다.
• 이 방법은 여덟 차원 행렬 표현을 가진 리만-실베르슈타인 벡터 표현 방식과 표준 전기장/자기장 표현 방식에 모두 적용된다.
• 공간 이산화를 위해 예이 스킴과 업윈드 스킴이 사용되며, 후자는 슈뢰딩거화 적용 시 일계 공간 수렴성을 보이며 양호한 안정성을 확보한다.
• 경계 조건(주기적, 완전 도체, 임피던스)은 매트릭스 표현에 유니터리 변환을 적용하여 구현되며, 물리적 제약 조건을 유지한다.
• 비균일 매질에서의 인터페이스 조건은 잠기고 있는 표면 방법(IIM)을 사용하여 유전율 및 투자율의 불연속성에 대해 정확하게 처리한다.
• 이 프레임워크는 연속 변수 양자 시스템으로 확장되어 공간적으로 변화하는 파동 속도로 인한 조밀한 해밀토니안 행렬을 피함으로써, 맥스웰 방정식의 더 효율적인 아날로그 양자 시뮬레이션을 가능하게 한다.

실험 결과

연구 질문

• RQ1 슈뢰딩거화 방법이 맥스웰 방정식에 효과적으로 적용되어 물리적 경계 조건을 갖춘 유니터리 양자 시뮬레이션을 가능하게 할 수 있는가?
• RQ2 슈뢰딩거화는 수치적 스킴에서 발산이 없는 자속밀도 및 총 전자기 에너지와 같은 기본 물리적 성질의 유지에 어떤 영향을 미치는가?
• RQ3 비균일 매질에서 유전율 및 투자율이 불연속적인 인터페이스를 통해 전자기파의 전파를 얼마나 정확하게 시뮬레이션할 수 있는가?
• RQ4 슈뢰딩거화와 결합했을 때 예이 스킴과 업윈드 스킴 간의 정밀도 및 안정성은 어떻게 비교되는가?
• RQ5 슈뢰딩거화 접근법은 연속 변수 양자 시스템에 어떻게 적응시켜 전자기장의 아날로그 양자 시뮬레이션을 가능하게 할 수 있는가?

주요 결과

• 수치 실험을 통해 슈뢰딩거화 방법은 예이 스킴과 업윈드 스킴 모두에서 자속밀도의 발산이 없는 조건과 전자기장의 총 에너지를 정확히 유지함을 입증하였다.
• 수치 시뮬레이션 결과, 동일한 메쉬 크기에서 슈뢰딩거화 방법을 적용한 예이 스킴은 업윈드 스킴보다 더 높은 정밀도를 보였으며, 후자는 일계 공간 수렴성을 보였다.
• 이 방법은 비균일 매질 인터페이스에서 전자기파의 반사 및 투과를 정확히 포착하였으며, ε₂=2, μ₂=2 인 유전체 인터페이스에 대해 정확한 해석적 해와 매우 유사한 결과를 도출하였다.
• 진공(ε₁=1, μ₁=1)에서 유전체 매질(ε₂=3, μ₂=1)로 전파되는 가우시안 펄스에 대해 슈뢰딩거화 방법은 미세 메쉬(M=2¹¹)를 사용한 양자 격자 알고리즘(QLA)의 결과와 잘 일치하는 해를 생성하였다.
• 완전 도체 및 임피던스 경계 조건은 슈뢰딩거화 프레임워크 내에서 유니터리 변환을 통해 구현되었으며, 물리적 일관성을 유지하고 기대되는 전자기장 행동과 일치하였다.
• 연속 변수 양자 프레임워크에서 이 방법은 공간적으로 변화하는 파동 속도로 인한 조밀한 해밀토니안 행렬을 피함으로써, 맥스웰 방정식의 더 효율적인 아날로그 양자 시뮬레이션을 가능하게 하였다.