Christiana E.
Udoh
,
Valeria
Garbin
and
João T.
Cabral
*
Department of Chemical Engineering, Imperial College London, London SW7 2AZ, UK. E-mail: j.cabral@imperial.ac.uk
First published on 4th June 2019
We report the formation of polymeric and nanocomposite capsules via droplet solvent extraction, focusing on the interplay between solvent exchange and removal, demixing and directional solidification kinetics. We investigate a model system of sodium poly(styrene sulfonate), NaPSS and silica nanoparticles in aqueous solution, whose phase behaviour is experimentally measured, and examine a series of selective extraction solvents (toluene, butyl acetate, ethyl acetate and methyl ethyl ketone), ranging from 0.04 to 11% v/v water solubility. Tuning the rate of solvent exchange is shown to provide an effective means of decoupling demixing and solidification timescales, and thereby tunes the internal microstructure of the capsule, including hollow, microporous, core–shell, and bicontinuous morphologies. In turn, these determine the capsule dissolution mechanism and kinetics, ranging from single to pulsed release profiles of nanoparticle clusters (at intermediate solubilities), to minimal dissolution (at either extremes). These findings provide facile design and assembly strategies for functional capsules with time-varying release profiles.
Common approaches for encapsulation include spray drying,8 emulsion polymerisation,9 photopolymerization,7 solvent displacement,10 and layer-by-layer assembly,11 as well as the more recent flash nanoprecipitation.12 The advent of droplet microfluidics, and single and multiple emulsification has provided exceptional capability for templating solute droplets, with relatively high throughput, enabling capsules to be formed by initiating polymerization, phase change, solvent evaporation/extraction or directed solidification.13–15 In solvent evaporation, solute-containing liquid droplets are concentrated by osmotic removal of the solvent in liquid or gas phase. In addition, their internal microstructure can be spontaneously generated by demixing, coarsening and eventual kinetic arrest upon solidification. This methodology has been employed in the fabrication of porous capsules with various architectures, including isolated micropores as well as interconnected domains,16,17 using multi-step processes. Single-step approaches able to spontaneously control capsule microstructure and triggered delivery of active agents are evidently attractive and have remained elusive.
We have previously demonstrated the formation of spherical and non-spherical polymeric particles with prescribed internal micropore structure and size distribution, via microfluidic extraction.18,19 The roles of polymer molecular mass Mw and water solubility (the good solvent), and starting composition (with respect to overlap and concentrated crossovers) and solution viscosity were investigated.20 In addition, we reported an approach for anisotropic polymer–nanocomposite capsule design, mapping accessible shapes and morphologies. Internally bicontinuous capsules were found to exhibit pulsed cargo release, in a series of bursts.21 Here, we seek to elucidate the role of the extraction solvent in modulating microcapsule architecture and, in turn, the release mechanism and kinetics of nanoparticle cargo. We therefore investigate a series of extraction solvents and carry out a comparative study of their effect in both capsule formation and dissolution. Our hypothesis is that the timescale for nanocomposite capsule formation can be decoupled into competitive solidification and demixing timescales, such that the microstructure of capsules could be tailored precisely by design of the extraction process.
With these ideas in mind, we first establish the solution thermodynamics for ternary systems, comprising solvent (water), non-solvent (MEK, EA, BA and T) and neat NaPSS or SiO2 by turbidimetry, as shown in Fig. 1A. The obtained binodal lines approximately follow the solubility limit of the neat solvent pairs (provided in ESI,† Table S1), providing a simple means to estimate the stability lines. The solubility of water in the various non-solvents ranges from ≃11 to 0.04% v/v, respectively from MEK to T. Demixing of the aqueous solutions by non-solvent addition is therefore expected to occur at progressively smaller concentrations of MEK, EA, BA and T as shown in Fig. 1A. Initial droplet compositions of interest are shown by the grey band (0.8–1) along the H2O:solute line on the phase diagram. The non-solvent quality appears thus to be an effective means of precisely selecting the location of the demixing trigger along the solvent extraction pathway of a droplet. However, given that solvent extraction will be correspondingly slower for poorer selective solvents, the interplay between demixing and extraction kinetics is a priori not obvious. Recent phase field simulations by Tree, Fredrickson and co-workers30 are able to qualitatively capture many features of the morphological evolution of non-solvent induced phase separation (NIPS) in a planar geometry, although not yet kinetic arrest upon coarsening, or two solute systems, such as those considered next, as relevant for encapsulation.
The phase behaviour of ternary mixtures of NaPSS/SiO2/H2O as a function of composition, at 20 °C, in terms of the polymer concentration Cp and nanoparticle volume fraction ϕNP is shown in Fig. 1B and C. Cross-sectional images of mixtures with constant ϕNP = 0.12 and increasing Cp (left to right) are provided in Fig. 1B, corresponding to the vertical line in Fig. 1C. Upon increasing Cp, initially homogeneous and stable mixtures become biphasic, exhibiting a bottom, liquid-phase, of higher ϕNP and a top, cloudy, gas phase. Coexisting three phases are observed upon further increasing Cp, corresponding to a bottom (solid), middle (liquid), and a top gas phase. At sufficiently high ϕNP and Cp, solid phases span the entire sample volume. The gas–liquid and solid–liquid interfaces were identified following the method by Faers and Luckham,31 by tilting demixed sample vials at a 30° angle and observing the interface tilt with respect to the vial, as well as its flow response; phase assignment was further confirmed by densitometry. Phase separation and coarsening in our system reached equilibrium after approximately two weeks, and the extent of the demixed regions did not change noticeably after 8 months. Crystal regions were not observed in our experiments probably due to the polydispersity of the silica dispersion and polymer chain lengths. It has been shown that even modest polydispersities suppress crystallization and fluid–solid transition, and introduce multiple solid phases,32,33 while gas–liquid demixing is less sensitive to polydispersity. We summarise our findings in a phase diagram, mapping out stable and immiscible compositions as shown in Fig. 1C.
The capsule formation process is expected to involve three simultaneous processes: solvent exchange and extraction, phase separation (via nucleation and growth or spinodal decomposition) and coarsening, and eventual phase inversion and solidification resulting in the formation of solid capsules, illustrated in the schematic presented in Fig. 2. Extraction is initiated by immersion of droplets in the non-solvent bath enabling solvent exchange at the interface, and a net osmotic removal of H2O, increasing both the polymer and/or nanoparticle concentration within the droplet. Non-solvent ingress within the droplet causes the solution to cross its phase boundaries, depicted in Fig. 1A, driving phase separation, beginning at the interface where the solvent composition gradient is highest, and continuing temporally as further extraction and internal recirculation flows occur. Eventually, kinetic arrest is expected when the solute concentration in one phase is sufficiently high to cause phase inversion and precipitation (as extensively employed in membrane fabrication30,34,35).
The pathway and kinetics of capsule formation are therefore expected to depend on an interplay between mutual affinity and diffusion between solvent and selective solvent, droplet concentration and dimensions, the rich phase behaviour of polymer–nanoparticle solutions, and eventual kinetic arrest. With these ideas in mind, we design drop shrinkage experiments to yield a range of effective Péclet numbers, estimated as Pe = [(dR/dt)R0]/D, where dR/dt is the initial shrinkage rate, R0 is the initial droplet radius, and D is the Stokes–Einstein diffusivity of the solute (detailed in ESI†). This dimensionless mass transport number characterises the relative rates of diffusion and convection, and ranges from Pe ≃ 0.9 to 60, for T to MEK, depending on solute content. Spray dried36–38 and solvent-extracted,39 single- and multi-component droplets can be expected to yield dimpled, hollow, or crumpled capsules under rapid extraction (Pe ≫ 1), while dense compact capsules have been reported for comparatively slow extraction (Pe ≪ 1).
At all conditions, droplet radii decrease monotonically with time and eventual solidification results in capsule formation which, under these conditions, remain largely spherical. Internal convection is observed within droplets, in particular for the faster extractions (ESI,† Supplementary movie S1). The external morphology of capsules obtained using MEK appears somewhat crumpled, which is expected for fast extractions (higher Pe) and shell formation at earlier times,21 followed by subsequent buckling upon continuous solvent removal. Buckling is driven by the unconstrained decrease in compressive stresses generated on the shell wall during droplet volume reduction.40
Notably, the timescales for extraction and solidification vary by more that two orders of magnitude from MEK to T, following the variation in H2O miscibility, depicted in Fig. 1A. The droplets exchanged with the solvents with intermediate miscibilities, EA and BA, exhibit clear demixed structures during extraction (ESI,† Supplementary movies S2 and S3). Fig. 3E shows a clear spinodal structure, which appears transiently and coarsens around t ∼ 140 s for an EA extraction. According to the phase diagram (Fig. 1C), the two internal phases are assigned to nanoparticle-rich and -poor, as expected for a quench into the unstable gas–liquid region, resulting in bicontinuous domains (Fig. 3E and F) as previously reported.29 In polymer–colloid demixing, the polymer component has been found to preferentially partition to the colloid-poor phase.27 We therefore assign the bicontinuous matrix phase and capsule skin as originating from the polymer-rich (colloid-poor) gas-phase, and the internal bicontinuous phase to the nanoparticle-rich liquid phase, validated by EDX measurements discussed below.
While the process is qualitatively similar in BA, spinodal demixing now occurs around t ∼ 1000 s (shown in Fig. 3F); further, given the lower H2O:BA miscibility, extraction takes considerably longer and thus coarsening of the structure can happen over longer timescales. Eventually, a well-defined skin layer (shell thickness of ∼10 μm) emerges upon approaching the solidification stage, depicted in Fig. 3G, which we interpret as resulting from the comparatively prolonged phase coarsening and stratification of the polymer-rich phase to the interface.
For the fastest extraction (MEK), a bicontinuous phase is not observed, which we interpret as due to the high Pe number resulting in rapid skin formation and precipitation, likely circumventing the spinodal pathway and leading to a hollow morphology examined below. At the other extreme of slow extraction, associated with the low solubility of H2O in T miscibility, and Pe ≃ 0.8, we obtain compact, dense particles with no visible demixed structure (ESI,† Supplementary movie S4), interpreted as due to the slow shrinkage aided by recirculation flows that enable gradual equilibration until kinetic arrest.
These observations are quantified by the evolution of droplet radius R(t) during the extraction process as shown in Fig. 4. The lines are empirical fits to , where R0 is the initial droplet radius, R∞ is the (final) capsule radius, τ is the extraction time (when R ceases to change), and α is a parameter accounting for non-Fickian diffusion, introduced previously.18,19 In agreement with the expected impact of solvent solubility on dissolution profile,41–43 extraction timescales ranging from tens to thousands of seconds can be readily tuned by the extraction medium. Further, the extraction timescale τ is found to depend linearly on R0 (ESI,† Fig. S3) with a prefactor which depends, as expected,18–20 on solute content.
Fig. 4 Droplet extraction kinetics in selected solvents, methyl ethyl ketone (MEK), ethyl acetate (EA), butyl acetate (BA), and toluene (T). Evolution of droplet radius with time for (A) pure H2O, (B) 1% (w/v) NaPSS, and (C) 1% w/v NaPSS + 10% w/v SiO2. Initial droplet radius R0 for data shown ranges between 250–270 μm. In (B and C), the radius is normalised by R0 to facilitate comparison. The raw data are provided in ESI,† Fig. S2. The open and closed symbols in (C) capture minor deviations from sphericity, by plotting both dimensions of the spheroid. Lines shown are empirical fits to the data, described in the text. |
Similar trends are observed with neat 1% w/v NaPSS (Fig. 4B) and 1% w/v NaPSS + 10% w/v SiO2 (Fig. 4C) droplets. Obviously, in these cases, instead of full dissolution and R(t) → 0, a final particle or capsule size R∞ is obtained. This size follows a dependence with the opposite sign to that of τ, and smaller (thus more compact) particles are attained with poorer extraction solvents, as demonstrated above.
A summary of these observations is presented in Fig. 5 in terms of final capsule size R∞ and extraction time τ. We find that droplet extraction kinetics τ are remarkably governed by solvent/non-solvent miscibility (summarised in ESI,† Table S1), regardless of solute (neat solvents, polymer or polymer/nanoparticle solutions) which has a comparably smaller effect. From droplets of comparable initial radius (selected in Fig. 5A for R0 ≃ 250–270 μm), the resulting capsule size R∞ also depends monotonically, approximately linearly, on H2O solubility, with minor variations likely due to surface and crowding droplet effects during extraction. Further, R∞ also scales monotonically with polymer and/or nanoparticle content.18,19 For reference, raw data are provided in ESI,† Fig. S2. Therefore, it emerges that R0, solute content and, to some extent solvent/non-solvent miscibility, largely control R∞, while the extraction time scale τ is largely governed by solvent/non-solvent miscibility, with solute content playing a comparatively minor role, and τ ∝ R0, the initial droplet size. The inverse relation between R∞ and τ can be qualitatively rationalized in the context that larger capsules form upon faster extractions, due to an early onset of skin formation. However, a theoretical framework to describe the phenomena requires the coupling of extraction (and mutual diffusion), demixing and directional solidification and is currently not available beyond the descriptive model introduced above and parameterized by these experiments. These parameters provide a facile strategy to design capsule size and formation timescale by osmotic solvent extraction, while the demixing processes described are expected to impact the internal microstructure.
Fig. 5 Dependence of (A) capsule radius, R∞, and (B) extraction time, τ, normalised by R0, on the solvent/nonsolvent solubility limit for neat polymer and polymer/nanoparticle droplets. Data shown in (A) include only droplets with initial droplet radius R0 ≃ 250–270 μm. Lines shown are guides to the eye. Selective solvents methyl ethyl ketone (MEK), ethyl acetate (EA), butyl acetate (BA), and toluene (T) were investigated. Additional data showing the linear dependence of τ on R0 for NaPSS and NaPSS:SiO2 droplets immersed in EA are provided in ESI,† S3. |
The counterpart internal structures of polymer–nanocomposite capsules are shown in Fig. 7 as a function of extraction solvent. In this case, the morphology is not only determined by the polymer-rich and -poor phase, but also by the demixing between polymer and nanoparticle phase in solution. Upon decreasing extraction solvent quality a succession of capsule microstructures is found, namely hollow, bicontinuous, core–shell, and compact, in good agreement with optical microscopy observations. SEM and EDX (show in ESI,† Fig. S4) measurements reveal that the internal spherical clusters, of ≈1–5 μm in diameter, are precipitated silica aggregates, bound by trace amounts of polymer. Varying extraction time is found to effectively tune structural coarsening time and thus internal architecture, demonstrated by the formation of a bicontinuous (Fig. 7B) or core–shell capsule (Fig. 7C) achieved using EA and BA respectively as non-solvents. At the extremes of solubility, hollow (larger) capsules and dense (compact) capsules are obtained for the fastest and slowest extraction rates, in MEK and T, respectively (Fig. 7A and D).
Finally, we examine the consequence of capsule morphology on nanoparticle release. The dissolution of nanocomposite capsules in DI H2O is shown in Fig. 8. Minimal release of nanoparticle clusters is observed from capsules extracted in MEK and T, as shown in Fig. 8A and D, within the timescales probed up to 104 s, and this intensity could not be resolved optically. Internally bicontinuous capsules formed by EA extraction (Fig. 8B), are found to release the encapsulated nanoparticle clusters in a series of bursts from specific sites around the capsule surface, over long timescales (ESI,† Supplementary movie S5), as shown in Fig. 8E. During release, the capsule's scaffold retains its integrity, thereby enabling the multiple bursts from internal compartments over time.21 Core–shell capsules, formed by BA extraction, possessing a large single compartment filled with nanoparticle clusters, releases from only one active site and in a single burst, at a comparatively shorter time (ESI,† Supplementary movie S6). We estimate nanoparticle cluster release from active sites of EA and BA capsules (indicated in Fig. 8B and C) and the overall release by computing the integrated release density over time, as shown in Fig. 8E and F. We interpret the minimal release observed in the most compact and least compact capsules as due to the formation of dense polymer–nanoparticle composites: in the case of fast extraction, at the capsule shell, while for slow extraction within the bulk volume. Short-range attractions between nanoparticle clusters are expected to hinder dissolution, in the absence of strong agitation. However, at intermediate solubilities, likely due to the gradually stratified structure and larger inter-particle distance, capsules exhibit tunable release profiles. These observations corroborate the proposed correlation between internal architecture and release of ‘cargo’, following from the simple design strategies outlined for capsule formation.
To examine the functional consequences of this design, capsules were subsequently immersed in good solvent (water). At the two extreme solubilities (MEK, T) no appreciable dissolution takes place over long times (≥24 h) under quiescent conditions (i.e. in the absence of agitation), interpreted as due to the formation of a compact skin or compact (bulk) particle in MEK and T, respectively. At the intermediate solubilities, with BA we find a single burst release at early times, compatible with the formation of a single large, particle reservoir, due to coarsening of the internal morphology; the bicontinuous structure formed by EA extraction (at faster timescales), releases over a series of pulses, and from multiple locations, reflecting its internal compartmentalized morphology. Polymer–nanoparticle thermodynamic and tuneable solvent extraction kinetics are therefore also shown to provide a means to modulate the dissolution of capsules. Further control can likely be exercised by employing solvent mixtures, across a large R0 range (typically 1–1000 μm) available via microfluidics.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9sm00708c |
This journal is © The Royal Society of Chemistry 2019 |