---
title: 'Quark-Meson Model: Chiral Dynamics'
url: https://www.emergentmind.com/topics/quark-meson-model
type: topic
---

# Quark-Meson Model: Chiral Dynamics

Searching arXiv for recent papers on the quark-meson model and closely related extensions to ground the article in current literature.
{"query":"all:\"quark-meson model\" OR title:\"quark-meson model\"","max_results": 10,"sort_by": "submittedDate","sort_order": "descending"}
{"query":"title:\"quark-meson model\" OR abs:\"quark-meson model\"","max_results": 10,"sort_by": "submittedDate","sort_order": "descending"}
The quark-meson model is a renormalizable low-energy effective model of QCD in which quark fields are coupled by a Yukawa interaction to a scalar \(\sigma\) field and a pseudoscalar pion triplet. In its standard two-flavor form it is the linear sigma model with quarks, designed to realize the chiral symmetry-breaking pattern of QCD and to study chiral restoration at finite temperature and baryon chemical potential. The \(\sigma\) field serves as the chiral order parameter, the pions are the Goldstone-like modes in the broken phase, and the quark sector provides the fermionic degrees of freedom that backreact on mesonic dynamics [1709.05991][1612.03668][2003.03270].

## 1. Canonical formulation and symmetry structure

In the two-flavor formulation, the field content consists of light quarks \(\psi=(u,d)^T\), a scalar meson \(\sigma\), and a pseudoscalar isotriplet \(\boldsymbol{\pi}=(\pi_1,\pi_2,\pi_3)\). A standard Minkowski-space Lagrangian is
\[
{\cal L}= \frac12\big[(\partial_\mu \sigma)^2+(\partial_\mu \boldsymbol{\pi})^2\big]
-\frac12 m^2(\sigma^2+\boldsymbol{\pi}^2)
-\frac{\lambda}{24}(\sigma^2+\boldsymbol{\pi}^2)^2
+h\sigma
+\bar\psi\left[ i\!\not\!\partial +\left(\mu+\frac12\tau_3\mu_I\right)\gamma^0 -g(\sigma+i\gamma^5\boldsymbol{\tau}\cdot\boldsymbol{\pi}) \right]\psi ,
\]
with \(\mu_B=3\mu\) and, in one common setup, \(\mu_I=0\) [1612.03668].

A frequently used FRG truncation is the local potential approximation (LPA), where the Euclidean effective action is taken as
\[
\Gamma_{k} = \int_x \Big\{ \bar{\psi} \left({\partial}\!\!\!\slash + h(\sigma+i\vec{\tau}\cdot\vec{\pi}\gamma_{5}) -\mu \gamma_0 \right)\psi +\frac{1}{2} (\partial_{\mu}\phi)^{2}+U_{k}(\phi^2) - c \sigma \Big\},
\]
with \(\phi\equiv(\sigma,\vec\pi)\) and chirally symmetric ultraviolet potential
\[
U_\Lambda(\phi^2)=\frac{1}{2}m_\Lambda^2\phi^2+\frac{1}{4}\lambda_\Lambda(\phi^2)^2 .
\]
In this formulation the term \(c\sigma\) explicitly breaks chiral symmetry [1709.05991].

At tree level, after shifting \(\sigma=\phi_0+\tilde\sigma\), the masses are
\[
m_\sigma^2 = m^2+\frac{\lambda}{2}\phi_0^2,\qquad
m_\pi^2 = m^2+\frac{\lambda}{6}\phi_0^2,\qquad
m_q=g\phi_0,
\]
and the vacuum minimum is commonly identified with \(f_\pi\), so \(\phi_0=f_\pi\) [1612.03668]. This makes the model technically simple while preserving the essential chiral structure required for equilibrium and nonequilibrium studies.

## 2. Parameter fixing and renormalization

A central technical issue is how the parameters \(m^2\), \(\lambda\), and \(g\) are fixed once loop effects are included. The familiar tree-level matching,
\[
m^2=-\frac12(m_\sigma^2-3m_\pi^2),\qquad
\lambda=3\frac{m_\sigma^2-m_\pi^2}{f_\pi^2},\qquad
g^2=\frac{m_q^2}{f_\pi^2},
\]
is widely used, but it is only a tree-level identity. Once one-loop vacuum fluctuations are included in the effective potential, these become renormalized quantities, and direct tree-level matching becomes inconsistent [1612.03668].

A consistent treatment uses on-shell and \(\overline{\mathrm{MS}}\) renormalization schemes to relate the physical inputs
\[
m_\sigma,\quad m_\pi,\quad m_q,\quad f_\pi
\]
to the running mass parameter and couplings. In the large-\(N_c\) limit employed in that construction, quark self-energy effects are subleading, \(Z_\psi=1\), \(\delta m_q=0\), and only fermion loops are kept in meson self-energies [1612.03668].

An important conceptual point is that the \(\sigma\) mass should be fixed by the physical pole mass, not by the curvature of the effective potential at zero external momentum. The curvature mass corresponds to the self-energy at zero momentum, whereas the pole mass is the physically correct definition used in the on-shell matching procedure [1612.03668]. This distinction is consequential because the phase structure obtained from the one-loop effective potential depends sensitively on whether parameter fixing is performed consistently.

In explicit FRG studies, a representative ultraviolet parameter set has been chosen as
\[
\Lambda=1~\text{GeV},\quad m_\Lambda/\Lambda=0.969,\quad \lambda_\Lambda=0.001,\quad c/\Lambda^3=0.00175,\quad h=4.2,
\]
with the infrared theory reproducing
\[
f_\pi \approx 92.5~\text{MeV},\quad m_\pi = 138~\text{MeV},\quad m_\sigma = 606~\text{MeV},\quad m_\psi = 388~\text{MeV}
\]
in vacuum [1709.05991].

## 3. Nonperturbative implementations

The quark-meson model is used with several complementary approximation schemes, each emphasizing different physics.

The functional renormalization group integrates fluctuations progressively from a UV scale \(\Lambda\) to the infrared through the Wetterich equation,
\[
\partial_k \Gamma_k = \frac{1}{2}\text{Tr}\left\{ \partial_k R_k\left(\Gamma^{(2)}_k+R_k\right)^{-1} \right\}.
\]
In LPA only the potential \(U_k\) flows, while wavefunction renormalizations and higher derivative terms are neglected. This is computationally efficient and often captures critical behavior reasonably well, but it can be insufficient in the low-\(T\), high-\(\mu\) regime where thermodynamic pathologies appear [1709.05991].

A related extension is the Polyakov-loop-extended quark-meson model (PQM), in which quarks couple to a static temporal background gauge field through the Polyakov loop. Its Lagrangian is
\[
{\cal L} = \bar{q}\,\left[i\sl{D} - g(\sigma + i\gamma_5 \vec{\tau}\cdot\vec{\pi})\right]q
+\frac{1}{2}(\partial_\mu\sigma)^2+\frac{1}{2}(\partial_\mu\vec{\pi})^2
-U(\sigma,\vec{\pi})-{\cal U}(\ell,\ell^*),
\]
and in the limit \(\ell,\ell^*\to 1\) it reduces to the standard quark-meson model [1004.2665]. In this formulation, Polyakov-loop-modified quark occupation numbers suppress low-temperature quark excitations and alter conserved-charge fluctuations.

For real-time dynamics, the model has also been formulated using the two-particle irreducible effective action on the Schwinger-Keldysh closed time path. The 2PIEA,
\[
\Gamma[\phi,G,\Delta] = S[\phi] +\frac{i}{2}\Tr\ln G^{-1} +\frac{i}{2}\Tr\!\left[G^{-1}_{\mathrm{cl}(\phi)}G\right] -i\Tr\ln \Delta^{-1} -i\Tr\!\left[\Delta^{-1}_{\mathrm{cl}(\phi)}\Delta\right] +\Gamma_2[\phi,G,\Delta]+\text{const.},
\]
yields causal initial-value equations for the order parameter, bosonic propagators, and fermionic propagators, including memory integrals and off-shell effects [2003.03270]. This framework makes it possible to study thermalization and dynamical spectral rearrangement rather than only equilibrium thermodynamics.

## 4. Thermodynamics and phase structure

The thermodynamic observables follow from the infrared grand potential. In the FRG treatment,
\[
\Omega(T,\mu)=U_{k\to 0},\qquad
p(T,\mu)=-\Omega(T,\mu)+\Omega(0,0),\qquad
s=\frac{\partial p}{\partial T},\qquad
n=\frac{\partial p}{\partial \mu}.
\]
The resulting phase diagram contains a chiral crossover at low \(\mu\) and higher \(T\), a first-order chiral transition at large \(\mu\) and low \(T\), and a critical end point near
\[
T \approx 50~\text{MeV},\qquad \mu \approx 345~\text{MeV}
\]
in the specific FRG-LPA setup quoted above [1709.05991].

A notable result of that FRG analysis is that the first-order line is not monotonic. Part of it has negative slope, but at low temperature it bends in the opposite direction and develops positive curvature. The paper analyzes this using the Clausius-Clapeyron relation,
\[
\frac{dT_c}{d\mu_c}=-\frac{\Delta n}{\Delta s},
\]
so that, if \(\Delta n>0\), a positive slope implies \(\Delta s<0\). The model then develops a region where the pressure decreases with increasing temperature and the entropy density becomes negative. This is described as thermodynamically unacceptable and is interpreted not as a physical state of quark matter but as evidence that the homogeneous FRG treatment in LPA is incomplete in that regime [1709.05991].

The dependence on approximation scheme is substantial. In the one-loop large-\(N_c\) analysis with consistent on-shell parameter fixing, including vacuum fluctuations shifts the phase boundary relative to the no-sea approximation. For \(m_\sigma=600\) MeV, the transition is first order at low \(T\) and high \(\mu\), terminating at a tricritical point located at
\[
(\mu,T)\approx (303.24~\text{MeV},\,55~\text{MeV}),
\]
and becoming second order on the low-\(\mu\) side. For \(m_\sigma=800\) MeV, the transition is second order throughout the entire \(\mu\)–\(T\) plane when vacuum fluctuations are included, whereas neglecting them restores a first-order line [1612.03668]. This contrast shows that the quark-meson phase diagram is not a unique output of the field content alone; it is strongly conditioned by renormalization, fluctuation content, and truncation.

In the PQM version at physical pion mass, the transition becomes a smooth crossover rather than a true phase transition. Mesonic fluctuations further smooth thermodynamic observables and broaden peaks in susceptibilities near the pseudocritical region [1004.2665].

## 5. Spectral functions, nonequilibrium evolution, and instability diagnostics

The quark-meson model is not restricted to static order parameters. In FRG-based real-time analyses, retarded meson two-point functions are obtained by analytic continuation of Euclidean correlators,
\[
\Gamma^{(2),R}(\omega,\vec p)
=
-\lim_{\epsilon\to 0}\Gamma^{(2),E}(p_0=-i(\omega+i\epsilon),\vec p),
\]
and the spectral function is constructed as
\[
\rho(\omega,\vec p)=\frac{1}{\pi}\frac{\text{Im}\,\Gamma^{(2),R}(\omega,\vec p)}{\left(\text{Re}\,\Gamma^{(2),R}(\omega,\vec p)\right)^2+\left(\text{Im}\,\Gamma^{(2),R}(\omega,\vec p)\right)^2}.
\]
These observables are used to diagnose possible instabilities of the homogeneous ground state [1709.05991].

One diagnostic is the static pion inverse propagator at finite momentum. The condition
\[
D_\pi^{-1}(\omega=0,|\vec p_c|)=0,
\]
equivalently
\[
\Gamma^{(2)}_{k,\pi}(\omega=0,|\vec p_c|)=0,
\]
at finite spatial momentum is interpreted as evidence for an instability toward a spatially modulated phase, such as a chiral spiral or chiral density wave. In the FRG study, regions where such zero crossings occur overlap with the low-\(T\), high-\(\mu\) domain in which the negative-entropy pathology appears, suggesting that inhomogeneous phases may be relevant [1709.05991].

The nonequilibrium 2PI treatment formulates the dynamics in terms of statistical and spectral functions for both bosons and fermions. The Wigner-transformed spectral function,
\[
\rho(\tau,\omega,|\mathbf p|) = \int_{-2\tau}^{2\tau}d(\Delta t)\,e^{i\omega\Delta t}\rho(\tau,\Delta t,|\mathbf p|),
\]
separates relative-time structure from macroscopic evolution [2003.03270]. In thermal equilibrium, the fluctuation-dissipation relation links \(F\) and \(\rho\); dynamically, the simulations show that this relation emerges around \(\tau\gtrsim 100\), with temperatures stable by \(\tau=130\) [2003.03270].

Spectrally, the bosonic channels are well described by relativistic Breit-Wigner forms, while the fermionic vector-zero component develops a second low-frequency mode in the infrared in the crossover region. This produces a double-peak structure interpreted as an additional light propagating fermionic degree of freedom near the crossover. At higher temperature the peaks overlap and the spectrum again resembles a single broadened quasiparticle peak [2003.03270]. These results make the model a laboratory for studying not only equilibrium symmetry restoration but also how thermal quasiparticles and collective modes emerge dynamically.

## 6. Extensions, related families, and nomenclature

The name “quark-meson model” is often used narrowly for the chiral quark-\(\sigma\)-\(\pi\) theory, but several related model families extend it by adding additional effective degrees of freedom or by changing the microscopic interpretation of the hadronic sector.

| Model | Additional ingredient | Main role |
|---|---|---|
| QM | \(\sigma\), \(\pi\), quarks | Chiral symmetry breaking and restoration |
| PQM | Polyakov loop background \(\ell,\ell^*\) | Approximate confinement/deconfinement physics |
| QMD | Explicit diquark fields \(\Delta\) | 2SC and CFL color-superconducting phases |
| QMC | Composite baryons with internal quark structure | Nuclear matter, finite nuclei, neutron stars |

The Polyakov-loop extension preserves the quark-meson chiral sector but supplements it with a background gluonic potential \({\cal U}(\ell,\ell^*)\). In this framework, the kurtosis \(R_{4,2}=c_4/c_2\) drops from about \(9\) at low temperature to below \(1\) around the crossover region, reflecting the change from three-quark baryonic excitations at low \(T\) to single-quark-like degrees of freedom at high \(T\). By contrast, the QM model without Polyakov loop gives \(R_{4,2}=1\) at low temperature because quark excitations are not suppressed in the same way [1004.2665].

The quark-meson diquark model (QMD) is a renormalizable extension of the QM model that includes explicit color-antitriplet diquark fields. It is designed to describe pion condensation at finite isospin chemical potential and the 2SC and CFL phases at finite baryon chemical potential. In this model all symmetries, including \(SU(N_c)\), are global, so symmetry breaking in the 2SC phase produces physical Goldstone bosons rather than Higgsed gauge fields; the paper classifies the resulting type-A and type-B Goldstone modes and shows that the BCS gaps approach constants for large chemical potentials while the speed of sound approaches the conformal value from above [2602.18256].

A different source of terminological confusion is the quark-meson coupling model (QMC). Despite the similar acronym, QMC is not the chiral quark-meson model of QCD thermodynamics. It is a relativistic mean-field theory for nuclei and dense matter in which baryons are composite objects made of confined quarks and the meson fields couple directly to the quarks inside the baryons. Its characteristic formula is a nonlinear effective baryon mass, for example
\[
M_N^*=M_N-g_\sigma\sigma+\frac{d}{2}(g_\sigma\sigma)^2,
\]
with scalar polarizability \(d\) encoding the internal quark response [1802.08368][1912.12461]. This distinction is essential: the standard QM model is a low-energy chiral model of quarks and mesons, whereas QMC is a quark-structure model of baryons in nuclear matter.

## 7. Limitations, open issues, and significance

The quark-meson model is valued because it isolates the interplay of quark dynamics and chiral symmetry in a renormalizable setting, but its limitations are equally prominent. The low-\(T\), high-\(\mu\) FRG study finds a broad region with negative entropy density, upward bending of the first-order line, and pressure decreasing with temperature. These features are interpreted as signs that the homogeneous LPA truncation is missing essential physics rather than discovering a stable new phase [1709.05991].

Several missing ingredients have been identified within the model literature. One possibility is a truncation or regulator artifact: higher-order derivative terms could alter the problematic thermodynamics. Another is color superconductivity: the attractive interaction induced by the mesonic sector supports the standard two-flavor color-superconducting channel
\[
\langle\psi^TCi\gamma_5\tau_2\lambda_2\psi\rangle ,
\]
so the negative-entropy region may be a symptom of a missing 2SC instability. A third is the formation of inhomogeneous chiral condensates, suggested by finite-momentum zero crossings of the pion inverse propagator. The same work also notes that low-temperature matter should probably include nucleonic degrees of freedom, which the homogeneous quark-meson model does not contain [1709.05991].

At a broader level, the model’s predictive content depends strongly on the fluctuation scheme. Mesonic fluctuations smooth the crossover and reduce the sharpness of higher susceptibilities, Polyakov-loop backgrounds are needed to reproduce hadronic suppression of low-temperature quark-number fluctuations, and consistent one-loop on-shell parameter fixing materially changes the phase diagram relative to the no-sea approximation [1004.2665][1612.03668]. This suggests that the quark-meson model is best understood not as a single immutable theory but as a structured family of effective descriptions whose reliability is regime-dependent.

Within that domain, it remains one of the standard controlled settings for studying chiral symmetry breaking, crossover dynamics, spectral rearrangement, conserved-charge fluctuations, and the interplay between thermodynamics and collective excitations in low-energy QCD [1709.05991][2003.03270].

Source: https://www.emergentmind.com/topics/quark-meson-model