Finite density QCD phase structure from strangeness fluctuations

Charting the phase diagram of quantum chromodynamics (QCD) at large density is a challenging task due to the complex action problem in lattice simulations. Through simulations at imaginary baryon chemical potential μ B we observe that, if the strangeness neutrality condition is imposed, both the strangeness chemical potential μ S / μ B and the strangeness susceptibility χ 2 S take on constant values at the chiral transition for varying μ B . We present new lattice data to extrapolate contours of constant μ S / μ B or χ 2 S to finite baryon chemical potential. We argue that they are good proxies for the QCD crossover because, as we show, they are only mildly influenced by criticality and by finite volume effects. We obtain continuum limits for these proxies up to μ B = 400 MeV , through a next-to-next-to-leading order ( N 2 LO ) Taylor expansion based on large-statistics data on 16 3 × 8 , 20 3 × 10 and 24 3 × 12 lattices with our 4HEX improved staggered action. We show that these are in excellent agreement with existing results for the chiral transition and, strikingly, also with analogous contours obtained with the hadron resonance gas (HRG) model. On the 16 3 × 8 lattice, we carry out the expansion up to next-to-next-to-next-to-next-to-leading order ( N 4 LO ), and extend the extrapolation beyond μ B = 500 MeV , again finding perfect agreement with the HRG model. This suggests that the crossover line constructed from this proxy starts deviating from the chemical freeze-out line near μ B ≈ 500 MeV , as expected but not yet observed.