A unified study of the \(B_c\) meson: from spectrum and form factors to weak and radiative decays

Abstract: The $B_c$ meson constitutes a unique system for investigating heavy-hadron dynamics, since it exhibits quarkonium-like bound-state structure while decaying predominantly through weak interactions. In this work, we present a unified study of $B_c$-meson spectroscopy, decay constants, weak decays, radiative transitions, and Regge trajectories within a single framework. The mass spectrum is computed in a screened potential model including relativistic kinetic-energy corrections and spin-dependent interactions, from which we obtain both spin-averaged and spin-resolved states together with the corresponding pseudoscalar and vector decay constants. Using these spectroscopic inputs, we then analyze weak decays within a mass-updated three-point QCD sum-rule framework for transitions to $S$-, $P$-, and $D$-wave charmonium states, as well as to final states containing charmonium and $D<sup>{(*)}_{(s)}$ mesons. In this part of the analysis, the hadronic thresholds, Borel-window prescriptions, Lorentz decompositions, and decay-width expressions of the underlying sum-rule formulations are retained, while the heavy-quark masses are updated to $m_c=1.48~\mathrm{GeV}$ and $m_b=4.90~\mathrm{GeV}$, leading to a controlled refit of the overall normalization rather than a full rederivation of all perturbative and condensate contributions. We further investigate purely leptonic decay widths and radiative $E1$ and $M1$ transitions, and examine the Regge behavior of the resulting spectrum. The present study therefore provides a coherent description of the $B_c$ meson in which the spectroscopic wave functions determine the short-distance couplings that enter both weak and electromagnetic observables, while the Regge analysis serves as a complementary global consistency test of the same dynamical picture.

• The paper introduces a single-parameter, screened-potential framework capturing the Bc meson’s spectrum, precise mass splittings, and state-dependent wavefunctions.
• It computes decay constants and weak transition form factors via QCD three-point sum rules, achieving results in line with lattice QCD and experimental benchmarks.
• Radiative, nonleptonic, and semileptonic decays are analyzed, revealing clear hierarchical channels and highlighting areas for improved QCD corrections.

Unified Heavy-Quarkonium Description of the $B_c$ Meson: Spectrum, Form Factors, and Decay Phenomenology

Framework for $B_c$ Meson Bound-State Dynamics

The study presents a comprehensive phenomenological analysis of the $B_c$ meson, modeling it as a heavy-heavy bound state within a semi-relativistic, screened-potential approach. The Hamiltonian formulation incorporates the full relativistic kinetic energy expanded in powers of momentum up to $\mathcal{O}(\mathbf{p}^{10})$ and a central potential, combining a short-distance Coulombic one-gluon-exchange term and a long-range linearly rising, but exponentially screened, confining interaction. Thereby, both spin-averaged masses and fine/hyperfine splittings are described through a consistent treatment incorporating spin-dependent interactions and relativistic corrections.

This single-parameter framework is variationally solved for all orbital and radial excitations. The resulting state-dependent wavefunctions underpin the calculation of spectroscopic masses, decay constants, overlap-sensitive short-distance observables, and wavefunction values at the origin. This approach, while reproducing quarkonium-like features for low-lying states, rationalizes softening at higher excitations, an essential aspect for interpreting Regge phenomenology and upper-level spectroscopy.

Mass Spectrum and Spin Structure

The computed spectrum reproduces ground-state and spin-averaged $B_c$ masses to within a few MeV of lattice-QCD and experimental results, establishing compatibility with existing benchmarks. Excited states, including $S$-, $P$-, $D$-, and $F$-wave levels, are consistently placed against quark model and relativistic potential model predictions (Table references in the paper). The incorporation of screening is especially evident in the moderate compression of levels for higher excitations, avoiding the overestimation of the mass gap observed in unscreened potentials.

Resolved fine and hyperfine splittings reflect the explicit mass asymmetry of the constituent $b$ and $B_c$0 quarks—contrasting with the near-degenerate pattern of equal-mass quarkonia—exposing a richer pattern of operator mixings. Parents and daughter Regge trajectories, both in $B_c$1 and $B_c$2 planes, demonstrate quasi-linear behavior for higher excitations, with deviations near the ground-state attributable to short-range effects.

Figure 1: Regge trajectories of the $B_c$3 meson in the $B_c$4 plane for natural and unnatural parity states; open and filled symbols distinguish theoretical predictions from observed levels.

Figure 2: Linear $B_c$5 Regge trajectories for spin-averaged and flavor-resolved $B_c$6 states, reflecting uniform radial excitation structure.

Decay Constants and Short-Distance Couplings

Wavefunctions at the origin are exploited, via a relativistically improved Van Royen-Weisskopf relation, to extract leptonic decay constants for both pseudoscalar and vector states. For the ground state, $B_c$7 converges with lattice and sum-rule estimates, $B_c$8–$B_c$9, while radial excitations display the expected monotonic suppression due to increased spatial extension, a direct prediction of the screened potential.

This decay-constant normalization coherently sets the scale for all short-distance observable calculations: weak-annihilation widths, purely leptonic branching ratios, and weak/EM transitions.

Weak Transition Form Factors and Semileptonic Decays

Transition amplitudes—pivotal for predicting exclusive weak decays—are calculated using QCD three-point sum rules. These are numerically evaluated in the spacelike domain and then analytically continued to the physical kinematic region via a truncated $B_c$0-series expansion. The renormalization procedure, including an explicit update of quark masses and normalization constants, is sector-specific to match the spectrum and decay framework.

A robust channel hierarchy emerges: $B_c$1-wave transitions $B_c$2, $B_c$3 dominate, with normalization sensitive to Coulomb and mass corrections. The analysis determines:

• Semileptonic widths: In the corrected scenario, $B_c$4 widths are the largest, followed by $B_c$5; $B_c$6- and $B_c$7-wave modes are strongly suppressed but hierarchically ordered (see Tables and Figures provided).
• LFU ratios: The results yield $B_c$8 and $B_c$9, compatible with the Standard Model and lattice QCD, and well below the original LHCb central value. No significant lepton flavor universality violation is observed given the uncertainties.

Figure 3: Differential decay widths in the $\mathcal{O}(\mathbf{p}^{10})$0-wave sector for bare and Coulomb-corrected calculations, showing the predominance of $\mathcal{O}(\mathbf{p}^{10})$1 over $\mathcal{O}(\mathbf{p}^{10})$2 transitions.

Figure 4: Bare and corrected $\mathcal{O}(\mathbf{p}^{10})$3 form factors $\mathcal{O}(\mathbf{p}^{10})$4 and $\mathcal{O}(\mathbf{p}^{10})$5; corrections systematically elevate normalization and improve theory/data consistency.

Figure 5: $\mathcal{O}(\mathbf{p}^{10})$6 form factors $\mathcal{O}(\mathbf{p}^{10})$7, $\mathcal{O}(\mathbf{p}^{10})$8, $\mathcal{O}(\mathbf{p}^{10})$9, $B_c$0 for both calculation variants, controlling distribution shapes for differential widths and angular observables.

For $B_c$1-wave final states, $B_c$2-rich channels are singled out as dominant over $B_c$3, a robust claim supported by sum-rule and Bethe-Salpeter comparanda. $B_c$4-wave transitions, still accessible only theoretically, present a consistent multiplet ordering with systematic suppression.

Nonleptonic Two-Body Modes

Factorizable nonleptonic decays involving charmonia and $B_c$5 are computed using the same form factors and modern decay constants. The flavor and spin hierarchy is unambiguous: decays to $B_c$6 are heavily favored over non-strange $B_c$7 via CKM and decay constant enhancement. However, the predicted normalization for leading channels is systematically lower than measured $B_c$8 branching fractions from LHCb/ATLAS, indicating a need for improved non-factorizable QCD effects or higher-order corrections.

Purely Leptonic and Inclusive Widths

Purely leptonic widths for $B_c$9 and $S$0 are computed, accounting for weak-annihilation and all-spectator processes. The predicted inclusive lifetime is on the longer side ($S$1 ps) relative to experimental determinations and recent OPE-based Standard Model calculations, signifying an underestimation of total width—potentially due to missing higher-order effects or unaccounted nonperturbative contributions.

Radiative Transitions

$S$2 radiative transitions dominate the EM sector, with computed widths matching the established magnitude ($S$3 keV)), validating the orbital structure of the model wavefunctions. $S$4 transitions, strictly allowed and hindered, show the anticipated suppression and are consistent with the spread among quark models.

Figure 6: $S$5-wave sector differential widths, highlighting the dominance of particular tensor and axial-vector channels after Coulomb correction.

Figure 7: $S$6 form factors $S$7 and $S$8, dictating exclusive width and $S$9 distributions for these suppressed final states.

Figure 8: $P$0, $P$1 form factors, controlling angular and polarization signatures in the $P$2-wave sector.

Global Phenomenology and Future Prospects

This work provides a unified and internally consistent phenomenology for the $P$3 system, where the screened-potential spectrum and state-dependent wavefunctions act as the organizing principle for all decay, radiative, and transition observables. The framework produces key results:

• Quantitative agreement with lattice and PDG for low-lying spectra.
• Channel-hierarchical weak decay predictions: S-wave semileptonic modes are dominant and numerically robust; P- and D-wave hierarchies align with recent sum-rule and lattice analyses.
• Strong numerical enhancement from Coulomb corrections, especially in pseudoscalar-mediated transitions.
• Regge trajectories manifesting approximately linear behavior—and hence consistent with the softening of the potential at high excitation.
• Discrepancies in nonleptonic decay normalization and inclusive lifetime, flagging the need for improved treatments of short-distance coefficients, non-factorizable corrections, and higher-order weak-interaction effects.

Implications for future research in AI-driven hadronic physics include:

• Enhanced modeling of heavy-heavy effective interactions with further systematic inclusion of coupled-channel and QCD corrections.
• Quantitative improvement in nonleptonic two-body modes through explicit QCD factorization and inclusion of power corrections.
• Utilization of Regge-trajectory fits as global consistency metrics for model-building.
• Provision of detailed, channel-differentiated observables as cross-validation targets for forthcoming LHC and Belle II measurements of excited $P$4 states and rare decay modes.

Conclusion

The screened-potential framework, incorporating relativistic and spin-dependent effects, provides a technically coherent and numerically viable description of the $P$5 meson from its spectrum up to its weak, radiative, and annihilation-induced decays. While the approach robustly captures the main features and hierarchies of the $P$6 system, persistent discrepancies in total width and selected exclusive modes highlight the importance of dynamic QCD effects and provide avenues for further theoretical refinement. As experimental data on excited states and rare transitions improves, frameworks such as this will underpin the systematic confrontation of strong- and weak-interaction dynamics in heavy-heavy mesons.