Interacting Boson System at Finite Temperature: The treatment of the lattice calculations

Technical Summary: Interacting Boson System at Finite Temperature

Problem Statement
The paper addresses the theoretical description of hot and dense hadronic matter, specifically a relativistic system of charged bosons (modeled as pions) at finite temperature and isospin density. A primary challenge in this regime is the occurrence of Bose-Einstein condensation (BEC) at sufficiently large isospin densities. In the condensed phase, the isospin chemical potential ($\mu_I$) becomes fixed by the condensation condition (reaching the lowest single-particle energy) and ceases to be an independent thermodynamic variable. Standard Grand Canonical Ensemble approaches, which treat $\mu_I$ as the independent variable, become inadequate for describing the thermodynamics in this phase. Furthermore, for interacting systems, the quasiparticle mass is a dynamical quantity determined self-consistently by the scalar density, complicating the formulation. The authors aim to test whether a thermodynamically consistent mean-field framework can reproduce finite-isospin lattice QCD thermodynamics, a regime accessible to Monte Carlo simulations due to the absence of the fermion sign problem.

Methodology
The authors develop a thermodynamically consistent mean-field treatment for repulsively interacting relativistic bosons using an Extended Canonical Ensemble (ECE).

• Formalism: Instead of using the isospin chemical potential ($\mu_I$) as the independent variable, the ECE formulation treats the conserved isospin density ($n_I$) as the independent variable. This is achieved via a Legendre transformation of the grand potential $\Omega$ to the free energy $F(T, n_I, V)$.
• Interactions: The study considers two self-interaction models for a complex scalar field: a repulsive $\phi^4$ interaction ($\lambda$) and a repulsive $\phi^6$ interaction ($b$). These interactions stabilize the condensate at finite isospin density.
• Mean-Field Approximation: The interaction terms are expanded around the mean scalar density $\sigma = \langle \hat{\phi}^\dagger \hat{\phi} \rangle$. This yields an effective quasiparticle mass $M(\sigma)$ and an excess pressure contribution $P_{ex}(\sigma)$. The condition $\sigma(\partial U/\partial \sigma) = \partial P_{ex}/\partial \sigma$ ensures thermodynamic consistency.
• Condensed Phase Treatment: In the condensed phase, the field is decomposed using the Bogolyubov method ($\hat{\phi} = \Phi_0 + \hat{\psi}$), separating the condensate density ($n_{cond}$) from thermal excitations. The chemical potential is constrained by the onset condition $\mu_I = M(\sigma)$. The system of equations for $n_I$, $\sigma$, and $n_{cond}$ is solved self-consistently without varying $\mu_I$ independently.
• Normalization: To compare with lattice QCD results, the authors adopt the lattice convention for isospin chemical potential normalization, where the condensation condition is $\mu_I = m_\pi/2$ (as opposed to $\mu_I = m_\pi$ in other conventions).

Key Contributions

1. Extended Canonical Ensemble (ECE) Implementation: The paper provides a rigorous framework for handling interacting bosonic systems in the condensed phase by shifting the independent variable from chemical potential to conserved density. This avoids the singularity issues associated with fixing $\mu_I$ during condensation.
2. Thermodynamic Consistency: The authors demonstrate that the mean-field approach, combined with the ECE, yields a self-consistent set of equations for the effective mass, scalar density, and condensate density that satisfies thermodynamic relations.
3. Lattice QCD Comparison: The formalism is applied to a pion-like system at $T = 122$ MeV. The interaction parameters ($\lambda$ or $b$) are fitted to the lattice pressure data using a single fitting parameter.

Results
The models were tested against lattice QCD results for pressure, isospin density, energy density, and trace anomaly at $T = 122$ MeV.

• Pressure and Density: Both the $\phi^4$ and $\phi^6$ models reproduce the general behavior of the lattice pressure and isospin density. The transition between the thermal phase (red segment in plots) and the condensed phase (blue segment) is continuous, with the phase transition point corresponding to the onset of Bose-Einstein condensation.
• Model Accuracy: While both models perform reasonably well, the $\phi^6$ model provides a visibly better overall description of the lattice results compared to the $\phi^4$ model.
• Predictive Power: Crucially, once the interaction constant is fitted to the pressure, the models successfully predict the isospin density, energy density, and trace anomaly without further adjustment. This non-trivial agreement validates the thermodynamic consistency of the framework.
• Effective Mass: The effective quasiparticle mass at the critical point is determined to be $M(\sigma_{cr}) = m_\pi/2$, consistent with the lattice condensation condition.

Significance and Claims
The authors claim that the Extended Canonical Ensemble formulation is a useful and necessary framework for studying interacting relativistic bosonic matter, particularly across the transition between thermal and Bose-condensed regimes. The success of the interacting pion-like models in reproducing lattice-QCD thermodynamics suggests that, at "hadron" temperatures (specifically $T = 122$ MeV) and the considered isospin densities, the vast majority of quark configurations reflect a pion structure (where $n_{\bar{d}} = n_u$ and $n_{\bar{u}} = n_d$). The work demonstrates that effective bosonic descriptions, when formulated with thermodynamic consistency and proper variable selection, can successfully connect with first-principles lattice QCD results in the finite isospin density regime.