Millisecond dynamics of colloidal suspension studied by X‑ray photon correlation spectroscopy at the Shanghai Synchrotron Radiation Facility

Chen‑Hui Cui · Zi‑Mu Zhou · Lin‑Feng Wei · Song‑Lin Li · Feng Tian · Xiu‑Hong Li · Zhi Guo · Yi‑Hui Xu ·Huai‑Dong Jiang2, · Ren‑Zhong Tai,2

Abstract X-ray photon correlation spectroscopy (XPCS) has emerged as a powerful tool for probing the nanoscale dynamics of soft condensed matter and strongly correlated materials owing to its high spatial resolution and penetration capabilities.This technique requires high brilliance and beam coherence, which are not directly available at modern synchrotron beamlines in China.To facilitate future XPCS experiments, we modified the optical setup of the newly commissioned BL10U1 USAXS beamline at the Shanghai Synchrotron Radiation Facility (SSRF).Subsequently, we performed XPCS measurements on silica suspensions in glycerol, which were opaque owing to their high concentrations.Images were collected using a high frame rate area detector.A comprehensive analysis was performed, yielding correlation functions and several key dynamic parameters.All the results were consistent with the theory of Brownian motion and demonstrated the feasibility of XPCS at SSRF.Finally, by carefully optimizing the setup and analyzing the algorithms, we achieved a time resolution of 2 ms, which enabled the characterization of millisecond dynamics in opaque systems.

Keywords XPCS · SSRF · Silica suspension · Glycerol · Brownian motion · Millisecond dynamics

1 Introduction

When coherent light impinges on a disordered sample, the resulting scattered signal forms a speckle pattern.Because the internal structure of the disordered sample evolves over time, the corresponding speckle pattern also fluctuates.Valuable insights into the internal dynamics of the sample can be obtained by analyzing the correlations between these fluctuating speckle patterns.Using this principle, dynamic light scattering (DLS) techniques that utilize coherent visible light have been developed for approximately half a century, facilitating the exploration of numerous intricate dynamics and relaxation processes in condensed systems [1–5].A fundamental assumption during the correlation of the operational scattering signals is that the photons detected by the detector are scattered only once by the sample, thereby excluding the effect of multiple scattering events.However, in the visible light spectrum, the probability of a single photon undergoing multiple scatterings becomes significant when the concentration of the sample under investigation is excessively high or the sample is optically opaque [6].Therefore, DLS is typically utilized only for analyzing materials with low concentrations or minimal absorption.

Owing to this limitation, XPCS techniques that use (partially) coherent X-rays as a substitute for visible light have emerged, thereby enabling the examination of highly concentrated or optically opaque materials [7].This is primarily because the scattering cross section of X-rays is smaller than that of visible light [8].This results in reduced susceptibility to multiple scattering during the experiment.In addition,the shorter wavelength of X-rays enables the detection of samples with smaller length scales [9].However, the smaller scattering cross section of the X-rays also results in a weaker scattering signal.When the number of coherent photons in the (partially) coherent X-ray beam is equivalent to that in a laser, XPCS experiments become more challenging than DLS experiments, resulting in poor signal-to-noise ratio(SNR) of the data obtained.Consequently, the coherent X-ray flux is a primary limiting factor for the development of XPCS.The solution can be derived from the SNR equation [10]:

whereIpixdenotes the average intensity per pixel, which is associated with the incident coherent flux.Npix,Nfr, andNrepdenote the number of pixels, number of frames, and number of repetitions, respectively.Ipixcan be enhanced by designing a dedicated XPCS beamline (for instance, focusing the beam using compound refractive lenses).Likewise, the coherent X-ray flux is anticipated to increase by several orders of magnitude by adopting more advanced fourth-generation X-ray sources, such as near diffraction-limited storage rings(DLSRs) [11–16] and free-electron lasers (FELs) [17–22].Npixcan be improved using area detectors along with the multi-speckle XPCS method [23].BothNrepandNfrenhance the SNR by increasing the number of sampled frames.

Another limitation of the XPCS is its relatively low temporal resolution that primarily arises from the lower (partially) coherent flux employed in XPCS experiments.For instance, at a constant coherent flux, increasing the detector frame rate requires a corresponding reduction in the exposure time, which in turn lowers the value ofIpix.According to Eq.1, this decrease inIpixreduces the SNR.Therefore,both high SNR and high temporal resolution are impossible to achieve simultaneously, as these two factors are inversely related to one another.The temporal resolution of XPCS was constrained to ensure that XPCS experimental data had a sufficiently high SNR.Furthermore, the maximum detection frame rate of the area detectors employed in XPCS experiments presented a notable constraint [7, 9].With the development of fourth-generation X-ray sources and the introduction of higher-frame-rate detectors, the range of XPCS detection times is expected to expand, particularly for the fastest dynamics [9, 24, 25].After more than two decades of development, XPCS has been extensively applied to condensed systems [7, 9], including colloidal suspensions [23,26], polymers [27, 28], metallic glasses [29, 30], proteins[31–33], and strongly correlated materials [10, 34, 35].

SSRF is the first third-generation synchrotron radiation light source in China.In recent years, its experimental capability has increased [36–40].More recently, the SSRF commissioned all its second-phase beamlines that employ state-of-the-art undulators to provide significantly higher beam quality.This endeavor will facilitate the development of coherent technology that has already achieved results in coherent diffraction imaging (CDI) and scanning coherent diffraction imaging (scanning CDI or ptychography) at the SSRF [41–43].However, significant progress is yet to be achieved in XPCS experiments.In addition, the absence of a dedicated XPCS beamline at SSRF presents challenges for conducting XPCS experiments.

In this work, we demonstrated the feasibility of XPCS at the SSRF with millisecond time resolution by optimizing the optics, samples, and analysis algorithm.The experiments employed standard silica colloidal suspensions as samples,which were ideal for XPCS experiments owing to their relatively simple and clear dynamic mechanisms.Furthermore,colloidal particles can serve as tracer particles to reveal structural changes in the sample through their dynamic alterations, such as glass transition [44, 45] and non-equilibrium dynamics of soft materials [46].Consequently, understanding the dynamic processes of dilute colloidal suspensions is beneficial for further investigations of related materials.

2 Sample preparation

An aqueous suspension of silica spheres (M814153, MACKLIN, 500 nm ) was used as the stock solution in this experiment.Characterized by scanning electron microscopy(SEM), the diameter of the silica spheres of the stock solution was determined by evaluating at least 1000 particles in the image.An average size ofRSEM=524.4 nm and dispersion ofp=7.2% were obtained [26].The volume fraction of silica in the stock solution was 1%.But the diffusion of silica in water was too rapid for the Eiger X 4 M detector in the range of the observable.To better characterize the dynamics, the stock solution was mixed with glycerol at a volume ratio of 1:1.Most of the water was removed using a rotary evaporator, yielding considerably slower diffusion[44, 47].The actual volume fraction of silica was slightly lower than 1% because the rotary evaporator could not completely remove the residual water.The interactions between silica particles can be neglected at this volume fraction [45].Before the XPCS experiment, the suspension was sonicated and injected into a 2 mm-thick sample cell, which was sealed with Kapton film.The experiment was conducted at a temperature of approximately ～298K.

Fig.1 (Color online) Schematic of the experimental setup: 1.Undulator; 2.Four-blade slits 1; 3.Horizontal deflection mirror; 4.Double multilayer monochromator; 5.Horizontal focusing mirror; 6.Vertical focusing mirror; 7.Four-blade slits 2; 8.Four-blade slits 3; 9.Sample; 10.Detector.[48]

3 Experimental configuration

The experiment was conducted at the BL10U1 USAXS beamline of SSRF, and the schematic diagram of the experimental setup is shown in Fig.1.The incident beam emitted from the undulator was monochromated by a double multilayer monochromator after passing through the four-blade slits 1.The photon energy and energy resolution ( ΔE∕E)were 10 keV (λ=0.124 nm ) and 5.8×10−3, respectively.The beam was subsequently focused using horizontal and vertical focusing mirrors with the focal point at the detector.The partially coherent light required for the experiment was selected using upstream four-blade slits 2 and 3.A pinhole with a diameter of 100 μm was placed upstream of the sample to reduce the beam size and increase the coherence of the incident beam on the sample.The scattering signal of the sample was recorded using an area detector (EIGER X 4 M)located 27.6m downstream from the sample.The effective pixel number of the detector was 2070×2167 with a pixel size of 75 μm×75 μm.Its frame rate could reach 750Hz ,i.e., 1.34 ms minimum exposure period.In this experiment,the availableq-range (scattering wave vector range) was approximately 1.2×10−2nm−1–2.0×10−1nm−1.The scattering wave vectorqwas determined from the scattering angle 2θusing the formulaq=4πsin(θ)∕λ.To ensure the repeatability of the experiment and accuracy of the results,the experiments were repeated three times under the same conditions, with each run collecting 1000 speckle patterns.Six exposure periods (2, 5, 10, 20, 50, and 100 ms) were used to evaluate the temporal resolving ability of the XPCS study.

4 DLS

Prior to the XPCS experiment, the silica suspension was pre-characterized using DLS technology.The coherent light source was a laser with a wavelength ofλlaser=633 nm , and the size of the laser beam was adjusted to 2 mm using a diaphragm.The laser was then directed into the sample contained in the cuvette, causing the light to scatter and form a speckle pattern.The scattered signal was collected by a detector (photomultiplier tube, PMT) after passing through an order-sorting aperture (OSA) with a diameter of 70 μm.The PMT was located at a scattering angle ofπ∕2 ( 90◦).The signal was displayed on a computer after correlation using a correlator (ALV-7004/USB-FAST).The entire experimental setup was placed in a dark room to avoid the influence of natural light.Additionally, considering the high sensitivity of the PMT to light, an extra layer of darkroom was added around the detector to prevent stray signals from the laser and allow only the OSA to receive light.The acquisition time for each signal was set to 10 s, and the process was repeated ten times to obtain the average value.

The stock solution was diluted with water to different concentrations, and then, these samples were placed in the cuvette for DLS measurement (Fig.2).The experimental temperature was maintained at 297 K.The measured results were fitted (Fig.2, solid line) using the equationto obtain the relaxation timeτ, whereβrepresents the speckle contrast.According to the experimental parameters, the hydrodynamic diameters of silica (Fig.3)can be obtained from the relaxation timesτof the samples.Figure 3 shows that in the case of a low dilution factor, the hydrodynamic diameter of silica first increases and then decreases owing to the influence of multiple scattering.The influence of multiple scattering becomes negligible with increasing dilution factor.The hydrodynamic diameter of silica is stable at approximately 550 nm , which is slightly larger than that obtained by SEM.This is consistent with the expected results [23].

Fig.2 (Color online) Correlation data (points) obtained using DLS measurements of the stock solution at different dilution factors or concentrations ( X10 ～0.1% , X20 ～0.05% , X50 ～0.02% ,X100 ～0.01% , X200 ～0.005% , X250 ～0.004% , X500 ～0.002% ,X1000 ～0.001%.).These solid lines are the fitting curves for the correlation data

Fig.3 Hydrodynamic diameters of silica in water suspensions of different concentrations obtained using the fitting results.The x-axis represents the dilution factor of the stock solution

These findings indicate that DLS is indeed constrained by the sample concentration and that multiple scattering significantly affects the experimental results.Therefore, determining the dynamic behavior of silica in the stock solution using the traditional DLS method is impossible because of the influence of multiple scattering.However, XPCS can easily overcome the influence of multiple scattering using hard X-rays.

5 Small angle X‑ray scattering (SAXS)

The two-dimensional (2D) small angle X-ray scattering(SAXS) data presented in Fig.4 were generated by overlaying 1000 frames of speckle patterns collected for each of the three sets of data.The results from the three sets of data were averaged.The exposure period was 50 ms , corresponding to a total exposure time of 50 s per frame.The 2D SAXS pattern was azimuthally integrated to obtain the one-dimensional (1D) SAXS profile (Fig.5 point), illustrating the scattering intensityI(q) with respect to the scattering wave vectorq.The signals of bad pixels, gaps, and beamstop were excluded from data processing.By fitting the 1D SAXS data (Fig.5 line), the diameter of the silica particles was determined to beRSAXS=528.2 nm (In the highqrange,the difference between the SAXS results and fitted curve is primarily because of the weak scattered signal at highq, which is considerably affected by noise.).This size was consistent with that obtained by SEM, indicating that the size of the silica particles produced via rotary evaporation did not change significantly, and substantial agglomeration was not observed.The agreement between the SAXS, SEM,and DLS measurements proves the reliability of the size and structure of the silica particles in the stock solutions.Therefore, a comprehensive understanding of the particle size and potential aggregation phenomena was obtained by employing the SAXS and DLS techniques.

Fig.4 (Color online) Two-dimensional (2D) SAXS from the sample overlaying 1000 frames collected continuously within 50 s.The result was further averaged from three repeated measurements to improve the statistics.The x- and y-axes are represented by the scattering vectors qx and qy , respectively.For a clear view, the logarithm of the scattering signal intensity is taken to the base 10.The illustration in the upper right corner is a speckle image

Fig.5 (Color online) One-dimensional (1D) SAXS data derived by using the azimuthal averaging of scattering signals from the sample shown in Fig.4.Point: SAXS scattering intensity I(q) with respect to scattering wave vector q; Line: fitting of silica size (considering the silica as homogeneous spherical particles) [49].For a clear view, the logarithm of the scattering signal intensity is taken to the base 10

6 X‑ray photon correlation spectroscopy(XPCS)

6.1 Theory

To investigate the dynamic properties of silica, the obtained scattering patterns were correlated using the intensity correlation equation.The influence of photon flux variation on the results was minimized by applying the two-time correlation (TTC) method to analyze the sample dynamics [50]:

whereI(q,t1,2) is the intensity value atqandt1,2; 〈…〉 represents the average of all pixels; and ||t1−t2||=τindicates the delay time.When the sample is in equilibrium during the observation time, the TTC can be expressed as [25, 51]:

Averaging the corresponding data with the same delay timeτimproves the SNR of the results.The intensity correlation functiong2(q,τ) can be related to the first-order correlation functiong1(q,τ) through the Siegert relation [1]:

Here,βrepresents the speckle contrast, which depends on the coherence of the incident beam, optical layout of the beamline and pixel size of the detector.In this study, the diluted colloidal suspension was expected to exhibit Brownian motion [26, 51, 52], such that

Combining Eqs.4 and 5, the correlation equation can be simplified to:

Γ(q)=1∕τc(q) represents the relaxation rate, which is the reciprocal of the characteristic relaxation timeτc(q):

Drepresents the Stokes–Einstein diffusion constant, which can be expressed as

wherekB,T,η, andRHrepresent the Boltzmann constant,thermodynamic temperature, viscosity coefficient of the solvent, and hydrodynamic radius of the silica particles,respectively.

6.2 Results and discussion

The two time correlation (TTC) maps in Fig.6 exhibit no significant change (narrowing) at the midridge line on the diagonal in the image during the entire experiment.This indicates that the silica in the sample is in equilibrium dynamics.Therefore, the effects of sample heating and damage by X-rays on the sample were negligible for the XPCS experiments conducted in this work.

The data collected from the same sample at different exposure periods were correlated and fitted using Eq.6(Fig.7).First, all the data were fitted well using a single exponential model.Second, data with shorter exposure periods exhibited lower SNR.From Eq.1, it is natural to obtain lower SNR when lower exposure periods reduced the intensityIpixwith all other parameters kept constant.Nonetheless,the SNR of the lowest exposure period (2 ms) was sufficient for the correlation analysis.The SNR of the low-exposure periods could be significantly increased by increasingNfrorNrep.Third, and most importantly, despite the significantly different exposure periods, all correlation curves exhibited similar decay behaviors.This is further demonstrated by the relaxation rate analysis below.

Fig.6 (Color online) TTC derived from the correlation calculation of 1000 frames of speckle patterns using Eq.2.t1,2 indicates the time of the experiment, where the exposure period is 50 ms,q=2.13×10−2 nm−1 ( Δq=2.06×10−3 nm−1 ).The inset figure shows a partially enlarged view of the results of the experimental time 25–35 s

Fig.7 (Color online) Correlation function of the sample at different exposure periods, where q=2.13×10−2 nm−1( Δq=2.07×10−3 nm−1 ), and g2(q,τ) is normalized.For clarity,g2(q,τ) was shifted vertically.Points and lines refer to correlation data and fitted curves using Eq.6, respectively.These results were averaged over three replicated measurements to improve the statistics

Fig.8 Relaxation rate Γ obtained by fitting the correlation function at different exposure periods (Fig.7), and the abscissa represents the different exposure periods (2, 5, 10, 20, 50, and 100 ms).The dashed line is the guide line for the eye.These results were averaged over three replicated measurements to improve the statistics

The relaxation rate Γ of the sample obtained from fitting is shown in Fig.8, and the relaxation rates obtained at different exposure periods are consistent.The dashed line is the guide line for the eye, and the error bars are the standard deviations of the three repetitions.The Γ values for different exposure periods are almost unchanged, indicating the feasibility of XPCS experiments with millisecond exposures in SSRF [47, 53–56].In addition, Fig.8 shows that the standard deviation decreases with increasing exposure time.From Eq.1, this can be attributed to an increase inIpix, which improves the SNR of the result.Therefore, we can improve the SNR by increasing the number of frames (Nfr) collected,thereby reducing the standard deviation of the millisecond exposure results.In summary, the standard deviation of our millisecond-exposure results is usable and can be reduced to ensure realistic dynamic sampling.

The relationship betweeng2(q,τ) and differentqvalues in the data was correlated and fitted using Eq.6 (Fig.9).The decay is slower at lowerqvalues because the characteristic decay time decreases withq.Figure 10 shows that the relaxation rate (points) of the sample varies linearly withq2.The Stokes–Einstein diffusion constantD=2425 nm2∕s is obtained by linear fitting with Γ=D×q2.This proves that the silica in the sample undergoes Brownian motion.

Because of the incomplete water removal by the rotary evaporator and the hygroscopicity of glycerol, the sample inevitably contains water.The water content of glycerol strongly affects the viscosity [57], and thus, the viscosity of the sample cannot be accurately determined.However, for a silica sphere,RHis slightly larger thanRSAXS[23].We assumed thatRH=RSAXS[47] and then estimated the viscosity of the sample using Eq.8, which yieldedη=0.34 Pa s.Based on the empirical formula [58],

Fig.9 (Color online) Correlation functions under different scattering wave vector q values ( Δq=6.88×10−4 nm−1 ), where the exposure period is 50 ms, and g2(q,τ) is normalized.For clarity, g2(q,τ) was shifted vertically.Points and lines refer to correlation data and fitted curves using Eq.6 at different q values, respectively.These results were averaged over three replicated measurements to improve the statistics

Fig.10 Relationship between the relaxation rate Γ and q2.The points refer to the relaxation rate Γ under different q2 values obtained from the fitting of Fig.9.The line refers to the linear fitting of the points,which has been forced through the origin.These results were averaged over three replicated measurements to improve the statistics

aandbare the coefficients of the temperatureTcin Celsius,wherea=0.705 −0.0017Tcandb=(4.9+0.036Tc)a2.5.The water mass concentration in the sample was estimated to be approximately 5.2%.This portion of water, which is difficult to remove, makes it challenging to determine the viscosity of the sample.Consequently, it is difficult to calculate the hydrodynamic diameter of silica using the Stokes–Einstein diffusion constantDand compare it with the results obtained from SEM or SAXS.Using another solvent that absorbs less water can help solve this problem.Nonetheless,we successfully demonstrated a systematic XPCS study of a highly concentrated silica suspension at SSRF.

7 Summary

XPCS is an emerging technique with great potential for characterizing sample dynamics at molecular and atomic scales.Its advantages over DLS include high penetration capability, high spatial resolution, and reduced multiscattering.However, the requirement for coherent X-ray sources limits its accessibility, as third-generation synchrotrons such as SSRF, cannot easily provide sufficient coherence.After the recent upgrade of SSRF, the newly commissioned BL10U1 USAXS beamline provides an ideal platform for XPCS experiments because of its unique properties such as relatively high coherence, high brilliance, and long sample-to-detector distance.

In this study, we successfully performed XPCS studies on a silica suspension model system.The linear Γ −q2dependence demonstrates that the correct Brownian motion dynamics was obtained.Careful collimation of micrometer optics ensured low parasitic scattering and high coherence, and a newly established framework for data analysis enabled the effective extraction of dynamic information from scattering images.Consequently,we achieved a millisecond time resolution, which is at the same level as most XPCS experiments at modern synchrotrons.

The establishment of the XPCS platform is expected to provide new options for researchers in China and abroad and facilitate the study of various nanoscale and microscale samples such as colloidal suspensions, polymers,glass, biomacromolecules, and strongly related materials.This powerful tool will enable researchers gain valuable insights into the dynamic properties of these systems in their natural environment and contribute to a deeper understanding of their behavior and properties.

AcknowledgementsThe authors thank the staffof USAXS beamline (BL10U1) at SSRF for their support during the beamtime.The authors thank the user experimental assistance system of SSRF for their assistance.

Author ContributionsHuai-Dong Jiang, and Ren-Zhong Tai contributed to the overall research concept and supervision of this study.Yi-Hui Xu, Chen-Hui Cui, Zi-Mu Zhou, Lin-Feng Wei, and Song-Lin Li participated in material preparation, data collection, data processing program writing, and data analysis.Feng Tian, Xiu-Hong Li provided experimental suggestions and optimized experimental conditions.Zhi Guo and Yi-Hui Xu provided fund support and reviewed the manuscript.The first draft of the manuscript was written by Chen-Hui Cui,and all authors commented on previous versions of the manuscript.All authors read and approved the final manuscript.

Data Availability StatementThe data that support the findings of this study are openly available in Science Data Bank at https:// doi.org/ 10.57760/ scien cedb.13559 and https:// cstr.cn/ 31253.11.scien cedb.13559.

Declarations

Conflict of interestRen-Zhong Tai is an editorial board member for Nuclear Science and Techniques and was not involved in the editorial review, or the decision to publish this article.All authors declare that there are no competing interests.