This article is part of the series INEC 2011.

Open Access Nano Express

Density functional theory calculations on graphene/α-SiO2(0001) interface

Zhimin Ao12*, Man Jiang3, Zi Wen1* and Sean Li2

Author affiliations

1 Key Laboratory of Automobile Materials, Ministry of Education and Department of Materials Science and Engineering, Jilin University, Changchun, 130022, China

2 School of Materials Science and Engineering, The University of New South Wales, Sydney, NSW 2052, Australia

3 Department of Chemistry, Jilin University, Changchun, 130022, China

For all author emails, please log on.

Citation and License

Nanoscale Research Letters 2012, 7:158  doi:10.1186/1556-276X-7-158


The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/7/1/158


Received:23 November 2011
Accepted:28 February 2012
Published:28 February 2012

© 2012 Ao et al; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

In this work, the graphene/α-SiO2(0001) interface is calculated using density functional theory. On the oxygen-terminated SiO2 surface, atomic structure reconstruction occurs at the graphene/SiO2 interface to eliminate the dangling bonds. The interface interaction is 77 meV/C atom, which indicates that van der Waals force dominates the interaction, but it is stronger than the force between the graphene layers in graphite. The distance between graphene and the SiO2 surface is 2.805 Å, which is smaller than the 3.4 Å interlayer distance of graphite. In addition, the SiO2 substrate induces p-type doping in graphene and opens a small gap of 0.13 eV at the Dirac point of graphene, which is desirable for electronic device applications.

Introduction

Graphene, a single two-dimensional layer of graphite in hexagonal structure, is the starting point for many nanographite devices with promising electrical properties [1]. After the theoretical prediction of the peculiar electronic properties of graphene in 1947 by Wallace [2] and the subsequent studies on its magnetic spectrum [3,4], it took half a century until the graphene could be first experimentally fabricated [1], and its anomalous quantum Hall effect has been measured [5-7], which encourages numerous works on it now [8-11].

However, the instability of a freestanding graphene (it has an intrinsic three-dimensional structure or ripples) [12] leads graphene to be used on devices by laying it on a substrate. For example, the insulating substrates made from silicon dioxide [SiO2] are widely used as a dielectric medium in electronic devices. The electrical properties of graphene can be modified using electrical gates, substrates, and chemical species such as atoms and molecules [13,14]. It has been reported that graphene grown on a SiC surface was n-type and was exhibiting a gap of about 0.26 eV [15]. A theoretical study on bilayer graphene grown on a SiC surface showed that the first carbon layer formed covalent bonds with SiC and acted as a buffer layer, and the graphene nature properties were recovered by the second carbon layer [16,17]. For graphene supported by SiO2, previous works studied the electrical properties of graphene and the doping effect of the SiO2 substrate [18-23]. They found that graphene was p doped when the graphene was weakly bonded to an O-terminated surface with hydrogen passivation. However, graphene would exhibit a finite bandgap if there is a strong interaction with an O-terminated surface while n doping took place on a Si-terminated surface with the active dangling bonds [20,21]. In these works, they considered that graphene would form C-O or C-Si bonds on the interface due to the dangling bonds associated with the absence of hydrogen passivation. However, when graphene is used in an electronic device, the graphene layer is usually deposited on a surface of the SiO2 substrate where the atomic structure of the surface is usually considered to reconstruct to eliminate the dangling bonds [24]. In addition, the atomic force microscopy [AFM] image indicated that the height of a single graphene layer on the SiO2 substrate is around 4 Å, which is a bit larger than the 3.4 Å of the graphite interlayer distance [25]. Therefore, the interaction between graphene and the SiO2 surface should be a weak van der Waals interaction without the formation of strong covalent bonds.

In this work, we investigate the graphene/α-SiO2(0001) interface through density functional theory [DFT] calculations with considerations on surface reconstruction of SiO2. The interfacial atomic structure and the effects of the substrate on corresponding electronic structure of graphene will also be studied.

Simulation details

There are two fundamental types of surfaces: the surface with the termination of the Si atoms or the surface with the termination of the O atoms. Based on the atomic structure study of α-SiO2(0001) surface using the first principles [24], the O-terminated surface was more stable where covalent bonds between two O atoms on a particular Si atom were formed. Therefore, the O-terminated surface is selected for study in this work.

All DFT calculations in this work are implemented by the DMOL3 code, which utilizes a norm-conserving pseudopotential to perform the first principle quantum mechanical calculations [26]. Local density approximation with the PWC function is employed as exchange correlation function [27]. The spin restriction is taken, which would not lead to a big error since the dangling oxygen bonds would be eliminated through the formation of covalent bonds between the two O atoms on the particular Si atom as shown in Figure 1b[24]. The K points are 6 × 6 × 6 for the bulk and 6 × 6 × 1 for all slabs, which make the convergence tolerance of energy of 1.0 × 10-5 eV/atom, maximal force of 0.03 eV/Å, and displacement of 1.0 × 10-3 Å.

thumbnailFigure 1. Atomic stacking sequence of graphene/α-SiO2(0001) interface without relaxation (a) and in structure A (b). The yellow, gray, and red spheres show Si, C, and O atoms. The numbers are the index of atoms in the simulation. The meanings of parameters in the figure are given in the text.

The graphene/α-SiO2(0001) interface is constructed based on the bulk structures. A vacuum thickness of 12 Å is added above the graphene layer to ensure the interaction between the repeated slabs in the normal direction of the surface to be small enough. The simulation cell parameters of the graphene and the SiO2 substrate are a = b = 4.920 Å, α = β = 90°, and γ = 120° and a = b = 4.910 Å, α = β = 90°, and γ = 120°, respectively, where a small lattice mismatch is within a range of 0.3%, as shown in Figure 1a. Based on the AFM result, it is known that the height of the single graphene layer on the SiO2 substrate is around 4 Å [25]. Thus, in the initial structure of graphene/α-SiO2 interface, as shown in Figure 1a, the distance between graphene and SiO2 [d0] is set as 4.0 Å. In the simulation process, a structural relaxation process is allowed for C atoms in the graphene and all Si and O atoms in the substrate, except that the atoms of the lowest O-Si-O monolayer are fixed.

Results and discussion

The atomic structure of the favorite graphene/α-SiO2(0001) interface configuration with the distance between the graphene and SiO2 surface, d0 = 2.820 Å, is shown in Figure 1b, which is named as structure A; d0 is in the van der Waals distance range, but it is smaller than the interlayer distance of graphite of 3.4 Å. From the figure, the surface reconstruction occurred after the relaxation which a covalent bond forms between the two O atoms on the particular Si atom to eliminate the dangling bonds. The related structure parameters are given in Table 1. As shown in the table, a bond length expansion between the surface or interface atoms is found with l0bulk < l0A < l0clean, where l0 denotes Si-O bond length, and the subscripts bulk, clean, and A denote bulk SiO2, clean SiO2 slab, and structure A, respectively. The result l0bulk = 1.596 Å is consistent with the reported result of 1.61 Å [20] while the extension of Si-O bonds at the surface or interface is found. We consider that the extension is caused by the one Si-O broken bond on the substrate surface or at the interface. In addition, the weak interaction between the graphene and the SiO2 surface slightly reduces the l0 at the interface compared with that on the surface.

Table 1. Atomic structure parameters of bulk SiO2, clean SiO2(0001) surface, and graphene/SiO2(0001) interface

In addition, Table 1 also shows the distance between the first and second Si layers [d1] and the Si-O-Si bond angle [α1] in different structures; d1bulk > d1clean > d1A and α1bulk > α1clean > α1A are found. In light of Figure 1, d1 is affected by both l0 and α1. The former brings out d1bulk < d1clean, and the latter leads to d1bulk > d1clean. Therefore, the effect of α1 on d1 is stronger than that of l0. Comparing structure A with the clean SiO2 surface slab, both l0 and α1 lead to the decrease of d1, where d1 and l0 of structure A shrink to about 4%.

Quantitative experimental data for the interaction strength of graphene/substrate interface is very limited. The first principle calculations showed that the binding energy of graphene on a Si-terminated SiO2 surface is around 20 meV/C atom with interface distance d0 = 3.29 Å and that on a hydrogen passivation O-terminated SiO2 surface is 0.13 eV/C atom with d0 = 2.58 Å [20]. In addition, the interlayer binding energy in graphite was reported to be 50 to 60 meV [28]. In the system of structure A, the binding energy calculated between graphene and the SiO2 substrate is about 77 meV/C atom, which is larger than that in graphite and also in the graphene on Si-terminated SiO2 surface, but it is smaller than that on the hydrogen-passivated O-terminated SiO2 surface. Note that d0 in structure A is 2.820 Å, which is smaller than 3.14 Å in graphite and graphene laid on the Si-terminated SiO2 surface, but it is larger than that on hydrogen-passivated O-terminated SiO2 surface. It is known that the binding energy is inversely proportional to d0. The system of graphene on hydrogen-passivated O-terminated SiO2 surface has the smallest d0 and the strongest binding energy. A previous study implied that the C-C bond is weakened through the strengthening of bonds to the substrate [29]. A similar phenomenon is found in structure A where the C-C bond length of graphene is 1.430 Å, which is longer than the 1.420 Å in graphite. Thus, the adsorption of graphene on SiO2 with structure A is stronger than that between the bulk graphite layers as shown previously.

On the other hand, in order to understand the effect of substrate on the graphene electrical properties, the band structure for three different systems in Figure 1b, namely free graphene monolayer, clean SiO2(0001) slab, and graphene/SiO2 interface, are represented in Figure 2. The band structure of free graphene monolayer is calculated on a periodic structure where the graphene monolayer is separated by the same vacuum distance as that for the graphene/SiO2 slabs. This band structure shows the crossing of π and π* bands at the K point and also at the Fermi level, which agrees with a well-known result that graphene is semimetallic with a 0 bandgap. The band structure of the interface shows that the bands of graphene layer are open with a 0.13-eV gap at the K point. It is interesting to note that the Fermi level is lowered with the amount of transferred charges. Thus, the charge should be transferred from graphene to the substrate in structure A. As a consequence of the charge transfer, the SiO2 substrate induces p-type doping in the graphene. Except for these variations, the band structure of interface is almost identical with the sum of band structure of the free graphene monolayer and SiO2 slab, as shown in Figure 2a, b.

thumbnailFigure 2. Calculated band structures for graphene monolayer (a), clean α-SiO2(0001) slab, (b) and graphene/α-SiO2 interface (c). The dash line at 0 value denotes the Fermi level.

The charge transfer and atomic charge can be obtained using the Mulliken analysis; it is shown in Table 2. Mulliken analysis is performed using a projection of the plane wave states onto a localized basis with the technique described by Sanchez-Portal et al. [30]. Subsequently, the resulting projected states are performed using the Mulliken formalism [31]. This technique has been widely used to analyze the electronic structures performed with linear combinations of atomic orbital basis sets. As shown in Table 2, all C atoms of graphene are positively charged in structure A, and the graphene transfers 0.144 e to the SiO2 substrate where e denotes one electron charge. This is consistent with the result of the band structure where the p-type-doped graphene induced by the substrate is found.

Table 2. Charge transfer and atomic charge obtained by Mulliken analysis

In addition, the electron charges of O4 and O8 atoms at the SiO2 surface are -0.448 and -0.424 e, respectively, for structure A. Those in the clean SiO2(0001) slab are -0.396 and -0.293 e, and in bulk SiO2 are both -0.781 e. On the other hand, the charges of Si15, which binds with O4 and O8, are 1.681, 1.510, and 1.562 e in structure A, clean SiO2 slab, and bulk SiO2, respectively. It is known that the bond strength is in proportion with the multiplication absolute value of changes of the O and Si atoms. Therefore, the interaction between O4 and Si15, and O8 and Si15 are the strongest in bulk SiO2, followed by that in structure A, and the weakest was that in the SiO2 slab. These also agree with the results in Table 1 where the Si-O l0 at the surface or interface is l0bulk < l0A < l0clean.

Conclusions

In conclusion, DFT calculations were employed to study the graphene/α-SiO2(0001) interface with an oxygen-terminated SiO2 surface. After geometry relaxation, a stable structure with a distance between graphene and the SiO2 surface, d0 = 2.805 Å, was observed. In addition, the structure reconstruction of the SiO2 surface took place where the two O atoms binding on a particular Si atom form a covalent bond on the surface to eliminate the surface dangling bonds. Simulation results indicate that the interface interaction is 77 meV/C atom, which indicates that the van der Waals force dominates the interaction. Furthermore, the SiO2 substrate induces p-type doping in graphene and open a small gap of 0.13 eV at the Dirac point of graphene.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

ZA did the calculation and drafted the manuscript. ZW, MJ, and SL co-drafted the manuscript. All authors read and approved the final manuscript.

Acknowledgements

The financial supports from the Vice-Chancellor's Postdoctoral Research Fellowship Program of the University of New South Wales (SIR50/PS19184), the ECR grant of the University of New South Wales (SIR30/PS24201), the Faculty Grant of the University of New South Wales (IR001/PS27218), ARC Discovery project of DP1096769 of Australia, and the National Key Basic Research and Development Program (grant no. 2010CB631001) of China are acknowledged.

References

  1. Novoselov KS, Geim AK, Morozov SV, Zhang D, Jiang Y, Dubonos SV, Grigorieva IV, Firsov AA: Electric filed effect in atomically thin carbon films.

    Science 2004, 306:666-669. PubMed Abstract | Publisher Full Text OpenURL

  2. Wallace PR: The band theory of graphite.

    Phys Rev 1947, 71:622-634. Publisher Full Text OpenURL

  3. McClure JW: Diamagnetism of graphite.

    Phys Rev 1956, 104:666-671. Publisher Full Text OpenURL

  4. Zheng Y, Ando T: Hall conductivity of a two-dimensional graphite system.

    Phys Rev B 2002, 65:245420. OpenURL

  5. Zhang Y, Tan Y-W, Strormer HL, Kim P: Experimental observation of the quantum Hall effect and Berry's phase in graphene.

    Nature 2005, 438:201-204. PubMed Abstract | Publisher Full Text OpenURL

  6. Novoselov KS, Geim AK, Morozov SV, Jiang D, Katsnelson MI, Grigorieva IV, Dubonos SV, Firsov AA: Two-dimensional gas of massless Dirac fermions in graphene.

    Nature 2005, 438:197-200. PubMed Abstract | Publisher Full Text OpenURL

  7. Zhang Y, Jiang Z, Small JP, Purewal MS, Tan Y-W, Fazlollahi M, Chudow JD, Jaszczak JA, Stormer HL, Kim P: Landau-level splitting in graphene in high magnetic fields.

    Phys Rev Lett 2006, 96:136806. PubMed Abstract | Publisher Full Text OpenURL

  8. Gusynin VP, Sharapov SG: Unconventional integer quantum Hall effect in graphene.

    Phys Rev Lett 2005, 95:146801. PubMed Abstract | Publisher Full Text OpenURL

  9. Kane CL, Mele EJ: Quantum spin Hall effect in graphene.

    Phys Rev Lett 2005, 95:226801. PubMed Abstract | Publisher Full Text OpenURL

  10. Guinea FA, Castro Neto H, Peres NMR: Electronic states and Laudau levels in graphene stacks.

    Phys Rev B 2006, 73:245426. OpenURL

  11. Hasegawa Y, Kohmoto M: Quantum Hall effect and the topological number in graphene.

    Phys Rev B 2006, 74:155415. OpenURL

  12. Meyer JC, Geim AK, Katsnelson MI, Novoselov KS, Booth TJ, Roth S: The structure of suspended graphene sheets.

    Nature 2007, 446:60. PubMed Abstract | Publisher Full Text OpenURL

  13. Ao ZM, Yang J, Li S, Jiang Q: Enhancement of CO detection in Al doped graphene.

    Chem Phys Lett 2007, 461:276-279. OpenURL

  14. Leenaerts O, Partoens B, Peeters FM: Adsorption of H2O, NH3, CO, NO2, and NO on graphene: a first-principles study.

    Phys Rev B 2008, 77:125416. OpenURL

  15. Zhou SY, Gweon G-H, Fedorov AV, First PN, De Heer WA, Lee D-H, Guinea F, Castro Neto AH, Lanzara A: Substrate-induced bandgap opening in epitaxial graphene.

    Nat Mater 2007, 6:770. PubMed Abstract | Publisher Full Text OpenURL

  16. Varchon F, Feng R, Hass J, Li X, Nguyen BN, Naud C, Mallet P, Veuillen J-Y, Berger C, Conrad EH, Magaud L: Electronic structure of epitaxial graphene layer on SiC: effect of the substrate.

    Phys Rev Lett 2007, 99:126805. PubMed Abstract | Publisher Full Text OpenURL

  17. Mattausch A, Pankratov O: Ab initio study of graphene on SiC.

    Phys Rev Lett 2007, 99:076802. PubMed Abstract | Publisher Full Text OpenURL

  18. Shi Y, Dong X, Chen P, Wang J, Li L-J: Effective doping of single-layer graphene from underlying SiO2 substrates.

    Phys Rev B 2009, 79:115402. OpenURL

  19. Dharma-wardana MWC: Electronic structure and the minimum conductance of a graphene layer on SiO2 from density functional methods.

    J Phys Condens Matter 2007, 19:386228. Publisher Full Text OpenURL

  20. Kang Y-J, Kang J, Chang KJ: Electronic structure of graphene and doping effect on SiO2.

    Phys Rev B 2008, 78:115404. OpenURL

  21. Shemella P, Nayak SK: Electronic structure and band-gap modulation of graphene via substrate surface chemistry.

    Appl Phys Lett 2009, 94:032101. Publisher Full Text OpenURL

  22. Ao ZM, Zheng WT, Jiang Q: The effects of electronic field on the atomic structure of the graphene/α-SiO2 interface.

    Nanotechnology 2008, 19:275710. PubMed Abstract | Publisher Full Text OpenURL

  23. Liu Y, Ao ZM, Wang T, Wang WB, Yu B: Transformation from AA to AB-stacked bilayer graphene on α-SiO2 under electric field.

    Chin Phys Lett 2011, 28:087303. Publisher Full Text OpenURL

  24. Rignanese G-M, De Vita A, Charlier J-C, Gonze X, Car R: First-principles molecular-dynamics study of the (0001) α-quartz surface.

    Phys Rev B 2000, 61:13250-13255. Publisher Full Text OpenURL

  25. Ishigami MJ, Chen H, Cullen WG, Fuhrer MS, Williams ED: Atomic structure of graphene on SiO2.

    Nano Lett 2007, 7:1643-1648. PubMed Abstract | Publisher Full Text OpenURL

  26. Hamann DR, Schluter M: Norm-conerving pseudopotentials.

    Phys Rev Lett 1979, 43:1494-1497. Publisher Full Text OpenURL

  27. Perdew JP, Wang Y: Accurate and simple analytic representation of the electron-gas correlation energy.

    Phys Rev B 1992, 45:13244-13249. Publisher Full Text OpenURL

  28. Hasegawa M, Nishidate K: Semiempirical approach to the energetic of interlayer binding in graphite.

    Phys Rev B 2004, 70:205431. OpenURL

  29. Aizawa T, Souda R, Otani S, Ishizawa Y, Oshima C: Anomalous bond of monolayer graphite on transition-metal carbide surfaces.

    Phys Rev Lett 1990, 64:768-771. PubMed Abstract | Publisher Full Text OpenURL

  30. Sanchez-Portal D, Artacho E, Soler JM: Projection of plane-wave calculations into atomic orbitals.

    Solid State Commun 1995, 95:685-689. Publisher Full Text OpenURL

  31. Mulliken RS: Electronic population analysis on LCAO-MO molecular wave functions. I.

    J Chem Phys 1955, 23:1833. Publisher Full Text OpenURL