Efficient, precise DFT calculations of NMR shieldings: Revisiting the finite field approach

化学の世界において、分子の正確な形状を知ることは、その仕組みを理解するための鍵となることが多い。ある物質が救命薬となるか、あるいは毒性のある汚染物質となるかは、原子の精密な配置によって決まる。これらの目に見えない構造を見るために、科学者は核磁気共鳴（NMR）と呼ばれる強力なツールに頼っている。この技術は、分子内部の原子の周囲にある磁気環境を検知する、分子の高解像度指紋スキャナーのような役割を果たす。分子が強い磁場の中に置かれると、その原子は独自の反応を示し、それによって隣接する原子や空間的な関係性が明らかになる。しかし、これらの信号を解釈することは必ずしも単純ではない。複雑な分子はしばしば重なり合った信号を生み出し、それらを解きほぐすことは困難であり、データに適合する可能性のある形状が複数残ってしまうこともある。このパズルを解くために、科学者はコンピュータシミュレーションを用いて、さまざまな候補構造に対してNMR信号がどのようになるかを予測する。コンピュータの予測が現実世界の実験と一致すれば、その構造が確定する。

長年の課題は、これらのコンピュータ予測を高い精度で計算することが極めて困難であることだった。磁場に対して分子がどのように反応するかをモデル化するために必要な数学は複雑であり、多くの高度な手法において、これらの計算を行う標準的な方法は、あまりに遅すぎるか、あるいは単に利用不可能である。これが、生物学や医学において極めて重要な、大きく複雑な分子の研究におけるボトルネックとなっている。研究チームは今回、これらの計算を実行するための、新しい効率的な方法を開発した。巧妙な数学的ショートカットと、作業を多くのコンピュータプロセッサに分割して行う戦略を組み合わせることで、彼らは高速かつ高精度な手法を作り出した。彼らのアプローチにより、科学者は、最も確立されてはいるがはるかに遅い手法に匹敵するレベルの精度で、大規模で現実的な化学系のNMR信号を予測することが可能になった。

研究者たちは、電子が分子内でどのように振る舞うかを予測するための、現代の計算化学における主力である「密度汎関数理論」と呼ばれる特定の種類の計算に焦点を当てた。伝統的に、磁気遮蔽（原子の信号がNMRスペクトル上のどこに現れるかを決定する特性）を計算するには、磁場が適用されたときに電子雲がどのようにシフトするかを記述する一連の困難な方程式を解く必要がある。これらの方程式は、特に最も正確なバージョンの理論を使用する場合、大規模なシステムに対して直接解くには複雑すぎることが多い。チームは、これらの複雑な方程式を直接解こうとする代わりに、磁場に小さな制御された変化を与え、システムのエネルギーがどのように変化するかを見ることで、答えを近似できることに気づいた。これは、地図から傾斜を計算しようとするのではなく、丘を2歩登ってみて、どれくらい高くなったかを見ることで、丘の傾斜を測定することに似ている。

この「ステップ・バイ・ステップ」のアプローチを必要な精度で機能させるために、チームは2つの大きな障害を克服しなければならなかった。第一に、計算には複素数が含まれる。複素数とは、虚数成分を含む実数の数学的な拡張である。標準的な化学計算は通常、実数を扱うが、磁場の存在によって数学は複素数の領域へと強制される。第二に、電子間の相互作用を計算するプロセスは計算コストが高く、膨大なコンピュータ時間を必要とすることが多い。研究者たちは、複素数を効率的に処理するバージョンの計算を実装することで、第一の問題を解決した。第二の問題については、「分解識別近似（resolution-of-identity approximation）」と呼ばれる手法を用いた。この手法は、電子間の相互作用の計算を、より小さく扱いやすい断片に分解することで簡略化するものであり、これは大規模な建設プロジェクトが、同時に実行できる小さなタスクに分割される様子に似ている。

この研究における最も重要な革新は、コンピュータ作業の組織化の仕方であった。彼らは、単一のプロセッサ上で計算を逐次実行するのではなく、タスクを複数のプロセッサ間で並列に分割して実行するシステムを設計した。磁場は異なる方向にも適用できるため、各方向に対する計算は互いに独立しているということに彼らは気づいた。これにより、各方向に異なるプロセッサを割り当てることができ、それらは互いに絶えず通信する必要なく、同時に作業を進めることができた。共有メモリと分散メモリコンピューティングの両方を使用するこのハイブリッドなアプローチにより、コンピュータのパワーを増強するにつれて、答えを得るまでに要する時間が大幅に減少した。対照的に、従来の手法では、プロセッサ同士が互いに待ち時間を要するため、プロセッサを追加しても計算が速くならないという壁に突き当たることがよくある。

チームは、単純なベンゼン環のような分子から、高コレステロール治療に使用されるアトルバスタビンのような大型で複雑な薬物に至るまで、多様な分子を用いてこの新手法をテストした。彼らは、正確ではあるものの、大規模なシステムに対して実行するのが極めて困難で遅いことで知られる、分析手法のゴールドスタンダードと比較した。結果は、彼らの新手法がゴールドスタンダードと実質的に同一の数値を生み出し、その誤差は実験測定の自然な不確かさに比べれば無視できるほど小さいことを示した。より小さな分子については、新手法は従来のアプローチよりも高速であった。数千の原子を含む最も大きな分子については、新手法はコンピュータクラスター上で実行した場合、従来のアプローチを大幅に上回るスピードで、劇的に高速化した。これは、彼らの手法が効率的にスケールすること、つまり、パフォーマンスの天井に突き当たるのではなく、計算能力が増すにつれて性能が向上することを証明した。

この研究の実用的な価値を証明するために、研究者たちは化学における有名な歴史的ミステリー、すなわち天然物であるヘキサシクリノール（hexacyclinol）の構造にこの手法を適用した。この分子はシベリアの白樺から発見された菌類から単離されたものだが、その正確な形状は長年激しい論争の対象となっていた。異なる研究グループが異なる構造を提案したが、2006年にこの分子の全合成によって正しい形状が確認されるまで、その論争は終わらなかった。研究者たちは、新しいアルゴリズムを使用して、最初に提案された構造と、修正された正しい構造の両方のNMR信号を計算した。その結果、修正された構造が、元の提案よりも実験データに格段によく一致していることが明確に示された。この成功した適用は、彼らの手法が単なる理論的な改善ではなく、かつて科学者を悩ませてきた現実世界の構造問題を解決できる実用的なツールであることを裏付けた。

また、研究チームは、新しい手法の速度が精度の犠_を伴わないことを確認するために、誤差の発生源を慎重に調査した。彼らは、計算が高速かつ精密である「スイートスポット」を見つけるために、磁場のステップのサイズやさまざまな数学的近似をテストした。その結果、特定の最適な磁場強度を使用することで誤差を最小限に抑え、実験的な不確かさの閾値を十分に下回るレベルに保てることを発見した。また、ソフトウェア内の一般的なショートカットのいくつかが重い原子に対して小さな誤差を導入する可能性があることも発見したが、それらの誤差は管理可能であり、より堅牢な設定を使用することで回避できることも分かった。この細部へのこだわりにより、この手法が軽い元素だけでなく、幅広い化学元素に対して信頼できるものであることが保証されている。

大規模で複雑な系に対して高精度なNMR計算を身近なものにすることで、本研究は、これまで分析が困難であった分子の研究への扉を開いた。分子が磁場に対してどのように振る舞うかを迅速かつ正確に予測できる能力は、研究者がより高い確信を持って構造を確認し、新薬を設計し、生物学的プロセスを理解することを可能にする。この手法は、厳密な誤差分析と効率的な並列コンピューティングの基盤の上に構築されており、現代の化学研究における増大する複雑さに対処できるものである。著者たちが述べているように、このアプローチはNMRに限定されるものではない。独立した計算をプロセッサ間で分割するという同じ戦略は、他の困難な化学の問題にも適用可能であり、科学者が分子の世界をモデル化する方法を変革する可能性がある。この研究は、計算タスクの組織化の仕方を再考することで、かつては遅すぎる、あるいは複雑すぎて対処できないと考えられていた問題を解決できることを示す実証となっている。