insulating states of non - trivial topological order have attracted much attention both theoretically and experimentally .
a topological insulator has a remarkable property of being metallic on the surface albeit insulating in the bulk .
recently much focus is on a specific type of topological insulators,@xcite which are said to be `` @xmath0-nontrivial '' .
@xcite the latter occurs as a consequence of interplay between a specific type of spin - orbit interaction and band structure .
@xcite such systems are invariant under time reversal and show kramers degeneracy . from the viewpoint of an experimental realization , the original idea of kane and mele ( km ) @xcite
was often criticized for being unrealistic , since it relies on a relatively weak spin - orbit coupling in graphene . in order to overcome such difficulty , bernevig , hugues and zhang ( bhz ) proposed an alternative system @xcite , which is also @xmath0-nontrivial but not based on graphene .
bhz model was intended to describe low - energy electronic properties of a two - dimensional ( 2d ) layer of hgte / cdte quantum well .
conductance measurement in a ribbon geometry @xcite showed that the system exhibits indeed a metallic surface state , which is also called _
helical _ edge modes .
this paper highlights the spectrum and wave function of such helical edge modes in the bhz model .
respecting appropriately the crystal structure of original hgte / cdte quantum well , one can safely implement it as a tight - binding model on a 2d square lattice.@xcite note , however , that in contrast to km model which can be represented as a purely lattice model , in bhz an internal spin-@xmath14 degree of freedom , stemming from the @xmath15-type @xmath16 and @xmath17-type @xmath18 orbitals , resides on each site of the square lattice in addition to the real electronic spin .
we also mention that in the continuum limit with vanishing topological mass term , km model has two valleys ( @xmath19 and @xmath20 ) , whereas bhz has a single valley ( at @xmath21 ) .
another idea which we can borrow from graphene study is the sensibility of edge spectrum on different types of edge structure , i.e. , either zigzag of armchair type .
@xcite this applies also to the helical edge modes of bhz topological insulator in a ribbon geometry , since the edge spectrum is predominantly determined by how the 2d bulk band structure is `` projected '' onto the 1d edge axis .
indeed , the structure of bhz helical edge modes has been extensively studied in ref .
@xcite in the @xmath7-edge geometry and in the tight - binding implementation .
however , in the practical experimental setup , @xcite this is certainly not the only one which is relevant to determine the transport characteristics at the edges . here , in this paper our main focus is on the other representative geometry , the @xmath11-edge . in the @xmath11-geometry , as a consequence of the specific way how `` hidden '' dirac cones ( or valleys ) in the 2d spectrum is projected onto the @xmath11-axis , edge modes show some unexpected behaviors . in order to motivate further the present study ,
let us first recall the importance of edge modes in the quantum hall state under magnetic field that exhibits a finite and quantized ( charge ) hall conductivity @xmath22 . in realistic samples with a boundary , quantization of hall conductivity is attributed to dissipationless transport due to a gapless chiral edge mode .
in contrast to charge hall effect , a finite _ spin _ hall current does not need breaking of the time reversal symmetry . in the quantized spin hall ( qsh )
effect , the chern number in the bulk takes an integral value , and correspondingly there appears integral pairs of gapless edge modes . on the other hand ,
the @xmath0 topological insulator is characterized by an odd number of gapless modes per edge that are robust against weak perturbations preserving the time - reversal symmetry .
the existence of such gapless edge states is generally guaranteed by a general theorem under the name of bulk / edge correspondence .
@xcite the bhz model has a convenient feature that the location of the minimum energy gap can be controlled by changing the parameters in the model .
in particular , the sign of spin hall conductivity changes discontinuously as the mass parameter of the model is varied .
hence , the corresponding change of edge spectrum should provide useful information on the bulk / edge correspondence .
furthermore , understanding of the nature of edge modes under specific edge geometries should serve as possible application of topologically protected phenomena in nano - architectonics .
we take the representative cases of the @xmath7- and @xmath11-edges , which we call also the straight and zigzag edges , respectively .
this paper is organized as follows : in sec .
ii , we clarify our motivations to study the lattice version of a @xmath0 topological insulator ( bhz model ) , implemented as a nearest - neighbor ( nn ) tight - binding model
. explicit form of the bhz lattice hamiltonian is introduced in sec .
it is demonstrated that by considering a lattice model , one naturally takes into account _ hidden _ dirac cones , the latter lacking in the analyses based on the continuum dirac model . in secs .
iv and v , we study the detailed structure of gapless edge modes under two different edge geometries : straight and zigzag edges .
v is the highlight of this paper , demonstrating that the zigzag edge modes of bhz lattice model show unique features .
we first point out that at @xmath23 , a pair of completely flat spectrum appears at @xmath24 ; besides the edge wave function can be trivially solved .
we then show , in a half - empirical way , that this analytic exact solution at @xmath23 can be generalized to the case of an arbitrary @xmath25 $ ] ( this idea is schematically represented in fig . [ zig_concept ] ) . using the exact solution
thus constructed , we highlight the nature of _ reentrant _ edge modes , another unique feature of the edge modes in the @xmath11-edge geometry .
the reentrant edge modes possess two contrasting characters in real and momentum spaces : though well distinguished in real space , they live in an extremely small energy scale in the spectrum .
vi is devoted to conclusions .
the gapless edge modes of bhz topological insulator are also treated in the framework of continuum dirac model in appendix a.
it has been well recognized that the quantized spin hall conductivity is determined by wave functions of bloch electrons over the entire brillouin zone . on the other hand ,
only the electronic states near the fermi level are relevant to the change of the hall conductivity when the fermi level is shifted .
simplified effective models are useful for the latter case since various theoretical techniques can be employed in the low - energy range . in this paper
, we work mainly with the lattice version of the bhz model , and make some comments in the low - momentum limit .      why do we have to go back to a lattice model ?
firstly , because we need to recover the _
correct _ absolute value of spin hall conductance @xmath26 .
the latter is defined as the difference of hall conductance for up and down ( pseudo ) spins multiplied by @xmath27 : @xmath28 in quantized spin hall ( qsh ) systems , the spin hall conductance is quantized to be an integer in units of @xmath29 . in the language of charge conductance . ]
this is completely in parallel with the quantization of _ charge _ hall conductance in units of @xmath30 in quantized hall system . in both cases ,
such integers are topological invariants and protected against weak perturbations .    on the other hand
, if one calculates , using kubo formula , the contribution of a single dirac fermion , e.g. , of the continuous dirac model at the @xmath21-point @xcite to spin hall conductance , then this gives _
half _ of the value expected from the topological quantization .
@xcite in order to be consistent , it is naturally assumed that there must be _ even _ number of dirac cones .
@xcite however , the low - energy effective theory with which we are starting contains obviously a single dirac cone .
@xcite as we will see explicitly , such trivial discrepancy is naturally resolved by considering a lattice version of the bhz model .
another aspect motivating us to employ the lattice version of bhz model is the fact the idea of an edge states is a real space concept , and we need _ a priori _ to go back to real space to give an unambiguous definition to it . in this paper
, we highlight the detailed structure of gapless edge modes under a specific edge geometry .
clearly , introduction of the latter needs a precise description in real space .
recall also here that edge modes of graphene nano - ribbon exhibit contrasting behaviors in zigzag and armchair edge geometries : e.g. , the system becomes either metallic or semi - conducting in the armchair geometry , depending on @xmath31 , with @xmath32 being the number of rows , whereas a completely flat edge mode appears in the zigzag geometry .
@xcite where does the difference comes from ? in momentum space , the question is how the bulk dirac cone structure look like when viewed from the edge .
note that in the zigzag geometry the flat edge mode connects the two dirac points : @xmath19 and @xmath20 , whereas in the armchair geometry these two points are projected onto the same point at the edge . in a topological insulator ,
this bulk to edge projection is even a more subtle issue , since not all the dirac cones are explicit ( see table i ) . in a sense , zigzag edge in the square lattice bhz model ( see fig .
[ geo_zig ] ) plays the following double role : it resembles the zigzag edge in graphene at @xmath33 , whereas it may rather resemble the armchair edge at @xmath34 and at @xmath35 .    [ cols="<,<,^,^,^,^,<,^,^",options="header " , ]     [ 4dpts ]      spin hall conductance is a topological quantity determined by the _ global _ structure of entire 2d brillouin zone . the helical edge modes , encoding the same information , is , therefore , _ a priori _ not derived from a local description in the 2d brillouin zone .
one exception to such a general idea is the study of ref .
@xcite ( see also appendix ) , in which gapless edge modes are derived from the continuous model in a strip geometry by simply applying the condition that all the pseudo spin components of wave function vanish at the boundary .
( @xmath3)-sublattice component of the wave function vanishes at one ( the other ) boundary . for details ,
see , ref . @xcite .
] this implies that information about the helical edge modes is actually encoded in the original ( single ) dirac cone .
this seems to be rather surprising , if one recalls that in the case of km model , the distinction between trivial and non - trivial phases is made by a relative sign of the mass gap at @xmath19- and @xmath20- points , which are , of course , macroscopically separated in momentum space . here
, in the bhz model , the same distinction is made by relative sign between the mass ( @xmath36- ) term and the @xmath37- term added to the dirac hamiltonian .
motivated by this observation , we investigate the structure of helical edge modes under different boundary conditions for the periodic bhz model , implemented as a square lattice and nearest - neighbor ( nn ) tight - binding model model .
in parallel with the arm chair and zigzag edges for graphene , we consider ( a ) usually considered @xmath7- ( straight ) boundary @xcite , as well as ( b ) @xmath11- ( zigzag ) boundary for the tight - binding bhz model .
the idea of focusing on dirac fermions in the description of quantized hall effect has appeared in the context of transitions between different plateaus . for describing the transitions ,
half - integer quantization is not a drawback , since the difference of hall conductance before and after the gap closing is quantized to be an integer in units of @xmath30 : @xmath38 or vice versa .
a discrete jump in the hall conductance across the transition is indeed consistent with counting based on the emergence of massless dirac fermions at the transition @xcite .
the absolute value of hall conductance is , on the other hand , a winding number , and determined by the vortices @xcite . here
, each vortex gives , in contrast to a dirac cone , an integral contribution to the hall conductance ( in units of @xmath30 ) .
an interesting question is whether the total hall conductance , often expressed as a topological invariant @xcite , can be also written as a sum of contributions from emergent dirac electrons in the spectrum .
our empirical answer is yes,@xcite indicating that the number of dirac electrons emergent in the spectrum is always even , reminiscent of the no - go theorem of nielsen and ninomiya in 3 + 1 dimensions @xcite .
it should be noted that here not only explicit dirac electrons ( gapless for a given set of parameters ) but also hidden dirac electrons ( massive for that set of parameters ) must be taken into account .
such massive dirac electrons are called `` spectators '' in ref .
@xcite , in the sense that they are inactive for the transition .
spectators are indispensable to ensure the correct integer quantization of the hall conductance .
( vertical axis ) of a bhz lattice model , here , implemented as a nn tight - binding model on a square lattice : cf .
( [ dk_nn ] ) .
the spectrum is shown over the entire brillouin zone : @xmath39 ( horizontal plane ) .
parameters are chosen such that @xmath40 and @xmath41 , i.e. , the system is in the topologically non - trivial phase : @xmath42 .
the spectrum shows a typical wine - bottle structure around the @xmath21-point.,width=302 ]
let us first consider the bhz model in the long - wave - length limit .
the low - energy effective hamiltonian , describing the vicinity of gap closing at @xmath43 , is the minimal model to capture the physics of a @xmath0-topological insulator .
this effective hamiltonian is also contrasting to the prototypical km model , in that the former describes only a single dirac cone .
the distinction between the @xmath0- trivial and non - trivial phases is , therefore , made by adding a @xmath37-term to the usual dirac hamiltonian .
the explicit form of bhz hamiltonian is implemented by the following @xmath44 matrix : @xmath45 , \label{htot}\end{aligned}\ ] ] where @xmath46 is a 2d crystal momentum , here measured from the @xmath21-point .
the lower - right block @xmath47 is deduced from @xmath48 by imposing time reversal symmetry .    the bulk energy spectrum
: @xmath49 is then given by solving the eigenvalue equation for the upper - left block @xmath48 of the @xmath44 bhz hamiltonian , i.e. , @xmath50 in order to represent @xmath48 in a compact form , we introduce a @xmath51-vector , @xmath52 , each component of which is either an even or an odd function of @xmath53 : @xmath54 , whose parity is determined by symmetry considerations .
@xcite as far as the low - energy universal properties in the vicinity of @xmath21-point is concerned , we need to keep only the lowest order terms of @xmath55 , and in this long - wave - length limit , they read explicitly , @xmath56 other parameters which appear in ref .
@xcite , i.e. , @xmath57 and @xmath58 are set to be zero , which , however , does not lose any generality .
the bulk energy spectrum is thus determined by diagonalizing the following `` spin hamiltonian '' , @xmath59 , where @xmath60 , are pauli matrices . using the standard representation for @xmath61 , @xmath62 reads explicitly as , @xmath63 .
\label{hk}\ ] ] each row and column of eq .
( [ hk ] ) represent an `` orbital spin '' associated with the @xmath15-type @xmath16 and the @xmath17-type @xmath18 orbitals of the original 3d band structure of hgte and cdte .
@xcite then , by choosing the `` spin quantization axis '' in the direction of @xmath51-vector , one can immediately diagonalize the hamiltonian @xmath48 , i.e. , @xmath64 where the eigenvalue @xmath65 given by , @xmath66 this implies , @xmath67 where @xmath68 .
when @xmath69 , the dispersion relation ( [ wine ] ) represents a wine - bottle structure ( fig .
2 ) , i.e. , @xmath70 shows a minimum at a finite value of @xmath71 . at the critical value @xmath72 , the density of states shows van hove singularity .     in units of @xmath73 , obtained by numerically evaluating the @xmath53-integral in eq .
( [ tknn ] ) , plotted as a function of @xmath8.,width=8 ]      lattice version of the bhz model is implemented as a tight - binding hamiltonian .
to construct such a hamiltonian explicitly , we replace linear and quadratic dependences in @xmath48 on @xmath74 and @xmath75 as in eqs .
( [ dk ] ) , by a function which has the right periodicity of the square lattice .
this can be implemented as , @xmath76 , \label{dk_nn}\end{aligned}\ ] ] where @xmath77 is the lattice constant .
( [ dk_nn ] ) corresponds to regularizing the effective dirac model on a square lattice with only nearest - neighbor ( nn ) hopping . in this setup , i.e. , eqs .
( [ htot ] ) , ( [ hk ] ) together with eqs .
( [ dk_nn ] ) , the lattice version of bhz model acquires four gap closing points shown in table i , if one allows the original mass parameter @xmath2 to vary beyond the vicinity of @xmath78 .
the new gap closing occurs at different points in the brillouin zone from the original dirac cone ( @xmath21-point ) , namely at @xmath79 , @xmath80 and @xmath81 .
the gap closing at @xmath82 occurs at @xmath83 , whereas the gap closing at @xmath84 and @xmath85 occurs simultaneously when @xmath86 .
each time a gap closing occurs , one can re - expand the lattice model with respect to small deviations of @xmath53 measured from the gap closing .
the new effective model in the vicinity of such hidden gap closing falls on the same dirac form as the original one at the @xmath21-point , up to the @xmath37-term . in order to quantify the emergence of such hidden dirac cones ,
one still needs the following two parameters : ( i ) the mass gap @xmath2 ( especially , its sign ) , and ( ii ) the chirality @xmath87 . the latter is associated with the homotopy in the mapping : @xmath88 .
note that in the gap closing at @xmath84 and @xmath85 at @xmath86 , the role of @xmath89 is interchanged compared with the original dirac cone at the @xmath21-point .
the former ( latter ) corresponds to @xmath90 ( @xmath91 ) .
the missing dirac partner , in the sense of ref .
@xcite , is found in this way .
once the explicit form of effective dirac hamiltonian in the continuum limit is given , one can determine its contribution to @xmath92 . in systems with trs ,
i.e. , of the form ( [ htot ] ) , the contribution from @xmath48 to @xmath93 cancels with that of @xmath94 . on the other hand ,
their contribution to @xmath95 , remains finite , takes a half - integral value , @xmath96 in units of @xmath73 ( @xmath97 in conductance ) .
contribution to @xmath98 from a dirac point with a mass gap @xmath2 and chirality @xmath87 is , @xmath99 provided that the fermi energy is in the gap .
this can be verified explicitly by applying the kubo formula to the continuum model .    as mentioned earlier , such counting based on the continuum dirac model ,
is known to describe correctly a discrete jump of @xmath100 in the quantum hall case . here , we apply the same logic to qsh case .
imagine that one observes the evolution of @xmath98 , starting with the trivial insulator phase , where @xmath101 , and varying the mass parameter @xmath2 .
the fermi energy is always kept in the gap unless there appears a dirac cone .
each time such a gap closing occurs , @xmath98 shows a discrete change , which can be attributed to the above dirac fermion argument . in the present model
, one can verify explicitly that this is indeed the case .
if one evaluates the spin hall conductance from tknn formula , @xcite @xmath102 allows for the following representation , in terms of berry curvature integrated over the entire brillouin zone : @xmath103 where @xmath104 and @xmath105 is given , e.g. , by eq .
( [ dk_nn ] ) .
for such an explicit choice of @xmath106 , eq . ( [ tknn ] ) is evaluated numerically , and plotted as a function of @xmath107 in fig .
[ sigma_xy ] . when @xmath106 is given by eq .
( [ dk_nn ] ) , the plotted curve ( the solid curve shown in blue in fig . [ sigma_xy ] , which looks practically like steps ) is comparable with the column @xmath108 of table i. note that the absolute value of @xmath102 is susceptible of the concrete implementation of @xmath106 over the entire brillouin zone , whereas its parity ( whether it is even or odd ) in units of @xmath109 remains the same . as well known , the latter determines system s @xmath0-property .
@xcite    we have seen that @xmath102 takes a finite value @xmath110 when @xmath111 , i.e. , the system is in the topological ( inverted gap ) phase .
this is also consistent with the gapless edge picture , in which the spin hall conductance of twice the unit of quantum conductance @xmath30 is attributed to two channels of edge modes , which form a pair of kramers partners .
the apparent half - integer quantization at the @xmath21-point , in the sense of eq .
( [ 1/2 ] ) , is compensated by the contribution from missing dirac partner(s ) , and as a result , is indeed shifted by one - half , replaced by an expected integral quantization .
edges ) are , here , chosen to be perpendicular to the @xmath112-direction . in this figure ,
the number of rows in the strip is @xmath113.,width=6 ]
let us first review the behavior of gapless edge modes in the straight edge geometry , the latter commensurate with the square lattice , and can be chosen either normal to the @xmath7- or to the @xmath112-direction ( as in fig .
[ geo_str ] ) . introducing
an edge leads to breaking of the translational invariance in the direction perpendicular to the edge , inducing a coupling between dirac cones .                  in the straight edge geometry shown in fig .
[ geo_str ] ) , electrons are confined inside a strip between the rows at @xmath114 and @xmath115 . the translational invariance along the @xmath116-axis is still maintained , allowing for constructing a 1d bloch state with a crystal momentum @xmath74 : @xmath117 where @xmath74 is measured in units of @xmath118 with @xmath77 being the lattice constant .
@xmath119 is a one - body electronic state localized on site @xmath120 , and @xmath121 is an operator creating such an electron .
it is also convenient to introduce @xmath122 , and express @xmath123 as @xmath124 .
naturally , the two creation operators are related by fourier transformation similarly to eq .
( [ bloch ] ) , i.e. , @xmath125 .    in order to introduce the edges , it is convenient to rewrite the bhz tight - binding hamiltonian in terms of the hopping between neighboring rows
. let us first consider the bhz tight - binding hamiltonian in real space : @xmath126 , \label{ham_bhz}\end{aligned}\ ] ] where @xmath127 hopping matrices , @xmath128 and @xmath129 , are given explicitly as , @xmath130 in order to rewrite it in terms of the bloch state , eq .
( [ bloch ] ) , or equivalently , in terms of the corresponding creation and annihiration operators , @xmath122 ( @xmath131 ) , we perform fourier transformation in the ( @xmath116- ) direction along the edge .
( [ ham_bhz ] ) thus rewrites , @xmath132 \label{ham_str}\ ] ] where @xmath133 s are diagonal ( on - row ) components , which read explicitly , @xmath134 @xmath135 represents a hopping amplitude in the @xmath136-direction , i.e. , between neighboring rows . inside the strip ,
i.e. , for @xmath137 , these amplitudes take the same value as in the bulk , given in eqs .
( [ gamma_xy ] ) , i.e. , @xmath138    in the tight - binding implementation , a strip geometry can be introduced by switching off all the hopping amplitudes connecting sites on the edge of the sample to the exterior of the sample . in our straight edge geometry ,
such outermost rows are located at @xmath139 and @xmath140 .
we turn off all the hopping amplitudes from @xmath139 to @xmath141 , and the ones from @xmath142 to @xmath143 , i.e. , @xmath144 this boundary condition , ( i ) breaks the translational invariance in the @xmath136-direction , and ( ii ) restricts the hamiltonian matrix into @xmath145 blocks .
and @xmath146 plotted as a function @xmath8 in the limit @xmath147.,width=8 ]      let us construct the eigenvector of the straight edge hamiltonian , eq .
( [ ham_str ] ) , with an eigenenergy @xmath148 .
( [ ham_str ] ) is already diagonal w.r.t .
@xmath74 , we diagonalize eq .
( [ ham_str ] ) for a given @xmath74 , to find the energy spectrum @xmath149 .
the corresponding eigenvector is thus specified by @xmath148 and @xmath74 , and takes generally the following form : @xmath150 where @xmath151 is a @xmath127 spinor specifying the amplitude and the pseudo spin state of eigenvector on row @xmath152 .
one might rather regard , @xmath153 , \label{wf}\ ] ] as the wave function of the corresponding eigenstate . the eigenvalue equation , @xmath154 can be rewritten , in terms of the @xmath155 s , in the form of a recursive equation : @xmath156    all the information on the spectrum and the wave function of both the extended bulk states and the localized edge states is encoded in eq .
( [ rec_str ] ) and the boundary condition which we will specify later . since the recursive relation , eq . ( [ rec_str ] ) ,
is linear , its eigenmodes take the form of a geometric series : @xmath157 where @xmath158 is a solution of the characteristic equation which we will derive later . if @xmath159 , eq . ( [ ser_rho ] ) may represent an edge mode .
since the recursive relation , eq . ( [ rec_str ] ) , is of second order , its characteristic equation becomes a quadratic equation , giving two solutions for @xmath158 , say , @xmath160 . on the other hand ,
our recursive equation has also a @xmath127 matrix form , we first have to solve a ( reduced ) eigenvalue equation for @xmath161 , assuming that @xmath158 is given .
the reduced eigen value equation for @xmath161 reads , @xmath162 \psi_0 = e \psi_0.\ ] ] using eqs .
( [ gamma_xy ] ) this can be also rewritten as , ) corresponds to the choice , @xmath163 , or @xmath164 which is consistent with the bloch theorem . in the strip geometry with the periodic boundary condition
, @xmath75 takes discrete values . with the open boundary condition relevant to the actual strip
, @xmath75 is no longer a good quantum number .
however , if the width of the strip is large , one may roughly interpret the 1d energy spectra in the strip geometry as composed of the many slices of bulk energy spectrum at different values of @xmath75 .
in addition , a pair of edge modes appear as a characteristic feature of the nontrivial topological property . ]
@xmath165 \psi_0 = e \psi_0 , \label{eeq_red}\ ] ] where @xmath166 is given in eq .
( [ dkx ] ) .
this is a @xmath127 eigenvalue equation , and there are generally two solutions for @xmath148 and two corresponding eigenvectors for a given @xmath158 .
recall here that for @xmath167 , our numerical data ( fig . [ spec_str ] ) show that the edge spectrum behaves as @xmath168 in the limit of @xmath147 .
hereafter , we will focus only on such edge solutions .
since in the same limit , @xmath169 , eq .
( [ eeq_red ] ) reduces to , @xmath170 \psi_0 = 0 . \label{eeq_0}\ ] ] note that we have multiplied both sides of eq .
( [ eeq_red ] ) by @xmath171 .
it is clear from this expression that @xmath161 can be chosen to be an eigenstate of @xmath172 , i.e. , @xmath173 , where @xmath174,\ \
\left [ \begin{array}{r } 1\\ -1 \end{array } \right ] .
\label{x+-}\ ] ] if one denotes the eigenvalue of @xmath172 by @xmath15 , as @xmath175 , then @xmath176 corresponds to @xmath177 , and @xmath178 to @xmath179 .
namely , @xmath15 specifies the eigenspinors given in eqs .
( [ x+- ] ) . based on these eigenspinors
, one can construct the total wave function .
of course , one still needs to do determine the allowed values of @xmath158 .
for an eigenstate specified by @xmath15 , @xmath158 must satisfy , @xmath180 for @xmath176 , @xmath181 in eq .
( [ quad ] ) , if @xmath158 is a solution of this quadratic equation for @xmath176 , then @xmath182 satisfies the same equation for @xmath178 . @xcite
thus the general solution becomes a linear combination of the following four basic solutions : @xmath183
|x+ \rangle   \nonumber \\ & + & % .\left [ c_{-1 } \left(1 \over \rho_1 ( 0 ) \right)^{j } +   \left [ c_{-1 } \rho_1 ( 0 ) ^{-j } + c_{-2 } \rho_2 ( 0 ) ^{-j }    \right ] |x- \rangle , \label{sol_0}\end{aligned}\ ] ] of course , at this point , this is just a solution at only one single @xmath184-point , @xmath185 . however , a solution in the form of eq .
( [ sol_0 ] ) , with @xmath186 given in eq .
( [ rho12_0 ] ) , can be easily generalized to satisfy eq .
( [ rec_str ] ) for an arbitrary , finite @xmath74 , after a simple replacement of parameters . as for the eigenmode of the form of eq .
( [ ser_rho ] ) , one has to solve a reduced eigenvalue equation for @xmath187 , and determine @xmath158 such that eq .
( [ eeq_red ] ) is satisfied .
however , since @xmath133 has the structure given in eq .
( [ dkx ] ) , the eigenmodes for @xmath187 given as eq .
( [ x+- ] ) remain to be valid for an arbitrary @xmath74 .
this might become clearer , if one decomposes eq .
( [ eeq_red ] ) into the following set of equations : @xmath188 \psi_0 & = & 0 , \label{eeq_2}\end{aligned}\ ] ] i.e. , eq .
( [ eeq_red ] ) is recovered by adding both sides of eqs .
( [ eeq_1 ] ) and ( [ eeq_2 ] ) .
the first equation gives the energy dispersion , @xmath189 , if @xmath187 is chosen to be an eigenstate of @xmath172 , i.e. , @xmath190 , where @xmath191 note that this is an exact edge spectrum valid over the entire range of @xmath74 , as far as the edge solution is possible ( see fig .
[ spec_str ] ) . on the other hand , eq .
( [ eeq_2 ] ) , analogous to eq .
( [ eeq_0 ] ) , justifies the previous conjecture : @xmath192 . while , in the characteristic equation for @xmath158 , one has to make the simple replacement : @xmath193 , i.e. , eqs . ( [ eeq_2 ] ) and ( [ eeq_0 ] )
are identical up to this replacement .
for @xmath176 , the solution for @xmath158 reads , @xmath194 correspondingly , a general solution for @xmath195 can be constructed as , @xmath196 |x+ \rangle   \nonumber \\ & + & \left [ c_{-1 } \rho_1 ( k_x ) ^{-j } + c_{-2 } \rho_2 ( k_x ) ^{-j } \right ] |x- \rangle , \label{sol_gen}\end{aligned}\ ] ] where the coefficients @xmath197 should be chosen to satisfy the boundary conditions .
( [ sol_gen ] ) is smoothly connected to eq .
( [ sol_0 ] ) in the limit : @xmath147 .
three panels of fig .
[ spec_str ] show the energy spectrum ( edge + bulk ) for different values of @xmath8 .
as for the edge part of the spectrum , only a part of eq .
( [ spec_th_str ] ) is realized . in order to determine which part of the spectrum in eq .
( [ spec_th_str ] ) is indeed activated , we discuss below the case of semi - infinite geometry in some detail .
[ spec_str ] also demonstrates one of another specific feature of straight edge mode that the main location of the mode moves from the zone center for @xmath198 , to the zone boundary for @xmath199 .
thus , the group velocity intersecting with the fermi level reverses its sign , reflecting the sign change of @xmath200 in the bulk .
this can be regarded as the concrete expression of bulk / edge correspondence in the present case .
@xcite    what kind of a boundary condition should we apply in eq .
( [ sol_gen ] ) ?
suppose that here our system is _ semi - infinite _
, for simplicity , extended from @xmath201 to @xmath202
. such a boundary condition can be applied , by formally requiring that the wave function ( [ sol_gen ] ) vanishes at @xmath203 , i.e. , @xmath204.\ ] ] this means that the coefficients @xmath197 in eq .
( [ sol_gen ] ) should be chosen to satisfy , @xmath205 this turns out to be rather an important requirement for determining the range of validity of the solution given in eq .
( [ sol_gen ] ) , since the wave function @xmath195 must be normalizable . in eq .
( [ sol_gen ] ) , only the eigenmodes of the form of eq . ( [ ser_rho ] ) with @xmath159 should be kept in the solution ( to be precise , _ both _ @xmath206 and @xmath146 must be smaller than 1 ) . in a strip geometry , another solution , consisting of both @xmath207 ,
describes the edge mode localized at the opposite end of the system .
and @xmath146 plotted as a function of @xmath74 at @xmath208 ( @xmath40 ) .
reference lines are at @xmath209 , and at @xmath164 and @xmath210  as for the latter , cf .
( [ rho_flat ] ) .
@xmath164 corresponds indeed to the point at which the edge modes merge with the bulk spectrum ( central panel , @xmath211 ) .
lower panel : an enlarged image of the central panel , plotted also for a system of larger size : @xmath212.,title="fig:",width=8 ]   and @xmath146 plotted as a function of @xmath74 at @xmath208 ( @xmath40 ) .
reference lines are at @xmath209 , and at @xmath164 and @xmath210  as for the latter , cf .
( [ rho_flat ] ) .
@xmath164 corresponds indeed to the point at which the edge modes merge with the bulk spectrum ( central panel , @xmath211 ) .
lower panel : an enlarged image of the central panel , plotted also for a system of larger size : @xmath212.,title="fig:",width=8 ]    and @xmath146 plotted as a function of @xmath74 at @xmath208 ( @xmath40 ) .
reference lines are at @xmath209 , and at @xmath164 and @xmath210  as for the latter , cf .
( [ rho_flat ] ) .
@xmath164 corresponds indeed to the point at which the edge modes merge with the bulk spectrum ( central panel , @xmath211 ) .
lower panel : an enlarged image of the central panel , plotted also for a system of larger size : @xmath212.,title="fig:",width=8 ]    in fig .
[ sol_str1 ] , @xmath206 and @xmath146 are plotted as a function of @xmath8 in the limit @xmath147 . when @xmath213 , both @xmath206 and @xmath146 are _ smaller _ than 1 , namely both @xmath214 and @xmath215 are _ larger _ than 1 .
this means that only the first two terms of eq .
( [ sol_gen ] ) , both corresponding to @xmath216 ( @xmath176 ) , should be kept in the solution , i.e. , @xmath217 outside this region , either @xmath206 or @xmath146 is larger than 1 . when @xmath207 , since this implies automatically @xmath218 , the eigenmode corresponding to the latter is still compatible with the boundary condition at @xmath219 .
however , because of the boundary condition at @xmath203 , i.e. , eq .
( [ req_bc ] ) , when @xmath220 and @xmath221 , or vice versa , the only possible choice for the coefficients @xmath197 is @xmath222 namely , a solution of the type of eq .
( [ sol_gen ] ) , or an edge mode crossing at @xmath185 is _
inexistent_. this is consistent with the fact that an edge mode crossing at @xmath185 exists only in the region , @xmath223 $ ] in the straight edge geometry .
( when @xmath224 $ ] , the edge modes cross at @xmath225 , i.e. , at the zone boundary , which is also time - reversal symmetric . )    coming back to the regime in which edge modes are existent , i.e. , @xmath226 $ ] , one can clearly see in fig .
[ sol_str1 ] that there are two different behaviors  a flat region where @xmath206 and @xmath146 are degenerate , and the remaining part with two branches .
this is due to the fact that the two solutions for @xmath158 could be either both real , or a pair of complex numbers conjugate to each other . in the latter case ,
the two solutions have the same absolute value , @xmath227 , whereas in the present case , one can verify that this degenerate value is independent of @xmath2 , i.e. , @xmath228 this explains the existence of a flat region in fig .
[ sol_str1 ] . from eq.([rho12_0 ] ) we see that the square root becomes pure imaginary for all @xmath184 provided @xmath229 , where appear in the regime : @xmath230 , and also at the other end . in the approximation that becomes valid in the small wave number
, this threshold value is given by @xmath231 .
for @xmath40 , @xmath232 , whereas , @xmath233 is , of course , @xmath234 . if one expands eq .
( [ delta_pm ] ) in powers of @xmath235 , then at leading order @xmath236 coincides with @xmath233 . ] @xmath237    at @xmath164 , the edge solution is expected to merge with the bulk spectrum .
this happens , when @xmath238 such a behavior becomes clearer by plotting @xmath239 s as a function of @xmath74 .
[ sol_str2 ] illustrates this feature at @xmath41 for @xmath40 .
one can indeed see that the edge spectra merge with the bulk at @xmath240 , satisfying eq .
( [ km ] ) .
the latter reduces , at this value of @xmath2 , to @xmath241 , i.e. , @xmath242 .    it is also instructive to investigate the nature of edge modes in _ real _ space , i.e. , the wave function , and compare it with the general solution ( [ sol_gen ] ) . in numerical experiments , one has to diagonalize the @xmath243 hamiltonian matrix , equivalent to eq .
( [ ham_str ] ) .
an eigen wave function is , therefore , obtained as a @xmath244-component vector ; here , in the straight edge geometry , the latter can be chosen to be real .
the edge wave function is easily identified if it exists , e.g. , by choosing the lowest - energy eigenmode @xmath161 in the upper band . by investigating the structure of such an edge wave function
, one can explicitly verify that the eigenmodes are spanned by two eigenspinors given in eq .
( [ x+- ] ) . in repeating such numerical experiments for different @xmath74 and @xmath8
, one can naturally distinguish an edge state from a bulk state by focusing on the spatial distribution of the wave function . here
, what deserves much attention is that one can recognize a one - to - one correspondence between localizability of the wave function @xmath161 and its spinor structure .
let us turn to the case of a different edge geometry , the zigzag edge geometry , shown schematically in fig .
[ geo_zig ] . as mentioned earlier
, the zigzag edge geometry considered here is , in a sense , analogous to a more popular edge geometry of graphene ribbon , named in the same way , but defined on a hexagonal lattice . here , on a square lattice ,
a zigzag edge is introduced , either normal to @xmath11- or @xmath245-direction ( as in fig .
[ geo_zig ] ) .
electrons in the zigzag edge geometry are , therefore , confined to a strip diagonal in the cartesian coordinates , say , along the @xmath11-direction as in fig .
[ geo_zig ] .
-edge geometry ( here , chosen to be normal to the @xmath245-axis ) , defined in terms of the original square lattice , on which the new @xmath246- and @xmath247-axes are superposed in blue ( online ) .
numbers on these axes are redefined indices @xmath248 and @xmath249 , i.e. , @xmath250 . the number of rows @xmath32 is here chosen to be , @xmath251 .
when @xmath32 is even ( odd ) , the two edges are inversion asymmetric ( symmetric ) w.r.t .
the center of strip.,width=6 ]    intuitively , say , because the zigzag surface is literally , `` rough '' , one expects that this edge geometry might have a stronger tendency to trap electrons in the vicinity of the boundary .
we show below , on one hand , that this intuition from the macroscopic world is still valid in the microscopic quantum mechanical world .
but just as a result of this stronger tendency to keep the electrons in its vicinity , on the other hand , the zigzag edge shows various curious phenomena , e.g. , completely flat edge modes , and the reentrance of edge modes in @xmath184-space , etc .
clearly , the translational invariance along the @xmath11-axis is maintained , on which @xmath246-axis is introduced , together with the conserved momentum @xmath252 in this direction .
accordingly , @xmath247-axis is chosen to be in the @xmath245-direction .
it may be also useful to _ redefine _ the indices @xmath248 and @xmath249 such that the lattice points in the original square lattice are located at @xmath250 , where @xmath248 is an even ( odd ) integer for @xmath249 : even ( odd ) . in the zigzag edge geometry , the spectrum @xmath253 and the wave function @xmath254 are determined by the following recursive equation : @xmath255 for @xmath195 , analogous to eq .
( [ rec_str ] ) in the straight edge geometry . in order to derive eq .
( [ rec_zig ] ) , we first rewrote the tight - binding hamiltonian ( 15 ) in the new labeling , and then considered a bloch state along the @xmath246-axis , analogous to eq .
( [ bloch ] ) , but with a crystal momentum @xmath184 conjugate to @xmath246 , i.e. , @xmath252 . as was the case in eq .
( [ rec_str ] ) , @xmath21 describes , here , in eq .
( [ rec_zig ] ) the hopping between adjacent rows , and reads explicitly as , @xmath256 note that here the crystal momentum @xmath184 is measured in units of @xmath257 so that the zone boundary is always given by @xmath258 . and @xmath259 .
] to find an edge solution , we first express the solution of eq .
( [ rec_zig ] ) , in the form of a geometric series , written formally in the same as eq .
( [ ser_rho ] ) . recall that @xmath158 is ( generally ) a ( complex ) number of , for an edge mode , amplitude smaller than unity ( with the understanding that the edge mode is localized in the vicinity of @xmath203 ) .
@xmath187 is a two component eigenvector of the following reduced eigenvalue equation : @xmath260 \psi_{0 } = e \psi_0 . \label{reduced}\ ] ] from the analogy to the straight edge case , one may express @xmath21 as @xmath261 and rewrite eq .
( [ reduced ] ) into the following explicit form : @xmath262 \psi_{0 } = e \psi_0 . \label{}\end{aligned}\ ] ] where we @xmath263 and @xmath264 are short - hand notations for , respectively , @xmath265 and @xmath266 . comparing this form with the straight edge case
, one can see that here one can not use the same recipe for solving the problem , i.e. , solving the problem at , say , @xmath267 and then extrapolate its solution to general @xmath184 .
it seems not impossible to proceed in that direction and solve the problem analytically , but here , we choose to take another route , which is much simpler , to find still an exact solution , but in a half - empirical way .
+                  some concrete examples of such energy spectrum are shown fig .
[ spec_zig ] .
a pair of gapless edges modes always appear iff @xmath42 .
in contrast to the straight edge case , however , they appear always in the vicinity of @xmath267 , and intersects at @xmath267 .
one of the last panels of fig .
[ spec_zig ] shows a unique feature of edge modes in the zigzag edge geometry . at @xmath86 ,
the edge modes become completely flat , apart from a small finite - size gap around the zone boundary .
we have already seen such flat edge modes in graphene in the case of zigzag edge geometry ( but on a hexagonal lattice ) .
@xcite in graphene , such flat edge modes connect 1d projection of @xmath19 and @xmath20 points via the 1d bz boundary .
this is a similar behavior to the present case , if one regards the former as the limit of vanishing intrinsic coupling ( or topological mass @xmath2 ) in the km model .
one of the differences between the two cases is that here the flat edge modes cover the entire 1d brillouin zone .    in order to elucidate the nature of flat edge mode at @xmath33 , first notice that at this value of @xmath2 the diagonal terms of ( the diagonal blocks of the hamiltonian matrix in ) eq . ( [ rec_zig ] ) vanish
this implies , as in graphene nano - ribbon in the zigzag edge geometry , the existence of an eigenstate of the form : @xmath268 \label{psi_2j=0}\ ] ] i.e. , an eigenvector satisfying @xmath269 ( @xmath270 ) . here ,
a semi - infinite geometry is implicit ( for approximating a ribbon of sufficiently large width or @xmath32 ; our system extended from @xmath201 to @xmath271 ) , with a boundary condition of @xmath272 . under this setup , and with the condition of vanishing diagonal matrix elements , eq .
( [ rec_zig ] ) implies , @xmath273 simultaneously , it should also have a vanishing eigenenergy @xmath24 for consistency .    clearly , eq . ( [ rec_d4b ] )
has a solution of the form of a geometric series , here , for @xmath274 ( @xmath275 ) : @xmath276 note that here @xmath277 plays , roughly , the role of @xmath278 , but their precise relation will become clearer when the entire problem is solved . in order to proceed , we recall that @xmath21 can be written explicitly as , @xmath279 . \label{gamma_zig_mat}\ ] ] then , by assuming a solution of the form of eq .
( [ ser_lambda ] ) , eq .
( [ rec_d4b ] ) can be reduced to the following eigenvalue problem for @xmath280 : @xmath281 with the eigenvalues , @xmath282 and the corresponding eigenvector , @xmath283 , i.e. , @xmath284 given explicitly as , @xmath285 .
\label{u+-}\ ] ] the coefficient @xmath286 is a function of @xmath184 , which takes precisely the following form : @xmath287 notice that @xmath288 , and the two eigenspinors @xmath283 are orthogonal : @xmath289 . a general solution in the form of eq .
( [ psi_2j=0 ] ) can be thus constructed by applying @xmath290 , recursively , to @xmath291 and the result is , @xmath292 + c_- \lambda_-^j \left [ \begin{array}{c } \alpha_- ( 1- i)\\ 1 \end{array } \right ] .
\label{sol_d4b}\ ] ] in this construction , the two eigenvectors @xmath283 always have the form of eq .
( [ u+- ] ) .
this feature remains when @xmath2 is away from @xmath293 at which the edge modes are no longer completely flat , or rather even in the regime in which the edge spectrum is not flat at all .
another remark , concerning the behavior of eq .
( [ sol_d4b ] ) is that under the choice of signs in eq .
( [ lambda+- ] ) , @xmath294 is always larger than 1 , whereas @xmath295 except at the zone boundary .
this can be easily verified either numerically , i.e. , by plotting @xmath296 as a function of @xmath184 , or by showing @xmath297 using directly the expression for @xmath296 in eq .
( [ lambda+- ] ) . in numerical simulation for systems of a finite number of rows ,
both of these two solutions play a role giving rise to a pair of edge solutions .
is odd , one solution , corresponding , say , to @xmath298 , is localized in the vicinity of @xmath201 , and indeed has the form of eq .
( [ psi_2j=0 ] ) with @xmath299 given as eq .
( [ sol_d4b ] ) and @xmath300 .
the other solution , corresponding to @xmath301 , has the same structure of eq .
( [ psi_2j=0 ] ) but with @xmath299 increasing with practically an equal geometric ratio of @xmath301 , and naturally localized in the vicinity of the other edge : @xmath271 . on the other hand ,
when @xmath32 is even , the eigenmode of the system becomes a linear combination of the above two types of solutions .
this even / odd feature occurs only at @xmath2 precisely equal to @xmath293 , since at this value of @xmath2 where the edge spectrum becomes completely flat , the two ( generally ) counter - propagating edge modes acquire the same ( zero ) group velocity , and get mixed . ]     plotted as a function of @xmath184 at @xmath302 but for different values of @xmath3 , i.e. , @xmath303 and @xmath304 , corresponding , respectively , the colors : red , green , blue , cyan , magenta and orange ( upper panel ) .
edge spectrum at @xmath86 and @xmath305 ( @xmath302 , @xmath212 , lower panel ) .
two reference lines are at @xmath306 , the value of @xmath184 at which @xmath298 vanishes at @xmath305 .
the axes are suppressed in the lower panel so as to highlight the completely flat edge modes.,title="fig:",width=8 ]   plotted as a function of @xmath184 at @xmath302 but for different values of @xmath3 , i.e. , @xmath303 and @xmath304 , corresponding , respectively , the colors : red , green , blue , cyan , magenta and orange ( upper panel ) .
edge spectrum at @xmath86 and @xmath305 ( @xmath302 , @xmath212 , lower panel ) .
two reference lines are at @xmath306 , the value of @xmath184 at which @xmath298 vanishes at @xmath305 .
the axes are suppressed in the lower panel so as to highlight the completely flat edge modes.,title="fig:",width=8 ]    in fig .
[ lambda ] , the upper panel shows @xmath298 plotted as a function of @xmath184 at @xmath302 but for different values of @xmath3 , naturally assuming @xmath86 .
when @xmath3 is smaller than a critical value @xmath307 @xmath298 changes its sign ( has a zero ) at intermediate @xmath184 .
the edge wave function is extremely localized , i.e. , onto a single row : @xmath201 or @xmath271 .
] this continues to be the case even in the limit @xmath3 vanishes .
the lower panel shows the edge spectrum at @xmath86 , @xmath305 and @xmath302 .
note that the zero of @xmath298 corresponds to the value of @xmath184 at which the bulk spectrum focuses onto a single point .      as we will describe in detail in the next subsection ,
our `` recipe '' for constructing the _ exact _ edge wave function , and simultaneously its spectrum , lies in `` extrapolating '' the exact solution available at @xmath86 to a general value of @xmath308 $ ] ( recall also fig .
[ zig_concept ] ) . to complete this program
, we need to refer to some results of the numerical experiments performed for a system of finite number of rows . we have already seen the spectrum of such systems in fig .
[ spec_zig ] ; here we focus on the behavior of wave function , i.e. , the behavior of @xmath195 as a function of @xmath152 .    in the zigzag edge geometry , it is remarkable that ( not only ) the edge wave function ( but also the bulk wave function ! ) has the following particular form : @xmath309 .
\label{evec_gen}\ ] ] when that eigenstate represents an edge mode , eq .
( [ evec ] ) further simplifies : @xmath310 \vspace{0.2 cm } \\ c_4 \left [ \begin{array}{c } \alpha ( 1-i ) \\ 1 \end{array } \right ] \vspace{0.2 cm } \\ c_6 \left [ \begin{array}{c } \alpha ( 1-i ) \\ 1 \end{array } \right ] \vspace{0.2 cm } \\ \vdots
\\ \end{array } \right ] , \label{evec_edge}\ ] ] i.e. , for a given set of parameters @xmath311 , @xmath3 ( and @xmath2 ) as well as for a fixed @xmath184 , @xmath312 is a constant ( @xmath313 ) .
the ratio , on the other hand , @xmath314 is a measure of , to what extent the edge mode is localized in the vicinity of a boundary , say , at @xmath201 .
( @xmath275 ) plotted ( red points with filling to the horizontal axis ) at @xmath315 for @xmath316 ( upper panel ) and @xmath317 ( lower panel ) .
@xmath40 , @xmath318 . in the upper panel
, the blue line corresponds to the `` theoretical '' value , @xmath319 , whereas in lower panel , the plots are fitted by a curve , @xmath320 , with the choice of parameters , @xmath321 , @xmath322,title="fig:",width=8 ]    ( @xmath275 ) plotted ( red points with filling to the horizontal axis ) at @xmath315 for @xmath316 ( upper panel ) and @xmath317 ( lower panel ) .
@xmath40 , @xmath318 . in the upper panel
, the blue line corresponds to the `` theoretical '' value , @xmath319 , whereas in lower panel , the plots are fitted by a curve , @xmath320 , with the choice of parameters , @xmath321 , @xmath322,title="fig:",width=8 ]     ( its magnitude , @xmath239 ) plotted as a function @xmath184 for @xmath323 ( @xmath40 ) . at @xmath324
there are two possible solutions for @xmath158 : @xmath325 and @xmath326 , the larger value of which determines large-@xmath152 behavior of @xmath327 .
merger with bulk occurs when @xmath328 , i.e. , at @xmath329 close to the zone boundary ( @xmath330 ) , reentrance of edge solution occurs ( see fig .
[ reent_k ] for details).,width=8 ]     for @xmath331 ( @xmath40 ) . on the reference line at @xmath332 , @xmath333 ( @xmath334 , @xmath335 = \pm 0.133499 ... $]).,width=8 ]
we have extensively studied such characteristic behaviors of the edge wave function in numerical experiments . in fig .
[ rho_j ] , results of such analyses are shown . for the choice of parameters such that @xmath336 , @xmath40 and two different values of @xmath107 : @xmath337 and @xmath338 for comparison . at this value of @xmath336
, we first verified that the wave function @xmath195 takes indeed the form of eq .
( [ evec_edge ] ) , with @xmath286 approximately given by , @xmath339 as in the straight edge case , for @xmath184 corresponding to an edge mode , i.e. , for such a state that are localized in the vicinity of either of the two boundaries , this value of @xmath286 is common practically to all @xmath152 in the strip , as far as a finite amplitude exists . for bulk states which are well extended into the interior of the sample , the wave function takes no longer the form of eq .
( [ evec_edge ] ) , but keeps still a characteristic form as eq .
( [ evec_gen ] ) .    in the two panels of fig .
[ rho_j ] , @xmath340 is plotted as a function of @xmath152 for two different values of @xmath107 . at @xmath337 ,
@xmath340 saturates at rows away enough from the boundary at @xmath201 ( but not too close to the other edge , either ) , to a value close to @xmath319 ( blue line ) , a value which is later `` derived '' ( see fig . [ rho_th_d025 ] ) .
let us assume , ) to a future publication .
instead , we take it here as a plausible hypothesis fully justified by numerical experiments . ] as in the straight edge case , that the wave function in the zigzag edge geometry takes the following form : @xmath341 u_+ \nonumber \\
+ \left [ c_{-1 } \rho_1^{-j } + c_{-2 } \rho_2^{-j } \right ] u_- , \label{sol_zig_gen}\ ] ] where the eigenspinors @xmath283 are always given by eq .
( [ u+- ] ) , the latter found analytically in the limit @xmath342 . we have extensively verified the validity of this hypothesis in numerical experiments . the coefficients @xmath197 are susceptible of system s geometry . here , in a strip geometry
, they satisfy , @xmath343 for one edge mode , and @xmath344 for the other . under this hypothesis , such a behavior as seen in the upper panel of fig .
[ rho_j ] ( @xmath337 ) is interpreted as a consequence of two `` real solutions '' for @xmath158 , which are also _ both _ smaller than unity ( cf .
[ rho_th_d025 ] ) .    on the other hand , at @xmath338 ( in the lower panel of fig .
[ rho_j ] ) @xmath340 shows an oscillatory behavior .
this implies , with the same hypothesis as above , i.e. , the wave function , @xmath195 , given as in eq .
( [ sol_zig_gen ] ) , with the choice of coefficients as eqs .
( [ cs_2 ] ) , a pair of complex solutions for @xmath158 .
indeed , the plots at @xmath338 are nicely fitted by a curve , of the form , @xmath345,\ ] ] with the choice of parameters , @xmath334 and @xmath346 , which will be also ( a posteriori ) justified ( see fig . [ rho_th_d030 ] ) .
( @xmath347 ) calculated in a strip geometry ( here , @xmath318 ) is plotted against the theoretical curve for @xmath348 ( shown in black ) given in eq .
( [ a+- ] ) up to @xmath349 for the lowest - energy eigenmode at different values of @xmath350 ( @xmath40 fixed ) , corresponding , respectively , to the colors : red , green , blue , cyan and magenta .
reference lines indicate the regime of @xmath184 in which the edge modes disappear at the given value of @xmath107.,width=8 ]    comparing these two contrasting cases , notice that the coefficients @xmath286 , estimated to be such as eq .
( [ a_pm_k005 ] ) , are common to the two cases , i.e. , independent of @xmath107 .
one can indeed verify ( by changing the parameter @xmath107 in numerical experiments ) that the eigenspinors @xmath283 remains always the same , as far as the state describes an edge mode ( see fig .
[ ak ] ) ; only @xmath340 changes as a function of @xmath107 .    in fig .
[ ak ] , @xmath312 ( @xmath351 ) is plotted for the lowest - energy eigenmode @xmath161 ( in the upper band ) at different values of @xmath107 .
one can see that in the range of @xmath184 at which @xmath161 is expected to represent an edge mode the plotted points fall roughly on the theoretical curve for @xmath348  cf .
( [ a+- ] )  apart from a small disagreement close to the zone boundary ( @xmath258 ) .
this is indeed a key discovery allowing us to proceed to the next step , of extrapolating the earlier exact solution at @xmath86 to an arbitrary value of @xmath107 .
let us reformulate the recipe for constructing the _ exact _ edge wave function and simultaneously its spectrum in the zigzag edge geometry , which has already been briefly outlined in the introduction ( recall also fig .
[ zig_concept ] ) .    1 .
we have seen in the previous subsection that the edge wave function @xmath254 in the zigzag edge geometry always takes , as far as it describes a _ localized edge mode _ , the form of eq .
[ evec_edge ] with a parameter @xmath286 depending only on @xmath184 ( and @xmath311 , @xmath3 ) .
all our numerical data agree with the hypothesis that @xmath286 , consequently the reduced two - component eigenvector @xmath283 , is independent of @xmath2 .
2 .   on the other hand ,
we know that the problem can be solved exactly at @xmath352 .
we have seen , in particular , that the wave function @xmath254 can be constructed from the same set of spinors @xmath283 with a choice of parameters @xmath286 given _ analytically _ as a function of @xmath184 in eq .
( [ a+- ] ) .    taking also into account the fact that the edge modes , gapless at @xmath267 and characterizing the topological insulator , evolves continuously to the completely flat edge mode at @xmath352 , one can deduce , from these two observations , that the solution of the eigenvalue equation for @xmath187 , i.e. , eq .
( [ reduced ] ) for an _ arbitrary _
@xmath2 should be given , indeed , by @xmath283 , defined as in eq .
( [ u+- ] ) , with the parameter @xmath353 obtained analytically in the limit : @xmath33 ( recall fig .
thus , for a general value of @xmath2 , only @xmath158 and @xmath148 are unknown ( recall that the edge spectrum is no longer flat for a general @xmath2 ) .
but , clearly , they are solutions of @xmath260 \left [ \begin{array}{c } \alpha ( 1- i)\\ 1 \end{array } \right ] = e   \left [ \begin{array}{c } \alpha ( 1- i)\\ 1 \end{array } \right ] , \label{mao - chan}\ ] ] where the @xmath127 matrix @xmath21 is given explicitly as eq . ( [ gamma_zig_mat ] ) . ) and ( [ gamma_zig_mat ] ) , notice that @xmath354 factors out .
so all the coefficients become real . ]
to find _ exactly _ the value of @xmath158 and @xmath148 , one has only to solve this set of equations , and at the end of the calculation , substitute the analytic expression for @xmath77 , i.e. , @xmath286 given in eq . ( [ a+- ] )
obtained in the limit of @xmath33 .
clearly , eqs .
( [ mao - chan ] ) are a set of coupled equations , linear in @xmath148 and quadratic in @xmath158 .
we expect , therefore , two sets of solutions for @xmath355 , which are given as a function of @xmath77 . to each of these two sets of solutions ,
one substitutes either @xmath356 or @xmath357 .
there exist , therefore , four sets of solutions , in general .    unfortunately , the analytic formula for these four sets of solutions are too lengthy to be shown here . instead , we plotted these four solutions in fig .
[ 4sols ] , for different values of @xmath8 .     plotted as a function of @xmath347 at @xmath40 and @xmath337 ( as in fig .
[ rho_j ] , upper panel ) , but for @xmath358 ( upper panel ) and @xmath359 ( lower panel ) . in the upper panel
, the blue line corresponds to the theoretical value , @xmath360 whereas , in the lower panel the plots ( in red ) are fitted by the curve , @xmath320 , with the choice of parameters , @xmath361 , @xmath362.,title="fig:",width=8 ]   plotted as a function of @xmath347 at @xmath40 and @xmath337 ( as in fig .
[ rho_j ] , upper panel ) , but for @xmath358 ( upper panel ) and @xmath359 ( lower panel ) . in the upper panel
, the blue line corresponds to the theoretical value , @xmath360 whereas , in the lower panel the plots ( in red ) are fitted by the curve , @xmath320 , with the choice of parameters , @xmath361 , @xmath362.,title="fig:",width=8 ]     plotted as a function @xmath184 at @xmath363 ( again , as in fig .
[ rho_j ] , upper panel ) . enlarged picture for @xmath184 close to the zone boundary : @xmath364 the reentrance of edge solution occurs at @xmath365 , whereas two real solutions for @xmath158 is possible when @xmath366 two reference lines are at @xmath358 and at @xmath367 , on which @xmath340 was plotted in fig .
[ reent_j ] .
the value of @xmath158 on these lines are , @xmath368 and @xmath369 on @xmath370 ( case of real solutions ) , whereas @xmath371 on @xmath372.,width=8 ]    ) .
this is an enlarged image of the spectrum in the zigzag edge geometry as shown in fig .
[ spec_zig ] . here
, the parameters are chosen such that @xmath373 , @xmath40 , @xmath212 .
two reference lines at @xmath374 and at @xmath375 introduce three different momentum regions : i ) @xmath376 $ ] , ii ) @xmath377 $ ] and iii ) @xmath378 $ ] , corresponding , respectively , to i ) the ordinary edge , ii ) the bulk and iii ) the reentrant regimes.,width=8 ]      _ reentrance _ of the edge mode is another characteristic feature of the edge mode of zigzag geometry , and occurs close to the zone boundary , @xmath379 , when @xmath107 is not too large : @xmath380 . very remarkably , the spectrum looks completely `` innocent '' when this occurs , i.e. , the edge mode , say , the lowest energy ( @xmath381 ) mode in the upper band looks almost completely degenerate with the bottom of the ( bulk ) spectrum ( @xmath382 ) , in this regime of @xmath184 ( see fig . [ spec_zig_enl ] ) .
existence of an edge mode of such specific character is , on the other hand , nothing exceptional in the zigzag edge geometry . at a value of @xmath2 ,
e.g. , @xmath383 or @xmath384 as in fig .
[ rho_j ] , such reentrant edge modes are indeed existent .
if one focuses on the wave function of , say , the lowest energy mode in the upper band , after touching the lower band at @xmath267 , it continues to be spatially localized when @xmath184 is small enough , but as the spectrum merges with the bulk continuum , the wave function also starts to penetrate into the bulk .
however , close to the zone boundary , it starts to be localized again .
this is what we call the reentrance of edge modes
.    figs .
[ reent_j ] and [ reent_k ] highlight the behavior of such reentrant edge modes , naturally in a different regime of @xmath184 from , say , fig .
[ rho_j ] .
[ reent_j ] shows the behavior of @xmath340 at @xmath385 and @xmath367 for @xmath363 . at this value of @xmath184
, the wave function @xmath195 takes always the form of eq .
( [ evec_edge ] ) , but with a different set of parameters for @xmath286 from the case of fig .
[ rho_j ] , upper panel , since it still depends on @xmath184 .
the two plots for @xmath340 in fig .
[ reent_j ] show two typical behaviors of the edge wave function , i.e. , one corresponding to real and the other to complex solutions for @xmath158 . as we have extensively studied in the case of ordinary edge modes ( appearing at @xmath386 ) , the two contrasting behaviors of @xmath340 as a function @xmath152 ( in fig .
[ rho_j ] ) are naturally understood by referring to the theoretical curve of @xmath158 as a function of @xmath184 , e.g. , such as the one shown in fig .
[ rho_th_d025 ] .    what is rather remarkable here , in the case of fig .
[ reent_j ] , is that this crossover between real and complex solutions occurs within a tiny change of @xmath184 , i.e. , from @xmath385 in the upper panel to @xmath367 in the lower panel .
this drastic change is , however , quite reasonable from the viewpoint of fig .
[ reent_k ] .
the upper panel of fig .
[ reent_j ] shows a monotonic decay , which converges asymptotically to a single exponential decay .
this is consistent with the behavior of theoretical curve for @xmath158 as a function of @xmath184 in fig .
[ reent_k ] .
the latter implies two real solutions for @xmath158 at @xmath385 : @xmath368 and @xmath369 the latter coincides with the value of @xmath340 in fig .
[ reent_j ] at which it saturates . at @xmath367 , on the other hand , the plots for @xmath340 are nicely fitted by the curve , @xmath320 , with the choice of parameters , @xmath387 and @xmath388 ( see the lower panel of fig .
[ reent_j ] ) .
this is a clear fingerprint that the reentrant edge mode at this value of @xmath184 corresponds to a pair of complex solutions for @xmath158 .
is plotted as a function of @xmath184 at @xmath389 , @xmath40 .
different curves correspond to different size ( width ) of the system : @xmath318 ( blue ) , 200 ( green ) and 300 ( red ) . the two reference lines are placed at @xmath390 and @xmath391 ( @xmath392 $ ] corresponds to the bulk regime ) .
the plots reveal an extremely small but still a finite binding energy of the reentrant edge modes.,width=8 ]    does the reentrant edge mode really have zero binding energy ? in order to address this question , we ( re)plotted the energy spectrum ( @xmath393 , to be precise ) but in an enlarged scale roughly by one thousand times in fig . [ en_bind ] .
first , for an `` ordinary '' edge state , occurring at @xmath394 the value @xmath393 is much above the threshold at this scale .
@xmath393 takes a value of order @xmath395 for such ordinary edge state . as for the reentrant edge mode , fig .
[ en_bind ] reveals that it has indeed an extremely small but still a finite binding energy .
notice different behaviors of @xmath393 as a function of @xmath32 in the bulk and reentrant regions of @xmath184 .
the former ( the latter ) corresponds to @xmath396 $ ] ( @xmath378 $ ] ) . in the bulk region
@xmath393 is expected to vanish in the thermodynamic limit .
[ en_bind ] shows indeed that the binding energy of reentrant edge mode , @xmath397 , is thousand times smaller than that of the ordinary edge state .
this implies the appearance of an extremely small energy scale which was not existing in the original hamiltonian ( cf .
kondo effect ) .
the reentrance of edge mode is indeed a unique feature , in its contrasting properties in real and momentum space , and in the appearance of an extremely small energy scale .
we have highlighted in this paper various unique properties of helical edge modes in @xmath0 topological insulator .
we have extensively investigated a lattice version of the bhz model , under different edge geometries .
one of the specific characters of bhz model is that the spin hall conductance in the bulk changes its sign in the middle of topological phase ( at @xmath23 ) , i.e. , @xmath398 , respectively , for @xmath213 and for @xmath399 , though both represent a non - trivial value . from the viewpoint of bulk - edge correspondence
, this information should be also encoded in the edge theory .
we have seen that the change of @xmath98 manifests in a very different way in the @xmath7- ( straight ) and @xmath11- ( zigzag ) edge geometries . in the @xmath7-edge case ,
the edge spectrum changes its global structure in the two parameter regimes , i.e. , the main location of the mode moves from the zone center for @xmath9 , to the zone boundary for @xmath10 . as a result
, the group velocity at the intersection with fermi level reverses its sign , leading to change of the sign in landauer conductance at the edge . in the @xmath11-edge case , on the other hand , the edge spectrum is symmetric w.r.t .
@xmath23 , i.e. , neither change of the position of gap closing , nor the reversal of group velocity at @xmath23 .
the change of @xmath98 is here encoded in the swapping of left- and right- going edge modes of the same spin .
much of our focuses has been on the analysis of the zigzag or @xmath11-edge geometry , the latter showing , as a consequence of specific way in which the bulk topological structure is projected onto the 1d edge , a number of unique features , such as the completely flat edge spectrum at @xmath23 , and the reentrance of edge modes .
we have also shown , here in a half - empirical way , that the exact edge wave function for zigzag edge geometry can be constructed , by extrapolating the solution at @xmath23 .
the reentrant edge mode , though sharing much of its characteristics with the usual edge mode in real space , introduces a new extremely small energy scale which was absent in the original bhz model .
ki , ay and ah have been much benefited from useful discussions with jun goryo on the bulk / edge correspondence .
ki also acknowledges christoph brne , hartmut buhmann and laurence molenkamp for their detailed explanation of the experimental situations in stcm ( kyoto ) , ngss-14 ( sendai ) and qhsyst10 ( dresden ) conferences .
ki and ay are supported by kakenhi ( ki : grant - in - aid for young scientists b-19740189 , ay : no .
08j56061 of mext , japan ) .
let us first recall that the eigenvector @xmath400 , which has appeared in eq .
( [ dk_pm ] ) , is a standard su(2 ) spinor , here chosen to be single - valued .
an eigenvector , corresponding to a positive energy eigenvalue @xmath401 , is @xmath402 with @xmath53 satisfying @xmath49 , and represented as , @xmath403 \\ \nonumber & = & { 1\over\sqrt{2d(d - d_z ) } } \left [ \begin{array}{c } d_x - i d_y \\ d - d_z \end{array } \right ] .
\label{evec}\end{aligned}\ ] ] @xmath404 , @xmath405 are polar coordinates in @xmath51-space , satisfying the relations such as , @xmath406    in order to find an edge solution , we focus on a solution of the form , @xcite @xmath407 say , in a semi - infinite plane : @xmath408 . the spatial dependence in the @xmath136-direction can be taken into account by applying pierls substitution : @xmath409 to @xmath75 s in @xmath48 .
the eigenenergy of such a solution is obtained by a simple replacement : @xmath410 in eq .
( [ wine ] ) , i.e. , @xmath411 this can be regarded as a quadratic equation w.r.t .
its two solutions are , @xmath413 for a given set of @xmath74 and @xmath148 , there are two possible values for @xmath412 , or equivalently , four possible values for @xmath414 . of course , in a semi - infinite plane ,
say , @xmath408 the edge solution of the form , eq . ( [ evec_kappa ] )
should decay as @xmath415 , so only two of such solutions are relevant .
we expect that the edge spectrum behaves as @xmath416 in the limit of@xmath147 .
let us parametrize the two solutions in this limit as @xmath417 ^ 2\equiv p\pm q , \label{kappa0}\ ] ] where @xmath418 when @xmath419 , @xmath420 and @xmath421 is real .
since @xmath422 as far as @xmath423 is real , eq . ( [ sol_kappa ] ) represents two positive solutions for @xmath412 , i.e. , the wave function represented by eq .
( [ evec_kappa ] ) shows simple exponential damping . on contrary
, we expect a pure imaginary solution for @xmath414 for an extended state in the bulk . on the other hand ,
when @xmath424 , @xmath425 is always positive but @xmath421 becomes purely imaginary .
thus two solutions for @xmath414 become complex numbers conjugate to each other . in this case , the wave function represented by eq .
( [ evec_kappa ] ) shows damped oscillation .
@xmath427 . for a given value of @xmath74 and @xmath401 ,
we thus have identified two solutions characterized by different values of @xmath414 . in order to construct a general solution in the presence of a boundary , we need to take a linear combination of these two solutions , i.e. , @xmath428 e^{\kappa_+ y }   +   c_-   \left [   \begin{array}{c }   u_-\\   v_-   \end{array }   \right ] e^{\kappa_- y},\end{aligned}\ ] ] where @xmath283 , @xmath429 are short - hand notations for @xmath430 and @xmath431 .
we now fix the boundary condition at @xmath432 , which we choose to be , @xmath433 \left [ \begin{array}{c }
c_+\\ c_- \end{array } \right]= \left [ \begin{array}{c } 0\\ 0 \end{array } \right ] , \label{bc}\ ] ] which implies the following secular equation : @xmath434=0\ ] ] this leads to , @xmath435      let us check whether this solution contains the edge modes .
we expect that the edge spectrum bahaves , as @xmath147 , @xmath416 . recall that in this limit , eq . ( [ sol_kappa ] ) reduces to eqs .
( [ kappa0 ] ) and ( [ ab ] ) .
( [ req_bc ] ) is also simplified in this limit , as @xmath436 focusing on the case , @xmath437 and @xmath438 , and using the parameterization in eqs .
( [ kappa0 ] ) and ( [ ab ] ) , one can readily verify , @xmath439 thus eq . ( [ req_bc0 ] ) is safely satisfied .    how about the first order corrections ? i.e. , contributions of order @xmath440 to the energy spectrum , @xmath149 .
first note that there is no @xmath440-correction to @xmath441 .
one can , therefore , safely replace , at this order , @xmath441 s in eq .
( [ req_bc ] ) with their values at @xmath147 , @xmath416 , i.e. , @xmath442 . \label{req_bc1}\ ] ] we have already seen that the first two terms cancel , whereas @xmath443 ^ 2= 2a+2\sqrt{p^2-q^2}={a^2\over b^2}.\ ] ] thus , the edge spectrum in the continuum limit is determined to be , @xmath444 remarkably , the slope of the edge spectrum depends only on a single parameter , @xmath311 .
an interesting question is to what extent this conclusion is general ?
if one calculates the edge spectrum , using a tight - binding model , generally the results depend on the way edges of the sample are introduced with respect to the lattice . in the case of zigzag edge ,
in particular , apparently the edge spectrum does not converge to eq .
( [ edge_cont ] ) even in the long - wave - length limit : @xmath147 .