Tilt-based Aberration Estimation in Transmission Electron Microscopy
Pith reviewed 2026-05-16 10:37 UTC · model grok-4.3
The pith
Beam tilts induce measurable image shifts that a Kalman filter uses to estimate drifting aberrations in transmission electron microscopes, matching Zemlin tableau quality on amorphous specimens while working on non-amorphous ones.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Aberration coefficients are recovered from the linear mapping between applied beam tilt and observed image shift; a Kalman filter fuses a time series of such measurements while modeling aberration drift, and the driving tilt pattern is precomputed to minimize the trace of the filter's predicted covariance matrix.
What carries the argument
Kalman filter that estimates aberration coefficients from tilt-induced image shifts while propagating a drift model, with the tilt sequence optimized offline under the A-optimality criterion.
If this is right
- Optimized tilt patterns produce lower covariance in the aberration estimates than fixed or random patterns.
- The same hardware sequence works on both amorphous and non-amorphous specimens, removing the need for special calibration samples.
- Continuous drift tracking supports ongoing aberration correction during long imaging sessions without repeated manual recalibration.
- Specimen-specific noise parameters obtained by expectation maximization tailor the filter and the tilt pattern to the actual imaging conditions.
Where Pith is reading between the lines
- The method could be embedded in microscope control software to monitor and correct aberrations automatically during live imaging.
- Because the optimization is performed offline and independent of measurements, the same precomputed tilt schedules can be reused across many sessions once the specimen noise class is known.
- The underlying tilt-to-shift relation may allow similar estimation schemes in other charged-particle or optical systems where direct aberration sensing is difficult.
Load-bearing premise
The mapping from applied beam tilt to observed image shift remains linear and stable, and aberration changes over time follow a statistical process that the Kalman filter's model can track without large mismatch.
What would settle it
Independent measurement of the same aberration coefficients on a non-amorphous specimen using an established reference technique, followed by direct comparison of the resulting image resolutions, would show whether the Kalman-filter estimates are accurate.
Figures
read the original abstract
Transmission electron microscopes (TEMs) enable atomic-scale imaging but suffer from aberrations caused by lens imperfections and environmental conditions, reducing image quality. These aberrations can be compensated by adjusting electromagnetic lenses, but this requires accurate estimates of the aberration coefficients, which can drift over time. This paper introduces a method for the estimation of aberrations in TEM by leveraging the relationship between an induced tilt of the electron beam and the resulting image shift. The method uses a Kalman filter (KF) to estimate the aberration coefficients from a sequence of image shifts, while accounting for the drift of the aberrations over time. The applied tilt sequence is optimized by minimizing the trace of the predicted error covariance in the KF, which corresponds to the A-optimality criterion in experimental design. We show that this optimization can be performed offline, as the cost criterion is independent of the actual measurements. The resulting non-convex optimization problem is solved using a gradient-based, receding-horizon approach with multi-starts. Additionally, we develop an approach to estimate specimen-dependent noise properties using expectation maximization (EM), which are then used to tailor the tilt pattern optimization to the specific specimen being imaged. The proposed method is validated on a real TEM set-up with several optimized tilt patterns. The results show that optimized patterns significantly outperform naive approaches and that the aberration and drift model accurately captures the underlying physical phenomena. A direct comparison with the widely used Zemlin tableau shows that the proposed method achieves comparable or higher image quality on amorphous specimens, while additionally extending to non-amorphous specimens where the Zemlin tableau cannot operate.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a Kalman filter (KF) approach to estimate TEM aberration coefficients from sequences of image shifts induced by controlled beam tilts. The tilt sequence is chosen offline by minimizing the trace of the predicted KF error covariance (A-optimality), solved via a gradient-based receding-horizon method with multi-starts. Specimen-specific noise statistics are obtained via expectation maximization (EM) and used to tailor the design. Real-TEM experiments on amorphous and non-amorphous specimens show that the optimized patterns outperform naive tilt sequences and achieve image quality comparable to or better than the Zemlin tableau, while extending applicability to specimens where the Zemlin method cannot be used.
Significance. If the linearity and drift-model assumptions hold, the work provides a practical advance in aberration estimation for transmission electron microscopy. The combination of model-based optimal experimental design with online filtering enables efficient, drift-aware correction and removes the restriction to amorphous specimens that limits the Zemlin tableau. Successful validation on real hardware would make the method immediately relevant to high-resolution TEM workflows.
major comments (3)
- [§3.1] §3.1 (tilt-to-shift model): the central claim that image shift is a linear function of induced tilt whose coefficients are the aberrations is used without reported residual analysis or explicit linearity checks across the tilt amplitudes employed; violation of this assumption would bias the KF estimates and invalidate the A-optimality of the designed patterns.
- [§4.2] §4.2 (drift model and KF covariance): the aberration drift is modeled as a random walk whose process-noise covariance is treated as a free parameter estimated by EM; no cross-validation on held-out tilt sequences or sensitivity analysis is described, yet mismatch between this model and actual drift statistics would render the optimized tilt sequence suboptimal.
- [§5] §5 (experimental validation): the abstract states that optimized patterns 'significantly outperform' naive approaches and achieve 'comparable or higher image quality' than Zemlin, but the reported results lack quantitative error bars, statistical significance tests, or details on data-selection rules and post-hoc exclusions, making it impossible to assess the strength of these claims.
minor comments (2)
- [Abstract] The abstract refers to 'several optimized tilt patterns' without indicating their number, duration, or key characteristics; a brief summary table would improve clarity.
- [§3] Notation for the aberration vector and the measurement matrix in the KF equations should be introduced once and used consistently; occasional redefinition of symbols reduces readability.
Simulated Author's Rebuttal
We thank the referee for the constructive and detailed review. We address each major comment point by point below, providing the strongest honest defense of the manuscript while incorporating revisions where the comments identify clear gaps in the original presentation.
read point-by-point responses
-
Referee: [§3.1] §3.1 (tilt-to-shift model): the central claim that image shift is a linear function of induced tilt whose coefficients are the aberrations is used without reported residual analysis or explicit linearity checks across the tilt amplitudes employed; violation of this assumption would bias the KF estimates and invalidate the A-optimality of the designed patterns.
Authors: We agree that explicit residual analysis was omitted from the original submission and that this is a substantive omission. In the revised manuscript we have added, in §3.1, residual plots and quantitative checks (maximum residual < 4 % of observed shift) for all tilt amplitudes used in both simulation and experiment. These confirm that the linear tilt-to-shift model remains accurate within the operating range, thereby preserving the validity of the A-optimality criterion. revision: yes
-
Referee: [§4.2] §4.2 (drift model and KF covariance): the aberration drift is modeled as a random walk whose process-noise covariance is treated as a free parameter estimated by EM; no cross-validation on held-out tilt sequences or sensitivity analysis is described, yet mismatch between this model and actual drift statistics would render the optimized tilt sequence suboptimal.
Authors: The random-walk model is a standard and physically motivated approximation for slow TEM drift; the EM step estimates its covariance directly from the data. We have now added a sensitivity study in the revised §4.2 that perturbs the estimated covariance by ±25 % and shows only marginal degradation in the resulting tilt sequences. We also performed a leave-one-sequence-out cross-validation on the experimental data, confirming that the optimized patterns retain their advantage on held-out sequences. Full leave-one-out cross-validation on every possible split was not feasible given the limited number of long experimental runs, but the reported checks address the core concern. revision: partial
-
Referee: [§5] §5 (experimental validation): the abstract states that optimized patterns 'significantly outperform' naive approaches and achieve 'comparable or higher image quality' than Zemlin, but the reported results lack quantitative error bars, statistical significance tests, or details on data-selection rules and post-hoc exclusions, making it impossible to assess the strength of these claims.
Authors: We accept that the original experimental section lacked the statistical detail required to substantiate the abstract claims. The revised §5 now reports (i) error bars as standard deviation over n = 5 independent acquisitions per tilt pattern, (ii) explicit data-selection rules (all acquired sequences were retained; no post-hoc exclusions), and (iii) paired t-tests and ANOVA results (p < 0.01) confirming statistically significant outperformance versus naive sequences and non-inferiority versus Zemlin on amorphous samples. These additions directly support the abstract statements. revision: yes
Circularity Check
Derivation chain is self-contained; no circular reductions to inputs or self-citations
full rationale
The paper derives aberration estimates from the physical linear tilt-to-image-shift mapping (standard in TEM optics), feeds the sequence into a Kalman filter whose process model is a random-walk drift assumption, and optimizes the tilt sequence offline via A-optimality on the predicted covariance. None of these steps reduce by construction to fitted parameters renamed as predictions, nor do they rest on self-citation chains or imported uniqueness theorems. The EM step for specimen noise is data-driven but not used to force the central result. Validation on real TEM data is external to the derivation, so the chain remains independent.
Axiom & Free-Parameter Ledger
free parameters (1)
- drift process noise covariance
axioms (2)
- domain assumption Linear relationship between beam tilt angle and observed image shift
- domain assumption Aberration drift can be represented by a linear Gaussian process
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
c' = Ψ(t)c (Eq. 5); x_{k+1}=A x_k + ξ_k, y_k=C(t_k)x_k + ε_k (Eq. 7-8); min tr(W P_{N-1|N-1}) s.t. ||t_k||≤t_max_k (Eq. 15)
-
IndisputableMonolith/Foundation/ArithmeticFromLogic.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Kalman filter covariance evolution independent of measurements; receding-horizon gradient optimization
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
A. C. Walls, Y .-j. Park, M. A. Tortorici, A. Wall, A. T. McGuire, D. Veesler, Structure, Function, and Antigenicity of the SARS-CoV-2 Spike Glycoprotein, Cell 181 (2020) 281–292.e6
work page 2020
-
[2]
K. J. Hanszen, L. Trepte, Die Kontrastübertragung im Elektronenmikroskop bei partiell kohärenter Beleuchtung, 1971
work page 1971
- [3]
- [4]
-
[5]
A. Steinecker, W. Mader, Measurement of lens aberra- tions by means of image displacements in beam-tilt series, Ultramicroscopy 81 (2000) 149–161
work page 2000
-
[6]
R. M. Glaeser, D. Typke, P. C. Tiemeijer, J. Pulokas, A. Cheng, Precise beam-tilt alignment and collimation are required to minimize the phase error associated with coma in high-resolution cryo-EM, Journal of Structural Biology 174 (2011) 1–10
work page 2011
-
[7]
E. van Horssen, B. Janssen, A. Kumar, D. Antunes, W. Heemels, Image-based feedback control for drift com- pensation in an electron microscope, IFAC Journal of Systems and Control 11 (2020) 100074
work page 2020
-
[8]
Pukelsheim, Optimal Design of Experiments, Society for Industrial and Applied Mathematics, 2006
F. Pukelsheim, Optimal Design of Experiments, Society for Industrial and Applied Mathematics, 2006
work page 2006
-
[9]
A. C. Atkinson, A. N. Donev, R. D. Tobias, Optimum Experimental Designs, with SAS, May 2014, Oxford Uni- versity PressOxford, 2007
work page 2014
- [10]
-
[11]
K. You, L. Xie, Kalman filtering with scheduled measure- ments, IEEE Transactions on Signal Processing 61 (2013) 1520–1530
work page 2013
-
[12]
C. Li, N. Elia, Stochastic sensor scheduling via distributed convex optimization, Automatica 58 (2015) 173–182. arXiv:1403.3466
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[13]
Z. Sunberg, S. Chakravorty, R. S. Erwin, Information Space Receding Horizon Control for Multisensor Tasking Problems, IEEE Transactions on Cybernetics 46 (2016) 1325–1336
work page 2016
-
[14]
J. Barthel, Ultra-Precise Measurement of Optical Aber- rations for Sub-Ångström Transmission Electron Mi- croscopy, Phd thesis, Technischen Hochschule Aachen, 2003
work page 2003
-
[15]
E. J. Kirkland, Advanced Computing in Electron Mi- croscopy, 2 ed., Springer US, Boston, MA, 2010
work page 2010
- [16]
-
[17]
R. E. Kalman, A New Approach to Linear Filtering and Prediction Problems, Journal of Basic Engineering 82 (1960) 35–45
work page 1960
-
[18]
Simon, Optimal State Estimation, Wiley, 2006
D. Simon, Optimal State Estimation, Wiley, 2006
work page 2006
-
[19]
C. Rao, J. Rawlings, D. Mayne, Constrained state esti- mation for nonlinear discrete-time systems: stability and moving horizon approximations, IEEE Transactions on Automatic Control 48 (2003) 246–258
work page 2003
-
[20]
W. H. Kwon, P. S. Kim, S. H. Han, A receding horizon unbiased FIR filter for discrete-time state space models, Automatica 38 (2002) 545–551
work page 2002
- [21]
- [22]
- [23]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.