[論文レビュー] Continual Model of Medium III: Calculation of analytical gradients of parameters of surface meshes on the molecular surfaces over atomic coordinates

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

Various Chemistry Research Topics参考文献 2被引用数 7

ひとこと要約

本稿では、分子の溶媒回避（SES）表面および溶媒可及（SAS）表面における表面メッシュパラメータの原子座標に関する解析的勾配を計算する手法を提示する。連続的媒質モデルを用いて、微分幾何学および陰関数微分を用いて、メッシュ幾何（位置、法線、面積）の導関数の閉形式表現を導出する。これにより、連続的溶媒化モデルにおけるエネルギー勾配計算が効率的かつ正確に可能になる。

ABSTRACT

The problem of finding of analytical gradients (derivatives over atoms coordinates) of solvation energies can be decomposed on two subtasks: at the first stage we search for parameters of the superficial devices (three coordinates, three components of normal vector and square) and their derivatives; at the second stage we differentiate energy and we express it through derivative of the matrixes featuring a problem. But these derivatives of matrix elements can be expressed through derivatives of parameters of surface meshes of SES (solvent excluded surface) or SAS (solvent accessible surface). The purpose of the given paper is finding of these analytical gradients of the parameters of the surface meshes.

研究の動機と目的

• 連続的溶媒化モデルにおける溶媒化エネルギー勾配の正確かつ効率的な計算を可能にするために、表面メッシュパラメータの解析的導関数を導出すること。
• 分子シミュレーションにおける原子座標に関する複雑な表面表現（SES/SAS）の微分の課題に対処すること。
• 溶媒化エネルギー勾配計算を2段階に分解すること：第一に、メッシュパラメータの導関数を計算すること；第二に、エネルギー導関数をこれらのメッシュ導関数を通じて表現すること。
• 微分幾何学および陰関数微分を用いた、表面メッシュ要素の勾配計算の数学的に厳密なフレームワークを提供すること。

提案手法

• 暗黙の曲面上での微分幾何学を用いて、頂点位置、法線ベクトル、面積といった表面メッシュパラメータの導関数の解析的表現を導出する。
• 連鎖則および陰関数微分を適用し、溶媒回避または溶媒可及表面定義を介して、原子座標に関するメッシュパラメータの勾配を計算する。
• 分子表面のレベルセット表現を用いて、表面幾何を原子座標およびその導関数の関数として表現する。
• 原子座標に関するメッシュパラメータのヤコビ行列を計算する形式的枠組みを構築し、表面メッシュを通じた勾配伝搬を可能にする。
• エネルギー勾配をメッシュパラメータの導関数の寄与に体系的に分解する。
• 数学的導出により手法の妥当性を検証し、SESおよびSAS表現の両方に適用可能なフレームワークを提供する。

実験結果

リサーチクエスチョン

• RQ1 分子表面における表面メッシュパラメータ（位置、法線、面積）の原子座標に関する導関数を、どのように解析的に計算できるか？
• RQ2 溶媒化エネルギー勾配の文脈において、SES や SAS のような複雑な表面表現を微分するための数学的枠組みは何か？
• RQ3 溶媒化エネルギーの勾配を、原子座標から直接ではなく、表面メッシュパラメータの導関数を通じてどのように表現できるか？
• RQ4 暗黙の分子表面におけるメッシュ要素の導関数の閉形式表現は何か？
• RQ5 提案手法は、一貫性のある数学的取り扱いのもとで、溶媒回避表面および溶媒可及表面の両方へ一般化可能か？

主な発見

• 本稿では、分子表面における表面メッシュパラメータ（頂点位置、法線ベクトル、面積）の原子座標に関する導関数の閉形式解析的表現を導出した。
• 微分幾何学を用いて表面メッシュを通じて導関数を伝搬させることで、連続的溶媒化モデルにおけるエネルギー勾配の正確な計算が可能になる。
• 導出された勾配は、溶媒回避（SES）表面および溶媒可及（SAS）表面の両方に対して有効であり、広範な適用可能性を有する。
• 手法は2段階に分ける：第一に、メッシュパラメータの導関数を計算すること；第二に、エネルギー導関数をこれらのメッシュ導関数を用いて表現すること。
• 形式的枠組みは数学的に整合性があり、分子モデリングおよび自由エネルギー計算における効率的な勾配ベース最適化を支援する。
• 本フレームワークは、連続的溶媒化モデルにおける高精度かつ高効率なアルゴリズムの基盤を提供する。