Miriam
Simon
*a,
Michael
Gradzielski
*a and
Ingo
Hoffmann
*b
aStranski-Laboratorium für Physikalische und Theoretische Chemie, Institut für Chemie, Technische Universität Berlin, Straße des 17. Juni 124, Sekr. TC 7, D-10623 Berlin, Germany. E-mail: miriam.simon@tu-berlin.de; michael.gradzielski@tu-berlin.de
bInstitut Max von Laue-Paul Langevin (ILL), 71 avenue des Martyrs, CS 20156, F-38042 Grenoble Cedex 9, France. E-mail: hoffmann@ill.fr
First published on 20th August 2020
Oil-in-water (O/W) microemulsion droplets are convenient carriers for hydrophobic molecules in an aqueous phase and are used for a wide range of applications. We studied weakly charged O/W microemulsion droplets complexed with oppositely charged polyacrylates that form long linear arrangements of droplets. All samples showed rather low viscosities, which is in contrast to similar systems of hydrophobically interconnected droplets. Here, we applied small-angle neutron scattering, dynamic light scattering and neutron spin-echo spectroscopy to characterise the dynamic properties of polyacrylate/microemulsion complexes in order to understand the origin of the low-viscous behaviour. We found that the electrostatic interactions lead to very dynamic complexes with high exchange rates of droplets and only a fraction of the droplets is contained within the transient complexes at a given time. These results were only accessible by the combination of different methods as one method alone would have given an incomplete picture.
Microemulsions have been known since the 1940s,5 they occur as finely dispersed droplets of water in oil (W/O) or oil in water (O/W) but also bicontinuous structures are possible.6 Microemulsions uniquely combine a number of properties, like thermodynamic stability, transparent appearance, low viscosity and spontaneous formation which makes them very interesting systems whenever it is necessary to disperse hydrophobic molecules in a hydrophilic solvent or vice versa as it is needed for example in drug delivery, nanoparticle synthesis, oil recovery or cosmetics and detergents.7–11 In order to tailor microemulsion properties to specific needs for applications (like enhanced viscosity), it might be necessary to add additives like polymers.
Complexes formed by microemulsion droplets and polymers,12,13 including neutral amphiphilic polymers,14,15 block copolymers16,17 or telechelic polymers18,19 have already been studied by several groups. In most cases a gelling of the sample due to network formation of the microemulsion droplets with the polymer was found. In contrast, little research has been done on complexes formed by ionic microemulsion droplets and oppositely charged polyelectrolytes, where the interaction is mainly of electrostatic nature.20 For cationic surfactant/anionic polymer complexes it was found that the solubilisation capacity of the micelles increases in presence of the polymer, due to the additional ionic strength.21 If charged microemulsion droplets bind to polymers through electrostatic and/or hydrogen bonds, a bridging of the droplets would be possible as well as a decoration of the polymers.
In a recent study we have therefore investigated complexes formed by cationic O/W microemulsion droplets with polyacrylates.22 The structurally most interesting complexes were found in polyelectrolyte excess close to the phase boundary, where small-angle neutron scattering (SANS) measurements showed long cylindrical arrangements of droplets on the polyacrylate chains with 3 to 23 droplets per aggregate, depending on the Mw of polyelectrolyte. When increasing the polyelectrolyte content, thus moving further away from the phase boundary, the complexes become smaller and less extended. Surprisingly, all of these samples showed rather low viscosities (see Fig. S1†). This is in contrast to previously studied systems of very similar (but uncharged) microemulsion droplets that formed highly viscous networks when interacting with telechelic polymers.18,19 The telechelic polymers bind to the droplets via hydrophobic forces with 12–18 kBT per endcap, depending on the length of the alkyl chain, while the electrostatic interactions due to opposite charges become comparatively weak already at moderate salt concentrations.23 However, previous static measurements were not able to explain the low-viscous nature of the electrostatically formed polyelectrolyte/microemulsion aggregates. For this reason we were now employing a high-resolution dynamic method, namely neutron spin-echo (NSE), to address this question in thorough detail. In that context, it might be noted that NSE was already successfully employed by Yearley et al.24 to study the relation of cluster formation and solution viscosity.
Sodium polyacrylate (NaPA) was employed as polyelectrolyte. It was either purchased directly or prepared by adding a stoichiometric amount of NaOH to the poly(acrylic acid) (PAA). Four different molecular weights (Mw) were used in this work: 5.1 kg mol−1 (sodium salt, used as received from Sigma Aldrich), 15 kg mol−1 (sodium salt solution from Sigma Aldrich, freeze dried before use), 60 kg mol−1 (sodium salt, used as received from Fluka Chemika) and 240 kg mol−1 (PAA; solution from Acros Organics, prepared by adding a stoichiometric amount of NaOH (≥98%, p.a., ISO, Carl ROTH) to the solution and freeze dried before use, which results in an average molecular weight of 315 kg mol−1). These molecular weights correspond to contour lengths (lengths of the stretched polyelectrolyte) of 140, 400, 1600, and 8500 Å respectively, the corresponding overlap concentrations for these polymers (assuming fully stretched chains) are 32, 3.9, 0.25, and 0.007 mM (in monomeric units).
Samples were prepared at a microemulsion concentration of 100 mM surfactant (with 95 mol% TDMAO and 5 mol% TTAB) with 12 mM of added NaPA (monomeric units) to obtain a charge ratio of z = [−]/([+] + [−]) = 0.7, where [+] and [−] are the nominal number of charges of employed TTAB molecules and NaPA monomer units respectively. A charge ratio of z = 0.7 corresponds to samples at polyelectrolyte charge excess but close to the phase boundary where long-time stable samples containing elongated aggregates were found in previous studies.22 The respective concentrations of surfactant and NaPA can be used to calculate the number of droplets and total length of PE chain. The ratio of both numbers results in an average length of 140 Å PE per ME droplet of R = 45 Å, so the space between droplets would be a bit smaller than one droplet diameter. Samples for NSE and SANS measurements were prepared in D2O (>99.5% D, Eurisotop), while samples for DLS were prepared in H2O from a Millipore system.
Neutron spin-echo (NSE) measurements were performed on the instrument IN1525,26 at Institut Laue Langevin (ILL, Grenoble, France) using wavelengths of 6, 10 and 14 Å reaching Fourier times of 42, 194 and 531 ns and covering a q range from 0.019 to 0.14 Å−1. Standard methods were applied for data reduction.
Small-angle neutron scattering (SANS) experiments were carried out on PAXY at Laboratoire Léon Brillouin (LLB, Saclay, France) and on V4 at Helmholtz-Zentrum Berlin (HZB, Germany). Data reduction was done using the software package BerSANS.27
All measurements were carried out at 25 °C. See ESI for more details.†
![]() | ||
Fig. 1 SANS curves for microemulsion droplets with polyacrylates of different Mw at z = 0.7. Solid lines are best fits of a model of coexisting spheres and wormlike chains, see the ESI for details.† The model was chosen after analysis of the NSE data. Pink areas indicate q regions of DLS and NSE measurements. |
NSE measurements were performed to gain information about the dynamics of the system at nanometre length scales and nanosecond timescales. At first, the measured intermediate scattering functions were fitted with a simple exponential S(q,t)/S(q,0) = exp(−Dappq2t), (magnitude of the scattering vector q, Fourier time t, see Fig. S9 to S12†) to obtain the apparent diffusion coefficients Dapp shown in Fig. 2. The peak at 0.1 Å−1 is due to undulation motions of the microemulsion membranes, which are mostly visible at the form factor minimum, located at qminR = tan(qmin × R), R ≈ 4.5/qmin for spheres. This means that the microemulsion droplets are about 45 Å in radius. This peak is in good agreement with the measured SANS data which shows a form factor minimum at the same q (see Fig. 1). The diffusion coefficient corresponding to a sphere with a radius of 45 Å in a non-interacting sample at our experimental conditions would be , which is significantly more than we observe at q < 0.1 Å−1 for the pure microemulsion meaning that there are relevant effects from interparticle interactions resulting in a static structure factor S(q) with a peak at about 0.04 Å−1 (see Fig. S13†). The q dependence of Dapp is given by Dapp = D0 × H(q)/S(q), where the hydrodynamic function H(q) ≈ 1 in our concentration range,29 and the strong increase of Dapp at the lowest q values is due to S(q) < 1.
The characteristic behaviour of microemulsion droplets with a peak in Dapp at the form factor minimum is also observed for the polyelectrolyte/microemulsion complexes but their Dapp at lower q is significantly lower than for the pure microemulsion, showing that the dynamics are slower for the complexes (see encircled region in Fig. 2). The Dapp of the different complexes is rather similar regardless of the molecular weight of the NaPA, despite the different structures that are observed using SANS. A closer inspection of the intermediate scattering functions at q values where only diffusion should be visible reveals that S(q,t) is not single exponential for the complexes as opposed to the pure microemulsion (see Fig. S2†). This is in contrast to our findings from DLS where the data of the complexes can be described well as a monomodal decay (see Fig. 3 and S3 to S5†).
Such a bimodal diffusion process as a function of length and timescale can be observed if particles are entrapped within a transient network (and one sees fast movement within the mesh size and slower one across the network) or if some of the droplets are bound within the complex while others are freely moving. Assuming the later scenario to be more likely in our case, this leaves two possible explanations for the monomodality of the DLS data. Either the complexes show much stronger scattering than the free microemulsion droplets at the low q values of the DLS due to their larger size, rendering the latter virtually invisible. In this case, the diffusion coefficient measured in DLS would be at least as low as the slow diffusion coefficient measured in NSE unless there were strong repulsive interactions between the complexes to which the SANS data do not give any evidence. Alternatively, the complexes could be highly dynamic with a lifetime longer than the nanosecond scale of NSE but shorter than the micro- to millisecond timescale of DLS. For a system where particles change between two different diffusion coefficients D1,2 (D1 < D2), on a timescale τ in the limit of fast exchange (τ ≪ 1/(Dq2)), a single averaged relaxation rate is observed:
S(q,t) = exp(−(xD1 + (1 − x)D2)q2t) | (1) |
S(q,t) = x![]() | (2) |
To address this question in more detail, we proceeded to fit the NSE data from the complexes with a combination of the Milner–Safran model31,32 to account for the membrane undulations at high q and the translational diffusion of the free microemulsion droplets and a second slow diffusive mode with q independent diffusion coefficient Dslow and q dependent amplitude xslow, see eqn (S14).† The contribution from the free microemulsion droplets (membrane undulations and translational diffusion) was fixed for all samples according to fits to the data from the pure microemulsion (see Fig. S6†), where we obtained a value for the bending rigidity κ of 4.5 kBT, which is in good agreement with previous results on similar systems.33 Different approaches to describe the interactions in this system were tried by employing different approximations for the structure factor but it turned out, that different treatments had a rather small effect on the outcome (see Fig. S7 and S8† for a comparison). Thus we decided to fit our data treating the NaPA as simple electrolyte and leaving the volume fraction at its nominal value. The fitting procedure is described in detail in the ESI and fits are shown in Fig. S9 to S12.†
Fit results for xslow are shown in Fig. 4 in a q range where the influence from the membrane undulations is weak. Even though xslow is an intensity weight and does not directly correspond to the fraction of microemulsion droplets bound to the complexes, the intensity ratio between a single microemulsion droplet and a single complex using the same droplets as building blocks should be roughly proportional to the volume ratio of the structures in the NSE q range, as the respective form factors show similar intensity. Therefore, xslow should be a reasonably good approximation for the fraction of the droplets bound in complexes, which is about 20 to 30% of the microemulsion droplets for all complexes. The results are quite similar for all samples and do not show a pronounced trend with q. The small increase of xslow at low q can be explained by the intensity of the complexes increasing relative to the intensity of the pure microemulsion droplets towards lower q.
![]() | ||
Fig. 4 x slow obtained from fitting NSE curves with eqn (S14) and (S15).† No significant differences between samples can be seen, while a slight increase with decreasing q can be observed. |
The diffusion coefficients obtained from DLS are systematically higher than Dslow obtained from NSE (see Fig. 5). At low q, the relative intensity of the complexes compared to the free microemulsion droplets, increases drastically which affects the comparison between DLS and NSE (compare Fig. 1). The fact that the single diffusion coefficient seen in DLS is faster than the slow diffusion coefficient seen in NSE indicates that the result obtained from DLS is an average of complexes and free microemulsion droplets as described by eqn (1). This means that the monomodality of the DLS curves is due to the averaged signal, which arises from the short lifetime of the complexes. With increasing Mw, the complexes become larger, resulting in a higher intensity at low q (Fig. 5, right axis). The relative intensity of the complexes in the DLS q range increases as well, causing them to preponderate in the measured average, which explains the decrease of the DLS diffusion coefficient with increasing Mw. However, the local binding and the binding probability are independent of the length of the NaPA (Fig. 4), confirming that the droplet dynamics in these complexes is exclusively governed by the local interaction between droplet and chain fragment.
![]() | ||
Fig. 5 D slow obtained from fitting NSE curves with eqn (S14),† diffusion coefficients obtained from DLS and I(0) values obtained from SANS as a function of the molecular weight of the PE. While the observed diffusion coefficient is constant for NSE, it decreases with increasing Mw for DLS where at high Mw the scattering signal is dominated by the complexes and the relative contribution of free ME droplets becomes negligible. |
With this knowledge it was possible to refit the SANS data using a two-component fit, including homogeneous wormlike chains and free spheres with the respective intensity fractions obtained from NSE. The new SANS model is fully consistent with the dynamic NSE results as opposed to the previous descriptions, where it was assumed that all material was contained in the complexes (see Fig. S14† for a comparison of the models). The new SANS model is shown as solid lines in Fig. 1 and describes the experimental data very accurately. See the ESI for details on the fit and Table SI† for the obtained fit parameters.
The values obtained for Dslow do not show a strong dependence on Mw (see Fig. 5). One likely explanation for this situation would be, that due to the soft and flexible nature of the complexes, the diffusion times Dslow measured in NSE belong to segment motions of the complex, rather than diffusion of large complexes. With a diffusion coefficient of 1 Å2 ns−1, this would correspond to segment lengths of about 200 Å, which is in good agreement with the values of the Kuhn length obtained from the new SANS model (see Table SI†). The effective segment lengths are similar for polyelectrolytes of the same nature but different Mw, consequently Dslow is also similar for all different Mw measured even though the static picture shows Mw dependent differences. This explanation could be verified by employing other polyelectrolytes with significantly different persistence lengths which should have a pronounced effect on Dslow and will be done in a consecutive study.
Reverting to the initial question concerning the viscosity of polyelectrolyte/microemulsion complexes, we can now explain the low viscosities observed with the very dynamic situation of these complexes. Even a very high shear modulus G0, arising from an interconnection of the droplets cannot compensate such small relaxation times τ, so no viscosity (η) increase can be observed in this system, as η = τ × G0. This conclusion could only be drawn from the comparison of different measurement techniques, especially the combination of static and dynamic methods. One method alone would have given an incomplete picture. This example shows the importance of always studying the static as well as the dynamic behaviour of a system for a thorough characterisation. The observed highly transient nature of the studied polyelectrolyte/microemulsion complexes (PEMECs) should not only be relevant for this particular system, but also be of importance for the understanding of related systems of oppositely charged colloids and polyelectrolytes, systems as they are abundantly occurring in biology, but also in applied formulations in pharmacy, cosmetics, etc.
Footnote |
† Electronic supplementary information (ESI) available: Viscosity measurements, DLS measurements, examples of NSE and DLS curves, modelling of NSE data and SANS (measurements, fitting models, discussion). See DOI: 10.1039/d0na00336k |
This journal is © The Royal Society of Chemistry 2020 |