Kathryn A.
Rosowski
a,
Estefania
Vidal-Henriquez
b,
David
Zwicker
b,
Robert W.
Style
*a and
Eric R.
Dufresne
*a
aDepartment of Materials, ETH Zürich, 8093 Zürich, Switzerland. E-mail: robert.style@mat.ethz.ch; eric.dufresne@mat.ethz.ch
bMax Planck Institute for Dynamics and Self-Organization, 37077, Göttingen, Germany
First published on 4th June 2020
When liquid droplets nucleate and grow in a polymer network, compressive stresses can significantly increase their internal pressure, reaching values that far exceed the Laplace pressure. When droplets have grown in a polymer network with a stiffness gradient, droplets in relatively stiff regions of the network tend to dissolve, favoring growth of droplets in softer regions. Here, we show that this elastic ripening can be strong enough to reverse the direction of Ostwald ripening: large droplets can shrink to feed the growth of smaller ones. To numerically model these experiments, we generalize the theory of elastic ripening to account for gradients in solubility alongside gradients in mechanical stiffness.
Fig. 1 Ostwald ripening in liquids and elastic ripening in polymer networks. (A) and (B) Schematic diagrams summarizing the key differences between Ostwald ripening and elastic ripening.3 Here, Pc ∼ E is the cavitation pressure. (C) In substrates with a gradient in stiffness (here a step change from E = 10–700 kPa), but with otherwise identical properties, droplets ripen from the stiffer side to the softer side (data from ref. 3). |
This picture changes significantly when droplets form by nucleation and growth in a polymer network. In homogeneous networks, these droplets tend to be monodisperse and very stable, with smaller droplets appearing in stiffer networks.4 When droplets exclude the polymer network, they push the network outward as they grow. In response, the network squeezes the droplets. This both suppresses droplet condensation,3 and increases the droplet's internal pressure by an amount comparable to the network's Young modulus, E.5–7 This increase in pressure can potentially far exceed the Laplace pressure. Thus, when the polymer network has heterogeneous mechanical properties, the elastic contribution to droplet pressure is heterogeneous and can drive transport of material from droplets in stiff regions to droplets in soft regions, summarized in Fig. 1B.3 Like Ostwald ripening, this ‘elastic ripening’ is mediated by the transport of material between droplets in the dilute phase. Related phenomena have been observed in the nucleus of living cells.8
Previous demonstrations of elastic ripening, like that in Fig. 1C, have superficially resembled Ostwald ripening, because large droplets (in soft regions of the network), grew at the expense of small droplets (in stiff regions of the network). Thus, the establishment of elastic stresses as the driving force for ripening relied on comparisons across experiments with varying elastic moduli, and comparisons to numerical models.
Here, we demonstrate that elastic ripening phenomena can qualitatively differ from classic Ostwald ripening. In particular, we demonstrate the growth of small droplets, fed by the dissolution of large ones. Experimentally, this is achieved by coupling two different families of silicone gels, which have different thermodynamic and transport properties at the same elastic modulus. These new results demand a generalization of the model for elastic ripening. In particular, gradients in solubility must be accounted for. With this generalization, numerical simulations capture essential features of experimental observations, using experimentally measured properties.
To create stiffness gradients, we make two different silicones side by side. We first cure a 3–5 mm layer of the stiffer silicone in a polystyrene Petri dish (Greiner). Half of this is cut out with a scalpel, pulled off the dish, and placed in half of a new, glass-bottomed dish (MatTek). The softer silicone is then poured into the other side of the dish and allowed to cure.3 When we drive phase separation by cooling, droplets on each side of the sample grow to a uniform distribution, with a sharp transition between them (Fig. 1C). After relatively fast droplet formation, we observe slow evolution of the droplets near the interface.
In previous experiments,3 we observed that smaller droplets, on the stiffer side, dissolved and fed the growth of larger droplets, on the softer side (e.g.Fig. 1C). While resembling familiar Ostwald ripening, due to surface tension, we concluded that ripening was driven by gradients in the stiffness of the polymer network. This conclusion was based on an observed increase of the ripening rate with the stiffness difference and comparison with numerical models. Here, we challenge this interpretation by considering ripening across a broader range of silicone samples.
We use two different families of silicone gels, as described in the Materials and Methods. ‘Gelest’ silicones are fabricated by mixing together functionalized silicone polymers, crosslinker, and a platinum catalyst following.3 ‘Sylgard’ silicones are fabricated with a popular commercial kit, Sylgard 184. Chemically, these are essentially the same as the Gelest silicones, but additionally contain a significant quantity of silica nanoparticle filler.9 With both types of silicone, we can achieve a range of Young's moduli, E, from tens to hundreds of kPa.
Despite some similarities, these two families of silicones show quantitative differences in their phase-separation behavior. In particular, condensed droplets of fluorinated oil have different sizes at the same network stiffness, as shown in Fig. 2. As the stiffness of Gelest networks increases from 10 to 580 kPa, the mean droplet radius reduces from 14.9 ± 1.4 μm to 5.7 ± 0.7 μm. As the stiffness of the Sylgard networks increases from 80 to 450 kPa, the mean droplet radius reduces slightly from 8.3 ± 0.6 μm to 6.0 ± 0.5 μm. Note that these radii are the final, stable sizes that droplets grow to when the homogeneous silicones fully cool to room temperature after saturation at 40 °C.
In samples made from only one of these silicone families, elastic ripening and Ostwald ripening proceed in the same direction. This is because the droplet size generically decreases with the stiffness. With two different families of silicone, we can remove this limitation, and explore combinations of materials where elastic and Ostwald ripening alternatively reinforce or oppose each other.
An example where Ostwald and elastic ripening reinforce each other is shown in Fig. 3A and Movie S1 (ESI†). Here, stiffer Sylgard silicone (E = 80 kPa, mean droplet size 8.4 μm) is in contact with softer Gelest silicone (E = 10 kPa, mean droplet size 14.9 μm). Consistent with previous results, smaller fluorinated oil droplets on the stiff side shrink while feeding the growth of larger droplets on the soft side. The direction of ripening is consistent with both Ostwald and elastic ripening.
An example where elastic forces oppose Ostwald ripening is shown in Fig. 3B and Movie S2 (ESI†). Here, stiff Gelest silicone (E = 140 kPa, mean droplet size 12.7 μm) is in contact with soft Sylgard silicone (E = 100 kPa, mean droplet size 7.3 μm). In this case, the larger droplets near the interface on the stiff side shrink while small droplets on the soft side grow. This is the opposite of the trend expected from Ostwald ripening, and provides a convincing visual case for elastic ripening.
An interesting complication in the interpretation of the data in Fig. 3 arises due to the differences in the solubility of fluorinated oil on the two sides. In both silicones, the saturation concentration, ϕsat, increases with the temperature, but does not change significantly with the elastic modulus (Fig. 4). However, fluorinated oil has significantly higher solubility in the Gelest than Sylgard silicones, Fig. 4. Additionally, the solubility of fluorinated oil is more temperature sensitive in Gelest than in Sylgard. In practice, this means that when the temperature drops, the Gelest samples become more strongly supersaturated than the Sylgard ones. As we will see below, these factors play important roles in understanding the observed ripening behaviour.
In our experiments, solute transport occurs in gels with a heterogeneous saturation concentration. In such cases, one needs to take a more general approach, where fluxes are driven by gradients in chemical potential, μ:
(1) |
μ = kBTlog(ϕ/ϕsat). | (2) |
In the droplet phase, the chemical potential is simply related to the pressure, μ ≈ P/nL.13 Here, nL is the number density of molecules in the droplet phase. Differences in chemical potential between droplets are thus equivalent to pressure differences, and eqn (1) matches our expectation that solute should be driven from high pressure to low pressure.
For droplets in an elastic network, the pressure within the droplet has the form,
(3) |
(4) |
With this simple extension of ripening theory, numerical simulations capture essential features of the experiments, as shown in Fig. 3. We randomly place droplets on either side of the sample, with density and size distributions matching the measured experimental results. Since the droplets equilibrate quickly with their surrounding, we use the chemical potential associated with the pressure P given by eqn (3) to calculate the concentration ϕ using eqn (2). Finally, we model diffusive flux between droplets by combining eqn (1) and (2) with a conservation law, as described in the Materials and Methods. Values of D, ϕsat and E are taken directly from the measured values.
The results, shown in Fig. 5 and Movies S3 and S4 (ESI†), are in qualitative agreement with the experiments. The left-hand column shows results of the simulation of the experiment of Fig. 3A. As in the experiment, we see depletion of the droplets in the stiffer Sylgard, and growth of droplets in the softer Gelest adjacent to the interface. Interestingly, we see that ϕ is always larger on the softer side than on the stiffer side, so diffusion progresses up concentration gradients. This is due to the differences in ϕsat between the two sides. When we account for this in the chemical potential, viaeqn (2), we verify that transport occurs down chemical potential gradients, as expected. The right-hand column of the figure simulates the experiment of Fig. 3B. Again, we see the ripening from the stiff side to the soft side. However, in this case, ripening occurs from larger droplets to smaller droplets – the opposite direction to what is usually expected.
Fig. 5 Numerical simulations of the experimental setups presented in Fig. 3. Top: Droplet radius. Middle: Concentration in the dilute phase, ϕ. Bottom: Chemical potential, μ. (A) For 80 kPa Sylgard next to 10 kPa Gelest, elastic ripening moves against concentration gradients. (B) For 140 kPa next to 100 kPa Sylgard, elastic ripening goes against classical Ostwald ripening. |
While the simulations capture the essential experimental trends, transport in experiments is approximately ten times faster than expected from simulations. The origin of this discrepancy remains elusive.
These results help to lay the foundations for the analysis of phase separation in complex, heterogeneous environments. Our experiments and theory are inspired by recent observations of phase separation in living cells (e.g.ref. 8 and 16). In that context, specific macromolecules, including nucleic acids and proteins, segregate into functional domains within the cytoplasm and nucleoplasm.17–20 The working model for interpreting these phenomena is the phase separation of two-component fluids. While this captures some of the basic phenomenology, recent theories aimed to account for activity,21,22 large-scale concentration gradients,23,24 complex compositions,25 and rheology.4,11,26–28 Our work further establishes a general framework for evaluating the role of a passive network in crowded systems (see also ref. 29). Note that we assume that our system remains near equilibrium and that the rheology of the continuous phase can be reduced to a single, static, inflation pressure. Thus, these results should not be directly applied to complex living systems without caution.
We also anticipate that these results may be relevant to several related systems. For example, this system has much in common with the volumetric phase separation that can occur when crosslinked hydrogels go through deswelling transitions,30–33 giving rise to phase-separated regions of collapsed and non-collapsed gel. In this case, our results suggest the possibility of interesting ripening behaviour in the presence of network heterogeneities. Similarly, there are strong ties to the stabilisation of foams by bulk or interfacial elasticity,34,35 temporary elastic stabilisation during viscoelastic phase separation,26 and in the design of other phase-separation processes in materials near equilibrium. These include the production of porous membranes and scaffolds,36 the formation of segregated ice during the processing of frozen foods (e.g.ref. 37), in cryopreservation,38 and other processes where freezing causes material damage to porous materials (e.g.ref. 39).
To create ‘Sylgard’ silicones, we use the commercial brand Sylgard 184 (Dow). The base is mixed with curing agent, using different ratios to attain different stiffnesses. The mixture is degassed thorougly and cured at 40–60 °C for at least one week.
For both silicones, the Young's modulus, E, is measured by indentation following the procedure in ref. 40. We slowly indent bulk samples with a cylindrical, flat-punch indenter with a 1 mm contact radius, and calculate E from the initial slope of the force–indentation curve. We never see creep, and there is no viscous dissipation and no frequency-dependence of E below indentation speeds of 0.1 mm s−1.
The diffusion coefficient, D, was calculated following.3 Briefly, for an infinitely thick layer of silicone, covered in a layer of fluorinated oil, we expect the concentration profile (in terms of mass per volume) to follow
(5) |
moil(t) = 2rdish2ceq(πDt)1/2 | (6) |
Fig. 6 Measuring the diffusivity of fluorinated oil in silicone. (A) Here, Sylgard silicone with E = 230 kPa is immersed in fluorinated oil and weighed at regular intervals to determine the uptake in mass, as shown. Initially the uptake follows moil ∼ t1/2, as expected from eqn (6). Subsequently, moil plateaus at the saturation concentration. The red dashed line shows a fit with slope 1/2. (B) Plots of diffusivity vs. E show how D is larger in Gelest silicones than in Sylgard silicones. D is also independent of stiffness for the two silicone families, and increases with temperature. |
(7) |
ϕcond = ϕsatexp(P/nLkBT). | (8) |
(9) |
We initialize the system by randomly placing droplets according to experimental density measurements.4 The droplet radii are randomly drawn from a Gaussian distribution following the data shown in Fig. 2. The dilute phase ϕ is initialized at ϕ = ϕcond to avoid initial droplet growth. In our simulations we assumed a densely packed droplet phase ϕin = 1, other parameters are α = 0.5, nLkBT = 11 MPa and ϒ = 4 mN m−2.
Footnote |
† Electronic supplementary information (ESI) available: Supplementary Movie 1: elastic ripening from small to big. Movie showing progression of the ripening experiment in Fig. 3A. Images are taken at 5 minutes intervals. Supplementary Movie 2: elastic ripening from big to small. Movie showing progression of the ripening experiment in Fig. 3B. Images are taken at 5 minutes intervals. Supplementary Movie 3: simulated elastic ripening from small to big. 2D projection of a numerical simulation with Gelest (ϕsat = 0.023 and E = 10 kPa) on the left, and Sylgard (ϕsat = 0.022 and E = 80 kPa) on the right. Chemical potential μ is shown as a density plot with colorbar at the right side. Supplementary Movie 4: simulated elastic ripening from big to small. 2D projection of a numerical simulation with Gelest (ϕsat = 0.023 and E = 140 kPa) on the left, and Sylgard (ϕsat = 0.022 and E = 100 kPa) on the right. Chemical potential μ is shown as a density plot with colorbar at the right side. See DOI: 10.1039/d0sm00628a |
This journal is © The Royal Society of Chemistry 2020 |