Bethe–Salpeter Analysis in Meson Spectroscopy

• Bethe–Salpeter analysis is a continuum QFT approach that couples Dyson–Schwinger and Bethe–Salpeter equations to determine meson bound state properties.
• The method employs symmetry-preserving truncations, such as rainbow–ladder and Ball–Chiu vertices, validated through Ward–Takahashi identities and effective interactions like the Maris–Tandy model.
• Numerical results verify the approach’s accuracy by reliably matching light-flavor meson masses and ensuring chiral symmetry is systematically maintained.

The Bethe–Salpeter analysis is a continuum quantum-field-theoretical approach for determining the properties of composite bound states, such as mesons, through the coupled solution of Dyson–Schwinger equations (DSE) for single-particle propagators and the Bethe–Salpeter equation (BSE) for two-body amplitudes. The implementation in "Masses of Light Flavor Mesons using Bethe–Salpeter Approach" formalizes this analysis for light-flavor meson masses, utilizing symmetry-preserving truncations anchored by the Maris–Tandy model for nonperturbative gluon dynamics and systematic control of chiral symmetry through Ward–Takahashi identities (Liaqat et al., 4 Nov 2025).

1. Bethe–Salpeter Equations: Homogeneous and Inhomogeneous Formalism

The central BSE formalism exploits both inhomogeneous and homogeneous structures. The total meson momentum is denoted $P$ ($P^2=-m^2$), with internal relative momentum $p$, and individual quark momenta in the Bethe–Salpeter loop are $k_\pm = k \pm \tfrac12 P$.

The inhomogeneous BSE for a given channel reads: $$\Gamma_{M}(p; P) = \Gamma_0(p; P) + \int \frac{d^4k}{(2\pi)^4} K(p, k; P)\left[S(k_+)\Gamma_M(k; P)S(k_-)\right]$$ where $\Gamma_0$ is the channel-dependent Dirac structure (e.g., $\gamma_5$ or $\gamma_5\gamma_\mu$), $S$ is the dressed quark propagator, and $K$ is the quark–antiquark interaction kernel.

The homogeneous (bound-state) BSE is an eigenvalue problem: $$\Gamma_M(p; P) = \lambda(P^2) \int \frac{d^4k}{(2\pi)^4} K(p, k; P) S(k_+) \Gamma_M(k; P) S(k_-)$$ Physical meson masses correspond to $P^2$ where $\lambda(P^2)=1$.

2. Rainbow–Ladder Truncation, Maris–Tandy Model, and Vertex Treatments

Quark Propagator DSE and Maris–Tandy Gluon Exchange

The DSE for the renormalized quark propagator is: $$S^{-1}(p) = Z_2(i\!\not p + m_0) + Z_1g^2 \int^\Lambda \frac{d^4k}{(2\pi)^4} D_{\mu\nu}(q) \frac{\lambda^c}{2}\gamma^\mu S(k) \Gamma_\nu^c(k, p)$$ where $D_{\mu\nu}(q)$ is the gluon propagator and $\Gamma_\nu^c$ the quark–gluon vertex. In rainbow–ladder, the vertex is set to bare, and the effective interaction is replaced by the Maris–Tandy model: $$Z_1g^2 D_{\mu\nu}(q) \to \frac{g(q^2)}{q^2}\left(\delta_{\mu\nu}-\frac{q_\mu q_\nu}{q^2}\right)$$ with

$$\frac{g(q^2)}{q^2} = 4\pi^2 D \frac{\omega^6 q^2}{\omega^2} e^{-q^2/\omega^2} + 4\pi \frac{\gamma_m \pi}{\tfrac12 \ln[\tau + (1 + q^2/\Lambda_{\text{QCD}}^2)^2]} f(q^2)$$

where $f(q^2) = \frac{1-e^{-q^2/(4 m_t^2)}}{q^2}$, $\gamma_m = \frac{12}{33-2N_f}$, $\tau = e^2-1$.

Vertex Ansätze

• Bare vertex (RL): $\Gamma_\mu(p,k) = \gamma_\mu$
• Ball–Chiu (BC) vertex: Minimal longitudinal solution to vector Ward–Green–Takahashi identity,

$$i\Gamma_\mu^{\rm BC}(k, p) = i \bar\Delta_A(k^2, p^2) \gamma_\mu + 2\ell_\mu\Big[i\gamma\!\cdot\!\ell \Delta_A + \Delta_B\Big]$$

where the dressing functions $A, B$ originate from $S^{-1}(p) = i\!\not p\,A(p^2)+B(p^2)$.

Both the RL (bare) and BC vertex choices guarantee vector WGTI preservation in the gap and Bethe–Salpeter equations.

3. Axial-Vector Ward–Takahashi Identity and Truncation Consistency

Chiral symmetry, encoded by the axial-vector WGTI,

$$P_\mu \Gamma_{5\mu}(p; P) = S^{-1}(p_+)i\gamma_5 + i\gamma_5 S^{-1}(p_-) - i[m_f(\zeta)+m_g(\zeta)]\Gamma_5(p; P)$$

must not be violated by truncation. The BSE kernel $K$ must cohere with the quark–gluon vertex and the DSE truncation to ensure this identity holds, which is achieved automatically in the RL scheme ($\Lambda_{5\mu\beta}=0$), but requires explicit nonzero corrections for BC-improved truncations.

4. Mass Extraction via Padé Approximation

Rather than directly solving the BSE in the timelike regime, the amplitude is computed for spacelike $P^2>0$. For pseudoscalar channels, the leading Dirac component $E(p^2, P^2)$ is extracted. The inverse $P_E(P^2) = 1/E(p^2=0, P^2)$ is fitted over a range of spacelike $P^2$ values to a diagonal Padé approximant: $$f_N(x) = \frac{c_0 + c_1x + \cdots + c_Nx^N}{1 + c_{N+1}x + \cdots + c_{2N}x^N}$$ The analytically continued pole at $P^2 = -m^2$ gives the ground-state meson mass.

5. Numerical Results: Meson Spectrum Across RL and BC Truncations

Parameters ($\omega, D$) are tuned so that chiral condensate and pion properties match experiment. Numerical outcomes for both truncation schemes are summarized below:

Vertex	ω	D	$(-\langle\bar{q}q\rangle)^{1/3}$	$m_\pi$	$m_K$	$m_\rho$	$m_\sigma$
RL	0.40	0.93	0.277	0.139	0.501	0.814	0.670
BC	0.40	0.406	0.307	0.140	0.474	—	—

The masses from the homogeneous RL (eigenvalue) calculation also corroborate: $$m_\pi = 0.138,\quad m_K = 0.494,\quad m_\rho = 0.802,\quad m_\sigma = 0.668$$

There is robust agreement between RL and BC vertex truncations regarding the pion mass, with minor shifts in the kaon and vector/scalar channels. This stability validates the approach's treatment of axial and vector symmetries in the truncation.

6. Systematic Bethe–Salpeter–Dyson–Schwinger Workflow

The calculation proceeds as:

• DSE for the quark propagator: Truncate using Maris–Tandy interaction, construct $S(p)$.
• Vertex ansatz selection: RL (bare) and BC solutions guarantee required WGTI symmetry.
• Homogeneous BSE: Solve for meson masses using eigenvalue cross-check ($\lambda(P^2)=1$).
• Inhomogeneous BSE: Spacelike calculation of leading amplitude component, Padé fit for mass extraction.
• Compare computed spectrum for both truncations, demonstrating mild quantitative sensitivity and systematic symmetry preservation.

This approach demonstrates that symmetry-respecting coupled DSE–BSE truncations, implemented with validated nonperturbative interaction models and controlled through vector and axial WGTIs, furnish a quantitatively accurate description of light-meson masses in continuum QCD (Liaqat et al., 4 Nov 2025).