[論文レビュー] Range-separated tensor representation of the discretized multidimensional Dirac delta and elliptic operator inverse

Boris N. Khoromskij|arXiv (Cornell University)|Dec 6, 2018

Tensor decomposition and applications参考文献 19被引用数 4

ひとこと要約

本稿では、グリーン核の低ランクテンソル分解を用いて、離散化された多次元デルタ関数および楕円型作用素の逆演算のための、作用素に依存する範囲分離（RS）テンソル近似を提案する。グリーン核を短距離および長距離成分に分割し、楕円型作用素を適用することで、$O(n)$の複雑性スケーリングを達成し、多数の特異的電荷を有するポテンシャル方程式の効率的かつ高精度な解法が可能になる。これは、ポアソン＝ボルツマン方程式によるバイオ分子電気的状態のモデル化にとって不可欠である。

ABSTRACT

In this paper, we introduce the operator dependent range-separated tensor approximation of the discretized Dirac delta in $\mathbb{R}^d$. It is constructed by application of the discrete elliptic operator to the range-separated decomposition of the associated Green kernel discretized on the Cartesian grid in $\mathbb{R}^d$. The presented operator dependent local-global splitting of the Dirac delta can be applied for solving the potential equations in non-homogeneous media when the density in the right-hand side is given by the large sum of pointwise singular charges. We show how the idea of the operator dependent RS splitting of the Dirac delta can be extended to the closely related problem on the range separated tensor representation of the elliptic resolvent. The numerical tests confirm the expected localization properties of the obtained operator dependent approximation of the Dirac delta represented on a tensor grid. As an example of application, we consider the regularization scheme for solving the Poisson-Boltzmann equation for modeling the electrostatics in bio-molecules.

研究の動機と目的

• 特異的デルタ関数源を伴う楕円型偏微分方程式を解く際、従来のグリッドベースの手法に生じる数値的不安定性と高い計算コストを解消すること。
• 数値的安定性と精度の向上を目的として、デルタ関数を短距離および長距離成分に分割する、作用素に依存するテンソルベースの分割法を開発すること。
• 非一様媒質に多数の点電荷を有するポテンシャル方程式の解法へと、この手法を楕円型作用素の逆演算に拡張すること。
• 低ランクテンソル形式を用いて、ポアソン＝ボルツマン方程式に従う大規模バイオ分子系のスケーラブルなシミュレーションを可能にすること。
• 集団的デルタ関数の長距離成分が粒子数に従って対数的に増加することを示し、大規模N粒子系において線形対数複雑性を達成できることを実証すること。

提案手法

• 楕円型作用素の離散化グリーン核を、範囲分離（RS）テンソル形式を用いて短距離および長距離成分に分解する。
• 各成分に楕円型作用素 $\mathcal{L}$ を適用し、デルタ関数の対応するRS分割を得る： $\delta = \mathcal{L}p_s + \mathcal{L}p_l = \delta_s + \delta_l$。
• 両方の $\delta_s$ および $\delta_l$ を低ランクテンソル形式で表現し、合計近似誤差が $\varepsilon$-切り捨てランクおよびグリッドサイズ $n$ によって制御されることを保証する。
• $d$ 重の $n \times \cdots \times n$ キャリー図形グリッド上でのニュートン核（$1/\|x\|$）の正準RSテンソル表現を、分解の基盤として用いる。
• 半径関数としての $1/\|x\|$ や $1/\|x\|^3$ が低ランクRSテンソル分解を有することを活用し、流体力学的および弾性ポテンシャルの効率的計算を可能にする。
• 長距離スムージング成分による特異的電荷分布の正則化を用いて、ポアソン＝ボルツマン方程式にこの手法を適用し、条件数の問題を軽減する。

実験結果

リサーチクエスチョン

• RQ1 $ℝ^d$ 内の離散化デルタ関数は、作用素に依存する範囲分離テンソル形式を用いて効率的に近似可能か？
• RQ2 グリーン核のRS分割は、特異的および滑らかな成分を分離する安定的かつ低ランクのデルタ関数表現をもたらすか？
• RQ3 $N$ 個の点電荷に対する集団的デルタ関数の長距離成分は、 $N$ に対してランクが対数的にしか増加しない形で表現可能か？
• RQ4 テンソルベースのアプローチは、多数の特異的源を有するポテンシャル方程式の条件数およびスケーラビリティをどのように改善するか？
• RQ5 定数係数を有する楕円型作用素の逆演算（例えば、バイハーモニック作用素やストークス作用素）へ、この手法を拡張可能か？

主な発見

• 作用素に依存するRS分割されたデルタ関数は、グリーン核のRS分解に楕円型作用素 $\mathcal{L}$ を適用することで構築され、短距離および長距離テンソル成分の和として得られる。
• 長距離成分 $\delta_l$ は滑らかであり、低ランクテンソル近似に適しており、短距離成分 $\delta_s$ は局所化された補正を特徴としている。
• テンソル形式のおかげでグリッドサイズ $n$ に対して $O(n)$ の複雑性スケーリングを達成し、大規模3次元グリッド上での効率的計算が可能になる。
• 集団的デルタ関数の長距離成分のランクは、粒子数 $N$ に対してのみ対数的に増加するため、大規模 $N$ 粒子系のスケーラブルな処理が保証される。
• 数値的実験により、近似の局在化および精度が確認され、誤差は $\varepsilon$-切り捨ておよびグリッドサイズに依存して制御される。
• バイオ分子電気的状態のポアソン＝ボルツマン方程式に、このアプローチを成功裏に適用し、点電荷特異性に起因する数値的不安定性を回避する正則化スキームを提供した。