Relative Stability of Graphene Nanoflakes Under Environmentally

Jul 11, 2013 - Ian K. Snook,. ‡ and Amanda S. Barnard*. ,†. †. CSIRO Materials Science & Engineering, Parkville, Australia. ‡. Applied Physics...
0 downloads 0 Views 2MB Size
Article pubs.acs.org/JPCC

Relative Stability of Graphene Nanoflakes Under Environmentally Relevant Conditions Hongqing Shi,†,‡ Lin Lai,† Ian K. Snook,‡ and Amanda S. Barnard*,† †

CSIRO Materials Science & Engineering, Parkville, Australia Applied Physics, RMIT University, Melbourne, Australia



S Supporting Information *

ABSTRACT: The stability of finite-sized graphene nanostructures under different conditions is contingent on the stability of the highly reactive undercoordinated atoms around the circumference. Chemical passivation is a convenient way of stabilizing the edges and corners, provided the conditions support adsorption. In this paper we examine the stabilities of hydrogen, oxygen, hydroxyl, and water functionalization of the edges of graphene nanoflakes using a combination of ab initio thermodynamics and self-consistent charge density functional tight-binding simulations. We find that the adsorption of hydrogen and oxygen on graphene nanoflakes is selective and sensitive to the edge and corner structure and the local environment. Under hydrogen-rich conditions (or in vacuo), we find that armchair edges are preferred, whereas under oxygen-rich conditions the stability of the zigzag edge is enhanced. In addition to this, we find that the types of corners play an important role, and the presence of acute 60° corners is thermodynamically preferred, particularly under humid conditions. In each case, the adsorption efficiency is related to the temperature and the pressure, which may be used to stabilize one type of structure or another.



INTRODUCTION Since the successful isolation in 2004,1 graphene has been the focus of intense research, and cost-effective methods to manufacture bulk quantities for commercial applications are already being developed. Graphene has extremely high carrier mobility and capacity, high electron transfer rate, excellent ability to quench fluorescence, high thermal conductivity, maximal surface-to-volume ratio, and outstanding robustness and flexibility. On the basis of these desirable properties, and despite its short history, graphene has already demonstrated great potentials in various novel sensor applications, such as physical sensors for detection of photons, magnetic fields, mass, and strain. However, graphene has one severe limitation from the point of view of electronics applications: it has no band gap and a vanishingly small density of states at the Fermi level, making it a semimetal.2,3 Decreasing the dimensionality is one of the methods to modify graphene4 and can open a band gap and tune it to broaden the range of applications, and in many cases this involves restricting the in-plane dimension and/or introducing undercoordinated atoms. The introduction of reactive undercoordinated atoms around the circumference of finite-sized graphene nanostructures is also critical to the development of sensing applications. Graphene has already been used as electrochemical sensors, electronic sensors, optical sensors, and nanopore sensors for biological or chemical detection, all of which draw upon graphene’s extraordinary electrical, optical, and chemical properties.5 In the latter case, knowledge of the interactions of graphene with the local chemical environment is imperative, and controlling these interactions is the topic of much research. Published 2013 by the American Chemical Society

For example, oxidation of carbon is one of the most common chemical reactions when hydrophobic graphene is exposed to ambient conditions, but the controlled oxidation of carbon to add functional O-groups is one way of tuning the chemical activity.6 When taken to the extreme, a new material (graphene oxide) is formed, which is hydrophilic.7−9 Graphene oxide formed by oxidizing crystalline graphite is similar to a graphene sheet with the in-plane carbon atoms being directly functionalized with oxygen-containing groups (such as carboxyl, hydroxyl, ketone, and ether).10−12 Graphene oxide is an insulator but has applications in transparent conductive films, paper-like and composite materials, energyrelated materials, and biological and medical applications. Reduction of graphite oxide is also a practical way to synthesize graphene sheets directly.13 Experimental data have shown that some hydroxyl functional groups exist also at reduced graphene sheets,14 but the configuration at higher coverage, where epoxy and hydroxyl groups coexist at the graphene surface,15 is not fully known. However, circumferential O-functionalized graphene is not the same as graphene oxide, as in the former case it is only the edges and corners that are decorated with Ogroups, whereas the latter also involves planar adsorption. Oxygen is also present to some degree in the vast majority of graphene synthesis methods.16 In contrast, the hydrogenated graphene has great potential in other kinds of applications, due to different characteristic Received: April 28, 2013 Revised: June 12, 2013 Published: July 11, 2013 15375

dx.doi.org/10.1021/jp404200r | J. Phys. Chem. C 2013, 117, 15375−15382

The Journal of Physical Chemistry C



properties. For example, patterning H atoms on graphene in a one-dimensional (1-D) manner leads to so-called graphene nanoroads, which have properties similar to graphene nanoribbons.17 Two-dimensional hydrogenation was found both theoretically and experimentally to be an effective way to convert graphene from a gapless semimetal to a semiconductor with a tunable band gap,18−20 and partial hydrogenation may induce interesting magnetic states that have potential applications in spintronics.21 Hydrogenation is usually a deliberate process, during which the absorption and diffusion of H atoms on graphene are of great importance to determine the hydrogenation pattern and related electronic properties. Hydrogen is also the most common terminal group used in computational modeling of graphene.22 Chemical reactions are even more complicated when we consider the zero-dimensional (0-D) form, known as graphene nanoflakes (GNFs) or graphene quantum dots (GQDs).23 The bottom-up synthesis of graphene using chemical methods often produces GNFs and GQDs, and there are already ways of producing GNFs from the top down. Substrate-free synthesis has been reported,24−26 which makes GNFs (and GQDs) a viable candidate for integration into future devices. The advantage of using these materials is that they have many more structural degrees of freedom than graphene membranes or graphene nanoribbons. They possess edge states and corner states and can range in size from molecular to semi-infinite 2-D structures.27 They can also be engineered with different geometries, each with unique (and potentially tunable) properties. Since GNFs may be cut into a much larger variety of different shapes,28−31 we have previously established a categorical scheme based on the geometric complexity.23 Class 1 of 0-D graphene is characterized by a single type of edges (arm-chair or zigzag), where the edges and corners are of the same type. Class 2 of 0-D graphene is also characterized by a single type of edge but with the opposite types of corners. Class 3 of 0-D graphene nanostructures is characterized by a combination of arm-chair- and zigzag-edges, and any may have both arm-chair- and zigzag-corners. In this paper, we have investigated the stability of the circumferential oxygen, hydrogen, and hydroxyl termination of Class 1, Class 2, and Class 3 GNFs with respect to varied temperatures and different O- and/or H-rich environments using a combination of ab initio thermodynamics and computer simulations based on the density functional tight-binding theory. Our results show that the thermostabilities of the O-, H-, and OH-passivation of GNFs are highly dependent upon the surrounding temperature and/or the chemical environment, as well as the geometric category of the shape. The O- and Hpassivation is more stable than the OH-passivation, unless the supersaturations of O and H are simultaneously high (humid conditions). We also find that the thermodynamically preferred edge type depends on the chemical environment and that since the adsorption efficiency is related to the temperature and the pressure these may be used along with chemical functionalization to support one type of structure or another. There has been some seemingly contradictory evidence in the literature as to which type of edges (armchair or zigzag) are most likely,32,33 and our results suggest that both types of edges may be stable, depending on the local conditions. These results are compared with recent experimental findings where available.

Article

COMPUTATIONAL METHODOLOGY

As mentioned above, this study uses the density functional based tight-binding method with self-consistent charges (SCCDFTB).34,35 In this approach, the Kohn−Sham density functional is expanded to second order around a reference electron density, which is obtained from self-consistent density functional calculations of weakly confined neutral atoms within the generalized gradient approximation (GGA). The confinement potential is optimized to anticipate the charge density and effective potential in molecules and solids and unrelated to charge transfer. A short-range element-specific diatomic repulsive potential accounts for double counting terms in the Coulomb and exchange-correlation contributions as well as the internuclear repulsion. A minimal valence basis is established, and one- and two-center tight-binding matrix elements are explicitly calculated within DFT. The PBC set of parameters developed by Köhler and Frauenheim36 is used to describe the contributions from diatomic interactions of hydrogen, oxygen, and carbon. The Fermi level is smeared by 300 K, which leads to the electronic temperature of 1 meV. Self-consistency is converged within 0.01 meV at the level of Mulliken charges.35 All structures were fully relaxed first with the final forces on the atoms being less than 0.5 meV/Å. DFTB is a highly transferable method that overcomes the main limitation of density functional theory (DFT) at the reduced computational cost of standard tight binding; a detailed comparison of DFTB with DFT, for graphene systems, is provided in ref 37. This approach has also proven successful in studying the electronic properties of GNFs in the past.23,28−31,38−40 In this study all calculations include spin polarization and full noncollinear spin relaxation, to ensure that each structure is relaxed with respect to the magnetic moment. Both ferromagnetic and ferrimagnetic initial spin states were tested to ensure that the spin relaxation converged, and the final magnetic configuration was independent of the starting configuration. A comparison with simple, self-consistent, ionic, and electronic relaxations (that included only spin polarization) confirmed that the thermodynamic results are largely unchanged with the addition of magnetic effects, as the energy difference falls within the uncertainties of the method. Nevertheless, they have been included for completeness and will be reported elsewhere as part of a large and detailed comparison of the size , shape, and chemical dependence of the magnetic properties of GNFs. To evaluate the stability of O-, H-, and OH-functionalized GNFs, we have used a standard thermodynamic formalism where the formation energy per atom, Ef, is defined as Ef = (μC,X − nCμC − nX μ X )/Nt

(1)

where μC,X denotes the chemical potential of each GNF combined with adsorbate (X), which is equivalent to its static total energy. The nC represents the number of C atoms; nX represents the number of adsorbed X atoms, molecules, or groups; and Nt represents the total number of atoms in the passivated GNF. The chemical potential of the carbon atom, μC, is equal to the total energy per C atom in the fully relaxed unterminated (radical) version of each structure. This means that we are calculating the relative formation energy of each system, to evaluate the stability of a particular combination, rather than a total formation energy (where the structures are 15376

dx.doi.org/10.1021/jp404200r | J. Phys. Chem. C 2013, 117, 15375−15382

The Journal of Physical Chemistry C

Article total where Etotal O2 and EH2 are the total energy of a free, isolated O2 and H2 molecule at T = 0 K, and hvO2/2 and hvH2/2 are the vibrational energy of an O2 and H2 molecule. This thermodynamic formalism is well established and routinely used throughout ab initio thermodynamics to study the environmental stability of materials and surfaces.43−54 As mentioned above, a detailed comparison of the computational method (with density functional theory) is provided, for graphene, in ref 37. On the basis of this combination of methodologies, it is estimated that the uncertainties in the following simulations are on the order of 0.05 eV.

products of synthetic reactions). It is worth pointing out that this definition of the formation energy has been published in work elsewhere and that the absolute value of the atomic energy cancels out.41,42 One can also easily include a comparison to an infinite graphene membrane by simply subtracting the constant bulk-graphene free energy (or any other frame-of-reference containing carbon, such as graphene, graphite, diamond, amorphous carbon, or a carbide), and this will not alter the relative stability of the different nanoflakes. It should be pointed out, however, that to obtain the correct absolute stability one should select the correct bulk frame-ofreference from the phase diagram of carbon, which will be graphite or diamond (depending on the pressure and temperature) and not graphene at all. The formation energy Ef is linked to the environmental conditions (such as temperature, pressure, and chemistry) via the chemical potential of adsorbates μX. The chemical potential μX can be derived in a standard way from the combination of chemical potentials of gaseous molecules at chemical equilibrium. Since we assume a gaseous chemical reservoir, the chemical potential of any kind of molecule can be defined as μX = E X +

⎛ pV ⎞ hνX + kBT ln⎜ X ⎟ 2 ⎝ kBT ⎠



RESULTS AND DISCUSSION The fully relaxed structures of unsaturated GNFs are presented in Figure 1, with representative structures from Classes 1, 2,

(2)

Here, EX represents the static total energy of the molecule obtained from calculations at the level of theory, and with the same computational settings, p, T, kB, h, and νX are the external pressure, temperature, Boltzmann’s constant, Planck’s constant, and sum of vibrational frequencies of the molecule (using experimental data) in the chemical reservoir, respectively. Finally, VX indicates the quantum volume of the molecule ⎛ ⎞3/2 h2 VX = ⎜ ⎟ ⎝ 2πmX kBT ⎠

(3)

where mX is the mass of adsorbate X. It is clearly shown in eq 2 that the chemical potential can be altered by changing the environmental temperature and/or pressure and by choosing suitable chemical reservoirs. At a given temperature, the pressure is the partial pressure (or supersaturation) of H2 (or O2); conversely, at a given supersaturation, the temperature is the surrounding temperature, which contributes to the vibrational entropy. The oxygen and hydrogen chemical potentials (μO and μH) are half of the chemical potential of O2 and H2 molecules, respectively. Most oxygen- and hydrogen-rich conditions can be defined as the point beyond which gaseous O and H would start to condense on the sample. However, in the temperature and pressure range we are interested in, a condensed O2- and H2-solid phase does not exist (the critical temperature of O2 and H2, i.e., above which the gas and liquid phase are in thermodynamic coexistence, is Tc ≈ 150 K and 30 K). Thus, an appropriate and well-defined estimate of the upper limit of the oxygen and hydrogen chemical potentials is !

max[μO(T , p)] = !

max[μH (T , p)] =

hνO2 ⎞ 1 ⎛ total ⎜EO2 + ⎟ 2⎝ 2 ⎠

(4)

hνH2 ⎞ 1 ⎛ total ⎜E H 2 + ⎟ 2⎝ 2 ⎠

(5)

Figure 1. Fully relaxed geometry of Class 1 trigonal GNFs: (a) with AC-edges and 468 C atoms and (b) with ZZ-edges and 481 C atoms; Class 2 rhombic GNFs: (c) with AC-edges and 448 C atoms; and (d) with ZZ-edges and 448 C atoms; Class 2 hexagonal GNFs: (e) with AC-edges and 480 C atoms, and (f) with ZZ-edges and 486 C atoms; and a Class 3 rectangular GNF: (g) with both AC- and ZZ-edges and 448 C atoms.

and 3.23 Class 1 is defined as having edges and corners of the same chirality (either arm-chair or zigzag); Class 2 is defined as having edges and corners of opposing chirality; and Class 3 is defined as having edges (and corners) with both arm-chair and zigzag chirality. We have selected model structures with between four and five hundred atoms, which provides a reasonable comparison while reducing complexity associated with size dependence.23,28,30,31 From a purely geometric perspective, four different shapes are included. The trigonal (Class 1), hexagonal (Class 2), and rhombic (Class 2) GNFs have a unique edge of arm-chair (AC) or zigzag (ZZ), while the rectangular flake (Class 3) has both 15377

dx.doi.org/10.1021/jp404200r | J. Phys. Chem. C 2013, 117, 15375−15382

The Journal of Physical Chemistry C

Article

combination of high temperature and/or P ≈ 0 will help to purify the GNFs, by encouraging H desorption. The latter is representative of the majority of electronic structure simulations of graphene reported in the literature, and using the same computational method as described above we have shown that the lowest energy ″radical″ structure is the ACedged Class 2 hexagonal nanoflake.28,30 As we increase the supersaturation of H, in the range μH = −2.91 to −3.08 eV, the ZZ-edged hexagonal (Class 2) GNF gains stability. As we increase the supersaturation of H a little further (μH = −2.86 to −2.91 eV), the ZZ-edged rhombic (Class 2) GNF becomes the most stable configuration. Under H-rich conditions, when μH = −0.11 to −2.86 eV, the ZZedged trigonal (Class 1) GNF with the hydrogen adsorption is the most stable. Across this entire range the ZZ-edge is preferred, but as we increase the hydrogen partial pressure, we see an increased preference for shapes with a higher edge-toplanar ratio (moving from Class 2 to Class 1). When H is fully saturated, or oversaturated (μH > −0.11 eV), the AC-edged trigonal (Class 1) GNF is the most stable configuration. At this point we see a change from ZZ- to ACedges, which again explains the common computational conclusion that AC-edges are most stable, as computational studies traditionally use H as a generic passivant. We can also see some trends associated with the temperature and/or pressure. When the ambient temperature is lower than room temperature (RT), Class 1 AC-edged GNFs are the most stable, with acute 60° corners. When the environmental temperature is close to, or higher than RT, Class 1 ZZ-edged GNFs dominant, which also have acute 60° corners. When the temperature is much higher, Class 1 ZZ-edged GNFs gradually give way to hexagonal Class 2 ZZ-edged GNFs, with obtuse 120° corners, via the rhombic GNF (has both acute and obtuse angles). Oxygen (O) Adsorption. The adsorption of oxygen atom on GNFs has also been studied using ab initio thermodynamics and the structures shown in Figure 1. The relative formation energies of these structures with respect to the oxygen gas reservoir are presented in Figure 3.

kinds of edges. As mentioned above, oxygen, hydrogen, and hydroxyl adsorption were considered with a complete circumferential monolayer (termed: 1 ML) with all the undercoordinated carbon atoms functionalized. In total, we have investigated 21 functionalized GNFs, and all the relaxed configurations with functionalization can be found in the Supporting Information. We also attempted to passivate the circumference with H2O and confirmed that water will not passivate the circumference of graphene in the range of temperatures under investigation, which is consistent with the fact that graphene is hydrophobic. During these simulations, in most cases, the H2O desorbed or deprotonated, neither of which resulted in stable H2O passivation. The final (unstable) structures are included in the Supporting Information. This property is useful in the application of graphene for a possible nanofiltration system (with water playing the role of a carrier).60 In contrast, an aqueous dispersion of chemically converted graphene (CCG) can be readily prepared by controlled reduction of graphene oxide in water.60 Compared to ″pristine″ graphene, CCG is more susceptible to corrugation due to the presence of some sp3 hybridized carbon atoms and topological defects. The amplitude of corrugation of CCG sheets in water can be readily controlled by hydrothermal treatment, leading to a new class of permeation-tunable nanofiltration membranes. Hydrogen (H) Adsorption. Using eq 1 the adsorption of the hydrogen atom on GNFs has been calculated for each of the 14 structure/chemical/coverage combinations (see Supporting Information). The formation energies of these structures with respect to the hydrogen gas reservoir are presented in Figure 2. We can see from Figure 2 when the chemical potential of H is low (H-poor conditions) the unterminated ZZ-edged Class 2 GNF is the most stable configuration. This indicates that a

Figure 2. Relative formation energies of hydrogen adsorption on Class 1, 2, and 3 GNFs. The dashed vertical lines indicate the upper limitation of the hydrogen chemical potential, μH(T,p), using half of Etotal H2 + (hvH2/2) as zero reference. The corresponding temperatures are given for two selected pressures: one corresponding to ultrahigh vacuum (UHV) conditions and the other to atmospheric pressure. The results are labeled in order of increasing stability at 1 ML hydrogen coverage: (A) AC-edged hexagonal (Class 2), (B) ZZ-edged hexagonal (Class 2), (C) AC-edged rhombic (Class 2), (D) ZZ-edged rhombic (Class 2), (E) AC/ZZ-edged rectangle (Class 3), (F) ACedged trigonal (Class 1), and (G) ZZ-edged trigonal (Class 1). The unterminated graphene reference system is represented by a solid horizonal line, which is set to zero, and the solid vertical lines highlight changes in the relative stability of different structures.

Figure 3. Relative formation energies of oxygen adsorption on GNFs with different shapes and edges. The dashed vertical lines indicate the upper limitation of the oxygen chemical potential, μO(T,p), using half of Etotal O2 + (hvO2/2) as zero reference. The corresponding temperatures are given for two selected pressures: one corresponding to ultrahigh vacuum (UHV) conditions and the other to atmospheric pressure. The results are labeled using the same method as in Figure 2. 15378

dx.doi.org/10.1021/jp404200r | J. Phys. Chem. C 2013, 117, 15375−15382

The Journal of Physical Chemistry C

Article

We can see from Figure 3 that at low oxygen partial pressure (under anaerobic conditions) oxygen dissociates from the circumference of GNFs, and clean ″radical structure″ is preferred. The corresponding transition temperature for this reaction is very high (even in vacuo), which is consistent with the known stability of graphene oxide at room temperature and supports the application of O-terminated graphene under the same conditions. For higher oxygen supersaturation, up to O-rich aerobic conditions (μO > −3.81 eV), we find that Class 1 ZZ-edged GNFs are the most stable configuration. Once again, we find that acute 60° corners (characteristic of Class 1 structures) are the most preferred, which suggests that acute corners influence the chemical activity of GNFs under saturated conditions. Recently graphene single crystals were grown using lowpressure chemical vapor deposition in copper-foil enclosures using methane as a precursor,61 and scanning electron microscopy (SEM) revealed the graphene domains tended to have high ″edge roughness″, with numerous acute corners. Hydroxyl (OH) Adsorption. The relative formation energy for GNFs terminated with hydroxyl was also investigated and is provided as a function of hydroxyl chemical potential as given in Figure 4, with respect to a gaseous O2 + H2 reservoir. In

Figure 5. Thermochemical mapping for the GNFs with hydrogen, oxygen, hydroxyl, and water adsorption in ″constrained equilibrium″ with a H2 and O2 environment on the left. The formation energy of water vapor is indicated by a black solid line.

most stable structure is an AC-edged radical (unterminated) graphene nanoflake of Class 2. Under anaerobic conditions (low oxygen and high hydrogen chemical potential) the most stable structures are the Class 1 AC-edged trigonal GNF with complete edge hydrogenation (1 ML coverage). Under aerobic conditions (high oxygen and low hydrogen chemical potential) the most stable structures are the Class 1 ZZ-edged trigonal GNF with complete edge oxidation (1 ML coverage). Finally, under humid conditions, with high hydrogen and high oxygen chemical potentials, the preferred GNFs are Class 1 AC-edged structures, with complete circumferential passivation. However, under the majority of environmentally relevant conditions, we expect that ZZ-edges and acute corners will prevail. Although not discussed in detail here, we also found that (in addition to the role that the environmental conditions play in the thermodynamic stability of different configurations) the complexity of the gaseous environment can affect the mechanical stability of the (most stable) GNFs. Both of the OH-functionalized configurations (including the energetically preferred Class 1 AC/OH combination) developed out-ofplane distortions during relaxation (as shown in the Supporting Information), which was necessary to accommodate and stabilize the steric effects between neighboring OH groups. Recently experiments have efficiently produced hydroxylfunctionalized graphene from exfoliation of graphite powder by ball milling the presence of potassium hydroxide (KOH).62 The participation of hydroxyl groups was confirmed by Fourier transform infrared (FTIR), nuclear magnetic resonance analysis (NMR), and Raman and X-ray photoelectron spectroscopic (XPS) measurements, and graphene sheets were imaged by atomic force microscopy (AFM), scanning electron microscopy (SEM), transmission electron microscopy (TEM), and UV−vis spectroscopy. These structures were produced at room temperature and show good hydrophilicity, excellent electroactivity, and biocompatibility with human RPE cells, which is consistent with an environment rich in OH or H2O and supported by the present results. Our results are also consistent with a recent study exploring the influence of reactive edges

Figure 4. Relative formation energies of hydroxyl adsorption on GNFs with different shapes and edges. The dashed vertical lines indicate the upper limitation of the hydroxyl chemical potential, μOH(T,p), using Etotal OH + (hvOH2/2) as zero reference. The corresponding temperatures are given for two selected pressures: one corresponding to ultrahigh vacuum (UHV) conditions and the other to atmospheric pressure. The results are labeled using the same method as in Figure 2.

general, we can see that hydroxyl functionalization of the edges and corners is stable under ambient conditions and up to high temperature, but under dry conditions the unterminated ACedged Class 2 nanoflake is thermodynamically preferred. For the μOH in the range of −0.82 and −5.62 eV, the Class 1 ZZ-edged GNF is preferred, and once again we see a change to the Class 1 AC-edged GNF when μOH > −0.82 eV. This can be summarized by saying that ZZ-edges are favored under dry conditions, and AC-edges are favored under humid conditions. Comparison and Map of Chemical Stability. By combining all of these results (above), a map of the chemical stability of different classes of GNFs with hydrogen, oxygen, and hydroxyl adsorption is presented in Figure 5. This mapping technique was also successfully used in ref 54. This map highlights that under UHV (and near UHV) conditions, with both low oxygen and low hydrogen chemical potentials, the 15379

dx.doi.org/10.1021/jp404200r | J. Phys. Chem. C 2013, 117, 15375−15382

The Journal of Physical Chemistry C

Article

temperature, pressure, or both. From a practical standpoint, this also means that in a realistic environmental chamber, or when developing graphene nanoflakes for sensor applications, the edge structure could potentially be engineered to provide greater selectivity when screening specific gases within gaseous mixtures, or more complex ″multifunctional″ graphene geometries could be developed that are sensitive to more than one gas at a time.

and defects of substrates for achieving a complete cycle of room-temperature molecular sensing using graphene.55 An important aspect of these results is how sensitively the chemical environment influences edge stability, which has been a topic of some discussion in recent years. It was established that the ZZ-edges are kinetically stable by considering the effect of ejecting an atom at the edge for each chirality under electron irradiation at this energy,32 but the influence of thermodynamics remained a point of contention. ZZ-edges have some intrinsic advantages as there are antibonding π-electron states (edge states) localized along edges in the ZZ direction but not in the AC direction.33,56,57 Since the antibonding edge states encourage chemical adsorption, a graphene gas sensor with ZZedges may be more sensitive to gaseous environments. In contrast, the AC-edges are well-known to be thermodynamically stable,28,30,57 due to the energetically favorable reconstructions involving sp1 bonds, and do not possess magnetic edge states.58 This can be rationalized with Clar’s Rule59 and is significant because a preference for one type of edge over another will have important implications for the electronic and magnetic properties. At this point it should be highlighted that, like any nanoscale material, kinetics will also play an important role in determining the structure of graphene nanoflakes under different conditions. The present work is restricted to thermodynamic considerations, and so we cannot speculate as to the kinetics; however, we can see that reversible functionalization will require a change in the thermodynamics environment to be thermodynamically driven, or to lower a kinetic energy barrier as required. In situations where the crossing points in the freeenergy diagrams are shallow, we expect kinetics to be important, as (within uncertainties) we approach a state of mutual thermodynamic coexistence between the competing structures. Although not discussed here, defects (such as vacancies with different shapes and sizes) can also impact the stability of GNFs under these conditions.40 The chemical activity of unsaturated carbon atoms along the boundary of a defect will be comparable to the carbon atoms around the circumference and can therefore be mapped in a similar way. There are far more configurations to be examined when we combine defects with the current geometric classes, but since this issue will have a profound impact on the application of graphene for gas sensing,55 this will be a topic of future work.



ASSOCIATED CONTENT

S Supporting Information *

Schematic representations of all of the fully relaxed, H, O, OH, and H2O terminated graphene nanoflake structures included in this study. 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 This research has been supported by the Australian Research Council under Discovery Grant DP110101362. Computational resources for this project have been supplied by the National Computing Infrastructure (NCI) national facility under MAS Grant e74. The authors would like to thank Hongxin Yang for useful discussions. This paper is dedicated to the memory of Ian Snook.



REFERENCES

(1) Novoselov, K. S.; Geim, A.; Morozov, S.; Jiang, D.; Zhang, Y.; Dubonos, S.; Grigorieva, I.; Firsov, A. Electric Field Effect in Atomically Thin Carbon Films. Science 2004, 306, 666−669. (2) Castro Neto, A.; Guinea, F.; Neres, N. M. R.; Novoselov, K.; Geim, A. The Electronic Properties of Graphene. Rev. Mod. Phys. 2009, 81, 109−162. (3) Abergel, D.; Apalkov, V.; Berashevich, J.; Zieler, K.; Chakraborty, T. Properties of Graphene: A Theoretical Perspective. Adv. Phys. 2010, 59, 261−482. (4) Park, S.; Ruoff, R. S. Chemical Methods for the Production of Graphenes. Nat. Nanotechnol. 2009, 4, 217−224. (5) Liu, Y.; Dong, X.; Chen, P. Biological and Chemical Sensors Based on Graphene Materials. Chem. Soc. Rev. 2012, 41, 2283−2307. (6) Mestl, G.; Maksimova, N. I.; Keller, N.; Roddatis, V. V.; Schlogl, R. Carbon Nanofilaments in Heterogeneous Catalysis: An Industrial Application for New Carbon Materials. Angew. Chem., Int. Ed. 2001, 40, 2066−2068. (7) Li, D.; Muller, M. B.; Gilje, S.; Kaner, R. B.; Wallace, G. G. Processable Aqueous Dispersions of Graphene Nanosheets. Nat. Nanotechnol. 2008, 3, 101−105. (8) Li, D.; Kaner, R. B. Graphene-based Materials. Science 2008, 320, 1170−1171. (9) Mao, S.; Pu, H.; Chen, J. Graphene Oxide and its Reduction: Modeling and Experimental Progress. RSC Adv. 2012, 2, 2643−2662. (10) Akbar Bagri, A.; Mattevi, C.; Acik, M.; Chabal, Y. J.; Chhowalla, M.; Shenoy, V. B. Structural Evolution During the Reduction of Chemically Derived Graphene Oxide. Nat. Chem. 2010, 2, 581−587. (11) Bagri, A.; Grantab, R.; Medhekar, N. V.; Shenoy, V. B. Stability and Formation Mechanisms of Carbonyl- and Hydroxyl-Decorated Holes in Graphene Oxide. J. Phys. Chem. C 2010, 114, 12053−12061. (12) Hernández Rosas, J. J.; Ramírez Gutiérrez, R. E.; EscobedoMorales, A.; Anota, E. C. First Principles Calculations of the Electronic



CONCLUSION In this paper we have used ab initio thermodynamics and selfconsistent charge density functional based tight-binding simulations to investigate the relative stability of GNF functionalization by hydrogen, oxygen, and hydroxyl. We find that the geometry (Class) and angle subtended at the corners play an important role in determining the stability, particularly when the partial pressure of H or O is low (in vacuo). In contrast, when the partial pressures are high, the relative stability depends on the types of edges, or (alternatively) the stability of different edges is moderated by the environment. Zigzag edges are preferred under a wider range of conditions, but a H-rich environment is necessary to support AC-edge stability. Since the surrounding environment affects the stability of H-, O-, and OH-adsorption on GNFs so dramatically, it may be possible to use this knowledge to probe the structure of samples of graphene nanostructure by simply adjusting the 15380

dx.doi.org/10.1021/jp404200r | J. Phys. Chem. C 2013, 117, 15375−15382

The Journal of Physical Chemistry C

Article

and Chemical Properties of Graphene, Graphane, and Graphene Oxide. J. Mol. Model. 2011, 17, 1133−1139. (13) Stankovich, S.; Dikin, D. A.; Dommett, C. H. B.; Kohlhaas, K. M.; Zimney, E. J.; Stach, E. A.; Piner, R. D.; Nguyen, S. T.; Ruoff, R. S. Graphene-based Composite Materials. Nature 2006, 442, 282−286. (14) Gilje, S.; Dubin, S.; Badakhshan, A.; Farrar, J.; Danczyk, S. A.; Kaner, R. B. Photothermal Deoxygenation of Graphene Oxide for Patterning and Distributed Ignition Applications. Adv. Mater. 2010, 22, 419−423. (15) Boukhvalov, D. W.; Katsnelson, M. I. J. Modeling of Graphite Oxide. J. Am. Chem. Soc. 2008, 130, 10697−10701. (16) Choi, W. B.; Lahiri, I.; Seelaboyina, R.; Kang, Y. S. Synthesis of Graphene and Its Applications: A Review. Crit. Rev. Solid State Mater. Sci. 2010, 35, 52−71. (17) 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.; et al. Control of Graphene’s Properties by Reversible Hydrogenation: Evidence for Graphane. Science 2009, 323, 610−613. (18) 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. (19) Haberer, D.; Vyalikh, D. V.; Taioli, S.; Dora, B.; Farjam, M.; Fink, J.; Marchenko, D.; Pichler, T.; Ziegler, K.; Simonucci, S.; et al. Tunable Band Gap in Hydrogenated Quasi-Free-Standing Graphene. Nano Lett. 2010, 10, 3360−3366. (20) Yang, M.; Nurbawono, A.; Zhang, C.; Feng, Y. P.; Ariando. Two-Dimensional Graphene Superlattice Made with Partial Hydrogenation. Appl. Phys. Lett. 2010, 96, 193115. (21) Zhou, J.; Wang, Q.; Chen, X. S.; Kawazoe, Y.; Jena, P. Ferromagnetism in Semihydrogenated Graphene Sheet. Nano Lett. 2009, 9, 3867−3870. (22) Snook, I. K.; Barnard, A. S. Graphene Nano-Flakes and NanoDots: Theory, Experiment and Applications. In Physics and Applications of Graphene - Theory; Mikhailov, S., Ed.; InTech: Croatia, 2011. (23) Shi, H. Q.; Barnard, A. S.; Snook, I. K. Quantum Mechanical Properties of Graphene Nano-flakes and Quantum Dots. Nanoscale 2012, 4, 6761−6767. (24) Dato, A.; Radmilovic, V.; Lee, Z.; Phillips, J.; Frenklach, M. Substrate-Free Gas-Phase Synthesis of Graphene Sheets. Nano Lett. 2008, 8, 2012−2016. (25) Dato, A.; Frenklach, M. Substrate-Free Microwave Synthesis of Graphene: Experimental Conditions and Hydrocarbon Precursors. New J. Phys. 2010, 12, 125013. (26) Zenasni, A.; Delamoreanu, A.; Rabot, C. Free-Suspended Graphene Synthesis via Carbon Diffusion Through Platinum-Based Metal. Appl. Phys. Lett. 2012, 100, 151907. (27) Gass, M.; Bangert, U.; Bleloch, A.; Wang, P.; Nair, R.; Geim, A. Free-Standing Graphene at Atomic Resolution. Nat. Nanotechnol. 2008, 3, 676−681. (28) Barnard, A. S.; Snook, I. K. Modelling the Role of Size, Edge Structure and Terminations on the Electronic Properties of Graphene Nano-flakes. Modell. Simul. Mater. Sci. Eng. 2011, 19, 054001. (29) Barnard, A. S.; Snook, I. K. Ripple Induced Changes in the Wavefunction of Graphene: An Example of a Fundamental Bymmetry Breaking. Nanoscale 2012, 4, 1167−1170. (30) Shi, H. Q.; Barnard, A. S.; Snook, I. K. Modelling the Role of Size, Edge Structure, and Terminations on the Electronic Properties of Trigonal Graphene Nano-flakes. Nanotechnology 2012, 23, 065707. (31) Shi, H. Q.; Barnard, A. S.; Snook, I. K. High Throughput Theory and Simulation of Nanomaterials: Exploring the Stability and Electronic Properties of Nanographene. J. Mater. Chem. 2012, 22, 18119−18123. (32) Ö . Girit, Ç .; Meyer, J. C.; Erni, R.; Rossell, M. D.; Kisielowski, C.; Yang, L.; Park, C.-H.; Crommie, M. F.; Cohen, M. L.; Louie, S. G.; Zettl, A. Graphene at the Edge: Stability and Dynamics. Science 2009, 323, 1705−1708.

(33) Enoki, T.; Kobayashi, Y.; Fukui, K.-I. Electronic Structures of Graphene Edges and Nanographene. Int. Rev. Phys. Chem. 2007, 26, 609−645. (34) Porezag, D.; Frauenheim, Th.; Köhler, Th.; Seifert, G.; Kaschner, R. Construction of Tight-Binding-Like Potentials on the Basis of Density-Functional Theory: Application to Carbon. Phys. Rev. B 1995, 51, 12947. (35) Frauenheim, Th.; Seifert, G.; Elstner, M.; Niehaus, Th.; Köhler, C.; Amkreutz, M.; Sternberg, M.; Hajnal, Z.; Di Carlo, A.; Suhai, S. Atomistic Simulations of Complex Materials: Ground-State and Excited-State Properties. J. Phys.: Condens. Matter 2002, 14, 3015. (36) Köhler, C.; Frauenheim, T. Molecular Dynamics Simulations of CFx (x= 2,3) Molecules at Si3N4 and SiO2 Surfaces. Surf. Sci. 2006, 600, 453−460. (37) Zobelli, A.; Ivanovskaya, V.; Wagner, P.; Suarez-Martinez, I.; Yaya, A.; Ewels, C. P. A Comparative Study of Density Functional and Density Functional Tight Binding Calculations of Defects in Graphene. Physica Status Solidi B 2012, 249, 276−282. (38) Barnard, A. S.; Snook, I. K. Transformation of Graphene into Graphane in the Absence of Hydrogen. Carbon 2010, 48, 981−986. (39) Barnard, A. S.; Snook, I. K. Size- and Shape-Dependence of the Graphene to Graphane Transformation in the Absence of Hydrogen. J. Mater. Chem. 2010, 20, 10459−10464. (40) Shi, H. Q.; Barnard, A. S.; Snook, I. K. Site-Dependent Stability and Electronic Structure of Single Vacancy Point Defects in Hexagonal Graphene Nano-Flakes. Phys. Chem. Chem. Phys. 2013, 15, 4897− 4905. (41) Lai, L.; Barnard, A. S. Modeling the Thermostability of Surface Functionalisation by Oxygen, Hydroxyl, and Water on Nanodiamonds. Nanoscale 2011, 3, 2566−2575. (42) Lai, L.; Barnard, A. S. Surface Phase Diagram and Thermodynamic Stability of Functionalisation of Nanodiamonds. J. Mater. Chem. 2012, 22, 16774−16780. (43) Reuter, K.; Stampfl, C.; Scheffler, M. Ab Initio Thermodynamics and Statistical Mechanics of Surface Properties and Functions. In Handbook of Materials Modeling, Part A. Methods; Yip, S., Ed.; Springer: Berlin, 2005. (44) Wang, X.-G.; Weiss, W.; Shaikhutdinov, Sh. K.; Ritter, M.; Petersen, M.; Wagner, F.; Schlögl, R.; Scheffler, M. The Hematite (Alpha-Fe2O3)(0001) surface: Evidence for Domains of Distinct Chemistry. Phys. Rev. Lett. 1998, 81, 1038−1041. (45) Wang, X.-G.; Chaka, A.; Scheffler, M. Effect of the Environment on Alpha-Al2O3(0001) Surface Structures. Phys. Rev. Lett. 2000, 84, 3650−3653. (46) Reuter, K.; Scheffler, M. First-Principles Atomistic Thermodynamics for Oxidation Catalysis: Surface Phase Diagrams and Catalytically Interesting Regions. Phys. Rev. Lett. 2003, 90, 046103. (47) Reuter, K.; Scheffler, M. Composition and Structure of the RuO2(110) Surface in an O2 and CO Environment: Implications for the Catalytic Formation of CO2. Phys. Rev. B 2003, 68, 045407. (48) Raty, J. P.; Galli, G. Ultradispersity of Diamond at the Nanoscale. Nat. Mater. 2003, 2, 792−795. (49) Lundgren, E.; Gustafson, J.; Mikkelsen, A.; Andersen, J. N.; Stierle, A.; Dosch, H.; Todorova, M.; Rogal, J.; Reuter, K.; Scheffler, M. Kinetic Hindrance During the Initial Oxidation of Pd(100) at Ambient Pressures. Phys. Rev. Lett. 2004, 92, 046101. (50) Barnard, A. S.; Zapol, P.; Curtiss, L. A. Anatase and Rutile Surfaces with Adsorbates Representative of Acidic and Basic Conditions. Surf. Sci. 2005, 582, 173−188. (51) Rohrer, J.; Hyldgaard, P. Ab-initio Thermodynamics of Deposition Growth: Surface Terminations of CVD Titanium Carbide and Nitride. Phys. Rev. B 2010, 82, 045415. (52) Paisley, E. A.; Losego, M. D.; Gaddy, B. E.; Tweedie, J. S.; Collazo, R.; Sitar, Z.; Irving, D. L.; Maria, J. P. Surfactant-Enabled Epitaxy Through Control of Growth Mode with Chemical Boundary Conditions. Nat. Commun. 2011, 2, 461. (53) Li, H.; Zhang, W.; Smith, J. R. Advances in Ab initio Thermodynamic Studies on Metal/Oxide Interfaces. Phys. Status Solidi A 2011, 208, 1166−1173. 15381

dx.doi.org/10.1021/jp404200r | J. Phys. Chem. C 2013, 117, 15375−15382

The Journal of Physical Chemistry C

Article

(54) Ijäs, M.; Havu, P.; Harju, A. Interaction of Chlorine with StoneWales Defects in Graphene and Carbon Nanotubes and Thermodynamical Prospects of Chlorine-Induced Nanotube Unzipping. Phys. Rev. B 2013, 87, 205430. (55) Randeniya, L. K.; Shi, H. Q.; Barnard, A. S.; Fang, J.; Martin, P. J.; Ostrikov, K. Harnessing the Influence of Reactive Edges and Defects of Substrates for Achieving Complete Cycle of RoomTemperature Molecular Sensing using Graphene. Small 2013, DOI: 10.1002/smll.201300689. (56) Enoki, T.; Takai, K. Unconventional Electronic and Magnetic Functions of Nanographene-Based Host−Guest Systems. Dalton Trans. 2008, 29, 3773−3781. (57) Enoki, T. Role of Edges in the Electronic and Magnetic Structures of Nanographene. Phys. Scr. 2012, 146, 014008. (58) Joly, V. L. J.; Kiguchi, M.; Hao, S.-J.; Takai, K.; Enoki, T.; Sumii, R.; Amemiya, K.; Muramatsu, H.; Hayashi, T.; Kim, Y. A.; et al. Observation of Magnetic Edge State in Graphene Nanoribbons. Phys. Rev. B 2010, 81, 245428. (59) Wassmann, T.; Seitsonen, A. P.; Marco Saitta, A.; Lazzeri, M.; Mauri, F. Structure, Stability, Edge States, and Aromaticity of Graphene Ribbons. Phys. Rev. Lett. 2008, 101, 096402. (60) Qiu, L.; Zhang, X.; Yang, W.; Wang, Y.; Simon, G. P.; Li, D. Controllable Corrugation of Chemically Converted Graphene Sheets in Water and Potential Application for Nanofiltration. Chem. Commun. 2011, 47, 5810−5812. (61) Li, X.; Magnuson, C. W.; Venugopal, A.; Tromp, R. M.; Hannon, J. B.; Vogel, E. M.; Colombo, L.; Ruoff, R. S. Large-Area Graphene Single Crystals Grown by Low-Pressure Chemical Vapor Deposition of Methane on Copper. J. Am. Chem. Soc. 2011, 133, 2816−2819. (62) Yan, L.; Lin, M.; Zeng, C.; Chen, Z.; Zhang, S.; Zhao, X.; Wu, A.; Wang, Y.; Dai, L.; Qu, J.; et al. Electroactive and Biocompatible Hydroxyl-Functionalized Graphene by Ball Milling. J. Mater. Chem. 2012, 22, 8367−8371.

15382

dx.doi.org/10.1021/jp404200r | J. Phys. Chem. C 2013, 117, 15375−15382