arxiv: v3 [hep-ph] 27 Jul 2012

Transcription

1 Revisiting the roduction of charmonium lus a light meson at PANDA Qing-Yong Lin,2,3 and Hu-Shan Xu,2 Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou , China 2 Research Center for Hadron and CSR Physics, Lanzhou University and Institute of Modern Physics of CAS, Lanzhou , China 3 Graduate University of Chinese Academy of Sciences, Beijing, 00049, China Xiang Liu,2 School of Physical Science and Technology, Lanzhou University, Lanzhou , China 2 Research Center for Hadron and CSR Physics, Lanzhou University and Institute of Modern Physics of CAS, Lanzhou , China (Dated: October, 208) arxiv: v3 [he-h] 27 Jul 202 In this work, we calculate the total cross sections and the center-of-mass frame angular distributions of the charmonium roduction lus a light meson by the low energy interaction. The results of π 0 Ψ with and without form factor (FF) indicate that the FF contribution in the calculation cannot be ignored. The obtained cross section of π 0 J/ψ with FF can fit the E760 data well. We also redict the total cross sections and the center-of-mass frame angular distributions of ωψ, which show that these hysical quantities are deendent on Pauli (g ω ) and Dirac (κ ω ) couling constants of the ω interaction. Thus, ωψ can be as the ideal channel to test the different theoretical values of g ω andκ ω. Alying the formulae of π 0 Ψ and ωψ, we redict the total cross sections of the ηψ and ρψ reactions. Our results show a common behavior of the charmonium roduction with a light meson by the interaction, where the total cross section of theη c roduction is the largest one among all discussed rocesses. The above observations can be directly tested at the forthcoming PANDA exeriment. PACS numbers: 4.40.Pq, 3.75.Cs, 3.75.Gx I. INTRODUCTION Since the discovery of J/ψ in 974 [, 2], more and more charmonia have been reorted by exeriment [3], which rovide an ideal latform to imrove our understanding of nonerturbative Quantum Chromodynamics (QCD) dynamics. In the ast decade, the observations of a series of charmoniumlike states named as XYZ state have stimulated extensive interest in studying higher charmonia for both theorists and exerimentalists. As one of the forthcoming exeriments relevant to the study of hadron hysics, Antiroton Annihilations at Darmstadt (PANDA) exeriment at the Facility for Antiroton and Ion Research (FAIR) can serve as the investigation of charmonium, which is also one of the main hysical aims of PANDA [4]. Thus, theoretically studying the charmonium roduction by the low energy interaction becomes an imortant and interesting research work, which can give valuable suggestions for the forthcoming PANDA exeriment. In Ref. [5], Gaillard and Maiani firstly calculated the differential cross section of the charmonium roduction accomanied by a soft ion in the low energy interaction, where two hadron-level diagrams were introduced by the Born aroximation. They indicated that the corresonding cross section is roortional to the artial decay width of charmonium decay into [5]. The authors of Ref. [6] further studied the cross sections of the charmonium (Ψ) roduction lus a light meson (m) by the Ψ+m rocesses, which can be related to the corresonding author Electronic address: qylin@imcas.ac.cn Electronic address: xiangliu@lzu.edu.cn measured artial decay width of charmonium decay into m [6]. By this aroach, the cross sections of J/ψ(ψ )m (m=π 0,η,ρ 0,ω,η,φ) are estimated. Among these redicted cross sections, the cross section of the π 0 J/ψ rocess can reach u to 300 b when taking s=3 GeV, which is close to the E760 data taking s= GeV [7]. Later, Barnes and Li develoed the initial state light meson emission model, which was alied to study the near threshold associated charmonium roduction rocess π 0 Ψ, whereψdenotesη c, J/ψ,ψ,χ c0,χ c [8]. By the initial state light meson emission model, they calculated the differential cross section and total cross sections of π 0 Ψ rocesses. In addition, dσ/dω, the center-of-mass frame unolarized angular distribution, was redicted for π 0 Ψ [8]. In Ref. [9], Barnes et al. further indicated that the cross section of π 0 J/ψ near threshold may be affected by the Pauli J/ψ couling, which will be an interesting research toic in PANDA [9]. By the initial state light meson emission model [8], π 0 Ψ occurs via a roton exchange between and. In Ref. [8], authors treated the m and J/ψ coulings as the oint-like interaction vertices. However, in reality we should consider the structure effect of the m and J/ψ interactions. Thus, for reflecting such structure effect, the form factor should be introduced in the m and J/ψ interaction vertices, which was listed as one of the future develoments of the initial state light meson emission model [8]. Along this way, in this work we revisit the roduction of charmonium lus a light meson in the low energy interaction by considering the contribution of form factor (FF) to these rocesses. The comarison of the results with and without including FF in the calculation can reveal the difference under two cases, which will be tested at the forthcoming PANDA exeriment. By this study, we can not only extract the hysical icture de-

2 2 icting the roduction of charmonium lus a light meson in the low energy, but also learn what kind of FF to be suitable to describe the structure effect of the m and J/ψ interactions. These rocesses discussed in this work also include the charmonium roduction with a light vector meson (ω, ρ), where the ω or ρ interaction relates to both Dirac and Pauli coulings. At resent, the couling constants of Dirac and Pauli coulings of ω and ρ are determined by some theoretical grous by different rocesses and different models [0 8]. Thus, we adot these determined couling constants in our calculation, which includes the total cross sections, the differential cross sections dσ/dω. These studies can serve as further exerimental test of these couling constants by the charmonium roduction with a light vector meson at PANDA. This work is organized as follows. After introduction, we resent the calculation of the roduction of charmonium lus a light meson in the low energy interaction. In Sec. III, the numerical results are given. The last section is the discussion and conclusion. II. THE PRODUCTION OF CHARMONIUM As deicted by Fig., the charmonium (Ψ) roduction lus a light meson (m) by the low energy interaction can occur via the transition of intoψ+m by exchanging a roton. Thus, there exist two hadron-level diagrams [5, 8] shown in Fig. if only considering the tree-level contributions. m Ψ + FIG. : The diagrams describing the roduction charmonium lus a light meson by the interaction. Here,Ψdenotes charmonium (η c, J/ψ,ψ,χ c0,χ c ) while m is light meson (π 0,η,ρ,ω). When deducing the corresonding roduction amlitude, we use effective Lagrangian aroach. The interaction vertices of m and Ψ include ig NNπ φγ 5 τ πφ, for m=π 0 ( L m = g NNω φγ µ φω µ κ ) ω φσ µν φf µν, for m=ω 4m and L Ψ = ig NNηc φγ 5 φη c, forψ=η c g NNχc0 φφχ c0, forψ=χ c0 g NNJ/ψ(ψ ) φγ µ φψ µ, forψ= J/ψ(ψ ) g NNχc φγ µ γ 5 φχ µ c, forψ=χ c whereφdenotes the field of nucleon and theτis Pauli matrix.σ µν = i[γ µ,γ ν ]/2 and F µν = µ ω ν ν ω µ. For the ω m Ψ interaction, there exists two indeendent couling constants, i.e, g ω g NNω andκ ω corresonding to Dirac and Pauli terms resectively. In the following, we illustrate the calculation of the charmonium roduction. Without introducing FF in each of interaction vertices in Fig., the resulting amlitude of π 0 Ψ is [8] (/ M π0 Ψ = g π g Ψ v ( 2, s 2 ) (Γ /k+m ) (t m 2 γ 5 +γ 5 (/k / 2 + m ) Γ (u m 2 )u (, s ), () where m is the mass of roton., 2 and k are the four momenta of roton, antiroton and the emitted light meson. s, t and u are the Mandelstam variables with the definitions s=( + 2 ) 2, t=( k) 2 and u=(k 2 ) 2. g π g π and g Ψ g Ψ denote the couling constants of the π and Ψ interactions, resectively. ν and u denote the sinors of antiroton and roton, resectively. In Eq. (),Γ is defined as γ 5, i, iγ µ ǫ µ J/ψ and iγ µγ 5 ǫ χ µ c corresonding to π 0 Ψ rocesses withψtaken asη c,χ c0, J/ψ(ψ ) andχ c, resectively. For the ωψ rocess discussed here, its roduction amlitude can be written as M ωψ = g ω g Ψ v ( 2, s 2 ) ( Γ 2 / /k m γ µ +γ µ u (, s )ǫµ+ i κ ωg ω g Ψ v ( 2, s 2 ) 2m σ µν k ν +σ µν k ν Γ 2 /k / 2 m (/ /k+m ) = g ω g Ψ v ( 2, s 2 ) Γ 2 (t m 2 ) γ µ +γ µ (/k / 2 + m ) (u m 2 ) Γ 2 u (, s )ǫµ +i κ ωg ω g Ψ v ( 2, s 2 ) 2m Γ 2 /k / 2 m ( Γ 2 / /k m ) u (, s )ǫ µ (/ Γ /k+m ) 2 σ µν k (t m 2 ν +σ µν (/k / k 2 + m ) ν Γ (u m 2 2 u (, s )ǫµ, (2) whereǫµ is the olarization four-vector of the emittedωmeson. The definition ofγ 2 is the same as that ofγ in calculating π 0 Ψ. Just discussed above, the FF contribution to π 0 Ψ and ωψ is not included in Eqs. () and (2). Since the m and Ψ vertices are not the oint-like interactions, we need to introduce FF in each of interaction vertices, where the FF not only reflects the structure effect of the interaction vertex but also lays an imortant role to comensate the off-shell effect of the exchanged roton. For comaring the results with and without FF, we calculate the π 0 Ψ and ωψ rocesses considering the FF contribution, where the corre- )

3 3 sonding amlitudes can be exressed as and M FF π 0 Ψ = g (/ πg Ψ v ( 2, s 2 ) [Γ /k+m ) γ 5 F 2 (q 2 (t m 2 t ) +γ 5 (/k / 2 + m ) ] Γ (u m 2 F 2 (q 2 u) u (, s ) (3) M FF ωψ (/ = g ω g Ψ v ( 2, s 2 ) Γ /k+m ) 2 (t m 2 ) γ µ F 2 (q 2 t )+γµ (/k / 2 + m ) (u m 2 ) Γ 2 F 2 (q 2 u ) u (, s )ǫµ + iκ ωg ω g Ψ 2m (/ v ( 2, s 2 ) [Γ /k+m ) 2 σ µν k (t m 2 ν F 2 (q 2 t )+σµν k ν (/k / 2 + m ) ] Γ (u m 2 2 F 2 (q 2 u) u (, s )ǫµ, (4) wheref 2 (q 2 i ) denotes the introduced monoole FF with the definitionf (q 2 i )=(Λ2 m 2 i )/(Λ2 q 2 i ). m i and q i are the mass and the four-momentum of the exchanged roton, resectively. Λ is a free arameter, which is exected to be around GeV. In the next section, we will discuss the value ofλadoted in our calculation. In the above exressions, q t and q u denote the four-momenta of the exchanged rotons of t-channel and u-channel of the mψ rocess just shown in Fig., resectively. In addition, the suerscrit FF is introduced for distinguishing the amlitudes with and without FF. As indicated in Ref. [20], the mψ rocesses may include transition distribution amlitudes [9], which generalize the form factors that we include in the m vertex [20]. We notice that the amlitudes listed in Eqs. () and (2) are indeed transverse. When introducing FF in these amlitudes, we cannot make these amlitudes kee transverse. Of course, finding a more suitable form of FF is an interesting research toic, where this FF not only reflects the realistic hysical icture but also can make the corresonding amlitudes still be transverse. In this work, we still choose the monoole FF in our calculation for reflecting the realistic hysical icture, and estimate the roduction of charmonium by the interaction rocesses. The general differential cross section of mψ is given by [3] dσ dt = 64π s cm 2 M 2, (5) where cm is the three-momentum of roton in the center-ofmass frame of. The overline indicates the average over the olarizations of / in the initial state and the sum over the olarization of m/ψ in the final state. With these obtained amlitudes listed in Eqs. ()-(4), we finally get the exressions of the differential cross sections of the ωψ rocesses without including the FF contribution, which are collected in aendix A. The corresonding total cross sections of ωψ are shown in aendix B. We also confirm the deductions of the differential and total cross sections of π 0 Ψ in Ref. [8]. Since the formulae of the differential and total cross sections of π 0 Ψ and ωψ with the FF contribution are very comlicated, we do not show their concrete exressions in detail, but directly aly their formulae to the numerical calculation. III. NUMERICAL RESULTS Before resenting the numerical result, we first introduce the couling constants adoted in our calculation. For the rocesses π 0 Ψ, the couling constant of the π interaction was given in many theoretical works, where we take g NNπ = 3.5 [2]. The couling constants of charmonium interacting with nucleons are not well established at resent. Thus, in this work we adoted the same values as those in Ref. [8], where these couling constants are g ηc = (9.0±3.2) 0 3, g J/ψ = (.62±0.03) 0 3, g ψ = (0.97±0.04) 0 3, g χc0 = (5.42±0.37) 0 3, g χc = (.03±0.07) 0 3, which were estimated by the measured artial width of Ψ [3]. For the discussed rocesses ωψ, there are two strong couling constants g ω andκ ω for the ω vertex. At resent, different models and different aroaches gave various theoretical values for these two couling constants, which are listed in Table. I. TABLE I: The estimated values of the couling constants g ω andκ ω. Mechanism/Model g ω κ ω Paris [2, 3] Nijmegen [4] Bonn [5] Pion hotoroduction [6] Nucleon EM form factors [7] 20.86± ±0.0 QCD sum rule [8] 8 ± ± P 0 quark model [] 3/2 Light meson emission model [0] 23± ± 2.0 3/2 In addition, the masses of the hadrons involved in our calculation include m π 0 = 35.0 MeV, m ω = MeV, m = MeV, m ηc = MeV, m J/ψ = MeV, m χc0 = MeV, m χc = MeV and m ψ = MeV, which are from Particle Data Grou [3]. A. The total cross sections and angular distributions of π 0 Ψ Firstly, in this work we confirm the results of the total cross sections of π 0 Ψ in Ref. [8], where the Ψ and π

4 4 vertices are as oint-like structures. Just shown in the lefthand diagram of Fig. 2, the deendence of the total cross sections of π 0 Ψ on the center-of-mass energy E cm is resented with a comarison with the exerimental data from the E760 exeriment [7], where the cross section of π 0 J/ψ near 3.5 GeV was measured at Fermilab. The theoretical result of π 0 J/ψ taking E cm = 3.5 GeV is larger than the exerimental value obviously. FIG. 3: (color online). The total cross sections of π 0 η c (lefthand) and π 0 J/ψ (right-hand) with differentλvalues. Here, the black oints with error bars are the E760 data [7] η c FF η c FIG. 2: (color online). The obtained total cross section of π 0 Ψ and the comarison of the exerimental and theoretical results. Here, the red oints with errors are the exerimental measurement from E760 [7]. The left-hand and right-hand diagrams corresond to the theoretical results without and with the FF contribution to π 0 Ψ, resectively. J/ψ FF J/ψ In Sec. II, we investigate the structure effect of the Ψ and ω vertices on the total cross sections of π 0 Ψ, which is resented in the right-hand diagram in Fig. 2. When takingλ =.9 GeV, we obtain the total cross section of π 0 J/ψ at E cm = 3.5 GeV consistent with the exerimental data, which indicates that the FF contribution cannot be ignored in studying π 0 Ψ. Adoting the sameλ value, we also obtain the total cross sections of other charmonium roductions, which are listed in the right-hand diagram of Fig. 2. By checking the results with and without considering the FF contribution, we find that there exist differences, where the obtained total cross sections of π 0 Ψ with FF are suressed comared with those without FF, which are shown in Fig. 2. In addition, with increasing E cm, the total cross section with FF goes down after reaching its maximum, while the total cross section without FF is goes u continuously. The forthcoming PANDA exeriment can test these theoretical results. FIG. 4: (color online). The center-of-mass frame angular distribution dσ/dω of π 0 η c and π 0 J/ψ. Here, the results are given by taking the range of E cm = GeV or GeV with ste of 0.2 GeV for theη c or J/ψ roduction. The diagrams in the first column are the results without FF [8] while the remaining diagrams are the results with FF, where we takeλ=.9 GeV. The results at the lower and uer limits of E cm range are highlighted with dot-dashed green and short-dashed red lines resectively, while long-dashed blue lines are the result at E cm = 4 GeV. All results are normalized to the forward intensity. With π 0 η c and π 0 J/ψ as examle, we also resent the variation of their total cross sections with different Λ values as shown in Fig. 3, where we take several tyical values in the range ofλ= with ste of 0.2 GeV. These results show that the total cross sections of π 0 η c and π 0 J/ψ deend on the value ofλ. Besides roviding the information of total cross section, we also give the result of the center-of-mass unolarized angular distributions of π 0 η c and π 0 J/ψ just suggested in Ref. [8]. Considering the FF contribution, the corresonding angular distributions dσ/dω are shown in the second column

5 5 χc0 χc ψ B. The obtained result for the ωψ reactions (a) (b) FF χc0 FF χc FF ψ FIG. 5: (color online). The center-of-mass frame angular distribution dσ/dω of π 0 χ c0, π 0 χ c and π 0 ψ corresonding to E cm = GeV, GeV, GeV with ste of 0.2 GeV, resectively. Here, the results without and with FF corresond to these diagrams in the first and the second rows, resectively. The style of the highlighted lines is arranged in the same way as in Fig. 4. (c) (d) of Fig. 4, where the results without FF are also listed. In addition, the center-of-mass unolarized angular distributions of π 0 χ c0, π 0 χ c and π 0 ψ are listed in Fig. 5. We can find that the angular distributions of π 0 Ψ with and without FF are anisotroic just shown in Figs. 4 and 5. Another common eculiarity of these results in Figs. 4 and 5 is that these obtained distributions become forwardand backward-eaked with increasing E cm, which is consistent with the conclusion in Ref. [8]. The comarison of the results in Figs. 4 and 5 also indicates the different behaviors of the center-of-mass frame angular distributions dσ/dω with and without FF, esecially for π 0 ψ, π 0 χ c0 and π 0 χ c. For π 0 η c, its dσ/dω distribution with FF become more forward- and backward-eaked than that without FF when taking the same E cm value. There obviously exists a node in its differential cross sections with and without FF when takingθ=90 in the center-of-mass frame or t=u, which is totally different from the corresonding result without FF, where θ is defined as the angle between the light meson and roton. In fact, this behavior is also suorted by analyzing the detailed analytic exressions of this rocess. For π 0 J/ψ, we notice that its differential cross section with FF is close to 0 when taking E cm = 5 GeV andθ=90 in the center-of-mass frame. The situations of π 0 ψ, π 0 χ c0 and π 0 χ c are similar to that of π 0 J/ψ. These results further indicate that the FF contribution should be considered in studying the charmonium roduction lus a light meson at the low energy interaction. The above investigation of the center-of-mass unolarized angular distributions of π 0 Ψ can rovide valuable information to the design of the PANDA detector and the analysis of the exerimental data. What is more imortant is that further exeriment can test the calculated total cross sections and the dσ/dω distributions of π 0 Ψ. (e) FIG. 6: The total cross sections of ωη c (left column) and ωj/ψ (right column) deendent on E cm. The diagrams in the first row are the result without FF, while the rest diagrams are the result with FF. Here, we take three different combinations of g ω andκ ω to illustrate the total cross sections deendent on (g ω,κ ω ) just shown in (a)-(d). Among these four diagrams, diagrams (c) and (d) are obtained by takingλ =.9 GeV. Diagrams (e) and (f) show the total cross sections with different Λ values and the tyical value (g ω,κ ω )=(2.2, -0.2). In the following, we illustrate the results of ωψ. Just shown in Table. I, there exist different values of g ω andκ ω for (f)

6 6 (gω,κω)=(2.2, -0.2) (gω,κω)=(2.5, 0.66) (gω,κω)=(5.9, 0) ηc ηc ηc FF ηc FF ηc FF ηc J/ψ J/ψ J/ψ FF J/ψ FF J/ψ FF J/ψ FIG. 7: The redicted total cross sections of ωψ corresonding to tyical valuesλ=.9 GeV and (g ω,κ ω )=(2.2, 0.2). the ω couling. With ωη c and ωj/ψ as examle, we list the variation of the total cross section of these two rocesses with E cm when taking three tyical combinations of g ω andκ ω. Our calculation indicates that the total cross sections of ωη c and ωj/ψ are deendent on the value of (g ω,κ ω ) whether we consider the FF contribution or not (see Fig. 6 (a)-(d) for more details). The results shown in Fig. 6 (a)-(d) indeed indicate that the ωψ rocesses can be alied to test the values of (g ω,κ ω ) listed in Table. I. It is obvious that the total cross sections of ωη c and ωj/ψ are suressed by FF. Since we cannot constrain theλvalue, in Fig. 6 (e)-(f) we discuss the deendence of the cross section of ωη c and ωj/ψ onλ. For other charmonium roductions withωmeson, the behavior of their cross sections is similar to that of ωη c and ωj/ψ. In Fig. 7, we further list the total cross sections for the ωψ rocesses taking tyical valuesλ=.9 GeV and (g ω,κ ω )=(2.2, 0.2), where we take the same cutoffλas that of π 0 J/ψ. Our calculation shows that the total cross section of ωη c is the largest one among all charmonium roductions discussed here. The roduction cross section of J/ψ is roughly 50 times smaller than that ofη c, while the roduction cross section ofχ c0 is 5 times larger than that of J/ψ. We also notice that the total cross sections of ωψ become stable after reaching u to their maximums with increasing E cm, which is different from the situation of π 0 Ψ discussed above. Just because of the considerable cross sections ofη c andχ c0 roductions, ωη c and ωχ c0 can be as the ideal channels to study the charmonium roduc- FIG. 8: (color online). The obtained center-of-mass frame angular distribution dσ/dω of ωη c and ωj/ψ. Here, the results are given by taking the range of E cm = GeV or GeV with ste of 0.2 GeV for theη c or J/ψ roduction, resectively. The diagrams in the first and third rows are the results without FF while the remaining diagrams are the results with FF, where we takeλ=.9 GeV. The diagrams with the same couling constants (g ω,κ ω ) are listed in the same column. The results at the lower and uer limits of E cm range are highlighted with dotted green and short-dashed red lines resectively, while long-dashed blue lines are the result at E cm = 5 GeV. tion at PANDA. Additionally, with ωη c and ωj/ψ as examle, we also give their center-of-mass frame angular distributions, which are shown in Fig. 8. There exists obvious difference among the results calculated with and without FF. We can do a cross check for our calculation, where we takeλ, which denotes that the ω and Ψ interactions are treated as the oint-like vertices. We find that the result with FF is gradually consistent with that without FF when Λ tends to infinity. We also find the node in the differential cross section with FF when takingθ = 90 in the center-of-mass frame, which is similar to π 0 Ψ. The investigation of ωψ further shows that the FF contribution cannot be ignored in the calculation.

7 7 C. Total cross sections for other rocesses We can easily extend the formulae of π 0 Ψ and ωψ to study ηψ and ρψ resectively, where we only need to relace the corresonding couling constants and masses. The couling constant of the η interaction is taken as g η =.5 []. In Ref. [22], the ρ couling constants (g ρ,κ ρ )=(3.249, 6.) obtained by the Bonn full model. In addition, the masses ofηandρmesons are M η = MeV and M ρ = GeV [3], resectively. In Fig. 9, the total cross sections of ηψ and ρψ are calculated by including the FF contribution, where we take the tyical cutoffλ=.9 GeV. The result also shows that the total cross section of theη c roduction is the largest one among all charmonium roductions by the ηψ or ρψ rocesses. FIG. 9: The total cross sections of the rocesses ηψ (lefthand) and ρψ (right-hand) deendent on E cm. Here, we take the tyical valueλ=.9 GeV. IV. DISCUSSIONS AND CONCLUSION As one of the main hysical aims of PANDA, studying the charmonium roduction at the low energy interaction is an interesting research toic not only for exerimentalists but also for theorists. For investigating the roduction of charmonium lus a ion by the annihilation, the roduction mechanism shown in Fig. was roosed in Ref. [5], where the π 0 J/ψ cross section was calculated. According to this mechanism, the authors in Refs. [6, 8 0] further develoed the model of the roduction of charmonium lus a light meson and alied it to calculate the π 0 Ψ rocesses, where the total cross sections and the unolarized angular distributions were obtained. Although there are many theoretical grous dedicated to the roduction of charmonium lus a ion by the interaction, further theoretical study of these rocesses is still valuable, where several further develoments of the model were indicated in Ref. [0]. Considering the resent research status of the roduction of charmonium lus a light meson by the low energy interaction, in this work we revisit this interesting research toic. Different from the former work in Ref. [8], in our calculation we consider the FF contribution to each interaction vertex, which was roosed as one of the ossible develoments of the model [8]. We calculate the total cross sections and the center-of-mass frame angular distributions of the π 0 Ψ rocesses. Since the E760 exeriment measured the cross section of π 0 J/ψ [7], we comare our result with the E760 data, which indicates that the calculated total cross section of π 0 Ψ overlas with the exerimental data, where the cutoffλ in FF is taken as.9 GeV. The inconsistence between the exerimental data and the result without FF in Ref. [8] are alleviated by considering the FF contribution, which indicates that the FF involved in each interaction vertex cannot be ignored. Adoting the sameλvalue, in this work we also give other charmonium roduction cross sections and the corresonding center-of-mass frame angular distributions of π 0 Ψ. The difference of the results with and without FF further shows that we should consider the FF contribution to our calculation. Thus, these results obtained with FF can be served as further exerimental investigation of charmonium roduction at PANDA. Besides studying the charmonium roduction lus a ion, in this work we also calculate the ωψ rocesses. The results of the total cross section and the center-of-mass frame angular distribution of ωψ also indicate the distinct difference of the calculations with and without FF, which is similar to the situation of π 0 Ψ. Different from the π couling, the ω interaction is related to two indeendent couling constants g ω andκ ω, which make ωψ an ideal channel to test different theoretical values of g ω andκ ω, where our results also show the total cross sections of ωψ are deendent on the values of g ω andκ ω. In order to reflect the comleteness of the study of the charmonium roduction lus a light meson, we also aly the formulae of π 0 Ψ and ωψ to calculate the ηψ and ρψ, resectively. The redicted total cross sections of ηψ and ρψ are accessible at the forthcoming PANDA. By the systematic investigation of the charmonium roduction lus a light meson at the low energy interaction, we notice a common behavior of the charamonium roduction, i.e., theη c roduction cross section is the largest one among all discussed rocesses. This fact shows that the low energy interaction is an ideal latform to roduceη c. In addition, we also find that the charmonium roduction cross sections satisfy the relationσ ηc >σ χc0 >σ J/ψ >σ χc >σ ψ, which does not deend on the associated light mesons. The authors in Ref. [6] calculated the total cross section of the J/ψ andψ roductions mj/ψ(ψ ) by relating these rocesses with J/ψ(ψ ) m. Thus, we also make a comarison of the results obtained by us and those listed in Ref. [6]. As shown in Table. II, our result with FF of ωψ,ρψ are consistent with those given in Ref. [6] while the cross sections of π 0 ψ,ηψ calculated by us are slightly smaller than those in Ref. [6]. The cross section of ρj/ψ shown in this work also falls into the range of the redicted cross section of ρj/ψ, where the u-

8 8 TABLE II: The comarison of the total cross sections of mj/ψ(ψ ) obtained in this work (the second column) and given in Ref. [6] (the third column) and the corresonding E cm value. Here, our results are the tyical values of the cross section with FF when takingλ=.9 GeV. Reaction σ FF (b) σ max CA (b) [6] E cm (GeV) π 0 J/ψ 6 420± ηj/ψ ± ωj/ψ ± ρj/ψ 27 < π 0 ψ 28 55±8 5.4 ηψ 9 33± ωψ 40 46± ρψ 32 38± er limit of the cross section of ρj/ψ was given in Ref. [6]. There also exist the differences of the results of π 0 J/ψ,ηJ/ψ,ωJ/ψ from this work and Ref. [6], i.e., total cross sections resented here are about 4 times, 42 times and 2 times smaller than those for π 0 J/ψ,ηJ/ψ,ωJ/ψ, resectively. In summary, in this work we systematically study the charmonium roduction lus a light meson in the interaction, and redict the corresonding total cross section and the center-of-mass frame unolarized angular distribution, which rovide valuable information to the exerimental investigation of this kind of reaction. As an ideal exeriment for studying the charmonium roduction, the forthcoming PANDA exeriment can directly verify the rediction given in this work and test the charmonium roduction mechanism at the low energy interaction adoted here. Acknowledgement This roject is suorted by the National Natural Science Foundation of China under Grant Nos , , , the Ministry of Education of China (FANEDD under Grant No , DPFIHE under Grant No , NCET, the Fundamental Research Funds for the Central Universities and the Fok Ying Tung Education Foundation (No. 3006). Aendix A: The differential cross section of ωψ By Eq. (5), we can easily deduce the exressions of the unolarized differential cross sections, i.e., dσ dt ωη c dσ dt ωχ c0 α ω α ηc (( = π R 2 (r 2 + 2) f 2 + ( f 2 + 2R 4 + 2r 2 R 2 + 2R 2 f )xy ) +κ ( 3r 2 R 2 f 2 (2 f 2 + 4r 2 R 2 )xy ) s(s 4m 2 )x2 y 2 8 κ2( (r 2 + 8)r 2 R 2 f 2 (r 2 f f 2 4r 2 R 4 4r 2 R 2 f+ 6r 2 R 2 )xy 4(R 2 + f )x 2 y 2) ), (A) α ω α χc0 (( = π (r 2 + 2)(R 2 4) f 2 + ( 2(r 2 + 2)(r 2 4)+8r 2 + 2R 4 2R (R 2 4) f+f 2) xy ) s(s 4m 2 )x2 y 2 +3κr 2 f [ (R 2 4) f+ 2xy ] 8 κ2[ r 2 (r 2 + 8)(R 2 4) f 2 ( 4r 4 R 2 + 4r 2 R 4 6r 2 R 2 + 4r 2 (R 2 8) f dσ dt ω(j/ψ ( ) ) +(r 2 4) f 2) xy 4(r 2 + R 2 + f )x 2 y 2]), (A2) α ω α J/ψ ( ) {[ = π 2 (r 2 + 2)(R 2 + 2) f 2 + ( 2(r 2 + 2)(R 2 + 2)+2(r 2 + R 2 + 2) f+ 2(r 2 + R 2 ) 2 8+ f 2) xy s(s 4m 2 )x2 y 2 2x 2 y 2] + 2κ [ 3r 2 (R 2 + 2) f 2 (6r 4 + 4r 2 R 2 + 6r 2 f+f 2 )xy ] + 4 κ2[ r 2 (r 2 + 8)(R 2 + 2) f 2 + (2r 6 + 6r 4 R 2 +6r 4 + 2r 4 f+ 32r 2 R 2 + 6r 2 f+ r 2 f 2 )xy+4(r 2 + f )x 2 y 2] }, (A3) dσ dt ωχ c α ω α χc { 2 [ = π R 2 s(s 4m 2 )x2 y 2 R 2 (r 2 + 2)(R 2 4) f 2 + ( 2R 6 + 6r 2 R 4 8R 4 6r 2 R 2 + 2r 4 R 2 + 2R 2 (r 2 + R 2 4) f +(R 2 + 2) f 2) xy 2R 2 x 2 y 2] 2κ R 2 [ 3r 2 R 2 (R 2 4) f 2 + ( 2r 4 R 2 + 6r 2 R 4 24r 2 R 2 2r 2 R 2 f+ (R 2 + 2) f 2) xy +2(r 2 3R 2 )x 2 y 2] κ2 4R 2 [ r 2 R 2 (r 2 + 8)(R 2 4) f 2 + ( 2r 6 R 2 + 2r 4 R 4 6r 4 R 2 24r 2 R r 2 R 2 + 2r 2 R 2 (r 2 +4) f+ (2r 2 r 2 R 2 8) f 2) xy 2 ( r 4 + r 2 R 2 + 4r 2 + 2R 4 2R 2 + (r 2 + 2R 2 4) f ) x 2 y 2 4x 3 y 3] }, (A4)

9 9 whereα ω g 2 ω /4π,α Ψ g 2 Ψ /4π, r m ω/m and R m Ψ /m. As the dimensionless variables, x and y are defined as x t/m 2 and y u/m 2, resectively. f denotes a dimensionless energy variable with definition f = (s m 2 ω m 2 Ψ )/m2 = (x+y). Aendix B: The total cross section of ωψ The detailed formulae of the total cross section of ωψ are σ ωηc = π α ωα ηc m 2 {[ R 2 (r 2 + 2) f 2 I s(s 4m ( f 2 + 2R 4 + 2r 2 R 2 + 2R 2 ] [ f )I +κ 3r 2 R 2 f 2 I 2 (2 f 2 + 4r 2 R 2 ] )I 8 κ2[ (r 2 + 8)r 2 R 2 f 2 I 2 (r 2 f f 2 4r 2 R 4 4r 2 R 2 f+ 6r 2 R 2 )I 4(R 2 + f )I 0 ]}, (B) σ ωχc0 σ ω(j/ψ ( ) ) σ ωχc = π α ωα χc0 m 2 s(s 4m 2 ) { [ (r 2 + 2)(R 2 4) f 2 I 2 + ( 2(r 2 + 2)(r 2 4)+8r 2 + 2R 4 2R (R 2 4) f + f 2) I ] + 3κr 2 f [ (R 2 4) fi 2 + 2I ] 8 κ2[ r 2 (r 2 + 8)(R 2 4) f 2 I 2 ( 4R 4 R 2 + 4r 2 R 4 6r 2 R 2 +4r 2 (R 2 8) f+ (r 2 4) f 2) I 4(r 2 + R 2 + f )I 0 ]}, (B2) = π α ωα (J/ψ ( ) )m 2 { [ 2 (r 2 + 2)(R 2 + 2) f 2 I s(s 4m ( 2(r 2 + 2)(R 2 + 2)+2(r 2 + R 2 + 2) f+ 2(r 2 + R 2 ) f 2) I 2I 0 ] + 2κ [ 3r 2 (R 2 + 2) f 2 I 2 (6r 4 + 4r 2 R 2 + 6r 2 f+f 2 )I ] + 4 κ2[ r 2 (r 2 + 8)(R 2 + 2) f 2 I 2 +(2r 6 + 6r 4 R 2 + 6r 4 + 2r 4 f+ 32r 2 R 2 + 6r 2 f+ r 2 f 2 )I + 4(R 2 + f )I 0 ]}, = π α ωα χc m 2 { 2 [ R 2 (r 2 + 2)(R 2 4) f 2 I s(s 4m 2 R ( 2R 6 + 6r 2 R 4 8R 4 6r 2 R 2 + 2r 4 R 2 + 2R 2 (r 2 + R 2 4) f +(R 2 + 2) f 2) I 2R 2 I 0 ] 2κ R 2 [ 3r 2 R 2 (R 2 4) f 2 I 2 + ( 2r 4 R 2 + 6r 2 R 4 24r 2 R 2 2r 2 R 2 f+ (R 2 + 2) f 2) I +2(r 2 3R 2 )I 0 ] κ2 4R 2 [r2 R 2 (r 2 + 8)(R 2 4) f 2 I 2 + ( 2r 6 R 2 + 2r 4 R 4 6r 4 R 2 24r 2 R r 2 R 2 + 2r 2 R 2 (r 2 +4) f+ (2r 2 r 2 R 2 8) f 2) I 2 ( r 4 + r 2 R 2 + 4r 2 + 2R 4 2R 2 + (r 2 + 2R 2 4) f ) I 0 4I ]}, (B4) (B3) where I m (m=, 0,, 2) is defined asi m = x x 0 dx (xy) m = I m (x ) I m (x 0 ) with y= x f, x 0 = m2 ω 2E E ω 2 ω m 2 and ( I (x)= 3 x3 + f ) 2 x2,, x = m2 ω 2E E ω + 2 ω, m 2 I 0 (x)= x, I (x)= ( x+ f f ln x I 2 (x)= 2 ( x+ f ) f ln f ( 3 x 2 ), x+ f + x ). When deducing these exressions of ωψ, the momenta and energies of roton andωmeson are related to the the Mandelstam variables by = ( ) s 4m 2 /2, 2 E = s/2 2, ω = (m 2 Ψ m2 ω )2 2(m 2 ω + m2 Ψ )s+ s2, 2 s ( ) E k = s m 2 2s /2 Ψ + m 2 ω. [] J. J. Aubert et al. [E598 Collaboration], Phys. Rev. Lett. 33, 404 (974). [2] J. E. Augustin et al. [SLAC-SP-07 Collaboration], Phys. Rev. Lett. 33, 406 (974).

10 0 [3] K. Nakamura et al. [Particle Data Grou Collaboration], J. Phys. G G 37, (200). [4] W. Erni et al. [PANDA Collaboration], arxiv: [heex]. [5] M. K. Gaillard, L. Maiani, R. Petronzio, Phys. Lett. B0, 489 (982). [6] A. Lundborg, T. Barnes, and U. Wiedner, Phys. Rev. D 73, (2006). [7] T. A. Armstrong, D. Bettoni, V. Bharadwaj, C. Biino, G. Borreani, D. R. Broemmelsiek, A. Buzzo, and R. Calabrese et al., Phys. Rev. Lett. 69, 2337 (992). [8] T. Barnes and X. Li, Phys. Rev. D 75, (2007). [9] T. Barnes, X. Li, and W. Roberts, Phys. Rev. D 77, (2008) [arxiv: [he-h]]. [0] T. Barnes, X. Li, and W. Roberts, Phys. Rev. D 8, (200). [] C. Downum, T. Barnes, J. R. Stone, and E. S. Swanson, Phys. Lett. B 638, 455 (2006). [2] W. N. Cottingham, M. Lacombe, B. Loiseau, J. M. Richard, and R. Vinh Mau, Phys. Rev. D 8, 800 (973). [3] M. Lacombe, B. Loiseau, J. M. Richard, R. Vinh Mau, J. Cote, P. Pires, and R. de Tourreil, Phys. Rev. C 2, 86 (980). [4] M. M. Nagels, T. A. Rijken, and J. J. de Swart, Phys. Rev. D 20, 633 (979). [5] R. Machleidt, Phys. Rev. C 63, (200). [6] T. Sato and T. S. H. Lee, Phys. Rev. C 54, 2660 (996). [7] P. Mergell, U. G. Meissner, and D. Drechsel, Nucl. Phys. A 596, 367 (996). [8] S. L. Zhu, Phys. Rev. C 59, 3455 (999). [9] B. Pire and L. Szymanowski, Phys. Lett. B 622, 83 (2005) [heh/ ]. [20] J. P. Lansberg, B. Pire, and L. Szymanowski, J. Phys. Conf. Ser. 295, (20) [arxiv: [he-h]]. [2] Z. W. Lin, C. M. Ko, and B. Zhang, Phys. Rev. C 6, (2000) [nucl-th/ ]. [22] R. Machleidt, K. Holinde, and C. Elster, Phys. Ret. 49, (987).