The Effective Potential of Reduced Models in (21)D+

1. Introduction

The experimental realization of two-dimensional materials [1] in condensed matter physics has attracted the interest of the community in high energy physics due to the emergence of Dirac cones and the possibility of observing some interacting effects [2], such as the Fermi velocity renormalization [3,4] and mass renormalization [5,6]. When considering the effects of electronic interactions in these materials, it is useful to consider a dimensional reduction of QED4, namely, PQED, which provides the physical Coulomb potential among static particles in the spatial plane [7]. Within this approach, several results, including strong and weak interactions, have been obtained and compared to experimental data, see Refs. [5,8,9,10,11,12], just to cite a few. Furthermore, the model has also been coined as reduced quantum electrodynamics in Refs. [13,14]. Although the Coulomb potential has been successful in describing these effects, we also expect that this electron-electron interaction become screened due to the microscopic interactions within the two-dimensional material [15]. These interactions should include impurities and phonons, which would make the model even more complicated and, for the best of our knowledge, a complete solution is yet to be known [9]. On the other hand, because the long-range Coulomb interaction is related to a massless photon, hence, a simplest method to generate a short-range potential is to obtain a massive bosonic field, which is straightforwardly derived from SSB. Indeed, when considering the Higgs mechanism in QED4, one would obtain a Yukawa interaction in (3+1)D, yielding a short-range interaction, as expected [16].

The Higgs mechanism for QED in (2+1)D (QED3) would provide the simplest solution for obtaining a mass term for the gauge field and, therefore, a short-range interaction. However, as discussed before, this mechanism must be considered in PQED rather than in QED3. In PQED, the power of the gauge-field propagator is $\propto {\left(p^2\right)}^{-1/2}$ while in QED it is $p^{-2}$, where p is the external momentum. Hence, in order to obtain a Yukawa potential in PQED, one must obtain a scalar propagator given by $1/\sqrt{p^2+m^2}$. This, on the other hand, is obtained from the dimensional reduction of QED4 in the presence of a mass term, as it has been shown in Ref. [17].

Here, we consider the coupling of the gauge field in PQED with a scalar field. This is meant to effectively describe the effects of the many-body system and, essentially, to generate a mass-like term for PQED and to reduce the range of the Coulomb interaction. Firstly, we review the dimensional reduction which generates reduced models in both (2+1)D and (1+1)D. In order to do so, we apply the dimensional reduction in the two-point Green functions of the classical fields, through the equation of motion. Thereafter, we calculate the effective potential for a reduced scalar model in (2+1)D and discuss SSB. This quantity shows a behavior similar to the effective potential in a self-interacting scalar theory in (3+1)D. However, the reduced model does not allow SSB due to its asymptotic limit, leading to a stable ground state in a symmetric phase. Thereafter, we consider effective models for describing the electronic interactions for Dirac-like quasiparticles in (2+1)D in terms of a gauge field interacting with a scalar field. The realization of SSB, however, may occur either in (3+1)D or (2+1)D, providing two different cases for investigation. In case (i), we conclude that the Yukawa interaction is obtained whenever SSB occurs in (3+1)D and derive an effective action for describing such interaction in (2+1)D. In case (ii), which we call as Abelian Higgs PQED (HPQED) in comparison to scalar QED3 [18], we show that the quantum correction does not provide a SSB, hence, the system remains in its symmetric phase. These results are obtained in the one-loop approximation, using the so-called background field method for calculating the effective potential [16].

This paper is organized as follows. In Sec. II, we review the concept of reduced models, using the classical equation of motion. In Sec. III, we calculate the effective potential for a reduced version of the Klein-Gordon theory in (2+1)D. In Secs. IV and V, we consider the effects of considering both SSB and the dimensional reduction in QED4 plus a Higgs field. In Appendixes A–C, we show some details of the calculations.

2. The Reduced Models

In this section, we derive a reduced model that describes the Yukawa interaction in (2+1) dimensions at classical level. In order to do so, we start with the Yukawa action in (3+1) dimensions, whose Euclidean action is given by

$${\mathcal{L}}_{4\mathrm{D}}=\frac{{\partial }_{\mu }\phi {\partial }^{\mu }\phi }{2}+\frac{m^2{\phi }^2}{2}+g\phi \overline{\psi }\psi +{\mathcal{L}}_M\left[\psi \right],$$

(1)

where g is a dimensionless coupling constant, $\phi$ is a real and massive Klein-Gordon field, $\psi$ is the Dirac field, and ${\mathcal{L}}_M\left[\psi \right]$ is the matter action. This model shall be relevant because it will work as a simple example for calculating the dimensional reduction in theories with SSB.

The equation of motion for $\phi$ is promptly obtained from Eq. (1) and reads

$$(-\square +m^2)\phi \left(x\right)=J\left(x\right),$$

(2)

where □ is the d’Alembertian operator and $J\left(x\right)\equiv -g\overline{\psi }\psi \left(x\right)$ works as an external source for the scalar field. The solutions of Eq. (2) are obtained by inverting the differential operator, i.e,

$$\phi \left(x\right)=\int d^{4}y{G}_{4D}(x-y)J\left(y\right),$$

(3)

where

$$G_{4D}(x-y)=\int \frac{d^4{k}_E}{{\left(2\pi \right)}^4}\frac{e^{ik(x-y)}}{k_E^2+m^2}$$

(4)

is the Fourier transform of the scalar-field propagator. Equation (4), within the static-limit, provides the Yukawa potential. Hence, after using $k_E^2=k_0^2+k_1^2+k_2^2+k_3^2$ with $k_0\to 0$ [16], we find

$$V\left(r\right)=\int \frac{d^3{k}_E}{{\left(2\pi \right)}^3}\frac{e^{ikr}}{k_E^2+m^2}=\frac{e^{-\left|m\right|r}}{4\pi r},$$

(5)

where the inverse of m is the interaction length of the model and $r^2={(x_1-y_1)}^2+{(x_2-y_2)}^2+{(x_3-y_3)}^2$, as expected.

2.1. From (3+1)D to (2+1)D

For calculating a reduced model, we consider that the dynamics of the matter field is constrained to the (2+1)-dimensional space-time. This is obtained by assuming that $J\left(y\right)=J(y_0,y_1,y_2)\delta \left(y_3\right)\equiv J_{3D}\left(y\right)\delta \left(y_3\right)$. Thereafter, we perform a dimensional reduction of the two-point Green function of the scalar field. Hence, we obtain the scalar-field propagator at $x_3=0$, i.e, $\phi \left(x\right)\to {\phi }_{3D}\left(x\right)\equiv \varphi (x_0,x_1,x_2)=\phi (x_0,x_1,x_2,x_3=0)$. Furthermore, we also adopt the notation, where the x in $\varphi \left(x\right)$ always do represent the coordinates in the spacetime where the field is defined, hence, for $\varphi \left(x\right)$, we have $x=(x_0,x_1,x_2)$. On the other hand, for $\phi \left(x\right)$, it follows that $x=(x_0,x_1,x_2,x_3)$. The same holds for loop integrals and propagators in momentum space. After using these conditions in Eq. (3), we find

$${\phi }_{3D}\left(x\right)=\int d^{3}y{G}_{3D}(x-y)J_{3D}\left(y\right),$$

(6)

where

$$G_{3D}(x-y)=\int \frac{d^3{k}_E}{{\left(2\pi \right)}^3}\frac{e^{ik(x-y)}}{2\sqrt{k_E^2+m^2}}$$

(7)

is the Fourier transform of the scalar-field propagator in (2+1) dimensions. This also procudes the Yukawa potential in Eq. (5). Indeed, after using $k_E^2=k_0^2+k_1^2+k_2^2$ with $k_0\to 0$ in Eq. (7), we find

$$V\left(r\right)=\int \frac{d^2{k}_E}{{\left(2\pi \right)}^2}\frac{e^{ikr}}{2\sqrt{k_E^2+m^2}}=\frac{e^{-\left|m\right|r}}{4\pi r},$$

(8)

where $r^2={(x_1-y_1)}^2+{(x_2-y_2)}^2$.

The main idea of a reduced model is to obtain a theory that yields the reduced propagator in Eq. (7). This may be obtained from

$${\mathcal{L}}_{3D}=\frac{{\partial }^{\mu }\varphi K_E[\square ]{\partial }_{\mu }\varphi }{2}+g\varphi \overline{\psi }\psi +{\mathcal{L}}_M\left[\psi \right]$$

(9)

with

$$K_E[\square ]\equiv \frac{2\sqrt{-\square +m^2}}{-\square }=\int \frac{d^3{k}_E}{{\left(2\pi \right)}^3}{e}^{ikx}\frac{2\sqrt{k_E^2+m^2}}{k_E^2}.$$

(10)

This model describes the dynamical effects of the Yukawa interaction in the plane and some of its quantum effects have been discussed in Ref. [17] when coupled to the Dirac field. Obviously, the Coulomb interaction is promptly obtained by using $m\to 0$ within these calculations. The case of dimensional reduction to (1+1)D is discussed in App. A.

The derivation of the reduced models is more convenient to be performed in the Euclidean spacetime [7]. However, for calculating the Coleman-Weinberg potential, we rather consider the action in the Minkowski space for the sake of comparison with well known results in the literature.

3. The One-Loop Effective Potential of a Reduced-Scalar Model in (2+1)D

In this section, we discuss the Coleman-Weinberg effective potential of a general scalar action in a D-dimensional spacetime with a Pseudo-diferential operator. When considering an interaction term $U\left[\varphi \right]$, our toy model action reads

$${\mathcal{L}}_D=\frac{{\partial }^{\mu }\varphi \mathcal{K}[\square ]{\partial }_{\mu }\varphi }{2}+U\left[\varphi \right],$$

(11)

where the kernel $\mathcal{K}[\square ]$ is an arbitrary pseudo-differential operator in the Minkowski spacetime. The calculations using Eq. (11) shall be useful for our purposes in the next sections, where we discuss the SSB for two-dimensional materials in the light of the PQED formalism.

The partition function of the model is defined as usual and, therefore, reads

$$\mathcal{Z}=\int D\phi exp\left\{i\int d^{D}x{\mathcal{L}}_D\left[\varphi \right]\right\}.$$

(12)

The main idea of the one-loop expansion is to expand the action in Eq. (12), using Eq. (11), up to order of ℏ, see App. B for more details. Thereafter, we replace the quantum field by its vacuum expectation value (VEV), i.e, $\varphi \left(x\right)\to \langle \varphi \left(x\right)\rangle =\rho$, which is acceptable for describing the ground state of the model through its effective potential $V_{\mathrm{eff}}\left(\rho \right)$. Therefore, we have

$$\mathcal{Z}=e^{-i{V}_{\mathrm{eff}}\left(\rho \right)\Omega },$$

(13)

where $\Omega =\int d^{3}x$ is a constant factor, namely, the spacetime volume. Having these steps in mind and, after some algebra, one finds

$$V_{\mathrm{eff}}\left(\rho \right)=-U\left(\rho \right)-\frac{i}{2}\int \frac{d^{D}k}{{\left(2\pi \right)}^D}ln[k^2\mathcal{K}\left(k^2\right)+U^{\prime \prime }\left(\rho \right)],$$

(14)

where $U^{\prime \prime }\left(\rho \right)$ means the second derivative of $U\left[\rho \right]$. Note that the first term in the rhs of Eq. (14) is the classical potential and the second term is the quantum correction (in order of ℏ). Eq. (14) is how far we may go for an arbitrary kernel $\mathcal{K}$.

3.1. The case $D=4$ and $\mathcal{K}=1$

This is the most standard case and the effective potential reads [16]

$$V_{\mathrm{eff}}^{\mathrm{KG}}(\rho ,\mu )=-U\left(\rho \right)-\frac{{\left(U^{\prime \prime }\right)}^2}{64{\pi }^2}ln\left[\frac{4\pi {\mu }^2}{-U^{\prime \prime }}\right],$$

(15)

which is derived through the dimensional regularization scheme. Note that $\mu$ is a renormalization point and that $V_{\mathrm{eff}}^{\mathrm{KG}}\left(\rho \right)$ is expected to be unchanged by a scale transformation in $\mu$. For a classical potential $U\left[\rho \right]=\lambda {\rho }^4/4!$ with $\lambda <0$, we conclude that the VEV of the field does not vanish and it is double degenerated. hence, the discrete symmetry $\varphi \to -\varphi$ is broken by the ground state [16]. Our main goal is to understand what are the effects of considering a pseudo-diferential operator in a scalar action, regarding the mechanism of SSB. Thereafter, we apply these results for PQED with a Higgs field.

3.2. The case $D=3$ and $\mathcal{K}=2/\sqrt{-k^2}$

This is the case where $\mathcal{K}[\square ]=2{(\square )}^{-1/2}$ in Eq. (11). This corresponds to a reduced Klein-Gordon model. Here, we calculate its effective potential, using both the dimensional regularization and the cutoff regularization. As it shall be clear later, both these methods provide the same logarithmic term for the effective potential.

3.2.1. The dimensional regularization scheme

Within this regularization scheme, Eq. (14) is written as

$$\begin{array}{ccc}\hfill V_{\mathrm{eff}}(\rho ,\mu )& =& -U\left(\rho \right)-\frac{i{\mu }^{3-D}}{2}\int \frac{d^{D}k}{{\left(2\pi \right)}^D}ln[k^2\mathcal{K}\left[k^2\right]+U^{\prime \prime }\left(\rho \right)]\hfill \\ & \equiv & -U\left(\rho \right)+V_1(\rho ,\mu ),\hfill \end{array}$$

(16)

where $\mu$ is an energy scale included in order to preseve the units of the integral and $V_1$ is the quantum correction. In App. C we calculate this correction and obtain

$$\begin{array}{ccc}\hfill V_{\mathrm{eff}}(\rho ,\mu )& =& -U\left(\rho \right)+\frac{{(1/2)}^{D-1}\mathsf{\Gamma }\left[D\right]2{\pi }^{D/2}{(-U^{\prime \prime })}^3}{{\left(2\pi \right)}^{D}4D\mathsf{\Gamma }[D/2]}\hfill \\ & \times & {\left(\frac{\mu }{-U^{\prime \prime }}\right)}^{(3-D)}\mathsf{\Gamma }[1-D].\hfill \end{array}$$

(17)

Expanding Eq. (17) around $D\to 3$, we have

$$\begin{array}{ccc}\hfill V_{\mathrm{eff}}(\rho ,\mu )& =& -U\left(\rho \right)-\frac{{\left(U^{\prime \prime }\right)}^3}{48{\pi }^2}\hfill \\ & \times & \left[\frac{1}{2(3-D)}+\frac{3}{4}-\frac{\gamma }{2}+\frac{1}{2}ln\left(\frac{\mu }{-U^{\prime \prime }}\right)\right],\hfill \end{array}$$

(18)

where $\gamma \approx 0.58$ is the Euler’s constant. Finally, we use the minimal subtraction scheme (where the pole of the Gamma function is neglected as well as some extra constants) in Eq. (18) and find the renormalized effective potential, namely,

$$V_{\mathrm{eff}}(\rho ,\mu )=-U\left(\rho \right)-\frac{{\left(U^{\prime \prime }\right)}^3}{96{\pi }^2}ln\left(\frac{\mu }{-U^{\prime \prime }}\right).$$

(19)

Note that Eq. (19) resembles the effective potential of the scalar model in (3+1)D, given by Eq. (15). In particular, it has an energy-scale term which depends on $\mu$. This is a surprisingly conclusion, because in an usual scalar theory in (2+1)D, such term does not exist. Obviously, the pseudo-differential operator plays a central role for this result. Next, let us consider the cutoff regularization scheme, which shall be useful when we include a finite lattice regulator, as it occurs for effective models in condensed matter physics [16].

3.2.2. The cutoff regularization scheme

In this case, we also perform the Wick rotation to the Euclidean spacetime, i.e, $d^{D}k\to i{d}^D{k}_E$ and $k^2\to -k_E^2$. Thereafter, we include an ultraviolet cutoff $\mathsf{\Lambda }$ in the integral, where $0\le k_E\le \mathsf{\Lambda }$ in Eq. (14). Furthermore, for solving this integral, we use that

$$\begin{array}{ccc}\hfill \int k_E^{2}d{k}_{E}ln[2k_E-U^{\prime \prime }]& =& -\frac{1}{12}\left(k_E^2{U}^{\prime \prime }+k_E{U}^{\prime \prime 2}\right)\hfill \\ & -& \frac{U^{\prime \prime 3}}{24}ln[2k_E-U^{\prime \prime }]\hfill \\ & +& \frac{k_E^3}{3}ln[2k_E-U^{\prime \prime }]-\frac{k_E^3}{9}.\hfill \end{array}$$

(20)

Next, using Eq. (20) in Eq. (14) and expanding this for $\mathsf{\Lambda }\gg U^{\prime \prime }$, we find

$$\begin{array}{ccc}\hfill V_{\mathrm{eff}}(\rho ,\mathsf{\Lambda })& =& -U\left(\rho \right)-\frac{U^{\prime \prime 3}}{288{\pi }^2}-\frac{U^{\prime \prime 2}\mathsf{\Lambda }}{32{\pi }^2}-\frac{U^{\prime \prime }{\mathsf{\Lambda }}^2}{16{\pi }^2}\hfill \\ & -& \frac{{\left(U^{\prime \prime }\right)}^3}{96{\pi }^2}ln\left(\frac{2\mathsf{\Lambda }}{-U^{\prime \prime }}\right).\hfill \end{array}$$

(21)

Note that the coefficient of the logarithmic term, i.e, $-{\left(U^{\prime \prime }\right)}^3/96{\pi }^2$ is the same as in Eq. (19). This is an relevant feature for calculating the beta functions of the theory in Eq. (11), when considering the continuum limit, i.e, $\mathsf{\Lambda }\to \infty$[16]. Here, because the dimensional reduction is related to two-dimensional materials, we consider the case when $\mathsf{\Lambda }\propto 1/a$ is a large but finite constant, where $a\approx {10}^{-10}$m is the lattice parameter for a typical crystal. For example, in graphene and other two-dimensional materials, we have $\mathsf{\Lambda }\approx 1$eV [1].

3.2.3. The vacuum stability of the reduced-scalar model with a ${\varphi }^4$-self-interaction

For a more straightforward application, let us consider a classical potential given by

$$U\left[\varphi \right]=\frac{M^2{\varphi }^2}{2}+\frac{\lambda {\varphi }^4}{4!},$$

(22)

where M is the bare mass term and $\lambda$ is the bare coupling constant. The analysis of the vacuum stability of Eq. (21) is similar to the stability given by the effective potential in Eq. (15), where we consider $\lambda <0$ in order to have a real-valued effective potential. In Figure 1, we plot Eq. (21), using Eq. (22). From this, we may conclude that the ground state of the system remains symmetric at $\rho =0$, hence, no SSB occurs.

Although it is interesting that the reduced model has a similar structure for the effective potential, in comparison to the scalar field in (3+1)D, a few comments have to be addressed. (i) Our example holds for the $Z_2$ symmetry $\varphi \to -\varphi$, however, SSB is defined for continuous symmetries, hence, our calculation works only as a first approximation for calculating SSB. Despite this, the generalization for a complex field is simple and the results are clearly the same; (ii) The effective potential of the scalar field in (2+1)D does not have the logarithmic term [18], which is a striking difference in comparison to both the higher dimensional model (see Eq. (15)) and our reduced-scalar model (see Eq. (21); (iii) Because we are assuming a finite $\mathsf{\Lambda }$, hence, our effective potential is also finite and we do not need to deal with divergencies.

Next, we shall consider a scalar version of PQED with a complex scalar field. This also provides a more physical situation, allowing us to make predictions regarding the interactions between quasiparticles in two-dimensional materials.

4. From SSB to the Dimensional Reduction

The static approximation of PQED describes a system of electrons interacting through the Coulomb potential. However, due to screening effects within the two-dimensional material, it is expected that the Coulomb potential becomes a short-range interaction [19]. This mechanism, however, must be considered either before or after the dimensional reduction of the Gauge field. This lead us to two important cases, which we shall consider next. First, let us assume that the SSB occurs in (3+1)D and then we perform the dimensional reduction. Thereafter, we consider the opposite case, where the Higgs field is included in PQED.

We start with QED4 and include a Higgs field, hence, our action in the Euclidean spacetime reads

$${\mathcal{L}}_{\mathrm{HQED}4}=\frac{1}{4}{\overline{F}}^{\mu \nu }{\overline{F}}_{\mu \nu }-\frac{1}{2\alpha }{\left({\partial }_{\mu }{\overline{A}}^{\mu }\right)}^2+|D_{\mu }{\phi |}^2+U\left[\right|\phi \left|\right]+\overline{\psi }(i\cancel{\partial }-m)\psi +e\overline{\psi }{\gamma }_{\mu }\psi {\overline{A}}^{\mu },$$

(23)

where $D_{\mu }\to {\partial }_{\mu }+ig{\overline{A}}_{\mu }$, with $\mu =0,1,2,3$, is the covariant derivative and g is the dimensionless charge of the complex scalar field $\phi$. ${\overline{A}}_{\mu }$ is the gauge field and ${\overline{F}}_{\mu \nu }={\partial }_{\mu }{\overline{A}}_{\nu }-{\partial }_{\nu }{\overline{A}}_{\mu }$ is its field strength tensor. Furthermore, $\alpha$ is the gauge fixing parameter. Because of charge conservation, this parameter is not relevant for our purposes, hence, we simply neglect this term. $\phi$ is a massive complex scalar field coined Higgs field and $U\left[\right|\phi \left|\right]$ represents its self-interacting potential. $\psi$ is the Dirac field, which is meant to describe the quasi-particles in the two-dimensional material [5,9], and m is its bare mass while e is the dimensionless electric charge.

The standard approach for SSB in Eq. (23) is to consider a polar representation of the scalar field, i.e, $\phi \left(x\right)=\xi \left(x\right)e^{i\theta \left(x\right)}$. Furthermore, we assume that $U\left[\right|\phi \left|\right]={\mu }^2{\xi }^2/2+\lambda {\xi }^4/4!$, where $(\mu ,\lambda )$ are known constants. After some algebra, we conclude that the scalar field $\theta \left(x\right)$ may be removed from the model by a gauge transformation, namely, ${\overline{A}}_{\mu }\to {\overline{A}}_{\mu }-{\partial }_{\mu }\theta /g$ [16]. Hence, we find

$${\mathcal{L}}_{\mathrm{HQED}4}=\frac{1}{4}{\overline{F}}^{\mu \nu }{\overline{F}}_{\mu \nu }+{\partial }_{\mu }\xi {\partial }^{\mu }\xi +U\left[\xi \right]+g^2{\xi }^2{\overline{A}}_{\mu }{\overline{A}}^{\mu }+\overline{\psi }(i\cancel{\partial }-m)\psi +e\overline{\psi }{\gamma }_{\mu }\psi {\overline{A}}^{\mu }.$$

(24)

Finally, we consider that the Higgs field is in the broken phase and take its lowest-order solution, given by $\xi \left(x\right)\to {\rho }_0$. Therefore,

$${\mathcal{L}}_{\mathrm{HQED}4}=\frac{1}{4}{\overline{F}}^{\mu \nu }{\overline{F}}_{\mu \nu }+g^2{\rho }_0^2{\overline{A}}_{\mu }{\overline{A}}^{\mu }+\overline{\psi }(i\cancel{\partial }-m)\psi +e\overline{\psi }{\gamma }_{\mu }\psi {\overline{A}}^{\mu }.$$

(25)

Next, we use Eq. (25) for calculating the dimensional reduction. Note that whenever ${\rho }_0\ne 0$ we may conclude that the Gauge field has acquired a mass-like term. The dimensional reduction, on the other hand, will follow exactly the same steps as we did in Sec. II. This also may be performed through the generating functional of the current-current correlation functions, as it has been done in Ref. [7] for PQED. Hence, it follows that the reduced version of Eq. (25) is

$$\begin{array}{ccc}\hfill {\mathcal{L}}_{\mathrm{HRQED}}& =& \frac{1}{4}{\overline{F}}^{\mu \nu }{K}_E[\square ]{\overline{F}}_{\mu \nu }+\overline{\psi }(i\cancel{\partial }-m)\psi \hfill \\ & +& \frac{1}{2\alpha }{\overline{A}}^{\mu }{\partial }_{\mu }{K}_E[\square ]{\partial }_{\nu }{\overline{A}}^{\nu }+e\overline{\psi }{\gamma }_{\mu }\psi {\overline{A}}^{\mu },\hfill \end{array}$$

(26)

which is fully definied in (2+1)D and

$$K_E[\square ]=\frac{2\sqrt{-\square +2g^2{\rho }_0^2}}{-\square }.$$

(27)

Despite our notation, the gauge field in Eq. (26) is not the same as in Eq. (23). This model describes the Yukawa interaction between Dirac-like electrons that are confined to move, for example, in a two-dimensional material [17]. Furthermore, because of the SSB, the interaction become short-ranged and given by

$$V\left(r\right)=\frac{e^{-\sqrt{2}g{\rho }_{0}r}}{4\pi r},$$

(28)

which is obtained similarly to the calculations in Sec. II. The mechanism we use in this section does not allow us to discuss further about SSB. Indeed, it occurs in the higher dimensional model and we simply performed the dimensional reduction. Nevertheless, we could also consider that the Higgs mechanism is driven in (2+1)D instead of in (3+1)D. In this new case, we should also ask whether the effective potential is minimized in the broken phase or not.

5. From Dimensional Reduction to SSB

In this section, we consider a different path for discussing SSB in PQED. Here, we start with PQED and then include a Higgs field, in the Minkowski spacetime, for calculating the SSB. In this case, we may use some results in Sec. III regarding the loop integrals. Therefore, we have

$${\mathcal{L}}_{\mathrm{HPQED}}=-\frac{1}{2}{F}^{\mu \nu }{(\square )}^{-1/2}{F}_{\mu \nu }+|D_{\mu }{\varphi |}^2-\frac{1}{2\alpha }\frac{{\left({\partial }^{\mu }{A}_{\mu }\right)}^2}{{(\square )}^{1/2}}+U\left[\right|\varphi \left|\right],$$

(29)

where the subscript HPQED stands for Higgs-PQED, $D_{\mu }\to {\partial }_{\mu }+ie{A}_{\mu }$, with $\mu =0,1,2$, is the covariant derivative, and e is the dimensionless electric charge of the complex scalar field $\varphi$. $A_{\mu }$ is the gauge field of PQED, $F_{\mu \nu }={\partial }_{\mu }{A}_{\nu }-{\partial }_{\nu }{A}_{\mu }$ is its field strength tensor, and $\alpha$ is the gauge-fixing parameter. $\varphi$ is a massive complex scalar field, coined Higgs field, and $U\left[\right|\varphi \left|\right]$ represents its symmetry-breaking potential. We neglect the Dirac term, because it does not play any role regarding the SSB and may be included later.

The model in Eq. (29) has been investigated in the perturbation theory in Ref. [19], within the symmetric phase. Here, we shall discuss its broken phase in the light of the effective potential. For the sake of simplicity, we assume the Feynman gauge, where the longitudinal part of the gauge-field propagator vanishes. As matter of fact, the effective potential in scalar QED4 is also gauge-dependent. Nevertheless, the physical predictions, such as the ratio between the masses of both gauge and scalar fields, are expected to independent on $\alpha$ [16]. Using the one-loop approximation, we may expand the action as

$$\begin{array}{ccc}\hfill {\mathcal{L}}_{\mathrm{HPQED}}& =& A^{\mu }{(\square )}^{1/2}{A}_{\mu }-ie{A}_{\mu }[{\varphi }^{*}\left({\partial }^{\mu }\varphi \right)-\left({\partial }^{\mu }{\varphi }^{*}\right)\varphi ]\hfill \\ & -& {\varphi }^{*}\square \varphi +e^2{A}_{\mu }{A}^{\mu }{\left|\varphi \right|}^2+U\left[\right|\varphi \left|\right].\hfill \end{array}$$

(30)

Next, we consider the one-loop approximation for calculating the effective potential. Similar to the calculation in Sec. III, we assume that the scalar field may have a nonzero VEV, such that $\langle \varphi \rangle =\rho$ and $\langle A_{\mu }\rangle =0$. Using these assumptions in Eq. (30), and after neglecting some irrelevant constants, we find

$$\begin{array}{ccc}\hfill V_{\mathrm{eff}}^{\mathrm{HPQED}}\left(\rho \right)& =& -U\left(\rho \right)-i\int \frac{d^{3}k}{{\left(2\pi \right)}^3}ln[k^2+U^{\prime \prime }]\hfill \\ & -& \frac{3i}{2}\int \frac{d^{3}k}{{\left(2\pi \right)}^3}ln[{(-k^2)}^{1/2}+e^2{\rho }^2].\hfill \end{array}$$

(31)

Eq. (31) is the main result of this section and should be compared to the standard scalar-QED in (2+1)D [18]. In order to do so, we must replace ${(-k^2)}^{1/2}\to {(-k)}^2$ for describing the Maxwell propagator. In this case, we may conclude that the quantum corrections to the effective potential are finite in the sense that they do not have the logarithmic term, which depends either on $\mu$ or $\mathsf{\Lambda }$, depending on the regularization scheme. Here, nevertheless, the third term in rhs of Eq. (31) provides a logarithmic term, such as in the previous case of the reduced-scalar model. This term is clearly related to the pseudo-differential operator.

5.1. The Higgs term

The second term in the rhs of Eq. (31) is the Higgs-field contribution to the effective potential. Within the cutoff regularization, after going to the Euclidean space, it yields

$$\frac{1}{2{\pi }^2}{\int }_0^{\mathsf{\Lambda }}{k}_E^{2}d{k}_{E}ln\left[k_E^2-U^{\prime \prime }\right].$$

(32)

After solving the integral over $k_E$, keeping only the $\rho$-dependent terms, and expanding for $\mathsf{\Lambda }\gg U^{\prime \prime }$, we find

$$-\frac{\left(U^{\prime \prime }\mathsf{\Lambda }\right)}{2{\pi }^2}-\frac{{(-U^{\prime \prime })}^{3/2}}{6\pi }.$$

(33)

A Maxwell-like propagator for the gauge field would essentially provide a similar result.

5.2. The PQED term

The third term in the rhs of Eq. (31) is the PQED contribution. Using the same assumptions as before, it yields

$$\begin{array}{c}\hfill \frac{3}{4{\pi }^2}{\int }_0^{\mathsf{\Lambda }}{k}_E^{2}d{k}_{E}ln\left[k_E+e^2{\rho }^2\right].\end{array}$$

(34)

Solving the integral over $k_E$ and expanding for $\mathsf{\Lambda }\gg e^2{\rho }^2$, see Eq. (20), we have

$$\frac{{\left(e^2{\rho }^2\right)}^3}{12{\pi }^2}-\frac{3\mathsf{\Lambda }{e}^4{\rho }^4}{8{\pi }^2}+\frac{3{\mathsf{\Lambda }}^2{e}^2{\rho }^2}{8{\pi }^2}+\frac{3{\left(e^2{\rho }^2\right)}^3}{24{\pi }^2}ln\left(\frac{{\mathsf{\Lambda }}^2}{{\left(e^2{\rho }^2\right)}^2}\right).$$

(35)

Note that due to the pseudo-differential operator we have obtained a logarithmic term that resembles the standard effective potential in the scalar QED in (3+1)D [16].

5.3. The Effective Potential

Finally, the effective potential, after taking the two terms, reads

$$\begin{array}{c}\hfill V_{\mathrm{eff}}^{\mathrm{HPQED}}(\rho ,\mathsf{\Lambda })=-U\left(\rho \right)-\frac{\left(U^{\prime \prime }\mathsf{\Lambda }\right)}{2{\pi }^2}-\frac{{(-U^{\prime \prime })}^{3/2}}{6\pi }+\frac{{\left(e^2{\rho }^2\right)}^3}{12{\pi }^2}-\frac{3\mathsf{\Lambda }{e}^4{\rho }^4}{8{\pi }^2}\\ \hfill +\frac{3{\mathsf{\Lambda }}^2{e}^2{\rho }^2}{8{\pi }^2}+\frac{3{\left(e^2{\rho }^2\right)}^3}{24{\pi }^2}ln\left(\frac{{\mathsf{\Lambda }}^2}{{\left(e^2{\rho }^2\right)}^2}\right).\end{array}$$

(36)

Next, let us discuss the possibility of finding a symmetry-breaking solution in Eq. (36). Let us consider a classical potential, given by Eq. (22) with $M^2=0$. In this case, we find a stable ground state at $\rho =0$ and the effective potential resembles Fig. Figure 1. Otherwise, when $\lambda >0$, we do not have a real-valued potential.

Although our results indicate that HPQED does not admits a SSB at one-loop approximation, we could also imagine that the scalar field has a nonzero VEV at tree level, i.e, an explicit breaking of symmetry. In this case, the interaction between the electrons is given by a two-dimensional Fourier transform of $1/(2\sqrt{p^2}+2e^2{\rho }_0^2)$. Indeed, such result has been considered in Ref. [20] in the light of superconductivity driven by a topological phase transition. Clearly, this is quite different in comparison to the Yukawa potential we have discussed in Sec. IV. Indeed, it is not hard to conclude that the resulting potential behaves as $ln\left(r\right)$ for $r\to 0$ and goes to zero as $r\to \infty$. This only proves the need of taking care between the order of the dimensional reduction and SSB.

6. Summary

The major relevance of SSB is in describing quantum states of matter, such as ferromagnetism and superconductivity. This is typically made by considering a Landau-Ginzburg theory [21]. For field theories, scalar QED4 is a relativistic generalization of the Landau-Ginzburg theory, where the Dirac field of QED4 is replaced by a scalar field with a self-interacting term [16]. However, when considering a (2+1)D version of scalar QED, one concludes that the effective potential behaves quite different in comparison with the higher-dimensional version of the model. For example, its effective potential trivially obeys the renormalization group equation [18]. Indeed, because the spatial dimension is reduced by one, hence, the ultraviolet divergences are reduced and all of the beta functions vanish.

In order to have a proper description of the electronic interactions between electrons in a condensed-matter system, we must consider the reduced version of QED4, namely, PQED [7]. This dimensional reduction also may be performed to preserve the Yukawa interaction in the plane [17]. Here, we show that the SSB does not occur in PQED, using the one-loop approximation within the background field method. Although we have considered a real scalar field, the generalization for the complex field is straightforward. Thereafter, we consider the Abelian version of scalar PQED in (2+1)D. In this case, we must be careful about the order in which we realize both the SSB and the dimensional reduction. When the SSB is considered in scalar QED4, hence, after the dimensional reduction, we obtain a version of PQED that describes the Yukawa interaction between electrons in (2+1)D. Nevertheless, starting with PQED and coupling it to a scalar field, we realize that SSB is not realized and the system remains in its symmetric phase. Obviously, our approach is not the only path for describing screening effects and it only has the virtue that it follows from a quite simple and well known method, namely, the coupling of the gauge and scalar fields. Indeed, screening effects due to the fermionic loop are also considered in literature of quantum field theory [16] and may also be relevant for two-dimensional materials [10].

It would be interesting to investigate the SSB in the Abelian PQED where both the gauge and scalar fields have a pseudo-differential operator as well as the effect of including a thermal bath. We shall discuss this elsewhere.