技術要約: QISKIT-QUDITS

問題提起

$d$ 個の計算レベルを持つ、2レベルの量子ビット（qubit）を一般化したクディット（qudit）は、量子情報処理において潜在的な利点を提供します。しかし、Qiskit フレームワークにはクディットに対するネイティブなサポートが欠けており、既知のサードパーティ製拡張によるギャップとして扱われています。既存のツールチェーンは二分化に直面しています。すなわち、ネイティブなクディット・シミュレータ（例：Cirq, QuTiP）は Qiskit の広範なコンパイル・スタック、トランスパイラ・パス、およびハードウェア・バックエンドを活用できない一方で、量子ビットを対象としたツールは、クディットレベルの操作のためのプログラミング・サーフェスを欠いています。ユーザーは現在、エンコーディング・スキーム、ゲート行列、分解、および古典的デコーディングを、物理的に正しくない状態空間の扱いに関するエラーが発生しやすいプロセスとして、手動で実装しなければなりません。

手法

本論文では、各クディットを $m = \lceil \log_2 d \rceil$ 個の量子ビットにエンコードすることで、$d$ レベルのクディットをシミュレートする Qiskit 拡張機能 qiskit-qudits を導入します。このソフトウェアは、既存のエコシステムとの互換性を確保するため、Qiskit の公開された Instruction、Gate、および ControlledGate クラスに基づいて構築されています。

コア設計決定

1. エンコーディングとアイデンティティ・パディング:
   $d$ が 2 の累乗であるとは限らないため、エンコードされたヒルベルト空間 ($D = 2^N$) は、通常、論理部分空間 ($L = \prod d_i$) よりも大きくなります。本ライブラリは アイデンティティ・パディング・コンベンション を採用しています。すなわち、すべてのゲートは、論理部分空間の直交補空間に対して恒等写像として作用します。形式的には、論理ゲート $G$ に対して、エンコードされたユニタリは $U = E G E^\dagger + (I - \Pi)$ となります。ここで、$E$ はエンコーディング・アイソメトリ（等長写像）であり、$\Pi$ は論理部分空間への射影です。この選択により、無雑音シミュレーション中のリークを防ぎ、置換特性を保持します。

2. ゲート階層と分解:
   本ライブラリは、3 つの基本クラスを定義します：QuditGate（単一クディット）、QuditControlledGate（マルチコントロール）、および QuditMultiGate。

• 2 の累乗パス: すべてのオペランドの次元が 2 の累乗である場合、ゲートは標準的な Qiskit 量子ビットゲート（例：mcx, cp, synth_qft_full）の固定されたカスケードへと分解されます。これにより、Qiskit トランスパイラによる最適化が可能になります。
• 高密度合成フォールバック: 2 の累乗ではない次元の場合、ゲートは正確な高密度ユニタリ合成へとフォールバックします。本ライブラリは、完全な $D \times D$ ユニタリ行列を生成し（NumPy の __array__ プロトコル経由）、それを単一の不透明なユニタリ命令として埋め込み、その後 UnitarySynthesis パスによって分解します。

3. 回路モデルと測定:
   QuditQuantumCircuit は、2 つの同期されたビューを維持します：クディットレベルの命令ログ（.data）と、エンコードされた量子ビットレベルの回路（.circuit）です。測定は QuditMeasure を通じて行われ、これはエンコーディング量子ビットを古典ビット（ClByte）にマッピングします。専用のデコーディング層が、バックエンドのカウントをクディットレベルに変換し、「リーク」（結果 $\ge d$）を明示的に検出し、これらは設定可能なポリシー（keep, drop, raise）によって処理されます。

主な貢献

• トランスパイラ可能なクディット層の実装: 本研究は、先行研究 [1] の理論的なゲート定義を、動作可能でトランスパイラ可能な Qiskit 拡張へと変換しました。これは次元 $2 \le d \le 16$ および混合次元レジスタをサポートしています。
• 2 の累乗ではない次元に対する正確なシミュレーション: $d=2^m$ に限定されていた従来の量子ビット・アプローチとは異なり、本ライブラリは高密度ユニタリ・フォールバックを使用して、サポート範囲内の任意の $d$ に対して正確な合成をサポートします。
• 理論と実装の検証:
  • チェック A: 高密度ユニタリ行列が、2 の累乗および 2 の累乗ではない両方の次元において、論理部分空間上の離散フーリエ変換（DFT）を再現することを検証しました。これは [1] の検証を拡張するものです。
  • チェック B: 分解された回路（カスケードおよび高密度パスの両方）が、138 個のゲート・インスタンスにわたって、高密度ユニタリ定義と機械精度（$\approx 10^{-15}$）以内で一致することを確認しました。
  • パディング不変性: アイデンティティ・パディング不変性が正確に保持されることを確認し、有効な部分空間と無効な部分空間の間の混合を防いでいます。
• リーク検出: シミュレーションにおけるリーク事象を検出し、扱うための明示的なメカニズムを提供し、理想的な挙動とノイズ誘発エラーを区別します。

結果

数値検証により、実装の正当性が確認されました：

• QFT vs. DFT: QFT 実装は、次元 $d \in \{2, 3, 4, 5, 8\}$ において、理論的な DFT 行列と $10^{-15}$ 以内で一致しています。
• ゲート分解: SUMX, SUMP, SWAP, QFT を含むすべてのゲートは、高密度ユニタリ定義と一致する分解を生成します。
• リソース分析:
  • 2 の累乗: 分解は効率的で、トランスパイラに適しています。例えば、3 クトリット QFT ($d=3$) は 1,783 個の CX ゲートを必要とするのに対し、3 ククオートット QFT ($d=4$) はわずか 56 個の CX ゲートを必要とします。
  • 2 の累乗ではない場合: 高密度フォールバックは大幅なオーバーヘッドを伴います。3 クトリット QFT ($d=3$) は 1,783 個の CX ゲートに合成されますが、3 ククイント QFT ($d=5$) は 119,383 個の CX ゲートを必要とします。これは、大きな非 2 の累乗次元における高密度合成のコストを浮き彫りにしています。
  • メモリ制約: 高密度パスはメモリによって制限されます。16 レベルの 8 個のクディット（$N=32$）に対するゲートは、$\approx 256$ EiB のメモリを必要とするため、実用的な限界は一般的なハードウェア上で $N \approx 13\text{--}14$ となります。

意義と限界

本論文は、qiskit-qudits を、ユーザーがクディット・アルゴリズムを論理レベルで記述しながら、Qiskit のコンパイル・スタック、バックエンド、およびプリミティブへのフルアクセスを維持できるようにするための架け橋として位置付けています。これは、手動の実装においてエラーが発生しやすい「物理的に正しくない部分」と測定デコーディングの問題を明示的に解決しています。

限界とトレードオフ:

• 次元の上限: 本ライブラリは $2 \le d \le 16$ およびゲートあたりの最大 8 個のクディットを強制します。これらは公称値であり、実用的な限界は高密度ユニタリ合成に必要なメモリ（$N \lesssim 13\text{--}14$）によって決まります。
• リークのセマンティクス: アイデンティティ・パディング・コンベンションは、物理的なリークをモデル化するのではなく、リークした振幅を凍結します。これはテンソル積と可換ではないため、回路の等価性は投影された部分空間（$\Pi U \Pi$）上でチェックする必要があります。
• 最適化: 現在の実装は、トランスパイラにおける「Don't-care」行/列を活用していません。また、非 2 の累乗次元に対する高密度フォールバックは、ネイティブなクディット・シミュレータや最適化された可逆論理合成と比較して計算コストが高くなります。

本研究は、高密度フォールバックが非 2 の累乗次元に対する正確なシミュレーションを可能にする一方で、将来の改善は、次元に依存しない可逆論理合成、およびリソース・オーバーヘッドを削減するための、非物理的な部分空間のトランスパイラへの露出に焦点を当てるべきであると結論付けています。