First-principles theoretical investigation of graphene layers for sensor applications: A review

Transcription

1 Special Collection: Graphene First-principles theoretical investigation of graphene layers for sensor applications: A review Nanomaterials and Nanotechnology Volume 7: 1 7 ª The Author(s) 2017 DOI: / journals.sagepub.com/home/nax Yoshitaka Fujimoto Abstract Graphene is expected to be a potential device material for sensor applications due to its high charge mobility and high sensitivity to adsorbates. This article reviews the first-principles density-functional study that clarifies gas adsorption effects on graphene layers doped with boron and nitrogen atoms. We show adsorption effects of not only common gas molecules but also environmentally polluting or toxic gas molecules on stabilities and structural properties of graphene layers and carbon nanotubes. We also show physical properties induced by the adsorption of the gas molecules and discuss the possibility to detect these gas molecules. Keywords Graphene, carbon nanotubes, doping, gas adsorption, sensor application, first-principles density-functional calculation Date received: 17 May 2017; accepted: 21 July 2017 Topic: Low Dimensional Semiconductor Nanostructures Topic Editor: Nuno Silvestre Associate Editor: Nuno Silvestre Introduction Graphene, a two-dimensional carbon (C) atomic layer with a hexagonal network, exhibits various unique and excellent properties. 1 7 Among them, one of the most remarkable properties is that graphene shows extremely high charge mobility. 3 Furthermore, graphene as well as carbon nanotubes (CNTs) is also sensitive to adsorbates For example, graphene and CNTs can detect nitrogen dioxide (NO 2 ) and ammonia (NH 3 ) gases with fast response time and high sensitivity at room temperatures. 8,9 Thus, graphene is a good candidate for device materials to be used in highperformance sensor applications due to the high carrier mobility as well as the high sensitivity to adsorbates. One of the most effective methods to tune electronic properties of graphene and CNTs is to dope with heteroatoms. Moreover, substitutional doping with heteroatoms can often enhance the reactivity of nanomaterials. The boron (B) and nitrogen (N) atoms are good dopants for C-based nanomaterials. There are so far a large number of reports on B- and N-doped graphenes and CNTs, since B, C, and N elements possess similar structural properties The B-doped graphene and CNTs have been successfully synthesized through chemical vapor deposition and electrothermal reactions The p-type field-effect transistors are realized using the B-doped graphene. 26 It is reported theoretically that CO molecules are adsorbed on B-doped (8,0) nanotubes, while those are not adsorbed on the undoped ones. 27 On the other hand, experimental observations by X-ray photoelectron spectroscopy revealed two major N-defect Department of Physics, Tokyo Institute of Technology, Meguro, Tokyo, Japan Corresponding author: Yoshitaka Fujimoto, Department of Physics, Tokyo Institute of Technology, , Oh-Okayama, Meguro, Tokyo , Japan. fujimoto@stat.phys.titech.ac.jp Creative Commons CC BY: This article is distributed under the terms of the Creative Commons Attribution 4.0 License ( which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages ( open-access-at-sage).

2 2 Nanomaterials and Nanotechnology Table 1. Adsorption energy E a (ev) and distance d (Å) between molecule and B(N) atom for gas molecules adsorbed on B(N)-doped MLG and BLG and CNTs. MLG BLG 39 CNT 27,41 B 40 N 40 PN 34,36 B N B N CO E a d CO 2 E a d NO E a d NO 2 E a d NH 3 E a d H 2 O E a d H 2 E a , 4.16 d CNT: carbon nanotube; B: boron; N: nitrogen; PN: pyridine-type N defects; MLG: monolayer graphene; BLG: bilayer graphene. bonding configurations in N-doped graphene as well as N-doped CNTs; one is the substitutional N defect where the C atom is simply replaced by the N atom, and the other is the pyridine-type N defects where the N atoms are located around atomic vacancies. 12,22,28,29 Theoretical calculations suggested that the structural variety of N defects in graphene and CNTs provides various electronic properties such as p-type and n-type conduction properties The N doping into graphene as well as CNTs improves the adsorption energies of hydrogen atoms. 33,34 It is also reported that the N-doped graphene with transition metals might provide hydrogen storages. 35 Furthermore, the N-dopant sites in the N-doped graphene behave as chemical reaction centers. 36 Thus, the graphene doped with B and N atoms not only offers rich electronic properties but also is expected to be useful device materials in sensor applications. In this article, we review the first-principles densityfunctional study of the adsorptions of environmentally polluting and common gas molecules on B-doped and N-doped graphene layers. The section on Structural properties and stabilities reports the adsorption effects of various gas molecules on the stabilities and the structural properties of B- and N-doped graphene layers and CNTs. The next section on Physical properties discusses several physical properties for detecting gas molecules. Finally, we summarize this article together with outlook in the Concluding remarks section. Structural properties and stabilities In this section, the adsorption effects on the stabilities and the structural properties of the environmentally polluting and common gas molecules (CO, CO 2, NO, NO 2,H 2,H 2 O, and NH 3 ) to the B- and N-doped graphene layers and CNTs are reported. Energetics Table 1 summarizes the adsorption energies (E a ) and the distances (d) between molecules and dopant atoms for the adsorption of various molecules on the B(N)-doped graphenes and CNTs. For adsorptions of CO and CO 2 molecules, the adsorption energies of B(N)-doped graphene monolayers as well as bilayers take relatively small values and the distances between the molecules and the dopant atoms are considerably long (d > 2.5 Å). On the other hand, in the case of the B-doped (8,0) CNT, the adsorption energy for the CO molecule (je a j¼approximately 0.8 ev) is much larger than the adsorption energies of the B(N)-doped graphene monolayers and bilayers, and the distance between the CO molecule and the B atom in the (8,0) CNT is relatively short (d ¼ approximately 1.5 Å). 27 Thus, the CO and CO 2 molecules are not chemically but rather physically adsorbed on the doped graphene monolayers and bilayers, while the B-doped (8,0) CNT can bind chemically with the CO molecule. It should be noticed that the adsorption energies for the monolayer graphenes are overall smaller than those of the bilayer graphenes. This is because the adsorption energies for the bilayer graphene and the CNTs with CO and H 2 O molecules are obtained within the localdensity approximation, 37 whereas those for the monolayer graphenes and the CNTs with NH 3 and NO 2 molecules are within the generalized gradient approximation, except for the pyridine-type defects. 38 In the case of NO and NO 2 molecule adsorptions, the adsorption energies for the B-doped graphene monolayers and bilayers are relatively larger compared with those for the CO and CO 2 molecules on the B(N)-doped graphene monolayers and bilayers. The distances between the NO (NO 2 ) molecule and the B-doped graphene layers are relatively short. 39,40 Therefore, the NO and NO 2 molecules can bind chemically on the B-doped graphene monolayers and

3 Fujimoto 3 Table 2. Charge transfer between molecules and B-doped MLG, BLG graphenes, and B(N)-doped CNTs. a MLG 40 BLG 39 CNT 41 B B B N NO 0.03 NO NH CNT: carbon nanotube; B: boron; N: nitrogen; MLG: monolayer graphene; BLG: bilayer graphene. a The positive and negative values mean electron transfers from molecules to doped graphene layers (CNTs) and from doped graphene layers (CNTs) to molecules, respectively. Figure 1. Optimized atomic geometries of the B-doped bilayer graphene with (a) NO molecule, (b) NO 2 molecule, (c) the trimerized, and (d) the tetramerized pyridine-type defected graphenes with H 2 molecules, and the pyridine-type defected graphene with (e) NH 3 and (f) H 2 O molecules. 36 Reproduced with permission from, 39 copyright 2016 the Elsevier, and from, 34 copyright 2014 the American Institute of Physics. bilayers. For the N-doped cases, these two molecules are physically adsorbed with the relatively small adsorption energies as well as the very long distances between the molecule and the dopant atom (d >2.6Å). On the other hand, the NO 2 molecule is chemically adsorbed on the N-doped (10,0) CNT. 41,42 It is interesting that the adsorption energies are enhanced with the curvature effects of nanotubes. The H 2,NH 3,andH 2 O molecules are physically adsorbed on the B-doped as well as the N-doped monolayer graphenes with small adsorption energies and long distances between the molecules and the dopant atoms. 40 In the case of the B-doped CNTs, NH 3 and H 2 O molecules can be adsorbed chemically, while, for the N-doped CNTs, these molecules are physically adsorbed with the small adsorption energies and the long distances. 33,41 The H 2, NH 3,andH 2 O molecules can chemically bind to the pyridine-type N defects in the monolayer graphenes with very large adsorption energies. 34,36 Structure The chemical adsorption of the gas molecules would give rise to the structural deformation of graphene layers. When the NO 2 molecule is adsorbed on the B-doped graphene bilayer, the B atom in the graphene moves outward from the planar sheet (Figure 1(b)). Thereby, the binding distance between the B atom and the O atom in the NO 2 molecule becomes very short. The structural deformation of the B atom is also reported in the B-doped monolayer graphene. 40 On the other hand, in the case of the NO molecule adsorption, the B atom still resides in the planar sheet (Figure 1(a)). 39 The short binding distance for the NO 2 adsorption gives rise to the strong chemical bonds between the NO 2 molecule and the B-doped graphene layer compared with the NO molecule adsorption. For the H 2 adsorptions on the pyridine-type defects, the H 2 molecule is dissociatively adsorbed on the trimerized as well as the tetramerized pyridine-type defects in the graphene (Figure 1(c) and (d)). In both the cases, the chemical dissociation of the H 2 molecule takes place without a reaction barrier, that is, the H 2 molecule near the pyridine-type defect can dissociate into two H atoms and eventually bind chemically with two different N atoms in the pyridine-type defect. 34 In the case of adsorptions of the NH 3 molecule on the pyridine-type defects, the NH 3 molecule dissociates into the H atom and the NH 2 molecule. For the H 2 O molecule adsorption, the dissociative adsorption takes place into the H atom and the OH molecule as well. On the other hand, for the adsorption of these three molecules on the substitutional N defects, the dissociative adsorption does not occur. 40 Unlike the substitutional N defect, the pyridine-type defect exhibits the high reactivity for the H 2,NH 3, and H 2 O molecules. Therefore, graphenes with the pyridine-type defects would be useful as not only sensors but also catalysts. Physical properties This section reports several physical properties arising from the chemical adsorption of the gas molecules. Charge transfer When the molecules are adsorbed chemically to the doped graphene layers and CNTs, the charge transfers are induced between the adsorbed molecules and the graphene layers (CNTs) due to the orbital hybridizations. Table 2 lists the charge transfers between adsorbed molecules and doped graphene layers (CNTs). By the NO molecule adsorption, the electrons move from the NO molecule to the B-doped

4 4 Nanomaterials and Nanotechnology bilayer graphene with approximately 0.03 electrons. 39 On the other hand, in the case of the NO 2 molecule adsorptions, it is interesting that the electrons also transfer from the B-doped graphene monolayer as well as bilayer to the NO 2 molecule. For the B-doped and N-doped CNTs, the electrons also move from doped CNTs to the NO 2 molecule. 41 For the adsorption of the NH 3 molecule, electron transfers occur from the NH 3 molecule to the B-doped CNT. However, in the case of the N-doped CNT, the electron transfer does not occur between the NH 3 molecule and the N-doped CNT, because the NH 3 molecule is not chemically but physically adsorbed to the N-doped CNT with the long binding distance (see Table 1). Thus, the NO and NH 3 molecules behave as electron donors, whereas the NO 2 molecule acts as electron acceptors. Work function The adsorption of the molecules would change the work functions of graphene layers. The work functions of the pristine graphene monolayer and bilayer are approximately 4.5 ev, which agree with experimentally observed results When the B atom is doped into bilayer graphene, the work function increases to approximately 5 ev, since the number of the valence electron of the B atom is less by one compared to the C atom. On the other hand, the work function of the N-doped bilayer graphene decreases below 4 ev, since the N atom has one extra electron than the C atom. When the NO 2 molecule is adsorbed to the B-doped bilayer graphene, its work function further increases up to beyond 5 ev (Figure 2(a)). 39 In the case of the adsorption of the NO 2 molecule, the work function diminishes below that of the pristine graphene. The variation in the work function induced by the adsorption of different types of adsorbates is mainly caused by the charge transfer. In the case of the NO 2 molecule, electrons move from the B-doped graphene layer into the NO 2 molecule. On the other hand, in the case of the NO molecule, there exists the electron transfer from the NO molecule into the B-doped layer. The charge transfer between the molecule and the graphene layers leads to the shift in the Fermi energy and should modify the electronic structures and the transport properties of the graphene layers. 41,46 51 The monolayer graphene with pyridine-type N defects exhibits unusual behaviors in work functions. 34 When the N atom is substitutionally doped to the graphene (SN), the work function decreases as discussed above. On the other hand, when the trimerized (PN 3 ) as well as the tetramerized pyridine-type (PN 4 ) defect is formed in the graphene, 52 it is interesting that the work functions increase, since the pyridine-type defects in the graphene could induce the acceptor states (see Figure 2(b)). When the H 2 molecules are adsorbed to the trimerized as well as the tetramerized pyridine-type defect in the graphene, their work functions decrease dramatically below that of the pristine graphene (H-PN 3 and H-PN 4 ). The work function varies depending Figure 2. Work functions (WFs) of (a) doped bilayer graphenes with and without NO and NO 2 molecules and (b) various N defects in monolayer graphene with and without the H 2 molecules. The horizontal lines denote WFs of (a) the pristine bilayer graphene and (b) the pristine monolayer graphene. Reproduced with permission from, 39 copyright 2016 the Elsevier and from, 34 copyright 2014 the American Institute of Physics. on the p-type and n-type conduction properties. 34,53,54 In addition, the work function depends on the dopant concentration. Scanning tunneling microscopy Scanning tunneling microscopy (STM) is a powerful tool to observe the surface electronic structures of semiconductors and metals at atomic levels Hence, not only the types of adsorbates but also the atomic configurations of adsorbed molecules could be identified using the STM methods. When the B atom is doped into the graphene layers, a triangle-shaped hillock is observed around the B atom in the STM image of the B-doped bilayer graphene (Figure 3(a)), 53,54 and it is also observed experimentally. 62 When the NO molecule is adsorbed, two ring-like corrugations appear above the NO molecule (Figure 3(b)). 39 For the case of the NO 2 molecule adsorption, the STM image has large and small spherical-shaped protrusions above the NO 2 molecule (Figure 3(c)). Thus, the B atom and the NO

5 Fujimoto 5 Figure 3. STM images of B-doped bilayer graphenes (a) without adsorbate, (b) with NO molecule, and (c) with NO 2 molecule. The STM images are generated at the bias voltage of þ0.5 ev. Reprinted with permission from, 39 copyright 2016 the Elsevier. STM: scanning tunneling microscopy; B: boron. molecule as well as the NO 2 molecule are observable using the STM methods. Concluding remarks This review reported the first-principles density-functional study that clarifies the gas adsorption effects on the atomic structures, the stabilities, and the physical properties of the graphenes doped with B and N atoms. The B-doped graphene can bind chemically with the NO and NO 2 molecules, while the CO and CO 2 molecules are not chemically adsorbed. The H 2,NH 3, and H 2 O molecules are not chemically bound to the substitutionally N-doped and B-doped graphenes as well, whereas those molecules can bind chemically with the pyridine-type N defective graphenes. The formation of chemical bonds between the adsorbed molecules and the graphenes produces the charge transfers, and the electronic structures of the graphenes are modified significantly. The work function varies depending on the types of the adsorbates, and therefore, the adsorbed molecules on the graphenes could be detected by gauging the work function. The STM method could identify the types of adsorbates as well as the atomic configurations of adsorbed molecules. From the viewpoint of the practical sensor applications, it is one of the effective ways to examine the variation in the electrical conductance of the graphenes induced by the adsorption of the gas molecules and it is an important issue of future work. The theoretical analysis based on the first-principles electronic-structure calculations suggests that the B-doped graphenes are candidate materials for detecting the NO and NO 2 molecules, and the pyridine-type N defects in the graphenes are also important for detecting various common gas molecules. Furthermore, the pyridine-type N defects would be useful as catalysts and chemical reaction centers. Author note Computations were partly done at the Institute for Solid State Physics, the University of Tokyo, and at the Cybermedia Center of Osaka University. Declaration of conflicting interests The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article. Funding The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was financially supported by MEXT Elements Strategy Initiative to Form Core Research Center through the Tokodai Institute for Element Strategy, JSPS KAKENHI grant number JP17K05053 and grant number JP References 1. Novoselov KS, Geim AK, Morozov SV, et al. Electric field effect in atomically thin carbon films. Science 2004; 306: Berger C, Song Z, Li X, et al. Electronic confinement and coherence in patterned epitaxial graphene. Science 2006; 312: Castro EV, Novoselov KS, Morozov SV, et al. Biased bilayer graphene: semiconductor with a gap tunable by the electric field effect. Phys Rev Lett 2007; 99: Zhang Y, Tan YW, Stormer HL, et al. Experimental observation of the quantum Hall effect and Berry s phase in graphene. Nature (London) 2005; 438: Zhang Y, Tang T, Girit C, et al. Direct observation of a widely tunable bandgap in bilayer graphene. Nature (London) 2009; 459: Williams JR, Dicalro L, and Marcus CM. Quantum hall effect in a gate-controlled p-n junction of graphene. Science 2007; 317: Young AF and Kim P. Quantum interference and Klein tunnelling in graphene heterojunctions. Nat Phys 2009; 5: Schedin F, Geim AK, Morozov SV, et al. Detection of individual gas molecules adsorbed on graphene. Nat Mater 2007; 6: Kong J, Franklin NR, Zhou C, et al. Nanotube molecular wires as chemical sensors. Science 2000; 287: Collins PG, Bradley K, Ishigami M, et al. Extreme oxygen sensitivity of electronic properties of carbon nanotubes. Science 2000; 287:

6 6 Nanomaterials and Nanotechnology 11. Goldoni A, Lareiprete R, Petaccia L, et al. Single-wall carbon nanotube interaction with gases: sample contaminants and environmental monitoring. JAmChemSoc2003; 125: Wei D, Liu Y, Wang Y, et al. Synthesis of N-doped graphene by chemical vapor deposition and its electrical properties. Nano Lett 2009; 9(5): Kang HS and Jeong S. Nitrogen doping and chirality of carbon nanotubes. Phys Rev B 2004; 70: Fujimoto Y and Saito S. Effects of strain on carbon donors and acceptors in hexagonal boron nitride monolayers. Phys Rev B 2016; 93: Fujimoto Y and Saito S. Band engineering and relative stabilities of hexagonal boron nitride bilayers under biaxial strain. Phys Rev B 2016; 94: Wang Y, Shao Y, Matson DW, et al. Nitrogen-doped graphene and its application in electrochemical biosensing. ACS Nano 2010; 4: Fujimoto Y and Saito S. Atomic geometries and electronic structures of hexagonal boron-nitride bilayers under strain. J Ceram Soc Jpn 2015; 123: Ayala P, Arenal R, Loiseau A, et al. The physical and chemical properties of heteronanotubes. Rev Mod Phys 2010; 82: Fujimoto Y, Koretsune T, and Saito S. Electronic structures of hexagonal boron-nitride monolayer: strain-induced effects. J Ceram Soc Jpn 2014; 122: Fujimoto Y and Saito S. Interlayer distances and band-gap tuning of hexagonal boron-nitride bilayers. J Ceram Soc Jpn 2016; 124: Cho YJ, Kim HS, Im H, et al. Nitrogen-doped graphitic layers deposited on silicon nanowires for efficient lithium-ion battery anodes. J Phys Chem C 2011; 115: Susi T, Kotakoski J, Arenal R, et al. Atomistic description of electron beam damage in nitrogen-doped graphene and single-walled carbon nanotubes. ACS Nano 2012; 6(10): Gopalakrishnan K, Govindaraj A, and Rao CNR. Extraordinary supercapacitor performance of heavily nitrogenated graphene oxide obtained by microwave synthesis. J Mater Chem A 2013; 1(26): Cattelan M, Agnoli S, Favaro M, et al. Microscopic view on a chemical vapor deposition route to boron-doped graphene nanostructures. Chem Mater 2013; 25: Gebhardt J, Koch RJ, Zhao W, et al. Growth and electronic structure of boron-doped graphene. Phys Rev B 2013; 87: Tang YB, Yin LC, Yang Y, et al. Tunable band gaps and p- type transport properties of boron-doped graphenes by controllable ion doping using reactive microwave plasma. ACS Nano 2012; 6(3): Peng S and Cho K. Ab initio study of doped carbon nanotube sensors. Nano Lett 2003; 3(4): Czerw R, Terrones M, Charlier JC, et al. Identification of electron donor states in N-doped carbon nanotubes. Nano Lett 2001; 1(9): Ayala P, Grüneis A, Gemming T, et al. Tailoring N-doped single and double wall carbon nanotubes from a nondiluted carbon/nitrogen feedstock. JPhysChemC2007; 111(7): Fujimoto Y and Saito S. Formation, stabilities, and electronic properties of nitrogen defects in graphene. Phys Rev B 2011; 84: Fujimoto Y and Saito S. Energetics and electronic structures of pyridine-type defects in nitrogen-doped carbon nanotubes. Physica E 2011; 43: Fujimoto Y. Formation, energetics, and electronic properties of graphene monolayer and bilayer doped with heteroatoms. Adv Condens Mat Phys 2015; 2015: Fujimoto Y and Saito S. Structure and stability of hydrogen atom adsorbed on nitrogen-doped carbon nanotubes. J Phys Conf Ser 2011; 302: Fujimoto Y and Saito S. Hydrogen adsorption and anomalous electronic properties of nitrogen-doped graphene. JAppl Phys 2014; 115: Kim G, Jhi SH, and Park N. Effective metal dispersion in pyridinelike nitrogen doped graphenes for hydrogen storage. Appl Phys Lett 2008; 92: Fujimoto Y and Saito S. Adsorption of molecules on nitrogen-doped graphene: a first-principles study. JPS Conf Proc 2015; 4: Perdew JP and Zunger A. Self-interaction correction to density-functional approximations for many-electron systems. Phys Rev B 1981; 23: Perdew JP, Burke K, and Ernzerhof M. Generalized gradient approximation made simple. Phys Rev Lett 1996; 77: Fujimoto Y and Saito S. Gas adsorption, energetics and electronic properties of boron- and nitrogen-doped bilayer graphenes. Chem Phys 2016; 478: Dai J, Yuan J, and Giannozzi P. Gas adsorption on graphene doped with B, N, Al, and S: a theoretical study. Appl Phys Lett 2009; 95: Bai L and Zhou Z. Computational study of B- or N-doped single-walled carbon nanotubes as NH 3 and NO 2 sensors. Carbon 2007; 45: Tang Q, Zhou Z, and Chen Z. Graphene-related nanomaterials: tuning properties by functionalization. Nanoscale 2013; 5: Kim H, Renault O, Tyurnina A, et al. Doping efficiency of single and randomly stacked bilayer graphene by iodine adsorption. Appl Phys Lett 2014; 105: Yu YJ, Zhao Y, Ryu S, et al. Tuning the graphene work function by electric field effect. Nano Lett 2009; 9: Yan R, Zhang Q, Li W, et al. Determination of graphene work function and graphene-insulator-semiconductor band alignment by internal photoemission spectroscopy. Appl Phys Lett 2012; 101: Fujimoto Y and Hirose K. First-principles treatments of electron transport properties for nanoscale junctions. Phys Rev B 2003; 67:

7 Fujimoto Fujimoto Y, Hirose K, and Ohno T. Calculations of surface electronic structures by the overbridging boundary-matching method. Surf Sci 2005; 586: Fujimoto Y, Asari Y, Kondo H, et al. First-principles study of transport properties of Al wires: comparison between crystalline and jellium electrodes. Phys Rev B 2005; 72: Fujimoto Y and Oshiyama A. Formation and stability of 90 degree dislocation cores in ge films on si(001). AIP Conf Proc 2011; 1399: Fujimoto Y and Hirose K. First-principles calculation method of electron-transport properties of metallic nanowires. Nanotechnology 2003; 14: Travlou NA, Ushay C, Seredych M, et al. Nitrogen-doped activated carbon-based ammonia sensors: effect of specific surface functional groups on carbon electronic properties. ACS Sens 2016; 1: Lin YC, Teng PY, Yeh CH, et al. Structural and chemical dynamics of pyridinic-nitrogen defects in graphene. Nano Lett 2015; 15: Fujimoto Y and Saito S. Energetics and scanning tunneling microscopy images of B and N defects in graphene bilayer. Springer Proc in Phys 2017; 186: Fujimoto Y and Saito S. Electronic structures and stabilities of bilayer graphene doped with boron and nitrogen. Surf Sci 2015; 634: Tersoff J and Hamann DR. Theory and application for the scanning tunneling microscope. Phys Rev Lett 1983; 50: Tersoff J and Hamann DR. Theory of scanning tunneling microscopy. Phys Rev B 1985; 31: Fujimoto Y, Okada H, Endo K, et al. Images of scanning tunneling microscopy on the Si(001)-p(22) reconstructed surface. Mater Trans 2001; 42: Okada H, Fujimoto Y, Endo K, et al. Detailed analysis of scanning tunneling microscopy images of the Si(001) reconstructed surface with buckled dimers. Phys Rev B 2001; 63: Fujimoto Y and Oshiyama A. Atomic structures and energetics of 90 dislocation cores in Ge films on Si (001). Phys Rev B 2010; 81: Fujimoto Y, Okada H, Inagaki K, et al. Theoretical study on the scanning tunneling microscopy image of Cl-adsorbed Si(001). Jpn J Appl Phys 2003; 42: Fujimoto Y and Oshiyama A. Structural stability and scanning tunneling microscopy images of strained Ge films on Si(001). Phys Rev B 2013; 87: Zhao L, Levendorf M, Goncher S, et al. Local atomic and electronic structure of boron chemical doping in monolayer graphene. Nano Lett 2013; 13(10):