A. Gama
Goicochea
*a,
E.
Mayoral
a,
J.
Klapp
ab and
C.
Pastorino
c
aInstituto Nacional de Investigaciones Nucleares, Carretera México-Toluca s/n, La Marquesa Ocoyoacac, Estado de México 52750, Mexico. E-mail: agama@alumni.stanford.edu
bDepartamento de Matemáticas, CINVESTAV del IPN, México D. F. 07360, Mexico
cDepartamento de Física, Centro Atómico Constituyentes, CNEA-CONICET, Av. General Paz 1499, Provincia de Buenos Aires 1650, Argentina
First published on 25th October 2013
We undertake the investigation of sheared polymer chains grafted onto flat surfaces to model liposomes covered with polyethylene glycol brushes as a case study for the mechanisms of efficient drug delivery in biologically relevant situations, for example, as carriers for topical treatments of illnesses in the human vasculature. For these applications, specific rheological properties are required, such as low viscosity at high shear rates, to improve the transport of the liposomes. Therefore, extensive non-equilibrium, coarse-grained dissipative particle dynamics simulations of polymer brushes of various lengths and shear rates are performed to obtain the average viscosity and the friction coefficient of the system as functions of the shear rate and polymerization degree under theta-solvent conditions, and we find that the brushes experience considerable shear thinning at large shear rates. The viscosity (η) and the friction coefficient (μ) are shown to obey the scaling laws η ∼−0.31 and μ ∼
0.69 at high shear rates (
) in a theta solvent, irrespective of the degree of polymerization of brushes. These results confirm recent scaling predictions and reproduce very well trends in measurements of the viscosity at a high shear rate (
) of red blood cells in a liposome containing medium.
Most of the studies have been carried out for polymers under good solvent conditions, with very few exceptions.20,22–24 However, it is known that some biopolymers found in an aqueous environment are in the borderline between good solvent and theta solvent conditions.22,26 An important example is that of polyethylene glycol (PEG) polymer brushes on liposomes, which are used as carriers under flow for efficient drug delivery.27,28 PEG is known to have a decreasing solubility as its molecular weight increases,29 and for molecular weights in the range of 2000 to 10000 g mol−1 in aqueous solution it has been reported to be slightly below its theta temperature, when placed at the human body temperature.30 A detailed understanding of drug delivery mechanisms where PEG is involved, which is of paramount importance for optimal pharmaceutical design, is still far from complete.31 Therefore, our purpose here is to study the viscosity and friction of PEG brushes on surfaces that experience an external shear force under theta solvent conditions, for the first time. To reach scales comparable to those of typical non-equilibrium experiments we have carried out particle-based, mesoscopic molecular dynamics simulations, using the method known as dissipative particle dynamics (DPD),32 which has successfully predicted properties of polymer brushes both in equilibrium and non-equilibrium situations.33
This work reports explicit-solvent DPD simulations of the rheological properties of biopolymer (PEG) brushes grafted onto parallel flat surfaces (liposomes). Including the solvent explicitly in the simulations is crucial to reproduce experiments measuring friction between polymer brushes.18,19 We want to model the behavior of tethered proteins and biopolymers in environments such as vessels, where the separation between plates is essentially constant but we allow for variations in the degree of polymerization, at a fixed grafting density and substrate separation. In particular, we would like to determine the optimal conditions under which PEG brushes on liposomes promote the flow of the latter as mechanisms for efficient drug delivery. We obtain the viscosity and the friction coefficient as functions of the shear rate and degree of polymerization. This is the first report of scaling behavior for the viscosity and the friction coefficient at high shear rates of polymer brushes in a theta solvent using soft DPD interactions. Understanding of these phenomena is useful for the interpretation of several recent experiments in biological and other mesoscopic systems of current academic and industrial interest.34
In Section II, called “Models and methods”, we provide details about the DPD interaction model, the simulation algorithm and the systems studied. In Section III the main results are presented, accompanied by their discussion. The final section is devoted to our major conclusions.
![]() | (1) |
![]() ![]() | (2) |
![]() | (3) |
In eqn (1)–(3), rij = ri − rj represents the relative position vector between particles i and j, êij is the unit vector in the direction of rij, while vij = vi − vj is the relative velocity between particles i and j, with ri and vi being the position and velocity of particle i, respectively. The random variable ξij is generated between 0 and 1 with a Gaussian distribution, zero mean and unit variance; Rc is the cut off distance, beyond which all forces are zero. The DPD beads are all of the same size, with radius Rc, which is set equal to 1. The constants in the dissipation and random forces are not independent and satisfy the relation:35
![]() | (4) |
All our simulations are performed in reduced units (marked with asterisks); distances are reduced with the cutoff radius, Rc, which for a coarse-graining degree equal to three water molecules per DPD particle is equal to Rc = 6.46 Å,38 hence r = r*Rc. The time step δt is reduced with δt = (mR2c/kBT)1/2δt*, where m is the mass of a DPD particle, while energy is reduced with kBT. All other quantities can be reduced through combinations of these relations, for example, the viscosity: η = η*(kBTδt/R3c). Using the mass of 3 water molecules per DPD particle, at room temperature, one obtains δt ≈ (6.3 × 10−12 s)δt* and η ≈ (9.64 cP)η*. All our simulations are performed for (kBT)* = m* = R*c = 1, δt* = 0.03. In addition to the fluid monomeric particles, our system contains soft parallel surfaces as models for biological membranes and polymer chains attached to them, forming brushes. The polymers are built as linear chains of DPD beads joined by freely rotating harmonic springs.40 The spring constant is chosen as κ = κ*(kBT/R2c), where κ* = 100 and the spring's equilibrium position is req = r*eqRc, with r*eq = 0.7 in all cases.41 In our simulations we have fixed the distance between the plates, which are placed perpendicular to the z-direction, at a distance of D* = 7. Periodic boundary conditions are applied on the xy-plane but not in the z-direction, to reinforce the confinement. The equation of motion is solved using the velocity Verlet algorithm, adapted to DPD.42 The parameters of the dissipative and random force intensities are, respectively, γ = 4.5 and σ = 3.0, so that kBT* = 1 (see eqn (4)). The value of the conservative force intensity (see eqn (1)) was set to aij = 78.0 for all cases, namely for the interaction between particles of the same type (solvent–solvent, monomer–monomer) and for particles of different types (solvent–monomer). The choice aij = 78.0 is obtained when one uses a coarse graining degree that groups 3 water molecules in a single DPD bead.38 Since all particle–particle interactions are equal, the polymers are in a theta solvent, as is the case for PEG-covered liposomes in the human circulatory system.30 It has been shown that this model for PEG on colloidal surfaces can successfully predict brush scaling laws,28 as well as adsorption isotherms on metallic oxides and disjoining pressure profiles.39 The soft membranes on which these biopolymers are tethered are modeled as linearly decaying forces that act on the ends of the simulation box (in the z-direction), given by the following expression:39
![]() | (5) |
In eqn (5), Fw(zi) symbolizes the force exerted by the effective wall on particle i, whose position component in the z-direction is zi. The intensity of the force is given by the constant aw, while zC is the cutoff distance also set equal to 1, which defines the reach of this force, i.e., Fw(zi) = 0 for zi > zC. Eqn (5) represents a soft membrane because the maximum repulsion between the beads and the surface is aw when zi = 0; harder walls can be used for other DPD applications.39 It is known that liposomes can be up to about 3 orders of magnitude larger than their brush thickness, and the solvent,2,43 therefore we consider planar walls only. To fix one end of the polymers on these substrates we chose their interaction equal to aw = 70.0 (reduced energy over length DPD units), while the rest of the polymer beads and the solvent interact with the walls with an intensity given by aw = 140.0. With this choice of parameters, the polymer end experiences a less repulsive force toward the walls than the rest of the polymer (and solvent) beads, and is therefore adsorbed on them. The tethered ends of the polymers remain free to move on the xy-plane, of course, as occurs also for biopolymers interacting with membranes. Although this model for membranes is soft, it remains impenetrable to polymer and solvent molecules, and we assume moreover that the drug molecules have already been incorporated into the structure.
To carry out non-equilibrium simulations of these systems, we apply a constant shear velocity to the polymer ends adsorbed on the walls, of equal magnitude but opposite sign for beads on different plates. This is equivalent to moving the plates in opposite directions under the influence of a fixed external force that keeps the plates moving with a constant speed, known as the Couette flow,20–24 see Fig. 1. The units of velocity are given by [v0] = v*0[Rc/δt] ≈ 1.01 × 102v*0 m s−1, using the values of Rc and δt given previously. Some of the solvent monomers are seen to penetrate the PEG brushes and carry them along in their flow to produce lubrication, while there is no interpenetration of the brushes. As pointed out before, our simulations are performed in the canonical ensemble, where particle density and temperature are kept constant; previous studies have shown that completely equivalent results are obtained if one performs the simulations in the grand canonical ensemble, at fixed chemical potential, volume and temperature.18 To keep the simulation results invariant to the particular choice of aij interaction parameters, we fix the global dimensionless density to 3 (ref. 38) in all the simulations reported here; this means that when the degree of polymerization of the PEG brushes (N) is increased, the number of solvent molecules is reduced.
The grafting density is defined as Γ* = Np/A*, where Np is the number of polymer chains tethered on the surface and A* = L*xL*y is its area. Our simulations are performed at grafting densities that lie within the brush regime (Γ*θ ∼ 1/N, for theta-solvent conditions, where N is the degree of polymerization; for polymer chains in a good solvent, Γ*good ∼ 1/N2(0.588), which means that the brush regime is quickly reached for large chains in a good solvent, see ref. 28), where the “stealth” properties of the PEG brushes work best.44 These properties are the increased lifetime in the bloodstream and the kinetic stability of nanoparticles, such as liposomes, due to the protective layer formed by the brushes.31 We obtain the viscosity (η) and the friction coefficient (μ) of the fluid through the relations.20
![]() | (6) |
![]() | (7) |
In eqn (6) and (7) above, 〈Fx()〉 and 〈Fz(
)〉 are the mean forces that the particles on the surfaces experience along the flow direction and perpendicularly to it, respectively; the brackets indicate time average over all particles. Eqn (6) can be understood as the local definition of viscosity, applied to the entire sample. Considering a liquid confined between two planar walls, a linear flow can be generated by moving for example the top wall at a constant velocity. This motion requires that a steady force be applied on the top wall and an equal in magnitude but opposite in direction force on the bottom wall, see Fig. 1. In linear flow, the velocity gradient (shear rate,
) is constant throughout the liquid. The shear rate
is defined as 2v*0/D*, where v*0 is the flow velocity exerted on the grafted monomers, and D* is the surface separation (see Fig. 1). To obtain first a local definition of viscosity, a small cubic volume in the fluid (fluid element) can be considered. This fluid element undergoes strain, being deformed from a cubic to a parallelepiped shape, with increasing deformation in time. The amount of strain is the added lateral displacement of top and bottom planes of the fluid volume divided by its height. This strain increases at a constant rate in time for a linear flow, such that the strain rate, i.e. the strain per unit time, is the velocity gradient in the limit of an infinitesimal fluid element. The strain rate of deformation is therefore the shear rate,
. On the other hand, a force acts on the top and bottom planes of the fluid element, with identical magnitude and opposite direction. The center of mass of the fluid element suffers no acceleration and therefore the total external force on it must vanish. These forces must be proportional to the area of the top and bottom planes of the fluid element. Also, the force-per-unit-area or stress, σ, must be uniform from top to bottom in the fluid element since the fluid is not accelerating. For a linear flow, this is true not only for a fluid element but for the whole fluid: every part of the fluid is under the same stress, which produces everywhere the same velocity gradient (
). The definition of viscosity assumes that the shear rate (
) produced in a fluid element is proportional to the shear stress (σ) exerted on it. The viscosity is defined as the constant of proportionality between them: σ = η
, or η = σ/
. The total shear stress on the sample (σ) can be determined from the average total force on the brush heads divided by the wall area (Fx/A) while the shear rate (
) is extracted from the slope of the linear fit of the average velocity profile. The viscosity is then expressed by eqn (6). As for the friction coefficient, shown in eqn (7), its calculation follows directly from the analogy with the friction coefficient between solid surfaces, namely Fx = μFz, where Fz is the force acting perpendicularly to the surface, and Fx is the force acting parallel to it, which is responsible for the shear stress (σ = Fx/A).25 Our results were obtained from averages of simulations of up to 4 × 103 blocks, with each block composed of 2 × 104 time steps, using the first 2 × 105 time steps for equilibration and the rest for the production phase; when properly dimensionalized this represents a time observation window of 0.12 ms. Typical density profiles of the solvent and polymer monomers are shown in Fig. 2 for two polymerization degrees, where it is clear that the solvent penetrates the PEG brushes and reaches the surfaces, at that given grafting density.
The profiles of the opposing polymer brushes show very little interpenetration and an almost bulk like concentration of solvent monomers at the center of the simulation box (see Fig. 2), which reduces the viscosity of the system, because there are more momentum – carrying solvent molecules in the gap between the surfaces than if there were longer brushes. The total density (solvent plus polymer monomers) of the system was fixed at three; hence at the center of the simulation box one finds mostly solvent particles. The structuring of the PEG brush monomers, represented by the maxima in Fig. 2, coincides with that of the solvent monomers, indicating that even for a relatively large grafting density as that shown in the figure, the solvent is able to reach the membrane. The maxima of both types of profiles (the solvent's and the brushes') occur at the same positions because the interactions between them are the same (recall aij = aii = 78.0 for all i and j), as it is befitting for theta-solvent conditions. Theta solvent conditions are obtained whenever the interactions between polymer chains are the same as the interactions between the solvent and the polymers.2 A strong layering is expected when the chains are relatively short, as for N = 7 in Fig. 2, because in such a case monomers are more easily arranged than when long chains form the brushes.17–21 In Fig. 3 we present the velocity profile for PEG brushes of polymerization degree N = 14 and grafting density Γ* = 0.30 under a shear velocity equal to v*0 = 0.1. Clearly, there appears a linear gradient in the center of the channel formed between liposomes, from which the shear rate can be obtained, through the relation
= 2v*0/D*, as pointed out before. The inset of Fig. 3 shows the temperature profile of the complex fluid in the z-direction, which indicates the brushes are at the temperature fixed by the thermostat (kBT* = 1), although at large shear rates one expects the brush to “heat up” somewhat,20 as a consequence of the increased dissipation rate.
Most simulations were performed for a brush grafting density Γ* = 0.30, except where indicated otherwise, because for this value the average distance between grafted heads on the surface is smaller than their radius of gyration in a theta solvent, for all values of N we studied.28 By doing so one makes sure that the grafted polymer chains are in the brush, rather than the mushroom regime,45 which is important for the situation we are interested in modeling. The polymerization degree was varied in the range N = 1 up to N = 25, so that scaling laws could be extracted from the data. Finally, the shear rate was chosen to vary from = 0 to
= 0.30, except where indicated otherwise, so that we could compare our results with those available in the literature.17–21 Molecular dynamics simulations like ours solve the equation of motion essentially exactly,36 although the use of soft forces in DPD (see, for example, eqn (1)) might limit its applicability to situations where atomistic details are not important, as is the case for the nanotribology studies that are the focus of this work.
Similar results have been obtained in equilibrium DPD simulations46 which proved that finite size effects in the interfacial tension between two model fluids were as small as those found in the viscosity and shown in Fig. 4, when the appropriate ensemble was used. This behavior is to be attributed to the soft, short-range repulsive nature of the DPD forces (see eqn (1)–(3)), although these effects may become non-negligible if a long-range interaction is included, such as the electrostatic one. In the remainder of this work we present results obtained using only the smallest simulation box in Fig. 4.
Let us now proceed to discuss the results obtained for the mean viscosity of a complex fluid made up of PEG brushes and an aqueous solvent, as a function of the applied shear rate on the surfaces, for different degrees of polymerization (N), as seen in Fig. 5. The simulations were performed with a cubic box of lateral size L* = 7 and for Γ* = 0.30. For very small values of (less than 10−4) the calculation of η* requires very long simulations due to poor signal-to-noise ratios,20,24 which makes it very difficult to obtain accurately the zero-shear viscosity, η0, from simulations. However, a good estimate of it can be obtained from the extrapolation of the rheology profiles in Fig. 5 as
approaches zero. For the shorter PEG brushes (N = 1 and N = 3) the fluid displays essentially Newtonian behavior, whereas for the case with N = 7 the fluid shows shear-thinning. One observes also that as N is increased so is the extrapolated η0 value, which is a consequence of the combined effects of an increasing brush thickness and a reduced number of solvent particles (to keep the global density ρ* = 3). Increasing N is equivalent to reducing the separation between the membranes, D*, when the grafting density is kept constant, and yields increasing values of η0.19 Hence, we have opted to vary N rather than D* in this work. The simulations shown in Fig. 5 were carried out for relatively short polymers with the purpose of establishing the minimum polymerization degree at which the fluid started to show shear thinning, which, from Fig. 5, is N = 7. For polymerization degrees smaller than 7, it is difficult for the polymer chains to interact with the solvent monomers enough to reduce the viscosity as the shear rate is increased, hence the fluid remains approximately Newtonian (N = 1 and N = 3 in Fig. 5).
It has been reported that increasing the polymerization degree of PEG brushes on liposomes increases their circulation longevity,29 which would be a desirable aspect to improve drug delivery mechanisms. However, as its degree of polymerization grows, PEG's solubility in aqueous environment is reduced,30 hence there would be a competition between these trends that can be investigated with DPD simulations, varying N. Also, it is important to find out if the viscosity behaves qualitatively the same as a function of shear rate, regardless of the value of N, so that general conclusions are obtained that can be used to interpret a variety of experiments.25 Before presenting our results for the effective mean viscosity of systems where PEG brushes have larger polymerization degrees we comment on the behavior of the mean shear stress on the membranes along the direction of the flow as a function of the shear, 〈Fx()〉. Following Galuschko and coworkers,47 we find a critical shear rate for our systems,
*, as the value of the shear rate when the behavior of 〈Fx(
)〉 changes from linear to sub linear, which depends on N, being
* ∼ 8.6 × 10−3 for N = 14 and
* ∼ 2.9 × 10−3 for N = 25. Finding
* is tantamount to finding the Weissenberg number, We, since We =
/
*, and it signals when shear-thinning behavior starts to set in (We ≥ 1). Fig. 6 shows the dependence of the normalized shear stress, u ≡ 〈Fx(
)〉/〈Fx(
*)〉, on We for fluids with brushes with polymerization degrees equal to N = 14, 20 and 25. The data for the three brushes are seen to collapse reasonably well, especially at low We; this is the so-called linear regime, u ∼ We, (dashed blue line in Fig. 6). For larger values of the Weissenberg number, we find that the power law u ∼ Weκ is obeyed with κ = 0.69 (solid black line in Fig. 6) for our systems under theta-solvent conditions. Galuschko et al.47 predicted κ = 9/13 ≈ 0.69 using scaling arguments for what they call “dry” polymer brushes, which are brushes where the hydrodynamic and excluded volume interactions are screened out, for fairly concentrated solutions and polymer melts, where the corresponding Flory exponent is ν = 0.5.48 Our simulations predict a κ value in excellent agreement with scaling arguments47 because under theta-solvent conditions and at the concentrations we modeled, the brushes are almost indistinguishable from the solvent and act as a melt.
Recent molecular dynamics simulations49 of strongly compressed polyelectrolyte brushes under shear at melt concentrations find that the shear force 〈Fx()〉 that appears in eqn (6) and (7) behaves at high shear rates (
) as 〈Fx(
)〉 ∼
0.69, in agreement with our predictions in Fig. 6, and with the scaling arguments of Galuschko et al.47 Although there are explicit electrostatic interactions in the work of Spirin and Kreer,49 they argue that their brush and solvent system behaves as a melt of neutral polymers at high We because of the immobilization of the counterions. Therefore, the scaling exponent κ = 0.69 we obtained in Fig. 6 is very robust.
The effective mean viscosity data of these systems can also be collapsed as a function of We if one defines s = η/η0,47 for We > 1, as seen in Fig. 7(a). The value of the zero shear viscosity η0 was obtained from extrapolation, using the same procedure as the one described in the discussion of Fig. 5. The data are scattered at low values of the applied shear rate because the signal-to-noise level in those cases is rather poor and one is required to perform very long simulations to get a good estimate of the extrapolated, zero-shear viscosity, as has been reported by others.20,24 All three cases experience considerable shear thinning as the Weissenberg number is increased, and in this regime the data follow the power law s ∼ Weζ, with ζ = −0.31, see the black line in Fig. 7(a). To our knowledge, this is the first time such scaling law has been obtained in a theta solvent.
![]() | ||
Fig. 7 (a) Normalized effective mean viscosity s = η/η0 as a function of the Weissenberg number, We, for fluids with PEG brushes of three different polymerization degrees, N. In all cases, L*x = L*y = D* = 7, and the grafting density is Γ* = 0.30. The line represents the power law, s ∼ Weζ, with ζ = −0.31. (b) Viscosity of a fluid made up of red blood cells (RBC) dispersed in a 40 kDa liposome–dextran medium, as a function of the shear rate.51 The blue line represents the power law η ∼ ![]() |
In ref. 47 the value ζ = −0.43 was obtained for linear Lennard-Jones brushes in a good solvent. A similar value (ζ = −0.42) was found in simulations of bottle-brush polyelectrolytes in a good solvent.19 Hence, the polymer architecture appears to be of secondary importance, compared with the quality of the solvent or the polymer concentration. Galuschko and collaborators47 proposed also the following relation between κ and ζ:
κ − ζ = 1 | (8) |
In Fig. 7(b) we have included the high shear rate values of the viscosity of human red blood cells (RBC) in a liposomal suspension,51 which shows evidently shear-thinning. This non-Newtonian behavior is interpreted to be the result of the dissociation of flocculated liposomes at the higher values of the shear rate, an interpretation that is consistent with our predictions in Fig. 7(a). The blue line in Fig. 7(b) represents the power law we have obtained for PEG-grafted surfaces in a theta solvent, which models accurately the results for RBC at high values of the shear rate, providing confirmation for our predictions. Experiments carried out for gels used as vehicles for drug delivering liposomes,52 such as hydroxyethylcellulose (HEC) and a mixture of HEC and an acrylic acid-based polymer (carbopol 974), yield values of the viscosity at high shear rates very close to that found in Fig. 7(a), namely −ζ = 0.30 ± 0.01 for HEC and −ζ = 0.310 ± 0.009 for the mixture. The fact that the exponent in these experiments is found to be very close to our prediction suggests that it is of paramount importance to take into account the properties of the solvent (theta conditions), as well as those of the brush (PEG), and lends additional support to these calculations.
Adding polymer brushes to sheared surfaces not only modifies the viscosity of the system but it has been shown to reduce its friction coefficient10 also. Since this is a parameter that can be directly compared with experiments, we have obtained the friction coefficient of PEG brushes in an explicitly included theta solvent. Using eqn (7) one calculates the mean friction coefficient between the two sliding surfaces. The data in Fig. 8(a) show that the friction coefficient μ at a constant PEG grafting density increases with the shear rate, in agreement with trends found by others;19,47,53 it increases also with N because the overlap between the brushes also grows. Moreover, the μ values for PEG brushes in a theta solvent are somewhat larger than those obtained in simulations of linear brushes in a good solvent;18 this is of course expected because the brush thickness is reduced when it is immersed in a theta solvent, and the contacts between the polymers grow when going from a good to a theta solvent.2 Experiments by Klein and collaborators54 on the friction of polymer brushes in a theta solvent compared to a good solvent show that the friction coefficient can be larger in the former by up to three orders of magnitude than in the latter. The reason relies on the fact that under theta solvent conditions the brushes can interdigitate and the solvent can penetrate the brushes (see Fig. 2), reducing the freely flowing solvent in the center of the channel, thereby increasing the shear force and the friction. Klein and coworkers54 obtain values for μ which are close to our predictions, shown in Fig. 8(a). When electrostatic interactions are included the values for μ obtained from computer simulations turn out to be slightly larger and closer to 1 for good solvent conditions at high values of the shear rate,19,53 due to the added osmotic pressure of the ions. However, in biological systems like those modeled here, the quality of the solvent plays a key role in providing a mechanism that promotes a lower viscosity environment.
![]() | ||
Fig. 8 (a) Friction coefficient for polymer brushes of various degrees of polymerization (N), obtained using eqn (7). In all cases, L*x = L*y = D* = 7, and the grafting density is Γ* = 0.30. Lines are only guides for the eye. (b) The friction coefficient data shown in (a), normalized by its value at ![]() ![]() ![]() ![]() |
In Fig. 8(b) we show the same data as in Fig. 8(a) on logarithmic scales to emphasize the power law behavior. For the shortest brush (see the inset of Fig. 8(b)), made up of chains with polymerization degree equal to N = 7, there is a linear relation between the friction coefficient and the shear rate, μ ∼ , for the entire range of shear rate values used in the simulations (blue line). To see sublinear behavior here one would have to impose larger shear. For the other polymerization degrees shown (N = 14 to 25) the main panel in Fig. 8(b) shows that the data collapse rather well on a single curve when the friction coefficient (μ) is normalized by its value (μ*) at the start of the shear-thinning regime (
*). The black line is the fit to the power law μ/μ* ∼ Weκ, with κ = 0.69. This is the same exponent as the one obtained from the normalized shear stress data, in Fig. 6, indicating that the force 〈Fz(
)〉 (see eqn (7)) is essentially independent of shear once the brushes have reached a certain length. At the same time, the fact that the exponent is κ = 0.69 at We > 1 in all three cases (N = 14, 20 and 25) shows that 〈Fx(
)〉 (see eqn (7)) is responsible for the behavior of the friction coefficient in the high shear rate regime. The same behavior is found for Lennard-Jones brushes,47 namely the friction coefficient dependence on We scales with the same exponent as that of the shear stress, although in such cases κ = 0.57 because that system is under good-solvent conditions. In ref. 55 we have included a table listing all the data shown in Fig. 8(a) and (b).
It is argued that our simulations are useful for understanding the behavior of biopolymer brushes coating drug-carrying liposomes in an aqueous environment that acts as a theta solvent, under shear. Increasing the thickness of the brush through the degree of polymerization or the shear velocity on the surfaces is shown to raise the friction coefficient, which is a deleterious effect for the transport properties of these liposomes. However, when liposomes are covered by PEG brushes in a theta solvent at high shear rates, their flowing characteristics make them optimal carriers for drug delivery, as the polymer brushes imprint them with efficient injectable characteristics (low viscosity at high shear rates) while at the same time providing them with thermodynamically stable (“stealth”) mechanisms. These simulations have the additional advantage of including hydrodynamic interactions, as well as the solvent explicitly, and being mesoscopic they reproduce length and time scales comparable with those of environments of biological interest. Therefore, we believe our work should be useful in the improved design of drug-carriers, the rheological characterization of sheared brushes, and in establishing general scaling laws for non-equilibrium polymer brushes.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c3sm52486h |
This journal is © The Royal Society of Chemistry 2014 |