Preprint
Article

Time Reversal Symmetry Broken Floquet State of SMB_6 - A Strong Topological Insulator

This version is not peer-reviewed.

Submitted:

15 August 2023

Posted:

16 August 2023

You are already at the latest version

Abstract
We show that the compound samarium hexaboride is a strong topological insulator using the eigenvalues of the space inversion operator in the low-energy limit of the periodic Anderson model. Additionally, we assume the presence of the ferromagnetic exchange interaction (M). A Dirac cone like feature in the surface state energy spectra is observed for M equals zero in a certain parameter range. For M not equal to zero, there is no Kramers degeneracy. We have been able to show that this phase corresponds to the quantum anomalous Hall state by calculating Berry curvature and the Chern number. Using Floquet theory, we further show that the access to a novel state with broken time reversal symmetry is possible due to the normal incidence of circularly polarized optical field on the surface of the compound despite M being absent.
Keywords: 
Periodic Anderson model; Dirac cone; Kramers degeneracy; Quantum anomalous Hall state; Chern number; Floquet theory
Subject: 
Physical Sciences  -   Condensed Matter Physics

1. Introduction

This communication is on a topical issue of the mixed valence compound SmB6.[1,2,3,4,5]—a narrow gap topological Kondo insulator (TKI) which, with a high-temperature metallic phase, transforms into a paramagnetic charge insulator below 45 K. It has been suggested [6,7], as well as it is increasingly apparent [8,9,10] during the past several years, that SmB6 is a non-trivial topological insulator. This has generated great deal of excitement in the condensed matter physics community and it still remains a matter of avid debate [6]. Despite the supporting evidence for the TKI scenario [1,2,3,4,5,6,7,8,9,10], there is no prognosis regarding the nature of the bulk and surface states of SmB6 [11,12,13]. In this paper our primary aim is to resolve this issue. We consider an extended periodic Anderson model (EPAM) [14] for the compound SmB6 for this purpose. We introduce here the exchange interaction (M) assuming the presence of the ferromagnetic magnetic impurities in the system. The slave boson (SB) mean-field-theoretic Anderson model of Legner [15] refers to a simple cubic lattice with one spin-degenerate orbital per lattice site each for d   and   f electrons. We consider the low-energy version of this model together with the exchange interaction. Our minimalistic Hamiltonian, based on the slave boson (SB)mean field theory of ref. [15], captures essential physics of TKI in the presence of the coulomb repulsion U f ( >> t d 1 ) between f electrons on the same site, and the spin-orbit hybridization V. The parameter V is the harbinger of a topological dispensation. The terms ( t d 1 ,   t f 1 ) are the nearest neighbor hopping parameters for d   and   f electrons. There are three other parameters (b, λ, ξ ) of our theory [14]. While the term ξ enforces the fact that there are equal number of d and f fermions, the parameter b represents a c-number slave-boson field. We note that the constraint U f   >> t d 1 imposes a non-holonomic constraint, viz. the exclusion of the double occupancy. The SB-protocol provides a platform to reformulate this nonholonomic constraint into a holonomic constraint that can be implemented with the Lagrange multiplier λ. We found that λ = 6 t f 1 + 6 b 2 t f 1 ,   and ξ = 3 t d 1 + 3 t f 1 . The admissible value of b 2 is 1. Since the method to obtain them is explained clearly in ref. [14] we will not reproduce the same here. We observe a Dirac cone like feature at k = 0   momentum in the surface state energy spectra for M = 0 [16], upon writing the SB Hamiltonian in the Dirac basis similar to the Bernevig–Hughes–Zhang (BHZ) model [17].  We obtain the Z2 invariant (Z2 = 1 ) using the eigenvalues of the space inversion operator in the Fu-Kane framework [18]. This is the conclusive evidence of the fact that SmB6 is a strong topological insulator for M = 0. By calculating Berry curvature and the Chern number we have been able to show that M 0 corresponds to the quantum anomalous Hall state. It must be mentioned here that our low-energy model has been found adequate enough to capture important features of SmB6, although the ground state of the compound SmB6 has been shown to be is a quartet state [1].
The exotic Floquet topological phases [19,20,21,22,23,24,25,26] with a high tunability could be realized using the circularly polarized optical field (CPOF). Analogous to the Bloch theory, here one can transform the time-dependent Hamiltonian problem to a time-independent one using the Floquet’s theorem [27,28,29,30]. The combinations of the Floquet theory with dynamical mean field theory [31], and slave boson protocol [32] were also formulated for strongly correlated systems. Our approach is in acquiescence to the latter. We use the Floquet theory[27,28,29,30,33] in section 3 in the high-frequency limit to examine the Co3 Sn2 S2 thin film system. Interestingly, the incidence of CPOF leads to the time reversal symmetry (TRS) broken phase despite M = 0. The conclusive evidence of this phase being quantum spin Hall (QSH) phase is presented by calculating the spin Chern number [34,35].
The paper is organized as follows: In section 2, we calculate Z2 invariant using the eigenvalues of the space inversion operator with in the Fu-Kane framework [18]. In section 3, we are able to show the emergence of a novel phase with broken-TRS by the normal incidence of tunable CPOF despite M = 0. For this purpose, we make use of the Floquet theory in the high-frequency limit to investigate the system. The paper ends with discussion and brief concluding remarks in section 4.

2. Surface State Hamiltonian in SB formalism and the Z2 invariant

A. Surface State Hamiltonian

We treat the model Hamiltonian ( H P A M b , λ , ξ   ) in Eq.(1) of the ref.[14] in the low-energy limit below. All energies in our calculation/graphical representation below are expressed in units of the first neighbor hopping t d 1 for d-electrons as this corresponds to the kinetic energy of these itinerant electrons and therefore the most dominant. In the above limit, the following replacements are necessary : sin a j k j   a j k j + O a j 3 k j 3 , cos a j k j ( 1 ( 1 2 a j 2 k j 2 ) ) where k j (j = ( x, y, z)) are momentum components, and a j is the lattice constant along j direction. Furthermore, in order to obtain surface state Hamiltonian, we make the replacement a k z i a z   and look for states localized within the surface z = 0 of the form e x p i χ z   . Further more , we seek such a value of the unknown wave number χ ( χ = i q ,   q > 0 ) for which the exponential e x p a q z a 1   for   z > 0 . For example, if we assume a q   ~ 1 the exponential e x p a q z a ~ exp 10 , i,e vanishingly small ,   for   z   ~ 50 n m given that SmB6 cubic crystal structure with lattice constant a = 0.413 nm. Therefore, χ = i q ,   q > 0 ensures a decaying term for z > 0 in the surface states. Upon including the exchange coupling M , we find that the surface state Hamiltonian h surface k ,   q , μ , b , M is
h surface = Ϧ + L 1   L 1 Ϧ , L 1 = 0   i A 1 a q i A 1 a q 0   ,   Ϧ + = Ϧ k x ,   k y ,   q , μ , b , M =   ϵ k ,   q , μ , b τ 0 + n k x ,   k y ,   q , b . σ + M σ z   , Ϧ = Ϧ * k x , k y   q , μ , b , M
The two blocks Ϧ + , Ϧ ,characterized by the pseudo-spin indices ( +, −), are related to each other by TRS for   M = 0 , and   n k x ,   k y ,   q , b = ( A 1 a k x ,   A 1 a   k y ,   ϑ k x ,   k y ,   q , b ) . The other parameters/functions in (1) are
A 1 =   2 V b ,   k = k x ,   k y ,   k 2 = ( k x 2 + k y 2 ) ,
ϵ k ,   q , μ , b = ϵ 0 μ , b D 1 ( b ) a 2 q 2 + D 1 ( b ) a 2 k 2 + O a 4 q 4 + O a 4 k 4
ϵ 0 μ ,   b = 𝜇 + [ b 2 2 ϵ f 3 t d 1 3 t f 1 6 t d 2 6 t f 2 b 2 4 t d 3 4 t f 3   b 2 ] ,
D 1 b = [   t d 1 + b 2 t f 1 2 + 2 ( t d 2 + b 2 t f 2 ) + 2 ( t d 3 + b 2 t f 3 ) ] ,
ϑ k , q , b   = ϑ 0 ( b ) B 1 ( b ) a   2 q 2   + B 1 ( b ) a 2 k 2 + O a 4 q 4 + O a 4 k 4 ,  
ϑ 0 b = [ b 2 ϵ f 2 6 t d 2 + 6 t f 2 b 2 4 t d 3 + 4 t f 3   b 2 ] ,
B 1 b = [   t d 1 b 2 t f 1 2 + 2 ( t d 2 b 2 t f 2 ) + 2 ( t d 3 b 2 t f 3 ) ] .
Here the terms t d 1   ,   t f 1 , ( t d 2 ,     t f 2 ) , and t d 3   ,   t f 3 , respectively, are the NN, NNN, and NNNN hopping parameters [14] for d   and   f electrons, ϵ f is the onsite energy of the f electrons, μ is the chemical potential of fermion number, σ x , y , z are the Pauli matrices, and   σ 0 is the 2×2 identity matrix.   On a quick side note , we observe   that h surface k ,   q , μ , b , M = 0 corresponds to Qi-Wu-Zhang (QWZ) model[35]. As shown by these authors, the situation corresponds to the QSH state, for the spin Hall conductance of Ϧ ( k x ,   k y ,   q , μ , b , M = 0 and Ϧ * k x , k y   q , μ , b , M = 0   are not zero but the net Hall conductance of the system described by the model is zero.
For calculating the Z2 invariant, the eigenvalues of the parity operator needs to be obtained. The objective could be accomplished with relative ease if the Hamiltonian in (1) is written down in the Dirac basis similar to the Bernevig–Hughes–Zhang (BHZ) model [17] presented over a decade and half ago for quantum wells. We now write the Hamiltonian in (1) in Dirac basis similar to the BHZ model:
h surface B H Z k   q , μ , b , M = ϵ d + ϵ c 2 I 4 × 4 + a = 0,12,3 , 5 d a k , q , μ , b γ a + M σ 0   τ z  
where d 0 = ϵ d ϵ c 2 , ϵ d = ( ϵ k , q , μ , b + ϑ k , q , b ) , ϵ c = ( ϵ k , q , μ , b ϑ k , q , b ) , d 1 = i A 1 a k y , d 2 = i A 1 a k x ,   d 3 = A 1 a χ ,   and d 5 = 0 . The Dirac matrices ( γ 0 , γ 1 , γ 2 , γ 3 , γ 5 ) in con-travariant notations are γ 0 =   σ z I 2 × 2 , γ j = i σ y τ j , j = 1,2 , 3 , and γ 5 = i γ 0 γ 1 γ 2 γ 3 . The Pauli matrices σ and τ are acting in the space of bands that give rise to Kramers degeneracy (see Figure 1).

B. Surface state spectrum

The eigenvalues j = j s , σ , k , b ,   M   of the matrix (3) are given by a quartic and, therefore, we use Ferrari’s solution of a quartic equation (see Appendix A). We have plotted the surface state energy spectra (SSES) given by Eq. (A.2) as function of the dimensionless wave vector a k in Figure 1(a) and Figure 1(b). Since the conduction bands are partially empty, the surface state will be metallic. These figures correspond to unbroken TRS (M = 0 ) . It will be shown below, calculating the Z2 invariant using the eigenvalues of the space inversion operator, that the figures correspond to QSH phase. There should be a surface Dirac cone (or at least a Kramers degeneracy) at k=0, as in ref. [16]. We indeed observe a Dirac cone like feature here (see Figure 1(a) and Figure 1(b)) in a certain parameter-window. The numerical values of the parameters used in the plots are t d 1 = 1, t f 1 = 0.8, t d 2 = 0.01, t f 2 = 0.01, t d 3 = 0.001, t f 3 = 0.001, ϵ f =   0.02, V = 0.1, b =0.98,μ = 0 ,       M = 0 ,   0.3 and   U f t d 1 . In both the figures, t f 1 is negative and, therefore, the figures correspond to the insulating bulk. It may be noted that one needs t f 1 ~   t d 1 to access the Dirac-cone feature. The Dirac-cone feature of SSES agrees with several experimental observations reported earlier, such as those by scanning-tunneling microscopy [36,37], angle-resolved photoemission spectroscopy (ARPES)[38,39] and the circular dichroism ARPES[40], and so on. Here, the red curve corresponds to the spin-up valence band 3 s = 1 , σ = + 1 , k , b ,   and the green curve to spin-down conduction band 2 s = + 1 , σ = 1 , k , b . The curves display the band-inversion close to the Fermi energy represented by the horizontal solid line. In Figure 1(c), though M = 0   ,   μ   0 .   We observe that the states corresponding to momenta ak = ( ± 2,0 )   or ( 0, ± 2 ) in Figure 1(a) are degenerate. Furthermore, they satisfy the condition a k + a G =     a k where a G is a reciprocal   lattice vector . For example, for ak = ( ± 2,0 )   and a G = 4,0 . Of course, there are other possibilities too, for example ak = ( ± 2 ,   ± 2 ) . These possibilities we are not taking into account for the simple reason that they do not satisfy the condition k + G =   k . In Figure 1(d), M 0 and therefore TRS is broken. There is no Kramers degeneracy as could be seen in this figure. The Figure 1(d) corresponds to QAH as is shown below.
Figure 1. The plots of surface state energy spectrum given by Equation (14) as a function of the dimensionless wave vector a k . The numerical values of the parameters used in the plots are t d 1 = 1, t f 1 = ( 0.8, 0.5), t d 2 = 0.01, t f   2 = 0.01, t d 3 = 0.001, t f 3 = 0.001, ϵ f =   0.02, V = 0.1, b =0.98,μ = (0, 0.5 ) ,       M = ( 0,0.5 )   and   U f t d 1 . The horizontal solid line represents the Fermi energy. Since the conduction bands are partially empty, the surface state will be metallic in all the cases. In (a)-(c), the system is TR symmetric. The energy bands of the system come in Kramers pairs. We observe a Dirac cone like feature in Figures (a) and (b). In Figure (d), TR symmetry is lacking and therefore no Kramers pair is possible.
Figure 1. The plots of surface state energy spectrum given by Equation (14) as a function of the dimensionless wave vector a k . The numerical values of the parameters used in the plots are t d 1 = 1, t f 1 = ( 0.8, 0.5), t d 2 = 0.01, t f   2 = 0.01, t d 3 = 0.001, t f 3 = 0.001, ϵ f =   0.02, V = 0.1, b =0.98,μ = (0, 0.5 ) ,       M = ( 0,0.5 )   and   U f t d 1 . The horizontal solid line represents the Fermi energy. Since the conduction bands are partially empty, the surface state will be metallic in all the cases. In (a)-(c), the system is TR symmetric. The energy bands of the system come in Kramers pairs. We observe a Dirac cone like feature in Figures (a) and (b). In Figure (d), TR symmetry is lacking and therefore no Kramers pair is possible.
Preprints 82469 g001
In Figure 2 we have shown the contour plots of the Berry-curvature(BC) in the z-direction for M 0 . The numerical values of the parameters used in the plots are t d 1 = 1, t d 2 = 0.01, t f   2 = 0.01, t d 3 = 0.001, t f 3 = 0.001, ϵ f =   0.02, V = 0.1, b = 0.98, μ = 0 ,     and   U f t d 1 . The parameter t f 1 = 0.8 in Figure (a) whereas t f 1 = 0.6 in Figure (b). Upon integrating BC on a k-mesh-grid of the Brillouin zone (BZ), we calculate the intrinsic anomalous Hall conductivity σ A H ( AHC). This yields the Chern number (C). A brief outline of the procedure followed is given below in Appendix D. We find that AHC is σ A H = 1.04 71 e 2   in the formercase ( C = 1.0471   is close to the integer value 1), while in the latter case it is   σ A H = 1.6607 ( e 2 ) . Furthermore, we found that while for t f 1 = 0.4, C = 3.3906, for t f 1 =     0.2, C = 3.2181. Our conjecture for not obtaining integer values of Chern number leans upon the following: The Hall conductivity σ A H   cannot be determined as such from the 2D Dirac model since Eq. (D.1) (see Appendix A) requires an integral over the whole BZ. The integral is outside the Dirac model’s range of validity. To circumvent the problem one may possibly choose a momentum space cut-off small compared to the size of BZ and large enough to capture nearly all the contributions to BC integral. This will be within the range of validity of the 2D Dirac model.
Figure 2. The contour plots of the Berry-curvature in the z-direction for M = 0.08   as a function of the dimensionless wave vector components a k x   and   a k y . The numerical values of the parameters used in the plots are t d 1 = 1, t d 2 = 0.01, t f   2 = 0.01, t d 3 = 0.001, t f 3 = 0.001, ϵ f =   0.02, V = 0.1, b =0.98, μ = 0 ,     and   U f t d 1 . The parameter t f 1 = 0.8 in Figure (a), whereas t f 1 = 0.6 in Figure (b).
Figure 2. The contour plots of the Berry-curvature in the z-direction for M = 0.08   as a function of the dimensionless wave vector components a k x   and   a k y . The numerical values of the parameters used in the plots are t d 1 = 1, t d 2 = 0.01, t f   2 = 0.01, t d 3 = 0.001, t f 3 = 0.001, ϵ f =   0.02, V = 0.1, b =0.98, μ = 0 ,     and   U f t d 1 . The parameter t f 1 = 0.8 in Figure (a), whereas t f 1 = 0.6 in Figure (b).
Preprints 82469 g002
Let us now note that the Fermi energy inside the gap intersects the surface state bands in the same BZ, in general, either an even or an odd pair number of times. If there are odd numbers of pair intersections, which guarantees the time reversal invariance, the surface state is topologically non-trivial (strong topological insulator). Furthermore, it is also evident that the number of TRIM pair involved in the surface-state crossings (SSC) is one (odd). However, when there are an even number of pair-surface-state crossings, the surface states are topologically trivial (weak TI or ordinary Bloch insulators that are topologically equivalent to the filled shell atomic insulator). The quantized topological numbers, the Kane–Mele index Z2 for QSH phase and the Chern number Ϲ for QAH phase, strongly support such topological states. The QSH band structures are characterized by the topological invariant ν 0 = 0 (Z2= +1) and ν 0 = 1 (Z2 = −1). The former corresponds to weak TI, while the latter to strong TI. In fact, materials with band structures with Z2= −1 are expected to exist in systems with strong spin-orbit coupling acting as an internal quantizing magnetic field on the electron system. The graphical representations in Figure 2 (except Figure 2(d)) indicate that the system considered here is a strong TI for M = 0. This is corroborated by analysis given below.

C. Z2 invariant

In a bid to to calculate the Z2 invariant using the eigenvalues of the space inversion operator, we note that the f- and d-states have different parities, the inversion symmetry(IS) operator in this band basis is constructed as Π = I 2 × 2 τ   z . The time reversal (TR) operator for a spin 1/2 particle is Θ =   I 2 × 2 τ y K .   The operator K   stands for the complex conjugation. Here the τ j are Pauli matrices on two-dimensional spin space. The Hamiltonian under consideration, for M = 0, preserves the time reversal (TRS)and inversion symmetries (IS). It can be easily shown that Θ ψ Θ φ =   φ ψ taking eigenstate of the z-component of the spin operator I 2 × 2 τ z as the basis. Also, Θ   γ 0 Θ −1 =   γ 0 , Θ   γ 5 Θ −1 =   γ 5 ,and Θ   γ j Θ −1 =   γ j , where j = 1,2 , and 3 . Similarly, Π   γ 0 Π −1 = γ 0 , Π   γ j Π −1 =   γ j ( j = 1,2 ) , and Π   γ k Π −1 =   γ k ( k = 3,5 ) .   Since only γ 0   a n d   γ 5 are even under time reversal and inversion (and d 5 = 0 ) , at a time reversal invariant momentum(TRIM) Ki where the system preserves both TR and IS , the Hamiltonian will have the form h surface B H Z k = K i   ,   q , μ , b , M = 0 = ϵ d + ϵ c 2 I 4 × 4 + d 0 k =   K i , q , μ , b γ 0 . The eigenvalues of γ 0 are ± 1 ( multiplicity 2). The corresponding eigenvectors are + = (1/ 2 ) ( 1 1 0 0)T and = (1/ 2 ) (0 0 1   1 )T. Here + γ 0 + = 1   and γ 0 =   1 . We obtain
+ h surface B H Z k = K i   ,   q , μ , b , M = 0 , W 0 + = ϵ d = E + ,   h surface B H Z k = K i   ,   q , μ , b , M = 0 ,   W 0 = ϵ c = E .
Here
d 0 = ϵ d ϵ c 2   ,   ϵ d = ϵ k ,   q , μ , b + ϑ k ,   q , b ,   ϵ c = ϵ k ,   q , μ , b ϑ k ,   q , b .
Obviously enough, if E < E + ,   the state is occupied and the parity of the state at TRIM Ki is 1 .   In the opposite case ( E > E + ) , the state + is occupied and the parity is + 1 . Therefore, the parity is given by ( s g n d 0 ) .
We now obtain the Z2 invariant simply by the parity eigenvalues at TRIMs. The surface states correspond to the eigenstates (or the Bloch states linked to the eigenvalues in j s , σ , k , b ,   M = 0 . These are presented in Appendix A. We consider a matrix representation of the time reversal (TR) operator Θ in the Bloch wave function basis. With α and β as the band indices we consider the representation is   ξ α β k   = ύ α k Θ ύ β k . This matrix relates the two Bloch states | ύ α / β k / k   and     via   | ύ α k =   β ξ α β * k Θ | ύ β k .   With the aid of this one can easily show that   ξ α β k is a unitary matrix ( ξ ξ = I ) . We also find that it has the property ξ α β k = ξ β α k . This implies that the matrix ξ α β K i at a TRIM becomes anti-symmetric, i.e. ξ α β K i =   ξ β α K i   0. Only when the bands α and β form a Kramers pair, such a non-zero ξ α β is obtained. Yet another which we need to consider is the Berry connection matrix defined as γ α β k = i ύ α k k ύ β k . In view of the results Θ ψ Θ φ =   φ ψ and | ύ α k =   β ξ α β * k Θ | ύ β k we arrive at the relation linking γ α β k and γ α β k :
γ k =   ξ k γ * k ξ k + i ξ k k ξ k .
Upon taking the trace we find tr γ k = t r ( γ * k ) + i t r ( ξ k k ξ k ) . Since γ β α =   γ α β * and ξ ξ = ( ξ ) ξ , upon replacing k by   k in the preceding equation one may write A = tr γ k = t r ( γ k ) + i t r ( ξ k k ξ k ) .We shall need this result below. The Berry curvature of tr γ k may be defined as Ω = curl A. Since the system preserves TRS and inversion symmetries (IS), one may select any gauge which renders A equal to zero. We now consider the anti-symmetric and unitary matrix   ζ α β k   = ύ α k Π Θ ύ β k (where Π 2 =1) to examine the consequence of setting A equal to zero. Since we find from ref.[18] that A = tr γ k = i 2 t r ζ k ζ = i 2 k t r l o g ζ = i k log det ζ ), it is clear that in order to make A = 0, one needs to adjust the phase of Bloch states | ύ α k such that Pf(ζ) = 1. Suppose now ρ (   K i trim ) = ± 1   are the eigenvalues of Π for band α at TRIM K i trim , one obtains the matrix
  ξ α β K i trim = ύ α K α trim Π Π Θ ύ β K i trim = ζ α β K i trim   ρ α (   K i trim ) .
Obviously enough, when ρα = ρβ, ξ α β K i =   ξ β α K i   0. Only when the bands α and β form a Kramers pair, such a non-zero ξ α β is obtained. It follows that if the bands α and β are the nth Kramers pair in the total of 2N bands, we may write ρ α = ρ β ρ 2 n . From Eq.(6), one can now see that
P f   ξ α β K i trim = P f   ζ α β K i trim   n = 1 N ρ 2 n   K i trim .  
Since   Pf ζ = 1 , in view of this result and Eq.(B.3) in the Appendix B, we find that that the Z2 invariant can be calculated simply by the parity eigenvalues ρ 2 n at TRIMs K i trim ,   that is 1 ν = i δ K trim i =   i n ρ 2 n K trim i . Upon getting back to Eq.(5) and taking into account the observations below this equation, the parity of the occupied state at K trim i , viz. δ K trim i , is given by δ K trim i =   s g n d 0 .   Besides, from Eq. (4), it is easy to infer ϵ c >   ϵ d   at a given momentum . This implies that E > E + . Regarding the topology of our band system, this in turn leads to the conclusion that ν0 = 1. Thus, indeed the system is a strong topological (non-trivial) insulator.

3. Floquet Theory

The polarized periodic optical field provides a potent modus operandi to carry out theoretical proposition and experimental realization, manipulation, and detection of diverse unconventional/novel optical and electronic properties of materials, such as the realization of novel quantum phases without static counterparts like light induced quantum anomalous phase (QAH) phase [19], the topological phase transitions in semi-metals [21,22,23,24], the Floquet engin-eering of magnetism in topological insulator thin films [41,42], and so on. The exotic Floquet topological phases with a high tunability could be realized using the polarized periodic optical field. In fact, there has been an upsurge on experimental front in the search for topological states, in solid state [42], cold-atom [43] and optical systems [44], which are driven periodically. The circularly polarized optical field (CPOF) is described by a time-periodic (time period = T = 2π/ ω where ω is the frequency of light) gauge field. Upon using the Peierls substitution, lattice electrons couple to the electromagnetic gauge field. In the presence of COPF, the thin film Hamiltonian H s u r f a c e , , apart from breaking the time reversal invariance (TRS), becomes periodic in time. One can now transform the time-dependent Hamiltonian problem to a time-independent one using the Floquet’s theorem [25,26,27,28]. Analogous to the Bloch theory, a solution for the time-dependent Schrodinger equation of the system is obtained here involving the Floquet quasi-energy and the time-periodic Floquet state with the periodicity T. The Floquet state could be expanded in a Fourier series which makes us arrive at an infinite dimensional eigenvalue equation in the Sambe space [27].The circularly polarized optical field incident upon the film may be described by a time-varying gauge field A t   = A t + T =   A 1 sin ω t , sin ω t   +   ψ , 0  through the relation:  E t = A t t = E cos ω t , cos ω t   +   ψ , 0 , E = A 1 ω . Here E t   i s   The optical field. In particular, when the phase ψ = 0 or π, the optical field is linearly polarized. When ψ = + π/2 ( ψ = −π/2), the optical field is left-handed ( right-handed) circularly polarized. Once we have included a gauge field, it is necessary that we make the Peierls substitution H s u r f a c e t = H s u r f a c e k e A t . The dimensionless quantity a A 0 = a e E ω corresponds to the frequency of the incident light.
We assume the normal incidence of CPOF on the surface SmB6 with the thickness d = 30 nm. Suppose the angular frequency of the optical field incident on the film is ω   10 15 r a d i a n s 1 and wavelength λ i n 1500   n m . Therefore the ratio d/ λ i n   0.02 1 Upon taking the field into consideration our Hamiltonian becomes time dependent. As stated above, the Floquet theory can be applied to our time-periodic Hamiltonian H s u r f a c e t =   H s u r f a c e t + T with the period T = 2π/ω. Analogous to the Bloch theory involving crystal quasi-momentum, a solution | η   t = exp ( i έ t ) | ξ   t involving the Floquet quasi-energy έ could be written down for the time-dependent Schrodinger equation of the system. The Floquet state satisfies ξ   t =   ξ   t + T and, therefore, could be expanded in a Fourier series | ξ   t =   r exp i r ω t   | ξ r   where r is an integer. Then the wave function, in terms of the quasi-energy έ has the form | η   t = r exp i έ + r ω t ξ r   . This makes us arrive at an infinite dimensional eigenvalue equation in the Sambe space (the extended Hilbert space)[27,28]:
s H s u r f a c e , r , s   | ξ n s   = ( s ω δ r , s + 1 T 0 T H surface t e i r s ω t   d t )   ξ n s   =   έ n   ξ n s   .
The matrix element of the Floquet state surface Hamiltonian H s u r f a c e ,   α , β   is given by     H s u r f a c e , α , β = α ω δ α , β + 1 T 0 T H surface t e i α β ω t d t ,   where ( α , β ) are integers. In view of the Floquet theory [29,30,31,32], in the high-frequency limit, a thin film system, irradiated by the circularly polarized radiation, can be described by an effective, static Hamiltonian. in the off-resonant regime using the Floquet-Magnus (high-frequency) expansion [30]:
H surface F l o q u e t ( k ) = H s u r f a c e , 0,0 +   H s u r f a c e , 0 , 1 , H s u r f a c e , 0,1 ω + O ω 2 ,
where H s u r f a c e , n , m = 1 T 0 T H surface t e i n m ω t d t with n m . For M << t d 1 , we can write
H s u r f a c e , 0,0 = h surface k , q , μ , b + [ t d 1 + b 2 t f 1 2 + 2 t d 2 + b 2 t f 2 + 2 t d 3 + b 2 t f 3 ]
  × ( a 2 A 0 2 ) σ 0 τ 0 +   [   t d 1 b 2 t f 1 2 + 2 ( t d 2 b 2 t f 2 ) + 2 ( t d 3 b 2 t f 3 ) ]   ( a 2 A 0 2 )   σ 0 τ z ,
H s u r f a c e , 0 , 1 = t d 1 + b 2 t f 1 2 + 2 t d 2 + b 2 t f 2 + 2 t d 3 + b 2 t f 3 i a 2 ( k x + e i   ψ k y ) A 0 σ 0 τ 0 [   t d 1 b 2 t f 1 2 + 2 ( t d 2 b 2 t f 2 ) + 2 ( t d 3 b 2 t f 3 ) ] i a 2   k x + e i   ψ k y A 0   σ 0 τ z ( i / 2 ) A 1 ( a A 0 ) σ z  ⊗ τ x ( i / 2 ) A 1 a A 0   e i   ψ σ 0   τ y   ,
H s u r f a c e , 0,1 = t d 1 + b 2 t f 1 2 + 2 t d 2 + b 2 t f 2 + 2 t d 3 + b 2 t f 3 i a 2   ( k x + e i   ψ k y ) A 0 σ 0 τ 0 + [   t d 1 b 2 t f 1 2 + 2 ( t d 2 b 2 t f 2 ) + 2 ( t d 3 b 2 t f 3 ) ] i a 2   k x + e i   ψ k y A 0   σ 0 τ z + ( i / 2 ) A 1 ( a A 0 ) σ z   τ x + ( i / 2 ) A 1 ( a A 0 )   e i   ψ σ 0   τ y .
From the action of the time reversal operator on the wave function we see, that it leads to a complex conjugation of the wave function. Thus, in the case of spin-less wave functions as Θ =K, where K is the operator for complex conjugation. More generally, we can write Θ = UK where U is a unitary operator. Furthermore, for a spin-1/2 particle, flipping the spin coincides with the time-reversal. This means Θ  S   ^ = S   ^ where   S   ^ = 1 2   σ ^   and σ ^ is the vector of Pauli matrices. In view of these, one may also choose Θ = i σ y τ 0 K . Upon making use of the results Θ   A   ^ Θ −1 = A   ^ , Θ   B   ^ Θ −1 = B   ^ , and so on, where A   ^ =   σ 0 τ 0 ,   σ z τ y   ,   and B   ^ =   σ 0 τ y   , ……., we find that
Θ   H surface F l o q u e t a k x , a k y   Θ 1 = H surface F l o q u e t a k x ,   a k y + ( 4 a 2 A 0 2 s i n ψ / ω ) { A 1 a k x σ 0   τ x + A 1 a k y σ z τ y }   × t d 1 + b 2 t f 1 2 + 2 t d 2 + b 2 t f 2 + 2 t d 3 + b 2 t f 3 + 4 a 2 A 0   2 A 1 2 s i n ψ ω σ z τ z ,
where Θ   H surface F l o q u e t a k x , a k y     Θ −1 = H surface F l o q u e t k a k x ,   a k y only when ψ = 0   or π ,   that is , when the optical field is linearly polarized. In this case, the time reversal symmetry (TRS) is not broken. However, when ψ 0   or π , TRS is broken. We now consider the particular cases where ψ = + π/2 and ψ = −π/2. For the former the optical field is left-handed circularly polarized, whereas for the latter it is right-handed. Thus, the (previously not known) consequence is that the incidence of the CPOF on the SmB6 surface will be able to create a novel state with the broken TRS.
The Hamiltonian to describe this broken TRS system, in the basis ( d k ,     b c k ,   d k , b c k , ) T, could be written as H surface F l o q u e t k =
E 1 A 1 O P + ( a k ) 0 i   A 1 O P + a q A 1 O P + a k + E   2 i   A 1 O P + a q 0 0 i   A 1 O P a q E 3 A 1 O P ( a k ) i   A 1 O P a q 0 A 1 O P a k + E 4
where k = k x , k y ,   a k = a k x i a k y , E 1 = ϵ O P + ϑ O P + , E 2 = ϵ O P ϑ O P + , E 3 = ϵ O P + ϑ O P , E 4 = ϵ O P ϑ O P , A 1 O P ± = A 1 1 ± 2 B 2 s i n ψ a 2 A 0 2 ω , and A 1 = 2 V b . The functions ϵ O P = ϵ O P k , q , μ , b and ϑ O P ± = ϑ O P ± k , q , b are defined below:
ϵ O P k ,   q , μ , b = ϵ 0 μ , b D 1 b a 2 q 2 +   D 2 b a 2 k 2 + a 2 A 0 2 D 2 b + O a 4 k 4 ,
ϑ O P ± k , q , b = ϑ 0 b B 1 b a   2 q 2 + B 2 b a 2 k 2 a 2 A 0 2 B 2 ± a 2 A 0 2 ω s i n ψ   A 1 2 .
The eigenvalues ( ε α ) of the matrix (15) is given by the quartic ε α 4 + γ 1 O P k , b ε α 3 + γ 2 O P k , b ε α 2   + γ 3 O P k , b   ε α + γ 4 O P k , b = 0 (α= 1,2,3,4 ) where the coefficients γ β O P ( k , b) (β = 1,2,3,4) are given in Appendix A (see Eqs.(A.21)—(A.23)). It may be noted that to denote these eigenvalues we have used the symbol var epsilon which is distinct from that in Eq. (A.1). Once again, in view of the Ferrari’s solution of a quartic equation, we find the roots as
  ε α s , σ , k , b = s η O P k 2 γ 1 O P k , b 4 + l b O P k η O P k 2 +   s   c O P k 2 η O P k 1 2 .
where α = 1,2,3,4, s = ± 1 is the spin index and l   = ± 1 is the band-index. The spin-down ( s = 1 ) conduction band ( l   = + 1 ) and the spin-up (down) ( s = ± 1 ) valence bands ( l   = 1 ) , denoted respectively by ε 2 l = + 1 , s = 1 , k , b , ε 3 l = 1 , s = + 1 , k , b ,   and   ε 4 l = 1 , s = 1 , k , b somewhat peculiar as will be shown below. The functions appearing in Eq. (18) are given by
η O P k = 2 b O P k 3 + O P k 0 O P k 1 3 O P k + 0 O P k 1 3 ,
  0 O P ( k ) = ( b O P 3 k 27 b O P k d O P k 3 c O P 2 ( k ) ) ,
O P k = ( 2 729 b O P 6 + 4 d O P 2 b O P 2 27 + c O P 4 d O P b O P 4 81 2 b O P 3 27 + 2 c O P 2 b O P d O P 3 + d O P 3 27 ) 1 / 2
b O P k = { 3 γ 1 O P k , b 2 8 γ 2 O P k , b 16   } ,   c O P k = {   γ 1 O P k , b 3 + 4 γ 1 O P k , b γ 2 O P k , b 8 γ 3 O P k , b   32   } ,
d O P k = 3 γ 1 O P k , b 4 + 256 γ 4 O P k , b 64 γ 1 O P k , b γ 3 O P k , b + 16 γ 1 O P k , b 2 γ 2 O P k , b 256 .
The eigenvectors corresponding to ε α s   could be calculated in a manner given in the Appendix A. The value of a 2 A 0 2   dimensionless intensity of the radiation is taken to be around 0.8 which is good enough for the radiation field of frequency ν ~ 3 × 10 14 H z   under consideration. Moreover, s i n ψ = +1 ( s i n ψ = − 1 sign) corresponds to the left-handed (right-handed) circularly polarized radiation above. We notice from above that CPOF not only renormalizes d and f electron hopping integrals but also does the renormalization of the hybridization parameter(HP) .   We have the renormalized hybridization parameters(HP) as A 1 O P ± = A 1 1 ± 2 B 2 s i n ψ a 2 A 0 2 ω , and A 1 =   2 V b . As B 2 b [ t d 1 b 2 t f 1 2 ] > 0 ,   w e find that the renormalized HP A 1 O P + > A 1   ( A 1 O P + < A 1 ) for the left-handed (right-handed) CPOF. However, the renormalized HP A 1 O P < A 1 ( A 1 O P > A 1 ) for the left-handed(right- handed ) CPOF .   We note that, in principle, when a renor-malized parameter is less than A 1 , it is possible that there is a critical intensity of the radiation a 2 A 0 2   ω     t d 1 b 2 t f 1 at a given frequency at which the RHP in question will be zero. This, however, may affect the topological nature of the material. Now the nearest neighbor hopping elements t d 1   and   t f 1 are related to the band masses by t d 1 , f 1 2 2 m d , f a 2   . If one takes for the band masses m d   m f = 1.5   m e   50   m e , where m e is the rest mass of the electron , then the corresponding values of the hopping matrix elements are t d 1 ≈ 150 meV and t f 1 ≈ 4.5 meV. This yields the critical intensity of the radiation a 2 A 0 2 2.25 which is roughly three times the intensity value assumed in the graphical representations in Figure 3.
Figure 3. The plots of energy eigenvalues   ε α s , σ , k , b in Eq.(31) as a function of ak for a given intensity of incident radiation a A 0 = (0.80, 0.50). The Figures 4(a) and 4(b), respectively, corresponds to the plots for the left handed and the right-handed CPOF. The same is true for the Figures 4(c) and 4(b). The numerical values of the parameters used in the plots in (a) and (b) are t d 1 = 1, t f 1 = 0.8, t d 2 = 0.01, t f   2 = 0.01, t d 3 = 0.001, t f 3 = 0.001, ϵ f =   0.02, V = 0.1, b =0.98,μ = 0 ,   a A 0 = 0.8 ,   M = 0.0001   0 ,   a n d   U f t d 1 . The numerical values corresponding to (c) and (d) are μ = 0.5 ,   and   a A 0 = 0.5 ; the rest of the values are the same as in (a) and (b).The horizontal lines represent the Fermi energy.
Figure 3. The plots of energy eigenvalues   ε α s , σ , k , b in Eq.(31) as a function of ak for a given intensity of incident radiation a A 0 = (0.80, 0.50). The Figures 4(a) and 4(b), respectively, corresponds to the plots for the left handed and the right-handed CPOF. The same is true for the Figures 4(c) and 4(b). The numerical values of the parameters used in the plots in (a) and (b) are t d 1 = 1, t f 1 = 0.8, t d 2 = 0.01, t f   2 = 0.01, t d 3 = 0.001, t f 3 = 0.001, ϵ f =   0.02, V = 0.1, b =0.98,μ = 0 ,   a A 0 = 0.8 ,   M = 0.0001   0 ,   a n d   U f t d 1 . The numerical values corresponding to (c) and (d) are μ = 0.5 ,   and   a A 0 = 0.5 ; the rest of the values are the same as in (a) and (b).The horizontal lines represent the Fermi energy.
Preprints 82469 g003
In Figure 3 (a),3(b), 3(c), and 3(d) we have plotted the energy eigenvalues   ε α s , σ , k , b as a function of ak for a A 0 = (0.80,0.50) for the circularly polarized light. Whereas Figures 3(a) and 3(c) correspond to left-handed CPOF, 3(b) and 3(d) to the right-handed CPOF. The numerical values of the parameters used in the plots are t d 1 = 1, t f 1 = 0.8, t d 2 = 0.01, t f   2 = 0.01, t d 3 = 0.001, t f 3 = 0.001, ϵ f =   0.02, V = 0.1, b =0.98,μ = ( 0 , 0.5 ) ,       M = 0.0001   0 ,   a n d   U f t d 1 . The horizontal lines represent the Fermi energy. The conduction and valence bands denoted by ε 2 l = + 1 , s = 1 , k , b , ε 3 l = 1 , s = + 1 , k , b ,   and   ε 4 l = 1 , s = 1 , k , b   represent-ed by differently colored curves, apart from the band-inversion, exhibit some peculiarities by way of the multiple avoided crossings and the near absence of a surface Dirac cone at k = 0 in 3(a) and 3(b) unlike that in Figure 1(a).This non-trivial feature could be ascribed to the interaction of the system with the incident radiation. The figures show that when TRS is broken despite M = 0, the fledgling novel phase of the system is very robust. The reason being in both the figures the Fermi energy intersects the band ε 3 l = 1 , μ = + 1 , k , b only in the same BZ an odd pair number of times. This pair of surface state crossings (SSC)corresponds to the momenta k = ( ± 2,0 ) or ( 0 , ± 2 )   in Figure 3(a) and Figure 3(b). However, in Figure 3(c) and Figure 3(d) the same happens at the momenta k = ( ± 1,0 ) or ( 0 , ± 1 ) . These momenta satisfy the condition a k + a G =   a k , where the reciprocal lattice vector G is ( 4,0 ) or ( 0 , 4   )   in Figure 3(a) and Figure 3(b) and ( 2,0 ) or ( 0 , 2   ) in Figure 3(c) and Figure 3(d). O ur graphical representation lead to the fact that, due to the light-matter interaction, the emergent unconventional phase possibly corresponds to QSH. However, as stated in section 1, the conclusive evidence of this TRS-broken phase being QSH/QAH phase will be obtained once we calculate the spin Chern number[34] and the Z2 invariant which are future tasks.

4. Discussion and concluding remarks

The strong correlation effects and diverse surface conditions make SmB6 extremely complicated and almost a Gordian knot. Despite this, as we have seen above, our low-energy model was able to show that the compound is a strong TI. Our low-energy model was also able to capture the fact that there should be a surface Dirac cone at k=0(as in ref.[16]) in Figure 1(a) and 1(b) for M = 0. For M 0   ( Figure 1(d)), since TRS is broken, there is no Kramers degeneracy. By calculating BC and the Chern number we have been able to show that the Figure 1(d) corresponds to QAH state. In the case of the light-matter interaction in section 3, however, we need to show that the novel TRS-broken phase (despite M = 0) corresponds to QSH state. The problem needs an extensive investigation introducing an additional term
h z = 1 2   α 0   s i n k y a   σ x ( τ z + τ 0 + 1 2   α 0   s i n k x a   σ y ( τ z + τ 0 ) ] ,
which is the Rashba spin-orbit coupling (RSOC) between the d-electrons, in Eq. (1). Here α 0 stands for the strength of RSOC . These are highlights of the present report.
The Rashba coupling can arise in the present system due to proximity of material lacking in the structural inversion symmetry. In view of the spin-polarized ARPES measurements which appear to confirm the surface helical spin texture [45,46], it would be interesting to see how does surface state react to Rashba splitting as there is evidence of for a massive surface state at the surface Brillouin zone center which can exhibit Rashba splitting[47]. The Rashba SOC is of particular importance as it is a crucial ingredient for several spintronics and topological phenomena [48].
There are many other complications to bring home the point that the system needs concerted investigations. In fact, we are presently working on three of the several issues to be discussed below in brief: (i) Since we have considered surface state quite extensively, the next step forward is an investigation on the Kondo break down[49,50]. (ii)In YbB12, a finite residual temperature-linear term in the thermal conductivity κ/T(T → 0) is observed demonstrating the presence of gapless and itinerant neutral fermions [51]. On the other hand, κ/T (T → 0) in SmB6 has been controversial [52,53]. While κ/T(T → 0) of SmB6 has been reported to be very small but finite [54,55,56], the absence of κ /T(T → 0) has been reported in references [55,57]. It is worth mentioning that in ref.[58] κ/T(T → 0) has been shown to be finite. In view of the fact that TKI are found to be exceptionally sensitive to impurities, the issue needs to be looked into. (iii) Finally, in the quantum oscillation (QO) experiments of Li et al. [59], the signatures of two-dimensional Fermi surfaces supporting the presence of topological surface states were obtained. Various theoretical models were put forward that invoke novel itinerant low-energy neutral excitations [60] within the charge gap. These excitations were proposed to produce magneto quantum oscillation (MQO) signals. The theoretical models which entered the fray are based on magnetoexcitons [61], scalar Majorana fermions [62], emergent fractionalized quasiparticles [63] and non-Hermitian states [64]. As has been reported earlier [61], in SmB6, the QOs are observed only in the magnetization (de Haas-van Alphen(dHvA)effect). The dHvA oscillations strongly deviate below 1 K from the Lifshitz-Kosevich theory (LKT) possibly due to the presence of magnetic impurities [65]. It must be mentioned that the QOs in YbB12 [66] are observed in both magnetization (the de Haas-van Alphen, dHvA, effect) and resistivity (the Shubnikov-de Haas, SdH, effect) at applied magnetic fields H where the hybridization gap is still finite. The temperature-dependence of the oscillation amplitude complies with the expectations of Fermi-liquid theory [66]. It is hoped that the details of the problematic issues given above, related/unrelated to the present communication, will motivate the condensed matter physics community to delve deeper into this problem.
In conclusion, looking at the controversies and the possibilities, it is anybody’s guess that there are many unsettled issues. Unless other TKI candidates are discovered and thoroughly studied, it is perhaps difficult to achieve enhancement in the current understanding of strongly correlated topological insulators. In this backdrop, it is pertinent to make an attempt to investigate thoroughly what exactly are the physical explanation of the issues involved. In a future communication, as we already stated, we undertake a part of this demanding task.

Appendix A

Eigenvalues and eigenvectors of the matrix in Eq. (3)
The eigenvalues j of the Hamiltonian matrix (12) are given by the quartic
j 4   + a   j 3 + b   j 2   + c   j + d = 0 .
In view of the Ferrari’s solution of a quartic equation, we find the roots as
  j s , σ , k , b ,   M = j , k s , σ b = σ η 0 k 2 a 4 + s b 0 k η 0 k 2 + σ   c 0 k 2 η 0 k   1 2   ,  
where j= 1,2,3,4, σ = ± 1 is the spin index and s   = ± 1 is the band-index. The coefficients (a, b, c, d) of the quartic are given by
a =   2   ( ϵ d + ϵ c ) ,   ϵ d = ( ϵ k ,   q , μ , b + ϑ k ,   q , b ) ,   ϵ c = ( ϵ k ,   q , μ , b ϑ k ,   q , b ) ,
b = A B + C D +   4   ϵ d ϵ c + ϵ c 2 + 2 ( d 1 2 + d 2 2 + d 3 2 )   ,  
c = 2 ( ϵ c   A B + ϵ d C D + ( ϵ d + ϵ c ) ( d 1 2 + d 2 2 + d 3 2 ) ) ,
d = [   A B C D + ( A D + B C ) ( d 1 2 + d 2 2 + A C + B D d 3 2 + ( d 1 2 + d 2 2 + d 3 2 ) 2 ,
A = ϵ d + M ,   B = ϵ d M ,   C = ϵ c + M ,   D = ϵ c M .
The following symbols are defined in the main text: d 0 = ϵ d ϵ c 2 , d 1 = i A 1 a k y ,   d 2 = i A 1 a k x ,   d 3 = A 1 a χ ,   and the hybridization parameter A 1 =   2 V b . The functions appearing in Eq. (A.2) are given by
  η 0 ( k ) = 2 b 0 ( k ) 3 + ( k ) 0 ( k ) 1 3 ( k ) + 0 ( k ) 1 3 ,   0 ( k ) = ( b   0 3 k 27 b 0 k d 0 k 3 c 0 2 ( k ) ) ,
k = ( 2 729 b 0 6 + 4 d 0 2 b 0 2 27 + c 0 4 d 0 b 0 4 81 2 b 0 3 27 + 2 c 0 2 b 0 d 0 3 + d 0 3 27 ) 1 / 2   ,   b 0 k = { 3 a 2 8 b 16   } ,
c 0 k = {   a 3 + 4 a b 8 c   32   } ,   d 0 k = 3 a 4 + 256 d 64 a c + 16 a 2 b 256 .
The eigenstates linked to the eigenvalues j s , σ , k   q , μ , b ,   M in (A.2) are given below. These are required for the calculation of the Chern number when M 0 .   These Bloch states are given by
  ψ j , k = ψ 1 j ( k ) ψ 2 j ( k ) ψ 3 j ( k ) ψ 4 j ( k ) , j = 1 , 2 , 3 , 4 ,
ψ   1 j k = g 1 j k , ψ 2 j k = ( i a k + )   f 2 j k g 1 j k ,   ψ 3 j k = i j ϵ d + M f 3 j k g 1 j k ,  
  ψ 4 j k = A 1 a k + f 4 j k g 1 j k ,
g 1 j k = ϒ j 1 / 2 k = (   1 + a 2 k 2 f 2 j k 2 + j ϵ d + M 2   f 3 j k 2 + A 1 2 a 2 k 2   f 4 j k 2 ) 1 2 ,
f 2 j k = ( a q ) 1 j ϵ d M j ϵ c + M + A 1 2 a 2 k 2 + A 1 2 a 2 q 2 × j ϵ d +   M j ϵ c + M + A 1 2 a 2 k 2 + A 1 2 a 2 q 2 1
f 3 j k = ( A 1 a q ) 1 j ϵ d M j ϵ c + M + A 1 2 a 2 k 2 + η A 1 2 a 2 q 2   × j ϵ d + M j ϵ c + M + A 1 2 a 2 k 2 + A 1 2 a 2 q 2 1 ,   η = j ϵ d M j ϵ d + M ,
f 4 j k = 2 M / j ϵ d + M j ϵ c + M + A 1 2 a 2 k 2 + A 1 2 a 2 q 2 .
In the special case M t d 1 ,   ψ 1 j k = ϒ j 1 / 2 k , ψ 2 j k = d 1 + i d 2 d 3 ϒ j 1 / 2 k , ψ 3 j k = ( j ϵ d ) d 3 ϒ j 1 / 2 k , and ψ 4 j k = 0 . We find that the function ϒ j 1 / 2 k is given by a q j ϵ d 2 A 1 2 + a k 2 + a q 2 1 / 2 . The eigenvectors in this special case are
ψ j , k a q j ϵ d 2 A 1 2 + a k 2 + a q 2 1 / 2 i a k + j ϵ d 2 A 1 2 + a k 2 + a q 2 1 / 2 i j ϵ d A 1 j ϵ d 2 A 1 2 + a k 2 + a q 2 1 / 2 0 .
Chern number calculation: We calculate here the intrinsic anomalous Hall conductivity (AHC)/Chern number to show that the system is in the quantum anomalous Hall (QAH)phase when M 0 . The expression of AHC is σ A H = ( e 2 h )   j B Z d 2 k 2 π 3   g E j k μ   j z ( k ) , where μ is the chemical potential of the fermion number, j is the occupied band index, g E j k μ   is the Fermi-Dirac distribution and j z k   is the z-component of the Berry curvature (BC) for the jth band. To obtain AHC, we calculate BC using the Kubo formula
j z k = 2   2 I m i j ( E j k E i k ) 2 j , k v x ^ i , k i , k v y ^ j , k .
Here k is the Bloch wave vector, E j k   is the band energy, | j , k are the Bloch functions of a single band. The operator v j ^ represents the velocity in the j direction. For a system in a periodic potential and its Bloch states as the eigenstates, in view of the Heisenberg equation of motion i   d x ^ d t = [ x ^ , H ^ ], the identity m , k v α ^ n , k     = 1   E j k E i k i , k k α j , k   is satisfied. Upon using this identity, we obtain AHC in the zero temperature limit as σ A H = C e 2   where C =   j C j , C j =   BZ Ω xy ( k ) d 2 k 2 π 2 . The z-component of the Berry-curvature (BC) is
Ω xy k   = j A j , y   k x A j , x   k y = 2   j I m   ψ j , k k x ψ j , k k y  
We use the results presented in (A.11) - (A.20) to calculate BC. Next, upon integrating BC on a k-mesh-grid of the Brillouin zone, we calculate the intrinsic AHC.
Eigenvalues of the matrix (15): The eigenvalues ( ε α ) of the matrix (15) is given by the quartic ε α 4 + γ 1 O P k , b ε α 3 + γ 2 O P k , b ε α 2   + γ 3 O P k , b   ε α + γ 4 O P k , b = 0 (α= 1,2,3,4 ) where the coefficients γ β O P ( k , b) (β = 1,2,3,4) are given by
γ 1 O P k , b = μ E μ , γ 2 O P k , b = ( 1 2 ) μ ν E μ E ν 2 ( A 1 O P + A 1 O P ) ( a q ) 2 ( A 1 O P + 2 + A 1 O P 2 ) ( a k ) 2 ,
γ 3 O P k , b = 1 6 μ ν σ E μ E ν E σ + ( A 1 O P + A 1 O P ) ( a q ) 2 μ E μ + A 1 O P 2 ( a k ) 2 E 1 + E 2 +   A 1 O P + 2 ( a k ) 2 E 3 + E 4 ,
γ 4 O P k = μ E μ A 1 O P 2 ( a k ) 2 E 1 E 2 A 1 O P + 2 ( a k ) 2 E 3 E 4 A 1 O P + A 1 O P ( a q ) 2 E 1 E 4 + E 2 E 3 + A 1 O P 2 A 1 O P + 2 ( a k ) 4 + A 1 O P 2 A 1 O P + 2 ( a q ) 4 2 A 1 O P 2 A 1 O P + 2 ( a q ) 2 ( a k x 2 a k y 2 ) .

Appendix B

Supporting result for the calculation of Z2 invariant of the system: We feel necessary to supplement Eq.(8) and note below this equation by presenting an outline of the reports in ref.[15] in an effort to make the present communication comprehensive. We consider the green and red bands in Figure 1(a) and 1(b), and denote their Bloch wave functions, respectively, as | ύ 1 k   and | ύ 2 k . Following Fu and Kane[15] we assume the system one dimensional, i.e. k = (k,0).We impart the time dependence assuming that the band parameters change with time and return to the original values at t = T. We also suppose that the Hamiltonian h surface B H Z   ( M = 0 ) satisfies the following condition h surface B H Z   ( M = 0 , t ) =   h surface B H Z   M = 0 , t 1   . It is well-known [15] that charge polarization P can be calculated by integrating the Berry connection of the occupied states over the BZ. In the present case of the two-band system, P may be written as P = P1 + P2 where the Berry connections i ύ α k k |   ύ α k   are given   by   c α α k   α = 1,2 and   P 1 = π π d k 2 π   c 11 k ,   P 2 = π π d k 2 π     c 22 k .   The Berry curvature is given by α(k) = ∇k ×   c α α k . The total polarization density C(k) = c 11 k + c 22 k . This yields the TR polarization which is defined by P t r = P1 – P3 = 2P1 − P. Here P t r gives the difference in charge polarization between spin-up and spin-down quasiparticle bands. We now go back the charge polarization P, calculated by integrating the Berry connection of the occupied states over the BZ. Furthermore, the time-reversed version of | ύ 2 ( k ) is equal to | ύ 1 k   except for a phase factor. Hence, at t = 0 and t = T/2 one may write Θ| ύ 2 k = e i ρ k | ύ 1 k and Θ| ύ 1 k =   e i ρ k | ύ 2 k .   It is not difficult to see that the matrix representation of the TR operator Θ in the Bloch wave function basis ξ α β k     ύ α k   |   ύ β k will now be given as   ξ k = 0 e i ρ k e i ρ k 0 . One can easily confirm that ξ α β k   is a unitary matrix. It has the property ξ α β   k =     ξ α β   k . This implies that ξ α β   k is anti-symmetric at a TRIM. Upon getting back to the connections, we note that the connections satisfy c 11 k = c 22 k   k ρ k . These lead us to the charge polarization between spin-up bands P1 as P 1 = 0 π d k 2 π   C k i 2 π ρ π ρ 0 . Since ρ k = i log ξ 12 ( k ) , and C(k) = t r ( c k ) , after a little algebra, we find   P t r   = 1 i π log ξ 12 0 2 ξ 12 0   ξ 12 π ξ 12 π 2 . Obviously enough, the argument of the logarithm is +1 or −1. Furthermore, since log(−1) = one can see that P t r is 0 or 1 (mod 2).Physically, of course, the two values of P t r corresponds to two different polarization states which the system can take at t = 0 and t = T/2.The Bloch functions | ύ α k , t introduced above could be visualized as maps from the 2D phase space (k, t) to the Hilbert space. As in refs.[15], the Hilbert space could be separated into two parts depending on the difference in P t r between t = 0 and t = T/2. This leads to introduction of a quantity ν0, specified only in mod 2, and defined as ( P t r T 2   P t r ( 0 ) ) : When P t r   does not changes between t = 0 to T/2, ν0 = 0, whereas if there is a change then ν0 = 1. The visualization mentioned above, thus, yields that the Hilbert space is trivial if ν0 = 0, while for ν0 = 1 it is nontrivial (twisted). Equivalently, the system band structures are characterized by Kane–Mele index[15] Z2= +1 (ν0 = 0) and Z2= −1 (ν0 = 1). We obtain
( 1 )   ν 0 = j ξ 12 k trim j ξ 12 k trim j 2  
We now consider the generalization of this result. For this purpose, we suppose that 2N bands are occupied and forming N Kramers pairs. For eachsuch pair n, at the TR symmetric times t = 0 and π the wave functions are related by write Θ| ύ n 2 k = e i ρ n k | ύ n 1 k and Θ| ύ n 1 k = e i ρ n k | ύ n 2 k and the matrix ξ k   =
0 e i ρ 1 ( k ) 0 0 e i ρ 1 ( k ) 0 0 0 0 0 0 e i ρ 2 ( k ) 0 0 e i ρ 2 ( k ) 0 .
This leads to ξ 12 k trim j ξ 34 k trim j .   ξ 2 N   2 N 1 k trim j = e i n = 1 N ρ n k trim j = Pf   ξ k trim j ,   where Pfaffian is defined for an antisymmetric matrix and is related to the determinant by Pf [ A ] 2 =   det A . The difference in charge polarization between spin-up and spin-down quasiparticle bands, in this general case, is now given by P t r   = 1 i π log det [ ξ 0 ] P f [ ξ 0 ]   P f { ξ π ] det [ ξ 0 ] . Thus, the Z2 topological invariant ν is given by ( 1 ) ν 0 = j δ k trim j where for each TRIM k trim j one defines
δ K trim j  ≡  Pf [ ξ K trim j ] det [ ξ K trim j ]
As we see in the section 2 of the main text, this leads to the classification of the Hilbert space into the twisted ( ν 0 =1) and the trivial one ( ν 0 = 0).

References

  1. Amorese, O. Stockert, K. Kummer, N. B. Brookes, Dae-Jeong Kim, Z. Fisk, M. W. Haverkort, P. Thalmeier, L. Hao Tjeng, and A. Severing, Phys. Rev. B 100, 241107(R) (2019). [CrossRef]
  2. S. R. Panday and M. Dzero, Interacting fermions in narrow-gap semiconductors with band inversion, J. Phys.: Condens. Matter 33 275601(2021). [CrossRef]
  3. A. Stern, M. Dzero, V.M. Galitski, Z. Fisk and J. Xia, Surface-dominated conduction up to 240 K in the Kondo insulator SmB6 under strain, Nat. Mater. 16 708 (2017). [CrossRef]
  4. P. P. Baruselli and M. Vojta, Kondo holes in topological Kondo insulators: Spectral properties and surface quasiparticle interference, Phys. Rev. B 89, 205105 (2014). [CrossRef]
  5. L. Miao, C.-H. Min, Y. Xu, Z. Huang, E. C. Kotta, R. Basak, M. S. Song, B. Y. Kang, B. K. Cho, K. Kißner, F. Reinert, T. Yilmaz, E. Vescovo, Yi-D. Chuang, W. Wu, J. D. Denlinger, and L. A. Wray, Phys. Rev. Lett. 126, 136401 (2021).
  6. M. Dzero, K Sun, V. Galitski, P. Coleman, Phys. Rev. Lett. 2010, 104, 106408(2010); ibid Rev. B 2012, 85, 045130 ( 2012).
  7. D. J. Kim, S.Thomas.; T. Grant, J. Botimer, Z. Fisk, J.Xia, Sci. Rep., 3, 3150(2013).
  8. S. Wolgast, C. Kurdak, K. Sun, J. W. Allen, D.-J. Kim, and Z. Fisk, Phys. Rev. B 88, 180405 (2013).
  9. X. Zhang, N. P. Butch, P. Syers, S. Ziemak, R. L. Greene, and J. Paglione, Phys. Rev. X 3, 011011 (2013).
  10. D. J. Kim, J. Xia, and Z. Fisk, Nat. Mater. 13, 466 (2014).
  11. P. K. Biswas, M. Legner, G. Balakrishnan, M. C. Hatnean, M. R. Lees, D. M. Paul, E. Pomjakushina, T. Prokscha, A. Suter, T. Neupert, et al., Phys. Rev. B 95, 020410 (2017).
  12. A. P. Sakhya and K.B. Maiti , Scientific Reports volume 10, Article number: 1262 (2020).
  13. N. Wakeham, P. F. S. Rosa, Y. Q. Wang, M. Kang, Z. Fisk, F. Ronning, and J. D. Thompson, Phys. Rev. B 94, 035127 (2016).
  14. Udai Prakash Tyagi, Kakoli Bera and Partha Goswami, On Strong f-Electron Localization Effect in a Topological Kondo Insulator, Symmetry,13(12), 2245 (2021); https://doi.org/10.3390/sym13122245. [CrossRef]
  15. M. Legner, Topological Kondo insulators: materials at the interface of topology and strong correlations(Doctoral Thesis),ETH Zurich Research Collection(2016).Link:https://www.research-collection.ethz.ch/bitstream/ handle/20.500. 11850/155932/eth-49918-02.pdf?isAllowed=y&sequence=2.
  16. F.Lou, J.-Z. Zhao, H.Weng, and X. Dai, Phys. Rev. Lett. 110, 096401(2013).
  17. B. A. Bernevig, T. L. Hughes, and S.-C. Zhang: Science 314, 1757 (2006).
  18. L. Fu, C.L. Kane.: Time reversal polarization and a Z2 adiabatic spin pump. Phys. Rev. B 74, 195312 (2006); L. Fu, C. L. Kane, and E.J. Mele: Topological insulators in three dimensions. Phys. Rev. Lett. 98, 106803 (2007); C. L. Kane, and E.J. Mele: Z2 Topological Order and the Quantum Spin Hall Effect. Phys. Rev. Lett. 95, 146802 (2005).
  19. S. S. Dabiri, H. Cheraghchi, and A. Sadeghi, Light-induced topological phases in thin films of magnetically doped topological insulators Physical Review B 103, 205130 (2021); H. Xu, J. Zhou, and J. Li, Light-induced quantum anomalous Hall effect on the 2D surfaces of 3D topological insulators, Advanced Science, 2101508 (2021).
  20. W. Zhu, M. Umer, and J. Gong, Floquet higher-order Weyl and nexus semimetals, Phys. Rev. Research 3, L032026 (2021). [CrossRef]
  21. L. Zhou, C. Chen, and J. Gong, Floquet semimetal with Floquet-band holonomy, Phys. Rev. B 94, 075443 (2016). [CrossRef]
  22. H. Hubener, M. A. Sentef, U. D. Giovannini, A. F. Kemper, and A. Rubio, Creating stable Floquet–Weyl semimetals by laser-driving of 3D Dirac materials, Nat. Commun. 8, 13940 (2017). [CrossRef]
  23. D. Zhang, H. Wang, J. Ruan, G. Yao, and H. Zhang, Engineering topological phases in the Luttinger semimetal α-Sn, Phys. Rev. B 97, 195139 (2018). [CrossRef]
  24. H. Liu, J.-T. Sun, and S. Meng, Engineering Dirac states in graphene: Coexisting type-I and type-II Floquet-Dirac fermions, Phys. Rev. B 99, 075121 (2019). [CrossRef]
  25. L. Li, C. H. Lee, and J. Gong, Realistic Floquet semimetal with exotic topological linkages between arbitrarily many nodal loops, Phys. Rev. Lett. 121, 036401 (2018). [CrossRef]
  26. X. Liu, P. Tang, H. H¨ubener, U. De Giovannini, W. Duan, and A. Rubio, arXiv preprint arXiv:2106.06977 (2021). 22. F. Qin, R. Chen, and H.-Z. Lu, J. Phys. Condens. Matter 34, 225001 (2022).
  27. H. Sambe, Steady states and quasi-energies of a quantum mechanical system in an oscillating field, Phys. Rev. A 7, 2203 (1973).
  28. A. A. Pervishko, D. Yudin, and I. A. Shelykh, Impact of high-frequency pumping on anomalous finite-size effects in three-dimensional topological insulators, Phys. Rev. B 97, 075420 (2018). [CrossRef]
  29. R. Chen, B. Zhou, and D.-H. Xu, Floquet Weyl semimetals in light-irradiated type-II and hybrid line node semimetals, Phys. Rev. B 97, 155152 (2018).; R. Chen, D.-H. Xu, and B. Zhou, Floquet topological insulator phase in a Weyl semimetal thin film with disorder, Phys. Rev. B 98, 235159 (2018). [CrossRef]
  30. N. Goldman and J. Dalibard, Periodically Driven Quantum Systems: Effective Hamiltonians and Engineered Gauge Fields, Phys. Rev. X 4, 031027 (2014); A. Eckardt and E. Anisimovas, High-frequency approximation for periodically driven quantum systems from a Floquet-space perspective, New J. Phys. 17, 093039 (2015). [CrossRef]
  31. N. Tsuji, T. Oka, and H. Aoki, Correlated electron systems periodically driven out of equilibrium: Floquet+DMFT formalism, Phys. Rev. B 78, 235124 (2008). [CrossRef]
  32. B. H. Wu and J. C. Cao, Noise of Kondo dot with ac gate: Floquet–Green’s function and noncrossing approximation approach,Phys. Rev. B 81, 085327 (2010). [CrossRef]
  33. S. Kohler, J. Lehmann, and P. H¨anggi, Driven quantum transport on the nanoscale ,Phys. Rep. 406, 379 (2005). [CrossRef]
  34. Y. Yang, Z. Xu , L. Sheng , B. Wang, D. Y. Xing , and D. N. Sheng, Time-reversal-symmetry-broken quantum spin Hall effect, Phys. Rev. Lett. 107, 066602(2011). [CrossRef]
  35. X.-L. Qi, Y.-S. Wu, and S.-C. Zhang, Topological quantization of the spin hall effect in two-dimensional paramagnetic semiconductors, Phys. Rev. B 74, 085308 (2006). [CrossRef]
  36. 36.W. Ruan, C. Ye, M. Guo, F. Chen, X.Chen, G.-Ming Zhang, and Y. Wang, Phys. Rev. Lett. 112, 136401( 2014).
  37. S. Rößler, T.-Hwan Jang, D.-Jeong Kim, , and S. Wirth , Proc. Natl. Acad. Sci. U.S.A., 111 (13) 4798 (2014).
  38. M. Neupane, N. Alidoust, S. Y. Xu, et al., Surface electronic structure of the topological Kondo-insulator candidate correlated electron system SmB6. Nat Commun 4, 2991 (2013). [CrossRef]
  39. J. Jiang, S. Li, T. Zhang, et al., Observation of possible topological in-gap surface states in the Kondo insulator SmB6 by photoemission. Nat Commun 4, 3010 (2013). [CrossRef]
  40. N. Xu, P. Biswas, J. Dil, et al., Direct observation of the spin texture in SmB6 as evidence of the topological Kondo insulator. Nat Commun 5, 4566 (2014). [CrossRef]
  41. F. Qin, R. Chen, and H.-Z. Lu, J. Phys. Condens. Matter 34, 225001 (2022).
  42. J. W. McIver, B. Schulte, F.-U. Stein, T. Matsuyama, G. Jotzu, G. Meier, and A. Cavalleri, Light-induced anomalous hall effect in graphene, Nature Physics 16, 38 (2019). [CrossRef]
  43. 43.B. K. Wintersperger, C. Braun, F. N. Unal, A. Eckardt, ¨ M. D. Liberto, N. Goldman, I. Bloch, and M. Aidelsburger, Realization of an anomalous floquet topological system with ultracold atoms, Nature Physics 16, 1058 (2020). [CrossRef]
  44. C. S. Afzal, T. J. Zimmerling, Y. Ren, D. Perron, and V. Van, Realization of anomalous Floquet insulators in strongly coupled nanophotonic lattices, Phys. Rev. Lett. 124, 253601 (2020). [CrossRef]
  45. N. Xu, P.K. Biswas, J.H. Dil, R.S. Dhaka, G. Landolt, S. Muff, C.E. Matt, X. Shi, N.C. Plumb, M. Radović, E. Pomjakushina, K. Conder, A. Amato, S.V. Borisenko, R. Yu, H.-M. Weng, Z. Fang, X. Dai, J. Mesot, H. Ding, M. Shi, Nat. Commun. 5, 4566 (2014).
  46. S. Suga et al., J. Phys. Soc. Jpn 83, 014705 (2014).
  47. Hlawenka, P., Siemensmeyer, K., Weschke, E. et al. Samarium hexaboride is a trivial surface conductor. Nat Commun 9, 517 (2018). https://doi.org/10.1038/s41467-018-02908-7. [CrossRef]
  48. A. Manchon, H. C. Koo, J. Nitta, S. M. Frolov, and R. A. Duine, New perspectives for rashba spin–orbit coupling, Nat. Mater. 871(2015). [CrossRef]
  49. Paul, C. P´epin, and M. Norman, Phys. Rev. Lett. 98, 026402 (2007).
  50. V. Alexandrov, P. Coleman, O. Erten, Phys. Rev. Lett. 114, 177202 (2015).
  51. Y. Sato et al., Nat. Phys. 15, 954 (2019).
  52. Fuhrman, W.T., Chamorro, J.R., Alekseev, P. et al. Screened moments and extrinsic in-gap states in samarium hexaboride. Nat Commun 9, 1539 (2018). https://doi.org/10.1038/s41467-018-04007-z. [CrossRef]
  53. In fact, we have the literature where it is reported that the material may show a linear T-term in the specific heat which has a coefficient that varies between 2 and 25 mJ/mole/K2 , depending upon the sample preparation [7, 16, 17]. Doping with 5% of magnetic impurities can lead to an order of magnitude increase in the heat capacity [18].
  54. M. E. Boulanger et al., Phys. Rev. B. 97, 245141 (2018).
  55. Y. Xu et al., Phys. Rev. Lett. 116, 246403 (2016).
  56. M. Orendáč, et al. Isosbestic points in doped SmB6 as features of universality and property tuning. Phys. Rev. B 96, 115101 (2017).
  57. S. Sen et al., Physical Review Research 2, 033370 (2020).
  58. Chowdhury, D., Sodemann, I. & Sentil, T., Nat. Commn. 9, 1766 (2018).
  59. G. Li et al., Science, 346, 1208–1212(2014). This study is the first report of the quantum oscillations in magnetization in Kondo insulators.
  60. G. Baskaran, arXiv: 1507.03477 v1 (2015).
  61. Knolle, J. & Cooper, N. R., Phys. Rev. Lett. 118, 096604 (2017).
  62. Erten, O., Chang, P.-Y., Coleman, P. & Tsvelik, A. M., Phys. Rev. Lett. 119, 057603 (2017).
  63. Chowdhury, D., Sodemann, I. & Sentil, T., Nat. Commn. 9, 1766 (2018).
  64. Shen, H. and Fu, L., Phys. Rev. Lett. 121, 026403 (2018).
  65. W. T. Fuhrman and P. Nikoli´c, Magnetic impurities in Kondo insulators: An application to samarium hexaboride, Phys. Rev. B 101, 245118 (2020).
  66. Z. Xiang, Y. Kasahara, T. Asaba, B. Lawson, C. Insman, Lu Chen, K. Sugimoto, S. Kawaguchi, Y. Sato, G. Li, S. Yao, Y.L. Chen, F. Iga, J. Singleton, Y. Matsuda, and L. Li, Quantum oscillations of electrical resistivity in an insulator, Science 69, 65 (2018). [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

Altmetrics

Downloads

153

Views

20

Comments

0

Subscription

Notify me about updates to this article or when a peer-reviewed version is published.

Email

Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2025 MDPI (Basel, Switzerland) unless otherwise stated