[論文レビュー] Unconditionally stable time splitting methods for the electrostatic analysis of solvated biomolecules

Leighton Wilson, Zhao, Shan|arXiv (Cornell University)|Oct 10, 2014

Electromagnetic Simulation and Numerical Methods参考文献 30被引用数 5

ひとこと要約

本稿では、バイオ分子電気双極子系における時間依存非線形ポアソン＝ボルツマン（NPB）方程式を解くために、特に局所的一様陰解法（LODIE）スキームを含む、無条件安定な作用素分割法を提案する。非線形項を解析的に統合し、有限差分を用いた陰解法時間刻みを採用することで、無条件安定性とタンパク質サイズに比例する線形スケーラビリティを達成し、リチャードソン補外と特化した回復スキームによる加速により、大きな時間刻みでも効率的かつ高精度な自由エネルギー計算が可能になる。

ABSTRACT

This work introduces novel unconditionally stable operator splitting methods for solving the time dependent nonlinear Poisson-Boltzmann (NPB) equation for the electrostatic analysis of solvated biomolecules. In a pseudo-transient continuation solution of the NPB equation, a long time integration is needed to reach the steady state. This calls for time stepping schemes that are stable and accurate for large time increments. The existing alternating direction implicit (ADI) methods for the NPB equation are known to be conditionally stable, although being fully implicit. To overcome this difficulty, we propose several new operator splitting schemes, in both multiplicative and additive styles, including locally one-dimensional (LOD) schemes and additive operator splitting (AOS) schemes. The proposed schemes become much more stable than the ADI methods, and some of them are indeed unconditionally stable in dealing with solvated proteins with source singularities and non-smooth solutions. Numerically, the orders of convergence in both space and time are found to be one. Nevertheless, the precision in calculating the electrostatic free energy is low, unless a small time increment is used. Further accuracy improvements are thus considered. After acceleration, the optimized LOD method can produce a reliable energy estimate by integrating for a small and fixed number of time steps. Since one only needs to solve a tridiagonal linear system in each independent one dimensional process, the overall computation is very efficient. The unconditionally stable LOD method scales linearly with respect to the number of atoms in the protein studies, and is over 20 times faster than the conditionally stable ADI methods.

研究の動機と目的

• 溶媒中バイオ分子における非線形ポアソン＝ボルツマン（NPB）方程式に適用可能な従来の条件付き安定な交替方向陰解法（ADI）の不安定性を解消すること。
• 非常に大きな時間刻みでも安定を保つ、無条件安定な時間積分スキームの開発。特に滑らかでない解や特異的源項を含む場合にも有効であること。
• 非常に小さな時間刻みを用いない限り、電気的自由エネルギー推定精度が低くなる問題の改善。
• 三重対角行列の反復解法とタンパク質サイズに比例する線形スケーリングを活用した計算効率の向上。
• リチャードソン補外と作用素分割に基づく回復スキームを用いて収束性と精度を向上させ、真空中ではFFTを置き換えること。

提案手法

• NPB方程式を一次元部分問題に分解するため、乗法的および加法的作用素分割スキーム（特に局所的一様（LOD）および加法的作用素分割（AOS））を提案する。
• NPB方程式に含まれる非線形双曲正弦項を解析的に統合することで、非線形不安定性を排除し、無条件安定性を達成する。
• 中心差分による空間離散化と、陰解法オイラーおよびクランク＝ニコルソン時間積分を組み合わせ、時間方向に1次精度を保証しながら安定性を確保する。
• 各方向の交互に解く際、トマス法を用いて三重対角行列系を効率的に解き、原子数に比例する線形スケーリングを実現する。
• リチャードソン補外と、FFTを作用素分割に置き換えた回復スキームを適用し、時間刻みサイズを小さくせずにエネルギーの精度を向上させる。
• 球状キャビティと実際のタンパク質系を用いた数値実験により、安定性と収束性を検証し、誘電率の不連続性や電荷特異項を含む状況でも有効性を確認した。

実験結果

リサーチクエスチョン

• RQ1 時間依存非線形ポアソン＝ボルツマン方程式に対して、滑らかでない解や特異的源項を含んでも、無条件安定性を達成できる作用素分割法を設計できるか？
• RQ2 実際のバイオ分子系において、本手法は従来の条件付き安定AD Iスキームと比較して、安定性と精度に優れているか？
• RQ3 非常に小さな時間刻みに依存せずに、電気的自由エネルギー推定の精度を顕著に向上させられるか？
• RQ4 本手法の計算複雑度と、大規模タンパク質系における原子数に対するスケーラビリティはいかほどか？
• RQ5 リチャードソン補外と特化した回復スキームは、効率を損なわせることなく、精度を効果的に向上させられるか？

主な発見

• 提案されたLODIE法は、すべてのテスト時間刻み（Δt = 0.4を含む）で無条件安定であり、実タンパク質系の滑らかでない解に対しても安定である。
• LODIE法は原子数に比例する線形スケーリングを達成しており、大規模タンパク質では従来の条件付き安定AD I法の20倍以上も高速である。
• 提案されたすべてのスキームは、中心差分と1次時間分割のため、空間および時間方向に1次収束を示すが、加速処理が施されない限り自由エネルギー推定の精度は低い。
• リチャードソン補外と回復スキームを適用した加速LODIE法は、わずかで固定された時間ステップ数のみで信頼性のある電気的自由エネルギー推定が可能になる。
• 全体の計算複雑度はO(N³)のままであるが、三重対角行列解法の活用により、大規模バイオ分子系に対しても高い効率と実用的妥当性を確保している。
• すべてのスキームの中で、LODIE法が安定性、精度、効率のバランスが最も優れており、実タンパク質の電気的自由エネルギー解析に推奨される。