[論文レビュー] Continual Model of Medium IV: Calculation of analytical gradients for solvation energy over atomic coordinates

Oleg Kupervasser, N. E. Wanner|arXiv (Cornell University)|Aug 8, 2011

Various Chemistry Research Topics参考文献 3被引用数 6

ひとこと要約

本稿では、連続媒質モデル（PCM）を用いて、原子座標に関する溶媒化エネルギーの解析的勾配を提示し、陽性溶媒モデルにおける効率的で高精度な自由エネルギー計算を可能にしている。非極性および極性溶媒化寄与の勾配は、PCM、COSMO、SGBを用いて導出され、溶媒を避ける（SES）および可及的（SAS）表面における表面メッシュパラメータの解析的微分を活用して、ローリングアルゴリズムによって計算されている。

ABSTRACT

In the given paper we describe methods finding analytical gradients (derivatives) of solvation energy over atomic coordinates. It is made both for not polar energy and for the polar energy found by methods PCM, COSMO and SGB. These gradients are searched with use of analytical gradients of parameters (coordinates, normal vectors and squares) the surface meshes. The surface meshes are located on both the optimal smooth molecular surfaces: solvent excluded surface (SES) and surface solvent accessible surface (SAS) received from SES. These surfaces are found by methods of primary and secondary rolling according to the algorithm described in the previous papers of the series.

研究の動機と目的

• 分子力学および幾何最適化に使用可能な、原子座標に関する溶媒化エネルギーの解析的勾配の開発を目的とする。
• 連続媒質モデルを拡張し、非極性および極性寄与の両方の溶媒化エネルギーの微分を計算することを目的とする。
• 座標、法線ベクトル、面積などの表面メッシュパラメータからの勾配を導出することにより、正確で効率的な自由エネルギー計算を可能とすることを目的とする。
• 主および二次のローリングアルゴリズムによって生成された溶媒を避ける（SES）および可及的（SAS）表面に適用することを目的とする。
• 自由エネルギープロバランスやコンformationalサンプリングなど、正確なエネルギー微分を必要とする高度な計算化学手法を支援することを目的とする。

提案手法

• 溶媒を連続的媒質として表現し、解析的表面表現を用いる連続媒質モデルを用いる。
• 主および二次のローリングアルゴリズムを適用して、原子座標から滑らかな溶媒を避ける（SES）および可及的（SAS）表面を計算する。
• 表面メッシュパラメータ（座標、法線ベクトル、表面積）の解析的微分を用いて、溶媒化エネルギーの解析的勾配を導出する。
• 非極性溶媒化エネルギーは表面積に基づき、極性溶媒化エネルギーはPCM、COSMO、SGBモデルを用いて計算する。
• 解析的微分を統合されたフレームワークに組み込み、分子力学および量子化学シミュレーションへの応用を可能とする。
• エネルギー最小化およびダイナミクスにおける数値的安定性と精度を保証するため、一貫性のある解析的微分による手法の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 陽性溶媒モデルにおける原子座標に関する溶媒化エネルギーの解析的勾配は、どのように計算可能か？
• RQ2 座標、法線、面積などの表面メッシュパラメータの微分は、溶媒化エネルギー勾配の計算において果たす役割は何か？
• RQ3 ローリングアルゴリズムは、解析的勾配計算に適した滑らかなSESおよびSAS表面をどのように生成するか？
• RQ4 同じフレームワークは、PCM、COSMO、SGBモデルにおける非極性および極性溶媒化寄与の両方を処理できるか？
• RQ5 解析的勾配は、分子シミュレーションにおける自由エネルギー計算の精度および効率にどのような影響を与えるか？

主な発見

• 本稿では、PCM、COSMO、SGBモデルにおいて、非極性および極性寄与の両方について、原子座標に関する溶媒化エネルギーの解析的勾配を成功裏に導出している。
• SESおよびSAS表面における表面メッシュパラメータの解析的微分を用いて勾配を計算しており、高い数値的精度を実現している。
• 主および二次のローリングアルゴリズムの使用により、滑らかで微分可能な溶媒可及的および避ける表面の構築が可能になっている。
• 本手法は一貫性があり安定したエネルギー微分を提供しており、幾何最適化および分子ダイナミクスシミュレーションに不可欠である。
• フレームワークは溶媒を避ける表面および可及的表面の両方に適用可能であり、溶媒モデルの実装における柔軟性を高めている。
• 導出された勾配は数学的整合性に基づいて検証されており、高度な自由エネルギー計算手法への統合に適している。