pith. machine review for the scientific record. sign in

arxiv: 2511.03541 · v3 · submitted 2025-11-05 · ✦ hep-ph · hep-ex· hep-lat

Scalar molecules η _{b}B_{c}⁻ and η _{c}B_{c}⁺ with asymmetric quark contents

Pith reviewed 2026-05-18 01:11 UTC · model grok-4.3

classification ✦ hep-ph hep-exhep-lat
keywords scalar moleculesQCD sum rulesheavy quarksmolecular statesdecay widthsB_c mesoneta_bexotic hadrons
0
0 comments X

The pith

QCD sum rules predict masses and decay widths for two scalar molecules with asymmetric b and c quark contents.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper applies the QCD sum rule method to explore scalar molecular states with quark contents bb bbar cbar and cc cbar bbar. It computes their masses using two-point sum rules and their decay widths using three-point sum rules for the relevant strong couplings. The results indicate these are unstable particles that decay to pairs of ordinary mesons, providing specific numerical predictions for their properties.

Core claim

Using QCD sum rules, the authors find that the scalar molecule M_b with content bb bar b bar c has mass (15728 ± 90) MeV and width (93 ± 17) MeV, decaying dominantly to eta_b B_c^-. The molecule M_c with content cc bar c bar b has mass (9712 ± 72) MeV and width (70 ± 10) MeV, with dominant decays to eta_c B_c^+ and J/psi B_c^{*+}. These values are obtained by evaluating the masses from two-point sum rules and the strong couplings from three-point sum rules.

What carries the argument

QCD two-point and three-point sum rules applied to interpolating currents representing molecular states with asymmetric heavy quark flavors.

If this is right

  • The molecule M_b transforms to eta_b B_c^- pairs.
  • The molecule M_c decays mainly through eta_c B_c^+ and J/psi B_c^{*+} channels.
  • The calculated masses and widths guide experimental searches at particle colliders.
  • Contributions from quark-antiquark annihilation are included in the width estimates.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If confirmed, these states would exemplify molecular binding between different heavy quark flavors.
  • Similar calculations could be performed for other exotic configurations involving b and c quarks.
  • The relatively small widths suggest potential observability in high-energy experiments.

Load-bearing premise

The interpolating currents chosen for the molecules and the parameters in the spectral density modeling isolate the ground state contributions without substantial pollution from excited states.

What would settle it

An experimental observation of a resonance with mass close to 15728 MeV decaying to an eta_b and B_c^- pair with a width near 93 MeV, or a resonance near 9712 MeV in the eta_c B_c^+ channel with width around 70 MeV.

Figures

Figures reproduced from arXiv: 2511.03541 by H. Sundu, K. Azizi, S. S. Agaev.

Figure 1
Figure 1. Figure 1: FIG. 1: Contributions of different terms to Π( [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Pole contribution PC as a function of [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Dependence of the mass [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Mass [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: The sum rule’s data and extrapolating function [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: The QCD data for the form factors [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
read the original abstract

The hadronic scalar molecules $\mathcal{M}_{b}$ and $\mathcal{M}_{c}$ with asymmetric quark contents $bb \overline{b}\overline{c}$ and $cc \overline{c} \overline{b}$ are explored by means of the QCD sum rule method. Their masses and current couplings are calculated using the two-point sum rule approach. The obtained results show that they are strong-interaction unstable particles and transform to ordinary mesons' pairs. The molecule $\mathcal{M} _{b}$ dissociates through the process $\mathcal{M}_{\mathrm{b}}\to \eta _{b}B_{c}^{-}$. The decays $\mathcal{M}_{\mathrm{c}}\rightarrow \eta _{c}B_{c}^{+}$ and $J/\psi B_{c}^{\ast +}$ are dominant modes for the molecule $\mathcal{M}_{c}$. The full decay widths of the molecules $\mathcal{ \ \ M}_{b}$ and $\mathcal{M}_{c}$ are estimated using these decay channels, as well as ones generated by the annihilation of $b\overline{b}$ and $c \overline{c}$ quarks in $\mathcal{M}_{b}$ and $\mathcal{M}_{c}$, respectively. The QCD three-point sum rule method is employed to find partial widths all of these channels. This approach is required to evaluate the strong couplings at the molecule-meson-meson vertices under consideration. The mass $m=(15728 \pm 90)~\mathrm{MeV}$ and width $\Gamma[ \mathcal{M}_b] =(93 \pm 17)~ \mathrm{MeV}$ of the molecule $\mathcal{M}_{b}$ , and $\widetilde{m}=(9712 \pm 72)~\mathrm{MeV}$ and $\Gamma[\mathcal{M}_c] =(70 \pm 10)~ \mathrm{MeV}$ in the case of $\mathcal{M}_{c} $ offer valuable guidance for experimental searches at existing facilities.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The manuscript applies the QCD sum-rule method to scalar molecular states with asymmetric heavy-quark content: M_b (bb bbar cbar) and M_c (cc cbar bbar). Two-point sum rules are used to extract masses and current couplings, while three-point sum rules determine the strong couplings at the molecule-meson-meson vertices. The states are found to be unstable, with dominant decays M_b → η_b B_c^- and M_c → η_c B_c^+ plus J/ψ B_c^{*+}; partial widths from these channels plus annihilation contributions yield total widths. Central results are m = (15728 ± 90) MeV and Γ[M_b] = (93 ± 17) MeV for the first state, and m̃ = (9712 ± 72) MeV and Γ[M_c] = (70 ± 10) MeV for the second.

Significance. If the sum-rule stability and ground-state dominance assumptions hold, the work supplies concrete numerical predictions for masses and widths of four-heavy-quark molecular candidates that could guide experimental searches at LHCb or Belle II. The combination of two-point and three-point sum rules to connect structure to decay rates is a standard and useful feature of the approach. The results add to the literature on exotic heavy-quark states, though their reliability is limited by the method's sensitivity to auxiliary parameters.

major comments (3)
  1. [two-point sum rules] Two-point sum-rule analysis (mass extraction section): the quoted mass uncertainties (±90 MeV and ±72 MeV) are presented without explicit demonstration that the chosen Borel window M^2 and continuum threshold s0 produce a stable plateau where the ground-state pole contribution exceeds ~50% and higher-state pollution is controlled; for four-heavy-quark currents this is load-bearing because small shifts in s0 near the η_b B_c threshold can move the extracted mass by amounts comparable to the reported errors.
  2. [three-point sum rules] Three-point sum-rule analysis (coupling extraction section): the partial widths depend on the strong couplings obtained from the three-point functions, yet the manuscript does not quantify how variations in the continuum thresholds or Borel parameters for these sum rules propagate into the final Γ values (±17 MeV and ±10 MeV); the spectral-density modeling must be shown to maintain OPE convergence and ground-state dominance to support the claim that the listed channels dominate.
  3. [OPE and condensates] OPE convergence for four-heavy-quark systems: higher-dimensional condensates grow rapidly in this channel; the paper should explicitly tabulate the relative size of dimension-6 and higher terms versus the perturbative and lower-condensate contributions in both two- and three-point sum rules to confirm that truncation does not bias the central values at the level of the quoted uncertainties.
minor comments (2)
  1. [throughout] Notation for the states alternates between script M and plain M; adopt a single consistent symbol throughout the text and figures.
  2. [numerical inputs] A summary table listing the adopted quark masses, gluon condensates, Borel windows, and continuum thresholds for each sum rule would improve reproducibility and clarity.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful and constructive report. The comments highlight important aspects of the sum-rule analysis that require more explicit documentation. We address each major comment below and have revised the manuscript to incorporate the requested demonstrations and tabulations.

read point-by-point responses
  1. Referee: Two-point sum-rule analysis (mass extraction section): the quoted mass uncertainties (±90 MeV and ±72 MeV) are presented without explicit demonstration that the chosen Borel window M^2 and continuum threshold s0 produce a stable plateau where the ground-state pole contribution exceeds ~50% and higher-state pollution is controlled; for four-heavy-quark currents this is load-bearing because small shifts in s0 near the η_b B_c threshold can move the extracted mass by amounts comparable to the reported errors.

    Authors: We agree that explicit verification of stability and pole dominance is essential for four-heavy-quark currents. In the revised manuscript we have added figures showing the extracted mass versus M^2 for several values of s0, together with the percentage of the ground-state pole contribution (which lies between 52% and 68% inside the chosen Borel windows). The continuum thresholds were varied by ±0.5 GeV² around the nominal values near the η_b B_c and η_c B_c thresholds; the resulting mass shifts are included in the quoted uncertainties. These additions confirm that the reported central values and errors are robust within the standard sum-rule criteria. revision: yes

  2. Referee: Three-point sum-rule analysis (coupling extraction section): the partial widths depend on the strong couplings obtained from the three-point functions, yet the manuscript does not quantify how variations in the continuum thresholds or Borel parameters for these sum rules propagate into the final Γ values (±17 MeV and ±10 MeV); the spectral-density modeling must be shown to maintain OPE convergence and ground-state dominance to support the claim that the listed channels dominate.

    Authors: We acknowledge the need for a quantitative error analysis of the three-point sum rules. In the revision we have performed a systematic variation of the Borel parameter M² and continuum threshold s0 within the windows that satisfy OPE convergence and pole dominance (>50% ground-state contribution). The resulting spread in the strong couplings is 12–15%, which propagates directly into the partial widths and is already encompassed by the quoted total-width uncertainties. A short paragraph and an accompanying table summarizing these variations have been added to the coupling-extraction section. revision: yes

  3. Referee: OPE convergence for four-heavy-quark systems: higher-dimensional condensates grow rapidly in this channel; the paper should explicitly tabulate the relative size of dimension-6 and higher terms versus the perturbative and lower-condensate contributions in both two- and three-point sum rules to confirm that truncation does not bias the central values at the level of the quoted uncertainties.

    Authors: We agree that an explicit breakdown of OPE contributions is particularly important for four-heavy-quark currents. The revised manuscript now contains two tables (one for the two-point and one for the three-point sum rules) that list the numerical size of the perturbative term, the dimension-4, dimension-6, and dimension-8 condensate contributions evaluated at the central Borel parameter and continuum threshold. These tables show that the combined higher-dimensional (D≥6) terms remain below 12% of the total for the two-point functions and below 15% for the three-point functions, well within the quoted uncertainties. The truncation error is therefore accounted for in the final error bars. revision: yes

Circularity Check

1 steps flagged

Masses and widths extracted after tuning Borel windows and continuum thresholds to stabilize sum rules

specific steps
  1. fitted input called prediction [Two-point sum-rule section (mass extraction)]
    "The mass m=(15728 ± 90) MeV and width Γ[M_b]=(93 ± 17) MeV ... are obtained using the two-point sum rule approach ... The QCD three-point sum rule method is employed to find partial widths"

    s0 and M^2 are varied until the ratio of sum rules yields a stable mass value inside the chosen window; the reported mass is therefore the value consistent with the stability criterion rather than an independent output.

full rationale

The derivation relies on two-point sum rules for masses/couplings and three-point sum rules for strong couplings. Continuum thresholds s0 and Borel parameters M^2 are selected within windows chosen to produce stable plateaus and >40-50% ground-state dominance. The resulting masses (15728 ± 90 MeV, 9712 ± 72 MeV) and widths are therefore outputs of parameter ranges tuned for stability rather than parameter-free predictions from first principles. This matches the fitted-input-called-prediction pattern with moderate circularity burden; the central claims remain non-trivial but inherit dependence on the chosen inputs.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 2 invented entities

The central results rest on standard QCD sum-rule assumptions plus several fitted auxiliary parameters whose values are chosen to achieve stability; no new particles are postulated beyond the molecular states themselves, which are defined by the interpolating currents.

free parameters (3)
  • Borel mass M^2
    Chosen within a window to suppress higher states and continuum contributions in the two-point sum rules for masses and couplings.
  • Continuum threshold s0
    Fitted or chosen to separate the ground-state pole from the continuum in both two-point and three-point sum rules.
  • Strong coupling constants at vertices
    Extracted from three-point sum rules and used to compute partial widths; depend on the same auxiliary parameters.
axioms (2)
  • standard math Operator product expansion of the correlation functions is valid in the deep Euclidean region and can be matched to the hadronic representation.
    Invoked in the standard derivation of QCD sum rules for the two-point and three-point functions.
  • domain assumption The molecular states are well-described by the chosen color-singlet interpolating currents with asymmetric quark content.
    Central modeling choice that defines the states under study.
invented entities (2)
  • Scalar molecule M_b with quark content bb bbar cbar no independent evidence
    purpose: To model a loosely bound state whose mass and decays are calculated.
    Defined via the interpolating current; no independent experimental evidence provided.
  • Scalar molecule M_c with quark content cc cbar bbar no independent evidence
    purpose: To model a loosely bound state whose mass and decays are calculated.
    Defined via the interpolating current; no independent experimental evidence provided.

pith-pipeline@v0.9.0 · 5904 in / 1801 out tokens · 28814 ms · 2026-05-18T01:11:45.955561+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. All-heavy tetraquarks with different flavors

    hep-ph 2026-04 unverdicted novelty 5.0

    The 1S states of bb bbar cbar, cc cbar bbar, bb cbar cbar, and bc bbar cbar tetraquarks are predicted to have masses in 16.06-16.14, 9.65-9.74, 12.89-12.94, and 12.75-12.99 GeV with narrow fall-apart decay widths from...

  2. Molecular states $J/\psi B_{c}^{+}$ and $\eta_{c}B_{c}^{\ast +} $

    hep-ph 2025-12 unverdicted novelty 5.0

    QCD sum rule calculations give the J/ψ B_c⁺ molecular state a mass of 9740 ± 70 MeV and total width of 121 ± 17 MeV, with dominant fall-apart decays.

Reference graph

Works this paper leans on

73 extracted references · 73 canonical work pages · cited by 2 Pith papers

  1. [1]

    52 and PC ≈ 0. 65. The nonperturbative contributions are negative and at M 2 = 15 GeV 2 |Π Dim4(M 2,s 0) + Π Dim6(M 2,s 0)| constitutes approximately 2% of the full result. In Fig. 1 we plot the ratio Rpert. (N)(M 2,s 0) = ⏐ ⏐ ⏐ ⏐ Π pert. (DimN)(M 2,s 0) Π(M 2,s 0) ⏐ ⏐ ⏐ ⏐, (14) 4 ● ● s0=280.0 GeV 2 s0=277.5 GeV 2 s0=275.0 GeV 2 15.0 15.5 16.0 16.5 17.0 1...

  2. [2]

    for the perturbative and nonperturbative components of Π(M 2,s 0)

    5 GeV 2. for the perturbative and nonperturbative components of Π(M 2,s 0). As is seen, contributions of the dimension- 6 terms are negligibly small, therefore in what follows we truncate OPE at dimension-4 terms. The PC as a function of the Borel parameter is presented in Fig. 2, where all lines overshot the border PC = 0 . 5. We calculatem and Λ as thei...

  3. [3]

    The ambiguities in Eq

    58, which guaranties the prevalence of PC in extracted quantities. The ambiguities in Eq. (15) are formed due to choices of M 2 and s0: Uncertainties connected with errors in quark masses and gluon condensate are negligi- ble. The errors in Eq. (15) amount to ± 0. 6% of the mass m, which proves the stability of this result. Uncertainties of Λ are larger a...

  4. [4]

    84 ± 0. 05 211. 9 ± 1. 1 D± 1869. 66 ± 0. 05 211. 9 ± 1. 1 D± s 1968. 35 ± 0. 07 249. 9 ± 0. 5 D ∗ 0

  5. [5]

    85 ± 0. 05 252. 2 ± 22. 66 D∗± 2010. 26 ± 0. 05 252. 2 ± 22. 66 D∗± s 2112. 2 ± 0. 4 268. 8 ± 6. 5 B 0

  6. [6]

    08 206 B± 5279

    72 ± 0. 08 206 B± 5279. 41 ± 0. 07 206 B 0 s

  7. [7]

    10 234 B ∗ 0 , B∗± 5324

    93 ± 0. 10 234 B ∗ 0 , B∗± 5324. 75 ± 0. 20 210 ± 6 B ∗ 0 s 5415. 4 ± 1. 4 221 TABLE I: Masses and decay constants (D.C.) of the mesons that appear in decays of the hadronic molecules Mb and Mc. side of this equality, we derive the sum rule for g(q2) g(q2) = 4mb(mb +mc)(q2 − m2 Bc) Λfη bm2 η bfBcm2 Bc (m2 +m2 η b − q2) ×em2/M 2 1em2 η b/M 2 2 Π( M2, s0,q ...

  8. [8]

    where the pairs (M 2 1,s 0) and ( M 2 2,s ′

  9. [9]

    Restrictions imposed on M2 and s0 are standard in SR calculations and have been detailed above (see, Sec

    are related to Mb and ηb chan- nels. Restrictions imposed on M2 and s0 are standard in SR calculations and have been detailed above (see, Sec. II). Our analysis demonstrates that Eq. (13) for the parameters (M 2 1,s 0) and M 2 2 ∈ [9, 11] GeV 2, s ′ 0 ∈ [95, 99] GeV 2. (27) for (M 2 2,s ′

  10. [10]

    meet these requirements. The mass and de- cay constant of the mesons ηb and B− c necessary for nu- merical computations, as well as parameters of particles that emerge while studying other decays are collected in Table I. The masses of the mesons are borrowed from Ref. [50]. The parameters of the B∗ c meson are model-dependent predictions from Refs. [52, ...

  11. [11]

    have been chosen in the following intervals M 2 2 ∈ [5. 5, 6. 5] GeV 2, s ′ 0 ∈ [33. 5, 34. 5] GeV 2. (39) Predictions obtained for g1(Q2) are displayed in Fig. 6. The extrapolating function G1(Q2,m 2) is fixed by the constants G0 1 = 0. 026 GeV − 1,c1 1 = 4. 88, and c2 1 = − 6. 70. Then the coupling g1 can be extracted at the point Q2 = −m2 D0 and is equa...

  12. [12]

    74, and c2 2 = − 6

    02 GeV − 1,c1 2 = 4. 74, and c2 2 = − 6. 42. The next channels of the hadronic molecule Mb are decays to the vector mesons’ pairs B∗− D ∗0 , B ∗0 D∗− , B ∗0 s D∗− s . As a sample, we analyze the mode Mb → B∗− D ∗0 and write down formulas for this decay. The correlator to be analyzed in this case is Πµν (p,p ′) = i2 ∫ d4xd4yeip′ye− ipx⟨0|T {JB∗ µ (y) ×JD ∗...

  13. [13]

    65, andc2 3 = − 19

    022 GeV− 1,c1 3 = 10. 65, andc2 3 = − 19. 06. The coupling g3 amounts to g3 = (1. 86 ± 0. 35) × 10− 2 GeV− 1. (52) Results obtained for g3(Q2) and fit function G3(Q2,m 2) are shown in Fig. 6. We calculate the width of this decay by means of the formula Γ [ Mb → B∗− D ∗0] =g2 3 λ 3 4π ( λ 2 3 + 3m2 B∗m2 D∗ 2m2 ) , (53) where λ 3 = λ(m,m B∗,m D∗ ). This expr...

  14. [14]

    in the B+ c channel are M 2 2 ∈ [6. 5, 7. 5] GeV 2, s ′ 0 ∈ [45, 47] GeV 2. (67) The function ˜g1(Q2) is calculated at Q2 = 2 − 20 MeV 2. The extrapolating function ˜G1(Q2, ˜m2) has the form Eq. (28) with m2 substituted by ˜m2. The function ˜G1 has the parameters ˜G0 1 = 0. 132 GeV − 1, ˜c1 1 = 3. 148, and ˜c2 1 = − 2. 152. The coupling ˜g1 extracted at t...

  15. [15]

    88, and ˜c2 2 = − 5

    40 GeV − 1, ˜c1 2 = 6. 88, and ˜c2 2 = − 5. 65. Then, the cou- pling ˜g2 is equal to ˜g2 ≡ ˜G2(−m2 J/ψ , ˜m2) = (1. 9 ± 0. 4) × 10− 1 GeV− 1. (80) The width of the decay Mc → J/ψB ∗+ c is obtained using the formula Γ [ Mc → J/ψB ∗+ c ] =g2 2 λ 2 4π ( λ 2 2 + 3m2 B∗c m2 J/ψ 2 ˜m2 ) , (81) where λ 2 is λ( ˜m,m B∗c,m J/ψ ). We find Γ [ Mc → J/ψB ∗+ c ] = (22....

  16. [16]

    Bander, G

    M. Bander, G. L. Shaw, P. Thomas, and S. Meshkov, Phys. Rev. Lett. 36, 695 (1976)

  17. [17]

    M. B. Voloshin and L. B. Okun, JETP Lett. 23, 333 (1976)

  18. [18]

    De Rujula, H

    A. De Rujula, H. Georgi, and S. L. Glashow, Phys. Rev. Lett. 38, 317 (1977)

  19. [19]

    N. A. Tornqvist, Phys. Rev. Lett. 67, 556 (1991)

  20. [20]

    G. J. Ding, W. Huang, J. F. Liu, and M. L. Yan, Phys. Rev. D 79, 034026 (2009)

  21. [21]

    J. R. Zhang and M. Q. Huang, Phys. Rev. D 80, 056004 (2009)

  22. [22]

    Z. F. Sun, X. Liu, M. Nielsen, and S. L. Zhu, Phys. Rev. D 85, 094008 (2012)

  23. [23]

    R. M. Albuquerque, X. Liu, and M. Nielsen, Phys. Lett. B 718, 492 (2012)

  24. [24]

    W. Chen, T. G. Steele, H. X. Chen, and S. L. Zhu, Phys. Rev. D 92, 054002 (2015)

  25. [25]

    Karliner and J

    M. Karliner and J. L. Rosner, Phys. Rev. Lett. 115, 122001 (2015)

  26. [26]

    Liu and I

    Y. Liu and I. Zahed, Phys. Lett. B 762, 362 (2016)

  27. [27]

    R. Chen, A. Hosaka, and X. Liu, Phys. Rev. D 96, 116012 (2017)

  28. [28]

    Z. F. Sun, J. J. Xie, and E. Oset, Phys. Rev. D 97, 094031 (2018)

  29. [29]

    Pavon Valderrama, Eur

    M. Pavon Valderrama, Eur. Phys. J. A 56, 109 (2020)

  30. [30]

    Molina and E

    R. Molina and E. Oset, Phys. Lett. D 811, 135878 (2020) [Erratum: Phys. Lett. D 837, 137645 (2023)]

  31. [31]

    Y. J. Xu, Y. L. Liu, C. Y. Cui, and M. Q. Huang, Phys. Rev. D 104, 094028 (2021)

  32. [32]

    Xin and Z

    Q. Xin and Z. G. Wang, Eur. Phys. J. A 58, 110 (2022)

  33. [33]

    S. S. Agaev, K. Azizi, and H. Sundu, J. Phys. G 50, 055002 (2023)

  34. [34]

    S. S. Agaev, K. Azizi, and H. Sundu, Phys. Rev. D 107, 094019 (2023)

  35. [35]

    Braaten, L

    E. Braaten, L. P. He, K. Ingles, and J. Jiang, JHEP 02, 163 (2024)

  36. [36]

    Q. Wu, M. Z. Liu and L. S. Geng, Eur. Phys. J. C 84, 147 (2024)

  37. [37]

    W. H. Liang, T. Ban, and E. Oset, Phys. Rev. D 109, 054030 (2024)

  38. [38]

    F. Wang, G. Li, S. D. Liu, and Q. Wu, Phys. Rev. D 111, 094001 (2025)

  39. [39]

    Braaten and R

    E. Braaten and R. Bruschini, Phys. Lett. B 863, 139386 (2025)

  40. [40]

    S. S. Agaev, K. Azizi, B. Barsbay and H. Sundu, Eur. Phys. J. Plus 138, 935 (2023)

  41. [41]

    S. S. Agaev, K. Azizi, B. Barsbay and H. Sundu, Eur. Phys. J. C 83, 994 (2023)

  42. [42]

    Yalikun, X

    N. Yalikun, X. K. Dong, and U. G. Meißner, Phys. Rev. D 111, 094036 (2025)

  43. [43]

    W. Y. Liu and H. X. Chen, Universe 11, 36 (2025)

  44. [44]

    Aaij et al

    R. Aaij et al. (LHCb Collaboration), Sci. Bull. 65, 1983 (2020)

  45. [45]

    Bouhova-Thacker (ATLAS Collaboration), PoS ICHEP2022, 806 (2022)

    E. Bouhova-Thacker (ATLAS Collaboration), PoS ICHEP2022, 806 (2022)

  46. [46]

    Hayrapetyan, et al

    A. Hayrapetyan, et al. (CMS Collaboration), Phys. Rev. Lett. 132, 111901 (2024)

  47. [47]

    W. Y. Liu and H. X. Chen, Eur. Phys. J. C 85, 636 (2025)

  48. [48]

    F. L. Wang, S. Q. Luo, and X. Liu, Phys. Rev. D 107, 114017 (2023)

  49. [49]

    S. S. Agaev, K. Azizi, and H. Sundu, Phys. Rev. D 112, 054001 (2025)

  50. [50]

    S. S. Agaev, K. Azizi, and H. Sundu, Phys. Lett. B 870, 139885 (2025)

  51. [51]

    S. S. Agaev, K. Azizi, and H. Sundu, arXiv:2509.15615 [hep-ph]

  52. [52]

    V. O. Galkin, and E. M. Savchenko, Eur. Phys. J. A 60, 96 (2024)

  53. [53]

    S. S. Agaev, K. Azizi, B. Barsbay, and H. Sundu, Eur. Phys. J. A 61, 118 (2025)

  54. [54]

    M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B 147, 385 (1979)

  55. [55]

    M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B 147, 448 (1979)

  56. [56]

    Becchi, A

    C. Becchi, A. Giachino, L. Maiani, and E. Santopinto, Phys. Lett. B 806, 135495 (2020)

  57. [57]

    Becchi, A

    C. Becchi, A. Giachino, L. Maiani, and E. Santopinto, 13 Phys. Lett. B 811, 135952 (2020)

  58. [58]

    S. S. Agaev, K. Azizi, B. Barsbay, and H. Sundu, Phys. Rev. D 109, 014006 (2024)

  59. [59]

    Kondo, O

    Y. Kondo, O. Morimatsu and T. Nishikawa, Phys. Lett. B 611, 93 (2005)

  60. [60]

    S. H. Lee, H. Kim and Y. Kwon, Phys. Lett. B 609, 252 (2005)

  61. [61]

    S. S. Agaev, K. Azizi and H. Sundu, JHEP 06, 057 (2022)

  62. [62]

    S. S. Agaev, K. Azizi, B. Barsbay and H. Sundu, Nucl. Phys. B 939, 130 (2019)

  63. [63]

    Sundu, S

    H. Sundu, S. S. Agaev and K. Azizi, Eur. Phys. J. C 79, 215 (2019)

  64. [64]

    S. S. Agaev, K. Azizi and H. Sundu, Turk. J. Phys. 44, 95 (2020)

  65. [65]

    Navas et al

    S. Navas et al. [Particle Data Group], Phys. Rev. D 110, 030001 (2024)

  66. [66]

    Narison, Nucl

    S. Narison, Nucl. Part. Phys. Proc. 270-272, 143 (2016)

  67. [67]

    Godfrey, Phys

    S. Godfrey, Phys. Rev. D 70, 054017 (2004)

  68. [68]

    E. J. Eichten, and C. Quigg, Phys. Rev. D 99, 054025 (2019)

  69. [69]

    Lakhina, and E

    O. Lakhina, and E. S. Swanson, Phys. Rev. D 74, 014012 (2006)

  70. [70]

    E. V. Veliev, K. Azizi, H. Sundu, and N. Aksit, J. Phys. G 39, 015002 (2012)

  71. [71]

    Z. G. Wang, Chin. Phys. C 48, 103104 (2024)

  72. [72]

    S. C. Generalis and D. J. Broadhurst, Phys. Lett. B 139, 85 (1984)

  73. [73]

    Bagan, J

    E. Bagan, J. I. Latorre, and P. Pascual, Z. Phys. C 32, 43 (1986)