Spectroscopic parameters of $B_c$ meson

Technical Summary: Spectroscopic Parameters of the $B_c$ Meson

Problem and Motivation
The $B_c$ meson, a bound state of a bottom quark and a charm antiquark ($c\bar{b}$), occupies a unique position in heavy-quark spectroscopy. Unlike charmonium ($c\bar{c}$) or bottomonium ($b\bar{b}$), the $B_c$ system lacks charge-conjugation symmetry due to its unequal constituent masses. This absence of $C$-symmetry leads to structurally distinct modifications in leptonic decay constants and radiative transition charges. While the ground state ($1S$) and its first radial excitation ($2S$) have been experimentally established, higher excitations (P-wave, D-wave) and the vector ground state ($B_c^*$) remain largely unobserved or unconfirmed. Existing theoretical predictions for decay observables, such as the pseudoscalar decay constant $f_{B_c}$ and radiative widths, exhibit significant model dependence, with values for $f_{B_c}$ spanning 383–580 MeV and radiative widths varying by orders of magnitude. This sensitivity arises because these observables probe the short-distance Coulombic region of the interaction, where treatments of relativistic corrections, running coupling, and constituent masses differ. Furthermore, the string tension parameter in the confining potential is often adopted from flavor-symmetric systems ($c\bar{c}$ or $b\bar{b}$) rather than determined specifically for the $c\bar{b}$ sector.

Methodology
The authors investigate the $B_c$ meson within a unified nonrelativistic quark model based on the Cornell potential, supplemented by Breit–Fermi spin-dependent interactions. The effective Hamiltonian includes the Cornell interaction (Coulombic one-gluon exchange plus linear confinement), a spin-spin hyperfine interaction ($V_{SS}$), and a spin-orbit interaction ($V_{LS}$).

Key methodological features include:

• Nonperturbative Treatment: Both $V_{SS}$ and $V_{LS}$ are incorporated directly into the radial Schrödinger equation rather than treated perturbatively. The contact term of the hyperfine interaction is regularized using a Gaussian smearing function. This ensures that singlet and triplet states acquire distinct radial wave functions from the outset, which is critical for calculating decay constants and radiative widths.
• Parameter Determination: The model parameters (quark masses $m_b, m_c$, strong coupling $\alpha_s$, string tension $b$, and smearing width $\sigma$) are taken from Ref. [55] and determined specifically from $B_c$ spectroscopic data (pseudoscalar masses), rather than being fitted to flavor-symmetric systems.
• Calculations: The study computes the mass spectrum for radially and orbitally excited states up to $n=5$ in S-, P-, and D-wave sectors. It calculates pseudoscalar and vector decay constants ($f_P, f_V$) including one-loop QCD radiative corrections. Radiative transition widths for Electric Dipole (E1) and Magnetic Dipole (M1) processes are computed using full recoil-corrected overlap integrals, avoiding the long-wavelength approximation. Finally, radial Regge trajectories are constructed in the $(n, M^2)$ plane to test the global consistency of the spectrum.

Key Contributions and Results

1. Mass Spectrum and Fine Structure:

   • The ground state mass $M(1^1S_0) = 6272$ MeV and the first radial excitation $M(2^1S_0) = 6873$ MeV agree with experimental measurements within 2–3 MeV.
   • The model predicts an inverted fine-structure ordering for P- and D-wave multiplets: $M(^3P_0) > M(^3P_1) > M(^3P_2)$ and $M(^3D_1) > M(^3D_2) > M(^3D_3)$. This inversion results from the Thomas-precession term in the spin-orbit interaction slightly dominating the one-gluon-exchange (OGE) term at the relevant interquark separation ($\approx 0.55$ fm).
   • The predicted hyperfine splitting for the ground state, $\Delta M(1S) \approx 3$ MeV, is significantly smaller than the theoretical consensus (50–70 MeV) and lattice QCD values. The authors attribute this to the fact that the model parameters were fixed using pseudoscalar masses, leaving the hyperfine scale unconstrained.

2. Decay Constants:

   • The predicted pseudoscalar decay constant is $f_{B_c} = 580$ MeV, which is higher than most nonrelativistic potential model predictions (433–530 MeV) but aligns closely with relativistic calculations.
   • The vector decay constant is $f_{B_c^*} = 546$ MeV, yielding a ratio $f_{B_c^*}/f_{B_c} \approx 0.94$. This suppression is identified as a direct one-loop QCD consequence of the unequal-mass system.
   • Using these constants, the leptonic branching fraction is predicted as $\mathcal{B}(B_c \to \tau\nu_\tau) \approx 2\%$.

3. Radiative Transitions:

   • E1 Transitions: The transition pattern is governed primarily by radial wave-function overlaps rather than phase space. Transitions between adjacent radial shells (e.g., $nS \to (n-1)P$) dominate over cross-shell transitions, even when the latter involve significantly higher photon energies. The dominant $1D \to 1P$ transitions (50–87 keV) suggest the $1D \xrightarrow{\gamma} 1P \xrightarrow{\gamma} 1S$ cascade as the primary experimental pathway to observe D-wave states.
   • M1 Transitions: The allowed M1 widths (e.g., $1^3S_1 \to 1^1S_0$) are strongly suppressed ($< 0.01$ eV) due to the small predicted hyperfine splitting ($E_\gamma \sim 1-3$ MeV). In contrast, hindered M1 transitions (e.g., $2^1S_0 \to 1^3S_1$) are enhanced by orders of magnitude (up to $\sim 434$ eV) due to the breaking of orthogonality by finite photon momentum and the distortion of singlet/triplet wave functions by the hyperfine interaction.

4. Regge Trajectories:

   • The calculated mass spectra form nearly linear radial Regge trajectories ($R^2 > 0.989$) across all ten spin-parity families.
   • The slopes ($\mu^2$) decrease monotonically with increasing orbital angular momentum ($L$), from $\approx 5.83$ GeV$^2$ (S-wave) to $\approx 4.73$–$4.84$GeV$^2$ (D-wave). This reflects the increasing dominance of the linear confining term at larger interquark separations sampled by higher $L$ states.

Significance and Claims
The paper claims that its results provide a self-consistent and quantitatively predictive description of the $B_c$ meson spectrum based on parameters derived specifically for the $c\bar{b}$ sector. The authors identify two critical areas for future experimental validation:

1. The $B_c^*$ Mass: A direct measurement of the $B_c^*(1S)$ mass is described as the "single most important missing input" in $B_c$ physics. Such a measurement would fix the hyperfine splitting, independently constrain the smearing parameter $\sigma$, and resolve the ambiguity regarding the allowed M1 widths.
2. Fine-Structure Ordering: Experimental confirmation of the P-wave multiplet members would test the predicted inverted ordering, thereby validating the competition between the Thomas-precession and OGE contributions to the spin-orbit interaction.

The authors conclude that while the current model successfully reproduces known S-wave states and offers specific predictions for unobserved states, the precision of these predictions regarding spin-dependent dynamics is contingent upon future high-luminosity measurements at the LHC and future $e^+e^-$ facilities. The study emphasizes that the small predicted hyperfine splitting is a structural consequence of the adopted parameter set, not necessarily a physical prediction, and highlights the necessity of direct mass measurements to distinguish between competing theoretical frameworks.