Collective electronic excitations (plasmons) in single-layer and few-layer graphene have been studied extensively in the past few years. In particular, the dispersion and nature of the π plasmon peak in free-standing single-layer graphene was investigated by means of momentum-resolved electron energy-loss spectroscopy (MREELS) [1-3]. Besides, it was also studied how the transition from mono- to multilayer graphene changes the shape of EEL spectra. Within a layered electron-gas (LEG) model, this transition can be modeled very precisely [4,5].
For other 2D materials such as transition metal dichalcogenides (TMDs), however, this approach may be highly inaccurate. We therefore evaluate more precise model calculations for multilayer systems. Our calculations are based on time-dependent DFT calculations for the individual monolayers. The plane-wave pseudopotential DFT code abinit  is used to simulate the ground state within local density approximation (LDA). The linear response is then calculated within random-phase approximation (RPA) using the dp-code . Compared to the LEG model where the monolayers are assumed to be perfectly 2-dimensional and homogeneous, we preserve the layers’ microscopic structure in our calculations.
We demonstrate that only by taking the finite thickness of the constituent monolayers into account, the spectra of multilayer MoS2 can be modeled correctly (see the figure). Our calculations can be also applied to arbitrary van-der-Waals heterostructures. Our results are directly compared to MREEL spectra recorded with a Zeiss Libra 200 based TEM prototype (“SALVE I”, [8,2]) equipped with a monochromator and an Ω type in-column energy filter. This allows for the verification of our calculations over a large range of energy losses (0-40 eV) and different momentum transfers within the Brillouin zone. 
 Kinyanjui et al., EPL 97(5), 2012
 Wachsmuth et al., Phys. Rev. B 88, 2013
 Liou et al., Phys. Rev. B 91, 2015
 Jovanović et al., Phys. Rev. B 84, 2011
 Wachsmuth et al., Phys. Rev. B 90, 2014
 Gonze et al., Comput. Phys. Commun. 180(12), 2009
 Olevano et al., www.dp-code.org
 Kaiser et al., Ultramicroscopy 111(8), 2011
 We acknowledge financial support by the German Research Foundation (DFG) and the Ministry of Science, Research and Arts (MWK) of the state Baden-Württemberg within the Sub-Angstrom Low-Voltage Electron Microscopy project (SALVE).
To cite this abstract:Michael Mohn, Ralf Hambach, Philipp Wachsmuth, Ute Kaiser; Model calculations for low-loss EEL spectra of 2D multilayer systems. The 16th European Microscopy Congress, Lyon, France. https://emc-proceedings.com/abstract/model-calculations-for-low-loss-eel-spectra-of-2d-multilayer-systems/. Accessed: September 18, 2020
EMC Abstracts - https://emc-proceedings.com/abstract/model-calculations-for-low-loss-eel-spectra-of-2d-multilayer-systems/