[논문 리뷰] Quark Mass Dependence of the QCD Critical End Point in the Strong Coupling Limit

Jangho Kim, Wolfgang Unger|arXiv (Cornell University)|2016. 11. 28.

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

한 줄 요약

이 논문은 강한 결합 근사에서 라티스 QCD의 이중 표현을 사용하여 몽테카를로 시뮬레이션을 통해 강력 상호작용의 임계 끝점(Critical End Point, CEP)의 쿼크 질량 의존성을 조사한다. 서명 문제를 다루기 위해 웜 알고리즘을 적용하여 $μ$-$T$ 상도를 그려내고 유한한 쿼크 질량에서 CEP의 위치를 결정하며, 이는 CEP가 $m_q$에 따라 이동하고 비임계 이방성 $\xi(\gamma)$가 상계 경계에서 물리적으로 비합리적인 뒤로 굽힘 현상을 제거함을 보여준다.

ABSTRACT

Strong coupling lattice QCD in the dual representation allows to study the full $μ$-$T$ phase diagram, due to the mildness of the finite density sign problem. Such simulations have been performed in the chiral limit, both at finite $N_t$ and in the continuous time limit. Here we extend the phase diagram to finite quark masses, with an emphasis on the low temperature first order transition. We present our results on the quark mass dependence of the critical end point and the first order line obtained by Monte Carlo via the worm algorithm.

연구 동기 및 목표

• 색 대칭 근사 이론을 넘어서 유한한 쿼크 질량에서 QCD의 $μ$-$T$ 상도를 확장하기 위해.
• 강한 결합 근사에서 다양한 쿼크 질량에 대해 임계 끝점(Critical End Point, CEP)의 위치를 결정하기 위해.
• 비임계 이방성 $\xi(\gamma)$를 사용하여 페르미온 이론에서의 $\xi' = \gamma^2$와 비교할 때 관측된 물리적으로 비합리적인 상계 경계의 뒤로 굽힘 현상을 제거하기 위해.
• 이전 몽테카를로 시뮬레이션과 평균장 이론 결과와의 비교를 통해 방법의 정확성을 검증하기 위해.

제안 방법

• 연구는 이중 표현에서의 강한 결합 라티스 QCD를 사용하며, 게이지 및 그라스만 변수를 통합하여 메손과 바리온의 월드라인 표현을 도출한다.
• 분할 함수는 유한한 쿼크 질량을 위한 메손 전이($k_b$), 단량체($n_x$), 바리온 고리($\ell$)와 그 래핑 수 $r_\ell$ 및 부호 $\sigma(\ell)$로 표현된다.
• 관측량인 쿼크-항성기대값 $\langle\bar{\chi}\chi\rangle$, 쿼크-항성기대값의 임계도수 $\chi_{ch}$, 바리온 밀도 $\langle n_B \rangle$, 바리온 임계도수 $\chi_B$ 등을 계산하기 위해 몽테카를로 샘플링에 웜 알고리즘을 사용한다.
• 크기 효과 보정을 위해 $O(2)$, $Z(2)$, 그리고 전이($\rho=1$) 임계 지수를 사용하여 전이를 식별하고 CEP의 위치를 특정한다.
• 시간적 및 공간적 척도를 정확히 매핑하기 위해 비임계적으로 유도된 이방성 매개변수 $\xi(\gamma)$를 사용하며, 물리적으로 비합리적인 $\xi' = \gamma^2$를 대체한다.
• 임계도수 및 기대값 데이터의 피크 분석을 통해 CEP의 위치를 특정하며, 다양한 $m_q$에 대해 $aT_E'$ 및 $a\mu_E'$를 추출한다.

실험 결과

연구 질문

• RQ1 강한 결합 근사에서 쿼크 질량이 증가함에 따라 $μ$-$T$ 평면의 임계 끝점(Critical End Point, CEP)은 어떻게 이동하는가?
• RQ2 비임계 이방성 $\xi(\gamma)$를 사용할 경우, $\xi' = \gamma^2$를 사용했을 때 관측된 제1종 전이 경계의 물리적으로 비합리적인 뒤로 굽힘 현상은 제거되는가?
• RQ3 이 몽테카를로 연구에서 도출된 CEP 위치는 이전 몽테카를로 결과와 평균장 이론 예측과 정량적으로 어떻게 비교되는가?
• RQ4 유한한 $m_q$에서 단량체의 변동은 쿼크-항성기대값의 임계도수에 어떤 역할을 하는가? 그리고 CEP 결정의 정밀도에 어떤 영향을 미치는가?
• RQ5 정확한 비임계 이방성 매개변수를 사용할 경우, CEP 궤적에 대해 $(am_q)^{2/5}$에서의 삼중임계 스케일링 행동이 유지되는가?

주요 결과

• 임계 끝점(Critical End Point, CEP)은 쿼크 질량에 따라 이동한다: $am_q = 0.05$일 때, CEP는 $aT_E' = 0.73(3)$ 및 $a\mu_E' = 0.7808(5)$에 위치하며, 물리적 단위에서는 $aT_E = 0.62(3)$ 및 $a\mu_E = 0.6653(4)$이다.
• 비임계 이방성 $\xi(\gamma)$를 사용함으로써, $\xi' = \gamma^2$를 사용했을 때 관측된 저온 영역에서의 물리적으로 비합리적인 상계 경계의 뒤로 굽힘 현상이 제거된다.
• 유한한 $m_q$에 대한 CEP 궤적은 평균장 이론과 일치하지 않으며, CEP 위치와 $(am_q)^{2/5}$에서의 기울기 모두 평균장 예측과 다름을 보인다.
• $am_q = 0.05$일 때, 쿼크-항성기대값의 임계도수 피크는 $aT' = 0.725$와 $0.75$ 사이에 위치하며, 이는 CEP가 이 범위 내에 있음을 시사하며 $Z(2)$ 임계 지수에 의해 확인된다.
• $am_q = 0.01$ 및 $0.02$에 대한 결과는 이전 몽테카를로 데이터와 일치하여, 방법의 정확성과 수치적 안정성을 검증한다.
• $aT' = 0.7$에서 $\rho = 3$를 사용한 크기 효과 보정은 전이의 제1종 성격을 확인하며, $aT' = 0.95$에서의 전이 스케일링에는 $\rho = 1$이 사용된다.