Formation Mechanisms of Graphitic-N: Oxygen Reduction and

Dae Sik Kim , Dong Jae Chung , Hyeong-Il Park , Mohammad Zahid Ansari , Taeseup Song , Hansu ... Moumita Kotal , Hyunjun Kim , Sandipan Roy , Il-Kwon ...
0 downloads 0 Views 2MB Size
Subscriber access provided by FLORIDA ATLANTIC UNIV

Article

Formation Mechanisms of Graphitic-N: Oxygen Reduction and Nitrogen Doping of Graphene Oxides Wei-Wei Wang, Jing-Shuang Dang, Xiang Zhao, and Shigeru Nagase J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.5b10607 • Publication Date (Web): 24 Feb 2016 Downloaded from http://pubs.acs.org on February 25, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Formation Mechanisms of Graphitic-N: Oxygen Reduction and Nitrogen Doping of Graphene Oxides Wei-Wei Wang, †,‡Jing-Shuang Dang, †Xiang Zhao,* †and Shigeru Nagase*‡



Institute for Chemical Physics & Department of Chemistry, School of Science, Xi’an Jiaotong

University, Xi’an 710049, China. ‡

Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan

* Fax: +86 29 82668559; Tel: +86 29 82665671; Email: [email protected]; [email protected] .

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 32

Abstract

Deoxygenation and nitrogen doping mechanisms of graphene oxides with participation of foreign NH3 molecules were investigated by density functional theory calculations. Firstly, reduction on perfect graphene oxide without any structural defect is proved to be facilitated at high temperature but the following doping process is impracticable because of the huge energy requirement for C-C cleavage. To elucidate the formation of hexagonal graphitic-N, oxygen reduction and subsequent nitrogen doping processes on defective graphene oxides with single vacancies were explored for the first time. All possible reaction pathways were taken into account and the results demonstrate that formation of graphitic-N from NH3 and defective graphene oxides with one carbonyl or two hydroxyl groups is feasible in energy. The dominant reaction route is found to be exothermic with a practical reaction rate of 2.26 × 106 s-1 at 900 ℃, which is in a good agreement with experimental observations.

Keywords Density functional theory, graphene oxides, nitrogen doping, oxygen reduction, structural defects

ACS Paragon Plus Environment

2

Page 3 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Introduction Ever since its isolation in 2004, the two-dimensional graphene has attracted tremendous attentions owing to its exotic properties.1 To date, great achievements have been made on graphene-based nanomaterials in the fields of chemistry, physics, and life science.2-19 Besides the well-recognized pristine graphene structure, both experimental and theoretical efforts were made to design and synthesize heteroatom-doped graphene materials to tailor the physicochemical properties for extending its applications. For example, owing to the unique electronic properties, nitrogen-doped graphene (NG, in which one or more carbon atoms are substituted by nitrogen atoms) is expected to be utilized in fuel cells and other electrochemical devices.20-26 In general, there are two approaches to obtain nitrogen-doped graphene in experiments. Firstly, NG can be yielded by directly mixing the carbon source with nitrogen-containing compounds (such as ammonia and hydrazine) as precursors.27,28 Beside such direct synthesis, another way is post-synthesis treatment. By thermal or plasma treatment on existing graphene or graphene oxide (GO), nitrogen atoms can be inserted into the network of the carbon skeleton (edges, defective sites, and also the interior of the carbon sheet).23-25,29-36Structurally, there are three types of bonding configurations for a doped atomic nitrogen in graphene, designated as graphitic-N, pyridinic-N, and pyrrolic-N (see Figure 1). Pyridinic-N and pyrrolic-N are formed at the edge or defective sites of graphene. Graphitic-N, which bonds with three carbon atoms, is formed inside the network of graphene. Figure 1 In 2009, Li et al reported a simple chemical method to obtain NG in bulk quantities. 32 By annealing GO in NH3, the N-doped graphene sheets were synthesized at elevated temperatures.

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 32

XPS spectra revealed that nitrogen atoms were embedded into carbon skeleton and the graphiticN exhibited a significant role over 900°C.32 Although such an effective method has been proposed in experiments, the oxygen reduction and following nitrogen doping processes of GO are still less known. Therefore, discussions on reaction pathways and corresponding influential factors in theory are necessary to understand the mechanisms in detail. Moreover, to our bestknowledge, almost all the theoretical studies on the reactivity of GO only focus on the oxygen reductions and there is no report on the nitrogen doping mechanisms so far. In this study, comprehensive density functional theory (DFT) computations were performed to explore the reduction and doping mechanisms of GO by using ammonia as the reducing agent and meanwhile the nitrogen source for doping. Our objectives are to elucidate the interactions between ammonia with pristine and defective graphene oxides with various oxygen-containing functional groups, and to uncover the formation mechanisms of graphitic-N from foreign NH3 and parental GO sheet. Computational details Experimental evidence has shown that the oxygen functionalities distribute on both the edge and basal plane of the defective GO sheet, and those oxygen groups lead to a separation of GO sheet into small in-plane aromatic domains.37-41 Accordingly, we modeled the structures of GO as finite fragmental carbon sheets with various oxygen groups (Figure 2), which have been extensively used for theoretical calculations of graphene oxides.39, 42-45The molecular radius of graphene flake is around 6.17 Å.As for the oxygen-containing functionalities, epoxide (-O-), hydroxyl (-OH), and carbonyl (-C=O) groups were established as the major components because their existence have been well identified.46-52 In addition, the minor component carboxyl (COOH) group was excluded in the present work since carboxyl is only distributed at the edge of

ACS Paragon Plus Environment

4

Page 5 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

GO whereas our present study are constrained on the formation of graphitic-N in basal plane of GO. Figure 2 All calculations were performed with the GAUSSIAN 09 program.53 The M06-2X density functional in conjunction with the basis set of 6-31G(d) were employed for all structural optimizations and vibrational frequency calculations. The correction of zero-point energy (ZPE) to the total energy of each molecule was taken into account. The electronic ground states for all GO models in Figure 2 were confirmed as closed-shell singlet with the singlet-triplet splitting energies more than 28.0 kcal mol-1. In the present study, the structures of transition states were located by using the Berny algorithm. Vibrational analyses were conducted to clarify the nature of stationary points as global minima or transition states with one imaginary frequency. On the basis of transition state theory (TST), the reaction rates were evaluated by the Arrhenius formula: νexp(−Eb/kBT),54 where ν is the attempt frequency, T is the reaction temperature, Eb is the computed energy barrier, and kB is the Boltzmann constant. Results and discussion We initially focused on the reaction of NH3to the perfect graphene sheet with singleepoxide group (GO-1 in Figure 2). Epoxide was chosen as the major oxygen-containing component because this type of functional group exists in abundance at both the edge and interior of an aromatic domain of GO. A previous report by Tang et al indicated that attack of NH3 on the inner epoxide will cause the reduction of graphene oxides, and the nitrogen-containing groups (NH2 and -NH) will externally attach onto the carbon surface.55 Therefore, the adsorption of NH3 groups can be considered as the prelude of nitrogen doping. Kinetically, the reaction of NH3 to

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 32

GO with single epoxide group was confirmed as a stepwise process. In the first instance, the foreign NH3 molecule is adsorbed on GO with a binding energy of -4.6 kcal mol-1. After that, a hydrogen atom in NH3 transfers from nitrogen to oxygen, resulting in a new C-NH2 bond and a – OH group on carbon surface. In the third step, the second N-H dissociation takes place and the H atom attacks on oxygen to yield a H2O molecule. According to the different binding sites of – NH2 in the second step, the reduction is divided into two reaction routes, as shown in Figure 3 and 4. Route-a in which -NH2 is located at the para position of oxygen is considered as a more favorable pathway because of a lower energy barrier for hydrogen transfer to form the H2O molecule. Based on our DFT calculations, the second N-H dissociation (TS3-a)actsas the ratedetermining step with an energy barrier of 42.0 kcal mol-1 (route-a), suggesting that the reduction of epoxide on graphene oxide is facilitated at elevated temperature. After elimination of H2O, the external nitrogen atom is attached onto the basal plane of the graphene to form a NH-containing reduced product (G-NH-1), as shown in Figure 3. Figure 3 Figure 4 Furthermore, the subsequent nitrogen doping process on pristine graphene (G-NH-1 in Figure 3 and 5) was explored here. In general, nitrogen doping includes two major steps: N-H dissociation and C-N formation. According to different sequence of the two steps and distinct positions of newly-formed C-N bonds, four reaction pathways were found to elucidate the doping process. As shown in Figure 5 and 6, the reaction is preferred to generate a sp 3 carbon in the first step (INT1-cd in route-c(d)) and then to form new C-N bonds (C4-N bond of NG-CH-1 in route-c and C3-N bond of NG-CH-2 in route-d, respectively). However, energetically, even

ACS Paragon Plus Environment

6

Page 7 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

assisted by the migrated hydrogen atom (route-d, in Figure 5 and 6), this substitution process is still endothermic with a huge activation barrier of 107.2 kcal mol-1, which indicates that the doping is unlike to take place even at elevated temperatures. Figure 5 Figure 6 Accordingly, direct doping from flawless graphene network seems to be unpractical. It should be mentioned that such a conclusion can also apply to the case of GO with hydroxyl groups located on the interior of the aromatic domain. It has been demonstrated that the reaction of NH3 to hydroxyl-participated GO will generate a water molecule as the reduction product and a -NH2 group which is adsorbed on carbon surface.55 However, similar to the case of epoxide, since dissociation of C-C bonds on perfect graphene skeleton needs huge energy requirements, the nitrogen substitution cannot take place in practice. Based on this reason, the reduction and doping reactions of GO with single hydroxyl group are not considered in this work. Different from the hypothetical GO model with structurally perfect graphene layer, experimental observations suggested that carbon surface is always interspersed with vacancies and the most reported defects are single vacancies (SV) and double vacancies (DV) which are thermodynamically stable.56-62 In the case of SV, the missing carbon atom from hexagonal lattice leads to a defective segment which includes a pentagon-nonagon pair and a sp-hybridized carbon with dangling bonds. Herein, the reactions of NH3 to defective GO with SV (GO-2, GO-3, and GO-4 in Figure 2) were explored, to determine whether it is feasible to dope nitrogen into the basal plane of carbon sheet to generate the experimentally observed graphitic-N. The carbonyls and hydroxyls were selected as the oxygen-containing groups because of the identified large

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 32

amounts of sp3 carbon with C-O and C=O in experiments.32 As for the double vacancies, since it is impossible to generate the graphitic-N after injecting only one nitrogen atom, investigations on GO with DV are beyond our present research. Therefore, in the following discussions our main goal is to study the reduction and doping processes of GO with SV to uncover the formation of graphitic-N and to clarify the role of defects in carbon skeleton on chemical doping. Herein, the reactions between ammonia and defective GO with a carbonyl group were studied kinetically. Firstly, optimizations of two different models with single carbonyl group were performed and the results show that the structure GO-2 shown in Figure 2 is 17.6 kcal/mol lower in energy than the other model (labeled as GO-2b, shown in Figure S1 in supporting information). Therefore GO-2 was employed as the precursor for further kinetic computations. The optimized structures and calculated reaction channels are shown in Figure 7-10. After adsorption of NH3 to generate the identical intermediate (INT1-ghi in Figure 7-10) in the first step, three distinct pathways based on different addition sites of NH3were discussed. The calculation results indicate that NH3 prefers to attack from the pentagonal site (C2, in route-g) to generate a stable doping product with graphitic-N (NG-H-1 in Figure 7) after N-H dissociations and water formation. Route-g is the only one exothermic reaction with a negative reaction energy of -55.5 kcal mol-1 and exhibits the lowest energy barrier of 55.3 kcal mol-1 among all possible reaction pathways. Such an activation barrier is much lower than that of the abovementioned nitrogen doping process from a perfect graphene layer (c. a. 110.0 kcal mol -1), implying that structural defect is essential for nitrogen doping to form graphitic-N. After reduction of carbonyl by transferring two hydrogen atoms from ammonia (NG-H-1), the single carbon vacancy on carbon surface is healed by the external nitrogen. Figure 7

ACS Paragon Plus Environment

8

Page 9 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 8 Figure 9 Figure 10 As depicted in Figure 10, in comparison with route-g, the energy barriers of two other channels are much higher. In the case of route-h, NH3firstly attacks the carbon next to the pentagon (C3) which exhibits the shortest distance to foreign nitrogen in INT1-ghi of 3.30 Å. Since the defective five-membered ring in GO-2 is not destroyed during nitrogen doping, the final doping network contains a defective pentagon-heptagon-pentagon adjacency (NG-H-2 in Figure 8). As seen in Figure 8 and 10, the much higher energy barrier (77.5 kcal mol-1) as well as the positive reaction energy (2.1 kcal mol-1) both suggest that the formation of this pentagon and heptagon-containing hetero-configuration is energetically unfavorable in practice. In the case of route-i, it is evident from Figure 9 and 10 that the reaction cannot embed the nitrogen inside the network of graphene to generate the graphitic-N. Moreover, the reaction needs to overcome a considerably huge barrier of 109.7 kcal mol-1, implying that route-i is kinetically unfavorable. Therefore, according to the abovementioned three reaction routes in Figures 7-9, route-g is considered as the dominant mechanism for C=O reduction and nitrogen doping, both thermodynamically and kinetically. Besides carbonyl, hydroxyl is another important oxygen-containing group in GO. It should be mentioned that if only one –OH is participated, the reduction products must be a water molecule and an adsorbed –NH2 group. The two N-H bonds in –NH2 are stable and obviously N-H dissociation cannot take place for further nitrogen doping. Therefore, the nitrogen doping with the presence of –OH should be considered as a multi-hydroxyl participated proceed. In the

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 32

present work, the reactions between NH3 and defective GO with two hydroxyl groups were investigated. As shown in Figure 2 and Figure S2 in supporting information, the two models GO3 and GO-4 which are proved as the most two stable GO with two hydroxylgroups among all possible configurations were established as the precursors. GO-3 is 0.2 kcal/mol lower than GO4 in energy and at least 1.2 kcal/mol lower than other species. In GO-3 and GO-4, one of the two hydroxyls is located on the dangling carbon and the other one is adsorbed around the vacancy. A hydrogen bond is formed between the two –OH groups for stabilization. In the first step, the foreign NH3 is adsorbed on GO. The NH3 molecule can both locate on the carbon surface and also the –OH group. As shown in Figure 11 and 12, the location of NH3 on the –OH group by hydrogen bonding interaction (INT1-j and INT1-k) is exothermic by -14.9 kcal mol-1 in GO-3 (-12.5 kcal mol-1 in GO-4). In contrast with the barrierless adsorption of NH3 to hydroxyl, the adsorption of NH3 to carbon surface in Figures 13 and 14 (INT1-l and INT1-m) is demonstrated as an endothermic reaction with an activation barrier of 27.6 kcal mol -1 in GO-3 (30.1 kcal mol-1 in GO-4), suggesting that the NH3molecule is preferred to adsorb to –OH rather than the carbon surface. Therefore, we further investigated the following reduction and doping reactions based on the structures of INT1-j and INT-k, as shown in Figure 11. In the presence of adsorbed NH3, the dehydroxylation from two –OH groups wascalculatedsubsequently. The reaction barrier of hydrogen transfer is 10.9 kcal mol-1 in GO-3 and 12.8 kcal mol-1 in GO-4, respectively. Interestingly, we also calculated the proton transfer between –OH groups without the participation of NH3 for comparison. As shown in Figures 15 and 16, the energy barrier is found to be slightly increased to 12.7 kcal mol-1 for GO-3 and 15.0 kcal mol-1 for GO-4, which indicates that the adsorbed NH3 can promote the hydrogen transfer to form an eliminated NH3H2O complex. After the removal of NH3-H2O (GO-2 in Figure 11), the –OH group which locates

ACS Paragon Plus Environment

10

Page 11 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

on the dangling C6 is reduced as a carbonyl group. In another word, GO-2 is obtained from GO3 and GO-4 by intramolecular dehydroxylation. As mentioned above, the C=O group in GO-2 can react with NH3 to dissociate the N-H bonds to generate a water molecule and -NH species (route-g in Figure 8). The nitrogen atom in NH3 is successfully inserted into the graphene surface to form a graphitic-N. Throughout the whole reaction processes (route-j for dehydroxylation and route-g for decarbonylation), dehydroxylation from two –OH is much easier and the ratedetermining step is the reduction of C=O and nitrogen doping in subsequent steps. Figure 11 Figure 12 Figure 13 Figure 14 Figure 15 Figure 16 Kinetically, on the basis of the Arrhenius formula, the calculated reaction rate of route-g at 900 ℃ is 2.26 × 106 s-1 (with a barrier of 55.3 kcal/mol and an attempt frequency of 429.5 cm-1 which is originated from the frequency analysis of rate-determining TS3-g), implying that the formation of graphitic-N from NH3 and defective GO with a carbonyl group at elevated temperature is available (the calculated reaction rates at other temperatures are listed in Table S1 in supporting information for comparison). Such a computational result is consistent with

ACS Paragon Plus Environment

11

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 32

experimental evidence32and successfully uncovers the nitrogen doping mechanisms of GO from foreign nitrogen sources for the first time. Conclusion In the present work, we studied the oxygen reduction and nitrogen doping mechanisms of GO in the presence of ammonia molecules. In the case of perfect graphene oxide without any defects (GO-1), de-epoxidation can be achieved by N-H dissociations but the following nitrogen doping is energetically unfavorable because of the endothermic character and the huge activation barrier of 107.2 kcal mol-1 for C-C cleavage. Moreover, defective graphene oxides with single vacancies were modeled to expose the formation mechanism of graphitic-N. Three distinct reaction pathways (route-g, h, i) were discovered for decarbonylation and nitrogen doping. Owing to the negative reaction energy of -55.5 kcal mol-1 and the lowest energy barrier of 55.3 kcal mol-1, route-g is predicted as the most favorable pathway for doping. The calculated reaction rate of route-g at 900 ℃ (2.26 × 106 s-1) indicates that the formation of graphitic-N from NH3 and defective GO with a carbonyl group (GO-2) is feasible at elevated temperature. Furthermore, we discussed the reaction of NH3 to defective GO with two hydroxyl groups (GO-3 and GO-4). Our calculation results suggest that the NH3 molecule is preferred to adsorbed with the hydroxyl group instead of the carbon surface and the reactivity of dehydroxylation between two hydroxyls is improved in the presence of adsorbed NH3 to generate the abovementioned defective GO-2 with a single carbonyl group. Similarly, subsequent nitrogen doping from NH3 can take place and the graphitic-N can be generated by route-g. Overall, we introduced the detailed formation mechanisms of nitrogen-doped graphene materials from various GO and foreign nitrogen sources for the first time. The present work

ACS Paragon Plus Environment

12

Page 13 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

suggests that not only the reducing agent but also the defects in carbon skeleton are essential for reduction and doping, which is useful to understand the chemical doping process of graphenebased materials for the design of hyperfine functionalized materials in the future. AUTHOR INFORMATION Corresponding Authors * Fax: +86-29-82668559; Tel: +86-29-82665671; E-mail: [email protected]; [email protected]. Author Contributions The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript. Notes Any additional relevant notes should be placed here. Acknowledgments This work has been financially supported by the National Natural Science Foundation of China (21171138, 21573172) and the Specially Promoted Research Grant (22000009) from the Ministry of Education, Culture, Sports, Science, and Technology of Japan. ASSOCIATED CONTENT Supporting Information

ACS Paragon Plus Environment

13

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 32

Structures and relative energies of single carbonyl and two-hydroxylGO configurations, calculated reaction ratesat different temperatures, and full citations of refs. 4, 11, 48, and 53. This material is available free of charge via the Internet at http://pubs.acs.org.

References 1.

Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Jiang, D.; Zhang, Y.; Dubonos, S. V.;

Grigorieva, I. V.; Firsov, A. A. Electric Field Effect in Atomically Thin Carbon Films. Science2004, 306, 666–669. 2.

Stankovich, S.; Dikin, D. A.; Dommett, G. H. B.; Kohlhaas, K. M.; Zimney, E. J.; Stach,

E. A.; Piner, R. D.; Nguyen S. T.; Ruoff, R. S. Graphene-based Composite Materials. Nature2006, 442, 282–286. 3.

Yoo, E.; Kim, J.; Hosono, E.; Zhou, H.; Kudo, T.; Honma, I. Large Reversible Li Storage

of Graphene Nanosheet Families for Use in Rechargeable Lithium Ion Batteries. Nano Lett.2008, 8, 2277–2282. 4.

Wang, D.; Choi, D.; Li, J.; Yang, Z.; Nie, Z.; Kou, R.; Hu, D.; Wang, C.; Saraf, L. V.;

Zhang, J.; et al. Self-Assembled TiO2–Graphene Hybrid Nanostructures for Enhanced Li-Ion Insertion. ACS Nano2009, 3, 907–914. 5.

Wang, G.; Shen, X.; Yao J.; Park, J. Graphene Nanosheets for Enhanced Lithium Storage

in Lithium Ion Batteries. Carbon2009, 47, 2049–2053.

ACS Paragon Plus Environment

14

Page 15 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

6.

Eda, G.; Lin, Y. Y.; Miller, S.; Chen, C. W.; Su W. F.; Chhowalla, M. Transparent and

Conducting Electrodes for Organic Electronics from Reduced Graphene Oxide. Appl. Phys. Lett. 2008, 92, 233305. 7.

Wang, X.; Zhi, L.; Mullen, K. Transparent, Conductive Graphene Electrodes for Dye-

Sensitized Solar Cells. Nano Lett. 2008, 8, 323–327. 8.

Wu, J.; Becerril, H. A.; Bao, Z.; Liu, Z.; Chen, Y.; Peumans, P. Organic Solar Cells with

Solution-processed Graphene Transparent Electrodes. Appl. Phys. Lett. 2008, 92, 263302. 9.

Yin, Z.; Wu, S.; Zhou, X.; Huang, X.; Zhang, Q.; Boey, F.; Zhang, H. Electrochemical

Deposition of ZnO Nanorods on Transparent Reduced Graphene Oxide Electrodes for Hybrid Solar Cells. Small2010, 6, 307–312. 10.

Ponomarenko, L. A.; Schedin, F.; Katsnelson, M. I.; Yang, R.; Hill, E. W.; Novoselov, K.

S.; Geim, A. K. Chaotic Dirac Billiard in Graphene Quantum Dots. Science2008, 320, 356–358. 11.

Wang, S.; Wang, J.; Miraldo, P.; Zhu, M.; Outlaw, R.; Hou, K.; Zhao, X.; Holloway, B.

C.; Manos, D.; Tyler, T.; et al. High Field Emission Reproducibility and Stability of Carbon Nanosheets and Nanosheet-based Backgated Triode Emission Devices. Appl. Phys. Lett. 2006, 89, 183103. 12.

Chen, H.; Muller, M. B.; Gilmore, K. J.; Wallace, G. G.; Li, D. Mechanically Strong,

Electrically Conductive, and Biocompatible Graphene Paper.Adv. Mater. 2008, 20, 3557–3561. 13.

Wang, L.; Lee, K.; Sun, Y. Y.; Lucking, M.; Chen, Z.; Zhao, J. J.; Zhang, S. B. Graphene

Oxide as an Ideal Substrate for Hydrogen Storage. ACS Nano2009, 3, 2995–3000.

ACS Paragon Plus Environment

15

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

14.

Page 16 of 32

Liu, C.; Yu, Z.; Neff, D.; Zhamu, A.; Jang, B. Z. Graphene-Based Supercapacitor with an

Ultrahigh Energy Density. Nano Lett. 2010, 10, 4863–4868. 15.

Mohanty, N.; Berry, V. Graphene-Based Single-Bacterium Resolution Biodevice and

DNA Transistor: Interfacing Graphene Derivatives with Nanoscale and Microscale Biocomponents. Nano Lett. 2008, 8, 4469–4476. 16.

Jiang, Z.; Wang, J.; Meng, L.; Huang, Y.; Liu, L. A Highly Efficient Chemical Sensor

Material for Ethanol: Al2O3/Graphene Nanocomposites Fabricated from Graphene Oxide. Chem. Commun. 2011, 47, 6350–6352. 17.

Mao, S.; Lu, G.; Yu, K.; Bo, Z.; Chen, J. Specific Protein Detection Using Thermally

Reduced Graphene Oxide Sheet Decorated with Gold Nanoparticle-Antibody Conjugates. Adv. Mater. 2010, 22, 3521–3526. 18.

Dong, X.; Shi, Y.; Huang, W.; Chen, P.; Li, L. J. Electrical Detection of DNA

Hybridization with Single-Base Specificity Using Transistors Based on CVD-Grown Graphene Sheets.Adv. Mater. 2010, 22, 1649. 19.

Nguyen,P.;Berry,V.Graphene Interfaced with Biological Cells: Opportunities and

Challenges. J. Phys. Chem. Lett. 2012, 3, 1024–1029. 20.

Wang, X.; Li, X.; Zhang, L.; Yoon, Y.; Weber, P. K.; Wang, H.; Guo, J.; Dai, H. N-

Doping of Graphene Through Electrothermal Reactions with Ammonia. Science2009, 324, 768– 771.

ACS Paragon Plus Environment

16

Page 17 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

21.

Reddy, A. L. M.; Srivastava, A.; Gowda, S. R.; Gullapalli, H.; Dubey, M.; Ajayan, P. M.

Synthesis of Nitrogen-doped Graphene Films for Lithium Battery Application. ACS Nano2010, 4, 6337–6342. 22.

Jin, Z.; Yao, J.; Kittrell, C.; Tour, J. M. Large-Scale Growth and Characterizations of

Nitrogen-Doped Monolayer Graphene Sheets. ACS Nano2011, 5, 4112–4117. 23.

Guo, B.; Liu, Q.; Chen, E.; Zhu, H.; Fang, L.; Gong, J. R. Controllable N-Doping of

Graphene. Nano Lett. 2010, 10, 4975–4980. 24.

Wang, Y.; Shao, Y.; Matson, D. W.; Li, J.; Lin, Y. Nitrogen-Doped Graphene and Its

Application in Electrochemical Biosensing. ACS Nano2010, 4, 1790–1798. 25.

Jeong, H. M.; Lee, J. W.; Shin, W. H.; Choi, Y. J.; Shin, H. J.; Kang, J. K.; Choi, J. W.

Nitrogen-Doped Graphene for High-Performance Ultracapacitors and the Importance of Nitrogen-Doped Sites at Basal Planes. Nano Lett. 2011, 11, 2472–2477. 26.

Wu, Z. S.; Ren, W.; Xu, L.; Li, F.; Cheng, H. M. Doped Graphene Sheets As Anode

Materials with Superhigh Rate and Large Capacity for Lithium Ion Batteries. ACS Nano2011, 5, 5463–5471. 27.

Wei, D.; Liu, Y.; Wang, Y.; Zhang, H.; Huang, L.; Yu, G. Synthesis of N-Doped

Graphene by Chemical Vapor Deposition and Its Electrical Properties. Nano Lett. 2009, 9, 1752– 1758. 28.

Qu, L.; Liu, Y.; Baek, J. B.; Dai, L. Nitrogen-Doped Graphene as Efficient Metal-Free

Electrocatalyst for Oxygen Reduction in Fuel Cells. ACS Nano2010, 4, 1321–1326.

ACS Paragon Plus Environment

17

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

29.

Page 18 of 32

Shao, Y.; Zhang, S.; Engelhard, M. H.; Li, G.; Shao, G.; Wang, Y.; Liu, J.; Aksay, I. A.;

Lin, Y. Nitrogen-doped Graphene and Its Electrochemical Applications. J. Mater. Chem. 2010, 20, 7491–7496. 30.

Geng, D.; Chen, Y.; Chen, Y.; Li, Y.; Li, R.; Sun, X.; Ye, S.; Knights, S. High Oxygen-

reduction Activity and Durability of Nitrogen-doped Graphene. Energy Environ. Sci. 2011, 4, 760–764. 31.

Sheng, Z. H.; Shao, L.; Chen, J. J.; Bao, W. J.; Wang, F. B.; Xia, X. H. Catalyst-free

Synthesis of Nitrogen-doped Graphene via Thermal Annealing Graphite Oxide with Melamine and Its Excellent Electrocatalysis. ACS Nano2011, 5, 4350–4358. 32.

Li, X.; Wang, H.; Robinson, J. T.; Sanchez, H.; Diankov, G.; Dai, H. Simultaneous

Nitrogen Doping and Reduction of Graphene Oxide. J. Am. Chem. Soc. 2009, 131, 15939– 15944. 33.

Imran

Jafri,

R.;

Rajalakshmi,

Nitrogen Doped Graphene Nanoplateletsas Catalyst Support

N.;

Ramaprabhu,

for Oxygen

S.

Reduction Reaction

in Proton Exchange Membrane Fuel Cell. J. Mater. Chem. 2010, 20, 7114–7117. 34.

Lin, Y. C.; Lin, C. Y.; Chiu, P. W. Controllable Graphene N-doping with Ammonia

Plasma. Appl. Phys. Lett. 2010, 96, 133110. 35.

Zhang, L. S.; Liang, X. Q.; Song, W. G.; Wu, Z. Y. Identification of the Nitrogen Species

on N-doped Graphene Layers and Pt/NG Composite Catalyst for Direct Methanol Fuel Cell. Phys. Chem. Chem. Phys. 2010, 12, 12055–12059.

ACS Paragon Plus Environment

18

Page 19 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

36.

Parvez, K.; Yang, S.; Hernandez, Y.; Winter, A.; Turchanin, A.; Feng, X.; Müllen, K.

Nitrogen-Doped Graphene and Its Iron-Based Composite As Efficient Electrocatalysts for Oxygen Reduction Reaction. ACS Nano2012, 6, 9541–9550. 37.

Liang, Y.; Wu, D.; Feng, X.; Müllen, K. Dispersion of Graphene Sheets in Organic

Solvent Supported by Ionic Interactions. Adv. Mater. 2009, 21, 1679–1683. 38.

Paredes, J. I.; Villar-Rodil, S.; Solís-Fernández, P.; Martínez-Alonso, A.; Tascón, J. M.

D. Atomic Force and Scanning Tunneling Microscopy Imaging of Graphene Nanosheets Derived from Graphite Oxide. Langmuir2009, 25, 5957–5968. 39.

Li, J. L.; Kudin, K. N.; McAllister, M. J.; Prud’homme, R. K.; Aksay, I. A.; Car, R.

Oxygen-Driven Unzipping of Graphitic Materials. Phys. Rev. Lett. 2006, 96, 176101. 40.

Gómez-Navarro, C.; Weitz, R. T.; Bittner, A. M.; Scolari, M.; Mews, A.; Burghard, M.;

Kern, K. Electronic Transport Properties of Individual Chemically Reduced Graphene Oxide Sheets. Nano Lett. 2007, 7, 3499–3503. 41.

Mkhoyan, K. A.; Contryman, A. W.; Silcox, J.; Stewart, D. A.; Eda, G.; Mattevi, C.;

Miller, S.; Chhowalla, M. Atomic and Electronic Structure of Graphene-Oxide. Nano Lett. 2009, 9, 1058–1063. 42.

Frankcombe, T. J.; Bhatia, S. K.; Smith, S. C. Ab initio Modelling of Basal Plane

Oxidation of Graphenes and Implications for Modelling Char Combustion. Carbon2002, 40, 2341–2349. 43.

Sendt, K.; Haynes, B. S. Density Functional Study of the Chemisorption of O2 Across

Two Rings of the Armchair Surface of Graphite. J. Phys. Chem. C2007, 111, 5465–5473.

ACS Paragon Plus Environment

19

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

44.

Page 20 of 32

Tachikawa, H.; Kawabata, H. Electronic States of Defect Sites of Graphene Model

Compounds: A DFT and Direct Molecular Orbital−Molecular Dynamics Study. J. Phys. Chem. C2009, 113, 7603–7609. 45.

Gao, X.; Jang, J.; Nagase, S. Hydrazine and Thermal Reduction of Graphene Oxide:

Reaction Mechanisms, Product Structures, and Reaction Design. J. Phys. Chem. C2010, 114, 832–842. 46.

Boukhvalov, D. W.; Katsnelson, M. I. Modeling of Graphite Oxide. J. Am. Chem. Soc.

2008, 130, 10697–10701. 47.

Li, Z.; Zhang, W.; Luo, Y.; Yang, J.; Hou, J. How Graphene is Cut upon Oxidation? J.

Am. Chem. Soc. 2009, 131, 6320–6321. 48.

Cai, W.; Piner, R. D.; Stadermann, F. J.; Park, S.; Shaibat, M. A.; Ishii, Y.; Yang, D.;

Velamakanni, A.; An, S. J.; Stoller, M.; et al. Synthesis and Solid-State NMR Structural Characterization of 13C-Labeled Graphite Oxide. Science2008, 321, 1815–1817. 49.

Szabó, T.; Tombácz, E.; Illés, E.; Dékány, I. Enhanced Acidity and pH-dependent

Surface Charge Characterization of Successively Oxidized Graphite Oxides. Carbon2006, 44, 537–545. 50.

He, H.; Klinowski, J.; Forster, M.; Lerf, A. A New Structural Model for Graphite Oxide.

Chem. Phys. Lett. 1998, 287, 53–56. 51.

He, H.; Riedl, T.; Lerf, A.; Klinowski, J. Solid-State NMR Studies of the Structure of

Graphite Oxide. J. Phys. Chem. 1996, 100, 19954–19958.

ACS Paragon Plus Environment

20

Page 21 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

52.

Szabó, T.; Berkesi, O.; Forgo, P.; Josepovits, K.; Sanakis, Y.; Petridis, D.; Dékány, I.

Evolution of Surface Functional Groups in a Series of Progressively Oxidized Graphite Oxides. Chem. Mater. 2006, 18, 2740–2749. 53.

Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.;

Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; et al. Gaussian 09 Revision A.02

(Gaussian, Inc., Wallingford CT, 2009). 54

Han, S.; Yoon, M.; Berber, S.; Park, N.; Osawa, E.; Ihm, J.; Tomanek, D.Microscopic

Mechanism of Fullerene Fusion. Phys. Rev. B, 2004, 70, 113402. 55.

Tang, S.; Cao, Z. Adsorption and Dissociation of Ammonia on Graphene Oxides: A

First-Principles Study. J. Phys. Chem. C 2012, 116, 8778–8791. 56.

Meyer, J. C.; Kisielowski, C.; Erni, R.; Rossell, M. D.; Crommie, M. F.; Zettl, A. Direct

Imaging of Lattice Atoms and Topological Defects in Graphene Membranes. Nano Lett. 2008, 8, 3582–3586. 57.

Warner, J. H.; Rümmeli, M. H.; Ge, L.; Gemming, T.; Montanari, B.; Harrison, N. M.;

Büchner, B.; Briggs, G. A. D. Structural Transformations in Graphene Studied with High Spatial and Temporal Resolution. Nat. Nanotechnol. 2009, 4, 500–504. 58.

Ugeda, M. M.; Brihuega, I.; Guinea, F.; Gómez-Rodríguez, J. M. Missing Atom as a

Source of Carbon Magnetism. Phys. Rev. Lett. 2010, 104, 096804. 59.

Krasheninnikov, A. V.; Lehtinen, P. O.; Foster, A. S.; Nieminen, R. M. Bending the

Rules: Contrasting Vacancy Energetics and Migration in Graphite and Carbon Nanotubes. Chem. Phys. Lett. 2006, 418, 132–136.

ACS Paragon Plus Environment

21

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

60.

Page 22 of 32

El-Barbary, A. A.; Telling, R. H.; Ewels, C. P.; Heggie, M. I.; Briddon, P. R. Structure

and Energetics of the Vacancy in Graphite. Phys. Rev. B2003, 68, 144107. 61.

Lee, G. -D.; Wang, C. Z.; Yoon, E.; Hwang, N.-M.; Kim, D.-Y.; Ho, K. M. Diffusion,

Coalescence, and Reconstruction of Vacancy Defects in Graphene Layers. Phys. Rev. Lett. 2005, 95, 205501. 62.

Banhart, F.; Kotakoski, J.; Krasheninnikov, A. V. Structural Defects in Graphene. ACS

Nano2011, 5, 26–41.

ACS Paragon Plus Environment

22

Page 23 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 1.Schematic representation of nitrogen-doped graphene.

Figure 2.Computational models of graphene oxides. Atoms representations are O (red), C (white) and H (pink).GO-1: perfect graphene with an epoxide group. GO-2: defective graphene with a carbonyl group in the center. GO-3 and GO-4: defective graphene with two hydroxyl groups.

ACS Paragon Plus Environment

23

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 32

Figure 3.Fragmental structures for stationary points involved in route-a and route-b. Atoms representations are O (red), N (blue), C (white) and H (pink). Other carbon atoms are omitted for clarity. The distances are represented in units of Å.Relative energies (in parenthesis) are given in kcal mol-1.

Figure 4.Reaction energy profiles for ammonia participated reduction of GO-1 (route-a in red and route-b in green). Relative energies referred to the total energy of reactants (GO-1+NH3) are calculated at the M06-2X/6-31G(d) level of theory, units in kcal mol-1.

ACS Paragon Plus Environment

24

Page 25 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 5.Fragmental structures for stationary points involved in route-c, d, e, and f. Atoms representations are O (red), N (blue), C (white) and H (pink).Other carbon atoms are omitted for clarity. The distances are represented in units of Å.Relative energies (in parenthesis) are given in kcal mol-1.

ACS Paragon Plus Environment

25

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 32

Figure 6.Reaction energy profiles for nitrogen doping of GO-1 (route-c in green, route-d in red, route-e in pink, and route-f in blue). Relative energies referred to the energy of reactant (G-NH1) are calculated at the M06-2X/6-31G(d) level of theory, units in kcal mol-1.

Figure 7.Fragmental structures for stationary points involved in route-g. Atoms representations are O (red), N (blue), C (white) and H (pink).Other carbon atoms are omitted for clarity. The distances are represented in units of Å.Relative energies (in parenthesis) are given in kcal mol-1.

ACS Paragon Plus Environment

26

Page 27 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 8.Fragmental structures for stationary points involved in route-h. Atoms representations are O (red), N (blue), C (white) and H (pink). Other carbon atoms are omitted for clarity. The distances are represented in units of Å.Relative energies (in parenthesis) are given in kcal mol-1.

Figure 9.Fragmental structures for stationary points involved in route-i. Atoms representations are O (red), N (blue), C (white) and H (pink).Other carbon atoms are omitted for clarity. The distances are represented in units of Å.Relative energies (in parenthesis) are given in kcal mol-1.

ACS Paragon Plus Environment

27

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 32

Figure 10.Reaction energy profiles for reactions of ammonia with GO-2 (route-g in red, route-h in blue, and route-i in green). Relative energies referred to the total energy of reactants (GO2+NH3) are calculated at the M06-2X/6-31G(d) level of theory, units in kcal mol-1.

Figure 11.Fragmental structures for stationary points involved in route-j and route-k. Atoms representations are O (red), N (blue), C (white) and H (pink).Other carbon atoms are omitted for clarity. The distances are represented in units of Å.Relative energies (in parenthesis) are given in kcal mol-1.

ACS Paragon Plus Environment

28

Page 29 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 12.Reaction energy profiles for ammonia participated dehydroxylation of GO-3 (route-j in red) and GO-4 (route-k in green). Relative energies referred to the total energy of reactants (GO-3+NH3) are calculated at the M06-2X/6-31G(d) level of theory, units in kcal mol-1.

Figure 13.Fragmental structures for stationary points involved in route-l and route-m. Atoms representations are O (red), N (blue), C (white) and H (pink).Other carbon atoms are omitted for clarity. The distances are represented in units of Å.Relative energies (in parenthesis) are given in kcal mol-1.

ACS Paragon Plus Environment

29

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 30 of 32

Figure 14.Reaction energy profiles for ammonia participated dehydroxylation of GO-3 (route-l in red) and GO-4 (route-m in green). Relative energies referred to the energy of reactant GO-3 are calculated at the M06-2X/6-31G(d) level of theory, units in kcal mol-1.

Figure 15.Fragmental structures for stationary points involved in route-n and route-o. Atoms representations are O (red), C (white) and H (pink).Other carbon atoms are omitted for clarity. The distances are represented in units of Å.Relative energies (in parenthesis) are given in kcal mol-1.

ACS Paragon Plus Environment

30

Page 31 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 16.Reaction energy profiles for intramolecular dehydroxylation of GO-3 (route-n in red) and GO-4 (route-o in green). Relative energies referred to the energy of reactant GO-3 are calculated at the M06-2X/6-31G(d) level of theory, units in kcal mol-1.

ACS Paragon Plus Environment

31

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 32

Insert Table of Contents Graphic and Synopsis Here

ACS Paragon Plus Environment

32