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

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775

Reconstruction of electron beam distribution using phase-retrieval algorithm

CROSSMARK_Color_square_no_text.svg

aDepartment of Physics, Pohang University of Science and Technology, Pohang 790-784, Republic of Korea, and bPohang Accelerator Labotratory, Pohang University of Science and Technology, Pohang 790-784, Republic of Korea
*Correspondence e-mail: chbkim@postech.ac.kr

(Received 12 May 2010; accepted 1 February 2011; online 11 March 2011)

The shapes of light sources such as electron beams can be reconstructed by inverse Fourier transformation of the complex degree of spatial coherency, which can be measured using Young's interferometer. The application of the phase-retrieval algorithm to reduce phase measurement errors in the complex degree of spatial coherency is numerically studied using an electron beam with an asymmetric distribution. This application is demonstrated with experimental data measured at the diagnostic beamline at the Pohang Accelerator Laboratory.

1. Introduction

Over the past several decades interferometer techniques have been successfully implemented to measure the size of an electron beam at various storage rings (Mitsuhashi, 1999[Mitsuhashi, T. (1999). Proceedings of the Joint US-CERN-Japan-Russia School on Particle Accelerator Beam Measurements, pp. 399-427. Singapore: World Scientific.]; Sakai et al., 2000[Sakai, I., Yamamoto, Y., Mitsuhashi, T., Amano, D. & Iwasaki, H. (2000). Rev. Sci. Instrum. 71, 1264-1267.]; Fisher et al., 2001[Fisher, A. S., Bong, E. L., Holtzapple, R. L. & Petree, M. (2001). Proceedings of the 2001 Particle Accelerator Conference, Chicago, USA, p. 547.]; Masaki & Takano, 2003[Masaki, M. & Takano, S. (2003). J. Synchrotron Rad. 10, 295-302.]). Recently, numerical and experimental studies have improved the measurement techniques (Naito & Mitsuhashi, 2006[Naito, T. & Mitsuhashi, T. (2006). Phys. Rev. ST Accel. Beams, 9, 122802.]; Parc et al., 2009[Parc, Y. W., Kim, C., Huang, J. Y. & Ko, I. S. (2009). J. Synchrotron Rad. 16, 642-646.]). The shape of an electron beam can be reconstructed by using the inverse Fourier transformation (IFT) of the complex degree of coherency, which can be measured by experiment. To determine the complex degree of coherency we need to measure the visibility and the phase of the interferogram (Mitsuhashi, 1999[Mitsuhashi, T. (1999). Proceedings of the Joint US-CERN-Japan-Russia School on Particle Accelerator Beam Measurements, pp. 399-427. Singapore: World Scientific.]). In general, the IFT result for the shape of the electron beam is distorted owing to the measurement error of the phase. Mitsuhashi (1999[Mitsuhashi, T. (1999). Proceedings of the Joint US-CERN-Japan-Russia School on Particle Accelerator Beam Measurements, pp. 399-427. Singapore: World Scientific.]) used the cosine Fourier transformation method to eliminate the phase measurement error. However, the application of this method is limited to symmetric distribution as in the case of the Gaussian shape of the electron beam (Mitsuhashi, 1999[Mitsuhashi, T. (1999). Proceedings of the Joint US-CERN-Japan-Russia School on Particle Accelerator Beam Measurements, pp. 399-427. Singapore: World Scientific.]).

Alternatively, the size of an electron beam with a Gaussian distribution can be determined by comparing the measured spatial coherency and the numerically calculated coherency (Sakai et al., 2000[Sakai, I., Yamamoto, Y., Mitsuhashi, T., Amano, D. & Iwasaki, H. (2000). Rev. Sci. Instrum. 71, 1264-1267.]; Parc et al., 2009[Parc, Y. W., Kim, C., Huang, J. Y. & Ko, I. S. (2009). J. Synchrotron Rad. 16, 642-646.]). However, this method requires several trials to match the numerically calculated result with the measured result under the assumption of a Gaussian distribution for the electron beam. If the electron beam shape is not Gaussian, it is impossible to estimate the source shape by using this method. Therefore, it is highly desirable to find a method that eliminates the phase error in reconstructing the non-symmetric distribution of the electron beam.

In this paper we show that the phase-retrieval algorithm introduced by Fienup (1980[Fienup, J. R. (1980). Opt. Eng. 19, 297-305.]) is appropriate for reconstructing the distribution of an electron beam that is asymmetric. The phase information required to reconstruct an image from the measured diffractogram can be retrieved by using the algorithm. To demonstrate the application of this method in the electron beam diagnostics, the complex degree of spatial coherence of light generated from the Pohang Light Source (PLS) is investigated using the diagnostics beamline at the Pohang Accelerator Laboratory. The theory of interference and spatial coherency is reviewed in §2[link]. A numerical study on the phase-retrieval algorithm and reconstructed image of the electron beam is described in §3[link]. The experimental application of the phase-retrieval algorithm is presented in §4[link]. Finally, a summary is provided in §5[link].

2. Theory

Using Young's interferometer the degree of correlation between two points in a wavefront generated from a light source can be investigated, as shown in Fig. 1[link]. The degree of correlation is called the spatial coherency of light. The interferogram constructed by an interferometer with two slits is calculated as (Goodman, 1985[Goodman, J. W. (1985). Statistical Optics, ch. 5. New York: John Wiley and Sons.])

[I\left(P\right)={\left\langle{{{\left|{{K_1}\tilde u\left({S_1},t - {{{r_1}} \over c}\right) + {K_2}\tilde u\left({S_2},t - {{{r_2}} \over c}\right)} \right|}^2}}\right\rangle_t},\eqno(1)]

where P is a point on the screen in Fig. 1[link], Ki is the amplitude of the wave at the slit Si for i = 1 and 2, and ri is the distance from the slit Si to the point P. If K1 and K2 are given by the same value K, equation (1)[link] can be written as

[I\left(P\right)= 2{K^2}{I_0}\left[{1 + \left|{\tilde \gamma \left({{{{r_2} - {r_1}} \over c}} \right)} \right|\cos \left({2\pi \bar \nu {{{r_2} - {r_1}} \over c} + \alpha } \right)}\right],\eqno(2)]

where I0 is the initial intensity of light as given by [|{\tilde u(t)}|^2], [\tilde \gamma] is the complex degree of spatial coherency, [\bar \nu] is the central frequency of light, c is the speed of light in the air and α is the phase of the interferogram. The visibility of the interferogram is defined as (Born & Wolf, 1999[Born, M. & Wolf, E. (1999). Principles of Optics, 7th ed. (extended), ch. 10. New York: Cambridge University Press.])

[V = {{{I_{\max }} - {I_{\min }}} \over {{I_{\max }} + {I_{\min }}}} = \left| {\tilde \gamma \left({{{{r_2} - {r_1}} \over c}} \right)} \right|.\eqno(3)]

The spatial coherency of light from a finite-size source can be estimated from the van Cittert–Zernike theorem (Cittert, 1934[Cittert, P. H. van (1934). Physica, 1, 201.]; Zernike, 1938[Zernike, F. (1938). Physica, 5, 785.]). When d is the distance between the slits, the van Cittert–Zernike theorem can be expressed as

[\tilde \gamma(d)=\int\limits_\sigma{I(\xi)\exp\left(-ik{d \over L}\xi\right)\,{\rm{d}}\xi },\eqno(4)]

where the wavenumber k = 2π/λ, ξ is the position of the infinitesimal light source, and L is the distance between the source and the slit. From (4)[link], the complex degree of spatial coherency is revealed as the Fourier transform (FT) of the intensity distribution I of the light source (Sakai et al., 2000[Sakai, I., Yamamoto, Y., Mitsuhashi, T., Amano, D. & Iwasaki, H. (2000). Rev. Sci. Instrum. 71, 1264-1267.]; Parc et al., 2009[Parc, Y. W., Kim, C., Huang, J. Y. & Ko, I. S. (2009). J. Synchrotron Rad. 16, 642-646.]).

[Figure 1]
Figure 1
Configuration of interference with finite-size light source. d is the distance between two slits. Si is the slit position and ri is the distance between points Si and P on the screen for i = 1, 2. L is the distance between the source and the slit mask.

3. Numerical study

In our numerical model all infinitesimal fragments of the light source emit light independently (Parc et al., 2009[Parc, Y. W., Kim, C., Huang, J. Y. & Ko, I. S. (2009). J. Synchrotron Rad. 16, 642-646.]). For the numerical study of the interferogram the intensity of the interferogram on the screen can be calculated by taking into consideration all of the phase changes owing to the path differences from an infinitesimal fragment of the light source to a point P on the screen.

A numerically calculated interferogram with light (wavelength 650 nm) generated from a Gaussian distribution source is shown in Fig. 2[link]. The points in Fig. 2[link] represent the numerical calculation results of the interferogram with a Gaussian light source. The inset of Fig. 2[link] shows the source distribution with a Gaussian shape, which is used in this numerical study. This Gaussian distribution will be used several times in this paper.

[Figure 2]
Figure 2
The numerically calculated result of the interferogram is shown by dots. The solid line is the plot of the theory in equation (5)[link]. The distance between the slits d is 1 cm and the half slit width a is 300 µm. The electron beam distribution used in this numerical study is shown in the inset. The standard deviation of the Gaussian distribution σ is 210 µm.

The formula for the interferogram from the Gaussian shape source with two slits is given by the relation

[I(P) = 2{\left [{{{\sin (u)} \over u}} \right]^2}\left [{1 + \exp \left( { - {{{\sigma ^2}} \over 2}{v^2}} \right)\cos \left({\delta + \alpha } \right)} \right], \eqno(5)]

with u = kax/R, v = kd/L and [\delta] = kdx/R, where the wavenumber k = 2π/λ, a is the half-width of the slit, R is the distance between the center (O) of the slits and a point (P) on the screen, x is the distance of the observation point from the screen center, σ is the standard deviation of the Gaussian source, L is the distance between the source and the slit mask, and d is the distance between the two slits.

The solid line in Fig. 2[link] is the plot of equation (5)[link] with the same parameters used in the simulation. The line in Fig. 2[link] is well matched with the interferogram from the numerical model. Thus, the theoretical and numerical models are in good agreement. In this numerical study the distance L between the source and the slit is 24.5 m, the slit distance d is set to 1.0 cm, the half-slit width a is 300 µm, and the distance D between the slit and the screen is 12 km (Parc et al., 2009[Parc, Y. W., Kim, C., Huang, J. Y. & Ko, I. S. (2009). J. Synchrotron Rad. 16, 642-646.]).

3.1. Amplitude and phase of FT

In this section the sum of two Gaussians is used to represent an asymmetric distribution of an electron beam, as shown in Fig. 3[link]. The standard deviations of the two Gaussian distributions are 210 µm (σ1) and 105 µm (σ2) in Fig. 3[link]. The separation of the two Gaussians (l) is 420 µm. Visibilities obtained from the numerical study are represented by circles in Fig. 4(a)[link]. The visibility is no longer Gaussian as shown in Fig. 4(a)[link].

[Figure 3]
Figure 3
Electron beam distribution used in this numerical study. The standard deviation of one Gaussian distribution is 210 µm and that of the other is 105 µm. The distance between the centers of the two Gaussian distributions is 420 µm.
[Figure 4]
Figure 4
(a) Visibility and amplitude of the FT of the distribution shown in Fig. 3[link]. The circles are the numerically calculated visibilities. The dots are the amplitudes of the FT. (b) Phase of the FT. The dots are the original phases obtained from the FT results with the distribution in Fig. 3[link]. The symbol × (+) is used for random phases with a maximum value of 0.4 (0.8) rad.

From equations (3)[link] and (4)[link] we can see that the amplitude of the FT of the electron beam distribution should be the same with this visibility. The amplitude of the FT of the same electron beam distribution used in the numerical study is plotted using dots in Fig. 4(a)[link] to compare the numerically calculated visibility. The numerically calculated visibility (circles) and the amplitude of the FT (dots) are closely matched.

The phases used in this section come from the FT of the electron beam distribution in the previous section. In the real experiment these will be measured from the interferograms. There will be some errors in the measurement result. These errors will provide distorted measurement results for a certain distribution of electron beams in a storage ring (Mitsuhashi, 1999[Mitsuhashi, T. (1999). Proceedings of the Joint US-CERN-Japan-Russia School on Particle Accelerator Beam Measurements, pp. 399-427. Singapore: World Scientific.]). In Fig. 4(b)[link] the original phases coming from the FT of Fig. 3[link] are shown by dots. To study the effects of the errors in the reconstruction of the electron beam, we assumed random errors in the original phase as shown in Fig. 4(b)[link]. Random errors within 0.4 (0.8) rad are added to the original phase, represented by the symbol × (+).

3.2. The effect on IFT owing to the phase error

The results of IFT with the phase including random errors are shown in Fig. 5[link]. The solid line is the result of IFT with original phases. The dash-dotted (dashed) line shows the IFT result with a random phase with a maximum error 0.4 (0.8) rad. The solid line is fitted by the two Gaussians; the standard deviation of one Gaussian is 210.1 µm and that of the other Gaussian is 105.1 µm. The separation of two Gaussians is 420.1 µm. The electron beam shape with larger phase errors is distorted more than that with smaller phase errors. As expected, the asymmetric shape is reconstructed when there is no phase error. However, the shape of the electron beam is distorted when there are random phase errors. Therefore, we need to find a method to reduce the phase error in this IFT process.

[Figure 5]
Figure 5
Electron beam shape obtained by IFT for the two Gaussians in Fig. 3[link]. The solid line represents the result of IFT with the original phase. The dash-dotted (dashed) line represents the IFT result with random phase of a maximum value of 0.4 (0.8) rad.

3.3. Phase-retrieval algorithm

The algorithm of the phase-retrieval method is shown in Fig. 6[link]. The initial starting function should be the result of the IFT shown in Fig. 5[link]. The function obtained by the IFT in Fig. 5[link] in the previous section is represented by g in Fig. 6[link]. We first make the FT of g, and the function obtained in the Fourier space is denoted by G in Fig. 6[link]. The amplitude of G is represented by |G|. This amplitude is replaced by the visibility V obtained in Fig. 4(a)[link] because this visibility will be measured in the experiment. The vindication of this replacement comes from the comparison between the visibility and the amplitude of the FT, as shown in Fig. 4(a)[link]. This replacement plays the key role in this algorithm, and this process is called the Fourier space constraint.

[Figure 6]
Figure 6
Scheme of the phase-retrieval algorithm.

We need one more strategy to complete the algorithm. The new function obtained after the replacement is denoted by [G^{\,\prime}] in Fig. 6[link]. By IFT we can obtain a function in real space denoted by [g^{\,\prime}] in Fig. 6[link]. This function usually has negative points. However, the electron beam distribution should not have a negative value because a negative distribution has no physical meaning. Thus, the negative points should be removed from [g^{\,\prime}] changing the negative values to zero. This process is called the real space constraint.

After changing the negative values to zero, the initial function g is replaced by [g^{\,\prime}]. We can iterate this process until the iteration does not show a meaningful difference from the previous step. This iteration process is called the phase-retrieval algorithm (Fienup, 1980[Fienup, J. R. (1980). Opt. Eng. 19, 297-305.]). Note that the negative value can be replaced with an arbitrary value if there is any need to try such a strategy. In this study the simple method of changing the negative value to zero is sufficient to reconstruct the electron beam shape.

Fig. 7[link] shows the recovery process of the electron beam distribution according to the phase-retrieval algorithm. The dashed line represents the initial distorted function from Fig. 5[link] with a maximum phase error of 0.8 rad. The dash-dotted line is the recovered function after the first iteration. The solid line is the result after 2000 iterations, which is sufficient to recover the electron beam shape shown in Fig. 3[link]. The asymmetric function is recovered with the phase-retrieval algorithm. The solid line is also fitted by two Gaussians. The standard deviation of one Gaussian is 216.3 µm and the other Gaussian is 100.2 µm. The separation of the two Gaussians is 426.59 µm. Even in the case of the maximum error of 0.8 rad the recovered result is reasonably close to the electron beam shape in Fig. 3[link].

[Figure 7]
Figure 7
Electron beam shape obtained by iterations. The dashed line represents the initial function with random phase errors with a maximum value of 0.8 rad, as shown in Fig. 4(b)[link]. The dash-dotted line is the reconstructed function [g^{\,\prime}] after one iteration. The solid line is the result after 2000 iterations.

4. Experimental study

A schematic diagram of the interferometer at the 1B1 beamline of PLS is shown in Fig. 8[link]. The radiation generated from the electron beam in the storage ring propagates to the end of the beamline. Three optical mirrors reflect only visible light in the transport path of the radiation to the interferometer (Huang & Ko, 1998[Huang, J. Y. & Ko, I. S. (1998). J. Synchrotron Rad. 5, 642-644.]). Note that these mirrors can be the source of phase measurement errors of the complex degree of coherency (Mitsuhashi, 1999[Mitsuhashi, T. (1999). Proceedings of the Joint US-CERN-Japan-Russia School on Particle Accelerator Beam Measurements, pp. 399-427. Singapore: World Scientific.]). A commercial band-pass filter [CVI Melles Griot, 03 FIV 022 (F10-650.0-50.0M)] is used between the slit and the detector in the interferometer. The FWHM bandwidth of the optical band-pass filter is 10 nm. In the experiment, 650 nm is selected by the band-pass filter. A focusing lens is used to overlap two lights passing through each slit on the CCD plane.

[Figure 8]
Figure 8
Schematic diagram of 1B1 beamline at PLS.

4.1. Measurement of visibility and phase

A measured interferogram is shown in Fig. 9[link] by squares. In the measurement the slit width is 500 µm, the slit length is 6 cm and the distance d between two slits is 6 mm. The distance L between the source and the slit is 24.5 m. The fitting formula on the basis of equation (5)[link] is given by

[I = 2{\left [{{{A\sin (Bx)} \over x}} \right]^2}\left [{1 + E\cos \left({Cx + D} \right)} \right], \eqno(6)]

where x is the pixel position. The dotted line in Fig. 9[link] is the fitting line with equation (6)[link]. The constants are given for the fitting line as A = 2.281, B = 0.4145, C = 7.995, D = 1.817 and E = 0.7943. D is the experimental phase of the complex degree of coherency. In this fitting result D, a fitting error and measurement error exist. Our objective is to reduce this error using the phase-retrieval algorithm.

[Figure 9]
Figure 9
Measured interferogram (squares) at PLS and the fitting result (lines) when the distance d between two slits is 6 mm. The constants for this fitting after equation (6)[link] are obtained as A = 2.281, B = 0.4145, C = 7.995, D = 1.817 and E = 0.7943.

The visibilities are measured by 25 slit masks with different slit distances d and the results are represented by circles in Fig. 10(a)[link]. Note that these visibilities are not obtained from the fitting results but from the maximum and minimum values from the measured interferograms. The dots in Fig. 10(a)[link] are the amplitudes of FT results with a single Gaussian distribution. Two results are different, which means that the electron beam distribution may be not Gaussian at PLS. The measured phase is shown by dots in Fig. 10(b)[link].

[Figure 10]
Figure 10
(a) The circles are measured visibilities. The dots are the amplitudes of the FT of a Gaussian distribution with σ = 210 µm. (b) Measured phase.

4.2. Reconstruction of electron beam shape by IFT and the phase-retrieval algorithm

The IFT result in Fig. 10[link] is shown by a dashed line in Fig. 11[link]. As expected, the result is distorted (Mitsuhashi, 1999[Mitsuhashi, T. (1999). Proceedings of the Joint US-CERN-Japan-Russia School on Particle Accelerator Beam Measurements, pp. 399-427. Singapore: World Scientific.]). The recovery of the electron beam is attempted by the phase-retrieval algorithm and the result is shown in Fig. 11[link]. The dash-dotted line is the first iteration result. The result after 100 iterations is plotted with a solid line. The distribution of the electron beam at PLS is revealed as not a Gaussian distribution. This fitting result tells us that the half width at half-maximum is 317 µm. The right-hand side of the recovered electron beam distribution in Fig. 11[link] can be fitted with a Gaussian function. The standard deviation of the fitted result is 273 µm.

[Figure 11]
Figure 11
Results obtained by IFT with measured visibility and phase. The dashed line represents the initial function obtained by IFT. The dash-dotted line is the reconstructed function after first iteration with the phase-retrieval algorithm. The solid line is the result after 100 iterations.

5. Summary

A theory of interference is reviewed in order to understand the relation between the source shape and the spatial coherency of light. A random phase error is added in the phase part of the FT of the given electron distribution to simulate the experimental situation. A distorted distribution is obtained by the IFT with a random phase error. The phase-retrieval algorithm is applied to recover the electron distribution used in the numerical study from the distorted IFT result. The phases are obtained from the fitting result of the measured interferograms. The phase errors from the fitting and experimental measure are decreased by the phase-retrieval algorithm, and the results show us that the electron distribution in PLS is not a Gaussian distribution. This technique can be applied to measure the distribution of the electron beam in fourth-generation light sources such as XFEL. This method is applicable to the reconstruction of any kinds of light sources such as stars.

Acknowledgements

This work was supported by a grant from the National Research Foundation of Korea, funded by the Korean government (Ministry of Education, Science and Technology, Korea) (grant No. 2008-0059842). This research was partly supported by the Basic Science Research Program of the National Research Foundation of Korea, funded by the Ministry of Education, Science and Technology, Korea (grant No. 2010-0008034). We would like to extend our gratitude to Dr Beom Hyun Kim (Department of Physics, Pohang University of Science and Technology) for his assistance and discussions.

References

First citationBorn, M. & Wolf, E. (1999). Principles of Optics, 7th ed. (extended), ch. 10. New York: Cambridge University Press.  Google Scholar
First citationCittert, P. H. van (1934). Physica, 1, 201.  Google Scholar
First citationFienup, J. R. (1980). Opt. Eng. 19, 297–305.  CrossRef Google Scholar
First citationFisher, A. S., Bong, E. L., Holtzapple, R. L. & Petree, M. (2001). Proceedings of the 2001 Particle Accelerator Conference, Chicago, USA, p. 547.  Google Scholar
First citationGoodman, J. W. (1985). Statistical Optics, ch. 5. New York: John Wiley and Sons.  Google Scholar
First citationHuang, J. Y. & Ko, I. S. (1998). J. Synchrotron Rad. 5, 642–644.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationMasaki, M. & Takano, S. (2003). J. Synchrotron Rad. 10, 295–302.  Web of Science CrossRef IUCr Journals Google Scholar
First citationMitsuhashi, T. (1999). Proceedings of the Joint US–CERN–Japan–Russia School on Particle Accelerator Beam Measurements, pp. 399–427. Singapore: World Scientific.  Google Scholar
First citationNaito, T. & Mitsuhashi, T. (2006). Phys. Rev. ST Accel. Beams, 9, 122802.  Web of Science CrossRef Google Scholar
First citationParc, Y. W., Kim, C., Huang, J. Y. & Ko, I. S. (2009). J. Synchrotron Rad. 16, 642–646.  Web of Science CrossRef IUCr Journals Google Scholar
First citationSakai, I., Yamamoto, Y., Mitsuhashi, T., Amano, D. & Iwasaki, H. (2000). Rev. Sci. Instrum. 71, 1264–1267.  Web of Science CrossRef CAS Google Scholar
First citationZernike, F. (1938). Physica, 5, 785.  CrossRef Google Scholar

© International Union of Crystallography. Prior permission is not required to reproduce short quotations, tables and figures from this article, provided the original authors and source are cited. For more information, click here.

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