[논문 리뷰] Numerical methods for the sign problem in Lattice Field Theory

Lorenzo Bongiovanni|arXiv (Cornell University)|2016. 03. 10.

Theoretical and Computational Physics참고 문헌 16인용 수 3

한 줄 요약

이 박사학위논문은 라티스 양자장이론에서의 부호 문제를 극복하기 위한 세 가지 주요 접근법을 조사한다: 복소 랑주뱅 동역학과 게이지 냉각, 레프셰츠 심부, 그리고 LLR 알고리즘을 사용한 상태 밀도 방법. 복소 랑주뱅 동역학이 SU(3) 양-밀스 이론에 θ-항이 있을 경우에도 타당함을 입증하고, 장난감 모형들에 대해 심부와 비교하며, 상대성 양자역학적 보즈 가스에 대해 유한한 화학적 포텐셜에서 상태 밀도 방법을 적용하여 해석적 결과와 양호한 일치를 보인다.

ABSTRACT

The great majority of algorithms employed in the study of lattice field theory are based on Monte Carlo's importance sampling method, i.e. on probability interpretation of the Boltzmann weight. Unfortunately in many theories of interest one cannot associated a real and positive weight to every configuration, that is because their action is explicitly complex or because the weight is multiplied by some non positive term. In this cases one says that the theory on the lattice is affected by the sign problem. An outstanding example of sign problem preventing a quantum field theory to be studied, is QCD at finite chemical potential. Whenever the sign problem is present, standard Monte Carlo methods are problematic to apply and, in general, new approaches are needed to explore the phase diagram of the complex theory. Here we will review three of the main candidate methods to deal with the sign problem, namely complex Langevin dynamics, Lefschetz thimbles and density of states method. We will first study complex Langevin dynamics, combined with the gauge cooling method, on the one-dimensional Polyakov line model, and then we will apply it to pure gauge Yang-Mills theory with a topological theta-term. It follows a comparison between complex Langevin dynamics and the Lefschetz thimbles method on three toy models, which are the quartic model, the U(1) one-link model with a mu dependent determinant, and the SU(2) non abelian one-link model with complex beta parameter. Lastly, we introduce the density of state method, based on the LLR algorithm, and we will employ it in the study of the relativistic Bose gas at finite chemical potential.

연구 동기 및 목표

• 행동이 복소수이거나 가중치가 정규형이 아닐 때 표준 몬테카를로 방법을 적용할 수 없게 되는 라티스 양자장이론에서의 부호 문제를 해결하기 위해.
• 부호 문제가 있는 이론을 다루는 데에 유망한 세 가지 수치적 방법—복소 랑주뱅 동역학, 레프셰츠 심부, 상태 밀도—의 평가 및 비교를 위해.
• 비아벨 게이지 이론, 특히 θ-항이 있는 SU(3) 양-밀스 이론에서 게이지 냉각을 사용한 복소 랑주뱅 동역학의 강건성과 정확도를 시험하기 위해.
• 상대성 양자역학적 보즈 가스에 대해 LLR 알고리즘을 통한 상태 밀도 방법을 적용하여 유한한 화학적 포텐셜에서의 부호 문제가 심각한 시스템에 대해 수행하기 위해.
• 복소 랑주뱅 동역학과 심부를 기준 모형들—4차 모형, U(1) 한점 모형, SU(2) 한점 모형—에 대해 체계적으로 비교하기 위해.

제안 방법

• 복소 랑주뱅 동역학은 장 변수를 복소수로 확장하고, Fokker-Planck 방정식에 따라 확률적으로 진화시키며, 비아벨 이론에서는 게이지 냉각을 사용해 동역학을 안정화시킨다.
• 레프셰츠 심부는 복소화된 다양체에서 행동의 임계점을 찾고, 행동의 허수부가 일정한 다양체를 따라 적분함으로써 실수이고 양수인 측도를 확보한다.
• 상태 밀도 방법은 LLR 알고리즘을 사용해 제약 조건이 있는 시뮬레이션의 연속에서 상태 밀도를 계산하고, 매우 진동하는 적분으로부터 분할함수를 재구성한다.
• 상대성 양자역학적 보즈 가스의 경우, 제어 변수에 대한 분할함수의 로그 도함수를 피팅하여 LLR 알고리즘을 통해 자유 에너지의 차이를 계산한다.
• 복소 랑주뱅 시뮬레이션에서는 기울기 흐름을 통해 위상 전하를 모니터링하여 위상 수반도를 추출하고 수렴성을 검증한다.
• 복소 랑주뱅과 심부 간의 비교는 장 구성의 분포 분석과 각 심부가 경로 적분에 기여하는 방식을 분석함으로써 수행된다.

실험 결과

연구 질문

• RQ1 복소 랑주뱅 동역학에 게이지 냉각을 적용할 경우, 부호 문제가 심각한 SU(3) 양-밀스 이론에 θ-항이 있을 때 신뢰할 수 있는 결과를 도출할 수 있는가?
• RQ2 복소 랑주뱅 동역학의 결과가 4차 모형, μ에 의존하는 행렬식이 있는 U(1) 한점 모형, 그리고 복소수 β를 가진 SU(2) 한점 모형에서 레프셰츠 심부로 얻은 결과와 어떻게 비교되는가?
• RQ3 LLR 알고리즘에 기반한 상태 밀도 방법이 상대성 양자역학적 보즈 가스에서 자유 에너지와 분할함수를 얼마나 정확하게 계산하는가?
• RQ4 게이지 냉각은 비아贝尔 게이지 이론에서 복소 랑주뱅 동역학을 안정화시키는 데 어떤 역할을 하는가? 그리고 위상 전하 분포의 수렴성에 어떤 영향을 미치는가?
• RQ5 다양한 레프셰츠 심부의 기여는 경로 적분에 어떻게 영향을 미치며, 복소 행동이 존재할 경우 어떤 고정점이 측도를 지배하는가?

주요 결과

• 게이지 냉각을 사용한 복소 랑주뱅 동역학은 θ-항이 있는 SU(3) 양-밀스 이론에서 기대되는 위상 전하 분포를 성공적으로 재현하며, 허수 θ일 경우 HMC 결과와 일치한다.
• 실수 θ와 허수 θ 값에 대해 위상 전하 분포는 반대 방향의 곡률을 보이며, 예상되는 CP-홀드 행동을 확인하고 방법의 일관성을 검증한다.
• 4차 모형에서 복소 랑주뱅 동역학은 장 구성의 분포를 올바르게 포착하며, 안정된 심부가 시뮬레이션에서 관측된 주요 기여와 일치한다.
• U(1) 한점 모형에서 이 방법은 올바른 기여하는 심부를 식별하며, μ에 의존하는 행렬식이 존재할 경우 복소 랑주뱅 결과가 심부 기반 예측과 일치한다.
• 복소수 β를 가진 SU(2) 한점 모형에서 게이지 냉각을 사용한 복소 랑주뱅 시뮬레이션은 정확한 분포로 수렴하며, 심부 분석은 단일 안정된 심부의 지배를 확인한다.
• LLR 알고리즘을 사용한 상태 밀도 방법은 상대성 양자역학적 보즈 가스에서 자유 에너지의 차이를 정확하게 계산하며, 통계 오차 범위 내에서 해석적 평균장 이론 예측과 일치한다.