量子計算のミラー回路、元へ戻る結果で確かめられること

量子コンピューターの計算結果を確かめるには、比較する答えが必要です。小さな回路なら、通常のコンピューターで同じ計算を再現できます。けれども、量子ビットが増えてその再現が重くなると、試したい回路だけでなく、答え合わせの計算にも大きな負担がかかります。

ミラー回路は、計算を進めたあとに、その操作を逆順に戻す回路を付ける方法です。理想的には出発点へ戻るので、途中の複雑な状態をすべて計算しなくても、最後に期待する答えを用意できます。QED-Cの共同研究は、この方法をハミルトニアン・シミュレーションの評価へ組み込んでいます。ただし、「戻れた」という結果と、途中で目的どおりの物理計算ができたことは、同じ意味ではありません。^[1]

物理モデルを回路へ移すときにも、誤差が入る

ハミルトニアンは、量子系のエネルギーや、状態が時間とともにどう変わるかを表す演算子です。ここでは時間によらないハミルトニアン H と、時間 t を指定し、状態の変化を表す U(t)=exp(−iHt) を量子回路で実現することを考えます。式ではプランク定数の一種である ℏ を1とする単位を使っています。^[1]

装置が一度に実行しやすい操作へ分ける方法の一つが、トロッター分解です。例えば H=A+B なら、Aによる短い変化とBによる短い変化を順に行い、それを何回も繰り返します。AとBを入れ替えると結果が変わる場合、この順番のある操作は、両方が同時に働く時間発展と一般には一致しません。これが近似に由来するトロッター誤差です。

時間を細かく分ければ、この近似誤差を抑えられます。一方、操作数は増えるため、実機では雑音の影響を受ける機会も増えます。回路が指定どおりに動いたかと、その回路が元の時間発展を十分よく表すかを、分けて調べる必要があります。^[1]

三つの答え合わせは、違う問いに答える

QED-Cの研究では、同じハミルトニアン・シミュレーションに対して、三種類の評価方法を整理しています。数値を見る前に「何を正しい分布として比べたか」を確かめると、結果の意味がはっきりします。ここでいう分布は、決めた測定方法で0と1の並びが出る確率です。量子状態のすべての情報を直接比較しているわけではありません。原論文で使う正規化ヘリンガー忠実度（normalized Hellinger fidelity）は、二つの測定結果の確率分布の近さを表す指標です。量子状態や量子操作そのものを比べる忠実度と、同じ数値になるとは限りません。^[1]

Method 1では、比較する側も同じトロッター近似を使います。実機と雑音なしの回路がよく一致しても、両方が同じ近似誤差を含んでいれば、その誤差は見えてきません。Method 2なら、回路に分割する前の時間発展と比べられます。ただし、ここでいう厳密な時間発展にも、古典計算の数値精度の限界があります。

Method 2を雑音なしの回路に対して行えば、近似そのものの影響を調べられます。まず小さい問題でこの差を確認し、次に実機での差を見ると、近似の刻み方と装置の雑音を考え分ける手掛かりになります。二つの指標の差をそのまま引いて「雑音だけの誤差」とする、という意味ではありません。

逆に戻すのは、作った回路そのもの

ある回路を C とすると、その操作を打ち消す逆回路は C⁻¹ です。先にC、続いてC⁻¹を実行すれば、理想的には C⁻¹C=I、つまり何もしない操作と同じになります。例えば0101から始めるなら、最後に0101が出るはずだと分かります。初期状態が複数の測定結果を取りうる場合は、その既知の分布と照合します。^[1]

大切なのは、逆に戻している相手が作成した近似回路Cである点です。Cが目的の時間発展Uと少し違っていても、Cの逆回路を正しく実行すれば戻れます。元へ戻れたことだけで、CとUの違いを判定することはできません。

これを一つの量子ビットで確かめてみます。説明用に H=X+Z、t=1、初期状態を0とします。XとZは異なる方向の状態変化を表すパウリ演算子で、順序を入れ替えると結果が変わります。短いXの操作、次に短いZの操作を行う組をK回繰り返す、一次のトロッター回路を考えます。

片道の結果は、Kを増やすとこの例の厳密な値へ近づいています。それでも、近似の粗いK=1を含め、逆回路の後はどれも元へ戻ります。「往復が合う」と「途中の答えが目標に合う」を分ける理由が、この二列の違いに表れています。実機でKを増やした場合に、雑音も含めた結果が同じように改善するとは限りません。

行きと帰りのずれが、打ち消し合うこともある

単純なミラー回路には、別の見落としもあります。例えば、ある軸の周りに90度回す指定が、装置の系統的なずれで100度になったとします。帰りの−90度という指定も−100度になれば、合計は0度です。往復だけを見ると元へ戻っても、片道では指定からずれています。これは、ずれが角度の符号に応じて同じ比率で生じると仮定した説明例です。

このように、量子操作の振幅や位相に系統的に生じるコヒーレント誤差は、往復の間で打ち消し合ったり、逆に強め合ったりすることがあります。原論文は、回路の中央へランダムなパウリ操作を入れ、後半もそれに合わせて変更する方法を扱っています。この乱択変種では、元のビット列がそのまま戻るとは限りません。初期状態に、追跡したパウリ操作を施した結果が、あらかじめ分かる比較対象になります。最後の答えを古典計算で追えるまま、単純な往復の対称性を崩して、誤差の見落としを減らす狙いです。ランダム操作を挟むだけで、後半をそのままにする方法ではありません。^[1]

固定した公開コードにも、通常の逆回路を使う分岐と、パウリ操作を考慮した回路を使う分岐があります。同じ「Method 3」という名前でも、どの変種を使い、何種類のランダム回路を何回ずつ測ったかを確認する必要があります。ここでは公開コードを読んだ範囲を扱い、QED-Cのプログラムや量子実機を実行した評価はしていません。^[2]

答え合わせが軽くなっても、試験の条件は残る

ミラー回路では、行きの回路に帰りの回路を加えるため、深さがおよそ2倍になります。原論文は往復で得た指標の平方根を使って片道と比較する整理も行います。その説明では、系全体にかかる脱分極雑音（global depolarizing noise）を仮定し、各ゲートとその逆は同じ程度の誤差を持つとしています。各ゲートの指標をFと置き、n個の操作の指標を Fⁿ と近似する考え方です。この条件で、往復の指標を片道の指標の二乗として扱い、平方根を比較に使います。一般の装置で、単に平方根を取れば片道の品質が正確に分かるわけではありません。状態の準備と測定に由来する誤差も、回路本体の誤差と分けて点検する課題です。^[1]

評価結果を選ぶ場面では、まず、物理モデルの近似精度を知りたいのか、用意した回路の実行品質を知りたいのかを決めます。前者なら、小さい問題で厳密な時間発展と比べる検証が手掛かりになります。後者で、途中の状態を古典計算で再現するのが難しければ、既知の答えを作れるミラー回路が候補になります。

そのうえで、初期状態、トロッター分解の順序と回数、測定方法、回路の変換条件、測定回数、往復の変種を揃えます。古典計算による答え合わせを軽くできることは、量子回路を任意の規模で正確に動かせることでも、実用計算で優位になったことでもありません。ミラー回路の結果は、何を確かめるために往復したのかと結び付けて読むと役立ちます。

参考資料

1. Avimita Chatterjeeほか（QED-C collaboration）「A Comprehensive Cross-Model Framework for Benchmarking the Performance of Quantum Hamiltonian Simulations」。arXiv:2409.06919v1、2024年9月11日投稿、PDF表記は9月12日。第II.B節、第III.A–D節、Fig.3、Table I・II、付録Aを参照。参照版PDF。本稿はこの投稿版に基づき、後の出版版の実験結果を検証した記事ではありません。
2. QED-C「QC-App-Oriented-Benchmarks」。commit 21dfc331ce61a2170b06f88259aa19f901c021bc、2026年10月18日確認。Hamiltonian Simulationの説明／回路生成と逆回路の分岐。公開実装は性能標準ではなく、参照用アプリケーションの試作群として提供されています。