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 QED3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_3$$\end{document}, 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 the chiral phase transition is apparently different from that of the deconfinement phase transition and in Nambu phase the fermion is confined only for small gauge-boson mass.


Introduction
The chiral and deconfinement phase transitions of nonperturbative 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 the 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 The two functions A( p 2 ) and B( p 2 ) in the above equation are related to the inverse fermion propagator The deconfinement phase transition is then related to the observation of the free particle and also the corresponding a e-mail: fenght@seu.edu.cn b e-mail: zonghs@chenwang.nju.edu.cn 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 as regards 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 nonperturbative features while being simpler. Three-dimensional quantum electrodynamics (QED 3 ) is just such a model which has many features similar to quantum chromodynamics (QCD), such as DCSB and confinement [2][3][4][5]. Moreover, its superrenormalization obviates the ultraviolet divergence which is present in QED 4 . For 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 [6][7][8] and high-T c cuprate superconductors [9][10][11][12].
The study of DCSB in QED 3 has been an active research subject for nearly 30 years, since Appelquist et al. [13] found that DCSB vanishes when the flavor of the massless fermions reaches a critical number N c ≈ 3.24. They gained 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 [14,15]. On the other hand, the achievement in the research of the mass singularity and confinement in QED 3 is caused by a paper of Maris [3] who found that the fermion is confined by the truncated DSE for the full fermion and boson propagators at N < N c where chiral symmetry is broken. This result might imply that the existence of confinement and DCSB depends on the same boundary conditions. Moreover, the authors of Refs. [2,16] 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 [17,18], or from a topological origin [19,20]. 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 is restored. Since DCSB and confinement are nonperturbative phenomena, both of them occur in the lowenergy region and might disappear with the increase of the boson mass. Therefore, it is very interesting to investigate whether or not both phase transitions occur at the same critical point in this case. One may argue that, with the introduction of the explicit photon mass term, the potential obtained from the full photon propagator will be attended at long distance, and the potential is no longer confining. Here, it is of heuristic interest to compare this case with that of QCD, where the quark interacts via the gluons. Lattice simulations illustrate that the gluon has a nonzero mass at zero temperature and zero density where, because of the occurrence of confinement, the free quarks cannot be observed [21]. This indicates that the nonzero mass of the boson does not destroy confinement. Based on the same spirit in this paper, we aim to 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.

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 nonperturbative effects. Just as mentioned in Sect. 1, the chiral symmetry is broken by the dynamical generation of the fermion mass (here N = 1).
Since the earlier work shows that the Schwinger function is suitable for the description of the fermion mass as regards the axiom of reflection positivity. Note that the Schwinger function with oscillatory behavior is not positive definite, which is an easily identifiable signal in (t) due to the pair of complex conjugate poles. This violates the axiom of reflection positivity. It follows from this that an oscillating behavior of the function describes a field with a complex mass spectrum and/or residues that are not positive. This is appropriate for particles that decay and forms the basis of the argument that such a propagator allows a state to exist only for a finite time before stabilization; i.e., that the propagator describes a confined fermion [22]. If one adopts the full boson propagator, the results of the 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 nonperturbative phenomena. The Schwinger function can be written as 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, if the system reveals a stable observable asymptotic state with a mass m for the fermion propagator, then In this way, the analysis of the 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 a powerful tool to study it, we shall use the coupled gap equations to calculate this function.

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 units of α, namely, we choose a kind of natural units in which α = 1. From Eqs. (2) and (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 [17,18] where (q 2 ) is the vacuum polarization for the gauge boson which is satisfied by the polarization tensor (11) and ζ is the gauge-boson mass which is acquired through the Higgs mechanism which happens when the gauge field interacts with a scalar field 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 as regards the Higgs mechanism in QED 3 can be found in Refs. [17,23]). 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 [4,5]. 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 transition and CPT in this Higgs model.

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 mass-less 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 the application of iterative methods, we can obtain A, B, and .

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][30], but none of them completely resolves the problem. However, the Ward-Takahashi identity, provides us an effectual tool to obtain a reasonable ansatz for the full vertex [31]. The portion of the dressed vertex which is free of kinematic singularities, i.e. the BC vertex, can be written as Since the numerical results obtained using the first part of the vertex coincide very well with earlier investigations [15], we choose this one as a suitable ansatz, to be used in our calculation. Following the procedure in the 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,

Numerical results
After solving the above coupled DSEs in the rainbow approximation by means of the iteration method, we can obtain the three functions 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 is 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 decreas-ing 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 gaugeboson 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, associated with the fermion propagator; thus the free particle can never be observed where the fermion is confined. On the increase of ζ , the oscillating behavior remains but it vanishes at another critical value ζ R dc ≈ 0.068 and around it 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.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 the 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 ζ BC 1 dc ≈ 0.038, which exhibits that the deconfinement phase transition occurs but here ψ ψ = 0. On the increase of ζ , the Schwinger function shows the real mass singularity of the propagator and chiral symmetry gets restored when the boson mass reaches ζ BC 1 c ≈ 0.071. Moreover, we investigate the critical behaviors with a range of N (< N c ). Our numerical results show that, whereas each of the two critical value changes with N , the qualitative conclusion that ζ c is larger than ζ dc remains.

Conclusions
The primary goal of this paper is to investigate the chiral and the 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 increase 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 understand 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 fact that the two critical numerical values are altered. The result indicates that, with the increasing gauge-boson mass in the chiral model, the occurrence of the deconfinement phase transition is apparently earlier than that of the chiral phase transition. Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited. Funded by SCOAP 3 / License Version CC BY 4.0.