⚛️ nuclear theory

Structure of Heavy Mesons in the Light-Front Quark Model

This paper demonstrates that while both the Single Gaussian Ansatz and Gaussian Expansion Method (GEM) effectively model static properties of heavy mesons, GEM offers a superior description of structural observables by correctly reproducing perturbative QCD asymptotic behaviors and endpoint distributions in wave functions and distribution amplitudes.

Original authors: Ahmad Jafar Arifi, Lucas Happ, Shuhei Ohno, Makoto Oka

Published 2026-10-27

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Ahmad Jafar Arifi, Lucas Happ, Shuhei Ohno, Makoto Oka

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

Inside the heart of every atom, protons and neutrons are built from smaller particles called quarks, which are held together by the strong force, the most powerful interaction in nature. To understand how these particles stick together to form the matter we see around us, physicists use a framework called Quantum Chromodynamics. However, solving the equations for this theory is incredibly difficult, especially when trying to describe how quarks move at speeds close to the speed of light. To bridge this gap, scientists often rely on simplified models that capture the essential behavior of these particles without getting lost in impossible calculations. One such approach is the light-front quark model, a tool that treats hadrons—particles like protons and mesons—as collections of quarks moving within a specific, simplified view of time and space. This method has been successful in describing many properties of matter, but researchers have long debated how best to describe the internal shape and movement of the quarks inside these particles, particularly the heavy ones containing charm or bottom quarks.

A team of researchers at RIKEN and other institutions in Japan recently set out to test two different ways of describing the internal structure of these heavy mesons. They focused on the ground state, which is the most stable and lowest-energy configuration of these particles. The scientists compared two mathematical approaches used to map out the wave function, a description of where the quarks are likely to be found and how they move. The first approach, known as the Single Gaussian Ansatz, uses a simple, bell-shaped curve to estimate the particle's shape. It is a straightforward method that has served the field well for decades. The second approach, called the Gaussian Expansion Method, is more complex. Instead of relying on a single curve, it builds the description by combining many different bell-shaped curves with varying widths, allowing it to adapt to a much wider range of shapes and behaviors.

The researchers applied both methods to a set of heavy mesons, including those made of charm and bottom quarks. They first tuned the parameters of their models to match known experimental data, such as the mass of the particles and how quickly they decay. In this initial test, both methods performed remarkably well. They were able to reproduce the known masses and decay rates with high accuracy, showing that the simple bell-curve approach is still a powerful tool for predicting basic properties. This result confirmed that for static characteristics—those that do not change with time or motion—the simple model is sufficient.

However, the story changed when the team looked deeper into the structural details of the particles. They examined how the quarks are distributed inside the meson, specifically looking at the edges of the particle where the momentum of the quarks becomes very high. Here, the two methods began to diverge. The simple bell-curve model failed to capture the correct behavior at these extreme edges. It showed a suppression, or a drop-off, that did not match what theory predicts should happen at high speeds. In contrast, the more complex expansion method revealed a different picture. It showed that the wave function extends further and behaves in a way that aligns perfectly with the predictions of perturbative Quantum Chromodynamics, the theory that describes how quarks interact at very short distances. Specifically, the complex method showed that the probability of finding a quark at the very edge of the particle follows a specific power-law pattern, a behavior that the simple model missed entirely.

This difference in the internal structure had real consequences for how the particles interact with light. The researchers calculated the electromagnetic form factors, which describe how the particle responds to an electric charge. They found that the complex method predicted a faster drop-off in this response at low energy levels compared to the simple method. This faster drop-off brought the predictions closer to the data obtained from lattice Quantum Chromodynamics, a highly detailed computer simulation that is considered a gold standard for these calculations. The simple model, while good for basic numbers, was less accurate in describing these subtle structural nuances.

The study concludes that while the simple, single-curve approach is adequate for calculating the mass and decay of heavy mesons, it falls short when describing their detailed internal architecture. The more sophisticated expansion method, by allowing for a richer variety of shapes, captures the correct behavior of quarks at high speeds and short distances. This suggests that to truly understand the structure of matter, especially in the realm of heavy particles, scientists need models that can adapt to the complex, power-law behaviors dictated by the fundamental laws of physics. The findings highlight that the choice of mathematical tool matters deeply: a simple approximation can get the big picture right, but only a more flexible approach can reveal the true nature of the particle's edge.

Technical Summary: Structure of Heavy Mesons in the Light-Front Quark Model

Problem Statement
The paper addresses the limitations of current Light-Front Quark Model (LFQM) approaches in describing the structural properties of hadrons, specifically ground-state heavy mesons. While the Single Gaussian Ansatz (SGA) has been successful in reproducing static observables like mass spectra and decay constants, it relies on a fixed functional form (a single Gaussian) that does not represent the true eigenstate of the Hamiltonian. Consequently, the SGA fails to accurately capture the asymptotic behaviors of wave functions (WFs) and distribution amplitudes (DAs), particularly near the endpoints ( $x \to 0, 1$ ) and in the high-momentum region. These regions are governed by short-range interactions and relativistic kinematics, which are critical for consistency with perturbative QCD (pQCD) predictions. The authors aim to assess these limitations by contrasting the SGA with the Gaussian Expansion Method (GEM), which aims to approximate the true eigenstate of the Hamiltonian more rigorously.

Methodology
The authors investigate heavy mesons (containing $c$ or $b$ quarks) using a relativized quark model within the instant form, which is subsequently mapped to light-front dynamics.

Hamiltonian: The study employs a QCD-motivated interquark potential consisting of a confining term, a color Coulomb term, and a hyperfine interaction. The kinetic energy is treated relativistically.
Two Approaches:
- Single Gaussian Ansatz (SGA): Assumes the meson wave function is a single Gaussian function. The width parameter is determined via the variational principle.
- Gaussian Expansion Method (GEM): Expands the wave function as a linear superposition of multiple Gaussian basis functions with varying range parameters. This method solves the generalized eigenvalue problem to find the coefficients that minimize the energy, effectively approximating the Hamiltonian's eigenstate.
Parameter Determination: For both methods, model parameters (quark masses, confinement parameters, strong coupling, and smearing parameters) are determined independently by performing a $\chi^2$ fit to experimental mass spectra and decay constants, as well as lattice QCD data.
Mapping to LFQM: The resulting instant-form wave functions are mapped to Light-Front Wave Functions (LFWFs) using the Melosh transformation. This allows for the calculation of structural observables.
Observables: The study calculates and compares static properties (mass spectra, decay constants) and structural properties, including:
- Twist-2 Distribution Amplitudes (DAs).
- Electromagnetic (EM) form factors in the low- and high- $Q^2$ regions.
- Asymptotic behaviors of the wave functions in position and momentum space.

Key Contributions and Results

Static Properties: Both SGA and GEM yield comparable accuracy in reproducing static properties. The $\chi^2/N_{dof}$ values indicate that both models fit the mass spectra and decay constants well, though the fitted model parameters (e.g., quark masses and potential strengths) differ significantly between the two approaches.
Wave Function Asymptotics: The study reveals a critical distinction in the short-distance (high-momentum) behavior of the wave functions.
- GEM: Accurately reproduces the power-law damping of the wave function in momentum space, $\psi(k_\perp \to \infty) \propto 1/k_\perp^2$ , which aligns with pQCD predictions derived from one-gluon exchange and relativistic kinematics.
- SGA: Fails to reproduce this power-law behavior, exhibiting a Gaussian fall-off instead.
Distribution Amplitudes (DAs): The DAs derived from GEM show enhanced amplitudes near the endpoints ( $x \to 1$ ) compared to SGA. Specifically, the GEM results exhibit the correct asymptotic behavior $\phi(x \to 1) \propto (1-x)$ , consistent with pQCD. In contrast, the SGA suppresses the DAs near the endpoints, deviating from lattice QCD data and theoretical expectations.
Electromagnetic Form Factors: In the low- $Q^2$ region, the EM form factors for pseudoscalar mesons calculated using GEM fall off faster than those from SGA. This faster fall-off leads to better agreement with lattice QCD data, particularly for the $\eta_c$ meson. The charge radii calculated via GEM are generally larger than those from SGA, reflecting the extended spatial distribution of the GEM wave function.
Model Parameters: The authors note that the parameters fitted for GEM and SGA are distinct. Using parameters optimized for SGA within a GEM framework (or vice versa) leads to significant deviations (up to 30%) in decay constants, underscoring that the two methods are conceptually different and require independent tuning.

Significance
The paper concludes that while the SGA is a simple and powerful tool for describing static hadron properties, it possesses inherent limitations in capturing the detailed structural features and asymptotic behaviors required for a complete description of heavy mesons. The GEM, by approximating the true eigenstate of the Hamiltonian, provides a more accurate description of the wave function's shape, particularly in the short-range and high-momentum regions. This leads to DAs and form factors that are more consistent with perturbative QCD expectations and lattice QCD data. The authors emphasize that the endpoint behaviors of DAs and the power-law damping of wave functions are intrinsically linked to the short-range color Coulomb interaction and relativistic kinematics, which the GEM captures more effectively than the fixed-form SGA. The study suggests that for structural observables sensitive to wave function details, methods that aim to reproduce the Hamiltonian's eigenstate (like GEM) are superior to simple Ansatz-based approaches.

More like this