量子コンピューターで材料のエネルギーを計算し、正解に近い値が出たら、電子の状態も正しく再現できたと考えてよいのでしょうか。グラフェンを扱う新しい研究は、近いエネルギーを持つ状態が並ぶとき、その確認だけでは不十分な場合を示しています。
2026年11月3日公開のプレプリントでは、大学の研究チームがグラフェン単層と細いナノリボンのバンド構造を計算し、一部をIBM Quantumの実機で測定しました。対象は電子間相互作用を含めないモデルです。実機での測定、雑音を含むシミュレーション、波動関数の評価を分けて読むと、今回の成果と残る課題が見えてきます。［1］
簡単にまとめると？
- 実機で測ったのは、シミュレーションで先に最適化した回路です。実機の測定を繰り返して最適化全体を回した結果とは区別します。
- 回路を深くすると理想計算は改善しましたが、雑音モデルでは測定誤差が増えました。表現力とゲート誤差の両方が深さを制限します。
- エネルギーの近い状態が混ざると、エネルギー誤差だけは小さく見える場合があります。狙った状態の空間にどれだけ入っているかという別の診断が役立ちます。
今回計算したグラフェンの範囲
グラフェンは炭素原子が蜂の巣状に並ぶ材料です。研究では、電子が原子間を移る様子を表すタイトバインディングモデルを使い、単層と、端の形がジグザグ型・アームチェア型のナノリボンを調べています。電子同士の強い相関を直接解いた計算ではなく、その研究へ向かうための、正解と比較しやすい出発点です。［1］
バンド構造は、電子の波数ごとに取り得るエネルギーを並べたものです。単層では各波数の問題が2×2の行列にまとまり、1量子ビットで表せます。リボンでは単位胞の原子数に応じて行列が大きくなります。ここで少数の量子ビットへ符号化しているのは、この一電子の行列であり、多数の電子が相互作用する全状態を同じ規模で計算したわけではありません。［1］
最も低いエネルギーの状態を探すVQE（変分量子固有値ソルバー）に加え、より高い状態にはVQD（変分量子デフレーション）を使います。VQDは、すでに見つけた状態と重なるとコストが増える項を加え、別の状態を順に探す方法です。そのため、先に求めた状態の誤差が後の探索へ伝わる可能性もあります。［1］
IBM実機は、最適化した回路の測定に使った
実機での最適化は、パラメーターを変えるたびに測定を繰り返す必要があります。今回の研究では、その負担を抑えるためにパラメーターの最適化をシミュレーション内で行い、得られた回路だけを実機へ送る方法を採りました。論文中のすべての実機結果が、この手順によるものです。［1］
単層は1量子ビットで、Heron系のibm_marrakeshを使用。リボンは最も細い2量子ビットの構成で、アームチェア型をibm_fez、ジグザグ型をibm_kingstonで測定しています。より幅広いリボンや、後述する3量子ビットの波動関数評価は、同じ実機実験の範囲ではありません。［1］
| 対象 | 量子ビット数 | 平均絶対誤差／|t| |
|---|---|---|
| グラフェン単層 | 1 | 約0.021 |
| アームチェア型リボン | 2 | 約0.044 |
| ジグザグ型リボン | 2 | 約0.065 |
tはモデル内で電子が隣の原子へ移る強さを表すホッピングのパラメーターです。上の値はtの大きさを尺度にした誤差で、単純な「正答率97.9%」などへ読み替えるものではありません。また、実機の結果には読出し誤差の緩和を適用しています。雑音下で最適化したシミュレーションの総誤差と比べるときは、最適化の誤差まで含むか、最後の測定だけかをそろえる必要があります。［1］
深い回路ほど正確、とは限らない
回路に調整可能な操作を増やすと、表せる状態の幅は広がります。4量子ビットのジグザグ型リボンで回路の繰返しを1回から4回へ増やした比較では、理想シミュレーションの平均誤差は約0.11から0.000026へ減りました。一方、装置の雑音モデルを使うと、測定したエネルギーの誤差は約0.10から0.162へ増えています。いずれも|t|を単位とする値です。［1］
このとき、2量子ビットゲート数は3個から12個へ増えました。雑音下で見つけたパラメーターを理想的に評価すると改善も見られるため、単に最適化がすべて失敗したという話ではありません。表現力を高めるために追加したゲートが、測定時には誤差を積み重ねるという関係です。この深さの比較は雑音シミュレーションであり、4量子ビットの実機で同じ数値を実測した結果ではありません。［1］
エネルギーが近くても、状態が違うことがある
別の落とし穴は、エネルギー準位が密に並ぶときです。狙った状態の代わりに、その近くの状態が混ざっても、エネルギーの平均は大きく変わらない場合があります。研究では、幅を広げたジグザグ型リボンでエネルギー誤差が小さくなる場面を、波動関数の再現が改善したとは扱っていません。準位間隔も狭くなるため、その間隔に対してどれくらい誤っているかを見る必要があります。［1］
そこで使うのが部分空間への包含度（subspace containment）です。計算で得た状態のうち、正解の状態が張る空間にどれくらいの重みが入っているかを測ります。同じエネルギーを持つ複数の状態がある「縮退」では、一つのベクトルだけでなく、それら全体が張る空間への重なりを見ることで、基底の選び方に左右されない評価ができます。［1］
今回の包含度は、雑音下で最適化したパラメーターを、雑音なしの状態ベクトルで評価しています。実機の出力状態をすべて測り直したトモグラフィーではありません。3量子ビットのリボンでは基底状態は高い包含度を保ちましたが、アームチェア型の第2励起状態で約0.809まで下がる点もありました。エネルギーだけでは見落としやすい、近隣の状態への混入を調べる役割があります。［1］
次の材料計算を比べるための基準
この研究は、古典計算で厳密に解けるモデルを使って量子回路を検査したもので、古典計算に対する優位性の実証ではありません。波動関数の診断にも、厳密対角化で得た参照状態を用いています。大きな未知の問題へ進んだ場合に、同じ診断を同じ費用で行えるとまでは言えません。［1］
材料計算の報告を比べる際には、相互作用を含むモデルか、最適化をどこで実行したか、雑音と測定回数は何か、エネルギー以外の状態の確認があるかを並べると、評価の範囲が明確になります。少ない量子ビットでバンドを描くことから、目的の電子状態を信頼して準備することへ進むには、この区別が役立ちます。
実験に使われた計算基盤の背景は、CyclewaveのIBMの企業解説でも紹介しています。
出典
- Medrano Millán et al., Band structure of two-dimensional graphene and of one-dimensional nanoribbons via variational quantum algorithms: from energies to wavefunction quality, arXiv:2611.04975v1（2026年11月3日提出、PDF表紙日付11月4日）。本文III〜VII節、表I、図5・9・11・12。