CIENCIAS NUCLEARES
Raman D-band in the irradiated graphene: Origin of the non-monotonous dependence of its intensity with defect concentration
La banda D de Raman del grafeno irradiado: Origen de la dependencia no monótona de su intensidad con la concentración de defectos
Daniel Codorniu Pujals
Instituto Superior de Tecnolog ías y Ciencias Aplicadas (INSTEC)
Avenida Salvador Allende esquinaa Luaces, Plaza de la Revolución, La Habana, Cuba
ABSTRACT
Raman spectroscopy is one of the experimental techniques more used in studying irradiated carbon nanostructures, in particular graphene, due to its high sensibility to the presence of defects in the crystalline lattice. Special attention has received the variation of the intensity of the Raman D-band of graphene with the concentration of defects produced by irradiation. Nowadays, there are enough experimental evidences about the non-monotonous character of that dependence, but the explanation of this behavior is still controversial. In the present work we developed a simplified mathematical model to obtain a functional relationship between these two magnitudes and showed that the non-monotonous dependence is intrinsic to the nature of the D-band and that it is not necessarily linked to amorphization processes. The obtained functional dependence was used to fit experimental data taken from other authors. The determination coefficient of the fitting was .
Key words: graphene, irradiation, Raman spectroscopy, nanostructures, mathematical models.
RESUMEN
La espectroscopía Raman es una de las técnicas experimentales más usadas en el estudio de las nano-estructuras de carbono irradiadas, en particular del grafeno, debido a su alta sensibilidad a la presencia de defectos en la red cristalina. Una atención especial ha recibido la variación de la intensidad de la banda D de los espectros Raman del grafeno con la concentración de defectos producida por la irradiación. Hoy día hay suficientes evidencias experimentales sobre el carácter no monótono de esa dependencia, pero la explicación de ese comportamiento todavía es polémica. En el presente trabajo se desarrolló un modelo matemático simplificado para obtener una relación funcional entre estas dos magnitudes y se demostró que la dependencia no-monótona es intrínseca a la naturaleza de la banda D y que no está ligada necesariamente a procesos de amorfización. La dependencia funcional obtenida fue utilizada para ajustar datos experimentales obtenidos por otros autores. Se obtuvo un ajuste con un coeficiente de determinación.
Palabras claves: grafeno, irradiación, espectroscopía Raman, nanoestructuras, modelos matemáticos.
INTRODUCTION
Graphene is a very interesting and promising material. Since its obtainment in 2004, a lot of attention has been devoted to the study of the exceptional electronic and physical properties of this two-dimensional nanostructure [1,2] and to the analysis of potential practical applications that cover many different fields from electronics to environmental protection and health [3]. Particularly interesting are the possible applications in biosciences , where there are a lot of ongoing research projects devoted to the use of graphene as the basis for develop DNA sequencing devices, biosensors, etc.[4].In many of the mentioned applications it is necessary or desirable to control the electronic structure and, in particular, the band gap. For this purpose it is possible to follow several approaches. One of them is the chemical functionalization of graphene, but another very important way is the controlled introduction of defects through its irradiation with different kind of particles. This alternative is being explored experimentally, as well as theoretically, by several groups around the world [5-8]. In the efforts devoted to understand the behavior of graphene under particle bombardment, an special role is played by Raman spectroscopy, a powerful technique that have been applied for several decades to the characterize of carbon bulk structures and nanostructures. The special usefulness of Raman spectroscopy in the study of irradiated graphene, graphite and carbon nanotubes is related with the so-called D-band, that only exists in the Raman spectra of these materials in case lattice defects or other kind of disorder [9]. The sensibility of this technique to the modifications to and defects in the lattice of carbon structures have been used to studied amorphization processes in graphite and other carbonaceous materials [10].An important experimental parameter in the study of the irradiated samples of graphite and carbon nanostructures is the ratio of the intensities of D-band and the G-band . The G-band is the main first order line in graphene spectra and it is not linked to the existence of defects [11]. Already in the seventies of the last century, it was clear that the ratio is strongly dependent of the concentration of defects in the crystal lattice. In 1970, Tuinstra and Konig [12] proposed for graphite the relationship
where
is the defect concentration1. Nevertheless, the experimental studies of irradiated graphene performed in the last years have shown that the Tuinstra-Konig relation is not applicable in those cases, and that the
ratio exhibits a non-monotonous behavior, i.e. for small concentration of defects the ratio increases with the concentration of defects, but after a critical concentration it begins to decrease [7,13,14]. The authors of those papers explained the decreasing of the ratio by the existence at high defect concentrations of very distorted zones in the lattice, where the material is no more crystalline but amorphous. They followed a theory developed in 2000 by Ferrari and Robertson [10] for the evolution of Raman spectra of graphite from pristine to amorphous samples, not directly related to irradiation processes. In that paper was reported by the first time the non-monotonous variation of
with the concentration.In spite of the fact that in many irradiation processes the sample gets some kind of amorphization, this is not always true, especially when the bombarding particles are electrons or protons with energies of several keV [15]. Besides that, in the zone of medium concentration of defects, when the ratio
begins to decrease with
the cited papers do not give enough evidence of the existence of amophization.Therefore, it is worth to obtain the theoretical dependence of the ratio
with
in order to elucidate whether the observed behavior is present there or it
requires additional mechanisms as ,for example, amorphization, for its
explanation. This was the main goal of the research reported in this paper.
Development of the model
In order to obtain a mathematical expression for the ratio, it is necessary
to consider the differences in the mechanisms that produce each band. As it
was mentioned above, the G-band is the only first order Raman process
allowed by the selection rule
=0 [9, 11]. On the other side, due to the
conical form of the energy bands in the graphene (Dirac cones) and to the
absence of gap, all the Raman lines in graphene are resonant. It is well
established in the literature [9,11] that the G band is a typical Raman resonant
line described by a Lorentzian function, as follows:
where is the frequency at the center of the G-band and
is the G-line
width.The nature of the D-band is totally different. In fact, it is linked with phonons
whose wavevectors are far from the center of the Brillouin zone and,
consequently, they are forbidden by the selection rule . Thomsen and
Reich [16] explained the existence of this band in the Raman spectra of
graphite by a double resonance process, in which the excited electron is
scattered first by a phonon and then by a lattice defect.The double resonance process is represented in the figure 1a. First , an
electron with wavevector k in the valence band (p), near of the Dirac point K, is
excited by the laser to the conduction band , i.e. an electron-hole pair is
created (step 1); then the excited electron is scattered by a phonon, acquiring a
wavevector k+q and passing to the Dirac cone near to the point
(step 2).
After that, the electron is scattered by a lattice defect, that absorbs the
wavevector q (step 3). Finally the electron goes back to the valence band and
recombines with the hole (step 4). The Feynman diagram corresponding to
that process is represented in figure 1b.
Due to the double resonance process, the mathematical expression that describes the intensity of the D band in graphene is rather complicated [17]. For the purposes of this work we can use a simplified expression taken from [11].
where W () is a factor that contains the matrix elements of the transitions
ocurring during the double resonance processes, L is the average distance
between defects,
is the wave vector corresponding to the maximum of the
D-line and
the line width. The integration is performed on the phonon
wavevectors in the first Brillouin Zone (BZ).As well as our interest is in the main dependence of
with
, it is enough to
consider the central part of the band, then we can substitute the gaussian function by a Dirac delta2 in the integrand of (2) , i.e
Making this change in (2), using the properties of the Dirac delta and changing by , we obtain:
where .
Besides the explicit dependence of with
that is showed in (4), there is
also an implicit dependence because the line width also depends of the
concentration of defects. In fact, for the excited electronic states, those that
participates in the processes originating the D-band, there are two basic
mechanisms of relaxation: scattering by phonons and scattering by defects.
We can then write the following:
The first term in the right side of (5) is the contribution to the line width of the electron-phonon interaction and the second one is the contribution of the electron-defect scattering. It is obvious that the last term should depend of the concentration of defects3. In the case of irradiated graphene the main kind of defects are single vacancies [18]. For this kind of defects, the expression:
Is frequently used [17] , where is the square mean deviation of the
nearest neighbor hopping energy calculated by the tight binding model. This variation is produced by the presence of defects, which change the positions of
the atomic orbitals. Using (6) we can rewrite (4) as follows:
In the case of the G-band, the relaxation of the electronic states that
participates in the Raman process is produced mainly through the scattering of
electrons by phonons. Then, in this case, the contribution to the line width of the
interaction electron-defects can be neglected and we can consider Then, the expression for the intensity of line G can be approximated to:
RESULTS AND DISCUSSION
Using (7) and (8) we can obtain the dependence of with
as follows:
Expression (9) shows clearly the non-monotonous dependence of with the
concentration of defects in irradiated graphene. In obtaining it ,we have not
used any kind of consideration about amorphization of another mechanism,
additional to the basic Raman processes. Then, the experimentally obtained dependence of
with the defect concentration is intrinsic to the nature of the
processes that produces the D-band in the Raman spectra of graphene.
Besides its theoretical value, expression (9) can be used as a function for fitting experimental data, using a, b, c as fitting parameters.
The fitting of the experimental data obtained by Luchesse et al4. [7 ] using the
expression (9), is showed in Figure 2. It was carried out with the fitting tools of
MATLAB that use the determination coefficient, for evaluating the quality of
the approximation. This coefficient is defined by:
where and
are the values of the independent variable (
in this
work ) obtained from the model and from the experimental data, respectively;
and y is the mean value of the experimental values. In our case
= 0,96.
Furthermore, the relative error in the determination of the parameters a, b, c is
less than 30%. This indicates that, although we have introduced several
approximations in the model, formula (9) is an acceptable representation of the
function
=
.
CONCLUSIONS
It has been demonstrated that the non-monotonous dependence of with
the defect concentration (
) is an intrinsic feature of the Raman spectra of
irradiated graphene and it is not necessarily related to the amorphization of the
samples. A simple mathematical relationship between these magnitudes was obtained.
This expression can be used to fit experimental data of the Raman spectra of
irradiated graphene in relation with the concentration of defects.
REFERENCES
1. GEIM AK. Graphene: Status and prospects. Science. 2009; 324(5934): 1531-1534.
2. CASTRO NETO AH, GUINEA F, PERES NMR, et. al. The electronic properties of grapheme. Rev. Mod. Physi. 2009; 81(1): 109-162.
3. COOPER DR, D'ANJOU B, GHATTAMANENI N, et. al. Experimental review of graphene. ISRN Condensed Matter Physics. 2012; 2012(article ID 501686): 1-56.
4. YAO J, SUN Y, YANG M, DUAN Y. Chemistry, physics and biology of graphene-based nanomaterials: new horizons for sensing, imaging and medicine. J. Mater. Chem. 2012; 22(29): 14313-14329.
5. TAPASZTO L, DOBRIK G, NEMES-INCZE P, et. al. Tuning the electronic structure of graphene by ion irradiation. Phys.Rev. B. 2008; 78(1-4): 233407.
6. LETHINEN O, KOTAKOSKI J, KRASHENINNIKOV AV, et. al. Effect of ion bombardment on a two-dimensional target. Phys. Rev. B. 2010; 81(1-4): 153401.
7. LUCCHESE MM, STAVALE F, MARTINS FERREIRA EH, et. al. Quantifying ion-induced defects and Raman relaxation length in grapheme. Carbon. 2010; 48(5): 1592-1597.
8. CODORNIU PUJALS D, AGUILERA CORRALES Y, BALDASARRE F. Calculation of the number of atoms displaced during the irradiation of monolayer grapheme. J. Radioanal. Nucl. Chem, 2011; 289(1): 167-172.
9. MALARDAL M, PIMENTAM A, DRESSELHAUS G, DRESSELHAUS MS. Raman spectroscopy in grapheme. Physics Reports. 2009; 473(5-6): 51-87.
10. FERRARI C. AND ROBERTSON J. Interpretation of Raman spectra of disordered and amorphous carbon. Physical Review B. 2000; 61(20): 14095-14101.
11. SAITO R, HOFMANN M, DRESSELHAUS G, et. al. Raman spectroscopy of graphene and carbon nanotubes. Advances in Physics. 2011; 60(3): 413-550.
12. TUINSTRA F, KOENIG JL. Raman spectrum of graphite. J Chem. Phys. 1970; 53(3): 1126–30.
13. TEWELDEBRHAN D, BALANDIN A A. Modification of graphene properties due to electron-beam irradiation. Appl. Phys. Lett. 2009; 94(1): 013101-013103.
14. BUCHOWICZ G, STONE P, ROBINSON JT, et. al. Correlation between structure and electrical transport in ion-irradiated graphene grown on Cu foils. [artículo en línea]. Appl. Phys. Letters. 2011; 98(3): 032102. <http: arxiv.org/pdf/1012.4060> [consulta: feb 2013]
15. KRASHENNINIKOV AV, NORDLUND K. Ion and electron irradiation-induced effects in nanostructured materials. J.Appl. Phys. 2010; 107(7): 071301.
16. THOMSEN C, REICH S. Double Resonant Raman Scattering in Graphite. Phys. Rev. Lett. 2000; 85(24): 5214-5217.
17. VENEZUELA P, LAZZERI M, MAURI F. Theory of double-resonant Raman spectra in graphene: intensity and line shape of defect-induced and two-phonon bands [artículo en línea]. Phys Rev B. 2011; 84(3): 035433. <http: arXiv: 1103.4582> [consulta: feb 2013]
18. LETHINEN O, KOTAKOSKI J, KRASHENNINIKOV AV, et. al. Effects of ion bombardenment in a two-dimensional target. Physical Review B. 2010; 81(15): 153401.
Recibido: 15 de noviembre de 2012
Aceptado: 25 de abril de 2013
1The so-called Tuinstra-Konig relationship was formulated as ,where
is a factor that depends of the wave length of the exciting laser and L is the distance between the crystallites, that is approximately the same as the average distance between defects. Then,
2The factor is introduced to assure that the integral of both functions is the same
3The anharmonic effects give an additional contribution to the line-width. As this contribution is not related to defects, we will include it in .
4The experimental results reported in [7 ] constitute one of the most complete set of data Raman intensities in irradiated graphene published so far.