Patterned Partially Hydrogenated Graphene (C4H) and Its One

Jan 30, 2012 - Department of Chemistry, Institute for Functional Nanomaterials, University of Puerto Rico, Rio Piedras Campus, San Juan, Puerto Rico 0...
0 downloads 0 Views 4MB Size
Article pubs.acs.org/JPCC

Patterned Partially Hydrogenated Graphene (C4H) and Its OneDimensional Analogues: A Computational Study Yafei Li and Zhongfang Chen* Department of Chemistry, Institute for Functional Nanomaterials, University of Puerto Rico, Rio Piedras Campus, San Juan, Puerto Rico 00931 S Supporting Information *

ABSTRACT: By means of density functional theory (DFT) computations, we systematically studied the structural and electronic properties of the experimentally just achieved new two-dimensional (2D) hydrocarbon  the patterned partially hydrogenated graphene with formula C4H (Adv. Mater 2011, 23, 4497), and in particular its one-dimensional (1D) analogues. The C4H layer is a stable 2D crystal featured with periodic Clar sextet aromatic rings and is semiconducting with a wide band gap; however, this single-sided patterned partially hydrogenated C4H layer can only be obtained when the possibility of double-sided hydrogenation is excluded, since the doublesided graphane-embedded structure is energetically more favorable. The 1D C4H nanotubes, rolled up by the C4H layer, exhibit excellent thermodynamic properties and all have a wide band gap regardless of the tube diameter and chirality. In contrast, cutting the C4H layer into 1D C4H nanoribbons can result in rich electronic characteristics: they can be metallic or semiconducting depending on the chirality and edge configuration.



INTRODUCTION Graphene, a one-atom-thick planar sheet of sp2-hybridized carbon atoms packed in a honeycomb crystal lattice, has been considered as the “miracle material” of the 21st century since its experimental synthesis in 2004.1 This thinnest two-dimensional (2D) carbon material has extraordinary properties over other allotropes.2 Among others, the carriers in graphene can behave as massless Dirac Fermions3 and have a particular high mobility (106 m/s, close to the speed of light);4,5 graphene has been established as the strongest material ever measured, with a Young’s modulus of >1000 GPa;6 the thermal conductivity of graphene is in excess of carbon nanotubes (CNTs) or diamond.7 These properties endow graphene many exciting applications in numerous fields, such as supercapacitors,8 transparent conducting electrodes,9 and biosensors.10−12 In particular, graphene is expected as a substitutional material for silicon to use in microelectronics.13 However, pristine graphene is a semimetal with zero band gap, which is a principle obstacle limiting the utilization of graphene in electronics,14 such as p-n junctions. Thus an intermediate band gap (e.g., 1−3 eV) must be opened before any graphene-based microprocessor can be realized. Both experimental15 and theoretical studies16,17 have confirmed that cutting 2D graphene into one-dimensional (1D) graphene nanoribbons (GNRs) can generate a nonzero energy gap, irrespective of the width and edge chirality. In particular, the zigzag-edged GNRs have localized edge states18−20 and can be tuned into half-metals under external electric field21−23 or chemical modification.24−26 However, the band gaps in GNRs diminish gradually with increasing ribbon width,16 and so far the large-scale fabrication of GNRs is still a big challenge. © 2012 American Chemical Society

In recent years, many other approaches have been proposed toward engineering a band gap in graphene, such as perpendicular electric field,27−29 periodic antidot lattices,30−34 uniaxial strain,35−37 and covalent38−44 or noncovalent45−47 functionalization. However, the above-mentioned approaches either can only open a tiny band gap ( 5, the stain energy becomes negative while the band gap exceeds that of the C4H layer; all (n,n) tubes have a negative strain energy, while their band gaps are all larger than that of the C4H layer. The above results are not simple coincidences. As we demonstrated above, the electronic properties of the C4H layer and nanotubes are mainly determined by the isolated aromatic rings. For these C4H systems, the aromatic rings are more localized in the energetically more favorable ones and correspondingly present a more pronounced band gap. The stability of C4HNTs can also be hinted at by examining their formation energies. In Figure 6b we present the variation of formation energy of (n,0) and (n,n) C4HNTs as a function of tube diameters. As a comparison, the formation energy of the C4H layer was also computed. The positive formation energy of the C4H layer (0.11 eV/atom) indicates that the reaction between the hydrogen molecules and graphene is endothermic; thus, the C4H layer can only be obtained by exposing graphene to H plasmas (but not molecular H2). The formation energies of (n,0) and (n,n) C4HNTs both increase with increasing tube diameters, and all the tubes examined here have a lower (thus more favorable) formation energy than that of the C4H layer. It is not surprising since CNTs are more reactive than graphene, and these thin tubes have a higher curvature than thick ones, thus possess higher reactivity.91 In particular, (3,0) and (4,0) C4HNTs have negative formation energies, indicating that these two hydrogenated tubes can be realized by direct exposure of (6,0) and (8,0) CNTs to H2 molecules, respectively. Electronic Properties of C4H Nanoribbons. Finally, we investigated the electronic properties of C4HNRs. C4HNRs with zigzag or armchair edges can be obtained by cutting the carbon skeleton of the C4H layer along different directions. Following the convention of GNRs, we define the ribbon parameter N as the number of zigzag chains (Nz) across the ribbon width for zigzag ribbons (Figure 7a), and the number of dimer lines (Na) across the ribbon width for armchair ribbons (Figure 7b). For both zigzag and armchair C4HNRs, due to the ordering distribution of Csp3 atoms in the C4H layer, there exists two kinds of edge configurations: the edge sites occupied with pure Csp2 atoms (I) or alternate Csp3 and Csp2 atoms (II). Obviously, according to our definition, the Nz can only be even

band gaps. This is in stark contrast to the conventional (n,m) CNTs, which are metallic only if (n − m)/3 is an integer; otherwise, they are semiconducting.86,90 The uniform semiconducting property of the C4HNTs is attributed to the fact that all the Csp2 atoms are in the Clar sextet aromatic rings. For (n,0) C4HNTs, the band gap increases monotonically with the increase of n, indicating that the quantum confinement effect is not pronounced here. In particular, when n > 5, the bands gaps of (n,0) tubes are even larger than that of the C4H layer. In comparison, the band gaps of (n,n) tubes initially increase as n increases, and then from (5,5) tube converge to a constant value of 3.62 eV. Interestingly, even the smallest (n,n) tube has a larger band gap than that of the C4H layer. Overall, the uniformity of electronic properties would be a key advantage for the application of C4HNTs in nanotechnology. Stability of C4H Nanotubes. The stability of C4HNTs is a very important issue since it indicates the feasibility of experimental realization. To estimate the stability, we first examined the strain energies of both (n,0) and (n,n) C4HNTs as a function of n and compared them with their corresponding (2n,0) and (2n,2n) CNTs. Note that the strain energy only reflects the relative stability, and to some extent, the synthesis difficulty, but not the absolute stability. The strain energy of (n,0) C4HNTs decreases as the increase of n (Figure 6a), similar to single-walled carbon nanotubes. In particular, the negative strain energies for the nanotubes with n > 5 indicate that these thick (n,0) C4HNTs are energetically even more favorable than the C4H layer. In comparison, all the (n,n) C4HNTs we considered herein, including the thinnest (3,3) tube, have negative strain energies. The strain energy of (n,n) C4HNTs initially decreases with increasing tube diameter and then approaches a value of −0.02 eV at (8,8) and is expected to converge to zero (as that of the C4H layer) with the further increase of the tube diameter. The negative strain energies of the above studied C4HNT nanotubes (except (3,0) and (4,0)) strongly indicate that the C4H monolayer may easily bend to form nanotubes under appropriate conditions. The above results can be understood by considering the fact the Csp3 atoms in C4HNTs are closer to ideal sp3 hybridization than those of the C4H layer. Generally speaking, rolling the C4H layer into C4HNTs could release the strain of basal plane, and thus make C4HNTs energetically more favorable than the C4H layer. However, this reasoning is not valid for (3,0) and (4,0) tubes: their Csp3 atoms are closer to ideal sp3 hybridization than the thick tubes, but these two very thin nanotubes still have positive strain energies. Why? This is also not difficult to understand. Although the strain due to the deviation from ideal sp3 hybridization is much more released in very thin tubes, 4530

dx.doi.org/10.1021/jp212499h | J. Phys. Chem. C 2012, 116, 4526−4534

The Journal of Physical Chemistry C

Article

I-8-ZC4HNR is metallic with two energy levels crossing over from conduction band to valence band across the Fermi level in the band structure. The two energy levels are mainly contributed by the pz orbitals of edge Csp2 atoms which are not in the Clar sextet ring. In contrast, II-8-ZC4HNR is semiconducting with a 3.48 eV indirect band gap. The VBM and CBM are both localized at the Csp2 atoms of Clar sextet aromatic rings throughout the ribbon (see Figure S3 in Supporting Information). Both having edge Csp2 atoms, why is I-8-ZC4HNR metallic while II-8-ZC4HNR is semiconducting? Careful examination of the locations of the unhydrogenated Csp2 atoms gives us the answer. In II-8-ZC4HNR, all the Csp2 atoms are in Clar sextet rings, the same as these of C4H layer. In contrast, the edge Csp2 atoms in I-8-ZC4HNR do not belong to any Clar sextet ring, thus, it is the residual pz electrons of the edge Csp2 atoms that contribute to the metallicity of I-8-ZC4HNR. Both I-9-AC4HNR and II-9-AC4HNR are semiconducting, with a 3.74 and 2.78 eV direct band gap, respectively. For I-9AC4HNR, the VBM is mainly contributed by the two inner aromatic rings, while the CBM is localized at the Clar sextet aromatic rings throughout the ribbon (Figure S3, Supporting Information). Therefore, the edge states are not found in I-9AC4HNR, though in each edge there are two Csp2 atoms which do not belong to any Clar sextet ring. In contrast, in the band structure of II-9-AC4HNR, several new levels are at the top of the valence band and at the bottom of the conduction band, which shrink the band gap correspondingly. Both VBM and

Figure 7. The optimized geometric structures of (a) I-8-ZC4HNR, (b) II-8-ZC4HNR, (c) I-9-AC4HNR, and (d) II-9-AC4HNR.

while the Na can only be odd. Each edge atom is passivated with a hydrogen atom to ensure that the atom coordination of nanoribbons is the same as that of the C4H layer. Figure 7 presents the optimized structures of two representative zigzag and armchair C4HNRs, and their corresponding band structures are shown in Figure 8. Our computations demonstrated that these C4HNRs all have a nonmagnetic ground state.

Figure 8. Electronic band structures of (a) I-8-ZC4HNR, (b) II-8-ZC4HNR, (c) I-9-AC4HNR, and (d) II-9-AC4HNR. Panels (e) and (f) are the variations of the band gaps of zigzag (6 ≤ Nz ≤ 18) and armchair (5 ≤ Na ≤ 17) C4HNRs as a function of ribbon widths, respectively. 4531

dx.doi.org/10.1021/jp212499h | J. Phys. Chem. C 2012, 116, 4526−4534

The Journal of Physical Chemistry C



CBM are localized at the edge Csp2 atoms, indicating that the electronic properties of II-9-AC4HNR are determined by these edge states. We also investigated the electronic properties of zigzag and armchair C4HNRs with two types of edge configurations as a function of ribbon widths (Figure 8e,f). All the I-ZC4HNRs are metallic regardless of the ribbon width, while the II-ZC4HNRs are all semiconducting, and due to the quantum confinement effect, their band gaps decrease slightly with increasing ribbon width, and gradually approach that of the C4H layer. In comparison, both I- and II- ACH4NRs are semiconducting, the band gap of I-ACH4NRs decreases from 4.06 eV at Na = 5 to 3.62 eV at Na = 15, which shows the quantum confinement effect, while due to the presence of edge states, II-AC4HNRs all have a lower band gap than that of the C4H layer. With the increase of ribbon width, the edge states get less and less pronounced and the band gap could be enlarged correspondingly; at the same time the quantum confinement effect essentially lowers the band gap. As a result, the band gap of IIAC4HNRs converges to a constant value of 2.80 eV as ribbon width increases.



Article

ASSOCIATED CONTENT

S Supporting Information *

The band structure of the C4H layer computed using HSE06 functional, electronic profiles of VBM and CBM for C4H nanotubes and C4H nanoribbons, and the structural properties of a series of C4H nanotubes. This material is available free of charge via the Internet at http://pubs.acs.org.



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS Support by NSF Grant EPS-1010094 is gratefully acknowledged.



REFERENCES

(1) (a) Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Jiang, D.; Zhang, Y.; Dubonos, S. V.; Gregorieva, I. V.; Firsov, A. A. Electric field effect in atomically thin carbon films. Science 2004, 306, 666−669. (b) 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. (2) For some very recent reviews, see (a) Allen, M. J.; Tung, V. C.; Kaner, R. B. Honeycomb carbon: a review of graphene. Chem. Rev. 2010, 110, 132−145. (b) Castro Neto, A. H.; Novoselov, K. New directions in science and technology: two-dimensional crystals. Rep. Prog. Phys. 2011, 74, 082501. (c) Rao, C. N. R.; Subrahmanyam, K. S.; Ramakrishna Matte, H. S. S.; Govindaraj, A. Graphene: synthesis, functionalization and properties. Mod. Phys. Lett. B 2011, 25, 427−451. (d) Guo, S.; Dong, S. Graphene nanosheet: synthesis, molecular engineering, thin film, hybrids, and energy and analytical applications. Chem. Soc. Rev. 2011, 40, 2644−2672. (e) Terrones, M.; BotelloMendez, A. R.; Campos-Delgado, J.; Lopez-Urias, F.; Vega-Cantu, Y. I.; Rodriguez-Macias, F. J.; Elias, A. L.; Munoz-Sandoval, E.; CanoMarquez, A. G.; Charlier, J.-C.; Terrones, H. Graphene and graphite nanoribbons: morphology, properties, synthesis, defects and applications. Nano Today 2010, 5, 351−372. (3) Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Jiang, D.; Khotkevich, I. V.; Gregorieva, I. V.; Dubonos, S. V.; Firsov, A. A. Twodimensional gas of massless dirac fermions in graphene. Nature 2005, 438, 197−200. (4) Morozov, S. V.; Novoselov, K. S.; Katsnelson, M. I.; Schedin, F.; Elias, D.; Jaszczak, J. A.; Geim, A. K. Giant intrinsic carrier mobilities in graphene and its bilayer. Phys. Rev. Lett. 2008, 100, 016602. (5) Zhang, Y.; Tan, Y.-W.; Stormer, H. L.; Kim, P. Experimental observation of the quantum hall effect and Berry’s phase in graphene. Nature 2005, 438, 201−204. (6) Lee, C. G.; Wei, X. D.; Kysar, J. W.; Hone, J. Measurement of the elastic properties and intrinsic strength of monolayer graphene. Science 2008, 321, 385−388. (7) Balandin, A. A.; Ghosh, S.; Bao, W.; Calizo, I.; Teweldebrhan, D.; Miao, F.; Lau, C. N. Superior thermal conductivity of single-layer graphene. Nano Lett. 2008, 8, 902−907. (8) Stoller, M. D.; Park, S.; Zhu, Y.; An, J.; Ruoff, R. S. Graphenebased ultracapacitors. Nano Lett. 2008, 8, 3498−3502. (9) Sordan, R.; Traversi, F.; Russo, V. Logic gates with a single graphene transistor. Appl. Phys. Lett. 2009, 94, 073305. (10) Garaj, W.; Hubbard, W.; Reina, A.; Kong, J.; Branton, D.; Golovchenko, J. A. Graphene as a subnanometre trans-electrode membrane. Nature 2010, 467, 190−194. (11) Merchant, C. A.; Healy, K.; Wanunu, M.; Ray, V.; Peterman, N.; Bartel, J.; Fischbein, M. D.; Venta, K.; Luo, Z.; Johnson, A. T. C.;

CONCLUSION

In summary, we systematically investigated the structural and electronic properties of the recently achieved 2D C4H layer and its 1D derivatives: C4HNTs and C4HNRs by means of DFT computations. Structurally, the C4H layer is rather unique: all its Csp2 atoms are benzenoid while the Csp3 atoms slightly deviate from the ideal sp3 hybridization. Electronically, the C4H monolayer is semiconducting with a sizable band gap. However, this singlesided patterned partially hydrogenated C4H layer can only be obtained when we prevent the double-sided hydrogenation, since the double-sided graphane-embedded structure is energetically more favorable. Rolling up the C4H layer into C4HNTs can make the Csp3 atoms even closer to ideal sp3 hybridization. Consequently, the majority of C4HNTs (except (3,0) and (4,0)) have negative strain energies, and all the C4HNTs have more favorable formation energies than the C4H layer. These energetic evaluations indicate that C4HNTs can be feasibly realized. In contrast to carbon nanotubes whose band gap depends on the chirality and diameter, all C4HNTs have wide energy gaps regardless of the tube chirality and diameter. C4HNRs present much richer electronic properties than C4HNTs. ZC4HNRs can be either metallic or semiconducting: those terminated with Csp2 atoms not associated in Clar sextet rings are metallic due to the residual pz electrons of edge atoms, while those terminated with alternate Csp2 and Csp3 atoms are wide-gap semiconductors. AC4HNRs are always semiconducting: those terminated with pure Csp2 atoms all have a larger band gap than that of the C4H layer; in contrast, due to the presence of edge states, those terminated with alternate Csp2 and Csp3 atoms all have a lower band gap than that of the C4H layer. The above unique structural and electronic properties endow the C4H layer and its 1D analogues many promising applications in optics and opto-electronics devices, besides others. We hope our theoretical studies will help promote more experimental and theoretical studies of these novel materials and exploration of their applications in various fields. 4532

dx.doi.org/10.1021/jp212499h | J. Phys. Chem. C 2012, 116, 4526−4534

The Journal of Physical Chemistry C

Article

Drndić, M. DNA translocation through graphene nanopores. Nano Lett. 2010, 10, 2915−2921. (12) Schneider, G. F.; Kowalczyk, S. W.; Calado, V. E.; Pandraud, G.; Zandbergen, H. W.; Vandersypen, L. M. K.; Dekker, C. DNA translocation through graphene nanopores. Nano Lett. 2010, 10, 3163−3167. (13) Lemme, M. C.; Echtermeyer, T. J.; Baus, M.; Kurz, H. A graphene field-effect device. IEEE Electron Device Lett. 2007, 28, 282− 284. (14) Novoselov, K. Graphene: mind the gap. Nat. Mater. 2007, 6, 720−721. (15) Li, X. L.; Wang, X. R.; Zhang, L.; Lee, S. W.; Dai, H. J. Chemically derived, ultrasmooth graphene nanoribbon semiconductors. Science 2008, 319, 1229−1232. (16) Son, Y.-W.; Cohen, M. L.; Louie, S. G. Energy gaps in graphene nanoribbons. Phys. Rev. Lett. 2006, 97, 216803. (17) Barone, V.; Hod, O.; Scuseria, G. E. Electronic structure and stability of semiconducting graphene nanoribbons. Nano Lett. 2006, 6, 2748−2754. (18) Fujita, M.; Wakabayashi, K.; Nakada, K.; Kusakabe, K. Peculiar localized state at zigzag graphite edge. J. Phys. Soc. Jpn. 1996, 65, 1920−1923. (19) Nakada, K.; Fujita, M.; Dresselhaus, G.; Dresselhaus, M. S. Edge state in graphene ribbons: nanometer size effect and edge shape dependence. Phys. Rev. B 1996, 54, 17954−17961. (20) Wakabayashi, K.; Sigrist, M.; Fujita, M. Spin wave mode of edgelocalized magnetic states in nanographite zigzag ribbons. J. Phys. Soc. Jpn. 1998, 67, 2089−2093. (21) Son, Y.-W.; Cohen, M. L.; Louie, S. G. Half-metallic graphene nanoribbons. Nature 2006, 444, 347−349. (22) Hod, O.; Barone, V.; Peralta, J. E.; Scuseria, G. E. Enhanced half-metallicity in edge-oxidized zigzag graphene nanoribbons. Nano Lett. 2007, 7, 2295. (23) Kan, E. J.; Li, Z. Y.; Yang, J. L.; Hou, J. G. Will zigzag graphene nanoribbon turn to half metal under electric field? Appl. Phys. Lett. 2007, 91, 243116. (24) Kan, E. J.; Li, Z. Y.; Yang, J. L.; Hou, J. G. Half-metallicity in edge-modified zigzag graphene nanoribbons. J. Am. Chem. Soc. 2008, 130, 4224−4225. (25) Li, Y. F.; Zhou, Z.; Shen, P. W.; Chen, Z. Spin gapless semiconductor-metal-half-metal properties in nitrogen-doped zigzag graphene nanoribbons. ACS Nano 2009, 3, 1952−1958. (26) 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. (27) Zhang, Y. B.; Tang, T. T.; Girit, C.; Zhao, H.; Martin, M. C.; Zettl, A.; Crommie, M. F.; Shen, Y. R.; Wang, F. Direct observation of a widely tunable bandgap in bilayer graphene. Nature 2009, 459, 820− 823. (28) Castro, E. V.; Novoselov, K. S.; Morozov, S. V.; Peres, N. M. R.; Dos Santos, J.; Nilsson, J.; Guinea, F.; Geim, A. K.; Neto, A. H. C. Biased bilayer graphene: semiconductor with a gap tunable by the electric field effect. Phys. Rev. Lett. 2007, 99, 216802. (29) Samarakoon, D. K.; Wang, X.-Q. Tunable band gap in hydrogenated bilayer graphene. ACS Nano 2010, 4, 4126−4130. (30) Ouyang, F. P.; Peng, S. L.; Liu, Z. R.; Liu, Z. F. Bandgap opening in graphene antidot lattices: the missing half. ACS Nano 2011, 5, 4023−4030. (31) Petersen, R.; Pedersen, T. G.; Jauho, A.-P. Clar sextet analysis of triangular rectangular and honeycomb graphene antidot lattice. ACS Nano 2011, 5, 523−529. (32) Li, Y. F.; Zhou, Z.; Shen, P. W.; Chen, Z. Two-dimensional polyphenylene: experimentally available porous graphene as a hydrogen purification membrane. Chem. Commun. 2010, 46, 3672− 3674. (33) Du, A. J.; Smith, S. C. Multifunctional porous graphene for nanoelectronics and hydrogen storage: new properties revealed by first principle calculations. J. Am. Chem. Soc. 2010, 132, 2876−2877.

(34) Zhang, A. H.; Teoh, H. F.; Dai, Z. X.; Feng, Y. P.; Zhang, C. Band gap engineering in graphene and hexagonal BN antidot lattices: a first principles study. Appl. Phys. Lett. 2011, 98, 023105. (35) Gui, G.; Li, J.; Zhong, J. X. Band structure engineering of graphene by strain: first-principles calculations. Phys. Rev. B 2008, 78, 075435. (36) Gyunea, F.; Katsnelson, M. I.; Geim, A. K. Energy gaps and a zero-field quantum Hall effect in graphene by strain engineering. Nat. Phys. 2010, 6, 30−33. (37) Sun, L.; Li, Q. X.; Ren, H.; Su, H. B.; Shi, Q. W.; Yang, J. L. Strain effect on electronic structures of graphene nanoribbons: a firstprinciples study. J. Chem. Phys. 2008, 129, 074704. (38) Liu, L.-H.; Yan, M.-D. Simple method for the covalent immobilization of graphene. Nano Lett. 2009, 9, 3375−3378. (39) He, H.; Gao, C. General approach to individually dispersed, highly soluble, and conductive graphene nanosheets functionalized by nitrene chemistry. Chem. Mater. 2010, 22, 5054−5064. (40) Quintana, M.; Spyrou, K.; MarekGrzelczak, M.; Browne, W. R.; Rudolf, P.; Prato, M. Functionalization of graphene via 1,3-dipolar cycloaddition. ACS Nano 2010, 4, 3527−3533. (41) Suggs, K.; Reuven, D.; Wang, X.-Q. Electronic properties of cycloaddition functionalized graphene. J. Phys. Chem. C 2011, 115, 3313−3315. (42) Zhang, H.; Bekyaova, E.; Huang, J.-W.; Zhao, Z.; Bao, W. Z.; Wang, F. L.; Haddon, R. C.; Lau, C. N. Aryl functionalization as a route to band gap engineering in single layer graphene devices. Nano Lett. 2011, 11, 4047−4051. (43) Niyogi, S.; Bekyarova, E.; Hong, J.; Khizroev, S.; Berger, C.; de Heer, W.; Haddon, R. C. Chemical modification of epotaxial graphene: spontaneous grafting of aryl groups. J. Am. Chem. Soc. 2009, 131, 1336−1337. (44) Niyogi, S.; Bekyariva, E.; Hong, J.; Khizroev, S.; Berger, C.; de Heer, W.; Haddon, R. C. Covalent chemistry for graphene electronics. J. Phys. Chem. Lett. 2011, 2, 2487−2498. (45) Lu, Y. H.; Chen, W.; Feng, Y. P.; He, P. M. Tuning the electronic structure of graphene by an organic molecule. J. Phys. Chem. B 2009, 113, 2−5. (46) Nduwimana, A.; Wang, X.-Q. Energy gaps in supramolecular functionalized graphene nanoribbons. ACS Nano 2009, 3, 1995−1999. (47) Zhang, Z. X.; Huang, H. L.; Yang, X. M.; Zang, L. Tailoring electronic properties of graphene by π−π stacking with aromatic molecules. J. Phys. Chem. Lett. 2011, 2, 2897−2905. (48) Sofo, J. O.; Chaudhari, A. S.; Barber, G. D. Graphane: a twodimensional hydrocarbon. Phys. Rev. B 2007, 75, 153401. (49) Elias, D. C.; Nair, R. R.; Mohiuddin, T. M. G.; Morozov, S. V.; Blake, P.; Halsall, M. P.; Ferrari, A. C.; Boukhvalov, D. W.; Katsnelson, M. I.; Geim, A. K.; Novoselov, K. S. Control of graphene’s properties by reversible hydrogenation: evidence for graphane. Science 2009, 323, 610−613. (50) Ryu, S. M.; Han, M. Y.; Maultzsch, J.; Heinz, T. F.; Kim, P.; Steigerwald, M. L.; Brus, L. E. Reversible basal plane hydrogenation of graphene. Nano Lett. 2008, 8, 4597−4602. (51) Balog, R.; Jorgensen, B.; Wells, J.; Laegsgaard, E.; Hofmann, P.; Besenbacher, F.; Hornekaer, L. Atomic hydrogen adsorbate structures on graphene. J. Am. Chem. Soc. 2009, 131, 8744−8745. (52) Sessi, P.; Guest, J. R.; Bode, M.; Guisinger, N. P. Patterning graphene at the nanometer scale via hydrogen desorption. Nano Lett. 2009, 9, 4343−4347. (53) Balog, R.; Jørgensen, B.; Nilsson, L.; Andersen, M.; Rienks, E.; Bianchi, M.; Fanetti, M.; Lægsgaard, E.; Baraldi, A.; Lizzit, S.; Sljivancanin, Z.; Besenbacher, F.; Hammer, B.; Pedersen, T. G.; Hofmann, P.; Hornekær, L. Bandgap opening in graphene induced by patterned hydrogen adsorption. Nat. Mater. 2010, 9, 315−319. (54) Grassi, R.; Low, T.; Lundstorm, M. Scaling of the energy gap in pattern-hydrogenated graphene. Nano Lett. 2011, 11, 4574−4578. (55) Boukhvalov, D. W.; Katsnelson, M. I. Tuning the gap in bilayer graphene using chemical functionalization: density functional calculations. Phys. Rev. B 2008, 78, 085413. 4533

dx.doi.org/10.1021/jp212499h | J. Phys. Chem. C 2012, 116, 4526−4534

The Journal of Physical Chemistry C

Article

(56) Samarakoon, D. K.; Wang, X.-Q. Chair and twist-boat membranes in hydrogenated graphene. ACS Nano 2010, 4, 4017− 4022. (57) Li, Y. F.; Zhou, Z.; Shen, P. W.; Chen, Z. Structural and electronic properties of graphane nanoribbons. J. Phys. Chem. C 2009, 113, 15043−15045. (58) Lu, N.; Li, Z. Y.; Yang, J. L. Electronic structure engineering via on-plane chemical functionalization: a comparison study on twodimensional polysilane and graphane. J. Phys. Chem. C 2009, 113, 16741−16746. (59) Lebegue, S.; Klintenberg, M.; Eriksson, O.; Katsnelson, M. I. Accurate electronic band gap of pure and functionalized graphane from GW calculations. Phys. Rev. B 2009, 79, 245117. (60) Zhou, J.; Wang, Q.; Sun, Q.; Chen, X. S.; Kawazoe, Y.; Jena, P. Ferromagnetism in semihydrogenated graphene sheet. Nano Lett. 2009, 9, 3867−3870. (61) Singh, A. K.; Penev, E. S.; Yakobson, B. I. Vacancy clusters in graphane as quantum dots. ACS Nano 2010, 4, 3510−3514. (62) Yang, M.; Nurbawono, A.; Zhang, C.; Feng, Y. P. Ariando. Twodimensional graphene superlattice made with partial hydrogenation. Appl. Phys. Lett. 2010, 96, 193115. (63) Gao, H. L.; Wang, L.; Zhao, J. J.; Ding, F.; Lu, J. P. Band gap tuning of hydrogenated graphene: H coverge and configuration dependence. J. Phys. Chem. C 2011, 115, 3236−3242. (64) Nair, R.; Ren, W.; Jalil, R.; Riaz, I.; Kravets, V.; Britnell, L.; Blake, P.; Schedin, F.; Mayorov, A. S.; Yuan, S.; et al. Fluorographene: a two-dimensional counterpart of Teflon. Small 2010, 6, 2877−2884. (65) Robinson, J. T.; Burgess, J. S.; Junkermeier, C. E.; Badescu, S. C.; Reinecke, T. L.; Perkins, F. K.; Zalautdniov, M. K.; Baldwin, J. W.; Culbertson, J. C.; Sheehan, P. E.; et al. Properties of fluorinated graphene films. Nano Lett. 2010, 10, 3001−3005. (66) Withers, F.; Bointon, T. H.; Dubois, M.; Russo, S.; Craciun, M. F. Nanopatterning of fluorinated graphene by electron beam irradiation. Nano Lett. 2011, 11, 3912−3916. (67) Lee, W. -K.; Robinson, J. T.; Gunlycke, D.; Stine, R. R.; Tamanaha, C. R.; King, W. P.; Sheehan, P. E. Chemically isolated graphene nanoribbons reversibly formed in fluorographene using polymer nanowire masks. Nano Lett. 2011, 11, 5461−5464. (68) Samarakoon, D. K.; Chen, Z. F.; Nicolas, C.; Wang, X. Q. Structural and electronic properties of fluorographene. Small 2011, 7, 965−969. (69) Singh., R.; Bester, G. Hydrofluorinated graphene: twodimensional analog of polyvinylidene fluoride. Phys. Rev. B 2011, 84, 155427. (70) Sahin, H.; Topsakal, M.; Ciraci, S. Structure of fluorinated graphene and their signatures. Phys. Rev. B 2011, 83, 115432. (71) Wu, X. S.; Sprinkle, M.; Li, X. B.; Ming, F.; Berger, C.; de Heer, W. A. epitaxial−graphene/graphene-oxide junction: an essential step towards epitaxial graphene electronics. Phys. Rev. Lett. 2008, 101, 026801. (72) Jung, I.; Dikin, D. A.; Piner, R. D.; Ruoff, R. S. Tunable electrical conductivity of individual graphene oxide sheets reduced at “low” temperatures. Nano Lett. 2008, 8, 4283−4287. (73) Jia-An, Y.; Lede, X.; Chou, M. Y. Structural and electronic properties of oxidized graphene. Phys. Rev. Lett. 2009, 103, 086802. (74) Xiang, H. J.; Wei, S. H.; Gong, X. G. Structural motifs in oxidized graphene: A genetic algorithm study based on density functional theory. Phys. Rev. B 2010, 82, 035416. (75) Samarakoon, D. K.; Wang, X.-Q. Twist-boat conformation in graphene oxides. Nanoscale 2011, 3, 192−195. (76) Haberer, D.; Guusca, C. E.; Wang, Y.; Sachdev, H.; Fedorov, A. V.; Farjam, M.; Jafari, S. A.; Vyalikh, D. V.; Usachov, D.; Liu, X. J.; et al. Evidence for a new two-dimensional C4H-type polymer based on hydrogenated graphene. Adv. Mater. 2011, 23, 4497−4503. (77) Kresse, G.; Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B 1993, 47, 558−561. (78) Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953−17979.

(79) Kresse, G.; Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 1999, 59, 1758−1775. (80) Perdew, J. P.; Burke, L.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (81) (a) Heyd, J.; Scuseria, G. E.; Ernzerhof, M. Hybrid functionals based on a screened Coulomb potential. J. Chem. Phys. 2003, 118, 8207−8215. (b) Heyd, J.; Scuseria, G. E.; Ernzerhof, M. Erratum: hybrid functionals based on a screened Coulomb potential. J. Chem. Phys. 2006, 124, 219906. (82) For a review, see Barone, V.; Hod, O.; Peralta, J. E.; Scuseria, G. E. Accurate prediction of the electronic properties of low-dimensional graphene derivatives using a screened hybrid density functional. Acc. Chem. Res. 2011, 44, 269−279. (83) (a) Clar, E. Polycyclic Hydrocarbons; Academic Press: New York, 1964. (b) Clar, E. The Aromatic Sextet; Wiley: London, 1972. (c) Randic, M. Aromaticity of polycyclic conjugated hydrocarbons. Chem. Rev. 2003, 103, 3449−3606. (84) (a) Wassmann, T.; Seitsonen, A. P.; Saitta, A. M.; Lazzeri, M.; Mauri, F. Structure, stability, edge states, and aromaticity of graphene ribbons. Phys. Rev. Lett. 2008, 101, 096402. (b) Wassmann, T.; Seitsonen, A. P.; Saitta, A. M.; Lazzeri, M.; Mauri, F. Clar’s theory, πelectron distribution, and geometry of graphene nanoribbons. J. Am. Chem. Soc. 2010, 132, 3440−3451. (c) Gao, X. F.; Zhao, Y. L.; Liu, B.; Xiang, H, J.; Zhang, S. B. π-Bond maximization of graphene in hydrogen addition reactions. Nanoscale 2011, DOI: 10.1039/ clnr11048a. (85) Moran, D.; Stahl, F.; Bettinger, H. F.; Schaefer, H. F. III; Schleyer, P. v. R. Toward graphite: magnetic properties of large polybenzenoid hydrocarbons. J. Am. Chem. Soc. 2003, 125, 6746− 6752. (86) (a) Ormsby, J. L.; King, B. T. Clar valence bond representation of π-bonding in carbon nanotubes. J. Org. Chem. 2004, 69, 4287−4291. (b) Matsuo, Y.; Tahara, K.; Nakamura, E. Theoretical studies on structures and aromaticity of finite-length armchair carbon nanotubes. Org. Lett. 2003, 5, 3181−3184. (c) Bettinger, H. F. Effects of finite carbon nanotube length on sidewall addition of fluorine atom and methylene. Org. Lett. 2004, 6, 731−734. (d) Martín-Martínez, F. J.; Melchor, S.; Dobado, J. A. Clar−Kekulé structuring in armchair carbon nanotubes. Org. Lett. 2008, 10, 1991−1994. (e) Baldoni, M.; Selli, D.; Sgamellotti, A.; Mercuri, F. Unraveling the reactivity of semiconducting chiral carbon nanotubes through finite-length models based on Clar Sextet theory. J. Phys. Chem. C 2009, 113, 862−866. (f) Balaban, A. T.; Klein, D. J. Claromatic carbon nanostructures. J. Phys. Chem. C 2009, 113, 19123−19133. (87) Manaa, M . R.; Sprehn, D. W.; Ichord, H. A. Prediction of extended aromaticity for a novel C48N12 azafullerene structure. J. Am. Chem. Soc. 2002, 123, 13990−13991. (88) Chen, Z .F.; Jiao, H. J.; Moran, D.; Hirsch, A.; Thiel, W.; Schleyer, P. v. R. Aromatic stabilization in heterofullerenes C48X12 (X = N, P, B, Si). J. Phys. Org. Chem. 2003, 16, 726−730. (89) Haddon, R. C. Hybridization and the orientation and alignment of .Pi. orbitals in nonplanar conjugated organic molecules: .Pi. orbitals axis vector analysis (POAV2). J. Am. Chem. Soc. 1986, 108, 2837− 2842. (90) Lu, X.; Chen, Z. F. Curved p-conjugation, aromaticity and the related chemistry of small fullerenes (