[논문 리뷰] Frequency-splitting estimators for single-propagator traces

Leonardo Giusti, T. H. Harris외 2명|arXiv (Cornell University)|2020. 01. 01.

Particle physics theoretical and experimental studies인용 수 1

한 줄 요약

이 논문은 라티스 QCD에서 단일 전파기 추적 계산의 통계적 노이즈를 줄이기 위해 주파수 분할 추정기법을 도입한다. 이 방법은 큰 쿼크 질량에서 히opping 전개를 사용하면서도 분할-짝수 기법을 결합하여, 복소도수 변동과 효율적인 스토하스틱 추정을 가능하게 하여, 특히 벡터 전류에서 매우 효과적으로 분산을 1–2계단 감소시킨다. 이는 강한 상호작용 진공 분극 및 뮤온 g−2 계산의 계산 비용을 크게 낮춘다.

ABSTRACT

In these proceedings we address the computation of quark-line disconnected diagrams in lattice QCD. The evaluation of these diagrams is required for many phenomenologically interesting observables, but suffers from large statistical errors due to the vacuum and random-noise contributions to their variances. Motivated by a theoretical analysis of the variances, we introduce a new family of stochastic estimators of single-propagator traces built upon a frequency splitting combined with a hopping expansion of the quark propagator, and test their efficiency in two-flavour QCD with pions as light as 190 MeV. The use of these estimators reduces the cost of the computation by one to two orders of magnitude over standard estimators depending on the fermion bilinear. As a concrete application, we show the impact of these findings on the computation of the hadronic vacuum polarization contribution to the muon anomalous magnetic moment.

연구 동기 및 목표

• 무작위 노이즈와 진공 기여로 인한 쿼크 라인 분리 도형에서의 큰 통계적 분산을 줄이기 위해.
• 기본 몬테카를로 시뮬레이션과 호환되면서도 노이즈를 억제하는 단일 전파기 추적을 위한 새로운 스토하스틱 추정기류를 개발하기 위해.
• 뮤온 비례자기모멘트의 강한 상호작용 진공 분극 기여를 계산하는 효율성을 향상시키기 위해.
• 전자기 전류 상관함수 및 이스표인버션 붕괴 보정과 같은 관측량에서 분리 기여를 비용 효율적으로 평가할 수 있도록 하기 위해.

제안 방법

• 다른 쿼크 질량에서의 전파기 간 차이를 위한 분할-짝수 기법을 사용하여 쿼크 전파기 추적을 분해하는 주파수 분할 추정기법을 제안한다.
• 일련의 가장 큰 쿼크 질량에서 히핑 전개(최대 차수 2n)를 적용하여 추적을 계산 가능한 부분(M2n,mr)과 나머지(R)로 분리하고, 나머지는 스토하스틱 추정을 위해 사용한다.
• 연속된 질량 간의 K−1 차이에 대해 분할-짝수 추정기법을 사용하며, 고질량에서 히핑 전개를 적용하여 고주파수 변동을 제어한다.
• 전체 주파수 분할 추정기법을 다음과 같은 합으로 구성한다: (1) 질량 차이에 대한 분할-짝수 추정기, (2) M2n,mr 부분의 정확한 계산, (3) 나머지 R에 대한 스토하스틱 추정기.
• 단위 분산을 가진 독립적인 가우시안 노이즈 장을 사용하여 나머지 항을 추정한다: ̄τRΓ,r(x0) = −1/(aL3Ns) Σi,x Re[aΓ (η†i Hn mr)(x)Γ (D−1mr Hn mr ηi)(x)].
• 최종 추정기법은 ̄τfsΓ,r1(x0) = Σk=1 to K−1 ̄τΓ,rkrk+1(x0) + ̄tMΓ,rK(x0) + ̄τRΓ,rK(x0)로 정의되며, 모든 구성 요소를 통합한다.

실험 결과

연구 질문

• RQ1 주파수(질량) 척도에 따라 기여를 분리함으로써 단일 전파기 추적의 분산을 줄일 수 있는가?
• RQ2 큰 쿼크 질량에서 히핑 전개와 분할-짝수 추정기법을 조합하여 무작위 노이즈를 얼마나 효과적으로 억제할 수 있는가?
• RQ3 물리적 점 근처의 벡터 전류 채널에서 이 방법이 계산 비용을 얼마나 줄일 수 있는가?
• RQ4 이 방법을 뮤온 비례자기모멘트의 분리 기여를 더 효율적으로 계산하는 데 적용할 수 있는가?

주요 결과

• 주파수 분할 추정기법은 F7 및 G8 엔셈블에서 무작위 노이즈 기여의 분산을 약 2계단 감소시켰다.
• 벡터 전류의 경우, 반복당 비용이 3.3–6배 증가했음에도 불구하고, 질량에 따라 10–30배의 계산 비용 절감을 달성하였다.
• 게이지 노이즈 기여는 큰 Ns에서 포화되며, 이는 작은 Ns에서 무작위 노이즈가 지배적임을 확인하고, 새로운 방법이 이를 효과적으로 억제함을 시사한다.
• 특히 히핑 전개를 고질량에서 적용할 경우, 표준 추정기법보다 엄격히 작은 분산을 달성하는 경우가 있었다.
• 수치적 테스트 결과, 1.0 fm까지 분리 전류 상관함수에서 양호한 신호를 확보하였고, 분할-짝수 추정기법을 사용해 aμ 기여를 효율적으로 계산할 수 있었다.
• 이 기법은 저모드 평균화 및 확산과 같은 다른 분산 감소 기법과도 호환되어, 향후 응용에서 추가적인 비용 절감이 가능하다.