The anomalous decay f_{1}(1285)\to\rho\gamma and related processes

The anomalous decay and related processes

A. A. Osipov111Email address:, A. A. Pivovarov222Email address:, M. K. Volkov333Email address: Joint Institute for Nuclear Research, Bogoliubov Laboratory of Theoretical Physics, 141980 Dubna, Russia

We work out the low-energy expansion of the anomalous decay amplitude by using the Nambu–Jona-Lasinio model with chiral symmetric four-quark interactions in the one-quark-loop approximation. The related processes , , and , are also considered. An effective meson Lagrangian responsible for , , and interactions is found. The predicted radiative decay widths, , , , , allow an experimental test of the hypothesis that and -mesons have a quark-antiquark nature. At present, only the decay has been measured. Our result is in remarkably good agreement with the recent data of CLAS Collaboration , but disagrees with the PDG-based estimate of . The calculations presented require a minimum of theoretical input, and are shown to be consistent with the non-renormalization theorems of QCD.

11.30.Rd, 11.30.Qc, 12.39.Fe, 13.40.Hq

I Introduction

Anomalies have important consequences for a wide range of issues in quantum field theory. This is why they are invariably under special attention of theoreticians and experimentalists. For instance, the Wess-Zumino effective Lagrangian Wess71 () summarizes the effects of anomalies in current algebra and finally relates uniquely the decay amplitude with other ones, such as , and a five pseudoscalar vertex. Thus, the anomaly based results are tightly restrictive and potentially very accurate. The latter is the consequence of the Adler-Bardeen theorem Adler69 () which states that chiral anomaly is not modified by higher order corrections.

The world average phenomenological data on the radiative decay PDG () selected through the study of the reaction , and new results of CLAS Collaboration at Jefferson Lab. on the photoproduction off a proton target Dickson16 () are especially important due to a presence of the anomaly: the underlying triangle quark loop diagrams describing the and vertices. These vertices determine the widths of the and decays and are basic elements in the description of the photoproduction data. In this respect the presently available phenomenological data allow a sensitive test of the anomaly. Yet some of the essential properties of the theoretical description of this vertex are only poorly understood.

The connections of vector and axial-vector mesons with the anomaly can be studied on general grounds, i.e without assuming the quark-antiquark structure for spin-1 states. For instance, the method, based on the massive Yang-Mills approach Schechter84 (), leads to the Bardeen’s form of the non-abelian anomaly. Unfortunately, this breaks explicitly the chiral symmetry and forbids the , , , and vertices.

On the opposite, if one starts from the most general anomalous action in terms of pseudoscalars and spin-1 states which is chirally (gauge) symmetric and embodies the chiral anomalies only through the Wess-Zumino-Witten action of the pseudoscalars Witten83 (), one gets a consistent scheme, and in this case there is the possibility for anomalous , , , and vertices Kugo85 (); Meissner90 (). However, the method fails to predict the couplings of the effective Lagrangian, and cannot be used to estimate the width of the decay.

The purpose of the paper is to clarify exactly this obscure aspect of the radiative decay of the axial-vector ) meson. For definiteness we will consider the Nambu–Jona-Lasinio (NJL) model with chiral symmetry spontaneously broken down to the diagonal subgroup (the quantum anomaly breaks the axial symmetry) Ebert83 (); Volkov84 (); Volkov86 (); Ebert86 (); Ebert94 (); Osipov94 (); Bernard96 (); Volkov06 (). The NJL model not only gives the structure of the vertex, but also fixes the values of the coupling constants involved.

There are at least three essential reasons for our calculations. First, we show that if one assumes that the meson is a bound quark-antiquark state one can obtain its radiative decay width by considering the anomalous quark triangle diagram. The result is restrictive because there are general statements about the longitudinal and transversal parts of the triangle Adler69 (); Vainshtein03 ().

Second, the phenomenological data on this decay presently are very contradictive. The PDG-based estimate is . The recent CLAS data give three times less this value . We argue that the quite low value reported by the CLASS collaboration can be perfectly understood if the -meson is approximately a quark-antiquark state.

Third, it has been noted in Wang17 () that known theoretical models Kochelev09 (); Domokos09 () failed to fit new CLAS data. We suppose that one of the reasons is related with the oversimplified expressions used for and vertices. In this respect a new attempt made in Wang17 () also suffers from a superficial approach to the anomalous structure of these vertices. We suggest a different effective Lagrangian which consistently describes the triangle anomaly and can be applied to fit CLAS data.

Ii Anomalous triangle diagram for the decay

The relevant vertices of the NJL quark-meson Lagrangian density are


where the are constituent quarks of mass (the color and flavor indices are suppressed); is an axial-vector field; is a neutral -meson field; is a photon field; is the matrix of the light quark’s charges, where is a diagonal Pauli-matrix; are the standard Dirac matrices in four dimensions. The couplings and are well established: is the proton charge with and is fixed from the decay, . The constituent quark mass is equal to Volkov86 ().

The amplitude of the decay can be written as follows


where is a number of colors, is the polarization vector of the -meson, is a polarization vector of -meson and is the polarization vector of the photon; and are 4-momenta of and photon correspondingly, is the 4-momentum of -meson. The tensor is a sum of two Feynman diagrams (see fig.1)

Figure 1: The one-quark-loop contribution to the radiative decay amplitude . The first diagram corresponds to the Lorentz tensor , the second one to the tensor . External lines represent the axial-vector field with polarization vector , the vector -meson field with polarization vector , and electromagnetic field with polarization vector .



It is obvious from (5) that the Lorentz tensor obeys the Bose symmetry requirement


Making a replacement of variables in one of the integrals in (5), and calculating traces, we obtain that .

After some mildly tedious calculations we find that the amplitude can be written in the form


where we have restricted ourselves up to the terms of the third power in momenta. (To describe correctly the low-energy limit, the amplitude must have the smallest possible number of momenta. One should not think about this truncation in terms of expansion which is not applicable here. We are following instead the idea of expansion. According to it, the meson physics in the large limit is described by the tree diagrams of an effective local Lagrangian, with local vertices and local meson fields Witten79 (). This is exactly what one obtains restricting to the leading in momenta terms of the constituent quark loops. The details of such description of low-energy meson physics in the framework of NJL model are given in Volkov84 (); Volkov86 (); Ebert86 (). In particular, our result (II) differs from the one obtained by Kaiser and Meissner Meissner90 () only by prescribing the definite values to the corresponding couplings of the effective -vertices in accord with the NJL model. In Sect. V we show that the truncated triangle diagrams of Fig.1 do reflect the QCD anomaly structure, in the infrared.)

The result (II) contains an ambiguous surface term, represented by the 4-vector . It is a well-known remain of superficial linear divergence of the quark-loop integral Jackiw72 (). Most generally , where and are free constants. The property (7) relates these constants . The requirement of gauge symmetry is . That totally fix the constants. Indeed, we get from (II)


and if one takes the Ward identity (9) is obviously fulfilled. We conclude that the surface term is completely fixed by Bose and gauge symmetry requirements, . As a result, the linear in momenta contribution in (II) is zero.

We can gain some deep understanding of this formula by considering the product . It is easy to find out that the result is zero. It means that if one considers the transition of in our amplitude (in accord with the idea of vector meson dominance) the formula (II) obeys the additional requirement of gauge symmetry.

Moreover, the celebrated Landau-Yang theorem Landau48 () states that a massive vector (i.e. spin-1) particle cannot decay into two on-shell massless photons. Let us show that our amplitude does not contradict this general result. Indeed, in a frame where meson is at rest, we can always choose the direction of the -axis along the spatial part of the photons momenta, i.e. and . The photon polarization vectors and are orthogonal to the photon momenta, and thus can be chosen as follows: and . The polarization vector of a massive meson is given by . It follows then that


Hence, the amplitude of decay is equal to zero in the boson rest frame. Since the amplitude is Lorentz invariant, it equals zero in any other frame as well.

Now, it is not difficult to relate the amplitude (4) (with given by (II)) with the Lagrangian density, describing this anomalous decay and leading to the same amplitude


Here the quantities stand for the field strengths associated with neutral vector -meson field , , electromagnetic field , , and the neutral axial-vector field, . This expression gives a definite meaning to our statement about the oversimplified form of the Lagrangian used in the literature for this vertex. To see the difference it is enough to compare (II) with the vertex used, for instance, in Wang17 (), where .

On the mass surface of and mesons from (4) and (II) we get


where . Then it follows that


and the radiative decay width is given by


This model estimate is in a perfect agreement with the experimental result, given by the CLAS Collaboration: , and about four times less the PDG estimate: PDG ().

Iii The related decays

Our purpose now is to describe the related processes, i.e. the radiative decays , , and . Their amplitudes are originated by the same quark-loop integrals as the amplitude (4). Therefore, the general factor which comes out from the isotopic trace calculations will be the only difference in the results. Let us remind that for the amplitude this factor is . Now, the corresponding factor in the amplitude is . It gives immediately the Lagrangian density


and the radiative decay width


The amplitude has the factor . It means that


and the radiative decay width is equal to


The radiative decay has three different channels: , and . Due to the property , we can sum the traces over Pauli matrices of these two contributions. That gives


We conclude that each of the three possible modes has a similar amplitude and the same expression for the decay width, i.e.


and , with


Iv Comparison with other approaches

No much efforts have been done till now for a theoretical description of the processes considered. Probably, this is related to the very poor experimental information available on these radiative decays. In Table 1 we collect the relatively old estimations made in the framework of the covariant oscillator quark model Ishida89 (). Our results, in general, have a tendency to be twice smaller of that predictions. Let us remind that the isoscalar member of the axial-vector -nonet (), -meson, has mostly content, but can mix with the mainly isosinglet state (and gluons). The authors of Ishida89 () considered two possible candidates for such partner: the usual candidate, , and the another promising state


In the case of the combination of two isoscalar members, and , the mass formula of their model gives the mixing angle . In the case of the other combination, and , the mass formula gives . One can see that data on and modes in that model slightly depend on the mixing angle . The tendency is the smaller the more radiative decay width (see Table 1).

In our work we consider the -meson as a pure non-strange state (). Presently there are some indications that such mixing in the axial-vector nonet is really small and is mostly made of and quarks. For instance, LHCb Collaboration Aaij14 () gives , assuming that is mixed with the state. This agrees with an earlier determination of in Gidal87 ().

311 34.3 238 26.8
Dickson16 ()
Amelin95 ()
PDG ()
Ishida89 () 509-565 48-57 537 62
Table 1: NJL-model predictions for anomalous radiative decays of and axial-vector mesons, . We give also some known empirical data, , and predictions of the covariant oscillator quark model Ishida89 (). All decay widths are given in keV.

Calculating the triangle diagrams we considered the lowest order terms in an external momenta expansion (minimal couplings). It means that we are only concerned with the part of the effective action having the smallest possible number of derivatives which is responsible for the intrinsic parity violating processes. This approximation corresponds to the standard counting of the spin-1 mesons in Resonance Chiral Theory. In that approach, they contribute at order of chiral counting in the effective meson Lagrangian (i.e. with terms kept up to four derivatives). For further arguments supporting this approximation, we refer to the original papers Ecker89 (); Gasser89 (); Donoghue89 (). Nonetheless, due to the importance of the question, we give some additional arguments in the following two sections.

It is worth to be mentioned here the quark model predictions by Lakhina and Swanson (see ref. [52] in Dickson16 ()). They have found that a nonrelativistic Coulomb-plus-linear quark potential model predicts , while a relativized version of the model gives much less value . One sees that the NJL model result based on the relativistic quantum field theory calculations is in agreement with the relativized version of the Lakhina and Swanson model. Both models nicely reproduce the CLAS Collaboration result. On the other hand, the PDG-based estimate favors the nonrelativistic result.

V Restrictions from QCD

The anomalous quark triangle diagrams considered here are the subject of special attention in the literature. The pioneering studies have been done by Rosenberg Rosenberg63 () and AdlerAdler69b (). Rosenberg got an explicit expression for the fermion triangle graph:


where we follow his notations


Here is a mass of the fermion in the triangle. The effective Lagrangian corresponding to (25) is, strictly speaking, local only if (28) can be well approximated by . In this case, we obtain from above


This is exactly the approximation that has been used in our estimates of (in our case is the constituent quark mass). Locality of the model requires us to work with form-factors and considered at and , i.e. off mass-shell of variables and .

It is, however, not well-understood what is the accuracy of such a step. Some interesting insight can be obtained by considering the amplitude in the limit of small external photon momentum Kuraev92 (). In this limit the tensor is linear in (one neglects quadratic and higher powers of ), and we obtain from (25)


Due to a special kinematics, which corresponds now to the decay of the axial-vector state in flight, the expressions for the form-factors and of the fermion triangle graph are considerably simplified


but still follow the pattern . In general these two form-factors are independent.

The equation (30) can be cast into the standard form with the aid of the Schouten’s identity Schouten38 ()


where . This identity allows us to write


where the invariant functions , and are the longitudinal and transversal parts of the quark triangle with respect to the axial-vector momentum . Both structures are transversal with respect to the photon momentum , .

Let us consider now the problem from a different angle, namely, by using the one-quark-loop QCD result for the triangle graph. In this case, the longitudinal and transversal form-factors are still given by the equations (V), where one should only replace the constituent quark mass by the current quark mass . The result is (for )


In particular, in the chiral limit, , one gets


The longitudinal part represents the axial anomaly associated with the divergence of the axial-vector current (where for the meson case, or for the meson case) constructed from the light quark fields. Indeed,


where , and . Note, that the imaginary part of is equal to zero since Dolgov71 ().

The Adler-Bardeen theorem Adler69 () implies that the one-loop result (35) for stays intact when the interaction with gluons is switched on. It is not corrected at the nonperturbative level too. Moreover, for the special kinematics considered (the photon carries a soft momentum) the transversal part is unambiguously fixed by the longitudinal one in the chiral limit of perturbation theory (up to the nonperturbative corrections to ) Vainshtein03 ().

The model result (II) taken in special kinematics can be compared with these general QCD statements. Indeed, one can see that (a) it follows the QCD pattern ; (b) it does not contain the imaginary part, which otherwise would contribute in the chiral limit to the divergence of the axial-vector current; (c) it correctly reflects the underlying anomaly. The latter needs some explanation. On the mass-shell of the -meson, we have


Numerically the factor agrees quite well with factor 1 in (V). (The difference in the overall sign is not essential here because the sign in (V) can be changed by the appropriate definition of the QCD AVV-currents correlator.) But it is valid to expect from (V) and (37) of the full agreement, because we are dealing with the anomaly. And this expectation is actually correct. To show this let us use the mass formulae of the model


It follows then, that the value would be in a total agreement with QCD requirement (V). Exactly this value of was obtained long ago on the basis of spectral-function sum rules Weinberg67 (), which are valid in QCD for , and of the KSFR relation for the coupling to the isospin current Suzuki66 (); Riazuddin66 (). Now, we come to the same conclusion on the basis of the anomaly consideration.

This reasoning, however, requires a caveat. From the QCD calculations, it follows that the axial anomaly has a pole at . This pole implies the presence of zero-mass bound states in the physical spectrum Frishman81 (); Achasov93 (). In the case of radiative decays, the pole in can be associated with pions. However, for the case, there is no such a light pseudoscalar state (the problem). For further progress with the singlet case, the gluon anomaly must be included. Some interesting attempts may be found in Efremov90 (); Teryaev90 (); Teryaev95 (), where the authors argue that the photon anomaly must be canceled by the gluon one. The question has been also addressed in the instanton liquid model of QCD vacuum Dorohov05 (), where it has been shown that the singlet longitudinal amplitude is renormalized at low momenta by the presence of the gluon anomaly. As a result, the product vanishes at , taking the value on the -meson mass-shell. On the contrary, the normalized nonsinglet amplitude follows the pattern (this value exactly corresponds to the factor 1 in eq.(V)). This shows that inclusion of the gluon anomaly will slightly change our estimates for the singlet decays. Nonetheless, the deviation from the equality used in our calculations is rather small for the physical region of considered here.

Vi The results

It is straightforward now to obtain in the approximation considered above the amplitude. The result is


where , are four-momenta of photons, is a four-momentum of the meson, , , are polarization vectors and


In particular, in the rest frame of , where , we obtain


Although the decay of a spin-one resonance is suppressed for real photons, according to the Landau-Yang theorem Landau48 (), the value of the matrix element (VI) is measured with a good accuracy Aihara88 (); Achard02 (). For instance, an estimate Kochelev17 () based on L3 Collaboration data gives


A similar result one obtains from PDG PDG () values


Our estimate,


agrees well with these phenomenological values and gives an additional argument in favour of approximations made in the previous sections.

Vii Conclusions

In this work, we evaluate the radiative decay widths of the and axial-vector mesons with the assumption that these states are mostly made of and quarks. The anomalous triangle diagrams are considered. We show that surface term of these diagrams are fixed by the vector Ward identity and Bose statistics. These leads to the total cancellation of the linear in external momenta contributions in the low-energy expansion of the amplitudes. As a result, the leading order contribution (in the external momenta expansion) is determined by the cubic in momenta terms. This finite contribution in the framework of NJL model is totally fixed by the coupling of the decay, , the finite structure constant , and masses of light quarks and mesons. Note, that the constituent quark mass in the NJL model can be related with the phenomenological parameters only. The relation is given by the formula


The effective Lagrangian describing the radiative decays of these mesons are obtained.

The decay width found,