Nano-Beam Precession Electron Diffraction (N-PED)

• N-PED is an advanced electron diffraction technique that combines a focused nanobeam with beam precession to average dynamical scattering and enhance data quality.
• It is applied for automated crystal orientation mapping, precise strain metrology, phase assignment, and thermal vibration analysis at nanometer resolution.
• Recent advances in detector technology and computational workflows have achieved misorientation resolution below 0.2° and enabled sub-pixel strain mapping.

Nano-Beam Precession Electron Diffraction (N-PED) is an electron diffraction methodology in which a focused nanobeam is precessed about the optic axis, usually while the probe is raster-scanned so that a diffraction pattern is recorded at each probe position. In contemporary practice it is implemented through scanned precession electron diffraction (SPED) and precession-assisted 4D-STEM workflows, where precession is used to suppress dynamical scattering effects, increase the number of excited reflections, and move the recorded patterns toward quasi-kinematic or pseudo-kinematical conditions. Across recent studies, N-PED has been used for automated crystal orientation mapping, phase assignment in multiphase alloys, multidirectional strain mapping in irradiated ceramics, in situ chemomechanical mapping in gas environments, and diffraction-intensity-based thermometry at nanometer resolution (Jeong et al., 2021, MacLaren et al., 2023, Daghbouj et al., 23 Sep 2025, Yang et al., 14 Apr 2025).

1. Physical principle and measurement geometry

The defining feature of N-PED is the combination of a nanometric probe with beam precession. In PED-assisted 4D-STEM, a focused, highly parallel nanobeam is scanned in a 2D grid across the specimen, and a diffraction pattern is recorded at each scan position, producing a 4D dataset with two spatial and two diffraction dimensions (Jeong et al., 2021). During acquisition, the incident beam is precessed at a constant small angle around the optic axis. This conical rocking combines diffraction intensities across an angular range, suppresses dynamical scattering effects, increases the number of excited reflections, and brings the pattern closer to quasi-kinematic conditions (Jeong et al., 2021).

Reported implementations show that the technique is not tied to a single operating regime. Published parameter sets include a beam size of approximately 1 nm with a 0.5° precession angle and step sizes of 2–2.7 nm for orientation mapping (Jeong et al., 2021), a 5–8 nm probe with 8.8 mrad precession and 1.59 nm step size for high-precision orientation analysis in InP nanowires (Corrêa et al., 2023), a spot size of approximately 2.5–3 nm with a 0.5° precession angle and 3 nm step size for mapping an $\alpha/\beta$ titanium alloy (MacLaren et al., 2023), a beam size of approximately 4 nm with a 1° precession angle and 5 nm mapping step for strain mapping in irradiated SiC (Daghbouj et al., 23 Sep 2025), and a nano-probe diameter of approximately 15 nm with 12–25 nm scan steps for in situ gas-solid reaction studies (Sun et al., 26 Apr 2025). For graphene thermometry, a scanning nanobeam around 1.39 nm in diameter was combined with a very small precession angle of 0.02° (Yang et al., 14 Apr 2025).

The principal analytical consequence of precession is not the complete elimination of dynamical effects, but their averaging. A common misconception is that PED produces strictly kinematical diffraction. The literature instead describes quasi-kinematic, pseudo-kinematical, or more uniform diffraction intensities, which are sufficiently regular to support template matching, strain extraction from Bragg-spot positions, and intensity-based refinement workflows (Jeong et al., 2021, Plana-Ruiz et al., 1 Apr 2025). This distinction is important because many of the most accurate implementations still rely on dynamical or two-beam dynamical models rather than purely kinematical ones (Corrêa et al., 2023).

2. Instrumentation, detectors, and acquisition environments

N-PED performance is strongly conditioned by detector architecture. A major methodological development has been the transition from conventional fluorescent-screen-plus-CCD systems to axial CMOS and direct electron detectors. In a PED-assisted 4D-STEM comparison, a scintillator-based CMOS detector provided higher dynamic range and spatial resolution, a wider acceptance angle of up to 143 mrad rather than about 88 mrad, greatly reduced artifacts, and signal-to-noise ratio up to approximately 80 times that of the conventional system (Jeong et al., 2021). The same study reported average dark noise values of 6.07 for the CMOS detector and 18.96 for the conventional detector, with the improved detector capturing approximately 1.6 times more reflections (Jeong et al., 2021).

This detector dependence is not incidental. High-angle reflections are weak but highly orientation-sensitive, and several studies identify their reliable detection as essential for accurate indexing and fine misorientation measurement. In the titanium alloy workflow, the use of a Quantum Detectors MerlinEM direct electron detector was described as crucial because its low noise enabled reliable detection of high-angle spots needed to resolve orientation ambiguities (MacLaren et al., 2023). In gas-environment 4D-STEM, integration of a CheeTah M3 hybrid pixel direct electron detector with PED improved the quantity and quality of detected Bragg peaks and allowed more precise strain measurements (Sun et al., 26 Apr 2025).

N-PED has also been adapted to increasingly complex acquisition environments. In situ gas-solid reaction studies combined a MEMS-based closed-cell TEM holder, a NanoMEGAS DigiStar P2010 precession unit, and a direct electron detector in a temperature-controlled gas environment (Sun et al., 26 Apr 2025). The gas-cell work showed that pressure must be optimized: it must remain high enough to permit reactions but low enough to maximize electron coherence and Bragg peak detection, with values below 300 mbar cited for 100–200 nm samples (Sun et al., 26 Apr 2025). A practical protocol was introduced in which the reaction is paused during 4D-STEM acquisition, typically by rapid temperature quench, so that large-field strain and orientation maps can be recorded without losing temporal alignment to the underlying reaction pathway (Sun et al., 26 Apr 2025).

3. Computational workflows and indexing strategies

The standard N-PED computational pipeline begins with diffraction-pattern acquisition at each scan position and proceeds through indexing, phase assignment, and crystallographic or strain analysis. One open workflow for a two-phase titanium alloy converted TopSpin .app5 files to HDF5 with the fpd Python package, indexed patterns with py4DSTEM against simulated $\alpha$-Ti and $\beta$-Ti templates derived from ICSD CIF files, exported orientation solutions as .ang files, and used orix for quantitative orientation and misorientation mapping (MacLaren et al., 2023). A custom Python script then selected the better-fit phase at each pixel according to the highest correlation and produced a composite orientation map (MacLaren et al., 2023). Variant grouping was carried out with DBSCAN in scikit-learn using a maximum cluster size of 15° in orientation space (MacLaren et al., 2023).

Template matching remains central, but it is no longer the only route to orientation extraction. In InP nanowires, quantitative analysis of diffracted intensities was used instead of conventional library matching: experimental PED intensities were compared with simulations based on the two-beam dynamical diffraction approximation, and the local orientation was obtained by minimizing a residue between experiment and simulation (Corrêa et al., 2023). This yielded a reported angular precision of approximately 0.03 degrees, substantially below the approximately 1 degree associated with traditional ACOM based on pattern matching algorithms (Corrêa et al., 2023).

For strain mapping, reference-based approaches dominate. In irradiated SiC, Topspin and ASTAR’s AutoSTRAIN module performed spot position determination at sub-pixel precision and calculated local strain through cross-correlation between each local diffraction pattern and a reference pattern acquired from an unirradiated, strain-free region (Daghbouj et al., 23 Sep 2025). The general strain measure was reported as

$$\epsilon = \frac{d_{\text{measured}} - d_{\text{ref}}}{d_{\text{ref}}}.$$

Specific reflections were used to separate out-of-plane, in-plane, shear, and oblique strain components (Daghbouj et al., 23 Sep 2025).

Several studies also emphasize that peak localization is itself a nontrivial algorithmic choice. In a small-convergence-angle SPED study on GaAs nanowires, center of mass, cross-correlation, and curve fitting were compared for Bragg-disk localization (Crout et al., 2023). Cross-correlation produced the most accurate and least noisy maps when the diffraction spots had been aggressively and correctly processed, center of mass was simple and fast but could over-smooth the result, and curve fitting produced noisy maps in practice (Crout et al., 2023). This work further stressed the importance of manual inspection, use of spot pairs at $\pm g$, and tuning of intensity thresholds (Crout et al., 2023).

Where weak minority-phase signals overlap with strong matrix reflections, multi-stage indexing becomes necessary. For irradiation-induced G-phase precipitates in austenite, a two-stage template matching scheme indexed the matrix first, masked or subtracted matrix reflections, and then indexed the precipitates from the residual signal (Cautaerts et al., 2021). This approach was enabled by a pixelated detector with sufficient dynamic range to record the weak precipitate reflections simultaneously with intense matrix spots (Cautaerts et al., 2021).

4. Quantitative observables: orientation, strain, and thermal vibration

Orientation mapping is one of the mature outputs of N-PED. In a detector-focused PED-assisted 4D-STEM study, orientation maps reconstructed by template matching showed reduced false indexing and higher reliability when acquired with a scintillator-based CMOS detector (Jeong et al., 2021). After masking reflections close to the direct beam and increasing the template library from 1,326 to 11,476 templates, the reported misorientation resolution reached 0.20° precision and 0.27° accuracy, with sub-degree angular resolution comparable to orientation mapping using Kikuchi diffraction patterns (Jeong et al., 2021). Fine structural features such as nanograins, nanotwins, and sub-grain boundaries were resolved (Jeong et al., 2021).

N-PED also supports quantitative crystallographic analysis in multiphase systems. In the $\alpha/\beta$ titanium alloy study, distinct $\alpha$-Ti laths down to below 10 nm width and the $\beta$-Ti matrix were clearly resolved at a 3 nm step size, with 10,817 pixels indexed as $\alpha$-Ti and 8,283 as $\beta$-Ti after cropping for quality (MacLaren et al., 2023). The resulting orientation relationships were consistent with the Burgers orientation relationship,

$\alpha$0

and variant misorientations near 60–66° about $\alpha$1 were reported (MacLaren et al., 2023).

Strain metrology is the second major N-PED output. In irradiated single-crystal 4H-SiC, N-PED measured a peak out-of-plane strain of approximately 1.0% at about 588 nm depth, while multidirectional maps revealed tensile out-of-plane strain, compressive in-plane strain, and localized shear up to 0.6% (Daghbouj et al., 23 Sep 2025). In polycrystalline $\alpha$2-SiC, it detected localized tensile strain amplification at grain boundaries up to approximately 2.5% together with strain relief in bubble-depleted zones adjacent to the boundaries (Daghbouj et al., 23 Sep 2025). In gas-cell 4D-STEM, strain tensors were expressed in terms of $\alpha$3, $\alpha$4, and $\alpha$5, and the hydrostatic in-plane strain was defined as

$\alpha$6

This framework was used for quantitative in situ chemomechanical mapping during zirconium oxidation (Sun et al., 26 Apr 2025).

A newer application is nanoscale thermometry through thermal vibration analysis. In graphene, scanning PED was used to obtain local kinematic diffraction intensities from which the Debye–Waller factor $\alpha$7 was extracted by a structure-factor-based correction and Wilson-plot-type linear fitting (Yang et al., 14 Apr 2025). The reported kinematical form was

$\alpha$8

The study reported precision of $\alpha$9 for the Debye–Waller factor as a function of temperature and further showed that the factor depends on graphene thickness, with bilayer and multilayer values approximately 2.6 times and 6.2 times the monolayer value, respectively (Yang et al., 14 Apr 2025).

5. Representative materials systems and scientific use cases

A prominent N-PED use case is phase-resolved mapping in engineering alloys. In TIMETAL 550, a two-phase $\beta$0 titanium alloy prepared by standard FIB liftout, N-PED resolved nanometric $\beta$1 laths, recovered expected Burgers orientation relationships, and identified eight main $\beta$2 orientation clusters that could be rationalized into five distinct orientations when mirror relations were taken into account (MacLaren et al., 2023). The study presented the workflow as a method for nanocrystalline two- or multi-phase materials and explicitly linked it to future analysis of microstructural evolution under thermomechanical treatment (MacLaren et al., 2023).

In radiation-damage research, N-PED has been used to separate bulk-average and grain-boundary-localized strain physics. For ion-irradiated SiC, N-PED strain profiles in single-crystal 4H-SiC quantitatively matched both the depth and amplitude of HR-XRD-derived strain from dynamical simulation, whereas XRD failed to detect distinct peak shifts in polycrystalline $\beta$3-SiC because of grain orientation averaging and depth convolution (Daghbouj et al., 23 Sep 2025). The N-PED maps directly visualized strain amplification at grain boundaries and strain relief in bubble-depleted zones, supporting the interpretation that grain boundaries act as active sinks for irradiation-induced defects (Daghbouj et al., 23 Sep 2025).

In situ gas-environment studies extend N-PED from static characterization to chemomechanical evolution. During initial oxidation of zirconium, precession-assisted 4D-STEM in a MEMS gas cell quantified tensile hydrostatic strain in the $\beta$4-Zr substrate increasing from 0 to 2.58% after 100 s oxidation (Sun et al., 26 Apr 2025). The same work emphasized that large fields of view on the order of microns can be combined with nanometer-level probe sizes, provided that reaction timing, gas pressure, and data-acquisition windows are carefully coordinated (Sun et al., 26 Apr 2025).

Precipitation and nanoscale orientation-relationship studies form another established domain. In austenitic steel, scanning nano-beam electron diffraction with a pixelated detector resolved weak reflections from 10–20 nm irradiation-induced G-phase precipitates and, using two-stage template matching, showed that they adopt orientation relationships characteristic of FCC-to-BCC transformations rather than a random relation to the matrix (Cautaerts et al., 2021). Pole-figure analysis was discussed in terms of Kurdjumov–Sachs-, Nishiyama–Wasserman-, Pitsch-, and Bain-type variants (Cautaerts et al., 2021).

A broader implication is that N-PED increasingly acts as a general nanoscale diffraction platform rather than a single-purpose orientation-mapping tool. The range of demonstrated outputs now includes phase mapping, variant analysis, strain tensor mapping, orientation precision below 0.1°, and local thermometry (Corrêa et al., 2023, Yang et al., 14 Apr 2025).

6. Limitations, artifacts, and methodological directions

Several recurring limitations structure current N-PED practice. One is symmetry-induced indexing ambiguity. In the titanium alloy study, a local 180° ambiguity arose because part of the scan was aligned too close to a zone axis of $\beta$5-Ti with two-fold projection symmetry in the zero-order Laue zone, even though the full 3D crystal symmetry does not support that equivalence (MacLaren et al., 2023). The explicit recommendation was to avoid acquisition near such high-symmetry zone axes, or alternatively to reduce camera length or use higher-angle diffraction spots if sufficient signal-to-noise is available (MacLaren et al., 2023).

A second limitation is probe wandering. Precession improves reciprocal-space sampling but can reduce spatial resolution because the illuminated probe follows a circular motion across the specimen. A correction based on precession path segmentation recorded multiple frames during each precession cycle and then aligned the reconstructed virtual images in post-processing (Nordahl et al., 2022). In the reported test cases, edge steepness in virtual bright-field images increased from 0.15 to 0.38 after correction, but the segmented dataset became $\beta$6 times larger and scan times increased by a factor of 2–3 (Nordahl et al., 2022).

Strain mapping introduces its own constraints. TEM lamella thinning can relax strain and thereby reduce the measured magnitude relative to bulk values, electron-beam exposure can induce local artifacts, the probe depth remains limited by TEM geometry, and careful reference-region selection is essential because errors in the reference propagate directly into the strain field (Daghbouj et al., 23 Sep 2025). In gas-cell measurements, insufficient Bragg reflections, high gas pressure, sample thickness, or dynamical effects can produce invalid pixels, and even when PED is used there remains a finite pseudo-strain contribution from out-of-plane tilt, though multislice simulation in one study found this to be below 0.4% for the low-order peaks used in analysis (Sun et al., 26 Apr 2025).

Current development directions are correspondingly clear. Open-source and notebook-driven workflows have already been demonstrated with py4DSTEM, orix, Pyxem, and related tools (MacLaren et al., 2023, Crout et al., 2023). Deep-learning approaches such as FCU-Net and whole-pattern fitting have been identified as emerging alternatives for diffraction analysis in complex environments (Sun et al., 26 Apr 2025). For three-dimensional strain metrology, a formal tensor-tomography framework has shown computationally that scanning precession electron diffraction data from multiple tilt series can, in principle, support reconstruction of the full non-symmetric strain tensor field, and that incorporating precession may be important because it improves the linearity and robustness of the forward model (Tovey et al., 2020). This suggests that N-PED is evolving from a 2D local-diffraction mapping method into a broader quantitative nanocrystallographic framework spanning orientation, strain, thermal vibration, phase relations, and eventually tensor tomography.