Skip to content
BY-NC-ND 3.0 license Open Access Published by De Gruyter Open Access December 30, 2016

The Interaction between Graphene and Oxygen Atom

Yifan Hao , Xuejun Zhao , Xuedan Song , Hongjiang Li , Xiaobing Zhu and Ce Hao EMAIL logo
From the journal Open Physics


Based on the density function theory (DFT) method, the interaction between the graphene and oxygen atom is simulated by the B3LYP functional with the 6-31G basis set. Due to the symmetry of graphene (C54H18, D6h), a representative patch is put forward to represent the whole graphene to simplify the description. The representative patch on the surface is considered to gain the potential energy surface (PES). By the calculation of the PES, four possible stable isomers of the C54H18-O radical can be obtained. Meanwhile, the structures and energies of the four possible stable isomers, are further investigated thermodynamically, kinetically, and chemically. According to the transition states, the possible reaction mechanism between the graphene and oxygen atom is given.

PACS: 31.50.Bc

1 Introduction

Since being discovered in 2004 [1], graphene, as the first truly two-dimensional crystal, with large surface area, high chemical stability, flexibility, superior electric and thermal conductivity [25], has drawn dramatic attention for use in such areas as sensors, energy conversion, photodetectors, and aerospace [610]. Graphene oxide (GO) as a significant raw material of mass-production graphene via reduction, is an important carbon nanomaterial in the graphene family [1114]. Because GO is rich for oxygen-containing groups in the structure, and the ratio of sp2 to sp3 hybridized carbon atoms is easy to control, GO possesses brilliant performance in electronics, optics, energy conversion, environmental science, and biotechnology. It is reported that GO nanoplatelets were self-assembled to form carbon films with tunable porous morphologies and mechanical flexibility, which incorporated good electrical properties [15]. Yeh and his co-workers used GO in the photocatalytic of H2 evolution activity. According to the Mott-Schottky equation, the results of electrochemical analysis showed that H2 could be generated stably from the pure water or methanol solution through the GO films under mercury vapor light irradiation, even under the condition of lack of cocatalyst Pt [16]. Meanwhile, compared with other benchmark material, i.e. C60, GO exhibits a good performance on the optical power limiting. In order to investigate optical power limiting behaviors of GO, Jiang et al. [17] prepared GO films on the plastic substrates and glass, which shows excellent broadband optical limiting. Besides, compounded GO, with other materials, can enhance the thermodynamic properties of the material [1820]. In addition, graphene quantum dots (GQDs) as an emerging fluorescent nanomaterial, with excellent optical performance [21], easy to surface modify [22], high chemical stability [2326], low toxicity, and good bio-compatibility [27, 28], have drawn dramatic attention in the scientific community. All in all, GO as well as GQDs have dramatically grabbed the attention of the scientific community due to their brilliant performances. Due to most of them being prepared by oxidation, it is of vital importance to understand the interaction between the graphene and oxygen atom.

In this paper, a systematic study of the interaction between the graphene (C54H18, D6h) and oxygen atom is performed by the density function theory. Based on the symmetry, a representative patch is put forward to represent the whole graphene (C54H18, D6h) to simplify the description. By the single point energy scan on the entire surface of the graphene, the PES between the graphene (C54H18, D6h) and oxygen atom was obtained. Meanwhile, all the possible isomers, as well as all the TSs, were obtained to explain the possible reaction mechanism.

2 Computational details

In this study, C54H18 (D6h) is chosen as the whole system. Based on the symmetry of graphene, a representative patch is put forward to represent this whole graphene (C54H18, D6h) to simplify the description. In calculation, we adopted the B3LYP hybrid functional with the Gaussian09 program [2931]. The STO-3G [32], as well as the standard 6-31G(d)[33], were set as the basis sets. Firstly, graphene (C54H18, D6h) and its corresponding symmetry group, as the initial geometries are optimized with the basis set of B3LYP/6-31G (d). Secondly, in terms of the calculation of the single point energy, the potential energy curve and the PES were obtained at the basic level of STO-3G. When calculating PES, we just need to consider the energy distribution in the whole graphene (C54H18, D6h). So, the basis set of STO-3G can satisfy our requirements according to our experience. Thirdly, the isomers of C54H18-O on the representative patch were optimized respectively at the standard 6-31G (d). In addition, the harmonic vibrational frequency was analyzed to make sure that the stationary points at the same basis set, were real.

3 Results and discussion

3.1 The representative patch

Based on the symmetry, we chose a triangular section of the graphene (C54H18, D6h) as a representative patch to make it more simple(as shown in Fig. 1). All the key points have their own symmetry. For instance, the key point of R1 is C6v symmetry; the key points of Rn (n = 2, 3, 4), bn (n = 3, 4, 5, 6), C3, C4 are C2v symmetry; the key points of bn (n = 1, 2, 7, 8) C1, C2 are Cs symmetry. The area of the total surface is equal to 24 and the representative patch is the 1/24.When calculated to the representative patch, it can contain the whole surface of the graphene (C54H18, D6h). Hence, we just need to calculate the smallest structure unit (the patch DEF) to replace the entire graphene. On the representative patch, there are eighteen key points including four six-membered rings’ centres (key points of R1, R2, R3, R4), eight bond centres (key points of b1, b2, b3, b4, b5, b6, b7, b8), six carbon atoms. Because the graphene (C54H18, D6h) edge is saturated with hydrogenation, we just need to take the four carbon atom (key points of C1, C2, C3, C4) into account.

Figure 1 The representative patch of graphene (C54H18, D6h).
Figure 1

The representative patch of graphene (C54H18, D6h).

To make the potential energy surface at a reasonable interaction distance between the graphene and oxygen atom, the potential energy curves of the representative patch at the different key points were discussed by the scanning of the single point energy. Fig. 2 is the potential energy curves of the representative patch, and it can be seen that two different trends are shown from the curves. The top curves show a relatively sharp trend, which belongs to the oxygen atom moving away from the bond centre of b1 and b2. Meanwhile, the bottom curves have potential wells, the height of which is about 1.4 Å. It expressed that there is a stable chemical adsorption site at the bond centre of the b3 and b4 when the oxygen atom moved away from the graphene at the normal directions.

Figure 2 The potential energy curves of the representative patch as function of the height from oxygen atom to the key points.
Figure 2

The potential energy curves of the representative patch as function of the height from oxygen atom to the key points.

According to the potential energy curves of the C54H18–O radical, the potential energy surfaces were scanned at the height of 1.4 Å, as shown in Fig. 3. From the representative patch, it can be clearly observed that the key points of b1, b2, b3, b4 are the minimum value of the potential energy in the wells, and the key points of R1, R2, R3, and R4 are the maximum value of the potential energy. Therefore, four possible isomers are proposed theoretically, which are named C54H18 - O1C54H18 - O2C54H18 - O3C54H18 - O4 corresponding to the bond centre of bn (n = 1, 2, 3, 4).

Figure 3 The potential energy surface of C54H18-O.
Figure 3

The potential energy surface of C54H18-O.

3.2 The thermodynamic stability

Table 1 is the geometric parameters, HOMO-LUMO gap energies, and relative energies for the four possible isomers of C54H18-O. However, the balanced distance of the C-O bond (including the bonds of C-O and C-O) is 1.52Å and the mean angel of C-O-C is 60.64°.

Table 1

Geometric parameters, HOMO-LUMO gap energies and relative energies of C54H18-O.

StructureC-C (Å)C-O (Å)C-O (Å)C-O-C(°)Gap Energy(eV)Rel Energ(eV)

Among the different isomers of C54H18-O, it can be noted that C54H18-O2 is the lowest energy in total, so we set it as the zero point of the relative energy. Then C54H18-O4 is only higher by 0.003 eV than C54H18-O2, which is the second most stable. And C54H18-O3 possesses the highest energy, which is the least stable. Based on the above analysis, the ordering of thermodynamic stability is received for C54H18-O2>C54H18-O4> C54H18-O1> C54H18-O3.

3.3 The chemical reaction activity

The HOMO–LUMO energy gap sequence enable us to describe the chemical reaction activity of the four C54H18-O isomers, which can also be shown in Table 1. Comparing with the HOMO–LUMO energy gaps, we can get that C54H18-O2 is the highest gap energy at 2.728 eV, and the C54H18-O3 is the lowest gap energy at 1.932 eV. According to the gap energy of HOMO–LUMO, the sequence of the chemical reaction activity is C54H18-O2 > C54H18-O4 > C54H18-O1 > C54H18-O3, which is consistent with the sequence of the thermodynamic stability.

3.4 The kinetics activity

TSs between different C54H18-O radicals are calculated at the 6-31G (d) level. The results show that each TS only has one single imaginary frequency. TS1 represents the transition state of C54H18-O1 and C54H18-O2, and the value of single imaginary frequency is 522.24 i cm−1. TS2 represents the transition state of C54H18-O2 and C54H18-O3 with the only imaginary frequency valu of 517.55 i cm−1. TS3 represents the transition state of C54H18-O3 and C54H18-O4 and the value of single imaginary frequency is 439.86 i cm−1.

Table 2 is the activation energies (Ea) as well as the relative energy of transition states (TS1, TS2, TS3). It shows that the Ea value of the oxygen atom to be transferred from the isomer of C54H18-O1 to the isomer of C54H18-O2 is 0.747 eV. On the reverse reaction path, the Ea value of the oxygen atom to be transferred from the isomer of C54H18-O2 to the isomer of C54H18-O1 is 0.902 eV. Another reaction path for the isomer of C54H18-O2 is the oxygen atom to be transferred from the isomer of C54H18-O2 to the isomer of C54H18-O3, the Ea value of which is 1.048 eV. So the minimum Ea value of the isomer C54H18-O2 is 0.902 eV. In addition, the Ea value of C54H18-O3 (including the isomer of C54H18-O3 to the isomer of C54H18-O2 and the isomer of C54H18-O3 to the isomer of C54H18-O4) is 0.703 eV and 0.604 eV. Thus, the minimum Ea value of the isomer C54H18-O3 is 0.604 eV. Meanwhile, the Ea value of the oxygen atom to be transferred from the isomer of C54H18-O4 to the isomer of C54H18-O3 is 0.946 eV. Therefore, the stability order of the kinetics activity is C54H18-O4 > C54H18-O2 > C54H18-O1 > C54H18-O3 for the four isomers of C54H18-O. This is a little different to the sequence of the thermodynamic stability and the chemical reaction activity.

Table 2

The calculated energies and activation energy (Ea) for radical isomers C54H18-O on the representative patch at B3LYP/6-31G (d) level.

E (eV)0.1550.90201.0480.3450.9490.003
Ea (eV)0.7470.9021.0480.7030.6040.946

4 Conclusions

According to the PES gained by the scanning of the single point energy, four possible isomers of C54H18-O on the representative patch are optimized by the density functional theory method. In short, the oxygen atom is easy to absorb above the bond to form the epoxy. C54H18-O2 is the most stable isomer, calculated by analyzing the the thermodynamics and chemical reaction activity. C54H18-O4 is the most stable isomer as calculated by analyzing the kinetics activity of radicals. Furthermore, investigation of the interaction between the oxygen atom and graphene is also crucial to further understanding the details of the GO and GQDs.

Tel./fax: +86 411 84986335


The financial support of the National Natural Science Foundation of China (Grant Nos. 21373042), the Fundamental Research Funds for the Central Universities of China (Grant No. DUT15RC(4)11) and Liaoning Provincial Doctor Startup Foundation of China (Grant No. 20141019). And the Project Sponsored by the Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education Ministry is greatly appreciated.


[1] Novoselov K.S., Geim A.K., Morozov S.V., Jiang D., Zhang Y., Dubonos S.V., et al., Electric field effect in atomically thin carbon films, Science, 2004, 306, 666-669.10.1126/science.1102896Search in Google Scholar PubMed

[2] Geim A.K., Novoselov K.S., The rise of graphene, Nat. Mater., 2007, 6, 183-191.10.1142/9789814287005_0002Search in Google Scholar

[3] Shahil K.M.F., Balandin A.A., Thermal properties of graphene and multilayer graphene: applications in thermal interface materials, Solid State Commun., 2012, 152, 1331-1340.10.1016/j.ssc.2012.04.034Search in Google Scholar

[4] Stoller M.D., Park S., Zhu Y., An J., Ruoff R.S., Graphene-based ultracapacitors, Nano Lett., 2008, 8, 3498-3502.10.1021/nl802558ySearch in Google Scholar PubMed

[5] Zhu Y., Murali S., Stoller M.D., Ganesh K.J., Cai W., Ferreira P.J., et al., Carbon-based supercapacitors produced by activation of graphene, Science, 2011, 332, 1537-1541.10.1126/science.1200770Search in Google Scholar PubMed

[6] Bandurin D.A., Torre I., Kumar R.K., Shalom M.B., Tomadin A., Principi A., et al., Negative local resistance caused by viscous electron backflow in graphene, Science, 2016, 351, 1055-1058.10.1126/science.aad0201Search in Google Scholar PubMed

[7] Anguita J.V., Ahmad M., Haq S., Allam J., Silva S.R P., Ultra-broadband light trapping using nanotextured decoupled graphene multilayersal, Sci. Adv., 2016, 2, e1501238.10.1126/sciadv.1501238Search in Google Scholar PubMed PubMed Central

[8] Katsnelson M.I., Graphene: carbon in two dimensions, Mater. Today, 2007, 10, 20-27.10.1017/CBO9781139031080Search in Google Scholar

[9] Zhang M.X., Zhao G.J., Modification of n-type organic semiconductor performance of perylene diimides by substitution in different positions: two-dimensional pi-stacking and hydrogen bonding, ChemSusChem, 2012, 5, 879-887.10.1002/cssc.201100510Search in Google Scholar PubMed

[10] Neto A.H.C., Guinea F., Peres N.M.R., Novoselov K.S., Geim A.K., The electronic properties of graphene, Rev. Mod. Phys., 2009, 81, 109-162.10.1103/RevModPhys.81.109Search in Google Scholar

[11] Wei Z., Wang D., Kim S., Kim S.Y., Hu Y., Yakes M.K., et al., Nanoscale tunable reduction of graphene oxide for graphene electronics, Science, 2010, 328, 1373-1376.10.1126/science.1188119Search in Google Scholar PubMed

[12] Li H., Song Z., Zhang X., Huang Y., Li S., Mao Y., et al., Ultrathin, molecular-sieving graphene oxide membranes for selective hydrogen separation, Science, 2013, 342, 95-98.10.1126/science.1236686Search in Google Scholar PubMed

[13] Binggeli M., Christoph R., Hintermann H.E., Colchero J., Marti O., Friction force measurements on potential controlled graphite in an electrolyte environment, Nanotechnology, 1993, 4, 59-63.10.1088/0957-4484/4/2/001Search in Google Scholar

[14] Nair R.R., Wu H.A., Jayaram P.N., Grigorieva I.V., Geim A.K., Unimpeded permeation of water through helium-leak-tight graphene-based membranes, Science, 2012, 335, 442-444.10.1126/science.1211694Search in Google Scholar PubMed

[15] Lee S.H., Kim H.W., Hwang J.O., Lee W.J., Kwon J., Bielawski C.W., et al., Three-dimensional self-assembly of graphene oxide platelets into mechanically flexible macroporous carbon films, Angew. Chem., 2010, 122, 10282-10286.10.1002/ange.201006240Search in Google Scholar

[16] Yeh T.F., Syu J.M., Cheng C., Chang T.H., Teng H., Graphite oxide as a photocatalyst for hydrogen production from water, Adv. Funct. Mater., 2010, 20, 2255-2262.10.1002/adfm.201000274Search in Google Scholar

[17] Jiang X.F., Polavarapu L., Neo S.T., Venkatesan T., Xu Q.H., Graphene oxides as tunable broadband nonlinear optical materials for femtosecond laser pulses, J. Phys. Chem. Lett., 2012, 3, 785-790.10.1021/jz300119tSearch in Google Scholar PubMed

[18] Zhao S., Li Y., Yin H., Liu Z., Luan E., Zhao F., et al., Three-dimensional graphene/Pt nanoparticle composites as freestanding anode for enhancing performance of microbial fuel cells, Sci. Adv., 2015, 1, e1500372.10.1126/sciadv.1500372Search in Google Scholar PubMed PubMed Central

[19] Stankovich S., Dikin D.A., Dommett G.H.B., Kohlhaas K.M., Zimney E.J., Stach E.A., et al., Graphene-based composite materials, Nature, 2006, 442, 282-286.10.1038/nature04969Search in Google Scholar PubMed

[20] Dikin D.A., Stankovich S., Zimney E.J., Piner R.D., Dommett G.H.B., Evmenenko G., et al., Preparation and characterization of graphene oxide paper, Nature, 2007, 448, 457-460.10.1038/nature06016Search in Google Scholar PubMed

[21] Zhu S., Meng Q., Wang L., Zhang J., Song Y., Jin H., et al., Highly photoluminescent carbon dots for multicolor patterning, sensors, and bioimaging, Chem. Int. Ed., 2013, 52, 3953-3957.10.1002/anie.201300519Search in Google Scholar PubMed

[22] Sun Y.P., Zhou B., Lin Y., Wang W., Fernando K.A.S., Pathak P., et al., Quantum-sized carbon dots for bright and colorful photoluminescence, J. Am. Chem. Soc., 2006, 128, 7756-7757.10.1021/ja062677dSearch in Google Scholar PubMed

[23] Pei S., Zhang J., Gao M., Wu D., Yang Y., Liu R., A facile hydrothermal approach towards photoluminescent carbon dots from amino acids, J. Colloid Interface Sci., 2015, 439, 129-133.10.1016/j.jcis.2014.10.030Search in Google Scholar PubMed

[24] Kong B., Zhu A., Ding C., Zhao X., Li B., Tian Y., Carbon Dot-Based Inorganic–organic nanosystem for two-photon imaging and biosensing of pH variation in living cells and tissues, Adv. Mater., 2012, 24, 5844-5848.10.1002/adma.201202599Search in Google Scholar PubMed

[25] Han B., Wang W., Wu H., Fang F., Wang N., Zhang X., Polyethyleneimine modified fluorescent carbon dots and their application in cell labeling, Colloids Surf. B, 2012, 100, 209-214.10.1016/j.colsurfb.2012.05.016Search in Google Scholar PubMed

[26] Ruan S., Zhu B., Zhang H., Chen J., Shen S., Qian J., et al., A simple one-step method for preparation of fluorescent carbon nanospheres and the potential application in cell organelles imaging, J. Colloid. Interface Sci., 2014, 422, 25-29.10.1016/j.jcis.2014.02.006Search in Google Scholar PubMed

[27] Yang S.T., Wang X., Wang H., Lu F., Luo P.G., Cao L., et al., Carbon dots as nontoxic and high-performance fluorescence imaging agents, J. Phys. Chem. C, 2009, 113, 18110-18114.10.1021/jp9085969Search in Google Scholar PubMed PubMed Central

[28] Ge S., Ge L., Yan M., Song X., Yu J., Liu S., A disposable immunosensor device for point-of-care test of tumor marker based on copper-mediated amplification, Biosens. Bioelectron., 2013, 43, 425-431.10.1016/j.bios.2012.12.047Search in Google Scholar PubMed

[29] Becke A.D., Density-functional thermochemistry. III. The role of exact exchange, J. Chem. Phys., 1993, 98, 5648-5652.10.1063/1.464913Search in Google Scholar

[30] Frisch M.J., Trucks G.W., Schlegel H.B., Scuseria G.E., Robb M.A., Cheeseman J.R., et al., Gaussian 09, revision A.01, Inc., Wallingford CT, 2009.Search in Google Scholar

[31] Lee C., Yang W., Parr R.G., Development of the colle-salvetti correlation-energy formula into a functional of the electron density, Phys. Rev. B, 1988, 37, 785-789.10.1103/PhysRevB.37.785Search in Google Scholar

[32] Pietro W.J., Hehre W.J., Molecular orbital theory of the properties of inorganic and organometallic compounds. 3. STO-3G basis sets for first- and second-row transition metals, J. Comput. Chem., 1983, 4, 241-251.10.1002/jcc.540040215Search in Google Scholar

[33] Hehre W.J., Ditchfield R., Pople J.A., Self-consistent molecular orbital methods. XII. Further extensions of gaussian-type basis sets for use in molecular orbital studies of organic molecules, J. Chem. Phys., 1972, 56, 2257-2261.10.1063/1.1677527Search in Google Scholar

Received: 2016-4-23
Accepted: 2016-5-24
Published Online: 2016-12-30
Published in Print: 2016-1-1

© 2016 Y. Hao et al., published by De Gruyter Open

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License.

Downloaded on 30.1.2023 from
Scroll Up Arrow