Chiral condensate and the equation of state at nonzero baryon density from the hadron resonance gas model with a repulsive mean field

We study the QCD equation of state and the chiral condensate using the hadron resonance gas model with repulsive mean-field interactions. We find that the repulsive interactions improve the agreement with the lattice results on the derivatives of the pressure with respect to the baryon chemical potential up to eighth order. From the temperature dependence of the chiral condensate we estimate the crossover temperature as a function of baryon chemical potential, <a:math xmlns:a="http://www.w3.org/1998/Math/MathML"><a:mrow><a:msub><a:mi>T</a:mi><a:mrow><a:mi>p</a:mi><a:mi>c</a:mi></a:mrow></a:msub><a:mrow><a:mo>(</a:mo><a:msub><a:mi>μ</a:mi><a:mi>B</a:mi></a:msub><a:mo>)</a:mo></a:mrow></a:mrow></a:math>. We find that the chiral crossover line starts to deviate significantly from the chemical freeze-out line already for <b:math xmlns:b="http://www.w3.org/1998/Math/MathML"><b:mrow><b:msub><b:mi>μ</b:mi><b:mi>B</b:mi></b:msub><b:mo>></b:mo><b:mn>400</b:mn><b:mspace width="4pt"/><b:mi>MeV</b:mi></b:mrow></b:math>. Furthermore, we find that the chiral pseudocritical line can be parametrized as <d:math xmlns:d="http://www.w3.org/1998/Math/MathML"><d:mrow><d:msub><d:mi>T</d:mi><d:mrow><d:mi>p</d:mi><d:mi>c</d:mi></d:mrow></d:msub><d:mrow><d:mo>(</d:mo><d:msub><d:mi>μ</d:mi><d:mi>B</d:mi></d:msub><d:mo>)</d:mo></d:mrow><d:mo>/</d:mo><d:msub><d:mi>T</d:mi><d:mrow><d:mi>p</d:mi><d:mi>c</d:mi></d:mrow></d:msub><d:mrow><d:mo>(</d:mo><d:mn>0</d:mn><d:mo>)</d:mo></d:mrow><d:mo>=</d:mo><d:mn>1</d:mn><d:mo>−</d:mo><d:msub><d:mi>κ</d:mi><d:mn>2</d:mn></d:msub><d:msup><d:mrow><d:mo>[</d:mo><d:msub><d:mi>μ</d:mi><d:mi>B</d:mi></d:msub><d:mo>/</d:mo><d:msub><d:mi>T</d:mi><d:mrow><d:mi>p</d:mi><d:mi>c</d:mi></d:mrow></d:msub><d:mrow><d:mo>(</d:mo><d:mn>0</d:mn><d:mo>)</d:mo></d:mrow><d:mo>]</d:mo></d:mrow><d:mn>2</d:mn></d:msup><d:mo>−</d:mo><d:msub><d:mi>κ</d:mi><d:mn>4</d:mn></d:msub><d:msup><d:mrow><d:mo>[</d:mo><d:msub><d:mi>μ</d:mi><d:mi>B</d:mi></d:msub><d:mo>/</d:mo><d:msub><d:mi>T</d:mi><d:mrow><d:mi>p</d:mi><d:mi>c</d:mi></d:mrow></d:msub><d:mrow><d:mo>(</d:mo><d:mn>0</d:mn><d:mo>)</d:mo></d:mrow><d:mo>]</d:mo></d:mrow><d:mn>4</d:mn></d:msup></d:mrow></d:math> with <e:math xmlns:e="http://www.w3.org/1998/Math/MathML"><e:mrow><e:msub><e:mi>κ</e:mi><e:mn>2</e:mn></e:msub><e:mo>=</e:mo><e:mn>0.0150</e:mn><e:mrow><e:mo>(</e:mo><e:mn>2</e:mn><e:mo>)</e:mo></e:mrow></e:mrow></e:math> and <f:math xmlns:f="http://www.w3.org/1998/Math/MathML"><f:mrow><f:msub><f:mi>κ</f:mi><f:mn>4</f:mn></f:msub><f:mo>=</f:mo><f:mn>3.1</f:mn><f:mrow><f:mo>(</f:mo><f:mn>6</f:mn><f:mo>)</f:mo></f:mrow><f:mo>×</f:mo><f:msup><f:mn>10</f:mn><f:mrow><f:mo>−</f:mo><f:mn>5</f:mn></f:mrow></f:msup></f:mrow></f:math>, which are in agreement with lattice QCD results for small values of <g:math xmlns:g="http://www.w3.org/1998/Math/MathML"><g:msub><g:mi>μ</g:mi><g:mi>B</g:mi></g:msub></g:math>. For the first time we find a tiny but nonzero value of <h:math xmlns:h="http://www.w3.org/1998/Math/MathML"><h:msub><h:mi>κ</h:mi><h:mn>4</h:mn></h:msub></h:math> in our study. Published by the American Physical Society 2024