SciELO - Scientific Electronic Library Online

 
vol.68 número4Cosmological-static metric correspondence and Kruskal type solutions from symmetry transformations índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Revista mexicana de física

versión impresa ISSN 0035-001X

Rev. mex. fis. vol.68 no.4 México jul./ago. 2022  Epub 19-Mayo-2023

https://doi.org/10.31349/revmexfis.68.040101 

Research

Non-perturbative field theoretical aspects of graphene and related systems

J. A. Casimiro Olivares1 

A. Julia Mizher2 

A. Raya3 

1Instituto de Física y Matemáticas, Universidad Michoacana de San Nicolás de Hidalgo, Edificio C-3, Ciudad Universitaria, Morelia, Michoacán 58040, México.

2Instituto de Física Teórica, Estadual Paulista, Rua Dr. Bento Teobaldo Ferraz, 271 - Bloco II, 01140-070 São Paulo, SP, Brazil. Centro de Ciencias Exactas, Universidad del Bío-Bío, Avda. Andrés Bello 720, Casilla 447, 3800708, Chillán, Chile.

3Instituto de Física y Matemáticas, Universidad Michoacana de San Nicolás de Hidalgo, Edificio C-3, Ciudad Universitaria, Morelia, Michoacán 58040, México. Centro de Ciencias Exactas, Universidad del Bío-Bío, Avda. Andrés Bello 720, Casilla 447, 3800708, Chillán, Chile.


Abstract

In this article, we review the dynamics of charge carriers in graphene and related 2D systems from a quantum field theoretical point of view. By allowing the electromagnetic fields to propagate throughout space and constraining fermions to move on a 2D manifold, the effective theory of such systems becomes a non-local version of quantum electrodynamics (QED) dubbed in literature pseudo or reduced QED. We review some aspects of the theory assuming the coupling arbitrary in strength. In particular, we focus on the chiral symmetry breaking scenarios and the analytical structure of the fermion propagator in vacuum and under the influence of external agents like a heat bath, in the presence of a Chern-Simons term, anisotropy and in curved space. We briefly discuss the major advances and some new results on this field.

Keywords: Quantum field theoretical; electromagnetic fields; constraining fermions

1. Introduction

Two dimensional relativistic Dirac fermions have been the subject of active research for many decades. These are by no means merely flat cousins of quarks and leptons on the high energy physics realm, but rather these particles offer an opportunity to explore vast and intriguing phenomenology in other branches of physics as well (see, for instance, Ref.1). The pioneer work of Wallace2 on the band structure of graphene paved the way to consider relativistic Dirac fermions in solid state physics realms. The realization of the quantum anomaly emerging from the electromagnetic dynamics of Dirac fermions and photons in two spatial dimensions was first simulated in the seminal work of Semenoff3 under condensed-matter considerations, whereas one the first glimpses into topological matter was achieved from the dynamics of relativistic Dirac excitations in a Quantum Hall set up without Landau levels, where the role of parity anomaly is seen in the quantization of the conductivity as discussed by Haldane4. High temperature superconductivity in layered cuprates5-9 has also been naturally explained in therms of the dynamics of planar Dirac fermions in these systems. All of these ideas have highlighted the intricacies and interesting features of the interactions among Dirac particles constrained to move on a plane. But the first isolation of graphene membranes by the Cambridge group10 and others11 with the overwhelming evidence of the relativistic nature of charge carriers in this material12 have indeed boosted the interest toward these systems.

The relativistic behavior of this material is rooted in the organization of the carbon atoms in honeycomb lattices (see, for instance,13). This structure is conveniently represented by two overlapped triangular sublattices (with a bipartite unit cell), as represented in Fig. 1.

Figure 1 Crystallographic structure of the honeycomb array. In the left panel, the two overlapping triangular lattices are rep- resented by the red and green points. The primitive vectors a1 and a2 connecting all points in the crystal structure are shown, and the unit cell of the honeycomb array highlighted in blue. In the right panel, the first Brillouin zone is shown. The reciprocal lattice also shows a hexagonal structure. K and K’ represent the inequivalent Dirac points. 

The index labelling these sublattices constitute a degree of freedom for the charge carriers. In the tight-binding approach, electrons are firmly bounded to the atom except for one electron per atom belonging to the p orbital. Considering that the charge carrier can only hop to the nearest neighbor (which means that only inter-lattice hopping is allowed) one gets an expression for the energy that vanishes in six points, two of them non-equivalent. Those are called the Dirac points K and K’, located precisely where the valence band touches the conduction band and where the Fermi level is found. One can easily verify that in the vicinity of each of the Dirac points, the dispersion relation becomes linear, a feature typically associated to relativistic fermions. Besides that, the operators associated to the sublattice degree of freedom acquire a Pauli matrices structure in such a way that the system can be represented by a Dirac equation, where the speed of the light is replaced by the Fermi velocity v F , which turns out to be v F ≃ 1×106 m/s13. For this reason, an extensive work to characterize and explore graphene has been developed in the framework of quantum field theory (QFT). Moreover, the analytic structure of the Dirac equation unveils the fact that electrons in graphene have a chirality degree of freedom, which is another feature associated to relativistic particles and which has important implications for the electronic transport phenomena in this material.

Although in graphene the charge carriers are constrained to the plane, the external electromagnetic fields fermions may interact with are not. Because of that, usual Quantum Electrodynamics fully defined in (2+1)D (QED3), does not provide a suitable description. Among other inaccuracies, QED3 yields to a logarithmic static interaction (see, for instance, Ref.14) rather than the expected Coulomb interaction15-17 between electrons in graphene. Instead of that, in a QFT approach, electromagnetic fields must be represented by Abelian gauge boson fields in (3+1)D and the interacting theory therefore must deal with the subtlety of combining particles moving in different dimensionality.

Mixed dimensional gauge theories have been proposed in the nineties aiming to explain the then recent discovery of the Quantum Hall effect18. After that, several aspects of these theories were developed, in both formal and applied contexts. The first appearance of theories that provide a suitable background to explore fermions in two-space dimension interacting with gauge fields in three-space dimension occurred almost simultaneously in Refs.16-17. The procedure adopted in both approaches was to perform a dimensional reduction of the gauge field, integrating out its third spatial component, obtaining an effective interaction in (2+1)D. Among other features, this dimensional reduction softens the infrared behavior of the gauge boson propagator, which behaves as ~ 1/q 2 rather than the usual pole-dependence 1/q 2. In Ref.16 the authors obtain an effective Lagrangian in (2+1)D which for all intents and purposes corresponds to the mixed dimensional theory described above, naming the theory pseudo-quantum Electrodynamics (PQED). In particular, the authors focus on tracing a correspondence between PQED and the Chern-Simons theory. Further formal aspects of PQED were developed in Refs.19-20. On the other side, in Ref.17, the discussion is focused on the renormalization group aspects of graphene, and it is shown that the mixed dimensional theory possess a fixed point when the Fermi velocity evolves to the speed of light. The renormalization group analysis of this theory applied to graphene was posteriorly extended in Refs. 21-24. Recently, renormalization group techniques applied to PQED could successfully describe experimental data on the renormalization of the band gap in other two-dimensional materials like diselenide (WSe2) and molybdenumm disulfide (MoS2)25.

A more general procedure that could also be applied to systems where the charged particles are constrained to a lower dimensional string was proposed a few years later in Ref.26. The theory resulting from this method was named reduced-quantum electrodynamics (RQED). Although the framework developed in Ref.26 could be used to construct theories in arbitrary dimensions for the fermions d e and gauge bosons d γ (which then allow to label the extended theory as RQED dγ,de ), the major goal in Ref.26 was to analyse the dynamical chiral symmetry breaking in systems where two-spatially dimensional fermions interact with three-spatially dimensional gauge fields, focusing on condensed matter systems i . For this purpose, the authors solve the Schwinger-Dyson equations within an improved rainbow-ladder approximation and found that the dynamics of chiral symmetry breaking is rich and nontrivial. Posteriorly, the renormalization of RQED to one- and two-loops was deeply investigated27-30. The scale invariance obtained in previous approaches was confirmed in the context of RQED.

Following Ref.26, some activity involving the chiral symmetry breaking of dimensionally reduced theories took place. In Ref.31 the Schwinger-Dyson formalism within the rainbow approximation was applied to conclude that there is a critical coupling α = π/8 (in the conventions adopted in this work) above which there is room for dynamical mass generation in the theory. The existence of a critical coupling was confirmed in Ref.32. Additionally, a critical number of fermion families was determined from the similar structure of the gap equation in QED3 and RQED.

Chiral symmetry breaking was also explored in different external conditions: in Refs.33,34 the dynamical symmetry breaking was studied at finite temperature. In Refs.35,37 the dynamical mass generation was studied in RQED and PQED coupled to a Chern-Simons term. In Ref.38, anisotropy associated to strained graphene was shown to slightly affect the critical coupling. Moreover, in Ref.39 perturbative aspects of RQED in curved space was analysed.

Additionally, analysis of the Landau-Khalatnikov-Fradkin transformations (LKFT)40,41 was performed in RQED. These are non-perturbative transformations that connect the Green’s functions in different gauges and give valuable information about the renormalization coefficients in multi-loop calculations. Using the knowledge previously obtained in QED4 and QED342-44 the LKFT for the fermion propagator in RQED were derived in Ref.45.

In this manuscript we review the foundations of Pseudo- and Reduced-QED and focus on several aspects of dynamical chiral symmetry breaking explored via Schwinger-Dyson equations. In Sec. 2 we deduce the Lagrangian of the theory reviewing formal development as locality, unitarity and scale invariance. In Sec. 3 we review the formalism of Schwinger-Dyson equations, analyse general aspects of the gap equations, deduce the gap equations for RQED in the vacuum, analysing the wave function renormalization and renormalization group results. Section 4 is dedicated to study chiral symmetry breaking in a medium, including finite temperature effects, interaction with a Chern-Simons term and possible effects of curved space. In Sec. 5 we discuss how LKF transformations can provide information about renormalizability of the theory. Finally in Sec. 6 we present an outlook of the subjects discussed in this work.

2. Reduced or Pseudo QED: A historical recount

2.1. The Pseudo- or Reduced-QED Lagrangian

The derivation of PQED in Ref.16 and RQED in Ref.26 follow a similar procedure, consisting of performing a dimensional reduction of standard QED, imposing the matter fields to be constrained to a lower dimensional space. Although the general case, where an Abelian gauge field in arbitrary dimension interacts with fermions constrained to a generic manifold, has been explored in Refs.26,27, here we focus on the particular situation where four-dimensional QED is reduced to a three-dimensional theory that represents fermions in a plane.

The starting point is the usual QED4 Lagrangian,

L=-14FμνFμν-ejμAμ+LM+LGF, (1)

where the first term is the Maxwell term, j u is a matter current that couples to the gauge field A, e is the electric charge, ℒ M is a generic matter kinetic lagrangian and ℒ GF is a gauge fixing term.

In order to describe matter confined to a plane, the current must be defined in such a way that there is no dynamics in the third spatial coordinate,

jμ(x)={jμ(x0,x1,x2)δ(x3)μ=0,1,20μ=3. (2)

An effective interaction can be obtained integrating Eq. (1) over the gauge field Aµ This results in the generating functional (in Euclidean space)

Zeff=exp[e22d4xd4x'j3+1μ(x)Gμν(x-x')j3+1ν(x')]=exp[e22d4x d4x'j3+1μ(x)1-j3+1ν(x')]exp[-Seff(j3+1μ)] (3)

where □ denotes the d’Alambertian operator.

Considering the definition in Eq. (2) and the usual photon propagator,

Gμν=[-δμν+1-1ξμν]1(-)2, (4)

the effective action can be written as

Seff=-e22d3xd3x'jμ(x)KE(x-x')jμ(x'), (5)

where the (Euclidean) kernel is given by

KE(x-x')=d4k(2π)4eik(x-x')k2|x3=x3'=0=18π2|x-x'|2, (6)

and now x and x’ are defined in three dimensions.

It is possible to construct a theory fully defined in (2+1)D that mimics this effective action. This can be done noticing that the kernel in Eq. (6) can be written as a three-dimensional integral,

18π2|x-x'|2=14d3k(2π)2ei(k(x-x'))k2141-E, (7)

where here the d’Alambertian operator is defined in (2+1)D and the label E denotes Euclidean space. This yields the following effective action,

Seff=-e28d3xd3x'jμ(x)1-Ejμ(x'). (8)

It is straightforward to verify that the action in Eq. (8) can be obtained from the Lagrangian

LRQED=-14Fμν2Fμν-ejμAμ+LM+LGF. (9)

This is the Pseudo or Reduced-QED Lagrangian.

A remarkable consequence of the dimensional reduction, is that the photon propagator of RQED becomes proportional to 1/q rather than the usual 1/q 2 found in QED4. This can be obtained writing the pure gauge sector of the Lagrangian Eq. (9), including a proper gauge fixing term

LG=-12Fμν1Fμν+12ξaAa2 (10)

Defining the most general form for the propagator,

Δμν=a(k2)gμν+b(k2)kμkν (11)

it is possible to find the coefficients a(k2) and b(k2) (see e.g.26),

Δμν(q2)=-iq2gμν-1-ξqμqνq2. (12)

Note that this propagator has a softer infrared behavior than the photon propagator in QED4 and QED3. We notice that several groups differ in conventions by a global factor of 1/2 in this propagator.

The Feynman rules for RQED are represented in Fig. 2.

Figure 2 Feynman rules for RQED or PQED. 

2.2. Coulomb interaction

A crucial feature of this theory is that it correctly reproduces the Coulomb interaction expected for charged fermions, no matter if they are constrained to the plane. Naively considering full QED3, where the gauge fields are also in lower dimension, yields to a logarithmic interaction that does not correspond to what is observed in graphene15-17.

This can be seen considering a pair of static point charges with associated current

jμr,t= δ(r-x)+δ(r-y)μ=00,μ=1,2 (13)

The interaction energy associated to the current defined above interacting with a gauge field is given by

E=e2d2rjurAμr (14)

Rewriting it in terms of the gauge field propagator,

E=e22d2rd2r'dt'jμ(r)d2k(2π)2×dω2πeik(r-r')-ω(t-t')Gμν(ω,k)jν(r'). (15)

Integrating over t’ and ω and taking μ = v = 0 to obtain the potential between static points,

E=e22d2rd2r'ρrxd2k(2π)2eik(r-r')G00(ω=0,k)ρ(r'). (16)

Considering the current (13) and inserting the propagator in (12), one obtains, apart from unphysical self-interaction terms, that

E=e24π|x-y|, (17)

namely, the expected Coulomb interaction.

2.3. Scale invariance

In opposition to QED3, in RQED the electromagnetic coupling e2 is dimensionless, since it originates from the four-dimensional theory. This can also be seen doing a classical power counting of dimensions in (9). The theory is therefore classically scale invariant.

In Ref.17, a 1-loop calculation in a graphene motivated model, similar to RQED, showed that the system flows towards a Lorentz covariant point in the infrared. Perturbative calculations performed in the MS¯ scheme at 1 and 2-loops have shown that at this order the beta function vanishes and consequently the coupling in RQED does not run27,30, suggesting the absence of charge renormalization. A generalization of these previous results demonstrated that the RQED beta function vanishes at all orders46.

The renormalized fields and parameters can be defined in terms of dimensionless renormalization constants,

ψ=Zψ1/2ψr,A~=ZA1/2A~r,a~=ZAa~r,e=Zα1/2er,Γμ=ZΓΓrμ, (18)

where ψ is the fermion field, Ã is the reduced gauge field, a is a gauge fix parameter, Γµ is the fermion-photon vertex and the index r refers to renormalized. Because of finiteness of the vertex and electron charge, there is a relation between the renormalization constants,

Zα=(ZΓZψ)-2ZA-1. (19)

In what follows we briefly review the 1-loop calculation of the β-function performed in Ref.27. The 1-loop corrections to RQED propagators and vertex is represented in Fig. 3. For fermion fields constrained to a space in de dimensions and gauge fields in dγ dimensions, the diagrams in Fig. 3 are defined as,

iΠ1μν(q) =-ddek(2π)deTr(-ieγμ)ik+q(k+q)2(-ieγν)ikk2,  Σ1(p)=ddek(2π)de(-ieγμ)ik+q(p+k)2(-ieγν)i(4π)εeΓ(1-εe)(-k2)1-εe(gμν-ξ~kμkνk2), -ieΛ1μ=ddek(2π)dei(4π)εeΓ(1-εe)(-k2)1-εe(gνρ-ξ~kνkρk2)(-ieγν)ikk2(-ieγμ)ikk2(-ieγρ), (20)

where εγ=(4-dγ)/2, εe=(dγ-de)/2, ξ̃=ξ(1-εe), ε = 1 - a and a is a gauge fixing parameter. Integrating out the expressions (20), after some algebra one obtains27,

Π1(q2)=-de2(-q2)-εγ-εe(4π)de/21-εγ-εe3-2εγ-2εeG(1,1),ΣV1(p2)=-e2Γ(1-εe)(-p2)-εγ(4π)dγ/22(1-εγ-εe)22-2εγ-εe,-ξ(1-εγ-εe)G(1,1-εe)Λ1μ=e2Γ(1-εe)m-2εγ(4π)dγ/2γμ2(1-εγ-εe)22-εγ-εe-ξ¯Γ(εγ)Γ(2εe), (21)

where G(v1, v2) is the massless one-loop propagator defined by

G(ν1,ν2)=Γ(-de/2+ν1+ν2)Γ(de/2-ν1)Γ(de/2-ν2)Γ(ν1)Γ(ν2)Γ(de-ν1-ν2). (22)

For the configuration relevant to graphene and further two-dimensional Dirac/Weyl materials, d γ = 4, d e = 3, d γ e = 1/2 and ℰ→0. The β-function is defined as

β(α(μ))=dlogα(μ)dlogμ, (23)

and the renormalized coupling constant α is defined as

α(μ)4π=μ-2εγe2(4π)dγ/2Zα-1(α(μ))e-γεγ, (24)

where μ -2ℰγ is a factor to compensate the dimension of the coupling.

Figure 3 RQED 1-loop corrections for a) the photon propagator, b) the fermion propagator and c) the vertex. 

At this order, it can be obtained from (21) that ZΓZψ = 1, which indicates that the Ward identity is satisfied. Therefore, according to (19), the charge renormalization constant Zα depends exclusively on ZA. The particular MS¯ scheme chosen by the author is such that the finite terms are absorbed in the coupling and the renormalization constants reduces to unity for finite theories. Conversely, divergent theories are written as a Laurent series in ℰΥ. Within this framework, it can be extracted from (21) that for dγ = 4 and de = 3, Π1 is independent of ℰγ, which implies that

ZA=1+O(α2). (25)

Together with Eqs. (23) and (24), it yields to

β(α(μ))=-2εγ+γA(α(μ)). (26)

Since this theory has no anomalous dimension associated to the gauge field, and ℰ γ →0, the beta function vanishes.

A similar but more troublesome calculation can be performed at 2-loops27-30, showing that at this order the β-function remains vanishing.

In order to generalize this result to all orders, it is argued in Ref.30 that there is no renormalization of the gauge field since it would come from a non-local term in the free part of the action while counterterms must be local polynomials in the fields and their derivatives. Reminding that the charge renormalization depends only on Z A , this would imply in scale invariance at all-order. It is pointed out in Ref.46 however that this argument does not hold for the present theory, since it is valid only for renormalizable theories, which does not happen to be the case. For a discussion on this topic we refer to46.

An alternative procedure was proposed in order to generalize the result at 1-loop and show that it is valid at all orders46. First of all an effective theory for the gauge fields is obtained integrating out the fermions in the action, namely,

Γ~A=d3-εx12ZA2Fμν1-2Fμν+lnDet(iD/)+Sgf. (27)

Invoking gauge symmetry, it is possible to show that Γ̃(A) depends solely on the transverse projection of the gauge field,

AμT=δμν-μν2Aν. (28)

Expressing AμT in terms of strength tensor,

AνT=μ2Fμν=d3r4π(x-r)μ|x-r|3Fμνr. (29)

it is possible to define the all-order expansion of the effective action as,

Γ~=n1d3x1d3xnA1TAnTjμ1x1jμnxn=n1d3r1d3rnFμ1ν1r1Fμnνnrnγμ1ν1,,μnνnr1,,rn (30)

with

γμ1ν1,,μnνnr1,,rn=d3x14πd3xn4π(x1-r1)μ1|x1-r1|3×(xn-rn)μn|xn-rn|3jν1x1jνnxn (31)

By power counting and taking into account that the Furry theorem holds for RQED, the only contribution where divergences may arise is for n=2. This contribution corresponds to the usual transverse photon self-energy that leads to the following finite correction26,

Πμν(q)=e28qδμν-qμqνq2. (32)

The effective action for the gauge field then reads

Γ~A=d3-ϵx12ZA2Fμν1-Fμν+e28Fμν1-Fμν+O(e4F4-5)+Sgf (33)

which shows that the contributions from higher orders in the field strength are sufficiently ultraviolet-suppressed and give only finite corrections. Given that, it is safe to say that the renormalization constant ZA = 1 in RQED to all orders.

2.4. Renormalization group - RQED as a fixed point of a nonrelativistic model

Although one of the main achievements of RQED is to correctly describe several aspects of graphene, one point that deserves special attention is that usually the calculations are performed in a relativistic framework, where the Fermi velocity is taken to be the speed of light, or in natural units vF →1. Since the actual value of the Fermi velocity vF/c ≈ 1/300, the standard procedure of recovering vF in the end of the calculation must be done with extreme caution.

Previous to27,30, a renormalization group analysis was performed considering a mixed dimensional non-relativistic model17, where the Fermi velocity is taken into account. This breaking of Lorentz invariance can be explicitly appreciated in the action,

S=d3rψ¯(-γ00+vγ)ψ-ied3r(-γ0A0+vγA)ψ. (34)

Here, a procedure similar to the one adopted in RQED was performed, where the gauge field was previously integrated out in the third space coordinate. As expected, this yields to the same dependence of the photon propagator on 1/q found before.

In order to fully renormalize this model, the Fermi velocity must be also renormalized, and one more renormalization constant must be introduced, v0 = ZVZR. This explicitly modifies the fermion self-energy as:

Σ(ω,k)=Zψ(ω,k)#ωγ0-Zv(ω,k)vγk (35)

The one-loop calculation within this model was the first to predict the absence of charge renormalization for this type of mixed dimensional theory. Besides that, the renormalization constant for the Fermi velocity was obtained17,24,

Zv=1-116πe2vlnΛ. (36)

Defining an effective coupling

g=e24πvF, (37)

where g plays the role of a fine structure constant with v replacing the speed of light. The Fermi velocities at different energies are related by

v(E)=v(E0)1-g4lnEE0. (38)

It is possible to determine the existence of fixed points, imposing

βv(v,e2)=0. (39)

The Callan-Symanzik equation allows for obtaining an explicit form for the β-function17, and (39) is satisfied for vF = 1.

From this, one can conclude that the Fermi velocity grows as the energy decreases. Furthermore, the β-function has a non-trivial zero when the Fermi velocity reaches the speed of light. This limit corresponds to a Lorentz invariant weak coupling model whose coupling corresponds to the fine structure constant of QED. RQED therefore can be seen as a fixed point of the renormalization group to which the system represented by (34) flows. A detailed extension of renormalization group aspects of this mixed dimensional theory was performed in Refs. 21-24.

2.5. Causality and unitarity

As it usually happens to theories where some degrees of freedom are integrated out and the resulting theory is represented by an effective Lagrangian, RQED is nonlocal. Because of this, it becomes important to check if the theory respects causality and unitarity.

In fact, it has been shown that this mixed dimension theory respects both conditions. In Ref.19, the behavior of the classical Green’s function, advanced and retarded, is investigated inside and outside the light-cone. Generalizing the theory, considering a generic power α for the d’Alambertian operator in the denominator of the pure gauge term in the Lagrangian (9), it is shown that the Green’s function vanishes outside the light-cone. This implies that this family of theories respects causality for any value of the parameter α, including RQED, for which α = 1/2. Furthermore, for this theory the Green’s function vanishes inside the light-cone as well, being finite only on its surface. This means that RQED is constrained by a more strict condition than causality, it obeys Huygens principle19,47. While causality is found for any α, the Huygens principle applies exclusively to RQED, not holding for its cousin QED3.

Working with a more restrict generalization, where the power of the d’Alambertian operator in the RQED Lagrangian α = ENT#091;0,1), the validity of unitarity was explored in Ref.20. To this end, the authors invoke the optical theorem as follows. The S-matrix, that relates the initial and the final state in a scattering process, must be unitary so that the theory is unitary. Representing the S-matrix as S = 1 + iT, this implies that

iT-T=TT, (40)

which is the optical theorem. Evaluating the operator in the above equation between initial and final states, it can be represented in terms of the propagator,

i|T|f=(2π)3δ3(ki-kf)Dif, (41)

where a phase space factor was introduced to ensure a correct dimensionality. For i→ƒ the propagator in (41) becomes the Feynman propagator and performing a Fourier transform, (40) yields to

DF*(ω,k)-DF(ω,k)=-iTγDF*(ω,k)DF(ω,k). (42)

Here Τ comes from integrating the phase space factor and is the characteristic time of the system, and γ = -2(1- α). The generalized Feynman propagator can be easily obtained from the generalized Lagrangian,

DF(t,r)=dω2πd2k(2π)2e-iωteikr(ω2-k2+iε)1-α. (43)

A straightforward algebra shows that the optical theorem (40) is obeyed exclusively for α = 0, which corresponds to QED3, and for RQED, where α = 1/2, as long as a particular relation between Τ and ℰ is respected. Therefore, unitarity holds for these two theories, while for other values of the parameter α it does not.

All this reasoning was obtained for the free theory. In the case of interacting theory a similar procedure can be performed considering a dressed propagator in the random phase approximation (RPA),

Gμν=Gμα(0)δα,ν+ΠαβGβν(0)+ΠαβGβσ(0)ΠσγGγν(0)+.... (44)

Inserting in the above equation the 1-loop expression for Πμν calculated in Ref.48, yields to a result similar to the one for the free theory. Results beyond RPA were also obtained using the 2-loop expression calculated in Ref.27. In this case, the functional form for the propagator remains the same, the only modification being the coefficients. Therefore, up to 2-loops the optical theorem still holds, fulfilling the unitarity condition.

3. Chiral symmetry breaking

3.1. Schwinger-Dyson equations

Schwinger-Dyson equations (SDEs) are the field equations of a given quantum field theory (see, for instance49). These conform an infinite tower of relations among the Green functions involved: n-point functions are related to other n-point and higher-point functions, each verifying its own SDE. In their formal derivation, no assumption is made regarding the strength of the coupling constant(s) of the theory. Therefore, these equations are non-perturbative in nature and provide a useful tool to understand phenomena like bound-states, chiral symmetry breaking and confinement. The only systematic scheme to truncate the infinite tower of SDEs is perturbation theory, but in this scheme none of the above mentioned phenomena can be addressed reliably. In gauge theories such as QCD and QED, symmetry-preserving truncations of SDEs have made a tremendous development in addressing several non-perturbative aspects of these theories, which would otherwise require other frameworks to be addressed, including lattice field theory and effective model considerations.

SDEs are formally derived from the observation of the vanishing of the functional derivatives of the connected Green functions generating functional with respect to the fields. An alternative diagrammatic derivation of these equations can be derived directly from the Feynman rules of the theory under study. In QED, for instance (see Ref.50), the perturbative expansion of the fermion propagator, depicted in Fig. 4 shows three types of corrections, those to the fermion propagator itself (first row of corrections), those to the gauge boson propagator (second row of corrections) and the third kind corresponds to vertex corrections (third row of corrections).

Figure 4 Radiative corrections to the fermion propagator in QED. These include corrections to the fermion propagator itself, to the gauge boson propagator and vertex corrections. 

The infinite resummation of diagrams is better carried out by defining the fermion self-energy Σ(p), which is represented in the diagram in Fig. 5 and upon which the perturbative expansion of the fermion propagator contains all the radiative corrections shown in Fig. 4. In terms of Σ(p), the perturbative expansion of the fermion propagator is shown in Fig. 6.

Figure 5 Fermion self-energy. Blobs over the different parts of the diagram indicate that all perturbative corrections to these pares are already taken into account. 

Figure 6 Perturbative expansion of the fermion propagator in terms of the fermion self-energy. 

It corresponds to the expansion

Sp=S0(p)+S0(p)Σ(p)S0(p)\na+S0(p)Σ(p)S0(p)Σ(p)S0(p)+=S0(p)+S0(p)Σ(p)S0(p)+S0(p)Σ(p)S0(p)+=S0(p)+S0(p)Σ(p)S(p) (45)

The last line corresponds to the SDE for the fermion propagator depicted in Fig. 7.

Figure 7 SDE for the fermion propagator. 

It is convenient to re-write the above equation in terms of the inverse fermion propagator, which then becomes

S-1(p)=(S0(p))-1-Σ(p), (46)

and diagrammatically can be depicted as in Fig. 8.

Figure 8 SDE for the inverse fermion propagator. 

This equation involves the full photon propagator and fermion-photon vertex. The former, by a similar reasoning, can be seen to obey its own SDE depicted in the diagram in Fig. 9.

Figure 9 SDE for the gauge boson propagator. 

This two-point function is coupled to the fermion propagator and the fermion-photon vertex which verifies the SDE shown in Fig. 10.

Figure 10 SDE for the fermion-boson vertex. 

This diagram shows that the three-point function is coupled to two-point and four-point functions, hence illustrating the structure of the infinite tower of SDE.

At first glance, it seems impossible to envisage a non-perturbative truncation of the tower of SDEs without compromising the reliability of the predictions hence derived. Nevertheless, because of the gauge symmetry, a number of relations exist in which a (n + 1)-point Green function can be written on terms of n-point function. That is the case of Ward identities in QED51-54. Thus, one can attempt to propose a symmetry-preserving truncation and explore general features of the non-perturbative solution to the infinite tower of equations.

Below we shall review a favorite truncation to the SDE for the fermion propagator or gap equation, the so-called rainbow-ladder truncation. We also review how to improve of this truncation by including vacuum polarization effects and vertex corrections.

3.2. Structure of the gap equation

Let us consider the parity preserving ordinary version of QED. The SDE for the fermion propagator in an arbitrary number space-time dimensions d is given by

S-1p=S0-1p+4πiα ×ddk(2π)dΓμ(k,p)S(k)γνΔμν(p-k) (47)

where α = e2 / (4π) is the (dimensionless) fine-structure constant and e is the electric charge. The most general form of this propagator is commonly expressed as

S(p)=F(p2)p-M(p2), (48)

where F(p2) is the wavefunction renormalization function and M(p2) is the mass function. The tree level values of these functions are F(p2) = 1 and M(p2) = m0 where m0 is the fermion bare mass. Furthermore, Δμν (p - k) represents the full gauge boson propagator and Γμ (k,p) the full fermion-boson vertex. By neglecting vacuum polarization effects, i.e., working within the quenched approximation, we replace the full photon propagator Δμν (q) by its bare counterpart

Δμν(0)(q)=1q2gμν+(ξ-1)qμqνq2, (49)

where ξ is the covariant gauge parameter. In terms of diagrams, it amounts to consider Fig. 11.

Figure 11 SDE for the fermion propagator in quenched approximation. 

Under this approximation, the gap equation can be solved with a suitable choice of the fermion-boson vertex. The rainbow-ladder approximation (see50) corresponds to the perturbation-theory inspired choice Γμ (k,p) = γμ, which allows to decouple the SDE for the fermion propagator from the infinite tower, as it corresponds to the diagram in Fig. 12.

Figure 12 SDE for the fermion propagator in rainbow-ladder approximation. 

Chiral symmetry breaking is usually approached starting with m0 = 0. Thus, the natural consequence of this phenomenon corresponds to the dynamical appearance of fermion mass where there was none. Upon inserting the fermion and gauge boson propagators and fermion-photon vertex in the gap equation, after multiplying it by p and 1, respectively, we obtain the coupled system of equations

1F(p2)=1-4πiαd p2×ddk(2π)dTrpγμS(k)γνΔμν(0)(p-k),M(p2)F(p2)=4πiαd×ddk(2π)dTrγμS(k)γνΔμν(0)(p-k). (50)

Let us explore the solutions to this coupled set of equations in different space-time dimensions.

3.2.1. Gap equation in d = 4

It is interesting to consider the scenario of chiral symmetry breaking in ordinary QED. After taking the traces and Wick rotate to Euclidean space Eqs. (50), by switching to hyper-spherical coordinates, angular integrals can be performed analytically. Radial integrations are divergent. However, regulating these integrals with a hard cut-off in the momentum, we reach to the following pair equations to be solved self-consistently,

1F(p2)=1+αξ4π0Λ2dk2F(k2)k2+M2(k2)×k4p2Θ(p2-k2)-Θ(k2-p2), M(p2)F(p2)=α4π(3+ξ)0Λ2dk2F(k2)M(k2)k2+M2(k2)×k2p2Θ(p2-k2)-Θ(k2-p2) (51)

where Θ(x) is the Heaviside step function. In Landau gauge, ξ = 0, we have that F(p2) = 1 and thus we have to solve the following non-linear integral equation for the mass function,

Mp2=3α4π1p20p2dk2k2M(k2)k2+M2(k2)+p2Λ2dk2M(k2)k2+M2(k2) (52)

This equation has no analytical solution. Nevertheless, we can gain some insight by transforming it into the differential equation

ddp2p4dM(p2)dp2=-3α4πp2M(p2)p2+M2(p2). (53)

For p2 >> M2(p2), the above linearizes to

ddp2p4dM(p2)dp2=-3α4πM(p2). (54)

which is restricted to the boundary conditions

ddp2p2M(p2) |p2=Λ2=0,dM(p2)dp2|p2=κ2=0. (55)

The infrared cut-off κ2 is introduced to preserve the non-invariance of the original under M(p2) → cM (p2) with c constant as p2 → 0. Physically, it serves to quantify the amount of mass that is dynamically generated. We look for solutions to Eq. (54) of the form

M(p2)=(p2)s (56)

which after substituting in Eq. (54) yields

s(s+1)=-3α4π. (57)

Solving for s, we have that

s±=-12±1-3απ2. (58)

Thus, depending on the value of α compared to the value of αc = π / 3, we could observe an oscillatory behavior or a power-law solution. Consistency with boundary conditions demands α > αc. Moreover, the infrared and ultraviolet cut-offs are constrained as

Λκ=expπααc-1-2. (59)

This behavior is called the Miransky scaling law which states that the chiral symmetry breaking in QED corresponds to a conformal phase transition55. Thus, we conclude that in order for chiral symmetry breaking to be broken in ordinary QED, the coupling α must exceed a critical value.

3.2.2. Gap equation in QED3

A similar reasoning can be followed to explore the scenario of chiral symmetry breaking in QED restricted to a plane, namely, QED3. The main difference in this case is that the electric charge e2 has mass dimension one and hence serves as a natural scale for all mass scales in the theory. Furthermore, QED3 is super-renormalizable, and thus there is no need to regulate the integrals. In the rainbow approximation and in Landau gauge the gap equation, reduces to (see Ref.56)

M(p)=e22π2p0dkkM(k)k2+M(k)2lnk+pk-p. (60)

At this point, although the integrals are finite, one can introduce an ultraviolet regulator Λ in the integral with Λ >> e2 and consider at the end the limit Λ → ∞ maintaining e2/Λ fixed. Proceeding in this way, after expanding the logarithm for k >> p and p >> k, the integral reduces to

Mp=e2π20dkM(k)k2+M2(k)+e2π2p2pΛdkk2M(k)k2+M2(k) (61)

This equation can be transformed into the differential equation

ddpp3dM(p)dp=-2e2π2p2M(p)p2+M2(p), (62)

subjected to the boundary conditions

ddpp3dM(p)dp|p=00,M(p)|p=Λ0. (63)

Again, in the regime where p2 >> M2(p), the differential equation (62) linearizes in the following form

ddpp3dM(p)dp=-2e2π2M(p), (64)

which admits the general solution

Mp=4e2π2pC1J222e2π2p+C2Y222e2π2p (65)

where Jv (x) and Yv (x) are Bessel functions of the first and second kind of order v, respectively. UV boundary condition demands that C2 = 0. The constant C1 cannot be fixed from the IR conditions. Nevertheless, noticing again that the nonlinear is non-invariant under M(p) → cM(p) with c a constant, we require to change that boundary condition by demanding that

M(κ)=κ, (66)

in the IR, where κ is a regulator quantifying the amount of mass being generated. Such boundary condition allows to write the solution to the gap equation as57

M̃(x)=nα^xJ22α^x, (67)

With

M̃x=M(p/Λ)Λ,α^=e24πΛ, n=κ2α^1J22α^κ (68)

It is interesting that this solution is positive definite for κα^ but when κα^/20, the solution develops a zero. Furthermore, when κα^/40 a new zero develops, as illustrated in Fig. 13.

Figure 13 Multiple solutions to the gap equation in QED3. The solid black curve corresponds to M~x=x. The red long-dashed curve corresponds to a positive definite solution with κ=α^2, the dashed blue curve to κ=α^20 and the dotted green curve to κ=α^40

This pattern continues ad infinitum and emerges from the analytical properties of the gap equation, which corresponds to a Hammerstein equation of the first kind. This structure is particular of the truncation and is modified when some of the assumptions are removed from consideration. See Ref.57 fir further discussion.

3.2.3. Gap equation in QED3 including vacuum polarization effects

Including vacuum polarization effects in this model turns out interesting. Let us consider the scenario where we include a large number N of massless fermion families circulating in loops. This amounts to ressuming an infinite number of planar diagrams as shown in Fig. 14 such that one considers the unquenched approximation depicted in Fig. 15.

Figure 14 Vacuum polarization in the leading 1=N approximation. 

Figure 15 SDE for the fermion propagator in the leading 1/N approximation. 

In this case, the non-pertubative gauge boson propagator can be calculated exactly. In an arbitrary covariant gauge, it reads58-60

Δμν(q)=1q2+α̃q8gμν-qμqνq2+ξqμqνq4, (69)

where α̃=e2N. We observe that the infrared behavior of this propagator softens as 1/q when q→0 as compared to the usual 1/q2 pole structure in Eq. (49). It turns out that this softening makes the theory to loose its confining properties. It also has consequences regarding the chiral transition, as we can see from the gap equation. In the leading order of the 1/N approximation, F(p)=1+O(1/N) and the full fermion photon vertex Γμ=γμ+O(1/N). Working in Landau gauge, the mass function obeys

Mp=α~2π2Np0dkkM(k)k2+M2(k)×α~/8+8|k+p|α~/8+8|k-p| (70)

Once more, by expanding the logarithm for κ >> p and p >> κ, we have

Mp=α~π2Np0pdkkM(k)k2+M2(k)kp+α~/8+α~π2NppdkkM(k)k2+M2(k)pk+α~/8. (71)

This expression is equivalent to the differential equation

ddpp2(p+α̃/8)22p+α̃/8dM(p)dp=-α̃π2Np2M(p)p2+M2(p), (72)

subjected to the boundary conditions

0M(0)<,pdM(p)dp+M(p);p=α̃=0. (73)

In the large-coupling regime α̃p, the differential equation (72) simplifies to

ddpp2dM(p)dp=-8π2Np2M(p)p2+M2(p). (74)

Moreover, in the momentum domain p2M2(p), the above equation linearizes and takes the simplified form

ddpp2dM(p)dp=-8π2NM(p). (75)

Equation (75) admits a power law solution of the form

M(p)=ps, (76)

which upon substitution into (75) yields

s±=-12±121-32π2N. (77)

Thus, the solution to the linearized equation is

Mp=C1p-12+121-32π2N+C2p-12-121-32π2N. (78)

Again, we notice the existence of a critical number of fermion families, Nc = 32/π2 that distinguishes an oscillatory from a decaying behavior of the mass function. Consistency with boundary conditions suggests that Nc > 32/π2 and the amount of dynamically generated mass, quantified as M(p = 0) also follows a Miransky scaling law of the form

M(0)=α̃exp-2πNcN-1+δ. (79)

Thus, this scenario has a similar scaling behavior as ordinary QED. Furthermore, vacuum polarization effects translate to the existence of a critical (large) number of fermion families above which chiral symmetry breaking is no longer possible and the chiral symmetry restoration corresponds to a conformal phase transition55.

It is interesting that the gauge dependence of the mass function and wavefunction renormalization have been considered in covariant gauges61 and up to the next-to-leading order of the approximation62-65. These ressumations carried out in an arbitrary non-local gauge do not spoil the existence of the critical value Nc for the restoration of chiral symmetry and only refine the value of this critical number Nc. The conclusion up to date is that such number is an integer Nc ≤ 3.

3.2.4 Gap equation in RQED

So far we have reviewed different scenarios for chiral symmetry breaking in ordinary QED. The gap equation in this theory is such that the wavefunction renormalization is trivial in Landau gauge. The scenario of criticality emerges in d = 4 as the need for the coupling to exceed a critical value for the theory to be able to break dynamically the chiral symmetry. In three dimensions, the quenched theory does not exhibit this feature. Chiral symmetry breaking can be broken for arbitrary values of the coupling. Nevertheless, including vacuum polarization effects, it is observed that if the number of fermion families circulating in loops exceeds a critical value, chiral symmetry cannot be broken. This statement is valid up to 1/N2 and in arbitrary non-local gauge. An interesting observation is that when fermions in loops remain massless, there is a softening of the infrared pole in the propagator 1/q2 as the gauge boson momentum q→0 to a 1/q softer behavior. This behavior is similar to the tree-level gauge boson propagator in RQED

Δμν(0)(q)=12qgμν-(1-ζ)qμqνq2. (80)

Notice that in this case, the gauge parameter ζ is not the same as ξ in the ordinary theory, but here it has to be renormalized because fermions and bosons live in different dimensions. In fact, ζ = ξ / 2, and hence ζ = 0 corresponds to Landau gauge.

Following the reasoning of the previous subsections, we truncate the gap equation in the rainbow approximation as in Ref.31. In an arbitrary covariant gauge, it is equivalent to the following coupled system of equations (see66)

1F(p)=1+α2πp20ΛdkF(k)k2k2+M(k)2θ(k-p)k-(2+ζ)p2+(1-ζ)p4k2+θ(p-k)p-(2+ζ)k2+(1-ζ)k4p2,M(p)F(p)=α(2+ζ)2π0Λdkk2F(k)M(k)k2+M2(k)θ(k-p)k+θ(p-k)p (81)

with the dimensionless coupling α = e2 / 4π as usual. It is evident from the above expressions that for ζ = 0, F(p) ≠ 1. Nevertheless, as a first approximation, let us take F(p) = 1 and explore the solution of the gap equation. Then, the gap equation becomes the non-linear integral equation for the mass function31,66

Mp=2απp0pdkk2M(k)k2+M2(k)+2απpΛdkk2M(k)&#091;k2+M2(k)&#093;k (82)

This expression can be straightforwardly converted into the following differential equation

p2M(p)+2pM'(p)+2απp2M(p)p2+M2(p)=0, (83)

restricted to the boundary conditions

limpΛpdM(p)dp+M(p)=0,limp0p2dM(p)dp=0. (84)

Upon linearizing the differential equation 83 when p >> M(p), we can again write the resulting equation in the Euler-Cauchy form

ddpp2dM(p)dp+απM(p)=0, (85)

which, admits a general solution of the form

M(p)=C1pn++C2pn-, (86)

where

n±=-12±121-8απ. (87)

As in previous cases, the non-invariance of the gap equation under scalings M(p) → cM(p) where c is a constant demands the introduction of an infrared cut-off κ and the modification of the infrared boundary condition to

M(κ)=κ. (88)

Thus, defining αc = π / 8, boundary conditions demand that α > αc. Furthermore, the dynamical mass obeys the Miransky scaling

Λκ=expAααc-1+δ, (89)

with A = -2π and ∂ = -4, hence indicating that the chiral symmetry is broken in a conformal phase transition provided the coupling exceeds a critical value.

In order to explore the effect of the wavefunction renormalization in the result89, particularly in relation with the gauge dependence, below we explore some variants of the rainbow-ladder truncation following66.

3.2.5. Effect of the wavefunction renormalization

To explore the sensitivity of the truncation to the effect of the wavefunction renormalization function in connection with the gauge parameter dependence of the critical coupling and the dynamical mass, we conduct the following exercise:

First, we impose F(p) = 1 in all covariant gauges and solve the equation

Mp=α2πp0pdkk2M(k)k2+M2(k)2k&#091;2+ζ&#093;k+α2πppΛdkk2M(k)&#091;k2+M2(k)&#093;k2p2+ζ, (90)

for various values of the gauge parameter ζ.

As a second variant, we solve the coupled system of equations in (81) for various values of the gauge fixing parameter.

Next we impose the Ward-Takahashi identity (WTI),

(k-p)μΓμ(k,p)=SF-1(k)-SF-1(p), (91)

which allow to split the vertex into its longitudinal and transverse pieces,

Γμ(k,p)=ΓLμ(k,p)+ΓTμ(k,p),

where (k-p)μΓTμ(k,p)=0. We exploit this identity into the gap equation in the following way. Assuming that all the dependence on the gauge parameter arises only from the photon propagator, we split this two-point function into its transverse component and the gauge parameter dependent longitudinal part,

Δμν(0)(q)=ΔμνT(q)+ζqμqνq3, (92)

the gap equation can be written as

S-1p=S0-1p+4iαπd3k(2π)3γμSkΓνΔμνTq+4iαπζd3k(2π)3γμSkΓνqνqμq3, (93)

with q = k - p. In the final term we replace the identity (91) and upon taking traces after multiplying by 1 and p, respectively, we find the gap equation to be equivalent to the coupled system of equations

1F(p)=1+α2πp20ΛdkF(k)k2k2+M(k)2θ(k-p)k-2p2-ζp4k2+θ(p-k)p-2k2-k4p2+αζ4πp20λdkF(k)k2F(p)(k2+M2(p)×θ(k-p)k2p2k2k2-p4k2-MkMpp2k2+2p2MkMpk2+θp-kp2p2k2p2-p2k2p2-MkMpk2p2+2p2M(k)M(pp2M(p)F(p)=απ0Λk2F(k)M(k)k2+M2(k)θ(k-p)k+θ(p-k)p+αζ2π0λdkF(k)k2F(p)(k2+M2(p)×θ(k-p)kM(p)(1-p2k2)-M(k)p22k2+θ(p-k)p-M(k)(1-k2p2)-M(p)k22p2 (94)

Then we solve the coupled system of equations for various values of ζ.

As a final variant of the gap equation, we use the Ball and Chiu vertex67, which is explicitly constructed to verify the Ward identity and is written as follows

ΓBCμ=γμ21F(k)+1F(p)+12(k+p)(k+p)μ(k2-p2)1F(k)-1F(p)+(k+p)μ(k2-p2)M(k)F(k)-M(p)F(p) (95)

This is a standard choice for the longitudinal piece of the vertex and it is explicitly constructed to satisfy these relations. Nevertheless, the additional terms to the central part exhibit spurious kinematic singularities as k2p2. These spurious singularities do not appear in perturbation theory and are unwanted in a non-perturbative construction of the vertex. These actually are cancelled by more educated ansatze of the vertex which include the transverse piece unconstrained by these identities. Thus, for our purposes, by keeping only the central part, which is

ΓCBCμ=γμ21F(k)+1F(p), (96)

we guarantee that the ward identities are satisfied up to spurious singular terms which are not expected to appear in a complete form of the vertex.

Then, the gap equation corresponds to the following coupled system of equations

1F(p)=1+α2πp20ΛdkF(k)k2k2+M(k)21F(k)+1F(p)θ(k-p)k-(2+ζ)p2+(1-ζ)p4k2+θ(p-k)p×-(2+ζ)k2+(1-ζ)k4p2,M(p)F(p)=α(2+ζ)2π0Λdkk2F(k)M(k)k2+M2(k)1F(k)+1F(p)θ(k-p)k+θ(p-k)p (97)

In Fig. 16 we compare the findings of the critical coupling in different gauges for the variants of the truncation described above.

Figure 16 Critical coupling as a function of the covariant gauge parameter for different models. 

At first glance, it might seem surprising that the simplest scenario of neglecting wavefunction renormalization effects and solving the mass function alone yields the least gauge parameter dependent results. This is better understood from a renormalization group analysis, that exploits the similarities of the 1/N approximation in QED3 and RQED itself, as discussed below.

3.3. Renormalization Group Analysis

It is remarkable that RQED becomes an infra-red Lorentz invariant fixed-point for a general theory of mixed dimensions for fermions and gauge bosons, RQED dγ,de 17, because it allows to analyze the critical behavior of the general theory, regarding chiral symmetry breaking, around such point. The important observation is that as we approach this fixed point, the Fermi velocity of fermion fields tends to the speed of light in vacuum and at the same time, the coupling constant approaches the fine-structure constant of QED. Therefore, one can exploit these facts to compare the critical structure of RQED and QED332,62-65,68 This is particularly relevant to understand the gauge dependence of the critical numbers associated to the chiral transition.

In order to make explicit the mapping between QED3 and RQED, we start by expressing the gauge boson propagator in RQED as32

Δμν(RQED)(q)=12qgμν+(ζ-1)qμqνq2, (98)

where we recall that the gauge fixing parameter ζ of RQED is half of the corresponding parameter of QED. On the other hand, the leading behavior of the gauge boson propagator (69) of QED3 in the large N approximation has the form

Δμν(QED3)(q)=8Ne2qgμν+(ξ-1)qμqνq2. (99)

Thus, identifying

1π2Nα4π,ζξ2, (100)

we can simply map the solutions of the gap equation in RQED to those of QED3 in the 1/N approximation.

For instance, the gauge parameter dependence of the critical number of fermion families for chiral symmetry restoration in QED3 was explored in Refs.32,68. The authors find that the critical condition can be written as

1=16(2+ξ)Lc. (101)

Then, from (100) we can straightforwardly find that for RQED, this result is translated to the critical coupling

1=16(2+ξ)αc4π, (102)

which shows a strong gauge dependence of the critical coupling, although it does not depend on the flavor number N. This result comes up as a consequence of a partial re-summation of diagrams in the approximation. It can be improved, for instance, by considering a RPA calculation of the vacuum polarization tensor. In the RPA, the gauge boson propagator has the leading behavior

Δμν(RQED)(q)=12q(1+Ne2/16)×gμν+(ζ-1)qμqνq2. (103)

Thus, by redefining the coupling as

αα'=α1+e2N/16, (104)

the gauge dependence of the critical coupling is now

αc=π2(5+ξ¯)-N2π/4, (105)

with ξ=(1+ξ¯)/2 and N = 2 for graphene. The gauge dependence of αc is milder in this case.

At next-to-leading order, one observes that the gauge dependence of the critical coupling has the form

αc=4π8(2+ξ)+d(x), (106)

where

d(ξ)=8S(ξ)-84-1123ξ+9ξ2-4Nπ2, (107)

and

Sξ=(1-ξ)R1-(1-ξ2)R28-(7+16ξ-3ξ2)P2128, (108)

with the coefficients

R1=163.7428,R2=209.175,P2=1260.720. (109)

This is a rather intricate form of the coupling which nevertheless improves upon the leading order gauge parameter dependence. The RPA still helps to reduce the resudial gauge dependence of the critical coupling by redefining

αα'=α1-αN4π2. (110)

Thus, we observe that in this framework, the appropriate re-interpretation of the coupling permits a less severe gauge dependence of the critical value of the coupling to trigger chiral symmetry breaking.

4. Chiral symmetry breaking in a medium

Mixed-dimensional theories allow natural extensions of the QFT formalism to incorporate the effects of (classical and quantum) external agents like a heat bath and/or external electromagnetic, strain or gravitational fields, among others. Generally speaking, these effects have the potential to enhance or inhibit phase transitions. In the case of the chiral transition, a heat bath is known to have the effect of restoring chiral symmetry. The imaginary-time formalism of thermal field theory (TFT) is a natural extension to consider in this case. On the contrary, electromagnetic fields promote the breaking this symmetry. In particular, in 3D materials, the configuration of parallel electric and magnetic fields gives rise to a special configuration that can be understood as the presence of a Chern-Simons (CS) term in the electromagnetic Lagrangian, with the added possibility of parity and time reversal symmetry breaking. In 2D systems, the presence of such a term for fermion fields is realized as the possibility of mass terms different from the ordinary Dirac mass.

A Haldane term of this type has interesting effects on the chiral phase transition in ordinary QED3. Curvature effects are also important in modifying the electric and optic properties of the materials and are generally expected to occur in materials with defects. In this section, we review the effect of a thermal bath, the role of a Chern-Simons term, the influence of strain and the formulation of RQED on curved space. We specialize the discussion in how these effects impact the scenario for chiral symmetry breaking.

4.1. Finite temperature effects

The scenario of chiral symmetry breaking in RQED has been considered at finite temperature. Let us revisit the calculations in Refs.26,31,33,34. We start our discussion from the imaginary time formalism of TFT. We introduce the fermionic Matsubara frequencies ω n = (2n +1) πT and replace of any integral over the temporal component of an arbitrary three-vector p = (p0, p) by the summation over frequencies,

dp0f(p0,p)Tn=-f(-iωn,p). (111)

Under this prescription, working in the Landau-like gauge ζ = 0 and neglecting wavefunction renormalization effects, the gap equation takes the form33

Mmp=n=-Λd2k(2π)2(4παT)Mn(k)(2n+1)2π2T2+k2+Mn(k)2×1&#091;4(m-n)2π2T2+(p-k)2&#093;1/2 (112)

where Mn (p) = M(ωn, p) is the mass function associated to every Matsubara frequency and the notation ?Λ indicates that divergent integrals are to be regularized with an ultraviolet cut-off Λ.

Several considerations are at hand. First of all, as customary in TFT, one would be tempted to perform the sum over Matsubara frequencies and introduce the cut-off in momentum integrals. Nevertheless, because these frequencies appear inside the square-root, the standard techniques to perform such sums are hard to implement. First of all, the convergence of the sum

U=n=-un=n=-1n2+A21(2n+1)2+B2=2n=1un(+)+u0 (113)

is ensured (see Ref.33) provided the sum in the last row of the above expression converges. In Eq. (33), it is assumed that A2, B2 > 0 and un(+)=(un+u-n)/2. The proof of convergence is established via the auxiliary sum

V=n=1vn=n=11(n+A)(n+B), (114)

which is convergent and satisfies that vnun(+)0 and

limn=un(+)vn=0. (115)

Then, because V is convergent, so it is nun(+) and thus, U is also convergent.

Alternatively, one can perform the momentum integrals first and then the sum over Matsubara frequencies. In this case, one can assume that the mass function is the same for all the frequencies such that the gap equation can be written as33

M(p)=4αTπ0dxM(k)I(k), (116)

where

Ik=n=-1(2n+1)2π2T2+k2+M(k)2×K(xn)(p-k)2+4n2π2T2. (117)

K(x) denoting an elliptic function and

xn=-4kp(k-p)2+4n2π2T2. (118)

Considering only the zeroth Matsubara frequency n = 0, an approximation that is commonly used in high-temperature regime studies, the gap equation reduces to33

M(p)=4αTπ0dkk|k-p|K(x0)M(k)π2T2+k2+M2(k). (119)

Expanding kK(x0) / |k - p| for k >> p and k << p, the gap equation becomes

Mp=2αTp0pdkkM(k)π2T2+k2+M2(k)+2αTpΛdkM(k)π2T2+k2+M2(k). (120)

Upon linearization, the resulting integral equation can be cast in the form of the differential equation

ddpp2dM(p)dp+2αTpM(p)=0, (121)

subject to the boundary conditions

limpπTcp2dM(p)dp=0,limpΛM(p)=0, (122)

where Tc would correspond to the critical temperature for chiral symmetry restoration, which will be defined from the analytical behavior of the mass function. The general solution to Eq. (121) is expressed as

Mp=C1αTpJ18αTp+C2αTpY18αTp, (123)

where J1(x) and Y1(x) are Bessel functions of the first and second kind. UV boundary condition imply C2 = 0. At the same time, the IR boundary condition imposes the following relation

2J1(ξ)+ξJ0(ξ)-ξJ2(ξ)=0, (124)

with ξ=8αT/πTc. The above relation has nontrivial solutions for a given n such that {ξ n } = {ξ0, ξ1, …} with the hierarchy ξ n < ξ n + 1 for any given n. The lowest solution ξ0 ≃ 2.4 fixes the value of the critical temperature T C = 0.44 αT. This result establishes that even at low temperatures, the value

αcTαcξ02, (125)

exceeds the critical value of the coupling in vacuum33. It has it s origin in the fact that there is no evidence that the mass function M(p) vanishes at the critical temperature. Therefore, a more refined treatment is called for.

One can improve on the above result recalling that the phenomenon of dynamical chiral symmetry breaking is infrared. Therefore, one can simply neglect the external momentum in the gap equation (119), which then becomes33

1=2αT0Λdkπ2T2+k2+m2(T), (126)

where we have introduced the shorthand notation m(T) = M(0, T). Retaining Λ as the largest scale of the problem, two solutions are found for the above expression, namely

m(T)=±πTc2-T2,Tc=α2T2Λ2(2αT+Λ)2. (127)

Keeping the positive solution alone, one can derive the critical coupling for chiral symmetry breaking of the form

αc(T)=ΛΛ-2T. (128)

Nevertheless, this expression is oblivious to the critical value of α at T = 0.

A different view of the gap equation can be achieved as follows34. Let us introduce dimensionless quantities,

p= Λσ,k=Λρ, T=ΛT~, Mmk=ΛM~mρ, (129)

such that the gap Eq. (112) can be written as

M~mσ=αT~n=-Nf-1Nf01dρπdθρ×M~m(ρ)(2n+1)2π2T~2+ρ2+M~n2(ρ)×1|4(m-n)2π2T~2+(σ-ρ)2|1/2, (130)

where θ is the angle between σ and ρ ii Notice that in this case, the cut-off is not introduced in momentum integrals, but in the number of frequencies Nƒ that are summed up. Whithin the so-called constant mass approximation69, all mass functions involved in the gap equation can be replaced with their values at zero momentum. Denoting M~m(σ)Sm, the gap equation under this approximation simplifies to

Sm=2αT~n=-Nf-1Nf01dρρSn(2n+1)2π2T~2+ρ2+Sn2×1|4(m-n)2π2T~2+ρ2|1/2 (131)

Fixing the the cut-off Nƒ such that

(2Nf+1)πT0=Λ(2Nf+1)πT~0=1\naT~0=1(2Nf+1)π, (132)

we are allowed to write all temperature scales as T=kT̃0,kN. Thus, performing momentum integration analytically and assuming a linear scaling of Smγm (T - TC) near criticality,

Sm2αT~n=-Nf-1NfSnπ2T~2-4m2+8mn+4n+1×tan-1π2T~2m-n2+1π2T~2-m-n4m2+8mn+4n+1-tan-12m-n4m2+8mn+4n+1 (133)

Specializing on the zeroth Matsubara frequency m = 0, observing that Sn/S0 ≤ 1, and recalling our assumption that near the critical point the temperature-dependent Sm (Ť)→0 we have that34

1αc=n=-Nf-1Nf2π4n+1tan-14π2T~c2n2+1π2T~c24n+1-tan-12n4n+1. (134)

The sum over n is finite for every value of temperature. This allows to obtain the behavior of the critical coupling, αc, as shown in Fig. 17.

Figure 17 Critical coupling as a function of the critical temperature. 

The behavior of the critical coupling as a function of temperature is consistent with

αc(T̃)=π8+aT̃+bT̃2, (135)

with a,b real numbers. This behavior exhibits the correct behavior of αc(T̃=0)=π/8 34.

It is remarkable that within the finite temperature framework, the gauge dependence of the critical temperature and coupling is a pending matter up to date.

4.2. Chiral symmetry breaking with a Chern-Simons term

Physics on two spatial dimensions is very interesting in its own right. For electromagnetic interactions, there exists the possibility to include a Chern-Simons (CS) term in the gauge sector through the Lagrangian

L=-iϑ4ϵμνρAμνAρ-AμJμ, (136)

where ϑ is the CS coefficient (see, for instance, Ref.70). Such a term induces a topological mass to the gauge bosons. It also induces fractional (anyon) statistics for matter particles as well as parity and time-reversal symmetry breaking of the theory, among other interesting effects. It is through this term that, for instance, the fractional quantum Hall effect can be explained and the emergence of the transverse conductivity without Landau levels can be explained.

The scenario of chiral symmetry breaking in QED3 with a CS has been explored by solving the gap equation in different truncation schemes71-74. The CS coefficient acts as an effective vacuum polarization such that when a large number of fermion families are considered, besides the existence of a critical number of families N c that restores chiral symmetry, there exists a critical value for the CS coefficient ϑ c that also induces chiral symmetry restoration. The order of the phase transition in this case changes to order one.

In RQED, the impact of the dimensionless parameter ϑ has been explored by our group (for a study exploring the effects of a modified Chern-Simons term containing a dimensionful parameter see75 and for other approaches involving time reversal symmetry breaking see76). The starting point is the addition of the CS term (136) in the RQED Lagrangian. The effective theory under consideration has the structure36

LRQEDCS=14Fμν2(-)1/2Fμν+ψ-(iγμμ+eγμAμ)ψ+12ζ(αAα)2+-iϑ4εμνρAμνAρ, (137)

where all the quantities carry their usual meaning. The tree-level gauge boson propagator derived from the above Lagrangian is

^μυq=12q1(1+ϑ2)(δμν-qμqνq2)+ζq2qμqνq2-12q2ϑ(1+ϑ2)ϵμνρqρ. (138)

From here it is immediate to notice the effect of the CS coefficient as an effective dielectric constant37 that has the potential to modify the conditions for chiral symmetry breaking. Moreover, for fermions, it is convenient to introduce their right- and left-handed projections, ψ± = Χ±ψ, defined from the chiral matrix τ = [Y3, Y5]/2. In this form, the chiral projectors are defined such that

χ±=121±τ, (139)

which fulfill the following properties

=χ±,χ+χ-=0,χ++χ-=1. (140)

The advantage of working with these projectors is the following: In order to explore the general scenario of chiral symmetry breaking and its consequence, the dynamical generation of fermion masses, we consider the possibility of emergence of an ordinary Dirac mass term meψ¯ψ, which breaks explicitly chiral symmetry, but is even under partity and time-reversal transformations. On top of that, there is another mass term that potentially can be generated, the Haldane mass term moψ¯τψ. Although this term does not break chiral symmetry, it violates parity (and time-reversal) and because the CS is parity violating we expect a relation between them. With these considerations, the most general Lagrangian we consider for our analysis is36,46

LRQEDCS=-14Fμν2(-)1/2Fμν+12ζ(αAα)2+ψ¯(iγμμ+eγμAμ+me+τmo)ψ+-iϑ4εμνρAμνAρ, (141)

These masses however do not correspond to poles in the fermion propagator. The operators χ± project the upper and lower two-component spinors, making explicit that the would-be poles of each fermion species correspond to combination of Dirac and Haldane mass terms. Explicitly, upon acting with the chiral projectors, the matter Lagrangian can be written as

LF=ψ¯+(i+M+)ψ++ψ¯-(i+M-)ψ-, (142)

where M+ = me + m0 and M- = me - m0.

As mentioned earlier, the presence of a CS term generates a mass for the photon in QED3, where ϑ has mass dimension one. Such a mass is known as topological mass (in this case no real topology is involved and the terminology comes from historical reasons). In QED3, radiative corrections may give a contribution to this mass. For the case of QED3 incremented with a CS term, the Coleman-Hill theorem states that only one-loop radiative corrections may be non-vanishing, with the other corrections being identically zero to all orders. We note that the argument to demonstrate this theorem is valid for a theory with massive fermions. Remarkably, it was shown that the Coleman and Hill theorem is also valid for RQED77 in spite of the fact that the CS coefficient is dimensionless. Therefore, the only possible correction to the photon propagator comes from the one-loop diagram Fig. 3a).

In order to address the issue of chiral symmetry and parity breaking, in Ref.35 a robust truncation of the SDE was incorporated. To this end, the possibility of Fermi velocity renormalization by splitting the fermion momentum P = (p0, p) into its temporal p0 and spatial p components such that the propagator reads

S(p0,p)=-F+0(p)γ0p0+F+(p)γp+M+(p)(F+0)2(p)p02+F+2(p)vFp2vF2+M+2(p)Χ+-F-0(p)γ0p0+F-(p)vFγp+M-(p)(F-0)2(p)p02+F-2(p)p2vF2+M-2(p)Χ- (143)

Tha gap equation is truncated including a gauge-boson-fermion vertex dressing model and vacuum polarization effects. The former is based in the central Ball-Chiu vertex, which in this formalism takes the explicit form

ΓνP,K=14Hνσ+P+Hνσ+Kγσ1+γ5+14Hνσ-P+Hνσ-Kγσ1-γ5 (144)

where

H±(P)=1F±0(P)0001vFF±(P)000M±(P)F±(P), (145)

whereas the latter are included modelling a Coulomb-like approximation. Noticing that the most general form of the gauge boson propagator can be written as

ΔμνQ=GLQPμν3+GTPμν1-Pμν3\na+GDPμν6+GEPμν10-Pμν11, (146)

where the involved tensors are

Pμν1=gμν-QμQνQ2, Pμν3=nμnνn2, Pμν6=ϵμναQα, Pμν10=-ϵμαβQαnβnνQ2q2, Pμν11=-ϵναβQαnβnμQ2q2 (147)

with

nμ=gμ0-q0QμQ2. (148)

It turns out that these vacuum polarization effects enter into the fermion gap equation only through the combinations GL and GDE = GD + GE. Then, the gap equation becomes a system of equation for the six unknown functions F±0(P), F±(P) and M±(P), which depend on the value of the coupling constant α and the CS coefficient ϑ. Different models of G L and G DE are proposed that are consistent with a Coulomb-like interaction. The main conclusions that are drawn from this work are

  • The gap equation supports solutions which break parity, associated to the generation of a Haldane mass.

  • Non-trivial parity-preserving solutions are found that, nevertheless break chiral symmetry.

  • The Dirac dynamical mass reduces as θ increases. This confirms that the CS coefficient behaves as an effective dielectric constant, as was proposed in39.

A point-wise evolution of the dynamical mass in terms of the coupling and the CS term was considered by some of us in Ref.36. Assuming that Fermi velocity renormalization has already been carried out and thus no need to split the spatial and temporal components of the momentum, a rainbow-ladder truncation of the gap equation was adopted with the additional simplification of neglecting the wavefunction renormalization effects setting F ±(p) = 1 in Landau gauge. The gap equation is found equivalent to the following uncoupled system of non-linear integral equations for the left and right handed mass functions

M±p=2παd3k(2π)32M±(k)k2+M±2(k)1q(1+ϑ2)1q2ϑ1+ϑ2kqk2+M±2(k). (149)

Inspired by the reasoning of the previous section, the gap equation can be linearized and reduced to the following differential equation

p2M±''p+2pM±'p+2απ1+δ2M±p=2aϑπ1+ϑ25p3-6p2κ+Λ39p2-p9 (150)

with the boundary conditions

p2M'±p|pκ=0, pM'±p+M±p|p=Λ=0, (151)

where κ is the IR regulator introduced as a scale that quantifies the amount of mass being dynamically generated. The general solution to Eq. (150) is

M±p=p-121-ααc-12×c±2p1-ααc+c±1+fϑ,p, (152)

where

αc=π8(1+ϑ2) (153)

and

fϑ,p=πϑ2+12aϑϑ2+1κα+3πϑ2+1×Λ3+14p3-54p4α+πϑ2+1×118κp2+α+πϑ2+1α+3πϑ2+1 (154)

We notice that when ϑ = 0, ƒ(ϑ, p) = 0 and we obtain precisely the solution discussed in (89) which demands α > αc = π/8. For ϑ ≠ 0, we observe the role of the CS coefficient ϑ as an effective dielectric constant. The function ƒ(ϑ, p) describes deviations from the Miransky scaling by virtue of the CS coefficient.

The numerical solution to the non-linear Eqs. (149) are shown in Figs. 18 and 19 for ϑ = ηϑc with varying 0 ≤ η ≤ 1 but keeping α = 1.07 αc fixed. These numbers are consistent with a critical value of ϑc = 6.6 × 10-14. For small values of p, we observe that both masses are finite and approximately flat as p → 0. For larger momentum, the height of the mass functions decreases. This happens for different values of ϑ provided α > αc. Furthermore, Figs. 20 and 21 show that the dynamical masses depend on the magnitude of the CS coefficient in opposite forms for right- and left-handed fermions. While M_(0) increases M+(0) decreases as ϑ grows bigger. At the critical value ϑc, M+(0) flips sign abruptly, signaling a first order phase transition.

Figure 18 Mass function M_(p) plotted as a function of momentum. The different curves correspond to different values of the CS parameter ϑ = ηϑc with 0 ≤ η ≤ 1 and fixed value of the coupling α = 1.07 αc

Figure 19 Mass function M+(p) plotted as a function of momentum. The different curves correspond to different values of the CS parameter ϑ = ηϑc with 0 ≤ η ≤ 1 and fixed value of the coupling α = 1.07 αc

Figure 20 The absolute value of the height of the right-handed fermion mass function M+(0) plotted as a function of the CS parameter ϑ. Above ϑc the mass M+(0) becomes negative and here we flipped its sign for a better visualization. The different curves correspond to different values of the coupling α = xαc

Figure 21 Height of the left-handed fermion mass function M_(0) plotted as a function of the CS parameter ϑ. The different curves correspond to different values of the coupling α = xαc

Notice that the heights M +(0) of the mass functions parametrize how much mass has been generated. These hights are depicted in Figs. 20 and 21 for a constant value of α = xαc. Above ϑc, M+(0) undergoes a discontinuity and jumps assuming negative values, as shown in Fig. 22, where we have flipped the sign after the discontinuity to lead the eye. The blue dots stand for the absolute value of M+ above ϑc, where we observe that the two masses become the mirror image of one another. We notice that since the curves are re-scaled, they coincide for any value of x > 1.

Figure 22 M+(0) re-scaled by their value at ϑ = 0. The solid line corresponds to M+(0) and dashed line M_(0). The red dots are the absolute value |M+(0)| above the critical ϑ. 

The Dirac and Haldane masses that appear in the Lagrangian (141) can be obtained summing up or subtracting M+(0), respectively. From Figs. 23 and 24, we notice that both these masses are non-vanishing for a large range of values of ϑ, clearly showing that CSB can occur within this model. For small ϑ, the Dirac mass me, associated to chiral symmetry breaking, presents a plateau, leaping at a critical ϑ to a value several orders of magnitude smaller. After that, me increases again and finally drops to zero. We interpret the discontinuity as an attempt to restore chiral symmetry, associated to the first piece in Eq. (149). That would be the critical point if ƒ(ϑ, p) as absent. However, as one moves towards the right of the discontinuity, the CS contribution ƒ(ϑ, p) becomes dominant and me starts to be significant again. For even larger ϑ the Dirac mass drops and the real critical ϑ, for which chiral symmetry is restored, is around ϑ = 1. In turn, parity breaking is encoded in m0. As the CS parameter becomes dominant, the Haldane mass increases, never restoring the symmetry.

As in the finite temperature case, the gauge dependence of these findings represent an open avenue to explore73,74. Furthermore, feedback effects of the dynamical mass and the CS coefficient ϑ need to be incorporated by solving the SDE for the fermion and photon propagators simultaneously.

Figure 23 Dynamical mass term, even under parity transformation, as a function of the CS parameter ϑ

Figure 24 Dynamical mass term, odd under parity transformation, as a function of the CS parameter ϑ

4.3. Anisotropy effects on chiral symmetry breaking

Straintronics has emerged as the field in which the electrical properties of graphene and other materials are manipulated through mechanical deformations (see, for instance, Ref.78). Effects of anisotropy (strain) in the gap equation have been recently considered in Ref.38. By writing the fermion propagator of RQED with strain as

S(p)=-γμMμνpν-1, (155)

with

Mμν=1000v1000v2, (156)

the gap equation is considered as a function of the anisotropy parameter η = ν1/ ν2. Including one-loop vacuum polarization effects and a suitable generalization of the Ball-Chiu vertex, the authors of Ref.38 solved the gap equation and searched for the effect of η on the critical coupling. They observe that the net effect is to increase the value of αc for chiral symmetry breaking.

4.4. Chiral symmetry breaking in curved space

It is well known that for 2D materials, electric and optical properties are modified by impurities and deformation. The latter can be introduced in a framework in which these mechanical deformations can be incorporated in terms of curved geometry of the underlying space-time in the equations of motion. In this regard, the dynamics of gauge bosons and electrons in mixed space-time dimensions has been considered in a curved-space background for the latter (and flat extra dimensions for the former)39. The starting point is the action

S=ddγx-q&#091;-14FμνFμν-12ξ(μAμ)2&#093;\na+ddex-Hψ¯iγμ(x)(μ+Ωμ+ieAμ)ψ, (157)

where Fμν=μAν-νAμ=μAν-νAμ, due to the cancellation of connection terms in the derivatives μ, γμ(x)=eνμγν, (Ωμ)βα=1/2ωμab(Jab)βα, where (Jab)βα represent the Lorentz generators in spinor space and ωμba=ebυ-δλνμ+Γλμνeλa stands for the spin connection, related to the Christoffel connection through the metricity condition namely μeνa=μeνa-Γμνλeaλ+(ωμ)abeνn=0. In all the above, the vielbeineλa verifies ηabeμaeνb=gμν. The induced metric on the brane (boundary space) is denoted by Hαβ and gμν is the bulk space metric. For applications in graphene, it is considered that

gμνdxμdxν=dt2-dz2-hijdxidxj, (158)

with i, j = 1,2 and thus

ddγx-g=dt dz dx1 dx2 h. (159)

Considering the representation

e~μa(x)={eμa(x)δ(xdγ-de),a,μ=μe0a,μ=de,dγ-1, (160)

where in the case of graphene xdγ-de corresponds to the third spatial dimension. In fact, assuming that all extra dimensions are flat, the action of this theory is very similar to the corresponding action of QED is curved space39, namely,

S=ddγx-g-14FϑνFϑν-12ξμAμ2ψ-γ-μxμ+Ωμ+ieAμψ (161)

where γ¯μ(x)=eμa(x)γa. In the local momentum representation, the propagators of the theory have the form

S0p=pp2-Me2,Δμν0q=1q2-Mγ2gμν+ξ-1qμqνq2-Mγ2, (162)

where

Me2=R(x')12,Mγ2=-R(x')6, (163)

where R(x’) is the Ricci scalar. For graphene, Mγ2=0 and Me is a constant.

The one-loop renormalization of the fermion propagator has been discused in Ref39. All divergent terms are the same as in curved-space ordinary QED. Curvature effects are only seen in the finite parts of the self-energies. As for the polarization tensor, such a term is found finite, hence making the beta function null at the one-loop level. Furthermore, the curvature effects in the case of graphene show up as an effective chemical potential of the form μ=Me2vF4, where vF is the effective Fermi velocity.

5. Landau-Khalatnokov-Fradkin transformations

Gauge symmetry lays at the very foundation of the description of fundamental interactions. It can manifest in different ways at many levels. For the Green functions of QED, different sets of relations among them can be derived from the fundamental symmetry of the theory, as for example Ward-,51 Ward-Green-Takahashi-51-53 and Transverse Ward-identities54,79-82 which relate (n + 1)-point to n-point functions in constructions resembling divergence and curl of currents. In a different setting, the so-called Landau-Khalatnikov-Fradkin transformations (LKFT)40,41 which have been analyzed in different versions of QED and extended to non-Abelian gauge theories like QCD, have been derived from different arguments in the past decades83-87. These transformations are non-perturbative in nature and hence have the nature to address the issue of gauge invariance in perturbative and non-perturbative studies of field theories.

A widely used example of LKFT has been in the case of the fermion propagator. It has been used to establish the multiplicative renormalizability of spinor and scalar QED in different space-time dimensions88-91. The main lesson is that this transformation fixes some of the coefficients of the perturbative expansion of the fermion propagator at all orders. In the case of RQED, given that WTI are satisfied, it turns interesting that the LKFT serves as a tool to establish general features of multi-loop calculations as well as non-perturbative calculations in connection with chiral symmetry breaking. In what follows of this section, we review the LKFT in RQED.

5.1. LKFT for the fermion propagator in equal dimensional QED

We start from the most general form of the Dirac fermion propagator in arbitrary space-time dimensions and in momentum space and relate it to its coordinate space representation

S(x;ξ)=xX(x;ξ)+Y(x;ξ) (164)

through a Fourier transformation. We have additionally added the gauge parameter label ξ as we are interested in the form of these propagators for different gauges.

For theories of fermions and bosons in the same space-time dimensions, the momentum space free gauge boson propagator Δμν(0)(q) has the general form (in analogy with (92))

Δμν(0)(q)=ΔμνT(q)+ξqμqν(q2)2, (165)

such that the longitudinal portion of the propagator, inversely proportional to q 4, describes how this Green function changes in different gauges, being proportional to the gauge fixing parameter. The LKFT emerges precisely from this part of the propagator40,41,83-87. These transformations are more clearly written in coordinate space. For the fermion propagator, the LKFT states that the fermion propagator in an arbitrary covariant gauge S(x; ξ) is related to its form in Landau gauge S(x; 0) through the transformation

S(x;ξ)=S(x;0)e-iΔd(0)-Δd(x), (166)

where

Δd(x)=-iξe2μ4-dddq(2π)de-iqxq4, (167)

e denoting the electric charge, and μ is an energy scale that renders e dimensionless in d = 4, but yields a dimensionful coupling e2 in QED3. Performing the momentum integration, Δd (x) is explicitly given by43

Δd(x)=-iξα4πd-22Γd-42(μx)4-d, (168)

where α = e2/(4 π) is the coupling constant, and Γ(z) is the Euler Gamma function.

To taste a feeling of the structure of LKFT, let us take F(p; 0) = 1 and M(p; 0) = 0, namely, the massless fermion propagator in Landau gauge as a seed of LKFT.

For d = 3, the transformation reveals that the wavefunction renormalization in an arbitrary covariant gauge is (see, for instance, Ref.42,43)

F(p;ξ)=1-α¯ξ2parctan2pα¯ξ, (169)

with α¯=e2/(4π) which has mass dimension one. We see directly that a weak coupling expansion fixes all terms of the form (α¯ξ)j at any given order in perturbation theory. This is a major asset of the LKFT. Below we discuss in detail the progress achieved so far regarding the LKFT in RQED.

5.2. The LKF transformation for RQED

Let us consider the free gauge boson propagator over bulk dimensions in RQED. It is of the same form as in Eq. eqn:phototond, but when considered reduced to the de -dimensional brane, it becomes27,29

Δμeνeq=Dq2gμeνe-qμeqνeq2+ξ~D(q2)qμeqνeq2, (170)

with,

D(q2)=i(4π)εeΓ(1-εe)(-q2)1-εe, (171)

where εe = ( - de)/2 and ξ̃=(1-εe)ξ. From here we see that the longitudinal part of the propagator changes the form of the LKFT for the fermion propagator from the equal dimensional theory. Considering that the propagator changes from gauge to gauge according to

Sde(x;ξ)=Sde(x;0)e-i&#091;Δ̃de(0;εe)-Δ̃de(x;εe)&#093;, (172)

we define

~dex; ε2=ifεeξe2μ4-dγddeq(2π)dee-iqxq4-2εe, (173)

=-if(εe)ξe2Γ(de-a2)2aπde/2Γ(a2)(μx)a-de (174)

With

f(εe)=Γ(1-εe)(1-εe)(4π)εe,a=4-2εe. (175)

This is the general form of LKFT for the fermion propagator in RQEDdγ - de first presented in REF.45. Note that when dγ - de = d, Eq. (173) matches the usual form of the LKFT Eq. (166) in any dimension d. For graphene and other 2D materials,

Δ̃3x,12=-iξe216π2Γ1-2ϵ2(μx)2ϵ-1, (176)

where the limit ∈→1/2 is understood. Expanding Eq. (176) around ∈ = 1/2, defining δ = ∈ = -1/2, we get

Δ̃3x,12=iξe216π21δ+γE+2ln(μx)+O(δ). (177)

Since the transformation function, Eq. (177), diverges at x = 0, we introduce a cutoff xmin and consider

-i~3xmin,12-~3x,12=lnxxmin-2ν, (178)

where we have defined v = ξα / (4 π), and the dimensionless coupling constant α=e2/(4π).

With Eq. (178), we can obtain the non-perturbative structure of the fermion propagator in RQED starting with the tree-level massless case in Landau gauge. In coordinate space, we have

Xx;ξ= Xx;0e-i~3xmin,12-~3x,12=-xmin2ν4πx-2ν-3. (179)

whereas in momentum space,

F(p,ξ)=π2Γ(1-ν)Γ(32+ν)p2Λ2ν. (180)

This is the non-perturbative form of the fermion propagator in RQED in any covariant gauge ξ. The power-law behavior of this two-point function is consistent with the multiplicative-renormalizable character of the theory. We notice that a similar behavior is found for the propagator in ordinary QED.

Next, expanding this result in powers of α, we get

F(p,ξ)=1+ξα4πF1+ξα4π2F2+O(α3), (181)

with

F1=ln(p2Λ2)-γE-ψ(32)=ln(p2Λ2)+2γE+ln(4)-2, (182)

F2=12lnp2Λ2-γE-ψ322-2ζ2+4=12lnp2Λ2+2γE+ln4-22-2ζ2+4, (183)

where ψ(z) is the digamma function, ζ(s) is the Riemann Zeta function, and we have made use of the identity ψ(3/2) = γE - ln(4) + 2.

We make a comparison against the exact one-loop calculation of the wavefunction renormalization28,29,

Fp;ξ=1+α4π49-1-3ξ3L¯+α4π2(1-3ξ)218L¯2-2ζ2+4+43ξ+7L¯+48ζ227-8ζ2L¯+2-ln4-28027, (184)

where

L¯=ln-p2μ2+ln(4)-2, (185)

and μ is the renormalization mass scale. This correction has a contribution that is proportional to the gauge parameter and one that remains finite in Landau gauge. Thus, we can only compare against the former from the LKFT result. F1 defined in is equivalent to L¯ in provided we identify

lnp2Λ2+2γEln-p2μ2. (186)

At O(α2), the propagator has terms independent, linear and quadratic ξ. From Eqs. (181) and (183), we observe the LKFT only gives terms proportional ξ2 at order α2. This means that we can only compare our α2 result, F 2, defined by with the coefficient of (α ξ / (4 π))2 in the perturbative result, Eq. (184). We observe that

F2918L¯2-2ζ(2)+4, (187)

with the identification Eq. (186), as expected. Thus, we have shown that there is full consistency between our LKFT result, Eq. (181), and the perturbative result, Eq. (184), up to order α2. We hence predict the form of all the coefficients of the form (αξ)j in the all order perturbative expansion from our LKFT result Eq. (180).

5.3. Extensions of the LKF transformation

It is interesting to notice that the LKF transformation can be represented directly in momentum space. In Refs.92-95 it is discussed the form in which the wavefunction renormalization is related in two gauges labeled by the parameters ξ and η, respectively. Assuming that in gauge η the multiloop expansion of the wavefunction renormalization is

F(p,η)=m=0am(η)αmμ̃2p2mϵ, (188)

with ∈ = 2 - d/2 and ũ is a renormalization scale, the LKFT implies that

F(p,ξ)=m=0am(ξ)αmμ̃2p2mϵ, (189)

where the new coefficients

amξ=am(η)Γ2-m+1ϵΓ1+mϵ×l=0Γ1+m+lϵΓl1-ϵl!Γ2-m-l+1ϵ×η-ξα-ϵlμ~2p2lϵ, (190)

With the appropriate choice of the renormalization scale ũ, the multiloop expansion of the fermion propagator can be obtained in an easier form.

Moreover, taking advantage of the above mentioned multiplicative renormalizability of the theory, by factorizing

F(p;ξ)=Zψ(α,ξ)Fr(p;ξ), (191)

one can see from LKFT that

logZψ(α,ξ)=logZψ(α,η)-(ξ-η)αϵ. (192)

Thus, for the fermion anomalous dimension

γψ(α,ξ)=-β(α)logZψ(α,ξ)α-ξαlogZψ(α,ξ)ξ, (193)

where β(α) is the beta function and γ(α) the gauge boson anomalous dimension, given explicitly as

β(α)=-2ϵα+l=0βlαl+1,γ(α)=-l=0βlαl, (194)

where β1 = 0 for graphene. From here its is evident that the fermion anomalous dimension in two gauges, namely,

γψ(α,ξ)=γψ(α,η)-2α(ξ-η), (195)

i.e., the gauge dependence is encoded in the one-loop correction alone, whereas higher loops are gauge invariant.

The two-loop correction to the mass anomalous dimension γm has been explored very recently in68. In this work, the authors consider a theory in which the gauge boson lives in 4 dimensions, but the fermion fields are allowed to live in a general dimension de. From the multiplicative renormalizability of the fermion mass,

m=Zmmr, (196)

where the renormalization constant Zm is obtained from the relation

1+ΣS1-ΣVZm=const. (197)

where Σs and Σv are the scalar and vector projections of the fermion self-energy. Furthermore, the mass anomalous dimension runs with the energy scale μ as

γm=dlogZmdlogμ. (198)

Using the MS¯-regularization scheme, it is observed that in the particular case of RQED,

γm=323αr4π-64αr4π2(Nζ2-827)+Oαr4π3, (199)

where N is the number of fermion families and ζ2 is the Euler function. This mass anomalous function is completely gauge invariant, and thus the critical coupling and/or critical number of fermion flavor derived from γm are also gauge invariant.

6. Outlook

In this review article we have revisited several aspects of a mixed-dimensional theory of electron and gauge boson interactions aiming to describe the new family of 2D Dirac material like graphene and its cousins. Quantum field theories in mixed dimensions have been discussed theoretically for over 4 decades in literature, but the isolation of graphene has indeed boosted the interest of these highly non-local yet well behaved theories. The fact that RQED, the theory we have specialized in, is a fixed point in the renormalization group allows to extend general conclusions from the behavior of the Green’s functions in RQED. From the historical point of view, we have revisited the formulation of the theory and the most fundamental aspects that render the theory suitable for phenomenological description of the properties of a variety of Dirac materials. The fact that the theory allows a standard 1/r fall-off for static Coulomb interactions to the observation of the scale invariance of the theory makes it already interesting to explore. Fundamental aspects of causality and unitarity have been revisited along with perturbation theory studies of the propagators and vertices.

An interesting aspect of this formulation is the possibility to carry out non-perturbative calculations in particular those connected with the chiral symmetry of the theory which might be broken when the coupling is strong enough. Schwinger-Dyson equations and renormaization group analysis exploit the structure of the gap equation as compared to other formulations of QED3 including vacuum polarization effects allow to extract critical properties of the chiral transition in RQED. Still a number of improvements and refining of the truncations to the infinite tower of SDE in particular in connection with the gauge invariance of the predictions still awaits for a formal answer.

The influence of external agents has also been considered in this theory. Finite temperature effects indicate a relation of the transition temperature competing against the strong coupling which by virtue of the heat bath diminishes, hence restoring chiral symmetry. Findings of different groups are still controversial, hence pointing toward the need of a careful revision of the working hypothesis. Moreover, the influence of a chemical potential is still required to have a more realistic description of the physical system and eventual test in the laboratory. In this sense, an effort was made in96 to address the role of chemical potential in a configuration with spin-orbit coupling. The phase diagram in the plane of temperature and chemical potential of RQED, similar to that of QCD, is a pending assignment for the interested community.

The possibility of parity and time reversal symmetry breaking has been addressed by different groups by including the CS term in the fundamental Lagrangian. The existence of parity and time reversal breaking solution is remarkable in this theory. The effect of the CS coefficient theta as an effective dielectric constant is seen to impact the chiral symmetry breaking pattern inasmuch as the effective screening might restore chiral symmetry but still break the discrete symmetries in the theory.

An interesting feature is the modeling of elastic deformation and anisotropy in relation to the dynamics of fermions in curved space. The effects of curvature are seen as density effects through a definition of a position dependent chemical potential.

LKF tranformations are a key ingredient from gauge invariance that have a fundamental role in perturbative and non-perturbative studies in RQED. In the former case, the closed connection between RQED and QED3 in the large N limit serves as a benchmark for perturbative calculations in improving the calculation of vacuum polarization effects.

Many features are still pending to address in this theory. Some of which include the influence of magnetic and electric field in different configurations and potential connections with other phenomenology are important features that perhaps the community will address in the future.

Acknowledgements

We benefited from discussions with L. A. Rangel, S. Henrnández-Ortiz, K. Raya, J. C. Rojas, J. Báez-Cuevas. JACO and AR acknowledge financial support from Consejo Nacional de Ciencia y Tecnología (México) under grant FORDECyT-PRONACES/61533/2020. AJM acknowledge financial support from FAPESP under grant 2016/12705-7.

References

1. V. A. Miransky and I. A. Shovkovy, Quantum field theory in a magnetic field: From quantum chromodynamics to graphene and Dirac semimetals, Phys. Rept. 576, 1 (2015) https://doi.org/10.1016/j.physrep.2015.02.003 [ Links ]

2. P. R. Wallace, The Band Theory of Graphite Phys. Rev. 71 (9), (1947), 622-634. https://doi.org/10.1103/PhysRev.71.622 [ Links ]

3. G. W. Semenoff, Condensed-Matter Simulation of a Three-Dimensional Anomaly Phys. Rev. Lett. 53, (1984) 2449. https://doi.org/10.1103/PhysRevLett.53.2449 [ Links ]

4. F. D. M. Haldane, Model for a Quantum Hall Effect without Landau Levels: Condensed-Matter Realization of the Parity Anomaly. Phys. Rev. Lett. 61, (1988) 2015-2018. https://doi.org/10.1103/PhysRevLett.61.2015 [ Links ]

5. N Dorey, and N. E. Mavromatos QED3 and two-dimensional superconductivity without parity violation Nucl Phys B 386 (1992) 614. https://doi.org/10.1016/0550-3213(92)90632-L [ Links ]

6. K. Farakos and N. E. Mavromatos Dynamical gauge symmetry breaking and superconductivity in three-dimensional systems Mod Phys Lett A 13 (1998) 1019. https://doi.org/10.1142/S0217732398001091 [ Links ]

7. M. Franz and Z. Tesanovic Algebraic Fermi Liquid from Phase Fluctuations: “Topological” Fermions, Vortex “Berryons,” and QED 3 Theory of Cuprate Superconductors Phys Rev Lett 87 (2001) 257003. https://doi.org/10.1103/PhysRevLett.87.257003. [ Links ]

8. I. F. Herbut, QED 3 theory of underdoped high-temperature superconductors Phys Rev B 66 (2002) 094504. https://doi.org/10.1103/PhysRevB.66.094504 [ Links ]

9. M. Franz, Z. Tesanovic, O Vafek QED 3 theory of pairing pseudogap in cuprates: From d-wave superconductor to antiferromagnet via an algebraic Fermi liquid Phys Rev B 66 (2002) 054535. https://doi.org/10.1103/PhysRevB.66.054535 [ Links ]

10. K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V Grigorieva and A. A. Firsov, Electric field effect in atomically thin carbon films, Science 306, (2004) 666. https://doi.org/10.1126/science.1102896 [ Links ]

11. Zhang Y, Tan Y-W, Horst L S and Kim P, Experimental observation of the quantum Hall effect and Berry’s phase in graphene Nature 438 (2005), 201. https://doi.org/10.1038/nature04235 [ Links ]

12. K. S. Novoselov, et al., Two-dimensional gas of massless Dirac fermions in graphene Nature 438, (2005) 197-200. https://doi.org/10.1038/nature04233 [ Links ]

13. A. Geim, and K. Novoselov, The rise of graphene, Nature Mater 6, (2007) 183. https://doi.org/10.1038/nmat1849 [ Links ]

14. C. J. Burden, J. Praschifka, C. D. Roberts, Photon polarization tensor and gauge dependence in three-dimensional quantum electrodynamics, Phys. Rev. D 46(6), (1992) 2695-2702. https://doi.org/10.1103/physrevd.46.2695 [ Links ]

15. E. Marino, Quantum Field Theory Approach to Condensed Matter Physics, (Cambridge University Press, Cambridge, UK, 2017). [ Links ]

16. E. Marino, Quantum electrodynamics of particles on a plane and the Chern-Simons theory, Nucl. Phys. B 408, (1993) 551-564. https://doi.org/10.1016/0550-3213(93)90379-4 [ Links ]

17. J. Gonzalez, F. Guinea and M. A. H. Vozmediano, Non-Fermi liquid behavior of electrons in the half filled honeycomb lattice (A Renormalization group approach), Nucl. Phys. B 424 (1994), 595-618. https://doi.org/10.1016/0550-3213(94)90410-3 [ Links ]

18. D. Tong, The Quantum Hall Effect, Lecture Notes, https://www.damtp.cam.ac.uk/user/tong/qhe/qhe.pdfLinks ]

19. R. L. P. G. do Amaral and E. C. Marino, Canonical quantization of theories containing fractional powers of the d’Alembertian operator, J. Phys. A 25 (1992), 5183-5200 https://doi.org/10.1088/0305-4470/25/19/026 [ Links ]

20. E. C. Marino, L. O. Nascimento, V. Alves and C. M. Smith, Unitarity of theories containing fractional powers of the d’Alembertian operator ,” Phys. Rev. D 90 (2014) no.10, 105003 https://doi.org/10.1103/PhysRevD.90.105003 [ Links ]

21. J. Gonzalez, F. Guinea and M. A. H. Vozmediano, Unconventional Quasiparticle Lifetime in Graphite, Phys. Rev. Lett. 77 (1996), 3589-3592 https://doi.org/10.1103/PhysRevLett.77.3589 [ Links ]

22. J. Gonzalez, F. Guinea and M. A. H. Vozmediano, Marginal-Fermi-liquid behavior from two-dimensional Coulomb interaction, Phys. Rev. B 59 (1999), 2474-2477 http://doi.org/10.1103/PhysRevB.59.R2474 [ Links ]

23. J. Gonzalez, F. Guinea and M. A. H. Vozmediano, Electron-electron interactions in graphene sheets Phys. Rev. B 63 (2001) 134421 http://doi.org/10.1103/PhysRevB.63.134421 [ Links ]

24. M. A. H. Vozmediano, Renormalization group aspects of graphene, Phil. Trans. Roy. Soc. Lond. A 369 (2011), 2625-2642 http://doi.org/10.1098/rsta.2010.0383 [ Links ]

25. L. Fernández, V. Alves, L. O. Nascimento, F. Peña, M. Gomes and E. C. Marino, Phys. Rev. D 102 (2020) no.1, 016020 doi:10.1103/PhysRevD.102.016020 ENT#091;arXiv:2002.10027 ENT#091;hep-thENT#093;ENT#093;. [ Links ]

26. E. V. Gorbar, V. P. Gusynin and V. A. Miransky, Dynamical chiral symmetry breaking on a brane in reduced QED, Phys. Rev. D 64, (2001) 105028 https://doi.org/10.1103/PhysRevD.64.105028 [ Links ]

27. S. Teber, Electromagnetic current correlations in reduced quantum electrodynamics, Phys. Rev. D 86 (2012) 025005, https://doi.org/10.1103/PhysRevD.86.025005 [ Links ]

28. A. V. Kotikov and S. Teber, Two-loop fermion self-energy in reduced quantum electrodynamics and application to the ultrarelativistic limit of graphene, Phys. Rev. D 89, no. 6, (2014) 065038 https://doi.org/10.1103/PhysRevD.89.065038 [ Links ]

29. S. Teber, Two-loop fermion self-energy and propagator in reduced QED 3,2 , Phys. Rev. D 89, no. 6, (2014) 067702, https://doi.org/10.1103/PhysRevD.89.067702 [ Links ]

30. S. Teber and A. Kotikov, Field theoretic renormalization study of reduced quantum electrodynamics and applications to the ultrarelativistic limit of Dirac liquids, Phys. Rev. D 97, no.7, (2018) 074004. https://doi.org/10.1103/PhysRevD.97.074004 [ Links ]

31. V. Alves, W. S. Elias, L. O. Nascimento, V. Juričić and F. Peña, “Chiral symmetry breaking in the pseudo-quantum electrodynamics ,” Phys. Rev. D 87, no.12, (2013) 125002. https://doi.org/10.1103/PhysRevD.87.125002 [ Links ]

32. A. Kotikov and S. Teber, Critical behaviour of reduced QED 4,3 and dynamical fermion gap generation in graphene, Phys. Rev. D 94, no.11, (2016) 114010. https://doi.org/10.1103/PhysRevD.94.114010 [ Links ]

33. L. O. Nascimento, V. Alves, F. Peña, C. M. Smith and E. Marino, Chiral-symmetry breaking in pseudoquantum electrodynamics at finite temperature, Phys. Rev. D 92, (2015) 025018. https://doi.org/10.1103/PhysRevD.92.025018 [ Links ]

34. J. Báez-Cuevas, A. Raya and J. C. Rojas, Chiral symmetry restoration in reduced QED at finite temperature in the supercritical coupling regime, Phys. Rev. D 102 (2020), 056020. https://doi.org/10.1103/PhysRevD.102.056020 [ Links ]

35. M. E. Carrington, Effect of a Chern-Simons term on dynamical gap generation in graphene, Phys. Rev. B 99 (2019) 115432. https://doi.org/10.1103/PhysRevB.99.115432 [ Links ]

36. J. A. C. Olivares, L. Albino, A. J. Mizher and A. Raya, Influence of a Chern-Simons term in the dynamical fermion masses in reduced or pseudo QED, Phys. Rev. D 102 (2020) no.9, 096023 https://doi.org/10.1103/PhysRevD.102.096023 ENT#091;arXiv:2009.08484 ENT#091;hep-phENT#093;ENT#093;. [ Links ]

37. G. C. Magalhães, V. Alves, E. C. Marino and L. O. Nascimento, Pseudo Quantum Electrodynamics and Chern-Simons theory Coupled to Two-dimensional Electrons, Phys. Rev. D 101 (2020) no.11, 116005 https://doi.org/10.1103/PhysRevD.101.116005. [ Links ]

38. M. E. Carrington, A. R. Frey and B. A. Meggison, Effect of anisotropy on phase transitions in graphene, Phys. Rev. B 102 (2020) no.12, 125427. https://doi.org/10.1103/PhysRevB.102.125427 [ Links ]

39. P. I. C. Caneda and G. Menezes, Reduced quantum electrodynamics in curved space, Phys. Rev. D 103, (2021) 065010. https://doi.org/10.1103/PhysRevD.103.065010 [ Links ]

40. L. D. Landau and I. M. Khalatnikov, The gauge transformation of the Green function for charged particles, Sov. Phys. JETP 2 (1956) 69 ENT#091;Zh. Eksp. Teor. Fiz. 29 (1955) 89ENT#093;. [ Links ]

41. E. S. Fradkin, Concerning some general relations of quantum electrodynamics, Zh. Eksp. Teor. Fiz. 29, 258 (1955) ENT#091;Sov. Phys. JETP 2, 361 (1956)ENT#093;. [ Links ]

42. A. Bashir, Non-perturbative Fermion Propagator for the Massless Quenched QED3 Phys. Lett. B 491, (2000) 280. https://doi.org/10.1016/S0370-2693(00)01043-1 [ Links ]

43. A. Bashir and A. Raya, Landau-Khalatnikov-Fradkin transformations and the fermion propagator in quantum electrodynamics, Phys. Rev. D 66, (2002) 105005. https://doi.org/10.1103/PhysRevD.66.105005 [ Links ]

44. A. Bashir and A. Raya, Fermion propagator in quenched QED3 in the light of the Landau-Khalatnikov-Fradkin transformation, Nucl. Phys. Proc. Suppl. 141, (2005) 259 https://doi.org/10.1016/j.nuclphysbps.2004.12.039 [ Links ]

45. A. Ahmad, J. J. Cobos-Martínez, Y. Concha-Sánchez, and A. Raya, Landau-Khalatnikov-Fradkin transformations in reduced quantum electrodynamics, Phys. Rev. D 93, (2016) 094035. https://doi.org/10.1103/PhysRevD.93.094035 [ Links ]

46. D. Dudal, A. J. Mizher and P. Pais, Exact quantum scale invariance of three-dimensional reduced QED theories, Phys. Rev. D 99 (2019) no.4, 045017- https://doi.org/PhysRevD.99.045017 [ Links ]

47. C. G. Bollini and J. J. Giambiagi, Arbitrary powers of D’Alembertians and the Huygens’ principle, J. Math. Phys. 34 (1993), 610-621 https://doi.org/10.1063/1.530263 [ Links ]

48. A. Coste and M. Luscher, Parity Anomaly and Fermion Boson Transmutation in Three-dimensional Lattice QED, Nucl. Phys. B 323 (1989), 631-659 https://doi.org/10.1016/0550-3213(89)90127-2 [ Links ]

49. C. D. Roberts and A. G. Williams Dyson-Schwinger Equations and the Application to Hadronic Physics Prog. Part. Nucl. Phys. 33 (1994) 477. https://doi.org/10.1016/0146-6410(94)90049-3 [ Links ]

50. M. R. Pennington Swimming with Quarks, J. Phys. Conf. Ser. 18 (2005) 1. https://doi.org/10.1088/1742-6596/18/1/001 [ Links ]

51. J. C. Ward, An Identity in Quantum Electrodynamics, Phys. Rev. 78 (1950) 182 https://doi.org/10.1103/PhysRev.78.182 [ Links ]

52. H. S. Green, A Pre-renormalized quantum electrodynamics, Proc. Phys. Soc. (London) A66, (1953) 873. https://doi.org/10.1088/0370-1298/66/10/303 [ Links ]

53. Y. Takahashi, On the generalized ward identity, Nuovo Cimento 6, (1957) 371. https://doi.org/10.1007/BF02832514 [ Links ]

54. Y. Takahashi, Canonical quantization and generalized Ward relations: Foundation of nonperturbative approach, Print-85-0421 (Alberta), (1985). [ Links ]

55. V. Miransky and K. Yamawaki, Conformal phase transition in gauge theories Phys. Rev. D 55, (1997) 5051. https://doi.org/10.1103/PhysRevD.55.5051. [ Links ]

56. A. Bashir, A. Huet, and A. Raya Gauge dependence of mass and condensate in chirally asymmetric phase of quenched three-dimensional QED Phys. Rev. D 66, (2002) 025029. https://doi.org/10.1103/PhysRevD.66.025029. [ Links ]

57. K. Raya, A. Bashir, S. Hernández-Ortiz, A. Raya, and C. D. Roberts, Multiple solutions for the fermion mass function in QED3, Phys. Rev. D 88, (2013) 096003. https://doi.org/10.1103/PhysRevD.88.096003 [ Links ]

58. T. Appelquist, M. J. Bowick, E. Cohler, and L. C. R. Wijewardhana Chiral-symmetry breaking in 2+1 dimensions Phys. Rev. Lett. 55 (1985) 1715. https://doi.org/10.1103/PhysRevLett.55.1715. [ Links ]

59. T. W. Appelquist, M. Bowick, D. Karabali, and L. C. R. Wijewardhana Spontaneous chiral-symmetry breaking in three-dimensional QED Phys. Rev. D 33 (1986) 3704. https://doi.org/10.1103/PhysRevD.33.3704. [ Links ]

60. T. Appelquist, D. Nash, and L. C. R. Wijewardhana, Critical Behavior in (2+1)-Dimensional QED Phys. Rev. Lett. 60 (1988) 2575. https://doi.org/10.1103/PhysRevLett.60.2575 [ Links ]

61. A. Bashir, A. Raya, S. Sánchez-Madrigal and C.D. Roberts, Gauge Invariance of a Critical Number of Flavours in QED3. Few-Body Syst 46 (2009) 229. https://doi.org/10.1007/s00601-009-0069-9. [ Links ]

62. A. V. Kotikov and S. Teber, Critical Behavior of (2+1)-Dimensional QED: 1/N Expansion, Particles 2020, 3, 345-354. https://doi.org/10.3390/particles3020026 [ Links ]

63. A. V. Kotikov and Sofian Teber, Critical behaviour of (2+1)-dimensional QED: 1/N-corrections, EPJ Web of Conferences 138 (2017), 06005. https://doi.org/10.1051/epjconf/201713806005 [ Links ]

64. A. V. Kotikov and S. Teber, Critical behavior of (2+1)-dimensional QED: 1/N f corrections in an arbitrary nonlocal gauge, Phys. Rev. D 94, (2016) 114011. https://doi.org/10.1103/PhysRevD.94.114011 [ Links ]

65. A. V. Kotikov, V.I. Shilin and S. Teber, Critical behavior of (2+1)-dimensional QED: 1/N f corrections in the Landau gauge, Phys.Rev.D 94 (2016) 5, 056009. https://doi.org/10.1103/PhysRevD.94.056009 [ Links ]

66. J. A. Casimiro-Olivares Generación dinámica de Masas en R(P)QED M. Sc. Thesis (In Spanish) UMSNH (2019). [ Links ]

67. J. S. Ball and T.-W. Chiu Analytic properties of the vertex function in gauge theories. II Phys. Rev. D 22 (1981) 2550. https://doi.org/10.1103/PhysRevD.22.2550 [ Links ]

68. S. Metayer and S. Teber Two-loop mass anomalous dimension in reduced quantum electrodynamics and application to dynamical fermion mass generation e-Print: 2107.07807 ENT#091;hep-thENT#093; https://arxiv.org/abs/2107.07807Links ]

69. A. Ahmad, A. Ayala, A. Bashir, E. Gutiérrez and A. Raya, QCD Phase Diagram and the Constant Mass Approximation, J. Phys. Conf. Ser. 651, no.1, (2015) 012018 https://doi.org/10.1088/1742-6596/651/1/012018 [ Links ]

70. A. Khare, Fractional Statistics and Quantum Theory (World Scientific, Singapore, 2005), 2nd ed. [ Links ]

71. K. I. Kondo and P. Maris, First-order phase transition in three-dimensional QED with Chern-Simons term Phys. Rev. Lett. 74(1995) 18. https://doi.org/10.1103/PhysRevLett.74.18 [ Links ]

72. K. I. Kondo and P. Maris, Spontaneous chiral symmetry breaking in three-dimensional QED with a Chern-Simons term, Phys. Rev. D 52, (1995) 1212. https://doi.org/10.1103/PhysRevD.52.1212 [ Links ]

73. C. P. Hofmann, A. Raya and Saúl Sánchez-Madrigal, Confinement in Maxwell-Chern-Simons planar quantum electrodynamics and the 1/N approximation Phys. Rev. D 82, (2010) 096011. https://doi.org/10.1103/PhysRevD.82.096011 [ Links ]

74. S. Sánchez Madrigal, C. P. Hofmann and A Raya, Dynamical Mass Generation and Confinement in Maxwell-Chern-Simons Planar Quantum Electrodynamics J. Phys.Conf. Ser. 287 (2011) 012028. https://doi.org/10.1088/1742-6596/287/1/012028 [ Links ]

75. R. F. Ozela, V. Alves, G. C. Magalhães and L. O. Nascimento, Effects of the pseudo-Chern-Simons action for strongly correlated electrons in a plane Phys. Rev. D 105 (2022) no.5, 056004 https://doi.org/10.1103/PhysRevD.105.056004 [ Links ]

76. E. C. Marino and Leandro O. Nascimento and Van Sérgio Alves and N. Menezes and C. Morais Smith, Quantum-electrodynamical approach to the exciton spectrum in transition-metal dichalcogenides, 2D Materials, 5, 4 (2018). https://doi.org/10.1088%2F2053-1583%2Faacc3f [ Links ]

77. D. Dudal, A. J. Mizher and P. Pais, Remarks on the Chern-Simons photon term in the QED description of graphene Phys. Rev. D 98 (2018) 065008. https://doi.org/10.1103/PhysRevD.98.065008 [ Links ]

78. G. García Naumis, Electronic properties of 2D materials and its heterostructures: a minimal review, Rev. Mex. Fis. 67 (5), (2021) 1. https://doi.org/10.31349/RevMexFis.67.050102 [ Links ]

79. K.-I. Kondo, Transverse Ward-Takahashi Identity, Anomaly and Schwinger-Dyson Equation, Int. J. Mod. Phys. A12, (1997) 5651. https://doi.org/10.1142/S0217751X97002978 [ Links ]

80. H.-X. He, F. Khanna and Y. Takahashi, Transverse Ward-Takahashi identity for the fermion-boson vertex in gauge theories, Phys. Lett. B480, (2000) 222. https://doi.org/10.1016/S0370-2693(00)00353-1 [ Links ]

81. H.-X. He, Transverse vector vertex function and transverse Ward-Takahashi relations in QED, Commun. Theor. Phys. 46, (2006) 109. https://doi.org/10.1088/0253-6102/46/1/025 [ Links ]

82. H.-X. He, Transverse ward-takahashi relation for the fermion-boson vertex function in four-dimensional abelian gauge theory, Int. J. Mod. Phys. A22, (2007) 2119. https://doi.org/10.1142/S0217751X07036257 [ Links ]

83. K. Johnson and B. Zumino, Gauge Dependence of the Wave-Function Renormalization Constant in Quantum Electrodynamics, Phys. Rev. Lett. 3, (1959) 351. https://doi.org/10.1103/PhysRevLett.3.351 [ Links ]

84. B. Zumino, Gauge properties of propagators in quantum electrodynamics, J. Math. Phys. 1, (1960) 1. https://doi.org/10.1063/1.1703632 [ Links ]

85. S. Okubo, The gauge properties of Green’s functions in quantum electrodynamics, Nuovo Cim. 15, (1960) 949. https://doi.org/10.1007/BF02860201 [ Links ]

86. I. Bialynicki-Birula, On the Gauge Properties of Green’s Functions, Nuovo Cim. 17, (1960) 951. https://doi.org/10.1007/BF02732140 [ Links ]

87. H. Sonoda, On the gauge parameter dependence of QED, Phys. Lett. B 499 (2001) 253. https://doi.org/10.1016/S0370-2693(01)00030-2 [ Links ]

88. L. A. Fernandez-Rangel, A. Bashir, L. X. Gutierrez-Guerrero and Y. Concha-Sanchez, Constructing Scalar-Photon Three Point Vertex in Massless Quenched Scalar QED, Phys. Rev. D 93, no. 6, (2016) 065022. https://doi.org/10.1103/PhysRevD.93.065022 [ Links ]

89. N. Ahmadiniaz, A. Bashir and C. Schubert, Multiphoton amplitudes and generalized LKF transformation in Scalar QED, Phys. Rev. D 93 (2016) 4, 045023 https://doi.org/10.1103/PhysRevD.93.045023 [ Links ]

90. A. Kizilersu and M. R. Pennington, Building the Full Fermion-Photon Vertex of QED by Imposing Multiplicative Renormalizability of the Schwinger-Dyson Equations for the Fermion and Photon Propagators, Phys. Rev. D 79, (2009) 125020. https://doi.org/10.1103/PhysRevD.79.125020 [ Links ]

91. A. Bashir and A. Raya, Constructing the fermion boson vertex in QED(3), Phys. Rev. D 64, (2001) 105001. https://doi.org/10.1103/PhysRevD.64.105001 [ Links ]

92. A.V. Kotikov and S. Teber Landau-Khalatnikov-Fradkin transformation and the mystery of even ζ-values in Euclidean massless correlators Phys. Rev. D 100, (2019) 105017. https://doi.org/10.1103/PhysRevD.100.105017 [ Links ]

93. A. James, A. V. Kotikov, and S. Teber Landau-Khalatnikov-Fradkin transformation of the fermion propagator in massless reduced QED Phys. Rev. D 101 (2020) 045011. https://doi.org/10.1103/PhysRevD.101.045011 [ Links ]

94. A. V. Kotikov and S. Teber, Landau-Khalatnikov-Fradkin Transformation and Hatted ζ-Values Phys. Part.Nucl. 51, (2020) 562. https://doi.org/10.1134/S1063779620040425 [ Links ]

95. V. P. Gusynin, A. V. Kotikov, and S. Teber, Landau-Khalatnikov-Fradkin transformation in three-dimensional quenched QED Phys. Rev. D 102, (2020) 025013 https://doi.org/10.1103/PhysRevD.102.025013 [ Links ]

96. D. Dudal, F. Matusalem, A. J. Mizher, A. R. Rocha and C. Villavicencio, Half-integer anomalous currents in 2D materials from a QFT viewpoint Sci. Rep. 12 (2022) no.1, 5439 https://doi.org/10.1038/s41598-022-09483-4 [ Links ]

i We refer to this scenario throughout the work simply as RQED.

ii In our notation, spatial vectors are written in boldface characters, but their magnitud is detoted by the character alone.

Received: September 21, 2021; Accepted: April 25, 2022

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License