Order-to-Disorder Transition in Ring-Shaped Colloidal Stains

Phys. Rev. Lett. 107, 085502
Order-to-disorder transition in the particle stain left by an evaporating drop. (a) A 3μL sessile water droplet evaporates from a glass substrate. (b) Ring-shaped stain of red particles (the coffee stain) left on the substrate after evaporation. (c) A close-up by an optical microscope of the bottom layer of the stain, taken from the white square in (b), shows that the outermost lines of the stain (left) have an ordered, crystalline structure. Towards the center of the drop (right), a transition to a disordered particle arrangement is observed. (d) A top view of the ring stain, taken from the red square in (c) with a scanning electron microscope (SEM), shows that the first lines of particles (left) are arranged in a hexagonal array, while the next lines (brighter in the image) are arranged as square, followed by again hexagonal array.Order-to-disorder transition in the particle stain left by an evaporating drop. (a) A 3μL sessile water droplet evaporates from a glass substrate. (b) Ring-shaped stain of red particles (the coffee stain) left on the substrate after evaporation. (c)... Show more

Ordered arrays of colloidal particles present important characteristics for fields as photonics [1] and biotechnology [2]. A relatively simple approach to deposit particles onto a substrate is by evaporation of a colloidal dispersion droplet. In droplets with pinned contact lines, evaporation gives rise to the so-called coffee-stain effect [3–8]: a ringlike stain remains on the substrate once the liquid has evaporated. The mechanism behind this ring-stain formation is explained in the pioneering work of Deegan et al. [3]: In an evaporating drop with an immobile contact line, a capillary flow is generated to replenish the liquid that has evaporated from the edges. This flow drags particles towards the contact line, forming the ring-shaped stain in particle suspensions such as coffee.

This deposition mechanism has been successfully used to generate colloidal crystals in three dimensions [9]. Furthermore, techniques have been developed to control the patterns by permitting evaporation through specially designed masks [10] or by controlling the wettability [8,11]. Also, large-scale colloidal crystals have been observed for particle sizes below 100 nm [12], but some displayed instabilities that gave rise to concentric rings [13], stripes, chevrons [14], and cracks [15]. At present, there is no way of predicting beforehand the structure of the deposit remaining after evaporation, due to a lack of understanding of the basic mechanisms. However, the processes driving colloidal crystallization are crucial for controlling the mechanical and the optical properties of the synthesized materials.

Analysis of the particle ordering in the stain. (a) Overview of the different patterns observed in the ring stain. The position of the contact line is on the left side of the figure (the first layers of particles are not shown). (b1)–(b3) A close-up reveals the sequence of patterns in the stain: square packing close to the contact line (b1), followed by hexagonal packing (b2), and finally disordered packing (b3). The Voronoi cells belonging to the particles are shown by blue lines. (c) The area of the Voronoi cells plotted versus the distance x from the contact line. The order-to-disorder transition at x=xc is determined by the increase in the mean Voronoi area and the dispersion, due to the lack of order. This transition is indicated by the gray (red) dotted line.Analysis of the particle ordering in the stain. (a) Overview of the different patterns observed in the ring stain. The position of the contact line is on the left side of the figure (the first layers of particles are not shown). (b1)–(b3) A close-up ... Show more

In this Letter, we reveal the conditions leading to colloidal crystallization in a paradigmatic system: a 3-μL sessile water droplet containing colloidal particles evaporating from a smooth glass slide [Fig. 1(a)]. The colloids consist of almost monodisperse red-fluorescent polystyrene particles with particle diameters in the range of 0.52μm. The droplet is gently deposited onto a glass, and is left to evaporate under controlled temperature and humidity. The contact line remains pinned during the entire experiment. Our aim is to analyze the structure of the particle arrangement in the remaining coffee stain. It turns out that the dynamics of particle deposition is crucial for the understanding of the final structure. Therefore, the entire evaporative process is filmed simultaneously from the side and from below [16]. This allows us to follow the evolution of the drop geometry, the formation of the deposit, and the dynamics of particle migration. The particles are filmed from below with a high sensitivity digital camera (PCO Imaging) connected to an inverted microscope. An integrated standard halogen light source is used to illuminate the particles; the emitted light is filtered before it reaches the camera. The camera is operated at a typical frame rate of 1 frame per second, which is sufficient to capture the movement of the particles. To determine the velocity field in the measurement area from the captured images, a custom-made MATLAB microparticle image velocimetry ( μPIV) code was used. The side-view images are used to calculate contact angle, volume, radius, and height of the drop at every instant with a custom-made MATLAB code. Side-view images are taken with a long-distance microscope and a CCD camera (Lumenera Corp.), synchronized with the visualization of the particles by a common trigger.

When the droplet has completely evaporated, most of the colloidal particles have aggregated into a ring-shaped coffee stain [Fig. 1(b)]. Zooming in on the stain, we see that the particle arrangement is not homogeneous [Fig. 1(c)]; there is a remarkable transition from a crystalline arrangement to a disordered phase. A top view of the stain [Fig. 1(d)] shows a sequence of particle arrangements within the ordered phase, starting from hexagonal packing, followed by square packing, and again hexagonal packing. The particle arrangement is analyzed in more detail in Fig. 2. Close to the contact line, we find a square arrangement of particles [Fig. 2(b1)], followed by a hexagonal array [Fig. 2(b2)], and finally a disordered phase [Fig. 2(b3)]. The same trend has been observed for different particle sizes and concentrations. A well-known method to quantitatively assess the different packing structures is to calculate the Voronoi areas [17] of the particles, as shown in Figs. 2(b1), 2(b2), and 2(b3). In Fig. 2(c), we identify the Voronoi areas corresponding to the different cell arrangements: the disordered phase exhibits a larger average cell area as well as a wider dispersion, compared to the ordered phase. This reveals the sharp transition from order to disorder.

Origin of the order-to-disorder transition. Plot of the radial particle velocity versus time, measured at a distance of 174μm from the contact line. The fluid velocity predicted by the model (1) is plotted as a solid line. Both model and experiment show the dramatic velocity increase at the end of the droplet’s life (the rush hour). The inset shows the direction of the radial velocity in the droplet. The time tc at which the order-to-disorder transition occurs is determined from the videos taken for the μPIV measurements; the corresponding critical velocity uc is read from this graph (indicated by the black line in the figure).Origin of the order-to-disorder transition. Plot of the radial particle velocity versus time, measured at a distance of 174μm from the contact line. The fluid velocity predicted by the model (1) is plotted as a solid line. Both model and experiment... Show more

The order-to-disorder transition can be traced back to the hydrodynamics within the drop. Figure 3 shows the radial velocity measured by μPIV of the layer of particles closest to the glass slide, in a region close to the contact line. The particle velocity increases dramatically in the last moments of the droplet’s life. We refer to this sudden change in speed as “rush hour.” The particles that arrive early, at a low deposition speed, form an ordered (square or hexagonal) structure. In contrast, particles that arrive during rush hour have a high speed and form a jammed, disordered phase.

The rush-hour behavior can be explained by simple mass-balance considerations. The volume flow towards the contact line inside the drop is driven by the evaporation from the drop surface [3]. It turns out that the rate of evaporation is approximately constant over time [6]. To replenish this evaporated liquid, a continuous volume flow towards the contact line is generated inside the drop. However, the drop height is vanishing during evaporation. This can be characterized by a contact angle θ(t)0 at the final stages of evaporation (see Fig. 3). Hence, the same amount of liquid has to be squeezed through an area which is vanishing, inducing a diverging radial velocity. Close to the contact line, the dependence of the height-averaged radial velocity u¯ on time t and distance from the drop center r is given by [4,6]

The experimentally determined critical velocity uc versus its theoretical prediction (3).

see [16] for details. Here, D*=22DvaΔc/(πρ) represents the driving of the flow, with Dva=24×10-6m2/s the diffusion constant for vapor in air, Δc=1.2×10-2kg/m3 the vapor concentration difference between the drop surface and the surroundings, ρ=998kg/m3 the liquid density, and R the drop radius. The spatial (R-r)-1/2 divergence originates from the diverging evaporative flux from the drop surface [3]. Although this spatial divergence generated a lot of interest [3–8], it is clear from (1) that there is also a temporal divergence because θ(t) vanishes. This temporal divergence provides the key to the rush-hour behavior: the radial velocity blows up towards the end of the droplet’s life.

The experimental velocity measurements are performed at a distance of 2μm above the glass slide, and therefore differ from the theoretical height-averaged velocity. To quantitatively compare the model with the experimental data, a description of the velocity profile is required. Given the low velocities measured inside the evaporating drop (low Reynolds number), the inertial terms in the Navier-Stokes equation can be neglected. Since the droplet height is much smaller than its radius, the thin-film (lubrication) approximation [18] can be used to describe the velocity field. Hence, the velocity u(r,z,t) [16] is given by

with h(r,t)=θ(t)(R2-r2)/(2R) the local droplet height. As shown in Fig. 3, this model quantitatively describes the rush-hour behavior. For the plot of (2) we took R-r=174μm (the distance from the contact line where the μPIV measurements are performed), z=2μm and te the droplet’s lifetime determined from the experiment. This means there are no adjustable parameters in the comparison between the model and the experiments.

We show that the structural change in the stain can be entirely attributed to the rush-hour behavior. The particles that form the ordered phase arrive at the contact line at an early time, when the deposition rate is still low. On the other hand, the particles that belong to the disordered phase arrive during rush hour, at a high deposition rate. The slow, early arriving particles have time to arrange themselves within the stain by Brownian motion, while those coming in rush hour have no time to find an appropriate spot and are quenched into a disordered phase. We can quantify this transition by comparing the time scale on which particles arrange by diffusion to the time scale of particle deposition. The diffusive time scale is tD=dp2/D, where D is the diffusivity of the particles in liquid, calculated by the Einstein relation. The hydrodynamic time scale is th=L/u [where L is the typical distance between the particles, dependent on their concentration in the solution np(t) as L=np-1/3, and u is the particle velocity]. When tD and th become of the same order, a transition is expected. We therefore define the critical transition velocity as

We validated this criterion experimentally for various particle concentrations and sizes. First, we determined the time tc when the ring-shaped stain undergoes the transition [the position where this occurs is xc, as indicated by the red line in Fig. 2(c)]. Next, we used the velocity-time plots (Fig. 3) to find the experimental value of u(tc)=uc. As demonstrated in Fig. 4, the experimental critical velocity perfectly follows the proposed scaling. The agreement suggests that the prefactor in (3) is actually of order unity. Our model also explains many results described in the literature [12–15], where large-scale crystalline deposits were observed for particle sizes ranging from 6 to almost 50 nm, in hindsight. In this range of particle sizes, the critical velocity required for a disordered phase (3) is unreachable within the droplet, which explains the absence of a disordered phase.

Apart from the order-to-disorder transition, there are also transitions within the ordered phase, from square to hexagonal packing and vice versa [Fig. 1(d)]. In an unconfined system, the most efficient packing lattice is the hcp, as conjectured by Kepler [19]. However, in the evaporating drop the particles are confined in a wedge formed by the glass slide on one side and the liquid-air interface with contact angle θ on the other side. It was shown by Pieranski et al. [20,21] that such a confinement indeed leads to a sequence of hexagonal and square packing structures, dependent on the most efficient packing for the available space. Indeed, when a new layer is formed, the confinement by the wedge favors the square packing; see Fig. 1. Away from the step edge, the available space increases and allows for the denser, hexagonal packing structure. We observe a sequence of this type of transitions each time a new layer is added. Similar transitions were also observed by Abkarian et al. [22] in a wedge formed by a moving contact line.

In conclusion, in this Letter we revealed the mechanism triggering crystallization in an evaporating colloidal solution drop with a pinned contact line. The temporal divergence (rush hour) observed in the velocity field inside the drop is responsible for the lack of order in the ring-shaped stains: when the deposition rate is low, particles have time to arrange by Brownian motion and form an ordered phase. During the rush hour, there is no time for such rearrangement and the particles form a disordered phase. Such effect resembles the popular computer game Tetris, where at the beginning of the game slowly falling objects are easily arranged in regular structures, while as the game evolves, the fast falling objects become jammed and disordered. Therefore, the key to generate large-scale colloidal crystals is to avoid the rush-hour behavior, i.e., to avoid small contact angles during evaporation. This can be achieved, for example, by modifying the hydrophobicity of the substrate or by continuously replenishing the evaporated liquid. Large-scale crystalline arrays can also be obtained by tuning the particle size, since nanometer-sized particles diffuse much more efficiently. As a rule of thumb, freely evaporating colloidal solutions with particles within the nanometric range will naturally crystallize, micrometer-sized particles will show both ordered and disordered phases, while larger particles lead to disordered stains


  1. A. Stein and R. Schroden, Curr. Opin. Solid State Mater. Sci. 5, 553 (2001).
  2. L. Zhang, S. Maheshwari, H. Chang, and Y. Zhu, Langmuir 24, 3911 (2008).
  3. R. D. Deegan, O. Bakajin, T. Dupont, G. Huber, S. Nagel, and T. Witten, Nature (London) 389, 827 (1997).
  4. R. D. Deegan, O. Bakajin, T. F. Dupont, G. Huber, S. R. Nagel, and T. A. Witten, Phys. Rev. E 62, 756 (2000).
  5. H. Hu and R. G. Larson, Langmuir 21, 3963 (2005).
  6. Y. O. Popov, Phys. Rev. E 71, 036313 (2005).
  7. H. Bodiguel and J. Leng, Soft Matter 6, 5451 (2010).
  8. H. B. Eral, D. M. Augustine, M. H. G. Duits, and F. Mugele, Soft Matter 7, 4954 (2011).
  9. K. P. Velikov, C. G. Christova, R. P. A. Dullens, and A. van Blaaderen, Science 296, 106 (2002).
  10. D. J. Harris, H. Hu, J. C. Conrad, and J. A. Lewis, Phys. Rev. Lett. 98, 148301 (2007).
  11. F. Fan and K. J. Stebe, Langmuir 20, 3062 (2004).
  12. T. P. Bigioni, X. M. Lin, T. T. Nguyen, E. I. Corwin, T. A. Witten, and H. M. Jaeger, Nature Mater. 5, 265 (2006).
  13. L. Pauchard, F. Parisse, and C. Allain, Phys. Rev. E 59, 3737 (1999).
  14. C. Allain and L. Limat, Phys. Rev. Lett. 74, 2981 (1995).
  15. E. R. Dufresne, E. I. Corwin, N. A. Greenblatt, J. Ashmore, D. Y. Wang, A. D. Dinsmore, J. X. Cheng, X. S. Xie, J. W. Hutchinson, and D. A. Weitz, Phys. Rev. Lett. 91, 224501 (2003).
  16. See Supplemental Material at http://link.aps.org/supplemental/10.1103/PhysRevLett.107.085502 for a video illustrating the rush-hour effect. A document is also included with further experimental details and a more detailed description of the theoretical model employed.
  17. G. Voronoi, J. Reine Angew. Math. 133, 97 (1908).
  18. A. Oron, S. H. Davis, and S. G. Bankoff, Rev. Mod. Phys. 69, 931 (1997).
  19. J. Kepler, Strena Seu de Nive Sexangula (Godefridum Tampach, Francofurti ad Moenum, 1611); [translated as J. Kepler, The Six-Cornered Snowflake (Clarendon, Oxford, 1966)].
  20. P. Pieranski, L. Strzelecki, and B. Pansu, Phys. Rev. Lett. 50, 900 (1983).
  21. B. Pansu, P. Pieranski, and P. Pieranski, J. Phys. (Paris) 45, 331 (1984).
  22. M. Abkarian, J. Nunes, and H. Stone, J. Am. Chem. Soc. 126, 5978 (2004).

About the Authors

Image of Álvaro G. Marín
Image of Hanneke Gelderblom
Image of Detlef Lohse

Detlef Lohse earned his Ph.D. on the theory of turbulence in Marburg, Germany, in 1992. As a postdoc in Chicago and later in Marburg and München he worked on single bubble sonoluminescence. In 1998 he was appointed Chair of Physics of Fluids at the University of Twente, The Netherlands (http://pof.tnw.utwente.nl/). Lohse’s present research subjects are turbulence and multiphase flow, granular matter, and micro- and nanofluidics, employing experimental, theoretical, and numerical methods. Lohse is Associate Editor of and several other journals. He is a Fellow of the American Physical Society and of the Institute of Physics. He was also elected a Member of the German Academy of Science and the Royal Dutch Academy of Science. In 2005 he received the Spinoza Prize.

Image of Jacco H. Snoeijer

Related Articles

Focus: 3D Structure Shrinks When Heated

Focus: 3D Structure Shrinks When Heated

The volume of a star-shaped structure decreases when baked. Combining this technology with more conventional structures could lead to materials that don't expand or contract with temperature changes. Read More »

Synopsis: A Crack in Earth’s Protective Shield

Synopsis: A Crack in Earth’s Protective Shield

Observations with India’s cosmic-ray telescope indicate that Earth’s magnetic field weakened during a 2015 geomagnetic storm, allowing cosmic rays to pass through. Read More »

Synopsis: Time Optimization in Quantum Computing  
Quantum Information

Synopsis: Time Optimization in Quantum Computing  

New experiments find the fastest way to manipulate logic gates with two qubits as inputs. Read More »

More Articles