[論文レビュー] Computing electrostatic potentials using regularization based on the range-separated tensor format

Peter Benner, Venera Khoromskaia+3 more|arXiv (Cornell University)|Jan 28, 2019

Tensor decomposition and applications参考文献 34被引用数 5

ひとこと要約

本稿では、生体分子内の電気的ポテンシャルを効率的に計算するために、範囲分離（RS）テンソル形式を用いたポアソン＝ボルツマン方程式（PBE）の新規正則化スキームを提案する。電気的ポテンシャルを低ランクカノニカルテンソル表現を用いて短距離成分と長距離成分に分解することで、特異的であるディラック・デルタ源を滑らかで局在化した長距離関数に置き換える。これにより、正則化された部分の数値解法を1回のFDM/FEM解法で行い、事前に計算された短距離寄与を直接加算することで、計算コストを低減しつつ、細かいグリッド上でも高い精度を達成できる。

ABSTRACT

In this paper, we apply the range-separated (RS) tensor format [6] for the construction of new regularization scheme for the Poisson-Boltzmann equation (PBE) describing the electrostatic potential in biomolecules. In our approach, we use the RS tensor representation to the discretized Dirac delta [21] to construct an efficient RS splitting of the PBE solution in the solute (molecular) region. The PBE then needs to be solved with a regularized source term, and thus black-box solvers can be applied. The main computational benefits are due to the localization of the modified right-hand side within the molecular region and automatic maintaining of the continuity in the Cauchy data on the interface. Moreover, this computational scheme only includes solving a single system of FDM/FEM equations for the smooth long-range (i.e., regularized) part of the collective potential represented by a low-rank RS-tensor with a controllable precision. The total potential is obtained by adding this solution to the directly precomputed rank-structured tensor representation for the short-range contribution. Enabling finer grids in PBE computations is another advantage of the proposed techniques. In the numerical experiments, we consider only the free space electrostatic potential for proof of concept. We illustrate that the classical Poisson equation (PE) model does not accurately capture the solution singularities in the numerical approximation as compared to the new approach by the RS tensor format.

研究の動機と目的

• 生体分子電気的ポテンシャルにおける特異的ディラック・デルタ源に起因する従来のPBEソルバーの数値不安定性と高コストを解消すること。
• 点電荷の近傍で粗いグリッド分解能を発揮しにくい古典的有限要素法の限界を克服すること。
• 界面処理を複雑にせず、溶質-溶媒界面におけるコーシー境界データの連続性を維持する正則化技術を開発すること。
• 長距離電気的相互作用に低ランクテンソル構造を活用することで、細かいグリッド上での正確かつ効率的な解法を可能にすること。
• 短距離成分と長距離成分を分離することで、分子の回転に伴う電気的ポテンシャルの再計算を効率的に行えるようにすること。

提案手法

• 3次元直交座標グリッド上において、範囲分離（RS）テンソル形式を用いてニュートンポテンシャルを局在化した短距離成分と滑らかな長距離成分に分解する。
• 長距離成分を、低ランクカノニカルテンソルとして、縮小された高階特異値分解（RHOSVD）を用いたカノニカルからタッカー（C2T）変換により表現する。
• PBEにおける特異的ディラック・デルタ分布を、長距離成分の滑らかで低ランクの近似に置き換え、標準的なFDM/FEMによる安定な数値解法を可能にする。
• 短距離寄与を局所的サポートを持つ低ランクテンソルとして事前に計算し、PBEの繰り返し解法を回避する。
• 数値的に解いた長距離成分と事前に計算された短距離テンソル表現を加算することで、全ポテンシャルを構築する。
• PBEの右辺にRSテンソル形式を適用し、修正された源項が分子領域内に局在化し、界面連続性が保持されることを保証する。

実験結果

リサーチクエスチョン

• RQ1 範囲分離テンソル形式を用いた正則化により、特異的源項を回避しつつ解の精度を維持できるか？
• RQ2 RSテンソルに基づく正則化は、古典的ポアソン方程式やPBEモデルと比較して、原子核付近の解の特異性をどの程度正確に捉えることができるか？
• RQ3 提案手法により、従来のFEM/FDM手法と比較して、どの程度細かいグリッド上で正確な電気的ポテンシャル計算が可能になるか？
• RQ4 長距離成分と短距離成分の分離により、分子の回転に伴う電気的ポテンシャルの再計算がどの程度効率的に行えるか？
• RQ5 RSテンソルに基づくPBEソルバーの実現可能な精度と計算効率は、低ランク近似誤差とシステム解法コストの観点からどの程度か？

主な発見

• RSテンソルに基づく正則化により、PBEは長距離成分の1次元線形系の解法に簡略化され、複数回の解法を要する手法と比較して計算コストが顕著に低減される。
• プロテインFasciculin 1において、$257^3$グリッド上での離散的$L_2$誤差が$4.3 \times 10^{-8}$未満に抑えられ、高い精度が達成された。
• 長距離成分は、$n=129^3$でRSランク10で$ ext{O}(10^{-5})$の誤差、$n=257^3$でより厳密な切り捨て閾値を用いることで$ ext{O}(10^{-6})$の誤差で近似可能である。
• 古典的ポアソン方程式は解の特異性を正確に捉えられないが、RSテンソル法は分解により正しく特異性を解明する。
• 長距離ポテンシャルの低ランク構造のおかげで、従来のFEMよりもはるかに細かいグリッド上での正確な計算が可能である。
• 非線形PBEに対しても自然に一般化可能であり、分子領域内の源項のみを変更すればよく、非線形領域はそのままである。