The Effect of Boron and Nitrogen Doping in ... - ACS Publications

Oct 27, 2016 - Department of Physics, Assam University, Silchar-788011, India. •S Supporting Information. ABSTRACT: The electronic, magnetic, and op...
0 downloads 0 Views 6MB Size
Article pubs.acs.org/JPCC

The Effect of Boron and Nitrogen Doping in Electronic, Magnetic, and Optical Properties of Graphyne Barnali Bhattacharya and Utpal Sarkar* Department of Physics, Assam University, Silchar-788011, India

Downloaded via IOWA STATE UNIV on January 27, 2019 at 06:20:36 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.

S Supporting Information *

ABSTRACT: The electronic, magnetic, and optical properties of boronand nitrogen-doped graphyne have been investigated with various doping positions and concentrations of boron and nitrogen atoms. We have explored how the presence of a single dopant atom changes the conductivity of doped graphyne from the semiconducting to metallic one. The boron atom at the chain site introduces spin polarization which is in the ferromagnetic (FM) ground state for minimal boron concentration and in the antiferromagnetic (AFM) ground state for an increasing number of boron atoms in the unit cell. We have examined the origin of spin polarization which increases with increasing dopant concentration. Our optical spectra show that the interband transition takes place in the low energy regime. Due to the presence of dopant atom, the absorption spectra extend from the infrared region to the UV region and exhibit a strong peak. The reflectivity and energy loss spectra derived the plasmon energy for these systems where the reflectivity displays a sharp decline.



INTRODUCTION Low dimensional structures of carbon allotropes have become the principal focus of the modern research community, due to their novel intriguing properties like field emission,1,2 quantum conductance,3,4 superconductivity,5−8 half-metallicity,9,10 high thermal conductivity,11−14 high capacitance,15−21 and transport.22,23 In this regard, graphene is the archetype of the twodimensional (2D) network of carbon having an sp2 hybridized carbon atom and it has been projected as a promising candidate for future nanoscale electronic devices.24,25 The pristine graphene is a semimetal and has some limitations for its use in the logic gate and other electronic circuits. Consequently, doped graphene offers an alternative pathway over their parent counterpart and also offers a significant potential as an electrocatalyst.26−30 Baughman et al.31 predicted a new type of carbon allotrope called graphyne having the same symmetry as graphene but being structurally different from graphene. Graphyne is a direct-band gap semiconductor,32−37 which enhances its applications in future nanoelectronics. It is encouraging that, the graphdiyne tube and film,38−42 flakes and substructures of graphyne43−45 are synthesized successfully using various chemical methods. However, it is in the initial stage toward the preparation of extended structures of graphyne. Additionally, the first principle calculations also established the stability of graphyne and its related derivatives.36,37,46 Chemical doping is a well-adopted and effective method for tailoring the electronic and optical properties of nanostructure.47−67 Doping helps to tune the work function and carrier concentration of nanostructure which broadens the range of © 2016 American Chemical Society

electronic and optical applications. Now, doping of graphene by adding an external atom or substitution has effectively tuned the electronic and optical properties of graphene, resulting in an application. Recent studies on N-doped graphene show excellent electrocatalytic activity,48−50 whereas Al and Si doping enhances the sensing activity51−53 of graphene by increasing the molecule graphene interaction. Nitrogen or silicon doping led to the appearance of a band gap in the graphene electronic spectrum.54,55 Current studies also indicate that nitrogen- and boron-doped graphene,56−61 nitrogen- and boron-doped carbon nanotubes,62,63 and carbon nitride64 could offer numerous optoelectronic properties61 including nonlinear optical properties59 as well as catalytic behavior over their traditional counterpart. Recently, many experiments have been performed to modify the properties of graphene with foreign atoms including nitrogen, boron, fluorine, and hydrogen.48,65,66 The successful preparation and characterization of B- and Ndoped graphene by Panchakarla et al.48 confirms that doped graphene is highly possible to synthesize and can be used to fabricate electronic devices, as the B- and N-doped graphene exhibit p- and n-type semiconducting properties, respectively, that can be tuned by varying the concentration of dopant atom. The advances on doped graphene gradually grow interest among the research community for the modification of mechanical and optical properties of graphyne with dopant atom. The previous study of BN-doped graphyne,67 graphyne Received: July 26, 2016 Revised: October 16, 2016 Published: October 27, 2016 26793

DOI: 10.1021/acs.jpcc.6b07478 J. Phys. Chem. C 2016, 120, 26793−26806

Article

The Journal of Physical Chemistry C nanoribbons,68 and graphyne nanotube69 exhibit an increase in band gap and shift the optical absorption spectra toward the UV region depending on the doping site and concentration. The possibility of using Ca-decorated graphyne70,71 and graphyne nanotube72 as promising hydrogen storage materials opens a new prospect for the study of doped graphyne. The applicability of Li-intercalated multilayer graphyne, as a lithium ion battery anode, is reported by Hwang et al.73 The Na decorated graphyne sheet37 has also been studied. Unlike B- or N-doped graphene, a little information is available in the literature about the effect of boron and nitrogen doping with different concentrations on the electronic properties of graphyne. Although the electronic and optical properties of BN-codoped graphyne have been explored, only boron- and nitrogen-doped graphyne74−76 with different concentrations has been less explored. Very recently, the B- and N-doped graphyne was studied by a few researchers to explore their electrocatalytic activity,75,77 oxygen reduction activity,78 and possible application as lithium-ion79 and hydrogen storage material.80 However, basic questions about the dopant distribution and how their electronic, optical, and magnetic properties change remain an open issue. To address these fundamental questions and to explore the possibility of using graphyne more widely in the field of electronics, optoelectronics, and spinotronics, modification in electronic structure and in optical property is necessary as a consequence of which a detailed investigation by doping is essential. Our study is complete in this sense and can serve the purpose suitably. In addition, to the best of our knowledge, an optical properties study of only boron- or nitrogen-doped graphyne has yet to be reported. Being motivated by the above, in the present study, we have explored systematically the electronic, magnetic, and optical properties of boron- and nitrogen-doped graphyne. Tailoring the above-mentioned properties with increasing dopant concentration is also discussed.

the number of carbon and dopant atoms present in the doped configuration, respectively. The optical properties of a system can be deduced from the frequency dependent complex dielectric function ε(ω) = ε1(ω) + iε2(ω).61,83,84 The imaginary part of the dielectric function, due to direct interband transition, can be described as ε2(ω)inter =

∑ ∑ Wk |ρija |2 δ(εkj − εki − ω) i ∈ VB, j ∈ CB

k

(1)

ρaij,

Here Ω, εkj, and εki represent the unit cell volume, dipole transition matrix, energy of occupied valence states, and energy of unoccupied conduction states, respectively. When the metallic intraband transition is added, the imaginary part of the dielectric function85,86 can be obtained using the relation ε2(ω)intra + inter =

Γωpl 2 ω(ω 2 + Γ 2)

+ ε2(ω)inter

(2)

in which ω, ωpl, and Γ stand for the photon energy, plasma frequency, and relaxation time in the Drude model, respectively. The relaxation time is taken as the same as alpha graphyne,87 as our considered doped graphyne structures are also metallic like alpha graphyne. The real part of the dielectric function depends on the imaginary part, and they are related to each other by the Kramer−Kronig relation.88,89 Using this, the real part of the dielectric function for interband transition can be obtained as ε1(ω)inter = 1 +

4 P π

∫0



dω′

ω′ε2(ω′) ω′2 − ω 2

(3)

Here P is the principal value of the integral. For the metallic system, the real part of the dielectric function involves the contribution of intraband and interband transitions and ε1(ω)intra can be expressed as



intra

ε1(ω)

COMPUTATIONAL DETAILS Spin polarized density functional theory (DFT) calculations were performed using the siesta package81 within the generalized gradient approximation (GGA) and Perdew− Burke−Ernzerhof (PBE) is employed to treat the exchangecorrelation part of the density functional. The DZP basis set has been applied, and the norm-conserving Troullier−Martins pseudopotential82 is utilized to describe the core−valence interactions. Structural relaxations are performed using the conjugate-gradient method until the maximum forces on each atom are less or equal to 0.010 eV/Å and the energy difference between the consecutive self-consistent field steps is less than 10−4 eV. The 11 × 11 × 1 Monkhorst−Pack set of k points is used for sampling the Brillouin zone, and the kinetic energy cutoff value was set at 300 Ry. To avoid the interaction between two neighboring images, periodic boundary condition along the z-axis is employed with a lateral vacuum region equal to 15 Å. The stability is determined through cohesive energy (Ecoh), which is calculated as Ecoh =

4π 2 Ωω 2

=1−

ωpl 2 ω(ω 2 + Γ 2)

(4)

From the real and imaginary parts of the dielectric function, the optical properties such as absorption coefficient (α),61,67,84 reflectivity (R),61,67,84 real part of optical conductivity (σ),89,90 and energy loss spectrum (L)61,67,83,84 can be obtained.



RESULTS AND DISCUSSION In the present study, we explore the influence of dopant atom, doping site, and concentration on the electronic and optical properties of graphyne. Our previous reports36,67,69 established the tunability of the electronic and optical properties of the graphyne family due to codoping with BN and motivated the present study. For substitution purposes, we have chosen two high symmetric sites of carbon atom having sp or sp2 hybridization (chain or ring). In addition, we have gradually increased the doping concentration step by step. The optimized 2D layered structures (Figure 1) of B- or N-doped graphyne corresponding to the favorable doping site ensure that substitutional doping of B or N brings a small distortion in local area only. A. Detailed Structures. Depending on dopant type and doping site, the bond length and lattice constant of B- or Ndoped graphyne get changed (Table S1). Due to doping, the deformation potential energy is developed by the lattice deformation and the hybrid orbitals are formed between

E D ‐ graphyne − [nEC + mEdopant] n+m

where ED‑graphyne, EC, and Edopant represent the total energies of the doped graphyne, isolated carbon, and isolated dopant (boron or nitrogen) atom, respectively. Here, n and m denote 26794

DOI: 10.1021/acs.jpcc.6b07478 J. Phys. Chem. C 2016, 120, 26793−26806

Article

The Journal of Physical Chemistry C

results in the same type of bond elongation (contraction), but the triple bond character in the BC (NC) bond is reduced. For “graphyne with 3B at chain”, the BC (NC) bond connecting the ring and chain gets contracted compare to “graphyne with 1B (N) at chain” and “graphyne with 2B (N) at chain”. The elongation of bond length is also noticed when the sp2 hybridized carbon atom is substituted by a boron atom, as found in “B-doped graphene” reported by Panchakarla et al.48 where they observed that the B−C bond is about 0.5% longer than the C−C bond. The B−C bond in the hexagon is the same as reported for the B−C single bond92 but slightly larger than that observed in “B-doped graphene”,95 and the B−C bond connecting the ring and chain is also nearly equal to the B−C single bond. The result gets changed when the sp2 hybridized C atom is substituted by a N atom. Same as “N-doped graphene”,48,95 the N−C bond length (1.406 Å) in hexagon is nearly the same as the C−C bond length which is much greater than the N−C double bond length (1.270 Å)93 but smaller than the reported N−C single bond length (1.474 Å).94 Hence, this bond has partial double bond character. Unlike the B−C bond connecting the ring and chain, the N−C bond connecting the N atom and the sp hybridized C atom at the chain is found to have strong π character, as the bond length (1.330 Å) is close to the reported N−C double bond (1.270 Å). Upon analyzing the optimized structure and bond length from Table S2, we can infer that boron doping expands the bond length, whereas the nitrogen doping reduces the bond length. Due to substitutional doping, lattice deformation arises and simultaneously a residual stress is generated. To release the generated residual stress, the B−C (N−C) bonds get stretched (squeezed). The calculated cohesive energies (Table S1) indicate the stability of our substituted graphyne systems and are comparable with the experimental cohesive energy of graphite (7.37 eV).96,97 The cohesive energy trend makes it clear that with the same doping concentration and doping site, substitutional B doping is comparatively more favorable than N doping. Similar results were reported by Jafari et al.68 for Bdoped graphyne nanoribbons (GNRs) where the lowest formation energy is assigned to B-doped GNR inferring their stability. The lowest formation energy concludes a greater stability of “B-doped graphene” than “N-doped graphene” and is reported by Ren et al.98 The most favorable site for B doping is the sp2 bonded C atom site instead of the sp hybridized C atom site as the former offers higher coordination number for bonding and thus facilitates better stability. However, for N doping, the most preferred site for substitution is the sp hybridized C atom at the chain. Generally, the breaking of two partial double bonds requires more energy than the breaking of one triple bond; therefore, the substitution of the sp-hybridized C atom by a dopant atom is energetically favorable which is observed for nitrogen doping too. However, boron doping breaks three partial double bonds to generate three single bonds. Actually, when the sp2-hybridized C atom is substituted by a B atom, a slight movement is sufficient to stretch the C−B bond and it needs a small energy to break the CC bond at the ring. Interestingly, the cohesive energy of B (N)-doped graphyne decreases with increasing doping concentration. B. Electronic Properties. The previous study of pristine graphyne and BN-codoped graphyne36,67 established their nonmagnetic nature. It is clear from Table S1 that only boron substituted at the chain position introduces spin

Figure 1. Optimized structure of doped graphyne with the most favorable doping site: (a) graphyne with 1B at ring; (b) graphyne with 2B at ring; (c) graphyne with 3B at ring; (d) graphyne with 1N at chain; (e) graphyne with 2N at chain; (f) graphyne with 3N at chain. The boron or nitrogen atoms are substituted in the unit cell.

dopant and carbon atom. Hence, the energy involved in the substitution process contains two parts; one part breaks the original C−C bond, and the other part is the deformation potential energy. We see that boron doping elongated the lattice vectors (both a and b) compared to pristine graphyne but nitrogen doping shortened the lattice vectors except for “graphyne with 1N at chain” and “graphyne with 2N at chain”. The elongation/shortening of the lattice vector is due to the larger/shorter atomic radius of B (85 pm)/N (65 pm) than that of the C (70 pm) atom. When sp hybridized C atom is substituted by B atom, two types of CC bonds have been broken; one is the strong triple bond, and the other has a double bond characteristic which is comparatively stronger than the CC bond at the hexagon of graphyne. As a result, the “graphyne with 1B at chain” contains two types of BC bonds with bond lengths of 1.361 and 1.498 Å, respectively. The former bond length is slightly less than the reported BC double bond (1.440 Å)91 but greater than the BC triple bond (1.230 Å) found in B-doped α-graphyne,92 and hence, this bond is expected to have mainly double bond character along with BC triple bond character. In contrast, the BC bond between the chain and ring is shorter than the BC single bond (1.522 Å)91 but larger than the BC double bond and thus contains some π character which matches well with other findings.91 Now the BC bond length at the chain site is slightly elongated compared to the pristine graphyne. To compensate this elongation, rearrangement of all bond lengths of the graphyne unit cell takes place. The nearest CC bond length connecting the ring and chain gets squeezed, but the nearest CC bond at the hexagon gets expanded (Table S2). When the C atom is substituted by a N atom, the entire situation gets changed. The substitution of the sp-hybridized carbon by nitrogen gives rise to two types of NC bonds having bond lengths of 1.186 and 1.341 Å, respectively, and consequently, the chain lengths get shortened compared to graphyne. The NC bond length in the chain is slightly shorter than the estimated NC triple bond length92 and considered to exhibit a triple bond character. The NC bond length connecting the chain and ring is slightly greater than the pure NC double bond (1.270 Å) found in dimethylglyoxime,93 and smaller than the NC single bond found in methylamine;94 thus, it behaves as a NC double bond. To balance the contraction of chain length due to N substitution, the nearest CC bond at the ring gets expanded. The gradual increase of boron (nitrogen) concentration at the chain site 26795

DOI: 10.1021/acs.jpcc.6b07478 J. Phys. Chem. C 2016, 120, 26793−26806

Article

The Journal of Physical Chemistry C polarization in graphyne. Thus, our “graphyne with B at chain” systems are magnetic with a noticeable magnetic moment which increases with increasing doping concentration. In contrast, “graphyne with B atom at ring” structures are nonmagnetic, the same as “B-doped graphene”.99 The difference of the total energy between ferromagnetic (FM) and antiferromagnetic (AFM) states indicates that “graphyne with 1B at chain” is in the FM ground state but “graphyne with 2B at chain” and “graphyne with 3B at chain” are in the AFM ground state. For systems with B substituted at the chain site, the spin density difference as depicted in Figure 2 is mainly located

densification of charge. Thus, localization of charge occurs around the B−C bond. Moreover, the asymmetrical π bonds (due to asymmetric dumbbell-like orbitals) between B and C also create a localized state near the Fermi level and break the spin degeneracy. The spin splitting takes place in the vicinity of the Fermi level. Substitutional doping with boron gives rise to impurity states in both spin-up and spin-down bands, and the impurity bands which cross the Fermi level are well dispersed. The well-dispersed nature of the B-induced impurity bands increases the electrical conductivity of “B-doped graphyne”. When the sp2 hybridized carbon atom is replaced by a boron atom, π bonds break and the boron atom forms three sigma bonds with three neighboring carbon atoms. Surprisingly, this type of π bond breaking in “graphyne with B at ring” does not break the spin degeneracy and no spin splitting has been observed near the Fermi level. However, the creation of holes due to boron doping shifts the Fermi level downward and the impurity states cross the Fermi level. As a result, “graphyne with B at ring” makes a contribution to electrical activity. Substitutional N doping increases the number of electrons irrespective of the doping site and creates an impurity state above the Fermi level except “graphyne with 3N at chain” (Figure S2g). Thus, the impurity states generated by nitrogen in N-doped graphyne make a contribution to electrical conductivity and this contribution increases with increasing N concentration. Spin splitting of the band has not been observed for N-doped graphyne, since the difference in atomic radius between N and C is less compared to B and C, a small amount of local strain is generated, which cannot break the spin degeneracy. To get detailed information on a particular band, fat band101 studies (Figures S3−S6) have been done. For “graphyne with 1B at chain”, the fat band plot shows (Figure S3) that both the C and B atoms contribute to the energy band around the Fermi level. The dispersion less band at the edge of the Fermi level is basically contributed by the px−py bonding orbital of B as well as the C atom (whose hybridization is in between sp and sp2 as discussed earlier) at the chain and can be considered as a px−py π band (as depicted in the −COHP plot, Figures S7 and S8). Similar to pristine graphyne, the bonding pz π band is mostly contributed by the carbon atom and ranging from −5.620 eV to Fermi energy in the VB. Besides, we found notable bands with px and py characters between −2.620 and −1.600 eV. These px and py states mainly originated from the Cchain−Cchain bond (as depicted in the −COHP plot, Figure S8) and can be considered as bonding px−py π bands, because the sp C atom has “one σ and two π” character. The bands below −7.000 eV are the bonding σ band that arises from the px and py orbital contribution. In CB, the antibonding pz π* and px−py π* bands expanded from 0.000 to 11.780 eV and from 2.280 to 6.780 eV, respectively. Above 12.000 eV, the bands are basically σ* bands contributed by px and py orbitals. For “graphyne with 1N at chain” (Figure S4), the bonding pz π bands are ranging from −7.820 eV to Fermi energy while the bonding px−py π bands arising from the Cchain−Cchain bond lie between −3.320 to −1.000 eV and −8.000 to −4.000 eV (as depicted in the −COHP plot, Figures S9 and S10). Below −8.500 eV, the bands exhibit bonding σ character. In CB, the antibonding pz π* bands are ranging from Fermi energy to 10.680 eV and the antibonding px−py π* bands appear in between 1.670 and 3.300 eV due to the contribution of sp hybridized C and N atoms. Just above 11.100 eV, the antibonding σ* bands originating from px and py orbitals are observed. When B (N) is substituted at the ring (Figures S5 and S6) in the VB, the bonding pz π and

Figure 2. Top view of the spin density difference ρ = ρ↑ − ρ↓ for Bdoped graphyne: (a) graphyne with 1B at chain; (b) graphyne with 2B at chain; (c) graphyne with 3B at chain configurations.

around the BC bonds at the chain and a very small contribution from the carbon atoms in the ring and chain, which indicates that the total magnetic moment of the system is essentially contributed by the BC bonds at the chain with a small contribution from the C atoms of the hexagon and chain. For N-substitution, similar to “N-doped graphene”,99 the total magnetic moment is zero irrespective of the substitutional site. Thus, only by varying the concentration of the boron atom at the chain site, the magnetism of doped graphyne can be tuned. C. Band Structures. The spin-resolved band diagrams of graphyne with one dopant atom and multiple dopant atoms are depicted in Figures 3 and S1, respectively. The bands around the Fermi energy are significantly affected due to the presence of a B (N) atom because the doping by B (N) is similar to hole (electron) doping which results in a gradual shift of bands in order to accommodate an extra hole (electron). A similar type of downward/upward shift of bands due to the presence of B and N is also visible in B- or N-doped graphenes,48 where some bands cross the Fermi level. Doping with boron or nitrogen switched the conducting nature of pristine graphyne36 from semiconducting to metallic even when the concentration of the dopant atom is minimum. No splitting of the band for up and down spin channels has been observed for nitrogen doping (Figure S1c, d, g, h) irrespective of the substitution site which indicates that spin degeneracy does not break for N-doped graphyne and remains the same as pristine graphyne. Similar spin degeneracy has been observed in the case of “graphyne with B at ring”. For “graphyne with B at chain”, the bands corresponding to up and down spins are asymmetrical, which implies that the spin degeneracy has been broken. When the sp hybridized C atom is substituted by a B atom, the symmetry of the π bond between the substituted C atom and its adjacent sp hybridized C atom is broken. Hence, B doping at the chain breaks the π bond and leaves unpaired the pz electron of the C atom at the chain close to the doped B atom. The local strain100 along with quantum entrapment gives rise to strong localization and densification of charge. In addition, the difference in radius between B and C generates local strain74,100 in B-doped graphyne, which also produces strong localization and 26796

DOI: 10.1021/acs.jpcc.6b07478 J. Phys. Chem. C 2016, 120, 26793−26806

Article

The Journal of Physical Chemistry C

Figure 3. Band structure, projected density of states (PDOS), and −COHP plot of (a) graphyne with 1B at chain; (b) graphyne with 1B at ring; (c) graphyne with 1N at chain; (d) graphyne with 1N at ring.

px−py π bands appear in between −6.180 eV to Fermi energy (−7.200 eV to Fermi energy) and −2.180 to −1.170 eV (−3.300 to −2.090 eV), respectively. Below −7.000 eV

(−8.000 eV), the bands are bonding σ bands. In CB, antibonding pz π* and px−py π* bands appear in between Fermi energy to ∼11.820 eV (Fermi energy to ∼10.900 eV) 26797

DOI: 10.1021/acs.jpcc.6b07478 J. Phys. Chem. C 2016, 120, 26793−26806

Article

The Journal of Physical Chemistry C

concentration increases in the chain (Figure S2a and e), the contribution of the pz orbital of the C and B atoms is no longer symmetric for both spins; a little bit of asymmetry is observed. In addition, the asymmetric contribution of px and py orbitals of B and C to DOS also increased. Hence, the spin polarization increases with increasing B concentration. The TDOS and PDOS plots indicate that substitution of an sp2 hybridized C atom by a B atom (Figure 3b, II and Figure S2b, d) makes “graphyne with B at ring” completely nonmagnetic, as the contributions of both spins are symmetric and this nature is also observed in “B-doped graphene”.104 Similar to “B-doped graphene”,102,103 for “graphyne with B at ring”, the contribution of pz is highest for both C and B atoms around the Fermi level than the px and py orbitals, which indicates stronger hybridization between pz orbitals of B and C than px−px, py−py, and px−py hybridizations. The contributions of px and py orbitals gradually come closer to the Fermi level with increasing number of B atoms in the ring site as the hybridization between px−px, py−py, and px−py increases with increasing concentration. The N-doped graphyne acts as an ntype material the same as “N-doped graphene”102,103 due to the electron rich character of N which results in an increase of electrons in the system and causes an impurity state. In the case of “graphyne with N at ring”, the impurity state lies at the Fermi level and spreads over the valence and conduction bands. However, for “graphyne with N at chain”, this impurity state is located in the conduction band away from the Fermi level (Figure 3c, II and Figure S2c, g) except “graphyne with 3N at chain”. Increasing N concentration activates the contribution of px, py orbitals of N and C atoms near the Fermi level in the CB, indicating various degrees of orbital hybridization between the p orbitals. The comparison of PDOS of both types of dopant shows that the B and N atoms create impurity states in the vicinity of the Fermi level and enhance the carrier density by creating an additional charge carrier independent of the doping position. Additionally, the comparison of doping position shows that the substitution of an sp-hybridized C atom by a dopant activates the contribution of px, py and s orbitals near the Fermi level more significantly than when an sp2 hybridized C atom is replaced by the dopant. It further implies that the substitution of an sp-hybridized C atom by a dopant changed the hybridization more. Thus, the change of impurity position has a considerable effect on the electronic properties of doped graphyne. To complete the study of B- and N-doped graphyne with different dopant concentrations, it is important to analyze their nature of chemical bonding for further realization of the electronic structure. The bonding nature is characterized by crystal orbital Hamilton population (−COHP) analysis which presents a clear-cut view of orbital-pair interactions (Figure 3 and Figure S11). This technique partitions the band structure energy into bonding, nonbonding, and antibonding energy regions within a specified energy range, based on which it is possible to interpret the bonding situation in doped graphyne and the nature of energy states at the Fermi level. A positive value of −COHP stands for the bonding state, whereas a negative value represents the antibonding state. The −COHP diagram of “graphyne with B at chain” (Figure 3a, III and Figure S11a, e) shows that the energy states at the Fermi level are basically bonding states. At the Fermi level, it appears that Cring−Bchain and Cring−Cchain interactions are antibonding states, while Cchain−Bchain, Cchain−Cchain, and Cring−Cing interactions are bonding states, indicating that the interactions between Cchain and Bchain in the unit cell are stronger than those of Cring and

and 3.820 to 6.500 eV (1.900 to 6.800 eV), respectively. Moreover, the px−py π and px−py π* bands in “graphyne with B or N at ring” are attributed to a C−C bond at the chain only and no contribution of dopant is visible in the px−py π and px− py π* bands, because now the dopant B atom is almost sp3 hybridized which mainly contributes to σ and σ* bands, whereas N is sp2 hybridized which has σ + π character and contributes to pz π, pz π* and σ, σ* bands. Increasing dopant concentration moves px−py π and px−py π* bands closer to the Fermi level. D. Density of States. In pristine graphyne,36 the energy states near the Fermi level arise from the pz orbital of the C atom, and at the Fermi level, no energy state is observed which confirms its semiconducting nature. However, for N- or Bdoped graphyne, a significant amount of energy states are present at the Fermi level (Figures 3 and S2), confirming their metallic character which is also supported by our fat band study (Figures S3−S6). The electronic density of states near the valence band is getting increased with increasing dopant atom concentration. This is due to the fact that B tends to act as an electron donor in the vicinity of the C−B bond and N tends to act as an acceptor in the vicinity of the C−N bond. Substitution of a single C atom by a B atom shifts the Fermi level downward due to the electron deficient character of the B atom, and this type of shift is also observed in B-doped graphene.102,103 The shift is a little bit more for “graphyne with B at ring” than “graphyne with B at chain”. The opposite effect has been observed for N doping where a N atom shifts the Fermi level toward high energy and the shift is higher for “graphyne with N at chain” than “graphyne with N at ring”. A similar shift of Fermi level toward CB is also reported for N-doped graphene102,103 to accommodate the extra electron. The downward (upward) shift of Fermi level for “B (N)-doped graphyne” increases with increasing concentration irrespective of doping site. It should be mentioned that the shift is more in the case of N doping compared to B doping for the same doping site and concentration (see Table S3). The PDOS of doped graphyne not only sheds insights on the contribution of orbitals of constituting atoms but also focuses on the origin of spin polarization. Notably, except “graphyne with B at chain”, PDOS graphs indicate the nonmagnetic nature by exhibiting symmetrical DOS for up and down spin. For “graphyne with B at chain” (Figure 3a, II), the asymmetric nature of TDOS is only visible around the Fermi level, which implies the breaking of spin degeneracy around the Fermi level. From the PDOS of “graphyne with 1B at chain” (Figure 3a, II, ii and iii), it is clear that the spin degeneracy originates from the px and py orbitals of the B atom and the sp hybridized C atom in the chain nearest to it. The contributions of the pz orbitals of C and B atoms are symmetric for up and down spin, and hence, spin splitting is attributed to the px and py orbitals at and near the Fermi level (see Figure 3a, II, ii and iii). At the Fermi level, the PDOS of the C atom (Figure 3a, II, ii) shows that the contribution of the pz (px) orbital is highest for the majority (minority) spin channel. However, the PDOS curve of the B atom indicates that the contribution of px is highest for both spin channels at the Fermi level which further implies the presence of various degrees of orbital hybridization between p orbitals of B and the nearest C atom. It is evident from −COHP calculation that the orbital hybridization takes place between px−px, py−py, and px−py orbitals of B and the nearest C atom and is comparatively stronger than pz−pz, px−pz, and py−pz orbital hybridization (Figure S7). When the B 26798

DOI: 10.1021/acs.jpcc.6b07478 J. Phys. Chem. C 2016, 120, 26793−26806

Article

The Journal of Physical Chemistry C Bchain, since the Cring−Bchain bond length is larger than the Cchain−Bchain bond length, as shown by the structural parameter analysis. On the contrary, for “graphyne with N at chain”, Cring−Cring contributes antibonding states and Cchain−Cring contributes bonding states at and near the Fermi level. The Cring−Nchain bond contributes bonding states at and near the Fermi level, whereas the Cchain−Nchain bond contributes a small antibonding state at the Fermi level and strong antibonding states at the CB in the energy range 1.200−2.000 eV. The bonding character gets entirely changed when the sp 2 hybridized carbon atom is substituted by a B or N atom. For “graphyne with B at ring”, the energy state at the Fermi level is bonding in nature, contributed by Cchain−Bring, Cring−Bring, Cchain−Cchain, and Cring−Cring bonds. Only the Cring−Cchain interaction appears as an antibonding state at the Fermi level. With increasing concentration of B atom in the unit cell, above and below the Fermi level, the −COHP analysis of the Cring− Cchain bond shows that the antibonding contribution increases. In the case of “graphyne with N at ring”, the Cring−Nring and Cring−Cring appear as antibonding states at the Fermi level and the whole specified energy range in the CB. However, Cchain− Cchain, Nring−Cchain, and Cring−Cchain interactions contribute bonding states at the Fermi level. With increasing N concentration, the antibonding contribution of Cring−Nring and Cring−Cring is expanding below the Fermi level. In addition, the Nring−Cchain interaction also appears as the antibonding state at and near the Fermi level, indicating the instability of the N−C bond due to increasing concentration. The change in electronic properties originates from the electron redistribution due to doping. The substitution of the C atom by a B or N atom introduces an unpaired electron in graphyne, which breaks the original electron distribution. Accordingly, the localized electron distribution is modified and will enhance its chemical reactivity as well as electrical conductivity. In pristine graphyne, most of the electron density is located around the carbon triple bond, whereas the charge density around the C− C bond at the ring shows the existence of a covalent bond between the carbon atoms. In B (N)-doped graphyne, the electron density is reduced (enhance) around the doping site because of the electronegativity difference of B (N) and C. Because of the electronegativity (χ) difference between boron (χ = 2.04), carbon (χ = 2.55), and nitrogen (χ = 3.04), B−C and N−C bonds get polarized and the bond electrons shift toward C and N for B−C and N−C bond, respectively, which interprets the considerable ionicity of the B−C and N−C bonds. The charge density plot (Figure 4) of “B-doped graphyne” shows that the asymmetricity of the electron cloud along the B−C bond is larger in the case of “substitution at ring site” than “substitution at chain site”. Hence, the B−C bond at the ring is more ionic than the B−C bond at the chain. The strong bonding between B and C atoms at the ring, due to mixed ionic and covalent interactions, clarifies the B atom’s affinity toward the sp2 bonded C atom over the sp bonded C atom. Similarly, in the case of nitrogen doping, asymmetric charge accumulation is also observed between the sp bonded N atom and the sp2 bonded N atom. Comparison between “N at chain” and “N at ring” shows that the charge accumulation is higher around “N at chain” than around “N at ring”, which reflects the N atom’s affinity toward the sp bonded C site. E. Optical Properties. The real and imaginary parts of the dielectric function for B- or N-substituted graphyne are calculated by considering the interband transition as well as intraband transition dominated by the free-carrier response

Figure 4. Electron density plot of (a) graphyne with 3B at chain; (b) graphyne with 3N at chain; (c) graphyne with 3B at ring; (d) graphyne with 3N at ring.

which confirms that the systems are optically metallic. Figure 5 shows the variation of real and imaginary parts of the dielectric

Figure 5. Comparison of the dielectric function between interband and intraband transition of “graphyne with 1B at chain” for average electric field: (a) real part of dielectric function; (b) imaginary part of dielectric function.

function with energy for “graphyne with 1B at chain” considering an average electric field in which polarization is in all the directions. The magnitude of the static dielectric tensor (Figure 5a) is many times larger for intraband transition than interband transition. Near zero energy (up to 0.400 eV), only intraband transition dominates. The imaginary part of the dielectric function (Figure 5b) also confirms the contribution of intraband transition near zero energy. When we compare the static dielectric tensor for interband transition of B- or N-doped graphyne and BN-codoped graphyne, we see that ε1(0) of “Bor N-doped graphyne” is much greater than BN-codoped graphyne and its BN analogues (Table S4), which indicates that B- or N-doped graphyne may have a higher conductivity and carrier mobility than BN-doped graphyne and its BN analogues. The inclusion of intraband transition shows that ε1(0) of B- or N-doped graphyne is much greater than pristine graphyne 26799

DOI: 10.1021/acs.jpcc.6b07478 J. Phys. Chem. C 2016, 120, 26793−26806

Article

The Journal of Physical Chemistry C (Figure S12). Hence, “B- or N-doped graphyne” has a significant future in optoelectronic device application. Figure 6 and Figure 7 demonstrate the frequency-dependent absorption spectra and the real part of optical conductivity of

optical response is due to intraband transition and it is arising from the free carrier. However, the interband absorption occurs from the direct optical transition between the valence and conduction bands and governs the spectra in the UV region. It is evident from Figure 6a and b that, for boron-doped structures, the strong absorption peak is shifted toward lower energy with increasing dopant concentration; i.e., a red shift of strong peak is observed with increasing dopant concentration. In contrast, for nitrogen-doped structures, the strong peak exhibits a more or less blue shift toward high energy with increasing nitrogen concentration irrespective of chain or ring position except for “graphyne with 2N at chain”. For “graphyne with 2N at chain”, two nearly equal strong peaks are observed in the energy range 19.000−22.000 eV. The first peak is redshifted, while the second peak is blue-shifted with respect to “graphyne with 1N at chain”. The absorption coefficient is pronounced in a wide UV region, and these peaks are originating from interband transitions. Similar to absorption spectra, with increasing B concentration, the strong peak of optical conductivity (Figure 7) shows a red shift. However, in the case of N doping, “graphyne with 2N” shows a red shift, whereas “graphyne with 3N” shows a blue shift with respect to “graphyne with 1N”. The optical conductivity also displays prominent peaks in the UV region. The interband transitions are liable for the well-defined peaks of conductivity spectra in the UV region (Figure S14). At very low energy (up to 0.400 eV), only the intraband transition takes place. Another important optical quantity, the reflectance of “B- or N-doped graphyne”, is presented in Figure 8 and Figure S15, considering the average electric field to explain the reflectivity from the surface of “B- or N-doped graphyne”. Main peaks of reflectivity are arising below 5 eV (where the absorption is comparatively low), and for high energy, the reflectivity vanishes and the system becomes transparent. Moreover, the free carriers have a pronounced effect on reflectance only at low frequencies, making the reflectance value high. At 0.000−0.380 eV energy, the reflectivity for all “B- or N-doped graphyne” is equal to unity, irrespective of the doping site and concentration, attributed by the free carrier (intraband transition). That means in this energy range the reflectivity is 100%, so an electromagnetic wave cannot pass through it. The energy loss function (ELF) spectrum arising from interband and intraband transitions is contributed by plasmon excitations; thus, by analyzing ELF in connection with the dielectric function, all excitations can be recognized. The energy loss spectra of “B- or N-doped graphyne” are presented in Figure 9, where the peaks may correspond to plasmons and are arising due to the collective oscillations of free electrons or may be due to interband transitions from VB to CB. Table S5 presents the plasmon energy at which the ELF curve exhibits the highest peak. In the case of “graphyne with B at chain”, except “graphyne with 3B at chain” (Figure 9c), the sharp resonance peak is situated in the infrared region which corresponds to plasmons that are contributed by free carriers. The highest peak of “graphyne with 3B at chain” is located in the visible region (1.851 eV) and is mainly contributed by interband transition with a small contribution of free carriers. It is obvious from Table S5 that, with increasing boron concentration at the chain site, the plasmon frequency increases. That means there is a blue shift in plasmonic frequencies with the increasing dopant concentration. In contrast, the “graphyne with N at chain” exhibits a red shift of plasma frequency (from UV to infrared) with increasing

Figure 6. Comparison of absorption spectra with increasing dopant concentration: (a) graphyne with B at chain; (b) graphyne with B at ring; (c) graphyne with N at chain; (d) graphyne with N at ring. “GY” stands for graphyne.

Figure 7. Comparison of the real part of optical conductivity with increasing dopant concentration: (a) graphyne with B at chain; (b) graphyne with B at ring; (c) graphyne with N at chain; (d) graphyne with N at ring. “GY” stands for graphyne.

doped graphyne, respectively, considering the average electric field. Similar to graphene, the absorption spectra of “B- or Ndoped graphyne” exhibit a wide spectral range, extended from the infrared to the ultraviolet region of the electromagnetic spectrum. There are two types of contributions (Figure S13), namely, intraband and interband optical transitions, in the absorption process. In the infrared region (up to 0.400 eV), the 26800

DOI: 10.1021/acs.jpcc.6b07478 J. Phys. Chem. C 2016, 120, 26793−26806

Article

The Journal of Physical Chemistry C

Figure 8. Comparison of reflectivity between interband and intraband transition: (a) graphyne with 1B at chain; (b) graphyne with 2B at chain; (c) graphyne with 3B at chain; (d) graphyne with 1N at chain; (e) graphyne with 2N at chain; (f) graphyne with 3N at chain. “GY” stands for graphyne.

Figure 9. Comparison of the energy loss function due to interband and intraband transitions: (a) graphyne with 1B at chain; (b) graphyne with 2B at chain; (c) graphyne with 3B at chain; (d) graphyne with 1N at chain; (e) graphyne with 2N at chain; (f) graphyne with 3N at chain. “GY” stands for graphyne.

dopant concentration. For “graphyne with 3N at chain”, the highest peak corresponds to the energy of volume plasmons and is located in the infrared region, contributed by the free carrier, while for “graphyne with 1N at chain” the major peak corresponds to interband transitions and is located in the UV region (7.471 eV). The strong peak of ELF spectra of “graphyne with 2N at chain” is also situated in the UV region (7.441 eV). Similar to “graphyne B at ring”, the plasmon frequency of “graphyne with N at ring” is located at the infrared and visible

regions, respectively. Here also the major peak can be tuned from infrared to visible region with increasing dopant concentration. The strong peaks of the ELF spectrum, which are located in the infrared region, and originate from the collective oscillations of free electrons, correspond to the same energies at which reflectance (Figure 8) exhibits sharp reductions and the real part of the dielectric function passes from negative to positive values. The highest peak contributed by interband transition also corresponds to energy where reflectance falls, and the real part of the dielectric function also 26801

DOI: 10.1021/acs.jpcc.6b07478 J. Phys. Chem. C 2016, 120, 26793−26806

Article

The Journal of Physical Chemistry C

The positions of the highest peak of ε2 for parallel field and the possible interband transition105 of “B- or N-doped graphyne” are tabulated in Table S6. Since the substitutional doping in graphyne with B or N does not hamper its hexagonal symmetry, we can identify the dominant interband transition in different energy ranges on the basis of symmetry arguments and a projected band diagram (in Figures S3−S6). Same as pristine graphyne,32 σ, σ*, px−py π, and px−py π* states have even parity, whereas pz π and pz π* states have odd parity. For inplane polarization, the allowed transitions are taking place between the states with the same parity, but for out-of-plane polarization, transitions take place between the states with different parity. In parallel polarization, for all systems except “graphyne with B at chain”, the strong interband absorption from ∼0.500 to ∼8.000 eV comes due to pz π → pz π* transition, since the bands around the Fermi level are pz π (in VB) and pz π* (in CB). However, in the energy range 8.000− 10.000 eV, the px−py π → px−py π* transition dominates, and between 10.000 and 14.000 eV, the absorption is arising due to the px−py π → σ* and σ → px−py π* transitions. In the case of perpendicular polarization, the weak absorption below 9.000 eV is arising from px−pyπ → pz π* and pz π → px−py π* transitions, whereas over 9.000 eV the absorption mainly originates from σ → pzπ* and pz π → σ* transitions. For “graphyne with B at chain”, the strong interband absorption from ∼0.500 to ∼8.300 eV arises due to both pz π → pz π* and px−pyπ → px−py π* transitions. The absorption spectra depend on the polarization of the electric field vector. In the case of a parallel field, the inclusions of intraband contribution significantly change the optical spectra of “B- or N-doped graphyne”. On the other hand, the optical spectra of “B- or N-doped graphyne” remain unaltered for the perpendicular field, since absorption starts at high energy in this case. The modification of conducting nature due to the different types of electric field is clearly visible from the optical conductivity plot presented in Figure S16. The optical conductivity starts at nearly zero for the parallel electric field but starts with a gap for the perpendicular electric field, which confirms that all of these structures are optical semiconductors for the perpendicular electric field. It is noticed (Figure S17) that the reflectivity spectra of all doped graphynes are characterized by strong peaks under parallel polarization in the low energy region; however, for perpendicular polarization, strong peaks are located in the UV region. In comparison with pristine graphyne,32,62 “B- or N-doped graphyne” is more sensitive in the low-frequency region but less sensitive in the high-frequency region. In the case of “N-doped graphyne” with parallel polarization (Figure S18), the reflectance shows a strong peak (Rmax) within the infrared to visible region which increases with increasing concentration but “graphyne with 3N at chain”. A similar observation has been made on “N-doped graphene”.61 The only difference is that for “N-doped graphyne” Rmax is more extended toward the visible region (even extended up to the UV region for “graphyne with N at chain”) than “N-doped graphene” which in turn indicates the illumination possibility of “N-doped graphyne” with visible light of a wide energy range. This feature is expected as N-doping increases electrons in the system. On the other hand, unlike “Bdoped graphene”,61 the Rmax value of “B-doped graphyne” also increases with increasing concentration and lies within the infrared to visible region.

falls but does not pass from negative to positive values except “graphyne with 3B at chain or ring”. To tune the opto-metallic properties of our B- or N-doped graphyne, we have considered three types of electric fields such as (a) average electric field, (b) parallel electric field, and (c) perpendicular electric field. As intraband transition due to free carrier dominates only in the low energy region, it does not affect the optical response for the perpendicular field. A comparison of dielectric functions for interband transition only with three types of electric fields is presented in Figure 10. We

Figure 10. Variation of the imaginary part of the dielectric function with respect to photon energy for average, parallel, and perpendicular electric field: (a) graphyne with 1B at chain; (b) graphyne with 1B at ring; (c) graphyne with 1N at chain; (d) graphyne with 1N at ring. “GY” stands for graphyne. The average, parallel, and perpendicular electric fields are represented by black, red, and green lines, respectively.

have found the same strong anisotropic nature in the dielectric function in the low-energy region for our B- or N-doped graphyne as we observed in pristine graphyne.67 For example, in “graphyne with 1B at chain” with the parallel electric field, the intensity of the dielectric function appears high in the energy range 0.020−10.000 eV which ultimately decreases with increasing photon energy. However, for the perpendicular electric field, the dielectric function is notable above 12.900 eV and below this value the absorption is negligible. At low energy, the parallel electric field dominates while at high energy the perpendicular electric field dominates. Furthermore, the direction of the electric field changes the value of the static dielectric tensor, as listed in Table S4. The static dielectric tensor (ε1(0)) has a larger value for the parallel electric field than the perpendicular electric field, the same as B- or N-doped graphene and greater than pristine graphene.61 Moreover, same as “B- or N-doped graphene”, the ε1(0) does not change appreciably with increasing concentration for perpendicular polarization and significantly changes for parallel polarization. Unlike “B-doped graphene”, ε1(0) increases with increasing concentration for “graphyne with B at chain”. For the perpendicular electric field, the optical response starts with a gap and optical responses show the optical semiconductor properties for these structures. 26802

DOI: 10.1021/acs.jpcc.6b07478 J. Phys. Chem. C 2016, 120, 26793−26806

The Journal of Physical Chemistry C





CONCLUSIONS

We have studied the structural, electronic, magnetic, and optical properties of “B- or N-doped graphyne” using density functional theory. The cohesive energies of these doped graphynes are comparable with the experimental cohesive energy of graphite, establishing their stability. Our study reveals that the B atom prefers to substitute the sp2-hybridized C atom in the hexagon while the favorable substitution site for the N atom is the sp-hybridized C atom in the chain. Both types of dopants produce a very small local structural deformation around the doping site. The presence of a dopant atom in the unit cell induced energy levels at or close to the Fermi level and converted the semiconductor graphyne into metallic in nature. All “N-doped graphyne” and “B-doped graphyne” are nonmagnetic, but in the case of “graphyne with B at chain”, a highly localized impurity state is created in the Fermi level which enhances the electrical activity and introduces the magnetism in “graphyne with B at chain”. The inclusion of a single boron atom at the chain site of graphyne introduces ferromagnetism, while the presence of two or three boron atoms at the chain site gives rise to antiferromagnetism. The spin density difference explores that the magnetism is basically attributed by the C−B bond at the chain. The PDOS gives the evidence that spin degeneracy is broken only at and near the Fermi level, not far away from the Fermi level, and the spin splitting is attributed by the px and py orbitals at and near the Fermi level for a single boron atom. Moreover, spin polarization increases with increasing boron concentration. Similar to graphyne, “B- or N-doped graphyne” shows the same type of optical anisotropy for parallel and perpendicular electric fields. It is found that, below 0.400 eV, the optical response is dominated by the intraband transition which arises from free carriers. The comparison of the static dielectric tensor predicts a higher conductivity and carrier mobility of B- or N-doped graphyne than pristine graphyne, and it’s BN analogues which in turn enhance the possibility of their use in optoelectronic devices. Like graphene, “B- or N-doped graphyne” shows strong absorption in the UV region and also exhibits absorption in a wide spectral range, extended from the infrared to the ultraviolet region. The optical conductivity for perpendicular polarization exhibits the optically semiconducting nature of these structures. From ELF curves and corresponding reflectance, we have derived plasmon energies for all B- or Ndoped structures. The reflectivity for all structures is 100% at extremely low energy (below 0.380 eV) irrespective of the doping site and concentration, showing their capability of reflecting electromagnetic waves in this range. In addition, for electric field polarized perpendicular to the plane, reflectivity starts with a gap and order of magnitude less than the parallel field, inferring the capability of transmitting electromagnetic waves.



Article

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Phone: +919401542687. Author Contributions

The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS U.S. wishes to acknowledge the assistance of computational resources from Sharcnet and Compute Canada. This research is supported by Assam University, Silchar, India.



REFERENCES

(1) Rinzler, A. G.; Hafner, J. H.; Nikolaev, P.; Lou, L.; Kim, S. G.; Tomanek, D.; Nordlander, P.; Colbert, D. T.; Smalley, R. E. Unraveling Nanotubes: Field Emission from an Atomic Wire. Science 1995, 269, 1550−1553. (2) Li, G.; Li, Y.; Qian, X.; Liu, H.; Lin, H.; Chen, N.; Li, Y. Construction of Tubular Molecule Aggregations of Graphdiyne for Highly Efficient Field Emission. J. Phys. Chem. C 2011, 115, 2611− 2615. (3) Li, T. C.; Lu, S. P. Quantum Conductance of Graphene Nanoribbons with Edge Defects. Phys. Rev. B: Condens. Matter Mater. Phys. 2008, 77, 085408. (4) Williams, J. R.; Abanin, D. A.; DiCarlo, L.; Levitov, L. S.; Marcus, C. M. Quantum Hall Conductance of Two-terminal Graphene Devices. Phys. Rev. B: Condens. Matter Mater. Phys. 2009, 80, 045408. (5) Zhou, J.; Sun, Q.; Wang, Q.; Jena, P. High-temperature Superconductivity in Heavily N- or B-doped Graphene. Phys. Rev. B: Condens. Matter Mater. Phys. 2015, 92, 064505. (6) Faye, J. P. L.; Sahebsara, P.; Sénéchal, D. Chiral Triplet Superconductivity on the Graphene Lattice. Phys. Rev. B: Condens. Matter Mater. Phys. 2015, 92, 085121. (7) Nandkishore, R.; Chubukov, A. V. Interplay of Superconductivity and Spin-density-wave Order in Doped Graphene. Phys. Rev. B: Condens. Matter Mater. Phys. 2012, 86, 115426. (8) Ludbrook, B. M.; Levy, G.; Nigge, P.; Zonno, M.; Schneider, M.; Dvorak, D. J.; Veenstra, C. N.; Zhdanovich, S.; Wong, D.; Dosanjh, P.; Straßer, C.; Stohr, A.; Forti, S.; Ast, C. R.; Starke, U.; Damascelli, A. Evidence for Superconductivity in Li-decorated Monolayer Graphene. Proc. Natl. Acad. Sci. U. S. A. 2015, 112, 11795−11799. (9) Gan, L. Y.; Zhang, Q.; Guo, C. S.; Schwingenschlög, U.; Zhao, Y. Two-Dimensional MnO2/Graphene Interface: Half-Metallicity and Quantum Anomalous Hall State. J. Phys. Chem. C 2016, 120, 2119− 2125. (10) Dutta, S.; Pati, S. K. Half-Metallicity in Undoped and Boron Doped Graphene Nanoribbons in the Presence of Semilocal ExchangeCorrelation Interactions. J. Phys. Chem. B 2008, 112, 1333−1335. (11) Park, W.; Guo, Y.; Li, X.; Hu, J.; Liu, L.; Ruan, X.; Chen, Y. P. High-performance Thermal Interface Material Based on Few-layer Graphene Composite. J. Phys. Chem. C 2015, 119, 26753−26759. (12) Tan, X.; Shao, H.; Hu, T.; Liu, G.; Jiang, J.; Jiang, H. High Thermoelectric Performance in Two-dimensional Graphyne Sheets Predicted by First-principles Calculations. Phys. Chem. Chem. Phys. 2015, 17, 22872−81. (13) Xiao, Y.; Wang, W.; Lin, T.; Chen, X.; Zhang, Y.; Yang, J.; Wang, Y.; Zhou, Z. Largely Enhanced Thermal Conductivity and High Dielectric Constant of Poly(vinylidene fluoride)/Boron Nitride Composites Achieved by Adding a Few Carbon Nanotubes. J. Phys. Chem. C 2016, 120, 6344−6355.

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpcc.6b07478. Detailed information on electronic and optical properties of B- or N-doped graphyne and lists of bond lengths, magnetic moments, lattice vectors, and possible interband transitions of B- or N-doped graphyne (PDF) 26803

DOI: 10.1021/acs.jpcc.6b07478 J. Phys. Chem. C 2016, 120, 26793−26806

Article

The Journal of Physical Chemistry C (14) Hu, M.; Jing, Y.; Zhang, X. Low Thermal Conductivity of Graphyne Nanotubes from Molecular Dynamics Study. Phys. Rev. B: Condens. Matter Mater. Phys. 2015, 91, 155408. (15) Ö zçelik, V. O.; Ciraci, S. Size Dependence in the Stabilities and Electronic Properties of α-Graphyne and Its Boron Nitride Analogue. J. Phys. Chem. C 2013, 117, 2175−2182. (16) Shi, G.; Hanlumyuang, Y.; Liu, Z.; Gong, Y.; Gao, W.; Li, B.; Kono, J.; Lou, J.; Vajtai, R.; Sharma, P.; Ajayan, P. M. Boron Nitride− Graphene Nanocapacitor and the Origins of Anomalous SizeDependent Increase of Capacitance. Nano Lett. 2014, 14, 1739−1744. (17) Peymanirad, F.; Neek-Amal, M.; Beheshtian, J.; Peeters, F. M. Graphene-Silicene Bilayer: A Nanocapacitor with Permanent Dipole and Piezoelectricity Effect. Phys. Rev. B: Condens. Matter Mater. Phys. 2015, 92, 155113. (18) Subramaniam, C. K.; Maiyalagan, T. Double Layer Energy Storage in Graphene - a Study. Micro Nanosyst. 2012, 4, 180−185. (19) Wang, Y.; Chen, X.; Ye, W.; Wu, Z.; Han, Y.; Han, T.; He, Y.; Cai, Y.; Wang, N. Side-Gate Modulation Effects on High-quality BNGraphene-BN Nanoribbon Capacitors. Appl. Phys. Lett. 2014, 105, 243507−4. (20) Zhan, C.; Neal, J.; Wu, J.; Jiang, D. Quantum Effects on the Capacitance of Graphene-Based Electrodes. J. Phys. Chem. C 2015, 119, 22297−22303. (21) Mousavi-Khoshdel, M.; Targholi, E.; Momeni, M. J. FirstPrinciples Calculation of Quantum Capacitance of Codoped Graphenes as Supercapacitor Electrodes. J. Phys. Chem. C 2015, 119, 26290−26295. (22) Berdiyorov, G.; Bahlouli, H.; Peeters, F. M. Effect of Substitutional Impurities on the Electronic Transport Properties of Graphene. Phys. E 2016, 84, 22−26. (23) Lherbier, V.; Liang, L.; Charlier, J. C.; Meunier, V. Charge Carrier Transport and Separation in Pristine and Nitrogen-Doped Graphene Nanowiggle Heterostructures. Carbon 2015, 95, 833−842. (24) Poetschke, M.; Rocha, C. G.; Foa Torres, L. E. F.; Roche, S.; Cuniberti, G. Modeling Graphene-Based Nanoelectromechanical Devices. Phys. Rev. B: Condens. Matter Mater. Phys. 2010, 81, 193404. (25) Nemec, N.; Tománek, D.; Cuniberti, G. Modeling Extended Contacts for Nanotube and Graphene Devices. Phys. Rev. B: Condens. Matter Mater. Phys. 2008, 77, 125420. (26) Zhang, L.; Xia, Z. Mechanisms of Oxygen Reduction Reaction on Nitrogen-Doped Graphene for Fuel Cells. J. Phys. Chem. C 2011, 115, 11170−11176. (27) 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. (28) Sen, D.; Thapa, R.; Chattopadhyay, K. K. Rules of Boron− Nitrogen Doping in Defect Graphene Sheets: A First-Principles Investigation of Band-Gap Tuning and Oxygen Reduction Reaction Catalysis Capabilities. ChemPhysChem 2014, 15, 2542−2549. (29) Lee, W. J.; Maiti, U. N.; Lee, J. M.; Lim, J.; Han, T. H.; Kim, S. O. Nitrogen-Doped Carbon Nanotubes and Graphene Composite Structures for Energy and Catalytic Applications. Chem. Commun. 2014, 50, 6818−6830. (30) Yang, L.; Jiang, S.; Zhao, Y.; Zhu, L.; Chen, S.; Wang, X.; Wu, Q.; Ma, J.; Ma, Y.; Hu, Z. Boron-Doped Carbon Nanotubes as MetalFree Electrocatalysts for the Oxygen Reduction Reaction. Angew. Chem. 2011, 123, 7270−7273. (31) Baughman, R. H.; Eckhardt, H.; Kertesz, M. Structure-Property Predictions for New Planar forms of Carbon: Layered Phases Containing sp2 and sp Atoms. J. Chem. Phys. 1987, 87, 6687−6699. (32) Kang, J.; Li, J.; Wu, F.; Li, S.; Xia, J. Elastic, Electronic, and Optical Properties of Two-Dimensional Graphyne Sheet. J. Phys. Chem. C 2011, 115, 20466−20470. (33) Narita, N.; Nagai, S.; Suzuki, S.; Nakao, K. Optimized Geometries and Electronic Structures of Graphyne and Its Family. Phys. Rev. B: Condens. Matter Mater. Phys. 1998, 58, 11009−11014. (34) Zhou, J.; Lv, K.; Wang, Q.; Chen, X. S.; Sun, Q.; Jena, P. Electronic Structures and Bonding of Graphyne Sheet and Its BN Analog. J. Chem. Phys. 2011, 134, 174701−5.

(35) Li, Y.; Xu, L.; Liua, H.; Li, Y. Graphdiyne and Graphyne: From Theoretical Predictions to Practical Construction. Chem. Soc. Rev. 2014, 43, 2572−2586. (36) Singh, N. B.; Bhattacharya, B.; Sarkar, U. A First Principle Study of Pristine and BN-Doped Graphyne Family. Struct. Chem. 2014, 25, 1695−1710. (37) Sarkar, U.; Bhattacharya, B.; Seriani, N. First Principle Study of Sodium Decorated Graphyne. Chem. Phys. 2015, 461, 74−80. (38) Li, G.; Li, Y.; Guo, H.; Li, Y.; Zhu, D. Architecture of Graphdiyne Nanoscale Films. Chem. Commun. 2010, 46, 3256−3258. (39) Qian, X.; Ning, Z.; Li, Y.; Liu, H.; Ouyang, C.; Chen, Q.; Li, Y. Construction of Graphdiyne Nanowires with High-Conductivity and Mobility. Dalton Trans 2012, 41, 730−733. (40) Du, H.; Deng, Z.; Lü, Z.; Yin, Y.; Yu, L. L.; Wu, H.; Chen, Z.; Zou, Y.; Wang, Y.; Liu, H.; Li, Y. The Effect of Graphdiyne Doping on the Performance of Polymer Solar Cells. Synth. Met. 2011, 161, 2055− 2057. (41) Wang, S.; Yi, L.; Halpert, J. E.; Lai, X.; Liu, Y.; Cao, H.; Yu, R.; Wang, D.; Li, Y. A Novel and Highly Efficient Photocatalyst Based on P25-Graphdiyne Nanocomposite. Small 2012, 8, 265−271. (42) Li, G. X.; Li, Y.; Qian, X.; Liu, H. B.; Lin, H. W.; Chen, N.; Li, Y. Construction of Tubular Molecule Aggregations of Graphdiyne for Highly Efficient Field Emission. J. Phys. Chem. C 2011, 115, 2611− 2615. (43) Kehoe, J. M.; Kiley, J. H.; English, J. J.; Johnson, C. A.; Petersen, R. C.; Haley, M. M. Carbon Networks Based on Dehydrobenzoannulenes. 3. Synthesis of Graphyne Substructures1. Org. Lett. 2000, 2, 969−972. (44) Haley, M. M. Synthesis and Properties of Annulenic Subunits of Graphyne and Graphdiyne Nanoarchitectures. Pure Appl. Chem. 2008, 80, 519−532. (45) Johnson, C. A.; Lu, Y.; Haley, M. M. Carbon Networks Based on Benzocyclynes. 6. Synthesis of Graphyne Substructures via Directed Alkyne Metathesis. Org. Lett. 2007, 9, 3725−3728. (46) Narita, N.; Nagai, S.; Suzuki, S. Potassium Intercalated Graphyne. Phys. Rev. B: Condens. Matter Mater. Phys. 2001, 64, 245408. (47) Mu, Y.; Li, S. D. Multiple Dirac Cones in BN Co-doped βGraphyne. J. Mater. Chem. C 2016, 4, 7339−7344. (48) 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. (49) Imran Jafri, R.; Rajalakshmi, N.; Ramaprabhu, S. Nitrogen Doped Graphene Nanoplatelets as Catalyst Support for Oxygen Reduction Reaction in Proton Exchange Membrane Fuel Cell. J. Mater. Chem. 2010, 20, 7114−7117. (50) Gao, Y.; Hu, G.; Zhong, J.; Shi, Z.; Zhu, Y.; Su, D. S.; Wang, J.; Bao, X.; Ma, D. Nitrogen-Doped sp2.2-Hybridized Carbon as a Superior Catalyst for Selective Oxidation. Angew. Chem., Int. Ed. 2013, 52, 2109−2113. (51) Lv, R.; dos Santos, M. C.; Antonelli, C.; Feng, S.; Fujisawa, K.; Berkdemir, A.; Cruz-Silva, R.; Elías, A. L.; Perea-Lopez, N.; LópezUrías, F.; Terrones, H.; Terrones, M. Large-Area Si-Doped Graphene: Controllable Synthesis and Enhanced Molecular Sensing. Adv. Mater. 2014, 26, 7593−7599. (52) Zou, Y.; Li, F.; Zhu, Z. H.; Zhao, M. W.; Xu, X. G.; Su, X. Y. An Ab Initio Study on Gas Sensing Properties of Graphene and Si-Doped Graphene. Eur. Phys. J. B 2011, 81, 475−479. (53) Ao, Z. M.; Yang, J.; Li, S.; Jiang, Q. Enhancement of CO Detection in Al Doped Graphene. Chem. Phys. Lett. 2008, 461, 276− 279. (54) Zhang, S. J.; Lin, S. S.; Li, X. Q.; Liu, X. Y.; Wu, H. A.; Xu, W. L.; Wang, P.; Wu, Z. Q.; Zhong, H. K.; Xu, Z. J. Opening the Band Gap of Graphene Through Silicon Doping for the Improved Performance of Graphene/GaAs Heterojunction Solar Cells. Nanoscale 2016, 8, 226−232. (55) Lherbier, A.; Botello-Mendez, A. R.; Charlier, J. C. Electronic and Transport Properties of Unbalanced Sublattice N-Doping in Graphene. Nano Lett. 2013, 13, 1446−1450. 26804

DOI: 10.1021/acs.jpcc.6b07478 J. Phys. Chem. C 2016, 120, 26793−26806

Article

The Journal of Physical Chemistry C (56) Zhang, L.; Xia, Z. Mechanisms of Oxygen Reduction Reaction on Nitrogen-Doped Graphene for Fuel Cells. J. Phys. Chem. C 2011, 115, 11170−11176. (57) 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. (58) Sforzini, J.; Hapala, P.; Franke, M.; Straaten, G. V.; Stöhr, A.; Link, S.; Soubatch, S.; Jelínek, P.; Lee, T. L.; Starke, U.; Švec, M.; Bocquet, F. C.; Tautz, F. S. Structural and Electronic Properties of Nitrogen-Doped Graphene. Phys. Rev. Lett. 2016, 116, 126805. (59) Zhang, F.; Wang, Z.; Wang, D.; Wu, Z.; Wang, S.; Xu, X. Nonlinear Optical Effects in Nitrogen−Doped Graphene. RSC Adv. 2016, 6, 3526−3531. (60) 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. (61) Nath, P.; Chowdhury, S.; Sanyal, D.; Jana, D. Ab-initio Calculation of Electronic and Optical Properties of Nitrogen and Boron Doped Graphene Nanosheet. Carbon 2014, 73, 275−282. (62) Lee, W. J.; Maiti, U. N.; Lee, J. M.; Lim, J.; Han, T. H.; Kim, S. O. Nitrogen-Doped Carbon Nanotubes and Graphene Composite Structures for Energy and Catalytic Applications. Chem. Commun. 2014, 50, 6818−6830. (63) Yang, L.; Jiang, S.; Zhao, Y.; Zhu, L.; Chen, S.; Wang, X.; Wu, Q.; Ma, J.; Ma, Y.; Hu, Z. Boron-Doped Carbon Nanotubes as MetalFree Electrocatalysts for the Oxygen Reduction Reaction. Angew. Chem. 2011, 123, 7270−7273. (64) Zheng, Y.; Liu, J.; Liang, J.; Jaroniec, M.; Qiao, S. Z. Graphitic Carbon Nitride Materials: Controllable Synthesis and Applications in Fuel Cells and Photocatalysis. Energy Environ. Sci. 2012, 5, 6717−6731. (65) Shen, B.; Chen, J.; Yan, X.; Xue, Q. Synthesis of Fluorine-Doped Multi-layered Graphene Sheets by Arc-discharge. RSC Adv. 2012, 2, 6761−6764. (66) Shaikjee, A.; Coville, N. J. The Synthesis, Properties and Uses of Carbon Materials with Helical Morphology. J. Adv. Res. 2012, 3, 195− 223. (67) Bhattacharya, B.; Singh, N. B.; Sarkar, U. Pristine and BN Doped Graphyne Derivatives for UV Light Protection. Int. J. Quantum Chem. 2015, 115, 820−829. (68) Jafari, M.; Asadpour, M.; Majelan, N. A.; Faghihnasiri, M. Effect of Boron and Nitrogen Doping on Electro-Optical Properties of Armchair and Zigzag Graphyne Nanoribbons. Comput. Mater. Sci. 2014, 82, 391−398. (69) Bhattacharya, B.; Singh, N. B.; Mondal, R.; Sarkar, U. Electronic and Optical Properties of Pristine and Boron−Nitrogen Doped Graphyne Nanotubes. Phys. Chem. Chem. Phys. 2015, 17, 19325− 19341. (70) Hwang, H. J.; Kwon, Y.; Lee, H. Thermodynamically Stable Calcium-Decorated Graphyne as a Hydrogen Storage Medium. J. Phys. Chem. C 2012, 116, 20220−20224. (71) Li, C.; Li, J.; Wu, F.; Li, S. S.; Xia, J. B.; Wang, L. W. High Capacity Hydrogen Storage in Ca Decorated Graphyne: A FirstPrinciples Study. J. Phys. Chem. C 2011, 115, 23221−23225. (72) Wang, Y. S.; Yuan, P. F.; Li, M.; Sun, Q.; Jia, Y. CalciumDecorated Graphyne Nanotubes as Promising Hydrogen Storage Media: A First-Principles Study. J. Solid State Chem. 2013, 197, 323− 328. (73) Hwang, H. J.; Koo, J.; Park, M.; Park, N.; Kwon, Y.; Lee, H. Multilayer Graphynes for Lithium Ion Battery Anode. J. Phys. Chem. C 2013, 117, 6919−6923. (74) Yun, J.; Zhang, Z.; Yan, J.; Zhao, W.; Xu, M. First-principles Study of B or Al-Doping Effect on the Structural, Electronic Structure and Magnetic Properties of γ-Graphyne. Comput. Mater. Sci. 2015, 108, 147−152. (75) Das, B. K.; Sen, D.; Chattopadhyay, K. K. Implications of Boron Doping on Electrocatalytic Activities of Graphyne and Graphdiyne Families: A First Principles Study. Phys. Chem. Chem. Phys. 2016, 18, 2949−2958.

(76) Kong, X.; Chen, Q.; Sun, Z. The Positive Influence of BoronDoped Graphyne on Surface Enhanced Raman Scattering with Pyridine as the Probe Molecule and Oxygen Reduction Reaction in Fuel Cells. RSC Adv. 2013, 3, 4074−4080. (77) Das, B. K.; Sen, D.; Chattopadhyay, K. K. Nitrogen Doping in Acetylene Bonded Two Dimensional Carbon Crystals: Ab-initio Forecast of Electrocatalytic Activities Vis-à-vis Boron Doping. Carbon 2016, 105, 330−339. (78) Chen, X.; Qiao, Q.; An, L.; Xia, D. Why Do Boron and Nitrogen Doped α- and γ-Graphyne Exhibit Different Oxygen Reduction Mechanism? A First-Principles Study. J. Phys. Chem. C 2015, 119, 11493−11498. (79) Zhang, S.; Du, H.; He, J.; Huang, C.; Liu, H.; Cui, G.; Li, Y. Nitrogen-Doped Graphdiyne Applied for Lithium-Ion Storage. ACS Appl. Mater. Interfaces 2016, 8, 8467−8473. (80) Lu, R.; Rao, D.; Meng, Z.; Zhang, X.; Xu, G.; Liu, Y.; Kan, E.; Xiao, C.; Deng, K. Boron-substituted Graphyne as a Versatile Material with High Storage Capacities of Li and H2: a Multiscale Theoretical Study. Phys. Chem. Chem. Phys. 2013, 15, 16120−16126. (81) Soler, J. M.; Artacho, E.; Gale, J. D.; García, A.; Junquera, J.; Ordejón, P.; Portal, D. S. The SIESTA Method for Ab Initio Order-N Materials Simulation. J. Phys.: Condens. Matter 2002, 14, 2745−2779. (82) Troullier, N.; Martins, J. A Straightforward Method for Generating Soft Transferable Pseudopotentials. Solid State Commun. 1990, 74, 613−616. (83) Das, R.; Chowdhury, S.; Jana, D. A First Principles Approach to Magnetic and Optical Properties in Single-layer Graphene Sandwiched Between Boron Nitride Monolayers. Mater. Res. Express 2015, 2, 075601. (84) Cao, X.; Li, Y.; Cheng, X.; Zhang, Y. Structural Analogues of Graphyne Family: New Types of Boron Nitride Sheets with Wide Band Gap and Strong UV Absorption. Chem. Phys. Lett. 2011, 502, 217−221. (85) Rashkeev, S. N.; Uspenskii, Y. A.; Mazin, I. I. Optical properties of transition metals at infrared frequencies. Sov. Phys. JETP 1985, 61, 5. (86) Shahrokhi, M.; Naderi, S.; Fathalian, A. A. Ab Initio Calculations of Optical Properties of B2C Graphene Sheet. Solid State Commun. 2012, 152, 1012−1017. (87) Shuai, Z.; Wang, D.; Peng, Q.; Geng, H. Computational Evaluation of Optoelectronic Properties for Organic/Carbon Materials. Acc. Chem. Res. 2014, 47, 3301−9. (88) Mahan, G. D. Many particle physics; Plenum Press: New York, 1990. (89) Dressel, M.; Gruner, G. Electrodynamics of solids; Cambridge University Press: Cambridge, U.K., 2002. (90) Putz, S.; Gmitra, M.; Fabian, J. Optical Conductivity of Hydrogenated Graphene from First Principles. Phys. Rev. B: Condens. Matter Mater. Phys. 2014, 89, 035437. (91) Olmstead, M. M.; Power, P. P.; Weese, K. J.; Doedens, R. J. Isolation and X-ray Crystal Structure of the Boron Methylidenide ion [Mes2BCH2]- (Mes = 2,4,6-Me3C6H2): A Boron-Carbon Double Bonded Alkene Analog. J. Am. Chem. Soc. 1987, 109, 2541−2542. (92) Majidi, R. Effect of Doping on the Electronic Properties of Graphyne. Nano 2013, 8, 1350060. (93) Merritt, L. L.; Lanterman, E. The Crystal Structure of Dimethylglyoxime. Acta Crystallogr. 1952, 5, 811−817. (94) Lide, D. R. Structure of the Methylamine Molecule. I. Microwave Spectrum of CD3ND2. J. Chem. Phys. 1957, 27, 343−352. (95) 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. (96) Jiang, X.; Arhammar, C.; Liu, P.; Zhao, J.; Ahuja, R. The R3Carbon Allotrope: A Pathway Towards Glassy Carbon Under High Pressure. Sci. Rep. 2013, 3, 1877. (97) Kittel, C. Introduction to Solid State Physics; Wiley: Singapore, 1996. (98) Ren, X.; Wang, B.; Zhu, J.; Liu, J.; Zhang, W.; Wen, Z. The Doping Effect on the Catalytic Activity of Graphene for Oxygen 26805

DOI: 10.1021/acs.jpcc.6b07478 J. Phys. Chem. C 2016, 120, 26793−26806

Article

The Journal of Physical Chemistry C Evolution Reaction in a Lithium−air Battery: a First-principles Study. Phys. Chem. Chem. Phys. 2015, 17, 14605−14612. (99) Wu, M.; Cao, C.; Jiang, J. Z. Light Non-metallic Atom (B, N, O and F)-doped Graphene: a First-principles Study. Nanotechnology 2010, 21, 505202. (100) Zheng, W. T.; Sun, C. Q. Underneath the Fascinations of Carbon Nanotubes and Graphene Nanoribbons. Energy Environ. Sci. 2011, 4, 627−655. (101) Flores, M.; Cisternas, E.; Correa, J. D.; Vargas, P. Moiré Patterns on STM Images of Graphite Induced by Rotations of Surface and Subsurface Layers. Chem. Phys. 2013, 423, 49−54. (102) Zheng, B.; Hermet, P.; Henrard, L. Scanning Tunneling Microscopy Simulations of Nitrogen- and Boron- Doped Graphene and Single-Walled Carbon Nanotubes. ACS Nano 2010, 4, 4165− 4173. (103) Mukherjee, S.; Kaloni, T. P. Electronic Properties of Boronand Nitrogen-doped Graphene: a First Principles Study. J. Nanopart. Res. 2012, 14, 1059. (104) Faccio, R.; Werner, L. F.; Pardo, H.; Goyenola, C.; Ventura, O. N.; Mombru, A. W. Electronic and Structural Distortions in Graphene Induced by Carbon Vacancies and Boron Doping. J. Phys. Chem. C 2010, 114, 18961−18971. (105) Kumar, A.; Kumar, A.; Ahluwalia, P. K. Ab Initio Study of Structural, Electronic and Dielectric Properties of Free Standing Ultrathin Nanowires of Noble Metals. Phys. E 2012, 46, 259−269.

26806

DOI: 10.1021/acs.jpcc.6b07478 J. Phys. Chem. C 2016, 120, 26793−26806