Inverted Cheerios effect: Liquid drops attract or repel by elasto-capillarity
Solid particles floating at a liquid interface exhibit a long-ranged attraction mediated by surface tension. In the absence of bulk elasticity, this is the dominant lateral interaction of mechanical origin. Here we show that an analogous long-range interaction occurs between adjacent droplets on solid substrates, which crucially relies on a combination of capillarity and bulk elasticity. We experimentally observe the interaction between droplets on soft gels and provide a theoretical framework that quantitatively predicts the migration velocity of the droplets. Remarkably, we find that while on thick substrates the interaction is purely attractive and leads to drop-drop coalescence, for relatively thin substrates a short-range repulsion occurs which prevents the two drops from coming into direct contact. This versatile, new interaction is the liquid-on-solid analogue of the “Cheerios effect”. The effect will strongly influence the condensation and coarsening of drop soft polymer films, and has potential implications for colloidal assembly and in mechanobiology.
The long-ranged interaction between particles trapped at a fluid interface is exploited for the fabrication of microstructured materials via self-assembly and self-patterning tessier2001 (); Furst11 (); Cavallaro13 (); bowden1997 (); Ershov13 () and occurs widely in the natural environment when living organisms or fine particles float on the surface of water loudet2011 (); peruzzo2013 (). In a certain class of capillary interactions the particles deform the interface because of their shape or chemical heterogeneity danov2005 (); botto2012 (); kumar2013 (). In this case the change in interfacial area upon particle-particle approach causes an attractive capillary interaction between the particles. In the so-called Cheerios effect, the interaction between floating objects is mainly due to the change in gravitational potential energy associated to the weight of the particles, which deform the interface while being supported by surface tension Vella2005 (), and the same principle applies when the interface is elastic Maha2015 (). The name “Cheerios effect” is reminiscent of breakfast cereals floating on milk and sticking to each other or to the walls of the breakfast bowl .
Here we consider a situation opposite to that of the Cheerios effect by considering liquid drops deposited on a solid. The solid is sufficiently soft to be deformed by the surface tension of the drops, resulting in a lateral interaction. Recent studies have provided a detailed view of statics of single-drop wetting on deformable surfaces JXWD11 (); Limat12 (); MDSA12b (); Style12 (); LUUKJFM14 (). The length scale over which the substrate is deformed is set by the ratio of the droplet surface tension and the shear modulus . The deformation can be seen as an elasto-capillary meniscus, or “wetting ridge”, around the drop (Fig. 1 A,B). Interestingly, the contact angles at the edge of the drop are governed by Neumann’s law, just as for oil drops floating on water. In contrast to the statics of soft wetting, its dynamics has only been explored recently. New effects such as stick-slip induced by substrate viscoelaticity limat14 (); karpitschka2015 () and droplet migration due to stiffness gradients durotaxisPNAS13 () have been revealed. The possibility that elasto-capillarity induces an interaction between neighboring drops is of major importance for applications such as drop condensation on polymer films Sokuler09 () and self-cleaning surfaces Schellenberger15 (); Rykaczewski14 (); Kim12 (); Subramanyam12 (). The interaction between drops on soft surfaces might also provide insights into the mechanics of cell locomotion Science2005 (); Lo:BPJ2000 (); Zhang:NatComm2014 () and cell-cell interaction Guo:BPJ2006 ().
Here we show experimentally that long-ranged elastic deformations lead to an interaction between neighboring liquid drops on a layer of cross-linked polydimethylsiloxane (PDMS). The layer is sufficiently soft for significant surface tension induced deformations to occur (Fig. 1). The interaction we observe can be thought of as the inversed Cheerios effect, since the roles of the solid and liquid phases are exchanged. Remarkably, the interaction can be either attractive or repulsive, depending on the geometry of the gel. We propose a theoretical derivation of the interaction force from a free energy calculation that self-consistently accounts for the deformability of both the liquid drop and the elastic solid.
I Experiment: attraction versus repulsion
Here, the inverted Cheerios effect is observed with sub-millimeter drops of ethylene glycol on a PDMS gel. The gel is a reticulated polymer formed by polymerizing small multifunctional prepolymers – contrary to hydrogels, there is no liquid phase trapped inside. The low shear modulus of the PDMS gel gives an elasto-capillary length sufficiently large to be measurable in the optical domain.
The interaction between two neighboring liquid drops is quantified by tracking their positions while they are sliding under the effect of gravity along a soft layer held vertically. The interaction can be either attractive (Fig. 1 A) or repulsive (Fig. 1 B): drops on relatively thick gel layers attract each other, while drops on relatively thin layers experience a repulsion.
The drop-drop interaction induces a lateral motion that can be quantified by the horizontal component of the relative droplet velocity, , with the convention that implies repulsion. In Fig. 1 C, we report as a function of the separation , defined as the shortest distance between the surfaces of the drops. The drops () exhibit attraction when sliding down a thick layer (, black curve), while they are repelled on a thin layer (, red curve). is larger at close proximity, signaling an increase in the interaction force. Spontaneous merging occurs where drops encounter direct contact. Importantly, these interactions provide a new mechanism for droplet coarsening (or ordering) by coalescence (or its suppression) that has no counterpart on rigid surfaces.
The interaction force can be inferred from the relative velocities between the drops, based on the effective “drag” due to sliding on a gel. We first calibrate this drag by considering drops that are sufficiently separated, so that they do not experience any mutual interaction. The motion is purely downward and driven only by gravitational force , and inertia is negligible. Figure 1 D shows that the droplet velocity is nonlinear, and approximately scales as . This force-velocity calibration curve is in good agreement with viscoelastic dissipation in the gel, based on which one expects the scaling law karpitschka2015 ():
Here is the rheological exponent that emerges from the scale invariance of the gel network ChambonWinter (); LAL96 (); deGennes1996 (), while is a characteristic timescale. The parameter values and are calibrated in a rheometer (see Methods). Equation (1) is valid for below the characteristic rheological speed, – this is justified here since for the silicone gel, while the reported speeds reach up to . The presence of a large viscoelastic dissipation in the gel exceeds the dissipation within the drop by orders of magnitude, and explains these extremely slow drop velocities observed experimentally CGS96 (); karpitschka2015 (). The force-distance relation for the Inverse Cheerios effect can now be measured directly using the independently calibrated force-velocity relation [Fig. 1 D and (1)]. By monitoring how the trajectories are deflected with respect to the downward motion of the drops, we obtain .
The key result is shown in Fig. 2, where we report the interaction force as a function of distance . Panel A shows experimental data for the attractive force () between drops on thick layers (black dots), together with the theoretical prediction outlined below. Supplementary Video 1 shows an example of attractive drop-drop interaction. The attractive force is of the order of N, which is comparable to both the capillary force-scale and the elastic force-scale . The force decreases for larger distance and its measurable influence was up to .
Figure 2 B shows the interaction force between drops on thin layers. The dominant interaction is now repulsive (), see Supplementary Video 2. Intriguingly, we find that the interaction is not purely repulsive, but also displays an attractive range at very small distance. It is possible to access this range experimentally in case the motion of the individual drops are sufficiently closely aligned (Supplementary Video 3). The “neutral” distance where the interaction force changes sign appears when the separation is comparable to the substrate thickness , suggesting that the key parameter governing whether the drops attract or repel is the thickness of the gel.
Ii Mechanism of interaction: rotation of elastic meniscus
We explain the attraction versus repulsion of neighboring drops by computing the total free energy of drops on gel layers of different thicknesses. In contrast to the normal Cheerios effect, which involves two rigid particles, in the current experiment an additional element of complexity is present: both the droplet and the elastic substrate are deformable, and their shapes will change upon varying the distance . Hence, the interaction force involves both the elastic and the surface tension contributions to the free energy that emerge from self-consistently computed shapes of the drops and elastic deformations.
To reveal the mechanism of interaction, we first consider two dimensional drops, for which the free energy can be written as
The geometry is sketched in Fig. 3 BC, and further details are given in the Supplementary Information. The elastic energy in the entire layer is a functional of the profile describing the shape of the elastic solid: the functional explicitly depends on the layer thickness, and is ultimately responsible for the change from attraction to repulsion. The function represents the shape of the liquid-vapor interface. The integrals in (2) represent the interfacial energies; they depend on the surface tensions , , associated with the liquid-vapor, solid-liquid and solid-vapor interfaces, respectively. Variation with respect to the contact line positions provide the relevant boundary conditions to the problem LUUKJFM14 (). In order to determine the force per unit length characterizing the lateral interaction between the two-dimensional drops, the free energy (2) is minimized under the constraint that the nearest contact lines are separated by a distance . This constraint is imposed by adding the term to the free energy, where is the inner contact line position (see Supplementary Information). With this convention attractive forces correspond to . Equivalently, is the external force necessary to prevent the drop from moving towards each other, yielding an external work when varying the spacing between the drops.
The energy minimization reveals the mechanism of drop-drop interaction: the interaction force appears in the boundary condition for the contact angles,
where the angles are defined in Fig. 3. Equation (3) can be thought of as an “imbalance” of the usual Neumann boundary condition. For a single droplet, the contact angles satisfy Neumann’s law, which is (3) with (Fig. 3 A). On a thick elastic layer, the overall shape of the wetting ridge is of the form Style12 (); LUUKJFM14 ()
where the horizontal scale is set by elasto-capillary length based on the typical solid surface tension . The origin of can be understood from the principle of superposition. Due to the substrate deformation of a single drop, a second drop approaching the first one will see a surface that is locally rotated by an angle . The elastic meniscus near the inner contact line of this approaching drop gets rotated by an angle (Fig. 3 B). Importantly, changes in the liquid angle exhibit a weaker dependence , which for large drops can be ignored. As a consequence, this meniscus rotation induces an imbalance of the surface tension forces according to (3), which for small rotations yields , where follows from the single drop deformation (4). There is no resultant interaction force from the stress below the drop, which, due to deformability of the drop, provides a uniform pressure.
The inverted Cheerios effect is therefore substantially different from the Cheerios effect between two particles floating at the surface of a liquid. Apart from the drop being deformable, we note that the energy driving the interaction is different for the two cases: while the liquid interface shape is determined by the balance between gravity and surface tension in the Cheerios effect, the solid shape is determined by elasto-capillarity in the inverted Cheerios effect. Another difference is the mechanism by which the interaction is mediated. The Cheerios effect is primarily driven by a change in gravitational potential energy which implies a vertical displacement of particles: a heavy particle slides downwards, like a bead on a string, along the deformation created by a neighboring particle Vella2005 (). A similar interaction was recently discussed for rigid cylinders that deform an elastic surface due to gravity Maha2015 (). In contrast, the inverted Cheerios effect discussed here does not involve gravity and can be totally ascribed to elasto-capillary tilting of the solid interfaces – as in Fig. 3 – manifesting the interaction as a force near the contact line.
The rotation of contact angles indeed explains why the drop-drop interaction can be either attractive or repulsive. On a thick substrate, the second drop experiences solid contact angles that are rotated counter clockwise, inducing an attractive force (Fig. 3 B). By contrast, on a thin substrate the elastic deformation induced by the second drop has a non-monotonic profile . This is due to volume conservation: the lifting of the gel near the contact line creates a depression at larger distances (Fig. 3 C). At large distance, the rotation of the contact angles thus changes sign, and, accordingly, the interaction force changes from attractive to repulsive. Naturally, the relevant length scale for this phenomenon is set by the layer thickness .
Iii Three-dimensional theory
The extension of the theory to three dimensions is straightforward and allows for a quantitative comparison with the experiments. For the three dimensional case we compute the shape of the solid numerically, by first numerically solving for the deformation field induced by a single using an axisymmetric elastic Green function Style12 (). Adding a second drop on this deformed surface gives an intricate deformation that is shown in Fig. 4. The imbalance of the Neumann law applies everywhere around the contact line: the background deformation induces a rotation of the solid contact angles around the drop. According to (3), these rotations result into a distribution of force per unit length contact line , where is the radial direction associated with the interacting drop and the azimuthal angle along the contact line (Fig. 4). The resultant interaction force is obtained by integration along the contact line, as (see Supplementary Information). By symmetry, this force is oriented along the line connecting the two drops.
The interaction force obtained by the three-dimensional theory is shown as the red curves in Fig. 2 AB. The theory gives an excellent description of the experimental data without adjustable parameters. The quantitative agreement indicates that the interaction mechanism is indeed caused by the rotation of the elastic meniscus.
In summary, we have shown that liquid drops can exhibit a mutual interaction when deposited on soft surfaces. The interaction is mediated by substrate deformations, and its direction (repulsive versus attractive) can be tuned by the thickness of the layer. The measured force/distance relation is in quantitative agreement with the proposed elasto-capillary theory. The current study reveals that multiple “pinchings” of an elastic layer by localized tractions lead to an interaction having a range comparable to . The key insight is that interaction emerges from the rotation of the elastic surface, providing a generic mechanism that should be applicable to a wide range of objects interacting on soft media.
In biological settings, elasto-capillary interactions may play a role in cell-cell interactions, which are known to be sensitive to substrate stiffness Guo:BPJ2006 (). In addition, the elastic interaction could also play a role in cell-extracellular matrix interactions, as a purely passive force promoting aggregation between anchor points on the surface of adhered biological cells. For example, it has been demonstrated that a characteristic distance of about nm between topographical features enables the clustering of integrins. These transmembrane proteins are responsible for cell adhesion to the surrounding matrix, mediating the formation of strong anchor points when cells adhere to substrates huang2009 (); dalby2014 (). Assuming that the topographical features “pinch” the cell with a force likely comparable to the cell’s cortical tension, which takes values in the range krieg2008 (); tinevez2009 (); fischer2014 (); sliogeryte2014 (), and an elastic modulus of in the physiological range of biological tissues swift2013 (), one predicts a range of interaction consistent with observations.
More generally, substrate-mediated interactions could be dynamically programmed using the responsiveness of many gels to external stimuli (pH, temperature, electric fields). Possible applications range from fog harvesting and cooling to self-cleaning or anti-fouling surfaces, which rely on controlling drop migration and coalescence. The physical mechanisms revealed here, in combination with the fully quantitative elasto-capillary theory, paves the way for new design strategies for smart soft surfaces.
V Material and Methods
The Supplementary Information provides further technical information, the variational derivation underlying Equation (3), and the numerical scheme that lead to the calculations of Fig. 2 and 4. Supplementary Videos 1, 2 and 3 show typical experiments of drop-drop interactions.
Substrate preparation The two prepolymer components (Dow Corning CY52-276 A and B) were mixed in a ratio of 1.3:1 (A:B). Thick elastic layers () were prepared in petri dishes (diameter ). Thin layers () were prepared by spin-coating the gel onto silicon wafers. Thicknesses were determined by color interferometry. See supplement for details on substrate curing & rheology.
Determining the interaction between drops Droplets of Ethylene Glycol () were pipetted onto a small region near the center of the cured substrate. The sample was then mounted vertically so that gravity acts along the surface ( direction, compare Fig. 1 A,B). The droplets were observed in transmission (thick layers) or reflection (thin layers) with collimated illumination, using a telecenric lens (JenMetar 1x) and a digital camera (pco 1200). Images were taken every 10 s. The contours of the droplets were determined by a standard correlation technique.
At large separation, droplets move downward due to gravity. The gravitational force for each droplet is proportional to its volume. The relation between force and velocity is follows the same power law as the rheology, which was explained recently karpitschka2015 ().
Due to their different volumes/velocities, distances between droplets change. Whenever two droplets come close, their trajectories change due to their interaction. Drops on thick substrates (Figure 1 C, black) attract and finally merge. On a rigid surface, these droplets would have passed by each other. Opposite holds for droplets on thin layers (red): the droplets repel each other, which prevents coalescence.
To determine the interaction forces, we first evaluate the velocity vector of each individual droplet. The droplets behave as quasi-stationary, and the total force vector acting on each droplet is aligned with its velocity vector. The magnitude of the total force is derived by the calibration shown in Fig. 1 D). The interaction force is obtained by subtracting the gravitational force vector from the total force vector. Figure 2 A shows data from nine individual droplet pairs, at different times and different locations on the substrate. Panel B) shows data from 18 different droplet pairs. The raw data has been averaged over distance bins, taking the standard deviation as error bar.
Acknowledgements.SK acknowledges financial support from NWO through VIDI Grant No. 11304. AP and JS acknowledge financial support from ERC (the European Research Council) Consolidator Grant No. 616918.
- (1) Tessier PM, Velev OD, Kalambur AT, Lenhoff AM, Rabolt JF, Kaler EW (2001) Structured metallic films for optical and spectroscopic applications via colloidal crystal templating Advanced Materials 13(6):396–400.
- (2) Furst E (2011) Directing colloidal assembly at fluid interfaces Proc. Natl. Acad. Sci. USA 108(52):20853–20854.
- (3) Cavallaro M, Gharbi M, Beller D, Čopar S, Shi Z, Baumgart T, Yang S, Kamien R, Stebe K (2013) Exploiting imperfections in the bulk to direct assembly of surface colloids Proc. Natl. Acad. Sci. USA 110(47):18804–18808.
- (4) Bowden N, Terfort A, Carbeck J, Whitesides G (1997) Self-assembly of mesoscale objects into ordered two-dimensional arrays Science 276(5310):233–235.
- (5) Ershov D, Sprakel J, Appel J, Cohen Stuart M, van der Gucht J (2013) Capillarity-induced ordering of spherical colloids on an interface with anisotropic curvature Proc. Natl. Acad. Sci. USA 110(23):9220–9224.
- (6) Loudet J Pouligny B (2011) How do mosquito eggs self-assemble on the water surface? The European Physical Journal E 34(8):1–17.
- (7) Peruzzo P, Defina A, Nepf HM, Stocker R (2013) Capillary interception of floating particles by surface-piercing vegetation Physical review letters 111(16):164501.
- (8) Danov K, Kralchevsky P, Naydenov B, Brenn G (2005) Interactions between particles with an undulated contact line at a fluid interface: Capillary multipoles of arbitrary order J. Colloid Interface Sci. 287:121–134.
- (9) Botto L, Lewandowski EP, Cavallaro M, Stebe KJ (2012) Capillary interactions between anisotropic particles Soft Matter 8(39):9957–9971.
- (10) Kumar A, Park BJ, Tu F, Lee D (2013) Amphiphilic janus particles at fluid interfaces Soft Matter 9(29):6604–6617.
- (11) Vella D Mahadevan L (2005) The cheerios effect Am J. Phys. 73(9):817–825.
- (12) A. C, L. R, M.K. C, Mahadevan L (2015) Elastic cheerios effect: Self-assembly of cylinders on a soft solid Europhys. Lett. 112:54001.
- (13) Jerison ER, Xu Y, Wilen LA, Dufresne ER (2011) Deformation of an elastic substrate by a three-phase contact line Phys. Rev. Lett. 106(18):186103.
- (14) Limat L (2012) Straight contact lines on a soft, incompressible solid Eur. Phys. J. E 35(12):134.
- (15) Marchand A, Das S, Snoeijer JH, Andreotti B (2012) Contact angles on a soft solid: From young’s law to neumann’s law Phys. Rev. Lett. 109(23):236101.
- (16) Style RW Dufresne ER (2012) Static wetting on deformable substrates, from liquids to soft solids Soft Matter 8(27):7177–7184.
- (17) Lubbers LA, Weijs J, Botto L, Das S, Andreotti B, Snoeijer JH (2014) Drops on soft solids: Free energy and double transition of contact angles J. Fluid. Mech. 747:R1.
- (18) Kajiya T, Brunet P, Royon L, Daerr A, Receveur M, Limat L (2014) A liquid contact line receding on a soft gel surface: dip-coating geometry investigation Soft Matter 10:8888–8895.
- (19) Karpitschka S, Das S, van Gorcum M, Perrin H, Andreotti B, Snoeijer J (2015) Droplets move over viscoelastic substrates by surfing a ridge Nature Commun. 6:7891.
- (20) Style RW, Che Y, Park SJ, Weon BM, Je JH, Hyland C, German GK, Power MP, Wilen LA, Wettlaufer JS, Dufresne ER (2013) Patterning droplets with durotaxis Proc. Natl. Acad. Sci. USA 110(31):12541–12544.
- (21) Sokuler M, Auernhammer G, Roth M, Liu C, Bonacurrso E, Butt HJ (2009) The softer the better: Fast condensation on soft surfaces Langmuir 26(3):1544–1547.
- (22) Schellenberger F, Xie J, Encinas N, Hardy A, Klapper M, et al. (2015) Direct observation of drops on slippery lubricant-infused surfaces Soft Matter 11:7617–7626.
- (23) Rykaczewski K, Paxson AT, Staymates M, Walker ML, Sun X, et al. (2014) Dropwise condensation of low surface tension fluids on omniphobic surfaces Sci. Rep. 4:4158.
- (24) Kim P, Wong TS, Alvarenga J, Kreder MJ, Adorno-Martinez WE, et al. (2012) Liquid-infused nanostructured surfaces with extreme anti-ice and anti-frost performance ACS Nano 6(8):6569–6577.
- (25) Subramanyam SB, Rykaczewski K, Varanasi KK (2012) Ice adhesion on lubricant-impregnated textured surfaces Langmuir 29(44):13414–13418.
- (26) Discher D, Janmey P, Wang Y (2005) Tissue cells feel and respond to the stiffness of their substrate Science 310:1139–1143.
- (27) Lo CM, Wang HB, Dembo M, Wang YL (2000) Cell movement is guided by the rigidity of the substrate Biophys. J. 79(1):144–152.
- (28) Zhang R, Ni L, Jin Z, Li J, Jin F (2014) Bacteria slingshot more on soft surfaces Nature Commun. 5:5541.
- (29) Guo WH, Frey MT, Burnham NA, Wang YL (2006) Substrate rigidity regulates the formation and maintenance of tissues Biophys. J. 90(6):2213–2220.
- (30) Winter H Chambon F (1986) Analysis of linear viscoelasticity of a cross-linking polymer at the gel point J. Rheology 30(2):367–382.
- (31) Long D, Ajdari A, Leibler L (1996) Static and dynamic wetting properties of thin rubber films Langmuir 12(21):5221–5230.
- (32) de Gennes PG (1996) Soft adhesives Langmuir 12:4497–4500.
- (33) Carre A, Gastel J, Shanahan M (1996) Viscoelastic effects in the spreading of liquids Nature 379(6564):432–434.
- (34) Huang J, Grater SV, Corbellini F, Rinck S, Bock E, Kemkemer R, Kessler H, Ding J, Spatz JP (2009) Impact of order and disorder in rgd nanopatterns on cell adhesion Nano Letters 9(3):1111–1116.
- (35) Dalby MJ, Gadegaard N, Oreffo RO (2014) Harnessing nanotopography and integrin-matrix interactions to influence stem cell fate Nature materials 13(6):558–569.
- (36) Krieg M, Arboleda-Estudillo Y, Puech PH, Käfer J, Graner F, Müller D, Heisenberg CP (2008) Tensile forces govern germ-layer organization in zebrafish Nature cell biology 10(4):429–436.
- (37) Tinevez JY, Schulze U, Salbreux G, Roensch J, Joanny JF, Paluch E (2009) Role of cortical tension in bleb growth Proc. Natl. Acad. Sci. USA 106(44):18581–18586.
- (38) Fischer-Friedrich E, Hyman AA, Jülicher F, Müller DJ, Helenius J (2014) Quantification of surface tension and internal pressure generated by single mitotic cells Scientific reports 4.
- (39) Sliogeryte K, Thorpe SD, Lee DA, Botto L, Knight MM (2014) Stem cell differentiation increases membrane-actin adhesion regulating cell blebability, migration and mechanics Scientific reports 4.
- (40) Swift J, Ivanovska IL, Buxboim A, Harada T, Dingal PDP, Pinter J, Pajerowski JD, Spinler KR, Shin JW, Tewari M, et al. (2013) Nuclear lamin-a scales with tissue stiffness and enhances matrix-directed differentiation Science 341(6149):1240104.