Different critical points of chiral and deconfinement phase transitions in (2+1)-dimensional fermion-gauge interacting model

Based on the truncated Dyson-Schwinger equations for fermion and massive boson propagators in QED$_3$, the fermion chiral condensate and the mass singularities of the fermion propagator via the Schwinger function are investigated. It is shown that the critical point of chiral phase transition is apparently different from that of deconfinement phase transition and in Nambu phase the fermion is confined only for small gauge boson mass.


I. INTRODUCTION
The chiral and deconfinement phase transitions of nonperaturbative systems are important issues of continuous interests both theoretically and experimentally. Although the mechanism is unknown, the originally chiral symmetric system may undergo chiral phase transition (CPT) into a phase with dynamical chiral symmetry breaking (DCSB) which explains the origin of constituent-quark masses in QCD and underlies the success of chiral effective field theory [1,2]. In the chiral limit, the order parameter of CPT is defined via the fermion propagator .
(1) The two functions A(p 2 ) and B(p 2 ) in the above equation are related to the inverse fermion propagator S −1 (p) = iγ · pA(p 2 ) + B(p 2 ). ( The deconfinement phase transition is then related to the observation of the free particle and also the corresponding propagator. If the full fermion propagator has no mass singularity in the timelike region, it can never be on mass shell and the free particle can never be observed where the confinement happens [3]. Accordingly, the appearance of the mass singularity in the system directly implies deconfinement. So in this way we can learn the deconfinement phase transition from the analytic structure of the fermion propagator. To indicate DCSB and confinement, it is very suggestive to study some model that reveals the general nonperaturbative features while being simpler. Threedimensional quantum electrodynamics (QED 3 ) is just such a model which has many features similar to quantum * Email:fenght@seu.edu.cn † Email:zonghs@chenwang.nju.edu.cn chromodynamics (QCD), such as DCSB and confinement [2][3][4][5][6][7][8]. Moreover, its superrenormalization obviate the ultraviolet divergence which is present in QED 4 . Due to these reasons, it can serve as a toy model of QCD. In parallel with its relevance as a tool through which to develop insight into aspects of QCD, QED 3 is also found to be equivalent to the low-energy effective theories of strongly correlated electronic systems. Recently, QED 3 has been widely studied in graphene [9][10][11] and high-T c cuprate superconductors [12][13][14][15].
The study of DCSB in QED 3 has been an active subject near 30 years since Appelquist et al. found that DCSB vanishes when the flavor of massless fermions reaches a critical number N c ≈ 3.24 [16]. They gain this conclusion by solving the truncated Dyson-Schwinger equation (DSE) for the fermion propagator in the chiral limit. Later, extensive analytical and numerical investigations showed that the existence of DCSB in QED 3 remains the same after including higher order corrections to the DSE [17,18]. On the other hand, the achievement in research of the mass singularity and confinement in QED 3 is caused by a paper of P. Maris who found that the fermion is confined by the truncated DSE for the full fermion and boson propagators at N < N c [3] where chiral symmetry is broken. This result might imply that the existence of confinement and DCSB depend on the same boundary conditions. Moreover, the authors of Ref. [2,19] pointed out that restoration of chiral symmetry and deconfinement are coincident owing to an abrupt change in the analytic properties of the fermion propagator when a nonzero scalar self-energy becomes insupportable.
Nevertheless, the above result will be altered when the gauge boson acquires a finite mass ζ through the Higgs mechanism [20,21]. For a fixed N (< N c ) and with the increasing boson mass, the fermion chiral condensate falls and diminishes at a critical value ζ c (which, of course, depends on N ) and then chiral symmetry restores. Since DCSB and confinement are nonperaturbative phenomena, both of them occur in the low energy region and might disappear with the rise of boson mass. There-fore, it is very interesting to investigate whether or not both phase transitions occur at the same critical point in this case. In this paper, we will adopt the truncated DSEs for the full propagators to study the behaviors of the mass singularity and the fermion chiral condensate with a range of gauge boson mass and try to answer this question.

II. SCHWINGER FUNCTION
The Lagrangian for massless QED 3 in a general covariant gauge in Euclidean space can be written as where the 4-component spinor ψ is the massless fermion field, ξ is the gauge parameter. This system has chiral symmetry and the symmetry group is U (2). The original U (2) symmetry reduces to U (1) × U (1) when the massless fermion acquires a nonzero mass due to nonperaturbative effects. Just as mentioned in Sec. I, the chiral symmetry is broken by the dynamical generation of the fermion mass (here N = 1). If one adopts the full boson propagator, the results of Euclidean-time Schwinger function reveal that the fermion propagator has a complex mass singularity and thus corresponds to a nonphysical observable state [3] which means the appearance of confinement. On the contrary, if the Schwinger function exhibits a real mass singularity of the propagator, the fermion is observable and the fermion is not confined [24,25]. Therefore, we also adopt this method to analyze those nonperaturbative phenomena. The Schwinger function can be written as with M (p 2 ) = B(p 2 )/A(p 2 ). If there are two complex conjugate mass singularities m * = a ± ib associated with the fermion propagator, the function will show an oscillating behavior for large (Euclidean) t. However, the system reveals a stable observable asymptotic state with a mass m for the fermion propagator, then By this way, the analysis of mass singularity can be used to determine whether or not the fermion is confined. Since the Schwinger function is determined by the fermion propagator and the DSEs provide us an powerful tool to study it, we shall use the coupled gap equations to calculate this function.

III. TRUNCATED DSE
Now let us turn to the calculation of A(p 2 ) and B(p 2 ). These functions can be obtained by solving DSEs for the fermion propagator, where Γ ν (p, k) is the full fermion-photon vertex and q = p − k. The coupling constant α = e 2 has dimension one and provides us with a mass scale. For simplicity, in this paper temperature, mass and momentum are all measured in unit of α, namely, we choose a kind of natural units in which α = 1. Form Eq. (2) and Eq. (7), we obtain the equation satisfied by A(p 2 ) and B(p 2 ) Another involved function D σν (q) is the full gauge boson propagator which is given by [20,21] where Π(q 2 ) is the vacuum polarization for the gauge boson which is satisfied by the polarization tensor and ζ is the gauge boson mass which is acquired though Higgs mechanism which happens when the gauge field interacts with a scalar filed in the phase with spontaneous gauge symmetry breaking (Here, we adopt the massive boson propagator to investigate the oscillation behavior of Schwinger function in DCSB phase, more details about Higgs mechanism in QED 3 can be found in Ref. [20,22]). Using the relation between the vacuum polarization Π(q 2 ) and Π σν (q 2 ), we can obtain an equation for Π(q 2 ) which has an ultraviolet divergence. Fortunately, it is present only in the longitudinal part and is proportional to δ σν . This divergence can be removed by the projection operator and then we obtain a finite vacuum polarization [5,23]. Finally, we choose to work in the Landau gauge, since the Landau gauge is the most convenient and commonly used one. Once the fermion-boson vertex is known, we immediately obtain the truncated DSEs for the fermion propagator and then analyze the deconfinement and chiral phase transitions in this Higgs model.

A. Rainbow approximation
The simplest and most commonly used truncated scheme for the DSEs is the rainbow approximation, since it gives us rainbow diagrams in the fermion DSE and ladder diagrams in the Bethe-salpeter equation for the fermion-antifermion bound state amplitude. In the framework of this approximation, the coupled equations for massless fermion and massive boson propagators reduce to the three coupled equations for A(p 2 ), B(p 2 ) and Π(q 2 ), with G(k 2 ) = A 2 (k 2 )k 2 + B 2 (k 2 ). By application of iterative methods, we can obtain A, B and Π.

B. Improved scheme for DSE
To improve the truncated scheme for DSE, there are several attempts to determine the functional form for the full fermion-gauge-boson vertex [26][27][28][29], but none of them completely resolve the problem. However, the Ward-Takahashi identity provides us an effectual tool to obtain a reasonable ansatze for the full vertex [30]. The portion of the dressed vertex which is free of kinematic singularities, i.e. BC vertex, can be written as, Since the numerical results obtained using the first part of the vertex coincide very well with earlier investigations [18], we choose this one as a suitable ansatze to be used in our calculation. Following the procedure in rainbow approximation, we also obtain the three coupled equations for A(p 2 ), B(p 2 ) and Π(q 2 ) in the improved truncated scheme for DSEs,

IV. NUMERICAL RESULTS
After solving the above coupled DSEs in rainbow approximation by means of the iteration method, we can obtain the three function A, B, Π for the propagator and plot them in Fig. 1. From Fig. 1 it can be seen that A(p 2 ) increases with increasing momenta but almost equal to one at large p 2 . In the range of small momenta, it decreases but does not vanish when p 2 → 0. Both of the other two functions B(p 2 ) and Π(q 2 ) decrease at large momenta but their rates of decreasing are different. B(p 2 ) decreases as rapidly as ∼ 1/p 2 , while Π(q 2 ) decreases as rapidly as ∼ 1/ q 2 . In addition, all the three functions are constant in the infrared region. Thus, we can obtain the values of the corresponding functions A, B and Π at zero momenta, which, as functions of the gauge boson mass ζ, are also shown in Fig. 1 Then, substituting the obtained A and B into Eq. (4), we immediately obtain the behavior of the Schwinger function with nonzero boson mass which is shown in Fig.  2. At small ζ, the Schwinger function reveals its typical oscillating behavior which illustrates the conjugate mass singularities like m * = a ± ib m * ∼ 0.043 ± 0.063i at ζ = 0.01, (24) m * ∼ 0.023 ± 0.025i at ζ = 0.06, (25) associated with the fermion propagator and thus the free particle can never be observed where the fermion is confined. As the rise of ζ, the oscillating behavior remains but it vanishes at another critical value ζ R dc ≈ 0.068 and around which both of the propagators do not exhibit any singularity.
Beyond ζ R dc , the function ln[Ω(t)] ∼ −mt where the stable asymptotic state of the fermion is observable m ≈ 0.021 at ζ = 0.07, (26) m ≈ 0.0041 at ζ = 0.09 (27) and hence the deconfinement phase transition happens, but the DCSB remains. With the enlargement of ζ, the absolute slope of ln[Ω(t)] decreases and m disappears at ζ R c . To validate the difference between ζ c and ζ dc , we also give the behavior of the Schwinger function beyond rainbow approximation in Fig. 3. In the BC 1 truncated scheme for DSE, the oscillation of the Schwinger function only appears at small ζ, which denotes the existence of confinement, but it disappears at ζ BC1 dc ≈ 0.038, which exhibits that deconfinement phase transition occurs but here ψ ψ = 0. As the rise of ζ, the Schwinger function shows the real mass singularity of the propagator and chiral symmetry gets restored when the boson mass reaches ζ BC1 c ≈ 0.071.

V. CONCLUSIONS
The primary goal of this paper is to investigate chiral and deconfinement phase transition by application of an Abelian Higgs model through a continuum study of the Schwinger function. Based on the rainbow approximation of the truncated DSEs for the fermion propagator and numerical model calculations, we study the behavior of the Schwinger function and the fermion chiral condensate. It is found that, with the rise of the gauge boson mass, the vanishing point (ζ dc ) of the oscillation behavior of the Schwinger function is apparently less than that of the fermion chiral condensate and each of the propagators does not reveal any singularity near ζ dc . To make know the difference between the two critical points, we also work in an improved scheme for the truncated DSEs and show that the above conclusion remains despite the two critical numerical values alter. The result indicates that, with the increasing gauge boson mass in the chiral model, the occurrence of de-confinement phase transition is apparently earlier than that of chiral phase transition.

VI. ACKNOWLEDGEMENTS
We would like to thank Prof. Wei-min Sun and Guozhu Liu for their helpful discussions. This work was supported by the National Natural Science Foundation