# Vacuum polarization at the boundary of a topological insulator

###### Abstract

In this paper we study the polarized vacuum energy on the conducting surface of a topological insulator characterized by both topological index and time reversal symmetry. This boundary is subject to the action of a static and spatially homogeneous magnetic field perpendicular to it as well as of an electric field that is uniform near the considered surface and produced by a biased voltage, at zero temperature. To do this, we consider modifications in the Gauss law that arise due to the nonzero gradient of the axionlike pseudoscalar factor coupled to the applied magnetic field, which accounts for the topological properties of the system. Such a term allows us to find a correction to the induced charges which modifies the quantum vacuum of the spinor field regarding an ordinary surface. The polarized vacuum energy is calculated in both the weak-field approximation and in the general case, and since the found energy depends on a length defined on the boundary, we show that there is a radial density of force or a surface shear stress that tends to shrink it.

Key words: Vacuum polarization, Axion electrodynamics, Topological insulator.

## I Introduction

Topological insulators (TIs) are electronic systems with a nonzero energy gap to the excited (conducting) states in the bulk, nevertheless possessing a gapless mode on their boundary surfaces (edges) Hasan (). The meaning of the term “topological” comes from the fact that the wave functions describing their electronic states span a Hilbert space with a nontrivial topology ryu (). Such systems are also characterized by a topological number that reveals how robust they are against external perturbations, including those ones related to the quantum vacuum.

A lot of theoretical work has been done since the original formulation by Casimir Casimir () concerning the macroscopic effect which is caused by modifications in the fluctuations of the quantum vacuum energy of the electromagnetic field and that bears his name. The effect originally described consists of an attraction force that arises between two ideal conducting plates closely spaced, parallel and uncharged, placed in a perfect vacuum at zero temperature. This phenomenon was experimentally confirmed ten years later with some level of precision. Numerous papers published since then have addressed various types of geometrical settings, quantum fields, materials, space-times, topologies and thermal conditions. These and the latest developments both theoretical and experimental, are described in detail in Bordag (); Miloni ().

Effects related to the quantum vacuum can also be present in empty spaces with non-Euclidean topology including cosmological models, in general relativity or alternative theories of gravity Celio1 (); Celio2 (). In these space-times there are no material borders, but the identification of conditions imposed on the quantum fields can play the same role as the one presented by these boundaries. When macroscopic fields are present, possible physical effects can arise due to the polarization of the quantum vacuum induced by these fields. Thus, it is interesting to verify how such effects operate over systems as topological insulators, in which the nontrivial topology is in their spaces of quantum states instead of in their physical spaces. In this investigation line, the literature presents some works Adolfo (); Rodriguez (); Alberto ().

In the present work we will study the polarized quantum vacuum associated to 2D charged fermions when an uniform electric field generated by a biased voltage and an equally uniform magnetic field act on the surface of a special class of TI’s, labeled by a topological number the index, which expresses the number of times that a 1D edge electronic state crosses the Fermi level between 0 and , where is the lattice parameter, besides exhibiting time reversal (TR) symmetry (e.g., BiSe). In these systems there is a term in the electromagnetic Lagrangian proportional to , derived from the axion electrodynamics wilck (), which, inside the TI, does not modify the usual Maxwell’s laws valid in material media by virtue of the factor being a constant, despite changing the electromagnetic constitutive relations. This leads, for example, to the appearance of magnetization (polarization) induced by the electric (magnetic) field ando (). It is interesting to point out that, according to Ref. Li (), it is possible to study the breakdown of TR symmetry via magnetic fluctuations in order to acquire dynamics, becoming a field variable.

Although in our approach is not a dynamical field, on the TI conducting surface the Maxwell’s laws are modified since . Thus, we will consider the variation in the Gauss’ law due to the nonzero gradient of which is coupled to a static and spatially homogeneous magnetic field perpendicular to the TI boundary in order to find an effective electric charge that modifies the polarized vacuum energy on that edge, firstly in the weak-field limit, where the energy gap is approximately zero. Next, we will analyze the more general situation by considering the nonzero gap opened through incidence of stronger magnetic fields. In this case, we will use the result found in Visser () in order to analyze the polarized quantum vacuum of a spinor field in 1+2 dimensions. Since this new vacuum energy depends on a length defined on the boundary, there is a tangential force or a shear stress associated to it.

The present paper unites frontier areas of research, namely polarized vacuum effects and topological insulators, investigating a phenomenon that can be relevant in the technological applications of TIs, since the work predicts the occurrence of a dynamical effect due to the quantum vacuum which is possibly non-negligible.

The paper is organized as follows: In Sec. II, we present the electrodynamics of a TI with topology and TR symmetry, obtaining a modification in the electric charges on its conducting boundary. In Sec. III, we reobtain the polarized vacuum energy density of a spinor field endowed with the found effective charge and interacting with a constant magnetic field on that boundary, in both the weak-field regimen and general case. Finally, in Sec. IV, we discuss the results.

## Ii Correcting the Induced Electric Charge

In this section, we will study how variations in the electrodynamics laws induced by the topology of TIs modify in turn the electric charge of the carriers. The electromagnetic Lagrangian density in a -type TI with TR symmetry can be written as wilck ()

(1) |

without source terms. The factor is the fine structure constant. The second term of Eq. (1) - the -term, which encodes the topological features of TIs, is quite similar to the one which occurs in the axion electrodynamics, where is a pseudoscalar field, for which its transformation in a temporal reversion, given by , is , and whose coupling with the electromagnetic field generates unusual effects as axion-photon oscillations raffelt () and shining light through walls Javier ().

The presence of the -term changes the electromagnetic constitutive relations, in the form

(2) | |||

(3) |

in a such way that within the TI bulk these quantities lead to the so-called magnetoelectric effect, in which the polarization is proportional to the magnetic field or the magnetization is proportional to the electric field. This effect can be exploited in order to obtain other exotic phenomena such as magnetic monopole charge-image qi (). In that region inside TI, is a constant of value equal to , and thus the electrodynamics laws in material media become unchanged because the extra term is just a total derivative Essin ().

On the boundary of the TIs under consideration, the Maxwell’s equations are modified besides the constitutive relations, because the axionlike factor presents a nonzero gradient, since in regions external to the TI (vacuum or ordinary insulators), , and inside it, . Thus, the Gauss’s law changes to

(4) |

The electric field is necessary in our analysis for two reasons: Firstly, the theta term does not exist without it. Secondly, it induces the electric charge concentration on the conducting surface. In fact, the experiments with BiSe described in Ref. Tong () used magnetic fields as well as biased electric fields in order to study the Landau quantization on the TI boundary, whose energy levels must be taken into account when one calculates the vacuum energy Beneventano (). Then the total charge density is a combination of that one induced by the external electric field, , with the one of topological origin due to the theta term, .

Thus, after the integration of Eq. (4) throughout a thin pill-box-type (cylindrical) volume encompassing the TI surface and within the limit of the fields action, we get

(5) |

where is the effective (topological) charge and is the magnetic flux on an area contained in the TI surface. Without the magnetic flux the induced charge is equal to the total one. This flux cannot be arbitrarily low, since the -term in the Lagrangian (1) is introduced after accounting for the proper quantization of the magnetic flux Hoyos (). If we take into account the quantum of magnetic flux, , then . In this case we would be dealing with effective electric charges which are fractionary ones, carried by quasiparticles.

## Iii Polarized Vacuum Energy on the TI Boundary

It is known that dynamical effects arise from the variations in the vacuum expectation value of observables associated with a quantum field, imposed by external boundary conditions or by nontrivial topologies of the physical space as well as by interactions with an external field generated from a particular source. The difference between the vacuum energies of these configurations is named Casimir energy, in the two first cases, and polarized vacuum energy, in the last one. In this section, we will consider the configuration of a static and spatially homogeneous magnetic field acting perpendicularly on the surface of a -type TI with TR symmetry in order to calculate the polarized vacuum energy of a spin 1/2 field situated on this edge, in both the weak-field approximation and general case.

### iii.1 Weak-field approximation

In the weak-field regimen, the energy gap opening caused by the incident magnetic field is neglected by considering the dispersion relation associated with the surface states of a TI around the Dirac points, given by Qien ()

(6) |

where is the Fermi velocity. Thus, we will deal with the polarized vacuum energy of a massless fermion in 2D space; the ordinary Lagrangian density of an electromagnetic field without sources plus a spinor field with charge interacting with the field is, in covariant language, given by

(7) |

where are the Dirac’s matrices (in 1+2 dimensions) and . The gauge choice is

(8) |

in cylindrical coordinates , compatible with a constant magnetic field pointing in the direction. The polarized vacuum energy density, , related to the spin 1/2 field in 1+2 dimensions for this particular gauge and for massive fermions, was found in Visser () by different methods, and it is given by

(9) |

in natural units (), with being the Hurwitz zeta function Milton (), which is the ordinary Riemann zeta function when the second entry vanishes. The last term in the expression above is the vacuum energy density of the spinor field with the magnetic field switched off. It is worth pointing out that this polarized vacuum energy did not need to pass by the renormalization of some parameters, since it is naturally finite, at least in 1+2 dimensions. We also remark that such energy does not generate any physical influence because it does not depend on a geometrical parameter.

By considering now the surface of a -type TI with TR symmetry, we have seen that one must add the axion-like term to the Lagrangian (7), given in covariant notation by

(10) |

where is the four-dimensional completely antisymmetric tensor. Therefore, the modified electromagnetic sector of the Lagrangian allows us to take into account the electric charge derived in Eq. (5), since this charge will correct the vacuum energy found in the literature, given by Eq. (9). Such a charge is responsible for the dynamical effect that we are looking for and it does not occur when ordinary charges alone are considered. Thus this topological charge must be taken into account in the coupling terms, and it will be assumed to possess fermionic nature because it is just a correction to the induced (ordinary) charges. Beside this, we must have , provided that we are dealing with carriers on a graphenelike conducting sheet that forms the TI surface, in which the linear dispersion for massless fermions is obeyed, at least around the Dirac’s points Castro (), as we have pointed out. We also must change , where this latter is the Fermi velocity. Joining all this and following Visser (), we arrive at

(11) | |||||

for the same gauge (8), after restoring the constants and .

The above result can also take into account the fact that the considered surface forms an interface between two TIs. In this case, the more general expression for the axionlike parameter is , Marcel (), and the expression for the polarized vacuum energy density becomes

(12) |

yielding a quantization of the polarized vacuum energy density. Notice that , too.

If we consider that the magnetic flux is given in a function of a linear dimension of the TI boundary, like , and if this flux is cylindrical with radius , then the polarized vacuum energy density depends on this length taken on the TI surface, which yields a density of radial force given by

(13) |

This Casimir force density is tangential to the TI surface, and can be seen as a radial shear stress that arises due to the quantum vacuum influence. At first instance, through the signal of this force we deduce that it tends to shrink the surface.

### iii.2 General case

Because a sufficiently strong magnetic field opens an interband gap on the boundary, Dirac fermions acquire an effective mass , which can be expressed as a function of the external magnetic field in the form

(14) |

with being the magnetic momentum of the carrier Shen (). Thus we must use Eq. (9) and adopt the same procedure made in the previous section; i.e., we must take into account the effective electric charge in order to get

(15) |

where, once more, . Explicitly, in terms of the geometric parameter , we have that

(16) |

where .

In Fig. 1 we depict the behavior of the polarized vacuum energy (designed here as Casimir energy) as a function of the radius of the region in which the magnetic field acts (i.e., the flux “tube” radius), taking all the non-numeric parameters equal to unity. We remark the monotonic energy increasing with this geometric parameter vanishes at for any and . For large values of , the polarized vacuum energy increases almost quadratically with this latter. We can still notice the cubic dependence of that energy on the magnetic field, such that it also grows monotonically with the field, for .

From Eq. (16) we calculate the radial force, which is given by

(17) | |||||

It is possible to show that this force diminishes with the radius for any positive value of and , becoming always more negative, which yields a growing shear stress that tends to shrink the surface.

## Iv Concluding Remarks

In this paper we have investigated the polarized vacuum on the surface of a class of TIs, labeled by the topological index and exhibiting time reversal (TR) symmetry, at zero temperature. Initially it was discussed how the Maxwell’s laws with source terms are modified in such systems, leading to the correction of the induced electric charge distribution by virtue of the presence of an axionlike pseudoscalar factor and when one considers both uniform magnetic and electric fields perpendicular to the boundary under analysis. Since the electric field does not influence the Landau levels on the TI surface Tong () and consequently the respective polarized vacuum energy, it was only indirectly taken into account in the electric charge induced by the field.

Thus, we have obtained an effective charge which is a correction to the induced ordinary charges due to the topological properties of the system and calculated the modifications in the polarized vacuum energy of a spin 1/2 field confined on one of the TI conducting edges, in both the weak-field approximation and general case. For the first case, we have also obtained the polarized vacuum energy density for the situation in which the surface forms an interface between two TIs, and found an expression that points out the quantization of this energy.

We have noticed that, differently from the ordinary polarized vacuum energy of 2-D fermions presented in the literature, the result found by us depends on a length parameter defined on the boundary, which yields a tangential force (in fact, a radial shear stress) associated to it. The negative signal of the found force indicates a tendency for it to shrink the surface. It is worth pointing out that the ordinary charges induced by the electric field do not generate alone any dynamical effects because they do not carry information on the topological properties of the system and thus the polarized vacuum energy is independent of its geometrical parameters Visser (); Beneventano ().

## Acknowledgements

M. O. Tahim and G. D. Saraiva would like to thank CNPq for partial financial support.

## References

- (1) M. Z. Hasan and C. L. Kane, Rev.Mod.Phys. 82, 3045 (2010).
- (2) A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. Ludwig: Phys. Rev. B 78, 195125 (2008).
- (3) H. B. G. Casimir, Proc. K. Ned. Akad. Wetensh. B 51, 793 (1948).
- (4) M. Bordag et al., Advances in the Casimir Effect (Oxford University Press, New York, 2009).
- (5) Casimir Physics, Series: Lecture Notes in Physics, Vol. 834, edited by D. Dalvit, P. Milonni, D. Roberts and F. da Rosa (Springer Heidelberg, 2011).
- (6) C. R. Muniz, V. B. Bezerra and M. S. Cunha, Phys. Rev. D 88, 104035 (2013).
- (7) V. B. Bezerra, H. F. Mota, and C. R. Muniz, Phys. Rev. D 89, 024015 (2014).
- (8) A. G. Grushin and A. Cortijo, Phys. Rev. Lett. 106, 020403 (2011).
- (9) P. Rodriguez-Lopez and A. G. Grushin, Phys. Rev. Lett. 112, 056804 (2014).
- (10) A. G. Grushin, P. Rodriguez-Lopez, and A. Cortijo, Phys. Rev. B 84, 045119 (2011).
- (11) F. Wilczek, Phys. Rev. Lett. 58, 1799 (1987).
- (12) Y. Ando, J. Phys. Soc. Jpn, 82, No. 10, 102001 (2013).
- (13) R. Li, J. Wang, X.-L. Qi, and S.-C. Zhang, Nat.Phys. 6, 284 (2010).
- (14) S. K. Blau, M. Visser, and A. Wipf, Int. J. Mod. Phys. A6, 5409 (1991).
- (15) G. Raffelt and L. Stodolsky, Phys. Rev. D 37, No. 5, 1237 (1988).
- (16) J. Redondo and A. Ringwald, Contemp. Phys. 52, 211 (2011).
- (17) X. L. Qi, R. Li, J. Zang, and S. C. Zhang, Science 323, 1184 (2009).
- (18) A. M. Essin, J. E. Moore, and David Vanderbilt, Phys. Rev. Lett. 102, 146805 (2009).
- (19) P. Cheng et al., Phys. Rev. Lett 105, 076801 (2010).
- (20) C. G. Beneventano and E. M. Santangelo, J.Phys. A 39, 7457 (2006).
- (21) C. Hoyos-Badajoz, K. Jensen, and A. Karch, Phys.Rev. D 82, 086001 (2010).
- (22) D. Hsieh et al., Nature (London) 452, 970 (2008).
- (23) M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions (Dover Publications, 1964).
- (24) A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 (2009).
- (25) X.-L. Qi, in Topological Insulators, edited by M. Franz and L. Molenkamp (Elsevier, New York, 2013).
- (26) S. Q. Shen, Topological Insulators: Dirac Equation in Condensed Matters (Springer-Verlag, Berlin-Heidelberg, 2012).