On the Electric-Field-Induced Responses of Charged Spherical Colloids in Uncharged Hydrogels and the Anomalous Bulk Viscosity of Polymer-Nanocomposite Melts Mu Wang Master of Engineering Department of Chemical Engineering McGill University Montreal,Quebec June 2008 A thesis submitted to McGill University in partial fulfilment of the requirements of the degree of Master of Engineering c Mu Wang 2008. All rights reserved.
148
Embed
On the Electric-Field-Induced Responses of Charged ...digitool.library.mcgill.ca/thesisfile21951.pdf · Colloids in Uncharged Hydrogels ... understanding and support during various
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
On the Electric-Field-InducedResponses of Charged SphericalColloids in Uncharged Hydrogels
and the Anomalous Bulk Viscosity ofPolymer-Nanocomposite Melts
Mu Wang
Master of Engineering
Department of Chemical Engineering
McGill University
Montreal,Quebec
June 2008
A thesis submitted to McGill University in partial fulfilment of therequirements of the degree of Master of Engineering
Dear Mu Wang The Royal Society of Chemistry (RSC) hereby grants permission for the use of your paper(s) specified below in theprinted and microfilm version of your thesis. You may also make available the PDF version of your paper(s) that theRSC sent to the corresponding author(s) of your paper(s) upon publication of the paper(s) in the following ways: in yourthesis via any website that your university may have for the deposition of theses, via your university’s Intranet or via yourown personal website. We are however unable to grant you permission to include the PDF version of the paper(s) onits own in your institutional repository. The Royal Society of Chemistry is a signatory to the STM Guidelines onPermissions (available on request). Please note that if the material specified below or any part of it appears with credit or acknowledgement to a third partythen you must also secure permission from that third party before reproducing that material. Please ensure that the published article states the following: Reproduced by permission of The Royal Society of Chemistry Regards Gill CockheadContracts & Copyright Executive Gill Cockhead (Mrs), Contracts & Copyright ExecutiveRoyal Society of Chemistry, Thomas Graham HouseScience Park, Milton Road, Cambridge CB4 0WF, UKTel +44 (0) 1223 432134, Fax +44 (0) 1223 423623http://www.rsc.org -----Original Message-----From: [email protected] [mailto:[email protected]] Sent: 30 May 2008 07:16To: CONTRACTS-COPYRIGHT (shared)Subject: Permission Request Form: Mu Wang Name : Mu WangAddress : Department of Chemical EngineeringMcGill UniversityRoom 3060, Wong Building, 3610 University StreetMontreal, Quebec H3A 2B2 Tel : 514-991-6000Fax :Email : [email protected] I am preparing the following work for publication: Article/Chapter Title : Electric-field-induced displacement of charged sphericalcolloids in compressible hydrogelsJournal/Book Title : On the Electric-field-induced Responses of Charged SphericalColloids in Uncharged Hydrogels and the Anomalous Bulk Viscosity ofPolymer-nanocomposite Melts Editor/Author(s) : Mu Wang Publisher : Thesis for Master of Engineering I would very much appreciate your permission to use the following material:
vi
RE: Permission Request Form: Mu Wang imap://exchange.mcgill.ca:993/fetch%3EUID%3E/INBOX%3E1997...
2 of 2 2008-6-9 0:43
Journal/Book Title : Soft MatterEditor/Author(s) : Mu Wang and Reghan J. HillVolume Number : 4Year of Publication : 2008 Description of Material : a research paper titled: Electric-field-induced displacementof charged spherical colloids in compressible hydrogels Page(s) : 1048-1058 Any Additional Comments : I am the first author of the above research paper (Soft Matter, 4, 2008, 1048-1058), andI would like to include a clearly duplicated version (not reprints) of the above paperin my thesis for the degree of Master of Engineering in McGill University, Montreal,Canada. Thank you very much in advance for your kind permission.
DISCLAIMER:
This communication (including any attachments) is intended for the use of the addressee only and may contain confidential, privileged or copyright material. It may not be relied upon or disclosed to any other person without the consent of the RSC. If you have received it in error, please contact usimmediately. Any advice given by the RSC has been carefully formulated but is necessarily based on the information available, and the RSC cannot be held responsible for accuracy or completeness. In this respect, the RSC owes no duty of care and shall not be liable for any resulting damage or loss. The RSC acknowledges that a disclaimer cannot restrict liability at law for personal injury or death arising through a finding of negligence. The RSC does not warrant that its emails or attachments are Virus-free: Please rely on your own screening.
vii
CONTRIBUTION OF AUTHORS
Contents of chapters 3–5 of this thesis are reproduced or adapted from the
papers that have been published or submitted for publication in scientific jour-
nals under the supervision of my research supervisor, Professor R. J. Hill, who
is also a co-author. Chapter 3 discussed the steady electric-field-induced dis-
placement of a colloidal particle in compressible hydrogels, chapter 4 extended
the steady displacement to the dynamic responses and connected the single
particle response to the bulk electroacoustic signals, and chapter 5 theoreti-
cally interpreted the recently discovered bulk viscosity reduction in polymer-
nanocomposite melts.
The research project was initiated by Professor Hill. In chapter 3, using a
displacement construction proposed by Professor Hill, I devised the computa-
tional methodology and boundary layer solution for the steady displacement,
and highlighted the importance of compressibility. In chapter 4, I numerically
calculated dynamic particle displacement and derived a boundary layer ap-
proximation. To overcome numerical difficulties, I also analytically solved a
so-called two-fluid model, and performed a far-field asymptotic analysis on the
governing equations. Furthermore, I have theoretically shown that electric mi-
crorheology and electroacoustic diagnostics can be applied to hydrogel-colloid
composites at low and high frequencies, respectively. In both chapters, I devel-
oped robust computer programs in C and FORTRAN to compute the steady
and dynamic particle responses. In chapter 5, I proposed that the bulk vis-
cosity reduction may arise from Rouse dynamics in entangled melts at small
scales. Professor Hill and I derived a hydrodynamic model to verify this idea,
and I have showed that the layer thickness is qualitatively independent of the
3–1 Poisson’s ratios of selected hydrogels ascertained from experi-ments under undrained and drained conditions . . . . . . . . 13
4–1 Parameters for the results shown in figure 4–1. . . . . . . . . . 76
5–1 Summary of the parameters that characterize the experimentsof Mackay et al. (2003) and Tuteja et al. (2005) with Rg/h &1 and φ = 0.005, and theoretical interpretations (providingfitted values for δ) based on (5.8) with χ = 0 . . . . . . . . . 118
5–2 Best-fit polymer correlation lengths ξ ascertained from experi-ments of Mackay et al. (2003) and Tuteja et al. (2005) and thetheoretical interpretation based on a continuous-layer-profilemodel with ka = 0 . . . . . . . . . . . . . . . . . . . . . . . . 120
analysis can greatly improve the accuracy and stability of the numerical ap-
proach. In this thesis, the solution of the Poisson-Boltzmann equation is based
on the methodology of MPEK, and the perturbed multi-phase electrokinetic
model is solved using TWPBVPL, a general-purpose boundary value problem
solver (Cash & Mazzia 2006). To help improve the computational accuracy,
asymptotic analysis of the governing equations is undertaken.
9
CHAPTER 3Electric-field-induced displacement of charged spherical colloids in
compressible hydrogels1
This chapter concerns the electric-field-induced displacement of a charged
spherical colloid embedded in an uncharged compressible hydrogel. Previ-
ous theoretical calculations for incompressible polymer skeletons predict sub-
nanometer particle displacements within the experimentally accessible param-
eter space (e.g., particle surface charge density, polymer shear modulus, and
electric field strength). Accordingly, the prevailing expectation is that an ex-
perimental test of the theory would be extraordinarily difficult. In this work,
however, we solved the electrokinetic model for compressible polymer skele-
tons with arbitrary Poisson’s ratio. The most striking result, obtained from
numerically exact solutions of the full model and an analytical boundary-layer
approximation, is that polymer compressibility admits particle displacements
that increase linearly with particle size when the radius is greater than the
Debye length. This scaling is qualitatively different than previously obtained
for incompressible skeletons, where the ratio of the particle displacement to
the electric field approaches a particle-size-independent constant. The dis-
placement is also much more sensitive to the hydrodynamic permeability of
the polymer skeleton. Therefore, when compressible hydrogels are deformed
at frequencies below their reciprocal draining time, our theory identifies the
parameter space where displacements could be registered using optical mi-
croscopy. In turn, this will help to establish a quantitative connection between
the electric-field-induced particle displacement and physicochemical character-
istics of the particle-polymer interface.
1 Reproduced by permission of The Royal Society of Chemistry
10
3.1 Introduction
Hydrogels are polymer networks that have found widespread use in tis-
sue engineering (Barndl et al. 2007), drug delivery (Qiu & Park 2001), and
molecular separations, e.g., gel-electrophoresis, isoelectric focusing, and iso-
tachophoresis (Westermeier 2005). The networks are often synthesized from
polymers such as poly(methyl methacrylate) (PMMA), poly(vinyl alcohol)
(PVA), and polyacrylamide (PA); as well as from macromolecules of biologi-
cal origin, such as collagen and agar.
Recently, several novel applications of hydrogel nano-composites have
been demonstrated where organic and inorganic nanoparticles are immobi-
lized in otherwise conventional hydrogel matrices. For example, wavelength-
selective light-induced swelling from gold and gold-coated silica nanoparticle
inclusions makes these intriguing materials useful as light-activated microflu-
idic valves (Sershen et al. 2005). Nanoparticles have also been introduced
into soft biological tissues to increase the sensitivity of ultrasound imaging for
early tumor detection (Liu et al. 2006; Dayton & Ferrara 2002), and to adsorb
optical energy for treating certain cancers (Loo et al. 2005). Larger colloidal
inclusions have been used to probe the local and bulk viscoelastic response
of polymer solutions and gels (Schnurr et al. 1997; MacKintosh & Schmidt
1999). Finally, in an effort to control the otherwise diffusion-limited transfer
of uncharged molecules across membranes in biosensing, silica nanoparticles
have been embedded in uncharged hydrogel gels to produce electroosmotic
flow (Matos et al. 2006).
As a first step toward understanding the coupling of electroosmotic flow
and polymer deformation in hydrogel composites, Hill & Ostoja-Starzewski
(2008) calculated the electric-field-induced displacement of particles embed-
ded in incompressible polymer skeletons. Their work demonstrates that a
simple balance between the bare Coulomb force and an elastic restoring force
on the particles prevails only when the particle radius a is smaller than the
11
Debye screening length κ−1. Otherwise, the theory quantifies how electroos-
motic flow—in the diffuse layer of countercharge that envelops each particle—
deforms the polymer skeleton and, therefore, influences the particle displace-
ment. For incompressible skeletons, the ratio of the particle displacement to
the electric field strength bears a striking resemblance to the electrophoretic
mobility (O’Brien & White 1978) at all values of κa. Accordingly, the electric-
field-induced particle displacement reflects the size and charge of the inclu-
sions, the viscosity and concentration of the electrolyte, and the shear modu-
lus and hydrodynamic permeability of the polymer skeleton. In principle, the
electric-field-induced displacement is an appealing diagnostic for probing the
physicochemical characteristics of the particle-polymer interface, in a similar
way that electrophoresis is routinely used to ascertain the surface charge of
colloidal particles dispersed in Newtonian electrolytes. However, in the exper-
imentally accessible parameter space, the particle displacements predicted by
Hill & Ostoja-Starzewski (2008) are extraordinarily small, making it difficult
to envision practical diagnostic applications.
In the absence of electroosmotic flow, the elastic restoring force of the
polymer skeleton varies by up to 25 percent over the experimentally accessible
range of Poisson’s ratios for hydrogels (Schnurr et al. 1997). Therefore, when
κa 1, the electric-field-induced displacement of a particle in a compressible
polymer network increases by only 25 percent over the value for an incom-
pressible skeleton. It is therefore unlikely that finite compressibility would
significantly influence the sensitivity of an experiment to test Hill and Ostoja-
Starzewski’s theory when κa 1. However, the situation when κa & 1, which
is generally achieved for particles larger than about one micron, is not as
straightforward to interpret. In this chapter, we show that the electric-field-
induced particle displacement of sufficiently large particles in compressible
matrices is qualitatively different than in incompressible skeletons. Rather
than tending to a size-independent value, the particle displacement increases
linearly with particle size when κa 1. For particles whose radius is greater
12
Table 3–1: Poisson’s ratios of selected hydrogels ascertained from experimentsunder undrained and drained conditions. Note that HEMA-AA represents2-hydroxyethyl methacrylate (HEMA) acrylic acid (AA) co-monomer gel.
Class I (undrained)Hydrogel νHEMA-AA (Johnson et al. 2004a,b) 0.42 – 0.45polyacrylamide (Boudou et al. 2006) 0.487± 0.013polyacrylamide (Takigawa et al. 1996) 0.457± 0.011polyacrylamide (Engler et al. 2004) 0.4 – 0.45poly(vinyl alcohol) (Urayama et al. 1993) 0.433
Class II (drained)Hydrogel νagarose gel (Freeman et al. 1994) 0.15± 0.09resorcinol-formaldehyde (Gross et al. 1997) 0.124 – 0.233silica gel (Scherer 1992) 0.216 – 0.244polyacrylamide (Li et al. 1993) 0.24 – 0.36polyacrylamide (Geissler & Hecht 1980, 1981) 0 – 0.25
than about one micron, our theory predicts displacements of tens to hundreds
of nanometers with modest electric-field strengths and electrolyte concentra-
tions.
While the thermodynamically admissible range for Poisson’s ratio is from
−1 to 0.5, Geissler & Hecht (1980, 1981) established Poisson’s ratio’s of 0 and
0.25 for polymer skeletons in poor and good solvents, respectively. This range
is corroborated to some extent by experiments, but only after two classes of
experiments are identified. As summarized in table 3–1, class I experiments,
which often involve measurements of strain immediately after the initial de-
formation, or with boundary conditions that prevent draining, yield Poisson’s
ratios greater than about 0.4. These experiments reflect the incompressibil-
ity of the solvent. In contrast, the Poisson’s ratios from class II experiments,
where the polymer is permitted to drain, are often in the range 0–0.25 pre-
dicted by Geissler & Hecht (1980, 1981).
The electrokinetic (multi-phase) model of Hill & Ostoja-Starzewski (2008)
generalizes a bi-phasic model (polymer and solvent) where the solvent is hy-
drodynamically coupled to a linearly elastic polymer skeleton. The bi-phasic
13
model can be traced to early works of Biot (1941) and Frenkel (1944) per-
taining, respectively, to the consolidation and seismoelectric behavior of soils,
as well as the propagation of sound waves in geological exploration (Biot
1956a,b). More recently, the bi-phasic model—also termed a two-fluid model—
has been adopted in the relatively new field of microrheology to interpret the
dynamics of entangled polymer solutions and gels (Brochard & de Gennes
figure 2) with ν = 0.5. As κa increases, the other independent dimensionless
parameter κ` = (`/a)κa = 0.01κa increases accordingly. As is customary,
we have plotted the results with constant ζ-potential, so the surface charge
density σ increases with κa according to the non-linear Poisson-Boltzmann
equation (Overbeek 1943). For example, when κa 1, σ ≈ εoεsa−1ζ, and
when κa 1, σ ≈ εoεsκζ. The ostensible increase in the displacement with
28
−15 −10 −5 0 5 10 15−15
−10
−5
0
5
10
15
κx
κy
(a) ψ = ψ0 − rE cos θ + ψ(r)E cos θ
−15 −10 −5 0 5 10 15−15
−10
−5
0
5
10
15
κx
κy
(b) −rE cos θ + ψ(r)E cos θ
−15 −10 −5 0 5 10 15−15
−10
−5
0
5
10
15
κx
κy
(c) ψ(r)E cos θ
Figure 3–1: Streamlines (black lines, with flow from left to right), polymerand particle displacement (vectors, not to scale), and electrostatic potentialisocontours (colored lines, with potential increasing from blue to red) on aplane that contains the polar axis (all fields are axisymmetric): a = 500 nm,κ−1 = 100 nm; ζ ≈ −100 mV; ` = 5 nm, E = 1 kPa, ν = 0.2; E ≈ 25 V cm−1.
29
increasing κa can therefore be attributed, in part, to the accompanying in-
crease in surface charge.
When κa = 0.1 (results not shown), the numerical calculations coincide
with (3.2), which, recall, reflects a direct balance between the bare Coulomb
force and elastic restoring force of the polymer skeleton. Under these con-
ditions, the particle displacement is extraordinarily small. When κa & 1,
however, the displacement is considerably larger. Comparing the numeri-
cal (solid lines) and boundary-layer theory (dashed lines) demonstrates that
(3.48) provides satisfactory predictions of the displacement when κa & 10 and
|ζ| < kT/e. It is noteworthy that the boundary-layer theory, which, recall,
also rests on the Debye-Huckel approximation (|ζ| < kT/e), is also reliable
when |ζ| > kT/e if κa is sufficiently large. Therefore, (3.48) furnishes accu-
rate predictions of the displacement for sufficiently large particles embedded
in any (uncharged) hydrogel at any reasonable ionic strength.
As identified in the introduction, negative Poisson’s ratios are not relevant
for hydrogel skeletons. Nevertheless, for completeness, it is interesting to note
that a negative Poisson’s ratio indicates an isotropic continuum that, under
homogeneous axial strain, adopts an equally signed transverse strain. More
specifically, the axial and transverse strain are equal when ν = −1, and,
therefore, the material changes density without changing shape. In general,
such materials can only support an isotropic state of strain, which is clearly
evident when writing the strain in terms of the stress (Landau & Lifshitz
1986). Therefore, because the strain tensor must be linear in the electric
field or displacement vector, the displacement must be zero when ν = −1, as
confirmed by the numerics and boundary-layer theory.
In the range of Poisson’s ratios appropriate for hydrogel skeletons (0 ≤
ν < 0.5), the displacement achieves a maximum that depends on κa. The
maximum displacement is achieved at a Poisson’s ratio that approaches zero
with increasing κa, and the sensitivity of the displacement to Poisson’s ratio—
as measured qualitatively by the ratio of the maximum displacement to the
30
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4
Poisson’s ratio ν
-0.018
-0.016
-0.014
-0.012
-0.010
-0.008
-0.006
-0.004
-0.002
0
Z/E
, [nm
/(V
/cm
)]
1
2
4
6
(a) κa = 1
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4
Poisson’s ratio ν
-0.05
-0.04
-0.03
-0.02
-0.01
0
Z/E
, [nm
/(V
/cm
)]
1
2
4
6
46
(b) κa = 10
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4
Poisson’s ratio ν
-0.28
-0.24
-0.20
-0.16
-0.12
-0.08
-0.04
0
Z/E
, [nm
/(V
/cm
)]
1
2
4
6
(c) κa = 100
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4
Poisson’s ratio ν
-0.6
-0.5
-0.4
-0.3
-0.2
-0.1
0
Z/E
, [nm
/(V
/cm
)]
1
2
4
6
(d) κa = 1000
Figure 3–2: The ratio Z/E as a function of Poisson’s ratio ν for various scaledreciprocal double-layer thickness κa = 1, 10, 100 and 1000, and scaled ζ-potentials −ζe/(kT ) = 1, 2, 4 and 6: NaCl at T = 298 K; a/` = 100; andE = 1 kPa. Solid lines are numerically exact calculations (a = 500 nm and` = 5 nm), and dashed lines are the boundary-layer theory in (3.48). Notethat the displacement is inversely proportional to E .
31
displacement when ν = 0.5—increases significantly with κa. This important
result, which is also captured by the boundary-layer theory in (3.47), produces
significantly larger particle displacements in compressible skeletons than in in-
compressible ones when κa is large. In other words, finite compressibility
permits large particles to undergo relatively large electric-field-induced dis-
placements. The displacement is plotted in figure 3–3 as a function of κa for
various ζ-potentials with ν = 0.2. As highlighted in the introduction, this
Poisson’s ratio is representative of the values ascertained by several indepen-
dent experiments reported in the literature involving hydrogels under drained
conditions. Again, note that the particle radius a = 500 nm and Brinkman
screening length ` = 5 nm are fixed, so κ` = (`/a)κa = 0.01κa. This way
of plotting the results clearly identifies the ranges of κa and ζ-potential over
which the analytical theories for small (dash-dotted lines) and large (dashed
lines) κa are accurate. Furthermore, in the parameter space where the particle
displacements are large and, therefore, most easily measured, the boundary-
layer theory in (3.47) is reliable.
To draw comparisons between the displacement and electrophoretic mo-
bility, which Hill & Ostoja-Starzewski (2008) demonstrated are very closely
connected when ν = 0.5, figure 3–4 presents the scaled particle displacement
−(Z/E)Ee/(2εoεskT ) as a function of the scaled ζ-potential for various re-
ciprocal double-layer thickness κa with `/a = 0.01 and ν = 0.2. In addi-
tion to testing the boundary-layer approximation (dashed lines, right panel),
this figure clearly identifies maximums in the displacement due to polariza-
tion and relaxation. These processes are well known from their influences on
the electrophoretic mobility (O’Brien & White 1978) and incremental pore
mobility (Hill 2006d). Note also that, due to finite polymer compressibility
(ν = 0.2), the scaled displacements in figure 3–4 (right panel) are significantly
larger than the corresponding scaled electrophoretic mobility (Hill & Ostoja-
Starzewski 2008, figure 5). Even the general shape of the collection of curves is
different from the incompressible limit. In fact, the overall trends bear a much
32
10-2
10-1
100
101
κa
-0.030
-0.025
-0.020
-0.015
-0.010
-0.005
0
Z/E
, [nm
/(V
/cm
)]
1
2
4
6
8
(a)
100
101
102
103
κa
-0.6
-0.5
-0.4
-0.3
-0.2
-0.1
0
Z/E
, [nm
/(V
/cm
)]
1
2
4
6
8
(b)
Figure 3–3: The ratio Z/E as a function of the scaled reciprocal double layer-layer thickness κa for various scaled ζ-potentials −ζe/(kT ) = 1, 2, 4, 6 and8: NaCl at T = 298 K; a/` = 100; E = 1 kPa; and ν = 0.2. Solid lines arenumerically exact calculations (a = 500 nm and ` = 5 nm), and the dashedand dash-dotted lines are, respectively, analytical theories for large and smallκa. Note that the displacement is inversely proportional to E .
33
closer resemblance to the incremental pore mobility (Hill 2006d), thereby pro-
viding a valuable clue toward understanding how the compressibility of the
polymer skeleton influences the particle displacement.
Note that the shape of the electrophoretic mobility versus ζ-potential and
κa relationship (and hence that of the scaled displacement when ν = 0.5) prin-
cipally reflects the variation of the electric-field-induced force on a fixed parti-
cle, since the balancing hydrodynamic drag is, to a first approximation (i.e., in
the absence of electroviscous effects), equal to the Stokes drag force −6πηaU ,
where U is the particle velocity. By direct analogy, the displacement of a par-
ticle embedded in an incompressible matrix reflects the electric-field-induced
force on a fixed particle, since the elastic restoring force is exactly the value
given by (3.1) with ν = 0.5. As demonstrated by Hill & Ostoja-Starzewski
(2008), the hydrodynamic coupling of the fluid and an incompressible polymer
skeleton leads to the same state of stress on the particle as in electrophoresis.
When the polymer is compressible, however, the net force due to the
electric field when the particle is fixed at the origin is evidently dominated by
the hydrodynamic and electrical contributions
f e,E + f d,E = 4π(η/`2)CEE → 6πεoεsζκa2
(κ`+ 1)E as κa→∞, (3.50)
not those arising from the electroosmotic-flow-induced distortion of the poly-
mer skeleton. Accordingly, Hill’s interpretation of (3.50), which does not in-
volve polymer deformation (Hill 2006b), also provides an appealing interpre-
tation of the particle displacement in compressible polymer skeletons. More
specifically, Hill (2006b) showed that the force represented by (3.50) exceeds
the bare electrical force when κa 1 and κ` 1. This is due to an ad-
verse pressure gradient (increasing pressure in the direction of electroosmotic
flow) that must develop to sustain a Darcy (pressure driven) flow in the far
field. Since this force increases with the square of the particle radius a, and
the elastic restoring force is linear in a as shown in (3.1), it follows that the
34
0 2 4 6 8 10
-ζe/(kT)
0
1
2
3
4
5
6
7
scal
ed d
ispl
acem
ent
0.01 0.10.2
0.4
1
2
(a)
0 2 4 6 8 10
-ζe/(kT)
0
20
40
60
80
100
scal
ed d
ispl
acem
ent
3
610
15
30
50
100
1000
(b)
Figure 3–4: The scaled displacement −(Z/E)Ee/(2εoεskT ) as a function of thescaled ζ-potential −ζe/(kT ) for various scaled reciprocal double-layer thick-ness κa = 0.01, 0.1, 0.2, 0.4, 1 and 2 (left panel); and κa = 3, 6, 10, 15, 30,50, 100 and 1000 (right panel): KCl at T = 298 K; a/` = 100; and ν = 0.2.Solid lines are numerically exact calculations (a = 500 nm and ` = 5 nm) andthe dashed lines (κa ≥ 10) are the boundary-layer theory in (3.47).
particle displacement should increase linearly with a when κa 1. This is
indeed verified by our numerics and boundary-layer theory.
Let us briefly address the influence of the polymer hydrodynamic perme-
figure 4) showed it has a very weak influence when ν = 0.5. The particle
displacement is plotted in figure 3–5 as a function of the scaled Brinkman
screening length κ` with ν = 0.2. In striking contrast to the incompress-
ible limit, but consistent with Hill’s interpretation of the force represented by
(3.50), the displacement undergoes a significant transition from a high plateau
when κ` is small (low permeability) to a much lower plateau when κ` is large
(high permeability). Again, when the permeability is low, the electrical force is
accentuated by an adverse pressure gradient. However, when the permeability
35
10-3
10-2
10-1
100
101
102
103
104
κl
10-3
10-2
10-1
100
-Z/E
, [nm
/(V
/cm
)]
10
100
1000
1
0.1 & 0.01
Figure 3–5: The ratio Z/E as a function of scaled Brinkman screening lengthκ` for various scaled reciprocal double-layer thicknesses κa = 0.01, 0.1, 1, 10,100 and 1000: NaCl at T = 298 K; ζe/(kT ) = −1; E = 1 kPa; and ν = 0.2.Solid lines are numerically exact calculations (a = 500 nm) and the dashedlines are the boundary-layer theory in (3.47). Note that the displacement isinversely proportional to E .
is high, the weak coupling of the polymer and fluid leads to a balance be-
tween the electrical and particle-displacement-induced elastic restoring force
as shown in (3.1). It should also be noted that the boundary-layer theory
in (3.47) is valid only when κa 1 and a/` 1, so, as expected, the ap-
proximation breaks down in figure 3–5 when κ` > κa. The small discrepancy
between the boundary-layer and numerical solutions when κa 1 and κ` 1
reflects a subtle but minor shortcoming of the boundary-layer approximation.
Specifically, the leading-order boundary-layer approximation for κa 1 and
` a in the inner region does not include a small but finite tangential pres-
sure gradient. A similar limitation was identified by Hill & Saville (2005)
when comparing numerically exact calculations of the electrophoretic mobility
of soft colloidal spheres with analytical theories in the literature for large κa.
36
3.5 Summary
We generalized the electrokinetic model of Hill & Ostoja-Starzewski (2008)
to calculate the electric-field-induced displacement of a charged, spherical col-
loid embedded in an electrolyte-saturated compressible polymer skeleton. The
fluid velocity and polymer displacement fields were calculated to linear order
in perturbations from an equilibrium state governed by the non-linear Poisson-
Boltzmann equation. Using linear superposition, we expressed the particle dis-
placement in terms of asymptotic coefficients that describe the far-field decays
of the fluid pressure and polymer displacement. Because the polymer skeleton
is compressible, two asymptotic coefficients are necessary to correctly quan-
tify the polymer distortion and, hence, to compute the electric-field-induced
particle displacement. In addition to numerically exact solutions of the full
model, we derived an analytical boundary-layer solution for the limit in which
the inclusion radius is larger than the Debye and Brinkman screening lengths.
Our theory reveals an electric-field-induced particle displacement that is
a sensitive and, in general, complicated function of the Poisson ratio and
hydrodynamic permeability of the polymer skeleton, the size and charge of
the inclusion, and the concentration of the electrolyte. However, the parti-
cle displacement remains inversely proportional to Young’s modulus (or shear
modulus) of the polymer. More importantly, polymer compressibility yields
a particle displacement that has a qualitatively different dependence on the
particle size than for incompressible polymer skeletons. Specifically, when
κa 1, the electric-field-induced displacement increases linearly with the
particle size, rather than tending to a size-independent value. Fortunately,
in the parameter space where experimentally measurable particle displace-
ments are expected, our boundary-layer approximation furnishes a reliable
and convenient alternative to numerical solutions of the full model. To our
best knowledge, experiments have not been reported in the literature, so we
hope our theory will stimulate future experimental work in this area.
37
CHAPTER 4Dynamic electric-field-induced response of charged spherical
colloids in uncharged hydrogels
Embedding colloidal particles in polymeric hydrogels often endows the
polymer skeleton with appealing characteristics for microfluidics and biosens-
ing applications. This chapter provides a rigorous foundation for interpreting
active electrical microrheology and electroacoustic experiments on such mate-
rials. In addition to viscoelastic properties of the composites, these techniques
sense physicochemical characteristics of the particle-polymer interface. We
extended the steady multi-phase electrokinetic model in the previous chapter
to calculate the dynamic response of charged spherical colloids embedded in
uncharged hydrogels when subjected to harmonically oscillating electric fields.
The frequency response depends on the particle size and charge, ionic strength
of the electrolyte, and elastic and hydrodynamic characteristics of the polymer
skeleton. Our calculations capture the transition from quasi-steady compress-
ible to quasi-steady incompressible dynamics as the frequency passes through
the reciprocal draining time of the gel. At higher frequencies, the dynamics
are dominated by hydrodynamic viscous and inertial forces, with the response
eventually becoming equal to the classical dynamic electrophoretic mobility
measured using modern electroacoustic instruments. We establish the connec-
tion between the electroacoustic signal for hydrogel composites and the single-
particle dynamic electrophoretic mobility. Finally, we provide an approximate
analytical theory that captures the elasticity of the hydrogel at the ultrasonic
frequencies used in commercially available electroacoustic instruments. This
agrees very well with our numerics over a wide range of the experimentally
accessible parameter space.
38
4.1 Introduction
Hydrogels are an important class of soft matter that have gained widespread
application in drug delivery (Qiu & Park 2001; Lin & Netters 2006; Peppas
et al. 2000), tissue engineering (Khademhosseini & Langer 2007; Barndl et al.
ωd = 2π(1 − 2ν)−1(E/η)(`/a)2, and ωv = 2πµ/η. The first two of these were
adopted by Levine & Lubensky (2001), and the other two are, respectively,
the reciprocal draining time and reciprocal viscous time introduced in § 4.1.
The results with ` = 1 and 100 nm are representative of tightly and loosely
coupled fluid and polymer. Other parameters are summarized in table 4–1.
Recall, the first and second Lame constants (µ and λ) of the hydrogel skeleton
are taken to be real constants.
The transition from quasi-steady compressible to incompressible dynam-
ics is evident from the plateau seen in figure 4–1(a) at intermediate frequencies
76
10-2
100
102
104
106
108
1010
ω, (rad/s)
-50
0
50
100
150
200
250
300
350R
e[α(
ω)]
, (m
/N)
103
104
105
106
107
108
10910
-3
10-2
10-1
100
101
102
103
|Re[
α(ω
)]|,
(m/N
)ω
B ωd
ωv ω∗
(a) ` = 1 nm, Re(α)
10-2
100
102
104
106
108
1010
ω, (rad/s)
0
40
80
120
160
Im[α
(ω)]
, (m
/N)
103
104
105
106
107
108
10910
-3
10-2
10-1
100
101
102
103
|Im[α
(ω)]
|, (m
/N)
ωB ω
d
ωv
ω∗
(b) ` = 1 nm, Im(α)
10-2
100
102
104
106
108
1010
ω, (rad/s)
-50
0
50
100
150
200
250
300
350
Re[
α(ω
)], (
m/N
)
ωB
ωd
ωv ω∗
(c) ` = 100 nm, Re(α)
10-2
100
102
104
106
108
1010
ω, (rad/s)
0
40
80
120
160
Im[α
(ω)]
, (m
/N)
ωB
ωd
ωv
ω∗
(d) ` = 100 nm, Im(α)
Figure 4–1: Response function α(ω) as a function of angular frequency ωfor different Brinkman screening lengths: ` = 1 nm [panels (a) and (b)] and` = 100 nm [panels (c) and (d)]. Other parameters are listed in table 4–1.Solid lines are exact solutions of the two-fluid model (4.20), dashed lines arethe approximation of Levine & Lubensky (2001), and the dash-dotted linesare the GSER. Several characteristic frequencies are identified (see text fordetails): ωv ≈ 2.9×106 rad s−1 and ω∗ ≈ 8.8×106 rad s−1 for all panels; ωB ≈12 rad s−1, ωd ≈ 47 rad s−1 [panels (a) and (b)]; and ωB ≈ 1.23× 105 rad s−1,ωd ≈ 4.70× 105 rad s−1 [panels (c) and (d)].
77
with ` = 1 nm. When ` = 100 nm, however, the draining time τd is compa-
rable to the viscous time τv, so the low-frequency plateau (compressible elas-
tic regime) in figure 4–1(c) transits to the high-frequency viscous dominated
regime without an intermediate (incompressible elastic) plateau.
The real parts of the quasi-steady compressible and quasi-steady incom-
pressible elastic plateaus differ by at most 25%, as predicted by Schnurr et al.
(1997). This relatively small change is often used to justify neglecting com-
pressibility when interpreting optical and magnetic tweezers microrheology
experiments (Schnurr et al. 1997; Ziemann et al. 1994). However, if the exter-
nal force is accompanied by electroosmotic flow, the effect of compressibility
is much more significant when κa & 1 (Wang & Hill 2008). This is explored
when we address the electric-field-induced response.
Levine and Lubensky’s approximation is valid when fluid inertia can be
neglected. Accordingly, it agrees well with the exact solution when ω ω∗.
As seen in the insets of figures 4–1(a) and 4–1(b), which have logarithmic
axes, Levine and Lubensky’s theory yields increasingly large relative errors
at higher frequencies. However, the absolute displacement is practically zero
at such high frequencies, so the errors are of minor concern for magnetic and
optical microrheology, but have important consequences for electroacoustics.
Earlier studies suggest the GSER is valid at frequencies between ωB and
ω∗ (Levine & Lubensky 2001, 2000). However, our exact results in figure 4–
1 show that the GSER is a good approximation only when the quasi-steady
incompressible elastic plateau is present. At higher frequencies, the GSER de-
viates from the exact solution at frequencies well below ω∗, because it does not
account for fluid or particle inertia. Consequently, our calculations reveal that
the frequency range of validity for the GSER is considerably narrower than
previously expected. Also, the characteristic frequencies ωB and ω∗ intro-
duced by Levine & Lubensky (2001) are similar here to ωd and ωv adopted in
this work; they are practically equivalent for mapping out the experimentally
accessible parameter space.
78
Our exact solution of the two-fluid model is essential for the following
numerical and analytical solutions of the multi-phase electrokinetic model. As
described in § 4.4, the two-fluid model provides far-field boundary conditions
for the fluid velocity and polymer displacement fields, and it is also the ba-
sis of our analytical approximation for the high-frequencies encountered in
electroacoustics.
4.7.2 Numerical solution of the multi-phase electrokinetic model
The model presented in § 4.3 is solved numerically by adopting κ−1, u∗ =
εsεo(kT/e)2/(ηa) and ηu∗/µ = εsεo(kT/e)
2/(µa) as the characteristic scales
for length, fluid velocity and polymer displacement, respectively; and similarly
to (4.66), the dimensional response Z/E is obtained from the dimensionless
Separate programs were written to compute the response for compressible and
incompressible hydrogels. Asymptotic coefficients are extracted from the far-
field decay of the perturbations, and the dynamic response Z/E is obtained
from the superposition in (4.66) or (4.128). Note that the asymptotic analysis
detailed in § 4.4 permits Z/E to be calculated over an extraordinarily wide
range of frequencies, from as low as 0.01 Hz to higher than 1 GHz. An algo-
rithmic description of the computational methodologies and external libraries
used in our programs is provided in Appendix 4.B.
A representative spectrum of Z/E for a colloidal particle in a compressible
hydrogel is presented in figure 4–2. The computations are validated, in part, by
the steady boundary-layer results for compressible and incompressible hydro-
gels shown in (4.1) (Hill & Ostoja-Starzewski 2008; Wang & Hill 2008). Note
that the response (real part) undergoes a distinct transition from quasi-steady
compressible to quasi-steady incompressible elastic plateaus as the frequency
passes through the reciprocal draining time τ−1d ≈ 150 Hz. The transition is
79
slow, however, spanning several frequency decades, with the steady compress-
ible asymptote realized at extremely low frequencies ∼ 10−2 Hz. Note that
the small discrepancy between the steady (horizontal line) and low-frequency
dynamic asymptote in figure 4–2 is due to small errors in the boundary-layer
approximation (4.1), which, in this example, is about 2% smaller than the
numerically exact value (Wang & Hill 2008).
The quasi-steady compressible and incompressible elastic responses differ
by an order of magnitude. At higher frequencies, |Z/E| decays to zero be-
cause the viscous and inertial stresses dominate the particle response. Note
that in active microrheology, most experimentally accessible frequencies are
in the transition from the quasi-steady compressible to incompressible elastic
regimes, so these dynamic calculations are essential for correctly interpreting
such experiments. For the hydrogel-colloid composite in figure 4–2, an applied
electric field E = 20 V/cm with frequency ω/(2π) = 1 Hz induces a parti-
cle displacement with amplitude Z ≈ 4 nm, which could be resolved using
back-focal-plane interferometry (Allersma et al. 1998). The sub-nanometer
displacements at higher frequencies (& 10 kHz) are clearly too small to detect
using any direct measurement of particle displacement. Instead, the parti-
cle velocity −iωZ should be measured, as routinely undertaken at ultrasonic
frequencies in electroacoustic experiments.
The high-frequency regime is examined in figure 4–3, where the absolute
values of the real and imaginary parts of Z/E for particles in a hydrogel (solid
and dashed lines) are compared with their counterparts (Z/E)∗ = −µd∗/(iω)
for the same particles dispersed in the electrolyte without polymer. Accord-
ingly, quantities with superscripts “∗” are from numerically exact solutions of
the standard electrokinetic model (Mangelsdorf & White 1992), as calculated
by the MPEK software package (Hill et al. 2003a). Important characteris-
tic frequencies in § 4.1, including the reciprocal draining time τ−1d , reciprocal
viscous time τ−1v , and reciprocal inertia time τ−1
f are identified. Again, the
80
10-2
100
102
104
106
108
frequency, (Hz)
0
0.05
0.1
0.15
0.2
0.25
0.3
-Z/E
, [nm
/(V
/cm
)]
10-2
100
102
104
106
10810
-6
10-5
10-4
10-3
10-2
10-1
100
| Z/E
|, [
nm/(
V/c
m)]
Figure 4–2: Representative frequency spectrum of Z/E for a charged colloidalsphere embedded in an uncharged, compressible electrolyte-saturated hydro-gel: NaCl at T = 298 K; a = 500 nm; κa = 500; −ζe/(kT ) = 3; ` = 5 nm;ρp = 1050 kg m−3, ν = 0.2; and E = 800 Pa. Solid and dashed lines are thereal and imaginary parts of Z/E from numerically exact solutions of the fullmulti-phase electrokinetic model; and dash-dotted lines are from (4.1) withν = 0.2 and 0.5.
transition from quasi-steady compressible to quasi-steady incompressible elas-
tic dynamics is clearly evident as the frequency passes through τ−1d ≈ 150 Hz.
More importantly, this figure highlights the transition from quasi-steady elas-
tic to viscous dynamics as the frequency passes through τ−1v ≈ 370 kHz. At
frequencies beyond τ−1f ≈ 3.6 MHz, both the real and imaginary parts of Z/E
equal their (Z/E)∗ counterparts. Therefore, it is only at these ∼MHz frequen-
cies that the electroacoustic response of the hydrogel composite is the same as
for its respective colloidal dispersion (Hunter 1998; O’Brien 1988). Moreover,
at lower frequencies, the electroacoustic response probes the shear modulus of
the polymer skeleton and the size and charge of the inclusions.
Having identified several qualitative features of a typical frequency spec-
trum, let us explore the influences of various parameters. First, figure 4–4
shows how Poisson’s ratio, increasing from 0 to 0.5 with fixed shear modulus
µ ≈ 0.333 kPa, affects the response. Recall, the spectrum for the incom-
pressible hydrogel (ν = 0.5) was calculated using an independently developed
81
10-2
100
102
104
106
108
frequency, (Hz)
10-6
10-5
10-4
10-3
10-2
10-1
100
| Z/E
|, [
nm/(
V/c
m)]
τv
-1
τd
-1
τf
-1
Figure 4–3: Comparison of Z/E for a compressible hydrogel with the response(Z/E)∗ = −µd∗/(iω) for a classical Newtonian dispersion. Parameters are thesame as in figure 4–2. Solid and dashed lines are the real and imaginary partsof Z/E from numerically exact solutions of the full multi-phase electrokineticmodel. Dash-dotted and dash-double-dotted lines are the real and imaginaryparts of (Z/E)∗ = −µd∗/(iω) calculated from the MPEK software package(Hill et al. 2003a).
program based on the theory in § 4.3.3 and § 4.4.3. Comparing the spectra for
ν = 0.5 and ν < 0.5 provides an important consistency check on our numeri-
cal computations, since the methodologies for compressible and incompressible
skeletons are distinct. In general, Z/E can vary with Poisson’s ratio by up
to an order of magnitude at frequencies below the reciprocal draining time.
The response is clearly very sensitive to Poisson’s ratio as ν → 0.5 at these
frequencies. At higher frequencies, however, Z/E is independent of Poisson’s
ratio with fixed shear modulus µ, because the compressible polymer skeleton
is hydrodynamically coupled to the incompressible fluid.
Next, figure 4–5 shows how Young’s modulus affects the response spec-
trum. As expected from the steady displacement (Hill & Ostoja-Starzewski
2008; Wang & Hill 2008), the response is indeed inversely proportional to
the elastic modulus at frequencies below the reciprocal viscous time τ−1v . In
addition, the elastic modulus changes both the draining and viscous times.
Accordingly, the spectra in figure 4–5 overlap at frequencies below τ−1v when
82
10-2
100
102
104
106
108
frequency, (Hz)
0
0.1
0.2
0.3
0.4
-Z/E
, [nm
/(V
/cm
)]
10-2
100
102
104
106
10810
-6
10-5
10-4
10-3
10-2
10-1
100
| Z/E
|, [
nm/(
V/c
m)]
0
0.1
0.2
0.3
0.4
0.5
Figure 4–4: Frequency spectra of Z/E for various Poisson ratios ν = 0, 0.1,0.2, 0.3, 0.4, and 0.5 with a fixed shear modulus µ = E/(2 + 2ν) ≈ 333 Pa.All other parameters are the same as in figure 4–2. Solid and dashed lines thereal and imaginary parts of Z/E from numerically exact solutions of the fullmulti-phase electrokinetic model.
multiplying Z/E by E and dividing the frequency ω/(2π) by E . Noteworthy
from the perspective of electroacoustics is that the real part of Z/E (solid
lines) is sensitive to E at ultrasonic frequencies, whereas the imaginary part
(dashed lines) is practically independent of E .
The influences of the scaled ζ-potential ζe/(kT ) and scaled reciprocal
double layer thickness κa on Z/E are examined in figure 4–6. Note that our
computational methodology is stable and accurate for all κa & 1. When κa .
1, however, electro-osmotic flow is extremely weak and the dynamic response
can be approximated by (4.21). Situations of practical significance most often
occur when κa & 1, so figure 4–6 presents spectra for six values of κa in the
range 1–500. When κa is large, electroosmotic flow significantly affects the
particle displacement in the quasi-steady compressible elastic regime in the
same manner as for steady electric fields (Wang & Hill 2008). Consequently,
the displacement at frequencies below the reciprocal draining time increases
with κa relative to the respective quasi-steady elastic plateaus (incompressible
regimes).
83
100
102
104
106
108
frequency, (Hz)
10-7
10-6
10-5
10-4
10-3
10-2
10-1
100
101
| Z/E
|, [
nm/(
V/c
m)]
0.1 kPa
1 kPa
10 kPa
Figure 4–5: Frequency spectra of Z/E for various Young’s moduli E = 100 Pa,1 kPa, and 10 kPa with Poisson ratio ν = 0.2. All other parameters are thesame as in figure 4–2. Solid and dashed lines the real and imaginary partsof Z/E from numerically exact solutions of the full multi-phase electrokineticmodel.
The response Z/E is also affected by polarization and relaxation of the
diffuse double layer. This is especially evident for particles with thick double
layers and high ζ-potentials, a situation where the back-field of the polarized
double layer is strong (Gibb & Hunter 2000). For example, in figure 4–6(a)
with κa = 1 and |ζ| & 2kT/e, the real part of Z/E increases with frequency,
and the imaginary part changes sign between approximately 2 and 30 kHz.
These changes occur at frequencies higher than the reciprocal diffusion re-
laxation time τ−1i , which represents the maximum frequency that the diffuse
double layer is capable of following the external field (DeLacey & White 1981).
The increase in the real part of Z/E indicates that the back-field decreases
with increasing frequency, thereby reducing the so-called retardation experi-
enced by the particle. The sign change of the imaginary part of Z/E indicates
that the double layer polarization lags the applied field when the frequency is
higher than τ−1i . For particles with large κa, the back-field is weak, because
relaxation via diffusion across a thin double layer is fast. In the panels with
κa = 5, 10, and 50, maximums in the real part of Z/E with respect to |ζ| are
84
100
102
104
106
108
frequency, (Hz)
0
0.005
0.01
0.015
0.02
-Z/E
, [nm
/(V
/cm
)]
-1
-2
-4
-6
-8
(a) κa = 1
100
102
104
106
108
frequency, (Hz)
0
0.005
0.01
0.015
0.02
0.025
0.03
-Z/E
, [nm
/(V
/cm
)]
100
102
104
106
10810
-7
10-6
10-5
10-4
10-3
10-2
10-1
| Z/E
|, [
nm/(
V/c
m)]
-1
-2
-4
-6 -8
(b) κa = 5
100
102
104
106
108
frequency, (Hz)
0
0.01
0.02
0.03
0.04
0.05
-Z/E
, [nm
/(V
/cm
)]
100
102
104
106
10810
-7
10-6
10-5
10-4
10-3
10-2
10-1
100
| Z/E
|, [
nm/(
V/c
m)]
-1
-2
-4
-6
-8
(c) κa = 10
100
102
104
106
108
frequency, (Hz)
0
0.05
0.1
0.15
0.2
-Z/E
, [nm
/(V
/cm
)]
100
102
104
106
10810
-7
10-6
10-5
10-4
10-3
10-2
10-1
100
| Z/E
|, [
nm/(
V/c
m)]
-1
-2
-4
-6
-8
(d) κa = 50
100
102
104
106
108
frequency, (Hz)
0
0.05
0.1
0.15
0.2
0.25
-Z/E
, [nm
/(V
/cm
)]
100
102
104
106
10810
-7
10-6
10-5
10-4
10-3
10-2
10-1
100
| Z/E
|, [
nm/(
V/c
m)]
-1
-2
-4
-6
-8
(e) κa = 100
100
102
104
106
108
frequency, (Hz)
0
0.1
0.2
0.3
0.4
0.5
-Z/E
, [nm
/(V
/cm
)]
100
102
104
106
10810
-7
10-6
10-5
10-4
10-3
10-2
10-1
100
| Z/E
|, [
nm/(
V/c
m)]
-1
-2
-4
-6
-8
(f) κa = 500
Figure 4–6: Frequency spectra of Z/E for various scaled ζ-potentials−ζe/(kT ) = 1, 2, 4, 6, and 8; and various scaled reciprocal double-layerthicknesses κa = 1, 5, 10, 50, 100, and 500: KCl at T = 298 K; a/` = 100(a = 500 nm and ` = 5 nm); ρp = 1050 kg m−3; ν = 0.2; and E = 1 kPa. Solidand dashed lines the real and imaginary parts of Z/E from numerically exactsolutions of the full multi-phase electrokinetic model.
85
evident. These can lead to ambiguity in determining the ζ-potential from the
steady response (Hill & Ostoja-Starzewski 2008; Wang & Hill 2008). However,
by measuring the frequency spectrum of Z/E (or the mobility µd = −iωZ/E),
it may be easier to unambiguously ascertain the correct ζ-potential. This ap-
proach has been useful for interpreting electroacoustic measurements of the
dynamic mobility (Hunter & O’Brien 1997).
Finally, figures 4–7 and 4–8, respectively, show the influence of hydrogel
permeability `2 for small and large values of κa. The Brinkman screening
length ` has a significant influence on the hydrodynamic coupling between
the fluid and polymer skeleton. Earlier studies (Hill & Ostoja-Starzewski
2008; Wang & Hill 2008) demonstrate that the polymer displacement at steady
state is practically independent of ` for incompressible hydrogels, but varies
significantly for compressible hydrogels due to an adverse electroosmotic-flow-
induced pressure gradient, particularly when κa is large. Note that the Brinkman
screening length also affects the draining time.
When κa is small (figure 4–7), the Brinkman screening length is most ef-
fective in changing reciprocal draining time τ−1d . Accordingly, as the Brinkman
screening length increases, the frequency range exhibiting a quasi-steady in-
compressible elastic response decreases, and eventually disappears, with Z/E
transferring directly from the quasi-steady compressible plateau to the vis-
cous and inertial stress dominated regimes. When κa is large (figure 4–8),
the Brinkman screening length also significantly changes the amplitude of the
quasi-steady compressible plateau. Similarly to the steady displacement (Hill
& Ostoja-Starzewski 2008; Wang & Hill 2008), decreasing the permeability
increases the magnitude of the adverse tangential pressure gradient, which, in
turn, increases the particle displacement.
4.7.3 High-frequency boundary-layer approximation and applica-tion to electroacoustics
The amplitude of the dynamic electrokinetic response Z/E was demon-
strated above to become extraordinarily small at high frequencies. The results
86
100
102
104
106
108
frequency, (Hz)
0
0.001
0.002
0.003
0.004
0.005
-Z/E
, [nm
/(V
/cm
)]
1 10
1 10
50
50
100
1005
5
(a) κa = 1, ζe/(kT ) = −1
100
102
104
106
108
frequency, (Hz)
0
0.004
0.008
0.012
0.016
-Z/E
, [nm
/(V
/cm
)]
1 10
1 10
50
50
100
100
5
5
(b) κa = 1, ζe/(kT ) = −4
100
102
104
106
108
frequency, (Hz)
0
0.002
0.004
0.006
0.008
0.01
0.012
-Z/E
, [nm
/(V
/cm
)]
100
102
104
106
10810
-710
-610
-510
-410
-310
-210
-1
| Z/E
|, [
nm/(
V/c
m)]
1 10
1 10
50
50
100
1005
5
(c) κa = 10, ζe/(kT ) = −1
100
102
104
106
108
frequency, (Hz)
0
0.01
0.02
0.03
0.04
-Z/E
, [nm
/(V
/cm
)]
100
102
104
106
10810
-710
-610
-510
-410
-310
-210
-1
| Z/E
|, [
nm/(
V/c
m)]
110
1 10
50
50
100
100
5
5
(d) κa = 10, ζe/(kT ) = −4
Figure 4–7: Frequency spectra of Z/E for various Brinkman screening lengths` = 1, 5, 10, 50, and 100 nm with scaled ζ-potentials −ζe/(kT ) = 1 and4, and scaled reciprocal double-layer thickness κa = 1 and 10. All otherparameters are the same as in figure 4–6. Solid and dashed lines the real andimaginary parts of Z/E from numerically exact solutions of the full multi-phase electrokinetic model.
87
100
102
104
106
108
frequency, (Hz)
0
0.01
0.02
0.03
0.04
0.05
0.06
-Z/E
, [nm
/(V
/cm
)]
100
102
104
106
10810
-7
10-6
10-5
10-4
10-3
10-2
10-1
| Z/E
|, [
nm/(
V/c
m)]1
10
50
100
5
(a) κa = 100, ζe/(kT ) = −1
100
102
104
106
108
frequency, (Hz)
0
0.05
0.1
0.15
0.2
0.25
0.3
-Z/E
, [nm
/(V
/cm
)]
100
102
104
106
10810
-8
10-6
10-4
10-2
100
| Z/E
|, [
nm/(
V/c
m)]
1
10
50100
5
(b) κa = 100, ζe/(kT ) = −4
100
102
104
106
108
frequency, (Hz)
0
0.05
0.1
0.15
0.2
0.25
-Z/E
, [nm
/(V
/cm
)]
100
102
104
106
10810
-7
10-6
10-5
10-4
10-3
10-2
10-1
100
| Z/E
|, [
nm/(
V/c
m)]
1
1050 100
5
(c) κa = 1000, ζe/(kT ) = −1
100
102
104
106
108
frequency, (Hz)
0
0.2
0.4
0.6
0.8
1
-Z/E
, [nm
/(V
/cm
)]
100
102
104
106
10810
-7
10-6
10-5
10-4
10-3
10-2
10-1
100
| Z/E
|, [
nm/(
V/c
m)]
1
10 50 100
5
(d) κa = 1000, ζe/(kT ) = −4
Figure 4–8: Frequency spectra of Z/E for various Brinkman screening lengths` = 1, 5, 10, 50, and 100 nm with scaled ζ-potentials −ζe/(kT ) = 1 and 4,and scaled reciprocal double-layer thickness κa = 100 and 1000. All otherparameters are the same as in figure 4–6. Solid and dashed lines the real andimaginary parts of Z/E from numerically exact solutions of the full multi-phase electrokinetic model.
88
below demonstrate that the dynamic mobility µd = −iωZ/E is large at the fre-
quencies used in commercial electroacoustic instruments. In § 4.5 we showed
that O’Brien’s macroscopic electroacoustic equations (O’Brien 1990) can be
applied to hydrogel-colloid composites. Accordingly, a close connection was
established between the electroacoustic signal (pressure fluctuations) in ESA
measurements and the dynamic mobility. This motivated the derivation in
§ 4.6 leading to the analytical approximation (4.126) for the high frequencies
encountered in electroacoustic experiments.
The dynamic electrophoretic mobilities for representative hydrogel-colloid
composites from (4.126) and numerically exact computations are presented in
figure 4–9. Because our calculations neglect interactions, they are suitable
for composites with low particle volume fractions. Spectra are shown with
Young’s modulus spanning three decades. Note that the spectrum with a finite
(real part) plateau at low frequencies is the mobility for the same particles
dispersed in a Newtonian electrolyte (without polymer); this was calculated
using the MPEK software package (Hill et al. 2003a). At this large value of
κa = 500, the analytical approximation (4.126) compares extremely well with
the numerically exact calculations. Note also that the real parts (dash-dotted
and solid lines) depart very slightly when the amplitude vanishes at lower
frequencies, and the imaginary parts are practically identical at all frequencies.
Noteworthy is that the real part of µd can be distinguished from µd∗
only at MHz frequencies when Young’s modulus of the skeleton is greater
than about 10 kPa. However, the imaginary part is very sensitive to the
elastic modulus at about 1 MHz, suggesting that, at a given fixed frequency,
or in a narrow range of frequencies, an electroacoustic experiment could probe
the kinetics of polymer gelation and aging on time scales less than a second.
Recall, commercial electroacoustic instruments operate at frequencies between
0.3 MHz and 11 MHz (Hunter 1998). Clearly, to probe the elastic modulus of
hydrogel skeletons with lower moduli, a wider frequency range—extending to
lower frequencies—is required. Nevertheless, many important hydrogels have
89
103
104
105
106
107
108
109
frequency, (Hz)
-8
-6
-4
-2
0
2
4
6
mob
ility
, [(µ
m/s
)/(V
/cm
)]
0.1 kPa 1 kPa 10 kPa
Figure 4–9: Frequency spectra of the dynamic electrophoretic mobility µd =−iωZ/E for Young’s moduli E = 100 Pa, 1 kPa, and 10 kPa with Poisson’sratio ν = 0.2. Other parameters are the same as in figure 4–2. Solid anddashed lines are the real and imaginary parts of µd = −iωZ/E from numeri-cally exact solutions of the full multi-phase electrokinetic model. Dash-dottedlines are the real parts of µd from the analytical approximation (4.126); theimaginary parts of µd from (4.126) are not shown because they overlap thenumerically exact solution. Note that the spectrum with a finite (real) low-frequency plateau is the mobility µd
∗ = −iω(Z/E)∗ for the same particlesdispersed in a Newtonian electrolyte (without polymer), calculated using theMEPK software package (Hill et al. 2003a).
an elastic modulus greater than 10 kPa, e.g., reverse thermoresponsive poly(N-
isopropylacrylamide) gels at 8 wt% and 40 C have E ≈ 170 kPa (Takigawa
et al. 1997), and polyacrylamide gels at 0.8% w/v have E ≈ 35 kPa (Takigawa
et al. 1996), so the real and imaginary parts of their electroacoustic responses
at MHz frequencies would be extremely sensitive to changes in the elastic
modulus. Note also that, because the electroacoustic response is prominent
at frequencies much higher than the reciprocal draining time, the response
only reflects changes in the shear modulus, not the accompanying changes in
hydrodynamic permeability.
Figure 4–10 compares numerically exact calculations of the dynamic elec-
trophoretic mobility with our analytical boundary-layer approximation (4.126)
90
with κa = 50 (top panels) and 1000 (bottom panels). The boundary-layer ap-
proximation is accurate when κa = 1000 for all the experimentally accessible
ζ-potentials and frequencies, and even when κa = 50 the approximation de-
viates only slightly from the exact calculations at the highest ζ-potentials
and frequencies. These dynamics are independent of the Brinkman screening
length and Poisson ratio (with fixed shear modulus). Note that (4.126) does
not capture the compressible dynamics at frequencies below the reciprocal
draining time, where the dynamic mobility is vanishingly small.
4.8 Summary
We extended the multi-phase electrokinetic model of Hill & Ostoja-Starzewski
(2008) and Wang & Hill (2008) to calculate the dynamic response of a charged,
spherical colloid embedded in uncharged hydrogels subjected to harmonically
oscillating electric fields. We began by solving the two-fluid model of Levine
& Lubensky (2001) exactly, and compared our analytical solution with two
approximations widely adopted in the microrheology literature. We then de-
veloped a powerful computational methodology to solve the full multi-phase
electrokinetic model by linearly perturbing an equilibrium base state governed
by the non-linear Poisson-Boltzmann equation. The particle response, defined
as the ratio of the displacement to the electric field strength, was obtained
by superposing two simpler sub-problems to satisfy the particle equation of
motion. Compressible and incompressible hydrogel skeletons had to be con-
sidered separately. By adopting an analytical solution in the far field, we
achieved accurate numerical solutions over an extraordinarily wide range of
frequencies, in a wide range of the experimentally accessible parameter space.
In addition, we examined the dynamic electrophoretic mobility, defined as the
ratio of the particle velocity to the electric field strength, and its connection
to electroacoustic diagnostics for characterizing hydrogel-colloid composites.
Noteworthy was an analytical boundary-layer approximation that compares
extremely well with the numerically exact results at the ultra-sonic frequen-
cies adopted in commercial electroacoustic instruments.
91
103
104
105
106
107
108
109
frequency, (Hz)
-8
-7
-6
-5
-4
-3
-2
-1
0
1
mob
ility
, [(µ
m/s
)/(V
/cm
)]
-2
-1
-4
-6
-8
(a) κa = 50, Re(µd)
103
104
105
106
107
108
109
frequency, (Hz)
-3
-2
-1
0
1
2
3
4
5
mob
ility
, [(µ
m/s
)/(V
/cm
)]
-2
-1
-4
-6
-8
(b) κa = 50, Im(µd)
103
104
105
106
107
108
109
frequency, (Hz)
-16
-14
-12
-10
-8
-6
-4
-2
0
2
mob
ility
, [(µ
m/s
)/(V
/cm
)]
-1
-2
-4
-6
-8
(c) κa = 1000, Re(µd)
103
104
105
106
107
108
109
frequency, (Hz)
-6
-4
-2
0
2
4
6
8
mob
ility
, [(µ
m/s
)/(V
/cm
)]
-1
-2
-4
-6
-8
(d) κa = 1000, Im(µd)
Figure 4–10: Frequency spectra of the dynamic electrophoretic mobility µd =−iωZ/E for various scaled ζ-potentials −ζe/(kT ) = 1, 2, 4, 6, and 8: κa = 50(top panels) and 1000 (bottom panels). All other parameters are the sameas in figure 4–6. Solid lines are the real (left panels) and imaginary (rightpanels) parts of µd from numerically exact solutions of the full multi-phaseelectrokinetic model, and dashed lines are µd from the analytical boundary-layer approximation (4.126).
92
The approximate solution of the two-fluid model by Levine & Lubensky
(2001) agrees well with our exact analytical solution when fluid and particle
inertia can be neglected. However, the range of applicability of the gener-
alized Stokes-Einstein relation (GSER) was found to be narrower than pre-
viously thought. The electric-field-induced dynamic response of a colloidal
particle in a hydrogel often exhibits an ostensible transition from quasi-steady
compressible to incompressible elastic dynamics—both characterized by dis-
tinct plateaus in the real part of the frequency spectrum—as the frequency
passes through the reciprocal draining time τ−1d of the hydrogel. At higher
frequencies, when the dynamics are dominated by viscous and inertial forces,
the response is similar to a particle in a Newtonian electrolyte. In general,
the response depends on Poisson’s ratio, Young’s modulus, and Brinkman
screening length of the hydrogel, as well as physicochemical characteristics of
the inclusions, including size and charge. At frequencies above the reciprocal
draining time, the response is practically independent of hydrogel permeabil-
ity and compressibility, since the the fluid and polymer skeleton are strongly
coupled by hydrodynamic drag forces. At frequencies below the reciprocal
draining time, hydrogel compressibility can increase the electric-field-induced
particle displacement by an order of magnitude relative to the displacement in
a perfectly incompressible skeleton with the same shear modulus. Accordingly,
the dynamics of compressible and incompressible hydrogels are qualitatively
different at low frequencies. Note that the response spectrum also reflects po-
larization and relaxation of the diffuse double layer, particularly for inclusions
with thick double layers (small κa) and high ζ-potentials.
The present theory provides a rigorous foundation for interpreting two
classes of electric-field-based diagnostic experiments involving hydrogel-colloid
composites. Such experiments probe both the physicochemical characteris-
tics of the charged inclusions, and the viscoelastic rheology of the hydrogel.
Our calculations demonstrate that the particle displacement at low frequen-
cies could be directly measured using active electrical microrheology. At
93
higher frequencies, however, the particle displacements are too small (sub-
nanometer) to measure directly, so electroacoustic techniques are necessary to
measure instead the dynamic electrophoretic mobility of the inclusions. Ac-
cordingly, we showed that the macroscopic relations for colloidal dispersions
developed by O’Brien (1990) can be directly applied to hydrogel-colloid com-
posites. Our calculations of the dynamic electrophoretic mobility demonstrate
that the strength of electroacoustic signals from hydrogel-colloid composites
are comparable to those from Newtonian electrolytes without a polymer skele-
ton. Accordingly, our calculations suggest that electroacoustic experiments
on hydrogel-colloid composites could be performed using presently available
commercial instruments.
This chapter and chapter 3 provide the complete solution of the multi-
phase electrokinetic model for dynamic and steady electric-field-induced re-
sponses of colloidal particles embedded in hydrogels. The particle responses
are sensitive to the hydrogel viscoelasticity, compressibility and hydrodynamic
permeability, and physicochemical properties of the inclusion. Therefore, char-
acterization techniques based on these responses can be developed. Clearly,
the continuous phase significantly affects particle dynamics. Effects of colloidal
particles on the bulk composite properties are revealed in the next chapter with
a different problem, where the nanoparticle-induced anomalous bulk viscosity
reductions found in polymer-nanocomposite melts are studied.
Appendicies
4.A Point-force representation of a particle in an uncharged hy-drogel matrix
Here we relate the net force on a spherical colloid in an uncharged hydrogel
to the strength of a point force that produces the same far-field disturbances.
The particle undergoes harmonic translation in an otherwise stationary hy-
drogel. The strength of the point force is obtained from reciprocal relations
similar to Hill et al. (2003a). However, in addition to the Lorentz reciprocal
94
relation for fluid in a domain S (Kim & Karrila 1991),∫∂S
(u′ · Tf − u · Tf ′) · ndA =
∫S
(u′ ·∇ · Tf − u ·∇ · Tf ′)dV, (4.129)
a similar reciprocal relation, known in solid mechanics as the Betty theo-
rem (Barber 2003) is required. This is∫∂S
(µ′v′ · Te − µv · Te′) · ndA =
∫S
(µ′v′ ·∇ · Te − µv ·∇ · Te′)dV, (4.130)
provided λ = λ′.
Consider a large domain Ω with boundary ∂Ω and outward unit normal
n that encloses an oscillating sphere centered at position r1 with radius a
(system 1) and a fixed point-force centered at position r2 (system 2). Note that
|r1−r2| a+κ−1. Furthermore, the sphere occupies volume Ω1 and undergoes
oscillatory translation with velocity −iωZ. The corresponding surface and
outward unit normal are denoted ∂Ω1 and n1, respectively.
The divergence of elastic and hydrodynamic stresses for system 1 (∇ ·Tf1
and ∇ · Te1) are given by (4.58b) and (4.58c), respectively. For system 2,
∇ · Tf2 = η∇2u2 −∇p2
= −iωρfu2 + (η/`2)(u2 + iωv2) + f fδ(r2), (4.131a)
∇ · Te2 = µ∇2v2 + (λ+ µ)∇(∇ · v2)
= f eδ(r2)− (η/`2)(u2 + iωv2), (4.131b)
where δ(r) is the Dirac-delta function, and f f and f e are the point forces
exerted on the fluid and elastic medium, respectively.
Applying the Lorentz reciprocal relation to the volume enclosed by ∂Ω1
and ∂Ω gives∫∂Ω
(u1 · Tf2 − u2 · Tf
1) · ndA−∫∂Ω1
(u1 · Tf2 − u2 · Tf
1) · n1dA
=
∫Ω−Ω1
(u1 ·∇ · Tf2 − u2 ·∇ · Tf
1)dV. (4.132)
95
Because u ∼ r−3 as r → ∞, the integral over ∂Ω on the left-hand side of
(4.132) vanishes when Ω is sufficiently large. Therefore, substituting (4.58b)
and (4.131a) into (4.132) gives
−∫∂Ω1
(u1 · Tf2 − u2 · Tf
1) · n1dA
= u1(r2) · f f +
∫Ω−Ω1
[iω(η/`2)(u1 · v2 − u2 · v1) + u2 ·∇ · Tm1 ]dV.(4.133)
Inside the particle, u1 = −iωZ, and since |r1 − r2| a + κ−1, u2(|x −
r1| ≤ a) can be considered constant. Therefore, applying Gauss’s divergence
theorem to the integral over ∂Ω1 on the left-hand side of (4.133) gives∫∂Ω1
(u1 · Tf2 − u2 · Tf
1) · n1dA
= −iωZ ·∫
Ω1
∇ · Tf2dV − u2(r1) ·
∫∂Ω1
Tf1 · n1dA
= −u2(r1) ·[ω2ρfVpZ +
∫∂Ω1
Tf1 · n1dA
], (4.134)
where Vp is the particle volume. Substituting (4.134) into (4.133) gives
f f · u1(r2) = u2(r1) ·[ω2ρfVpZ +
∫∂Ω1
Tf1 · n1dA
]−∫
Ω−Ω1
[iω(η/`2)(u1 · v2 − u2 · v1) + u2 ·∇ · Tm1 ]dV. (4.135)
Similarly, for the elastic displacements, applying the same procedure as above,
but with the Betty theorem, yields
f e ·v1(r2) = v2(r1) ·∫∂Ω1
Te1 ·n1dA−
∫Ω−Ω1
(η/`2)(u1 ·v2−u2 ·v1)dV. (4.136)
where v2(|x− r1| ≤ a) may be considered constant.
Again, since u1(r2) = −iωv1(r2) and u2(r1) = −iωv2(r1), multiplying
(4.136) by −iω and adding (4.135) gives
f ·u1(r2) = u2(r1)·[ω2ρfVpZ +
∫∂Ω1
(Tf1 + Te
1) · n1dA
]−∫
Ω−Ω1
u2 ·∇·Tm1 dV,
(4.137)
96
where f = f e + f f is the total point force. Note that the volume integral on
the right-hand side of (4.137) can be factored to give∫Ω−Ω1
u2 ·∇ · Tm1 dV = −u2(r1) ·
∫∂Ω1
Tm1 · n1dA, (4.138)
because ∇ · Tm is exponentially small when |x − r1| a + κ−1, and u2 can
be considered constant where ∇ · Tm is finite.
Substituting (4.138) into (4.137) yields
f = ω2ρfVpZ +
∫∂Ω1
(Te1 + Tf
1 + Tm1 ) · n1dA. (4.139)
The integral over ∂Ω1 on the left-hand side of (4.139) is the total force on the
sphere, which according to Newton’s second law must equal the acceleration
of its mass, −Vpρpω2Z, so the strength of the point force is
f = ω2VpZ(ρf − ρp). (4.140)
In other words, similarly to bare particles (Mangelsdorf & White 1992) and
particles with polymer coatings (Hill et al. 2003a) dispersed in Newtonian
electrolytes, the acceleration of the mass of fluid displaced by a finite sized
inclusion in a hydrogel must be added to the force on a point particle producing
the same far-field fluid velocity and polymer displacement disturbances. Note
that the foregoing analysis neglects the mass of the polymer.
4.B Numerical solution of the field equations
The field equations are solved according to the outline presented in § 4.3.
First, the Poisson-Boltzmann equation is solved efficiently using the adaptive
mesh algorithm developed by Hill et al. (2003a), and then various equilibrium
quantities and their derivatives are computed. Next, the linearly perturbed
equations are solved. Before this calculation, the equations are transformed to
simplify the numerical methods discussed in § 4.3 and § 4.4. The matrix alge-
bra and eigenvalue calculations involved in the transformations are performed
97
using BLAS and LAPACK routines (Anderson et al. 1999). Completely differ-
ent computational strategies are adopted for calculating the linearized pertur-
bations for incompressible and compressible hydrogel skeletons. Asymptotic
coefficients and physical quantities are then constructed from the numerical so-
lutions. The following discusses in further detail how the perturbed problems
are solved.
For incompressible skeletons, the perturbed equations are transformed to
the decoupled forms outlined in § 4.3, and the resulting differential equations
are discretized using a second-order central difference scheme, and solved using
a banded matrix solver. Solutions are then improved iteratively using a mov-
ing mesh method based on the methodology of Hill et al. (2003a). When the
solution has converged, asymptotic coefficients based on the far-field asymp-
totic analysis are obtained. The far-field solution is calculated using LAPACK.
Our program to compute the response for incompressible skeletons is written
entirely in C.
For compressible skeletons, the perturbed solutions oscillate in space with
several wave lengths, e.g., the construction of the fluid velocity and polymer
displacement in § 4.2 involves three wave lengths. The second-order central
difference scheme with the moving mesh method of Hill et al. (2003a) does not
converge. We therefore modified a general-purpose boundary value problem
software package TWPBVPL (Cash & Mazzia 2006), which solves the differ-
ential equations using fourth-, sixth- and eighth-order methods with hybrid
mesh selection, to solve the linearly perturbed problem. The second-order or-
sion times a2/Ds ∼ 0.1 and 10 s when a = 10 nm, which are much shorter
than ω−1 and γ−1 in the experiments.
5.3 Intrinsic viscosity from the single-layer model
The intrinsic viscosity from (5.8) depends on three dimensionless param-
eters: the scaled layer thickness δ/a, the scaled reciprocal slipping length
ka, and a layer property parameter χ = (ηi/ηo)/(ρi/ρo). In polymer solu-
tions (de Gennes 1979; Larson 1999), the bulk viscosity increases with the
polymer concentration according to a power law. Considering the similarity
between polymer melts and solutions, it is natural to assume that the viscos-
ity in the layer varies with the polymer segment density as η ∼ ρn (n > 0).
Consequently, χ = (ρi/ρo)n−1, and when ρi < ρo, it follows that χ < 1 if n > 1.
Figures 5–1 and 5–2 present the effect of scaled layer thickness δ/a on the
intrinsic viscosity [η] for χ ≤ 1 and χ ≥ 1, respectively. Evidently, when the
layer properties are identical to those of the bulk polymer melt, i.e., χ = 1, [η]
is independent of the layer thickness. Note that ka = 0 gives [η] = 1, which
mimics the intrinsic viscosity of an inviscid spherical bubble, and ka→∞ gives
[η] → 5/2, which recovers the hard-sphere result from Einstein (1906). With
increasing layer thickness δ/a, the intrinsic viscosity decreases from positive
112
10-3
10-2
10-1
100
101
Scaled layer thickness, δ/a
10-2
10-1
100
101
102
Intr
insi
c vi
scos
ity, [
η]0
10-4
1
0.01
0.5
0.1
(a) ka = 0, perfect slip
10-3
10-2
10-1
100
101
Scaled layer thickness, δ/a
10-2
10-1
100
101
102
Intr
insi
c vi
scos
ity, [
η]
0
10-4
1
0.01
0.5
0.1
(b) ka→∞, no slip
Figure 5–1: Intrinsic viscosity [η] as a function of the scaled layer thicknessδ/a according to (5.8) for polymer-nanocomposite melts with χ = 0, 10−4,0.01, 0.1, 0.5 and 1; ka = 0 (left panel) and ka→∞ (right panel). Solid anddashed lines indicate positive and negative values of [η], respectively.
to negative values when χ < 1, and increases when χ > 1. In both cases, [η]
asymptotes to the results with χ = 1 for small δ/a, and its absolute magnitude
increases with increasing layer thickness. Moreover, comparing results for
ka = 0 and ka→∞ reveals that the intrinsic viscosity is relatively insensitive
to the degree of slip when δ/a is large. Here, the case with χ ≤ 1 shown in
figure 5–1 is more important, because it exhibits negative intrinsic viscosities.
When the magnitude of the intrinsic viscosity is small, the corresponding
layer thickness δ/a is sensitive to χ, but this sensitivity vanishes when the
intrinsic viscosity is large. For example, in figure 5–1(a), the layer thicknesses
corresponding to [η] = −1 with χ = 10−4 and 0.5 are δ/a = 0.03 and 0.5,
respectively, but the layer thicknesses corresponding to [η] = −100, for the
same χ, are δ/a = 3 and 4, indicating a reduced sensitivity to χ at large
by Mackay et al. (2003) and Tuteja et al. (2005) correspond to large layer
thicknesses, and are therefore relatively insensitive to χ.
Two limiting cases, χ → 0 and χ → ∞, with large δ/a are identified
in figures 5–1 and 5–2, respectively. As χ → ∞, the layer becomes rigid
113
10-3
10-2
10-1
100
101
Scaled layer thickness, δ/a
100
101
102
Intr
insi
c vi
scos
ity, [
η]
11.1
25100
∞
Figure 5–2: Intrinsic viscosity [η] as a function of the scaled layer thicknessδ/a according to (5.8) for polymer-nanocomposite melts with ka = 0 (solidlines) and ka→∞ (dashed lines); χ = 1, 1.1, 2, 5, 100 and ∞.
(infinitely viscous), so nanoparticles behave as rigid spheres with radius a+ δ
in the melt, giving
[η]→ (5/2)(1 + δ/a)3 as χ→∞. (5.11)
Similarly, as χ → 0, nanoparticles behave as instantaneously spherical, de-
formable bubbles with radius a+ δ, giving
[η]→ −(5/3)(1 + δ/a)3 as χ→ 0. (5.12)
Note that these limiting cases are independent of slip at the nanoparticle-
polymer interface.
The influence of the property parameter χ on the intrinsic viscosity is
shown in figure 5–3. These results are for a no-slip nanoparticle-polymer
interface, i.e., ka→∞. As expected, [η] = 5/2 when χ = 1, and the intrinsic
viscosity asymptotes to the values given by (5.12) and (5.11) when χ→ 0 and
∞, respectively.
Figure 5–4 demonstrates how the scaled reciprocal slipping length ka
affects [η] at χ = 0.1 for several values of δ/a. The slipping parameter is
114
10-3
10-2
10-1
100
101
102
103
Dimensionless parameter, χ
10-2
10-1
100
101
102
103
104
Intr
insi
c vi
scos
ity, [
η]
0.1
0.3
1
3
10
Figure 5–3: Intrinsic viscosity [η] as a function of the parameter χ =(ηi/ηo)/(ρi/ρo) according to (5.8) for polymer-nanocomposite melts withka → ∞ (no-slip) and δ/a = 0.1, 0.3, 1, 3, and 10. Solid and dashed linesindicate positive and negative values of [η], respectively.
clearly most influential with thin layers, where the magnitude of [η] is small,
and has negligible effect when δ/a & 1. For the experiments of Mackay et al.
(2003) and Tuteja et al. (2005) where [η] ∼ −100, slipping at the nanoparticle-
polymer interface does not significantly affect [η].
From the results above, the model demonstrates that the layer thickness
must be comparable to or larger than the nanoparticle radius to achieve the
negative intrinsic viscosities measured experimentally. Under these conditions,
the intrinsic viscosity is almost independent of ka and χ when χ . 0.1. It
should therefore be reasonable to interpret the experiments by focusing on the
influence of layer thickness, as captured by the single parameter δ/a.
5.4 Theoretical interpretation of experiments
Here we interpret the experiments of Mackay et al. (2003) and Tuteja
et al. (2005) using the model presented in § 5.2 and § 5.3. To ensure the
polymer molecules are entangled and confined (Tuteja et al. 2005), only the
experiments with entangled polymer melts and Rg/h & 1 are considered. Here
h is the interparticle half gap defined by Tuteja et al. (2005), and Rg/h & 1
indicates the polymers are confined by nanoparticles. Also, since our model
115
10-2
10-1
100
101
102
Scaled reciprocal slipping length, ka
10-3
10-2
10-1
100
101
102
103
104
Intr
insi
c vi
scos
ity, [
η]
0.01
0.1
0.15
1
10
Figure 5–4: Intrinsic viscosity [η] as a function of the scaled reciprocal slippinglength ka according to (5.8) for polymer-nanocomposite melts with χ = 0.1and δ/a = 0.01, 0.1, 0.15, 1, and 10. Solid and dashed lines indicate positiveand negative values of [η], respectively.
does not account for hydrodynamic (and other) interactions, we restrict our
analysis to the lowest nanoparticle volume fraction, i.e., φ = 0.005. From
the estimated viscosity and density in the layers, the layer thicknesses are
obtained by fitting the theory to experiments. Such comparisons shed light
on the mechanisms of bulk viscosity reduction.
We propose two possible origins for the layers. Firstly, these may arise
from changes in free volumes induced by nanoparticles. In this case, polymer
segments are depleted at the nanoparticle surfaces due to excluded volume
effects, and the melt viscosity and density recover their bulk values over a
length scale of the polymer correlation length ξ. For polymer melts, ξ ≈
b ∼ 1 nm (Fleer et al. 1993), where, recall, b is the segment size. Using
this idea, Hill (2006a,c) successfully explained the enhanced permeability and
reverse selectivity of nanocomposite membranes (Merkel et al. 2002). The
segment density profile can be calculated by mean-field theories (de Gennes
1979; Wu et al. 1995) and computer simulations (Daoulas et al. 2005). Here,
as a first approximation, we assume that the polymer segment density inside
a uniform layer is half of the bulk value, i.e., ρi/ρo = 1/2. This is equivalent
116
to assuming a linear polymer density profile. The simple scaling of de Gennes
(1976b) for entangled polymer solutions connects the viscosity and segment
density in the layer by η ∼ ρ15/4. Therefore, the layer property parameter
χ ≈ (ρi/ρo)11/4 ≈ 0.15.
Another possible origin for the layers is the local Rouse viscosity of the
melt in close proximity of the nanoparticles. In this case, the layer has the
Rouse viscosity and bulk polymer density, and the layer thickness is charac-
terized by the tube diameter dt ∼ 10 nm (Dealy & Larson 2006). This is
advanced from the idea of Brochard-Wyart & de Gennes (2000), where they
proposed that nanoparticles smaller than the entanglement tube diameter dt
experience Rouse dynamics in entangled melts, since the nanoparticle motion
only involves simple chain rearrangements. The concept partially explained
the unusually high nanoparticle diffusivity in polymer-nanocomposite melts
(Tuteja et al. 2007b). The Rouse viscosity can be calculated according to
Tuteja et al. (2007b) as ηRouse = ηc(M/Mc)aT , where ηc = 292 Pa s is the PS
melt viscosity at the entanglement critical molecular weight Mc = 32.7 kDa,
and aT is the shift factor given by log(aT ) = −7.65[T (C)−170]/[T (C)−28.1].
Note that the actual PS melt viscosity can be calculated according to ηmelt =
ηc(M/Mc)3.68aT (Tuteja et al. 2005). Evidently, ηRouse ηmelt for entangled
melts at high molecular weight, and this yields χ ≈ 0.11 and 0.0015 for the
the 75 and 396 kDa polymer melts, respectively.
The preceding analysis shows that χ 1 regardless of whether the layer is
attributed to an increase in free volume or a decrease in viscosity due to Rouse
dynamics. Recall from § 5.3, the depletion layer thickness δ/a is insensitive to
χ and ka when the magnitude of the intrinsic viscosity is large. Consequently,
as a first approximation, we accept the surface-slip independent case when
χ = 0 to approximate the layer thickness with a uniform layer.
Table 5–1 presents the effective layer thickness δ and other important
quantities obtained by fitting the analytical model to the experiments of
Mackay et al. (2003) and Tuteja et al. (2005). The best-fit values of δ are: (1)
117
Table 5–1: Summary of the parameters that characterize the experiments ofMackay et al. (2003) and Tuteja et al. (2005) with Rg/h & 1 and φ = 0.005,and theoretical interpretations (providing fitted values for δ) based on (5.8)with χ = 0. Radii of gyration Rg are calculated from Cotton et al. (1974)and values of Rg/h are from Tuteja et al. (2005). Note that with propertemperature adjustment, the PS tube diameter dt ≈ 9.4 nm at 170C (Dealy& Larson 2006).
actions, because the thickness of the layer exceeds the entanglement tube
diameter for a pure melt.
Evidently, the experimental data giving (δ + a)/dt > 1/2 are always ac-
companied by small φm(1 + δ/a)−3 close to the actual nanoparticle volume
fraction φ = 0.005. This implies that distorted polymer configurations may
arise from nanoparticle interactions. Moreover, careful examination of data
from Tuteja et al. (2005) reveals that interactions also introduce strong free-
volume effects on the composite, as is evidenced by the large change in Tg.
Note that this is different from the free volumes introduced by excluded vol-
ume effects described above, which give rise to a layer thickness δ ∼ 1 nm,
and consequently a moderate negative intrinsic viscosity [η] ≈ −5. Clearly,
the large negative intrinsic viscosities observed by Mackay et al. (2003) and
Tuteja et al. (2005) are the result of strong polymer-mediated interactions
between nanoparticles and the polymer. Accordingly, we envision no clear
distinction between the bulk polymer and polymer in layers, so extending our
notion of perturbed polymer layers surrounding individual nanoparticles in an
unperturbed bulk polymer, we hypothesize that the polymer is everywhere
perturbed by the nanoparticles.
With weak interactions, the nano-scale Rouse dynamics in the polymer
entanglement tubes give rise to negative intrinsic viscosities much smaller in
magnitude than those with strong interactions. This is evidenced from the PS
75 kDa/25 kDa NP polymer nanocomposite in table 5–1, which shows (δ +
a)/dt < 1/2, φm(1 + δ/a)−3 φ, and the glass transition temperature change
∆Tg = 0.1C (Tuteja et al. 2005). Clearly, weak interactions explain the
thinnest depletion layer and moderate negative intrinsic viscosity achieved by
this sample, which is representative of a single nanoparticle in an unperturbed
entangled polymer melt.
Our uniform layer model can be extended to handle continuous polymer
density and viscosity layers by extending the uniform layer model to multiple
119
Table 5–2: Best-fit polymer correlation lengths ξ ascertained from experi-ments of Mackay et al. (2003) and Tuteja et al. (2005) and the theoreticalinterpretation based on a continuous-layer-profile model with ka = 0. Thepolymer segment densities are from Tuinier & Lekkerkerker (2002), and theviscosity-density relationship is from Colby et al. (1994).
layers, and discretizing the continuous profiles accordingly. However, a con-
tinuous change in the layer density and viscosity does not significantly modify
the qualitative pictures emerging from the foregoing analytical theory for uni-
form layers. This is demonstrated in table 5–2, which presents the best-fit
correlation lengths ξ for a continuous polymer segment density described by
the Edwards-de Gennes equation (de Gennes 1979), and the polymer viscos-
ity from the two-parameter scaling theory of Colby et al. (1994) for good
solvents. Since the spatial dependence of nano-scale Rouse viscosity around
nanoparticles in an entanglement tube is not available, we assume, as a first
approximation, that the layer arises from excluded volume effects. Note that
Edwards-de Gennes equation is the ground state dominance approximation of
Doi-Edwards theory (Doi & Edwards 1987), valid for infinitely long chains,
as well as long, but finite, chains near a hard wall (Wu et al. 1995). Here,
the segment density profile is described using the analytical approximate so-
lution from Tuinier & Lekkerkerker (2002). Moreover, the polymer melt in
the layer is considered a polymer solution, due to the similar polymer dynam-
ics (de Gennes 1979). Interestingly, even when the layers are attributed to
excluded volume effects, the resulting best-fit correlation lengths ξ are still
on the order of the entanglement tube diameter dt. This indicates the char-
acteristic layer thickness is qualitatively independent of the layer profile, and
supports the foregoing interpretation based on nano-scale Rouse dynamics.
120
5.5 Summary and conclusions
We proposed a simple hydrodynamic model to interpret the reduced
bulk viscosity of polymer-nanocomposite melts observed by Mackay et al.
(2003) and Tuteja et al. (2005). The model adopts a uniform viscosity and
density layer around nanoparticles, and can be easily extended to continu-
ous layer profiles (as demonstrated by the results in table 5–2). Comparing
the theory with experiments suggests that the reduced bulk viscosity arises
from nano-scale Rouse dynamics experienced by nanoparticles in entanglement
tubes (Brochard-Wyart & de Gennes 2000) when the nanoparticle-polymer
interactions are weak. Our model also confirms that polymer-mediated inter-
actions are indeed crucial for the large bulk viscosity reductions found with
high molecular weight polymer, as first pointed out by Tuteja et al. (2005).
The comparison also suggests that the polymers in samples that exhibit a large
bulk viscosity reduction are strongly perturbed. We further demonstrated that
a continuous variation of polymer density and viscosity around nanoparticles
does not change the qualitative picture emerging from our analytical expres-
sion for a single uniform layer. Our model showed how the microstructure
affects the bulk viscosity in a manner that is consistent with experiments, and
it highlighted how sensitive the bulk viscosity is to interactions between the
nanoparticles and polymer.
121
CHAPTER 6Summary and conclusions
This thesis showed (i) the influences of the continuous phase on the
single-particle and bulk-composite responses by investigating the electric-field-
induced steady and dynamic responses of charged colloidal particles embedded
in uncharged hydrogel matrices and (ii) the influences of colloid-induced mi-
crostructural changes on composite transport properties by theoretically inter-
preting the recently discovered anomalous bulk viscosity reduction in polymer-
nanocomposite melts. The objectives outlined in § 1.2 have been accomplished,
and the theories and results presented in chapters 3–5 contribute to our un-
derstanding of the dynamics of interacting dispersed and continuous phases in
colloid composites. These dynamics are crucial for developing new composite
materials and diagnostic techniques to probe their microstructures.
Chapters 3 and 4 investigated the steady and dynamic electric-field-induced
responses of a charged colloidal particle in an uncharged hydrogel matrix.
These theoretical studies are not only important for understanding the in-
fluences of complex fluids, but also serve as rigorous foundations for electric
microrheology and electroacoustic techniques. A multi-phase electrokinetic
model, generalized from the standard electrokinetic model, was used to de-
scribe the polymer displacement, fluid velocity, ion fluxes and electrostatic
potentials. Computer programs based on MPEK (Hill et al. 2003a) were de-
veloped to solve the multi-phase electrokinetic model accurately with steady
and dynamic electrical forcing. Generally, the electric-field-induced particle
responses depend on both the hydrogel viscoelastic and the particle physico-
chemical properties, making the electric-field-based characterizations distinct
from existing techniques. Chapter 3 showed that the particle displacements
in compressible hydrogels with Poisson’s ratio close to zero can be an order of
122
magnitude larger than those for incompressible hydrogels. This is confirmed
by dynamic studies in chapter 4, where the response spectra for compress-
ible hydrogels present an evident transition from quasi-steady compressible
to quasi-steady incompressible elastic regimes. These results imply that elec-
tric microrheology could be successfully applied to hydrogel-colloid composites
at frequencies below the reciprocal draining time. On the other hand, high
frequency responses can be measured using electroacoustics, and chapter 4
also showed that the electroacoustic signals and the single particle response
are connected by the theory of O’Brien (1988, 1990). Similarly to other elec-
trokinetic phenomena (Hunter 2001; Lyklema 1995), both steady and dynamic
responses exhibit double layer polarization and relaxation at large ζ-potentials.
Moreover, boundary-layer analytical approximations—valuable for interpret-
ing experiments—were derived in chapters 3 and 4 for both the steady and
dynamic forcing.
After studying the influences of a hydrogel matrix on particle responses
in chapters 3 and 4, chapter 5 focused on the effect of colloid-induced mi-
crostructural changes on the bulk viscosity of polymer-nanocomposite melts.
This chapter developed a hydrodynamic model that incorporates a polymer
layer with properties different from the bulk. The model suggested that the
bulk viscosity reductions observed in the experiments of Mackay et al. (2003)
and Tuteja et al. (2005) can be attributed to the Rouse dynamics in poly-
mer entanglement tubes experienced by the nanoparticles, as proposed by
Brochard-Wyart & de Gennes (2000). The model also revealed that the inter-
actions are strong in polymer-nanocomposite melts, and that polymer chains
are likely everywhere disturbed, as evidenced from the large fitted layer thick-
ness and free volume changes. This chapter shows the effect of the inhomo-
geneous layer as a result of colloid-induced microstructural changes, and is an
important first step toward more quantitative theories to interpret intriguing
nano-scale phenomena, such as bulk viscosity reduction.
123
Experimental verifications of the theoretical treatment presented in chap-
ters 3 and 4 using electric microrheology and electroacoustics are recommended
for future investigations. Other possibilities include studying the dielectric re-
laxation spectra of dilute hydrogel-colloid composites theoretically and exper-
imentally, and developing new theories that take particle-particle interactions
into account. For the colloid-induced bulk viscosity reduction, theoretical,
simulation, and experimental studies focused on resolving the polymer-particle
and particle-particle interactions that influence these anomalous behaviors are
recommended.
124
REFERENCES
Abetz, V., Brinkmann, T., Dijkstra, M., Ebert, K., Fritsch, D.,Ohlrogge, K., Paul, D., Peinemann, K. V., Nunes, S. P.,Scharnagl, N. & Schossig, M. 2006 Developments in membrane re-search: from material via process design to industrial application. Ad-vanced Engineering Materials 8, 328–358.
Ahualli, S., Delgado, A. V., Miklavcic, S. J. & White, L. R. 2006Dynamic electrophoretic mobility of concentrated dispersions of sphericalcolloidal particle. On the consistent use of the cell model. Langmuir 22,7041–7051.
Allersma, M. W., Gittes, F., deCastro, M. J., Stewart, R. J. &Schmidt, C. F. 1998 Two-dimensional tracking of ncd motility by backfocal plane interferometry. Biophysical Journal 74, 1074–1085.
Anderson, E., Bai, Z., Bischof, C., Blackford, S., Demmel, J.,Dongarra, J., Du Croz, J., Greenbaum, A., Hammarling, S.,McKenney, A. & Sorensen, D. 1999 LAPACK Users’ Guide, 3rdedn. Society for Industrial and Applied Mathematics.
Ascher, U. M., Mattheij, R. M. M. & Russell, R. D. 1988 Numeri-cal Solution of Boundary Value Problems for Ordinary Differential Equa-tions . Englewood Cliffs: Prentice Hall.
Barber, J. R. 2003 Elasticity , 2nd edn. Dordrecht: Springer.
Barndl, F., Sommer, F. & Goepferich, A. 2007 Rational design of hy-drogels for tissue engineering: impact of physical factors on cell behavior.Biomaterials 28 (2), 134–146.
Barrire, B. & Leibler, L. 2003 Kinetics of solvent absorption and perme-ation through a highly swellable elastomeric network. Journal of PolymerScience Part B: Polymer Physics 41, 166–182.
Batchelor, G. K. 1967 An Introduction to Fluid Dynamics . Cambridge:Cambridge University Press.
Baumgaertel, M., Schausberger, A. & Winter, H. H. 1990 The relax-ation of polymers with linear flexible chains of uniform length. RheologicaActa 29, 400–408.
Biot, M. A. 1941 General theory of three-dimensional consolidation. Journalof Applied Physics 12, 155–164.
Biot, M. A. 1956a Theory of propagation of elastic waves in a fluid-saturated
125
porous solid. I. low-frequency range. The Journal of The Acoustical Soci-ety of America 28, 168–178.
Biot, M. A. 1956b Theory of propagation of elastic waves in a fluid-saturatedporous solid. II. higher frequency range. The Journal of The AcousticalSociety of America 28, 179–191.
Booth, F. & Enderby, J. A. 1952 On electrical effects due to sound wavesin colloidal suspensions. Proceedings of the Physical Society. Section A65, 321–324.
Boudou, T., Ohayon, J., Arntz, Y., Finet, G., Picart, C. & Trac-qui, P. 2006 An extended modeling of the micropipette aspiration exper-iment for the characterization of the Young’s modulus and Poisson’s ratioof adherent thin biological samples: Numerical and experimental studies.Journal of Biomechanics 39, 1677–1685.
Brennen, C. 1975 A concentrated suspension model for the Couette rheologyof blood. Canadian Journal of Chemical Engineering 53, 126–133.
Breuer, K., ed. 2005 Microscale Diagnostic Techniques . Berlin: Springer.
Brinkman, H. C. 1947 A calculation of the viscous force exerted by a flowingfluid on a dense swarm of particles. Applied Scientific Research Section A1, 27–34.
Brochard, F. & de Gennes, P.-G. 1977 Dynamical scaling for polymersin theta solvents. Macromolecules 10, 1157–1161.
Brochard, F. & de Gennes, P.-G. 1992 Shear dependent slippage at apolymer/solid interface. Langmuir 8, 3033–3037.
Brochard-Wyart, F. & de Gennes, P.-G. 2000 Viscosity at small scalesin polymer melts. European Physical Journal E 1, 93–97.
Buscall, R. & Ettelaie, R. 2006 Colloidal dispersions in polymer melts.Industrial and Engineering Chemistry Research 45, 6915–6922.
Candau, S., Bastide, J. & Delsanti, M. 1982 Structural, elastic, anddynamic properties of swollen polymer networks. Advances in PolymerScience 44, 27–71.
Cash, J. R. & Mazzia, F. 2006 Hybrid mesh selection algorithms based onconditioning for two-point boundary value problems. Journal of Numeri-cal Analysis, Industrial and Applied Mathematics 1, 81–90.
Chaterji, S., Kwon, I. K. & Park, K. 2007 Smart polymeric gels: Re-defining the limits of biomedical devices. Progress in Polymer Science 32,1083–1122.
Chung, Y.-I., Ahn, K.-M., Jeon, S.-H., Lee, S.-Y., Lee, J.-H. &Tae, G. 2007 Enhanced bone regeneration with BMP-2 loaded functionnanoparticle-hydrogel complex. Journal of Controlled Release 121, 91–99.
126
Cicuta, P. & Donald, A. M. 2007 Microrheology: a review of the methodand applications. Soft Matter 3, 1449–1455.
Cohen, M. H. & Turnbull, D. 1959 Molecular transport in liquids andglasses. Journal of Chemical Physics 31, 1164–1169.
Colby, R. H., Rubinstein, M. & Daoud, M. 1994 Hydrodynamics ofpolymer solutions via two-parameter scaling. Journal de Physique II 4,1299–1310.
Cotton, J. P., Decker, D., Benoit, H., Farnoux, B., Higgins, J.,Jannink, G., Ober, R., Picot, C. & Cloizeau, J. D. 1974 Confor-mation of polymer chain in the bulk. Macromolecules 7, 863–872.
Coussy, O. 2004 Poromechanics . Chichester: Wiley.
Daoulas, K. C., Theodorou, D. N., Harmandaris, V. A., Karayian-nis, N. C. & Mavrantzas, V. G. 2005 Self-consistent-field study ofcompressible semiflexible melts adsorbed on a solid substrate and com-parison with atomistic simulations. Macromolecules 38, 7134–7149.
Dayton, P. A. & Ferrara, K. W. 2002 Targeted imaging using ultra-sound. Journal of Magnetic Resonance Imaging 16, 362–377.
Dealy, D. M. & Larson, R. G. 2006 Structure and Rheology of MoltenPolymers. From Structure to Flow Behavior and Back Again. Munich:Hanser.
Debye, P. 1933 A method for the determination of the mass of electrolyticions. Journal of Chemical Physics 1, 13–16.
DeLacey, E. H. B. & White, L. R. 1981 Dielectric response and conduc-tivity of dilute suspensions of colloidal particles. Journal of the ChemicalSociety, Faraday Transactions II 77, 2007–2039.
Doi, M. & Edwards, S. F. 1987 The Theory of Polymer Dynamics . NewYork: Oxford University Press.
Drury, J. L. & Mooney, D. J. 2003 Hydrogels for tissue engineering:scaffold design variables and applications. Biomaterials 24, 4337–4351.
Dukhin, A. S. & Goetz, P. J. 2002 Ultrasound for Characterizing ColloidsParticle Sizing, Zeta Potential, Rheology . Amsterdam: Elsevier.
Dukhin, A. S., Goetz, P. J., Wines, T. H. & Somasundaran, P.2000 Acoustic and electroacoustic spectroscopy. Colloids and Surfaces A:Physicochemical and Engineering Aspects 173, 127–158.
Dukhin, A. S., Ohshima, H., Shilov, V. N. & Goetz, P. J. 1999aElectroacoustics for concentrated dispersions. Langmuir 15, 3445–3451.
Dukhin, A. S., Shilov, V. N., Ohshima, H. & Goetz, P. J. 1999b
127
Electroacoustic phenomena in concentrated dispersions: new theory andCVI experiment. Langmuir 15, 6992–6706.
Eddington, D. T. & Beebe, D. J. 2004 Flow control with hydrogels.Advanced Drug Delivery Reviews 56, 199–210.
Einstein, A. 1906 The theory of the Brownian motion. Annalen Der Physik19, 371–381.
Enderby, J. A. 1951 On electrical effects due to sound wave in colloidalsuspensions. Proceedings of the Royal Society of London Series A 207,329–342.
Enge, A., Pelissier, P. & Zimmermann, P. 2007 MPC: Multiple Preci-sion Complex Library . INRIA.
Engler, A., Bacakova, L., Newman, C., Hategan, A., Griffin, M.& Discher, D. 2004 Substrate compliance versus ligand density in cellon gel responses. Biophysical Journal 86, 617–628.
Ferry, J. D. 1980 Viscoelastic Properties of Polymers . New York: JohnWiley and Sons.
Fleer, G. J., Cohen-Stuart, M. A., Scheutjens, J. M. H. M., Cos-grove, T. & Vincent, B. 1993 Polymers at Interfaces . London: Chap-man and Hall.
Fousse, L., Hanrot, G., Lefevre, V., Pelissier, P. & Zimmermann,P. 2007 MPFR: A multiple-precision binary floating-point library withcorrect rounding. ACM Transactions on Mathematical Software 33, 13.
Freeman, P. M., Natarajan, R. N., Kimura, J. H. & Andriacchi,T. P. 1994 Chondrocyte cells respond mechanically to compressive loads.Journal of Orthopaedic Research 12, 311–320.
Frenkel, J. 1944 On the theory of seismic and seismoelectric phenomena ina moist soil. Journal of Physics-USSR 3, 230–241.
Fung, Y. C. & Tong, P. 2001 Classical and Computational Solid Mechan-ics . Singapore: World Scientific.
Galneder, R., Kahl, V., Arbuzova, A., Rebecchi, M., Radler,J. O. & McLaughlin, S. 2001 Microelectrophoresis of a bilayer-coatedsilica bead in an optical trap: Application to enzymology. BiophysicalJournal 80, 2298–2309.
Ganesan, V., Pryamitsyn, V., Surve, M. & Narayanan, B. 2006Noncontinuum effects in nanoparticle dynamics in polymers. Journal ofChemical Physics 124 (22), 221102.
Geissler, E. & Hecht, A. M. 1980 The Poisson ratio in polymer gels.Macromolecules 13, 1276–1280.
128
Geissler, E. & Hecht, A. M. 1981 The Poisson ratio in polymer gels. 2.Macromolecules 14, 185–188.
de Gennes, P.-G. 1976a Dynamics of entangled polymer solutions: I. theRouse model. Macromolecules 9, 587–593.
de Gennes, P.-G. 1976b Dynamics of entangled polymer solutions. II. in-clusion of hydrodynamic interactions. Macromolecules 9, 594–598.
de Gennes, P.-G. 1979 Scaling Concepts in Polymer Physics . Cornell Uni-versity Press.
Gibb, S. E. & Hunter, R. J. 2000 Dynamic mobility of colloidal particleswith thick double layers. Journal of Colloid and Interface Science 224,99–111.
Granlund, T. 2007 GNU MP: The GNU Multiple Precision Arithmetic Li-brary , 4th edn.
Green, P. G. & Kramer, E. J. 1986 Matrix effects on the diffusion of longpolymer chains. Macromolecules 19, 1108–1114.
Gross, J., Scherer, G. W., Alviso, C. T. & Pekala, R. W. 1997Elastic properties of crosslinked resorcinol-formaldehyde gels and aerogels.Journal of Non-Crystalline Solids 211, 132–142.
Haraguchi, K. & Takehisa, T. 2002 Nanocomposite hydrogels: a uniqueorganic-inorganic network structure with extraordinary mechanical, op-tical, and swelling/de-swelling properties. Advanced Materials 14, 1120–1124.
Haraguchi, K., Takehisa, T. & Fan, S. 2002 Effects of clay con-tent on the properties of nanocomposite hydrogels composed of poly(N-isopropylacrylamide) and clay. Macromolecules 35, 10162–10171.
Hill, R. J. 2006a Diffusive permeability and selectivity of nanocompositemembranes. Industrial and Engineering Chemistry Research 45, 6890–6898.
Hill, R. J. 2006b Electric-field-induced force on a charged spherical colloidembedded in an electrolyte-saturated Brinkman medium. Physics of Flu-ids 18, 043103.
Hill, R. J. 2006c Reverse-selective diffusion in nanocomposite membranes.Physical Review Letters 96, 216001.
Hill, R. J. 2006d Transport in polymer-gel composites: theoretical method-ology and response to an electric field. Journal of Fluid Mechanics 551,405–433.
Hill, R. J. 2007 Electric-field-enhanced transport in polyacrylamide hydrogelnanocomposites. Journal of Colloid and Interface Science 316, 635–644.
129
Hill, R. J. & Ostoja-Starzewski, M. 2008 Electric-field-induced dis-placement of a charged spherical colloid embedded in an elastic Brinkmanmedium. Physical Review E 77, 011404.
Hill, R. J. & Saville, D. A. 2005 ‘Exact’ solutions of the full electrokineticmodel for soft spherical colloids: electrophoretic mobility. Colloids andSurfaces A: Physicochemical and Engineering Aspects 267, 31–49.
Hill, R. J., Saville, D. A. & Russel, W. B. 2003a Electrophoresisof spherical polymer-coated colloidal particles. Journal of Colloid andInterface Science 258, 56–74.
Hill, R. J., Saville, D. A. & Russel, W. B. 2003b High-frequencydielectric relaxation of spherical colloidal particles. Physical ChemistryChemical Physics 5, 911–915.
Hoffman, B. D., Massiera, G., Van Citters, K. M. & Crocker,J. C. 2006 The consensus mechanics of cultured mammalian cells. Pro-ceedings of the National Academy of Sciences of the United States ofAmerica 103, 10259–10264.
Hunter, R. J. 1998 Recent developments in the electroacoustic character-isation of colloidal suspensions and emulsions. Colloids and Surfaces A:Physicochemical and Engineering Aspects 141, 37–65.
Hunter, R. J. 2001 Foundations of Colloid Science, 2nd edn. Oxford: OxfordUniversity Press.
Hunter, R. J. & O’Brien, R. W. 1997 Electroacoustic characterization ofcolloids with unusual particle properties. Colloids and Surfaces A: Physic-ochemical and Engineering Aspects 126, 123–128.
Johnson, B., Bauer, J. M., Niedermaier, D. J., Crone, W. C. &Beebe, D. J. 2004a Experimental techniques for mechanical character-ization of hydrogels at the microscale. Experimental Mechanics 44 (1),21–28.
Johnson, B. D., Beebe, D. J. & Crone, W. C. 2004b Effects of swellingon the mechanical properties of a pH-sensitive hydrogel for use in mi-crofluidic devices. Materials Science and Engineering: C 24, 575–581.
Kang, S. W., Kim, J. H., Char, K., Won, J. & Kang, Y. S. 2006Nanocomposite silver polymer electrolytes as facilitated olefin transportmembranes. Journal of Membrane Science 285, 102–107.
Khademhosseini, A & Langer, R. 2007 Microengineered hydrogels fortissue engineering. Biomaterials 28, 5087–5092.
Kickelbick, G. 2003 Concepts for incorporation of inorganic building blocksinto organic polymers on a nanoscale. Progress in Polymer Science 28,83–114.
130
Kim, J. J. & Park, K. 1998 Smart hydrogels for bioseparation. Bioseparation7, 177–184.
Kim, S. & Karrila, S. J. 1991 Microhydrodynamics: Principles and SelectedApplications . Boston: Butterworth-Heinemann.
Kimura, Y. & Mizuno, D. 2007 Microrheology of a swollen lyotropic lamel-lar phase. Molecular Crystals and Liquid Crystals 478, 3–13.
Kopesky, E. T., Haddad, T. S., Cohen, R. E. & McKinley, G. H.2004 Thermomechanical properties of poly(methly methacrylate)s con-taining tethered and untethered polyhedral oligomeric silsesquioxanes.Macromolecules 37, 8992–9004.
Lamb, H. 1945 Hydrodynamics . New York: Dover Publications.
Landau, L. D. & Lifshitz, E. M. 1986 Theory of Elasticity , 3rd edn.Oxford: Pergamon Press.
Landau, L. D. & Lifshitz, E. M. 1987 Fluid Mechanics . Oxford: Perga-mon Press.
Larson, R. G. 1999 The Structure and Rheology of Complex Fluids . NewYork: Oxford University Press.
Levine, A. J. & Lubensky, T. C. 2000 One- and two-particle microrheol-ogy. Physical Review Letters 85, 1774–1777.
Levine, A. J. & Lubensky, T. C. 2001 Response function of a sphere in aviscoelastic two-fluid medium. Physical Review E 63, 041510.
Li, Y., Hu, Z. & Li, C. 1993 New method for measuring Poisson ratio inpolymer gels. Journal of Applied Polymer Science 50, 1107–1111.
Lin, C.-C. & Netters, A. T. 2006 Hydrogels in controlled release for-mulations: Network design and mathematical modeling. Advanced DrugDelivery Reviews 58, 1379–1408.
Lin, D. C., Yurke, B. & Langrana, N. A. 2004 Mechanical proper-ties of a reversible, DNA-crosslinked polyacrylamide hydrogel. Journal ofBiomechanical Engineering-Transactions of ASME 26, 104–110.
Lin, D. C., Yurke, B. & Langrana, N. A. 2005 Use of rigid spher-ical inclusions in Young’s moduli determination: application to DNA-crosslinked gels. Journal of Biomechanical Engineering-Transactions ofASME 127, 571–579.
Liu, J., Levine, A. L., Mattoon, J. S., Yamaguchi, M., Lee, R. J.,Pan, X. L. & Rosol, T. J. 2006 Nanoparticles as image enhancingagents for ultrasonography. Physics in Medicine and Biology 51, 2179–2189.
131
Loo, C., Lowery, A., Halas, N., West, J. & Drezek, R. 2005 Im-munotargeted nanoshells for integrated cancer imaging and therapy. NanoLetters 5, 709–711.
Lyklema, J. 1995 Fundamentals of Interface and Colloid Science. II. Solid-Liquid Interfaces . London: Academic Press.
Mackay, M. E., Dao, T. T., Tuteja, A., Ho, D. L., Van Horn,B., Kim, H.-C. & Hawker, C. J. 2003 Nanoscale effects leading tonon-Einstein-like decrease in viscosity. Nature Materials 2, 762–766.
MacKintosh, F. C. & Schmidt, C. F. 1999 Microrheology. Current Opin-ion in Colloid and Interface Science 4, 300–307.
MacRobert, T. M. 1967 Spherical Harmonics: An elementary treatise onharmonic functions, with applications., 3rd edn. Oxford: Pergamon Press.
Mangelsdorf, C. S. & White, L. R. 1992 Electrophoretic mobility of aspherical colloidal particle in an oscillating electric field. Journal of theChemical Society, Faraday Transactions 88, 3567–3581.
Mangelsdorf, C. S. & White, L. R. 1997 Dielectric response of a dilutesuspension of spherical colloidal particles to an oscillating electric field.Journal of the Chemical Society, Faraday Transactions 93, 3145–3154.
Markov, M. G. 2005 Propagation of longitudinal elastic waves in a fluid-saturated porous medium with spherical inclusions. Acoustical Physics51, S115–S121.
Masliyah, J. H. & Bhattacharjee, S. 2006 Electrokinetic and ColloidTransport Phenomena. Hoboken: Wiley-Interscience.
Mason, T. G. & Weitz, D. A. 1995 Optical measurements of frequency-dependent linear viscoelastic moduli of complex fluids. Physical ReviewLetters 74 (7), 1250–1253.
Matos, M. A., White, L. R. & Tilton, R. D. 2006 Electroosmoticallyenhanced mass transfer through polyacrylamide gels. Journal of Colloidand Interface Science 300, 429–436.
Matos, M. A., White, L. R. & Tilton, R. D. 2008 Enhanced mixing inpolyacrylamide gels containing embedded silica nanoparticles as internalelectroosmotic pumps. Colloids and Surfaces B: Biointerfaces 61, 262–269.
Maxwell, J. C. 1873 Treatise on Electricity and Magnetism, vol. 1 . London:Oxford University Press.
Merkel, T. C., Freeman, B. D., Spontak, R. J., He, Z., Pinnau,I., Meakin, P. & Hill, A. J. 2002 Ultrapermeable, reverse-selectivenanocomposite membranes. Science 296, 519–522.
Migler, K. B., Hervet, H. & Leger, L. 1993 Slip transition of a polymermelt under shear stress. Physical Review Letters 70, 287–290.
132
Milner, S. T. 1993 Dynamical theory of concentration fluctuations in poly-mer solutions under shear. Physical Review E 48, 3674–3691.
Mizuno, D., Kimura, Y. & Hayakawa, R. 2000 Dynamic electrophoreticmobility of colloidal particles measured by the newly developed methodof quasi-elastic light scattering in a sinusoidal electric field. Langmuir 16,9547–9554.
Mizuno, D., Kimura, Y. & Hayakawa, R. 2001 Electrophoretic microrhe-ology in a dilute lamellar phase of a nonionic surfactant. Physical ReviewLetters 87, 088104.
Mizuno, D., Kimura, Y. & Hayakawa, R. 2004 Electrophoretic microrhe-ology of a dilute lamellar phase: Relaxation mechanisms in frequency-dependent mobility of nanometer-sized particles between soft membranes.Physical Review E 70, 011509.
O’Brien, R. W. 1979 A method for the calculation of the effective trans-port properties of suspensions of interacting particles. Journal of FluidMechanics 91, 17–39.
O’Brien, R. W. 1986 The high frequency dielectric dispersion of a colloid.Journal of Colloid and Interface Science 113, 81–93.
O’Brien, R. W. 1988 Electro-acoustic effects in a dilute suspension of spher-ical particles. Journal of Fluid Mechanics 190, 71–86.
O’Brien, R. W. 1990 The electroacoustic equations for a colloidal suspen-sion. Journal of Fluid Mechanics 212, 81–93.
O’Brien, R. W., Jones, A. & Rowlands, W. N. 2003 A new formulafor the dynamic mobility in a concentrated colloid. Colloids and SurfacesA: Physicochemical and Engineering Aspects 218, 89–101.
O’Brien, R. W. & White, L. R. 1978 Electrophoretic mobility of a spher-ical colloidal particle. Journal of the Chemical Society, Faraday Transac-tions II 74, 1607–1626.
Oestreicher, H. L. 1951 Field and impedance of an oscillating sphere ina viscoelastic medium with an application to biophysics. The Journal ofThe Acoustical Society of America 23, 707–714.
O’Konski, C. T. 1960 Electric properties of macromolecules V. Theory ofionic polarization in polyelectrolytes. Journal of Physical Chemistry 64,605–619.
Overbeek, J. Th. G. 1943 Theorie der elektrophorese. Kolloid-Beih 54,287–364.
Pal, R. 2007 Rheology of Particulate Dispersions and Composites . Boca Ra-ton: CRC Press.
Peppas, N. A., Bures, P., Leobandung, W. & Ichikawa, H. 2000
133
Hydrogels in pharmaceutical formulations. European Journal of Pharma-ceutics and Biopharmaceutics 50, 27–46.
Peppas, N. A., Hilt, J. Z., Khademhosseini, A. & Langer, R. 2006Hydrogels in biology and medicine: from molecular principles to bionan-otechnology. Advanced Materials 18, 1345–1360.
Pine, D. J., Weitz, D. A., Chaikin, P. M. & Herbolzheimer, E. 1988Diffusing wave spectroscopy. Physical Review Letters 61, 1134–1137.
Pozrikidis, C. 1996 Introduction to Theoretical and Computational FluidDynamics . New York: Oxford University Press.
Preston, M. A., Kornbrekke, R. & White, L. R. 2005 Determinationof the dynamic electrophoretic mobility of a spherical colloidal particlethrough a novel numerical solution of the electrokinetic equations. Lang-muir 21 (22), 9832–9842.
Qiu, Y. & Park, K. 2001 Environment-sensitive hydrogels for drug delivery.Advanced Drug Delivery Reviews 53, 321–339.
Rider, P. F. & O’Brien, R. W. 1993 The dynamic mobility of particlesin a non-dilute suspension. Journal of Fluid Mechanics 257, 607–636.
Roberts, C., Cosgrove, T., Schmidt, R. G. & Gordon, G. V. 2001Diffusion of poly(dimethyl siloxane) mixtures with silicate nanoparticles.Macromolecules 34, 538–543.
Rotello, V. M., ed. 2004 Nanoparticles: Building Blocks for Nanotechnol-ogy . New York: Kluwer Academic/Plenum Publishers.
Russel, W. B., Schowalter, W. R. & Saville, D. A. 1989 ColloidalDispersions . Cambridge: Cambridge University Press.
Rutgers, Ir. R. 1962 Relative viscosity of suspensions of rigid spheres inNewtonian liquids. Rheologica Acta 2, 202–210.
Scherer, G. W. 1992 Bending of gel beams: method for characterizing elas-tic properties and permeability. Journal of Non-Crystalline Solids 142,18–35.
Schmidt, G. & Malwitz, M. M. 2003 Properties of polymer-nanoparticlecomposites. Current Opinion in Colloid and Interface Science 8, 103–108.
Schnurr, B., Gittes, F., MacKintosh, F. C. & Schmidt, C. F. 1997Determining microscopic viscoelasticity in flexible and semiflexible poly-mer networks from thermal fluctuations. Macromolecules 30, 7781–7792.
Sershen, S. R., Mensing, G. A., Ng, M., Halas, N. J., Beebe, D. J.& West, J. L. 2005 Independent optical control of microfluidic valvesformed from optomechanically responsive nanocomposite hydrogels. Ad-vanced Materials 17, 1366–1368.
Shkel, I. A, Tsodikov, O. V. & Record Jr, M. T. 2000 Complete
134
asymptotic solution of cylindrical and spherical Poisson-Boltzmann equa-tions at experimental salt concentrations. Journal of Physical ChemistryB 104, 5161–5170.
Snieder, R. & Page, J. 2007 Multiple scattering in evolving media. PhysicsToday May 2007, 49–55.
Soskey, P. R. & Winter, H. H. 1984 Large step shear strain experimentswith parallel-disk rotational rheometers. Journal of Rheology 28, 625–645.
Speight, J. G. 2005 Lange’s Handbook of Chemistry , 16th edn. New York:McGraw-Hill.
Strubulevych, A., Leroy, V., Scanlon, M. G. & Page, J. H. 2007Characterizing a model food gel containing bubbles and solid inclusionsusing ultrasound. Soft Matter 3, 1388–1394.
Takigawa, T., Morino, Y., Urayama, K. & Masuda, T. 1996 Poisson’sratio of polyacrylamide (PAAm) gels. Polymer Gels and Networks 4, 1–5.
Takigawa, T., Yamawaki, T., Takahashi, K. & Masuda, T. 1997Change in Young’s modulus of poly(N-isopropylacrylamide) gels by vol-ume phase transition. Polymer Gels and Networks 5, 585–589.
Tanaka, T., Hocker, L. O. & Benedek, G. B. 1973 Spectrum of lightscattered from a viscoelastic gel. Journal of Chemical Physics 59, 5151–5159.
Taylor, G. I. 1932 The viscosity of a fluid containing small drops of anotherfluid. Proceedings of the Royal Society of London Series A 138, 41–48.
Temkin, L. & Leung, C. M. 1976 On the velocity of a rigid sphere in asound wave. Journal of Sound and Vibration 49, 75–92.
Thevenot, C., Khoukh, A., Reynaud, S., Desbrieres, J. & Grassl,B. 2007 Kinetic aspects, rheological properties and mechanoelectrical ef-fects of hydrogels composed of polyacrylamide and polystyrene nanopar-ticles. Soft Matter 3, 437–447.
Torquato, S. 2002 Random Heterogeneous Materials: Microstructure andMacroscopic Properties . New York: Springer-Verlag.
Tuinier, R. & Lekkerkerker, H. N. W. 2002 Polymer density around asphere. Macromolecules 35, 3312–3313.
Tuteja, A., Duxbury, P. M. & Mackay, M. E. 2007a Multifunctionalnanocomposites with reduced viscosity. Macromolecules 40, 9427–9434.
Tuteja, A., Duxbury, P. M. & Mackay, M. E. 2008 Polymer chainswelling induced by dispersed nanoparticles. Physical Review Letters 100,077801.
Tuteja, A., Mackay, M. E., Hawker, C. J. & Van Horn, B. 2005
135
Effect of ideal, organic nanoparticles on the flow properties of linear poly-mers: Non-Einstein-like behavior. Macromolecules 38, 8000–8011.
Tuteja, A., Mackay, M. E., Narayanan, S., Asokan, S. & Wong,M. S. 2007b Breakdown of the continuum Stokes-Einstein relation fornanoparticle diffusion. Nano Letters 7, 1276–1281.
Urayama, K., Takigawa, T. & Masuda, T. 1993 Poisson ratio ofpoly(vinyl alcohol) gels. Macromolecules 26, 3092–9096.
Valentine, M. T., Dewalt, L. E. & Ou-Yang, H. D. 1996 Forces ona colloidal particle in a polymer solution: a study using optical tweezers.Journal of Physics: Condensed Matter 8, 9477–9482.
Verwey, E. J. W. & Overbeek, J. Th. G. 1948 Theory of Stability ofLyophobic Colloids . Elsevier, Amsterdam.
Wang, K. L., Burban, J. H. & Cussler, E. L. 1993 Hydrogels as sepa-ration agents. Advances in Polymer Science 110, 67–79.
Wang, M. & Hill, R. J. 2008 Electric-field-induced displacement of chargedspherical colloids in compressible hydrogels. Soft Matter 4, 1048–1058.
Westermeier, R. 2005 Electrophoresis in Practice. Weinheim: Wiley-VCH.
Wu, D. T., Fredrickson, G. H., Carton, J. P., Ajdari, A. &Leibler, L. 1995 Distribution of chain-ends at the surface of a poly-mer melt: compensation effects and surface-tension. Journal of PolymerScience Part B: Polymer Physics 33, 2373–2389.
Yamaguchi, N., Chae, B.-S., Zhang, L., Kiick, K. L. & Furst, E. M.2005 Rheological characterization of polysaccharide-poly(ethylene glycol)star copolymer hydrogels. Biomacromolecules 6 (4), 1931–1940.
Yao, S. H., Hertzog, D. E., Zeng, S. L., Mikkelsen, J. C. & Santi-ago, J. G. 2003 Porous glass electroosmotic pumps: design and experi-ments. Journal of Colloid and Interface Science 268, 143–153.
Yao, S. H. & Santiago, J. G. 2003 Porous glass electroosmotic pumps:theory. Journal of Colloid and Interface Science 268, 133–142.
Ziemann, F., Radler, J. & Sackmann, E. 1994 Local measurement ofviscoelastic moduli of entangled actin networks using an oscillating mag-netic bead micro-rheometer. Biophysical Journal 66, 2210–2216.