QCD phase diagram for large Nf : analysis from contact interaction effective potential

In this paper, we discuss the impact of a higher number of light quark flavors, Nf , on the QCD phase diagram under extreme conditions. Our formalism is based on the Schwinger-Dyson equation, employing a specific symmetry-preserving vector-vector flavor-dressed contact interaction model of quarks in Landau gauge, utilizing the rainbow-Ladder truncation. We derive expressions for the dressed quark mass Mf and effective potential Ω f at zero, at finite temperature T and the quark chemical potential µ. The transition between chiral symmetry breaking and restoration is triggered by the effective potential of the contact interaction, whereas the confinement and deconfinement transition is approximated from the confinement length scale τ˜ir. Our analysis reveals that at (T = µ = 0), increasing Nf leads to the restoration of chiral symmetry and the deconfinement of quarks when Nf reaches its critical value, N c f ≈ 8. At this critical value, in the chiral limit (mf = 0), the global minimum of the effective potential occurs at the point where the dressed quark mass approaches zero (Mf → 0). However, when a bare quark mass of mf = 7 MeV is introduced, the global minimum shifts slightly to a nonzero value, approaching Mf → mf . At finite T and µ, we illustrate the QCD phase diagram in the (T χ,C c − µ) plane, for various numbers of light quark flavors, noting that both the critical temperature Tc and the critical chemical potential µc for chiral symmetry restoration and deconfinement decrease as Nf increases. Moreover, the critical endpoint (TEP , µEP ) also shifts to lower values with increasing Nf . Our findings are consistent with other low-energy QCD approaches.