^a) For the predefined gases except H₂, data are inaccurate below 100 K. For lower temperatures, better select 'OTHER' and enter the molecular structure. Dissociation is only supported for molecules dissociating into two atoms of equal type (H₂, D₂, N₂, and O₂).
^b) E_ads⁰ is the adsorption energy for T = 0, i.e., without vibrational energies of the adsorbate. Energy per atom when dissociating. Each line should contain two numbers, |E_ads⁰|, f, where f is the fraction or relative number of adsorption sites with this adsorption energy. (You need not care about normalization; the program will normalize it such that the sum of all f values becomes unity.) |E_ads⁰| values must be given in ascending order. Further data columns will be ignored.
^c) Comma-separated list, leave empty if none. Enter all modes for the adsorbate if the gas is one of the pre-defined ones. For 'OTHER' gases, omit vibration modes that are equal for the gas phase and the adsorbed molecule. Enter vibration modes present only in the adsorbed molecules as positive numbers and modes present only in the gas phase as negative numbers. Modes with high E_vib = hν (above a few k_BT) can be ignored.
^d) Adsorbate entropy not related to site occupation and vibrations, divided by R. E.g., ln(3) for three distinguishable azimuthal orientations. See ^f) for how to enter a formula.
^e) Enter the step Δθ₀, or a list of initial coverages θ₀₍₁₎, θ₀₍₂₎,... (best in descending sequence). These coverages should be less than the θ₀ value of the main calculation.
^f) Variables: θ or th, T, and s (output; optional if a simple expression is given). Declare own variables with var. Predefined constants h (Planck), k (Boltzmann), R (gas), and e (elementary charge, eV/J). Functions: sqrt, sqr, ln or log, exp, pow(base, exponent), and JavaScript Math functions. Examples: sqr(1-θ) or var x=10*(1-th); s=x/(1+x).
^g) Assumes that the ionization probability in the mass spectrometer is inversely proportional to the velocity, which is proportional to √T. This correction should be applied only if the desorbing molecules directly enter the mass spectrometer, without being reflected at any surfaces. The correction is also inaccurate in case of a barrier on adsorption, causing temperature-dependent sticking.
^h) Select 'Expert mode' to control the fitting process used for determination of the energy distribution (Richardson–Lucy method). The raw energy distribution is iteratively refined to better fit the experimental data. A high number of iterations (default 50) can lead to a better fit, but takes more computing time and may also result in isolated desorption energies even in case of a continuous distribution of energies, to better fit the experimental noise. If you get at least one rather sharp peak in the energy distribution and the calculated TPD peak corresponding to it is not as sharp as the experimental one, consider increasing the number of iterations. The noise suppression (default value 5) reduces sharpening broad distributions of the adsorption energy.
In expert mode, the energy distribution plot contains additional information: The raw correction factor is the ratio of experimental and calculated peak intensity in the last iteration of the Richardson–Lucy refinement, and the correction factor actually applied is obtained by smoothing this raw factor. The plot also shows the raw energy distribution ρ_raw and an intermediate result for ρ(|E_ads⁰|), to visualize how it evolves from ρ_raw to the final result. Hide footnotes