2013/07/27 by Xiao-Yong Jin, Yoshinobu Kuramashi, Y. Kuramashi +4 · 11 citations
Mathematics · Physics and Astronomy · #Combinatorics #Condensed matter physics #Exponent #High-Energy Particle Collisions Research #Lattice QCD #Mathematical physics #Mathematics #Order (exchange) #Particle physics #Phase transition #Physics #Quantum Chromodynamics and Particle Interactions #Quantum chromodynamics #Quark #Theoretical and Computational Physics #hep-lat
paper · pdf · doi:10.1103/physrevd.88.094508
published in Physical review. D. Particles, fields, gravitation, and cosmology/Physical review. D. Particles and fields 88(9) (American Physical Society) · 35 pages
arxiv created 2013/07/27 · openalex publication_date 2013/11/22 · arxiv updated 2013/11/27 · openalex created_date 2016/06/24 · openalex updated_date 2026/08/05
We explore the phase space spanned by the temperature and the chemical potential for four-flavor lattice QCD using the Wilson-clover quark action. In order to determine the order of the phase transition, we apply finite-size scaling analyses to gluonic and quark observables, including plaquette, Polyakov loop, and quark number density, and examine their susceptibility, skewness, kurtosis, and Challa-Landau-Binder cumulant. Simulations were carried out on lattices of a temporal size fixed at Nt=4 and spatial sizes chosen from 63 up to 103. Configurations were generated using the phase-reweighting approach, while the value of the phase of the quark determinant was carefully monitored. The \ensuremathμ-parameter reweighting technique is employed to precisely locate the point of the phase transition. Among various approximation schemes for calculating the ratio of quark determinants needed for \ensuremathμ reweighting, we found the Taylor expansion of the logarithm of the quark determinant to be the most reliable. Our finite-size analyses show that the transition is first order at (\ensuremathβ,\ensuremathκ,\ensuremathμ/T)=(1.58,0.1385,0.584\ifmmode±\else\textpm\fi0.008), where (m_\ensuremathπ/m_\ensuremathρ,T/m_\ensuremathρ)=(0.822,0.154). It weakens considerably at (\ensuremathβ,\ensuremathκ,\ensuremathμ/T)=(1.60,0.1371,0.821\ifmmode±\else\textpm\fi0.008), where (m_\ensuremathπ/m_\ensuremathρ,T/m_\ensuremathρ)=(0.839,0.150), and a crossover rather than a first-order phase transition cannot be ruled out.