Single shot coherence properties of the free-electron laser SACLA in the hard X-ray regime.

We measured the coherence properties of the free-electron laser SACLA on a single shot basis at an X-ray energy of 8 keV. By analysing small-angle X-ray scattering speckle patterns from colloidal dispersions we found a degree of transverse coherence of βt = 0.79 ± 0.09. Taking detector properties into account, we developed a simulation model in oder to determine the degree of coherence from intensity histograms. Finally we calculated a coherence time of τc = 0.1 fs and a pulse duration of 5.2 fs which corresponds with previous predictions.

The rise of X-ray free electron lasers (XFEL) allows the study of physical and chemical properties that could not be accessed experimentally before. In particular, the ultrashort pulse duration on the order of 10 fs at hard x-ray wavelengths of around 1.5 Å offers the possibility to access physical and chemical processes on molecular time and length scales by single shot x-ray diffraction patterns.

Many of such studies rely on the coherence properties of the X-ray beam. This includes, e.g., the study of dynamic properties via X-ray photon correlation spectroscopy123, imaging single particles and biological specimens by coherent diffractive imaging45, cross correlation techniques to reveal a single particle's shape6 and orientational order in liquids and glasses7, or the imaging of ultrafast magnetization processes8. XFEL sources such as the Linac Coherent Light Source (LCLS) at SLAC9 and SACLA at SPring-810 are based on self-amplified spontaneous emission (SASE)1112. The SASE process originates from random fluctuation in the electron beam. Therefore, the coherence properties of such an X-ray source differ from optical lasers. In particular, every single X-ray shot is expected to show certain random features because of the stochastic nature of the SASE process. This was observed recently in experiments at the LCLS and in simulations on the longitudinal coherence properties, where a significant shot-to-shot variation has been observed13. In contrast, a SASE beam is expected to show a superior transverse coherence12.

Recently, several studies have been performed on the different FEL sources' coherence properties. In the soft x-ray regime, pinhole and double slit experiments have been performed at FLASH1415 and LCLS16 in order to study the coherence lengths and times. In further studies, interferometry methods have been used1718 to determine longitudinal and transverse coherence lengths. In contrast, coherence properties in the range of hard X-ray have been examined by analysing coherent diffraction patterns, so-called speckle patterns, from disordered samples on a shot-to-shot basis1319. A similar study has been performed at FLASH20. In this study, the degree of coherence is obtained from the speckle contrast. By analysing single shot diffraction patterns from colloidal particles in small angle X-ray scattering (SAXS) geometry, the degree of transverse coherence was found to be almost 11319. Furthermore, the longitudinal coherence was determined by wide angle scattering (WAXS) experiments on amorphous metals and it also appeared to be very sensitive on the SASE fluctuations319.

Here we present a study of the coherence properties of the free-electron laser SACLA in the hard X-ray regime. In contrast to the LCLS, three bunch compressors and short-period in-vacuum undulators are used at SACLA that result in a compact design compared to LCLS10. Furthermore, initial experiments reported pulse durations well below 10 fs21. The question arises if and how this design combined with the short pulse duration have an impact on the coherence properties of the single x-ray pulses. We calculated the degree of coherence for single shots from speckle patterns of colloidal particles at a hard X-ray energy of 8 keV. The number of transverse modes was determined by the need of correcting for the energy bandwidth of the x-ray pulses. The results are supported by simulations of the impact of the charge sharing process of CCD detectors on intensity histograms of speckle patterns. Finally, we estimate the pulse duration from the intensity fluctuations of the beam and compare it to the expected values in the order of 5 fs.

Results

A typical speckle pattern from an amorphous sample together with the azimuthally averaged intensity I(q) is shown in Fig. 1. The experimental details are given in the Methods section. The size of individual speckles depends on set-up parameters only and is given by s = λd/l, with wavelength λ, sample-detector distance d and beam size at the sample l. The experimental speckle size can be determined by calculating the spatial intensity correlation function g(R) = 〈I(r + R)I(r)〉/〈I(r)〉2, that is connected to the coherence properties of the beam via the mutual coherence function γ(R, τ) by g(R) = |γ(R, τ = 0)|2 + 1 at the spatial distance R and time difference τ22. γ(R, τ) is a complex function that contains all information of the coherence properties. The degree of coherence for a single speckle pattern is given by , with M denoting the mode number of the beam1923. By modelling g(R) with a Lorentzian function for both the vertical and horizontal directions, we obtain a beam size of (1.9 ± 0.1) × (1.6 ± 0.1) μm2 (h × v) which corresponds with the directly measured beam size of 1.8 × 1.5 μm2 excellently.

Figure 1

(a) Single shot speckle pattern and corresponding I(q) averaged over 500 shots (inset). (b) Typical intensity autocorrelation g(R) and Lorentzian fit to determine the speckle size. Due to the finite size of the extracted region of interest no full decorrelation (g(R) − 1 = 0) can be observed.

By analysing the speckle pattern the most direct quantity that can be accessed is the speckle contrast or visibility β. It is given by the normalized standard deviation of the intensity yielding222324 Due to the equality , the speckle contrast directly reflects the degree of coherence of the x-ray beam. A fully coherent beam would have a mode number of M = 1 and thus a contrast of β = 1. The experimentally measured contrast is affected by both the longitudinal and transverse coherence of the x-ray beam yielding β = β(q), with the transverse coherence β and a q-dependent correction β(q) that depends on the energy bandwidth and the scattering geometry192526.

In Fig. 2 (a) the contrast β is shown for 2000 single shots at q = 0.041 ± 0.001 nm−1. The average scattering intensity in this q-range was 35 photons/pixel. We find an averaged contrast of 〈β〉 = 0.70 ± 0.08 which can be taken from the histogram of the contrast values as shown in Fig. 2 (b). This corresponds to a mode number of M = 2.0 ± 0.5. The single shot contrast varies from slightly above 0.5 up to some fully coherent shots with contrast values of 1 following an asymmetric distribution. This is demonstrated by the solid line in Fig. 2 (b) reflecting a lognormal distribution. In addition, we did not find any significant variation as a function of q, see Fig. 2 (c). Taking into account the natural bandwidth of the FEL radiation of Δλ/λ = 5 × 10−3, the correction accounting for the effect of the longitudinal coherence can be calculated to1925 with beam size L, sample thickness W, and the coefficients , , and k0 = 2π/λ. For the q-range studied here, we find a correction factor of β = 0.89, i.e. the transverse coherence is found to be β = β/β = 0.79 ± 0.09. This corresponds to a transverse mode number of M = 1.6 ± 0.4.

Figure 2

(a) Single shot speckle contrast β from 2000 consecutive shots calculated at q = 0.041 nm−1. (b) Histogram of β. The solid line represents a lognormal distribution using a mean value of 0.7 and a width of 0.105. (c) Speckle contrast β as function of q averaged over 2000 consecutive single shots. The error bars represent the standard deviation of the single shot values from the averaged value.

Typically, a certain minimum intensity is necessary for contrast calculation using Eq. 1. At low scattering intensities the contrast is usually obtained by describing the intensity distribution by intensity probability distribution functions22023. This method was shown to sufficiently describe the contrast of low-intensity speckle patterns in particular at large q 3 and by using split-and-delay devices27. Due to charge sharing of neighbouring pixels in the CCD detector and read out noise a continuous histogram is measured with peaks at intensity values corresponding to multiples of single photon events rather than discrete photon events. Typically, experimental intensity values are discretized by rounding up or down to integer values. This may introduce some artifacts. As the ratio of single, double, triple, … photon events is crucial for obtaining the mode number13, the interpretation can have a strong impact on the resulting contrast. At low intensities, such problems can be avoided by using so-called droplet algorithms to identify the location of the photons impacting the pixels of the detector3. However, those algorithms are limited to very low intensities, in particular work best for non-overlapping single photon events28.

Therefore, here we use the experimental, non-discretized intensity histograms and model those by simulations of intensity patterns following a negative binomal distribution and take into account charge sharing effects of the CCD detector2930 and detector noise, see Methods for details. We used M = 2, which corresponds to the experimental value discussed above, as input parameter for the mode numbers. The results are shown as red line in Fig. 3. For comparison, the black circles represent the discrete values before blurring. In particular, at low intensities the results correspond to the experimental data. At higher intensities (see Fig. 3 (b)), the model describes the general shape of the experimentally obtained histogram very well, however, small differences can be found for and also at higher intensities. We think this is due to the simplicity of our model and due to a better statistic in the simulation model. Because of blurring, real scattering intensities can appear as lower intensities. For example, a single event of intensity I0 hitting the CCD sensor in between two pixels, will appear as two independent events of intensity I0/2. These can not be distinguished from other events of intensity I0/2 in the same pixels. As a result, experimentally obtained mode numbers from intensity histograms may be overestimated.

Figure 3

Intensity histograms at q = 0.071 ± 0.003 nm−1 (a) and q = 0.053 ± 0.003 nm−1 (b) averaged over 600 consecutive shots and simulation results taking blurring and noise into account.

Finally we want to discuss the impact of the shot-to-shot intensity distribution. Thanks to the energy bandwidth of 40 eV, we estimate the coherence time τ = ξ/c = 0.1 fs, with the speed of light c. The coherence time is proportional to the average x-ray pulse duration 〈t〉 via 〈t〉 = M19, with the average number of temporal modes M. We determined M from the shot-to-shot intensity distribution via a Gamma-distribution1123 This equation is valid at least for the linear FEL regime32. The pulse intensity was measured by a diode located in front of the optical components of the beamline for more than 10000 single shots. The shot-to-shot variation of the beam intensity is displayed in Fig. 4 (a). By fitting its distribution with Eq. 3 we obtain M = 52.2 ± 2.0, see Fig. 4 (b). This results in an average x-ray pulse length of t = 5.2 ± 0.2 fs which perfectly corresponds to other recent results from SACLA21.

Figure 4

(a) Shot to shot variation of the beam intensity. (b) Intensity distribution. The solid line is a fit of Eq. 3 yielding M = 52.2 ± 2, the mean intensity was 〈I〉 = 300 μJ.

Discussion

In conclusion, we present a study on the coherence properties of the hard x-ray free-electron laser SACLA based on the analysis of single shot speckle patterns from colloidal particles. The size of the speckles in horizontal and vertical directions are in perfect agreement with expected values. We found a degree of transverse coherence of 0.79 ± 0.09 corresponding to a mode number of M = 1.6. The degree of transverse coherence is affected by set-up and machine parameters such as bunch compression schemes, saturation length, and photon wavelength resulting in observed values below 1. This is in agreement with expectations from FEL theory1231 and previous experiments161920.

By using a simple simulation model, we were able to model experimental intensity histograms. Because of the model's simplicity, the illustration of the values in the histogram is vague, i.e. they are blurred. This blurring's origin is related to the charge sharing and the noise level which is present in CCD detectors used in modern x-ray scattering experiments. Such simulation methods will help to determine contrasts more efficiently, in particular for weak scattering intensities where time-consuming droplet algorithms are used to overcome the effects of charge sharing3. However, we want to note that the simple simulation model discussed here does not give a full description of such influences. Further knowledge of detector properties is needed for a full understanding of our observations.

Methods

Sample

We synthesized the silica colloidal particles by using a modified Stöber method3334. The particles were dispersed in water at a volume fraction of about 5%. They had a mean radius of 41 nm with a polydispersity of 10% which was determined by a form factor analysis (see following section). After the preparation the dispersion was filled into a quartz capillary (diameter of about 0.7 mm) which has been sealed by melting the glass.

Experiment

The experiment was performed at SACLA's EH3 of BL3 at a photon energy of 8 keV. With a bandwidth of 5 × 10−3 at a repetition rate of 20 Hz35, the pulse energy was (300 ± 40) μJ. The beam was focalized with KB mirrors to a size of 1.8 × 1.5 μm2 (h × v) at the sample position. The beam size was measured in the vertical and horizontal directions by knife edge scans. The sample capillaries were placed into a holder that was mounted in the MAXIC sample chamber of BL336. To resolve the structure factor of the samples, the detector (dual multiport charge coupled device (MPCCD)38, i.e. two sensors with 1024 × 512 pixels each with individual pixel size of 50 μm × 50 μm) was placed at a distance of 3 m behind the sample. The accessible q-range was limited by the beam stop and the size of the detector to 0.02 nm−1 < q < 0.34 nm−1. For each sample at least 10000 single shot speckle patterns were taken. These were averaged and azimuthally integrated to obtain I(q) as shown in Fig. 1. We measured several 1000 speckle patterns from a diluted dispersion (volume fraction ) in order to extract the particle form factor. This was fitted by a form factor of a sphere convoluted with a size distribution to account for the polydispersity39. At the volume fraction of 5% used for the coherence analysis a structure factor reflecting particle-particle correlations was visible39, see Fig. 1.

Charge sharing simulations

Depending on the detector used, the photons impinging on the CDD sensor generate a charge cloud of different sizes, typically in the order of 10 μm37. This may result in shortcomings, e.g. an intensity of 1.4 can either be achieved by one photon with negligible charge sharing and less charge sharing contribution from neighbouring pixels, or e.g. by two photons (or more) with non-negligible charge sharing contribution. The impact of charge sharing of the CCD pixels and noise was estimated by computer simulations. Therefore, we used several 500 × 500 pixel grids with values corresponding a negative binomal distribution at an average intensity of 0.36 and 1.54 photons per pixel (identical to the experimental data), respectively. The negative binomal distribution function p(I) is given by23 where I is the intensity in number of photons, 〈I〉 its average, and Γ the Gamma-function. The initial mode number was set to M = 2, matching the result discussed above. Each pixel pattern was blurred by a random Gaussian distribution to model the charge sharing of pixels given by Here, x, y denote the pixel coordinate and x, y the random center of the Gaussian. After that, Gaussian noise was added to the pattern. The amplitude of the charge sharing parameter σ and the noise were adapted to match the experimentally obtained histograms. Therefore, we calculated an averaged histogram from 600 single shots showing low intensities at q = 0.071 ± 0.003 nm−1, see Fig. 3 (a). The intensities in histograms are scaled for the ADU (Analog-Digital Unit) value for 1 photon of I(1 photon) = 124.7 adu. By comparing the simulated histogram with the averaged experimental one we find σ ≈ 0.15 for the strength of charge sharing and a width of Gaussian noise of σnoise ≈ 0.1 photons. The effect of charge sharing is only slightly above the expected value for the MPCCD detector of pixels38.

Author Contributions

F.L., C.G., W.R., M.Y., T.I. and G.G. designed research. T.K., K.T. and M.Y. prepared the experimental setup. B.F. synthesized and prepared the samples. F.L., C.G., B.F., M.A.S., M.S., S.S., J.G., M.C., S.N., K.T., T.K., A.R. and G.G. performed the experiment. F.L., C.G., W.R. and G.G. analyzed and modeled the data. F.L. and C.G. wrote the manuscript text. All authors reviewed and discussed the manuscript.