[論文レビュー] GPU accelerated Trotter-Suzuki solver for quantum spin dynamics

Axel D. Dente, Carlos Bederián+2 more|arXiv (Cornell University)|Apr 30, 2013

Quantum many-body systems参考文献 16被引用数 6

ひとこと要約

この論文では、多数のスピン系における非平衡量子スピンダイナミクスをシミュレートするためのGPUアクセラレートド4次Trotter-Suzukiアルゴリズムを提示する。スピン相互作用と局所的回転の並列処理を活用することで、CPUベースのOpenMPコードに比べ最大30倍の高速化を達成し、実用的なメモリ制限内で最大27スピンのシミュレーションが可能になった。

ABSTRACT

The resolution of dynamics in out of equilibrium quantum spin systems lies at the heart of fundamental questions among Quantum Information Processing, Statistical Mechanics and Nano-Technologies. Efficient computational simulations of interacting many-spin systems are extremely valuable tools for tackling such questions. Here, we use the Trotter-Suzuki (TS) algorithm, a well-known strategy that provides the evolution of quantum systems, to address the spin dynamics. We present a GPU implementation of a particular TS version, which has been previously implemented on single cores in CPUs. We develop a massive parallel version of this algorithm and compare the efficiency between CPU and GPU implementations. This method reduces the execution time considerably and is capable of dealing with systems of up to 27 spins (only limited by GPU memory).

研究の動機と目的

• 相互作用する多数スピン量子系における非平衡ダイナミクスをシミュレートする計算的課題に対処すること。
• 従来のCPUベースのシミュレーションに制限をもたらす指数的増加するヒルベルト空間スケーリングを克服すること。
• 量子状態の精度を保つために、ユニタリで断片化しない数値的手法を構築すること。
• 対称性の仮定なしに、スピンハミルトニアンのTrotter-Suzuki分解をGPUの並列処理で高速化すること。
• GPUメモリの制限のみに依存して、最大27スピンのスピン系のスケールアップシミュレーションを可能にすること。

提案手法

• スピン1/2ハミルトニアンの時間発展演算子に対して4次Trotter-Suzuki分解を実装する。
• XYZ分解を用いてハミルトニアンをX、Y、Z成分に分割し、Z項は対角で、X/Y項は非対角となるようにする。
• 非対角のXおよびY項を効率的な位相操作が可能な有効なZ項に変換するための局所的回転を適用する。
• GPUカーネルにおけるスレッドレベルの独立性を活用し、位相更新と局所的回転を大規模並列処理で実装する。
• 断片化なしに、$2^N$個のイジング状態からなる計算基底を用いて多体波動関数を表現する。
• 時間発展演算中における状態ベクトルのノルムを保存することでユニタリティを保証し、ルンゲ・クッタのような非ユニタリな積分法を避ける。

実験結果

リサーチクエスチョン

• RQ1 GPUアクセラレーションにより、Trotter-Suzuki法を用いた多数スピン量子系の時間発展演算の実行時間を著しく短縮できるか？
• RQ2 非自明なスピンハミルトニアンに対して、GPUアーキテクチャ上でTrotter-Suzukiアルゴリズムをどの程度並列化できるか？
• RQ3 GPUアクセラレートドTSを用いて、完全なヒルベルト空間表現で効率的にシミュレート可能な最大系サイズ（スピン数）は何か？
• RQ4 GPU最適化されたTS実装の精度は、正確なまたは基準となるCPUベースのシミュレーションと比較してどの程度か？
• RQ5 このGPUアクセラレートドTSソルバは、閉じた量子系における熱化や行列積状態の時間発展といった物理的に意味のある状況に効果的に応用可能か？

主な発見

• GPU実装は、6コアプロセッサ上で最適化されたOpenMP CPU実装に比べ、最大で約30倍の高速化を達成した。
• 本手法により、GPUメモリ容量（テスト環境では約6 GB）の制限のみに依存して最大27スピンのシミュレーションが可能になった。
• Trotter-Suzuki近似の妥当性を検証するためのロシュミルトエコーを用いた時間発展演算の相対誤差は、許容範囲内に保たれた。
• ロシュミルトエコーは、GPUおよびCPU実装の両方における計算誤差を検出する信頼できる指標となった。
• 本実装は、ランダム重ね合わせ状態や行列積状態といった高度な量子状態表現とも互換性があり、量子多体物理学における広範な応用を可能にした。
• アルゴリズムは時間発展演算の全過程でユニタリティを維持しており、確率の保存と物理的整合性を保証した。