# Influence of Carrier-Carrier Scattering on Electron Transport in Monolayer Graphene

## Abstract

The influence of electron-electron scattering on the distribution function and transport characteristics of intrinsic monolayer graphene is investigated via an ensemble Monte Carlo simulation. Due to the linear dispersion relation in the vicinity of the Dirac points, it is found that pair-wise collisions in graphene do not conserve the ensemble average velocity in contrast to conventional semiconductors with parabolic energy bands. Numerical results indicate that electron-electron scattering can lead to a decrease in the low field mobility by more than 80 % for moderate electron densities. At high densities, the impact gradually diminishes due to increased degeneracy.

Since graphene was first realized experimentally in 2004, (1) it has attracted significant interest due to its unique electronic properties. At low electron energies near the inequivalent Dirac points, there is no gap between the valence and conduction bands and the dispersion of the energy bands is linear. (2) Extremely high electrical mobilities have been reported in suspended graphene, exceeding 10 cm/Vs near room temperature, (3) suggesting potential applications to ultrahigh frequency electronic devices. (4) The carrier density, which is controlled by the gate voltage, can be expected to vary by orders of magnitude. An interesting consequence of the linear energy dispersion is that the ensemble average velocity is not necessarily conserved upon an electron-electron (e-e) scattering event. Accordingly, inter-carrier collisions deserve careful consideration in the determination of the transport properties.

Das Sarma et al. (5) found the inelastic e-e scattering rate and mean free path in graphene through the analysis of the quasiparticle self-energies. The scattering rate, calculated for electron densities from cm, was found to be of the same order of magnitude for electron-phonon scattering rates evaluated in the deformation potential approximation (DPA). (6) Several authors have considered the electronic transport properties of graphene based on approaches such as the Monte Carlo simulation (see, for example, Ref. (7)); however, the effects of e-e scattering have yet to be addressed to the best of our knowledge. In the present study, we examine the influence of this interaction mechanism in intrinsic monolayer graphene at room temperature. A full-band ensemble Monte Carlo method is used for accurate analysis of the distribution function as well as its macroscopic manifestations, particularly, the electron low-field mobilities and drift velocities.

In both bulk and two-dimensional (2D) conventional semiconducting systems, e-e scattering has been well studied as is documented in the literature. (8) During an e-e scattering event, both the energy and momentum are conserved. In a parabolic band structure, common in conventional semiconductors, momentum conservation directly leads to the conservation of velocity. This can be readily shown by multiplying the momentum conservation equation by , which gives . It is then clear that, in a material with a parabolic band structure, e-e scattering has no direct effect on the drift velocity (which is an ensemble averaged quantity). The situation becomes different in graphene, where the energy dispersion in the region of the inequivalent Dirac points is linear, therefore removing the constraint that momentum conservation leads directly to velocity conservation.

The problem of e-e scattering in nonparabolic bands was previously treated by Bonno et al. (9) Here we use a similar approach, making it specific to the linear band structure of graphene. The transition probability for two electrons, with wavevectors and , is given by Fermi’s golden rule

(1) | |||||

where is the electron energy dispersion relation and and are the wavevectors of the final states after scattering. The matrix element of interaction , when accounting for exchange scattering (i.e., the indistinguishability of electron pairs with like spin after collision), is given by (10)

(2) |

where the Coulomb scattering matrix between two electrons (; ) is

(3) |

with , is the static dielectric function, is the angle between the wavevectors and , and is the normalization area. The corresponding element (; ) can be expressed likewise with . The dielectric function, which is valid for densities larger than cm, is considered in the random-phase approximation. (11)

The scattering rate, including the effect of degeneracy, is then obtained as the double sum over all the final states and the partner electrons (with determined by momentum conservation)

(4) | |||||

where the ’s are the occupation probabilities. Within the Monte Carlo simulation, the Pauli exclusion principle is treated using the standard rejection technique.

In graphene, the energy dispersion of the conduction band is linear in the vicinity of the Dirac points; i.e., , with the Fermi velocity ( cm/s). Therefore, the electrons behave as massless particles and can not be analyzed by the center of mass concept. (12) However, by considering the equations of momentum and energy conservation, it can be shown that all the possible final states after scattering lie on an ellipse defined by and . As shown in Fig. 1, the initial point (i.e., the coordinate centrum at the or point in the Brillouin zone) and the terminal point of the vector specify the two foci of the ellipse, with the length of long axis of the ellipse determined by .

Our Monte Carlo model applied in the simulation includes the full electronic band structure as calculated in the tight binding approximation. The electron-phonon scattering rates and phonon dispersion relations, previously obtained from a first principles approach based on density functional perturbation theory, are included for all six branches. (13) All electron-phonon transition possibilities, including intravalley and intervalley, are accounted for as well as the e-e scattering. Finally, since our simulations include large electron densities we account for degeneracy, where Pauli exclusion is non-negligible, by implementation of the rejection technique in the selection of the final state after scattering. (14)

In Fig. 2, we show the calculated e-e scattering rate, assuming a Fermi-Dirac distribution of initial states and that all final states are all available. The scattering rate is found to have a maximum of approximately 40-80 ps, for the energy range and electron densities considered, which is at least an order of magnitude larger than the total electron-phonon scattering rate. (13) In Fig. 3, we examine the effect of e-e scattering on the drift velocity from the low-field to saturation regimes. It is evident that there is a more than 80 degradation of the low field mobility, which can be estimated by taking the slope of the drift velocity in the low field region, for an electron concentration of cm. Additionally, there is a degradation of the saturation drift velocity by 16 . For low fields, when electrons rarely gain enough energy to emit an optical phonon, acoustic phonon and e-e scattering events dominate the electron dynamics and determine the mobility. (5) As the field is raised, optical phonon scattering becomes more prevalent and the effect of e-e scattering is somewhat mitigated. The dominance of e-e scattering at low energies, relative to the electron-phonon scattering rate, suggests that other scattering mechanisms (i.e., substrate dependent surface polar phonons, remote impurities, etc.) may tend to wash out the resulting degradation.

For larger electron densities, cm and cm, this effect is suppressed by the reduction of available states with high and low energies. This can be seen from the electron energy distribution functions, given in Fig. 4, at an electric field of kV/cm for the highest and lowest carrier densities considered. Just as in conventional semiconductors, e-e scattering results in an extension of the thermal tail of the distribution function, while simultaneously shifting its peak backward. Consequently, though an increase in electron density results in an increased e-e scattering rate, the distribution function is less affected than in the moderate density case. This results in a larger degradation of the low-field mobility and saturation drift velocity for lower concentrations, as is evident in Fig. 3.

Figure 5 presents how the selection of final state pairs effects the velocity of the ensemble which, as stated previously, is not necessarily conserved. On average, the momentum of an initial electron pair will be directed along the electric field which corresponds to the long axis of the ellipse, in Fig. 1. The curve corresponding to , where is the energy difference of the final state electron pair, and , has maximal values when is equal to or . Also shown in Fig. 5, is the value of the final state drift velocity in the direction, , normalized by the fermi velocity . Opposite to , the velocity in the direction of the electric field has minima at is equal to or . This results in a decrease of the the ensemble drift velocity for moderate electron densities, on the order of = 110 cm, since there are available high and low energy states to scatter into. At larger densities, the lowest energy states are already mostly occupied, thereby preventing transition to the energy states which most greatly reduce the drift velocity. As a result, the effect of e-e scattering is most pronounced for = 110 cm, and largely suppressed for = 110 cm.

In summary, our results clearly indicate that due to e-e scattering, there is an electron density dependent degradation of the low-field mobility and saturation drift velocity in intrinsic monolayer graphene. This effect is most evident at moderate electron densities, while being suppressed at higher densities due to the unavailability of low energy states. For moderate electron densities, it is predicted that the low-field mobility and saturation drift velocity are degraded on the order of 80 % and 16 %, respectively.

This work was supported, in part, by the DARPA/HRL CERA, US ARO, and SRC/FCRP FENA programs.

### References

- K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, Science 306, 666 (2004).
- A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 (2009).
- K.I. Bolotin, K.J. Sikes, Z. Jiang, M. Klima, G. Fudenberg, J. Hone, P. Kim and H. L. Stormer, Solid State Commun. 146, 351 (2008).
- Y. M. Lin, C. Dimitrakopoulos, K. A. Jenkins, D. B. Farmer, H. Y. Chiu, A. Grill, and Ph. Avouris, Science 327, 622 (2010).
- S. Das Sarma, E. H. Hwang, and W. K Tse, Phys. Rev. B 75, 121406 (2007); E. H. Hwang, B. Y.-K. Hu, and S. Das Sarma, Phys. Rev. B 76, 115434 (2007);Wang-Kong Tse, E. H. Hwang, and S. Das Sarma, Appl. Phys. Lett. 93, 023128 (2008).
- R. S. Shishir, F. Chen, J. Xia, N. J. Tao, and D. K. Ferry, J. Comput. Electron. 8, 43 (2009)
- A. Akturk and N. Goldsman, J. Appl. Phys. 103, 053702 (2008); R. S. Shishir and D. K Ferry, J. Phys. Condens. Matter, 21 344201 (2009); J. Chauhan and J. Guo, Appl. Phys. Lett. 95, 023120 (2009).
- O. Bonno, J. L. Thobel, and F. Dessenne, J. Appl. Phys. 97, 043702 (2005); U. Bhapkar and M. S. Shur, J. Appl. Phys. 82, 1649 (1997); M. V. Fischetti and S. E. Laux, Phys. Rev. B 38, 9721 (1988); S. M. Goodnick and P. Lugli, Phys. Rev. B 37, 2578 (1988).
- O. Bonno and J. L. Thobel, J. Appl. Phys. 104, 053719 (2008).
- A. Mošková and M. Moško, Phys. Rev. B 49, 7443 (1994).
- E. H. Hwang and S. Das Sarma, Phys. Rev. B 75, 205418 (2007).
- J. Sabio, F. Sols, and F. Guinea, Phys. Rev. B 81, 045428 (2010).
- K. M. Borysenko, J. T. Mullen, E. A. Barry, S. Paul, Y. G. Semenov, J. M. Zavada, M. Buongiorno Nardelli, and K. W. Kim, Phys. Rev. B 81, 121412(R) (2010).
- P. Lugli and D. K. Ferry, IEEE Trans. Electron Devices 32, 2431 (1985).