-
The radiative decay mode of the
f1(1285) resonance is interesting because it is the basic element in the description of thef1(1285) photoproduction data [1, 2]. It is also advocated as one of the observables most suitable for learning about the nature of thef1(1285) state [3-7]. Using the chiral unitary approach,f1(1285) appears as a pole in the complex plane of the scattering amplitude of theK∗ˉK+c.c. interaction in the isospinI=0 andJPC=1++ channel [8-10]. In other words, the axial-vector mesonf1(1285) can be taken as aK∗ˉK molecular state. For brevity, we useK∗ˉK to represent the positive C-parity combination ofK∗ˉK andˉK∗K in what follows.The experimental decay width of
f1(1285) is22.7±1.1 MeV [7], quite small compared with its mass. This is naturally explained in Ref. [8] using the molecular picture, implying thatf1(1285) is a dynamically generated state. TheK∗ˉK channel is the only allowed and considered pseudoscalar-vector channel in the chiral unitary approach, and the pole off1(1285) is below theK∗ˉK threshold; therefore, the total width of thef1(1285) resonance was not obtained in Ref. [8]. If the convolution of theK∗ width was taken into account, the partial decay width of theK∗ˉK channel would be approximately0.3 MeV (see more details in Ref. [8]). In fact, the dominant decay modes contributing to the width are peculiar. For example, theηππ channel accounts for 52% of the width, and the branching ratio ofπa0(980) channel is 38%. The decay off1(1285)→πa0(980) has been well investigated in Ref. [11] within theK∗ˉK molecular state picture forf1(1285) . These theoretical calculations in Ref. [11] have been confirmed in a recent BESIII experiment [12].There is another important decay channel, i.e., the
KˉKπ channel, the branching ratio of which is(9.1±0.4) % [7]. This decay mode was investigated in Ref. [13] with the same picture as in Ref. [11], and the theoretical predictions agree with existing experimental data. One could posit that the decay off1(1285)→ˉKK∗→KˉKπ should be much enhanced, owing to the strong coupling off1(1285) to theˉKK∗ channel. Actually, the mass off1(1285) is below the mass threshold ofˉKK∗ ; hence, it is easy to see that the above mechanism is much suppressed owing to the highly off-shell effect of theK∗ propagator, which was already found and discussed in Ref. [13] (see more details in that reference). Yet, all of the above tests have been performed for hadronic decay modes and not for radiative decays. In this work, we study the radiative decays of thef1(1285) resonance, assuming that it is aK∗ˉK state.On the experimental side, the particle data group (PDG) averaged values for the radiative decays of
f1(1285) are [7]①Br(f1→γρ0)=(5.3±1.2)%,
(1) Br(f1→γϕ)=(7.5±2.7)×10−4,
(2) which leads to the partial decay width
Γf1→γρ0=1.2±0.3 MeV and a ratioR1=Br(f1→γρ0)/Br(f1→γϕ)=71±30 . There is currently no experimental data on thef1(1285)→γω decay. On the other hand, the recent value ofΓf1→γρ0 obtained by the CLAS collaboration at Jafferson Lab, utilizing the analysis of theγp→pf1(1285) reaction, is much smaller, at0.45±0.18 MeV [1]. These values were obtained withBr(f1→ηππ)=0.52±0.02 [7]; the measured branching ratio wasBr(f1→γρ0)/ Br(f1→ηππ)=0.047±0.018 and the width wasΓf1= 18.4±1.4 MeV in Ref. [1]. The measured mass of thef1(1285) state wasMf1=1281.0±0.8 MeV, compatible with the known properties [7] of thef1(1285) resonance. On the theoretical side, the authors in Ref. [2] reportΓf1→γρ0=0.311 MeV andΓf1→γω=0.0343 MeV under the assumption thatf1(1285) has a quark-antiquark nature. ThisΓf1→γρ0 value is compatible with that obtained by the CLAS collaboration, within the error range, but is much smaller than the above PDG averaged value. Within the picture off1(1285) being a quark-antiquark state, another theoretical prediction for thef1(1285) radiative decay was reported in Ref. [14] using a covariant oscillator quark model. It predictedΓf1(1285)→γρ0 in the range of 0.509~0.565 MeV,Γf1(1285)→γω in the range of 0.048~0.057 MeV, andΓf1(1285)→γϕ in the range of 0.0056~0.02 MeV; these predictions depend on a particular mixing angle between the(uˉu+dˉd)/√2 andsˉs components. Note thatf1(1285) andf1(1420) are the members of the pseudovector nonet in theqˉq quark model [2, 14], wheref1(1285) is a mostlyuˉu+dˉd state andf1(1420) is ansˉs state. However, the study in Ref. [15] shows thatf1(1420) is not a genuine resonance and it shows up as a peak because of theK∗ˉK andπa0(980) decay modes off1(1285) around1420 MeV. In fact, as discussed by the PDG [7], although these two states are well known, their nature remains to be established. Thus, further investigations about them are needed [16].Here, we extend the work in Refs. [11, 13] for the hadronic decays of
f1(1285) to the case of radiative decays. In the molecular state scenario,f1(1285) decays intoγV (V=ρ0 ,ω , andϕ ) via kaon loop diagrams, and we can evaluate simultaneously these processes. It is shown that the theoretical results are in a good agreement with experimental data, hence supporting the strong coupling of thef1(1285) state to theˉKK∗ channel.The present paper is organized as follows. In Sec. 2, we discuss the formalism and the main ingredients of the model. In Sec. 3 we present our numerical results and conclusions. A short summary is given in the last section.
-
We study the
f1(1285)→γV decays under the assumption thatf1(1285) is dynamically generated from theK∗ˉK+c.c. interaction; thus, this decay can proceed viaf1(1285)→K∗ˉK→γV through triangle loop diagrams, which are shown in Fig. 1. In this mechanism,f1(1285) first decays intoK∗ˉK , thenK∗ decays intoKγ , andKˉK interacts to produce the vector meson V in the final state. We use p, k, and q for the momentum off1(1285) ,γ andK− andˉK0 in Figs. 1 (a, b) , respectively. Then, one can easily obtain that the momentum of the final vector meson isp−k , and the momenta ofK∗ and K arep−q andp−q−k , respectively. On the other hand, the decay off1(1285)→γV can also go withK∗ exchange, where one needs aK∗K∗γ vertex; then,K∗ˉK interacts to produce the vector meson V. However, it is easy to see that, compared with the mechanism shown in Fig. 1, this mechanism is strongly suppressed owing to the highly off-shell effect of the exchangedK∗ propagator when theK∗ˉK invariant mass is the mass of the vector meson V. In fact, as shown in Ref. [17], for the case ofa1/b1→γπ decays, the contribution of theK∗ exchange is rather small, on the order of 0.5%, compared with the one from the K exchange. Therefore, it is expected that the contributions from theK∗ exchange will be also small for thef1→γV decays, as studied here, and those contributions can be safely neglected. -
To evaluate the radiative decay of
f1(1285)→γV , we need the decay amplitudes of these diagrams, shown in Fig. 1. As mentioned above, thef1(1285) resonance is dynamically generated from the interaction ofK∗ˉK . For the charge conjugate transformation, we take the phase conventionsCK∗=−ˉK∗ andCK=ˉK , which are consistent with the standard chiral Lagrangians, and write|f1(1285)>=1√2(K∗ˉK−ˉK∗K)=−12(K∗+K−+K∗0ˉK0−K∗−K+−ˉK∗0K0) .
(3) Then we can write the
f1(1285)ˉKK∗ vertex as−itf1→ˉKK∗=−igf1C1ϵμ(f1)ϵμ(K∗),
(4) where
ϵμ(f1) andϵμ(K∗) stand for the polarization vector off1(1285) andK∗ (ˉK∗ ), respectively. We will take the value of the coupling constant ofgf1ˉKK∗(≡gf1=7555MeV) as obtained in the chiral unitary approach [8]. The factorsC1 account for the weight of eachˉKK∗ (KˉK∗ ) component off1(1285) , corresponding to thef1ˉKK∗ vertex for each diagram shown in Fig. 1, and can be easily obtained from Eq. (3) as,CA,B1=−12; CC,D1=12.
(5) For the
ˉKKV vertices, we take the effective Lagrangian describing the pseudoscalar-pseudoscalar-vector (PPV ) interaction as [18-21],LPPV=−ig<Vμ[P,∂μP]> ,
(6) where
g=M/2f=4.2 withM≈(mρ+mω)/2 andf=93 MeV. The symbol<> denotes the trace, while the pseudoscalar- and vector-nonets are collected in the P and V matrices, respectively. We can write them asVμ=(ω+ρ0√2ρ+K∗+ρ−ω−ρ0√2K∗0K∗−ˉK∗0ϕ)μ,
(7) and
P=(ξ1π+K+π−ξ2K0K−ˉK0ξ3),
(8) with
ξ1=1√2π0+1√3η+1√6η′ ,ξ2=−1√2π0+1√3η+1√6η′ , andξ3=−1√3η+2√6η′ .Thus, the
ˉKKV vertex can be written as−itˉKK→V=igC2(2q+k−p)μεμ(p−k,λV),
(9) where
εμ(p−k,λV) is the polarization vector of the vector meson. From Eq. (6) and from the explicit expressions for the V and P matrices as shown in Eqns. (7) and (8), the factorsC2 for each diagram shown in Fig. 1 can be obtained,CA,C2=−1√2;CB,D2=1√2;forρproduction,CA,C2=−1√2;CB,D2=−1√2;forωproduction,CA,C2=1;CB,D2=1;forϕproduction.
(10) In terms of Eqns. (5) and (10), it is easy to see that Figs. 1 (a, c) give the same contribution and Figs. 1 (b, d) also give the same contribution. We hence only consider Figs. 1 (a, b) in the following calculation.
In addition, according to the Lagrangian in Eq. (6), the
ϕ→KˉK decay width is given byΓϕ→KˉK=g2mϕ48π(1−4m2Km2ϕ)3/2,
and we can obtain the coupling
g≃4.5 with the averaged experimental value ofΓϕ→KˉK=1.77±0.02 MeV,mϕ=1019.46 MeV, andmK=(mK++mˉK0)/2=495.6 MeV as quoted by the PDG [7]. Hence, in this work, we will takeg=4.2 as in Eq. (6).For the electromagnetic vertex
K∗Kγ , the effective interaction Lagrangian takes the form as in Refs. [22-25]LK∗Kγ=egK∗KγmK∗εμναβ∂μK∗ν∂αAβK,
(11) where
K∗ν ,Aβ and K denote theK∗ vector meson, photon, and the K pseudoscalar meson, respectively. The partial decay width ofK∗→Kγ is given byΓK∗→Kγ=e2g2K∗Kγ96π(m2K∗−m2K)3m5K∗.
(12) The values of the coupling constants
gK∗Kγ can be determined from the experimental data [7],ΓK∗+→K+γ=50.3±4.6 keV andΓK∗0→K0γ=116.4±10.2 keV, which lead togK∗+K+γ=0.75±0.03,gK∗0K0γ=−1.14±0.05,
(13) where the small errors are determined with the uncertainties of
ΓK∗→Kγ as above. In addition, we fix the relative phase between the above two couplings, taking into account the quark model expectation [26]. -
The partial decay width of the
f1(1285)→γρ0 decay is given byΓf1(1285)→γρ0=Eγ12πM2f1∑λf1,λγ,λρ|MA+MB|2,
(14) where
MA andMB are the decay amplitudes in Figs. 1 (a, b), respectively, and the energy of the photon isEγ=|→k|=(M2f1−m2ρ0)/2Mf1 . In the cases ofω andϕ production, these can be obtained in a straightforward manner.The above amplitudes,
MA andMB , can be easily obtained with effective interactions. Here, we give explicitly the amplitudeMA for theρ0 production,MA=−eggf1gK∗+K+γ2√2mK∗+∫d4q(2π)41q2−m2K−+iϵ×12ω∗(q)D1Mf1−q0−ω∗(q)+iΓK∗+/2×D2(p−q−k)2−m2K++iϵ,
(15) where
ω∗(q)=√|→q|2+m2K∗+ is theK∗+ energy, and we have taken the positive energy part of theK∗ propagator into account, which is a good approximation, given the large mass ofK∗ (see more details in Ref. [11]). In Eq. (15), the factorsD1 andD2 read②D1=εμναβ(p−q)μεν(p,λf1)kαε∗β(k,λγ),
(16) D2=(2q+k−p)σε∗σ(p−k,λρ) ,
(17) with
λf1 ,λγ , andλρ the spin polarizations off1(1285) , photon, andρ0 meson, respectively. The amplitudeMB corresponding to Fig. 1 (b) can be easily obtained through the substitutionsmK∗+→mK∗0 ,mK+→mK0 , andmK−→mˉK0 intoMA . The decay amplitudes off1(1285)→γϕ andf1(1285)→γω share the similar formalism as in Eq. (15).To calculate
MA in Eq. (15), we first integrate overq0 using Cauchy's theorem. For doing this, we take the rest frame off1(1285) , in which one can writep=(Mf1,0,0,0),k=(Eγ,0,0,Eγ),
(18) q=(q0,|→q|sinθcosϕ,|→q|sinθsinϕ,|→q|cosθ),
(19) with
θ andϕ as the polar and azimuthal angles of→q along the→k direction. The energy of the final vector meson isEV=(M2f1+m2V)/2Mf1 . Then, we haveV1=D1D2=∓iEγ|→q|2sin2θ,
(20) for
λf1=0 ,λγ=±1 , andλρ=∓1 , andV2=D1D2=±i2E2γmρ0(q0−Mf1−|→q|cosθ)×(q0+EVEγ|→q|cosθ),
(21) for
λf1=±1 ,λγ=±1 , andλρ=0 . Notice that we have dropped those terms containingsinϕ orcosϕ , because after the integration over the azimuthal angleϕ , they do not yield contributions.After integrating over
q0 in Eq. (15), we haveFA1=|→q|4(1−cos2θ)ωω′ω∗(XA1+XA2+XA3),
(22) FA2=|→q|2ωω′ω∗[(Mf1−ω∗−EVEγ|→q|cosθ)(ω∗+|→q|cosθ)XA1+(ω−Mf1−|→q|cosθ)(ω+EVEγ|→q|cosθ)XA2+(ω′−Eγ−|→q|cosθ)(EV+ω′+EVEγ|→q|cosθ)XA3],
(23) where
XA1=1(Mf1−ω∗−ω+iΓK∗+2)(Eγ−ω∗−ω′+iΓK∗+2),XA2=1(Mf1−ω∗−ω+iΓK∗+2)(EV−ω−ω′+iϵ),XA3=1(ω+ω∗−Eγ−iΓK∗+2)(EV+ω+ω′−iϵ),
with
ω′=√|→q|2+E2γ+2Eγ|→q|cosθ+m2K+ andω=√|→q|2+m2K− the energies ofK− andK+ in the diagram of Fig. 1 (a).FB1 andFB2 will be obtained just by applying the substitution toFA1 andFA2 withmK∗+→mK∗0 ,mK−→mˉK0 , andmK+→mK0 .Finally, the partial decay width takes the form
Γf1→γV=e2g2g2f1E5γ192π2M2f1m2V∑i=1,2|∫Λ0d|→q|∫1−1dcosθ×(CAFAi+CBFBi)|2,
(24) with
CA=−√24gK∗+K+γmK∗+,forV=ρ0,ω,
(25) CA=12gK∗+K+γmK∗+,forV=ϕ,
(26) CB=√24gK∗0K0γmK∗0,forV=ρ0,
(27) CB=−√24gK∗0K0γmK∗0,forV=ω,
(28) CB=−12gK∗0K0γmK∗0,forV=ϕ.
(29) For
ρ0 production, the relative minus sign betweenCA andCB combined with the minus sign between the couplingsgK∗+K+γ andgK∗0K0γ is positive, and hence the interference of the two diagrams(a) and(b) shown in Fig. 1 is constructive. However, it is destructive forω andϕ production, which makeΓf1(1285)→γρ0 much larger compared with the other two partial decay widths.In Eq. (24), we have introduced a momentum cutoff
Λ for preventing the ultraviolet divergence and for compensating the off-shell effects that appear in the triangle loop integral. It can also be done by introducing form factors to the intermediate particles, as shown in Refs. [27-32].Again, we want to stress that, in this work, those contributions of the
K∗ exchange via diagrams containing anomalous vector-vector-pseudoscalar (VVP) vertices are not taken into account.③ Such contributions were extensively studied in Refs. [17, 33-35] for the low-lying scalar, axial vector, and tensor meson radiative decays. As discussed in Refs. [33, 34], these contributions are very sensible to the exact value of the VVP coupling. Furthermore, including such diagrams, the decay amplitudes would become more complex, owing to additional model parameters, which cannot be exactly determined. Hence, we leave these contributions to further studies when more precise experimental measurements become available. -
In this section, we explain how the large
ρ0 width contributions are implemented. We studyf1(1285)→γρ0 with theρ0→π+π− decay. For this purpose we replaceΓf1→γρ0 in Eq. (24) by¯Γf1→γρ0 :¯Γf1→γρ0=∫(mρ0+2Γ0ρ0)2(mρ0−2Γ0ρ0)2d˜m2S(˜m)Γf1→γρ0(mρ0→˜m),
(30) where
˜m is the invariant mass of theπ+π− system. Then,S(˜m) has the formS(˜m)=−1πIm(1˜m2−m2ρ0+imρ0Γρ(˜m)),
(31) where
Γρ(˜m) is energy-dependent, and it can be written as [36-42],ρ(˜m)=Γ0ρ0(˜m2−4m2πm2ρ0−4m2π)3/2,
(32) with
mρ0=775.26 MeV,Γ0ρ0=149.1 MeV andmπ=mπ+=mπ−=139.57 MeV. -
The partial decay width of the
f1(1285)→γV decay as a function ofΛ from 800 to 1500 MeV is illustrated in Fig. 2, where the black solid, dashed, and dotted curves stand for the theoretical results of theρ0 ,ω , andϕ production. It is worth mentioning that the results forω andϕ are multiplied by a factor of 100, while the red solid line stands for the results for theρ0 production but with the contributions of theρ0 mass as in Eq. (30). One can see that, from Fig. 2, the theoretical results have the same order of magnitude within the given range of the cutoff parameterΛ values. In the considered range of cutoffs,Γf1→γρ0 varies from0.4 to0.9 MeV, which is consistent with the experimental result within the error range [1, 7]. In addition, the contribution of theρ0 width is also important and it will reduce the numerical results ofΓf1→γρ0 by a factor of 18%.Figure 2. (color online) Partial decay width of the
f1(1285)→γV decay as a function of the cutoff parameterΛ . The black solid, dashed, and dotted curves denote the results for theρ0 ,ω , andϕ production, while the results forω andϕ are multiplied by a factor of 100. The red solid line denotes the results for theρ0 production but with the contributions of theρ0 mass as in Eq. (30).In Table 1 we show explicitly the numerical results for the
f1(1285)→γV decays with some particular cutoff parameters. We show also the theoretical calculations of Refs. [2, 14] and the experimental results [1, 7], for comparison.Λ f1→γρ0 Γ (¯Γ )f1→γω [×10−2 ]f1→γϕ [×10−2 ]R1 R2 800 0.42 (0.34 )1.36 0.71 59 31 1000 0.56 (0.46 )1.87 0.93 60 30 1500 0.88 (0.72 )3.01 1.41 62 29 Ref. [2] 0.311 3.43 — — 9 Ref. [14] (set I)④ 0.509 4.8 2.0 25 11 Ref. [14] (set II) 0.565 5.7 0.56 101 10 Exp. [7] 1.2±0.3 — 1.7±0.6 71±30 — Exp. [1]a 0.45±0.18 — — — — aThe measured width of f1(1285) is ~6 MeV smaller than the previous world average [7].Table 1. Partial decay width for
f1(1285)→γV . All units are in MeV.In general, we cannot provide the value of the cutoff parameter; however, if we divide
Γf1(1285)→γρ0 byΓf1(1285)→γω orΓf1(1285)→γϕ , the dependence of these ratios on the cutoff will be smoothed. Two ratios are defined asR1=Γf1(1285)→γρ0Γf1(1285)→γϕ,R2=Γf1(1285)→γρ0Γf1(1285)→γω.
(33) These two ratios are correlated with each other. With
R1 measured experimentally, one can fix the cutoff in the model and predict the ratioR2 . We also show, in Table 1 , the explicit numerical results forR1 andR2 , for some particular cutoff parameters.In Fig. 3, we show the numerical results for the above two ratios, where the solid line denotes the results for
R1 , while the dashed line denotes the results forR2 . Indeed, one can see that the dependence of both ratios on the cutoffΛ is rather weak. The ratioR1≃60 is in agreement with the experimental result71±30 [7]. On the other hand, the result forR2 is approximately30 . We can conclude firmly that the partial decay width off1(1285)→γρ0 is much larger than the ones toγω andγϕ channels. This is owing to the destructive interference between Figs. 1 (a, b) forω andϕ production. Our present conclusion agrees wtih quark model calculations [2, 14]. However, from Table 1 one can see that the presently obtained ratiosR1 andR2 are much different from the values obtained by the quark models, especially forR2 . In the quark model calculations,R2 is always around9 , which is owing to the isospin difference ofρ0 andω mesons. We hope that future experimental measurements will help to clarify this issue.Figure 3. (color online) The
Λ dependence of the ratiosR1 (solid line) andR2 (dashed line) defined in Eq. (33). The error band corresponds to the experimental result forR1. It is worth mentioning that there is only one free parameter
Λ in the present work (all the other parameters were fixed in previous works). In addition, the dependence ofR1 andR2 on the cutoffΛ is rather weak; thus, these can be the model predictions, and they would be compared with future experimental measurements.In addition, we want to note that, although we have assumed that
f1(1285) is a dynamically generated state, the numerical results here are not tied to the assumed nature off1(1285) . The crucial point is that it couples strongly to theˉKK∗ channel, whatever its origin. -
We have evaluated the partial decay rates of the radiative decays
f1(1285)→γV with the assumption thatf1(1285) is a dynamically generated state from the strongˉK∗K interaction, and in this picture thef1(1285) state has a strong coupling to theˉKK∗ channel. The theoretical results we obtained for the partial widths are sensitive to the free parameterΛ , but they are compatible with experimental data within the error range. Furthermore, the ratiosR1=Γf1→γρ0Γf1→γϕ andR2=Γf1→γρ0Γf1→γω , which are not sensitive to the only free parameterΛ , are predicted. It is found that the values ofR1 andR2 obtained here are different from other theoretical predictions using quark models. The precise experimental observations of those radiative decays would then provide very valuable information about the relevance of the strong coupling off1(1285) to theˉKK∗ channel.
