CC BY
42
THE ELECTRON-PHONON MATRIX ELEMENT IN THE DIRAC POINT OF GRAPHENE
!'2S. V. Koniakhin, 1'2'SE. D. Eidelman
1 Ioffe Physical-Technical Institute of the Russian Academy of Sciences - 194021
St. Petersburg, Russia 2 St. Petersburg Academic University - Nanotechnology Research and Education Centre of the Russian Academy of Sciences - 194021 St. Petersburg, Russia 3 St. Petersburg Chemical Pharmaceutical Academy - 197022 St. Petersburg, Russia
kon@mail.ioffe.ru
PACS 63.22.Rc,63.20.kd
The chief aim of this paper is to derive the matrix element of electron-phonon interaction in graphene as a function of phonon wave vector. The tight binding model, harmonic crystal approximation, and deformation potential approximation were employed for obtaining the matrix element. Required microscopic parameters are available in the current literature. This technique allows the most precise derivation of the electron-phonon matrix element in graphene based on the semiempirical models. Scattering of electrons from the Dirac point is considered as most important. The 2D plots of the e-ph matrix element absolute value as a function of the phonon wave vector for in-plane modes are given as a result. These plots show the high anisotropy of the e-ph matrix element and singularities in high symmetry points. The results are in agreement with the long-wavelength approximation.
Keywords: graphene, electron-phonon interaction.
1. Introduction
Electron-phonon interaction can affect transport properties of extremely pure graphene [1,2] and consequently its investigation is important for creating electronic devices based on graphene layers, particularly thermoelectric devices [3,4]. Raman spectroscopy is one of the most comprehensive techniques for studies of graphene [5, 6] and double resonant Raman peaks are footprints of electronic and vibrational properties of the crystal. This encourages the research of electron-phonon interaction in graphene. The knowledge of vibrational properties, electron wave functions and deformation potential in total can lead to derivation of the electron-phonon matrix element.
2. Harmonic crystal approximation
Graphene is a two dimensional crystal with a honeycomb lattice of carbon atoms. One of the alternatives for the unitary cell is rhombic one. It contains two atoms in basis: A and B types. They have different directions of the chemical bonds, which causes their inequivalence. Reciprocal lattice for the two-dimensional hexagonal crystal is also hexagonal lattice, which is in addition rotated over 60 degrees. Brillouin zone carries symmetrical properties of the crystal and consequently it must be hexagon-shaped. To specify the position of each atom, one can write R(w, s) = RM + rs, where RM is the position of unitary cell in crystal and rs is the position of atom of s type in the unitary cell.
The vibrational properties of graphene can be described on the basis of the harmonic crystal approximation [7]. This method starts from the Newton's second law for each atom in crystal:
mxa (u, s) = - ^ (u - U, s, s')xp (u', s1), (1)
u's'P
where xa(u,s) is a component of displacement of the s type atom in the u cell, and are interatomic force constants. The infinite system of the classical motion equations after substitution are:
xa (q, u, s) = —Uwa(q, s) exp (iqR (u, s) - iu(q)t) (2)
m
reduces to the eigenvector problem:
u2(q)wa(q,s) = Cal3{ q )wl3(q,s'), (3)
s'P ^ '
for the dynamical matrix Cap that is defined as:
Caß\ q ' ) = ~ y^ (u " U, s, S) ■ exp(«q • (R(u, s) - R(u', s'))). V ss I m \ / -.,/
q ) = ~ y J (u - u''s's') ■ exp(iq • (R(u>s) - R(u'> s')))- (4)
1 °° ' m / ^
It is useful to introduce the vector Ws(q) = (wx(q,s),wy(q,s))T. Eigenvalues for the dynamical matrix give squared frequencies for phonons of different phonon modes and eigenvectors represent corresponding directions and phases of atomic displacements. For a crystal with 2 atoms in its basis, the dynamical matrix dimensions are 6 by 6. But out of plane modes can be separated, and therefore, the in-plane modes can be described by a 4 by 4 dynamical matrix. Interatomic force constants utilized to construct the dynamical matrix are available in current literature [8-12]. The density-functional theory (DFT) is used in [10] to obtain the interatomic force constants for neighbors, up to twentieth, and we use here these data as the most complete.
3. Derivation of electron-phonon matrix element
The wave functions of electron with wave vector k in graphene can be written as
^k(r) = —t^Yj Cs(k)^(r - Ru - rs) ■ exp(ik(Ru + rs)). (5)
u,s
Coefficients Cs(k) define the true Hamiltonian eigenstates with wave vector k, and ^(r) is the atomic wave function. Most of important optical and electronic properties of graphene are defined by the bands that appear as a result of interaction between electrons from n orbitals of carbon atoms. Consequently, ^(r) here is close to the n orbitals.
The phonon is a wave of atomic displacements from their equilibrium positions and consequently it causes a deviation potential from the potential of unperturbed crystal. It leads to the possibility of electron scattering on corresponding perturbation. For small atomic displacements, the perturbation potential is directly proportional to these displacements. Thereby, one can write the scalar product of atomic displacement and gradient of the certain potential in the following expression:
ÔVq(r) = J] W (q) • VV(r - R u'' — rs«)) • exp(ïq(Ru» + rs«)).
(6)
In the hard ion approach the potential of an atom is assumed to be rigid and consequently in relation (6) V is directly the potential of the atom that is used to calculate band structure in the material, e.g. to construct tight-binding matrix elements in the tight binding model. Indeed, the form of atomic potential changes when the atom is displaced and it leads to the concept of deformation potential. In [13], the deformation potentials in nanotubes and graphene were studied on the basis of ab initio calculations.
The Electron-phonon matrix element defines the process of electron scattering with wavevector k into the state with wavevector k' on introduced potential:
Me_ph(k ^ k' = k + q) = Me_ph(k, q)
dr , (r)#Vq(r)-0k(r).
(7)
Basic formulae for calculation of the electron-phonon matrix element are given in [13] in clear form, and where possible, designations from [13] are used. However, in [13], the authors mainly investigate nanotube transport properties and do not give the electron-phonon matrix element as a function of the wave vector or represent 2D plots of Me_ph . For instance, the precise Raman spectra calculation requires such information.
After substitution of relations (2), (5), and (6) to (7) one obtains:
Me-ph(k, q) = ^ J] J] Cs>t(k')Cs(k)x
u,s s' u",s"
exp(ïk(RAu + rs) — ¿(k + q)Ru' + zq(RAu'' + )) • (Wy (q) • g(s'; Au, s; Au'', s'')) , (8) where
g(u';Au,s;Au'',s'')
dr [^(r)VV(r — (Rau'' + rs» — rs')Mr — (Rau + rs — Rs'))] . (9)
Fig. 1 represents relative positions of atoms that figures in the integral in relation (9).
Fig. 1. Diagram illustrating the structure of terms in the expression for the deformation potential.
As was mentioned previously, the relevant electron states in graphene are the result of n orbital interactions. Due to the geometry of the n-n bond, the energy of interaction
u'' ,s
between carbon atoms changes significantly only if atoms are displaced along the bond axis. Mathematically it means that the direction of the potential gradient lies along vectors between atom s' and its neighbors. Also, only interactions between nearest neighbors are dominant. Consequently, the two types of terms appear. The first one is so-called off-site term, where Au = 0, s = s':
g(s';0,s';Au'',s'') = ,RAu'' + ^ - ^ ■ A(|Ra«« + ^ - r„|), (10)
|Ra«" + rs" - rs' |
This corresponds to perturbation of the diagonal terms of Hamiltonian. The second type are on-site terms, where Au'' = 0, s'' = s' or Au'' = Au, s'' = s':
RAu+rs-rs'
u+rs-V | ^ a(|RAu + is (tt)
R.A«+rs-rs, , „ r IN i11;
g(s'; Au, s; 0, s') = |RA:+r:-r:'| ■ a(|Rau + rs - rs'|),
g(s'; Au, s; Au, s) = - |RAU+rS—r:'| ■ a(|Rau + rs - rs'|)
They define the off-diagonal perturbation of the Hamiltonian. Numeric values of a and A are given in [13].
4. Results and discussion
Electrons near the Fermi level make the main contribution to the transport phenomena in graphene. Due to the shape of electron dispersion, these electrons have a wave vector lying near the Dirac point (K point or conical point) of the Brilloiun zone. Visible light can excite electron hole pairs also in the region of the Dirac point only. Consequently, scattering of electrons from the Dirac point of the Brillouin zone attracts most interest.
Fig 2. shows the absolute value of electron-phonon matrix element as a function of the phonon wave vector q for TA,LA,TO and LO modes. Initial electron wave vector k corresponds to the Dirac point (k = (4^/3a0, 0)). High anisotropy of the electron-phonon matrix element is detected. Patterns obey D3h symmetry, which is strictly necessary for the Dirac point of Brillouin zone in graphene. The directions in Brillouin zone and segments of hexagonal Brilloiun zone with high and low values of the electron-phonon matrix element tolerably match to the ones from [14]. In [14], the authors used density-functional perturbation theory (DFPT) to calculate the electron-phonon matrix element, but the details of their calculations are unfortunately very brief.
Next, there are estimations [15] of the matrix element of scattering of electron in Dirac point within the long wave length approximation. For small phonon wave vectors q, the following behavior was indicated:
ii jTA,LA
M h - q
j.JTOLO , (12)
Me—ph - const.
It can be seen that the obtained results agree well with estimations from [15].
Introduced derivation of the electron-phonon matrix element is based on the most complete and full sets of microscopic material parameters, and thus, it provides the most precise semiempirical approach at this point. The results obtained can be used for improvement of Raman spectra calculation. Taking into account the high anisotropy of the electron-phonon matrix element is also important for investigation of the electron drag by phonons in graphene.
TA LA
Fig. 2. Absolute value of electron-phonon matrix element for TA,LA,TO, and LO modes as a function of wave vector. The initial state of the electron corresponds to the Dirac point. Lighter areas correspond to larger values of matrix elements.
Acknowledgements
We would like to thank A. Ya. Vul' and A. T. Dideykin for their attention. This work was supported by the Dynasty foundation Grant for PhD students and by the St. Petersburg government grant for PhD students. The work was partly supported by programs of the Presidium of the Russian Academy of Sciences "Fundamental Foundations of the Technologies of Nanostructures and Nanomaterials" and "Quantum Mesoscopic and Unordered Systems".
References
[1] Tse Wang-Kong, Das Sarma S. Energy relaxation of hot Dirac fermions in graphene. Phys. Rev. B., 79, P. 235406 (2009, Jun.).
[2] Mariani Eros, von Oppen Felix. Temperature-dependent resistivity of suspended graphene. Phys. Rev. B, 82, P. 195403 (2010, Nov.).
[3] Wang Deqi, Shi Jing. Effect of charged impurities on the thermoelectric power of graphene near the Dirac point. Phys. Rev. B, 83, P. 113403 (2011, Mar.).
[4] Checkelsky Joseph G., Ong N. P. Thermopower and Nernst effect in graphene in a magnetic field. Phys. Rev. B, 80. P. 081413 (2009, Aug.).
[5] A. C. Ferrari, J. C. Meyer, V. Scardaci et al. Raman Spectrum of Graphene and Graphene Layers. Phys. Rev. Lett., 97, P. 187401 (2006, Oct.).
[6] Venezuela Pedro, Lazzeri Michele, Mauri Francesco. Theory of double-resonant Raman spectra in graphene: Intensity and line shape of defect-induced and two-phonon bands. Phys. Rev. B., 84, P. 035433 (2011, Jul.).
[7] Born M., Huang K. Dynamical Theory of Crystal Lattices. Oxford classic texts in the physical sciences. Clarendon Press (1996).
[8] Al-Jishi R., Dresselhaus G. Lattice-dynamical model for graphite. Phys. Rev. B., 26, P. 4514-4522 (1982, Oct.).
[9] A. Grüneis, R. Saito, T. Kimura et al. Determination of two-dimensional phonon dispersion relation of graphite by Raman spectroscopy. Phys. Rev. B., 65, P. 155405 (2002, Mar.).
[10] Dubay O., Kresse G. Accurate density functional calculations for the phonon dispersion relations of graphite layer and carbon nanotubes. Phys. Rev. B., 67, P. 035401 (2003, Jan.).
[11] Wirtz Ludger, Rubio Angel. The phonon dispersion of graphite revisited. Solid State Communications, 131(3-4), P. 141 (2004).
[12] Hui Wang, Yufang Wang, Xuewei Cao et al. Vibrational properties of graphene and graphene layers. Journal of Raman Spectroscopy, 40(12), P. 1791-1796 (2009).
[13] J. Jiang, R. Saito, Ge. G. Samsonidze et al. Electron-phonon matrix elements in single-wall carbon nanotubes. Phys. Rev. B, 72, P. 235408 (2005, Dec.).
[14] K. M. Borysenko, J. T. Mullen, E. A. Barry et al. First-principles analysis of electron-phonon interactions in graphene. Phys. Rev. B, 81, P. 121412 (2010, Mar.).
[15] Manes J. L. Symmetry-based approach to electron-phonon interactions in graphene. Phys. Rev. B., 76, P. 045430 (2007, Jul.).
