The purpose of this paper is to simulate charge transport in monolayer graphene on a substrate made of hexagonal boron nitride (h-BN). This choice is motivated by the fact that h-BN is one of the most promising substrates on account of the reduced degradation of the velocity due to the remote impurities.
The semiclassical Boltzmann equations for electrons in the monolayer graphene are numerically solved by an approach based on a discontinuous Galerkin (DG) method. Both the conduction and valence bands are included, and the inter-band scatterings are taken into account as well.
The importance of the inter-band scatterings is accurately evaluated for several values of the Fermi energy, addressing the issue related to the validity of neglecting the generation-recombination terms. It is found out that the inclusion of the inter-band scatterings produces sizable variations in the average values, like the current density, at zero Fermi energy, whereas, as expected, the effect of the inter-band scattering becomes negligible by increasing the absolute value of the Fermi energy.
The correct evaluation of the influence of the inter-band scatterings on the electronic performances is deeply important not only from a theoretical point of view but also for the applications. In particular, it will be shown that the time necessary to reach the steady state is greatly affected by the inter-band scatterings, with not negligible consequences on the switching on/off processes of realistic devices. As a limitation of the present work, the proposed approach refers to the spatially homogeneous case. For the simulation of electron devices, non-homogenous numerical solutions are required. This last case will be tackled in a forthcoming paper.
As observed in Majorana et al. (2019), the use of a Direct Simulation Monte Carlo (DSMC) approach, which properly describes the inter-band scatterings, is computationally very expensive because the valence band is highly populated and a huge number of particles is needed. Even by simulating holes instead of electrons does not overcome the problem because there is a certain degree of ambiguity in the generation and recombination terms of electron-hole pairs. The DG approach, used in this paper, does not suffer from the previous drawbacks and requires a reasonable computing effort.
