Band gap opening in graphene: a short theoretical study

Abstract

Graphene, being a gapless semiconductor, cannot be used in pristine form for nano-electronic applications. Therefore, it is essential to generate a finite gap in the energy dispersion at Dirac point. We present here the tight-binding model Hamiltonian taking into account of various interactions for tuning band gap in graphene. The model Hamiltonian describes the hopping of the π-electrons up to third nearest-neighbours, substrate effects, Coulomb interaction at two sub-lattices, electron–phonon interaction in graphene-on-substrates and high phonon frequency vibrations, besides the bi-layer graphene. We have solved the Hamiltonian using Zubarev’s double time single particle Green’s function technique. The quasi-particle energies, electron band dispersions, the expression for effective band gap and the density of states (DOS) are calculated numerically. The results are discussed by varying different model parameters of the system. It is observed that the electron DOS and band dispersion exhibit linear energy dependence near Dirac point for nearest-neighbour hopping integral. However, the second and third nearest-neighbour hoppings provide asymmetry in DOS. The band dispersions exhibit wider band gaps with stronger substrate effect. The modified gap in graphene-on-substrate attains its maximum value for Coulomb interaction energy UC=1.7t1 . The critical Coulomb interaction is enhanced to UC=2.5t1 to produce maximum band gap in the presence of electron–phonon interaction and phonon vibration. The bi-layer graphene exhibits Mexican hat type band gap near Dirac point for transverse gating potential. The other conclusions for the present work are described in the text.

Introduction

Graphene is one-atom-thick two-dimensional structure with carbon atoms packed in a honeycomb lattice. Its recent experimental discovery has stimulated extensive investigations on every aspect of this novel material [ 1 , 2 ]. The tight-binding calculation for graphene shows that its conduction and valence bands touch at six Dirac points in the Brillouin zone [ 3 ] where energy dispersions are linear with respect to momentum. This unique band dispersion in graphene leads to graphene’s novel physical and electronic properties such as room temperature quantum Hall effects and high charge carrier mobility [ 4 – 6 ]. Graphene, being a gapless semi-metal, cannot be used in pristine form for nano-electronic applications. Therefore, it is necessary to open a finite gap in the energy dispersions at K point by various mechanisms [ 7 , 8 ]. When graphene lies on born nitride (BN) substrate, a small gap (∼100 meV) is observed [ 9 , 10 ]. Similarly, a band gap of 250 meV is observed for silicon carbide substrate [ 11 , 12 ]. There is a great effort for producing a gap (∼1 eV) which is observed in germanium/silicon for the application in digital electronics.

The recent work indicates that the carbon layer is co-covalently bonded to the SiC sub-lattice [ 13 – 17 ]. The band gap opening for graphene-on-ruthenium metal also arises due to symmetry breaking in the system [ 18 , 19 ]. McCann and Falko [ 8 , 20 ] have proposed that bi-layer graphene can develop a gap, when transverse electric field is applied between two surfaces of the system acting as a gating. Band gap in bi-layer graphene has been observed experimentally by infrared spectroscopy [ 21 , 22 ] and angle resolved photo emission spectroscopy (ARPES) [ 23 ]. Hague [ 24 , 25 ] has proposed a theoretical model calculation taking attractive interaction modified through phonons in polarisable substrate for the strong enhancement of the band gap in graphene substrates. The band gap of several orders in eV can be prepared in Graphane [ 26 ] and Fluorographene [ 27 , 28 ] by chemical modification of the system with hydrogen and fluorine, respectively. Atomic thick boron nitride (BN) forms a honeycomb lattice where the π orbitals on N sites are shifted up in energy by +∆ and decreased in energy of −∆ on B site causing a gap of 2∆ [ 29 ]. The band gap of 5.56 eV is observed experimentally on monolayer BN systems.

The role of Coulomb interaction in low-dimensional systems provides unique opportunity for theoretical as well as experimental studies. The graphene [ 30 , 31 ], varieties of semiconductor surfaces like Si, Ce, Sn, Pb [ 32 ], Bechgard salt [ 33 ] and doped polymers [ 34 , 35 ] display strong on-site as well as inter-site Coulomb interactions. It is observed that on-site Coulomb interaction in graphene is U  = 3.3 t 1 and the nearest-neighbour Coulomb interaction is V=2.0t1 where the nearest-neighbour hopping integral is t 1  = 2.8 eV [ 36 ]. It has been reported that the effective Hubbard interaction is U=3.3t1 in the close vicinity of the separation between conducting graphene and insulating phases [ 37 , 38 ]. The Coulomb interaction between massless fermions appears to be unscreened in pristine graphene. At present, it is not clear whether this type of Coulomb interaction would lead to weak correlation in graphene or strongly correlated electronic phases like an insulator [ 39 , 40 ]. The long ranged Coulomb interaction in pristine graphene leads to unusual behaviour [ 41 , 42 ]. The on-site Coulomb interaction plays a crucial role for the understanding of the defect induced magnetisation [ 43 – 45 ].

Earlier, we have reported the study of band gap opening in graphene by a single impurity taking the tight binding model up to the third nearest-neighbour hoppings in the absence of Coulomb interaction [ 46 ], in presence of Coulomb interaction [ 47 ], electron phonon interaction [ 48 ] and bi-layer graphene [ 49 ]. In the present brief review, we study the effect of all interactions on the band gap opening of graphene. We propose a tight-binding model Hamiltonian consisting of the site energy of the carbon atoms and nearest- neighbour hopping of π electrons of carbon atoms taking into account the substrate effects, Coulomb interaction, electron–phonon interaction effect on monolayer graphene and finally the effect of gate potential on the band gap of bi-layer graphene. The rest of the work is as follows. We describe the tight-binding model Hamiltonian for different types of interaction. We briefly outline the single particle Green’s function to calculate correlation functions, physical parameters, quasi-particle energies and electron density of states in the next section. Then we present the results and discussion and finally the conclusion.

Tight-binding model

The graphene sheet is formed by carbon atoms arranged in a non-Bravais honeycomb lattice with nearest-neighbour C–C distance of a0=1.43 Å where the lattice constant is a=3a0 . The s, px and py orbitals hybridise to form sp2 bonds leading to high energy sigma bonds. The pz orbitals in graphene form the π bond which is responsible for the electronic properties of graphene. The 2D character of the honeycomb lattice does not allow the overlap of the pz orbital of carbon atom and s, px and py orbitals of its nearest-neighbour carbon atoms. One can construct a simple tight-binding model incorporating only the nearest-neighbour hopping of electrons between the adjacent A and B sub-lattices of honeycomb lattice to study the low energy electron excitations and doping effect in pristine graphene.

Assuming the electron hopping up to the third nearest-neighbours, the tight-binding Hamiltonian for electron in graphene is written as

H0=∑i,σϵaai,σ†ai,σ+ϵbbi,σ†bi,σ-t1∑<i,j>,σai,σ†bj,σ+bj,σ†ai,σ-t2∑≪i,j≫,σai,σ†aj,σ+bj,σ†bi,σ-t3∑≪<i,j≫>,σai,σ†bj,σ+bj,σ†ai,σ,\sigma }} \left( {a_{i,\sigma }^{\dag } b_{j,\sigma } + b_{j,\sigma }^{\dag } a_{i,\sigma } } \right) - t_{2} \mathop \sum \limits_{ \ll i,j \gg,\sigma } \left( {a_{i,\sigma }^{\dag } a_{j,\sigma } + b_{j,\sigma }^{\dag } b_{i,\sigma } } \right) - t_{3} \mathop \sum \limits_{{ \ll < i,{\text{j}} \gg { > },\sigma }} \left( {a_{i,\sigma }^{\dag } b_{j,\sigma } + b_{j,\sigma }^{\dag } a_{i,\sigma } } \right)$$\end{document}]]>
where ai,σ†(ai,σ) creates (annihilates) an electron with spin σ ( σ  = ↑,↓) on site R→i on sub-lattice A. Similarly, bi,σ†(bi,σ) creates (annihilates) an electron on sub-lattice B. Here, t1=2.5to3.0eV is the nearest-neighbour hopping energy, t 2 with 0.02t1≤t2≤0.2t1 [ 36 , 37 ] is the next-nearest-neighbour hopping energy, t3 is the third nearest-neighbour hopping integral with t3=0.024t1eV and εa(εb) is the site energy at the sub-lattice site A (B). Further, <i,j>, ,$$\end{document}]]>  ≪  i ,  j  ≫,  and ≪ <  i ,  j  > ≫  stand for nearest, next-nearest and next–next-neighbour hoppings from site R→i to R→j .

The Fourier transformed dispersion γ 1 ( k )for the nearest-neighbour hopping is

γ1k=eikxa0+2e-i2kxa0·cos32kya0
and dispersions for the next-and next-to-next-nearest-neighbour hoppings are γ2k=∑δ2eik→.δ→2 and γ3k=∑δ3eik→.δ→3 . where δ2→ and δ3→ are the corresponding lattice vectors. Graphene deposited on silicon dioxide (SiO 2 ) substrate can be described by the 2D massless Dirac fermions [ 3 ]. On the other hand, graphene deposited on silicon carbide (SiC) substrate can be described by massive 2D Dirac electron. A substrate-induced potential thus generated can break the symmetry of the honeycomb lattice and generate gap in electron system. In this case, the electron interacts with the static potential induced by the substrates and hence A sub-lattice site acquires an energy +Δ and the B site with energy -Δ leading to the symmetry breaking between A and B sites. Such a symmetry breaking Hamiltonian is written as
Hsub=Δ∑i,σai,σ†ai,σ-Δ∑i,σbi,σ†bi,σ.

The system exhibits an insulating ground state with a band gap 2Δ , while Fermi level lies in the gap in the un-doped system. The Coulomb interaction forbids both the electron occupancies at the same site. The Hubbard interaction representing the two sub-lattices with an effective Coulomb energy U is written as

HU=U∑ini,↑ani,↓a+ni,↑bni,↓b
where n i ,↑ α ( n i ,↓ α ) with α  ∊  A ,  B sub-lattices, represents the occupation number with up(down) spin. For weak coupling, the Hamiltonian can be decoupled by Hartree–Fock mean-field approximation decoupling scheme, i.e. Uni,↑αni,↓α≈U<ni,↑α>ni,↓α+U<ni,↓α>ni,↑α-U<ni,↑α><ni,↓α> n_{i, \downarrow }^{\alpha } + U < n_{i, \downarrow }^{\alpha } > n_{i, \uparrow }^{\alpha } - U < n_{i, \uparrow }^{\alpha } > < n_{i, \downarrow }^{\alpha } >$$\end{document}]]> , where α  ≡  a ,  b corresponding to A and B site interactions. The mean-field solutions are taken as <ni,↑a>+<ni↑b>2=n + < n_{i \uparrow }^{b} > }}{2} = n$$\end{document}]]> and <ni,↑a>-<ni↑b>2=d - < n_{i \uparrow }^{b} > }}{2} = d$$\end{document}]]> and this leads to the condition, <ni,↑a>=n+d = n + d$$\end{document}]]> and <ni,↑b>=n-d = n - d$$\end{document}]]> where n represents the mean electron occupation and d , the deviation from the mean occupation. Similar expression can be formed for the down spin electron. The electron–phonon interaction and lattice vibration are written as
He - ph=∑α,i,jfi-jni,σα(bj†+bj)
Hp=∑jω0bj†bj.
Here H e-ph describes electron–phonon interaction where the electron density n i , σ α at sub-lattice site is coupled to the phonon displacement uj=bj†+bj with phonon creation (annihilation) operator bj†(bj) at site j and f ( i − j ) is the electron–phonon coupling. For high frequency phonons, Eq. ( 5 ) represents the Holstein interaction. The Hamiltonian H p represents the free phonon energy with phonon frequency ω 0 .

Calculation of Green’s functions and quasi-particle bands

The Green’s functions for the electrons of A and B sub-lattices are calculated by Zubarev’s Green’s function technique [ 50 ]. The Zubarev’s double time single particle retarded Green’s function is defined

Grt,t′=≪AT;B(t′)≫r=-iΘt-t′<At;Bt′> $$\end{document}]]>
where <…> indicates the average over a grand canonical ensemble. AtandB(t′) are Heisenberg representations of the operators, while Θ = 1 for t  >  t ′ and Θ = 0 for t  <  t ′ . Here, [ A ( t );  B ( t ′ )] indicates commutation or anti-commutation depending on Boson or fermions operators. The poles of the Green’s functions provide quasi-particle energies.

The correlation function corresponding to Green’s function G r ( t , t ′ ) is defined as

<Bt′At>=∫-∞∞Jk(ω)e-iω(t-t′)dω = \mathop \int \limits_{ - \infty }^{\infty } J_{k} (\omega ) e^{{ - i\omega (t - t^{'} )}} {\text{d}}\omega$$\end{document}]]>
where the spectral density function J k ( ω ) is written as
Jkω=ilimη→01eωθ+1Grω+iη-Grω-iη
where θ  =  k B T with η as a small spectral width. Other physical parameters can be calculated using the correlation functions of the corresponding Green’s functions.

The density of states which is proportional to tunneling conductance is calculated from the imaginary part of Green’s function. The electron density of states for the 2D graphene is given by the formula

DOS=A(2π)2∫∫dk→xdk→yρk,σ(ω)
where A is the area of the honeycomb lattice in real space and spectral function ρ k , σ ( ω ) is given by
ρk,σω=-2π∑σImGrω+iη
with η as a small spectral width. Finally, the occupations and their difference for sub-lattice electrons for different spin orientations are calculated to study the magnetic effect of Coulomb interaction in the gap formation in graphene. All the energy parameters are scaled by the hopping integral t1 .

Results and discussion

The band gaps can be induced in graphene near Dirac point by the following techniques.

Effect of electron hopping

The electron density of states (DOS) and band dispersion are numerically computed [ 46 ] and are shown in Figs.  1 and 2 . The tight-binding calculations give the first-nearest-neighbour hopping integral t 1  = 2.5–3.0 eV [ 30 , 31 ]. In the present calculation, we have taken ɛ a  =  ɛ b  = 0, t1=-2.78eV , t2=-0.12eV , t3=-0.068eV . Here, the scaled hopping integrals become t~1=-1 , t~2=-0.043 and t~3=-0.024 . The density of states (DOS) for electrons for the graphene is plotted for different band energies (c) (see Fig.  1 ). The DOS exhibits V-shaped gap structure at K-point (Dirac point) for nearest-neighbour hopping energy t~1=-1 , i.e. DOS shows linear dependence of band energy. When second nearest-neighbour hopping ( t~2=-0.043 ) is included, it still retains the V-shape, but shifts to lower energies becoming asymmetric in nature with respect to Fermi level ( εF=0 ) at Dirac point (Fig.  1 ).

Fig. 1

Plot of the electronic density of states (DOS) of ideal graphene for different band energy with different hopping integrals, t~1=-1 ( solid line ) for NN; t~1=-1 , t~2=-0.043 ( dotted line ) for NNN, t~1=-1 , t~2=-0.043 , t~3=-0.024 ( dashed line ) for NNNN

Fig. 2

Plot of the electronic band dispersion of ideal graphene for different band energy with different hopping integrals, t~1=-1 ( solid line ) for NN; t~1=-1 , t~2=-0.043 ( dotted line ) for NNN, t~1=-1 , t~2=-0.043 , t~3=-0.024 ( dashed line ) for NNNN

The inclusion of third neighbour interaction induces a gap near K-point and the gap shifts further to lower energies. Figure  2 shows the electron energy dispersion for the graphene and is plotted for different band hopping energy. The energy band dispersion exhibits V-shaped nature at K-point (Dirac point) for nearest-neighbour hopping energy t~1=-1 , i.e. band dispersion shows linear dependence of band energy. When second nearest-neighbour hopping ( t~2=-0.043 ) is included, it still retains the V-shape, but shifts to lower energies becoming asymmetric in nature with respect to Fermi level ( ɛ F  = 0) at Dirac point (Fig.  2 ). The inclusion of third neighbour interaction induces a gap near K points and the gap shifts further to lower energies.

Substrate effect

To study the substrate effect in band gap opening in graphene, we compute DOS and electron band dispersion near Dirac point using the expressions given in Eqs. ( 1 ) and ( 3 ) of our earlier publications [ 46 ] and show the plots in Figs.  3 and 4 . Figure  3 shows the DOS for different values of band gaps d 1  = 0.035–0.107 developing in graphene sheet due to substrate effect [ 9 , 11 , 13 , 18 ] and [ 32 – 34 ]. With the onset of substrate effect, the Dirac point moves towards the valence band occupied by electrons and introduces a band gap below Fermi level ɛ F  = 0. The gap is further enhanced with the increase of d 1 . Figure  4 shows the enhancement of band gap between lower electron occupied valence band and upper hole band due to the increase of d 1 . The middle of the band gap also moves down the Fermi level ( ɛ F  = 0).

Fig. 3

Plot of the variation of density of states (DOS) with energy ( c ) for different substrates d 1  = 0.035, 0.071, 0.089, 0.0107

Fig. 4

Plot of the variation of energy band dispersion with momentum ( k ) for different substrates d 1  = 0.035, 0.071, 0.089, 0.0107

Effect of Coulomb interaction

The sub-lattice Coulomb interaction is treated within a mean-field approximation. Using the Hamiltonian given in Eqs. ( 1 ), ( 3 ) and ( 4 ) of our earlier calculations [ 47 ], we calculate the difference ( d ) in occupation numbers as

d=∑kΔ¯fβω1k↑-fβω2k↑(ω1k↑-ω2k↑)
where f is the Fermi distribution function and ω1k,↑ and ω 2 k ,↑ are the electron band dispersions. The modified gap d2T=Δ~t1 is plotted in Figs.  5 , 6 for different Coulomb energies and hopping parameters. The ferromagnetic magnetizations and spin polarizations are reported [ 51 ].

Fig. 5

Plot of the variation of modified band gap ( d 2 ) with temperature ( t ) for different values of Coulomb energy u  = 1.0, 1.5, 1.7, 2.0, 2.5, 3.5, 4.5 and 5.0 for fixed substrate-induced band gap d 1  = 0.1

Fig. 6

Plot of the variation of modified band gap ( d 2 ) with temperature ( t ) at u  = 1.7 for different values of electron hopping t~1=-1 , t~2=-0.043 , t~3=-0.024 for fixed substrate-induced band gap d 1  = 0.1. The band dispersions for different hopping are also given in inset of the figure

The graphene acquires a band gap of ∆ due to substrate effect. The effective band gap then becomes Δ¯=Δ+Ud2 due to Coulomb interaction between electrons. We scale the physical parameters by nearest-neighbour hopping integral ( t 1 ) and hence the modified gap appears as d2=d1+ud2 . The temperature dependent difference ( d ) in electron occupancies of A and B sub-lattices for up spin electrons is computed numerically and self consistently and hence the modified band gap ( d2=Δ¯t1 ) is computed for different Coulomb correlation energies as shown in Fig.  5 . For lower Coulomb energies, the modified gap ( d 2 ) at t  = 0 gradually increases with increase of u from 1 to 1.7 and attains the maximum of d2=0.182 . On further increasing to higher Coulomb energies, the modified gap ( d 2 ) at t  = 0 decreases and attains the bare gap ( d 1 ) arising due to only substrate effect for Coulomb energy u  = 5.0. It is observed that the magnetic gap vanishes for u  = 5.0 indicating that n↑a=n↑b in paramagnetic phase. The temperature dependence of modified gap shows that the magnitude of gap for lower ‘ u ’ gradually decreases with temperature. For higher u , the modified gap increases with temperature, attains its maximum and then decreases with temperature. However, the maximum of modified gap nearly remains same for higher value of ‘ u ’ at higher temperatures indicating that d=<n↑a>-<n↑b> - < n_{ \uparrow }^{b} >$$\end{document}]]> remains unchanged at very high temperatures and high Coulomb interactions.

Figure  6 shows the effect of different hopping integrals on temperature-dependent modified gap ( d 2 ). In the absence of Coulomb interaction and electron hopping, the band gap in graphene is d 1  = 0.1 which arises due to substrate effect only. This band gap at temperature 0 K is enhanced to d 2  = 0.167 due to the onset of a critical Coulomb energy u c  = 1.7 for the nearest-neighbour hopping energy of t 1  = −1.0. This gap at temperature 0 K is further enhanced to d 2  = 0.181 for the same critical Coulomb energy u c  = 1.7 and for hopping integrals taken up to third nearest neighbours. It appears that the contributions of third nearest-neighbour hopping integrals and beyond have little effect on the band gap. The temperature-dependent modified gap shows that the gap is the highest at very low temperatures and gradually decreases with increase of temperature. This effect of hopping integrals on the modified gap is also seen in the band dispersions shown in the inset of the Fig.  6 .

Effect of electron–phonon interaction

The electrons on graphene sheet interact strongly with the phonons on the polarized surface of substrates. Applying Lang–Firsov canonical transformation [ 48 ] to Eqs. ( 1 ), ( 3 ), ( 4 ), ( 5 ) and ( 6 ) in high frequency limit of localized phonons, we calculate the effective Coulomb energy U~=(U-2λt~1) and effective hopping ( t~1 ) which is a function of phonon frequency ω 0 , electron–phonon coupling ( λ ) and NN hopping integral ( t~1 ). The temperature-dependent modified gaps are plotted in Figs.  7 and 8 .

Fig. 7

Plot of the variation of modified band gap ( d 2 ) vs. phonon frequency ω 0 for different values of Coulomb potential u=0.5,1.0,1.5,2.0,2.5 for fixed electron–phonon coupling constant al=0.2 , temperature (t)=0.01 and substrate-induced gap d 1  = 0.1

Fig. 8

Plot of the variation of modified band gap ( d 2 ) vs. temperature ( t ) for different values of phonon frequency ω0=0.2,0.4,0.6,0.8,1.0 for fixed electron–phonon coupling constant al  = 0.2, Coulomb interaction (u)=1.7 and substrate-induced gap ( d 1 ) = 0.1

The effect of Coulomb interaction ( u ) on the modified gap d 2 for high phonon frequency ( ω 0 ) vibration is shown in Fig.  7 . For given value of lower Coulomb interaction, the modified gap gradually increases with phonon frequency. With further increase of Coulomb energy, the modified gap gradually increases from d 1  = 0.1 and attains maximum value d 2  = 0.165 for given Coulomb interaction u  = 2.0. With further increase in Coulomb interaction, the modified gap gradually increases with phonon frequency ω 0 attains a maximum flat peak with d 2  = 0.165 for critical Coulomb interaction u c  = 2.5 for critical phonon frequency ω 0c  ≈  0.35. Thus, it is clear that modulated gap is maximized for critical Coulomb interaction u c  = 2.5 and critical phonon frequency ω 0c  = 0.35 for given values of electron–phonon coupling al  = 0.2. The critical Coulomb interaction U C  = 2.5 t 1 for producing maximum modified gap is slightly higher than the critical Coulomb interaction Uc=uct1=1.7t1 obtained earlier in the absence of electron–phonon interaction and phonon frequency [ 47 ].

The effect of phonon frequency ( ω 0 ) on gap is shown in Fig.  8 . For a given low electron–phonon coupling al=0.2 and relatively low phonon frequency ω 0  = 0.2, the modified gap becomes d 2  = 0.13 gap at temperature t  = 0. With increase of phonon vibrational frequency, the modified gap is enhanced to the higher value, i.e. d 2  = 0.165 at temperature t  = 0 for vibration frequency ω 0  = 1. However, for a given phonon vibration frequency, the modified gap decreases with temperature and remains nearly constant at higher temperatures. Thus, it is clear that phonon vibration frequency enhances the gap near room temperature.

Effect of bi-layer graphene (BLG)

Experiments [ 8 , 20 , 23 ] show the evidence of band gap opening in graphene in BLG by gating between two layers. We have proposed model Hamiltonian consisting of intra-layer and inter-layer hopping integrals t 1 and t ⊥ , respectively, in the presence of gating potential V [ 49 , 52 , 53 ]. The four bands for BLG are written as

ωα,s,k=-μ-(-1)s[V22+2εk2+εk,⊥2-(-1)αR]2
R=V2/2+2εk2+εk⊥22-4V2/16+εk2V2/2+εk⊥2V2/4+εk4
where s,α=1,2 . The band dispersion is plotted in Fig.  9 which shows four bands ωα,sk of which ω11k and ω 22 k are high energy bands and ω 12 k (conduction) and ω 21 k (valence) bands are the low lying dispersion bands near Dirac point with zero energy. In the absence of gate potential v=Vt1=0 , low energy bands touch at Dirac point. For low electric potential ( v=0.054 ), the low energy bands show a Mexican hat shape with a gap energy of 0.054 as observed experimentally [ 8 , 23 ].

Fig. 9

Plot of energy band dispersions vs. momentum (k) for pristine bi-layer graphene for different values of external biasing v  = 0 ( dashed line ) and v  = 0.054 ( solid line )

Conclusions

The pristine graphene is a gapless semiconductor. It is necessary to induce a band gap in graphene near the Dirac point for its application in several electronic devices. We present here different model calculations for band gap opening in monolayer and bi-layer graphenes taking electron hoppings up to third nearest-neighbours, substrate effect, on-site Coulomb correlation effect and finally electron–phonon interaction along with lattice vibration effect. All the calculations are done using Zubarev’s double time single particle Green’s function. The expression for modified band gap, dispersions and density of states (DOS) are calculated and computed numerically. The results are discussed varying different model parameters. The DOS exhibits a V-shaped gap near Dirac point with linear energy dependence for nearest-neighbour hopping t 1 . The DOS retains its V-shape for second and third nearest-neighbour hoppings except that the DOS shifts towards low energy valence band. Similar result is concluded from band dispersions for different hopping integrals.

The DOS exhibits wider gap near Dirac point with the increase of band gap arising due to substrate effect. Similar conclusion is derived from the band dispersions also. The substrate-induced gap is modified due to Coulomb interaction energy. The modified gap is enhanced at temperature T  = 0 0 K with the increase of Coulomb interaction. The effective gap attains its maximum value, d 2  = 0.182 at t  = 0 for critical Coulomb interaction U C  = 1.7 t 1 . On further increase in the Coulomb energy to U=5.0t1 , we observe that the modified gap d 2 = d 1 indicating that the electron occupancies of both sub-lattices becoming same (i.e. n↑a=n↑b ) in paramagnetic phase. We have calculated the modified band gap taking the electron–phonon interaction between the graphene layer and the substrate in high frequency limit of the localized phonon. In the presence of electron–phonon interaction and phonon vibration, the modified gap attains its maximum value ( i.e.d2=0.165 ) for a critical Coulomb interaction Uc=2.5t1 which is higher than the critical Coulomb interaction Uc=1.7t1 obtained earlier in the absence of electron–phonon interaction and phonon vibration [ 47 ]. Here, we have proposed tight-binding model calculation for AB-stacked bi-layer graphene taking intra-layer and inter-layer hopping in the presence of transverse gating potential. We have obtained four quasi-particle bands of which the inner two bands touch each other at Dirac point. A band gap opens near Dirac point due to application of the gating potential. These two inner bands exhibit Mexican hat shape for gating potential V=0.054t1 . Since we have considered dispersion for full Brillouin zone in both the layers, we can take any gating potential higher than the gating potential V=0.054t1 applied to bi-layer systems experimentally [ 8 , 23 ].

Acknowledgements

The authors gracefully acknowledge the research facilities offered by the Institute of Physics, Bhubaneswar, India.

References

  1. Novoselov and Geim (2004) Electric field effect in atomically thin carbon films 10.1126/science.1102896
  2. Novoselov et al. (2005) Two-dimensional gas of massless Dirac fermions in graphene 10.1038/nature04233
  3. Wallace (1947) The Band Theory of Graphite 10.1103/PhysRev.71.622
  4. Miller and Kubista (2009) Observing the Quantization of ZeroMass Carriers in Graphene 10.1126/science.1171810
  5. Stander et al. (2009) Evidence of Klein tunneling in graphene pn junctions 10.1103/PhysRevLett.102.026807
  6. Bolotin et al. (2008) Temperature-Dependent Transport in Suspended Graphene 10.1103/PhysRevLett.101.096802
  7. Han and Özyilmaz (2007) Energy Band-Gap Engineering of Graphene Nanoribbons 10.1103/PhysRevLett.98.206805
  8. McCann and Falko (2006) Landau level degeneracy and quantum Hall effect in a graphite bilayer 10.1103/PhysRevLett.96.086805
  9. Giovannetti and Khomyakov (2007) Substrate-induced band gap in graphene on hexagonal boron nitride: Ab initio density functional calculations 10.1103/PhysRevB.76.073103
  10. Jung et al. (2015) Origin of band gaps in graphene on hexagonal boron nitride 10.1038/ncomms7308
  11. Zhou et al. (2007) Substrate-induced bandgap opening in epitaxial graphene 10.1038/nmat2003
  12. Nevius et al. (2015) Semiconducting Graphene from Highly Ordered Substrate Interactions 10.1103/PhysRevLett.115.136802
  13. Varchon et al. (2007) Electronic Structure of Epitaxial Graphene Layers on SiC: Effect of the Substrate 10.1103/PhysRevLett.99.126805
  14. Mattausch and Pankratov (2007) Ab Initio Study of Graphene on SiC 10.1103/PhysRevLett.99.076802
  15. Emtsev et al. (2008) Interaction, growth, and ordering of epitaxial graphene on SiC{0001} surfaces: A comparative photoelectron spectroscopy study 10.1103/PhysRevB.77.155303
  16. Varchon et al. (2008) Ripples in epitaxial graphene on the Si-terminated SiC(0001) surface 10.1103/PhysRevB.77.235412
  17. Varchon et al. (2008) Why Multilayer Graphene on 4H−SiC(000–1) Behaves Like a Single Sheet of Graphene 10.1103/PhysRevLett.100.125504
  18. Fratini and Guinea (2008) Substrate-limited electron dynamics in graphene 10.1103/PhysRevB.77.195415
  19. Zhao et al. (2015) One-step green synthesis of a ruthenium/graphene composite as a highly efficient catalyst 10.1039/C4RA11397G
  20. McCann and Abergel (2007) Electrons in bilayer graphene 10.1016/j.ssc.2007.03.054
  21. Zhang et al. (2009) Direct observation of a widely tunable bandgap in bilayer graphene 10.1038/nature08105
  22. Abergel (2015) Mucha-Kruczyński, Marcin.: Infrared absorption of closely aligned heterostructures of monolayer and bilayer graphene with hexagonal boron nitride 10.1103/PhysRevB.92.115430
  23. Ohta et al. (2006) Controlling the electronic structure of bilayer graphene 10.1126/science.1130681
  24. Hague (2011) Tunable graphene band gaps from superstrate-mediated interactions 10.1103/PhysRevB.84.155438
  25. Hague (2012) Polarons in highly doped atomically thin graphitic materials 10.1103/PhysRevB.86.064302
  26. Halasall and Ferrari (2009) Control of graphene's properties by reversible hydrogenation: evidence for graphane 10.1126/science.1167130
  27. Eklund et al. (2010) Reversible fluorination of graphene: Evidence of a two-dimensional wide bandgap semiconductor 10.1103/PhysRevB.81.205435
  28. Lin et al. (2013) Ab initio optical study of graphene on hexagonal boron nitride and fluorographene substrates 10.1039/c2tc00570k
  29. Song et al. (2010) Large scale growth and characterization of atomic hexagonal boron nitride layers 10.1021/nl1022139
  30. Wehling et al. (2011) Strength of Effective Coulomb Interactions in Graphene and Graphite 10.1103/PhysRevLett.106.236805
  31. Kotov and Uchoa (2012) Electron-Electron Interactions in Graphene: Current Status and Perspectives 10.1103/RevModPhys.84.1067
  32. Hansmann et al. (2013) Long-Range Coulomb Interactions in Surface Systems: A First-Principles Description within Self-Consistently Combined GW and Dynamical Mean-Field Theory 10.1103/PhysRevLett.110.166401
  33. Pariser and Parr (1953) A Semi-Empirical Theory of the Electronic Spectra and Electronic Structure of Complex Unsaturated Molecules 10.1063/1.1699030
  34. Soos et al. (1993) Band to correlated crossover in alternating Hubbard and Pariser-Parr-Pople chains: Nature of the lowest singlet excitation of conjugated polymers 10.1103/PhysRevLett.71.1609
  35. Friend et al. (1999) Electroluminescence in conjugated polymers 10.1038/16393
  36. Schuler et al. (2013) Optimal Hubbard Models for Materials with Nonlocal Coulomb Interactions: Graphene, Silicene, and Benzene 10.1103/PhysRevLett.111.036601
  37. Meng et al. (2010) Quantum spin liquid emerging in two-dimensional correlated Dirac fermions 10.1038/nature08942
  38. Osofsky et al. (2016) Functionalized graphene as a model system for the two-dimensional metal-insulator transition 10.1038/srep19939
  39. Igor (2006) Herbut.: Interactions and Phase Transitions on Graphene’s Honeycomb Lattice 10.1103/PhysRevLett.97.146401
  40. Drut and Lähde (2009) Is Graphene in Vacuum an Insulator? 10.1103/PhysRevLett.102.026802
  41. Castro Neto and Guinea (2009) The electronic properties of graphene 10.1103/RevModPhys.81.109
  42. Das Sarma and Shaffique (2011) Electronic transport in two-dimensional graphene 10.1103/RevModPhys.83.407
  43. Parr et al. (1950) Molecular orbstal calculations of the lower excited electronic levels of benzene, configuration interaction included 10.1063/1.1747540
  44. Wang et al. (2015) Visualizing chemical states and defects induced magnetism of graphene oxide by spatially-resolved-X-ray microscopy and spectroscopy 10.1038/srep15439
  45. Midtvedt and Alexander (2016) Strain–displacement relations for strain engineering in single-layer 2d materials 10.1088/0953-8984/28/4/045302
  46. Sahu and Rout (2014) Study of Band Gap Opening in Graphene by Impurity and Substrate-Mediated Interactions 10.1166/asl.2014.5374
  47. Sahu and Rout (2015) Model study of the effect of coulomb interaction on band gap of graphene-on-substrates 10.1016/j.physb.2014.12.014
  48. Sahu et al. (2016) Model study of band gap opening in graphene by electron–electron and electron–phonon interaction in high frequency range 10.1166/asl.2016.6894
  49. Sahu et al. (2016) Theoretical study of band gap opening in AB-stacked bi-layer graphene by impurity and electric field effects 10.1016/j.matpr.2016.01.119
  50. Zubarev (1960) Double-time green functions in statistical physics 10.1070/PU1960v003n03ABEH003275
  51. Sahu and Rout (2016) Tight-binding model study of substrate induced pseudo-spin polarization and magnetism in mono-layer graphene 10.1016/j.jmmm.2016.01.062
  52. Rout et al. (2016) Effect of impurity doping on tunnelling conductance in AB-stacked bi-layer graphene: A tight-binding study 10.1063/1.4945245
  53. Sahu et al. (2016) Microscopic theoretical model study of band gap opening in AA-stacked bi-layer graphene 10.1063/1.4946087