X-ray imaging inside the focal spot of polycapillary optics using the coded aperture concept

We present a method for x-ray imaging of objects inside the focal spot of polycapillary optics that resolves details smaller than the focal spot dimensions. This method employs coded aperture imaging, in which the micro-structure of polycapillary optics is treated as the coding pattern. Projection of the object is decoded from a magnified x-ray image of the polycapillary structure which is specifically sharpened by the object. Field of view can be extended by scanning the object across the focal spot. © 2013 Optical Society of America OCIS codes: (340.0340) X-ray optics; (180.7460) X-ray microscopy. References and links 1. M. Kumakhov, “Channeling of photons and new x-ray optics,” Nucl. Instrum. Methods B 48, 283–286 (1990). 2. M. A. Kumakhov and F. F. Komarov, “Multiple reflection from surface x-ray optics,” Phys. Rep. 191, 289–350 (1990). 3. S. Dabagov, “Channeling of neutral particles in micro-and nanocapillaries,” Phys.-Usp. 46, 1053– 1075 (2003). 4. S. A. Hoffman, D. J. Thiel, and D. H. Bilderback, “Developments in tapered monocapillary and polycapillary glass X-ray concentrators ,” Nucl. Instrum. Meth. A 347, 384–389 (1994). 5. L. Vincze, K. Janssens, F. Adams, A. Rindby, and P .Engström, “Interpretation of capillary generated spatial and angular distributions of x rays: Theoretical modeling and experimental verification using the European Synchrotron Radiation Facility Optical beam line,” Rev. Sci. Instr. 69, 3494–3503 (1998). 6. A. Bjeoumikhov, “Observation of peculiarities in angular distributions of x-ray radiation after propagation through polycapillary structures,” Phys. Lett. A 360, 405–410 (2007). 7. S. B. Dabagov, M. A. Kumakhov, and S. V. Nikitina, “On the interference of X-rays in multiple reflection optics,” Phys. Lett. A 203, 279–282 (1995). 8. S. B. Dabagov, M. A. Kumakhov, S. V. Nikitina, V. A. Murashova, R. V. Fedorchuk, and M. N. Yakimenko, “Observation of Interference Effects at the Focus of an X-ray Lens,” J. Synchrotron Radiat. 2, 132–135 (1995). 9. C. MacDonald, “Focusing polycapillary optics and their applications,” X-ray Optics and Instr. 2010, 867049 (2010). 10. N. Gao, I. Y. Ponomarev, Q. F. Xiao, W. M. Gibson, and D. A. Carpenter, “Monolithic polycapillary focusing optics and their applications in microbeam x-ray fluorescence,” Appl. Phys. Lett. 69, 1529–1531 (1996). 11. L. Vincze, B. Vekemans, F. Brenker, G. Falkenberg, K. Rickers, A. Somogyi, M. Kersten, and F. Adams, “Three-dimensional trace element analysis by confocal X-ray microfluorescence imaging,” Anal. Chem. 76, 6786–6791 (2004). 12. B. Kanngiesser, W. Malzer, A. Rodriguez, and I. Reiche, “Three-dimensional micro-XRF investigations of paint layers with a tabletop setup,” Spectroc. Acta Pt. B-Atom. Spectr. 60, 41–47 (2005). 13. V. Chernik and A. Romanov, “X-ray polycapillary imaging microscopy,” Proc. SPIE 5943, 158– 166 (2005). 14. J. G. Ables, “Fourier transform photography: a new method for X-ray astronomy,” Proc. ASA 1, 172–173 (1968). 15. E. E. Fenimore and T. M. Cannon, “Coded aperture imaging with uniformly redundant arrays,”, Appl. Opt. 17, 337–347 (1978). 16. K. A. Nugent, “Coded aperture imaging: a Fourier space analysis ,”Appl. Opt. 26, 563-569 (1987). 17. A. Haboub, A. A. MacDowell, S. Marchesini, and D. Y. Parkinson, “Coded aperture imaging for fluorescent x-rays,” Proc. SPIE 8502, 850209 (2012). 18. L. D. Caro, C. Giannini, A. Cedola, S. Lagomarsino, and I. Bukreeva, “X-ray pointand lineprojection microscopy and diffraction,” Opt. Commun. 265, 8–28 (2006). 19. D. Sayre, “Some implications of a theorem due to Shannon,” Acta Cryst. 5, 843–843 (1952). 20. P. Thibault, M. Dierolf, A. Menzel, O. Bunk, C. David, and F. Pfeiffer, “High-resolution scanning x-ray diffraction microscopy,” Science 321, 379–382 (2008). 21. A. M. Zysk, R. W. Schoonover, Q. Xu, and M. A. Anastasio, “Framework for computing the spatial coherence effects of polycapillary x-ray optics,” Opt. Express 20, 3975–3982 (2012). 22. T. Sun, M. Zhang, Z. Liu, Z. Zhang, G. Li, Y. Ma, X. Du, Q. Jia, Y. Chen, Q. Yuan, W. Huang, P. Zhu, and X. Ding, “Focusing synchrotron radiation using a polycapillary half-focusing X-ray lens for imaging,” J. Synchrotron Radiat. 16, 116–118 (2009). 23. D. Hampai, S. B. Dabagov, G. D. Ventura, F. Bellatreccia, M. Magi, F. Bonfigli, and R. M. Montereali, “High-resolution x-ray imaging by polycapillary optics and lithium fluoride detectors combination,” EPL (Europhys. Lett.) 96, 60010 (2011). 24. A. Kuehn, O. Scharf, I. Ordavo, H. Riesemeier, U. Reinholz, M. Radtke, A. Berger, M. Ostermann, and U. Panne, “Pushing the limits for fast spatially resolved elemental distribution patterns,” J. Anal. At. Spectrom. 26, 1986–1989 (2011). 25. A. G. Peele, K. A. Nugent, A. V. Rode, K. Gabel, M. C. Richardson, R. Strack, and W. Siegmund, “X-ray focusing with lobster-eye optics: a comparison of theory with experiment,” Appl. Opt. 35, 4420–4425 (1996). 26. A. Jarre, C. Fuhse, C. Ollinger, J. Seeger, R. Tucoulou, and T. Salditt, “Two-dimensional hard x-ray beam compression by combined focusing and waveguide optics,”Phys. Rev. Lett. 94, 074801 (2005).


Introduction
Focusing polycapillary devices are reflective x-ray optical elements composed of arrays of glass micro-capillaries [1,2,3].Radiation in individual capillaries propagates by means of multiple total external reflections.A special bending of capillaries enables to concentrate x-rays into an intense micro-spot, hereafter referred to as the focal spot of a polycapillary optics.It is worth to emphasize, that while polycapillary elements are frequently referred to as "lenses" they cannot be described in terms of classical imaging optics [4].For polycapillary optics, the formation of the focal spot can be usually described using geometrical optics as a result of the overlap of beams emitted from individual capillaries [5,6].The lateral dimensions of the focal spot are thus of the order of ∼2 f θ c , where θ c is the critical angle for the total external reflection (θ c [mrad] ≈ 30/E[keV]) and f is the focal length determined by the bending radius of capillaries.For low x-ray energies and/or small capillary diameters, interference effects can introduce fine structures to the shape of the focal spot [7,8].
For a recent review of polycapillary optics discussing applications in spectrometry, imaging and diffraction see [9].Most frequently, focusing polycapillary optics is used in scanning geometry for element specific imaging [10].The use of two crossed lenses permits a 3D confocal x-ray microscopy [11] with a table-top setup [12].In such geometries, imaging resolution is usually defined as the size of the finest focal spot that can be formed using the optics and is typically at the level of 5-50 μm.In contact geometries the resolution is limited to several widths of the capillary channels and also has a value at the level of several microns [13].
In this work we present a method to resolve details of an object located inside the focal spot of polycapillary x-ray optics, which is based on the coded aperture principle [14,15,16,17].The proposed approach works for hard and polychromatic x-rays.In the proof-of-principle experiment we obtained a five-fold improvement in the resolution as compared to the size of the focal spot.

Principle of x-ray imaging inside the focal spot of polycapillary optics
The principle of x-ray imaging inside the focal spot of polycapillary optics is presented in Fig. 1(a).For its description we introduce a simplified model based on geometrical optics.For brevity we restrict the calculation to the xz plane.Consider an extended source at z = 0 (exit surface of lens) having a brightness S(x)G(α), an object located at z = d and x-ray camera placed at z = D. S(x) describes the internal structure of the polycapillary optics [c.f.Fig. 1(b)] and the transmission of the lens.G(α) defines the angular distribution of the radiation emitted from each point of optics.Since each capillary emits radiation in a shallow cone directed towards the focal spot, we assume that G(α) = exp(−0.5α 2/σ 2 ), where σ ∼ θ c , and α is the angle between vectors r − r s and f − r s , where r = (x, D) is the observation point, r s = (x s , 0) runs over the exit surface of the lens and f = (0, f ) is the position of the focal spot.In the paraxial approximation α ≈ (x − Mx s )/D, where M = ( f − D)/ f .The intensity at the detector can be written as: where T (x) = exp[−P(x)] describes the transmission of the object and P(x) = μ(x, z)dz is the projection of the object's absorption coefficient.Convolution of I(x) with the point-spread function of the the detector should be also taken into account in Eq. ( 1).An image recorded without an object will be denoted as I 0 .I 0 evaluated for D = f gives the intensity distribution at the focal spot: a Gaussian with a FWHM of 2.35 f σ .While the information about the structure of S(x) is lost in the focal plane [c.f.Fig. 1(c)], an xray camera placed at position D > f records a magnified (by a factor of M) and blurred image of the internal structure of the polycapillary element.If an object is placed at d > f then, in a crude approximation, the focal spot can be treated as a secondary source producing a magnified point-projection of the object.The imaging resolution in such a case is limited by penumbra blur resulting from the size of the focal spot [18].
The most interesting situation arises when an object is placed inside the focal spot (d = f ).In such a case Eq. ( 1) becomes: From Eq. ( 2) it follows that the shadow of the object vanishes from the camera.Instead of this, the image of the capillary structure changes drastically.Since the object's details introduce high frequencies to the convolution, the image I is unblurred as compared to I 0 in a similar way as in a pinhole camera.What's most important is that Eq. ( 2) describes coded aperture imaging [15,16,17], where the capillary pattern S defines the coded aperture.In fact, in our notation, S combines information about geometrical micro-structure of the polycapillary and its x-ray transmission characteristics.

Experimental setup
Imaging inside the focal spot of polycapillary optics was realized using polychromatic radiation from an x-ray tube with a W anode and a 75 μm spot (Oxford Instruments, XTF 5011) operating at 1 mA and 50 keV, which was focused by a polycapillary lens (IfG).The optics had an input focal length of 40 mm, output focal length of 14 mm, length of 100 mm and exit surface diameter of 2.4 mm.The acceptance cone of the optics had an opening angle of 7 • .The optical microscope image of the lens exit surface [Fig.1(b)] reveals a hexagonal mesh of individual capillaries (diameter 3.5 μm) as well as a coarser one corresponding to bundles of capillaries.The focal spot of the lens was featureless and had an approximate Gaussian shape with FWHM of (40 ± 1) μm as demonstrated in Fig. 1(c).This image was obtained by scanning the focal spot with a 5 μm pinhole and measuring the transmitted intensity.The total flux in the focal spot was approx.10 8 photon/s and the mean energy of the beam was 22 keV.The intensity gain in the focal spot was at the level of 5×10 3 as compared to the intensity of the direct beam through 5 μm pinhole.X-ray images were acquired using a CMOS camera (Rad-Icon, RadEye1).The sensor had 1024×512 pixels separated by 48 μm and was optically coupled to a scintillator.The camera was placed at a distance D = 171 mm from the lens.Since θ c changes with energy (a slow 1/E dependence), for polychromatic radiation the effective focal spot size depends on energetic response of the detector.Therefore, to maintain constant experimental conditions, auxiliary x-ray transmission scans were performed with the same x-ray camera but used as a point-detector just by integrating the intensity over all camera pixels.

X-ray images of simple objects placed inside the focal spot
Let us first demonstrate experimental images obtained for simple objects placed inside the focal spot of a polycapillary lens.Figure 2  Conversely, for x-ray images recorded with objects inside the focal spot, the hexagonal polycapillary mesh is clearly visible despite normalization by I 0 .This means that the object influences this structure or, in other words, the information about the object is encoded in the entire image.Let us discuss the most important features of images from Fig. 2. First, the visible hexagonal structure corresponds to bundles of capillaries [c.f.Fig. 1(b)] and individual capillaries are not resolved.Hence, in this pilot experiment, capillaries determine the behavior of the G function, but the S function is determined by the structure of capillary bundles.Second, the object inside the focal spot produces a sharpening of the images.The degree of sharpening is higher for the smaller object.For the wire, sharpening takes place only along the horizontal direction.These observations are confirmed by the analysis of Fourier transforms presented in Fig. 2(b).Furthermore, a clear contrast reversal in images recorded for the wire and for the pinhole can be observed.

Reconstruction of the object's projection
Next, let us describe the process of decoding or in other words reconstruction of the object's projection.Since this procedure is based on Eq. ( 2), the knowledge of the S(x) function characterizing the capillary structure i.e. the coded aperture is essential.This function has to be known only up to a maximal frequency of the magnified image that can be resolved by the detector.Therefore, S(x) can be approximated by an image I P recorded for a pinhole placed inside the focal spot.In our case, x-ray camera ultimately cuts-off frequencies above ∼ (10μm) −1 .Therefore I P recorded for a 5 μm pinhole, that is shown in Fig. 2, is a sufficient approximation of S(x).
The decoding was performed by a direct implementation of the convolution theorem: For thin objects, U(x) approximates P(x)G(x).Fourier transforms from Fig. 2(b) show discrete peaks corresponding to the periodicity of the capillary bundles.The deconvolution was limited to this discrete set of values which extremely reduced noise.However, let us show that the reconstruction using such a discrete set is at all possible.Shannon sampling theorem [19] says that if the object has a finite support (−a/2, a/2) then its shape can be completely reconstructed from Fourier transform sampled at least at a critical frequency ν C = 1/a.In our case, the signal from an arbitrary large object originates from the area of the focal spot.In practice it is sufficient to assume that the focal spot is confined to a diameter of 2 × FWHM.The corresponding critical frequency ν C is marked in Fig. 2(b) and is higher than the spatial frequency of the 1-st order peaks.This denotes a proper sampling.
In order to test the decoding procedure, coded aperture images were measured for two golden grids with different meshes (600 and 1000 mesh grids).Results of the reconstruction are presented in Fig. 3(a).The mesh structure is clearly resolved inside the area of the focal spot.This proves the possibility of imaging objects inside the focal spot of a polycapillary lens.However, the quality of images is only moderate.Moreover, the increase of the acquisition time did not significantly improve the quality.
The improvement of quality and, simultaneously, the possibility of imaging of larger objects can be realized by an extension of the reconstruction procedure.X-ray coded aperture images can be recorded for the object scanned across the focal spot in such a way that there is an overlap between neighboring decoded images i.e: where Δx and Δy are smaller than 4 f σ .This procedure is very loosely related to the x-ray ptychographic procedure [20].This simple operation improves the data quality.
It corrects the image for the shape of the focal spot and for the asymmetry of the deconvolution kernel as shown in Figs.3(b).A reconstruction of a more complex object (number 12 from a copper electron microscopy finder grid) is presented in Fig. 4(a).The reconstruction was performed from  4 show a drastic improvement in the resolution of the coded aperture imaging as compared to conventional scanning.Features separated by 10 μm are very well resolved.Features separated by 8 μm are at the limit of resolution.The corresponding frequency of (16 μm) −1 is marked in Fig. 2(b).This denotes a five-fold improvement as compared to the transmission scan for which the resolution is limited by focal spot size to 40 μm.
The resolution of the projection from Fig. 4(c) is slightly better than for the scan from Fig. 4(b).In projection geometry the penumbral blur is reduced by a factor of (D − d)/(D − f ) ≈ 0.86 as compared to the focal spot dimensions.In fact, for d → D this factor tends to zero.However, the magnification of the point projection is given by (D − f )/(d − f ) and tends to unity.Hence, for large d the resolution of the detector becomes the dominant factor.This is exemplified in Fig. 4(c) where individual pixels of the detector can be observed.It means that coded aperture imaging gives optimal spatial resolution for a given polycapillary and detector combination.
Finally, we will discuss the limits of the resolution.The fundamental limit is the highest spatial frequency of the S(x) function describing the coded aperture.In practice, this frequency is set by the highest resolved harmonics of the capillary structure.Commercially available lenses have capillary channels with diameter smaller than 1 μm.In order to resolve such a frequency one has to obtain a sufficiently high magnification M. For example, a standard polycapillary optics having f = 5 mm and a 1k×1k pixel detector with a resolution of 50 μm placed at D = 300 mm could provide resolution better than 1 μm.Precise measurement of S(x) with a correction for directly transmitted beam seems to be possible for pinholes with such a diameter.

Conclusions
We presented a method for resolving details of objects inside the focal spot of polycapillary optics which treats the micro-structure of polycapillary optics as a coded aperture.The proof-of-principle experiment was performed using a low-resolution equipment and the formalism used to model the image formation and reconstruction was largely simplified.After improvements, resolution at the level of 1 μm could be possible with commercially available optics.Further developments in x-ray optics fabrication and use of more exact theories of x-ray propagation in capillaries [3,21] could help to obtain sub-micron resolution.The presented results can have impact on other imaging geometries using polycapillaries [13,22,23,24].Presumably, the proposed approach could be used with other kinds of micro or nano-structured optics [25,26].

Fig. 1 .
Fig. 1.Scheme of an experiment for x-ray imaging inside the focal spot of a polycapillary optics using the coded aperture principle.(a) The internal structure of the polycapillary optics is treated as a coding aperture.Projection of the object is decoded from a magnified x-ray image of the polycapillary structure which is specifically sharpened by the object placed at the focal spot.(b) Optical microscope image of the exit surface of the lens.The marked hexagon corresponds to a single capillary bundle.(c) Transmission scan of the focal spot with a 5 μm pinhole.Points show a profile through the center and the line is a Gaussian fit with a FWHM of 40 μm.
(a) shows x-ray images recorded without an object (I 0 ) and for two simple objects: a 5 μm pinhole in a tungsten foil and a 25 μm tungsten wire.For a better presentation ratios I/I 0 are shown.For the pinhole, the image was corrected by a subtraction of an image recorded with the pinhole moved aside of the focal spot.This eliminated the influence of high-energy radiation passing through the foil.The upper insets in the images show projections recorded for objects placed at d = f + 20 mm.The regions in images outside the shadow of the objects are almost uniform which is due to normalization by I 0 .

Fig. 2 .
Fig. 2. X-ray imaging of simple objects placed inside the focal spot.(a) Left: image recorded without object (I 0 ).Middle: data for a 5 μm pinhole.Right: data for a 25 μm wire.Pinhole and wire images were divided by I 0 .Insets at the bottom show zooms of the data.The zoom of I 0 was flattened for a better comparison.The upper insets show projection-like images obtained for objects displaced from focal spot along z axis.(b) Corresponding logarithms of the Fourier transform intensities.Numeric labels provide spatial frequencies in the focal plane.

Fig. 3 .
Fig. 3. Reconstruction of coded aperture images recorded for golden grids with different meshes placed inside the focal spot of a polycapillary optics.(a) Reconstruction from a single image.Left: pitch 42 μm, hole 37 μm.Right: pitch 25 μm, hole 19 μm.(b) Corresponding reconstructions from a set of images.The number of images and their acquisition times are given at the bottom of the plots.

Fig. 4 .
Fig. 4. Comparison of coded aperture imaging with a standard x-ray transmission scan and a projection image.(a) Reconstruction from a set of coded aperture images.The line plot is a profile at the position marked with a dashed line.(b) Transmission scan for d = f .(c) Projection image (I/I 0 ) for d = f + 20mm.(d) Optical microscope image of the object.