[논문 리뷰] Improving DWF Simulations: the Force Gradient Integrator and the Möbius Accelerated DWF Solver

Hantao Yin, Robert D. Mawhinney|arXiv (Cornell University)|2011. 11. 21.

Quantum Chromodynamics and Particle Interactions참고 문헌 5인용 수 5

한 줄 요약

이 논문은 라티스 QCD 시뮬레이션의 효율성을 향상시키기 위해 힘 기울기 적분기(FGI)와 모비우스 가속 도메인 월 페르미온(DWF) 해법기를 제안한다. FGI는 통합 오차를 감소시키고 더 큰 스텝 크기를 허용하며, 모비우스 해법기는 더 작은 다섯 번째 차원을 사용해 디랙 연산자 적용 횟수를 약 2배 감소시켜 32³×64 격자와 180 MeV의 파이온 질량을 가진 큰 격자에서 2배의 속도 향상을 달성한다.

ABSTRACT

We have implemented a variant of the force gradient integrator proposed by Kennedy et.al. and are using it in our production 2+1 flavor DWF simulations with pion masses of 180 MeV in (4.5fm)3 volumes. We find modest speed-ups (\sim 20%) from using the force gradient integrator, compared to our previously used Omelyan integrator. On other ensembles, primarily finite temperature 2+1 flavor DWF QCD, we have extensively tuned the Hasenbusch preconditioning masses and achieved speed-ups of 2-3x. Here we have also switched to the force gradient integrator, but this change has not had any impact on the speed. We also report on an improved solver for DWF, which uses Möbius fermions, with a smaller fifth dimension than the original DWF fermions, as an intermediate step in the generation of solutions of the Dirac equation. This approach cuts the number of effective Dirac applications by approximately a factor of 2 when the conjugate gradient iteration count is large.

연구 동기 및 목표

• 분자역학 통합 기법을 최적화하여 2+1 편성 도메인 월 페르미온(DWF) 라티스 QCD 시뮬레이션의 효율성을 향상시키는 것.
• HMC 알고리즘에서 시뮬레이션 시간의 대부분을 차지하는 DWF 디랙 연산자 해법의 계산 비용을 감소시키는 것.
• 모비우스 페르미온과 표준 DWF 간의 수학적 유사성을 활용하여 DWF 디랙 연산자 해법을 가속화하는 것.
• 안정성과 스텝 크기를 향상시키면서도 심플렉틱성과 가역성을 유지하는 새로운 힘 기울기 적분기(FGI)를 구현하고 테스트하는 것.
• FGI와 하제누스비히 질량 분할 및 몫 작용을 결합하여 생산 규모의 HMC 시뮬레이션에서 상당한 속도 향상을 달성하는 것.

제안 방법

• 대칭적이고 심플렉틱적인 힘 기울기 적분기(FGI)를 구현하기 위해 두 단계의 업데이트를 사용: 먼저 힘 기울기의 작은 지수를 사용해 게이지 장을 변형한 후, 업데이트된 장을 바탕으로 다시 힘을 재계산한다.
• 2차 힘 기울기 항을 타일러 전개 근사로 계산하여, $ U' = e^{-\frac{\tau^2}{24}F^j T_j}U $를 통해 게이지 장 이동으로 대체함으로써 효율적인 구현을 가능하게 한다.
• 첫 번째 단계의 해를 두 번째 단계의 초기 추측으로 재사용함으로써 '힘 기울기 예측'을 적용하여, 힘 평가 시 디랙 해법 비용을 감소시킨다.
• 더 작은 다섯 번째 차원 $L'$을 사용하는 모비우스 페르미온 작용을 활용해 전체 DWF 해법을 근사하는 모비우스 가속 DWF 해법기를 도입하여 4차원 윌슨 디랙 연산자 적용 횟수를 감소시킨다.
• 잔여 오차를 모비우스 연산자를 사용해 반복적으로 해법함으로써 정확도를 향상시키는 결함 보정 방법을 사용한다.
• HMC 비용을 줄이기 위해 FGI를 하제누스비히 질량 분할 및 몫 작용과 결합하여, 힘과 관련된 디랙 방정식의 해법 정확도를 낮추는 데 사용한다.

실험 결과

연구 질문

• RQ1 대규모 DWF 시뮬레이션에서 오멜리안 적분기 대비 힘 기울기 적분기(FGI)가 측정 가능한 성능 향상을 제공하는가?
• RQ2 모비우스 페르미온 형식은 효과적인 다섯 번째 차원 크기를 줄여 DWF 디랙 방정식의 해법 속도를 높일 수 있는가?
• RQ3 DWF 해법기에서 정확도와 계산 비용을 균형 잡기 위해 모비우스 근사의 최적 파라미터 조합은 무엇인가?
• RQ4 FGI와 하제누스비히 질량 분할 및 몫 작용의 조합이 생산 수준의 HMC 시뮬레이션에서 상당한 속도 향상을 이끌 수 있는가?
• RQ5 힘 기울기 예측 기법은 가역성과 심플렉틱성을 해치지 않으면서도 힘 평가 비용을 어떻게 감소시키는가?

주요 결과

• 4.5\textrm{fm}^3$ 격자에서 180 MeV의 파이온을 가진 경우 FGI는 약 20%의 소량의 속도 향상을 달성했으며, 더 큰 격자에서는 더 큰 성과가 기대된다.
• $32^3 \times 64$ 격자에서 모비우스 가속 DWF 해법기는 4차원 윌슨 디랙 연산자 적용 횟수를 3.6×10⁵에서 2.0×10⁵로 줄여 약 1.8배 감소시켰다.
• 모비우스 가속 DWF 해법기는 직접 해법 대비 벽 시계 시간에서 2배의 속도 향상을 달성했으며, 이는 해법기가 총 시간의 42%만을 차지했다.
• 결함 보정 방법은 모비우스 근사의 정확도를 향상시켜 최소한의 추가 비용으로 고정밀도 해법을 가능하게 하였다.
• 하제누스비히 질량 분할과 몫 작용을 통해 힘과 관련된 디랙 방정식의 해법 정확도를 낮추어 유한 온도 계열에서 2–3배의 속도 향상을 기여하였다.
• 모비우스 파라미터는 한 번의 조정만으로도 앙세블리 전체에 걸쳐 구성과 소스 벡터에 관계없이 효과를 유지하여 생산 시뮬레이션에서 널리 재사용 가능하였다.