a)
For the predefined gases except H2, 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 (H2, D2, N2, and O2).
b) Eads0 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, |Eads0|, 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.) |Eads0| 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 Evib = hν (above a few kBT) 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 Δθ0, or a list of initial coverages θ0(1), θ0(2),... (best in descending sequence).
These coverages should be less than the θ0 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 interaction 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
ρ(|Eads0|), to visualize how it evolves from ρraw to the final result.
Hide footnotes
Author: Michael Schmid, IAP/TU Wien.
This work is licensed under
GNU General Public License (GPLv3) and
Creative Commons Attribution-ShareAlike 4.0 (CC-BY-SA 4.0).
The plot library Plotly is licensed under the MIT license.
When using and/or modifying this program for scientific work please cite our paper:
Michael Schmid, Gareth S. Parkinson, and Ulrike Diebold, Analysis of temperature-programmed desorption via equilibrium thermodynamics,
ACS Phys. Chem. Au (2022), doi: 10.1021/acsphyschemau.2c00031