Density-Dependent Quark Mass Framework

• Density-Dependent Quark Mass (DDQM) framework is a model for deconfined quark matter that employs effective quark masses dependent on density to mimic confinement and asymptotic freedom.
• It uses a canonical ensemble formulation to derive energy density, pressure, and chemical potential rigorously, incorporating a density-dependent bag term to enforce stability.
• The framework advances astrophysical research by providing a tool for predicting self-bound strange stars, hybrid star phases, and stability windows in dense quark matter.

The density-dependent quark mass (DDQM) framework, also referred to as the quark mass density-dependent (QMDD) model, is a phenomenological approach for modeling deconfined quark matter at high baryon density and low temperature. Originally formulated to mimic key features of QCD—such as confinement and asymptotic freedom—via a density dependence in the quark effective mass, the framework has evolved to incorporate a rigorous thermodynamic structure, parameter-space diagnostics, astrophysical applications, and various refinements to address longstanding inconsistencies. Below is a comprehensive technical overview of the DDQM formalism, emphasizing the self-consistent canonical-ensemble construction and its implications for the phase structure and stability of dense quark matter (Lugones et al., 2022).

1. Canonical Ensemble Formulation and Key Equations

The foundational advance in modern DDQM is the adoption of the canonical, rather than grand canonical, ensemble. For a system at zero temperature ($T=0$) with $N$ particles in volume $V$, one posits an effective mass ansatz

$M(n) = m_0 + D\, n^{-\alpha},$

where $n = N/V$ is the number density, $m_0$ the current (bare) mass, and $D,\alpha>0$ are parameters set by phenomenological or stability constraints. This form ensures $M(n)\rightarrow\infty$ as $n\rightarrow 0$ (confinement) and $M(n)\rightarrow m_0$ at asymptotically large densities (asymptotic freedom).

The canonical free energy per unit volume is

$\epsilon(n) = g\, M(n)^4\, \chi[x(n)],$

where $g$ is the degeneracy and

$$x(n) = \frac{k_F}{M(n)},\qquad k_F = \left(\frac{6\pi^2 n}{g}\right)^{1/3},$$

with

$$\chi(x) = \frac{1}{16\pi^2} \left[x\sqrt{x^2+1}(2x^2 + 1) - \mathrm{arcsinh}\,x\right].$$

The pressure is not given by the naive Fermi-gas result due to the density dependence of $M(n)$. It is instead obtained as

$$P(n) = n^2 \frac{d}{dn}\left(\frac{\epsilon(n)}{n}\right) = g\, M^4 \phi(x) - B(n),$$

with

$$\phi(x) = \frac{1}{48\pi^2}\left[x\sqrt{x^2+1}(2x^2-3) + 3\mathrm{arcsinh}\,x\right],$$

$B(n) = -g\, n\, M^3 \frac{dM}{dn}\, \beta(x),$

$$\beta(x) = 4\chi(x) - x\chi'(x) = \frac{1}{4\pi^2}\left[x\sqrt{x^2+1} - \mathrm{arcsinh}\,x\right].$$

The term $B(n)>0$ functions as a density-dependent “bag constant,” automatically enforcing confinement at low density.

The chemical potential exhibits a similar splitting: $$\mu(n) = \frac{d\epsilon}{dn} = M \sqrt{x^2+1} - \frac{B(n)}{n}.$$

The number density, energy density, and other observables can subsequently be constructed.

2. Thermodynamic Consistency and the Energy/Pressure Structure

The canonical ensemble construction ensures rigorous thermodynamic consistency. The Euler relation,

$n\, \frac{d\epsilon}{dn} = \epsilon + P,$

is identically satisfied due to the explicit structure of $\epsilon(n)$ and $P(n)$. Critically, the stationarity of the energy-per-baryon occurs precisely at zero pressure: $$\frac{d}{dn_B}\left(\frac{\epsilon}{n_B}\right) = 0 \iff P=0,$$ where $n_B = n/3$ is the baryon number density. This property, missing in earlier grand-canonical treatments, guarantees a well-defined minimum for $E/A$ at $P=0$, a prerequisite for meaningful stability and phase separation analyses.

The absolute-stability condition for quark matter is set by comparison to iron: $$\left.\frac{E}{A}\right|_{P=0} < 930\,\textrm{MeV},$$ ensuring self-bound quark phases suitable for strange star scenarios.

3. Extension to β-Equilibrated, Charge-Neutral Three-Flavor Matter

For astrophysically relevant applications, the construction must incorporate fermion flavor structure, electric neutrality, and weak equilibrium. Densities $n_i$ for $i \in \{u, d, s, e\}$ and $n_B = (n_u + n_d + n_s)/3$ are considered, with

$$\epsilon = \sum_{i=u,d,s} g\, M_i^4\, \chi[x_i] + g_e\, m_e^4\, \chi[x_e].$$

Weak equilibrium conditions: $\mu_d = \mu_u + \mu_e,\quad \mu_s = \mu_d,$ and electric neutrality: $$\frac{2}{3}n_u - \frac{1}{3}n_d - \frac{1}{3}n_s - n_e = 0,$$ are imposed. Solving these four coupled equations at a trial $n_B$ yields all chemical potentials and partial densities, after which the pressure is evaluated by summing all species contributions.

Two mass formula architectures are often used:

• Flavor-dependent masses (Model 1): $M_i(n_i) = m_i + D n_i^{-\alpha}$,
• Flavor-blind masses (Model 2): $M_i(n_B) = m_i + D n_B^{-\alpha}$.

The pressure and chemical potential structure for each species is

$$P_i = g M_i^4 \phi(x_i) - B_i,\qquad \mu_i = M_i \sqrt{x_i^2+1} - \Delta_i,$$

with $B_i$ and $\Delta_i$ defined according to the chosen mass formula.

4. Parameter Space and Phase Structure: Stability, Self-Bound, and Hybrid Windows

The $(D, \alpha)$ parameter space can be analytically partitioned using

$$e_{ud}^0 = \left.\frac{\epsilon+P}{n_B}\right|_{\textrm{2f}}, \qquad e_{uds}^0 = \left.\frac{\epsilon+P}{n_B}\right|_{\textrm{3f}},$$

defining critical curves at $e_{ud}^0=930$ MeV and $e_{uds}^0=930$ MeV.

• Below both curves: self-bound quark matter (two- or three-flavor).
• Between: "2f-only" self-bound matter.
• Above both: hybrid matter (quark only at high pressure).
• The $\alpha$ value at which $c_s(P=0)=c$ marks the causality boundary; self-bound regions above this are excluded.

These sharp boundaries allow identification of windows for stable strange stars, hybrid stars, and causal EoSs.

5. Analytical Summary and Comparison with Alternative Models

The DDQM framework in canonical ensemble yields the following analytical results:

• Free energy per volume:

$\epsilon(n) = g\, M(n)^4\, \chi[x(n)].$

• Pressure and chemical potential:

$P(n) = g\, M^4\, \phi(x) - B(n),$

$\mu(n) = M \sqrt{x^2+1} - \frac{B(n)}{n},$

$B(n) = -g\, n\, M^3\, \frac{dM}{dn}\, \beta(x).$

• Thermodynamic consistency: $n\mu = \epsilon + P$.
• Minimum $E/A$ at $P=0$.

Inclusion of a density-dependent mass modifies the pressure and chemical potential so as to naturally generate a "bag-like" term and incorporate the effects of repulsive interactions. Importantly, as shown in both the canonical and grand-canonical consistent frameworks, the resulting EoS and stability window for strange quark matter nearly coincide with those of the MIT bag model once vacuum corrections are included, but are now rooted in dynamical, density-dependent mass generation (0801.2813).

6. Physical Implications and Applications

The rigorously defined DDQM model provides a tractable equation of state for deconfined quark matter, suitable for:

• Predicting properties and the very existence of self-bound strange stars,
• Delineating the boundaries between hadronic, hybrid, and quark-matter phases,
• Systematically exploring parameter regimes where absolutely stable strange matter can exist,
• Ensuring causality in the EoS, with transparent diagnostic criteria.

In realistic applications, simultaneous beta equilibrium and charge neutrality are treated self-consistently, and the full EoS $P(\epsilon)$ can be used as input for Tolman–Oppenheimer–Volkoff integration in compact-star models, and for comparative analysis against astrophysical mass–radius observations, tidal deformability constraints, and multi-messenger signals.

7. Advances, Limitations, and Extensions

• The canonical ensemble solution removes longstanding thermodynamic inconsistencies present in earlier grand-canonical implementations.
• The density-dependent "bag" term $B(n)$, emergent from explicit mass differentiation, is an intrinsic realization of confinement physics; it is not merely a phenomenological constant.
• The model supports both flavor-dependent and flavor-blind mass formulae, enabling exploration of a variety of microphysical scenarios.
• Extensions incorporate isospin effects, magnetic fields, excluded volume corrections, and order parameter approaches (e.g., Polyakov-loop–driven deconfinement), further broadening the physical scope and astrophysical relevance of the framework.
• The phenomenological nature of the mass-density relation remains a limitation; the underlying microscopic QCD derivation is absent, and matching to lattice QCD or perturbative calculations at high density requires additional developments.

The current DDQM framework thus offers a fully self-consistent, thermodynamically robust approach to modeling dense quark matter in compact astrophysical objects and provides the technical basis for systematic phenomenological and observational exploration of deconfined, self-bound matter in the universe (Lugones et al., 2022).