B–N@Graphene: Highly Sensitive and Selective Gas Sensor - The

We have performed density functional theory (DFT) calculations to study the gas (CO, CO2, NO, and NO2) sensing mechanism of pure and doped (B@, N@, an...
0 downloads 0 Views 6MB Size
Article pubs.acs.org/JPCC

B−N@Graphene: Highly Sensitive and Selective Gas Sensor Indrani Choudhuri,† Nandini Patra,‡,§ Arup Mahata,† Rajeev Ahuja,‡,∥ and Biswarup Pathak*,†,§ †

Discipline of Chemistry, School of Basic Sciences, Indian Institute of Technology (IIT) Indore, 452020 M.P., India Condensed Matter Theory Group, Department of Physics and Astronomy, Uppsala University, 751 20 Uppsala, Sweden § Centre for Material Science and Engineering, Indian Institute of Technology (IIT) Indore, 452020 M.P., India ∥ Applied Materials Physics, Department of Material Science and Engineering, Royal Institute of Technology (KTH), 100 44 Stockholm, Sweden ‡

Downloaded via UNIV OF SOUTH DAKOTA on June 30, 2018 at 16:07:13 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.

S Supporting Information *

ABSTRACT: We have performed density functional theory (DFT) calculations to study the gas (CO, CO2, NO, and NO2) sensing mechanism of pure and doped (B@, N@, and B−N@) graphene surfaces. The calculated adsorption energies of the various toxic gases (CO, CO2, NO, and NO2) on the pure and doped graphene surfaces show, doping improves adsorption energy and selectivity. The electronic properties of the B−N@graphene surfaces change significantly compared to pure and B@ and N@graphene surfaces, while selective gas molecules are adsorbed. So, we report B−N codoping on graphene can be highly sensitive and selective for semiconductor-based gas sensor.

1. INTRODUCTION Gas sensors are very much important to detect the presence of various toxic gases, which are very much harmful to organic life.1 Gaseous pollutants, which release in a large amount from the industries and automobiles, spread in our environments very rapidly and affect the human health, ozone layer, and global climates. Semiconductor based gas sensors are highly demanding because of their exceptional stability and surface reactivity in the ambient conditions. The primary challenge of the semiconductor based gas sensor is to meet the demand for higher sensitivity and selectivity so that even the trace amounts of toxic gases can be detected. The metal oxide based semiconductors are quite efficient,2−5 but they are not environment friendly.6 So, the two-dimensional metal free semiconducting material with excellent surface, optical, and electrical properties could be the most promising gas sensing material. Graphene is the thinnest possible material,7−9 with the largest surface area (3000 m2 per gram) and got remarkable strength and stiffness, which can be stretched up to 20%, elastically. Theoretically, monolayer of graphene and its electronic properties were proposed as early as in 1947.10,11 However, 40 years later the material is synthesized and showed similar electronic properties as proposed theoretically.9,12 Thus, it clearly shows the importance of theoretical studies in the prediction of new structures and their properties. There is a number of researchers that believes graphene can be an excellent sensor owing to its two-dimensional structure and excellent electrical, optical, and mechanical properties.13−15 Recently, Schedin et al.16 showed graphene is a unique and © 2015 American Chemical Society

attractive sensing material to detect all the individual events when some gaseous molecules attaches or detaches from its surface. Even its local carrier concentration changes when gaseous molecules are adsorbed, which leads to changes in resistance. Graphene is an exceptionally low noise material thus very sensible toward the charge fluctuations and so a suitable candidate for chemical detection.17,18 Graphene based devices are reported to be a promising candidate for DNA sequencing19 as DNA can interact with graphene through π−π stacking, which diminishes the ballistic conductance of graphene based materials.20−22 Despite of all, the zero band gap and absence of dangling bonds in graphene limits its device-based applications. Various approaches have been taken to fabricate the high-performance graphene devices by engineering their band gaps for improving their semiconducting properties.23 Doping is an efficient approach, which tunes the electronic properties of graphene. Substitutional doping of graphene with different atoms (e.g., B, N, B−N) affects the sp2 hybridization of carbon atoms.24 The reason for choosing the N and B atoms to replace the C atoms is the atomic masses of these dopants are closer to carbon atom. The electron rich and deficient nature of N and B atoms helps to alter significantly the electronic properties of the host material. The researches of Lherbie et al.25 studied the charge mobility and conductivity of the system by doping graphene with different concentrations of B and N impurities. CervantesReceived: July 30, 2015 Revised: October 16, 2015 Published: October 16, 2015 24827

DOI: 10.1021/acs.jpcc.5b07359 J. Phys. Chem. C 2015, 119, 24827−24836

Article

The Journal of Physical Chemistry C Sodi et al.26 theoretically predicted the band gap opening in graphene by elemental (B and N) doping. Since then, there are many experimental evidence of band gap opening in graphene by B/N doping.27−29 Panchakarla et al.27 experimentally synthesized B/N doped graphene and confirmed the Fermi energy shifting (p/n doping) as theoretically proposed. Fan et al.30 theoretically predicted band gap opening of graphene by doping small boron-nitride domains. Later Bepete et al. experimentally31 confirmed such predictions. So, lately, theoretical studies becoming an essential tool for predicting new materials and their properties. There are several theoretical and experimental studies on adsorption of gaseous molecules on pure and doped graphene surfaces.32−36 Ko et al.32 recently demonstrated that graphene can selectively sense NOx molecules at room temperature. The sensitivity was up to 9% in the presence of 100 ppm of NO2 gas. So, theoretical calculations are needed to show the changes in electronic properties of graphene in presence of gases. Zhang et al.33 showed how B/N doped and defected graphene interacts with various gas molecules. They have concluded that the gas sensing property of graphene can be drastically improved by introducing dopants and defects. Dai et al.34 theoretically reported how elemental (B/N/Al and S) doping on graphene improves the gas sensing properties. Such calculations are helpful in predicting the changes in electrical properties of graphene, which in turn are helpful for gas sensing. However, most of the previous research mainly focused on the perfect or mono doped graphene32−34 or metal decorated35 graphene for gas sensing properties. Those studies rather predict the selectivity on the basis of adsorption energies, and incorporation of metals limits their applicability. However, it is experimentally proven that codoped graphene has many excellent applications toward sensor,36−38 photocatalysis,39 ORR,40,41 and supercapacitor. 42 Recently, Niu et al.36 experimentally reported that nitrogen and silica codoped graphene nanosheet is a very good sensor for NO2. However, whether such codoping could improve the sensitivity and selectivity of the gas sensing properties is not well studied theoretically. In reality, the graphene sheets, which are prepared by the available fabrication methods, have many defects. Deliberately or accidentally graphene could be doped by non-carbon elements. So, our primary aim is to codope graphene by environmental friendly nonmetal (like B or N) elements using the state of the art theoretical techniques. This is done to improve the selectivity and sensitivity of the gas sensor, as such doping not only improves the interaction between the graphene sheet and the gaseous molecule but also changes their electronic properties.

correction46 as the long-range interaction present between the surface and the gaseous molecules. The Brillion zone is integrated using Monkhorst−Pack generated sets of 11 × 11 × 1 k-points.47 For the density of states (DOS) calculations, the Monkhorst−Pack generated 45 × 45 × 1 sets of k-points are used with a Gaussian smearing of 0.001 eV. The adsorption energy (EA) of the toxic greenhouse gases (CO, CO2, NO, and NO2) on the pure and doped graphene surfaces is calculated using the following equation: EA = (Egraphene + Egas) − Egas − graphene

(1)

Here Egas‑graphene represents the energy of the optimized structure of the gaseous molecules adsorbed on the pure graphene surfaces. The Egraphene and Egas represent the singlepoint energy of the pure graphene and gas molecules, respectively. The single-point energy is calculated using the optimized coordinates of the respective systems. The charge density difference (CDD)48 is plotted to understand the nature of bonding between the pure/doped graphene and the gas molecules. The CDD is calculated using the following equation: ρCDD = ρ total −

∑ ρifragments i

(2)

Here the ρtotal is the total charge density of the system and ρifragments is the charge density of the individual fragments by which the system is made of. Here the charge density of the fragments (ρifragments) is calculated by a pseudo structure in which fragment part retains the same structure as in the total system but other parts are deleted. Here, we have mainly two different fragments, the pure graphene sheet and adsorbed gases. In CDD the positive and negative charge density is plotted by blue and yellow colors, respectively.

3. RESULTS AND DISCUSSION 3.1. Pure Graphene for Gas Sensor. A two-atom unit cell is considered for the graphene hexagonal structure. The structure is fully relaxed and the carbon−carbon bond length (1.42 Å) of the optimized structure matches with the experimental value of 1.42 Å.49 A two-atom unit cell is extended to 4 × 4 supercell of 32 carbon atoms to study the gaseous molecules interacting with the graphene surface. The total charge density of the pure graphene is plotted, and we find the charges are highly delocalized over the graphene. Bader charge analysis is done50 using the Henkelman program51,52 with near-grid algorithm refine-edge method. It shows, all the carbon atoms of graphene have the same charges (∼4e), thus confirming the delocalization. The electronic structure of pure graphene is one of the most exciting findings of this century. The total density (TDOS) and partial density of states (PDOS) of graphene is plotted in Figure S1 in Supporting Information. It shows graphene has 0 eV band gap, and carbon 2p orbitals are appearing on both sides of the Fermi energy (0 to −6.0 eV and 0 to 6.0 eV) agrees well with the previous experimental and theoretical studies.53 The change in the semiconductor’s electronic properties due to the gaseous molecular adsorption is the most important criteria for semiconductor based gas sensor.54 Band theory suggests, when gaseous molecules are adsorbed, the valence band (VB) electrons interact with the gaseous molecules and thus tunes the band position. Such tuning controls the electron flow from VB to CB, and thus the current. So, the main

2. COMPUTATIONAL DETAILS We have used the Vienna ab initio simulation package (VASP)43 to do all the calculations. The exchange-correlation interaction is treated in the level of the GGA using the Perdew−Burke−Ernzerhof (GGA-PBE) exchange-correlation functional.44 The projected augmented wave (PAW) method45 is employed using an energy cutoff of 470 eV to describe the electronic wave function. Thus, generalized gradient approximation (GGA-PBE) within projector-augmented wave (PAW) methodology is adopted to investigate the electronic properties of such materials. In all the calculations, self-consistency is achieved with a tolerance in the total energy set to 10−4 eV. We have included semiempirical DFT-D3 type of dispersion 24828

DOI: 10.1021/acs.jpcc.5b07359 J. Phys. Chem. C 2015, 119, 24827−24836

Article

The Journal of Physical Chemistry C

Figure 1. Optimized structures of (a) CO, (b) CO2, (c) NO, and (d) NO2 molecules adsorbed on graphene surface (brown, blue, and red colored balls denote C, N, and O, respectively). Charge density plots (Isosurface value 0.2322 e·Å−3) of (e) CO, (f) CO2, (g) NO, and (h) NO2 molecules adsorbed on graphene surface. The CDD (charge density difference) plots of (i) CO, (j) CO2, (k) NO, and (l) NO2 molecules adsorbed on graphene surface (isosurface value = 0.0036 e·Å−3).

Figure 2. Total and partial density of states (TDOS/PDOS) are plotted for (a) CO, (b) CO2, (c) NO, and (d) NO2 molecules adsorbed on pure graphene surfaces (EF is set to zero and indicated by the black dashed line).

objective of this study is to find out the changes in the semiconductor’s electronic properties due to the adsorption of the gases. Now, the question is whether toxic/greenhouse gases (CO, CO2, NO, NO2) could be adsorbed on the pure graphene surfaces or not? The CO/CO2 and NO/NO2 molecules are adsorbed on the 32-atom unit cell of graphene through C-/Oand N-/O- center, respectively. The optimized structures of the adsorbed gaseous molecules (CO, CO2, NO, and NO2) are presented in Figure 1a−d. The CO/CO2 molecules adsorbed through C-center is thermodynamically 0.02/0.02 eV [Figure 1a,b] more stable than through the O-center. The distance between the graphene and the adsorbed CO molecule is 3.08 Å, thus interacts weakly. Similarly NO/NO2 molecules are more stable (by 0.01/0.01 eV) while adsorbed through N-center [Figure 1c,d]. Here again the interaction between the gaseous molecules (NO and NO2) and graphene surface is very weak as the distance between the surface and the gaseous molecules are more than 2.98 Å. The calculated adsorption energies for CO, CO2, NO, and NO2 molecules are 0.12, 0.11, 0.10, and 0.26 eV, respectively. Our results show NO2 interacts better among all the four gases, and this is fairly in agreement with previous experimental study, where they show graphene could selectively sense NO2.32 The total charge density and CDD [Figure 1e−l] plots confirm the electrostatic nature of the interaction between

the graphene and the gas molecules. Therefore, we can say the molecules are physisorbed on the graphene surface. Bader charges are calculated for the gas molecules adsorbed on the graphene surface. We find that there are no changes in the total Bader charges on the C1 as well on the neighboring carbon atoms (C2, C3, and C4) [Table S1, Supporting Information] when gaseous molecules are adsorbed on the graphene. In fact, the total number of valence electrons on the CO (10e), CO2 (16e), NO (11e), and NO2 (17e) molecules remains constant even after adsorption on the graphene surface. Thus, there are no charge transfer between the graphene and the adsorbed gas molecules. Our adsorption energy, charge density, and Bader charge analysis study shows that the interaction between the pure graphene and toxic gases are very weak. We have plotted the TDOS and PDOS of the various gases adsorbed on the graphene surface [Figure 2a−d]. The total density of states of graphene does not change much due to the adsorption of CO and CO2 molecules [Figure 2a,b] but gives a sharp peak at Fermi energy while NO/NO2 molecules [Figure 2c,d] are adsorbed. These sharp peaks are [Figure 2c,d] coming from the 2p orbitals of the N and O atoms of the NO and NO2. However, electronic properties did not change much due to the gas adsorption. Thus, we believe pure graphene cannot be a promising material for semiconductor based gas sensor. 24829

DOI: 10.1021/acs.jpcc.5b07359 J. Phys. Chem. C 2015, 119, 24827−24836

Article

The Journal of Physical Chemistry C The B-, N-, and B−N doped graphene is studied to improve the interaction with the gaseous molecules. Moreover, B-doped graphene (B@graphene) is a p-type semiconductor so hole is the major carrier. Now, when oxidizing gases (NO2, CO2) come in contact with the semiconductor it increases the number of holes in the VB. As a result, the conductivity of the semiconductor increases, but in the case of reducing gases (NO, CO) the effect is opposite. Similarly, N-doped graphene (N@graphene) is an n-type semiconductor. So, when it comes in contact with the reducing gases the number of electron increases in valence band. So the electrical conductivity increases. It also shows opposite effects for oxidizing gases. Therefore, such doping is very important for improving the sensitivity and selectivity of the gas sensors. 3.2. B@-, N@-, B−N@Graphene for Gas Sensor. 3.2.1. Formation Energy. One boron/nitrogen atom is doped in the 32 atomic unit cell of graphene, which corresponds to 3.12% of boron/nitrogen doping. The doping concentration is varied from 3.12% to 50.00% by doping 3 (9.38%), 5 (15.60%), 7 (21.80%), and 16 (50.00%) boron/ nitrogen atoms in the 32 atomic unit cell of graphene. The doping is done such a way that the B/N−B/N distances are the longest to avoid severe structural distortion. Similarly, the B−N codoping is done by doping 1, 3, 5, 7, and 8 pairs of B−N atoms, which correspond to 6.25%, 18.75%, 31.25%, 43.75%, and 50.00% doping. The formation energy (Ef) of the B-, N-, and B−N doped graphene system is calculated using the following equations: Ef = Edoped − graphene − Egraphene − x(E B/N − EC)

Figure 3. Graphical representation of the change in doping formation energy of B-, N-, and B−N codoped graphene with increasing dopant concentration.

The TDOS and PDOS of 3.12% B@graphene is plotted to understand the electronic structure of B-doped graphene. We find some orbital density at the Fermi energy region and such orbital density coming from the dopants interacting with the neighboring carbon atoms to form degenerate states at the Fermi energy. Such doping shifts the Fermi energy (by 0.76 eV) below the Dirac point and thus makes B-doped graphene a p-type semiconductor as theoretically predicted by CervantesSodi and co-workers.26 Though such B-doping does not open up the band gap at the Fermi energy but creates a gap (0.20 eV) above the Fermi energy [Figure S2, Supporting Information], it behaves like a degenerate semiconductor.58,59 Thus, B-doped graphene’s electronic properties can be controlled by tuning the position of the Dirac point using the gate voltage.26,59 The boron doping concentrations are increased to open up the band gap at the Fermi energy. Therefore, we have doped 3, 5, 7, and 16 boron atoms in the 32-atom sheet of carbons, which corresponds to 9.38%, 15.60%, 21.80%, and 50.00% of boron doping. However, the increase of B doping concentration leads to a systematic opening up of band gap (0.20 to 2.56 eV) above the Fermi energy. After 50.00% of boron doping, 2.56 eV of band gap is produced above the Fermi energy. We assume as the doping concentration increases (hole concentration increases), the amount of degenerate states appearing at the Fermi energy level increases; thus, the effective band gap is shifted toward the more positive side [Figures S3− S6 and Table S3, Supporting Information]. The interaction of the gaseous molecules on the B-dopedgraphene is studied for low (3.12%) and high doping (50.00%) concentrations to get the minimum and maximum effect on the electronic structure of the doped graphene systems. The concentrations of the gas molecules are varied along with their doping concentrations. For the 3.12% doped-system, the dopant and gaseous molecular ratio is 1:1, whereas it is 16:4 for the 50.00% doping concentration. The dopant vs gaseous molecular ratio (16:4) is increased to minimize the interaction between the gaseous molecules. Gaseous molecules are most stable while adsorbed at the Bsite of the B-doped graphene (3.12%) [Figure 4a−d]. They are adsorbed via C-/O- and N-/O-center for CO/CO2 and NO/ NO2 molecules, respectively. Our relative energetic study shows CO and CO2 molecules are preferred to be adsorbing via C-

(3)

Ef = E B − N@graphene − Egraphene] − x(E B + E N) − 2EC (4)

Here, Edoped‑graphene and Egraphene are the energy of the optimized structure of doped-graphene and pure graphene sheet, respectively, and x is the number of dopant atoms. EB/N and EC is the atomic energy of doped (B/N) and carbon atom, respectively. EB and EC are calculated from the α-rhombohedral boron and graphene structure, respectively. EN is calculated from the total energy of the N2 molecule. We have varied the doping concentration from 3.12 to 50.00% for all the cases as B-, N-, and B−N doped graphene is experimentally achieved up to 13.00%, 15.60%, and 27.00%, respectively.55,56 So, we assume the B-, N-, and B−N doped graphene can be synthesized up to 50.00% doping concentration. The formation energy value for various doping concentrations is given in Supporting Information (Table S2). Our graphical representation of doping concentration vs formation energy [Figure 3] clearly shows that B−N codoping is energetically more favorable than N-doping and B-doping. Such trend also is observed experimentally,57 and this could be due to their (C and N) similar sizes (C = 0.73 Å and N = 0.71 Å) compared to the B atom (0.84 Å). 3.2.2. B@Graphene. As boron has one less valence electron than carbon, therefore such type of doping will lead into a ptype semiconductor. We have optimized the 3.12% B@ graphene, and the structure is slightly distorted as the B−C bond length is increased to 1.48 Å and thus deviated from the C−C bond length of 1.42 Å. The charge density plot show that there is an electron deficient area created around the B center, and the nature of the B−C bond is not purely covalent. Bader charge analysis also shows that B is positively charged (+1.97) while the neighboring C (−0.57) atoms are not. 24830

DOI: 10.1021/acs.jpcc.5b07359 J. Phys. Chem. C 2015, 119, 24827−24836

Article

The Journal of Physical Chemistry C

Figure 4. Optimized structures of (a) CO, (b) CO2, (c) NO, and (d) NO2 molecules adsorbed on 3.12% B@graphene surface (brown, green, blue, and red colored balls denote C, B, N, and O, respectively). Charge density plot of (e) CO, (f) CO2, (g) NO, and (h) NO2 molecules adsorbed on 3.12% B@graphene surface (isosurface value = 0.2322 e·Å−3). The CDD plots of (i) CO, (j) CO2, (k) NO, and (l) NO2 molecules adsorbed 3.12% B@graphene surface (isosurface value = 0.0036 e·Å−3).

Figure 5. Total and partial density of states (TDOS/PDOS) are plotted for (a) CO, (b) CO2, (c) NO, and (d) NO2 molecules adsorbed on 3.12% B@graphene surfaces (EF is set to zero and indicated by the black dashed line).

demonstrated a similar trend in gas adsorptions on B doped graphene surfaces.33,34 Zhang et al.33 reported that NO and NO2 bind more strongly than CO. Our calculated gaseous molecular adsorption energy trends are very much similar to previous theoretical reports on B@graphene surfaces.33,34 The calculated adsorption energy of the gases on the 50.00% B@graphene system shows similar trend as observed in 3.12% B@graphene. It ranges from 0.17 to 0.73 eV per gas molecules. We have also plotted the charge density and CDD for the 50.00% B-doped + gas systems to understand the extent of interactions. The adsorption energy and charge density analysis shows a similar trend in bonding as observed in 3.12% B@ graphene surface. The TDOS and PDOS are plotted [Figure 5a−d] for the 3.12% B@graphene + gas molecules. In the case of the 3.12% B@graphene + CO2 system, neither band gap opening nor any CO2 orbitals are present at the Fermi energy. It concludes that the presence of CO2 on B@graphene does not make any difference in the electronic structure of B@graphene. However, adsorption of CO and NO creates some impurity states at the Fermi energy level [Figure 5a,c]. PDOS analysis confirms the impurity states appearing at the EF mainly come from the gaseous C (CO) and N (NO) 2p orbitals. However, NO2 gas

center by 0.60 and 0.09 eV, respectively. The distance between the B and C of CO and CO2 gas molecules are 1.59 and 3.14 Å, respectively. The adsorption energies calculated for the most stable CO and CO2 structures are 0.66 and 0.12 eV, respectively. Therefore, CO interacts with the B-doped graphene strongly, whereas CO2 interacts weakly. This is because CO2 is a linear molecule; thus, it prefers to be adsorbed in a parallel fashion. As a result, the electrostatic interaction between the negatively charged O (of CO2) and the graphene π-electrons might be strong enough to repel the molecule. Similarly, NO and NO2 molecules are preferred to be adsorbing via N-center than the O-center [Figure 4c,d]. The relative energetics study show structures adsorbed via N-center are stable by 0.38 and 0.48 eV, respectively. The B−N distances are 2.22 and 1.65 Å for the NO and NO2 gases, respectively. The calculated adsorption energies for the most stable NO and NO2 structures are 0.58 and 1.23 eV, respectively. The charges on the N of NO2 are higher than on the NO; thus, the interaction between the B and the bent NO2 molecule is stronger than the NO. So, the adsorption energy and charge density analysis shows that CO/CO2 molecules are physisorbed, whereas NO/NO2 molecules are chemisorbed on the 3.12% B@graphene surface. Previous theoretical studies 24831

DOI: 10.1021/acs.jpcc.5b07359 J. Phys. Chem. C 2015, 119, 24827−24836

Article

The Journal of Physical Chemistry C

Figure 6. Optimized structures of (a) CO, (b) CO2, (c) NO, and (d) NO2 molecules adsorbed on 3.12% N@graphene surface (brown, blue, and red colored balls denote C, N, and O, respectively). Charge density plot of (e) CO, (f) CO2, (g) NO, and (h) NO2 molecules adsorbed on 3.12% N@ graphene surface (isosurface value 0.2322 e·Å−3). The CDD plots of (i) CO, (j) CO2, (k) NO, and (l) NO2 molecules adsorbed on 3.12% N@ graphene surface (isosurface value 0.0036 e·Å−3).

atomic unit cell of graphene, corresponding to 9.38%, 15.60%, 21.80%, to 50.00% of N doping. As the N doping concentration increases, the effective band gap of N@graphene system increases from 0.15 to 2.82 eV [Figures S7−S9, Supporting Information] below the Fermi energy. Here also some impurity states appear at the Fermi energy. This can also be explained as degenerate states coming from the N and C 2p orbitals. As N has more electrons than C, the degenerate states are coupled with the valence band (VB) side. These types of systems are called degenerate semiconductors.58 The gaseous molecular adsorptions are studied at the Ncenter of 3.12% N@graphene to see whether the interaction between the N and gaseous molecules could be improved or not. CO adsorbed through C-center is thermodynamically more stable by 0.01 eV than through the O-center [Figure 6a− d]. The B−C bond distance between the B and adsorbed CO gas molecule is 3.21 Å and thus weakly interacting. Similarly CO2 is 0.01 eV more stable when adsorbed via C-center than through the O-center, and the B−C distance is 3.08 Å. Thus, both the gases (CO and CO2) interact weakly with the Ndoped graphene. The adsorption energies calculated to be 0.11 and 0.12 eV for the most stable CO and CO2 structures, respectively. Similarly, NO and NO2 gases are preferred to be adsorbed via N-center by 0.02 and 0.16 eV, respectively. Their adsorption energies are calculated to be 0.18 and 0.26 eV, respectively. The distance between the N and adsorbed NO and NO2 gases are 3.05 and 3.02 Å, respectively. Therefore, the interaction between the N-doped graphene and the NO/NO2 gas molecules are comparably stronger. Our calculated adsorption energy trends are in agreement on previous studies on N-doped graphene surfaces reported by Zheng and coworkers.33 When the gaseous molecules adsorbed on the N@graphene surface, there is not much change in the Bader charges on the N and the surrounding C atoms, thus proving that there are only weak interactions between the gases and N@graphene. Only in the case of NO2 adsorption, the charge on the N atom changes (6.18) due to the stronger interaction with the NO2 molecule. Thus, NO2 has far higher adsorption energy in comparison to

adsorption opens up the band gap [0.50 eV, Figure 5d] at the Fermi energy. Hence, we can say B-doped graphene can be a selective and sensitive gas sensor for the NO2 molecule. Similarly, we have studied the electronic structure of the 50.00% B@graphene + gas systems. There is no band gap opening at the Fermi energy due to the gaseous (CO, CO2, NO, and NO2) molecular adsorption. For the higher doping concentration, the degenerate orbitals of B appear at the Fermi energy and thus not good for sensors. 3.2.3. N@Graphene. Similarly, we thought if we dope graphene with a more electro-negative atoms such as nitrogen (N), then the interaction between the nitrogen doped graphene and the gaseous molecules can be improved. However, nitrogen has one more valence electron than the carbon atom; therefore, such doping could lead to an n-type semiconductor. So, if we dope one N atom in the 32-atom unit cell of graphene, which corresponds to 3.12% of nitrogen doping. We have relaxed the N@graphene, and the optimized structure retains the similar geometry as pure graphene. This is because nitrogen and carbon are the same size; thus, the N−C bond distances (1.41 Å) are close to C−C bond distances (1.42 Å) of graphene. Charge density analysis shows more electron density around the N atom of N@graphene. Our Bader charge calculation also shows that N is more (−1.22) negatively charged than the surrounding C atoms (0.30). Thus, we assume the extra electron is localized mainly on the nitrogen atom. The TDOS and PDOS are plotted to understand the electronic structure of the N-doped graphene. From the TDOS analysis, we find band gap opens up by 0.20 eV below the Fermi energy. Here Fermi energy position shifted 0.76 eV above the Dirac point; thus, N@graphene is an n-type semiconductor. Some orbital density appearing just below the Fermi energy (−0.15 eV) and PDOS analysis confirms that this density is mainly coming from the N 2p. Thus, N-doped graphene’s electronic properties can be controlled by tuning the position of the Dirac point using gate voltage.26 We assumed if we increase the N doping concentration then we would open up the graphene’s band gap significantly. Hence, we have doped 3, 5, 7, and 16 nitrogen atoms in the 32 24832

DOI: 10.1021/acs.jpcc.5b07359 J. Phys. Chem. C 2015, 119, 24827−24836

Article

The Journal of Physical Chemistry C

Figure 7. Total and partial density of states (TDOS/PDOS) are plotted for (a) CO, (b) CO2, (c) NO, and (d) NO2 molecules adsorbed on 3.12% N@graphene surfaces (EF is set to zero and indicated by the black dashed line).

Figure 8. Optimized structures of (a) 6.25% B−N@graphene (front view) and (b) charge density plots of 6.25% B−N@graphene (isosurface value = 0.2322 e·Å−3). (c) TDOS and PDOS of 6.25% B−N@graphene (EF is set to zero and indicated by the black dashed line).

Figure 9. Optimized structures of (a) CO, (b) CO2, (c) NO, and (d) NO2 molecules adsorbed on 6.25% B−N@graphene surface. Charge density plots of (e) CO, (f) CO2, (g) NO, and (h) NO2 molecules absorbed on 6.25% B−N@graphene (isosurface value = 0.2322 e·Å−3). CDD of (i) CO, (j) CO2, (k) NO, and (l) NO2 (isosurface value = 0.0036 e·Å−3).

S4, Supporting Information] to understand the effect on the electronic structure of the N@graphene due to the toxic gases adsorbed. The electronic structure of the 3.12% and 50.00% N@graphene + gas system does not open up any band gap at the Fermi energy. It means adsorption of the toxic gases did not alter the electrical current flow through the 3.12% N@ graphene. Thus, N@graphene (3.12%) cannot be a promising semiconductor based gas sensor. 3.2.4. B−N@Graphene. We have considered B−N codoping to check whether these types of doping could open up the band gap at the Fermi energy and such doping might improve the selectivity and sensitivity of the gas sensor. Therefore, we have doped one boron and nitrogen atom into the six-membered

the other three gases. The charge density and CDD [Figure 6e−l] is plotted to understand whether there is any orbital mixing present between the N- and gaseous molecules. The CDD analysis shows there is not enough orbital mixing between the N-doped graphene and the gas molecules. The calculated adsorption energy of the gases on the 50.00% N@graphene system ranges from 0.16 to 0.52 eV per gas molecule. The charge density and CDD plot confirms that there is only electrostatic attraction between the N atom and the adsorbed gases. The electronic structure (TDOS and PDOS) [Figure 7a−d] of the N@graphene + gas molecules is studied for 3.12 and 50.00% N@graphene systems [Figures S10 and S11 and Table 24833

DOI: 10.1021/acs.jpcc.5b07359 J. Phys. Chem. C 2015, 119, 24827−24836

Article

The Journal of Physical Chemistry C

Figure 10. Total and partial density of states (TDOS/PDOS) are plotted for (a) CO, (b) CO2, (c) NO, and (d) NO2 molecules adsorbed on 6.25% B−N@graphene surface (EF is set to zero and indicated by the black dashed line].

Similarly, we have studied the adsorption property of these gas molecules on 50.00% B−N@graphene as well. The calculated adsorption energy of the gases on the 50.00% B− N@graphene system ranges from 0.24 to 1.23 eV per gas molecule. The charge density and CDD is plotted for the 50.00% B−N@graphene system, and it confirms that the gas molecules are chemisorbed on the 50.00% B−N@graphene surface. TDOS and PDOS of 6.25% B−N@graphene + gas molecules are plotted in Figure 10a−d. The electronic structure of 6.25% B−N@graphene+CO/CO2 molecules [Figure 10a,b] shows that there is not much change in the electronic structure when CO and CO2 molecules are adsorbed. As the band gap changes from 0.24 to 0.43 for CO and 0.24 to 0.35 for CO2. In the case of CO, the N 2p (loan pair) orbital appears just below the Fermi energy and the vacant O 2p orbital just above the Fermi energy. However, there are impurity states appearing at the Fermi energy when NO/NO2 gases are adsorbed on the B− N@graphene surface [Figure 10c,d]. We have found the same type of interaction of NO with the B−N@graphene surface as we have discussed in B@graphene. So, there is no band gap at the Fermi energy, and we find 6.25% B−N@graphene is highly sensitive and selective for N-based (NO and NO2) gas molecules. Similarly, the 50.00% B−N@graphene shows a similar trend. So, the B−N@graphene could be a highly sensitive and selective semiconducting based gas sensor.

ring of the 32 atomic unit cell of graphene [Figure 8a] in such a way that the interaction between the B and N is minimal. The charge density study [Figure 8b] shows that there is less electron density around the B atom but more on the N atom as found in their respective B@graphene and N@graphene systems. Bader charge analysis shows a similar trend as B and N bears +1.97 and −1.19 charges, respectively [Table S5, Supporting Information]. Interestingly, in the monodoped case (B@graphene and N@graphene), the 2p orbitals of the dopants are mainly present at the Fermi energy; thus, we do not get clean band gap at the Fermi energy. Even, if we increase the concentration of B and N doping, the effective band gap shifts further and further from the Fermi energy. However, if we codope (6.25% B−N doping), then the band gap opens up (0.24 eV) at the Fermi energy [Figure 8c]. This is because the total number of electrons in the system remains constant and thus gives a clean band gap at the Fermi energy. We have also increased the B−N concentration in the B−N@graphene from 6.25% to 18.75%, 31.25%, 43.75%, and 50.00%, and the calculated band gaps are 0.24, 0.90, 0.30, 0.19, and 0.11 eV, respectively [Figures S12 and S13 and Table S6, Supporting Information]. As the B−N codoping pattern is different for the higher concentrations (31.25% to 50.00%), the band gap openings are different. First, we have considered 6.25% B−N@graphene for our further study. The gaseous molecules (CO, CO2, NO, and NO2) are adsorbed at the B-site of the 6.25% B−N@graphene as B interacts more strongly with the gas molecules in comparison to N. The adsorption energies calculated for all these gaseous molecules were found to be stronger (0.25 to 1.35 eV) with respect to their monodoped systems. Among these gases, NO2 (1.36 eV), NO (0.85 eV), and CO (0.77 eV) bind very strongly, whereas CO2 (0.25 eV) weakly binds on the B−N@ graphene surface [Figure 9a−d]. We have plotted the charge density and CDD of the gases adsorbed on the B−N@ graphene surface [Figure 9e−l] to see the interaction between the B−N@graphene and gas molecules. The charge density and CDD also shows that NO2, NO, and CO have considerable orbital mixing with the B of the B−N@graphene. Bader charge analysis shows that B gains 0.12e charge (+1.85) when CO is adsorbed and loses 0.11e (+2.08) when NO2 is adsorbed. As we know CO is a reducing gas so it transfers electrons to the B, whereas NO2 is an oxidizing gas and thus oxidizes B. Such considerable amount of electron transfer is only possible when they are bonded, thus proving that CO and NO2 gases are chemisorbed on B−N@graphene surface.

4. CONCLUSION Here, we have performed density functional theory (DFT) calculations to study the gas (CO, CO2, NO, and NO2) sensing mechanism of pure and doped (B@, N@ ,and B−N@) graphene surfaces. The adsorption energy, charge density, and CDD of the pure graphene + gas molecules are studied to understand their bonding and selectivity. The p- and n-type of doping certainly improves the interaction between the surface and the gaseous molecules but such doping does not change the band gap properties at the Fermi energy. Therefore, semiconductor-based gas sensing technique is not promising for such type of doping. However, the B−N codoping changes the graphene’s electronic properties at the Fermi energy. Once the gaseous molecules are adsorbed on the B−N@graphene surface, the electronic properties change maximum for NO and NO2 but not for CO and CO2. Therefore, B−N@graphene could be a promising semiconductor based gas sensor. 24834

DOI: 10.1021/acs.jpcc.5b07359 J. Phys. Chem. C 2015, 119, 24827−24836

Article

The Journal of Physical Chemistry C



(15) Rumyantsev, S.; Liu, G.; Shur, M. S.; Potyrailo, R. A.; Balandin, A. A. Selective Gas Sensing with a Single Pristine Graphene Transistor. Nano Lett. 2012, 12, 2294−2298. (16) Schedin, F.; Geim, A. K.; Morozov, S. V.; Hill, E. W.; Blake, P.; Katsnelson, M. I.; Novoselov, K. S. Detection of Individual Gas Molecules Adsorbed on Grapheme. Nat. Mater. 2007, 6, 652−655. (17) Yuan, W.; Shi, G. Graphene-Based Gas Sensors. J. Mater. Chem. A 2013, 1, 10078−10091. (18) Ratinac, K. R.; Yang, W.; Ringer, S. P.; Braet, F. Toward Ubiquitous Environmental Gas Sensors-Capitalizing on the Promise of Graphene. Environ. Sci. Technol. 2010, 44, 1167−1176. (19) Rajan, A. C.; Rezapour, R. M.; Yun, J.; Cho, Y.; Cho, W. J.; Min, S. K.; Lee, G.; Kim, K. S. Two Dimensional Molecular Electronics Spectroscopy for Molecular Fingerprinting, DNA Sequencing and Cancerous DNA Recognition. ACS Nano 2014, 8, 1827−1833. (20) Min, S. K.; Kim, W. Y.; Cho, Y.; Kim, K. S. Fast DNA Sequencing with a Graphene-Based Nanochannel Device. Nat. Nanotechnol. 2011, 6, 162−165. (21) Prasongkit, J.; Grigoriev, A.; Pathak, B.; Ahuja, R.; Scheicher, R. H. Transverse Conductance of DNA Nucleotides in a Graphene Nanogap from First Principles. Nano Lett. 2011, 11, 1941−1945. (22) Prasongkit, J.; Grigoriev, A.; Pathak, B.; Ahuja, R.; Scheicher, R. H. Theoretical Study of Electronic Transport through DNA Nucleotides in a Double-Functionalized Graphene Nanogap. J. Phys. Chem. C 2013, 117, 15421−15428. (23) Cai, J.; Pignedoli, C. A.; Talirz, L.; Ruffieux, P.; Söde, H.; Liang, L.; Meunier, V.; Berger, R.; Li, R.; Feng, X.; et al. Graphene Nanoribbon Heterojunctions. Nat. Nanotechnol. 2014, 9, 896−900. (24) Wang, X.; Sun, G.; Routh, P.; Kim, D.-H.; Huang, W.; Chen, P. Heteroatom-Doped Graphene Materials: Syntheses, Properties and Applications. Chem. Soc. Rev. 2014, 43, 7067−7098. (25) Lherbier, A.; Blase, X.; Niquet, Y.-M.; Triozon, F.; Roche, S. Charge Transport in Chemically Doped 2D Graphene. Phys. Rev. Lett. 2008, 101, 036808. (26) Cervantes-Sodi, F.; Csányi, G.; Piscanec, S.; Ferrari, A. C. Edgefunctionalized and Substitutionally Doped Graphene Nanoribbons: Electronic and Spin Properties. Phys. Rev. B: Condens. Matter Mater. Phys. 2008, 77, 165427. (27) Panchakarla, L. S.; Subrahmanyam, K. S.; Saha, S. K.; Govindaraj, A.; Krishnamurthy, H. R.; Waghmare, U. V.; Rao, C. N. R. Synthesis, Structure and Properties of Boron and Nitrogen Doped Graphene. Adv. Mater. 2009, 21, 4726−4730. (28) Sheng, Z.-H.; Gao, H.-L.; Bao, W.-J.; Wang, F.-B.; Xia, X.-H. Synthesis of Boron Doped Graphene for Oxygen Reduction Reaction in Fuel Cells. J. Mater. Chem. 2012, 22, 390−395. (29) Wei, D.; Liu, Y.; Wang, Y.; Zhang, H.; Huang, L.; Yu, G. Synthesis of N-Doped Graphene by Chemical Vapor Deposition and Its Electrical Properties. Nano Lett. 2009, 9, 1752−1758. (30) Fan, X.; Shen, Z.; Liu, A. Q.; Kuo, J.-L. Band Gap Opening of Graphene by Doping Small Boron Nitride Domains. Nanoscale 2012, 4, 2157−2165. (31) Bepete, G.; Voiry, D.; Chhowalla, M.; Chiguvare, Z.; Coville, N. J. Incorporation of Small BN Domains in Graphene During CVD Using Methane, Boric Acid and Nitrogen Gas. Nanoscale 2013, 5, 6552−6557. (32) Ko, G.; Kim, H.-Y.; Ahn, J.; Park, Y.-M.; Lee, K.-Y.; Kim, J. Graphene-Based Nitrogen Dioxide Gas Sensors. Curr. Appl. Phys. 2010, 10, 1002−1004. (33) Zhang, Y.-H.; Chen, Y.-B.; Zhou, K.-G.; Liu, C.-H.; Zeng, J.; Zhang, H.-L.; Peng, Y. Improving Gas Sensing Properties of Graphene by Introducing Dopants and Defects: A First-Principles Study. Nanotechnology 2009, 20, 185504. (34) Dai, J.; Yuan, J.; Giannozzi, P. Gas Adsorption on Graphene Doped with B, N, Al, and S: A Theoretical Study. Appl. Phys. Lett. 2009, 95, 232105. (35) Lv, Y.-a; Zhuang, G.-l; Wang, J.-g; Jia, Y.-b.; Xie, Q. Enhanced Role of Al or Ga-Doped Graphene on the Adsorption and Dissociation of N2O Under Electric Field. Phys. Chem. Chem. Phys. 2011, 13, 12472−12477.

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpcc.5b07359. Optimized structures (graphene, graphene + gases, B@ graphene, 50% B@graphene + gases, N@graphene, 50.00% N@graphene + gases, and 50.00% B−N@ graphene) and their respective TDOS/PDOS, charge densities, and Bader charges (PDF)



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS We thank IIT Indore for the lab and computing facilities. This work is supported by Council of Scientific and Industrial Research [CSIR, Grant number: 01(2723)/13/EMR(II)], New Delhi. A.M. thanks CSIR for the Senior Research Fellowship. R.A. thanks Swedish Research Council (VR), Swedish Institute, and Swedish Energy Agency for financial support and SNIC for computing time.



REFERENCES

(1) Wang, X.-d.; Wolfbeis, O. S. Optical Methods for Sensing and Imaging Oxygen: Materials, Spectroscopies and Applications. Chem. Soc. Rev. 2014, 43, 3666−3761. (2) Bai, J.; Zhou, B. Titanium Dioxide Nanomaterials for Sensor Applications. Chem. Rev. 2014, 114, 10131−10176. (3) Wang, C.; Yin, L.; Zhang, L.; Xiang, D.; Gao, R. Metal Oxide Gas Sensors: Sensitivity and Influencing Factors. Sensors 2010, 10, 2088− 2106. (4) Nisar, J.; Topalian, Z.; Sarkar, A. D.; Ö sterlund, L.; Ahuja, R. TiO2-Based Gas Sensor: A Possible Application to SO2. ACS Appl. Mater. Interfaces 2013, 5, 8516−8522. (5) Liu, S.; Hu, F.; Zhang, J.; Tang, H.; Shao, M. Surface-Doping Effect of InVO4 Nanoribbons and the Distinctive Behaviour as Gas Sensors. ACS Appl. Mater. Interfaces 2013, 5, 3208−3211. (6) Karlsson, H. L.; Cronholm, P.; Gustafsson, J.; Möller, L. Copper Oxide Nanoparticles are Highly Toxic: A Comparison between Metal Oxide Nanoparticles and Carbon Nanotubes. Chem. Res. Toxicol. 2008, 21, 1726−1732. (7) Huang, X.; Qi, X.; Boey, F.; Zhang, H. Graphene-Based Composites. Chem. Soc. Rev. 2012, 41, 666−686. (8) Huang, X.; Yin, Z.; Wu, S.; Qi, X.; He, Q.; Zhang, Q.; Yan, Q.; Boey, F.; Zhang, H. Graphene-Based Materials: Synthesis, Characterization, Properties, and Applications. Small 2011, 7, 1876−1902. (9) Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Jiang, D.; Zhang, Y.; Dubonos, S. V.; Grigorieva, I. V.; Firsov, A. A. Electric Field Effect in Atomically Thin Carbon Films. Science 2004, 306, 666−669. (10) Wallace, P. R. The Band Theory of Graphite. Phys. Rev. 1947, 71, 622−634. (11) McClure, J. W. Diamagnetism of Graphite. Phys. Rev. 1956, 104, 666−671. (12) Novoselov, K. S.; Jiang, D.; Schedin, F.; Booth, T. J.; Khotkevich, V. V.; Morozov, S. V.; Geim, A. K. Two-Dimensional Atomic Crystals. Proc. Natl. Acad. Sci. U. S. A. 2005, 102, 10451− 10453. (13) Huang, X.; Zeng, Z.; Fan, Z.; Liu, J.; Zhang, H. Graphene-Based Electrodes. Adv. Mater. 2012, 24, 5979−6004. (14) He, Q.; Wu, S.; Yin, Z.; Zhang, H. Graphene-Based Electronic Sensors. Chem. Sci. 2012, 3, 1764−1772. 24835

DOI: 10.1021/acs.jpcc.5b07359 J. Phys. Chem. C 2015, 119, 24827−24836

Article

The Journal of Physical Chemistry C

(57) Rao, C. N. R.; Gopalakrishnana, K.; Govindaraj, A. Synthesis, Properties and Applications of Graphene Doped with Boron, Nitrogen and Other Elements. Nano Today 2014, 9, 324−343. (58) Lee, K.-W.; Pickett, W. E. Superconductivity in Boron-Doped Diamond. Phys. Rev. Lett. 2004, 93, 237003. (59) Novoselov, K. S. Graphene: Mind the Gap. Nat. Mater. 2007, 6, 720−721.

(36) Niu, F.; Liu, J.-M.; Tao, L.-M.; Wang, W.; Song, W.-G. Nitrogen and Silica Co-Doped Graphene Nanosheets for NO2 Gas Sensing. J. Mater. Chem. A 2013, 1, 6130−6133. (37) Guo, P.; Xiao, F.; Liu, Q.; Liu, H.; Guo, Y.; Gong, J. R.; Wang, S.; Liu, Y. One-Pot Microbial Method to Synthesize Dual-Doped Graphene and Its Use as High-Performance Electrocatalyst. Sci. Rep. 2013, 3, 3499. (38) Yang, G.-H.; Zhou, Y.-H.; Wu, J.-J.; Cao, J.-T.; Li, L.-L.; Liu, H.Y.; Zhu, J.-J. Microwave-Assisted Synthesis of Nitrogen and Boron CoDoped Graphene and Its Application for Enhanced Electrochemical Detection of Hydrogen Peroxide. RSC Adv. 2013, 3, 22597−22604. (39) Qu, D.; Zheng, M.; Du, P.; Zhou, Y.; Zhang, L.; Li, D.; Tan, H.; Zhao, Z.; Xie, Z.; Sun, Z. Highly Luminescent S, N Co-Doped Graphene Quantum Dots with Broad Visible Absorption Bands for Visible Light Photocatalysts. Nanoscale 2013, 5, 12272−12277. (40) Dong, Q.; Zhuang, X.; Li, Z.; Li, B.; Fang, B.; Yang, C.; Xie, H.; Zhang, F.; Feng, X. Efficient Approach to Iron/Nitrogen Co-Doped Graphene Materials as Efficient Electrochemical Catalysts for the Oxygen Reduction Reaction. J. Mater. Chem. A 2015, 3, 7767−7772. (41) Duan, X.; O’Donnell, K.; Sun, H.; Wang, Y.; Wang, S. Sulfur and Nitrogen Co-Doped Graphene for Metal-Free Catalytic Oxidation Reactions. Small 2015, 11, 3036−3044. (42) Wu, Z. S.; Winter, A.; Chen, L.; Sun, Y.; Turchanin, A.; Feng, X.; Müllen, K. Three-Dimensional Nitrogen and Boron Co-Doped Graphene for High-Performance All-Solid-State Supercapacitors. Adv. Mater. 2012, 24, 5130−5135. (43) Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for ab initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B: Condens. Matter Mater. Phys. 1996, 54, 11169−11186. (44) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (45) Kresse, G.; Joubert, D. From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method. Phys. Rev. B: Condens. Matter Mater. Phys. 1999, 59, 1758−1775. (46) Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate ab initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104. (47) Monkhorst, H. J.; Pack, J. D. Special Points for Brillouin-Zone Integrations. Phys. Rev. B 1976, 13, 5188−5192. (48) Mahata, A.; Rai, R. K.; Choudhuri, I.; Singh, S. K.; Pathak, B. Direct vs. Indirect Pathway for Nitrobenzene Reduction Reaction on a Ni Catalyst Surface: A Density Functional Study. Phys. Chem. Chem. Phys. 2014, 16, 26365−26374. (49) Heyrovska, R. Atomic Structures of Graphene, Benzene and Methane with Bond Lengths as Sums of the Single, Double and Resonance Bond Radii of Carbon, arXiv:0804.4086, 2008. (50) Bader, R. F. W. A Quantum Theory of Molecular Structure and Its Applications. Chem. Rev. 1991, 91, 893−928. (51) Henkelman, G.; Arnaldsson, A.; Jónsson, H. A Fast and Robust Algorithm for Bader Decomposition of Charge Density. Comput. Mater. Sci. 2006, 36, 354−360. (52) Posternak, M.; Baldereschi, A.; Freeman, A. J.; Wimmer, E.; Weinert, M. Prediction of Electronic Interlayer States in Graphite and Reinterpretation of Alkali Bands in Graphite Intercalation Compounds. Phys. Rev. Lett. 1983, 50, 761−764. (53) Drissi, L. B.; Saidi, E. H.; Bousmina, M.; Fassi-Fehri, O. DFT Investigations of the Hydrogenation Effect on Silicene/Graphene Hybrids. J. Phys.: Condens. Matter 2012, 24, 485502. (54) Yavari, F.; Koratkar, N. Graphene-Based Chemical Sensors. J. Phys. Chem. Lett. 2012, 3, 1746−1753. (55) Xu, J.; Jang, S. K.; Lee, J.; Song, Y. J.; Lee, S. The Preparation of BN-Doped Atomic Layer Graphene via Plasma Treatment and Thermal Annealing. J. Phys. Chem. C 2014, 118, 22268−22273. (56) Chang, C.-K.; Kataria, S.; Kuo, C.-C.; Ganguly, A.; Wang, B.-Y.; Hwang, J.-Y.; Huang, K.-J.; Yang, W.-H.; Wang, S.-B.; Chuang, C.-H.; et al. Band Gap Engineering of Chemical Vapor Deposited Graphene by in Situ BN Doping. ACS Nano 2013, 7, 1333−1341. 24836

DOI: 10.1021/acs.jpcc.5b07359 J. Phys. Chem. C 2015, 119, 24827−24836