text:
# The gluonic 𝐵 and 𝐽/𝜓 decays into the 𝜂' meson

## Abstract

The inclusive and exclusive $`B`$ decays into the $`\eta ^{}`$ meson plus others are investigated in a model based on the QCD anomaly. The invariant mass distribution is discussed. The constraint of the effective coupling is obtained from the data of the exclusive decay modes. The branching ratio of $`J/\psi \eta ^{}\eta \gamma `$ is predicted to be as large as $`5.4\times 10^5`$, which can be tested in the forthcoming CLEO-c experiments.

preprint: LMU 20-03preprint: hep-ph/0301038

1. Introduction.
A few years ago, unexpected large branching ratios of $`B`$ decaying into final states with an $`\eta ^{}`$ meson such as $`B\eta ^{}X_s`$ and $`B\eta ^{}K`$ were observed by the CLEO collaborationBrowder et al. (1998); Behrens et al. (1998) and recently confirmed by BaBarAubert et al. (2001) and BelleBel . This stimulated many theoretical activities in understanding the special role of the $`\eta ^{}`$ meson in B decays. As the contribution of traditional four quark operators from the effective Hamiltonian in the standard model (SM) is far below the dataDatta et al. (1998); Ali and Greub (1998), various exotic mechanisms were introduced such as a large coupling between the gluon and $`\eta ^{}`$ through the QCD anomalyAtwood and Soni (1997); Hou and Tseng (1998); Fritzsch (1997); Gronau and Rosner (1996); Dighe et al. (1997), intrinsic $`\overline{c}c`$ content inside $`\eta ^{}`$Halperin and Zhitnitsky (1998); Yuan and Chao (1997) and positive interference between several contributions in the SMLipkin (1991); Beneke and Neubert (2002) et.al. The large contribution arising from new physics was also discussedKagan and Petrov (1997); Hou and Tseng (1998). Among those theoretical efforts, the possible enhancement from the QCD anomaly is of particular interest since it is well known that the $`\eta ^{}`$ meson plays an very special role in the dynamics of low energy QCD Fritzsch and Minkowski (1975).

The mass eigenstates $`\eta ^{}`$ and $`\eta `$ are related to flavor octet and singlet states $`\eta _8`$, $`\eta _0`$ through a mixing matrix:

$`\eta =\eta __8\mathrm{cos}\theta \eta __0\mathrm{sin}\theta ,\eta ^{}=\eta __8\mathrm{cos}\theta +\eta __0\mathrm{sin}\theta ,`$ (1)
where $`\theta `$ is the mixing angle and $`\eta _8`$, $`\eta _0`$ have the flavor content: $`\eta _8=\frac{1}{\sqrt{6}}(\overline{u}u+\overline{d}d2\overline{s}s)`$ and $`\eta _0=\frac{1}{\sqrt{3}}(\overline{u}u+\overline{d}d+\overline{s}s)`$. The associated axial currents are $`j_5^{\mu 8}=\frac{1}{\sqrt{6}}(\overline{u}\gamma ^\mu \gamma _5u+\overline{d}\gamma ^\mu \gamma _5d2\overline{s}\gamma ^\mu \gamma _5s)`$ and $`j_5^{\mu 0}=\frac{1}{\sqrt{3}}(\overline{u}\gamma ^\mu \gamma _5u+\overline{d}\gamma ^\mu \gamma _5d+\overline{s}\gamma ^\mu \gamma _5s)`$ respectively. Through the QCD anomaly, the divergence of the flavor singlet axial current is non-zero and is given by

$`_\mu j_5^{\mu 0}={\displaystyle \frac{1}{\sqrt{3}}}\left(2i{\displaystyle \underset{q=u,d,s}{}}m_q\overline{q}\gamma _5q+{\displaystyle \frac{3\alpha _s}{4\pi }}G_{\mu \nu }\stackrel{~}{G}^{\mu \nu }\right)`$ (2)
where $`\stackrel{~}{G}^{\mu \nu }=\frac{1}{2}\epsilon ^{\mu \nu \rho \sigma }G_{\rho \sigma }`$ is the dual of $`G_{\mu \nu }`$. This breaks the global $`U(1)_A`$ symmetry for massless particles and makes the flavor singlet state $`\eta _0`$ evade to be a Goldstone boson of chiral $`SU(3)_LSU(3)_R`$ symmetry. The QCD anomaly gives its main contribution to the large mass of $`\eta ^{}(m_\eta ^{}=0.958\text{GeV})`$ which is much heavier than the other flavor octet states such as $`\pi `$, $`K`$s and suggests a large gluon content in $`\eta ^{}`$.

The QCD anomaly indicates a strong coupling between $`\eta ^{}`$ and the gluon field. It is then natural to understand the large branching ratio of $`B\eta ^{}X_s`$ through the QCD anomaly. In the literatures there are two different ways to handle this problem. The one is through two-body decay process $`bs\eta ^{}`$ from some effective Hamiltonian due to QCD anomalyFritzsch (1997). The other one is through three-body process $`bsg\eta ^{}`$Atwood and Soni (1997); Hou and Tseng (1998). In the first step decay, the $`b`$ quark decays into the $`s`$ quark and a virtual gluon $`g^{}`$, then $`g^{}`$ decays into $`\eta ^{}`$ and a on shell gluon $`g`$. This model has some advantages in explaining the spectrum of invariant mass distribution of recoiling hadrons. However, the effective $`gg^{}\eta ^{}`$ vertex seems to be too small from various approaches Hou and Tseng (1998); Muta and Yang (2000); Ali and Parkhomenko (2002); Artuso et al. (2002). In both of the approaches, the effective coupling between $`\eta ^{}`$ and gluon may contain complicate non-perturpative quark-gluon interactions. It is then better to treat it as a free phenomenological parameter rather than to evaluate it from perturbative calculationsFritzsch (1997).

In this note we focus on the phenomenological analysis of the first mechanism. The effective Lagrangian in this model is given by Fritzsch (1997)

$`H_{eff}=a\alpha _sG_F\overline{s}_Lb_RG_{\mu \nu }\stackrel{~}{G}^{\mu \nu }+\text{h.c}`$ (3)
where $`\alpha _s`$ and $`G_F`$ are the strong coupling constant and Fermi constant, $`a`$ is the effective coupling. From this effective Hamiltonian, the decay $`B\eta ^{}X_s`$ arises from the subprocess $`bs\eta ^{}`$. The evaluation of matrix elements is straightforward:

$`s\eta ^{}\left|H_{eff}\right|b=a\alpha _sG_Fs\left|\overline{s}_Lb_R\right|b\eta ^{}\left|G_{\mu \nu }\stackrel{~}{G}^{\mu \nu }\right|0`$ (4)
Applying relation $`0\left|j_5^{\mu 8(0)}\right|\eta _{8(0)}=if_{8(0)}P^\mu `$ to the divergences of both flavor singlet and octet axial currents and ignoring small $`u,d`$ quarks masses, the matrix elements $`0\left|G_{\mu \nu }\stackrel{~}{G}^{\mu \nu }\right|\eta ^{}`$ can be rewritten as

$`0\left|\alpha _sG_{\mu \nu }\stackrel{~}{G}^{\mu \nu }\right|\eta ^{}`$ $`={\displaystyle \frac{4\pi }{3}}\sqrt{{\displaystyle \frac{3}{2}}}m_\eta ^{}^2(f_8\mathrm{sin}\theta +\sqrt{2}f_0\mathrm{cos}\theta )`$
$`0\left|\alpha _sG_{\mu \nu }\stackrel{~}{G}^{\mu \nu }\right|\eta `$ $`={\displaystyle \frac{4\pi }{3}}\sqrt{{\displaystyle \frac{3}{2}}}m_\eta ^2(f_8\mathrm{cos}\theta \sqrt{2}f_0\mathrm{sin}\theta )`$ (5)
In the $`b`$ quark rest frame, the decay branching ratio is given by

$`Br(B\eta ^{}X_s)={\displaystyle \frac{\pi }{12}}\tau _Ba^2G_F^2m_\eta ^{}^4(f_8\mathrm{sin}\theta +\sqrt{2}f_0\mathrm{cos}\theta )^2{\displaystyle \frac{(m_b^2m_\eta ^{}^2)^2}{m_b^3}}`$ (6)
$`Br(B\eta X_s)={\displaystyle \frac{\pi }{12}}\tau _Ba^2G_F^2m_\eta ^4(f_8\mathrm{cos}\theta \sqrt{2}f_0\mathrm{sin}\theta )^2{\displaystyle \frac{(m_b^2m_\eta ^2)^2}{m_b^3}}`$ (7)
where $`\tau _B`$ is the lifetime of $`B`$ meson.

2. Recoil mass distribution

For two-body like subprocess such as $`b\eta ^{}s`$, the invariant mass is directly related to the energy $`E_\eta ^{}`$ of the $`\eta ^{}`$ meson through the relation: $`m_X^2=m_B^2+m_\eta ^{}^22m_BE_\eta ^{}`$, where $`m_B`$ and $`m_\eta ^{}`$ are the masses of $`B`$ and $`\eta ^{}`$ meson. The small $`s`$ quark mass has been ignored. In the two-body decay of $`b\eta ^{}s`$, the energy of $`\eta ^{}`$ is fixed from energy-momentum conservation. A typical value of the pole mass $`m_b=4.8`$ GeV will lead to a narrow peak with the invariant mass of $`m_{X_s}1.5`$ GeV. This seems to be disfavored by the current data since the experiment reported a peak at about 2 GeV with a relative wide width in the recoil mass distribution Browder et al. (1998); Aubert et al. (2001).

However, the above estimation may be too naive. Note that in the two-body decay process, the exact distribution of the recoil mass strongly depends on the wave function of the $`B`$ meson which is theoretically hard to estimate. It is too early to draw the conclusion that the current data already disfavored all the two-body models.

To illustrate the non-perturbative bound state effects here we adopt a simple model proposed by Altarelli et.al. Altarelli et al. (1982) a number of years ago which is based on the Fermi motion of the $`b`$ quark inside $`B`$ meson. The basic idea of this model is that the Fermi motion of the $`b`$ quark and the spectator quark $`q`$ in the B meson make them have back-to-back relative three-momenta $`𝐩`$ in the $`B`$ rest frame. The momentum is assumed to obey a Gaussian distribution as follows

$`\varphi (p)={\displaystyle \frac{4}{\sqrt{\pi }p__F^3}}e^{p^2/p__F^2},p=\left|𝐩\right|`$ (8)
where $`\varphi (p)`$ is normalized as $`_0^{\mathrm{}}\varphi (p)p^2𝑑p=1`$. The mean value of $`p`$ is $`<p>=\frac{3}{2}p__F`$. In this model the spectator quark $`q`$ is on always handled as on shell while the $`b`$ quark is treated as off-shell. Through energy-momentum conservation, the effective mass $`W`$ of the b quark is determined as

$`W^2`$ $`=m_B^2+m_q^22m_B\sqrt{m_q^2+p^2}`$ (9)
and the energy of the $`b`$ quark is $`E_W=\sqrt{W^2+p^2}`$. Here one parameter $`p__F`$ is introduced which specifies both the distribution width and the mean value. As $`p__F`$ is linked to the average energy of $`b`$ quark inside the $`B`$ meson, in principle it can be calculated from theories based on non-perturbative methods or from some models. For example, the calculations from QCD sum rule Bagan et al. (1995) give $`p__F=0.58\pm 0.06`$ GeV, and the value from relativistic quark model Hwang et al. (1995) is $`0.54\pm 0.16`$ GeV. The value of $`p__F`$ can also be extracted directly from the data. A fit to $`BX_s\gamma `$ photon energy spectrum give a value of about $`0.45`$ GeVAli and Greub (1995) while the fits to semi leptonic $`B`$ decays and $`BJ/\psi X`$ give a value of 0.57 GeV Palmer et al. (1997); Hwang et al. (1996). Thus the value of $`p__F`$ is likely to lie in the range $`0.4p__F0.6`$ GeV. After including the Fermi motion, the differential decay width $`d\mathrm{\Gamma }(m_b)/dm_X`$, should be replaced by

$`{\displaystyle \frac{d\mathrm{\Gamma }}{dm_X}}={\displaystyle _0^{p_{max}}}𝑑p\varphi (p)p^2{\displaystyle \frac{d\mathrm{\Gamma }(W)}{dm_X}}`$ (10)
where $`p_{max}`$ is the allowed maximum value of $`p`$, $`d\mathrm{\Gamma }(W)/dm_X`$ is the differential decay rate in the $`B`$ meson rest frame, which is linked to the one in $`b`$ quark rest frame through a Lorentz boostPalmer et al. (1997); Du and Yang (1998) In Fig.1 the invariant mass distribution is generated in this model with different values of $`p__F`$. Here we use $`\frac{1}{\mathrm{\Gamma }}\frac{d\mathrm{\Gamma }}{dm_X}`$ which is normalized to unity and independent on the value of $`a`$.

It can be clearly seen that the $`p__F`$ dependence is rather strong. The peak of the distributions significantly shifts from around $`1.4`$ GeV ( for $`p__F=0.4`$ GeV) to $`1.8`$ GeV (for $`p__F=0.6`$ GeV ). Considering the considerable uncertainties in both theory and experiment data, there is no significant disagreement in the recoil mass distribution of $`B\eta ^{}X_s`$.

3. Bound on $`a`$ from inclusive and exclusive $`B`$ decays
The value of $`a`$ could be constrained from the exclusive decay modes $`B\eta ^{}(\eta )K^{()}`$. Note that although predictions of the standard effective Hamiltonian approach are too low to account for the data of inclusive decay modes, the disagreement in the exclusive decay modes are smallerAli and Greub (1998); Beneke and Neubert (2002). Furthermore, the effective Hamiltonian approach can reproduce the correct patterns of $`Br(B\eta ^{}K)Br(B\eta K)`$ and $`Br(B\eta ^{}K^{})Br(B\eta K^{})`$ which is observed in the experiments. It implies that in exclusive decays modes, it may still play an important role, and the interference between different contributions may also be importantChiang and Rosner (2002); Gronau and Rosner (2000).

Nevertheless, by saturating the current data of exclusive decays, the upper bound of the parameter $`a`$ can still be obtained. The decay amplitudes of decay modes $`B\eta ^{}K^{()}`$ and $`B\eta K^{()}`$ in this model read

$`(B^{\pm ,0}\eta ^{}K^{\pm ,0})`$ $`=aG_F{\displaystyle \frac{4\pi }{3}}\sqrt{{\displaystyle \frac{3}{2}}}m_\eta ^{}^2(f_8\mathrm{sin}\theta +\sqrt{2}f_0\mathrm{cos}\theta ){\displaystyle \frac{m_B^2m_K^2}{2(m_bm_s)}}F_0^{BK}(m_\eta ^{}^2)`$
$`(B^{\pm (0)}\eta K^{\pm (0)})`$ $`=aG_F{\displaystyle \frac{4\pi }{3}}\sqrt{{\displaystyle \frac{3}{2}}}m_\eta ^2(f_8\mathrm{cos}\theta \sqrt{2}f_0\mathrm{sin}\theta ){\displaystyle \frac{m_B^2m_K^2}{2(m_bm_s)}}F_0^{BK}(m_\eta ^2)`$
$`(B^{\pm (0)}\eta ^{}K^{\pm (0)})`$ $`=aG_F{\displaystyle \frac{4\pi }{3}}\sqrt{{\displaystyle \frac{3}{2}}}m_\eta ^{}^2(f_8\mathrm{sin}\theta +\sqrt{2}f_0\mathrm{cos}\theta ){\displaystyle \frac{\left|P_{\eta ^{}K^{}}\right|m_B}{(m_b+m_s)}}A_0(m_\eta ^{}^2)`$
$`(B^{\pm (0)}\eta K^{\pm (0)})`$ $`=aG_F{\displaystyle \frac{4\pi }{3}}\sqrt{{\displaystyle \frac{3}{2}}}m_\eta ^2(f_8\mathrm{cos}\theta \sqrt{2}f_0\mathrm{sin}\theta ){\displaystyle \frac{\left|P_{\eta K^{}}\right|m_B}{(m_b+m_s)}}A_0(m_\eta ^2)`$ (11)
where $`\left|P_{\eta ^{}K^{}}\right|\left|P_{\eta K^{}}\right|\frac{1}{2}m_B`$. $`F_0^{BK}(q^2)`$ and $`A_0(q^2)`$ are the form factors for $`BK`$ and $`BK^{}`$ transition with momentum transfer $`q^2`$. The value of $`m_b`$ is taken to be the effective one. i.e $`m_b^2W^2\varphi (p)p^2𝑑p`$. In the calculations, we take $`m_b=4.65`$ GeV which corresponds to $`p__F=0.5`$GeV and $`m_q=0.15`$GeV.

The corresponding branching ratio can be evaluated through the following relation

$`Br={\displaystyle \frac{\tau _B\left|P\right|}{8\pi m_B^2}}\left|\right|^2`$ (12)
where $`\left|P\right|`$ is the momentum of one of the final state mesons in $`B`$ rest frame.

It is useful to define two kind of ratios which are independent of the parameter $`a`$:

1) The ratio between $`B\eta ^{}X`$ and $`B\eta X(X=X_s,KorK^{})`$. This ratio is independent of the value of $`a`$ and only sensitive to the $`\eta ^{}\eta `$ mixing. In this model one findsHe et al. (1998)

$`R`$ $`{\displaystyle \frac{Br(B\eta X_s)}{Br(B\eta ^{}X_s)}}={\displaystyle \frac{Br(B\eta K)}{Br(B\eta ^{}K)}}={\displaystyle \frac{Br(B\eta K^{})}{Br(B\eta ^{}K^{})}}`$
$`={\displaystyle \frac{m_\eta ^4}{m_\eta ^{}^4}}\left({\displaystyle \frac{f__8\mathrm{cos}\theta \sqrt{2}f__0\mathrm{sin}\theta }{f__8\mathrm{sin}\theta +\sqrt{2}f__0\mathrm{cos}\theta }}\right)^2`$ (13)
In the following numerical calculations we take $`\theta =17^{}`$ and $`f_8=f_0=1.06f_\pi `$ Ball et al. (1996) as an illustration. This leads to a value of $`R0.16`$. Considering the CLEO data of $`R(0.10.8)`$ Behrens et al. (1998), it follows that with the constraints from $`B\eta K`$, this model can account for at most $`60\%`$ of the observed $`B\eta ^{}K`$ branching ratio. Note that the exact value of $`R`$ may vary with different sets of parameters $`\theta `$, $`f_8`$ and $`f_0`$; the constraints from $`R`$ are only an order of magnitude estimate.

2) The ratio between $`BPK^{}`$ and $`BPK(P=\eta ^{}or\eta )`$. In this model it is independent of $`both`$ the value of $`a`$ and details of $`\eta ^{}\eta `$ mixing.

$`R^{}{\displaystyle \frac{Br(B\eta ^{}K^{})}{Br(B\eta ^{}K)}}={\displaystyle \frac{(m_B^2+m_K^{}^2m_\eta ^2)^24m_B^2m_K^{}^2}{(m_B^2m_K^2)^2}}\left({\displaystyle \frac{A_0(m_\eta ^{}^2)}{F_0^{BK}(m_\eta ^{}^2)}}\right)^2`$ (14)
The values of $`F_0`$ and $`A_0`$ in the BSW model Bauer et al. (1987) are $`F_0=0.38,A_0=0.32`$ which corresponds to $`R=0.84`$, while from light cone QCD sum rule Ali et al. (1994); Ball (1998) $`F_0=0.35\pm 0.05,A_0=0.39\pm 0.1`$ and $`R=1.1\pm 0.3`$. Thus if this model gives the dominant contribution to these modes, the value of $`R`$ should be around 1. However, the current data gives a value of $`R^{}(0.50.4)`$ Behrens et al. (1998). This is a more clear and stronger constraint than the one from $`R`$. With the observed small value of $`R^{}`$, this model can explain at most half of the branching ratio of $`B\eta ^{}(\eta )K`$ and therefore is not the dominant mechanism of these processes. In Fig.2(c-f) the numerical results of branching ratios as a function of the effective coupling $`a`$ is given and compared with the data. As some inclusive decay modes have not yet been observed by the Babar and Belle collaborations, only the CLEO data are used in the numerical evaluations. It can be seen from the figure that the data of exclusive decay modes $`B\eta ^{}K^{}`$ and $`B\eta K`$ impose strong constraints on the effective coupling. With these constraints, the maximum value of $`a`$ lies in the range:

$`a(89)\times 10^3\text{GeV}\text{-1}`$ (15)
From Eq.(6) and (7), the branching ratio of inclusive decays $`B\eta ^{}(\eta )X_s`$ as a function of $`a`$ is plotted in Fig.2(a-b) and compared with the data. In the decay $`B\eta ^{}X_s`$ the acceptance cut effects is taken into account, which leads to a 19$`\%`$ reduction from the calculation in Eq.(6). Given the upper bound of $`a`$ in Eq.(15) this model can still successfully reproduce the $`B\eta ^{}X_s`$ branching ratio within 1$`\sigma `$ range.

4. Prediction of radiative decay $`J/\psi \eta ^{}\eta \gamma `$
From the effective Hamiltonian in Eq.(3), this model can also contribute to the radiative $`J/\psi `$ decays into $`\eta ^{}`$. Using relation Eq.(The gluonic $`B`$ and $`J/\psi `$ decays into the $`\eta ^{}`$ meson), the ratio between $`J/\psi \eta ^{}\gamma `$ and $`J/\psi \eta \gamma `$ can be predicted and found to be the same as in Ref.Ball et al. (1996)

$`{\displaystyle \frac{\mathrm{\Gamma }(J/\psi \eta ^{}\gamma )}{\mathrm{\Gamma }(J/\psi \eta \gamma )}}=\left|{\displaystyle \frac{0\left|G\stackrel{~}{G}\right|\eta ^{}}{0\left|G\stackrel{~}{G}\right|\eta }}\right|^2{\displaystyle \frac{(1m_\eta ^{}^2)^3}{(1m_\eta ^2)^3}}`$ (16)
which is in good agreement with the data.

Furthermore, given the value of the effective coupling $`a`$ the decay rate of $`J/\psi \eta ^{}\eta \gamma `$ can be predicted. To this end let us first define the ratio

$`r(\eta ^{})`$ $`=`$ $`\mathrm{\Gamma }(B\eta ^{}X_s)/\mathrm{\Gamma }(Bg^{}X_s)`$ (17)
which can be understood as the size of $`bs\eta ^{}`$ relative to $`bsg`$. Taking $`\mathrm{\Gamma }(Bg^{}X_s)1\%`$ and $`a0.008`$ GeV<sup>-1</sup> which comes from the bounds from exclusive decays as an example, one finds

$`r(\eta ^{})0.045`$ (18)
Note that the strong coupling constant in the effective Hamiltonian has been separated from the effective coupling $`a`$ and absorbed in the matrix element of $`0\left|\alpha _sG\stackrel{~}{G}\right|\eta ^{}(\eta )`$. It is expected that there is no significant running on the value of $`a`$ from energy scale $`m_B`$ to $`m_{J/\psi }`$.

Since the radiative decay of $`J/\psi \gamma X`$ is dominated by the process $`J/\psi g^{}g^{}\gamma `$, the branching ratio of $`J/\psi \eta ^{}\eta ^{}\gamma `$ can be simply estimated as

$`{\displaystyle \frac{Br(J/\psi \eta ^{}\eta ^{}\gamma )}{Br(J/\psi \gamma X)}}r(\eta ^{})^2.`$ (19)
The observation of the process $`J/\psi \gamma X`$ gives $`Br(J/\psi \gamma X)=(17.0\pm 2.0)\times 10^2`$. Thus taking $`r(\eta ^{})=0.045`$ the maximum branching ratio of $`J/\psi \eta ^{}\eta ^{}\gamma `$ is estimated as

$`Br(J/\psi \eta ^{}\eta ^{}\gamma )3.4\times 10^4`$ (20)
The decay rate of $`J/\psi \eta ^{}\eta \gamma `$ can be estimated as follows

$`{\displaystyle \frac{Br(J/\psi \eta ^{}\eta \gamma )}{Br(J/\psi \eta ^{}\eta ^{}\gamma )}}=R`$ (21)
Using the value of $`R=0.16`$ from Eq.(The gluonic $`B`$ and $`J/\psi `$ decays into the $`\eta ^{}`$ meson) one finds for the maximum branching ratio for $`J/\psi \eta ^{}\eta \gamma `$

$`Br(J/\psi \eta ^{}\eta \gamma )5.4\times 10^5`$ (22)
Considering the detection efficiency of $`\eta ^{}`$ is about a few percent ( through $`\eta ^{}\eta \gamma \gamma `$), it may be hard to find a signal of such a decay mode in BES due to limited statistics ( in BES $`5\times 10^7J/\psi `$ samples are collected). But in the forthcoming CLEO-c project $`1\times 10^9J/\psi `$ samples are planned to be produced. It is then promising to search for the signal and test the predictions from this model in the CLEO-c experiment.

###### Acknowledgements.

Y.F. Zhou acknowledges the support by the Alexander von Humboldt Foundation.