Meson spectrum and low-energy constants in large-$N$ QCD

This paper presents new non-perturbative results on the meson spectrum and low-energy constants in large-$N$ QCD, derived from lattice Monte Carlo simulations of the Twisted Eguchi-Kawai model up to $N=841$ combined with standard finite-$N$ data.

The Gist

Imagine the universe is built from tiny, invisible Lego bricks. In the world of particle physics, the most important bricks are quarks and gluons. When they stick together, they form particles called mesons (like the pion or the rho meson).

For decades, physicists have tried to understand exactly how these bricks snap together. The rules of the game are written in a theory called Quantum Chromodynamics (QCD). But here's the problem: calculating how these bricks interact is incredibly hard. It's like trying to predict the weather in a hurricane by looking at a single raindrop. The math gets so messy that even the world's fastest supercomputers struggle to solve it.

The "Magic Trick": The Large-N Limit

To make the math manageable, physicists use a clever trick. They imagine a version of the universe where the number of "colors" (a property of quarks, not actual colors) is not 3 (like in our real world), but infinity.

Think of it like this:

• Our Real World: Imagine a small, chaotic dance floor with 3 dancers. It's hard to predict who bumps into whom.
• The "Large-N" World: Imagine a massive stadium filled with thousands of dancers. When the crowd is huge, the chaos smooths out into a predictable, flowing pattern.

By studying this "infinite crowd" version first, physicists can find the underlying rules of the dance. Then, they can work backward to understand our real world with just 3 dancers.

The Problem: The "One-Point" Box

Usually, to simulate this on a computer, you need a huge grid (a lattice) to represent space and time. But if you want to simulate a universe with thousands of "colors," the computer needs to be impossibly large. It's like trying to build a model of the entire stadium on a single postage stamp.

Enter the Twisted Eguchi-Kawai (TEK) model. This is the paper's main innovation.

The authors discovered a way to shrink the entire stadium down to a single point without losing any of the physics.

• The Analogy: Imagine you have a giant, complex tapestry. Usually, you need a huge loom to weave it. The TEK model is like a magical loom that can weave the entire tapestry using just one thread, provided you twist that thread in a very specific, intricate way.
• By using this "twist," they could simulate a universe with 841 "colors" (a massive number) on a computer that would normally only handle a tiny fraction of that.

What They Found

Using this super-efficient "one-point" simulation, the authors calculated two main things:

1. The "Musical Notes" of the Universe (Meson Spectrum)
Just as a guitar string vibrates at specific frequencies to create notes, particles vibrate at specific masses. The authors mapped out the "notes" (masses) of these particles in their infinite-color universe.

• They found that the heavier, excited particles follow a beautiful, predictable pattern called a Regge trajectory.
• The Metaphor: Imagine a ladder. The rungs represent different particle masses. The authors found that in the infinite-color world, the rungs are perfectly evenly spaced, like a well-tuned musical scale. This helps them understand how the "ladder" looks in our real world.

2. The "Glue" Constants (Low-Energy Constants)
To describe how these particles interact, physicists use a set of numbers called "Low-Energy Constants." Think of these as the recipe ingredients for the strong force.

• How sticky are the quarks? (The Chiral Condensate)
• How fast do they decay? (The Pion Decay Constant)
• The authors calculated these ingredients with extreme precision.
• The Surprise: When they compared their "infinite-color" results with real-world data (where there are only 3 colors), they found that the "infinite" version is actually a very good approximation, but there are some subtle differences. It's like tasting a soup made with a million spices versus a soup made with just three; the flavor profile is similar, but the nuances are different.

Why This Matters

This paper is a breakthrough because it proves that we can simulate the "infinite" universe on a computer with incredible accuracy.

• Before: We had to guess how the infinite universe behaved by looking at small, imperfect simulations.
• Now: We have a direct window into the "perfect" version of the theory.

By understanding the "perfect" infinite version, we can better understand the messy, real world. It's like studying the physics of a perfect, frictionless sphere to understand how a bumpy, real-world ball rolls down a hill.

In a Nutshell

The authors built a mathematical time machine that shrinks the entire universe into a single point, allowing them to simulate a version of reality with thousands of particle types. They used this to map out the "musical notes" of the subatomic world and measure the "ingredients" of the strong force, giving us a clearer picture of how the universe is built from the bottom up.

Technical Summary

1. Problem Statement

The paper addresses the challenge of obtaining quantitative, non-perturbative results for Quantum Chromodynamics (QCD) in the 't Hooft large-$N$ limit (where the number of colors $N \to \infty$ and the 't Hooft coupling $\lambda = Ng^2$ is fixed).

• Theoretical Context: In the large-$N$ limit with $N_f/N \to 0$, QCD simplifies significantly. Gluon dynamics dominate at leading order, while quark contributions are sub-leading ($1/N$ suppressed). This limit allows for a reorganization of the diagrammatic expansion and offers insights into non-perturbative aspects of strong interactions.
• Limitations of Standard Methods: Analytical methods in large-$N$ QCD are limited. Standard lattice QCD simulations typically extrapolate results from small $N$ (usually $N \leq 10$) to $N \to \infty$ assuming a $1/N$ expansion. However, finite-$N$ corrections can be large, making such extrapolations potentially unreliable or misleading.
• Goal: The authors aim to provide direct, non-perturbative data for meson spectra and Low-Energy Constants (LECs) at very large $N$ values (up to $N=841$) to benchmark theoretical predictions and improve the understanding of the $1/N$ expansion.

2. Methodology

The study utilizes Lattice Monte Carlo simulations based on the Twisted Eguchi–Kawai (TEK) model, which relies on the concept of large-$N$ volume independence.

• Volume Reduction: Instead of simulating on a large space-time lattice, the TEK model reduces the space-time volume to a single point (a one-point box). This is valid if center symmetry is unbroken.
• Twisted Boundary Conditions: To enforce center symmetry in a single-point box, the authors employ twisted boundary conditions in all directions. The action is the TEK Wilson plaquette action:
  $$S_W[U] = -Nb \sum_{\mu \neq \nu} z_{\nu\mu} \text{Tr}\{U_\mu U_\nu U_\mu^\dagger U_\nu^\dagger\}$$
  where $z_{\nu\mu}$ are twist factors and $b = 1/(Ng^2)$ is the inverse 't Hooft coupling.
• Effective Volume: In this setup, the effective physical volume scales as $V_{\text{eff}} \propto N^2$. The effective size is $\ell = a\sqrt{N}$, where $a$ is the lattice spacing. This allows simulations at $N=841$ to probe physics equivalent to much larger volumes in standard lattice setups.
• Quenched Approximation: The simulations are performed in the quenched approximation ($N_f=0$), meaning quark loops are neglected (the quark determinant is dropped). Quarks are treated as dynamical probes on fixed gluon backgrounds.
• Fermion Discretization: The authors use the Wilson Dirac operator adapted for the TEK formulation. Fermionic observables are computed using the inverse of this operator.
• Data Analysis:
  • Meson Masses: Extracted via the Generalized Eigenvalue Problem (GEVP) analysis of correlation matrices. Correlators are computed in momentum space and transformed to coordinate space.
  • Extrapolation: Results are extrapolated to the continuum limit ($a \to 0$) and the chiral limit ($m_\pi \to 0$) using a fit function that accounts for $O(a)$ lattice artifacts and leading $O(m_\pi^2)$ chiral behavior.
  • $1/N$ Expansion: The TEK results (at $N=\infty$ effectively) are combined with standard finite-$N$ results from literature to determine the coefficients of the $1/N$ expansion for various observables.

3. Key Contributions

• Record-Breaking $N$ Values: The study presents results for $N$ up to 841, significantly surpassing previous lattice studies which typically stopped around $N=361$.
• Direct Access to Large-$N$ Limit: By using the TEK model, the authors bypass the need for uncertain extrapolations from small $N$, providing a direct anchor point for the $N \to \infty$ limit.
• Comprehensive Spectrum: The paper provides a detailed spectrum of mesons (pions, rhos, and their excited states) and determines the coefficients of the $1/N$ expansion for key Low-Energy Constants (LECs) up to $O(1/N^3)$.
• Validation of Volume Independence: The work demonstrates the efficacy of the TEK model for studying meson dynamics and chiral properties in large-$N$ QCD.

4. Key Results

A. Meson Spectrum and Regge Trajectories

• Spectrum: The authors computed masses for ground and excited states ($\pi, \rho, a_0, a_1, b_1$, etc.) in the chiral limit.
• Finite-$N$ Effects: Comparing large-$N$ results with experimental data and finite-$N$ simulations reveals that finite-$N$ corrections are small for low-lying states but grow significantly for higher excited states in the meson tower.
• Radial Regge Trajectories: The squared masses of excited mesons ($m_n^2$) follow a linear Regge trajectory:
  $$\left(\frac{m_{A_n}^{(\chi)}}{\sqrt{\sigma}}\right)^2 = C + \frac{\mu_r^2}{\sigma} n$$
  The extracted Regge slopes ($\mu_r$) are:
  • $\pi$-channel: $\mu_r/\sqrt{\sigma} = 3.65(21)$
  • $\rho$-channel: $\mu_r/\sqrt{\sigma} = 3.95(24)$
  • These values are in excellent agreement with large-$N$ Polyakov effective model predictions ($\approx 3.55$) and are larger than phenomenological values derived from experimental masses ($\approx 2.65$).

B. Low-Energy Constants (LECs)

The authors determined the $N \to \infty$ values for several LECs (renormalized in the $\overline{\text{MS}}$ scheme at 2 GeV):

• Chiral Condensate: $\frac{\Sigma_R}{N\sqrt{\sigma}^3} = 0.0889(23)$
• Pion Decay Constant: $\frac{F_\pi}{\sqrt{N}\sqrt{\sigma}} = 0.1262(34)$
• Chiral Slope: $\frac{B_R}{\sqrt{\sigma}} = 5.58(26)$
• NLO Coupling: $\frac{\bar{\ell}_4}{N} = 0.446(55)$

C. The $1/N$ Expansion

By combining their $N=\infty$ results with finite-$N$ data (specifically for $N_f=2$), the authors performed polynomial fits to extract the coefficients of the $1/N$ expansion.

• Critical Finding: They found that sub-leading terms in the $1/N$ expansion are surprisingly large when dynamical fermions are present.
• Implication: Extrapolating from small finite-$N$ data (e.g., $N=3, 4, 5$) to $N=\infty$ without including the large-$N$ anchor point can be misleading. The inclusion of the TEK results significantly alters the extrapolated values, particularly for the chiral condensate and $\bar{\ell}_4$.
• Convergence: The results highlight the convergence of $SU(N_f)$ and $U(N_f)$ chiral effective theories in the large-$N$ limit, where the $\eta'$ meson becomes light and degenerate with pions, reducing sub-leading corrections in the $U(N_f)$ case.

5. Significance

This work represents a major step forward in the non-perturbative understanding of QCD in the large-$N$ limit.

1. Benchmarking Theory: The precise determination of meson masses and LECs at $N=841$ provides a rigorous benchmark for analytical models (like holographic QCD or effective string theories) that rely on the large-$N$ limit.
2. Correcting Extrapolations: The study demonstrates that relying solely on small-$N$ lattice data to infer large-$N$ physics is risky due to large sub-leading corrections. The TEK model is essential for obtaining the true large-$N$ baseline.
3. Phenomenological Insight: The results clarify the behavior of the chiral condensate and pion decay constant in the large-$N$ limit, offering insights into the nature of confinement and chiral symmetry breaking.
4. Future Directions: The success of the TEK model for meson dynamics paves the way for studying more complex phenomena, such as $\pi$-$\pi$ scattering, in the large-$N$ limit, which is currently a frontier in lattice QCD.