We present temperature programmed desorption (TPD) measurements of CO, CH$_4$, O$_2$ and CO$_2$ from the forsterite(010) surface in the sub-monolayer and multilayer coverage regimes. In the case of CO, CH$_4$ and O$_2$, multilayer growth begins prior to saturation of the monolayer peak, resulting in two clearly distinguishable desorption peaks. On the other hand a single peak for CO$_2$ is observed which shifts from high temperature at low coverage to low temperature at high coverages, sharpening upon multilayer formation. The leading edges are aligned for all the molecules in the multilayer coverage regime indicating zero order desorption. We have extracted multilayer desorption energies for these molecules using an Arrhenius analysis. For sub-monolayer coverages, we observe an extended desorption tail to higher temperature. Inversion analysis has been used to extract the coverage dependent desorption energies in the sub-monolayer coverage regime, from which we obtain the desorption energy distribution. We found that owing to the presence of multiple adsorption energy sites on the crystalline surface the typical desorption energies of these small molecules are significantly larger than obtained in previous measurements for several other substrates. Therefore molecules bound to crystalline silicate surfaces may remain locked in the solid state for a longer period of time before desorption into the gas phase.

The thermal desorption of ammonia (NH$_3$) from single crystal forsterite (010) has been investigated using temperature-programmed desorption. The effect of defects on the desorption process has been probed by the use of a rough cut forsterite surface prepared from the cleaved forsterite sample. Several approaches have been used to extract the desorption energy and pre-exponential factor describing the desorption kinetics. In the sub-monolayer coverage regime, the NH$_3$ desorption shows a broad distribution of desorption energies, indicating the presence of different adsorption sites, which results in an apparent coverage-dependent desorption energy. This distribution is sensitive to the surface roughness with the cut forsterite surface displaying a significantly broader distribution of desorption energies compared to the cleaved forsterite surface. The cut forsterite surface exhibits sites with desorption energies up to 62.5 kJ mol$^{-1} $ in comparison to a desorption energy of up to 58.0 kJ mol$^{-1} $ for the cleaved surface. Multilayer desorption is independent of the nature of the forsterite surface used, with a desorption energy of ($25.8\pm0.9$) kJ mol$^{-1} $. On astrophysically relevant heating time-scales, the presence of a coverage dependent desorption energy distribution results in a lengthening of the NH$_3$ desorption time-scale by $5.9\times 10^4$ yr compared to that expected for a single desorption energy. In addition, the presence of a larger number of high-energy adsorption sites on the rougher cut forsterite surface leads to a further lengthening of ca. 7000 yr.