Molecular simulation of temperature-programmed desorption

We use accelerated molecular dynamics to probe the thermal desorption of n-alkanes from the Au(1 1 1) and C(0 0 0 1) surfaces. Studying an alkane series ranging from CH₄ to C₁₆H₃₄, we find that the desorption prefactor increases with increasing chain length for the short chains until a certain chain length is reached when it becomes essentially constant and independent of chain length. We can understand the dependence of the preexponential factor on alkane chain length in terms of conformational changes within the alkane molecules. For the shorter molecules, the desorption temperatures probed in experimental temperature-programmed desorption studies lie below those for which torsional motion is activated. Thus, the short alkanes can be treated as rigid rods, and the loss in translational and rotational entropy upon adsorption leads to a preexponential factor that increases with increasing chain length. As the alkane chain length increases, the binding energy and the experimental desorption temperatures also increase. Thus, torsional motion is activated to an extent that increases with increasing chain length. This torsional activation increases the entropy of adsorption and counteracts the entropy decrease due to a loss of translation and rotation. This leads to a virtually constant prefactor for sufficiently long chains. Our findings are consistent with experimental data for the thermal desorption of alkanes from the Au(1 1 1), Pt(1 1 1), MgO(1 0 0) and C(0 0 0 1) surfaces.