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

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775

Partially coherent light propagation through a kinoform lens

crossmark logo

aShanghai Institute of Applied Physics, Chinese Academy of Sciences, Jialuo Road 2019, Jiading District, Shanghai 201800, People's Republic of China, bUniversity of Chinese Academy of Sciences, Yuquan Road 19, Shijingshan District, Beijing 100049, People's Republic of China, and cShanghai Advanced Research Institute, Chinese Academy of Sciences, Zhangheng Road 239, Pudong District, Shanghai 201204, People's Republic of China
*Correspondence e-mail: wangyong@zjlab.org.cn, mengxiangyu@zjlab.org.cn, tairenzhong@zjlab.org.cn

Edited by Y. Amemiya, University of Tokyo, Japan (Received 5 November 2022; accepted 31 January 2023; online 22 March 2023)

Combining wave optics propagation and geometric ray tracing, the mutual optical intensity (MOI) model is extended to quantitatively simulate the propagation of partially coherent light through a kinoform lens at high speed. The MOI model can provide both a high accuracy and a high efficiency simulation. The intensity and coherence degree distributions at the focal plane are calculated using the MOI model. It is beneficial to improve the focusing capability of the kinoform lens by reducing the coherence or increasing the number of lens steps. In addition, increasing the number of steps is also beneficial to increase the photon flux and reduce the depth of focus.

1. Introduction

With the development of free-electron lasers (Emma et al., 2010[Emma, P., Akre, R., Arthur, J., Bionta, R., Bostedt, C., Bozek, J., Brachmann, A., Bucksbaum, P., Coffee, R., Decker, F. J., Ding, Y., Dowell, D., Edstrom, S., Fisher, A., Frisch, J., Gilevich, S., Hastings, J., Hays, G., Hering, P., Huang, Z., Iverson, R., Loos, H., Messerschmidt, M., Miahnahri, A., Moeller, S., Nuhn, H. D., Pile, G., Ratner, D., Rzepiela, J., Schultz, D., Smith, T., Stefan, P., Tompkins, H., Turner, J., Welch, J., White, W., Wu, J., Yocky, G. & Galayda, J. (2010). Nat. Photon. 4, 641-647.]; Ishikawa et al., 2012[Ishikawa, T., Aoyagi, H., Asaka, T., Asano, Y., Azumi, N., Bizen, T., Ego, H., Fukami, K., Fukui, T., Furukawa, Y., Goto, S., Hanaki, H., Hara, T., Hasegawa, T., Hatsui, T., Higashiya, A., Hirono, T., Hosoda, N., Ishii, M., Inagaki, T., Inubushi, Y., Itoga, T., Joti, Y., Kago, M., Kameshima, T., Kimura, H., Kirihara, Y., Kiyomichi, A., Kobayashi, T., Kondo, C., Kudo, T., Maesaka, H., Maréchal, X. M., Masuda, T., Matsubara, S., Matsumoto, T., Matsushita, T., Matsui, S., Nagasono, M., Nariyama, N., Ohashi, H., Ohata, T., Ohshima, T., Ono, S., Otake, Y., Saji, C., Sakurai, T., Sato, T., Sawada, K., Seike, T., Shirasawa, K., Sugimoto, T., Suzuki, S., Takahashi, S., Takebe, H., Takeshita, K., Tamasaku, K., Tanaka, H., Tanaka, R., Tanaka, T., Togashi, T., Togawa, K., Tokuhisa, A., Tomizawa, H., Tono, K., Wu, S. K., Yabashi, M., Yamaga, M., Yamashita, A., Yanagida, K., Zhang, C., Shintake, T., Kitamura, H. & Kumagai, N. (2012). Nat. Photon. 6, 540-544.]; Hettel, 2014[Hettel, R. (2014). J. Synchrotron Rad. 21, 843-855.]) and diffraction-limited storage rings (Allaria et al., 2012[Allaria, E., Appio, R., Badano, L., Barletta, W. A., Bassanese, S., Biedron, S. G., Borga, A., Busetto, E., Castronovo, D., Cinquegrana, P., Cleva, S., Cocco, D., Cornacchia, M., Craievich, P., Cudin, I., D'Auria, G., Dal Forno, M., Danailov, M. B., De Monte, R., De Ninno, G., Delgiusto, P., Demidovich, A., Di Mitri, S., Diviacco, B., Fabris, A., Fabris, R., Fawley, W., Ferianis, M., Ferrari, E., Ferry, S., Froehlich, L., Furlan, P., Gaio, G., Gelmetti, F., Giannessi, L., Giannini, M., Gobessi, R., Ivanov, R., Karantzoulis, E., Lonza, M., Lutman, A., Mahieu, B., Milloch, M., Milton, S. V., Musardo, M., Nikolov, I., Noe, S., Parmigiani, F., Penco, G., Petronio, M., Pivetta, L., Predonzani, M., Rossi, F., Rumiz, L., Salom, A., Scafuri, C., Serpico, C., Sigalotti, P., Spampinati, S., Spezzani, C., Svandrlik, M., Svetina, C., Tazzari, S., Trovo, M., Umer, R., Vascotto, A., Veronese, M., Visintini, R., Zaccaria, M., Zangrando, D. & Zangrando, M. (2012). Nat. Photon. 6, 699-704.]; Eriksson et al., 2014[Eriksson, M., van der Veen, J. F. & Quitmann, C. (2014). J. Synchrotron Rad. 21, 837-842.]; de Jonge et al., 2014[Jonge, M. D. de, Ryan, C. G. & Jacobsen, C. J. (2014). J. Synchrotron Rad. 21, 1031-1047.]), the coherence of synchrotron radiation X-rays has greatly improved (Qi et al., 2014[Qi, J.-C., Ye, L.-L., Chen, R.-C., Xie, H.-L., Ren, Y.-Q., Du, G.-H., Deng, B. & Xiao, T.-Q. (2014). Acta Phys. Sin. 63, 104202.]), which can advance the performance of coherent diffraction imaging (Liang et al., 2015[Liang, M., Williams, G. J., Messerschmidt, M., Seibert, M. M., Montanez, P. A., Hayes, M., Milathianaki, D., Aquila, A., Hunter, M. S., Koglin, J. E., Schafer, D. W., Guillet, S., Busse, A., Bergan, R., Olson, W., Fox, K., Stewart, N., Curtis, R., Miahnahri, A. A. & Boutet, S. (2015). J. Synchrotron Rad. 22, 514-519.]), X-ray photon correlation spectroscopy (Stephenson et al., 2009[Stephenson, G. B., Robert, A. & Grübel, G. (2009). Nat. Mater. 8, 702-703.]), coherent X-ray scatter imaging (Ice et al., 2011[Ice, G. E., Budai, J. D. & Pang, J. W. L. (2011). Science, 334, 1234-1239.]) and so on. These techniques not only rely on advanced X-ray sources but also on high efficiency focusing optical elements. Refractive focusing elements can perform submicrometre focusing in the hard X-ray energy range while keeping the characteristics of high theoretical working energy and flexibility (Snigirev & Snigireva, 2008[Snigirev, A. & Snigireva, I. (2008). C. R. Phys. 9, 507-516.]). However, due to the strong absorption of X-rays by the edges of conventional X-ray refractive focusing elements such as compound refractive lenses, the photon flux and focusing capability cannot be improved by increasing the geometric aperture. Based on the plano-concave lens, the kinoform lens is realised by removal of passive material where the path length is an integer multiple of the X-ray wavelength, which effectively increases the geometric aperture (Aristov et al., 2000[Aristov, V., Grigoriev, M., Kuznetsov, S., Shabelnikov, L., Yunkin, V., Weitkamp, T., Rau, C., Snigireva, I., Snigirev, A., Hoffmann, M. & Voges, E. (2000). Appl. Phys. Lett. 77, 4058-4060.]). At present many models can be used to optimize the design of kinoform lenses, but with some disadvantages. For example, the thin-lens approximation method (Buralli et al., 1989[Buralli, D. A., Morris, G. M. & Rogers, J. R. (1989). Appl. Opt. 28, 976-983.]) is a simple model with low calculation accuracy. The Takagi–Taupin description (TTD) (Yan et al., 2007[Yan, H. F., Maser, J., Macrander, A., Shen, Q., Vogt, S., Stephenson, G. B. & Kang, H. C. (2007). Phys. Rev. B, 76, 115438.]) of X-ray dynamical diffraction theory is a rigorously accurate method that is only suitable for the analysis of short kinoform lenses (Yan, 2010[Yan, H. F. (2010). Phys. Rev. B, 81, 075402.]). The beam propagation method (Van Roey et al., 1981[Van Roey, J., van der Donk, J. & Lagasse, P. E. (1981). J. Opt. Soc. Am. 71, 803-810.]) is based on dividing the object into slices, and angular spectral theory has high calculation accuracy and low calculation efficiency (Yan, 2010[Yan, H. F. (2010). Phys. Rev. B, 81, 075402.]). However, the existing model is only suitable for fully coherent light, and cannot simulate the propagation of partially coherent light through a kinoform lens.

Based on statistical optics, we have established the MOI model (Meng et al., 2015[Meng, X., Xue, C., Yu, H., Wang, Y., Wu, Y. & Tai, R. (2015). Opt. Express, 23, 29675-29686.]; Ren et al., 2019[Ren, J., Wang, Y., Meng, X., Shi, X., Assoufid, L. & Tai, R. (2019). J. Synchrotron Rad. 26, 1198-1207.]) that simulates the propagation of partially coherent light through synchrotron radiation beamlines. The model can provide the intensity, coherence and phase information of the partially coherent light at a specified position on the beamline. In this work, combining wave optics propagation and geometric ray tracing, the MOI model is extended to simulate the propagation of partially coherent light through a kinoform lens. The effects of the coherence and the number of kinoform lens steps on the focusing capability are quantitatively analyzed.

2. Model description

2.1. Establishing the MOI model for a kinoform lens

Since the long kinoform lens is widely applied in hard X-ray submicrometre focusing, we take the long kinoform lens as an example to introduce the MOI model. The refractive index of a plano-concave lens is given by n = 1 − δ + iβ, where δ is the refractive coefficient of the medium and β is the attenuation coefficient of the medium. X-rays can be focused by the plano-concave lens. Based on the plano-concave lens, the kinoform lens is realised by removal of passive material where the path length is an integer multiple of the X-ray wavelength, which not only maintains the focusing characteristics but also improves the X-ray transmittance. The step length, l = mλ/δ, generates an m-integer-2π phase difference between X-ray propagation inside and outside the step, where λ is the wavelength and m is a positive integer. The curved surface of the long kinoform lens is elliptical and can be expressed as (Evans-Lutterodt et al., 2003[Evans-Lutterodt, K., Ablett, J. M., Stein, A., Kao, C. C., Tennant, D. M., Klemens, F., Taylor, A., Jacobsen, C., Gammel, P. L., Huggins, H., Bogart, G., Ustin, S. & Ocola, L. (2003). Opt. Express, 11, 919-926.])

[{{x^2\,(2-\delta)}\over{f^{2}\delta}} \ + \ {{\big[z\,(2-\delta)-f\big]^{\,2}}\over{f^{2}}} \ = \ 1, \eqno(1)]

where x is the transverse coordinate of the lens, z is the coordinate along the optical axis, and f is the focal length. The MOI model uses the mutual optical intensity to describe the partially coherent light. The X-ray transmission from the upstream source to the incident surface of the kinoform lens and from the exit surface to the focal plane are free-space propagation. The mutual optical intensity propagation through free space can be written as (Born et al., 1959[Born, M., Clemmow, P. C., Gabor, D., Stokes, A. R., Taylor, A. M., Wayman, P. A. & Wilcock, W. L. (1959). Principles of Optics - Electromagnetic Theory of Propagation, Interference and Diffraction of Light, 1st ed. Pergamon Press.]; Goodman, 2015[Goodman, J. W. (2015). Statistical Optics. John Wiley & Sons.])

[\eqalignno{ J\,\left({{w_1},{w_2}}\right) \ = {}& \int\!\!\!\int J\left({{x_1},{x_2}}\right) \exp\left[-\,i\,\,{{2\pi}\over{\lambda}}\,\left(r_2-r_1\right)\right] \cr& \times {{\chi\left(\gamma_1\right)}\over{\left({\lambda}r_1\right)^{1/2}}} {{\chi\left(\gamma_2\right)}\over{\left({\lambda}r_2\right)^{1/2}}} \,{\rm{d}}x_1\,{\rm{d}}x_2, & (2)}]

where λ is the wavelength, x1 and x2 are any two points on the object plane, w1 and w2 are any two points on the image plane, r1 and r2 are the x1-to-w1 and x2-to-w2 distances, respectively, χ(γ1) and (γ2) are the inclination factors for the inclination angle γ1 and γ2, respectively, and J(x1, x2) and J(w1, w2) are the mutual optical intensity at the object and the image planes, respectively. A schematic diagram of the mutual optical intensity propagation through free space is shown in Fig. 1[link].

[Figure 1]
Figure 1
Schematic diagram of mutual optical intensity propagation through free space.

The mutual optical intensity propagation through free space based on formula (2)[link] is numerically calculated. The MOI model calculation procedure is built as follows. Firstly, the source plane is divided into many small elements. Each element is small enough to be considered to have full coherence and constant intensity. Secondly, the propagation of mutual optical intensity in each element is carried out using the Fresnel approximation or the Fraunhofer approximation (Born et al., 1959[Born, M., Clemmow, P. C., Gabor, D., Stokes, A. R., Taylor, A. M., Wayman, P. A. & Wilcock, W. L. (1959). Principles of Optics - Electromagnetic Theory of Propagation, Interference and Diffraction of Light, 1st ed. Pergamon Press.]). Finally, the mutual optical intensity at the image plane can be obtained by summing the contributions of all elements. The mutual optical intensity propagation through free space can be expressed as (Ren et al., 2019[Ren, J., Wang, Y., Meng, X., Shi, X., Assoufid, L. & Tai, R. (2019). J. Synchrotron Rad. 26, 1198-1207.])

[J\,\left(w_1,w_2\right) \ =\, \sum \left\{ A\left(x_2,w_2\right)^{*}\, \left[\sum J\left(x_1,x_2\right) \,A\left(x_1,w_1\right)\right] \right\}, \eqno(3)]

where A(x1, w1) is the integral of x1 on the object plane as follows,

[\eqalignno{ A\left(x_1,w_1\right) \ = {}& \int \exp\left[ i\,{{2\pi}\over{\lambda}} \left(r_{01}-{{w_1x}\over{r_{01}}}+{{x^2}\over{2r_{01}}}\right) + i\varphi(x) \right] \cr& \times {{\chi\left(\gamma_1\right)}\over{\left({\lambda}r_{01}\right)^{1/2}}} \,{\rm{d}}x, & (4)}]

where r01 is the distance between the center points of the two elements x1 and w1, and φ(x) is the phase distribution within each element.

The incident surface (red lines) is defined as the upstream surface of the kinoform lens, as shown in Fig. 2[link]. The exit surface (blue curve) is defined as the back surface of the kinoform lens. The incident and exit surfaces are described by the P and Q planes. The geometrical tracing method is used to analyze the wavefront propagation through the incident surface to the exit surface. The wavefront at the incident surface is divided into many small surface elements. The rays travel along the direction defined by the phase gradient within each element at the incident plane. Following Fermat's principle (Born & Wolf, 1999[Born, M. & Wolf, E. (1999). Principles of Optics - Electromagnetic Theory of Propagation, Interference and Diffraction of Light, 7th ed. Pergamon Press.]), the rays enter the kinoform lens medium, transmit through the medium, and hit the exit plane. The path length of the ray traveling inside the kinoform lens is used to describe the phase change. Assume that a certain ray enters the medium from point P at the incident surface and leaves the medium from point Q at the exit surface. J(P1, P2) is the mutual intensity of the incident surface, and J(Q1, Q2) is the mutual intensity of the exit surface, which can be expressed as

[\eqalignno{ J\left(Q_1,Q_2\right) = {}& \,J\left(P_1,P_2\right) \exp\left\{i\,{{2\pi}\over{\lambda}}\big[\Gamma\left(P_1,Q_1\right)-\Gamma\left(P_2,Q_2\right)\big]\right\} \cr& \times t\,\left(P_1,Q_1\right)\,\,t\,\left(P_2,Q_2\right), & (5)}]

[t\,\left(P,Q\right) = \left[ {{{\rm{d}}x\left(P\right)}\over{{\rm{d}}x\left(Q\right)}} \right]^{1/2}, \eqno(6)]

[\Gamma\left(P,Q\right) = \tau\left(P,Q\right) \, \left(1-\delta+i\beta\right), \eqno(7)]

where t(P, Q) is the complex amplitude transmission function; x(P) and x(Q) are the coordinates of P and Q, respectively; τ(P, Q) is the distance from point P to point Q; [\Gamma\left(P,Q\right)] is the path length from point P to point Q, which describes the phase change and attenuation of the wavefront by the medium. Since the wavefront at the incident surface of the kinoform lens is not a plane wave, some rays may escape from the step edge of the kinoform lens. For simplicity, the escaped rays are ignored in the subsequent propagation. The total lost flux in the whole kinoform is less than 0.1%. Therefore, the MOI model still has high accuracy. The real component 1−δ of the refractive index for the silicon medium is very close to 1. The phase shift [\varphi] = [(2\pi/\lambda)\,\delta\Delta{h}] is very small from the figure error or roughness Δh. Unlike reflective optical elements such as Kirk­patrick–Baez mirrors, the figure error and roughness on the kinoform lens have negligible effect on damage of the wavefront. For other refractive optical elements such as compound refractive lenses, similar conclusions have also been given (Pantell et al., 2001[Pantell, R. H., Feinstein, J., Beguiristain, H. R., Piestrup, M. A., Gary, C. K. & Cremer, J. T. (2001). Rev. Sci. Instrum. 72, 48-52.]). In this paper, the effect from optics-related errors is not considered.

[Figure 2]
Figure 2
Schematic diagram of a long kinoform lens. The x and z axes denote the transversal and longitudinal directions, L denotes the total length of the lens, R denotes the aperture of the lens, l denotes the length of each step, and P and Q denote points at the incident and exit surfaces.

Unlike the long kinoform lens, the steps for the short kinoform lens are designed for one plane, as shown in Fig. 3[link]. The thickness of the short kinoform lens is l = mλ/δ. The left plane (red line) is the incident plane, and the right curved surface (blue curve) is the exit surface. The exit surface function of different steps for the short kinoform lens is hyperbolic, while the surface function for the long kinoform is elliptical. The formula for the short kinoform lens can be expressed as follows (Cao et al., 2016[Cao, J., Wang, Y., Zou, Y., Zhang, X., Wu, Y. & Tai, R. (2016). J. Synchrotron Rad. 23, 436-442.]),

[z_k\left(x\right) \,=\, {{\left(x^2+f^{\,2}\,\right)^{1/2}\,-\,\,f-\left(k-1\right)\,m}\over{\delta}} \qquad \left(x_{k-1}\,\lt\,\,x\,\lt\,\,x_k\right), \eqno(8)]

where k is the step number. The wavefront propagation through the short kinoform lens is also defined by formula (5)[link].

[Figure 3]
Figure 3
Schematic diagram of a short kinoform lens. The x and z axes denote the transversal and longitudinal directions, R denotes the aperture of the lens, l denotes the total length of the lens, P and Q denote points at the incident and exit surfaces.

The propagation of mutual optical intensity through the lens can be performed in three steps: firstly, free-space propagation from the source plane to the incident plane using formula (3)[link]; secondly, propagation from incident surface to the exit surface using formula (5)[link]; and, finally, propagation from the exit surface to the focal plane using formula (3)[link].

2.2. Gaussian Schell model

The mutual optical intensity at the source plane can be described by the Gaussian Schell model (GSM) (Gori et al., 2001[Gori, F., Santarsiero, M., Piquero, G., Borghi, R., Mondello, A. & Simon, R. (2001). J. Opt. A: Pure Appl. Opt. 3, 1-9.]; Starikov & Wolf, 1982[Starikov, A. & Wolf, E. (1982). J. Opt. Soc. Am. 72, 923-928.]) and is given as

[J\left(x_1,x_2\right) = I_0 \exp\left(-{{\left(x_1-x_2\right)^2}\over{2\chi^2}}\right) \exp\left(- {{x_1^2+x_2^2}\over{4\sigma^2}}\right), \eqno(9)]

where J(x1, x2) is the mutual optical intensity between two points x1 and x2 at the source plane, χ represents the coherence length, σ represents the root-mean-square (r.m.s.) beam size and I0 represents the optical intensity at the central point. Note that the GSM is not necessary for the MOI code, and any mutual optical intensity distribution can be used for the calculation. The global degree of coherence of the source is described by the following formula (Vartanyants & Singer, 2010[Vartanyants, I. A. & Singer, A. (2010). New J. Phys. 12, 035004.]),

[G = {{ \textstyle\int\!\!\int|\,J\left(x_1,x_2\right)\,|^2\,{\rm{d}}x_1\,{\rm{d}}x_2 }\over{ \left[\textstyle\int{I}\left(x\right)\,{\rm{d}}x\right]^2 }}, \eqno(10)]

where I(x) denotes the intensity distribution at the source plane.

3. Simulation results

3.1. Simulation of partially coherent light propagation through a long kinoform lens

The MOI model is used to simulate the partially coherent light propagation through a long kinoform lens. The optical setup is shown in Fig. 4[link]. The mutual optical intensity at the source plane is obtained by using formula (9)[link]. The GSM source energy is 10 keV with beam r.m.s. size of σ = 10 µm and coherence length of χ = 11.5 µm, and the global degree of coherence G is 0.50. The source-to-kinoform and kinoform-to-image distances are 100.013 m and 11.111 m, respectively. The aperture of the kinoform lens R is 3.18 mm. The total length of the lens L is 25.96 mm with each step length 203 µm; the step number of the lens is 128, and the material of the lens is silicon. When the energy is 10 keV, δ and β are 1.2398 × 10−6 and 7.3841 × 10−8, respectively. The focal length of the lens is 10 m.

[Figure 4]
Figure 4
Optical setup for the long kinoform lens.

The incident surface is divided into 5000 elements with each element of size 0.64 µm. The intensity and coherence degree distribution at the focal plane are shown in Figs. 5[link](a) and 5[link](b). The focus spot r.m.s. size is 1.17 µm, calculated using the MOI model, which is very close to the spot size of 1.11 µm estimated by the geometrical optics. The reason for the slight difference in the spot size is that the geometrical optics cannot simulate the diffraction effect of the kinoform lens structure on the focus spot. The X-ray transmission varies from 100% to 31% due to the stepped structure of the kinoform lens. The diffraction effect from the stepped structure results in a slight enhancement of the global degree of coherence to 0.55 at the focal plane. The apparent oscillation peaks on the edges of the coherence degree profile are caused by the diffraction effect of the kinoform lens stepped structure. The MOI model has high calculation efficiency. It can be used to calculate the intensity distribution along the optical axis, as shown in Fig. 5[link](c). The focus depth is approximately 240 mm. The minimal focus spot located at −18 mm is due to the flat plane wave of the source and the diffraction effect of the kinoform lens stepped structure. In our previous paper (Sun et al., 2022[Sun, W., Wang, Y., Meng, X., Ren, J., Cao, J., Li, J. & Tai, R. (2022). J. Synchrotron Rad. 29, 1338-1343.]), the ladder wavefront (LWF) model can be seen as a special MOI model under the full coherence condition. The spot intensity at the focal plane calculated from the LWF model is in good agreement with that calculated from the beam propagation method model. Therefore, the MOI model has high accuracy.

[Figure 5]
Figure 5
Partially coherent light propagation through the long kinoform lens. (a) Intensity distribution at the focal plane, (b) coherence degree distribution at the focal plane, (c) intensity distribution along the focus depth. The coherence degree denotes the correlation between any point and the center point.

The source spot size is σ = 10 µm. The coherence length is chosen to be 5.2 µm, 11.5 µm, 22.7 µm and ∞, with corresponding global degree of coherence of 0.25, 0.5, 0.75 and 1. The intensity distribution at the focal plane corresponding to the minimum focus spot is shown in Fig. 6[link](a). When the global degree of coherence is 0.25, the spot r.m.s. size at the focal plane is 1.12 µm. With increasing coherence length, the spot size at the focal plane increases gradually. When the global degree of coherence is 1, the spot r.m.s. size at the focal plane is 1.21 µm. The source divergence angle is inversely proportional to the global degree of coherence (Vartanyants & Singer, 2010[Vartanyants, I. A. & Singer, A. (2010). New J. Phys. 12, 035004.]). To make the full beam accepted by the kinoform lens with aperture of 3.18 mm, the minimum global degree of coherence should be greater than 0.2. The spot size as a function of the global degree of coherence is shown in Fig. 6[link](b). The greater the global degree of coherence, the stronger the diffraction effect from the stepped structure. Therefore, the spot size at the focal plane gradually rises with increasing global degree of coherence.

[Figure 6]
Figure 6
(a) Intensity distribution at the focal plane with different global degrees of coherence. (b) Spot size at the focal plane as a function of global degree of coherence.

The kinoform lens is composed of many steps. The number of steps is related to the thickness of the steps. A fewer number of steps reduces the light transmission of the lens, while a higher number of steps greatly increases the difficulty of kinoform lens fabrication. At present the shortest processing length of kinoform lens steps is 0.4 µm (Xu et al., 2018[Xu, J., Liu, G., Huang, Q., Liu, M., Zhou, X., Wu, H., Li, N., Li, Y., Xu, X., Liang, D., Jiang, H., Zhang, L., Wang, J., Cao, J. & Lin, H. (2018). J. Instrum. 13, C07005.]); therefore, the effect of the number of steps on the focusing capability of the long kinoform lens needs to be analyzed.

The GSM source parameters are σ = 10 µm and χ = 11.5 µm, and the global degree of coherence is G = 0.50. The source energy is 10 keV. The kinoform lens has the same specific parameters as in the previous paragraph. We analyze the effect of the number of steps k varying from 1 to 1024 on the focusing capability. Peak intensity, spot r.m.s. size, focal depth, global degree of coherence and normalized photon flux at the position of the minimum focus spot are shown in Figs. 7[link](a)–7[link](e). The normalized photon flux is the ratio between the photon flux at the focal plane and the source plane. As the number of steps increases and the length of the steps decreases, the light transmission difference for various positions of each step gradually decreases. As the number of steps k varies from 1 to 1024, the minimum light transmittance of the kinoform lens stepped structure increases from 0% to 86%. Therefore, the diffraction effect of the steps decreases with increasing step number, which is beneficial to increasing the focusing capability and reducing the focus depth. Figs. 7[link](a), 7[link](b) and 7[link](e) show that the peak intensity increases, the spot size decreases and the normalized photon flux increases with increasing step number. When the step number is 128, the light transmission changes with the position of one step from 31% to 100%. When the step number is greater than 128, the spot size r.m.s. of 1.14 µm tends to be stable. At the same time, the focus depth as a function of step number is shown in Fig. 7[link](c). When the number of steps is 128, the focal depth is 240 mm, which tends to be stable. Since the diffraction effect of the steps decreases with increasing step number, the global degree of coherence decreases with increasing step number as shown in Fig. 7[link](d). When the number of steps is 128, the global degree of coherence at the focal plane is 0.55, which is close to stable. Increasing the number of steps is beneficial to improve the focusing capability, increase the photon flux and reduce the focal depth, but it also increases the manufacturing difficulty. Selecting the appropriate number of steps is very important for the design of a kinoform lens. In this case, our analysis shows that 128 steps are good enough to achieve a high-focusing performance. All MOI simulations were performed on a laptop with i5-9300h CPU and 16 GB RAM. The element numbers of the wavefront for different planes are chosen to be 5000, which is large enough to achieve high accuracy calculation. The total calculation time from the source to focal planes is 4 s. Therefore, the MOI model can provide both a high accuracy and a high efficiency simulation.

[Figure 7]
Figure 7
Focusing performance of a long kinoform lens as a function of the number of steps. (a) Peak intensity, (b) spot r.m.s. size, (c) focal depth, (d) global degree of coherence and (e) normalized photon flux with the various numbers of steps.

3.2. Simulation of partially coherent light propagation through a short kinoform lens

The MOI model is used to simulate the propagation of partially coherent light through a short kinoform lens. The parameters of the source are the same as those given in Section 3.1[link]. The optical setup is shown in Fig. 8[link]. The source-to-lens and lens-to-focal distances are 100 m and 11.111 m, respectively. The aperture of the kinoform lens, R, is 3.19 mm, the total length of the lens, l, is 203 µm and the step number of the lens is 128. The focal length of the lens is 10 m. The intensity and coherence degree distribution at the focal plane are shown in Figs. 9[link](a) and 9[link](b). The focus spot r.m.s. size is 1.14 µm which is close to the spot size of 1.11 µm estimated by the geometrical optics. The spot sizes are slightly different for the long and short kinoform lens. This is due to the difference between the long and short kinoform lens. The short kinoform lens is more like an ideal lens than the long one because its length along the beam is shorter than that of the long one. Therefore, the actual focus point of the long kinoform lens is farther away from the theoretical focus than the short one. We calculated the intensity profile at the theoretical position, so the focal spot size of the short kinoform lens is smaller than that of the long one. Like the long kinoform lens, the stepped structure of the kinoform lens causes a diffraction effect. This diffraction effect results in a slight enhancement of the global degree of coherence of 0.55 at the focal plane. The intensity distribution along the focus depth is shown in Fig. 9[link](c). The focus depth is 240 mm. There is a 3.6 mm position difference between the simulated and geometrical focal planes.

[Figure 8]
Figure 8
Optical setup for the short kinoform lens.
[Figure 9]
Figure 9
Partially coherent light propagation through the short kinoform lens. (a) Intensity distribution at the focal plane, (b) coherence degree distribution at the focal plane and (c) intensity distribution along the focus depth.

The source spot size is σ = 10 µm. The coherence length is chosen to be 5.2 µm, 11.5 µm, 22.7 µm and ∞, with corresponding global degree of coherence of 0.25, 0.5, 0.75 and 1. The intensity distribution at the plane for the minimum focus spot is shown in Fig. 10[link](a). The focus spot size with various global degrees of coherence is shown in Fig. 10[link](b). Like the long kinoform lens, the focus spot size gradually increases as the global degree of coherence increases.

[Figure 10]
Figure 10
(a) Intensity distribution at the focal plane with different global degrees of coherence. (b) Focus spot size for various global degrees of coherence.

We simulate the effect of the number of steps on the focusing capability of the short kinoform lens. The parameters of the source are the same as those given in Section 3.1[link]. The aperture of the lens is R = 3.19 mm, the focal length is f = 10 m and the total length of the lens is L = 203 µm. Peak intensity, spot r.m.s. size, focal depth, global degree of coherence and normalized photon flux at the focal plane with various step numbers k are shown in Figs. 11[link](a)–11(e). Like the long kinoform lens, increasing the number of steps is beneficial to improve the focusing capability, increase the photon flux and reduce the focal depth. When the number of steps is larger than 128, the spot r.m.s. size, focal depth and global degree of coherence at the focal plane tend to be stable. Therefore, 128 steps are sufficient to achieve high focusing performance.

[Figure 11]
Figure 11
Focusing performance of the short kinoform lens with various numbers of steps. (a) Peak intensity, (b) spot r.m.s. size, (c) focal depth, (d) global degree of coherence and (e) normalized photon flux as a function of the number of steps at the focal plane.

4. Conclusion

In this paper, the MOI model is developed to simulate the propagation of partially coherent light through the long and short kinoform lenses with high speed. The simulation can be performed within several seconds which is practically useful for the optical design of kinoform lenses. Reducing the coherence is beneficial to improve the focusing capability of the kinoform lens. Increasing the number of steps is beneficial to improve the focusing capability of the kinoform lens, increase the photon flux and reduce the focus depth. The MOI model is a useful tool for the design of kinoform lenses in synchrotron radiation beamlines.

Funding information

The following funding is acknowledged: the Nation Key R&D Program of China (grant No. 2021YFA1601003 and 2017YFA0403400); the National Natural Science Foundation of China (grant No. 11875314, 11805260, 11775291 and U1932116).

References

First citationAllaria, E., Appio, R., Badano, L., Barletta, W. A., Bassanese, S., Biedron, S. G., Borga, A., Busetto, E., Castronovo, D., Cinquegrana, P., Cleva, S., Cocco, D., Cornacchia, M., Craievich, P., Cudin, I., D'Auria, G., Dal Forno, M., Danailov, M. B., De Monte, R., De Ninno, G., Delgiusto, P., Demidovich, A., Di Mitri, S., Diviacco, B., Fabris, A., Fabris, R., Fawley, W., Ferianis, M., Ferrari, E., Ferry, S., Froehlich, L., Furlan, P., Gaio, G., Gelmetti, F., Giannessi, L., Giannini, M., Gobessi, R., Ivanov, R., Karantzoulis, E., Lonza, M., Lutman, A., Mahieu, B., Milloch, M., Milton, S. V., Musardo, M., Nikolov, I., Noe, S., Parmigiani, F., Penco, G., Petronio, M., Pivetta, L., Predonzani, M., Rossi, F., Rumiz, L., Salom, A., Scafuri, C., Serpico, C., Sigalotti, P., Spampinati, S., Spezzani, C., Svandrlik, M., Svetina, C., Tazzari, S., Trovo, M., Umer, R., Vascotto, A., Veronese, M., Visintini, R., Zaccaria, M., Zangrando, D. & Zangrando, M. (2012). Nat. Photon. 6, 699–704.  Web of Science CrossRef CAS Google Scholar
First citationAristov, V., Grigoriev, M., Kuznetsov, S., Shabelnikov, L., Yunkin, V., Weitkamp, T., Rau, C., Snigireva, I., Snigirev, A., Hoffmann, M. & Voges, E. (2000). Appl. Phys. Lett. 77, 4058–4060.  Web of Science CrossRef CAS Google Scholar
First citationBorn, M., Clemmow, P. C., Gabor, D., Stokes, A. R., Taylor, A. M., Wayman, P. A. & Wilcock, W. L. (1959). Principles of Optics – Electromagnetic Theory of Propagation, Interference and Diffraction of Light, 1st ed. Pergamon Press.  Google Scholar
First citationBorn, M. & Wolf, E. (1999). Principles of Optics – Electromagnetic Theory of Propagation, Interference and Diffraction of Light, 7th ed. Pergamon Press.  Google Scholar
First citationBuralli, D. A., Morris, G. M. & Rogers, J. R. (1989). Appl. Opt. 28, 976–983.  CrossRef CAS PubMed Web of Science Google Scholar
First citationCao, J., Wang, Y., Zou, Y., Zhang, X., Wu, Y. & Tai, R. (2016). J. Synchrotron Rad. 23, 436–442.  Web of Science CrossRef IUCr Journals Google Scholar
First citationEmma, P., Akre, R., Arthur, J., Bionta, R., Bostedt, C., Bozek, J., Brachmann, A., Bucksbaum, P., Coffee, R., Decker, F. J., Ding, Y., Dowell, D., Edstrom, S., Fisher, A., Frisch, J., Gilevich, S., Hastings, J., Hays, G., Hering, P., Huang, Z., Iverson, R., Loos, H., Messerschmidt, M., Miahnahri, A., Moeller, S., Nuhn, H. D., Pile, G., Ratner, D., Rzepiela, J., Schultz, D., Smith, T., Stefan, P., Tompkins, H., Turner, J., Welch, J., White, W., Wu, J., Yocky, G. & Galayda, J. (2010). Nat. Photon. 4, 641–647.  Web of Science CrossRef CAS Google Scholar
First citationEriksson, M., van der Veen, J. F. & Quitmann, C. (2014). J. Synchrotron Rad. 21, 837–842.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationEvans-Lutterodt, K., Ablett, J. M., Stein, A., Kao, C. C., Tennant, D. M., Klemens, F., Taylor, A., Jacobsen, C., Gammel, P. L., Huggins, H., Bogart, G., Ustin, S. & Ocola, L. (2003). Opt. Express, 11, 919–926.  Web of Science PubMed Google Scholar
First citationGoodman, J. W. (2015). Statistical Optics. John Wiley & Sons.  Google Scholar
First citationGori, F., Santarsiero, M., Piquero, G., Borghi, R., Mondello, A. & Simon, R. (2001). J. Opt. A: Pure Appl. Opt. 3, 1–9.  Google Scholar
First citationHettel, R. (2014). J. Synchrotron Rad. 21, 843–855.  Web of Science CrossRef IUCr Journals Google Scholar
First citationIce, G. E., Budai, J. D. & Pang, J. W. L. (2011). Science, 334, 1234–1239.  Web of Science CrossRef CAS PubMed Google Scholar
First citationIshikawa, T., Aoyagi, H., Asaka, T., Asano, Y., Azumi, N., Bizen, T., Ego, H., Fukami, K., Fukui, T., Furukawa, Y., Goto, S., Hanaki, H., Hara, T., Hasegawa, T., Hatsui, T., Higashiya, A., Hirono, T., Hosoda, N., Ishii, M., Inagaki, T., Inubushi, Y., Itoga, T., Joti, Y., Kago, M., Kameshima, T., Kimura, H., Kirihara, Y., Kiyomichi, A., Kobayashi, T., Kondo, C., Kudo, T., Maesaka, H., Maréchal, X. M., Masuda, T., Matsubara, S., Matsumoto, T., Matsushita, T., Matsui, S., Nagasono, M., Nariyama, N., Ohashi, H., Ohata, T., Ohshima, T., Ono, S., Otake, Y., Saji, C., Sakurai, T., Sato, T., Sawada, K., Seike, T., Shirasawa, K., Sugimoto, T., Suzuki, S., Takahashi, S., Takebe, H., Takeshita, K., Tamasaku, K., Tanaka, H., Tanaka, R., Tanaka, T., Togashi, T., Togawa, K., Tokuhisa, A., Tomizawa, H., Tono, K., Wu, S. K., Yabashi, M., Yamaga, M., Yamashita, A., Yanagida, K., Zhang, C., Shintake, T., Kitamura, H. & Kumagai, N. (2012). Nat. Photon. 6, 540–544.  Web of Science CrossRef CAS Google Scholar
First citationJonge, M. D. de, Ryan, C. G. & Jacobsen, C. J. (2014). J. Synchrotron Rad. 21, 1031–1047.  Web of Science CrossRef IUCr Journals Google Scholar
First citationLiang, M., Williams, G. J., Messerschmidt, M., Seibert, M. M., Montanez, P. A., Hayes, M., Milathianaki, D., Aquila, A., Hunter, M. S., Koglin, J. E., Schafer, D. W., Guillet, S., Busse, A., Bergan, R., Olson, W., Fox, K., Stewart, N., Curtis, R., Miahnahri, A. A. & Boutet, S. (2015). J. Synchrotron Rad. 22, 514–519.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationMeng, X., Xue, C., Yu, H., Wang, Y., Wu, Y. & Tai, R. (2015). Opt. Express, 23, 29675–29686.  Web of Science CrossRef CAS PubMed Google Scholar
First citationPantell, R. H., Feinstein, J., Beguiristain, H. R., Piestrup, M. A., Gary, C. K. & Cremer, J. T. (2001). Rev. Sci. Instrum. 72, 48–52.  Web of Science CrossRef CAS Google Scholar
First citationQi, J.-C., Ye, L.-L., Chen, R.-C., Xie, H.-L., Ren, Y.-Q., Du, G.-H., Deng, B. & Xiao, T.-Q. (2014). Acta Phys. Sin. 63, 104202.  Google Scholar
First citationRen, J., Wang, Y., Meng, X., Shi, X., Assoufid, L. & Tai, R. (2019). J. Synchrotron Rad. 26, 1198–1207.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationSnigirev, A. & Snigireva, I. (2008). C. R. Phys. 9, 507–516.  Web of Science CrossRef CAS Google Scholar
First citationStarikov, A. & Wolf, E. (1982). J. Opt. Soc. Am. 72, 923–928.  CrossRef Web of Science Google Scholar
First citationStephenson, G. B., Robert, A. & Grübel, G. (2009). Nat. Mater. 8, 702–703.  Web of Science CrossRef PubMed CAS Google Scholar
First citationSun, W., Wang, Y., Meng, X., Ren, J., Cao, J., Li, J. & Tai, R. (2022). J. Synchrotron Rad. 29, 1338–1343.  Web of Science CrossRef IUCr Journals Google Scholar
First citationVan Roey, J., van der Donk, J. & Lagasse, P. E. (1981). J. Opt. Soc. Am. 71, 803–810.  CrossRef Web of Science Google Scholar
First citationVartanyants, I. A. & Singer, A. (2010). New J. Phys. 12, 035004.  Web of Science CrossRef Google Scholar
First citationXu, J., Liu, G., Huang, Q., Liu, M., Zhou, X., Wu, H., Li, N., Li, Y., Xu, X., Liang, D., Jiang, H., Zhang, L., Wang, J., Cao, J. & Lin, H. (2018). J. Instrum. 13, C07005.  Web of Science CrossRef Google Scholar
First citationYan, H. F. (2010). Phys. Rev. B, 81, 075402.  Web of Science CrossRef Google Scholar
First citationYan, H. F., Maser, J., Macrander, A., Shen, Q., Vogt, S., Stephenson, G. B. & Kang, H. C. (2007). Phys. Rev. B, 76, 115438.  Web of Science CrossRef 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