pith. sign in

arxiv: 2606.13278 · v1 · pith:UNIGAKF3new · submitted 2026-06-11 · ✦ hep-ph

Resolving the CP Asymmetry Puzzle in B Decays with Unitarized Final-State Interactions

Pith reviewed 2026-06-27 06:27 UTC · model grok-4.3

classification ✦ hep-ph
keywords B meson decaysCP asymmetryfinal-state interactionsunitarizationchiral effective field theoryheavy-quark spin symmetryannihilation decaysD Dbar system
0
0 comments X

The pith

A parameter-free unitarized FSI framework using gamma-gamma data predicts CP asymmetries in B decays including large effects in pure-annihilation channels.

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

The paper develops a unitarized description of final-state interactions in hadronic B decays by solving the Lippmann-Schwinger equation to all orders. The interaction kernel comes from chiral effective field theory constrained by heavy-quark spin symmetry, with all low-energy constants fixed exclusively by BaBar measurements of the gamma-gamma to D Dbar cross section. This leaves no free parameters when the framework is applied to B decays. The resulting predictions match existing data on CP asymmetries and branching fractions. Pure-annihilation modes that give zero asymmetry in short-distance treatments acquire sizable asymmetries from the long-distance dynamics.

Core claim

Using the coupled D Dbar system, the interaction kernel can be derived from chiral effective field theory constrained by heavy-quark spin symmetry, and low-energy constants are fixed by the gamma gamma to D Dbar cross-section data from BaBar, leaving no adjustable parameter in the FSI sector when applied to B decays. The resulting predictions for CP asymmetries and branching fractions show excellent agreement with the available experimental data. A defining consequence is that CP asymmetries in the pure-annihilation channels B0bar to D0 D0bar and B0bar to Ds+ Ds-, identically zero in any short-distance treatment, are dramatically enhanced by FSI.

What carries the argument

The unitarized final-state interaction framework based on the Lippmann-Schwinger equation solved to all orders in the coupled D Dbar system with a chiral EFT kernel fixed by gamma-gamma data.

If this is right

  • CP asymmetries in the pure-annihilation channels B0bar to D0 D0bar and B0bar to Ds+ Ds- are dramatically enhanced by FSI.
  • Predictions for CP asymmetries and branching fractions agree with available experimental data.
  • Definite predictions are given for partial widths and CP asymmetries of all yet-unmeasured channels.
  • A measured nonzero asymmetry in pure-annihilation modes serves as a direct signature of long-distance dynamics.
  • The framework carries implications for CKM parameter extraction and searches for physics beyond the Standard Model.

Where Pith is reading between the lines

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

  • The same unitarized approach may resolve analogous strong-phase puzzles in other nonleptonic heavy-meson decays.
  • Future precision measurements in the pure-annihilation channels can distinguish the size of long-distance contributions from short-distance expectations.
  • Incorporating these FSI effects could tighten constraints on CKM angles extracted from multiple B decay modes simultaneously.

Load-bearing premise

The low-energy constants determined from gamma-gamma to D Dbar data can be transferred without modification or additional intermediate states to describe final-state interactions in B decays.

What would settle it

A measurement of the CP asymmetry in B0bar to D0 D0bar that is consistent with zero within experimental precision would contradict the predicted dramatic enhancement from final-state interactions.

Figures

Figures reproduced from arXiv: 2606.13278 by Pengyu Niu, Qian Wang, Wei Wang, Xin-Yue Hu.

Figure 1
Figure 1. Figure 1: FIG. 1: The schematic diagram of [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Comparison of the predicted partial widths [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: The [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: The topologies considered in this work. [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: From left to right, the [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: The distribution of the parameters approximates a normal distribution. The average of the distribution [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: Distributions of weak decays parameters. The upper, middle and lower panels correspond to scheme I, [PITH_FULL_IMAGE:figures/full_fig_p014_8.png] view at source ↗
read the original abstract

A reliable description of direct CP asymmetries in hadronic $B$ decays, among the most sensitive probes of the CKM mechanism and of physics beyond the Standard Model, remains unresolved. Conventional factorization approaches adopt distinct treatments for the strong phases originating from long-distance QCD interactions, leading to predictions that differ by factors of several for identical decay channels. Existing final-state interaction (FSI) models are limited to one-loop approximations that break unitarity and rely on channel-specific phenomenological cutoffs, severely restricting the predictive capability. We introduce a unitarized FSI framework based on the Lippmann-Schwinger equation solved to all orders, restoring unitarity that is manifestly broken in one-loop treatments. Using the coupled $D\bar{D}$ system, we show that the interaction kernel can be derived from chiral effective field theory constrained by heavy-quark spin symmetry, and low-energy constants are fixed by the $\gamma\gamma\to D\bar{D}$ cross-section data from BaBar, leaving no adjustable parameter in the FSI sector when applied to $B$ decays. The resulting predictions for CP asymmetries and branching fractions show excellent agreement with the available experimental data. A defining consequence is that CP asymmetries in the pure-annihilation channels $\bar{B}^0\to D^0\bar{D}^0$ and $\bar{B}^0\to D_s^+D_s^-$, identically zero in any short-distance treatment, are dramatically enhanced by FSI, and a measured nonzero asymmetry in these modes is therefore a direct experimental signature of long-distance dynamics. We present definite predictions for the partial widths and CP asymmetries of all yet-unmeasured channels, establishing final-state interactions as a predictive, data-driven ingredient of the Standard Model with direct implications for CKM parameter extraction and BSM searches.

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

2 major / 1 minor

Summary. The manuscript introduces a unitarized final-state interaction (FSI) framework for hadronic B decays based on the Lippmann-Schwinger equation solved to all orders. The interaction kernel for the coupled D Dbar system is derived from chiral effective field theory constrained by heavy-quark spin symmetry; its low-energy constants are fixed exclusively by BaBar γγ → D Dbar cross-section data near threshold. No additional parameters are introduced when the resulting T-matrix is applied to B → D Dbar amplitudes at the B-meson mass shell. The approach is claimed to restore unitarity (broken in one-loop treatments), yield predictions for CP asymmetries and branching fractions in agreement with existing data, and generate large CP asymmetries in pure-annihilation channels (B0bar → D0 D0bar and B0bar → Ds+ Ds-) that vanish in short-distance treatments.

Significance. If the central results hold, the work would provide a parameter-free, data-driven method for incorporating long-distance strong phases in B decays, addressing long-standing discrepancies among factorization approaches. The explicit predictions for yet-unmeasured channels and the identification of nonzero CP asymmetries in annihilation modes as direct signatures of FSI would have implications for CKM extractions and BSM searches. The use of external γγ data to fix the kernel and the all-orders unitarization are notable strengths.

major comments (2)
  1. [framework application to B decays] The transfer of the interaction kernel (fixed by BaBar γγ → D Dbar data at √s ≲ 4.5 GeV) to B-decay kinematics at s = m_B² ≈ 28 GeV² is load-bearing for the 'no adjustable parameter' claim and the predicted CP asymmetries. The manuscript does not demonstrate that higher-order momentum dependence or additional coupled channels remain negligible over this factor-of-six increase in energy scale; uncontrolled corrections could alter the strong phases extracted from the Lippmann-Schwinger solution.
  2. [results and comparison with data] The abstract states that the resulting predictions show 'excellent agreement' with data, yet the manuscript provides no quantitative assessment (e.g., χ² per degree of freedom or explicit comparison tables) of how sensitive this agreement is to variations in the kernel at high s. Without such a test, the agreement cannot be distinguished from possible fine-tuning within the extrapolation uncertainty.
minor comments (1)
  1. [formalism] Notation for the coupled-channel T-matrix and the precise definition of the on-shell projection in the Lippmann-Schwinger equation should be clarified with an explicit equation reference.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough review and valuable comments on our manuscript. We address each major comment below in detail. Our responses aim to clarify the robustness of the framework while acknowledging areas where additional clarification or revision can strengthen the presentation.

read point-by-point responses
  1. Referee: [framework application to B decays] The transfer of the interaction kernel (fixed by BaBar γγ → D Dbar data at √s ≲ 4.5 GeV) to B-decay kinematics at s = m_B² ≈ 28 GeV² is load-bearing for the 'no adjustable parameter' claim and the predicted CP asymmetries. The manuscript does not demonstrate that higher-order momentum dependence or additional coupled channels remain negligible over this factor-of-six increase in energy scale; uncontrolled corrections could alter the strong phases extracted from the Lippmann-Schwinger solution.

    Authors: We appreciate the referee pointing out this critical aspect of the extrapolation. The interaction kernel is constructed from chiral effective field theory with constraints from heavy-quark spin symmetry, which provides a systematic expansion valid in the relevant kinematic regime. The low-energy constants are determined solely from the BaBar data near threshold, and the all-orders solution of the Lippmann-Schwinger equation generates the T-matrix at higher energies without introducing new parameters. While we recognize that an explicit demonstration of the suppression of higher-order terms over this energy range is not provided in the current manuscript, the framework's success in reproducing data in the D Dbar system and the resulting agreement in B decays support its applicability. We will add a paragraph discussing the expected magnitude of corrections based on the power counting of the EFT in the revision. revision: partial

  2. Referee: [results and comparison with data] The abstract states that the resulting predictions show 'excellent agreement' with data, yet the manuscript provides no quantitative assessment (e.g., χ² per degree of freedom or explicit comparison tables) of how sensitive this agreement is to variations in the kernel at high s. Without such a test, the agreement cannot be distinguished from possible fine-tuning within the extrapolation uncertainty.

    Authors: We agree that including a quantitative assessment would improve the clarity of the results section. Although the model has no adjustable parameters in the FSI sector, making fine-tuning unlikely, we will revise the manuscript to include explicit comparison tables of predicted versus measured branching fractions and CP asymmetries, along with a discussion of the sensitivity to reasonable variations in the kernel. This will include estimates of theoretical uncertainties arising from the extrapolation. revision: yes

Circularity Check

0 steps flagged

No significant circularity; LECs fixed by external BaBar data

full rationale

The paper states that the interaction kernel is derived from chiral EFT + HQSS and its low-energy constants are fixed solely by independent BaBar γγ→D Dbar cross-section data, with the resulting unitarized amplitudes then applied to B decays without further parameters. No self-citation load-bearing steps, self-definitional relations, or fitted inputs renamed as predictions appear in the abstract or described framework. The central claim rests on external experimental input rather than internal reduction, making the derivation self-contained against the stated benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 3 axioms · 0 invented entities

No free parameters are introduced for the B-decay application; constants are fixed by external BaBar data. The framework rests on standard domain assumptions of chiral EFT and heavy-quark symmetry plus the mathematical property that the Lippmann-Schwinger equation restores unitarity.

axioms (3)
  • domain assumption Heavy-quark spin symmetry constrains the interaction kernel in the coupled DDbar system
    Invoked to derive the kernel from chiral effective field theory.
  • standard math Solution of the Lippmann-Schwinger equation to all orders restores unitarity
    Mathematical basis for replacing one-loop approximations.
  • domain assumption The coupled DDbar system suffices to model the dominant FSI relevant to the B decays under study
    The paper restricts attention to this channel.

pith-pipeline@v0.9.1-grok · 5862 in / 1789 out tokens · 33667 ms · 2026-06-27T06:27:42.365706+00:00 · methodology

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

54 extracted references · 30 linked inside Pith

  1. [1]

    A. D. Sakharov, Violation of CP Invariance, C asymme- try, and baryon asymmetry of the universe, Pisma Zh. Eksp. Teor. Fiz.5, 32 (1967)

  2. [2]

    Bona et al

    M. Bona et al. (UTfit), New UTfit Analysis of the Unitarity Triangle in the Cabibbo-Kobayashi-Maskawa scheme, Rend. Lincei Sci. Fis. Nat.34, 37 (2023), arX- iv:2212.03894 [hep-ph]

  3. [3]

    Beneke, G

    M. Beneke, G. Buchalla, M. Neubert, and C. T. Sachrajda, QCD factorization for B —>pi pi decays: Strong phases and CP violation in the heavy quark limit, Phys. Rev. Lett.83, 1914 (1999), arXiv:hep-ph/9905312

  4. [4]

    Beneke, G

    M. Beneke, G. Buchalla, M. Neubert, and C. T. Sachrajda, QCD factorization for exclusive, nonlepton- ic B meson decays: General arguments and the case of heavy light final states, Nucl. Phys. B591, 313 (2000), arXiv:hep-ph/0006124

  5. [5]

    Y. Y. Keum, H.-N. Li, and A. I. Sanda, Penguin enhance- ment andB→Kπdecays in perturbative QCD, Phys. Rev. D63, 054008 (2001), arXiv:hep-ph/0004173

  6. [6]

    C.-D. Lu, K. Ukai, and M.-Z. Yang, Branching ratio and CP violation of B —>pi pi decays in perturbative QCD approach, Phys. Rev. D63, 074009 (2001), arXiv:hep- ph/0004213

  7. [7]

    C. W. Bauer, S. Fleming, D. Pirjol, and I. W. Stewart, An Effective field theory for collinear and soft gluons: Heavy to light decays, Phys. Rev. D63, 114020 (2001), arXiv:hep-ph/0011336

  8. [8]

    C. W. Bauer, D. Pirjol, I. Z. Rothstein, and I. W. Stewart,B→M 1M2: Factorization, charming penguins, strong phases, and polarization, Phys. Rev. D70, 054015 (2004), arXiv:hep-ph/0401188

  9. [9]

    Beneke and M

    M. Beneke and M. Neubert, QCD factorization forB→ P PandB→P Vdecays, Nucl. Phys. B675, 333 (2003), arXiv:hep-ph/0308039

  10. [10]

    Cheng and C.-K

    H.-Y. Cheng and C.-K. Chua, QCD Factorization for Charmless HadronicB s Decays Revisited, Phys. Rev. D 80, 114026 (2009), arXiv:0910.5237 [hep-ph]

  11. [11]

    Cheng and C.-K

    H.-Y. Cheng and C.-K. Chua, Revisiting Charmless HadronicB u,d Decays in QCD Factorization, Phys. Rev. D80, 114008 (2009), arXiv:0909.5229 [hep-ph]

  12. [12]

    J. Chai, S. Cheng, Y.-h. Ju, D.-C. Yan, C.-D. L¨ u, and Z.-J. Xiao, Charmless two-body B meson decays in the perturbative QCD factorization approach, Chin. Phys. C 46, 123103 (2022), arXiv:2207.04190 [hep-ph]

  13. [13]

    Z.-T. Zou, A. Ali, C.-D. Lu, X. Liu, and Y. Li, Improved Estimates of TheB (s) →V VDecays in Perturbative QCD Approach, Phys. Rev. D91, 054033 (2015), arX- iv:1501.00784 [hep-ph]

  14. [14]

    A. R. Williamson and J. Zupan, Two body B decays with isosinglet final states in SCET, Phys. Rev. D74, 014003 (2006), [Erratum: Phys.Rev.D 74, 03901 (2006)], arXiv:hep-ph/0601214

  15. [15]

    Wang, Y.-M

    W. Wang, Y.-M. Wang, D.-S. Yang, and C.-D. Lu, Charmless Two-bodyB (s) →V Pdecays In Soft- Collinear-Effective-Theory, Phys. Rev. D78, 034011 (2008), arXiv:0801.3123 [hep-ph]

  16. [16]

    Cheng, C.-K

    H.-Y. Cheng, C.-K. Chua, and A. Soni, Final state inter- actions in hadronic B decays, Phys. Rev. D71, 014030 (2005), arXiv:hep-ph/0409317

  17. [17]

    Bediaga and C

    I. Bediaga and C. G¨ obel, DirectCPviolation in beauty and charm hadron decays, Prog. Part. Nucl. Phys.114, 103808 (2020), arXiv:2009.07037 [hep-ex]

  18. [18]

    A. Pich, E. Solomonidi, and L. Vale Silva, Final-state interactions in the CP asymmetries of charm-meson two- body decays, Phys. Rev. D108, 036026 (2023), arX- iv:2305.11951 [hep-ph]

  19. [19]

    Bediaga, T

    I. Bediaga, T. Frederico, and P. C. Magalh˜ aes, Enhanced Charm CP Asymmetries from Final State Interactions, Phys. Rev. Lett.131, 051802 (2023), arXiv:2203.04056 [hep-ph]

  20. [20]

    A. F. Falk, A. L. Kagan, Y. Nir, and A. A. Petrov, Final state interactions and new physics inB→piKdecays, Phys. Rev. D57, 4290 (1998), arXiv:hep-ph/9712225

  21. [21]

    J. F. Donoghue, E. Golowich, A. A. Petrov, and J. M. Soares, Systematics of soft final state interactions in B decay, Phys. Rev. Lett.77, 2178 (1996), arXiv:hep- ph/9604283

  22. [22]

    Mohammadi and H

    B. Mohammadi and H. Mehraban, Final state interaction inB 0 →D 0 ¯D0, JHEP07, 089

  23. [23]

    Aubert et al

    B. Aubert et al. (BaBar), Observation of theχ c2(2p) Meson in the Reactionγγ→D ¯Dat BaBar, Phys. Rev. D81, 092003 (2010), arXiv:1002.0281 [hep-ex]

  24. [24]

    A. J. Buras and L. Silvestrini, Nonleptonic two-body B decays beyond factorization, Nucl. Phys. B569, 3 (2000), arXiv:hep-ph/9812392

  25. [25]

    Ciuchini, E

    M. Ciuchini, E. Franco, G. Martinelli, and L. Silvestrini, Charming penguins in B decays, Nucl. Phys. B501, 271 (1997), arXiv:hep-ph/9703353

  26. [26]

    Geng, X.-N

    C.-Q. Geng, X.-N. Jin, and C.-W. Liu, Large CP asym- metries from final-state interactions in charmful baryonic decays of B0→Ξc+Ξ¯c- and Bs0→Λc+Λ¯c-, Phys. Rev. D112, 013001 (2025), arXiv:2502.12770 [hep-ph]

  27. [27]

    Duan, J.-P

    Z.-D. Duan, J.-P. Wang, R.-H. Li, C.-D. L¨ u, and F.-S. Yu, Final-state rescattering mechanism in bottom-baryon de- cays, JHEP09, 160, arXiv:2412.20458 [hep-ph]

  28. [28]

    Bediaga, T

    I. Bediaga, T. Frederico, and P. C. Magalh˜ aes, Charm rescattering contribution in rareB + c →K +K −π+ decay, Phys. Lett. B785, 581 (2018), arXiv:1808.02945 [hep- ph]

  29. [29]

    Zhou, Y.-B

    S.-H. Zhou, Y.-B. Wei, Q. Qin, Y. Li, F.-S. Yu, and C.-D. Lu, Analysis of Two-body CharmedBMeson Decays in Factorization-Assisted Topological-Amplitude Approach, Phys. Rev. D92, 094016 (2015), arXiv:1509.04060 [hep- ph]

  30. [30]

    Zhou, R.-H

    S.-H. Zhou, R.-H. Li, and X.-Y. L¨ u, Analysis of three- body decays B→D(V→)PP under the factorization- assisted topological-amplitude approach, Phys. Rev. D 110, 056001 (2024), arXiv:2406.00373 [hep-ph]

  31. [31]

    Ou-Yang, R.-H

    J. Ou-Yang, R.-H. Li, and S.-H. Zhou, Analysis of three- body charmed B meson decays B→D(V*→)VP, Phys. Rev. D112, 056005 (2025), arXiv:2506.14675 [hep-ph]

  32. [32]

    Fu-Sheng, X.-X

    Y. Fu-Sheng, X.-X. Wang, and C.-D. Lu, Nonleptonic Two Body Decays of Charmed Mesons, Phys. Rev. D 84, 074019 (2011), arXiv:1101.4714 [hep-ph]

  33. [33]

    Colangelo, F

    P. Colangelo, F. De Fazio, and T. N. Pham, Nonfactorizable contributions in B decays to charmoni- um: The Case ofB − →K −hc, Phys. Rev. D69, 054023 (2004), arXiv:hep-ph/0310084

  34. [34]

    Colangelo, F

    P. Colangelo, F. De Fazio, and T. N. Pham,B − → K −χc0 decay from charmed meson rescattering, Phys. Lett. B542, 71 (2002), arXiv:hep-ph/0207061

  35. [35]

    Rohrken et al

    M. Rohrken et al. (Belle), Measurements of Branching 6 Fractions and Time-dependent CP Violating Asymmetries inB 0 →D (∗)±D∓ Decays, Phys. Rev. D 85, 091106 (2012), arXiv:1203.6647 [hep-ex]

  36. [36]

    Zupanc et al., Improved measurement of ¯B0 →D − s D+ and search for ¯B0 →D + s D− s at Belle, Phys

    A. Zupanc et al., Improved measurement of ¯B0 →D − s D+ and search for ¯B0 →D + s D− s at Belle, Phys. Rev. D75, 091102 (2007), arXiv:hep-ex/0703040

  37. [37]

    Aaij et al

    R. Aaij et al. (LHCb), First observations of ¯B0 s → D+D−,D + s D− andD 0 ¯D0 decays, Phys. Rev. D87, 092007 (2013), arXiv:1302.5854 [hep-ex]

  38. [38]

    Abudin´ enet al

    F. Abudin´ enet al. (Belle-II), Measurement of theB 0 lifetime and flavor-oscillation frequency using hadronic decays reconstructed in 2019–2021 Belle II data, Phys. Rev. D107, L091102 (2023), arXiv:2302.12791 [hep-ex]

  39. [39]

    Aaij et al

    R. Aaij et al. (LHCb), Measurement ofB 0 s andD − s me- son lifetimes, Phys. Rev. Lett.119, 101801 (2017), arX- iv:1705.03475 [hep-ex]

  40. [40]

    Aaij et al

    R. Aaij et al. (LHCb), Measurements of theB +, B0, B0 s meson and Λ 0 b baryon lifetimes, JHEP04, 114, arX- iv:1402.2554 [hep-ex]

  41. [41]

    Adachi et al

    I. Adachi et al. (Belle-II), Measurement of the branching fractions of B→D (∗)K −K (∗)0 (S) and B→D (∗)D− s decays at Belle II, JHEP08, 206, arXiv:2406.06277 [hep-ex]

  42. [42]

    Adachi et al

    I. Adachi et al. (Belle), Measurement of the branch- ing fraction and charge asymmetry of the decayB + → D+ ¯D0 and search forB 0 →D 0 ¯D0, Phys. Rev. D77, 091101 (2008), arXiv:0802.2988 [hep-ex]

  43. [43]

    Aaij et al

    R. Aaij et al. (LHCb), Measurement of theCPasymme- try inB − →D − s D0 andB − →D −D0 decays, JHEP05, 160, arXiv:1803.10990 [hep-ex]

  44. [44]

    Aaij et al

    R. Aaij et al. (LHCb), Measurement of CP asymmetries in B 0 →D − s D+ and B 0 s →D + s D− decays, (2026), arX- iv:2603.28132 [hep-ex]

  45. [45]

    Aaij et al

    R. Aaij et al. (LHCb), Measurement of CP violation in B0 →D +D− andB 0 s →D + s D− s decays, JHEP01, 061, arXiv:2409.03009 [hep-ex]

  46. [46]

    M. B. Wise, Chiral perturbation theory for hadrons con- taining a heavy quark, Phys. Rev. D45, R2188 (1992)

  47. [47]

    M. T. AlFiky, F. Gabbiani, and A. A. Petrov, X(3872): Hadronic molecules in effective field theory, Phys. Lett. B640, 238 (2006), arXiv:hep-ph/0506141

  48. [48]

    Ji, X.-K

    T. Ji, X.-K. Dong, M. Albaladejo, M.-L. Du, F.-K. Guo, and J. Nieves, Establishing the heavy quark spin and light flavor molecular multiplets of the X(3872), Zc(3900), and X(3960), Phys. Rev. D106, 094002 (2022), arXiv:2207.08563 [hep-ph]

  49. [49]

    Guo and U.-G

    F.-K. Guo and U.-G. Meissner, Where is theχ c0(2P)?, Phys. Rev. D86, 091501 (2012), arXiv:1208.1134 [hep- ph]

  50. [50]

    Gomez Nicola and J

    A. Gomez Nicola and J. R. Pelaez, Meson meson scat- tering within one loop chiral perturbation theory and its unitarization, Phys. Rev. D65, 054009 (2002), arXiv:hep-ph/0109056

  51. [51]

    Zhang, X.-W

    L. Zhang, X.-W. Kang, and X.-H. Guo, Composite nature ofZ b states from data analysis, Eur. Phys. J. C82, 375 (2022), arXiv:2203.02301 [hep-ph]. 7 SUPPLEMENTAL MATERIALS Appendix A: Topology Diagrams The current-current operators,O 1 andO 2, have the same chiral structure, but different color structures, which will be reflected in corresponding Wilson coe...

  52. [52]

    +O(· · ·), ∆2 =f(∆ ′ 12 + ∆′

  53. [53]

    +O(· · ·), ∆3 =f(δ p + ∆3 12 + ∆3 32 −∆ 3

  54. [54]

    Appendix D: Bootstrap The errors of the fitted parameters are derived by making bootstrap from three fit schemes, respectively

    +O(· · ·), (C8) for each channel. Appendix D: Bootstrap The errors of the fitted parameters are derived by making bootstrap from three fit schemes, respectively. From the experimental data within the experimental error range, 10,000 sample sets are generated for each fit scheme and three fit schemes are performed to derive 10,000 sets of parameters for th...