Charmonium-Glueball spectroscopy with improved hadron creation operators

This paper presents a lattice QCD study of scalar glueball and charmonium mixing that employs improved hadron creation operators, including derivative-based meson and chromo-magnetic field glueball operators, to successfully resolve the low-lying spectrum and identify the lightest iso-scalar state as glueball-dominated.

The Gist

Imagine the universe is built out of tiny, invisible Lego bricks. Most of the bricks we know are "quarks," which snap together to form protons and neutrons (the stuff inside atoms). But there's another kind of brick called a "gluon." Gluons are the super-strong glue that holds the quarks together.

Usually, gluons are just the glue. But physicists think that sometimes, if you have enough of them, they can snap together to form their own little structures, completely without any quarks. These ghostly, glue-only structures are called Glueballs.

The problem? We've never actually seen a Glueball in real life. It's like trying to find a specific, invisible ghost in a crowded room. The math says they should exist, but they are incredibly hard to spot because they look very similar to normal particles made of quarks (like mesons).

This paper is a report from a team of scientists who decided to build a better "flashlight" to find these ghosts. Here is how they did it, explained simply:

1. The Problem: The "Noisy Room"

In the world of quantum physics, scientists use giant supercomputers to simulate the universe. They try to create a particle and see how long it lasts to figure out what it is.

Think of it like trying to hear a whisper in a very noisy room.

• The Whisper: The signal from the particle you are looking for.
• The Noise: Random statistical errors that happen in the computer simulation.
• The Issue: For a long time, the "flashlights" (mathematical tools called operators) the scientists used were too dim. They couldn't distinguish the whisper from the noise, or they confused a Glueball ghost with a normal particle.

2. The Solution: Building Better Flashlights

The authors of this paper decided to upgrade their flashlights. They created two new types of tools to "create" particles in the computer simulation:

A. The "Meson" Flashlight (For normal particles)

Normal particles are like a dance between two partners (quarks). The scientists realized that to see the dance clearly, they needed to look at it from different angles and with different levels of detail.

• The Old Way: They used a standard, blurry lens.
• The New Way: They used a technique called Distillation. Imagine taking a photo of a dancer and applying a "smart filter" that only keeps the most important parts of the image while removing the blur. They combined this with looking at the dancers moving in specific patterns (using derivatives). This made the signal much clearer and the "whisper" much louder.

B. The "Glueball" Flashlight (For the ghost particles)

This was the big innovation. Glueballs are made of magnetic fields (gluons).

• The Old Way: Scientists used "Wilson Loops." Imagine trying to map a complex 3D shape by drawing a single, giant, tangled string around it. It's messy, and if you smooth out the string (a process called smearing), all the different shapes start to look exactly the same. It's like trying to tell a square from a circle by looking at a blurry blob.
• The New Way: The team built their tools using the actual "magnetic field" components and their mathematical derivatives.
  • The Analogy: Instead of drawing a tangled string, they built a 3D sculpture out of specific magnetic bricks. Because they built it with specific shapes and directions (using math called derivatives), the sculpture has a unique "fingerprint." Even if you smooth out the edges, the sculpture still looks distinct. This allowed them to see the Glueball clearly without it getting lost in the noise.

3. The Discovery: Mixing the Signals

The scientists put these two new flashlights together. They created a simulation where normal particles and Glueballs could mix (since they look similar).

• The Result: They successfully identified the lightest Glueball.
• The Character: They found a particle that is mostly made of "glue" (gluons) with a mass of about 1.9 GeV (which is roughly twice the mass of a proton).
• The Excitation: They also found a "sibling" particle that is mostly made of normal quarks.

Why This Matters

Think of this like tuning a radio. Before, the radio was full of static, and you could only hear the loudest stations (normal particles). The Glueball station was so quiet and mixed with the static that no one could find it.

This paper is like the engineers who finally built a noise-canceling headset and a better antenna.

1. They proved that Glueballs can be found in these simulations.
2. They showed that using "smart" mathematical shapes (derivatives) is much better than using old, messy loops.
3. They confirmed that the lightest Glueball is a stable, distinct object, mostly made of glue.

In a nutshell: The team built a better microscope to see the invisible glue that holds the universe together. They proved that "glue-only" particles exist in their simulations and gave us a new, clearer way to find them in the future. This brings us one step closer to solving one of the biggest mysteries in physics: What happens when you have a ball made entirely of glue?

Technical Summary

1. Problem Statement

Lattice Quantum Chromodynamics (QCD) spectroscopy faces two primary challenges when studying states with mixed compositions (e.g., mesons and glueballs):

• Signal-to-Noise Ratio: Glueball correlators suffer from severe statistical noise, often obscuring the signal before the asymptotic regime (where the ground state dominates) is reached.
• Operator Overlap: Extracting the energy spectrum relies on the Generalized Eigenvalue Problem (GEVP). If the set of creation operators used does not sufficiently overlap with the true energy eigenstates, the GEVP fails to resolve excited states or requires large temporal separations where noise dominates.
• Basis Degeneracy: Traditional glueball operators, constructed from spatial Wilson loops with link smearing, often become nearly degenerate (linearly dependent) after smearing, leading to ill-conditioned correlation matrices.
• Flavor Mixing: In the iso-scalar channel ($I=0$), charmonium ($c\bar{c}$) and glueball states mix. Standard approaches often treat these sectors separately or use suboptimal operator bases, making it difficult to disentangle the dominant composition of the resulting eigenstates.

2. Methodology

The authors performed a study in $N_f=2$ QCD using two degenerate clover-improved Wilson quarks at a mass equal to half the physical charm mass. The lattice spacing is $a \approx 0.0658$ fm, and the ensemble size is $48 \times 24^3$.

A. Improved Meson Operators

• Strategy: Combined optimal distillation profiles with derivative-based operators.
• Construction: Used three types of Dirac structures ($\Gamma$) for charmonium:
  1. $\Gamma = I$ (Standard)
  2. $\Gamma = \gamma_i \nabla_i$ (Spatial derivative)
  3. $\Gamma = \gamma_0 \gamma_5 \gamma_i B_i$ (Coupling to chromo-magnetic field)
• Distillation: Instead of a single profile, they applied 7 different Gaussian widths to each operator type, creating a basis of 21 operators. This avoids the need for trial-and-error tuning of the number of distillation vectors.

B. Improved Glueball Operators

• Strategy: Replaced traditional Wilson loops with continuum-like operators constructed from the chromo-magnetic field ($B_i$) and its gauge-covariant derivatives ($D_i$).
• Construction: Operators are built using Clebsch-Gordan coefficients of the cubic group to ensure specific lattice irreducible representations (irreps) and continuum-like angular momentum ($J$) properties.
  • Example ($0^{++}$): $O(t) = \sum_{\vec{x}, i} \text{Tr}(B_i(\vec{x}, t)^2)$.
  • Derivative operators (e.g., involving $D^2_k B_i$) introduce non-trivial spatial structure and relative minus signs, enhancing linear independence.
• Implementation: Used the clover definition for $F_{\mu\nu}$ and symmetric finite differences for derivatives. This approach is computationally efficient and parallelizes well compared to complex Wilson loop shapes.

C. Mixed Analysis

• Constructed a full correlation matrix mixing meson and glueball operators.
• Applied a partial pruning technique: Singular vectors were calculated separately for the meson and glueball blocks to maintain numerical stability while preserving the distinction between the two operator types.
• Solved the GEVP to extract the iso-scalar ($I=0$) and iso-vector ($I=1$) spectra.

3. Key Contributions

1. Novel Glueball Operator Basis: Introduced a systematic construction of glueball operators using $B_i$ and $D_i$ fields. This approach retains continuum angular momentum information, ensures better linear independence (avoiding the degeneracy issues of Wilson loops), and simplifies implementation for parallel computing.
2. Optimized Distillation for Mesons: Demonstrated that combining derivative-based operators with multiple optimal distillation profiles yields faster convergence to ground state plateaus than standard distillation or derivative-only approaches.
3. Flavor-Dependent Operator Optimization: Showed that the optimal operator choice depends on the flavor channel. For the iso-scalar channel, the derivative operator ($\gamma_i \nabla_i$) outperformed the standard operator ($I$), whereas the opposite was true for the iso-vector channel.
4. Systematic Mixing Analysis: Successfully resolved the mixing between charmonium and glueball states in the iso-scalar sector, quantifying the overlap of specific operators with the resulting energy eigenstates.

4. Results

• Iso-Vector Channel ($I=1$): The combined approach (derivatives + profiles) resolved the ground state and first excitation with high precision and faster convergence than previous methods.
• Iso-Scalar Channel ($I=0$):
  • Ground State: Identified as a glueball-dominated state ($0^{++}$). The glueball operators (particularly those coupling to $J=0$) showed large overlaps with this state. The mass was found to be $\approx 1900$ MeV (consistent with quenched predictions of $\sim 1.7-1.8$ GeV, adjusted for the specific quark mass setup).
  • First Excitation: Identified as a meson-dominated state. Charmonium operators showed significantly larger overlaps with this state than the glueball operators did.
  • Stability: Due to the heavy quark mass used ($m_c/2$), the scalar glueball lies below the $2\pi$ threshold, making it stable and allowing for a clean extraction of the spectrum without decay width complications.
• Correlation Matrix Structure: The new glueball operators produced a correlation matrix with a near-block-diagonal structure (grouping operators by their dominant $J$ content), analogous to the success of Hadron Spectrum Collaboration (HadSpec) meson operators. This contrasts with the dense, degenerate matrices produced by Wilson loops.

5. Significance

• Resolution of the Mixing Problem: The study provides a robust framework for disentangling states with the same quantum numbers but different internal compositions (gluonic vs. quark-antiquark). This is crucial for identifying candidates for the scalar glueball in the physical spectrum (e.g., $f_0(1500)$, $f_0(1710)$).
• Methodological Advancement: The proposed glueball operators offer a superior alternative to Wilson loops, addressing issues of linear independence, computational cost, and parallelization. This paves the way for more complex studies, including multi-level sampling techniques in the presence of dynamical quarks.
• Future Applications: The techniques (optimal profiles + derivative operators for mesons; $B_i/D_i$ operators for glueballs) are directly applicable to scattering studies, tetraquark spectroscopy, and the determination of the full scalar meson spectrum in physical QCD. The authors plan to use these tools to study scalar glueball scattering in future work.

In conclusion, the paper demonstrates that constructing operators with continuum-like structures and optimizing their overlap via distillation profiles significantly improves the reliability and precision of lattice QCD spectroscopy, particularly in the challenging regime of glueball-meson mixing.