Rovibrational energy levels of H$_2$O by quantum computing

この論文は、イオントラップ型量子コンピュータを用いて水分子のワトソンハミルトニアンの量子ビット表現を導出し、量子選択構成相互作用法と古典計算を組み合わせるハイブリッド手法により、数 cm$^{-1}$の精度で低励起回転振動エネルギー準位を計算することを報告しています。

要約

この論文は、**「量子コンピュータを使って、水分子（H2O）がどう振動し、どう回転しているかを、これまでになく正確に計算しようとした」**という研究報告です。

専門用語を並べると難しく聞こえますが、実はとても面白い「料理」と「迷路」の話に例えることができます。

1. 物語の舞台：水分子の「ダンス」

まず、水分子（H2O）を想像してください。これは 3 つの原子（酸素 1 つ、水素 2 つ）でできています。
この分子は止まっているのではなく、常に**「ダンス」**をしています。

• 振動（Vibration）： 原子同士がバネでつながれているように、伸び縮みしたり曲がったりする動き。
• 回転（Rotation）： 分子全体がくるくる回る動き。

この「振動」と「回転」が組み合わさった状態を**「回転振動状態（ロビブレーション）」**と呼びます。この状態のエネルギー（ダンスの激しさ）を正確に知ることは、大気中の水蒸気が太陽光をどう吸収するか、あるいは宇宙の星の温度を測るなど、科学にとって非常に重要です。

2. 従来の「古典コンピュータ」の限界

これまで、このダンスのエネルギーを計算するには、スーパーコンピュータ（古典コンピュータ）を使っていました。
しかし、分子が大きくなったり、回転が複雑になったりすると、計算すべきパターンが**「天文学的な数」**に膨れ上がります。

• 例え話： 1 人のダンサーの動きを追うのは簡単ですが、100 人のダンサーが複雑に絡み合うダンスを、1 人ずつ順番にシミュレーションしようとしたら、計算が終わる前に宇宙が滅んでしまうかもしれません。
• 水分子は小さいですが、それでも「回転」と「振動」を同時に高精度で計算するのは、古典コンピュータには非常にハードルが高いのです。

3. 新しいアプローチ：量子コンピュータという「魔法の鏡」

そこで、この研究では**「量子コンピュータ」**という新しい道具を使いました。
量子コンピュータは、計算のルールそのものが分子の動きと似ているため、この問題を解くのに適しています。

しかし、現在の量子コンピュータは**「ノイズ（雑音）」**が多く、完全には信頼できません。完全な計算をするには、まだ技術が追いついていないのです。

4. 解決策：「QSCI」という賢い選び方

研究チームは、完全な計算をしようとして失敗するのではなく、**「QSCI（量子選択コンフィギュレーション相互作用）」**という賢い戦略を取りました。

これを**「料理の味見」**に例えてみましょう。

• 従来の方法（VQE など）：
  巨大な鍋（全エネルギー）を、すべて一口ずつ味わって、一番美味しいものを探す。
  → 鍋が巨大すぎると、味見するだけで疲れてしまいます（計算リソース不足）。
• この論文の方法（QSCI）：
  1. まず、量子コンピュータを使って「どんな味がしそうな料理か」を**「味見（サンプリング）」**します。
  2. 量子コンピュータは「ここが美味しそうだ」という**「重要な味（重要な状態）」**だけを抽出して、古典コンピュータに渡します。
  3. 古典コンピュータは、その「重要な味」だけを組み合わせて、最終的な「完璧なレシピ（エネルギー値）」を計算します。

つまり、**「全部を計算するのではなく、量子コンピュータが『ここが重要だ！』と選りすぐった部分だけを、人間（古典コンピュータ）が詳しく調べる」**という、チームワークの勝利です。

5. 実験の結果：水分子の「ダンス」を捉えた

研究チームは、日本の理化学研究所（RIKEN）にある**「Quantinuum Reimei（レイメイ）」**という量子コンピュータを使って実験を行いました。

• 結果： 水分子の回転振動エネルギーを、**「1 cm⁻¹（波数）」**という非常に高い精度で計算することに成功しました。
• 意味： これは、実験室で実際に測定された値と、ほぼ同じ精度が出たということです。
• 驚き： 量子コンピュータには「ノイズ（雑音）」がありましたが、それをうまく補正する手法（パリティのチェックなど）を使うことで、ノイズに負けない結果を出せました。

6. この研究の未来

この研究は、**「量子コンピュータが、化学の分野で実際に役立つ（量子優位性を示せる）最初のステップ」**になる可能性があります。

• なぜ重要か？
  これまで「電子の動き」を計算する研究は多かったですが、「分子全体の振動・回転」を量子コンピュータで計算したのは、この研究が先駆けの一つです。
• 将来：
  もしこの技術がさらに進歩すれば、**「人間には計算しきれない巨大な分子」や「新しい薬の候補」**の動きを、量子コンピュータを使って正確に予測できるようになります。

まとめ

この論文は、**「不完全な量子コンピュータでも、賢い選び方（QSCI）を使えば、水分子の複雑なダンスを、実験とほぼ同じ精度で再現できる！」**と証明した画期的な研究です。

まるで、**「騒がしい部屋（ノイズのある量子コンピュータ）で、重要な会話（エネルギー値）だけを聞き取り、それを整理して完璧な記録にする」**ような、非常にクリエイティブで実用的な成果と言えます。

技術概要

以下は、RIKEN-iTHEMS-Report-26「量子コンピューティングによる H2O の回転振動エネルギー準位の計算」の技術的サマリーです。

1. 研究の背景と課題 (Problem)

分子のエネルギー準位を計算する際、電子構造問題（電子の波動関数とエネルギーの計算）は量子コンピューティングの主要な応用分野として長年研究されてきましたが、**回転・振動構造（rovibrational structure）**の問題については、現在の量子コンピュータを用いた研究は限定的でした。

• 既存の課題:
  • 従来の量子化学計算（VQE など）は、電子構造問題に焦点が当てられており、分子の回転・振動運動を同時に扱う「回転振動シュレーディンガー方程式」の解法は十分に発展していませんでした。
  • 回転・振動波動関数は原子数に比例して変数が増加し（$N_{var} = 3N_a - 3$）、古典コンピュータでは多原子分子の高精度計算が困難です。
  • 既存の量子コンピュータを用いた振動構造の研究（O3, H2O, CO2 など）では、分子の回転運動が考慮されておらず、回転・振動結合項（Coriolis 結合など）を無視しているものが大半でした。
  • 現在の量子コンピュータはノイズ（NISQ 時代）にさらされており、大規模なハミルトニアンの直接対角化や VQE のような変分法は、項数が膨大になるため実行が困難です。

2. 手法とアプローチ (Methodology)

本研究では、現在のノイズあり量子コンピュータ（Quantinuum のトラップドイオン量子コンピュータ「Reimei」）を用いて、水分子（H2O）の回転振動エネルギー準位を計算する新しいハイブリッド手法を提案・実装しました。

• ハミルトニアンの導出:

  • Watson のハミルトニアンを量子ビット形式に変換しました。これには、剛体回転項、調和振動子項、非調和項、振動コリオリ結合項、そして回転・振動結合項が含まれます。
  • 水分子のポテンシャルエネルギー曲面は CCSD(T)/aug-cc-pVQZ レベルで計算し、4 次までのテイラー展開を用いて非調和結合定数を評価しました。
  • 量子ビットへの符号化には「バイナリマッピング」を採用し、振動量子数と回転量子数を量子ビット列にマッピングしました。

• 量子選択配置相互作用 (QSCI) 法の適用:

  • 従来の VQE 法はハミルトニアンの項数が膨大（水分子で 20 万項以上）なため実行不可能でした。
  • 代わりに、QSCI (Quantum-Selected Configuration Interaction) 法（またはサンプルベース量子対角化法）を採用しました。
  • 手順:
    1. 試行波動関数の生成: 量子コンピュータ上で、ハミルトニアンの時間発展（Suzuki-Trotter 分解）を用いて試行波動関数 $|\psi(\tau, N_{ST})\rangle$ を準備します。ここで、時間ステップ $\tau$ やステップ数 $N_{ST}$、カットオフパラメータ $\lambda$ を変数として最適化します。
    2. 基底状態のサンプリング: 量子コンピュータ上でこの波動関数を測定し、計算基底における確率分布 $p(x)$ を取得します。
    3. 基底集合の選択: 確率が高い基底状態を選択し、コンパクトな基底集合 $\Omega$ を構築します。
    4. 古典的な対角化: 選択された基底集合 $\Omega$ 上でハミルトニアン行列を構成し、古典コンピュータで対角化してエネルギー固有値を求めます。
  • 励起状態の計算には、異なる初期状態から異なる基底集合を生成し、それらを直交化して最終的な行列を対角化するという戦略を用いました。

• 誤差対策:

  • 量子コンピュータのノイズにより、パリティ（対称性）が破れた状態が測定される可能性があります。本研究では、正しいパリティを持つ状態のみを事後選択（post-selection）し、分布を再規格化する簡易的な誤差軽減手法を適用しました。

3. 主要な貢献と結果 (Key Contributions & Results)

• 回転・振動結合の重要性の証明:

  • 水分子のエネルギー準位計算において、非調和項だけでなく、振動コリオリ結合項や回転・振動結合項を含めることが、高精度な結果を得るために不可欠であることを示しました。
  • 特に回転量子数 $J > 0$ の場合、回転・振動結合項を無視すると、エネルギー準位が最大 30 cm$^{-1}$ 程度シフトし、状態の順序さえ入れ替わることが確認されました。

• 高精度なエネルギー準位の計算:

  • Quantinuum の Reimei 量子コンピュータ（20 量子ビット、全結合、2 量子ビットゲート誤差 $\sim 10^{-3}$）を用いた実験により、低励起状態の回転振動エネルギー準位を数 cm$^{-1}$ の精度で再現することに成功しました。
  • 特定の条件下（$v_{max}=3, J=0$ など）では、実験値との誤差が1 cm$^{-1}$ 未満に収まるケースも確認されました。これは「分光学的精度（spectroscopic accuracy）」の基準を満たすレベルです。

• ノイズ耐性とサンプリングの有効性:

  • 量子ノイズが存在する環境下でも、サンプリングされた基底集合が有効な波動関数の主要部分を捉えており、古典的な摂動論に基づくサンプリングよりも優れた、あるいは同等の精度で基底集合を構築できることを示しました。
  • 場合によっては、適度なノイズが適切な基底関数の探索に寄与する可能性も示唆されました。

• スケーラビリティの示唆:

  • 基底集合のサイズ $|\Omega|$ を増やすことでエネルギーが収束することを確認し、この手法がより大きな分子やクラスターへ拡張可能であることを示唆しました。

4. 意義と将来展望 (Significance)

• 量子優越性の検証可能性:
  • 回転振動構造の計算は、高精密分光実験（HITRAN データベースなど）と比較可能な高精度な実験データが存在するため、量子コンピュータの計算結果を直接検証できます。これは、化学分野における「検証可能な量子優越性（verifiable quantum advantage）」の実現に向けた重要なステップです。
• NISQ 時代の有効活用:
  • 完全誤り耐性量子コンピュータが実現する前の、現在のノイズあり量子コンピュータ（NISQ）でも、適切なアルゴリズム（QSCI）とハイブリッド手法を用いることで、化学的に意味のある高精度計算が可能であることを実証しました。
• 今後の展開:
  • 本研究はベンゼンや水のクラスターなど、より大規模な分子系への適用、およびより効率的な量子ビット符号化や基底集合選択法の開発へと発展させることが期待されます。

結論:
本論文は、量子コンピュータを用いて分子の回転・振動運動を高精度にシミュレーションする最初の体系的な実証の一つであり、特に回転・振動結合を考慮した Watson ハミルトニアンの量子ビット化と、QSCI 法によるノイズ耐性のある計算手法の確立において画期的な成果を挙げています。