Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Quantitative characterization of x-ray differential interference contrast microscopy using modulation transfer function

Open Access Open Access

Abstract

Performance of two types of differential interference contrast objectives, i.e., the XOR pattern and the zone-plate doublet, is quantitatively characterized and compared using modulation transfer function. Effects of partial coherence, finite absorption and phase in a complex object, as well as bias retardation are also examined.

©2011 Optical Society of America

1. Introduction

High resolution soft x-ray microscopy is a useful analytic tool in nanosciences [13]. Its nanoscale resolving power, combined with available x-ray phase imaging techniques, has been successfully used to image otherwise transparent samples with nanometer resolution [49]. Comparing with other phase contrast methods available at the x-ray wavelengths such as ptychography [10, 11], interferometry [1214], and Zernike phase plates [79], x-ray differential interference contrast (DIC) microscopy is advantageous in that DIC objectives can be readily implemented in existing x-ray microscopes, thus forgoing dedicated experimental apparatus [46]. In addition, x-ray DIC microscopy has been demonstrated with both synchrotron and compact plasma x-ray sources, making it accessible to a wider range of users [15, 16]. Currently, two types of DIC objectives, i.e., the exclusive OR (XOR) patterns [4, 16, 17] and the zone-plate doublets (ZPD) [5, 18, 19], have been demonstrated. In both objective designs, arbitrary bias retardation and lateral shear can be incorporated to improve image contrast. The XOR pattern’s fabrication method is similar to that of a Fresnel zone-plate (FZP) and as a result the same finest outermost zone width, thus spatial resolution, can be achieved. ZPD on the other hand can be used for non-iterative quantitative x-ray phase imaging [19], although its nanofabrication process requires a more stringent placement accuracy. Figure 1 schematically illustrates the difference in imaging performance between conventional FZP, XOR, and ZPD objectives for a complex sinusoidal grating object. A conventional FZP reproduces the object’s intensity, defined as the absolute square of the object transmittance function; and the effect of object’s phase is not reflected in its image. In contrast, XOR and ZPD use the object’s phase, in addition to its absorption, for their image formation and therefore higher contrast can be obtained. Clearly, XOR and ZPD based imaging systems are more effective in generating image contrast for complex objects.

 figure: Fig. 1

Fig. 1 Performance comparison between FZP (conventional), XOR, and ZPD based x-ray microscopes. The complex object (left column) is illuminated coherently and imaged by FZP, XOR, and ZPD objectives, respectively (center column). Note that the object’s intensity is defined as the absolute square of the object transmittance function. The resultant image intensity (right column) shows that XOR and ZPD objectives generate higher image contrast for the complex object. It is clearly seen that image visibility produced by the two DIC objectives exceeds the object’s visibility while the image visibility produced by FZP is at best equal to that of the object. The value of MTF for DIC objectives can therefore exceed 1 (see text for details). The same bias retardation value 2Δθ = π/2 is used for the XOR and ZPD objectives in this figure.

Download Full Size | PDF

Despite the usefulness of x-ray DIC microscopy, little effort has been devoted to quantify its imaging performance. Modulation transfer function (MTF) is the standard evaluation criterion for imaging system responses, and it has been widely used to quantify absorption-based microscopes [20, 21]. An absorption grating is typically used as the object for MTF analyses in the case of absorption contrast microscopes. Evaluation of phase contrast microscopes requires either a pure phase object with no absorption [22, 23] or a complex object which introduces both phase shift and absorption to the transmitted wave-front. Since most biological and magnetic samples have complex refractive indices at the soft x-ray wavelengths, analyses using complex grating objects render more realistic performance characterization of x-ray phase contrast microscopes. Theoretical analyses using complex objects exist only for conventional absorption contrast and Zernike phase contrast microscopes [24, 25], but not for DIC microscopes. Following the classical approach [20, 26, 27], here we systematically quantify and compare the performance of XOR and ZPD based x-ray DIC microscopes using MTF. The DIC microscopes studied here are configured to be space invariant in order to allow for the exact transfer function analysis without approximation [28, 29]. Performance of the DIC microscopes are evaluated with complex grating objects using varying levels of partial coherence and bias retardation.

2. Modulation Transfer Function for X-ray Microscopy

2.1. Mathematical Description of the Image Formation Process on an X-ray Microscope

A typical partially coherent cascade imaging system is illustrated in Fig. 2. Coordinate systems for the source, object, objective lens, and the image planes are denoted by (α, β), (ξ˜, η˜), (x,y), and (u, v) and the distances between the source, condenser, object, objective lens, and image are z 1, z 2, zo, and zi, respectively. When the relative positions of the condenser, object, and objective lens are properly set, this partially coherent imaging system becomes exactly space invariant and the resultant image intensity distribution can be analytically expressed in terms of transfer functions without low coherence approximation and additional phase correcting lenses [28, 29]. This illumination apparatus is especially useful for x-ray DIC microscopy where a high degree of illumination coherence is generally required. Image intensity distribution in this case is given by [27, 28],

Ii(u,v)=+dpdq𝕁o(p,q)×+dw1dw2𝕋o(w1,w2)𝕂(w1p,w2q)×+dνUdνVexp{j2π(uνU+vνV)}×𝕋o*(w1νU,w2νV)𝕂*(w1pνU,w2qνV),
where 𝕁o, 𝕂, and 𝕋o are the Fourier spectra of the illumination mutual intensity J o, the space invariant amplitude spread function K, and the object transmittance function t o, respectively. Note that following prior classical MTF analyses [20, 27], a circular finite extent incoherent source Is(α, β) is used here to illuminate the microscope, i.e.,
Is(α,β)={Io,α2+β2<Rs0,otherwise,
where Io is a constant and Rs is the source radius. Fourier spectrum 𝕁o of this source is thus,
𝕁o(p,q)=κM2Is(λfcMp,λfcMq),
where λ is the wavelength, M = zi/zo is the magnification, fc is the focal length of the condenser, and κλ 2/π. Similarly, we have
𝕂(p,q)=MP(λzoMp,λzoMq),
where P is the pupil function of the objective.

 figure: Fig. 2

Fig. 2 Schematic diagram of a partially coherent imaging system.

Download Full Size | PDF

2.2. Pupil Functions of the Objective Lenses

The objective lens plays a critical role in the image forming process of the microscope. Three types of objective lenses, i.e., Fresnel zone-plates (FZP) [21, 30], exclusive OR (XOR) patterns [4, 31], and zone-plate doublets (ZPD) [19], are examined here and their pupil functions P(x,y) are summarized in Table 1. Note that 2Δθ is the bias retardation, Δs is the lateral shear, and the parameter γ is defined as π/[Δr(2Ro – Δr)] where Δr and Ro are the outermost zone width and radius, respectively. These pupil functions are used in Eq. (4) for simulations performed in the following sections. Figure 3 shows examples of XOR and ZPD objectives with bias retardation values of 2Δθ = 0, π/2, and π, and Fig. 4 shows their corresponding pupil functions. For comparison, a conventional FZP objective and its pupil function are also shown in the figures. Note that the pupil function of a visible-light DIC microscope is identical to that of ZPD [19, 32], and the MTF analyses presented here for ZPD are also valid for visible-light DIC microscopy.

 figure: Fig. 3

Fig. 3 XOR and ZPD based DIC objectives with various bias retardation values. (a) XOR with Δθ = 0, (b) XOR with Δθ = π/4, (c) XOR with Δθ = π/2, (d) ZPD with Δθ = 0, (e) ZPD with Δθ = π/4, and (f) ZPD with Δθ = π/2.

Download Full Size | PDF

 figure: Fig. 4

Fig. 4 Pupil functions corresponding to the objectives in Fig. 3. (a) XOR with Δθ = 0, (b) XOR with Δθ = π/4, (c) XOR with Δθ = π/2, (d) ZPD with Δθ = 0, (e) ZPD with Δθ = π/4, and (f) ZPD with Δθ = π/2.

Download Full Size | PDF

Tables Icon

Table 1. Pupil Functions P(x,y) of FZP, XOR, and ZPD Objectives

2.3. Periodic Objects: Gratings

Gratings that introduce both absorption and phase variations are needed in order to evaluate the performance of DIC microscopes. Here three types of gratings, i.e., absorption sinusoidal, phase sinusoidal, and complex sinusoidal (Fig. 5), are used and their transmittance functions t o as well as Fourier coefficients an are summarized in Table 2 (Detailed derivations in Appendix A). Note that the refractive index ñ at the soft x-ray spectral region is in general complex and typically denoted by ñ = 1 – δ + where (1 – δ) determines the phase shift introduced and β the absorption. Denoting the maximum thickness variation of the complex sinusoidal grating as to and the attenuation length as ta = 1/(2) where k = 2π/λ is the vacuum wavenumber, the parameter η is defined as the ratio of to to ta. Characteristics of the complex object can then be modeled by varying the values of δ/β and η.

 figure: Fig. 5

Fig. 5 Absorption (top) and phase (bottom) variations of (a) absorption sinusoidal, (b) phase sinusoidal, and (d) complex sinusoidal grating objects. Period of the grating is denoted by T.

Download Full Size | PDF

Tables Icon

Table 2. Transmittance Functions and Fourier Coefficients of the Gratings Used in the Text*

2.4. Implementing the Numerical Evaluation of the Image Integral

Given a grating with period T, its object transmittance function t o(ξ′, η′) can be expressed in a Fourier series as,

to(ξ,η)=n=+anexp(j2πnνoξ),
where νo ≜ 1/T is the fundamental spatial frequency of the grating and an is the nth order Fourier coefficient (Table 2). Note that Eq. (5) is expressed in the magnified coordinate (−Mξ˜, −Mη˜) = (ξ′, η′), representing the geometric image of the object. Fourier transform of the object transmittance function t o(ξ′, η′) is given by,
𝕋o(z1,z2)=n=+anδ(z1nνo,z2),
where δ is the Dirac delta function.

Image intensity distribution of this periodic grating formed by a microscope can therefore be obtained by substituting 𝕋o into Eq. (1), i.e.,

Ii(u,v)=m=+n=+anam*exp{j2π[(nm)νou]}×+dpdq𝕁o(p,q)𝕂(p+nνo,q)𝕂*(p+mνo,q).

The integral in Eq. (7) depends solely on the properties of the imaging system, i.e., the illumination 𝕁o and the amplitude spread function 𝕂 of the microscope, while am and an are the Fourier coefficients of the object being imaged. Since 𝕁o and 𝕂 are given by Eq. (3) and Eq. (4) as scaled versions of Is and P, respectively, the integral can be numerically evaluated by overlaying scaled Is, P, and P*, i.e.,

Ii(u,v)=m=+n=+anam*exp{j2π[(nm)νou]}×κM2+Is(λfcMp,λfcMq)P(λzo(Mp+Mnνo),λzoMq)×P*(λzo(Mp+Mmνo),λzoMq)d(Mp)d(Mq),
and this formulation is used for numerical simulations presented in Sec. 3.

To simplify the notations, we define

νcRoλzo
νsRsλfc,
where νc is the cutoff frequency and νs is the extent of source angular distribution [33]. The effect of partial coherence can be characterized by the σ value, defined as [27, 34],
σνsνc.
As σ increases, the degree of illumination coherence on the object decreases. Coherent and incoherent imaging processes are therefore denoted by σ = 0 and σ = ∞, respectively. The normalized object spatial frequency α is defined by normalizing the object’s fundamental spatial frequency νo with respect to the cutoff frequency νc, i.e.,
αMνoνc,
where M is the magnification. Note that νo is the fundamental spatial frequency of the grating as seen in the image plane while the actual grating spatial frequency is o.

2.5. Modulation Transfer Function and Visibility

Modulation transfer function (MTF) of an optical imaging system is defined by,

MTF(α)Vimage(α)Vobject(α),
where Vimage(α) and Vobject (α) are the visibilities of the image and object, respectively [27]. Note that the object’s intensity is defined as the absolute square of the object transmittance function, i.e., |t o(ξ′, η′)|2, and the visibility for a phase object is therefore zero by definition (see Appendix B). As a result, MTF in Eq. (13) cannot be used in this case, and instead image visibility Vimage is used directly to evaluate the performance of the imaging systems when imaging pure phase objects.

3. Quantitative Characterization using Modulation Transfer Function

Numerical evaluations of the performance of XOR and ZPD based x-ray DIC microscopes are presented in this section. Bias retardation values 2Δθ = 0, π/6, π/2, 2π/3, and π are used for both DIC objectives. Lateral shear Δs is fixed at Δs = 2Δr, where Δr is the outermost zone width, in all cases. Effect of partial coherence is examined by varying the value of σ. Performance of the FZP objective is also shown for comparison.

Each grating is characterized by its periodicity, cross-sectional shape, and its material properties, i.e., the real and imaginary parts of the refractive index ñ. A sinusoidal absorption grating is imaged with degree of partial coherence varying from σ= 0 to 2. Phase sinusoidal gratings with phase variations from φ = 0 to φ = π/3 are examined under σ = 0, 0.1, 0.2, and ∞. Note again that for phase objects, image visibility Vimage is used instead of MTF since the object’s visibility Vobject is zero by definition (see Appendix B). Complex objects with both phase and absorption effects are also examined, noting that such complex objects more realistically resemble samples at the x-ray wavelengths. Performance of the imaging system is examined by varying the magnitude of the object’s phase variation while maintaining a constant absorption, i.e., varying the value of δ/β while keeping the value of η constant. Specifically, the δ/β value is increased from 0 to 10 while the value of η is fixed at 0.1 in order to isolate the effect of phase on image contrast.

3.1. Absorption Sinusoidal Grating/Absorption Objects

We first examine the sinusoidal absorption grating which has only three diffraction orders, i.e., 0 and ±1. Under coherent illumination, the ±1 diffraction orders are allowed to pass the pupil if α < 1. With FZP, the pupil acts as a low-pass filter allowing all harmonics below the cut-off frequency νc to pass without modification. Therefore, within α < 1, the object’s waveform is reconstructed without loss of information on the image plane (Fig. 6(a)) [20]. On the other hand, performance of XOR and ZPD depends on their bias retardation since the ±1 diffraction orders are modified by the patterns within their pupils. Figures 6(b)–6(f) and Figs. 7(b)–7(f) respectively show MTFs obtained using XOR and ZPD under coherent illumination with various bias retardation values. These results show that XOR and ZPD are not optimized for absorption objects since their MTF is in general lower than that of FZP. For XOR, the low MTF region gradually widens as Δθ increases from 0 to π/2, while the low MTF region for ZPD remains fixed for 0 < Δθ < π/2. Note that XOR and ZPD with Δθ = π/2 always result in MTF = 1 for α < 1 because in these cases value of the pupil function at the origin is zero, and the zeroth order is consequently eliminated. As a result, minimum intensity in the image is suppressed to zero and MTF becomes one for α < 1. Since the pupil function of XOR consists of a binary grating while that of ZPD has a cosine structure (Fig. 4), the transitions between low and high MTF regions are steep for XOR and smooth for ZPD.

 figure: Fig. 6

Fig. 6 MTF of the absorption sinusoidal grating object for (a) FZP, (b) XOR with Δθ = 0, (c) XOR Δθ = π/12, (d) XOR Δθ = π/4, (e) XOR Δθ = π/3, and (f) XOR Δθ = π/2. The coherent cases, i.e., σ = 0, are plotted separately and located at the bottom of each panel.

Download Full Size | PDF

 figure: Fig. 7

Fig. 7 MTF of the absorption sinusoidal grating object for (a) FZP, (b) ZPD Δθ = 0, (c) ZPD Δθ = π/12, (d) ZPD Δθ = π/4, (e) ZPD Δθ = π/3, and (f) ZPD Δθ = π/2. The results for Δθ = 0 and Δθ = π/12 are very similar. The coherent cases, i.e., σ = 0, are plotted separately and located at the bottom of each panel.

Download Full Size | PDF

For partially coherent illumination, i.e., σ > 0, non-zero MTF values are observed beyond α = 1 as expected. Again, for absorptive objects performance of the XOR and ZPD are inferior to FZP especially for XOR with large Δθ values (Figs. 6 and 7). Effect of bias retardation on MTF is relatively small when Δθ is smaller than π/12 for ZPD, as can be seen by comparing Figs. 7(b) and 7(c). Overall, DIC microscopy is not suited for imaging absorption objects and the use of a conventional x-ray microscope is recommended for such objects.

3.2. Phase Sinusoidal Grating/Phase Objects

In this section, we use image visibility Vimage to evaluate the performance of DIC microscopes for pure phase objects. The first columns of Figs. 8 and 9 show the image visibilities of a phase sinusoidal grating under coherent illumination. In general, a phase object is nearly invisible with FZP when its phase variation magnitude φ is small. However, such phase grating objects can become visible even with FZP when φ is large, as shown in Fig. 8(a). This is due to the fact that, unlike the absorption sinusoidal grating which has only three diffraction orders, the phase sinusoidal grating has an infinite number of diffraction orders (Table 2). When α is small, i.e., α < 0.5, a large number of diffraction orders pass though the pupil in the image plane. This allows the pure phase object to be reconstructed without loss of information. Consequently, the image, as a higher fidelity representation of the phase object, lacks contrast in intensity. As α increases beyond 0.5, the finite extent of the pupil prevents the passing of higher order harmonics. As a result, only lower order harmonics contribute to the reconstructed image and the image becomes visible at the cost of missing higher order harmonics. However, image visibility in this case remains relatively low (Vimage = 0.43) even for φ = π/3.

 figure: Fig. 8

Fig. 8 Image visibility Vimage of the phase sinusoidal grating object with φ = 0 ∼ π/3 for various σ values. First column shows the coherent case, i.e., σ = 0. (a) FZP, (b) XOR Δθ = 0, (c) XOR Δθ = π/12, (d) XOR Δθ = π/4, (e) XOR Δθ = π/3, and (f) XOR Δθ = π/2.

Download Full Size | PDF

 figure: Fig. 9

Fig. 9 Image visibility Vimage of the phase sinusoidal grating object with φ = 0 ∼ π/3 for various σ values. First column shows the coherent case, i.e., σ = 0. (a) FZP, (b) ZPD Δθ = 0, (c) ZPD Δθ = π/12, (d) ZPD Δθ = π/4, (e) ZPD Δθ = π/3, and (f) ZPD Δθ = π/2.

Download Full Size | PDF

As seen in Figs. 8(b)–8(f) and Figs. 9(b)–9(f), XOR and ZPD under coherent illumination are able to produce large image contrast for phase objects. Bias retardation is essential for the performance of XOR and ZPD. Certainly, given the same phase object, the resultant image visibility varies significantly with different Δθ values and larger Δθ generally results in higher image visibility over a wider range of α. As in the case of absorption objects, the transitions between low and high image visibility regions are steep for XOR and smooth for ZPD. Also note that XOR and ZPD with Δθ = π/2 always yield Vimage = 1 below the cutoff frequency due to the elimination of the zeroth order component.

Image visibility values Vimage of a phase sinusoidal grating under partially coherent illuminations for σ = 0.1, σ = 0.2, and σ = ∞ are also shown in Figs. 8 and 9. The image visibility gradually decreases with increasing σ and completely fades at σ = ∞. The range of high Vimage region depends on bias retardation and in general larger Δθ values produce higher Vimage for a wider range of α. Note that the XOR with Δθ = π/2 has lower visibility overall as shown in Fig. 8(f).

3.3. Complex Sinusoidal Grating

A complex sinusoidal grating is used to examine the effect of absorption and phase concurrently. Results presented in the previous two sections show that XOR and ZPD with proper bias retardation are able to produce relatively large image contrast for sinusoidal phase gratings, whereas FZP is suited for absorption objects. By varying the value of δ/β from 0 to 10 while fixing η at 0.1, MTFs of the XOR and ZPD are examined and compared. Note that for δ/β = 0, the object introduces no phase effect and as a result characteristics of an absorption sinusoidal grating is observed as expected. As the δ/β value increases, the resultant images reflect features of both absorption and phase gratings.

MTFs under coherent illumination for complex sinusoidal objects are shown in the first columns of Figs. 10 and 11. Observed image intensity modulation in the FZP imaging system is mostly due to the absorption of the complex sinusoidal grating and the effect of the phase variation is negligible. As a result, the MTF value of FZP remains around unity even for δ/β = 10 (Fig. 10(a)). Similarly, image intensities obtained using XOR and ZPD with Δθ = 0 barely reflect the presence of phase variations in the complex object, as shown in Figs. 10(b) and 11(b). These observations agree with the results presented in the previous section, i.e., phase effect on the image intensity is generally negligible for XOR and ZPD with Δθ = 0. For 0 < Δθ < π/2, image intensity of the complex object with large δ/β is visibly enhanced by XOR and ZPD. MTF values as high as 8.45 for XOR and 8.15 for ZPD can be obtained for the complex object with δ/β = 10 (Figs. 10(d)–10(f) and Figs. 11(c)–11(f)). Note that the values of MTF can be greater than 1 for complex objects imaged by DIC objectives because DIC image contrast in this case is enhanced by the object’s phase variations. As Δθ increases, high MTF region widens and eventually extends up to the theoretical resolution limit, i.e., α = 1, for both XOR and ZPD. Note that both XOR and ZPD with Δθ = π/2 produce MTF = 10 for all α < 1 (Figs. 10(f) and 11(f)). This is because in these cases the zeroth order Fourier coefficients are eliminated by the pupil function, and the resultant visibility in image intensity becomes 1 for all δ/β values. Since the object’s intensity visibility is Vobject = 0.1 for the complex sinusoidal object with η = 0.1 (see Appendix B), MTF = 1/0.1 = 10 in this case.

 figure: Fig. 10

Fig. 10 MTF of the complex sinusoidal object with η = 0.1 and δ/β = 0 ∼ 10 for various σ values. (a) FZP, (b) XOR Δθ = 0, (c) XOR Δθ = π/12, (d) XOR Δθ = π/4, (e) XOR Δθ = π/3, and (f) XOR Δθ = π/2.

Download Full Size | PDF

 figure: Fig. 11

Fig. 11 MTF of the complex sinusoidal object with η = 0.1 and δ/β = 0 ∼ 10 for various σ values. (a) FZP, (b) ZPD Δθ = 0, (c) ZPD Δθ = π/12, (d) ZPD Δθ = π/4, (e) ZPD Δθ = π/3, and (f) ZPD Δθ = π/2.

Download Full Size | PDF

Modulation transfer functions of complex objects under partially coherent illuminations are shown in Figs. 10 and 11 with σ = 0.1, 0.2, and ∞. As expected, the effect of phase in a complex object is minor for FZP but significant for XOR and ZPD. ZPD and XOR with Δθ = 0 again show similar characteristics as FZP, as seen in Figs. 10(b) and 11(b). Increasing the value of Δθ widens the high MTF region for both DIC objectives.

3.4. Effect of Rotation on MTF

In the previous sections, we have examined the performance of DIC objectives whose orientation along the shear direction is perpendicular to the grating structure, i.e., the relative rotational angle between the DIC objective and the grating object θr = 0° (Fig. 12). Here we examine the effect of rotation on MTF by orienting each DIC objective with respect to the grating objects. Figures 13 and 14 show the results with rotation angles θr = 0° ∼ 90° for various σ values. Note that a complex sinusoidal grating with η = 0.1 and δ/β = 0 ∼ 10 is used as the object and the object spatial frequency α is fixed at 0.6 for all cases. In general, MTF is largest at θr = 0° and minimizes at θr = 90°. When σ = 0, the transitions between low and high MTF regions are again steep for XOR and smooth for ZPD.

 figure: Fig. 12

Fig. 12 DIC objectives, (a) XOR and (b) ZPD, with rotation angle θr. In both examples, Δθ = π/2 and θr = 40° are used.

Download Full Size | PDF

 figure: Fig. 13

Fig. 13 Effect of XOR’s orientation on MTF for various σ values. A complex sinusoidal object with η = 0.1 and δ/β = 0 ∼ 10 is used. XOR is rotated an angle of θr = 0° ∼ 90° with respect to the grating. The normalized object spatial frequency α is fixed at 0.6. (a) XOR Δθ = 0, (b) XOR Δθ = π/12, (c) XOR Δθ = π/4, (d) XOR Δθ = π/3, and (e) XOR Δθ = π/2.

Download Full Size | PDF

 figure: Fig. 14

Fig. 14 Effect of ZPD’s orientation on MTF for various σ values. A complex sinusoidal object with η = 0.1 and δ/β = 0 ∼ 10 is used. ZPD is rotated an angle of θr = 0° ∼ 90° with respect to the grating. The normalized object spatial frequency α is fixed at 0.6. (a) ZPD Δθ = 0, (b) ZPD Δθ = π/12, (c) ZPD Δθ = π/4, (d) ZPD Δθ = π/3, and (e) ZPD Δθ = π/2.

Download Full Size | PDF

4. Conclusion

Performance of x-ray DIC microscopes employing XOR and ZPD objectives was quantitatively examined using MTF. It has been demonstrated that x-ray DIC microscopy can utilize the phase information in the sample to produce large image contrast for both phase and complex objects. Effect of partial coherence was also quantitatively examined and shown that illumination coherence plays an important role in DIC image contrast.

Appendix A Fourier Coefficients of Grating Objects

A.1 Absorption Sinusoidal Grating

Object transmittance function t o(ξ′, η′) of the absorption sinusoidal grating, illustrated in Fig. 5(a), is defined as,

to(ξ,η)=1+cos(2πνoξ).
Its Fourier spectrum is of the form of Eq. (6) with the nth order Fourier coefficients,
an={1,n=01/2,n=±10,otherwise.

A.2 Phase Sinusoidal Grating

Object transmittance function t o(ξ′, η′) of a phase sinusoidal grating shown in Fig. 5(d) is given by [29],

to(ξ,η)=exp[jϕsin(2πνoξ)],
where φ is the magnitude of the phase variation. Using the identity [29]
exp(±jzsinθ)=n=+(±1)nJn(z)ejnθ,
where Jn is the nth order Bessel function of the first kind and denoting θ = 2πνoξ′, Eq. (A.3) can be rewritten as,
to(ξ,η)=n=+Jn(ϕ)exp(j2πnνoξ).
Fourier spectrum of Eq. (A.5) is therefore,
𝕋o(z1,z2)=n=+Jn(ϕ)δ(z1nνo,z2),
where δ is the Dirac delta function and the nth order Fourier coefficient for a phase sinusoidal grating is given by,
an=Jn(ϕ).

A.3 Complex Sinusoidal Grating

Refractive index ñ at x-ray wavelengths is in general complex and has the form ñ = 1 – δ + j β where (1 – δ) determines the phase shift introduced by the material and β is the absorption. Here we assumed that the grating is sinusoidally shaped and its material has a refractive index ñ = 1 – δ + j β. Assuming that the object is placed in vacuum (refractive index no = 1), the object transmittance function t o(ξ′, η′) is given by,

to(ξ,η)=exp[jk(nno)t(ξ,η)],
where k is the vacuum wavenumber, and t(ξ′, η′) describes the grating thickness variation across the object plane (ξ′, η′), i.e.,
t(ξ,η)=to[1+sin(2πνoξ)].
Note that to is the maximum thickness variation in the grating. Denoting the attenuation length as ta = 1/(2), η as the ratio of to to ta, and considering only the spatially varying part of the grating, Eq. (A.8) becomes,
to(ξ,η)=exp[jk(δ+jβ)ηtasin(2πνoξ)]=exp[j(δβ+j)η2sin(2πνoξ)].

Using the identity in Eq. (A.4) and

z=(δβ+j)η2,
Equation (A.10) can be expressed as,
to(ξ,η)=n=+Jn(z)exp(j2πnνoξ).
Note that complex sinusoidal grating is characterized by η and the ratio of δ to β. Figure 5(c) shows the complex sinusoidal grating with η = 1 and δ/β = 1.

Fourier spectrum of the complex sinusoidal transmittance function is therefore,

𝕋o(z1,z2)=n=+Jn(z)δ(z1nνo,z2),
where δ is the Dirac delta function. Fourier coefficients of the complex sinusoidal grating are
an=Jn(z),
where z is given by Eq. (A.11).

Appendix B Visibility of Grating Objects

Visibility V of an intensity pattern is defined as [27],

VImaxIminImax+Imin,
where Imax and Imin are the maximum and the minimum of the intensity distribution, respectively. Here we define the object’s intensity distribution as the absolute square of its transmittance function, i.e., |t o(ξ′, η′)|2. This definition follows the classical paper by Becherer and Parrent [20] and is widely used in MTF definition [29].

B.1 Absorption Sinusoidal Grating

Intensity distribution of the absorption sinusoidal object is given by,

|to(ξ,η)|2=|1+cos(2πνoξ)|2=32+2cos(2πνoξ)+12cos(4πνoξ).
Object’s visibility of the absorption sinusoidal object is therefore 1.

B.2 Phase Sinusoidal Grating

Intensity distribution of the phase sinusoidal object is given by,

|to(ξ,η)|2=|exp[jϕsin(2πνoξ)]|2=1,
and its visibility is therefore zero.

B.3 Complex Sinusoidal Grating

Intensity distribution of the complex sinusoidal object is given by,

|to(ξ,η)|2=|exp[j(δβ+j)η2sin(2πνoξ)]|2=exp[ηsin(2πνoξ)].
Accordingly, this complex sinusoidal object’s visibility is given by,
Vobject=eηeηeη+eη=tanhη.

References and links

1. A. Sakdinawat and D. Attwood, “Nanoscale X-ray imaging,” Nat. Photonics 4, 840–848 (2010). [CrossRef]  

2. W. Chao, J. Kim, S. Rekawa, P. Fischer, and E. H. Anderson, “Demonstration of 12 nm resolution Fresnel zone plate lens based soft X-ray microscopy,” Opt. Express 17, 17669–17677 (2009). [CrossRef]   [PubMed]  

3. S. Rehbein, S. Heim, P. Guttmann, S. Werner, and G. Schneider, “Ultrahigh-resolution soft-X-ray microscopy with zone plates in high orders of diffraction,” Phys. Rev. Lett. 103, 110801 (2009). [CrossRef]   [PubMed]  

4. C. Chang, A. Sakdinawat, P. Fischer, E. Anderson, and D. Attwood, “Single-element objective lens for soft X-ray differential interference contrast microscopy,” Opt. Lett. 31, 1564–1566 (2006). [CrossRef]   [PubMed]  

5. T. Wilhein, B. Kaulich, E. Di Fabrizio, F. Romanato, S. Cabrini, and J. Susini, “Differential interference contrast X-ray microscopy with submicron resolution,” Appl. Phys. Lett. 78, 2082–2084 (2001). [CrossRef]  

6. E. Di Fabrizio, D. Cojoc, S. Cabrini, B. Kaulich, J. Susini, P. Facci, and T. Wilhein, “Diffractive optical elements for differential interference contrast X-ray microscopy,” Opt. Express 11, 2278–2288 (2003). [CrossRef]   [PubMed]  

7. A. Sakdinawat and Y. Liu, “Phase contrast soft X-ray microscopy using Zernike zone plates,” Opt. Express 16, 1559–1564 (2008). [CrossRef]   [PubMed]  

8. O. von Hofsten, M. Bertilson, M. Lindblom, A. Holmberg, and U. Vogt, “Compact Zernike phase contrast X-ray microscopy using a single-element optic,” Opt. Lett. 33, 932–934 (2008). [CrossRef]   [PubMed]  

9. C. Holzner, M. Feser, S. Vogt, B. Hornberger, S. B. Baines, and C. Jacobsen, “Zernike phase contrast in scanning microscopy with X-rays,” Nat. Phys. 6, 883–887 (2010). [CrossRef]  

10. M. Dierolf, A. Menzel, P. Schneider, C. M. Kewish, R. Wepf, O. Bunk, and F. Pfeiffer, “Ptychographic X-ray computed tomography at the nanoscale,” Nature (London) 467, 436–439 (2010). [CrossRef]  

11. K. Giewekemeyer, P. Thibault, S. Kalbfleisch, A. Beerlink, C. M. Kewish, M. Dierolf, F. Pfeiffer, and T. Salditt, “Quantitative biological imaging by ptychographic X-ray diffraction microscopy,” Proc. Natl. Acad. Sci. U.S.A. 107, 529–534 (2010). [CrossRef]  

12. F. Pfeiffer, T. Weitkamp, O. Bunk, and C. David, “Phase retrieval and differential phase-contrast imaging with low-brilliance X-ray sources,” Nat. Phys. 2, 258–261 (2006). [CrossRef]  

13. A. Momose, W. Yashiro, H. Maikusa, and Y. Takeda, “High-speed X-ray phase imaging and X-ray phase tomography with talbot interferometer and white synchrotron radiation,” Opt. Express 17, 12540–12545 (2009). [CrossRef]   [PubMed]  

14. P. Zhu, K. Zhang, Z. Wang, Y. Liu, X. Liu, Z. Wu, S. A. McDonald, F. Marone, and M. Stampanoni, “Low-dose, simple, and fast grating-based X-ray phase-contrast imaging,” Proc. Natl. Acad. Sci. U.S.A. 107, 13576–13581 (2010). [CrossRef]   [PubMed]  

15. U. Vogt, M. Lindblom, P. A. C. Jansson, T. T. Tuohimaa, A. Holmberg, H. M. Hertz, M. Wieland, and T. Wilhein, “Single-optical-element soft-X-ray interferometry with a laser-plasma x-ray source,” Opt. Lett. 30, 2167–2169 (2005). [CrossRef]   [PubMed]  

16. M. C. Bertilson, O. von Hofsten, M. Lindblom, T. Wilhein, H. M. Hertz, and U. Vogt, “Compact high-resolution differential interference contrast soft X-ray microscopy,” Appl. Phys. Lett. 92, 064104 (2008). [CrossRef]  

17. C. Chang, P. Naulleau, E. Anderson, K. Rosfjord, and D. Attwood, “Diffractive optical elements based on Fourier optical techniques: a new class of optics for extreme ultraviolet and soft X-ray wavelengths,” Appl. Opt. 41, 7384–7389 (2002). [CrossRef]   [PubMed]  

18. B. Kaulich, T. Wilhein, E. Di Fabrizio, F. Romanato, M. Altissimo, S. Cabrini, B. Fayard, and J. Susini, “Differential interference contrast X-ray microscopy with twin zone plates,” J. Opt. Soc. Am. A 19, 797–806 (2002). [CrossRef]  

19. T. Nakamura and C. Chang, “Quantitative X-ray differential interference contrast microscopy with independently adjustable bias and shear,” Phys. Rev. A 83, 043808 (2011). [CrossRef]  

20. R. J. Becherer and G. B. Parrent Jr., “Nonlinearity in optical imaging systems,” J. Opt. Soc. Am. 57, 1479–1482 (1967). [CrossRef]  

21. W. Chao, B. D. Harteneck, J. A. Liddle, E. H. Anderson, and D. T. Attwood, “Soft X-ray microscopy at a spatial resolution better than 15nm,” Nature (London) 435, 1210–1213 (2005). [CrossRef]  

22. O. von Hofsten, M. Bertilson, and U. Vogt, “Theoretical development of a high-resolution differential-interference-contrast optic for X-ray microscopy,” Opt. Express 16, 1132–1141 (2008). [CrossRef]   [PubMed]  

23. A. Noguchi, H. Ishiwata, M. Itoh, and T. Yatagai, “Optical sectioning in differential interference contrast microscopy,” Opt. Commun. 282, 3223–3230 (2009). [CrossRef]  

24. Y. Ichioka, K. Yamamoto, and T. Suzuki, “Image of a sinusoidal complex object in a partially coherent optical system,” J. Opt. Soc. Am. 65, 892–902 (1975). [CrossRef]  

25. Y. Ichioka and T. Suzuki, “Image of a periodic complex object in an optical system under partially coherent illumination,” J. Opt. Soc. Am. 66, 921–932 (1976). [CrossRef]  

26. S. Inoué and K. R. Spring, Video Microscopy: The Fundamentals, 2nd ed. (Plenum Press, 1997). [CrossRef]  

27. J. W. Goodman, Statistical Optics (Wiley, 1985).

28. T. Nakamura and C. Chang, “Exact space invariant illumination for partially coherent imaging systems,” J. Opt. Soc. Am. A 27, 1953–1961 (2010). [CrossRef]  

29. J. W. Goodman, Introduction to Fourier Optics, 3rd ed. (Roberts & Company Publishers, 2005).

30. D. Attwood, Soft X-rays and Extreme Ultraviolet Radiation: Principles and Applications (Cambridge University Press, 1999).

31. C. Chang, E. Anderson, P. Naulleau, E. Gullikson, K. Goldberg, and D. Attwood, “Direct measurement of index of refraction in the extreme-ultraviolet wavelength region with a novel interferometer,” Opt. Lett. 27, 1028–1030 (2002). [CrossRef]  

32. C. Preza, D. L. Snyder, and J. A. Conchello, “Theoretical development and experimental evaluation of imaging models for differential-interference-contrast microscopy,” J. Opt. Soc. Am. A 16, 2185–2199 (1999). [CrossRef]  

33. C. Chang and T. Nakamura, “Partially coherent image formation theory for X-ray microscopy,” in Microscopy: Science, Technology, Applications and Education, 4th ed., A. Méndez-Vilas and J. Díaz, eds. (Formatex Research Center, 2010), Vol. 3, pp. 1897–1904.

34. M. Born and E. Wolf, Principles of Optics, 6th ed. (Cambridge University Press, 1997).

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (14)

Fig. 1
Fig. 1 Performance comparison between FZP (conventional), XOR, and ZPD based x-ray microscopes. The complex object (left column) is illuminated coherently and imaged by FZP, XOR, and ZPD objectives, respectively (center column). Note that the object’s intensity is defined as the absolute square of the object transmittance function. The resultant image intensity (right column) shows that XOR and ZPD objectives generate higher image contrast for the complex object. It is clearly seen that image visibility produced by the two DIC objectives exceeds the object’s visibility while the image visibility produced by FZP is at best equal to that of the object. The value of MTF for DIC objectives can therefore exceed 1 (see text for details). The same bias retardation value 2Δθ = π/2 is used for the XOR and ZPD objectives in this figure.
Fig. 2
Fig. 2 Schematic diagram of a partially coherent imaging system.
Fig. 3
Fig. 3 XOR and ZPD based DIC objectives with various bias retardation values. (a) XOR with Δθ = 0, (b) XOR with Δθ = π/4, (c) XOR with Δθ = π/2, (d) ZPD with Δθ = 0, (e) ZPD with Δθ = π/4, and (f) ZPD with Δθ = π/2.
Fig. 4
Fig. 4 Pupil functions corresponding to the objectives in Fig. 3. (a) XOR with Δθ = 0, (b) XOR with Δθ = π/4, (c) XOR with Δθ = π/2, (d) ZPD with Δθ = 0, (e) ZPD with Δθ = π/4, and (f) ZPD with Δθ = π/2.
Fig. 5
Fig. 5 Absorption (top) and phase (bottom) variations of (a) absorption sinusoidal, (b) phase sinusoidal, and (d) complex sinusoidal grating objects. Period of the grating is denoted by T.
Fig. 6
Fig. 6 MTF of the absorption sinusoidal grating object for (a) FZP, (b) XOR with Δθ = 0, (c) XOR Δθ = π/12, (d) XOR Δθ = π/4, (e) XOR Δθ = π/3, and (f) XOR Δθ = π/2. The coherent cases, i.e., σ = 0, are plotted separately and located at the bottom of each panel.
Fig. 7
Fig. 7 MTF of the absorption sinusoidal grating object for (a) FZP, (b) ZPD Δθ = 0, (c) ZPD Δθ = π/12, (d) ZPD Δθ = π/4, (e) ZPD Δθ = π/3, and (f) ZPD Δθ = π/2. The results for Δθ = 0 and Δθ = π/12 are very similar. The coherent cases, i.e., σ = 0, are plotted separately and located at the bottom of each panel.
Fig. 8
Fig. 8 Image visibility Vimage of the phase sinusoidal grating object with φ = 0 ∼ π/3 for various σ values. First column shows the coherent case, i.e., σ = 0. (a) FZP, (b) XOR Δθ = 0, (c) XOR Δθ = π/12, (d) XOR Δθ = π/4, (e) XOR Δθ = π/3, and (f) XOR Δθ = π/2.
Fig. 9
Fig. 9 Image visibility Vimage of the phase sinusoidal grating object with φ = 0 ∼ π/3 for various σ values. First column shows the coherent case, i.e., σ = 0. (a) FZP, (b) ZPD Δθ = 0, (c) ZPD Δθ = π/12, (d) ZPD Δθ = π/4, (e) ZPD Δθ = π/3, and (f) ZPD Δθ = π/2.
Fig. 10
Fig. 10 MTF of the complex sinusoidal object with η = 0.1 and δ/β = 0 ∼ 10 for various σ values. (a) FZP, (b) XOR Δθ = 0, (c) XOR Δθ = π/12, (d) XOR Δθ = π/4, (e) XOR Δθ = π/3, and (f) XOR Δθ = π/2.
Fig. 11
Fig. 11 MTF of the complex sinusoidal object with η = 0.1 and δ/β = 0 ∼ 10 for various σ values. (a) FZP, (b) ZPD Δθ = 0, (c) ZPD Δθ = π/12, (d) ZPD Δθ = π/4, (e) ZPD Δθ = π/3, and (f) ZPD Δθ = π/2.
Fig. 12
Fig. 12 DIC objectives, (a) XOR and (b) ZPD, with rotation angle θr . In both examples, Δθ = π/2 and θr = 40° are used.
Fig. 13
Fig. 13 Effect of XOR’s orientation on MTF for various σ values. A complex sinusoidal object with η = 0.1 and δ/β = 0 ∼ 10 is used. XOR is rotated an angle of θr = 0° ∼ 90° with respect to the grating. The normalized object spatial frequency α is fixed at 0.6. (a) XOR Δθ = 0, (b) XOR Δθ = π/12, (c) XOR Δθ = π/4, (d) XOR Δθ = π/3, and (e) XOR Δθ = π/2.
Fig. 14
Fig. 14 Effect of ZPD’s orientation on MTF for various σ values. A complex sinusoidal object with η = 0.1 and δ/β = 0 ∼ 10 is used. ZPD is rotated an angle of θr = 0° ∼ 90° with respect to the grating. The normalized object spatial frequency α is fixed at 0.6. (a) ZPD Δθ = 0, (b) ZPD Δθ = π/12, (c) ZPD Δθ = π/4, (d) ZPD Δθ = π/3, and (e) ZPD Δθ = π/2.

Tables (2)

Tables Icon

Table 1 Pupil Functions P(x,y) of FZP, XOR, and ZPD Objectives

Tables Icon

Table 2 Transmittance Functions and Fourier Coefficients of the Gratings Used in the Text *

Equations (32)

Equations on this page are rendered with MathJax. Learn more.

I i ( u , v ) = + d p d q 𝕁 o ( p , q ) × + d w 1 d w 2 𝕋 o ( w 1 , w 2 ) 𝕂 ( w 1 p , w 2 q ) × + d ν U d ν V exp { j 2 π ( u ν U + v ν V ) } × 𝕋 o * ( w 1 ν U , w 2 ν V ) 𝕂 * ( w 1 p ν U , w 2 q ν V ) ,
I s ( α , β ) = { I o , α 2 + β 2 < R s 0 , otherwise ,
𝕁 o ( p , q ) = κ M 2 I s ( λ f c M p , λ f c M q ) ,
𝕂 ( p , q ) = M P ( λ z o M p , λ z o M q ) ,
t o ( ξ , η ) = n = + a n exp ( j 2 π n ν o ξ ) ,
𝕋 o ( z 1 , z 2 ) = n = + a n δ ( z 1 n ν o , z 2 ) ,
I i ( u , v ) = m = + n = + a n a m * exp { j 2 π [ ( n m ) ν o u ] } × + d p d q 𝕁 o ( p , q ) 𝕂 ( p + n ν o , q ) 𝕂 * ( p + m ν o , q ) .
I i ( u , v ) = m = + n = + a n a m * exp { j 2 π [ ( n m ) ν o u ] } × κ M 2 + I s ( λ f c M p , λ f c M q ) P ( λ z o ( M p + M n ν o ) , λ z o M q ) × P * ( λ z o ( M p + M m ν o ) , λ z o M q ) d ( M p ) d ( M q ) ,
ν c R o λ z o
ν s R s λ f c ,
σ ν s ν c .
α M ν o ν c ,
M T F ( α ) V image ( α ) V object ( α ) ,
t o ( ξ , η ) = 1 + cos ( 2 π ν o ξ ) .
a n = { 1 , n = 0 1 / 2 , n = ± 1 0 , otherwise .
t o ( ξ , η ) = exp [ j ϕ sin ( 2 π ν o ξ ) ] ,
exp ( ± j z sin θ ) = n = + ( ± 1 ) n J n ( z ) e j n θ ,
t o ( ξ , η ) = n = + J n ( ϕ ) exp ( j 2 π n ν o ξ ) .
𝕋 o ( z 1 , z 2 ) = n = + J n ( ϕ ) δ ( z 1 n ν o , z 2 ) ,
a n = J n ( ϕ ) .
t o ( ξ , η ) = exp [ j k ( n n o ) t ( ξ , η ) ] ,
t ( ξ , η ) = t o [ 1 + sin ( 2 π ν o ξ ) ] .
t o ( ξ , η ) = exp [ j k ( δ + j β ) η t a sin ( 2 π ν o ξ ) ] = exp [ j ( δ β + j ) η 2 sin ( 2 π ν o ξ ) ] .
z = ( δ β + j ) η 2 ,
t o ( ξ , η ) = n = + J n ( z ) exp ( j 2 π n ν o ξ ) .
𝕋 o ( z 1 , z 2 ) = n = + J n ( z ) δ ( z 1 n ν o , z 2 ) ,
a n = J n ( z ) ,
V I m a x I m i n I m a x + I m i n ,
| t o ( ξ , η ) | 2 = | 1 + cos ( 2 π ν o ξ ) | 2 = 3 2 + 2 cos ( 2 π ν o ξ ) + 1 2 cos ( 4 π ν o ξ ) .
| t o ( ξ , η ) | 2 = | exp [ j ϕ sin ( 2 π ν o ξ ) ] | 2 = 1 ,
| t o ( ξ , η ) | 2 = | exp [ j ( δ β + j ) η 2 sin ( 2 π ν o ξ ) ] | 2 = exp [ η sin ( 2 π ν o ξ ) ] .
V object = e η e η e η + e η = tanh η .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.