原子は物質の構成要素ですが、学校の図解でよく描かれているような、固定的で静止した球体ではありません。あらゆる原子の中心には原子核があり、その周囲を電子の雲が取り囲んでいます。ほとんどの化学計算において、科学者は重い原子核を固定された点として扱い、より軽い電子がその周囲を駆け巡るものとして扱います。これは多くの物質に対してうまく機能する簡略化です。しかし、最も軽い元素である水素の場合、この簡略化は通用しません。水素の原子核は非常に軽いため、小さな静止したビー玉のように振る舞うのではなく、量子物理学によってのみ記述できる方法で、エネルギー障壁の間を震えたり、トンネル効果によって通り抜けたりします。この挙動は「核量子効果」として知られており、光合成のような生物学的プロセスにおけるプロトンの移動や、水分子がどのように結合するかを理解する上で極めて重要です。しかし、こうした原子核の震えるような量子的な性質を従来のスーパーコンピュータでシミュレーションすることは、数学的な複雑さが手に負えないほど急速に増大するため、しばしば不可能です。
ドイツ航空宇宙センターとヘルムホルツ・インスティテュート・ウルムの研究チームは、シミュレーションを量子コンピュータへと移行させることで、この問題の解決に向けた重要な一歩を踏み出しました。彼らは、量子プロセッサの力と古典的コンピュータの信頼性を組み合わせた、新しいハイブリッド手法を開発しました。彼らのアプローチでは、量子コンピュータは、ある瞬間に原子に作用するエネルギーと力を決定する高度に特化した計算機として機能し、一方で古典的コンピュータは原子自体の動きを処理します。決定的なのは、この手法が特定の原子核を固定された点としてではなく、一度に複数の場所に存在し得る量子オブジェクトとして扱うことであり、これにより分子がどのように振動し反応するかについて、より現実的な描写が可能になります。
アルゴリズムをテストするために、研究者たちは水素分子、水分子、そして2つの水分子が1つのプロトンを共有しているズンデルイオンという3つの特定の系をシミュレートしました。ズンデルイオンは、化学や生物学における根本的なプロセスであるプロトン移動の最小モデルとなるため、特に重要です。チームは、原子核の量子的な性質を無視した場合、シミュレーションの結果がしばしば不正確になることを発見しました。しかし、これらの核量子効果を含めると、分子の挙動は劇的に変化しました。ズンデルイオンの場合、これらの効果を含めることで、プロトンが2つの水分子間を自由に移動することを妨げていたエネルギー障壁が取り除かれました。谷に閉じ込められる代わりに、プロトンの動きは実質的に障壁のないものとなり、容易に前後へと往来できるようになりました。この変化は実験データと非常によく一致しており、分子の予測振動数における誤差を大幅に減少させました。
研究者たちはまた、現在利用可能な（まだ初期段階でノイズの多い）量子コンピュータにおいて、これらのシミュレーションをいかに実用的なものにするかについても探求しました。彼らは、問題のエンコーディング（符号化）の方法をいくつかテストし、非常に正確だが複雑な手法と、現在のハードウェア向けに設計されたより単純でコンパクトな手法を比較しました。その結果、複雑な手法が最も精密な結果をもたらす一方で、簡略化されたコンパクトな手法も驚くほど効果的であることが分かりました。水分子や水素分子については、これらの合理化されたアプローチは、はるかに少ない計算リソースで、複雑な手法とほぼ同等の結果を生み出しました。これは、このアルゴリズムが、限られた能力しか持たない近未来の量子デバイスでも実行可能であることを示唆しています。
しかし、完璧な精度への道のりは決して平坦ではありませんでした。チームは、シミュレーションの成否が、原子核の量子状態をどのように表現するかという選択に大きく依存することを発見しました。水分子については単純な表現で十分でしたが、より複雑なズンデルイオンについては、正しい相互作用を捉えるために、より詳細な核軌道の記述を含める必要がありました。より詳細な記述を用いたとき、シミュレーションは実世界の実験で見られるものと一致する、特定の二峰性のパターンを再生することに成功しました。この追加の詳細がなければ、シミュレーションはこれらの重要な特徴を見落としていました。本研究は、現在の量子コンピュータは、まだ大規模な化学系のあらゆる細部をシミュレートできるほど強力ではないものの、このハイブリッドなアプローチが有望な道筋を提供していると結論付けています。どの部分を量子コンピュータでシミュレートし、どの部分を古典的に処理するかを慎重に選択することで、科学者たちは、以前は到達不可能であった精度で原子核の量子的な振る舞いをモデル化し始めることができるのです。
技術要約：量子コンピュータにおける核効果を伴う分子動力学
問題提起
プロトン・トンネル効果や零点エネルギーといった核量子効果（NQE）は、プロトン移動や水素結合などのプロセスを正確に記述するために極めて重要である。しかし、これらの効果を古典的ハードウェア上の量子化学計算に組み込むことは、フル構成相互作用（FCI）空間がシステムサイズに対して指数関数的に増大するため、計算量的に困難な場合が多い。量子コンピュータは、システムサイズに対する量子ビット数の要求が線形スケーリングするという利点があるものの、NQEを伴う動的なシミュレーションにおいて実用的な量子優位性を実証することは依然として課題となっている。既存のアルゴリズムは、核量子効果の欠如、あるいは動力学のシミュレーションの欠如、または既知の軌道への依存といった問題を抱えている。
手法
著者らは、核量子効果（NQE）を核電子軌道（NEO）フレームワークを通じて統合した、新規のハイブリッド量子・古典アルゴリズムによる第一原理分子動力学（AIMD）を提案している。
- NEOフレームワーク: 選択された軽い原子核（本研究ではプロトン）は、電子と同様に核軌道を用いて量子力学的に扱われる。これは、これらの特定の原子核に対して標準的なボルン・オッペンハイマー近似による分離を回避するものである。
- ハイブリッド動力学アルゴリズム: アルゴリズムは、古典的なボルン・オッペンハイマー動力学に基づく第二量子化のアプローチに従う。
- 量子コンポーネント: 変分量子固有値ソルバー（VQE）を用いて、量子コンピュータ上で瞬時基底状態エネルギー（E0）および原子核に作用する力を決定する。力は、基底状態 ∣ψ0(R)⟩ における力演算子の期待値 ∇IH(R) を評価するヘルマン・ファインマンの定理を用いて計算される。
- 古典コンポーネント: 量子プロセッサから計算された力に基づき、速度ベルレ法を用いて分子の幾何構造を更新する。
- 活性空間手法: 量子ノイズあり中規模（NISQ）デバイスの限られた量子ビット数を対処するため、著者らはハートリー・フォック埋め込みを実装している。量子コンピュータのために活性空間が選択され、不活性軌道はハートリー・フォックレベルで扱われる。選択基準には、単一軌道エントロピー、反応経路に沿った変動性、および対称性の保持が含まれる。
- アンザッツの選択: 本研究では、化学的に動機付けられたユニタリー結合クラスター（UCCSD/NEOUCCSD）アンザッツと、ハードウェア効率的なアンザッツ（Real Amplitudes および Tiled Unitary Product Ansatz - tUPS）の両方を用いて、近未来のデバイスにおける実現可能性を評価している。
- 基底関数系: シミュレーションには、6-31G電子基底関数系と、特定の核基底関数系（DZSPDN：s型のみのDZSNB、およびs型とp型の両方のDZSPNB）を用いている。
主な貢献
- アルゴリズムの実装: 本研究は、量子コンピュータ上でNEO法を用いて核量子効果を明示的に組み込んだ、ハイブリッド量子・古典AIMDアルゴリズムの初の実装を提示している。
- 活性空間戦略: 著者らは、動力学的プロセス中に物理的一貫性を確保するために、軌道の変動性と対称性を考慮したNEO動力学のための特定の活性空間選択プロトコルを開発した。
- ハードウェア効率: 本研究は、コンパクトなハードウェア効率的アンザッツが、より正確なUCCSDアンザッツと同等の結果を達成できることを示しており、これは近未来の量子ハードウェアでの実行可能性を示唆している。
- 基底関数の感度: シミュレーションの精度が、核基底関数系の選択（特に複雑なモード結合のためのp型軌道の必要性）および活性空間に決定的に依存することを強調している。
結果
アルゴリズムは、水素分子（H2）、水（H2O）、およびズンデルイオン（H5O2+）の3つの系を用いたシミュレーションを通じて検証された。
- 水素（H2）: NEOによるNQEの導入は、純粋な電子的手法で見られる振動数の過大評価を補正したが、わずかな過小評価も招いた。ハードウェア効率的アンザッツ（tUPSおよびReal Amplitudes）は、NEOの場合において1%未満の誤差を達成した。これは、アンザッツによる過大評価とNEOによる過小評価の間の幸運な誤差の相殺に起因する。
- 水（H2O）: NEOシミュレーションは、すべての正規モード（曲げ、対称伸縮、反対称伸縮）において、純粋な電子的手法よりも低い相対誤差を示した。プロトン間の相互作用が近い曲げモードは依然として困難であったが、NEOによる改善が見られた。
- ズンデルイオン（H5O2+）:
- ポテンシャルエネルギー面: 純粋な電子的なケースでは、プロトン移動のポテンシャルは障壁を持つ二重井戸構造を示す。NQEを導入すると、プロトンの往来運動は実質的に障壁のないものとなり、ポテンシャルを単一井戸へと変容させる。これは、共有プロトンの零点運動が顕著な影響を及ぼし、振動基底状態が純粋な電子的な障壁よりも高い位置にあることを示している。
- スペクトル: 純粋な電子シミュレーションは、1000 cm−1付近の特徴的なスペクトルの二重線や、プロトン移動と水の揺らぎモードとの結合を捉えることができなかった。
- NEOの成功: p型核軌道を用いたNEOシミュレーションは、プロトン移動モードと水の揺らぎ／伸縮モードとの間の結合およびスペクトルの二重線を正常に再現した。より大きな活性空間の最適化問題を解くために、縮退したNEOUCCSDアンザッツ（63個の励起のうち17個を選択）が必要であった。
- 精度: NEOの結果は、活性空間内の厳密なNEOCASCI解と比較して、ピークあたりの最大誤差が40 cm−1未満であった。また、実験的な解離前スペクトルとも定性的に一致した。
意義と主張
本論文は、彼らのハイブリッドアルゴリズムが、量子コンピュータ上で核量子効果を伴う分子動力学をシミュレートすることの実現可能性を正常に示したと主張している。著者らは以下の点を強調している：
- プロトン移動を正しくモデル化するためには、人工的なエネルギー障壁を取り除き、正しいモード結合を可能にするために、ズンデルイオンにおけるNQEが不可欠である。
- UCCSDは最高の精度を提供するが、ハードウェア効率的なアンザッツは、活性空間と核基底関数系を慎重に選択すれば、近未来のデバイスへの実行可能な道筋を提供する。
- 活性空間と核基底関数系（特にズンデルイオンにおけるp軌道の包含）の選択は、特定のアンザッツの使用よりも精度の向上において重要である。
- 本研究は、より大きなシステムへの将来の研究のための基礎を築くものであり、リソース使用量や精度をさらに向上させるために、カノニカル基底から自然軌道基底へ移行できる可能性があることを述べている。
著者らは、現在の技術の状態に対して謙虚であり、結果がノイズのない状態ベクトルシミュレータに依存していること、および、まだ物理的なハードウェア上で顕著な量子優位性が実証されていないことを認めている。本研究は、量子化学の手法を動的に核自由度へと拡張するための概念実証として機能している。
毎週最高の quantum physics 論文をお届け。
スタンフォード、ケンブリッジ、フランス科学アカデミーの研究者に信頼されています。
受信トレイを確認して登録を完了してください。
問題が発生しました。もう一度お試しください。
スパムなし、いつでも解除可能。
週刊ダイジェスト — 最新の研究をわかりやすく。登録