Pith. sign in

REVIEW 3 major objections 5 minor 8 references

General approaches for shear-correcting coordinate transformations in Bragg coherent diffraction imaging: Part 1

T0 review · 3 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Bragg coherent diffraction images carry a predictable shear that one linear-algebra identity removes.

desk verdict The math holds and the paper is honest about its novelty, but the SEM-based validation is circular and needs to be fixed before publication. read the letter →

arxiv 1909.05353 v1 pith:HTGQYRCV submitted 2019-08-15 physics.ins-det

Siddharth Maddali , Peng Li , Anastasios Pateras , Daniel Timbie , Nazar Delegan , Alex Crook , Hope Lee , Irene Calvo-Almazan
show 7 more authors
Dina Sheyfer Wonsuk Cha F. Joseph Heremans David D. Awschalom Virginie Chamard Marc Allain Stephan O.Hruszkewycz
This is my paper · ORCID
classification physics.ins-det
keywords Braggcoherentdiffractionimagingphaseretrievalshearcorrectioncoordinatetransformationscatteringgeometryptychographyconjugatespacesnanocrystalstrain
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper establishes that the three-dimensional shear seen in Bragg coherent diffraction imaging (BCDI) reconstructions is a predictable geometric consequence of how Fourier space is sampled, not a flaw of phase retrieval, and can be removed by a general post-processing transformation. The argument reduces to a matrix identity linking the Fourier-space sampling basis to the real-space sampling basis of the reconstructed object, so the goniometer settings alone determine how to render an undistorted image. This matters because BCDI encodes a single component of the lattice strain tensor, a frame-dependent quantity, so displaying the recovered complex field on the wrong grid can corrupt both morphology and strain readings. The paper generalizes earlier facility-specific prescriptions into a frame-agnostic relation and demonstrates it on a silicon carbide nanoparticle checked against scanning electron microscopy.

What carries the argument

The load-bearing object is the identity $B_{\mathrm{real}} = B_{\mathrm{recip}}^{-T} D$ (Eq. 12), the discrete analogue of the continuous reciprocal-basis relation $B_q = B_r^{-T}$. $B_{\mathrm{recip}}$ is the matrix whose columns are the non-orthogonal Fourier-space step vectors $\mathbf{q}_i, \mathbf{q}_j, \mathbf{q}_k$ fixed by pixel size, object-detector distance, wavelength, Bragg angles, and rocking increment; $D = \mathrm{diag}(1/N_1, 1/N_2, 1/N_3)$ accounts for the array size. The identity converts the easily measured Fourier-space geometry into the correct real-space grid and reveals the shear as unavoidable: in Bragg geometry $\mathbf{q}_i \perp \mathbf{q}_j$, while $\mathbf{q}_k$ cannot be perpendicular to both. The paper also supplies explicit laboratory-frame expressions for $\mathbf{q}_i, \mathbf{q}_j, \mathbf{q}_k$ from two detector rotations and the rocking step, together with a Rodrigues rotation formula so the same computation can be adapted to other goniometer arrangements.

What would settle it

Simulate a faceted nanocrystal, sample its noiseless Fourier transform on the non-orthogonal grid defined by the paper's $\mathbf{q}_i, \mathbf{q}_j, \mathbf{q}_k$, reconstruct by phase retrieval, apply Eq. (12), and compare rendered facet angles and edge lengths to the ground truth; any systematic discrepancy beyond pixel or angular resolution would show the conjugacy relation is incomplete.

Watch

Extended reading notes

Core claim

On its own terms, this paper establishes that the sheared appearance of a phase-retrieved BCDI object is not an arbitrary artifact but the exact discrete consequence of sampling Fourier space on a non-orthogonal grid. The central identity is Eq. (12): $B_{\mathrm{real}} = B_{\mathrm{recip}}^{-T} D$, with $D=\mathrm{diag}(1/N_1, 1/N_2, 1/N_3)$, which follows from requiring the discrete phase $q^T r = il/N_1 + jm/N_2 + kn/N_3$ to hold on every array index. Here $B_{\mathrm{recip}} = [\mathbf{q}_i\ \mathbf{q}_j\ \mathbf{q}_k]$ lists the Fourier-space step vectors fixed by detector pixel pitch, sample-detector distance, wavelength, Bragg angle, and rocking step, so once those are known the correct real-space sampling basis $B_{\mathrm{real}}$ is determined. Rendering the retrieved array on this sheared basis yields the undistorted scatterer, and the paper verifies the prescription on a SiC nanoparticle by matching the corrected rendering to SEM images.

Load-bearing premise

The whole result rests on the assumption that the measured diffraction pattern is a well-sampled, un-aliased digital Fourier transform of the object; if the fringes are too coarse or the detector window clips them, the computed real-space grid is unreliable no matter how the shear is corrected.

Editorial extensions

If this is right

  • After phase retrieval, applying Eq. (12) with the measured goniometer parameters converts any BCDI reconstruction into a physically accurate, shear-corrected 3D image without modifying the reconstruction algorithm.
  • The corrected grid fixes the scatterer's orientation in the laboratory frame, allowing morphological features to be compared directly with scanning electron micrographs.
  • The linear dependence of the columns of $B_{\mathrm{recip}}$ serves as a degeneracy test: a rocking axis in the scattering plane produces parallel sampling vectors and cannot yield a 3D BCDI dataset.
  • The strain component encoded in BCDI can be computed on the sheared grid by finite differences of the phase along the columns of $B_{\mathrm{real}}$, avoiding interpolation to an orthogonal grid.
  • The same conjugacy relation provides the foundation for Part II's in-algorithm modification of the discrete Fourier transform, enabling reconstruction directly on an orthogonal grid with physical constraints.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The identity is the finite-array analogue of converting a unit cell to its reciprocal cell; if the Fourier aperture is marginal, an aperture-weighted version of Eq. (12) would also predict small shape errors at the object boundary, a testable refinement.
  • A sensitivity analysis of Eq. (12) would quantify which experimental parameter dominates the distortion, something the paper does not do; this could guide the design of new BCDI beamlines.
  • Since the derivation only assumes a fixed Fourier sampling matrix, the shear correction should transfer unchanged to Bragg ptychography and non-standard rocking trajectories; a synthetic-data test with an arbitrary rocking axis would confirm the transfer.
  • If phase retrieval enforces orthogonal-grid real-space support, the post-processing shear correction can only be exact away from constraint boundaries, which is a reason to prefer Part II's in-algorithm approach in constrained reconstructions.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. This paper presents a general derivation of the coordinate transformation needed to correct shear distortion in three-dimensional Bragg coherent diffraction imaging (BCDI) reconstructions. The central result is Eq. (12), which states that the real-space sampling basis B_real is related to the Fourier-space sampling basis B_recip by B_real = B_recip^{-T} D, where D = diag(1/N1, 1/N2, 1/N3). The derivation starts from the continuous Fourier relationship between conjugate bases and adapts it to the discrete Fourier transform, assuming adequate oversampling and no cyclic aliasing. Section 3 develops explicit expressions for the Fourier-space sampling vectors q_i, q_j, q_k for a specific goniometer geometry, and Section 4 applies the correction to a SiC nanocrystal reconstructed from BCDI data, comparing the shear-corrected object with an SEM image.

Significance. If the central derivation is correct, this paper provides a useful, parameter-free unification of previously scattered prescriptions for BCDI shear correction. The derivation of Eq. (12) is clean, self-contained, and appears mathematically sound, and the paper is careful to note the standard prerequisite of adequate fringe sampling. The main weakness is the experimental demonstration: the SEM image is used both to disambiguate the twin-image degeneracy of phase retrieval and then cited as corroboration of the shear-corrected morphology, making the validation circular. The comparison is also qualitative and against a nominally identical but different particle. Because the theoretical contribution is strong, the circularity is a validation issue rather than a flaw in the derivation, but it prevents the paper from claiming an independent experimental confirmation.

major comments (3)
  1. [Section 4 and Fig. 4 caption] The experimental validation is circular. The caption of Fig. 4 explicitly states that the twin-image degeneracy (psi(r) versus psi*(-r)) was resolved by choosing the solution that best matches the SEM image after the shear correction. The text then uses the agreement with that same SEM image as evidence that the correction is valid. This does not independently test Eq. (12), because the SEM image has already influenced which phase-retrieval solution is presented. The paper should either use an independent criterion for twin selection (for example, known support constraints or a second particle not used for selection) or explicitly reframe the example as illustrative rather than confirmatory. In addition, the comparison is purely visual; a quantitative metric (for instance, measured facet angles or dimensions with uncertainties) would materially strengthen the demonstration.
  2. [Section 4, Fig. 4] The SEM comparison is made against a different particle from the same batch, not the particle that was imaged by BCDI. The text notes that the SEM image is of the batch of pillars prior to release, and the BCDI particle was one of many nominally identical pillars. This introduces particle-to-particle variability into the comparison, which is not discussed. The paper should state this limitation explicitly and, if possible, compare the shear-corrected reconstruction to an SEM image of the same particle or otherwise quantify the expected variability.
  3. [Section 5, Summary] The summary states that the shear correction was 'demonstrated' and 'corroborated with SEM images.' Given the circularity described above, this overstates the evidence. The paper's lasting contribution is the theoretical derivation and the general computational prescription, not the independent experimental validation. The summary and abstract should be revised to separate the theoretical result from the illustrative example, and the claim of experimental corroboration should be either removed or qualified.
minor comments (5)
  1. [Abstract] There are typographical errors in the abstract: 'realiable' should be 'reliable' and 'communitcy' should be 'community.'
  2. [Section 2, around Eq. (8)] The phrase 'orthonormal bases defined in this manner are self-congujate' contains a typo: 'self-congujate' should be 'self-conjugate.'
  3. [Section 4, first paragraph] The phrase 'With the the theoretical' contains a duplicated 'the.'
  4. [Section 2, Eq. (11)] The derivation assumes the phase relation in Eq. (11) holds exactly for all integer indices, which requires sufficient oversampling and no cyclic aliasing. The paper states this, but it would be helpful to state the oversampling condition quantitatively or to estimate the error when the condition is only approximately satisfied.
  5. [Section 3.3, Eqs. (14)-(16)] The derivation of the columns of B_real from Eq. (12) is terse; a brief intermediate step showing the reciprocal-lattice-like relation would improve readability for readers not familiar with the analogous Bravais-lattice conversion.

Circularity Check

1 steps flagged · score 2.0 of 10

Derivation of the shear correction is self-contained; only the SEM validation is mildly circular because the same SEM image selects the twin-image solution and then serves as corroboration.

  1. other [Fig. 4 caption (Section 4)]
    "Here, the inherent mathematical degeneracy in the phase retrieval inverse problem (i.e. if ψ(r) is a real-space solution for an observed BCDI diffraction pattern, then so is ψ∗(−r)) was resolved by choosing the solution that most closely reproduced the asymmetric morphological features of the SiC particle in the SEM image, after application of the shear correction."

    The SEM image is first used to disambiguate the twin-image degeneracy by selecting the phase-retrieval solution that best matches it after the shear correction, and then the same SEM image is cited in Section 4 as corroboration: 'The essential morphological features in the SEM image are seen to be reproduced faithfully with the appropriate shear correction.' The agreement is therefore partly by construction: of the two candidate solutions, the one that already resembles the SEM image is chosen, and that choice is then used as evidence. This weakens the experimental demonstration but does not affect the derivation of Eq. (12), which is obtained algebraically from the DFT phase condition with no fitted parameters.

full rationale

The central mathematical derivation is self-contained and not circular. Starting from the continuous Fourier transform, the paper derives the conjugate-basis relation Bq = Br^{-T} in Eq. (8) via a change of variables. It then assumes the discrete phase condition in Eq. (11), q^T r = il/N1 + jm/N2 + kn/N3, which is the standard compatibility condition for a sampled BCDI measurement to be representable as a DFT on an N1 by N2 by N3 grid, and algebraically obtains Breal = Brecip^{-T} D in Eq. (12). No empirical parameter is fitted, and the cited prior prescriptions (Pfeifer, Pateras, Berenguer, Yang) are noted only as equivalent formulations after the derivation is completed, so the self-citations are not load-bearing. The only circularity is in the validation: the SEM image is used once to resolve the intrinsic psi(r) versus psi*(-r) phase-retrieval degeneracy by selecting the solution that matches the SEM, and it is then used again as the corroborating reference for the shear-corrected morphology. That makes the experimental agreement partially self-referential, but it does not reduce the derivation itself to its inputs. Overall the paper's core first-principles result stands independently, with a minor validation caveat.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

No free parameters are fitted: all inputs are measured experimental quantities (wavelength, pixel size, distance, angles, rocking step). The derivation rests on standard Fourier analysis plus domain assumptions about BCDI sampling: sufficient oversampling, parallel slices, small rotation steps, and removable origin offsets. No new physical entities are introduced.

assumptions (5)
  • standard math Continuous and discrete Fourier transform conventions with phase factor exp(-i 2 pi q^T r), giving conjugate bases Bq = Br^{-T}.
    Used to derive Eq. (8) and its discrete analog Eq. (12); this is textbook Fourier analysis, not a domain-specific postulate.
  • domain assumption Discrete phase approximation q^T r = il/N1 + jm/N2 + kn/N3 (Eq. 11) holds for all integer index pairs.
    Bridges the continuous Fourier transform to the DFT used by phase retrieval; requires adequate oversampling and no cyclic aliasing, as the paper itself states before Eq. (11).
  • domain assumption Detector slices of the 3D diffraction pattern can be treated as parallel planes in Fourier space.
    Invoked in Section 3.1 to model the 3D data array as stacked 2D detector images; valid for small rocking angles typical of BCDI.
  • domain assumption The rocking step is small enough that qk can be linearized to first order in Delta theta (Eq. 31).
    Used to compute the third Fourier sampling vector; the full expression Eq. (30) is also provided, so this is an approximation rather than a required exact assumption.
  • domain assumption Origin offsets q0 and r0 can be set to zero without loss of generality.
    Section 2 shows the corresponding phase terms do not affect measured intensity, enabling direct application of Eq. (12).

how reviews work

0 comments
Cite this review

Pith. "Pith review of General approaches for shear-correcting coordinate transformations in Bragg coherent diffraction imaging: Part 1." pith.science (2026). https://pith.science/paper/HTGQYRCV

@misc{pith2026190905353,
  author       = {Pith},
  title        = {Pith review of: General approaches for shear-correcting coordinate transformations in Bragg coherent diffraction imaging: Part 1},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/HTGQYRCV}},
  note         = {Machine review of arXiv:1909.05353}
}
read the original abstract

In this two-part article series we provide a generalized description of the scattering geometry of Bragg coherent diffraction imaging (BCDI) experiments, the shear distortion effects inherent to the resulting three-dimensional (3D) image from current phase retrieval methods and strategies to mitigate this distortion. In this Part I, we derive in general terms the real-space coordinate transformation to correct this shear, which originates in the more fundamental relationship between the representations of mutually conjugate 3D spaces. Such a transformation, applied as a final post-processing step following phase retrieval, is crucial for arriving at an un-distorted and physically meaningful image of the 3D scatterer. As the relevance of BCDI grows in the field of materials characterization, we take this opportunity to generalize the available sparse literature that addresses the geometric theory of BCDI and the subsequent analysis methods. This aspect, specific to coherent Bragg diffraction and absent in two-dimensional transmission CDI experiments, gains particular importance concerning spatially-resolved characterization of 3D crystalline materials in a realiable, non-destructive manner. These articles describe this theory, from the diffraction in Bragg geometry, to the corrections needed to obtain a properly rendered digital image of the 3D scatterer. Part I provides the experimental BCDI communitcy with the theoretical underpinnings of the 3D real-space distortions in the phase-retrieved object, along with the necessary post-retrieval correction method. Part II builds upon the geometric theory developed in Part I with the formalism to correct the shear distortions directly on an orthogonal grid within the phase retrieval algorithm itself, allowing more physically realistic constraints to be applied.

Figures

Figures reproduced from arXiv: 1909.05353 by the authors.

Figure 1
Figure 1. (a) Basic anatomy of a BCDI measurement of an isolated crystalline nano￾particle. Rotating the scatterer in small increments (for instance about the θ￾direction) causes the reciprocal lattice point q0 of the scatterer’s crystal structure to sweep an incremental angle in Fourier space. The ‘rocking’ of the crystal’s posi￾tion about the Bragg condition effectively causes the measurement plane of the area detector to q… view at source ↗
Figure 2
Figure 2. Basic geometry of a BCDI measurement. Also shown are the laboratory frame [PITH_FULL_IMAGE:figures/full_fig_p030_2.png] view at source ↗
Figure 3
Figure 3. Isosurface plots of the reconstructed object ( [PITH_FULL_IMAGE:figures/full_fig_p030_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: Comparison of the shear-corrected BCDI rendering of the SiC nano-particle [PITH_FULL_IMAGE:figures/full_fig_p031_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

8 extracted references · 8 canonical work pages

  1. [1]

    & Chamard, V

    Berenguer, F., Godard, P., Allain, M., Belloir, J.-M., Talneau, A., Ravy, S. & Chamard, V. (2013). Phys. Rev. B , 88, 144101. URL: https://link.aps.org/doi/10.1103/PhysRevB.88.144101 Calvo-Almaz´ an, I., Allain, M., Maddali, S., Chamard, V. & Hruszkewycz, S. O. (2019). Sci- entific Reports, 9(1),

    work page doi:10.1103/physrevb.88.144101 2013

  2. [3]

    Top row: Direct isosurface plot of the scatterer from the phase retrieval solution array, without the required shear correction

    Isosurface plots of the reconstructed object ( XY , YZ and XZ views), with the color scale depicting complex phase in radians. Top row: Direct isosurface plot of the scatterer from the phase retrieval solution array, without the required shear correction. Axis units are in pixels. Bottom row: Isosurface plots after the shear correction has been applied (r...

    work page 2014

  3. [6]

    (a) Basic anatomy of a BCDI measurement of an isolated crystalline nano- particle. Rotating the scatterer in small increments (for instance about the θ- direction) causes the reciprocal lattice point q0 of the scatterer’s crystal structure to sweep an incremental angle in Fourier space. The ‘rocking’ of the crystal’s posi- tion about the Bragg condition e...

    work page 2014

  4. [100]

    Pixel array dimensions IUCr macros version 2.1.6: 2014/01/16

    work page 2014

  5. [535]

    C., Leake, S

    URL: http://science.sciencemag.org/content/348/6234/530 Newton, M. C., Leake, S. J., Harder, R. & Robinson, I. K. (2009). Nature Materials, 9, 120 EP –. URL: https://doi.org/10.1038/nmat2607 Pateras, A. (2015). Three dimensional X-ray Bragg ptychography of an extended semiconduc- tor heterostructure. Ph.D. thesis, Aix Marseille University. Thse de doctora...

    work page doi:10.1038/nmat2607 2009

  6. [2015]

    I., Allain, M., Godard, P., Largeau, L., Patriarche, G., Talneau, A., Pantzas, K., Burghammer, M., Minkevich, A

    URL: http://www.theses.fr/2015AIXM4366 Pateras, A. I., Allain, M., Godard, P., Largeau, L., Patriarche, G., Talneau, A., Pantzas, K., Burghammer, M., Minkevich, A. A. & Chamard, V. (2015). Phys. Rev. B , 92, 205305. URL: https://link.aps.org/doi/10.1103/PhysRevB.92.205305 Pfeifer, M. (2005). Structural Studies of Lead Nanocrystals Using Coherent X-ray Diff...

    work page doi:10.1103/physrevb.92.205305 2015

  7. [3776]

    & Verdier, M

    URL: https://doi.org/10.1038/s41467-018-06166-5 Dupraz, M., Beutier, G., Rodney, D., Mordehai, D. & Verdier, M. (2015). Journal of Applied Crystallography, 48(3), 621–644. URL: https://doi.org/10.1107/S1600576715005324 Fienup, J. R. (1982). Appl. Opt. 21(15), 2758–2769. URL: http://ao.osa.org/abstract.cfm?URI=ao-21-15-2758 Fienup, J. R. (1987). JOSA A, 4(...

    work page doi:10.1038/s41467-018-06166-5 2015

  8. [6386]

    J., Maser, J., Fuoss, P

    URL: https://doi.org/10.1038/s41598-019-42797-4 Cha, W., Ulvestad, A., Allain, M., Chamard, V., Harder, R., Leake, S. J., Maser, J., Fuoss, P. H. & Hruszkewycz, S. O. (2016). Phys. Rev. Lett. 117, 225501. URL: http://link.aps.org/doi/10.1103/PhysRevLett.117.225501 Cherukara, M. J., Pokharel, R., O’Leary, T. S., Baldwin, J. K., Maxey, E., Cha, W., Maser, J...

    work page doi:10.1038/s41598-019-42797-4 2016

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.