Calculation of electronic properties of multilayer graphene with Monte Carlo method

In this study, the electronic transport properties of bilayer graphene is investigated by an ensemble Monte Carlo method. The bilayer graphene has a quadratic energy dependence on wave vector near the points known as Dirac points in the reciprocal lattice. For bilayer graphene the scatterings due to acoustic and optic phonons and ionized impurities are taken into account. Velocity-time and steady state velocity-applied field curves are obtained and from the slope of velocity-field curves at low fields, the low field mobility of bilayer graphene is obtained. The dependence of mobility of bilayer graphene on temperature, electron concentration, impurity concentration, acoustic and optic deformation constants is investigated and it is observed that the most important mechanism limiting the mobility is the phonon scattering.