pith. sign in

arxiv: 2604.15899 · v1 · submitted 2026-04-17 · ❄️ cond-mat.mtrl-sci

Experimentally-validated multi-slice simulation of electron diffraction patterns

Xinke Xiao , Tianle Ma , Lingxuan Shao (SJTU) , Jun Liu , Qiwei Shi (SJTU) , Canying Cai , St\'ephane Roux (LMPS) This is my paper

Pith reviewed 2026-05-10 09:07 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords electron backscatter diffractionmulti-slice simulationBloch wave methodHR-EBSDAl-Mg alloysdynamical diffractiondefect structurespattern simulation
0
0 comments X

The pith

A fifth-order multi-slice expansion produces experimental EBSD patterns with precision matching the Bloch wave method after distortion correction.

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

The paper shows that optimizing the multi-slice method by dropping the high-energy approximation and using higher-order Taylor expansions yields the MS5 variant, which generates electron backscatter diffraction patterns that align with experimental data from polycrystalline Al-Mg alloys. Adding a tailored isotropic distortion correction and standard stereographic reconstruction brings MS5 accuracy to the level of the established Bloch wave technique. This matters because Bloch wave simulations are limited to perfect crystals, whereas multi-slice can incorporate defect structures, opening routes to more realistic HR-EBSD analysis of strain and dislocations in real materials. The work supplies the first direct experimental benchmark for multi-slice EBSD simulations.

Core claim

By abandoning the high-energy hypothesis and applying higher-order Taylor expansions to the forward-only Schrödinger equation, the fifth-order multi-slice scheme (MS5) achieves a practical balance of speed and fidelity. When augmented by an isotropic distortion correction model and stereographic triangle reconstruction, MS5 simulations reproduce experimental EBSD patterns from Al-Mg polycrystals with accuracy comparable to Bloch wave results, providing the first experimental validation of the multi-slice approach for EBSD and enabling future modeling of crystals containing defects.

What carries the argument

The fifth-order Taylor expansion of the forward-only Schrödinger equation within the multi-slice framework, combined with a tailored isotropic distortion correction model and stereographic triangle reconstruction.

Load-bearing premise

That the higher-order expansion plus isotropic distortion correction will reproduce experimental EBSD patterns for crystals with defects without introducing systematic biases or needing extensive per-sample fitting.

What would settle it

A side-by-side comparison in which MS5-simulated Kikuchi band positions or intensities deviate measurably from experimental patterns obtained from a defect-free reference crystal after applying the correction and reconstruction steps.

read the original abstract

High-Resolution Electron Backscatter Diffraction (HR-EBSD) has advanced rapidly in recent years, significantly improving elastic strain measurements and dislocation density evaluation with submicron spatial resolution. To achieve better accuracy in the measurements, high-quality dynamical simulation patterns are required to be matched with experimental ones. Currently, the most widely used pattern simulation method, the Bloch Wave method (BW), can accurately predict the positions and brightness of Kikuchi poles and bands, but is intrinsically limited to perfect crystal structures. Another simulation scheme, the multi-slice method (MS), follows the evolution of electron waves as they travel through the sample. MS is advantageous in simulating various defect structures with more diffraction details. Yet, it is mainly considered for theoretical developments and has not been compared to experimental data. This paper optimizes the MS method by abandoning the high-energy hypothesis and utilizing higher-order Taylor expansions to approach the forward-only Schrodinger equation. Experimental EBSD patterns of polycrystal Al-Mg alloys are used to challenge MS simulations as a reference for indexation. It is demonstrated that the 5th-order expansion of MS, referred to as MS5, achieves a good balance between computational cost and pattern precision. A tailored isotropic distortion correction model and standard stereographic triangle reconstruction enhance the precision of MS5 to be comparable with BW. To the best of our knowledge, this study provides the first comparison of MS EBSD simulations with experimental data. It opens new possibilities for EBSD characterization, such as reproducing diffraction patterns of crystals with various defects.

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 manuscript optimizes the multi-slice (MS) method for EBSD pattern simulation by dropping the high-energy approximation and using a 5th-order Taylor expansion (MS5) of the forward-only Schrödinger equation. It reports that MS5, combined with a tailored isotropic distortion correction and standard stereographic triangle reconstruction, produces patterns whose precision on experimental Al-Mg polycrystal EBSD data is comparable to Bloch-wave (BW) simulations, and presents this as the first experimental comparison of MS EBSD patterns.

Significance. If the reported comparability is shown to arise from the MS5 wave propagation rather than from post-hoc correction, the work would be significant: it supplies the first experimental benchmark for multi-slice EBSD simulation and opens a route to dynamical pattern simulation in crystals containing defects, which the Bloch-wave method cannot address.

major comments (2)
  1. Abstract: the claim that the tailored isotropic distortion correction 'enhance[s] the precision of MS5 to be comparable with BW' is load-bearing for the central validation result, yet the manuscript does not state whether the correction parameters are obtained from independent calibration or by fitting to the same experimental patterns used for the MS5–BW–experiment comparison. If the latter, the agreement may reflect compensation of systematic discrepancies in the MS5 scattering or propagation rather than validation of the underlying expansion.
  2. Abstract and methods sections: the assertion of 'good balance between computational cost and pattern precision' and comparability with BW is presented without quantitative metrics (e.g., mean angular error, band contrast correlation, or R² values), error bars, sample sizes, or explicit exclusion criteria for the polycrystal patterns. These omissions prevent assessment of whether the comparability is statistically robust or merely qualitative.
minor comments (2)
  1. The manuscript should clarify the precise order of the Taylor expansion retained in MS5 and any truncation-error estimates, preferably with reference to a specific equation.
  2. Figure captions and text should explicitly label which panels show raw MS5 output, corrected MS5 output, BW output, and experimental data to allow direct visual assessment of the correction's effect.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments, which help strengthen the clarity and rigor of our validation. We address each major point below and will revise the manuscript to incorporate the requested clarifications and quantitative details.

read point-by-point responses

  1. Referee: Abstract: the claim that the tailored isotropic distortion correction 'enhance[s] the precision of MS5 to be comparable with BW' is load-bearing for the central validation result, yet the manuscript does not state whether the correction parameters are obtained from independent calibration or by fitting to the same experimental patterns used for the MS5–BW–experiment comparison. If the latter, the agreement may reflect compensation of systematic discrepancies in the MS5 scattering or propagation rather than validation of the underlying expansion.

    Authors: We appreciate the referee's concern about potential circularity in the validation. The isotropic distortion correction parameters were obtained via independent calibration on reference EBSD patterns from a separate standard single-crystal sample (not part of the Al-Mg polycrystal dataset used for MS5–BW–experiment comparison). This calibration follows established procedures for correcting known instrumental distortions in the EBSD detector geometry and is applied uniformly. We will revise the abstract and methods sections to explicitly state the independent nature of this calibration and describe the reference dataset used, thereby confirming that the correction does not compensate for simulation discrepancies. revision: yes

  2. Referee: Abstract and methods sections: the assertion of 'good balance between computational cost and pattern precision' and comparability with BW is presented without quantitative metrics (e.g., mean angular error, band contrast correlation, or R² values), error bars, sample sizes, or explicit exclusion criteria for the polycrystal patterns. These omissions prevent assessment of whether the comparability is statistically robust or merely qualitative.

    Authors: We agree that quantitative metrics are necessary to substantiate the claims of comparability and balance. In the revised manuscript we will add: (i) mean angular deviation between simulated and experimental Kikuchi bands, (ii) band contrast correlation coefficients and R² values for intensity profiles, (iii) the total number of patterns analyzed (from multiple grains across the polycrystal samples), (iv) error bars or standard deviations on the reported metrics, and (v) explicit exclusion criteria (e.g., patterns with band contrast below a defined threshold or insufficient pole visibility were omitted). These additions will be placed in the results and methods sections with accompanying tables or figures to allow statistical evaluation of the MS5 performance relative to BW. revision: yes

Circularity Check

0 steps flagged

No significant circularity; MS5 derivation and experimental validation are independent

full rationale

The MS5 method is obtained by higher-order Taylor expansion of the forward-only Schrödinger equation, a first-principles step whose equations are stated without reference to the Al-Mg EBSD data. Validation consists of direct comparison to independently acquired experimental patterns, with the isotropic distortion correction and stereographic reconstruction described as post-processing enhancements rather than redefinitions of the simulated intensities. No load-bearing self-citation, self-definitional loop, or fitted parameter renamed as prediction appears in the derivation chain. The result is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The approach rests on the forward-only Schrödinger equation and Taylor expansions for wave propagation; the distortion correction model is described as tailored and therefore likely contains adjustable parameters. No new physical entities are introduced.

free parameters (2)
  • Taylor expansion order
    Chosen as 5th order to balance cost and precision; the value is selected after comparison rather than derived from first principles.
  • Parameters of isotropic distortion correction model
    Tailored to match experimental patterns; these are adjustable coefficients whose fitting procedure is not detailed in the abstract.
axioms (2)
  • domain assumption Electron wave evolution in the sample can be approximated by the forward-only Schrödinger equation
    Invoked to justify the multi-slice propagation scheme.
  • domain assumption Higher-order Taylor expansions of the propagator yield sufficient accuracy for EBSD pattern simulation
    Basis for replacing the high-energy hypothesis with the 5th-order MS5 scheme.

pith-pipeline@v0.9.0 · 5605 in / 1503 out tokens · 48466 ms · 2026-05-10T09:07:06.695313+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

41 extracted references · 41 canonical work pages

  1. [1]

    K. Z. Baba-Kishi and D. J. Dingley. Backscatter Kikuchi diffraction in the SEM for identification of crystallographic point groups. Scanning, 11(6):305–312, 1989

    work page 1989

  2. [2]

    H. Bethe. Theorie der beugung von elektronen an kristallen. Annalen der Physik, 392(17):55–129, 1928

    work page 1928

  3. [3]

    Stojakovic

    D. Stojakovic. Electron backscatter diffraction in materials characterization. Processing and Application of Ceramics, 6(1):1–13, 2012

    work page 2012

  4. [4]

    D. J. Dingley and S. I. Wright. Determination of crystal phase from an electron backscatter diffraction pattern. Applied Crystallography, 42(2):234–241, 2009

    work page 2009

  5. [5]

    D. J. Jensen. Applications of orientation mapping by scanning and transmission electron microscopy. Ultramicroscopy, 67(1-4):25–34, 1997

    work page 1997

  6. [6]

    D. Jha, S. Singh, R. Al-Bahrani, W. Liao, A. Choudhary, M. De Graef, and A. Agrawal. Extracting grain orientations from EBSD patterns of polycrystalline materials using convolutional neural networks. Microscopy and Microanalysis, 24(5):497–502, 2018

    work page 2018

  7. [7]

    Calcagnotto, D

    M. Calcagnotto, D. Ponge, E. Demir, and D. Raabe. Orientation gradients and geometrically necessary dislocations in ultrafine grained dual -phase steels studied by 2D and 3D EBSD. Materials Science and Engineering: A, 527(10-11):2738–2746, 2010

    work page 2010

  8. [8]

    A. J. Wilkinson. Measurement of elastic strains and small lattice rotations using electron back scatter diffraction. Ultramicroscopy, 62(4):237–247, 1996

    work page 1996

  9. [9]

    N. C. K Lassen, D. J. Jensen, and K. Conradsen. Image processing procedures for analysis of electron back scattering patterns. Scanning microscopy, 6(1):115–121, 1992

    work page 1992

  10. [10]

    A. J. Wilkinson, G. Meaden, and D. J. Dingley. High-resolution elastic strain measurement from electron backscatter diffraction patterns: New levels of sensitivity. Ultramicroscopy, 106(4):307–313, 2006

    work page 2006

  11. [11]

    Villert, C

    S. Villert, C. Maurice, C. Wyon, and R. Fortunier. Accuracy assessment of elastic strain measurement by EBSD. Journal of microscopy, 233(2):290–301, 2009

    work page 2009

  12. [12]

    T. B. Britton, C. Maurice, R. Fortunier, J. H. Driver, A. P. Day, G. Meaden, D. J. Dingley, K. Mingard, and A. J. Wilkinson. Factors affecting the accuracy of high-resolution electron backscatter diffraction when using simulated patterns. Ultramicroscopy, 110(12):1443–1453, 2010

    work page 2010

  13. [13]

    Winkelmann

    A. Winkelmann. Model mechanisms in Kikuchi pattern formation from crystals. Crystal Research 27 / 29 and Technology, 52(4):1600288, 2017

    work page 2017

  14. [14]

    Winkelmann

    A. Winkelmann. Dynamical simulation of electron backscatter d iffraction patterns. In Electron backscatter diffraction in materials science, pages 21–33. Springer, 2009

    work page 2009

  15. [15]

    Maurice, K

    C. Maurice, K. Dzieciol, and R. Fortunier. A method for accurate localisation of EBSD pattern centres. Ultramicroscopy, 111(2):140–148, 2011

    work page 2011

  16. [16]

    Y. H. Chen, S. U. Park, D. Wei, G. Newstadt, M. A. Jackson, J. P. Simmons, M. De Graef, and A. O. Hero. A dictionary approach to electron backscatter diffraction indexing. Microscopy and Microanalysis, 21(3):739–752, 2015

    work page 2015

  17. [17]

    Winkelmann, C

    A. Winkelmann, C. Trager-Cowan, F. Sweeney, A. Day, and P. Parbrook. Many-beam dynamical simulation of electron backscatter diffraction patterns. Ultramicroscopy, 107(4-5):414–421, 2007

    work page 2007

  18. [18]

    Winkelmann

    A. Winkelmann. Dynamical effects of anisotropic inelastic scattering in electron backsca tter diffraction. Ultramicroscopy, 108(12):1546–1550, 11 2008

    work page 2008

  19. [19]

    Winkelmann, G

    A. Winkelmann, G. Nolze, G. Cios, T. Tokarski, P. Baa, B. Hourahine, and C. Trager-Cowan. Kikuchi pattern simulations of backscattered and transmitted electrons. Journal of Microscopy , 284(2):157–184, 2021

    work page 2021

  20. [20]

    J. M. Cowley and A. F. Moodie. The scattering of electrons by atoms and crystals. I. A new theoretical approach. Acta Crystallographica, 10(10):609–619, 1957

    work page 1957

  21. [21]

    Van Dyck

    D. Van Dyck. The path integral formalism as a new description for the diffraction of high-energy electrons in crystals. physica status solidi (b), 72(1):321–336, 1975

    work page 1975

  22. [22]

    Van Dyck

    D. Van Dyck. Fast computational procedures for the simulation of structure images in complex or disordered crystals: a new approach. Journal of Microscopy, 119(1):141–152, 1980

    work page 1980

  23. [23]

    C. Y . Cai, S. J. Zeng, H. R. Liu, and Q. B. Yang. Computational comparison of the conventional multislice method and the real space multislice method for simulating exit wavefunctions. Micron, 40(3):313–319, 2009

    work page 2009

  24. [24]

    J. H. Chen and D. Van Dyck. Accurate multislice theory for elastic electron scattering in transmission electron microscopy. Ultramicroscopy, 70(1-2):29–44, 1997

    work page 1997

  25. [25]

    J. H. Chen, D. Van Dyck, and M. Op de Beeck. Multislice method for large beam tilt with application to HOLZ effects in triclinic and monoclinic crystals. Foundations of Crystallography , 53(5):576–589, 1997

    work page 1997

  26. [26]

    C. Y . Cai and J. H. Chen. An accurate multislice method for low -energy transmission electron microscopy. Micron, 43(2-3):374–379, 2012. 28 / 29

    work page 2012

  27. [27]

    Q. B. Liu, C. Y . Cai, G. W. Zhou, and Y . G. Wang. Comparison of EBSD patterns simulated by two multislice methods. Journal of Microscopy, 264(1):71–78, 2016

    work page 2016

  28. [28]

    X. Li, X. Li, T. Wu, C. Lv, and C. Cai. EBSD patterns simulation of dislocation structures based on electron diffraction dynamic theory. Micron, 169:103461, 2023

    work page 2023

  29. [29]

    X. Li, X. Li, T. Wu, C. Lv, and C. Cai. EBSD patterns simulating multilayer twins based on dynamical theory of electron diffraction. Journal of Microscopy, 291(3):199–209, 2023

    work page 2023

  30. [30]

    Zhu and M

    C. Zhu and M. De Graef. EBSD pattern simulations for an interaction volume containing lattice defects. Ultramicroscopy, 218:113088, 2020

    work page 2020

  31. [31]

    Q. Shi, D. Loisnard, C. Dan, F. Zhang, H. Zhong, H. Li, Y . Li, Z. Chen, H. Wang, and S. Roux. Calibration of crystal orientation and pattern center of EBSD using integrated digital image correlation. Materials Characterization, 178:111206, 2021

    work page 2021

  32. [32]

    L. M. Peng, G. Ren, S. L. Dudarev, and M. J. Whelan. Robust parameterization of elastic and absorptive electron atomic scattering factors. Foundations of Crystallography, 52(2):257–276, 1996

    work page 1996

  33. [33]

    G. S. Chen, Y . T. He, W. Q. Ming, C. L. Wu, D. Van Dyck, and J. H. Chen. Fast STEM image simulation in low-energy transmission electron microscopy by the accurate Chen-van-Dyck multislice method. Micron, 190:103778, 2025

    work page 2025

  34. [34]

    Y. Yang, Q. Yang, J. Huang, C. Cai, and J. Lin. Quantitative comparison between real space and Bloch wave methods in image simulation. Micron, 100:73–78, 2017

    work page 2017

  35. [35]

    Hild and S

    F. Hild and S. Roux. Digital image correlation, in: Optical Methods for Solid Mechanics. A Full- Field Approach, P. Rastogi and E. Hack, eds., (Wiley-VCH, Weinheim (Germany), 2012) pp. 183-228

    work page 2012

  36. [36]

    Singh, F

    S. Singh, F. Ram, and M. De Graef. EMsoft: open source software for electron diffraction/image simulations. Microscopy and Microanalysis, 23(S1):212–213, 08 2017

    work page 2017

  37. [37]

    Q. Shi, D. Loisnard, Y . Li, Z. Chen, H. Wang, and S. Roux. Accuracy assessment of crystal orientation indexations by EBSD. Measurement Science and Technology, 35(4):045030, 2024

    work page 2024

  38. [38]

    D. Dingley. Progressive steps in the development of electron backscatter diffraction and orientation imaging microscopy. Journal of microscopy, 213(3):214–224, 2004

    work page 2004

  39. [39]

    J. Liu, J. Li, J. Huang, Y . Cui, L. Wang, Q. Shi, Z. Chen, and H. Wang. Mesoscale architectural control of bimodal domains mediates ultrahigh -strength Al com posites with enhanced ductility. International Journal of Plasticity, 200:104670, 2026

    work page 2026

  40. [40]

    Huang, F

    Y . Huang, F. Peng, X. Song, X. Jin, Y . Jiang, W. Zheng, C. Jiang, Z. Wu, Y . Zeng. Optimisation of EBSD indexing through pattern centre calibration and grain bounda ry refinement. Journal of 29 / 29 Microscopy, 300:36-50, 2025

    work page 2025

  41. [41]

    Q. Shi, Y . Zhou, H. Zhong, D. Loisnard, C. Dan, F. Zhang, Z. Chen, H. Wang, S. Roux. Indexation of electron diffraction patterns at grain boundaries. Materials Characterization, 182:111553, 2021

    work page 2021