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

Journal logoJOURNAL OF
APPLIED
CRYSTALLOGRAPHY
ISSN: 1600-5767

On the analysis of two-time correlation functions: equilibrium versus non-equilibrium systems

crossmark logo

aInstitute of Applied Physics, University of Tübingen, Auf der Morgenstelle 10, 72076, Tübingen, Germany, bCluster of Excellence Machine Learning: New Perspectives for Science, University of Tübingen, Maria-von-Linden-Straße 6, 72076, Tübingen, Germany, and cDepartment of Physics, University of Siegen, Emmy-Noether-Campus, Walter-Flex-Straße 3, 57076, Siegen, Germany
*Correspondence e-mail: fajun.zhang@uni-tuebingen.de, frank.schreiber@uni-tuebingen.de

Edited by A. Barty, DESY, Hamburg, Germany (Received 16 February 2024; accepted 16 May 2024; online 4 July 2024)

X-ray photon correlation spectroscopy (XPCS) is a powerful tool for the investigation of dynamics covering a broad range of timescales and length scales. The two-time correlation function (TTC) is commonly used to track non-equilibrium dynamical evolution in XPCS measurements, with subsequent extraction of one-time correlations. While the theoretical foundation for the quantitative analysis of TTCs is primarily established for equilibrium systems, where key parameters such as the diffusion coefficient remain constant, non-equilibrium systems pose a unique challenge. In such systems, different projections (`cuts') of the TTC may lead to divergent results if the underlying fundamental parameters themselves are subject to temporal variations. This article explores widely used approaches for TTC calculations and common methods for extracting relevant information from correlation functions, particularly in the light of comparing dynamics in equilibrium and non-equilibrium systems.

1. Introduction

X-ray photon correlation spectroscopy (XPCS) is a highly versatile experimental technique that is widely used to study the dynamics of both soft and hard condensed matter. The current state of XPCS and light sources enables the investigation of the dynamics across an unprecedented range of timescales, from femtoseconds to hours, and length scales that span from micrometres down to ångströms (Shpyrko, 2014[Shpyrko, O. G. (2014). J. Synchrotron Rad. 21, 1057-1064.]; Lehmkühler et al., 2021[Lehmkühler, F., Roseker, W. & Grübel, G. (2021). Appl. Sci. 11, 6179.]). XPCS is employed in various areas of condensed matter research to explore the dynamics of colloids (Westermeier et al., 2012[Westermeier, F., Fischer, B., Roseker, W., Grübel, G., Nägele, G. & Heinen, M. (2012). J. Chem. Phys. 137, 114504.]; Kwaśniewski et al., 2014[Kwaśniewski, P., Fluerasu, A. & Madsen, A. (2014). Soft Matter, 10, 8698-8704.]; Angelini et al., 2014[Angelini, R., Madsen, A., Fluerasu, A., Ruocco, G. & Ruzicka, B. (2014). Colloids Surf. A Physicochem. Eng. Asp. 460, 118-122.]; Angelini & Ruzicka, 2015[Angelini, R. & Ruzicka, B. (2015). Colloids Surf. A Physicochem. Eng. Asp. 483, 316-320.]; Liu et al., 2021[Liu, W., Zheng, B., Yin, X., Yu, X., Zhang, Y., Wiegart, L., Fluerasu, A., Armstrong, B. L., Veith, G. M. & Bhatia, S. R. (2021). Appl. Mater. Interfaces, 13, 14267-14274.]), liquids and liquid crystals (Seydel et al., 2001[Seydel, T., Madsen, A., Tolan, M., Grübel, G. & Press, W. (2001). Phys. Rev. B, 63, 073409.]; Lu et al., 2008[Lu, X., Mochrie, S. G. J., Narayanan, S., Sandy, A. R. & Sprung, M. (2008). Phys. Rev. Lett. 100, 045701.]; Madsen et al., 2003[Madsen, A., Als-Nielsen, J. & Grübel, G. (2003). Phys. Rev. Lett. 90, 085701.]; van 't Zand et al., 2012[Zand, D. D. van 't, Chushkin, Y., Belkoura, L., Lobo, C. V., Strey, R., Lyakhova, K. & Clegg, P. S. (2012). Soft Matter, 8, 4062-4066.]), polymers (Narayanan et al., 2007[Narayanan, S., Lee, D. R., Hagman, A., Li, X. & Wang, J. (2007). Phys. Rev. Lett. 98, 185506.]; Conrad et al., 2015[Conrad, H., Lehmkühler, F., Fischer, B., Westermeier, F., Schroer, M. A., Chushkin, Y., Gutt, C., Sprung, M. & Grübel, G. (2015). Phys. Rev. E, 91, 042309.]; Nogales & Fluerasu, 2016[Nogales, A. & Fluerasu, A. (2016). Eur. Polym. J. 81, 494-504.]), metallic and molecular glasses (Leitner et al., 2012[Leitner, M., Sepiol, B., Stadler, L.-M. & Pfau, B. (2012). Phys. Rev. B, 86, 064202.]; Ruta et al., 2012[Ruta, B., Chushkin, Y., Monaco, G., Cipelletti, L., Pineda, E., Bruna, P., Giordano, V. M. & Gonzalez-Silveira, M. (2012). Phys. Rev. Lett. 109, 165701.], 2013[Ruta, B., Baldi, G., Monaco, G. & Chushkin, Y. (2013). J. Chem. Phys. 138, 054508.]), proteins (Begam et al., 2021[Begam, N., Ragulskaya, A., Girelli, A., Rahmann, H., Chandran, S., Westermeier, F., Reiser, M., Sprung, M., Zhang, F., Gutt, C. & Schreiber, F. (2021). Phys. Rev. Lett. 126, 098001.]; Girelli et al., 2021[Girelli, A., Rahmann, H., Begam, N., Ragulskaya, A., Reiser, M., Chandran, S., Westermeier, F., Sprung, M., Zhang, F., Gutt, C. & Schreiber, F. (2021). Phys. Rev. Lett. 126, 138004.]; Ragulskaya et al., 2021[Ragulskaya, A., Begam, N., Girelli, A., Rahmann, H., Reiser, M., Westermeier, F., Sprung, M., Zhang, F., Gutt, C. & Schreiber, F. (2021). J. Phys. Chem. Lett. 12, 7085-7090.]; Reiser et al., 2022[Reiser, M., Girelli, A., Ragulskaya, A., Das, S., Berkowicz, S., Bin, M., Ladd-Parada, M., Filianina, M., Poggemann, H.-F., Begam, N., Akhundzadeh, M. S., Timmermann, S., Randolph, L., Chushkin, Y., Seydel, T., Boesenberg, U., Hallmann, J., Möller, J., Rodriguez-Fernandez, A., Rosca, R., Schaffer, R., Scholz, M., Shayduk, R., Zozulya, A., Madsen, A., Schreiber, F., Zhang, F., Perakis, F. & Gutt, C. (2022). Nat. Commun. 13, 5528.]; Chushkin et al., 2022[Chushkin, Y., Gulotta, A., Roosen-Runge, F., Pal, A., Stradner, A. & Schurtenberger, P. (2022). Phys. Rev. Lett. 129, 238001.]), magnetic systems (Shpyrko, 2014[Shpyrko, O. G. (2014). J. Synchrotron Rad. 21, 1057-1064.]; Zhang et al., 2017[Zhang, Q., Dufresne, E. M., Chen, P., Park, J., Cosgriff, M. P., Yusuf, M., Dong, Y., Fong, D. D., Zhou, H., Cai, Z., Harder, R. J., Callori, S. J., Dawber, M., Evans, P. G. & Sandy, A. R. (2017). Phys. Rev. Lett. 118, 097601.]), and clays (Bandyopadhyay et al., 2004[Bandyopadhyay, R., Liang, D., Yardimci, H., Sessoms, D. A., Borthwick, M. A., Mochrie, S. G. J., Harden, J. L. & Leheny, R. L. (2004). Phys. Rev. Lett. 93, 228302.]).

The dynamics of a system under investigation are revealed by analyzing the temporal correlations of the scattered intensity. Under coherent illumination, the resulting far-field pattern of the scattered intensity exhibits spots of constructive and destructive interference known as speckles. The dynamics of the sample lead to changes in the speckle pattern. The XPCS technique exploits the fluctuations of these speckles to extract information about the dynamic behavior of the sample. Comprehensive summaries of XPCS research studies and future prospects can be found in several reviews (e.g. Shpyrko, 2014[Shpyrko, O. G. (2014). J. Synchrotron Rad. 21, 1057-1064.]; Grübel et al., 2008[Grübel, G., Madsen, A. & Robert, A. (2008). X-ray Photon Correlation Spectroscopy (XPCS), pp. 953-995. Dordrecht: Springer Netherlands.]; Sutton, 2008[Sutton, M. (2008). C. R. Phys. 9, 657-667.]; Lehmkühler et al., 2021[Lehmkühler, F., Roseker, W. & Grübel, G. (2021). Appl. Sci. 11, 6179.]; Perakis & Gutt, 2020[Perakis, F. & Gutt, C. (2020). Phys. Chem. Chem. Phys. 22, 19443-19453.]; Sinha et al., 2014[Sinha, S. K., Jiang, Z. & Lurio, L. B. (2014). Adv. Mater. 26, 7764-7785.]). An overview on the qualitative analysis was presented by Bikondoa (2017[Bikondoa, O. (2017). J. Appl. Cryst. 50, 357-368.]). Nevertheless, a satisfying link between the theory derived for equilibrium systems (including the estimation of physical parameters of the system such as diffusion, viscosity etc.) and the analogous quantitative analysis of non-equilibrium systems, where the underlying physical parameters evolve with time, is still missing and indeed difficult to achieve.

The subject of non-equilibrium is, of course, not limited to a specific technique but increasingly relevant in various fields, for example, glass physics and mode coupling theory (Götze, 1999[Götze, W. (1999). J. Phys. Condens. Matter, 11, A1-A45.]; Martinez et al., 2010[Martinez, V. A., Bryant, G. & van Megen, W. (2010). J. Chem. Phys. 133, 114906.]), as well as growth phenomena (Headrick et al., 2019[Headrick, R. L., Ulbrandt, J. G., Myint, P., Wan, J., Li, Y., Fluerasu, A., Zhang, Y., Wiegart, L. & Ludwig, K. F. (2019). Nat. Commun. 10, 2638.]; Ju et al., 2019[Ju, G., Xu, D., Highland, M. J., Thompson, C., Zhou, H., Eastman, J. A., Fuoss, P. H., Zapol, P., Kim, H. & Stephenson, G. B. (2019). Nat. Phys. 15, 589-594.]; Dax et al., 2023[Dax, I., Zaluzhnyy, I. A., Pylypenko, A., Russegger, N., Starostin, V., Rysov, R., Westermeier, F., Sprung, M., Hinderhofer, A., Pithan, L. & Schreiber, F. (2023). New J. Phys. 25, 103033.]). Time-resolved correlations have also been pioneered in dynamic light scattering with important insights into temporal heterogeneities, higher-order correlations and spatial–temporal correlations (Cipelletti et al., 2003[Cipelletti, L., Bissig, H., Trappe, V., Ballesta, P. & Mazoyer, S. (2003). J. Phys. Condens. Matter, 15, S257-S262.]; Duri & Cipelletti, 2006[Duri, A. & Cipelletti, L. (2006). Europhys. Lett. 76, 972-978.]; Duri et al., 2005[Duri, A., Bissig, H., Trappe, V. & Cipelletti, L. (2005). Phys. Rev. E, 72, 051401.]; Cipelletti & Weitz, 1999[Cipelletti, L. & Weitz, D. A. (1999). Rev. Sci. Instrum. 70, 3214-3221.]; and references therein). The present paper attempts first to provide an overview of the conventional data analysis of XPCS studies and then to complement previous studies by discussing the quantitative analysis in the light of the connection between equilibrium and non-equilibrium, including a comparison based on specific case studies. Although the findings are limited to specific conditions and may not be readily generalized, we hope to inspire further theoretical and numerical investigations to shed light on this issue in a broader context, which is currently underrepresented in the literature.

2. Two-time correlation function

2.1. Conventional calculations of two-time correlations and their connection

Data analysis is a crucial step in XPCS, as it involves extracting the relevant information from the correlation functions of the measured time-resolved 2D speckle patterns [I(t)]. To follow non-equilibrium dynamical evolution during the XPCS measurement, it is customary to use the two-time correlation function (TTC) (Madsen et al., 2010[Madsen, A., Leheny, R. L., Guo, H., Sprung, M. & Czakkel, O. (2010). New J. Phys. 12, 055001.]; Bikondoa, 2017[Bikondoa, O. (2017). J. Appl. Cryst. 50, 357-368.]; Sutton et al., 2003[Sutton, M., Laaziri, K., Livet, F. & Bley, F. (2003). Opt. Express, 11, 2268.]):

[{Corr}({\bf q},t_{1},t_{2}) = {{\langle I({\bf q},t_{1})I({\bf q},t_{2}) \rangle} \over {\langle I({\bf q},t_{1})\rangle\langle I({\bf q},t_{2})\rangle}}. \eqno (1)]

This `Corr-TTC' calculates the correlation between intensities at times t1 and t2 averaged over all pixels at the same q ring. q is the scattering vector [q = |q| = (4π/λ)sinθ, where θ is half the scattering angle and λ is the wavelength of the incident radiation]. For simplicity, we shall assume an isotropic sample. Here, the normalization is performed by the mean intensity.

Another possibility for the calculation of the TTC is G-TTC – the autocovariance of the intensity normalized by its standard deviation σ (Bikondoa, 2017[Bikondoa, O. (2017). J. Appl. Cryst. 50, 357-368.]; Brown et al., 1997[Brown, G., Rikvold, P. A., Sutton, M. & Grant, M. (1997). Phys. Rev. E, 56, 6601-6612.]):

[\eqalignno{ G({\bf q},t_{1},t_{2}) & = {{\overline{I(t_{1})I(t_{2})}-\overline{I(t_{1})} \ \overline{I(t_{2})}} \over {[\overline{I^{2}(t_{1})}-\overline{I(t_{1})}^{2}]^{ {{1} / {2}}}\,[\overline{I^{2}(t_{2})}-\overline{I(t_{2})}^{2}]^{{{1}/ {2 }}}}} \cr & = {{\overline{I(t_{1})I(t_{2})}-\overline{I(t_{1})}\ \overline{I(t_{2 })}} \over {\sigma(t_{1})\ \sigma(t_{2})}},& (2)}]

where [\overline{I} = \langle I\rangle].

If the scattered intensity fluctuates around a stable mean and has a negative exponential distribution (i.e. fully coherent) (Goodman, 2020[Goodman, J. W. (2020). Speckle Phenomena in Optics: Theory and Applications, 2nd ed. Bellingham: SPIE Press.]; Pusey & Van Megen, 1989[Pusey, P. & Van Megen, W. (1989). Physica A, 157, 705-741.]), it can be shown that the average speckle intensity equals the standard deviation of speckle intensities (Brown et al., 1997[Brown, G., Rikvold, P. A., Sutton, M. & Grant, M. (1997). Phys. Rev. E, 56, 6601-6612.]; Loudon, 1983[Loudon, R. (1983). The Quantum Theory of Light. Oxford University Press.]):

[\overline{I} = \sigma = \left( {\overline{I^{2}}-\overline{I}^{2}} \right)^{1/2}.\eqno (3)]

Thus, G(q, t1, t2) = Corr(q, t1, t2) − 1, and the use of G-TTC and Corr-TTC is physically equivalent under these conditions. In the more general case of partial coherence and non-stable mean intensity, the normalized standard deviation β = σ2/〈I2 = βsourceβsample is a product of speckle contrast due to the properties of the X-ray source/experimental setup (Möller et al., 2021[Möller, J., Reiser, M., Hallmann, J., Boesenberg, U., Zozulya, A., Rahmann, H., Becker, A.-L., Westermeier, F., Zinn, T., Sprung, M., Narayanan, T., Gutt, C. & Madsen, A. (2021). New J. Phys. 23, 093041.]) and fluctuations from the non-constant mean scattering intensity from the sample (Goodman, 1985[Goodman, J. W. (1985). Statistical Optics. New York: Wiley-Interscience.]; Pusey & Jakeman, 1975[Pusey, P. N. & Jakeman, E. (1975). J. Phys. A Math. Gen. 8, 392-410.]; Dainty, 1977[Dainty, J. (1977). Prog. Opt. 14, 1-46.]). Using this relation, equation (2)[link] can be rewritten as follows:

[G = {{\overline{I(t_{1})I(t_{2})}-\overline{I(t_{1})}\ \overline{I(t_{2}) }} \over {{ {\beta} ^{1/2}\ \overline{I(t_{1})}}\ {\beta} ^{1/2}\ \overline{I(t_{ 2})}}} = {{1} \over {\beta}}\ {{\overline{I(t_{1})I(t_{2})}} \over {\overline{I(t_{1})}\ \overline{I(t_{2})}}}-{{1} \over {\beta}} = {{1} \over {\beta}}({ Corr}-1).\eqno (4)]

This equation gives the general relation between the two conventional ways of calculating the TTC. A list of functions used in the manuscript, as well as their relations (including the general scheme for XPCS data analysis), can be found in Section S2 in the supporting information. In the case of low-intensity statistics, while equation (4)[link] with its definition of contrast and the definitions of Corr-TTC and G-TTC are mathematically true, they are not operational. In such instances, the contrast needs to be determined from the binomial distribution, and alternative data analysis methods, as well as experimental procedures, are employed instead of TTCs to obtain the intermediate scattering function (Roseker et al., 2018[Roseker, W., Hruszkewycz, S. O., Lehmkühler, F., Walther, M., Schulte-Schrepping, H., Lee, S., Osaka, T., Strüder, L., Hartmann, R., Sikorski, M., Song, S., Robert, A., Fuoss, P. H., Sutton, M., Stephenson, G. B. & Grübel, G. (2018). Nat Commun, 9, 1704.]; Hua et al., 2020[Hua, N., Zaluzhnyy, I. A., Hrkac, S. B., Shabalin, A. G. & Shpyrko, O. G. (2020). J. Synchrotron Rad. 27, 1626-1632.]; Hruszkewycz et al., 2012[Hruszkewycz, S. O., Sutton, M., Fuoss, P. H., Adams, B., Rosenkranz, S., Ludwig, K. F., Roseker, W., Fritz, D., Cammarata, M., Zhu, D., Lee, S., Lemke, H., Gutt, C., Robert, A., Grübel, G. & Stephenson, G. B. (2012). Phys. Rev. Lett. 109, 185502.]). These approaches are beyond the scope of the current article.

2.2. TTC for a non-equilibrium system with evolution of the intensity distribution.

Classical data analysis assumes the extraction of the one-time correlation function, g2, from the TTC by different coordinate systems (see Section 3[link]) (Bikondoa, 2017[Bikondoa, O. (2017). J. Appl. Cryst. 50, 357-368.]). For the correlation map defined by equation (1)[link], the g2 function may be determined via the generalized Siegert relation:

[g_{2_{ Corr}}({\bf q},\Delta t) = 1+\beta({\bf q})|g_{1}(t,\Delta t)|^{2}.\eqno (5)]

Depending on the sample dynamics, the g1 function may have different forms. Nevertheless, the standard approach which is valid for the majority of systems is a Kohlrausch–Williams–Watts (KWW) relation: [g_{1} = \exp\left[-(\Delta t/\tau)^{\gamma}\right]], with the KWW exponent γ and the relaxation time of the system τ (Williams & Watts, 1970[Williams, G. & Watts, D. C. (1970). Trans. Faraday Soc. 66, 80-85.]) and Δt = |t2t1|. Substituting equation (5)[link] into equation (4)[link], it is possible to obtain the g2 function for the correlation map G-TTC:

[\eqalignno{ g_{2_{G}}({\bf q},\Delta t) & = {{1} \over {\beta}}[g_{2_{ Corr}}({\bf q},\Delta t)-1] = {{1} \over {\beta}}\left[1+\beta({\bf q})|g_{1}(t,\Delta t)|^{2}-1\right] \cr & = |g_{1}(t,\Delta t)|^{2} = \exp\left[-2(\Delta t/\tau)^{\gamma} \right].&(6)}]

Thus, we can conclude that if under the investigated conditions the Siegert relation is applicable [equation (5[link])], g2G does not depend on the contrast β and is a function of the parameters of dynamics of the system only. This feature will be used in Section 3[link]. While g2G continues to capture all physical parameters accessible by g2Corr, the normalization already performed by the contrast should be taken into account when extracting the parameters of dynamics that depend on it (e.g. the non-ergodicity parameter).

The above derivations lead to an important point. The dynamics of a system in equilibrium can be described using the fluctuations of the static structure factor, where the average density is constant and the density fluctuates around this mean value with time. In this equilibrium case, [\overline{I(q)}] is constant and the dynamics can be directly extracted. Thus, there is no physical difference between the use of G-TTC or Corr-TTC.

If the system is non-equilibrium and the fluctuations from the sample cannot be described by zero-mean Gaussian statistics, the contrast β may evolve with time. Consider a system that exhibits not only dynamics but also kinetic evolution, e.g. a system undergoing phase separation. In contrast to equilibrium systems, there is a change in the behavior of the mean intensity [\overline{I(q)}] due to its kinetic evolution. This change may lead to a variation of β. Therefore, for kinetically evolving systems, the calculation of TTCs with Corr [equation (1[link])] followed by the extraction of g2Corr [equation (5[link])] may lead to ambiguous data analysis, while the correlation map G-TTC [equation (4[link])] followed by extraction of g2G [equation (6[link])] is not influenced by the contrast evolution.

In order to illustrate the effects for a kinetically evolving system, Fig. 1[link] shows the example of simulations based on the Cahn–Hilliard equation for spinodal decomposition (Ragul­skaya et al., 2022[Ragulskaya, A., Starostin, V., Begam, N., Girelli, A., Rahmann, H., Reiser, M., Westermeier, F., Sprung, M., Zhang, F., Gutt, C. & Schreiber, F. (2022). IUCrJ, 9, 439-448.]; Girelli et al., 2021[Girelli, A., Rahmann, H., Begam, N., Ragulskaya, A., Reiser, M., Chandran, S., Westermeier, F., Sprung, M., Zhang, F., Gutt, C. & Schreiber, F. (2021). Phys. Rev. Lett. 126, 138004.]; Cahn & Hilliard, 1958[Cahn, J. W. & Hilliard, J. E. (1958). J. Chem. Phys. 28, 258-267.], 1959[Cahn, J. W. & Hilliard, J. E. (1959). J. Chem. Phys. 31, 688-699.]) under fully coherent light (βsource = 1). The system evolves kinetically, resulting in a significant variation of the mean intensity [Fig. 1[link](c)]. G-TTC is stable, while Corr-TTC shows fluctuations along the diagonal [compare Fig. 1[link](a) and Fig. 1[link](b)] with values that sometimes exceed 2 owing to the limited number of scatterers (Dainty, 1977[Dainty, J. (1977). Prog. Opt. 14, 1-46.]; Pusey & Jakeman, 1975[Pusey, P. N. & Jakeman, E. (1975). J. Phys. A Math. Gen. 8, 392-410.]). These fluctuations are the same as the contrast of the system, β = σ2/〈I2 = βsourceβsample = βsample, and reflect the change in the scattering intensity distribution [Fig. 1[link](d)].

[Figure 1]
Figure 1
Examples of Corr-TTC (a) and G-TTC (b) based on Cahn–Hilliard simulations for the same q value. The different options for the g2 cuts are also illustrated in (b): the CCS (purple dashed line) represents the correlation between the present (at time t0) and future dynamics, the ACS (gray dashed line) that between the past and future. The black dashed line represents another possibility – the correlation between the present and past. (c) The evolution of the mean intensity [\overline{I} = \langle I\rangle] of the system. (d) A comparison between β = σ2/〈I2 (purple), [g_{2_{Corr}}(t_{1} = t_{2})-1] (orange) and [g_{2_{G}}(t_{1} = t_{2})] (green). G-TTC is less sensitive to the kinetic changes of the system than Corr-TTC.

3. Analysis of TTC functions using different time coordinate systems

3.1. Quantitative analysis of TTCs

Quantitative description of the evolution of the correlation function typically requires slicing (`cutting') the TTC at distinct observation times t0 and extracting the parameters of the dynamics of the system, such as relaxation time, τ(t0), and Kohlrausch–Williams–Watt exponent, γ(t0).

The evolution of the relaxation time of a non-equilibrium system is frequently used for the calculation of the evolution of the diffusion coefficient, velocity or other macroscopic observables of the system (Chushkin et al., 2022[Chushkin, Y., Gulotta, A., Roosen-Runge, F., Pal, A., Stradner, A. & Schurtenberger, P. (2022). Phys. Rev. Lett. 129, 238001.]; Reiser et al., 2022[Reiser, M., Girelli, A., Ragulskaya, A., Das, S., Berkowicz, S., Bin, M., Ladd-Parada, M., Filianina, M., Poggemann, H.-F., Begam, N., Akhundzadeh, M. S., Timmermann, S., Randolph, L., Chushkin, Y., Seydel, T., Boesenberg, U., Hallmann, J., Möller, J., Rodriguez-Fernandez, A., Rosca, R., Schaffer, R., Scholz, M., Shayduk, R., Zozulya, A., Madsen, A., Schreiber, F., Zhang, F., Perakis, F. & Gutt, C. (2022). Nat. Commun. 13, 5528.]; Czakkel & Madsen, 2011[Czakkel, O. & Madsen, A. (2011). Europhys. Lett. 95, 28001.]; Lehmkühler et al., 2020[Lehmkühler, F., Dallari, F., Jain, A., Sikorski, M., Möller, J., Frenzel, L., Lokteva, I., Mills, G., Walther, M., Sinn, H., Schulz, F., Dartsch, M., Markmann, V., Bean, R., Kim, Y., Vagovic, P., Madsen, A., Mancuso, A. P. & Grübel, G. (2020). Proc. Natl Acad. Sci. USA, 117, 24110-24116.]). Through this approach, a time series of cuts and corresponding parameters [τ(t0), γ(t0)] obtained from the TTC results in a time-resolved evolution of macroscopic observables of the system under study.

This quantitative analysis is based on derivations made for equilibrium systems. Therefore, such treatment of the experimental data requires that τ(t0) represents an `instantaneous' description of the system, i.e. it is obtained for a specific moment in time t0 (Ladd et al., 1995[Ladd, A. J. C., Gang, H., Zhu, J. X. & Weitz, D. A. (1995). Phys. Rev. Lett. 74, 318-321.]). We will refer to such dynamics as effective dynamics. For instance, τ(t0) is considered not to be affected by any possible future (t > t0) perturbations of the system (e.g. beam damage) (Ruta et al., 2017[Ruta, B., Zontone, F., Chushkin, Y., Baldi, G., Pintori, G., Monaco, G., Rufflé, B. & Kob, W. (2017). Sci. Rep. 7, 3962.]; Reiser et al., 2022[Reiser, M., Girelli, A., Ragulskaya, A., Das, S., Berkowicz, S., Bin, M., Ladd-Parada, M., Filianina, M., Poggemann, H.-F., Begam, N., Akhundzadeh, M. S., Timmermann, S., Randolph, L., Chushkin, Y., Seydel, T., Boesenberg, U., Hallmann, J., Möller, J., Rodriguez-Fernandez, A., Rosca, R., Schaffer, R., Scholz, M., Shayduk, R., Zozulya, A., Madsen, A., Schreiber, F., Zhang, F., Perakis, F. & Gutt, C. (2022). Nat. Commun. 13, 5528.]; Chushkin et al., 2022[Chushkin, Y., Gulotta, A., Roosen-Runge, F., Pal, A., Stradner, A. & Schurtenberger, P. (2022). Phys. Rev. Lett. 129, 238001.]; Timmermann et al., 2023[Timmermann, S., Anthuparambil, N. D., Girelli, A., Begam, N., Kowalski, M., Retzbach, S., Senft, M. D., Akhundzadeh, M. S., Pog­gemann, H., Moron, M., Hiremath, A., Gutmüller, D., Dargasz, M., Öztürk, Ö., Paulus, M., Westermeier, F., Sprung, M., Ragulskaya, A., Zhang, F., Schreiber, F. & Gutt, C. (2023). Sci. Rep. 13, 11048.]). This can be interpreted as the relaxation time of the corresponding equilibrium system, wherein the macro-parameters align with the investigated moment of the non-equilibrium system.

To this end, one needs to estimate the relaxation time of the corresponding equilibrium system based on the TTC from a non-equilibrium system. This is typically done by extracting g2 functions in the form of one-dimensional cuts. There are several ways to obtain cuts from TTCs, which are discussed below.

3.2. Frequently employed coordinate systems

The most frequently employed methods for extracting g2 functions are horizontal and diagonal cuts, as illustrated in Fig. 1[link](b).

Diagonal cuts were introduced alongside TTCs by Brown et al. (1997[Brown, G., Rikvold, P. A., Sutton, M. & Grant, M. (1997). Phys. Rev. E, 56, 6601-6612.]). In the context of XPCS, diagonal cuts have a longstanding history and were commonly used in the past (Malik et al., 1998[Malik, A., Sandy, A. R., Lurio, L. B., Stephenson, G. B., Mochrie, S. G. J., McNulty, I. & Sutton, M. (1998). Phys. Rev. Lett. 81, 5832-5835.]; Brown et al., 1999[Brown, G., Rikvold, P. A., Sutton, M. & Grant, M. (1999). Phys. Rev. E, 60, 5151-5162.]; Livet et al., 2001[Livet, F., Bley, F., Caudron, R., Geissler, E., Abernathy, D., Detlefs, C., Grübel, G. & Sutton, M. (2001). Phys. Rev. E, 63, 036108.]; Sutton et al., 2003[Sutton, M., Laaziri, K., Livet, F. & Bley, F. (2003). Opt. Express, 11, 2268.]; Fluerasu et al., 2005[Fluerasu, A., Sutton, M. & Dufresne, E. M. (2005). Phys. Rev. Lett. 94, 055501.]; Müller et al., 2011[Müller, L., Waldorf, M., Gutt, C., Grübel, G., Madsen, A., Finlayson, T. R. & Klemradt, U. (2011). Phys. Rev. Lett. 107, 105701.]; Orsi et al., 2010[Orsi, D., Cristofolini, L., Fontana, M. P., Pontecorvo, E., Caronna, C., Fluerasu, A., Zontone, F. & Madsen, A. (2010). Phys. Rev. E, 82, 031804.]; Bikondoa et al., 2012[Bikondoa, O., Carbone, D., Chamard, V. & Metzger, T. H. (2012). J. Phys. Condens. Matter, 24, 445006.]; Ruta et al., 2012[Ruta, B., Chushkin, Y., Monaco, G., Cipelletti, L., Pineda, E., Bruna, P., Giordano, V. M. & Gonzalez-Silveira, M. (2012). Phys. Rev. Lett. 109, 165701.]). In 2017, Bikondoa introduced horizontal cuts as an alternative to diagonal cuts (Bikondoa, 2017[Bikondoa, O. (2017). J. Appl. Cryst. 50, 357-368.]). Since the horizontal cuts were argued to be more intuitive and consistent with the standard calculations in statistical mechanics, they were termed a `conventional coordinate system' (CCS). Moreover, it was discussed that the use of diagonal cuts could lead to interpretation issues when external forces or perturbations are present [e.g. such as the system presented by Ruta et al. (2017[Ruta, B., Zontone, F., Chushkin, Y., Baldi, G., Pintori, G., Monaco, G., Rufflé, B. & Kob, W. (2017). Sci. Rep. 7, 3962.])], since it might be argued that they mix events prior and subsequent to the perturbation. Consequently, diagonal cuts were referred to as the `alternative coordinate system' (ACS) according to Bikondoa (2017[Bikondoa, O. (2017). J. Appl. Cryst. 50, 357-368.]). As a result, the CCS approach has gained popularity and is increasingly used (Zhang et al., 2021[Zhang, Q., Dufresne, E. M., Nakaye, Y., Jemian, P. R., Sakumura, T., Sakuma, Y., Ferrara, J. D., Maj, P., Hassan, A., Bahadur, D., Ramakrishnan, S., Khan, F., Veseli, S., Sandy, A. R., Schwarz, N. & Narayanan, S. (2021). J. Synchrotron Rad. 28, 259-265.]; Girelli et al., 2021[Girelli, A., Rahmann, H., Begam, N., Ragulskaya, A., Reiser, M., Chandran, S., Westermeier, F., Sprung, M., Zhang, F., Gutt, C. & Schreiber, F. (2021). Phys. Rev. Lett. 126, 138004.]; Lehmkühler et al., 2021[Lehmkühler, F., Roseker, W. & Grübel, G. (2021). Appl. Sci. 11, 6179.]). Nonetheless, the traditional ACS approach is still being used today, inter alia, for facilitating comparison with earlier XPCS studies.

As the name suggests, diagonal cuts [g_{2_{\rm ACS}}(t_{0},\Delta t)] at different sample ages t0 = tage = (t1 + t2)/2 are obtained by taking the line perpendicular to the t1 = t2 diagonal with the delay time Δt = |t2t1|. Horizontal cuts at different waiting times t0 = tw are performed by extraction of lines with t1 = constant (or t2 = constant) and the delay time defined as Δt := |t2t1|. In this manner, diagonal cuts (ACS) represent correlations between the past and the future states of the system relative to the considered moment t0, whereas horizontal cuts represent correlations of the system at moment t0 with its future states. Therefore, in fact, the horizontal cuts may also mix events prior and subsequent to the perturbation if the latter happens in the future. We note that, for the complete description, one may also consider vertical cuts [see Fig. 1[link](b)], i.e. correlations of the system at moment t0 with its past states. In the following, we mainly focus on the ACS and CCS, but this third option can be straightforwardly derived similarly to the CCS.

Importantly, for non-equilibrium systems, the ACS and CCS can produce different evolutions of τ(t0) and γ(t0). These values play a vital role in interpreting the results, and as such, their correct determination is crucial for characterizing the dynamics of the system.

3.3. Connection between ACS and CCS for equilibrium and non-equilibrium systems

First, we discuss how the ACS and CCS correlation functions are connected. They can be defined as

[g_{1_{\rm CCS}}(t,\Delta t) = \langle E^{*}(t)E(t+\Delta t)\rangle_{N}\eqno (7)]

and

[g_{1_{\rm ACS}}(t,\Delta t) = \langle E^{*}(t-\Delta t/2)E(t+\Delta t/2)\rangle_{N},\eqno (8)]

where 〈〉N is the average over the ensemble, E is the electric field and the asterisk represents complex conjugation.

We note that, for each definition of the correlation function, the Siegert relation holds [equation (5[link])]:

[g_{2_{Corr,{\rm CCS}}}(t,\Delta t) = 1+\beta|g_{1_{\rm CCS}}(t,\Delta t)|^{2},\eqno (9)]

[g_{2_{Corr,{\rm ACS}}}(t,\Delta t) = 1+\beta|g_{1_{\rm ACS}}(t,\Delta t)|^{2}.\eqno (10)]

In the case of equilibrium systems, correlation functions are translation invariant in time, i.e.

[g_{1_{\rm CCS}}(t,\Delta t) = g_{1_{\rm CCS}}(t+T,\Delta t) = g_{1_{\rm CCS}}(0,\Delta t),\eqno (11)]

[g_{1_{\rm ACS}}(t,\Delta t) = g_{1_{\rm ACS}}(t+T,\Delta t) = g_{1_{\rm ACS}}(0,\Delta t),\eqno (12)]

so they only depend on Δt. Therefore, [g_{1_{\rm CCS}}] and [g_{1_{\rm ACS}}] are equivalent in the case of equilibrium systems, since

[g_{1_{\rm ACS}}(t,\Delta t) = g_{1_{\rm ACS}}(t+\Delta t/2,\Delta t) = g_{1_{\rm CCS}}(t, \Delta t).\eqno (13)]

For non-equilibrium systems, this is not the case, as the first equality in equation (13)[link] generally does not hold. However, the second equality still holds by definition, so that for non-equilibrium we expect generally

[g_{1_{\rm ACS}}(t,\Delta t)\neq g_{1_{\rm ACS}}(t+\Delta t/2,\Delta t) = g_{1_{\rm CCS}}(t, \Delta t).\eqno (14)]

This relation only reveals a trivial transformation of the variables. Nevertheless, it shows how [g_{1_{\rm CCS}}] and [g_{1_{\rm ACS}}] are related in non-equilibrium systems and that, in general, they are not identical.

3.4. Geometrical illustration of ACS and CCS

The disparity between the ACS and CCS extends beyond mere mathematical expressions. Instead, it lies in their fundamental approach to correlation. The CCS correlates the system with itself in the future, whereas the ACS reveals the correlation between the system in the future and itself in the past, with equal distance from the time under investigation. These approaches are illustrated in Fig. 2[link]. As was discussed before, it is assumed that the correlation function 〈 E*(t′)E(t′′)〉N between times t′ and t′′ is equal to the correlation function of the corresponding equilibrium system with some fixed dynamical property τinstantaneous = τ(tcen). Furthermore, it is natural to assume that the time moment tcen is the average tcen = (t′ + t′′)/2 (see Fig. 2[link]). If this is the case, at any moment t0 in time, we shall consider the ACS definition of the correlation function (diagonal cut) as the one that more closely resembles the effective dynamics of the system τ(t0), while the horizontal cut will correspond to τ(t0 + Δt/2). Therefore, the g2G functions for the ACS and CCS can be represented as follows:

[g_{2_{G,{\rm ACS}}}(t_{0},\Delta t) = \exp\left\{-2\left[{{\Delta t} \over {\tau\left(t_{0 }\right)}}\right]^{\gamma\left(t_{0}\right)}\right\},\eqno (15)]

[g_{2_{G,{\rm CCS}}}(t_{0},\Delta t) = \exp\left\{-2\left[{{\Delta t} \over {\tau\left(t_{0 }+{{\Delta t} /{2}}\right)}}\right]^{\gamma\left(t_{0}+{{\Delta t}/ {2}} \right)}\right\}.\eqno (16)]

[Figure 2]
Figure 2
Schematic description of obtaining dynamics in a non-equilibrium system through the ACS (left panel) and CCS (right panel) in the (t, Δt) plane [t′ = t′(t, Δt), t′′ = t′′(t, Δt)]. The time dependency is illustrated by different colors, gradually changing from red to deep blue. For ACS: t′ = t − Δt/2 and t′′ = t + Δt/2. Therefore, at any time moment t0, tcen = (t′ + t′′)/2 = t0 and, thus, the calculated τ can be approximated by the effective dynamics with τ(t0) (see examples for t01 and t02). For CCS: t′ = t and t′′ = t + Δt. Therefore, at any time moment t0, tcen = (t′ + t′′)/2 = t0 + Δt/2. Thus, the relaxation time τ = τ(t0 + Δt/2) cannot be approximated by effective dynamics anymore, except when the system is in quasi-equilibrium and τ(t0 + Δt/2) ≈ τ(t0).

Obviously, for an equilibrium process, τ(t0) = τ(t0 + Δt/2) = constant and γ(t0) = γ(t0 + Δt/2) = constant, which leads to the same results for the ACS and CCS. However, in the case of a non-equilibrium process, only the ACS has the conventional form of the g2 function, which is typically used for data analysis for the extraction of the relaxation time τ(t0): [g_{2_{G,{\rm ACS}}}(t_{0},\Delta t) = g_{2_{G}}(t_{0},\Delta t)] [compare equations (15[link]) and (6[link])]. In contrast, for each Δt, [g_{2_{\rm CCS}}] represents a different set of τ(t0 + Δt/2) and γ(t0 + Δt/2) pairs. Therefore, while the ACS cuts may be used to describe the effective dynamics, followed by extraction of the momentary properties of the system, the CCS cuts provide some averaged properties of the correlation between the considered moment and the future. These fundamental differences should be taken into account for the interpretation of the results of the XPCS data analysis.

We note that, while the ACS may correspond to the dynamics of the corresponding equilibrium of the non-equilibrium system, that is not necessarily guaranteed. In the following, we demonstrate the system for which this assumption holds true.

3.5. Simulation example

In this section, we demonstrate the use case of ACS cuts for the extraction of effective dynamics. We employ the simulations for a model example of a non-equilibrium system and compare the ACS and CCS with the corresponding equilibrium for each sample age.

The model system consists of a set of particles. The position of a particle is represented by the vector [{\bf r} = (x,y)], for which the movement along the Cartesian coordinates x and y is statistically independent. The probability density function [P({\b r},t) = P(x,t)P(y,t)] is set to

[\eqalign{ P(x,t) & = {{1} \over {\left[ {4\pi D(t)t} \right]^{1/2}}}\exp\left[-{{\left(x-x_{0} \right)^{2}} \over {4D(t)t}}\right], \cr P(y,t) & = {{1} \over {\left[ {4\pi D(t)t} \right]^{1/2}}}\exp\left[-{{\left(y-y_{0} \right)^{2}} \over {4D(t)t}}\right],}\eqno (17)]

where t is time and x0 and y0 are the initial positions of the particle. At each discrete time (sample age) of a simulated observation, we modify the D(t) parameter of the system to simulate non-equilibrium dynamics. This approach allows us to obtain both the TTC of the non-equilibrium system and the corresponding equilibrium parameters for each discrete moment t0. Parameters of the simulations can be found in Section S1.

While, in the general case, such a non-equilibrium system does not correspond to any known example of anomalous diffusion, and the corresponding equation of motion remains unclear, for each specific t0, equation (17)[link] describes `momentarily' Brownian motion, which serves as a corresponding equilibrium scenario of this non-equilibrium system. In this case, D(t0) is the diffusion coefficient of this corresponding equilibrium and the relaxation time τ(t0)cor.eq. ∝ 1/D(t0). Furthermore, γ(t0)cor.eq. is equal to 1 and Γ(t0)cor.eq. = 1/τ(t0) ∝ q2. The values of the τ(t0)cor.eq. and γ(t0)cor.eq. parameters for the equilibrium were double-checked by performing a simulation with constant D = D(t0).

For the sake of illustration, we demonstrate three cases of the evolution of the 1/D = 1/D(t) parameter: linear (Case 1), exponential (Case 2) and sinusoidal (Case 3). In Case 1, we assume that 1/D of the model system changes linearly for each discrete time step j:

[1/D(j+1) = {\rm constant} +1/D(j),\quad j = [1,N_{\rm t}]\rm {,} \eqno (19)]

where Nt is the total number of time steps of the simulation. In this case, a linear change in the relaxation time of the corresponding equilibrium scenario τ(t0)cor.eq. = τ(j)cor.eq. ∝ 1/D(j) with time is expected. Depending on the sign of the constant in equation (18)[link], the system slows down (Case 1a) or accelerates (Case 1b).

The results of the Case 1a simulation of a linear slowdown (constant > 0) are presented in Fig. 3[link]. Both the CCS and ACS qualitatively capture the linear behavior of the relaxation time as well as the Γ(t) = 1/τ(t) ∝ q2 behavior. Furthermore, in this case, the ACS manages to capture the momentary description and, thus, the effective dynamics. Remarkably, the relaxation time and γ have values that are similar to the corresponding equilibrium equivalents. On the other hand, the CCS does not provide a momentary description of the system. Instead, the obtained relaxation time for all q values is larger than that for the corresponding equilibrium scenarios and also suggests a different slope [see Figs. 3[link](b) and 3[link](c)]. Furthermore, γ is around 0.75 [Fig. 3[link](d)] in contrast to 1 for the equilibrium system. The CCS results come as a direct reflection of the slowdown of the system with time. As was discussed earlier, [g_{2_{G,{\rm CCS}}}] for each Δt represents a different set of τ(t0 + Δt/2) and γ(t0 + Δt/2) pairs [see schematic representation in Fig. 4[link](a) and equation (16)[link]]. For our model system, for any positive Δt, τ(t0) < τ(t0 + Δt/2) and γ(t)cor.eq. = 1. Therefore, [g_{2_{G,{\rm CCS}}}] is stretched in comparison to [g_{2_{G,{\rm ACS}}}] [see equations (16[link]) and (15[link])]. Similarly, if this system is accelerated (Case 1b), [g_{2_{G,{\rm ACS}}}] will catch the effective dynamics, while [g_{2_{G,{\rm CCS}}}] will be compressed in comparison with [g_{2_{G,{\rm ACS}}}]. This conclusion is demonstrated in the simulation of Case 1b in supplementary Fig. S2.

[Figure 3]
Figure 3
Data analysis for a model system for Case 1a – a linear increase in 1/D. (a) G-TTC for q = 74 pixels. (b) Relaxation rate Γ as a function of q2 at time = 25. (c) and (d) represent relaxation time τ and γ as functions of time, correspondingly. Orange circles display results for CCS analysis and light-blue squares those for the ACS. The dashed blue line shows results from corresponding equilibrium systems. Results for other q and time values are similar.
[Figure 4]
Figure 4
(a) Schematic representation of [g_{2_{G,{\rm CCS}}}], which follows equation (16)[link] and Fig. 2[link]. For each Δt, [g_{2_{G,{\rm CCS}}}] represents a different set of τ(t0 + Δt/2) and γ(t0 + Δt/2) pairs. (b) Comparison of g2 cuts for the simulated system for Case 1a, presented in Fig. 3[link], for q = 74 pixels and time = 50. The CCS and ACS cuts from the simulated TTC are represented with orange circles and blue squares, respectively. The black dashed line shows [g_{2_{G,{\rm CCS}}}], calculated via equation (16)[link]. The blue dashed line shows the conventional g2G function, calculated via equation (2)[link].

This conclusion can be further supported by the comparison of the ACS and CCS g2 cuts of our model system with the conventional g2 function [equation (6[link])] and [g_{2_{\rm CCS}}] [equation (16[link])], both calculated from the instantaneous τ, obtained from the corresponding equilibrium equivalents and presented in Fig. 4[link](b). The CCS g2 cuts align well with the estimation via equation (16)[link] and the ACS g2 cuts overlap with the conventional representation [equation (6[link])].

These general conclusions also remain valid for a possible nonlinear evolution of 1/D. Fig. 5[link] (Case 2) and Fig. 6[link] (Case 3) demonstrate exponential and sinusoidal changes in 1/D of the investigated system, correspondingly. In both cases, the ACS data analysis results overlap with the corresponding equilibrium scenario and capture the effective dynamics, while the CCS results represent some averaged dynamics and reflect the acceleration or slowdown of the system. For example, γ and τ in Case 3 correspond to γCCS < γcor.eq. = 1 and τCCS > τcor.eq. in the accent of the sinus, while γCCS > 1 and τCCS < τcor.eq. for the descent part of the sinus, which is consistent with the previous discussions.

[Figure 5]
Figure 5
Data analysis for a model system for Case 2 – an exponential increase in 1/D [to be compared with Fig. 3[link] (linear increase)]. (a) G-TTC for q = 86 pixels. (b) Relaxation rate Γ as a function of q2 at time = 70. (c) and (d) represent relaxation time τ and γ as functions of time, correspondingly. Orange circles display results for CCS analysis and light-blue squares those for the ACS. The dashed blue line shows theoretical behavior. Results for other q and time values are similar.
[Figure 6]
Figure 6
Data analysis for a model system for Case 3 – a sinusoidal change in 1/D [to be compared with Figs. 3[link] and 5[link] (linear and exponential cases)]. (a) G-TTC for q = 83 pixels. (b) Relaxation rate Γ as a function of q2 at time = 25. (c) and (d) represent relaxation time τ and γ as functions of time, correspondingly. Orange circles display results for CCS analysis and light-blue squares those for the ACS. The dashed blue line shows results from corresponding equilibrium systems. Results for other q and time values are similar.

Therefore, the ACS and CCS analyses of the system give generally different descriptions with different quantities of key parameters, which leads to a different qualitative interpretation of the dynamics. The ACS analysis of various simulations of the model system suggests that the system under investigation can be described by Brownian motion momentarily at each measurement time, capturing the effective dynamics. The CCS analysis suggests subdiffusion or superdiffusion behavior with a relaxation time evolution distinct from the corresponding equilibrium scenario. The observed discrepancies arise as a natural outcome of the functional form of the corresponding g2 functions. Nevertheless, the prevailing XPCS data analysis does not encompass an assessment of the distinctions between the two types of cuts and their resulting quantitative effects on the relaxation time and KWW parameters, which may lead to ambiguous interpretation.

4. Summary and conclusions

In this article, we discussed the analysis of XPCS data, focusing specifically on the comparison of two widely used TTC calculation methods: normalization by the mean (Corr-TTC) and normalization by the standard deviation (G-TTC). We demonstrate that for kinetically evolving systems Corr-TTC is susceptible to intensity variations, potentially leading to inconclusive data interpretation. In contrast, G-TTC is generally robust against these fluctuations. Therefore, we recommend using G-TTC for analyzing processes such as film growth, coarsening, phase separation etc.

We then compared the two widely used methods for extracting one-time correlation and relevant information from the TTC: the ACS introduced by Brown et al. (1997[Brown, G., Rikvold, P. A., Sutton, M. & Grant, M. (1997). Phys. Rev. E, 56, 6601-6612.]) and the CCS recommended by Bikondoa (2017[Bikondoa, O. (2017). J. Appl. Cryst. 50, 357-368.]). While for equilibrium systems these methods produce consistent results, for non-equilibrium systems the ACS and CCS can produce distinct evolutions of the relaxation time (τ) and Kohlrausch–Williams–Watt exponent (γ), which are crucial for interpreting experimental data.

On the basis of a geometrical representation, we derived the functional forms of g2 for ACS and CCS cuts. Our analysis revealed that for non-equilibrium systems the ACS yields the conventional (theoretical) form of the g2 function, which is typically used for data analysis to extract the relaxation time τ(t0): [g_{2_{G,{\rm ACS}}}(t_{0},\Delta t) = g_{2_{G}}(t_{0},\Delta t)]. In contrast, [g_{2_{\rm CCS}}] for each Δt corresponds to a distinct set of τ(t0 + Δt/2) and γ(t0 + Δt/2) pairs. We demonstrated these dependencies through simulations of case studies. Therefore, while ACS cuts may be used to extract the effective properties of the system, the CCS reflects the average properties of correlations between the considered time point and the future. These distinctions between the two analysis methods and their consequential quantitative impacts on the relaxation time and KWW parameters should be carefully considered when interpreting XPCS data. For instance, if the experimental results are compared with theoretical predictions, it is important to first identify the analysis method that aligns with the assumptions of the theoretical model used. As we are lacking general analytical expressions for non-equilibrium two-time correlation functions, the TTC analysis would greatly benefit from guidance through further simulation work.

5. Related literature

The following additional reference is cited in the supporting information: Barton et al. (1998[Barton, B. F., Graham, P. D. & McHugh, A. J. (1998). Macromolecules, 31, 1672-1679.]).

Supporting information


Acknowledgements

The authors would like to thank Martin Oettel for productive discussions. Open access funding enabled and organized by Projekt DEAL.

Conflict of interest

There are no conflicts of interest to declare.

Funding information

Funded by EXC number 2064/1 – Project number 390727645. This work was supported by the DFG and BMBF (05K19PS1, 05K20PSA, 05K22PS1, 05K20VTA, 05K22VTA).

References

First citationAngelini, R., Madsen, A., Fluerasu, A., Ruocco, G. & Ruzicka, B. (2014). Colloids Surf. A Physicochem. Eng. Asp. 460, 118–122.  CrossRef CAS Google Scholar
First citationAngelini, R. & Ruzicka, B. (2015). Colloids Surf. A Physicochem. Eng. Asp. 483, 316–320.  CrossRef CAS Google Scholar
First citationBandyopadhyay, R., Liang, D., Yardimci, H., Sessoms, D. A., Borthwick, M. A., Mochrie, S. G. J., Harden, J. L. & Leheny, R. L. (2004). Phys. Rev. Lett. 93, 228302.  Web of Science CrossRef PubMed Google Scholar
First citationBarton, B. F., Graham, P. D. & McHugh, A. J. (1998). Macromolecules, 31, 1672–1679.  Web of Science CrossRef CAS Google Scholar
First citationBegam, N., Ragulskaya, A., Girelli, A., Rahmann, H., Chandran, S., Westermeier, F., Reiser, M., Sprung, M., Zhang, F., Gutt, C. & Schreiber, F. (2021). Phys. Rev. Lett. 126, 098001.  Web of Science CrossRef PubMed Google Scholar
First citationBikondoa, O. (2017). J. Appl. Cryst. 50, 357–368.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationBikondoa, O., Carbone, D., Chamard, V. & Metzger, T. H. (2012). J. Phys. Condens. Matter, 24, 445006.  Web of Science CrossRef PubMed Google Scholar
First citationBrown, G., Rikvold, P. A., Sutton, M. & Grant, M. (1997). Phys. Rev. E, 56, 6601–6612.  Web of Science CrossRef CAS Google Scholar
First citationBrown, G., Rikvold, P. A., Sutton, M. & Grant, M. (1999). Phys. Rev. E, 60, 5151–5162.  Web of Science CrossRef CAS Google Scholar
First citationCahn, J. W. & Hilliard, J. E. (1958). J. Chem. Phys. 28, 258–267.  CrossRef CAS Web of Science Google Scholar
First citationCahn, J. W. & Hilliard, J. E. (1959). J. Chem. Phys. 31, 688–699.  CrossRef CAS Web of Science Google Scholar
First citationChushkin, Y., Gulotta, A., Roosen-Runge, F., Pal, A., Stradner, A. & Schurtenberger, P. (2022). Phys. Rev. Lett. 129, 238001.  Web of Science CrossRef PubMed Google Scholar
First citationCipelletti, L., Bissig, H., Trappe, V., Ballesta, P. & Mazoyer, S. (2003). J. Phys. Condens. Matter, 15, S257–S262.  Web of Science CrossRef CAS Google Scholar
First citationCipelletti, L. & Weitz, D. A. (1999). Rev. Sci. Instrum. 70, 3214–3221.  Web of Science CrossRef CAS Google Scholar
First citationConrad, H., Lehmkühler, F., Fischer, B., Westermeier, F., Schroer, M. A., Chushkin, Y., Gutt, C., Sprung, M. & Grübel, G. (2015). Phys. Rev. E, 91, 042309.  Web of Science CrossRef Google Scholar
First citationCzakkel, O. & Madsen, A. (2011). Europhys. Lett. 95, 28001.  Web of Science CrossRef Google Scholar
First citationDainty, J. (1977). Prog. Opt. 14, 1–46.  CrossRef Google Scholar
First citationDax, I., Zaluzhnyy, I. A., Pylypenko, A., Russegger, N., Starostin, V., Rysov, R., Westermeier, F., Sprung, M., Hinderhofer, A., Pithan, L. & Schreiber, F. (2023). New J. Phys. 25, 103033.  CrossRef Google Scholar
First citationDuri, A., Bissig, H., Trappe, V. & Cipelletti, L. (2005). Phys. Rev. E, 72, 051401.  Web of Science CrossRef Google Scholar
First citationDuri, A. & Cipelletti, L. (2006). Europhys. Lett. 76, 972–978.  Web of Science CrossRef CAS Google Scholar
First citationFluerasu, A., Sutton, M. & Dufresne, E. M. (2005). Phys. Rev. Lett. 94, 055501.  Web of Science CrossRef PubMed Google Scholar
First citationGirelli, A., Rahmann, H., Begam, N., Ragulskaya, A., Reiser, M., Chandran, S., Westermeier, F., Sprung, M., Zhang, F., Gutt, C. & Schreiber, F. (2021). Phys. Rev. Lett. 126, 138004.  Web of Science CrossRef PubMed Google Scholar
First citationGoodman, J. W. (1985). Statistical Optics. New York: Wiley-Interscience.  Google Scholar
First citationGoodman, J. W. (2020). Speckle Phenomena in Optics: Theory and Applications, 2nd ed. Bellingham: SPIE Press.  Google Scholar
First citationGötze, W. (1999). J. Phys. Condens. Matter, 11, A1–A45.  Google Scholar
First citationGrübel, G., Madsen, A. & Robert, A. (2008). X-ray Photon Correlation Spectroscopy (XPCS), pp. 953–995. Dordrecht: Springer Netherlands.  Google Scholar
First citationHeadrick, R. L., Ulbrandt, J. G., Myint, P., Wan, J., Li, Y., Fluerasu, A., Zhang, Y., Wiegart, L. & Ludwig, K. F. (2019). Nat. Commun. 10, 2638.  Web of Science CrossRef PubMed Google Scholar
First citationHruszkewycz, S. O., Sutton, M., Fuoss, P. H., Adams, B., Rosenkranz, S., Ludwig, K. F., Roseker, W., Fritz, D., Cammarata, M., Zhu, D., Lee, S., Lemke, H., Gutt, C., Robert, A., Grübel, G. & Stephenson, G. B. (2012). Phys. Rev. Lett. 109, 185502.  Web of Science CrossRef PubMed Google Scholar
First citationHua, N., Zaluzhnyy, I. A., Hrkac, S. B., Shabalin, A. G. & Shpyrko, O. G. (2020). J. Synchrotron Rad. 27, 1626–1632.  CrossRef IUCr Journals Google Scholar
First citationJu, G., Xu, D., Highland, M. J., Thompson, C., Zhou, H., Eastman, J. A., Fuoss, P. H., Zapol, P., Kim, H. & Stephenson, G. B. (2019). Nat. Phys. 15, 589–594.  Web of Science CrossRef CAS Google Scholar
First citationKwaśniewski, P., Fluerasu, A. & Madsen, A. (2014). Soft Matter, 10, 8698–8704.  Web of Science PubMed Google Scholar
First citationLadd, A. J. C., Gang, H., Zhu, J. X. & Weitz, D. A. (1995). Phys. Rev. Lett. 74, 318–321.  CrossRef PubMed CAS Google Scholar
First citationLehmkühler, F., Dallari, F., Jain, A., Sikorski, M., Möller, J., Frenzel, L., Lokteva, I., Mills, G., Walther, M., Sinn, H., Schulz, F., Dartsch, M., Markmann, V., Bean, R., Kim, Y., Vagovic, P., Madsen, A., Mancuso, A. P. & Grübel, G. (2020). Proc. Natl Acad. Sci. USA, 117, 24110–24116.  Web of Science PubMed Google Scholar
First citationLehmkühler, F., Roseker, W. & Grübel, G. (2021). Appl. Sci. 11, 6179.  Google Scholar
First citationLeitner, M., Sepiol, B., Stadler, L.-M. & Pfau, B. (2012). Phys. Rev. B, 86, 064202.  Web of Science CrossRef Google Scholar
First citationLiu, W., Zheng, B., Yin, X., Yu, X., Zhang, Y., Wiegart, L., Fluerasu, A., Armstrong, B. L., Veith, G. M. & Bhatia, S. R. (2021). Appl. Mater. Interfaces, 13, 14267–14274.  CrossRef CAS Google Scholar
First citationLivet, F., Bley, F., Caudron, R., Geissler, E., Abernathy, D., Detlefs, C., Grübel, G. & Sutton, M. (2001). Phys. Rev. E, 63, 036108.  Web of Science CrossRef Google Scholar
First citationLoudon, R. (1983). The Quantum Theory of Light. Oxford University Press.  Google Scholar
First citationLu, X., Mochrie, S. G. J., Narayanan, S., Sandy, A. R. & Sprung, M. (2008). Phys. Rev. Lett. 100, 045701.  Web of Science CrossRef PubMed Google Scholar
First citationMadsen, A., Als-Nielsen, J. & Grübel, G. (2003). Phys. Rev. Lett. 90, 085701.  Web of Science CrossRef PubMed Google Scholar
First citationMadsen, A., Leheny, R. L., Guo, H., Sprung, M. & Czakkel, O. (2010). New J. Phys. 12, 055001.  Web of Science CrossRef Google Scholar
First citationMalik, A., Sandy, A. R., Lurio, L. B., Stephenson, G. B., Mochrie, S. G. J., McNulty, I. & Sutton, M. (1998). Phys. Rev. Lett. 81, 5832–5835.  Web of Science CrossRef CAS Google Scholar
First citationMartinez, V. A., Bryant, G. & van Megen, W. (2010). J. Chem. Phys. 133, 114906.  CrossRef PubMed Google Scholar
First citationMöller, J., Reiser, M., Hallmann, J., Boesenberg, U., Zozulya, A., Rahmann, H., Becker, A.-L., Westermeier, F., Zinn, T., Sprung, M., Narayanan, T., Gutt, C. & Madsen, A. (2021). New J. Phys. 23, 093041.  Google Scholar
First citationMüller, L., Waldorf, M., Gutt, C., Grübel, G., Madsen, A., Finlayson, T. R. & Klemradt, U. (2011). Phys. Rev. Lett. 107, 105701.  Web of Science PubMed Google Scholar
First citationNarayanan, S., Lee, D. R., Hagman, A., Li, X. & Wang, J. (2007). Phys. Rev. Lett. 98, 185506.  Web of Science CrossRef PubMed Google Scholar
First citationNogales, A. & Fluerasu, A. (2016). Eur. Polym. J. 81, 494–504.  Web of Science CrossRef CAS Google Scholar
First citationOrsi, D., Cristofolini, L., Fontana, M. P., Pontecorvo, E., Caronna, C., Fluerasu, A., Zontone, F. & Madsen, A. (2010). Phys. Rev. E, 82, 031804.  Web of Science CrossRef Google Scholar
First citationPerakis, F. & Gutt, C. (2020). Phys. Chem. Chem. Phys. 22, 19443–19453.  Web of Science CrossRef CAS PubMed Google Scholar
First citationPusey, P. & Van Megen, W. (1989). Physica A, 157, 705–741.  CrossRef CAS Google Scholar
First citationPusey, P. N. & Jakeman, E. (1975). J. Phys. A Math. Gen. 8, 392–410.  CrossRef CAS Google Scholar
First citationRagulskaya, A., Begam, N., Girelli, A., Rahmann, H., Reiser, M., Westermeier, F., Sprung, M., Zhang, F., Gutt, C. & Schreiber, F. (2021). J. Phys. Chem. Lett. 12, 7085–7090.  Web of Science CrossRef CAS PubMed Google Scholar
First citationRagulskaya, A., Starostin, V., Begam, N., Girelli, A., Rahmann, H., Reiser, M., Westermeier, F., Sprung, M., Zhang, F., Gutt, C. & Schreiber, F. (2022). IUCrJ, 9, 439–448.  Web of Science CrossRef CAS PubMed IUCr Journals Google Scholar
First citationReiser, M., Girelli, A., Ragulskaya, A., Das, S., Berkowicz, S., Bin, M., Ladd-Parada, M., Filianina, M., Poggemann, H.-F., Begam, N., Akhundzadeh, M. S., Timmermann, S., Randolph, L., Chushkin, Y., Seydel, T., Boesenberg, U., Hallmann, J., Möller, J., Rodriguez-Fernandez, A., Rosca, R., Schaffer, R., Scholz, M., Shayduk, R., Zozulya, A., Madsen, A., Schreiber, F., Zhang, F., Perakis, F. & Gutt, C. (2022). Nat. Commun. 13, 5528.  CrossRef PubMed Google Scholar
First citationRoseker, W., Hruszkewycz, S. O., Lehmkühler, F., Walther, M., Schulte-Schrepping, H., Lee, S., Osaka, T., Strüder, L., Hartmann, R., Sikorski, M., Song, S., Robert, A., Fuoss, P. H., Sutton, M., Stephenson, G. B. & Grübel, G. (2018). Nat Commun, 9, 1704.  CrossRef PubMed Google Scholar
First citationRuta, B., Baldi, G., Monaco, G. & Chushkin, Y. (2013). J. Chem. Phys. 138, 054508.  Web of Science CrossRef PubMed Google Scholar
First citationRuta, B., Chushkin, Y., Monaco, G., Cipelletti, L., Pineda, E., Bruna, P., Giordano, V. M. & Gonzalez-Silveira, M. (2012). Phys. Rev. Lett. 109, 165701.  Web of Science CrossRef PubMed Google Scholar
First citationRuta, B., Zontone, F., Chushkin, Y., Baldi, G., Pintori, G., Monaco, G., Rufflé, B. & Kob, W. (2017). Sci. Rep. 7, 3962.  Web of Science CrossRef PubMed Google Scholar
First citationSeydel, T., Madsen, A., Tolan, M., Grübel, G. & Press, W. (2001). Phys. Rev. B, 63, 073409.  Web of Science CrossRef Google Scholar
First citationShpyrko, O. G. (2014). J. Synchrotron Rad. 21, 1057–1064.  Web of Science CrossRef CAS IUCr Journals Google Scholar
First citationSinha, S. K., Jiang, Z. & Lurio, L. B. (2014). Adv. Mater. 26, 7764–7785.  Web of Science CrossRef CAS PubMed Google Scholar
First citationSutton, M. (2008). C. R. Phys. 9, 657–667.  Web of Science CrossRef CAS Google Scholar
First citationSutton, M., Laaziri, K., Livet, F. & Bley, F. (2003). Opt. Express, 11, 2268.  Web of Science CrossRef PubMed Google Scholar
First citationTimmermann, S., Anthuparambil, N. D., Girelli, A., Begam, N., Kowalski, M., Retzbach, S., Senft, M. D., Akhundzadeh, M. S., Pog­gemann, H., Moron, M., Hiremath, A., Gutmüller, D., Dargasz, M., Öztürk, Ö., Paulus, M., Westermeier, F., Sprung, M., Ragulskaya, A., Zhang, F., Schreiber, F. & Gutt, C. (2023). Sci. Rep. 13, 11048.  CrossRef PubMed Google Scholar
First citationWestermeier, F., Fischer, B., Roseker, W., Grübel, G., Nägele, G. & Heinen, M. (2012). J. Chem. Phys. 137, 114504.  Web of Science CrossRef PubMed Google Scholar
First citationWilliams, G. & Watts, D. C. (1970). Trans. Faraday Soc. 66, 80–85.  CrossRef CAS Web of Science Google Scholar
First citationZand, D. D. van 't, Chushkin, Y., Belkoura, L., Lobo, C. V., Strey, R., Lyakhova, K. & Clegg, P. S. (2012). Soft Matter, 8, 4062–4066.  Google Scholar
First citationZhang, Q., Dufresne, E. M., Chen, P., Park, J., Cosgriff, M. P., Yusuf, M., Dong, Y., Fong, D. D., Zhou, H., Cai, Z., Harder, R. J., Callori, S. J., Dawber, M., Evans, P. G. & Sandy, A. R. (2017). Phys. Rev. Lett. 118, 097601.  CrossRef PubMed Google Scholar
First citationZhang, Q., Dufresne, E. M., Nakaye, Y., Jemian, P. R., Sakumura, T., Sakuma, Y., Ferrara, J. D., Maj, P., Hassan, A., Bahadur, D., Ramakrishnan, S., Khan, F., Veseli, S., Sandy, A. R., Schwarz, N. & Narayanan, S. (2021). J. Synchrotron Rad. 28, 259–265.  Web of Science CrossRef IUCr Journals Google Scholar

This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.

Journal logoJOURNAL OF
APPLIED
CRYSTALLOGRAPHY
ISSN: 1600-5767
Follow J. Appl. Cryst.
Sign up for e-alerts
Follow J. Appl. Cryst. on Twitter
Follow us on facebook
Sign up for RSS feeds