Towards more accurate natural orbital functional approximations: including 4-index cumulant contributions

あなたは、分子の完璧なデジタルツインを構築しようとしていると想像してください。それは、電子というレゴブロックで作られた小さな構造物です。何十年もの間、科学者たちは「還元密度行列汎関数理論（RDMFT）」と呼ばれるツールを使用してきました。RDMFTは、個々の電子を一つずつ追跡する必要なく、電子がどのように振る舞うかを教えてくれる、非常にスマートなレシピ本のようなものです。個々の電子を追跡することは、風が吹いている中でビーチの砂粒を一つ一つ数えようとするようなものです。目標は、分子がどのように結合し、そして極めて重要なことに、それらがどのように壊れるかを予測することです。

分子が壊れるとき、それは2つの別々の破片に分かれます。ゴムバンドをパチンと切るようなものです。現実の世界では、一度切れてしまえば、その2つの破片は独立した、穏やかな島のように振る舞うはずです。彼らはもう魔法のように繋がっているわけではなく、その内部の「スピン」（量子的な性質であり、小さな方位磁石のようなもの）は、特定の静かな状態に落ち着くはずです。しかし、現在の私たちのレシピの多くは、これを間違えてしまいます。シミュレーションの中で分子が壊れる際、破片が奇妙に繋がり続けていたり、内部のコンパスが制御不能に激しく回転したりすることがあり、それがエネルギーや安定性に関する不正確な予測につながります。これは大きな問題です。なぜなら、結合がどのように壊れるかを理解することは、新しい電池の設計から命を救う薬の開発に至るまで、あらゆることに不可ло欠だからです。もしシミュレーションが「破断」を誤れば、モデル全体の信頼性が揺らいでしまいます。

この論文は、この特定の悩みに取り組んでいます。著者たちは、PNOF5と呼ばれる人気のレシピを修正しました。著者たちは、PNOF5が多くのことには優れているものの、分子が壊れるために引き伸ばされている最中の、電子のペアとしての相互作用に関する決定的な細部を見落としていることに気づきました。それはまるで、レシピが電子は直接のパートナーとだけ会話していると仮定しており、近くにある他のペアともささやき合っている事実を無視しているかのようです。著者たちは、レシピが「破断のルール」を尊重するように強制する新しい手法を導入しました。つまり、分子が分裂するとき、破片が真に分離し、そのスピンが正しく落ち着くようにすることです。これらの厳格なルールを加え、さらに物理法則に従うように結果を磨き上げることで、彼らは改良された新しいバージョンを作り上げました。彼らのシミュレーションによれば、窒素（N2）や酸素（O2）のような分子において、この新しい手法は最も高度で高価なシミュレーションと全く同じ挙動を示し、ついに結合の「スナップ（破断）」を正しく捉えることに成功しました。

折れたゴムバンドの物語

電子の世界へ飛び込んでみましょう。そこでは物事が少し奇妙になります。分子を、電子がダンサーであるダンスフロアだと想像してみてください。安定した分子の中では、彼らは密接に、かつ調和のとれたペアで踊っています。しかし、分子を引き離し、結合が切れるまで引き伸ばしたらどうなるでしょうか？これは「結合解離」と呼ばれます。現実の世界では、結合が壊れた後、2つの半分は独立した、幸せな原子として漂っていくはずです。彼らはもう手をつないでいるわけではなく、彼らの内部の「スピン」（これらを回転するコマだと考えてください）は、穏やかで予測可能なリズムに落ち着くはずです。

問題は、分子のこの挙動を予測しようとする現在のコンピュータモデルが、しばれてこの「破断シーン」を間違えてしまうことです。それらは、幕が降りた後でも俳優たちがステージを去るのを拒んでいる悪い映画のようです。分子の2つの半分は奇妙に繋がり続け、スピンは混乱し、その結果、結合を切るのにどれだけのエネルギーが必要かについて誤った答えを導き出します。これは、新しい材料を設計したり、化学反応を理解しようとしたりしている科学者にとって、大きな悩みの種です。

ヴァレリー・チュイコ、ポール・W・エアーズ、エドゥアルド・マティトのチームが登場しました。彼らは、PNOF5と呼ばれる特定のレシピを修正することに決めました。PNOF5を、電子がどのようにペアを作るかを示す指示書と考えてください。これは巧妙なレシピですが、盲点があります。それは、電子は自分のダンスパートナーのことだけを気にかけ、フロア上の他の誰もを無視するという仮定です。これは分子が完全な状態であれば問題ありませんが、分子を引き離し始めると、その仮定は崩れ去ります。このレシピは、分子が解離する際に必要な、異なる電子ペア間の微妙な遠距離の「ささやき」を考慮できていないのです。

解決策：足りないルールを加える

著者たちの解決策は、レシピに新しい「交通ルール」を加えることでした。彼らは、分子が壊れるとき、破片のスピンがどのように振る舞うべきかについて、具体的で譲れないルールがあることに気づきました。もし分子が2つの破片に分かれるなら、各破片の総スピンは、単独で孤立した原子が持つものと一致しなければなりません。それは、おもちゃの車を真っ二つに割ったとき、それぞれの半分の車輪がランダムに回ったり止まったりするのではなく、正しい速度で回転していなければならないと言うようなものです。

PNOF5を修正するために、彼らは単に推測したのではなく、これらの新しいスピンのルールに従う「可能な限り最も近い」電子マップを見つけるために数学を用いました。割れた花瓶のぼやけた写真があり、その破片が「どうあるべきか」を正確に知っている状況を想像してください。彼らは数学的なアルゴリズムを使用して、その写真を鮮明にし、既知の破断ルールに一致するように強制しました。しかし、注意点がありました。写真を特定の見た目に強制しようとすると、誤って物理法則を壊してしまう（例えば、花瓶の破片を浮かせたりしてしまう）ことがあります。そのため、彼らは第2のステップとして「精製（purification）」プロセスを追加しました。このステップは、新しい、より鮮明になった写真が、依然として宇宙の基本法則（具体的には「N-表現可能性」と呼ばれるもの。これは、その電子配置が実際に現実として存在し得るか、という意味の難しい言葉です）に従っているかどうかをチェックするものです。

彼らは、この「スピンのルールを修正し、次に物理を精製する」という2段階のプロセスを、数値が変化しなくなるまで何度も繰り返しました。それはギターの弦をチューニングするようなものです。締め、ピッチを確認し、テンションを調整し、完璧な音が出るまで繰り返します。

結果：スナップを正しく捉える

彼らが窒素（$N_2$）、酸素（$O_2$）、一酸化窒素（$NO^+$）のような分子を用いてこの新しい改良手法をテストしたところ、驚くべき結果が得られました。「破断ゾーン」において（例えば、原子が3.5 Å以上の距離まで引き離されたとき）、彼らの新しい手法は、コンピュータ化学のゴールドスタンダードとして知られるCASSCFとほぼ完璧に一致するエネルギー値を出しました。

例えば、窒素（$N_2$）のシミュレーションでは、古いPNOF5のレシピは -107.315 ハルティー（エネルギーの単位）を与えましたが、新しい修正版は -107.438 ハルティーを与えました。ゴールドスタンダードであるCASSCFの計算も -107.438 ハルティーでした。これは完璧な一致です。酸素や一酸化窒素についても同様でした。新しい手法はエネルギーを正しく捉えただけでなく、電子が2つの破片間でどの程度共有されているかの尺度である「非局在化指数」も修正しました。古いレシピでは、破片は本来離れているべき時でも、まだ電子を共有していました。新しいレシピでは、その共有は消滅しました。まさに現実の世界と同じように。

しかし、著者たちはこれがすべてに対する魔法の杖であるとは主張していません。彼らは、この手法は分子が実際に壊れているときに最もよく機能することを指摘しています。もし分子がまだ完全な状態で心地よく座っているなら、古いレシピでもすでに十分な仕事をしており、新しいレシピは大きな優位性をもたらしません。また、一酸化炭素（$CO$）のような非常にトリッキーな分子については、数学を機能させるためにセットアップをわずかに微調整する必要があったことも述べています。これは、正しい答えを得るために、原子が特定のやり方で電子を共有するように強制したことを意味します。

なぜこれが重要なのか

この研究は、レシピ全体を捨て去ることなく、コンピュータモデルをより正確にできることを示しているため、大きな前進です。電子がどのように振る舞うかについての、適切に選ばれたルールをいくつか加えるだけで、彼らは長年の問題を解決しました。これは、化学反応を正確にシミュレートしたいのであれば、破断を特別な、乱れた例外として扱うのではなく、システム全体を形作る厳格なルールとして扱う必要があることを示唆しています。

著者たちは、まだやるべきことはあると認めています。彼らの手法は計算量が多く、特に大きな分子に対しては多くのコンピュータパワーを必要とします。また、非常に複雑な状況においては、現在の彼らのアプローチではすべての詳細を捉えきれない可能性があることも述べています。しかし、彼らが切り開いた道は明確です。もし分子がどのように壊れるかを予測したいのであれば、モデルが破片の物理学を尊重するようにしなければなりません。そして、この新しい手法によって、彼らは結合の「スナップ」を正確に捉えることが可能であることを証明したのです。