Electric polarization in magnetic topological nodal semimetal thin films

We theoretically study the electric polarization in magnetic topological nodal semimetal thin films. In magnetically doped topological insulators, topological nodal semimetal phases emerge once the exchange coupling overcomes the band gap. Changing the magnetization direction, nodal structure is modulated and the system becomes topological nodal point or line semimetals. We find that nodal line semimetals are characterized by non-linear electric polarization, which is not observed in nodal point semimetals. The non-linear response originates from the existence of the surface states. Screening effect is self consistently included within a mean field approximation and the non-linear electric polarization is observed even in the presence of screening effect.

We theoretically study the electric polarization in magnetic topological nodal semimetal thin films. In magnetically doped topological insulators, topological nodal semimetal phases emerge once the exchange coupling overcomes the band gap. Changing the magnetization direction, nodal structure is modulated and the system becomes topological nodal point or line semimetals. We find that nodal line semimetals are characterized by non-linear electric polarization, which is not observed in nodal point semimetals. The non-linear response originates from the existence of the surface states. Screening effect is self consistently included within a mean field approximation and the non-linear electric polarization is observed even in the presence of screening effect.

I. INTRODUCTION
Topological semimetals are characterized by topological numbers and topological responses [1][2][3][4]. One of famous examples is a magnetic Weyl semimetal which is characterized by a Chern number defined on a closed surface enclosing a nodal point called Weyl point. The Weyl points behave as a sink or source of Berry curvature, which leads to intrinsic anomalous Hall effect [5]. Magnetic Weyl semimetals exhibit a semiquantized anomalous Hall effect as a topological response [1,6,7]. The intrinsic anomalous Hall effect in a magnetic Weyl semimetal is observed in recent experiments [8,9]. On the other hand, nodal line semimetals are characterized by a Zak phase [10]. The Zak phase is closely related to the electric polarization [11][12][13], so that we expect non-trivial electric polarization in the nodal line semimetals [3,14]. However, the electric polarization is well defined only in insulators and the applicability of the theory of polarization based on the Zak phase is ambiguous in the nodal semimetals.
In this work, we theoretically study electric polarization in magnetic topological nodal semimetal thin films. We consider magnetically doped topological insulators [15][16][17][18][19]. We assume ferromagnetic ordering of magnetic moments and introduce the exchange coupling between the band electrons and the magnetization. Once the exchange coupling overcomes the band gap, the system is brought into the magnetic topological semimetal phases [1,20]. The topological property of the band structure depends on the direction of the magnetization and the system becomes the nodal point and line semimetals by changing the direction of the magnetization [2,21,22]. There are other candidate materials of the magnetic topological semimetals [23,24], where the topology of the electronic structure depends on the magnetic structure. We numerically calculate the band structure of the topological semimetal thin films in the presence of the external potential. The screening effect is self-consistently included. We find that nodal line semimetals exhibit nonlinear electric polarization, which originates from the existence of the surface states. In the nodal point semimetals, the non-linear electric polarization is not observed. Therefore, the nodal line semimetal phase is characterized by the non-linear electric polarization.
This work is organized as follows. In Section II, we introduce model Hamiltonian which describes electronic structures of topological insulators and exchange coupling between band electrons and local magnetic moments. The screened external potential is also added. We explain the band structure of bulk systems. In Section III, we formulate a theoretical procedure to calculate a screened potential and a self-consistently determined electric polarization. In Section IV, we show numerically calculated band structures of the bulk and slab systems. We discuss the polarization of each state in the presence of the external potential and the numerically calculated electric polarization. Conclusions are given in Section V.

II. MODEL HAMILTONIAN
We start with the lattice model for three dimensional topological insulators (Bi 2 Se 3 , Bi 2 Te 3 , and Sb 2 Te 3 ) [25,26], where t is a hopping parameter, a is a lattice constant, and m k is a mass term, σ and τ are the Pauli matrices acting on the real spin and the pseudo spin (orbital) degrees of freedoms. The exchange coupling between the band electrons and the local magnetic moments is written as [27] where J 0 , J 3 are exchange coupling constants andM is a unit vector representing the direction of magnetization. In this work,M is set on the x-z plane as we see in Figure.1 (a). The first term is an usual exchange coupling and the second term originates from nonequality of the exchange coupling constants between two orbitals considered here, i.e., p-orbitals of (Bi,Sb) and (Te,Se). In magnetic topological insulators [16,18,19], Cr atoms are substituted for Bi or Sb atoms. As a result, there is the nonequality of exchange coupling between p-orbitals of (Bi,Sb) and (Te,Se). The total Hamiltonian is given as where U (r) is the screened external potential. In the following calculation, we self-consistently include the screening effect and the detail of the numerical procedure is given in the next section. The topological phase diagram of bulk system with no-external electric field is shown in Figure (1) (b) θ = 0 and (c) θ = π/2. The system becomes topological semimetal phase when the exchange coupling overcomes the energy gap. In the present system, both of topological nodal line and point semimetal phases emerge. At θ = 0, the system becomes a nodal line semimetal in J 0 < J 3 and a nodal point semimetal in J 0 > J 3 . At θ = π/2, on the other hand, the system becomes a nodal line semimetal in J 0 > J 3 and a nodal point semimetal in J 0 < J 3 . This means that the topological property of the electronic structure can be modulated by manipulating the magnetization direction. The topologically non-trivial band structure in nodal line semimetals is characterized by the Zak phase where the summation is over the occupied states and |u nk is an eigenstate of the bulk Hamiltonian. In Section IV, we discuss the relation between the Zak phase and existence of surface states. With a surface boundary for x direction as represented in Figure 1 (a), the wave function is written as where ψ(r) and φ(x) are four component spinors, and L y , L z are the system sizes along y and z directions. The Schrödinger equation becomes We numerically diagonalize the above Hamiltonian and derive the eigenstates and the energy bands.

III. ELECTRIC POLARIZATION INDUCED BY THE SCREENED POTENTIAL
We calculate electronic structure in the presence of the external electric field. The screening effect is self-consistently included within a mean field approximation by simultaneously solving the Schrödinger equation and the Poisson's equation. Our numerical procedure is as follows. We set an initial value for the external potential as where N x a is the width of the slab. We derive eigenstates in the presence of the external potential. We consider the slab as multilayers of zero-thickness planes with charge density ρ 2D (x). Using the derived eigenstates, ρ 2D (x) is calculated as where f (ε nk ) is the Fermi distribution function. The electric field is calculated as where κ is the permitivity and n x = ⌊ x a ⌋. ∆ρ 2D (ja) is the induced charge density, ∆ρ 2D (ja) = ρ 2D (ja) − ρ 2D 0 (ja), where ρ 2D 0 (ja) is the charge density in the absence of the external potential. Here, we assume charge neutrality of the slab, E(x) = E 0 (x < 0, L x < x) or j ∆ρ 2D (ja) = 0. Integrating the above electric field, the screened potential is calculated as where we set a boundary condition for the screened potential so that the average of the screened potential is zero, j U (ja)=0. We perform numerical iteration so that the screened potential converges. The strength of the screening is characterized by a dimensionless parameter α defined as Using the eigenstates with the screened potential, we calculate the electric polarization of each eigenstate, The total electric polarization from the occupied eigenstates is written as In the following section, we numerically diagonalize Equation (7) and self-consistently determine the electric polarization.

IV. NUMERICAL RESULTS
As represented in Figure 1 (a), we consider the open boundary condition for the x direction, and the periodic boundary condition for the y and z directions. We set the parameters as m 2 /t = 1, m 0 /t = −0.5, J 0 /t = 1, and N x = 30. First, we consider a nodal line semimetal phase, θ = π/2. Figure 2 shows the band structure of nodal line semimetals for several parameters. At J 3 = 0, there is a nodal line on the k x = 0 plane in the bulk system. In Figure 3, we see that the Zak phase is quantized as π when the integration path go through the area enclosed by the nodal line. Otherwise, the Zak phase becomes zero. In the slab system, the wave number k x is discretized and we obtain subband structures. In Figure 2 (b), we show the lowest subbands. There are degenerate flat bands in the wave numbers (k y , k z ) enclosed by the nodal line where the Zak phase is quantized as π. The flat band states are localized on the left and right surfaces. In the presence of the external potential, the degenerate flat bands are gapped as shown in Figure 2 (c). The colored bar represents the polarization of each state. We see that the surface states localized on the left surface are occupied and on the right surface are empty. It means that the surface states are completely polarized and the surface electric polarization occurs. The occupied bulk states are polarized to the opposite direction so that the surface electric polarization is screened. At finite J 3 , the qualitative behavior is the same as J 3 = 0 but the band structure is slightly modified. The nodal line is not bounded on ε = 0 and the surface bands are bended. As a result, a larger external potential is required in order to achieve complete surface polarization compared with the flat surface bands i.e., J 3 = 0. Figure 4 shows the band structure of nodal point (Weyl) semimetal phase, θ = 2π/5. When the magnetization has a finite component along z direction, the nodal line is gapped except for a pair of nodal points and the system becomes nodal point (Weyl) semimetals. In the slab system, there are chiral surface states localized on the left or right surfaces and they are seamlessly connected to the bulk states. The velocity of the chiral surface states increases with the decrease of θ and it becomes maximum at θ = 0. In the presence of the external potential, the left surface states are pushed down and the right surface states are pushed up. Consequently, the surface electric polarization occurs but it is weak compared with that in the nodal line semimetals. In the nodal line semimetals, we observe the surface electric polarization which is proportional to the number of surface states in the sufficiently strong external potential. On the other hand, the required external potential to achieve the surface electric polarization which is proportional to the number of the surface states in the nodal point semimetals is much larger than that in the nodal line semimetals.
Using the self-consistently derived electronic states, we calculate the electric polarization in the slab system. Figure  5 shows the electric polarization as a function of the external electric field for several parameters. In the nodal line 3: (a) The black closed curve represents nodal line in the Brillouin zone. Two paths C1 and C2 are depicted by red arrows and give the Zak phase π and 0, respectively. (b) The Zak phase calculated at fixed (ky, kz). The Zak phase is quantized as π in the area enclosed by the nodal line projectoed on ky-kz plane and it is zero in the other area.
Left Right semimetal, Figure 5 (a), we see that the system exhibits non-linear electric polarization. In the absence of the screening effect (α = 0), the electric polarization abruptly changes at E 0 = 0. In the bulk limit, the electric polarization is induced by infinitesimal electric field and discontinuously changes at E 0 = 0. However, in our numerical calculation, the electric polarization continuously changes at E 0 = 0 because of the finite size effect. In the presence of the screening effect, the step structure is broadened and the electric polarization vanishes at E 0 = 0. However, the kink structure of the electric polarization is preserved. Figure 5 (b) shows the electric polarization for finite chemical potential µ/t = 0.03. Even in the finite chemical potential, the kink structure is observed, though the kink appears at a finite electric field. This kink structure is broadened by the screening effect. In finite J 3 term, Figure 5 (c), the surface bands are bended and the situation is similar to the finite chemical potential case, Figure 5 (b), so that the qualitative behavior is the same. In the nodal point semimetal, Figure 5  to the external electric field and there is no kink structure. Therefore, the kink structure is characteristic to the nodal line semimetals and is closely related to the topological property of the system. Figure 6 (a) shows the electric polarization as a function of the angle θ. There are peak structure of the electric polarization at θ = π/2 because the system becomes the nodal line semimetal at this angle and exhibits large electric polarization originating from the surface states. With the increase of J 3 , the peak value becomes smaller. This is because the surface bands are bended and the kink structure is broadened. Figure 6 (b) shows variation of the total energy defined as where ε nk (θ) is an eigenenergy of the slab system with the magnetization angle θ. The variation of the total energy also becomes maximum at θ = π/2 and decreases away from θ = π/2. Unlike the electric polarization, ∆ε(θ) weakly depends on J 3 and varies slowly compared with the electric polarization. This is because the electric polarization is almost determined by the property of the surface states but the variation of the total energy is almost determined by the bulk band structure. Figure 6 (b) means that the magnetization easy axis is parallel to the surface.

V. CONCLUSIONS
We investigated the electronic structure in the magnetically doped topological insulator thin films. We found that the system becomes topological nodal semimetal phases when the exchange coupling exceeds the band gap and the topological property of the band structure is modulated by changing the magnetization direction. We calculated the electric polarization in the topological nodal semimetal thin films. We self-consistently include the screening effect and find that the nodal line semimetals exhibit non-linear electric polarization even in a finite chemical potential. This is a topological response and the nodal line semimetals are characterized by the non-linear electric polarization, which is not observed in the nodal point (Weyl) semimetals.