The recoil energy
is given by
where
represents the mean recoil energy
and
q is the transferred momentum for a target atom with mass
M.
is the spread in the recoil energy due to atomic vibrations. A Gaussian distribution can therefore be observed, whose standard deviation
is given by [
5,
11,
26]
where we have specified the mean kinetic energy of the target atoms with
.
As far as the Monte Carlo method is concerned, it is well known [
12,
27,
28] and therefore we limit ourselves here to the description of its main features. An electron beam irradiates the target surface (which lies in the plane
) with the primary energy
and the angle of incidence
. The elastic scattering cross-section
is calculated by
. In this equation,
stands for the elastic scattering cross-sections of carbon,
for the elastic scattering cross-sections of hydrogen and
and
for their respective atomic concentrations. The inelastic scattering cross section
is calculated by
. In this equation,
is the inelastic mean free path of the electrons and
N is the number of molecules per unit volume. The probabilities of elastic and inelastic scattering are given by
and
, respectively. Having calculated the mean free path of the electrons as
, we can obtain the step length
between the collisions by
where
is a random number that is sampled with a uniform distribution between 0 and 1. The scattering angle
is determined by
where
is the cumulative probability of elastic scattering. The azimuth angle
is sampled uniformly between 0 and 2
. For the choice between elastic or inelastic collisions, a random number
is generated, which is sampled with a uniform distribution between 0 and 1. If
, the collision is inelastic. In this case, the simulation of the electron’s trajectory is finished, as it is no longer of interest to follow the trajectory of this electron. If, on the other hand,
, then the collision is elastic. If this is the case, another random number
is sampled uniformly between 0 and 1 to determine the type of elastic collision and the recoil energy. If, in particular,
then an electron-carbon collision takes place and the recoil energy is given by
The quantities
and
describe the Doppler broadening. Since they are determined using random numbers derived from a Gaussian distribution with the standard deviation calculated according to Equation (
5), they can be positive or negative. Please note that
and
represent the atomic masses of carbon and hydrogen, respectively.
m stands for the electron mass and
E for the electron energy.