arxiv: v1 [hep-ph] 7 Mar 2011

Transcription

1 YITP-11-8 Magnetic Susceptibility o the Quark Condensate and Polarization rom Chiral Models Marco Frasca 1, and Marco Ruggieri, 1 via Erasmo Gattamelata, 3, Roma, Italy Yukawa Institute or Theoretical Physics, Kyoto University, Kitashirakawa Oiwake-cho, Sakyo-ku, Kyoto , Japan We compute the magnetic susceptibility o the quark condensate and the polarization o quarks at zero temperature and in a uniorm magnetic background. Our theoretical ramework consists o two chiral models that allow to treat sel-consistently the spontaneous breaking o chiral symmetry: the linear σ model coupled to quarks, dubbed quark-meson model, and the Nambu-Jona-Lasinio model. We also perorm analytic estimates o the same quantities within the renormalized quarkmeson model, both in the regimes o weak and strong ields. Our numerical results are in agreement with the recent literature; moreover, we conirm previous Lattice indings, related to the saturation o the polarization at large ields. arxiv: v1 [hep-ph] 7 Mar 011 PACS numbers: 1.38.Aw,1.38.Mh Keywords: Magnetic Susceptibility o the Quark Condensate, Eective Chiral Models, QCD in Strong Magnetic Background. I. INTRODUCTION One o the most attractive aspects o the vacuum o Quantum Chromodynamics (QCD), the theory o strong interactions, is its non-perturbative nature. Mainly by means o Lattice QCD simulations at zero quark chemical potential [1 5] it is established that two crossovers take place in a narrow range o temperature; one or quark deconinement, and another one or the (approximate) restoration o chiral symmetry. Besides, powerul analytic and semi-analytic techniques have been developed to understand the coupling between chiral symmetry restoration and deconinement, see [6, 7] and reerences therein. Moreover, Lattice QCD and Operator Product Expansion (OPE in the ollowing) o the correlators o hadronic currents show that the QCD vacuum can be characterized by several quark, gluon and mixed condensates [8]. A ruitul theoretical approach to the physics o strong interactions, which is capable to capture some o the nonperturbative properties o the QCD vacuum, is the use o eective chiral models. One o them is the Nambu- Jona-Lasinio (NJL) model [9] (see Res. [10] or reviews), in which the QCD gluon-mediated interactions are replaced by eective interactions among quarks, which are built in order to respect the global symmetries o QCD. Under some approximations, it is possible to derive the NJL model eective interaction kernel rom irst principles QCD, see [11, 1]. Besides, a linear σ model augmented with a Yukawa-type coupling to quarks, named Quark-Meson model (QM model in the ollowing), has been developed as an eective model o QCD [13, 14]. In this model, quartic meson sel-couplings allow to absorb the cuto dependence o the coupling constants; moreover, ollowing an idea by Weinberg [15], tree level prop- marcorasca@mclink.it ruggieri@yukawa.kyoto-u.ac.jp agating mesons are suicient to avoid triviality o the NJL model in 3+1 dimensions in the one-loop approximation [16] (see [14] or an excellent discussion about these points). The chiral models are widely used to map qualitatively, and to some extent also quantitatively, the phase diagram o strongly interacting matter along several directions like temperature, chemical potential, isospin chemical potential and external ields [17 30]. It is thus interesting to use these models to compute other quantities, which are related to the QCD vacuum condensates. As a matter o act, the chiral models allow or a sel-consistent treatment o the spontaneous chiral symmetry breaking in the vacuum. Thereore, ater the ground state properties (i.e., the chiral condensate in the NJL model, or the expectation value o the σ-ield in the quark-meson model) under external actors are computed, it is straightorward to estimate numerical values o other condensates, using the one-loop propagators o the theory. It has been realized that external ields can induce QCD condensates that are absent otherwise [31]. O particular interest or this article is magnetic moment, Σ µν where denotes the ermion ield o the lavor th, and Σ µν = i(γ µ γ ν γ ν γ µ )/. At small ields one can write, according to [31], Σ µν = χ Q eb, (1) and χ is a constant independent on lavor, which is dubbed magnetic susceptibility o the quark condensate. In [31] it is proved that the role o the condensate (1) to QCD sum rules in external ields is signiicant, and it cannot be ignored. The quantity χ has been computed by means o special sum rules [31 35], OPE combined with Pion Dominance [36], holography [37, 38], instanton vacuum model [39], analytically rom the zero mode o the Dirac operator in the background o a SU() instanton [40], and on the Lattice in two color quenched simulations at zero and inite temperature [41]. It has also been suggested that in the photoproduction o lep-

2 ton pairs, the intererence o the Drell-Yan amplitude with the amplitude o a process where the photon couples to quarks through its chiral-odd distribution amplitude, which is normalized to the magnetic susceptibility o the QCD vacuum, is possible [4]. This intererence allows in principle to access the chiral odd transversity parton distribution in the proton. Thereore, this quantity is interesting both theoretically and phenomenologically. The several estimates, that we briely review in Section III, lead to the numerical value o χ as ollows: χ = MeV. () A second quantity, which embeds non-linear eects at large ields, is the polarization, µ, deined as µ = Σ, Σ = Σ 1, (3) which has been computed on the Lattice in [41] or a wide range o magnetic ields, in the ramework o twocolor QCD with quenched ermions. At small ields µ = χq eb naturally; at large ields, non-linear eects dominate and an interesting saturation o µ to the asymptotic value µ = 1 is measured. According to [41] the behavior o the polarization as a unction o eb in the whole range examined, can be described by a simple inverse tangent unction. Besides, magnetization o the QCD vacuum has been computed in the strong ield limit in [43] using perturbative QCD, where it is ound it grows as BlogB. In this article, we compute the magnetic susceptibility o the quark condensate by means o the NJL and the QM models. This study is interesting because in the chiral models, it is possible to compute sel-consistently the numerical values o the condensates as a unction o eb, once the parameters are ixed to reproduce some characteristic o the QCD vacuum. We irstly perorm a numerical study o the problem, which is then complemented by some analytic estimate o the same quantity within the renormalized QM model. Moreover, we compute the polarization o quarks at small as well as large ields, both numerically and analytically. In agreement with the Lattice results [41], we also measure a saturation oµ to one at largeields, in the caseothe eective models. Our results push towards the interpretation o the saturation as a non-artiact o the Lattice. On the contrary, we can oer a simple physical understanding o this behavior, in terms o lowest Landau level dominance o the chiral condensate. As a matter o act, using the simple equations o the models or the chiral condensate and or the magnetic moment, we can show that at large magnetic ield µ has to saturate to one, because in this limit the higher Landau levels are expelled rom the chiral condensate; as a consequence, the ratio o the two approaches one asymptotically. We also obtain a saturation o the polarization within the renormalized QM model. There are some dierences, however, in comparison with the results o the nonrenormalized models. In the ormer case, the asymptotic value o µ is charge-dependent; moreover, the interpretation o the saturation as a lowest Landau level (LLL) dominance is not straightorward, because the renormalized contribution o the higher Landau levels is important in the chiral condensate, even in the limit o very strong ields. It is possible that the results obtained within the renormalized model are a little bit ar rom true QCD. As a matter o act, in the renormalized model we assume that the quark sel-energy is independent on momentum; thus, when we take the limit o ininite quark momentum in the gap equation, and absorb the ultraviolet divergences by means o counterterms and renormalization conditions, we implicitly assume that that quark mass at large momenta is equal to its value at zero momentum. We know that this is not true, see or example [44, 45]: even in the renormalized theory, the quark sel-energy naturally cuts o the large momenta, leading to LLL dominance in the traces o quark propagator which are relevant or our study. Nevertheless, it is worth to study this problem within the renormalized QM model in its simplest version, because it helps to understand the structure o this theory under the inluence o a strong magnetic ield. In our calculations we neglect, or simplicity, the possible condensation o ρ mesons at strong ields [46, 47]. Vector meson dominance [46] and Sakai-Sugimoto model [48] suggest or the condensation a critical value o eb c m ρ 0.57 GeV, where m ρ is the ρ meson mass in the vacuum. Beside these, a NJL-based calculation within the lowest Landau level (LLL) approximation [47] predicts ρ meson condensation at strong ields as well, even i in the latter case it is hard to estimate exactly eb c, mainly because o the uncertainty o the parameters o the model. It would certainly be interesting to address this problem within our calculations, in which not only the LLL but also the higher Landau levels are considered, and in which the spontaneous breaking o chiral symmetry is kept into account sel-consistently. However, this would complicate signiicantly the calculational setup. Thereore, or simplicity we leave this issue to a uture project. The plan o the article is as ollows. In Section II we describe the QM and NJL models and ix our notation. In Section III we discuss our numerical results or the polarization o quarks, and compute the magnetic susceptibility o the quark condensate. In Section IV we compute the renormalized Quantum Eective Potential (QEP) o the QM model in a magnetic background, in the one-loop approximation, and compute analytically the solution o the gap equation in the weak ield case, and semi-analytically in the strong ield limit. We then use the results to estimate µ and χ. Finally, in Section V we draw our conclusions.

3 3 II. CHIRAL MODELS COUPLED TO A MAGNETIC FIELD In the irst part o this article, we derive numerical results or the spin polarization o quarks in a magnetic ield, and or the magnetic susceptibility o the quark condensate, using two eective chiral models: the Nambu- Jona-Lasinio (NJL) model, and the Quark-Meson (QM) model. We describe the two models in some detail in this Section. We work in the Landau gauge and take the magnetic ield along the z axis, B = (0,0,B). A. Quark-Meson model In the QM model, a meson sector described by the linear sigma model lagrangian, is coupled to quarks via a Yukawa-type interaction. The model is renormalizable in D = 3+1dimensions. However, since we adopt the point o view o it as an eective description o QCD, it is not necessary to use the renormalized version o the model itsel. On the contrary, it is enough to ix an ultraviolet scale to cuto the divergent expectation values; the UV scale is then chosen phenomenologically, by requiring that the numerical value o the chiral condensate in the vacuum obtained within the model, is consistent with the results obtained rom the sum rules [49]. This is a rough approximation o the QCD eective quark mass, which smoothly decays at large momenta [44, 45]. In Section IV we will use a renormalized version o the model, to derive semi-analytically some results in the two regimes o weak and strong ields. The lagrangian density o the model is given by L = q[id µ γ µ g(σ +iγ 5 τ π)]q + 1 ( µσ) + 1 ( µπ) U(σ,π). (4) In the above equation, q corresponds to a quark ield in the undamental representation o color group SU(3) and lavor group SU(); the covariant derivative, D µ = µ Q ea µ, describes the coupling to the background magnetic ield, where Q denotes the charge o the lavor. Besides, σ, π correspond to the scalar singlet and the pseudo-scalar iso-triplet ields, respectively. The potential U describes tree-level interactions among the meson ields. In this article, we take its analytic orm as U(σ,π) = λ 4 ( σ +π v ) hσ, (5) where the irst addendum is chiral invariant; the second one describes a sot explicit breaking o chiral symmetry, and it is thus responsible or the non-zero value o the pion mass. For h = 0, the interaction terms o the model are invariant under SU() V SU() A U(1) V. This group is broken explicitly to U(1) 3 V U(1)3 A U(1) V i the magnetic ield is coupled to the quarks, because o the dierent electric charge o u and d quarks. Here, the superscript 3 in the V and A groups denotes the transormations generated by τ 3, τ 3 γ 5 respectively. Thereore, the chiral group in presence o a magnetic ield is U(1) 3 V U(1)3 A. This group is then explicitly broken by h-term to U(1) 3 V. In this article, we restrict ourselves to the large-n c (that is, one-loop) approximation, which amounts to consider mesons as classical ields, and integrate only over ermions in the generating unctional o the theory to obtain the Quantum Eective Potential(QEP). As a matter o act, quantum corrections arising rom meson bubbles are suppressed o a actor 1/N c with respect to case o the ermion bubble. In the integration process, the meson ields are ixed to their classical expectation value, π = 0 and σ 0 (in particular, σ has the quantum numbers o the chiral condensate, qq ). The physical value o σ will be then determined by minimization o the QEP. To compute QEP in presence o a magnetic background, we use the Leung-Ritus-Wang method [50] which allows to write down the quark propagator or the lavor in terms o Landau levels, S (x,y) = k=0 dp0 dp dp 3 (π) 4 E P (x)λ k i P γ MĒP(y), (6) where E P (x) corresponds to the eigenunction o a charged ermion in magnetic ield, and ĒP(x) γ 0 (E P (x)) γ 0. In the above equation, P = (p 0,0,Q k Q eb,p 3 ), (7) where k = 0,1,,... labels the k th Landau level, and Q sign(q ), with Q denoting the charge o the lavor ; Λ k is a projector in Dirac space which keeps into account the degeneracy o the Landau levels; it is given by Λ k = δ k0 [P + δ Q,+1 +P δ Q, 1 ]+(1 δ k0 )I, (8) where P ± are spin projectorsand I is the identity matrix in Dirac spinor indices. At the one-loop level, the QEP then reads V = λ ( σ v ) hσ 4 N c Q eb π k β k dp3 π ω k(p 3 ). (9) In the above equation we have deined ω k (p 3 ) = p 3 +k Q eb +m q, (10) with m q = gσ; β k = δ k0 counts the degeneracy o the Landau levels. The one-loop ermion contribution, which corresponds to the last addendum in the r.h.s. o Eq. (9), is divergent in the ultraviolet. In order to regularize it, we adopt a

4 4 smooth regulator U Λ as in [6], which is more suitable, rom the numerical point o view, in our model calculation with respect to the hard-cuto which is used in analogous calculations without magnetic ield. In this article we chose U Λ = Λ N Λ N +(p z + Q eb k) N, N = 5. (11) The (more usual) 3-momentum cuto regularization scheme is recovered in the limit N ; we notice that, even i the choice N = 5 may seem arbitrary to some extent, it is not more arbitrary than the choice o the hard cuto scheme, that is, o a regularization scheme. In eective models, the choice o a regularization scheme is a part o the deinition o the model itsel. Momentum integrals are understood as ollows: n β n dp π n β n dp π U Λ. (1) Once the expectation value o σ is computed as a unction o eb by a procedure o minimization o V, we compute the expectation values that are relevant to the context. To begin with, we consider the chiral condensate or the lavor, = Tr[S (x,x)], with S given by Eq. (6). It is straightorward to derive the relation = N c Q eb π k=0 β k dp3 π m q ω k (p 3 ), (13) where the divergent integral on the r.h.s. o the above equation has to be understood regularized as in (1). From Eq. (13) we notice that the prescription (1) is almost equivalent to the introduction o a running eective quark mass, m q = gσθ ( Λ p 3 k Q eb ), (14) that can be considered as a rough approximation to the eective running quark mass in QCD [45] which decays at large quark momenta, see also the discussion in [44]. Once the scale Λ is ixed, the Landau levels with n 1 are removed rom the chiral condensate i eb Λ. Next we turn to the magnetic moment or the lavor, where Σ µν = Tr[Σ µν S (x,x)], (15) Σ µν = 1 i (γµ γ ν γ ν γ µ ), (16) is the relativistic spin operator. We take B = (0,0,B); in this case, only Σ 1 Σ is non-vanishing. Using the properties o γ matrices it is easy to show that only the Lowest Landau Level (LLL) gives a non-vanishing contribution to the trace: Σ = N c Q eb π dp3 π m q ω 0 (p 3 ), (17) where ω 0 = ω k=0. Once again, the divergent integral on the r.h.s. o the above equation has to be understood regularized ollowing the prescription in Eq. (1). For the QM model, the parameters are tuned as ollows. We ixthe constituent quarkmassin the vacuumto a phenomenological value, m q = 335 MeV; urthermore, we require that in the vacuum the ollowing condition holds: V(σ,B = 0) σ = 0, (18) σ=π which implies that σ = π in the vacuum, thus m q = g π. As a consequence we ind g = 3.6. The parameter h is ixed by the condition h = π m π, where m π is the average value o the pion mass in the vacuum and π = 9.4 MeV. We have then h = 0.047m 3 q. To determine v and λ we solve simultaneously Eq. (18) and m σ = V σ, (19) σ=π with m σ = 700 MeV. The divergences in these two equations are cured with the prescription (1). In the UV regulator we chose Λ = 560 MeV which implies ūu = ( 31 MeV) 3. This results in the values λ = 4.67 and v = 1.8m q. Beore going ahead, we notice that in the non-renormalized model, the 1-loop ermion contribution regularized at the scale Λ is included into the conditions (18) and (19). B. Nambu-Jona-Lasinio model ThequarklagrangiandensityotheNJLmodelisgiven by [ L = q(iγ µ D µ m 0 )q +G ( qq) +(i qγ 5 τq) ] ; (0) hereq isthequarkdiracspinorin theundamentalrepresentation o the lavor SU() and the color group; τ correspond to the Paulimatrices in lavorspace. Asum over color and lavor is understood. Once again, the covariant derivative embeds the QED coupling o the quarks with the external magnetic ield. The explicit sot breaking o chiral symmetry in this model is achieved by introducing a current quark mass, m 0, whose numerical value is ixed by itting the vacuum pion mass. The NJL model is not renormalizable in the usual sense or D > space-time dimensions; it is renormalizable in D < 4 in the one-loop approximation. At D = 4 it might represent a trivial theory o non-interacting bosons i the ultraviolet cuto is allowed to be ininite [14, 16]. However, adopting the point o view o the NJL model as an eective (and rough) description o low-energy QCD, it is not necessary torequirethecutotobeininite: the cutocanbe interpreted as a physical quantity which makes the eective quark mass to be constant at small quark momenta, and vanishing at large momenta, thus mimicking roughly the

5 5 running eective mass o real QCD [45]. It is worth to notice that a non-local interaction among quarks with the same quantum numbers o the our-ermion local term in Eq. (0) can be derived by the QCD action under some approximation, taking the low energy limit o the latter [7, 11]. For simplicity, we treat here only the case o local interaction. Once again, the one-loop ermion contribution can be obtained within the Leung-Ritus-Wang method: V = GΣ N c =u,d Q eb π + β n n dp π ω n(p), (1) with ω n (p) = p +n Q eb +Mq, and M q = m 0 GΣ with Σ = ūu + dd. Equations (13) and (15) are still valid in the NJL model, with the replacement m q M q. A comparison o the QEPs o the QM model and o the NJL model shows that at the one-loop level, the two models dier only or the classic part o the eective meson potential, and or the deinition o the constituent quark mass. The parameters in the NJL model are ixed as ollows. The bare quark mass is computed by virtue o the Gell- Mann-Oakes-Renner relation, m π π = m 0 ūu, which is satisied in the NJL model [10]. Then, two equations are solved simultaneously: one or π and one or the chiral condensate in the vacuum at zero magnetic ield strength. Requiring that π = 9.4 MeV and ūu = ( 53 MeV) 3, this procedure gives the numerical values o G and Λ. Thus we ind Λ = MeV, G =.0/Λ and m 0 = 5 MeV. Condensates Π eb GeV eb GeV FIG. 1. Upper panel. Chiral condensates o u-quarks (red) and d-quarks (blue), in units o the same quantities at zero magnetic ield, as a unction o the magnetic ield. Lower panel. Expectation value o the magnetic moment operator, in units o 3 π, as a unction o eb. Data correspond to the QM model. III. NUMERICAL RESULTS In this Section, we collect our numerical results which are relevant or the computation o the magnetic susceptibility o the chiral condensate, and o the spin polarization. From the numerical point o view, it is more convenient to compute irstly the latter; then, a it o the polarization data at small ields will enable to extract the value o the magnetic susceptibility o the quark condensate. For what concerns data about chiral condensate, magnetic moment and polarization, we plot results only or the QM model, since the results obtained within the NJL model are qualitatively very similar to those obtained within the ormer model. Then, we will give the inal result or the chiral magnetization or the two models considered here. A. Polarization The physical value o the (total) chiral condensate Σ or the NJL model, and o σ or the QM model, are obtained numerically by a minimization procedure o the QEP (1) or any value o eb. Then, we make use o Eqs. (13) and (15), with the replacement (1), to compute the chiral condensate and the magnetic moment or each lavor. In the upper panel o Fig. 1, we plot the chiral condensates or u and d quarks, as a unction o eb, or the QM model. The magnetic ield splits the two quantities because o the dierent charge or the two quarks. The small oscillations, which are more evident or the case o the u-quark, are an artiact o the regularization scheme, and disappear i smoother regulators are used, see the discussion in [9]. In the regime o weak ields, our data are consistent with the scaling eb /M where M denotes some mass scale; in the strong ield limit we ind instead eb 3/. The behavior o the quark condensate as a unction o magnetic ield is in agreement with the magnetic catalysis scenario [7, 51]. In the lower panel o Fig. 1 we plot our data or the expectation value o the magnetic moment. Data correspond to the QM model (or the NJL model we obtain similar results). At weak ields, Σ eb as expected rom Eq. (15). In the strong ield limit, non-linearity arise because o the scaling o quark mass (or chiral con-

6 u quark Μ u quark d quark Μ Trig Clas Quan eb GeV FIG.. Polarization o u-quarks (red) and d quarks (blue) as a unction the magnetic ield strength, or the QM model. densate); we ind Σ eb 3/ in this limit. In Fig. we plot our results or the polarization. Data are obtained by the previous ones, using the deinition (3). At small ields, the polarization clearly grows linearly with the magnetic ield. This is a natural consequence o the linear behavior o the magnetic moment as a unction o eb or small ields, see Fig. 1. On the other hand, within the chiral models we measure a saturation o µ at large values o eb, to an asymptotic value µ = 1. This conclusion remains unchanged i we consider the NJL model, and it is in agreement with the recent Lattice indings [41]. It should be noticed that, at least or the u quark, saturation is achieved beore the expected threshold or ρ meson condensation [46 48]. Thereore, our expectation is that our result is stable also i vector meson condensation is considered. The saturation to the asymptotic value µ = 1 o polarization is naturally understood within the models we investigate, as a LLL dominance in the chiral condensate (i.e., ull polarization). As a matter o act, Σ and turn to be proportional in the strong ield limit, since only the LLL gives a contribution to the the latter, compare Eq. (13) and (17) which imply µ = 1 HLL, () where HLL corresponds to the higher Landau levels contribution to the chiral condensate. In the strong ield limit HLL 0 because o Eq. (1) which is a rough approximation to the QCD eective quark mass, as depicted in Eq. (14); hence, µ has to approach the asymptotic value µ = 1. On the other hand, in the weak ield limit HLL and the proportionality among Σ and is lost. In Re. [41], data o polarization are it by means o Μ d quark eb GeV Trig Clas Quan eb GeV FIG. 3. Upper panel. Polarization data or u quark (grey dots) and itting curves. See the text or details. Lower panel. Polarization data and itting unctions or d quark. Results correspond to the QM model. three dierent unctions, namely µ clas = µ coth 3χQ eb µ µ 3χQ eb, (3) µ quan = µ coth χq eb coth χq eb µ µ, (4) µ trig = µ π arctan πχq eb µ. (5) The three unctions above share the behavior at the origin, µ χq eb, and the asymptotic one, µ µ. In [41] a two-parameterit on χ and µ is perormed; the model with the lowest chi-squared per degrees o reedom is represented by the trigonometric one, Eq. (5). Inspired by these results we have tried to it our data using the same unctions (3)-(5). In our model calculation, the asymptotic value µ = 1 is achieved straightorwardly, thereore it is enough to perorm a one-parameter it leaving χ as a ree parameter. The results o this procedure are collected in Fig. 3, where we plot our data o quark s polarization as a unction o eb (gray dots) and the three itting unctions with µ = 1. From Fig. 3 we read that both the trigonometric and the classic unctionsdonotadaptwelltoourdata(orthenjlmodelwe obtain similar results). In particular, both o the aore-

7 7 mentioned unctions overestimate χ, and reach slowly the asymptotic value. On the other hand, the quantum unction in Eq. (4) oers a better description o our data in the whole range o eb examined here. We have checked that the value o χ obtained within this it overestimates the value obtained by a weak-ield linear it only o the 10%, see below. Moreover, the it unctions smoothly ollow the large eb data to the asymptotic value µ = 1. We conclude that within our model, the quantum itting unction is a more aithul representation o data. The dierence with Lattice data is probably due to the act that the latter ones are obtained with quenched (hence, non-dynamical) ermions, while in our case dynamical ermions are considered. It will be interesting, thereore, to compare the model predictions with Lattice data obtained with other kinds o ermions in the near uture. Beore going ahead, it is interesting that our data on polarization, and our interpretation o the saturation o the latter at large ields, gives a quantitative estimate o the goodness o the LLL approximation or the models at hand. In particular, at eb c 0.6 GeV, which in turns is the estimated critical value or vector meson condensation [46, 48], the LLL almost dominates the chiral condensate o the u quark; on the other hand, at this value o eb the LLL dominance has not yet been achieved or the d quark condensate; but the higher Landau levels at this value o eb give a contribution almost to the 0% o the chiral condensate. This is an indirect check o the consistence o the LLL approximation used in [47], which in turn should be good within a 0% accuracy. B. Magnetic susceptibility o the quark condensate At small ields µ = χq eb rom Eq. (1). Hence, we use the data on polarization at small ields, to obtain the numerical value o the magnetic susceptibility o the chiral condensate. Our results are as ollows: χ 4.3 GeV, NJL (6) χ 5.5 GeV, QM (7) respectively or the NJL model and the QM model. To obtain the numerical values above we have used data or eb up to 5m π 0.1 GeV, which are then it using a linear law. Using the numerical values o the chiral condensate in the two models, we obtain χ 69 MeV, NJL (8) χ 65 MeV, QM (9) The numerical values o χ that we obtain within the eective models are in air agreement with recent results, see Table I. In our model calculations, the role o the renormalization scale is played approximately by the ultraviolet cuto, that is Λ in Eq. (11), which is equal to GeV in the QM model, and 0.67 GeV in the NJL model. To acilitate the comparison with previous estimates, we review briely the rameworks in which the results in Table I are obtained. In [36] the ollowingresult is ound, within OPE combined with Pion Dominance (we ollow when possible the notation used in [41]): χ PD = c χ N c 8π F π, Pion Dominance (30) with F π = π = MeV and c χ = ; the estimate o [36] is done at a renormalization point M = 0.5 GeV. It is remarkable that Eq. (30) has been reproduced recently within AdS/QCD approach in [38]. Probably, this is the result more comparable to our estimate, because the reerence scales in [36] and in this article are very close. Withinourmodelcalculationsweindc NJL χ = 1.93 and c QM χ =.36. Using the numerical value o F π and c χ we get χ PD = 4.45 GeV, which agrees within 3% with our NJL model result, and within 18% with our QM model result. In [37] the authors ind c χ =.15 within hard-wall holographic approach, at the scale M 1 GeV. The results o [37] are thus in very good parametric agreement with [36]; on the other hand, the numerical value o F π in the holographic model is smaller than the one used in [36], pushing the holographic prediction or χ to slightly higher values than in [36]. However, the scale at which the result o [37] is valid should be much smaller than M = 1 GeV, thus some quantitative disagreement with [36] is expected. As the authors have explained, it might be possible to tune the parameters o the holographic model, mainly the chiral condensate, to reproduce the correct value o F π ; their numerical tests suggest that by changing the ratio /m ρ o a actor o 8, then the numerical value o c χ is inluenced only by a 5%. It is thereore plausible that a best tuning makes the quantitative prediction o [37] much closer to the estimate o [36]. In [39] an estimate o χ within the instanton vacuum model has been perormed beyond the chiral limit, both or light and or strange quarks (the result quoted in Table I corresponds to the light quarks; or the strange quark, χ s /χ u,d 0.15 is ound). Taking into account the numerical value o the chiral condensate in the instanton vacuum, the numerical estimate o [39] leads to χ =.5 ± 0.15 GeV at the scale M = 1 GeV. An analytic estimate within a similar ramework has been obtained in [40], in which the zero-mode o the Dirac operator in the background o a SU() instanton is used to compute the relevant expectation values. The result o [31] gives χ = 3.5 GeV at M 1 GeV. In [41] the result χ = GeV is achieved within a two-color simulation with quenched ermions. It is interesting that in [41] the same quantity has been computed also at inite temperature in the coninement phase, at T = 0.8T c, and the result seems to be independent on temperature. The reerence scale o [41], determined by the inverse lattice spacing, is M GeV. Thereore the lattice results are not quantitatively comparable with our model calculation. However, they share

8 8 an important eature with the results presented here, namely the saturation o the polarization at large values o the magnetic ield. Finally, estimates o the magnetic susceptibility o the chiral condensate by means o several QCD sum rules there exist [31 35]. The results are collected in Table I. IV. RENORMALIZED QM MODEL In this Section, we make semi-analytic estimates o the polarization and the magnetic susceptibility o the quark condensate, as well as or the chiral condensate in magnetic background, within the renormalized QM model. This is done with the scope to compare the predictions o the renormalized model with those o the eective models, in which an ultraviolet cuto is introduced to mimic the QCD eective quark mass. In the renormalized model, we allow the eective quark mass to be a constant in the whole range o momenta, which is dierent rom what happens in QCD [45]. Thus, the higher Landau levels give a inite contribution to the vacuum chiral condensate even at very strong ields. This is easy to understand: the ultraviolet cuto, Λ, in the renormalized model can be taken larger than any other mass scale, in particular Λ eb 1/ ; as a consequence, the condition p 3 + n eb < Λ is satisied taking into account many Landau levels even at very large eb. The contribution o the higher Landau levels, once renormalized, appears in the physical quantities to which we are interested here, in particular in the chiral condensate. Since the computation is a little bit lengthy, it is useul to anticipate its several steps: irstly we perorm regularization, and then renormalization, o the QEP at zero magnetic ield (the corrections due to the magnetic ield turn out to be ree o ultraviolet divergences). Secondly, we solve analytically the gap equation or the σ condensate in the limit o weak ields, and semi-analytically in the opposite limit. The ield-induced corrections to the QEP and to the solution o the gap equation are divergence-ree in agreement with [7], and are thereore independent on the renormalization scheme adopted. Then, we compute the renormalized and sel-consistent values o the chiral condensate and o the magnetic moment, as a unction o eb, using the results or the gap equation. Within this theoretical ramework, it is much more convenient to compute and Σ by taking derivatives o the renormalized potential; in act, the computation o the traces o the propagator in the renormalized model is much more involved i compared to the situation o the non-renormalized models, since in the ormer a non-perturbative(and non-trivial) renormalization procedure o composite local operators is required [54]. Finally, we estimate χ, as well as the behavior o the polarization as a unction o eb. A. Renormalization o the QEP To begin with, we need to regularize the one-loop ermion contribution in Eq. (9) namely V ermion 1 loop = N c n=0 Q eb π + dk ( β n k +n Q eb +m ) 1/ q. π (31) To this end, we deine the unction, V(s), o a complex variable, s, as V(s) = N c n=0 Q eb π + dk ( β n k +n Q eb +m q)1 s. π (3) The unction V(s) can be analytically continued to s = 0. We deine then V1 loop ermion = lim s 0 + V(s). Ater elementary integration over k, summation over n and taking the limit s 0 +, we obtain the result V1 loop ermion = N (Q eb) ( ) c 4π s log( Q eb )+a B (q) (Q eb) N c π ζ ( 1,q) N c Q eb m q 8π ( ) s log(m q )+a, (33) where we have subtracted terms which do not depend explicitly on the condensate. In the above equation, ζ(t,q) isthehurwitzzetaunction; orre(t) > 1andRe(q) > 0, it is deined by the series ζ(t,q) = n=0 (n +q) t ; the series can be analytically continued to a meromorphic unctiondeinedin thecomplexplanet 1. Moreoverwe have deined q = (m q + Q eb )/ Q eb ; urthermore, a = 1 γ E ψ( 1/),where γ E isthe Eulero-Mascheroni number and ψ is the digamma unction. The derivative ζ ( 1,q) = dζ(t,q)/dt is understood to be computed at t = 1. The irst two addenda in Eq. (33) arise rom the higher Landau levels; on the other hand, the last addendum is the contribution o the LLL. The unction B is the second Bernoulli polynomial; using its explicit orm, it is easy to show that the divergence in the LLL term in Eq. (33) is canceled by the analogous divergence in the irst addendum o the same equation. It is interesting that the LLL contribution, which is in principle divergent, combines with a part o the contribution o the higher Landau levels, leading to a inite result. This can be interpreted as a renormalization o the LLL contribution. On the other hand, the remaining part arising

9 9 Method χ (GeV ) Ren. Point (GeV) Re. Sum rules 8.6 ± [31] Sum rules [3] Sum rules 4.4 ± [33] Sum rules 3.15 ± [34] Sum rules.85 ± [35] OPE + Pion Dominance N c/(4π F π) 0.5 [36] Holography 1.075N c/(4π F π) 1 [37] Holography N c/(4π F π) 1 [38] Instanton vacuum.5 ± [39] Zero mode o Dirac Operator [40] Lattice 1.547(3) [41] NJL model This work QM model This work TABLE I. Magnetic susceptibility o the quark condensate obtained within several theoretical approaches. In the table, F π = MeV. See the text or more details. rom the higher Landau levels is still divergent; this divergence survives in the B 0 limit, and is due to the usual divergence o the vacuum contribution. We then have V1 loop ermion = N m 4 ( ) q c 16π s log( Q eb )+a +N c Q eb m q 8π log m q Q eb (Q eb) N c π ζ ( 1,q). (34) The divergence o the vacuum energy is made explicit in the expression in the square brackets in Eq. (34). The mathematical structure o the divergence, namely a pole in s = 0 and a logarithm with a dimensional argument, is similar to that obtained within the dimensional regularization scheme. The scale o the logarithm is hidden in the 1/s term, and appears explicitly when the divergence is subtracted. Such a divergence aects only the B = 0eective potential; the correctionsdue tothe magnetic ield are either inite or independent on the condensate. As a matter o act, in the zero magnetic ield limit q, the ermion bubble becomes m 4 ( q N c N 16π s logm q +a+ 1 ) V 0, (35) where we have used the relation ( ) 1 ζ ( 1,q) 1 q +log(q) B (q) 4, (36) with B corresponding to the second Bernoulli polynomial. Using Eqs. (35) and (34) we notice that the pole /s is cancelled in the dierence V 1 V1 loop ermion V 0, ( m 4 q V 1 = N c 16π + Q ) eb m q 8π log Q eb m q Q eb N c π ζ m 4 q ( 1,q) N c N 3π. (37) To remove the divergence o the vacuum term, we ollow Re. [7] and adopt the renormalization conditions that at zero magnetic ield, the quantum corrections do not shit the classical expectation value o the σ ield in the vacuum, σ = π, as well as the classical value m σ. To ulill these conditions, two counterterms have to be added to the eective potential, V c.t. = δλ σ 4 /4+δv σ /. This amounts to require (V 0 +V c.t. ) σ = (V 0 +V c.t. ) σ=π σ = 0. (38) σ=π Taking into account the one-loop divergent contribution and the condition (38), we can write the renormalized potential at zero ield as V = λ ( σ ṽ ) N c N m 4 q hσ 4 16π log m q g π, (39) where λ and ṽ are the renormalized parameters, whose expression is not needed here. Equation (39) is in agreement with the textbook result or the one-loop eective potential o a linear sigma model coupled to ermions by means o a Yuwaka-type interaction [5]. This is a consistency check o our regularization and renormalization procedures. B. Approximate solutions o the gap equation Weak ields. We now restore the magnetic ield, and

10 10 irstly take the weak ield limit. In this context, small ield means eb m q. Using again Eq. (36) we notice that the derivative o the ζ unction cancels the other addenda in Eq. (37), and the remaining contribution is V 1 N c = N c (Q eb) 4π (Q eb) 4π log m q Q eb log m q µ, (40) whichisin agreementwiththe resulto[7]. Inthe above equation we have ollowed the notation o [43] introducing an inrared scale µ, isolating and then subtracting the term which does not depend on the condensate. The scale µ is arbitrary, and we cannot determine it rom irst principles; on the other hand, it is irrelevant or the determination o the σ condensate. We expect µ π since this is the typical scale o chiral symmetry breaking in the model or the σ ield. The correction (40) lowers the eective potential o the broken phase, thus avoring the spontaneous breaking o chiral symmetry. The result in Eq. (40) is UV divergence-ree as anticipated; thus, it is independent on the particular procedure used to regularize the QEP. In this limit, it is easyto obtainanalyticallythe behavior o the constituent quark mass as a unction o eb. As amatteroact, wecanexpandthederivativeothe QEP with respect to σ, around the solution at B = 0, writing σ = π +δσ. Then, a straightorward evaluation leads to ( m q = g π 1+ 5 ) N c 91π πm (eb). (41) σ As anticipated, the scale µ is absent in the solution o the gap equation. Strong ields. In the limit eb m q, we can ind an asymptotic representation o V 1 by using the expansion ζ ( 1,q) = c 0 +c 1 (q 1) valid or q 1, with c 0 = 0.17 and c 1 = 0.4. Then we ind V 1 N c N c m q 8π ( m q + Q eb ) log Q eb m q Q eb m q π c 1, (4) where we have subtracted condensate-independent terms. In the strong ield limit it is not easy to ind analytically an asymptotic representation or the sigma condensate as a unction o eb; thereore we solve the gap equation numerically, and then it data with a convenient analytic orm as ollows: m q = b eb 1/ + c3 π eb, (43) where b = 0.3 and c = At large ields the quark mass grows as eb 1/ as expected by dimensional analysis; this is a check o the equations that we use. C. Evaluation o chiral condensate and magnetic moment Chiral condensate. To compute the chiral condensate we ollow a standard procedure: we introduce source term or, namely a bare quark mass m, then take derivative o the eective potential with respect to m evaluated at m = 0. This amounts to derive only, since both the classic potential and counterterms do not have a dependence on m. According to Eq. (37) we ind V1 loop ermion = V 0 m m =0 + V 1 m m =0. (44) The irst and second addenda on the r.h.s. o Eq. (44) represent the vacuum and the ield-induced contributions to the chiral condensate respectively. The vacuum contribution is divergent, and can be renormalized according to the procedure o renormalization o composite local operators [54]. It is not necessary to perorm this procedure here, because the numerical value o the vacuum condensate is not necessary in our discussion 1. Thereore, it is enough to compute only the contribution at B 0 arising rom V 1, which is inite. In particular, or the weak ield case we obtain = 0 N c 1π Q eb m q. (45) On the other hand, in the strong ield limit we have = N cm q ( Q 4π eb +m ) Q eb q log m q. (46) Using Equations (41) and (43), we show that the chiral condensate scales as a + b(eb) or small ields, and as eb 3/ or large ields. Magnetic moment. Next we turn to the computation o the expectation value o the magnetic moment. The expression in terms o Landau levels is given by Eq. (17), which clearly shows that this quantity has a log-type divergence. In order to avoid a complicated renormalization procedure o a local composite operator, we notice that it is enough to take the minus derivative o V 1 with respect to B to get magnetization, M [43], then multiply by m/q to get the magnetic moment. This procedure is very cheap, since the B dependent contributions to the eective potential are inite, and the resulting expectation value will turn out to be inite as well (that is, already renormalized). 1 In principle, in the renormalization program o the composite operator qq one can impose, as a renormalization condition, that the chiral condensate in the vacuum is consistent with the results o [49] at the renormalization point M = 1 GeV.

11 11 In the case o weak ields, rom Eq. (40) we ind Σ = N c Q eb m q 6π log m q µ. (47) On the other hand, in the strong ield limit we get rom Eq. (4) Σ = N c m 3 q 4π log Q eb m q. (48) Μ MeV The above result is in parametric agreement with the estimate o magnetization in [43]. In act, m q eb in the strong ield limit, which leads to a magnetization M BlogB. Using the expansions or the sigma condensate at small and large values o the magnetic ield strength, we argue that Σ eb in a weak ield, and Σ eb 3/ in a strong ield. D. Computation o chiral magnetization and polarization We can now estimate the magnetic susceptibility o the quark condensate and the polarization as a unction o eb. For the ormer, we need to know the behavior o the magnetic moment or weak ields. From Eq. (47) and rom the deinition (1) we read χ = N cm q 6π log m q (µ). (49) µ The presence o the inrared scale µ makes the numerical estimateoχuncertain; however,takingoritavalueµ π, whichisthetypicalscaleochiralsymmetrybreaking, we have χ 44 MeV, which is in agreement with the expected value, see Eq. (). In Fig. 4 we plot (µ) as a unction o µ. The interval on the µ-axis delimited by the green and the blue vertical lines is the range in which we obtain a value o χ which is consistent with phenomenology. Next we turn to the polarization. For weak ields we ind trivially a linear dependence o µ on Q eb, with slope given by the absolute value o χ in Eq. (49). On the other hand, in the strong ield limit we ind, according to Eq. (48), µ m q m q + Q eb 1 Q b+ Q, (50) where we have used Eq. (43). This result shows that the polarization saturates at large values o eb, but the asymptotic value depends on the lavor charge. It is interesting to compare the result o the renormalized model with that o the eective models considered in the previous Section. In the ormer, the asymptotic value o µ is lavor-dependent; in the latter, µ 1 independently on the value o the electric charge. Our Μ Π FIG. 4. Magnetic susceptibility o the quark condensate multiplied by the chiral condensate, in MeV, as a unction o the inrared scale µ, as given by Eq. (49). The interval on the µ-axis delimited by the green and the blue vertical lines is the range in which we obtain a value o χ which is consistent with phenomenology. interpretation o this dierence is as ollows: comparing Eq.(50) with the general model expectation, Eq.(), we recognize in the actor Q /(b + Q ) the contribution o the higher Landau levels at zero temperature, which turns out to be inite and non-zero ater the renormalization procedure. This contribution is then transmitted to the physical quantities that we have computed. The trace o the higher Landau levels is implicit in the solution o the gap equation in the strong ield limit, namely the actor b in Eq. (43), and explicit in the additional Q dependence in Eq. (50). A posteriori, this conclusion seems quite natural, because in the renormalization procedure we assume that the eective quark mass is independent on quark momentum, thus there is no cut o the large momenta in the gap equation (and in the equation or polarization as well). In the eective models considered in the irst part o this article, on the other hand, the cuto procedure is equivalent to have a momentum-dependent eective quark mass, m q = gσθ(λ p 3 n Q eb ), which naturally cuts o higher Landau levels when eb Λ. At the end o the days, the expulsion o the higher Landau levels rom the chiral condensate makes µ 1 in the strong ield limit. Our expectation is that i we allow the quark mass to run with momentum and decay rapidly at large momenta, mimicking the eective quark mass o QCD, higher Landau levels would be suppressed in the strong ield limit, and the result (50) would tend to the result in Fig.. V. CONCLUSIONS In this article, we have computed the magnetic susceptibility o the quark condensate, χ, and the polarization, µ, in a background o a magnetic ield, B, by means o

12 1 two chiral models o QCD: the quark-meson model and the Nambu-Jona-Lasinio model. The knowledge o these quantities is relevant both theoretically and phenomenologically. Indeed, the magnetic susceptibility o the quark condensate might lead, in the photoproduction o lepton pairs, to an intererence between the Drell-Yan process and the photon-quark coupling, the latter induced by the chiral-odd distrubution amplitude o the quark [4]. The two models are widely used to study the phase diagram o QCD in several regimes; it has been proved in dierent contexts, that they oer a good theoretical tool to compute low-energy QCD properties. It is thus interesting to compute, within these models, quantities which characterize the QCD vacuum in a magnetic background, and compare the results with those obtained within dierent theoretical rameworks, both analytically and numerically. In the irst part o this article, we have reported our results about χ and µ obtained within a numerical selconsistent solution o the model, in the one-loop approximation. Our results on χ are summarized in Table I, and are in air agreement with previous estimates. Besides, our data on polarization are collected in Fig. 3. For the latter, we obtain a saturation to the asymptotic value µ = 1 at large values o eb, which is understood as a lowest Landau level dominance in the quark condensate. This is in agreement with recent results obtained within Lattice QCD simulations with two colors and quenched ermions [41]. In the second part o the article, we have estimated χ and µ within a renormalized version o the QM model. The numerical value o χ in this case is quite uncertain becauseothe presence,inthe inalresult, oanunknown inrared scale µ. However, as shown in Fig. 4, taking or µ a numerical value around π, which is the typical scale o chiral symmetry breaking in this model, we obtain a value o χ which is consistent with phenomenology. Within the renormalized model we ind a saturation o µ at large B, in qualitative agreement with our indings within the eective models in Section III and quenched QCD [41]. In the renormalized model, the asymptotic value o µ is lavor-dependent; in the eective models and quenched QCD, µ 1 independently on the value o the electric charge. We attribute this dierence to the presence, in the theory, o the renormalized contribution o the higher Landau levels: in the renormalized model, the higher Landau levels give a inite contribution to the chiral condensate at zero temperature, even in the case o very strong ields. Ater renormalization o the QEP, this contribution is inite and non-zero, and is transmitted to the physical quantities that we have computed. The results we obtained are quite encouraging, and suggest that a systematic study o external ield induced condensates at zero, as well as inite, temperature within chiral models is worth to be done. As a natural continuation o this study, it would be interesting to add the strange quark, ollowing the work[39]. Moreover, Lattice simulations have shown that χ is almost insensitive to the temperature, at least in the coninement phase [41]. This result is achieved with quenched ermions. On the other hand, the chiral model handle with dynamical ermions; it would be thereore o interest to address the question o the behavior o χ as a unction o temperature using chiral models. Furthermore, it has been shown [53] that local luctuations o topological charge induce a quark electric dipole moment along the direction o a strong magnetic ield. This problem can be studied easily within the chiral models, introducing a pseudo-chemical potential µ 5 conjugated to chirality imbalance as already done in [4, 30]. Acknowledgements. We acknowledge M. Chernodub, T. Cohen, A. Flachi, R. Gatto, A. Ohnishi or several discussions and correspondence. Moreover, we thank H. Suganuma or valuable and numerous discussions on the renormalization o the quark-meson model and or a careul reading o the manuscript. The Laboratoire de Matematique et Physique Teorique o Tours University, as well as the Physics Department o Jyvaskyla University, where part o this work was completed, are acknowledged or the kind hospitality. The work o M. R. is supported by JSPS under contract number P0908. [1] P.deForcrandandO.Philipsen, JHEP0701, 077(007); JHEP 0811, 01 (008); PoS LATTICE008, 08 (008). [] Y. Aoki et al., Nature 443, 675 (006); Y. Aoki, S. Borsanyi, S. Durr, Z. Fodor, S. D. Katz, S. Krieg and K. K. Szabo, JHEP 0906, 088 (009); S. Borsanyi et al., S. Borsanyi et al., JHEP 1009 (010) 073; [3] A. Bazavov et al., Phys. Rev. D 80, (009). [4] M. Cheng et al., Phys. Rev. D 81, (010). [5] F. Karsch, E. Laermann and A. Peikert, Nucl. Phys. B 605, 579 (001); F. Karsch, Lect. Notes Phys. 583, 09 (00); O. Kaczmarek and F. Zantow, Phys. Rev. D 71, (005). [6] C. S. Fischer, Phys. Rev. Lett. 103, (009); C. S. Fischer and J. A. Mueller, Phys. Rev. D 80, (009); C. S. Fischer, A. Maas and J. A. Muller; Eur. Phys. J. C 68, 165 (010). A. C. Aguilar and J. Papavassiliou, arxiv: [hep-ph]; C. Feuchter and H. Reinhardt, Phys. Rev. D 70, (004); Phys. Rev. D 71, (005). M. Leder, J. M. Pawlowski, H. Reinhardt and A. Weber, arxiv: [hep-th]; J. Braun, Eur. Phys. J. C 64, 459 (009); J. Braun, L. M. Haas, F. Marhauser and J. M. Pawlowski, arxiv: [hep-ph]; J. Braun and A. Janot, arxiv: [hep-ph]. [7] M. Frasca, Phys. Lett. B 670, 73 (008); Mod. Phys. Lett. A 4 (009) 45; arxiv: [hep-ph]. [8] M. A.Shiman, A.I.VainshteinandV.I.Zakharov, Nucl. Phys. B 147, 385 (1979); Nucl. Phys. B 147, 448 (1979). [9] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 1, 345