Next Article in Journal
Some Classical Models of Particles and Quantum Gauge Theories
Previous Article in Journal
Learning Mixed Strategies in Quantum Games with Imperfect Information
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Topological Properties of the 2D 2-Band System with Generalized W-Shaped Band Inversion

Department of Physics, Faculty of Science, University of Zagreb, Bijenička 32, 10000 Zagreb, Croatia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Quantum Rep. 2022, 4(4), 476-485; https://doi.org/10.3390/quantum4040034
Submission received: 26 September 2022 / Revised: 25 October 2022 / Accepted: 29 October 2022 / Published: 2 November 2022

Abstract

:
We report the topological properties, in terms of the Berry phase, of the 2D noninteracting system with electron–hole band inversion, described by the two-band generalized analogue of the low-energy Bernevig–Hughes–Zhang Hamiltonian, yielding the W-shaped energy bands in the form of two intersecting cones with the gap along the closed continuous loop. We identify the range of parameters where the Berry phase attains qualitatively different values: (a) the integer multiplier of 2 π , (b) the integer multiplier of π , and (c) the nontrivial value between the latter two, which depends on the system parameters. The system thus exhibits the anomalous quantum Hall effect associated with the nontrivial geometric phase, which is presumably tunable through the choice of parameters at hand.

1. Introduction

The Berry phase [1] was introduced at least conceptually for the first time most likely in the 1950s in D. Bohm’s Quantum Theory [2], Ch. 20, Sec. 1 in equation 8, as the geometric phase accumulated in the wave function during the cyclic adiabatic change of parameters in the Hamiltonian; today, it still grasps the focus of interest of the modern physics, particularly in the fields of condensed matter and optics [3]. The concept of nontrivial topological properties of the wave function is important not only from the fundamental point of view of the general quantum mechanics but also from the experimental and technological point of view. The topological insulators [4] and quantum Hall effect (QHE) [5,6] are just some of the most well-known examples of real systems with the above-mentioned effects experimentally observed. Furthermore, the anyons, firstly suggested by F. Wilcek as the composite particles consisting of the electric charge moving in the 2D plane around the infinite magnetic flux tube perpendicular to it [7], are objects satisfying neither fermionic nor bosonic statistics. In turn, the phase in the wave function, accumulated by braiding the two of them, is neither zero (bosons) nor π (fermions) but some real number between the two instead, reflecting the nontrivial topology of the wave function. The fundamental requirement to be able to achieve such a nontrivial braiding is the two-dimensionality (2D) of the system. So far, the only experimentally realized system exhibiting anyons in the field of solid-state physics is the fractional QHE [8,9], which is emerging due to electron–electron interactions. The extra efforts to realize anyons in different systems, especially the so-called non-Abelian ones, is also motivated by their potential application for quantum computation [10]. Some of these efforts are directed at finding the systems in which anyons may emerge in the absence of interactions (e.g., the electron–electron interaction) but solely due to the special topology of the wave function, resulting from diagonalization of the Hamiltonian containing parameters that can be controlled from the outside [11], the realization of which is still an open question.
We analyze the topological properties of the 2D system with inverted electron and hole bands. The band inversion and so-called “conical intersections” [12] are practically paradigmatic prerequisites for the nontrivial topological properties and onset of the topological insulators. Although the full scale of mechanisms leading to the band inversion is not fully known yet, usually, the spin–orbit interaction in the solid state is used as its main generator. The quest to generate the intrinsic QHE state, as an intrinsic property of the system without the external magnetic field (so-called anomalous QHE) [13], has been continued since Haldane’s proposal [14], its experimental realizations [15], as the important topic of modern solid-state physics and material science through the design of new topological materials. Here, we shall not discuss the mechanisms leading to the band inversion. We shall start with the Hamiltonian already possessing that feature and explore its properties. The system under consideration features the W-shaped bands with rotational symmetry (the so-called “Witch’s hat” form created by the two oppositely oriented intersecting cones with a gap created along the circular line of intersection). It is in essence the two-band derivative of the low-energy Bernevig–Hughes–Zhang Hamiltonian (BHZ) [16] for the 2D system, which is sometimes called the “half BHZ model”, leading to the modified Dirac equation [17], which we generalize in terms of the number of parameters presumably tunable by the external influence such as electric field, mechanical strain, etc. The results indicate regions of parameters for which the Berry phase can achieve different nontrivial values; i.e., we find regions in the parameter space with the nontrivial Berry phase characterized by the integer multiplier of π and of 2 π as well as the region where the Berry phase attains the nontrivial value between the latter two depending on the parameters in the Hamiltonian. Likewise, the anomalous QHE determined by it emerges, possessing the corresponding nonstandard properties.

2. Methods

2.1. The Model

The low-energy Hamiltonian of the 2D two-band system under consideration, in the basis of the plane waves characterizes by the 2D wave vector k = ( k x , k y ) , reads
H ^ = ( Δ α 1 | k | m ) α 2 ( k x i k y ) n α 2 ( k x + i k y ) n ( Δ α 1 | k | m ) ,
where Δ is one of the parameters determining the band gap, and parameters ( n , m ) { 1 , 2 , 3 , } are positive integers. In particular, m describes different possible initial electron dispersions such as linear, parabolic etc., while n describes the possibility of introducing more general topology beyond simple conical band intersections. The case n = 0 (and m = 0 ) is also covered by the result, but it leads to the trivial topology; thus, we shall comment on it at the end only for the sake of completeness. Here, α 1 ( α 2 ) stand for positive parameters giving, multiplied with k m ( k n ), dimension of energy, e.g., for m = n = 1 , α 1 , 2 are the Fermi velocity parameters α 1 , 2 = v 1 , 2 such as for graphene, for example; for m = n = 2 , α 1 , 2 are the effective mass ( m ) parameters α 1 , 2 = 2 / 2 m 1 , 2 such as for the free electron gas, etc.
Using the in-plane rotational symmetry of the problem, i.e., using the polar coordinates k = | k | > 0 and φ = arctan ( k y / k x ) , in which k x ± i k y = k exp ( ± i φ ) , and expressing energy E in units of α 1 , i.e., E E E / α 1 consequently scaling the Hamiltonian in the same way, i.e., H ^ H ^ H ^ / α 1 , we can write the Hamiltonian (1) in the convenient form for analytical analysis
H ^ = ( δ k m ) γ k n e i n φ γ k n e i n φ ( δ k m ) ,
where δ Δ / α 1 and γ α 2 / α 1 are the control parameters of the model, together with the corresponding exponents n and m. It is straightforward to diagonalize the Hamiltonian (2). The eigen-states, i.e., the electronic energy bands, are
E ± ( k ) = ± ( δ k m ) 2 + γ 2 k 2 n ,
where the sign “+” stands for the conduction (or electron) band, and sign “−” stands for the valence (or hole) band. The bands are shown schematically in Figure 1 for several choices of parameters.
The (normalized) eigen-vectors in the plane wave basis in polar coordinates, corresponding to the eigen-states (3), are
ψ ± ( k , φ ) = 1 2 ( δ k m ) 2 + γ 2 k 2 n 1 4 ± γ k n ( δ k m ) + ( δ k m ) 2 + γ 2 k 2 n 1 2 e i n φ ( δ k m ) + ( δ k m ) 2 + γ 2 k 2 n 1 2 .

2.2. The Berry Phase

Knowing the eigen-vectors, we can calculate the Berry connection A ( ± ) = ( A k ( ± ) , A φ ( ± ) ) in polar coordinates, where
A k ( ± ) ( k , φ ) = i ψ ± ( k , φ ) | k | ψ ± ( k , φ ) , A φ ( ± ) ( k , φ ) = i k ψ ± ( k , φ ) | φ | ψ ± ( k , φ ) ,
playing the role of the “Berry (vector) potential”, or synthetic gauge field generated by the “topological charge” located in the origin. Performing the partial differentiation k , φ and the straightforward algebra, we obtain
A k ( ± ) ( k , φ ) = 0 A φ ( ± ) ( k , φ ) = n 2 1 k ± δ k m k ( δ k m ) 2 + γ 2 k 2 n .
The Berry phase for each band is then calculated as the contour integral of A along the path encircling the “topological charge” in the origin. In the case of standard solid-state systems at the lattice, that path encircles the Brillouin zone, being the natural boundary along which the derivative of the wave function vanishes due to the periodicity condition. The topology of the 2D Brillouin zone is the one of the torus, granting the integer (or trivial, i.e., zero) TKNN invariant, also called the Chern number, in the standard noninteracting systems, thus reflecting its topological nature [18]. The presented model is approximate in the sense that it is continuous; we modeled just the vicinity of the conical intersection without full periodicity conditions of the lattice (so-called low-energy Hamiltonian). Thus, the path encircling the “topological charge” in the origin should be taken in the limit k . This reflects its geometric instead of pure topological nature (such as for the pure Dirac cone in graphene, for example). Therefore, the Berry phase B ( ± ) is
B ( ± ) = lim | k | A ( ± ) ( k ) d k = lim k 0 2 π A φ ( ± ) ( k ) k d φ
which transforms, after inserting the expression for the Berry connection (6) and taking the integral, into
B ( ± ) = n π 1 ± lim k δ k m ( δ k m ) 2 + γ 2 k 2 n ,
and, taking the limit, we finally obtain
B ( ± ) = π n 1 1 1 + γ 2 , n = m π n , n > m 0 or 2 n π , n < m .

2.3. The Anomalous QHE

Another way to address the intrinsic QHE originating solely from the band topology is the calculation of the electrical conductivity tensor within the adiabatic approximation, namely its off-diagonal element that accounts for the Hall conductivity. The contribution in the case of the full band accounts exactly for the QHE as illustrated in detail in early works of Streda, Kohomoto and others [19,20,21]. The simplest way is to calculate the average of the, say, x-component of the current operator as the response to the perpendicular electric field introduced through the time-dependent gauge, i.e., the vector potential of the form A = ( 0 , E y t , 0 ) , where E y is the perpendicular component of electric field, and t is time.
The current operator is defined by the expression
J ^ = e 1 ( 2 π ) 2 d k x d k y k H ^ ( k , t ) ,
where H ^ is the Hamiltonian of the system, and the Hall conductivity σ x y is defined as the response function, i.e., J ^ x = σ x y E y . The Hall conductivity in this case, per state | η with energy E η in the, by assumption, fully occupied band, is given by the well-known Kubo formula
σ x y η = i η η η | J x | η η | J y | η η | J y | η η | J x | η ( E η E η ) 2 .
In our case with two bands ν and ν , Equation (11), after summing up the contributions per states | k in the band, reduces to the expression for the Hall conductivity per band ν [22]
σ x y ( ν ) = k 2 Im J x ( ν ν ) ( k ) J y ( ν ν ) ( k ) ( E ν ( k ) E ν ( k ) ) 2 .
where we calculated it per unit area. To pursue this quest, one needs to calculate the current vertices
J κ ( ν ν ) ( k ) = e l l H l l k κ U l ν U l ν ,
where κ { x , y } accounts for the Cartesian component and l , l , ν , ν are the indices denoting (row, column) the corresponding matrix elements. U l ν are elements of the k -dependent unitary matrix
U ^ = e i ϕ k cos ( ϑ k / 2 ) e i ϕ k sin ( ϑ k / 2 ) sin ( ϑ k / 2 ) cos ( ϑ k / 2 ) ,
determined by the condition U ^ H ^ U ^ = E ^ , where E ^ is the diagonal eigenvalue matrix and
H 12 = | H 12 | e i ϕ k , tan ϕ k = Im [ H 12 ] Re [ H 12 ] , tan ϑ k = 2 H 12 H 11 H 22 .
Taking it into account, the current vertices (13) read
J κ ( ν ν ) ( k ) = e 1 2 cos ϑ k ( H 11 H 22 ) k κ + sin ϑ k | H 12 | k κ , J κ ( ν ν ) ( k ) = e 1 2 sin ϑ k ( H 11 H 22 ) k κ + i | H 12 | ϕ k k κ + cos ϑ k | H 12 | k κ ,
for the intraband ( ν ν ) and interband ( ν ν ) case, respectfully.
We illustrate the procedure and the result for the particular problem, using the Hamiltonian (1) for the case n = m = 1 . Having H 11 = H 22 = Δ α 1 | k | , H 12 = α 2 | k | exp ( i ϕ k ) , tan ϕ k = k y / k x and tan ϑ k = α 2 | k | / ( Δ α 1 | k | ) , we calculate the interband current vertices for our two-band system with the valence band ( ν = ) and conduction band ( ν = + )
J x ( + ) ( k ) = e α 2 | k | Δ E k k x i k y , J y ( + ) ( k ) = e α 2 | k | Δ E k k y + i k x ,
for which J κ ( + ) = [ J κ ( + ) ] , and where the spectrum, appearing in the expression, is E k = ( Δ α 1 | k | ) 2 + α 2 2 | k | 2 .
The (quantum) Hall conductivity of the system with a fully occupied valence band ( ν = ) and empty conduction band ( ν = + ), with E k E + ( k ) = E ( k ) , is given by the expression
σ x y = 2 d k x 2 π d k y 2 π 1 E k 2 Im J x + ( k ) J y + ( k ) = e 2 h α 2 2 Δ 2 0 k d k ( Δ α 1 k ) 2 + α 2 2 k 2 3 / 2 = e 2 h 1 2 1 + 1 1 + γ 2 .

3. Results

The most natural way to analyze the obtained results is to start with the properties of the electron energy bands (3) which also reflect the topological properties, depending on parameters. The band structure, depending on the characteristic choices of parameters, is shown in Figure 2. In principle, the energy bands have two intersecting cones with a gap opened between them along the circular intersecting line (we call it “the nodal ring”). The tips of the protruding cones are located symmetrically, beyond and beneath the origin, at the distance equal to 2 δ . The circular nodal ring, giving the W-shape (or “Mexican hat-shape”, depending on n, m), shrinks in radius as γ increases, but it is preserved up to the limit γ .
For the case n = m , the radius of the nodal ring is k 0 = ( δ / ( 1 + γ 2 ) ) 1 / n , and the real band gap, i.e., the vertical distance between the lower and upper band, along the nodal ring is equal to E + E = 2 δ γ / 1 + γ 2 (or E + E = 2 Δ γ / 1 + γ 2 in the original units of energy). In the general case n m , one needs to solve the equation n γ 2 k 0 2 n m + m k 0 m m δ = 0 for k 0 and insert it into expression (3) to find the value of the gap E + E .
The Berry phase of the system is a sum of the contributions of all occupied bands. We consider the case with only the lower band (“−”) fully occupied and upper band (“+”) empty. Then, the Berry phase of the system is equal to B = B ( ) , i.e.,
B = π n 1 + 1 1 + γ 2 , n = m π n , n > m 2 π n , n < m .
opening the possibility of different classes of nontrivial results dependent on the system parameters (but not on the scaling of energy). It is essential to stress that this result holds as long as we have two distinct bands, i.e., until the gap closes. Basically, the result exhibits three main classes depending on the relation of n and m: (a) if n < m , B attains the value 2 π n ; if n > m , B attains the value π n ; if n = m , B attains the value between 2 π n and π n depending on γ . We show the Berry phase for the case n = m in Figure 3. Apparently, as γ 0 , the band gap closes and B 2 n π . In the opposite case γ , the W-shape of the band is lost, yielding the value B n π .
The anomalous QHE is confirmed also by the “electrodynamical approach”, i.e., by the calculation of the Hall conductivity in the gapped state using the interband current vertices. The result (18), obtained for the case n = m = 1 , can be easily generalized using the analogous straightforward procedure for the case of general ( n , m ) { 1 , 2 , 3 , } . It yields the result σ x y = ( e 2 / h ) · C , where
C = 1 2 1 + 1 1 + γ 2 n , n = m 1 2 n , n > m n , n < m .
The conductivity appears in the units of conductivity quanta e 2 / h multiplied by the dimensionless quantity C which appears as the “winding number” equal to the Berry phase divided by 2 π , therefore being in accordance with Equation (19). In the next section, we shall discuss when it may have the meaning of the topological invariant known as the Chern number.

4. Discussion

Physical systems featuring different topological properties are nowadays well categorized in terms of classes of corresponding topological invariants depending on the symmetries that the system possesses [23]. The general topological properties of conical intersections, as generators of nontrivial synthetic gauge fields, and inverted bands are rather well-understood. The interesting upgrade to that established picture is to find a way how to tune or violate certain assumptions, upon which the symmetry classes are built, through the change of parameters, and then track the transitions between the topological phases. The model that we analyze is symmetric with respect to the time reversal; it is the low-energy approximation which reflects topological properties of the single (generalized) conical intersection, which is expressed in terms of the “winding number” C given by Equation (20). The character of the conical intersection, expressed by C , is determined by the system parameters m, n and γ = α 2 / α 1 , yielding the topological phase diagram with two phases and the unusual boundary between them—Figure 4.
The result reproduces, for example, the result for the modified Dirac equation, C = 1 [17], taking parameters m = 2 and n = 1 , or the K point of simple “massive graphene” [24] with C = 1 / 2 for m = 0 and n = 1 , or trivial insulator with C = 0 for n = 0 , m (e.g., n = 0 , m = 2 case is considered in Ref. [25]), etc. However, the analyzed system goes qualitatively beyond that. Apart from the trivial phase ( C = 0 ) for n = 0 , in the case of the nonzero n, the “winding number” C = 1 · n remains such for all m > n . It attains the value C = 1 / 2 · n for all m < n . Finally, for all m = n , C = n 1 + 1 / 1 + γ 2 / 2 . As a result of the lack of other effective means of visualization, we illustrate the effect of topology, “encoded” into the wave function (4), on the spin operator ( σ x , σ y , σ z ) , by calculating the vector field Σ ( k x , k y ) having the expected values of σ x and σ z operators as its components, i.e.,
Σ ( k x , k y ) = ( Ψ ( k x , k y ) | σ x | Ψ ( k x , k y ) , Ψ ( k x , k y ) | σ z | Ψ ( k x , k y ) ) ,
as depicted in Figure 5 (the topologically trivial case C = 0 would have all spins directed “upwards”).
To apply the presented results to the real systems requires taking into account all symmetries that it imposes. For example, the system on the lattice, such as graphene, is limited with symmetry yielding the necessity of an even number of Dirac points in the Brillouin zone. Those may cancel each other’s Berry phase if the Berry curvatures are of the opposite signs; thus, some time-reversal symmetry-lifting mechanism (for example, such as the Haldane model [14]) is necessary to observe the finite effect. In this quantum report, we do not go into such details but rather focus on reporting the unusual topological property of the quantum system (1) possessing the anomalous QHE and the non-integer “winding number” for the certain set of parameters.

Author Contributions

Both Z.R. and D.R. contributed to this paper in terms of conceptualization, methodology, formal analysis and draft preparation; D.R. also handled the project administration and funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the QuantiXLie Centre of Excellence, a project co-financed by the Croatian Government and European Union through the European Regional Development Fund—the Competitiveness and Cohesion Operational Programme (Grant KK.01.1.1.01.0004).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest whatsoever.

References

  1. Berry, M.V. Quantal phase factor accompanying adiabatic changes. Proc. R. Soc. Lond. Ser. A 1984, 392, 45–57. [Google Scholar]
  2. Bohm, D. Quantm Theory; Prentice-Hall, Inc.: New York, NY, USA, 1951; Chapter 20, Section 1. [Google Scholar]
  3. Cohen, E.; Larocque, H.; Bouchard, F.; Nejadsattari, F.; Gefen, Y.; Karimi, E. Geometric phase from Aharonov–Bohm to Pancharatnam–Berry and beyond. Nat. Rev. Phys. 2019, 1, 437–449. [Google Scholar] [CrossRef] [Green Version]
  4. Moore, J. The birth of topological insulators. Nature 2010, 464, 194–198. [Google Scholar] [CrossRef] [PubMed]
  5. Klitzing, K.V.; Dorda, G.; Pepper, M. New Method for High-Accuracy Determination of the Fine-Structure Constant Based on Quantized Hall Resistance. Phys. Rev. Lett. 1980, 45, 494–497. [Google Scholar] [CrossRef] [Green Version]
  6. Klitzing, K.V. The quantized Hall effect. Rev. Mod. Phys. 1986, 58, 519–531. [Google Scholar] [CrossRef]
  7. Wilczek, F. Magnetic Flux, Angular Momentum, and Statistics. Phys. Rev. Lett. 1982, 48, 1144–1146. [Google Scholar] [CrossRef]
  8. Tsui, D.C.; Stormer, H.L.; Gossard, A.C. Two-Dimensional Magnetotransport in the Extreme Quantum Limit. Phys. Rev. Lett. 1982, 48, 1559–1562. [Google Scholar] [CrossRef] [Green Version]
  9. Laughlin, R.B. Anomalous Quantum Hall Effect: An Incompressible Quantum Fluid with Fractionally Charged Excitations. Phys. Rev. Lett. 1983, 50, 1395–1398. [Google Scholar] [CrossRef] [Green Version]
  10. Kitaev, A. Fault-tolerant quantum computation by anyons. Ann. Phys. 1997, 303, 2–30. [Google Scholar] [CrossRef] [Green Version]
  11. Todorić, M.; Jukić, D.; Radić, D.; Soljačić, M.; Buljan, H. Quantum Hall Effect with Composites of Magnetic Flux Tubes and Charged Particles. Phys. Rev. Lett. 2018, 120, 267201. [Google Scholar] [CrossRef] [Green Version]
  12. Larson, J.; Sjöquist, E.; Öhberg, P. Conical Intersections in Physics, an Introduction to Synthetic Gauge Theories; Springer Nature: Cham, Switzerland, 2020; Chapter 4; pp. 55–91. [Google Scholar]
  13. Chang, C.-Z.; Zhang, J.; Feng, X.; Shen, J.; Zhang, Z.; Guo, M.; Li, K.; Ou, Y.; Wei, P.; Wang, L.-L.; et al. Experimental observation of the quantum anomalous Hall effect in a magnetic topological insulator. Science 2013, 340, 167–170. [Google Scholar] [CrossRef] [PubMed]
  14. Haldane, F.D.M. Model for a quantum Hall Effect without Landau levels: Condensed-matter realization of the ‘parity’ anomaly. Phys. Rev. Lett. 1988, 61, 2015–2018. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Jotzu, G.; Messer, M.; Desbuquois, R.; Lebrat, M.; Uehlinger, T.; Greif, D.; Esslinger, T. Experimental realization of the topological Haldane model with ultracold fermions. Nature 2014, 515, 237–240. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Bernevig, B.A.; Hughes, T.L.; Zhang, S.-C. Quantum Spin Hall Effect and Topological Phase Transition in HgTe Quantum Wells. Science 2006, 314, 1757–1761. [Google Scholar] [CrossRef] [Green Version]
  17. Peres, N.M.R.; Santos, J.E. Strong light—Matter interaction in systems described by a modified Dirac equation. J. Phys. Condens. Matter 2013, 25, 305801. [Google Scholar] [CrossRef] [Green Version]
  18. Thouless, D.J.; Kohmoto, M.; Nightingale, M.P.; den Nijs, M. Quantized Hall conductance in a two-dimensional periodic potential. Phys. Rev. Lett. 1982, 49, 405–408. [Google Scholar] [CrossRef] [Green Version]
  19. Streda, P. Theory of quantised Hall conductivity in two dimensions. J. Phys. C Solid State Phys. 1982, 15, L717–L721. [Google Scholar] [CrossRef]
  20. Kohmoto, M. Topological invariant and the quantization of the Hall conductance. Ann. Phys. 1985, 160, 343–354. [Google Scholar] [CrossRef]
  21. Kohomoto, M. Zero modes and the quantized Hall conductance of the two-dimensional lattice in a magnetic field. Phys. Rev. B 1989, 39, 11943–11949. [Google Scholar] [CrossRef]
  22. Sinitsyn, N.A.; Hill, J.E.; Min, H.; Sinova, J.; MacDonald, A.H. Charge and Spin Hall Conductivity in Metallic Graphene. Phys. Rev. Lett. 2006, 97, 106804. [Google Scholar] [CrossRef] [Green Version]
  23. Kitaev, A. Periodic table for topological insulators and superconductors. AIP Conf. Proc. 2009, 1134, 22–30. [Google Scholar]
  24. Xiao, D.; Yao, W.; Niu, Q. Valley-Contrasting Physics in Graphene: Magnetic Moment and Topological Transport. Phys. Rev. Lett. 2007, 99, 236809. [Google Scholar] [CrossRef] [PubMed]
  25. Rukelj, Z.; Akrap, A. Carrier concentrations and optical conductivity of a band-inverted semimetal in two and three dimensions. Phys. Rev. B 2021, 104, 075108. [Google Scholar] [CrossRef]
Figure 1. Schematics of the energy bands (3), conduction E + and valence E as a function of k assuming the rotational symmetry for several choices of system parameters. The vertical distance of the cone tips, along the energy axis, is 2 δ ( δ = 1 in all four panels). The other parameters are: (a) n = m = 1 , γ = 0 (the gap is closed); (b) n = m = 1 , γ = 0.5 ; (c) n = m = 2 , γ = 0.5 ; (d) n = 4 , m = 7 , γ = 0.5 .
Figure 1. Schematics of the energy bands (3), conduction E + and valence E as a function of k assuming the rotational symmetry for several choices of system parameters. The vertical distance of the cone tips, along the energy axis, is 2 δ ( δ = 1 in all four panels). The other parameters are: (a) n = m = 1 , γ = 0 (the gap is closed); (b) n = m = 1 , γ = 0.5 ; (c) n = m = 2 , γ = 0.5 ; (d) n = 4 , m = 7 , γ = 0.5 .
Quantumrep 04 00034 g001
Figure 2. The energy band E + ( k ) , assuming the angular rotational symmetry, for several choices of system parameters. The other energy band is simply E ( k ) = E + ( k ) . In both panels, δ = 0.5 , the energy gap between E + ( k ) and E ( k ) closes for γ = 0 and remains opened for finite γ varied between values zero and 20. The other parameters are: (a) n = m = 1 ; (b) n = m = 2 . The minimum of E + ( k ) determines the radius of the nodal ring (see Figure 1) k 0 which shrinks with increasing γ and becomes zero in the limit γ , turning the picture of two intersecting cones (with the gap opened along the intersection) into the picture of two distant cones, which are separated “vertically” by 2 δ from each other.
Figure 2. The energy band E + ( k ) , assuming the angular rotational symmetry, for several choices of system parameters. The other energy band is simply E ( k ) = E + ( k ) . In both panels, δ = 0.5 , the energy gap between E + ( k ) and E ( k ) closes for γ = 0 and remains opened for finite γ varied between values zero and 20. The other parameters are: (a) n = m = 1 ; (b) n = m = 2 . The minimum of E + ( k ) determines the radius of the nodal ring (see Figure 1) k 0 which shrinks with increasing γ and becomes zero in the limit γ , turning the picture of two intersecting cones (with the gap opened along the intersection) into the picture of two distant cones, which are separated “vertically” by 2 δ from each other.
Quantumrep 04 00034 g002
Figure 3. The Berry phase B (19) as a function of parameter γ for the case n = m . It attains a value of B = 2 π n for γ = 0 and tends to the value B π n as γ .
Figure 3. The Berry phase B (19) as a function of parameter γ for the case n = m . It attains a value of B = 2 π n for γ = 0 and tends to the value B π n as γ .
Quantumrep 04 00034 g003
Figure 4. Topological phase diagram of the system described by the Hamiltonian (1) characterized by the “winding number” C , given by Equation (20), depending on parameters ( n , m ) { 0 , 1 , 2 , 3 , } and γ R + . (0) n = 0 , m (orange): the trivial phase, C = 0 ; (I) m < n (green): C = n / 2 ; (II) m > n (blue): C = n ; (III) m = n 0 (red): C = n 1 + 1 / 1 + γ 2 / 2 . Along the m = n > 0 line, C is a function of γ which interpolates between values of C in the regions I and II, being equal to the value in I in the limit γ , and equal to the value in II in the limit γ = 0 .
Figure 4. Topological phase diagram of the system described by the Hamiltonian (1) characterized by the “winding number” C , given by Equation (20), depending on parameters ( n , m ) { 0 , 1 , 2 , 3 , } and γ R + . (0) n = 0 , m (orange): the trivial phase, C = 0 ; (I) m < n (green): C = n / 2 ; (II) m > n (blue): C = n ; (III) m = n 0 (red): C = n 1 + 1 / 1 + γ 2 / 2 . Along the m = n > 0 line, C is a function of γ which interpolates between values of C in the regions I and II, being equal to the value in I in the limit γ , and equal to the value in II in the limit γ = 0 .
Quantumrep 04 00034 g004
Figure 5. The vector field representing the expected value of ( σ x , σ z ) operator with respect to the state Ψ ( k x , k y ) (4) on the ( k x , k y ) domain. The horizontal labels (1), (2), (3) represent values of parameter γ = 0.1, 1, 10, respectively, while the vertical ones are: (a) ( m , n ) = ( 1 , 1 ) , C = −0.998, −0.854, and −0.651, respectively; (b) ( m , n ) = ( 1 , 2 ) , C = 1 / 2 , (c) ( m , n ) = ( 2 , 1 ) , C = 1 . The topologically trivial case for which the bands touch ( γ = 0 , C = 0 ) would have all “spins” directed upwards.
Figure 5. The vector field representing the expected value of ( σ x , σ z ) operator with respect to the state Ψ ( k x , k y ) (4) on the ( k x , k y ) domain. The horizontal labels (1), (2), (3) represent values of parameter γ = 0.1, 1, 10, respectively, while the vertical ones are: (a) ( m , n ) = ( 1 , 1 ) , C = −0.998, −0.854, and −0.651, respectively; (b) ( m , n ) = ( 1 , 2 ) , C = 1 / 2 , (c) ( m , n ) = ( 2 , 1 ) , C = 1 . The topologically trivial case for which the bands touch ( γ = 0 , C = 0 ) would have all “spins” directed upwards.
Quantumrep 04 00034 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rukelj, Z.; Radić, D. Topological Properties of the 2D 2-Band System with Generalized W-Shaped Band Inversion. Quantum Rep. 2022, 4, 476-485. https://doi.org/10.3390/quantum4040034

AMA Style

Rukelj Z, Radić D. Topological Properties of the 2D 2-Band System with Generalized W-Shaped Band Inversion. Quantum Reports. 2022; 4(4):476-485. https://doi.org/10.3390/quantum4040034

Chicago/Turabian Style

Rukelj, Zoran, and Danko Radić. 2022. "Topological Properties of the 2D 2-Band System with Generalized W-Shaped Band Inversion" Quantum Reports 4, no. 4: 476-485. https://doi.org/10.3390/quantum4040034

Article Metrics

Back to TopTop