[論文レビュー] Quantum simulation of Maxwell's equations via Schrödingersation

Shi Jin, Nana Liu+1 more|arXiv (Cornell University)|Aug 16, 2023

Quantum Computing Algorithms and ArchitectureComputer Science被引用数 3

ひとこと要約

本稿では、非ユニタリな電磁気的ダイナミクスを拡張空間における歪みのある位相変換によりユニタリな時間発展に変換する『シュレーディンガー化』法を用いて、マクスウェル方程式の量子アルゴリズムを提示する。この手法により、物理的境界条件および界面条件を正確に再現できる電磁界のシミュレーションが可能となり、離散的（Yee法および上流化スキーム）および連続変数量子フレームワークの両方で高い忠実度が達成された。数値結果は、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.

研究の動機と目的

• 標準的な量子シミュレーションでは非ユニタリなダイナミクスのため困難である物理的境界条件および界面条件を組み込んだマクスウェル方程式のための量子アルゴリズムを開発すること。
• 従来、ODEおよびPDEに適用されてきたシュレーディンガー化法を、マクスウェル方程式に従う電磁界系へと拡張すること。
• 量子シミュレーションフレームワーク内で、磁束密度の発散なしの性質および全電磁エネルギーの保存といった、重要な物理的性質が保持されることを示すこと。
• 近い将来の装置に適したアクセシビリティを高めるために、離散的（キュービットベース）および連続変数量子コンピューティングプラットフォームの両方で量子シミュレーションを可能にすること。
• 完全導体、インピーダンス、不均質媒質を含む界面を有する系について、解析解および基準数値解と比較して、数値実験を通じて手法の妥当性を検証すること。

提案手法

• シュレーディンガー化法により、非ユニタリなダイナミクスをユニタリな時間発展問題に変換するため、歪みのある位相変換を用いて補助的な空間次元を導入する。
• この手法は、リーマン＝シルバーマンベクトル表現（8次元行列表現）および標準的な電場・磁場表現の両方へ適用される。
• 空間離散化にはYeeスキームおよび上流化スキームが用いられ、後者はシュレーディンガー化においても第一階層の空間精度を示しながら良好な安定性を示す。
• 境界条件（周期的、完全導体、インピーダンス）は、行列表現上でのユニタリ変換により実装され、物理的制約が保持される。
• 不均質媒質における界面条件は、浸漬界面法（IIM）を用いて処理され、誘電率および透磁率の不連続性が正確に扱われる。
• 連続変数量子系へとフレームワークを拡張し、空間的に変化する波動速度に起因する密度の高いハミルトニアン行列を回避することで、アナログ量子シミュレーションがより効率的に行えるようになった。

実験結果

リサーチクエスチョン

• RQ1 シュレーディンガー化法は、マクスウェル方程式に効果的に適用可能であり、物理的境界条件を満たすユニタリな量子シミュレーションを可能にするか？
• RQ2 シュレーディンガー化は、数値スキームにおいて、磁束密度の発散なしの性質や全電磁エネルギーの保存といった、基本的な物理的性質をどのように保持するか？
• RQ3 不連続な誘電率および透磁率を有する界面を有する媒質を通過する電磁波の伝搬を、この手法はどの程度正確にシミュレートできるか？
• RQ4 Yeeスキームと上流化スキームをシュレーディンガー化と組み合わせた場合、量子シミュレーションにおける精度および安定性において、両者の相違は何か？
• RQ5 シュレーディンガー化アプローチは、連続変数量子系に適応可能であり、電磁界のアナログ量子シミュレーションを可能にするか？

主な発見

• 数値実験により、Yeeスキームおよび上流化スキームの両方において、シュレーディンガー化法が磁束密度の発散なしの条件および電磁界全エネルギーの保存を正確に保持していることが確認された。
• 数値シミュレーションの結果、同じメッシュサイズにおいて、シュレーディンガー化法とYeeスキームを組み合わせた場合、上流化スキームよりも高い精度が得られ、後者は第一階層の空間収束を示した。
• この手法は、誘電率ε₂=2、透磁率μ₂=2の誘電体界面における電磁波の反射・透過を正確に捉えており、正確な解析解と密接に一致した。
• 真空中（ε₁=1、μ₁=1）から誘電体媒質（ε₂=3、μ₂=1）へと進行するガウスパルスについて、シュレーディンガー化法により得られた解は、細かいメッシュ（M=2¹¹）を用いた量子格子アルゴリズム（QLA）の結果とよく一致した。
• シュレーディンガー化フレームワーク内での完全導体およびインピーダンス境界条件のユニタリ変換による実装は、物理的整合性を保ち、期待される場の挙動と一致した。
• 連続変数量子フレームワークにおいて、空間的に変化する波動速度に起因する密度の高いハミルトニアン行列を回避でき、マクスウェル方程式のより効率的なアナログ量子シミュレーションが可能になった。