⚛️ phenomenology

Hyperfine Structure of $B$ and $D$ Mesons in a QCD-Inspired Potential Model

This paper investigates the hyperfine structure of $B$ , $D$ , and $B_c$ mesons using a QCD-inspired Cornell potential model with Gaussian-smeared spin-spin interactions and Dalgarno-Lewis analytic wavefunctions to provide stable, testable benchmarks for upcoming spectroscopy measurements.

Original authors: Krishna Kingkar Pathak

Published 2026-09-25

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Krishna Kingkar Pathak

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

Deep within the fabric of the universe, matter is held together by a force so powerful it binds the smallest known building blocks into the protons and neutrons that make up our world. This force, known as the strong interaction, is governed by a set of rules called Quantum Chromodynamics. While scientists have a complete mathematical description of this force, solving the equations to understand how particles behave inside a particle is like trying to predict the weather by calculating the motion of every single air molecule; it is computationally impossible with current technology. To bridge this gap, physicists use simplified models that capture the essential behavior of these interactions, treating particles like tiny balls connected by invisible springs. These models allow researchers to predict the properties of exotic particles called mesons, which are made of a heavy quark and a lighter antiquark, and to understand why they have the specific masses and energy levels they do.

One of the most intriguing features of these particles is a phenomenon called hyperfine splitting. Imagine two tiny magnets inside a particle; depending on whether they point in the same direction or opposite directions, the particle's total energy changes slightly. This tiny difference in energy, caused by the spin of the quarks, creates a measurable gap between two otherwise identical particles. For decades, scientists have used potential models to calculate these gaps, but they have faced a persistent problem: the math breaks down when the particles get too close, predicting infinite energy where there should be a finite value. Furthermore, while these models work reasonably well for the lowest energy states of particles, they often struggle to predict the behavior of excited states, which are like higher notes on a musical instrument, or to handle systems where the two constituent particles have very different masses.

In a recent study, Krishna Kingkar Pathak from Arya Vidyapeeth College in India tackled these challenges by refining the mathematical tools used to describe these particles. The researcher focused on a specific type of particle called the Cornell potential, which combines a short-range attraction with a long-range stretching force, much like a rubber band that pulls back harder the more you stretch it. To fix the mathematical breakdown at the center of the particle, Pathak introduced a smoothing technique. Instead of treating the interaction between the quarks as a sharp, singular point that causes the equations to explode, the model spreads this interaction out over a tiny, finite region. This adjustment removes the infinite values and allows the equations to produce stable, realistic results. By using a specific mathematical method to solve the equations, the researcher was able to derive clear, analytic descriptions of the particle's internal structure without needing to rely on massive, brute-force computer simulations.

The study applied this refined model to a wide range of mesons, including those containing charm and bottom quarks, which are significantly heavier than the up and down quarks found in ordinary matter. The results showed that the model could reproduce the known masses and energy gaps of these particles with remarkable precision, often matching experimental measurements within a few million electron volts. This level of accuracy is significant because it confirms that the underlying physics of the model correctly captures the dynamics of the strong force, even when the particles involved have very different masses. The model successfully predicted the energy gaps for the ground states of these particles, aligning closely with data collected by major particle physics experiments around the world.

Beyond just confirming what is already known, the study ventured into uncharted territory by making predictions for particles in higher energy states, specifically the first excited states. These excited states are harder to observe experimentally and have been a source of uncertainty in previous theories. The model predicted that the energy gaps for these excited states would be smaller than those for the ground states, a trend that matches the physical intuition that particles in higher energy orbits are less tightly bound. The study also extended these predictions to the $B_c$ meson, a rare particle made of a bottom quark and a charm quark, for which experimental data is scarce. The predictions for this system followed the same consistent patterns observed in other mesons, suggesting that the model's parameters are robust and reliable across different types of quark combinations.

The research also examined how sensitive these predictions are to the specific values chosen for the strength of the strong force within the model. By varying this strength slightly, the researcher found that the predicted energy gaps remained stable, changing by less than three percent. This stability indicates that the model is not fragile; it does not collapse if the input numbers are tweaked slightly, which gives scientists confidence in its reliability. The study also clarified the physical meaning of the smoothing parameter used to fix the mathematical singularity, linking it to the effective size of the quarks themselves. This connection helps ground the mathematical adjustments in physical reality, suggesting that the "smearing" is not merely a mathematical technique to make the numbers work, but a reflection of the finite size of the particles involved.

Ultimately, this work provides a unified and accurate framework for understanding the internal structure of heavy mesons. By combining a smoothed interaction with a precise mathematical solution, the model offers a clear picture of how quarks bind together and how their spins influence the particle's energy. The predictions for the excited states and the rare $B_c$ meson serve as valuable benchmarks for future experiments at facilities like the LHCb and Belle II. As physicists continue to map the spectrum of these particles, this study offers a reliable guide, ensuring that the theoretical tools used to interpret new data are as solid and consistent as the experimental measurements themselves.

Technical Summary: Hyperfine Structure of B and D Mesons in a QCD-Inspired Potential Model

Problem Statement
Quantum Chromodynamics (QCD) remains difficult to solve directly in the non-perturbative low-energy regime. While potential models, particularly the Cornell potential, offer intuitive descriptions of quark-antiquark dynamics, they face specific challenges in accurately reproducing hyperfine splittings (the mass difference between vector and pseudoscalar mesons). Notable discrepancies exist in the literature: the charmonium $2S$ splitting is significantly smaller than standard potential model predictions, and light quark systems (like $\rho-\pi$ ) exhibit splittings far exceeding standard expectations. Furthermore, the conventional contact spin-spin interaction term in the Fermi-Breit Hamiltonian diverges at the origin ( $r=0$ ), requiring regularization. There is a need for a consistent framework that addresses these singularities, incorporates relativistic corrections, and provides stable predictions for both ground states ( $1S$ ) and radially excited states ( $2S$ ) across heavy-light ( $D, B$ ) and heavy-heavy ( $c\bar{c}, b\bar{b}, B_c$ ) systems.

Methodology
The study employs a QCD-inspired potential model based on the Cornell ansatz:
$V(r) = -\frac{4\alpha_s}{3r} + br + c$
where the short-range Coulombic term represents one-gluon exchange and the long-range linear term represents confinement.

Key methodological components include:

Analytic Wavefunctions: The author utilizes the Dalgarno–Lewis perturbative method, treating the Coulombic part as the parent Hamiltonian and the linear confinement term as a perturbation. This yields an analytic wavefunction $\psi(r)$ that includes a factor $(r/a_0)^{-\epsilon}$ to account for relativistic corrections.
Regularization of the Contact Term: To resolve the singularity of the Fermi-Breit spin-spin interaction $V_{SS}(r) \propto \delta^{(3)}(\vec{r})$ , the author replaces the Dirac delta function with a Gaussian smearing function:
$\delta^{(3)}(\vec{r}) \rightarrow \frac{\sigma^3}{\pi^{3/2}} e^{-\sigma^2 r^2}$
Here, $\sigma$ is a phenomenological smearing parameter (set to 1.8 GeV) representing the effective size of constituent quarks. This ensures a finite, physically meaningful hyperfine integral.
Mass Spectrum Calculation: The physical mass is calculated as the sum of the spin-independent central mass and spin-dependent corrections (spin-orbit, tensor, and hyperfine terms). For $S$ -wave ground states, the mass difference is dominated by the hyperfine term $\Delta E_{HF}$ .
Parameter Space: The model uses fixed constituent quark masses ( $m_u=m_d=0.33$ GeV, $m_s=0.50$ GeV, $m_c=1.50$ GeV, $m_b=4.80$ GeV) and an effective strong coupling constant $\alpha_s = 0.36$ . The sensitivity of results to variations in $\alpha_s$ (within 0.34–0.38) is explicitly tested.

Key Contributions

Extension to Radial Excitations: Unlike many previous studies focusing solely on ground states, this work explicitly calculates hyperfine splittings for the first radially excited states ( $2S$ ) for $D$ , $D_s$ , $B$ , and $B_s$ mesons.
Unified Framework for Unequal Masses: The model successfully applies the same parameter set (effective $\alpha_s$ , string tension $b$ , and smearing width $\sigma$ ) to both heavy-light ( $D, B$ ) and heavy-heavy ( $c\bar{c}, b\bar{b}$ ) systems, demonstrating the stability of the parameter space across different mass asymmetries.
Analytic Approach with Regularization: The combination of the Dalgarno–Lewis analytic wavefunction with Gaussian smearing provides a compact, testable framework that avoids numerical divergences while maintaining physical consistency at short distances.

Results

Ground State ( $1S$ ) Accuracy: The model reproduces experimental ground-state meson masses with deviations typically below 3 MeV for heavy-light systems. The calculated hyperfine splittings for $1S$ states (e.g., $D^*-D \approx 140$ MeV, $B^*-B \approx 45$ MeV, $\Upsilon(1S)-\eta_b(1S) \approx 60$ MeV) align closely with experimental values (PDG 2024) and other theoretical models (e.g., Godfrey-Isgur).
Radial Excitation ( $2S$ ) Predictions: The model predicts reduced hyperfine splittings for $2S$ states compared to $1S$ states, consistent with the expectation of lower wavefunction density at the origin. Predicted values include $\sim 70$ MeV for $D(2S)$ , $\sim 36$ MeV for $B(2S)$ , and $\sim 57$ MeV for charmonium $2S$ .
$B_c$ System: Phenomenological estimates for the $B_c$ system are provided, showing consistency with the $1/(m_q m_{\bar{q}})$ scaling law. Specifically, the $B_c$ splittings are obtained through Fermi–Breit scaling, with the $B_c(2S)$ splitting estimated at 37 MeV.
Sensitivity Analysis: Variations in the effective coupling $\alpha_s$ within the range 0.34–0.38 result in changes to the hyperfine splitting of less than 3%, indicating the robustness of the results against small coupling fluctuations.

Significance and Claims
The paper claims to provide an accurate and unified description of meson hyperfine structure that bridges the gap between ground-state phenomenology and excited-state predictions. By incorporating Gaussian smearing, the model resolves the short-distance singularity inherent in the contact interaction without arbitrary parameterization. The author asserts that their framework offers "compact, testable benchmarks" for upcoming spectroscopy measurements at facilities such as LHCb and Belle II. The inclusion of $2S$ states and the $B_c$ system serves to validate the stability of the potential parameter space for unequal-mass systems and radial excitations. The study concludes that while perturbative treatment of the linear term limits accuracy for higher excitations, the current approach successfully captures the essential dynamics of meson binding for the $1S$ and $2S$ sectors.

Technical Summary: Hyperfine Structure of B and D Mesons in a QCD-Inspired Potential Model

More like this