arxiv: v2 [hep-ph] 24 Apr 2008

Transcription

1 Octet baryon magnetic moments from QCD sum rules Lai Wang and Frank X. Lee Physics Department, he George Washington University, Washington, DC 005, USA arxiv: v [hep-ph] 4 Apr 008 A comprehensive study is made for the magnetic moments of octet baryons in the method of QCD sum rules. A complete set of QCD sum rules is derived using the external field method and generalized interpolating fields. For each member, three sum rules are constructed from three independent tensor structures. hey are analyzed in conjunction with the corresponding mass sum rules. he performance of each of the sum rules is examined using the criteria of OPE convergence and ground-state dominance, along with the role of the transitions in intermediate states. Individual contributions from the u, d and s quarks are isolated and their implications in the underlying dynamics are explored. Valid sum rules are identified and their predictions are obtained. he results are compared with experiment and previous calculations. PA numbers: 3.40.Em,.38.-t,.38.Lg,.55.Hx, 4.0.Gk, 4.0.Jn I. INRODUCION he QCD sum rule method is a nonperturbative analytic formalism firmly entrenched in QCD with minimal modeling. he field remains active judging by the 3000 and growing references to the seminal paper of Shifman, Vainshtein and Zakharov[] that introduced the method. he approach provides a general way of linking hadron phenomenology with the interactions of quarks and gluons via only a few parameters: the QCD vacuum condensates and susceptibilities. he studies give an unique perspective on how the properties of hadrons arise from nonperturbative interactions in the QCD vacuum and how QCD works in this context. It has been successfully applied in almost every aspect of strong-interaction physics. Calculations of the magnetic moment were carried out soon after the method was introduced for the proton, neutron [, 3] and hyperons [4] in the external field method. In this method, a static magnetic field is introduced that couples to the quarks and polarizes the QCD vacuum. Magnetic moments can be extracted from the linearresponsetothis field. he resultsofthe studies validated the external field method as a way of probing hadron properties other than the mass, such as magnetic moments, form factors, axial charge, isospin-breakings. Later, a more systematic study was made for the magnetic moments of octet baryons [5, 6, 7, 8]. Calculations were also carried out for decuplet baryons [9, 0,, ] and the rho meson [3]. here are other studies of magnetic moments using the light-cone QCD sum rule method [4, 5, 6, 7] which will not discuss here. In this work, we carry out a comprehensive, independent calculation of the magnetic moments of the octet baryons in the external field method. It can be considered as an update over the previous calculations [5, 6, 7] which were done more than 0 years ago. here are a number of things we do differently. First, we employ generalized interpolating fields which allow us to use the optimal mixing of interpolating fields to achieve the best match. Second, we derive a new, complete set of QCD sum rules at all three tensor structures and analyze all of them. he previous sum rules, which were mostly limited to one of the tensor structures, correspond to a special case of the mixing in our sum rules. In this way, we provide an independent check of the previous sum rules. hird, we perform a Monte-Carlo analysis which has become standard nowadays. he advantage of such an analysis is explained later. Fourth, we use a different procedure to extract the magnetic moments and to treat the transition terms in the intermediate states. Our results show that these transitions cannot be simply ignored. Fifth, we isolate the individual quark contributions to the magnetic moments and discuss their implications in the underlying quark-gluon dynamics in the baryons. he paper is organized as follows. In Section II, the method of QCD sum rules is introduced. Using the interpolating fields for the octet baryons, master formula are calculated. hen both the phenomenological representation and QCD side are derived. Section III will list the sum rules we derived for the octet baryon family, followed by the analysis to extract the magnetic moments in Section IV. Section V summarizes the results and gives a comparison of our results with experiment and previous calculation, followed by an in-depth discussion of our findings. Our conclusions are given in Section VI. II. MEHOD he starting point is the time-ordered correlation function in the QCD vacuum in the presence of a constant background electromagnetic field F µν : Π(p) = i d 4 x e ip x 0 { η(x) η(0) } 0 F. () he QCD sum rule approach is to evaluate this correlation at two different levels. On the quark level, it describes a hadron as quarks and gluons interacting in the QCD vacuum. On the phenomenological level, It is saturated by a tower of hadronic intermediate states

2 with the same quantum numbers. his way, a connection can be established between a description in terms of hadronic degrees of freedom and one based on the underlying quark and gluon degrees of freedom governed by QCD. Here η is the interpolating field (or hadron current) with the quantum numbers of the hadron under consideration. he subscript F means that the correlation function is to be evaluated with an electromagnetic interaction term added to the QCD Lagrangian: L I = A µ J µ, () where A µ is the external electromagnetic potential and J µ = e q qγ µ q is the quark electromagnetic current. Since the external field can be made arbitrarily small, one can expand the correlation function Π(p) = Π (0) (p)+π () (p)+, (3) where Π (0) (p) is the correlation function in the absence of the field, and gives rise to the mass sum rules of the baryons. he magnetic moments will be extracted from the QCD sum rules obtained from the linear response function Π () (p). he action of the external electromagnetic field is twofold: it couples directly to the quarks in the baryon interpolating fields, and it also polarizes the QCD vacuum. he latter can be described by introducing new parameters called vacuum susceptibilities. he interpolating field is constructed from quark fields with the quantum number of baryon under consideration and it is not unique. We consider a linear combination of the two standard local interpolating fields. hey read for the baryon octet family: η p (uud) = ǫ abc [(u a Cγ 5 d b )u c +β(u a Cd b )γ 5 u c ], η n (ddu) = ǫ abc [(d a Cγ 5 u b )d c +β(d a Cu b )γ 5 d c ], η Λ (uds) = 6 ǫabc [(u a Cγ 5 d b )s c +(u a Cγ 5 s b )d c (d a Cγ 5 s b )u c +β((u a Cd b )γ 5 s c +(u a Cs b )γ 5 d c (d a Cs b )γ 5 u c )], η Σ (dds) = ǫ abc [(d a Cγ 5 s b )d c +β(d a Cs b )γ 5 d c ], η Σ0 (uds) = ǫ abc [(u a Cγ 5 s b )d c +(d a Cγ 5 s b )u c +β((u a Cs b )γ 5 d c +(d a Cs b )γ 5 u c )], η Σ+ (uus) = ǫ abc [(u a Cγ 5 s b )u c +β(u a Cs b )γ 5 u c ], η Ξ (ssd) = ǫ abc [(s a Cγ 5 d b )s c +β(s a Cd b )γ 5 s c ], η Ξ0 (ssu) = ǫ abc [(s a Cγ 5 u b )s c +β(s a Cu b )γ 5 s c ]. (4) Here u and d are up-quark and down-quark field operators, C is the charge conjugation operator, the superscript means transpose, and ǫ abc makes it color-singlet. he normalization factors are chosen so that correlation functions of these interpolating fields coincide with each other under SU(3)-flavor symmetry. he real parameter β allows for the mixture of the two independent currents. hechoiceadvocatedbyioffe[8]andoften usedin QCD sum rules studies corresponds to β =. We will take advantage of this freedom to achieve optimal matching in the sum rule analysis. A. Phenomenological Representation We start with the structure of the two-point correlation function in the presence of the electromagnetic vertex to first order Π(p) = i d 4 xe ipx < 0 η(x)[ i d 4 ya µ (y)j µ (y)] η(0) 0 >. (5) Inserting two complete sets of physical intermediate states, we restrict our attention only to the positive energy ones and write Π(p) = d 4 xd 4 y d4 k d 4 k i i (π) 4 (π) 4 k M N N s s N iε k MN iε (6) e ipx A µ (y) < 0 η(x) N k s >< N k s J µ (y) Nks >< Nks η(0) 0 >. We can use the translation invariance to express η(x) in terms of η(0) < 0 η(x) N k s >=< 0 η(0) N k s > e ik x. (7) he interpolating field excites (or annihilates) the ground state as well as the excited states of the baryon from the QCD vacuum. he ability to do so is described by a phenomenological parameter λ N (called current coupling or pole residue), defined by the overlap for the ground state where u is the Dirac spinor. ranslation invariance on J µ (y) gives < 0 η(0) Nks >= λ N u(k,s), (8) < N k s J µ (y) Nks >= e iqy < N k s J µ (0) Nks >, (9)

3 3 where q = k k is the momentum transfer and Q = q. he matrix element of the electromagnetic current has the general form < k s J µ (0) ks >= ū(k,s )[F (Q )γ µ +F (Q µν qν )iσ M N ]u(k,s), (0) where the Dirac form factors F and F are related to the Sachs form factors by G E (Q ) = F (Q ) Q (M N) F (Q ) G M (Q ) = F (Q )+F (Q ). () At Q = 0, F (0) =, F (0) = µ a which is the anomalous magnetic moment, and G M (0) = F (0)+F (0) = µ which is the total magnetic moment. Writing out explicitly only the contribution of the ground-state nucleon and denoting the excited state contribution by ESC, we have Π(p) = λ N d 4 xd 4 y d4 k d 4 k (π) 4 (π) [k M 4 N iε] [k MN iε] A µ (y)e i(p k )x e iqy u(k,s )ū(k,s )[F (Q )γ µ +F (Q µν qν )iσ M N ] u(k,s)ū(k,s)+esc. s s () he spin sums are of the form u(k,s)ū(k,s) = ˆk +M N. (3) s QCD sum rule calculations are most conveniently done in the fixed-point gauge. For electromagnetic field, it is defined by x µ A µ (x) = 0. In this gauge, the electromagnetic potential is given by Changing variables from k to q = k k, then d 4 k = d 4 q, we have A µ (y) = F µνy ν. (4) Π(p) = λ N F µν d 4 xd 4 y d4 k d 4 q (π) 4 (π) [k M 4 N iε] [(q k ) MN iε] e i(p k )x ( i q ν e iq y )(ˆp+M N )[F (Q )γ µ +F (Q µν qν )iσ M N ](ˆk ˆq +M N )+ESC. (5) Integrating over x, we get a delta function d 4 x (π) 4ei(p k )x = δ 4 (p k ). (6) Integrating / q ν by parts, then doing d 4 y, we can get another delta function δ(q). Since we have a δ(q), when doing / q ν only terms linearin q ν contribute. We havef (Q ) q=0 =,F (Q ) q=0 = µ a = µ, and / q ν F, (Q ) q=0 = 0, so derivatives of the structure functions do not enter. Finally, we arrive at Π(p) = i λ N F µν[p M N iε] (ˆp+M N ){ i(µn ) M N σ µν (ˆp+M N )+iσ µν (p µ γ ν p ν γ µ )(ˆp+M N )[p M N iε] }+ESC. (7) Examination of its tensor structure reveals that its has 3 independent combinations: F µν (ˆpσ µν +σ µνˆp), F µν i(p µ γ ν p ν γ µ )ˆp and F µν σ µν. he momentum-space correlation function in the above equation can be written in terms of these three structures λ N Fµν Π(p) = 4 {[σ µν ](M [p MN N µ N )+ µn iε] M N (p MN )+µ N[ˆpσ µν +σ µνˆp] [i(p µ γ ν p ν γ µ )]ˆp}+ESC, + (µn ) M N (8) where we have used the following identities (ˆp+M N )σ µν = M N σ µν + (ˆpσµν +σ µνˆp)+i(p µ γ ν p ν γ µ ), (ˆp+M N )σ µν (ˆp+M N ) = σ µν (p +M N )+M N(ˆpσ µν +σ µνˆp)+i(p µ γ ν p ν γ µ )ˆp, (ˆp+M N )(p µ γ ν p ν γ µ )(ˆp+M N ) = (p M N )(pµ γ ν p ν γ µ ). (9)

4 4 FIG. : he three kinds of contributions to the spectral function in the presence of an external field: ground state, transitions between ground state and excited states, and pure excited states Next step is to perform the Borel transform defined by ˆB[f(p )] = lim p,n p /n=m Upon Borel transform the ground state takes the form n! ( p ) n+ ( d dp )n f(p ). (0) ˆB[Π(p)] = λ N 4M e M N /M { M N (MN µ N M (µ N ))[F µν σ µν ]+ µ N [F µν (ˆpσ µν (µn ) +σ µνˆp)]+ [F µν i(p µ γ ν p ν γ µ )ˆp]}, M N () where M is the Borel mass, not to be confused with the nucleon mass M N. Here we must treat the excited states with care. For a generic invariant function, the pole structure can be written as C N N (p M N ) + N C N N (p M N )(p M N ) + N C N N (p M N ), () where C N N, C N N and C N N are constants. he first term is the ground state pole which contains the desired magnetic moment µ N. he second term represents the non-diagonal transitions between the ground state and the excited states caused by the external field. he third term is pure excited state contributions. hese different contributions can be represented by the diagrams in Fig.. Upon Borel transform, it takes the form [ λ N ( µn M e M N /M +e M N /M C N N e (M M N N MN N M N)/M )] + C N N M e M N /M. (3) N he important point is that the transitions give rise to a contribution that is not exponentially suppressed relative to the ground state. his is a general feature of the external-field technique. he strength of such transitions at each structure is a priori unknown and is an additional source of contamination in the determination of the magnetic moment µ N. he standard treatment of the transitions is to approximate the quantity in the square brackets by a constant, which is to be extracted from the sum rule alongwith the ground state propertyof interest. Inclusion of such contributions is necessary for the correct extraction of the magnetic moments. he pure excited state contributions are exponentially suppressed relative to the ground state and can be modeled in the usual way by introducing a continuum model and threshold parameter. B. Calculation of the QCD Side We start by contracting out the quark pairs in Eq. () using Wick s theorem, resulting in the so-called master formula in terms of quark propagators. he master formula for the proton (with uud quark content) is Ω {η N (x) η N (0) Ω = 4ǫ abc ǫ a b c { +Su aa γ 5d Cγ5 S bb u +Saa u r( d Cγ5 Su bb γ 5) +βγ 5 Su aa γ 5 d u +βγ 5 Su aa r(u d γ 5 ) +βsu aa d Cγ5 Su bb γ 5 +βsu aa γ 5r(u Cγ5 Sd bb ) +β γ 5 Su aa d u γ 5 +β γ 5 Su aa γ 5r(d u ) }. (4)

5 he master formula for neutron (with ddu quark content) can be obtained by exchanging the d quark with a u quark from Eq. (??). By replacing the d quark with a s quark, one can get the master formula for Σ + (with uus quark content). While the master formula for Σ (with uus quark content) can be obtained by replacing the u quark with a d quark from Σ + master formula. By exchanging the u quarks and s quarks in Σ + master formula, the master formula for Ξ 0 (with ssu quark content) can be obtained. Likewise, by replacing the u quark with a d quark, one can get the master formula for Ξ + (with ssd quark content). he master formulae for Σ 0 and Λ (uds quark content) have a more complicated structure. hey can be written in a combined way as Ω { η(x) η(0) } Ω = fǫ abc ǫ a b c { f f Ss aa γ 5 u Cγ5 Sd bb f f 3 Ss aa γ 5 d Cγ5 S bb +f f S aa d γ 5 u f 3 f Su aa γ 5 d +f f 5 βs aa s u +f f 4 βsd aa u Cγ5 S bb Cγ5 S bb Cγ5 S bb s +f f 3 S aa d γ 5 s +f 3 f Su aa γ 5 s d γ 5 f f 6 βs aa s s γ 5 +f f 6 βsd aa s f 3 f 4 βsu aa d Cγ5 Ss bb γ 5 +f 3 f 5 βs aa +f 4 f βγ 5 Ss aa γ 5 u +f 5 f βγ 5 Sd aa γ 5 u f 6 f βγ 5 Su aa γ 5 d +f 4 f 5 β γ 5 Ss aa u +f 5 f 4 β γ 5 Sd aa u f 6 f 4 β γ 5 S aa u d s u s d f 4 f 3 βγ 5 Ss aa γ 5 d s +f 5 f 3 βγ 5 Sd aa γ 5s s +f 6 f βγ 5 Su aa γ 5 s d γ 5 f 4 f 6 β γ 5 Ss aa d s γ 5 +f 5 f 6 β γ 5 Sd aa s s γ 5 +f 6 f 5 β γ 5 S aa u s Cγ5 Su bb +f f Ss aa r(u Cγ5 Sd bb γ 5 ) Cγ5 Su bb +f f Sd aa r(s Cγ5 Su bb γ 5) Cγ5 Sd bb +f 3 f 3 Su aa r(s Cγ5 Sd bb γ 5 ) d Cγ5 Su bb γ 5 +f f 4 βss aa γ 5 r(d Cγ5 Su bb ) Cγ5 Su bb γ 5 +f f 5 βsd aa γ 5 r(s Cγ5 Su bb ) Cγ5 Sd bb γ 5 +f 3 f 6 βsu aa γ 5r(s Cγ5 Sd bb ) u +f 4 f βγ 5 Ss aa r(d u +f 5f βγ 5 Sd aa r( s d +f 6 f 3 βγ 5 Su aa r(s u γ 5 +f 4 f 4 β γ 5 Ss aa γ 5 r(u u γ 5 +f 5 f 5 β γ 5 Sd aa γ 5 r(s d γ 5 +4f 6 f 6 β γ 5 S aa u γ 5 r( s u γ 5 ) u γ 5) d γ 5 ) d ) u ) 5 d )}, (5) where the various factors are: for Σ 0, f = /3,f = 0,f =,f 3 =,f 4 = 0,f 5 =,f 6 = ; and for Λ, f =,f =,f =,f 3 =,f 4 =,f 5 =,f 6 =. In the above equations, S ab q (x,0;f) 0 { qa (x) q b (0) } 0 F, q = u,d,s, (6) is the fully interacting quark propagator in the presence of the electromagnetic field. o first order in F µν and m q (assume m u = m d = 0,m s 0), and order x 4, it is given by the operator product expansion (OPE) [3, 6, 7]: S ab q (x,0;z) i π ˆx x 4δab m q 4π x δab qq δab + im q 48 qq ˆxδab + 9 qg cσ Gq x δ ab im q 5 qg cσ Gq ˆxx δ ab qq g cg x 4 δ ab ( ) λ n ab + i 483π g c G ˆxσαβ +σ αβˆx ( λ n x ) ab ( ) λ 9 qg cσ Gq σ αβ n ab + im q 768 qg cσ Gq (ˆxσ αβ +σ αβˆx ) ( ) λ n ab + ie q ˆxσ αβ +σ αβˆx 3π F αβ x δ ab e q 4 χ qq F αβσ αβ δ ab + ie qm q 96 χ qq F αβ (ˆxσ αβ +σ αβˆx ) δ ab + e q 88 qq F ( αβ x σ αβ x ρ x β σ βα) δ ab + e q [ x (κ+ξ)σ αβ x ρ x β (κ ξ)σ βα] δ ab + i 3π (g cg n )ˆxσαβ +σ αβˆx αβ x + ( λ 3 0 qq g cg x σ αβ n 576 qq F αβ e q (κf 6 qq αβ i4 )( ) λ ξǫ αβµνf µν n ab +higher order terms. (7) ) ab

6 6 We use the convention ǫ 03 = + in this work. In addition to the standard vacuum condensates, the vacuum susceptibilities induced by the external field are defined by qσ µν q F e q χ qq F µν, qg c G µν q F e q κ qq F µν, qg c ǫ µνρλ G ρλ γ 5 q F ie q ξ qq F µν. (8) Note that χ has the dimension of GeV, while κ and ξ are dimensionless. With the above elements in hand, it is straightforward to evaluate the correlation function by substituting the quark propagator into the various master formulae. We keep terms to first order in the external field and in the strange quark mass. erms up to dimension 8 are considered. he algebra is extremely tedious. Each term in the master formula is a product of three copies of the quark propagator. here are hundreds of such terms over various color permutations. he calculation can be organized by diagrams (similar to Feynmann diagrams) in Fig. andfig.3. Notethat eachdiagramisonlygenericandall possible color permutations are understood. he QCD side has the same tensor structure as the phenomenological side and the results can be organized according to the same 3 independent structures. III. QCD SUM RULES Once we have both the QCD side (LHS) and the phenomenological side(rhs), we can derive the sum rules by matching both sides. Since there are three independent tensor structures, three sum rules can be constructed. We denote these tensor structures by the following shorthand notation WE = F µν (ˆpσ µν +σ µνˆp), (9) WO = F µν σ µν, (30) WO = F µν i(p µ γ ν p ν γ µ )ˆp (3) he sum rule from WE involves only dimension-even condensates, so we call this sum rule chiral-even. he sum rule from both WO and WO involves only dimension-odd condensates, so we call them chiral-odd. Note that in previous works [3, 5] the dimension of the tensor structures, rather than the dimension of the condensates, was used to refer to the sum rules. he two names are opposite. Now we are ready to collect all of the QCD sum rules. At the structure WE, all of the sum rules can be expressed in the following form c L 4/9 E (w)m 4 +c m s χal 6/7 E (w)m +c 3 χa L 4/7 E 0 (w)+c 4 bl 4/9 E 0 (w)+(c 5 +c 6 )m s al 4/9 E 0 (w) +(c 7 +c 8 )a L 4/9 M +c 9 χm 0 a L 8/7 M +c 0 m s m 0 al 6/7 M = λ N [ µn M +A]e M N /M, (3) where the coefficients differ from member to member. he quark condensate, gluon condensate, and the mixed condensate are a = (π) ūu, b = g c G, ūg c σ Gu = m 0 ūu. (33) he quark charge factors e q are given in units of electric charge e u = /3, e d = /3, e s = /3. (34) Note that we choose to keep the quark charge factors explicit in the sum rules. he advantage is that it can facilitate the study of individual quark contribution to the magnetic moments. he parameters f and φ account for the flavor-symmetry breaking of the strange quark in the condensates and susceptibilities: f = ss ūu = sg cσ Gs ūg c σ Gu, φ = χ s χ = κ s κ = ξ s ξ. (35) he anomalous dimension corrections of the interpolating fields and the various operators are taken into account in the leading logarithmic approximation via the factor [ L γ αs (µ ] γ [ ) ln(m /Λ QCD = α s (M = ) ] γ ) ln(µ /Λ QCD ), (36) where µ = 500 MeV is the renormalization scale and Λ QCD is the QCD scale parameter. As usual, the pure excited state contributions are modeled using terms on the OPE side surviving M under the assumption of duality, and are represented by the factors E n (w) = e w /M (w /M ) n, (37) n! where w is an effective continuum threshold and it is in principle different for different sum rules and we will treat it as a free parameter in the analysis. Also, λ N is the rescaled current coupling λ N λ /4. We only need to carry out four separate calculations for Σ + (uus), Σ 0 (uds), Ξ + (ssd) and Λ(uds). he QCD n

7 7 ( d ) ( d ) ( d3 ) ( d4 ) ( d5 ) ( d6 ) ( d7 ) ( d8 ) ( d9 ) ( d0 ) ( d ) ( d ) ( d3 ) ( d4 ) ( d5 ) ( d6 ) ( d7 ) ( d8 ) ( d9 ) ( d0 ) ( d ) ( d ) FIG. : Non-mass diagrams considered for the octet baryon magnetic moments. ( d ) ( d ) ( d3 ) ( d4 ) ( d5 ) ( d6 ) ( d7 ) ( d8 ) ( d9 ) ( d0 ) ( d ) ( d ) ( d3 ) ( d4 ) ( d5 ) ( d6 ) ( d7 ) ( d8 ) ( d9 ) ( d0 ) FIG. 3: Diagrams considered for the strange quark mass corrections to the octet baryon magnetic moments.

8 8 sum rules for other members can be obtained from them by making appropriate substitutions specified below. For Σ + at WE : c = 6 [(+β) c s (3+β +3β )c u ]; c = 6 (+β)[( +3β)c sf s φ (+3β)c u ]; c 3 = ( +β)[(+β)c sf s φ c u ( β +f s +βf s )]; c 4 = 384 [(3 β +3β )c s +(+6β +β )c u ]; c 5 = 4 [(3 8β +3β f s +β f s )c s +( 5+β +7β 9f s +6βf s +7β f s )c u ]; c 6 = 4 ( +β)[(+β)c sf s φ(8κ ξ)+c u (4κ+6βκ+ξ)]; c 7 = 7 ( +β)[c s( 3+3β +f s +βf s )+c u ( +β +8f s +8βf s )]; c 8 = 7 ( +β)[(+β)c sf s φ(8κ ξ) +c u ( κ+βκ+0f s κ+0βf s κ+ξ βξ +f s ξ +βf s ξ)]; c 9 = 88 [(7(+β)c sf s φ c u (3 3β +f s +βf s )]; c 0 = 96 [( 4β +β )c s +( 3 β +3β )c u (+f s )]. (38) For Ξ 0 at WE : c = 6 [(+β) c u (3+β +3β )c s ]; c = 6 (+β)[( +3β)c u (+3β)c s f s φ]; c 3 = ( +β)[(+β)c uf s c s f s φ(f s βf s ++β)]; c 4 = 384 [(3 β +3β )c u +(+6β +β )c s ]; c 5 = 4 [((3 8β +3β )f s +β )c s +(( 5+β +7β )f s 9+6β +7β )c s ]; c 6 = 4 ( +β)[(+β)c u(8κ ξ)+c s f s φ(4κ+6βκ+ξ)]; c 7 = 7 ( +β)[c uf s (( 3+3β)f s ++β)+c s f s (( +β)f s +8+8β)]; c 8 = 7 ( +β)[(+β)c uf s (8κ ξ) +c s f s φ( κf s +βκf s +0κ+0βκ+ξ βf s ξ +f s ξ +βξ)]; c 9 = 88 [(7(+β)c uf s c s f s φ((3 3β)f s ++β)]; c 0 = 96 [( 4β +β )c u +( 3 β +3β )c s (+f s )]. (39) For Λ at WE : c = 48 (( +β) (c d +c u )+(3+0β +3β )c s ); c = 6 (+β)(( +β)(c d +c u ) (+5β)c s f s φ); c 3 = 7 ( +β)((c d +c u )( β +f s +5βf s ) (5+β)c s f s φ); c 4 = 5 ((7+β +7β )(c d +c u )+(4+6β +4β )c s ); c 5 = 44 (c s( 78+β(84 6f s )+β (0 4f s )+0f s ) +( +β)(c d +c u )(7+3f s +β(4+33f s ))); c 6 = 7 ( +β)(c sf s φ(8(+β)κ+(5+β)ξ)+(c d +c u )((0+β)κ (+β)ξ)); c 7 = 6 ( +β)(c s(39+33β 0f s βf s )+(c d +c u )( β +4f s +0βf s )); c 8 = 43 ( +β)(c sf s φ(8(+β)κ+(5+β)ξ) +(c d +c u )((9+7β +f s +5βf s )κ ( β +f s +5βf s )ξ)); c 9 = 576 ( +β)((c d +c u )(5+3β +3f s +7βf s ) (5+β)c s f s φ); c 0 = 9 ((6 4β 34β )c s ( 7 0β +5β )(c d +c u )(+f s )). (40) At the structure WO, the sum rules can be expressed in the following form c m s L 8/9 E (w)m 4 +c χal 6/7 E (w)m 4 +(c 3 +c 4 )ae (w)m +c 5 m s χa L 6/7 E 0 (w) +c 6 χabl 6/7 [ E 0 (w)+(c 7 +c ] 8 )m s a M +c 9 ab M = λ N M µn N M + µn +A e MN M N /M. (4)

9 9 For Σ + at WO : c = ( +β)((+β)c s c u ); c = 4 ( +β)[( +β)c sf s φ+8(+β)c u ]; c 3 = 4 ( +β)[c s( 6 6β f s +βf s ) 6c u (+β f s +βf s )]; c 4 = 96 ( +β)[( +β)c sf s φ(4κ 3ξ) 8(+β)c u (κ+ξ)]; c 5 = [(3+β +3β )c s f s φ (3+4β +9β )c u (+f s )]; c 6 = 576 ( +β)[( +β)c sf s φ 8(+β)c u ]; c 7 = 36 [c s( 3+3β +5f s βf s +5β f s ) c u (3+4β +9β 3f s +0βf s +9β f s )]; c 8 = 44 [c sf s φ(4κ+8βκ+4β κ+ξ 0βξ +β ξ) c u ( (3+4β +9β )(+f s )ξ +4(β +6β +3f s +βf s +3β f s )κ)]; c 9 = 576 [3c s c u ( 3+3β f s +βf s )]. (4) For Ξ 0 at WO : For Λ at WO : c = ( +β)((+β)c u c s ); c = 4 ( +β)[( +β)c u +8(+β)c s f s φ]; c 3 = 4 ( +β)[c u(( 6 6β)f s +β) 6c s ((+β)f s +β)]; c 4 = 96 ( +β)[( +β)c u(4κ 3ξ) 8(+β)c s f s φ(κ+ξ)]; c 5 = [(3+β +3β )c u (3+4β +9β )c s f s φ(+f s )]; c 6 = 576 ( +β)[( +β)c u 8(+β)c s f s φ]; c 7 = 36 [c uf s (( 3+3β )f s +5 β +5β ) c s f s ((3+4β +9β )f s 3+0β +9β )]; c 8 = 44 [c uf s (4κ+8βκ+4β κ+ξ 0βξ +β ξ) c s f s φ( (3+4β +9β )(+f s )ξ +4((β +6β )f s +3+β +3β )κ)]; c 9 = 576 [3c u c s (( 3+3β )f s +β)]. c = 6 ( +β)((c d +c u )(+β) (5+β)c s ); c = 7 ( +β)((7+β)(c d +c u )+(37+35β)c s f s φ); c 3 = 7 ( +β)((c d +c u )(4+β( 3f s )+3f s )+c s ( 30 6β +35f s +37βf s )); c 4 = 88 ( +β)((c d +c u )(( 3+5β)κ+(7 35β)ξ) c s f s φ(58κ+86βκ+49ξ +3βξ)); c 5 = 4 (( +β )(c d +c u )(+f s )+(5+6β +3β )c s f s φ); c 6 = 78 (( +β)((+7β)(c d +c u )+(35+37β)c s f s φ)); c 7 = 6 ( ( +β)(c d +c u )( 7+β f s +3βf s )+c s ( 30+34f s +4β(6+7f s )+β (6+8f s ))); c 8 = 864 (c sf s φ(4(7+0β +9β )κ (+6β +35β )ξ) +( +β)(c d +c u )(4( 4 β +5f s +7βf s )κ ( 5+β)(+f s )ξ)); c 9 = 3456 (( +β)((5+β)c s +(c d +c u )(5+β +f s +3βf s ))). At the structure WO, the sum rules can be expressed in the following form c m s L 8/9 E (w)m +c χal 6/7 E (w)m +(c 3 +c 4 )ae 0 (w)+c 5 m s χa L 6/7 E 0 (w)+c 6 m 0 al 4/9 M +c 7 χabl 6/7 M +(c 8 +c 9 )m s a M 4 +c 0 ab M 4 = λ N M N [ (µn ) M +A]e M N /M. (43) (44) (45) For Σ + at WO : c = ( +β)[(+β)c s c u ]; c = 6 ( +β) c s f s φ; c 3 = 4 ( +β)[c s(+β f s +βf s )+c u ( β f s +βf s )]; c 4 = 6 ( +β)[c sf s φ+(+β)c u ](κ ξ); c 5 = 3 ( +β)[c sf s φ+c u ( β f s +βf s )]; c 6 = 4 ( +β)[7(+β)c s +c u ( β 3f s +3βf s )]; c 7 = 44 ( +β) c s f s φ; c 8 = 8 [c s( 3+3β +5f s βf s +5β f s ) c u (3+4β +9β 3f s +0βf s +9β f s )]; c 9 = 7 [c sf s φ(( 0β +β )ξ +4(+β) κ) c u (4(β +6β +3f s +βf s +3β f s )κ (3+4β +9β )(+f s )ξ)]; c 0 = 88 ( +β)((+β)c s c u ( 3 3β f s +βf s )]. (46)

10 0 For Ξ 0 at WO : For Λ at WO : c = ( +β)[(+β)c u c s ]; c = 6 ( +β) c u ; c 3 = 4 ( +β)[c uf s ((+β)f s +β)+c s f s ( (+β)f s +β)]; c 4 = 6 ( +β)[c sf s φ+(+β)c u ](κ ξ); c 5 = 3 ( +β)[c uf s +c s f s φ( (+β)f s +β)]; c 6 = 4 ( +β)[7(+β)c uf s +c s ( (+β)f s 3+3β)]; c 7 = 44 ( +β) c u ; c 8 = 8 [c uf s (( 3+3β )f s +5 β +5β ) c s f s ((3+4β +9β )f s 3+0β +9β )]; c 9 = 7 [c uf s (( 0β +β )ξ +4(+β) κ) c s f s (4((β +6β )f s +3+β +3β )κ (3+4β +9β )(+f s )ξ)]; c 0 = 88 ( +β)((+β)c u c s ( (3+3β)f s +β)]. c = ( +β)((c d +c u )(+β) (5+β)c s ); c = 36 ( +β) ((c d +c u ) c s f s φ); c 3 = 4 ( +β)((c d +c u )(+β( 6+f s ) f s )++c s (0+3f s +β(+5f s ))); c 4 = 96 ( +β)(( +3β)(c d +c u )+(5+3β)c s f s φ)(κ ξ); c 5 = ( +β)((c d +c u )( +β( +f s ) f s )+4c s f s φ); c 6 = 96 ( +β)( (5+β)c s +(c d +c u )(3+7β +5f s +3βf s )); c 7 = 864 ( +β) ((c d +c u ) c s f s φ); c 8 = 6 (( +β)(c d +c u )( 7+β f s +3βf s ) c s ( 5+7f s +β(6+7f s )+β (3+4f s ))); c 9 = 864 (c sf s φ( 4(7+0β +9β )κ+(+6β +35β )ξ) ( +β)(c d +c u )(4( 4 β +5f s +7βf s )κ ( 5+β)(+f s )ξ)); c 0 = 3456 (( +β)((5+β)c s +(c d +c u )(5+β +f s +3βf s )). (47) (48) It is not necessary to list all the coefficients for all 8 members since the coefficients for other members of the octet family can be obtained by appropriate replacements of quark contents in the following way: for proton p, replace s quark by d quark in Σ +. for neutron n, exchange d quark with u quark in proton p, for Σ, replace u quark by d quark in Σ +, for Ξ, replace u quark by d quark in Ξ 0, for Σ 0, change corresponding factors in Ξ 0. Here the conversions between u and d quarks are achieved bysimply switchingtheirchargefactorse u and e d. Moreover, he conversions from s quark to u or d quarks involve setting m s = 0, f = φ =, in addition to the switching of charge factors. his way, they can also provide additional checks for our calculations. Atthis point, wecanmakesomecomparisonswith previous calculations in Ref. [5, 6, 7]. First, we use general interpolating fields where we can vary β to achieve the best match in sum rules. he previous calculations correspond to a fixed value of β=-. his effect was studied in detail in Ref. [9], and it was found that β =. is the optimal value. Most of our results are at β =.. Second, we have checked that our sum rules agree with those in the previous calculations for the most part. For example, for the proton at WE, we completely agree except for the κ ξ term in Eq. (.6) in Ref. [5]. In all of our sum rules, we have the combination κ ξ instead of κ ξ. For the strange members (Σ, Ξ, Λ), they have 8 terms in the OPE, while we have 0 terms. For the sum rules at structure WO, they have only 3 terms, while we have 9 terms. For the sum rules at structure WO, they have only 4 terms, while we have 0 terms. hird, they only analyzed the sum rules at structure WE, while we will examine all the structures. Fourth, we use a completely different analysis method. Before going into the analysis, we would like to point out some relations among the correlation functions (or OPE) based on symmetries, which lead to the same relations in the magnetic moments. In exact SU(3)-flavor symmetry, it is known that the magnetic moments of the octet family are related by (see, for example, Ref. [0]) µ Σ + = µ p, µ Λ = µ n, µ Σ +µ n = µ p, µ Ξ = µ Σ, µ Ξ 0 = µ n. (49) heserelationsareborneoutintheopeofoursumrules if SU(3)-flavor symmetry is enforced. hey are only approximately true since SU(3)-flavor symmetry is broken by the strange quark. Here we have the advantage of

11 studying the symmetry-breaking effects since the terms are explicit in our QCD sum rules. IV. SUM RULE ANALYSIS he sum rules for magnetic moments have the generic form of OPE - ESC = Pole + ransition, or Π mag (QCD,β,w,M ) = λ ( µn ) N M +A e M N /M, (50) where QCD represents all the QCD input parameters. he task then becomes: given the function Π mag with known QCD input parameters and the ability to vary β, find the phenomenological parameters (magnetic moment µ N, transition strength A, coupling strength λ N, and continuum threshold w) by matching the two sides over some region in the Borel mass M. A χ minimization is best suited for this purpose. It turns out that there are too many fit parameters for this procedure to be successful in general. o alleviate the situation, we employ the corresponding mass sum rules which have a similar generic form of OPE - ESC = Pole, or Π mass (QCD,β,w,M ) = λ Ne M N /M, (5) which shares some of the common parameters. Note that the continuum threshold may not be the same in the two sum rules. By taking the ratio of the two equations, we are left with Π mag (QCD,β,w,M ) Π mass (QCD,β,w,M ) = µ N +A. (5) M his is the form we are going to implement. By plotting the two sides as a function of /M, the slope will be the magnetic moment and the intercept the transition strength. he linearity (or deviation from it) of the lefthand side gives an indication of OPE convergence and the role of excited states. he two sides are expected to match for a good sum rule. his way of matching the sum rules has two advantages. First, the slope, which is the magnetic moment of interest, is usually better determined than the intercept. Second, by allowing the possibility of different continuum thresholds, we ensure that both sum rules stay in their valid regimes. We use the chiral-even mass sum rules in Ref. [] which are listed here in the same notation, p L 4/9 E 3 (w )M 6 +p bl 4/9 E (w )M +p 3 m s al 4/9 +p 4 a L 4/9 +p 5 a k v L 4/9 +p 6 m 0a L /7 M = λ N e M N /M. (53) he coefficients for N are: p = 64 (5+β +5β ); p = 56 (5+β +5β ); p 3 = 0; p 4 = 4 (7 β 5β); p 5 = 0; p 6 = 96 (3 β β ). (54) For Λ: p = 64 (5+β +5β ); p = 56 (5+β +5β ); p 3 = 96 ((0 5f s) (6+6f s )β (4+5f s )β ); p 4 = 96 ((4f s 5 6t)+(4+4f s )β +(4f s ++6t)β ); p 5 = 7 ((0f s +)+( 8f s )β (f s +3)β ); p 6 = 88 (( 6f s 3)+(8f s )β +(8f s +5)β ). (55) For Σ: p = 64 (5+β +5β ); p = 56 (5+β +5β ); p 3 = 3 (( 5f s) f s β (+5f s )β ); p 4 = 94 ((4f s ++8t)+4f s β +(4f s 8t)β ); p 5 = 4 ((6f s +) β (6f s )β ); p 6 = 96 ((f s +) β (f s )β ). (56) For Ξ: p = 64 (5+β +5β ); p = 56 (5+β +5β ); p 3 = 3 6 (( f s) f s β (+f s )β ); p 4 = 96 ((5 f s +8t) 0f s β (5+f s +8t)β ); p 5 = 4 f s((f s +6) f s β +(f s 6)β ); p 6 = 94 f s((f s +) f s β +(f s )β ). (57) he function t is defined as t ln M µ γ EM with γ EM the Euler-Mascheroni constant. Use this mass sum rule for the coupling constant for magnetic moment sum rules, and move to the left side, then only the magnetic moment µ N and the transition strength A dependence are on the right hand side. We use the Monte-Carlo procedure first introduced in Ref. [9] to carry out the search which allows a rigorous error analysis. In this method, the entire phase-space of the input QCD parameters is explored simultaneously, and is mapped into uncertainties in the phenomenological parameters. his lead to more realistic uncertainty estimates than traditional approaches. First, a set of randomly-selected, Gaussianlydistributed condensates are generated with a assigned uncertainties. Here we give 0% for the uncertainties of input parameters, and this number can be adjusted to test the sensitivity of the QCD parameters. hen the OPE is constructed in the Borel window with evenly distributed points M j. Note that the uncertainties in the OPE are not uniform throughout the Borel window. hey are larger at the lower end where uncertainties in the higher-dimensional condensates dominate. hus, it is crucial that the appropriate weight is used in the calculation of χ. For the OPE obtained from the k th set of QCD parameters, the χ per degree of freedom is χ k N DF = n B j= [Π OPE k (Mj,β,w,w ) Π Phen k (n B n p ) σope(m j ) (Mj,µ,A)], (58) where Π OPE refers to the LHS of Eq. (5) and Π Phen its RHS. he integer n p is the number of phenomenologi-

12 cal search parameters. In this work, n B =5 points were used along the Borel axis. he procedure is repeated for many QCD parameter sets, resulting in distributions for phenomenological fit parameters, from which errors are derived. In practice, 00 configurations are sufficient for getting stable uncertainties. We used about 000 sets to resolve more subtle correlations among the QCD parameters and the phenomenological fit parameters. his means that each sum rule is fitted 000 times to arrive at the final results. he QCD input parameters are given as follows. he condensates are taken as a = 0.5 GeV 3,b =. GeV 4,m 0 = 0.7 GeV. For the factorization violation parameter, we use κ v =.0. he QCD scale parameter is restricted to Λ QCD = 0.5 GeV. he vacuum susceptibilities have been estimated in studies of nucleon magnetic moments [3, 5, 9, ], but the values vary depending on the method used. We use χ = 6.0 GeV and κ = 0.75, ξ =.5. Note that χ is almost an orderof magnitudelargerthanκandξ, andisthemostimportant of the three. he strange quark parameters are placed at m s = 0.5 GeV, f = 3, φ = 0 [6, 9]. hese input parameters are just central values. We will explore sensitivity to these parameters by assigning uncertainties to them in the Monte-Carlo analysis. Counts Counts u B (u N ) 0.4. ransition Amplitude A V. RESUL AND DISCUSSION We have 4 sum rules in total to analyze: 3 for each member of the octet. For each sum rule, we have in principle 5 parameters to determine: µ, A, w, w, β. But a search treating all five parameters as free does not work because there is not enough information in the OPE. In fact, the freedom to vary β can be used as an advantage to yield the optimal match. We find that β =. gives the best match in most cases. his agrees with the value suggested in Ref. [9]. One exception is the proton: we found a better solution at β = 0.4 than at β =.. Another parameter that can be used to our advantage is the continuum threshold w for the corresponding mass sum rule. We fix it to the value that gives the best solution to the mass sum rule independently. he following values for w are used: for the nucleon, w =.44 GeV; for Λ, w =.60 GeV; for Σ, w =.66 GeV; for Ξ, w =.8 GeV. In this way the magnetic moment sum rule and the mass sum rule can stay in their respective valid Borel regimes. his leaves us with three parameters: µ, A, w. Unfortunately, a three-parameter search is either unstable or returns values for w smaller than the particle mass, an unphysical situation. Again we think this is a symptom of insufficient information in the OPE. So we are forced to fix the continuum threshold w that corresponds to the best match for the central values of the QCD parameters. FIG. 4: Histogram for the proton magnetic moment (top) and transition amplitude(bottom) obtained from Monte-Carlo fits of Eq. 3) at WE for 000 QCD parameter sets. hey are based on 0% uncertainty given to all the QCD input parameters. A. he Sum Rule at WE he results determined this way at the WE structure are displayed in able I. he Borel window is determined by the following two criteria: OPE convergence which gives the lower bound, and ground-state dominance which gives the upper bound. It is done iteratively. For each value of β, we adjust the Borel window until the best solution is found. We see that our calculated magnetic moments agree with experiment fairly well within error bars. We stress that the errors are derived from Monte-Carlo distributions which give the most realistic estimation of the uncertainties. An example of such distributions is given in Fig. 4. We see that they are roughly Gaussian distributions. he central value is taken as the average, and the error is one standard deviation of the distribution. We foundabout 0%accuracyforthe magneticmoments in our Monte-Carlo analysis, resulting from 0% uniform uncertainty in all the QCD input parameters. Of course, the uncertainties in the QCD parameters can be non-uniform. For example, we tried the uncertainty assignments (which are quite conservative) in Ref. [9],

13 3 ABLE I: Results for the magnetic moment of octet baryons from the QCD sum rule in Eq. (3) (structure WE ). he seven columns correspondto, from left toright: particle, β value, Borel region inwhichthetwosides oftheqcdsumruleare matched, continuum threshold, transition strength, extracted magnetic moment in unit of nuclear magnetons, and experimental value. he errors are derived from 000 samples in the Monte-Carlo analysis with 0% uncertainty on all QCD input parameters. β Region w A µ B Exp. (GeV) (GeV) (GeV ) (µ N) (µ N) p to ± ± to ± ±.79 n to ± ± Λ -.. to ± ± Σ to ± ± Σ to ± ± Σ -..3 to ± ± Ξ to ± ± Ξ -.. to ± ± and found about 30% uncertainties in our output. o gain a better appreciation on how the QCD sum rules produce the results, we show Fig. 5, using the proton as an example. here are three graphs in this figure to give three different aspects of the analysis. he first graph shows how the two sides of Eq. (5) match over the Borel window, which should be linear as a function of /M according to the right-hand side of this equation. Indeed, we observe excellent linear behavior from the OPE side (LHS). he match is almost perfect (barely distinguishable between the solid and dotted lines). he slope gives the magnetic moment µ, and the intercept give the transition contribution A. We find that the inclusion of A is important in producing the best match. Also plotted are the individual contributions from u and d quarks. We see that for the proton, the u-quark contribution is the dominant one, which is expected because it is doubled-represented in the proton (uud). We define the slope from an individual quark contribution as the effective magnetic moment of that quark in the particle. he second graph in Fig. 5 shows how the various terms in the OPE contribute to determination of magnetic moments. he M 0 term, which contains the contributions from the condensates χa and b, plays an important role. It is the leading contribution in the region below M <. GeV. For this reason, the sum rule at WE is expected to have good spectral properties. Indeed this is confirmed in the third graph where we plot the three terms in the phenomenological side (pole, transition, and excited) as a function of M. he groundstate pole is dominant (over 70% of the RHS at the low end of the Borel window). he excited-state contribution starts small, then grows with M, as expected from the continuum model. he transition contribution is small in this sum rule. It is consistently smaller than the excitedstate contribution and has a weak dependence on the Borel mass. B. he Sum Rule at WO Next, we analyze the sum rule in Eq. (4) at the structure WO, using the same procedure. able II displays the results extracted from this sum rule. he magnetic moments have larger errors than those from WE : about 5% as opposed to 0%. he agreement with experiment is reasonable(with the exception of Σ ), but not as good as those from WE. We had to searcha wider regionin β to find the best match. he transition contribution (A) are larger for the strange particles, as well as their errors. Fig. 6 shows the details of the analysis in the case of the proton. he matching is very good, as indicated in the top graph. he middle graphs shows that the leading contribution in the OPE (M 4 term) is χa, followed by the quark condensate a (M term). he condensate χab (M 0 term) and ab (M term) are very small in this sum rule. he bottom graph reveals a surprising result: the excited-state dominates over the pole and the transition. As a result, this sum rule is less reliable. his is the reason why the results from this sum rule are not as good as those from WE. his sum rule also shows the importance of checking the individual terms in the phenomenological side, in addition to looking at the best matchofthetwosides. Inthiscase,thereisnopoledominance, even though the leading terms is non-perturbative and the match is almost perfect. C. he Sum Rule at WO Finally, we present the results from the sum rule in Eq. (45) at the structure WO in able II. he agreement with experiment is not as good as the other two sum rules. For example, Σ and Ξ have the wrong sign. Fig. 7 shows the details of the analysis for the proton. he matching is very good, as indicated in the top graph. he middle graph shows that the leading contri-

14 4 ABLE II: Similar to able I, but for the QCD sum rule in Eq. (4) (structure WO ). β Region w A µ B Exp. (GeV) (GeV) (GeV ) (µ N) (µ N) p -.4 to ±.67 ±.79 n -.. to ± ± Λ -..0 to ±0.3 - ±0.7 - Σ to ± ± Σ to ± ± Σ -. to ± ± -.6 Ξ to ± ± Ξ to ± ±0.9-5 butioninthe OPEisM withacoefficientofχa, followed by the M 0 term. he /M term is slightly negative, while the /M 4 term is very small. he bottom graph shows that the excited-state dominates over the pole and the transition, like the sum rule from WO, but the relative size of the pole is much larger. Since the WO sum rule has power corrections up to /M 4, it is expected to be more reliable than the WO sum rule. But our analysis shows that this advantage is offset by the smallness of the /M and /M 4 terms. As a result, the reliability of the WO sum rule is about the same as the WO sum rule. We have performed the same analysis for all the members and all three structures. Fig. 8 shows the graphs for the neutron. In this case, the slope is negative. Again, the sum rule at WE has excellent convergence properties. he WO sum rule has a good match, but the pole is less than the excited-state. he WO sum rule does not have a good match. Fig. 9 shows the case for the Σ 0, which has all three quark contributions (u, d and s). Overall, based on the quality of the match, the broadness ofthe Borelwindow and its reachinto the lowerend, the size of the continuum contribution, and the OPE convergence, we find that the sum rule at WE is the most reliable of the three sum rules. D. Some physics discussions Based on the results of our comprehensive analysis, we conclude that the QCD sum rules at WE in Eq. (5) are the most reliable. Here we discuss some physics implications of the results extracted from them (able I). First we look at some ratios of magnetic moments listed in able IV. Here we compare our results with those from the SU(6) symmetry, lattice calculations [] and experiment. From the table, we see that the QCD sum rule results compare well against other approaches and experiment. hey agree a little better with experiment than the lattice results. Furthermore, our QCD rum rule results are an improvement over the ratios from previous QCD sum rule calculations in Ref. [5]. Next, we consider a few sum rules among magnetic moments that have been discussed in literature. hey reveal interesting quark dynamics in the baryons. hey are mostly based on SU(6)-symmetry considerations in the quark model. We begin with the sum rule [3] 3(p Σ + ) Ξ Ξ 0 = (p+3λ). (59) p It assumed baryon independence of quark moments: the independence of which baryon the same quark is in, a concept first mentioned by Franklin [4]. In other words, eachquarkisnotsensitivetotheenvironmentitresidesin and quarks in different spin states had the same effective moments. In this sum rule, it also assumed that the s quark moment in the Λ is the same as in the Σ and Ξ, even though the spin states are different. his rum rule isviolatedthemost (byafactorof5)usingpresentvalues of magnetic moments. his large violation is mostly due tothe smalldifferenceinthe denominatorξ Ξ 0 which magnifies the apparent discrepancy. So this is not a good way of testing baryon independence. Our determination of the ratio of the left hand side over the right hand side (LHS/RHS) is 7.4(.0), compared to the experiment value of 5.7(4) and the lattice calculation of 4.(.5). he error are added in quadrature in forming these ratios. One interesting sum rule [9] which is fairly accurately satisfied is p+n = 3Λ+ (Σ+ +Σ ) (Ξ 0 +Ξ ). (60) It is derived using SU(3)-flavor symmetry to characterize the Λ wave function Our determination of the LHS/RHS ratio is.4(79), compared to the experimental values of.(5) and the lattice calculation of.45(7). Another sum rule, first derived by Franklin [4], p n = Σ + Σ +Ξ Ξ 0, (6) is another test of baryon independence of the quark moments. he strange quarks approximately cancel, leaving only u and d quarks. We give the ratio of.0(8),

15 5 ABLE III: Similar to able I, but for the QCD sum rule in Eq. (45) (structure WO ). β Region w A µ B Exp. (GeV) (GeV) (GeV ) (µ N) (µ N) p -..3 to ± ± n -..7 to.9.60 ± ± Λ -..3 to ± ±0.7 - Σ to ± ±0..45 Σ to ± ± Σ -..4 to ±0.0 ± Ξ to ± ± Ξ -..4 to ± ±0. -5 while the experiment measurements gives.33(0) and lattice results is.07(8) for this ratio. Finally, the Sachs sum rule [6] 3(p+n) = Σ + Σ Ξ +Ξ 0 (6) is satisfied by the more general extension of SU(6) symmetry and is another test of baryon independence of quark moments. his sum rule is just the sum of the two separate sum rules proposed earlier by Franklin [5]. Our result gives the ratio of 6(35) which is in agreement with the experiment measurements of 8(). Whereas it yields an opposite violation with a ratio of.9(0) from lattice moments. his contradiction is possibly because the Sachs sum rule may be sensitive to dynamics not included in the lattice calculation. E. Individual quark contributions function of /M. We call this the raw individual quark contributions to the magnetic moments. able V gives the result of raw individual u, d and s quark sector contributions to the magnetic moments from the QCD sum rules at WE. It is compared with the lattice QCD result in []. he lattice results were rescaled later in [7] and we use the rescaled results for the comparison. Our results agree with the lattice results reasonably well. he biggest discrepancy is that our light quark moments in Ξ 0 and Ξ are small. In Σ 0, and Λ the total of u and d quark moments agree very well with lattice data for light quark moments. Note that only a combined number for light quark (µ l = µ u + µ d ) was given on the lattice []. o compare with our separated u and d quark numbers, we break up the lattice number by using the relation µ u = µ d. o gain a deeper understanding of the dynamics, it is useful to consider the individual quark sector contributions to the magnetic moment. In our approach, we can easily dial individual quark contributions to the QCD sum rules. For example, to turn off all u-quark contributions, we set the charge factor c u = 0. o turn off all s-quark contributions, we set c s = 0, m s = 0, f =, and φ =. We can extract a number corresponding to each quark contribution from the slope of Eq. (5) as a ABLE IV: A comparison of selected ratios of magnetic moments. Ratio Sum Rule SU(6) Lattice Exp. result symmetry results n/p -0.70(9) /3-3(5) -8 Σ /Σ (7) /3-0.37(3) Ξ /Ξ (6) / 0.58(5) 0.5 Ξ 0 /Λ.05(56).4(5).04 Ξ /Λ.4(3).0().3 In the simple quark model [8], the magnetic moment of the proton is given by p = 4 3 µ u 3 µ d. In the SU() limit, there is the relation µ u = µ d. Our results for Σ 0 and Λ, the two uds structure baryons, agree with this very well. his relation also suggests that the ratio of the u quark and d quark contributions in proton is 4 3 µ u/( 3 4 µ d) = 8. Our QCD sum rule ratio is for β =., which agrees with the lattice ratio of 0.3(7) at κ. For neutron the quark model ratio is 4 3 µ d/( 3 µ u) =, and our result is 3.64(40) for this ratio, which differs from the SU(6) spin-flavor symmetry prediction. he lattice ratio is.6. he bigger ratio represents an enhancement of the doubly represented quark contribution. hese results suggest that although the total magnetic moments ratios agree with SU(6) ratios (see able IV), the underlying quark dynamics are really quite different from different individual quark contributions. o include the possible nonstatic effect, Franklin [9] proposed a generalization of SU(6) results so that non-

16 6 ABLE V: Individual quark contributions and total magnetic moments in unit of nuclear magnetons extracted from the QCD sum rules at WE (denoted by the subscript SR) are compared with those from a lattice QCD calculation (denoted by the subscript LA) [, 7]. µ u SR µ u LA µ d SR µ d LA µ s SR µ s LA µ otal SR µ otal LA p.46(5).59(4) 0.38(4) 0.5(9) 0 0.8(6).79 n -0.4(9) -0.3(0) -.53(3) -.3(4) (5) -.60() Λ (7) 0.09(6) -0.09(4) -0.03() -3(5) -0.54(7) -0.56(5) -0.50(7) Σ +.06(5).33(9) () 0.08(5).3(5).37(8) Σ 0.09() 7(9) -0.54(6) -0.9(3) 0.4() 0.08(5) 9(7) 5(6) Σ (0) -.4(4) 0.0() 0.08(5) -.6(0) -.08(0) Ξ 0 0.0() -0.44(6) (7) -0.73(6) -.5(5) -.7(9) Ξ () 0.(3) -7(6) -0.73(6) -4(6) -0.5(7) static components are the same for each octet baryon: p = 4 3 µ u 3 µ d, n = 4 3 µ d 3 µ u, Σ + = 4 3 µ u 3 µ s, Σ = 4 3 µ d 3 µ s, Ξ 0 = 4 3 µ s 3 µ u, Ξ = 4 3 µ s 3 µ d, Σ 0 = 3 µ u + 3 µ d 3 µ s, Λ = µ s (63) where quark symbols refer to quark moment contributions including nonstatic effect. he magnetic moment contribution of the unlike quark in the baryon is primed. In order to compare with the effective quark moments defined in the quark model, we convert our raw quark moments in able quark in a similar fashion. ake the proton for example, we define effective quark moments µ u and µ d from QCD sum rules by u SR = 4 3 µ u, d SR = 3 µ d. One can define effective moments for quarks in other baryons in a similar manner. In the proton, µ u is.84(8) while in Σ +.54(8) and in Σ 0 is.63(8). he magnetic moments list here are all in unit of µ N. In the neutron, µ u is.5(7), which is smaller than µ u. his is an example that effective quark magnetic moment is sensitive to the environment the quark resides in. We find the following relation from our results µ p u < µσ0 u < µσ+ u. (64) For µ d, in n it is.5() while for Σ 0 is (8) and for Σ is 0.93(8). From p we can get µ d =.5(). It has a similar relation for the absolute value µ n d < µσ0 d < µ Σ d. (65) For µ s, it is (6) in Ξ 0 and 0.50(5) in Ξ. For µ s, it is 0.45(5), 0.4(5) and 0.30(5) in Σ +, Σ 0 and Σ, respectively. FromΛwecangetµ s about 3(5). For these individual quark effective moments, we notice that in general we have the following relation µ s < µ d < µ u. (66) It is expected because of the quark mass effects, which is analogous to those seen in the electric properties. Another way of looking at the individual effective quark moments is by expressing them in terms of baryon magnetic moments using Eq.(63) and isospin symmetry. For example, the d-quark effective moment can be expressed as µ d = 4 (p+n) = 4 (Σ Σ + ). (67) Our result indicates µ d = 0.9(4) < 6(7) using the baryon moments in able V or in able I. It agrees well with the experimental moments µ d = 0.98 < 94(7) and the lattice result of µ d =.00(5) < 6(6). Similarly the d quark effective moment is µ d = p+n = Ξ0 Ξ. (68) We have µ d =.(39) < 0.5(8). It also agree very well with the experimental moments µ d =.003 < 0.578(6) although the two sides of this sum don t agree well. he strange quark can be isolated as µ s = 4 (Ξ0 +Ξ ), µ s = Σ+ Σ, (69) µ s = Λ δ. Our results are compared with experiment and lattice calculation in able VI. Not only the results agree with the experiment and lattice calculation, they also roughly agree with the strange quark moments from individual quark moments. Now, we can look at the quark moment difference such as µ s µ d and µ s µ d. Accordingto simple quarkmodel, the difference should be the same and approximately as In fact, from µ s µ d = 3(p Σ+ ), (70)