[論文レビュー] QSystem: bitwise representation for quantum circuit simulations

Evandro C. R. Rosa, Bruno G. Taketani|arXiv (Cornell University)|Apr 7, 2020

Quantum Computing Algorithms and Architecture参考文献 3被引用数 4

ひとこと要約

本稿では、限られた重ね合わせを示す量子状態のシミュレーションにおいて指数的高速化を達成する、新しいビット単位の表現に基づくオープンソースの量子回路シミュレータであるQSystemを紹介する。ハッシュマップとビット演算を用いた新規表現により、非ゼロの基底状態のみを効率的に格納・操作することで、Qiskit、Cirq、Forest SDKを凌駆し、特にスパースな量子状態のシナリオにおいて顕著な性能向上を示した。事前に状態構造を把握していなくても、性能向上が達成可能であることを実証した。

ABSTRACT

We present QSystem, an open-source platform for the simulation of quantum circuits focused on bitwise operations on a Hashmap data structure storing quantum states and gates. QSystem is implemented in C++ and delivered as a Python module, taking advantage of the C++ performance and the Python dynamism. The simulators API is designed to be simple and intuitive, thus streamlining the simulation of a quantum circuit in Python. The current release has three distinct ways to represent the quantum state: vector, matrix, and the proposed bitwise. The latter constitutes our main results and is a new way to store and manipulate both states and operations which shows an exponential advantage with the amount of superposition in the systems state. We benchmark the bitwise representation against other simulators, namely Qiskit, Forest SDK QVM, and Cirq.

研究の動機と目的

• 限られた重ね合わせを示す量子状態を効率的に扱う量子回路シミュレータの開発を目的とし、メモリおよび時間コストの低減を図る。
• 従来の量子シミュレータが示す指数的リソーススケーリングの問題を、ビット演算に基づく新規データ表現を導入することで解決することを目的とする。
• C++の高速性とPythonの使いやすさを統合した、高性能でオープンソースのシミュレーションプラットフォームを提供することを目的とする。
• さまざまな量子回路ワークロード下で、新規ビット表現を既存のシミュレータと比較してベンチマークすることを目的とする。
• 関連する基底状態の事前知識がなくても、指数的高速化が達成可能であることを実証することを目的とする。

提案手法

• ビット単位の表現は、量子状態をハッシュマップとして格納し、キーを計算基底状態のインデックス（整数）とし、値を複素振幅とする。これにより、スパースな格納が可能となる。
• 量子ゲートは、XORなどのビット論理演算（CNOT用）およびビットシフトを用いて実装され、行列演算を一切行わずに状態インデックスを直接操作できる。
• 性能を確保するためC++コアを採用し、使いやすさを考慮してPython APIを提供することで、量子コンputation用語を用いた直感的な回路構築を可能にする。
• 3種類の表現形式をサポートする：ベクトル（密度状態）、行列（誤差モデル用密度行列）、および提案されたビット単位（ビット演算を用いたスパース状態）。
• 測定プロセスは、各キュービットに対してアルゴリズム4を逐次適用し、ハッシュマップからの確率的サンプリングにより状態の崩壊と測定結果の抽出を実行する。
• 任意の基底状態をサポートするが、現在の実装ではパウリZ演算子によって定義される計算基底を用いているが、行列積状態で用いられる他の基底へも拡張可能である。

実験結果

リサーチクエスチョン

• RQ1 ハッシュマップとビット演算を用いたビット単位の表現が、限られた重ね合わせを示す状態において、量子回路シミュレーションで指数的高速化を達成できるか？
• RQ2 回路の複雑さや状態のスパarsityに応じて、Qiskit、Cirq、Forest SDKといった既存のシミュレータと比較して、ビット単位の表現の性能はどのように変化するか？
• RQ3 観察された高速化は、基底状態の事前知識に依存するのか、それとも関連する状態が未知であっても有効に機能するのか？
• RQ4 高エンタングルメントを示す状態（例：GHZ状態）の準備においても、ビット単位のアプローチが指数的スケーリングの利点を維持できるか？
• RQ5 この手法は、計算基底を超えた他のデータ構造や基底集合へも一般化可能か、その範囲はどの程度か？

主な発見

• 理論的モデルが予測した通り、わずかで未知の基底状態集合を必要とする量子状態において、ビット単位の表現がシミュレーション時間において指数的高速化を達成した。
• GHZ状態の準備に関するベンチマークでは、QSystemのビット単位表現がQiskit、Cirq、Forest SDKを大きく上回り、キュービット数の増加に伴う時間の増加が著しく抑制された。
• 測定処理において、同時に測定する最適化を施していなくても、Cirqに比べてオーバーヘッドを1桁以上削減した。
• 密度状態のシミュレーションでは高速化が予想されない場合でも、指数的スケーリングの利点が観察されたことから、あらゆる状態タイプにわたる頑健性が示された。
• 関連する基底状態の知識がなくても、状態表現のスパarsityにのみ依存して性能向上が達成された。
• エンタングルド状態（例：GHZ状態）のシミュレーションにおいても、高い効率性を示した。これにより、エンタングルメント自体が高速化を制限するのではなく、状態のスパarsityがその鍵であることが実証された。