Lattice quantum chromodynamics at large isospin density

INSPIRE: 2682305
arXiv: 2307.15014
DOI: 10.1103/PhysRevD.108.114506

Authors: Abbott, Ryan, Detmold, William, Romero-López, Fernando, Davoudi, Zohreh, Illa, Marc, Parreño, Assumpta, Perry, Robert J., Shanahan, Phiala E., Wagman, Michael L.

Submitted: 27 July 2023

Subjects:

Journal reference: Phys.Rev.D 108 (11) 114506 (2023)

Citations: 75

Abstract

We present an algorithm to compute correlation functions for systems with the quantum numbers of many identical mesons from lattice quantum chromodynamics (QCD). The algorithm is numerically stable and allows for the computation of <math display="inline"><mi>n</mi></math>-pion correlation functions for <math display="inline"><mi>n</mi><mo>∈</mo><mo stretchy="false">{</mo><mn>1</mn><mo>,</mo><mo>…</mo><mo>,</mo><mi>N</mi><mo stretchy="false">}</mo></math> using a single <math display="inline"><mi>N</mi><mo>×</mo><mi>N</mi></math> matrix decomposition, improving on previous algorithms. We apply the algorithm to calculations of correlation functions with up to 6144 charged pions using two ensembles of gauge field configurations generated with quark masses corresponding to a pion mass <math display="inline"><msub><mi>m</mi><mi>π</mi></msub><mo>=</mo><mn>170</mn><mtext> </mtext><mtext> </mtext><mi>MeV</mi></math> and spacetime volumes of <math display="inline"><mo stretchy="false">(</mo><msup><mn>4.4</mn><mn>3</mn></msup><mo>×</mo><mn>8.8</mn><mo stretchy="false">)</mo><mtext> </mtext><mtext> </mtext><msup><mi>fm</mi><mn>4</mn></msup></math> and <math display="inline"><mo stretchy="false">(</mo><msup><mn>5.8</mn><mn>3</mn></msup><mo>×</mo><mn>11.6</mn><mo stretchy="false">)</mo><mtext> </mtext><mtext> </mtext><msup><mi>fm</mi><mn>4</mn></msup></math>. We also discuss statistical techniques for the analysis of such systems, in which the correlation functions vary over many orders of magnitude. In particular, we observe that the many-pion correlation functions are well-approximated by log-normal distributions, allowing the extraction of the energies of these systems. Using these energies, the large-isospin-density, zero-baryon-density region of the QCD phase diagram is explored. A peak is observed in the energy density at an isospin chemical potential <math display="inline"><msub><mi>μ</mi><mi>I</mi></msub><mo>∼</mo><mn>1.5</mn><msub><mi>m</mi><mi>π</mi></msub></math>, signaling the transition into a Bose-Einstein condensed phase. The isentropic speed of sound, <math display="inline"><msub><mi>c</mi><mi>s</mi></msub></math>, in the medium is seen to exceed the ideal-gas (conformal) limit (<math display="inline"><msubsup><mi>c</mi><mi>s</mi><mn>2</mn></msubsup><mo>≤</mo><mn>1</mn><mo>/</mo><mn>3</mn></math>) over a wide range of chemical potential before falling towards the asymptotic expectation at <math display="inline"><msub><mi>μ</mi><mi>I</mi></msub><mo>∼</mo><mn>15</mn><msub><mi>m</mi><mi>π</mi></msub></math>. These, and other thermodynamic observables, indicate that the isospin chemical potential must be large for the system to be well described by an ideal gas or perturbative QCD.