Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Thermal effects of beam profiles on X-ray photon correlation spectroscopy at megahertz X-ray free-electron lasers

Open Access Open Access

Abstract

X-ray free-electron lasers (XFELs) with megahertz repetition rates enable X-ray photon correlation spectroscopy (XPCS) studies of fast dynamics on microsecond and sub-microsecond time scales. Beam-induced sample heating is one of the central concerns in these studies, as the interval time is often insufficient for heat dissipation. Despite the great efforts devoted to this issue, few have evaluated the thermal effects of X-ray beam profiles. This work compares the effective dynamics of three common beam profiles using numerical methods. Results show that under the same fluence, the effective temperatures increase with the nonuniformity of the beam, such that the Gaussian beam profile yields a higher effective temperature than the donut-like and uniform profiles. Moreover, decreasing the beam sizes is found to reduce beam-induced thermal effects, in particular the effects of beam profiles.

© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

X-ray photon correlation spectroscopy (XPCS) is an emerging technique to probe nanoscale and microscale dynamics [14]. It has been applied to various sample systems including amorphous materials [58], colloidal dispersions [911], and strongly correlated materials [12,13]. XPCS is typically performed in the sequential mode, where the sample is continuously illuminated by a coherent X-ray beam to produce a series of speckle patterns at desired delays. The dynamic properties of the sample are then extracted from the patterns via correlation analysis and model fitting.

High repetition-rate X-ray free-electron lasers (XFELs) such as the European XFEL, LCLS-II, and SHINE can deliver XFEL pulses at megahertz (MHz) repetition rates, enabling XPCS studies on microsecond and sub-microsecond time scales. In such studies, the time interval is often insufficient for heat dissipation [14], such that beam-induced sample heating becomes one of the central concerns. Pioneer studies conducted by Lehmkühler et al., Dallari et al., and Reiser et al. have carefully evaluated the thermal effects of fluences and dose rates on MHz XPCS [1517]. Meanwhile, although the beam profiles of XFEL pulses can vary [1820], few have evaluated the thermal effects of beam profiles. It is possibly because of the difficulty in facilitating beams of different profiles and comparable fluences. However, as the thermal patterns of sample heating at MHz XFELs largely depend on beam profiles, they would presumably affect the effective sample dynamics.

In this work, we performed a numerical study to evaluate the thermal effects of beam profiles on XPCS at MHz XFELs, using parameters in a recent experiment performed at the European XFEL [15]. We first simulated the Brownian motion of particles dispersed in water, which were illuminated by MHz XFEL pulses of three common beam profiles. We then retrieved the dynamics via XPCS correlation analysis and model fitting. Results show that nonuniform beam profiles lead to faster effective dynamics, such that the Gaussian profile yields higher effective temperatures than the donut-like and uniform ones. Moreover, decreasing the beam sizes is found to reduce beam-induced thermal effects, in particular the effects of beam profiles.

2. Methods

2.1 Initial conditions

Brownian motion simulations were performed on spherical particles dispersed in water, following the recent experiment performed at the European XFEL [15]. Key parameters were almost identical to the ones therein, including the hydrodynamic radius (69 nm), the capillary diameter (0.7 mm), the capillary wall thickness (10 µm), the photon energy (9.3 keV), the maximum pulse number (120), the repetition rate (1.128 MHz), and the size of the Gaussian beam (4 µm). A transverse circular region of 8 µm diameter was taken as the probing region, which corresponds to 4.71σ acceptance of the Gaussian beam. Four fluences therein (H = 3.9, 10.5, 27.7, and 56.8 mJ/mm2) as well as H = 0 were simulated.

Initially, 400 particles were randomly distributed in the probing region (the dashed circle). Figure 1(a) shows the typical initial particle distribution, where the X- and Z-positions were separately generated using the numpy.random.random function of Python. X and Z are the two transverse directions perpendicular to the beam and orthogonal to each other. The environmental temperature was set to 293 K for all conditions. Three beam profiles were used, including the uniform profile (flat-top), the donut-like profile (TEM01*), and the Gaussian profile (TEM00).

 figure: Fig. 1.

Fig. 1. (a) The typical particle distribution in the transverse probing region of 8 µm diameter, as indicated by the red dashed circle. Every blue sphere denotes one particle. X and Z are the two transverse directions perpendicular to the beam. (b)(c)(d) The adiabatic temperature increases of water induced by each uniform, donut-like, and Gaussian pulse. The fluence H is 10.5 mJ/mm2. The color bar indicates the temperature increase in K.

Download Full Size | PDF

2.2 Modeling beam-induced temperature increases

In this work, photon absorption by particles is ignored. Week-absorbing materials like water mainly dissipate the heat in the radial directions (transversely from the beam center outward). We thus integrated the heat in the beam direction (Y) and obtained the temperature increase on the XZ plane.

The effective fluences Heff at the sample depend on the original fluences H and the transmission of the beam through quartz and water, yielding

$${H_{\textrm{eff}}} = H\textrm{ }Ts\textrm{ }[1 - \exp ( - W/L)], $$
where Ts = 0.95 is the transmission of the front capillary wall, W is the capillary diameter, and L is the attenuation length of water at 9.3 keV.

The adiabatic temperature increase on the XZ plane induced by each uniform pulse is shown in Fig. 1(b). It was calculated with

$$\Delta {T_U} = \frac{{{H_{\textrm{eff}}}}}{{4W\rho {c_p}}},$$
where the area term gets canceled, ρ is the water density, and cp is the specific heat capacity of water. The reduction factor 4 results from the fact that the probing region diameter is twice the beam size in the Ref. [15]. It ensures that the applied pulse energy is the same as the one of the same H therein.

The adiabatic temperature increase on the XZ plane induced by each donut-like pulse is shown in Fig. 1(c). It was calculated with

$$\Delta {T_D}(r) = \frac{{{H_{\textrm{eff}}}}}{{4W\rho {c_p}}}N{F_D} \exp [ - \frac{{{{(r - \mu )}^2}}}{{2{\sigma ^2}}}], $$
where NFD = 1.70 is the normalization factor of the pulse energy for donut-like pulses, $r = \sqrt {{X^2} + {Z^2}}$ is the radial distance from the beam center, µ = d/4 is the offset, $\sigma = {\raise0.7ex\hbox{$d$} \!\mathord{\left/ {\vphantom {d {2\sqrt {2\ln 2} }}}\right.\kern-\nulldelimiterspace}\!\lower0.7ex\hbox{${2\sqrt {2\ln 2} }$}}$ is the standard deviation of the Gaussian beam size d. After normalization, the pulse energies of different profiles are the same at the same H.

The adiabatic temperature increase on the XZ plane induced by each Gaussian pulse is shown in Fig. 1(d). It was calculated with

$$\Delta {T_G}(r) = \frac{{{H_{\textrm{eff}}}}}{{4W\rho {c_p}}}N{F_G}\exp ( - \frac{{{r^2}}}{{2{\sigma ^2}}}), $$
where NFG= 2.96 is the normalization factor of the pulse energy for Gaussian pulses. The resulting peak temperature agrees well with the literature value [15].

The temperature increase after N pulses was obtained by multiplying the right sides of Eqs. (2) (3) and (4) by an accumulation factor $\sum\nolimits_{n = 0}^N {\{ 1 - \exp [\frac{{ - 1}}{{2(1 + n\Delta t/{t_0})}}]} \}$, where $\Delta t = \; 0.886\;\mathrm{\mu s}$ is the time interval and t0 is the thermal relaxation time. The t0 for Gaussian beams was adopted from the literature [15], whereas the ones for non-Gaussian beams were calculated according to the proportions of deposited pulse energies as shown in the first section of the Supplement.

2.3 Brownian motion simulation

Brownian motion simulations were performed in 2D on the XZ plane. The particles motion was driven by diffusion with

$$D = \frac{{{k_b}T}}{{6\pi \eta R}}, $$
where D is the diffusion coefficient, kb is the Boltzmann constant, T is the solution temperature, and η is the T-dependent dynamic viscosity [21]. After each time interval Δt, every particle moved by two random normal steps, one in the X-direction and one in the Z-direction. The step sizes were generated independently for each particle and each direction with a standard deviation of $\sqrt {2D\Delta t} $ and a mean value of 0, using the numpy.random.randn function of Python. If some of the particles were to escape the probing region, the moving directions were reversed. As the particle movement was minimal compared with the size of the probing region, these events were rather rare, and the dynamics were almost unchanged by reversals.

2.4 Correlation analysis and model fitting

In this section, intermediate scattering functions (ISFs) were obtained by autocorrelating the scattering amplitudes, and the dynamics parameters (the effective temperature Teff and the compression factor γ) were then obtained by fitting the ISFs.

First, the scattering amplitudes were calculated with

$$E(q) = \sum\nolimits_{i = 1}^{Np} {A({r_i})\exp ( - iq{r_i})}, $$
where the form factor of particles is ignored, Np = 400 is the total particle number in the probing region, q is the scattering vector, and ri is the position of the ith particle. A(ri) is the position-dependent amplitude factor, which equals unity for uniform pulses, $\sqrt {\exp [ - \frac{{{{({r_i} - \mu )}^2}}}{{2{\sigma ^2}}}]} $ for donut-like pulses, and $\sqrt {\exp ( - \frac{{{r_i}^2}}{{2{\sigma ^2}}})} $ for Gaussian pulses.

Then, the ISFs were extracted from E(q) via [22]

$$f(q,\tau ) = \frac{{ < E(q,t)E(q,t + \tau ) > }}{{ < E(q,t){ > ^2}}}. $$
For each condition, 30’000 simulations and autocorrelations were performed, yielding 30’000 ISFs of rather high statistic errors. They were first grouped into 3 datasets (3 × 10’000). Then, the ISFs of each dataset were averaged to yield 3 statistically reliable ISFs. They were all fitted by the Kohlrausch–Williams–Watts (KWW) function
$$f(q,\tau ) = \exp ( - |\Gamma (q)\tau {|^\gamma }), $$
where Γ(q) is the relaxation rate and γ is the compression factor. For diffusion, Γ(q) depends on the effective temperature Teff, with
$$\Gamma (q) = \frac{{{k_b}{T_{\textrm{eff}}}}}{{6\pi \eta R}}{q^2}. $$
 Finally, Teff and γ were obtained by averaging the fitting results of 3 ISFs. All Teff were normalized, such that the uniform Teff of H = 0 equals 293 K.

3. Results and discussion

3.1 Intermediate scattering functions and fits

Figure 2 shows the simulated intermediate scattering functions (symbols) and fitting curves (solid lines) of all three beam profiles at different fluences H. Firstly, at H = 0, the ISFs of different beam profiles almost overlap. Despite different statistic weights applied by nonuniform profiles to particles at different positions, no differences are observable when no heating takes place. It validates our simulation protocol. Secondly, at H > 0, ISFs of different profiles (colored symbols and lines) separate from each other. At the same H, the ISFs of donut-like beams decay faster than the ones of uniform beams and slower than the ones of Gaussian beams. It confirms our presumption about the thermal effects of beam profiles. As the heat accumulates, nonuniform beams apply higher statistic weights to particles in the hotter regions (see Fig. 1(b) and Fig. 1(c)). The higher the nonuniformity is, the faster the decays of the ISFs become. Moreover, the ISFs of H = 27.7 mJ/mm2 decay slightly slower than the ones in the reference because of the neglect of photon absorption by particles [15]. Below, we will further analyze the results by inspecting the effective temperatures and the compression factors obtained from model fitting.

 figure: Fig. 2.

Fig. 2. Simulated intermediate scattering functions (symbols) and their fitting curves (solid lines) at q = 0.125 nm-1. The fluences and beam profiles are shown in the legend.

Download Full Size | PDF

Figure 3(a) shows the effective temperatures Teff as a function of the fluences H. At H = 0, Teff of different beam profiles are almost the same, in line with the overlapping ISFs above. At H > 0, the Gaussian profile always yields the highest Teff, whereas the uniform profile always yields the lowest Teff, again indicating the significant effects of beam nonuniformity on the effective dynamics. Moreover, all curves are highly linear, suggesting that increasing H only enhances these effects. As expected, our simulations result in lower Teff than the experimental ones because of the neglect of photon absorption by particles [15]. Figure 3(b) shows the compression factors γ plotted as a function of H. At H = 0, all γ are close to unity. At H > 0, γ decrease with the increase of H because of the nonequilibrium induced by continuous heating and averaging over pulses [15], which is confirmed by a pulse-dependent analysis shown in Fig. S1 of the Supplement. More importantly, the decrease of γ is more profound for the nonuniform profiles, again suggesting strong effects of beam profiles on the effective dynamics. All Teff and γ are summarized in Table S1 of the Supplement.

 figure: Fig. 3.

Fig. 3. (a) Effective temperatures at q = 0.125 nm-1 obtained from KWW model fitting plotted as a function of fluences. Dashed lines serve as visual indicators. All temperatures are normalized such that the uniform ones equal 293 K. All standard deviations are less than 0.10 K. (b) Compression factors at q = 0.125 nm-1 obtained from KWW model fitting plotted as a function of fluences. Dashed lines serve as visual indicators. All standard deviations are less than or equal to 0.01.

Download Full Size | PDF

The above results confirm our presumption about the thermal effects of beam profiles. Nonuniform profiles apply different local fluences to different regions, such that the statistical weights of particles in hotter regions are higher in the ISFs. Therefore, results of nonuniform beams show faster dynamics and higher degrees of nonequilibrium. The dynamics parameters Teff and γ scale not only with H but also with the nonuniformity of the beams. The Gaussian profile brings about the highest Teff and the lowest γ compared with the uniform and donut-like profiles.

3.2 Γ-q2 dependence

Figure 4 shows the q2-dependence of the relaxation rates Γ for all beam profiles at different fluences H. Firstly and expectedly, Γ increase with H owing to stronger heating. They are slightly smaller than the experimental ones of the same H because of the neglect of photon absorption by particles [15]. In addition, they also increase with the nonuniformity of the beam profiles, such that the Gaussian ones are always the largest at the same H. These results agree with the findings above that the effective temperatures increase with H and the beam nonuniformity. Moreover, small deviations from Eq. (9) appear at small and large q for all fluences including H = 0. The result of H = 0 suggests that these deviations stem neither from heating-induced nonequilibrium nor from averaging over pulses. Further simulations demonstrate that they result from the spread of particle radial positions on the XZ plane (see Fig. S2 in the Supplement). Similar deviations were also observed experimentally [15]. However, the experimental Γ of H = 56.8 mJ/mm2 therein deviates significantly from Eq. (9). It could be attributed to the high degree of nonequilibrium induced by the strong photon absorption of particles, which is not included here. The numerical studies in this work are thus more instructive for soft-matter particle dispersions, where the photon absorption by particles is on the same level as water.

 figure: Fig. 4.

Fig. 4. The q2-dependence of the relaxation rates Γ for all beam profiles at different fluences H. The fluences and beam profiles are listed in the legend. Dashed lines are the Γ = Dq2 fits using the effective temperatures of q = 0.125 nm-1 (q2 = 0.016 nm-2). The highest and lowest effective temperatures are given in black. Standard deviations of all data points are less than 100 s-1.

Download Full Size | PDF

3.3 Thermal effects at different beam sizes

The above simulations and analysis were performed at the Gaussian beam size of d = 4 µm, where beam nonuniformity induces faster dynamics and higher degrees of nonequilibrium. It would be interesting to see how they change with the beam sizes. To this end, the Gaussian beam size d was varied. The other sizes were also varied proportionally to d. The thermal relaxation times of different beam sizes and profiles were calculated using the value and the dependence in the Refs. [15,23]. The pulse energies were kept constant, which is the common situation in experiments.

Figure 5 shows the effective temperatures Teff of different beam profiles at H = 57.7 mJ/mm2 plotted as a function of d. In the tested d range, Teff increase with d for all beam profiles. The increase is significant at smaller d but minimal at larger d. It indicates that the thermal effects induced by photon absorption could be reduced by using smaller beams. Regarding different beam profiles, the finding that nonuniform beams yield higher Teff remains valid. Teff of the Gaussian beams are always the highest, whereas Teff of the uniform beams are always the lowest. However, for smaller beams, the differences in Teff of different beam profiles are also smaller. For $d = \; 8\;\mathrm{\mu m}$, the differences are about twice as large as those of $d = \; 4\;\mathrm{\mu m}$. For $d = 1\;\mathrm{\mu m}$, the differences are almost negligible. Therefore, decreasing the beam size could be a practical strategy for MHz XPCS experiments to reduce beam-induced thermal effects, in particular the effects of beam profiles.

 figure: Fig. 5.

Fig. 5. The effective temperatures of different beam profiles plotted as a function of the Gaussian beam size d. All the other sizes are proportional to d. Simulations and analysis were performed at q = 0.125 nm-1. When varying d, the pulse energies were kept constant using the value of H = 57.7 mJ/mm2 in [15]. Dashed lines serve as visual indicators. The thermal relaxation time was calculated using the one in [15] and its dependence on d in [23]. All standard deviations are smaller than 0.15 K. Corresponding compression factors are between 0.89 and 0.96.

Download Full Size | PDF

4. Conclusion

XPCS at high repetition-rate XFELs such as the European XFEL offers unique capabilities and is thus under rapid development. However, critical issues such as the heating problem must first be better understood. The motivation of this work is to assess the thermal effects of different beam profiles, which can be altered by optics. The heat dissipation of water is relatively slow at high repetition rates. For the uniform beam, all particles are subject to nearly the same temperature, such that the statistic weights are equal for all particles. For nonuniform beams, however, heating is stronger in the regions with higher local fluences, such that ‘hotter’ particles obtain higher statistic weights in the correlation functions. Therefore, nonuniform beams could presumably facilitate faster effective dynamics. To verify this presumption, we performed Brownian motion simulations of particle dispersions under the illumination of differently-shaped beams at several fluences, where photon absorption by particles was neglected. Particle dynamics was revealed by correlation analysis and model fitting.

Results show that the effective temperatures of nonuniform beam profiles are indeed significantly higher than the uniform ones. Further analyses suggest that the dynamics of all beam profiles follows the Γ-q2 dependence. Additional studies of various beam sizes show that decreasing the beam size reduces the thermal effects induced by photon absorption of the beam profiles. Moreover, it diminishes the differences in the effective dynamics of different beam profiles.

For most of the simulations in this work, the differences in the effective dynamics of different profiles are only significant at high H, where radiation damage and particle heating are of greater concern. The thermal effects of nonuniform profiles would be critical only if significant differences occur at low H. Solutions with higher absorption or higher viscosity-temperature dependence are likely to facilitate such significant differences. Should proper thermal relaxation time be available, similar simulations will be performed for these two scenarios as well.

The numerical studies in this work cannot fully describe experimental conditions. The main limitation of this work is the neglect of photon absorption by particles. It depends on many factors such as absorbance, heat capacity, particle geometry, and convective heat transfer between particles and solvent [15,16]. While considering the photon absorption by particles would better describe the experiments, it would also include a great amount of complexity. In this first study, we have focused on the thermal effects of photon absorption by the solvent. Another limitation of our numerical studies is that heat mainly dissipates in the radial directions, such that the temperature distribution of the uniform profile is not entirely uniform. Nonetheless, this work offers a glimpse into the thermal effects of beam profiles on XPCS at high repetition-rate XFELs. Moreover, similar simulations could also be performed for XPCS studies at 4th generation synchrotrons. As they provide highly brilliant and coherent X-ray beams, they are also attractive sources for XPCS studies on short time scales. The average brilliance of 4th generation synchrotrons is only a few orders of magnitude lower than MHz XFELs [24], such that the heating problem is also worth noting. This work investigates the thermal effects of three common beam profiles. The findings could help select proper optics and beam profiles to reduce beam-induced sample heating in high repetition-rate XPCS experiments at these two highly brilliant X-ray facilities.

Funding

Natural Science Foundation of Shanghai (22ZR1442100).

Acknowledgments

The authors thank Dr. Felix Lehmkühler, Dr. Yajun Tong, Dr. Fan Yang, and Dr. Menglu Hu for their helpful discussions. This work was sponsored by Natural Science Foundation of Shanghai (Grant No. 22ZR1442100).

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Supplemental document

See Supplement 1 for supporting content.

References

1. B. Chu, Q. C. Ying, F. J. Yeh, A. Patkowski, W. Steffen, and E. W. Fischer, “An X-ray Photon Correlation Experiment,” Langmuir 11(5), 1419–1421 (1995). [CrossRef]  

2. S. B. Dierker, R. Pindak, R. M. Fleming, I. K. Robinson, and L. Berman, “X-Ray Photon Correlation Spectroscopy Study of Brownian Motion of Gold Colloids in Glycerol,” Phys. Rev. Lett. 75(3), 449–452 (1995). [CrossRef]  

3. O. G. Shpyrko, “X-ray photon correlation spectroscopy,” J. Synchrotron Radiat. 21(5), 1057–1064 (2014). [CrossRef]  

4. A. Madsen, A. Fluerasu, and B. Ruta, “Structural Dynamics of Materials Probed by X-Ray Photon Correlation Spectroscopy,” in Synchrotron Light Sources and Free-Electron Lasers, E. Jaeschke, S. Khan, J. R. Schneider, and J. B. Hastings, eds. (Springer International Publishing, 2015), pp. 1–21.

5. F. Perakis, G. Camisasca, T. J. Lane, A. Späh, K. T. Wikfeldt, J. A. Sellberg, F. Lehmkühler, H. Pathak, K. H. Kim, K. Amann-Winkel, S. Schreck, S. Song, T. Sato, M. Sikorski, A. Eilert, T. McQueen, H. Ogasawara, D. Nordlund, W. Roseker, J. Koralek, S. Nelson, P. Hart, R. Alonso-Mori, Y. Feng, D. Zhu, A. Robert, G. Grübel, L. G. M. Pettersson, and A. Nilsson, “Coherent X-rays reveal the influence of cage effects on ultrafast water dynamics,” Nat. Commun. 9(1), 1917 (2018). [CrossRef]  

6. V. M. Giordano and B. Ruta, “Unveiling the structural arrangements responsible for the atomic dynamics in metallic glasses during physical aging,” Nat. Commun. 7(1), 10344 (2016). [CrossRef]  

7. H. Zhou, S. Hilke, E. Pineda, M. Peterlechner, Y. Chushkin, S. Shanmugam, and G. Wilde, “X-ray photon correlation spectroscopy revealing the change of relaxation dynamics of a severely deformed Pd-based bulk metallic glass,” Acta Mater. 195, 446–453 (2020). [CrossRef]  

8. J. J. Song, R. Bhattacharya, H. Kim, J. Chang, T. Y. Tang, H. Guo, S. K. Ghosh, Y. Yang, Z. Jiang, H. Kim, T. P. Russell, G. Arya, S. Narayanan, and S. K. Sinha, “One-Dimensional Anomalous Diffusion of Gold Nanoparticles in a Polymer Melt,” Phys. Rev. Lett. 122(10), 107802 (2019). [CrossRef]  

9. F. Lehmkühler, J. Valerio, D. Sheyfer, W. Roseker, M. A. Schroer, B. Fischer, K. Tono, M. Yabashi, T. Ishikawa, and G. Grübel, “Dynamics of soft nanoparticle suspensions at hard X-ray FEL sources below the radiation-damage threshold,” IUCrJ 5(6), 801–807 (2018). [CrossRef]  

10. J. Möller, M. Sprung, A. Madsen, and C. Gutt, “X-ray photon correlation spectroscopy of protein dynamics at nearly diffraction-limited storage rings,” IUCrJ 6(5), 794–803 (2019). [CrossRef]  

11. A. Girelli, H. Rahmann, N. Begam, A. Ragulskaya, M. Reiser, S. Chandran, F. Westermeier, M. Sprung, F. Zhang, C. Gutt, and F. Schreiber, “Microscopic Dynamics of Liquid-Liquid Phase Separation and Domain Coarsening in a Protein Solution Revealed by X-Ray Photon Correlation Spectroscopy,” Phys. Rev. Lett. 126(13), 138004 (2021). [CrossRef]  

12. O. G. Shpyrko, E. D. Isaacs, J. M. Logan, Y. Feng, G. Aeppli, R. Jaramillo, H. C. Kim, T. F. Rosenbaum, P. Zschack, M. Sprung, S. Narayanan, and A. R. Sandy, “Direct measurement of antiferromagnetic domain fluctuations,” Nature 447(7140), 68–71 (2007). [CrossRef]  

13. S. Gorfman, A. A. Bokov, A. Davtyan, M. Reiser, Y. Xie, Z.-G. Ye, A. V. Zozulya, M. Sprung, U. Pietsch, and C. Gutt, “Ferroelectric domain wall dynamics characterized with X-ray photon correlation spectroscopy,” Proc. Natl. Acad. Sci. U.S.A. 115(29), E6680–E6689 (2018). [CrossRef]  

14. G. Grübel, G. B. Stephenson, C. Gutt, H. Sinn, and T. Tschentscher, “XPCS at the European X-ray free electron laser facility,” Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 262(2), 357–367 (2007). [CrossRef]  

15. F. Lehmkühler, F. Dallari, A. Jain, M. Sikorski, J. Möller, L. Frenzel, I. Lokteva, G. Mills, M. Walther, H. Sinn, F. Schulz, M. Dartsch, V. Markmann, R. Bean, Y. Kim, P. Vagovic, A. Madsen, A. P. Mancuso, and G. Grübel, “Emergence of anomalous dynamics in soft matter probed at the European XFEL,” Proc. Natl. Acad. Sci. USA 117(39), 24110–24116 (2020). [CrossRef]  

16. F. Dallari, A. Jain, M. Sikorski, J. Möller, R. Bean, U. Boesenberg, L. Frenzel, C. Goy, J. Hallmann, Y. Kim, I. Lokteva, V. Markmann, G. Mills, A. Rodriguez-Fernandez, W. Roseker, M. Scholz, R. Shayduk, P. Vagovic, M. Walther, F. Westermeier, A. Madsen, A. P. Mancuso, G. Grübel, and F. Lehmkühler, “Microsecond hydrodynamic interactions in dense colloidal dispersions probed at the European XFEL,” IUCrJ 8(5), 775–783 (2021). [CrossRef]  

17. M. Reiser, A. Girelli, A. Ragulskaya, S. Das, S. Berkowicz, M. Bin, M. Ladd-Parada, M. Filianina, H.-F. Poggemann, N. Begam, M. S. Akhundzadeh, S. Timmermann, L. Randolph, Y. Chushkin, T. Seydel, U. Boesenberg, J. Hallmann, J. Möller, A. Rodriguez-Fernandez, R. Rosca, R. Schaffer, M. Scholz, R. Shayduk, A. Zozulya, A. Madsen, F. Schreiber, F. Zhang, F. Perakis, and C. Gutt, “Resolving molecular diffusion and aggregation of antibody proteins with megahertz X-ray free-electron laser pulses,” (2022).

18. A. P. Mancuso, A. Aquila, L. Batchelor, R. J. Bean, J. Bielecki, G. Borchers, K. Doerner, K. Giewekemeyer, R. Graceffa, O. D. Kelsey, Y. Kim, H. J. Kirkwood, A. Legrand, R. Letrun, B. Manning, L. Lopez Morillo, M. Messerschmidt, G. Mills, S. Raabe, N. Reimers, A. Round, T. Sato, J. Schulz, C. Signe Takem, M. Sikorski, S. Stern, P. Thute, P. Vagovič, B. Weinhausen, and T. Tschentscher, “The Single Particles, Clusters and Biomolecules and Serial Femtosecond Crystallography instrument of the European XFEL: initial installation,” J Synchrotron Rad 26(3), 660–676 (2019). [CrossRef]  

19. A. Madsen, J. Hallmann, G. Ansaldi, T. Roth, W. Lu, C. Kim, U. Boesenberg, A. Zozulya, J. Möller, R. Shayduk, M. Scholz, A. Bartmann, A. Schmidt, I. Lobato, K. Sukharnikov, M. Reiser, K. Kazarian, and I. Petrov, “Materials Imaging and Dynamics (MID) instrument at the European X-ray Free-Electron Laser Facility,” J Synchrotron Rad 28(2), 637–649 (2021). [CrossRef]  

20. M. Vassholz, H. P. Hoeppe, J. Hagemann, J. M. Rosselló, M. Osterhoff, R. Mettin, T. Kurz, A. Schropp, F. Seiboth, C. G. Schroer, M. Scholz, J. Möller, J. Hallmann, U. Boesenberg, C. Kim, A. Zozulya, W. Lu, R. Shayduk, R. Schaffer, A. Madsen, and T. Salditt, “Pump-probe X-ray holographic imaging of laser-induced cavitation bubbles with femtosecond FEL pulses,” Nat. Commun. 12(1), 3468 (2021). [CrossRef]  

21. J. Kestin, M. Sokolov, and W. A. Wakeham, “Viscosity of liquid water in the range −8 °C to 150 °C,” J. Phys. Chem. Ref. Data 7(3), 941–948 (1978). [CrossRef]  

22. J. Stetefeld, S. A. McKenna, and T. R. Patel, “Dynamic light scattering: a practical guide and applications in biomedical sciences,” Biophys. Rev. 8(4), 409–427 (2016). [CrossRef]  

23. J. Carnis, W. Cha, J. Wingert, J. Kang, Z. Jiang, S. Song, M. Sikorski, A. Robert, C. Gutt, S.-W. Chen, Y. Dai, Y. Ma, H. Guo, L. B. Lurio, O. Shpyrko, S. Narayanan, M. Cui, I. Kosif, T. Emrick, T. P. Russell, H. C. Lee, C.-J. Yu, G. Grübel, S. K. Sinha, and H. Kim, “Demonstration of Feasibility of X-Ray Free Electron Laser Studies of Dynamics of Nanoparticles in Entangled Polymer Melts,” Sci. Rep. 4(1), 6017 (2015). [CrossRef]  

24. T. S. Plivelic, “X-ray tracing, design and construction of an optimized optics scheme for CoSAXS, the small angle x-ray scattering beamline at MAX IV laboratory: AIP Conference Proceedings: Vol 2054, No 1,” https://aip.scitation.org/doi/abs/10.1063/1.5084576.

Supplementary Material (1)

NameDescription
Supplement 1       Supplement

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (5)

Fig. 1.
Fig. 1. (a) The typical particle distribution in the transverse probing region of 8 µm diameter, as indicated by the red dashed circle. Every blue sphere denotes one particle. X and Z are the two transverse directions perpendicular to the beam. (b)(c)(d) The adiabatic temperature increases of water induced by each uniform, donut-like, and Gaussian pulse. The fluence H is 10.5 mJ/mm2. The color bar indicates the temperature increase in K.
Fig. 2.
Fig. 2. Simulated intermediate scattering functions (symbols) and their fitting curves (solid lines) at q = 0.125 nm-1. The fluences and beam profiles are shown in the legend.
Fig. 3.
Fig. 3. (a) Effective temperatures at q = 0.125 nm-1 obtained from KWW model fitting plotted as a function of fluences. Dashed lines serve as visual indicators. All temperatures are normalized such that the uniform ones equal 293 K. All standard deviations are less than 0.10 K. (b) Compression factors at q = 0.125 nm-1 obtained from KWW model fitting plotted as a function of fluences. Dashed lines serve as visual indicators. All standard deviations are less than or equal to 0.01.
Fig. 4.
Fig. 4. The q2-dependence of the relaxation rates Γ for all beam profiles at different fluences H. The fluences and beam profiles are listed in the legend. Dashed lines are the Γ = Dq2 fits using the effective temperatures of q = 0.125 nm-1 (q2 = 0.016 nm-2). The highest and lowest effective temperatures are given in black. Standard deviations of all data points are less than 100 s-1.
Fig. 5.
Fig. 5. The effective temperatures of different beam profiles plotted as a function of the Gaussian beam size d. All the other sizes are proportional to d. Simulations and analysis were performed at q = 0.125 nm-1. When varying d, the pulse energies were kept constant using the value of H = 57.7 mJ/mm2 in [15]. Dashed lines serve as visual indicators. The thermal relaxation time was calculated using the one in [15] and its dependence on d in [23]. All standard deviations are smaller than 0.15 K. Corresponding compression factors are between 0.89 and 0.96.

Equations (9)

Equations on this page are rendered with MathJax. Learn more.

$${H_{\textrm{eff}}} = H\textrm{ }Ts\textrm{ }[1 - \exp ( - W/L)], $$
$$\Delta {T_U} = \frac{{{H_{\textrm{eff}}}}}{{4W\rho {c_p}}},$$
$$\Delta {T_D}(r) = \frac{{{H_{\textrm{eff}}}}}{{4W\rho {c_p}}}N{F_D} \exp [ - \frac{{{{(r - \mu )}^2}}}{{2{\sigma ^2}}}], $$
$$\Delta {T_G}(r) = \frac{{{H_{\textrm{eff}}}}}{{4W\rho {c_p}}}N{F_G}\exp ( - \frac{{{r^2}}}{{2{\sigma ^2}}}), $$
$$D = \frac{{{k_b}T}}{{6\pi \eta R}}, $$
$$E(q) = \sum\nolimits_{i = 1}^{Np} {A({r_i})\exp ( - iq{r_i})}, $$
$$f(q,\tau ) = \frac{{ < E(q,t)E(q,t + \tau ) > }}{{ < E(q,t){ > ^2}}}. $$
$$f(q,\tau ) = \exp ( - |\Gamma (q)\tau {|^\gamma }), $$
$$\Gamma (q) = \frac{{{k_b}{T_{\textrm{eff}}}}}{{6\pi \eta R}}{q^2}. $$
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.