[논문 리뷰] Electric polarizability of hadrons with overlap fermions on multi-GPUs

Michael Lujan, Andrei Alexandru외 1명|arXiv (Cornell University)|2011. 11. 27.

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

한 줄 요약

이 논문은 양자 색역학 격자에서 강입자, 특히 중성 π 보존의 전기 분극율을 계산하기 위해 고성능 GPU 구현을 제시한다. 다중 GPU 아키텍처와 변동형 공액 기울기 해법기(Deflation)를 활용하여, 각 GPU에서 약 30개의 CPU 코어에 해당하는 성능를 달성하였으며, 이는 카이랄 대칭성이 중요한 저질량 π 보존에서의 시뮬레이션을 가능하게 한다. 초보적 결과는 중성 π 보존의 분극율이 약 500 MeV 근처에서 부호를 바꾸는 것으로 나타났으며, 이는 이론적 예측과 일치한다.

ABSTRACT

Electric polarizability is an important parameter for the internal structure of hadrons. Previous studies of polarizabilities have been done at relatively heavy pion masses, leaving the chiral region largely unexplored. In this report, we use overlap fermions which are known to be computationally demanding to properly capture the chiral dynamics. We present an implementation strategy to construct overlap on multi-GPUs. We find that our GPU code has an equivalent of \sim30 CPU cores to 1 GPU. We also present preliminary results for the polarizability of the neutral pion.

연구 동기 및 목표

• 카이랄 대칭성을 유지하고 특이 구성을 피하기 위해 오버랩 페르미온을 사용하여 강입자의 전기 분극율을 계산하기 위해.
• 오버랩 페르미온의 계산 블로킹 문제를 다중 GPU 아키텍처에 효율적으로 구현하여 해결하기 위해.
• 분극율이 비정상적인 행동을 보일 것으로 예상되는 카이랄 영역에서의 QCD를 탐색하기 위해 저질량 π 보존(약 240 MeV까지)에서 시뮬레이션을 수행하기 위해.
• 대규모 격자 QCD 시뮬레이션을 위해 적응형 다중질량 CG 해법기와 고유계수계 탈출 기법을 결합한 확장 가능하고 메모리 효율적인 GPU 전략을 개발하기 위해.

제안 방법

• 행렬 부호 함수의 Zolotarev 유리 근사법을 사용하여 오버랩 페르미온 작용을 구현하였으며, 다항식 근사를 소·대특이값 하위공간으로 분할하였다.
• 낮은 고유모드를 계산하기 위해 암시적 재시작 아르노르드 분해를 사용하는 고유계수 해법기를 적용하여, 수렴을 가속화하는 탈출 기법을 가능하게 하였다.
• 다중질량, 적응형 공액 기울기(CG) 방법을 사용하여 오버랩 디랙 연산자를 역행렬화하였으며, 성능 최적화를 위해 다항식 차수와 정밀도(단일 정밀도 산술 포함)를 동적으로 조정하였다.
• 배경장 방법을 적용하기 위해, 게이지 링크에 단위 복소수 인자(Phase factors)를 도입하여 일정한 전기장을 도입하였으며, 이로 인해 공변 미분을 $ D_\mu = \partial_\mu - iG_\mu - iqA_\mu $ 로 수정하였다.
• 소스 위치를 다양화하여 통계적 정밀도를 향상시키기 위해, 각 구성에서 5개의 쿼크 전파함수(별도의 $ u $ 및 $ d $ 전파함수 포함)를 사용하여 전기장 유무에 따른 상관 함수를 계산하였다.
• 메모리 제약을 해결하기 위해 데이터를 다수의 GPU에 분산시키고, GPU 메모리가 부족할 경우 하이브리드 CPU-GPU 메모리 전략을 사용하였다.

실험 결과

연구 질문

• RQ1 특히 카이랄 한계 근처에서, 저질량 π 보존의 전기 분극율은 얼마인가?
• RQ2 물리적 π 보존 질량에서의 시뮬레이션을 가능하게 하기 위해, 계산 비용이 큰 오버랩 페르미온 형식을 다중 GPU 아키텍처에 어떻게 효율적으로 구현할 수 있는가?
• RQ3 카이랄 효과 이론에 의해 예측된 것처럼, $ m_\pi \approx 500 \, \text{MeV} $ 영역에서 중성 π 보존의 분극율의 부호가 바뀌는가?
• RQ4 탈출 기법과 적응형 정밀도 CG 방법이 GPU에서 오버랩 페르미온 전파함수의 계산 비용을 어느 정도 감소시킬 수 있는가?

주요 결과

• GPU 구현은 각 GPU당 약 30개의 CPU 코어에 해당하는 성능를 달성하여, 오버랩 페르미온 계산을 크게 가속화하였다.
• 적응형 다중질량 CG와 동적 정밀도 조정을 사용함으로써, $ m_\pi \approx 200 \, \text{MeV} $ 에서 단일 전파함수 계산 시간을 비적응형일 경우 1.6시간에서 0.8시간으로 단축시켰다.
• 중성 π 보존의 분극율은 $ m_\pi \approx 500 \, \text{MeV} $ 에서 부호를 바꾸며, 이는 이론적 예측과 이전 연구 결과와 일치한다.
• 통계적 자료 수가 제한되어 오차 구간이 크므로, 저질량 영역에서 신호 대 잡음 비율을 향상시키기 위해 4배 더 많은 구성 수를 생성하고 있다.
• 16개의 GPU에 고유계수계 및 전파함수 계산을 분산하여 메모리 제약을 성공적으로 처리하였으며, 必要한 경우 하이브리드 CPU-GPU 메모리 사용 전략을 적용하였다.
• 이 방법은 π 보존 질량이 300 MeV 이하인 조건에서 오버랩 페르미온을 사용한 중성 π 보존 분극율의 첫 번째 격자 QCD 계산을 가능하게 하였다. 물리적 점에 가까워지고 있다.