LETTER pubs.acs.org/NanoLett
Pressure-Mediated Doping in Graphene Jimmy Nicolle, Denis Machon, Philippe Poncharal, Olivier Pierre-Louis, and Alfonso San-Miguel* Universite de Lyon, F-69000, France; Universite Lyon 1, Laboratoire PMCN; CNRS, UMR 5586; F-69622 Villeurbanne Cedex
bS Supporting Information ABSTRACT: Exfoliated graphene and few layer graphene samples supported on SiO2 have been studied by Raman spectroscopy at high pressure. For samples immersed on a alcohol mixture, an electron transfer of ∂n/∂P ∼ 8 1012 cm 2 GPa 1 is observed for monolayer and bilayer graphene, leading to giant doping values of n ∼ 6 1013 cm 2 at the maximum pressure of 7 GPa. Three independent and consistent proofs of the doping process are obtained from (i) the evolution of the Raman G-band to 2D-band intensity ratio, (ii) the pressure coefficient of the G-band frequency, and (iii) the 2D band components splitting in the case of the bilayer sample. The charge transfer phenomena is absent for trilayer samples and for samples immersed in argon or nitrogen. We also show that a phase transition from a 2D biaxial strain response, resulting from the substrate drag upon volume reduction, to a 3D hydrostatic compression takes place when going from the bilayer to the trilayer sample. By model calculations we relate this transition to the unbinding of the graphene SiO2 system when increasing the number of graphene layers and as function of the surface roughness parameters. We propose that the formation of silanol groups on the SiO2 substrate allows for a capacitance-induced substrate-mediated charge transfer. KEYWORDS: Carbon, graphene, Raman spectroscopy, strain sensors, high pressure, piezoelectric actuator raphene1 4 raises an increasing interest due to its unusual electronic properties,2 6 for its extraordinary membranetype mechanical response,4,7 12 or its electromechanical response.13 Pressure application continuously modifies the interaction of low-dimensional systems with its environment,14 which in the case of supported graphene under hydrostatic pressure15 are the pressure transmitting medium (PTM) and the substrate on which graphene lays. Because of the weakness of the interlayer van der Waals interactions, the in-plane elastic properties of graphene have been proposed to be directly transposed from the graphite ones as confirmed by a number of experiments and calculations.16 18 Moreover, the interaction of graphene with a rough surface probes also the out-of-plane elastic behavior. It has been shown that graphene can conform to high rough surfaces down to subnanometric features with nearly 99% fidelity.19 This situation allows for very strong substrate graphene interactions that can be additionally increased by pressure through a double mechanism: substrate graphene distance reduction and substrate shrinkage hauling of the carbon system.15 With an increasing number of layers, the carbon-based membrane is expected to rigidify and a modification of this scenario seems plausible. Exploring the high-pressure behavior of different number of graphene layers supported by SiO2-covered substrates appears then as an appealing method to probe both the graphene substrate interactions and the elastic properties of the graphene system. Using Raman spectroscopy, a powerful tool for graphene and carbon nanotube studies,20 we have studied the pressure-induced evolution of exfoliated samples with a number of graphene layers,
G
r 2011 American Chemical Society
n = 1, 2, and 3 supported on silica-coated silicon substrates. Samples were loaded in a diamond anvil cell (Figure 1a) up to a maximum pressure of 7.0 GPa using three different pressure transmitting media (PTM): a standard mixture of (4:1) methanol/ ethanol, nitrogen, and argon. In Figure 2 is shown the pressure evolution of the graphene and bilayer graphene Raman spectra under pressure in the alcohol mixture PTM. Spectra obtained on trilayer samples or for the other PTM were of similar quality. We appreciate at ambient pressure the G and 2D bands of these two samples. The 2D band allows the identification of the number of layers21,22 (Figure 1b). These Raman bands can be distinctly followed under pressure, but the additional signal from the diamond anvils or from the PTM need to be taken into account in the fitting protocol as shown in Figure 2. The distinct pressure evolution of these two bands as function of the number of graphene layers or the PTM will constitute the key experimental information. With increasing pressure, we observe in all cases an expected blue shift of the G band which is linear with the applied pressure (Figure 3a). The pressure-induced evolution is reversible except for a partial peak broadening. Values of the linear coefficients, ∂ωG/∂P, are plotted (Figure 3b) as a function of the number of layers. The most striking result is that the monolayer and bilayer samples show similar values of ∂ωG/∂P and that a sudden drop takes place at n = 3, which adopts a value close to graphite. Received: April 14, 2011 Revised: July 20, 2011 Published: August 01, 2011 3564
dx.doi.org/10.1021/nl201243c | Nano Lett. 2011, 11, 3564–3568
Nano Letters
LETTER
Figure 1. Experimental setup. (a) A diamond anvil cell (DAC) is the pressure generator. The SiO2/Si substrate (b) is placed in between the two diamond anvils, confined by a metallic gasket, and immersed in pressure transmitting media. The graphene and few-layer graphene samples on the substrate (c) can be identified and followed under pressure by Raman spectroscopy. (d) The corresponding Raman signal is shown, performed using a red laser (λ = 647.1 nm) as excitation light, for the n = 1, 2, and 3 layered samples seen in (c).
To understand the origin of such differences, let us consider the mechanical response of graphene. The pressure induced biaxial hydrostatic strain (ε) in the plane, which will modify the G-band frequencies, can be related to the applied stress (σ) by the compliance tensor Sij, as ε2D = 2(S11+S12)σ if we consider a pure biaxial compression or ε3D = 2(S11 + S12 + S13)σ if we consider the effect in the biaxial compression of the c-axis interaction. Furthermore, it is possible to calculate the G band shift as a function of applied biaxial stress using the relation23 ωG ω0 = Aε/ω0 with ω0 and ωG the position of the G band, respectively, at ambient pressure and at a given pressure and the A parameter derived from the data obtained by Blakslee et al.24 We can then calculate the G band pressure derivative in case of pure biaxial compression (ε = ε2D). We obtain, ∂ωG/∂P = 7.4 cm 1 3 GPa 1 in close agreement with our obtained values for n = 1, 2 in argon and nitrogen PTM. We then apply this model for graphite by placing ε = ε3D. The calculated pressure slope is then ∂ωG/∂P = 4.4 cm 1 3 GPa 1, which is close to graphite values in literature and to our values measured for graphite and trilayer graphene samples (see Figure 3b). The observed transition at n = 2f3 is then well explained as a transition in the mechanical response of graphene from biaxial compression to hydrostatic compression. We note here that for graphene (n = 1) in argon PTM we find (∂ω2D/∂P)/(∂ωG/∂P) = 2.8, that is, exactly the same value that can be obtained from the data of Ding et al.10 for biaxially strained graphene using piezoelectric actuators. In order to investigate more precisely the observed change in graphene mechanical response with the number of layers, we can consider the adhesion potential of a membrane (the graphene plane) on a rough surface that can be quantified by comparing surface roughness, adhesion energy, and bending rigidity modulus of the membrane. The morphology of supported graphene on a rough substrate results from a competition between adhesion
forces, which are minimized when the graphene shape perfectly follows that of the substrate, and the elastic cost of graphene shape changes, which are minimized when graphene is perfectly flat. This competition, which has been addressed in detail in ref 25, results in a complex sequence of partial unbinding transitions. Following ref 25, the unbinding transition depends on a single dimensionless parameter (R), which compares the typical curvature of our substrate with the adhesion equilibrium curvatures of the graphene layers. We consider that the main contribution to the elastic energy of graphene is bending, and we neglect the contributions related to in-plane strain. The unbinding transition then depends on a single dimensionless parameter Rn = (keq/kg)1/2, where kg is the typical substrate curvature, and keq = (2γn/Cn)1/2 is the adhesion equilibrium curvature, with γn and Cn, respectively, the multilayer graphene adhesion energy on SiO2 and the bending rigidity. Considering typical values given for the accurate determination of the roughness of an SiO2 surface,19 perfect adhesion is expected for R > 0.8 0.86, and partial unbinding starts for R < 0.8 0.86 (0.8 for sharp needlelike patterns, and 0.86 for smooth sinusoidal patterns). We plot the evolution of (R) parameter versus the number of graphene layers (see Figure 4), and we obtain that the transition between perfect adhesion everywhere along the surface and partial or total unbinding appears to take place precisely between the bilayer and the trilayer sample. As a consequence, the carbon-layered system goes from the biaxial compression through the substrate for n e 2 (surface adhesion) to the hydrostatic compression by the PTM for n g 3 (insufficient adhesion). This transition represented graphically in Figure 2b and Figure 4 could then be dependent on the type of surface and of its roughness characteristics (for more information about calculation see Supporting Informations). Let us turn now to the particular values of ∂ωG/∂P in the alcohol experiment (Figure 3b). Strong n-doping effects of systems 3565
dx.doi.org/10.1021/nl201243c |Nano Lett. 2011, 11, 3564–3568
Nano Letters
LETTER
Figure 2. Evolution with pressure of the Raman signal of (a) monolayer and (b) bilayer samples immersed in a mixture of (4:1) methanol/ethanol using a red laser (λ = 647.1 nm) as excitation light. The blue areas show the G band (ωG = 1580 cm 1 at ambient pressure) and the 2D band (ω2D = 2680 cm 1 at ambient pressure), and the red areas show the additional signal due either to the PTM (i) or to the second order Raman active band of diamond (ii) diamond (not visible) due to the sample environment. These bands have also been taking into account on fitting as presented under each spectra with the blue lines corresponding to the sample Raman signal. The fitting process is shown for both G and the 2D graphene bands, using for the 2D-band one component for the single layer sample and four components for the bilayer sample. All sample peaks shift to higher energies with pressure application.
having alcohol radicals interacting with sp2 carbon materials have been already reported.26 The relation between doping and both the G-band position and the ratio of intensities of the 2D- and G-bands have been well characterized for n = 1 and n = 2 samples.27 29 We plot the relative pressure evolution of the I(2D)/I(G) ratio for n = 1, 2 in Figure 5a both for argon and alcohol PTM. For n = 2 for which 4 peaks are needed to describe the 2D band, we take the full integrated intensity.28 Argon is not expected to induce any doping and in fact I(2D)/I(G) does not evolve significantly. We have used the measured values of I(2D)/I(G) and their pressure evolution to quantify both the initial doping values at the pressure cell loading (some kbars) and to calculate the doping rate with pressure.27 The origin of the initial doping could be associated to the samples alcohol cleaning protocol. At these doping levels, following Das et al.,27,28 I(2D)/I(G) evolves
linearly with pressure both with increasing p- or n-type doping. We obtain much more consistent results considering n-type doping, which is in good correspondence with other observations.26 We find an electron initial doping at loading conditions of ∼7 1012 cm 2 both for n = 1, 2 independently of the used PTM. We obtain a doping rate with pressure of ∂n/∂P ∼ 8 1012 electrons cm 2 3 GPa 1, both for the mono- and bilayer samples in alcohol. The value of ∂n∂P can be neglected within error bars for samples immersed in argon or nitrogen. In the doping domain here observed, the G-band frequency linearly increases with charge injection,27,28 and we can then evaluate the associated contribution of doping to ∂ωG/∂P. As it is shown in Table 1, this contribution is of ∼3.5 cm 1 3 GPa 1 and corresponds perfectly to the observed difference between ∂ωG/∂P in alcohol and argon PTM (Figure 3b), providing a cross-check of the pressure doping effect. 3566
dx.doi.org/10.1021/nl201243c |Nano Lett. 2011, 11, 3564–3568
Nano Letters
LETTER
Figure 3. Evolution of the Raman G-band position with pressure for different number of graphene layers(n) and different pressure transmitting media. (a) G-band frequency for n = 1, 2, and 3 immersed in alcohol. The solid line is the measured evolution under pressure for graphite, and the dotted lines are the linear fit of graphene data. The insert shows the evolution of the monolayer G band loaded in alcohol mixture during a pressure cycle. (b) Summary of the G band pressure linear coefficients obtained in alcohol (red squares), argon (blue circles), and nitrogen (green triangle) depending on the number of layers (n). Two effects are illustrated in the figure and discussed in the article: (i) a change of behavior between n = 2 and n = 3 due to a modification of the sample substrate binding as illustrated by the inserts and (ii) an electron doping process observed in the n = 1, 2 samples immersed in alcohol, which is highlighted by the red arrows.
Figure 4. Dimensionless parameter R as a function of n. The red (empty) and the blue (filled symbols) symbols respectively represent the results with kg = 0.3 nm 1, and kg = 0.4 nm 1. Triangles and circles correspond19 to γ1 = 0.6 eV 3 nm 2, and γ2 = 0.5 eV 3 nm 2. Unbinding is expected to occur in the vicinity of the shaded area between R = 0.8 and R = 0.86.
For the n = 2 case, the doping effect can even be confirmed when considering the relative pressure evolution of the four 2Dband Raman components (Figures 2 and 5b). There is no additional pressure-induced splitting of these modes for argon PTM, while there is a progressive splitting in alcohol PTM, which has been predicted as having its origin on doping.30 For the n = 3 sample, a detail analysis of the small differences between the measured values of ∂ωG/∂P in the different PTM and their proximity to the graphite ones, allow us to bound the doping effect to values at least two times smaller than for n = 1, 2. The mechanical transition between n = 2 and 3 from a biaxial to a hydrostatic response associated to graphene substrate unbinding is then concomitant with a reduction of charge injection. All these observations allow us to propose that the electron charge injection observed in alcohol-immersed samples takes place through the SiO2 substrate as it has been observed for O2 or H2O.31 This is reinforced by the fact that both the initial doping
Figure 5. Doping effect on the Raman signal. (a) Relative evolution with pressure of the ratio of the intensities of the 2D and G Raman peaks. The closed symbols are used for the monolayer samples and the opened symbols for the bilayer samples, loaded in mixture of alcohol (black squares) and in argon (red triangles). These ratios have been quantified as doping dependent.28 (b) Pressure evolution of the bilayer 2D-band components in alcohol (b1) and in argon (b2). The different colors correspond to the four different components of the bilayer 2D band.21 The observed splitting in alcohol has been predicted as originated by doping.30
values of the n = 1,2 samples and their pressure derivatives are the same. This is consistent with a one-layer charge transfer process. 3567
dx.doi.org/10.1021/nl201243c |Nano Lett. 2011, 11, 3564–3568
Nano Letters
LETTER
Table 1. Pressure Induced Doping Dn/DP Obtained from the Evolution of the G-Band to 2D-Band Intensity Ratio and Its Contribution to the G-Band Position under Pressure ([DωG/DP]dop) for Graphene (n = 1) and Bilayer Samples (n = 2)a ∂n/∂P (1013) n PTM (cm 1 2
2
∂ωG/∂P [∂ωG/∂P]dop
GPa 1) (cm
1
(measured)
GPa 1) (cm
1
[∂ωG/∂P]mech
GPa 1) (cm
1
GPa 1)
alc
0.7 ( 0.2
3.6 ( 1.1
10.5 ( 0.2
6.9 ( 1.4
Ar
0.2 ( 0.2
1.0 ( 1.1
7.6 ( 1.0
6.6 ( 2.0
alc
0.8 ( 0.2
3.4 ( 1.1
10.4 ( 0.3
7.0 ( 1.4
Ar
0.1 ( 0.2
0.3 ( 0.6
6.9 ( 1.0
6.6 ( 1.6
N2
0.1 ( 0.2
0.3 ( 0.6
6.9 ( 0.2
6.6 ( 0.8
a
The mechanical contribution to the G-band pressure coefficient (last column) is then obtained by subtracting the doping contribution from the measured one. Within error bars all mechanical contributions coincide with the calculated value for biaxial compression of a graphene system (7.4 cm-1 3 GPa-1, see text) and with independence of the used pressure transmitting media (PTM): argon, nitrogen and a 4:1 ethanol/ methanol mixture (alc).
The formation of silanol groups (Si O H) on the surface of oxidized silicon substrates, which has been proposed as a potential source of charge injected through the surface in carbon nanotubes,32 could then be here at the origin of the observed doping in graphene. The pressure linear evolution of such charge injection could be associated to a capacitive mechanism, which is in correspondence with the expected approximate linear evolution of the van der Waals driven evolution of the graphene SiO2 distance.33 At our maximum pressure of 7 GPa, the observed doping corresponds to a Fermi level shift close to 1 eV (1.4 eV in the Dirac cone approximation).27,34 On the basis of the observed mechanism, it is possible to conjecture about the development of strain sensors based on contacted graphene immersed in selected fluids.
’ ASSOCIATED CONTENT
bS
Supporting Information. Description of experimental methods, samples preparation, high pressure apparatus, Raman spectra under high pressure of the different samples, and the development of the theoretical model. This material is available free of charge via the Internet at http://pubs.acs.org.
’ AUTHOR INFORMATION Corresponding Author
*E-mail:
[email protected].
’ ACKNOWLEDGMENT We sincerely acknowledge G. Montagnac from LGL (ENS Lyon, CNRS and University Lyon 1) and R. Mascart from LPMCN (CNRS and University Lyon 1) for technical support. ’ REFERENCES
(3) Zhang, Y.; Tan, Y.; Stormer, H.; Kim, P. Nature 2005, 438, 201–204. (4) Meyer, J.; Geim, A.; Katsnelson, M.; Novoselov, K.; Booth, T.; Roth, S. Nature 2007, 446, 60–63. (5) Gomez-Navarro, C.; Weitz, R.; Bittner, A.; Scolari, M.; Mews, A.; Burghard, M.; Kern, K. Nano Lett. 2007, 7, 3499–3503. (6) Eda, G.; Fanchini, G.; Chhowalla, M. Nat. Nanotechnol. 2008, 3, 270–274. (7) Dikin, D.; Stankovich, S.; Zimney, E.; Piner, R.; Dommett, G.; Evmenenko, G.; Nguyen, S.; Ruoff, R. Nature 2007, 448, 457–460. (8) Lee, C.; Wei, X.; Kysar, J.; Hone, J. Science 2008, 321, 385–388. (9) Ferralis, N. J. Mater. Sci. 2010, 45, 5135–5149. (10) Ding, F.; Ji, H.; Chen, Y.; Herklotz, A.; D€orr, K.; Mei, Y.; Rastelli, A.; Schmidt, O. G. Nano Lett. 2010, 10, 3453–3458. (11) Huang, M.; Yan, H.; Heinz, T. F.; Hone, J. Nano Letters 2010, 10 (10), 4074–4079. (12) Metzger, C.; Remi, S.; Liu, M.; Kusminskiy, S.; Neto, A.; Swan, A.; Goldberg, B. Nano Lett. 2010, 10, 6–10. (13) Huang, M.; Pascal, T.; Kim, H.; Goddard, W.; Greer, J. Nano Lett. 2011, 11, 1241–1246. (14) San-Miguel, A. Chem. Soc. Rev. 2006, 35, 876–889. (15) Proctor, J.; Gregoryanz, E.; Novoselov, K.; Lotya, M.; Coleman, J.; Halsall, M. Phys. Rev. B 2009, 80, 073408. (16) Thomsen, C.; Reich, S.; Ordejon, P. Phys. Rev. B 2002, 65, 155411. (17) Ni, Z.; Wang, H.; Ma, Y.; Kasim, J.; Wu, Y.; Shen, Z. ACS Nano 2008, 2, 1033–1039. (18) Mohiuddin, T. M. G.; Lombardo, A.; Nair, R. R.; Bonetti, A.; Savini, G.; Jalil, R.; Bonini, N.; Basko, D. M.; Galiotis, C.; Marzari, N.; Novoselov, K.; Geim, A. K.; Ferrari, A. C. Phys. Rev. B 2009, 79, 205433. (19) Cullen, W.; Yamamoto, M.; Burson, K.; Chen, J.; Jang, C.; Li, L.; Fuhrer, M.; Williams, E. Phys. Rev. Lett. 2010, 105, 215504. (20) Dresselhaus, M.; Jorio, A.; Hofmann, M.; Dresselhaus, G.; Saito, R. Nano Lett. 2010, 10, 751–758. (21) Ferrari, A.; Meyer, J.; Scardaci, V.; Casiraghi, C.; Lazzeri, M.; Mauri, F.; Piscanec, S.; Jiang, D.; Novoselov, K.; Roth, S.; Geim, A. Phys. Rev. Lett. 2006, 97, 187401. (22) Poncharal, P.; Ayari, A.; Michel, T.; Sauvajol, J.-L. Phys. Rev. B 2009, 79, 195417. (23) Sakata, H.; Dresselhaus, G.; Dresselhaus, M.; Endo, M. J. Appl. Phys. 1988, 63, 2769–2772. (24) Blakslee, O. J. Appl. Phys. 1970, 41, 3373–&. (25) Pierre-Louis, O. Phys. Rev. E 2008, 78, 021603. (26) Blake, P.; Brimicombe, P. D.; Nair, R. R.; Booth, T. J.; Jiang, D.; Schedin, F.; Ponomarenko, L. A.; Morozov, S. V.; Gleeson, H. F.; Hill, E. W.; Geim, A. K.; Novoselov, K. S. Nano Lett. 2008, 8, 1704–1708. (27) Das, A.; Pisana, S.; Chakraborty, B.; Piscanec, S.; Saha, S. K.; Waghmare, U. V.; Novoselov, K. S.; Krishnamurthy, H. R.; Geim, A. K.; Ferrari, A. C.; Sood, A. K. Nat. Nanotechnol. 2008, 3, 210–215. (28) Das, A.; Chakraborty, B.; Piscanec, S.; Pisana, S.; Sood, A. K.; Ferrari, A. C. Phys. Rev. B 2009, 79, 155417. (29) Gava, P.; Lazzeri, M.; Saitta, A.; Mauri, F. Phys. Rev. B 2009, 80, 155422. (30) Attaccalite, C.; Wirtz, L.; Lazzeri, M.; Mauri, F.; Rubio, A. Nano Lett. 2010, 10, 1172–1176. (31) Ryu, S.; Liu, L.; Berciaud, S.; Yu, Y.; Liu, H.; Kim, P.; Flynn, G.; Brus, L. Nano Lett. 2010, 10, 4944–4951. (32) Lee, J.; Ryu, S.; Yoo, K.; Choi, I.; Yun, W.; Kim, J. J. Phys. Chem. C 2007, 111, 12504–12507. (33) Caillier, C.; Ayari, A.; Gouttenoire, V.; San Miguel, A.; Jourdain, V.; Picher, M.; Sauvajol, J. Appl. Phys. Lett. 2010, 97, 173111. (34) Giovannetti, G.; Khomyakov, P.; Brocks, G.; Karpan, V.; van den Brink, J.; Kelly, P. Phys. Rev. Lett. 2008, 101, 026803.
(1) Novoselov, K.; Geim, A.; Morozov, S.; Jiang, D.; Zhang, Y.; Dubonos, S.; Grigorieva, I.; Firsov, A. Science 2004, 306, 666–669. (2) Novoselov, K.; Jiang, D.; Schedin, F.; Booth, T.; Khotkevich, V.; Morozov, S.; Geim, A. Proc. Natl. Acad. Sci. U.S.A. 2005, 102, 10451–10453. 3568
dx.doi.org/10.1021/nl201243c |Nano Lett. 2011, 11, 3564–3568