www.rsc.org/pccp | Physical Chemistry Chemical Physics
PAPER
Scaling form of viscosity at all length-scales in poly(ethylene glycol)
solutions studied by fluorescence correlation spectroscopy and
capillary electrophoresis
Robert Holyst,*ac Anna Bielejewska,a Je˛drzej Szymański,a Agnieszka Wilk,bd
Adam Patkowski,b Jacek Gapiński,b Andrzej Żywociński,a Tomasz Kalwarczyk,a
Ewelina Kalwarczyk,a Marcin Tabaka,a Natalia Zie˛bacza and
Stefan A. Wieczoreka
Received 27th April 2009, Accepted 9th July 2009
First published as an Advance Article on the web 3rd August 2009
DOI: 10.1039/b908386c
We measured the viscosity of poly(ethylene glycol) (PEG 6000, 12 000, 20 000) in water using
capillary electrophoresis and fluorescence correlation spectroscopy with nanoscopic probes of
different diameters (from 1.7 to 114 nm). For a probe of diameter smaller than the radius of
gyration of PEG (e.g. rhodamine B or lyzozyme) the measured nanoviscosity was orders of
magnitude smaller than the macroviscosity. For sizes equal to (or larger than) the polymer radius
of gyration, macroscopic value of viscosity was measured. A mathematical relation for macro and
nanoviscosity was found as a function of PEG radius of gyration, Rg, correlation length in
semi-dilute solution, x, and probe size, R. For R o Rg, the nanoviscosity (normalized by water
viscosity) is given by exp(b(R/x)a), and for R 4 Rg, both nano and macroviscosity follow the
same curve, exp(b(R/x)a), where a and b are two constants close to unity. This mathematical
relation was shown to equally well describe rhodamine (of size 1.7 nm) in PEG 20 000 and the
macroviscosity of PEG 8 000 000, whose radius of gyration exceeds 200 nm. Additionally, for the
smallest probes (rhodamine B and lysozyme) we have verified, using capillary electrophoresis and
fluorescence correlation spectroscopy, that the Stokes–Einstein (SE) relation holds, providing that
we use a size-dependent viscosity in the formula. The SE relation is correct even in PEG solutions
of very high viscosity (three orders of magnitude larger than that of water).
Introduction
Water-based polymer solutions appear in a wide variety of
systems and industrial processes and products (e.g. biological
cells, food production processes, paints, and personal care
products). Their viscosity can be larger by several orders of
magnitude than that of water. However, because water is
the dominant fraction in these solutions, small objects of
sub-nanometer size, d, should experience only the viscosity
of water while diffusing in the solutions. On the other hand,
large objects, much greater than the polymer size, should
experience the large macro-viscosity of the solution. It follows
immediately that the coefficient of viscosity depends on the
length-scale at which it is probed: viscosity should change
from the value for water, Zsolvent, at the nano-scale to a large
macro-viscosity, Zmacro, at the macroscale. A number of
questions are still open in this context: what is the length-scale
for which we observe a crossover from Zsolvent to Zmacro? What
is the relation between nanoviscosity determined from the
a
Institute of Physical Chemistry PAS, Kasprzaka 44/52,
01-224 Warsaw, Poland
b
Institute of Physics, A. Mickiewicz University, Umultowska 85,
61-614 Poznań, Poland
c
Cardinal Stefan Wyszyński University, WMP-SNS´, Dewajtis 5,
01-805 Warsaw, Poland
d
Laboratory for Neutron Scattering, ETH Zurich & Paul Scherrer
Institut, 5232 Villigen PSI, Switzerland
This journal is
c
the Owner Societies 2009
diffusion of nanoprobes and macroviscosity measured by
standard rheometers? Is the Stokes–Einstein relation valid
for nanoprobe diffusion in polymer solutions? Partial answers,
summarized below, have been given to these questions.1–27
The problem of nano and macroviscosity of polymer solutions
received a lot of attention in experimental and theoretical
studies. In the early fifties, Schachman et al.1 used an
ultracentrifuge to study the dependence of viscosity of DNA
solutions on size of sedimenting nanoprobes and observed that
small and large probes exhibited dramatically different
sedimentation coefficients. Later, the dependence of the
viscosity on the size of probes and concentration of polymer
was found to be a stretched exponential function.2–5 Laurent
et al.2 had investigated the sedimentation of bovine serum
albumin (BSA—radius 3.55 nm) in hyaluronic acid solution.
Rodbard and Chrambach3 applied capillary electrophoresis to
the motion of various proteins and dyes of sizes from 0.51 to
6 nm in polyacrylamide gel. A thorough account of all the
works in this direction based on capillary electrophoresis,3,6,7
sedimentation experiments in the ultracentrifuge1,2,5,8 or
diffusion9–12 was given by Odijk.13 This exponential dependence
of nanoviscosity on probe size and polymer concentration
prompted scientists to analyze the problem theoretically.13–21
Theories developed by Odijk,13 following Ogston,5,14 Philles,15–18
Cukier,19 Altenberger,20 Amsden21 and de Gennes22,23
predicted the same exponential dependence of viscosity on
Phys. Chem. Chem. Phys., 2009, 11, 9025–9032 | 9025
probe size, but differed in the physical mechanisms which led
to such dependence. A more recent theoretical approach which
follows Odijk13 is based on the nonuniform viscosity in
polymer solutions around a nanoprobe.24,25 All experiments
and theoretical calculations showed the same universal
stretched exponential dependence of viscosity on the size of
nanoprobes and concentration of polymer solutions.
An important breakthrough in the study of the relation
between nanoviscosity, probe size and polymer concentration
came in 1978 with the experiments of Langevin and Rondelez.8
They showed experimentally that, instead of two variables
(probe size and concentration of a polymer), the nanoviscosity
felt by a nanoprobe is an exponential function of only a single
variable, R/x, where R is the size of a nanoprobe and x is the
correlation length (distance between entaglement points of
polymer chains), which is dependent on the concentration, in
the semi-dilute polymer solution. This important experimental
work8 followed the suggestion of de Gennes22,23 who developed
a theoretical approach to the viscosity of polymer solutions. In
his approach, a dense polymer solution can be viewed as a
transient statistical network of mesh size (correlation length), x,
in a solvent. This size also corresponds to the size of a ‘‘blob’’
inside which all monomers belong to the same polymer chain.
De Gennes postulated that the viscosity should depend on R as
Z(R/x), as verified by Langevin and Rondelez. Moreover, he
suggested that the crossover from nanoviscosity to macroviscosity occurs when the size of a nanoprobe exceeds the blob
size. Thus, in this approach, the viscosity experienced by an
object of size R c x should have a constant value equal to the
macroscopic viscosity Zmacro, while for R o x, the viscosity
should depend on R as Z(R/x). Thus, in this model, the crossover length scale, L, is equal to x. Since x decreases with
polymer concentration, for any size, R, there is a well defined
concentration, x(x), for which a crossover occurs to the macroscopic viscosity. The latter prediction was not confirmed by
Langevin and Rondelez,8 who concluded their paper with the
following remark ‘‘Also the experimental observation that at
high polymer concentration the sedimentation coefficients do
not converge towards a constant value proportional to the
macroscopic viscosity, is still unclear’’. Interestingly, de Gennes’
prediction concerning the crossover length scale was
confirmed26 in sedimentation experiments in polyethylene–
propylene in decane (as a solvent) for calcium carbonate
core-shell nano-particles of size R = 4 nm. Clearly, these two
experiments8,26 gave conflicting results. Further studies in this
direction have been performed by Michelman–Ribeiro et al.27
for probes from 2 to 44 nm in poly(vinyl alcohol) solutions
using fluorescence correlation spectroscopy. They noted that
44 nm spheres move in a macroviscosity environment, and that
smaller probes feel nano-viscosity, but they did not correlate
their results with the blob size in PVA and, consequently, did
not obtain a clear-cut crossover between nano and macroviscosity. Finally, in recent computer simulations, Liu et al.28
found that when the size of a probe is larger than the gyration
radius of polymer chain, probes experience macroviscosity.
They also observed that the crossover is not sharp and occurs
when the radius of gyration of a polymer is between the
diameter and radius of a nanoprobe. Summarizing, they found
a crossover at a size comparable not to the correlation length, x,
9026 | Phys. Chem. Chem. Phys., 2009, 11, 9025–9032
but to the much larger radius of gyration, Rg. They also found
that the Stokes–Einstein relation is valid only for particles
larger than the radius of gyration. In view of these conflicting
results, the problem of crossover between nano and macroviscosity as well as the test of the Stokes–Einstein relation is still
an open problem, requiring further experimental tests.
The Stokes–Einstein (SE) equation belongs to the larger
class of relations known under the common name fluctuation–
dissipation relations.28 The SE relation was partially addressed
by Ogston,5 but there is still a lot of controversy as to how to
approach the SE relation at the nanoscale. Liu et al.29 and also
Tuteja et al.30 concluded that the Stokes–Einstein relation
holds only when the size of a particle is larger than the radius
of gyration of a polymer. In both cases, the diffusion of
particles was compared to the macroviscosity of a polymer
solution. Also, Michelman–Ribeiro et al.27 verified the SE
relation for large probes (44 nm) in poly(vinyl alcohol) solutions using fluorescence correlation spectroscopy. Similar
conclusions were reached by Zanten et al.31 using diffusive
wave spectroscopy. The list of papers validating the SE
relation for large probe sizes using diffusion measurements
in polymer solutions is long.32–41 A similar problem of the
relation between macroviscosity and diffusion exists in the
dynamics of concentrated colloidal solutions.42,43 Poon et al.42
using two-color dynamic light scattering (to avoid multiple
scattering) showed that collective long-time diffusion
coefficients measured for the q vector corresponding to the
maximum in the structure factor satisfies the generalized
Stokes–Einstein relation. Also, the viscosity of colloidal
systems was investigated by Brady and co-workers.44 The
question which still remains is as follows: can we still use the
SE relation at the nanoscale? In this case, a typical test of
the SE relation is not a simple comparison of diffusion with the
macroviscosity, but rather a comparison of hydrodynamic
drag exerted on a moving object (measured by e.g. capillary
electrophoresis or sedimenation experiments in the ultracentifuge) to its diffusion coefficent (measured by dynamic
light scattering at vanishing concentration or by fluorescence
correlation spectroscopy).
In this manuscript, we study the nanoviscosity and
macroviscosity of PEG solutions covering a wide range of
macroscopic viscosity from Zsolvent to Zmacro E 3000Zsolvent for
PEG 20 000 and Zmacro E 35 000Zsolvent for PEG 8 000 000. It is
worth mentioning that in all previous studies, the solution
viscosity rarely exceeded an order of magnitude beyond that of
the solvent8,26 (with the exception of gels) and here we cover
almost 5 orders of magnitude. Our studies are concentrated
on establishing a clear relation between nano and macroviscosity and the crossover length scale from nanoviscosity
to macroviscosity.
The article is organized as follows: in the next section we
will present the materials and describe methods used in
experiments. In section 3 we will show and discuss our results.
Conclusions are contained in section 4.
Materials and methods
We studied the viscosity of PEG (molecular mass 6000, 12 000,
20 000) solutions in water using fluorescence correlation
This journal is
c
the Owner Societies 2009
spectroscopy (FCS) (for diffusion), capillary electrophoresis
(CE) (for flow) and standard rheology (for macroscopic
viscosity) over a wide range of polymer concentration
(from 0 to 50% with macroscopic viscosity from Zsolvent to
Zmacro E 3000Zsolvent) and for objects ranging in size from
1.7 nm (rhodamine B) to 114 nm (silica spheres). Fluorescence
correlation spectroscopy, although an old technique, was not
used on a regular basis for the first 20–30 years after its
invention.45 It became a standard technique for diffusion
measurements only 10 years ago, after its combination with
confocal microscopy.46,47
PEG solutions
We used aqueous solutions (appropriate buffers) of PEG
(poly(ethylene glycol)). The molecular weights of PEG,
purchased from Fluka, were, on average, 6000, 12 000,
20 000. The measured polydispersity (using mass spectroscopy
and size exclusion chromatography GPC) was typically
PDI = 1.09 to 1.13 for PEG 20 000 with Mn = 16 300 for
PEG 20 000 (average mass between 16 000 and 24 000). For
comparison, we have also used a high molecular weight PEG
8 000 000.
Sample preparation
We purchased highly purified (by crystallization repeated three
times, and by dialysis) lysozyme protein from Sigma-Aldrich.
The tertiary structure of lysozyme is a prolate ellipsoid with
short and long axes of 30 and 50 Å, respectively (based on the
crystallographic structure). For fluorescence correlation
spectroscopy, the lysozyme was labeled with the TAMRA
(i.e. 5(6)-carboxy-tetramethylrhodamine) fluorescent dye with
the absorption and emission in the region of 555 and 580 nm,
respectively. We purchased TAMRA from Sigma-Aldrich.
TAMRA fluorescent dye can be attached to the amino group,
either at the N-terminus of the protein or on the side chain
of a lysine amino acid. The amino group has to be in the
non-protonated form to react with the fluorescent dye, thus
the reaction must be performed at a pH higher than the pK
value for the protonation of the amino group. The protein
solution was dialyzed to the 0.1 M phosphate buffer pH 8.4 in
order to label the lysine’s amino groups. A 10-fold molar
excess of the dye dissolved in DMSO was added to the
10 mg ml1 protein solution while slowly vortexing. After 24
h, the reaction was terminated by the addition of 10 mM TRIS
buffer, pH 7.5. The solution was then incubated for 12 h at
room temperature. To remove unbound dye molecules,
exhaustive filtration with the 20 mM phosphate buffer
(pH 7, 0.154 M NaCl) was performed. For the filtration we
used Millipore Amicon Ultra, centrifugal filter devices with a
molecular weight cut-off of 3 kD. The latter step was necessary
because the fluorescence signal originating from the free dye
could seriously diminish the signal from the dye bound to the
protein. The filtration procedure was repeated about 50 times.
Experimental solutions were made in a following way. First
we prepared different PEG solutions; PEG concentrations
were 5, 10, 15, 20, 30, 40 and 50 wt%. To avoid aggregation
of proteins, and to keep a constant pH, we prepared solutions
using a phosphate buffer (pH 7). Buffer solutions were
This journal is
c
the Owner Societies 2009
dialyzed. We added about 2 mL of protein/buffer solution to
these PEG/buffer solutions. The final protein concentration
was about 4 108 M. For capillary electrophoresis, the
charge ladder of lysozyme was prepared. The preparation of
the protein charge ladder is described in detail elsewhere.48
Capillary electrophoresis (CE) set up for determination of
hydrodynamic drag
Electrophoretic experiments were performed on a custom-build
apparatus consisting of high-voltage supply model CZE
1000 R (Spellman,. Plainview, NY, USA) and Spectra 100
variable-wavelength UV-Vis CE detector from Thermo
Separation Products (Fremont, CA, USA). The capillary
was inserted in the plexiglas box with fan. The temperature
was controlled inside the box. The fused-silica 75 mm i.d.
capillaries with total length of 55 cm and effective length of
43 cm to the detector window for lysozyme and 38–26 cm for
rhodamine were used. New capillaries were conditioned by
passing a 0.1 M NaOH solution for 20 min and then washing
with distilled water for 10 min, and finally equilibrating with
an appropriate electrolyte. The electrolyte solutions were
prepared by dissolving the appropriate amount of polyethylene glycol (PEG) 6000, 12 000 or 20 000 (concentration
range 0–50% v/v) with water for HPLC adjusted to pH = 2.9
with H3PO4 for rhodamine or PEG 20 000 with tris-glycine
buffer of pH 8.4 (25 mM tris base and 192 mM glycine)
for lysozyme. A new capillary was used for each PEG
(6000, 12 000 and 20 000) and sample. One capillary was used
for all concentrations of each PEG and sample. The capillary
was filled with new solution by applied pressure until a drop of
the solution appeared on the opposite end of the capillary. The
capillary was not washed between runs.
The sample of rhodamine was prepared by dissolving in
water adjusted to pH = 2.9 with H3PO4. The protein charge
ladder was prepared as described elsewhere48 and the reaction
was terminated by the addition of tris-glycine buffer and this
solution was introduced into the capillary. The samples were
introduced electrokinetically (depending on PEG concentration,
for 0.5 to 5 min and at 5 to 15 kV). The applied voltage for
analysis was 20 kV. During the experiment, the temperature in
the capillary box was maintained at 23 1C. The electroosmotic
flow (EOF) was determined using DMSO as the neutral
marker.
On electropherograms, for PEG solution of concentration
0–20%, peaks of rhodamine and DMSO (electroosmotic flow
marker) were visible. For higher concentrations of PEG, only
the rhodamine peak was observed. The electrophoretic mobility
of the sample was calculated according to the formula:
1
1
LD LT
mep ¼
tP tEOF
V
for PEG concentrations 0 to 20% and
mep ¼
1 LD L
tP V
for PEG concentrations 20%–42% (for concentrations higher
than 20% of PEG the input of electroosmotic mobility was not
higher than 5%), where LT, LD are the total length of the
capillary and the length of the capillary from the inlet to the
Phys. Chem. Chem. Phys., 2009, 11, 9025–9032 | 9027
detector, respectively, and V is the applied voltage during the
separation.
Fluorescence correlation spectroscopy (FCS) setup for diffusion
measurements
The FCS setup used in the experiments was a commercial
monochromatic FCS ConfoCor II (Carl Zeiss, Jena,
Germany). The experiments were conducted at 23 1C using a
543 nm He–Ne laser for illumination. The objective used was
C-Apochromat 40 /1.2 (N.A. 1.2) with the pinhole diameter
set to 78 mm (1 Airy unit for 543 nm laser). An avalanche
photo diode was used for detection. The laser intensities were
about 40 mW corresponding to the power density at the focal
point B30 kW cm2. The accessible time range for the
measurements of the autocorrelation function was 1 ms–10 s.
In the FCS experiment, the fluorescence intensity emitted
from a small volume element of solution, optically defined as
the focal volume of the confocal microscope objective, is
recorded as a function of time. The recorded fluorescence
fluctuates as the fluorescently labeled proteins diffuse in
and out of the focal volume or change their photo-physical
properties. The focal volume is of the order of 1 fl
(i.e. 1015 liter) and the concentration of the fluorescent
molecules is in the nanomolar range, making FCS an example
of a single molecule detection technique. The fluorescence
intensity fluctuation time series can be analyzed by means of
the autocorrelation function which contains the information
about the average number of fluorescent molecules in the focal
volume and their average residence time in the focal volume.
Distribution of the intensity (I) of the laser light in the focal
volume is often approximated as a three dimensional
Gaussian: I(x,y,z) = I0 exp(2(x2 + y2)/F2 2z2/P2), where
F is the cross-sectional length in the x–y plane, and P is the
height of the illuminated element of volume. We recorded the
intensity, S(t), of fluorescence emitted from this volume
as a function of time. This intensity fluctuated as single
fluorescently labeled molecules diffused in and out of the focal
volume. The signal, S(t), allowed the extraction of the
distribution of residence times (tD) of the tracers in the focal
volume by
R analyzing the autocorrelation function of S(t):
g(Dt) = dt S(t) S(t + Dt).
The autocorrelation function adopts different analytical
forms depending on the assumed physicochemical processes
taking place in the studied solution. For three dimensional
multicomponent diffusion, with the triplet state correction, the
autocorrelation function has the following form:
GðtÞ ¼ 1 þ
n
s t=t 1 X
e
1s
N i¼1
Ai
1=2 ;
1 þ o12 tDt
1 þ tD
t
i
i
ð1Þ
where s is the fraction of dye molecules in the triplet state, t is
the triplet lifetime (usually in the range of nanoseconds up to
microsecond), N is the average number of molecules in the
focal volume (usually of the order of 1, thus we register a
signal from a single molecule at a given time), tDi is
the residence time of the molecule in the focal volume,
o = P/F is the structure parameter describing the ratio
9028 | Phys. Chem. Chem. Phys., 2009, 11, 9025–9032
Fig. 1 Semi log-plot of the normalized autocorrelation function for
rhodamine 6G in PEG 20 000. PEG weight concentrations increase
from 0 to 50%. The average time that the fluorescent probe spends in
the focal volume is given by the characteristic decay time of the
autocorrelation function. For FCS, this autocorrelation function is
algebraic in time (eqn (1)) and only one mode of diffusion was present.
between longitudinal and transverse size of the focal volume,
Ai is the fraction of the ith component, n is the number of
diffusing species (e.g. a free dye and a protein-bound dye are
the two objects (n = 2) with different diffusion coefficients).
The residence time, tD, is related to the diffusion coefficient, D,
by D = F2/(4t) for the transverse direction (in the x–y plane).
The transverse radius of the focal volume, F, is taken from the
calibration measurement of the residence time, tcal, of the
fluorescent dye with a known diffusion coefficient, Dcal.
Finally, we can write the simple expression relating the
diffusion coefficient, D, of the studied molecule to the diffusion
coefficient of the fluorescent dye, Dcal, used for calibration,
together with their residence time
D ¼ Dcal
tcal
:
tD
ð2Þ
To exclude the possibility of a non-Gaussian shape of the
confocal volume, free diffusion of rhodamine 6G in water was
measured prior to each experiment and obtained data fitted
well with a single component autocorrelation function for
normal diffusion. The value of the fitted structure parameter
o = P/F for rhodamine 6G amounted to about 5. Typical
residence time for rhodamine was about 31 ms in our FCS
setup. It gives the following dimensions of the focal volume:
F = 0.186(3) mm, P = 0.932(15) mm. At least 20 correlation
functions accumulated for 30 s were recorded for each solution
and their average was analyzed with Origin in a nonlinear least
squares fitting module and by the MEMFCS program49,50
based on the maximum entropy method. The standard
deviation of the values of the autocorrelation function was
calculated as the standard deviation of the average and used
for the calculation of the w2 test. Fig. 1 shows typical
autocorrelation functions for rhodamine in different PEG
20 000 solutions. Details concerning the algebraic decay of
the autocorrelation function are contained elsewhere.45–47,49,50
Macroviscosity in PEG solutions
The shear viscosity was measured using a standard
TA Instruments AR2000N. The macroscopic values of the
viscosity for PEG 6000, 12 000, 20 000 and 8 000 000 (8M)
solutions are shown in Fig. 2.
This journal is
c
the Owner Societies 2009
Fig. 2 Macro-viscosity is shown for different molecular weights
(PEG 6000, 12 000, 20 000, 8 000 000) versus polymer weight fraction.
that xlimit = 0.18 g cm3 (15% w/w) for rhodamine B. Above
this concentration of PEG, rhodamine B should diffuse or flow
in all PEG–water solutions and experience the
macroscopic viscosity of the solution. Fig. 4 contradicts this
theory. Rhodamine B diffuses in PEG 6000, 12 000 and 20 000
experiencing a viscosity which is much smaller than the
macroscopic viscosity (at high concentrations of PEG 20 000,
the nanoviscosity ‘‘felt’’ by rhodamine is almost 100 smaller
than the macro-viscosity of the solution). Moreover, this
nanoviscosity weakly depends on the molecular mass
of PEG. The viscosity was determined from the diffusion
coefficient, D, measured in FCS and from the electrophoretic
mobility, m, measured in CE. The diffusion coefficient, D, was
normalized by its value for pure water, D0, similarly to the
electrophoretic mobility. The analysis of CE and FCS
Using four different polymers, we covered almost 5 orders
of magnitude of solution viscosity and two orders of
magnitude in concentration.
Results and discussion
PEG (or PEO) is a flexible polymer and its radius of gyration
in water, as a function of molecular mass, is given51 by
Rg = 0.02 M0.58
[nm]. We find Rg = 3.1 nm (for PEG
p
6000), Rg = 4.6 nm (for PEG 12 000) and Rg = 6.2 nm (for
PEG 20 000). The size of a ‘‘blob’’ inside which all monomers
belong to the same polymer chain, is a function of polymer
concentration, x, x = Rg(x/x*)0.75, where x* is the polymer
concentration at which polymer chains start to overlap. The
overlap concentration is given by x* = Mp/(4/3 pRg3 NA),
where Mp is the molar mass of polymer and NA is Avogadro’s
number.
The overlap concentration depends on the molecular mass
of PEG: x* = 0.08 g cm3 (for PEG 6000), x* = 0.05 g cm3
(for PEG 12 000), x* = 0.03 g cm3 (for PEG 20 000). Thus
for a 3% solution of PEG 20 000 in water, the polymer chains
start to overlap. The ‘‘blob size’’, x, only weakly depends on
the molecular mass (Fig. 3).
According to de Gennes’ theory,22,23 above a certain
concentration xlimit = x*(R/Rg)4/3, a small probe (e.g.
rhodamine B, R = 1.7 nm), of size comparable to the ‘‘blob’’
size, x, should experience, during its motion, the macroscopic
viscosity of the solution. For PEG 6000, 12 000, 20 000 we find
Fig. 3 Calculated (see text) blob size (correlation length) x as a
function of the polymer weight fraction for PEG 6000, 12 000,
20 000 and 8M. The blob size does not exceed the radius of gyration.
We consider concentrations larger than the overlap concentration.
This journal is
c
the Owner Societies 2009
Fig. 4 Nano-viscosity, normalized by the water viscosity, as a
function of the polymer weight fraction, determined from the flow
of rhodamine B under the influence of the electric field measured in
capillary electrophoresis (CE) (open circles) and from its diffusion
measured in fluorescence correlation spectroscopy (FCS) (filled circles)
compared to the macroviscosity (see also Fig. 2) (shown as filled
squares) (a) for solution of PEG 6000 (b) 12 000 and (c) 20 000. Please
note that rhodamine B, according to current theories, should follow
the macroviscosity line for concentrations above 15% w/w.
Phys. Chem. Chem. Phys., 2009, 11, 9025–9032 | 9029
experiments for rhodamine leads to the following equation
(Fig. 4):
D0 m0
Z
¼
¼
m
D
Z0
ð3Þ
and therefore we identified Z as the nanoviscosity (normalized
by the water viscosity, Z0). Eqn (3) and Fig. 4 give direct
experimental proof of the validity of the Stokes–Einstein
relation at the nanoscale. Thus, the coefficient of nanoviscosity,
Z, at the nanoscale can be determined using the Stokes–
Einstein formula, Z = kT/3pDR, where k is the Boltzmann
constant, T is the temperature, and R is the diameter of the
probe. In fact, for each system we have to check that the SE
relation is valid before we can use it to determine the viscosity.
One can also look at the SE relation in a slightly different
way. For a probe moving in a complex liquid there are two
important time scales: a typical time, tP, needed for the probe
to change the local structure of the fluid and the time, tR,
required for the relaxation of this perturbation. If tR { tP, the
structured fluid relaxes rapidly during the motion of the object
and the SE relation is satisfied. We stress that, in order to have
reasonable and comparable results from the flow and diffusion
(each provides a drag coefficient proportional to viscosity), the
Stokes–Einstein relation must be valid. The same is not true
for solutions of stiff polymers, semi-dilute networks of F-actin
(with a persistence length 1000 the diameter of the actin
fiber), solutions of fd virus52,53 (of size 800 nm and persistence
length 2200 nm) or in solutions of supramolecular polymers
that are hydrogen bonded (of persistence length well exceeding
100 nm54,55). In these cases, several characteristics are
observed, including: different diffusion coefficients at short
and long time scales, anomalous diffusion, and the lack of
clear crossover between the nano- and macro-scale. The
diffusing tracers do not probe the average structure of
the fluid—instead, they probe the local non-equilibrium
configurations of the complex liquid. The long-time and
short-time diffusion coefficients differ significantly52,55 from
each other and the Stokes–Einstein equation is not valid. In
such systems, visco-elasticity changes with the scale of motion
i.e. at the nanoscale we observe a viscous flow, while at the
macroscale we have a fully elastic system (e.g. like in a gel56).
In our system, satisfying the Stokes–Einstein formula, we can
use it to determine the nanoviscosity at all length scales.
Therefore, in our case, a sphere of diameter R experiences
hydrodynamic drag in the PEG solution of 3pZ(R)R, with a
size-dependent viscosity, Z(R). In Fig. 5, another test of the
Stokes–Einstein relation for a lysozyme protein moving in a
concentrated PEG 20 000 solution is given.
For the determination of the crossover length scale we used
the following probes: rhodamine B (R = 1.7 nm), lysozyme
(and lysozyme charge ladder for electrophoretic mobility48)
(R = 3.1 nm), apoferritin (R = 13.8 nm) quantum dots
(R = 25 nm) and silica spheres (R = 114 nm). The results
are summarized in Fig. 6. The rhodamine B and lysozyme
experience a viscosity very different from the macro-viscosity,
while appoferitin, quantum dots and silica spheres move
according to the macro-viscosity in PEG 20 000 solution.
All data for these nanoprobes, and also for macro-viscosity
of PEG 6000, 12 000, 20 000 and, additionally, for PEG
9030 | Phys. Chem. Chem. Phys., 2009, 11, 9025–9032
Fig. 5 Nano-viscosity, normalized by the water viscosity, as a
function of PEG 20 000 weight fraction, determined from the flow of
the lysozyme charge ladders48 measured in capillary electrophoresis
(open circles) and from diffusion measured in fluorescence correlation
spectroscopy (filled circles) in the PEG 20 000 solution with macroviscosity shown as filled squares. Fig. 4 and this figure provide an
experimental test of the Stokes–Einstein relation (see also eqn (3)).
Fig. 6 Nano-viscosity for rhodamine (circles), lysozyme (triangles),
apoferritin (squares) in PEG 20 000 solution and macroviscosity for
PEG 20 000 (filled squares). Apoferritin, silica spheres (crosses) and
quantum dots (rhombi) follow the macroviscosity. Thus, the crossover
length scale in PEG20000 is larger than 3.1 nm and smaller than
13 nm. Fig. 7 and eqn (4) and (5) leave no doubt that it is the radius of
gyration (here 6.1.nm) which is the crossover length scale from nano to
macroviscosity.
8 000 000, were collected and put on a single plot shown in
Fig. 7. From this plot we found the viscosity in the form:
a
Z
R
ð4Þ
¼ exp b
Z0
x
for the probes of size smaller than the radius of gyration i.e.
R o Rg. For R 4 Rg, the viscosity given by eqn (4) attains its
macroscopic value given by:
a
Rg
Zmacro
:
ð5Þ
¼ exp b
x
Z0
All our results (Fig. 7) follow the master curve given by eqn (4)
and (5) with a = 0.70 0.05 and b = 1.45 0.15. Fig. 3 and
eqn (4) and (5) shown in Fig. 7 combine all our results from
macrorheology and nanorheology and demonstrate (together
with Fig. 6) that the crossover length-scale, L, is indeed given
by the radius of gyration, Rg. These results are consistent with
our previous studies where we unambiguously demonstrated57
a sharp crossover from nano- to macro-viscosity at length
This journal is
c
the Owner Societies 2009
Fig. 7 The normalized plot for macroviscosity of PEG 6000, 12 000,
20 000, 8 000 000 solutions and nanoviscosity of rhodamine B and
lysozyme in PEG 20 000 solution. On this plot we present all data
using the scaling form given by eqn (4) and (5). We plot ln(Z/Z0)
(viscosity divided by the viscosity of water—see eqn (3)) versus R/x,
where x is the ‘‘blob’’ size in the polymer solution and R is the size of
the probe (for macroviscosity we use eqn (5) and instead of R we use
the radius of gyration). All data points collapse on a single master
curve showing that the crossover length scale is given by the radius of
gyration. On this plot, we present data for rhodamine of size 1.7 nm
and PEG 8 000 000 of radius of gyration 201 nm—showing an intrinsic
connection between nano and macroviscosity. Please note that for the
polymer itself Rg/x = (x/x*)0.75 is related to concentration, x, divided
by the overlap concentration, x*.
scale L E 17 nm in aqueous surfactant solutions (hexaethylene
glycol monododecyl ether , C12E6 ). This length-scale did not
depend on surfactant concentration. In this experiment,57
the concentration of surfactant varied over two orders of
magnitude (from 0.1 to 35% by weight in water, w/w) and
the diameters, R, of the probes varied over three orders of
magnitude (ranging from 0.28 nm for heavy water to 190 nm
for fluorescent polystyrene spheres). We have recently
performed SANS scattering studies and identified 17 nm in
our previous study57 with the length of semi-flexible
(practically rigid) micelles of C12E6. Our results for the
macroviscosity (eqn (5) and Fig. 3) are also consistent
with the macroviscosity measurements performed for linear
polymers over a wide range of concentrations by Takahashi
et al.58 In previous works,1–12 summarized by Odijk,13 the
following form of the nanoviscosity was proposed:
Z
¼ expðKRm xn Þ:
ð6Þ
Z0
In many experimental works,2–12 values of two exponents
m and n were determined. For different systems, they varied from
0.69 to 1.0 and 0.5 to 1.1, respectively. If we transform our
eqn (4) to the form of eqn (6), we obtain the following relation
between the parameters of K, m and n and our parameters:
0
1
3= 5 a
4
M
=
P
K ¼ b@
Rg 4 A , m = a, and n = 0.75 a.
4=3pNA
fluids,59 where the confining surfaces at the nanoscale very
strongly affect rheological behavior. It is also complementary
to the bulk rheology of polymers.58 We found a connection
between the motion of nanoscopic probes (nanoviscosity) in
polymer solutions and their macroviscosity (Fig. 7 and eqn (4)
and (5)). Our results for aqueous PEG solutions demonstrated
that the crossover length scale, L, in polymer solutions is not
related to the ‘‘blob’’ size, x, contrary to the theoretical
assumption.8,22,23 We found experimentally that the crossover
length scale, L, is given by Rg, the radius of gyration (in
accordance with recent computer simulations of Liu et al.28),
Our results explain the lack of a crossover in the experiments
of Langevin and Rondelez8 (their probes BSA and Ludox
spheres were smaller than the radius of gyration of PEG in the
experiment). Finally, we observed that de Gennes’ scaling
form, Z(R/x), is obeyed for all probes and even for macroviscosity, if we identify R with the radius of gyration, Rg. In
this approach we identified R as the diameter of a probe (and
not its radius). In fact, from computer simulations, we have
learned that the crossover is not28 sharp (see Fig. 1 of ref. 28)
but occurs for R between Rg and 2Rg. We may expect some
correction to our results which comes from the specific
structure of a polymer solution near the surface of a
nano-object.60 At the crossover length scale, i.e. R comparable
to Rg, the distribution of polymer chains near the surface is
very different from the bulk i.e. a depletion layer forms around
the object in the non-adsorbing case. Therefore, the local
viscosity is smaller near the surface than in the bulk.24,25,60
Thus, the local motion of an object in the depletion layer is
faster (larger diffusion coefficient) than the long time motion.
In our case, at the crossover, the depletion layer is, at most, of
nanometer size, whereas the size of the focal volume in FCS
confocal microscopy is around 1000 nm i.e. three orders of
magnitude larger. Therefore FCS determines the long time and
large scale motion and the depletion layer does not influence
our results. In FCS, we are not able to observe the local
motion of a protein at the nanometer scale.
The results of the present paper, together with the results of
our previous study,57 give a clear picture of the crossover
length scale in complex fluids. This length scale is related to the
size of the object which forms a complex liquid (apart from the
solvent). A consequence of this observation is the hypothesis
that a polymer chain in a polymer solution should move
experiencing the macroviscosity of the solution and a
micelle from a micellar solution should also experience the
macroviscosity of the solution during its motion. Using FCS
we can trace the motion of a single object (which forms the
complex liquid) in this complex liquid and determine the
diffusion coefficient Dmacro. Additionally, we can easily
measure the macroviscosity of such a solution, Zmacro. Using
the Stokes–Einstein formula, Dmacro = kT/3pZmacroR (or its
analog for elongated molecules) allows the determination of
the hydrodynamic diameter of the object, R.
Conclusions
Acknowledgements
Nanorheology is an emerging field of physical chemistry. Our
present studies using CE and FCS for bulk aqueous polymer
solutions are complementary to those performed for confined
This work was supported by the Project operated within the
Foundation for Polish Science Team Programe co-financed
by the EU ‘‘European Regional Development Fund’’
This journal is
c
the Owner Societies 2009
Phys. Chem. Chem. Phys., 2009, 11, 9025–9032 | 9031
TEAM/2008-2/2 and from the budget of the Ministry of Science
and Higher Education 2007–2010. R.H., E.K., T.K. and M.T.
acknowledge support from the Foundation for Polish Science
(‘‘Mistrz’’ and TEAM stipendships). J.S. acknowledges support
from the Unilever company. N.Z. acknowledges a PhD
scholarship from the President of the Polish Academy of
Sciences. This work was also supported by a research grant
from the Human Frontier Science Program Organization. We
acknowledge fruitful discussion with Piotr Garstecki.
Notes and references
1 H. K. Schachman and W. F. Harrington, J. Am. Chem. Soc., 1952,
74, 3965–39661.
2 T. C. Laurent and A. Pietruszkiewicz, Biochim. Biophys. Acta,
1961, 49, 258–264.
3 D. Rodbard and A. Chrambach, Anal. Biochem., 1971, 40, 95–134;
D. Rodbard and A. Chrambach, Science, 1971, 172, 440–451.
4 T. C. Laurent, A. Björk, A. Pietruszkiewicz and H. Perrson,
Biochim. Biophys. Acta, 1963, 78, 351–359.
5 A. G. Ogston, B. N. Preston and J. D. Wells, Proc. R. Soc. London,
Ser. A, 1973, 333, 297–316.
6 S. P. Radko and A. Chrambach, Electrophoresis, 1996, 17,
1094–1102.
7 S. P. Radko and A. Chrambach, Biopolymers, 1997, 42, 183–189.
8 D. Langevin and F. Rondelez, Polymer, 1978, 19, 875–882.
9 G. D. J. Phillies, G. S. Ullmann and K. Ullmann, J. Chem. Phys.,
1985, 82, 5242–5246.
10 M. R. Wattenberger, V. A. Bloomfield, Z. Bu and P. S. Russo,
Macromolecules, 1992, 25, 5263–5265.
11 M. Tokita, M. Toshikazu, T. Kiyonori and H. Kunio, Phys. Rev.
E., 1996, 53, 1823–1827.
12 T. F. Kossar and R. J. Phillips, AIChE J., 1995, 41, 701–711.
13 T. Odijk, Biophys. J., 2000, 79, 2314–2321.
14 A. G. Ogston, Trans. Faraday Soc., 1958, 54, 1754–1757.
15 G. D. J. Phillies, W. Brown and P. Zhou, Macromolecules, 1992,
25, 4948.
16 G. D. J. Phillies and D. Clomenil, Macromolecules, 1993, 26, 167.
17 G. D. J. Phillies, Macromolecules, 1987, 20, 558.
18 G. D. J. Phillies, J. Non-Cryst. Solids, 1991, 612, 131–133.
19 R. I. Cukier, Macromolecules, 1984, 17, 252.
20 A. R. Altenberger and M. Tirrell, J. Chem. Phys., 1984, 80, 2208.
21 B. Amsden, Polymer, 2002, 43, 1623.
22 P. G. De Gennes, Scaling concept in polymer physics, Cornell
University Press, 1979.
23 F. Brochard Wyart and P. G. de Gennes, Eur. Phys. J. E, 2000, 1,
93–97.
24 T. H. Fan, B. Xie and R. Tuinier, Phys. Rev. E, 2007, 76, 051405.
25 R. Tuinier and T. H. Fan, Soft Matter, 2008, 4, 254–257.
26 P. Tong, X. Ye, B. J. Ackerson and L. J. Fetters, Phys. Rev. Lett.,
1997, 79, 2363–2366.
27 A. Michelman-Ribeiro, F. Horkay, R. Nossal and H. Boukari,
Biomacromolecules, 2007, 8, 1595–1600.
28 J. Liu, D. Cao and L. Zhang, J. Phys. Chem. C, 2008, 112,
6653–6661.
29 E. Frey and K. Kroy, Ann. Phys., 2005, 14, 20–50.
30 A. Tuteja, M. E. Mackay, S. Narayanan, S. Asokan and
M. S. Wong, Nano Lett., 2007, 7, 1276–1281.
31 J. H. van Zanten, S. Amin and A. A. Abdala, Macromolecules,
2004, 37, 3874–3880.
32 K. A. Streletzky and G. D. J. Phillies, ACS Symp. Ser., 2000, 739,
297.
33 J. Won, C. Onyenemezu, W. G. Miller and T. P. Lodge, Macromolecules, 1994, 27, 7389.
34 K. E. Bremmell, N. Wissenden and D. E. Dunstan, Adv. Colloid
Interface Sci., 2001, 89–90, 141.
9032 | Phys. Chem. Chem. Phys., 2009, 11, 9025–9032
35 D. Gold, C. Onyenemezu and W. G. Miller, Macromolecules, 1996,
29, 5700.
36 G. D. J. Phillies, J. Gong, L. Li, L. P. Rau and J. Rollings, J. Phys.
Chem., 1989, 93, 6219.
37 W. Brown and R. Rymden, Macromolecules, 1988, 21, 840.
38 C. N. Onyenemezu, D. Gold, M. Roman and W. G. Miller,
Macromolecules, 1993, 26, 3833.
39 G. D. J. Phillies, C. Richardson, C. A. Quinlan and S. Z. Ren,
Macromolecules, 1993, 26, 6849.
40 T. H. Lin and G. D. J. Phillies, Macromolecules, 1984, 17, 1686.
41 M. Mustafa and P. S. Russo, J. Colloid Interface Sci., 1989, 129,
240.
42 W. C. K. Poon, S. P. Meeker, P. N. Pusey and P. N. Segré,
J. Non-Newtonian Fluid Mech., 1996, 67, 179–189; P. N. Segré,
S. P. Meeker, P. N. Pusey and W. C. K. Poon, Phys. Rev. Lett.,
1995, 75, 958; R. Sigel, S. Pispas, D. Vlassopoulos,
N. Hadjichristidis and G. Fytas, Phys. Rev. Lett., 1999, 83, 4666.
43 P. N. Pusey, P. N. Segré, O. P. Bekrend, S. P. Meeker and W. C.
K. Poon, Physica A, 1997, 235, 1–8.
44 D. R. Foss and J. F. Brady, J. Fluid Mech., 2000, 407, 167–200;
A. M. Leshansky, J. F. Morris and J. F. Brady, J. Fluid Mech.,
2008, 597, 305–341; J. F. Morris and J. F. Brady, J. Fluid Mech.,
1997, 348, 103–139.
45 D. Magde, E. L. Elson and W. W. Webb, Phys. Rev. Lett., 1972,
29, 705.
46 M. Eigen and R. Rigler, Proc. Natl. Acad. Sci. U. S. A., 1994, 91,
5740–5747.
47 P. Schwille, J. Korlach and W. W. Webb, Cytometry, 1999, 36,
176–182; S. Hess and W. W. Webb, Biophys. J., 2002, 83,
2300–2317.
48 J. Gao, M. Mammen and G. M. Whitesides, Science, 1996, 272,
535; J. D. Carbeck and R. S. Negin, J. Am. Chem. Soc., 2001, 123,
1252; I. Gitlin, J. D. Carbeck and G. M. Whitesides, Angew.
Chem., Int. Ed., 2006, 45, 3022; J. Szymanski, E. Pobozy,
M. Trojanowicz, A. Wilk, P. Garstecki and R. Holyst, J. Phys.
Chem. B, 2007, 111, 5503–5510.
49 P. Sengupta, K. Garai, J. Balaji, N. Periasamy and S. Maiti,
Biophys. J., 2003, 84, 1977–1984.
50 J. Szymanski, A. Patkowski, J. Gapinski, A. Wilk and R. Holyst,
J. Phys. Chem. B, 2006, 110, 7367–7373.
51 K. Devanand and J. C. Selser, Macromolecules, 1991, 24,
5943–5947.
52 I. Y. Wong, M. L. Gardel, D. R. Reichman, E. R. Weeks,
M. T. Valentine, A. R. Bausch and D. A. Weitz, Phys. Rev. Lett.,
2004, 92, 178101.
53 K. Kang, J. Gapiński, M. P. Lettinga, J. Buitenhuis, G. Meier,
M. Ratajczyk, J. K. G. Dhont and A. Patkowski, J. Chem. Phys.,
2005, 122, 044905.
54 J. Gucht, N. A. M. Besseling, W. Knoben, L. Bouteiller and
M. A. Cohen Stuart, Phys. Rev. E., 2003, 67, 051106.
55 K. Kang, A. Wilk, J. Buitenhuis, A. Patkowski and J. K. G.
Dhont, J. Chem. Phys., 2006, 124, 044907.
56 N. Fatin-Rouge, K. Starchev and J. Buffle, Biophys. J., 2004, 86,
2710–2719.
57 J. Szymański, A. Patkowski, A. Wilk, P. Garstecki and R. Holyst,
J. Phys. Chem. B, 2006, 110, 25593–25597.
58 Y. Takahashi, Y. Isono, I. Noda and M. Nagasawa, Macromolecules,
1985, 18, 1002–1008.
59 A. L. Demirel and S. Granick, Phys. Rev. Lett., 1996, 77,
2261–2264; C. Cottin-Bizonne, S. Jurine, J. Baudry, J. Crassous,
F. Restagno and E. Charlaix, European J. Phys. E, 2002, 9, 47–53;
A. Mukhopadhyay and S. Granick, Curr. Opin. Colloid Interface
Sci., 2001, 6, 423–429.
60 T. H. Fan, J. K. G. Dhont and R. Tuinier, Phys. Rev. E, 2007, 75,
011803; J. Sprakel, J. van der Gucht, M. A. Cohen Stuart and N.
A. M. Besseling, Phys. Rev. E, 2008, 77, 061502; Q. Lu and
M. J. Solomon, Phys. Rev. E, 2002, 66, 061504; M. J.
Solomon and Q. Lu, Curr. Opin. Colloid Interface Sci., 2001, 6,
430–437.
This journal is
c
the Owner Societies 2009