Axial-anomaly effects and chiral phase structure in holographic QCD

1. Problem Statement

Understanding the chiral phase structure of Quantum Chromodynamics (QCD), particularly the nature of the chiral phase transition at finite temperature, remains a central challenge. A key tool for this analysis is the Columbia plot, which maps the order of the phase transition (first-order, second-order, or crossover) as a function of light ($m_{u,d}$) and strange ($m_s$) quark masses.

A critical theoretical uncertainty in this domain is the role of the axial $U(1)_A$ anomaly. While the anomaly explicitly breaks $U(1)_A$ symmetry in QCD, its effective implementation in phenomenological models varies significantly. Previous holographic QCD studies (specifically soft-wall models) often restricted the flavor sector to $SU(2)$ or $SU(3)$, neglecting the pseudoscalar singlet sector ($\eta'$ meson) or treating the anomaly interaction as a constant. This limitation prevents a dynamical description of the $\eta-\eta'$ system and obscures how different realizations of the anomaly strength affect the finite-temperature phase diagram.

Core Question: How do different functional forms (profiles) of the axial-anomaly interaction, constrained by vacuum observables, influence the predicted chiral phase structure and the topology of the Columbia plot?

2. Methodology

The authors employ a bottom-up holographic QCD approach based on the AdS/CFT correspondence, specifically extending a previous soft-wall model to a $U(3)$ flavor framework.

• Model Construction:

  • Geometry: 5D AdS space with a soft-wall dilaton profile $\Phi(z)$ to ensure linear Regge trajectories.
  • Fields: A bifundamental scalar field $X$ representing chiral symmetry breaking, coupled to $U(3)_L \times U(3)_R$ gauge fields.
  • Anomaly Implementation: The authors introduce a holographic-coordinate-dependent determinant interaction term: $-\gamma(z) \text{Re}\{\det X\}$. This term dynamically generates the $\eta'$ mass and induces mixing between the flavor-singlet ($\eta_0$) and flavor-octet ($\eta_8$) components.
  • Anomaly Profiles: Three distinct ansätze for the strength function $\gamma(z)$ are tested:
    • Type A: $\gamma_A(z) = \gamma_0 z^2$ (monotonically increasing).
    • Type B: $\gamma_B(z) = \frac{\gamma_0 z^2}{1 + \gamma_1 z^2}$ (saturates in the infrared).
    • Type C: $\gamma_C(z) = \gamma_0 z^2 e^{-\gamma_1 z^2}$ (peaked at intermediate scales).

• Vacuum Constraints:

  • The model parameters are fitted to reproduce the vacuum pseudoscalar spectrum (masses of $\pi, K, \eta, \eta'$ and their excitations) and mixing patterns.
  • A simultaneous fit is attempted for both scalar and pseudoscalar sectors to further constrain the anomaly profile, specifically testing the Type-B profile.

• Finite-Temperature Analysis:

  • The background geometry is replaced by an AdS black-hole metric to introduce temperature $T$.
  • The order of the chiral phase transition is determined by analyzing the behavior of the chiral condensates ($\sigma_u, \sigma_s$) as functions of temperature. Discontinuities indicate first-order transitions; smooth variations indicate crossovers.
  • Columbia Plots are constructed by scanning the light and strange quark masses to map the regions of different transition orders.

3. Key Contributions

1. $U(3)$ Extension: The paper successfully extends soft-wall holographic QCD to include the full $U(3)$ flavor group, enabling a dynamical description of the $\eta-\eta'$ system without imposing mixing angles by hand.
2. Profile Sensitivity: It demonstrates that while vacuum pseudoscalar observables (masses and mixing) constrain the overall magnitude of the anomaly, they do not uniquely determine the functional form of the anomaly profile $\gamma(z)$. Multiple profiles can yield equally good vacuum fits.
3. Scalar Sector Constraints: The study reveals that including scalar meson observables imposes stricter constraints. Only the Type-B profile allows for a simultaneous reasonable fit of both scalar and pseudoscalar spectra, though this introduces tension in decay constants.
4. Phase Structure Divergence: The most significant finding is that different anomaly profiles, which are indistinguishable in vacuum spectroscopy, lead to qualitatively different chiral phase structures at finite temperature.

4. Key Results

• Vacuum Sector:

  • All three profiles (A, B, C) successfully reproduce the $\eta-\eta'$ mass splitting and mixing patterns within experimental uncertainties.
  • The Type-B profile is the only one capable of providing a simultaneous fit to both scalar and pseudoscalar meson spectra, although it results in decay constants ($f_\pi, f_K$) that are lower than experimental values.
  • The scalar sector exhibits a reversed singlet-octet ordering compared to the pseudoscalar sector, a non-trivial consequence of the determinant interaction.

• Finite-Temperature Phase Structure (Columbia Plot):

  • Type-A Profile: Predicts a crossover or second-order transition across the entire quark-mass plane. The anomaly strength effectively suppresses near the critical temperature, eliminating the first-order region.
  • Type-B and Type-C Profiles: Predict a first-order phase transition region in the light-quark mass corner (small $m_{u,d}$), consistent with the conventional Columbia plot scenario.
  • Sensitivity: The size of the first-order region varies significantly between profiles. The Type-B profile (with scalar constraints) yields the largest first-order region among the viable cases.

• Lattice Comparison:

  • The results align with recent lattice QCD discussions suggesting that the existence and size of the first-order region depend on the persistence of $U(1)_A$ breaking near the critical temperature. The Type-A profile mimics scenarios where the anomaly is effectively restored/suppressed at $T_c$.

5. Significance

• Theoretical Insight: The paper highlights a critical limitation in effective QCD models: vacuum spectroscopy alone is insufficient to fix the anomaly sector's behavior at finite temperature. The "Columbia plot" is highly sensitive to the phenomenological implementation of the anomaly.
• Model Validation: It underscores the necessity of incorporating finite-temperature and topological observables (like topological susceptibility) alongside vacuum data to constrain holographic models effectively.
• Phenomenological Implications: The findings suggest that the nature of the chiral phase transition in QCD (specifically the existence of a first-order corner) may hinge on the detailed energy-scale dependence of the axial anomaly, a feature that current lattice simulations are still refining.
• Future Directions: The authors propose that future work should focus on observables directly sensitive to the anomaly (e.g., topological susceptibility) and extend the analysis to finite baryon chemical potential to further test the robustness of these phase structures.

In summary, this work provides a controlled holographic framework demonstrating that the modeling of axial-anomaly effects is the decisive factor in determining the global chiral phase structure of QCD, even when vacuum constraints are satisfied.