Droplets move over viscoelastic substrates by surfing a ridge

Abstract

Liquid drops on soft solids generate strong deformations below the contact line, resulting from a balance of capillary and elastic forces. The movement of these drops may cause strong, potentially singular dissipation in the soft solid. Here we show that a drop on a soft substrate moves by surfing a ridge: the initially flat solid surface is deformed into a sharp ridge whose orientation angle depends on the contact line velocity. We measure this angle for water on a silicone gel and develop a theory based on the substrate rheology. We quantitatively recover the dynamic contact angle and provide a mechanism for stick-slip motion when a drop is forced strongly: the contact line depins and slides down the wetting ridge, forming a new one after a transient. We anticipate that our theory will have implications in problems such as self-organization of cell tissues or the design of capillarity-based microrheometers.

ARTICLE

Received 20 Feb 2015 | Accepted 23 Jun 2015 | Published 4 Aug 2015

Liquid drops on soft solids generate strong deformations below the contact line, resulting from a balance of capillary and elastic forces. The movement of these drops may cause strong, potentially singular dissipation in the soft solid. Here we show that a drop on a soft substrate moves by surng a ridge: the initially at solid surface is deformed into a sharp ridge whose orientation angle depends on the contact line velocity. We measure this angle for water on a silicone gel and develop a theory based on the substrate rheology. We quantitatively recover the dynamic contact angle and provide a mechanism for stickslip motion when a drop is forced strongly: the contact line depins and slides down the wetting ridge, forming a new one after a transient. We anticipate that our theory will have implications in problems such as self-organization of cell tissues or the design of capillarity-based microrheometers.

DOI: 10.1038/ncomms8891 OPEN

Droplets move over viscoelastic substrates by surng a ridge

S. Karpitschka1, S. Das2, M. van Gorcum1, H. Perrin3, B. Andreotti3 & J.H. Snoeijer1,4

1 Physics of Fluids Group, Faculty of Science and Technology, Mesa Institute, University of Twente, 7500 AE Enschede, The Netherlands. 2 Department of

Mechanical Engineering, University of Maryland, College Park, Maryland 20742, USA. 3 Physique et Mcanique des Milieux Htrognes, UMR 7636 ESPCI-CNRS, Universit Paris-Diderot, 10 rue Vauquelin, 75005 Paris, France. 4 Department of Applied Physics, Eindhoven University of Technology, PO Box 513, 5600 MB Eindhoven, The Netherlands. Correspondence and requests for materials should be addressed to SK(email: mailto:[email protected]

Web End [email protected] ) or to B.A. (mailto:[email protected]

Web End [email protected]) or to J.H.S. (mailto:[email protected]

Web End [email protected]).

NATURE COMMUNICATIONS | 6:7891 | DOI: 10.1038/ncomms8891 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

& 2015 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8891

Capillary interactions of soft materials are ubiquitous in nature and play a major role in self-organization of cell tissues1, in embryotic development2,3, in wound healing4,

and in controlling the spreading of cancer cells5,6. Not least motivated by this, soft wetting79 recently came to the attention of both physicists and biologists. Despite its potential for applications, such as the patterning of cells10 or droplets11 onto soft surfaces, or the optimization of condensation processes12, our fundamental understanding of the dynamics of soft wetting lags behind by far of what is known about rigid surfaces13.

Partial wetting of a liquid on a rigid (smooth) substrate is controlled by intermolecular interactions, whose strength is characterized by surface energies13. The motion of the three-phase contact line is governed by the viscous dissipation in the liquid. A dissipation singularity arises at the moving contact line14 and its regularization at the nanoscopic scale can result from various processes13,15. When the substrate is deformable, a sharp ridge forms below the contact line at the edge of the droplet79. The ridge geometry (Fig. 1a) originates from the coupling between elasticity and surface energy1620. The problem is inherently multi-scale and non-local, even at equilibrium, due to the long range of elastic interactions2123.

Pioneering experiments have shown that the softness drastically slows down the wetting dynamics24,25 in comparison to rigid solids. This viscous braking in the wetting of soft solids has been attributed to a viscoelastic force, as discussed in several recent experimental articles11,26. The theoretical description of moving contact lines over soft solids is so far limited to global dissipation arguments17, which, at least for wetting of rigid solids, are known not to capture the entire physics behind the process.

In this paper, the physical mechanism that governs soft wetting dynamics is revealed. We measure the dynamical wetting of small water droplets on a rheologically characterized silicone gel and discover a saturation of the dynamical contact angle for large speeds, associated with a maximum friction force. Driving the contact line motion beyond this maximal force eventually leads to a dynamical depinning where the contact line surfs down the wetting ridge, providing a mechanism for recently observed stickslip motion26. We develop a theoretical framework for dynamical soft wetting, suitable for any substrate rheology. The dynamic wetting angle is calculated from the velocity-dependent shape of the deformed solid (Fig. 1c). The experimental results

are matched quantitatively, including saturation/dynamical depinning. The latter arises from an upper limit of the viscoelastic braking effect, which, by exploring different rheologies, is found to be a robust, universal feature of soft wetting and should thus be relevant far beyond droplets on silicone gels. In addition, the analysis captures recent X-ray measurements on the slow growth of wetting ridges when a drop is deposited on a substrate27 (Fig. 1b).

ResultsExperiments. Experiments were performed using water drops on a silicone gel (cf. Methods section for details). This system was previously used in static28 and transient27 experiments. Figure 2a shows the rheology of this gel; similar data were reported in ref. 29. The storage G0 and loss G00 moduli are related by

KramersKronig relation: they originate from the same relaxation function C(t). More precisely, the complex shear modulus obey the relation m o

G0 iG00 io R

1

0 dt C t

exp iot.

A silicone gel is a reticulated polymer formed by polymerizing small multifunctional prepolymers: contrarily to other types of gels, there is no liquid phase trapped inside. Such cross-linked polymer networks exhibit scale invariance that yields power-law response of the form 17,30,31:

C t

G 1 G 1 n

1

t t

n

h i

; 1 where G is a static shear modulus and G is the gamma function. The associated complex modulus reads m G0 iG00 G[1 (ito)n].

The value of n is not universal but depends on the stoichiometric ratio r between reticulant and prepolymer, with n typically between 2/3 and 1/2 (ref. 32). The best t in Fig. 2a gives an

a

104

G, G(Pa)

103

102

1/2

1

1

10 102 101 1 101 102

(rad s1)

a b

b

102

Stationary Deposition

Moving Depinning

t

101

Time

max crit

1/2

1

s

()

h0

h(x,t)

1

c

d

10 106 105 104 103 102 101 1

1

v

v (mm s1)

Figure 1 | Dynamics of wetting ridges. (a) Equilibrium deformation by a three-phase contact line, inducing a solid contact angle ys. The liquid contact angle is denoted y. (b) Growth of the wetting ridge after a drop is deposited. (c) Contact line moving at a velocity V. The motion induces a rotation j of the wetting ridge and the liquid contact angle, while ys remains constant. (d) Dynamical depinning occurs at a critical angle jcrit.

Figure 2 | Rheology of the substrate and dynamic contact angle. (a) Storage modulus G0(o) (open symbols) and loss modulus G00(o)

(closed symbols) of the silicone gel. Lines are best t of m G(1 (ito)n),

giving n 0.55, G 1.2 kPa and t 0.13 s. (b) Dynamic angle j y yeq

for water on the silicone gel (symbols). Data are averaged over 10 independent experiments, error bars represent the s.d. The small-v behaviour exhibits the same power-law as G00. Dashed line is the bestt of the asymptote equation (8). Solid line corresponds to equation (7), describing the full range of velocities. The critical angle of depinningjcrit 39 3 , measured separately, is plotted at an arbitrary velocity.

max

2 NATURE COMMUNICATIONS | 6:7891 | DOI: 10.1038/ncomms8891 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2015 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8891 ARTICLE

r ^

s 0, the constitutive relation (equation (3)) and the traction

at the free surface is identical to the static problem, but features dependences on the frequency o. The time-dependent traction can therefore be solved analogously to (refs 9,22,23), by an additional spatial Fourier transformation (noted by B):

^~

h q; o

^~G q; o

exponent n 0.55. Note that the associated effective viscosity

G00/oBG(t/o)n is very large, beyond 10 Pa.s, over the entire frequency domain. This will imply that the dissipation mainly occurs in the solid, not in the liquid.

Figure 2b shows the dynamical angle j as a function of velocity, both of which are measured while the droplet slowly relaxes over time towards its equilibrium shape (after the injection phase). The resulting j versus v is not sensitive to the history of the relaxation process, and the dynamics can thus be considered quasi-steady. The loglog plot reveals a power-law relation between j and v at small velocities, with an exponent equal to the rheological value of n 0.55, within error bars. Such

a power-law dependence is similar to previously obtained results24,25. For large velocities, we report a striking saturation of the dynamical contact angle. Neither the small-v power law nor the saturation can be explained by dissipation mechanisms in the liquid, and one needs to account for the dissipation within the solid. Long et al.17 addressed this using a global dissipation argument based on the solid rheology, but this argument fails to capture key features such as the saturation.

When the drop is kept inating with a large over-pressure, we observe a depinning of the wetting front, with a sudden increase of its velocity, as the dynamical angle reaches a value jcrit 39 3 . This compares well to the saturation of 37

observed during relaxation (after injection), as indicated in Fig. 2. When forcing the contact angle beyond this angle, the contact line dynamically depins from the wetting ridge, surng it, until a new ridge forms. Note that such a depinning, leading to stickslip motions, had already been reported in refs 26,33. Our current measurements show that this is a direct consequence of the saturation of the dynamic contact angle.

Theoretical framework. A liquid drop deposited on a substrate exerts a capillary traction on the surface7,3437. While the resulting elastic deformation has been computed and measured for static situations8,9,18,20,23,28, the traction becomes time dependent in the case of dynamical wetting. Here we consider a single straight contact line, for which the elastic problem is two-dimensional. Our goal is to nd the deformation of the solid, h(x, t), resulting from the time-dependent capillary traction, T(x, t). For simplicity, we consider the same surface energy gs for the wet and the dry parts of the substrate, and assume that there is no Shuttleworth effect, i.e. that gs does not depend on strain16,19. The resulting traction on the solid is then purely normal, and reads T(x, t) gs qxxh, the latter term being the solid

Laplace pressure9,23,38,39. The theory is rigorous for small slopes (qxh)2oo1, but can be extended in a semi-quantitative way to nite slopes.

The shape of the deformed substrate h(x, t) follows from the normal substrate displacements. Inside a purely elastic material, the displacements adapt instantaneously to changes in the capillary traction; the problem is therefore essentially static. For realistic soft materials, however, the displacements are delayed with respect to the imposed forcing. For small deformations, this is captured by a linear stressstrain rate relation

sij x; t

Z

t

1

^~

T q; o

1 gsq2

^~G q; o

; 4

where q is the wavenumber. The Greens function ^~

Gq; o is the

product of the time kernel m(o) 1 by the space kernel K(q). For an incompressible layer of thickness h0 (ref. 40) it is

K q

sinh 2qh0

2qh0

cosh 2qh0

1

2q 5

Leftright symmetry and volume conservation are reected by K(q) K( q) and K(0) 0. Sharp features in the solid prole,

like the solid contact angle, are found in the large-q asymptotics for which K(q)C(2|q|) 1.

The moving contact line. We now apply our theory to a contact line moving at a constant velocity v, which induces a traction

T x; t

g sinyd x vt

6 This reects the normal force per unit contact line that is exerted by the liquid on the solid, while y is the liquid angle at the location of the cusp. For simplicity we consider that the drop size is much larger than the substrate thickness, in which case the Laplace pressure inside the liquid can be neglected28. We indeed veried that the nite drop size has a negligible inuence on the resulting motion: the relevant scale for the dynamics is the size of the ridge gs/G, which is much smaller than the drop size.

This also justies a two-dimensional model. Another important simplication comes from the quasi-steady nature of the droplet relaxation: temporal changes of contact angle and contact line velocity are small in our experiments (dy/dtoot 1), so that the process can be modelled by a constant velocity.

According to equation (4), the capillary traction induces a wetting ridge moving at a velocity v (see Methods section):

q

g sinygs q2

2 qh0

2 1

m qv

=gs

Kq

1 ; 7

in the co-moving frame. In real space this gives proles such as shown in Fig. 1c. The motion induces a leftright symmetry breaking: the asymmetric deformation of the solid results in a tilt angle jv of the cusp. Since the liquid is close to equilibrium, the

change in the solid angle j induced by the dissipation taking place inside the solid directly yields a change in the liquid angle y (Fig. 1d). Hence y yeq j, where yeq is the equilibrium liquid

angle by Neumanns law. The liquid contact angle y gets deviated away from yeq due to the viscoelastic forces in the substrate.

Figure 2b shows the calculated tilt curve for the gel used in experiments. It predicts not only the power-law behaviour at low velocity but also presents a saturation of the tilt angle. The tilt angle quanties the velocity-dependent viscoelastic force between the solid and liquid phases. For a well-established moving ridge, it behaves as a resistive force increasing with the velocity. When the drop is forced to inate with a driving force larger than the maximal braking force, the contact line can no longer remain pinned to the steadily moving solid ridge and surfs the gel wave. To investigate further the relation between the tilt j and the substrate constitutive relation, we use the gel rheology (equation (1)), and expand (equation (7)) in the small-v asymptotics (and hence small j, that is, sinyEsinyeq), which

dij; 2

where C is the relaxation function previously introduced (see equation (1)) and p is the pressure. Like ref. 17, we apply a Fourier transform with respect to time (noted by ^):

sij x; o

^ m o

^Eij x; o

^p x; o

dij; 3

The mathematical problem dened by mechanical equilibrium,

dt0 C t t0

@t

0 Eij x; t0

p x; t

NATURE COMMUNICATIONS | 6:7891 | DOI: 10.1038/ncomms8891 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

& 2015 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8891

gives (see Methods section):

j

2n 1n

cos np=2

g siny gs

v v

n

; 8

where the characteristic velocity scale emerges as v* gs/(Gt).

Note that the outer length scale (thickness of substrate) does not appear. This expression can be simply interpreted. At vanishing response time t, a deformation matching the static ridge would propagate at a velocity v, pushing the substrate material up and down at a characteristic frequency o equal to the velocity v divided by the characteristic width of the ridge Bgs/G. The

perturbation introduced by a nite t is encoded by the loss modulus G00(o). As the characteristic strain is set by the slope of the ridge Bgsiny/gs, one obtains dimensionally equation (8). The scaling law j / v = v n thus simply carries over the low

frequency behaviour G00(o)pon, which is a robust mechanism valid beyond the small-slope approximation of our theory. At small v (small o), dissipation will dominantly occur in the solid because no1, while the loss modulus of a Newtonian liquid

G00liq / o.

In contrast to the tilt angle j, we nd that that the solid angle ys does not depend on v. This result can be derived analytically from the large-q asymptotics of equation (7), gsinygs q 2, valid for all v and arbitrary m(o). In real space, this implies a slope discontinuity

ys

p g siny=gs; 9

which is the small-slope limit (g/gsoo1) of Neumanns contact angle law. Physically, equation (9) reects that ys is determined by surface tensions only9,18,41: bulk viscoelastic stresses are not singular enough to contribute to the contact angle selection. This feature remains true for arbitrary angles23. To match Neumanns law quantitatively, our theory must be corrected at large slopes to take geometric nonlinearities into account. This can be achieved phenomenologically in equation (4) by replacing

gs !

q g2s g2 4.

Indeed, the Neumann condition for y 90

reads 2gs sin a g, where a is the angle of the solid interface with

the horizontal. Small-slope theory gives 2 h0

j j 2tana ggs, and

hence lacks a factor cos a

1 g=2gs

2

q .

For the rst time, we reveal that the exponent of the dynamical contact angle directly originates from the gel rheology. The rheological parameters being calibrated independently, the dynamic contact angle can be tted to the model to extract the solid surface tension. Using equation (8), which is valid for small slopes, we nd gs 16 mN m 1. This is a reasonable value,

though a bit lower than the value previously derived from Neumanns law28. We think this difference can be attributed to the small-slope nature of our theory: condering the phenomenological correction for large slopes gives a value gs

39 mN m 1, in close agreement with ref. 28. The solid line in Fig. 2b shows the prediction from equation (7), providing an excellent description over the full range of velocities. The model captures also the saturation, though the value for jcrit is slightly overestimated.

Depinning and growth of a new wetting ridge. How can the contact line escape pinning, without dragging the capillary wedge along with it? To answer this question, let us consider the recent experiments investigating the growth of a wetting ridge after depositing a droplet on a silicone gel27. The substrate was observed to only very slowly establish the nal shape of the wetting ridgesuch a delay in growth (or decay) of wetting ridges would explain how a sufciently rapid contact line could escape

from the ridge. However, the solid angle ys (cf. Fig. 1a) appeared very quickly and remained constant during the entire growth of the ridge27.

These features of ridge growth can all be explained by considering our theory for a traction that is suddenly imposed at the time of deposition (t 0), so that

T x; t

g siny d x

Y t

10 where y is the liquid contact angle and Y(t) is the Heaviside step function. Combining this traction with equation (4), one can compute the resulting h(x, t) for any rheology m(o). An example of the evolution of the wetting ridge is shown in Fig. 1b. A movie is given in Supplementary Movie 1.

First, the theory recovers the experimental nding that ys is

constant at all times. As for the moving contact line, this result can be derived analytically from the large-q asymptotics of equation (4), which again results in equation (9): the asymptotics are independent of the rheology and the history of the traction, but entirely governed by the surface tensions.

Second, the theory explains why, contrarily to the rapid appearance of ys, the global shape of the ridge evolves much more slowly. Figure 3 shows the evolution of the central height of the ridge, h(x 0), towards its static value hN, for the two idealized

rheological models. The relaxation towards the equilibrium height is algebraic for the gel model, with an exponent directly following that of rheological relaxations (Fig. 3a, as t 1/2 for n 1/2). This claries the complex evolution of the wetting ridge

of the silicone gel in ref. 27: small-scale characteristics such as ys are dominated by surface tension and relax quickly, while large-scale features inherit the relaxation dynamics from the bulk rheology. This means that immediately after depinning, where the contact line exhibits a rapid motion, the solid cusp cannot adapt quickly. The liquid will slide down the wetting ridge, which appears frozen on the timescale of the depinning. During this phase it is clear that the liquid dynamics will be importantstill, the onset of the depinning can be explained quantitatively without invoking the uid dynamics inside the liquid, because the saturation angle jmax coincides with the observed depinning angle jcrit (Fig. 2b).

Robustness and interpretation. The theory can be applied to any viscoelastic substrate, assuming that it is probed in the linear regime. Generic reticulated polymer networks possess a long-time entropic elasticity42,43, that is characterized by a static shear modulus G. Such networks become viscoelastic when excited over timescales shorter than a certain response time t. To investigate the robustness of the phenomenology that was observed experimentally and reproduced quantitatively by our model, we

a

1

b

(h h) /(/G)

101

102

1/2

1

10 105 103 101 101 103

-3

0 1 2 3 4 5 t / ( )

t /

Figure 3 | Relaxation of the central height of the wetting ridge after drop deposition. The curves show the approach to equilibrium height, hN h at x 0, for two rheologies: (a) power-law rheology (equation (1) with

n 1/2) and (b) standard linear model, b 300 (equation (11), see inset).

The dimensionless substrate thickness for these plots was gs/(Gh0) 0.5.

4 NATURE COMMUNICATIONS | 6:7891 | DOI: 10.1038/ncomms8891 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2015 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8891 ARTICLE

101

0.5

tan( ) / ( sin( )/ s)

1

= 1,000 100

10

101

xh / ( / s)

1.0

102

10 104 102 1 102

3

v / ( /(G )) = 0.1

v / ( / (G ))

Figure 4 | Cusp tilt for standard linear model. Solid lines: numerical results; dashed lines: analytical approximation (equation (14)). The upper bound for the viscous braking force is robust with respect to the details of the rheology.

= 102

= 104

1.5 106 105 104 103 102 101

x / ( /G)

will consider a different rheological limit. When cross-linking long polymer chains, one forms an elastomer rather than a gel. Assuming a single (Rouse) timescale t to characterize the onset of entanglements, the rheology can be idealized as17,31,43:

C t

G 1 be t =t

Figure 5 | Logarithmic variation of the prole slope. qxh is plotted as a function of the distance to the contact line. Different curves correspond to different instantaneous relaxation moduli bE, for identical dimensionless velocity v 1. In the limit b-N (the KelvinVoigt solid), the slope diverges

logarithmically at the contact line. For nite b, the slope saturates at a regularization length given by equation (13).

: 11 This single timescale response is also referred to as standard linear model and has frequently been used to describe the transition from rubber to glass behaviours. In general, several relaxation times must be introduced to capture quantitatively the rheology of actual elastomers.

The wetting ridge relaxation, which follows the rheological relaxation, becomes exponential for the standard linear model (Fig. 3b). The case of a moving contact line with the rheology (equation (11)) is given by the solid lines in Fig. 4, showing the tilt curves for various parameters b. As for the gel case, the tilt has an upper bound and this therefore appears a robust feature of soft wetting. Note that the maximum depends on g/gs and b (or n for the gel case), illustrating that the value of jmax will in general depend on the details of the rheology (for example, jmax in ref. 26 is much smaller than that in the presented experiments). The dynamic contact angle for a standard linear solid can actually be captured in a simple analytical form. For this we consider the limit b-N while keeping t0 bt constantthis corresponds to

the KelvinVoigt model with a frequency-independent effective viscosity Z Gt0. Intriguingly, this limit turns out to be singular:

the high-frequency behaviour of equation (11) becomes purely viscous and gives a non-integrable singularity of the dissipation. This singularity could already be anticipated from equation (8), since the KelvinVoigt rheology has G00Bo1, while equation (8)

presents a divergence for n 1.

In fact, this viscoelastic singularity is the soft-solid analogue of the classical Huh and Scriven paradox for viscous contact line motion14. This is demonstrated from the large-q asymptotics of equation (7) in the KelvinVoigt limit, giving a slope close to the contact line (see the Methods section):

@xh

g siny 2pgs

: 14

This result is analogous to the CoxVoinov law44,45 in partial wetting of viscous uids. In that case, a similar logarithmic factor linking microscopic and macroscopic scales appears, for arbitrary contact angles45, and the resulting expression for small Ca is of the form equation (14). Interestingly, the analysis reveals that the relevant dimensionless velocity for soft wetting is not the classical liquid capillary number Ca vZ = g, based on the liquid viscosity

Z, but the solid capillary number

Cas

vGbtg g2s

h0 egE

; 12

with gE Eulers constant and v* gs/(Gt0). This expression

reveals a logarithmic divergence of the slope, in perfect analogy to the CoxVoinov result for liquid contact lines44,45.

Contrarily to the viscous-liquid singularity, the presence of an instantaneous elastic response, that is, a nite value of b, is sufcient to regularize the divergence. This is illustrated in Fig. 5, which shows the slope ahead of the moving contact line as a function of the distance to the contact line. The dashed line

corresponds to the KelvinVoigt limit equation (12), showing the logarithmic steepening of the slope. For nite b, the slope saturates on approaching the moving contact line. The saturation wavenumber is found q0B(vt) 1, which corresponds to a length

tv: 13 where is a dynamical regularization length that depends linearly on the velocity of the contact line; for v O1 this scale is still

much smaller than the substrate thickness, by a factor b 1. The physical origin of the regularization lies in the instantaneous elasticity in the high-frequency limit, which applies at frequencies beyond Bbt.

Inserting the regularization length into the KelvinVoigt limit equation (12), we identify the tilt and get an analytical expression of the dynamic liquid contact angle (the strict validity of the analysis requires small slopes, that is, small j; we therefore replaced siny sinp = 2 j 1):

y yeq

2p Cas ln

: 15

Equation (14) closely follows the numerical results (Fig. 5, dashed and solid lines, respectively). The moving contact line singularity is avoided altogether when G00 has an exponent no1, as was the case for the power-law gel, in perfect analogy to shear-thinning uids moving on a rigid substrate.

DiscussionWe have shown how contact lines can surf on a wetting ridge, and that this governs the remarkable spreading of drops on viscoelastic substrates. We have quantied this dynamics by

v=v

1 4 v=v

2

x

j j=x 4 gE ln x

j j=h0

NATURE COMMUNICATIONS | 6:7891 | DOI: 10.1038/ncomms8891 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

& 2015 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8891

measuring the dynamic contact angle of water on a silicone gel for a wide range of velocities and described, for the rst time, a saturation of the dynamical angle for large velocities. This saturation is in harsh contrast to wetting dynamics of rigid solids and leads to depinning, where the contact line slides down the ridge until a new wetting ridge has had time to grow and sustain a steady motionhence, explaining the remarkable stickslip motion26 found recently on soft solids. We develop a theory that identies a robust maximum in viscous braking force that correctly predicts the onset of dynamical depinning. In addition, our theory captures the unsteady growth of a wetting ridge27. This work provides a framework for viscoelasto-capillary dynamics valid beyond droplets, and should be applicable, for example, within a biological context. It also opens a new perspective where droplets can be used as microrheometers, since the length scales probed by the droplets are given by the elastocapillary length (that is, a few microns), and the tilt saturation occurs at velocities that are directly related to the relaxation timescale.

Methods

Wetting experiments. The silicone gels (Dow Corning CY52-276) are prepared by curing the mixed components onto glass slides, yielding 0.8-mm thick substrates. The rheology was determined using a MCR 501 rheometer (Anton Paar). Dynamic contact angles were measured using droplets of MilliQ water dispensed from a clean Hamilton syringe. First, a small droplet (220 ml)

was placed onto the substrate, leaving the syringe needle attached to the droplet. Then, the contact angle of the droplet was increased by quickly injecting water (3 20 ml with 28 ml s 1) into it. After the injection phase the drop relaxes quasi-statically, causing the contact line velocity to decay slowly. The advancing motion of the contact line and the relaxation of the contact angle were imaged at 50 Hz with a long-distance video microscope. The droplet contour was extracted with sub-pixel resolution, and velocities down to Bnm s 1 could be detected.

The measured contact angles were translated to tilt angles j by subtracting yeqE1061.

The moving contact line. Fourier transforming equation (6) from x to q and from t to o preserves the d-shape of the traction:

^~

T q; o

2p g sin y d o vq

: 16 Inserting the above into equation (4), the inverse transform to the timedomain yields

~h q; t

g sin ygs q2

of rheology:

lim

x!0

h0x lim

x!0

h0x

21

The rotation of the wetting ridge is given by the symmetric part of h0(x) and

thus obtained by the backward transform of the real part, evaluated at x 0:

tanj lim

x!0

1

2 h0 x

h0 x

1p lim

x!0

0

@

Z1

1

ig sinygsq e iqxdq

1

A

g siny gs :

1

2p

Z1

1

22

< iq

~h q

h i

dq

With the symmetry property K(q) K( q) and with small, positive v,

equation (22) simplies to (primes omitted):

j

g siny G

sinnp=2

p

Z1

0

qqvtnKq

gs = Gq

dq: 23

In the limit of thick elastic layers, K(q) (2|q|) 1. After non-dimensionalizing the

integration variable as q0 gG q, one obtains (primes omitted):

j

g siny gs

2Kq 1

2

nsinnp=22

p

Z1

0

v gs= Gt

qn q=2 1

2

dq

24

2n 1n

cosnp = 2

g siny gs

v v

n

;

where v* gs/(Gt) is the characteristic velocity.

Growth of a wetting ridge. Here we give the full derivation of the time-dependent wetting ridge shape after the deposition of a droplet. We only discuss the result for the exponential relaxation model. An analogous calculation can be performed for the power-law relaxation.

In the following, we non-dimensionalize x with h0, q with h 10, t with bt, o with (bt) 1 and h with g siny/G. With this scaling, the Fourier transform of the time kernel for exponential relaxation (equation (11)) reads

m o

G 1

o o=b i

: 25

The space kernel in scaled variables is

k q

K q

h0

1

2q : 26

sinh 2q

2q

cosh 2q

2q

2 1

The traction equation (10) is transformed to

^~

T q; o

g siny

h0

1eiqvt: 17

The only explicit time dependence appears in the phase factor that shifts the prole in x-direction linearly with time. The transformation to the co-moving frameis done by multiplication with e iqvt, which cancels the only explicit time dependence, and one obtains equation (7).

The slopes are evaluated by multiplication with iq before the inverse

transform (in the co-moving frame):

h0 x

1

2p

m vq

K q

: 27

equations (16)(18) are inserted into the general expression equation (3), which yields:

^~h q; o

i o

io pd o

o o=b i

k q

1

pd o

2 ; 28

with the dimensionless parameter as gs/(Gh0). as compares the elastocapillary

length for the solid surface tension to the layer thickness h0. The inverse Fourier transform to the time domain yields

q; t

1

1 asq

Z1

iq~h q

e iqxdq: 18

h0(x) is a real function because < iq

1 ~h q <iq

q

and I iq

~h q

b exp bt

h i

1 b a q k q

k q

1 asq

2 : 29

Fourier transformation to real space is performed numerically.

References

1. Manning, M. L., Foty, R. A., Steinberg, M. S. & Schoetz, E.-M. Coaction of intercellular adhesion and cortical tension species tissue surface tension. Proc. Natl Acad. Sci. USA 107, 1251712522 (2010).

2. Trinkhaus, J.-P. & Groves, P. Differentiation in cultures of mixed aggregates of dissociated tissue cells. Proc. Natl Acad. Sci. USA 41, 787795 (1955).

3. Steinberg, M. Reconstruction of tissues by dissociated cells. Science 141,

401408 (1963).

4. Armstrong, M. & Armstrong, P. Mechanisms of epibolic tissue spreading analyzed in a model morphogenetic systemroles for cell-migration and tissue contractility. J. Cell. Sci. 102, 373385 (1992).

5. Douezan, S., Dumond, J. & Brochard-Wyart, F. Wetting transitions of cellular aggregates induced by substrate rigidity. Soft Matter 8, 45784583 (2012).

Iiq

q

. h0(x) can be split into a symmetric and an antisymmetric part,

where the symmetric part is given by the inverse transform of the real part of iq~hq:

1

2 h0 x

h0 x

1

2p

Z1

1

< iq

q

h i

e iqxdq; 19

The antisymmetric part is obtained form the imaginary part:

1

2 h0 x

h0 x

1

2p

Z1

1

I iq

q

h i

e iqxdq: 20

The solid angle ys is given by the (antisymmetric) slope discontinuity at x 0 and is thus encoded in the backward transform of the imaginary

part. The discontinuity is caused by the large-q asymptotics alone. If O m vq

o O q

, which is the case for the exponential- and power-law (no1)

relaxation (but not for the KelvinVoigt model), it is independent

6 NATURE COMMUNICATIONS | 6:7891 | DOI: 10.1038/ncomms8891 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2015 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms8891 ARTICLE

6. Sabari, J. et al. Fibronectin matrix assembly suppresses dispersal of glioblastoma cells. PLoS ONE 6, e24810 (2011).

7. Shanahan, M. The inuence of solid micro-deformation on contact-angle equilibrium. J. Phys. D: Appl. Phys. 20, 945950 (1987).

8. Pericet-Camara, R., Best, A., Butt, H.-J. & Bonaccurso, E. Effect of capillary pressure and surface tension on the deformation of elastic surfaces by sessile liquid microdrops: an experimental investigation. Langmuir 24, 1056510568 (2008).

9. Jerison, E. R., Xu, Y., Wilen, L. A. & Dufresne, E. R. Deformation of an elastic substrate by a three-phase contact line. Phys. Rev. Lett. 106, 186103 (2011).

10. Discher, D., Janmey, P. & Wang, Y. Tissue cells feel and respond to the stiffness of their substrate. Science 310, 11391143 (2005).

11. Style, R. W. et al. Patterning droplets with durotaxis. Proc. Natl Acad. Sci. USA 110, 1254112544 (2013).

12. Sokuler, M. et al. The softer the better: fast condensation on soft surfaces. Langmuir 26, 15441547 (2010).

13. Bonn, D., Eggers, J., Indekeu, J., Meunier, J. & Rolley, E. Wetting and spreading. Rev. Mod. Phys. 81, 739805 (2009).

14. Huh, C. & Scriven, L. E. Hydrodynamic model of steady movement of a solid/ liquid/uid contact line. J. Coll. Int. Sci. 35, 85101 (1971).

15. Snoeijer, J. H. & Andreotti, B. Moving contact lines: scales, regimes, and dynamical transitions. Ann. Rev. Fluid Mech. 45, 269292 (2013).

16. Shuttleworth, R. The surface tension of solids. Proc. Phys. Soc. A 63, 444457 (1950).

17. Long, D., Ajdari, A. & Leibler, L. Static and dynamic wetting properties of thin rubber lms. Langmuir 12, 52215230 (1996).

18. Limat, L. Straight contact lines on a soft, incompressible solid. Eur. Phys. J. E 35, 134 (2012).

19. Weijs, J. H., Andreotti, B. & Snoeijer, J. Elasto-capillarity at the nanoscale: on the coupling between elasticity and surface energy in soft solids. Soft Matter 9, 84948503 (2013).

20. Bostwick, J., Shearer, M. & Daniels, K. Elastocapillary deformations on partially-wetting substrates: rival contact-line models. Soft Matter 10, 73617369 (2014).

21. White, L. The contact angle on an elastic substrate. 1. The role of disjoining pressure in the surface mechanics. J. Colloid Interface Sci. 258, 8296 (2003).

22. Style, R. W. & Dufresne, E. R. Static wetting on deformable substrates, from liquids to soft solids. Soft Matter 8, 71777184 (2012).

23. Lubbers, L. A. et al. Drops on soft solids: free energy and double transition of contact angles. J. Fluid. Mech. 747, R1 (2014).

24. Shanahan, M. & Carre, A. Viscoelastic dissipation in wetting and adhesion phenomena. Langmuir 11, 13961402 (1995).

25. Carre, A., Gastel, J. & Shanahan, M. Viscoelastic effects in the spreading of liquids. Nature 379, 432434 (1996).

26. Kajiya, T. et al. Advancing liquid contact line on visco-elastic gel substrates: stick-slip vs. continuous motions. Soft Matter 9, 454461 (2013).

27. Park, S. et al. Visualization of asymmetric wetting ridges on soft solids with x-ray microscopy. Nat. Commun. 5, 4369 (2014).

28. Style, R. W. et al. Universal deformation of soft substrates near a contact line and the direct measurement of solid surface stresses. Phys. Rev. Lett. 110, 066103 (2013).

29. Style, R. W. et al. Traction force microscopy in physics and biology. Soft Matter 10, 40474055 (2014).

30. Winter, H. & Chambon, F. Analysis of linear viscoelasticity of a cross-linking polymer at the gel point. J. Rheol. 30, 367382 (1986).

31. de Gennes, P.-G. Soft adhesives. Langmuir 12, 44974500 (1996).32. Scanlan, J. C. & Winter, H. H. Composition dependence of the viscoelasticity of end-linked poly(dimethylsiloxane) at the gel point. Macromolecules 24, 4754 (1991).

33. Kajiya, T. et al. A liquid contact line receding on a soft gel surface: dip-coating geometry investigation. Soft Matter 10, 88888895 (2014).

34. Lester, G. Contact angles of liquids at deformable solid surfaces. J. Colloid Sci. 16, 315326 (1961).

35. Rusanov, A. Theory of wetting of elastically deformed bodies.1. Deformation with a nite contact-angle. Colloid J. USSR 37, 614622 (1975).

36. Rusanov, A. Thermodynamics of deformable solid-surfaces. J. Colloid Interface Sci. 63, 330345 (1978).

37. de Gennes, P.-G., Brochard-Wyart, F. & Quere, D. Capillarity and Wetting Phenomena: Drops, Bubbles, Pearls, Waves (Springer, 2004).

38. Mora, S., Phou, T., Fromental, J.-M., Pismen, L. M. & Pomeau, Y. Capillarity driven instability of a soft solid. Phys. Rev. Lett. 105, 214301 (2010).

39. Marchand, A., Das, S., Snoeijer, J. H. & Andreotti, B. Capillary pressure and contact line force on a soft solid. Phys. Rev. Lett. 108, 094301 (2012).

40. Xu, Y. et al. Imaging in-plane and normal stresses near an interface crack using traction force microscopy. Proc. Natl Acad. Sci. USA 107, 14964 (2010).

41. Marchand, A., Das, S., Snoeijer, J. H. & Andreotti, B. Contact angles on a soft solid: from Youngs law to Neumanns law. Phys. Rev. Lett. 109, 236101 (2012).

42. Flory, P. & Rehner, J. Statistical mechanics of cross-linked polymer networks I: Rubberlike elasticity. J. Chem. Phys. 11, 512520 (1943).

43. Rubinstein, S. & Panyukov, M. Elasticity of polymer networks. Macromolecules 35, 66706686 (2002).

44. Voinov, O. V. Hydrodynamics of wetting [English translation]. Fluid Dyn. 11, 714721 (1976).

45. Cox, R. G. The dynamics of the spreading of liquids on a solid surface. Part 1. Viscous ow. J. Fluid Mech. 168, 169194 (1986).

Acknowledgements

We thank T. Baumberger, A. Eddi and E. Raphael for fruitful discussions.

S.K. acknowledges nancial support from NWO through VIDI Grant No. 11304. J.S. acknowledges nancial support from ERC (the European Research Council) Consolidator Grant No. 616918.

Author contributions

J.H.S. and B.A. designed the research; S.K., J.H.S. and S.D. developed the theory; M.G. and S.K. performed and analysed the wetting experiments; H.P. performed the rheological measurements; B.A. and J.H.S. wrote the manuscript. All authors discussed the results and commented for the nal manuscript.

Additional information

Supplementary information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions

Web End =reprintsandpermissions/

How to cite this article: Karpitschka, S. et al. Droplets move over viscoelastic substrates by surng a ridge. Nat. Commun. 6:7891 doi: 10.1038/ncomms8891 (2015).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

NATURE COMMUNICATIONS | 6:7891 | DOI: 10.1038/ncomms8891 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

& 2015 Macmillan Publishers Limited. All rights reserved.