Diracmode analysis for quark number density and its application for deconfinement transition
Abstract
The quark number density at finite imaginary chemical potential is investigated in the lattice QCD using the Diracmode expansion. We find the analytical formula of the quark number density in terms of the Polyakov loop in the large quark mass regime. On the other hand, in the small quark mass region, the quark number density is investigated by using the quenched lattice QCD simulation. The quark number density is found to strongly depend on the lowlying Dirac modes while its sign does not change. This result leads to that the quark number holonomy is not sensitive to the lowlying Dirac modes. We discuss the confinementdeconfinement transition from the property of the quark number density and the quark number holonomy.
Diracmode analysis for quark number density and its application for deconfinement transition
Takahiro M. Doi^{}^{}Speaker. T.M.D. is supported by the GrantinAid for JSPS fellows (No.15J02108) and the RIKEN Special Postdoctoral Researchers Program. email: takahiro.doi.gj@riken.jp and Kouji Kashiwa
1Theoretical Research Division, Nishina Center, RIKEN, Wako 3510198, Japan
\@textsuperscript2Yukawa Institute for Theoretical Physics,
Kyoto University, Kyoto 6068502, Japan

Abstract.
1 Introduction
One of the important task in nuclear and particle physics is the nonperturbative understanding of the phase structure of Quantum Chromodynamics (QCD) at finite temperature () and real chemical potential . However, the confinementdeconfinement transition has not been fully understood yet. For example, we cannot find any order parameter of the deconfinement transition in the case with dynamical quarks
The chiral and confinementdeconfinement transitions are key phenomena for this purpose, but the confinementdeconfinement transition is not yet fully understood comparing with the chiral transition. Although the chiral transition can be described by the spontaneous breaking of the chiral symmetry, but we can not find any classical orderparameters of the confinementdeconfinement transition in the presence of dynamical quarks.
The Polyakov loop, which respect the gaugeinvariant holonomy, is the exact orderparameter of the confinementdeconfinement transition in the infinite quark mass limit. However, it is no longer the order parameter in the presence of dynamical quarks. Although other candidates for the order parameter of the deconfinement transition have been proposed [1, 2, 3, 4, 5, 6], they are also not the exact order parameter. Therefore, we need some extension of ordinary determinations to clearly discuss and investigate the confinementdeconfinement transition in the system with dynamical quarks.
Recently, in Refs. [7, 8, 9], it is found that the topological change of QCD thermodynamics at finite can be used to determine the confinementdeconfinement transition, and based on the nontrivial freeenergy degeneracy, the quark number holonomy which is defined by the contour integral of the quarknumber susceptibility of has been proposed as the quantum orderparameter for the confinementdeconfinement transition. This argument is based on the analogy of the topological order discussed in Refs. [10] and QCD at [11] in the condensed matter physics. The quark number holonomy counts the gapped points of the quark number density along . As a results, it becomes nonzero/zero in the deconfined/confined phase. In particular, the quark number density at is important because it gives the property of the quark number holonomy. In order to investigate the quark number density at , we use the Diracmode expansion [12]. And we will see the behavior of the quark number density by removal of the lowlying Dirac modes.
This paper is organized as follows. In the section 2, we show the heavy quark mass expansion of the quark number density. In the section 3, we discuss the Diracmode expansion of the quark number density in both large and small quark mass regime. Section 4 is devoted to summary and discussions.
2 largemass expansion of quark number density
In this study, we consider the lattice QCD on the standard square lattice. We denote each sites as and linkvariables as . We impose the temporal periodic boundary condition for linkvariables to generate configurations in the quenched calculation to manifest the imaginarytime formalism.
On the lattice, the quark number density is defined as
(1) 
where denotes the functional trace and is taken over spinor and color indices. The Dirac operator in this article is taken as the WilsonDirac operator with quark mass and the chemical potential in the lattice unit as
(2) 
where is the identity matrix. The linkvariable operator is defined by the matrix element
(3) 
with and with . The chemical potential is the dimensionless on the lattice and we define .
We impose the temporal antiperiodicity and spatial periodicity for . To that end, we add a minus sign to the matrix element of the temporal linkvariable operator at the temporal boundary of :
(4) 
In this notation, the Polyakov loop is expressed as
(5) 
The minus sign stems from the additional minus on in Eq.(4).
In the heavy quark mass region, the quark number density (1) can be expressed by using the quark mass expansion as
(6) 
where we define the effective mass and the operator . In the case of the large quark mass, with smaller are dominant. The th order contribution has the terms
(7) 
and is a product of linkvariables. The examples of the paths of the products are shown in Fig. 1. Note that many of them become exactly zero because of Elitzur’s theorem [13]; only the gaugeinvariant terms corresponding to closed loops are nonzero. Moreover, spatially closed loops which do not wind the temporal length are canceled out each other and have no contribution to the quark number density in total in each order . Noting these important facts, it is confirmed that the th order contribution is constituted of the gaugeinvariant loops with the length which winds the temporal direction. In particular, the nonzero leading term in the expansion (6) is the th order term which relates to the Polyakov loop () and its complex conjugate (). The leading term, , to the quark number density at finite is written as
(8) 
In the heavy quarkmass expansion of the quark number density (6), there are higher order terms beyond the leading terms, which are the Polyakov loop and its conjugate. By similar analysis in the case of the leading term, the higherorder terms can be expressed in terms of the quantities which correspond to loops on the lattice. For example, the subleading terms include a term proportional to the quantity
(9) 
These terms correspond to the closed paths which wind the temporal length and bypasses the spatial direction. As another example, loops winding the temporal direction twice or more can be contribute to the expansion (6). For example, a loop winding the temporal direction twice
(10) 
is a possible contribution to Eq. (6) as the ()th order term.
3 Diracmode expansion of the quark number density
In the following, we consider the quark number density in terms of the Dirac eigenmode. In large quark mass region, we analytically investigate it in all order of the large quark mass expansion (6). In small quark mass region, we perform the quenched lattice QCD simulation.
The WilsonDirac eigenvalues are obtained from the eigenvalue equation as
(11) 
where is the WilsonDirac eigenstate. Considering the WilsonDirac mode expansion of the chiral condensate, the lowlying eigenmodes of the operators have dominant contribution to the chiral condensate known as BanksCasher relation [14, 15].
3.1 large quark mass region
We start the leading term to express it in terms of the WilsonDirac modes. The leading contribution of the quark number density in large quark mass region (8) is expressed by the Polyakov loop and its complex conjugate. It is already known that the Polyakov loop can be expressed in terms of the eigenmodes of the naiveDirac operator which corresponds to the case of [16, 17] and the WilsonDirac operator [18, 19]. In the following, we derive a different form of the Dirac spectral representation of the Polyakov loop using the operator on the square lattice with the normal nontwisted periodic boundary condition for linkvariables, in both temporal and spatial directions. We define a key quantity,
(12) 
This quantity is defined as the slightly changed quantity from the Polyakov loop by replacing a temporal linkvariable to the WilsonDirac operator . This quantity can be calculated as
(13) 
Thus, the quantity is proportional to . Here, other terms vanish because of the Elitzur’s theorem and the trace over the Dirac indecies. On the other hand, since in Eq. (12) is defined through the functional trace, it can be expressed in the basis of Dirac eigenmodes as
(14) 
The term arises because WilsonDirac operator is not normal due to the Wilson term and the completeness of the WilsonDirac eigenstates has the error:
(15) 
However, this error is controllable and can be ignored in close to the continuum limit.
Combining Eqs. (13) and (14), one can derive the relation between and the Dirac modes as
(16) 
with the error. From the formula (16), it is analytically found that the lowlying Dirac modes with have negligible contribution to the Polyakov loop because the eigenvalue plays as the damping factor. It is also numerically shown that there is no dominant contribution in the Dirac modes to the Polyakov loop [17]. Thus, one can find the Dirac spectrum representation of the quark number density in the leading order and the lowlying Dirac modes have little contribution to the quark number density.
The above discussion on the leading term can be applicable to the Dirac spectrum representation of the higher order terms. The detailed discussion is shown in our recent paper [20]. In the same way in the case of the leading term, all the terms in the expansion (6) can be expressed in terms of the WilsonDirac modes. Thus, it is analytically found that the quark number density does not depend on the density of the lowlying WilsonDirac modes in the allorder. However, this fact is only valid in the sufficiently large quark mass region since other contributions which can not be expressed by and can appear in the small region.
3.2 Small quark mass region
Next, we consider the small mass regime. In this regime, we perform the lattice QCD simulation to investigate the quark number density. In this study, we perform the quenched calculation with the ordinary plaquette action and then fermionic observables are evaluated by using the WilsonDirac operator (2) with the imaginary chemical potential . Our calculation is performed in both the confinement phase and the deconfinement phase. In the confinement phase, we consider lattice with and which corresponds to fm and MeV. In the deconfinement phase, we consider lattice with and which corresponds to fm and MeV. Both values of correspond to . In both cases, we set the quark mass as in the lattice unit, which is equivalent to the hopping parameter , for the calculation of the eigenmodes of the WilsonDirac operator in the small quark mass region [21].
We calculate the quark number density as
(17) 
This form trivially takes pure imaginary value up to the error. Each contribution, , to the quark number density of the Dirac mode with can be defined as
(18) 
and then the quark number density becomes
(19) 
In Fig. 3, is shown in the cases with and . We here only show results with one particular configuration.
To investigate the quark number density in terms of Dirac modes, we define the infrared (IR) cutted quark number density with the cutoff as
(20) 
In Fig. 4, we shown the result at , where the system is in the deconfinement phase. In general, we can perform the configuration average in the evaluation of Eq. (20). However, since the averaging well works after summing over all Diracmodes, it misses a physical meaning. Then, we show in one particular configuration.
Note that the sign of the quark number density does not changed by removing the lowlying Dirac modes while the absolute value of the quark number density seems to drastically change. This tendency has been found in almost all our configurations. This result means that the absolute value of the quark number density shares a same property in terms of Dirac modes with the chiral condensate, while its sign shares the property with the Polyakov loop.
Finally, the topological confinementdeconfinement transition can be discussed in terms of the Diracmode analysis. The orderparameter of the topological confinementdeconfinement transition can be expressed [8] as
(21) 
where is the dimensionless quark number density such as . It counts gapped points of the quark number density along direction and thus it becomes zero/nonzero in the confined/deconfined phase. The quark number holonomy (21) can be expressed as
(22) 
when the RW endpoint which is the endpoint of the RW transition line becomes the secondorder point at . In Eq. (22), characterizes and thus shares the same property about the Dirac modes with . The important point here is that the absolute value of does not have so much meaning even if it is nonzero, but its sign is important since the sign flipping at characterizes the gapped points along direction. From our quenched lattice QCD data, the sign of the quark number density are insensitive to the lowlying Diracmodes and this behavior is similar to the Polyakov loop.
4 Summary and discussion
In this paper, we have discussed properties of the quark number density at finite temperature () and imaginary chemical potential () by using the Diracmode expansion.
From the heavy quark mass expansion with the Dirac mode expansion, we found that lowlying Dirac modes do not dominantly contribute to the quark number density in all order of the heavy quark mass expansion.
In the small quark mass regime, we found that the absolute value of the quark number density strongly depends on lowlying Dirac modes, but its sign does not. Our result shows that the quark number holonomy is sensitive to the confinement properties of QCD and it is the good quantum order parameter for the confinementdeconfinement transition.
In order to avoid the incompleteness of the WilsonDirac modes, one can consider the hermitian WilsonDirac operator . In fact, the actual calculation is shown in our recent paper [20]. Moreover, the more detailed calculation on the higher order terms in the largemass expansion is also included in it.
References
 [1] E. Bilgici, F. Bruckmann, C. Gattringer, C. Hagen, Phys.Rev. D77, 094007 (2008), 0801.4051
 [2] C.S. Fischer, Phys.Rev.Lett. 103, 052003 (2009), 0904.2700
 [3] K. Kashiwa, H. Kouno, M. Yahiro, Phys.Rev. D80, 117901 (2009), 0908.1213
 [4] S. Benič, Phys.Rev. D88, 077501 (2013), 1305.6567
 [5] P.M. Lo, B. Friman, O. Kaczmarek, K. Redlich, C. Sasaki, Phys. Rev. D88, 074502 (2013), 1307.5958
 [6] T.M. Doi, K. Redlich, C. Sasaki, H. Suganuma, Phys. Rev. D92, 094004 (2015), 1505.05752
 [7] K. Kashiwa, A. Ohnishi, Phys. Lett. B750, 282 (2015), 1505.06799
 [8] K. Kashiwa, A. Ohnishi, Phys. Rev. D93, 116002 (2016), 1602.06037
 [9] K. Kashiwa, A. Ohnishi (2017), 1701.04953
 [10] X.G. Wen, Int.J.Mod.Phys. B4, 239 (1990)
 [11] M. Sato, Phys.Rev. D77, 045013 (2008), 0705.2476
 [12] S. Gongyo, T. Iritani, H. Suganuma, Phys. Rev. D86, 034510 (2012), 1202.4130
 [13] S. Elitzur, Phys. Rev. D12, 3978 (1975)
 [14] T. Banks, A. Casher, Nucl. Phys. B169, 103 (1980)
 [15] L. Giusti, M. Luscher, JHEP 03, 013 (2009), 0812.3638
 [16] H. Suganuma, T.M. Doi, T. Iritani, PTEP 2016, 013B06 (2016), 1404.6494
 [17] T.M. Doi, H. Suganuma, T. Iritani, Phys. Rev. D90, 094505 (2014), 1405.1289
 [18] H. Suganuma, T.M. Doi, K. Redlich, C. Sasaki (2016), 1610.02999
 [19] H. Suganuma, T.M. Doi, K. Redlich, C. Sasaki, EPJ Web Conf. 137, 04003 (2017), 1611.08742
 [20] T.M. Doi, K. Kashiwa (2017), 1706.00614
 [21] S. Aoki et al. (CPPACS), Phys. Rev. Lett. 84, 238 (2000), heplat/9904012