Asymmetry-enriched electronic and optical properties of bilayer graphene

The electronic and optical response of Bernal stacked bilayer graphene with geometry modulation and gate voltage are studied. The broken symmetry in sublattices, one dimensional periodicity perpendicular to the domain wall and out-of-plane axis introduces substantial changes of wavefunctions, such as gapless topological protected states, standing waves with bonding and anti-bonding characteristics, rich structures in density of states and optical spectra. The wavefunctions present well-behaved standing waves in pure system and complicated node structures in geometry-modulated system. The optical absorption spectra show forbidden optical excitation channels, prominent asymmetric absorption peaks, and dramatic variations in absorption structures. These results provide that the geometry-modulated structure with tunable gate voltage could be used for electronic and optical manipulation in future graphene-based devices.

Scientific RepoRts | (2019) 9:859 | https://doi.org/10.1038/s41598-018-37058-9 highly distorted energy dispersions with an eye-shaped stateless region accompanied by saddle points. In twisted BLGs, the moiré superlattice enlarges the primitive unit cell and thus creates a lot of 2D energy subbands. The geometry-modulated BLG is expected to generate many characters different from 2D bulk properties. Many works have studied band structures of AB-BA BLG by effective models 19,28,30 or tight-binding model with sharp DW 18,31 . Finite but short width has been studied in ref. 18 by first-principle calculation. However, the width of domain should be large enough to avoid interaction between different DW states. Reference 32 has studied DW states in large domain and DW width.
In this work, we investigate electronic and optical properties of geometry-and gate-modulated BLGs in tight-binding model. We present energy subbands, tight-binding functions on distinct sublattices, and DOS by exact diagonalization. Optical absorption spectra are calculated within gradient approximation. We study the effects caused by geometry modulation and gate voltage. The band structures show diverse 1D phenomena, including 1D energy subbands with various band-edge states, band splitting through geometry modulation, and metallic behavior in geometry-modulated system with a gate voltage. The wavefunctions present well-behaved standing waves in system without geometry modulation, while their node structures become complicated as applying geometry modulation and/or gate voltage. We discuss the wavefunctions of the topological protected DW states. For optical absorption spectra, we find forbidden optical excitation channels under specific linear relations between layer-dependent sublattices, prominent asymmetric absorption peaks in absence of selection rule, and DW-and gate-induced dramatic variations in optical absorption structures. We also discuss the connection to experiments. Our predicted results could be verified by the experimental measurements. Figure 1 shows schematic diagram for stacking DWs, also known as stacking solitons, between AB and BA domains in Bernal-stacked BLG. Such topological line defects involve both tensile and shear tilt boundaries, as a consequence of the competition between strain energy and misalignment energy cost within the DW, spanning in space a few to ten nanometers 20 . In this paper, we consider tensile soliton, which has zigzag boundary along the DW, created by uniformly stretching the C-C bonds of one layer along the armchair direction by an atomic bond distance. We assume the system has translational invariance along the tangent (ŷ) direction. For pristine Bernal-stacked BLG, all carbon atoms in a primitive unite cell are labeled by top-and bottom-layer indices, and A-and B-sublattice indices, (A 1 , B 1 , A 2 , B 2 ). The positions of these atoms can further be classified into even and odd for each sublattice. Because the wavefunctions with even and odd position indices only have a π phase difference, the considered bilayer system will be described by the components of wavefunctions only with odd position indices. In our notation, A 1 and A 2 (B 1 and B 2 ) are on top of each other for AB (BA) domain, as shown in Fig. 1(b).

Methods
The low-energy Hamiltonian is described by the single orbital p z tight-binding model, where c + (r) and c(r) are creation and annihilation operators at position r, respectively. The hopping amplitude is characterized by the empirical formula 9,28,33 where d = |r − r′| is the distance connecting two lattice points, b 0 = 1.42 Å the in-plane C-C bond length, d 0 = 3.35 Å the interlayer distance, and ρ = 0.184b 0 the characteristic decay length. γ 0 = −2.7 eV is the intralayer nearest-neighbor hopping integral and γ 1 = 0.48 eV is the interlayer interaction. We consider a system of 100  unite cells for each domain with various width of DWs (N dw strained unit cells with a modified lattice constant a′ = a(3N dw + 4)/(3N dw + 2) for geometry-modulated layer in each DW region, where = a b 3 0 is the lattice constant of single layer graphene). The cut-off of the distance for hopping amplitude is set to be d ∥ ≤ 2a. Because of translational invariance along zigzag direction, the quasi-1D Hamiltonian is a real symmetric matrix, which can possess real eigenfunction for each eigenvalue (If Hψ = Eψ, Eψ* = Eψ with real E for real symmetric H. When H does not have degeneracy, ψ = Eψ* real. When H have degeneracy, two of eigenfunctions can be chosen as ψ + ψ* ∈ real and i(ψ − ψ*) real. Therefore, all eigenfunctions of H can be real). Opposite on-site energy between two layers is included to describe layer-dependent Coulomb potential to simulate the effect of gate voltage 34 .
When BLG is subject to an electromagnetic field at zero temperature, occupied valence states are excited to unoccupied conduction ones by incident photons. In general, only vertical transitions can occur because of negligible momenta carried by photons. Based on the Kubo formula for linear response of an electromagnetic field, the optical absorption function is given by [35][36][37][38][39]  where → P is the momentum operator, Ê the unit vector of an electric polarization along ŷ, m e the bare electron mass, and Γ the broadening factor due to various de-excitation mechanisms. For a clean sample, physical properties can be observed under a sufficiently low temperature and Γ (=2 meV) will be small enough for observing the fine structures. The superscripts c and v represent the indices of conduction and valence subbands, respectively. The velocity matrix element, , for optical properties of carbon-related sp 2 bonding is evaluated from the gradient approximation 37,38 . The velocity operator is calculated by ∂H(k)/∂k y . The γ 0 -dependent velocity matrix elements dominate the optical excitations. Whether vertical excitation channels could survive is mainly related to the eigenfunction of the A-sublattice of the initial state to that of the B-sublattice of the final state within each layer, or vice versa. E c (k) − E v (k) is the excitation energy. The joint density of states (JDOS), which is defined by setting the velocity matrix elements in Eq. (3) equal to one, determine the available number of vertical optical excitation channels. While JDOS directly links with the band-edge states, it will exhibit the prominent structures.

Results and Discussion
Without geometry modulation. Low-lying band structures of pure and geometry-modulated BLGs exhibit unusual features. For pristine AB stacking graphene, the first conduction and valence bands are 2D parabolic bands and slightly overlap with each other around π π = ± ± K a a (2 / 3 , 2 /3 ). We first study the energy bands of AB stacking graphene without geometry modulation around K (k y a = 2π/3) in Fig. 2(a). 1D parabolic dispersions is due to zone-folding effects. Because of translational invariance along x direction, each band, except the first conduction and valence bands, in Fig. 2(a) involves two degenerate states, which can be referred to ±k x of the pristine AB stacking graphene. (The first conduction and valence bands are corresponding to high-symmetry line, therefore they are not degenerate.) Therefore, the band indices are closed related to discrete k x , provided applying Fourier transformation along the x direction. The energy dispersions along k x is almost negligible due to the sufficiently large unit cell. The particle-hole symmetry is broken because of non-vanished hopping integrals between atoms with distance longer than the spacing of the nearest neighbors. According to state energies measured from E F , the first valence and conduction bands (v 1 , c 1 ) touch with each other around K and have a very  Fig. 3(a), these two states are degenerate at K and can have other ratio between the weights of B 1 -and B 2 -sublattice through linear combination. Therefore, the wavefunctions at K are localized only at B 1 -and B 2 -sublattices, which are on top of each other. Notice that, because of armchair-sharp in the x direction, the position of atoms can further be classified by even or odd index. Within the same sublattice, the eigenfunctions with even or odd position index only has a π phase difference. For clear presentation, we only show the components with odd position indices. With increase of subband indices, the wavefunctions become well-behaved standing waves instead of uniform spatial distributions. For (v 2 , c 2 ) or (v′ 2 , c′ 2 ), the states have dominant components in the B 1 -and B 2 -sublattices and minor weights in A 1 -and A 2 -sublattices, which means that finite momentum away from K smear the wavefunction through finite weights in the sublattices with dangling atoms. The tight-binding functions show standing-wave behaviors with two nodes for each component in the spatial distributions, which is due to the linear combination between ±k x as mentioned before and the constrain to have real wavefunctions for a real symmetric matrix (If Hψ = Eψ, Hψ* = Eψ with real E for real symmetric H. When H does not have degeneracy, ψ = ψ* real. When H have degeneracy, two of eigenfunctions can be chosen as ψ + ψ* ∈ real and i(ψ − ψ*) ∈ real. Therefore, all eigenfunctions of H can be real). The π/2 phase for each sublattice between v 2 and v′ 2 or between c 2 and c′ 2 is consistent with uniform distribution for the bulk system. Our results also show that, the eigenfunctions between B 1 -and B 2 -sublattices are almost in-phase for conduction bands and out-of-phase for valence bands, and those between A 1 -and A 2 -sublattices are almost out-of-phase for conduction bands and in-phase for valence bands. The slight phase shift between A 1 and A 2 or between B 1 -and B 2 is intrinsic, depending on the material parameters. The phase shift becomes less obvious for modes with larger band indices. For the next band index, the tight-binding functions exhibit four-nodes standing waves (as shown in Fig. 3(g and h)). The number of nodes, Geometry modulation. Dramatic changes in electronic structure come to exist for the geometry-modulated BLGs (Fig. 1). Because of more complicated interlayer hopping integrals, asymmetry of the energy spectrum about Fermi energy is greatly enhanced by the geometric modulation, as clearly displayed in Fig. 2(b). The overlap of valence and conduction bands is getting larger, and so do the free electron and hole densities. The destruction of the inversion symmetry leads to the splitting of doubly degenerate states. More pairs of neighboring energy subbands are created. The splitting electronic states evolve more excitation channels and additional structures in optical absorption spectrum. Most of the energy bands have parabolic dispersions, while (v 1 ,c 1 ) and (v 2 ,c 2 ) around half filling exhibit oscillating and crossing behaviors. (The indices of valence and conduction bands are starting from the Fermi energy, as representatively shown in Fig. 2(b).) In general, the band-edge states in various energy subbands seriously deviate from K . They are responsible for the vHSs in DOS and thus the number of vertical optical transition channels. Spatial distributions of the wavefunctions belong to unusual standing waves. Symmetric and antisymmetric standing waves for each tight binding function thoroughly disappear under the modulation of stacking configuration, as clearly indicated in Fig. 4. The tight-binding functions on four sublattices roughly have a linear superposition relationship within the AB and BA domains. When the weight is large for one domain, it becomes small for the other. While the modulation of interlayer hopping amplitudes and relaxed intralayer hopping amplitudes act as scattering centers for the extended solutions of infinite system, the weight within DWs is small comparing with that in domains. Each component of a wavefunction still has fixed number of nodes. It grows with the 1D subband indices, e.g., two and four nodes in (v 1 /v′ 1 , c 1 /c′ 1 ) and (v 2 /v′ 2 , c 2 /c′ 2 ) subbands, respectively. Wavefunction with zero node disappears because of no translational invariance in the x direction. The position of nodes can be located within domains or DWs. The left (right) domain, which corresponding to BA (AB) stacking BLG, contains the major weight from the eigenfunctions of A 1 -and A 2 -(B 1 -and B 2 -) sublattices, consistent with the case without geometry modulation. For conduction (valence) bands, the eigenfunctions with major weight are in-phase (out-of-phase) and those with minor weight are out-of-phase (in-phase). On the other hand, for conduction (valence) bands, the eigenfunctions between A 1 and B 1 are symmetric (anti-symmetric) and anti-symmetric (symmetric) for c 1 and c 2 (v 1 and v 2 ), respectively. Similar behaviors of the eigenfunctions can be obtained for larger energy subband indices. Due to the major weight confined within one domain, the phase shift becomes less obvious, which is similar to the cases with nodes more than 2 in the uniform system.  Gate voltage. Electronic energy spectra are greatly affected by an external electric field. For pristine AB stacking BLG, the gate voltage (V z ) creates a band gap 8 . Roughly speaking, for AB-BA BLG, each domain still prefers to have a gap when applying a gate voltage. However, the characteristics of the wavefunctions in the AB domain is different from that in the BA domain, as discussed in the previous section. In order to bend in with these two subsystems, the geometry-modulation BLG with gate voltage will have metallic states within the gap. The reason is likely to be the two subsystems carry different valley Chern number 18 . Figure 2(c) shows rich and unique energy dispersions around the Fermi level. Most of energy levels repelled away from the Fermi energy by V z , where the energy dispersions become weak within a certain range of k y . Moreover, two Dirac dispersions around K lead to a finite contribution in DOS and provide metallic behavior near the Fermi energy, in great contrast with semiconducting behavior in the gated AB stacking BLGs. The valence and conduction wavefunctions become highly complicated, as displayed in Fig. 5, when the system sustains geometry modulation and gate voltage. In AB stacking BLG with a positive gate voltage in the z direction, which breaks inversion symmetry, the tight-binding functions of B 2 -sublattice are mainly responsible for the valence band states and those of B 1 -sublattice for conduction band states. This is different from the un-gated case, where the top layer and bottom layer components are distributed evenly between layers, as shown in Fig. 4. In other words, electron states are driven to the top layer, and hole states stay at bottom layer, and the BLG is therefore polarized under the electric field. This property of applying gate voltage is still maintained in the geometry-modulation BLGs. The left (right) column of the plots in Fig. 5, which corresponding to conduction (valence) bands, shows major contribution from the bottom (top) layer except the first row of the plots. We shall discuss these two states with more details later. As shown in Fig. 5(c-h), the abnormal standing waves show irregular features in oscillatory forms, amplitudes, numbers of nodes, and relationship among four sublattices. In general, there are no analytic sine or cosine waves suitable for the spatial distributions of the wavefunctions. Because the gated AB BLG has the band crossovers in both conduction and valence bands, the number of nodes does not grow with state energies monotonously. Furthermore, it might be identical or different for the top and bottom layers. For example, the v 1 (v 2 ) subband at K exhibits the 4-and 4-zero-point (4-and 6-zero-point) subenvelope functions on the top and bottom layers, respectively (Fig. 5(a and c)). If examining the wavefunctions from higher energy to lower energy (from Fig. 5 (h to b) of the conduction bands through (a) up to (g) of the valence bands) with fixed k y , we find that the major weights of the wavefunctions in corresponding domains seem to inject into the opposite domains through both DWs and have some rebounds. Similar behavior can be found in the opposite way. As a result, simple linear combination of the (A 1 , A 2 )-and (B 1 , B 2 )-related tight-binding functions disappear. These features of wavefunctions near Fermi energy are expected to induce very complex optical excitation channels for small energy frequency.
To compromise different characteristics in different DWs, the major components of the wavefunctions will leak into the other sublattices with energy outside the corresponding quasi-1D bulk bands. It is only possible to create this non-uniform distribution near DWs. A pair of special wavefunctions is shown in Fig. 5(a,b), which has significant enhancement for the weight within the DWs. These DW states are robust once the system has a gate voltage. The localization in Fig. 5(b) is weak, where the energy level at K is close to the band edges created by the gate voltage. When the k y is set to be near the crossing points (K L and K R ), the DW states become apparent, as shown in Fig. 6. In Fig. 6, instead of tight-binding functions, we present the probability distribution for each component in real space, which can be directly measured in experiments (discussed later). Note that, even we name these localized states DW states, the weight outside DW is still large and the corresponding decay length is around several nanometers for V z = 0.1 eV. The decay length becomes smaller as V z is larger. Therefore, large size of domains, as considered in this paper, is crucial to have well-define DW states when the applied gate voltage is limited. Two DW states, localized at each DW, belong to two different linear dispersive bands but with group velocities of the same sign. The DW states localized at the first DW have negative group velocities in ŷ and dominant weights for top (bottom) layer when k is around K L (K R ), as shown in Fig. 6(a,d,e,h)). On the other hand, the DW states localized at the second DW have positive group velocities in ŷ and dominant weights for bottom (top) layer when k is around K L (K R ), as shown in Fig. 6(b,c,f,g)). The dominant component of these DW states can be different from that of the other states which may be referred to bulk states in this gated geometry-modulated system. For example, probability distribution of A 1 component of the DW states in Fig. 6(f and g) appears in the right (AB) domain with countable weight in the DW, while that of bulk states is used to be located at the left (BA) domain. These unusual properties of DW states, which is presented in real space wavefunctions, might be used to distinguish topological protected states from the others.

Density of state.
Main features of electronic structures are directly reflected in the DOS, as shown in Fig. 7.
The pristine system has a finite DOS at E F (the dot line in Fig. 7(a)), indicating the gapless property. Weak overlap between conduction and valence bands leads to a very small shoulder structure at E F . In addition to a pair of shoulder structures around γ 1 , which comes from the valence and conduction bands of the second group, the absence of special structures near E F is consistent with the first-principles calculations 40 . In this work, we focus on the fine structure within small energy region and do not show the shoulder structures around γ 1 in Fig. 7. For a finite system without geometry modulation, Fig. 7(a) presents a lot of asymmetric peaks of DOS divergent in square-root form, due to 1D parabolic energy subbands as shown in Fig. 2(a). On the other side, the geometry modulation can create additional peaks ( Fig. 7(b-c)), which corresponds to the splitting of doubly degenerate states ( Fig. 2(b-c)). A finite DOS at E F in the geometry-modulated system shows metallic behavior, combined with a pair of rather strong peak structures around ±0.01 eV. The latter are related to the weakly dispersive energy bands near E F . With a gate voltage, the DOS near E F becomes a plateau structure, which is created by the linear energy dispersions across E F , as shown in Fig. 2(c). Two very prominent peaks appear around ±V z mainly due to the weak dispersions generated by the gate voltage. There are also some sub-structures for |E| < V z , which corresponds to band edge states with major weights around DWs. These accidental sub-structures come from large size of DWs.
optical absorption. The geometry-and gate-modulated BLG exhibits rich and unique optical properties.
The JDOS gives the weight of vertical optical excitations between different bands. For pristine AB stacking BLG, JDOS shows monotonically increasing within the considered frequency region 41 , which has the same reason to have monotonically increasing in DOS for small energy region away from E F . The inset of Fig. 8 shows peak structures, which indicate many channels between different 1D subbands, in the JDOS of BLG without and with geometry modulation (the dashed curves). Peaks in the JDOS are strong when the valence and conduction band-edge states have the same wave vector, e.g. (v 1 , c 1 , v 2 , c 2 ) around K and (v 2 , c 3 ) at k y a ≈ 2.06 for BLG with DWs, as shown in Fig. 2(b). There are some weak but observable structures due to non-vertical relation between different subbands. The peak structures in the JDOS are greatly reduced when system is under a gate voltage, as shown in the inset of Fig. 9, because most states with |E| < V z are repelled away from Fermi energy (Fig. 2(c)). For the system without geometry modulation, the threshold vertical channels is around V z = 0.1 eV. However, some weak structures appear in the geometry-modulated BLG for ω < V z because of linear dispersions and band edge states with major weights around DWs (Fig. 2(c)). This is one of major differences in optical absorption when geometry modulation is included. Structures in the optical absorption spectrum is determined by the available optical transition channels and the velocity matrix elements between channels. From Fig. 8, we can find both optical gap and energy gap are identical and equal to zero. The optical absorption spectrum also shows 1D absorption peaks. However, many peaks in JDOS do not appear in the optical absorption spectrum, as shown in the insets of Figs 8 and 9. The relations between wavefunctions of valence and conduction bands strongly affect the appearance of absorption peaks.
The v n to c n vertical excitations, which are due to the same pair of valence and conduction bands, are forbidden. The main mechanism is the linear symmetry or anti-symmetric superposition of the (A 1 , A 2 ) and (B 1 , B 2 ) sublattices, as discussed in Fig. 3. When the geometry modulation is introduced, the peaks are weaker as compared to uniform system and the curve for the optical absorption spectrum becomes smoother. The intensity, frequency, and structures of optical special absorption are very sensitive to the changes in the width of DW (Fig. 8) because of the shifting of band-edge states of the valence and conduction subbands. We find red-shift phenomena when increasing width of DW and suppressed strength for large ω. The first effect is mainly due to more atoms involved in the calculation and more energy subbands as a result. The second effect is related to larger band splitting for larger width of DW. Also, the first absorption structure might be replaced by another excitation channel during the variation of DW width.
The changes of the structures of optical absorption spectra become clear when applying various gate voltages, as shown in Fig. 9. The reduced intensity and the enhanced number of absorption structures occur in the increment of V z . This directly reflects gradual changes in energy dispersions, band-edge states and wavefunctions with the gate voltages. The first absorption peak, corresponding to (v 2 , c 3 ) at k y a = 2.06, for V z = 0 appears at 0.049 eV. First few absorption peaks are indicated by arrows above the curves with corresponding absorption channels in Fig. 9. By comparing the band structures and corresponding optical absorption spectra carefully, we find that those peaks involving v 2 (the corresponding band edges are sensitive to V z ) are shifted to higher frequency and become weaker when the gate voltage becomes larger. We find strong peaks which are contributed from multioptical excitation channels. Some weak absorption structures, which corresponds to new optical excitation channels, emerge in between as indicated by arrows below the curves in Fig. 9. New shoulder structures for frequency around 2V z become clearer for lager V z because of the band edges induced by the gate voltage. We note that the optical absorptions get greatly enhanced when the system has geometry modulation, as shown in the inset of Fig. 9. Most of the contribution is from the linear dispersions and band edge states with major weights around DWs (Fig. 2(c)).
Connection to experiments. The periodic boundary condition in asymmetry-enriched BLGs is responsible for 1D electronic and optical properties. The energy dispersions might belong to linear, parabolic and oscillatory forms, in which the first ones are the well-known vertical/non-vertical Dirac cones in the AA/AA′ stacking 9,27,42 . The band-edge states, which are the critical points in energy dispersion, correspond to the extreme points, saddle points, and effective 1D constant-energy loops. DOS presents various structures, including V-shape structures from linear energy dispersions, shoulder structures from 2D band edges, logarithmically symmetric peaks from the vHSs, and square-root asymmetric peaks from 1D band edges. These four kinds of special structures can further be revealed in the optical absorption spectra 43 . An electric field would create optical and energy gaps in most of BLG systems 7,8,12 . The experimental measurements can verify the predicted band structures, DOS, wavefunctions, and optical absorption spectra. The high-resolution angle-resolved photoemission spectroscopy (ARPES) can be used to directly examine the energy dispersion. The measured results have confirmed feature-rich band structures of the carbon-related sp 2 -bonding systems. Graphene nanoribbons are identified to possess 1D parabolic energy subbands centered at the high-symmetry point, accompanied with an energy gap and non-uniform energy spacings 44 . Recently, a lot of ARPES measurements are conducted on few-layer graphenes, covering the linear Dirac cone in monolayer system [45][46][47] , two pairs of parabolic bands in AB stacking BLG 45,48 , the coexistent linear and parabolic dispersions in symmetry-destroyed bilayer systems 49 , one linear and two parabolic bands in tri-layer ABA stacking 45,50 , the linear, partially flat and sombrero-shaped bands in tri-layer ABC stacking 50 . The Bernal-stacked graphite possesses 3D band structure, with the bilayer-and monolayer-like energy dispersions, respectively, at K and H points of the first Brillouin zone. The ARPES examinations on the geometry-modulated BLGs could provide unusual band structures, such as, the split of electronic states, diverse energy dispersions near E F , band-edge states, and metallic properties. These directly reflect the strong competition between stacking symmetries, interlayer hopping integrals and Coulomb potential.
The scanning tunneling spectroscopy (STS) measurements, in which the differential conductance (dI/dV) is proportional to DOS, are very powerful in exploring the vHSs due to the band-edge states and the metallic/semiconducting/semi-metallic behaviors. They have successfully identified diverse electronic properties in graphene nanoribbons [51][52][53] , carbon nanotubes 54,55 , few-layer graphenes [56][57][58][59][60][61][62] , and graphite 63,64 . Concerning graphene nanoribbons, the width-and edge-dominated energy gaps and the asymmetric peaks due to 1D parabolic bands are confirmed from the precisely defined crystal structures [51][52][53][54][55]65 . The similar prominent peaks obviously appear in carbon nanotubes, where they present chirality-and radius-dependent band gaps and energy spacings between two neighboring subbands 54,55 . A plenty of STS measurements on few-layer graphenes clearly reveal diverse low-lying DOS, including a V-shape dependence initiated from the Dirac point in monolayer system 64 , the peak structures closely related to saddle points in asymetric BLGs [57][58][59] , an electric-field-induced gap in bilayer AB stacking and tri-layer ABC stacking 60,61 , a pronounced peak at E F due to partially flat bands in tri-layer and penta-layer ABC stackings 54,55 , and a dip structure at E F accompanied with a pair of asymmetric peaks arising from constant-energy loops in tri-layer AAB stacking 54 . The measured DOS of the AB-stacked graphite is finite near E F characteristic of semi-metallic property 64 and exhibits the splitting π and π * strong peaks at deeper/ higher energy 63 . The focuses of the STS examinations on the geometry-modulated 22 and gated BLGs should be the square-root asymmetric peaks, the single-or double-peak structures, finite DOS at E F , and strong valence and conduction peaks caused by gate voltage. The STS can also be used to measure the two-dimensional structure of individual wavefunctions in metallic single-walled carbon nanotubes 66,67 . Therefore, the wavefunctions found in the paper can be verified by STS measurements.
Up to date, four kinds of optical spectroscopies, including absorption, transmission, reflection, and Raman scattering spectroscopies, are frequently utilized to accurately explore vertical optical excitations 43 . Concerning the AB-stacked BLG, the experiments have confirmed the 0.3-0.4 eV shoulder structure under zero field, the V z -created semimetal-semiconductor transition and two low-frequency asymmetric peaks, the two strong π-electronic absorption peaks at the middle frequency, specific magneto-optical selection rule for the first group of Landau levels (LLs), and linear magnetic-field-strength dependence of the inter-LL excitation energies. Similar verifications performed on trilayer ABA stacking cover one shoulder around 0.5 eV, the gapless behavior unaffected by gate voltage, the V z -induced low-frequency multi-peak structures, several π-electronic absorption peaks, and monolayer-and bilayer-like inter-LL absorption frequencies. Moreover, the identified spectral features in trilayer ABC stacking are two low-frequency characteristic peaks and gap opening under an electric field. The above-mentioned optical spectroscopies are available in examining the vanishing optical gaps for any metallic systems, prominent asymmetric absorption peaks, absence of selection rule, forbidden optical excitations associated with linear relations in the (A 1 , A 2 ) and (B 1 , B 2 ) sublattices, and the variations in absorption structures due to the modulation of DW width and gate voltage.

Concluding Remarks
We have studied the electronic and optical properties of AB stacking BLG with geometry modulation and gate voltage in tight-binding model. We present energy subbands, tight-binding functions on distinct sublattices, and DOS by exact diagonalization. Effects of geometry modulation in the presence or absence of a gate voltage are discussed. The metallic system exhibits a plenty of 1D energy subbands, accompanied by well-behaved or irregular standing waves. Specifically, the layer-dependent Coulomb potential destroys a simple relation between the (A 1 , A 2 ) and (B 1 , B 2 ) subenvelope functions and gives complicated node structures. With a gate voltage, the system is still metallic due to the existence of the DW states. The wavefunctions of the topological protected DW states present unusual space distributions. DOS shows various vHSs, including single-and double-peak structures with the square-root divergent forms, a pair of prominent peaks caused by gate voltage, and a plateau structure across E F .
Optical absorption spectra are calculated within gradient approximation. We find forbidden optical excitation channels under specific linear relations between layer-dependent sublattices, prominent asymmetric absorption peaks in absence of selection rule, and DW-and gate-induced dramatic variations in optical absorption structures. Concerning the geometry-modulated systems, the observable absorption peaks could survive only under the destruction of symmetric or anti-symmetric linear superposition of the (A 1 , A 2 ) and (B 1 , B 2 ) sublattices. Simple dependence of absorption structures on the modulation width of DW is absent. However, the reduced intensity and the enhanced number are regularly revealed in the increment of V z . The frequency, number, and form intensity of optical absorption peaks strongly depend on the modulation period and electric-field strength. Our predicted results could be verified by the experimental measurements. This geometry-modulated system is suitable for studying various physical phenomena. For example, magneto-electronic properties, which exist in a uniform perpendicular magnetic field, might be dominated by the 1D Landau subbands, being in sharp contrast with the Landau levels in 2D systems. This problem is under current investigation. How to generalize the manipulations of geometric structures to emergent layered materials is the near-future focus, e.g., important studies on essential properties of bilayer silicene 68 , germanene 69 , phosphorene 70 , and bismuthene 71 .