A study of the gamma* - f0(980) transition form factors

The spin wave function of the f0(980)-meson is constructed under the assumption that the meson is dominantly a quark-antiquark state. This wave function is used in a calculation of the gamma* - f0 transition form factors. In the real-photon limit the results for the transverse form factor are compared to the large momentum transfer data measured by the BELLE collaboration recently. It turns out that, for the momentum-transfer range explored by BELLE, the collinear approximation does not suffer, power corrections to it, modeled as quark transverse moment effects, seem to be needed. Mixing of the f0 with the sigma(500) is also briefly discussed.


Introduction
Recently the BELLE collaboration [1] has measured the cross section for γ * γ → π 0 π 0 for large photon virtuality, Q 2 1 , and small energy in the γ * γ center-of-mass system. From these data the photon-meson transition form factors have been extracted for the scalar, f 0 (980), and tensor, f 2 (1270), mesons for Q 2 1 < ∼ 30 GeV 2 . These transition form factors are similar to those for the pseudoscalar mesons which have been extensively studied by both experimentalists and theoreticians. In Ref. [2] the γ − f 0 and the γ − f 2 form factors have been investigated within the NRQCD factorization framework [3], in which relativistic corrections and higher Fock state contributions are suppressed by powers of the relativistic velocity of the quarks in the meson, i.e. up to some minor modifications, the light mesons are treated like heavy Quarkonia. Super-convergence relations have been derived in [4] and shown to provide constraints on the γ − f 2 transition form factor. The latter form factor has also been studied within the framework of collinear factorization [5]. A phenomenological model for this form factor is discussed in [6]. The process γ * γ → ππ has been discussed in the framework of generalized distribution amplitudes, time-like versions of generalized parton distributions [7]. In this paper the interest is focused on the γ − f 0 transition form factor.
The f 0 (980) meson is a complicated system whose nature is not yet fully understood. Its peculiar properties have led to many speculations about its quark content. A comparison of the partial widths for the f 0 decays into pairs of pions and Kaons [8] under regard of the respective phase spaces reveals that the matrix element for f 0 → K + K − is much larger than that for f 0 → π + π − . Thus, if the f 0 is viewed as a quark-antiquark state, it is dominantly an ss state. The comparison of the branching ratios for the radiative decays of the φ-meson into the f 0 and π 0 leads to the same conclusion. However, the f 0 -meson is not a pure ss state as is, for instance, obvious from the decay widths for J/Ψ → f 0 ω and J/Ψ → f 0 φ. This fact is interpreted as f 0 −σ(500) mixing. Detailed phenomenological analyses of f 0 −σ mixing in various decay processes [9,10,11,12] revealed two ranges for the mixing angle, ϕ, A light scalar glueball may affect this result [9].
As an alternative to the quark-antiquark interpretation other authors [13,14] have suggested a tetraquark configuration for the f 0 -meson. This appears as a natural explanation for the fact that the a 0 (980) and the f 0 mesons are degenerate in mass and are the heaviest particles of the lightest scalar-meson nonet. For the tetraquark interpretation there seems to be no f 0 − σ mixing [11]. The drawback of this picture is that the two-pion decay of the f 0 is too small as compared to experiment whereas the a 0 → ηπ is too large. In [15] it has been suggested that the lightest scalar-meson nonet, considered as tetraquarks states, mixes with the scalar-meson nonet with masses around 1200 MeV under the effect of the instanton force. The latter nonet is believed to have a predominant qq structure. This mixing leads to a better description of the light scalar-meson decays. The f 0 may also have a substantial KK molecule component [16]. It goes without saying that the real f 0 -meson is a superposition of all these configurations.
The goal of the present paper is the calculation of the γ * − f 0 transition form factors at large photon virtualities. For this calculation the pQCD framework developed by Brodsky and Lepage [17] is utilized in which the process is factorized in a perturbatively calculable hard subprocess (here γ * γ * → qq) and a soft hadronic matrix element, parametrized as a light-cone wave function, which is under control of soft, long-distance QCD. As any hadron the f 0 -meson possesses a Fock decomposition [18] starting with the simple quark-antiquark components where Ψ 2,β is the light-cone wave function of the qq Fock state; the index β labels its decomposition in flavor, color and helicity. The integration measures are defined by In the photon-photon interactions at large photon virtualities the f 0 -meson is generated through its lowest Fock components, mainly the ss one. As can be shown [17] the hard generation of the f 0 through higher Fock components is suppressed by inverse powers of the photon virtuality and is therefore neglected. Once the meson is produced it gets dressed by fluctuations into higher Fock components under the effect of long-distance QCD. The calculation of the γ * − f 0 transition form factors is similar to the one of the photon-pseudoscalar-meson form factors [17]. The latter calculation is to be generalized in such a way that also hadrons with non-zero orbital angular between their constituents can be treated.
The paper is organized as follows: In the next section the spin part of the light-cone wave function, termed the spin wave function, of the f 0 is constructed assuming that this mesons is an ss state. In Sect. 2.1 the collinear reduction of the spin wave function is discussed and, in Sect. 2.2, an example of a light-cone wave function of the f 0 is introduced and compared to the twist-2 and 3 distribution amplitudes. The γ * − f 0 transition form factors are defined in Sect. 3.1, followed by a LO perturbative calculation within the modified perturbative approach in which quark transverse degrees of freedom are retained (Sect 3.2). Numerical results for the form factors in the realphoton limit are given in Sect. 4.1 and compared to the BELLE data. Some comments on the behavior of the γ * − f 0 form factors are presented in Sect. 4.2. Finally, the summary will be given in Sect. 5.

The spin wave function of the f 0 -meson
For the description of the hadron the light-cone approach is used which enables one to completely separate the dynamical and kinematical features of the Poincaré invariance [19,20]. The overall motion of the hadron is decoupled from the internal motion of the constituents, i.e. the light-cone wave function of the hadron, Ψ, is independent of the hadron's momentum and is invariant under the kinematical Poincaré transformations ( boosts along and rotations around the 3-directions as well as transverse boosts). Hence, Ψ is determined if it is known at rest. The ss Fock component given in (2), is split in a spin part (hereafter denoted as spin wave function) and a reduced light-cone wave function, Ψ 0 , which represents the full, soft wave function, Ψ, with a factor K µ removed from it. As discussed in detail in Ref. [21] the covariant spin wave function can be constructed starting from the observation [22] that, in zero binding energy approximation, an equal-time hadron state (in the spin basis) in the constituent center-of-mass frame equals the (helicity) light-cone state at rest. Consequently, one can use the standard ls coupling scheme in order to couple quark and antiquark to a state of given spin and parity. On boosting the results to a frame with arbitrary hadron momentum one easily reads off the covariant spin wave function 2 .
Since the f 0 (980)-meson is a J P C = 0 ++ state the quark and antiquark have to couple in a spin-1 state and one unit of orbital angular momenta is required 3 . The ls coupling scheme leads to the following ansatz for the spin wave function of a final state meson in its rest frame [21,23] (4) Note that µ 1 , µ 2 denote spin components and v,ū are equal-t spinors here. In the meson's rest frame the meson and the constituent momenta read where k is the three-momentum part of the relative momentum of quark and antiquark In order to retain a covariant formulation, the four-vector K = (0, k) is introduced 4 . As is customary in the parton model, the binding energy is neglected and the constituents are considered as quasi on-shell particles. That possibly crude approximation can be achieved by putting the minus components of the constituents to zero. Hence, k 3 = 0 and our relative vector reduces to In this case the spin wave function (4) reads where k ⊥± = k ⊥1 ± ik ⊥2 . This spin wave function is of the same type as is discussed in [24] for the l = 1 Fock components of ρ and π-mesons. In the infinite momentum frame (IMF), obtained by boosting the meson rest frame momenta along the 3-direction, p · K = 0 holds and the quark and 3 In spectroscopy notation the valence Fock component of f 0 -meson is a 3 P 0 state. 4 As discussed in [21] each unit of orbital angular momentum will be represented by is the velocity 4-vector. In the rest frame clearly K ⊥ → (0, k) and one has the appropriate object transforming as a 3-vector under O(3). Thus, K µ introduced in the line after (6), is strictly speaking K µ ⊥ .
antiquark momenta are parametrized as 5 whereτ = 1 − τ and The boost to the IMF leads to: For convenience the variable ξ = 1−2τ is introduced. For ξ = 0 this covariant spin wave function coincides with the one employed for the χ c0 in [26]. The normalization of the spin wave function is chosen such that where E is the meson's energy. The meson's ss Fock state (2) explicitly reads The number of colors is denoted by N c and c,c are color labels. Proper state normalization requires the condition 1 2 where P f 0 is the probability of the ss Fock component. 5 In [25] the parton momenta are parametrized as wherep is a light-like vector whose 3-component points in the opposite direction of p. For this parametrization momentum conservation only holds up to corrections of order k 2 ⊥ /p. It however also leads to the spin wave function (11) up to corrections of order k 3 ⊥ /M 0 .

Collinear reduction
In collinear approximation the limit k ⊥ → 0 in the hard subprocess is to be taken in general. However, terms ∝ K α in it combine with terms linear in K in the spin wave function and therefore survive the k ⊥ -integration of the wave function. These terms are in general of the same order as the other terms in the spin wave function and it is therefore unjustified to neglect these terms 6 . Consider the expansion of the subprocess amplitude with respect to K: where A 0 is of order 1 while A 1 is of order 1/p + for dimensional reason. The k ⊥ -integration yields Formally this is equivalent to the replacement The transverse metric tensor is defined by g 11 ⊥ = g 22 ⊥ = −1 while all other components are zero in a frame where the meson moves along the 3-direction. In the collinear limit the spin wave function becomes Multiplying the spin wave function with the reduced wave function and integrating over k ⊥ , one arrives at the associated distribution amplitudes. The first term in (18) generates the twist-2 distribution amplitudē Because of charge conjugation invariance the twist-2 distribution amplitude is antisymmetric in ξ. It possesses a Gegenbauer expansion and depends on the factorization scale, µ F , [10,27,28] Evidently, the reduced wave function must be symmetric in ξ. The Gegenbauer coefficients in (20) which encode the soft, non-perturbative QCD, evolve with the factorization scale as where Here, β 0 = (11N c − 2n f )/3, C F = 4/3 and n f denotes the number of active flavors. For the initial scale, µ 0 , the value 1.41 GeV is chosen in this article.
The decay constantf 0 depends on the scale too [10] f The other terms in (18) are of twist-3 nature although they do not correspond to the full twist-3 contributions since, in general, they also receive contributions from a second reduced wave function. This is however of no relevance for the purpose of the present paper, namely the calculation of the γ * − f 0 transition form factors. As we shall see in the following there is no twist-3 contribution to it. Anyway the k ⊥ -integration of the other terms leads to two further distribution amplitudes which are related to Φ 0 in the case at hand: Both these distribution amplitudes are symmetric in ξ and only the even terms appear in their Gegenbauer expansions.
With the help of these distribution amplitudes one can transform the product of wave function and collinear spin wave function (18), integrated over k ⊥ , into the form This expression resembles the corresponding pion spin wave function to twist-3 accuracy [29,30].

A wave function for the f 0 meson
For the evaluation of the transition form factors the light-cone wave function is to be specified. It is modeled as a Gaussian in k 2 with This wave function is similar to the one for the pion advocated for in [18]. It has been used for instance in the calculation of the photon-pseudoscalar transition form factors [31] or in the analysis of pion electroproduction [32]. Insertion of the wave function into Eq. (19) leads to the associated distribution amplitude (20) with the Gegenbauer coefficients (m is an odd integer) As a consequence of charge conjugation invariance which forces Ψ 0 to be symmetric in ξ, the matrix element vanishes in accord with the result quoted in [10]. On the other hand, the scalar density provides Evaluation of the integral leads toB 0 ≃ −1. This estimate is to be taken with caution since Φ 0s in (24) is likely not the full twist-3 distribution amplitude, but it provides orientation. As is obvious from the vacuum-particle matrix element of quark field operators given in (30), the decay constant is a short-distance quantity; it represents the wave function at the origin of the configuration space. It is also clear that only the ss Fock component of the f 0 -meson contributes to this matrix element.
For the numerical evaluation of the γ * − f 0 transition form factors the wave function will be restricted to the first Gegenbauer term, all others are neglected. with In this case the twist-2 distribution amplitude reads For the transverse-size parameter, a 0 , the value 0.8 GeV −1 is taken in the following. This value is very close to the corresponding value for the pion, see [31]. The r.m.s. k ⊥ is related to the transverse-size parameter by is adopted which has been derived by De Fazio and Pennington [33] from radiative φ → f 0 γ decays with the help of QCD sum rules (see also [34]). In [33] the f 0 -meson is considered as a (dominantly) ss state. The value (35) is extracted from the stability window for the Borel parameter between 1.2 and 2 GeV 2 . This is consistent with the initial scale chosen in this article. 3 The γ * − f 0 transition form factors

The definition of the form factors
Let us consider the general case of two virtual photons where q i and p denote the momenta of the photons and the mesons while λ i are the helicities of the photons. One has It is convenient to introduce the following variables [35] where, obviously, −1 ≤ ω ≤ 1.
The transition vertex is defined by the matrix element of the time-ordered product of two electromagnetic currents where and e i are the quark charges in units of the positron charge, e 0 . Following [4] the vertex is covariantly decomposed as where Current conversation is manifest: As one sees from (41) there are two form factors, one for transverse photon polarization, F T , and another one for longitudinal polarization, F L . By definition the form factors are dimensionless. p q 1 q 2 Figure 1: The equal-energy brick wall frame.
Contracting the vertex function with the polarization vectors of the photons and using transversality (ǫ i q i = 0), one arrives at One can show, most easily in the equal-energy brick wall frame (see Fig. 1), defined by q 1 = (ν 0 0 a 1 ) , q 2 = (ν 0 0 a 2 ) , p = (2ν 0 0 a 1 + a 2 ) , that the contraction with transverse photon polarization vectors with the same helicity projects out the form factor F T and with longitudinal ones F L : If the photons have different helicities the vertex function is zero.

The LO perturbative calculation
In the perturbative calculation of the form factors, performed at large Q 2 , the mass of the f 0 -meson is neglected whenever this is possible. From the Feynman graphs shown in Fig. 2 one finds for the vertex function (39) Figure 2: LO Feynman graphs for the γ * → f 0 transition form factors. The momenta of the virtual partons are denoted by g 1 and g 2 .
where the parton virtualities read (see also Fig. 2) Taking into consideration that the traces are only non-zero for even numbers of γ matrices, one notices that only the first term of the spin wave function (11), i.e. the leading-twist piece, contributes to the traces. The twist-3 terms lead to an odd number of γ matrices in the traces, the fourth term is neglected. With the help of (46) one finally arrives at the following expressions for the form factors: Because of the variation of ω 2 between 0 and 1 the mass dependent terms in front of the integral are kept. For ω → ±1 they exactly cancel whereas for ω → 0 F T ∝ Q −4 . The wave function (26) generates a factor (1 − ξ 2 ) 2 in the k 2 ⊥ integration. Hence, there is no singularity at the end points ξ → ±1 for all ω. One also notices from (49) that For (49) can be neglected and For ω → 0, on the other hand, only the term ∼ M 2 0 /(2Q 2 ) remains and Explicitly, for ω → 1 (i.e. Q 2 2 = 0) For a wave function of the type (26) one can write Eq. (53) as and For this type of wave functions the transition form factor is given by the collinear result multiplied by a universal reduction factor ̺. The latter function is shown in Fig. 3. It is interesting that, in the collinear approximation, the LO perturbative result for the form factor is related to the sum over all Gegenbauer coefficients. The γ − π transition form factor possesses this property too. This makes it clear that it is impossible to extract more than one Gegenbauer coefficient from the γ − f 0 transition form factor data. This coefficient is to be regarded as an effective one. NLO corrections may allow one to fix a second coefficient [35]. The situation improves for |ω| < 1 as will be discussed in Sect. 4.2.

The real-photon limit
The BELLE collaboration [1] extracted the γ − f 0 transition form factor from the cross sections on γ * γ → π 0 π 0 . In order to fix the normalization of that form factor the couplings of the f 0 to both the two-photon and the ππ channels are required. Both these couplings are not well known [8]. Hence, the normalization of the γ−f 0 transition form factor is subject to considerable uncertainties. The published data on the transition form factor, F T (Q 2 1 ), are scaled by the value of the form factor at Q 2 1 = 0 obtained from the width of the two-photon decay of the f 0 -meson (M 0 = (990 ± 20) MeV [8]) From the average decay width quoted in [8], one obtains The BELLE collaboration uses the slightly different value |F T (0)| BELLE = 0.0832 ± 0.0136. In a first step the BELLE data are compared to the collinear result  (32). For these wave function parameters the probability of the ss Fock component of the f 0 -meson is (see (14)): The results of the fit are shown in Fig. 4. Reasonable agreement with experiment is to be seen within rather large errors although the shape of the fit is opposite to that of the data: the collinear result for the scaled form factor, Q 2 1 F coll T slightly decreases with increasing Q 1 due to the evolution of the decay constant and the Gegenbauer coefficient, B 1 , whereas the data increase in tendency.
An increasing scaled form factor can be generated by quark transverse momenta in the hard scattering kernel and in the wave function, see Fig. 3. Retaining the quark transverse momenta implies that quarks and antiquarks are pulled apart in the transverse configuration or impact-parameter space. The separation of color sources is accompanied by the radiation of gluons. These radiative corrections have been calculated in Ref. [36] in the form of a Sudakov factor in the impact parameter plane. The Sudakov factor, e −S , comprises resummed leading and next-to-leading logarithms which are not taken into account by the usual QCD evolution. The k ⊥ -factorization combined with the Sudakov factor is termed the modified perturbative approach (mpa) [36]. It has been used, for instance, in calculations of the pion electromagnetic form factor [36] or the π − γ transition form factor [31] and will be used here as well. In the impact-parameter plane the transition form factor (53) reads The integrand is completed by the Sudakov factor, exp (−S), its explicit form can be found for instance in [31]. The Sudakov factor provides the sharp cut-off of the b-integral at 1/Λ QCD . Since 1/b in the Sudakov factor marks the interface between the non-perturbative soft momenta which are implicitly accounted for in the meson wave function, and the contributions from soft gluons, incorporated in a perturbative way in the Sudakov factor [36,31], it naturally acts as the factorization scale. The Bessel function K 0 is the Fourier transform of the hard scattering kernel and k 2 ⊥ Ψ 0 is the Fourier transform of the wave function (31) multiplied by k 2 ⊥ . It reads Evaluating the form factor within the modified perturbative approach and fitting B 1 to the BELLE data [1] one arrives at the results shown in Fig. 4. The fit provides the following value for the Gegenbauer coefficient 7 and χ 2 = 5.9 for 9 data points. The normalization uncertainty of the theoretical result follows from the errors of B 1 and F T (0), see (58). The agreement of the result obtained within the modified perturbative approach, with experiment is somewhat better than for the collinear approximation -the scaled form factor increases with Q 2 1 as the data do. This increase is the effect of the k ⊥ corrections shown in Fig. 3, the Sudakov factor plays a minor role in this context 8 . In passing it is noted that the predictions presented in [2] lie markedly below experiment for Q 2 1 > ∼ 10 GeV 2 . The value (62) of the Gegenbauer coefficient B 1 is not far from the QCD sum result [10]: The coefficent B 3 is found to be zero within errors. However, the value off 0 is substantially larger than the value (35) used in the form factor calculation. More precisely, the fit to the BELLE data fixes the product off 0 and B 1 for which the following results exist The product off 0 and B 1 derived in [10] is substantially larger than the BELLE data [1] on the γ − f 0 transition form factors allow. This product of f 0 and B 1 is also in conflict with a light-cone wave function interpretation since it leads to a probability larger than 1. Of course a smaller value of the transverse-size parameter would cure this problem for the prize of an implausible compact valence Fock component. For instance, if one halves a 0 the probability is about 0.12 but k 2 ⊥ ≃ 1.2 GeV. The last issue to be discussed is the contribution from the non-strange qq Fock state to the γ − f 0 transition form factor. This is usually considered as f 0 −σ mixing [9]- [15]. As for the η−η ′ system [37] this mixing is treated in the quark-flavor basis. As a consequence of the smallness of OZI-rule violations η − η ′ mixing is particularly simple in that basis -there is a common mixing angle for the states and the decay constants. It is assumed that this mixing scheme also holds for the case of interest here. Let σ n and σ s be states with the lowest Fock components nn = (uū + dd)/ √ 2 and ss, respectively. In analogy to (30) the corresponding decay constants are defined by the σ ivacuum matrix elements of the quark field operators: Since in hard processes only small spatial quark-antiquark separations are of relevance it seems plausible to embed the particle dependence and the mixing behavior of the qq Fock components solely into the decay constants 9 (for a detailed discussion of this procedure in the η − η ′ case see [38]). In generalization of (30) one may also define the decay constantsf q These decay constants mix according tō According to [10,39]f n ≃f s . Since the decay constant quoted in (35) is to be identified withf s 0 and since | cos ϕ| is close to 1, see (1), it suffices to assumef n ≃f s ≃f 0 for a rough estimate. For the range ϕ = (25 − 40) • of the mixing angle quoted in (1) one finds Clearly, this leads to a transition form factor which is in conflict with the BELLE data [1]. Using the second range of mixing angles in (1) one obtains reasonable agreement with experiment. Particularly favored is the range ϕ = (145 − 151) • for which the form factor stays within the uncertainty band displayed in Fig. 4. An exact determination of the mixing angle is not possible at present given the poor information available for the basic decay constants,f n ,f s , and the assumption on the explicit form of the light-cone wave function.

The case of two virtual photons
Here, in this subsection, it will be commented on the γ * − f 0 transition form factor. As is the case for ω = 1, the Sudakov factor plays a minor role. In order to estimate the importance of the power corrections taken into account in the modified perturbative approach the ratio of the form factors evaluated from (49) (transformed to the impact parameter plane and with the Sudakov factor included) and from the collinear approximation is displayed in Fig. 5. As expected the power corrections become smaller with increasing Q 2 and their importance decreases if ω deviates from 1. The same observation has been made in [35] in case of the γ * − π transition form factor. As noticed in [35] the reason for this effect is the term 1 − ξ 2 ω 2 in the hard scattering kernel which controls to which extent the form factor is sensitive to contributions from the end-point regions ξ → ±1 where soft effects can be important.
Since the power corrections are small at small ω, it is of interest to look at the transition form factor (71) in this region. Using the Gegenbauer expansion of the distribution amplitude the integral can be carried out term by term. The full result is a power series in ω 2 leaving aside the ω-dependence of the prefactor. The first terms of this series read As one notices the m-th Gegenbauer coefficient comes with the power ω m−1 first. For Q 2 larger than 4 GeV 2 the difference between the modified perturbative approach and the collinear result is smaller than 10%. Hence, the result in the modified perturbative approach evaluated from the wave function (26), is not far from the collinear result (72). Thus, as is the case for the γ * − π transition form factor [35], a measurement of the γ * − f 0 transition form factors for a range of small ω would therefore provide valuable constraints on the f 0 distribution amplitude. In Fig. 6 the γ * − f 0 transition form factor, evaluated from the wave function (31) within the modified perturbative approach, is shown for several small values of ω. It is clearly seen that the form factor drops with Q 2 increasingly stronger than 1/Q 2 with decreasing ω. At ω = 0 it decreases as 1/Q 4 (aside from evolution logarithms).

Summary
In this article the spin wave function of the f 0 (980) meson is constructed under the assumption that the meson is dominantly a strange-antistrange quark state. The collinear limit of the spin wave function is also discussed and the connection to the twist-2 and twist-3 distribution amplitudes is made. The spin wave function is applied in a calculation of the γ * − f 0 transition form factors. In the real-photon limit the results for the transverse form factor are compared to the large-Q 2 data measured by the BELLE collaboration recently. It turns out that, for the Q 2 range explored by BELLE, the collinear approximation does not suffice, power corrections to it, modeled as quark transverse moment effects, seem to be needed. The parameters required in this calculation in order to achieve agreement with BELLE form factor data, the transverse-size parameter, a 0 , the decay constant,f 0 , and the lowest (effective) Gegenbauer coefficient, B 1 , have plausible values. However, Cheng et al [10] in their analysis of charmless B-meson decays, adopt a much larger value forf 0 than (35). It remains to be seen whether the B-meson decays can be reconciled with the decay constant (35). The implications of σ − f 0 mixing for the transition form factors are also briefly discussed. A mixing angle of about 150 • seems to be favored. The paper is completed by presenting results on the γ * − f 0 form factors and on their collinear limits. It turns out that, in many aspects, the photon -f 0 form factors have properties similar to the form factors for the transition from a photon to the π 0 or other pseudoscalar mesons. However, the limits for Q 2 1 → ∞ are different. Whereas for the pseudoscalar mesons the limits of the scaled form factors are finite (e.g. Q 2 1 F γπ 0 → √ 2f π ) the γ − f 0 form factor F T tends to zero ∼ f 0 (µ 0 )B 1 (µ 0 )(α s (µ 0 )/α s (Q 2 1 )) −4/25 . The γ * − f 0 transition form factors also play a role in the calculation of the hadronic light-by-light contribution to the muon anomalous magnetic moment [40] - [43]. In particular, the results presented in this article clarify the asymptotic behavior of the γ * − f 0 form factors.