[논문 리뷰] The strong convergence of operator-splitting methods for the Langevin dynamics model

Adam Telatovich, Xiantao Li|arXiv (Cornell University)|2017. 06. 13.

Theoretical and Computational Physics참고 문헌 34인용 수 6

한 줄 요약

이 논문은 추가 노이즈가 있는 랭지비안 역학 모델에서 연산자 분할 방법에 대한 엄밀한 강수렴 분석을 제시한다. 키투나의 해 표현식과 이토-테일러 전개를 활용하여, 강수렴 순서가 최대 3에 이르는 새로운 고차수 연산자 분할 스킴을 유도하였으며, 표준 분할 방법인 SVV와 난이도 없는 분할 방법이 오직 순서 1에만 도달한다는 것을 입증한다.

ABSTRACT

We study the strong convergence of some operator-splitting methods for the Langevin dynamics model with additive noise. It will be shown that a direct splitting of deterministic and random terms, including the symmetric splitting methods, only offers strong convergence of order 1. To improve the order of strong convergence, a new class of operator-splitting methods based on Kunita's solution representation are proposed. We present stochastic algorithms with strong orders up to 3. Both mathematical analysis and numerical evidence are provided to verify the desired order of accuracy.

연구 동기 및 목표

• 추가 노이즈가 있는 랭지비안 역학 모델에 널리 사용되는 연산자 분할 방법의 강수렴 행동을 분석하기 위해.
• 기존 분할 방법의 근본적 한계를 규명하기 위해 — 예를 들어, 스토케스틱 비어티스 Verlet(SVV) 및 난이도 없는 분할 방법은 오직 강수렴 순서 1에만 도달한다는 것을 보여주기 위해.
• 스토크스틱 미분 방정식에 대한 키투나의 해 표현식을 기반으로 한 새로운 고차수 연산자 분할 스킴의 클래스를 개발하기 위해.
• 수학적 분석과 수치 실험을 통해 이론적 수렴 순서를 검증하기 위해.
• 제안된 고차수 스킴의 과다진동 및 과다감쇠 극한에서의 거동을 검토하기 위해.

제안 방법

• 저자들은 이토 SDE에 대한 키투나의 해 표현식을 사용하여, 해를 드리프트 및 확산 미분 연산자의 교환자들을 포함하는 연산자 지수로 표현한다.
• 키투나 전개를 점진적으로 단계별로 잘라내어 고차수 수치 스킴을 구성하며, 계수들은 드리프트 및 확산 항에 대해 반복적으로 미분 연산자를 적용하여 유도된다.
• 다변수 이토 적분을 최대 차수 3까지 계산하며, 예를 들어 $ \bigint_{0}^{\tau} \bigint_{0}^{s_3} \bigint_{0}^{s_2} dW_{s_1} ds_2 ds_3 $ 와 같은 항들은 정확도를 확보하기 위해 해석적으로 평가된다.
• 스킴은 다중지표 $ \boldsymbol{\nu} $ 에 따라 항들을 그룹화하며, 총 수준 $ l(\boldsymbol{\nu}) + n(\boldsymbol{\nu}) $ 에 의해 순서가 결정된다. 여기서 $ l $ 은 위너 적분의 수를 세며, $ n $ 은 미분 연산자의 적용 횟수를 의미한다.
• 추가 노이즈의 경우, 강수렴 순서 3과 순서 2.5 스킴이 동일하다. 이는 모든 추가 계수 함수 $ g_{\boldsymbol{\nu}} $ 가 $ \boldsymbol{\nu} \notin \tilde{\boldsymbol{\nu}}_{2.5} $ 일 때, 상수 또는 0 도함수이기 때문에 사라지기 때문이다.
• 결과로 도출된 스킴은 명시적이고 적분 가능한 부분단계로 구현되며, 기존 분할 방법의 구현 용이성은 유지하면서도 더 높은 강수렴 순서를 달성한다.

실험 결과

연구 질문

• RQ1 왜 SVV 및 난이도 없는 분할과 같은 표준 연산자 분할 방법은 분자의 역학 시뮬레이션에서 널리 사용되지만도 강수렴 순서 1에만 도달하는가?
• RQ2 추가 노이즈가 있는 랭지비안 역학 모델에 대해 연산자 분할 스킴에서 고차수 강수렴을 달성할 수 있는가?
• RQ3 키투나의 해 표현식과 이토-테일러 전개를 사용한 고차수 스킴의 이론적 기반은 무엇인가?
• RQ4 강수렴 오차 노름을 기준으로 새로운 스킴이 고전적 방법에 비해 정확도에서 어떻게 비교되는가?
• RQ5 새로운 고차수 스킴은 과다진동 및 과다감쇠 극한에서 정확한 거동을 유지하는가?

주요 결과

• 스토크스틱 비어티스 버전(SVV) 및 난이도 없는 분할 방법은 부분단계에서 정확함에도 불구하고, 강수렴 순서 1에만 도달한다.
• 키투나의 해 표현식을 적용하고 이토-테일러 전개를 고차수 수준에서 잘라내어, 강수렴 순서가 최대 3에 이르는 새로운 연산자 분할 스킴을 유도하였다.
• 추가 노이즈의 경우, 강수렴 순서 3과 순서 2.5 스킴은 동일하다. 이는 모든 추가 계수 함수 $ g_{\boldsymbol{\nu}} $ 가 $ \boldsymbol{\nu} \notin \tilde{\boldsymbol{\nu}}_{2.5} $ 일 때, $ L^1 $ 및 $ L^0 $ 작용에 의해 $ g_{(1,0,0)} $ 와 $ g_{(1,0,0,0)} $ 가 상수이기 때문에 사라지기 때문이다.
• 새로운 스킴은 수학적으로 필요한 강수렴 순서를 달성함을 증명하였으며, 수치적 증거는 이론적 예측을 확인한다.
• 이 방법들은 과다진동 및 과다감쇠 극한 모두에 적용 가능하며, 두 영역 모두에서 일관된 수렴 행동을 보였다.
• 이 방법은 이토-테일러 방법에서 요구하는 力의 고차수 도함수 계산이 필요로 하지 않으며, 룬게-쿠타 방법과는 달리 다수의 힘 평가도 필요로 하지 않아 분자의 역학에 실용적이다.