Quantum Computing for Electronic Circular Dichroism Spectrum Prediction of Chiral Molecules

技術要約：キラル分子の電子円二色性スペクトル予測のための量子コンピューティング

問題提起
電子円二色性（ECD）分光法は、安全な創薬およびエナンチオ選択的合成に不可欠な、キラル分子の絶対配置を決定するための重要なツールである。しかし、理論計算の計算コストが、ECDの予測モデリングにおける実用的な適用を妨げている。時間依存密度汎関数理論（TD-DFT）や高レベルの波動関数法（例：結合クラスター法）などの従来の第一原理的手法は、分子のサイズや複雑さに比例して指数関数的に増大する広範なコンフォメーション・サンプリングと励起状態計算を必要とする。機械学習は高速化を実現する一方で、複雑な系や強相関系において、信頼できる絶対配置の割り当てに必要な物理的厳密さや転移性を欠くことが多い。したがって、古典的なアルゴリズムの指数関数的なスケーリングを回避しつつ、第一原理の精度を維持できる、スケーラブルで物理学に基づいた計算手法への切実なニーズが存在する。

手法
著者らは、ECDスペクトルを予測するために、量子方程式論（qEOM）形式と統合された変分量子フレームワークを導入している。このワークフローは、マルチGPUおよび量子プロセッシングユニット（QPU）によって加速されたハイブリッド量子・古典環境内で動作する。主要な手法ステップは以下の通りである：

1. 活性空間近似： 量子ビット数の要求を管理するため、著者らは活性空間アプローチを採用している。静的および動的相関を捉えるために不可欠な、分数占有数を示す軌道のサブセットを選択するために、古典的な相関計算（MP2/cc-pVDZ）から得られた自然軌道占有数（NOON）を使用する。
2. ハミルトニアンのマッピング： 活性空間の電子ハミルトニアンは、第二量子化演算子を用いて構築され、ジョルダン・ウィグナー変換を介して量子ビット表現にマッピングされる。これにより、約20〜24量子ビット（CAS(10,10)およびCAS(12,12)活性空間に対応）に作用する量子ビット・ハミルトニアンが得られる。
3. 基底状態の最適化（VQE）： 変分量子固有値ソルバー（VQE）を用いて、相関のある基底状態の波動関数（$|\Psi_0\rangle$）およびエネルギーを決定する。ユニタリ結合クラスター（UCC）に着想を得たアンザッツを用い、ハミルトニアンの期待値を最小化するために古典的最適化アルゴリズム（COBYLA）によって最適化を行う。
4. 励起状態の計算（qEOM）： 量子方程式論（qEOM）形式を用いて励起状態を生成する。励起演算子（単一および二重励起）が、VQEによって準備された基底状態に作用する。励起エネルギーと波動関数の係数は、量子デバイス上で評価された交換子の行列要素を含む一般化固有値問題を解くことで得られる。
5. カイロプティカル特性の計算： 電気双極子モーメント（$\mu_{0\mu}$）および磁気双極子モーメント（$m_{0\mu}$）は、qEOM励起状態を用いて計算される。これらを組み合わせて回転強度（$R_{0\mu}$）を決定し、これが励起エネルギーとともにECDスペクトルを定義する。
6. 実装： シミュレーションには、CUDA-Qフレームワーク内のNVIDIAマルチGPU状態ベクトルシミュレータを利用する。ハミルトニアンの期待値は、パウリ文字列の評価を並列化するために4つの仮想QPUに分散され、大きな活性空間に伴う実行時間のボトルネックを軽減している。

主な貢献

• フレームワークの開発： 著者らは、基底状態および励起状態の両方の性質に対して完全に量子整合的な取り扱いを用いて、電子構造モデリングからECDスペクトルの予測までを繋ぐ、エンドツーエンドの量子・古典ワークフローを提示している。
• スケーラビリティのデモンストレーション： 本手法は、20〜24量子ビットにマッピングされた活性空間を利用して、単一または複数のキラル中心を持つものを含む、12種類の臨床的に関連性の高いキラル薬物分子に適用されている。
• 検証： 量子計算の結果は、確立された古典的手法に対して系統的に検証されている。基底状態エネルギーは、結合クラスター単一励起・二重励起（CCSD）および完全活性空間配置間相互作用法（CASCI）に対してベンチマークされている。励起状態エネルギーおよびECDスペクトルは、同一の活性軌道空間内でのCASCIと比較されている。

結果
本研究は、VQE–qEOMフレームワークが、古典的なCASCI参照計算と定量的レベルで一致することを示している：

• スペクトルの精度： 予測されたECDスペクトルは、単一および多キラル中心分子の両方において、スペクトル線形、コットン効果の符号、および相対的なピーク強度を正確に再現している。
• エナンチオマーの識別： 本フレームワークは、R体とS体の間の鏡像対称性を正しく捉え、カイロプティカル応答の符号反転を正確に再現している。
• パフォーマンス： D-リボース（CAS(12,12)、24量子ビット）やオメプラゾールのような分子において、量子計算の結果は古典的参照値と定量的に一致している。トレオニンのような柔軟な分子では、相関効果への感度の高さに起因して、ピーク強度にわずかな偏差が見られたが、全体的なスペクトル傾向およびエナンチオマーの割り当ては一貫していた。
• 効率性： マルチGPU並列化戦略により、単一QPU実行と比較して実行時間が大幅に短縮され、より大きな活性空間でのシミュレーションが可能となった。

意義および主張
本論文は、本研究が、第一原理に基づくカイロプティカルスペクトルの予測のための実用的かつスケーラブルな経路として、量子アルゴリズムを確立したものであると主張している。VQEの変分的な堅牢性と、方程式論（equation-of-motion）の手法の概念的構造を組み合わせることで、本フレームワークは、データ駆動型モデルや計算負荷の高い古典的手法に代わる、物理学的に根ざした選択肢を提供する。著者らは、本アプローチが、現在のノイズあり中間規模量子（NISQ）ハードウェアの能力と、製薬研究における正確なキラル識別へのニーズとの間の溝を埋めることに成功したと断言している。これらの結果は、量子を用いた分子分光学が、化学的および製薬学的に関連のあるシステムにおいて実現可能であることを示唆しており、量子ハードウェアの進展に伴う、より大きな分子やより複雑な立体化学環境への将来的な拡張のための基礎を提供するものである。