<h1>Band gap opening in graphene: a short theoretical study</h1>
			<ul class="item-list">
	        	        <li>
	            Sivabrata Sahu	            	            	            	            <sup aria-label="Affiliated with School of Applied Sciences (Physics), KIIT University, Bhubaneswar, Odisha, 751024, IN">
	                1	            </sup>
	            	        </li>
	        	        <li>
	            G. C. Rout	            	            <abbr title="This is the corresponding author for this article">*</abbr>
	            	            	                <a href="mailto:gcr@iopb.res.in" class="tiny-icon email-link mx-1" title="Email G. C. Rout">
	                    Email
	                </a>
	            	            	            <sup aria-label="Affiliated with Condensed Matter Physics Group, Physics Enclave, Bhubaneswar, Odisha, 751031, IN">
	                2	            </sup>
	            	        </li>
	        	    </ul>
	    	    <ol class="affiliations" aria-hidden="true">
	        <li>School of Applied Sciences (Physics), KIIT University, Bhubaneswar, Odisha, 751024, IN</li><li>Condensed Matter Physics Group, Physics Enclave, Bhubaneswar, Odisha, 751031, IN</li>	    </ol>
	    
<h2>Abstract</h2>
<p>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 
<span id="IEq1"><mml:math id="IEq1_Math"><mml:mrow><mml:msub><mml:mi>U</mml:mi><mml:mtext>C</mml:mtext></mml:msub><mml:mo>=</mml:mo><mml:mn>1.7</mml:mn><mml:msub><mml:mi>t</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:math><tex-math id="IEq1_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$U_{\text{C}} = 1.7 t_{1}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq1.gif"/></span>
. The critical Coulomb interaction is enhanced to 
<span id="IEq2"><mml:math id="IEq2_Math"><mml:mrow><mml:msub><mml:mi>U</mml:mi><mml:mtext>C</mml:mtext></mml:msub><mml:mo>=</mml:mo><mml:mn>2.5</mml:mn><mml:msub><mml:mi>t</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:math><tex-math id="IEq2_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$U_{\text{C}} = 2.5 t_{1}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq2.gif"/></span>
 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.</p><hr/><section><h2>Introduction</h2>
<p>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 [
<a href="#CR1"><sup>1</sup></a>
, 
<a href="#CR2"><sup>2</sup></a>
]. The tight-binding calculation for graphene shows that its conduction and valence bands touch at six Dirac points in the Brillouin zone [
<a href="#CR3"><sup>3</sup></a>
] 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 [
<a href="#CR4"><sup>4</sup></a>
–
<a href="#CR6"><sup>6</sup></a>
]. 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 [
<a href="#CR7"><sup>7</sup></a>
, 
<a href="#CR8"><sup>8</sup></a>
]. When graphene lies on born nitride (BN) substrate, a small gap (∼100 meV) is observed [
<a href="#CR9"><sup>9</sup></a>
, 
<a href="#CR10"><sup>10</sup></a>
]. Similarly, a band gap of 250 meV is observed for silicon carbide substrate [
<a href="#CR11"><sup>11</sup></a>
, 
<a href="#CR12"><sup>12</sup></a>
]. There is a great effort for producing a gap (∼1 eV) which is observed in germanium/silicon for the application in digital electronics.</p>
<p>The recent work indicates that the carbon layer is co-covalently bonded to the SiC sub-lattice [
<a href="#CR13"><sup>13</sup></a>
–
<a href="#CR17"><sup>17</sup></a>
]. The band gap opening for graphene-on-ruthenium metal also arises due to symmetry breaking in the system [
<a href="#CR18"><sup>18</sup></a>
, 
<a href="#CR19"><sup>19</sup></a>
]. McCann and Falko [
<a href="#CR8"><sup>8</sup></a>
, 
<a href="#CR20"><sup>20</sup></a>
] 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 [
<a href="#CR21"><sup>21</sup></a>
, 
<a href="#CR22"><sup>22</sup></a>
] and angle resolved photo emission spectroscopy (ARPES) [
<a href="#CR23"><sup>23</sup></a>
]. Hague [
<a href="#CR24"><sup>24</sup></a>
, 
<a href="#CR25"><sup>25</sup></a>
] 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 [
<a href="#CR26"><sup>26</sup></a>
] and Fluorographene [
<a href="#CR27"><sup>27</sup></a>
, 
<a href="#CR28"><sup>28</sup></a>
] by chemical modification of the system with hydrogen and fluorine, respectively. Atomic thick boron nitride (BN) forms a honeycomb lattice where the 
<em>π</em>
 orbitals on N sites are shifted up in energy by +∆ and decreased in energy of −∆ on B site causing a gap of 2∆ [
<a href="#CR29"><sup>29</sup></a>
]. The band gap of 5.56 eV is observed experimentally on monolayer BN systems.</p>
<p>The role of Coulomb interaction in low-dimensional systems provides unique opportunity for theoretical as well as experimental studies. The graphene [
<a href="#CR30"><sup>30</sup></a>
, 
<a href="#CR31"><sup>31</sup></a>
], varieties of semiconductor surfaces like Si, Ce, Sn, Pb [
<a href="#CR32"><sup>32</sup></a>
], Bechgard salt [
<a href="#CR33"><sup>33</sup></a>
] and doped polymers [
<a href="#CR34"><sup>34</sup></a>
, 
<a href="#CR35"><sup>35</sup></a>
] display strong on-site as well as inter-site Coulomb interactions. It is observed that on-site Coulomb interaction in graphene is 
<em>U</em>
 = 3.3
<em>t</em>
<sub>1</sub>
 and the nearest-neighbour Coulomb interaction is 
<span id="IEq3"><mml:math id="IEq3_Math"><mml:mrow><mml:mi>V</mml:mi><mml:mo>=</mml:mo><mml:mn>2.0</mml:mn><mml:msub><mml:mi>t</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:math><tex-math id="IEq3_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$V = 2.0 t_{1}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq3.gif"/></span>
 where the nearest-neighbour hopping integral is 
<em>t</em>
<sub>1</sub>
 = 2.8 eV [
<a href="#CR36"><sup>36</sup></a>
]. It has been reported that the effective Hubbard interaction is 
<span id="IEq4"><mml:math id="IEq4_Math"><mml:mrow><mml:mi>U</mml:mi><mml:mo>=</mml:mo><mml:mn>3.3</mml:mn><mml:msub><mml:mi>t</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:math><tex-math id="IEq4_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$U = 3.3 t_{1}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq4.gif"/></span>
 in the close vicinity of the separation between conducting graphene and insulating phases [
<a href="#CR37"><sup>37</sup></a>
, 
<a href="#CR38"><sup>38</sup></a>
]. 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 [
<a href="#CR39"><sup>39</sup></a>
, 
<a href="#CR40"><sup>40</sup></a>
]. The long ranged Coulomb interaction in pristine graphene leads to unusual behaviour [
<a href="#CR41"><sup>41</sup></a>
, 
<a href="#CR42"><sup>42</sup></a>
]. The on-site Coulomb interaction plays a crucial role for the understanding of the defect induced magnetisation [
<a href="#CR43"><sup>43</sup></a>
–
<a href="#CR45"><sup>45</sup></a>
].</p>
<p>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 [
<a href="#CR46"><sup>46</sup></a>
], in presence of Coulomb interaction [
<a href="#CR47"><sup>47</sup></a>
], electron phonon interaction [
<a href="#CR48"><sup>48</sup></a>
] and bi-layer graphene [
<a href="#CR49"><sup>49</sup></a>
]. 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.</p></section>
<section><h2>Tight-binding model</h2>
<p>The graphene sheet is formed by carbon atoms arranged in a non-Bravais honeycomb lattice with nearest-neighbour C–C distance of 
<span id="IEq5"><mml:math id="IEq5_Math"><mml:mrow><mml:msub><mml:mi>a</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mn>1.43</mml:mn></mml:mrow></mml:math><tex-math id="IEq5_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$a_{0} = 1.43$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq5.gif"/></span>
Å where the lattice constant is 
<span id="IEq6"><mml:math id="IEq6_Math"><mml:mrow><mml:mi>a</mml:mi><mml:mo>=</mml:mo><mml:msqrt><mml:mn>3</mml:mn></mml:msqrt><mml:msub><mml:mi>a</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow></mml:math><tex-math id="IEq6_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$a = \sqrt 3 a_{0}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq6.gif"/></span>
. The s, 
<span id="IEq7"><mml:math id="IEq7_Math"><mml:msub><mml:mi>p</mml:mi><mml:mi>x</mml:mi></mml:msub></mml:math><tex-math id="IEq7_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$p_{x}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq7.gif"/></span>
 and 
<span id="IEq8"><mml:math id="IEq8_Math"><mml:msub><mml:mi>p</mml:mi><mml:mi>y</mml:mi></mml:msub></mml:math><tex-math id="IEq8_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$p_{y}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq8.gif"/></span>
 orbitals hybridise to form 
<span id="IEq9"><mml:math id="IEq9_Math"><mml:mrow><mml:mi>s</mml:mi><mml:msup><mml:mi>p</mml:mi><mml:mn>2</mml:mn></mml:msup></mml:mrow></mml:math><tex-math id="IEq9_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$sp^{2}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq9.gif"/></span>
 bonds leading to high energy sigma bonds. The 
<span id="IEq10"><mml:math id="IEq10_Math"><mml:msub><mml:mi>p</mml:mi><mml:mi>z</mml:mi></mml:msub></mml:math><tex-math id="IEq10_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$p_{z}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq10.gif"/></span>
 orbitals in graphene form the 
<span id="IEq11"><mml:math id="IEq11_Math"><mml:mi mathvariant="italic">π</mml:mi></mml:math><tex-math id="IEq11_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\pi$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq11.gif"/></span>
 bond which is responsible for the electronic properties of graphene. The 2D character of the honeycomb lattice does not allow the overlap of the 
<span id="IEq12"><mml:math id="IEq12_Math"><mml:msub><mml:mi>p</mml:mi><mml:mi>z</mml:mi></mml:msub></mml:math><tex-math id="IEq12_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$p_{z}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq12.gif"/></span>
 orbital of carbon atom and s, 
<span id="IEq13"><mml:math id="IEq13_Math"><mml:msub><mml:mi>p</mml:mi><mml:mi>x</mml:mi></mml:msub></mml:math><tex-math id="IEq13_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$p_{x}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq13.gif"/></span>
 and 
<span id="IEq14"><mml:math id="IEq14_Math"><mml:msub><mml:mi>p</mml:mi><mml:mi>y</mml:mi></mml:msub></mml:math><tex-math id="IEq14_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$p_{y}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq14.gif"/></span>
 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.</p>
<p>Assuming the electron hopping up to the third nearest-neighbours, the tight-binding Hamiltonian for electron in graphene is written as
<section id="Equ1"><mml:math display="block" id="Equ1_Math"><mml:mrow><mml:msub><mml:mi>H</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:munder><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow></mml:munder><mml:mfenced close=")" open="(" separators=""><mml:msub><mml:mi mathvariant="italic">ϵ</mml:mi><mml:mi>a</mml:mi></mml:msub><mml:msubsup><mml:mi>a</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow><mml:mo>†</mml:mo></mml:msubsup><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msub><mml:mi mathvariant="italic">ϵ</mml:mi><mml:mi>b</mml:mi></mml:msub><mml:msubsup><mml:mi>b</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow><mml:mo>†</mml:mo></mml:msubsup><mml:msub><mml:mi>b</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow></mml:msub></mml:mfenced><mml:mo>-</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:munder><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:mo>&lt;</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mtext>j</mml:mtext><mml:mo>&gt;</mml:mo><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow></mml:munder><mml:mfenced close=")" open="(" separators=""><mml:mrow><mml:msubsup><mml:mi>a</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow><mml:mo>†</mml:mo></mml:msubsup><mml:msub><mml:mi>b</mml:mi><mml:mrow><mml:mi>j</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msubsup><mml:mi>b</mml:mi><mml:mrow><mml:mi>j</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow><mml:mo>†</mml:mo></mml:msubsup><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfenced><mml:mo>-</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:munder><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:mo>≪</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi><mml:mo>≫</mml:mo><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow></mml:munder><mml:mfenced close=")" open="(" separators=""><mml:mrow><mml:msubsup><mml:mi>a</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow><mml:mo>†</mml:mo></mml:msubsup><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>j</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msubsup><mml:mi>b</mml:mi><mml:mrow><mml:mi>j</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow><mml:mo>†</mml:mo></mml:msubsup><mml:msub><mml:mi>b</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfenced><mml:mo>-</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mn>3</mml:mn></mml:msub><mml:munder><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:mo>≪</mml:mo><mml:mo>&lt;</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mtext>j</mml:mtext><mml:mo>≫</mml:mo><mml:mo>&gt;</mml:mo><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow></mml:munder><mml:mfenced close=")" open="(" separators=""><mml:mrow><mml:msubsup><mml:mi>a</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow><mml:mo>†</mml:mo></mml:msubsup><mml:msub><mml:mi>b</mml:mi><mml:mrow><mml:mi>j</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:msubsup><mml:mi>b</mml:mi><mml:mrow><mml:mi>j</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow><mml:mo>†</mml:mo></mml:msubsup><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mfenced></mml:mrow></mml:math><tex-math id="Equ1_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$H_{0} = \mathop \sum \limits_{i,\sigma } \left( \epsilon_{a} a_{i,\sigma }^{\dag } a_{i,\sigma } +\epsilon_{b } b_{i,\sigma }^{\dag } b_{i,\sigma } \right) - t_{1} \mathop \sum \limits_{{ < i, {\text{j}} >,\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}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_Equ1.gif"/></section>
where 
<span id="IEq15"><mml:math id="IEq15_Math"><mml:mrow><mml:msubsup><mml:mi>a</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow><mml:mo>†</mml:mo></mml:msubsup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:math><tex-math id="IEq15_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$a_{i,\sigma }^{\dag } (a_{i,\sigma } )$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq15.gif"/></span>
 creates (annihilates) an electron with spin 
<em>σ</em>
 (
<em>σ</em>
 = ↑,↓) on site 
<span id="IEq16"><mml:math id="IEq16_Math"><mml:msub><mml:mover accent="true"><mml:mtext>R</mml:mtext><mml:mo stretchy="false">→</mml:mo></mml:mover><mml:mtext>i</mml:mtext></mml:msub></mml:math><tex-math id="IEq16_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$${\vec{\text{R}}}_{\text{i}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq16.gif"/></span>
 on sub-lattice A. Similarly, 
<span id="IEq17"><mml:math id="IEq17_Math"><mml:mrow><mml:msubsup><mml:mi>b</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow><mml:mo>†</mml:mo></mml:msubsup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>b</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:math><tex-math id="IEq17_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$b_{i,\sigma }^{\dag } (b_{i,\sigma } )$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq17.gif"/></span>
 creates (annihilates) an electron on sub-lattice B. Here, 
<span id="IEq18"><mml:math id="IEq18_Math"><mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mfenced close=")" open="(" separators=""><mml:mrow><mml:mo>=</mml:mo><mml:mn>2.5</mml:mn><mml:mrow><mml:mspace width="0.277778em"/><mml:mrow><mml:mi>t</mml:mi><mml:mi>o</mml:mi></mml:mrow></mml:mrow><mml:mspace width="0.277778em"/><mml:mn>3.0</mml:mn><mml:mtext>eV</mml:mtext></mml:mrow></mml:mfenced></mml:mrow></mml:math><tex-math id="IEq18_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$t_{1} \left( { = 2.5 {\text\;{to }}\;3.0 {\text{eV}}} \right)$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq18.gif"/></span>
 is the nearest-neighbour hopping energy, 
<em>t</em>
<sub>2</sub>
 with 
<span id="IEq19"><mml:math id="IEq19_Math"><mml:mrow><mml:mn>0.02</mml:mn><mml:msub><mml:mi>t</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mo>≤</mml:mo><mml:msub><mml:mi>t</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mo>≤</mml:mo><mml:mn>0.2</mml:mn><mml:msub><mml:mi>t</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:math><tex-math id="IEq19_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$0.02t_{1} \le t_{2} \le 0.2 t_{1}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq19.gif"/></span>
 [
<a href="#CR36"><sup>36</sup></a>
, 
<a href="#CR37"><sup>37</sup></a>
] is the next-nearest-neighbour hopping energy, 
<span id="IEq20"><mml:math id="IEq20_Math"><mml:msub><mml:mi>t</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:math><tex-math id="IEq20_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$t_{3}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq20.gif"/></span>
 is the third nearest-neighbour hopping integral with 
<span id="IEq21"><mml:math id="IEq21_Math"><mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mn>3</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mn>0.024</mml:mn><mml:mspace width="0.166667em"/><mml:msub><mml:mi>t</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mspace width="0.277778em"/><mml:mtext>eV</mml:mtext></mml:mrow></mml:math><tex-math id="IEq21_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$t_{3} = 0.024\, t_{1} \;{\text{eV}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq21.gif"/></span>
 and 
<span id="IEq22"><mml:math id="IEq22_Math"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mtext>a</mml:mtext></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mtext>b</mml:mtext><mml:mspace width="0.333333em"/></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:math><tex-math id="IEq22_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\varepsilon_{\text{a}} (\varepsilon_{\text{b }} )$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq22.gif"/></span>
 is the site energy at the sub-lattice site A (B). Further, 
<span id="IEq23"><mml:math id="IEq23_Math"><mml:mrow><mml:mo>&lt;</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi><mml:mo>&gt;</mml:mo><mml:mo>,</mml:mo></mml:mrow></mml:math><tex-math id="IEq23_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$< i, j > ,$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq23.gif"/></span>
 ≪ 
<em>i</em>
, 
<em>j</em>
 ≫,  and ≪ &lt; 
<em>i</em>
, 
<em>j</em>
 &gt; ≫  stand for nearest, next-nearest and next–next-neighbour hoppings from site 
<span id="IEq24"><mml:math id="IEq24_Math"><mml:msub><mml:mover accent="true"><mml:mtext>R</mml:mtext><mml:mo stretchy="false">→</mml:mo></mml:mover><mml:mi>i</mml:mi></mml:msub></mml:math><tex-math id="IEq24_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$${\vec{\text{R}}}_{i}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq24.gif"/></span>
 to 
<span id="IEq25"><mml:math id="IEq25_Math"><mml:msub><mml:mover accent="true"><mml:mtext>R</mml:mtext><mml:mo stretchy="false">→</mml:mo></mml:mover><mml:mi>j</mml:mi></mml:msub></mml:math><tex-math id="IEq25_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$${\vec{\text{R}}}_{j}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq25.gif"/></span>
.</p>
<p>The Fourier transformed dispersion 
<em>γ</em>
<sub>1</sub>
(
<em>k</em>
)for the nearest-neighbour hopping is
<section id="Equ2"><mml:math display="block" id="Equ2_Math"><mml:mrow><mml:msub><mml:mi mathvariant="italic">γ</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mfenced close=")" open="("><mml:mi>k</mml:mi></mml:mfenced><mml:mo>=</mml:mo><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:msub><mml:mi>k</mml:mi><mml:mi>x</mml:mi></mml:msub><mml:msub><mml:mi>a</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow></mml:msup><mml:mo>+</mml:mo><mml:mn>2</mml:mn><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mfrac><mml:mi>i</mml:mi><mml:mn>2</mml:mn></mml:mfrac><mml:msub><mml:mi>k</mml:mi><mml:mi>x</mml:mi></mml:msub><mml:msub><mml:mi>a</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow></mml:msup><mml:mo>·</mml:mo><mml:mo>cos</mml:mo><mml:mfrac><mml:msqrt><mml:mn>3</mml:mn></mml:msqrt><mml:mn>2</mml:mn></mml:mfrac><mml:msub><mml:mi>k</mml:mi><mml:mi>y</mml:mi></mml:msub><mml:msub><mml:mi>a</mml:mi><mml:mn>0</mml:mn></mml:msub></mml:mrow></mml:math><tex-math id="Equ2_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\gamma_{1} \left( k \right) = e^{{ik_{x} a_{0} }} + 2 e^{{ - \frac{i}{2}k_{x} a_{0} }} \cdot { \cos }\frac{\sqrt 3 }{2}k_{y} a_{0}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_Equ2.gif"/></section>
and dispersions for the next-and next-to-next-nearest-neighbour hoppings are 
<span id="IEq26"><mml:math id="IEq26_Math"><mml:mrow><mml:msub><mml:mi mathvariant="italic">γ</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mfenced close=")" open="("><mml:mi>k</mml:mi></mml:mfenced><mml:mo>=</mml:mo><mml:msub><mml:mo>∑</mml:mo><mml:msub><mml:mi mathvariant="italic">δ</mml:mi><mml:mn>2</mml:mn></mml:msub></mml:msub><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mover accent="true"><mml:mi>k</mml:mi><mml:mo stretchy="false">→</mml:mo></mml:mover><mml:mo>.</mml:mo><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">δ</mml:mi><mml:mo stretchy="false">→</mml:mo></mml:mover><mml:mn>2</mml:mn></mml:msub></mml:mrow></mml:msup></mml:mrow></mml:math><tex-math id="IEq26_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\gamma_{2} \left( k \right) = \sum\nolimits_{{\delta_{2} }} {e^{{i\vec{k}.\vec{\delta }_{2} }} }$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq26.gif"/></span>
 and 
<span id="IEq27"><mml:math id="IEq27_Math"><mml:mrow><mml:msub><mml:mi mathvariant="italic">γ</mml:mi><mml:mn>3</mml:mn></mml:msub><mml:mfenced close=")" open="("><mml:mi>k</mml:mi></mml:mfenced><mml:mo>=</mml:mo><mml:munder><mml:mo movablelimits="false">∑</mml:mo><mml:msub><mml:mi mathvariant="italic">δ</mml:mi><mml:mn>3</mml:mn></mml:msub></mml:munder><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mover accent="true"><mml:mi>k</mml:mi><mml:mo stretchy="false">→</mml:mo></mml:mover><mml:mo>.</mml:mo><mml:msub><mml:mover accent="true"><mml:mi mathvariant="italic">δ</mml:mi><mml:mo stretchy="false">→</mml:mo></mml:mover><mml:mn>3</mml:mn></mml:msub></mml:mrow></mml:msup></mml:mrow></mml:math><tex-math id="IEq27_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\gamma_{3} \left( k \right) = \mathop \sum \limits_{{\delta_{3} }} e^{{i\vec{k}.\vec{\delta }_{3} }}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq27.gif"/></span>
. where 
<span id="IEq28"><mml:math id="IEq28_Math"><mml:mover accent="true"><mml:msub><mml:mi mathvariant="italic">δ</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mo stretchy="false">→</mml:mo></mml:mover></mml:math><tex-math id="IEq28_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\overrightarrow {{\delta_{2} }}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq28.gif"/></span>
 and 
<span id="IEq29"><mml:math id="IEq29_Math"><mml:mover accent="true"><mml:msub><mml:mi mathvariant="italic">δ</mml:mi><mml:mn>3</mml:mn></mml:msub><mml:mo stretchy="false">→</mml:mo></mml:mover></mml:math><tex-math id="IEq29_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\overrightarrow {{\delta_{3} }}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq29.gif"/></span>
 are the corresponding lattice vectors. Graphene deposited on silicon dioxide (SiO
<sub>2</sub>
) substrate can be described by the 2D massless Dirac fermions [
<a href="#CR3"><sup>3</sup></a>
]. 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 
<span id="IEq30"><mml:math id="IEq30_Math"><mml:mrow><mml:mo>+</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi></mml:mrow></mml:math><tex-math id="IEq30_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$+ \Delta$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq30.gif"/></span>
 and the B site with energy 
<span id="IEq31"><mml:math id="IEq31_Math"><mml:mrow><mml:mo>-</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi></mml:mrow></mml:math><tex-math id="IEq31_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$- \Delta$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq31.gif"/></span>
 leading to the symmetry breaking between A and B sites. Such a symmetry breaking Hamiltonian is written as
<section id="Equ3"><mml:math display="block" id="Equ3_Math"><mml:mrow><mml:msub><mml:mi>H</mml:mi><mml:mtext>sub</mml:mtext></mml:msub><mml:mo>=</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:munder><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow></mml:munder><mml:msubsup><mml:mi>a</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow><mml:mo>†</mml:mo></mml:msubsup><mml:msub><mml:mi>a</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:munder><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow></mml:munder><mml:msubsup><mml:mi>b</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow><mml:mo>†</mml:mo></mml:msubsup><mml:msub><mml:mi>b</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow></mml:msub><mml:mo>.</mml:mo></mml:mrow></mml:math><tex-math id="Equ3_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$H_{\text{sub}} = \Delta \mathop \sum \limits_{i,\sigma } a_{i,\sigma }^{\dag } a_{i,\sigma } - \Delta \mathop \sum \limits_{i,\sigma } b_{i,\sigma }^{\dag } b_{i,\sigma } .$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_Equ3.gif"/></section></p>
<p>The system exhibits an insulating ground state with a band gap 
<span id="IEq32"><mml:math id="IEq32_Math"><mml:mrow><mml:mn>2</mml:mn><mml:mi mathvariant="normal">Δ</mml:mi></mml:mrow></mml:math><tex-math id="IEq32_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$2\Delta$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq32.gif"/></span>
, 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
<section id="Equ4"><mml:math display="block" id="Equ4_Math"><mml:mrow><mml:msub><mml:mi>H</mml:mi><mml:mtext>U</mml:mtext></mml:msub><mml:mo>=</mml:mo><mml:mi>U</mml:mi><mml:munder><mml:mo movablelimits="false">∑</mml:mo><mml:mi>i</mml:mi></mml:munder><mml:mfenced close="]" open="[" separators=""><mml:mrow><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mo stretchy="false">↑</mml:mo></mml:mrow><mml:mi>a</mml:mi></mml:msubsup><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mo stretchy="false">↓</mml:mo></mml:mrow><mml:mi>a</mml:mi></mml:msubsup><mml:mo>+</mml:mo><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mo stretchy="false">↑</mml:mo></mml:mrow><mml:mi>b</mml:mi></mml:msubsup><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mo stretchy="false">↓</mml:mo></mml:mrow><mml:mi>b</mml:mi></mml:msubsup></mml:mrow></mml:mfenced></mml:mrow></mml:math><tex-math id="Equ4_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$H_{\text{U}} = U\mathop \sum \limits_{i} \left[ {n_{i, \uparrow }^{a} n_{i, \downarrow }^{a} + n_{i, \uparrow }^{b} n_{i, \downarrow }^{b} } \right]$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_Equ4.gif"/></section>
where 
<em>n</em>
<sub><em>i</em>
,↑</sub>
<sup><em>α</em></sup>
(
<em>n</em>
<sub><em>i</em>
,↓</sub>
<sup><em>α</em></sup>
) with 
<em>α</em>
 ∊ 
<em>A</em>
, 
<em>B</em>
 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. 
<span id="IEq33"><mml:math id="IEq33_Math"><mml:mrow><mml:mi>U</mml:mi><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mo stretchy="false">↑</mml:mo></mml:mrow><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mo stretchy="false">↓</mml:mo></mml:mrow><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mo>≈</mml:mo><mml:mi>U</mml:mi><mml:mo>&lt;</mml:mo><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mo stretchy="false">↑</mml:mo></mml:mrow><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mo>&gt;</mml:mo><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mo stretchy="false">↓</mml:mo></mml:mrow><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mo>+</mml:mo><mml:mi>U</mml:mi><mml:mo>&lt;</mml:mo><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mo stretchy="false">↓</mml:mo></mml:mrow><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mo>&gt;</mml:mo><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mo stretchy="false">↑</mml:mo></mml:mrow><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mo>-</mml:mo><mml:mi>U</mml:mi><mml:mo>&lt;</mml:mo><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mo stretchy="false">↑</mml:mo></mml:mrow><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mo>&gt;</mml:mo><mml:mo>&lt;</mml:mo><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mo stretchy="false">↓</mml:mo></mml:mrow><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mo>&gt;</mml:mo></mml:mrow></mml:math><tex-math id="IEq33_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$Un_{i, \uparrow }^{\alpha } n_{i, \downarrow }^{\alpha } \approx U < n_{i, \uparrow }^{\alpha } > n_{i, \downarrow }^{\alpha } + U < n_{i, \downarrow }^{\alpha } > n_{i, \uparrow }^{\alpha } - U < n_{i, \uparrow }^{\alpha } > < n_{i, \downarrow }^{\alpha } >$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq33.gif"/></span>
, where 
<em>α</em>
 ≡ 
<em>a</em>
, 
<em>b</em>
 corresponding to A and B site interactions. The mean-field solutions are taken as 
<span id="IEq34"><mml:math id="IEq34_Math"><mml:mrow><mml:mfrac><mml:mrow><mml:mo>&lt;</mml:mo><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mo stretchy="false">↑</mml:mo></mml:mrow><mml:mi>a</mml:mi></mml:msubsup><mml:mo>&gt;</mml:mo><mml:mo>+</mml:mo><mml:mo>&lt;</mml:mo><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo stretchy="false">↑</mml:mo></mml:mrow><mml:mi>b</mml:mi></mml:msubsup><mml:mo>&gt;</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:mfrac><mml:mo>=</mml:mo><mml:mi>n</mml:mi></mml:mrow></mml:math><tex-math id="IEq34_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\frac{{ < n_{i, \uparrow }^{a} > + < n_{i \uparrow }^{b} > }}{2} = n$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq34.gif"/></span>
 and 
<span id="IEq35"><mml:math id="IEq35_Math"><mml:mrow><mml:mfrac><mml:mrow><mml:mo>&lt;</mml:mo><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mo stretchy="false">↑</mml:mo></mml:mrow><mml:mi>a</mml:mi></mml:msubsup><mml:mo>&gt;</mml:mo><mml:mo>-</mml:mo><mml:mo>&lt;</mml:mo><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo stretchy="false">↑</mml:mo></mml:mrow><mml:mi>b</mml:mi></mml:msubsup><mml:mo>&gt;</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:mfrac><mml:mo>=</mml:mo><mml:mi>d</mml:mi></mml:mrow></mml:math><tex-math id="IEq35_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\frac{{ < n_{i, \uparrow }^{a} > - < n_{i \uparrow }^{b} > }}{2} = d$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq35.gif"/></span>
 and this leads to the condition,
<span id="IEq36"><mml:math id="IEq36_Math"><mml:mrow><mml:mo>&lt;</mml:mo><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mo stretchy="false">↑</mml:mo></mml:mrow><mml:mi>a</mml:mi></mml:msubsup><mml:mo>&gt;</mml:mo><mml:mo>=</mml:mo><mml:mi>n</mml:mi><mml:mo>+</mml:mo><mml:mi>d</mml:mi></mml:mrow></mml:math><tex-math id="IEq36_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$< n_{i, \uparrow }^{a} > = n + d$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq36.gif"/></span>
 and 
<span id="IEq37"><mml:math id="IEq37_Math"><mml:mrow><mml:mo>&lt;</mml:mo><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mo stretchy="false">↑</mml:mo></mml:mrow><mml:mi>b</mml:mi></mml:msubsup><mml:mo>&gt;</mml:mo><mml:mo>=</mml:mo><mml:mi>n</mml:mi><mml:mo>-</mml:mo><mml:mi>d</mml:mi></mml:mrow></mml:math><tex-math id="IEq37_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$< n_{i, \uparrow }^{b} > = n - d$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq37.gif"/></span>
 where 
<em>n</em>
 represents the mean electron occupation and 
<em>d</em>
, 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
<section id="Equ5"><mml:math display="block" id="Equ5_Math"><mml:mrow><mml:msub><mml:mi>H</mml:mi><mml:mrow><mml:mtext>e - ph</mml:mtext></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:munder><mml:mo movablelimits="false">∑</mml:mo><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mo>,</mml:mo><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi></mml:mrow></mml:munder><mml:mi>f</mml:mi><mml:mfenced close=")" open="(" separators=""><mml:mrow><mml:mi>i</mml:mi><mml:mo>-</mml:mo><mml:mi>j</mml:mi></mml:mrow></mml:mfenced><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow><mml:mi mathvariant="italic">α</mml:mi></mml:msubsup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msubsup><mml:mi>b</mml:mi><mml:mrow><mml:mi>j</mml:mi></mml:mrow><mml:mo>†</mml:mo></mml:msubsup><mml:mo>+</mml:mo><mml:msub><mml:mi>b</mml:mi><mml:mi>j</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:math><tex-math id="Equ5_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$H_{\text{e - ph}} = \mathop \sum \limits_{\alpha ,i,j} f\left( {i - j} \right)n_{i,\sigma }^{\alpha } (b_{j}^{\dag } + b_{j} )$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_Equ5.gif"/></section>
<section id="Equ6"><mml:math display="block" id="Equ6_Math"><mml:mrow><mml:msub><mml:mi>H</mml:mi><mml:mtext>p</mml:mtext></mml:msub><mml:mo>=</mml:mo><mml:munder><mml:mo movablelimits="false">∑</mml:mo><mml:mi>j</mml:mi></mml:munder><mml:msub><mml:mi mathvariant="italic">ω</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mspace width="0.277778em"/><mml:msubsup><mml:mi>b</mml:mi><mml:mrow><mml:mi>j</mml:mi></mml:mrow><mml:mo>†</mml:mo></mml:msubsup><mml:msub><mml:mi>b</mml:mi><mml:mi>j</mml:mi></mml:msub><mml:mo>.</mml:mo></mml:mrow></mml:math><tex-math id="Equ6_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$H_{\text{p}} = \mathop \sum \limits_{j} \omega_{0}\; b_{j}^{\dag } b_{j} .$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_Equ6.gif"/></section>
Here 
<em>H</em>
<sub>e-ph</sub>
 describes electron–phonon interaction where the electron density 
<em>n</em>
<sub><em>i</em>
,
<em>σ</em></sub>
<sup><em>α</em></sup>
 at sub-lattice site is coupled to the phonon displacement 
<span id="IEq38"><mml:math id="IEq38_Math"><mml:mrow><mml:msub><mml:mi>u</mml:mi><mml:mi>j</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:msubsup><mml:mi>b</mml:mi><mml:mrow><mml:mi>j</mml:mi></mml:mrow><mml:mo>†</mml:mo></mml:msubsup><mml:mo>+</mml:mo><mml:msub><mml:mi>b</mml:mi><mml:mi>j</mml:mi></mml:msub></mml:mrow></mml:math><tex-math id="IEq38_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$u_{j} = b_{j}^{\dag } + b_{j}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq38.gif"/></span>
 with phonon creation (annihilation) operator 
<span id="IEq39"><mml:math id="IEq39_Math"><mml:mrow><mml:msubsup><mml:mi>b</mml:mi><mml:mrow><mml:mi>j</mml:mi></mml:mrow><mml:mo>†</mml:mo></mml:msubsup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi>b</mml:mi><mml:mi>j</mml:mi></mml:msub><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:math><tex-math id="IEq39_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$b_{j}^{\dag } (b_{j} )$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq39.gif"/></span>
 at site 
<em>j</em>
 and 
<em>f</em>
(
<em>i</em>
−
<em>j</em>
) is the electron–phonon coupling. For high frequency phonons, Eq. (
<a href="#Equ5"><sup>5</sup></a>
) represents the Holstein interaction. The Hamiltonian 
<em>H</em>
<sub>p</sub>
 represents the free phonon energy with phonon frequency 
<em>ω</em>
<sub>0</sub>
.</p></section>
<section><h2>Calculation of Green’s functions and quasi-particle bands</h2>
<p>The Green’s functions for the electrons of A and B sub-lattices are calculated by Zubarev’s Green’s function technique [
<a href="#CR50"><sup>50</sup></a>
]. The Zubarev’s double time single particle retarded Green’s function is defined
<section id="Equ7"><mml:math display="block" id="Equ7_Math"><mml:mrow><mml:msub><mml:mi>G</mml:mi><mml:mtext>r</mml:mtext></mml:msub><mml:mfenced close=")" open="(" separators=""><mml:mrow><mml:mi>t</mml:mi><mml:mo>,</mml:mo><mml:msup><mml:mi>t</mml:mi><mml:msup><mml:mrow/><mml:mo>′</mml:mo></mml:msup></mml:msup></mml:mrow></mml:mfenced><mml:mo>=</mml:mo><mml:mo>≪</mml:mo><mml:mi>A</mml:mi><mml:mfenced close=")" open="("><mml:mi>T</mml:mi></mml:mfenced><mml:mo>;</mml:mo><mml:mi>B</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msup><mml:mi>t</mml:mi><mml:msup><mml:mrow/><mml:mo>′</mml:mo></mml:msup></mml:msup><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:msub><mml:mo>≫</mml:mo><mml:mtext>r</mml:mtext></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mi>i</mml:mi><mml:mi mathvariant="italic">Θ</mml:mi><mml:mfenced close=")" open="(" separators=""><mml:mrow><mml:mi mathvariant="normal">t</mml:mi><mml:mo>-</mml:mo><mml:msup><mml:mi mathvariant="normal">t</mml:mi><mml:msup><mml:mrow/><mml:mo>′</mml:mo></mml:msup></mml:msup></mml:mrow></mml:mfenced><mml:mo>&lt;</mml:mo><mml:mfenced close="]" open="[" separators=""><mml:mrow><mml:mi mathvariant="normal">A</mml:mi><mml:mfenced close=")" open="("><mml:mi mathvariant="normal">t</mml:mi></mml:mfenced><mml:mo>;</mml:mo><mml:mi mathvariant="normal">B</mml:mi><mml:mfenced close=")" open="(" separators=""><mml:msup><mml:mi mathvariant="normal">t</mml:mi><mml:msup><mml:mrow/><mml:mo>′</mml:mo></mml:msup></mml:msup></mml:mfenced></mml:mrow></mml:mfenced><mml:mo>&gt;</mml:mo></mml:mrow></mml:math><tex-math id="Equ7_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$ G_{\text{r}} \left( {t,t^{'} } \right) = \ll A\left( T \right);B(t^{'} ) \gg_{\text{r}} = - i\rm\varTheta \left( {t - t^{'} } \right) < \left[ {A\left( t \right);B\left( {t^{'} } \right)} \right] > $$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_Equ7.gif"/></section>
where &lt;…&gt; indicates the average over a grand canonical ensemble. 
<span id="IEq40"><mml:math id="IEq40_Math"><mml:mrow><mml:mi>A</mml:mi><mml:mfenced close=")" open="("><mml:mi>t</mml:mi></mml:mfenced><mml:mspace width="0.277778em"/><mml:mtext>and</mml:mtext><mml:mspace width="0.277778em"/><mml:mi>B</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msup><mml:mi>t</mml:mi><mml:msup><mml:mrow/><mml:mo>′</mml:mo></mml:msup></mml:msup><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:math><tex-math id="IEq40_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$A\left( t \right)\; {\text{and}}\;B(t^{'} )$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq40.gif"/></span>
 are Heisenberg representations of the operators, while Θ = 1 for 
<em>t</em>
 &gt; 
<em>t</em>
<sup>′</sup>
 and Θ = 0 for 
<em>t</em>
 &lt; 
<em>t</em>
<sup>′</sup>
. Here, [
<em>A</em>
(
<em>t</em>
); 
<em>B</em>
(
<em>t</em>
<sup>′</sup>
)] indicates commutation or anti-commutation depending on Boson or fermions operators. The poles of the Green’s functions provide quasi-particle energies.</p>
<p>The correlation function corresponding to Green’s function 
<em>G</em>
<sub>r</sub>
(
<em>t</em>
, 
<em>t</em>
<sup>′</sup>
) is defined as
<section id="Equ8"><mml:math display="block" id="Equ8_Math"><mml:mrow><mml:mo>&lt;</mml:mo><mml:mi>B</mml:mi><mml:mfenced close=")" open="(" separators=""><mml:msup><mml:mi>t</mml:mi><mml:msup><mml:mrow/><mml:mo>′</mml:mo></mml:msup></mml:msup></mml:mfenced><mml:mi>A</mml:mi><mml:mfenced close=")" open="("><mml:mi>t</mml:mi></mml:mfenced><mml:mo>&gt;</mml:mo><mml:mo>=</mml:mo><mml:munderover><mml:mo movablelimits="false">∫</mml:mo><mml:mrow><mml:mo>-</mml:mo><mml:mi>∞</mml:mi></mml:mrow><mml:mi>∞</mml:mi></mml:munderover><mml:msub><mml:mi>J</mml:mi><mml:mi>k</mml:mi></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:msup><mml:mi>e</mml:mi><mml:mrow><mml:mo>-</mml:mo><mml:mi>i</mml:mi><mml:mi mathvariant="italic">ω</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo>-</mml:mo><mml:msup><mml:mi>t</mml:mi><mml:msup><mml:mrow/><mml:mo>′</mml:mo></mml:msup></mml:msup><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:msup><mml:mtext>d</mml:mtext><mml:mi mathvariant="italic">ω</mml:mi></mml:mrow></mml:math><tex-math id="Equ8_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$< B\left( {t^{'} } \right)A\left( t \right) > = \mathop \int \limits_{ - \infty }^{\infty } J_{k} (\omega ) e^{{ - i\omega (t - t^{'} )}} {\text{d}}\omega$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_Equ8.gif"/></section>
where the spectral density function 
<em>J</em>
<sub><em>k</em></sub>
(
<em>ω</em>
) is written as
<section id="Equ9"><mml:math display="block" id="Equ9_Math"><mml:mrow><mml:msub><mml:mi>J</mml:mi><mml:mi>k</mml:mi></mml:msub><mml:mfenced close=")" open="("><mml:mi mathvariant="italic">ω</mml:mi></mml:mfenced><mml:mo>=</mml:mo><mml:mi>i</mml:mi><mml:munder><mml:mo movablelimits="false">lim</mml:mo><mml:mrow><mml:mi mathvariant="italic">η</mml:mi><mml:mo stretchy="false">→</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:munder><mml:mfrac><mml:mn>1</mml:mn><mml:mrow><mml:msup><mml:mi>e</mml:mi><mml:mfrac><mml:mi mathvariant="italic">ω</mml:mi><mml:mi mathvariant="italic">θ</mml:mi></mml:mfrac></mml:msup><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:mfrac><mml:mfenced close="]" open="[" separators=""><mml:mrow><mml:msub><mml:mi>G</mml:mi><mml:mtext>r</mml:mtext></mml:msub><mml:mfenced close=")" open="(" separators=""><mml:mrow><mml:mi mathvariant="italic">ω</mml:mi><mml:mo>+</mml:mo><mml:mi>i</mml:mi><mml:mi mathvariant="italic">η</mml:mi></mml:mrow></mml:mfenced><mml:mo>-</mml:mo><mml:msub><mml:mi>G</mml:mi><mml:mtext>r</mml:mtext></mml:msub><mml:mfenced close=")" open="(" separators=""><mml:mrow><mml:mi mathvariant="italic">ω</mml:mi><mml:mo>-</mml:mo><mml:mi>i</mml:mi><mml:mi mathvariant="italic">η</mml:mi></mml:mrow></mml:mfenced></mml:mrow></mml:mfenced></mml:mrow></mml:math><tex-math id="Equ9_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$J_{k} \left( \omega \right) = i \mathop {\lim }\limits_{\eta \to 0} \frac{1}{{e^{{\frac{\omega }{\theta }}} + 1}}\left[ { G_{\text{r}} \left( {\omega + i\eta } \right) - G_{\text{r}} \left( {\omega - i\eta } \right)} \right]$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_Equ9.gif"/></section>
where 
<em>θ</em>
 = 
<em>k</em>
<sub>B</sub>
<em>T</em>
 with 
<em>η</em>
 as a small spectral width. Other physical parameters can be calculated using the correlation functions of the corresponding Green’s functions.</p>
<p>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
<section id="Equ10"><mml:math display="block" id="Equ10_Math"><mml:mrow><mml:mtext>DOS</mml:mtext><mml:mo>=</mml:mo><mml:mfrac><mml:mi>A</mml:mi><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mn>2</mml:mn><mml:mi mathvariant="italic">π</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:msup></mml:mfrac><mml:mrow><mml:mspace width="0.277778em"/><mml:mpadded width="-3pt"><mml:mo>∫</mml:mo></mml:mpadded><mml:mpadded width="-3pt"><mml:mo>∫</mml:mo></mml:mpadded><mml:mspace width="0.277778em"/></mml:mrow><mml:msub><mml:mover accent="true"><mml:mrow><mml:mtext>d</mml:mtext><mml:mi>k</mml:mi></mml:mrow><mml:mo stretchy="false">→</mml:mo></mml:mover><mml:mi>x</mml:mi></mml:msub><mml:msub><mml:mover accent="true"><mml:mrow><mml:mtext>d</mml:mtext><mml:mi>k</mml:mi></mml:mrow><mml:mo stretchy="false">→</mml:mo></mml:mover><mml:mi>y</mml:mi></mml:msub><mml:msub><mml:mi mathvariant="italic">ρ</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi mathvariant="italic">ω</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:math><tex-math id="Equ10_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$${\text{DOS}} = \frac{A}{{(2\pi )^{2} }}\mathop {\iint }\nolimits \overrightarrow {{{\text{d}}k}}_{x} \overrightarrow {{{\text{d}}k}}_{y} \rho_{k,\sigma } (\omega )$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_Equ10.gif"/></section>
where 
<em>A</em>
 is the area of the honeycomb lattice in real space and spectral function 
<em>ρ</em>
<sub><em>k</em>
,
<em>σ</em></sub>
(
<em>ω</em>
) is given by
<section id="Equ11"><mml:math display="block" id="Equ11_Math"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ρ</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:mrow></mml:msub><mml:mfenced close=")" open="("><mml:mi mathvariant="italic">ω</mml:mi></mml:mfenced><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>2</mml:mn><mml:mi mathvariant="italic">π</mml:mi><mml:munder><mml:mo movablelimits="false">∑</mml:mo><mml:mi mathvariant="italic">σ</mml:mi></mml:munder><mml:mfenced close="]" open="[" separators=""><mml:mrow><mml:mtext>Im</mml:mtext><mml:mfenced close=")" open="(" separators=""><mml:mrow><mml:msub><mml:mi>G</mml:mi><mml:mtext>r</mml:mtext></mml:msub><mml:mfenced close=")" open="(" separators=""><mml:mrow><mml:mi mathvariant="italic">ω</mml:mi><mml:mo>+</mml:mo><mml:mi>i</mml:mi><mml:mi mathvariant="italic">η</mml:mi></mml:mrow></mml:mfenced></mml:mrow></mml:mfenced></mml:mrow></mml:mfenced></mml:mrow></mml:math><tex-math id="Equ11_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\rho_{k,\sigma } \left( \omega \right) = - 2\pi \mathop \sum \limits_{\sigma } \left[ {{\text{Im}}\left( {G_{\text{r}} \left( {\omega + i\eta } \right)} \right)} \right]$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_Equ11.gif"/></section>
with 
<em>η</em>
 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 
<span id="IEq41"><mml:math id="IEq41_Math"><mml:msub><mml:mi>t</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:math><tex-math id="IEq41_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$t_{1}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq41.gif"/></span>
.</p></section>
<section><h2>Results and discussion</h2>
<p>The band gaps can be induced in graphene near Dirac point by the following techniques.</p>
<section><h2>Effect of electron hopping</h2>
<p>The electron density of states (DOS) and band dispersion are numerically computed [
<a href="#CR46"><sup>46</sup></a>
] and are shown in Figs. 
<a href="#Fig1"><sup>1</sup></a>
 and 
<a href="#Fig2"><sup>2</sup></a>
. The tight-binding calculations give the first-nearest-neighbour hopping integral 
<em>t</em>
<sub>1</sub>
 = 2.5–3.0 eV [
<a href="#CR30"><sup>30</sup></a>
, 
<a href="#CR31"><sup>31</sup></a>
]. In the present calculation, we have taken 
<em>ɛ</em>
<sub><em>a</em></sub>
 = 
<em>ɛ</em>
<sub><em>b</em></sub>
 = 0,
<span id="IEq42"><mml:math id="IEq42_Math"><mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>2.78</mml:mn><mml:mspace width="0.277778em"/><mml:mtext>eV</mml:mtext></mml:mrow></mml:math><tex-math id="IEq42_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$t_{1} = - 2.78 \;{\text{eV}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq42.gif"/></span>
, 
<span id="IEq43"><mml:math id="IEq43_Math"><mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>0.12</mml:mn><mml:mspace width="0.277778em"/><mml:mtext>eV</mml:mtext></mml:mrow></mml:math><tex-math id="IEq43_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$t_{2} = - 0.12\; {\text{eV}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq43.gif"/></span>
, 
<span id="IEq44"><mml:math id="IEq44_Math"><mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mn>3</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>0.068</mml:mn><mml:mspace width="0.277778em"/><mml:mtext>eV</mml:mtext></mml:mrow></mml:math><tex-math id="IEq44_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$t_{3} = - 0.068\; {\text{eV}}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq44.gif"/></span>
. Here, the scaled hopping integrals become 
<span id="IEq45"><mml:math id="IEq45_Math"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>1</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:math><tex-math id="IEq45_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{1} = - 1$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq45.gif"/></span>
, 
<span id="IEq46"><mml:math id="IEq46_Math"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>2</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>0.043</mml:mn></mml:mrow></mml:math><tex-math id="IEq46_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{2} = - 0.043$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq46.gif"/></span>
 and 
<span id="IEq47"><mml:math id="IEq47_Math"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>3</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>0.024</mml:mn></mml:mrow></mml:math><tex-math id="IEq47_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{3} = - 0.024$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq47.gif"/></span>
. The density of states (DOS) for electrons for the graphene is plotted for different band energies (c) (see Fig. 
<a href="#Fig1"><sup>1</sup></a>
). The DOS exhibits V-shaped gap structure at K-point (Dirac point) for nearest-neighbour hopping energy 
<span id="IEq48"><mml:math id="IEq48_Math"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>1</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:math><tex-math id="IEq48_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{1} = - 1$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq48.gif"/></span>
, i.e. DOS shows linear dependence of band energy. When second nearest-neighbour hopping (
<span id="IEq49"><mml:math id="IEq49_Math"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>2</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>0.043</mml:mn></mml:mrow></mml:math><tex-math id="IEq49_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{2} = - 0.043$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq49.gif"/></span>
) is included, it still retains the V-shape, but shifts to lower energies becoming asymmetric in nature with respect to Fermi level (
<span id="IEq50"><mml:math id="IEq50_Math"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ε</mml:mi><mml:mtext>F</mml:mtext></mml:msub><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:math><tex-math id="IEq50_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\varepsilon_{\text{F}} = 0$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq50.gif"/></span>
) at Dirac point (Fig. 
<a href="#Fig1"><sup>1</sup></a>
).
<figure id="Fig1"><h3>Fig. 1</h3>
<figcaption><p>Plot of the electronic density of states (DOS) of ideal graphene for different band energy with different hopping integrals, 
<span id="IEq78"><mml:math id="IEq78_Math"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>1</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:math><tex-math id="IEq78_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{1} = - 1$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq78.gif"/></span>
 (
<em>solid line</em>
) for NN; 
<span id="IEq79"><mml:math id="IEq79_Math"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>1</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:math><tex-math id="IEq79_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{1} = - 1$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq79.gif"/></span>
, 
<span id="IEq80"><mml:math id="IEq80_Math"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>2</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>0.043</mml:mn></mml:mrow></mml:math><tex-math id="IEq80_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{2} = - 0.043$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq80.gif"/></span>
 (
<em>dotted line</em>
) for NNN, 
<span id="IEq81"><mml:math id="IEq81_Math"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>1</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:math><tex-math id="IEq81_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{1} = - 1$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq81.gif"/></span>
, 
<span id="IEq82"><mml:math id="IEq82_Math"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>2</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>0.043</mml:mn></mml:mrow></mml:math><tex-math id="IEq82_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{2} = - 0.043$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq82.gif"/></span>
, 
<span id="IEq83"><mml:math id="IEq83_Math"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>3</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>0.024</mml:mn></mml:mrow></mml:math><tex-math id="IEq83_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{3} = - 0.024$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq83.gif"/></span>
 (
<em>dashed line</em>
) for NNNN</p></figcaption>
<img src="40089_2017_203_Fig1_HTML.gif" /></figure>
<figure id="Fig2"><h3>Fig. 2</h3>
<figcaption><p>Plot of the electronic band dispersion of ideal graphene for different band energy with different hopping integrals, 
<span id="IEq84"><mml:math id="IEq84_Math"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>1</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:math><tex-math id="IEq84_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{1} = - 1$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq84.gif"/></span>
 (
<em>solid line</em>
) for NN; 
<span id="IEq85"><mml:math id="IEq85_Math"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>1</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:math><tex-math id="IEq85_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{1} = - 1$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq85.gif"/></span>
, 
<span id="IEq86"><mml:math id="IEq86_Math"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>2</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>0.043</mml:mn></mml:mrow></mml:math><tex-math id="IEq86_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{2} = - 0.043$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq86.gif"/></span>
 (
<em>dotted line</em>
) for NNN, 
<span id="IEq87"><mml:math id="IEq87_Math"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>1</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:math><tex-math id="IEq87_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{1} = - 1$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq87.gif"/></span>
, 
<span id="IEq88"><mml:math id="IEq88_Math"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>2</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>0.043</mml:mn></mml:mrow></mml:math><tex-math id="IEq88_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{2} = - 0.043$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq88.gif"/></span>
,
<span id="IEq89"><mml:math id="IEq89_Math"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>3</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>0.024</mml:mn></mml:mrow></mml:math><tex-math id="IEq89_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{3} = - 0.024$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq89.gif"/></span>
 (
<em>dashed line</em>
) for NNNN</p></figcaption>
<img src="40089_2017_203_Fig2_HTML.gif" /></figure></p>
<p>The inclusion of third neighbour interaction induces a gap near K-point and the gap shifts further to lower energies. Figure 
<a href="#Fig2"><sup>2</sup></a>
 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 
<span id="IEq51"><mml:math id="IEq51_Math"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>1</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:math><tex-math id="IEq51_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{1} = - 1$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq51.gif"/></span>
, i.e. band dispersion shows linear dependence of band energy. When second nearest-neighbour hopping (
<span id="IEq52"><mml:math id="IEq52_Math"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>2</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>0.043</mml:mn></mml:mrow></mml:math><tex-math id="IEq52_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{2} = - 0.043$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq52.gif"/></span>
) is included, it still retains the V-shape, but shifts to lower energies becoming asymmetric in nature with respect to Fermi level (
<em>ɛ</em>
<sub>F</sub>
 = 0) at Dirac point (Fig. 
<a href="#Fig2"><sup>2</sup></a>
). The inclusion of third neighbour interaction induces a gap near K points and the gap shifts further to lower energies.</p></section>
<section><h2>Substrate effect</h2>
<p>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. (
<a href="#Equ1"><sup>1</sup></a>
) and (
<a href="#Equ3"><sup>3</sup></a>
) of our earlier publications [
<a href="#CR46"><sup>46</sup></a>
] and show the plots in Figs. 
<a href="#Fig3"><sup>3</sup></a>
 and 
<a href="#Fig4"><sup>4</sup></a>
. Figure 
<a href="#Fig3"><sup>3</sup></a>
 shows the DOS for different values of band gaps 
<em>d</em>
<sub>1</sub>
 = 0.035–0.107 developing in graphene sheet due to substrate effect [
<a href="#CR9"><sup>9</sup></a>
, 
<a href="#CR11"><sup>11</sup></a>
, 
<a href="#CR13"><sup>13</sup></a>
, 
<a href="#CR18"><sup>18</sup></a>
] and [
<a href="#CR32"><sup>32</sup></a>
–
<a href="#CR34"><sup>34</sup></a>
]. 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 
<em>ɛ</em>
<sub>F</sub>
 = 0. The gap is further enhanced with the increase of 
<em>d</em>
<sub>1</sub>
. Figure 
<a href="#Fig4"><sup>4</sup></a>
 shows the enhancement of band gap between lower electron occupied valence band and upper hole band due to the increase of 
<em>d</em>
<sub>1</sub>
. The middle of the band gap also moves down the Fermi level (
<em>ɛ</em>
<sub>F</sub>
 = 0).
<figure id="Fig3"><h3>Fig. 3</h3>
<figcaption><p>Plot of the variation of density of states (DOS) with energy (
<em>c</em>
) for different substrates 
<em>d</em>
<sub>1</sub>
 = 0.035, 0.071, 0.089, 0.0107</p></figcaption>
<img src="40089_2017_203_Fig3_HTML.gif" /></figure>
<figure id="Fig4"><h3>Fig. 4</h3>
<figcaption><p>Plot of the variation of energy band dispersion with momentum (
<em>k</em>
) for different substrates 
<em>d</em>
<sub>1</sub>
 = 0.035, 0.071, 0.089, 0.0107</p></figcaption>
<img src="40089_2017_203_Fig4_HTML.gif" /></figure></p></section>
<section><h2>Effect of Coulomb interaction</h2>
<p>The sub-lattice Coulomb interaction is treated within a mean-field approximation. Using the Hamiltonian given in Eqs. (
<a href="#Equ1"><sup>1</sup></a>
), (
<a href="#Equ3"><sup>3</sup></a>
) and (
<a href="#Equ4"><sup>4</sup></a>
) of our earlier calculations [
<a href="#CR47"><sup>47</sup></a>
], we calculate the difference (
<em>d</em>
) in occupation numbers as
<section id="Equ12"><mml:math display="block" id="Equ12_Math"><mml:mrow><mml:mi>d</mml:mi><mml:mo>=</mml:mo><mml:munder><mml:mo movablelimits="false">∑</mml:mo><mml:mi>k</mml:mi></mml:munder><mml:mfrac><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false">¯</mml:mo></mml:mrow></mml:mover><mml:mfenced close="]" open="[" separators=""><mml:mrow><mml:mi>f</mml:mi><mml:mfenced close=")" open="(" separators=""><mml:mrow><mml:mi mathvariant="italic">β</mml:mi><mml:msub><mml:mi mathvariant="italic">ω</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mi>k</mml:mi><mml:mo stretchy="false">↑</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:mfenced><mml:mo>-</mml:mo><mml:mi>f</mml:mi><mml:mfenced close=")" open="(" separators=""><mml:mrow><mml:mi mathvariant="italic">β</mml:mi><mml:msub><mml:mi mathvariant="italic">ω</mml:mi><mml:mrow><mml:mn>2</mml:mn><mml:mi>k</mml:mi><mml:mo stretchy="false">↑</mml:mo></mml:mrow></mml:msub></mml:mrow></mml:mfenced></mml:mrow></mml:mfenced></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:msub><mml:mi mathvariant="italic">ω</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mi>k</mml:mi><mml:mo stretchy="false">↑</mml:mo></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mi mathvariant="italic">ω</mml:mi><mml:mrow><mml:mn>2</mml:mn><mml:mi>k</mml:mi><mml:mo stretchy="false">↑</mml:mo></mml:mrow></mml:msub><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mfrac></mml:mrow></mml:math><tex-math id="Equ12_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$d = \mathop \sum \limits_{k} \frac{{\bar{\Delta } \left[ {f\left( {\beta \omega_{1k \uparrow } } \right) - f\left( {\beta \omega_{2k \uparrow } } \right)} \right]}}{{(\omega_{1k \uparrow } - \omega_{2k \uparrow } )}}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_Equ12.gif"/></section>
where 
<em>f</em>
 is the Fermi distribution function and 
<span id="IEq53"><mml:math id="IEq53_Math"><mml:msub><mml:mi mathvariant="italic">ω</mml:mi><mml:mrow><mml:mn>1</mml:mn><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mo stretchy="false">↑</mml:mo></mml:mrow></mml:msub></mml:math><tex-math id="IEq53_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\omega_{1k, \uparrow }$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq53.gif"/></span>
 and 
<em>ω</em>
<sub>2
<em>k</em>
,↑</sub>
 are the electron band dispersions. The modified gap 
<span id="IEq54"><mml:math id="IEq54_Math"><mml:mrow><mml:msub><mml:mi>d</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mfenced close=")" open="("><mml:mi>T</mml:mi></mml:mfenced><mml:mo>=</mml:mo><mml:mfrac><mml:mover accent="true"><mml:mi mathvariant="normal">Δ</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:msub><mml:mi>t</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mfrac></mml:mrow></mml:math><tex-math id="IEq54_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$d_{2} \left( T \right) = \frac{{\tilde{\Delta }}}{{t_{1} }}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq54.gif"/></span>
 is plotted in Figs. 
<a href="#Fig5"><sup>5</sup></a>
, 
<a href="#Fig6"><sup>6</sup></a>
 for different Coulomb energies and hopping parameters. The ferromagnetic magnetizations and spin polarizations are reported [
<a href="#CR51"><sup>51</sup></a>
].
<figure id="Fig5"><h3>Fig. 5</h3>
<figcaption><p>Plot of the variation of modified band gap (
<em>d</em>
<sub>2</sub>
) with temperature (
<em>t</em>
) for different values of Coulomb energy 
<em>u</em>
 = 1.0, 1.5, 1.7, 2.0, 2.5, 3.5, 4.5 and 5.0 for fixed substrate-induced band gap 
<em>d</em>
<sub>1</sub>
 = 0.1</p></figcaption>
<img src="40089_2017_203_Fig5_HTML.gif" /></figure>
<figure id="Fig6"><h3>Fig. 6</h3>
<figcaption><p>Plot of the variation of modified band gap (
<em>d</em>
<sub>2</sub>
) with temperature (
<em>t</em>
) at 
<em>u</em>
 = 1.7 for different values of electron hopping 
<span id="IEq90"><mml:math id="IEq90_Math"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>1</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:math><tex-math id="IEq90_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{1} = - 1$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq90.gif"/></span>
, 
<span id="IEq91"><mml:math id="IEq91_Math"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>2</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>0.043</mml:mn></mml:mrow></mml:math><tex-math id="IEq91_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{2} = - 0.043$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq91.gif"/></span>
, 
<span id="IEq92"><mml:math id="IEq92_Math"><mml:mrow><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>3</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mn>0.024</mml:mn></mml:mrow></mml:math><tex-math id="IEq92_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{3} = - 0.024$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq92.gif"/></span>
 for fixed substrate-induced band gap 
<em>d</em>
<sub>1</sub>
 = 0.1. The band dispersions for different hopping are also given in inset of the figure</p></figcaption>
<img src="40089_2017_203_Fig6_HTML.gif" /></figure></p>
<p>The graphene acquires a band gap of ∆ due to substrate effect. The effective band gap then becomes 
<span id="IEq55"><mml:math id="IEq55_Math"><mml:mrow><mml:mover accent="true"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false">¯</mml:mo></mml:mrow></mml:mover><mml:mo>=</mml:mo><mml:mi mathvariant="normal">Δ</mml:mi><mml:mo>+</mml:mo><mml:mi>U</mml:mi><mml:mfrac><mml:mi>d</mml:mi><mml:mn>2</mml:mn></mml:mfrac></mml:mrow></mml:math><tex-math id="IEq55_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\bar{\Delta } = \Delta + U\frac{d}{2}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq55.gif"/></span>
 due to Coulomb interaction between electrons. We scale the physical parameters by nearest-neighbour hopping integral (
<em>t</em>
<sub>1</sub>
) and hence the modified gap appears as 
<span id="IEq56"><mml:math id="IEq56_Math"><mml:mrow><mml:msub><mml:mi>d</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>d</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mo>+</mml:mo><mml:mi>u</mml:mi><mml:mfrac><mml:mi>d</mml:mi><mml:mn>2</mml:mn></mml:mfrac></mml:mrow></mml:math><tex-math id="IEq56_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$d_{2} = d_{1} + u\frac{d}{2}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq56.gif"/></span>
. The temperature dependent difference (
<em>d</em>
) 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 (
<span id="IEq57"><mml:math id="IEq57_Math"><mml:mrow><mml:msub><mml:mi>d</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mover accent="true"><mml:mrow><mml:mi mathvariant="normal">Δ</mml:mi></mml:mrow><mml:mrow><mml:mo stretchy="false">¯</mml:mo></mml:mrow></mml:mover><mml:msub><mml:mi>t</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mfrac></mml:mrow></mml:math><tex-math id="IEq57_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$d_{2} = \frac{{\bar{\Delta }}}{{t_{1} }}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq57.gif"/></span>
) is computed for different Coulomb correlation energies as shown in Fig. 
<a href="#Fig5"><sup>5</sup></a>
. For lower Coulomb energies, the modified gap (
<em>d</em>
<sub>2</sub>
) at 
<em>t</em>
 = 0 gradually increases with increase of 
<em>u</em>
 from 1 to 1.7 and attains the maximum of 
<span id="IEq58"><mml:math id="IEq58_Math"><mml:mrow><mml:msub><mml:mi>d</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mn>0.182</mml:mn></mml:mrow></mml:math><tex-math id="IEq58_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$d_{2} = 0.182$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq58.gif"/></span>
. On further increasing to higher Coulomb energies, the modified gap (
<em>d</em>
<sub>2</sub>
) at 
<em>t</em>
 = 0 decreases and attains the bare gap (
<em>d</em>
<sub>1</sub>
) arising due to only substrate effect for Coulomb energy 
<em>u</em>
 = 5.0. It is observed that the magnetic gap vanishes for 
<em>u</em>
 = 5.0 indicating that 
<span id="IEq59"><mml:math id="IEq59_Math"><mml:mrow><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mo stretchy="false">↑</mml:mo></mml:mrow><mml:mi>a</mml:mi></mml:msubsup><mml:mo>=</mml:mo><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mo stretchy="false">↑</mml:mo></mml:mrow><mml:mi>b</mml:mi></mml:msubsup></mml:mrow></mml:math><tex-math id="IEq59_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$n_{ \uparrow }^{a} = n_{ \uparrow }^{b}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq59.gif"/></span>
 in paramagnetic phase. The temperature dependence of modified gap shows that the magnitude of gap for lower ‘
<em>u</em>
’ gradually decreases with temperature. For higher 
<em>u</em>
, 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 ‘
<em>u</em>
’ at higher temperatures indicating that 
<span id="IEq60"><mml:math id="IEq60_Math"><mml:mrow><mml:mi>d</mml:mi><mml:mo>=</mml:mo><mml:mo>&lt;</mml:mo><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mo stretchy="false">↑</mml:mo></mml:mrow><mml:mi>a</mml:mi></mml:msubsup><mml:mo>&gt;</mml:mo><mml:mo>-</mml:mo><mml:mo>&lt;</mml:mo><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mo stretchy="false">↑</mml:mo></mml:mrow><mml:mi>b</mml:mi></mml:msubsup><mml:mo>&gt;</mml:mo></mml:mrow></mml:math><tex-math id="IEq60_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$d = < n_{ \uparrow }^{a} > - < n_{ \uparrow }^{b} >$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq60.gif"/></span>
 remains unchanged at very high temperatures and high Coulomb interactions.</p>
<p>Figure 
<a href="#Fig6"><sup>6</sup></a>
 shows the effect of different hopping integrals on temperature-dependent modified gap (
<em>d</em>
<sub>2</sub>
). In the absence of Coulomb interaction and electron hopping, the band gap in graphene is 
<em>d</em>
<sub>1</sub>
 = 0.1 which arises due to substrate effect only. This band gap at temperature 0 K is enhanced to 
<em>d</em>
<sub>2</sub>
 = 0.167 due to the onset of a critical Coulomb energy 
<em>u</em>
<sub>c</sub>
 = 1.7 for the nearest-neighbour hopping energy of 
<em>t</em>
<sub>1</sub>
 = −1.0. This gap at temperature 0 K is further enhanced to 
<em>d</em>
<sub>2</sub>
 = 0.181 for the same critical Coulomb energy 
<em>u</em>
<sub>c</sub>
 = 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. 
<a href="#Fig6"><sup>6</sup></a>
.</p></section>
<section><h2>Effect of electron–phonon interaction</h2>
<p>The electrons on graphene sheet interact strongly with the phonons on the polarized surface of substrates. Applying Lang–Firsov canonical transformation [
<a href="#CR48"><sup>48</sup></a>
] to Eqs. (
<a href="#Equ1"><sup>1</sup></a>
), (
<a href="#Equ3"><sup>3</sup></a>
), (
<a href="#Equ4"><sup>4</sup></a>
), (
<a href="#Equ5"><sup>5</sup></a>
) and (
<a href="#Equ6"><sup>6</sup></a>
) in high frequency limit of localized phonons, we calculate the effective Coulomb energy 
<span id="IEq61"><mml:math id="IEq61_Math"><mml:mrow><mml:mover accent="true"><mml:mi>U</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mo>=</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>U</mml:mi><mml:mo>-</mml:mo><mml:mn>2</mml:mn><mml:mi mathvariant="italic">λ</mml:mi><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>1</mml:mn></mml:msub><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:math><tex-math id="IEq61_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{U} = (U - 2\lambda \tilde{t}_{1} )$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq61.gif"/></span>
 and effective hopping (
<span id="IEq62"><mml:math id="IEq62_Math"><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>1</mml:mn></mml:msub></mml:math><tex-math id="IEq62_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{1}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq62.gif"/></span>
) which is a function of phonon frequency 
<em>ω</em>
<sub>0</sub>
, electron–phonon coupling (
<em>λ</em>
) and NN hopping integral (
<span id="IEq63"><mml:math id="IEq63_Math"><mml:msub><mml:mover accent="true"><mml:mi>t</mml:mi><mml:mo stretchy="false">~</mml:mo></mml:mover><mml:mn>1</mml:mn></mml:msub></mml:math><tex-math id="IEq63_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\tilde{t}_{1}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq63.gif"/></span>
). The temperature-dependent modified gaps are plotted in Figs. 
<a href="#Fig7"><sup>7</sup></a>
 and 
<a href="#Fig8"><sup>8</sup></a>
.
<figure id="Fig7"><h3>Fig. 7</h3>
<figcaption><p>Plot of the variation of modified band gap (
<em>d</em>
<sub>2</sub>
) vs. phonon frequency 
<em>ω</em>
<sub>0</sub>
 for different values of Coulomb potential 
<span id="IEq93"><mml:math id="IEq93_Math"><mml:mrow><mml:mi>u</mml:mi><mml:mo>=</mml:mo><mml:mn>0.5</mml:mn><mml:mo>,</mml:mo><mml:mn>1.0</mml:mn><mml:mo>,</mml:mo><mml:mn>1.5</mml:mn><mml:mo>,</mml:mo><mml:mn>2.0</mml:mn><mml:mo>,</mml:mo><mml:mn>2.5</mml:mn></mml:mrow></mml:math><tex-math id="IEq93_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$u = 0.5, 1.0, 1.5, 2.0, 2.5$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq93.gif"/></span>
 for fixed electron–phonon coupling constant 
<span id="IEq94"><mml:math id="IEq94_Math"><mml:mrow><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:mo>=</mml:mo><mml:mn>0.2</mml:mn></mml:mrow></mml:math><tex-math id="IEq94_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$al = 0.2$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq94.gif"/></span>
, temperature 
<span id="IEq95"><mml:math id="IEq95_Math"><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mn>0.01</mml:mn></mml:mrow></mml:math><tex-math id="IEq95_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$(t) = 0.01$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq95.gif"/></span>
 and substrate-induced gap 
<em>d</em>
<sub>1</sub>
 = 0.1</p></figcaption>
<img src="40089_2017_203_Fig7_HTML.gif" /></figure>
<figure id="Fig8"><h3>Fig. 8</h3>
<figcaption><p>Plot of the variation of modified band gap (
<em>d</em>
<sub>2</sub>
) vs. temperature (
<em>t</em>
) for different values of phonon frequency 
<span id="IEq96"><mml:math id="IEq96_Math"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ω</mml:mi><mml:mn>0</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mn>0.2</mml:mn><mml:mo>,</mml:mo><mml:mn>0.4</mml:mn><mml:mo>,</mml:mo><mml:mn>0.6</mml:mn><mml:mo>,</mml:mo><mml:mn>0.8</mml:mn><mml:mo>,</mml:mo><mml:mn>1.0</mml:mn></mml:mrow></mml:math><tex-math id="IEq96_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\omega_{0} = 0.2, 0.4, 0.6, 0.8, 1.0$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq96.gif"/></span>
 for fixed electron–phonon coupling constant 
<em>al</em>
 = 0.2, Coulomb interaction 
<span id="IEq97"><mml:math id="IEq97_Math"><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mi>u</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>=</mml:mo><mml:mn>1.7</mml:mn></mml:mrow></mml:math><tex-math id="IEq97_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$(u) = 1.7$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq97.gif"/></span>
 and substrate-induced gap (
<em>d</em>
<sub>1</sub>
) = 0.1</p></figcaption>
<img src="40089_2017_203_Fig8_HTML.gif" /></figure></p>
<p>The effect of Coulomb interaction (
<em>u</em>
) on the modified gap 
<em>d</em>
<sub>2</sub>
 for high phonon frequency (
<em>ω</em>
<sub>0</sub>
) vibration is shown in Fig. 
<a href="#Fig7"><sup>7</sup></a>
. 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 
<em>d</em>
<sub>1</sub>
 = 0.1 and attains maximum value 
<em>d</em>
<sub>2</sub>
 = 0.165 for given Coulomb interaction 
<em>u</em>
 = 2.0. With further increase in Coulomb interaction, the modified gap gradually increases with phonon frequency 
<em>ω</em>
<sub>0</sub>
 attains a maximum flat peak with 
<em>d</em>
<sub>2</sub>
 = 0.165 for critical Coulomb interaction 
<em>u</em>
<sub>c</sub>
 = 2.5 for critical phonon frequency 
<em>ω</em>
<sub>0c</sub>
 ≈  0.35. Thus, it is clear that modulated gap is maximized for critical Coulomb interaction 
<em>u</em>
<sub>c</sub>
 = 2.5 and critical phonon frequency 
<em>ω</em>
<sub>0c</sub>
 = 0.35 for given values of electron–phonon coupling 
<em>al</em>
 = 0.2. The critical Coulomb interaction 
<em>U</em>
<sub>C</sub>
 = 2.5
<em>t</em>
<sub>1</sub>
 for producing maximum modified gap is slightly higher than the critical Coulomb interaction 
<span id="IEq64"><mml:math id="IEq64_Math"><mml:mrow><mml:msub><mml:mi>U</mml:mi><mml:mtext>c</mml:mtext></mml:msub><mml:mo>=</mml:mo><mml:msub><mml:mi>u</mml:mi><mml:mtext>c</mml:mtext></mml:msub><mml:msub><mml:mi>t</mml:mi><mml:mn>1</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mn>1.7</mml:mn><mml:msub><mml:mi>t</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:math><tex-math id="IEq64_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$U_{\text{c}} = u_{\text{c}} t_{1} = 1.7 t_{1}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq64.gif"/></span>
 obtained earlier in the absence of electron–phonon interaction and phonon frequency [
<a href="#CR47"><sup>47</sup></a>
].</p>
<p>The effect of phonon frequency (
<em>ω</em>
<sub>0</sub>
) on gap is shown in Fig. 
<a href="#Fig8"><sup>8</sup></a>
. For a given low electron–phonon coupling 
<span id="IEq65"><mml:math id="IEq65_Math"><mml:mrow><mml:mi>a</mml:mi><mml:mi>l</mml:mi><mml:mo>=</mml:mo><mml:mn>0.2</mml:mn></mml:mrow></mml:math><tex-math id="IEq65_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$al = 0.2$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq65.gif"/></span>
 and relatively low phonon frequency 
<em>ω</em>
<sub>0</sub>
 = 0.2, the modified gap becomes 
<em>d</em>
<sub>2</sub>
 = 0.13 gap at temperature 
<em>t</em>
 = 0. With increase of phonon vibrational frequency, the modified gap is enhanced to the higher value, i.e. 
<em>d</em>
<sub>2</sub>
 = 0.165 at temperature 
<em>t</em>
 = 0 for vibration frequency 
<em>ω</em>
<sub>0</sub>
 = 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.</p></section>
<section><h2>Effect of bi-layer graphene (BLG)</h2>
<p>Experiments [
<a href="#CR8"><sup>8</sup></a>
, 
<a href="#CR20"><sup>20</sup></a>
, 
<a href="#CR23"><sup>23</sup></a>
] 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 
<em>t</em>
<sub>1</sub>
 and 
<em>t</em>
<sub>⊥</sub>
, respectively, in the presence of gating potential V [
<a href="#CR49"><sup>49</sup></a>
, 
<a href="#CR52"><sup>52</sup></a>
, 
<a href="#CR53"><sup>53</sup></a>
]. The four bands for BLG are written as
<section id="Equ13"><mml:math display="block" id="Equ13_Math"><mml:mrow><mml:msub><mml:mi mathvariant="italic">ω</mml:mi><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mo>,</mml:mo><mml:mi>s</mml:mi><mml:mo>,</mml:mo></mml:mrow></mml:msub><mml:mfenced close=")" open="("><mml:mi>k</mml:mi></mml:mfenced><mml:mo>=</mml:mo><mml:mo>-</mml:mo><mml:mi mathvariant="italic">μ</mml:mi><mml:mo>-</mml:mo><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mo>-</mml:mo><mml:mn>1</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mi>s</mml:mi></mml:msup><mml:msqrt><mml:mfrac><mml:mrow><mml:mo stretchy="false">[</mml:mo><mml:mfrac><mml:msup><mml:mi>V</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mn>2</mml:mn></mml:mfrac><mml:mo>+</mml:mo><mml:mn>2</mml:mn><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup><mml:mo>+</mml:mo><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mo>,</mml:mo><mml:mi>⊥</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup><mml:mo>-</mml:mo><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mo>-</mml:mo><mml:mn>1</mml:mn><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mi mathvariant="italic">α</mml:mi></mml:msup><mml:msqrt><mml:mi>R</mml:mi></mml:msqrt><mml:mo stretchy="false">]</mml:mo></mml:mrow><mml:mn>2</mml:mn></mml:mfrac></mml:msqrt></mml:mrow></mml:math><tex-math id="Equ13_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\omega_{\alpha ,s,} \left( k \right) = - \mu - ( - 1)^{s} \sqrt {\frac{{[ \frac{{V^{2} }}{2} + 2\varepsilon_{k}^{2} + \varepsilon_{k, \bot }^{2} - ( - 1)^{\alpha } \sqrt R ]}}{2}}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_Equ13.gif"/></section>
<section id="Equ14"><mml:math display="block" id="Equ14_Math"><mml:mrow><mml:mi>R</mml:mi><mml:mo>=</mml:mo><mml:msqrt><mml:mrow><mml:msup><mml:mfenced close=")" open="(" separators=""><mml:mrow><mml:msup><mml:mi>V</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mo stretchy="false">/</mml:mo><mml:mn>2</mml:mn><mml:mo>+</mml:mo><mml:mn>2</mml:mn><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup><mml:mo>+</mml:mo><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mi>⊥</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup></mml:mrow></mml:mfenced><mml:mn>2</mml:mn></mml:msup><mml:mo>-</mml:mo><mml:mn>4</mml:mn><mml:mfenced close=")" open="(" separators=""><mml:mrow><mml:msup><mml:mi>V</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mo stretchy="false">/</mml:mo><mml:mn>16</mml:mn><mml:mo>+</mml:mo><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup><mml:msup><mml:mi>V</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mo stretchy="false">/</mml:mo><mml:mn>2</mml:mn><mml:mo>+</mml:mo><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi>k</mml:mi><mml:mi>⊥</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msubsup><mml:msup><mml:mi>V</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:mo stretchy="false">/</mml:mo><mml:mn>4</mml:mn><mml:mo>+</mml:mo><mml:msubsup><mml:mi mathvariant="italic">ε</mml:mi><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mn>4</mml:mn></mml:msubsup></mml:mrow></mml:mfenced></mml:mrow></mml:msqrt></mml:mrow></mml:math><tex-math id="Equ14_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$R = \sqrt {\left( {V^{2} /2 + 2\varepsilon_{k}^{2} + \varepsilon_{k \bot }^{2} } \right)^{2} - 4\left( {V^{2} /16 + \varepsilon_{k}^{2} V^{2} /2 + \varepsilon_{k \bot }^{2} V^{2} /4 + \varepsilon_{k}^{4} } \right)}$$\end{document}]]></tex-math><graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_Equ14.gif"/></section>
where 
<span id="IEq66"><mml:math id="IEq66_Math"><mml:mrow><mml:mi>s</mml:mi><mml:mo>,</mml:mo><mml:mi mathvariant="italic">α</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:math><tex-math id="IEq66_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$s,\alpha = 1 , 2$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq66.gif"/></span>
. The band dispersion is plotted in Fig. 
<a href="#Fig9"><sup>9</sup></a>
 which shows four bands 
<span id="IEq67"><mml:math id="IEq67_Math"><mml:msub><mml:mi mathvariant="italic">ω</mml:mi><mml:mrow><mml:mi mathvariant="italic">α</mml:mi><mml:mo>,</mml:mo><mml:mi>s</mml:mi><mml:mi>k</mml:mi></mml:mrow></mml:msub></mml:math><tex-math id="IEq67_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\omega_{\alpha ,sk}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq67.gif"/></span>
 of which 
<span id="IEq68"><mml:math id="IEq68_Math"><mml:msub><mml:mi mathvariant="italic">ω</mml:mi><mml:mrow><mml:mn>11</mml:mn><mml:mi>k</mml:mi></mml:mrow></mml:msub></mml:math><tex-math id="IEq68_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\omega_{11k}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq68.gif"/></span>
 and 
<em>ω</em>
<sub>22
<em>k</em></sub>
 are high energy bands and 
<em>ω</em>
<sub>12
<em>k</em></sub>
 (conduction) and 
<em>ω</em>
<sub>21
<em>k</em></sub>
 (valence) bands are the low lying dispersion bands near Dirac point with zero energy. In the absence of gate potential 
<span id="IEq69"><mml:math id="IEq69_Math"><mml:mrow><mml:mi>v</mml:mi><mml:mo>=</mml:mo><mml:mfrac><mml:mi>V</mml:mi><mml:msub><mml:mi>t</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mfrac><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:math><tex-math id="IEq69_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$v = \frac{V}{{t_{1} }} = 0$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq69.gif"/></span>
, low energy bands touch at Dirac point. For low electric potential (
<span id="IEq70"><mml:math id="IEq70_Math"><mml:mrow><mml:mi>v</mml:mi><mml:mo>=</mml:mo><mml:mn>0.054</mml:mn></mml:mrow></mml:math><tex-math id="IEq70_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$v = 0.054$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq70.gif"/></span>
), the low energy bands show a Mexican hat shape with a gap energy of 0.054 as observed experimentally [
<a href="#CR8"><sup>8</sup></a>
, 
<a href="#CR23"><sup>23</sup></a>
].
<figure id="Fig9"><h3>Fig. 9</h3>
<figcaption><p>Plot of energy band dispersions vs. momentum (k) for pristine bi-layer graphene for different values of external biasing 
<em>v</em>
 = 0 (
<em>dashed line</em>
) and 
<em>v</em>
 = 0.054 (
<em>solid line</em>
)</p></figcaption>
<img src="40089_2017_203_Fig9_HTML.gif" /></figure></p></section></section>
<section><h2>Conclusions</h2>
<p>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 
<em>t</em>
<sub>1</sub>
. 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.</p>
<p>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 
<em>T</em>
 = 0
<sup>0</sup>
<em>K</em>
 with the increase of Coulomb interaction. The effective gap attains its maximum value, 
<em>d</em>
<sub>2</sub>
 = 0.182 at 
<em>t</em>
 = 0 for critical Coulomb interaction 
<em>U</em>
<sub>C</sub>
 = 1.7
<em>t</em>
<sub>1</sub>
. On further increase in the Coulomb energy to 
<span id="IEq71"><mml:math id="IEq71_Math"><mml:mrow><mml:mi>U</mml:mi><mml:mo>=</mml:mo><mml:mn>5.0</mml:mn><mml:msub><mml:mi>t</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:math><tex-math id="IEq71_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$U = 5.0 t_{1}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq71.gif"/></span>
, we observe that the modified gap d
<sub>2</sub>
 = d
<sub>1</sub>
 indicating that the electron occupancies of both sub-lattices becoming same (i.e. 
<span id="IEq72"><mml:math id="IEq72_Math"><mml:mrow><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mo stretchy="false">↑</mml:mo></mml:mrow><mml:mi>a</mml:mi></mml:msubsup><mml:mo>=</mml:mo><mml:msubsup><mml:mi>n</mml:mi><mml:mrow><mml:mo stretchy="false">↑</mml:mo></mml:mrow><mml:mi>b</mml:mi></mml:msubsup></mml:mrow></mml:math><tex-math id="IEq72_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$n_{ \uparrow }^{a} = n_{ \uparrow }^{b}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq72.gif"/></span>
) 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 (
<span id="IEq73"><mml:math id="IEq73_Math"><mml:mrow><mml:mtext>i</mml:mtext><mml:mo>.</mml:mo><mml:mtext>e</mml:mtext><mml:mo>.</mml:mo><mml:mspace width="0.277778em"/><mml:msub><mml:mi>d</mml:mi><mml:mn>2</mml:mn></mml:msub><mml:mo>=</mml:mo><mml:mn>0.165</mml:mn></mml:mrow></mml:math><tex-math id="IEq73_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$${\text{i}} . {\text{e}} . \;d_{2} = 0.165$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq73.gif"/></span>
) for a critical Coulomb interaction 
<span id="IEq74"><mml:math id="IEq74_Math"><mml:mrow><mml:msub><mml:mi>U</mml:mi><mml:mtext>c</mml:mtext></mml:msub><mml:mo>=</mml:mo><mml:mn>2.5</mml:mn><mml:msub><mml:mi>t</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:math><tex-math id="IEq74_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$U_{\text{c}} = 2.5 t_{1}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq74.gif"/></span>
 which is higher than the critical Coulomb interaction 
<span id="IEq75"><mml:math id="IEq75_Math"><mml:mrow><mml:msub><mml:mi>U</mml:mi><mml:mtext>c</mml:mtext></mml:msub><mml:mo>=</mml:mo><mml:mn>1.7</mml:mn><mml:msub><mml:mi>t</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:math><tex-math id="IEq75_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$U_{\text{c}} = 1.7 t_{1}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq75.gif"/></span>
 obtained earlier in the absence of electron–phonon interaction and phonon vibration [
<a href="#CR47"><sup>47</sup></a>
]. 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 
<span id="IEq76"><mml:math id="IEq76_Math"><mml:mrow><mml:mi>V</mml:mi><mml:mo>=</mml:mo><mml:mn>0.054</mml:mn><mml:msub><mml:mi>t</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:math><tex-math id="IEq76_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$V = 0.054 t_{1}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq76.gif"/></span>
. Since we have considered dispersion for full Brillouin zone in both the layers, we can take any gating potential higher than the gating potential 
<span id="IEq77"><mml:math id="IEq77_Math"><mml:mrow><mml:mi>V</mml:mi><mml:mo>=</mml:mo><mml:mn>0.054</mml:mn><mml:msub><mml:mi>t</mml:mi><mml:mn>1</mml:mn></mml:msub></mml:mrow></mml:math><tex-math id="IEq77_TeX"><![CDATA[\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym}
				\usepackage{amsfonts}
				\usepackage{amssymb}
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$V = 0.054 t_{1}$$\end{document}]]></tex-math><inline-graphic specific-use="web" mime-subtype="GIF" xlink:href="40089_2017_203_Article_IEq77.gif"/></span>
 applied to bi-layer systems experimentally [
<a href="#CR8"><sup>8</sup></a>
, 
<a href="#CR23"><sup>23</sup></a>
].</p></section><hr/><h2>Acknowledgements</h2>
<p>The authors gracefully acknowledge the research facilities offered by the Institute of Physics, Bhubaneswar, India.</p>
<hr/><h2>References</h2>
<ol><li>Novoselov and Geim (2004) <em>Electric field effect in atomically thin carbon films</em> <a href="https://doi.org/10.1126/science.1102896" target="_blank">10.1126/science.1102896</a></li><li>Novoselov et al. (2005) <em>Two-dimensional gas of massless Dirac fermions in graphene</em> <a href="https://doi.org/10.1038/nature04233" target="_blank">10.1038/nature04233</a></li><li>Wallace (1947) <em>The Band Theory of Graphite</em> <a href="https://doi.org/10.1103/PhysRev.71.622" target="_blank">10.1103/PhysRev.71.622</a></li><li>Miller and Kubista (2009) <em>Observing the Quantization of ZeroMass Carriers in Graphene</em> <a href="https://doi.org/10.1126/science.1171810" target="_blank">10.1126/science.1171810</a></li><li>Stander et al. (2009) <em>Evidence of Klein tunneling in graphene pn junctions</em> <a href="https://doi.org/10.1103/PhysRevLett.102.026807" target="_blank">10.1103/PhysRevLett.102.026807</a></li><li>Bolotin et al. (2008) <em>Temperature-Dependent Transport in Suspended Graphene</em> <a href="https://doi.org/10.1103/PhysRevLett.101.096802" target="_blank">10.1103/PhysRevLett.101.096802</a></li><li>Han and Özyilmaz (2007) <em>Energy Band-Gap Engineering of Graphene Nanoribbons</em> <a href="https://doi.org/10.1103/PhysRevLett.98.206805" target="_blank">10.1103/PhysRevLett.98.206805</a></li><li>McCann and Falko (2006) <em>Landau level degeneracy and quantum Hall effect in a graphite bilayer</em> <a href="https://doi.org/10.1103/PhysRevLett.96.086805" target="_blank">10.1103/PhysRevLett.96.086805</a></li><li>Giovannetti and Khomyakov (2007) <em>Substrate-induced band gap in graphene on hexagonal boron nitride: Ab initio density functional calculations</em> <a href="https://doi.org/10.1103/PhysRevB.76.073103" target="_blank">10.1103/PhysRevB.76.073103</a></li><li>Jung et al. (2015) <em>Origin of band gaps in graphene on hexagonal boron nitride</em> <a href="https://doi.org/10.1038/ncomms7308" target="_blank">10.1038/ncomms7308</a></li><li>Zhou et al. (2007) <em>Substrate-induced bandgap opening in epitaxial graphene</em> <a href="https://doi.org/10.1038/nmat2003" target="_blank">10.1038/nmat2003</a></li><li>Nevius et al. (2015) <em>Semiconducting Graphene from Highly Ordered Substrate Interactions</em> <a href="https://doi.org/10.1103/PhysRevLett.115.136802" target="_blank">10.1103/PhysRevLett.115.136802</a></li><li>Varchon et al. (2007) <em>Electronic Structure of Epitaxial Graphene Layers on SiC: Effect of the Substrate</em> <a href="https://doi.org/10.1103/PhysRevLett.99.126805" target="_blank">10.1103/PhysRevLett.99.126805</a></li><li>Mattausch and Pankratov (2007) <em>Ab Initio Study of Graphene on SiC</em> <a href="https://doi.org/10.1103/PhysRevLett.99.076802" target="_blank">10.1103/PhysRevLett.99.076802</a></li><li>Emtsev et al. (2008) <em>Interaction, growth, and ordering of epitaxial graphene on SiC{0001} surfaces: A comparative photoelectron spectroscopy study</em> <a href="https://doi.org/10.1103/PhysRevB.77.155303" target="_blank">10.1103/PhysRevB.77.155303</a></li><li>Varchon et al. (2008) <em>Ripples in epitaxial graphene on the Si-terminated SiC(0001) surface</em> <a href="https://doi.org/10.1103/PhysRevB.77.235412" target="_blank">10.1103/PhysRevB.77.235412</a></li><li>Varchon et al. (2008) <em>Why Multilayer Graphene on 4H−SiC(000–1) Behaves Like a Single Sheet of Graphene</em> <a href="https://doi.org/10.1103/PhysRevLett.100.125504" target="_blank">10.1103/PhysRevLett.100.125504</a></li><li>Fratini and Guinea (2008) <em>Substrate-limited electron dynamics in graphene</em> <a href="https://doi.org/10.1103/PhysRevB.77.195415" target="_blank">10.1103/PhysRevB.77.195415</a></li><li>Zhao et al. (2015) <em>One-step green synthesis of a ruthenium/graphene composite as a highly efficient catalyst</em> <a href="https://doi.org/10.1039/C4RA11397G" target="_blank">10.1039/C4RA11397G</a></li><li>McCann and Abergel (2007) <em>Electrons in bilayer graphene</em> <a href="https://doi.org/10.1016/j.ssc.2007.03.054" target="_blank">10.1016/j.ssc.2007.03.054</a></li><li>Zhang et al. (2009) <em>Direct observation of a widely tunable bandgap in bilayer graphene</em> <a href="https://doi.org/10.1038/nature08105" target="_blank">10.1038/nature08105</a></li><li>Abergel (2015) <em>Mucha-Kruczyński, Marcin.: Infrared absorption of closely aligned heterostructures of monolayer and bilayer graphene with hexagonal boron nitride</em> <a href="https://doi.org/10.1103/PhysRevB.92.115430" target="_blank">10.1103/PhysRevB.92.115430</a></li><li>Ohta et al. (2006) <em>Controlling the electronic structure of bilayer graphene</em> <a href="https://doi.org/10.1126/science.1130681" target="_blank">10.1126/science.1130681</a></li><li>Hague (2011) <em>Tunable graphene band gaps from superstrate-mediated interactions</em> <a href="https://doi.org/10.1103/PhysRevB.84.155438" target="_blank">10.1103/PhysRevB.84.155438</a></li><li>Hague (2012) <em>Polarons in highly doped atomically thin graphitic materials</em> <a href="https://doi.org/10.1103/PhysRevB.86.064302" target="_blank">10.1103/PhysRevB.86.064302</a></li><li>Halasall and Ferrari (2009) <em>Control of graphene's properties by reversible hydrogenation: evidence for graphane</em> <a href="https://doi.org/10.1126/science.1167130" target="_blank">10.1126/science.1167130</a></li><li>Eklund et al. (2010) <em>Reversible fluorination of graphene: Evidence of a two-dimensional wide bandgap semiconductor</em> <a href="https://doi.org/10.1103/PhysRevB.81.205435" target="_blank">10.1103/PhysRevB.81.205435</a></li><li>Lin et al. (2013) <em>Ab initio optical study of graphene on hexagonal boron nitride and fluorographene substrates</em> <a href="https://doi.org/10.1039/c2tc00570k" target="_blank">10.1039/c2tc00570k</a></li><li>Song et al. (2010) <em>Large scale growth and characterization of atomic hexagonal boron nitride layers</em> <a href="https://doi.org/10.1021/nl1022139" target="_blank">10.1021/nl1022139</a></li><li>Wehling et al. (2011) <em>Strength of Effective Coulomb Interactions in Graphene and Graphite</em> <a href="https://doi.org/10.1103/PhysRevLett.106.236805" target="_blank">10.1103/PhysRevLett.106.236805</a></li><li>Kotov and Uchoa (2012) <em>Electron-Electron Interactions in Graphene: Current Status and Perspectives</em> <a href="https://doi.org/10.1103/RevModPhys.84.1067" target="_blank">10.1103/RevModPhys.84.1067</a></li><li>Hansmann et al. (2013) <em>Long-Range Coulomb Interactions in Surface Systems: A First-Principles Description within Self-Consistently Combined GW and Dynamical Mean-Field Theory</em> <a href="https://doi.org/10.1103/PhysRevLett.110.166401" target="_blank">10.1103/PhysRevLett.110.166401</a></li><li>Pariser and Parr (1953) <em>A Semi-Empirical Theory of the Electronic Spectra and Electronic Structure of Complex Unsaturated Molecules</em> <a href="https://doi.org/10.1063/1.1699030" target="_blank">10.1063/1.1699030</a></li><li>Soos et al. (1993) <em>Band to correlated crossover in alternating Hubbard and Pariser-Parr-Pople chains: Nature of the lowest singlet excitation of conjugated polymers</em> <a href="https://doi.org/10.1103/PhysRevLett.71.1609" target="_blank">10.1103/PhysRevLett.71.1609</a></li><li>Friend et al. (1999) <em>Electroluminescence in conjugated polymers</em> <a href="https://doi.org/10.1038/16393" target="_blank">10.1038/16393</a></li><li>Schuler et al. (2013) <em>Optimal Hubbard Models for Materials with Nonlocal Coulomb Interactions: Graphene, Silicene, and Benzene</em> <a href="https://doi.org/10.1103/PhysRevLett.111.036601" target="_blank">10.1103/PhysRevLett.111.036601</a></li><li>Meng et al. (2010) <em>Quantum spin liquid emerging in two-dimensional correlated Dirac fermions</em> <a href="https://doi.org/10.1038/nature08942" target="_blank">10.1038/nature08942</a></li><li>Osofsky et al. (2016) <em>Functionalized graphene as a model system for the two-dimensional metal-insulator transition</em> <a href="https://doi.org/10.1038/srep19939" target="_blank">10.1038/srep19939</a></li><li>Igor (2006) <em>Herbut.: Interactions and Phase Transitions on Graphene’s Honeycomb Lattice</em> <a href="https://doi.org/10.1103/PhysRevLett.97.146401" target="_blank">10.1103/PhysRevLett.97.146401</a></li><li>Drut and Lähde (2009) <em>Is Graphene in Vacuum an Insulator?</em> <a href="https://doi.org/10.1103/PhysRevLett.102.026802" target="_blank">10.1103/PhysRevLett.102.026802</a></li><li>Castro Neto and Guinea (2009) <em>The electronic properties of graphene</em> <a href="https://doi.org/10.1103/RevModPhys.81.109" target="_blank">10.1103/RevModPhys.81.109</a></li><li>Das Sarma and Shaffique (2011) <em>Electronic transport in two-dimensional graphene</em> <a href="https://doi.org/10.1103/RevModPhys.83.407" target="_blank">10.1103/RevModPhys.83.407</a></li><li>Parr et al. (1950) <em>Molecular orbstal calculations of the lower excited electronic levels of benzene, configuration interaction included</em> <a href="https://doi.org/10.1063/1.1747540" target="_blank">10.1063/1.1747540</a></li><li>Wang et al. (2015) <em>Visualizing chemical states and defects induced magnetism of graphene oxide by spatially-resolved-X-ray microscopy and spectroscopy</em> <a href="https://doi.org/10.1038/srep15439" target="_blank">10.1038/srep15439</a></li><li>Midtvedt and Alexander (2016) <em>Strain–displacement relations for strain engineering in single-layer 2d materials</em> <a href="https://doi.org/10.1088/0953-8984/28/4/045302" target="_blank">10.1088/0953-8984/28/4/045302</a></li><li>Sahu and Rout (2014) <em>Study of Band Gap Opening in Graphene by Impurity and Substrate-Mediated Interactions</em> <a href="https://doi.org/10.1166/asl.2014.5374" target="_blank">10.1166/asl.2014.5374</a></li><li>Sahu and Rout (2015) <em>Model study of the effect of coulomb interaction on band gap of graphene-on-substrates</em> <a href="https://doi.org/10.1016/j.physb.2014.12.014" target="_blank">10.1016/j.physb.2014.12.014</a></li><li>Sahu et al. (2016) <em>Model study of band gap opening in graphene by electron–electron and electron–phonon interaction in high frequency range</em> <a href="https://doi.org/10.1166/asl.2016.6894" target="_blank">10.1166/asl.2016.6894</a></li><li>Sahu et al. (2016) <em>Theoretical study of band gap opening in AB-stacked bi-layer graphene by impurity and electric field effects</em> <a href="https://doi.org/10.1016/j.matpr.2016.01.119" target="_blank">10.1016/j.matpr.2016.01.119</a></li><li>Zubarev (1960) <em>Double-time green functions in statistical physics</em> <a href="https://doi.org/10.1070/PU1960v003n03ABEH003275" target="_blank">10.1070/PU1960v003n03ABEH003275</a></li><li>Sahu and Rout (2016) <em>Tight-binding model study of substrate induced pseudo-spin polarization and magnetism in mono-layer graphene</em> <a href="https://doi.org/10.1016/j.jmmm.2016.01.062" target="_blank">10.1016/j.jmmm.2016.01.062</a></li><li>Rout et al. (2016) <em>Effect of impurity doping on tunnelling conductance in AB-stacked bi-layer graphene: A tight-binding study</em> <a href="https://doi.org/10.1063/1.4945245" target="_blank">10.1063/1.4945245</a></li><li>Sahu et al. (2016) <em>Microscopic theoretical model study of band gap opening in AA-stacked bi-layer graphene</em> <a href="https://doi.org/10.1063/1.4946087" target="_blank">10.1063/1.4946087</a></li></ol>