Metal–Semiconductor Behavior along the Line of Stacking Order Change in Gated Multilayer Graphene

We investigated gated multilayer graphene with stacking order changes along the armchair direction. We consider that some layers cracked to release shear strain at the stacking domain wall. The energy cones of graphene overlap along the corresponding direction in the k-space, so the topological gapless states from different valleys also overlap. However, these states strongly interact and split due to atomic-scale defects caused by the broken layers, yielding an effective energy gap. We find that for some gate voltages, the gap states cross and the metallic behavior along the stacking domain wall can be restored. In particular cases, a flat band appears at the Fermi energy. We show that for small variations in the gate voltage, the charge occupying this band oscillates between the outer layers.


I. INTRODUCTION
Stacking order domain walls in multilayer graphene are most usually considered along the zigzag direction.They can be created when some layers are stretched, corrugated, or delaminated along the armchair direction [1][2][3][4].Stacking order change causes the appearance of topological states.When the multilayer is gated and the gap opens, the one-dimensional topological gapless states appear at valleys K and K ′ separated by 4  3 π/a, [1,[5][6][7][8][9].The gapless states are then valley-protected as long as there is no valley mixing caused, for example, by atomicscale defects.It was however shown that even under strong defects like bonds breaking in some layers at the domain wall, some gap states that connect the valence band (VB) and conduction band (CB) continua persist, and the systems preserve metallic character along the stacking domain wall [10][11][12].
The situation is different if we allow stacking domain wall along the armchair direction.Firstly, because along the armchair direction, the valleys appear at the same K = 0 [13].Therefore, the gapless states, which have opposite slopes E(k) at each valley, overlap and may no longer be valley-protected.Secondly, because such stacking change causes strong shear strain at the interface between different stackings, at least for narrow domain walls.
In this work, we investigate the behavior of gapless states in such systems.We consider multilayers in which the interface bonds are broken in some layers to release any shear strain in the system.Fig. 1 shows how the unstrained bilayer (BLG), trilayer (TLG), and four-layer graphene (FLG) with stacking domain walls along the armchair direction can be formed.One can easily imagine that such cracks of layers may be present in realistic multilayers [14,15].The bonds broken in some layers * Electronic address: wj@fizyka.umk.plact as atomic-scale defects.Therefore, the gapless states strongly interact and split yielding the appearance of the effective energy gap and thus semiconducting properties of such multilayers.However, we show that for some gate voltages, some gap states cross and the energy gap disappears.The metallic behavior along the stacking domain wall can be restored.Transitions from semiconductor to metal can be achieved by changing the gate voltage.Moreover, for some gate voltages, the separated gap states in FLG connect at the Fermi level (E F ) lead-ing to a kind of flat band at E F , with charge occupying this band localized mainly in the outer layers.

II. SYSTEM AND METHOD
We investigate Bernal stacked bilayer and rhombohedral four-layer graphene with stacking order domain walls.Fig. 1 demonstrates how the consecutive layers are deposited.Panels (a) and (c) show the investigated BLG and FLG systems, and panel (b) shows an intermediate system of the TLG.Only small fragments of the multilayers are shown, while the investigated systems are infinite in both in-plane directions.The bottom layer is always a uniform graphene, as well as the top layer in the case of FLG.The nodes of the FLG top layer have the same in-plane positions as the nodes of the bottom layer.The remaining layers are broken, they consist of two nonconnected half-infinite graphene layers shifted relative to each other by a C−C along the armchair direction.This allows for the formation of a strainless interface with different stacking orders of rhombohedral multilayers on both sides.Note that due to the presence of the stacking domain wall, the systems studied are periodic only in the direction defined by the domain wall, i.e., in the armchair direction.Fig. 1 (d) shows the stacking order of layers on both sides of the domain wall.For BLG they are simply Bernal stackings AB and BA, respectively, while for FLG we have ABCA and CBAC stackings.
To investigate the electronic structure of the defined multilayers we calculate the local density of states (LDOS) at the interface.The calculations are performed using the Green function matching technique for a system composed of a conductor connected to two semiinfinite leads [16].In our case, the conductor consists of the interface region, as defined in Fig. 1   The voltages we consider are available in the experiment and are of the same order as typically used for multilayer graphene systems [20][21][22][23].
In Fig. 2 we show LDOS calculated in the BLG system for four values of the gate voltage V .Of the four gapless states that may appear at all at K = 0 in the BLG with a stacking domain wall [13], there is little left in the energy gap.For V = 0.4 eV, four separated flat bands are visible as a remnant of interacting topological gapless states (two of them are situated very close to the band continua).This is caused by the atomic-scale crack in the upper layer.For smaller V , fewer bands are visible in the energy gap.Finally, for V = 0.025 eV, the band that "grows up" from the VB continuum connects with the conduction band and the system takes on a metallic character along the domain wall.
Figure 3 shows LDOS for FLG system.Now, the situation is more complex.The stacking order change may, in principle, lead to even eight gapless states at K = 0, but here the number of bands in the energy gap depends strongly on V .For very small V = 0.025 eV the overall picture resembles two pairs of intersecting topological states with opposite slopes.They slightly repeel at the crossing points, so a small gap opens at the E F .For higher values, e.g., V = 0.075 eV, this small gap closes, the band crossing points approach the cone center and two new bands detach from the VB and CB continua.Increasing farther the gate, the gap opens (for V ≈ 0.1 eV), but closes again for V higher than 0.3 eV and the system restores again its metallic character.Especially interesting is the case when V ≈ 0.092 eV.The corresponding LDOS is presented in Fig. 4 (a).The four bands in the energy gap show flat areas and a short flat (in the range of k) band is formed also at the Fermi level where two bands touch.The increased density of this flat band at E = 0.14 eV is well seen in Fig. 4, where the projected density of states, i.e., PDOS(E)= k∈IBZ LDOS(E, k)dk is shown.The distribution of band densities in the layers is visualized in Fig. 4 (c) and (d).Of the two bands that touch at the Fermi level, the band below the E F is localized in two lower layers, while the band above the E F localizes in the two upper layers.Slightly away from the cone center, i.e., for k = 0.1, the distribution over layers is uniform.However, for k = 0, i.e., just at the flat band, a fraction of the charge moves to the outer layers.The similar effect of leakage of charge to the outer layers was recently observed also in [23] for a system of double twisted trilayers.
In low temperatures the flat band at E F is only halffilled.Note, that in graphene systems such a band can give rise to magnetic effects [24][25][26][27] and superconductivity [28][29][30][31].The width of the flat band at E F is about 5 meV.In Fig. 5 we show the LDOS(E,k = 0) of this band resolved into the layers, for two values of the gate voltage, V = 0.091 eV and V = 0.093 eV, for which the falt band is still present at E F .This result confirms that at k = 0 the density is localized mainly in the outer layers (blue and black).More importantly, we find that a slight increase of V by only a few mV, moves the charge occupying the flat band at E F , from the top to the bottom layer.The calculations have been also performed for trilayer (Fig. 1 (b)), but since the results do not differ substantially from those for BLG, they are not presented.We have also checked that wider separation of the halfgraphene layers does not essentially change the results.The band crossing points and the separation of bands can occur for different voltages, but the overall picture of the gap states and their behavior with the changing V is similar.

IV. CONCLUSIONS
We have studied the electronic structure of gated bilayer and four-layer graphene with stacking order domain walls created along the armchair direction.Bernal and rhombohedral stackings have been considered for BLG and FLG, respectively.Since the change of stacking across the armchair direction causes shear strain at the domain wall, we allow some layers to crack, i.e., to break C-C bonds along the domain wall to release the strain.
We have investigated the local density of states along the domain wall, i.e., at the interface of the adjacent stacking orders.The graphene energy cones E(k) overlap in the corresponding direction in the k-space.Thus, all the possible topological gapless states appear at the same valley, K = 0. Since the broken bonds act as atomic-scale defects and lead to the removal of valley protection of gapless topological states, these states interact and split yielding an effective energy gap.We have, however, found that for some gate voltages, the metallic character of the multilayer along the domain wall is restored.In the case of BLG, it happens for very small gates.In the case of the four-layer graphene, the gap states cross or anticross, i.e., they close or open the energy gap, depending on the gate.
In the case of four-layer graphene and for a specific value of the gate voltage, some bands in the gap connect forming a kind of a short flat band at the Fermi energy.We have demonstrated that the density of this band is moved to the outer layers.Finally, we have also shown that for very small changes in the gate voltage, the charge that fills the flat band at the E F , can change its localization between the bottom and the top layers.This gives us a tool for changing the layer localization of one-dimensional currents that at E F can flow along the stacking domain wall.

FIG. 1 :
FIG. 1: Schematic representations of the investigated systems.(a) BLG with stacking order change AB/BA.(b) TLG with stacking orders ABC/BAC.(c) FLG with stacking order change ABCA/BACB.(d) Lateral view of sequences of carbon atoms in the consecutive layers along the armchair direction.Vertical broken lines mark the interface between different stackings.The layers from bottom to top are marked in colors: blue, red, green, and dark grey.The bottom (blue) layer in (c) is covered and not seen.
(a)-(c), connected on the left-and right-hand sides to two semiinfinite multilayers.The local density of states is calculated as LDOS(E) = − 1 π Im[TrG C (E)], where G C is the interface Green function determined by the Hamiltonians of the conductor and the leads.The Hamiltonians are calculated using a tight-binding approximation.The values of intra-layer and inter-layer hopping parameters are t i = 2.7 eV and t e = 0.27 eV, respectively, following the[17][18][19][20].Since the system is periodic in the armchair direction, the local density of states is k-dependent, i.e., LDOS(E,k), where k is the wave vector corresponding to this periodicity.III.RESULTS AND DISCUSSIONWe perform calculations for defined graphene multilayers in a perpendicular external electric field.Several values of the field are considered, namely, E=7, 22, 44, and 120 mV/ Å.Assuming the interlayer distance d = 3.35 Å the values of E translate into gate voltages between the neighbor layers V = 0.025, 0.075, 0.15, and 0.4 eV.

FIG. 2 :
FIG. 2: LDOS for BLG with stacking order domain wall along the armchair direction, as visualized in Fig. 1 (a), for different gate voltages V applied between the layers.(a) V = 0.025 eV, (b) V = 0.075 eV, (c) V = 0.15 eV, (d) V = 0.4 eV.The Fermi level lies in the middle of the gap defined by the energy continua.

FIG. 3 :
FIG. 3: LDOS for FLG with stacking order domain wall along the armchair direction, as visualized in Fig. 1 (c), for different gate voltages V .(a) V = 0.025 eV, (b) V = 0.075 eV, (c) V = 0.15 eV, (d) V = 0.4 eV.The Fermi level lies in the middle of the energy gap.

FIG. 4 :
FIG. 4: (a) LDOS for FLG and V = 0.092 eV.(b) PDOS, i.e., LDOS integrated over k in the IBZ.(c) and (d) LDOS distribution in the layers for k = 0.1 and k = 0, respectively.The density in consecutive layers, from bottom to top, is shown in colors: blue, red, green, and black.