A phenomenological analysis on isospin-violating decay of X(3872)

In a molecular scenario, we investigate the isospin-breaking hidden charm decay processes of X(3872), i.e., X(3872)→π+π-J/ψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$X(3872) \rightarrow \pi ^+ \pi ^- J/\psi $$\end{document}, X(3872)→π+π-π0J/ψ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$X(3872) \rightarrow \pi ^+ \pi ^- \pi ^0 J/\psi $$\end{document}, and X(3872)→π0χcJ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$X(3872)\rightarrow \pi ^0\chi _{cJ}$$\end{document}. We assume that the source of the strong isospin violation comes from the different coupling strengths of X(3872) to its charged components D∗+D-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^{*+} D^-$$\end{document} and neutral components D∗0D¯0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^{*0 } {\bar{D}}^0$$\end{document} as well as the interference between the charged meson loops and neutral meson loops. The former effect could fix our parameters by using the measurement of the ratio Γ[X(3872)→π+π-π0J/ψ]/Γ[X(3872)→π+π-J/ψ]\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma [X(3872) \rightarrow \pi ^+ \pi ^- \pi ^0 J/\psi ]/\Gamma [X(3872) \rightarrow \pi ^+ \pi ^- J/\psi ]$$\end{document}. With the determined parameter range, we find that the estimated ratio Γ[X(3872)→π0χc1/Γ[X(3872)→π+π-J/ψ]\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma [X(3872) \rightarrow \pi ^0 \chi _{c1}/\Gamma [X(3872) \rightarrow \pi ^+ \pi ^- J/\psi ]$$\end{document} is well consistent with the experimental measurement from the BESIII collaboration. Moreover, the partial width ratio of π0χcJ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pi ^0 \chi _{cJ}$$\end{document} for J=0,1,2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J=0,1,2$$\end{document} is estimated to be 1.77-1.65:1:1.09-1.43\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.77{-}1.65:1:1.09{-}1.43$$\end{document}, which could be tested by further precise measurements of BESIII and Belle II.

Beside the resonance parameters, the experimental measurements indicated that the π + π − invariant mass for X (3872) → π + π − J/ψ concentrates near the upper kinematic boundary, which corresponds to the ρ meson mass [1]. As for X (3872) → π + π − π 0 J/ψ, the π + π − π 0 invariant mass distribution has a strong peak between 750 MeV and the kinematic limit of 775 MeV, suggesting that the process is dominated by the sub-threshold decay X (3872) → ω J/ψ. The ratio of the branching fractions of π + π − J/ψ and π + π − π 0 J/ψ is determined to be, The large isospin violation implied by the almost equality of the branching fractions of ω J/ψ and ρ J/ψ channels further makes the nature of X (3872) complicated and confusing.
To date, the nature of X (3872) still remains unclear. Besides the mass spectrum, the investigations of decay behaviors are also crucial to understand the properties of X (3872). The isospin breaking effects of X (3872) were studied in Ref. [120], where X (3872) was considered as a dynamically generated state and the coupling strengths of X (3872)D * + D − and X (3872)D * 0D0 were assumed to be the same. Under 2 3 P 1 assignment, the decay channel X (3872) → ρ/ω + J/ψ was estimated via intermediate charmed meson loops [101]. Using a phenomenological Lagrangian approach, the authors studied radiative decays to J/ψ/ψ(2S) with the X (3872) being a composite state containing both D 0D * 0 molecule and a cc component [78,121], and hidden charm and radiative decays of X (3872) were investigated with the X (3872) being a composite state comprised of the dominant molecular D 0D * 0 component and other hadronic pairs, which could be D ± D * ∓ and J/ψω/ρ [122]. The final state interaction effects of the hidden charm decay of X (3872) were examined in Ref. [75], they found that the FSI contribution to X (3872) → J/ψρ is tiny. Assuming the decays of X (3872) through ρ J/ψ and ω J/ψ, the authors in Ref. [68] calculated the decay rates of X (3872) → π + π − J/ψ and X (3872) → π + π − π 0 J/ψ.
Besides π + π − J/ψ and π + π − π 0 J/ψ, the pionic transition from X (3872) to χ cJ was predicted in Refs. [77,123], in which it is found that the ratio of different transitions with different angular momentum J was sensitive to the inner structure of X (3872). In 2019, the BESIII Collaboration searched for the process e + e − → γ X (3872) by using the collision data with center-of-mass energies between 4.15 and 4.30 GeV and a new decay mode, χ c1 π 0 , of X (3872) was observed with a statistical significance of more than 5σ but no significant X (3872) signal was observed in the invariant mass distributions of π 0 χ c0,2 . The ratios of the branching ratios of X (3872) → π 0 χ cJ for J = 0, 1, 2 and X (3872) → π + π − J/ψ were measured to be [38], where the numbers in parentheses for J = 0 and 2 are the upper limits in 90% C.L. [38]. Later, the Belle collaboration searched for X (3872) in B + → χ c1 π 0 K + decay, and the ratio was measured to be B[X → π 0 χ c1 ]/B[X → π + π − J/ψ] < 0.97 at 90% C.L. [124]. The upper limit of the ratio measured by the Belle collaboration does not contradict the BESIII data [38]. The experimental measurements of the ratios in Eqs. (1) and (2) imply a strong isospin violation. The explanation of this fact is important for revealing the nature of X (3872). In the present work, we attempt to hunt for the source of the isospin violation in the molecular scenario by assuming that X (3872) is an S-wave molecule with J PC = 1 ++ given by the superposition of D 0D * 0 and D ± D * ∓ hadronic configurations. The fundamental source of the isospin violation is the mass difference of up and down quarks. Specific to the present discussed issue, the concrete manifestation is the mass difference of charged and neutral charmed mesons, which leads to the different coupling strengths of X (3872)D * 0D0 and X (3872)D * + D − . This coupling strength difference in part provides the source of the isospin violation in the decays of X (3872). As a molecular state, the hidden charm decays of X (3872) occur via the charmed meson loops, where the interferences between the charged and neutral meson loops provide another important source of the isospin violation. In the present work, we consider these two sources of isospin violation, the uncertainties of the former one, i.e. the different coupling strengths, can be determined by the ratio B[X → π + π − π 0 J/ψ]/B[X → π + π − J/ψ], and then with the fixed parameters, we can further estimate the ratios Comparing the present estimation for J = 1 with the BESIII data can also check the present model's reasonability. As for J = 0 and J = 2, the present estimations can narrow down the ratios' range, which could be tested by further measurements.
This paper organized as follows: After introduction, we present the model used in the present estimations of X (3872) → ρ/ω J/ψ and X (3872) → χ cJ π 0 . The numerical results and discussions are presented in Sects. 3, and 4 is devoted to a short summary.

hidden charm decay of X (3872)
As discussed above, we assume that the coupling strengths of X (3872)D * 0D0 and X (3872)D * + D − are different. The effective coupling of X (3872) with its components can be, where g X is the coupling constant, θ is a phase angle describing the proportion of neutral and charged constituents. It should be mention that from a more fundamental quark level point of view, the difference in coupling strength should be dynamically generated from u, d quark-antiquark pair generating processes and also from the wave function difference of the charged and neutral DD * . Furthermore, from the dispersion discussion in Refs. [101,125], the different coupling strength comes also from the dispersion integrals where the thresholds are different for the charged and neutral channels. From the phenomenological point of view, we can parameterize the coupling strength and the isospin breaking effects into a common factor g X and a phase angle θ in the effective Lagrangian as shown in Eq. (3).
Moreover, the distributions of the components, i.e., DD * , in the molecular state could be described by a wave function, which would then be integrated in the Feynman diagram calculations and affect the magnitude of partial widths. In the present work, we mainly focus on the ratios of the partial widths as given in Eqs. (1), (2), the form factor appears in both the numerators and denominators. As discussed in Ref. [78], the estimated ratio in the nonlocal case may not be too much different from the local one. Then the simple parameterization in Eq. (3) could be a reasonable approximation in estimating the order of magnitude of the ratio.
In the present work, the hidden charm decay processes of X (3872) occur via charmed meson loops, i.e., the charmo- In the present work, all these diagrams are estimated in hadronic level and all the involved interactions are depicted by effective Lagrangians. In heavy quark limit, one can con-struct the effective Lagrangian for charmonium and charmed mesons, which are [126][127][128] The coupling between light meson and charmed mesons could be obtained based on the heavy quark limit and chiral symmetry, which are [128][129][130] With the Lagrangian listed above, we can obtain the decay amplitude corresponding to X (3872) → ρ J/ψ, ω J/ψ, and π 0 χ cJ with J = 0, 1, 2. For brevity, we collect all the amplitudes corresponding to diagrams in Figs. 1 and 2 in Appendix. A and leave the coupling constants to be discussed in the following section.

Numerical results and discussion
Since mass difference of X (3872) and D * 0D0 is very tiny, the coupling constants g X are very sensitive to the mass of X (3872). Thus, in the present work, we mainly focus on the ratios of the hidden charm decay channels, which are independent on the coupling constants g X . Moreover, the involved charmonia in the present estimation are J/ψ and χ cJ . In the heavy quark limit, the coupling constants of the involved charmonia and charmed mesons can be related to the gauge couplings g 1 and g 2 by, f ψ = 426MeV and f χ c0 = 510MeV are the J/ψ and χ c0 decay constants [127], respectively. In the heavy quark and chiral limits, the charmed meson couplings to the light vector and pesudoscalar mesons have the following relationship [128,130], where the parameter g V = m ρ / f π with f π = 132 MeV being the pion decay constant and β = 0.9 [128]. By matching the form factor obtains from the light cone sum rule with that calculated from lattice QCD, one obtained the parameters λ = 0.56 GeV −1 and g = 0.59 [131].
In the amplitudes, the form factors should be considered to depict the inner structures and off shell effects of the charmed mesons in the loop. However, the mass of X (3872) is very close to the thresholds of D * D, which indicates that the components of X (3872), i.e., the charmed mesons connected to X (3872) in Figs. 1 and 2, are almost on shell. Therefore, we introduce only one form factor in a monopole form to depict the inner structure and the off-shell effects of the exchanged charmed meson, which is [47,130,[132][133][134], where the parameter can be further reparameterized as D ( * ) = m D ( * ) + α QCD with QCD = 0.22 GeV and m D ( * ) is the mass of the exchanged meson. The model parameter α should be of order of unity [47, [132][133][134], but its concrete value cannot be estimated by the first principle. In practice, the value of α is usually determined by comparing theoretical estimates with the corresponding experimental measurements.

Summary
In the present work, we have investigated the decay behaviors of X (3872) → π + π − J/ψ, π + π − π 0 J/ψ, and π 0 χ cJ in a molecular scenario and tried to understand isospin violations in X (3872) hidden charm decays. The fundamental source of the isospin violation has been shown in two dif-ferent aspects. The mass difference of the charged and neutral charmed mesons leads to different coupling strengths of X (3872)D * 0D0 and X (3872)D * + D − , which in part provides the source of isospin violation in the decay of X (3872). Another important source of isospin violation is the interference between the charged and neutral meson loops.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: This is a theoretical study and no external data are associated with this work.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .