Loading [MathJax]/jax/output/HTML-CSS/jax.js

Chiral symmetry restoration and deconfinement in the contact interaction model of quarks with parallel electric and magnetic fields

  • We study the impact of steady, homogeneous, and external parallel electric and magnetic field strengths (eEeB) on the chiral symmetry breaking-restoration and confinement-deconfinement phase transition. We also sketch the phase diagram of quantum chromodynamics (QCD) at a finite temperature T and in the presence of background fields. The unified formalism for this study is based on the Schwinger-Dyson equations, symmetry preserving vector-vector contact interaction model of quarks, and an optimal time regularization scheme. At T=0, in the purely magnetic case (i.e., eE0), we observe the well-known magnetic catalysis effect. However, in a pure electric field background (eB0), the electric field tends to restore the chiral symmetry and deconfinement above the pseudo-critical electric field eEχ,Cc. In the presence of both eE and eB, we determine the magnetic catalysis effect in the particular region where eB dominates over eE, whereas we observe the chiral inhibition (or electric chiral rotation) effect when eE overshadows eB. At finite T, in the pure electric field case, the phenomenon of inverse electric catalysis appears to exist in the proposed model. Conversely, for a pure magnetic field background, we observe the magnetic catalysis effect in the mean-field approximation and inverse magnetic catalysis with eB-dependent coupling. The combined effects of eE and eB on the pseudo-critical Tχ,Cc yields an inverse electromagnetic catalysis, with and without an eB-dependent effective coupling of the model. The findings of this study agree well with the already predicted results obtained via lattice simulations and other reliable effective models of QCD.
  • Dynamical chiral symmetry breaking and confinement are the two fundamental properties of non-perturbative quantum chromodynamics (QCD). At zero or low temperatures T, the fundamental degrees of freedom of quantum chromodynamics (QCD) are the low energy hadrons, whereas at high T, the interaction gets increasingly screened and thus becomes weak, causing hadrons to melt down to a new phase in which the predominant degrees of freedom are quarks and gluons. This phenomenon is referred to as the confinement–deconfinement phase transition. As the strength of the QCD interaction decreases with increasing T, only the current quark mass survives when the parameter T exceeds a critical value. This is termed as the chiral symmetry breaking–chiral symmetry restoration phase transition. Such a phase transition is expected to have occured in the early universe, a few microseconds after the Big Bang, and it is experimentally observed in heavy-ion collisions at the Large Hadron Collider (LHC) in CERN and Relativistic Heavy Ion Collider (RHIC) at Brook Heaven National Laboratory (BNL). In addition, when the hadronic matter is subjected to an external electromagnetic field background, it yields a significant impact on phase transition. It is well known that at T=0, in the pure magnetic case, the strong magnetic field tends to strengthen the formation of quark anti-quark condensate, and the system remains in the chiral symmetry broken phase, even at a high magnetic field strength eB. This phenomenon is known as magnetic catalysis (MC) [1-10]. It was verified in earlier studies that at a finite T, the pseudo-critical temperature Tχ,Cc of the chiral symmetry restoration and deconfinement increases with an increase in eB; hence, magnetic catalysis is also observed at finite T [1, 2, 5, 11, 12]. Recently, the lattice QCD simulation [13-15] predicted that at finite T, a magnetic field suppresses the formation of quark-antiquark condensates and tends to restore the chiral symmetry at approximately the pseudo-critical temperature Tχ,Cc. Consequently, Tχ,Cc decreases with an increase in eB, and such a phenomenon is known as the inverse magnetic catalysis (IMC). This phenomenon is validated and supported by effective models of low energy QCD [16-23], as well as in holographic QCD models [24].

    In a pure electric case, and at T=0, the situation is relatively different from that of a pure magnetic field background. The strong electric field suppresses the formation of a quark-antiquark condensate, and thus tends to restore the chiral symmetry, i.e., the electric field anti-screens the strong interaction. Such a phenomenon is known as the chiral electric inhibition effect [1, 2, 5, 11, 12, 25-28] or the chiral electric rotation effect [29]. The nature of chiral phase transitions is of the second order in the chiral limit but cross-over when the bare quark mass is considered. At finite temperature, it is well understood that the pseudo-critical temperature Tχ,Cc decreases with an increase in electric field strength eE, and this phenomenon is known as the inverse electric catalysis (IEC) [28, 30]. The study on the influence of electric fields on the chiral phase transitions, as well as the effect of the magnetic field, is equally important from theoretical and experimental perspectives. Experimentally, in heavy-ion collisions, the electric and magnetic fields are generated with the same order of magnitude (1018 to1020Gauss) [31-34] in the event-by-event collisions using Au + Au at RHIC-BNL and in a non-central heavy-ion collision of Pb + Pb in ALICE-LHC. Moreover, some interesting anomalous effects, such as the chiral magnetic [35, 36], chiral electric separation [37, 38], and particle polarization [39-41] effects, which may arise owing to the generation of vector and/or axial currents in the presence of strong electromagnetic fields, are also required for theoretical exploration. Recently, a special case of considering the electric field strength parallel to the magnetic field strength (eEeB) significantly focused on exploring the above-mentioned phenomenon in the effective models of QCD [28, 29, 42, 43]. A major reason behind this approach is that the parallel electric and magnetic fields play important and prominent roles in several heavy-ion collision experiments [44, 45].

    Considering the aforementioned facts and findings, the objective of this study is to elucidate the traits of the dynamical chiral symmetry breaking-restoration transition in the presence of parallel electric and magnetic fields. The unified framework of this study is based on Schwinger-Dyson equations (SDE) in the rainbow-ladder truncation, Landau gauge, symmetry preserving the confining vector-vector contact interaction model (CI) [46], and Schwinger optimal time regularization scheme [47]. We adopt the quark-antiquark condensate ˉqq as an order parameter for the chiral symmetry breaking-restoration, whereas for the confinement-deconfinement transition, we use the confinement length scale [22, 48]. It should be noted that the chiral symmetry restoration and deconfinement occur simultaneously in this model [22, 49].

    This article is organized as follows. In Sec. II, we present the general formalism and contact interaction model at zero temperature, in the absence of background fields. In Sec. III, we discuss the gap equation at zero temperature, in the presence of parallel electric and magnetic fields presented in Sec. IV. Next, in Sec. V, we present the phase diagram at finite temperature, in the presence of parallel electric and magnetic fields. Finally, in Sec. VI, we present the summary and perspectives of this study.

    We begin with the Schwinger-Dyson's equations (SDE) for a dressed-quark propagator Sf, which is given by

    S1f(p)=iγp+mf+d4k(2π)4g2Δμν(pk)λa2γμSf(k)λa2Γν(p,k),

    (1)

    the subscript f represents the two light quark flavors, i.e., up (u) and down (d) quarks, g is the coupling constant, and mf is the current quark mass, which can be set to zero in the chiral limit. λa represents the conventional Gell-Mann matrices, Γν is the dressed quark-gluon vertex, and Δμν depicts the gluon propagator.

    In literature, it is well known that the properties of low energy hadrons can be reproduced by assuming that the interaction among the quarks occurs not via a mass-less vector boson exchange but by the symmetry preserving four-fermions vector-vector CI with a finite gluon mass [46, 50-53]:

    g2Δμν(k)=δμν4παirm2Gδμναeff,

    (2)

    where αir=0.93π represents the infrared enhanced interaction strength parameter, and mG=800 MeV is the gluon mass scale [54]. In the CI model, for a small value of αir and a large value of the gluon mass scale mG, there must be a critical value αeff; above this critical value, the chiral symmetry is broken, and below this value, there is a lesser chance for the generation of the dynamical mass. The d-dimensional (arbitrary space-time dimensions) dependence of this effective coupling and its critical value on the chiral symmetry breaking, using an iterative method in the superstrong regime, has been comprehensively studied in Ref. [55].

    The CI model Eq. (2), together with the choice of rainbow-ladder truncation Γν(p,k)=γν, forms the kernel of the quark SDE, Eq. (1), which brings the dressed-quark propagator into a very simple form [56]:

    S1f(p)=iγp+Mf.

    (3)

    This is possible because the wave function renormalization trivially tends to unity in this case, and the quark mass function Mf becomes momentum independent:

    Mf=mf+4αeff3Λd4k(2π)4Tr[Sf(k)].

    (4)

    In this truncation, the quark-anitquark condensate is given by

    ˉqqf=NcΛd4k(2π)4Tr[Sf(k)],

    (5)

    with Nc=3 representing the number of colors. The form of the proposed gap equation, Eq. (4), is significantly similar to the NJL model gap equation, except the coupling parameter αeff=92G [49].

    Further simplification of Eq. (4) yields the following gap equation.

    Mf=mf+16αeff3Λd4k(2π)4Mfk2+M2f,

    (6)

    where Mf represents the dynamical mass, and the symbol Λ stresses the need to regularize the integrals. Using d4k=(1/2)k2dk2sin2θdθsinϕdϕdψ, and performing trivial regular integrations with the variable s=k2, the above expression reduces to:

    Mf=mf+αeffMf8π20dsss+M2f.

    (7)

    The integral in Eq. (7) is not convergent; therefore, we need to regularize it via the proper-time regularization scheme [47]. In this scheme, we take the exponent of the integrand's denominator and then introduce an additional infrared cutoff, in addition to the conventional ultraviolet cut-off, which is widely adopted in NJL model studies. Accordingly, the confinement is implemented using an infrared cut-off [57]. By adopting this scheme, the quadratic and logarithmic divergences are eliminated, and the axial-vector Ward-Takahashi identity [58, 59] is satisfied. From Eq. (7), the denominator of the integrand is given by

    1s+M2f=0dτeτ(s+M2f)τ2irτ2uvdτeτ(s+M2f)=eτ2uv(s+M2f)eτ2ir(s+M2f)s+M2f.

    (8)

    Here, τ1uv=Λuv is an ultra-violet regulator that plays the dynamical role and sets the scale for all dimensional quantities. τ1ir=Λir represents the infrared regulator whose non zero value implements confinement by ensuring the absence of quarks production thresholds [60]. Hence, τir corresponds to the confinement scale [22]. From Eq. (8), it is clear that the location of the original pole is at s=M2, which is canceled by the numerator. Accordingly, we discarded the singularities, and the propagator is free from real and complex poles, which is consistent with the definition of confinement: "an excitation described by a pole-less propagator would never reach its mass-shell" [57].

    After performing integration over "s", the gap equation is given by

    Mf=mf+M3fαeff8π2Γ(1,τuvM2f,τirM2f),

    (9)

    where

    Γ(a,x1,x2)=Γ(a,x1)Γ(a,x2),

    (10)

    and Γ(a,x)=xtα1etdt, which is the incomplete Gamma function. By using the parameters of Ref. [50], i.e., τir=(0.24GeV)1 and τuv=(0.905GeV)1, with the bare quark mu=md=0.007 GeV, we obtain the dynamical mass Mu=Md=0.367 GeV and the condensate ˉqqu=ˉqqd=(0.0143) GeV3.

    In this section, we study the gap equation in the presence of a uniform and homogeneous electromagnetic field with eEeB and at zero temperature. In the QCD Lagrangian, the interaction with the parallel electromagnetic field Aextμ embedded in the covariant derivative is expressed as

    Dμ=μiQfAextμ,

    (11)

    where Qf=(qu=+2/3,qd=1/3)e refers to the electric charges of u and d-quarks, respectively. We adopt the symmetric gauge vector potential Aextμ=(iEz,0,Bx,0) in Euclidean space, where both the electric and magnetic fields are selected along the z-axis. The gap equation in the presence of the parallel electromagnetic field remains in the form of Eq. (4), where Sf(k) is dressed with parallel background fields, i.e., Sf(k)~Sf(k). ~Sf(k), in the Schwinger proper time representation [1, 2, 5, 47], in the presence of the parallel magnetic field in Euclidean space [26, 29], is given by

    ~Sf(k)=0dτeτ(M2f+(k24+k23)tan(|QfE|τ)|QfE|τ+(k21+k22)tanh(|QfB|τ)|qfB|τ)×[γ4k4+Mf+tan(|QfE|τ)(γ4k3γ3k4)itanh(|QfBτ|)(γ1k2γ2k1)]×[1itanh(|QfB|τ)tan(|QfE|τ)γ5itanh(|QfB|τ)γ1γ2+tan(|QfE|τ)γ4γ3],

    (12)

    where the magnetic field couples with the coordinate indices 1 and 2, while the electric field couples with 3 and 4. Now, taking the trace of Eq. (12) and introducing both infrared and ultraviolet cut-offs, the gap equation in the presence of parallel electric and magnetic fields is given by

    ~Mf=mf+16αeff3f=u,d˜τ2irτ2uvdτ˜Mfeτ˜M2f×dk1(2π)dk2(2π)dk3(2π)dk4(2π)eτ((k24+k23)tan(|QfE|τ)|QfE|τ+(k21+k22)tanh(|QfB|τ)|QfB|τ.

    (13)

    The infrared cu-off τir is introduced in the model to mimic confinement by ensuring the absence of quarks production thresholds. In the presence of both eE and eB, it is required to vary slightly with both eE and eB. Therefore, the entanglement between dynamical chiral symmetry breaking and confinement is expressed using an explicit eE and eB-dependent regulator in the infrared [22]:

    ˜τir=τirMf˜Mf,

    (14)

    where Mf is the dynamical mass in the absence of background fields, and ˜Mf represents the eE and eB dependent dynamical mass. For chiral quarks (i.e., mf=0), the confining scale ˜τir diverges at the chiral symmetry restoration parameter near its critical value, which ensures the simultaneity of deconfinement and chiral symmetry restoration [22]. After the integration over k, the gap equation, Eq. (13), can be written as

    ˜Mf=mf+αeff3π2f=u,dˉτ2irτ2uvdτ˜Mfeτ˜M2f[|QfE|tan(|QfE|τ)|QfB|tanh(|QfB|τ)].

    (15)

    When both fields eE and eB0, Eq. (15) reduces to Eq. (9). The gap equation for the pure electric field can be obtained by setting eB0, while for pure magnetic field, eE0. The quark-antiquark condensate in the presence of background fields is in this form.

    In this study, we adopt two flavors f=2, i.e., u- and d-quarks. We use the same current quark mass for the up and down quarks, i.e., mu=md=0.007 GeV, such that the iso-spin symmetry is preserved. As is well-known, the response to the electromagnetic field is different for u and d-quarks because they have different electric charges [1, 5].

    The numerical solution of the gap equation Eq. (15) as a function of eB for fixed values of eE, is plotted inFig. 1. In the pure magnetic case (eE0), the dynamical mass ~Mf monotonically increases with an increase in eB , which ensures the magnetic catalysis phenomenon. The increase in ~Mf as a function of eB reduces in magnitude upon varying eE from its smaller to larger given values. We noticed that at eE0.33 GeV2, the dynamical mass ~Mf exhibits the de Haas-van Alphen oscillatory type behavior [61]. In other words, it remains constant for small eB values, then monotonically decreases with an increase in eB, and abruptly declines to lower values in the region eB[0.30.54] GeV2, where the chiral symmetry is partially restored via first-order phase transition; above it, the value increases again. This phenomenon can be attributed to the strong competition that occurs between parallel eB and eE, i.e., on the one hand, eB enhances the mass function, while on the other hand, eE suppresses it. We also sketch a zoom-in plot of the mass function, which shows the de Haas-van Alphen oscillatory behavior in the region eB[0.30.54] GeV2 with eE=0.33 GeV2 , as illustrated in Fig. 2. Such a behavior type is also explored and discussed in the other effective model of QCD, for example, refer to Refs. [29, 62].

    Figure 1

    Figure 1.  (color online) Behavior of the dynamically generated quark mass, Eq. (15), as a function of the magnetic field strength eB at several given values of eE.

    Figure 2

    Figure 2.  Zoom-in plot of the mass function that shows the de Haas-van Alphen oscillatory behavior in the magnetic field region eB=[0.30.5] GeV2 with the electric field eE=0.33 GeV2

    The behavior of the quark-antiquark condensate, Eq. (16), as a function of eB at various given values of eE, is illustrated in the Fig. 3. For a pure magnetic case (eE0), the magnetic field strength eB facilitates the formation of the quark-antiquark condensate. For given non zero values of eE, particularly, at eE0.33 GeV2 the evolution in the condensate is suppressed in the region eB=[0.30.54] GeV2. The nature of the transition in this region suddenly changes to first-order; however, above this region, it is enhanced with eB. The pseudo-critical values of the fields at which the chiral symmetry partially is restored and first-order phase transition occurres are eBc0.3 GeV2 and eEc=0.33 GeV2. Although such values of electric or magnetic field strength are sufficiently large in terms of what is typically generated during heavy-ion collisions, they may be relevant to astronomical objects, such as neutron stars, magnetars, etc. The confinement parameter ˜τ1ir as a function of eB for different fixed values of eE is presented in Fig. 4. We observed the same pseudo-critical fields eBc0.3 GeV2 and eEc=0.33 GeV2 for the confinement transition, similar to the case of chiral symmetry breaking.

    Figure 3

    Figure 3.  (color online) Quark-antiquark condensate, Eq. (16), as a function of eB for several values of eE.

    Figure 4

    Figure 4.  (color online) Confinement scale ˜τ1ir as a function of eB for several values of eE.

    In the following, we discuss the variation of ~Mf and ˜τ1ir as a function of eE, in the pure electric field, as well as for several non-zero values of eB. In Figs. 5, 6, and 7, we illustrate the behaviors of all three parameters as a function of eE, for various values of eB0.

    Figure 5

    Figure 5.  (color online) Dynamical mass as a function of eE for several values of eB.

    Figure 6

    Figure 6.  (color online) Quark-antiquark condensate, Eq. (16), as a function of eE for several values of eB.

    Figure 7

    Figure 7.  (color online) Confinement scale ˜τ1ir as a function of eE for a given several values of eB.

    In the pure electric field case (eB0), the three-parameters monotonically decrease with an increase in eE, and at a pseudo-critical field strength eEχ,Cc , the chiral symmetry is partially restored and deconfinement transitions occurred; hence, the nature of the transitions here becomes smooth cross-over. Accordingly, we observed the chiral rotation or chiral electric inhibition effect in the contact interaction model, as already predicted by other effective models of QCD [1, 2, 25-29]. For non-zero eB values, we determined an interesting behavior for all the three parameters, as a function of eE. All three parameters are enhanced for given smaller to larger values of eB, except the chiral symmetry restoration and deconfinement regions, where all the parameters are suppressed by higher eB values. The pseudo-critical field strength Eχ,Cc decreases with an increase in eB, and at a pseudo-critical eBc0.3 GeV2, the transition changes from cross-over to first order. The nontrivial behaviors of all three parameters represent the competition between the magnetic catalysis and electric inhibition effects, both induced by eB in the presence of parallel eE [29]. As elucidated and validated in Ref. [25], the cause of the chiral inhibition effect is the second Lorentz invariant of the electromagnetic field, EB. The magnitude of eEχ,Cc at which the chiral symmetry restoration and deconfinement occurred are triggered from the inflection point of the electric gradients of andeE˜τ1ir, as presented in Figs. 8 and 9, respectively. In the pure electric case, the chiral symmetry restoration and deconfinement transition occur at eEχ,Cc0.34 GeV2. For several given values of eB0, eEχ,Cc decreases with an increase in eB. We determined a smooth cross-over phase transition up to the pseudo-critical magnetic field strength eBc0.3 GeV2, where the transitions take the first-order nature. Here, eEχ,Cc remains constant in the region eB[0.30.54] GeV2 and then increases with larger values of eB , as shown in Fig. 10. A similar behavior has already been demonstrated in [29]. The boundary point where the cross-over phase transition ends and the first order phase transition starts is known as a critical end point, and its co-ordinates are at (eBp0.3,eEp0.33) GeV2.

    Figure 8

    Figure 8.  (color online) Electric gradient of the quark-antiquark condensate as a function of eE for fixed values of eB. At a particular fixed value of eB=0.3 GeV2, the electric gradient diverges at eEχc,0.33 GeV2 and above this value, the first order phase transition occurs.

    Figure 9

    Figure 9.  (color online) Electric gradient of the confinement scale eE˜τ1ir, as a function of eE for different fixed values of eB. At a particular fixed value of eB=0.3 GeV2, the electric gradient of the confinement scale diverges at eECc0.33 GeV2.

    Figure 10

    Figure 10.  (color online) Phase diagram of chiral symmetry and confinement for eEC,χc vs eB. eE=eEC,χc is obtained from the inflection points of the electric gradient of the condensate and the confining length scale eE˜τ1ir.

    In this section, we elucidate the behaviors of the dynamical mass, condensate, and confinement length scale at a finite temperature and in the presence of parallel electric and magnetic fields. We also explore the IEC, MC, and IMC phenomena, as well as the competition among them. Finally, we sketch the QCD phase diagram.

    The finite temperature version of the gap equation, Eq. (13), in the presence of parallel electric and magnetic fields can be obtained by adopting the standard convention for momentum integration:

    d4k(2π)4Tnd3k(2π)3,

    (17)

    and the four momenta k(ωn,k), with ωn=(2n+1)πT, represent the fermionic Matsubara frequencies. The Lorentz structure does not preserve anymore at finite temperature. By making the following replacements in Eq. (13),

    dk1(2π)dk2(2π)dk3(2π)dk4(2π)Tndk1(2π)dk2(2π)dk3(2π),

    and k4ωn, we have

    ^Mf=mf+16αeff3Tnf=u,dˆτ2irτ2uvdτ^MfeτˆM2f×dk1(2π)dk2(2π)dk3(2π)eτ((ω2n+k23)tan(|QfE|τ)|QfE|τ+(k21+k22)tanh(|QfB|τ)|QfB|τ),

    (18)

    with ˆMf=Mf(eE,eB,T). Performing sum over Matsubara frequencies and integrating over ks, the gap equation can be written as

    ˆMf=mf+αeff3π2f=u,dˆτ2irτ2uvdτˆMfeτˆM2fΘ3(π2,e|QfE|4T2tan(|QfE|τ))×[|QfE|tan(|QfE|τ)|QfB|tanh(|QfB|τ)],

    (19)

    where Θ3(π2,ex) is the third Jacobi's theta function. The confinement scale ˆτir here slightly varies with T, eE , and eB, and it takes the form

    ˆτir=τirMfˆMf.

    (20)

    Here, Mf=Mf(0,0,0) is the dynamical mass at eB=eE=T=0, whereas ˆMf=Mf(eE,eB,T) is the variation of the dynamical mass at finite eB, eE, and T. The quark-antiquark condensate is given by the following.

    First, we consider the case of a pure electric field at a finite temperature. The numerical solution of Eq. (19) as a function of T, for different given values of eE, is presented in Fig. 11. The dynamical mass ˆMf monotonically decreases with the increase of T until the dynamical chiral symmetry is partially restored. The response of eE is to suppress the dynamical mass. The quark-antiquark condensate, Eq. (14), and confinement length scale, Eq. (15), as a function of T for various given values of eE are depicted in Figs. 12 and 13, respectively. We infer that both parameters decrease with temperature T, and at a pseudo-critical temperature Tχ,Cc, the chiral symmetry is partially restored, and deconfinement occurs. We note that the electric field eE suppresses both parameters in the low-temperature regions and also reduces the pseudo-critical temperature Tχ,Cc. In Figs. 14 and 15, we plotted the thermal gradients of the condensate and the confinement length scale, respectively. The peaks in the thermal gradients of both parameters shift toward the low-temperature regions upon increasing the value of eE. Consequently, the critical temperature Tχ,Cc decreases with an increase in electric field strength eE; hence, inverse electric catalysis is observed at a finite temperature in the proposed contact interaction model. Therefore, the observations of this study are consistent with other effective models of QCD [28-30]. The magnitude of the critical temperature Tχ,Cc0.22 GeV2 is obtained from the inflection points of the thermal gradients of and Tˆτ1ir.

    Figure 11

    Figure 11.  (color online) Dynamical mass, Eq. (19), as a function of temperature for various given values of electric field strength eE.

    Figure 12

    Figure 12.  (color online) Quark-antiquark condensate, Eq. (21), plotted as a function of temperature for various given values of electric field eE.

    Figure 13

    Figure 13.  (color online) Behavior of confinement scale, Eq. (20), plotted as a function of temperature for different given values eE.

    Figure 14

    Figure 14.  (color online) Thermal gradient of the quark-antiquark condensate plotted as a function of temperature T for various values of the electric field eE. The peaks in the derivatives shift toward low temperature regions from small to large values of eE.

    Figure 15

    Figure 15.  (color online) Thermal gradient of the confinement scale Tˆτ1ir plotted as a function of temperature T for various values of eE.

    Second, we consider the case of a pure magnetic field at finite temperature T. We plot the thermal gradient of the quark-antiquark condensate and the confinement length scale, as presented in Figs. 16 and 17, respectively. We note that the inflection points in both parameters shift toward the higher temperatures. In other words, eB enhances the pseudo-critical temperature Tχ,Cc of chiral symmetry restoration and deconfinement. Hence, in a pure magnetic case and at a finite temperature, we observe the magnetic catalysis phenomenon, which has already been observed in the CI model [22, 63].

    Figure 16

    Figure 16.  (color online) Thermal gradient of the quark-antiquark condensate plotted as a function of temperature T for various values of magnetic field eB. From the plots, it is evident that the peaks shift toward higher temperatures.

    Figure 17

    Figure 17.  (color online) Thermal gradient of the confinement scale Tˆτ1ir plotted as a function of temperature for several values of the magnetic field strength field eB.

    In most effective model calculations of QCD, it is well demonstrated that to reproduce the inverse magnetic catalysis effect as predicted by the lattice QCD [13, 14], the effective coupling must be taken as magnetic field dependent [17, 22, 28, 63] or both temperature and magnetic field dependent [16, 19, 23, 63, 64]. In the case of this study, we just employed the following functional form of the eB-dependent effective coupling αeff(eB) [22], where the coupling decreases with the magnetic field strength as

    αeff(κ)=αeff(1+aκ2+bκ31+cκ2+dκ4),

    (22)

    here κ=eB/Λ2QCD, with ΛQCD=0.24 GeV. The parameters a, b, c, and d were extracted to reproduce the behavior of critical temperature Tχ,Cc for the chiral symmetry restoration and deconfinement in the presence of magnetic field strength, obtained by lattice QCD simulations [13, 14]. The thermal gradients of the condensate and the confinement length scale, with the magnetic field dependent coupling, Eq. (22), is plotted in Figs. 18 and 19, respectively. It can be observed that the critical Tχ,Cc decreases with an increase in eB ; hence, the inverse magnetic catalysis can be observed at a finite T.

    Figure 18

    Figure 18.  (color online) Thermal gradient of the condensate with magnetic field dependent coupling, Eq. (22), plotted as a function of temperature for several values of magnetic field strength. The IMC effect is depicted in this figure.

    Figure 19

    Figure 19.  (color online) Behavior of Tˆτ1ir with the magnetic field dependent coupling, Eq. (22), as a function of temperature for several values of magnetic field strength eB.

    In Fig. 20, we sketch the combined phase diagram in the Tχ,CceE,eB planes, where we demonstrated the inverse electric catalysis, magnetic catalysis, and inverse magnetic catalysis (with magnetic dependent coupling). In the pure electric background, the solid-black curve indicates that the critical temperature Tχ,Cc decreases with an increase in eE; hence, the electric field strength inhibits the chiral symmetry breaking and confinement. In the pure magnetic limit (without magnetic field dependent coupling), the magnetic field eB enhances the critical temperature Tχ,Cc, and thus eB acts as a facilitator of chiral symmetry breaking and confinement. If we use the magnetic field dependent coupling, it can be observed that the critical temperature Tχ,Cc decreases as the magnetic field strength eB increases (red dotted-dashed curve). In this case, eB acts as an inhibitor of the chiral symmetry and confinement.

    Figure 20

    Figure 20.  (color online) Combined phase diagram in the Tχ,Cc vs eE,eB planes of chiral symmetry breaking and confinement. The plot shows the evalution of IEC as well as MC without and IMC with eB-dependent coupling, Eq. (22). The Tχ,Cc values are obtained from the inflection points in Tˆσ and Tˆτ1ir.

    Third, we adopt both the non-zero values of eE and eB , and draw the phase diagram in the Tχ,CceE plane for various given values of eB , as shown in Fig. 21. At this point, the competition between IEC vs MC starts: the eB tends to catalyze the chiral symmetry breaking and confinement; consequently, Tχ,Cc is enhanced. In contrast, eE tries to inhibit the chiral phase transition and tends to suppres Tχ,Cc; finally, it is concluded with a combined inverse electromagnetic catalysis (IEMC). This may be different for a very strong eB, where the eB dominates over the eE. Next, we consider the eB-dependent coupling, Eq. (22), and sketch the phase diagram Tχ,Cc vs eE for various given values of eB in Fig. 22. It can be observed that the effect of parallel eE and eB with eB-dependent coupling tends to reduce Tχ,Cc. Therefore, both eE and eB produce IEC and IMC simultaneously.

    Figure 21

    Figure 21.  (color online) Phase diagram in the Tχ,CceE plane of chiral symmetry breaking and confinement for various given values of eB. Here, eE inhibits the chiral phase transition, whereas eB facilitates it; consequently, the IEMC effect is observed.

    Figure 22

    Figure 22.  (color online) Phase diagram in the Tχ,Cc vs eE of chiral symmetry breaking and confinement for various given values of eE with eB-dependent coupling Eq. (22). Here, both fields eE and eB inhibit the chiral phase transition; hence, the IEMC is observed.

    In this work, we studied the influence of uniform, homogeneous, and external parallel electric and magnetic fields on the chiral phase transitions. In this context, we implemented the Schwinger-Dyson formulation of QCD, with a gap equation kernel comprising a symmetry preserving vector-vector contact interaction model of quarks in a rainbow-ladder truncation. Subsequently, we adopted the well-known Schwinger proper-time regularization procedure. The findings of this study are presented as follows.

    At zero temperature, in the pure magnetic background, the magnetic field facilitated the dynamical chiral symmetry breaking and confinement; hence, we observed the magnetic catalysis phenomenon. However, the electric field tends to restore chiral symmetry and deconfinement above the pseudo-critical electric field eEχ,Cc0.34 GeV2, i.e., the chiral rotation effect demonstrated in the proposed model. We observed that the electric field acted as an inhibitor of the chiral symmetry breaking and confinement. When both eE and eB were considered, we determined the magnetic catalysis effect for small given values of eE up to and above the critical electric field strength eEc0.33 GeV2, where the other parameters exhibited the de Haas-van Alphen oscillatory type behavior in the region eB=[0.30.54] GeV2. In this particular region,eE dominated over eB , and we observed the electric chiral inhibition effect, while above that region where eB was superior to eE, we observed magnetic catalysis again. We also realized that the pseudo-critical strength eEχ,Cc was suppressed with the increase in eB , and the transition nature was determined to be the smooth cross-over nature up to and above the pseudo-critical magnetic field strength eBχ,Cc0.3 GeV2 , where the transition suddenly changed to the first order transition. Subsequently, we located the position of the critical endpoint at (eBp0.3,eEp0.33) GeV2. We further determined that eEχ,Cc initially remains constant in the limited region eB=[0.30.54] GeV2 and then increases with larger values of eB.

    Finally, we sketched the phase diagram at a finite temperature and in the presence of parallel electric and magnetic fields. At a finite T, in the pure electric field limit, we determined that the pseudo-critical temperature decreased as we increasedeE; hence, the inverse electric catalysis was observed. However, for the pure magnetic field background, we observed the magnetic catalysis effect in the mean-field approximation and inverse magnetic catalysis with eB-dependent coupling. The combined effect of both eE and eB on Tχ,Cc yielded the inverse electromagnetic catalysis, with and without eBdependent effective coupling of the model.

    Qualitatively and quantitatively, the predictions of the presented CI-model agree well with results obtained from other effective QCD models, as well as modern Lattice QCD results. In the near future, we plan to extend this work to study the Schwinger pair production rate, including the dynamical chiral symmetry breaking for a higher number of colors, flavors, and in a parallel electromagnetic field. We are also interested in extending this work to study the IEMC phenomenon with the electric field and temperature-dependent coupling. Regarding the dependence of the model coupling on the considered electric field, there is no first-principle simulation for comparisons because lattice simulations in the presence of real electric background fields are not feasible to date. However, progress in this case has commenced, and results will be published somewhere else. The next strategy will be to study the properties of light hadrons in the background of parallel electric and magnetic fields.

    The author acknowledges A. Bashir, A. Raya, R.L.S. Farias, and E. V. Gorbar for their valuable suggestions in the process of completing this work. The author also thanks his colleagues at the Institute of Physics, Gomal University, Pakistan.

    [1] S. P. Klevansky and R. H. Lemmer, Phys. Rev. D 38, 3559-3565 (1988)
    [2] H. Suganuma and T. Tatsumi, Annals. Phys. 208, 470-508 (1991) doi: 10.1016/0003-4916(91)90304-Q
    [3] K. G. Klimenko , Theor. Math. Phys. 89, 1161-1168 (1992)
    [4] K. G. Klimenko, Teor. Mat. Fiz. 89, 211 (1991)
    [5] S. P. Klevansky, Rev. Mod. Phys 64, 649-708 (1992) doi: 10.1103/RevModPhys.64.649
    [6] I. V. Krive and S. A. Naftulin, Phys. Rev. D 46, 2737-2740 (1992)
    [7] V. P. Gusynin, V. A. Miransky, and I. A. Shovkovy, Phys. Rev. Lett. 73, 3499-3502 (1994) doi: 10.1103/PhysRevLett.73.3499
    [8] V. P. Gusynin, V. A. Miransky, and I. A. Shovkovy, Erratum: Phys. Rev. Lett. 76, 1005 (1996), arXiv:hep-ph/9405262
    [9] V. P. Gusynin, V.A. Miransky, and I. A. Shovkovy, Phys. Rev. D 52, 4718-4735 (1995), arXiv:hep-th/9407168
    [10] V. P. Gusynin, V. A. Miransky, and I. A.Shovkovy, Phys. Lett. B 349, 477-483 (1995), arXiv:hep-ph/9412257
    [11] K. G. Klimenko, Theor. Math. Phys. 90, 1-6 (1992) doi: 10.1007/BF01018812
    [12] K. G. Klimenko, Teor. Mat. Fiz. 90, 3 (1992)
    [13] G. S. Bali, F. Bruckmann, G. Endrodi et al., PoS LATTICE2011, 192 (2011), arXiv:1111.5155
    [14] G. S. Bali, F. Bruckmann, G. Endrodi et al., JHEP 04, 130 (2013), arXiv:1303.1328
    [15] G. Endrodi, M. Giordano, S. D. Katz et al., JHEP 07, 007 (2019), arXiv:1904.10296
    [16] R.L S. Farias, K. P. Gomes, G. I. Krein et al., Phys. Rev. C 90, 025203 (2014), arXiv:1404.3931
    [17] P. Costa, M. Ferreira, D. P. Menezes et al., Phys. Rev. D 92, 036012 (2015), arXiv:1508.07870
    [18] E. J. Ferrer, V. de la Incera, and X. J. Wen, Phys. Rev. D 91, 054006 (2015), arXiv:1407.3503
    [19] A. Ayala, C. A. Dominguez, L. A. Hernandez et al., Phys. Lett. B 759, 99-103 (2016), arXiv:1510.09134
    [20] J. O. Andersen, W. R. Naylor, and A. Tranberg, JHEP 02, 042 (2015), arXiv:1410.5247
    [21] N. Mueller and J. M. Pawlowski, Phys. Rev. D 91, 116010 (2015), arXiv:1502.08011
    [22] A. Ahmad and A. Raya, J. Phys. G 43, 065002 (2016), arXiv:1602.06448
    [23] R. L. S. Farias, V. S. Timoteo, S. S. Avancini et al., Eur. Phys. J. A 53, 101 (2017), arXiv:1603.03847
    [24] S. He, Y. Yang, and P. H. Yuan, Analytic Study of Magnetic Catalysis in Holographic QCD, (2020), arXiv:2004.01965
    [25] A. Yu Babansky, E. V. Gorbar, and G. V. Shchepanyuk, Phys. Lett. B 419, 272-278 (1998), arXiv:hep-th/9705218
    [26] G. Cao and X. G. Huang, Phys. Rev. D 93, 016007 (2016), arXiv:1510.05125
    [27] W. R. Tavares and S. S. Avancini, Phys. Rev. D 97, 094001 (2018), arXiv:1801.10566
    [28] M. Ruggieri and G. X. Peng, Phys. Rev. D 93, 094021 (2016), arXiv:1602.08994
    [29] L. Wang and G. Cao, Phys. Rev. D 97, 034014 (2018), arXiv:1712.09780
    [30] W. R. Tavares, R. L. S. Farias, and S. S. Avancini, Phys. Rev. D 101, 016017 (2020), arXiv:1912.00305
    [31] A. Bzdak and V. Skokov, Phys. Lett. B 710, 171-174 (2012), arXiv:1111.1949
    [32] W. T. Deng and X. G. Huang, Phys. Rev. C 85, 044907 (2012), arXiv:1201.5108
    [33] J. Bloczynski, X. G. Huang, X. Zhang et al., Phys. Lett. B 718, 1529-1535 (2013), arXiv:1209.6594
    [34] J. Bloczynski, X. G. Huang, X. Zhang et al., Nucl. Phys. A 939, 85-100 (2015), arXiv:1311.5451
    [35] D. E. Kharzeev, L. D. McLerran, and H. J. Warringa, Nucl. Phys. A 803, 227-253 (2008), arXiv:0711.0950
    [36] K. Fukushima, D. E. Kharzeev, and H. J. Warringa, Phys. Rev. D 78, 074033 (2008), arXiv:0808.3382
    [37] X. G. Huang, and J. Liao, Phys. Rev. Lett. 110, 232302 (2013), arXiv:1303.7192 doi: 10.1103/PhysRevLett.110.232302
    [38] Y. Jiang, X. G. Huang, and J. Liao, Phys. Rev. D 91, 045001 (2015), arXiv:1409.6395
    [39] I. Karpenko and F. Becattini, Eur. Phys. J. C 77, 213 (2017), arXiv:1610.04717
    [40] X. L. Xia, H. Li, Z. B. Tang et al., Phys. Rev. C 98, 024905 (2018), arXiv:1803.00867
    [41] D. X. Wei, W. T. Deng, and X. G. Huang, Phys. Rev. C 99, 014905 (2019), arXiv:1810.00151
    [42] L. Wang, G. Cao, X. G. Huang et al., Phys. Lett. B 780, 273-282 (2018), arXiv:1801.01682
    [43] G. Cao, Phys. Rev. D 101, 094027 (2020), arXiv:2003.09209
    [44] Q. Li, D. E. Kharzeev, C. Zhang et al., Nature Phys. 12, 550-554 (2016), arXiv:1412.6543 doi: 10.1038/nphys3648
    [45] X. L. Zhao, G. L. Ma, and Y. G. Ma, Phys. Lett. B 792, 413-418 (2019), arXiv:1901.04156
    [46] H. L. L. Roberts, A. Bashir, L. X. Gutierrez-Guerrero et al., Phys. Rev. C 83, 065206 (2011), arXiv:1102.4376
    [47] J. S. Schwinger, Phys. Rev. 82, 664-679 (1951) doi: 10.1103/PhysRev.82.664
    [48] K. L. Wang, Y. X. Liu, L. Chang et al., Phys. Rev. D 87, 074038 (2013), arXiv:1301.6762
    [49] F. Marquez, A. Ahmad, M. Buballa et al., Phys. Lett. B 747, 529-535 (2015), arXiv:1504.06730
    [50] L. X. Gutierrez-Guerrero, V. Bashir, I. C. Cloet et al., Phys. Rev. C 81, 065202 (2010), arXiv:1002.1968
    [51] H. L. L. Roberts, L. Chang, I. C. Cloet et al., Few Body Syst. 51, 1-25 (2011), arXiv:1101.4244 doi: 10.1007/s00601-011-0225-x
    [52] H. L. L. Roberts, C. D. Roberts, A. Bashir et al., Phys. Rev. C 82, 065202 (2010), arXiv:1009.0067
    [53] C. Chen, L. Chang, C. D. Roberts et al., Few Body Syst. 53, 293-326 (2012), arXiv:1204.2553 doi: 10.1007/s00601-012-0466-3
    [54] P. Boucaud, J. P. Leroy, A. L. Yaouanc et al., Few Body Syst. 53, 387-436 (2012), arXiv:1109.1936 doi: 10.1007/s00601-011-0301-2
    [55] A. Ahmad, A. Martínez and A. Raya, Phys. Rev. D 98, 054027 (2018), arXiv:1809.05545
    [56] A. Ahmad, A. Ayala, A. Bashir et al., J. Phys. Conf. Ser. 651, 012018 (2015) doi: 10.1088/1742-6596/651/1/012018
    [57] D. Ebert, T. Feldmann, and H. Reinhardt, Phys. Lett. B 388, 154-160 (1996), arXiv:hep-ph/9608223
    [58] J. C. Ward, Phys. Rev. 78, 182 (1950)
    [59] Y. Takahashi, Nuovo Cim. 6, 371 (1957) doi: 10.1007/BF02832514
    [60] C. D. Roberts, M. S. Bhagwat, A. Holl et al., Eur. Phys. J. ST 140, 53-116 (2007), arXiv:0802.0217 doi: 10.1140/epjst/e2007-00003-5
    [61] W. de Haas and P. van Alphen, Leiden Comm. 208d 212a, (1930)
    [62] G. Cao and P. Zhuang, Phys. Rev. D 92, 105030 (2015), arXiv:1505.05307 doi: 10.1103/PhysRevD.92.105030
    [63] A. Ahmad, A. Bashir, M. A. Bedolla et al., J. Phys. G 48(7), 075002 (2021), arXiv:2008.03847
    [64] A. Bandyopadhyay and R. L. S. Farias, Inverse magnetic catalysis -- how much do we know about?, (2020), arXiv:2003.11054
  • [1] S. P. Klevansky and R. H. Lemmer, Phys. Rev. D 38, 3559-3565 (1988)
    [2] H. Suganuma and T. Tatsumi, Annals. Phys. 208, 470-508 (1991) doi: 10.1016/0003-4916(91)90304-Q
    [3] K. G. Klimenko , Theor. Math. Phys. 89, 1161-1168 (1992)
    [4] K. G. Klimenko, Teor. Mat. Fiz. 89, 211 (1991)
    [5] S. P. Klevansky, Rev. Mod. Phys 64, 649-708 (1992) doi: 10.1103/RevModPhys.64.649
    [6] I. V. Krive and S. A. Naftulin, Phys. Rev. D 46, 2737-2740 (1992)
    [7] V. P. Gusynin, V. A. Miransky, and I. A. Shovkovy, Phys. Rev. Lett. 73, 3499-3502 (1994) doi: 10.1103/PhysRevLett.73.3499
    [8] V. P. Gusynin, V. A. Miransky, and I. A. Shovkovy, Erratum: Phys. Rev. Lett. 76, 1005 (1996), arXiv:hep-ph/9405262
    [9] V. P. Gusynin, V.A. Miransky, and I. A. Shovkovy, Phys. Rev. D 52, 4718-4735 (1995), arXiv:hep-th/9407168
    [10] V. P. Gusynin, V. A. Miransky, and I. A.Shovkovy, Phys. Lett. B 349, 477-483 (1995), arXiv:hep-ph/9412257
    [11] K. G. Klimenko, Theor. Math. Phys. 90, 1-6 (1992) doi: 10.1007/BF01018812
    [12] K. G. Klimenko, Teor. Mat. Fiz. 90, 3 (1992)
    [13] G. S. Bali, F. Bruckmann, G. Endrodi et al., PoS LATTICE2011, 192 (2011), arXiv:1111.5155
    [14] G. S. Bali, F. Bruckmann, G. Endrodi et al., JHEP 04, 130 (2013), arXiv:1303.1328
    [15] G. Endrodi, M. Giordano, S. D. Katz et al., JHEP 07, 007 (2019), arXiv:1904.10296
    [16] R.L S. Farias, K. P. Gomes, G. I. Krein et al., Phys. Rev. C 90, 025203 (2014), arXiv:1404.3931
    [17] P. Costa, M. Ferreira, D. P. Menezes et al., Phys. Rev. D 92, 036012 (2015), arXiv:1508.07870
    [18] E. J. Ferrer, V. de la Incera, and X. J. Wen, Phys. Rev. D 91, 054006 (2015), arXiv:1407.3503
    [19] A. Ayala, C. A. Dominguez, L. A. Hernandez et al., Phys. Lett. B 759, 99-103 (2016), arXiv:1510.09134
    [20] J. O. Andersen, W. R. Naylor, and A. Tranberg, JHEP 02, 042 (2015), arXiv:1410.5247
    [21] N. Mueller and J. M. Pawlowski, Phys. Rev. D 91, 116010 (2015), arXiv:1502.08011
    [22] A. Ahmad and A. Raya, J. Phys. G 43, 065002 (2016), arXiv:1602.06448
    [23] R. L. S. Farias, V. S. Timoteo, S. S. Avancini et al., Eur. Phys. J. A 53, 101 (2017), arXiv:1603.03847
    [24] S. He, Y. Yang, and P. H. Yuan, Analytic Study of Magnetic Catalysis in Holographic QCD, (2020), arXiv:2004.01965
    [25] A. Yu Babansky, E. V. Gorbar, and G. V. Shchepanyuk, Phys. Lett. B 419, 272-278 (1998), arXiv:hep-th/9705218
    [26] G. Cao and X. G. Huang, Phys. Rev. D 93, 016007 (2016), arXiv:1510.05125
    [27] W. R. Tavares and S. S. Avancini, Phys. Rev. D 97, 094001 (2018), arXiv:1801.10566
    [28] M. Ruggieri and G. X. Peng, Phys. Rev. D 93, 094021 (2016), arXiv:1602.08994
    [29] L. Wang and G. Cao, Phys. Rev. D 97, 034014 (2018), arXiv:1712.09780
    [30] W. R. Tavares, R. L. S. Farias, and S. S. Avancini, Phys. Rev. D 101, 016017 (2020), arXiv:1912.00305
    [31] A. Bzdak and V. Skokov, Phys. Lett. B 710, 171-174 (2012), arXiv:1111.1949
    [32] W. T. Deng and X. G. Huang, Phys. Rev. C 85, 044907 (2012), arXiv:1201.5108
    [33] J. Bloczynski, X. G. Huang, X. Zhang et al., Phys. Lett. B 718, 1529-1535 (2013), arXiv:1209.6594
    [34] J. Bloczynski, X. G. Huang, X. Zhang et al., Nucl. Phys. A 939, 85-100 (2015), arXiv:1311.5451
    [35] D. E. Kharzeev, L. D. McLerran, and H. J. Warringa, Nucl. Phys. A 803, 227-253 (2008), arXiv:0711.0950
    [36] K. Fukushima, D. E. Kharzeev, and H. J. Warringa, Phys. Rev. D 78, 074033 (2008), arXiv:0808.3382
    [37] X. G. Huang, and J. Liao, Phys. Rev. Lett. 110, 232302 (2013), arXiv:1303.7192 doi: 10.1103/PhysRevLett.110.232302
    [38] Y. Jiang, X. G. Huang, and J. Liao, Phys. Rev. D 91, 045001 (2015), arXiv:1409.6395
    [39] I. Karpenko and F. Becattini, Eur. Phys. J. C 77, 213 (2017), arXiv:1610.04717
    [40] X. L. Xia, H. Li, Z. B. Tang et al., Phys. Rev. C 98, 024905 (2018), arXiv:1803.00867
    [41] D. X. Wei, W. T. Deng, and X. G. Huang, Phys. Rev. C 99, 014905 (2019), arXiv:1810.00151
    [42] L. Wang, G. Cao, X. G. Huang et al., Phys. Lett. B 780, 273-282 (2018), arXiv:1801.01682
    [43] G. Cao, Phys. Rev. D 101, 094027 (2020), arXiv:2003.09209
    [44] Q. Li, D. E. Kharzeev, C. Zhang et al., Nature Phys. 12, 550-554 (2016), arXiv:1412.6543 doi: 10.1038/nphys3648
    [45] X. L. Zhao, G. L. Ma, and Y. G. Ma, Phys. Lett. B 792, 413-418 (2019), arXiv:1901.04156
    [46] H. L. L. Roberts, A. Bashir, L. X. Gutierrez-Guerrero et al., Phys. Rev. C 83, 065206 (2011), arXiv:1102.4376
    [47] J. S. Schwinger, Phys. Rev. 82, 664-679 (1951) doi: 10.1103/PhysRev.82.664
    [48] K. L. Wang, Y. X. Liu, L. Chang et al., Phys. Rev. D 87, 074038 (2013), arXiv:1301.6762
    [49] F. Marquez, A. Ahmad, M. Buballa et al., Phys. Lett. B 747, 529-535 (2015), arXiv:1504.06730
    [50] L. X. Gutierrez-Guerrero, V. Bashir, I. C. Cloet et al., Phys. Rev. C 81, 065202 (2010), arXiv:1002.1968
    [51] H. L. L. Roberts, L. Chang, I. C. Cloet et al., Few Body Syst. 51, 1-25 (2011), arXiv:1101.4244 doi: 10.1007/s00601-011-0225-x
    [52] H. L. L. Roberts, C. D. Roberts, A. Bashir et al., Phys. Rev. C 82, 065202 (2010), arXiv:1009.0067
    [53] C. Chen, L. Chang, C. D. Roberts et al., Few Body Syst. 53, 293-326 (2012), arXiv:1204.2553 doi: 10.1007/s00601-012-0466-3
    [54] P. Boucaud, J. P. Leroy, A. L. Yaouanc et al., Few Body Syst. 53, 387-436 (2012), arXiv:1109.1936 doi: 10.1007/s00601-011-0301-2
    [55] A. Ahmad, A. Martínez and A. Raya, Phys. Rev. D 98, 054027 (2018), arXiv:1809.05545
    [56] A. Ahmad, A. Ayala, A. Bashir et al., J. Phys. Conf. Ser. 651, 012018 (2015) doi: 10.1088/1742-6596/651/1/012018
    [57] D. Ebert, T. Feldmann, and H. Reinhardt, Phys. Lett. B 388, 154-160 (1996), arXiv:hep-ph/9608223
    [58] J. C. Ward, Phys. Rev. 78, 182 (1950)
    [59] Y. Takahashi, Nuovo Cim. 6, 371 (1957) doi: 10.1007/BF02832514
    [60] C. D. Roberts, M. S. Bhagwat, A. Holl et al., Eur. Phys. J. ST 140, 53-116 (2007), arXiv:0802.0217 doi: 10.1140/epjst/e2007-00003-5
    [61] W. de Haas and P. van Alphen, Leiden Comm. 208d 212a, (1930)
    [62] G. Cao and P. Zhuang, Phys. Rev. D 92, 105030 (2015), arXiv:1505.05307 doi: 10.1103/PhysRevD.92.105030
    [63] A. Ahmad, A. Bashir, M. A. Bedolla et al., J. Phys. G 48(7), 075002 (2021), arXiv:2008.03847
    [64] A. Bandyopadhyay and R. L. S. Farias, Inverse magnetic catalysis -- how much do we know about?, (2020), arXiv:2003.11054
  • 加载中

Cited by

1. Tavares, W.R., Ramos, R.O., Farias, R.L.S. et al. Charged scalars at finite electric field and temperature in the optimized perturbation theory[J]. Physical Review D, 2024, 110(11): 116024. doi: 10.1103/PhysRevD.110.116024
2. Upadhaya, S.. Fluctuations and Correlations of Conserved Charges Serving as Signals for QGP Production: An Overview from Polyakov Loop Enhanced Nambu–Jona-Lasinio Model[J]. Universe, 2024, 10(8): 332. doi: 10.3390/universe10080332
3. Ahmad, A., Azhar, M., Raya, A. Robust features of a QCD phase diagram through a contact interaction model for quarks: a view from the effective potential[J]. European Physical Journal A, 2023, 59(11): 252. doi: 10.1140/epja/s10050-023-01169-3
4. Corrêa, E.B.S.. Phase transition in a four-fermion interaction model under boundary conditions and electromagnetic effects[J]. Physical Review D, 2023, 108(7): 076002. doi: 10.1103/PhysRevD.108.076002
5. Tavares, W.R., Avancini, S.S., Farias, R.L.S. Quark matter under strong electric fields in the linear sigma model coupled with quarks[J]. Physical Review D, 2023, 108(1): 016017. doi: 10.1103/PhysRevD.108.016017
6. Ahmad, A., Murad, A. Color-flavor dependence of the Nambu-Jona-Lasinio model and QCD phase diagram[J]. Chinese Physics C, 2022, 46(8): 083109. doi: 10.1088/1674-1137/ac6cd8

Figures(22)

Get Citation
Aftab Ahmad. Chiral symmetry restoration and deconfinement in the contact interaction model of quarks with a parallel electric and magnetic fields[J]. Chinese Physics C. doi: 10.1088/1674-1137/abfb5f
Aftab Ahmad. Chiral symmetry restoration and deconfinement in the contact interaction model of quarks with a parallel electric and magnetic fields[J]. Chinese Physics C.  doi: 10.1088/1674-1137/abfb5f shu
Milestone
Received: 2021-03-13
Article Metric

Article Views(1926)
PDF Downloads(27)
Cited by(6)
Policy on re-use
To reuse of Open Access content published by CPC, for content published under the terms of the Creative Commons Attribution 3.0 license (“CC CY”), the users don’t need to request permission to copy, distribute and display the final published version of the article and to create derivative works, subject to appropriate attribution.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Chiral symmetry restoration and deconfinement in the contact interaction model of quarks with parallel electric and magnetic fields

    Corresponding author: Aftab Ahmad, aftabahmad@gu.edu.pk
  • Institute of Physics, Gomal University, 29220, D.I. Khan, Khyber Pakhtunkhaw, Pakistan

Abstract: We study the impact of steady, homogeneous, and external parallel electric and magnetic field strengths (eEeB) on the chiral symmetry breaking-restoration and confinement-deconfinement phase transition. We also sketch the phase diagram of quantum chromodynamics (QCD) at a finite temperature T and in the presence of background fields. The unified formalism for this study is based on the Schwinger-Dyson equations, symmetry preserving vector-vector contact interaction model of quarks, and an optimal time regularization scheme. At T=0, in the purely magnetic case (i.e., eE0), we observe the well-known magnetic catalysis effect. However, in a pure electric field background (eB0), the electric field tends to restore the chiral symmetry and deconfinement above the pseudo-critical electric field eEχ,Cc. In the presence of both eE and eB, we determine the magnetic catalysis effect in the particular region where eB dominates over eE, whereas we observe the chiral inhibition (or electric chiral rotation) effect when eE overshadows eB. At finite T, in the pure electric field case, the phenomenon of inverse electric catalysis appears to exist in the proposed model. Conversely, for a pure magnetic field background, we observe the magnetic catalysis effect in the mean-field approximation and inverse magnetic catalysis with eB-dependent coupling. The combined effects of eE and eB on the pseudo-critical Tχ,Cc yields an inverse electromagnetic catalysis, with and without an eB-dependent effective coupling of the model. The findings of this study agree well with the already predicted results obtained via lattice simulations and other reliable effective models of QCD.

    HTML

    I.   INTRODUCTION
    • Dynamical chiral symmetry breaking and confinement are the two fundamental properties of non-perturbative quantum chromodynamics (QCD). At zero or low temperatures T, the fundamental degrees of freedom of quantum chromodynamics (QCD) are the low energy hadrons, whereas at high T, the interaction gets increasingly screened and thus becomes weak, causing hadrons to melt down to a new phase in which the predominant degrees of freedom are quarks and gluons. This phenomenon is referred to as the confinement–deconfinement phase transition. As the strength of the QCD interaction decreases with increasing T, only the current quark mass survives when the parameter T exceeds a critical value. This is termed as the chiral symmetry breaking–chiral symmetry restoration phase transition. Such a phase transition is expected to have occured in the early universe, a few microseconds after the Big Bang, and it is experimentally observed in heavy-ion collisions at the Large Hadron Collider (LHC) in CERN and Relativistic Heavy Ion Collider (RHIC) at Brook Heaven National Laboratory (BNL). In addition, when the hadronic matter is subjected to an external electromagnetic field background, it yields a significant impact on phase transition. It is well known that at T=0, in the pure magnetic case, the strong magnetic field tends to strengthen the formation of quark anti-quark condensate, and the system remains in the chiral symmetry broken phase, even at a high magnetic field strength eB. This phenomenon is known as magnetic catalysis (MC) [1-10]. It was verified in earlier studies that at a finite T, the pseudo-critical temperature Tχ,Cc of the chiral symmetry restoration and deconfinement increases with an increase in eB; hence, magnetic catalysis is also observed at finite T [1, 2, 5, 11, 12]. Recently, the lattice QCD simulation [13-15] predicted that at finite T, a magnetic field suppresses the formation of quark-antiquark condensates and tends to restore the chiral symmetry at approximately the pseudo-critical temperature Tχ,Cc. Consequently, Tχ,Cc decreases with an increase in eB, and such a phenomenon is known as the inverse magnetic catalysis (IMC). This phenomenon is validated and supported by effective models of low energy QCD [16-23], as well as in holographic QCD models [24].

      In a pure electric case, and at T=0, the situation is relatively different from that of a pure magnetic field background. The strong electric field suppresses the formation of a quark-antiquark condensate, and thus tends to restore the chiral symmetry, i.e., the electric field anti-screens the strong interaction. Such a phenomenon is known as the chiral electric inhibition effect [1, 2, 5, 11, 12, 25-28] or the chiral electric rotation effect [29]. The nature of chiral phase transitions is of the second order in the chiral limit but cross-over when the bare quark mass is considered. At finite temperature, it is well understood that the pseudo-critical temperature Tχ,Cc decreases with an increase in electric field strength eE, and this phenomenon is known as the inverse electric catalysis (IEC) [28, 30]. The study on the influence of electric fields on the chiral phase transitions, as well as the effect of the magnetic field, is equally important from theoretical and experimental perspectives. Experimentally, in heavy-ion collisions, the electric and magnetic fields are generated with the same order of magnitude (1018 to1020Gauss) [31-34] in the event-by-event collisions using Au + Au at RHIC-BNL and in a non-central heavy-ion collision of Pb + Pb in ALICE-LHC. Moreover, some interesting anomalous effects, such as the chiral magnetic [35, 36], chiral electric separation [37, 38], and particle polarization [39-41] effects, which may arise owing to the generation of vector and/or axial currents in the presence of strong electromagnetic fields, are also required for theoretical exploration. Recently, a special case of considering the electric field strength parallel to the magnetic field strength (eEeB) significantly focused on exploring the above-mentioned phenomenon in the effective models of QCD [28, 29, 42, 43]. A major reason behind this approach is that the parallel electric and magnetic fields play important and prominent roles in several heavy-ion collision experiments [44, 45].

      Considering the aforementioned facts and findings, the objective of this study is to elucidate the traits of the dynamical chiral symmetry breaking-restoration transition in the presence of parallel electric and magnetic fields. The unified framework of this study is based on Schwinger-Dyson equations (SDE) in the rainbow-ladder truncation, Landau gauge, symmetry preserving the confining vector-vector contact interaction model (CI) [46], and Schwinger optimal time regularization scheme [47]. We adopt the quark-antiquark condensate ˉqq as an order parameter for the chiral symmetry breaking-restoration, whereas for the confinement-deconfinement transition, we use the confinement length scale [22, 48]. It should be noted that the chiral symmetry restoration and deconfinement occur simultaneously in this model [22, 49].

      This article is organized as follows. In Sec. II, we present the general formalism and contact interaction model at zero temperature, in the absence of background fields. In Sec. III, we discuss the gap equation at zero temperature, in the presence of parallel electric and magnetic fields presented in Sec. IV. Next, in Sec. V, we present the phase diagram at finite temperature, in the presence of parallel electric and magnetic fields. Finally, in Sec. VI, we present the summary and perspectives of this study.

    II.   GENERAL FORMALISM AND CONTACT INTERACTION MODEL
    • We begin with the Schwinger-Dyson's equations (SDE) for a dressed-quark propagator Sf, which is given by

      S1f(p)=iγp+mf+d4k(2π)4g2Δμν(pk)λa2γμSf(k)λa2Γν(p,k),

      (1)

      the subscript f represents the two light quark flavors, i.e., up (u) and down (d) quarks, g is the coupling constant, and mf is the current quark mass, which can be set to zero in the chiral limit. λa represents the conventional Gell-Mann matrices, Γν is the dressed quark-gluon vertex, and Δμν depicts the gluon propagator.

      In literature, it is well known that the properties of low energy hadrons can be reproduced by assuming that the interaction among the quarks occurs not via a mass-less vector boson exchange but by the symmetry preserving four-fermions vector-vector CI with a finite gluon mass [46, 50-53]:

      g2Δμν(k)=δμν4παirm2Gδμναeff,

      (2)

      where αir=0.93π represents the infrared enhanced interaction strength parameter, and mG=800 MeV is the gluon mass scale [54]. In the CI model, for a small value of αir and a large value of the gluon mass scale mG, there must be a critical value αeff; above this critical value, the chiral symmetry is broken, and below this value, there is a lesser chance for the generation of the dynamical mass. The d-dimensional (arbitrary space-time dimensions) dependence of this effective coupling and its critical value on the chiral symmetry breaking, using an iterative method in the superstrong regime, has been comprehensively studied in Ref. [55].

      The CI model Eq. (2), together with the choice of rainbow-ladder truncation Γν(p,k)=γν, forms the kernel of the quark SDE, Eq. (1), which brings the dressed-quark propagator into a very simple form [56]:

      S1f(p)=iγp+Mf.

      (3)

      This is possible because the wave function renormalization trivially tends to unity in this case, and the quark mass function Mf becomes momentum independent:

      Mf=mf+4αeff3Λd4k(2π)4Tr[Sf(k)].

      (4)

      In this truncation, the quark-anitquark condensate is given by

      ˉqqf=NcΛd4k(2π)4Tr[Sf(k)],

      (5)

      with Nc=3 representing the number of colors. The form of the proposed gap equation, Eq. (4), is significantly similar to the NJL model gap equation, except the coupling parameter αeff=92G [49].

      Further simplification of Eq. (4) yields the following gap equation.

      Mf=mf+16αeff3Λd4k(2π)4Mfk2+M2f,

      (6)

      where Mf represents the dynamical mass, and the symbol Λ stresses the need to regularize the integrals. Using d4k=(1/2)k2dk2sin2θdθsinϕdϕdψ, and performing trivial regular integrations with the variable s=k2, the above expression reduces to:

      Mf=mf+αeffMf8π20dsss+M2f.

      (7)

      The integral in Eq. (7) is not convergent; therefore, we need to regularize it via the proper-time regularization scheme [47]. In this scheme, we take the exponent of the integrand's denominator and then introduce an additional infrared cutoff, in addition to the conventional ultraviolet cut-off, which is widely adopted in NJL model studies. Accordingly, the confinement is implemented using an infrared cut-off [57]. By adopting this scheme, the quadratic and logarithmic divergences are eliminated, and the axial-vector Ward-Takahashi identity [58, 59] is satisfied. From Eq. (7), the denominator of the integrand is given by

      1s+M2f=0dτeτ(s+M2f)τ2irτ2uvdτeτ(s+M2f)=eτ2uv(s+M2f)eτ2ir(s+M2f)s+M2f.

      (8)

      Here, τ1uv=Λuv is an ultra-violet regulator that plays the dynamical role and sets the scale for all dimensional quantities. τ1ir=Λir represents the infrared regulator whose non zero value implements confinement by ensuring the absence of quarks production thresholds [60]. Hence, τir corresponds to the confinement scale [22]. From Eq. (8), it is clear that the location of the original pole is at s=M2, which is canceled by the numerator. Accordingly, we discarded the singularities, and the propagator is free from real and complex poles, which is consistent with the definition of confinement: "an excitation described by a pole-less propagator would never reach its mass-shell" [57].

      After performing integration over "s", the gap equation is given by

      Mf=mf+M3fαeff8π2Γ(1,τuvM2f,τirM2f),

      (9)

      where

      Γ(a,x1,x2)=Γ(a,x1)Γ(a,x2),

      (10)

      and Γ(a,x)=xtα1etdt, which is the incomplete Gamma function. By using the parameters of Ref. [50], i.e., τir=(0.24GeV)1 and τuv=(0.905GeV)1, with the bare quark mu=md=0.007 GeV, we obtain the dynamical mass Mu=Md=0.367 GeV and the condensate ˉqqu=ˉqqd=(0.0143) GeV3.

    III.   GAP EQUATION AT T=0 AND IN THE BACKGROUND OF PARALLEL eE AND eB
    • In this section, we study the gap equation in the presence of a uniform and homogeneous electromagnetic field with eEeB and at zero temperature. In the QCD Lagrangian, the interaction with the parallel electromagnetic field Aextμ embedded in the covariant derivative is expressed as

      Dμ=μiQfAextμ,

      (11)

      where Qf=(qu=+2/3,qd=1/3)e refers to the electric charges of u and d-quarks, respectively. We adopt the symmetric gauge vector potential Aextμ=(iEz,0,Bx,0) in Euclidean space, where both the electric and magnetic fields are selected along the z-axis. The gap equation in the presence of the parallel electromagnetic field remains in the form of Eq. (4), where Sf(k) is dressed with parallel background fields, i.e., Sf(k)~Sf(k). ~Sf(k), in the Schwinger proper time representation [1, 2, 5, 47], in the presence of the parallel magnetic field in Euclidean space [26, 29], is given by

      ~Sf(k)=0dτeτ(M2f+(k24+k23)tan(|QfE|τ)|QfE|τ+(k21+k22)tanh(|QfB|τ)|qfB|τ)×[γ4k4+Mf+tan(|QfE|τ)(γ4k3γ3k4)itanh(|QfBτ|)(γ1k2γ2k1)]×[1itanh(|QfB|τ)tan(|QfE|τ)γ5itanh(|QfB|τ)γ1γ2+tan(|QfE|τ)γ4γ3],

      (12)

      where the magnetic field couples with the coordinate indices 1 and 2, while the electric field couples with 3 and 4. Now, taking the trace of Eq. (12) and introducing both infrared and ultraviolet cut-offs, the gap equation in the presence of parallel electric and magnetic fields is given by

      ~Mf=mf+16αeff3f=u,d˜τ2irτ2uvdτ˜Mfeτ˜M2f×dk1(2π)dk2(2π)dk3(2π)dk4(2π)eτ((k24+k23)tan(|QfE|τ)|QfE|τ+(k21+k22)tanh(|QfB|τ)|QfB|τ.

      (13)

      The infrared cu-off τir is introduced in the model to mimic confinement by ensuring the absence of quarks production thresholds. In the presence of both eE and eB, it is required to vary slightly with both eE and eB. Therefore, the entanglement between dynamical chiral symmetry breaking and confinement is expressed using an explicit eE and eB-dependent regulator in the infrared [22]:

      ˜τir=τirMf˜Mf,

      (14)

      where Mf is the dynamical mass in the absence of background fields, and ˜Mf represents the eE and eB dependent dynamical mass. For chiral quarks (i.e., mf=0), the confining scale ˜τir diverges at the chiral symmetry restoration parameter near its critical value, which ensures the simultaneity of deconfinement and chiral symmetry restoration [22]. After the integration over k, the gap equation, Eq. (13), can be written as

      ˜Mf=mf+αeff3π2f=u,dˉτ2irτ2uvdτ˜Mfeτ˜M2f[|QfE|tan(|QfE|τ)|QfB|tanh(|QfB|τ)].

      (15)

      When both fields eE and eB0, Eq. (15) reduces to Eq. (9). The gap equation for the pure electric field can be obtained by setting eB0, while for pure magnetic field, eE0. The quark-antiquark condensate in the presence of background fields is in this form.

      In this study, we adopt two flavors f=2, i.e., u- and d-quarks. We use the same current quark mass for the up and down quarks, i.e., mu=md=0.007 GeV, such that the iso-spin symmetry is preserved. As is well-known, the response to the electromagnetic field is different for u and d-quarks because they have different electric charges [1, 5].

      The numerical solution of the gap equation Eq. (15) as a function of eB for fixed values of eE, is plotted inFig. 1. In the pure magnetic case (eE0), the dynamical mass ~Mf monotonically increases with an increase in eB , which ensures the magnetic catalysis phenomenon. The increase in ~Mf as a function of eB reduces in magnitude upon varying eE from its smaller to larger given values. We noticed that at eE0.33 GeV2, the dynamical mass ~Mf exhibits the de Haas-van Alphen oscillatory type behavior [61]. In other words, it remains constant for small eB values, then monotonically decreases with an increase in eB, and abruptly declines to lower values in the region eB[0.30.54] GeV2, where the chiral symmetry is partially restored via first-order phase transition; above it, the value increases again. This phenomenon can be attributed to the strong competition that occurs between parallel eB and eE, i.e., on the one hand, eB enhances the mass function, while on the other hand, eE suppresses it. We also sketch a zoom-in plot of the mass function, which shows the de Haas-van Alphen oscillatory behavior in the region eB[0.30.54] GeV2 with eE=0.33 GeV2 , as illustrated in Fig. 2. Such a behavior type is also explored and discussed in the other effective model of QCD, for example, refer to Refs. [29, 62].

      Figure 1.  (color online) Behavior of the dynamically generated quark mass, Eq. (15), as a function of the magnetic field strength eB at several given values of eE.

      Figure 2.  Zoom-in plot of the mass function that shows the de Haas-van Alphen oscillatory behavior in the magnetic field region eB=[0.30.5] GeV2 with the electric field eE=0.33 GeV2

      The behavior of the quark-antiquark condensate, Eq. (16), as a function of eB at various given values of eE, is illustrated in the Fig. 3. For a pure magnetic case (eE0), the magnetic field strength eB facilitates the formation of the quark-antiquark condensate. For given non zero values of eE, particularly, at eE0.33 GeV2 the evolution in the condensate is suppressed in the region eB=[0.30.54] GeV2. The nature of the transition in this region suddenly changes to first-order; however, above this region, it is enhanced with eB. The pseudo-critical values of the fields at which the chiral symmetry partially is restored and first-order phase transition occurres are eBc0.3 GeV2 and eEc=0.33 GeV2. Although such values of electric or magnetic field strength are sufficiently large in terms of what is typically generated during heavy-ion collisions, they may be relevant to astronomical objects, such as neutron stars, magnetars, etc. The confinement parameter ˜τ1ir as a function of eB for different fixed values of eE is presented in Fig. 4. We observed the same pseudo-critical fields eBc0.3 GeV2 and eEc=0.33 GeV2 for the confinement transition, similar to the case of chiral symmetry breaking.

      Figure 3.  (color online) Quark-antiquark condensate, Eq. (16), as a function of eB for several values of eE.

      Figure 4.  (color online) Confinement scale ˜τ1ir as a function of eB for several values of eE.

      In the following, we discuss the variation of ~Mf and ˜τ1ir as a function of eE, in the pure electric field, as well as for several non-zero values of eB. In Figs. 5, 6, and 7, we illustrate the behaviors of all three parameters as a function of eE, for various values of eB0.

      Figure 5.  (color online) Dynamical mass as a function of eE for several values of eB.

      Figure 6.  (color online) Quark-antiquark condensate, Eq. (16), as a function of eE for several values of eB.

      Figure 7.  (color online) Confinement scale ˜τ1ir as a function of eE for a given several values of eB.

      In the pure electric field case (eB0), the three-parameters monotonically decrease with an increase in eE, and at a pseudo-critical field strength eEχ,Cc , the chiral symmetry is partially restored and deconfinement transitions occurred; hence, the nature of the transitions here becomes smooth cross-over. Accordingly, we observed the chiral rotation or chiral electric inhibition effect in the contact interaction model, as already predicted by other effective models of QCD [1, 2, 25-29]. For non-zero eB values, we determined an interesting behavior for all the three parameters, as a function of eE. All three parameters are enhanced for given smaller to larger values of eB, except the chiral symmetry restoration and deconfinement regions, where all the parameters are suppressed by higher eB values. The pseudo-critical field strength Eχ,Cc decreases with an increase in eB, and at a pseudo-critical eBc0.3 GeV2, the transition changes from cross-over to first order. The nontrivial behaviors of all three parameters represent the competition between the magnetic catalysis and electric inhibition effects, both induced by eB in the presence of parallel eE [29]. As elucidated and validated in Ref. [25], the cause of the chiral inhibition effect is the second Lorentz invariant of the electromagnetic field, EB. The magnitude of eEχ,Cc at which the chiral symmetry restoration and deconfinement occurred are triggered from the inflection point of the electric gradients of andeE˜τ1ir, as presented in Figs. 8 and 9, respectively. In the pure electric case, the chiral symmetry restoration and deconfinement transition occur at eEχ,Cc0.34 GeV2. For several given values of eB0, eEχ,Cc decreases with an increase in eB. We determined a smooth cross-over phase transition up to the pseudo-critical magnetic field strength eBc0.3 GeV2, where the transitions take the first-order nature. Here, eEχ,Cc remains constant in the region eB[0.30.54] GeV2 and then increases with larger values of eB , as shown in Fig. 10. A similar behavior has already been demonstrated in [29]. The boundary point where the cross-over phase transition ends and the first order phase transition starts is known as a critical end point, and its co-ordinates are at (eBp0.3,eEp0.33) GeV2.

      Figure 8.  (color online) Electric gradient of the quark-antiquark condensate as a function of eE for fixed values of eB. At a particular fixed value of eB=0.3 GeV2, the electric gradient diverges at eEχc,0.33 GeV2 and above this value, the first order phase transition occurs.

      Figure 9.  (color online) Electric gradient of the confinement scale eE˜τ1ir, as a function of eE for different fixed values of eB. At a particular fixed value of eB=0.3 GeV2, the electric gradient of the confinement scale diverges at eECc0.33 GeV2.

      Figure 10.  (color online) Phase diagram of chiral symmetry and confinement for eEC,χc vs eB. eE=eEC,χc is obtained from the inflection points of the electric gradient of the condensate and the confining length scale eE˜τ1ir.

    IV.   QCD PHASE DIAGRAM AT T0 AND IN THE PRESENCE OF PARALLEL eE AND eB
    • In this section, we elucidate the behaviors of the dynamical mass, condensate, and confinement length scale at a finite temperature and in the presence of parallel electric and magnetic fields. We also explore the IEC, MC, and IMC phenomena, as well as the competition among them. Finally, we sketch the QCD phase diagram.

      The finite temperature version of the gap equation, Eq. (13), in the presence of parallel electric and magnetic fields can be obtained by adopting the standard convention for momentum integration:

      d4k(2π)4Tnd3k(2π)3,

      (17)

      and the four momenta k(ωn,k), with ωn=(2n+1)πT, represent the fermionic Matsubara frequencies. The Lorentz structure does not preserve anymore at finite temperature. By making the following replacements in Eq. (13),

      dk1(2π)dk2(2π)dk3(2π)dk4(2π)Tndk1(2π)dk2(2π)dk3(2π),

      and k4ωn, we have

      ^Mf=mf+16αeff3Tnf=u,dˆτ2irτ2uvdτ^MfeτˆM2f×dk1(2π)dk2(2π)dk3(2π)eτ((ω2n+k23)tan(|QfE|τ)|QfE|τ+(k21+k22)tanh(|QfB|τ)|QfB|τ),

      (18)

      with ˆMf=Mf(eE,eB,T). Performing sum over Matsubara frequencies and integrating over ks, the gap equation can be written as

      ˆMf=mf+αeff3π2f=u,dˆτ2irτ2uvdτˆMfeτˆM2fΘ3(π2,e|QfE|4T2tan(|QfE|τ))×[|QfE|tan(|QfE|τ)|QfB|tanh(|QfB|τ)],

      (19)

      where Θ3(π2,ex) is the third Jacobi's theta function. The confinement scale ˆτir here slightly varies with T, eE , and eB, and it takes the form

      ˆτir=τirMfˆMf.

      (20)

      Here, Mf=Mf(0,0,0) is the dynamical mass at eB=eE=T=0, whereas ˆMf=Mf(eE,eB,T) is the variation of the dynamical mass at finite eB, eE, and T. The quark-antiquark condensate is given by the following.

      First, we consider the case of a pure electric field at a finite temperature. The numerical solution of Eq. (19) as a function of T, for different given values of eE, is presented in Fig. 11. The dynamical mass ˆMf monotonically decreases with the increase of T until the dynamical chiral symmetry is partially restored. The response of eE is to suppress the dynamical mass. The quark-antiquark condensate, Eq. (14), and confinement length scale, Eq. (15), as a function of T for various given values of eE are depicted in Figs. 12 and 13, respectively. We infer that both parameters decrease with temperature T, and at a pseudo-critical temperature Tχ,Cc, the chiral symmetry is partially restored, and deconfinement occurs. We note that the electric field eE suppresses both parameters in the low-temperature regions and also reduces the pseudo-critical temperature Tχ,Cc. In Figs. 14 and 15, we plotted the thermal gradients of the condensate and the confinement length scale, respectively. The peaks in the thermal gradients of both parameters shift toward the low-temperature regions upon increasing the value of eE. Consequently, the critical temperature Tχ,Cc decreases with an increase in electric field strength eE; hence, inverse electric catalysis is observed at a finite temperature in the proposed contact interaction model. Therefore, the observations of this study are consistent with other effective models of QCD [28-30]. The magnitude of the critical temperature Tχ,Cc0.22 GeV2 is obtained from the inflection points of the thermal gradients of and Tˆτ1ir.

      Figure 11.  (color online) Dynamical mass, Eq. (19), as a function of temperature for various given values of electric field strength eE.

      Figure 12.  (color online) Quark-antiquark condensate, Eq. (21), plotted as a function of temperature for various given values of electric field eE.

      Figure 13.  (color online) Behavior of confinement scale, Eq. (20), plotted as a function of temperature for different given values eE.

      Figure 14.  (color online) Thermal gradient of the quark-antiquark condensate plotted as a function of temperature T for various values of the electric field eE. The peaks in the derivatives shift toward low temperature regions from small to large values of eE.

      Figure 15.  (color online) Thermal gradient of the confinement scale Tˆτ1ir plotted as a function of temperature T for various values of eE.

      Second, we consider the case of a pure magnetic field at finite temperature T. We plot the thermal gradient of the quark-antiquark condensate and the confinement length scale, as presented in Figs. 16 and 17, respectively. We note that the inflection points in both parameters shift toward the higher temperatures. In other words, eB enhances the pseudo-critical temperature Tχ,Cc of chiral symmetry restoration and deconfinement. Hence, in a pure magnetic case and at a finite temperature, we observe the magnetic catalysis phenomenon, which has already been observed in the CI model [22, 63].

      Figure 16.  (color online) Thermal gradient of the quark-antiquark condensate plotted as a function of temperature T for various values of magnetic field eB. From the plots, it is evident that the peaks shift toward higher temperatures.

      Figure 17.  (color online) Thermal gradient of the confinement scale Tˆτ1ir plotted as a function of temperature for several values of the magnetic field strength field eB.

      In most effective model calculations of QCD, it is well demonstrated that to reproduce the inverse magnetic catalysis effect as predicted by the lattice QCD [13, 14], the effective coupling must be taken as magnetic field dependent [17, 22, 28, 63] or both temperature and magnetic field dependent [16, 19, 23, 63, 64]. In the case of this study, we just employed the following functional form of the eB-dependent effective coupling αeff(eB) [22], where the coupling decreases with the magnetic field strength as

      αeff(κ)=αeff(1+aκ2+bκ31+cκ2+dκ4),

      (22)

      here κ=eB/Λ2QCD, with ΛQCD=0.24 GeV. The parameters a, b, c, and d were extracted to reproduce the behavior of critical temperature Tχ,Cc for the chiral symmetry restoration and deconfinement in the presence of magnetic field strength, obtained by lattice QCD simulations [13, 14]. The thermal gradients of the condensate and the confinement length scale, with the magnetic field dependent coupling, Eq. (22), is plotted in Figs. 18 and 19, respectively. It can be observed that the critical Tχ,Cc decreases with an increase in eB ; hence, the inverse magnetic catalysis can be observed at a finite T.

      Figure 18.  (color online) Thermal gradient of the condensate with magnetic field dependent coupling, Eq. (22), plotted as a function of temperature for several values of magnetic field strength. The IMC effect is depicted in this figure.

      Figure 19.  (color online) Behavior of Tˆτ1ir with the magnetic field dependent coupling, Eq. (22), as a function of temperature for several values of magnetic field strength eB.

      In Fig. 20, we sketch the combined phase diagram in the Tχ,CceE,eB planes, where we demonstrated the inverse electric catalysis, magnetic catalysis, and inverse magnetic catalysis (with magnetic dependent coupling). In the pure electric background, the solid-black curve indicates that the critical temperature Tχ,Cc decreases with an increase in eE; hence, the electric field strength inhibits the chiral symmetry breaking and confinement. In the pure magnetic limit (without magnetic field dependent coupling), the magnetic field eB enhances the critical temperature Tχ,Cc, and thus eB acts as a facilitator of chiral symmetry breaking and confinement. If we use the magnetic field dependent coupling, it can be observed that the critical temperature Tχ,Cc decreases as the magnetic field strength eB increases (red dotted-dashed curve). In this case, eB acts as an inhibitor of the chiral symmetry and confinement.

      Figure 20.  (color online) Combined phase diagram in the Tχ,Cc vs eE,eB planes of chiral symmetry breaking and confinement. The plot shows the evalution of IEC as well as MC without and IMC with eB-dependent coupling, Eq. (22). The Tχ,Cc values are obtained from the inflection points in Tˆσ and Tˆτ1ir.

      Third, we adopt both the non-zero values of eE and eB , and draw the phase diagram in the Tχ,CceE plane for various given values of eB , as shown in Fig. 21. At this point, the competition between IEC vs MC starts: the eB tends to catalyze the chiral symmetry breaking and confinement; consequently, Tχ,Cc is enhanced. In contrast, eE tries to inhibit the chiral phase transition and tends to suppres Tχ,Cc; finally, it is concluded with a combined inverse electromagnetic catalysis (IEMC). This may be different for a very strong eB, where the eB dominates over the eE. Next, we consider the eB-dependent coupling, Eq. (22), and sketch the phase diagram Tχ,Cc vs eE for various given values of eB in Fig. 22. It can be observed that the effect of parallel eE and eB with eB-dependent coupling tends to reduce Tχ,Cc. Therefore, both eE and eB produce IEC and IMC simultaneously.

      Figure 21.  (color online) Phase diagram in the Tχ,CceE plane of chiral symmetry breaking and confinement for various given values of eB. Here, eE inhibits the chiral phase transition, whereas eB facilitates it; consequently, the IEMC effect is observed.

      Figure 22.  (color online) Phase diagram in the Tχ,Cc vs eE of chiral symmetry breaking and confinement for various given values of eE with eB-dependent coupling Eq. (22). Here, both fields eE and eB inhibit the chiral phase transition; hence, the IEMC is observed.

    V.   SUMMERY AND PERSPECTIVES
    • In this work, we studied the influence of uniform, homogeneous, and external parallel electric and magnetic fields on the chiral phase transitions. In this context, we implemented the Schwinger-Dyson formulation of QCD, with a gap equation kernel comprising a symmetry preserving vector-vector contact interaction model of quarks in a rainbow-ladder truncation. Subsequently, we adopted the well-known Schwinger proper-time regularization procedure. The findings of this study are presented as follows.

      At zero temperature, in the pure magnetic background, the magnetic field facilitated the dynamical chiral symmetry breaking and confinement; hence, we observed the magnetic catalysis phenomenon. However, the electric field tends to restore chiral symmetry and deconfinement above the pseudo-critical electric field eEχ,Cc0.34 GeV2, i.e., the chiral rotation effect demonstrated in the proposed model. We observed that the electric field acted as an inhibitor of the chiral symmetry breaking and confinement. When both eE and eB were considered, we determined the magnetic catalysis effect for small given values of eE up to and above the critical electric field strength eEc0.33 GeV2, where the other parameters exhibited the de Haas-van Alphen oscillatory type behavior in the region eB=[0.30.54] GeV2. In this particular region,eE dominated over eB , and we observed the electric chiral inhibition effect, while above that region where eB was superior to eE, we observed magnetic catalysis again. We also realized that the pseudo-critical strength eEχ,Cc was suppressed with the increase in eB , and the transition nature was determined to be the smooth cross-over nature up to and above the pseudo-critical magnetic field strength eBχ,Cc0.3 GeV2 , where the transition suddenly changed to the first order transition. Subsequently, we located the position of the critical endpoint at (eBp0.3,eEp0.33) GeV2. We further determined that eEχ,Cc initially remains constant in the limited region eB=[0.30.54] GeV2 and then increases with larger values of eB.

      Finally, we sketched the phase diagram at a finite temperature and in the presence of parallel electric and magnetic fields. At a finite T, in the pure electric field limit, we determined that the pseudo-critical temperature decreased as we increasedeE; hence, the inverse electric catalysis was observed. However, for the pure magnetic field background, we observed the magnetic catalysis effect in the mean-field approximation and inverse magnetic catalysis with eB-dependent coupling. The combined effect of both eE and eB on Tχ,Cc yielded the inverse electromagnetic catalysis, with and without eBdependent effective coupling of the model.

      Qualitatively and quantitatively, the predictions of the presented CI-model agree well with results obtained from other effective QCD models, as well as modern Lattice QCD results. In the near future, we plan to extend this work to study the Schwinger pair production rate, including the dynamical chiral symmetry breaking for a higher number of colors, flavors, and in a parallel electromagnetic field. We are also interested in extending this work to study the IEMC phenomenon with the electric field and temperature-dependent coupling. Regarding the dependence of the model coupling on the considered electric field, there is no first-principle simulation for comparisons because lattice simulations in the presence of real electric background fields are not feasible to date. However, progress in this case has commenced, and results will be published somewhere else. The next strategy will be to study the properties of light hadrons in the background of parallel electric and magnetic fields.

    ACKNOWLEDGMENTS
    • The author acknowledges A. Bashir, A. Raya, R.L.S. Farias, and E. V. Gorbar for their valuable suggestions in the process of completing this work. The author also thanks his colleagues at the Institute of Physics, Gomal University, Pakistan.

Reference (64)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return