@Article{CiCP-26-114, author = {Armando Majorana, Giovanni Nastasi and Vittorio Romano}, title = {Simulation of Bipolar Charge Transport in Graphene by Using a Discontinuous Galerkin Method}, journal = {Communications in Computational Physics}, year = {2019}, volume = {26}, number = {1}, pages = {114--134}, abstract = {
Charge transport in suspended monolayer graphene is simulated by a numerical deterministic approach, based on a discontinuous Galerkin (DG) method, for
solving the semiclassical Boltzmann equation for electrons. Both the conduction and
valence bands are included and the interband scatterings are taken into account.
The use of a Direct Simulation Monte Carlo (DSMC) approach, which properly
describes the interband scatterings, is computationally very expensive because the valence band is very populated and a huge number of particles are needed. Also the choice
of 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. Often, direct solutions of the Boltzmann equations with a DSMC neglect
the interband scatterings on the basis of physical arguments. The DG approach does
not suffer from the previous drawbacks and requires a reasonable computing effort.
In the present paper the importance of the interband 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 interband scatterings produces huge variations in the average values, as the current, with zero Fermi energy while, as expected, the effect of the interband scattering
becomes negligible by increasing the absolute value of the Fermi energy.