Home | About | Journals | Submit | Contact Us | Français |

**|**Opt Express**|**PMC3482909

Formats

Article sections

Authors

Related links

Optics Express

Opt Express. 2012 February 13; 20(4): 3975–3982.

Published online 2012 February 2. doi: 10.1364/OE.20.003975

PMCID: PMC3482909

Received 2011 October 17; Revised 2012 January 13; Accepted 2012 January 19.

Copyright © 2012 Optical Society of America

This is an open-access article distributed under the terms of the Creative Commons Attribution-Noncommercial-No Derivative Works 3.0 Unported License, which permits download and redistribution, provided that the original work is properly cited. This license restricts the article from being modified or used commercially.

Despite the extensive use of polycapillary x-ray optics for focusing and collimating applications, there remains a significant need for characterization of the coherence properties of the output wavefield. In this work, we present the first quantitative computational method for calculation of the spatial coherence effects of polycapillary x-ray optical devices. This method employs the coherent mode decomposition of an extended x-ray source, geometric optical propagation of individual wavefield modes through a polycapillary device, output wavefield calculation by ray data resampling onto a uniform grid, and the calculation of spatial coherence properties by way of the spectral degree of coherence.

The use of x-ray optics in imaging studies has taken on a vital role in the synchrotron and laboratory communities [1, 2]. Polycapillary devices, also known as Kumakhov lenses, consist of collections of glass capillaries that direct x-rays by total external reflection to focus or collimate the output. Coherence effects are well known to exist in capillary structures, as a number of wave theory-based investigations have demonstrated the simulation of x-ray propagation through single capillary devices, usually with synchrotron or soft x-ray laser illumination, and the presence of diffraction and interference effects in the output wavefield [3–6]. As the capillary device size is typically large relative to the wavelength, these simulations have generally relied on geometrical optics or Fresnel-Kirchoff approximations [7–9].

Many simulation studies of polycapillary devices [10] are also performed with ray-tracing toolboxes such as SHADOW [11, 12] and PolyCAD [13, 14]. These simulations function by computing the paths of x-rays through a device and evaluating the number of rays incident on each pixel in the output plane as an estimate of the output intensity. Importantly, the intensity estimate found with these techniques is not directly derived from the underlying wavefield and therefore neglects intensity modulations due to interference effects.

Although diffraction and interference effects in single capillary devices have been extensively investigated experimentally [15–18] and in simulation, relatively little is known of the changes in spatial coherence brought about by propagation through polycapillary devices, especially in the presence of extended laboratory sources [19, 20]. Access to the wavefield data after propagation through the device is necessary for analysis of these effects. In this paper, we present the first demonstration of a method to simulate the wavefield data at the output of a polycapillary device due to an extended source and analyze the resulting spatial coherence properties without the need for an analytically tractable wave-based representation of the optic.

Statistical optics, or coherence theory, describes observable phenomena for non-deterministic optical wavefields [21]. The coherence properties of wavefields generated by a statistically stationary optical source, *Ũ _{s}*(

$${W}_{s}({\mathbf{r}}_{1},{\mathbf{r}}_{2},\omega )=\sum _{n=1}^{\infty}{\beta}_{n}\left(\omega \right){\varphi}_{n}^{*}\left({\mathbf{r}}_{1},\omega \right){\varphi}_{n}\left({\mathbf{r}}_{2},\omega \right),$$

(1)

where ^{*} indicates complex conjugation. Here, the cross-spectral density is equal to the incoherent sum of statistically independent modes, *ϕ _{n}*(

It can be readily shown that the cross-spectral density for a stochastic wavefield *Ũ*(**r**,*ω*), which is the propagated wavefield *Ũ _{s}*(

$$\begin{array}{lll}W\left({\mathbf{r}}_{1},{\mathbf{r}}_{2},\omega \right)& =& {\tilde{U}}^{*}\left({\mathbf{r}}_{1},\omega \right)\tilde{U}\left({\mathbf{r}}_{2},\omega \right)& =& \sum _{n=1}^{\infty}{\beta}_{n}\left(\omega \right){\psi}_{n}^{*}\left({\mathbf{r}}_{1},\omega \right){\psi}_{n}\left({\mathbf{r}}_{2},\omega \right),\end{array}$$

(2)

where *ψ _{n}*(

Finally, the spectral degree of coherence, which quantifies the degree of statistical similarity between the wavefield at two points in space, is defined

$$\mu \left({\mathbf{r}}_{1},{\mathbf{r}}_{2},\omega \right)=\frac{W\left({\mathbf{r}}_{1},{\mathbf{r}}_{2},\omega \right)}{\sqrt{W\left({\mathbf{r}}_{1},{\mathbf{r}}_{1},\omega \right)W\left({\mathbf{r}}_{2},{\mathbf{r}}_{2},\omega \right)}}.$$

(3)

When |*μ*(**r**_{1}, **r**_{2}, *ω*)| = 1, the wavefield is fully coherent between the two points. When |*μ*(**r**_{1}, **r**_{2}, *ω*)| = 0, the wavefield is completely incoherent between the two points.

In this work, we introduce and demonstrate a method for the calculation of the spectral degree of coherence for a wavefield that has propagated through a polycapillary device. This method consists of the following steps: (1) The spatially partially-coherent anode source is decomposed into a collection of statistically-independent point sources that are interpreted as the modes *ϕ _{n}*(

In this section, we show an implementation of this method in which an x-ray tube source and collimating polycapillary optic are modeled. In step 1, an x-ray tube source anode is modeled by a collection of incoherent point emitters with a Gaussian spatial profile. In the CMD representation, the modes and weights are of the form *ϕ _{n}*(

In step 2, interaction of individual wavefields, which correspond to coherent modes, with the polycapillary optic is modeled using geometrical optical propagation. In geometrical optics (also known as the high-frequency limit), the propagating wavefield is assumed to be of the form *U*(**r**,*ω*) = *A*(**r**)exp[*iωS*(**r**)/*c*], where *A*(**r**), the amplitude, and *S*(**r**), the eikonal or optical path length, are independent of frequency, and *c* is the speed of light.

The polycapillary optic modeled in this work is a collimating optic consisting of 19 curved glass capillaries arranged in a hexagonal pattern on the input and output surfaces. The interior wall of a given capillary is defined by
${\left[x-{x}_{l}{g}_{l}(z)\right]}^{2}+{\left[y-{y}_{l}{g}_{l}\left(z\right)\right]}^{2}={r}_{0}^{2}{f}_{l}^{2}\left(z\right)$, where (*x _{l}*,

A ray-tracing algorithm finds the capillary, if any, into which each ray enters. Each entering ray has a straight-line trajectory until it intersects a capillary wall at which point Snell’s law is used to calculate the new trajectory and the amplitude reflection coefficient *R*. After each reflection, the ray amplitude is updated *A _{m}* →

In step 3, the wavefield at the polycapillary output is computed by combining and resampling the discrete ray data onto a uniform grid. A number of methods exist to accomplish this task, including Gaussian beam weighting methods from the optical physics community, which have been used to account for the non-uniformity of ray data [29–31] in simulations of a number of system geometries [32,33]. We have chosen to implement a computationally-efficient weighted beam gridding method, variants of which have been used in the computed tomography community for treatment of non-uniformly sampled spatial Fourier-domain data [34]. In this simulation, the gridding operation is performed by convolving ray point data from each capillary with a Kaiser-Bessel kernel and then sampling on a uniform grid [34].

Specifically, the amplitude and eikonal ray data from each capillary are grouped by the total number of reflections. The gridding routine by Hargreaves and Beatty [36] is applied to each group with a 1024 × 1024 pixel grid size, a kernel width of 5, and an overgrid factor of 3 as described by Jackson *et al.* [34]. The resulting wavefields from the reflection groups are summed to arrive at a final output wavefield from each capillary and then combined to form a finely-sampled final output wavefield for the entire device due to one anode source point (50 nm pixel spacing). These data, when computed for each anode source point, are the coherent modes, *ψ _{n}*(

To validate the wavefield generation method used here, the output of a parallel-plane waveguide illuminated by a soft x-ray laser source has been computed for comparison to results by Kukhlevsky, *et al.* [8]. In the previous work, the intensity output of a SiO_{2} device (length = 21 cm, plane separation = 300 *μ*m) illuminated by a quasi-monochromatic source (*λ* = 20 nm, diameter = 50 *μ*m) [35] at a distance of 78.5 cm was measured at a distance of 33.4 cm from the output and simulated using an approximation method, namely the Fresnel-Kirchoff theory. The field at the output of the device is described as the superposition of *m* “field modes,” where each mode corresponds to the initial field having been reflected *m* times.

In Fig. 1, the simulated (labeled C) and experimental (labeled A) results reported in [8, see Fig. 3(a)] and the simulation results generated with the method described in this work (labeled B) are shown. One notes that both sets of simulation results predict the same peak and trough structure as well as the same approximate width for the output field. In the work by Kukhlevsky *et al.*, the simulated field is reported to be consistent with one in which *m _{max}* = 1. The field produced with the method in this paper also consists of

The transverse intensity response from a parallel-plane waveguide illuminated by a soft x-ray laser source. The simulation results using the techniques in this work (B) are compared to simulated results (C) and experimental measurements (A) from the literature **...**

The two sets of simulated results in Fig. 1 differ in terms of the predcted depth of the troughs. One notes, however, that neither the Fresnel-Kirchoff theory or geometrical optical approximation are exact solutions to Maxwell’s equations. The strong agreement between the two simulation results in terms of the peak-trough structure and the number of relevant reflections indicate that the method used in this work provides both qualitative and quantitative information about wavefields propagating through capillary devices.

To demonstrate the influence of the extended source on the device response, the ray locations at the device output are shown in Fig. 2 for different anode point positions. The ray output locations are strongly influenced by the anode point location, demonstrating that an experimentally significant anode size has been chosen. The ray locations are similar to those shown in previous simulation studies [27]. It is noteworthy that the ray locations alone do not provide the wavefield intensity and do not enable the evaluation of spatial coherence properties.

Simulation results showing the spatial distribution of ray locations at the polycapillary device output plane with varying input anode source point locations. The panels show (from left to right) the anode source point locations (−6.4 *μ* **...**

The intensity of the simulated output wavefield, computed in step 3, containing responses from all anode points is shown in Fig. 3. Additionally, the intensity of wavefields propagated away from the output plane, as computed in step 4 with a conventional Fresnel propagation operator, are shown. As the wavefields propagate away from the device, the intensity from the individual capillaries begins to diverge and overlap somewhat while the response from the overall device remains relatively collimated. This general behavior corresponds well with experimental results in the literature [37].

Finally, the spectral degree of coherence, computed in step 5 from the wavefield responses of individual anode source points, or modes, is shown in Fig. 4. The two-point spectral degree of coherence is shown in each frame as an image with reference to a point **r**_{1} in the output plane. The simulated optic has a highly structured coherence response that, in the output plane, is characterized by several “bright spots” corresponding to a high spatial degree of correlation with **r**_{1}. As the wavefields propagate away from the device, the structured spatial response begins to disperse, causing a more uniform effective spectral degree of coherence of approximately *μ*_{eff} = 0.3, especially when referenced to **r**_{1} = (0,0).

This framework for the analysis of polycapillary device spatial coherence properties opens the door for more comprehensive investigations of the influence of device parameters on the coherence response and optimization of these parameters for specific applications. One potential application for collimating x-ray optics generally is in phase-contrast imaging techniques. In-line phase contrast methods, specifically, suffer from a trade-off between the source coherence necessary for generation of phase-contrast features and the source flux. Collimating optics can potentially alter this relationship by capturing a larger solid angle of an x-ray tube’s output than would be possible otherwise. A key question, however, is how the wavefield coherence properties are changed due to propagation through the polycapillary device and the effect on image quality. The framework described here could potentially be used to evaluate the suitability of these devices for phase-contrast imaging. In addition, the simulation methods presented here may be useful for investigation of the wavefield output from polycapillary devices. Experimental measurements of these devices have demonstrated that the output intensity of certain devices contains ringing features that are strongly suggestive of the Fraunhofer diffraction patterns found due to light propagation through an aperture [17]. These effects are best characterized when the phase of propagating rays are accounted for in the output, especially as capillary devices are manufactured at smaller diameters where these patterns are most evident [18]. Interference effects, which are well known to occur in multiple-reflection optics [8,16], may also be investigated with access to simulations of the wavefield without the need for an analytic solution to the wave equation.

The framework presented in this work consists of an adaptable five-step process. In future work, these steps may be altered to address specific topics of interest. For example, enhanced polycapillary device simulation methods may be incorporated into step 2. Specifically, the effects of capillary waviness and surface roughness have not been included here and may potentially have a substantial effect on the coherence properties of a device. In addition, a number of methods exist to perform the wavefield construction in step 3, the accuracy of which could form the basis of future studies. The methods presented here form the basis for improved understanding of the properties of polycapillary devices and provide a framework for future investigations into the effects of polycapillary optical devices on spatial coherence properties.

We acknowledge the helpful suggestions from the anonymous reviewer that served to significantly improve this paper. We are also grateful to Drs. Carolyn A. MacDonald and P. Scott Carney for their helpful advice and guidance. PSC advised on the use of the CMD and geometrical optics. AMZ acknowledges support from the NIH ( CA136102) and MAA acknowledges support from the NSF ( CBET-0546113 and CBET-0854430) and NIH ( EB010049 and EB009715).

1. Schroer C., Lengeler B., “X-ray optics,” in Springer Handbook of Lasers and Optics, Träger F., editor. , ed., (Springer-Verlag, 2007).10.1007/978-0-387-30420-5_18 [Cross Ref]

2. MacDonald C. A., Gibson W. M., “Applications and advances in polycapillary optics,” X-ray Spectrom. 32, 258–268 (2003).10.1002/xrs.648 [Cross Ref]

3. Alexandrov Yu. M., Dabagov S. B., Kumakhov M. A., Murashova V. A., Fedin D. A., Fedorchuk R. V., Yakimenko M. N., “Peculiarities of photon transmission through capillary systems,” Nucl. Instrum. Methods Phys. Res., Sect. B 134, 174–180 (1998).10.1016/S0168-583X(98)00554-0 [Cross Ref]

4. Dabagov S. B., Marcelli A., “Single-reflection regime of x rays that travel into a monocapillary,” Appl. Opt. 38, 7494–7497 (1999).10.1364/AO.38.007494 [PubMed] [Cross Ref]

5. Kukhlevsky S. V., Flora F., Marinai A., Nyitray G., Kozma Zs., Ritucci A., Palladino L., Reale A., Tomassetti G., “Wave-optics treatment of x-rays passing through tapered capillary guides,” X-Ray Spectrom. 29, 354–359 (2000).10.1002/1097-4539(200009/10)29:5<354::AID-XRS436>3.0.CO;2-G [Cross Ref]

6. Dabagov S. B., “Wave theory of x-ray scattering in capillary structures,” X-Ray Spectrom. 32, 223–228 (2003).10.1002/xrs.593 [Cross Ref]

7. Vincze L., Janssens K., Adams F., Rindby A., “Detained ray-tracing code for capillary optics,” X-Ray Spectrom. 24, 27–37 (1995).10.1002/xrs.1300240107 [Cross Ref]

8. Kukhlevsky S. V., Flora F., Marinai A., Nyitray G., Ritucci A., Palladino L., Reale A., Tomassetti G., “Diffraction of X-ray beams in capillary waveguides,” Nucl. Instrum. Methods Phys. Res., Sect. B 168, 276–282 (2000).10.1016/S0168-583X(99)00874-5 [Cross Ref]

9. Kukhlevsky S. V., “Interference and diffraction in capillary x-ray optics,” X-Ray Spectrom. 32, 223–228 (2003).10.1002/xrs.644 [Cross Ref]

10. Xiao Q. F., Poturaev S. V., “Polycapillary-based X-ray optics,” Nucl. Instrum. Methods Phys. Res., Sect. A 347, 376–383 (1994).10.1016/0168-9002(94)91913-5 [Cross Ref]

11. Liu A., “The X-ray distribution after a focussing polycapillary a shadow simulation,” Nucl. Instrum. Methods Phys. Res., Sect. B 243, 223–226 (2006).10.1016/j.nimb.2005.07.209 [Cross Ref]

12. Welnak C., Chen G. J., Cerrina F., “Shadow: a synchrotron radiation and X-ray optics simulation tool,” Nucl. Instrum. Methods Phys. Res., Sect. A 347, 344–347 (1994).10.1016/0168-9002(94)91906-2 [Cross Ref]

13. Hampai D., Dabagov S. B., Cappuccio G., Cibin G., “X-ray propagation through hollow channel: PolyCAD - a ray tracing code,” Nucl. Instrum. Methods Phys. Res., Sect. B 244, 481–488 (2006).10.1016/j.nimb.2005.10.009 [Cross Ref]

14. Hampai D., Dabagov S. B., Cappuccio G., Cibin G., “X-ray propagation through polycapillary optics studied through a ray tracing approach,” Spectrochim. Acta, Part B 62, 608–614 (2007).10.1016/j.sab.2007.03.002 [Cross Ref]

15. Dabagov S. B., Kumakhov M. A., Nikitina S. V., Murashova V. A., Fedorchuk R. V., Yakimenko M. N., “Observation of interference effects at the focus of an x-ray lens,” J. Synchrotron Radiat. 2, 132–135 (1995).10.1107/S0909049595003864 [PubMed] [Cross Ref]

16. Dabagov S. B., Kumakhov M. A., Nikitina S. V., “On the interference of X-rays in multiple reflection optics,” Phys. Lett. A 203, 279–282 (1995).10.1016/0375-9601(95)00393-H [Cross Ref]

17. Bjeoumikhov A., “Observation of peculiarities in angular distributions of X-ray radiation after propagation through polycapillary structures,” Phys. Lett. A 360, 405–410 (2007).10.1016/j.physleta.2006.08.046 [Cross Ref]

18. Bjeoumikhov A., Bjeoumikhova S., Riesemeier H., Radtke M., Wedell R., “Propagation of synchrotron radiation through nanocapillary structures,” Phys. Lett. A 366, 283–288 (2007).10.1016/j.physleta.2007.01.092 [Cross Ref]

19. Dabagov S. B., Fedorchuk R. V., Murashova V. A., Nikitina S. V., Yakimenko M. N., “Interference phenomenon under focusing of synchrotron radiation by a Kumakhov lens,” Nucl. Instrum. Methods Phys. Res., Sect. B 108, 213–218 (1996).10.1016/0168-583X(95)00835-7 [Cross Ref]

20. Vincze L., Janssens K., Kukhlevsky S. V., “Simulation of polycapillary lenses for coherent and partially coherent x-rays,” Proc. SPIE 5536, 81–85 (2004).10.1117/12.560740 [Cross Ref]

21. Mandel L., Wolf E., Optical Coherence and Quantum Optics (Cambridge University Press, 1995).

22. Wolf E., “New spectral representation of random sources and of the partially coherent field that they generate,” Opt. Commun. 38, 3–6 (1981).10.1016/0030-4018(81)90295-9 [Cross Ref]

23. Liu H., Mu G., Lin L., “Propagation theories of partially coherent electromagnetic fields based on coherent or separated-coordinate mode decomposition,” J. Opt. Soc. Am. A 23, 2208–2218 (2006).10.1364/JOSAA.23.002208 [PubMed] [Cross Ref]

24. Zysk A. M., Carney P. S., Schotland J. C., “Eikonal method for calculation of coherence functions,” Phys. Rev. Lett. 95, 043904 (2005).10.1103/PhysRevLett.95.043904 [PubMed] [Cross Ref]

25. Schoonover R. W., Zysk A. M., Carney P. S., “Geometrical optics limit of stochastic electromagnetic fields,” Phys. Rev. A 77, 043831 (2008).10.1103/PhysRevA.77.043831 [Cross Ref]

26. Liu A., “Simulation of x-ray propagation in a straight capillary,” Math. Comput. Simulat. 65, 251–256 (2004).10.1016/j.matcom.2004.01.001 [Cross Ref]

27. Liu A., “Simulation of x-ray beam collimation by polycapillaries,” Nucl. Instrum. Methods Phys. Res., Sect. B 234, 555–562 (2005).10.1016/j.nimb.2005.02.001 [Cross Ref]

28. Xiao Q., Ponomarev I., Kolomitsev A., Kimball J., “Numerical simulations for capillary-based x-ray optics,” Proc. SPIE 1736, 227–238 (1992).10.1117/12.140477 [Cross Ref]

29. Popov M., “A new method of computation of wave fields using gaussian beams,” Wave Motion 4, 85–97 (1982).10.1016/0165-2125(82)90016-6 [Cross Ref]

30. Norris A., “Complex point-source representation of real point sources and the gaussian beam summation method,” J. Opt. Soc. Am. A 3, 2005–2010 (1986).10.1364/JOSAA.3.002005 [Cross Ref]

31. Forbes G., Alonso M., “Using rays better. I. theory for smoothly varying media,” J. Opt. Soc. Am. A 18, 1132–1145 (2001).10.1364/JOSAA.18.001132 [PubMed] [Cross Ref]

32. Heilpern T., Heyman E., Timchenko V., “A beam summation algorithm for wave radiation and guidance in stratified media,” J. Acoust. Soc. Am. 121, 1856–1864 (2007).10.1121/1.2537221 [PubMed] [Cross Ref]

33. Svensson E., “Gaussian beam summation in shallow waveguides,” Wave Motion 45, 445–456 (2008).10.1016/j.wavemoti.2007.09.003 [Cross Ref]

34. Jackson J., Meyer C., Nishimura D., Macovski A., “Selection of a convolution function for fourier inversion using gridding [computerised tomography application],” IEEE Trans. Med. Imag. 10, 473–478 (1991).10.1109/42.97598 [PubMed] [Cross Ref]

35. Bollanti S., Albertano P., Belli M., Di Lazzaro P., Faenov A. Ya., Flora F., Giordano G., Grilli A., Ianzinni F., Kukhlevsky S. V., Letardi T., Nottola A., Palladino L., Pikuz T., Reale A., Reale L., Scafati A., Tabocchini M. A., Turcu I. C. E., Vigli-Papadaki K., Schina G., “Soft X-ray plasma source for atmospheric-pressure microscopy, radiobiology and other applications,” Il Nuovo Cimento 20, 1685–1701 (1998).

36. Hargreaves B., Beatty P., “Gridding functions,” http://mrsrl.stanford.edu/brian/gridding/

37. Hampai D., Dabagov S. B., Cappuccio G., Cibin G., Sessa V., “X-ray micro-imaging by capillary optics,” Spectrochim. Acta, Part B 64, 1180–1184 (2009).10.1016/j.sab.2009.08.006 [Cross Ref]

Articles from Optics Express are provided here courtesy of **Optical Society of America**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |