---
title: Operando Synchrotron Nano-holo-Tomography
url: https://www.emergentmind.com/topics/operando-synchrotron-x-ray-nano-holo-tomography
type: topic
---

# Operando Synchrotron Nano-holo-Tomography

Searching arXiv for recent papers on synchrotron X-ray holo-tomography, operando holotomography, and reconstruction methods.
Operando synchrotron X-ray nano-holo-tomography is a propagation-based, phase-contrast tomographic imaging methodology that combines coherent hard X-ray illumination, multi-distance holographic projection acquisition, and tomographic reconstruction to recover three-dimensional maps of refractive-index decrement and, in suitable regimes, electron density at nanoscale resolution during or near device-relevant conditions. In the literature considered here, the method is established through non-destructive synchrotron X-ray holographic tomography of buried three-dimensional silicon photonic crystals [1808.01392], extended algorithmically by software and inverse methods for holographic and tomographic X-ray imaging [2506.11567; 2407.20304], and demonstrated in fully operando 4D form for a silicon-graphite battery anode, where repeated multi-distance scans are coupled with Digital Volume Correlation and transport analysis to connect evolving nanostructure to chemomechanical function [2508.06413].

## 1. Definition and methodological scope

Synchrotron X-ray holographic tomography, also referred to as nano-holo-tomography, is described as propagation-based (inline) phase-contrast imaging with multi-distance phase retrieval and tomographic rotation to reconstruct a 3D map of the refractive index decrement $\delta$ and, via $\delta$, the electron density [1808.01392]. The acquisition uses a focused, coherent hard X-ray beam to illuminate the sample with a spherical wave; radiographs are recorded at multiple propagation conditions while rotating the sample, after which phase retrieval is performed projection-wise and the recovered phase projections are reconstructed tomographically [1808.01392].

The relevant imaging physics is expressed through the complex refractive index
$$
n(x,y,z) = 1 - \delta(x,y,z) + i \beta(x,y,z),
$$
with absorption coefficient
$$
\mu = \frac{4\pi \beta}{\lambda},
$$
exit-wave amplitude
$$
A(x,y) = \exp\!\left[-\frac{2\pi}{\lambda}\int \beta(x,y,z)\,dz\right],
$$
and projected phase
$$
\phi(x,y) = -\frac{2\pi}{\lambda}\int \delta(x,y,z)\,dz = -k\int \delta(x,y,z)\,dz,
$$
where $k=2\pi/\lambda$ [1808.01392]. As a general phase-contrast framework, the Transport-of-Intensity Equation is
$$
\frac{\partial I}{\partial z} = -k \nabla \cdot (I \nabla \phi),
$$
although the cited silicon photonic-crystal study employed a Paganin-type single-material phase-retrieval approach extended to multiple propagation distances rather than directly solving the TIE [1808.01392].

In the operando context, the method becomes 4D when successive 3D reconstructions are acquired during a realistic process sequence. In the silicon-graphite anode study, operando synchrotron X-ray nano-holo-tomography is defined as a 4D $(3\mathrm{D}+\mathrm{time})$ methodology using propagation-based phase contrast in the near-field Fresnel regime to recover quantitative phase images of a low-$Z$, multiphase battery electrode during its formation cycle, then reconstruct them tomographically at nanoscale resolution with minutes temporal sampling [2508.06413]. This establishes the distinction between non-destructive in situ tomography, in which no irreversible preparation is required, and true operando imaging, in which the sample is imaged under an active thermal, electrochemical, fluidic, optical, or electrical stimulus [1808.01392; 2508.06413].

## 2. Imaging geometry, beamline configurations, and acquisition protocols

The method has been implemented at high-coherence synchrotron nanoprobe beamlines using focused hard X-ray illumination and geometric magnification. In the ESRF ID16A-NI study on silicon photonic crystals, the source was the synchrotron ID16A-NI nano-imaging beamline, operated at photon energy $E = 17$ keV with wavelength $\lambda \approx 0.729$ Å, using multilayer-coated Kirkpatrick–Baez mirrors to produce a diffraction-limited focused spot of $23 \times 37$ nm$^2$ [1808.01392]. The sample was placed a short distance $z_s$ downstream of the focus and the detector at distance $z_d$ downstream of the sample, with spherical-wave illumination and Fresnel scaling relations
$$
D = \frac{z_s z_d}{z_s + z_d}, \qquad M = \frac{z_s + z_d}{z_s}.
$$
Four different sample-to-focus distances were used, with 1500 radiographs acquired at each distance while rotating the sample from $0^\circ$ to $180^\circ$; exposure time was $0.3$ s per projection, for a total of 6000 projections [1808.01392].

A later nano-holotomography reconstruction study at ESRF ID16A used the Projection X-ray Microscope with KB mirrors at $17.05$ keV and four sample planes with $z_{1,j} = 4.026$, $4.199$, $4.890$, and $6.325$ mm, total focus-to-detector distance $1.208$ m, geometrical magnifications approximately $300$, $287.7$, $247$, and $191$, and $1500$ angles over $180^\circ$ at $0.1$ s exposure per projection [2407.20304]. The detector had $4096 \times 4096$ pixels with $15\,\mu$m pitch, objective magnification $10\times$, and an effective pixel of approximately $10$ nm; data were typically binned $2\times2$ for processing [2407.20304].

The 2025 battery study demonstrates explicitly operando acquisition. It used the ESRF ID16B hard X-ray nanoprobe at $29.4$ keV with photon flux approximately $3.4\times10^{11}$ ph s$^{-1}$ and a focused nanobeam of about $55$ nm (horizontal) by $60$ nm (vertical) in cone-beam geometry [2508.06413]. For each projection angle, images were recorded at four sample-detector distances. The operando geometry placed the sample approximately $52$ mm from the focal spot and approximately $280$ mm from the detector, giving voxel size $50$ nm, field of view approximately $102 \times 102\,\mu$m$^2$, and reconstructed volume approximately $102 \times 102 \times 102\,\mu$m$^3$ [2508.06413]. Each scan used $2205$ projections over $360^\circ$ with exposure time $35$ ms per projection, plus $20$ flats and $21$ darks per distance, for a total operando scan time of approximately $9$ minutes per four-distance dataset [2508.06413].

These configurations illustrate the defining acquisition pattern of nano-holo-tomography: coherent nanofocused or quasi-point-source illumination, geometric magnification, multi-distance Fresnel imaging, and dense angular sampling. A plausible implication is that the principal design variable is not a single propagation distance but a coupled geometry that must simultaneously satisfy phase-contrast transfer, magnification, field-of-view coverage, and mechanical stability constraints.

## 3. Phase retrieval, probe treatment, and tomographic inversion

The standard reconstruction pipeline in the photonic-crystal study consists of scaling radiographs to a common magnification, registering them across distances and angles, performing initial phase retrieval using the Paganin et al. single-material method assuming a homogeneous $\delta/\beta$ ratio, refining the phase map with 15 iterations of non-linear least-squares optimization, and reconstructing $\delta(x,y,z)$ by filtered backprojection in ESRF PyHST2 [1808.01392]. For silicon at $17$ keV, the single-material ratio used was $\delta/\beta = 174$ [1808.01392]. Electron density was then obtained through
$$
\rho_e(x,y,z) = \frac{2\pi}{r_e \lambda^2}\,\delta(x,y,z),
$$
equivalently,
$$
\delta = \frac{r_e \lambda^2}{2\pi}\rho_e
$$
[1808.01392].

The 2025 battery work follows the same broad logic but in a multiphase porous composite. After multi-distance stack alignment and flat-field/dark-field correction, a Paganin-like phase-retrieval filter with $\delta/\beta = 2744$ was applied to the four-distance stack for each angle, followed by FBP in PyHST2, ring artifact correction, and conversion to 16-bit [2508.06413]. The study states explicitly that the constant $\delta/\beta$ assumption is a source of systematic bias in multiphase objects [2508.06413]. This is a central methodological caveat: the phase retrieval is quantitative within the adopted prior, but compositional heterogeneity violates the single-material model.

A major algorithmic development addresses probe-related artifacts in coherent nano-holotomography. The conventional flat-field step divides each sample projection by a reference image taken without the sample, but when the illumination is coherent this approach overlooks the phase component of the probe, and focusing optics can introduce high-frequency probe structure [2407.20304]. The cited work formulates a conic-beam holotomography forward model in which the complex probe and the angular sequence of sample transmissions are jointly retrieved from multi-distance data [2407.20304]. The preferred objective is the amplitude-space least-squares functional
$$
F_1(\{T_{\theta_k}\},q)
=
\sum_{k,j}\left\|\,|L_{k,j}(T_{\theta_k})|-\sqrt{\tilde d_{k,j}}\,\right\|^2
+
\left\|\,|\mathcal{P}_{\zeta_0}\mathcal{S}^{r}_{s} q|-\sqrt{\tilde d^{r}}\,\right\|^2,
$$
with Wirtinger-calculus gradients and nonlinear conjugate gradients using the Dai–Yuan update [2407.20304]. Experimentally, this scheme resolved twice thinner layers in a 3D ALD standard sample measured using nano-holotomography [2407.20304]. In the reported ALD case, layers as thin as approximately $40$ nm were resolved, versus approximately $80$ nm with conventional TXM full-field imaging, and synthetic-data SSIM improved from $0.49$ to $0.81$ in the noise-free case and from $0.44$ to $0.61$ with Poisson noise [2407.20304].

The software dimension is represented by HoToPy, a Python toolbox for holographic and tomographic X-ray imaging [2506.11567]. It supports both a direct contrast regime and a holographic regime, with explicit implementations of CTF, ICT, nonlinear Tikhonov, TikhonovTV, AP, Paganin, Generalized Paganin, Bronnikov-aided correction, and Modified Bronnikov [2506.11567]. HoToPy is implemented in PyTorch for GPU acceleration and automatic differentiation, uses ASTRA for high-performance tomographic primitives, and provides modules for phase retrieval, preprocessing, alignment, ring removal, and trajectory modeling [2506.11567]. The toolbox therefore systematizes the algorithmic ecosystem around nano-holo-tomography rather than defining a new imaging physics.

## 4. Operando implementation: environmental integration, stability, and throughput

The 2018 photonic-crystal demonstration is explicitly non-destructive and in situ in the sense that the samples were mounted “as is” without cutting, thinning, or heavy-element doping, but no operando stimuli were applied during X-ray imaging [1808.01392]. The same source states that to make the technique operando one would integrate sample environments such as heating, biasing stages, microfluidic cells, or optical excitation compatible with near-focus geometry and rotational tomography [1808.01392]. It also notes that current acquisition of four distances with 1500 projections and $0.3$ s exposures implies tens of minutes per tomogram, so time-resolved operando imaging would require reducing projections or distances, faster detectors, continuous rotation, and potentially compressive or iterative reconstruction strategies [1808.01392].

The battery study realizes this transition in practice. It used a miniaturized operando half-cell with approximately $1.5$ mm diameter housing and $0.7$ mm sample, PFA body, stainless-steel connectors, constant compression of approximately $23$ kPa, shortened upper connector to minimize wobble during rotation, and UV glue sealing [2508.06413]. Electrochemistry was run galvanostatically at C/6 between $5$ mV and $1$ V vs Li$^+$/Li, with a potentiostatic hold at $5$ mV and tomography performed at open circuit during scans; OCV drift was approximately $50$ mV over $12$ min [2508.06413]. Imaging during OCV, small rotating mass, rigid compression, short exposures, and dose control were all part of the stability strategy [2508.06413].

Dose and damage considerations become much more stringent in operando settings than in the silicon-photonics case. The battery study reports cumulative dose approximately $0.52$ MGy in electrolyte and approximately $0.39$ MGy in active material over the entire experiment, below a previously established threshold of approximately $1.5$ MGy for observable beam-induced artifacts [2508.06413]. This contrasts with the photonic-crystal paper, which did not report dose and notes only contextual values for other material classes [1808.01392]. A plausible implication is that operando nano-holo-tomography is constrained not only by temporal resolution and mechanical stability but also by dose budgeting across an entire time series.

Algorithmic support for operando deployment is supplied by HoToPy, which provides center-of-rotation estimation, tilt estimation, registration-based reprojection alignment, bandpass filtering, directional Fourier filters, wavelet ring removal, and per-projection geometry modeling in ASTRA [2506.11567]. These capabilities are explicitly described as relevant for operando conditions with mechanical or thermal drifts [2506.11567]. The simultaneous-probe-retrieval framework is similarly positioned as operando-relevant because it reduces reliance on frequent empty-beam measurements, allows chunked processing with periodic probe re-estimation, and is compatible with streaming reconstruction workflows, although time-resolved 4D reconstruction is not claimed in HoToPy itself [2407.20304; 2506.11567].

## 5. Quantitative performance and representative case studies

### Representative implementations

| Study | Imaging context | Key reported performance |
|---|---|---|
| [1808.01392] | Buried 3D silicon photonic band-gap crystals | Real-space 3D density distributions with 20 nanometer resolution |
| [2407.20304] | ALD standard in nano-holotomography with simultaneous probe retrieval | Resolved layers as thin as ~40 nm; SSIM 0.81 vs 0.49 noise-free |
| [2508.06413] | 4D operando silicon-graphite anode | Voxel size 50 nm; effective spatial resolution typically ~150–200 nm; ~9 min per scan |

In the silicon photonic-crystal study, the central scientific result is the structural discrimination of samples that appear similar in scanning electron microscopy but differ radically in photonic function [1808.01392]. The “Good” crystal showed the designed 3D inverse woodpile periodicity, with $Z$-directed pores of depth $D = 6280 \pm 20$ nm and radius $R = 183 \pm 10$ nm, aspect ratio $17.15 \pm 0.04$, and $X$-directed pores of maximum depth $D = 9460 \pm 20$ nm and aspect ratio $25.8 \pm 0.1$; tomography also revealed a slight $2^\circ$ shear between pore arrays, within a typical range $0$–$6^\circ$, and optical reflectivity showed a broad band gap consistent with theory and tomography-derived geometry with $R/a \approx 0.267$ [1808.01392]. The “Bad” crystal lacked a band gap because tomography revealed a buried stiction-induced internal void formed during drying after colloidal quantum dot infiltration. The “Ugly” crystal had very shallow pores of about $70$ nm owing to erroneous etching parameters and therefore lacked the required 3D periodic depth [1808.01392]. The significance is methodological as much as substantive: surface-sensitive SEM could not distinguish these buried failure modes, whereas non-destructive 3D electron-density mapping could [1808.01392].

The simultaneous-probe-retrieval study provides a different performance benchmark. For an ALD standard sample, the method improved reconstruction fidelity by accounting for the complex probe rather than assuming flat-field division was sufficient [2407.20304]. The reported reconstruction of the full ALD volume, including phase retrieval and tomography, took approximately four hours on two nodes with $4\times$A100 GPUs each, while tomography of the full volume completed in approximately $15$ minutes for $500$ conjugate-gradient iterations [2407.20304]. The result that “3D reconstruction resolving twice thinner layers” was achieved is concretized experimentally by clean separation of approximately $40$ nm ZnO/Al$_2$O$_3$ layers that were not reliably separated by conventional holotomography [2407.20304].

The HoToPy paper reports computational benchmarks rather than new physical application claims. For phase retrieval on a single angle with two distances at $2160\times2560$ pixels, timings range from approximately $4.1$ ms for unconstrained CTF on an A6000 Ada to approximately $439.5$ ms for Tikhonov with non-positivity and disk-shaped support; for the full catalytic-particle dataset, timings range from approximately $9$ s to approximately $11.8$ min depending on algorithm and constraints [2506.11567]. Tomographic reprojection alignment with $100$ iterations, $8\times8$ binning, and bandpass filtering took approximately $9.3$ min on a 24-core CPU plus A6000 Ada GPU, and FDK reconstruction of the full volume took approximately $2$ min [2506.11567]. These figures indicate that advanced holographic phase retrieval and alignment are already computationally practical at synchrotron scale, though not yet necessarily real-time.

## 6. Multiphysics analysis and scientific reach in true operando studies

The 2025 silicon-graphite anode study extends nano-holo-tomography from structural diagnosis to multiscale process science [2508.06413]. Quantitative phase maps reconstructed at successive time points were combined with Digital Volume Correlation to extract displacement and strain fields and with diffusion/tortuosity analysis to relate evolving pore topology to ionic transport [2508.06413]. For DVC, local affine transforms were computed over cubic correlation windows at a pyramidal sequence of $20 \rightarrow 8 \rightarrow 4\,\mu$m$^3$ with $50\%$ overlap, and strain invariants were derived through the deformation gradient $\mathbf{F}$ and its polar decomposition:
$$
\varepsilon_V = \det(\mathbf{F}) - 1,
$$
$$
\mathbf{U}^{\mathrm{dev}} = \det(\mathbf{F})^{-1/3}\mathbf{U},
\qquad
\varepsilon_D = \|\mathbf{U}^{\mathrm{dev}} - \mathbf{I}\|.
$$
For reference, the small-strain tensor is
$$
\varepsilon_{ij} = \tfrac{1}{2}\left(\partial u_i/\partial x_j + \partial u_j/\partial x_i\right)
$$
[2508.06413].

The electrode was observed to have initial macro-porosity of approximately $38$ vol%, partitioned into approximately $13\%$ inner-particle and approximately $25\%$ interparticle porosity [2508.06413]. During lithiation, a two-step sequence occurred: above about $50$–$60$ mV vs Li$^+$/Li, graphite and porosity increased together, whereas below about $35$ mV, in the silicon activity regime, macro-porosity dropped abruptly, with relative decreases of $-7.5$ to $-17\%$, inner-particle porosity decreasing by about $-22\%$ and interparticle porosity by about $-15\%$ [2508.06413]. During delithiation, interparticle porosity recovered while inner-particle porosity increased beyond the fresh state by $+18\%$ relative, with local pore size in graphite increasing from approximately $280$ nm to approximately $420$ nm [2508.06413]. Electrode thickness changed from about $60\,\mu$m in the fresh state to about $75\,\mu$m at end of lithiation and about $64\,\mu$m at end of cycle [2508.06413].

At the mechanics level, mean orthogonal strain components during lithiation were reported as $\varepsilon_{zz} \approx 8.8\%$, $\varepsilon_{xx} \approx 1.24\%$, and $\varepsilon_{yy} \approx 2.43\%$, with deviatoric strain showing staircase increases aligned with voltage plateaus near $70$, $35$, and $20$ mV vs Li$^+$/Li and about $74\%$ irreversibility by end of formation [2508.06413]. Hot spots of volumetric and shear strain localized near large CBD-Si clusters and at graphite-graphite interfaces, and local maximum tensile strain approached approximately $7\%$, which the study links to risk of binder failure [2508.06413].

Transport analysis was performed by solving steady-state Laplace’s equation across the segmented pore network,
$$
\nabla^2 C = 0,
$$
with Dirichlet conditions at current collector and separator and no-flux lateral boundaries, relating effective diffusivity to porosity and tortuosity by
$$
D_{\mathrm{eff}} = D\frac{\epsilon}{\tau}.
$$
Geodesic tortuosity was estimated from
$$
\tau_{\mathrm{geo}} = \frac{L_{\mathrm{geodesic}}}{L_{\mathrm{Euclidean}}},
$$
and “fast geometric paths” were defined by $L_{\mathrm{geodesic}}/L_{\mathrm{Euclidean}} \le 1.2$ across the electrode thickness [2508.06413]. In-plane tortuosity remained $\lesssim 2.6$, through-plane tortuosity was about $3.3$–$3.8$ during cycling and stabilized near approximately $3.5$ at the end, while the number of persistent fast paths increased from $9$ initially to $13$ at the end [2508.06413]. These paths correlated with local flux maxima and were identified as “highways” for ionic transport [2508.06413].

A scale-bridging analysis grouped 14 local parameters into morphology, diffusion, geometry/position, and mechanics, then correlated them with particle activity and irreversibility using Pearson correlation and an “impact factor” [2508.06413]. The study reports that up to approximately $90\%$ of the particle-level deviation from ensemble behavior can be rationalized by these local parameters, compared with approximately $45\%$ if only macro-scale descriptors such as porosity, tortuosity, and thickness are used [2508.06413]. Activity drivers were ranked, with distance to fast path, particle/electrolyte surface fraction, and inner porosity positive, while average local geometric tortuosity and depth position were negative; irreversibility drivers included number of neighbors, maximum volumetric strain at lithiation, irreversible deviatoric strain, and proximity to CBD-Si cluster [2508.06413]. This demonstrates that operando nano-holo-tomography can move beyond morphology tracking into multiscale structure-function inference.

## 7. Relation to complementary methods, misconceptions, and current limitations

The method differs fundamentally from absorption CT because weakly absorbing materials such as silicon at hard X-ray energies yield poor attenuation contrast, whereas holographic tomography derives contrast from phase shifts and is therefore sensitive to electron-density variations in weakly absorbing matter [1808.01392]. It also differs from ptychographic X-ray computed tomography, which uses scanning coherent diffraction imaging with overlapping probe positions and iterative far-field phase retrieval; PXCT can achieve very high resolution but often requires smaller samples or destructive thinning and heavier computation [1808.01392]. SEM reveals surfaces only, FIB-SEM serial slicing is destructive and volume-limited, TEM tomography requires sub-micrometre lamellae, and SAXS or coherent scattering provide reciprocal-space information rather than local real-space 3D maps [1808.01392]. The battery study adds that ptychography and scanning nano-XCT can reach higher spatial resolution but are slower and can be less compatible with minutes-scale operando full-electrode imaging, whereas STXM and BCDI offer single-particle chemical or strain sensitivity but limited representativeness [2508.06413].

A common misconception is that nano-holo-tomography is inherently quantitative simply because it reconstructs phase. The cited literature qualifies this strongly. Quantitativeness depends on the validity of assumptions such as homogeneous $\delta/\beta$ in Paganin-type retrieval [1808.01392; 2508.06413], adequacy of probe treatment under coherent illumination [2407.20304], and stability of geometry, magnification, and alignment [2506.11567]. Another misconception is that flat-field correction is sufficient preprocessing in coherent cone-beam holography. The simultaneous-probe-retrieval study shows explicitly that this is not the case because $\mathcal{P}\{q\psi\} \neq \mathcal{P}\{q\}\mathcal{P}\{\psi\}$; ignoring probe phase injects line artifacts that later amplify in CT [2407.20304].

Several limitations remain explicit in the record. The photonic-crystal paper does not report exact $z_s$ and $z_d$ values, detector native pixel size, field of view, dose, or the method used to estimate resolution [1808.01392]. The battery study states that the effective spatial resolution is typically about $150$–$200$ nm despite $50$ nm voxels, sufficient for graphite microporosity and large CBD-Si clusters but not for individual approximately $200$ nm silicon nanoparticles [2508.06413]. It further notes that only about $4.1$ vol% of CBD-Si clusters are visible versus about $18$ vol% expected, implying that approximately $75\%$ of Si is mixed at sub-resolution [2508.06413]. The simultaneous-probe-retrieval framework assumes a stationary probe within an angular chunk, so stronger time-varying probe behavior would require more frequent probe updates [2407.20304]. HoToPy does not claim real-time streaming, dedicated sparse-view recovery, explicit detector PSF modeling, or automated dose optimization [2506.11567].

Taken together, these studies define operando synchrotron X-ray nano-holo-tomography as a method family rather than a single fixed protocol. Its core is propagation-based multi-distance phase-contrast tomography under coherent synchrotron illumination; its present frontier lies in robust inverse methods, stability-aware acquisition, and coupling to downstream analyses such as DVC, segmentation, and transport modeling. The available evidence shows that the method can non-destructively image buried nanostructures in millimetre-scale supports [1808.01392], can be substantially improved by explicit treatment of the complex probe [2407.20304], is now supported by modular GPU-enabled reconstruction software [2506.11567], and has matured into a genuinely operando 4D platform capable of revealing multiscale chemomechanical dynamics in realistic electrochemical materials [2508.06413].

Source: https://www.emergentmind.com/topics/operando-synchrotron-x-ray-nano-holo-tomography