document: graphene has garnered interest from broad spectrum of communities , ranging from those aiming at atomic scale circuit devices to those searching for new topological phases . both communities sought after ways to gap the massless dirac spectrum . 
   the realization of a gate - induced band - gap in the bernal stacked bi - layer graphene  @xcite following the prediction in ref . 
   @xcite brought the holy grail of graphene based transistor one step closer to reality . 
   however , the sub - gap conductance measured by  @xcite with weak temperature dependence well below the optically measured gap as large as 250 mev@xcite introduced a new puzzle and obstacle : the gapped bilayer is not as insulating as it should be . 
   dominant transport along physical edge of the samples proposed earlier by @xcite have been ruled out by corbino geometry measurements@xcite , which observed two - dimensional variable range hopping type temperature dependence , independent of geometry . in this paper 
   we predict existence of topological gapless channel of transport along recently imaged ab - ba tilt boundary network @xcite which solves the puzzle . 
   the predicted topological edge state holds the promise of the first realization of topological surface(edge ) state hosted by structural topological defect . 
   though there has been much theoretical interest in topological gapless modes hosted by structural topological defects@xcite no such topological gapless mode has been observed so far . 
   the lattice dislocations in three dimensional crystals previously discussed occur deep in the sample that is not directly accessible . 
   however , the tilt boundary of interest have recently been observed@xcite . 
   the tilt boundary is a structural topological line defect along which each neighboring layer is displaced by one inter - atomic spacing 
   . such defect can occur due to the third dimension added by the stacking of the graphene layers ; it forms a boundary between two inequivalent stacking structures frequently referred to as ab and ba . 
   here we show that the tilt boundaries host gapless modes of topological origin and form the first example of a naked structural defect hosting topological electronic states . 
   ( @xmath2 ) sublattice sites . 
   @xmath3 represent hopping matrices for a tight - binding model . 
   , scaledwidth=50.0% ]    topological aspects of gapped multi - layer graphene have been previously discussed@xcite and it was pointed out that they should exhibit quantum valley hall effect with corresponding edge states . however , to this date there has been no experimental detection of proposed edge state  @xcite . 
   moreover , little is known about how the topological aspects of gapped multi - layer graphene relates to topological insulators  @xcite . 
   the idea of classifying different topological insulator ( superconductor ) candidates based on symmetries  @xcite have played a key role in the field of topological insulators . in particular the observation that additional symmetries such as the crystalline symmetries can enlarge the possibilities of topological phases@xcite led to the discovery of three - dimensional topological crystalline insulators@xcite . 
   on one hand we propose feasible experiments to detect topological edge states at naturally occurring tilt boundaries . at the same time , we make first concrete application of the spt approach  @xcite for two dimensional ( 2d ) system and study a large class of gapped graphene systems placing the quantum valley hall insulator in the larger context and predicting conditions for topological superconductors . 
   the rest of the paper is organized as follows . in section 
   [ sec : micro ] we show that a ab - ba tilt boundary in gated bilayer graphene supports gapless edge states through explicit microscopic calculations . 
   specifically we consider an abrupt boundary in tight - binding model and then investigate the effect of strain using ab - initio calculation . in section [ sec : topo ] we show that these edge - states are protected by no valley mixing , electron number conservation , and time reversal ( @xmath1 ) symmetries within the framework of spt . 
   hence we identify chirally stacked gated @xmath4-layer graphene layers as _ time - reversal symmetric _ 
   @xmath0-type spt . in section [ sec : expt ] 
   we discuss experimental implications . 
   finally in section [ sec : summary ] we summarize the results and comment on practical implications . 
   +    c +    fig . 
   [ fig : bilayer ] and fig . 
   [ fig : domainwall](a ) show tilt boundaries of interest in gapped bernal stacked bi - layer graphene . in the case sketched , strain is concentrated at the tilt boundary with the top layer stretched by one inter - atomic spacing with respect to the bottom layer . for a general orientation , tilt boundaries can involve both strain and shear . as the tilt boundaries in layered graphene form a type of topological line defects in structure , they can be characterized using the tangent vector @xmath5 and the burger s vector @xmath6 . 
   the tangent vector @xmath5 points along the tilt boundary which can point along any direction with respect to the burger s vector @xmath6 . 
   when the tilt boundary only involves strain as in the case depicted in fig . 
   [ fig : bilayer ] and fig . 
   [ fig : domainwall](a ) , the @xmath6 is perpendicular to @xmath5 . in the opposite extreme limit of @xmath7 , 
   shear is concentrated at the boundary . 
   independent of the angle between @xmath6 and @xmath5 , the burger s vector magnitude is the inter - atomic spacing i.e. @xmath8 for a bilayer system , as it is shown explicitly for the strain tilt boundary in fig . 
   [ fig : domainwall](a ) . 
   since @xmath9 is a fraction of the bravis lattice primitive vector magnitude @xmath10 , the bilayer domain boundary is a partial dislocation from quasi two - dimensional view . in a general mult - layer a vertical array of these partial dislocations form a tilt - boundary . in typical samples , the domain wall separating the ab and ba stacked domains has substantial width spanning 5 - 20 nm and the angle between @xmath6 and @xmath5 ranges between @xmath11 and @xmath12  @xcite . 
   [ fig : bilayer ] allow us to makes two important microscopic observations about tilt - boundaries . 
   ( 1 ) as the boundary requires a shift of one - layer with respect to the other by one inter - atomic spacing along the bond direction , there are three natural directions for the tilt - boundary to run for each fixed angle between @xmath6 and @xmath5 . 
   ( 2 ) the boundary is arm - chair for @xmath13 ( pure shear ) whereas it is zigzag for @xmath14 ( pure strain case shown in fig . 
   [ fig : bilayer ] ) . based on observation ( 1 ) 
   we expect a given type of tilt boundary to possibly form a triangular network seen in experiments@xcite . 
   the observation ( 2 ) combined with earlier microscopic studies of boundary condition effects on edge states in ref . 
   @xcite implies that electronic spectrum at tilt - boundaries with @xmath13 will be gapped though the gap magnitude will be small when the tilt boundary is spread over finite width . 
   in this section , we consider the electronic structure of tilt boundaries with @xmath14 and return to more general case in the section [ sec : expt ] . as it is shown schematically in fig . 
   [ fig : bilayer ] , the bulk of each domain is gapped in the presence of inter - layer hopping and the external electric field . the latter is important for breaking the inversion symmetry between the layers and gapping otherwise touching bands  @xcite . 
   below we present two separate microscopic calculations of the ab - ba tilt boundary electronic structure for @xmath14 , which shows gapless edge states . 
   we consider a tight - binding hamiltonian with nearest - neighbor intra- and inter - layer hopping . for the ab - stacking region ( see fig . 
   [ fig : domainwall](a ) ) , @xmath15 where @xmath16 is a layer index , @xmath17 and @xmath18 are intra - layer and interlayer nearest neighbor hopping respectively , and @xmath19 is the chemical potential difference between two layers due to the gate voltage . 
   ( m , n ) labels the position of the two site unit - cell with @xmath20 and @xmath21 annihilating electrons at the two sites of layer @xmath22 . for the ba stacked region , 
   the only change is in the inter - layer term with @xmath23 replacing the inter - layer term in @xmath24 . 
   as we address the effect of strain through ab - initio simulation , we focus here on a sharp tilt boundary as shown in fig . 
   [ fig : domainwall ] . 
   we plot the energy spectrum in fig . 
   [ fig : domainwall](b ) with the model parameters set to be @xmath25ev , @xmath26ev , and @xmath27ev . 
   the size of the system was 200 unit - cell in each direction under periodic boundary condition with each domain spanning 100 unit - cell width separated by two sharp tilt boundaries . from the spectra it is clear that @xmath28 and @xmath29 valleys each have two edge states per spin . 
   further investigation of the wave function shows right and left moving edge states associated with the given valley are spatially separated between the two edges : the edges offer valley filtering ( see fig.[fig : domainwall](c ) ) .      . 
   the fermi level is indicated with the blue dashed line and the bands with marked linear dispersion relation intersecting near the fermi level are outlined by the purple ( light ) curves . 
   ( b ) top and side views of the charge distribution for a state near the @xmath30 point of ( a ) . 
   the yellow rectangles indicate the ab - ba domain wall structure . 
   , scaledwidth=48.0% ]    the electronic structure in realistic tilt boundary will be affected by both the span over which the lattice structure transition from ab to ba stacking , and the strain concentrated at the tilt boundary . to address these issues we further carried out first - principles calculation using density functional theory ( dft ) within the local density approximation ( lda )  @xcite . 
   we constructed a periodic supercell with two tilt boundaries characterized by burgers vectors identical to the ones considered within tight binding calculation , but the domain wall is set to have finite width over which the one unit - cell mismatch is spread . 
   we have tried several configurations with domains and tilt boundaries of different widths to find qualitatively similar results . in the rest of this paper we focus on a representative example with 3.1 nm wide domain wall between 1.5 nm wide domains . 
   the choice of narrow width for the domains was due to the limitation in the simulation capability . 
   to this 2-d system , we then applied a slightly exaggerative perpendicular electric field of 5 v / nm . by relaxing the domain boundaries until the forces reaches below 0.5 ev / nm 
   , we took the effect of both strain and the width into account . 
   we have tried several configurations with domains and tilt boundaries of different width but otherwise similar setting to find qualitatively similar results . 
   the dft - lda simulation results are presented in fig . 
   [ fig : ldaband ] , which confirms the existence of the gapless edge states predicted by the tight - binding model . the electronic structure along the extended direction of edge 
   is shown in fig .  [ 
   fig : ldaband](a ) , focusing on the region near @xmath28 point . in this figure two gapless 1d dirac dispersion 
   is clearly resolved from the gapped bulk specta , with two distinct dirac points @xmath30 and @xmath31 in the vicinity of the @xmath28 point . 
   the strain concentrated at the tilt boundary causes energy splitting of @xmath30 and @xmath31 states ; the energy of 1d dirac points is increased by compressive strain and decreased by tensile strain . 
   this energy splitting will become negligibly small in realistic tilt boundaries with much wider span , as the strain will become smaller . 
   [ fig : ldaband](a ) shows that the existence of gapless edge states found in our tight - binding calculation are robust against long range perturbations such as tilt boundary width or subtle bond - length variation inside the domain wall , as well as the strain at the tilt boundary . 
   we now turn to the spatial distribution of the gapless edge states . for illustration , 
   [ fig : ldaband ] ( b ) shows the wavefunction amplitude of the state slightly above the @xmath30 point . 
   as expected from the tight - binding results , the charge distribution is prominent inside the domain wall but rapidly decays away from the tilt boundary . 
   the dirac points of edge states and the decaying feature of edge states always exist , making us believe it is a universal feature . 
   meanwhile , the charge for the states near @xmath30 is highly localized on the layer subject to compressive strain , regardless of the direction of the applied electric field . 
   the situation for the states in vicinity of @xmath31 is similar except that the charge prefer to highly localized on the layer under tensile strain . hence forming a layer selective contact to an isolated edge state 
   could be a mechanism for valley filtering@xcite .    finally , we comment on the so far ignored effect of interaction . 
   if the edge boundaries have substantial width , forward scattering part of the coulomb interaction will be the dominant correlation effect to the edge states and lead to luttinger liquid behavior @xcite . 
   in order to address the robustness of the edge states , we investigate topological aspects of the low energy effective theory in the continuum limit . we first show that the ab - ba tilt boundary can be mapped to a gate - polarity boundary of uniform bi - layer . based on this mapping and results of refs . 
   @xcite on the gate - polarity boundary , we discuss the valley chern number of the tilt boundary edge states . 
   we than apply the notion of spt  @xcite and identify chirally - stacked multi - layer graphene as a realization of @xmath32-type spt , protected by time - reversal ( @xmath1 ) , absence of valley mixing , charge conservation symmetries . 
   this identification enables us to address effects of symmetry changes : topological quantum phase transitions . 
   there are recent studies of such perturbations for specific cases such as rashba spin - orbit coupling  @xcite and magnetic ordering  @xcite . through our first application of spt classification scheme by @xcite to a concrete physical system of multi - layer graphene 
   , we obtain an exhaustive systematic study of topological quantum phase transition possibilities . 
   the low energy effective hamiltonian near the @xmath28 valley for uniformly ab or ba stacked bi - layer is @xmath33 where @xmath34 or @xmath35 sign should be used for ab or ba stacking respectively . in eq .  , @xmath36 , @xmath37 s are pauli matrices acting on the sub - lattice indices , and @xmath38 s are pauli matrices acting on the layer indices . 
   the effective hamiltonian near @xmath29 is @xmath39 . 
   now it is straight forward to show that ba stacking is equivalent to ab stacking subject to the opposite gate polarity . at zero field 
   the hamiltonian of the ab stacked bilayer can be transformed to that of the ba stacking by interchanging the two layers via the following unitary transformation : @xmath40 , with @xmath41 . for a gated bi - layer however , the gate polarity has to flip since @xmath42 hence at the level of low energy effective theory , the tilt boundary between ab and ba stacking under uniform external field is equivalent to the gate polarity domain wall of structurally uniform bi - layer proposed by @xcite .    the above equivalence combined with earlier results on valley chern number of gated chirally stacked multi - layer offers the topological origin of the helical edge states observed in the microscopic calculation of section  [ sec : micro ] . 
   first for bi - layer and then for general @xmath4-layers , it was shown that low energy effective theory of chiral stacked @xmath4-layer under uniform vertical electric field @xmath19 has finite chern number per spin for each valley of equal magnitude and opposite sign  @xcite : @xmath43 the chern numbers in eq . 
   can be obtained by integrating the berry curvature over momenta @xmath44 continuing the linearized dispersion to infinity . 
   this combined with the equivalence relation of eq . 
   means the valley chern numbers change sign at the tilt boundary . 
   such sign change leads to @xmath45 across the tilt boundary and @xmath4 branches of valley helical edge states  @xcite , as long as the two valleys @xmath28 and @xmath29 remain distinct . 
   hence the two valley helical edge states per spin observed in section [ sec : micro ] originate from the valley chern number change across the tilt boundary , as in the gate polarity boundary edge states  @xcite . hence , our prediction is the tilt boundaries will be the first experimentally observed crystalline topological defects to host topological gapless mode due to change in the chern number . 
   we now apply the procedure for identifying the class of spt based on symmetries of free fermion hamiltonian developed by @xcite , which predicts possible number of protected edge ( surface ) states . 
   this procedure allows us to consider additional symmetries in the multi - layer graphene in addition to the @xmath46 , @xmath1 , and @xmath47 taken into account in the pioneering work by @xcite , and by @xcite . 
   the procedure consists of three steps : ( 1 ) find a gapless dirac hamiltonian ( by keeping the kinetic term only ) with the same symmetries . 
   then we find all the symmetry preserving mass terms that can gap out the gapless part and are amenable to classification using clifford algebra . 
   this is based on the assumption that the spt order is robust as long as the energy gap stays finite and the symmetries remain the same and hence any gapped hamiltonian can be adiabatically transformed into a gapped dirac hamiltonian . 
   ( 2 ) express the hamiltonian and the conserved quantities associated with symmetries in the majorana basis . 
   this leads to the clifford algebra ( i.e. real representation of the dirac algebra ) associated with the gapless part of the hamiltonian . 
   ( 3 ) find the space of mass matrices that anti - commutes with all the generators of this clifford algebra . 
   the resulting space may have disconnected pieces , the number of which gives the classification of the spt . 
   two mass matrices are topologically distinct if and only if they belong to two different pieces . applying this procedure to chiral multi - layer graphene will enable us to study phase transitions into different spts upon symmetry changes .    for chiral multi - layer graphene system , we assume no valley mixing , electron number conservation @xmath48 , and time reversal ( @xmath1 ) symmetries . 
   the relevant gapless dirac hamiltonian is : @xmath49 where @xmath50 , and @xmath51 , and @xmath52 , in which @xmath53 is given by the number of layers . 
   this hamiltonian can be written in the majorana fermion basis using the following decomposition @xmath54 where @xmath55 denotes the majorana fermion satisfying @xmath56 where @xmath57 ( @xmath58 ) denotes the a , or b sublattice ( @xmath28 or @xmath29 valley ) indices , and @xmath59 denotes the flavor of the majorana fermions ( + or - ) . in the majorana fermion basis the hamiltonian is represented as follows :    @xmath60    where @xmath61 is a real anti - symmetric matrix ( differential operator ) , and @xmath62 is an eight component vector whose components are @xmath63 . 
   now we express the conserved quantities associated with the symmetries of the hamiltonian in the majorana fermion basis . 
   first , no valley mixing combined with total electron number conservation symmetry leads to separate conservation of the electron number at each valley @xmath64 and @xmath65 . 
   hence the total electron number @xmath66 , and the valley polarization @xmath67 are conserved . in the majorana fermion basis , 
   @xmath68 where @xmath69 ( @xmath70 ) pauli matrices act on the majorana flavors ( valley indices ) . 
   so defined @xmath71 and @xmath72 satisfy @xmath73 .    under time 
   reversal symmetry @xmath1 , @xmath28 and @xmath29 valley indices are exchanged i.e. , @xmath74 . hence @xmath1 acts like @xmath75 in the valley basis with the matrix part of the time reversal operator satisfying @xmath76 . 
   on the other hand , no valley mixing implies the hamiltonian is invariant under @xmath77 transformation which acts like @xmath78 in the valley basis . 
   hence , in the presence of no valley mixing symmetry , we can define a new time reversal operator @xmath79 , which acts like @xmath80 . in terms of majorana fermions 
   @xmath81    in order to find the relevant clifford algebra , we need to form anti - commuting generators in terms of @xmath61 in eq . combined with symmetries , @xmath82 , @xmath83 , and @xmath84 
   however , the symmetries require @xmath85=\left[\mathcal{a},\hat{q}_{v}\right]=0 $ ] and @xmath86=0 . moreover , symmetry operators satisfy @xmath87 , @xmath88 , and @xmath89 relations . using these relations 
   , we find the full set of generators of the relevant clifford algebra as @xmath90 , @xmath91 , and @xmath92 as @xmath93 for @xmath94 . 
   the resulting full set of anti - commutation relations is @xmath95 and it defines a clifford algebra cliff(3,2 ) . 
   now , we will find the space of mass matrices , @xmath96 , that can gap out dirac hamiltonian associated with this clifford algebra in order to obtain spt classification . 
   the mass term with matrix representation @xmath97 should satisfy the following algebra : @xmath98 where we have normalized @xmath99 . solving the above equation yields the allowed space for the mass matrix , @xmath100 . 
   it has been shown @xcite that @xmath96 which solves eq . 
   for the case of cliff@xmath101 is @xmath102 the spt classification is then given by the number of disconnected pieces in the space of mass matrix @xmath100 i.e. its zeroth homotopy group : @xmath103 . 
   using eq . , it can be verified that @xmath104  @xcite . 
   consequently , each class of the time reversal invariant multilayer graphene in the absence of intervalley scattering is indexed by a @xmath32-valued number : in this case the valley chern number . 
   [ tab1 ]    [ cols="^,^,^,^,^",options="header " , ]     now we are in a good position to consider symmetry changes . important to note here that spontaneously ordered phases 
   can be considered alongside systems under external field , as once a system is deep inside the ordered phase it can be treated within mean - field theory . 
   we first consider the symmetry reduction possibilities while maintaining spin degeneracy ( see table [ table : spinless ] ) . 
   if we only break the time reversal symmetry , the system is characterized by two independent topological indices @xmath105 , hence the classification is given by @xmath106 . 
   we refer to these spt phases as intra - valley quantum ( anomalous ) hall ( qah ) states . 
   such phases may be realized by placing trigonally - strained graphene @xcite under an external magnetic field , as the sign of the pseudo - magnetic field caused by strain is opposite for the two valleys . 
   further reducing the symmetry by introducing inter - valley scattering leads to inter - valley qah state indexed by a single integer @xmath32 : the total chern number  . breaking 
   electron number conservation turns the above insulators into superconductors . following the procedure above 
   , we obtain the same classification for the topological superconductors in 2d , resulting in topological valley superconductor ( tvsc ) and intra- or inter - valley topological superconductors . 
   table [ table : spinless ] summarizes all symmetry reduction possibilities and their classifications starting form gated multi - layer graphene . 
   now we consider extending our classification to take the electron spin into account as a dynamical degree of freedom . 
   this will allow us to consider interaction effects at the level of spin ordering . with both spin and valley degrees of freedom , 
   there are three symmetry operators related to time reversal : @xmath107 exchanges two valleys , @xmath108 which acts on the spin indices flips spin , and @xmath109 does both . 
   if both @xmath110 and @xmath111 are broken , the classification reduces to that of spinless electrons ( see upper part of table [ table : spinful ] ) . 
   however , taking any of these symmetry operators into account leads to new classes ( see lower part of table [ table : spinful ] ) . 
   when all @xmath112 , @xmath110 ( and as a result @xmath111 ) and no valley - mixing symmetries are imposed the state corresponds to the so - called `` layer antiferromagnetic '' ( laf ) phase predicted in refs . 
   @xcite and possibly occurring as a ground state in bi - layer graphene at neutrality point  @xcite . 
   in this phase the product of spin and valley of edge quasi - particles is locked to their momentum . 
   therefore , one may index this state by its _ spin - valley chern number _ 
   @xmath113  @xcite . 
   in this phase , for each valley , the edge states associated with two spins are counter - propagating . as long as 
   these two counter - propagating modes do not couple , there can be any number of them per each valley , leading to @xmath114 spin - valley hall conductivity  @xcite . 
   however there is a form of symmetry allowed coupling between a pair of sets of counter propagating edges ; this coupling can gap out the edge modes  @xcite . 
   therefore , the number of symmetry protected edge modes for each valley is @xmath113 mod @xmath115 . 
   hence , we obtain a @xmath116 classification for this phase labeled by @xmath117 as supposed to @xmath32 classification which would be implied by refs . 
   @xcite . 
   another interesting possibility is breaking @xmath112 , and @xmath110 , while respecting their product @xmath111 and no - valley mixing symmetry . 
   this leads to the quantum spin hall ( qsh ) phase , in which @xmath118 , and @xmath119 due to @xmath111 symmetry , while there is no constraint on the @xmath120 . 
   hence , unlike usual 2d qsh states  @xcite , we obtain @xmath32 . 
   this is because multiple type of time reversal operators can be defined  @xcite in the presence of no valley mixing symmetry . 
   this and other spt possibilities are summarized in table [ table : spinful ] .    with the classification at hand , we return to what it implies for the fate of the gapless edge states at the tilt boundary of chirally stacked and gated n - layer graphene . 
   in general , a topological phase transition between two phases within the same class requires the bulk gap to close and reopen ( as in inter - plateaux transitions in integer quantum hall effect ) . 
   gapless edge states are guranteed at a physical boundary between such two phases . the edge states at the ab - ba tilt boundary we established in the section ii are examples of such edge states . 
   hence the edge states will remain gapless as long as the time ( valley ) reversal and charge conservation symmetry are maintained . 
   on the other hand , symmetry change can either yield a trivial phase which does not support an edge state or a different type of spt with different type of edge states . 
   according to tables i and ii , ruining the no valley mixing symmetry is the only way to render the system trivial and gap the edge states . 
   however this requires large momentum transfer which generally requires fine tuning unless the unit - cell becomes enlarged either for the entire system  @xcite or for the edge through arm - chair edge , or short - range disorder such as a vacancy breaks @xmath121-@xmath122 sub - lattice symmetry  @xcite . 
   as both rarely occur , we anticipate gapless edge states at most tilt boundaries . 
   in particular , the natural zigzag boundary formation for the tilt boundaries make such edge states more robust than the edge states in gate polarity boundaries  @xcite .    among symmetry change possibilities leading to another spt , transition from qvh with spin degeneracy to laf phase where spin is 
   a dynamic degree of freedom is of particular interest as laf is suspected to be the ground state of bi - layer graphene near neutrality point  @xcite . upon this phase transition 
   the nature of edge states change from spin degenerate valley helical states ( qvh ) to spin - valley hall edge states ( laf ) . 
   in such transitions , the bulk gap has to close and reopen ; this is indeed seen in the experiment of  @xcite . 
   in this section , we discuss the experimental implications of our findings . specifically we propose the transport through the network of tilt boundaries as a solution to the long standing mystery of sub - gap transport@xcite . 
   further we propose feasible experiments to test the proposal .    in order to discuss the topological transport through the network of tilt boundaries observed in refs . 
   @xcite we should first discuss the effect of the arbitrary angle between the burger s vector @xmath6 and the tangent vector @xmath5 . 
   microscopic study of tilt boundaries at various angles between @xmath6 and @xmath5 will be presented in the future@xcite . 
   however , it has been known that a single domain bernal stacked gapped bilayer ribbon should support gapless edge state for zig - zag edges , but not for arm - chair edges  @xcite . 
   this is because the arm - chair edge enlarges the unit - cell along the direction parallel to the boundary and makes the projection of k and k@xmath123 valley identical . 
   however , @xcite showed that the polarity boundary edge states only develop barely visible gap which is orders of magnitude less then the bulk gap even for a sharp boundary , and the gap decreases quickly when the polarity boundary becomes smooth . 
   these arguments apply to our tilt boundary and the edge state will develop a small gap when @xmath13 . 
   however , given the large width of the observed tilt boundaries we expect that all straight tilt boundaries will have nearly gapless edge states except those with small angles between @xmath6 and @xmath5 . 
   when the tilt boundary meanders and changes directions , likely there will be portions with small gap segmenting the gapless regions and the transport will occur through hopping between the gapless regions . the observed 2d network of such tilt boundaries would yield 2d variable range hopping temperature dependence @xmath124@xcite at low temperatures governed by the 2d connectivity and the small characteristic gaps of gapped regions . 
   this explains the observed temperature @xmath1 dependence of resistance at low temperatures  @xcite . 
   we propose following experiments to test our proposal . 
   ( 1 ) four terminal transport measurements with two of the contacts , say contacts 1 and 3 , at two ends of a tilt boundary . 
   this would yield highly anisotropic transport proving dominant transport along the tilt boundary i.e. , @xmath125 . 
   ( 2 ) scanning tunneling spectroscopy measurements of local density of states . 
   this should measure a gapless spectrum at the tilt boundary but exhibit a gapped spectrum with the gap magnitude of the optical gap away from the tilt boundary . 
   ( 3 ) thermoelectric imaging . 
   the mid - gap density of state at tilt boundaries would appear in scanning thermopower images . 
   unpublished thermopower imaging data by @xcite indeed show a network with local density of state near fermi energy , that is reminiscent of the tilt boundary network . 
   ( 4 ) edge current imaging using scanning squid which can detect magnetic field generated by edge currents . 
   we showed that spin - degenerate tilt boundaries of gated multi - layer graphene support topological gapless edge states protected by three symmetries : time ( valley ) reversal , no - valley - mixing , and electron number conservation . 
   we demonstrated the existence of gapless edge states through a tight - binding model calculation and a first principal calculation , where the latter took strain effects into account . 
   we then addressed the symmetry protection of the edge states and consequences of symmetry changes within the framework of spt  @xcite . 
   the framework of spt allowed us to place the 2d topological phase supporting the edge states , namely qvh , among various topological insulator / superconductor phases alongside previously postulated qah , laf and qsh . 
   while previous literature postulated qvh , qah , laf and qsh to be all supporting number of edge states growing with the number of layers @xmath4 ( i.e. @xmath0-type in the language of classification ) , we found that the symmetry of laf only protects odd number of edge modes for each valley . 
   hence laf is a @xmath116-topological insulator much like quantum spin hall insulator@xcite . 
   transition between these different spt s require closing and re - opening of the bulks gap as already been observed in ref . 
   @xcite . 
   we predict the naturally occurring tilt boundary  @xcite to be the first topological structural defect hosting topologically protected gapless mode of transport , most importantly , our findings on tilt boundaries combined with the recent observations @xcite solve the long standing mystery of sub - gap transport  @xcite . 
   our explanation can be tested through proposed transport , scanning tunneling spectroscopy and thermopower imaging experiments , and scanning squid experiments . 
   experimental confirmation of the tilt - boundary transport origin of the sub - gap transport will open doors to control the sub - gap transport and enable device application of gated multi - layer graphene systems . 
   : we thank paul mceuen for numerous discussions and sharing his unpublished data . 
   we thank joe stroscio for sharing his unpublished data on thermopower imaging . 
   we thank e.j . 
   mele for useful discussions and sharing the preprint @xcite and j. sethna for useful discussions regarding structural topological defect aspect of the tilt boundary . 
   e .- a.k . and d.n . 
   were supported by nsf award eec-0646547 through cornell center for nanoscience . 
   e .- a.k . and a.v . 
   were supported in part by nsf career grant dmr-0955822 . 
   was also supported in part by nsf grant dmr-1120296 through cornell center for materials research . 
   y.l . and l.y . were supported by nsf grant no . 
   the computational resources have been provided by lonestar of teragrid at the texas advanced computing center ( tacc ) . 
   45ifxundefined [ 1 ] ifx#1 ifnum [ 1 ] # 1firstoftwo secondoftwo ifx [ 1 ] # 1firstoftwo secondoftwo `` `` # 1''''@noop [ 0]secondoftwosanitize@url [ 0 ] 
   + 12$12  & 12#1212_12%12@startlink[1]@endlink[0]@bib@innerbibempty link:\doibase 10.1038/nmat2082 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.74.161403 [ * * ,   ( ) ] http://dx.doi.org/10.1038/nature08105 [ * * ,   ( ) ] http://dx.doi.org/10.1038/nphys1822 [ * * ,   ( ) ] link:\doibase 
   10.1021/nl102459 t [ * * , ( ) ] ,   @noop `` , ''   @noop ( ) ,   @noop * * ,   ( ) http://dx.doi.org/10.1038/nphys1220 [ * * ,   ( ) ] link:\doibase 10.1103/physrevx.2.031013 [ * * , ( ) ] link:\doibase 10.1103/physrevb.82.115120 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.99.236809 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.100.036804 [ * * , ( ) ] link:\doibase 10.1103/physrevlett.95.226801 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.96.106802 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.78.195125 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.85.085103 [ * * , ( ) ] link:\doibase 10.1103/physrevlett.106.106802 [ * * ,   ( ) ] http://dx.doi.org/10.1038/nmat3449 [ * * , ( ) ] link:\doibase    10.1103/physrevb.84.075418 [ * * ,   ( ) ] link:\doibase 10.1103/physrev.136.b864 [ * * , ( ) ] link:\doibase 10.1103/physrev.140.a1133 [ * * , ( ) ] http://dx.doi.org/10.1038/nphys547 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.104.216406 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.107.256801 [ * * ,   ( ) ] link:\doibase    10.1103/physrevb.83.155447 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.82.195438 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.106.156801 [ * * ,   ( ) ] in  link:\doibase 10.1063/1.3149495 [ _ _ ] ,  , vol .  ,  ( ) pp .  ,   link:\doibase 10.1126/science.1191700 [ * * ,   ( ) ] http://dx.doi.org/10.1038/nphys1420 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.82.205106 [ * * , ( ) ] http://dx.doi.org/10.1038/nnano.2011.251 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.95.146802 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.86.205116 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.98.186809 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.98.266402 [ * * , ( ) ] link:\doibase 10.1103/physrevlett.89.266603 [ * * , ( ) ] @noop `` , ''   link:\doibase    10.1103/physrevlett.100.026802 [ * * ,   ( ) ] @noop * * ,   ( ) @noop * * @noop `` , ''   @noop `` , ''
summary: despite much interest in engineering new topological surface(edge ) states using structural defects , such topological surface states have not been observed yet . 
   we show that recently imaged tilt boundaries in gated multi - layer graphene should support topologically protected gapless edge states . 
   we approach the problem from two perspectives : the microscopic perspective of a tight - binding model and an ab - initio calculation on a bilayer , and the symmetry protected topological ( spt ) states perspective for a general multi - layer . 
   hence we establish the tilt boundary edge states as the first concrete example of edge states of symmetry enriched @xmath0-type spt , protected by no valley mixing , electron number conservation , and time reversal @xmath1 symmetries . 
   further we discuss possible phase transitions between distinct spt s upon symmetry changes . 
   combined with recently imaged tilt boundary network , our findings offer a natural explanation for the long standing puzzle of sub - gap conductance in gated bilayer graphene , which can be tested through future transport experiments on tilt boundaries . in particular , the tilt boundaries offer an opportunity for in - situ imaging of topological edge transport