research papers
Measurable structure factors of dense dispersions containing polydisperse optically inhomogeneous particles
aInstitut für Chemie, Universität Rostock, 18051 Rostock, Germany
*Correspondence e-mail: joachim.wagner@uni-rostock.de
Here, it is investigated how optical properties of single scatterers in interacting multi-particle systems influence measurable structure factors. Both particles with linear gradients of their scattering length density and core–shell structures evoke characteristic deviations between the weighted sum 〈S(Q)〉 of factors in a multi-component system and experimentally accessible measurable structure factors SM(Q). While 〈S(Q)〉 contains only the structural information of self-organizing systems, SM(Q) is additionally influenced by the optical properties of their constituents, resulting in features such as changing amplitudes, additional peaks in the low-wavevector region or splitting of higher-order maxima, which are not related to structural reasons. It is shown that these effects can be systematically categorized according to the qualitative behaviour of the form factor in the Guinier region, which enables assessing the suitability of experimentally obtained structure factors to genuinely represent the microstructure of complex systems free from any particular model assumption. Hence, a careful data analysis regarding size distribution and optical properties of single scatterers is mandatory to avoid a misinterpretation of measurable structure factors.
Keywords: small-angle scattering; contrast variation; structure factors; hard-sphere interaction; dense dispersions; polydisperse particles; optical inhomogeneities; size distributions.
1. Introduction
Colloidal dispersions attract wide interest in condensed matter physics as highly tunable model systems, mimicking atoms and molecules on the much larger mesoscopic scale with typical length scales between 10 and 1000 nm. Studying these systems enabled major advances in the comprehension of the characteristics of simple fluids and solids, which stimulated the progress of significant theoretical developments towards the understanding of complex systems and materials (Lu & Weitz, 2013).
Scattering experiments serve as essential methods for structural and dynamical investigations in colloidal many-particle systems (Li et al., 2016). Small-angle scattering (SANS with neutrons or SAXS with X-rays as a probe) enables the characterization of colloidal suspensions across the entire range of relevant scattering vector magnitudes Q (Glatter, 2018). By employing visible light, which is also a natural choice since its wavelength is of the same order of magnitude as the typical size of a colloidal particle, the same type of analysis is in principle also possible in a simpler laboratory setup. This is however connected with the cost of a limited resolution and, as a consequence thereof, the restriction to comparatively large structures (Bohren & Huffmann, 2008).
In non-interacting systems, the positions and orientations of the colloidal particles are completely uncorrelated. Thus, the scattered intensity results solely from the superposition of the scattering functions of the single constituents. However, when the particles do interact, higher-level structures emerge from inherent self-organization due to interparticle forces, such as electrostatic and steric interactions or van der Waals attractions. The intensity is then influenced both by the optical properties of the scatterers themselves and by the spatial correlations between them. For idealized radially symmetric and monodisperse systems, where all particles are assumed to be identical, the two contributions can be rigorously separated into the form factor P(Q), containing the single-particle properties, and the S(Q), which encodes the structural correlations, employing the well known factorization I(Q) ∝ P(Q) S(Q) (Hansen et al., 1991).
Realistic dispersions typically exhibit a distribution of characteristics, prominently through particle size. In polydisperse interacting systems, the characterization via scattering experiments is in general significantly more complicated, as the factorization of the intensity into form factor and ). Additionally, the observed diffraction patterns become increasingly featureless for broader size distributions, further obstructing the interpretation of experimental intensities. The analysis of multi-component systems thus requires a thorough understanding of the underlying distributions of scattering properties and particle interactions. Insights can be gained through contrast-variation techniques (Ballauff, 2001); selectively altering the contrast between specific particle types or between particles and the surrounding medium allows for the isolation and probing of distinct species, aiding the validation of theoretical models.
can no longer be employed in a straightforward way (Salgi & Rajagopalan, 1993From an experimental standpoint, it is useful to analyse the measurable SM(Q), which is defined in such a way that the factorization property I(Q) ∝ P(Q) SM(Q) is also recovered in the polydisperse case (Hansen et al., 1991). SM(Q) is comparatively easy to access experimentally from the ratio between the intensities of an interacting suspension and a highly diluted non-interacting one. It is as such also widely used as a measure for structural correlations in polydisperse systems, where the height of the principal peak is especially well established as an (Banchio et al., 1998). Under specific circumstances, this type of analysis can however turn into a serious pitfall: SM(Q) is also fundamentally affected by optical properties of the particles and not only by their interactions (Salgi & Rajagopalan, 1993).
For certain types of dispersions, some simplifying assumptions can be employed. In dilute suspensions of strongly interacting charged particles, for example, the interparticle distances typically are about an order of magnitude larger than the particle sizes because of the large electrostatic repulsion (Hayter & Penfold, 1981). In such a case, the correlation between the particle positions and the scattering amplitudes can be neglected. This neglect of correlations leads to the `decoupling approximation', under which SM(Q) can be decomposed into a that genuinely represents the averaged structural correlations and weighting factors solely dependent on the scattering amplitudes (Pusey et al., 1982; Kotlarchyk & Chen, 1983). This type of analysis was also recently used in a SANS study of moderately concentrated poly(N-isopropylacrylamide) microgels (Zhou et al., 2023), where again the importance of an accurate treatment of polydispersity was stressed. Non-spherical particles cause effects in scattering patterns that appear quite similar to those introduced by size dispersity in systems of spherical particles (Pusey et al., 1982). These effects can to a certain degree also be treated within the decoupling approximation (Kotlarchyk & Chen, 1983), whose range of validity has been extensively characterized (Greene et al., 2016). In highly concentrated suspensions, where particles are in close contact, the accuracy of the decoupling approximation is strongly diminished, as in these systems the correlation lengths of the particles' centres of mass are comparable to the correlation lengths inside the particles themselves (Pedersen, 1997). At such high particle volume fractions, excluded volume effects are the predominant contribution to the total interaction potential.
The fundamental interactions in dense colloidal dispersions consisting of spherical particles can to a good approximation be theoretically described with the hard-sphere model (Kirkwood & Boggs, 1942; Widom, 1967). The radial distribution functions of an n-component mixture of hard spheres can be calculated within the Percus–Yevick closure of the Ornstein–Zernike equation (Percus & Yevick, 1958) using Baxter's technique (Baxter, 1970), giving access to the corresponding factors (Vrij, 1978, 1979; Blum & Stell, 1979, 1980). Building on Vrij's work (Vrij, 1979; van Beurten & Vrij, 1981), Griffith et al. (1987) presented an analytical scattering function of a polydisperse hard-sphere fluid with a Schulz–Flory distribution (Flory, 1936; Schulz, 1939) of particle diameters. Despite their helpfulness, these expressions are not widely used because of their perceived complexity. Nayeri et al. (2009) later extended this approach to core–shell structured hard spheres and used their expressions to describe experimental intensities of a hard-sphere-like microemulsion system. Only recently, Botet et al. (2020) provided expressions for SM(Q) in a simple accessible form and for a number of commonly encountered size distributions. Their analytical expressions are valid for hard optically homogeneous spheres.
This resurgence of interest is an incentive to systematically examine how different form-factor models affect the characteristics of measurable structure factors. It is a well known fact that, especially in dense dispersions, SM(Q) is generally not equal to the 〈S(Q)〉 representing the averaged spatial correlations of the entire system (Salgi & Rajagopalan, 1993; Frenkel et al., 1986). The mismatch between SM(Q) and 〈S(Q)〉 is precisely the reason why in many past studies an extraction of SM(Q) is deliberately not attempted. By describing the scattered intensity of concentrated suspensions with sophisticated models, a thorough characterization of particle properties is possible without calculating SM(Q), as performed, for example, by Stieger, Pedersen et al. (2004); Stieger, Richtering et al. (2004); Zackrisson et al. (2005); Balogh et al. (2007); or Scotti (2021).
Certainly, however, for many applications in condensed matter science, gaining an accurate approximation of the et al., 2020; Peláez-Fernández et al., 2011; Anta & Madden, 1999; Krause et al., 1991; Stellbrink et al., 2002). Beyond using structure factors to deduct structural patterns from experimental scattering data (Mohanty et al., 2017; Phalakornkul et al., 1996; Scheffold & Mason, 2009; Mason et al., 2006), these quantities are used in theoretical approaches to calculate short- and long-time dynamics in many-particle systems. Two notable examples in this context are the δγ-expansion by Beenakker & Mazur (1983, 1984) used to model hydrodynamic effects during short-time diffusion (Genz & Klein, 1991) and the widely known mode-coupling theory of the (Janssen, 2018), which both need the static as an input for calculations. In these theories, a popular method is to circumvent computationally expensive multi-component calculations by considering an effective one-component analysis using the experimentally obtained SM(Q) directly, assuming this quantity to be an accurate representation of the average 〈S(Q)〉 (Robert et al., 2008; Di Cola et al., 2009).
is still highly desired. This is especially the case when employing computer simulations or many-body theory to model structural phenomena and wanting to compare detailed facets of particle self-organization directly with experimental outcomes (DekkerFor a number of specialized cases, deviations between 〈S(Q)〉 and SM(Q) and the general influence of varying optical properties on SM(Q) have already been assessed (Banchio et al., 1998; Pedersen, 2001). The purpose of this contribution is to raise further awareness on how particles' optical properties influence the shape of SM(Q) while the underlying interactions remain unchanged and to show that these observations can be model-independently systematized on the basis of quite universal principles. This enables practitioners to make informed judgements under which circumstances such an experimentally obtained can still serve as a valid We show typical examples of shapes that can be realistically encountered during contrast-variation experiments, so even without explicitly employing established theoretical models [see e.g. Pedersen (1997) for a large collection of scattering functions], a qualitative assessment of experimental findings is possible. We also analyse two simplified models for optically inhomogeneous particles: those with a linear gradient of the scattering contrast and spheres with a core–shell structure. Nevertheless, the approach is readily adaptable to any model and provides a toolbox for the modelling of measurable structure factors for hard-sphere suspensions with arbitrary form factors, as demonstrated in several past studies (Vrij, 1979; Frenkel et al., 1986; Pedersen, 2001; Nayeri et al., 2009; Botet et al., 2020; Diaz Maier & Wagner, 2024).
2. Scattering of hard-sphere mixtures
We consider a mixture of spherical particles, where each particle can be categorized into one of n species. The composition of the mixture is specified by the number fractions xα = Nα/N, where N is the total number of particles and Nα is the number of particles belonging to species α. We further restrict ourselves to elastic single scattering events where the Born approximation is applicable. In such a case, the mean intensity
is proportional to the weighted sum of the single-particle scattering amplitudes fα(Q) and the factors Sαβ(Q) (Salgi & Rajagopalan, 1993). Herein, the scattering amplitude
is the Fourier–Bessel transform of the scattering contrast ρα(r), whereas the factors Sαβ(Q) are obtained from the solution of the multi-component Ornstein–Zernike equation. Expressions for Sαβ(Q) of the hard-sphere fluid within the Percus–Yevick closure are given by Vrij (1979) but, for the convenience of the reader, the solution is re-articulated in Appendix A, and presented in a manner that is accessible and easily applicable. This genuine multi-component approach involving factors is a general and versatile formalism, from which both the decoupling approximation (Kotlarchyk & Chen, 1983) and the local monodisperse approximation (Pedersen, 1994), serving as a localized effective one-component approach suited for very polydisperse systems, can be derived.
For non-interacting particles, the Sαβ(Q) = δαβ, where δαβ denotes the Kronecker symbol. Equation (1) then reduces to the size average of the squared scattering amplitudes,
factors are simplyThe average form factor
is familiarly obtained from the normalization to I(Q) ∝ P(Q) SM(Q), the expression
As the measurable should satisfy the relationresults from the combination of equations (1) and (3). The averaged
provides information about the total spatial correlations between all present particles, regardless of their species labels. It represents a true thermodynamic average independent of any optical properties. Any deviation between SM(Q) and 〈S(Q)〉 is thus a measure for the perturbation of 〈S(Q)〉 caused by the scattering amplitudes.
We now want to explore the influence of the scattering amplitudes on the shape of SM(Q). The aim is to gain a qualitative understanding of generic patterns; so to keep the analysis tractable, only a single representative size distribution is considered. For this purpose, the Schulz–Flory distribution with probability density
is chosen. Here, R is the particle radius with mean 〈R〉 and Γ(x) represents the gamma function. The polydispersity p of the system is specified by the shape parameter Z via p2 = (〈R2〉 − 〈R〉2)/〈R〉2 = 1/(Z + 1). The idea is now to discretize the distribution to a representative n-component mixture. For the Schulz–Flory distribution, an efficient way to achieve this is by exploiting the generalized Gauss–Laguerre quadrature rule, specifically used to calculate integrals with a weighting function like equation (7) (D'Aguanno & Klein, 1992; D'Aguanno, 1993; Olver et al., 2010). The nodes and weights generated by such a procedure are equivalent to the particle radii and number fractions of a discrete mixture which shares the first 2n − 1 moments 〈Rn〉 with the original continuous distribution. For each calculated scattering function, we carefully checked that the number of nodes necessary for convergence was reached. The numerical scheme was further tested against the analytical SM(Q) for homogeneous spheres provided by Botet et al. (2020), where excellent agreement was found.
3. Measurable structure factors of polydisperse systems
3.1. General remarks
Fig. 1 provides a general introductory overview of the influence of polydispersity on P(Q), SM(Q) and 〈S(Q)〉, discussed for a dense suspension of homogeneous spheres, serving as a reiteration of well known phenomenology (Botet et al., 2020). Concerning the form factors, only those corresponding to polydispersities of less than 10% appear structured. Familiarly, the characteristic minima in P(Q) become increasingly smeared out for broader size distributions.
Polydispersity also causes a change in the initial slope of P(Q) in the Guinier region. Reflecting the distribution of particle sizes when calculating the Taylor expansion of P(Q), the slope is now given by , where the familiar RG is substituted by an apparent (Glatter, 2018; Tomchuk et al., 2014). For homogeneous spheres,
is obtained, which reduces to the well known result of for monodisperse systems.
Similarly to P(Q), both the measurable SM(Q) and the average 〈S(Q)〉 become increasingly featureless at high polydispersities, which is distinctively noticeable as the principal peak's amplitude decreases and the secondary oscillations gradually disappear. Shifting the focus to direct comparison between the two structure factors SM(Q) and 〈S(Q)〉, multiple observations are apparent. While the amplitude of the principal peak is similar for both functions, differences appear at larger wavevectors, where secondary peaks in SM(Q) appear at roughly the locations of the form-factor minima, as similarly noticed by Ginoza & Yasutomi (1999). With increasing polydispersity, these maxima evolve into broad shoulders that get smeared out eventually. As also noted by Ginoza & Yasutomi (1999), sharp secondary maxima are hard to observe experimentally because a very narrow size distribution in combination with a homogenous distribution of the scattering length density (SLD) inside the particles is required. On the other hand, shoulder-like features in experimentally determined structure factors are well documented [see, as an example, Di Cola et al. (2009)]. In the low-Q region, a striking observation is the significant increase of 〈S(0)〉 at elevated polydispersities in comparison with SM(0). According to the fluctuation-dissipation theorem from statistical mechanics, the isothermal compressibility κT is for monodisperse systems connected to the zero-wavevector limit of S(Q) via S(0) = ρkBTκT, where ρ denotes the and kBT is the thermal energy. The extension of this concept to mixtures must however be treated with caution because, for multi-component systems, the connection between structure and thermodynamics is not simply given by the size average 〈S(0)〉. According to the Kirkwood–Buff theory of solutions, it is instead given by the relation , where is the αβ element of the inverse structure-factor matrix (Hansen & McDonald, 2013).
3.2. Linear contrast gradient
As a prototypical example for particles with inhomogeneous scattering strength, particles with a linear gradient of the SLD are investigated. This is particularly relevant for swellable particles into which the suspension medium can diffuse. This can occur with microgel particles (Karg et al., 2019), for example. Under certain reaction conditions, an inhomogeneous degree of cross-linking arises, which also leads to inhomogeneous scattering properties. Particles with intrinsic material gradients are also plausible, obtained for example by continuously changing the monomer composition in a feed process during synthesis. Then, in principle, a suspension in which the contrast within a particle changes its sign can also be realized. The form of a linear gradient is assumed for the sake of simplicity in order to investigate the phenomenology of continuous contrasts as an example.
The scattering contrast as a function of the distance r from the centre can for a single particle be parametrized as
where R is the particle radius, ρ0 is the contrast in the centre and ρR is the contrast at the interface to the surrounding medium. Accordingly, the resulting single-particle scattering amplitude is given by
which reduces to
in the forward-scattering limit. A closer look at equation (11) reveals that the forward-scattering contribution disappears if the condition ρR/ρ0 = −1/3 is fulfilled. In particular, when the maximum accessible scattering vector is limited, as in the case of contributes significantly to the total scattering If the is zero, the sample appears almost optically transparent. Refractive-index matching can be achieved for particles with a homogeneous scattering capacity if the SLD of the suspension medium is adapted to that of the particles. If the scattering capacity is inhomogeneous, index matching can only minimize the total scattering which is often achieved by making the almost zero. In the following, the condition when the power is minimal is referred to as the index match point.
To gain a systematic understanding of the behaviour of the measurable SM(Q) as a function of the contrast ratio ρR/ρ0, it will prove advantageous to investigate the Guinier region of the form factor. Using the contrast profile from equation (9), for a single particle with radius R,
is obtained for the effective squared ρ0 and ρR. For polydisperse suspensions, a similar expression emerges:
which depends not only on the particle's radius but also on the two contrast parametersAs such, the contrast dependence of the prefactor is not altered by polydispersity and the qualitative discussion can instead be based on monodisperse suspensions. We will thus refer to the prefactor simply as , even in the polydisperse case.
Inspecting equation (12), several characteristic ratios ρR/ρ0 are apparent: becomes zero for ρR/ρ0 = −1/5; exhibits a pole at ρR/ρ0 = −1/3, incident with the index match point; and has an asymptotic limit of for ρR/ρ0 → ±∞. Here, it will be shown that the behaviour of the scattering functions can be divided into three qualitatively distinct classes, and that form factors and measurable structure factors within each domain share unique features. The classification based on the behaviour of , together with form factors P(Q) and measurable structure factors SM(Q) representative of each region, is visualized in Fig. 2. The regions are characterized as follows:
(I) For ρR/ρ0 > −1/5, is positive and the form factors have the familiar decaying shape known from homogeneous spheres. With decreasing contrast ratio, the decay becomes increasingly gradual until is reached for ρR/ρ0 = −1/5. Around the principal peak of SM(Q) and for lower wavevectors, changes in the contrast have a negligible influence on the measurable structure factors. However, at wavevectors beyond the principal peak's location, SM(Q) is greatly affected by contrast variation. Depending on the specific location of the first form-factor minimum, which shifts to larger wavevectors with lower contrast ratios, the shoulder-like artefact also visible in Fig. 1 moves through SM(Q) towards larger wavevectors and therein most prominently affects the shape of the first local minimum and the following secondary maximum.
(II) For the contrast ratios −1/3 < ρR/ρ0 < −1/5, becomes negative, which implies an imaginary RG leading to a positive initial slope of P(Q). Form factors in this region therefore initially increase from P(0) = 1 until a global maximum is reached at QR ≃ 4, after which they decay. The height of the maximum increases as the contrast ratio moves towards the index match point at ρR/ρ0 = −1/3. Curiously, the measurable structure factors in this domain are almost indistinguishable, even though the variation of is much more pronounced in comparison with region (I), where the span of is small but SM(Q) shows a much more diverse behaviour. Also, the distorting artefacts from region (I) disappear almost completely.
(III) Contrast ratios of ρR/ρ0 < −1/3 again result in positive and negative initial slopes. Close to the index match point, where is comparatively large, P(Q) exhibits an intriguing shape. At small wavevectors, a pronounced minimum occurs even in very polydisperse suspensions. Beyond the minimum, P(Q) rises to a global maximum reminiscent of region (II). For contrasts in this range, an additional local maximum in SM(Q) appears at low wavevectors, caused by the presence of the first form-factor minimum. Such secondary maxima are often discussed in the literature as an indication of self-organization on length scales beyond the distance of nearest neighbours, i.e. the formation of correlated clusters (Sciortino et al., 2004; Liu et al., 2005). The secondary maxima occurring here are exclusively caused by the scattering amplitudes and cannot be attributed to structural properties of the sample. This constitutes a valuable example of a situation where a careless inspection of experimentally determined SM(Q) can in the worst case lead to unjustified assumptions about the structure of a system. Moving further away from the index match point, the first form-factor minimum moves towards larger wavevectors and gets shallower. At the same time, the following maximum declines and, as such, the shape of P(Q) morphs back into the familiar decaying shape from region (I). Simultaneously, the location of the secondary maximum in SM(Q) drifts towards higher wavevectors. Fig. 2 also displays a situation where the form-factor minimum exactly coincides with the location where the principal peak of SM(Q) would normally occur. In this case, the main peak is drastically diminished, which is again not an indicator for a less pronounced short-range order in this particular instance, but can certainly be mistaken as such.
The principal-peak height of a ) noticed in their investigation of highly concentrated nanoemulsions that the peak amplitude in SM(Q) is deeply affected by polydispersity. As such, the evolution of this height during contrast variation is also of special interest. Fig. 3 compares the peak height of the average 〈S(Qmax)〉 with the value of SM(Qmax) at the same wavevector as a function of the contrast ratio ρR/ρ0 and for different degrees of polydispersity. Overall, it is clearly shown that SM(Qmax) is deeply affected by changes in the contrast. There exist two contrast ratios where SM(Qmax) and 〈S(Qmax)〉 coincide. One of them is to a good approximation given by ρR/ρ0 ≃ −1/5, the location where the apparent disappears and P(Q) decays very slowly. The other location is at a positive contrast ratio and drifts towards higher ρR/ρ0 with increasing polydispersity. Bounded by those two ratios is a regime where SM(Qmax) exceeds 〈S(Qmax)〉, while for all other contrast ratios, the peak height from SM(Qmax) underestimates the actual height. For comparatively small polydispersities around 5%, the deviation from 〈S(Qmax)〉 is small and only amounts to a few per cent, as long as the contrast ratio is larger than ρR/ρ0 ≃ −1/5. For lower ratios, SM(Qmax) is strongly diminished, most pronouncedly at contrast ratios of ρR/ρ0 ≃ −1. For higher polydispersities, the deviations become even more severe, as best visualized in Fig. 3(b), where the relative deviation between 〈S(Qmax)〉 and SM(Qmax) is depicted. Even in the immediate vicinity of 〈S(Qmax)〉 = SM(Qmax), already deviations of the order of 5–10% appear for the highest shown polydispersities. This demonstrates that, no matter what the actual degree of polydispersity, SM(Qmax) can only serve as a reliable for very specific contrast ratios.
is an often employed structural Scheffold & Mason (20093.3. Core–shell particles
Core–shell models are commonly employed to describe particles consisting of different layers of material, e.g. nanoparticles with grafted shells (Hallett et al., 2020; Diaz Maier & Wagner, 2024) or micellar structures (Szymusiak et al., 2017). As core and shell naturally differ in their material properties, in principle both positive and negative contrast differences with respect to the surrounding medium can occur, similarly to particles with continuous material gradients. For Schulz–Flory-distributed core–shell particles, analytical expressions for the form factor P(Q) exist in the case of a polydisperse core and a shell of constant thickness (Bartlett & Ottewill, 1992), for a polydisperse total diameter and a constant core-to-shell ratio (Wagner, 2004), and for both core radius and shell thickness independently distributed (Wagner, 2012). Moreover, an analytical solution for the problem of correlated hard-sphere core–shell systems was provided by Nayeri et al. (2009).
The scattering amplitude of a single core–shell particle,
is the sum of the amplitudes of a sphere and a spherical shell, weighted by their respective contrasts, ρc and ρs. Rc and R are the core radius and the total radius of the particle, respectively, and we specifically consider the case where the core radius and the total radius are connected by a constant species-independent size ratio δ = Rc/R.
Similar to the gradient model, the forward-scattering contribution
disappears for specific contrast combinations of the ratio of contrasts ρs/ρc = δ3/(δ3 − 1), which now additionally depends on the size ratio δ. For the effective of a polydisperse system, an expression with similar structure to equation (13) emerges:
That, again, a prefactor containing the contrasts can be decoupled from the size average is a peculiarity of this model with constant size ratio and a key reason why this assumption was made for this investigation.
In Fig. 4, the contrast dependence of is visualized for different size ratios δ. As in the case of spheres with a linear gradient of the SLD, this results in hyperbola-like curves, where the location of the pole is now influenced by δ; an increasing ratio of core diameter to total diameter shifts the location of the pole to more negative contrast ratios (ρs/ρc). The contrast ratio where is, in comparison, only slightly altered by δ. This leads to a larger range of contrast ratios with negative as the shell thickness decreases.
This shows that core–shell particles exhibit qualitatively comparable optical characteristics to particles with a linear density gradient. As such, the form factors P(Q) of core–shell systems can likewise be categorized into three classes based on their behaviour at low wavevectors. Example form factors for each class are also visualized in Fig. 4.
Because of these similarities, we focus the remainder of the discussion on aspects that are unique to core–shell particles, i.e. how measurable structure factors are influenced by different core-to-shell ratios. For this purpose, structure factors corresponding to two important edge cases, particles with a small core and particles with a thin shell, are compared in Fig. 5 for different degrees of polydispersity and chosen contrast ratios ρs/ρc. Core–shell models with thin shells are often encountered when characterizing particles stabilized by a grafted polymer layer, which are prototypical colloidal model particles displaying hard-sphere behaviour (Royall et al., 2013). The case of hard spheres with a strongly scattering small core and a weakly scattering comparatively large shell is equally of interest. Under these conditions, essentially, the behaviour of highly charged strongly repelling particles whose interparticle distance is several times larger than their diameter is artificially mimicked. For these systems, the measurable SM(Q) should in theory to a good approximation coincide with the average 〈S(Q)〉. To reasonably compare models with different size ratios (δ), two specific contrast ratios (ρs/ρc) are depicted: the ratio ρs/ρc = δ3/(δ3 − 1) at the index matching point, where forward scattering is minimized; and the ratio ρs/ρc = δ5/(δ5 − 1), where and P(Q) shows the weakest decay. In the case of , both conditions basically lead to the same result: the shell is virtually hidden with ρs ≃ 0.
As can be observed in Fig. 5, for moderate polydispersities of 5–10%, the small core-to-total ratio δ = 0.1 indeed yields measurable structure factors SM(Q) that are indistinguishable from 〈S(Q)〉 for both depicted contrast ratios. For particles with thin shells (δ = 0.9), SM(Q) and 〈S(Q)〉 also agree well in the vicinity of the principal peak. However, differences arise around the secondary maxima, where the peak amplitudes in SM(Q) are diminished because of the interference of the scattering amplitudes. With increasing polydispersity, this deviation becomes more pronounced. Still, even for particles that are seemingly quite close to homogeneous spheres, artefacts in SM(Q) can be significantly reduced by careful contrast variation.
Looking at highly polydisperse systems, it is evident that, even for rather small cores with δ = 0.1, 〈S(Q)〉 cannot be accurately represented by any SM(Q). Only the height of the principle peak is correctly estimated. This stresses again the importance of an accurate treatment of very broad size distributions, where any kind of approximation must be carefully checked for validity.
4. Conclusions
Colloidal dispersions generally exhibit a particle size distribution, which needs to be taken into account when interpreting results from scattering experiments. The measurable SM(Q) is an experimental, comparatively easily accessible measure for the interparticle structure in interacting systems. However, in polydisperse systems, SM(Q) is, beyond the structural correlations, also decisively affected by the optical properties of the individual particles. To this end, we systematically investigated the influence of different form-factor models on the shape of SM(Q) of dense dispersions with hard-sphere interactions. The characterization of measurable structure factors was extended to two classes of spherical particles with inhomogeneous scattering capacity: first, spheres with a linear SLD profile as a general model for particles with continuous contrast gradients and, second, a core–shell system as a prototype for particles with layered structures.
For both models, we find that the structure factors can be categorized into three distinctive classes of shared qualitative features, based on the behaviour of the form factor P(Q) in the Guinier region. SM(Q) can, for these optically inhomogeneous model particles, be significantly influenced by the variation of the scattering contrasts relative to the surrounding medium. Depending on the specific contrast combination, shoulder-like features emerge, maxima are diminished or split, and even secondary maxima in the low-wavevector region, reminiscent of cluster peaks, can be observed. These effects are solely due to the optical properties of the particles and are not caused by structural changes in the sample. We further showed that the height of the principal peak of SM(Q) can only be regarded as a representative in a very restricted range of contrasts, especially for broad size distributions.
These observations emphasize the need to properly address the distribution of particle size (and possibly also other characteristics) in the interpretation of static scattering experiments. Actually, for many applications, deliberately broad size distributions are a desired feature; an academically relevant example is studies of deeply supercooled glass-forming systems (Ninarello et al., 2017), where crystallization needs to be suppressed and where polydispersity effects in any form certainly cannot be neglected (Zaccarelli et al., 2015; Pihlajamaa et al., 2023).
Beyond providing an enhanced qualitative understanding of features that can possibly be encountered when analysing experimentally extracted measurable structure factors, the numerical scheme presented in this contribution in principle provides a means to model the scattered intensity of any polydisperse hard-sphere system, provided a model for the single-particle scattering amplitude and an appropriate size distribution is available. Performing fits with such advanced models directly on experimentally observed intensities gives access to the underlying et al., 2022).
factors, enabling a characterization and possible further theoretical analysis on a genuine multi-component foundation, rather than employing effective one-component approaches. The current restriction to hard-sphere interactions is a major incentive to promote advancements in the analytical evaluation of factors for other interaction potentials, since numerically solving integral equations or employing computer simulations with reasonable statistics are currently only realistically feasible for a restricted number of components, especially in mixtures with large size disparities (AllahyarovAPPENDIX A
Percus–Yevick structure factors for hard-sphere mixtures
The analytical solution of the Ornstein–Zernike equation for the hard-sphere potential within the Percus–Yevick closure in terms of the Sαβ(Q), presented by Vrij (1979) and reformulated by Voigtmann (2003), is restated here. In short, an expression for the partial direct correlation functions cαβ(r) in real space can be found using Baxter's factorization technique (Baxter, 1970). The transformed solution in wavevector space cαβ(Q) can subsequently be used to obtain the factors Sαβ(Q).
factorsLet φ be the total of all spheres, dα be the diameter and xα be the of the spheres of species α. The total ρ of the system is related to the by . With the abbreviations
and
the set of coefficients
and
can be determined. Furthermore, by introducing and , the terms
and
can be calculated, which finally leads to
The partial direct correlation functions form the matrix C with elements Cαβ = (xαxβ)1/2cαβ, which is related to the matrix of factors S by the Ornstein–Zernike relation
The δαβ is the Kronecker delta.
factors here are defined within the convention , whereAcknowledgements
Open access funding enabled and organized by Projekt DEAL.
References
Allahyarov, E., Löwen, H. & Denton, A. R. (2022). Phys. Chem. Chem. Phys. 24, 15439–15451. CrossRef CAS PubMed Google Scholar
Anta, J. A. & Madden, P. A. (1999). J. Phys. Condens. Matter, 11, 6099–6111. CrossRef CAS Google Scholar
Ballauff, M. (2001). Curr. Opin. Colloid Interface Sci. 6, 132–139. Web of Science CrossRef CAS Google Scholar
Balogh, J., Olsson, U. & Pedersen, J. S. (2007). J. Phys. Chem. B, 111, 682–689. CrossRef PubMed CAS Google Scholar
Banchio, A. J., Nägele, G. & Ferrante, A. (1998). J. Colloid Interface Sci. 208, 487–499. CrossRef CAS PubMed Google Scholar
Bartlett, P. & Ottewill, R. H. (1992). J. Chem. Phys. 96, 3306–3318. CrossRef CAS Web of Science Google Scholar
Baxter, R. J. (1970). J. Chem. Phys. 52, 4559–4562. CrossRef CAS Web of Science Google Scholar
Beenakker, C. W. J. & Mazur, P. (1983). Physica A, 120, 388–410. CrossRef Google Scholar
Beenakker, C. W. J. & Mazur, P. (1984). Physica A, 126, 349–370. CrossRef Web of Science Google Scholar
Beurten, P. van & Vrij, A. (1981). J. Chem. Phys. 74, 2744–2748. Google Scholar
Blum, L. & Stell, G. (1979). J. Chem. Phys. 71, 42–46. CrossRef CAS Web of Science Google Scholar
Blum, L. & Stell, G. (1980). J. Chem. Phys. 72, 2212. CrossRef Google Scholar
Bohren, C. F. & Huffman, D. R. (2008). Absorption and Scattering of Light by Small Particles. New York: John Wiley. Google Scholar
Botet, R., Kwok, S. & Cabane, B. (2020). J. Appl. Cryst. 53, 1570–1582. CrossRef CAS IUCr Journals Google Scholar
D'Aguanno, B. (1993). Phys. Scr. T49A, 84–88. CAS Google Scholar
D'Aguanno, B. & Klein, R. (1992). Phys. Rev. A, 46, 7652–7656. CAS PubMed Google Scholar
Dekker, F., Kuipers, B. W. M., González García, Á., Tuinier, R. & Philipse, A. P. (2020). J. Colloid Interface Sci. 571, 267–274. CrossRef CAS PubMed Google Scholar
Diaz Maier, J. & Wagner, J. (2024). Soft Matter, 20, 1309–1319. CrossRef CAS PubMed Google Scholar
Di Cola, E., Moussaïd, A., Sztucki, M., Narayanan, T. & Zaccarelli, E. (2009). J. Chem. Phys. 131, 144903. CrossRef PubMed Google Scholar
Flory, P. J. (1936). J. Am. Chem. Soc. 58, 1877–1885. CrossRef CAS Google Scholar
Frenkel, D., Vos, R. J., de Kruif, C. G. & Vrij, A. (1986). J. Chem. Phys. 84, 4625–4630. CrossRef CAS Web of Science Google Scholar
Genz, U. & Klein, R. (1991). Physica A, 171, 26–42. CrossRef Google Scholar
Ginoza, M. & Yasutomi, M. (1999). J. Phys. Soc. Jpn, 68, 2292–2297. CrossRef CAS Google Scholar
Glatter, O. (2018). Scattering Methods and their Application in Colloid and Interface Science. Amsterdam: Elsevier. Google Scholar
Greene, D. G., Ferraro, D. V., Lenhoff, A. M. & Wagner, N. J. (2016). J. Appl. Cryst. 49, 1734–1739. Web of Science CrossRef CAS IUCr Journals Google Scholar
Griffith, W. L., Triolo, R. & Compere, A. L. (1987). Phys. Rev. A, 35, 2200–2206. CrossRef CAS PubMed Web of Science Google Scholar
Hallett, J. E., Grillo, I. & Smith, G. N. (2020). Langmuir, 36, 2071–2081. Web of Science CrossRef CAS PubMed Google Scholar
Hansen, J.-P., Levesque, D. & Zinn-Justin, J. (1991). Liquids, Freezing and Glass Transition. Amsterdam: North Holland. Google Scholar
Hansen, J.-P. & McDonald, I. R. (2013). Theory of Simple Liquids with Applications to Soft Matter. Amsterdam: Elsevier. Google Scholar
Hayter, J. B. & Penfold, J. (1981). Mol. Phys. 42, 109–118. CrossRef CAS Web of Science Google Scholar
Janssen, L. M. C. (2018). Front. Phys. 6, 97. CrossRef Google Scholar
Karg, M., Pich, A., Hellweg, T., Hoare, T., Lyon, L. A., Crassous, J. J., Suzuki, D., Gumerov, R. A., Schneider, S., Potemkin, I. I. & Richtering, W. (2019). Langmuir, 35, 6231–6255. CrossRef CAS PubMed Google Scholar
Kirkwood, J. G. & Boggs, E. M. (1942). J. Chem. Phys. 10, 394–402. CrossRef CAS Google Scholar
Kotlarchyk, M. & Chen, S.-H. (1983). J. Chem. Phys. 79, 2461–2469. CrossRef CAS Web of Science Google Scholar
Krause, R., D'Aguanno, B., Mendez-Alcaraz, J. M., Nägele, G., Klein, R. & Weber, R. (1991). J. Phys. Condens. Matter, 3, 4459–4475. CrossRef CAS Google Scholar
Li, T., Senesi, A. J. & Lee, B. (2016). Chem. Rev. 116, 11128–11180. Web of Science CrossRef CAS PubMed Google Scholar
Liu, Y., Chen, W.-R. & Chen, S.-H. (2005). J. Chem. Phys. 122, 044507. Web of Science CrossRef Google Scholar
Lu, P. J. & Weitz, D. A. (2013). Annu. Rev. Condens. Matter Phys. 4, 217–233. CAS Google Scholar
Mason, T. G., Graves, S. M., Wilking, J. N. & Lin, M. Y. (2006). J. Phys. Chem. B, 110, 22097–22102. CrossRef PubMed CAS Google Scholar
Mohanty, P. S., Nöjd, S., van Gruijthuijsen, K., Crassous, J. J., Obiols-Rabasa, M., Schweins, R., Stradner, A. & Schurtenberger, P. (2017). Sci. Rep. 7, 1487. CrossRef PubMed Google Scholar
Nayeri, M., Zackrisson, M. & Bergenholtz, J. (2009). J. Phys. Chem. B, 113, 8296–8302. CrossRef PubMed CAS Google Scholar
Ninarello, A., Berthier, L. & Coslovich, D. (2017). Phys. Rev. X, 7, 021039. Google Scholar
Olver, F. W. J., Lozier, D. W., Boisvert, R. F. & Clark, C. W. (2010). The NIST Handbook of Mathematical Functions. Cambridge University Press. Google Scholar
Pedersen, J. S. (1994). J. Appl. Cryst. 27, 595–608. CrossRef Web of Science IUCr Journals Google Scholar
Pedersen, J. S. (1997). Adv. Colloid Interface Sci. 70, 171–210. CrossRef CAS Web of Science Google Scholar
Pedersen, J. S. (2001). J. Chem. Phys. 114, 2839–2846. Web of Science CrossRef CAS Google Scholar
Peláez-Fernández, M., Moncho-Jordá, A. & Callejas-Fernández, J. (2011). J. Chem. Phys. 134, 054905. PubMed Google Scholar
Percus, J. K. & Yevick, G. J. (1958). Phys. Rev. 110, 1–13. CrossRef CAS Web of Science Google Scholar
Phalakornkul, J. K., Gast, A. P., Pecora, R., Nägele, G., Ferrante, A., Mandl-Steininger, B. & Klein, R. (1996). Phys. Rev. E, 54, 661–675. CrossRef CAS Web of Science Google Scholar
Pihlajamaa, I., Laudicina, C. C. L. & Janssen, L. M. C. (2023). Phys. Rev. Res. 5, 033120. CrossRef Google Scholar
Pusey, P. N., Fijnaut, H. M. & Vrij, A. (1982). J. Chem. Phys. 77, 4270–4281. CrossRef CAS Web of Science Google Scholar
Robert, A., Wagner, J., Härtl, W., Autenrieth, T. & Grübel, G. (2008). Eur. Phys. J. E, 25, 77–81. Web of Science CrossRef PubMed CAS Google Scholar
Royall, C. P., Poon, W. C. K. & Weeks, E. R. (2013). Soft Matter, 9, 17–27. CrossRef CAS Google Scholar
Salgi, P. & Rajagopalan, R. (1993). Adv. Colloid Interface Sci. 43, 169–288. CrossRef CAS Google Scholar
Scheffold, F. & Mason, T. G. (2009). J. Phys. Condens. Matter, 21, 332102. CrossRef PubMed Google Scholar
Schulz, G. V. (1939). Z. Phys. Chem. 43B, 25–46. CrossRef Google Scholar
Sciortino, F., Mossa, S., Zaccarelli, E. & Tartaglia, P. (2004). Phys. Rev. Lett. 93, 055701. Web of Science CrossRef PubMed Google Scholar
Scotti, A. (2021). Soft Matter, 17, 5548–5559. CrossRef CAS PubMed Google Scholar
Stellbrink, J., Allgaier, J., Richter, D., Moussaid, A., Schofield, A. B., Poon, W. C. K., Pusey, P. N., Lindner, P., Dzubiella, J., Likos, C. N. & Löwen, H. (2002). Appl. Phys. Mater. Sci. Process. 74, s355–s357. CrossRef CAS Google Scholar
Stieger, M., Pedersen, J. S., Lindner, P. & Richtering, W. (2004). Langmuir, 20, 7283–7292. Web of Science CrossRef PubMed CAS Google Scholar
Stieger, M., Richtering, W., Pedersen, J. S. & Lindner, P. (2004). J. Chem. Phys. 120, 6197–6206. Web of Science CrossRef PubMed CAS Google Scholar
Szymusiak, M., Kalkowski, J., Luo, H., Donovan, A. J., Zhang, P., Liu, C., Shang, W., Irving, T., Herrera-Alonso, M. & Liu, Y. (2017). ACS Macro Lett., 6, 1005–1012. CrossRef CAS PubMed Google Scholar
Tomchuk, O. V., Bulavin, L. A., Aksenov, V. L., Garamus, V. M., Ivankov, O. I., Vul', A. Y., Dideikin, A. T. & Avdeev, M. V. (2014). J. Appl. Cryst. 47, 642–653. CrossRef CAS IUCr Journals Google Scholar
Voigtmann, T. (2003). Mode Coupling Theory of the Glass Transition in Binary Mixtures, Technische Universität München. Google Scholar
Vrij, A. (1978). J. Chem. Phys. 69, 1742–1747. CrossRef CAS Web of Science Google Scholar
Vrij, A. (1979). J. Chem. Phys. 71, 3267–3270. CrossRef CAS Web of Science Google Scholar
Wagner, J. (2004). J. Appl. Cryst. 37, 750–756. Web of Science CrossRef CAS IUCr Journals Google Scholar
Wagner, J. (2012). J. Appl. Cryst. 45, 513–516. CrossRef CAS IUCr Journals Google Scholar
Widom, B. (1967). Science, 157, 375–382. CrossRef PubMed CAS Google Scholar
Zaccarelli, E., Liddle, S. M. & Poon, W. C. K. (2015). Soft Matter, 11, 324–330. CrossRef CAS PubMed Google Scholar
Zackrisson, M., Stradner, A., Schurtenberger, P. & Bergenholtz, J. (2005). Langmuir, 21, 10835–10845. Web of Science CrossRef PubMed CAS Google Scholar
Zhou, B., Gasser, U. & Fernandez-Nieves, A. (2023). Phys. Rev. E, 108, 054604. CrossRef PubMed 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.