IBMの104量子ビットで粒子衝突をシミュレーション、重い励起の生成を示す証拠

何を実行したのか――チップ内で粒子をぶつけたわけではない

Roland C. Farrell、Nikita A. Zemlevskiy、Marc Illa、John Preskillの研究チームは、IBMの104量子ビット量子プロセッサーibm_marrakeshを用い、1次元イジング場理論での衝突をデジタルかつリアルタイムにシミュレーションしました。

準備したのは、離れた位置から互いに向かって進む二つの軽い粒子の波束です。系を時間発展させ、波束が重なって衝突し、再び離れた後の局所エネルギー密度を測定しました。衝突後のダイナミクスに到達するまで、最大5,589個の2量子ビットゲートを用いています。7

重要なのは、これは半導体チップの中で現実の素粒子を衝突させた実験ではない、という点です。量子ビットに格子化した場の理論を符号化し、その量子力学的な時間発展を模擬したものです。実時間の散乱では多体系の量子ダイナミクスを扱う必要があり、一般には古典計算で難しくなるため、量子ハードウェアで直接表現する意義があります。

シミュレーションの流れ

計算は大きく3段階です。

1. 入射粒子を準備する：局在した二つの軽い粒子波束を、反対向きの運動量を持つように遠く離して初期化する。
2. 衝突をまたいで時間発展させる：トロッター分解した量子回路で時間発展を近似し、波束の接近、重なり、分離を追う。
3. 衝突後を測定する：局所エネルギー密度を測り、エネルギーがどのように分布したかを調べる。7

中心的な観測結果は、非弾性の粒子生成を示す証拠です。二つの軽い粒子がそのまま二つの軽い粒子として飛び去るだけでなく、最終状態が「軽い粒子1個と、より遅く重い励起1個」から成ることと整合的だったと報告されています。7

最大の難所は、衝突前の波束を正しく作ること

散乱計算では、入射状態の質が結果を左右します。研究チームはまず、波束の空間的な形、運動量成分、量子数を定めるよう設計した (|W(k_0)\rangle) 型の状態を作りました。ただし、この出発状態には不要な多粒子成分が混じり得ます。そこで次に、対称性を保つ変分回路でエネルギーを下げ、狙った一粒子固有状態へ近づけました。7

大規模格子向けのイジング波束回路は、104量子ビットの実機に投入する前に、行列積状態（MPS）による古典シミュレーションで設計されています。7

W状態、途中測定、フィードフォワードが果たした役割

通常、波束が多くの格子点に広がるほど、その準備回路は深くなりがちです。研究チームはW状態準備の技法を応用し、途中測定（mid-circuit measurement）と古典フィードフォワードを活用しました。

• 途中測定は、量子プログラムが終わる前に量子ビットを測定し、古典的な測定結果を取り出す操作です。
• フィードフォワードは、その測定結果に応じて後続の操作を条件付きで選ぶ仕組みです。
• この組み合わせにより、補助量子ビットを使わず、物理的な相関長への依存を除けば、波束の大きさや空間次元に依存しない回路深さで波束を準備するプロトコルを可能にします。7

著者らはこのスケーリングを、古典前処理のコストが大きい、あるいは波束サイズに対して回路深さが多項式的に増える従来手法に比べて、超指数的な改善だと説明しています。27

「歪度」で重い粒子を見分けた

研究チームは、通常の粒子検出器のように飛跡を撮影して出射粒子を同定したのではありません。手掛かりにしたのは、測定されたエネルギー密度プロファイルの形です。

比較的弾性的な散乱なら、衝突後に出ていく二つの軽い粒子は、より対称的なエネルギー分布を作ると予想されます。一方、非弾性の経路では、重い粒子は遅く動くため、衝突後のエネルギー密度に検出可能な非対称性が現れます。7

その非対称性を数値化したのが**歪度（skewness）**です。歪度は、分布のどちら側の裾がより長い、または重いかを表す統計量です。予想される重粒子の軌跡に対応する非ゼロの歪度が、非弾性過程のシグネチャーとなりました。7

これは、単にエネルギーの配置が変化したことを意味するだけではありません。入射した軽粒子の運動エネルギーが、軽い粒子と新しい重い励起を含む最終状態へ変換されたことと整合する結果です。7

なぜQCDやハドロン生成への道筋として注目されるのか

高エネルギー衝突では、粒子が新たに生み出されることが本質的な現象です。長期的には、研究者たちは強結合の量子場理論で、衝突がどのように粒子や集団的な振る舞いを生むのかを予測したいと考えています。これはハドロン生成、さらに究極的には量子色力学（QCD）につながる課題です。

一般的な実時間量子ダイナミクスは古典的手法では特に扱いにくく、将来の量子コンピュータが力を発揮し得る領域とみなされています。7

今回の波束準備法は、イジング理論の実証だけに限られません。著者らは、1次元スカラー場理論、シュウィンガー模型、2次元イジング場理論についても、この手法を構成または実演しています。27

ただし、QCDの衝突を解いたわけではない

成果の範囲は明確に区別する必要があります。

• 対象は1次元イジング場理論であり、QCD、現実の3次元ハドロン衝突、あるいは加速器実験のジェットをシミュレーションしたものではありません。7
• 得られたのは非弾性経路の証拠であり、誤り耐性量子計算による系統的に高精度な散乱予測ではありません。7
• 実験はノイズのある量子ハードウェアと誤り緩和技術に依存しており、完全な量子誤り訂正を備えた計算とは根本的に異なります。78
• より高次元のゲージ理論へ拡張するには、符号化、長時間の実時間発展、スケーラブルな状態準備、多粒子観測量の精密な抽出などで、はるかに高性能な量子系が必要です。

当面の成果は限定的ながら具体的です。量子プロセッサー上で準備した二粒子状態を衝突後まで追跡し、より重い出射励起の存在と整合する測定シグナルを得たことです。これにより、将来の量子ハードウェアでスケーラブルな散乱状態準備と実時間の粒子生成研究を進められる見通しが、より現実的になりました。7