Robust midgap states in band-inverted junctions under electric and magnetic fields

  1. Álvaro Díaz-FernándezORCID Logo,
  2. Natalia del ValleORCID Logo and
  3. Francisco Domínguez-AdameORCID Logo

GISC, Departamento de Física de Materiales, Universidad Complutense, E-28040 Madrid, Spain

  1. Corresponding author email

This article is part of the Thematic Series "Topological materials".

Guest Editor: J. J. Palacios
Beilstein J. Nanotechnol. 2018, 9, 1405–1413.
Received 22 Sep 2017, Accepted 16 Apr 2018, Published 14 May 2018


Several IV–VI semiconductor compounds made of heavy atoms, such as Pb1−xSnxTe, may undergo band-inversion at the L point of the Brillouin zone upon variation of their chemical composition. This inversion gives rise to topologically distinct phases, characterized by a change in a topological invariant. In the framework of the k·p theory, band-inversion can be viewed as a change of sign of the fundamental gap. A two-band model within the envelope-function approximation predicts the appearance of midgap interface states with Dirac cone dispersions in band-inverted junctions, namely, when the gap changes sign along the growth direction. We present a thorough study of these interface electron states in the presence of crossed electric and magnetic fields, the electric field being applied along the growth direction of a band-inverted junction. We show that the Dirac cone is robust and persists even if the fields are strong. In addition, we point out that Landau levels of electron states lying in the semiconductor bands can be tailored by the electric field. Tunable devices are thus likely to be realizable, exploiting the properties studied herein.

Keywords: crystalline topological insulators; electric and magnetic fields; Landau levels; midgap states


In 1982, Thouless et al. [1] made a connection between the quantum Hall effect and a topological invariant, the so-called first Chern number [2]. The fact that a quantum Hall system was insulating in the bulk but had a quantized conductivity on the edge could be related to the non-trivial topology of the band structure. In 2006, topology came up to stage once again with the theoretical prediction by Bernevig et al. [3] of a topological insulating behaviour in a HgTe/CdTe quantum well. The difference between the latter and the quantum Hall system lies in the fact that the quantum well required no magnetic field at all, but just relativistic corrections (Darwin and mass–velocity interactions) large enough so as to invert the Γ6 and Γ8 bands [4]. The HgTe/CdTe quantum well possesses non-trivial edge states when a certain width is exceeded. In 2007, experiments verified this remarkable result and established the existence of the quantum spin Hall effect [5]. However, no clear signatures of conductance quantization have been observed yet [6,7].

Besides II–VI compound semiconductors, such as HgTe, IV–VI semiconductors support non-trivial edges states as well [8]. In this regard, Dziawa et al. reported evidence of topological crystalline insulator states in Pb1−xSnxSe [9]. High-resolution scanning tunneling microscopy studies of these topological crystalline insulators provided strong evidence of the coexistence of massless Dirac fermions, protected by crystal symmetry, with massive Dirac fermions consistent with crystal-symmetry breaking [10]. Recently, these results have received further support with the aid of Dirac Landau level spectroscopy [11,12].

Band-inverted structures were already studied back in the 1980’s and 1990’s under the name of band-inverted junctions, in which the fundamental gap has opposite sign on each semiconductor. A salient feature is the existence of interface states lying within the gap, provided that the two gaps overlap [13-17]. These states are protected by symmetry, and are responsible for the conducting properties of the surface. In IV–VI heterojunctions, such as PbTe/SnTe, interface states are accurately described by means of a two-band model using the effective k·p approximation [18,19]. The equation governing the conduction- and valence-band envelope functions reduces to a Dirac-like equation after neglecting far-band corrections. In view of this analogy, exact solutions can be then straightforwardly found by means of supersymmetric [16] or Green’s function approaches [20]. In the context of symmetry-protected topological phases, our model can be applied not only to topological crystalline insulators, like the ones mentioned above [8], but also to more general three-dimensional topological insulators, such as Bi2Se3, in contact with a trivial insulator [21,22]. In the former case, mirror symmetry makes it possible to define mirror Chern numbers, which determine the topological crystalline phase [8]. In the latter, time-reversal symmetry, parity and particle–hole symmetry allow us to define a topological index given by the sign of the Dirac mass [21].

In 1994, Agassi studied the case of a band-inverted junction with a magnetic field applied parallel to the junction [23]. This author showed that the Dirac point remains robust upon the application of a magnetic field of arbitrary strengths and that the Landau levels in the continuum split for non-zero values of the in-plane momentum in the direction perpendicular to the magnetic field. By means of the modern theory of symmetry-protected topological phases, the protection of the Dirac point can be understood in the case of topological crystalline insulators from the fact that a magnetic field perpendicular to a mirror plane renders a system that is still symmetric about that plane [8]. This is not the case in a magnetic field parallel to the mirror plane, where the Dirac cone turns into the usual relativistic Landau levels [13,15,24]. Going back to the parallel magnetic field, Agassi demonstrated that for large values of this in-plane momentum, the states evolve to the bulk Landau states and the midgap state becomes the zero Landau level, usual of these Dirac systems. The reason is that the in-plane momentum perpendicular to the magnetic field is proportional to the position of the Landau orbits. If it becomes very large and the magnetic length is at the same time small, which happens for large magnetic fields, then the orbits do not intersect the junction and they might not notice that boundary. Hence, the case of most interest is in the vicinity of low in-plane momentum perpendicular to the field, where the states differ the most from the Landau levels of the bulk and we can see the effects of the interface.

In this same topic of external fields applied to band-inverted junctions, we have recently studied band-inverted junctions based on IV–VI compounds using a two-band model when an electric field is applied along the growth direction [25]. We have demonstrated that the Dirac cone of midgap states is robust against moderate values of the electric field but Fermi’s velocity decreases quadratically with the applied field. The aim of this paper is to characterize electron states in band-inverted junctions using a two-band model in the presence of crossed magnetic and electric fields, the former parallel to the junction, the latter perpendicular to it. We show that the Dirac cone of midgap states arising in the single-junction configuration is robust against crossed electric and magnetic fields. In addition, Landau levels of electron states lying in the semiconductor bands can be tailored by the electric field. Finally, the electronic structure of band-inverted junctions when the magnetic field is applied along the growth direction, parallel to the electric field, will also be briefly discussed for comparison.

Theoretical model

We consider heterojunctions of IV–VI compound semiconductors, such as Pb1−xSnxTe and Pb1−xSnxSe. The latter are known to shift from being semiconductors to topological crystalline insulators due to the band inversion at the L points of the Brillouin zone as the Sn fraction increases [8,26,27]. In order to keep the algebra as simple as possible, we restrict ourselves to the symmetric heterojunction with same-sized and aligned gaps, as depicted in Figure 1a. This assumption simplifies the calculations while keeping the underlying physics [28]. Thus, a single and abrupt interface presents the following profile for the magnitude of the gap


where sgn(z) = |z|/z is the sign function. Here the Z-axis is parallel to the growth direction [111].


Figure 1: (a) [Graphic 1] and [Graphic 2] band-edge profile of an abrupt band-inverted junction with aligned and same-sized gaps, located at the XY-plane, and b) schematic view of the applied electric and magnetic fields.

The envelope functions of the electron states near the band extrema [Graphic 3] and [Graphic 4] in IV–VI compounds are determined from the following Dirac-like Hamiltonian [15,16,19]


Here α = (αx, αy, αz) and β denote the usual 4 × 4 Dirac matrices, [Graphic 5] and [Graphic 6], σi and [Graphic 7] being the Pauli matrices and n × n identity matrix, respectively. Moreover, v is an interband matrix element having dimensions of velocity and it is assumed scalar, corresponding to isotropic bands around the L point. It is worth mentioning that the bands of IV–VI compounds around the L points are actually anisotropic. Nevertheless, this anisotropy can be absorbed in the definition of the dimensionless parameters defined below. That is, it is possible to consider a direction-dependent velocity, but it will not change the results shown below, except for a proportionality constant in the definition of the dimensionless in-plane momenta (see [19,28] for further details). In addition, we focus on states close to one of the L points of the Brillouin zone [8] and neglect other valleys in what follows since midgap states are stable against gap opening by valley mixing. Also notice that only linear momentum terms are taken into account in Equation 2 but quadratic momentum terms could have an impact of the electronic levels [29,30]. However, the two-band model Hamiltonian (Equation 2) successfully describes the hybridization of interface states in band-inverted quantum wells [31], in perfect agreement with more elaborated models including quadratic momentum terms [30]. The Hamiltonian (Equation 2) acts upon the envelope function χ(r), which is a four-component vector composed of the two-component spinors χ+(r) and χ(r) belonging to the [Graphic 8] and [Graphic 9] bands. The interface momentum is conserved and the envelope function can be expressed as [Graphic 10], where it is understood that the subscript “[Graphic 11]” in a vector indicates the nullification of its z-component. In the case of aligned and same-sized gaps, it is found that [Graphic 12], with [Graphic 13] and the interface dispersion relation is a single Dirac cone [Graphic 14], where the origin of energy is taken at the center of the gaps [20]. v is the group velocity at the Fermi level in undoped samples and it will be referred to as Fermi velocity hereafter.

Electron states under crossed electric and magnetic fields

We now turn to the electronic states of a single band-inverted junction exposed to a perpendicular electric field [Graphic 15] and a parallel magnetic field [Graphic 16], as shown schematically in Figure 1b. By choosing the Landau gauge, the vector potential is given as [Graphic 17].

The electrostatic potential eFz and the vector potential A(z) only depend on the z-coordinate. Therefore, [Graphic 18] is a constant of motion and the envelope function can still be factorized to the form [Graphic 19]. Now the longitudinal envelope function [Graphic 20] satisfies the following Dirac equation:


where [Graphic 21] is given by Equation 2. To address this problem we shall follow the Feynman–Gell-Mann ansatz [32] and define a new four-component vector ψ(z) as


It is convenient to introduce the following dimensionless quantities [Graphic 22], ξ = z/d, ε = E/Δ, f = eFd/Δ, and [Graphic 23]. Notice that f/2 is the ratio between the electric potential drop across the spatial extent of the midgap states [Graphic 24] in the absence of fields and the magnitude of the fundamental gap 2Δ. Similarly, b is the square of the ratio between d and the magnetic length [Graphic 25]. Hereafter we shall consider b > f ≥ 0 for reasons that will become clear shortly. Let us define


where μ = (b2f 2)1/4 is real. Then, inserting the ansatz (Equation 4) in Equation 3 and taking into account Equation 5, we get


where s0s(ξ = 0). Here [Graphic 26] and [Graphic 27] are 4 × 4 matrices given by


Let us diagonalize the left-hand side of the equation by introducing a unitary matrix U such that [Graphic 28]. Doing so and defining [Graphic 29] and [Graphic 30] we obtain


In order to solve Equation 8 we shall use the Green’s function method. The solution to Equation 8 will be given by


where the retarded Green’s function G(s,s′) satisfies


and G(s,s′)→0 as |s|,|s|→∞. Note that G(s,s′) is continuous on the line s = s′. Equation 9 can be particularized for s = s0, leading to a homogeneous system of equations with non-trivial solutions existing for energies satisfying the vanishing of the determinant


Since G(s,s′) can be considered as the inverse of the operator that acts upon it and the latter is diagonal, we may consider G(s,s′) to be block diagonal. Hence,


where [Graphic 31] is the 2 × 2 null matrix and the scalar functions g±(s,s′) satisfy


with [Graphic 32]. Since s is real because we have chosen μ to be so, then s2 > 0 and this equation corresponds to a harmonic oscillator. Notice that this would not be the case if μ were imaginary as in that case s2 < 0 and we would not have the positive parabola required for a harmonic oscillator. The solution to this problem is known to be [33,34]


where Γ(z) is the Gamma function, Dγ(z) is the parabolic-cylinder function, s> = max(s,s′) and s< = min(s,s′). Now that we have G(s,s′), it is straightforward to obtain from Equation 11 that g+(s0,s0)g(s0,s0) = μ2/2. Equivalently


Equation 15 determines the dispersion relation ε(κ) of midgap interface states as well as Landau levels lying in the semiconductor bands. It reduces to the result found by Agassi when the electric field vanishes [23].

Results and Discussion

Energy levels in the absence of electric field

This section is added for completeness and essentially reproduces the results found by Agassi [23] for small values of κx. However, we will be able to give approximate dispersion relations for the midgap state and the Landau levels that will provide us with a clearer view of the effect of the magnetic field in our case of interest. This section then corresponds to the f = 0 case, where approximate results can be obtained. In fact, these results are exact when κx = 0, where s0 = 0. Let us explore this last case. Using Γ(1 + z) = zΓ(z) and the Legendre duplication formula [Graphic 33], it is straightforward to obtain from Equation 15


There are now two possibilities, either the numerator goes to zero or the denominator goes to infinity. If p < 0, it is necessary to have a numerator equal to zero, which amounts to having,


This is nothing but a Dirac linear dispersion in the y-direction. It is remarkable that the Dirac point remains robust for any strength of b. Taking into account the definition of p, the case where p < 0 corresponds to |ε| < 1 at κx = κy = 0, meaning that these states lie within the gap.

Let us explore other possibilities. If p = 0, then both the numerator and the denominator are finite, which implies that p = 0 is not a solution. The other option where p > 0 is only achieved if the denominator goes to infinity because the numerator is always positive in this case. For this to happen, p must be a positive integer. The corresponding energies are the usual Landau levels of a relativistic particle


There is no zero Landau level because the requirement of p > 0 implies |ε| > 1 at κx = κy = 0, that is, Landau levels exist outside the gap. With this results in mind, we can now turn to the case where κx ≠ 0, but s0→0. After some tedious algebra we arrive at the following expression:




If s0 = 0 we obtain again Equation 16, corresponding to κx = 0. Now if κx ≠ 0, then either the term in curly brackets is zero or the prefactor multiplying this term is zero. As before, if the prefactor is zero then p is a positive integer. However, that would imply two possible energies for each integer, but numerically we will show briefly that this is not the case. Thus, we must consider the term in curly brackets to be equal to zero. If we consider b→0, but at the same time κx→0 sufficiently fast so that s0→0, then it is not difficult to obtain for the states inside the gap


whereas for the Landau levels we obtain to lowest order in κx


where c(n) results from the expansion around integer values of p of η(p) + η−1(p) + 2 ≈ c(n)(pn)−2. For instance, c(1) = 2/π, c(2) = 1/π, c(3) = 3/2π,… Before we consider each case, it is important to mention that the approximation of low b corresponds to the range of interest in experiments since typically d ≈ 4.5 nm and as a result b = 0.5 corresponds to a very large magnetic field of about 16 T.

Let us now consider each case separately. On the one hand, Equation 21 corresponds to an elliptic cone and for b = 0 we recover the original Dirac cone. It is not only remarkable, as we mentioned above, that the Dirac point is robust, but also that the shape of the dispersion relation remains a cone but slightly widened in the x-direction, as shown in Figure 2. Hence, the Fermi velocity becomes anisotropic and can actually be modulated with the magnetic field. It is expected that the application of an electric field will lead to further reduction of the Fermi velocity [25]. We will prove later that this is actually the case.


Figure 2: Dirac cones with, b≠ 0, and without, b = 0, a magnetic field applied. The original cone is distorted along the x-direction and the Fermi velocity, i.e., the slope, becomes anisotropic.

In Figure 3a we show a comparison between the Fermi velocity in the x-direction (recall that it does not change in the y-direction) given by the numerical evaluation of Equation 15 and the approximation in Equation 21. The agreement is noteworthy for low values of b.


Figure 3: Comparison between exact and approximate results given by (a) Equation 21 and (b) Equation 22. In panel (a) the Fermi velocity along the x-direction, calculated as the slope of the dispersion relation, is substantially reduced and the agreement between the exact and approximate results is noteworthy up to b ≤ 0.2. In panel (b), the Landau level splitting in the x-direction is very well predicted even for b = 0.5, as shown for the first level.

We can now focus on the Landau levels given by Equation 22. As it can be seen, for non-zero values of κx, each Landau level at κx = 0 splits into two Landau levels at κx ≠ 0 due to the occurrence of a “±”-sign inside the square root. The comparison for the first Landau level, n = 1, between the approximate result and the numerical calculations from Equation 15 are shown in Figure 3b. In contrast to Figure 3a, there is still agreement between both approaches for a large field of b = 0.5.

Energy levels at finite electric field

Let us now bring our attention to the case of f ≠ 0. In contrast to the case of f = 0, we have been unable to obtain explicit expressions of the dispersion relation, but the numerics shows remarkable results. Let us focus first on the midgap states. Since the magnetic field did not erase the Dirac point and based on known results of a band-inverted junction under an electric field [25,28], it seems plausible to argue that the effect of the electric field will be to enhance the reduction of the Fermi velocity in the x-direction and to introduce a reduction in the y-direction as well, leaving however the Dirac point untouched. This is indeed what we observe and we show our results in Figure 4. The insets show the Fermi velocity reduction as a function of the electric field for a fixed value of b = 0.5. It is remarkable how the Fermi velocity along the x-direction is substantially decreased in band-inverted junctions subject to crossed magnetic and electric fields.


Figure 4: The additional effect of the electric field leads to a further reduction of the Fermi velocity in the x-direction and to a reduction along the y-direction as well. The Dirac point, however, remains robust. The insets show the Fermi velocity reduction as a function the electric field for a fixed magnetic field of b = 0.5.

We may now turn to the evolution of the Landau levels as a function of the electric field. For simplicity, we shall consider only the first Landau level. It is illustrative to consider first the evolution of the lowest point of the Landau bands, that is, [Graphic 34]. If the electric field is zero, we already know what the energy will be from the discussion above. However, as we turn on the electric field, a splitting similar to the one we had with κx begins to develop. This splitting increases with electric field, up to a point where it starts decreasing again as f approaches b. In the limiting case where fb, the splitting goes to zero, as we show in Figure 5 for b = 0.5.


Figure 5: Splitting of the Landau levels at [Graphic 35] and b = 0.5 as a function of the electric field. It is important to note that the Landau levels move below the band edge due to the bending by the electric field (see main text for details).

In Figure 5 it may be surprising to see that the Landau bands shift below the band edge, leading to the apparent and erroneous belief that the latter enter the band gap. The effect of the electric field is to bend the constant band edges shown in Figure 1a upwards due to the presence of the electrostatic potential eFz, and so the Landau levels of the conduction band can move towards lower energies as long as the corresponding wave functions are not inside the band gap in position space.

Finally, it deserves consideration the previous discussion for low values of [Graphic 36]. As we can see in Figure 6, the parabolic dispersion that we obtained in the y-direction in the absence of an electric field splits into two parabolic bands. However, it is more remarkable to see that, instead of obtaining a splitting similar to that in Figure 3b, the dispersion goes downwards.


Figure 6: First Landau level dispersions for b = 0.5 and f =0.499. In panel (a), the original parabolic dispersion along the y-direction splits into two parabolic dispersions with energies below the band edge for the chosen fields, whereas in panel (b), the previously obtained splitting in the x-direction is now exclusively downwards.

Electron states under perpendicular electric and magnetic fields

In the previous sections we considered electron states when the magnetic field is parallel to the band-inverted junction, as depicted in Figure 1. For completeness, we now briefly discuss the salient features of the energy spectrum when the electric and magnetic fields are both perpendicular to the junction. The vector potential is then given as [Graphic 37] in the Landau gauge and thus [Graphic 38]. Starting from the Dirac Equation 3 with this vector potential and using the Feynman–Gell-Mann ansatz (Equation 4), one is led to a two-dimensional Schrödinger equation in the XZ plane. The resulting equation turns out to be separable in the x and z coordinates and can be straightforwardly solved by Green’s function techniques. At low or moderate electric and magnetic fields (f < b < 1), the energy levels within the gap are found to be


where n = 0, 1,… The above expression resembles the Landau levels of relativistic particles for an effective dimensionless magnetic field beffb(1 − 5f/8)2b(1 − 5f/4). Therefore, the electric field decreases the Landau level spacing as in the previous field configuration. There is yet another way of interpreting this result. If we undo the change of variables, we obtain for the energy the usual expression for the Landau levels that develop from a Dirac cone, the same as in graphene,


but with a renormalized Fermi velocity,


where FC = Δ/ed. In [25], it was anticipated that this renormalization of the Fermi velocity in a band-inverted junction with a perpendicular electric field could be measured by means of magnetotransport experiments, a prediction that is confirmed here.


We have studied band-inverted junctions under crossed electric and magnetic fields, the electric field being applied along the growth direction. Electron states are described by a spinful two-band model that is equivalent to the Dirac model for relativistic electrons. The mass term is half the bandgap and changes its sign across the junction. For the sake of algebraic simplicity, we assumed same-sized and aligned gaps, although this is not a serious limitation to the validity of the results [28].

In the absence of external fields, it is well known that band-inverted junctions support topologically protected states located at the interface. Their energy lies within the common gap of the two semiconductors and the dispersion relation is a Dirac cone [13,15,16,20]. The Dirac cone remains even if an electric field perpendicular to the junction is applied, but it widens and the Fermi velocity is quadratically reduced with the electric field [25,28]. In this paper we have proved that electrons with energy within the gap still behave as massless fermions when an additional magnetic field parallel to the band-inverted junction is applied. The original Dirac cone widens only in the direction perpendicular to the magnetic field but remarkably the dispersion relation remains gapless. Hence, the Fermi velocity becomes anisotropic and the combination of both electric and magnetic fields allows the Fermi velocity to be finely tuned. In addition, states lying within the semiconductor bands display relativistic-like Landau levels that split upon the application of the magnetic and electric fields. Interestingly, if both fields are parallel to the growth direction, the Landau level spacing can be further reduced by the electric field. We expect that the control of the Fermi velocity of topologically protected states will have applications for the design of novel electronic devices based on topological materials.


The authors thank L. Chico and J. W. González for helpful discussions. This work was supported by the Spanish MINECO under grant MAT2016-75955.


  1. Thouless, D. J.; Kohmoto, M.; Nightingale, M. P.; den Nijs, M. Phys. Rev. Lett. 1982, 49, 405. doi:10.1103/PhysRevLett.49.405
    Return to citation in text: [1]
  2. Kane, C. L.; Mele, E. J. Phys. Rev. Lett. 2005, 95, 146802. doi:10.1103/PhysRevLett.95.146802
    Return to citation in text: [1]
  3. Bernevig, B. A.; Hughes, T. L.; Zhang, S.-C. Science 2006, 314, 1757. doi:10.1126/science.1133734
    Return to citation in text: [1]
  4. Chu, J.; Sher, A. Physics and Properties of Narrow Gap Semiconductors; Springer: Berlin, Germany, 2008; p 157.
    Return to citation in text: [1]
  5. König, M.; Wiedmann, S.; Brüne, C.; Roth, A.; Buhmann, H.; Molenkamp, L. W.; Qi, X.-L.; Zhang, S.-C. Science 2007, 318, 766. doi:10.1126/science.1148047
    Return to citation in text: [1]
  6. Gusev, G. M.; Kvon, Z. D.; Shegai, O. A.; Mikhailov, N. N.; Dvoretsky, S. A.; Portal, J. C. Phys. Rev. B 2011, 84, 121302. doi:10.1103/PhysRevB.84.121302
    Return to citation in text: [1]
  7. Grabecki, G.; Wróbel, J.; Czapkiewicz, M.; Cywiński, Ł.; Gierałtowska, S.; Guziewicz, E.; Zholudev, M.; Gavrilenko, V.; Mikhailov, N. N.; Dvoretski, S. A.; Teppe, F.; Knap, W.; Dietl, T. Phys. Rev. B 2013, 88, 165309. doi:10.1103/PhysRevB.88.165309
    Return to citation in text: [1]
  8. Hsieh, T. H.; Lin, H.; Li, J.; Duan, W.; Bansil, A.; Fu, L. Nat. Commun. 2012, 3, 982. doi:10.1038/ncomms1969
    Return to citation in text: [1] [2] [3] [4] [5] [6]
  9. Dziawa, P.; Kowalski, B. J.; Dybko, K.; Buczko, R.; Szczerbakow, A.; Szot, M.; Łusakowska, E.; Balasubramanian, T.; Wojek, B. M.; Berntsen, M. H.; Tjernberg, O.; Story, T. Nat. Mater. 2012, 11, 1023. doi:10.1038/nmat3449
    Return to citation in text: [1]
  10. Okada, Y.; Serbyn, M.; Lin, H.; Walkup, D.; Zhou, W.; Dhital, C.; Neupane, M.; Xu, S.; Wang, Y. J.; Sankar, R.; Chou, F.; Bansil, A.; Hasan, M. Z.; Wilson, S. D.; Fu, L.; Madhavan, V. Science 2013, 341, 1496. doi:10.1126/science.1239451
    Return to citation in text: [1]
  11. Serbyn, M.; Fu, L. Phys. Rev. B 2014, 90, 035402. doi:10.1103/PhysRevB.90.035402
    Return to citation in text: [1]
  12. Phuphachong, T.; Assaf, B. A.; Volobuev, V. V.; Bauer, G.; Springholz, G.; de Vaulchier, L.-A.; Guldner, Y. Crystals 2017, 7, 29. doi:10.3390/cryst7010029
    Return to citation in text: [1]
  13. Volkov, B. A.; Pankratov, O. A. Sov. Phys. - JETP 1985, 42, 178.
    Return to citation in text: [1] [2] [3]
  14. Korenman, V.; Drew, H. D. Phys. Rev. B 1987, 35, 6446. doi:10.1103/PhysRevB.35.6446
    Return to citation in text: [1]
  15. Agassi, D.; Korenman, V. Phys. Rev. B 1988, 37, 10095. doi:10.1103/PhysRevB.37.10095
    Return to citation in text: [1] [2] [3] [4]
  16. Pankratov, O. A. Semicond. Sci. Technol. 1990, 5, S204. doi:10.1088/0268-1242/5/3S/045
    Return to citation in text: [1] [2] [3] [4]
  17. Kolesnikov, A. V.; Silin, A. P. J. Phys.: Condens. Matter 1997, 9, 10929. doi:10.1088/0953-8984/9/49/012
    Return to citation in text: [1]
  18. Kriechbaum, M. Envelope Function Calculations for Superlattices. In Two-Dimensional Systems: Physics and New Devices; Bauer, G.; Kuchar, F.; Heinrich, H., Eds.; Springer: Berlin, Germany, 1986; pp 120–129.
    Return to citation in text: [1]
  19. Ando, Y.; Fu, L. Annu. Rev. Condens. Matter Phys. 2015, 6, 361. doi:10.1146/annurev-conmatphys-031214-014501
    Return to citation in text: [1] [2] [3]
  20. Domínguez-Adame, F. Phys. Status Solidi B 1994, 186, K49. doi:10.1002/pssb.2221860231
    Return to citation in text: [1] [2] [3]
  21. Zhang, F.; Kane, C. L.; Mele, E. J. Phys. Rev. B 2012, 86, 081303. doi:10.1103/PhysRevB.86.081303
    Return to citation in text: [1] [2]
  22. Tchoumakov, S.; Jouffrey, V.; Inhofer, A.; Bocquillon, E.; Plaçais, B.; Carpentier, D.; Goerbig, M. O. Phys. Rev. B 2017, 96, 201302. doi:10.1103/PhysRevB.96.201302
    Return to citation in text: [1]
  23. Agassi, D. Phys. Rev. B 1994, 49, 10393. doi:10.1103/PhysRevB.49.10393
    Return to citation in text: [1] [2] [3]
  24. Pankratov, O. A.; Pakhomov, S. V.; Volkov, B. A. Solid State Commun. 1987, 61, 93. doi:10.1016/0038-1098(87)90934-3
    Return to citation in text: [1]
  25. Díaz-Fernández, A.; Chico, L.; González, J. W.; Domínguez-Adame, F. Sci. Rep. 2017, 8, 8058. doi:10.1038/s41598-017-08188-3
    Return to citation in text: [1] [2] [3] [4] [5]
  26. Assaf, B. A.; Phuphachong, T.; Volobuev, V. V.; Inhofer, A.; Bauer, G.; Springholz, G.; de Vaulchier, L. A.; Guldner, Y. Sci. Rep. 2016, 6, 20323. doi:10.1038/srep20323
    Return to citation in text: [1]
  27. Xu, S.-Y.; Liu, C.; Alidoust, N.; Neupane, M.; Qian, D.; Belopolski, I.; Denlinger, J. D.; Wang, Y. J.; Lin, H.; Wray, L. A.; Landolt, G.; Slomski, B.; Dil, J. H.; Marcinkova, A.; Morosan, E.; Gibson, Q.; Sankar, R.; Chou, F. C.; Cava, R. J.; Bansil, A.; Hasan, M. Z. Nat. Commun. 2016, 7, 12505. doi:10.1038/ncomms12505
    Return to citation in text: [1]
  28. Díaz-Fernández, A.; Domínguez-Adame, F. Phys. E (Amsterdam, Neth.) 2017, 93, 230. doi:10.1016/j.physe.2017.06.026
    Return to citation in text: [1] [2] [3] [4] [5]
  29. Kriechbaum, M.; Ambrosch, K. E.; Fantner, E. J.; Clemens, H.; Bauer, G. Phys. Rev. B 1984, 30, 3394. doi:10.1103/PhysRevB.30.3394
    Return to citation in text: [1]
  30. Buczko, R.; Cywiński, Ł. Phys. Rev. B 2012, 85, 205319. doi:10.1103/PhysRevB.85.205319
    Return to citation in text: [1] [2]
  31. Díaz-Fernández, A.; Chico, L.; Domínguez-Adame, F. J. Phys.: Condens. Matter 2017, 29, 475301. doi:10.1088/1361-648X/aa91a6
    Return to citation in text: [1]
  32. Feynman, R. P.; Gell-Mann, M. Phys. Rev. 1958, 109, 193. doi:10.1103/PhysRev.109.193
    Return to citation in text: [1]
  33. Domínguez-Adame, F. Europhys. Lett. 1991, 15, 569. doi:10.1209/0295-5075/15/6/001
    Return to citation in text: [1]
  34. Glasser, M. L.; Nieto, L. M. Can. J. Phys. 2015, 93, 1588. doi:10.1139/cjp-2015-0356
    Return to citation in text: [1]

© 2018 Díaz-Fernández et al.; licensee Beilstein-Institut.
This is an Open Access article under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
The license is subject to the Beilstein Journal of Nanotechnology terms and conditions: (

Back to Article List

Other Beilstein-Institut Open Science Activities

Keep Informed

RSS Feed

Subscribe to our Latest Articles RSS Feed.


Follow the Beilstein-Institut


Twitter: @BeilsteinInst