Abstract
Anisotropic dynamics on the colloidal length scale is ubiquitous in nature. Of particular interest is the dynamics of systems approaching a kinetically arrested state. The failure of classical techniques for investigating the dynamics of highly turbid suspensions has contributed toward the limited experimental information available up until now. Exploiting the recent developments in the technique of differential dynamic microscopy (DDM), we report the first experimental study of the anisotropic collective dynamics of colloidal ellipsoids with a magnetic hematite core over a wide concentration range approaching kinetic arrest. In addition, we have investigated the effect of an external magnetic field on the resulting anisotropic collective diffusion. We combine DDM with small-angle x-ray scattering and rheological measurements to locate the glass transition and to relate the collective short- and long-time diffusion coefficients to the structural correlations and the evolution of the zero shear viscosity as the system approaches an arrested state.
INTRODUCTION
Unlike their spherical counterparts, anisotropic particles exhibit a richer phase behavior with enhanced complexity due to their orientational degrees of freedom. In addition to the usual gas, liquid, crystal, and glass phases, anisotropic particles such as rods are known to exhibit liquid crystalline phases. While most of the existing studies address structural properties, less is known about the dynamic behavior of anisotropic particles. The most well-studied systems in this regard are perhaps colloidal rods and ellipsoids (1–9). However, because most of these studies involved particle tracking, they are therefore limited to quasi-2D (two dimension) and to self-diffusion. In bulk, several studies combining polarized and depolarized light scattering report the average translational and rotational diffusion coefficients (10–15) of different anisotropic particles. In these systems, multiple scattering and absorption of light at high concentrations make the measurement not only challenging but also sometimes impossible. As a result, most of the studies were carried out in the dilute regime, and a systematic description of the anisotropic dynamics approaching the ordered or arrested glassy phase is still lacking for bulk systems.
Theoretical studies on ellipsoidal particles predict a kinetic phase diagram as a function of packing fraction ϕ and aspect ratio ρ (16, 17). Depending on the aspect ratio, three types of glass transitions have been predicted, with the first one being the conventional cage driven one that is also found for spherical particles, while the other two have their origins related to the orientational degrees of freedom. For nonspherical particles, a glass phase is predicted consisting of nematic domains, where the interdomain orientations build an orientational glass. Experimentally, such an orientational glass has only been reported for a system of very long and thin rods (fd viruses) by Kang and Dhont (18) and Kang (19). Contrary to the prediction by Schilling et al. (16, 17), particle arrest and the freezing of the nematic domain texture were found to occur at the same concentration. The experiments thus reveal that for these long rods, there is a single glass transition, and both particle dynamics and the nematic domain texture dynamics arrest simultaneously. It is also important to point out that while the particles used by Kang et al. resemble rods rather than ellipsoids, at these large axial ratios, we would not expect notable differences in their behavior. The third type of glass transition is predicted for nearly spherical ellipsoids where the orientational degrees of freedom with odd parity flip and freeze independently from the positions. Recent simulations explored the dynamics around the glass transitions in 3D using ellipsoids of smaller aspect ratios (ρ = 1.25 prolate and 0.8 oblate). Consequently, the two-step glass transition with the orientational glass was not observed (20).
Despite a number of theoretical and simulation studies, opacity or turbidity apparently acts as a road block in the direct experimental investigations of concentrated anisotropic particles in 3D. This obstacle can sometimes be overcome by changing the solvent and matching the index of refraction of the particles, and in turn, it affects the interparticle interactions, which is not desirable. Differential dynamic microscopy (DDM) is a recently developed high-throughput technique that can overcome this limitation (21). In DDM, a time series of digital video images is acquired in bright-field or phase-contrast mode using a fast camera fitted to an ordinary optical microscope. The averaged power spectrum of the difference images is used to calculate the intermediate scattering function (ISF), which describes the dynamics of the system. DDM thus accesses the same dynamical quantity as dynamic light scattering (DLS), but in contrast to DLS, can also be applied to absorbing samples and is less affected by multiple scattering.
DDM has already been successfully used (22) to characterize the anisotropic dynamics of dilute colloidal ellipsoids over a wide range of scattering vectors q. We have exploited the advantage of DDM of not being limited by the opacity of the samples to explore the anisotropic dynamics of a bulk colloidal system in the collective regime. In this article, we report a detailed experimental study of the evolution of collective anisotropic dynamics of ellipsoidal colloids of aspect ratio ρ = 3.76 as a function of volume fraction. Through a combination of DDM and shear rheometry, we demonstrate the existence of a glass transition at a volume fraction of ϕ ≈ 0.52. DDM reveals several regimes for the collective dynamics in our ellipsoidal model system. At low concentrations, density fluctuations relax via a single decay process, which can be characterized by a collective diffusion coefficient
Being made up of hematite cores coated with silica shells, these particles are responsive toward an external magnetic field, which helps them to align along a particular direction. We can thus use a magnetic field to investigate the effects of a reduction in the rotational degrees of freedom on collective diffusion and the location of the glass transition. Although both
Theoretical background: Anisotropic diffusion of oriented prolate particles
For a dilute dispersion of noninteracting prolate ellipsoids with semi-long and short axes, a and b, the isotropic translational diffusion coefficient, Diso, can be expressed as the average over the diffusion coefficients along the different axes (26)
In the presence of an external magnetic field B, the individual silica-coated hematite (SCH) particles align with their magnetic moments μ being parallel to B. Because for hematite particles μ is perpendicular to the long axis of the ellipsoid, the particle orientation a is restricted to the plane perpendicular to B (27). However, a can still rotate around the direction of B (Fig. 1A) (27–29). As a result, the diffusion of the particles becomes anisotropic in nature (22). For noninteracting particles, one would expect that with increasing field, an alignment of the individual particles along with a slowing down of the dynamics measured parallel to B develops. The diffusion coefficient decreases from Diso at zero field to a lower value Db in the limit of high field when the particles are fully aligned. In contrast, the dynamics perpendicular to B is expected to increase from Diso at zero field to (Da + Db)/2 at high field. The field-dependent diffusion coefficients, measured parallel and perpendicular to B, are given by (30)
(A) Schematic of arrangement of ellipsoids in the presence of an external magnetic field with their shorter axes being parallel to the field direction. (B) Assortment of varied packing fractions of SCH suspension, which are confined between two glass coverslips separated by a spacer of width 120 μm. The opacity of the samples increases as the volume fraction ϕ increases, as illustrated with the black lines on a paper background positioned behind the samples, which become invisible for volume fractions above ϕ = 0.004 (photo credit: Antara Pal, Lund University).
RESULTS AND DISCUSSION
In DDM, the intensity fluctuations are analyzed from the spatial Fourier transform of a sequence of images to obtain the power spectrum, g(q, τ). g(q, τ) is directly related to the ISF, f (q, τ), where q is the wave vector of the fluctuations being probed and τ is the delay time (21, 31–33). The ISF, which represents the density-density time correlation function, is a measure of the dynamics on the spatial scale 2π/q. By analyzing the ISF in different azimuthal directions with respect to the applied field B, we can determine the effective translational diffusion coefficients, D(B, θ), of the particles as a function of B. In the version of DDM we implement (22, 32, 33), the incident light is completely incoherent, and the Fourier transform of the intensity correlation function of the images is directly related to the ISF by
In the current study, we have explored the dynamics of SCH particles over a wide range of volume fraction, ϕ, starting from ϕ = 0.0002 to ϕ = 0.51 (details of particles synthesis and characterization as well as the volume fraction calculation are provided in Materials and Methods). Figure 1B illustrates the increased turbidity of the dispersion with ϕ, leading to almost opaque samples above ϕ = 0.04. As a result, none of the traditional techniques could be used to study their dynamics but DDM. Depending on the nature of g(q, τ), the whole concentration range can be divided into three regimes: dilute, semi-dilute, and concentrated.
Dilute regime: 0.002 ≤ ϕ < 0.04
A dilute SCH dispersion (ϕ = 0.0002) was first examined following the methodology described by Reufer et al. (22). At these concentrations, the particles are freely diffusive in nature, and the corresponding ISF along a particular direction can be described by a single exponential relaxation mode as
(Top) Dilute sample of volume fraction ϕ = 0.002: (A) Variation of g(q, τ) as a function of delay time, τ, at different wave vectors q (2 μm−1 < q < 5 μm−1) at an azimuthal angle θ = −10°. Symbols show the experimental data, while the lines represent the fits with Eq. 5. q increases from blue to red. (B) Variation of diffusion coefficients, D(B, θ), as a function of azimuthal angle θ at different magnetic fields, which increases from blue (7.9 mT) to red (340 mT). The external magnetic field was applied at an angle θ ∼ 25°. (C) Variation of diffusion coefficients along, D∥, and perpendicular, D⊥, to the applied magnetic field (symbols) along with the fitting with Eq. 2 (solid line). Diso is also shown. (Bottom) Semi-dilute sample of volume fraction ϕ = 0.041: (D) Variation of g(q, τ) as a function of delay time, τ, at different wave vectors q (2 μm−1 < q < 5 μm−1). Symbols show the experimental data, while the lines represent the fit with Eq. 6. q increases from blue to red. (E) Variation of diffusion coefficients, D(θ), as a function of azimuthal angle θ at different magnetic fields, which increases from blue to red. (F) Variation of diffusion coefficients along, D∥, and perpendicular, D⊥, to the applied magnetic field.
Upon increasing the amplitude of the magnetic field, the alignment of the particles perpendicular to the field is promoted, which is manifested by the stronger decoupling of the diffusion coefficient along and perpendicular to the field (Fig. 2C). The particles were found to start to align at about 15 mT and were almost fully aligned for B > 300 mT. A sinusoidal shape was observed for D(B, θ) as a function of θ that can be well described by Eq. 3 (Fig. 2B), which allows us to obtain the values of D∥ and D⊥ as a function of B. According to the theoretical predictions, D⊥(B) increases and D∥(B) decreases with field satisfying Eq. 2, which can be observed in Fig. 2C.
In addition, the measurements with decreasing fields superimposed within experimental error with those obtained for increasing fields, confirming that the particles get reversibly aligned at these concentrations (fig. S2). By fitting Eq. 2 to the observed D∥ and D⊥ (Fig. 2C), Da and Db were determined to be Da = 2.94 μm2/s and Db = 1.77 μm2/s. The expected value for Da and Db based on the size of the particles are 2.95 and 2.32 μm2/s, respectively. The origin of the discrepancy between the observed and expected value for Db is discussed by Reufer et al. (22). These authors have argued that for comparable systems, particle sizes obtained by transmission electron microscopy (TEM) often suffer from the chosen threshold to find the particle contours and that this introduces a systematic error in the calculation of the expected D values. Moreover, the calculations assume perfect prolate shape, and the small geometrical imperfections at the tips of the ellipsoids also contribute to this discrepancy between measured and calculated values. Similar results were found up to ϕ = 0.04.
Semi-dilute regime: 0.04 < ϕ < 0.21
As the volume fraction reaches a value around ϕ = 0.04, the particles start to noticeably interact with each other, which, in turn, affects their dynamics. The ISF along a particular direction, f(q, τ), can no longer be described by a single exponential but with a stretched exponential
Concentrated regime: 0.21 < ϕ < 0.51
In the concentrated regime, a second slow relaxation mode starts to develop in g(q, τ) (Fig. 3A). This new relaxation mode is associated with the structural rearrangements of the cages formed due to the crowding of the particles, while the fast mode represents the rattling inside the cages (17). At this concentration range, f (q, τ) can be well expressed by a double stretched exponential function representing two different relaxation modes with two different characteristic time scales, τ1 and τ2, as
Concentrated sample of volume fraction ϕ = 0.48: (A) Variation of g(q, τ) as a function of delay time, τ, at different wave vectors q (2 μm−1 < q < 5 μm−1). Symbols show the experimental data, while the lines represent the fit with the equation (Eq. 7). q increases from blue to red. (B) Variation of short-time diffusion coefficients,
With increasing field,
The average diffusion coefficients can be obtained by taking an azimuthal average of the diffusion coefficients. At the lowest field, the normalized average short-time diffusion coefficient,
(A) Variation of normalized average short time diffusion coefficient (red circles) and 1/S(qDDM) at qDDM = 10 μm−1 (black squares) with ϕ at the lowest filed. (B) Variation of normalized average long-time diffusion coefficient with ϕ (black squares) at the lowest field; magenta circles show experimentally observed variation of inverse of normalized zero shear viscosity with ϕ, while the dashed orange and blue lines represent the fit with Quemada and VFT models, respectively. Variation of normalized long-time diffusion coefficients with ϕ, along parallel (blue diamond) and perpendicular (green triangle) to the field at high (340 mT) magnetic field.
The normalized average long-time collective diffusion coefficient,
In the presence of an external field, the decoupling of the diffusion coefficients between parallel and perpendicular directions increases with ϕ (Fig. 5). The decoupling behavior for the short-time dynamics can be quantified with the help of a decoupling parameter
Observed variation of the decoupling parameter
In the absence of a magnetic field, simulation (17) has predicted the formation of domains with local nematic order to explain the slowing down mechanism for aspect ratios ρ > 2. However, at low magnetic fields, the SAXS experiments provide no evidence for any nematic ordering as demonstrated with the 2D diffraction patterns shown in (Fig. 6A). This changes markedly when going to larger field strength, where a corresponding nematic order becomes apparent (Fig. 6B). The fact that we find no evidence for the formation of nematic domains in the absence of a field indicates that the slowing down of long-time diffusion and the occurrence of an arrest transition are the results of cage formation only. The question then remains whether this scenario changes in the presence of a high field, where the SAXS experiments reveal the formation of nematic domains. At this point, it is worth looking more into the volume fraction dependence of the anisotropic diffusion coefficients perpendicular and parallel to the field. Their normalized values
2D scattering pattern for ϕ = 0.3 at (A) 7 mT and (B) 300 mT. Panels (C) and (D) represent two models that are based on the diffraction patterns to understand the local structures at low and high fields.
CONCLUSION
We have studied the field-induced anisotropic collective dynamics of colloidal ellipsoids over a broad concentration range by exploiting the recently developed technique DDM. The chosen model system, hematite-silica core-shell particles, has a number of advantages. It not only provides us with a simple synthesis route toward quite monodisperse prolate particles with tunable axial ratios but also allows us to modify the rotational degrees of freedom using external magnetic fields. However, the very high scattering cross section for visible light combined with a strong absorption has previously made it impossible to use the traditional optical tools such as DLS or confocal microscopy to study diffusion at high concentrations and thus characterize the arrest transition expected to occur at high packing fractions. Our work now demonstrates that DDM can be extended to study samples with anisotropic dynamics in the limit of high turbidity and high absorption.
Although the literature is replete with simulations concerning the glass transition for ellipsoidal particles, there has been an almost complete lack of experimental information on these systems. The only study we are aware of for the ellipsoidal glass transition in bulk is based on rheological measurements (45) for an axial ratio ρ = 4.8, and no information on collective dynamics and complementary structural information that would provide further insight into the underlying arrest scenario has so far been presented. It is for this reason that we have performed a systematic study of the structural and dynamic properties of a model ellipsoid system using a combination of DDM, SAXS, and rheology. This has allowed us to study collective diffusion in different concentration regimes and particularly look for the typical signatures of an arrest transition in the ISF measured by DDM. We have characterized both short- and long-time collective diffusion upon approaching kinetic arrest. In the absence of an applied magnetic field, we found that the inverse reduced zero shear viscosity and the long-time collective diffusion coefficient are following the same characteristic scaling relation that has previously been used to reproduce these quantities in hard sphere systems. Moreover, we have been able to demonstrate that the arrest transition for these ellipsoids is not linked to the formation of nematic domains, quite in contrast to previous computer simulations of ellipsoids with comparable axial ratios (17). Furthermore, we have also shown that the arrest mechanism appears to be similar in the presence and absence of the external field and linked to caging and not influenced by the field-induced formation of nematic domains at volume fractions significantly smaller than ϕmax ≈ 0.52.
The combination of DDM and SAXS in the presence of an external magnetic field has shown high potential for future studies. While more work will be required using additional complementary techniques, our results on the field dependence of the anisotropic short- and long-time collective diffusion and the average structure indicate no change in the dominating arrest mechanism with increasing field strength, despite the fact that nematic domains form at high field strengths already before arrest. Here, future experiments using x-ray photon correlation spectroscopy as a function of field strength would allow us to probe the anisotropic short- and long-time collective dynamics at much shorter characteristic length scales of order of the nearest-neighbor distance and using very small scattering volumes (47). This should then allow us to investigate the dynamics of cage relaxation in the vicinity of the arrest transition and its dependence on a field-induced local order and the formation of nematic domains. Moreover, a systematic investigation of the structural and dynamic properties of concentrated suspensions of prolate ellipsoids when approaching kinetic arrest and the underlying arrest mechanisms as a function of axial ratio is needed and will allow us to test the existing literature on simulations.
Until recently, such an endeavor would have been perceived as very difficult experimentally due to a lack of available model systems when using the traditional experimental tools such as DLS or confocal microscopy. However, our study has shown that the type of anisotropic systems that can be used to study the colloidal glass transition can now be greatly extended, as turbidity becomes a minor issue when using DDM. This opens new avenues for exploring the dynamics of strongly interacting colloidal systems with a large variety of shapes and chemical compositions. We believe that our results will inspire future experimental studies on concentrated anisotropic particles and provide the basis for comparison with theory and simulation.
MATERIALS AND METHODS
Synthesis of the SCH
Spindle Fe2O3 hematites were first synthesized in water following the approach described by Ocaña et al. (48) and further coated with a thin silica layer in ethanol using the method developed by Graf et al. (49). After purification by repeated centrifugation/redispersion cycles in water, the SCHs were kept in water as a stock dispersion of 2.5 weight %. Details of the synthesis and characterization of similar particles can be found elsewhere (28, 29).
Characterization of SCHs
The size and shape characterization of the ellipsoidal colloids were carried out using TEM (TEM-CM100 microscope from Philips operating at 100 keV). Particle size distributions were calculated by measuring at least 100 particles from TEM images using ImageJ. The average particle semi-long, a, and semi-short, b, axes are found to be a = 148.5 nm and b = 39.5 nm, leading to an aspect ratio of ρ = 3.76.
Sample preparation
The concentrated SCH dispersions were prepared by centrifuging the stock dispersion to the desired weight fraction, whereas more dilute dispersions were obtained by diluting more concentrated samples. To ensure the proper redispersion, the dense samples were sonicated for 24 hours in a thermostated sonication bath. DDM samples were prepared by placing a 9-μl drop of the sample between two coverslips separated by a 120-μm spacer. Each sample was sonicated for about 2 hours before the preparation of the DDM sample cell. To avoid an eventual drift of the denser and more viscous dispersions, cells filled with dense samples were let to relax typically for 12 to 24 hours in a water-saturated atmosphere.
To convert the weight fraction into volume fraction, the density of a single core-shell ellipsoidal particle was determined with two pair of semiaxes, a and b, representing the long and short axis, respectively. One pair is related with the core dimensions, acore and bcore, while the other is related with the whole particle dimension, apart and bpart. The core is composed of hematite with a density, ρhem, and the shell is composed of silica with a density, ρsil. Because both silica and hematite are porous materials, it is necessary to include a porosity contribution in terms of volume fraction with respect to the overall particle volume. The porosity in the core
Considering the densities, ρsil = 2.06 g cm−3, ρhem = 5.26 g cm−3, and ρwtr = 1.00 g cm−3 (28), and the dimensions from TEM as acore = 131 nm, bcore = 24.5 nm, apart = 148.5 nm, and bpart = 39.5 nm, a particle density (ρpart) of 2.64 g cm−3 can be observed. The following expression can be used to calculate the volume fraction based on the densities of the particles and the solvent, and the weight fraction (χpart) for each sample
Experimental technique I: DDM
The DDM setup consists of an inverted optical microscope (Nikon Eclipse TE200) equipped with a fast camera (Mikrotron). The samples were imaged with a 40× objective in phase contrast mode. The sample stage was modified to install a solid core electromagnet connected to a dc power source (0 to 2.7 A). To prevent an excessive heating of the coils, the electromagnet was water-cooled using an external thermostat. The magnetic field was measured at the sample position using a tesla meter (TL Atomic SMS 102), and the applied magnetic field in the sample was derived from the calibration of the measured magnetic field to the applied current. The measurements were performed with increasing and then decreasing magnetic fields from 7.9 to 380.0 mT, 7.9 mT being the remanent magnetic field for the electromagnet. The temperature of the sample was maintained with an accuracy of 0. 1oC using a homemade sample holder connected to another thermostat. Depending on the concentration of the samples, the data acquisition mode has been varied from 200 to 500 frames/s.
Variation of diffusion coefficient as a function of q
In the range 2 μm−1 < q < 5 μm−1, the diffusion coefficient is independent of q. In fig. S1, we have plotted f(q, τ) as a function of τ. q2 for different q values, and the collapsing of all of them indicates the same fact.
Reversibility of diffusion coefficients
In the dilute regime, the diffusion coefficients measured during increasing the field almost superimposed with that measured while decreasing the field, which is shown in fig. S2.
Variation of stretching exponent β as a function of ϕ
In the semi-dilute and concentrated regimes, we have used single and double stretched exponential function(s) to describe the ISF. We find the stretching parameter β to be less than 1 (fig. S3). This indicates the presence of more than one single relaxation processes, although there is no visual appearance of that in g(q, τ) (Fig. 2D). The presence of more than one single relaxation processes becomes evident once the need to use double stretched exponential functions for concentrated regime arises. In this scenario, we observed the value of β1 abruptly acquiring a higher value as can be observed in fig. S3.
Experimental technique II: SAXS
SAXS experiments were carried out at the cSAXS beamline of the Swiss Light Source (Paul Scherrer Institute, Switzerland) with a sample-detector distance of 7.160 m and an energy of 8.5 keV (corresponding to a wavelength of λ = 0.1459 nm). The scattering patterns were recorded on a single-photon counting PILATUS detector (1475 pixels × 1679 pixels). The samples with different concentrations were sealed in glass capillaries of 1 mm diameter and 0.1 mm wall thickness and investigated at room temperature. A horizontal magnetic field was produced with a water-thermostated electromagnet such that the field lines were always normal to the incident beam. The structure factor was calculated by dividing the scattered intensity I with the experimentally obtained form factor F(q) (fig. S4). The variation of 1/S(q) for different concentrations (shown in Fig. 4A) has been calculated at qDDM = 10 μm−1 (shown by dashed vertical line).
Experimental technique III: Rheology
The rheological measurement was done by an ARES-RFS (rheometrics fluid spectrometer). The shear rate (
SUPPLEMENTARY MATERIALS
Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/6/3/eaaw9733/DC1
Fig. S1. Variation of f(q, τ) as a function of τ. q2 for different q’s, which are collapsing together (ϕ = 0.0002 and azimuthal angle θ = −10∘) as the diffusion coefficient is independent of q in 2 μm−1 < q < 5 μm−1 range.
Fig. S2. Variation of diffusion coefficients at ϕ = 0.0002.
Fig. S3. Variation of stretching exponent β as a function of ϕ from semi-dilute and concentrated regime.
Fig. S4. Variation of the structure factor at various ϕ (0.004 < ϕ < 0.3).
Fig. S5. Symbols represent the variation of viscosity as a function of shear rates for different volume fractions, whereas the continuous lines show the fits of the experimental data with the Cross model.
This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.
REFERENCES AND NOTES
- Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).