[논문 리뷰] Computing electrostatic potentials using regularization based on the range-separated tensor format

Peter Benner, Venera Khoromskaia외 3명|arXiv (Cornell University)|2019. 01. 28.

Tensor decomposition and applications참고 문헌 34인용 수 5

한 줄 요약

이 논문은 생물분자의 전기적 잠재력을 효율적으로 계산하기 위해 범위분리(RS) 텐서 형식을 사용하여 포아송-볼츠만 방정식(PBE)에 대한 새로운 정규화 기법을 제안한다. 전기적 잠재력을 단거리 및 장거리 성분으로 분리하여 저질서 캐논리컬 텐서 표현을 통해, 특이 딜타 원천을 부드럽고 국소화된 장거리 함수로 대체함으로써, 정규화된 부분에 대한 단일 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 범위분리 텐서 형식을 사용하여 PBE를 정규화할 수 있는가? 이는 특이 원천 항을 피하면서도 해의 정확도를 유지하는가?
• RQ2 RS 텐서 기반 정규화는 고전적 포아송 방정식 및 PBE 모델에 비해 원자핵 근처의 해의 특이성을 얼마나 잘 포착하는가?
• RQ3 제안된 방법은 기존의 FEM/FDM 접근법에 비해 얼마나 더 높은 해상도의 격자에서 정확한 전기적 상호작용 계산을 가능하게 하는가?
• RQ4 장거리 및 단거리 성분의 분리가 분자의 회전에 따른 전기적 잠재력의 효율적 재계산을 가능하게 하는가?
• RQ5 RS 텐서 기반 PBE 해법의 달성 가능한 정확도와 계산 효율성은 저질서 근사 오차와 시스템 해법 비용 측면에서 어떻게 평가되는가?

주요 결과

• RS 텐서 기반 정규화는 PBE를 장거리 성분에 대한 단일 선형 시스템 해법으로 줄여 기존의 다중 해법 접근법에 비해 계산 비용을 크게 감소시킨다.
• 단백질 패시큘린 1에 대해 $257^3$ 격자에서 이산 $L_2$ 오차가 $4.3 \times 10^{-8}$ 이하로 유지되어 높은 정확도를 입증한다.
• RS-랭크 10일 때 $n=129^3$에서 장거리 성분의 근사 오차는 $\text{O}(10^{-5})$이며, $n=257^3$에서 더 강한 절단 임계값을 적용하면 오차는 $\text{O}(10^{-6})$ 수준이 된다.
• 고전적 포아송 방정식은 해의 특이성을 정확히 포착하지 못하지만, RS 텐서 방법은 분해를 통해 이를 정확히 해결한다.
• 장거리 잠재력의 저질서 구조 덕분에 기존 FEM보다 훨씬 더 미세한 격자에서 정확한 계산이 가능하다.
• 비선형 PBE로의 일반화가 자연스럽게 가능하며, 분자의 영역 내 원천 항만 수정되고 비선형 영역은 그대로 유지되기 때문이다.