Pith. sign in

REVIEW 4 major objections 6 minor 39 references

The T_cbar{s} state is generated by the off-diagonal DK–D_sπ coupling, not by diagonal interactions.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · deepseek-v4-flash

2026-08-04 12:51 UTC pith:5SS7QLRX

load-bearing objection Promising unified mechanism for T_cbar{s} with a testable D_s1(2536) prediction, but the central pole's existence rests on a channel-dependent cutoff that needs a robustness check. the 4 major comments →

arxiv 2510.01564 v2 pith:5SS7QLRX submitted 2025-10-02 hep-ph hep-ex

A Paradigm for the Coupled-Channel Origin of Resonances: the Exotic T_{cbar{s}} in D_(s1)(2460/2536)to D_sππ

classification hep-ph hep-ex
keywords T_cbar{s}tetraquarkcoupled-channeloff-diagonal potentialD_s1(2460)D_s1(2536)triangle singularityexotic hadron
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

This letter aims to show that the exotic T_cbar{s} tetraquark candidate seen in D_s1(2460)+ → D_s+ π+ π− is not a kinematic triangle-singularity artifact but a genuine resonance generated by off-diagonal coupling between the DK and D_sπ channels. The same framework, combining triangle loops with DK–D_sπ rescattering, reproduces the two-peak structure in D_s1(2460) decay and predicts a single broad peak in D_s1(2536) decay. That contrast is traced to the S-wave versus D-wave character of the two D_s1 states' couplings to D*K, making the D_s1(2536) lineshape a decisive experimental test. A sympathetic reader would care because it turns a single exotic-state observation into a template for how hadronic structure and decay encode the same coupled-channel dynamics.

Core claim

On the paper's own terms, the T_cbar{s} resonance is dynamically generated by the off-diagonal coupled-channel potential between DK and D_sπ, while the diagonal DK→DK and D_sπ→D_sπ interactions vanish in the isovector channel. Solving the Lippmann-Schwinger equation yields a pole on the second Riemann sheet at about 2288 − 90i MeV. The fit to LHCb efficiency-corrected lineshapes works only when the D_sπ cutoff (Λ1≈2.18 GeV) is much larger than the DK cutoff (Λ2=0.5 GeV), an asymmetry the paper interprets as a momentum-scale or SU(3)-breaking effect. The same parameter set predicts D_s1(2536)+ → D_s+π+π− to be dominated by rescattering, giving a single broad enhancement rather than the two pe

What carries the argument

The central machinery is a coupled-channel Lippmann-Schwinger equation for the DK–D_sπ system with purely off-diagonal potentials (mediated by K* exchange), regularized by a dipole form factor with channel-dependent cutoffs. Triangle-loop diagrams feed the D_s1D*K vertex into the final D_sππ state, and the D_s1→D*K couplings are fixed from residues of the T-matrix, not fitted to the decay. The interference between the f0(500) triangle diagram and the rescattering diagram is what shapes the two-peak versus one-peak lineshapes.

Load-bearing premise

The claimed resonance depends on a hand-set channel-dependent cutoff asymmetry (Λ1≈2.18 GeV for D_sπ vs Λ2=0.5 GeV for DK); if that asymmetry is not physically justified, the T_cbar{s} pole could be an artifact of regularization rather than a dynamical prediction.

What would settle it

Measure the D_s1(2536)+ → D_s+ π+ π− invariant-mass distribution with high statistics. The framework predicts a single broad peak dominated by rescattering; observation of a clear two-peak structure, or a distinctly different peak position, would falsify the mechanism. Alternatively, a lattice QCD calculation of the I=1 DK–D_sπ scattering amplitude could check whether the off-diagonal potential alone produces a second-sheet pole without the cutoff asymmetry.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • If correct, T_cbar{s} is a genuine dynamically generated resonance on the second Riemann sheet, not merely a triangle singularity.
  • The D_s1(2536)+ → D_s+ π+ π− invariant-mass spectrum is predicted to show a single broad peak; LHCb or Belle II can test this directly.
  • The two-peak structure in D_s1(2460) decay is explained by interference between f0(500) and rescattering, distinguishing the S-wave vs D-wave coupling nature of the two D_s1 states.
  • The channel-dependent cutoff asymmetry becomes a physical feature (momentum scale or SU(3) breaking), not a technical nuisance.
  • The framework unifies spectroscopy and decay of D_s1 states, providing a template for analyzing other exotic hadrons.

Where Pith is reading between the lines

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

  • If the single-peak prediction is confirmed, the triangle-singularity explanation for T_cbar{s} would be strongly disfavored, and the off-diagonal coupled-channel mechanism would be established as the operative one.
  • The same off-diagonal mechanism may apply to other flavored tetraquark candidates where diagonal channels are weak, e.g., analogous states in bottom or charm-strange sectors.
  • The cutoff asymmetry suggests a concrete microscopic picture: the light pion in D_sπ permits larger virtual momenta, so future lattice QCD calculations could check whether the effective DK→DK interaction indeed strengthens near threshold due to D_sπ loops.
  • If the D_s1(2536) lineshape instead shows a two-peak pattern, the present framework would need revision, and the role of D-wave couplings in delaying resonance formation would be the first suspect.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 6 minor

Summary. The paper develops a unified model for the LHCb T_cbar{s} exotic state, combining triangle-loop diagrams with a coupled-channel DK–D_sπ rescattering T-matrix whose interactions are purely off-diagonal in the isovector channel. The eight-parameter model is fitted to efficiency-corrected LHCb lineshapes of D_s1(2460)+ → D_s+π+π−, giving χ²/d.o.f. = 1.43 and a second-sheet pole at E_p = 2288.4 − 89.6i MeV, identified as T_cbar{s}. Pole generation is attributed to a channel-dependent cutoff asymmetry (Λ1 = 2.18 GeV for D_sπ, Λ2 = 0.5 GeV fixed for DK). The same framework predicts a single-peak structure for D_s1(2536)+ → D_s+π+π−, proposed as the decisive experimental test.

Significance. If the central claim holds, the paper supplies a concrete dynamical mechanism for T_cbar{s} (off-diagonal DK–D_sπ generation) and a sharp, testable fingerprint distinguishing D_s1(2460) from D_s1(2536). The manuscript has real strengths: the D_s1(2536) lineshape is predicted with no additional free parameters (the overall coupling λ cancels); the D_s1D*K couplings are imported from the spectroscopy fit of Ref. [28] rather than re-fitted; the Dalitz plot is compared without direct fitting; and the fit quality is reported. However, the second-sheet pole — the linchpin of the 'dynamically generated' claim — is generated by a regularization asymmetry, is a post-diction of a fit to the same data that exhibits the structure, and is not tested against the triangle-singularity-only alternative. The robustness checks needed to support the headline claim are not presented.

major comments (4)
  1. [Sec. IV; Eq. (6); Table I] The existence of the pole is driven by the regularization asymmetry: Λ1=2.18 GeV (D_sπ) while Λ2=0.5 GeV is fixed by hand, and the text itself states that the larger D_sπ cutoff 'enhances the effective DK interaction, thereby generating a resonance near the DK threshold.' No test is shown that a common-cutoff regularization (Λ1=Λ2) or a moderate variation of Λ2 preserves the second-sheet pole and the fit quality. Since the off-diagonal potential at relevant momenta is set by the dipole form factors, the pole could be an artifact; a (Λ1,Λ2) scan reporting pole positions and χ²/d.o.f. is required to substantiate 'dynamically generated.'
  2. [Pole-extraction paragraph; Conclusions] The off-diagonal coupling g_K* and the cutoff Λ1 are fitted to the same LHCb lineshape data in which the T_cbar{s} enhancement is observed; the pole is therefore a post-diction, not an independent prediction. As the paper notes, 'By fitting the LHCb lineshapes, we determine the coupled-channel interactions ... and extract a pole.' The only genuinely independent content is the D_s1(2536) prediction. To support the claim that the pole is dynamically generated rather than a fit artifact, the authors should (i) report the fit quality with the coupled-channel pole removed (triangle diagrams only), and/or (ii) fix g_K*, Λ1, Λ2 from the spectroscopy analysis of Ref. [28] and test whether the pole survives.
  3. [Table I; D_s1(2460) fit] The identified pole, E_p = 2288.4 − 89.6i MeV, implies a pole mass about 39 MeV below and a width (~179 MeV) nearly twice the LHCb-reported values (mass 2327±13 MeV, width 96±16, +170/−23 MeV). Since the pole is presented as 'corresponding to T_cbar{s},' the mapping between the second-sheet pole and the experimental line-shape parameters needs explicit discussion (threshold/Flatté-type distortion, line-shape shift); otherwise the identification is not quantitatively established.
  4. [Sec. II; Eq. (2); Fig. 2] Refs. [6,7] interpret the same LHCb structure as a triangle singularity without a genuine resonance. Although the present model includes both triangle loops and rescattering, no quantitative comparison with a triangle-only fit (diagrams (a)+(b) without the T-matrix pole) is presented. The central message that the two-peak structure requires the coupled-channel pole would be strengthened by reporting Δχ²/d.o.f. for the pole-less alternative; as written, the fitted data cannot disambiguate the two mechanisms.
minor comments (6)
  1. [Abstract] 'Reimann Sheet' should be 'Riemann Sheet'; 'combing' should be 'combining.'
  2. [Sec. III; Table I] The text says the model has eight free parameters, but Table I lists seven fitted parameters with Λ2 fixed. Please reconcile the counting.
  3. [Sec. II and Fig. 2 caption] The caption calls the panels 'the fitted lineshapes of T_cbar{s}'; the panels are actually the D_s+π+ and π+π− spectra of D_s1(2460)+ decay. Please clarify the caption and state that the middle panel applies the m(π+π−)>0.39 GeV selection.
  4. [Sec. III] Please clarify whether χ²/d.o.f. = 1.43 is computed against the pseudo-data generated from the LHCb model or against the measured efficiency-corrected data. Fitting pseudo-data generated from LHCb's amplitude model may partially mask genuine discrepancies between the model and the actual measurements.
  5. [Sec. V] The statement that the theoretical and experimental g_S/g_D ratios give 'very similar invariant-mass distributions' despite a phase difference of about π is asserted without a comparison plot. Please quantify the difference or show the corresponding lineshapes.
  6. [Sec. IV] The asserted cancellation of the ρ/ω-exchange diagonal potentials in the I=1 channel is load-bearing for the 'rather than by diagonal interactions' conclusion but is presented in one sentence. A one-line isospin derivation or an explicit reference should be added.

Circularity Check

0 steps flagged

No significant circularity: the T_cbar{s} pole is a fitted output and the D_s1(2536) lineshape is an independent parameter-free prediction, though cutoff sensitivity is a caveat.

full rationale

The derivation chain is not circular in the sense defined here. The T_cbar{s} pole is not an input fitted as a parameter; it is obtained by solving the coupled-channel Lippmann-Schwinger equation (Eq. 6) with parameters (Λ1, gK*, r1, r2, mf0, Γf0, phase) that are fitted to the LHCb efficiency-corrected lineshapes, and then locating the pole by complex scaling. That is a fitted extraction, not a prediction from a quantity defined in terms of the pole. The assumption that diagonal DK and Dsπ potentials vanish in the I=1 system is a symmetry-based input; the summary statement that the pole is generated by the off-diagonal potential is an implication of that input, not a tautological recycle of the conclusion. The genuinely predictive element is the D_s1(2536)+ → D_s+ π+π− lineshape, computed with the same fitted T-matrix and the D_s1D*K couplings from Ref. [28], with no new free parameters (λ cancels, gD/gS can be replaced by the experimental value). Ref. [28] is a self-citation, but it is anchored to lattice-QCD spectra of the D(*)K system and is therefore external support rather than a self-justifying uniqueness argument. The main caveat—channel-dependent cutoffs Λ1=2.18 GeV (fitted) vs Λ2=0.5 GeV (fixed)—is a robustness concern; the paper itself admits 'the location and nature of the resonance pole are sensitive to the details of the coupled-channel dynamics.' However, model-dependence of a fitted pole is not circularity, because the pole-producing effect of the larger Dsπ cutoff is an output of the fit and is separately testable through the D_s1(2536) lineshape prediction.

Axiom & Free-Parameter Ledger

8 free parameters · 5 axioms · 0 invented entities

The paper's central result is model-dependent: eight parameters are fitted to the LHCb lineshapes, and the pole and its Riemann-sheet assignment are outputs of that fit, not consequences of first principles. The dominant free-parameter burden is the cutoff pair (Λ1, Λ2), whose asymmetry produces the resonance.

free parameters (8)
  • Λ1 (D_sπ cutoff) = 2.18^{+0.24}_{-0.04} GeV
    Dipole form-factor cutoff for the D_sπ channel; drives the effective DK interaction and is load-bearing for the T_cbar{s} pole (Table I).
  • Λ2 (DK cutoff) = 0.5 GeV (fixed)
    Second dipole cutoff, fixed by hand; the asymmetry with Λ1 is essential for generating the pole.
  • g_K* = 55.3^{+0.8}_{-2.5}
    Overall off-diagonal coupling constant in the effective potential; controls the DK–D_sπ transition strength.
  • ϕ = 3.78^{+0.38}_{-0.26} rad
    Relative phase between diagram (a) and diagrams (b,c); essential for the two-peak interference pattern.
  • r1 = 215^{+51}_{-86}
    Overall coupling of triangle diagram (a).
  • r2 = -9.0^{+4.3}_{-0.8}
    Overall coupling of rescattering diagrams (b,c).
  • m_f0 = 519^{+31}_{-89} MeV
    f0(500) Breit-Wigner mass, fitted to the data.
  • Γ_f0 = 242^{+90}_{-88} MeV
    f0(500) Breit-Wigner width, fitted to the data.
axioms (5)
  • domain assumption In I=1, diagonal DK→DK and D_sπ→D_sπ potentials vanish due to cancellation of ρ/ω exchange; only the off-diagonal potential remains.
    Stated after Eq. (6): 'in the isovector configuration, these contributions cancel, leaving vanishing diagonal interactions.' This zeroing is central; any SU(3) breaking would alter the coupled-channel dynamics.
  • domain assumption D_s1(2460) and D_s1(2536) couple predominantly to D*K; direct c̄s decays are OZI-suppressed.
    Section II: 'Since both D_s1(2460) and D_s1(2536) predominantly couple to D*K, the triangle diagram mechanism naturally dominates.' Underpins the entire diagram set.
  • domain assumption Only DK and D_sπ coupled channels are relevant for T_cbar{s} dynamics; other thresholds are neglected.
    The T-matrix in Eq. (6) is defined with α,β,γ=1,2 for D_sπ and DK only.
  • domain assumption The f0(500) is represented by a Breit-Wigner propagator; f0(980) is neglected because it lies outside the accessible ππ range.
    Footnote 1 and Eq. (2) diagram (a): the f0(980) contribution 'can safely be neglected'.
  • domain assumption The D_s1D*K S- and D-wave couplings are taken from the T-matrix residues of Ref. [28].
    Section II: 'the D_s1D*K couplings used here are fixed from the pole residues of the coupled-channel T-matrix analysis in Ref. [28].' This is a self-cited phenomenological input, not independently derived here.

pith-pipeline@v1.3.0-alltime-deepseek · 9722 in / 11593 out tokens · 110265 ms · 2026-08-04T12:51:02.439084+00:00 · methodology

0 comments
read the original abstract

The $T_{c\bar{s}}$ state observed in the decay $D_{s1}(2460)^+ \to D_s^+\pi^+\pi^-$ provides direct evidence for an isovector open-charm tetraquark state with strangeness--a discovery that demands a systematic framework connecting its origin to the nature of the parent $D_{s1}$. We successfully achieve this connection by two mechanisms, triangle loops and the coupled channel of $DK$-$D_s\pi$ with pure off-diagonal potentials. We first point out the behavior of propagator of $D_s\pi$ will influence the effective potential of $DK\to DK$, then we can successfully obtain the pole of $T_{c\bar{s}}$ on the second Reimann Sheet. By combing with the $\pi\pi$-$KK$ rescattering, not only the two-peak structure in $D_{s1}(2460)$ decay is well reproduced, but also a single-peak structure is predicted in $D_{s1}(2536)$ decay. The marked difference, testable at LHCb and Belle II, is driven by the $S$-wave versus $D$-wave nature of their $D^*K$ couplings, revealing the underlying structural distinction between the two $D_{s1}$ states. By directly linking hadronic structure to decay patterns, this work provides a template for deciphering the nature of such exotic states. More broadly, by revealing how non-perturbative coupled-channel effects manifest in exotic hadrons, our analysis connects to a universal mechanism shared by systems ranging from halo nuclei to atomic Feshbach resonances, offering a unified perspective across these fields.

Figures

Figures reproduced from arXiv: 2510.01564 by Guang-Juan Wang, Jia-Jun Wu, Makoto Oka, Zhi Yang.

Figure 1
Figure 1. Figure 1: FIG. 1. The illustrative Feynman diagrams for [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The fitted lineshapes of the [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. The Dalitz plot of the [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. The predicted lineshape of the [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. The predicted dalitz plot of the [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗
Figure 2
Figure 2. Figure 2: Meanwhile, the π +π − spectrum shows the same features, because the contribution from [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗

discussion (0)

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

Reference graph

Works this paper leans on

39 extracted references · 33 linked inside Pith

  1. [1]

    Aaij et al

    R. Aaij et al. (LHCb), (2024), arXiv:2411.03399 [hep-ex]

  2. [2]

    E. E. Kolomeitsev and M. F. M. Lutz, Phys. Lett. B582, 39 (2004), arXiv:hep-ph/0307133

  3. [3]

    F.-K. Guo, C. Hanhart, and U.-G. Meissner, Eur. Phys. J. A40, 171 (2009), arXiv:0901.1597 [hep-ph]

  4. [4]

    Guo, U.-G

    Z.-H. Guo, U.-G. Meißner, and D.-L. Yao, Phys. Rev. D 92, 094008 (2015), arXiv:1507.03123 [hep-ph]

  5. [5]

    Wang, Y.-S

    Z.-Y. Wang, Y.-S. Li, and S.-Q. Luo, (2024), arXiv:2412.06446 [hep-ph]

  6. [6]

    L. Roca, J. M. Dias, and E. Oset, (2025), arXiv:2502.18401 [hep-ph]

  7. [7]

    J. M. Dias, Y.-Y. Li, and E. Oset, (2025), arXiv:2507.20425 [hep-ph]

  8. [8]

    Godfrey and N

    S. Godfrey and N. Isgur, Phys. Rev. D32, 189 (1985)

  9. [9]

    A. P. Szczepaniak, Phys. Lett. B567, 23 (2003), arXiv:hep-ph/0305060

  10. [10]

    Hofmann and M

    J. Hofmann and M. F. M. Lutz, Nucl. Phys. A733, 142 (2004), arXiv:hep-ph/0308263

  11. [11]

    van Beveren and G

    E. van Beveren and G. Rupp, Phys. Rev. Lett.91, 012003 (2003), arXiv:hep-ph/0305035

  12. [12]

    Barnes, F

    T. Barnes, F. E. Close, and H. J. Lipkin, Phys. Rev. D 68, 054006 (2003), arXiv:hep-ph/0305025

  13. [13]

    Guo, P.-N

    F.-K. Guo, P.-N. Shen, H.-C. Chiang, R.-G. Ping, and B.-S. Zou, Phys. Lett. B641, 278 (2006), arXiv:hep- ph/0603072

  14. [14]

    Albaladejo, P

    M. Albaladejo, P. Fernandez-Soler, J. Nieves, and P. G. Ortega, Eur. Phys. J. C78, 722 (2018), arXiv:1805.07104 [hep-ph]

  15. [15]

    Wu, M.-Z

    T.-W. Wu, M.-Z. Liu, L.-S. Geng, E. Hiyama, and M. P. Valderrama, Phys. Rev. D100, 034029 (2019), arXiv:1906.11995 [hep-ph]

  16. [16]

    Kong, J.-T

    S.-Y. Kong, J.-T. Zhu, D. Song, and J. He, (2021), arXiv:2106.07272 [hep-ph]

  17. [17]

    E. B. Gregory, F.-K. Guo, C. Hanhart, S. Krieg, and T. Luu, (2021), arXiv:2106.15391 [hep-ph]

  18. [18]

    Wang and X

    P. Wang and X. G. Wang, Phys. Rev. D86, 014030 (2012), arXiv:1204.5553 [hep-ph]

  19. [19]

    Huang, Z.-Y

    B.-L. Huang, Z.-Y. Lin, and S.-L. Zhu, Phys. Rev. D 105, 036016 (2022), arXiv:2112.13702 [hep-ph]

  20. [20]

    Maiani, F

    L. Maiani, F. Piccinini, A. D. Polosa, and V. Riquer, Phys. Rev. D71, 014028 (2005), arXiv:hep-ph/0412098

  21. [21]

    Dai, X.-Q

    Y.-B. Dai, X.-Q. Li, S.-L. Zhu, and Y.-B. Zuo, Eur. Phys. J. C55, 249 (2008), arXiv:hep-ph/0610327

  22. [22]

    H.-X. Chen, W. Chen, X. Liu, Y.-R. Liu, and S.-L. Zhu, Rept. Prog. Phys.80, 076201 (2017), arXiv:1609.08928 [hep-ph]

  23. [23]

    F.-K. Guo, C. Hanhart, U.-G. Meißner, Q. Wang, Q. Zhao, and B.-S. Zou, Rev. Mod. Phys.90, 015004 (2018), [Erratum: Rev.Mod.Phys. 94, 029901 (2022)], arXiv:1705.00141 [hep-ph]

  24. [24]

    H.-X. Chen, W. Chen, X. Liu, and S.-L. Zhu, Phys. Rept.639, 1 (2016), arXiv:1601.02092 [hep-ph]

  25. [25]

    Esposito, A

    A. Esposito, A. Pilloni, and A. D. Polosa, Phys. Rept. 668, 1 (2017), arXiv:1611.07920 [hep-ph]

  26. [26]

    Brambilla, S

    N. Brambilla, S. Eidelman, C. Hanhart, A. Nefediev, C.- P. Shen, C. E. Thomas, A. Vairo, and C.-Z. Yuan, Phys. Rept.873, 1 (2020), arXiv:1907.07583 [hep-ex]. 7

  27. [27]

    Y. S. Kalashnikova and A. V. Nefediev, Phys. Usp.62, 568 (2019), arXiv:1811.01324 [hep-ph]

  28. [28]

    Yang, G.-J

    Z. Yang, G.-J. Wang, J.-J. Wu, M. Oka, and S.-L. Zhu, Phys. Rev. Lett.128, 112001 (2022), arXiv:2107.04860 [hep-ph]

  29. [29]

    C. B. Lang, L. Leskovec, D. Mohler, S. Prelovsek, and R. M. Woloshyn, Phys. Rev. D90, 034510 (2014), arXiv:1403.8103 [hep-lat]

  30. [30]

    G. S. Bali, S. Collins, A. Cox, and A. Sch¨ afer, Phys. Rev. D96, 074501 (2017), arXiv:1706.01247 [hep-lat]

  31. [31]

    Navas et al

    S. Navas et al. (Particle Data Group), Phys. Rev. D110, 030001 (2024)

  32. [32]

    Matsuyama, T

    A. Matsuyama, T. Sato, and T. S. H. Lee, Phys. Rept. 439, 193 (2007), arXiv:nucl-th/0608051

  33. [33]

    J.-J. Wu, T. S. H. Lee, and B. S. Zou, Phys. Rev. C85, 044002 (2012), arXiv:1202.1036 [nucl-th]

  34. [34]

    Shtabovenko, R

    V. Shtabovenko, R. Mertig, and F. Orellana, Comput. Phys. Commun.306, 109357 (2025), arXiv:2312.14089 [hep-ph]

  35. [35]

    Hahn and M

    T. Hahn and M. Perez-Victoria, Comput. Phys. Com- mun.118, 153 (1999), arXiv:hep-ph/9807565

  36. [36]

    Kamano, S

    H. Kamano, S. X. Nakamura, T. S. H. Lee, and T. Sato, Phys. Rev. D84, 114019 (2011), arXiv:1106.4523 [hep- ph]

  37. [37]

    Ikeno, R

    N. Ikeno, R. Molina, and E. Oset, Phys. Rev. D 105, 014012 (2022), [Erratum: Phys.Rev.D 106, 099905 (2022)], arXiv:2111.05024 [hep-ph]

  38. [38]

    Wang and X

    F.-L. Wang and X. Liu, Phys. Lett. B835, 137583 (2022), arXiv:2207.10493 [hep-ph]

  39. [39]

    Aaij et al

    R. Aaij et al. (LHCb), JHEP10, 106 (2023), arXiv:2308.00587 [hep-ex]