π0–η–η' Mixing in the Theory with Four-Quark Interactions

Mass formulas, mixing angles, and decay constants of the pseudoscalar π0, η, and η' mesons have been obtained in the theory with four-quark interactions with an accuracy to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal{O}(1{\text{/}}N_{{\text{c}}}^{2})$$\end{document} terms inclusively. Unlike the standard Nambu–Jona-Lasinio model, the meson Lagrangian has been obtained using the Volterra series, which allows a more detailed description of effects caused by the inequality of quark masses. The results have been compared to similar calculations in the 1/Nc chiral perturbation theory. The first corrections to the interaction violating the Zweig formulas have been calculated, and it has been shown that they are important for the description of the mass spectrum of the η–η' mesons.


INTRODUCTION
To study the properties of the pseudoscalar nonet of mesons, Leutwyler [1,2] used the Lagrangian whose effective vertices are classified in terms of the powers of momenta, masses of light quarks, and 1/N c , where N c is the number of color degrees of freedom. The method called the 1/N c chiral perturbation theory [3] made it possible to calculate the first correction to the main result of the current algebra for mass formulas of the charged pseudoscalar mesons with nonzero strangeness, to show that the correction is small (at that time, contradictions existed in the description of the η → 3π decay), and to obtain constraints on the masses of light quarks. Later, the method was extended to spin-1 states [4]. Among recently results obtained with this effective Lagrangian, we mention the analysis of the two-photon decays of the π 0 , η, and η' mesons including the first-and second-order corrections in 1/N c [5].
I have recently shown [6] that the theory with fourquark interactions of the Nambu-Jona-Lasinio type [7,8] gives the same mass formulas for the pseudoscalar , , , and mesons as in [1]. This theory allows one to relate the parameters of the effective Leutwyler theory to the parameters of the dynamic Nambu-Jona-Lasinio model. This is achieved by expanding the effective action of the Nambu-Jona-Lasinio model in the Volterra series in inverse powers of the masses of the constituent quarks rather than in the Taylor series in powers of the proper time [9][10][11]. In addition, the Leutwyler hypothesis on the behavior of the masses of light quarks ( ) in the limit of large N c values is used.
The aim of this work is to study the physical characteristics of the remaining terms of the nonet-π 0 , η, and η' mesons-for which the masses and mixing angles are calculated with an accuracy to the first correction in 1/N c inclusively and the constant of the anomalous violation of the symmetry, as well as the degree of violation of the Okubo-Zweig-Iizuka rule, is estimated. The latter estimates are of particular interest in the context of the current active study of the gluon structure of the η and η' mesons from different points of view (see, e.g., [12], where the dispersion approach was applied to study the axial anomaly). The formalism used in this work, as well as the numerical parameters of the model, is the same as in [6].

MODIFIED NAMBU-JONA-LASINIO
MODEL The initial four-quark Lagrangian used in this work is the same as in the standard quark version of the Nambu-Jona-Lasinio model [13][14][15][16] (1) Here, are the Dirac matrices, q are the quark fields, and m = diag(m u , m d , m s ) is the diagonal matrix consisting of the current masses. The Lagrangian density has the form , where the first and second terms are U(3) L × U(3) R chirally symmetric combinations describing four-quark interactions with spins zero and unity, respectively, given by the expressions

FIELDS, PARTICLES, AND NUCLEI
Here, , are the Gell-Mann matrices, and the constants G S and G V in the limit N c → ∞ have the order ; their values determined in [6] are G S = 6.6 GeV -2 and G V = 6.8 GeV -2 .
The functional integral method allows the transformation of the Lagrangian density (1) to the form (4) Here, the functional freedom of choice of dynamic variables was used in favor of the nonlinear realization of chiral symmetry; , , , and are the Hermitian matrices describing the vector, axial vector, scalar, and pseudoscalar fields, respectively; , where P L = (1 -γ 5 )/2 and P R = (1 + γ 5 )/2, are the quark fields belonging to the fundamental representation; The pseudoscalar field φ is dimensionless and acquires the necessary mass dimension at the transition to the field functions of physical states. The matrix M = diag(M u , M d , M s ) is the diagonal matrix consisting of the masses of constituent quarks Q. These masses appear through the dynamical breaking of symmetry and are related to the masses of light quarks by the gap equation (8) where (9) Here, Λ is the cutoff parameter characterizing the scale at which the effective theory under consideration is studied. In this case, this is the scale of hadron masses, which is Λ ≈ 1.1 GeV.
The gap equation is a mathematical minimum condition for the effective potential; this potential and kinetic part of the effective action are obtained after integration over the quark fields Q. This integration gives the quark determinant with the local part described by the first terms of its asymptotic expansion in powers of the proper time. The difference between the quark masses leads to a problem concerning the inclusion of difference effects. To solve this problem, an asymptotic expansion in the inverse powers of large masses was proposed in [11] on the basis of the Volterra series. This expansion distinguishes our calculations from the standard approach involving the Taylor series. This difference is significant because the Volterra series contains numerous finite (in the limit ) vertices that vanish in the limit of equal quark masses. These vertices carry important additional information on violations of the isospin and flavor symmetries that is absent in the standard meson Lagrangian of the Nambu-Jona-Lasinio model.
If the quarks are massive, pseudoscalar fields are mixed with axial vector fields. To remove this mixing, it is necessary to redefine axial vector fields as [17] (10) where is a matrix and stands for the Hadamard product of matrices [18] defined as without summation over repeated indices.
Only the diagonal elements of the matrix that are required here have the form (11) where (12) The quark determinant includes the kinetic terms of the Lagrangian of free meson fields. They have the canonical form after the corresponding redefinition of the variables (13) where the relation of the components to the components , as well as the relation between the components and , is standard. As easily seen, the new field variables marked by the superscript R have the dimension of mass because the constants are given by the expression As a result, The mass part of the Lagrangian (4) is diagonal in the flavor basis However, the physical π 0 , η, and η' mesons are not pure flavor states. For this reason, it is convenient to pass to the singlet-octet components 0, 3, 8. In this basis, the off-diagonal elements of the mass matrix are zero only in the case of the exact SU(3) f symmetry. Thus, we obtain (17) where is the symmetric matrix with the elements (18) Here, it is necessary to take into account two important circumstances-anomaly and violation of the Okubo-Zweig-Iizuka rule-both explained within the 1/N c expansion [19][20][21]. Their leading term has the order coinciding with the order of the leading term in Eqs. (18), where . The Lagrangians corresponding to these processes have the form of the product of two traces. At the quarkgluon level, such a contribution comes from diagrams with quark loops coupled through the gluon exchange.
The Lagrangian breaking the symmetry was obtained in [22] (see also [23]). Using this result, we set (19) The dimensional constant is fixed from the physical masses of the pseudoscalar mesons, and it does not change in the limit N c → ∞, having the order where is the dimensional constant and χ = 2BM is the matrix having the form The quadratic part of the Lagrangian (20) leads to mixing (22) Nonphysical fields and should be transformed to physical ones according to Eq. (13). In particular, Consequently, beyond the leading approximation in Eq. (22), additional mixing between the neutral components is induced by breakings of the isospin and SU(3) f symmetries.
Further, we consider the first terms of the 1/N c expansions of the resulting formulas. Here, we do only the first two steps, i.e., represent the elements of the mass matrix of the physical states in the form of the sum of the leading order (LO) term having the order and the next-to-leading order (NLO) term with the order , which is the first correction to the LO term. The calculation of the second correction (NNLO term) requires additional consideration of meson one-loop contributions, which is beyond the scope of this work.

1/N c EXPANSION
The six model parameters Λ, , , and have the orders and . The first parameter specifies the characteristic energy scale and the other parameters are much smaller than the first parameter. This property allows the systematic 1/N c expansion of the effective theory.
Let us seek the solution of the gap equation (8) in the form (24) where is the solution of the equation at large N c values. We find (25) Here and below, the index 0 means that this function of the quark masses is calculated in the limit . In particular, and . In turn, the constants given by Eq. (14) have the form (26) where (27) The elements of the mass matrix including all above contributions can be represented in the form = , where the first term is the leading contribution and the second term is the first correction to it. The leading contribution has the order and is described by the formulas (28) where (29) This result coincides with the formulas obtained by Leutwyler [1], but in the case under consideration, all parameters except for and are related to the main constants of the four-quark dynamics.
Mixing of with and is due to the breaking of isospin symmetry. In the first order in the mass difference , this mixing is removed by rotating to small angles and , respectively. Mixing of the components and is due to the breaking of SU(3) f symmetry and can be removed by rotating to the angle θ. With an accuracy to the first order in the breaking of isospin symmetry, the transformation of the neutral components to the physical π 0 , η, and η' states has the form  Since Eqs. (28) contain only two unknown parameters and , they can be determined from the known masses of the η and η' mesons. To this end, we first determine the masses of the quarks using the physical masses of the , , and mesons, which in this approximation are given by the current algebra formulas and Here, the overline means that these expressions were obtained disregarding the electromagnetic interaction. The inclusion of this interaction increases the masses of the charged states: It is known that the difference between the masses of the charged and neutral pions is due primarily to the electromagnetic interaction. The contribution of the strong interaction is proportional to and is thereby negligibly small. Using the Dashen theorem [24], we obtain m u = 2.6 MeV, m d = 4.7 MeV, and m s = 95 MeV from the above mass formulas. Then, the parameters GeV 2 and are determined from the system of two equations for the masses of the and mesons. The mixing angles are , , and . The estimates presented above show that the Zweig rule is strongly violated in the leading approximation, as also mentioned in [1,25]. Only paying this price can one satisfactorily describe the mass spectrum of the η-η' mesons. In the absence of the Lagrangian (20), the mass of the η meson is much smaller than its phenomenological value and the angle θ is . The interaction violating the Zweig rule allows one to reproduce the mass spectrum but reduces the absolute value of the angle θ. The result for the angle θ is close to [26] and exactly coincides with the  result obtained in [27] with the inclusion of the first correction in 1/N c chiral perturbation theory. Now, we do the next step and calculate the first correction to the leading term. This correction includes the contributions from Eqs. (18) and (22). As a result, Corrections caused by Eq. (18) with an accuracy to a common factor coincide with the result of 1/N c chiral perturbation theory [27]. The correspondence between factors has the form (36)

Corrections
were not considered in [27]. As before, the second-order terms in the breaking of isospin symmetry were neglected when deriving Eq. (35). We recall that the mass matrix in this approximation is still diagonalized by the orthogonal transformation (31).
To obtain numerical values, we first determine the masses of light quarks using the following mass formulas for the , , and mesons, which also include the first correction to the current algebra result:

CONCLUSIONS
To summarize, the main characteristics of the π 0η-η' system in the model with four-quark interactions have been calculated. A new method for the lowenergy expansion of the quark determinant with virtual particles with different masses has been used to bosonize the quark vertices. Assuming that the masses of light quarks vanish in the limit N c → ∞, we have η → π π π 0 Δ ± = (44. OSIPOV obtained mass formulas and decay constants of these pseudoscalar states in the form of the first two terms of their expansion in series in 1/N c . As a result, we have described the mass spectrum of the η-η' mesons and have analyzed the degree of the violation of the Zweig rule. The estimate obtained in this work indicates the moderate violation of the Okubo-Zweig-Iizuka rule. We emphasize that the considered corrections to the interaction violating the Zweig rule have not yet been studied. In particular, at least the effective Lagrangian in [32] includes interactions violating the Zweig rule but only in the leading order in 1/N c .