1. Introduction

X-ray microscopy has unique characteristics of short wavelength and high penetration, and is especially well suited to high resolution imaging of several micrometer thick biological specimens [1, 2]. As a result, a number of groups have developed high resolution transmission x-ray microscopes (TXMs) using Fresnel zone plate optics for routine, high resolution imaging of biological specimens in both 2D [3, 4, 5] and more recently in 3D via tomography [6, 7], and commercial lab-based [8] and synchrotron-based [9] systems have begun to appear. These TXM systems complement Fresnel zone plate based scanning transmission x-ray microscopes (STXMs) [10, 11, 12] by delivering higher throughput for a given source, which is especially important when collecting a tilt series of images for tomography.

In any full-field transmission imaging system, one must collect light scattered by the specimen and use it to obtain a projection image (in 2D) or a volume reconstruction (in 3D). In a TXM operating in brightfield incoherent imaging mode, this is done by using a condenser lens to illuminate the specimen, and an objective lens to collect the scattered light and deliver a phased, magnified, real-space image to a detector. Fresnel zone plates are usually used for the high resolution objective lens, with the image typically recorded on an x-ray sensitive charge-coupled device (CCD) detector.

An alternative, lensless scheme was proposed in 1980 by Sayre [13], and demonstrated in 1999 by Miao et al. [14]: illuminate an isolated specimen with a coherent wave, collect the far-field diffraction pattern, and obtain a reconstructed image by computationally phasing and then inverting the diffraction pattern (such as by using iterative finite-support phase retrieval methods pioneered by Feinup [15]). X-ray diffraction microscopy (XDM; also called x-ray coherent diffraction imaging or CDI) has the advantage that no optics-imposed losses on either the efficiency, or finite numerical aperture, are imposed upon the experiment. However, in x-ray diffraction microscopy the scattered signal must be phased to obtain an image; since the scattered signal drops off significantly at higher spatial frequencies, one must be able to computationally phase weak scattering signals.

In this work we address the following question through simulations: which approach yields more information on an example biological specimen for a given radiation dose? This is an important consideration for radiation-sensitive specimens such as frozen hydrated cells, since in x-ray microscopy (as in electron microscopy) radiation damage is the ultimate resolution-limiting factor.

To answer this question, we simulate x-ray imaging of two biological-cell-like computed objects consisting of lipid membrane, low concentration protein solute, and protein bars of various diameters and orientations inside the “cells.” For diffraction microscopy, we simulate the illumination of the cells by a coherent x-ray beam, calculate the resulting diffraction intensity distributions, add shot noise, and reconstruct the images using an iterative phase retrieval algorithm. For conventional microscopy, we calculate the absorption profiles of the cells, convolve them with the point spread function of a representative soft x-ray zone plate, reduce the image intensities according to the efficiency of the zone plate, and add shot noise to the resulting intensity images. We will discuss results of these comparative simulations and the implications for x-ray microscopy of radiation-sensitive specimens.

2. Specimen illumination, dose, and exposure

Consider a specimen that is illuminated by an average of n̄ photons per pixel, with a pixel dimension of Δ (and thus an areal exposure of n̄/Δ²). If an incident beam is attenuated according to the Lambert-Beer law of I=I ₀ exp[-µz], where µ is the linear absorption coefficient, then the fraction of energy absorbed per thickness is given by -dI/dz=µ I ₀ exp[-µz], where z is the distance into which the beam has already penetrated the specimen. We can therefore describe the “skin dose” D at z=0 by

where ρ is the density of the specimen material and E_photon is the photon energy. The skin dose _D calculation above assumes no escaped energy, which is usually the case for soft x-ray microscopy because of the short range of secondary electrons produced following x-ray absorption and the very low fluorescence yield of the low-Z elements that make up the main part of the mass of biological specimens. The relative biological effectiveness or RBE factor is considered to be equal to unity as usually the case with x-ray photons. Based on the above, we assume a direct proportionality between photons per pixel n̄, areal exposure n̄/Δ², and absorbed dose D in what follows.

3. Signal and noise in images

In Sec. 3.1, we consider the expected relationship between photon exposure n̄ and the signal-to-noise ratio SNR for detecting an object of a given contrast given signal-dependent noise. This provides a basis for the expectation that SNR should scale with n̄. In Sec. 3.2, we consider the problem of estimating the signal-to-noise ratio SNR for imaging of an object based on image correlations in the presence of signal-dependent noise. This provides a method which is used in subsequent SNR calculations based on simulated images.

3.1. Signal-to-noise ratio, and detection of objects

We wish to consider the signal-to-noise ratio for the detection of an object with a specified contrast. We follow here the treatments of Glaeser [16] and Sayre et al. [17]. Our goal is to distinguish measurements (such as image pixels) in the case that a unit-incident-flux measurement would lead to a measured intensity of I. To detect an object, assume that the signal is given by the difference between an intensity with a feature present I_f, and the background intensity I_b; that is, the average signal is n̄|I_f-I_b| when each measurement is made with n̄ incident photons. In the Gaussian approximation to the Poisson distribution for event counting, the noise can be expressed as the square root variance of each measurement √n̄I, and if the noise is uncorrelated from measurement to measurement (which is the case for photon noise) the total noise will be given by the root mean square sum. As a result, the expected signal-to-noise ratio (SNR) is

where Θ of

is a contrast parameter [17] which differs slightly from the usual definition of contrast C as

Based on Eq. (2), a log-log plot of SNR versus incident photon number n̄ should follow a straight line with a slope of 1/2 and an ordinate intercept of log₁₀Θ:

3.2. Estimation of the signal-to-noise ratio from image correlations

The easiest way to measure signal-to-noise ratio is to compare measured images of an object against the known object itself. While this approach can be used in simulations, we wish to employ a method that can be used in experiments where the object is not known except through images that have been obtained. An intuitive approach is to calculate the cross correlation coefficient of two independent images of the same object, so as to compare features that are reproduced against those that fluctuate. We summarize here an approach described by Bershad and Rockmore [18] and first used in electron microscopy by Frank and Al-Ali [19]. We extend this method to signal-dependent noise cases, and cases where the signal is not required to have a mean of zero (such as in intensity rather than amplitude measurements).

Let us consider two separate 2D intensity measurements I ₁ and I ₂ of the same signal array S with stochastic noise arrays N ₁ and N ₂, or

where the mean value of S over all 2D pixels is 〈S〉, and N ₁ and N ₂ come from the same source obeying the same Gaussian distribution with zero mean value (〈N ₁〉=〈N ₂〉=0). As a result, 〈I ₁〉=〈I ₂〉=〈S〉.

The imaging of a 2D object depends both on how the signal changes across the object due to illumination and the object’s intrinsic contrast, and on fluctuations due to noise. We characterize both of these by defining the total signal and noise for the entire 2D image with their variances:

where 〈N _1,2〉=0 has been used in the final equality of Eq. (8). Again, the average is done over all pixel indices of one 2D image, which differs from a variance calculation for a particular pixel in a set of separately measured images. The variances of images I ₁ and I ₂ can be calculated as

where we used the fact that the cross terms 〈SN _1,2〉 are in practice negligible compared to 〈S ²〉. This would be expected for low contrast objects, but in practice it also applies to cases with high contrast objects; in simulations, we find that 〈SN _1,2〉 are several orders of magnitude lower than 〈S ²〉 even with full contrast objects. As a result, we are able to drop the terms 〈SN _1,2〉 that would otherwise appear in Eq. (9). Since N ₁ and N ₂ obey the same Gaussian distribution, they have identical variances (〈N ² ₁〉=〈N ² ₂ 〉). So, we can denote σ₁=σ₂=σ. The covariance between I ₁ and I ₂ can be derived following analogous steps:

where r is the correlation coefficient. It can be calculated as

Combining Eqs. (7) and (8) with Eq. (9), we have

From Eqs. (7) and (10), we have

From Eqs. (12) and (13), we see that the signal-to-noise ratio can be calculated from the correlation coefficient r as

which is the square root of the expression α=r/(1-r) used by Frank and Al-Ali [19].

 

Fig. 1. (a) Images with noise simulated according to different mean photon numbers per pixel n̄. Without noise, each pixel p has a value I_p of between 0 and 1. With illumination of n̄ photons per pixel, each pixel has a starting value of n̄I_p photons which was used on a pixel-by-pixel basis to generate an image using signal-dependent noise calculated using the positive-integer-truncated Gaussian approximation. (b) signal-to-noise ratio (SNR) calculated from the images versus incident photon number per pixel n̄. The signal-to-noise ratio was calculated both by comparison of the noisy image with the noise-free original image described in Eqs. (7) and (8), and by using the two-noisy-image correlation method described in Eqs. (11) and (14). Since a fit of the data on a log-log plot of SNR versus n̄ shows a slope of 0.5 as expected from Eq. (5), both image SNR methods give the expected scaling of SNR versus exposure n̄.

Download Full Size | PDF

To confirm that Eq. (14) provides a good measure of the SNR from multiple noisy images of identical objects, we used simulations to explore its scaling with photon number. As an “object” we used the well-known “Lena” image with transparency I normalized between 0 and 1. For each calculational run we multiplied the image by a value of n̄, and then added to each pixel p a pseudo-random noise value based on a positive-integer-truncated Gaussian approximation of the true Poisson distribution for each pixel’s value of n̄I_p. By repeating this process, we obtained two images of an identical object but with different signal-dependent noise added, allowing us to measure the image SNR using the cross-correlation result of Eqs. (11) and (14). Figure 1(a) shows the resulting images with noise corresponding to the indicated number n̄ of incident photons per pixel. Figure 1(b) shows that a log-log plot of the calculated SNR values as a function of n̄ has a slope of 1/2 as expected from Eq. (5), reflecting the √n̄ dependence of SNR on incident exposure expected from Eq. (2). Finally, we also show in Fig. 1(b) the signal-to-noise ratio calculated by comparison of the original, noise-free image with a noise-included version at each photon exposure n̄ using Eqs. (7) and (8). The result shows that these SNR values are identical to those calculated from the cross-correlation method. However, in subsequent calculations we use the two-noisy-image correlation method of Eqs. (11) and (14) because it is applicable to a broader range of cases including experiments where the noise-free object is not known.

4. Simulations with defined “cells”

Having established the image correlation method of Eqs. (11) and (14) for estimating SNR from image pairs, we now turn to simulations of the two experimental approaches to be compared (Fig. 4): lens-based imaging (incoherent brightfield imaging in a TXM), and lensless imaging (x-ray diffraction microscopy or XDM). We discuss the “cells” to be “imaged” in Sec. 4.1, the TXM imaging process in Sec. 4.2, and the XDM imaging process in Sec. 4.3.

4.1. Defined “cells” A and B

We generated two different types of biological-cell-like defined objects for our simulations, which we will call cell A and cell B (see Fig. 2). In both cases we assumed an x-ray energy of 520 eV (within the “water window” [17, 20] between the carbon and oxygen K absorption edges), and values of the refractive indices calculated according to the tabulation of Henke et al. [21] using an assumed stoichiometric composition of H_48.6C_32.9N_8.9O_8.9S_0.3 and density of ρ=1.35 g/cm³ for protein, and H_62.5C_31.5O_6.3 with ρ=1.0 for lipid [22]. The two objects were defined as follows:

 

Fig. 2. Defined objects used for image simulations. Shown here are the magnitudes of the simulated exit waves resulting from plane wave illumination of the objects. Cell A has random protein thicknesses within an irregular boundary, while Cell B has a lipid membrane and several protein bars and ellipses inside.

Download Full Size | PDF

Cell A: The first defined object is modeled as a 2D continuous region with irregular boundary and random protein thicknesses over a range of 0 to 500 nm inside (Fig. 2(a)). The entire image is a 256×256 pixel array with a pixel size of 15 nm.

Cell B: The second defined object is embedded in the center of a 400×400×400 pixel ice cube, also with a pixel size of 15 nm. The cell-like object has a diameter of 200 pixels, or 3 microns, with a 3 pixel or 45 nm thick lipid membrane as boundary. Inside the object, several 225 nm diameter protein rods were placed (1.8 µm long for the vertical bars, and 1.35 µm long for the horizontal bars), along with some protein ellipsoids (the larger ones are 0.9×0.45 µm, and the smaller ones are 0.68×0.45 µm in size). A “bud” with a diameter of 1.2 µm was added to the top right shoulder of this fake cell, both in order to break the rotational symmetry for eliminating potential problems with enantiomorphs in XDM, and to approximate the appearance of a budding yeast cell. This object was then assumed to be illuminated by a plane wave, and a multislice propagation process [23, 24] was used to generate an exit wave (400×400 pixels across) leaving the cube.

In the case of cell A, we have a defined object with broad spatial frequency content due to the random thickness variations inside. In the case of cell B, we have an object with readily recognized structures. The two defined objects are shown in real space in Fig. 2, and in Fourier space in Fig. 3. To compare zone plate imaging with diffraction microscopy, we simulated both imaging techniques on these two different defined objects.

 

Fig. 3. Diffraction patterns of the exit waves from the two defined objects, cell A and cell B, and their azimuthally averaged, unit-integral-normalized power spectral densities. In both cases the exit wave amplitudes were Fourier transformed and then squared to yield the diffraction intensity. Both defined objects have the signal decline with spatial frequency in a power law relationship with a slope of about 3.3 over most frequencies. For cell A, this trend then levels off at a spatial frequency of about 10 µm^-1 where the pixel-by-pixel uncorrelated protein thickness dominate the diffraction pattern; this yields a flat power spectrum corresponding to a delta (δ) function in real space. For cell B, the overall round shape of the object gives rise to Airy rings in the diffraction pattern which show up at spatial frequencies above 10 µm^-1.

Download Full Size | PDF

4.2. Zone plate imaging process

X-ray microscopes using synchrotron radiation and Fresnel zone plates for full-field imaging have been in existence for some time [3], and are finding considerable success for 25–40 nm resolution imaging applications, tomographic imaging of frozen hydrated cells [6, 7, 25], with commercial laboratory source versions now becoming available [8]. We model here a representative microscope with a Fresnel zone plate with 30 nm outermost zone width and 10% diffraction efficiency [26], used in incoherent brightfield mode with a 100% efficient detector. In this case the recorded image intensity is a convolution of the intensity transmitted through the object with the intensity point spread function of the imaging system, or

Using the convolution theorem of Fourier transforms ℱ{}, the image can be represented as

where we used the fact that the modulation transfer function or MTF is the Fourier transform of the intensity point spread function. The MTF can be calculated[27] as

where f represents spatial frequency, and f ₀ is the spatial frequency cutoff. For a zone plate with outermost zone width of dr_N (set to be 30 nm in simulation), the cutoff of the incoherent MTF is at a spatial frequency of f ₀=1/dr_N. The integral of Eq. (17) has a numerical value of 0.2 compared to the integral of an MTF of 1 for all frequencies up to the cutoff at f ₀. The resulting image intensity was then multiplied by 0.1 to account for a typical zone plate focusing efficiency of 10% (while theoretical efficiencies can approach 20%, 10% is representative of the best experimental measurements). Finally, the resulting image was multiplied by an exposure of n̄ photons per pixel, and simulated photon noise was added using the positive-integer-truncated Gaussian approach to yield a final image. An image pixel size of 15 nm was used in these simulations.

 

Fig. 4. Schematic of the two x-ray imaging systems considered in simulations. In the TXM or transmission x-ray microscope at left, incoherent brightfield imaging is assumed where the numerical aperture of the condenser is 1.5 times the numerical aperture of the objective lens; a magnified image is recorded on a detector such as a CCD. The TXM objective is a zone plate with 30 nm outermost zone width and 10% diffraction efficiency. In XDM or x-ray diffraction microscopy at right, the specimen is assumed to be illuminated by a fully coherent beam and the far-field x-ray diffraction pattern is recorded on a detector such as a CCD. A reconstructed image is obtained by computational phasing of the coherent diffraction pattern.

Download Full Size | PDF

4.3. Reconstruction from diffraction pattern

In x-ray diffraction microscopy, the image is reconstructed from the far-field diffraction intensity of a coherently illuminated object which is assumed to be within a finite support (area of non-zero optical interaction). To simulate the imaging process, the complex exit waves leaving the defined objects were multiplied by √n̄ to account for an exposure of n̄ photons per pixel, after which they were Fourier transformed to yield the far-field diffraction amplitudes. These diffraction amplitudes were then squared to yield the diffraction intensities, and simulated photon noise was added to produce the simulated data recordings. The diffraction patterns were assumed to be recorded in the far field, with proper Shannon sampling between the real-space object array and the Fourier plane detector array (that is, 256² detector pixels for Cell A, and 400² detector pixels for Cell B). As with the TXM case, a real space image pixel size of 15 nm was used in these simulations. To reconstruct the image, we took the square root of the noisy diffraction intensity to revert to Fourier amplitude, and then reconstructed the image using the difference map algorithm developed by Elser [28] with β=1.15.

5. Simulation results

5.1. Image simulation and SNR calculation

We carried out simulation runs for each of the two defined objects (cell A and cell B), using both TXM and XDM, and using incident photon per pixel values of n̄={1,2,5}×10^{1,2,3,4,5} and 1×10⁶ (that is, exposures of 1×10¹, 2×10¹, …, 1×10⁶). Figure 5 shows some examples of simulated zone plate and reconstructed diffraction images. For each incident photon number n̄, we calculated 10 TXM images and 10 XDM images, each with separately simulated signal-dependent photon noise in the intensity recording. From each set of 10 images, one can form 45 different two-image pairs for calculating the SNR using the method of Eqs. (11) and (14) over the region where the object is located (that is, inside the object’s support for both the XDM and TXM cases). The final SNR value was calculated from the average of these 45 measurements.

 

Fig. 5. Images resulting from the simulations. Shown here are both TXM (top row) and XDM (bottom row) images calculated for the two defined objects, cells A and B, with exposures of n̄=1×10³ and n̄=1×10⁵ photons per pixel. The lack of “salt and pepper” shot noise in background region outside of the cell in the XDM reconstructions is a result of the imposition of a finite support constraint in the reconstruction process.

Download Full Size | PDF

Figure 6 shows plots of log₁₀(SNR) versus log₁₀(n̄) for the set of image simulations of both defined objects. As can be seen, in each case the slope of the fitted line is near 1/2, as expected from Eq. (5); in other words, the signal-to-noise ratio is proportional to √n̄ as expected from Eq. (2). An important difference is that x-ray diffraction microscopy gives a higher SNR for a given photon exposure n̄ in these simulations. For zone plates of 10% efficiency and an incoherent brightfield MTF integral about 0.2 as described in Sec. 4.2, one might expect a factor of 1/(0.1×0.2)=50 in signal loss, or a factor of √50⋍7 in the signal-to-noise ratio decrease for TXM. The fitted lines of Fig. 6 are indeed higher for x-ray diffraction microscopy versus transmission x-ray microscopy, with an improvement of XDM/TXM of 7.2±1.6 for cell A, and 6.3±2.1 for cell B. This means that the iterative reconstruction algorithm used in x-ray diffraction microscopy is able to phase even weak diffracted signals with significant signal-dependent photon noise present. We also note an early observation by Fienup [15] who found that finite support phase retrieval algorithms appeared to be quite robust in the presence of signal-independent noise added to the Fourier magnitudes.

As noted in Fig. 6, the log₁₀(SNR) versus log₁₀(n̄) dependence of cell B shows some oscillations about the linear trend. This oscillation corresponds to successive Airy rings in the diffraction pattern of the overall spherical shape of cell B, as shown in Fig. 3. In the dark bands between Airy rings corresponding to the far-field diffraction pattern of a disk, the signal is lower than the ring-free trend would suggest that the addition of more photons is less helpful until the next Airy ring begins to be filled in. A different effect can be seen in the log₁₀(SNR) versus log₁₀(n̄) curve for cell A: as shown in Fig. 3, the normalized power spectrum for cell A “levels off” at a constant value of about 10^-2 at spatial frequencies f above about 10 µm^-1. This suggests that when the incident photon number n̄ is larger than 100, all the diffraction intensity pixel values approach 1 so that all spatial frequencies above 10 µm^-1 have measurable, single-photon or larger intensities, and start to contribute to increase reconstruction quality. This may explain the slight flattening of the log₁₀(SNR) versus log₁₀(n̄) curve at the lowest photon exposures for cell A in x-ray diffraction microscopy (XDM).

 

Fig. 6. Plots of the signal-to-noise ratio SNR as a function of incident photons per pixel n̄ for our simulation set. The results at left are for the defined object “cell A”, while those at right are for “cell B” (see Fig. 2). The results for both x-ray diffraction microscopy (XDM) and transmission x-ray microscopy (TXM) are shown. In each case a slope of about 1/2 is observed in the log₁₀(SNR) versus log₁₀(n̄) plot (Eq. (5)), and for both objects the SNR for x-ray diffraction microscopy is about 7 times higher than it is for transmission x-ray microscopy. For x-ray diffraction microscopy of cell B, the SNR curve oscillates around the straight line fit. This oscillation corresponds to providing enough signal to phase data in successive Airy rings in the diffraction pattern of the overall spherical shape of cell B, as shown in Fig. 3.

Download Full Size | PDF

5.2. Resolution estimation

As shown in Fig. 3, the azimuthally-averaged diffracted signal P(f) in our defined objects declines in a power law relationship of log10 P(f)=mlog₁₀ f+b or

with m=-3.3 (azimuthal averaging smooths out the sharp variations due to individual speckles to reveal an overall scattering trend). Theoretical estimates for signal decline in x-ray diffraction microscopy have ranged from m=-4 [29], to m=-3, -4, and -6 for various specimen models [30]. In small angle x-ray scattering, Porod’s law suggests that the diffraction signal should decline with spatial frequency as f^m with m=-4, while measured power spectra tend to have positive deviations from this dependence (i.e., values of m of -3 to -4).

As the diffracted signal for a particular length scale declines, the available number of photons for reconstructing structure at that length scale decreases; the signal-to-noise ratio then worsens as a result (as illustrated in Fig. 1). Since images become unrecognizable when a minimum signal-to-noise ratio is not satisfied [31], this suggests a direct relationship between achievable resolution Δ_min and some minimum signal level S _min. Since the signal scales with incident photons per pixel n̄ yet declines with spatial frequency f in the Fourier domain, we will assume that the minimum reconstructable signal level is reached at a cutoff spatial frequency f_c, or

Associating the minimum object size Δ_min with the half-period of a grating at the maximum resolvable spatial frequency f_c, we obtain

Inserting the diffracted signal trend of Eq. (18) evaluated at frequency f_c into Eq. (19) and then employing Eq. (20) leads to

where we have used m_r to denote the power law dependence in reconstructed images even though we expect it to be the same as the specimen diffracted power scaling m of Eq. (18). The result of Eq. (21) can then be rearranged to give

where we have chosen to place Δ_min on the abscissa and n̄ on the ordinate for comparison with the dose-versus-resolution plots estimated by Howells et al. [29] and Shen et al. [30].

 

Fig. 7. Investigation of the dose versus resolution trend for x-ray diffraction microscopy of defined object B. (a) The ratio 〈I _recon(f)〉/〈I _data(f)〉 for the example case of incident photon number n̄=1000. This ratio measures how consistent the reconstructed intensity I_recon averaged over many iterates compared to the recorded intensity I_data. We chose a cutoff value of 0.7 as providing an estimate of the resolution of the reconstructed image. (b) The resulting trend of incident photon number n̄ versus resolution Δ_min, along with a straight line fit to determine the dose versus resolution scaling parameter m_r. The value of m_r=-3.2±0.2 is consistent with the scaling of m=-3.30±0.03 shown in Fig. 3 for the diffracted signal from this object.

Download Full Size | PDF

The above arguments suggest that the required dose (proportional to exposure n̄) for achieving a desired resolution Δ_min should have the same log-log slope m_r as the scaling m in object’s diffracted signal (because of the use of Eq. (19) in deriving Eq. (22)). We therefore explored the scaling parameter m_r in our simulations of imaging defined cell B using x-ray diffraction microscopy. The spatial frequency dependent ratio 〈I _recon(f)〉/〈I _data(f)〉 provides a good measure of the resolution in diffraction microscopy reconstructions [24, 32] (similar to a phase retrieval transfer function [33]), and in our simulations we used a value of 0.7 for this ratio as a way to identify the cutoff spatial frequency f_c and thus the resolution Δ_min=1/(2f_c). The resulting data of signal level n̄ required to achieve a given resolution Δ_min are shown in Fig. 7; a linear fit of the data gives a value of m_r=-3.2±0.2. Since defined cell B has a signal that scales like P(f) ∝ f^m with m=-3.30±0.03 as shown in Fig. 3, these simulation results are consistent with the expected result: the scattering strength of the specimen determines the dose versus resolution trend for imaging the specimen. The same conclusion is also obtained from cell A.

While Fig. 7 shows a good overall agreement between the dose-versus-resolution slope m_r and the scattering-versus-spatial-frequency slope m, there are “local” departures from this “global” trend. That is, the power law fit of Fig. 7 is done over a very large, 10⁵:1 dynamic range in incident photons per pixel n̄; this is the “global” trend. If one were instead to carry out the fit over a smaller exposure dynamic range, a different, “local” power law dependance might be observed. Given that the Airy ring characteristics of our defined object B produce departures from a simple power law relationship as seen in Fig. 3 (due to concentration of extra diffraction signal into certain spatial frequency ranges), we expect that both computer simulations and experimental observations can have somewhat different dose-versus-resolution slopes m_r over particular exposure ranges.

6. Conclusions

The goal of this work was to better understand the tradeoffs between two different x-ray microscopy methods for high resolution imaging of dose-sensitive specimens. Although in these simulations we had perfect knowledge of the objects that were “imaged,” we chose to adopt a cross-correlation method for quantitative comparisons that is better suited to experimental work with unknown objects. Following the work of Bershad and Rockmore [18], we modified the signal-to-noise measure of Frank and Al-Ali [19] with the result that our measure gives the expected scaling with incident photon number n̄. We then used this method to compare incoherent brightfield imaging in a transmission x-ray microscope with a zone plate objective lens of specified efficiency versus x-ray diffraction microscopy using an iterative phase retrieval algorithm for image reconstruction. In these simulations, x-ray diffraction microscopy gave a higher signal-to-noise ratio for equivalent dose, with a gain consistent with the losses imposed by the modulation transfer function and overall efficiency of the zone plate objective.

Of course these simulations are idealized; real imaging experiments involve a number of factors not accounted for here. One of them is that transmission x-ray microscopy has the significant advantage of providing images immediately, unlike x-ray diffraction microscopy where one must use sophisticated image reconstruction algorithms, and specimens that satisfy a finite-support constraint (though we note that ptychography can remove that latter limitation [34]). A limitation of transmission x-ray microscopy as well as its scanning “cousin” is that the resolution-determining outermost zones on zone plates usually have a lower diffraction efficiency than the coarser, inner zones. In x-ray diffraction microscopy, complications include undesired scattering from nearby high-contrast objects such as specimen support grid bars, and the detrimental effects of partial coherence which can contribute noise to reconstructed images [35]. Nevertheless, the simulations shown here illustrate how x-ray diffraction microscopy has the potential to deliver images with higher resolution when studying dose-sensitive specimens.