Radiative decays of f1(1285) as the KˉK molecular state

  • With f1(1285) as a dynamically generated resonance from KˉK interactions, we estimate the rates of the radiative transitions of the f1(1285) meson to the vector mesons ρ0, ω and ϕ. These radiative decays proceed via the kaon loop diagrams. The calculated results are in a fair agreement with experimental measurements. Some predictions can be tested experimentally; their analysis will be valuable for decoding the strong coupling of the f1(1285) state to the ˉKK channel.
  • The radiative decay mode of the f1(1285) resonance is interesting because it is the basic element in the description of the f1(1285) photoproduction data [1, 2]. It is also advocated as one of the observables most suitable for learning about the nature of the f1(1285) state [3-7]. Using the chiral unitary approach, f1(1285) appears as a pole in the complex plane of the scattering amplitude of the KˉK+c.c. interaction in the isospin I=0 and JPC=1++ channel [8-10]. In other words, the axial-vector meson f1(1285) can be taken as a KˉK molecular state. For brevity, we use KˉK to represent the positive C-parity combination of KˉK and ˉKK in what follows.

    The experimental decay width of f1(1285) is 22.7±1.1 MeV [7], quite small compared with its mass. This is naturally explained in Ref. [8] using the molecular picture, implying that f1(1285) is a dynamically generated state. The KˉK channel is the only allowed and considered pseudoscalar-vector channel in the chiral unitary approach, and the pole of f1(1285) is below the KˉK threshold; therefore, the total width of the f1(1285) resonance was not obtained in Ref. [8]. If the convolution of the K width was taken into account, the partial decay width of the KˉK channel would be approximately 0.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 of f1(1285)πa0(980) has been well investigated in Ref. [11] within the KˉK molecular state picture for f1(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 of f1(1285)ˉKKKˉKπ should be much enhanced, owing to the strong coupling of f1(1285) to the ˉKK channel. Actually, the mass of f1(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 the K 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 the f1(1285) resonance, assuming that it is a Kˉ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)×104,

    (2)

    which leads to the partial decay width Γf1γρ0=1.2±0.3 MeV and a ratio R1=Br(f1γρ0)/Br(f1γϕ)=71±30. There is currently no experimental data on the f1(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 γppf1(1285) reaction, is much smaller, at 0.45±0.18 MeV [1]. These values were obtained with Br(f1ηππ)=0.52±0.02 [7]; the measured branching ratio was Br(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 the f1(1285) state was Mf1=1281.0±0.8 MeV, compatible with the known properties [7] of the f1(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 that f1(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 of f1(1285) being a quark-antiquark state, another theoretical prediction for the f1(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 and sˉs components. Note that f1(1285) and f1(1420) are the members of the pseudovector nonet in the qˉq quark model [2, 14], where f1(1285) is a mostly uˉu+dˉd state and f1(1420) is an sˉs state. However, the study in Ref. [15] shows that f1(1420) is not a genuine resonance and it shows up as a peak because of the KˉK and πa0(980) decay modes of f1(1285) around 1420 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 the f1(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 that f1(1285) is dynamically generated from the KˉK+c.c. interaction; thus, this decay can proceed via f1(1285)KˉKγV through triangle loop diagrams, which are shown in Fig. 1. In this mechanism, f1(1285) first decays into KˉK, then K decays into Kγ, and KˉK interacts to produce the vector meson V in the final state. We use p, k, and q for the momentum of f1(1285), γ and K and ˉK0 in Figs. 1 (a, b) , respectively. Then, one can easily obtain that the momentum of the final vector meson is pk, and the momenta of K and K are pq and pqk, respectively. On the other hand, the decay of f1(1285)γV can also go with K exchange, where one needs a KKγ 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 exchanged K propagator when the KˉK invariant mass is the mass of the vector meson V. In fact, as shown in Ref. [17], for the case of a1/b1γπ decays, the contribution of the K 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 the K exchange will be also small for the f1γV decays, as studied here, and those contributions can be safely neglected.

    Figure 1

    Figure 1.  Triangle loop diagrams representing the process f1(1285)γV , with V being the ρ0, ω, or ϕ meson.

    To evaluate the radiative decay of f1(1285)γV, we need the decay amplitudes of these diagrams, shown in Fig. 1. As mentioned above, the f1(1285) resonance is dynamically generated from the interaction of KˉK. For the charge conjugate transformation, we take the phase conventions CK=ˉK and CK=ˉK, which are consistent with the standard chiral Lagrangians, and write

    |f1(1285)>=12(KˉKˉKK)=12(K+K+K0ˉK0KK+ˉK0K0) .

    (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 of f1(1285) and K (ˉK), respectively. We will take the value of the coupling constant of gf1ˉKK(gf1=7555MeV) as obtained in the chiral unitary approach [8]. The factors C1 account for the weight of each ˉKK (KˉK) component of f1(1285), corresponding to the f1ˉ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 with M(mρ+mω)/2 and f=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 as

    Vμ=(ω+ρ02ρ+K+ρωρ02K0KˉK0ϕ)μ,

    (7)

    and

    P=(ξ1π+K+πξ2K0KˉK0ξ3),

    (8)

    with ξ1=12π0+13η+16η, ξ2=12π0+13η+16η, and ξ3=13η+26η.

    Thus, the ˉKKV vertex can be written as

    itˉKKV=igC2(2q+kp)μεμ(pk,λV),

    (9)

    where εμ(pk,λ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 factors C2 for each diagram shown in Fig. 1 can be obtained,

    CA,C2=12;CB,D2=12;forρproduction,CA,C2=12;CB,D2=12;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π(14m2Km2ϕ)3/2,

    and we can obtain the coupling g4.5 with the averaged experimental value of ΓϕKˉK=1.77±0.02 MeV, mϕ=1019.46 MeV, and mK=(mK++mˉK0)/2=495.6 MeV as quoted by the PDG [7]. Hence, in this work, we will take g=4.2 as in Eq. (6).

    For the electromagnetic vertex KKγ, the effective interaction Lagrangian takes the form as in Refs. [22-25]

    LKKγ=egKKγmKεμναβμKναAβK,

    (11)

    where Kν, Aβ and K denote the K vector meson, photon, and the K pseudoscalar meson, respectively. The partial decay width of KKγ is given by

    ΓKKγ=e2g2KKγ96π(m2Km2K)3m5K.

    (12)

    The values of the coupling constants gKKγ can be determined from the experimental data [7], ΓK+K+γ=50.3±4.6 keV and ΓK0K0γ=116.4±10.2 keV, which lead to

    gK+K+γ=0.75±0.03,gK0K0γ=1.14±0.05,

    (13)

    where the small errors are determined with the uncertainties of ΓKKγ 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 and MB are the decay amplitudes in Figs. 1 (a, b), respectively, and the energy of the photon is Eγ=|k|=(M2f1m2ρ0)/2Mf1. In the cases of ω and ϕ production, these can be obtained in a straightforward manner.

    The above amplitudes, MA and MB, can be easily obtained with effective interactions. Here, we give explicitly the amplitude MA for the ρ0 production,

    MA=eggf1gK+K+γ22mK+d4q(2π)41q2m2K+iϵ×12ω(q)D1Mf1q0ω(q)+iΓK+/2×D2(pqk)2m2K++iϵ,

    (15)

    where ω(q)=|q|2+m2K+ is the K+ energy, and we have taken the positive energy part of the K propagator into account, which is a good approximation, given the large mass of K (see more details in Ref. [11]). In Eq. (15), the factors D1 and D2 read

    D1=εμναβ(pq)μεν(p,λf1)kαεβ(k,λγ),

    (16)

    D2=(2q+kp)σεσ(pk,λρ) ,

    (17)

    with λf1, λγ, and λρ the spin polarizations of f1(1285), photon, and ρ0 meson, respectively. The amplitude MB corresponding to Fig. 1 (b) can be easily obtained through the substitutions mK+mK0, mK+mK0, and mKmˉK0 into MA. The decay amplitudes of f1(1285)γϕ and f1(1285)γω share the similar formalism as in Eq. (15).

    To calculate MA in Eq. (15), we first integrate over q0 using Cauchy's theorem. For doing this, we take the rest frame of f1(1285), in which one can write

    p=(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 is EV=(M2f1+m2V)/2Mf1. Then, we have

    V1=D1D2=iEγ|q|2sin2θ,

    (20)

    for λf1=0, λγ=±1, and λρ=1, and

    V2=D1D2=±i2E2γmρ0(q0Mf1|q|cosθ)×(q0+EVEγ|q|cosθ),

    (21)

    for λf1=±1, λγ=±1, and λρ=0. Notice that we have dropped those terms containing sinϕ or cosϕ, because after the integration over the azimuthal angle ϕ, they do not yield contributions.

    After integrating over q0 in Eq. (15), we have

    FA1=|q|4(1cos2θ)ωωω(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 of K and K+ in the diagram of Fig. 1 (a). FB1 and FB2 will be obtained just by applying the substitution to FA1 and FA2 with mK+mK0, mKmˉK0, and mK+mK0.

    Finally, the partial decay width takes the form

    Γf1γV=e2g2g2f1E5γ192π2M2f1m2Vi=1,2|Λ0d|q|11dcosθ×(CAFAi+CBFBi)|2,

    (24)

    with

    CA=24gK+K+γmK+,forV=ρ0,ω,

    (25)

    CA=12gK+K+γmK+,forV=ϕ,

    (26)

    CB=24gK0K0γmK0,forV=ρ0,

    (27)

    CB=24gK0K0γmK0,forV=ω,

    (28)

    CB=12gK0K0γmK0,forV=ϕ.

    (29)

    For ρ0 production, the relative minus sign between CA and CB combined with the minus sign between the couplings gK+K+γ and gK0K0γ 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 study f1(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ρ02Γ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 form

    S(˜m)=1πIm(1˜m2m2ρ0+imρ0Γρ(˜m)),

    (31)

    where Γρ(˜m) is energy-dependent, and it can be written as [36-42],

    ρ(˜m)=Γ0ρ0(˜m24m2πm2ρ04m2π)3/2,

    (32)

    with mρ0=775.26 MeV, Γ0ρ0=149.1 MeV and mπ=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 from 0.4 to 0.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

    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.

    Table 1

    Table 1.  Partial decay width for f1(1285)γV. All units are in MeV.
    Λ f1γρ0 Γ (¯Γ) f1γω [×102] f1γϕ [×102] 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].
    DownLoad: CSV
    Show Table

    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 as

    R1=Γ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 ratio R2. We also show, in Table 1 , the explicit numerical results for R1 and R2 , 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 for R2. Indeed, one can see that the dependence of both ratios on the cutoff Λ is rather weak. The ratio R160 is in agreement with the experimental result 71±30 [7]. On the other hand, the result for R2 is approximately 30. We can conclude firmly that the partial decay width of f1(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 ratios R1 and R2 are much different from the values obtained by the quark models, especially for R2. In the quark model calculations, R2 is always around 9, 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

    Figure 3.  (color online) The Λ dependence of the ratios R1 (solid line) and R2 (dashed line) defined in Eq. (33). The error band corresponds to the experimental result for R1.

    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 of R1 and R2 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 of f1(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 that f1(1285) is a dynamically generated state from the strong ˉKK interaction, and in this picture the f1(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 ratios R1=Γf1γρ0Γf1γϕ and R2=Γf1γρ0Γf1γω, which are not sensitive to the only free parameter Λ, are predicted. It is found that the values of R1 and R2 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 of f1(1285) to the ˉKK channel.

    [1] R. Dickson et al. [CLAS Collaboration], Phys. Rev. C, 93: 065202 (2016) doi: 10.1103/PhysRevC.93.065202
    [2] A. A. Osipov, A. A. Pivovarov, and M. K. Volkov, Phys. Rev. D, 96: 054012 (2017) doi: 10.1103/PhysRevD.96.054012
    [3] E. Oset et al., Int. J. Mod. Phys. E, 25: 1630001 (2016)
    [4] L. S. Geng and E. Oset, Phys. Rev. D, 94: 014018 (2016) doi: 10.1103/PhysRevD.94.014018
    [5] J. M. Dias, V. R. Debastiani, J. J. Xie et al., Chin. Phys. C, 42: 043106 (2018) doi: 10.1088/1674-1137/42/4/043106
    [6] Y. Kalashnikova, A. E. Kudryavtsev, A. V. Nefediev et al., Phys. Rev. C, 73: 045203 (2006)
    [7] M. Tanabashi et al., Phys. Rev. D, 98: 030001 (2018)
    [8] L. Roca, E. Oset, and J. Singh, Phys. Rev. D, 72: 014002 (2005) doi: 10.1103/PhysRevD.72.014002
    [9] Y. Zhou, X. L. Ren, H. X. Chen et al., Phys. Rev. D, 90: 014020 (2014) doi: 10.1103/PhysRevD.90.014020
    [10] L. S. Geng, X. L. Ren, Y. Zhou et al., Phys. Rev. D, 92: 014029 (2015) doi: 10.1103/PhysRevD.92.014029
    [11] F. Aceti, J. M. Dias and E. Oset, Eur. Phys. J. A, 51: 48 (2015) doi: 10.1140/epja/i2015-15048-5
    [12] M. Ablikim et al. [BESIII Collaboration], Phys. Rev. D, 92: 012007 (2015) doi: 10.1103/PhysRevD.92.012007
    [13] F. Aceti, J. J. Xie, and E. Oset, Phys. Lett. B, 750: 609 (2015) doi: 10.1016/j.physletb.2015.09.068
    [14] S. Ishida, K. Yamada, and M. Oda, Phys. Rev. D, 40: 1497 (1989)
    [15] V. R. Debastiani, F. Aceti, W. H. Liang et al., Phys. Rev. D, 95: 034015 (2017) doi: 10.1103/PhysRevD.95.034015
    [16] W. H. Liang and E. Oset, Eur. Phys. J. C, 80: 407 (2020) doi: 10.1140/epjc/s10052-020-7966-y
    [17] L. Roca, A. Hosaka, and E. Oset, Phys. Lett. B, 658: 17 (2007) doi: 10.1016/j.physletb.2007.10.035
    [18] M. Bando, T. Kugo, S. Uehara et al., Phys. Rev. Lett., 54: 1215 (1985) doi: 10.1103/PhysRevLett.54.1215
    [19] M. Bando, T. Kugo, and K. Yamawaki, Phys. Rept., 164: 217 (1988) doi: 10.1016/0370-1573(88)90019-1
    [20] M. Harada and K. Yamawaki, Phys. Rept., 381: 1 (2003) doi: 10.1016/S0370-1573(03)00139-X
    [21] U. G. Meissner, Phys. Rept., 161: 213 (1988) doi: 10.1016/0370-1573(88)90090-7
    [22] Y. Oh and H. Kim, Phys. Rev. C, 73: 065202 (2006) doi: 10.1103/PhysRevC.73.065202
    [23] S. H. Kim, S. i. Nam, Y. Oh et al., Phys. Rev. D, 84: 114023 (2011) doi: 10.1103/PhysRevD.84.114023
    [24] S. H. Kim, S. i. Nam, A. Hosaka et al., Phys. Rev. D, 88: 054012 (2013) doi: 10.1103/PhysRevD.88.054012
    [25] A. C. Wang, W. L. Wang, F. Huang et al., Phys. Rev. C, 96: 035206 (2017)
    [26] F. E. Close, Academic Press/london 1979, 481p
    [27] X. H. Liu and Q. Zhao, Phys. Rev. D, 81: 014017 (2010) doi: 10.1103/PhysRevD.81.014017
    [28] G. Li, X. h. Liu, Q. Wang et al., Phys. Rev. D, 88: 014010 (2013) doi: 10.1103/PhysRevD.88.014010
    [29] W. Qin, C. An, G. Li, C. Wang et al., Eur. Phys. J. C, 79: 757 (2019) doi: 10.1140/epjc/s10052-019-7274-6
    [30] Y. Zhang and G. Li, Phys. Rev. D, 97: 014018 (2018) doi: 10.1103/PhysRevD.97.014018
    [31] X. H. Liu, G. Li, J. J. Xie et al., Phys. Rev. D, 100: 054006 (2019) doi: 10.1103/PhysRevD.100.054006
    [32] Q. Wu, G. Li, F. Shao et al., Phys. Rev. D, 94: 014015 (2016) doi: 10.1103/PhysRevD.94.014015
    [33] H. Nagahiro, L. Roca, and E. Oset, Eur. Phys. J. A, 36: 73 (2008) doi: 10.1140/epja/i2008-10567-8
    [34] H. Nagahiro, L. Roca, A. Hosaka et al., Phys. Rev. D, 79: 014015 (2009) doi: 10.1103/PhysRevD.79.014015
    [35] R. Molina, H. Nagahiro, A. Hosaka et al., Phys. Rev. D, 83: 094030 (2011) doi: 10.1103/PhysRevD.83.094030
    [36] H. C. Chiang, E. Oset, and L. C. Liu, Phys. Rev. C, 44: 738 (1991) doi: 10.1103/PhysRevC.44.738
    [37] L. S. Geng and E. Oset, Phys. Rev. D, 79: 074009 (2009) doi: 10.1103/PhysRevD.79.074009
    [38] C. Hanhart, Y. S. Kalashnikova, and A. V. Nefediev, Phys. Rev. D, 81: 094028 (2010) doi: 10.1103/PhysRevD.81.094028
    [39] J. J. Xie and E. Oset, Eur. Phys. J. A, 51: 111 (2015) doi: 10.1140/epja/i2015-15111-3
    [40] J. J. Xie and G. Li, Eur. Phys. J. C, 78: 861 (2018) doi: 10.1140/epjc/s10052-018-6353-4
    [41] X. Zhang and J. J. Xie, Commun. Theor. Phys., 70: 060 (2018) doi: 10.1088/0253-6102/70/1/60
    [42] X. Zhang and J. J. Xie, Chin. Phys. C, 43: 064104 (2019) doi: 10.1088/1674-1137/43/6/064104
  • [1] R. Dickson et al. [CLAS Collaboration], Phys. Rev. C, 93: 065202 (2016) doi: 10.1103/PhysRevC.93.065202
    [2] A. A. Osipov, A. A. Pivovarov, and M. K. Volkov, Phys. Rev. D, 96: 054012 (2017) doi: 10.1103/PhysRevD.96.054012
    [3] E. Oset et al., Int. J. Mod. Phys. E, 25: 1630001 (2016)
    [4] L. S. Geng and E. Oset, Phys. Rev. D, 94: 014018 (2016) doi: 10.1103/PhysRevD.94.014018
    [5] J. M. Dias, V. R. Debastiani, J. J. Xie et al., Chin. Phys. C, 42: 043106 (2018) doi: 10.1088/1674-1137/42/4/043106
    [6] Y. Kalashnikova, A. E. Kudryavtsev, A. V. Nefediev et al., Phys. Rev. C, 73: 045203 (2006)
    [7] M. Tanabashi et al., Phys. Rev. D, 98: 030001 (2018)
    [8] L. Roca, E. Oset, and J. Singh, Phys. Rev. D, 72: 014002 (2005) doi: 10.1103/PhysRevD.72.014002
    [9] Y. Zhou, X. L. Ren, H. X. Chen et al., Phys. Rev. D, 90: 014020 (2014) doi: 10.1103/PhysRevD.90.014020
    [10] L. S. Geng, X. L. Ren, Y. Zhou et al., Phys. Rev. D, 92: 014029 (2015) doi: 10.1103/PhysRevD.92.014029
    [11] F. Aceti, J. M. Dias and E. Oset, Eur. Phys. J. A, 51: 48 (2015) doi: 10.1140/epja/i2015-15048-5
    [12] M. Ablikim et al. [BESIII Collaboration], Phys. Rev. D, 92: 012007 (2015) doi: 10.1103/PhysRevD.92.012007
    [13] F. Aceti, J. J. Xie, and E. Oset, Phys. Lett. B, 750: 609 (2015) doi: 10.1016/j.physletb.2015.09.068
    [14] S. Ishida, K. Yamada, and M. Oda, Phys. Rev. D, 40: 1497 (1989)
    [15] V. R. Debastiani, F. Aceti, W. H. Liang et al., Phys. Rev. D, 95: 034015 (2017) doi: 10.1103/PhysRevD.95.034015
    [16] W. H. Liang and E. Oset, Eur. Phys. J. C, 80: 407 (2020) doi: 10.1140/epjc/s10052-020-7966-y
    [17] L. Roca, A. Hosaka, and E. Oset, Phys. Lett. B, 658: 17 (2007) doi: 10.1016/j.physletb.2007.10.035
    [18] M. Bando, T. Kugo, S. Uehara et al., Phys. Rev. Lett., 54: 1215 (1985) doi: 10.1103/PhysRevLett.54.1215
    [19] M. Bando, T. Kugo, and K. Yamawaki, Phys. Rept., 164: 217 (1988) doi: 10.1016/0370-1573(88)90019-1
    [20] M. Harada and K. Yamawaki, Phys. Rept., 381: 1 (2003) doi: 10.1016/S0370-1573(03)00139-X
    [21] U. G. Meissner, Phys. Rept., 161: 213 (1988) doi: 10.1016/0370-1573(88)90090-7
    [22] Y. Oh and H. Kim, Phys. Rev. C, 73: 065202 (2006) doi: 10.1103/PhysRevC.73.065202
    [23] S. H. Kim, S. i. Nam, Y. Oh et al., Phys. Rev. D, 84: 114023 (2011) doi: 10.1103/PhysRevD.84.114023
    [24] S. H. Kim, S. i. Nam, A. Hosaka et al., Phys. Rev. D, 88: 054012 (2013) doi: 10.1103/PhysRevD.88.054012
    [25] A. C. Wang, W. L. Wang, F. Huang et al., Phys. Rev. C, 96: 035206 (2017)
    [26] F. E. Close, Academic Press/london 1979, 481p
    [27] X. H. Liu and Q. Zhao, Phys. Rev. D, 81: 014017 (2010) doi: 10.1103/PhysRevD.81.014017
    [28] G. Li, X. h. Liu, Q. Wang et al., Phys. Rev. D, 88: 014010 (2013) doi: 10.1103/PhysRevD.88.014010
    [29] W. Qin, C. An, G. Li, C. Wang et al., Eur. Phys. J. C, 79: 757 (2019) doi: 10.1140/epjc/s10052-019-7274-6
    [30] Y. Zhang and G. Li, Phys. Rev. D, 97: 014018 (2018) doi: 10.1103/PhysRevD.97.014018
    [31] X. H. Liu, G. Li, J. J. Xie et al., Phys. Rev. D, 100: 054006 (2019) doi: 10.1103/PhysRevD.100.054006
    [32] Q. Wu, G. Li, F. Shao et al., Phys. Rev. D, 94: 014015 (2016) doi: 10.1103/PhysRevD.94.014015
    [33] H. Nagahiro, L. Roca, and E. Oset, Eur. Phys. J. A, 36: 73 (2008) doi: 10.1140/epja/i2008-10567-8
    [34] H. Nagahiro, L. Roca, A. Hosaka et al., Phys. Rev. D, 79: 014015 (2009) doi: 10.1103/PhysRevD.79.014015
    [35] R. Molina, H. Nagahiro, A. Hosaka et al., Phys. Rev. D, 83: 094030 (2011) doi: 10.1103/PhysRevD.83.094030
    [36] H. C. Chiang, E. Oset, and L. C. Liu, Phys. Rev. C, 44: 738 (1991) doi: 10.1103/PhysRevC.44.738
    [37] L. S. Geng and E. Oset, Phys. Rev. D, 79: 074009 (2009) doi: 10.1103/PhysRevD.79.074009
    [38] C. Hanhart, Y. S. Kalashnikova, and A. V. Nefediev, Phys. Rev. D, 81: 094028 (2010) doi: 10.1103/PhysRevD.81.094028
    [39] J. J. Xie and E. Oset, Eur. Phys. J. A, 51: 111 (2015) doi: 10.1140/epja/i2015-15111-3
    [40] J. J. Xie and G. Li, Eur. Phys. J. C, 78: 861 (2018) doi: 10.1140/epjc/s10052-018-6353-4
    [41] X. Zhang and J. J. Xie, Commun. Theor. Phys., 70: 060 (2018) doi: 10.1088/0253-6102/70/1/60
    [42] X. Zhang and J. J. Xie, Chin. Phys. C, 43: 064104 (2019) doi: 10.1088/1674-1137/43/6/064104
  • 加载中

Cited by

1. Lin, J.-X., Li, J.-T., Liang, W.-H. et al. J/ψ decays into ω(ϕ) f1(1285) and ω(ϕ) ` ` f1(1420) ”[J]. European Physical Journal C, 2024, 84(1): 52. doi: 10.1140/epjc/s10052-024-12405-5
2. He, D., Xie, Y., Sun, H. Role of triangle singularity in the decay process D0 →π+π-f0 (980), f0 →π+π-[J]. Physical Review D, 2023, 107(9): 094040. doi: 10.1103/PhysRevD.107.094040
3. Wang, K., Liu, B.-C. The π-p→ a0- ηp reaction in an effective Lagrangian model[J]. Physical Review C, 2023, 107(2): 025203. doi: 10.1103/PhysRevC.107.025203
4. Shen, Q.-H., Xie, J.-J. Faddeev fixed-center approximation to the ηk∗ K ¯ ∗, πk∗ K ¯ ∗, and KK∗ K ¯ ∗ systems[J]. Physical Review D, 2023, 107(3): 034019. doi: 10.1103/PhysRevD.107.034019
5. Yang, F., Zhu, H.Q., Huang, Y. Analysis of the η1(1855) as a KK¯1(1400) molecular state[J]. Nuclear Physics A, 2023. doi: 10.1016/j.nuclphysa.2022.122571
6. Molina, R., Döring, M., Liang, W.H. et al. The πf(500) decay of the a1(1260)[J]. European Physical Journal C, 2021, 81(9): 782. doi: 10.1140/epjc/s10052-021-09574-y
7. Lebiedowicz, P., Nachtmann, O., Salabura, P. et al. Exclusive f1 (1285) meson production for energy ranges available at the GSI-FAIR with HADES and PANDA[J]. Physical Review D, 2021, 104(3): 034031. doi: 10.1103/PhysRevD.104.034031
8. Ling, X.-Z., Liu, M.-Z., Lu, J.-X. et al. Can the nature of a0 (980) be tested in the Ds+ →π+π0η decay?[J]. Physical Review D, 2021, 103(11): 116016. doi: 10.1103/PhysRevD.103.116016

Figures(3) / Tables(1)

Get Citation
Ju-Jun Xie, Gang Li and Xiao-Hai Liu. Radiative decays of f1(1285) as the KˉK molecular state[J]. Chinese Physics C. doi: 10.1088/1674-1137/abae51
Ju-Jun Xie, Gang Li and Xiao-Hai Liu. Radiative decays of f1(1285) as the KˉK molecular state[J]. Chinese Physics C.  doi: 10.1088/1674-1137/abae51 shu
Milestone
Received: 2020-02-20
Revised: 2020-07-08
Article Metric

Article Views(2025)
PDF Downloads(36)
Cited by(8)
Policy on re-use
To reuse of Open Access content published by CPC, for content published under the terms of the Creative Commons Attribution 3.0 license (“CC CY”), the users don’t need to request permission to copy, distribute and display the final published version of the article and to create derivative works, subject to appropriate attribution.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Radiative decays of f1(1285) as the KˉK molecular state

  • 1. Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China
  • 2. School of Physics and Engineering, Qufu Normal University, Shandong 273165, China
  • 3. Center for Joint Quantum Studies and Department of Physics, School of Science, Tianjin University, Tianjin 300350, China
  • 4. School of Nuclear Science and Technology, University of Chinese Academy of Sciences, Beijing 101408, China
  • 5. School of Physics and Microelectronics, Zhengzhou University, Zhengzhou, Henan 450001, China

Abstract: With f1(1285) as a dynamically generated resonance from KˉK interactions, we estimate the rates of the radiative transitions of the f1(1285) meson to the vector mesons ρ0, ω and ϕ. These radiative decays proceed via the kaon loop diagrams. The calculated results are in a fair agreement with experimental measurements. Some predictions can be tested experimentally; their analysis will be valuable for decoding the strong coupling of the f1(1285) state to the ˉKK channel.

    HTML

    1.   Introduction
    • The radiative decay mode of the f1(1285) resonance is interesting because it is the basic element in the description of the f1(1285) photoproduction data [1, 2]. It is also advocated as one of the observables most suitable for learning about the nature of the f1(1285) state [3-7]. Using the chiral unitary approach, f1(1285) appears as a pole in the complex plane of the scattering amplitude of the KˉK+c.c. interaction in the isospin I=0 and JPC=1++ channel [8-10]. In other words, the axial-vector meson f1(1285) can be taken as a KˉK molecular state. For brevity, we use KˉK to represent the positive C-parity combination of KˉK and ˉKK in what follows.

      The experimental decay width of f1(1285) is 22.7±1.1 MeV [7], quite small compared with its mass. This is naturally explained in Ref. [8] using the molecular picture, implying that f1(1285) is a dynamically generated state. The KˉK channel is the only allowed and considered pseudoscalar-vector channel in the chiral unitary approach, and the pole of f1(1285) is below the KˉK threshold; therefore, the total width of the f1(1285) resonance was not obtained in Ref. [8]. If the convolution of the K width was taken into account, the partial decay width of the KˉK channel would be approximately 0.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 of f1(1285)πa0(980) has been well investigated in Ref. [11] within the KˉK molecular state picture for f1(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 of f1(1285)ˉKKKˉKπ should be much enhanced, owing to the strong coupling of f1(1285) to the ˉKK channel. Actually, the mass of f1(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 the K 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 the f1(1285) resonance, assuming that it is a Kˉ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)×104,

      (2)

      which leads to the partial decay width Γf1γρ0=1.2±0.3 MeV and a ratio R1=Br(f1γρ0)/Br(f1γϕ)=71±30. There is currently no experimental data on the f1(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 γppf1(1285) reaction, is much smaller, at 0.45±0.18 MeV [1]. These values were obtained with Br(f1ηππ)=0.52±0.02 [7]; the measured branching ratio was Br(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 the f1(1285) state was Mf1=1281.0±0.8 MeV, compatible with the known properties [7] of the f1(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 that f1(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 of f1(1285) being a quark-antiquark state, another theoretical prediction for the f1(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 and sˉs components. Note that f1(1285) and f1(1420) are the members of the pseudovector nonet in the qˉq quark model [2, 14], where f1(1285) is a mostly uˉu+dˉd state and f1(1420) is an sˉs state. However, the study in Ref. [15] shows that f1(1420) is not a genuine resonance and it shows up as a peak because of the KˉK and πa0(980) decay modes of f1(1285) around 1420 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 the f1(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.

    2.   Formalism
    • We study the f1(1285)γV decays under the assumption that f1(1285) is dynamically generated from the KˉK+c.c. interaction; thus, this decay can proceed via f1(1285)KˉKγV through triangle loop diagrams, which are shown in Fig. 1. In this mechanism, f1(1285) first decays into KˉK, then K decays into Kγ, and KˉK interacts to produce the vector meson V in the final state. We use p, k, and q for the momentum of f1(1285), γ and K and ˉK0 in Figs. 1 (a, b) , respectively. Then, one can easily obtain that the momentum of the final vector meson is pk, and the momenta of K and K are pq and pqk, respectively. On the other hand, the decay of f1(1285)γV can also go with K exchange, where one needs a KKγ 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 exchanged K propagator when the KˉK invariant mass is the mass of the vector meson V. In fact, as shown in Ref. [17], for the case of a1/b1γπ decays, the contribution of the K 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 the K exchange will be also small for the f1γV decays, as studied here, and those contributions can be safely neglected.

      Figure 1.  Triangle loop diagrams representing the process f1(1285)γV , with V being the ρ0, ω, or ϕ meson.

    • 2.1.   Effective interactions and coupling constants

    • To evaluate the radiative decay of f1(1285)γV, we need the decay amplitudes of these diagrams, shown in Fig. 1. As mentioned above, the f1(1285) resonance is dynamically generated from the interaction of KˉK. For the charge conjugate transformation, we take the phase conventions CK=ˉK and CK=ˉK, which are consistent with the standard chiral Lagrangians, and write

      |f1(1285)>=12(KˉKˉKK)=12(K+K+K0ˉK0KK+ˉK0K0) .

      (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 of f1(1285) and K (ˉK), respectively. We will take the value of the coupling constant of gf1ˉKK(gf1=7555MeV) as obtained in the chiral unitary approach [8]. The factors C1 account for the weight of each ˉKK (KˉK) component of f1(1285), corresponding to the f1ˉ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 with M(mρ+mω)/2 and f=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 as

      Vμ=(ω+ρ02ρ+K+ρωρ02K0KˉK0ϕ)μ,

      (7)

      and

      P=(ξ1π+K+πξ2K0KˉK0ξ3),

      (8)

      with ξ1=12π0+13η+16η, ξ2=12π0+13η+16η, and ξ3=13η+26η.

      Thus, the ˉKKV vertex can be written as

      itˉKKV=igC2(2q+kp)μεμ(pk,λV),

      (9)

      where εμ(pk,λ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 factors C2 for each diagram shown in Fig. 1 can be obtained,

      CA,C2=12;CB,D2=12;forρproduction,CA,C2=12;CB,D2=12;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π(14m2Km2ϕ)3/2,

      and we can obtain the coupling g4.5 with the averaged experimental value of ΓϕKˉK=1.77±0.02 MeV, mϕ=1019.46 MeV, and mK=(mK++mˉK0)/2=495.6 MeV as quoted by the PDG [7]. Hence, in this work, we will take g=4.2 as in Eq. (6).

      For the electromagnetic vertex KKγ, the effective interaction Lagrangian takes the form as in Refs. [22-25]

      LKKγ=egKKγmKεμναβμKναAβK,

      (11)

      where Kν, Aβ and K denote the K vector meson, photon, and the K pseudoscalar meson, respectively. The partial decay width of KKγ is given by

      ΓKKγ=e2g2KKγ96π(m2Km2K)3m5K.

      (12)

      The values of the coupling constants gKKγ can be determined from the experimental data [7], ΓK+K+γ=50.3±4.6 keV and ΓK0K0γ=116.4±10.2 keV, which lead to

      gK+K+γ=0.75±0.03,gK0K0γ=1.14±0.05,

      (13)

      where the small errors are determined with the uncertainties of ΓKKγ as above. In addition, we fix the relative phase between the above two couplings, taking into account the quark model expectation [26].

    • 2.2.   Decay amplitudes

    • 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 and MB are the decay amplitudes in Figs. 1 (a, b), respectively, and the energy of the photon is Eγ=|k|=(M2f1m2ρ0)/2Mf1. In the cases of ω and ϕ production, these can be obtained in a straightforward manner.

      The above amplitudes, MA and MB, can be easily obtained with effective interactions. Here, we give explicitly the amplitude MA for the ρ0 production,

      MA=eggf1gK+K+γ22mK+d4q(2π)41q2m2K+iϵ×12ω(q)D1Mf1q0ω(q)+iΓK+/2×D2(pqk)2m2K++iϵ,

      (15)

      where ω(q)=|q|2+m2K+ is the K+ energy, and we have taken the positive energy part of the K propagator into account, which is a good approximation, given the large mass of K (see more details in Ref. [11]). In Eq. (15), the factors D1 and D2 read

      D1=εμναβ(pq)μεν(p,λf1)kαεβ(k,λγ),

      (16)

      D2=(2q+kp)σεσ(pk,λρ) ,

      (17)

      with λf1, λγ, and λρ the spin polarizations of f1(1285), photon, and ρ0 meson, respectively. The amplitude MB corresponding to Fig. 1 (b) can be easily obtained through the substitutions mK+mK0, mK+mK0, and mKmˉK0 into MA. The decay amplitudes of f1(1285)γϕ and f1(1285)γω share the similar formalism as in Eq. (15).

      To calculate MA in Eq. (15), we first integrate over q0 using Cauchy's theorem. For doing this, we take the rest frame of f1(1285), in which one can write

      p=(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 is EV=(M2f1+m2V)/2Mf1. Then, we have

      V1=D1D2=iEγ|q|2sin2θ,

      (20)

      for λf1=0, λγ=±1, and λρ=1, and

      V2=D1D2=±i2E2γmρ0(q0Mf1|q|cosθ)×(q0+EVEγ|q|cosθ),

      (21)

      for λf1=±1, λγ=±1, and λρ=0. Notice that we have dropped those terms containing sinϕ or cosϕ, because after the integration over the azimuthal angle ϕ, they do not yield contributions.

      After integrating over q0 in Eq. (15), we have

      FA1=|q|4(1cos2θ)ωωω(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 of K and K+ in the diagram of Fig. 1 (a). FB1 and FB2 will be obtained just by applying the substitution to FA1 and FA2 with mK+mK0, mKmˉK0, and mK+mK0.

      Finally, the partial decay width takes the form

      Γf1γV=e2g2g2f1E5γ192π2M2f1m2Vi=1,2|Λ0d|q|11dcosθ×(CAFAi+CBFBi)|2,

      (24)

      with

      CA=24gK+K+γmK+,forV=ρ0,ω,

      (25)

      CA=12gK+K+γmK+,forV=ϕ,

      (26)

      CB=24gK0K0γmK0,forV=ρ0,

      (27)

      CB=24gK0K0γmK0,forV=ω,

      (28)

      CB=12gK0K0γmK0,forV=ϕ.

      (29)

      For ρ0 production, the relative minus sign between CA and CB combined with the minus sign between the couplings gK+K+γ and gK0K0γ 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.

    • 2.3.   The ρ0 width contributions

    • In this section, we explain how the large ρ0 width contributions are implemented. We study f1(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ρ02Γ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 form

      S(˜m)=1πIm(1˜m2m2ρ0+imρ0Γρ(˜m)),

      (31)

      where Γρ(˜m) is energy-dependent, and it can be written as [36-42],

      ρ(˜m)=Γ0ρ0(˜m24m2πm2ρ04m2π)3/2,

      (32)

      with mρ0=775.26 MeV, Γ0ρ0=149.1 MeV and mπ=mπ+=mπ=139.57 MeV.

    3.   Numerical results and discussion
    • 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 from 0.4 to 0.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γω [×102] f1γϕ [×102] 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 as

      R1=Γ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 ratio R2. We also show, in Table 1 , the explicit numerical results for R1 and R2 , 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 for R2. Indeed, one can see that the dependence of both ratios on the cutoff Λ is rather weak. The ratio R160 is in agreement with the experimental result 71±30 [7]. On the other hand, the result for R2 is approximately 30. We can conclude firmly that the partial decay width of f1(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 ratios R1 and R2 are much different from the values obtained by the quark models, especially for R2. In the quark model calculations, R2 is always around 9, 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 ratios R1 (solid line) and R2 (dashed line) defined in Eq. (33). The error band corresponds to the experimental result for R1.

      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 of R1 and R2 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 of f1(1285). The crucial point is that it couples strongly to the ˉKK channel, whatever its origin.

    4.   Summary
    • We have evaluated the partial decay rates of the radiative decays f1(1285)γV with the assumption that f1(1285) is a dynamically generated state from the strong ˉKK interaction, and in this picture the f1(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 ratios R1=Γf1γρ0Γf1γϕ and R2=Γf1γρ0Γf1γω, which are not sensitive to the only free parameter Λ, are predicted. It is found that the values of R1 and R2 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 of f1(1285) to the ˉKK channel.

Reference (42)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return