[논문 리뷰] Speeding up HMC with better integrators

M. A. Clark, A.D. Kennedy|ArXiv.org|2007. 10. 18.

Numerical methods for differential equations참고 문헌 7인용 수 10

한 줄 요약

이 논문은 동적 페르미온 계산을 위한 해밀턴 몬테카를로(HMC) 시뮬레이션에서 에너지 보존을 크게 향상시키기 위해 고차수 심플렉틱 적분기—특히 力量 기울기 및 헤시안 기반 적분기—를 사용하는 것을 제안한다. 포아송 괄호를 활용해 영향 Hamiltonian를 유도하고 적분기 파arameter를 최적화하여 단계 크기 의존성을 줄이고 안정성을 유지함으로써, 격자 양성장 이론 시뮬레이션에서 더 낮은 계산 비용과 더 높은 수용률을 달성한다.

ABSTRACT

We discuss how dynamical fermion computations may be made yet cheaper by using symplectic integrators that conserve energy much more accurately without decreasing the integration step size. We first explain why symplectic integrators exactly conserve a ``shadow'' Hamiltonian close to the desired one, and how this Hamiltonian may be computed in terms of Poisson brackets. We then discuss how classical mechanics may be implemented on Lie groups and derive the form of the Poisson brackets and force terms for some interesting integrators such as those making use of second derivatives of the action (Hessian or force gradient integrators). We hope that these will be seen to greatly improve energy conservation for only a small additional cost and that their use will significantly reduce the cost of dynamical fermion computations.

연구 동기 및 목표

• 격자 양성장역학(QCD)에서 동적 페르미온 시뮬레이션의 계산 비용을 줄이기 위해.
• 통합 단계 크기를 줄이지 않고도 HMC에서 에너지 보존을 향상시키기 위해 고차수 심플렉틱 적분기를 사용하기 위해.
• 영향 Hamiltonian의 측정 가능한 포아송 괄호 항을 바탕으로 한 적분기 최적화 프레임워크를 개발하기 위해.
• 격자 양성장 이론에 대한 HMC에서 力量 기울기 및 헤시안 기반 적분기의 성능 향상 잠재력을 탐색하기 위해.
• 파arameter 최적화를 통해 진정한 Hamiltonian와 영향 Hamiltonian 사이의 차이를 최소화하여 더 효율적인 HMC를 가능하게 하기 위해.

제안 방법

• 베이커–캠프벨–하우스도르프(BCH) 전개를 통해 유도된 진짜 Hamiltonian에 가까운 '영향' Hamiltonian를 정확히 보존하는 심플렉틱 적분기를 사용한다.
• 포아송 괄호를 적용하여 영향 Hamiltonian를 계산하고, 고차수 적분기의 경우 O(δτ⁴)까지 보정을 적용한다.
• 행동의 두 번째 도함수를 포함시켜 에너지 보존을 향상시키는 力量 기울기 및 헤시안 기반 적분기를 유도한다.
• 영향 Hamiltonian의 O(δτ²) 및 O(δτ⁴) 항을 최소화함으로써 적분기 파arameter(예: PQPQP 계획에서의 α)를 최적화한다.
• 측정 가능한 포아송 괄호 평균을 사용하여 적분기 파arameter를 사전 최적화하여 δH = H̃ − H를 최소화하는 것을 제안한다.
• 리 군 방법을 활용하여 다양체 위에서 고전 역학을 구현함으로써, 복잡한 적분기의 경우 힘 항과 포아송 괄호의 체계적인 유도를 가능하게 한다.

실험 결과

연구 질문

• RQ1 고차수 심플렉틱 적분기는 단계 크기를 줄이지 않고도 HMC에서 에너지 이탈을 줄일 수 있는가?
• RQ2 BCH 전개에서의 포아송 괄호 구조는 영향 Hamiltonian와 시뮬레이션 안정성에 어떻게 영향을 미치는가?
• RQ3 력 기울기 및 헤시안 기반 적분기는 동적 페르미온 시뮬레이션에서 에너지 보존을 어느 정도 향상시킬 수 있는가?
• RQ4 포아송 괄호 측정을 바탕으로 한 적분기 파arameter의 사전 최적화는 계산 비용을 줄일 수 있는가?
• RQ5 영향 Hamiltonian의 고차수 항은 수용률과 시뮬레이션 효율성에 어떤 영향을 미치는가?

주요 결과

• 최적화된 α = (1 - √3)/6를 가진 PQPQP 적분기는 영향 Hamiltonian의 O(δτ²) 항을 최소화하여 에너지 이탈을 감소시킨다.
• 영향 Hamiltonian에 O(δτ⁴) 보정을 포함하는 力量 기울기 적분기는 추가 비용이 거의 들지 않으면서도 에너지 보존 향상 잠재력이 뚜렷하다.
• 포아송 괄호에서 유도된 계수를 가진 헤시안 기반 적분기는 영향 Hamiltonian 보정을 O(δτ⁴)로 달성하여 표준 적분기 대비 δH를 수개월 정도로 감소시킨다.
• 력 기울기 적분기의 영향 Hamiltonian는 복잡한 중첩된 포아송 괄호의 조합을 포함하며, BCH 전개를 O(δτ⁴)까지 적용하여 유도된 계수를 가진다.
• 형태 e^{(1/6)δτ T}e^{(3/8)δτ S}e^{(1/3)δτ T}...를 가진 적분기는 영향 Hamiltonian 보정 δH ∝ δτ⁴ × (2259{S,{S,{S,{S,T}}}} + ...)를 달성하며, 계수는 BCH 전개에서 유도된다.
• 측정 가능한 포아송 괄호 평균을 직접 측정하여 적분기 파arameter를 최적화함으로써, δH를 감소시키고 HMC 시뮬레이션의 수용률을 향상시켜 성능 향상을 이룬다.