[논문 리뷰] Stable and accurate schemes for Langevin dynamics with general kinetic energies

Gabriel Stoltz, Zofia Trstanova|arXiv (Cornell University)|2016. 09. 09.

Markov Chains and Monte Carlo Methods참고 문헌 3인용 수 5

한 줄 요약

이 논문은 비표준 운동 에너지(예: 적응형으로 제약된 또는 전역 리프시츠 조건을 만족하지 않는 형태)를 갖는 랑주뱅 동역학에 대해 안정적이고 약한 두 번째 차수의 수치적 방법을 제안한다. 이는 분할 적분기와 메트로폴리스 보정을 조합하여 정확성과 안정성을 확보한다. 이 방법은 보른츠만-기블즈 분포의 효율적 샘플링을 가능하게 하며, 복잡한 시스템에서의 비정상 상태를 감소시킨다.

ABSTRACT

We study Langevin dynamics with a kinetic energy different from the standard, quadratic one in order to accelerate the sampling of the Boltzmann-Gibbs distribution. We consider two cases: kinetic energies which are local perturbations of the standard kinetic energy around the origin, where they vanish (this corresponds to the so-called adaptively restrained Langevin dynamics); and more general non-globally Lipschitz energies. We develop numerical schemes which are stable and of weak order two, by considering splitting strategies where the discretizations of the fluctuation/dissipation are corrected by a Metropolis procedure. We use the newly developped schemes for two applications: optimizing the shape of the kinetic energy for the adaptively restrained Langevin dynamics, and reducing the metastability of some toy models with non-globally Lipschitz kinetic energies.

연구 동기 및 목표

• 비표준 운동 에너지를 사용할 때 보른츠만-기블즈 분포를 효율적으로 샘플링하는 데 도전하는 것.
• 전역 리프시츠 조건을 만족하지 않거나 표준 제곱형 형태의 국소적 변형을 갖는 운동 에너지에 대해서도 안정적이고 정확한 수치적 방법을 개발하는 것.
• 적응형으로 제약된 랑주뱅 동역학에서 운동 에너지의 형태를 최적화하여 샘플링 성능을 향상시키는 것.
• 비전역 리프시츠 조건을 만족하지 않는 운동 에너지를 갖는 시스템에서 비정상 상태를 줄이기 위해 개선된 수치적 통합 기법을 적용하는 것.

제안 방법

• 역학을 결정론적 및 확률적 성분으로 분리하는 분할 적분기를 설계하며, 변동/소산 항에 대해 보정을 적용한다.
• 각 통합 단계 이후 메트로폴리스-해스팅스 보정을 적용하여 상세 균형을 확보하고 약한 두 번째 차수 정확성을 확보한다.
• 분할 및 보정 단계를 신중히 조율하여 특이점이나 비리프시츠 행동을 보이는 운동 에너지에 대해서도 안정성을 확보하는 방법을 구성한다.
• 제안된 방법을 사용하여 수치 실험을 통해 적응형으로 제약된 랑주뱅 동역학에서 운동 에너지의 함수 형태를 최적화한다.
• 비전역 리프시츠 조건을 만족하지 않는 운동 에너지를 갖는 만능 모델에 방법을 적용하여 비정상 상태 감소 및 더 빠른 수렴을 입증한다.

실험 결과

연구 질문

• RQ1 운동 에너지가 전역 리프시츠 조건을 만족하지 않을 경우, 랑주뱅 동역학의 수치적 방법을 어떻게 안정적이고 약한 두 번째 차수로 유지할 수 있는가?
• RQ2 비표준 운동 에너지 환경에서 분할 적분기의 안정성과 정확성에 메트로폴리스 보정이 미치는 영향은 무엇인가?
• RQ3 적응형으로 제약된 랑주뱅 동역학에서 운동 에너지의 형태를 최적화하여 보른츠만-기블즈 분포의 샘플링 속도를 높일 수 있는가?
• RQ4 비전역 리프시츠 조건을 만족하는 운동 에너지를 갖는 모델에서 제안된 방법이 비정상 상태를 어느 정도 감소시키는가?

주요 결과

• 제안된 메트로폴리스 보정이 적용된 분할 적분기 방법은 약한 두 번째 차수 정확성을 확보하며, 비전역 리프시츠 조건을 만족하지 않는 운동 에너지에 대해서도 안정성을 유지한다.
• 이러한 방법은 병리적인 운동 에너지 함수를 갖는 만능 모델에서 비정상 상태를 성공적으로 감소시켜 평형 상태에 더 빠르게 수렴하도록 한다.
• 적응형으로 제약된 랑주뱅 동역학에서 운동 에너지의 형태를 최적화하면 수치적으로 입증된 바와 같이 샘플링 효율성이 향상된다.
• 이 방법은 목표 보른츠만-기블즈 분포의 불변 측도를 유지하여 장기적인 통계적 행동이 정확하게 유지됨을 보장한다.