In this work, we propose a theoretical and computational model for taking into account the anisotropic structure of highly oriented pyrolitic graphite (HOPG) in the Monte Carlo simulations of charge transport. In particular, the dielectric characteristics, such as the inelastic mean free path and energy losses, are treated by linearly combining the contributions to these observables along the two main orthogonal directions identifying the layered crystalline structure of HOPG (along the layer plane and perpendicular to it).
Energy losses are evaluated from ab initio calculations of the dielectric function of the system along these two perpendicular directions. Monte Carlo simulated spectra, obtained with our anisotropic approach, are compared with acquired experimental data of reflection electron energy loss and secondary electron spectra, showing a good agreement.
These findings validate the idea of the importance of considering properly weighted interplanar and intraplanar interactions in the simulation of electron transport in layered materials.