QCD phase diagram in a constant magnetic background. Inverse magnetic catalysis: where models meet the lattice

Magnetic catalysis is the enhancement of a condensate due to the presence of an external magnetic field. Magnetic catalysis at $T=0$ is a robust phenomenon in low-energy theories and models of QCD as well as in lattice simulations. We review the underlying physics of magnetic catalysis from both perspectives. The quark-meson model is used as a specific example of a model that exhibits magnetic catalysis. Regularization and renormalization are discussed and we pay particular attention to a consistent and correct determination of the parameters of the Lagrangian using the on-shell renormalization scheme. A straightforward application of the quark-meson model and the NJL model leads to the prediction that the chiral transition temperature $T_{\chi}$ is increasing as a function of the magnetic field $B$. This is in disagreement with lattice results, which show that $T_{\chi}$ is a decreasing function of $B$, independent of the pion mass. The behavior can be understood in terms of the so-called valence and sea contributions to the quark condensate and the competition between them. We critically examine these ideas as well recent attempts to improve low-energy models using lattice input.


Introduction
The phase diagram of QCD has received a lot of attention since the first ideas appeared in the 1970s. At that time, it was thought that QCD has two phases, a hadronic phase at low temperatures and a deconfined phase of quarks and gluons at high temperatures. In 1984, Bailin and Love [1] suggested that at high density, quark matter should be a color superconductor. The ideas are analogous to those of ordinary superconductivity and BCS theory [2], namely the instability of the Fermi surface to form Cooper pairs under an attractive interaction. In QCD, an attractive interaction a e-mail: andersen@tf.phys.ntnu.no is provided by one-gluon exchange in the triplet channel. Since then, there have been large efforts to map out the phase diagram of QCD and study the properties of its different phases [3][4][5]. The phase diagram has shown to be surprisingly rich at high baryon density and low temperatures. It includes a quarkyonic phase [6] as well as a number of superconducting phases, some of them being inhomogeneous. Most of these results have been obtained using low-energy models of QCD, notably the quark-meson (QM) model and the Nambu-Jona-Lasinio (NJL) model, with or without coupling to the Polyakov loop. The reason is that lattice simulations are notoriously difficult to perform at finite baryon chemical potential µ B due to the sign problem, so that one cannot use techniques involving importance sampling.
The temperature T and baryon chemical potential µ B are not the only relevant parameters of QCD. For example, one can introduce a separate chemical potential µ f for each quark flavor f . For two flavors, this leads to another independent chemical potential besides µ B , namely the isospin chemical potential µ I . For the three flavors, µ S = 1 2 (µ u + µ d − 2µ s ) is added. The addition of these chemical potentials gives rise to pion and kaon condensation. At T = 0, pion condensation occurs for µ I > m π , while kaon condensation takes place for | ± 1 2 µ I + µ S | > m K (upper sign for charged kaons and lower sign for neutral kaons). The former is particularly interesting since finite µ I and vanishing µ B has no sign problem and is therefore amenable to lattice simulations.
The final example of an external parameter, which is the topic of this review, is a (constant) magnetic background. There are several areas of high-energy physics, where such a background is relevant. One is non-central heavy-ion collisions, where large, time-dependent fields are generated. These fields are short-lived and have a maximum value of approximately eB = 6m 2 π [7]. The basis mechanism is simply that (in the center-of mass frame) the two nuclei repre-sent electric currents that according to Maxwell's equations generate a magnetic field. Another example where strong magnetic fields appear, are magnetars [8]. This is a special class of neutron stars with relatively low rotation frequencies. It is believed that the magnetic fields on the surface are 10 14 − 10 15 Gauss, while in the interior they can be as strong as 10 16 − 10 19 Gauss.
We consider QCD with an SU(3) gauge group, a global SU(N f ) vector symmetry and quark masses m f . The QCD Lagrangian is where the gluon field strength tensor is F a µν = ∂ µ A a ν −∂ ν A a µ − g f abc A b µ A c ν , f abc are the structure constants the covariant derivative in the presence of an abelian background field Moreover, m f is the mass of a quark of flavor f and there is a sum of flavors in Eq. (1). The nonabelian gauge field is A µ = t a A a µ , t a = 1 2 λ a , and λ a are the Gell-Mann matrices. Finally L gf and L ghost are the gauge-fixing and ghost part of the Lagrangian, respectively.
The partition function in QCD can be written as where S QCD is the Euclidean action for QCD. In the second line, we have integrated over the fermions which can be done exactly since L QCD is bilinear in the quark fields. Moreover, S g is the Euclidean action for the gluons and This yields The last equation shows that the fermion determinant is manifestly positive. As in the case of finite isospin chemical potential, QCD in a magnetic field is also free of the sign problem, and one can therefore carry out lattice simulations. Interestingly, the combination of finite isospin and magnetic field is free of the sign problem only if the charges of the u and d-quark are the same. This is of course not real QCD, but it offers the possibility to compare lattice predictions with those of low-energy effective theories and models.
In this review, we will discuss (inverse) magnetic catalysis and the phase diagram of QCD in a strong magnetic background, paying attention to recent developments. There are other reviews [9][10][11][12] focusing on different aspects of the field. The paper is organized as follows. In the next section, we discuss the physics of magnetic catalysis at T = 0. In Sec. 3, we introduce the Polyakov loop and discuss magnetic catalysis in model calculations at nonzero temperature. In Sec. 4, we review inverse magnetic catalysis on the lattice focusing on the competing sea and valence effects. In Sec. 5, the improvement of models to incorporate inverse magnetic catalysis is discussed and in Sec. 6, we summarize. In Appendix A, we discuss renormalization of the quarkmeson model in the on-shell scheme, while in Appendix B, we show how the parameters of the model are fixed.
2 Magnetic catalysis at zero temperature Magnetic catalysis can be defined as 1. The magnitude of a condensate is enhanced by the presence of an external magnetic field B if the condensate already is present at vanishing field.
2. An external magnetic field induces symmetry breaking and a nonzero value of a condensate when the symmetry is intact for B = 0.
The condensate is the expectation value of a field, which can be either fundamental or composite. The expectation of a scalar field φ in low-energy models is an example of the former, whileψψ is the chiral condensate in e.g. the NJL model or QCD is an example of the latter. One refers to the second case as dynamical symmetry breaking by a magnetic field. We will discuss both cases below. The first papers on magnetic catalysis at T = 0 appeared three decades ago in the study of the NJL model in three dimensions [13]. Shortly thereafter in the linear sigma model [14] and the NJL model in two dimensions [15][16][17]. Since then it has been demonstrated in QED [18], chiral perturbation theory [19][20][21], in the Walecka model in nuclear physics [22], and also on the lattice, see e.g. [23][24][25][26][27].
After symmetry breaking, the sigma field has a nonzero expectation value φ 0 . The classical potential is The tree-level relations between the parameters of the Lagrangian m 2 , λ , g, and h and the physical masses m σ and m π , the pion decay constant f π , and the quark mass m q are Using the relations (10)-(11), we obtain where we have introduced ∆ = gφ 0 . The minimum of the classical potential is given by ∆ = g f π .
The classical potential has by construction its minimum at ∆ = m q or φ 0 = f π . In the large-N c limit, the mesons are included at tree level, while we include the Gaussian fluctuations of the fermions. Including the one-loop corrections from the fermions using a minimal subtraction scheme, leaves a renormalized one-loop effective potential that depends on the renormalization scale Λ . The minimum of the effective potential therefore depends on Λ . In order to ensure that the one-loop effective potential has its minimum at φ 0 = f π for zero magnetic field B, several methods have been used in the literature. One method is simply to subtract the one-loop contribution to the effective potential for B = 0. Then the renormalization scale dependence drops out and the correction to Eq. (12) is a finite B-dependent term that vanishes for B = 0. However this is inconsistent since one includes fermion fluctuations in the effective potential at finite magnetic field, but not for B = 0. Moreover, it is also incorrect since Eqs. (10)-(11) are tree-level relations that receive radiative corrections. One can also choose a specific value for Λ such that the one-loop correction to the position of the minimum of the effective potential vanishes. In this case, one has included quantum fluctuations also for B = 0, but again, the tree-level relations between the parameters of the Lagrangian and physical quantities receive loop corrections. In order to be consistent, the parameters of the Lagrangian must be determined to the same order in the loop expansion as one calculates the effective potential. The solution to the problem is to combine the minimal subtraction scheme with the on-shell scheme [28][29][30][31]. In this way one includes loop corrections to Eqs. (10)- (11), while at the same time ensures that the effective potential has its minimum at ∆ = g f π . Details of the renormalization of the oneloop effective potential in the large-N c limit can be found in Appendix A and the parameter fixing in Appendix B. It reads where x f = ∆ 2 2|q f B| and ζ (a, x) is the Hurwitz zeta-function. Here and in the remainder of the paper, ζ (1,0) (n, The first four lines of the one-loop effective potential are independent of the magnetic field and this part was first calculated in Ref. [32]. The last line is the B-dependent correction to V 0+1 . Note also that final result is independent of the renormalization scale Λ . In Fig. 1, we show the effective potential divided by f 4 π at T = 0. The black line is the tree-level potential Eq. (12), while the green and blue lines are the one-loop effective potential Eq. (13) for |eB| = 0 and |eB| = 10m 2 π , respectively. We have used m σ = 600 MeV, m π = 140 MeV, f π = 93 MeV, and m q = 300 MeV. The classical potential as well the oneloop effective potential with |eB| = 0 both have a minimum at ∆ = g f π by construction. Notice, however, that the latter is significantly deeper. The blue line with |eB| = 10m 2 π , shows that the minimum of the effective moves to a larger value, i. e. the system exhibits magnetic catalysis. While the above clearly demonstrates magnetic catalysis numerically, we would like to gain insight in the mechanism behind the effect. Instead of analyzing Eq. (13), we will discuss the gap equation in the NJL model.
In order to simplify the discussion, we will consider the NJL model with a single quark flavor and color, N f = N c = 1 with electric charge q f . In the chiral limit, the Lagrangian is [33] L =ψiγ µ ∂ µ ψ + G 2 (ψψ) 2 + (ψiγ 5 ψ) 2 .
This Lagrangian has a U(1) V × U(1) A symmetry. We introduce the gap M = −G ψψ and linearize the interaction terms, writing (ψψ) 2 ≈ ψψ 2 +2 ψψ ψψ and (ψiγ 5 ψ) 2 ≈ 0. M is now an effective quark mass arising after breaking the axial symmetry spontaneously, i.e. when ψψ = 0. In the mean-field approximation, we perform the Gaussian integral over the fermion field giving rise to the following oneloop effective potential, The gap equation for M is found by extremizing V 0+1 which yields Conventionally, since the NJL model is non-renormalizable, one has used a three-dimensional or a four-dimensional momentum cutoff Λ to regulate divergences. If Λ is a fourdimensional cutoff, the gap equation (16) M = 0 is always a solution, however for G > G c = 4π 2 Λ 2 there is also a nontrivial solution. Thus for G larger than the critical value 4π 2 Λ 2 , quantum fluctuations induce symmetry breaking in the model.
At finite magnetic field, the partial derivative in Eq. (14) is replaced by the covariant derivative and we add a term 1 2 B 2 to the effective potential. The gap equation becomes where M 2 B = M 2 + |q f B|(2k + 1 − s), p 2 = p 2 0 + p 2 z and q f is the charge. The divergences in Eq. (18) can be isolated by adding and subtracting the right-hand side of Eq. (16). The right-hand side of Eq. (18) minus the subtracted term is finite and is conveniently evaluated using dimensional regularization in the same way as done appendix A [34,35]. Finally, we impose a four-dimensional cutoff on the added term as in Eq. (16). Factoring out the trivial solution M = 0, this yields the regularized gap equation where x f = M 2 2|q f B| . This equation has only a nonzero M as solution. For G < G c , the solution is [33,36] In the limit |q f B| → 0, this solution connects to the trivial solution M = 0. In the lowest Landau level approximation, the gap equation has solution M 2 = Λ 2 e −4π 2 /G|q f B| , which is reminiscent of Eq. (20) if we identify the cutoff Λ with |q f B|. We can then think of magnetic catalysis as a 1+1 dimensional phenomenon, i.e. a dimensional reduction from 3+1 dimensions has taken place. The functional form of the gap equation is the same as for the gap equation in BCS theory of superconductivity as well as the gap equation found in the large-N limit of O(N)-symmetric nonlinear sigma model in 1+1 dimensions. The 1+1 dimensional nature of magnetic catalysis raises the question of whether this phenomenon is in conflict with the Coleman theorem, which forbids spontaneous symmetry breaking in less than two spatial dimensions at zero temperature [37]. As pointed out in Ref. [36], the fieldψψ is neutral with respect to the magnetic field. The neutral pion is the associated Goldstone boson that appears after breaking the U(1) symmetry. The charged pions are now massive even in the chiral limit.
There are other ways of regularizing the gap equation (18) or the fermion contribution to the one-loop effective potential (A.1), for example Schwinger's proper time method [38]. Let us illustrate this by computing the corresponding bosonic functional determinant, which shows up in chiral perturbation theory. It is based on the representation in Euclidean space where the sum over Landau levels k as well the momentum integral over p is convergent. The result is For ε = 0, the integral is divergent for small s, i.e. for large momentum. By adding and subtracting the divergent terms, we can isolate the divergences. One finds The integrals in the first line are divergent for ε = 0, while the last integral is convergent. The divergences show up as poles in ε. The first term in Eq. (23) is a vacuum energy counterterm while the second term corresponds to charge and wavefunction renormalization [38]. The last integral in Eq. (23) can be calculating exactly and involves the Hurwitz zeta function. Using the proper time method with the momentum integrals evaluated in d = 2 − 2ε dimensions yields the same results as those obtained by combining dimensional regularization and zeta-function regularization, as done in Appendix A. Alternatively, one can evaluate Eq. (23) with ε = 0 using a cutoff 1/Λ 2 as the lower limit of the sintegration in the divergent integrals. The regularization methods discussed so far separates in clean way the B-independent divergences from the Bdependent terms whether they are finite or divergent. 1 There are other regularization methods that do not separate this contributions, for example a sharp cutoff imposed directly on the integral in Eq. (22) or a form factor that is a function of e.g. p 2 z + 2k|qB|(2k + 1 − s). One has to be careful choosing such regulators since nonphysical oscillations may result [39,40].
The on-shell scheme used to obtain the final result Eq. (13) has two important virtues as first pointed out in Ref. [41]. By considering the small-B (large-x f ) behavior of the Hurwitz zeta-function, one finds that the only contributions at order B 2 comes from the renormalized magnetic field term 1 2 B 2 , It also ensures that the magnetic-field contribution to the effective potential and the magnetization vanishes in the limit m q → ∞. Both properties are expected from a physical point of view. It leads to a paramagnetic vacuum, in agreement with the hadron resonance gas model calculations [41] and lattice QCD lattice simulations [42]. Other renormalization schemes, such as the (modified) minimal subtraction scheme are connected to the above by a finite renormalization. However, the effective potential and the magnetization then grow logarithmically in the limit m q → ∞. Similar remarks apply to the (P)NJL model. Regularizing the model by imposing a UV cutoff Λ on the divergent integrals in the fermionic version of Eq. (23), the authors of [43] show that it predicts diamagnetic behavior for low values of B and paramagnetic behavior for large magnetic fields. By defining a subtraction procedure that resembles the renormalization of the magnetic field in the on-shell scheme, their effective potential leads to paramagnetic behavior as seen on the lattice. After having discussed magnetic catalysis in models, we now turn to lattice gauge theory. The first lattice simulations were carried out for an SU(2) gauge group for magnetic field strengths up to √ eB ∼ 3 GeV in the quenched approximation, i.e. setting the quark determinant to unity [23]. The simulations confirmed that the quark condensates is enhanced by the magnetic field and that the enhancement is qualitative linear with eB. The quark condensate itself was calculated using the Banks-Casher relation [44], which relates the density of eigenvalues close to zero of the Dirac operator and the condensate. Their calculations showed a monotonic increase of the spectral density for typical gauge field configurations. This enhancement induced by the magnetic field can be considered the basic mechanism behind magnetic catalysis. Below, we will discuss this mechanism further, here it suffices to add that the enhancement of the spectral density as a function of B for typical gauge configuration is also seen in full QCD [45]. Even in the free case, there is a proliferation of small eigenvalues due to the degeneracy of states, which in a constant magnetic field is proportional to |qB| [19].

Magnetic catalysis at nonzero temperature
In the previous section, we reviewed magnetic catalysis at T = 0 in some detail. A survey of the literature shows that it is a robust feature of lattice simulations as well model calculations: Magnetic catalysis does not depend on particular values of the masses or couplings. Since the condensate increases as a function of the magnetic field, it should raise the transition temperature for the chiral transition. This expectation was made explicit a long time ago in Ref. [19]. A number of model calculations have confirmed this expectation, e.g. [13,[46][47][48][49][50][51][52], although some of them suggest that the chiral and deconfinement transition split for larger values of eB.
The effective potential of the quark-meson model in the large-N c approximation at finite temperature is where E f = p 2 z + ∆ 2 + |q f B|(2k + 1 − s) and V 0+1 is given by Eq. (13). In Fig. 2, we show the results of a typical calculation where the quark-meson model was used. The curves show the transition temperature for the chiral transition as a function of |qB|/m 2 π in the chiral limit (green points) and at the physical point (red points). At B = 0, the gap between the two critical temperatures is approximately 10 MeV, which decreases as |qB| grows. In both cases, it is clear that the transition temperature increases with the magnetic field. Here the transition temperature was defined as the inflection point of the curve φ 0 (T ) at the physical point and φ 0 (T ) = 0 in the chiral limit. An alternative definition of the critical temperature is the peak of the chiral susceptibility 2 2 Using the peak of ∂ Φ ∂ T , where Φ is the Polyakov loop, yields a transition temperature for deconfinement, which is very close to the chiral transition temperature. In QCD, two transitions take place as one increases the temperature, namely the chiral transition and the deconfinement transition. Lattice calculations suggest that chiral symmetry is "restored" at a temperature of approximately T χ c = 155 MeV [53][54][55][56][57] though strictly speaking the transition is only a crossover. The crossover temperature is defined by the peak of the chiral susceptibility. This temperature is slightly less than the crossover temperature for the deconfinement transition, T dec c =170 MeV. However this temperature difference is observable dependent. In most cases, T dec c has been determined by the behavior of the Polyakov loop. Recently, it has been defined by the behavior of the quark entropy and in this case the two crossover temperatures agree within errors [57].
We will next discuss the Polyakov loop and how it can be incorporated in model calculations. The Wilson line is defined as where P denotes path ordering, A 4 = iA 0 and A 0 = t a A a 0 . The Polyakov loop operator l is the trace of the Wilson line (27). Together with its Hermitian conjugate, it is defined as where N c is the number of colors. The expectation values of l and l † are denoted by Φ andΦ. Under the center group Z N c of the gauge group SU(N c ), the Polyakov loop transforms as Nc Φ with n = 0, 1, 2...N c − 1. In pure-glue QCD it is an order parameter for confinement, while for QCD with dynamical fermions it is only an approximate order parameter [58]. Note also that Φ =Φ at zero density, i.e. for µ f = 0.
For N c = 3 and in the Polyakov gauge, one can write a nonabelian background gauge field as Introducing the fields φ 1 = 1 2 β A 3 4 and φ 2 = 1 2 √ 3 β A 8 4 , the thermal Wilson line reads for constant gauge fields Since the Polyakov loop is an approximate order parameter for deconfinement, the strategy put forward in Ref. [59] is to write down a phenomenological effective potential for Φ,Φ and the chiral condensate that describes the thermodynamics of the system. This potential consists of a gluonic part U(Φ,Φ) as well as a matter part. The term U(Φ,Φ) is constructed such that it reproduces the pure-glue pressure calculated on the lattice [60]. A number of different forms of U(Φ,Φ) have been proposed [61][62][63][64]. In Ref. [61], they used a Polynomial expansion incorporating the Z 3 center symmetry, Here the coefficients are and T 0 = 270 MeV, the transition temperature for pure-glue QCD [60]. A drawback of the proposed pure-glue potentials is that they are independent of the number of flavors n f . The transition temperature for B = 0 depends on the number of flavors and one should incorporate the back-reaction from the fermions to the gluonic sector [62]. This is done by using an n f -dependent T 0 . Once the coupling between the gluonic sector and the matter sector has been implemented, the two transitions take place at approximately the same temperature: The chiral transition moves to larger temperatures, while the deconfinement transition moves to lower temperatures. Finally, the Polyakov-loop potential is coupled to the matter sector via replacing the partial derivatives in the fermionic part of the Lagrangian by covariant ones including the constant background gauge field. This is implemented by making the substitution in Eq. (25). In the same way, the Fermi-Dirac distribution function is generalized, For small values of the Polyakov loop, Φ ≈ 0,Φ ≈ 0 Eq. (36) reduces to a Fermi-Dirac distribution with excitation energy 3E f , i.e. that of three quarks. For large temperatures, when Φ ≈ 1,Φ ≈ 1, the excitation energy is E f , which is the distribution function of deconfined quarks.
To the best of our knowledge, there are no systematic studies of the transition temperature as a function of the magnetic field B in various approximations. However, some interesting results using the quark-meson model exist.   3 shows the normalized transition temperature from Ref. [50] in two approximations using the functional renormalization group (FRG) [65]. In this approach, one solves a flow equation for the effective potential numerically by lowering a sliding scale k from an initial UV cutoff k = Λ (where the effective potential is equal to the classical potential) down to k = 0. The bare parameters at k = Λ are tuned such that one obtains the physical values of the masses and the pion decay constant in the vacuum (i. e. for k = 0). In this way, all quantum and thermal fluctuations are included. The black solid line is the mean-field result, i.e. the bosons are excluded from the flow equation, whereas the brown line is the result using the functional renormalization group. Clearly, the addition of bosonic fluctuations increases the transition temperature significantly.
In Fig. 4, we show the transition temperature at the physical point using the functional renormalization group [52]. The green points are the results without the Polyakov loop, whereas the blue points are the results including it. Clearly, the Polyakov loop lowers the transition temperature for fixed B, but it is still increasing as we increase the magnetic field. The above FRG results are obtained in the so-called localpotential approximation. In Ref. [66], the authors added the effects of wavefunction renormalization and the curve for the critical temperature lies between mean-field and the local- potential approximation. Thus the coupling of the Polyakov loop to the chiral sector is not sufficient to reproduce (qualitatively) the results seen on the lattice.

(Inverse) Magnetic catalysis on the lattice
After having discussed magnetic catalysis in low-energy models and theories of QCD, we next consider QCD lattice simulations. In the past decade, there have been a number of lattice calculations of QCD in a magnetic field [23-27, 45, 67-73], which have improved our understanding of QCD in a magnetic background.
In order to discuss (inverse) magnetic catalysis as seen on the lattice, it is advantageous to take a look at the pathintegral representation of a number of expectation values. The QCD Lagrangian is bilinear in the quark fields ψ f and so one can integrate over them, giving for the partition function as a path integral over gauge configurations A µ We can think of P = 1 Z (B) e −S g det(D /(B) + m) as a measure that depends on the gauge-field configuration A µ , the magnetic field, and the quark masses. Note that the B-dependence is in the functional determinant as well as the the trace of the propagator. In order to study the contributions to the quark condensate coming separately from the change of the operator and the change of the measure, it is convenient to introduce the valence and sea contributions defined as This can be thought of as an expansion of the quark condensate around B = 0. A priori, the sum of the two contributions needs not add up to the total quark condensate unless we are at small fields. However, it turns out that writing the condensate as a sum of the valence and sea contribution is remarkably good. This is clearly demonstrated in Fig. 5 from Ref. [26], which shows the relative increment r of the valence and sea contributions, their sum as well as the complete results for the quark condensate as a function of a dimensionless quantity b. The relative increment is defined as where ψψ is the average of the u and d quark condensates. Within error, the additivity is confirmed for values of b up to 8, which corresponds to magnetic fields up to eB = (500MeV) 2 [26]. It is also of interest to notice that both contributions work in the same direction, namely to increase the quark condensate as B grows. This is unlike what happens at temperatures around the critical temperature T c , as we shall see below. As pointed out in Ref. [45], ψψ sea can be thought of as the quark condensate of an electrically neutral fermion flavor coupled to an electrically charged fermion flavor, since the magnetic field only appears in the functional determinant and not in the propagator. On the other hand, ψψ val is reminiscent of the expression of the quark condensate in model calculations, except in models one does not integrate over gauge-field configurations.
Let us now turn to finite temperature. Inverse magnetic catalysis seems to have two somewhat different meanings in the literature. The first meaning corresponds directly to the concept magnetic catalysis discussed above: it simply means that a condensate, for example ψψ , decreases with the magnetic field at a fixed temperature. The second meaning is that the transition temperature itself is a decreasing function of the magnetic field.
The first finite-temperature lattice simulations were carried out in [23,24] for SU (2) gauge theory in the quenched approximation, focusing on the B-dependence of the chiral condensate for temperatures below the transition. In twoflavor QCD, simulations at finite temperature were carried out for pion masses in the range 200 − 480 MeV in Ref. [26] and it was concluded that the chiral and deconfinement transitions take place at the same temperature and that they increase slightly with the external magnetic field. The increase of the transition temperature with B is, at least qualitatively, in agreement with model calculations. Bali et al [27,68] carried out lattice simulations at the physical point with 2+1 flavors, i.e. for quark masses that correspond to m π = 140 MeV, and the result was somewhat surprising: The transition temperature turned out to be decreasing as B increases. The different behavior of T c is not a consequence of the different pion masses, rather it results from lattice artefacts and that the results of [26] were not continuum extrapolated. Today there is consensus that the chiral transition temperature is a decreasing function of the magnetic field. This behavior is illustrated in Fig. 6, which shows the results of a recent lattice simulation [70], namely the transition temperature in MeV as a function of the magnetic field eB in GeV for three different pion masses. The pion mass is 343 MeV (red points), 440 MeV (blue points), and 664 MeV (green points), which is much larger than the physical pion mass of 140 MeV. The transition temperature increases as a function of the pion mass for fixed value of B, which is also known from B = 0 calculations. In Ref. [45], the authors carried out a thorough analysis of the quark condensate around the critical temperature to understand the behavior of the transition temperature, focusing on disentangling the valence and sea effects. The valence contribution Eq. (39) can also be written as where the subscript indicates that the quark determinant is without a magnetic field. The spectral density of the quark operator for different values of the magnetic field is shown in Fig. 7. From the figure, it is evident that there is an increase in the spectral density around zero with increasing magnetic field. The corresponding ensemble was generated at finite temperature, T = 142 MeV and for vanishing magnetic background [45]. The Banks-Casher relation [44] then implies an increase of the valence contribution.
Note that Eq. (44) reduces to the valence contribution Eq. (42) if one replaces ∆ S f (B) by unity. Fig. 8 from Ref. [45] shows a scatter plot of the condensate as a function of the change in the action ∆ S f (B) due to the magnetic field. In this plot, the magnetic field strength is eB ≈ 5 GeV 2 and T close to the transition temperature. Each point represents a gauge configuration and they were generated at vanishing magnetic field. The plot suggests that larger values of the condensates correspond to larger values of the weight e −∆ S f (B) and therefore suppresses the weight of the associated gauge configuration. As a result, this counteracts the valence effect, and leads to a decrease in the critical temperature. For pion masses that are not too large, it also leads to a decrease of the condensate itself (see discussion below). In Refs. [70,71], the effects of varying the pion mass on the quark condensate as a function of the temperature have been studied in detail (again 2+1 flavors). Fig. 9 from [70] shows the difference between the quark condensates as a function of the temperature at B = 0 and eB = 0.425 GeV 2 (blue data points) and eB = 0.85 GeV 2 (red points) for three values of the pion mass. The authors find that the sea contribution is a decreasing function of B around T c for the different values of the pion masses, while the valence contribution is on the other hand an increasing function of the magnetic field for all temperatures and pion masses. In the upper panel it is clear that the sea contribution wins the competition around the transition temperature implying inverse magnetic catalysis in the strict sense of the word. This effect can barely been seen in the middle panel and is completely absent in the lower panel. In other words, a decreasing function of the transition temperature does not imply that the chiral condensate decreases as a function of temperature and it is therefore not clear that the latter is the driving mechanism of the former [70]. This nontrivial behavior was also demon-strated in [71], where the authors fixed the magnetic field to eB = 0.6 GeV 2 and varied the pion mass. In QCD, there is inverse magnetic catalysis for pion masses up to 500 MeV, and magnetic catalysis for larger values 3 We finally comment on the nature of the chiral transition and the temperature as a function of B. The simulations have been done with magnetic fields up to eB = 1 GeV 2 . They all show an analytic crossover and that the transition temperature is a decreasing function of B. However, it has been conjectured that the transition would start increasing again for sufficiently large temperatures, a phenomenon dubbed delayed magnetic catalysis. In Ref. [69] the author went as high as eB = 3.25 GeV in the simulations. The transition remains a crossover (albeit sharper), there is no sign of delayed magnetic catalysis, and the chiral and deconfinement transitions coincide. The sharper crossover suggests that there may be a critical point for even larger values of the magnetic field. For asymptotically large fields, QCD can be mapped onto an anistropic pure-glue theory [74]. This theory was simulated on the lattice and strong evidence for a first-order transition was found [69]. This implies the existence of a critical point and its position was estimated to be at eB 10 GeV 2 .

Improvements of models
The failure of models to correctly describe the behavior of QCD around the critical temperature, even after the introduction of the Polyakov loop, has lead to significant efforts to improve them, see e.g. [75][76][77][78][79][80][81][82][83][84][85]. The temperature T 0 that enters the Polyakov-loop potential depends on the number of flavors, it is therefore a reasonable assumption that it also depends on the magnetic field. In Ref. [75], the authors fitted the strange-quark susceptibilities from their calculations in the entangled PNJL model to the lattice results of Ref. [68]. Their ansatz for the field dependence was a simple polynomial in (eB) 2 up to quadratic order, giving two fitting parameters, An interesting feature here is that the model predicts a firstorder transition for magnetic fields larger than approximately eB = 0.25 GeV 2 . As mentioned above, such a critical point is expected in QCD, albeit at much larger magnetic fields [69]. The crossover nature of the chiral transition was a guide for the authors of Ref. [76], trying to incorporate the decreasing behavior of the transition temperature in the (Polyakovloop extended) two-flavor quark-meson model. In their meanfield analysis, they allowed the Yukawa coupling to vary with the magnetic field, g = g(B) using the boundary value g(0) = 3.3. This value is indicated by the vertical dotted line in Fig. 10 and corresponds to a fixed quark mass in the vacuum. The critical temperature as a function of the Yukawa coupling for three values of the magnetic field is also shown in Fig. 10. The grey shaded region indicates the values of g for which the transition is first order.
However, to obtain a transition temperature which is decreasing with the magnetic field, any curve g(B) must start at g(0) = 3.3 and successively cross the dashed (red) and solid (black) curves. One therefore soon enters the shaded region which indicates a first-order transition in the QM model. Thus a function g(B) can simply not describe the correct Bdependence of the transition temperature, while at the same time having a crossover transition.
Similar approaches have been used in the NJL model allowing the coupling G to depend on both T and B. For example, in Ref. [77], the authors fix a set of parameters to get a reasonable fit to the lattice data for the sum of the light quark condensates. The form of the B-dependent coupling was motivated by the running of the coupling in QCD for strong magnetic fields [86]. One particular appealing idea was recently put forward in Ref. [87] (see also Ref. [88]). Only T = 0 physics from the lattice is used as input to improve the model. In other words, there is no fitting to lattice data for T c or a coupling that depends both upon B and T . The authors first performed a determination of the baryon spectrum at the physical point as a function of magnetic field using lattice simulations. The authors focused on strong magnetic fields, which are relevant for the phase diagram. Making the simple assumption that the baryon masses can be written as the sum of the masses of their constituents, they derived Bdependent constituent quark masses. This was used as input in the PNJL model at zero temperature: Using the gap equation with the B-dependent quark masses, a B-dependent four-fermion coupling was obtained. The B-dependent constituent quark masses as well as G(B) are decreasing functions of the magnetic field. For B = 0, M 2 = 0.097 MeV 2 and G(0) = 12.8 GeV −2 and for the largest magnetic field used (|eB| = 0.6 GeV 2 ), M 2 = 0.079 MeV 2 and G(B) = 6.7 GeV −2 . Fig. 11 shows the normalized average quark condensate at T = 0 as a function of the magnetic field. The dashed black line is obtained from a next-to-leading order (oneloop) calculation in chiral perturbation theory [19,20], the light-blue points are obtained from the standard PNJL model, while the green points are obtained from the lattice-improved PNJL model. Finally, the red band shows the results from lattice simulations including errors. The plot has several interesting features. Firstly, all low-energy approaches are in good agreement with lattice results for low values of the magnetic field. For larger values, both χPT and PNJL underestimate the quark condensate. This is in contrast with the lattice-improved PNJL model, which is in quantitative We next consider the finite-temperature calculations of Ref. [87]. Fig. 12 shows the quark condensate as a function of T for different strengths of the magnetic field. The dashed lines are the predictions from the standard PNJL model without error bands, while the solid bands are the lattice-improved PNJL model predictions. We first note that quark condensate at T = 0 for fixed magnetic field is higher for the lattice-improved PNJL model in agreement with Fig. 11. This behavior persists for low temperature, for larger temperatures, however, the condensate drops faster compared to that calculated using the standard PNJL model. Thus for fixed magnetic field, the transition temperature defined as the inflection point of the quark condensate moves to a lower value of T compared to the PNJL model. This in itself is not enough to conclude that we have inverse magnetic catalysis, but the effect is so strong that the inflection point moves to the left as a function of B so the transition temperature decreases. Fig. 13 shows the normalized transition temperature for the chiral transition, as defined by the inflection point of the quark condensate, as a function of eB in units of GeV 2 . The pink band is from the lattice results of Refs. [68]. The width of the band indicates the errors of the simulations. The lightblue points are the results of a calculation from the standard PNJL model showing that the transition temperature increases as the magnetic field grows. This is in sharp contrast to the lattice-improved PNJL model, where the results are shown by the green points including uncertainties coming from the lattice determination of the baryon masses. Given the large uncertainties, the results for the transition temperature are in good agreement with the simulations. Similarly, the analytic crossover found also agrees with lattice results.

Summary and final remarks
The idea of magnetic catalysis at zero temperature has been around for three decades after its discovery in the NJL model. It is a robust phenomenon. For large values of the magnetic field, it can be understood in terms of dimensional reduction from 3 + 1 to 1 + 1 dimensions (or from 2+1 to 0+1 dimensions). Lattice simulations in the last decade have improved our understanding of the effect significantly by showing that both the valence and sea effect contribute in a nontrivial way.
Since the discovery of the the decrease of the chiral transition temperature with increasing magnetic field, a lot of work has been devoted to incorporate this feature in models. This includes T -and B-dependent couplings and Polyakovloop potentials. Fitting parameters can be considered an indirect way of incorporating the sea effect and requires input from lattice simulations at finite temperature. In our opinion, a cleaner approach is provided by the work [87], which uses lattice input at T = 0 only. In much the same way as one uses experimentally measured meson masses in the vacuum, they use B-dependent baryon masses measured on the lattice as input in their PNJL-model calculations, although assuming that the baryon mass is the sum of its constituents perhaps is somewhat simplistic.
The mass of a particle is given by the pole of the propagator.
In the on-shell scheme, the sum of the self-energy evaluated on-shell and the counterterms vanishes, In addition, the residue of the propagator evaluated on shell is unity. This implies The self-energies are where the integrals A(m 2 ) and B(p 2 ) are defined in Eqs. (B.27)-(B.28). We have omitted the tadpole diagram which is oneparticle reducible. It is cancelled by a counterterm, when we impose the condition that φ 0 = f π . The counterterms are given by the expressions . . where we have added the counterterm δt for the one-point function. We next need to relate the above counterterms to the counterterms of the parameters of the Lagrangian. These relations follow immediately from Eqs. (10)- (11), (B.13) In the large-N c limit, δ m 2 q = 0 implying that δ g 2 OS = −g 2 δ f 2 π f 2 π . There is also no loop correction to the quark-pion vertex and δ Z ψ = 1. This implies that the associated counterterms must cancel as well, leading to δ g 2 OS = −g 2 δ Z OS π . We can therefore write (B.14) The counterterm δ h OS is found from the one-point function.
At tree level, we have t = h − m 2 π f π = 0, which yields δt = δ h OS −δ m 2 π f π −m 2 π δ f π = or δ h OS = δt +δ m 2 π f π + 1 2 m 2 π f π δ Z OS π . Finally, we need the counterterm for the electromagnetic field, with the integral B(p 2 ) defined in Eq. (B.28). Renormalization is then carried out by making the substitution B 2 → B 2 (1 + δ Z OS A ). We note that the on-shell scheme is not welldefined when the fermions are massless. In that case, the (modified) minimal subtraction scheme may be used. Since the bare parameters are independent of the renormalization scheme, we can immediately write down the relations between the renormalized parameters in the on-shell and MS schemes. For example g 2 MS +δ g 2 MS = g 2 OS +δ g 2 OS . From Eq. (B.9), we find ε . From this, one finds the running coupling g 2 MS using g 2 MS = g 2 OS + δ g 2 OS − δ g 2 MS and given by Eq. (B.19). The running parameters are (B.20) The running parameters satisfy renormalization group equations that follow from Eqs.