The QCD phase diagram for external magnetic fields
The effect of an external (electro)magnetic field on the finite temperature transition of QCD is studied. We generate configurations at various values of the quantized magnetic flux with flavors of stout smeared staggered quarks, with physical masses. Thermodynamic observables including the chiral condensate and susceptibility, and the strange quark number susceptibility are measured as functions of the field strength. We perform the renormalization of the studied observables and extrapolate the results to the continuum limit using and lattices. We also check for finite volume effects using various lattice volumes. We find from all of our observables that the transition temperature significantly decreases with increasing magnetic field. This is in conflict with various model calculations that predict an increasing . From a finite volume scaling analysis we find that the analytic crossover that is present at persists up to our largest magnetic fields , and that the transition strength increases mildly up to this .
Keywords:Lattice QCD, finite temperature, external magnetic field
The properties of QCD in strong magnetic111Throughout the paper ‘magnetic’ refers to electromagnetic i.e. not chromomagnetic. fields are relevant for at least three important physical situations. First, cosmological models suggest that extremely strong magnetic fields ( GeV) could be produced during the electroweak phase transition of the early universe. This effect might also have an impact on subsequent strong interaction processes Vachaspati:1991nm (). Second, large magnetic fields ( MeV) are present in the interior of dense neutron stars called magnetars Duncan:1992hi (). Finally, in a noncentral heavy ion collision the spectators – being two beams of positive charges moving in opposite directions – also create an intense magnetic field which, depending on the centrality and the beam momentum, reaches up to GeV for RHIC and GeV for the LHC Skokov:2009qp (). This magnetic field is external since it is produced by the spectators, and though it has a very short lifetime (of the order of fm/), the magnetic ‘impulse’ coincides with the generation of the quark-gluon plasma and thus may have a significant effect on the properties of the transition.
For noncentral heavy ion collisions, an exciting consequence of the interplay between the strong magnetic field and the nontrivial topological structure of the quark-gluon plasma is the so-called chiral magnetic effect Kharzeev:2007jp (); Fukushima:2008xe (). This effect creates an electric current of quarks (anti)parallel to the external magnetic field, which may result in a preferential emission of charged particles perpendicular to the reaction plane, leading to event-by-event CP-violation Kharzeev:2004ey (). Recent measurements from the STAR experiment at RHIC Voloshin:2008jx (); :2009uh () and the ALICE experiment at the LHC Selyuzhenkov:2011xq () are in qualitative agreement with this picture, however, the interpretation of these results is still under discussion Wang:2009kd (); Muller:2010jd (); Voronyuk:2011jd ().
Because of this high phenomenological relevance, the effect of a finite magnetic field on the strong interactions has been studied extensively in the last years, both using model calculations and lattice simulations. In particular, the structure of the QCD phase diagram in the plane has received increasing attention recently. Calculations have been carried out within various low energy effective models of QCD. In the linear sigma model coupled to quarks and the Polyakov loop it was observed that the transition temperature increases with Mizher:2010zb (). Furthermore, a splitting between the deconfinement and chiral transitions was predicted to take place for large external fields. The strength of the transition was also observed to increase, which eventually results in a first-order phase transition Fraga:2008um (). Similar conclusions with respect to the increase in and in the strength of the transitions were drawn from studies of the Nambu-Jona-Lasinio model and extended versions thereof, like the EPNJL and PNJL models Gatto:2010pt (); Gatto:2010qs (), see also Osipov:2007je (), and the nonlocal PNJL model Kashiwa:2011js ().
The presence of the external magnetic field was also shown to increase the transition temperature within other types of models like in the Sakai-Sugimoto model of large gauge theories Johnson:2008vna (), in the Gross-Neveu model in lower dimensions Kanemura:1997vi (); Klimenko:1992ch (), in 2+1 dimensional QED when described by Schwinger-Dyson equations Alexandre:2000yf () and within the holographic approach Evans:2010xs (). However, the opposite effect of a decreasing deconfinement transition temperature was predicted using chiral perturbation theory for two quark flavors Agasian:2008tb (). A decrease in was also observed in the linear sigma model if the quark vacuum contributions are neglected Mizher:2010zb () and in the Sakai-Sugimoto model with nonzero chemical potential Preis:2010cq (). We mention that lattice simulations indicate a reduction of the transition temperature of QCD in an external chromomagnetic field Cea:2005td (); Cea:2002wx (); Cea:2007yv ().
The phase diagram is in most cases predicted by studying chiral symmetry breaking, i.e. the behavior of the chiral condensate or the dynamical quark mass as a function of . Most of the low energy models agree that chiral symmetry breaking is enhanced as the magnetic field grows Gusynin:1995nb (); Nam:2011vn (); Boomsma:2009yk (); in particular the value of the chiral condensate was found to increase linearly with in leading order Shushpanov:1997sf (); Cohen:2007bt (); Agasian:2001hv (). The condensate also increases with – although with a quadratic leading order – within the AdS/CFT duality picture Zayakin:2008cy (), and with in holography Evans:2010xs (). On the other hand, it was also conjectured that the running of the strong coupling in the presence of magnetic fields may modify this magnetic catalysis, and even turn the effect around to make the dynamical mass decrease with in some regions Miransky:2002rp ().
In recent lattice simulations with flavors of staggered quarks D'Elia:2010nq () the chiral condensate was observed to grow with the external field for any temperature in the transition region. The size of this effect was however found to be different for different values of , resulting in an increase in both the pseudocritical temperature and the strength of the transition. Furthermore, according to the findings of D'Elia:2010nq (), the relative change in is of the order of one percent for several larger-than-physical pion masses.
In this paper our aim is to perform a similar lattice study, but with improved gauge and smeared fermionic actions and with flavors of quarks, at the physical pion mass, and extrapolate the results to the continuum limit. We include the magnetic field in the fermion determinant and study its effect dynamically to investigate how the strength and the pseudocritical temperature of the QCD transition change as the external magnetic field is switched on. We explore a wide temperature region around the zero-field pseudocritical temperature , for various values of the magnetic field, ranging from MeV to GeV, i.e. covering the regions that are phenomenologically interesting for noncentral heavy ion collisions and for the evolution of the early universe.
This paper is structured as follows: first the implementation of the magnetic field on the lattice is described. Then we define the observables of interest, including the chiral condensate, chiral susceptibility and strange quark number susceptibility, and discuss their renormalization at zero and nonzero . After presenting the simulation setup and the details of the analysis we show our results for the transition temperature and the width of the transition.
2 Magnetic field on the lattice
Let us consider the case of a constant external magnetic field , that is pointing in the direction. In the continuum such a magnetic field can be realized by, e.g., the following vector potential,
Any other vector potential satisfying corresponds to the same physical system, and is connected to the above choice by an appropriate gauge transformation.
It is well known that in a finite box with periodic boundary conditions the magnetic flux cannot be arbitrary, but is quantized in terms of the area of the system in the plane orthogonal to the external field 'tHooft:1979uj (); AlHashimi:2008hr (). This leads to the quantization condition,
where is the charge of the particle. In turn, on the lattice the area is also quantized as , where is the lattice spacing (we restrict the discussion to isotropic lattices) and is the number of lattice points in the direction . This implies that the lattice discretization also imposes an upper bound on the magnetic flux. To see this explicitly, let us write down how the continuum vector potential (1) can be represented by complex phases that multiply the links of the lattice,
where the sites are labeled by integers , with . Constant magnetic background fields were first used in lattice studies of nucleon magnetic moments at zero temperature Martinelli:1982cb (); Bernard:1982yu (); Zhou:2002km (). At finite temperature this approach was first realized in D'Elia:2010nq (); Roberts:2010cz ().
On the lattice the quantization condition (2) thus takes the same form as in the continuum with the area discretized in terms of the lattice spacing,
In this formulation we have periodic boundary conditions in all spatial directions and the magnetic flux going through any plaquette in the plane is constant. Furthermore, this implementation of the magnetic field is periodic in the flux quantum with a period of , so its value is effectively constrained to . This prescription is discussed in more detail in appendix A. We note that the periodicity of the field in implies a stronger constraint for the flux,
where the correspondence between the implementation and the actual value of is unambiguous. Namely, at larger values of the periodicity is expected to introduce saturation effects, like it was observed in D'Elia:2011zu (); Bruckmann:2011zx (). The largest possible magnetic field is therefore .
On the lattice the temperature of the system is given by the inverse temporal extension as . It is therefore clear from equation (4) that the minimal value of the magnetic field is . Thus, to increase the maximal field one has to decrease , and to decrease the minimal magnetic field one has to increase . Taking these considerations into account, with reasonable lattice spacings and lattice extensions, the lattice magnetic field covers the region GeV.
We remark that if there are particles with different charges in the system, then the quantization condition for has to be fulfilled for the greatest common divisor – in our case this is the down quark charge, in equation (4) with the charge of the electron. This will then determine the minimal field. Fortunately in nature the ratio of quark charges is a small natural number so that the up and down quarks can be studied together. We note furthermore that the above implementation of the magnetic field leads to no sign problem, in contrast to a finite chemical potential, or a (Minkowskian) electric field.
3 Observables at finite
Let us consider the staggered partition function with three flavors and . Each quark flavor has to be treated separately since the charges/masses are different: , and we assume . The partition function reads, after taking the fourth roots of the fermion determinants,
where the fermion matrix is . (Here we do not address problems arising from the rooting trick Durr:2005ax ()). The dependence on the chemical potential is only made explicit to define derivatives of the partition function with respect to , see equation (8), and later we set all chemical potentials to zero. Since we are only concerned with a constant external field, the dynamics of the field introduced above does not have to be taken into account; in the gauge sector we only have the kinetic term with inverse gauge coupling .
To study thermodynamics in a nonzero external field we analyze the chiral condensates and chiral susceptibilities for the light flavors ,
and the strange quark number susceptibility,
where we defined the spatial volume of the system as with . The condensate for a particular flavor will be denoted in the following by the first letter of the flavor name, e.g. .
To take the continuum limit, the renormalization of these observables has to be carried out. The logarithm of the partition function (i.e. the free energy) at contains additive divergences of the forms , and Leutwyler:1992yt (). In section 4 we will show – based on the behavior of the beta function measured at zero temperature – that there are no additional -dependent divergences. Therefore the additive divergences of the observables derived from the free energy can be eliminated by subtracting the , contribution. In the chiral quantities there are also multiplicative divergences caused by the derivative with respect to the quark mass. To eliminate this multiplicative divergence in the chiral condensate (susceptibility), we multiply by the first (second) power of the bare quark mass Endrodi:2011gv (). Finally, to obtain a dimensionless combination we divide by the fourth power of the pion mass ,
Note that this procedure leads to a renormalized condensate that, for , is zero at and approaches a negative value as is increased.
Considering the strange quark number susceptibility, needs no renormalization (neither at nor at ) since it is connected to a conserved current.
4 Renormalization at finite
We expect that a nonzero external magnetic field does not introduce new divergences in the free energy density, since the external field – just like a chemical potential – is coupled to the current which is conserved. This expectation is also supported by the fact that in the presence of the external field there are no additional, divergent Feynman-diagrams due to the absence of internal photon lines.222Consider e.g. the gluon self-energy diagram (with one quark loop) which is – in a gauge invariant regularization – logarithmically divergent. The coupling to the external magnetic field in the lowest order in is given by two external photon legs attached to the quark loop. This diagram is clearly finite since it contains two extra quark propagators. Moreover, the vacuum energy was also calculated in the effective potential approach Salam1975203 () and its divergent part was found to be independent of . The absence (or presence) of -dependent divergences in the free energy density is closely related to the non-renormalization (or renormalization) of itself. In fact, in a gauge invariant renormalization scheme the product needs no renormalization because of the Ward-Takahashi identity (see appendix B). In our case the magnetic field always appears in the combination and, therefore, we expect that it is not subject to renormalization333Note that for a dynamical theory, would appear separately in the photon gauge action. This is, however, not the case for the present study.. We now check these expectations numerically.
From the point of view of renormalization theory it may be instructive to draw a parallel between the magnetic field and the quark mass. If the magnetic field were to induce new divergences (like the mass does) then itself (like ) would be subject to renormalization as with a corresponding renormalization constant. However, in the lattice approach the magnetic field – unlike the mass – fulfills a quantization condition, as in equation (4). Therefore the renormalization of can only amount to a shift in the lattice spacing , such that = is satisfied. This implies that the lattice scale has to change if the renormalization of the magnetic field is nontrivial, i.e. if there are -dependent divergences. This is expected since the lattice scale is determined by the beta-function of the theory which is given in terms of the renormalization scale-dependent, divergent Feynman-diagrams (see e.g. Peskin:1995ev ()). For the magnetic field however, due to the quantization condition, the only possible effect of such divergent diagrams is to alter the lattice scale.
Therefore we propose to measure a physical quantity at as a function of the magnetic flux for different lattice spacings. We take the lattice scale and the line of constant physics (LCP) which are measured at (see section 5), and assume that they are also valid at . We use the scale to determine the magnetic field from the flux according to equation (4) and the quantity in physical units. Then we compare the running of with the magnetic field for different lattice spacings. If this quantity has a meaningful continuum limit,
i.e. if for small enough lattice spacings the dependence of on is suppressed, then our assumption was valid and the lattice scale and the LCP are also correct for . In view of the discussion in the previous paragraph, this suggests that there are no -dependent divergences. In the opposite case the lattice scale does depend on , which in turn would indicate the presence of -dependent divergences.
For physical quantities we choose the charged pion mass and the Sommer parameter . For their definition and measurement details see, e.g. Aoki:2009sc (). For the charged pion mass we expect a strong dependence on the magnetic field, in the form Chernodub:2010qx (),
which can be deduced from the dispersion relation for spin-1 mesons444We note that this expression ought to receive corrections for large due to pair production.. On the other hand, for the Sommer parameter, which is defined using the potential between static color charges, we expect the dependence on to be suppressed.
In figure 1 the mass of the charged pion and the Sommer parameter are plotted as functions of the external field for various lattice spacings, at , , and fm (for the Sommer parameter the coarsest lattice spacing is not shown since here has large systematic errors). The lattice geometries and simulation parameters for these runs are tabulated in appendix C. At we used the measurements presented in Aoki:2009sc (). We observe a nice scaling with for both quantities, in the region , where we have data for the three finest lattices. Results are consistent with a constant behavior for the Sommer parameter which indicates that the lattice spacing is not modified by the external field beyond our statistical accuracy.
Data for the charged pion mass are also as expected and agree with the analytic prediction (11) within . For large we see deviations from the continuum scaling only for the coarsest lattice which is most probably due to lattice artefacts stemming from the periodicity of the lattice magnetic field (3). Based on theoretical arguments (see appendix B) and on these observations we conclude that it is safe to use the lattice scale and LCP at nonzero external fields, and to exclude the possibility of -dependent divergences.
5 Simulation details
We use the tree-level improved Symanzik gauge action and stout smeared staggered fermions; details about the action can be found in Aoki:2005vt (). We generate lattice configurations both at and with an exact RHMC algorithm, for various values of the gauge coupling and the magnetic flux. (To discretize the external magnetic field, the smeared links are multiplied by the links of equation (3).) For our zero temperature measurements we simulate , and lattices, while for the finite temperature runs we have lattice configurations with and . Finite volume effects are studied on the ensemble using sets of , and lattices. The masses of the up, down and strange quarks are set to their physical values along the line of constant physics (LCP) by fixing the ratios and to their experimental values. The lattice spacing is determined by . Details of the determination of the LCP and the lattice scale can be found in, e.g. Borsanyi:2010cj (). Based on the reasoning presented in section 4, we use the lattice spacing measurements at and to set the scale also at and . The nonzero value of the magnetic field may of course modify e.g. the pion decay constant, just as the temperature can, but this is not important from the aspect of matching the lattice quantities at , to their experimental values (which are also measured at and ).
Altogether we generated several hundred to few thousand thermalized trajectories for each and (see list of simulation parameters in appendix C), and performed the measurements on every fifth one to decrease autocorrelations. The observables presented in section 3 were measured using the random estimator method, with 40 random vectors. The production of configurations and the measurements were performed on CUDA-capable GPU clusters at the Eötvös University in Budapest and on the Bluegene/P at FZ Jülich.
We mention here that the staggered formulation of fermions introduces lattice artefacts due to the splitting of hadron states into multiplets with different masses Ishizuka:1993mt (). We keep the lowest lying pion state at the physical pion mass, while the other members of the pion multiplet are heavier. In the continuum limit this mass splitting between the tastes vanishes. However, at finite lattice spacing it can distort thermodynamic quantities. To reduce this splitting we apply stout smearing in the fermionic action, which is known to significantly reduce taste symmetry violation Borsanyi:2010cj ().
6 Analysis details
To study the -dependence of the observables of section 3 we scan a wide interval in both the temperature and the flux quantum . The latter is proportional to , see equation (4), so, since the transition spreads over a wide temperature region, the physical magnetic field also changes by up to a factor of two along an line between MeV and MeV. To correct for this change one can simulate at parameters tuned such that the physical magnetic field remains the same. However, since cannot be varied continuously, here we follow a slightly different approach.
We measure our observables along a grid of points in the plane, as depicted in figure 2. The simulation points are denoted by the blue crosses, while the curves are shown by the red dashed lines. To perform the interpolation of the measurements along these lines in a systematic and effective way, we fit a two-dimensional spline function to the data points. A similar approach is described in Endrodi:2010ai () for the fitting of the gradient of a two-dimensional function. In figure 3 we show the observables as functions of and for our lattices. We obtain reliable results with good fit qualities; being in the range .
We perform simulations over the same physical temperature and magnetic field range for two smaller lattice spacings at and , with very similar values for the spline fits as above. We use these three lattice spacings (around they correspond roughly to , 0.15 and 0.12 fm) to extrapolate our results to the continuum limit.
7 Behavior of the condensate
We remark already at this point that the pseudocritical temperature – as probably best visible in the upper right panel of figure 3 for the chiral susceptibility – apparently decreases with increasing , thereby contradicting a vast number of model calculations, see the summary given in the introduction. Furthermore this observation also disagrees with the lattice result of D'Elia:2010nq (). First of all, to check our simulation code we reproduced the results of D'Elia:2010nq () at a couple of points, see appendix D. Since we find a perfect agreement, we conclude that we are left with three possible reasons for the discrepancy. First, the lattice spacing of D'Elia:2010nq () is larger, fm, and also an unimproved action is used, so lattice discretization errors may be significant. Second, the present study uses flavors as opposed to the of D'Elia:2010nq (), and the pseudocritical temperature is known to depend on the number of flavors Karsch:2000kv (), which may also introduce systematic differences in the dependence on the external field. Third, the quark masses of D'Elia:2010nq () are larger than in the present study, which can also cause drastic changes in thermodynamics – for example the nature of the transition at depends very strongly (and non-monotonically) on the quark masses.
On closer inspection, the differences between our results and those of D'Elia:2010nq () can actually be traced back to the behavior of the chiral condensate as a function of for a given temperature. While the authors of D'Elia:2010nq () observed that at any temperature the condensate increases with , we find that this dependence is more complex, see the left panel of figure 4 for our results. At MeV, which is just above the zero-field pseudocritical temperature, the bare condensate decreases by a factor of 2 between and . As the temperature is reduced the function starts to develop a maximum, clearly visible for MeV and MeV. This non-monotonic behavior is not due to the saturation effects caused by the periodic implementation of the magnetic field on the lattice, since this maximum is located at very different values of for temperatures differing only by a few percent. Furthermore, for high temperatures the decrease is already visible at which is in the first 5 percent of the period, even for the up quark. To better illustrate this effect and to show that renormalization and conversion from to does not change the picture qualitatively, in the right panel of figure 4 we plot the value of the external field where the renormalized chiral condensate takes its maximum, as a function of the temperature. At high temperatures this maximum is located at , while below MeV it shifts to a nonzero magnetic field, in accordance with the left panel of figure 4. As already mentioned in the introduction, the possibility of such a decrease in the condensate with was also raised in low energy model calculations Miransky:2002rp ().
We summarize our findings as a) the dependence of the condensate on the external field is non-monotonic and varies strongly with temperature, and b) as a result the pseudocritical temperature shifts to lower values at large as compared to the case. The latter observation is supported by a similar dependence deduced from the chiral susceptibility or the strange quark number susceptibility, see section 9. We investigate the reason for this behavior further by increasing our light quark masses up to the physical strange quark mass, studying the theory. As a first approximation we apply the same lattice scale and line of constant physics as was used for the flavor analysis. This clearly introduces a systematic error, but most probably does not affect the qualitative behavior. In figure 5 we show the results for the chiral condensate and susceptibility (left panels), compared to the data (right panels). In the upper panels we plot the unrenormalized chiral susceptibility for the case of a vanishing external field (red bands) and a large field of GeV (blue bands). Furthermore, in order to see the change in the condensate due to the presence of the external field, we plot the difference between the condensate at GeV and at in the lower panels. For the case of we plot the slice of the 2-dimensional surfaces we obtained as described in section 6, while for we fit the data to a simple spline function (in the latter analysis we keep the physical value of the magnetic field constant by tuning as a function of to keep fixed).
As is clearly visible in the lower left panel of figure 5, the magnetic field reduces the chiral condensate for temperatures MeV, thus pushing the inflection point of the condensate towards the left and causing a decrease in . This decrease is also visible from the behavior of the corresponding susceptibility, shown in the upper left panel of the figure. On the other hand, for larger quark masses the situation drastically changes: the condensate increases with the magnetic field for all temperatures (see the lower right panel of figure 5), similarly as was observed in D'Elia:2010nq (). Moreover, there is no clear change in : the chiral susceptibility (upper right panel) gives consistent pseudocritical temperatures for both and GeV. This observation supports our explanation number three, namely that the difference regarding the change in between the present work and the study of D'Elia:2010nq () stems at least partially from the larger-than-physical quark masses of the latter.555Note that while the mass of the Goldstone pion of D'Elia:2010nq () is below 200 MeV, due to the larger taste splitting the higher lying pion tastes may have a larger impact on the response to the magnetic field.
To illustrate the behavior of the condensate from yet another aspect, we show in figure 6 the contour plot of the renormalized chiral condensate as a function of and . The color of the curves encodes the value of the condensate, ranging from to (blue towards red) in steps of . A similar plot of the results of D'Elia:2010nq () would consist of curves having positive slopes, indicating that each point of moves towards the right as a result of a finite . Here we find that for example the curve clearly has segments with negative slope, which once again reflects the complex behavior of the condensate as a function of and .
8 Nature of the transition - finite size effects
Here we address the question of how the strength of the transition changes as the external field is switched on. At the transition is known to be a broad crossover Aoki:2006we (), where the approximate order parameters like the chiral condensate change smoothly with the temperature, and no finite volume scaling is visible in the observables. Furthermore, the crossover nature of the transition implies that – as we will also observe, see figure 9 – different observables give different pseudocritical temperatures Aoki:2006br ().
As can be seen from figure 3, not just the transition temperature changes with , but also the shapes of our observables as functions of are altered by a finite magnetic field. Specifically, we find that the maximum value of the chiral susceptibility increases with , which may suggest the transition to become stronger for large magnetic fields, as was also reported in D'Elia:2010nq (). To properly determine the nature of the transition we search for finite volume scaling in our observables.
To this end we perform simulations at our largest magnetic field on the lattices with and . The largest lattice in the transition region corresponds to a box of linear size fm. Here we keep fixed (and not itself) as we are only interested in differences between the various volumes. In figure 7 the results for the chiral susceptibility (left panel) and for the chiral condensate (right panel) are shown as functions of the temperature for . The figure shows that our results agree within statistical errors with the and data, indicating that finite size errors are small, compared to statistical errors. This observation also implies that the transition at this high magnetic field is still an analytic crossover.
To further study how the strength of the transition changes we investigate the width of the renormalized chiral susceptibility. In figure 8 we plot the susceptibility divided by its maximum value as a function of for three different values of the magnetic field for the lattices. We find that, although the height of the peak in grows significantly (by almost a factor of 2 between and the largest , see also upper left panel of figure 5), the width of the peak is only mildly affected by the magnetic field. In particular, the width of the peak at half maximum decreases from MeV to MeV as the external field is increased from zero to . We find a very similar behavior on the and lattices. From this analysis our final conclusions are that the width of the transition decreases only mildly with increasing magnetic field, and as the finite size scaling analysis has shown, the transition remains an analytic crossover at least up to GeV.
9 The phase diagram
Finally, using the fitted two-dimensional surfaces of section 6, we study the observables as functions of the temperature, along the lines of constant magnetic field. In particular we analyze the renormalized chiral susceptibility , the renormalized chiral condensate and the strange quark number susceptibility . For the latter two observables we determine the pseudocritical temperature as the inflection points of the curves, while for the former we calculate the position of the maximum value of the observable. The results are shown in figure 9.
To carry out the continuum extrapolation, we fit the results for for all three lattice spacings ( and ) together with an -dependent polynomial function of order four of the form . This ensures the scaling of the final results with . We obtain indicating good fit qualities. In order not to make the plots overcrowded, we only show error bars for the continuum curves. The error coming from the continuum extrapolation is estimated to be 2 MeV and is added to the statistical error in quadrature. The error in the lattice scale determination Borsanyi:2010cj () propagates in the function and amounts to an additional MeV systematic error, which is not added to the errors for figure 9 since we find that it does not influence the shape of the curves.
As is clearly visible in figure 9, all three observables show that the pseudocritical temperature decreases with growing external field . Preliminary results for the Polyakov loop at one lattice spacing show a very similar decrease in , see appendix E. We observe that the strange susceptibility (which can be viewed as a quantity signaling the deconfinement transition) is the observable most sensitive to the external field. changes most drastically in this case, by almost MeV between and . We note that our results at are all consistent with earlier determinations of the pseudocritical temperature where the stout smeared staggered lattice action was used Aoki:2006br (); Aoki:2009sc (); Borsanyi:2010bp (). We mention that one would expect effects to become more pronounced as the magnetic field grows. However, the numerical data of the pseudocritical temperatures seem to scale well, even up to our maximum value of .
In this paper we studied the finite temperature transition of QCD in the presence of external (electro)magnetic fields via lattice simulations at physical quark masses. The extrapolation to the continuum limit is carried out, and finite size effects are under control. The results are relevant for the description of both the evolution of the early universe and of noncentral heavy ion collisions.
We obtained the phase diagram of QCD in the plane using three observables in the phenomenologically interesting region of . Performing a finite volume scaling study we found that the transition remains an analytic crossover up to our largest magnetic fields, with the transition width decreasing only mildly. This rules out the existence of a critical endpoint in the phase diagram below . Moreover, our results indicate that the transition temperature significantly decreases with increasing . This result contradicts several model calculations present in the literature which predict an increase in as grows (see the summary in section 1). We presented indications that the response of to the external field can be traced back to the behavior of the chiral condensate as a function of and . We showed that this behavior is more complex than is predicted by most model calculations (where the condensate increases with for any temperature), and that it depends very strongly on the quark masses.
We summarize our results in figure 10, which shows the QCD phase diagram in the plane as defined using the renormalized chiral condensate and the strange quark number susceptibility in the continuum limit. By comparing our magnetic fields to the maximal fields that may be produced in noncentral heavy ion collisions we conclude that the decrease in is negligible for RHIC and may be up to MeV for the LHC. Moreover, the effect grows with the magnetic field, exceeding 20 for at . This may have a significant impact on the description of the QCD transition during the evolution of the early universe.
This work has been supported by DFG grants SFB-TR 55, FO 502/1-2 and BR 2872/4-2, the EU grants (FP7/2007-2013)/ERC no 208740 and PITN-GA-2009-238353 (ITN STRONGnet). Computations were carried out on the GPU cluster at the Eötvös University in Budapest and on the Bluegene/P at FZ Jülich. We thank Ferenc Niedermayer for useful discussions, interesting ideas and for careful reading of the manuscript. G. E. would like to thank Massimo D’Elia, Swagato Mukherjee, Dániel Nógrádi, Tamás Kovács and Igor Shovkovy for useful discussions.
Appendix A Lattice vector potential and periodic boundary conditions
In this appendix we show that the lattice prescription for the links, as in equation (3), is indeed equivalent to the continuum vector potential up to a local gauge transformation. The direct lattice discretized version of the continuum vector potential (1) can be written as
The periodic boundary conditions are now only satisfied up to a local gauge transformation (transition function),
where is the abelian gauge field and the gauge transformation that acts in space. However, on the lattice it is more convenient to have exactly periodic boundary conditions. Hence we perform the inverse gauge transformation on the last -slice of the lattice, which changes fermions as
and the links in both the and directions as
resulting in periodic boundary conditions and the “twisted” links that we presented in equation (3).
Appendix B Renormalization properties of from gauge invariance
QCD with an external magnetic field has a local gauge invariance. Let be the bare quark field and and the bare external electromagnetic field and electromagnetic coupling. The renormalization of these are given as,
If both the regularization and the renormalization prescriptions are gauge invariant, then so are the renormalization constants and . The gauge transformation for the bare and renormalized quark fields is of the form
which shows that , due to the gauge invariance of . The same transformations for the external electromagnetic field are then
Dividing the second equation by and equating it to the first we obtain (inserting ),
This is a well known result in QED which therefore also applies for the case of QCD with an external magnetic field. Since is the particular combination which renormalizes the product , our conclusion is that does not need renormalization.
Appendix C Simulation parameters
In this appendix we tabulate the simulation parameters for the and the runs. The number of thermalized trajectories generated for each set of parameters ranges from several hundred to a few thousand.
|3.45, 3.465, 3.48, 3.488, 3.492, 3.495, 3.497, 3.5,|
|3.505, 3.507, 3.51, 3.514, 3.518, 3.525, 3.54|
|3.45, 3.465, 3.48, 3.495, 3.51, 3.525,|
|3.54, 3.555, 3.57, 3.585, 3.6, 3.625|
|0, 1, 2, 3, 4, 8, 16, 24, 32, 50, 70|
|3.465, 3.48, 3.488, 3.492, 3.495, 3.497, 3.5, 3.503,|
|3.505, 3.507, 3.51, 3.514, 3.518, 3.525, 3.54|
|3.525, 3.55, 3.575, 3.6, 3.625, 3.64,|
|3.65, 3.675, 3.7, 3.725, 3.75, 3.775|
|0, 1, 2, 3, 4, 5, 7, 9, 13, 18, 29, 40|
|3.6, 3.625, 3.65, 3.675, 3.687, 3.7, 3.712, 3.725,|
|3.738, 3.75, 3.762, 3.775, 3.8, 3.825, 3.85, 3.875|
|0, 1, 2, 4, 8, 12, 16, 25, 34, 40|
|1, 2, 3, 4, 8, 12, 16, 20, 24, 32, 70|
|1, 2, 3, 4, 5, 6, 12, 24|
|1, 2, 3, 4, 6, 12|
Appendix D Code check
To check our code and simulation algorithm we reproduced the results of D'Elia:2010nq () at one temperature. We employ exactly the same simulation setup, i.e. we use the Wilson gauge action and flavors of unsmeared naive staggered quarks on a lattice. We measure the up and down quark condensates and the Polyakov loop (see definition in appendix E) at gauge coupling and mass . We plot the average condensate and the Polyakov loop in figure 11 (to conform to the notation of D'Elia:2010nq () we divide our Polyakov loops by 3). We observe that results for both and agree within statistical errors, with the exception of one point for where the values differ by , as expected for 8 points on statistical grounds. Therefore we confirm that there is no discrepancy between results from the two algorithms/implementations.
Appendix E Polyakov loop
We carry out the same analysis as presented in section 6 for the Polyakov loop,
We note that while the quark condensates and susceptibilities and the quark number susceptibility depend explicitly on the magnetic field, the Polyakov loop, as a purely gluonic operator, is only affected by the modified spatial links indirectly; its expectation value at is influenced by the magnetic factors of equation (3) appearing in the fermion determinant. To cancel the multiplicative divergences of , we define the renormalized Polyakov loop Aoki:2006br () using the static quark-antiquark potential
We measure this observable and perform the spline fitting, see figure 12 for the results. The inflection point moves to smaller temperatures as we increase the magnetic field. This behaviour is similar to the decreasing transition temperature observed for the condensate or susceptibility (c.f. figure 3).
- (1) T. Vachaspati, Magnetic fields from cosmological phase transitions, Phys. Lett. B265 (1991) 258–261.
- (2) R. C. Duncan and C. Thompson, Formation of very strongly magnetized neutron stars - implications for gamma-ray bursts, Astrophys. J. 392 (1992) L9.
- (3) V. Skokov, A. Y. Illarionov, and V. Toneev, Estimate of the magnetic field strength in heavy-ion collisions, Int. J. Mod. Phys. A24 (2009) 5925–5932, [arXiv:0907.1396].
- (4) D. E. Kharzeev, L. D. McLerran, and H. J. Warringa, The Effects of topological charge change in heavy ion collisions: ’Event by event P and CP violation’, Nucl.Phys. A803 (2008) 227–253, [arXiv:0711.0950].
- (5) K. Fukushima, D. E. Kharzeev, and H. J. Warringa, The Chiral Magnetic Effect, Phys. Rev. D78 (2008) 074033, [arXiv:0808.3382].
- (6) D. Kharzeev, Parity violation in hot QCD: Why it can happen, and how to look for it, Phys.Lett. B633 (2006) 260–264, [hep-ph/0406125].
- (7) STAR Collaboration Collaboration, S. A. Voloshin, Probe for the strong parity violation effects at RHIC with three particle correlations, Indian J.Phys. 85 (2011) 1103–1107, [arXiv:0806.0029].
- (8) STAR Collaboration Collaboration, B. Abelev et. al., Azimuthal Charged-Particle Correlations and Possible Local Strong Parity Violation, Phys.Rev.Lett. 103 (2009) 251601, [arXiv:0909.1739].
- (9) ALICE Collaboration Collaboration, I. Selyuzhenkov, Anisotropic flow and other collective phenomena measured in Pb-Pb collisions with ALICE at the LHC, arXiv:1111.1875.
- (10) F. Wang, Effects of Cluster Particle Correlations on Local Parity Violation Observables, Phys.Rev. C81 (2010) 064902, [arXiv:0911.1482].
- (11) B. Müller and A. Schäfer, Charge Fluctuations from the Chiral Magnetic Effect in Nuclear Collisions, Phys.Rev. C82 (2010) 057902, [arXiv:1009.1053].
- (12) V. Voronyuk, V. Toneev, W. Cassing, E. Bratkovskaya, V. Konchakovski, et. al., (Electro-)Magnetic field evolution in relativistic heavy-ion collisions, Phys.Rev. C83 (2011) 054911, [arXiv:1103.4239].
- (13) A. J. Mizher, M. N. 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].
- (14) E. S. Fraga and A. J. Mizher, Can a strong magnetic background modify the nature of the chiral transition in QCD?, Nucl. Phys. A820 (2009) 103c–106c, [arXiv:0810.3693].
- (15) R. Gatto and M. Ruggieri, Deconfinement and Chiral Symmetry Restoration in a Strong Magnetic Background, Phys. Rev. D83 (2011) 034016, [arXiv:1012.1291].
- (16) R. Gatto and M. Ruggieri, Dressed Polyakov loop and phase diagram of hot quark matter under magnetic field, Phys. Rev. D82 (2010) 054027, [arXiv:1007.0790].
- (17) A. Osipov, B. Hiller, A. Blin, and J. da Providencia, Dynamical chiral symmetry breaking by a magnetic field and multi-quark interactions, Phys.Lett. B650 (2007) 262–267, [hep-ph/0701090].
- (18) K. Kashiwa, Entanglement between chiral and deconfinement transitions under strong uniform magnetic background field, Phys. Rev. D83 (2011) 117901, [arXiv:1104.5167].
- (19) C. V. Johnson and A. Kundu, External Fields and Chiral Symmetry Breaking in the Sakai- Sugimoto Model, JHEP 12 (2008) 053, [arXiv:0803.0038].
- (20) S. Kanemura, H.-T. Sato, and H. Tochimura, Thermodynamic Gross-Neveu model under constant electromagnetic field, Nucl. Phys. B517 (1998) 567–598, [hep-ph/9707285].
- (21) K. G. Klimenko, Three-dimensional Gross-Neveu model at nonzero temperature and in an external magnetic field, Theor. Math. Phys. 90 (1992) 1–6.
- (22) J. Alexandre, K. Farakos, and G. Koutsoumbas, Magnetic catalysis in QED(3) at finite temperature: Beyond the constant mass approximation, Phys. Rev. D63 (2001) 065015, [hep-th/0010211].
- (23) N. Evans, T. Kalaydzhyan, K.-y. Kim, and I. Kirsch, Non-equilibrium physics at a holographic chiral phase transition, JHEP 1101 (2011) 050, [arXiv:1011.2519].
- (24) N. O. Agasian and S. M. Fedorov, Quark-hadron phase transition in a magnetic field, Phys. Lett. B663 (2008) 445–449, [arXiv:0803.3156].
- (25) F. Preis, A. Rebhan, and A. Schmitt, Inverse magnetic catalysis in dense holographic matter, JHEP 1103 (2011) 033, [arXiv:1012.4785].
- (26) P. Cea and L. Cosmai, Color dynamics in external fields, JHEP 08 (2005) 079, [hep-lat/0505007].
- (27) P. Cea and L. Cosmai, Abelian chromomagnetic fields and confinement, JHEP 0302 (2003) 031, [hep-lat/0204023].
- (28) P. Cea, L. Cosmai, and M. D’Elia, QCD dynamics in a constant chromomagnetic field, JHEP 0712 (2007) 097, [arXiv:0707.1149].
- (29) V. P. Gusynin, V. A. Miransky, and I. A. Shovkovy, Dimensional reduction and catalysis of dynamical symmetry breaking by a magnetic field, Nucl. Phys. B462 (1996) 249–290, [hep-ph/9509320].
- (30) S.-i. Nam and C.-W. Kao, Chiral restoration at finite T under the magnetic field with the meson-loop corrections, Phys. Rev. D83 (2011) 096009, [arXiv:1103.6057].
- (31) 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].
- (32) I. A. Shushpanov and A. V. Smilga, Quark condensate in a magnetic field, Phys. Lett. B402 (1997) 351–358, [hep-ph/9703201].
- (33) T. D. Cohen, D. A. McGady, and E. S. Werbos, The chiral condensate in a constant electromagnetic field, Phys. Rev. C76 (2007) 055201, [arXiv:0706.3208].
- (34) N. O. Agasian, Chiral thermodynamics in a magnetic field, Phys. Atom. Nucl. 64 (2001) 554–560, [hep-ph/0112341].
- (35) A. V. Zayakin, QCD Vacuum Properties in a Magnetic Field from AdS/CFT: Chiral Condensate and Goldstone Mass, JHEP 07 (2008) 116, [arXiv:0807.2917].
- (36) V. A. Miransky and I. A. Shovkovy, Magnetic catalysis and anisotropic confinement in QCD, Phys. Rev. D66 (2002) 045006, [hep-ph/0205348].
- (37) M. D’Elia, S. Mukherjee, and F. Sanfilippo, QCD Phase Transition in a Strong Magnetic Background, Phys. Rev. D82 (2010) 051501, [arXiv:1005.5365].
- (38) G. ’t Hooft, A Property of Electric and Magnetic Flux in Nonabelian Gauge Theories, Nucl. Phys. B153 (1979) 141.
- (39) M. H. Al-Hashimi and U. J. Wiese, Discrete Accidental Symmetry for a Particle in a Constant Magnetic Field on a Torus, Annals Phys. 324 (2009) 343–360, [arXiv:0807.0630].
- (40) G. Martinelli, G. Parisi, R. Petronzio, and F. Rapuano, The proton and neutron magnetic moments in lattice QCD, Phys. Lett. B116 (1982) 434.
- (41) C. W. Bernard, T. Draper, K. Olynyk, and M. Rushton, Lattice QCD calculation of some baryon magnetic moments, Phys. Rev. Lett. 49 (1982) 1076.
- (42) L. Zhou, F. X. Lee, W. Wilcox, and J. C. Christensen, Magnetic polarizability of hadrons from lattice QCD, Nucl. Phys. Proc. Suppl. 119 (2003) 272–274, [hep-lat/0209128].
- (43) D. S. Roberts, P. O. Bowman, W. Kamleh, and D. B. Leinweber, Wave Functions of the Proton Ground State in the Presence of a Uniform Background Magnetic Field in Lattice QCD, Phys. Rev. D83 (2011) 094504, [arXiv:1011.1975].
- (44) M. D’Elia and F. Negro, Chiral Properties of Strong Interactions in a Magnetic Background, Phys. Rev. D83 (2011) 114028, [arXiv:1103.2080].
- (45) F. Bruckmann and G. Endrődi, Dressed Wilson loops as dual condensates in response to magnetic and electric fields, Phys.Rev. D84 (2011) 074506, [arXiv:1104.5664].
- (46) Dürr, Stephan, Theoretical issues with staggered fermion simulations, PoS LAT2005 (2006) 021, [hep-lat/0509026].
- (47) H. Leutwyler and A. V. Smilga, Spectrum of Dirac operator and role of winding number in QCD, Phys. Rev. D46 (1992) 5607–5632.
- (48) G. Endrődi, Z. Fodor, S. D. Katz, and K. K. Szabó, The QCD phase diagram at nonzero quark density, JHEP 04 (2011) 001, [arXiv:1102.1356].
- (49) A. Salam and J. Strathdee, Transition electromagnetic fields in particle physics, Nuclear Physics B 90 (1975), no. 0 203–220.
- (50) M. E. Peskin and D. V. Schroeder, An Introduction to quantum field theory, . Section 12.2, paragraph “Computation of and ”. Reading, USA: Addison-Wesley (1995) 842 p.
- (51) Y. Aoki et. al., The QCD transition temperature: results with physical masses in the continuum limit II, JHEP 06 (2009) 088, [arXiv:0903.4155].
- (52) M. N. Chernodub, Superconductivity of QCD vacuum in strong magnetic field, Phys. Rev. D82 (2010) 085011, [arXiv:1008.1055].
- (53) Y. Aoki, Z. Fodor, S. D. Katz, and K. K. Szabó, The equation of state in lattice QCD: With physical quark masses towards the continuum limit, JHEP 01 (2006) 089, [hep-lat/0510084].
- (54) S. Borsányi et. al., The QCD equation of state with dynamical quarks, JHEP 11 (2010) 077, [arXiv:1007.2580].
- (55) N. Ishizuka, M. Fukugita, H. Mino, M. Okawa, and A. Ukawa, Operator dependence of hadron masses for Kogut-Susskind quarks on the lattice, Nucl. Phys. B411 (1994) 875–902.
- (56) G. Endrődi, Multidimensional spline integration of scattered data, Comput. Phys. Commun. 182 (2011) 1307–1314, [arXiv:1010.2952].
- (57) F. Karsch, E. Laermann, and A. Peikert, Quark mass and flavor dependence of the QCD phase transition, Nucl. Phys. B605 (2001) 579–599, [hep-lat/0012023].
- (58) Y. Aoki, G. Endrődi, Z. Fodor, S. D. Katz, and K. K. Szabó, The order of the quantum chromodynamics transition predicted by the standard model of particle physics, Nature 443 (2006) 675–678, [hep-lat/0611014].
- (59) Y. Aoki, Z. Fodor, S. D. Katz, and K. K. Szabó, The QCD transition temperature: Results with physical masses in the continuum limit, Phys. Lett. B643 (2006) 46–54, [hep-lat/0609068].
- (60) S. Borsányi et. al., Is there still any Tc mystery in lattice QCD? Results with physical masses in the continuum limit III, JHEP 09 (2010) 073, [arXiv:1005.3508].