research papers\(\def\hfill{\hskip 5em}\def\hfil{\hskip 3em}\def\eqno#1{\hfil {#1}}\)

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775

Synchrotron radiation transmission by two coupled flat microchannel plates: new opportunities to control the focal spot characteristics

crossmark logo

aPhysics Department, Southern Federal University, Sorge Str. 5, 344090 Rostov-on-Don, Russian Federation, bINFN – Laboratori Nazionali di Frascati, Via Enrico Fermi 54, 00044 Frascati, Italy, cRICMASS, Rome International Center for Materials Science Superstripes, 00185 Rome, Italy, dIstituto Struttura della Materia, CNR, Via del Fosso del Cavaliere 100, 00133 Rome, Italy, eRAS P. N. Lebedev Physical Institute, Leninsky Pr. 53, 119991 Moscow, Russian Federation, and fNR Nuclear University MEPhI, Kashirskoe Sh. 31, 115409 Moscow, Russian Federation
*Correspondence e-mail: mazurmik@gmail.com, augusto.marcelli@lnf.infn.it

Edited by S. Svensson, Uppsala University, Sweden (Received 2 October 2021; accepted 3 December 2021; online 19 January 2022)

An improved theoretical model to calculate the focal spot properties of coherent synchrotron radiation (SR) soft X-ray beams by combining and aligning two microchannel plates (MCPs) is presented. The diffraction patterns of the radiation behind the MCP system are simulated in the framework of the electrodynamical model of the radiation emission from two-dimensional finite antenna arrays. Simulations show that this particular optical device focuses the soft X-ray radiation in a circular central spot with a radius of ∼4 µm. The study points out that such MCP-based devices may achieve micrometre and sub-micrometre spot sizes as required by many applications in the soft X-ray range. Finally, based on experimental and theoretical results of the radiation transmission by this MCP-based device, a new method to characterize the spatial properties of brilliant SR sources is discussed.

1. Introduction

Third-generation synchrotron radiation facilities and free-electron lasers (FELs) are brilliant sources of coherent X-rays open to a interdisciplinary scientific community that allow many new experiments to be performed. In addition to undulator devices installed at third-generation synchrotron radiation (SR) facilities, FELs are unique radiation sources that deliver narrow-band coherent radiation beams. Operational SR facilities are characterized by wide spectral ranges delivering radiation from the IR/THz range up to the `nm domain' (Tsuji et al., 2004[Tsuji, K., Injuk, J. & Van Grieken, R. (2004). Editors. X-ray Spectrometry: Recent Technological Advances. John Wiley & Sons Ltd.]).

In addition to the brilliance, an important property of the X-ray beam is the coherence, which characterizes the correlation between two electromagnetic waves at different points in space in both transverse and longitudinal directions. As known, at high spatial coherence the beam divergence tends to 2λ/d, where λ and d are the wavelength and diameter of the radiation beam, respectively (Hiroyuki, 2002[Hiroyuki, D. (2002). Rep. Prog. Phys. 65, 1513-1576.]).

The technologies to obtain a micrometre-size soft X-ray spot are described by Tsuji et al. (2004[Tsuji, K., Injuk, J. & Van Grieken, R. (2004). Editors. X-ray Spectrometry: Recent Technological Advances. John Wiley & Sons Ltd.]). There, a comparison of various combinations of slits and collimators with polycapillary-based optical systems as grazing-incidence optics is discussed. The main characteristics of polycapillary spectroscopic devices are size and divergence of the transmitted beam. These quite compact diffractive optics are applied in X-ray spectroscopy and in many other analytical techniques obtaining small spots (from several micrometres to tens of nanometres) or a high photon density at the focal point or both, i.e. in the micrometre and sub-micrometre ranges (Pfeiffer et al., 2002[Pfeiffer, F., David, C., Burghammer, M., Riekel, C. & Salditt, T. (2002). Science, 297, 230-234.]; Dabagov, 2003a[Dabagov, S. B. (2003a). X-ray Spectrom. 32, 179-185.]; Bukreeva et al., 2006[Bukreeva, I., Popov, A., Pelliccia, D., Cedola, A., Dabagov, S. & Lagomarsino, S. (2006). Phys. Rev. Lett. 97, 184801.], 2010[Bukreeva, I., Pelliccia, D., Cedola, A., Scarinci, F., Ilie, M., Giannini, C., De Caro, L. & Lagomarsino, S. (2010). J. Synchrotron Rad. 17, 61-68.]; Sun et al., 2009[Sun, T., Zhang, M., Liu, Z., Zhang, Z., Li, G., Ma, Y., Du, X., Jia, Q., Chen, Y., Yuan, Q., Huang, W., Zhu, P. & Ding, X. (2009). J. Synchrotron Rad. 16, 116-118.]; MacDonald, 2010[MacDonald, C. A. (2010). X-ray Opt. Instrum. 10, 867049.]; Dabagov & Gladkikh, 2019[Dabagov, S. B. & Gladkikh, Yu. P. (2019). Radiat. Phys. Chem. 154, 3-16.]).

X-ray beams shaped by various polycapillary devices have been characterized since the 1990s (Kumakhov & Komarov, 1990[Kumakhov, M. A. & Komarov, F. F. (1990). Phys. Rep. 191, 289-350.]; Bilderback et al., 1994[Bilderback, D. H., Hoffman, S. A. & Thiel, D. J. (1994). Science, 263, 201-203.]; Dabagov et al., 1995a[Dabagov, S. B., Kumakhov, M. A., Nikitina, S. V., Murashova, V. A., Fedorchuk, R. V. & Yakimenko, M. N. (1995a). J. Synchrotron Rad. 2, 132-135.]; Dabagov, 2003b[Dabagov, S. B. (2003b). Phys.-Usp. 46, 1053-1075.]; MacDonald & Gibson, 2000[MacDonald, C. A. & Gibson, W. M. (2000). Handbook of Optics, Vol. III, ch. 30, pp. 30.1-30.12. McGraw-Hill.]). Among polycapillary-based optical systems, microchannel plates (MCPs) are different compact devices allowing beam focusing and shaping. Studies of MCPs date back to the 1990s (Chapman et al., 1990[Chapman, H. N., Nugent, K. A., Wilkins, S. W. & Davis, T. J. (1990). J. X-ray Sci. Tech. 2, 117-126.], 1993[Chapman, H. N., Nugent, K. A. & Wilkins, S. W. (1993). Appl. Opt. 32, 6316-6332.]; Kaaret et al., 1992[Kaaret, P., Geissbühler, P., Chen, A. & Glavinas, E. (1992). Appl. Opt. 31, 7339-7343.]; Nussey, 2005[Nussey, P. (2005). PhD Thesis, University Leicester, UK.]). The radiation propagation through MCPs is typically described within classical (ray) optics as multiple reflections from their channel's internal walls. This approach, however, neglects the wave nature of the radiation and, in particular, the phase evolution of the propagating wave. This approach also cannot explain interference phenomena occurring while crossing different components (Dabagov et al., 1995b[Dabagov, S. B., Kumakhov, M. A. & Nikitina, S. V. (1995b). Phys. Lett. A, 203, 279-282.], 2000[Dabagov, S. B., Marcelli, A., Murashova, V. A., Svyatoslavsky, N. L., Fedorchuk, R. V. & Yakimenko, M. N. (2000). Appl. Opt. 39, 3338-3343.]). Indeed, surface and bulk channeling of radiation in micro- and nano-guides in capillary-based systems have been described and observed (Dabagov & Uberall, 2007[Dabagov, S. B. & Uberall, H. (2007). Nucl. Instrum. Methods Phys. Res. A, 580, 756-763.], 2008[Dabagov, S. B. & Uberall, H. (2008). Nucl. Instrum. Methods Phys. Res. B, 266, 3881-3887.]). In this work we analyze the radiation propagation in thin MCPs with narrow channels taking into account the phase evolution of the transmitted waves. The approach allows interference and focusing of X-ray beams behind different multichannel MCP configurations.

2. MCP optical systems

We have characterized several different MCPs, made by Si–Pb glass with composition (PbO)0.7(SiO2)0.3 manufactured by the Vladikavkaz Technological Center `Baspik' (https://www.baspik.com/eng/products/nauka/) (Mazuritskiy et al., 2019b[Mazuritskiy, M. I., Lerer, A. M., Kulov, S. K. & Samkanashvili, D. G. (2019b). J. Synch. Investig. 13, 499-507.]). We used relatively thin MCPs, e.g. 0.3–1.5 mm thick, with 104–107 hollow hexagonal-structured microchannels of equal length, extending through the whole thickness. The length-to-diameter ratio of a typical hollow microchannel spans from 40 to 150. Exhibiting rather high transmission efficiency (up to ∼60%) they are efficient waveguide optical systems. Moreover, MCP optical properties allow increasing the radiation density as well as shaping the beam (Mazuritskiy et al., 2016a[Mazuritskiy, M. I., Dabagov, S. B., Marcelli, A., Lerer, A. M. & Dziedzic Kocurek, K. (2016a). J. Synchrotron Rad. 23, 274-280.],b[Mazuritskiy, M. I., Lerer, A. M. & Makhno, P. V. (2016b). J. Exp. Theor. Phys. 123, 942-949.], 2017a[Mazuritskiy, M. I. & Lerer, A. M. (2017a). JETP Lett. 105, 572-576.],b[Mazuritskiy, M. I., Dabagov, S. B., Lerer, A. M., Dziedzic Kocurek, K., Sokolov, A., Coreno, M., Turchini, S., D'Elia, A., Sacchi, M. & Marcelli, A. (2017b). Nucl. Instrum. Methods Phys. Res. B, 402, 282-286.], 2018[Mazuritskiy, M. I., Lerer, A. M., Marcelli, A., Dabagov, S. B., Coreno, M. & D'Elia, A. (2018). JETP Lett. 107, 600-605.], 2019a[Mazuritskiy, M. I., Lerer, A. M., Dabagov, S. B., Marcelli, A., Hampai, D. & Dziedzic Kocurek, K. (2019a). J. Synch. Investig. 13, 1005-1013.]). Although many experimental and theoretical investigations have been published on single flat or bent MCPs, no data have been published on assembled and aligned MCPs.

In this work we discuss an original layout named Microchannel Aligned Assembled System (MAAS), which is formed by two different flat MCPs assembled such that their axis of the channels are maintained parallel. Fig. 1[link] shows the layout with the parallel coherent beam (1) propagating along the z-axis and two coaxial plates at the distance d1. In this layout, the primary beam – a monochromatic coherent plane wave with a Gaussian amplitude distribution in the transverse cross-section – is focused by the hollow microchannels of both MAAS MCPs. The distance d1 between two MCPs and the detector position at the distance d2 from the second MCP are the main parameters of the simulation. The d2 distance also determines the 2D-detector offset position (5).

[Figure 1]
Figure 1
Scheme of two flat parallel coaxial MCPs for focusing of the X-ray beam; d1 is the distance between the two MCPs and d2 the distance between the detector and second MCP. 1 – primary radiation; 2 – first MCP; 3 (inset) – image of the capillary structure of both MCPs; 4 – second MCP; 5 – 2D detector of radiation; 6 (insert) – focal spot profile; 7 (insert) – focal spot cross-section.

3. The theoretical model

In the simulations of the MAAS device (Fig. 1[link]), microchannel diameters as well as pitch sizes and periods may change to set the optimal parameters for focusing, interference, diffraction, etc. The model considers two MCPs of ideal hexagonal symmetry in the transverse cross-sections and with the symmetry axis along the z-axis. In the simulations, each MCP is characterized by its dispersion equation, optimized by suitable microchannels. The first MCP is used to shape the primary parallel beam (coherent or partially coherent). Behind the first MCP the beam becomes conical in shape, and this `secondary beam' defines the radiation footprint on the surface of the second MCP, chosen thicker and with smaller microchannels.

Focusing depends on the energy, the refractive index of the MCP medium, the channel diameter, the number of illuminated microchannels, and the thickness. By changing the distance between the two MCPs, the parameters of the antenna array of the second MCP change the diffraction pattern. The radiation reaches the detector plane forming a pattern with an intense narrow center peak (`zero position') and weaker side peaks. Let us assume that the original radiation is represented by a monochromatic plane wave with a Gaussian intensity distribution in the transverse cross-section,

[U_0(x,y,z) = \exp \left[ -\left( {{x-x_0}\over{\Delta_x}}\right)^2 -\left( {{x-y_0}\over{\Delta_y}}\right)^2 \right] \exp(-ikz), \eqno(1)]

where x0, y0 are the coordinates of the beam center; Δx, Δy are the parameters of the Gaussian function and k is the wavevector. We consider now the radiation propagating along the z-axis (Fig. 1[link]) normal to the surfaces of two flat MCPs of the MAAS device. The diffraction field formed behind the first MCP is the radiation source (`a secondary source') of the second MCP. Moreover,

(i) Both devices have a hexagonal packing of hollow cylindrical microchannels characterized by microchannel radii and pitch sizes: R1 and D1 and R2 and D2.

(ii) The coaxial MCPs are set at the distance d1, while d2 (Fig. 1[link]) is the distance between the second MCP and the detector.

(iii) Considering the hexagonal symmetry, the coordinates of the microchannel centers for both MCPs can be calculated using the following expression,

[x_{mn}^{\left(\,j\,\right)} = m{D_j} + n{D_j}\sin \psi, \qquad y_{mn}^{\left(\,j\,\right)} = n{D_j}\cos \psi, \qquad \psi = \pi /6.]

The microchannels are defined by `m,n' in the xy-plate and j = 1 or 2 for the first and second MCP, respectively.

We evaluate the diffraction at the far zone in the frame of the Kirchhoff–Huygens approximation. The model rules out any interaction among MCP microchannels, and the wavefield emitted by each cylindrical channel of the jth MCP can be described as

[\eqalign{ &V_{mn}^{\,\left(\,j\,\right)}\left({x,y,z} \right) = {{ \exp \left({ - ikr_{mn}^{\left(\,j\,\right)}} \right) }\over{ r_{mn}^{\left(\,j\,\right)} }} \, \Phi _{mn}^{\left(\,j\,\right)}\left({x,y,z} \right), \cr& r_{mn}^{\,\left(\,j\,\right)} = \left[ {{{\left({x - x_{mn}^{\left(\,j\,\right)}} \right)}^2} + {{\left({y - y_{mn}^{\left(\,j\,\right)}} \right)}^2} + {{\left({z - {z_j}} \right)}^2}} \right]^{1/2},} \eqno(2)]

defined by the space distribution [\Phi_{mn}^{\left(\,j\,\right)}\left({x,y,z} \right)] of the radiation pattern at the exit of each channel (m,n) of the jth MCP. The number of illuminated channels of the first MCP equals (2M1 + 1) × (2N1 + 1), while of the second equals (2M2 + 1) × (2N2 + 1). M1 and N1 are determined by the profile of the primary beam and the footprint on the front surface of the first MCP. At variance, due to the conical shape of the beam emerging from the first MCP (see Fig. 1[link]), the values M2, N2 depend also on the distance d1. Within the MAAS framework, at a large d1 distance we have many illuminated channels on the front side of the second MCP. The footprint dimension, obviously, affects the interference at the second MCP exit. The electromagnetic field at this plane is due to the interference of the beams emerging by the channels of the first MCP,

[\eqalignno{ U\left(x,y,d_1\right) \,\simeq\, {}& {\bar V}^{\,\left(1\right)} \left(x,y,d_1\right) \sum\limits_{m\,=\,-M_1}^{M_1} \sum\limits_{n\,=\,-N_1}^{N_1} U_0\left(x_{mn}^{\left(1\right)},y_{mn}^{\left(1\right)}\right) \cr& \times \exp \left[i\left(\chi_x^{\left(1\right)}x_{mn}^{\left(1\right)} + \chi_y^{\left(1\right)}y_{mn}^{\left(1\right)} \right) \right], &(3)}]

[\chi_x^{\left(\,j\,\right)} = k{{x}\over{d_j}}, \qquad \chi_y^{\left(\,j\,\right)} = k{{y}\over{d_j}}. \eqno(4)]

In our model we assume the radiation field constant inside each microchannel of the second MCP (multiple reflection mode),

[U_1(x,y,d_1) \simeq U_1\left(x_{mn}^{(2)},y_{mn}^{(2)},d_1\right), \eqno(5)]

so that the final diffraction pattern at the detector plane can be calculated as follows,

[\eqalignno{ U_2\left(x,y,d_2\right) \, \simeq \, {}& {\bar V}^{\,\left(2\right)} \left(x,y,d_2\right) \sum\limits_{m\,=\,-M_2}^{M_2} \sum\limits_{n\,=\,-N_2}^{N_2} U_1\left(x_{mn}^{\left(2\right)},y_{mn}^{\left(2\right)},d_1\right) \cr& \times \exp\left[i\left(\chi_x^{(2)}x_{mn}^{\left(2\right)} + \chi_y^{(2)}y_{mn}^{\left(2\right)} \right) \right]. &(6)}]

While all equations have been written for parallel flat MCPs with the central axes (m = n = 0) set along the z-axes, the model enables simulations for any ([\Delta{x},\Delta{y}]) displacement of the first MCP as well as for any angular rotation respect to the initial position of ([{\bar x}_{mn}^{\,(1)},{\bar y}_{mn}^{\,(1)}]),

[\eqalignno{ x_{mn}^{\left(1\right)} & = \Delta x + {\bar x}_{mn}^{\,\left(1\right)}C - {\bar y}_{mn}^{\,\left(1\right)}S, \qquad y_{mn}^{\left(1\right)} = \Delta y + {\bar x}_{mn}^{\,\left(1\right)}S + {\bar y}_{mn}^{\left(1\right)}C, \cr C & = \cos\varphi_M, \qquad S = \sin\varphi_M.&(7)}]

For the sake of simplicity the field at the exit of a single microchannel can be derived for a circular hole within an opaque screen. This model of propagation inside a circular hole well approximates a waveguide and the radiation field at the exit can be calculated using the following equation (Mazuritskiy & Lerer, 2015[Mazuritskiy, M. I. & Lerer, A. M. (2015). JETP Lett. 102, 483-486.], 2016b[Mazuritskiy, M. I., Lerer, A. M. & Makhno, P. V. (2016b). J. Exp. Theor. Phys. 123, 942-949.]),

[\eqalign{ {\bar V}_{\!j\,}\left(x,y,z\right) & = {{\exp\left(-ikr\right)}\over{r}} \,{\hat R}_j^2\, {{J_1\left(\rho{\hat R}_j\right)}\over{\rho}}, \cr r_{mn}^{\left(\,j\,\right)} & = \left[x^2+y^2+\left(z-z_j\right)^2\right]^{1/2}, \cr \rho & = k\left(x^2+y^2\right)^{1/2} \big/ r,} \eqno(8)]

where [J_1(\rho{\hat R}_j)] is the Bessel function and [{\hat R_j}], the model parameter, is the effective radius of the microchannel hole. For a single MCP the value [{\hat R_j}] can be evaluated by comparing experimental and simulated diffraction patterns. The decrease of the MCP thickness, i.e. the channel length, can be described by a smaller radius [{\hat R_j}] because this parameter defines the solid angle of the radiation spatial distribution at the microcapillary exit. Let us now examine the pattern of the radiation emerged by a single microchannel that can be applied also to any selected channel or set of channels within these assumptions:

(i) The validity of the complex propagation constants for the waveguide modes inside a circular hollow dielectric guide inhomogeneous in depth (see Mazuritskiy et al., 2014[Mazuritskiy, M. I., Dabagov, S. B., Marcelli, A., Lerer, A., Novakovich, A. & Dziedzic Kocurek, K. (2014). J. Opt. Soc. Am. B, 31, 2182-2187.]).

(ii) The validity of the amplitudes of the waveguide modes within the Kirchhoff–Huygens approach.

(iii) The validity of the mode amplitudes at the waveguide output using the complex propagation constants.

(iv) The far-field angular distributions of the radiation from each waveguide using the Kirchhoff–Huygens method.

The electric field of the plane wave inside the waveguide (capillary) is

[E^{\,{\rm{ext}}}(x,y,z) = \exp\left[i\left(k_x^{\prime}x + k_y^{\prime}y + k_z^{\prime}z\right)\right], \eqno(9)]

where [k_{x,y,z}^{\prime}] are the wavevector projections. The primary radiation field at the aperture of the waveguide and inside the waveguide can be written in a similar way,

[\eqalignno{ E^{\,{\rm{ext}}}(x,y,0) & = \exp\left[i\left(k_x^{\prime}x + k_y^{\prime}y\right)\right] \cr & = \sum\limits_{m\,=\,0}^\infty v_m\cos (m\varphi) \sum\limits_{n\,=\,1}^\infty U_{mn}\,R_n^{\,(m)}(r), & (10)}]

where

[v_m = \left\{ \matrix{ 1/2,\hfill & m=0,\hfill \cr 1,\hfill & m\ne0,\hfill}\right.]

and Umn Rn (m)(r) is the amplitude that depends on the radius of the wavefield distribution characterized by the indexes m,n. To calculate equation (10)[link] we use the M and N maxima corresponding to the m and n indexes. M is usually not greater than 3, while N is the maximum number of the wave modes accepted by the capillary. Taking into account the orthogonal properties of the function we can rewrite equation (10)[link] as

[U_{mn} = {{1}\over{\pi N_{mn}}} \int\limits_0^{2\pi} \cos (m\varphi)\,{\rm{d}}\varphi \int\limits_0^\infty r\, \exp\left[i\left(k_x^{\prime}x + k_y^{\prime}y\right)\right] \, R_n^{\,(m)}(r)\,{\rm{d}}r, \eqno(11)]

[N_{mn} = \int\limits_0^\infty r \left[ R_n^{\,(m)}(r)\right]^2\,{\rm{d}}r, \eqno(12)]

where [k_x^{\prime}] = [\kappa_0\cos\vartheta], [k_y^{\prime}] = [\kappa_0\sin\vartheta], x = [r\cos\phi], y = [r\sin\phi] and κ0 is the perpendicular component of the radiation wavevector to the capillary axis at the entrance. In addition, we assume the wavefield to be confined inside the microcapillary. In this approximation, integrating equations (10)[link]–(12)[link] over the segment and assuming for the transverse wavenumber, we obtain the following expressions,

[\eqalignno{ U_{mn} & = {{1}\over{\pi N_{mn}}} \int\limits_0^{2\pi} \cos(m \varphi)\,{\rm{d}}\varphi \int\limits_0^\rho r \exp \big[i\kappa_0\,r\cos(\varphi-\vartheta)\big] \cr& \quad \times J_m\left(\kappa_n^{(m)}r\right)\,{\rm{d}}r_{\vphantom{\big|}} \cr& = {{2i^{\,m}\cos(m\vartheta)}\over{N_{mn}}} \, F\left(\kappa_0,\kappa_n^{(m)},\rho\right), &(13)}]

[\eqalignno{ N_{mn} &\simeq \int\limits_0^\rho r\, \left[J_m\left(\kappa_n^{(m)}r\right)\right]^2{\rm{d}}r \cr& =\, {{\rho^2}\over{2}}\Big\{\left[J_m\left(\kappa_n^{(m)}\rho\right)\right]^2 \,-\,\, J_{m-1}\left(\kappa_n^{(m)}\rho\right) \,J_{m+1}\left(\kappa_n^{(m)}\rho\right)\Big\}, &(14)}]

[\eqalignno{ F(\alpha,\beta,\rho) & = \int\limits_0^\rho r \,J_m(\alpha r)\, J_m(\beta r)\,{\rm{d}}r \cr& = {{\rho}\over{\alpha^2-\beta^2}} \,\Big[ \kappa_n^{(m)} J_m(\alpha\rho)\, J_{m-1}(\beta\rho) \cr& \quad - \kappa_0\,J_{m-1}(\alpha\rho)\,J_m(\beta\rho)\Big].&(15)}]

At the capillary exit the field is calculated as the sum of the waves,

[E(x,y,h) = \sum\limits_{m\,=\,0}^M v_m\cos(m\varphi) \sum\limits_{n\,=\,1}^N U_{mn}\,R_n^{(m)}(r) \exp\left[-i\gamma_n^{(m)}h\right], \eqno(16)]

where [\gamma_n^{(m)}] is the longitudinal wavenumber and h is the capillary length. The radiation field is defined as

[\eqalignno{ E^{\,{\rm{rad}}}(x,y,z) & = {{1}\over{2\pi}} \int\limits_{-\infty}^{\infty}\int\limits_{-\infty}^{\infty} E(x^{\prime},y^{\prime},h) \, {{\exp(-ik\rho)}\over{\rho}} \, {\rm{d}}x^{\prime}\,{\rm{d}}y^{\prime}, \cr \rho & = \left[(x-x^{\prime\,})^2+(y-y^{\prime})^2+(z-h)^2\,\right]^{1/2},&(17)}]

assuming

[\eqalign{ & \rho \simeq r-\tau, \qquad r=\left[x^2+y^2+(z-h)^2\right]^{1/2}, \cr& \tau={{xx^{\prime}+yy^{\prime}}\over{r}} = \sin\theta\,(\cos\varphi\,x^{\prime}+\sin\varphi\,y^{\prime}).} \eqno(18)]

Using the spherical coordinates r, θ, φ, the radiation field can be rewritten as

[\eqalignno{ E^{\,{\rm{rad}}}(x,y,z) = {}& {{\exp(-ikr)}\over{2\pi{r}}}\, \Phi(x,y,z) \cr& \times \int\limits_{-\infty}^\infty\int\limits_{-\infty}^\infty E(x^{\prime},y^{\prime},h) \exp(ik\tau)\,{\rm{d}}x^{\prime}\,{\rm{d}}y^{\prime}, &(19)}]

[\Phi(x,y,z) = \int\limits_{-\infty}^\infty\int\limits_{-\infty}^\infty E(x^{\prime},y^{\prime},h) \exp(ik\tau)\,{\rm{d}}x^{\prime}\,{\rm{d}}y^{\prime}, \eqno(20)]

and finally, replacing the field in equations (19)[link]–(20)[link] with that of equation (16)[link], we obtain

[\Phi(\theta,\varphi) = \sum\limits_{m\,=\,0}^Mv_m\sum\limits_{n\,=\,1}^NU_{mn} \exp\left[-i\gamma_n^{(m)}h\right]\, V_n^{\,(m)}(\theta,\varphi), \eqno(21)]

where

[\eqalignno{ V_n^{\,(m)}(\theta,\varphi) & = \int\limits_{-\infty}^\infty\int\limits_{-\infty}^\infty \cos(m\varphi)\,R_n^{(m)}(r) \exp(ik\tau) \,{\rm{d}}x^{\prime}\,{\rm{d}}y^{\prime} \cr & \simeq \int\limits_0^{2\pi}\cos(m\varphi^{\prime})\,{\rm{d}}\varphi^{\prime} \int\limits_0^{R_k} r^{\,\prime}\,J_m\left(\kappa_n^{(m)}r^{\,\prime}\right) \cr& \quad \times \exp\left[ik\sin\theta\cos(\varphi^{\prime}-\varphi)\right]\,{\rm{d}}r, &(22)}]

and this expression is equal to

[V_n^{\,(m)}(\theta,\varphi) = 2\pi i^{\,m} \cos(m\varphi)\,\,F\left(k\sin\theta,\kappa_n^{(m)},\rho\right). \eqno(23)]

4. Results and discussion

Using the above model, simulations were performed for a primary parallel beam with a cross section of 100 µm × 100 µm at an energy of 94 eV. The parameters of the calculations correspond to the energy and the characteristics associated with the BESSY II (Berlin) and ELETTRA (Trieste) experimental parameters.

In the experiments, the primary beam hits the first MCP with hexagonal packing, cylindrical channels of 10 µm diameter and 12 µm pitch, and with a ratio thickness to microchannel diameter of 40. The second MCP was thicker (0.5 mm) with a smaller microchannel diameter (3.4 µm), 4.0 µm pitch, and ratio thickness to microchannel diameter of 150. The distance d1 between MCPs varied between 100 and 300 mm, while d2 was in the range 10–30 mm (see Fig. 1[link]).

For these simulations the 2D-detector was placed perpendicularly to the direction of the primary beam (CiPo-Beamline; Elettra Synchrotron https://www.elettra.trieste.it/elettra-beamlines/cipo.html) and characterized by 10 µm spatial resolution (Marcelli et al., 2018[Marcelli, A., Mazuritskiy, M. I., Dabagov, S. B., Hampai, D., Lerer, A. M., Izotova, E. A., D'Elia, A., Turchini, S., Zema, N., Zuccaro, F., de Simone, M., Rezvani, S. J. & Coreno, M. (2018). J. Instrum. 13, C03035.]). Fig. 2[link] shows the simulated diffraction patterns for this layout. Panel (a) in Fig. 2[link] shows the radiation spatial distribution at the second MCP exit. The simulations of the radiation transmitted by the MAAS device show a minimum circular focal spot of ∼4 µm.

[Figure 2]
Figure 2
Spatial SR intensity distributions behind the device: (a) at the exit of the second MCP (just beyond the second plate); (b) at the entrance of the second MCP with 3.4 µm channel diameter and 4.0 µm pitch.

The diffraction pattern at the far-zone in Fig. 2[link](a) (at the exit of the second MCP) shows a highly intense and small divergent peak demonstrating that the SR beam transmitted by the two MCPs is concentrated at the center of the transverse space. In an ideal layout with two hexagonal MCPs perfectly aligned along the z-axis, this artificial photonic device works as a condenser of the primary parallel beam, but at the same time as an efficient diffractive X-ray optical device forming a constructive interference of the transmitted radiation.

Regarding its focusing efficiency, at 94 eV the field at the first MCP exit [Fig. 2[link](b)] exhibits a wide conical distribution. The secondary radiation source generated by this MCP changes propagating towards the second becoming the primary source at the entrance plane of the second MCP. The radiation change between MCPs redistributes the intensity and reshapes the beam. Actually, the primary parallel SR beam is different from the beam hitting the second MCP. The difference becomes more evident for a thicker MCP, since longer microchannels reduce the angular and spatial distributions of the radiation at the exit of the second MCP.

Simulations for a circular-shaped primary beam of 100–300 µm radius in cross section have revealed an optimal interpolated distance d1, at which the illuminated area of the second MCP is larger. For such a layout we maximize the number of microchannels transmitting radiation. Obviously, the MAAS device is more efficient in generating a brilliant radiation source either in a condenser or in a diffractive regime. Within the wave approach, this device exhibits a strong correlation of the diffraction pattern at the second MCP exit with the illuminated area of the first MCP. In Fig. 3[link](a), looking at the radiation intensity after the second MCP exit as a function of the primary radiation, the dashed line is the result of the calculations, while the red line is the fit.

[Figure 3]
Figure 3
(a) The radiation density at the exit of the second MCP as a function of the primary beam cross-section (squared radius R2). Inset 1 – for 100 µm; inset 2 – for 300 µm. (b) Comparison of beam profiles for different MCPs distances in the range 100–300 mm.

The pattern central maximum can be fit by a 2D-Gaussian function, i.e. insets 1 and 2 for X = 1 (R = 100 µm) and X = 9 (R = 300 µm), respectively. The radiation densities on the Y-axis are calculated as the ratio of the radiation intensity to the Gaussian profile area. As shown in Fig. 3[link](a), the fit, which represents the radiation density behind the optical system versus the profile area of a primary beam, is well described by a second-order polynomial curve.

The pattern is the result of the interference of waves emitted by each pair of symmetrical microchannels of the first MCP having the same phases at the center of the second MCP. Thus, the first plate shapes the primary beam, i.e. the radiation hitting the surface of the second MCP. For a defined layout and symmetry properties the observed interference pattern also depends by the spatial coherence of the incoming beam.

Devices working as diffractive optical elements allow controlling the radiation density, redistributing the radiation in a different way compared with the condensing regime. Actually, the presence of two MCPs makes possible the redistribution of the condensed radiation maximum at the focal position because of the interference of the radiation. A more concentrated and narrower central maximum surrounded by high-order contributions occurs, the latter reflecting the symmetry of the device, that in our case is hexagonal. Considering a coherent primary beam, this MCP system should be optimized to obtain an interference pattern at the focal plane with a well defined intense central peak.

Fig. 3[link](b) shows the profiles of the central peak for different distances between the MCPs. We changed the illuminated area at the entrance surface of the second MCP changing the distance d1 (see Fig. 1[link]). At large distances, i.e. at high d1 value, the illuminated area on the second MCP contains more microchannels [Fig. 3[link](b)]. The image shows the profiles of the intensity distribution of the radiation at the exit of the second MCP for different distances d1 in the range 100–300 mm. The flux increases of about one order of magnitude, while the FWHM of the Gaussian distributions remains constant as expected for a diffractive optics [Fig. 3[link](b)].

In our case, the interference of waves has a maximum at large distance, i.e. the radiation intensity at the center exhibits also a density gain. In these simulations we assumed a fully coherent beam (i.e. 100% degree of coherence) becoming partially coherent while propagating along the channels of the MCPs. With our model we could simulate various beam conditions including the coherence fraction.

The model allows also to simulate patterns at resonance detuning, considering small shifts of the MCP in the transversal plane (xy) and small angular rotations of the second MCP around the z-axis of the system. At any non-symmetrical MCPs-position the model shows a reduced contrast of the diffraction pattern as well as a broadening of the main peak. The latter exhibits the maximum density of the radiation on the z-axis only for a full axial symmetry of the field because of the constructive interference of waves in the far-zone after the second MCP.

Simulations for ideal layouts showed a clear reduction of the intensity (roughly a factor of 2–3) along the z-axis for a linear translation of MCPs of the order of tenths of a micrometre in the transversal plane (xy). We may hypothesize a decrease of the radiation density of the same order, which should take place for an (ideal) angular rotation of 2–3 mrad of the MCPs out of the symmetry position around the z-axis. However, in a `real' device a full decoupling of rotations and translations is never straightforward, and any `action' may result in much larger effects in terms of contrast and broadening of the experimental features. Moreover, in real MAAS devices the distance between the two MCPs cannot be changed, and a precise alignment depends also on this parameter. At smaller distances between the two MCPs, due to their hexagonal symmetry, the device could be sensitive even to `small' rotations within the 60° rotation range.

In the energy range investigated, the transmission efficiency of this focusing device is about 32%. This evaluation is based on experimental data collected using a single flat MCP at the same energy that exhibits an efficiency of 57% (Mazuritskiy et al., 2019a[Mazuritskiy, M. I., Lerer, A. M., Dabagov, S. B., Marcelli, A., Hampai, D. & Dziedzic Kocurek, K. (2019a). J. Synch. Investig. 13, 1005-1013.]). The results of this model point out that a device based on a pair of MCPs may be highly efficient in the soft X-ray range because of the low intensity loss caused by radiation absorption of the microchannel walls when precise alignment conditions for both MPCs are fulfilled.

To test the theoretical model and to compare these results with previous experiments performed using a single MCP, we also carried out simulations of X-ray diffraction at 94 eV. Simulations of an MCP, 0.3 mm thick with microchannels of diameter 3.4 µm and 4.2 µm pitch, are compared with experimental data in Fig. 4[link]. These data at high spatial resolution were collected at the Reflectometer UHV station of the SR facility of BESSY II (Sokolov et al., 2014[Sokolov, A. A., Eggenstein, F., Erko, A., Follath, R., Künstner, S., Mast, M., Schmidt, J. S., Senf, F., Siewert, F., Zeschke, Th. & Schäfers, F. (2014). Proc. SPIE, 9206, 92060J.], 2018[Sokolov, A., Sertsu, M. G., Gaupp, A., Lüttecke, M. & Schäfers, F. (2018). J. Synchrotron Rad. 25, 100-107.]). In this experiment the detector was set in the perpendicular plane at 310 mm distance from the MCP, whose front surface was set perpendicular to the primary SR beam. In this layout, the pattern of the diffracted radiation is characterized by an intense central peak, a set of first-order contributions and weak contributions of the second diffraction order. The comparison confirms the reliability of the model to describe the behavior of MAAS devices.

[Figure 4]
Figure 4
Single MCP diffraction pattern and beam profiles in the selected plane at 94 eV. (a, b) Experimental data; (c, d) simulated distributions.

In our device the field at the exit of the first MCP has a hexagonal symmetry. Due to this geometry for defined parameters and dimensions of the MCPs, estimations provide the validity of the plane wave approximation behind the first MCP that propagates toward the second MCP. Considering a distance of several centimetres (e.g. from 10 to 30 cm) the divergent radiation exiting from one channel of the first MCP interferes with the radiation emerging from the adjacent channels. If we consider an entrance cross-section of 100 µm × 100 µm for which many channels are illuminated, the structured radiation field created by the first MCP illuminates the front surface of the second MCP. The structured pattern at the exit is quasi-homogeneous at the entrance of the second MCP, i.e. each channel produces a beam that at a large distance (e.g. from 10 to 30 cm) illuminates hundreds of channels of the second MCP.

At variance, for a device with much smaller distances between plates (e.g. a few mm), the structured illumination produced by the first MCP may generate different characteristics patterns.

5. Conclusion

This work presents the behavior of a new optical layout named Microchannel Aligned Assembled System (MAAS) based on a couple of two flat MCPs simulating the transmission of synchrotron radiation inside this device. We introduced a theoretical model built from a two-dimensional finite antenna array consisting of non-interacting emitters, which in our case are the MCP cylindrical microchannels acting as waveguides. The diffraction patterns of the primary synchrotron radiation transmitted through this device were calculated. In the soft X-ray range this waveguide model reproduces the focusing of the radiation from this assembled couple of MCPs. The device we characterized combines two different flat MCPs: the first with microchannels with a larger diameter that guarantees a high transmission while shaping the beam; and the second with narrower and longer microchannels that focuses the radiation at the detector plane. In this way it may behave as an optical condenser or a diffractive optics.

The analysis of the radiation distribution at its exit highlights the possibility to control macroscopic properties of the parallel or quasi-parallel primary beam. Moreover, using the theoretical model we introduced, the analysis of patterns generated by properly aligned couples of MCPs could be used to evaluate the coherent fraction of a SR source or a FEL.

Acknowledgements

We are grateful to S. K. Kulov and D. G. Samkanashvili (OOO Vladikavkaz Technology Center BASPIK, Vladikavkaz, Russia) for manufacturing the device used in this research. We are deeply grateful to A. A. Sokolov for the assistance during the experimental runs at BESSY. Special thanks are due to D. Catone, S. Rinaldi, S. Turchini, N. Zema and F. Zuccaro of CNR-ISM for technical support during the experiments performed at Elettra. Acknowledgments are also due to M. Coreno, A. D'Elia, M. De Simone and S. J. Rezvani for many useful discussions. AM thanks Italian MUR, which supports activities related to Progetto Roadmap ESFRI/DFM.AD006.003 EUROFEL-ISM. AM acknowledges also the Italian Presidenza del Consiglio dei Ministri for its support through the DARA Department, within the framework of the MIAMI project managed in cooperation with the University of Milano Bicocca and the University of Roma Tre.

Funding information

The following funding is acknowledged: Helmholtz Zentrum Berlin, BESSY II (grant No. Project n. 181-06366ST; grant No. Project n. 191-07912ST); Elettra-Sincrotrone Trieste (grant No. Grant No. 20180083; grant No. Grant No. 20180241); CALIPSOplus, EU Framework Program for Research (grant No. Grant No. 730872).

References

First citationBilderback, D. H., Hoffman, S. A. & Thiel, D. J. (1994). Science, 263, 201–203.  CrossRef CAS PubMed Web of Science Google Scholar
First citationBukreeva, I., Pelliccia, D., Cedola, A., Scarinci, F., Ilie, M., Giannini, C., De Caro, L. & Lagomarsino, S. (2010). J. Synchrotron Rad. 17, 61–68.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationBukreeva, I., Popov, A., Pelliccia, D., Cedola, A., Dabagov, S. & Lagomarsino, S. (2006). Phys. Rev. Lett. 97, 184801.  Web of Science CrossRef PubMed Google Scholar
First citationChapman, H. N., Nugent, K. A. & Wilkins, S. W. (1993). Appl. Opt. 32, 6316–6332.  CrossRef CAS PubMed Web of Science Google Scholar
First citationChapman, H. N., Nugent, K. A., Wilkins, S. W. & Davis, T. J. (1990). J. X-ray Sci. Tech. 2, 117–126.  CAS Google Scholar
First citationDabagov, S. B. (2003a). X-ray Spectrom. 32, 179–185.  Web of Science CrossRef CAS Google Scholar
First citationDabagov, S. B. (2003b). Phys.-Usp. 46, 1053–1075.  CAS Google Scholar
First citationDabagov, S. B. & Gladkikh, Yu. P. (2019). Radiat. Phys. Chem. 154, 3–16.  Web of Science CrossRef CAS Google Scholar
First citationDabagov, S. B., Kumakhov, M. A. & Nikitina, S. V. (1995b). Phys. Lett. A, 203, 279–282.  CrossRef CAS Web of Science Google Scholar
First citationDabagov, S. B., Kumakhov, M. A., Nikitina, S. V., Murashova, V. A., Fedorchuk, R. V. & Yakimenko, M. N. (1995a). J. Synchrotron Rad. 2, 132–135.  CrossRef CAS Web of Science IUCr Journals Google Scholar
First citationDabagov, S. B., Marcelli, A., Murashova, V. A., Svyatoslavsky, N. L., Fedorchuk, R. V. & Yakimenko, M. N. (2000). Appl. Opt. 39, 3338–3343.  Web of Science CrossRef PubMed CAS Google Scholar
First citationDabagov, S. B. & Uberall, H. (2007). Nucl. Instrum. Methods Phys. Res. A, 580, 756–763.  Web of Science CrossRef CAS Google Scholar
First citationDabagov, S. B. & Uberall, H. (2008). Nucl. Instrum. Methods Phys. Res. B, 266, 3881–3887.  Web of Science CrossRef CAS Google Scholar
First citationHiroyuki, D. (2002). Rep. Prog. Phys. 65, 1513–1576.  Google Scholar
First citationKaaret, P., Geissbühler, P., Chen, A. & Glavinas, E. (1992). Appl. Opt. 31, 7339–7343.  CrossRef CAS PubMed Web of Science Google Scholar
First citationKumakhov, M. A. & Komarov, F. F. (1990). Phys. Rep. 191, 289–350.  CrossRef CAS Web of Science Google Scholar
First citationMacDonald, C. A. (2010). X-ray Opt. Instrum. 10, 867049.  Google Scholar
First citationMacDonald, C. A. & Gibson, W. M. (2000). Handbook of Optics, Vol. III, ch. 30, pp. 30.1–30.12. McGraw-Hill.  Google Scholar
First citationMarcelli, A., Mazuritskiy, M. I., Dabagov, S. B., Hampai, D., Lerer, A. M., Izotova, E. A., D'Elia, A., Turchini, S., Zema, N., Zuccaro, F., de Simone, M., Rezvani, S. J. & Coreno, M. (2018). J. Instrum. 13, C03035.  Web of Science CrossRef Google Scholar
First citationMazuritskiy, M. I., Dabagov, S. B., Lerer, A. M., Dziedzic Kocurek, K., Sokolov, A., Coreno, M., Turchini, S., D'Elia, A., Sacchi, M. & Marcelli, A. (2017b). Nucl. Instrum. Methods Phys. Res. B, 402, 282–286.  Web of Science CrossRef CAS Google Scholar
First citationMazuritskiy, M. I., Dabagov, S. B., Marcelli, A., Lerer, A., Novakovich, A. & Dziedzic Kocurek, K. (2014). J. Opt. Soc. Am. B, 31, 2182–2187.  Web of Science CrossRef CAS Google Scholar
First citationMazuritskiy, M. I., Dabagov, S. B., Marcelli, A., Lerer, A. M. & Dziedzic Kocurek, K. (2016a). J. Synchrotron Rad. 23, 274–280.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationMazuritskiy, M. I. & Lerer, A. M. (2015). JETP Lett. 102, 483–486.  Web of Science CrossRef CAS Google Scholar
First citationMazuritskiy, M. I. & Lerer, A. M. (2017a). JETP Lett. 105, 572–576.  Web of Science CrossRef CAS Google Scholar
First citationMazuritskiy, M. I., Lerer, A. M., Dabagov, S. B., Marcelli, A., Hampai, D. & Dziedzic Kocurek, K. (2019a). J. Synch. Investig. 13, 1005–1013.  Web of Science CrossRef CAS Google Scholar
First citationMazuritskiy, M. I., Lerer, A. M., Kulov, S. K. & Samkanashvili, D. G. (2019b). J. Synch. Investig. 13, 499–507.  Web of Science CrossRef CAS Google Scholar
First citationMazuritskiy, M. I., Lerer, A. M. & Makhno, P. V. (2016b). J. Exp. Theor. Phys. 123, 942–949.  Web of Science CrossRef CAS Google Scholar
First citationMazuritskiy, M. I., Lerer, A. M., Marcelli, A., Dabagov, S. B., Coreno, M. & D'Elia, A. (2018). JETP Lett. 107, 600–605.  Web of Science CrossRef CAS Google Scholar
First citationNussey, P. (2005). PhD Thesis, University Leicester, UK.  Google Scholar
First citationPfeiffer, F., David, C., Burghammer, M., Riekel, C. & Salditt, T. (2002). Science, 297, 230–234.  Web of Science CrossRef PubMed CAS Google Scholar
First citationSokolov, A. A., Eggenstein, F., Erko, A., Follath, R., Künstner, S., Mast, M., Schmidt, J. S., Senf, F., Siewert, F., Zeschke, Th. & Schäfers, F. (2014). Proc. SPIE, 9206, 92060J.  Google Scholar
First citationSokolov, A., Sertsu, M. G., Gaupp, A., Lüttecke, M. & Schäfers, F. (2018). J. Synchrotron Rad. 25, 100–107.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationSun, T., Zhang, M., Liu, Z., Zhang, Z., Li, G., Ma, Y., Du, X., Jia, Q., Chen, Y., Yuan, Q., Huang, W., Zhu, P. & Ding, X. (2009). J. Synchrotron Rad. 16, 116–118.  Web of Science CrossRef IUCr Journals Google Scholar
First citationTsuji, K., Injuk, J. & Van Grieken, R. (2004). Editors. X-ray Spectrometry: Recent Technological Advances. John Wiley & Sons Ltd.  Google Scholar

This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775
Follow J. Synchrotron Rad.
Sign up for e-alerts
Follow J. Synchrotron Rad. on Twitter
Follow us on facebook
Sign up for RSS feeds