pith. sign in

arxiv: 2511.05448 · v1 · submitted 2025-11-07 · ❄️ cond-mat.mtrl-sci · physics.ins-det

An improved reliability factor for quantitative low-energy electron diffraction

Alexander M. Imre , Lutz Hammer , Ulrike Diebold , Michele Riva , Michael Schmid This is my paper

Pith reviewed 2026-05-17 23:40 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci physics.ins-det
keywords low-energy electron diffractionLEED I(E)reliability factorR factorsurface structure determinationPendry R_Poptimizationelectron diffraction
0
0 comments X

The pith

A modified R factor R_S replaces Pendry's R_P for LEED while fixing noise, offset sensitivity, and false-zero problems.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Quantitative low-energy electron diffraction determines surface atomic positions by optimizing the match between experimental and calculated I(E) intensity curves. The standard Pendry R_P factor suffers from noise during optimization, extreme sensitivity to tiny intensity offsets, and the ability to reach zero for qualitatively dissimilar curves. The paper introduces R_S as a direct replacement that removes these defects. In tests with imperfect data, R_S steers the search to the correct structure at least as reliably as R_P and better than the Zanazzi-Jona R factor. A more stable agreement metric supports more trustworthy surface-structure solutions.

Core claim

The central claim is that the new reliability factor R_S can be used as a direct replacement for R_P but avoids its shortcomings as a noisy target function, its high sensitivity to small intensity offsets, and the fact that R_P equals zero can be achieved by qualitatively very different curves; R_S performs as well as or better than R_P in steering optimizations to the correct result under imperfections in experimental data, while R_ZJ performs worse.

What carries the argument

The modified reliability factor R_S, a quantitative measure of agreement between two I(E) curves that reduces noise and intensity-offset sensitivity compared with prior factors.

If this is right

  • R_S can be substituted directly for R_P inside existing LEED analysis codes.
  • Optimizations become more robust to common experimental imperfections when R_S is used.
  • R_ZJ should be avoided relative to R_S or R_P when data quality is imperfect.
  • A zero value of R_S corresponds more reliably to qualitatively similar curves than a zero value of R_P.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Routine adoption of R_S could raise the accuracy of published surface structures in materials science.
  • The same style of refinement might benefit agreement metrics used in other electron or X-ray diffraction methods.
  • LEED software developers can implement R_S immediately and validate it on their own benchmark datasets.

Load-bearing premise

The specific imperfections tested are representative of those present in typical real-world LEED data.

What would settle it

Add realistic levels of noise and small intensity offsets to a high-quality experimental I(E) set for a known surface structure, then run optimizations with R_S, R_P, and R_ZJ and check which factor recovers the accepted geometry.

Figures

Figures reproduced from arXiv: 2511.05448 by Alexander M. Imre, Lutz Hammer, Michael Schmid, Michele Riva, Ulrike Diebold.

Figure 1
Figure 1. Figure 1: (a) The Y (L) function of Pendry’s R factor (blue). The gray curve shows an invertible function, YM, which forms the basis for the R factor RS presented in this work. The horizontal axis is the scaled logarithmic derivative L of the intensity [equation (3)]. (b) The I(E) plot shows two curves (black, red) that share the same YP(E) function, shown in the bottom. This implies that Pendry’s R factor between t… view at source ↗
Figure 2
Figure 2. Figure 2: Intensities and various Y functions at a minimum with different intensity offsets. Note that all Y functions are sensitive to relative changes only, but invariant to a scale factor of the intensity. Thus, the red intensity curve in (a) would be equivalent to a curve with a peak-valley height of 1.0 (like the other curves), but an intensity of 2.6 at the minimum. (b) For small, but non-zero intensity offset… view at source ↗
Figure 3
Figure 3. Figure 3: (a) R factors and (b) their derivatives near the minimum. The horizontal axis is the displacement ∆z of the two symmetry-equivalent uppermost Fe atoms of the α-Fe2O3(1¯102)-(1 × 1) structure [9]. The vertical axis for the Zanazzi–Jona R factor has been scaled by a factor of 7 since the curvature of RZJ vs. ∆z is lower than that of the two other R factors by roughly this factor. The dashed lines indicate th… view at source ↗
Figure 4
Figure 4. Figure 4: Contour plots of (a) RP and (b) RS as a function of the x and y displacement of the upper two symmetry-equivalent Fe atoms of α-Fe2O3(1¯102)-(1 × 1), based on tensor-LEED calculations. The red circles mark a few points where the gradients of RP are grossly misleading, even far from the minimum. Minima are marked by small, red crosses. 7 [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Experimental I(E) curve and its Y functions. The black curve at the bottom is I(E) of the (1,2) beam of α-Fe2O3(1¯102)-(1 × 1). The upper curves are the Y functions of Pendry’s R factor RP (blue) and of the R factor RS introduced in this work (orange). Blue arrows mark cusps of YP, where the YP(L) function reaches an extremum and folds back (cf. figure 1a). The inset around E = 350 eV shows the Y functions… view at source ↗
Figure 6
Figure 6. Figure 6: Correlation between RS and RP. Each data point gives the R factors between the calculated and experimental I(E) curve of a particular beam; also the overall R factors for all beams are included (marked by black circles). The data are from three systems, α-Fe2O3(1¯102)-(1 × 1) [9] (blue squares), Ir(100)-p(3 × 1)-MnO2 [25] (pink diamonds), and Pt(111)-(10 × 10)-49Te [20] (brown crosses). YS functions are es… view at source ↗
Figure 7
Figure 7. Figure 7: Deviations of the best-fit parameters for data sets with imperfect I(E) curves. The plot shows statistics of the normalized deviations of 30 fit parameters obtained using imperfect experimental data sets, with respect to the parameters derived from the best data set, when using a given R factor for both. It therefore provides a measure for how sensitive structure optimization with a given R factor is with … view at source ↗
read the original abstract

Quantitative low-energy electron diffraction [LEED $I(V)$ or LEED $I(E)$, the evaluation of diffraction intensities $I$ as a function of the electron energy] is a versatile technique for the study of surface structures. The technique is based on optimizing the agreement between experimental and calculated intensities. Today, the most commonly used measure of agreement is Pendry's $R$ factor $R_\mathrm{P}$. While $R_\mathrm{P}$ has many advantages, it also has severe shortcomings, as it is a noisy target function for optimization and very sensitive to small offsets of the intensity. Furthermore, $R_\mathrm{P} = 0$, which is meant to imply perfect agreement between two $I(E)$ curves can also be achieved by qualitatively very different curves. We present a modified $R$ factor $R_\mathrm{S}$, which can be used as a direct replacement for $R_\mathrm{P}$, but avoids these shortcomings. We also demonstrate that $R_\mathrm{S}$ is as good as $R_\mathrm{P}$ or better in steering the optimization to the correct result in the case of imperfections of the experimental data, while another common $R$ factor, $R_\mathrm{ZJ}$ (suggested by Zanazzi and Jona) is worse in this respect.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper proposes a modified reliability factor R_S for quantitative LEED I(V) analysis as a direct replacement for Pendry's R_P. It claims that R_S avoids key shortcomings of R_P, including being a noisy target for optimization, high sensitivity to small intensity offsets, and the possibility that R_P=0 for qualitatively dissimilar I(E) curves. The authors further demonstrate through tests on imperfect experimental data that R_S steers structural optimization to the correct result at least as reliably as R_P and better than the Zanazzi-Jona R_ZJ factor.

Significance. If the performance claims hold under representative conditions, R_S could provide a more stable and reliable metric for LEED structure refinement, potentially improving the robustness of surface crystallography results in materials science. The work directly addresses practical limitations in a widely used technique and includes comparative optimization tests, which is a strength.

major comments (2)
  1. [§4] §4 (Optimization tests under data imperfections): The specific imperfections applied to the experimental I(E) curves (e.g., energy calibration offsets, background artifacts, intensity scaling variations, or instrument-specific noise) are not described in sufficient detail. This is load-bearing for the central claim that R_S performs as well as or better than R_P while R_ZJ is worse, because without knowing the test protocol it is impossible to judge whether the result generalizes to typical LEED data quality.
  2. [Eq. (definition of R_S)] Eq. (definition of R_S): The manuscript should explicitly show that R_S does not reduce to a post-hoc scaling or fitted parameter derived from the same dataset, to confirm the circularity concern raised in the review process is avoided.
minor comments (2)
  1. [Abstract] The abstract states the claims without equations or quantitative metrics; moving a brief definition or key formula to the abstract or early introduction would improve accessibility.
  2. [Figures in §4] Figure captions for the optimization steering plots should include the exact number of test cases and the precise error models used.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough review and constructive comments on our manuscript. We address each major comment below and will revise the manuscript to incorporate clarifications that strengthen the presentation of our results.

read point-by-point responses

  1. Referee: [§4] §4 (Optimization tests under data imperfections): The specific imperfections applied to the experimental I(E) curves (e.g., energy calibration offsets, background artifacts, intensity scaling variations, or instrument-specific noise) are not described in sufficient detail. This is load-bearing for the central claim that R_S performs as well as or better than R_P while R_ZJ is worse, because without knowing the test protocol it is impossible to judge whether the result generalizes to typical LEED data quality.

    Authors: We agree that more explicit detail on the imperfections is warranted to allow readers to assess the generalizability of the optimization tests. In the revised manuscript we will expand the description in Section 4 to specify the exact energy calibration offsets, background levels, intensity scaling factors, and noise characteristics applied to the experimental I(E) curves, together with the rationale for choosing these representative values. revision: yes

  2. Referee: [Eq. (definition of R_S)] Eq. (definition of R_S): The manuscript should explicitly show that R_S does not reduce to a post-hoc scaling or fitted parameter derived from the same dataset, to confirm the circularity concern raised in the review process is avoided.

    Authors: R_S is defined using a fixed smoothing kernel applied to the normalized intensity curves; the kernel width is a constant chosen independently of any particular dataset. No parameters are fitted to the experimental or theoretical curves being compared. In the revised manuscript we will add an explicit demonstration (both algebraically and with numerical examples) that R_S is invariant under uniform intensity scaling and does not involve any post-hoc adjustment derived from the data under comparison, thereby confirming the absence of circularity. revision: yes

Circularity Check

0 steps flagged

No circularity: R_S definition and validation stand independently

full rationale

The paper proposes R_S as a modified reliability factor with an explicit mathematical definition intended as a direct replacement for R_P. This definition does not reduce by construction to any fitted parameter or normalization drawn from the same dataset under analysis. The performance demonstration—that R_S steers optimization to the correct structure at least as reliably as R_P under data imperfections—is an empirical comparison against external test cases and is not required for the definition itself. No self-citation chain, uniqueness theorem imported from prior author work, or ansatz smuggled via citation appears in the abstract or reader's summary. The derivation chain is therefore self-contained: a new functional form plus separate empirical checks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract provides no explicit equations, so no free parameters, axioms, or invented entities can be identified; the new R_S is described only at the level of its intended properties rather than its mathematical construction.

pith-pipeline@v0.9.0 · 5540 in / 1243 out tokens · 37030 ms · 2026-05-17T23:40:19.047922+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

26 extracted references · 26 canonical work pages

  1. [1]

    Pendry J B 1974Low Energy Electron Diffraction: The Theory and its Application to Determination of Surface Structure(London: Academic Press)

    work page

  2. [2]

    ed (Springer Series in Surface Sciencesno 6) (Berlin: Springer)

    Van Hove M A, Weinberg W H and Chan C M 1986Low-Energy Electron Diffraction: Experiment, Theory and Surface Structure Determinationsoftcover reprint of the hardcover 1. ed (Springer Series in Surface Sciencesno 6) (Berlin: Springer)

    work page

  3. [3]

    Van Hove M A, Moritz W, Over H, Rous P J, Wander A, Barbieri A, Materer N, Starke U and Somorjai G A 1993 Automated determination of complex surface structures by LEEDSurf. Sci. Rep.19191–229

    work page 1993

  4. [4]

    Kristallogr.213615–634

    Heinz K and Hammer L 1998 Surface crystallography by low energy electron diffractionZ. Kristallogr.213615–634

    work page 1998

  5. [5]

    Heinz K 2012 Electron based methods: Low-energy electron diffraction (LEED)Surface and Interface Sciencevol 1 ed Wandelt K (Weinheim: Wiley-VCH) pp 93–150

    work page 2012

  6. [6]

    Fauster T, Hammer L, Heinz K and Schneider M A 2020Surface Physics: Fundamentals and Methods(De Gruyter Oldenbourg)

    work page

  7. [7]

    Moritz W and Van Hove M A 2022Surface Structure Determination by LEED and X-rays (Cambridge: Cambridge University Press)

    work page

  8. [8]

    Sci.754122696

    Held G 2025 Structure determination by low-energy electron diffraction–A roadmap to the futureSurf. Sci.754122696

    work page 2025

  9. [9]

    Kraushofer F, Imre A M, Franceschi G, Kißlinger T, Rheinfrank E, Schmid M, Diebold U, Hammer L and Riva M 2025 ViPErLEED package I: Calculation ofI(V) curves and structural optimizationPhys. Rev. Res.7013005

    work page 2025

  10. [10]

    Schmid M, Kraushofer F, Imre A M, Kißlinger T, Hammer L, Diebold U and Riva M 2025 ViPErLEED package II: Spot tracking, extraction, and processing ofI(V) curvesPhys. Rev. Res.7013006 15 Imreet al

    work page 2025

  11. [11]

    D¨ orr F, Schmid M, Kraushofer F, Hammer L, Diebold U and Riva M ViPErLEED package III: Data acquisition for quantitative low-energy electron diffractionin preparation

    work page

  12. [12]

    Imre A M, Haidegger P, Kraushofer F, Wanzenb¨ ock R, Hable T, Tobisch S, Kienzer M, Buchner F, Carrete J, Madsen G K H, Schmid M, Diebold U and Riva M 2025 Structural optimization in tensor LEED using a parameter tree andR-factor gradientssubmitted

    work page 2025

  13. [13]

    Sci.6261–80

    Zanazzi E and Jona F 1977 A reliability factor for surface structure determinations by low-energy electron diffractionSurf. Sci.6261–80

    work page 1977

  14. [14]

    Pendry J B 1980 Reliability factors for LEED calculationsJ. Phys. C: Solid State Phys.13937

    work page 1980

  15. [15]

    Sci.416423–429

    Sporn M, Platzgummer E, Forsthuber S, Schmid M, Hofer W and Varga P 1998 The accuracy of quantitative LEED in determining chemical composition profiles of substitutionally disordered alloys: a case studySurf. Sci.416423–429

    work page 1998

  16. [16]

    Sci.488219–232

    Blum V, Hammer L, Meier W and Heinz K 2001 Quantification of substitutional disorder and atomic vibrations by LEED – the role of parameter correlationsSurf. Sci.488219–232

    work page 2001

  17. [17]

    Rouˇ cka R, Jiruˇ se J andˇSikola T 2002 Spot intensity processing in LEED imagesVacuum65 121–126 ISSN 0042-207X

    work page 2002

  18. [18]

    Mayer A, Salopaasi H, Pussi K and Diehl R D 2012 A novel method for the extraction of intensity–energy spectra from low-energy electron diffraction patternsComputer Physics Communications1831443–1447

    work page 2012

  19. [19]

    Rous P J 1992 The tensor LEED approximation and surface crystallography by low-energy electron diffractionProg. Surf. Sci.393–63

    work page 1992

  20. [20]

    Kißlinger T, Schewski A, Raabgrund A, Loh H, Hammer L and Schneider M A 2023 Surface telluride phases on Pt(111): Reconstructive formation of unusual adsorption sites and well-ordered domain wallsPhys. Rev. B108205412

    work page 2023

  21. [21]

    Intell.42514–526

    Chen S, Montgomery J and Boluf´ e-R¨ ohler A 2015 Measuring the curse of dimensionality and its effects on particle swarm optimization and differential evolutionAppl. Intell.42514–526

    work page 2015

  22. [22]

    Nocedal J and Wright S J 2006Numerical Optimizationsecond edition ed Springer series in operations research and financial engineering (New York, NY: Springer)

    work page

  23. [23]

    Schmid M, Rath D and Diebold U 2022 Why and how Savitzky–Golay filters should be replacedACS Meas. Sci. Au2185–196

    work page 2022

  24. [24]

    Dahlquist G and Bj¨ orck˚A 2008Numerical Methods in Scientific Computing, Volume I (Philadelphia: Society for Industrial and Applied Mathematics)

    work page

  25. [25]

    Ferstl P, Hammer L, Sobel C, Gubo M, Heinz K, Schneider M A, Mittendorfer F and Redinger J 2016 Self-organized growth, structure, and magnetism of monatomic transition-metal oxide chainsPhys. Rev. Lett.117046101

    work page 2016

  26. [26]

    Tear S P, Roll K and Prutton M 1981 A comparison of reliability (R) factors in a LEED structural analysis of the copper (111) surfaceJ. Phys. C: Solid State Phys.143297 16

    work page 1981