分子の電子の状態を、少ない量子ビットで調べたい。そのためには、同じ問題を小さく詰め込むだけでなく、何を計算するかを組み替える方法があります。電子一つひとつの複雑な組み合わせではなく、電子の分布を手がかりにする量子密度汎関数理論（Q-DFT）も、その一つです。
2026年11月3日公開のプレプリントは、水素化リチウム（LiH）を対象に、Q-DFTで用いる計算を二つの光量子ビットで実行しました。ただし、実機で確かめたのは、あらかじめ用意した問題の最適化です。電子密度から問題を作り直す循環全体と、分子の結合を正確に描くことには課題が残りました。［1］
簡単にまとめると？
- 四つの軌道を基底として使うコーン・シャム問題を、2量子ビットへ符号化しました。電子間相互作用をすべて同じ精度で扱う多体計算とは条件が違います。
- 光回路上で最適化は進みましたが、密度を更新する自己無撞着計算は実機では閉じていません。各段階の問題は、雑音のない古典シミュレーションから供給しました。
- 測定数の不足と、光子を完全には区別不能にできないことが誤差に効きます。量子ビット数だけでなく、測定費用・結合の再現・循環の収束を分けて見る研究です。
四つの軌道を、二つの量子ビットに載せる
密度汎関数理論（DFT）は、電子の密度を使って物質の性質を求める方法です。よく使われるコーン・シャム法では、相互作用する多数の電子の問題を、密度に依存する有効な場の中で動く電子の問題に置き換えます。電子間の効果は有効な場に含めますが、実際には交換・相関を表す汎関数の近似が必要です。Q-DFTは、この枠組みの一部を量子計算で解きます。［1］
今回使う基底は、リチウムの1s・2s・2pzと、水素の1sに対応する四つの直交化された軌道です。これらを2量子ビットの四つの状態に対応させます。四つの軌道を一度に四つの別々の量子ビットへ割り当てるのではなく、各軌道の重ね合わせを2量子ビットの状態として表すため、少ない量子ビットで扱えます。
この縮約には前提があります。対称性から寄与しないリチウムの2px・2pyを除き、さらにコーン・シャムの一体問題へ移しています。論文が比較する多体波動関数の表現は、最小基底全体なら12量子ビット、同じ軌道の除外後でもJordan–Wigner変換では8量子ビットです。12から2という数字だけで、同じ多体計算をそのまま6分の1にしたとは読めません。［1］
二つの占有軌道を、一緒に最適化する
LiHには四つの電子があり、この計算では二つのコーン・シャム軌道を二個ずつ占有します。そのため、一番低い状態だけを探せば十分ではありません。研究では、二つの直交する初期状態へ同じ可変回路をかけ、二つのエネルギーの重み付き和を小さくするアンサンブルVQEを使いました。
同じ重みで最適化する場合、重要なのは二つの状態が張る空間です。その中で状態の組み合わせが変わっても、和と電子密度を正しく求められます。各状態がそれぞれ厳密な固有状態に一致することを、最適化の必須条件にはしていません。ここでいう第一励起状態も、有効な一体問題の状態であり、LiH分子全体の励起をそのまま測ったものではありません。［1］
光回路では、一つの光子が二つの経路のどちらを通るかで1量子ビットを表します。装置は8モードの集積回路で、実際の計算には4計算モードと2補助モードを使いました。二つの光子を入力し、回転とC-NOTに対応する操作を行い、出力を測って古典側で四つのパラメーターを更新します。
ただし、このC-NOTは所定の検出結果を選び出す事後選択を必要とします。二つの量子ビットを表す各経路の組で、光子が一個ずつ検出された事象を採用します。入力した光子がすべて有効な計算結果になるわけではなく、将来の規模拡大では成功率と損失も数える必要があります。［1］
実機で閉じていない、密度の更新ループ
コーン・シャム法では、仮の密度から有効な問題を作り、それを解いて新しい密度を得ます。この更新を、密度とエネルギーが十分変わらなくなるまで繰り返します。これが自己無撞着計算、またはSCFループです。
雑音のない状態ベクトルシミュレーションでは、この循環が6回の更新で収束しました。ところが、測定回数を有限にしたシミュレーションでは、10億ショットでも、同じ6回の範囲で完全には収束しません。後半の問題同士のエネルギー差が小さく、測定の揺らぎの中で見分けにくくなるためです。［1］
そこで実験では、雑音のない古典シミュレーションから得た六つのハミルトニアン、つまり解くべき有効な問題を、個別に光回路へ渡しました。実機の測定で得た密度を次の問題へ戻して、循環全体を収束させたわけではありません。
| 検証の場 | 確かめた範囲 | 残った条件 |
|---|---|---|
| 雑音なしの古典シミュレーション | 6回のSCF更新で収束し、参照計算へ近づく | 交換・相関汎関数などの近似は残る |
| 有限測定のシミュレーション | 測定の揺らぎで密度更新が停滞する影響 | 10億ショットでも6回以内の完全収束なし |
| 光量子回路の実験 | 事前に用意した各問題のVQE最適化と電子密度の推定 | 閉じたSCF循環と定量的に正確な結合描写は未達 |
光子を増やすだけでは、消えない誤差
測定回数が少ないと、光子を数える統計の揺らぎが目立ちます。一方、十分な数を集めても、二つの光子の区別不能性が完全でなければ誤差が残ります。干渉させたい光子に、見分けられる違いが残ると、意図した光回路の働きからずれるためです。
実験のHong–Ou–Mandel干渉の可視度は約89％でした。これは光子の区別不能性を評価する指標で、計算全体の正答率ではありません。雑音モデルでは、測定数を増やすと誤差の支配要因が統計から区別不能性へ移ります。さらに、実際の装置にはモデルに含めていない位相制御などの誤差もあります。［1］
電子密度を見ると、リチウム側から水素側への電荷の移動は実験でも表れました。しかし、原子間の結合に関わる密度の形は参照結果とずれ、イオン結合を定量的に正確に描くには足りませんでした。「最適化の値が下がった」と「化学的に必要な答えが得られた」を分ける意味が、ここにあります。
この研究から読み取れるのは、少ない量子ビットへ問題を写す工夫と、それを測って更新する費用は別だということです。次の比較では、符号化する量子ビット数に加え、採用できる検出事象の割合、密度の精度、閉じた循環の収束、全体の計算時間をそろえる必要があります。大きな分子での量子優位を示した段階ではありません。シミュレーションには公開ソフトウェアが使われていますが、この研究の結果を支えるコードとデータは請求に応じて提供すると記されています。公開実装があることと、結果一式が公開され再現確認済みであることも分けて捉えたいところです。
出典
- Giovanni Adrianoほか「Quantum Simulation of Electronic Structure using Quantum Density Functional Theory on a Programmable Integrated Photonic Circuit」、arXiv:2611.04694v1、2026年11月3日投稿。論文・書誌情報。