🔬 physics

Meson mass spectrum in QCD $_2$ 't Hooft's model

This paper generalizes analytical methods to derive the meson mass spectrum in large- $N_c$ QCD $_2$ with arbitrary equal quark masses, providing analytical expressions for spectral sums and a systematic large- $n$ expansion that aligns with known results in both chiral and heavy quark limits.

Original authors: Alexey Litvinov, Pavel Meshcheriakov

Published 2026-10-30

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Alexey Litvinov, Pavel Meshcheriakov

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the universe of particle physics, there is a persistent mystery known as confinement. It is the reason we never see a single quark floating freely in nature. Quarks are the tiny building blocks of protons and neutrons, but they are bound together by a force so powerful that trying to pull them apart is like stretching a rubber band that never snaps; instead, the energy you put in creates new particles before the connection breaks. Understanding exactly how this binding works is one of the most difficult challenges in modern science. While our best theories describe the behavior of these particles in the three dimensions of space plus time we inhabit, solving the equations for these interactions is often impossible with current mathematical tools. To get around this, physicists often turn to simplified versions of reality, known as toy models. These are stripped-down universes where the rules are similar to our own but the math is manageable enough to solve. By studying these simpler worlds, scientists hope to uncover the fundamental mechanisms that govern the complex, real world.

One such simplified universe is a version of quantum chromodynamics, the theory of the strong nuclear force, reduced to just one dimension of space and one of time. In this two-dimensional world, the behavior of quarks and the mesons they form is governed by a famous equation developed by physicist Gerard 't Hooft. This equation acts as a map, describing the possible energy levels, or masses, of the particles that can exist in this model. For decades, scientists have been able to calculate these masses using powerful computers, but finding a direct mathematical formula to predict them has remained elusive. The difficulty lies in the fact that the equation involves a complex relationship between the masses of the quarks and the energy of the resulting particle, a relationship that changes depending on how heavy the quarks are.

In a recent study, researchers Alexey Litvinov and Pavel Meshcheriakov tackled this problem by focusing on a specific, symmetrical case where the two quarks making up a meson have the exact same mass. They sought to find a way to describe the entire spectrum of possible meson masses using pure mathematics, rather than relying solely on computer simulations. Their approach involved translating the original, difficult integral equation into a different mathematical language known as integrability. This translation allowed them to use a set of powerful tools developed for other areas of physics to analyze the problem. By doing so, they were able to derive new, precise formulas that describe the masses of these particles, particularly for the lighter, lower-energy states and the heavier, higher-energy states.

The researchers found that they could express the sum of all possible meson masses in a specific way, creating a kind of mathematical inventory of the particle spectrum. They developed detailed formulas that work for any mass of the quarks, provided the two are equal. These formulas are not just rough estimates; they are systematic expansions that become more accurate as you look at higher energy levels. The team verified their results by comparing their new mathematical predictions against high-precision computer calculations. The agreement was striking, with the new formulas matching the computer data to a very high degree of accuracy, confirming that their mathematical approach correctly captures the physics of the model.

One of the most significant achievements of this work was extending the understanding of these particles to two extreme physical situations. The first is the "chiral limit," a scenario where the quarks become effectively massless. In this state, the theory predicts the existence of a special, massless particle similar to the pion in our own universe, which plays a crucial role in the structure of matter. The researchers showed that their formulas naturally lead to this result, correctly predicting how the mass of this particle behaves as the quark mass approaches zero. The second extreme is the "heavy quark limit," where the quarks are incredibly massive. In this regime, the mesons behave differently, and the researchers' formulas successfully described how the energy levels shift and settle as the quark mass increases.

The study does not claim to have solved the confinement problem for our real, four-dimensional universe. Instead, it provides a complete and rigorous solution for a specific, simplified model that shares many features with the real world. The authors have demonstrated that by using advanced mathematical techniques, it is possible to move beyond numerical approximations and find exact analytical expressions for the properties of these bound states. Their work offers a deeper understanding of how mass emerges from the interaction of fundamental particles in a confined system. By generalizing previous methods to cover a wider range of conditions, they have provided a more robust toolkit for physicists to explore the mechanics of confinement. The results stand as a verified bridge between abstract mathematical theory and the concrete numerical behavior of particles, offering a clear view of how these systems function when the variables are controlled and the masses are equal.

Technical Summary: Meson mass spectrum in QCD2 't Hooft's model

Problem Statement
The paper addresses the determination of the meson mass spectrum in two-dimensional Quantum Chromodynamics (QCD $_2$ ) in the large- $N_c$ limit, governed by 't Hooft's integral equation. While the equation is solvable numerically, the authors aim to extend analytical methods to the case of equal but arbitrary quark masses ( $m_1 = m_2$ , corresponding to the parameter $\alpha_1 = \alpha_2 = \alpha$ ). Previous analytical work by Fateev, Lukyanov, and Zamolodchikov (FLZ) [2] was restricted to the massless case ( $\alpha = 0$ ). The challenge lies in generalizing the integrability-based approach to arbitrary $\alpha$ to derive analytical expressions for spectral sums and the large- $n$ expansion of the spectrum.

Methodology
The authors adopt the integrability framework established by FLZ, reformulating the 't Hooft integral equation into Baxter's TQ equation. The core strategy involves:

Fourier Transform: Transforming the integral equation into momentum space using the variable $\frac{1}{2} \log \frac{x}{1-x}$ , leading to a functional difference equation for the wave function $\Psi(\nu)$ .
Baxter's TQ Equation: Defining a function $Q(\nu)$ such that the problem reduces to solving the functional difference equation $Q(\nu + 2i) + Q(\nu - 2i) - 2Q(\nu) = -\frac{2z}{\nu + \alpha x}Q(\nu)$ , where $z$ and $x$ depend on $\nu$ and $\alpha$ .
Limiting Cases: Solving the TQ equation in two regimes:
- Small $\lambda$ (Low energy): Expanding solutions in powers of $\lambda$ using confluent hypergeometric functions. A key technical innovation involves introducing a modified digamma-like function $\psi_\alpha(\nu)$ to handle branch points at $\nu = \pm i$ that arise when $\alpha \neq 0$ .
- Large $\lambda$ (High energy): Constructing an asymptotic expansion in powers of $\lambda^{-1}$ using a series related to the Barnes function and an auxiliary function $S_0(\nu)$ defined via an integral representation.
Spectral Determinants: Utilizing relations between the spectral determinants $D_\pm(\lambda)$ and the solutions $Q_\pm(\nu)$ to extract spectral data. The authors propose generalized relations (2.32) and (2.33) connecting the logarithmic derivatives of the spectral determinants to the derivatives of the $Q$ -functions at specific points, extending the FLZ formula (2.31) to $\alpha \neq 0$ .

Key Contributions and Results

Generalization of FLZ Methods: The paper successfully generalizes the analytical machinery of FLZ to the case of equal non-zero quark masses ( $\alpha \neq 0$ ). This includes deriving the specific form of the TQ equation solutions and the necessary auxiliary functions.
Spectral Sums: The authors derive analytical expressions for the spectral sums $G_\pm^{(s)} = \sum \lambda_n^{-s}$ for $s \le 7$ . These expressions are given in terms of the parameter $\alpha$ , Riemann zeta values, and specific integrals $u_{2k-1}(\alpha)$ and $i_{2k}(\alpha)$ .
Large- $n$ Expansion: A systematic expansion for the eigenvalues $\lambda_n$ (and thus meson masses $M_n^2$ ) is derived up to high orders in $1/n$ . The expansion includes terms involving logarithms of $n$ and depends on the quark mass parameter $\alpha$ .
Limiting Cases:
- Chiral Limit ( $\alpha \to -1$ ): The authors analyze the behavior of spectral sums as the quark mass vanishes. They demonstrate that even spectral sums diverge (indicating a massless ground state, the pion) while odd sums remain finite. The derived ground state mass satisfies the Gell-Mann-Oakes-Renner (GMOR) relation.
- Heavy Quark Limit ( $\alpha \to \infty$ ): The asymptotic behavior of spectral sums is derived, showing agreement with results obtained from simplified heavy-quark approximations in the literature.
Numerical Verification: The analytical results are rigorously tested against numerical solutions of the 't Hooft equation using a basis of Chebyshev polynomials. The relative error between the analytical large- $n$ expansion (truncated at $O(\lambda^{-6})$ ) and numerical eigenvalues is found to be extremely small (often $< 10^{-4}$ for $n \ge 2$ ). The spectral sums and the ratio of normalization constants $d_-/d_+$ are also verified numerically with high precision.

Significance
The paper claims to provide a complete analytical framework for the meson spectrum in 't Hooft's model with equal quark masses, bridging the gap between the massless case and the general massive case. By verifying the results numerically and recovering known limits (chiral and heavy quark), the authors establish the validity of their generalized integrability approach. The derived formulas for spectral sums and the large- $n$ expansion offer precise tools for studying the bound state spectrum without relying solely on numerical diagonalization. The authors note that the analytic structure of these sums as functions of the complex parameter $\alpha$ presents a new field for theoretical investigation, potentially linking to non-unitary Conformal Field Theories (CFTs). Future work is suggested to extend these results to the case of unequal quark masses ( $\alpha_1 \neq \alpha_2$ ).

More like this