<<

Effective Chiral Theory for Pseudoscalar and Vector

E.Gedalin∗, A.Moalem† and L.Razdolskaya‡ Department of Physics, Ben Gurion University, 84105,Beer Sheva, Israel We consider the vector mixing scheme and mass splitting within the framework of an extended U(3)L U(3)R chiral effective field theory based on the hidden local symmetry approach where, the pseudoscalar and nonets play the role of dynamic variables. Unlike other variants of thisN model, we show that the diagonalization of the vector meson mass matrix and the assumption that its eigenvalues are identical with the physical meson masses, determines the mixing scheme as well as the model free parameters. We show that VMD can be derived from our lagrangian and that for electromagnetic processes at low momenta the VMD model is a good first approximation. The model reproduces nicely the radii of the charged and .

I. INTRODUCTION

In a recent contribution [1] it was indicated that, for any effective field theory (EFT) of colorless meson fields, the mixing schemes of states and decay constants are determined in a unique way by the kinetic and mass lagrangian densities. In a general case, these densities are bilinear in terms of the intrinsic fields, involving nondiagonal kinetic and mass matrices. However, they can be reduced into a standard quadratic form by transforming the intrinsic fields into the physical ones in three consecutive steps. These steps include : (i) the diagonalization of the kinetic matrix, (ii) rescalling of intrinsic fields to restore the standard normalization of the kinetic term, and (iii) the diagonalization of the resulting mass matrix. In such case where the dimensions of the nondiagonal kinetic and mass submatrices are respectively, k k and n n, this procedure leads to schemes which involve [k(k 1)/2]+ [n(n 1)/2] mixing angles and k field rescalling× parameters.× The commonly used mixing schemes correspond− to specific− choice of the kinetic and mass matrices. In particular, η η0 mixing requires one mixing angle scheme, if and only if, the kinetic term for the intrinsic fields has a quadratic− form. Such are the traditional scheme [2] and the one proposed by Feldmann et al. [3–5], the so called flavor basis (QFB) scheme. For bilinear kinetic and mass lagrangian densities (i.e. with nondiagonal kinetic and mass matrices) mixing schemes involve two mixing angles similar to the ones proposed by Escribano and Fr`ere [6]. There have been several attempts to incorporate within a single effective lagrangian a wide range of different electroweak and strong processes [7–10]. Particularly interesting are the Bando, Kugo and Yamawaki (BKY) model [7], based on the hidden local symmetry (HLS) approach, and variants of this model [8,9], where the vector mesons ¯ (ρ, ω, K∗, K∗,φ) are identified with the dynamical gauge of the HLS nonlinear chiral lagrangian. It was indicated by Benayoun et al. [9], that the original BKY model neither reproduces the observed vector mass splitting nor preserves current conservation for the pseudoscalar octet-singlet sector. A variant of this model proposed by Bramon et al. [8] also does not predict the masses correctly, although it guarantees current conservation, explicitly. Yet another variant of this model [9] includes the η0 meson also, and to some extent cures these deficiencies. However, once the HLS lagrangian is extended to include the η0 meson, the trace of the pseudoscalar (and vector) nonet field matrix no longer vanishes, and therefore, a general expression for the lagrangian must include additional terms which are quadratic in the traces of the nonet field matrices. More recently, a most general lagrangian based on local U(3)L U(3)R symmetry was constructed [11,12], by combining the HLS approach with a general procedure of including the η0 meson into a chiral theory [13–16]. Both the symmetric parts of the lagrangian, as well as the symmetry breaking companions,N include terms quadratic in the traces of the field matrices. Likewise, the lagrangian includes an explicit U(3)L U(3)R symmetry violating mass term for the pseudoscalar mesons [14,16]. These additional terms influence the meson mass spectra and more importantly, by means of the diagonalization procedure of Ref. [1], not only the mixingN scheme is defined uniquely but also, taking the eigenvalues of the resulting mass matrix to be equal to the physical masses, it was possible to fix a number of the model free parameters. All observables are calculated from the lagrangian straightforwardly without additional assumptions. Albeit, the lagrangian constructed reflects the

[email protected][email protected][email protected]

1 fundamental symmetries of the QCD lagrangian and provides a coherent theoretical framework where the pseudoscalar and vector meson nonets are treated on equal footing. Our main interest in the present work is to further develop and explore the U(3)L U(3)R chiral theory of Ref. [12] by considering the vector meson sector. It is demonstrated that in this case the relevant model parameters can all be fixed through the diagonalization procedure of Ref. [1]. It has been illustrated alreadyN [7] that the vector meson dominance (VMD) model emerges as a natural consequence of the basic assumption of the HLS approach where vector mesons play the role of gauge bosons. Particularly, the VMD model becomes a dynamical consequence for a specific choice of the HLS model parameter a = 2. Detailed analyses of data argue for departures from this value [9,17], and that in addition to the vector meson poles in the charged meson form factors, a direct - coupling may also contribute significantly. More recently, Harada and Yamawaki [18] have shown that in an effective field theory based on HLS, VMD may be badly violated. We demonstrate in this work that VMD derived from the lagrangian constructed in Ref. [12] is a good first approximation for electromagnetic processes at low momenta. Our paper is organized as follows. In Section II we present briefly the lagrangian. In section III we discuss vector mixing scheme and mass splitting corresponding to the lagrangian presumed in the simple case of exact isospin symmetry and in the more evolved case of isospin symmetry breaking. In Section IV we consider the photon interaction with hadronic matter and discuss the relation between the lagrangian constructed and the traditional VMD model. We summarize and conclude in Section V.

II. THE EFFECTIVE LAGRANGIAN

As in Ref. [12], the lagrangian is written in the form, 1 L = L + L¯ + L + a(L + L¯ ) Tr(V V µν )+L + ... , (1) A A m V V − 4 µν WZW ¯ ¯ where LA(LA)andLV (LV ) are the symmetric parts (symmetry breaking companions) for the pseudoscalar and vector mesons, Lm is a most general U(3)L U(3)R symmetry violating mass term [14,16], LWZW stands for the well known Wess-Zumino-Witten term [19,20], and the ellipsis ”...” represents terms accounting for the regularization of one loop N contributions [22,23,13]. All terms in the expression above are constructed from gauge fields Vµ, and vector and axial-vector covariants defined as [21,24,25], i 1 Γ = ξ ,∂ ξ + ξ r ξ + ξl ξ , (2) µ 2 † µ 2 † µ µ † i   1  ∆ = ξ†,∂ ξ + ξ†r ξ ξl ξ† . (3) µ 2 µ 2 µ − µ Here r and l represent the  external gauge fields; r =v + a and l = v a ,withv and µ µ µ µ µ µ µ − µ µ aµ being the vector and axial vector external electroweak fields, respectively. The non-linear representation of the pseudoscalar nonet fields is taken as [13,14,16],

2 √2 2 1 U(P, η0 + F0ϑ) ξ (P, η0 + F0ϑ) exp i P + i η01 , (4) ≡ ≡ ( F8 r3 F0 ) with η0(x) being the pseudoscalar singlet and P (x) the pseudoscalar Goldstone octet matrix,

1 π0 + 1 η π+ K+ √2 √6 8 1 0 1 0 P = π− π + η K , (5)  √2 √6 8  − ¯ 0 2 K− K η8  − √6    with obvious notations. The vacuum angle ϑ(x) is an auxiliary field which renders the variable X(x) √6η (x)/F + ≡ 0 0 ϑ(x) to be invariant under U(3)L U(3)R transformations [21]. In this representation, the unimodular part of U contains the octet degrees of freedom while the phase detU =exp(i√6η0/F0) involves the singlet only. To lowest N µ order (i.e., smallest number of derivatives) the lagrangian is built from traces of the covariants ∆µ,Γµ gVµ,∆µ∆ , µ µ µ − (Γµ gVµ)(Γ gV ), D ϑ, and arbitrary functions of the variable X(x), all being invariant under U(3)L U(3)R transformations.− − Explicit expressions of various lagrangian parts are listed for convenience in the Appendix. Their derivation can be found in Ref. [12]. Although the lagrangian, Eqn.1 appears similar to that of Bando et al.N [7] and

2 ¯ other variants of this model [8,9], the terms Li and their symmetry breaking companions Li(i = A, V ) are different, including additional terms proportional to Tr(∆ )Tr(∆µ), Tr(∆ )Dµϑ and Tr(Γµ gV µ)Tr(Γ gV ). Clearly, µ µ − µ − µ such terms must be included in any U(3)L U(3)R symmetry based HLS theory, where the traces of ∆µ and Γµ gVµ do not vanish. It is demonstrated below that these terms are needed in order to reproduce the meson masses− from the lagrangian without imposing additionalN constraints. We add also that the matrix used to generate the symmetry breaking companions is a universal Hermitian matrix defined as,

B = ξbξ + ξ†bξ† ; b = diag(bu,bd,bs) , (6) where bi = mi/(mu + md + ms). This form guarantees that the ratio of SU(3)F symmetry breaking scales to isospin symmetry breaking scales be identical to the once of the underlying QCD. Note that, the lagrangian1 involves coefficient functions Wi(X) which are absent in the SU(3) symmetry limit. Again such functions must be included in the process of extending the SU(3)L SU(3)R based model into a general U(3)L U(3)R theory. As a final remark we note that the symmetry violating mass term in Eqn. 1 has the exact form derived by Herera-Sikl´ody et al. [16]. ¯N N µ As in Refs. [7–9], the sum a(LV + LV ) contains, amongst other contributions, a vector meson mass term VµV , µ ∼ a vector-photon conversion factor VµA and the coupling of pseudoscalar pairs to both vectors and . The latter coupling can be eliminated∼ by choosing a = 2 allowing to incorporate the strict VMD lagrangian [7,26]. We shall attempt to determine the value of this parameter from data in order to assess the validity of the traditional VMD model. Altogether, the lagrangian Eqn. 1 (see Appendix), involves symmetry breaking scales ci,di (i=A,V), a general coupling constant g and a free parameter a. The symmetry breaking scales cA,dA and r are determined from global fit analysis of radiative decay widths of pseudoscalar and vector mesons [12]. The other symmetry breaking 2 scales cV ,dV and ag are determined through the diagonalization of the vector mass matrix.

III. VECTOR MESON MASSES AND MIXING

¯ We now turn to consider mixing and masses of vector mesons. The vector part a(LV + LV ) of the lagrangian 62 (See Appendix) involves a vector meson mass term V V µ which determines as already indicated the vector mixing ∼ µ scheme, uniquely. Generally speaking, the vector nonet matrix V is the sum of an octet matrix V8, involving the 0 0 ¯ 0 ρ ,ρ±,K∗ , K∗ ,K∗± and ω8 meson fields, and an SU(3)F singlet ω0 matrix, e.g., 1 V = V8 + rV ω0 1 . (7) √3 ×

Here a departure of the parameter rV from one accounts for possible nonet symmetry breaking. However, unlike the case of pseudoscalar mesons, the light vector meson spectrum shows no clear evidence for such a departure. In the discussion to follow we take rV = 1, assuming exact nonet symmetry. Then, the vector meson nonet matrix has the form (for brevity the Lorentz indices of the vector fields are omitted),

0 ρ + ω8 + ω0 ρ+ K + √2 √6 √3 ∗ 0 ρ ω8 ω0 0 V = ρ− + + K∗ . (8)  − √2 √6 √3  ¯ 0 2ω8 ω0 K∗− K∗ +  − √6 √3    As was already indicated in the introduction, the kinetic and mass terms of the lagrangian can be reduced to a standard quadratic form by first diagonalizing the kinetic matrix via a unitary transformation Υ, rescale the fields through a transformation R to restore the standard normalization of the kinetic term, and finally diagonalizing the resulting mass matrix through another unitary transformation Ω. These three steps are sufficient to define the state mixing scheme [1]. From Eqns. 8, 56 and 62 the vector meson kinetic and mass terms amounts to, 1 1 L = (∂ V ∂ V )(∂µV ν ∂ν V µ)+ (V 2 V µ) , (9) km −4 µ ν − ν µ − 2 µMV with, 1 V 2 V µ = aF 2g2 [Tr(V V µ)+w ˜ Tr(V )Tr(V µ)+ 2 µMV 8 µ 4 µ µ µ µ µ cV (Tr(BVµV )+w ˜4Tr(BVµ)Tr(V )) + dV (Tr(BVµBV )+w ˜4Tr(BVµ)Tr(BV ))] . (10)

3 Thus the kinetic term has the standard normalized quadratic form, and following [1] the mixing scheme is determined through the diagonalization of the mass matrix only and involves one mixing angle. The nonvanishing elements of the vector meson mass matrix read, c d µ2(ρ0ρ0)=µ2 1+ V (b + b )+ V [b2 + b2 +˜w (b b )2] , (11a) V 2 u d 2 u d 4 u − d  c  µ2(ρ+ρ+)=µ2 1+ V (b + b )+d b b , (11b) V 2 u d V u d  c d µ2(ω ω )=µ2 1+ V (b + b +4b )+ V [(b2 + b2 +4b2)+w ˜ (b + b 2b )2] , (11c) 8 8 V 6 u d s 6 u d s 4 u d − s   c d µ2(ω ω )=µ2 1+3˜w + V (b + b + b )(1 + 3w ˜ )+ V [b2 + b2 + b2 +˜w (b + b + b )2] , (11d) 0 0 V 4 3 u d s 4 3 u d s 4 u d s   2 0 2 1 µ (ρ ω8)=µ (bu bd)(cV + dV [bu + bd +˜w4(bu + bd 2bs)]) , (11e) V 2√3 − − 1 3 µ2(ρ0ω )=µ2 (b b ) c (1 + w˜ )+d [b + b +˜w (b + b + b )] , (11f) 0 V 6 u − d V 2 4 V u d 4 u d s r   2 2 1 3 2 2 2 µ (ω8ω0)=µ cV (bu + bd 2bs)(1 + w˜4)+dV [b + b 2b +˜w4(bu + bd 2bs)(bu + bd + bs)] , (11g) V 3√2 − 2 u d − s −   c d µ2(K K )=µ2 1+ V (b + b )+ V b b , (11h) ∗± ± V 2 u s 2 u s   c d µ2(K 0K 0)=µ2(K¯ 0K¯ 0)=µ2 1+ V (b + b )+ V b b , (11i) ∗ ∗ V 2 d s 2 d s   with,

2 2 2 µV =2ag F8 , (12)

We stress here that even in the limits of exact chiral symmetry and nonet symmetry, i.e. cV = dV =0,rV =1, the singlet vector meson mass, mω0 = µV √1+3˜w4, differs from that of the octet mω8 = µV . This distinct property ˜ µ µ occurs due to the term W4(X)Tr(Γµ gVµ)Tr(Γ gV ) in the lagrangian Eqn. 1. Only forw ˜4 = 0 (corresponding to lagrangians previously used) all nonet− meson masses− are equal. 1. Clearly, the mass matrix is nondiagonal due to the ω ω , ρ0 ω and ρ0 ω mixing. Notice though that the ρ0 ω and ρ0 ω admixtures are extremely small 8 − 0 − 8 − 0 − 8 − 0 and and vanish in the limit of exact isospin symmetry (bu = bd). First we diagonalize the mass matrix in the limit of exact isospin symmetry. In this limit, taking bs =1and neglecting terms of the order of cV bi,dV bi(i = u, d), the matrix mass assumes a particularly simple form with elements,

2 2 µ (ρρ)=µV , (13a) 2c 2d µ2(ω ω )=µ2 1+ V + V (1 +w ˜ ) , (13b) 8 8 V 3 3 4   c d µ2(ω ω )=µ2 (1 + V )(1 + 3w ˜ )+ V (1 +w ˜ ) , (13c) 0 0 V 3 4 3 4   √2 3 µ2(ω ω )= µ2 c (1 + w˜ )+d (1 +w ˜ ) , (13d) 8 0 − V 3 V 2 4 V 4 c   µ2(K )=µ2 1+ V , (13e) ∗ V 2   which we diagonalize using the unitary transformation, ω cos θ sin θ ω 8 = V V . (14) ω0 sin θV cos θV φ    −  

1 In case of broken nonet symmetry it is still possible to maintain mω0 = µV but with the nonet symmetry breaking parameter 2 6 fixed arbitrary to be rV =1/(1 + 3w ˜4)

4 It is straightforward to show that the resulting physical meson masses and mixing angle are, 1 m2 = µ2(ω ω )+µ2(ω ω ) (µ2(ω ω ) µ2(ω ω ))2 +4(µ2(ω ω ))2 , (15) ω 2 8 8 0 0 − 0 0 − 8 8 0 8 1  p  m2 = µ2(ω ω )+µ2(ω ω )+ (µ2(ω ω ) µ2(ω ω ))2 +4(µ2(ω ω ))2 , (16) φ 2 8 8 0 0 0 0 − 8 8 0 8   2µ2(ω ω ) p (m2 m2 )2 tan 2θ = 0 8 = φ − ω 1 . (17) V −(µ2(ω ω ) µ2(ω ω ) − (m2 + m2 2µ2(ω ω )2 − 0 0 − 8 8 s φ ω − 8 8 From Eqns. 13a-13e, withw ˜ = 0 the mixing angle reduces to the ideal mixing angle value θ = 54.7o,and 4 Videal − more importantly does not depend on the symmetry breaking scales cV and dV . Correspondingly, the ω, K and φ meson masses become,

2 2 2 2 1 2 2 mω = mρ ,mK = mρ(1 + 2 cV ) ,mφ = mρ(1 + cV + dV ) . (18)

2 Clearly, assuming now dV = cV /4 leads to the BKY [7] mass relations,

m2 m2 1 K = φ =1+ c . (19) 2 2 V mρ mK 2

Otherwise, taking dV = 0 leads to the mass relations of Bramon et al. [8]

m2 1 m2 K =1+ c , φ =1+c . (20) 2 V 2 V mρ 2 mρ

It should be stress though that in all of these private cases, the ρ and ω masses are equal and θV = θVideal.To describe the actual vector meson masses we should keep d as a free parameter andw ˜ = 0. Based on the proximity V 4 6 of the ρ and ω masses, the parameterw ˜4 should be very small, and to first order the ω and φ meson masses are,

2 2 2 2 mω = mρ(1 + 2w ˜4) ,mφ = mρ(1 + cV + dV )(1 +w ˜4) . (21)

2 2 2 2 Thus, the parameterw ˜4 determines the ρ ω mass splitting while (cV + dV )=2(mφ mρ)/(mρ + mω). At this − 2 − order, using the expressions above and Eq.18 for mK, along with the experimental meson mass values [2] one obtains, cV =0.69 0.01,dV 0.03 andw ˜4 0.01. Indeed,w ˜4 is very small but it plays an important role in reproducing ± ≈ ≈ 2 the ρ ω mass splitting. Note also that the actual value of the parameter dV is far smaller than cV /4=0.12 favoring the Bramon− et al. [8] relations Eqn.20. More accurate values of these parameters can be deduced from the exact expressions for the vector meson masses. These are,w ˜4 =0.009 0.0001 and dV =0.025 0.01. Correspondingly, ±o ± the mixing angle for the vector meson states is θV = (52 1) , rather close to ideal mixing. In turn, this value constrains the field admixtures due to nonideal mixing− to be± very small. With nonideal mixing taken into account, the matrix of the physical vector fields reads,

0 ρ +ω+φ ρ+ K + √2 ∗ ρ0+ω+φ 0 V = ρ − K . (22)  − √2 ∗  K K¯ 0 φ ω  ∗− ∗   −  where  =0.047 0.001, being a measure of a non-strange (strange) admixture in the physical φ and ω) fields. We ± stress again that this parameter like the symmetry breaking scales cV ,dV andω ˜4 are all determined through the diagonalization of the mass matrix and the requirement that the eigenvalues of this matrix should reproduce the physical masses. We now consider the effects of isospin symmetry breaking on the mass splitting and mixing angle. The nondiagonal part of the mass matrix now reads

2 0 0 2 0 2 0 µ (ρ ρ ) µ (ρ ω8) µ (ρ ω0) 2 0 2 0 2 2 (ρ ,ω8,ω0)= µ (ρ ω8) µ (ω8ω8) µ (ω8ω0) (23) M  2 0 2 2  µ (ρ ω0) µ (ω8ω0) µ (ω0ω0)   It is convenient to write the transformation Ω which diagonalize this submatrix in the form

5 0 Ω=Ω(ω8,ω0)Ω(ρ ,ω,φ) (24) where the transformation Ω(ω8,ω0) diagonalize the ω8,ω0 submatrix

10 0 Ω(ω8,ω0)= 0cosθV sin θV (25)  0 sin θ cos θ  − V V   and the transformation Ω(ρ0,ω,φ) diagonalize the reduced mass matrix

2 2 2 mρ mρω mρφ 1 2 0 2 2 Ω− (ω8,ω0) (ρ ,ω8,ω0)Ω(ω8,ω0)= mρω mω 0 , (26) M  2 2  mρφ 0 mφ   where,

2 1 2 2 1 2 m = m (bu bd)cV ,m= m (bu bd)cV w˜4 . (27) ρω 2 ρ − ρφ 2√2 ρ − From the Eq.26it is clear that effects of isospin symmetry breaking on the ρ and ω masses are negligibly small 0.1% 2 ≈ Takingw ˜4 =0.009 as above, it is rather easy to show that the matrix element mρφ is negligibly small and consequently the ρ0 φ mixing angle is practically zero. Therefore we may write the Ω(ρ0,ω,φ) transformation in the simple form −

cos θρω sin θρω 0 0 Ω(ρ ,ω,φ)= sin θρω cos θρω 0 (28)  − 001   where the ρ ω mixing angle is given by the expression, − m2 tan 2θ = 2 ρω . (29) ρω − m2 m2 ρ − ω With this in mind, we may now write the vector meson matrix in the form,

0 X1ρ +X2ω+φ ρ+ K + √2 ∗ 0 X2ρ +X1ω+φ 0 V = ρ − K , (30)  − √2 ∗  K K¯ 0 φ + (ρ0 sin θ ω cos θ )  ∗− ∗ ρω ρω   −  with,

X =cosθ sin θ , (31) 1 ρω − ρω X2 =cosθρω +sinθρω. (32)

To evaluate the ρ0 ω mixing angle, we need an estimate of the quark masses ratios b , (i = u, d, s). We recall that, − i

2 1 2 1 mK+ = mπ(bu + bs)(1 + cA) , (33) 2(bu + bd) 2

2 1 2 1 mK0 = mπ(bd + bs)(1 + cA) . (34) 2(bu + bd) 2 Then, we may estimate (b b ) from the masses [2] according to, u − d 2 2 b b m + m 0 u − d =2 K − K , (35) 1 (b + b )/2 m2 + m2 − u d K+ K0 which yields,

b b = 0.0160 0.0005 . (36) u − d − ±

6 With this value of b b the calculated ρ0 ω mixing angle reads, u − d − θ =( 8.9 0.3)o . (37) ρω − ±

For the parameter bs we have [12],

m 2m 2b 1 m 2 s = s = s =2 1+ c K 1 , (38) m m + m 1 b 2 A m − u d − s   π  which yields for c =0.64 0.06 (Alternative 1 of Ref. [12]) A ± b =0.940 0.01 , (39) s ± and consequently b =0.022 0.005,b=0.038 0.005. For the Alternative II (c =0.2 0.05) an QFB mixing u ± d ± A ± scheme (cA =1.4 0.1) we obtain bs =0.928,bu =0.028,bd =0.044 and bs =0.953,bu =0.0155,bd =0.0315 respectively. This± complete defining the mixing scheme in the case of isospin symmetry breaking. Here as in the analysis of the pseudoscalar meson sector [12], the exact form of the matrix B used to generate flavor symmetry breaking (see Appendix), guarantees that the ratios of SU(3) symmetry breaking scales to isospin symmetry breaking scales are identical to the ones of the underlying QCD.

IV. VECTOR MESON DOMINANCE MODEL AND PSEUDOSCALAR MESON FORM FACTORS

A particularly interesting aspect of the theory presented above is the interaction between the photon and hadronic matter. It was already illustrated by Bando et al. [7] and O’Connell et al. [26] that the HLS lagrangian is closely related to the well known VMD type 2 lagrangian and that by suitable transformation it can be reduced into the VMD type 1 lagrangian. (See also Ref. [9] on this matter). The lagrangian, Eqn.1, contains terms proportional µ µ µ to VµV and AµA and VµA , corresponding respectively to a vector meson and photon mass terms and photon-vector∼ meson interactions.∼ likewise, it contains the coupling of both vector mesons and photons to pairs of pseudoscalar mesons. By following arguments similar to those of Ref. [26] we derive in this section explicit expressions for the photon, vector, and pseudoscalar parts of the lagrangian, Eqn.1. We recall first that the covariants Γµ and ∆µ, Eqns.2,3 can be expanded as, i ie Γ = eA Q + [ ,∂ ] A 2,Q Q + ... , (40) µ µ 2 µ 2 µ − 4Fπ P P − 2Fπ {P }−P P 1 e  ∆µ = ∂µ + i Aµ[Q, ]+... , (41) −Fπ√2 P Fπ√2 P where Aµ stands for the photon field, Q = diag(2/3, 1/3, 1/3) the quark charge operator, and is the pseudoscalar nonet matrix expressed in terms of the physical fields− (see Ref.− [12]). The ellipsis ”...” denotes termsP involving higher 2 powers of the fields. The covariant Γµ gVµ contains then terms proportional to (eAµQ gVµ). By substituting − ¯ − the expression, Eqn.40 into the lagrangian part LV + LV we obtain the following photon and vector meson terms : (i) a direct photon-vector meson coupling,

aL(V,γ)=aF 2geAµ [Tr( V ,Q (1 + c B)) + 8 { µ } V c Tr(V )Tr(BQ)+d Tr( BV ,BQ )+2d Tr(BQ)Tr(bV )] , (42) V µ V { µ } V µ (ii) a photon mass term,

F 2e2 1 aL = a 8 A2 2+ (c +2d ) , (43) mγ 3 µ 3 V V   (iii) a photon-vector meson-pseudoscalar meson interactions,

2 0 e 0 e Strictly speaking, Γµ gVµ involves the terms e(Aµ tan θW Zµ)Q 2sinϑ cos ϑ TzZµ Wµ gVµ which give − − − − W W − 2√2sinϑW − rise to direct V A, V Z0 and V W couplings. − − −

7 a ie aL(V,γ, )= z2 + A Tr( Q, [ ,∂µ ] )+c Tr(B Q, [ ,∂µ ] )+ P − 2 2 π µ { { P P } V { P P } d (Tr( BQ,B[ ,∂µ ] )+2Tr(BQ)Tr(B[ ,∂µ ])) + ... , , (44) V { P P } P P } and (iv) a vector-pseudoscalar meson interaction term,

a gz2 aL(V, )= π Tr( V , [ ,∂µ ] )+c [Tr(B Q, [ ,∂µ ] )+ P − 2 2 { { µ P P } V { P P } Tr(V )Tr(B[ ,∂µ ])] µ P P d [Tr( BV ,B[ ,∂µcalP] )+2Tr(BV )Tr(B[ ,∂µ ])] + ... . (45) V { µ P } µ P P } ¯ Similarly, substituting the expansion of ∆µ, Eqn.41, into the pseudoscalar lagrangian part, LA + LA, leads to the photon-pseudoscalar meson interactions, ie L(A, γ, )= z2 A Tr( [Q, ],∂µ )+c Tr(B [Q, ],∂µ ) + ... . . (46) P − 2 π µ { { P P} A { P P} } Here the ellipsis ... denote terms involving ∂ n,n 2. To summarize the photon-meson part corresponding to the lagrangian, Eqn.1 is given by, PP ≥ 1 L = Tr(V V µν )+aL(V,γ)+aL + aL(V,γ, )+aL(V, )+L(A, γ, ) . (47) VMDII −4 µν mγ P P P

This last expression, with its photon mass Lmγ term is rather similar to the popular VMDII model [26]. By suitable transformation of the fields this term can be removed but prior to doing that some comments are in order. First, in complete analogy with the discussion of Ref. [26], we state here that the dressed photon propagator has the proper behavior of i/q2 at small photon momenta but significantly modified away from q2 = 0. This certainly may be important in− describing processes at sufficiently high momenta. Secondly, it is clear that the neutral ρ0,ωand φ meson mix with the photon, spontaneously breaking symmetry down to the U(1)em. We may exploit the diagonalization procedure described in section III to introduce physical photon and mesons. In terms of these physical fields, the lagrangian has no explicit coupling between the photon and neutral mesons, whilst there is a direct coupling between the photon to hadronic currents in the limit of exact U(3)F . Thirdly, so far the value of the parameter a is left free and in order to assign its value one needs considering the dynamics underlying QCD. Taking a = 2 at the instance of exact symmetry, i.e.cA = dA = cV = dV = 0, leads successfully to the phenomenological KSFR relations [27,28] 2 2 2 2 mρ =4F8 g ,gρ =2gρππFπ with gρ and gρππ are the ργ and ρππ effective coupling constants, respectively. Furthermore, with a = 2 the symmetric lagrangian part L(A, γ, )+aL(V,γ, ) = 0 so that the effective γππ coupling constant vanishes. Thus in the limit of exact symmetry the pionP form factor isP saturated by the ρ meson contribution. All departures from the KSFR relations, nonvanishing effective γππ coupling gγππ =0,andtheρ-meson dominance for the form factor are due to contributions from the symmetry breaking companions.6 The value a =2.13 0.06 deduced from the ρ meson mass and its decay width into two pions, indicates that VMD provides, a good± first approximation for the electromagnetic part of our lagrangian. With this value of a our lagrangian includes in addition to the usual VMD corrections due to direct photon-pseudoscalar meson coupling. We now turn to derive a more elegant form of the lagrangian Eqn.47 with manifestly a massless photon. We do that via transforming the vector meson fields into [26], 3 e = V A Q, (48) Vµ µ − g µ ˆ ˆ and define Γµ gVµ = Γµ g µ where Γµ does not involve the pure photon field contributions. In terms of − − V ¯ these new fields, the lagrangian part LV + LV neither contains a photon mass term nor an explicit coupling of the photon to pseudoscalar currents, jµ [∂µ , ]. Thus the terms Lmγ and aL(V,γ, ) are removed from LV + ¯ ∼ P P P LV . An explicit photon coupling to pseudoscalar currents, is now confined to the L(A, γ, ) term resulting from ¯ P ˆ the pseudoscalar LA + LA part which is left unaffected by the transformation Eqn.48. Clearly, since the Γµ as 2 2 ¯ before involves terms like ieAµ( Q ,Q /2)/2F8 etc. the lagrangian LV + LV still contains the photon vector meson pseudoscalar− mesonP P−{Pcouplings. Correspondingly,} since − −

3 0 Note that the same way the momentum independent direct Vµ Zµ and Vµ Wµ couplings may be transformed into momentum dependent couplings − −

8 e V = A Q, (49) µν Vµν − g µν the vector meson kinetic term transforms into, 1 1 e 2 e2 Tr(V V µν )= Tr( µν ) 2 A Tr(Q µν )+ A Aµν , (50) −4 µν −4 Vµν V − g µν V 3 g2 µν   where as usual Aµν = ∂µAν ∂ν Aµ. The first term on rhs of Eqn. 50 is the common kinetic term of the vector fields while the second term describes− the direct photon-vector meson coupling. The last term is associated with the V µν kinetic term of the photon Aµν A /4 and can be removed by rescalling the charge e and the photon field according to [26], −

2 2 2 2 Aµν0 = 1+2e /3g Aµν ; e0 = e/ 1+2e /3g , (51) Inserting all of these modifications intop the lagrangian Eqn.47 leadsp to, 1 e L = Tr( µν ) 2 A Tr(Q µν ) + VMDI −4 Vµν V − g µν V   aF 2g2 [Tr( )+w ˜ Tr( )Tr( )+c (Tr(B )+w ˜ Tr(B )Tr( )) + d (Tr(B B )+w ˜ Tr(B )Tr(B ))] 8 VV 4 V V V VV 4 V V V V V 4 V V a gz2 π Tr( , [ ,∂µ ] )+c [Tr(B Q, [ ,∂µ ] )+Tr( )Tr(B[ ,∂µ ])] − 2 2 { {Vµ P P } V { P P } Vµ P P d [Tr( B ,B[ ,∂µ ] )+2Tr(B )Tr(B[ ,∂µ ])] V { Vµ P P } Vµ P P } ie z2 A Tr( [Q, ],∂µ )+c Tr(B [Q, ],∂µ ) , (52) − 2 π µ { { P P} A { P P} } where for convenience we have omitted the 0 from Aµ0 and e0. The last expression is closely related to the VMDI model with the characteristics of (i) a massless photon, (ii) a momentum dependent direct photon-vector meson coupling µν µν A , and (iii) a direct photon-pseudoscalar meson coupling as well as the coupling through direct photon-vector Vmeson transitions. To conclude this section we now use the lagrangian LVMDI, Eq.52, to calculate the pion and kaon form factors. In the three diagram approximation (see Fig. 1) one obtains, a 1 F (q2)=1+q2 , (53a) π 2 m2 q2 ρ − 2 2 3a 2 1 1 cV 1 2(1 + cV + dV ) 1 FK (q )=1+q zK + − + , (53b) ± 2 m2 q2 3 m2 q2 3 m2 q2 " ρ − ω − φ − #

2 2 3a 2 1 1 cV 1 2(1 + cV + dV ) 1 F 0 (q )=q z − , (53c) K 2 K m2 q2 − 3 m2 q2 − 3 m2 q2 " ρ − ω − φ − # where zK =1/ 1+cA/2 is a rescalling parameter defined from applying the diagonalization procedure to the pseudoscalar meson sector. It must be stressed that the three diagrams involving V A transition vertex contribute to those terms proportionalp to Q2, only. The meson charge radii are completely determined− by diagrams involving the conversion of ρ, ω, φ mesons into a photon (diagram 1b). We have 3a 1 = , (54a) π 2 2 mρ 3a 1 1 c 1 2(1 + c + d ) 1 = z2 + V + V V , (54b) K K 2 − 2 2 ± 2 "mρ 3 mω 3 mφ # 3a 1 1 c 1 2(1 + c + d ) 1 = z2 + V + V V . (54c) K0 K 2 − 2 2 2 "−mρ 3 mω 3 mφ # The calculated radii of the pions and kaons along with data are listed in Table I. In the first three columns we list the results from our model for the Alternative I (c =0.64 0.06), Alternative II (c =0.20 0.05) and QFB scheme A ± A ± (cA =1.4 0.1) [12]. The results from the Alternative I and Alternative II explain equally well the pion and kaon radii although± the χ2 criteria is slightly better for the former. The results from the QFB scheme are by far inferior. For comparison we show the radii calculated with a =2.4 as determined by Benayoun and O’Connell [9]. Clearly, this value can not explain the charged pion radius.

9 V. SUMMARY AND CONCLUSIONS

We have considered the vector meson sector within the framework of an effective field theory based on local U(3)L U(3)R symmetry, where the vector mesons play the role of the dynamical gauge bosons of the HLS nonlinear chiral lagrangian. This supplement a previous study of the pseudoscalar sector [12] within this same framework. The lagrangianN is written in the most general way, including additional terms quadratic in the traces of the field matrices, which no longer vanish once singlet contributions are included. With these terms included, particularly ˜ the contribution from the small W4(X)Tr(Γµ gVµ)Tr(Γµ gVµ) term and its symmetry breaking companion, it is possible to reproduce the ω ρ mass splitting.− The explicit form− of the lagrangian, through its kinetic and mass terms corresponding to vector mesons,− determine the vector meson mixing scheme unambiguously. This observation which has been overlooked in previous studies allows us to calculate the model symmetry breaking scales and the value of ag2 through the diagonalization of the mass matrix. Once the value of g and a are disentangled, say from using another data point like the decay width of ρ ππ, the lagrangian is completely defined, and can be used to calculate all observables without any additional assumption.→ As an example we have calculated the form factors and radii of the pions and kaons. Our calculations reproduce straightforwardly the proper normalization of the form factors and explain nicely the charged radii. We finally indicate that VMD model emerges naturally as a first approximation of our lagrangian. This stands in marked difference with Ref. [9,18] where VMD seems to be badly violated.

Acknowledgment This work was supported in part by the Israel Ministry of Absorption.

VI. APPENDIX

For convenience we quote here the axial and vector parts of the lagrangian 1, explicitly. First the symmetric parts are,

µ µ LA = W1(X)Tr(∆µ∆ )+W4(X)Tr(∆µ)Tr(∆ )+ µ µ W5(X)Tr(∆µ)D ϑ + W6(X)DµϑD ϑ, (55) L = W˜ (X)Tr([Γ gV ][Γµ gV µ]) + V 1 µ − µ − W˜ (X)Tr(Γ gV )Tr(Γµ gV µ) , (56) 4 µ − µ − where

D ϑ = ∂ ϑ + Tr(r l ) , (57) µ µ µ − µ V = ∂ V ∂ V ig[V ,V ] . (58) µν µ ν − ν µ − µ ν ¯ The symmetry breaking companions LA are constructed in two alternative ways. The first (referred to as Alternative I) breaks SU(3)F symmetry, and corresponds to a quadratic form of the Goldstone meson kinetic energy term,

¯ ¯ ¯ µ ¯ ¯ µ LA = W1(X) cATr(B∆µ∆ )+dATr(B∆µB∆ ) + ¯ ¯ µ ¯ µ W4(X)dATr(B ∆µ)Tr(B∆ )+W5(X)cATr(B∆µ)D ϑ. (59) where the axial vector covariant, i 1 ∆¯ = ξ†,∂ ξ + ξ†r ξ ξ l ξ† , (60) µ 2 8 µ 8 2 8 µ 8 − 8 µ 8 n o   involves the pseudoscalar octet field matrix only. The second alternative (Alternative II) breaks U(3)F symmetry and corresponds to bilinear kinetic energy term. Namely,

¯ µ µ LA = W1(X)(cATr(B∆µ∆ )+dATr(B∆µB∆ )) µ µ +W4(X)(cATr(B∆µ)Tr(∆ )+dATr(B∆µ)Tr(B∆ )) + µ W5(X)cATr(B∆µ)D ϑ, (61) with ∆µ involves now the nonet pseudoscalar field matrix. Similarly the asymmetric companion of LV is,

10 L¯ = W˜ (X)(c Tr(B[Γ gV ][Γµ gV µ]) + V 1 V µ − µ − d Tr(B[Γ gV ]B[Γµ gV µ])) + V µ − µ − W˜ (X)(c Tr(Γ gV )Tr(B[Γµ gV µ]) 4 V µ − µ − + d Tr(B[Γ gV ])Tr(B[Γµ gV µ])) . (62) V µ − µ −

In the expressions above, cA,dA,cV and dV are the symmetry breaking parameters determined from data. In Ref. [12] cA,dA and r = F8/F0 were deduced from global fit of calculated radiative decay widths. Their values as determined for the Alternative I , Alternative II and QFB mixing schemes are listed in Table II.

TABLE I. Pseudoscalar meson squared charge radii in fm2 Alternative I Alternative II QFB g =2.4Data 2 0.43 0.01 0.43 0.01 0.43 0.01 0.49 0.01 0.439 0.03 2 ± ± ± ± ± ± 0.29 0.02 0.35 0.02 0.23 0.01 0.33 0.02 0.31 0.05 2± ± ± ± ± ± 0.040 0.003 0.048 0.003 0.031 0.003 (0.045 0.003) 0.054 0.026 χ2 − 0.±50− .±73− .±32− .±1 − ±

TABLE II. Symmetry breaking scales and χ2/dof from global fit to data. Values marked with an asterisk were kept fixed.

cW cA dA r Alternative I (0.20 0.05) (0.64 0.06) 0.25 0.04 0.91 0.04 Alternative II −(0.27 ± 0.05) 0.2 ±0.05− 0.1 ±0.02 0.94 ± 0.05 QFB −(0.19 ± 0.05) (1.4± 0.1) 1.1± 0.1*1± − ± ± − ±

11 [1] E.Gedalin, A.Moalem and L.Razdolskaya, Phys. Rev. 64 D, 076007 (2001). [2] D.E.Groom et al. (Particle Data Group) , Eur. Phys. J. C15, 1 (2000) . [3] T.Feldmann et al., Phys. Rev.D58, 114006 (1998). [4] T.Feldmann et al., Phys. Lett., B449, 339 (1999); hep-ph/9812269. [5] T.Feldmann, Nucl. Phys. Suppl. 74, 1512 (1999); hep-ph/9807367. [6] R.Escribano, J.-M.Fr`ere, Phys.Lett., B459, 288 (1999). [7] M.Bando, T.Kugo and K.Yamawaki, Nucl.Phys.B259, 493 (1985); Prog. Theor. Phys. 73, 1541 (1985); Phys.Rep.164, 218 (1988). [8] A.Bramon, A.Grau and G.Pancheri,The second DAΦNE physics handbook,v.II, p.477, Eds. L.Maiani et al.,INFN, Frascati, May 1995. [9] M.Benayoun and H.B.O’Connell, Phys. Rev. D58, 074006 (1998). [10] M.Benayoun et al. Eur. Phys. J. C2,269 (1998). [11] E.Gedalin, A.Moalem and L.Razdolskaya, nucl-th/0006073. [12] E.Gedalin, A.Moalem and L.Razdolskaya, hep-ph/0108124. [13] J.Gasser, H.Leutwyler, Nucl.Phys. B250, 465 (1985). [14] H.Leutwyler, Phys.Lett.B374, 163 (1996). [15] H.Leutwyler, Proc.QCD 97, Montpellier, France, July 1997,Ed.S.Narison,Nucl.Phys. B(Proc. Suppl.) 64,223 (1998); hep- ph/9709408. [16] P.Herera-Sikl´ody et al., Nucl.Phys. B497, 345 (1997) . [17] Jun Gao and Bing An Li, Phys. Rev. D61, 113006 (2000). [18] M.Harada and K.Yamawaki, Phys. Rev. Lett. 87, 152001 (2001). [19] E.Witten, Nucl. Phys.B223, 422 (1983). [20] C.G.Callan, E Witten, Nucl.Phys.,B239, 161 (1984). [21] J. Gasser, M. E. Sainio and A.Svarc,ˆ Nucl. Phys. B 307, 779 (1988). [22] J.Bijnens, A.Bramon and F.Cornet, Z.f¨ur Phys.C46, 595 (1990). [23] J.Bijnens, A.Bramon and F.Cornet, Phys. Lett. B237, 488 (1990). [24] A.Krause, Helv. Phys. Acta, 63, 3 (1990). [25] V. Bernard, N. Kaiser and Ulf-G.Meissner, Int. J. Mod. Phys. E4, 193 (1995). [26] H.B.O’Connell et al., Prog. Part. Nucl. Phys. 39, 201 (1997). [27] K.Kawarabayashi and M.Suzuki, Phys. Rev. Lett. 16, 255 (1966). [28] and Fayazuddin, Phys. Rev. 147, 1071 (1966).

~ ~ ~~P ~~P ~~ ~~ ~~~ ~~~ ~~ ~~~~~ ~~ ~~ ~~P ~~P a b

FIG. 1. Diagram contributing into pseudoscalar meson form factors. Diagram a corresponds to direct photon-pseudoscalar meson coupling, diagram b involves the photon-vector meson conversion.

12