# A search for inverse magnetic catalysis in thermal quark-meson models

## Abstract

We explore the parameter space of the two-flavor thermal quark-meson model and its Polyakov loop-extended version under the influence of a constant external magnetic field . We investigate the behavior of the pseudo critical temperature for chiral symmetry breaking taking into account the likely dependence of two parameters on the magnetic field: the Yukawa quark-meson coupling and the parameter of the Polyakov loop potential. Under the constraints that magnetic catalysis is realized at zero temperature and the chiral transition at is a crossover, we find that the quark-meson model leads to thermal magnetic catalysis for the whole allowed parameter space, in contrast to the present picture stemming from lattice QCD.

## I Introduction

The phase diagram of magnetic quantum chromodynamics QCD, i.e. for strong interactions under the influence of an external classical Abelian magnetic field, is currently under construction. This phase diagram corresponds to a special case, as it does not suffer from the Sign Problem, and can be easily simulated on the lattice. So, from the theoretical point of view, it serves as a crucial check for effective models of QCD extended to regions not easily accessible by lattice simulations. From the experimental standpoint, this setup is also quite remarkable owing to the fact that strong magnetic fields are relevant in non-central heavy ion collisions, and play a major role in the possibility of observing the chiral magnetic effect (for a comprehensive review, see Ref. Kharzeev:2013jha ().)

The chiral and deconfining transitions under the effect of a magnetic background are, of course, amenable also to effective model descriptions Fraga:2012rr (); Gatto:2012sp (). Those models have predicted several outstanding new features to the thermodynamics of strong interactions, from shifting the chiral and the deconfinement crossover lines in the phase diagram Agasian:2008tb (); Fraga:2008qn (); Menezes:2008qt (); Boomsma:2009yk (); Mizher:2010zb (); Fukushima:2010fe (); Gatto:2010pt (); Kashiwa:2011js (); Chatterjee:2011ry (); Andersen:2011ip (); Skokov:2011ib (); Andersen:2012bq (); Fukushima:2012xw (); Fraga:2012fs (); Strickland:2012vu (); Fukushima:2012kc (); Ferrari:2012yw (); Allen:2013eha (); Fu:2013ica (); Andersen:2012jf (); Andersen:2012zc (); Avancini:2012ee (); Andersen:2012dz (); Costa:2013zca (); Allen:2013lda (); Ballon-Bayona:2013cta (); Callebaut:2013ria (); Ferreira:2013tba (); Ruggieri:2013cya () to transforming the vacuum into a superconducting medium via -meson condensation Chernodub:2010qx (); Chernodub:2011mc (), for high enough magnetic fields, i.e. a few times . Nevertheless, the available lattice data D’Elia:2010nq (); D’Elia:2011zu (); Bali:2011qj (); Bali:2012zg (); Bruckmann:2013oba () contradicted essentially all predictions regarding the behavior of the pseudo critical lines for deconfinement and chiral symmetry restoration coming from chiral models (including their Polyakov loop extensions). The reason for this failure is unclear, but the fact that confinement is not properly captured in such chiral models might play a role Fraga:2012fs (); Fraga:2012ev (), as well as an apparently nontrivial chiral limit Blaizot:2012sd ().

In this paper, we explore the parameter space of the two-flavor thermal quark-meson (QM) model and its Polyakov loop-extended version (PQM) under the influence of a constant external magnetic field . We investigate the behavior of the pseudo critical temperature for chiral symmetry breaking taking into account the likely dependence of two parameters on the magnetic field: the Yukawa quark-meson coupling and the parameter of the Polyakov loop potential. We scan an important part of the parameter spaces of these models, in order to check whether they can accommodate, at least qualitatively, the trend of inverse magnetic catalysis found in Ref. Bali:2011qj (). In doing so, we keep two important constraints: (i) that magnetic catalysis is realized at zero temperature, and (ii) that the chiral transition at is a crossover. In these two limits, i.e. zero temperature or zero magnetic field, chiral effective models usually produce results that are, at least qualitatively, in line with lattice QCD. We find, nevertheless, that the extensions considered by introducing a dependence in the Yukawa coupling and in the parameter are not enough to account for the behavior of the critical temperature as found in lattice QCD simulations. We believe this should be also the case for other chiral models, which would signal the lack of some fundamental ingredient, possibly quark confinement.

The paper is organized as follows. In Section II we describe the thermal effective potential in the presence of a magnetic background. In Section III we consider the effects from the running of the Yukawa coupling, having the magnetic field as the momentum scale, in the QM model. In Section IV we investigate the consequences of making dependent, in a fashion usually made with chemical potentials. Section V contains our summary.

## Ii Thermal effective potential in the presence of a magnetic field

Let us first briefly review the Polyakov-Quark-Meson (PQM) model while we introduce our notation. The effective potential of the PQM model at finite temperature and magnetic field was calculated in Mizher:2010zb (). At the mean-field level, it is a sum of four contributions,

(1) |

where

(2) |

is the tree-level potential for the field. The Polyakov loop potential (in the logarithmic parametrization) is given by Roessner_08 ()

(3) |

The contribution to the effective potential (with the normalization ) has the form

(4) |

whereas is given by

(5) |

and represents the paramagnetic (thermal) contribution to the effective potential.

In the previous formulas, is the magnetic field, is the temperature, is the expectation value of the sigma meson field (the approximate order parameter for the chiral transition), , , and are the (approximate) deconfinement order parameters related to the Polyakov loop (see, e.g., Mizher:2010zb () for more details). The quantities , and are parameters of the mesonic self interaction (2), which are adjusted according to the constants on Table 1, and the functions and are defined as

(6) |

where the values of the parameters , , , and are listed in Table 2. In eq. (4), we used , where is the meson-quark coupling and is the charge of the quark of flavor .

Constant | [MeV] | [MeV] | [MeV] | ||||
---|---|---|---|---|---|---|---|

Value | 93 | 138 | 20 | 3.3 | 307 |

Constant | [MeV] | ||||
---|---|---|---|---|---|

Value |

Notice that the Quark-Meson (QM) model is recovered, as a particular case of the PQM model, when one sets (or, equivalently, ) and neglects the Polyakov loop potential (). In the PQM model, it can be shown that the reality of the minimum of the effective potential implies , so that there are actually only two order parameters in the model, and .

The state of thermodynamic equilibrium is the minimum of the effective potential with respect to the order parameters and at fixed external conditions ( and ). One defines, then, the chiral condensate , as well as the chiral susceptibility . The pseudo critical temperature for chiral symmetry breaking is identified as the peak of the chiral susceptibility . As discussed in Ref. Mizher:2010zb (), the pseudo critical temperature for the chiral and the deconfinement transitions split as the magnetic field increases sufficiently, a fact that is not in line with lattice results Bali:2011qj (). To keep our discussion simpler, we will focus our analysis on the chiral transition.

In the next two sections we promote two parameters to functions of the external magnetic field: the Yukawa coupling and the parameter of the Polyakov loop potential.

## Iii B-dependent Yukawa coupling in the QM model

Let us first consider the modifications brought about by introducing a dependence in the Yukawa coupling, i.e. by considering . For that purpose, we adopt the QM model (which is equivalent to setting and in the PQM model). This will change our formulas (4) and (5) for the effective potential simply by the replacement . In order to describe the parameter space of , we show the pseudo critical temperature as a function of the Yukawa coupling for different values of the applied field , as seen in Fig. 1.

Fig. 1 should be understood in the following way. Each continuous function will correspond to some continuous path on the plot. All these paths must start on the continuous line (), then proceed to some point on the dashed line (), then to the dotted line (), and so on, as increases. It is clear from Fig. 1 that any choice (a constant) leads to an increasing , i.e. magnetic catalysis. The same behavior will set in if is a decreasing function of . The only possible way to find a decreasing would be to let increase with .

However, as can be seen from Fig. 1, such a function would have to grow very strongly. Thus, the coupling would very soon reach the region , where the transition becomes of first order. Lattice results Bali:2011qj () clearly show that the transition is a crossover for finite and, therefore, we do not consider first-order transitions as acceptable solutions in this paper. Besides, if , as required by the parameter fit in the vacuum ( and ), there is clearly no continuous that would give rise to a decreasing (at least not for a crossover, see Fig. 1). Therefore, we can conclude that, if one takes the usual parameter fixing in the vacuum, , there is no continuous function that could lead to inverse magnetic catalysis in the QM model at finite temperature and zero quark chemical potential, unless the chiral transition is of first order.

To illustrate the general trend discussed above, we consider two explicit examples. We show in Fig. 2 the evolution of the chiral order parameter as a function of the temperature for different values of . We take the (increasing) ansatz . For small values of the applied field, the transition is a crossover. However, as soon as reaches a moderate value (), the chiral crossover turns into a first-order transition.

## Iv B-dependent in the PQM model

Let us now consider the PQM model with the standard (fixed) Yukawa coupling , but with a varying . As discussed in Ref. Schaefer_07 (); Roessner_08 (), the model parameter is equal to the deconfinement critical temperature in a pure gauge Polyakov loop model, MeV. However, due to the back reaction from the matter fields, should be lower than its pure gauge value when dynamical quarks are present. More specifically, one takes as the running scale in the perturbative renormalization group running of , the strong coupling constant. Once depends also on the number of quark flavors, quark masses, and possibly on external constraints (e.g. on the chemical potential, as in Ref. Schaefer_07 (), or on the magnetic field, as considered here), a function may be implemented, where is a given energy scale used to fix (e.g., the tau lepton mass GeV PDG_12 ()).

According to the analysis of Ref. Bruckmann:2013oba (); Bruckmann:2013ufa (), the form of the Polyakov loop effective potential could be important for the mechanism of inverse magnetic catalysis around the transition temperature. A natural way to parametrize the dependence of the Polyakov loop potential is through the parameter, described above. We investigate here whether or not this possibility can lead to inverse magnetic catalysis. We find that, within the PQM model (with all other parameters kept at their vacuum values), no function can lead to a decreasing for sufficiently high magnetic fields, such as those investigated in Bali:2011qj (). Some choices for may, at most, make the pseudo critical temperature decrease for small values of . However, even for magnetic field intensities , the catalysis-inducing vacuum term (4) dominates the whole picture and inevitably rises.

In order to see that the pseudo critical temperature tends to rise for every parametrization , we proceed in the same fashion as in the analysis of in the previous section. Fig. 3 shows, for different intensities of magnetic field, how the pseudo critical temperature depends on the value of . For reference, we also show the value MeV, obtained in the two-flavor lattice calculation of Karsch:2000kv (). In this diagram, each continuous function corresponds to some line connecting the curves shown. One can immediately realize that the pseudo critical temperature may go down only if decreases sufficiently fast. Such behavior, however, can not be sustained even at moderate intensities of the magnetic field, regardless of how fast may decrease as a function of .

Let us provide explicit examples. The choices
and , where ,
lead to the functions shown in Fig. 4(a).
The resulting normalized pseudo critical temperatures are
shown^{1}

## V Conclusions

We examined the behavior of the pseudo critical temperature for the chiral transition in the QM and in the PQM model in the presence of a constant magnetic field background by allowing for two natural extensions of the parameter space: a -dependent Yukawa coupling, , in the quark-meson model and a -dependent Polyakov loop potential parameter, , in the PQM model.

We found that even strongly varying functions and are not able to reproduce (not even qualitatively) the phase diagram from lattice simulations in the presence of a magnetic field up to fields GeV. This robustness of the (P)QM model can be attributed to the vacuum term, Eq. (4), which is responsible for magnetic catalysis in the vacuum and still dominates the effective potential at finite temperature, even for moderate values of the magnetic field intensity.

Even though our analysis has been performed at mean-field level, we believe that our conclusions will not be changed if one considers fluctuations (as studied, e.g., in Refs. Fukushima:2012xw () and Andersen:2013swa () for fixed values of the model parameters). Likewise, there is no reason to believe that other chiral models, such as the PNJL model, will behave differently.

The physical mechanism for inverse magnetic catalysis recently proposed in Bruckmann:2013oba (); Bruckmann:2013ufa () relies on a competition between valence and sea quarks. While the coupling of the magnetic field to valence quarks enhances magnetic catalysis, its interaction with sea quarks tends to decrease the chiral condensate, leading to inverse magnetic catalysis. In the context of chiral models, however, quarks are introduced as free quasiparticles, so that no degrees of freedom corresponding to sea quarks are present. We believe this can be an important, if not decisive, factor for the inability of chiral models to correctly describe the phase diagram.

It seems clear that further comparisons between predictions from these effective models and lattice data in systems containing other control parameters besides the temperature are crucial to test the reliability of those models. Besides the case with a magnetic field, these models seem to face difficulties describing lattice data also when one considers the dependence of the critical temperature on isospin and quark masses Dumitru:2003cf (); Fraga:2008be (); Stiele:2013pma (). Apparently, the presently used chiral models miss some important physics to describe the properties of bulk QCD matter as extracted from lattice simulations when introducing new control parameters, here the magnetic field, even on a qualitative level. In view of this circumstance, it seems that extrapolations of these chiral models to regions of the QCD phase diagram not easily accessible by lattice simulations have to be taken with caution before these basic discrepancies are resolved.

###### Acknowledgements.

The authors thank Rainer Stiele for many valuable discussions during the initial stage of this work. The lattice data points of Fig. 4(b) were kindly provided by Gergely Endrődi. BWM thanks the Institut für Theoretische Physik of the Goethe Universität, Frankfurt, for their kind hospitality during a visit in which this project was initiated. The work of ESF was financially supported by the Helmholtz International Center for FAIR within the framework of the LOEWE program (Landesoffensive zur Entwicklung Wissenschaftlich-Ökonomischer Exzellenz) launched by the State of Hesse. The work of JSB was supported by the Hessian LOEWE initiative through the Helmholtz International Center for FAIR (HIC for FAIR) and the Alliance Program of the Helmholtz Association (HA216/EMMI).### Footnotes

- Our calculation ceases at GeV due to numerical instabilities which arise when is sufficiently small. Such instabilities are related to the Polyakov loop potential (3), which is ill-defined at .

### References

- D. Kharzeev, K. Landsteiner, A. Schmitt, and H.-U. Yee, “Strongly Interacting Matter in Magnetic Fields,” Lect.Notes Phys. 871 (2013) 1–624.
- E. S. Fraga, “Thermal chiral and deconfining transitions in the presence of a magnetic background,” Lect.Notes Phys. 871 (2013) 121–141, arXiv:1208.0917 [hep-ph].
- R. Gatto and M. Ruggieri, “Quark Matter in a Strong Magnetic Background,” Lect.Notes Phys. 871 (2013) 87–119, arXiv:1207.3190 [hep-ph].
- N. Agasian and S. Fedorov, “Quark-hadron phase transition in a magnetic field,” Phys.Lett. B663 (2008) 445–449, arXiv:0803.3156 [hep-ph].
- E. S. Fraga and A. J. Mizher, “Chiral transition in a strong magnetic background,” Phys.Rev. D78 (2008) 025016, arXiv:0804.1452 [hep-ph].
- D. Menezes, M. Benghi Pinto, S. Avancini, A. Perez Martinez, and C. Providencia, “Quark matter under strong magnetic fields in the Nambu-Jona-Lasinio Model,” Phys.Rev. C79 (2009) 035807, arXiv:0811.3361 [nucl-th].
- J. K. Boomsma and D. Boer, “The Influence of strong magnetic fields and instantons on the phase structure of the two-flavor NJL model,” Phys.Rev. D81 (2010) 074005, arXiv:0911.2164 [hep-ph].
- A. J. Mizher, M. Chernodub, and E. S. Fraga, “Phase diagram of hot QCD in an external magnetic field: possible splitting of deconfinement and chiral transitions,” Phys.Rev. D82 (2010) 105016, arXiv:1004.2712 [hep-ph].
- K. Fukushima, M. Ruggieri, and R. Gatto, “Chiral magnetic effect in the PNJL model,” Phys.Rev. D81 (2010) 114031, arXiv:1003.0047 [hep-ph].
- R. Gatto and M. Ruggieri, “Deconfinement and Chiral Symmetry Restoration in a Strong Magnetic Background,” Phys.Rev. D83 (2011) 034016, arXiv:1012.1291 [hep-ph].
- K. Kashiwa, “Entanglement between chiral and deconfinement transitions under strong uniform magnetic background field,” Phys.Rev. D83 (2011) 117901, arXiv:1104.5167 [hep-ph].
- B. Chatterjee, H. Mishra, and A. Mishra, “Vacuum structure and chiral symmetry breaking in strong magnetic fields for hot and dense quark matter,” Phys.Rev. D84 (2011) 014016, arXiv:1101.0498 [hep-ph].
- J. O. Andersen and R. Khan, “Chiral transition in a magnetic field and at finite baryon density,” Phys.Rev. D85 (2012) 065026, arXiv:1105.1290 [hep-ph].
- V. Skokov, “Phase diagram in an external magnetic field beyond a mean-field approximation,” Phys.Rev. D85 (2012) 034026, arXiv:1112.5137 [hep-ph].
- J. O. Andersen and A. Tranberg, “The Chiral transition in a magnetic background: Finite density effects and the functional renormalization group,” JHEP 1208 (2012) 002, arXiv:1204.3360 [hep-ph].
- K. Fukushima and J. M. Pawlowski, “Magnetic catalysis in hot and dense quark matter and quantum fluctuations,” Phys.Rev. D86 (2012) 076013, arXiv:1203.4330 [hep-ph].
- E. S. Fraga and L. F. Palhares, “Deconfinement in the presence of a strong magnetic background: an exercise within the MIT bag model,” Phys.Rev. D86 (2012) 016008, arXiv:1201.5881 [hep-ph].
- M. Strickland, V. Dexheimer, and D. Menezes, “Bulk Properties of a Fermi Gas in a Magnetic Field,” Phys.Rev. D86 (2012) 125032, arXiv:1209.3276 [nucl-th].
- K. Fukushima and Y. Hidaka, “Magnetic Catalysis vs Magnetic Inhibition,” Phys.Rev.Lett. 110 (2013) 031601, arXiv:1209.1319 [hep-ph].
- G. N. Ferrari, A. F. Garcia, and M. B. Pinto, “Chiral Transition Within Effective Quark Models Under Magnetic Fields,” Phys.Rev. D86 (2012) 096005, arXiv:1207.3714 [hep-ph].
- P. G. Allen and N. N. Scoccola, “Phase diagram of strongly interacting matter under strong magnetic fields,” arXiv:1307.4070 [hep-ph].
- W.-j. Fu, “Fluctuations and correlations of hot QCD matter in an external magnetic field,” arXiv:1306.5804 [hep-ph].
- J. O. Andersen and A. A. Cruz, “Two-color QCD in a strong magnetic field: The role of the Polyakov loop,” Phys.Rev. D88 (2013) 025016, arXiv:1211.7293 [hep-ph].
- J. O. Andersen, “Chiral perturbation theory in a magnetic background - finite-temperature effects,” JHEP 1210 (2012) 005, arXiv:1205.6978 [hep-ph].
- S. S. Avancini, D. P. Menezes, M. B. Pinto, and C. Providencia, “The QCD Critical End Point Under Strong Magnetic Fields,” Phys.Rev. D85 (2012) 091901, arXiv:1202.5641 [hep-ph].
- J. O. Andersen, “Thermal pions in a magnetic background,” Phys.Rev. D86 (2012) 025020, arXiv:1202.2051 [hep-ph].
- P. Costa, M. Ferreira, H. Hansen, D. P. Menezes, and C. ProvidÃªncia, “Phase transition and CEP driven by an external magnetic field in asymmetric quark matter,” arXiv:1307.7894 [hep-ph].
- P. G. Allen and N. N. Scoccola, “Quark matter under strong magnetic fields in SU(2) NJL-type models: parameter dependence of the cold dense matter phase diagram,” arXiv:1309.2258 [hep-ph].
- A. Ballon-Bayona, “Holographic deconfinement transition in the presence of a magnetic field,” arXiv:1307.6498 [hep-th].
- N. Callebaut and D. Dudal, “On the transition temperature(s) of magnetized two-flavour holographic QCD,” Phys.Rev. D87 (2013) 106002, arXiv:1303.5674 [hep-th].
- M. Ferreira, P. Costa, D. P. Menezes, C. ProvidÃªncia, and N. Scoccola, “Deconfinement and chiral restoration within the SU(3) PNJL and EPNJL models in an external magnetic field,” arXiv:1305.4751 [hep-ph].
- M. Ruggieri, M. Tachibana, and V. Greco, “Renormalized vs Nonrenormalized Chiral Transition in a Magnetic Background,” arXiv:1305.0137 [hep-ph].
- M. Chernodub, “Superconductivity of QCD vacuum in strong magnetic field,” Phys.Rev. D82 (2010) 085011, arXiv:1008.1055 [hep-ph].
- M. Chernodub, “Spontaneous electromagnetic superconductivity of vacuum in strong magnetic field: evidence from the Nambu–Jona-Lasinio model,” Phys.Rev.Lett. 106 (2011) 142003, arXiv:1101.0117 [hep-ph].
- M. D’Elia, S. Mukherjee, and F. Sanfilippo, “QCD Phase Transition in a Strong Magnetic Background,” Phys.Rev. D82 (2010) 051501, arXiv:1005.5365 [hep-lat].
- M. D’Elia and F. Negro, “Chiral Properties of Strong Interactions in a Magnetic Background,” Phys.Rev. D83 (2011) 114028, arXiv:1103.2080 [hep-lat].
- G. Bali, F. Bruckmann, G. Endrodi, Z. Fodor, S. Katz, et al., “The QCD phase diagram for external magnetic fields,” JHEP 1202 (2012) 044, arXiv:1111.4956 [hep-lat].
- G. Bali, F. Bruckmann, G. Endrodi, Z. Fodor, S. Katz, et al., “QCD quark condensate in external magnetic fields,” Phys.Rev. D86 (2012) 071502, arXiv:1206.4205 [hep-lat].
- F. Bruckmann, G. Endrodi, and T. G. Kovacs, “Inverse magnetic catalysis and the Polyakov loop,” JHEP 1304 (2013) 112, arXiv:1303.3972 [hep-lat].
- E. S. Fraga, J. Noronha, and L. F. Palhares, “Large Nc Deconfinement Transition in the Presence of a Magnetic Field,” Phys.Rev. D87 (2013) 114014, arXiv:1207.7094 [hep-ph].
- J.-P. Blaizot, E. S. Fraga, and L. F. Palhares, “Effect of quark masses on the QCD presssure in a strong magnetic background,” Phys.Lett. B722 (2013) 167–171, arXiv:1211.6412 [hep-ph].
- S. Rößner, T. Hell, C. Ratti, and W. Weise, “The chiral and deconfinement crossover transitions: PNJL model beyond mean field,” Nucl. Phys. A 814 (2008) 118–143, arXiv:0712.3152 [hep-ph].
- J. Beringer and Particle Data Group, “Review of Particle Physics,” Phys. Rev. D 86 no. 1, (2012) 010001.
- B. Schaefer, J. M. Pawlowski, and J. Wambach, “Phase structure of the Polyakov-quark-meson model,” Phys. Rev. D 76 no. 7, (2007) 074023–+, arXiv:0704.3234 [hep-ph].
- F. Karsch, E. Laermann, and A. Peikert, “Quark mass and flavor dependence of the QCD phase transition,” Nucl.Phys. B605 (2001) 579–599, arXiv:hep-lat/0012023 [hep-lat].
- F. Bruckmann, G. Endrodi, and T. G. Kovacs, “Inverse magnetic catalysis in QCD,” arXiv:1311.3178 [hep-lat].
- J. O. Andersen, W. R. Naylor, and A. Tranberg, “Chiral and deconfinement transitions in a magnetic background using the functional renormalization group with the Polyakov loop,” arXiv:1311.2093 [hep-ph].
- A. Dumitru, D. Roder, and J. Ruppert, “The Quark mass dependence of T(c) in QCD: Working up from m = 0 or down from m = infinity?,” Phys.Rev. D70 (2004) 074001, arXiv:hep-ph/0311119 [hep-ph].
- E. Fraga, L. Palhares, and C. Villavicencio, “Quark mass and isospin dependence of the deconfining critical temperature,” Phys.Rev. D79 (2009) 014021, arXiv:0810.1060 [hep-ph].
- R. Stiele, E. S. Fraga, and J. Schaffner-Bielich, “Thermodynamics of (2+1)-flavor strongly interacting matter at nonzero isospin,” arXiv:1307.2851 [hep-ph].