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

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775

Mango wiggler as a novel insertion device providing a large and symmetrical imaging field of view

crossmark logo

aInstitute of High Energy Physics, Chinese Academy of Sciences, 19B Yuquan Road, Beijing 100049, People's Republic of China, and bUniversity of Chinese Academy of Sciences, Yuquan Road, Beijing 100049, People's Republic of China
*Correspondence e-mail: lim@ihep.ac.cn, lixiaoyu@ihep.ac.cn

Edited by Y. Amemiya, University of Tokyo, Japan (Received 28 February 2024; accepted 11 May 2024; online 21 June 2024)

A novel insertion device is introduced, designated as the Mango wiggler, designed for synchrotron radiation (SR) imaging that provides a large field of view. This innovative device is constructed from two orthogonal planar wigglers with a small difference in their period lengths, eliciting the phase difference of the magnetic fields to incrementally transitions from 0 to π/2. Such a configuration enlarges the vertical divergence of the light source, as with the horizontal divergence. The appellation `Mango wiggler' derives from the distinctive mango-shaped contour of its radiation field. A comprehensive suite of theoretical analyses and simulations has been executed to elucidate the radiation properties of the Mango wiggler, employing SPECTRA and Mathematica as calculation tools. In conjunction with the ongoing construction of the High Energy Photon Source in Beijing a practical Mango wiggler device has been fabricated for utilization in SR imaging applications. Theoretical analyses were applied to this particular Mango wiggler to yield several theoretical conclusions, and several simulations were performed according to the measured magnetic field results.

1. Introduction

The construction of the High Energy Photon Source (HEPS) (Xiaoming et al., 2014[Xiaoming, J., Jiuqing, W., Qing, Q., Yuhui, D., Weifan, S., Jian, C., Gang, X., Tiandou, H., Hu, D., Fusan, C., Fengli, L., Huihua, L., Junhui, Y., Chunhua, L., Yi, S., Jinhui, C., Haiyi, D., Quan, C., Wei, X., Ming, L., Guangcai, C., Hongwei, Z., Ye, T., Peng, L., Jing, L., Dongbai, S., Gang, L., Hong, S. & Jianshe, C. (2014). Synchrotron Radiat. News, 27(6), 27-31.]), a fourth-generation synchrotron radiation (SR) source, is currently underway. This advanced facility will feature a beamline dedicated to X-ray in-line imaging (Snigirev et al., 1995[Snigirev, A., Snigireva, I., Kohn, V., Kuznetsov, S. & Schelokov, I. (1995). Rev. Sci. Instrum. 66, 5486-5492.]) as demonstrated in Fig. 1[link]. This beamline will support the in situ, multiscale, multicontrast and nondestructive study of the structure of advanced engineering materials and components. It will also support the precise, mesoscopic, nondestructive and high-sensitivity three-dimensional imaging of complete organs. Moreover, it will support the high-sensitivity, high-penetration and nondestructive three-dimensional imaging of palaeontology samples and human remains. The typical sample sizes involved in these studies vary greatly, as does the X-ray absorption; thus, the photon energy required for ideal SR X-ray in-line imaging also greatly varies. Many important studies, such as precise mesoscopic imaging of intact large organs (Landhuis, 2017[Landhuis, E. (2017). Nature, 541, 559-561.]), require a field of view (FOV) of the order of centimetres or even tens of centimetres and a spatial resolution in the micrometre or even submicrometre range. A large FOV with a high resolution is desirable for all imaging systems. Conventional planar wigglers, which are commonly used as X-ray sources for in-line imaging, have limitations in terms of their narrow radiation FOV due to the two-dimensional deflection of electrons within the wiggler. This results in a smaller vertical divergence compared with horizontal divergence, necessitating the scanning and stitching of images for large samples (Xing et al., 2016[Xing, L., McKellar, R. C., Xu, X., Li, G., Bai, M., Persons, W. S., IV, Miyashita, T., Benton, M. J., Zhang, J., Wolfe, A. P., Yi, Q., Tseng, K., Ran, H. & Currie, P. J. (2016). Curr. Biol. 26, 3352-3360.]), which can be inefficient and complex. From Fig. 1[link]. the FOV of the in-line imaging depends on the divergence of the light source, and the spatial resolution depends on the source size observed at a point on the sample, which approximates the source size at a given angle. To address these issues, there is a clear need for a light source that can provide a larger divergence while maintaining a small enough source size to achieve real-time dynamic imaging with a large FOV.

[Figure 1]
Figure 1
Typical schematic display of SR X-ray in-line imaging.

To address the limitations of conventional planar wigglers and the need for a larger vertical divergence in SR imaging, various methods have been discussed in this article to modify the design of planar wigglers. These methods aim to expand the vertical divergence of the light source to achieve a larger FOV. One of the methods designated as the Mango wiggler with decisive advantages will be discussed in detail. The article commences with the genesis of the Mango wiggler, elucidating its conception. Subsequently, the radiation properties of the Mango wiggler are discussed, particularly its uneven angular distribution of radiation at both low and high energies. Both theoretical analyses and simulations are presented. The Mango wiggler device applied at HEPS is then presented. Theoretical analyses are applied to this practical Mango wiggler, and field measurement data of this actual device are used to perform some simulations. All the theoretical analyses and simulations are completed using SPECTRA (Tanaka, 2021[Tanaka, T. (2021). J. Synchrotron Rad. 28, 1267-1272.]) and Mathematica (Wolfram Research, 2018[Wolfram Research (2018). Mathematica, Version 11.3. Wolfram Research, Champaign, IL, USA.]).

Two distinct Mango wigglers are utilized in the study. The first one is the ideal Mango wiggler (Bx0 = By0 = 1.8 T, λux = 60.85 mm, λuy = 60 mm and Nu = 18), and the second one is the practical Mango wiggler at HEPS (Bx0 = By0 = 1 T, λux = 50.70 mm, λuy = 50 mm, Nux = 17.75 and Nuy = 18). The ideal Mango wiggler, which differs from the HEPS Mango wiggler in fundamental parameters, is used due to its stronger magnetic field. This allows it to generate more photons with higher energy, thus providing a clearer demonstration of radiation properties when using SPECTRA. The ideal Mango wiggler serves to illustrate the conceptual design and to derive theoretical analyses complemented by simulations in the discussion of radiation properties. These theoretical analyses are then applied to the HEPS Mango wiggler to yield theoretical outcomes. Additionally, several simulations are performed based on the measured magnetic field results of this practical Mango wiggler.

2. Concept of the Mango wiggler

Four methods to expand the vertical divergence of the light source are listed here. The first method is to add a bending magnet with a uniform horizontal field to a planar wiggler with a vertical field, abbreviated as the BP wiggler. The magnetic field strength of the added bending magnet should be much smaller than the peak field of the planar wiggler. Hence, the SR is mainly generated by the planar wiggler with a vertical field. Compared with the Mango wiggler presented in this section for which SR is produced by two wigglers in two directions, both the critical energy of the radiation and the flux generated in this first method are lower.

The second method is to use two orthogonal planar wigglers with the same peak magnetic field strength increasing from zero and a constant phase difference of π/2. The motion of electrons is gradually deflected from small circles to large circles along the beam axis. Thus, we call this a spiral wiggler. For instance, the magnetic fields of the two planar wigglers are defined as follows: Bx(z) = Bx0(z)sin(2πz/λux + π/2) and By(z) = By0(z)sin(2πz/λuy). The peak magnetic field strengths Bx0(z) and By0(z) are the same, and both increase from 0 to 1.8 T. The period lengths λux and λuy are both 60 mm, and the number of periods Nu is 18. The magnetic field and the angular trajectory of electrons are plotted in Fig. 2[link].

[Figure 2]
Figure 2
Magnetic field and angular trajectory of electrons in a spiral wiggler. (a) Magnetic field. (b) Magnetic field along z. (c) Angular trajectory of electrons. (d) Angular trajectory of electrons along z.

The third method is to use two orthogonal planar wigglers with a constant phase difference of π/2 but different peak magnetic field strengths. The horizontal peak magnetic field strength increases from zero while the vertical peak magnetic field strength is a constant. For instance, the magnetic fields for the two planar wigglers are Bx(z) = Bx0(z)sin(2πz/λux + π/2) and By(z) = By0sin(2πz/λuy). Bx0(z) increases from 0 to 1.8 T, and By0 is 1.8 T. The period lengths λux and λuy are both 60 mm, and the number of periods Nu is 18. The magnetic field and the angular trajectory of electrons are plotted in Fig. 3[link].

[Figure 3]
Figure 3
Magnetic field and angular trajectory of electrons in two combined wigglers with a constant phase difference of π/2 but different peak magnetic field strengths. (a) Magnetic field. (b) Magnetic field along z. (c) Angular trajectory of electrons. (d) Angular trajectory of electrons along z.

However, the magnetic field strength is too weak in the central area where x′ = 0 and y′ = 0 in Fig. 2[link](c). Therefore, both the power of the radiation and the flux of high-energy photons are very low due to the weak magnetic field strength in the centre of the FOV for the spiral wiggler. Moreover, realizing a magnetic field strength starting from a value of zero is very difficult in the design, processing, installing and adjusting.

The fourth method is still to use two orthogonal planar wigglers but for both directions the peak magnetic field strengths are the same constant value. There is a small difference in the period length for two directions, which causes the electron motion to vary from linear motion to circular motion. The magnetic fields of these combined wigglers are illustrated as equations (1)[link] and (2)[link],

[B_x\left(z\right) = B_0\sin \left({{2\pi{z}}\over{\lambda_{ux}}}\right), \qquad \lambda_{ux} = {{L}\over{N_u\pm(1/4)}}, \eqno(1)]

[B_y\left(z\right) = B_0\sin\left({{2\pi{z}}\over{\lambda_{uy}}}\right), \qquad \lambda_{uy} = {{L}\over{N_u}}. \eqno(2)]

where B0 is the constant peak magnetic field strength and L is the length of the whole magnet. For instance, the peak magnetic field strength B0 is 1.8 T, the period length λux is 60.85 mm, λuy is 60 mm and the number of periods Nu is 18 (Fig. 4[link]). The strongest combined magnetic field strength is [\sqrt{2}B_0]. Angular obits are evenly distributed. So, most of the angular space (x′–y′) is evenly covered with longitude lines (demonstrated in Fig. 12 of Section 3.2[link]). And the expanded FOV is centrally symmetrical.

[Figure 4]
Figure 4
Magnetic field and angular trajectory of electrons in the ideal Mango wiggler. (a) Magnetic field. (b) Magnetic field along z. (c) Angular trajectory of electrons. (d) Angular trajectory of electrons along z.

The Mango wiggler has the following advantages. First, the Mango wiggler has a large integral radiation power. Second, the Mango wiggler has a high critical energy. Third, the high-energy SR of the Mango wiggler is produced at the centre of the FOV. Fourth, the Mango wiggler can not only enlarge the FOV but also maintain good imaging resolution which will be discussed in detail in Section 3.4[link]. Fifth, the magnetic field of the Mango wiggler is feasible in the design, processing, installing and adjusting.

3. Radiation properties of the Mango wiggler

The Mango wiggler appears to be capable of providing a large and uniform FOV [Fig. 4[link](c)]. The radiation distribution is still uneven at low energy and high energy for different reasons. In this section, the properties of the radiation at low energy and high energy will be analysed and discussed. In addition, the photon flux, the critical energy and the source size will be discussed.

3.1. Uneven radiation distribution at low energy

For periodic wigglers, including planar wigglers and the Mango wiggler, the spectral distribution at the beam axis is discrete at low energy due to the undulator effect. Considering the coupling between the spectrum and deviation angle, the discretization of the spectrum at low energy leads to an uneven angular distribution of radiation.

First, the SR of a single electron is considered. As derived from Kim (1989[Kim, K. J. (1989). AIP Conf. Proc. 184, 565-632.]), the relative spectral line width for the nth harmonic is 1/(nNu). This line width is the distance between the main peak and closest valley according to the Rayleigh criterion (Rayleigh, 1879[Rayleigh, J. W. S. (1879). London, Edinb. Dubl. Philos. Mag. J. Sci. 8, 261-274.]). In order to unify the width criterion during analysis, the standard deviation is used to describe the width. Therefore, the relative line width for the nth harmonic, 1/(nNu), needs to be transformed into a standard deviation form. The radiation mechanism from different magnetic poles is similar to multislit diffraction. The spectral distribution is regarded as a [sin(mx)/sin(x)]2-term function based on multislit diffraction theory, where m has a physical meaning of 2Nu for an insertion device. The Rayleigh line width of [sin(mx)/sin(x)]2 is π/m [Fig. 5[link]a], and its standard deviation σi is defined in equation (3)[link],

[\sigma_i = \left( \! \Bigg\{ \int\limits_{-\pi/m}^{\,\pi/m} \! \left[{{\sin\left(mx\right)}\over{\sin\left(x\right)}}\right]^2 \! x^2\,{\rm{d}}x \Bigg\} \Bigg/ \Bigg\{\int\limits_{-\pi/m}^{\,\pi/m} \! \left[{{\sin\left(mx\right)}\over{\sin\left(x\right)}}\right]^2 {\rm{d}}x \Bigg\}\, \right)^{1/2}. \eqno(3)]

The ratio between the standard deviation and the Rayleigh line width is plotted in Fig. 5[link](b). The ratio is a constant value of 0.335 when m is large. Considering the reality of a wiggler, m is usually of the order of dozens. Thus, for the Rayleigh line width of the nth harmonic 1/(nNu), the corresponding standard deviation is 0.335/(nNu).

[Figure 5]
Figure 5
(a) The [sin(mx)/sin(x)]2 function when m = 5, for instance. (b) Ratio between the standard deviation and Rayleigh line width of [sin(mx)/sin(x)]2.

Second, the divergence of the electron bunch must be taken into account. The wavelength for the nth harmonic of a single electron observed at any deviation angle θ off the direction of the velocity of the single electron is illustrated as equation (4)[link] (Kim, 1989[Kim, K. J. (1989). AIP Conf. Proc. 184, 565-632.]), where K is the deflection parameter and γ is the particle energy in units of its rest mass. When observed along a single electron direction, where θ is zero, the wavelength for the nth harmonic is illustrated as equation (5)[link]. Therefore, the relative redshift for the nth harmonic of an electron with a deviation angle of θ when observed along the beam axis is illustrated as equation (6)[link]. The electrons distribution approximately follows a Gaussian distribution, and σθ is the standard deviation of the deviation angle θ of electrons in the electron bunch, with σθ = [(\sigma_{x^{\prime}}^{2}+\sigma_{y^{\prime}}^{2})^{1/2}], where σx is known as the horizontal divergence of the electron bunch, and σy is known as the vertical divergence of the electron bunch. Each deviating electron has a relative redshift for the nth harmonic when observed along the beam axis. The average relative redshift of electrons in the electron bunch is shown in equation (7)[link]. The degree of dispersion for all relative redshifts compared with the average relative redshift can be described by equation (8)[link]. For the Mango wiggler, the deflection parameter K is composed of deflection parameters from both the horizontal and vertical magnets, K = (Kx2+Ky2)1/2, where Kx = eBx0λux/2πmec, Ky = eBy0λuy/2πmec, e is the electron charge, me is the static mass of the electron and c is the velocity of light.

[\lambda_{n,\theta} = {{ 1+(K^2/2)+\gamma^2\theta^2 }\over{ n\,2\gamma^2 }} \, \lambda_u, \eqno(4)]

[\lambda_{n,0} = {{ 1+(K^2/2) }\over{ n\,2\gamma^2}} \, \lambda_u, \eqno(5)]

[{{\Delta\lambda_n}\over{\lambda_n}} = {{\lambda_{n,\theta}-\lambda_{n,0}}\over{\lambda_{n,0}}} = {{2\gamma^2\theta^2}\over{2+K^2}}, \eqno(6)]

[\eqalignno{ {{\overline{\Delta\lambda_n}}\over{\lambda_n}} & = {{ \int\limits_{-\infty}^{\,\infty} \left\{\left[1/\sqrt{2\pi}\sigma_{\theta}\right]\exp\left(-\theta^2/2\sigma_{\theta}^2\right)\right\} \, {{\Delta\lambda_n}\over{\lambda_n}}\,{\rm{d}}\theta }\over{ \int\limits_{-\infty}^{\,\infty} \left\{\left[1/\sqrt{2\pi}\sigma_{\theta}\right]\exp\left(-\theta^2/2\sigma_{\theta}^2\right)\right\}\,{\rm{d}}\theta }} \cr& = {{2\gamma^2\sigma_{\theta}^2}\over{2+K^2}}, & (7)}]

[\eqalignno{ \sigma_d & = \left( {{ \int\limits_{-\infty}^{\,\infty} \left\{\left[1/\sqrt{2\pi}\sigma_{\theta}\right]\exp\left(-\theta^2/2\sigma_{\theta}^2\right)\right\} \, \left({{\Delta\lambda_n}\over{\lambda_n}}-{{\overline{\Delta\lambda_n}}\over{\lambda_n}}\right)^2\,{\rm{d}}\theta } \over { \int\limits_{-\infty}^{\,\infty} \left\{\left[1/\sqrt{2\pi}\sigma_{\theta}\right]\exp\left(-\theta^2/2\sigma_{\theta}^2\right)\right\}\,{\rm{d}}\theta }} \right)^{1/2} \cr& = \,{{2\sqrt2\gamma^2\sigma_{\theta}^2}\over{2+K^2}}. & (8)}]

Third, the energy spread of the electron bunch δe is taken into account. When observed along the beam axis, the SR energy for the nth harmonic ɛn,0 is illustrated as equation (9)[link], where E is the electron energy (Kim, 1989[Kim, K. J. (1989). AIP Conf. Proc. 184, 565-632.]). Then, equations (10)[link] and (11) are derived from equation (9)[link],

[\varepsilon_{n,0}\ [{\rm{keV}}] \,=\, n \times 0.950 {{ E^{\,2}\ [{\rm{GeV}}] }\over{ \left[1+(K^2/2)\right]\,\lambda_{u}\ [{\rm{cm}}]}}, \eqno(9)]

[{{d\varepsilon_{n,0}}\over{\varepsilon_{n,0}}} \,=\, 2\,{{dE}\over{E}}, \eqno(10)]

[{{\sigma_{\varepsilon_{n,0}}}\over{\varepsilon_{n,0}}} \,=\, 2\,{{\sigma_{E}}\over{E}} \,=\, 2\delta_e. \eqno(11)]

In summary, the spectral width δn,0 for the nth harmonic along the beam axis is given by equation (12)[link], ignoring the phase errors of the magnets here,

[\delta_{n,0} = \left[ \left({{0.335}\over{nN_u}}\right)^2 + \,\left({{2\sqrt2\gamma^2\sigma_{\theta}^2}\over{2+K^2}}\right)^2 + \,\left(2\delta_{\rm{e}}\right)^2 \right]^{1/2}. \eqno(12)]

Whether the spectrum is uniformly distributed depends on the relationship between the spectral width and the distance between any two adjacent harmonics in the spectrum. When adjacent harmonics are much closer to each other, their spectral widths will overlap much more, and the spectrum is more uniform. The relative distance between any two adjacent harmonics in the spectrum is 1/n. A parameter is defined here to determine the uniformity of the spectrum, χ1, as presented in equation (13)[link]. The ideal Mango wiggler with the parameters in Section 2[link] is utilized to simulate the flux density with an energy dependence along the beam axis and an angular dependence via SPECTRA (Figs. 6[link], 7[link] and 8[link]). The spectrum and angular distribution of radiation are both uniform for those harmonics for which χ1 is less than 1. The spectrum and angular distribution of radiation are uneven for those harmonics for which χ1 is greater than 3. 1 < χ1 < 3 is the transition region from uniformity to unevenness,

[\chi_1 = {{1/n}\over{\delta_{n,0}}} = {{ 1/n }\over{ \left[ \left({{0.335}\over{nN_u}}\right)^2 +\, \left({{2\sqrt2\gamma^2\sigma_{\theta}^2}\over{2\,+\,K^2}}\right)^2 +\, \left(2\delta_{\rm{e}}\right)^2\right]^{1/2} }}. \eqno(13)]

To make the discrete spectrum uniform so that the angular distribution of radiation is uniform, breaking up the magnetic periodicity is one approach, such as through tapering (Kroll et al., 1981[Kroll, N., Morton, P. & Rosenbluth, M. (1981). IEEE J. Quantum Electron. 17, 1436-1468.]; Orzechowski et al., 1986[Orzechowski, T. J., Anderson, B. R., Clark, J. C., Fawley, W. M., Paul, A. C., Prosnitz, D., Scharlemann, E. T., Yarema, S. M., Hopkins, D. B., Sessler, A. M. & Wurtele, J. S. (1986). Phys. Rev. Lett. 57, 2172-2175.]; Boyanov et al., 1994[Boyanov, B. I., Bunker, G., Lee, J. M. & Morrison, T. I. (1994). Nucl. Instrum. Methods Phys. Res. A, 339, 596-603.]). For instance, the magnetic fields of the ideal Mango wiggler presented in Section 2[link] can be tapered. Bx0 and By0 linearly increase from 1.62 T to 1.8 T (Fig. 9[link]). Then, the spectrum and angular distribution of this tapered Mango wiggler at 8 keV are simulated (Figs. 10[link] and 11[link]). The angular distribution of radiation at 8 keV becomes more uniform based on a comparison of Fig. 8[link](b) and Fig. 11[link](b).

[Figure 6]
Figure 6
Flux density spectrum (5–30 keV) of the ideal Mango wiggler along the beam axis.
[Figure 7]
Figure 7
Flux density spectra of the ideal Mango wiggler along the beam axis: (a) SR energy around 8 keV, χ1 = 3.15. (b) SR energy around 11.36 keV, χ1 = 2.22. (c) SR energy around 16 keV, χ1 = 1.57. (d) SR energy around 23 keV, χ1 = 1.10.
[Figure 8]
Figure 8
Angular distribution of the flux density at low energy of the ideal Mango wiggler. (a) SR energy at 8 keV, χ1 = 3.15. (b) One-dimensional distribution of the flux density along the 45° white dashed line in (a). (c) SR energy at 11.36 keV, χ1 = 2.22. (d) One-dimensional distribution of the flux density along the 45° white dashed line in (c). (e) SR energy at 16 keV, χ1 = 1.57. (f) One-dimensional distribution of the flux density along the 45° white dashed line in (e). (g) SR energy at 23 keV, χ1 = 1.10. (h) One-dimensional distribution of the flux density along the 45° white dashed line in (g).
[Figure 9]
Figure 9
Magnetic fields of the tapered Mango wiggler.
[Figure 10]
Figure 10
Flux density spectrum at 8 keV of the tapered Mango wiggler.
[Figure 11]
Figure 11
Angular distribution of the flux density at 8 keV of the tapered Mango wiggler. (a) Angular distribution of radiation. (b) One-dimensional distribution of the flux density along the 45° white dashed line in (a).

3.2. Uneven distribution of radiation at high energy

As analysed above, the spectral harmonics can overlap with each other much more when the harmonics become higher for a planar wiggler. This property can ensure a uniform angular distribution of radiation for a planar wiggler at high energy but cannot guarantee that the Mango wiggler has a uniform angular distribution of radiation. In the radiation view of the Mango wiggler, the angular distribution of radiation produced by the divergence of the electron bunch and movement trajectory should be considered. Electrons in the Mango wiggler move in three-dimensional space rather than in a two-dimensional plane. There are angular gaps among the longitude lines in x′–y′ space (Fig. 12[link]). The SR divergence decreases as the radiation energy increases. At high energy, the radiation divergence produced by the electron bunch is relatively too small for the gaps to form a uniform angular distribution of radiation.

[Figure 12]
Figure 12
Demonstration of longitude lines, angular gaps among adjacent longitude lines, and divergence of the electron bunch in the 45° direction: σt.

The divergence of the SR produced by an electron bunch is the convolution between the SR produced by a single electron and the divergence of the electron bunch. The divergence of the SR produced by a single electron, σr, is illustrated in equation (14)[link] (Kim, 1989[Kim, K. J. (1989). AIP Conf. Proc. 184, 565-632.]),

[\sigma _{r^{\prime}} = \sqrt{{{2\pi}\over{3}}} \, {{1}\over{\gamma}} \, y^{-1} {{ \int_y^{\,\infty}K_{5/3}\left(y^{\prime}\right)\,{\rm{d}}y^{\prime} }\over{ K_{2/3}^2\left(y/2\right) }}. \eqno(14)]

where y equals ε/εc. εc is the critical energy of radiation, with εc = hωc/2π. ωc is the corresponding critical frequency, with ωc = 3γ3c/2ρ. ρ is the radius of instantaneous curvature of the electron trajectory, with ρ = mecγ/eB. h is the Planck constant. K5/3(y′) and K2/3(y/2) are modified Bessel functions. According to this equation, σr decreases as the SR energy increases. The divergence of the electron bunch in the 45° direction is σt = [[(\sigma_{x^{\prime}}^{2}+\sigma_{y^{\prime}}^{2})/2]^{1/2}](Fig. 12[link]), and the divergence of the SR produced by the electron bunch is σT = [(\sigma_{r^{\,\prime}}^{2}+\sigma_{t^{\,\prime}}^{2})^{1/2}]. Therefore, σT decreases as the SR energy increases.

The central angular gap Ag0 is considered representative because it is the largest among all the gaps. Ag is used to describe the gap between any adjacent longitude lines, and Ag0 is a subset of it. Since the period lengths in the two directions differ little, the phase difference of oscillation for electrons in the two directions Δφ after undergoing an average period of |λuxλuy|/[(λux + λuy)/2] is illustrated as equation (15)[link]. The angular motion formulas in the two directions are x′ = Kx/γcosφ and y′ = Ky/γcos(φ + Δφ), which are parametric functions. Thus, the central angular gap Ag0 is derived based on analytical geometry and illustrated as equation (16)[link]. Whether the angular distribution of radiation at high energy is uniform depends on the relationship between Ag0 and σT. Similar to the analysis at low energy, a parameter χ2 is defined here, which is illustrated in equation (17)[link], to determine the uniformity of the angular distribution of radiation at high energy. From the simulation results of the ideal Mango wiggler, the following conclusions can be drawn. The angular distribution of radiation can be uniform for those energies with χ2 < 1, the angular distribution of radiation can be uneven for those energies with χ2 > 3, and 1 < χ2 < 3 is the transition region from a uniform angular distribution of radiation to an uneven angular distribution of radiation. Four simulation examples of the ideal Mango wiggler at high energy are shown in Fig. 13[link]. Therefore, factors influencing uniformity at high energy can be analysed from equations (16)[link] and (17)[link]. If SR with higher critical energy is required, the magnetic field strength Bx0 and By0 must be stronger. Then this will result in a larger Ag0 in equation (16)[link] as Bx0 and By0 become larger. SR with higher energy will result in a smaller σT as mentioned before. Thus, χ2 becomes larger in equation (17)[link]. Considering this situation, in order to decrease χ2, the period length needs to be shortened. The challenge is the realization of strong magnetic fields with short periods length.

[\Delta\varphi = {{ 2\pi\left|\lambda_{ux}-\lambda_{uy}\right| }\over{ (\lambda_{ux}+\lambda_{uy})/2 }} \,=\, {{4\pi\left|\lambda_{ux}-\lambda_{uy}\right| }\over{ \lambda_{ux}+\lambda_{uy}}}, \eqno(15)]

[\eqalignno{ A_{g0} = & \Big\{ \,e\,\Big[B_{x0}^2\lambda_{ux}^2+B_{y0}^2\lambda_{uy}^2-\big(B_{x0}^4\lambda_{ux}^4+B_{y0}^4\lambda_{uy}^4\cr& +2B_{x0}^2B_{y0}^2\lambda_{ux}^2\lambda_{uy}^2\cos2\Delta\varphi\big)^{1/2}\Big]^{1/2} \Big\} \cr& \ \Big/ \big(2\sqrt2cm_{\rm{e}}\pi\gamma\big), & (16)}]

[\chi_2 = {{A_{g0}}\over{\sigma_{T^{\prime}}}}. \eqno(17)]

[Figure 13]
Figure 13
Angular distribution of the flux density at high energy of the ideal Mango wiggler. (a) SR energy at 100 keV, χ2 = 1.27. (b) One-dimensional distribution of the flux density along the 45° white dashed line in (a). (c) SR energy at 150 keV, χ2 = 1.58. (d) One-dimensional distribution of the flux density along the 45° white dashed line in (c). (e) SR energy at 267 keV, χ2 = 2.15. (f) One-dimensional distribution of the flux density along the 45° white dashed line in (e). (g) SR energy at 300 keV, χ2 = 2.29. (h) One-dimensional distribution of the flux density along the 45° white dashed line in (g).

3.3. Photon flux and critical energy

The simulated photon fluxes of the ideal Mango wiggler and a planar wiggler [By(z) = By0sin(2πz/λuy), By0 = 1.8 T, λuy = 60 mm and Nu = 18] are shown in Fig. 14[link] [curves (b) and (a)]. Another planar wiggler [By(z) = By0sin(2πz/λuy), By0 = [1.8\sqrt{2}] T, λuy = 60 mm and Nu = 18] [Fig. 14[link], curve (c)] and a helical wiggler [the magnetic field resembles a helical undulator (Alferov et al., 1974[Alferov, D. F., Bashmakov, Y. A. & Besonov, E. G. (1974). Sov. Phys. Tech. Phys. 18, 1336.])] [Bx(z) = Bx0cos(2πz/λux), By(z) = By0sin(2πz/λuy), Bx0 = By0 = 1.8 T, λux = λuy = 60 mm and Nu = 18] [Fig. 14[link], curve (d)] are also simulated. The critical energy of the planar wiggler is 3B0hγ2e/4πme (Thompson & Vaughan, 2001[Thompson, A. C. & Vaughan, D. (2001). X-ray Data Booklet. Lawrence Berkeley National Laboratory, University of California Berkeley, CA, USA.]). The critical energy of the Mango wiggler changes gradually from [3\sqrt{2}\,]B0hγ2e/4πme to 3B0hγ2e/4πme along z, and from Fig. 4[link](c) the higher critical energy appears at the central area of SR angular distribution. From the comparison, the Mango wiggler provides not only a larger FOV (see Fig. 15[link], angular distribution of the power density) but also a higher photon flux and a higher critical energy than the planar wiggler, as shown in Fig. 14[link] [curve (a)]. Moreover, the flux of the Mango wiggler is between those of the planar wiggler with a higher magnetic field, shown in Fig. 14[link] [curve (c)], and the helical wiggler, shown in Fig. 14[link] [curve (d)].

[Figure 14]
Figure 14
Photon fluxes within a 3 mrad acceptance angle for different wigglers. (a) Planar wiggler [By(z) = By0sin(2πz/λuy), By0 = 1.8 T, λuy = 60 mm and Nu = 18]. (b) The ideal Mango wiggler. (c) Planar wiggler [By(z) = By0sin(2πz/λuy), By0 = [1.8\sqrt{2}] T, λuy = 60 mm and Nu = 18]. (d) Helical wiggler [Bx(z) = Bx0cos(2πz/λux), By(z) = By0sin(2πz/λuy), Bx0 = By0 = 1.8 T, λux = λuy = 60 mm and Nu = 18].
[Figure 15]
Figure 15
Angular distribution of the power density. (a) Planar wiggler [By(z) = By0sin(2πz/λuy), By0 = 1.8 T, λuy = 60 mm and Nu = 18]. (b) The ideal Mango wiggler.

3.4. The source size at a given angle

For the applications of X-ray in-line imaging, the source size at a given angle influencing imaging resolution refers to the source size at the end of the wiggler observed from a point in the sample, shown in Fig. 16[link]. For a planar wiggler, the source size at a given angle θx has been derived from Walker (2003[Walker, R. P. (2003). Undulators, Wigglers and their Applications, edited by H. Onuki & P. Elleaume, pp. 108-147. London: Taylor & Francis.], 2017[Walker, R. P. (2017). Phys. Rev. Accel. Beams, 20, 020703.]) and illustrated as equation (18)[link]. The transverse size of the electron bunch, the transverse oscillation and the oblique observation angle contribute to the source size at a given angle. Σx is the r.m.s. of the distribution at a given angle θx, σx is the size of the electron bunch, σx is the divergence of the electron bunch, a is the amplitude of sinusoidal oscillation of the motion of the electron bunch, and L is the length of the wiggler. For a Mango wiggler, the magnet length contributing to the formation of radiation emitted at a given angle is much shorter than the total wiggler length L. Therefore, [L/2\sqrt{3}] in the third and the fourth terms on the right side of equation (18)[link] needs to be adjusted to λuσT/Ag where λu = (λux + λuy)/2 is the average period length. The source size of a Mango wiggler at a given angle (x′, y′) is illustrated in equation (19)[link], only around −x′ = y′, ax = aζx, ay = aζy, a = δλu/2π, δ = K/γ. ζx and ζy are illustrated as equations (20)[link] and (21)[link] separately. When |x′+ y′| is much larger than σT, ax and ay are both zero. Comparing equations (18)[link] and (19)[link], the source size at a given angle of the Mango wiggler is smaller than that of the planar wiggler. Thus, the Mango wiggler has more advantages in imaging resolution,

[\Sigma_x = \left(\sigma_x^2+a^2+{{L^2\sigma_{x^{\prime}}^2}\over{12}} + {{L^2\theta_x^2}\over{12}} \right)^{1/2}, \eqno(18)]

[\Sigma_{x\hbox{-}{\rm{Mango}}} = \left[ \sigma_x^2+a_x^2\left(x^{\prime},y^{\prime}\right) + {{\lambda_u^2\sigma_{T^{\prime}}^2\sigma_{x^{\prime}}^2}\over{A_g^2}} + {{\lambda_u^2\sigma_{T^{\prime}}^2x^{\prime\,2}}\over{A_g^2}} \right]^{1/2}, \eqno(19)]

[\zeta_x \simeq \left[1-\left({{x^{\prime}}\over{\delta}}\right)^2\right]^{1/2} - {{x^{\prime}}\over{\delta}}\arccos\left({{x^{\prime}}\over{\delta}}\right), \eqno(20)]

[\zeta_y \simeq \left[1-\left({{y^{\prime}}\over{\delta}}\right)^2\right]^{1/2} - {{y^{\prime}}\over{\delta}}\arccos\left({{y^{\prime}}\over{\delta}}\right). \eqno(21)]

[Figure 16]
Figure 16
The source size at the end of wiggler observed from a point in the sample.

3.5. Some comparisons with the BP wiggler

The BP wiggler has the same effect of expanding FOV. Fig. 17[link] shows the angular trajectory of the BP wiggler [Bx0 = 0.03 T, By(z) = By0sin(2πz/λuy), By0 = 1.8 T, λuy = 60 mm and Nu = 18].

[Figure 17]
Figure 17
The angular trajectory of the BP wiggler: Bx0 = 0.03 T, By(z) = By0sin(2πz/λuy), By0 = 1.8 T, λuy = 60 mm and Nu = 18.

The Mango wiggler has advantages in the uniformity of the FOV compared with the BP wiggler. If the vertical divergence is the same as that in the horizontal direction, which is 2δ, the angular gap in the central FOV of the BP wiggler will be 2δ/2Nu and the angular gap at the edges will be 2δ/Nu (Fig. 18[link]). Even if the central FOV is uniform, it does not guarantee uniformity across the entire FOV. As for the Mango wiggler, since its trajectory of the electron bunch unfolds uniformly in space, the angular gap along x′ = y′ is the same, which is 2δ/2Nu (Fig. 18[link]). If the central FOV is uniform, then the entire FOV must also be uniform.

[Figure 18]
Figure 18
A comparison of the uniformity of the FOV between the BP wiggler and Mango wiggler.

The Mango wiggler has another advantage over the BP wiggler. Due to the vertical deflection in the BP wiggler, the orbit excursion gives rise to lifetime degradation (due to effective loss of the vertical physical aperture). The Mango wiggler can avoid such an issue and thus the minimum gap can be narrower than the BP wiggler.

4. The Mango wiggler at the HEPS

The HEPS is a fourth-generation SR source that is currently under construction. We developed a Mango wiggler for X-ray in-line imaging at HEPS. The parameters of the Mango wiggler at HEPS are listed in Table 1[link]. The NdFeB magnet arrays of the delta type were adopted. This Mango wiggler works outside of vacuum. A practical Mango wiggler was constructed, as shown in Fig. 19[link], and measurements and adjustments of its magnetic field were carried out. The measured magnetic field and the trajectory at the smallest gap (12.1 mm) are shown in Fig. 20[link]. The field errors and effects on the electron trajectory can be found in Fig. 20[link].

Table 1
Parameters of the Mango wiggler designed for HEPS

Peak magnetic field strength, B0 1.0 T (0.2–1.0 T adjustable for the tapered mode)
Period lengths, λux/λuy 50.70 mm/50 mm
Period, Nux/Nuy 17.75/18
[Figure 19]
Figure 19
Photograph of the HEPS Mango wiggler.
[Figure 20]
Figure 20
Magnetic field and angular trajectory of electrons in the HEPS Mango wiggler according to the data measured at the smallest gap (12.1 mm). (a) Magnetic field. (b) Magnetic field along z. (c) Angular trajectory of electrons. (d) Angular trajectory of electrons along z.

The energy of the electrons in the storage ring is 6 GeV; the divergences of the electron bunch σx and σy are 3.14 µrad and 1.20 µrad, respectively, in the middle of the insertion device. The energy spread δe is 0.0011. An SR energy range no higher than 300 keV is chosen. According to the analyses discussed in Section 3[link], the angular distribution of radiation is uneven when the SR energy is less than 48.8 keV (χ1 > 3); the angular distribution of radiation is very uniform from 147 keV to 151 keV (χ1 < 1 and χ2 < 1); and the angular distribution of radiation is in a transition state that can be used in many experiments from 48.8 keV to 147 keV (1 < χ1 < 3) and from 151 keV to 300 keV (1 < χ2 < 1.43) (see Fig. 21[link]). Field measurements were carried out under a series of gaps. The data at the smallest gap (12.1 mm) were used to simulate the radiation at several SR energies (Figs. 22[link], 23[link], 24[link] and 25[link]). The dark stripe in the centre of the FOV is caused by the longitudinal positioning deviation between the horizontal magnet rows and the vertical magnet rows.

[Figure 21]
Figure 21
The designed parameters of the HEPS Mango wiggler and the parameters of the storage ring were substituted into equations (13)[link] and (17)[link]. (a) χ1 with an energy dependence. (b) χ2 with an energy dependence.
[Figure 22]
Figure 22
Simulated angular distribution of the power density of the HEPS Mango wiggler based on the magnetic field measured at the smallest gap (12.1 mm).
[Figure 23]
Figure 23
Flux density spectrum (0.8–300 keV, along the beam axis) of the HEPS Mango wiggler, simulated based on data measured at the smallest gap (12.1 mm).
[Figure 24]
Figure 24
Photon flux spectrum (0.8–300 keV, within 1.5 mrad) of the HEPS Mango wiggler, simulated based on data measured at the smallest gap (12.1 mm).
[Figure 25]
Figure 25
Angular distribution of the flux density for the HEPS Mango wiggler, simulated based on data measured at the smallest gap (12.1 mm). (a) SR energy at 48.8 keV. (b) SR energy at 147 keV. (c) SR energy at 151 keV. (d) SR energy at 300 keV.

5. Summary

A novel insertion device referred to as the Mango wiggler is introduced. It is constructed from two orthogonal planar wigglers with a small difference in their period lengths, eliciting the phase difference of the magnetic fields to incrementally transitions from 0 to π/2. The Mango wiggler can provide a large and symmetric FOV with a high flux and a high critical energy by expanding its vertical divergence to be the same as the horizontal divergence. As a result, the efficiency of SR X-ray in-line imaging experiments is improved while ensuring good imaging resolution. To study its radiation properties, detailed theoretical analyses and simulations were carried out via SPECTRA and Mathematica. The uneven angular distribution of radiation occurs at both low and high photon energies for different reasons. Mainly, the discretization of the spectrum results in the unevenness of the angular distribution of radiation at low energy. A tapered magnetic field can be applied to make the angular distribution of radiation at low energy more uniform. At high energy, a decrease in the SR divergence of the electron bunch results in longitude lines that cannot be uniformly filled. These theoretical analyses were then applied to the HEPS Mango wiggler and several simulations were conducted based on the measured magnetic field. Although the simulation results cannot provide final conclusions as the magnetic field requires further adjustment before the installation of the HEPS Mango wiggler, presenting the idea of the Mango wiggler is still meaningful.

Acknowledgements

We acknowledge the support in the processing, manufacturing, and field measurements of the Mango wiggler at HEPS from the optics group and the insertion device group.

Funding information

The authors acknowledge the support from the High Energy Photon Source (HEPS), a major national science and technology infrastructure in China (Grant No. (2017)2173).

References

First citationAlferov, D. F., Bashmakov, Y. A. & Besonov, E. G. (1974). Sov. Phys. Tech. Phys. 18, 1336.  Google Scholar
First citationBoyanov, B. I., Bunker, G., Lee, J. M. & Morrison, T. I. (1994). Nucl. Instrum. Methods Phys. Res. A, 339, 596–603.  CrossRef CAS Web of Science Google Scholar
First citationKim, K. J. (1989). AIP Conf. Proc. 184, 565–632.  CrossRef Google Scholar
First citationKroll, N., Morton, P. & Rosenbluth, M. (1981). IEEE J. Quantum Electron. 17, 1436–1468.  CrossRef Web of Science Google Scholar
First citationLandhuis, E. (2017). Nature, 541, 559–561.  CrossRef CAS PubMed Google Scholar
First citationOrzechowski, T. J., Anderson, B. R., Clark, J. C., Fawley, W. M., Paul, A. C., Prosnitz, D., Scharlemann, E. T., Yarema, S. M., Hopkins, D. B., Sessler, A. M. & Wurtele, J. S. (1986). Phys. Rev. Lett. 57, 2172–2175.  CrossRef PubMed CAS Web of Science Google Scholar
First citationRayleigh, J. W. S. (1879). London, Edinb. Dubl. Philos. Mag. J. Sci. 8, 261–274.  Google Scholar
First citationSnigirev, A., Snigireva, I., Kohn, V., Kuznetsov, S. & Schelokov, I. (1995). Rev. Sci. Instrum. 66, 5486–5492.  CrossRef CAS Web of Science Google Scholar
First citationTanaka, T. (2021). J. Synchrotron Rad. 28, 1267–1272.  Web of Science CrossRef IUCr Journals Google Scholar
First citationThompson, A. C. & Vaughan, D. (2001). X-ray Data Booklet. Lawrence Berkeley National Laboratory, University of California Berkeley, CA, USA.  Google Scholar
First citationWalker, R. P. (2003). Undulators, Wigglers and their Applications, edited by H. Onuki & P. Elleaume, pp. 108–147. London: Taylor & Francis.  Google Scholar
First citationWalker, R. P. (2017). Phys. Rev. Accel. Beams, 20, 020703.  CrossRef Google Scholar
First citationWolfram Research (2018). Mathematica, Version 11.3. Wolfram Research, Champaign, IL, USA.  Google Scholar
First citationXiaoming, J., Jiuqing, W., Qing, Q., Yuhui, D., Weifan, S., Jian, C., Gang, X., Tiandou, H., Hu, D., Fusan, C., Fengli, L., Huihua, L., Junhui, Y., Chunhua, L., Yi, S., Jinhui, C., Haiyi, D., Quan, C., Wei, X., Ming, L., Guangcai, C., Hongwei, Z., Ye, T., Peng, L., Jing, L., Dongbai, S., Gang, L., Hong, S. & Jianshe, C. (2014). Synchrotron Radiat. News, 27(6), 27–31.  CrossRef Google Scholar
First citationXing, L., McKellar, R. C., Xu, X., Li, G., Bai, M., Persons, W. S., IV, Miyashita, T., Benton, M. J., Zhang, J., Wolfe, A. P., Yi, Q., Tseng, K., Ran, H. & Currie, P. J. (2016). Curr. Biol. 26, 3352–3360.  CrossRef CAS PubMed 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