Chiral phase transition of three flavor QCD with nonzero magnetic field using standard staggered fermions
Abstract
Lattice simulations for (2+1)flavor QCD with external magnetic field demonstrated that the quark mass is one of the important parameters responsible for the (inverse) magnetic catalysis. We discuss the dependences of chiral condensates and susceptibilities, the Polyakov loop on the magnetic field and quark mass in three degenerate flavor QCD. The lattice simulations are performed using standard staggered fermions and the plaquette action with spatial sizes = 16 and 24 and a fixed temporal size = 4. The value of the quark masses are chosen such that the system undergoes a first order chiral phase transition and crossover with zero magnetic field. We find that in light mass regime, the quark chiral condensate undergoes magnetic catalysis in the whole temperature region and the phase transition tend to become stronger as the magnetic field increases. In crossover regime, deconfinement transition temperature is shifted by the magnetic field when quark mass is less than . The lattice cutoff effects are also discussed.
EPJ Web of Conferences \woctitleLattice2017 \selectlanguageenglish
1 Introduction
A strong magnetic field is expected to be produced in the early stage of the peripheral heavy ion collisions [1]. This interesting phenomenon drives a lot of theoretical studies since the external magnetic field affects the system through the electric charge of quarks especially phase and phase transition of QCD matter [2, 3].
From a view point of the field theory, the external magnetic field breaks several symmetries: 1) Rotational symmetry, 2) Time reversal symmetry and 3) Isospin symmetry. These breakings allow to mix channels which are prohibited by the symmetries of the original theory. For example, 1) allows to mix (axial)vector and (pseudo)scalar channels. Breaking of 3) can be obtained by introducing an isospin chemical potential but the external magnetic field also breaks 1) and 2) , thus the phase diagram is expected to be different from a system.
Recently QCD phase structure with magnetic field has been intensively studied on the lattice [2]. Earlier studies, that employed fourth rooted standard staggered fermions with found that the critical temperature increases as a function of magnitude of the magnetic field and the strength of the phase transition becomes stronger with larger magnetic field [4]. Afterwards by using stout staggered fermions with a physical pion mass a totally different result was found. There the inverse catalysis and the decreasing of with was observed [5, 6]. In the simulation with standard staggered fermions the rootmeansquare pion mass is much larger than the physical pion mass and this could be a cause of the observation of magnetic catalysis [7].
In this proceedings, we intend to study the influence of the value of quark mass to the behavior of as a function of the magnetic field , namely we investigate with QCD, which corresponds to a diagonal line for the extended Columbia plot (Right panel in Fig. 1). As a starting point of this project we will use the standard staggered fermions in our simulation. The numerical setups will be presented in the next section.
2 Setup
Here we introduce our numerical setup. We employ mass degenerated three flavor standard staggered fermions, , with the plaquette gauge action. Parameters for light mass regime are taken from [9], with which QCD has a first order phase transition at vanishing magnetic field . This simulation is performed with forth rooting technique as in previous studies and rational hybrid MonteCarlo algorithm.
The external U(1)\xspacemagnetic field is implemented in the following way.
The magnetic field only couples to quarks thus implementation is done just by changing SU(3)\xspacelinks to
. Here represents U(1)\xspacelinks
which contribute to the Dirac operator.
Finiteness of lattice size introduces an infrared cutoff to the U(1)\xspacefield [10].
Let us denote the lattice size and coordinate as ().
The external magnetic field in direction is described by the link variable of the U(1)\xspacefield and is expressed as follows,
{align}
u_x(n_x,n_y,n_z,n_t)&=
{exp[iq B N_x n_y] &(n_x= N_x1)
1 &(\textOtherwise)
\notag
u_y(n_x,n_y,n_z,n_t)&=exp[iq B n_x],
u_z(n_x,n_y,n_z,n_t)&=u_t(n_x,n_y,n_z,n_t)=1.\notag
where is the electric charge of each quark. Onevaluedness of the one particle wave function requires the Dirac quantization to be similar as in superconductors,
{align}
q^B = 2πNbNxNy ,
where is the number of magnetic flux through unit area for  plane
and .
The ultraviolet cutoff introduces also a periodicity of the magnetic field along with .
Namely, a range
represents an independent magnitude of the magnetic field similar to momentum in the Brillouin zone for crystals.
In order to examine effects on the phase transition coming from mass and the external field, we introduce a “threshold mass”. The continuum Dirac operator with a SU(3)\xspacevalued gluon field and external magnetic field is . By inserting the external magnetic field in the Landau gauge we obtain, {align} M_\textcon. ^†M_\textcon. &= D_μ^2 + (qxB)^2+m^2. This suggests that the term plays a similar role as the mass in this gauge. We introduce the “threshold mass” with by a dimensional analysis, {align} aeB = 2πNbNxNy—q— ≡a m_N_b^\textth. This mass is considered to be a threshold for magnetic field dominated regime to a mass dominated regime. We take in our study. By setting , and if taking larger than that value, the system thus may be considered to be dominated by the magnetic field rather than the quark mass. Details of our parameters are summarized in Tab.\xspace1.
range  # Conf.  Note  

5.128 – 5.160  0 – 56  First order for  
5.130 – 5.170  0 – 56  Crossover for  
5.10 – 5.65  0 – 56  Around for up quark  
5.35 – 5.65  0 – 56  Around for up quark  
5.35 – 5.85  0 – 56  Around for up quark 
3 Results
Here we summarize our preliminary results, which are listed in Tab.\xspace2. We measure the chiral condensate and the Polyakov loop and their susceptibilities. The critical temperature (critical ) for the chiral phase transition is determined by the susceptibility of chiral condensates in the light mass regime. The confinement/deconfinement transition is determined by one of the Polyakov loop in the heavy mass regime. In the light mass regime, we calculate the Binder cumulant [11] for the chiral condensate as a function of , {align} B_4 (β) = ⟨(δ¯ψψ)4⟩⟨(δ¯ψψ)2⟩2, where . The minimum of the Binder cumulant indicates the order of the phase transition: corresponds to a crossover, to a second order phase transition with the Ising universality class, to a first order phase transition [12]. However, contains forth order fluctuations, thus it is difficult to obtain accurate result compared to the susceptibility or the condensate itself.
along with  Chiral  Confinement  

Increase*  1st  1st  
Increase*  Crossover to critical  Crossover to critical  
Increase**  Not critical  Crossoverlike  
Increase**  Not critical  Crossoverlike  
Not clear  Not critical  Phase transitionlike*** 
Light mass regime
Here we show results for and . The former belongs to the first order regime for [9] and results are summarized in Fig.\xspace2 for the up quark. Down and strange quarks show similar behavior except for the absolute magnitude. The critical shifts as a function of monotonically, showing the magnetic catalysis instead of the inverse one.
The latter parameter corresponds to to the regime just above the critical mass for and the results are summarized in Fig.\xspace3. Standard fermions with , and case, it shows a crossoverlike behavior. This is consistent with the Binder cumulant (Bottom of Fig.\xspace3, black line). The dependence of shows similar tendencies in the case of .
The first order phase transition tends to become stronger as magnetic field increases in contrast to a model prediction in [13]. We can see from the chiral condensate that the critical temperature increases with the magnetic field.
Heavy mass regime
Here we show results with , and . In all cases here chiral condensates do not show any critical behavior as expected.
In Fig.\xspace4 and Fig.\xspace5 we show results of the Polyakov loop and its susceptibility for and , respectively. The confinement transition is slightly affected by the magnetic field. This is natural because gluons do not couple to the external magnetic field directly and effects of the magnetic field here come through dynamical but very heavy quarks (almost quenched). Fig.\xspace6 shows results of the Polyakov loop and its susceptibility for . For this mass, effects of the magnetic field are tiny. In this mass regime, the response of the Polyakov loop is similar to predicted by the PNJL model [14]. In addition the magnetic field increases the critical temperature in the same way as in the lighter mass regime except for the heaviest mass case.
By comparing the values of the threshold mass and the quark mass it seems that in the current simulation the QCD thermodynamics is more dominated by effects introduced by the quark mass.



4 Summary and discussion
In this work, we investigate the phase structure of QCD with three degenerate flavor and a U(1)\xspaceexternal magnetic field for various masses using standard staggered fermions. We have observed a tendency of strengthening of the phase transition in the light mass regime. Except for the heaviest case, the critical temperature increases with the magnetic field instead of decrease.
There are several issues which must be addressed in our forthcoming study. Firstly, we need to increase the statistics and improve resolution of to obtain more accurate results especially the determination of the order of phase transition from the Binder cumulant. Secondly, we find the increase of instead of decrease. That might come from cutoff effects (lattice artifacts) of our setup. This is because, for the case, the pion mass is relatively heavy (around 290 MeV [15, 16]), even in lighter quark mass regime. We have used parameters as in [9], which correspond to a first order chiral phase transition regime. However, the existence of a first order phase transition regime for massless three flavor QCD is not clear yet (see [12] and the references therein). In order to confirm or refine our results, we need to improve our lattice action to reduce the root mean square pion mass and lattice artifacts. This is important because the chiral features strongly depends on the lowlaying modes of the Dirac operator even in the case of vanishing magnetic field, which are tightly connected to discretization of the Dirac operator. Even more, lowlaying modes with magnetic field, namely the lowest Landau level, plays essential role when magnetic field is applied [17]. Thus, we have to use an improved action in the next step and clarify which phenomena are coming from continuum physics and which are just lattice artifacts.
Acknowledgement
We thank M. D’Elia for giving us exact numbers of their original work. AT thank the members of HotQCD collaboration for fruitful discussions and Y. Nakamura for providing information of current status for three flavor massless QCD in the critical regime. The simulations are performed at Fermilab Lattice Gauge Theory Computational Facility. The work of AT was supported in part by NSFC under grant no. 11535012.
References
 D.E. Kharzeev, L.D. McLerran, H.J. Warringa, Nucl. Phys. A803, 227 (2008), 0711.0950
 G. Endrodi, PoS LATTICE2014, 018 (2014), 1410.8028
 D.E. Kharzeev, J. Liao, S.A. Voloshin, G. Wang, Prog. Part. Nucl. Phys. 88, 1 (2016), 1511.04050
 M. D’Elia, S. Mukherjee, F. Sanfilippo, Physical Review D 82, 051501 (2010)
 G. Bali, F. Bruckmann, G. Endrődi, Z. Fodor, S. Katz, S. Krieg, A. Schäfer, K. Szabo, Journal of high energy physics 2012, 44 (2012)
 G. Bali, F. Bruckmann, G. Endrődi, Z. Fodor, S. Katz, A. Schäfer, Physical Review D 86, 071502 (2012)
 F. Bruckmann, G. Endrodi, T.G. Kovacs, JHEP 04, 112 (2013), 1303.3972
 F.R. Brown, F.P. Butler, H. Chen, N.H. Christ, Z.h. Dong, W. Schaffer, L.I. Unger, A. Vaccarino, Phys. Rev. Lett. 65, 2491 (1990)
 D. Smith, C. Schmidt, PoS LATTICE2011, 216 (2011), 1109.6729
 M. AlHashimi, U.J. Wiese, Annals of physics 324, 343 (2009)
 K. Binder, Physical Review Letters 47, 693 (1981)
 A. Bazavov, H.T. Ding, P. Hegde, F. Karsch, E. Laermann, S. Mukherjee, P. Petreczky, C. Schmidt, Phys. Rev. D95, 074505 (2017), 1701.03548
 N. Mueller, J.M. Pawlowski, Physical Review D 91, 116010 (2015)
 K. Fukushima, M. Ruggieri, R. Gatto, Phys. Rev. D81, 114031 (2010), 1003.0047
 F. Karsch, E. Laermann, C. Schmidt, Phys. Lett. B520, 41 (2001), heplat/0107020
 F. Karsch, C.R. Allton, S. Ejiri, S.J. Hands, O. Kaczmarek, E. Laermann, C. Schmidt, Nucl. Phys. Proc. Suppl. 129, 614 (2004), [,614(2003)], heplat/0309116
 T. Kojo, N. Su, Nuclear Physics A 931, 763 (2014)