[논문 리뷰] Unconditionally stable time splitting methods for the electrostatic analysis of solvated biomolecules

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

Electromagnetic Simulation and Numerical Methods참고 문헌 30인용 수 5

한 줄 요약

이 논문은 생물분자의 전기적 성질을 분석하기 위해 시간에 따라 변하는 비선형 푸아송-볼츠만(NPB) 방정식을 해결하기 위해 조건부 안정성의 연속된 방향 은폐(ADI) 방법의 문제점을 해결하고자, 특히 국소 1차원 은밀한 오일러(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 방정식을 1차원 하위문제로 분리하기 위해 곱셈형 및 덧셈형 연산자 분할 기법, 특히 국소 1차원(LOD) 및 덧셈형 연산자 분할(AOS)을 제안한다.
• NPB 방정식의 비선형 hyperbolic sine 항을 해석적으로 통합하여 비선형 불안정성을 제거하고 조건부 안정성을 달성한다.
• 중앙 유한 차분을 사용한 은밀한 오일러 및 크랭크-니콜슨 시간 통합 기법을 사용하여 안정성과 시간에 대한 1차 정확도를 확보한다.
• 각 교차 방향에서 삼중 대각선 시스템을 효율적으로 해결하기 위해 토머스 알고리즘을 적용하여 원자 수에 비례하는 선형적 확장성을 확보한다.
• 리치어드슨 보간법과 FFT를 대체하는 연산자 분할 기반 복구 기법을 적용하여 시간 스텝 크기를 줄이지 않아도 에너지 정확도를 향상시킨다.
• 구형 공동체와 다양한 비극성 불연속성 및 전하 특이성을 가진 실제 단백질 시스템에서의 수치 실험을 통해 안정성과 수렴성의 유효성을 검증한다.

실험 결과

연구 질문

• RQ1 비선형 푸아송-볼츠만(NPB) 방정정의 시간에 따라 변하는 방정식에 대해, 부드럽지 않은 해와 특이한 소스 항이 존재할 경우에도 조건부 안정성 없는 안정성을 확보할 수 있는 연산자 분할 방법을 설계할 수 있는가?
• RQ2 실제 생물분자 시스템에서 기존의 조건부 안정성 있는 ADI 방법과 비교해 볼 때, 제안된 방법의 안정성과 정확도는 어떠한가?
• RQ3 매우 작은 시간 스텝에 의존하지 않고도 전기적 자유 에너지 추정의 정밀도를 크게 향상시킬 수 있는가?
• RQ4 제안된 방법의 계산 복잡도와 대규모 단백질 시스템에서 원자 수에 비례하는 확장성은 어떠한가?
• RQ5 리치어드슨 보간법과 특화된 복구 기법은 정확도를 향상시키면서도 효율성을 유지하는 데 효과적인가?

주요 결과

• 제안된 LODIE 방법은 테스트된 모든 시간 스텝, 특히 큰 Δt = 0.4에서도 조건부 안정성을 확보하며, 실제 단백질 시스템의 비스듬한 해에서도 안정성을 유지한다.
• LODIE 방법은 원자 수에 대해 선형적 확장성을 확보하여 대규모 단백질에서는 조건부 안정성 있는 ADI 방법보다 20배 이상 빠르게 작동한다.
• 모든 제안된 방법은 중심 차분과 1차 시간 분할 기법으로 인해 공간과 시간 모두에서 1차 수렴성을 보이며, 가속화 없이서는 자유 에너지 추정의 정확도가 낮다.
• 리치어드슨 보간법과 복구 기법을 적용한 후, 가속화된 LODIE 방법은 고정된 소수의 시간 스텝만으로도 신뢰할 수 있는 전기적 자유 에너지 추정을 제공한다.
• 전체 복잡도는 여전히 O(N³)이지만, 삼중 대각선 해법의 사용으로 높은 효율성과 대규모 생물분자 시스템에 대한 실용적 타당성을 확보한다.
• 모든 방법 중에서 LODIE 방법은 안정성, 정확도, 효율성의 최적의 균형을 보이며, 실제 단백질의 전기적 자유 에너지 분석에 권장된다.