Structured Fluids - Polymers, Colloids, Surfactants (PDFDrive)
Structured Fluids - Polymers, Colloids, Surfactants (PDFDrive)
T. W I T T E N
The University of Chicago,
Chicago, Illinois, USA
with
P. P I N C U S
University of California,
Santa Barbara, USA
1
3
Great Clarendon Street, Oxford OX2 6DP
Oxford University Press is a department of the University of Oxford.
It furthers the University’s objective of excellence in research, scholarship,
and education by publishing worldwide in
Oxford New York
Auckland Cape Town Dar es Salaam Hong Kong Karachi
Kuala Lumpur Madrid Melbourne Mexico City Nairobi New Delhi
Shanghai Taipei Toronto
With offices in
Argentina Austria Brazil Chile Czech Republic France Greece
Guatemala Hungary Italy Japan Poland Portugal
Singapore South Korea Switzerland Thailand Turkey Ukraine Vietnam
You must not circulate this book in any other binding or cover
and you must impose this same condition on any acquirer
ISBN 978-0-19-958382-9
10 9 8 7 6 5 4 3 2 1
For Molly
This page intentionally left blank
Preface
This book is about fluids containing polyatomic structures such as polymer
molecules or colloidal grains. Such structured fluids have come to be known
as soft condensed matter. The book takes a scaling approach, following
de Gennes’s important monograph Scaling Concepts in Polymer Physics. My
purpose is to provide a unified, pedagogical introduction to soft-matter phenom-
ena that embodies this scaling approach. The book focuses on how to account
in the simplest way for the distinctive length, time, and energy scales that char-
acterize each phenomenon. This approach allows unity and simplicity, but it
provides only an initial glimpse into the rich phenomena and science to be found
in soft condensed matter. Notably, I omit an important aspect of soft matter: the
broken symmetry of the liquid crystalline state and the distinctive responses
and structures that arise from it. This realm of soft matter has benefited greatly
from the modern tools of differential geometry and field theory, as discussed
elegantly in the advanced text Principles of Condensed Matter Physics by Paul
Chaikin and Tom Lubensky (Cambridge, 1995).
The book began life in the late eighties as a joint venture with Phil Pincus
of the University of California, Santa Barbara. The chapters developed over
several cycles of teaching a course on structured fluids at Chicago. The intended
audience is the advanced undergraduate in physical science or engineering
who is comfortable with elementary physics and has seen elementary statistical
mechanics. The needed knowledge of physics is at the level of e.g., D. Halliday
and R. Resnick’s Fundamentals of Physics (New York, Wiley 2000). I assume
the student has seen statistical physics at the level of Reif’s book, the second
item of the list below. A working knowledge of Fourier transformation will also
be helpful.
Important conceptual points are illustrated by problems interspersed through
the text. The serious student should work these problems since the logical
development remains incomplete without them. Almost all of the problems
have been assigned, graded, and refined at least once. The book also includes
a number of suggested experimental projects using household materials. These
are intended to show concretely the principles discussed. They also give the
student a sense of how well the idealizations treated in the book apply to real
liquids. Students in the structured fluids course were required to do one of these
projects or some other project of their own devising. The projects are meant to
open a line of inquiry; only some have been tested, and the student may well find
that great modifications are needed to make them work. To guide the student to
further information, the book includes references to more detailed work. I have
tried to suggest sources in each major subject area. In some cases I have cited
primary journal articles, but these references are far from being thorough or
balanced. Corrections and additions to the text will be available via the Oxford
University Press web site, http://www.oup.com/.
viii Preface
In teaching the course, I found that the text contains too much material to be
covered in a one-quarter term with 27 class hours. It could probably be covered
rapidly in a one-semester course.
In learning from this book, there are several useful reference works to keep
in mind. The first and most important is de Gennes’s monograph on polymer
physics, Item 1 on the list below. Polymers serve as the paradigmatic system in
the book. Most of the ideas have been developed using polymers as an example,
then applied to other structured fluids. For further depth or clarification of these
ideas, de Gennes’s book is the best single source of which we are aware. The
subject depends heavily on ideas of statistical physics. The book develops all the
needed ideas within the text and aims to keep their number to a minimum. Still,
more explanation or depth on statistical physics may be helpful; a good source
is Reif’s book cited in Item 2. For the discussion of colloids and surfactants,
Item 3 is useful. This practical book by Jacob Israelachvili is about liquids
generally but contains much information about colloids and surfactants. For
collective properties of surfactant assemblies, Sam Safran’s book (Item 4) is
useful.
1. P.-G. de Gennes, Scaling Concepts in Polymer Physics (Ithaca NY: Cornell,
1979).
2. F. Reif, Fundamentals of Statistical and Thermal Physics (New York:
McGraw-Hill, 1965).
3. Jacob N. Israelachvili, Intermolecular and Surface Forces, 2nd ed. (London:
San Diego, CA: Academic Press, 1991)
4. Samuel A. Safran, Statistical Thermodynamics of Surfaces, Interfaces, and
Membranes (Reading, MA: Addison-Wesley Publication, 1994).
The reader should know about other classic and recent works that cover
the same material. The polymer physics material in Chapters 3 and 4 is
well covered by other textbooks and monographs. Two of the best known
are Paul J. Flory’s Statistical Mechanics of Chain Molecules (New York,
Interscience Publishers, 1969) and Charles Tanford’s Physical Chemistry of
Macromolecules (New York, Wiley 1961). There is also an important property
handbook: Polymer Handbook, 4th edition, by J. Brandrup, E. H. Immer-
gut and E. A. Grulke (New York, Wiley, 1999). An important monograph
written from a statistical physicist’s point of view is Polymers in Solution by
G. Jannink and J. Des Cloizeaux (Oxford, Clarendon Press, 1992). Another
monograph gives an authoritative treatment of polymer motions. It is The
Theory of Polymer Dynamics by Masao Doi and S. F. Edwards (New York,
Oxford University Press, 1986). Two recent textbooks offer a pedagogical
treatment of polymers: Statistical physics of Macromolecules by A. Yu. Gross-
berg and A. R. Khokhlov (New York, AIP Press, 1994), and Polymer Physics
by M. Rubinstein, R. H. Colby (New York, Oxford University Press, 2003).
A recent text covering a variety of structured fluids is Ronald G. Larson’s The
Structure and Rheology of Complex Fluids (New York, Oxford University Press,
1999). It is more advanced and more comprehensive than this book. Richard
A. L. Jones’s new text Soft Condensed Matter (Oxford, 2002), is similar to the
present book in scope and level.
Preface ix
By now the efforts of several people have improved the text greatly. The
book owes much to the generations of students who took the course, read early
drafts, and forced me to be more clear. Arlette Baljon, Alfred Liu, Keith Brad-
ley, Joe Plewa, and other students weeded out errors and confusing statements
in the text. Jung-ren Huang improved every chapter by his close reading and
thoughtful suggestions. He also compiled the data for the big table of semi-
dilute universal ratios in that chapter. Olgica Bakajin did the simulation to
make the dilation symmetry figure in the appendix of Chapter 3. Sidney Nagel
and other colleagues at the James Franck Institute took a lively interest in the
course, and the book benefited from their discussions and support. Oxford
editor Bob Rodgers provided early encouragement and obtained several very
useful reviews of an early draft of the book. Seven other reviewers looked at
the completed book last year and many improvements were made in response
to their wise suggestions. Phil Pincus inspired this entire project and wrote the
first draft of Chapter 5. He also spent many hours over several collaborative
visits to bring the book along. The book would not exist without the support
and encouragement of my wife Molly.
Chicago T. A. Witten
May 2003
This page intentionally left blank
Contents
1 Overview 1
1.1 Introduction 1
1.2 A gallery of structured fluids 2
1.2.1 Self-organization 2
1.2.2 Rheology 4
1.2.3 Scaling 5
1.3 Types of structured fluids 6
1.3.1 Colloids 6
1.3.2 Aggregates 7
1.3.3 Polymers 7
1.3.4 Surfactant assemblies 9
1.3.5 Association 10
1.4 The chapters to follow 11
References 11
2 Fundamentals 13
2.1 Statistical physics 13
2.1.1 Thermal equilibrium 13
2.1.2 Probability and work 17
2.1.3 Lattice gas 22
2.1.4 Approach to equilibrium 24
2.2 Magnitude of a liquid’s response 27
2.3 Experimental probes of structured fluids 31
2.3.1 Macroscopic responses 31
2.3.2 Probes of spatial structure 34
2.3.3 Probes of atomic environment 38
Solution to Problem 2.1 39
References 40
3 Polymer molecules 41
3.1 Types of polymers 41
3.1.1 Monomers 42
3.1.2 Architecture 43
3.1.3 Polymerization 45
3.2 Random-walk polymer 47
3.2.1 End-to-end probability 48
3.3 Interior structure 54
3.3.1 Scattering 56
3.4 Self-avoidance and self-interaction 61
3.4.1 Local and global avoidance 62
xii Contents
3.4.2 Estimating D 64
3.4.3 Self-interaction and solvent quality 67
3.4.4 Universal ratios 72
3.4.5 Polyelectrolytes 73
Appendix A: Dilation symmetry 75
Appendix B: Polymeric solvents and screening 79
References 82
4 Polymer solutions 83
4.1 Dilute solutions 83
4.2 Semidilute solutions 85
4.2.1 Structure 86
4.2.2 Energy 87
4.2.3 Concentrated solutions and melts 90
4.3 Motion in a polymer solution 90
4.3.1 Brownian motion of a sphere 90
4.3.2 Intrinsic viscosity 94
4.3.3 Polymer in dilute solution: hydrodynamic
opacity 96
4.3.4 Internal f luctuations 98
4.3.5 Hydrodynamic screening 98
4.3.6 Semidilute diffusion 99
4.3.7 Semidilute self-diffusion without
entanglement 102
4.3.8 Motion with entanglements 103
4.3.9 Stress relaxation and viscosity 105
4.4 Conclusion 109
Appendix A: Origin of the Oseen tensor 109
Solution to Problem 4.5 (Deriving permeability) 110
References 111
5 Colloids 113
5.1 Attractive forces: why colloids are sticky 114
5.1.1 Induced-dipole interactions 114
5.1.2 Solid bodies 116
5.1.3 Perturbation-Attraction Theorem 117
5.1.4 Depletion forces 120
5.2 Repulsive forces 122
5.2.1 Steric stabilization 123
5.2.2 Electrostatic stabilization 127
5.3 Organized states 132
5.3.1 Colloidal crystals 132
5.3.2 Lyotropic liquid crystals 133
5.3.3 Fractal aggregates 134
5.3.4 Anisotropic interactions 134
Contents xiii
5.4 Colloidal motion 135
5.4.1 Electrophoresis 136
5.4.2 Soret effect 137
Appendix A: Perturbation attraction in a square-gradient
medium 137
Appendix B: Colloidal aggregates 139
References 149
6 Interfaces 151
6.1 Probes of an interface 151
6.2 Simple fluids 153
6.2.1 Interfacial energy 153
6.2.2 Contact angle 156
6.2.3 Wetting dynamics 158
6.2.4 Surface heterogeneity 159
6.2.5 Other interfacial flows 160
6.3 Solutes and interfacial tension 161
6.3.1 Fluid mixtures 162
6.4 Polyatomic solutes 163
6.4.1 Polymer adsorption 163
6.4.2 Concentration profile 165
6.4.3 Hard wall 167
6.4.4 Kinetics of adsorption 168
6.4.5 Surface interaction 168
6.4.6 Flow 168
6.5 Conclusion 171
References 171
7 Surfactants 173
7.1 Introduction 173
7.2 Mixing principles 174
7.2.1 Positivity 175
7.2.2 Additivity 175
7.2.3 Ordering: like dissolves like 176
7.2.4 Reciprocity 176
7.2.5 Transitivity 177
7.2.6 Effect of permanent dipoles: water 177
7.2.7 Effect of charges: ionic separation 177
7.3 Surfactant molecules 178
7.4 Surfactants in solution: micelles 180
7.4.1 Open aggregation: wormlike micelles 182
7.4.2 Open aggregation: two-dimensional micelles 185
7.4.3 Aggregation kinetics 185
7.5 Micelle interaction 186
7.5.1 Energy of two-dimensional micelles 188
7.5.2 Energy to confine a fluctuating membrane 191
xiv Contents
7.6 Mixing immiscible liquids: microemulsions 192
7.6.1 Interfacial tension 195
7.6.2 Emulsions and foams 197
7.7 Amphiphilic polymers 198
7.7.1 Micelle size 199
7.7.2 Other copolymers 201
7.7.3 Polymeric amphiphiles in solution 202
7.8 Dynamics and rheology 203
7.8.1 Wormlike entanglement and relaxation 204
7.8.2 Rheology of lamellar solution 206
7.8.3 Shear-induced restructuring 207
Appendix: Gauss–Bonnet Theorem 210
References 211
Index 213
1
Overview
1.1 Introduction
As physical scientists, we are concerned with the behavior of matter in all
its forms. We want to know what matter does and why. This is our goal
in studying the primordial universe, the tenuous interstellar medium, the
gaseous atmosphere of the Earth, the ionized plasma of the Sun’s corona,
the mundane liquids and solids of our human surroundings, and the exotic
dense matter within a molecule, an atomic nucleus, or a proton. This book
is about a tiny subset of this vast range of forms of matter: structured fluids.
Structured fluids are liquids, i.e., condensed matter in which the atoms are
adjacent but freely mobile on a local scale. These liquids contain connected
polyatomic structures, such as solid grains or large molecules† . Though † By “polyatomic structure” we mean any
this is a small sampling of all forms of matter, it is a significant one. These assembly of many atoms that retains its
integrity over experimental time scales.
materials are important to our basic goal in studying matter: to explore the
This includes macromolecules and also
limits of what matter is capable of doing. Though they form a narrow range more weakly-bound structures.
of materials, structured fluids show a wide range of behavior. Indeed, they
show behavior not seen in other forms of matter, however exotic.
We may glimpse this distinctive behavior by recalling a few structured
fluids from everyday experience: pancake syrup, egg white, Silly Putty,
corn-starch-and-water paste, or toothpaste. Though each of these gooey
materials resembles a liquid, all respond to forces in a distinctive way that
distinguishes them qualitatively from simple fluids and from one another.
They owe their distinctive properties to connected structures much larger
than an atom, but much smaller than a macroscopic, classical body. They
belong to an intermediate, mesoscopic size regime. These structures give
raw egg white its springy consistency. They cause flowing corn-starch-
and-water paste to shatter like a brittle solid when struck. And they cause
toothpaste to flow out of its tube as a plug. These structures individually
confer distinctive properties to a liquid. More importantly, they interact
to create new forms of cooperative behavior and self-organization, as we
shall illustrate below. Some of these structured fluids have an additional
significance: they are important constituents of living cells or are important
technologically.
It was not possible until recently to study the molecular basis of
structured-fluid behavior. Now however experimental probes such as elec-
tron microscopy and dynamic scattering allow us to study structured fluids
with a resolution suited to the distinctive mesoscopic structures in them.
In addition present-day sophistication in chemical synthesis has made it
2 Overview
Fig. 1.1
Transmission electron micrograph of the
bi-continuous copolymer domain structure,
reproduced from [1]. © 1987 American
Chemical Society. Repeat distance is about
100 nm. The minority species in black
forms two disjoint domains, labeled 1 and
2. The sketch at right, after [2], shows the
three-dimensional structure inferred from
such micrographs and from x-ray
diffraction studies. © American Physical
Society. The vertical dashed line gives the
line of sight for the micrograph. The
thickness of the two minority domains is
reduced to make the illustration clearer.
Domains 1 and 2 are shown with different
shadings to distinguish them. The
magnified sketch at bottom shows the
orientation of individual polymers in the
domains.
A gallery of structured fluids 3
Water
Octane
Fig. 1.2
Freeze-fracture micrograph of a
water–oil–surfactant microemulsion
reprinted with permission from [4].
50
© 1988 American Chemical Society. The
smooth hills and valleys are thought to
represent the boundaries of oil- and
water-containing regions.
4 Overview
oil and water into these mesoscopic domains. In later chapters we shall be
interested in how small surfactant molecules can induce oil and water to
mix at this coarse scale. We shall consider what features of the molecules
govern the size of the domains and discuss how the evident randomness of
the interfaces may be characterized quantitatively. Such questions are of
broad significance, because surfactant-mediated interfaces between liquids
are ubiquitous. In everyday life and in technology surfactants at interfaces
are used to control the mixing of oily and watery liquids. And living cells
use these surfactant interfaces as a fundamental structural building block.
Mesoscopic structures can create macroscopic self-organization. The
spatial scale of the stripes in Fig. 1.3 is a few millimeters. This is a pho-
tograph of a black droplet of “ferrofluid” sandwiched between two glass
plates a few millimeters apart. The plates and fluid are in a magnetic field
of several hundred Gauss, pointing out of the figure. As the field increases
from zero the ferrofluid droplet goes from a circular shape to the labyrinthine
Fig. 1.3 shape seen here. Evidently this fluid interacts strongly with magnetic fields.
Labyrinth pattern of a ferrofluid drop It consists of 10-nm bits of the magnetic mineral magnetite dispersed in
(courtesy of R. E. Rosensweig) [5]. A drop
of ferrofluid was put between two
kerosene. In this fluid the polyatomic structures of interest are these bits of
horizontal glass plates about a millimeter dispersed solid. Solids dispersed in this “colloidal” form often have strong
apart. The plates were then placed between interactions with external forces such as the magnetic field used here. As we
the poles of a magnet. As the (vertical) proceed, we shall want to account for these strong interactions and explore
field increases from zero, the circular their limits.
droplet of fluid develops a fringe of spikes,
which elongate and branch to fill the area
Just as colloidal particles interact strongly with external fields, they
between the plates with the pattern shown. can interact strongly with each other, leading to further forms of
The thickness of the channels is self-organization. They can repel each other so strongly that they form
comparable to the spacing of the plates. a crystalline arrangement like that of Fig. 1.4. Figure 1.5 shows what can
happen when dispersed colloidal particles are suddenly made to attract each
other strongly. The particles assemble into a complex branched aggregate.
The figure is a transmission electron micrograph of colloidal silica particles
in water that have aggregated and then settled on a carbon film. A drop of
the water was placed on a carbon film, dried, and then placed in the micro-
scope. The aggregate looks different from a bulk precipitate: it is wispy and
tenuous. Aggregates like this are used commercially to toughen rubber in
tires and to thicken liquids like ice cream. In the chapters to follow, we shall
consider what makes aggregates take this form and we shall explore ways
to characterize its wispiness quantitatively. We shall also investigate how
such wispy objects interact with the fluid and with other aggregates near
them.
10 μm
1.2.2 Rheology
Fig. 1.4
An emulsion of water, decane oil, and The self-organizing behavior of structured fluids illustrated above is only
surfactant in which the oil droplets were part of their distinctive behavior. Further striking behavior is revealed when
purified to have uniform size and then the equilibrium of a structured fluid is perturbed. A classic example is shown
concentrated [6], © American Physical in on the right side of Fig. 1.6, a photograph of a polymer solution in a beaker.
Society, courtesy J. Bibette. The droplet
diameter is 0.93 μm. The long parallel
A rotating rod has been inserted into the center. The polymer liquid climbs
rows of particles and the crosshatched this rod like bread dough in a mixer. A simple liquid behaves oppositely;
texture indicate an incipient periodic the rotating liquid drops in the center and rises up the outer wall because
lattice. of centrifugal force. In the sequel we shall identify the forces that cause
A gallery of structured fluids 5
Fig. 1.5
Electron micrograph of aggregated
3.5-nm-radius particles of colloidal silica
from Lin et al. [7], courtesy of D. Weitz,
reprinted by permission from Nature
© 1989 Macmillan Publishers Ltd.
Aggregate was formed in water under
diffusion-limited conditions. The fraction
of the volume occupied by particles is
expected to vary as the sphere radius to
the −1.3 power whenever aggregation is
diffusion limited. This “universal” scaling
has been tested for a variety of materials.
Fig. 1.6
Rod-climbing liquid, reproduced from [8]
reproduced by permission of John Wiley &
Sons, Inc. (© 1977 John Wiley & Sons).
The control fluid in the left has the same
viscosity as the polymer solution on the
right, but is depressed in the center rather
than climbing the rod.
this rod-climbing behavior. We shall investigate what rate of stirring is (a) (b)
required to produce this effect. A less classical effect is sketched in Fig. 1.7.
This sketch depicts another polymer liquid. But the polymers in this liquid
contain surfactants on some of their monomers. When the liquid is properly
formulated, the act of quickly inverting its container produces a qualitative
change. The quiescent liquid has the consistency of motor oil; the disturbed
liquid is like half-cooled Jello. Phenomena like this show the potentialities
of these liquids for surprising behavior. We can only begin to sketch its
origins in what follows.
1.3.1 Colloids
Colloids are fluids containing compact, polyatomic particles suspended
in a liquid solvent. A familiar example is black ink, which is made from
colloidal carbon. Colloidal particles give distinctive physical properties of
fluids. Thus the colloidal form of the carbon in black ink is what makes it
absorb light efficiently. These distinctive physical properties are incidental
to many applications. Instead the colloidal particles are often present to
give the fluid some chemical property of interest. Examples are the cells
in blood, manufactured particles which remove molecules from a liquid
by adsorbing them, the light-sensitive grains in photographic film, or the
pH-buffer particles in detergent motor oil. In other cases, such as paint and
rubber cement, the colloidal particles produce the desired property not in the
fluid but in the solid that forms when the solvent dries. Powder processing of
ceramics utilizes the fluid properties of suspensions of ceramic powders to
optimize packing in order to achieve low defect bulk ceramics for structural
materials.
There are colloids that are important for the physical properties they
impart to fluids. One such colloid is the dispersion of 10-nm particles in
the ferrofluid of Fig. 1.3. In moderate magnetic fields the magnetic energy
of the particles is strong enough to alter the fluid’s energy significantly. In
a similar way, colloidal particles, because of their bulk, respond strongly
to electric fields or to flow. In emulsions and foams the dispersed particles
are liquid or gaseous droplets. These can be concentrated so much that
the dispersed particles push against each other and give the dispersion the
solid-like consistency of mayonnaise or whipped cream. Figure 1.4 shows
an example of a concentrated emulsion.
The interaction energy of two colloidal particles in a given solvent is
also magnified because of their bulk. Consequently, small changes in the
solvent can have a large effect on the interaction energy. This makes it
possible to change the interaction between two colloidal particles abruptly
from an effective hard-core repulsion to an attraction whose strength is
Types of structured fluids 7
many times the thermal energy kB T † . With such an attraction the particles † Here T is the absolute temperature and kB
must stick together when they encounter each other. The particles flocculate is the Boltzmann constant. The next chap-
ter reviews the significance of the thermal
or precipitate. This effect is exploited to sense small changes in a solution or
energy kB T .
to determine the presence or absence of a biological antigen—as in certain
early pregnancy tests.
1.3.2 Aggregates
The enhancement of particle–particle interaction between colloids makes
possible a form of self-organization not seen in simple molecular liquids—
namely colloidal aggregates like that of Fig. 1.5. Such aggregates form
when the attraction between two particles in contact is so strong that they
must stick together permanently. They cannot slide or roll around in order
to maximize the amount of contact. The result is the tenuous structure
shown in the figure. It is qualitatively unlike the dense-particle phase that
forms when small molecules precipitate from solution. The average particle
density within a radius r of a given particle decreases as r increases. Thus the
average density in an arbitrarily large aggregate becomes arbitrarily small:
the aggregates are “fractal” structures [9]. The origins and the consequences
of fractal structure will be a large concern in the chapters to follow.
Aggregated colloids show the enhancement effects discussed above for
dispersed colloids. In addition they have properties arising from their frac-
tal structure. Even though a tenuous aggregate occupies an arbitrarily
small fraction of the volume it pervades†† , it transmits forces efficiently †† The pervaded volume of an object is the
throughout the pervaded volume. An aggregate in a shear flow screens the spherical region of space that contains the
object.
surrounding solvent: the fluid is obliged to flow around rather than through
the aggregate. Since this screening inhibits flow, it enhances viscosity. Thus
the tenuous property of these aggregates makes them particularly effective
in increasing viscosity. Their fractal structure confers many other distinctive
properties to be explored in the course.
1.3.3 Polymers
Another way of producing a tenuous, polyatomic structure is to link small
molecules together into a flexible chain to form a polymer, as shown in
Fig. 1.8. The successive bonds between the monomers making up a flexible
chain have some randomness in their relative directions. Thus the directional
correlation between bonds more than a few bond lengths apart becomes
negligible; accordingly, a long polymer has the statistical properties of a
random walk. In some situations the statistical properties of a long polymer
are instead those of a self-repelling random walk. Like colloidal aggregates,
such polymers have the spatial scaling properties of a fractal: the average
monomer density of a self-repelling polymer of size R scales as R −4/3 .
This power for a simple random walk is −1. Chains whose pervaded-volume
fraction is a tenth of a percent or less can be readily produced. In a solution of
such polymers, there is room for a thousand chains in the volume pervaded
by one chain.
The flexibility of a tenuous polymer gives it properties that a tenuous
aggregate does not have. Unlike an aggregate, a polymer is not quenched
8 Overview
Fig. 1.8
Detail of a large polystyrene molecule as it
might appear in a good solvent such as
toluene. Each sphere represents a carbon
atom and one or two small hydrogen atoms
attached to it. The chemical bonds are
superimposed on one repeating unit or
monomer, and on a section of the chain
backbone. The backbone bonds may rotate
freely. A few successive backbone repeat
units are labeled 0, 1, 2, . . . . The vector a 1
for a four-monomer segment is shown.
This structure was generated by a Monte 2
Carlo computer simulation [Simulation by
M. Mondello, H.-J. Yang and R. J. Roe at 3
a
the University of Cincinnati using the Cray 4
Y-MP at the Ohio Supercomputer Center,
circa 1990, private communication], which
simulates the random rotations of the
bonds as they might occur in solution.
The simulated molecule has about 1/20 the
mass of those in a typical polystyrene
foam cup.
into some particular configuration but is free to explore the ensemble of
random bond directions. The randomness in each chain’s configuration thus
amounts to thermodynamic entropy, which may serve as a reservoir for heat
and work. The chain may be dramatically deformed by mild perturbations
without permanent effects. Thus the spontaneous thermal fluctuations in the
end-to-end distance of a flexible polymer are about as large as the average
end-to-end distance. Consequently, externally imposed energies as small as
the thermal energy kB T are sufficient to distort the shape of a polymer by
factors of order unity.
Unlike rigid colloidal aggregates, polymers may be concentrated to vol-
ume fractions up to unity. In this solvent-free limit called the melt state,
chains interpenetrate and entangle strongly. Each chain interacts directly
with hundreds of others and the forces thus communicated can produce
large, reversible deformation in each chain. The deformation of these ran-
dom chains is what produces the restoring force in a rubber band or a
plate of Jello. When a polymer liquid is abruptly deformed, it too responds
elastically over short times, like a rubber band. But over longer times the
chains disentangle and forget their initial distortions. The time scale for
disentanglement can easily be as long as seconds. By mixing polymers of
different lengths and architectures, one may produce liquids that behave
like a tough rubber on short time scales, like a weak rubber on longer
time scales and like a flowing latex at yet longer times [10]. This power
to control the storage of energy over time allows one to adapt polymer
solutions to the needs of a many-step manufacturing process, such as the
assembly of a car tire. Similarly, many everyday structured fluids such as
pancake syrup, shampoo, and paint are deliberately thickened to keep them
from flowing too much during application. The energy stored in a polymer
when it is deformed has striking effects on the flow properties of the liquid.
Elongating a chain produces spring tension along its length. This “normal
Types of structured fluids 9
stress” combines with other applied forces to accelerate each small volume
of the fluid. The elastic energy stored in the polymers can easily exceed
the kinetic energy in the flowing fluid. This produces nonintuitive flow
properties like the rod-climbing behavior of Fig. 1.6 [8]. One such effect
of commercial importance is turbulent drag reduction. Trace amounts of
a polymer can substantially reduce the power required to push a turbulent
fluid through a pipe, in spite of the small increase in viscosity due to the
polymer [11].
The deformability of a polymer has dramatic consequences when elec-
trically charged species are attached along the polymer chain at a given
density. The electrostatic repulsion in a long enough chain is sufficient to
stretch the chain from a random-walk configuration to that of a rigid rod,
for which the end-to-end distance is proportional to the molecular weight.
Polyelectrolytes, as such polymers are called, can be controlled in a fashion
not possible for neutral polymers. When the interaction along the chain is
screened, by the addition of some salt or the presence of other chains, the
electrostatic repulsion is reduced and the polyelectrolyte chain shrinks in
size. This alters fluid properties such as the viscosity.
Another important consequence of the deformability of polymers is seen
in their behavior near an interface. Even when the binding energy of a
monomer to a surface or interface is much smaller than the thermal energy
kB T , the total binding energy of a polymer chain made up of such monomers
may be several times the thermal energy. Then the polymer chain can
increase its binding by flattening itself closer to the surface, with little
cost in deformation energy. The distinctive features of the adsorbed state
will be considered later.
1.3.5 Association
The broad range of behavior discussed above is augmented even further
when we consider “association” of mesoscopic structures. By association
we mean a temporary joining together of the structures. The structures thus
joined can transmit forces and thus alter mechanical properties strongly. But
still these junctions are weak enough to break and reform over the time of an
experiment. Thus the associations alter themselves in response to the local
stress or flow in the liquid. We have already met the most prevalent form of
association: the joining together of surfactant molecules to make a micelle
fits the definition of association. Because these micelles are temporary,
the number and the average size of the micelles relaxes over time when
the micellar solution is subjected to e.g., a temperature jump. Association
behavior is common in colloidal aggregates. For example, silica aggregates
in water may be made to associate with one another via hydrogen bonds.
The result is a type of gel network whose links may be broken and readily
re-formed. Such a “network fluid” holds its shape in quiescent conditions
but flows under sufficiently strong shear. Upon removal of this shear, the
network is reestablished and the flow stops. These associating aggregates
are used to keep paint from running before it dries.
One can make polymers associate by attaching immiscible chemical
groups sparsely along the chains. Like surfactants, these groups assem-
ble themselves into micelles and form temporary crosslinks between the
polymers. Such polymer solutions can show a new form of response to
shear, namely reversible, shear-induced gelation. This is the phenomenon
sketched in Fig. 1.7. Another form of polymer association occurs in the
block copolymers of Fig. 1.1. Block copolymers, being amphiphilic, have
properties similar to those of surfactants. But the polymers’ entanglement
and deformability cause their resulting micellar microdomains to have some
distinctive properties not seen in surfactants. We discuss these in Chapter 7.
An example of commercial importance is Kraton, a rubbery polymer tipped
at each end with a small section of an immiscible, glass-forming polymer.
The ends of the Kraton polymer congregate into spherical micelles, each
containing many chain ends. The midsections of such a copolymer, being
attached at both ends, cannot disentangle themselves, and a strong, rubbery
material results. But when the material is heated above about 100◦ C the
spherical micelles at the ends melt and become more miscible, so that the
material flows. It can then be molded and processed. Such polymers, called
References 11
thermoplastic elastomers, are used e.g., in adhesive coatings and for the
elastic stripes painted on disposable diapers.
References
1. H. Hasegawa, H. Tanaka, K. Yamazaki, and T. Hashimoto, Macromolecules 20
1651 (1987).
2. M. W. Matsen, Phys. Rev. Lett. 80 4470 (1998).
3. D. A. Hajduk, P. E. Harper, S. M. Gruner, and C. C. Honeker, Macromolecules
27 4063 (1994); M. W. Matsen and F. S. Bates, Macromolecules 29 1091 (1996).
4. W. Jahn and R. Strey, J. Phys. Chem. 92 2294 (1988).
5. See e.g. R. E. Rosensweig, in Physics of Complex and Supermolecular Fluids,
ed. S. A. Safran and N. A. Clark (New York: Wiley-Interscience, 1987), p. 699.
6. J. Bibette, D. Roux, and F. Nallet, Phys. Rev. Lett. 65 2470 (1990).
7. See e.g. M. Y. Lin, H. N Lindsay, D. A. Weitz, R. C. Ball, and R. Klein, Nature
339 360 (1989); R. Jullien and R. Botet, Aggregation and Fractal Aggregates
World Scientific (1987).
12 Overview
This chapter assembles some underlying concepts that we will need through-
out our study of structured fluids. The first section reviews the main needed
ideas from statistical mechanics. For systems in equilibrium, we recall how
probabilities of states are determined, and the connection between probabil-
ity and the work to alter a system’s state. We also describe ways to estimate
the time for a system to reach equilibrium. The second section is about the
orders of magnitude common to the liquid environment of structured flu-
ids. In this section we justify the magnitude of a simple liquid’s viscosity.
The final section surveys the experimental probes that determine the useful
questions we can ask about structured fluids.
with e0 = e0 . Now we calculate the rel-
probabilities depend only on the respective energies: The energy E(c1 ) for
the first subsystem, E(c2 ) for the second subsystem and the total energy of ative probability for some configuration c̃1
of subsystem 1 using the f1+2 (c1 , c2 ) of
the two for the combined system: E(c1 , c2 ) = E(c1 ) + E(c2 ). The relative the combined system. Evidently we must
probabilities f must satisfy f1+2 (E1 + E2 ) = f1 (E1 ) × f2 (E2 ). sum the probabilities of all combined con-
One way to satisfy this requirement is f (E) = (constant)e−E/e0 , where figurations for which the simple subsystem
e0 is some constant. This f (E) in fact gives the equilibrium probability is in configuration c̃1 . Thus our desired
of a configuration c of any small subsystem whose energy is E(c). This probability f (c̃1 ) = c2 f (c̃1 , c2 ). Using
our Boltzmann formula for f1+2 (c1 , c2 ) this
probability distribution is called the Boltzmann distribution; deriving it is means
a major topic in statistical mechanics courses. By assuming this form one
readily verifies that f1+2 (E1 + E2 ) = f1 (E1 ) × f2 (E2 ); moreover, the f1 (c̃1 ) = exp(−E(c̃1 )/e0 )
constant e0 must be the same for subsystems 1, 2, and (1 + 2)† . Thus any
× exp(−E(c2 )/e0 ).
two (small) subsystems of the same overall isolated system have the same c2
value for the coefficient e0 . Problem 2.1 gives another way of understanding
how the Boltzmann distribution arises. The sum factor is a mere constant independ-
ent of c̃1 and has no significance in our
2.1. Justifying the Boltzmann distribution It was stated above that the probability relative probability. Comparing this f1 (c̃1 )
that a system in thermal equilibrium is in a configuration c of energy E(c) varies with our original formula, exp(−E(c̃1 )/e0 ),
as exp(−E(c)/e0 ). This is true for any conserved quantity in a small subsys- we see that the two expressions will be the
tem of a closed random system. To illustrate this point, consider the following same only if e0 = e0 .
16 Fundamentals
Express your answer in atomic mass units (6.023×(1023 )−1 g). (b) What would
the particle radius be in this case?
An important part of statistical mechanics is the relationship between
the temperature of a subsystem and its ability to do work. To recall this
relationship, we consider a specific example: an atom in a container in a
uniform gravitational field. This field exerts a constant force mg on the atom
and gives it a potential energy mgz where z is the vertical distance above
(say) the bottom of the container. According to Boltzmann’s principle, we
may take this atom as our subsystem. Its relative probability to be at height
z is then given by f (z) = exp(−mgz/T )† . From this one may immediately † Here our subsystem consists of only the
find the mean height height coordinate z of the atom. Since
we know the energy associated with a
∞ ∞
dzf (z)z given height without knowing the other
z = dzz p(z) = ∞ . configuration variables (such as the atom’s
0 0 0 dzf (z) momentum), we are allowed to use the
height coordinate as a subsystem. Of course
Simple integration yields z = T /(mg). it is also correct to find the probability f (z)
The mean height is evidently influenced by the temperature and the grav- starting from the full configuration prob-
itational field g. By increasing g slightly, we would do work −mg δz on ability f (x, y, z, px , py , pz ). But as in the
the system. This fluctuating work is difficult to give experimental meaning, (c1 , c2 ) note above, the effect of the other
variables is simply to multiply f (z) by an
but the average work δW ≡ −mgδz is well defined. unimportant constant.
gThus the work W done by an external agent to change g from g1 to g2 is
g1
2
δW . Expressing δz as (dz /dg)δg, or −(T /mg)δg/g, we conclude
g2
T dg g2 z1
W = mg = T log = T log . (2.1)
g1 mg g g1 z2
The point of this example is to remind you that (a) a system in thermal equi-
librium exerts forces on its surroundings due to its thermal fluctuations,
(b) it can be made to do work, and (c) the amount of work can be calcu-
lated using the Boltzmann probabilities f (c). Since the work is evidently
recoverable (by returning g to its initial value), the work done on the system
represents a form of stored energy.
liquid containing two colloidal particles, we often wish to know the prob-
ability that these are separated by a displacement r. In general the liquid has
many configurations in which the two particles have this separation. And
the number of configurations typically depends on the distance r. Solvent
molecules are nestled and packed around each colloidal particle, altering the
molecular arrangements near each particle. When two particles approach,
the “spheres of influence” around them overlap. In the overlap region the
solvent molecules must accommodate to two colloidal particles instead
of one. Thus the number of arrangements or configurations of the liquid
particles is altered. Figure 2.1 shows a schematic example of this effect. To
find the probability that the particles have separation r, we must add the
probabilities for all those configurations in the liquid with this separation.
We denote this subset of configurations by cr . Thus except for normaliza-
tion this probability is cr f (cr ). When the two particles are very far away,
the liquid near one is unaffected by the other, so that the probability is a
constant. It is convenient to normalize our probability so that this constant
is unity, the resulting relative
probability
is called the “pair distribution
function” g(r): g(r) ≡ cr f (cr )/ c∞ f (c∞ ).
Depending on the nature of the colloidal particles and the solvent, the pair
distribution function g(r) may increase or may decrease from unity as the
particles approach each other from infinity. To understand the meaning of
Fig. 2.1
these changes, we compare our fluid with the simplest possible one-particle
Schematic picture of two colloidal spheres system: an imaginary particle alone in an empty space and experiencing a
surrounded by small solvent molecules. potential U (r). The particle is weakly coupled to some thermal reservoir at
Only the solvent molecules closest to the temperature T . We shall take a potential that vanishes at infinity. To find
spheres are shown. In the upper picture, the g(r) for this simple system is an easy task. This is because there is only
the colloidal spheres are far apart, and one
sphere has little effect on the solvent
one configuration at position r. (We take the subsystem to be the coordinate
molecules surrounding the other. In the of the particle; so its momentum need not be considered.) Thus the only
lower picture, the colloidal spheres are r-dependent quantity in f (cr ) is e−U (r)/T . And thus cr f (cr ) = e−U /T
close together and the arrangement of the
and g(r) = e−U (r)/T .
solvent molecules is disrupted.
One way to interpret the probabilities g(r) for our colloidal system is to
compare with the g(r) of the simple system. In particular, one can find a
potential U (r) such that the g(r) of the simple system matches that of the
real system. Comparing the g(r)’s of the two systems, it is clear that they
are equal if U (r) is given by
U (r) ≡ −T log(g(r)) = −T log f (cr ) + T log f (c∞ )
cr c∞
of the real system. With this choice of U the simple system is equivalent to
the real one as far as probabilities are concerned. Now we may give meaning
to large or small values of g(r). When g(r) increases from one, the U of
the system is negative: the pair acts as though there were an attractive
interaction. Conversely, if g(r) is smaller than 1, the two particles act as
though they repel each other. In Chapter 5 to follow, we shall see how such
effective interactions behave in practice.
Remarkably, this effective attraction or repulsion describes not only the
probabilities of the two systems, but also their ability to do work. In the
Statistical physics 19
simple system the work required to bring the particle from infinity to r is
evidently U (r)† . The analogous work in the real system is not obvious to † It may not be so clear how this work is
deduce, since many parts of the system change when r changes. To evaluate measured. One way is to subject our particle
to an external potential Vx (r) with a deep
this work, we resort to an indirect method. We apply an external force F
minimum that holds the particle at a spe-
in the x direction to pull the two particles apart. (This indirect method is cific r = x. Then we may move our particle
easiest to imagine for a case where U (r) is attractive and much larger than from infinity to r0 by manipulating V so
T , so that there is a large probability that the particles are close together.) that its minimum moves from infinity to
We may do recoverable work on the system by increasing F slightly and r0 . We shall move slowly enough that the
thus increasing x . We shall apply the same force to the real system and particle stays in equilibrium with the res-
ervoir. The work required is the work done
the simple system and show that the work done on the two systems is the by the external force F = −dV /dr. At
same. Thus we can express the work done on the real system in terms of each point x en route, the particle is sit-
the U of the simple system. ting arbitrarily near the minimum of U +
When the force F changes, the resulting work δW is simply F δx , as V , so that dU /dr = −dV r0/dr = F .
in the gravity example above. To know the work for a given shift of F it is r0 work W done by F is ∞ F (r) dr =
The
∞ dU /dr dr = U (r0 ) − U (∞).
sufficient to know the change in x . The force affects the system through
the microscopic energy E(c). This E(c) has the form E0 (c) − F x, where
E0 (c) was the energy of configuration c before the force was applied.
To compute x we must average over all c. It is convenient to perform this
average in two stages, first summing over all c’s for a given displacement
r and then summing (integrating) over r:
3
d r c f (cr )x
x = 3 r . (2.2)
d r cr f (cr )
Using these facts we can readily express the x in terms of the pair
distribution function g(r):
3
d rxg(r) exp(F x/T )
x = 3 .
d r g(r) exp(F x/T )
This U (r) is called the potential of mean force. It is a form of “free energy.”
20 Fundamentals
F ≡ −T log f (c) . (2.3)
c
In our colloidal example, the external variable was the coordinate x and the
parameter b was −F . Thus −bδB → F δx = δW , as anticipated.
The free energy F is much used to discuss systems in thermal equilib-
rium. It is worthwhile to note how our discussion is related to standard
treatments.
1. One might worry that F is ill-defined, since it is based on the relative
probabilities f (c), whose normalization is arbitrary. If this normalization
were changed by multiplying all f (c) by a constant factor ea , one readily
sees from Eq. (2.3) that F is increased by an amount a. This ambiguity
is harmless; as with any energy, only changes in F are significant. Still,
it is conventional to fix the normalization by using f (c) = e−E(c)/T
with no other prefactor. Then, for a system with a single configuration
c, F is simply equal to E(c).
2. In macroscopic systems one is often interested in that part of the work
done on a
system which does not increase its internal energy E ≡
E f
c c c / c fc . In an isolated system the work done must be equal to
the change in E , in order to conserve energy. However a subsystem
which can exchange energy with a reservoir may transfer energy to the
reservoir when work is done on it. The energy transferred to an energy
reservoir in this way is heat. The amount of heat transferred, denoted
δQ, is evidently given by δQ = δW − δE . It can be shown that this
heat, like the free energy, can be expressed in terms of the normalized
probabilities pc . Specifically, δQ = −T δS, where the entropy S has
the form
S=− pc log pc ,
c
and pc is the normalized probability fc /( c fc ). The entropy is a meas-
ure of the randomness of the system. Many of the concepts in this section
can be expressed in terms of maximizing this entropy. However, in this
book, issues of heat are not essential and thus the notion of entropy will
play a minor role.
22 Fundamentals
counting the unwanted contribution from the μφ energy) must be given
by the change of (Fμ − (μ)d Fμ /dμ) = (Fμ − μφ ).
We now calculate these quantities explicitly:
d Fμ e−μ/T 1
φ = = −μ/T
= . (2.4)
dμ 1+e 1 + eμ/T
Thus the free energy is given by
Clearly the lattice gas differs significantly from an ideal gas. In an ideal gas
each particle may occupy any position independent of the other particles.
In a lattice gas the particles are obliged to avoid sites occupied by others. In
the dilute limit of a single particle in a large volume V the volume fraction
φ becomes 1/V and the lattice gas free energy of Eq. (2.5) has the limiting
form
F = −T log(V ),
like the ideal gas treated in the preceding Problem. In this limit the second
term in Eq. (2.5) is negligible. But it becomes increasingly important as the
volume fraction increases, making the lattice gas free energy larger than
the ideal gas free energy. In the opposite limit, as the φ → 1, the lattice
gas free energy becomes small again. A nearly full lattice is equivalent to a
dilute gas of vacancies, since F (1 − φ) = F (φ).
24 Fundamentals
say, for all x < x ∗ . (Thus in the fictitious system, U (x) = U (0).) Without
a biasing potential the random moves of x can easily take it from its initial
position to x ∗ . We denote the time to reach x ∗ with no potential as τ0 .
Specifically, τ0 is that time for which x has reached x ∗ with probability 12 .
The inverse of this τ0 is called the attempt rate. The Eyring–Kramers theory
describes how the barrier increases the escape time from τ0 .
If the barrier is high, the objects must make many “attempts”—returns to
x = 0—before reaching x ∗ . As time passes, the (un-normalized) probability
f (x) that the separation is x, approaches its equilibrium value for a greater
and greater range of x around its starting point x = 0. After sufficient time,
this range reaches x ∗ . Then f (x ∗ )/f (0) approaches its equilibrium value
of exp[−(U ∗ −U0 )/T ]. Whenever x does reach x ∗ , it may move to the right
or to the left without bias. Thus with probability of roughly 12 , it crosses the
barrier.
To see the effect of U , we first imagine the motion from the initial x = 0
state for a time τ0 in the fictitious system with U nullified by the artificial
force F . In a population of such molecules one encounters many possible
motions. As noted above, this time τ0 is sufficient for typical motions
to reach x ∗ . One way to account for the effect of U is to form a new
ensemble of motions starting from the motions of the fictitious system. One
removes some of these motions from the fictitious ensemble to account for
the time they spend in regions of high energy. These removals must assure
that the overall probability for the separation to be x is the equilibrium
probability f (x), in the potential U . In the fictitious system all x’s have
the same equilibrium probability f0 . But in the real system, x occurs with
equilibrium probability f (x) ∝ exp(−U (x)/T ). Evidently f (x)/f0 =
exp[−(U (x) − U (0))/T ]. In equilibrium one may thus recover the proper
probabilities by simply throwing away all but a fraction f (x)/f0 of all
fictitious systems that end up at x. This same prescription is approximately
correct even for our ensemble of motions in a finite time t. The probability
that x has reached x ∗ during the motion is that of the fictitious system (viz. 12 )
times the probability that the motion has survived the removal process.
This probability is f (x ∗ )/f (0) = exp[−(U (x ∗ ) − U (0))/T ]. Since the
separation has reached x ∗ in time τ0 with this probability, the average rate
at which the separation crosses x ∗ is thus 12 τ0−1 exp[(−U (x ∗ ) + U (0))/T ].
The pair of objects stay together for an average time that is the inverse of
this rate. Over times much shorter than this escape time, the two objects
may be considered as permanently bound together.
Clearly the most important influence on the escape time is the activation
energy U ∗ − U0 . Escape over a barrier is exponentially slow if its height
is much more than T . Increasing this height by an amount T increases the
escape time by a factor e. The escape time also depends on the attempt
time, but much more weakly. The Eyring–Kramers formula is applicable
whenever the barrier height is much greater than T . Remarkably, the shape
† To check for activated kinetics one con-
of the U (x) function is immaterial in this limit. The formula has been
ventionally plots the log of the measured
derived under rather general conditions, as discussed in [2]. It has also been
rate of e.g., dissociation against inverse tem-
experimentally verified for a wide range of processes† . Though we have perature 1/T . Activated kinetics gives a
illustrated this behavior in the context of dissociation, it is applicable to any straight line on this “Arrhenius plot;” the
quantity x with a region of high free energy. Activated hopping is not the slope is the activation energy U ∗ − U (0).
26 Fundamentals
only way a quantity like x can pass through a high-energy region. Another
process is quantum tunneling. This requires no thermal fluctuations, but
only the positional uncertainty all particles have by virtue of the wave
nature of matter. Such tunneling phenomena are not generally important
for the liquids we will consider, and we will ignore this and other quantum
effects from here on.
Using the Eyring–Kramers theory, one may get a rough idea of what parts
of the system can or cannot relax over the timescale of a given experiment.
For example, the energy required to break a carbon–carbon bond is 2 eV, i.e.,
80 T . The attempt time τ0 is in the range of 10−12 s. or more, as discussed
in the next section. From this one obtains an estimated escape time of some
1015 years. We shall need the notion of activated lifetimes when we discuss
colloidal stability in Chapter 5.
2.4. Metropolis dynamics: a concrete kinetic model It is useful to have a specific
model to examine phenomena like dissociation and to gauge the validity of
the Eyring–Kramers kinetic picture in specific cases. The Metropolis algorithm
[3] used for Monte-Carlo simulations provides a simple example that is easy
to implement on a programmable calculator or computer. Like the motion
of real molecules, the Metropolis dynamics consists of a sequence of small
steps. And after a long time, the Metropolis algorithm leads to the proper
equilibrium Boltzmann distribution f (x). Applied to our dissociation problem,
the Metropolis algorithm works as follows. Time is divided into discrete steps
of length t. The reaction coordinate changes by discrete, random increments
of a fixed magnitude x x ∗ . In the absence of a potential U , the motion is a
sequence of random steps. At timestep t, the current x(t) is changed by ±x,
with the + or − sign being chosen at random. In the presence of a potential U (x),
a given step results in a change of potential U . The dynamics must be modified
to account for the potential if it is to lead to the proper equilibrium distribution
f (x) ∝ exp(−U (x)/T ). The modification is that after each random step the
change of energy U is computed. If the step is found to decrease U , the
stepping process proceeds: the step is “accepted.” But if the step is found to
increase U , the step is “rejected” with probability 1 − exp(−U /T ). That is,
the position x is reset to its value before the last step. To implement this on a
computer, one chooses a random number between zero and one. If this number
is larger than exp(−U /T ), one resets the x coordinate; otherwise one retains
the last increment of x and proceeds to the next random step.
If the range of x is limited, the entire range is eventually visited, and after
a long time the probability f (x) that the coordinate is at x reaches a constant
distribution, independent of time t. One can show what this f (x) must be,
using the dynamical rule above. In general, one can express f (x, t) in terms of
f ’s at the previous timestep. In the absence of the potential U , if the coordinate
is x at timestep t, it must have been x − 1 or x + 1 at the previous timestep. If
it was at x − 1, then it moves to x with probability 12 . This 12 is the conditional
probability that the coordinate goes to x given that it was at x − 1. Thus f (x, t)
is the sum of two contributions for the two distinct states from the previous
timestep that can reach x:
q− = 12 , q+ = 1
2 exp(−U /T ), q0 = 1−q− −q+ = 12 [1−exp(−U /T )].
Magnitude of a liquid’s response 27
ẇ = ηγ̇ 2 . (2.7)
Evidently, the dimensions of η are energy per unit volume times time.
Why should a fluid’s dissipation be governed by a law like Eq. (2.7)?
What determines the size of the parameter η? To discuss these questions,
it is useful to imagine applying the flow between our plates in a different
way. We imagine not a smooth motion at a fixed velocity, but a motion in
sudden, tiny steps. We apply a small step strain, wait for a moment, and
then apply another. These steps may be made so small that the dissipation
they produce is the same as in the steady flow. But this stepwise picture
gives us a way to examine the molecular basis of viscosity. For simplicity,
we shall also imagine the simplest of fluids, a monatomic liquid like liquid
Argon. (We should not need to treat too-specific features of the fluid, since
Fig. 2.3
The isotropic liquid at the left is suddenly
viscosities depend but little on the molecular nature of a simple fluid.)
subjected to a 20% shear on the right, thus Figure 2.3 shows a fluid element before and immediately after a step
pushing the molecules together from strain. The sketch attempts to show an affine deformation of the centers—
upper-left to lower-right. that which would occur if these were embedded in rubber and then distorted.
Magnitude of a liquid’s response 29
typical speed v of the atoms. We may focus on one component of the throughout this book. Our use of these sym-
velocity, say vx . The energy E(c) of our system depends on vx in a simple bols is exactly that of DeGennes’ book on
way: E(c) = E(c ) + E(vx ), where the variables c are independent of polymers, op. cit. These weak forms of
equality are very useful for describing the
vx and the E(vx ) is independent of the c . In such cases, as we have seen scaling behavior we shall explore. The state-
above, the vx variable may be considered as an independent subsystem. In ment A
B means that (in the regime or
our case the (kinetic) energy E(vx ) is quadratic in the variable. Then one limit under discussion) A/B is a numerical
can readily calculate that E(vx ) = 12 T , as shown in Problem 2.5. This constant. This constant need not be close
equipartition theorem is true for any variable whose energy is quadratic. to 1. For example, a rod of length L and
mass m has a moment of inertia I
mL2 .
Using the equipartition theorem, the typical speeds of atoms in our simple The statement A ∼ B is weaker; it means
liquid are evidently given by mv 2
T or v
(T /m)1/2 , where m is the that A/B remains finite in the regime or
atomic mass. For the light atoms found in air or common liquids, it is† limit of interest. For example, rods of differ-
a few hundred m/s. The distance between collisions is roughly an atomic ent length L and width a have moments of
diameter a—i.e., roughly 0.2 nm. We thus estimate that the relaxation time inertia I
ML2 ∼ L3 when L a. Gen-
erally we use ∼ rather than the stronger ∝;
τ is of order a/v
10−12 s †† . A ∝ B indicates strict proportionality, with
For the tiny step strains in our experiment the stress σ must be pro- no implication of an asymptotic regime or
portional to the strain, as in any weakly deformed material. Accordingly, limit.
we define the “step-strain modulus” G0 as the ratio of initial stress to the
magnitude of the step strain: σ ≡ G0 γ .
The stored energy per unit volume w is the work done per unit volume
by the stress force. We may compute it by integrating the force σ xy
onγ the top face with the displacement dx = γ z: wxyz =
[σ (γ )xy]d[γ z]. Expressing σ in terms of the modulus G we
0 0
have w = 2 G0 γ . We see that G0 has dimensions of energy per unit
1 2
volume, like the dimensions of w, so that G0 is roughly the free energy cost
extrapolated to a strain of unity. This free energy is at least partly interaction
potential energy between the atoms; the free energy per atom is roughly the
cost of removing an atom. Such energies can be found from thermodynamic
measurements on the liquid. But they can be crudely estimated on funda-
mental grounds. Our atoms condense into a liquid phase (in equilibrium
with gas) because they attract each other. To form a phase, surface atoms
moving away from the surface under thermal motion must be opposed by
this attraction. The potential energy must thus be at least of order of the
kinetic energy, which in turn is of the order of T . Were this is not true,
the atoms could freely leave the liquid and the liquid would evaporate. But
the interaction energy cannot be too much greater than T ; otherwise the
30 Fundamentals
2.6. Rubber-band modulus The postal service rubber bands you find on the side-
walk have a cross-sectional area of 5 mm2 . They consist of hydrocarbon poly-
mers with a molecular weight per atom of about 5. The density is about
1000 kg/m3 (1 g/cm3 ). Under a 0.200 kg load, such a rubber band elongates by
a factor of 2. (a) What is the modulus of this rubber band? (b) What is the modu-
lus expressed in atmospheres of pressure? (c) Under a factor-of-two elongation
(unit strain) using the density, Avogadro’s number and the average mass per
atom, how many atoms store energy T ? Order-of magnitude estimates are fine.
There is a major effect that influences the viscosity of small-molecule
liquids, which we have not discussed. This is the so-called glass transition,
and is most important in molecules of irregular shape. As these liquids are
cooled, the viscosity grows too large to be measurable as the liquid cools
towards a characteristic temperature Tg . One can certainly not account for
this divergent viscosity by the arguments discussed above. However, one
may still ask whether the large viscosity is due to a large modulus or a long
relaxation time. The answer is the latter. A liquid near its glass transition
cannot equilibrate from an imposed distortion in a few atomic collision
times. The molecules have so little room to move and are so interlocked
that stress takes exceptionally long to relax. This phenomenon is important
for us in two ways. First, many of the background solvents of our structured
fluids have experimentally important glass transitions. This often provides
a useful means of controlling their viscosity. Second, large fluid structures
often themselves lead to high viscosities. There are some parallels between
this type of high viscosity and the glass transition, as we will discuss in due
course.
In accounting for the viscosity we have defined the fundamental response
of liquids. We’ve also identified the lengths, times, and energies that char-
acterize liquid structure and dynamics on the atomic scale. These conditions
set the backdrop against which the distinctive behavior of structured fluids
occurs.
107
G⬘
shear rate of 1/s. Comment on whether this estimate accounts for the damping
of the Jello.
In solutions, the modulus generalizes in an interesting way. A solution
containing dispersed objects resists compression of these objects into a
smaller volume of solvent. This osmotic pressure is evident if one encloses
the concentrated solute in a membrane which can pass the solvent, but not
the solute: solvent flows into the enclosure increasing the pressure there.
Present day osmometers use more sensitive methods, such as e.g. the reduc-
tion in vapor pressure of the solvent or the depression of the freezing point
due to the presence of the dispersed objects [5]. This measurement provides
a way of inferring how many solute particles are present, and of verifying
that they are well dispersed. Other systems, like polymer solutions, have
distinctive predicted behavior of their osmotic pressure, which can be tested
by this means.
2.8. Osmotic sensitivity Many commercial osmometers for water solutions work
by lowering the temperature until ice crystals begin to form. The solute has
the effect of depressing the freezing point, and this depression allows one to
measure the osmotic pressure. The solute particles are virtually all excluded
from the ice crystals. Thus when ice forms, the fluid volume accessible to the
solute is reduced. Each molecule in the water, displaces a volume vw . When it
transforms into ice, it reduces the water volume by vw . The work per molecule
done against the osmotic pressure vw . In order for the molecule to be in
equilibrium with ice, its free energy in the ice must increase by this amount.
For this purpose we may estimate that the free energy is of order T , so that
the needed shift in temperature T required to provide the needed increase
in temperature is given roughly by T
vw . (a) Find vw in nm3 given its
mass density and molecular weight. (A useful benchmark is that 1 g/cm3 =
3
0.6023 amu/Å .) (b) If an instrument is capable of measuring T = 0.001◦ ,
what is the largest volume per solute particle V that it can detect?
One of the most important macroscopic properties in a structured fluid
is almost trivial to measure. One simply determines whether or not the
system is stable as a single phase. Phase stability is by no means a foregone
conclusion in a solvent with polyatomic solutes. In some microemulsion and
colloid systems the conditions for single-phase stability and the nature of
the coexisting phases when a single phase is not stable are central questions
for study. We shall see that phase stability reflects distinctive conditions in
the large-scale spatial structure.
Recently the measurement of mechanical response of a fluid has been
extended in space as well as time. Forces can be measured across just a
tenths of a nanometer of a fluid using the Tabor/Israelachvili surface forces
apparatus (Fig. 2.5). This bold and simple device relies on the flatness of
cleaved crystals of mica. Such a surface can be atomically flat over square
millimeters of surface. Thus two such surfaces may be brought to e.g. a
nanometer spacing over many square microns. The surface forces appar-
atus is simply a pair of spring-mounted surfaces immersed in a fluid with
optical means for measuring the distance between them. With this device it
has proven possible to measure normal and shear forces at distances from
1 to 100 nm separations. The technique is sensitive: it routinely measures
surface interaction energies amounting to less than T per (30 nm)2 area.
34 Fundamentals
Force/radius (mN/m)
1 Linear alkanes
Branched alkanes
–1
vdW
–2
Fig. 2.5
Data from Jacob Israelachvili et al. from –3
the surface forces apparatus [6]. The liquid
between the mica plates is tetradecane, an
–4
unbranched hydrocarbon chain 14 carbons
long. Straight, rising segments of the curve
0 1 2 3 4 5 6 7
are unstable regions of the spring balance,
and contain no measured data. Distance, D (nm)
Fig. 2.6
AFM image (20 nm by 20 nm) of a
three-monolayer-thick film of the
surfactant cadmium aracidate (CdA2 ) on
mica, courtesy of J. Zasadzinski. Lighter
colors correspond to higher areas, and the
peak-to-valley height modulations about
0.2 nm. The two-dimensional fourier
transform is inset. See [8].
36 Fundamentals
Fig. 2.7
Scattering data obtained at National
Synchrotron Light source by Didier Roux 10–1
and Cyrus Safinya [9] on a lamellar phase
of 20% surfactant, oil, and water. The
ordinate I (q) is the intensity of scattered
radiation, in arbitrary units. The quantity q
on the abcissa means 2 sin(θ/2)/λ, where
10–2
θ is the deflection angle of the x-rays.
Three sets of data are for three different
amounts of water. The solid curves show –4 –2 0 2 4
the predicted power-law dependence on q. (q – G) 1000 (Å–1)
Experimental probes of structured fluids 37
2.9. Scattering assay A HeNe laser beam (wavelength about 500 nm) passes
through a glass container of unknown solution and onto the wall. We observe
that there is a faint glow of laser light around the central spot. The angular width
of this glowing region is about 5◦ . From this we infer that there are objects in
the solution which are much larger than atoms. (a) About how large are they?
(b) A second solution, similar to the first, shows the same type of scattering,
◦
except that the angular width is now only 2 12 . By what factor are these objects
larger or smaller than the ones in (a)?
a light bulb. It is far from the angular or energetic purity of a laser beam or source.
a synchrotron source. Thus the neutron beam must be collimated and made
mono-energetic by great sacrifices in intensity. Only a few facilities in the
world have enough intensity to do adequate scattering studies. Despite these
obstacles, neutrons offer some unique advantages. The main one is that they
permit isotopic labeling. One may significantly alter the scattering power
of an atom by neutrons if one substitutes a different isotope of the atom.
Such a change has virtually no effect on its behavior in a fluid. Thus, by
judicious attachment of different isotopes to different parts of a structured
38 Fundamentals
fluid, e.g. a colloidal core, the surfactants around it, and the solvent, one
may get much information about the spatial correlations of each component
and about cross-correlations between components.
Neutrons also give dynamic information. If a neutron scatters off a mov-
ing object, its kinetic energy and speed are altered. The distribution of
speeds in the scattered beam can be elegantly and sensitively measured by
sending this beam through a magnetic field. The magnetic field causes the
magnetic moment of the neutrons to oscillate with a specific frequency.
The neutron spin echo technique uses these oscillations to infer motions of
the scatterers in this frequency range. Frequencies between 106 and 109 /s
or so can be measured this way.
Waves and particles are used in another way to probe spatial structure
near surfaces. Specifically, they can sense the profile of some atomic spe-
cies with depth. The simplest of such method is to sense the reflectivity of
a wave incident on the surface. The profile of scatterers results in a profile
of index of refraction seen by the wave. This in turn shapes the reflectivity
as a function of angle, especially near the total-internal-reflection angle.
The method senses depths as small as a tenth of a wavelength and depths
as large as 100 wavelengths. It is not trivial to infer the index-of-refraction
profile from the reflectivity spectrum. One must assume various profiles
and calculate their spectra until a spectrum matching the measured one is
obtained. Sometimes two rather different profiles can give similar spec-
tra so that interpretation of the data is ambiguous. A less ambiguous but
more brute-force probe of surfaces is called Rutherford backscattering [11].
There the probe is a beam of massive particles like alpha particles, whose
wavelength is negligible. When a mono-energetic beam is directed onto a
surface, the particles backscatter elastically from the nuclei in the substrate,
losing a specific amount of energy that depends on the mass of the nucleus
hit. One may infer the depth of a backscattering event by the further energy
loss of the emerging particle on its way in and out. This method requires a
small room-sized particle accelerator. The method is routinely used to sense
the profile of composition in solid samples to depths of the order of 10 nm.
line, one may sense the average orientation of that part of the chain in space. Free chains ( = 0.92)
As Fig. 2.8 shows, the changes in the resonance can reveal subtle changes 100 Hz
in the atomic orientations due to overall deformation of the fluid. One can
sense the effect of a 10% stretch of a rubber in this way; likewise, one = 1.0
can sense the alignment of specific carbon–deuterium bonds in a surfactant
monolayer.
Another class of local probes involves optical fluorescence of small dye
molecules [14]. The dye molecule is chemically attached to the macro-
= 1.22
molecule to be probed. For example, fluorescence may be excited with
a pulse of light and its intensity or polarization are monitored with time.
The change of polarization senses the dye molecule’s rotation from its ini-
tial orientation. The time dependence of intensity is of interest because
the fluorescent decay can be strongly affected by the presence of quench- = 1.40
ing molecules. From the intensity decay one may infer the rate at which
excited molecules encounter quenching molecules. By attaching the dye
and quencher at suitable places, one may use this as a probe of mobility, and
of where the dye molecule likes to sit. A variant of the quencher molecule = 1.55
is an excimer-former—usually a second dye molecule. Excimers are bound
states of two dye molecules. The proximity of the second molecule shifts Fig. 2.8
the fluorescence in a specific way. Fluorescence probes of this kind are Nuclear magnetic resonance spectrum of a
very sensitive and flexible in their time scale. Temporal responses have rubber after [13]. The rubber is a polymer
been measured on scales ranging from fractions of nanoseconds to many made of hydrogen, carbon, oxygen, and
silicon (poly dimethyl siloxane) in which
seconds. the hydrogen atoms have been replaced by
deuterium atoms. The plot shows
absorption versus frequency of oscillation.
A strip of the rubber was stretched in the
Solution to Problem 2.1 direction of the static field by various
The outer part of the system (the “reservoir”) has N − k ≡
m balls in U − V ≡ v elongation factors λ shown at the left. The
sites. The number of distinct configurations with m balls {dm } 1 will be called splitting of the resonant peak indicates a
resulting slight anisotropy in the
the “partition function” zm . We wish to infer zm+1 given zm . Any m + 1-ball
neighborhood of the deuterons. The
configuration may be made from some m-ball configuration by adding one ball.
splitting is proportional to the amount of
Thus to count all m + 1 configurations we may consider each m-ball configuration
stretching. The frequency of the unsplit
and ask how many m + 1-ball configurations can be made from it.
peak is 13 MHz. The bottom curve shows a
(a) Whatever m-ball configuration we start from, there are exactly v m + 1-ball splitting of about 100 Hz.
configurations that can be made from it, since every possible site to add the
new ball results in a distinct configuration.
(b) Each of these m + 1-ball configurations could have been made from m + 1
different m-ball configurations. If we remove any one of the m + 1 balls, we
get one m-ball configuration. If v m, the chance of two balls at a site is
negligible and so essentially all these m-ball configurations are distinct.
(c) Not all the configurations found in (a) are distinct: one can arrive at the same
m + 1-configuration starting from many different m configurations. In fact, if
we generate all the m + 1-configurations separately from each m configuration,
we will end up generating each m + 1 configuration m + 1 times—once for
every m configuration that could have given the m+1-configuration in question.
Thus of vzm configurations generated, a factor m + 1 are redundant: they are
duplicates of others. The number of distinct m + 1 configurations zm+1 is thus
(v/(m + 1))zm . We recall that the desired probability p(ck ) = (const)zN−k .
Thus we have found that p(ck+1 )/p(ck ) = zN −k /zN −k−1 = (N − k)/v. Since
N k, this ratio becomes N /v and is independent of k, as we wished to show.
(If the reservoir is not dilute, the ratio p(ck+1 )/p(ck ) changes. It becomes a
function of the volume fraction N /U of the reservoir, because the number of
40 Fundamentals
References
1. See e.g. N. G. van Kampen, Stochastic Processes in Physics and Chemistry
(Amsterdam, New York: North-Holland, 1992).
2. S. Glasstone, K. J. Laidler, and H. Eyring, The Theory of Rate Processes,
(New York: McGraw Hill, 1941); H. A. Kramers, Physica 7 284 (1940);
S. Chandrasakhar, Rev. Mod. Phys. 15 1 (1943).
3. See e.g. H. Gould and J. Tobochnik, An Introduction to Computer Simulation
Methods (New York: Addison Wesley, 1988).
4. D. S. Pearson, Rubber Chem. Technol. 60 439 (1987).
5. C. Morris and H. Coll, Determination of Molecular Weight (New York: Wiley
and Sons, 1989).
6. Reprinted from Jacob N. Israelachivili, Intermolecular and surface forces, 2nd
ed. (London; San Diego, CA: Academic Press, 1991). p272. © 1991 with
permission from Elsevier.
7. R. Weisendanger and H. J. Guntherodt, Theory of STM and Related Scanning
Probe Methods, Springer Series in Surface Sciences, V 29 (1993).
8. J. A. Zasadzinski, R. Viswanathan, L. Madsen, J. Garnaes, and D. K. Schwartz,
Science 263 1726 (1994).
9. B. J. Berne and R. Pecora, Dynamic Light Scattering: With Applications to
Chemistry, Biology, and Physics (New York: Wiley, 1976); Dynamic Light
Scattering: Applications of Photon Correlation Spectroscopy, ed. R. Pecora
(New York: Plenum Press, 1985).
10. A. Q. R. Baron, H. Franz, A. Meyer, R. Ruffer, A. I. Chumakov, E. Burkel, and
W. Petry, Phys. Rev. Lett. 79 2823 (1997).
11. For this and related methods, see T. P. Russell, Ann. Rev. Materials Sci. 21 249
(1991).
12. See e.g. V. D. Fedotov and H. Schneider, Structure and Dynamics of Bulk
Polymers by NMR-Methods, NMR Basic Principles and Progress series, Vol. 21
(Heidelberg: Springer-Verlag, 1989).
13. B. Deloche and A. Dubault, Europhys. Lett. 1 629 (1986).
14. See e.g. C. E. Hoyle and J. M. Torkelson, eds., Photophysics of Polymers, ACS
Symposium Series Vol. 358 (Columbus, Ohio: American Chemical Society,
1987).
3
Polymer molecules
In this chapter we have our first in-depth encounter with a specific fluid
structure: the flexible chain molecules known as polymers. The basis for
much of their distinctive behavior lies in their flexibility and randomness.
This chapter focus on the properties of individual molecules that embody
this flexibility and randomness. The next chapter focuses on the further
properties that emerge when large numbers of polymers interact in a liquid.
There we will also treat the basic motions of a polymer solution: diffusion
and flow. The Preface notes a number of other texts covering these polymer
phenomena.
We begin this chapter with an informal introduction to some common
monomers and the ways they can be joined to make a polymer chain. Next
we use the idealized example of a random-walk polymer to demonstrate
scaling and renormalization. The following section describes the interior
structure of polymers, emphasizing their fractal properties and dilation
symmetry, and showing how these features can be revealed by scattering
experiments. Next, we study the impact of self-avoidance on the structure
of a polymer. We find that self-avoiding polymers have fractal structure
but the scaling exponents are altered by the self-avoidance constraint. We
discuss another form of self-interaction that occurs when the polymer is
electrically charged; this leads to dramatic changes in its structure.
3.1.1 Monomers
†For reference, carbon in organic molecules The simplest and most common polymers are hydrocarbons, i.e., molecules
has a valence of four; i.e., it makes four consisting entirely of carbon and hydrogen. It is useful to classify the differ-
bonds to adjacent atoms. Two of these may
ent types according to their polarity, i.e., the degree of electric polarizability.
go to the same atom, forming a double bond.
Silicon also has a valence of four. Nitrogen The least polarizable hydrocarbons are the saturated hydrocarbons, where
has a valence of three. Oxygen and sulfur each carbon is bonded to four different atoms (i.e., there are no double
have a valence of two. Hydrogen has a bonds)† . The bonds in a saturated hydrocarbon are very stable, symmetrical,
valence of one. and difficult to deform. Small saturated hydrocarbons are oils and waxes.
All these species mix well with each other but less well with the more polar
†† We return in Chapter 7 to a more quant- species to be described below. They mix poorly with water†† . Polypropyl-
itative discussion of molecular mixing. ene is the main example besides polyethylene. It has a CH 3 group replacing
one of the hydrogens on every second CH 2 group of polyethylene††† . Poly-
††† A polymer is generally named after the
ethylene and polypropylene are structurally simple enough to crystallize
monomer from which it is made. For readily. This contrasts with another very common polymer, polystyrene,
example the ethylene monomer has four
the material of cheap plastic spoons and styrofoam cups. Figure 1.8 shows
hydrogens surrounding two double-bonded
carbons. In polymerization, the second a typical configuration of polystyrene. It is related to polypropylene by
carbon–carbon bond of ethylene becomes replacing each CH 3 group of polypropylene by a benzene C6 H5 ring. The
attached to the end of the polyethylene bulky rings make crystallization difficult; thus polystyrene solidifies in a
chain. Sometimes a polymer has two differ- glassy rather than a crystalline state.
ent names, because either of two monomers
Polystyrene is an unsaturated hydrocarbon. That is, it contains double
can be used to make it.
bonds (between some carbon atoms and their neighbors). Some of these
redundant bonds could be used to attach additional hydrogen atoms, thus
“saturating” the molecule with hydrogen. The redundant bonding in poly-
styrene resides in the delocalized electrons in the benzene rings. These, like
other double bonds, are more polarizable than simple bonds. This means
that polystyrene (at high molecular weight) is not miscible in saturated
polymers. Two other important unsaturated polymers are natural rubber
or polyisoprene, and the synthetic rubber polybutadiene. Several common
hydrocarbon polymers are sketched in Fig. 3.1.
Polyethylene The inclusion of elements beyond carbon and hydrogen makes for asym-
Polyisoprene
metric bonds with carbon, and often increases polarity. For example, one
may make a carbon chain with oxygen at every third position. The simplest
case is [−−CH 2 −−CH 2 − −O −−]n —polyethylene oxide§ . Polyethylene oxide
is sufficiently polar to be soluble in water. Polymer backbones can be made
Polypropylene Polystyrene with no carbon at all. The most common case is polydimethyl siloxane.
This is the material of silicone rubber, sealants, and lubricants. It has
Fig. 3.1 a backbone of alternating silicon and oxygen atoms. Each silicon has in
Four common hydrocarbon polymers. addition two methyl (CH 3 ) groups. Polyesters such as DacronTM have the
Black dots represent carbon atoms. Solid
lines represent chemical bonds. A
generic structure [− −R −−COO− −]n , where the R is an arbitrary molecular
hydrogen atom is at the end of each bond substructure called a residue and COO is a carbon with a double-bonded
where no carbon atom is shown. Dashed oxygen and a single-bonded oxygen taking three of its four available bonds.
lines represent bonds from the preceding The last bond is with the R group. The double-bonded oxygen may form
monomer in the chain. Variants of weak hydrogen bonds with any covalently-bonded hydrogen present, e.g.,
polybutadiene are shown in Fig. 3.3.
in the R group. Another example is the polypeptides; all proteins are
polypeptides. Their backbone has the form [− −N −−C −−C − −]n , with an H
§ also known as poly-oxy-ethylene, or poly- taking the third bond of the nitrogen, an H and a variable R taking the
ethylene glycol (PEG). remaining two bonds of the middle carbon, and a double-bonded O taking
the two bonds of the last carbon. In many polypeptides (with favorable
Types of polymers 43
R groups) this oxygen forms a hydrogen bond with the hydrogen sev-
eral monomers down the chain, thus forming a helical structure called the
alpha helix. The alpha helix is one example of weak, cooperative bonding
within a chain which favors a certain configuration or crosslinking topo-
logy. This type of bonding is called secondary structure (in proteins) or
association (in synthetic polymers). Secondary structure can often be dis-
rupted reversibly without disrupting the primary, covalent bonds of the
polymer.
Other types of monomers respond more strongly to electric fields than can
be described by a polarizability. Conjugated polymers like polyacetylene
have alternating single and double bonds along the backbone. Moreover,
translating all the double bonds forward by one unit along the structure res-
ults in an identical structure. The position of the double bonds represents a
broken symmetry. By doing this translation, one has transported an electron
in each double bond one step along the chain. This means that such poly-
mers may readily be made into electrical conductors. This requires doping,
the removal or addition of a small fraction of the electrons. Conjugated
polymers can have electron mobilities comparable to those found in metals.
This conductivity can be made to occur in the liquid state and in crosslinked
solutions [1].
Polymers can be synthesized with ionic side groups. A common example
is polystyrene sulfonate. This derivative of polystyrene has an SO− 4 ion
attached to each benzene ring. The neutralizing counterion is typically a
metal like Na+ . In a polar solvent like water, the electrostatic binding energy
of the counterions is smaller than the thermal energy T , and accordingly
most of the ions dissociate and move freely in the solvent. The polymer is
left with a large net charge. Such polymers are called polyelectrolytes. We
will see later in this chapter that the charge has a dramatic effect on their
shape and properties.
3.1.2 Architecture
In polymerization monomers like those catalogued above are joined end
to end. The result is generally a flexible chain molecule, in which each
single bond can rotate. That is, the bond acts like a rigid rod, with the
atoms at the two ends attached as if on freely rotating bearings. The bonds
of a carbon atom point symmetrically outward at fixed angles of about
109◦ . A double bond may be regarded as a pair of rods; it does not allow
rotation. Even in a flexible polymer chain certain aspects of the structure
are locked in and unchangeable, over and above the fixed structure of the
constituent monomers. This structure is determined at the moment each
monomer joins the chain. Two types of locked-in structure are important:
isomerism and tacticity. Isomerism concerns the orientation of the two ends
of the monomer along the chain. Each monomeric link is attached to the
chain at its two ends. In general the two are not equivalent, so that monomers
may be attached “head to tail” in the same orientation, or “head to head”
in opposite orientations. How much of this isomeric irregularity occurs
depends on the details of the polymerization reaction. Of the polymers
shown in Fig. 3.1, only polyisoprene has isomeric variation of this type.
44 Polymer molecules
Fig. 3.2 R R R
H H H
A generic polymer chain with a C2 repeat
unit along the backbone. Bonds directed
out of or into the backbone plane are shown
as triangles. The residues R are at the
isotactic positions; they are all on the same H H
side (the front) of the backbone plane. H H
3.1.3 Polymerization
To create a polymer, one must find a chemical reaction mechanism that
joins the monomers together. The organic chemistry of these reactions
[3] is similar to that of small molecules, and we shall not discuss it here.
But some important features of polymerization can be appreciated without
organic chemistry. For example, statistical fluctuations are important in
the many-body process of polymerization. These produce variations in
isomerization and tacticity along the chain. More centrally, they produce
variations of molecular weight from chain to chain. We discuss three types
of polymerization to illustrate the possibilities.
The least variability occurs in biosynthesis, the natural process by which
proteins are synthesized in living organisms. By recombinant DNA techno-
logy the natural processes can now be commandeered to produce peptide
copolymers with the precise sequence desired [4]. These are virtually
identical in structure and molecular weight from chain to chain. This kind of
synthesis is still experimental as a means of producing synthetic polymers
and is restricted to polypeptides and DNA-like polysaccharides.
Some of this molecular uniformity occurs in the important type of poly-
merization called addition polymerization. In this method a catalyst initiates
polymerization in a solution of monomers. This turns a small fraction of
the monomers into a chemically active form such as a free radical. An
example is the polymerization of styrene monomers CH2 = CH R. In
the active form the double bond is broken and made reactive. This act-
ive form is able to combine with a second monomer. The activity of the
first monomer propagates to the second. Thus it is able to react with
a third, and so forth. Ideally many monomers are joined in sequence
before the activity of the chain end is lost. Finally the reaction is termin-
ated by deliberately adding a chemical that inactivates the active ends.
Another example of addition polymerization is shown in more detail
in Fig. 3.3.
The successive addition of monomers to the chains is a statistical process;
still, it results in relatively uniform chains when these chains are long. The
distribution of lengths P (n) becomes relatively narrower as the mean chain
length n gets larger. This results from the law of averages. We may see how
this works by considering the time T required to make an n-mer. Evidently
T is the time t2 required to attach the second monomer to the first, plus the
time t3 for the third monomer, and so on up to tn . Of course each of these
ti varies widely from chain to chain, owing to the statistical fluctuations in
the system when the ith monomer was being added. These variations are
reflected in the total time T , by the familiar law of propagation of errors:
T 2 = ti2 . The variations of the ti are roughly the same for all i,
so that T 2
ntn2 . Thus T
T n−1/2 ; the relative variation in T
becomes progressively smaller as n increases.
This increasing uniformity contrasts with the results of condensation
polymerization. Here instead of each chain having a single active end that
reacts only with monomers, many chain ends may react with one another.
An example is the polyamide nylon, in which any two chain ends may join,
producing a water molecule as a by-product (Fig. 3.4).
46 Polymer molecules
Fig. 3.3
M M
Addition polymerization of butadiene. *
Proper bond angles are not shown. The
initiation of the reaction by a metal atom *
M produces an unpaired electron, denoted
by a star. Under certain reaction
conditions, this unpaired electron resides
preferentially on the second carbon atom.
This unpaired electron then initiates M M
bonding with the end carbon of a second *
butadiene. It then moves to the second
*
carbon, as before. Additional monomers
are incorporated in the same way. The
resulting “1,2 polybutadiene” has
two-carbon side groups. Each may extend
either out of or into the plane of the
backbone, so that the chain has tacticity as
illustrated in Fig. 3.2. Alternatively, under
Cis Trans
different reaction conditions the unpaired
electron moves to the carbon atom farthest
from the metal initiator. This electron
induces bonding with the end carbon of
The kinetics of this reaction involves many processes contributing to the
another monomer, leaving an unpaired
electron that moves again to the farthest change of the n-mer concentration [n]. For simplicity we assume that all
carbon atom. Successive addition of n-mers are equally reactive with one another. Each reaction of i-mer with
monomers in this way produces “1,4 j -mer to produce n-mer is a second-order reaction, whose rate is propor-
polybutadiene.” Double bonds do not tional to the amount of each reactant present, i.e., to [i][j ]. Considering all
rotate freely like single bonds. Thus, the
these reactions that affect the number of n-mers, we have:
four backbone atoms around any double
bond can be in either the planar “cis” or
1
“trans” isomeric states shown in ovals. n−1
d[n]
=K [k][n − k] − [n] [k] . (3.1)
dt 2
k=1
H H H H Fig. 3.4
C C C C O H A polycondensation reaction. An amine
H O C
NH 2 end reacts with an alcohol COH end
O H H H H by eliminating the circled atoms to form a
H H H H H H H H
water molecule, thus joining the two
N C C C C C C N
H H chains together. The polymers thus
H H H H H H produced are called nylon 6-6.
considering all the reactions which increase or decrease the number of n-mers,
the rate of increase of the concentration [n] obeys Eq. (3.1) in the text, where K
is an overall rate constant. In analyzing this equation, you may assume that n, k,
and l are everywhere large enough that they may be considered as continuous
variables running from 0 to ∞ instead of discrete variables running from 1 to
∞. (a) Show that the total mass n[n] is independent of time, as it must be. (b)
[n](t) = [0](t) exp(−nc(t)), where2 c
Show that there is a solution of the form
and [0] are some function of t. (Since n[n] = [0]/c2 , [0] must be (const)c .)
How does c vary with t? (c) How is c related to n ≡ ( n[n])/( [n])?
angles as the third one did. Here the last bond may point in any direction
relative to the first, and these directions are nearly equally likely. With
the addition of another bond or two, the directional randomness of the last
bond is nearly complete. To discuss the randomness in a long polymer, it
is convenient to imagine that its building blocks are segments which are
several atoms long, so that the directions of the beginning and the end of
† In some cases this directional randomness each segment may be assumed statistically independent† .
cannot be assumed. An important one is the Another feature appears when our hydrocarbon chain exceeds five bonds
case when the backbone is rigid rather than
in length: it may hit itself. There are some rotational configurations of the
flexible. Important examples are DNA and
RNA helices and straight-chain hydrocar- bond angles that would place one carbon atom on top of another. Clearly
bons like polyacetylene. Here each segment these configurations are absent in real polymers: the polymer must avoid
has nearly the same direction as its pre- itself. This requirement is clearly separate from the requirements on the
decessor. The direction can only become successive bond angles or local structure. Local structure can be expressed
random by the addition of many of these ran-
in terms of random variables (e.g. our bond angles) that are statistically
dom small increments over a long length of
chain. Chains shorter than this persistence independent for each bond or small segment. The self-avoidance constraint
length but much longer than a monomer cannot be expressed in this way. A major topic of this chapter is to analyze
have a universal statistical behavior more the effects of self-avoidance. For the moment, we shall ignore it. We shall
complicated than that of the random-walk depart from the real world and imagine phantom polymers which are free
polymers treated in the text. Rigid chains in
to intersect themselves. After their behavior is understood, we will be in a
this regime are called wormlike chains. [5].
position to attack the constraint of self-avoidance.
3.2. Law of Averages for r 2 A polymer in one dimension consists of n vectors
a1 , a2 , . . . , an of length 1 that point randomly left or right. Each vector is
independent of the others. (a) Find a5 , · a5 , and
a5 a4 · a5 . (b) Using
averages such as these, and noting that r = i ai , find r and (r )2 as a
function of n.
We note that the first integral is the average of r for a one-segment chain:
r 1 and the second is r 2 1 .
We have taken a case where the step distribution p0 (r) is isotropic; thus
the first integral must vanish. The same is true with all odd powers and their
associated odd-order derivatives. Thus our equation for p reduces to
Here we have indicated the form of the first nonvanishing term beyond the
∇ 2 term; the “∇ 4 ” denotes various fourth-order derivatives.
50 Polymer molecules
We can deduce much about this equation without solving it, by invest-
igating its asymptotic behavior as n → ∞. But it will not do simply to
evaluate p(∞, r). We would find that this is zero for all r. We are in the
paradoxical situation of taking n arbitrarily large, but not infinite. Since
it is arbitrarily large, the specific value of n cannot be important. In the
regime of interest the behavior of p is the same whether we consider n,
n/1000, or 1000n. That is, p should show the same behavior when n is
multiplied by an arbitrary factor λ. By assuming such a limiting behavior,
we are implicitly postulating some kind of symmetry; namely, invariance
under multiplicative scaling of n.
We anticipate that a change of n by some factor λ preserves the shape
of p(n, r), but not its specific value. That is, p(λn, r) should differ from
p(n, r) by a mere change of scale on the r axis and the p axis. We express
the r scale factor as μ(λ) and the p scale factor as η(λ). Then our anticipated
λ variation takes the form p(n, r) = η(λ)p(λn, μ(λ)r).
We may use Eq. (3.6) to decide whether p(n, r) indeed has an asymptotic
behavior of this form. Substituting the proposed form,
dηp(λn, μr) 1
= r 2 1 ∇ 2 ηp(λn, μr) + const. “∇ 4 ”ηp(λn, μr) + · · · .
dn 6
(3.7)
To see the effect of λ and μ on the equations, we define the rescaled λn as
ñ, and the rescaled μr as r̃. Similarly, ηp ≡ p̃. We may readily express
Eq. (3.7) in terms of these renormalized variables. Evidently, d/dn =
λd/d ñ and ∇r2 = μ2 ∇r̃2 . Likewise, any fourth-order derivative satisfies
∇r4 = μ4 ∇r̃4 . In these new co-ordinates,
d p̃ 1 μ2 μ4 4
= r 2 1 ∇r̃2 p̃(ñ, r̃) + const. “∇r̃ ”p̃ + · · · . (3.8)
d ñ 6 λ λ
We may now freely take n to be arbitrarily large; we simply select a
compensating scale factor λ so that ñ remains finite. Evidently, λ → 0
in this limit. In order for Eq. (3.8) to remain finite in this limit, μ must
evidently change with λ. If the ∇ 2 term is to remain finite, we must have
† Technically, the n derivative on the left
μ2 /λ = constant. With this choice, the ∇ 4 term has a factor μ4 /λ = λ.
side of our equation is also an approxima- The ∇ 4 term is arbitrarily small, so that it can be neglected. Such quantities
tion. Like the spatial derivatives on the right
that become arbitrarily small when an asymptotic limit is taken are termed
side, it should be replaced by a full Taylor
expansion. But the higher terms not shown irrelevant variables. We may clearly extend this reasoning to show that all
here may be readily shown to be irrelevant higher orders in ∇ are irrelevant in this same way† . Likewise, just as the ∇ 2
like the higher spatial derivatives. term involved a moment r 2 1 , one can readily show that the ∇ 4 terms have
factors of the form r 4 1 . By our reasoning these moments of p0 (r) are also
†† One may ask why we chose μ(λ) to make irrelevant for the asymptotic p(n, r). Indeed, the only feature of p0 (r) that is
the first term in our equation finite. By relevant is r 2 1 †† . By similar reasoning we
may determine the scale factor
choosing μ ∝ λ1/4 we can make the second η for p. A simple method is to note that d 3 rp(n, r) = 1: the probability
term finite instead. But this choice does not that the end monomer is somewhere is a certainty. Inserting the scaling form
lead to a consistent equation. Now the first
and noting that d 3 r = μ−3 d 3 r̃, we infer η = μ3 = λ3/2 .
term has a factor λ−1/2 . Instead of being
irrelevant, the first term diverges, domin- The requirement that n and r must scale together means that the asymp-
ating the one we have rendered finite. We totic p(n, r) is really only a function of one variable. For any given n
must thus reject this choice of μ. and r, we are free to set λ = 1/n. Then recalling that μ ∝ λ1/2 ,
Random-walk polymer 51
p(n, r) = n−3/2 p(ñ, r̃)|ñ=1,r̃=rn−1/2 † . By use of the scaling symmetry, † Here p(1, r̃) means the asymptotic func-
we have reduced our unknown function of two variables, p(n, r), to an tion p evaluated for ñ = 1. It does not mean
the monomeric p0 (r) from which p(n, r)
unknown function f of a single variable:
was constructed.
= d 3 (rn−1/2 )[(rn−1/2 )p np/2 ]f (rn−1/2 ),
= np/2 d 3 r̂ r̂ p f (r̂), (3.9)
where the dummy variable r̂ stands for rn−1/2 . Since the integral is inde-
pendent of n, we infer that any measure of the polymer size of the form
[r p ]1/p scales the same with n:
[r p n ]1/p = Kp n1/2 . (3.10)
The constants Kp are numerical factors related to the r̂ integral in Eq. (3.9).
This same type of scaling appears whenever we do mathematics with
dimensioned quantities like mass, length, and time. Here we know that the
number representing, say, the time in our formulas has no significance in
itself. The formulas must give the same physical result whether we use one
unit of time or some other unit λ times smaller. Thus the number repres-
enting time can be multiplied by an arbitrary factor λ without changing the
physical result of the formula. This is the familiar requirement of dimen-
sional consistency: the answer must be unchanged when all quantities are
scaled in accordance with their time dimensions by an arbitrary factor λ.
Mathematically, this is a requirement of homogeneity. If we express our
formula in the form f (t, a, b, c, . . .) = 0, where t is the time variable, then
dimensional analysis amounts to the following invariance property: for each
physical quantity a, b, c, etc. in our formula there are exponents α, β, γ ,
etc. such that
f (λt, λα a, λβ b, λγ c, . . .) = 0. (3.11)
The exponents, usually integers or simple fractions, are the time dimensions
of the quantities a, b, c, etc. Though n in our polymer has no physical dimen-
sions, it has something in common with a dimensioned variable: namely,
the invariance of the results under changing the scale factor λ. The resulting
symmetry is the same homogeneity familiar in dimensional analysis: thus
asymptotic quantities like n behave as though they had dimensions. This
will prove true in the more subtle forms of scaling to be treated later.
52 Polymer molecules
Returning to our polymer, the asymptotic equation now has a simple and
familiar form. It is the heat equation, the equation that describes how a
† The polymer scaling is familiar in the con- diffusing substance spreads out with time† .
text of diffusion. It is the counterpart of the
fact that the distance a diffusing substance dp/dn = 1/6r 2 1 ∇ 2 p(n, r), (3.12)
spreads is proportional to the square root of
the spreading time. Here we have set λ = μ = 1 now that they have served their purpose. One
may readily verify that the solution is a Gaussian:
2 −3/2 3 2 2
p(n, r) = [2π nr 1 /3] exp − r /(nr 1 ) . (3.13)
2
U 3.3. Toy renormalization In the text we deduced a way of inferring polymer beha-
vior by analyzing the scaling properties of the equation without solving it.
This problem is meant to show the same procedures in the familiar context
of classical mechanics. Here we wish to find the limiting properties of a sys-
x tem’s motion as the energy of a particle goes to zero. The potential energy
0 function U (x)
of a point particle has the form of two semicircles side by side:
U (x) = U0 1 − (1 − |x|)2 . A particle with energy E U0 is trapped in the
region near x = 0. It is obliged to oscillate with some period τ . The object
is to find how the period τ varies as the energy E → 0. (a) Write conserva-
tion of energy √for the particle in the region x 1. Notice that in this region
U (x)
U0 2|x|(1 − 14 |x|). (b) In order to find the limiting behavior as
E → 0, define Ẽ = λE so that Ẽ will stay finite when E → 0 when expressed
in terms of the ˜ variables. We wish to find a corresponding time scale t˜ = μt
and distance scale x̃ = ηx such that the equation of motion will stay finite
when λ → ∞. (c) Restate the conservation of energy from (a) in terms of the
˜ variables. (d) How must μ and η vary with λ in order for the terms in the
equation to stay finite? (e) How do the maximum displacement in x and the
temporal period τ vary with E as E → 0? ( f) If μ and η vary with λ as found
in (b), does the correction in 14 |x| become small or become large?
This function describes the end-to-end distance of long polymers.
Remarkably, the microscopic structure of the chain, given by p0 (r), has
little impact on p(n, r); the only relevant feature of this p0 (r) was its
mean-square average r 2 1 . Since virtually any p0 (r) leads to the same
asymptotic p(n, r), we say that the Gaussian p(n, r) is universal.
It is worth noting that this Gaussian is not the asymptotic p(n, r) in every
sense. To see this, we might consider our hydrocarbon polymer for very
large extensions r. The polymer has a maximal extension where the back-
bone bonds form a coplanar, zigzag line. Beyond this distance p(n, r) is
strictly zero, so that the Gaussian form is qualitatively wrong for such val-
ues of r. This rmax is evidently proportional to n, so that it scales entirely
differently from r 2 n . This illustrates the point that asymptotic scal-
ing behavior depends on the question asked. In determining the Gaussian
scaling, we asked for the behavior of p(n, r) for “typical” r values with
appreciable probability. The rmax is extremely atypical; the probability that
r is any finite fraction of rmax is exponentially small in n, as one may
verify using Eq. (3.13). For our purposes in this book, we shall be inter-
ested in typical polymer configurations, so that the Gaussian scaling is the
right regime to consider. But there are some questions of physical import-
ance where the Gaussian treatment is not sufficient. The exercise on the
coil-stretch transition gives an example.
Random-walk polymer 53
3.4. Short lattice chains A four-step polymer in one dimension consists of left-
and right- directed bonds of unit length. (a) How many configurations of the
chain have the end four units to the right of the beginning? How many have
the end three steps to the right? Two? One? Zero? You may find it easiest to
draw a picture of all 24 configurations. (b) Find r 4 /r 2 2 for this polymer
and compare with 3, the corresponding ratio for a one-dimensional Gaussian
chain.
Knowing the probability p(n, r) gives us information about the energy
needed to distort the polymer. We recall that any system with a given relative
probability f (r) had a free energy given by −T log f . We saw that this free
energy is equivalent to a potential energy: if r is changed from one definite
value to another the work required is the change in this free energy. For our
polymer, the free energy
3 T r2
U (r) = −T log(p(n, r)) = . (3.14)
2 r 2
The polymer stores energy like an ideal spring with zero unstretched length.
As with the spring, the force required to elongate is proportional to the
distance. In magnitude, the energy U is of order T for typical extensions
r
r 2 1/2 . This spring energy appears in many forms in polymer liquids. It
is responsible for the rod-climbing behavior of Fig. 1.6. It is also responsible
for the microphase separation patterns seen in block copolymers, like that
of Fig. 1.1.
3.5. Coil-Stretch Transition An elongational flow field may be made by aiming
two round nozzles at each other underwater, then sucking water out through
both nozzles at the same rate. Any polymers at the center of this contraption
are pulled towards both nozzles at once. The effect is roughly equivalent to a
force pulling on each end, of magnitude bz2 , where z is the distance from the
center along the nozzle axis. The coefficient b is proportional to the elongation
rate γ̇ . One can detect the amount of elongation z induced in the polymers by
the flow. (a) What energy W (z) would a free particle starting at the origin gain
if it moved to z under the action of this force? This energy is part of the energy
of any polymer ending at z in the presence of the flow. (b) Denote the z2 of
the polymer without flow as Z 2 . Treating the polymer as a spring, as explained
in this section, sketch the work U (z) required to elongate the chain to length
z as a function of z for small b, counting both the elastic free energy of the
polymer and energy in (a), and write its functional form. (c) This energy U (z)
attains a maximum for some elongation z∗ and energy U ∗ (relative to its z = 0
value). How does the height of this maximum vary with b, when b is small? If
the height U ∗ is adjusted to be the thermal energy T , what is the corresponding
position z∗ of this maximum in terms of Z? (We have seen that this maximum
amounts to an activation barrier. When it is higher than about T , passage over
the barrier becomes very slow.) (d) If this polymer were in thermal equilibrium
at temperature T in the presence of this energy U (z), find z to lowest order
in b. In this case U ∗ T and |W (z)| T for all z’s that matter. According to
this U (z), once z exceeds z∗ , it will increase to infinity. In reality z increases
only to the maximum extension rmax of the polymer. This maximum extension
is not accounted for in the Gaussian p(n, r), as the text points out.
3.6. Series and parallel polymers Two identical random-walk polymers are
fastened together end to end, to make a chain of double the length of each.
(a) What is the root-mean-squared end-to-end distance r 2 of the two
free ends relative to that of a single polymer? (b) The same two polymers
54 Polymer molecules
are fastened together at both ends to make a double chain. What is the
root-mean-square end-to-end distance between the two ends in this parallel
configuration? Hint: One may construct the ensemble of double chains by
letting the two chains end anywhere, and then discarding all configurations
except those in which the two chains end at the same point. By this proced-
ure, one can relate the p(n, r) of the double chain to that of the two single
chains.
∞
−3/2 r2
ρ(r) 0 = 2 di i f . (3.15)
0 i
Interior structure 55
The integral is finite: at small ĩ, the integrand is exponentially small because
f is small for large argument; at large ĩ f is a constant, but the integrand
still falls off as an integrable power of ĩ. The integral is simply a constant
independent of r, so that ρ(r) 0 = (const)/r.
This power-law dependence of the local density means that our random-
walk polymer is a fractal set† in the sense of Mandelbrot [6]. That is, †That is, the polymer approaches the fractal
the average mass (number of monomers) M(r) within distance r of an behavior as one approaches asymptotic con-
ditions of large chains and long distances.
arbitrary point of the set varies as some power r D . Then this power is
called the fractal dimension. This M(r) can be found for the polymer by
integrating the local density: M(r) = r <r d 3 r ρ(r ) 0 = (const)r 2 .
A random-walk polymer has a fractal dimension of two, like a uniform
surface. Clearly, this property generalizes: any fractal with dimension D
has a local density ρ(R) 0 that falls off as R D−d in any spatial dimension
d > D. This is discussed in Problem 3.7.
3.7. Why the fractal dimension D is called a dimension It was said in this section
that in a fractal object, the number of points of the object within a given distance
R of a given point grows as R D . Show that this property holds for a straight
line or a plane in three-dimensional space. (a) What is the fractal dimension D
for these two objects? (b) What happens to the two D values when the line and
the plane are embedded in four dimensions?
3.3.1 Scattering
The behavior of the local density ρ(r) 0 is directly observable by
† The subject of scattering from condensed scattering† . We saw in Chapter 1 that scattering senses the size of an object:
matter is discussed in detail in e.g. Stephen small objects diffract the incoming waves like a slit diffracts light. Scat-
W. Lovesey Theory of neutron scattering
tering senses not only the overall size of a structure but also the internal
from condensed matter (Oxford, Clarendon
Press, 1984). The book by Jannink and Des distribution of matter within it, as we now show.
Cloizeaux mentioned at the beginning of In a scattering experiment a plane wave or beam of scattering particles
this chapter also has an extensive discussion (photons, neutrons, etc.) impinges on an unknown object, such as liquid con-
of neutron and light scattering. taining polymers. For definiteness we will imagine that our beam particles
are neutrons: the waves are then complex scalar Schrödinger waves. The
beam wave has the form A exp(i k · r − iωt). The constant A is the amp-
litude of the wave; its square is the number of neutrons per unit volume in
the beam. The angular frequency ω is proportional to the kinetic energy of
the neutrons. The wave vector k is proportional to their momentum. Evid-
ently, one can change k and ω together by varying the speed of the neutron
beam. But one cannot vary k and ω independently.
When this wave encounters an atomic nucleus, such as a carbon nucleus in
our solvent liquid, the great majority of it continues unperturbed. But a tiny
fraction radiates outward from the nucleus in a spherical wave of the form
B exp(i |k| r)/r. Of course doubling the number of neutrons in the beam
doubles the number scattered in this spherical wave. Thus B is proportional
to A; the proportionality factor is the (complex) scattering amplitude f . A
detector receives the wave far away in some particular direction (Fig. 3.6).
The part of this spherical wave entering the detector is a plane wave with
some wave vector k of magnitude k. The detector measures the intensity of
the wave as a function of the direction of k . The r −1 factor in the scattered
wave is virtually the same distance for all the scatterers, since the detector
is far away from the sample.
The detector senses the scattered waves from all the M scatterers j
in a given volume. The amplitude from a particular scatterer at rj is
B/r exp(k ·(r −rj )). As we have noted, the prefactor B must be proportional
to the wave function of the incident wave: B = f A exp(i k · rj − iωt). The
wave entering the detector is the sum of these contributions. The intensity I
is proportional to the absolute square of this wave. Since the scatterers in our
Interior structure 57
k⬘
Fig. 3.6
Schematic view of neutron scattering from
a solution containing a polymer. The
detector is shown and the wavevectors k
and k are indicated.
The quantity ρρ is called the density correlation function. This ρρ is
simply related to the local density defined in the previous section. We may
think of it as the joint probability that a monomer is at the origin and that
another is at r. This is the overall probability that a monomer is at the origin
times the conditional probability that a second monomer is at r given that
one is at 0. The probability that a monomer
is at the origin (or any other
point in the sample) is ρ = M/( d 3 r). The conditional probability is
† The density correlation function also con- the same as the local density ρ(r) 0 . Thus ρ(r )ρ(0) = ρ ρ(r) 0 † . The
tains the same information as the pair correl- relation between S(q) and the local density is simple:
ation function g(r) introduced in Chapter 2.
We recall that g(r) means the average dens-
ity of particles at displacement r from S(q) = d 3 rρ(r) 0 exp(i q · r). (3.21)
a particle, relative to the average dens-
ity ρ . Thus the connection is g(r) =
ρ(r) 0 /ρ = ρ(r)ρ(0) /ρ 2 . Of the The reduced scattering intensity S(q) is related to the local density by a
three quantities, the local density is most simple integral transform—the Fourier transform. This means that if one
convenient for discussing fractal structure,
expresses the scattering density ρ(r ) as a sum of plane waves of all different
because it avoids unimportant factors of
average density. wavevectors, the scattering intensity and S(q) arise from those waves that
have wavevector q. If the local density is constant the S(q) is simply zero
for q > 0.
Interior structure 59
That is, the internal structure is completely invisible to S(q) and the scat-
tering is about the same as though all the scatterers were concentrated at
a point. Our polymers show this behavior like any other large object. It is
straightforward to see how S(q) begins to vary when q grows comparable
to 1/R; this is discussed in an exercise.
In the complementary limit where qR 1 the S(q) senses the internal
structure of the object. If this object is a fractal made of units of size a,
as our polymer is, we may readily infer how it scatters in this limit. As
discussed above, a fractal of dimension D in three-dimensional space has a
local density of the form ρ(r) 0 = c r D−3 for a r R. For such power
laws it is easy to perform the integral in Eq. (3.21) by defining r̃ ≡ |q|r .
This yields
|q|R
−D
S(q) = c|q| d 3 r̃ r̃ D−3 exp(i q · r̃/|q|). (3.23)
The scattering is the same as though the fractal were chopped up into pieces
the size of a scattering wavelength.
1.0 Figure 3.7 summarizes the expected scattering from a fractal. As q
S (q) increases from 1/R to the inverse scatterer size 1/a, the scattering intensity
M drops by a factor of M, the number of scatterers in the fractal. Remarkably,
the S(q) from a random-walk polymer can be calculated explicitly as a
0.1 –D
check on these ideas.
3.8. Scattering from a random-walk polymer It was shown in the text that
the structure factor S(q) of a set of n scatterers at positions rj is
0.01 1/n nj,k=1 exp[iq · (rj − rk )] . Recall that the probability distribution p(i, r)
0.1 1.0 10.0
q RG
for a segment of length i to have its ends at displacement r can be expressed
in the product form with three factors for the x, y, and z directions:
p(i, r) = A(i) exp[−x 2 /(2ix 2 1)] × [x → y] × [x → z], where the pref-
Fig. 3.7 actor A assures normalization: r p(i, r) = 1. (a) From this and the fact
Schematic log–log plot of S(q) versus q
that dx exp(−x 2 ) exp(iqx) = π 1/2 exp(−q 2 /4), derive an expression for
for a fractal set of M 1 scatterers with S(q) in closed form. Note that the number of monomer pairs separated by
fractal dimension D. (In fact, this is the a distance m monomers is not the same for all m. There is only one such
exact structure function for a random-walk pair for m = n, but there are n − 1 pairs when m = 1. If one has a suf-
polymer, as requested in Problem 3.7. The ficient range of q accessible experimentally, one may discern two scaling
wave vector q is given in the combination
regimes in this S(q): for q 2 x 2 n 1, S(q) ∼ q 0 ; for q 2 x 2 n 1,
qRG , where RG is the radius of gyration
discussed in Problem 3.11.) S(q) ∼ q −2 . To get a feeling of what range of q is needed, plot log(S(q))
1/2 1/2
versus log(q) for (b) 0.3 < qx 2 n < 3, (c) 0.03 < qx 2 n < 30,
1/2
(d) 0.001 < qx n < 1000. The scales should be adjusted so that the three
2
plots are of roughly equal size. What range of q would be needed in your judg-
ment to discriminate between the expected q −2 scaling and a hypothetical q −1.8
scaling?
3.9. Scattering from a line The text shows that the scattering intensity at large
wave vector q from a fractal object was proportional to k −D . (a) Find explicitly
the structure factor S(q) from a line of length 2L oriented along the z-axis
between −L and L. The scattering wave vector q makes an arbitrary angle θ
with the line. (b) For a line at random orientation, show that the scattering is
proportional to q −1 for large q, as it should be for a one-dimensional fractal.
Note
1 that in three dimensions an angular average of a function f (θ) amounts
to 0 d(cos θ )f (θ )/(4π ).
3.10. Scattering from a Gaussian Cloud Since the monomers in a polymer chain
are randomly arranged throughout a volume of limited size R, it is sugges-
tive to approximate the polymer as a mere cloud of n monomeric scatterers
confined to a region of this size. A simple realization would be a Gaussian
cloud, whose local density ρ(r) is given by nCR −3 exp(−(r/R)2 ). How
does the scattering intensity from such a cloud vary with wave vector
q? Does it have the proper fractal behavior at large q? A useful integral
∞
is −∞ dx exp(iqx − x 2 ) = π 1/2 exp(−q 2 /4).
3.11. Radius of gyration The structure factor S(q) of any object composed of
n scatterers may be expanded in powers of the wave vector q: S(q) =
n(1 + aq + bq 2 + · · · ). (a) What are the coefficients a and b, expressed in
terms of the positions rj ? Assume that the scatterers are isotropically arranged
so that the scattering is the same in all directions. Use the power series expan-
sion for ex : ex = 1 + x + 12 x 2 + · · · . (b) The radius of gyration of any set
of n particles is defined as
the root-mean-square distance between an arbitrary
pair of them: RG 2 ≡ n−2 1 rj − rk )2 . Find an expression for RG in terms
2 j k (
of the scattering coefficients a and b. (c) For a random-walk polymer, how is
the radius of gyration related to the mean-squared end-to-end distance r 2 n ?
Self-avoidance and self-interaction 61
Since we know the gaussian form of p(n, r), this integral can be worked
out explicitly (along with the needed normalization). Since p(n, r) falls
off faster than any power for large r, all moments r a 1/a are of order
1/2
r 2 RW . (The reasoning is the same as for the P (M) distribution discussed
in Appendix A.) Thus any power-law dependence of the remaining frac-
tion t would leave the average size of the self-avoiding chains basically
unaffected. They would increase by only a finite factor. Conversely, if
the self-avoiding chains are to be arbitrarily larger than the random-walk
chains, t must vary more strongly than a power law. Clearly, the overall
fraction remaining from the initial random walks must be indefinitely small.
62 Polymer molecules
Our aim in this section is to understand what types of contacts are important
and how much they alter the shape of the polymer.
dimensions. Thus,
MAB ∼ d d rr DA −d (r − x)DB −d . (3.25)
When all the blocks of a chain swell by a factor b the overall size of the
chain must increase by the same factor, since the chain is still a random
walk of the expanded blocks. We have argued above that the expansion
factor b must be finite irrespective of the number of monomers involved.
Thus the expansion factors for all stages k of our process are essentially the
† The initial blocks, since they involve short same† . Therefore the self-avoiding size RSA expands by a factor of order
segments, may have different expansion bk relative to the random-walk size, RRW . Recalling that RRW ∼ n1/2
factors b from the majority. But the longer
and noting that bk = 2k log2 b = nlog2 b , we conclude that RSA ∼ nν ,
the chain is and the more stages k of subdivi-
sion are present, the less these initial blocks where ν = 12 + log2 b is a power larger than 12 . This argument says that
can matter. self-avoidance alters the scaling relationship between chain length and
size. The exponent ν describing the new scaling is called the Flory swelling
exponent.
In random-walk polymers, the scaling of the overall size R reflected a
pervasive dilation symmetry for the internal mass distribution. The same
thing proves to be true of self-avoiding polymers. First, the end-to-end
distribution function p(n, r) should have a scaling property analogous to
that of the random walk, since it represents a limit of indefinitely large n.
††The radius of gyration grows by a factor The scaling form must give average end-to-end distances proportional to
of about 10%, compared with the same nν , and must preserve normalization. Thus p(n, r) = n−νd f (r/nν ). If we
chain without tails: cf. [8], Chapter XV.4,
consider a finite piece of length n within an infinite self-avoiding polymer,
Table 8.
its mean size is of the same order as that same piece in isolation. That
is, the attachment of an infinite tail to each end of the piece extends it
†††It is natural to ask whether a self- by only a finite factor [8]†† , as with the finite tails discussed above. This
avoiding polymer has the full dilation sym- scaling law can be used to infer the local density ρ(r) 0 , in the same
metry of the random walk, as discussed in way done for the simple random-walk polymer in Eq. (3.16). The result
Appendix A. The answer is yes, though the is ρ(r) 0 = (constant) r 1/ν−3 . Since ρ(r) 0 falls off as a power of r,
proof is too complicated to give here [8]:
M(r) is also a power of r. Repeating our reasoning for the random walk,
ρ(r1 ), . . . , ρ(rk ) 0 we find M(r) = (constant) r 1/ν . Thus the self-avoiding polymer is a
fractal whose fractal dimension D = 1/ν ††† .
= λ−k(D−d) ρ(λr1 ), . . . , ρ(λrk ) 0 .
the needed expansion by pulling on the chain ends. Then the work S
required is the free energy of extension, calculated in the last chapter:
S = 32 T R 2 /RRW
2 = 32 R 2 /(a 2 n), where a is the size of a monomer. Such
expansion is favored because it reduces the number of contacts. To account
for these contacts, we imagine that contacts are allowed, but at a cost of a
small energy U per contact. That is, we replace the mutually self-avoiding
monomers by weakly self-repelling ones. As we have seen, local contacts,
despite their large number, are not important for the chain expansion. The
number of nonlocal contacts may be found by imagining that the chain is
a dispersed cloud of n monomers confined to the region of size R. The
number of contacts is n times the probability that a given monomer has
a contact. This probability is the monomer volume v times the monomer
density n/R d . Combining these facts, the overall energy V due to these
contacts is given by V
vU n2 /R d . Asymptotically (if d < 4 and the
chain is opaque to itself) there are indefinitely many contacts and the con-
tact energy grows much larger than the thermal energy T ; this is enough
free energy to account for substantial elastic stretching.
self-repulsion. (The invisible beads are freely allowed to sit at occupied sites. This
10–1
reduces the average b2 by a factor of 4.) Observe the resulting changes in M(r)/r 2 .
Determine how much the asymptotic D changes as a result. Compare the radius R I –1(q)
10–2 1.67
with the case of full self-avoidance and compare with the prediction of the Flory 2.00
formula.
10–3
The polymer is free to choose its size R. By the above estimates, the A
potential of mean force associated with size R is S + V . We thus expect
10–4 B
that the typical size will be that which minimizes S + V . Since S and V are
both powers of R, the minimum occurs when S
V : C
10–5
10–2 2 5
10–1
R 2 /n ∼ U n2 /R d , or R 2+d ∼ U n3 , or n ∼ R (2+d)/3 . (3.27) q Å–1
We are led to the conclusion that for sufficiently large n, the fractal Fig. 3.9
dimension Log–log plot of the inverse scattered
intensity as a function of the scattered wave
D = (2 + d)/3. (3.28) vector q, confirming the internal fractal
structure, after [11, Figure 1], © American
Physical Society, courtesy G. Jannink. The
In three dimensions D = 53 . We note that the D value is independent of the full points are experimental data for a
strength U assumed for the contact interaction: for sufficiently long chains, polymer in a good solvent, curve A, a theta
even an arbitrarily small U still leads to a strong expansion. solvent, curve B, and the bulk polymer,
This prediction for D is consistent with estimates based on systematic curve C. The behavior of curves B and C
will be discussed later in the chapter. The
calculations, experiments and computer simulations. One experimental polymers used were polystyrene and
confirmation is shown in Fig. 3.9. Moreover, in two and in four dimen- deuterated polystyrene with molecular
sions, D is known exactly and it agrees with the Flory prediction. Despite weight about 106 . The solvents were
the success of the Flory argument, it has not proved possible to justify it† . carbon disulfide (curve A), deuterated
Because of the success of the Flory formula, efforts continue to find a simple cyclohexane at 40◦ C (curve B), and
undeuterated polystyrene (curve C). The
way to understand the scaling of self-avoiding polymers [12]. dashed lines are reference straight lines
3.13. A lower bound on D for a self-avoiding polymer According to the Flory with slopes of 2.00 and 53 showing the
argument, the constraint of self-avoidance lowers the fractal dimension D. power-law behavior expected for D = 2
But it must not be too much lower: if D were too low, the expanded polymer and D = 53 . The full lines are calculated
would not be opaque to itself, and the contact energy would be insignificant. curves incorporating the expected scaling.
(a) From this consideration, give a lower limit on D in d dimensions [13]. Is
this bound consistent with the Flory formula for D? †
3.14. Triple contacts In some situations polymers obey a modified type of The Flory reasoning can easily lead one to
self-avoidance: intersections of two monomers (of radius a) are freely wrong conclusions. For example, the argu-
allowed, but intersections involving three monomers are forbidden. (a) For ment appears to be equally justified for a
n monomers in a region of radius R how does the number of triple intersec- chain where only the first and last third repel
tions vary with n and R. Your estimate should be in the same spirit as the each other. Yet we have shown that such
Flory estimate for double contacts cited above. (b) Assuming that each triple avoidance can only expand a chain by a
contact costs a given energy U , apply the Flory reasoning to estimate the factor of order unity.
fractal dimension of the chain in three dimensions, and in two dimensions.
osmotic pressure encountered in the last chapter. The volume within r0 will
be denoted v0 . To exclude segments completely from volume v0 would
require work v0 ; to exclude a fraction (1 − g(r0 )) requires a proportion-
ate amount of work, since the work is proportional to the number excluded:
(1−g(r0 ))v0 . This process has produced the correct g(r) for one segment
out to a distance r0 . We may extend it to produce a correct g(r) out to some
slightly larger radius r1 . The new volume affected will be denoted v1 . The
associated work is evidently (1 − g(r1 ))v1 . Continuing in this way, we
may impose the full g(r) between the chosen segment and the others, with
an associated work
[T (n − 1)/ ] vi (1 − g(ri )) = [T (n − 1)/ ] d 3 r(1 − g(r)). (3.29)
i
(a) Use this expression for F (R) to find the work to confine as a specific
function of R. This estimate gives the correct scaling of the confinement
energy, agreeing with more rigorous arguments. It also can be extended
to the confinement of self-avoiding polymers, as the rigorous arguments
cannot.
(b) Compare the work of confinement found in (a) to the interaction energy
of a collapsed chain discussed in the text. Is it important or not for long
chains?
3.4.5 Polyelectrolytes
It was mentioned in the first section that polymers in solution can be made
with an ionic charge every few bonds along the chain. The electrostatic
repulsion between these ions tends to expand the chain, and the degree of
expansion proves to be more extreme than that of neutral chains. To see
this, we imagine turning on the electrostatic repulsion from the most global
progressively down the most local scales, inverting the procedure of the last
section. We begin by considering the repulsion only between the first half of
the chain and the second half. When this is done for a neutral chain, it causes
the two halves to swing away from each other; this reduces the number of
contacts between the halves, and thus its interaction energy, dramatically.
The result of this global interaction by itself is to increase the radius by a
factor of order unity. At that point the energy of interaction between the two
halves is of order of the thermal energy T † . †For this statement to be strictly true, we
When we do the same thing with a polyelectrolyte chain, we must must also turn off the interactions involving
monomers near the center one.
consider the electrostatic repulsion between the first and last halves. If
the entire chain has n monomers, each bearing, say, one charge, the repul-
sion is roughly that between two charges of magnitude en/2 at a separation
R as large as the polymer. This energy is roughly e2 (n/2)2 /(R), where
is the dielectric constant of the solvent. It may be expressed in units of
the thermal energy T as T (n/2)2 (/R). The Bjerrum length ≡ e2 /(T )
is roughly 0.7 nm in water at room temperature. If R were the size of a
random-walk chain, so that R ∼ n1/2 , this energy clearly grows indefin-
itely large with n. Unlike the short-ranged energy, it cannot be much reduced
by simply making the two halves swing out of each other’s way. This global
energy would cause a large expansion of the chain—not just a finite-factor
expansion as in the neutral case. Balancing this global Coulomb repulsion
against the elastic energy T R 2 /(a 2 n), the two energies are comparable
when R ∼ n. We are led to the conclusion that the swelling exponent
ν = 1; this is as large as ν can be for a connected chain. The repulsion
is strong enough to stretch the chain to within a finite factor of its full
extension, as though the chain were a rigid rod.
To confirm this picture and to see how rigid this rod is, we are led to model
the chain as a row of evenly spaced charges at separation R/n, connected
1/2
by random-walk polymers of unstretched length r 2 1 ≡ a. The preferred
distance between charges can be found by considering the tension in the
polymer near the center. This tension F is equal to the repulsive force
between the left and the right halves. Adding up the forces on each charge
on the left half, one finds†† †† Here we have approximated the sum over
random-walk chains store energy like an ideal spring, they exert a force
proportional to their elongation. Specifically, F = 3T (R/n)/a 2 . Equating
these tensions leads to the preferred spacing: R/n
a( log(n)/a)1/3 . If the
charges were placed sparsely along the chain, a , and the spacing R/n
would be much smaller than the unperturbed size a. Thus locally, the effect
of the charges can be weak. Still, the overall effect is strong. Even a weakly
charged chain, if long enough, would stretch the monomers out to their
full extension, because of the logarithmic factor. We shall assume in what
follows, though, that the charges are sparse enough so that the monomers
are not strongly distorted: R/n a.
The real chain must fluctuate around this regular rod state by stretching
and by bending. If the chain stretches by a small amount γ R, the asso-
ciated energy must be proportional to the square of this small change. It
is convenient to express the constant of proportionality as ER so that the
stretching energy is ERγ 2 . To estimate the modulus E we consider strains γ
of order unity, so that the length R is doubled. Such a stretching must change
both Coulomb and elastic energies by factors of order unity; hence the total
energy is changed by such a factor. This means the stretching energy is of
the order of the original Coulomb energy: ERγ 2 |γ =1
nT ( log(n)/a)2/3 ,
or E
a −1 ( log(n)/a)1/3 . Except for the weakly varying log(n) factor,
the modulus E is independent of n, as in a solid, macroscopic rod. One
may now estimate the amount of strain γT typically present because of
thermal fluctuations by setting ERγT2
T . Thus γT ∼ 1/R 1/2 if we
again neglect the log n so that E is independent of size. The log(n) factor
only reduces γT further. As R becomes large, the thermal fluctuations in R
become increasingly insignificant.
R A polyelectrolyte is also rigid with respect to bending, e.g., into an arc
of radius B. One can readily find the Coulomb energy needed to bend a
uniformly charged rod in this way. Any small bending must cost an energy
of the form ẼR(R/B)2 . But a bend of, say, 90◦ with B
R changes
the distance between two arbitrary charges by a factor of order unity as
B illustrated in Fig. 3.11. Thus the Coulomb energy increases by a similar
factor: ẼR
n2 T /R. The coefficient Ẽ is of the same order of mag-
nitude (except for the log(n)) as the modulus E. The amount of thermal
Fig. 3.11 bending R/BT , like the thermal strain γT , falls off as the square root of the
A charged rod of length R is bent into an chain length n. We conclude that a polyelectrolyte chain has little of the
arc of radius of curvature B
R, thus fluctuations that random-walk or self-avoiding chains have. The assumed
shortening the distance between a typical state of a uniformly charged line is affected but little by thermal fluctuations
pair of charges (heavily drawn) by a finite
factor.
in stretching or bending.
In this chapter we have seen that connecting arbitrary small molecules to
form a flexible chain leads to remarkable structures. These structures are
very different from their compact constituents. They are indefinitely tenu-
ous and deformable. In order to think about these objects quantitatively, a
new concept is needed. This is the notion of fractal scaling. Using it, we
found that chains made of different constituents nevertheless converge to a
common quantitative behavior, provided they are long enough. This com-
mon asymptotic behavior is embodied in scaling laws and scaling functions
like the reduced scattering function S̃(q) defined above.
Appendix A: Dilation symmetry 75
In the next chapter we exploit and extend these scaling laws to explain the
striking properties of polymer solutions: their springiness, their viscosity,
and their tendency to absorb further solvent.
3.16. Screened polyelectrolyte: the Odijk length In our discussion of polyelec-
trolytes we ignored the countercharge that neutralized the ions on the chain
before it was put into solution. The interaction between the polyelectrolyte
charges is modified by the countercharges. In addition, the solvent always
contains further ions, not related to the polymers. We shall see in Chapter 5
that the effect of these ions is to change the interaction potential between two
unit charges to the form (e2 /r) exp(−κr): the potential falls off exponentially
for distances larger than the screening length 1/κ. This κ depends only on the
concentration of ions and is independent of the chain length n. Thus when
n is large, the interaction is relatively short-ranged: κR 1. Still, even such
screened chains are remarkably rigid. By extending the methods sketched in
the text, Odijk [16] found that that screened chains can be rigid on length
scales much greater than the screening length 1/κ.
Consider a length R of chain much longer than 1/κ, bent into a gentle arc of
radius B R. As in the text, the chain is treated as a line of uniformly spaced
charges separated by a distance a. (a) Find the change in Coulomb energy of a
given charge in the limit of weak bending. It must be of order 1/B 2 . The total
energy cost of bending the rod is evidently this energy per charge times the
number of charges n. For sufficiently long R, the cost of bending the chain
through a unit angle, with B = R, becomes of order T . This length R is called
the electrostatic persistence length, or Odijk length L0 . Chains shorter than the
Odijk length behave as rigid rods, even with screening; chains much longer
than the Odijk length coil randomly. (b) The quantity L0 κ is the persistence
length relative to the screening length. In a limit where the distance between
charges a is much smaller than the screening length 1/κ, how does L0 κ scale
with κa? Is the chain flexible or rigid on the scale of 1/κ?
Fig. 3.12
Eight randomly sampled images of a long
self-avoiding random-walk polymer. (We
show self-avoiding polymers because
simple random-walk polymers are difficult
to illustrate in two dimensions.) Half of the
images are at a 2.5-fold magnification
relative to the other half. The monomer
positions are blurred to a fixed fraction of
the picture size. In each case the picture
was centered on a monomer. The two
magnifications are interspersed randomly.
The reader is invited to guess which are the
magnified images. (Polymers were
generated via a Monte-Carlo simulation
originally written by Betsy Weatherhead
and drawn by Olgica Bakajin, 1995) For
the solution see the end of this section.
blurry line of the same darkness passing through the center. There is no
way to distinguish which has the larger magnification. An example lacking
dilation symmetry would be a soap froth in a pan, sampled by preparing the
froth repeatedly and then photographing it. Pictures at different magnifica-
tions are clearly distinguishable because the typical bubbles have different
sizes at different magnifications.
To state this invariance mathematically, we may define a set of statistics
that describe the random density ρ(r). One way to characterize the vari-
ations of ρ over the images is by means of the correlation functions
ρ(r1 )ρ(r2 ) . . . ρ(rk ) 0 . The 0 subscript indicates that the pictures in our
sample are selected to have the origin occupied. The local density discussed
in the main text is the simplest example of these correlation functions. The
points r1 , . . . , rk are specific points on each image. For example, r1 could
be 1 cm inward from the upper left corner. Any of these correlation func-
tions can be measured to arbitrary precision using a large enough sample of
images. The same process performed on images magnified by a factor λ−1
gives ρ(λr1 )ρ(λr2 ) · · · ρ(λrk ) 0 . Dilation invariance implies that the two
statistics are the same:
each segment between these points. Of these cases, we consider first the one
where the points are visited in order: the chain goes from the origin to r1 ,
then to r2 , etc. We further suppose that the number of monomers between
0 and r1 is i1 , and that the succeeding segments have i2 , . . . , ik monomers.
Once all this has been specified, the probability of such a chain can easily
be stated: it is the product p(i1 , r1 )p(i2 , r2 − r1 ) · · · p(ik , rk − rk−1 ). The
probability takes this simple form because each of the chain segments is
statistically independent. To find the joint probability for all chains that pass
through the points in this order, we must simply add all the contributions
from the different possible i1 , etc.:
⎡ ⎤ ⎡ ⎤
⎣ p(i1 , r1 )⎦ [ ] · · · ⎣ p(ik , rk − rk−1 )⎦ .
i1 ik
This sum is itself dilation invariant. To see this we use the scaling form
p(i, r) = i −3/2 f (r 2 /i), and replace the monomer sums by integrals. The
same sum dilated by a factor λ is
−3/2 2 −3/2 2
di1 i1 f ((λr1 ) /i1 ) [ ] · · · dik ik f (ik , (λrk − λrk−1 ) /ik ) .
(3.35)
Except for the λ−k prefactor, this is the same as the undilated sum. Thus
the sum is dilation invariant. We may repeat the reasoning for any order
of passage among the points 0, . . . , rk . All the corresponding sums have
dilation symmetry with the same prefactor. To obtain the overall probability
of passage through these points, we must simply add these sums together
for all possible orders of passage. The total evidently has the same dilation
symmetry as each part. Finally the density of monomers at r is evidently the
same as the probability density for a monomer to be there. That means that
our joint probability is simply ρ(r1 ) · · · ρ(rk ) 0 . We have thus shown that
Comparing this result with the general definition in Eq. (3.34), we see that
the μ prefactor takes a simple and suggestive form: it is a power of λ which
is proportional to the number of density factors ρ. The density behaves as
though it had dimensions of length−1 rather than its standard dimensions
of length−3 . This latter is the scaling we would find if we multiplied all
lengths, including the size of a monomer, by a factor λ. By contrast, our
78 Polymer molecules
where M̃ ≡ M/M . Since the moments scale the same, they are all finite
for a given r, however large. Thus the P (M̃) must fall off faster than any
Appendix B: Polymeric solvents and screening 79
log(具M典 P (M ))
off exponentially with M̃. The distribution also falls off for small M̃, as
exp(−1/M̃). It thus has the form sketched in Fig. 3.13. The function P (M̃)
must also be universal, since it can be calculated from the universal p(n, r).
Compared to a solid body, the fluctuations in M(R) for our random-
walk polymer are large as Fig. 3.13 suggests. The relative width of the
–2 0 2
distribution, e.g., (M − M )2 /M 2 , does not grow with r, but neither log M/具M典
does it go to zero† . Real, good-solvent polymers are dilation-symmetric
[8],with μk (λ) = (λ3−D )k . Their P (M) distributions are thought to be Fig. 3.13
universal and to resemble Fig. 3.13 [17]. Other fractal objects, such as Qualitative sketch of the mass distribution
colloidal aggregates are also thought to have this form of dilation symmetry. P (M) as a log–log plot. The M variable
has been scaled with M ; the P variable
has been scaled with M −1 . With this
Appendix B: Polymeric solvents and screening scaling the curve is universal. The
downward curvature for large and for small
In this Appendix we discuss how linking solvent molecules together results M indicates that P (M) falls off faster than
in a weakening of interaction between segments of a polymer dissolved any power for large or for small M.
in that solvent. This is easiest to understand if we place the solvent and
polymer on a lattice, with each segment occupying a lattice site. The solvent † Some known sets of points called multi-
molecules, being larger than the segments, occupy k > 1 sites within some fractals have scale invariance of a more
radius Rs . For concreteness, we can consider the k-site solvent molecules to general type. An example of a multifractal
be small polymers of the same species as the large polymer being studied. is the charge distribution on a conduct-
Thus in principle all the issues of fractal structure arise for both the solvent ing fractal that has been charged to some
potential V . Multifractal densities have μ
and the polymer, as discussed in the chapter. However, the weakening of powers that are not simply proportional to
interactions we wish to explain is independent of the details of this fractal the number of ρ factors. Thus multifractals
structure. Accordingly, we will ignore questions of fractal structure here. have mass moments that diverge at differ-
At first glance, our change to k-site solvent molecules seems unimport- ent rates as r → ∞. For a multifractal one
ant. For example, if these solvent molecules occupied a compact cube of must determine individually the scaling of
all these moments in order to characterize
k sites, our system would amount to a change to a coarser lattice whose the fluctuations of M.
sites had a linear size of k 1/3 instead of 1. But we wish to consider tenuous
solvent molecules that can interpenetrate. The k sites of each molecule per-
vade a volume shared by many other such molecules. This is the situation
when the solvent itself consists of polymers. But it is equally true when the
solvent consists of any tenuous structures, whose volume fraction decreases
indefinitely as k increases.
150x150 Polymer 2
150x150 Polymer 3
60x60 Polymer 1
60x60 Polymer 4
150x150 Polymer 4
150x150 Polymer 1
60x60 Polymer 2
60x60 Polymer 3
Fig. 3.14
The polymers of Fig. 3.12 without
blurring.
80 Polymer molecules
reduction factor applies roughly for typical points r in the pervaded volume:
1 − g(r)|r∈V
W /T .
Now we can find the mutual excluded volume b̃2 of our pair of molecules:
1
b̃2 = (1 − g(r))
W /T V . (3.41)
2 r
Using our estimate for W , we find b̃2
k. Using this result, we can infer
the self-interaction energy of an n-site polymer, containing n/k of the k-site
molecules. Repeating the reasoning of Eq. (3.31), the self-interaction energy
V of the polymer is given by
References
1. See e.g., D. W. Van Krevelen, Properties of Polymers: Their Correlation with
Chemical Structure; Their Numerical Estimation and Prediction from Additive
Group Contributions, 3rd ed. (Amsterdam: Elsevier, 1990) p. 333.
2. P. Rempp, E. Franta, and J.-E. Herz, in Polysiloxane Copolymers/Anionic Poly-
merization, Advances in Polymers Science series, no. 86 (Heidelberg: Springer,
1986), p. 147
3. G. G. Odian, Principles of Polymerization, 3rd ed. (New York: Wiley, 1991).
4. K. P. McGrath, M. J. Fournier, T. L. Mason, and D. A. Tirrell, J. Am. Chem.
Soc. 114 727 (1992).
5. H-G. Elias, Macromolecules, 2nd ed. (New York: Plenum Press, 1984), p. 124.
6. B. B. Mandelbrot, The Fractal Geometry of Nature (San Francisco: Freeman,
1982).
7. T. A. Witten, Rev. Mod. Phys. 70 1531 (1998).
8. G. Jannink and J. Des Cloizeaux, Polymers in Solution (Oxford: Oxford, 1992).
9. J. C. Le Guillou and J. Zinn-Justin, Phys. Rev. Lett. 39 95 (1977).
10. P. Flory, Principles of Polymer Chemistry (Ithaca, New York: Cornell, 1971),
Chap. XII.
11. J. P. Cotton, D. Decker, B. Farnoux, G. Jannink, R. Ober, and C. Picot, Phys.
Rev. Lett. 32 1170 (1974)
12. M. A. Moore and A. J. Bray, J. Phys. A: Math. Gen. 11 1353 (1978); S. F.
Edwards and P. Singh, J. Chem. Soc. Faraday Trans. 2, 75 1001 (1979); M. K.
Kosmas and K. F. Freed, J. Chem. Phys. 68 4878 (1978).
13. J. Des Cloizeaux, J. Phys. 31 715 (1970).
14. R. K. Pathria, Statistical mechanics (Oxford: Pergamon Press, 1972).
15. B. DuPlantier, Europhys. Lett 1 491 (1986); G. Jannink and J. Des Cloizeaux
op. cit., Chap. XIV.
16. T. Odijk, J. Polym. Sci. 15 477 (1977).
17. B. M. Sterner, to be published, demonstrates numerically the universality
of P (M) for two-dimensional self-avoiding polymers and obtains values for
several of the numerical coefficients Ck for positive and negative k.
18. M. Huggins, J. Phys. Chem. 46 151 (1942); M. Huggins, Ann. NY. Acad. Sci.
41 1 (1942); M. Huggins, J. Am. Chem. Soc. 64 1712 (1942).
19. S. F. Edwards, Proc. Phys. Soc. 88 265 (1966).
4
Polymer solutions
The previous chapter has shown how a polymer chain behaves when it is
by itself in a solvent. But the properties seen in a macroscopic fluid result
from many polymers interacting in the solvent. In this chapter we turn to
a survey of how polymer structure, energy, and motion are altered by their
mutual interactions. We first describe the energy associated with interac-
tion in dilute solution. The following section discusses the new structural
and energetic features that appear when the polymers interpenetrate. Then
we turn to dynamics. The first subsection discusses the Brownian motion
of a sphere, the starting point for understanding how objects move in a
liquid. We show how spontaneous thermal motions are related to the work
needed to create motion. We discuss the viscosity added to a liquid by
inserting particles or polymers into it. Next we consider how these proper-
ties are altered when polymers interpenetrate. We describe how flow fields
are modified by hydrodynamic screening. We distinguish between self-
diffusion and cooperative diffusion. Finally we discuss how entanglement
affects both diffusion and the relaxation of stress.
preceding chapter:
1
B2 = d 3 r(1 − gp (r)). (4.1)
2
The presence of V means that each increment of compression d requires
an extra work −(dV/d) d. Thus V adds an amount −dV/d to the
osmotic pressure :
2
dV NT N N
= 0 − = 1 + B2 + O . (4.2)
d
This formula for osmotic pressure is valid irrespective of the type of solute
molecule; the distinctive polymer features show up in how the excluded
volume B2 depends on molecular weight n and size R. We may infer how
B2 behaves by considering the pair correlation function gp (r). For r R,
two polymers are out of each others’ range of influence, so that gp (r)
should be nearly 1. But for r < R there is the possibility that gp will be
strongly influenced. To gauge the magnitude of such an influence, we con-
sider our polymers to be in a good solvent, i.e., one in which the monomers
have b2 > 0. We have seen that the effect of the b2 may be reproduced
by excluding all configurations in which two monomers approach closer
1/3
than some small distance of order b2 . To see the effect of this exclusion,
we take two polymer configurations at random, move them to a separation
r < R, and then discard the result if there are any intersections. The gp (r)
is the probability that the result is retained, i.e., the probability that there
are no intersections. As we saw in Eq. (3.26), the average number of inter-
† Here and in the following we shall use sections MAB (r) ∼ r 2D−d ∼ r 1/3† , and the probability of no contact is
the Flory value for D, keeping in mind significantly less than 1 for r < R. For example, gp (r) < 12 for all r smaller
that the actual number may well be slightly
than some fixed fraction αR as R → ∞. (This was a consequence of two
different.
polymers being mutually opaque.) This sets a lower bound on B2 :
1 3 1 1 4π 3
B2 = d r(1 − gp (r)) > (αR) . (4.3)
2 2 2 3
On the other hand B2 cannot grow faster than R 3 : polymers cannot exclude
each other from regions arbitrarily larger than their size. There must thus
be some constant C for which B2 < CR 3 ; this is an upper bound. In
view of these two bounds, B2 /R 3 must remain finite (and greater than 0)
as R → ∞. The polymers increase the pressure as though they were
mutually impenetrable spheres whose radius is roughly the geometric radius
R. This is a consequence of their mutual opacity.
4.1. Osmotic pressure of hard-sphere gas A dilute gas of hard spheres of radius
R has a g(r) which is zero for r < 2R and 1 beyond. (a) What is the excluded
volume B2 for such a gas? (b) For what volume fraction φ = N 43 π R 3 / of
spheres is the B2 term in the osmotic pressure equal to the ideal-gas term?
(Some amusing points of comparison: a close-packed crystal has a volume
fraction of 74%; random close packing has a volume fraction of 64%.)
An important further consequence of this opacity is implicit in the finite-
ness of B2 /R 3 as R → ∞. We have argued that it is finite whenever the
Semidilute solutions 85
this is called the volume fraction and is denoted φ. The overlap volume
fraction φ ∗ is clearly much smaller than the melt volume fraction—viz. 1.
Volume fraction is an unambiguous way to measure concentration, since it
does not depend on our definition of the monomer building block for the
chain. For typical hydrocarbon polymers, with molecular weight of 105 , φ ∗
is roughly 1%.
4.2.1 Structure
As the concentration rises above φ ∗ the spatial structure of the poly-
mers is significantly altered. We may express this structure in terms of
the local volume fraction analogous to the local density ρ 0 : φ(r) 0 ≡
ρ(r) 0 /ρmax . The r dependence of φ(r) 0 for r R can be written
(A/r)d−D ; the length A depends on the type of polymer and solvent but
not on its length n. Thus e.g., for polystyrene in benzene A
1 nm. If
φ φ ∗ , there is some distance r R for which φ(r) 0 = φ. This dis-
tance, called the correlation length ξφ = Aφ −1/(d−D) , plays an important
role in describing the structure of the solution. If we consider a region
around a given monomer much closer than this ξφ , the local density is
much greater than the average density. The average density can thus have
little effect on the self-avoidance interaction. We infer that the portion of
the polymer in this region is but little affected by the other chains. However,
for distances larger than ξφ , the original φ(r) 0 is much smaller than the
average volume fraction φ: the chain passing through the origin contributes
only a tiny fraction of the overall density. Thus the other chains may alter
the self-avoidance effects strongly. The segment of chain within a distance
ξφ is often called a blob. Each chain thus may be thought of a sequence of
blobs, each having an internal structure that is relatively unaffected by the
other chains.
To see how self-avoidance works on scales much larger than ξφ , we con-
sider global contacts between two blobs that are far apart along their chain.
As shown in the preceding section, the swelling of a chain depends on
the mutual excluded volume b2 of its constituent monomers. Here we may
take our monomers to be blobs. To find the b2 of two blobs, we imagine a
dilute admixture one-blob chains in the semidilute solution of k-blob chains.
Because interpenetration of the blobs is significantly suppressed, this situ-
ation is similar to a gas of monomers in a lattice, with k-site molecules
filling all the other sites. We discussed the effect of such multi-site solvents
in Appendix B of the preceding chapter. There we found that the extended
solvent molecules reduce the b2 between single sites by a factor of k.
To gauge the importance of this reduced b2 , we consider the number
of effective contacts within a given chain in the absence of swelling. This
chain has a radius R of order ξ k 1/2 and a volume of order (ξ k 1/2 )3 . Each
of the k blobs has an effective intersection with roughly kb2 /R 3 other
blobs of the same chain. Thus the number of effective intersections is
k 2 b2 /R 3 . Accounting for the k dependence of R and b2 , this gives an
effective number of contacts of k 2−1−3/2 ∼ k −1/2 . The effective interac-
tion energy between the monomers is thus of order T k −1/2 —much smaller
than the thermal energy. The interaction energy is not sufficient to expand
Semidilute solutions 87
the polymer significantly, since such expansion requires work of the order 1
of T . We infer that a chain has the structure of an unperturbed random
walk beyond the scale ξφ : it is like a random walk of blobs. The local
log ((r))
– 4/3
density thus behaves as shown in Fig. 4.2. The density from the chain at
the origin falls off as in a self-avoiding walk within the blob radius, and
like a simple random walk beyond. Monomers from other chains have a
+0.8
suppressed local density near the origin. These dominate at distances bey- –1
ond the blob radius. The suppression as a function of distance also obeys a a log(r)
power law ([1], chapter 13). Calculations suggest that its numerical value
is about 0.8.
By scattering one can verify the behavior anticipated above. One may
isotopically label a dilute fraction of the chains to measure the density profile Fig. 4.2
of individual chains in the solution. As anticipated, on scales smaller than ξφ Local volume fractions of monomer φ(r)
they have the fractal dimension of a self-avoiding chain (viz. 53 ), while over at distance r from an arbitrary monomer in
a semidilute solution of very low
distances much larger than ξφ they have the fractal dimension of a simple concentration φ. Solid curve: monomers
random walk (viz. 2) [1], as shown in Fig. 3.9. Each random-walk chain from the chain going through the origin.
has a self density of order R −1 . It is arbitrarily small as the chains are made Dot-dashed curve: volume fraction of
longer, with fixed total density. To achieve the given density with chains other chains. Scales are logarithmic, so
of increasing length there must be more and more chains in the volume R 3 that power laws appear as straight lines.
Within a blob size ξ the total volume
pervaded by one chain. The chains must interpenetrate. They interpenetrate fraction is dominated by the chain passing
more and more as φ increases above the overlap concentration φ ∗ . This through the origin. Beyond the distance ξ
interpenetration will become important when we come to consider how this the total volume fraction is dominated by
solution can flow. other chains. Simple dashed curve shows
excess density φ(r) 0 − φ(∞) 0 , which
decays exponentially to zero with a decay
4.2.2 Energy length ξs .
The blob picture is also useful for describing the interaction energy and the
osmotic pressure in the semidilute state. We saw above that interactions
between pieces of a chain larger than ξφ were unimportant relative to inter-
actions between different chains. Chemically distant parts of a given chain
might as well be different chains as far as the interactions are concerned.
This means that the interaction energy would be not much changed if the
chains were cut into pieces the size of a blob. But we can readily estimate
the energy of such cut-up chains, for these chains are at the overlap con-
centration φ ∗ . The distance between them is of the order of their size ξφ .
As we have noted at the beginning of this section, the interaction energy at
φ ∗ is of order T per cut-up chain piece—i.e., T per blob. The same is true
for the osmotic pressure:
T ξφ−3 . The interaction energy and pressure
are oblivious to the chain length in the semidilute regime.
According to this reasoning the quantity ξφ3 /T is predicted to be a
finite number of order unity regardless of the quantitative goodness of the
solvent, e.g., b2 /a 3 . In this sense it is like the quantity B2 /R 3 encountered
in the last section. And for similar reasons, this ratio too is universal. It has
been measured via scattering experiments for several polymers and several
solvents [3]. One finds /T = (3.2ξφ )−3 . Recalling that φ = φ(ξφ ) 0 =
(A/ξφ )−4/3 , one finds
Fig. 4.3
Reduced osmotic compressibility d /d ˜ φ̃
versus reduced volume fraction φ̃ = φ/φ ∗ Polyisoprene in cyclohexane
for various polymers and various good ~ ~ M = 640,000
d Π/d
solvents. For each polymer sample the ~
˜ φ̃ values are normalized so that they
d /d 5/4
approach 1 at low concentration. The φ̃
˜ φ̃ 2
values are normalized so that d 2 /d
approaches 2 at low concentration. Data
10
are shown for one sample of polyisoprene Polystyrene
in cyclohexane from [3] and three [Noda]
molecular weights of polystyrene from [4].
Theory [5] predicts that the pressure thus
normalized should be a universal function
of the normalized concentration. Theory
˜ φ̃ should scale as
also predicts that d /d
φ̃ p at high concentration, where the
exponent p is close to 54 . Thus
1
˜ φ̃)/φ̃ 5/4 should approach a (nearly)
(d /d 0.01 0.1 1 10 100
constant value at large φ̃. Both predictions ~
are well confirmed.
d 3r
S(q → 0) = lim d 3 r ρ(0)ρ(r ) exp(i q · r ) .
M q→0
We consider the scattering from some large volume within the sample. The
number of scatterers N in this volume is evidently d 3 rρ(r); its average is
Semidilute solutions 89
N̄ = ρ̄. In terms of N ,
S(q → 0) = 1/ρ̄ (ρ(r) − ρ̄)(N − N̄ ) = (N − N̄ )2 /N̄ .
given projected area. That is, the intensity seen within the beam goes as 1/ sin(θ),
where θ = 0 is the beam direction. You could find a way to illuminate the two
samples and view the two beams at different angles until their intensities looked
the same.
b3 p(t, r) giving the probability that the walker is at position r on the lat-
tice. If the walker is just above our imaginary plane, it has a probability 16
of passing down through the plane. The overall probability that this step
happens is 16 b3 p(t, r). On the other hand if the walker is just below the
plane, at r − b, the net probability that a particle moves up through the
plane is 16 b3 p(t, r − b). The probability current b2 jz moving through the
plane at this point is the difference of these two probabilities. Since p is
assumed to vary smoothly on the scale of a lattice step, we can express
b2 jz = − 16 (b3 /t)b∂p/∂z. Considering the current in arbitrary directions,
we have
† A comment
about the connection between
j = − 16 b2 /t ∇p.
(4.5) mean-squared displacement and ζ is in
order. Our random walker takes indepen-
dent steps of length b at every time interval
The coefficient in [ ] is called the diffusion constant ζ † . t. Thus the mean-squared displacement
To characterize the brownian motion further, we now ask what controls r 2 (t) = b2 t/(t). But a diffusing particle
the size of ζ . We can now consider the brownian motion of a simple particle, with diffusion constant ζ by definition has
a colloidal sphere of radius R. The solvent around such a sphere may be a mean-squared displacement in a given
direction of x 2 (t) = 2ζ t. We may also
treated as a continuum fluid with a velocity field v(r , t). The fluid velocity
recover the familiar diffusion equation for
on the boundary of the sphere must equal that of the sphere, but elsewhere the probability p(t, r), by noting that the
it may fluctuate because of thermal excitation: any configuration v(r ) has change of p at r during a time step is the sum
a probability given by the Boltzmann principle discussed in Chapter 2. of probability currents j onto site r during
Fortunately, one may determine the diffusion constant ζ without treat- that timestep: p = t∇ · j = tζ ∇ 2p.
This equation has the same form as the
ing these fluctuations explicitly, by means of a trick discovered by Albert
random-walk probability Eq. (3.12). See
Einstein [7]. He noted that there is a rigorous connection between ζ of Table 4.1.
any object and the drag coefficient describing the force required to pull
it along at some small speed v †† . To understand the Einstein relation, we †† This Einstein relation is a special case
imagine that our brownian particle is subjected to a weak gravitational of the fluctuation-dissipation theorem relat-
ing any small perturbation requiring work
field that exerts a force of magnitude F . As a result of this field we antic-
to a fluctuating quantity in the unperturbed
ipate that the particle gradually moves downward with an average speed system in thermal equilibrium [8].
v proportional to F ; this proportionality defines the drag coefficient :
F = v . Thus if the particle is at r with probability p(r), this motion
gives rise to a downward current density j (r) = p(r)v = −p(r)F / .
After the gravitational force has acted long enough for the colloidal parti-
cles to come to equilibrium, p(r) becomes a constant distribution given by
the Boltzmann Principle: p(r) = Z −1 exp(−F z/T ). Here F represents the
magnitude of the force, and z the height; the Z is a normalization constant.
Since this p(r) is not uniform, it gives rise to a diffusing upward current
jζ = −ζ ∇p, as discussed above. If gravity were suddenly turned off, one
could immediately observe this current: jζ = ζ F /Tp. But in the equi-
librium state with gravity acting, there must be no net current (otherwise
p(r) would not be constant in time): jζ + j = 0, or ζ F /Tp = pF / .
We conclude that (a) the induced speed v is indeed proportional to F as
anticipated, and (b) its magnitude is given by = T /ζ .
One may calculate the drag coefficient using conventional hydrody-
namics [9] and thus account for both drag and diffusion. But the basic
behavior of can be understood without hydrodynamic formalism. We
first recall the relationship between force and velocity v(r ) in a simple shear
flow, as discussed in Chapter 2. Here the flow was in the x direction but
varied only in the z direction. As we discussed, such flows transmit forces
92 Polymer solutions
The F̂ and r̂ are unit vectors in the F and r directions. The factor in ( )
depends only on the angle the vector r makes with v and with F . It is a
consequence of removing the longitudinal part of the v field. If we consider
only the downward velocity in response to a downward force, it is twice as
strong directly above and below the source as it is on the sides. There are
also velocities directed perpendicular to the force F . The fluid flows around
the particle. We shall need the Oseen tensor to analyze how the solvent
flows near tenuous objects like polymers. A mathematical explanation of
the Oseen tensor is given in Appendix A.
To find the drag force, we may equate the rate of work done by the force
with
3 the power P dissipated in the viscous flow. This latter has the form
d rηs γ̇ 2 . Here γ̇ is the shear rate introduced in Chapter 2; in a simple
shear flow, it is the derivative of the velocity. In more general flows it is a
combination of derivatives of the velocity. In the Oseen field γ̇
dv/dr
F /r 2 . Thus the power density
ẇ(r) = ηs γ̇ 2
F 2 /(ηs r 4 ).
The power dissipated beyond some radius R in the Oseen field is of order
∞
d 3 r ẇ(r)
F 2 /(ηs R).
R
Notes:
a ρ is mass density of the fluid.
m
composed of wave vectors
b The symbol v has to be interpreted using the Fourier transform method discussed in Appendix A. It is that part of v
⊥
transverse to v.
Fig. 4.4
Extra viscous dissipation from a sphere in
a liquid. Arrows show unperturbed
velocity field, with shear rate γ̇0 . Graph at
right shows dissipation rate along the
vertical dotted line through the sphere. Far
from the sphere the rate ẇ(z) has the
0
unperturbed value ẇ0 . Near the sphere the z
shear rate and dissipation rate are 2R
enhanced. Within the sphere, they vanish.
The
extra dissipation due to the sphere,
(ẇ(r) − ẇ0 ) d 3 r is suggested by the ·
shading.
w(z) w· 0
Motion in a polymer solution 95
solute present, the shear has a uniform value γ̇0 everywhere and a uniform
dissipation rate ẇ0 = ηs γ̇02 . A solid object such as a sphere perturbs this
flow and increases the dissipation. With such objects present, greater dissip-
ation occurs for a given macroscopic flow rate: the macroscopic viscosity is
increased. This increase is a major way of sensing the presence of structures
in a fluid. To understand how this perturbation works, we consider a single † In a viscous liquid, and many other sys-
colloidal sphere of radius R in a shear flow between two horizontal sliding tems where the forces are proportional to
plates. The solvent is obliged to flow around the sphere. This increases the the velocities, the energy dissipation rate
shear rate γ̇ outside the sphere. The sphere alters the shear flow from its ini- must be a quadratic function of the velocit-
tial uniform state of minimum dissipation; thus it increases the dissipation† . ies. For such systems the steady-state velo-
The dissipation per unit volume ẇ(r) is zero inside the sphere. But ẇ(r) city field is that which minimizes this energy
dissipation rate. This can be proven from the
must be increased just outside the sphere. For example, the integral of the equations of hydrodynamics, which amount
shear rate γ̇ along the vertical line through the center is the difference in to Newton’s laws of motion applied to these
speed of the two plates, as for any other vertical line. Since the shear rate linear forces ([10], Section 4.6).
vanishes within the sphere, it must be greater than γ0 above and below in
order to compensate. The perturbation in velocity v(r) caused by the sphere
resides mostly near the sphere. It can be thought of as the sum of many
Oseen tensors centered at each point of the sphere. Since these produce a
velocity field that falls off as 1/r, the resulting shear rate γ̇ (r) = dv/dr
falls off as 1/r 2 as discussed above in connection with the drag force. In
order to make up the deficit in shear rate within the sphere, γ̇ /γ̇0 must be
significantly greater than 1 for distances r
R. In this same region the
dissipation ẇ(r)/ẇ0 is also significantly greater than 1. Thus the effect of
the sphere is to increase the dissipation by a finite factor
in a region of the
order of the sphere size. That is, the extra dissipation d 3 r(ẇ(r) − ẇ0 ) is
roughly ẇ0 R 3 .
By accounting for the excess dissipation, we may find the increased
viscosity. On the one hand, the viscosity is defined by ẇ = ηγ̇ 2 , where γ̇
is the average velocity gradient. On the other hand if the fluid has volume ,
3 ẇ(r) 3
ẇ = ẇ0 + (ẇ(r) − ẇ0 ) d r = ẇ0 + −1 d r ,
ẇ0
Here the first term is the dissipation due to the unperturbed flow throughout
the fluid. The second term is the excess dissipation near the sphere. The
integrand is significant over a volume comparable to the sphere volume
V = 43 π R 3 . The integral can be found explicitly in terms of the known
velocity field around the sphere [9,11]. The result is (ẇ(r)/ẇ0 −1) d 3 r =
5
2 V . The dissipation from a sphere is the same as if the dissipation were
doubled over a volume equal to 52 that of the sphere. Comparing the two
expressions for the dissipation, we infer
5V
η = ηs 1 + .
2
If the solution contains N spheres far from each other, the dissipation is the
sum of that due to each, so that
2
5 NV N
η = ηs 1 + +O . (4.8)
2
96 Polymer solutions
F
v∞ − vi (r)
v∞ − d 3 r0 ρ(r0 ) .
ηs (r − r0 )
i
Motion in a polymer solution 97
Table 4.2 Universal ratios of length dimensions of a dilute polymer solutiona Notes:
a Adapted from [14], Table VI, using their
Ratiosb RG /Rh /R
RG Rt /Rh Rv /Rh
h original data and earlier data they compiled
from various sources.
Hard sphere (3/5)1/2 = 0.775 1 1 1 b R is the hydrodynamic radius—the radius
h
Ideal chain (theory)c 1.48 1.91 0 1.23 of a hard sphere with the same diffusion
Self-avoiding chain (theory)c 1.56 2.01 1.04 1.14 constant as the polymer. RG is the radius
Polyisoprene in cyclohexaned 1.39 1.79 0.95 1.11 of gyration, as defined in Problem 3.11.
Polystyrene in benzened 1.51 1.94 1.01 1.03 RG ≡ ( 3 )−1/2 R is the radius of a uniform
5 G
Poly(α-methylstyrene)d 1.55 2.0 1.05 1.11 sphere with the same RG as the polymer.
Rt is the thermodynamic radius, that of hard
spheres having the same excluded volume
B2 as the polymers. Rv is the viscometric
If we take as our origin some bead of the polymer, then on average ρ(r0 ) = radius, that of hard spheres having the same
ρ(r0 ) 0 ∼ r0 D−3 , where D is the fractal dimension of the polymer. Thus intrinsic viscosity as the polymers.
c These theories estimate the given ratios by
R treating ≡ (d − 4) as a small parameter
v(r) − v∞ ∼ d 3 r0 r0 D−3 (r − r0 )−1 . and calculating the given ratio to first or to
second order in .
d These experiments used scattering and
The integration is now of the same form as Eq. (3.26) for MAB , the viscometry to measure the various R values
for different molecular weights. The lim-
number of intersections of two fractals with dimensions D and 2 in three- iting ratios for high molecular weights are
dimensional space. For all fractals with D > 1, including our polymer, the reported.
integral goes as a positive power of R. We say that such fractals are opaque
to the flow. We note that opacity is attained for large enough R no matter
how small the individual beads were or how weakly they interacted with
the fluid.
The result of this flow opacity is to invalidate our supposition that each
bead feels the unperturbed speed v∞ . For this leads the conclusion that
v(r) − v∞ diverges with R. In that case the backflow would be stronger
than the asymptotic speed v∞ . To avoid this absurd result, we conclude that
the fluid speed within the fractal is actually much smaller than v∞ , so that
the forces Fi on each bead are also smaller than an isolated bead would feel.
Thus, the fluid cannot move transparently through the fractal; the flow must
go around it. The flow speed throughout the fractal volume is suppressed
by at least a finite factor.
Because a fractal screens out the flow within itself qualitatively, it also
alters the flow far away. The far-field velocity must extrapolate to a finite
fraction of v∞ at r = R. This means that the far field is that of a hard sphere
whose radius is of order R † . Thus a fractal with D > 1 in a uniform flow †This hard-sphere behavior was first pro-
causes dissipation like that of a hard sphere whose radius is of order R. posed by Zimm, and is called the Zimm
model [13].
Accordingly, its drag coefficient or hydrodynamic radius Rh are equal to
R up to a finite factor independent of the microscopic structure. A polymer
must also perturb a shear flow like a hard sphere of comparable radius.
Thus its effect on the viscosity is like that of a hard sphere whose radius Rv
is of order R. Like other asymptotically finite ratios we have encountered
above, the Rh , Rv as well as the thermodynamic Rt defined above are uni-
versal multiples of the radius of gyration RG for any self-avoiding polymer.
Table 4.2 summarizes the experimental data on these ratios. Remarkably
these universal ratios are not far from 1.
4.4. Diffusive analog of hydrodynamic opacity We saw in this section that tenuous
fractals can interact strongly with each other by causing the surrounding solvent
to move. We suggested that this interaction could be understood in terms of
the mutual opacity of two fractals. This analogy is clearest if we consider the
98 Polymer solutions
expressed as −b3 ρ2 v(r), where ρ2 is the density of spheres. The other
force on this volume element is the viscous force ηs b2 dv/dz on the bottom
and the top of the element. In steady state, if the fluid is not accelerating,
these forces must balance: ηd 2 v/dz2 = ρ2 v. This familiar equation has
solutions exponentially growing with z or exponentially decreasing. The
physical solution is the decreasing one† : v(z) = v(0) exp(−z/ξh ), where † Strictly, some of the exponentially grow-
ξh−2 = ρ2 /ηs . The flow is now confined to a region of size ξh near the ing solution is needed to satisfy the bound-
ary condition v(h) = 0. Here we restrict
bottom plate. Most of the momentum is absorbed before reaching the top ourselves to a thick layer with h ξh , so
plate. The same hydrodynamic screening length ξh appears when many that this part is negligible.
such one-dimensional flows are combined to make more general flows. For
example the Oseen tensor describing flow around a moving sphere acquires
a damping factor exp(−r/ξh ). Thus if a small sphere is moved through a
fluid containing stationary spheres, the flow around it is not much perturbed
for distances r ξh . But for distances r
ξh , the velocity is reduced by
a substantial factor.
If a small sphere is moved through a semidilute polymer solution, hydro- †† You may wonder what the polymers do
dynamic screening also occurs, since the polymers absorb momentum being with the momentum they have absorbed.
The polymers are not attached to the con-
supplied by the moving sphere†† . If the small sphere is inserted at random
tainer, so they cannot absorb momentum
into the fluid, the typical distance to a monomer is the correlation length without starting to move themselves. We
ξφ introduced above. Thus if one examines the flow at distances r ξφ , will consider the fate of the momentum
there are essentially no monomers there to impede it; accordingly, the flow absorbed by the polymers in Problem 4.8.
is unscreened. However, at distances r
ξφ the polymers alter the flow.
If our sphere were near an isolated polymer of radius ξφ , the flow would
be damped appreciably within the coil volume because of its opacity, as
discussed above. The same is true for a solution of polymers of size ξφ at
their overlap concentration φ ∗ . It is still true if these polymers are connected
together to form a semidilute solution with the same φ and ξφ . At distances
of order ξφ from the sphere, there is an appreciable chance that monomers
are present to impede the flow speed by a finite factor. Since screening
reduces the flow by a finite factor at r
ξφ , we infer that the screening
length is of order ξφ . At distances much larger than ξφ , the polymer solution
may be regarded as a more-or-less uniform mass of momentum-absorbing
blobs. Thus exponential damping of the flow is expected, as in the last
paragraph, with screening length ξh
ξφ .
From this fact we can infer the behavior of ζc far above the concen-
tration φ ∗ . We may attain this semidilute regime starting from a solution
near φ ∗ by joining chains together, as we have done previously. Then each
of the original chains becomes a blob of the solution, and each original
chain’s size R is now of the order of the blob radius ξφ . The joining of these
polymers does not change the distribution of monomers in space very much;
accordingly the solvent is free to flow around and through the blobs about as
much as it could with the original solution before joining. Thus a gradient of
density produces a current of monomers of about the same size whether the
polymers are joined or unjoined. This means that the cooperative diffusion
constant ζc is about the same in the two cases:
plex amplitudes ρ̃q . Applying the diffusion law of Eq. (4.12) to ρ(r) we
find an equation for the time derivative of ρ̃(q): ∂/∂t ρ̃(q) = −q 2 ζc ρ̃(q).
The amplitude evidently decays exponentially with a decay time τ (q) =
(q 2 ζc )−1 . As shown in the last chapter, waves scattered at wavevector q
have an intensity I (q) proportional to the square of this amplitude; thus,
Motion in a polymer solution 101
Table 4.3 Universal semidilute length ratios inferred from polystyrene solutionsa Notes:
a These ratios were inferred by Huang
Quantity Theta solventb Carbon disulfidec Toluened and Witten [18] using scattering and
cooperative-diffusion data from the litera-
ture.
b These values were inferred from data on
A (nm)e 0.32 (± 7%) 0.28 (± 7%) 0.32(± 7%)
ξφ /ξs f 0.61(± 9%) 1.23 (± 8%) 1.23 (± 9%) polystyrene in the theta solvent cyclohex-
ξ /ξs g 2.97 (± 5%) – 3.81 (± 6%) ane at 36 ◦ C [19–22].
c These values were inferred from data on
ξf /ξs h – 1.27(± 1%) –
perdeuterated polystyrene in the good sol-
ξζ /ξs i 4.29 (± 10%) – 1.65 (± 7%) vent carbon disulfide at 20 ◦ C [23–25].
r2 j 3 (± 22%) 1.2 (± 14%) 1.1 (± 23%) d These values were inferred from data on
r4 48 (± 28%) 4.5 (± 21%) 4.3 (± 29%) ordinary polystyrene in the good solvent
toluene at 25 ◦ C [26–28].
e Fractal amplitude defined by φ(r) =
0
(A/r)D−d , as in the subsection “Semidilute
I (q, t) ∝ |ρ(q, t)|2 = I (q, 0)e−2t/τ (q) . The decay time is roughly the time Solutions” above.
f ξ is the correlation length inferred from
for the density to diffuse a wavelength 2π/q. s
scattering: S(q) = const.(1 − (qξs )2 +
If many initial states were prepared at random with proper Boltzmann- O(q 4 )). ξφ is the distance from a monomer
weighted probabilities, and the intensity I (q, t) monitored for each, each at which the local volume fraction φ(r) 0
experiment would show a different initial intensity, but the same decay of a large single chain inferred from high-
q scattering is equal to the overall volume
time τ (q)/2. One may measure this decay conveniently by measuring the fraction φ: φ(r) 0 |r=ξφ = φ, as described
correlation function I (q, t + t )I (q, t ) t , where t is the average over in the subsection “Semidilute Solutions.”
initial times t . Under the diffusive law, this correlation function has the g ξ
is the length scale implicit in the
−3
form [16] osmotic pressure : = T ξ .
h ξ is defined by S(q) = S(0)(qξ )−D for
f f
I (q, t + t )I (q, t ) t = (I (q) − I (q) )2 e−2t/τ (q) + I (q) 2 . (4.15) q in the fractal regime.
i ξ is the radius of a sphere whose stokes
ζ
diffusion constant is the same as the coop-
The scattering from an equilibrium fluid behaves in this way, as well. The erative diffusion constant. It is rigorously
random thermal fluctuations have the effect of preparing initial states of related to ξ and the permeability length
ξp , as discussed in the text.
random amplitudes, which then decay according to the diffusion equation. j r is the reduced moment of the con-
n
At any given moment the intensity is the cumulative effect of many partially centration profile defined by rn ≡
∞
decayed random fluctuations. 0 dxx (φ(xξφ ) 0 /φ − 1). The r2 and r4
n
Careful scattering measurements have been done to test the predicted determine the ratios ξs : ξφ : ξ : ξf , as
explained in [18].
behavior of τ (q). The data [17] show good consistency with the φ depen-
dence predicted above. They lead to the estimate ζc = T/[6π ηs (αξφ )],
where α = 1.3 ± 10%, as inferred from Table 4.3.
4.5. Cooperative diffusion and permeability The cooperative diffusion constant
ζc tells the flow of concentration in response to a concentration gradient. One
can readily infer from this the flow of fluid in response to an osmotic pres-
sure gradient. The ratio of fluid velocity to a pressure gradient is defined as
the permeability P : v = j/ρ = −(P /η)∇p. The permeability is defined for
flow of any fluid; no polymers or other solute need be involved. Here η is the
viscosity of the fluid and p is the pressure, j is the current density of fluid and
ρ is the particle density of the fluid. Evidently the permeability has the dimen-
sions of a length squared. (a) Find the permeability P relating the average
velocity to the pressure gradient in a long circular pipe of radius R. Recall that
the velocity profile in the pipe is parabolic while the pressure is constant over
the cross-section of the pipe. You can find the pressure gradient by equating
the power required to maintain the pressure with the viscous energy dissipation
in the fluid. (b) If monomers with number density ρn are flowing to the right
with a current density jn , what is the average speed of monomers with respect
to solvent? Assume monomers have the same volume as solvent molecules.
Considering the monomers as fixed, what is the solvent current js in terms
of jn and volume fraction φ = ρn /(ρn + ρs )? (c) A semidilute solution has
a monomer density ρn and solvent viscosity ηs . Relate its permeability P to
its cooperative diffusion constant ζc and the osmotic compressibility d/dρn .
102 Polymer solutions
k
= ri (t)rj (t) .
i,j
see whether we can account better for its damping. If a block of Jello were
stretched to double its natural length, it would eventually expel a finite fraction
of the water inside. The stretched water wants to conserve its volume when
stretched, while the polymer network wants to increase its volume. Thus when
the Jello is distorted, there is a current of water through the polymer network.
From this we can estimate the dissipation of Jello. From the oscillation fre-
quency and density of Jello, we know its modulus G is roughly 4000 J/M3 (or
40,000 erg/cm3 ). This modulus is of the same order as the osmotic pressure.
(a) From this information, estimate the blob size ξ . A block of Jello is rapidly
subjected to a shear strain of magnitude γ . This strain produces an osmotic
pressure difference of order Gγ between the center and the outside. (b) Esti-
mate the rate of energy dissipation per unit volume caused by the resulting
flow in a cube of side b. (c) Assuming this dissipation has a small effect on
the motion, estimate the damping rate for a jiggling cube of Jello 5 cm on a
side.
simple model of motion is called the Rouse model [5]. The Rouse model
describes self-diffusion reasonably well when the concentration is not too
far above φ ∗ . But for solutions far above φ ∗ we must consider another
difference between our original solution of subchains and the final solution
of connected k-chains: entanglement. The forbidding task of accounting for
entanglements quantitatively was reduced to a simple and intuitive model
by Edwards and by deGennes [5, 10] in the 1970s. Our discussion below
recounts these ideas.
Intuitively, entanglements between chains impose constraints that go
beyond mere hydrodynamic drag. We may begin to anticipate how they
influence a given chain by replacing the other chains with fixed, regular
obstacles, such as a regular lattice of bars of spacing ξ . If a random walk
polymer threads through such a “jungle-gym” lattice, it certainly cannot
diffuse according to the Rouse model discussed above. If a random force
pushes the chain against one of the bars, the bar exerts a restoring force
that prevents the chain from crossing it. The chain is confined. If each
end of the chain were to be fixed, this confinement would be permanent;
the chain must forever follow the same path through the jungle gym that
it followed initially. But if the ends are free to move, there is one type
of motion it can perform without a restoring force from the bars. This
is for each monomer to move in the direction of its successor along the
chain. The result of such a motion is equivalent to removing a bit of the
chain at one end and adding it on to the other end. This motion is called
reptation [5].
Of all the different small random motions of our chain only a reptational
motion results in no restoring force from the bars. Thus for this motion
the random thermal forces are free to act unopposed. This means that the
chain is free to do random small steps along its own path: it performs
curvilinear brownian motion. The displacement along the path s(t) obeys
s 2 (t) ∝ t. The constant of proportionality is the curvilinear or tube dif-
fusion constant ζt . This diffusion motion is the same as it would be if the
contorted sequence of blobs were straightened out into a straight path of
length L
ξ k. Each section of the chain of size ξ is subjected to independ-
ent random forces, as in the Rouse model above. Accordingly the diffusion
constant for the overall motion of a chain containing k such sections is
reduced by a factor of k. The curvilinear diffusion constant is simply the
Rouse diffusion constant derived above, viz. ζc (M(ξ )/M). Under this dif-
fusion, the chain ultimately moves a curvilinear distance equal to its path
length L. At this point the chain is no longer constrained by its original
entanglements. The time τrep required is given by L2 /τrep
ζc (M(ξ )/M).
This τrep plays a central role in many phenomena of entangled polymers.
Using L = ξ(M/M(ξ )), we may infer the scaling of τrep :
pressure . The polymers also store energy when the solution is distorted,
even with no change of volume. We discussed how this energy is measured
in the Experimental Probes section of Chapter 2. There we told how rheo-
meters apply a known, oscillating shear strain and measure the stress that
oscillates in phase or in quadrature with that strain.
In strongly interpenetrating polymer solutions, this stress occurs because
of entanglement. To understand the effect, we return to our jungle-gym
model. We impose a step shear strain like that shown in Fig. 2.3. This
amounts to tilting the jungle-gym lattice so that each square of the lattice
distorts (affinely) into a parallelogram, like the sample as a whole. This
distorts the shape of the polymer and its internal density profile. If we con-
sider the end-to-end vectors of each blob of the chain, they are initially
isotropically distributed. But after the distortion, the vectors become aniso-
tropic. They are more concentrated along the long diagonal of the sheared
square and less concentrated along the short diagonal. The anisotropy has
a strength of order unity under one unit of shear. Each of these blob distor-
tions costs energy. To estimate this energy, we imagine that the anisotropic
distribution of vectors were produced by an external potential U (θ ). The
distribution of angles f (θ) = exp(−U (θ)/T ). In order to produce an f (θ )
that varies over a factor of order unity, U (θ ) must be of order T . The work
required to turn on this potential must also be of order T , since the energy
U affects most of the possible angles θ. This means that the work required
to make the shear distortion is of order T for every cell of the lattice, or for
every blob of the chain.
4.C Suggested experiment: Rubber elasticity and temperature The elastic
restoring force of rubber is essentially the restoring force of an elongated random
walk. This means the force should be proportional to the absolute temperature,
provided the basic local structure of the rubber polymers does not depend on tem-
perature. If this is true, one should observe a roughly 30% increase in the spring
constant as the temperature is changed from 0 ◦ C (273 K) to 100 ◦ C (373 K). Devise
an experiment with rubber bands to test this prediction. You might use weights, a
ruler, boiling water, and ice water.
Since a unit of shear strain stores energy at a density T /ξ 3 , the solution
has a step-strain modulus G0 of this order, owing to entanglements. This
G0 is independent of molecular weight and is of the order of the osmotic
† The ratio of G0 to is predicted to be pressure † . The solution retains this stored energy until the entangle-
universal in the semidilute regime in good ment constraints are released. If each chain is prevented from reptating, by
solvents. But this universality has not been
immobilizing its ends or by crosslinking, the stored energy lasts indefinitely
verified experimentally. Instead, appears
to vary more strongly with solvent quality and the sample is an elastic solid—a gel. Similarly a polymer melt immob-
than G0 does. As concentration increases ilized by a few crosslinks makes an elastic solid—a rubber. The order of
above the semidilute regime or solvent qual- magnitude of this modulus has been mentioned in Problems 2.6 and 2.7.
ity diminishes, the predicted universality is The highest-modulus rubbers are those with the highest density of entan-
expected to be compromised. These limita-
glements. Empirically the largest moduli for a weakly crosslinked rubber
tions may account for why the universality
has not been observed. are of the order of 1 atm, or 105 J/m3 .
4.8. Stokes diffusion in semidilute solution In treating the self-diffusion of a poly-
mer in semidilute solution, we have pretended that the neighboring chains
completely remove any momentum in the fluid, so that the flow around a point
source of momentum is completely screened out. But this can not be completely
true, since a chain that absorbs the momentum can only give it to the adjacent
solution; the chains do not remove the momentum entirely. This means that if
Motion in a polymer solution 107
Fig. 4.5
How reptation leads to stress relaxation.
Top picture shows a chain confined in a
jungle-gym representing the other chains.
Bottom left picture shows the effect of a
horizontal stretching: the whole structure
is distorted and each section of the chain
becomes anisotropic. Bottom right picture
shows the chain after some time has
passed. The heavy section on the left
shows the part of the chain that has
reptated out of the initial constraints. This
part is isotropic and holds no stress. The
dashed section on the right shows the part
of the initial chain that has reptated away.
This section is equal in length to the heavy
section on the left.
If the chains are free to reptate, the stored energy can relax (Fig. 4.5). In
our jungle-gym model, a given cell of the lattice continues to store energy
until the chain end has reached that cell and relaxed the constraint. The
same is true for a real entangled solution. The time required to relax a finite
fraction of these constraints is evidently the reptation time τrep . By extending
this reasoning, one can make detailed predictions of the fraction of the
initial stress relaxed as a function of time. These predictions agree well
with experiment [10] (Fig. 4.6).
4.9. Partial stress relaxation function in polymer melt If the volume fraction of
polymer is increased above the semidilute regime, the correlation length shrinks
to the size of a monomer. This is the concentrated solution regime. When all
the solvent is removed, the polymer liquid is called a melt. In this regime most
of what we derived about the semidilute regime remains valid. The point of
this problem is to calculate how the stress G(t)γ begins to fall below its initial
value G(0)γ in a strongly entangled solution or melt. The stress falls over time
because part of each chain has reptated out of its initial (distorted) tube and into
a new (isotropic and stress-free) set of entanglements, as shown in Fig. 4.5.
In effect, part of the initial distorted tube has disappeared. The lost stress is
just the initial stress times the fraction of initial tube which has disappeared.
This disappearance occurs as the chain executes a random walk forward and
backward in its tube. At a given time t the amount of tube removed from the
left is the maximum curvilinear distance the chain has traveled to the right in
that time. This maximum is roughly the root mean square curvilinear distance
108 Polymer solutions
Fig. 4.6
Storage modulus G versus oscillation PS 160 °C
8
frequency ω for narrow-distribution
polystyrene melts, reprinted with
permission from [30, Fig. 2]. © 1987 7
American Chemical Society. Molecular
weight ranges from 8.9×103 (curve L9) to 6
5.8×105 (curve L18). The modulus drops
log G⬘ dynes/cm2
off below a characteristic frequency on
each curve. This frequency is the inverse 5
of the stress relaxation time. The high
molecular weight samples have a large 4
range of frequencies where the modulus is
independent of frequency. Over this
plateau region, the liquid behaves like an 3 L34
L27 L37
elastic solid. The associated modulus is the L18 L15
step-strain modulus discussed in the text. L14
2 L19
(To achieve the huge reported range of L5 L22
oscillation frequencies, the experimenters L16 L9
L12
use a trick. They cool the short polymers to 1
slow their motion and heat the long ones to –6 –5 –4 –3 –2 –1 0 1 2 3 4
speed up their motion. The temperature log ωaT s–1
effect can be separately quantified. The
reported frequencies have been multiplied
by a factor aT , determined separately for
each sample, to show what the moduli
would be at a fixed temperature of 160 ◦ C.)
traveled. In a given time, about the same amount is lost on each end. (a) Derive
a function for the fractional stress lost as a function of time t after the initial
strain in terms of the curvilinear diffusion constant ζt . Do not worry about late
times when the amount of tube lost is nearly the entire tube. And do not worry
about numerical prefactors.
Knowing the step-strain modulus and the relaxation time, we can determ-
ine the viscosity η of the entangled solution. As we saw in Chapter 2, the
viscosity scales as the modulus times the relaxation time. If one increases
molecular weight M at a given concentration φ, the reptation time increases
as M 3 , while the modulus remains constant. Thus the viscosity should
increase as M 3 . Experimentally, the viscosity agrees better with M 3.4 than
with M 3 . After many years of puzzlement it appears that this disagreement
is not fundamental. It appears that in the case of viscosity it is particu-
larly difficult to attain the asymptotic M M(ξ ) behavior. Several effects
alter the viscosity when M is not asymptotically large. One may estimate
these effects quantitatively and they seem to explain why an apparent M 3.4
behavior is observed in the experimental range of M [10].
If a liquid undergoes step strains that are more frequent than the relaxa-
tion rate, the stress from each step progressively builds up. Thus, if a steady
shear rate exceeds the inverse relaxation time, substantial elastic stress
can build up. In particular, normal stress builds up as the polymer chains
stretch along the streamlines of the flow. This stress from stretched poly-
mer chains is the origin of the unusual flow behavior pictured in Fig. 1.6.
The appearance of steady-state deformation and its attendant change in
the stress at shear rates exceeding the relaxation rate is known as the
Cox–Merz rule [31]. Cox–Merz behavior is observed in many structured
Appendix A: Origin of the Oseen tensor 109
fluids; this is natural since the large structures found in such fluids make for
long relaxation times. Thus the Cox–Merz criterion for shear-rate dependent
viscosity can be met at moderate shear rates.
4.4 Conclusion
In this chapter we have recounted the basic behavior of large, flexible
chain molecules in solution. The qualitative elements of randomness,
self-repulsion, and entanglement combine to make distinctive spatial struc-
ture, energy storage, and dynamic response. We have seen that many of
these properties can be understood using primitive geometric notions such
as the number of contacts between two fractals. These geometric properties
allow one to predict all the asymptotic scaling behavior of long polymers
as a function of molecular weight and concentration. In many cases the
common geometric basis of these polymer properties allows one to predict
one property like the diffusion constant in terms of a completely different
property like the osmotic pressure, as exemplified in e.g., the universal
ratios of Table 4.2. In the coming chapters we explore to what degree such
generality can be extended to other forms of structured fluid.
Thus
ηṽα + ikα p̃/k 2 = f˜α /k 2 .
110 Polymer solutions
0 + i p̃ = kα f˜α /k 2 .
Thus
1
(fα − kα (kβ f˜β /k 2 )).
ηṽα =
k2
For a localized force F at the origin f̃k = F /(4π ). The desired Oseen
tensor is thus given by
|F | 1
v(r) = 2
F̂ − k̂(k̂ · F̂ ) eik·r .
4π η k
k
The first term must integrate to F /(ηr) up to a numerical factor. The second
term must be in the radial r̂ direction and must be such that ∂β vβ = 0. The
result is the Oseen tensor given in the text.
so that
dρn P 1 dρn 1 dφ
ζc = ρn ; P = ηs ζc = ηs ζc .
d ηs ρn d φ d
References 111
Thus the solvent flows through the polymer network as though it were a bunch
of pipes of radius R
ξ .
(e) Specifically, R 2 /8 = (3 − D/18π )ξ 2 or using D
53 , R
0.4ξ . We can get
a more refined estimate using the data in Table 4.3: ξ = 3.81ξs ; ξζ = 1.65ξs .
Then
R
(2.5 ± 10%)ξs ,
where ξs is the scattering correlation length defined by S(q) = S(0)
(1 − (qξs )2 + · · · ).
References
1. G. Jannink and J. Des Cloizeaux, Polymers in Solution (Oxford, UK: Oxford
University Press, 1992).
2. N. Nemoto, Y. Makita, Y. Tsunashima, and M. Kurata, Macromolecules 17 425
(1984).
3. M. Adam, L. J. Fetters, W. W. Graessley, and T. A. Witten, Macromolecules 24
2434 (1991).
4. I. Noda, N. Kato, T. Kitano, and M. Nagasawa, Macromolecules 14 668 (1981).
5. P. G. deGennes, Scaling Concepts in Polymer Physics (Ithaca, NY: Cornell,
1979).
6. M. Daoud and G. Jannink, J. Physique 37 973 (1976).
7. A. Einstein, Ann. Physik 17 549 (1905); 19 371 (1906).
8. R. K. Pathria, Statistical Mechanics (Oxford: Pergamon Press, 1972), chap. 13.
9. J. Happel and H. Brenner, Low Reynolds Number Hydrodynamics With Spe-
cial Applications to Particulate Media, 2nd rev. ed. (Leiden: Noordhoff
International Publishing, 1973).
10. M. Doi and S. F. Edwards, The Theory of Polymer Dynamics (Oxford: Oxford
University Press, 1986).
11. A. Einstein, Ann. Phys. 34 591 (1911).
12. P. E. Rouse, J. Chem. Phys. 21 1273 (1953).
13. B. Zimm, J. Chem. Phys. 24 269 (1956).
14. N. S. Davidson, L. J. Fetters, W.G. Funk, N. Hadjichristidis, and
W. W. Graessley, Macromolecules 20 2614 (1987).
15. B. Ewen and D. Richter, Adv. Polymer. Sci. 134 1–129 (1997).
16. See e.g. B. J. Berne and R. Pecora, Dynamic Light Scattering: With Applications
to Chemistry, Biology, and Physics (New York: Wiley, 1976).
17. A. Z. Akcasu, G. C. Summerfield, C. C. Han, C. Y. Kim, and H. Yu, J. of Polym.
Sci. Part B Polym. Phys. 18 863 (1980).
18. J. R. Huang and T. A. Witten, Macromolecules 35 10225 (2002).
19. M. Adam and M. Delsanti, Macromolecules 18 1760 (1985).
20. P. Stepanek, R. Perzynski, M. Delsanti, and M. Adam, Macromolecules 17 2340
(1984).
21. J. P. Cotton, M. Nierlich, F. Boué, M. Daoud, B. Farnoux, G. Jannink,
R. Duplessix, and C. Picot, J. Chem. Phys. 65 1101 (1976).
22. J. Roots and B. Nyström, Macromolecules 13 1595 (1980).
112 Polymer solutions
Colloidal dispersions† , are homogeneous suspensions of solid particles in a †Historically colloids referred to any sticky,
fluid. Most inks and paints are colloidal dispersions. The particle size ranges gelatinous material. Here, we follow the
more recent usage.
from nanometers to microns. Such a particle is large enough to approach
its bulk properties on the atomic scale but small enough that its thermal
energy dominates its gravitational energy†† . As discussed in Chapter 1, †† This means that gravity is not strong
colloidal suspensions or dispersions have a variety of direct applications as enough to drive the particles to the bottom or
top of the container, as discussed in Problem
well as serving as intermediates in various processing technologies, such as
2.2.
ceramics. Colloidal particles tend to attract each other and thus aggregate
together. If this aggregation grows unchecked the aggregates eventually
migrate to the top or bottom of the container, thus destroying the desired
dispersed state. Much of colloid science is concerned with maintaining the
dispersed state. This is the problem of colloidal stabilization. Stabilization
is achieved by modifying the particles’ surfaces to prevent aggregation.
Thus understanding and control of interfacial forces is central to colloid
science and to this chapter.
In this respect, colloidal dispersions resemble emulsions and foams. As
noted in Chapter 1, emulsions differ from colloids in that the dispersed
objects are liquid droplets rather than solid particles. In a foam the dispersed
objects are in the gas phase. Most of our reasoning about colloidal stability
applies equally to emulsions and foams.
In the next section we discuss the main causes of the attractive interac-
tions that compromise the stability of colloidal dispersions. The following
section deals with means of combating these attractions via repulsive forces.
The last section discusses consequences of strong colloidal interaction: the
particles take on various forms of spatial organization: crystals, liquid crys-
tals, fractal aggregates, and magnetic states. Appendix A treats one of the
major mechanisms of attraction in detail. Appendix B explains the various
forms of fractal aggregation mentioned in Chapter 1.
How are colloids made? There are two general approaches to the pro-
duction of colloidal-sized particles: dispersion††† and condensation. These ††† “Dispersion” is used to denote both a
are reviewed in Everett’s [2] book on colloid science. Dispersion means the process for making colloids and the result-
ing state in which particles are dispersed in
breakup of larger-sized particles into smaller ones, usually by some type of
a liquid.
rough mechanical treatment, such as grinding or shaking. It is difficult by
such means to apply sufficiently strong stresses to achieve the lower range
of colloidal dimensions. More importantly, dispersion generally results in
* Thischapter is heavily based on an earlier draft by Phillip A. Pincus, itself based on [1].
The organization of that draft, some text, and several figures survive in the present version.
114 Colloids
5.1. Fluctuating dipoles We assumed above that the instantaneous dipole was
directed parallel to r. Recalculate U1,2 , relaxing this assumption.
This result merits some discussion. The stability of the unpolarized atom
implies α ≥ 0. Then the interaction energy U1,2 is negative definite; the
interaction is attractive independently of whether the fluctuating moment
1 is pointing toward or away from atom 2. We can now average U1,2 over
μ
all fluctuations of atom 1; in Eq. (5.3), the factor μ21 is then replaced by
its average μ21 , which is nonzero even though μ1 = 0. The r −6 force
law is characteristic of the unretarded van der Waals interaction [4]. Let us
try to estimate the order of magnitude of U1,2 . The mean square value of
the fluctuating dipole moment, μ21 , is related to the elementary atomic
dipole given by the electronic charge e multiplied by an atomic length, a,
i.e., μ21 ∼ (ea)2 ; the atomic polarizability scales with the atomic volume,
α ∼ a 3 [3]. The interatomic coupling may be generally expressed as
For simplicity we have discussed the van der Waals interaction in the
context of single atoms and atomic media. But this mechanism of attraction
is much more general. It arises simply because materials are polarizable
and the charge within them fluctuates. Since all condensed matter has
these properties, van der Waals attractions operate in all condensed matter,
atomic or molecular, insulating or metallic. Because the van der Waals
energy has the same distance dependence for all materials, all like surfaces
separated by another medium have characteristic Hamaker constants, and
attractive interactions of the form of Eq. (5.6). But these attractions are
just a special case of an even more general attractive mechanism, to be
discussed next.
5.3. Derjaguin approximation The Derjaguin approximation relates the force
between curved bodies to the interactions between flat surfaces composed of
the same materials. This approximation can be used to find the force between a
sphere and a plane, equivalent to the interaction between two spheres. Consider
a sphere of radius b at a distance of closest approach h from a half space of
the same material, with b h. Suppose that the interaction energy per unit
area between to half spaces separated by a distance z is known to be v(z). (a)
Show that if v(z) decays sufficiently rapidly with z, then the force between the
sphere and the plane is F (h) = 2π bv(h). (b) Using the results of Problem 5.2,
show that the van der Waals interaction potential energy between a sphere and
a plane is of the form Vb (h) = −kH (b/h). What is the constant k?
When we insert the particle, ψ(0) adjusts itself to minimize W1 . One readily
verifies that ψ(0) = −W /χ and W1 = W0 − 12 W 2 /χ.
Our first particle has perturbed the quantity ψ in the whole vicinity of
the origin. The perturbed ψ at x can be expressed f(x)ψ(0). The simplest
and most common behavior is for f(x) to decay monotonically to zero
as x increases, as in our foam-rubber example above. We shall assume
that f(x) behaves in this way. Later we shall identify conditions when
this monotonic decrease occurs. In Problem 5.4 we will note a contrasting
Attractive forces: why colloids are sticky 119
case where f(x) oscillates instead. This f(x) can be related to the net
work W1 (ψ(0), ψ(x)) to change both ψ(0) and ψ(x). This W must be
minimal when ψ(x) = f(x)ψ(0). Also, since ψ is assumed small, we take
W to be quadratic in both ψ’s. One can easily deduce a W1 (ψ(0), ψ(x))
that is consistent with Eq. (5.8) and which gives ψ(x) = f(x)ψ(0) when
∂W1 /∂ψ(x) = 0:
W1 (ψ(0), ψ(x))
!
= W0 + W ψ(0) + 12 χ ψ(0)2 + ψ(x)2 − 2f(x)ψ(0)ψ(x) .
and
∂W
0= = W + χ(ψ(x) − f(x)ψ(0)).
∂ψ(x)
Since the two ψ’s play equivalent roles in these linear equations, they have
equal values:
−W
ψ(0) = ψ(x) = .
χ (1 − f (x))
Since f was positive, both ψ’s increase owing to their proximity. Moreover,
at this equilibrium value of the ψ’s, The free energy W is given by
W 2
W = Wmin = 2W0 − . (5.10)
χ (1 − f )
The nonlocal response f lowers the energy and thus induces an attractive
interaction between the two particles, as claimed.
Our proof covers the case where the perturbation on the medium is
weak. But perturbation leads to attraction in some cases of strong perturb-
ation, provided further conditions are met. Appendix A treats an important
category of strong perturbations.
One further category of perturbation not treated in our discussion is the
case where the perturbed quantity ψ is not a scalar. An important example is
the case where the particles induce a vector dipole moment in the molecules
around the particle. Such polarization often occurs when the solvent is water.
We suppose that the polarization vectors tend to point towards particle 1
and then estimate the energy to add particle 2. The dipoles nearest particle 1
120 Colloids
† 2,6-Lutidine or 2,6-dimethyl pyridine is a are now oriented away from where particle 2 is to be inserted. This adds to
small organic molecule whose chemical the W2 rather than diminishing it. Thus non-scalar perturbations need not
formula is C7 H9 N. Lutidine and water mix-
produce attraction.
tures undergo a continuous phase transition
into a lutidine-rich and a water-rich phase An example of the Perturbation-Attraction Theorem occurs when col-
at a convenient composition and temper- loidal particles are suspended in a mixed solvent, composed of a miscible
ature. At such a transition fluctuations of binary mixture such as lutidine and water† . Normally the particle surface
composition occur that have spatial extents prefers one of the solvents; this preference alters the concentration near
arbitrarily larger than the size of a molecule.
the particle. The concentration gradients extend a correlation length ξ from
Lutidine–water mixtures are often used for
studying the effects of phase separation in a the particle; when the composition and temperature approach the conditions
solvent. of phase separation, ξ becomes large. We refer the reader to a text on phase
transitions such as the book by H. E. Stanley [5]. Other particles are then
attracted to the regions of high favorable solvent concentration. As the The-
orem suggests, colloidal particles in mixed solvents tend to aggregate when
the solvent is near the threshold of demixing [6].
5.4. When perturbation causes repulsion Generally two like objects that perturb
a medium experience a mutual attraction. But sometimes they do not. Here
is a simple and important counterexample. The medium is a membrane with
a bending rigidity: the membrane resists curvature. Such a membrane whose
heightabove some reference plane is ψ has a free energy W [ψ] of the form
x W = dx Aψ 2 + B(d 2 ψ/dx 2 )2 provided ψ varies only with x. If dψ/dx is
h not too large, the second term is simply the square of the curvature. This is a
quadratic free energy and can be readily solved, though there is an awkward
integral to do. One finds that a perturbation ψ0 induced by a phantom vertical
plane at x = 0 with interaction energy V = λψ(0) has the form ψ0 (x) =
(a) Cλe−κz [cos(kx) + δ sin(kx)] for positive x. Evidently, ψ0 (−x) = ψ0 (x).
(a) Find the prefactor C and the coefficients δ, κ and k, that minimize W [ψ] +
V (ψ(0)). Note that d/dxψ|x=0 = 0 by symmetry.
(b) Now a second plane is inserted at distance x = h from the first. Using the
discussion in the text, find the energy W of the medium with ψ(0) and
ψ(h) fixed.
(b) (c) How does the total energy V + W vary with separation h? You can ignore
the overall magnitude of this energy. Are there separations where W + V
is larger than at h = ∞? If so, where?
excluded from the interplanar region, the pressures on each wall are bal-
anced. However, when the particles are excluded by their size from entering
the channel between the plates, the osmotic pressure on the outer surfaces
which is associated with the suspended particles is not counterbalanced,
leading to an effective attractive force per unit area between the plates of
magnitude . The range of this interaction is approximately the particle
diameter.
What happens when the colloidal particles of Fig. 5.2 are replaced by
flexible polymer coils? There are two distinct cases: the polymers are either
attracted to or repelled from the solid substrate. Which case is operative in h
any particular situation depends upon the three media—polymer, substrate,
and solvent—and the interactions between them. Polymer adsorption due
to attraction to the surfaces often leads to repulsive colloidal forces and Fig. 5.3
will be discussed in the next section. The non-adsorbing situation leads to Two dimensional representation of a
depletion attractions. Indeed, when the solution is dilute and the polymer polymer coil confined to the region
coils are well separated from one another, the system is quite analogous to between two plates for h < RG . The
that depicted in Fig. 5.2. In fact the osmotic pressure which drives the plates circles denote blobs of diameter h.
together at separations smaller than the polymer radius of gyration, RG , is
simply = (c/M)T , where c is the monomer concentration of chains
having degree of polymerization M. Of course, polymer coils are different
from colloidal particles in that polymers are deformable. Thus the chains
may distort and squeeze into the interplanar space even if the separation,
h, is small, h ≤ RG , as in Fig. 5.3. The free energy associated with this
deformation can be estimated (in a good solvent) as T per blob, as studied
in Problem 3.15. We discuss this confinement energy more carefully in
Chapter 6 to follow. Totalling over all the blobs in the channel yields a free
energy per polymer F given by
10,000
Force/radius (μNm–1)
1000
Fig. 5.7
Measured forces between two surfaces
100
coated with grafted polystyrene chains
about 1400 monomers long in toluene,
after [9]. The curves show predicted force
using free-solution properties of
10 polystyrene in toluene, together with the
0 200 400 600 800 1000 1200 1400 parabolic-profile theory described in the
Separation (Å) text, after [8].
arising from, e.g., a Na+ ion at the origin? The Na+ ion attracts Cl− ions
and repels other Na+ ions. Therefore if the test charge is at a large dis-
tance (compared to the average interionic spacing d ≡ c−1/3 ) from the
origin, Gauss’ Theorem would yield a reduction in the effective electro-
static potential: (r ). This reduction is called Debye–Hückel screening.
The form of the Debye–Hückel screening may be derived directly from
Maxwell’s electrostatic equation (written in cgs units)
where is the dielectric constant (80 for water), is the electrostatic poten-
tial and ρ(r ) = ρ+ (r ) − ρ− (r ) is the charge density at r; ρ± (r ) are the
anion and cation concentrations. In a solution, the charges are mobile and
may adjust to the local forces that they experience; the position of each ion
is governed by its Boltzmann probability owing to the local potential that
it experiences.
0
ρ± (r )
ρ± exp(∓ e(r )/T ), (5.16)
where ρ± 0 are constants that adjust the total number of ions of each species
to correspond to the concentration of ions dissolved, i.e., c = (ρ+ +
ρ− )(r ) d 3 r, where is the total volume of solution. We now suppose that
the potential seen by an ion is that arising from the mean charge density
ρ(r) of Eq. (5.15). Combining Eqs. (5.15) and (5.16) generates a nonlinear
second order partial differential equation for the self-consistent electrostatic
potential, (r ):
2 0 −e(r) 0 e(r)
∇ = 4π −ρ+ exp + ρ− exp . (5.17)
T T
κ 2 = (4πe/T )(ρ+
0 0
+ ρ− ). (5.19)
which is of the same form as the Yukawa potential that is familiar from nuc-
lear physics. This is a reasonable approximation if there are many screening
ions in a screening length, which corresponds to /d 1, as may be verified
using the formula for κ.
5.6. Debye–Hückel screening (a) Show that the condition that there are many
screening ions in a Debye length κ −1 is equivalent to /d 1 for monovalent
salts, where d is the mean spacing between ions in solution and is the Bjerrum
length. (b) What is the corresponding result for multivalent salts? (c) For a
CaCl2 solution in water, what is the range of concentrations for which the
Debye–Hückel approximation is valid?
A natural extension of Eq. (5.20) to the case of a colloidal particle of size
b and having Z ionizable groups firmly anchored to the surface is
−κ(r−b) 1
e(r) = ZT e . (5.21)
r κb + 1
However this is typically an overestimate of the strength of the repulsive
potential between charged colloidal particles. This can be seen by consid-
ering the potential energy of a negative counterion at the surface of the
particle, e(b) = −4π bσ T , where σ is the number of anchored charges
per unit area. For b
1 μ and σ
10−4 / nm2 , this gives |e(b)|
104 T .
There is therefore a strong attractive force driving the counterions to recon-
dense on the particles, so that the particle is incompletely ionized. This arises
because the strong, generally covalent anchoring of the surface charges
leads to very high local electric fields. This results in a reduction of the
effective charge, Z ∗ , of the particle as experienced by a distant test charge.
We can readily estimate Z ∗ for colloidal particles whose radius b is much
larger than the screening length κ −1 , the typical case. The surface charge
gives an enhanced density of counterions near the surface. The Boltzmann
enhancement factor is exp(e(b)/T ), where e(b)
Z ∗ T /b is the sur-
face potential. The enhancement acts over a volume of order b2 /κ. The
probability that the charge is in the enhanced region is the Boltzmann-
weighted ratio of the enhanced volume to the volume per ion 1/c. This
ratio is about 1 when eφ(b)
−T ln(cκ/(b2 )) or
d 2φ
− = (2λ2 )−1 e−φ , (5.23)
dz2
With this definition of λ, ρ 0 is the free charge density adjacent to the sur-
face. One readily verifies this by evaluating Eq. (5.23) at the surface, and
recognizing Poisson’s equation with charge density ρ 0 . We must solve the
equation in order to find this ρ 0 . This one dimensional Boltzmann–Poisson
equation has the solution,
This may easily be verified by inserting (5.24) into (5.23). We note that
φ → ∞ as z → ∞. The boundary condition which is satisfied at the surface
is given by Gauss’ Theorem and is ∂φ/∂z|z=0 = 4π σ . The corresponding
charge distribution is
ρ = ρ 0 [1 + (z/2λ)]−2 . (5.25)
two surfaces is 2h, then by symmetry there is no electric field on the mid-
plane, i.e., ∂φ/∂z|z=h = 0. Then the only force that can be transmitted is
associated with the osmotic pressure of the counterion “gas”:
= T ρ(h)/e. (5.26)
5.10. Rod excluded volume A slender rod of length L and radius b L is centered
at the origin and pointed along the z axis. A second identical rod lies along
the x axis at a distance x from the first. Its orientation is arbitrary.
(a) What is the probability that the second rod will intersect the first? Ignore
numerical prefactors and consider the dependence on L, b, and x. How
small must x become in order for the probability to be of order unity?
134 Colloids
(b) Using the results of (a), estimate the mutual excluded volume of two such
rods if they are dispersed in a solution. Again, ignore numerical prefactors.
(c) Repeat parts (a) and (b) with the rods replaced by disks of diameter L
and thickness b. Compare with the result for the rods. Compare with the
result for hard spheres of diameter L found in Chapter 4.
U1,2 = r −3 [μ
1 · μ
2 − 3(μ 2 · r)r −2 ].
1 · r)(μ (5.28)
5.4.1 Electrophoresis
Charged colloidal particles move in a distinctive way when an electric field
acts on them. The contrast with conventional forcing is greatest when there
are many free ions and the screening length is short. In Chapter 4, we
analyzed the response to external forces by tracing the added momentum as
it flowed outward from the forced particle. The long-ranged flow field of the
surrounding fluid leads to hydrodynamic drag. However, when a charged
particle is forced by an electric field, the situation is radically different.
In any unit of time, the momentum added to the particle is balanced by
opposite momentum added to the counterions in the oppositely charged
screening layer. For a uniformly charged body that is smooth on the scale of
a screening length, the two momenta cancel nearly completely, and virtually
no momentum is transferred to the fluid. Thus no long-range flow occurs.
The particle crawls through the liquid rather than swimming through it. All
the dissipation leading to the retarding force occurs in the narrow screening
zone. Thus the velocity of each bit of surface is determined independently of
the rest. The electrophoretic mobility depends on the surface charge density
and the screening length, but not on the size of the particle.
When the charge on a particle is nonuniform, the phenomenon of electro-
phoresis becomes startlingly rich and complex, as recent discoveries have
shown [25]. If two regions on the particle have different charge densities,
the electrophoretic force wants to pull them at different speeds. Since the
two regions are obliged to move at the same speed, one region must exert a
force on the other. This means that each region experiences a net unbalanced
force, which it must give to the surrounding fluid, producing hydrodynamic
backflow. An opposing momentum is injected at the other region. A net
momentum current comparable to the internal force is injected into the
fluid, and long-range Stokes flow occurs. One may control the nature of
the induced motion in striking ways by controlling the shape of the particle
Appendix A: Perturbation attraction in a square-gradient medium 137
and the placement of charge on it. One may induce motion at right angles
to the electric field or rotation with no translation, for example [25].
h 1 dψ 2
W [ψ] = 2 dz w(ψ(z)) + m(ψ(z)) . (5.29)
0 2 dz
Here the w(ψ) is the free energy per unit volume required to make a
uniform change from ψ0 to ψ . We have restricted the integration to the
left half of the system, noting that the contribution from the right half must
be identical. We shall suppose that the variation of ψ is gentle, so that we
† There is no piece linear in this gradient. may neglect higher powers of the gradient dψ/dz† . Its coefficient m may
Such a piece would not be invariant under depend in general on the local value of ψ †† .
a change of coordinate system that replaces
Two identical surfaces immersed in such a square-gradient medium must
z by −z; it is thus ruled out by the require-
ment that the free energy density must be attract each other at long distances. Further, the attractive force per unit
independent of coordinate system. area is simply the local energy density at the midplane wm . To see this, we
††
separate the surfaces by a small amount 2h and examine the change in free
As in the polymer solution, any fluid
with a free energy of this form has spatial
energy F . We perform the separation in two stages. First we separate the
correlations in the fluctuating equilibrium surfaces without allowing a change in the ψ profile. We simply extend its
ψ(z). These correlations die off exponen- midpoint ψm over the extra interval h on each side of the system. Second
tially in space with a decay length ξ given we allow the ψ field to relax to minimize F .
by ξ 2 = d 2 w/dψ 2 |0 /m(ψ0 ). In the first stage F changes by an amount 2hL2 wm . Since ψ is constant
here, the gradient energy is zero. In the second stage, ψ changes in the pre-
existing region and in the newly created gap region. The resulting change
of F in the gap region is second order; it is the product of the small change
ψm and the small interval 2h. The profile in the old region from 0 to h
was that which minimized the free energy F with the constraint that dψ/dz
vanish at h. This derivative vanished because h was at the midpoint of the
full system. To examine how this profile changes, we first determine the
optimal profile between 0 and h without constraints. We allow dψ/dz(h) to
vary freely. But the optimal value of dψ/dz(h) is zero in a square-gradient
medium. We may check this by dividing the interval up into many equal,
discrete intervals ending at z1 , z2 , . . . , zk = h. The corresponding ψ values
are ψ1 , . . . , ψk . The free energy F is the sum of contributions from each
interval. To minimize F in the i’th interval we need only adjust ψi and
ψi−1 . The only interval that involves ψk is the last one. The free energy
contribution to it is
−2
h h h
w(ψk−1 ) + m(ψk−1 )(ψk − ψk−1 )2 .
k k k
Since ψk appears only in the gradient term, the optimal value for it is
that which makes the gradient vanish. Minimizing F automatically makes
the gradient vanish at h. The vanishing of the gradient is not an additional
††† The author is grateful to Alexei constraint††† . When we allow ψ(z) to relax to the new, wider interval, there
Tkachenko for pointing this redundancy out is no proportionate change in F . The original ψ(z) already minimized F
to him.
in the old region 0–h with respect to arbitrary variations. This means any
small change ψ can make at most a second-order change in F in that
region. Thus the only change in F is the contribution from the gap region,
Appendix B: Colloidal aggregates 139
d(F /L2 )
= wh .
d(2h)
Though we have assumed our surfaces to be identical and flat, this force
law holds more generally. If two surfaces are sufficiently similar, there must
be a point between them where the gradient of ψ vanishes. Then the above
argument can be extended to show that the force per unit area at that point
is attractive and equal to w.
While this perturbation attraction occurs quite broadly, there are import-
ant cases where the mechanism does not apply. One is the case of surfaces
with grafted polymers treated in the chapter. Here polymers are chosen
which are attracted to the surfaces only at one end. The attraction results
in a polymer-enriched layer near each surface. When the two surfaces are
brought together so that the layers interpenetrate, a repulsion occurs, not
the attraction implied by the Perturbation-Attraction Theorem. Why is the
theorem not applicable? The reason is that the free energy does not have the
square-gradient form, as in Eq. (5.29). Instead, the energy cost of a nonuni-
form concentration is nonlocal. The cost at position z does not depend
merely on ψ at z and on ψ at neighboring z. It depends on ψ throughout
the layer. Thus e.g., if some chains were removed from the layer, the free
energy at the midpoint h would be affected—even if some external agent
forced ψ to be fixed in the intervening layer. Thus the free energy cost of
nonuniformity cannot be found by adding up gradient contributions from
the whole layer. Nonlocal differences such as ψs − ψ(h) also contribute
directly to the free energy.
Another case where perturbation does not lead to attraction is the
deformation of a membrane with bending rigidity, as treated in Problem 5.4.
nature of the dense phase and the amount of attraction or colloid concen-
tration required to produce it are described by the statistical mechanics of
phase transitions [5]. This well-developed branch of theoretical physics
gives powerful, specific, and well-verified predictions about how mutually
attracting particles undergo phase separation. We encountered this phe-
nomenon briefly in Chapter 4, in discussing the “collapse” of a polymer in
a poor solvent. For the most part the theory of phase transitions applies as
well to colloidal particles as it does to small molecules. Thus the nature of
precipitation and creaming is virtually the same as, e.g., the precipitation of
steam or the liberation of carbon dioxide from a freshly opened soft drink.
But in some cases the distinctive attraction of colloidal particles results in
qualitatively new phenomena.
Phase transition theory applies only when the particles have come to
thermodynamic equilibrium, at least locally. As we saw in Chapter 2, equi-
librium is the state attained after a sufficiently long time. For any change
of configuration the system has made, it must also have time to make the
reverse of that change. For example, any given pair of particles that are
together should have had enough time to separate and recombine several
times in order to attain equilibrium. With small molecules, phase separation
is usually slow enough that each local region (containing a few molecules) is
close to thermal equilibrium. But in colloids, one encounters a new regime
(a) of irreversible attraction, far from equilibrium.
(b)
As we have noted, most forms of colloidal attraction grow in strength with
the particle size. For large particles, the attractive energy may far exceed
the thermal energy T . Then the attraction becomes for practical purposes
irreversible, as discussed in Chapter 2. Two particles that stick together
(c) (d) have virtually no chance to unstick in the time the system is observed:
the strength of the attraction prevents even local equilibration. Thus phase
transition theories cannot describe the assembly process; a new approach
is needed.
The assembly process in this irreversible regime is called kinetic aggrega-
tion (Fig. 5.12). Kinetic aggregation of a colloidal solution ultimately
Fig. 5.12 produces large aggregates like the one shown in Fig. 1.5. These are strik-
Nonequilibrium aggregation. (a) Typical ingly different from the compact droplets or grains of normal precipitation:
contacting configuration of three solid the aggregates are open, tenuous assemblies. The nature of these assem-
colloidal particles. The angle between the blies has been much studied in the last two decades and has been treated in
three is arbitrary. (b) Configuration of
comprehensive reviews [29, 30].
lower attraction energy. In thermal
equilibrium this state is much more The colloidal particles of Fig. 1.5 were initially dispersed in a very dilute
probable than (a), but in irreversible state. Four-nanometer charge-stabilized silica particles were suspended in
aggregation it is no more probable. water at a volume fraction of roughly 10−6 [31]. Then a large amount
(c) Typical contacting configuration of of salt was added, screening out the repulsive Coulomb barrier between
three triangular clusters of liquid droplets.
the particles nearly completely, and exposing the particles’ strong van der
(d) Compact configuration. Like
configuration (b), this state is unreachable Waals attraction. A few minutes after the salt was added, the aggregate had
in the available time. Liquid drops can go formed. When two particles encountered each other in the course of their
from state (a) to state (b) with no energy ordinary Brownian motion, they stuck together permanently. The resulting
barrier, but not from state (c) to state (d) pair continued its Brownian motion until it encountered another particle or
[28]. To pass from (c) to (d) the shaded
cluster. These two clusters stuck together permanently in their contacting
ball would have to break the contact
indicated by the arrow and overcome the configuration and continued to move. In this way the clusters grew to the
associated energy barrier. size shown. One may readily simulate this process on a computer [29], to
Appendix B: Colloidal aggregates 141
1
Gold
Fig. 5.13
Silica Structure factors S(q) from colloidal
Polystyrene aggregates like that pictured in Chapter 1,
0.1
obtained from light scattering, after [31].
The upper three data sets marked DLCA
DLCA were from aggregates grown under
diffusion-limited conditions using 7.5 nm
I (q)
0.01
gold particles, 3.5 nm silica particles, and
RLCA 20 nm polystyrene spheres. Straight lines
have the slopes expected for a fractal
0.001 structure with D
1.84. The lower three
data sets marked RLCA were from
aggregates grown under reaction-limited
conditions made from the same three
0.0001 dispersed colloids. Straight lines have
0.002 0.005 0.01 0.02 0.05 slopes expected for a fractal structure with
q (nm–1) D
2.1.
larger cluster. The two joined clusters are typically unequal, and they may be
greatly unequal. The second drastic simplification is that self-intersections
are ignored. We may see the effect of both of these simplifications by simple
modifications of the model.
It is easy to modify the Ghost model to account in a primitive way for
unequal clusters. This modification is called the “fixed-ratio Sutherland’s
Ghost” model [32]. This model is like the original model, except that the
reacting clusters are forced to have a particular mass ratio, such as 3 : 1.
One may readily modify the above reasoning to determine the average
bond distance, B4n , in terms of those of its constituents, Bn and B3n . Now
when two particles are chosen at random to compute the combined bond
distance, B4n , the probability that both particles are on the smaller cluster is
now reduced (to 1/16); the other probabilities are also altered. For a general
ratio r : 1, the probability that an arbitrarily picked point belongs to the
small cluster is 1/(r + 1). The probability for the larger cluster is r/(r + 1).
Eq. (5.30) becomes
2 2
1 r r
B(r+1)n = Bn + Brn + 2 (Bn + Brn + 1).
r +1 r +1 (r + 1)2 † The probability P (n, t) that a given cluster
has mass n at time t is called the mass
This equation also allows power-law solutions of the form Bn ∼ nx . distribution. To account for the change of
Substituting into the above equation, one finds an implicit equation for x. this distribution as the aggregation pro-
ceeds, one must consider all processes
2 2 which may increase or decrease the num-
1 r r
(r + 1) =
x
+ rx + 2 (1 + r x ). (5.31) ber of n-clusters. The joining of an i cluster
r +1 r +1 (r + 1)2 with an (n − i)-cluster increases P (n, t).
The joining of an n-cluster with any other
The exponent x decreases with increasing ratio r. For small r, Eq. (5.31) cluster decreases P (n, t). The rate of join-
reduces to ing of i-clusters with j -clusters is evidently
proportional to their numbers, P (i, t) and
x = 2(1 + r x ).
to P (j , t). But two different mass pairs
As in the simple Ghost model, the geometric path between any two sites on a with the same cluster numbers P need not
cluster is still a random walk, so that n ∼ R D with D = 2/x. Evidently the have the same joining rate. A given pair
of clusters with masses i and j has a spe-
nominal fractal dimension D goes to infinity as the ratio r → 0. It is possible cific joining rate, denoted K(i, j ). Adding
to extend this treatment to include a complete distribution of aggregating together all the processes that change the
sizes, not merely a fixed ratio [33]. For the expected mass distributions† , D number of n-clusters, one obtains the so-
increases moderately from the 3.4 of the simple Sutherland’s Ghost model. called Smoluchowski equation [34]:
Clearly, to account for D of kinetic aggregates quantitatively requires a ∂P (n, t) 1
realistic treatment of the relative masses of the aggregating clusters. = K(i, j )P (i, t)P (j , t)
∂t 2
i,j
The second glaring defect of the model is its neglect of self avoidance. We
may evaluate the importance of this defect in the same way we have done ×[δi+j ,n −δi,n −δj ,n ].
for polymers in Chapter 3. We modify the process to assure self-avoidance.
The Smoluchowski equation gives several
Pairs of clusters are joined as in the original model. Then the combined different types of P distributions, depend-
cluster is checked for self-intersections. If any are found, the combined ing on the specific joining rate K. The
cluster is discarded. If this discarding probability approaches unity for large equation itself is an approximation, for
clusters, the average properties of the remaining ones may be qualitatively it accounts for only the masses of the
affected. But if the discarding probability does not approach unity, the clusters and ignores other variables, such
as the cluster shapes and their positions
average properties cannot be qualitatively affected and scaling exponents in space. Nevertheless, it is believed to
such as D must be unchanged. be well-justified for three-dimensional col-
By this reasoning, we may see that self-avoidance has no impact on D loidal solutions. These subjects are treated
in sufficiently high spatial dimensions d. We saw in Chapters 3 and 4, that at length in [34].
144 Colloids
when two fractals of size R are placed at random in the same volume, the
average number of intersections between them within distance R goes as
R D1 +D2 −d . When this exponent is negative, the number of intersections
goes to zero as R → ∞. Thus our two aggregates with dimension D have
no intersections in spatial dimension d > 2D. Even if the two aggregates
are connected at the origin, the number of intersections is limited and does
not grow indefinitely with R. Thus in these high spatial dimensions the
probability of self-avoidance is finite and does not go to zero as the cluster
sizes go to infinity. Then the aggregation process is not qualitatively affected
by the self-avoidance constraint, and the D of the aggregates is unchanged.
For Sutherland’s Ghost aggregates with D
3.4, self-avoidance has no
affect on D above d = 6.8. If the aggregation process were carried out on
a computer in a virtual space of seven or more dimensions, the resulting D
should be 3.4 even with self-avoidance imposed. We have encountered the
analogous property for polymers: their self-avoidance does not affect their
fractal dimension in more than four spatial dimensions.
In lower spatial dimensions, self-avoidance becomes relevant. In poly-
mers the discarding process has the greatest impact on the more compact
configurations. The remaining, self-avoiding, ones are generally increased
in size, and their D is reduced. The same qualitative behavior is expected
for aggregates. D should decrease progressively as the spatial dimension d
is reduced.
does not affect the conclusions.) The resulting object is a fractal whose
dimension is two greater than that of the cluster itself. We now choose a
random point on this cluster and on another cluster. We then join these
according to the Ghost model. However, for a valid joining configuration,
the history must not intersect the other cluster. If it does, it must be discarded.
The remaining configurations are valid joining configurations, in which the
one cluster has avoided the history of the other† . † In making this construction we have
This constraint is again that of the mutual avoidance of two fractals. treated one cluster as stationary and the
other one as moving in a random walk. This
The Brownian-motion history of a fractal of dimension D is a fractal of
amounts to working in the frame of refer-
dimension D + 2. The history of each particle of the fractal lying in a ence of one cluster. This is equivalent to
sphere of radius R is a random walk, with roughly R 2 particles in the sphere. the actual situation in which both clusters
The number of particles leaving these random-walk trails is the number of execute random walks.
particles of the fractal—roughly R D . The total number of particles in the
history within the sphere is thus roughly R 2 R D = R 2+D . As before, it must
be irrelevant in sufficiently high spatial dimensions d. If the clusters have
dimension D, then the avoidance has no effect on D provided (D + 2) +
D − d < 0. For our Ghost aggregates, this condition holds for d 8.8.
Evidently this random-walk avoidance has a stronger effect than the simple
avoidance discussed previously. D begins to decrease from its asymptotic
value for higher d (8.8 versus 6.8), and at a given d we would expect it
to decrease further. When clusters move by Brownian motion and stick
on contact, the aggregation is called diffusion-limited. This form is very
prevalent in colloidal solutions. The aggregate of Fig. 1.5 was produced in
diffusion-limited conditions.
Sometimes diffusion is not the rate-limiting aspect of the aggregation
process. It can happen that particles stick irreversibly with very low prob-
ability even when they are adjacent. Indeed this is the case for stabilized
colloidal particles. If the stability is diminished somewhat, e.g. by adding a
little salt, the sticking occurs fast enough to observe. Still the sticking rate
may be much slower than the rate of encounters via Brownian motion. This
regime is called reaction-limited aggregation. In the reaction-limited regime,
two clusters which join at a particular moment have had ample opportunity
to visit each other’s neighborhood. Thus all configurations that do not inter-
sect are equally likely, as in equilibrium. (Here we ignore the local repulsion
that inhibits very close approaches; we implicitly count such approaches
as intersections.) Now the joining proceeds as in the self-avoiding Suther-
land model discussed above: the joined configurations are taken at random
from all contacting configurations. As we have seen, this model is less
constrained than the history-avoiding diffusion-controlled case. Thus we
expect the resulting aggregates to have a fractal dimension D closer to the
large D = 3.4 of the unconstrained model.
The diffusion-limited and reaction-limited regimes are the most import-
ant types of aggregation seen in colloidal suspensions. But other types of
motion are also of interest. In ballistic aggregation, the joining clusters
follow straight-line paths instead of random walks. One form of ballistic
aggregation is sedimentation aggregation, in which the joining clusters are
drifting downwards in a gravitational field. They join when a heavier cluster
overtakes a lighter one [35]. Large, composite snowflakes are made by this
form of aggregation.
146 Colloids
4
D=d D = d/2
Not Transparent
tenuous to self
D
×o + + + + +
o? +◊
◊?
+
1 +
o Reaction-limited aggregates, after [31]
◊ Diffusion-limited aggregates, after [31]
+ Good-solvent polymer
Fig. 5.16
× Randomly branched polymers
Fractal dimension D of colloidal
aggregates and other fractals versus spatial 2 4 6 8
dimension d. d
Properties of aggregates
The fractal structure of colloidal aggregates makes them broadly useful, as
noted in Chapter 1. The aggregates are used like polymers to thicken liquids
and solidify them by producing a gel network. They are used to reinforce
rubber and make it tougher. The origin of these special mechanical prop-
erties can be readily understood, adapting the reasoning used in Chapter 3
and 4.
A small volume fraction of colloidal aggregates in a liquid increases its
viscosity markedly. The reason is the same as in a polymer solution. The
aggregates, like polymers, are opaque to flow since they have D > 1: a
shear flow goes around rather than through them. Thus they impede flow as
though their pervaded volume were essentially filled in with solid material
instead of being nearly pure solvent.
Aggregates like polymers produce osmotic pressure in dilute solution.
Since these fractals are opaque to each other, their mutual excluded volume
Appendix B: Colloidal aggregates 147
i E i . Each θi , acting on its own, would produce a displacement ui = θi ri . Local bending of two particles in a
spanning arm caused by displacement of
Here ri is the distance from the bead i to the displaced end. Clearly, beads
one end. The undistorted configuration is
further from the end produce the displacement more efficiently. The bulk of shown in light shading. The two lines,
the bending will involve beads at distances of order R from the end. All such drawn transverse to each bead, were
beads will have a comparable share of the bending and comparable angles parallel in the unbent configuration.
148 Colloids
θi . This θ will be such that u
i θ R. The number of beads involved
is roughly the number of beads Ma in the arm. Thus θ
u/(Ma R). The
resulting energy E is given by
E= Ei
kMa θ 2
ku2 /(Ma R 2 ).
i
This expression defines the desired spring constant K(R), the coefficient of
u2 . Using the fractal law for Ma , we find
K(R)
kR −C−2 .
From this basic elastic constant K(R), we may infer the bulk elasticity of
a mass of aggregates at volume fraction φ. When an aggregate suspension is
compressed to produce this elasticity, the individual aggregates have begun
to press against each other. The aggregates are at their overlap concentration:
φ = φi
(R/a)D−3 . The longest spanning arms of a typical aggregate
are deformed by these external contacts. If a small strain γ is now applied,
then these longest arms are displaced by an amount u
γ R. Each stores
an energy
E = Ku2
kR −C−2 γ 2 R 2 .
Each aggregate has only a few of these spanning arms: their number does
not increase with R. Thus the energy stored per aggregate volume is of order
E/R 3 . Combining, we find that the overall strain γ produces an energy per
unit volume of order kR −C−3 γ 2 . The coefficient of γ 2 defines an elastic
modulus G for the gel. This G, expressed in terms of volume fraction, obeys
G
kφ (C+3)/(3−D) a −3 . (5.32)
Figure 5.18 shows a comparison of the predicted power and
experimental data.
0.1
Modulus (GPa)
0.01
Fig. 5.18
Modulus versus concentration in two
different aggregate materials, after [40]
(silica aerogel, filled dots) and [41], 0.001
(fumed silica, open rectangles). These
experiments approach the percolation
threshold at low concentrations; this is
0.0001
expected to depress the moduli below the
predicted power law. The straight lines,
with slopes 3.2 and 3.55, indicate the
power law predicted by Eq. (5.32), 0.00001
including the experimental uncertainties in 0.01 0.1 1
D and C [38, 42]. Volume fraction
References 149
References
1. P. Pincus, in Lectures on Thermodynamics and Statistical Mechanics: XVII
Winter Meeting on Statistical Physics, ed. A. E. Gonzalez and C. Varea
(Singapore: World Publishing Co., 1988).
2. D. H. Everett, Basic Principles of Colloid Science (London: Royal Society of
Chemistry, 1988).
3. J. N. Israelachvili, Intermolecular and surface forces, 2nd ed. (London: San
Diego, CA: Academic Press, 1991).
4. J. Mohanty and B. W. Ninham, Dispersion Forces (London: Academic Press,
1976).
5. H. E. Stanley, Introduction to Phase Transitions and Critical Phenomena (New
York: Oxford University Press, 1971).
6. F. Brochard and P. G. de Gennes, Ferroelectrics 30 33 (1980).
7. For an excellent review of polymer brushes see A. Halperin, M. Tirrell and
T. Lodge, Adv. Polym. Sci. 100 (1991).
8. S. T. Milner, Europhys. Letts. 7, 695 (1988).
9. H. J. Taunton, C. Toprakcioglu, L. J. Fetters, and J. Klein, Nature 332 712
(1988).
10. A. N. Semenov, J.-F. Joanny, A. Johner, and J. Bonet-Avalos, Macromolecules
30 1479 (1997).
11. P. G. de Gennes, Scaling Concepts in Polymer Physics, (Ithaca: Cornell
University Press, 1979).
12. S. Bucci, C. Fagotti, V. Degiorgio, and R. Piazza, Langmuir, 7 824, (1991).
13. S. H. Behrens and D. G. Grier, in Electrostatic Effects in Soft Matter and Bio-
physics, ed. C. Holm, P. Kekicheff, and R. Podgornik (Dordrecht: Kluwer,
2001).
14. M. Alexander, L. F. Rojas-Ochoa, M. Leser, and P. Schurtenberger, J. Colloid
Interface Sci. 253 35 (2002).
15. R. Kjellander, T. Akesson, B. Jonsson, and S. Marcelja, J. Chem. Phys. 97 1424
(1992).
16. L. Guldbrand, B. Jonsson, H. Wennerstrom, and P. Linse, J. Chem. Phys. 80
2221 (1984).
17. P. Pieranski, Contemp. Phy. 24 25 (1983).
18. See e.g., M. Constantinos, Paleos ed., Polymerization in Organized Media
(London: Taylor and Francis, 1992).
19. M. J. Stevens, and M. O. Robbins, J. Chem. Phys. 98 2319 (1993).
20. L. Onsager, Ann. N.Y. Acad. Sci. 51 627 (1949).
21. P.G. de Gennes, and J. Prost The Physics of Liquid Crystals, 2nd ed. (Oxford :
Clarendon Press, 1993).
22. R. E. Rosensweig, M. Zahn, and R. Shumovich, J. Magn. Magn. Mater. 39
127(1983); R. E. Rosensweig, Ferrohydrodynamics (New York: Cambridge
Univ. Press, 1985).
23. A. P. Gast and C. F. Zukoski, Adv. Colloid Interface Sci. 30 153 (1989).
24. E. J. Hinch and J. D. Sherwood, J. Fluid Mech. 132 337 (1983).
25. D. Long and A. Ajdari, Phys. Rev. Lett. 81 1529 (1998).
26. J. L. Anderson, Ann. Rev. Fluid. Mech. 21 61 (1989).
27. This treatment extends P. G. de Gennes, Macromolecules 15 492 (1982); 14
1637 (1981).
28. P. Poulin, J. Bibette, and D. A. Weitz, Eur. Phys. J. B 7 277 (1999).
150 Colloids
29. P. Meakin, in Phase Transitions and Critical Phenomena, Vol. 12 C, eds. Domb,
J. L. Lebowitz (New York: Academic, 1988) p. 335.
30. See e.g. R. Jullien and R. Botet, Aggregation and Fractal Aggregates
(Singapore: World Scientific, 1987).
31. M. Y. Lin, H. M. Lindsay, D. A. Weitz, R. C. Ball, R. Klein, and P. Meakin,
Nature 339 360 (1989).
32. R. C. Ball and T. A. Witten, J. Stat. Phys. 36 873 (1984).
33. R. Botet, J. Phys. A: Math. Gen. 18 847 (1985).
34. P. G. J. VanDongen and M. H. Ernst, J. Phys. A: Math. Gen. 18 2779 (1985).
35. P. Meakin, Rev. Geophys. 29 317 (1991).
36. R. Buscall, P. D. A. Mills, J. W. Goodwin, and D. W. Lawson, J. Chem. Soc.,
Faraday Trans. I 84 4249 (1988).
37. T. A. Witten, M. Rubinstein, and R. H. Colby, J. Phys. (Paris) II 3 367 (1993).
38. P. Meakin, I. Majid, S. Havlin, and H. E. Stanley, J. Phys. A 17 L975 (1984).
39. Y. Kantor and I. Webman, Phys. Rev. Lett. 52 1891 (1984).
40. T. Woignier, J. Phalippou, Revue de Phys. Appliquee 24 C4-179 (1989).
41. J. Forsman, J. P. Harrison, and A. Rutenberg, Can. J. Phys. 65 767 (1987).
42. P. Dimon, S. K. Sinha, D. A. Weitz, C. R. Safinya, G. S. Smith, W. A. Varady,
and H. M. Lindsay, Phys. Rev. Lett. 57 595 (1986).
6
Interfaces
Most of our experience with liquids arises not from the bulk of the liquid but
from its interfaces with other media. As we gaze onto the ocean or into our
cup of coffee, it is the surface of the liquid that catches our eye. The property
that distinguishes a liquid from a gas is the existence of this surface. This
chapter deals with such interfaces—how they are structured, how much
energy they have, and how fast they move and relax. Most importantly, it
deals with how solute molecules, both simple and complex, affect these
properties.
A great deal of the importance of structured molecules in liquids arises
from their affect on interfaces. Molecules that migrate to an interface and
affect its properties are called surfactants. Often such molecules are added
to fluids like motor oil, cosmetics, cleaning products, and food, in order to
achieve desired interfacial properties. In rubber and plastic composites like
automobile tires and appliance cases, the strength of the material hinges
on the properties of the interface between the components making up the
composite. Discoveries involving the interface between water and fabric or
solids and internal body tissues have had recent technological impact. But
the most dramatic property of surfactants is their ability to cause spontan-
eous production of interfaces. Such surfactant-generated interfaces are the
subject of the next chapter.
We begin this chapter by describing what the quantitative properties mean
in experimental terms. We survey ways of measuring energy, structure, and
motion associated with interfaces. Next we examine how interfacial energy
comes about for simple liquids. We discuss wetting of an interface by a
liquid, including the motion of the liquid in the process of wetting. Next we
turn to the effect of added solute molecules on an interface. Such molecules
can alter the energy of an interface and create structure there. Finally we turn
to structured fluids like polymers. These create structures at the interface
with distinctive scaling properties extending over distances as large as the
polymers.
Plate
Air
Tube of
Light source Drop of
heavy liquid
light liquid
Liquid
Fig. 6.1
Devices for measuring interfacial properties. (a) contact angle micrometer after [1]. Schematic light bulb in foreground illuminates the black
droplet on the white sample plate from the side. Lens system behind the sample plate projects the shadow of the drop. Contact angle is measured
from this shadow. (b) Spinning drop tensiometer, after [2], used for measuring small interfacial tensions. Horizontal tube at center rotates about
horizontal axis. Centrifugal forces elongate the black drop of lower-density fluid. The microscope at the top views the distorted drop shown and
measures its dimensions. (c) Wilhemeny plate, pulling upwards on a liquid surface. The support F is attached to an analytical balance to measure
the force on the plate. A motor, not shown, raises or lowers the fluid so that the fluid film meets the plate vertically.
The idea is sketched in Fig. 6.4. The incident light is arranged to send an
evanescent wave along the surface in question. The light intensity decays
exponentially with distance from the interface, and the decay length may be
controlled within limits. This light excites the adsorbate molecules, which
have to be labeled=† with a fluorescent dye. By measuring the intensity
of the fluorescence as the evanescent decay length is varied, one may †That is, a fluorescent chemical group must
obtain information about how the fluorescent molecules are arranged in be attached to the molecule to be probed.
space [4].
Often one wants to know about structure at the single nanometer
scale. Here x-ray and neutron scattering are used. To infer the profile
of adsorbate density with distance from the surface, one may meas-
ure reflectivity near conditions of total internal (or external) reflection.
As the incident beam angle is increased, the fraction reflected begins
to drop from 100% as transmission becomes possible. If the sample is
a thin film, the reflectance variation with angle can be used to infer
the composition as a function of depth. The inference is not straight-
forward, though. One must guess the profile and keep improving the
Fig. 6.4
guess until the calculated reflectivity versus angle matches the measure- Selective detection of molecules near a
ments. surface via evanescent wave fluorescence,
Quite often molecules are not spread uniformly over an interface, but after [4]. Beam of light hits the interface
are heterogeneous or patterned in some way. This kind of spatial struc- at a low angle and is totally internally
reflected, entering the liquid only in the
ture can make a big difference in the properties of the interface, yet
form of an evanescent wave near the
this in-plane structure is not apparent from the depth profiles treated interface. Any fluorescent molecules in the
above. To find such information from scattering, one must probe scattered evanescent layer emit light of a lower
wave vectors with components parallel to the surface. One must detect frequency. The intensity of fluorescent
not just the reflected waves, but waves scattered at an azimuthal angle light is measured by a photodetector, not
shown.
to the incident beam Fig. 6.5. Such experiments can detect moving
inhomogeneity as well as static patterns, using the methods of dynamic
scattering sketched in Chapter 4. Naturally, since the scattering occurs
from only a small number of molecules at the surface, the scattering is
weak and hard to detect. Thus these measurements require large-scale
neutron or x-ray facilities that are based on nuclear reactors and particle
accelerators.
Another way to see irregular coverage of an interface is to use microscopy,
as we previously discussed in Chapter 2.
Fig. 6.6
A snapshot of a fluid near a hard wall. Left
picture shows a “phantom” wall that does
not alter the particle arrangements. In the
center picture, local configurations in
which particles intersect the wall have
been removed. In the right hand picture a
hard slab has been inserted into the liquid,
distorting the particle configurations.
the interfacial energy per molecule in a liquid is also of order T or less. But
interfacial energy is significant even when there is no attraction between the
atoms. The boundaries of any fluid have a spatial structure and an energetic
cost. If the fluid fills a box, then these fluid boundaries are at the walls of
the box. We may account for the effect of these walls using the statistical
mechanics principles introduced in Chapters 2, 3, and 4. We consider the
simplest possible wall, a plane which adds an infinite energy cost to any
particle intersecting it. We call this a hard wall. The effect of such a wall
is to remove any configurations in which wall intersects particles. The
result is something like the middle of Fig. 6.6. Clearly the wall disturbs the
arrangement of particles in its vicinity. The density near the wall is reduced.
The disturbance extends roughly a particle diameter into the liquid.
To add a hard wall to a liquid requires work. As we know from Chapter 2,
the work can be found by considering properties of the pure liquid. In a pure
liquid there is probability pA that a given region A has no particles passing
through it. As told in Chapter 2, the work required to remove the particles
from such a region is −T log pA . If the region happens to be the slab shown
on the right side of Fig. 6.6, then this is the work required to insert the slab.
If two slabs were inserted at separate places in the fluid, the probability that
each region is empty is pA , independent of the other region. The probability
that both regions are empty is pA 2 , and the work required is just double that
for one of the slabs. That is, the work is proportional to area of the slabs. The
same is true if two large slabs are joined to form a single slab of twice the
area. Then almost all of each slab is far from the other slab, and independent
of it. Thus the work to insert the slab has the form α(area), where (area)
means the area (on both sides) of the slab. The coefficient α is called the
interfacial energy or interfacial tension. If the liquid is a hard sphere liquid
like the one pictured, the probability pA is some number of order unity if the
slab is no larger than a particle radius a. That means the interfacial tension
α
T /a 2 for such a liquid.
Real liquids generally have interfacial energies of this order, as well. The
fact that the molecules attract or repel each other somewhat does not change
the picture above qualitatively. Likewise, if the interface has some attrac-
tion or repulsion for the molecules, or some mild curvature or roughness,
this does not change the interfacial energy greatly. Even the free surface
between a liquid and its vapor is qualitatively similar. Like a hard wall, a free
surface is a surface where the density must be a small fraction of the liquid
density. Thus the disturbance it causes is qualitatively as great as that of a
hard wall.
Simple fluids 155
3 Fig. 6.7
Surface energies for the liquids tabulated
in [5]. Horizontal scale is the surface
energy, in millijoules per square meter (or
2.5
Scaled dynes per centimeter). Molecules with NH
surface or OH groups are denoted by squares;
energy others are denoted by diamonds. Vertical
2 axis is the scaled surface energy
α/(T /δA). This δA is the area per flexible
segment of the molecule. It is calculated by
1.5 first determining the volume of a flexible
segment from its mass m and the liquid
density. Then δA is the projected area of a
sphere containing that volume. For small,
1
rigid molecules, the flexible segment is the
whole molecule. For flexible molecules, it
is the smallest freely jointed unit larger
0.5 than a C2 H2 unit. For cases where the
flexible unit is taken as smaller than the
Surface energy whole molecule, an arrow extends upward
0 to indicate the range of values resulting
0 20 40 60 80 from choosing larger pieces of the
molecule as the flexible unit. Upper limit
of the range is indicated by a ceiling on its
arrow. Arrows with no ceiling are
Figure 6.7 shows some representative interfacial energies. As argued polymers. The lowest α’s shown are for
above, liquid interfacial energies are of the order of T per surface degree liquid helium at 4◦ K, m = 4 and liquid
of freedom. For small rigid molecules, this amounts to T per molecular helium at 1.6◦ K, m = 4. The remaining
area. For large flexible molecules, there is a characteristic area for each liquids, at room temperature, are, from left
to right, n-pentane, C5 H12 , m = 72.15,
independently moving piece† , which displaces a characteristic area at the polytetrafluoroethylene (Teflon),
surface. Then the interfacial energy is roughly T for each such area. Within (C2 F4 )n , m = 100, n-octane,
this generalization, polar molecules such as water tend to have more energy C8 H18 , m = 29, ethanol,
than nonpolar ones. Solids as well as liquids have interfacial energy. To C2 H5 OH, m = 46, methanol,
create new surface (e.g. by breaking a block of the solid in half ) requires CH3 OH, m = 32, acetone,
CH3 COCH3 , m = 58, n-dodecane,
work. One must overcome the attraction that binds the molecules of the CH3 (CH2 )10 CH3 , m = 28,
solid together. Thus the surface energy is of order of the binding energy cyclohexane, C6 H11 OH, m = 100,
of a molecule per area occupied by a molecule on the surface. The bind- n-hexadecane, CH3 (CH2 )14 CH3 , m = 28,
ing energy of a solid is generally much larger than the thermal energy T . benzene, C6 H6 , m = 78, carbon
Thermal fluctuations are minor and each atom is fixed at a position that tetrachloride, CCl4 , m = 153, polystyrene,
(CH2 CH(C6 H5 ))n , m = 52, polyvinyl
minimizes its potential energy (cf. the argument in Chapter 2 for viscosity). chloride, C2 H3 Cl, m = 62, ethanediol,
Thus the surface energy of a solid is often much more than T per surface HOCH2 CH2 OH, m = 62, formamide,
atom. Among solids, metals and ionic crystals have high surface energies HCONH2 , m = 45, glycerol,
(which moreover depends on the orientation of the surface with respect to HOCH2 CH(OH)CH2 OH, m = 31, water,
the crystal axes). Hydrocarbons have low surface energies. Rare-gas solids H2 O, m = 18, hydrogen peroxide,
H2 O2 , m = 34. For this wide range of
have the lowest ones liquids, the scaled surface energy remains
The interfacial energy shows up concretely whenever a liquid can change near 1.
its interfacial area. If a drop of liquid is suspended in space, it minimizes its
interfacial energy by adopting a spherical shape. If an external agent alters †For a polymer this would be a few back-
this shape, the work it does is α times the increase in area. If the fluid is dis- bone bonds and their attached atoms.
torted and then released, it begins to move so as to decrease its interfacial
area, converting interfacial energy into kinetic energy. The kinetic energy
is equal to α times the reduction in area (less any loss from viscous dissipa-
tion). If a spherical drop of radius R could reduce its volume by an amount
156 Interfaces
known as complete wetting. The liquid spreads to cover all the exposed area
of the solid surface (as much as the supply of liquid will allow). Complete
wetting tends to occur on high-energy solid surfaces, where αs0 is large. For
example, hydrocarbon oils typically wet metal surfaces completely. Often
a liquid is attracted to such a high-energy surface by enough to justify
complete wetting [6].
At the other extreme, it can happen that any contact between the solid
and the liquid costs energy. This occurs when the energy of an air–liquid
plus an air–solid interface is lower than that the of the solid–liquid interface:
α + αs0 < αs . Clearly this can only happen when αs > αs0 , so that the
surface repels the liquid. In this situation, a spherical drop in grazing contact
with the surface would separate from the surface (in the absence of gravity).
Such a surface is said to be completely nonwetting for the liquid.
For both complete wetting and complete nonwetting, α is smaller than
|αs − αs0 |. If instead α is larger than this difference, partial wetting occurs.
The liquid surface meets the solid surface at contact angle θ. Figure 6.8
shows how interfacial energy is related to forces in this familiar case. If
a section of the liquid of length L advances, a work −Lδxαs0 is done in
decreasing the bare surface. This is just the work that would be done if the
solid exerted a tensile force Lαs0 on the liquid trying to advance it. The
interfacial energy thus can be viewed as a force per unit length. This is why
interfacial energies are often called interfacial tensions. The liquid-covered
surface exerts an opposing tension Lαs . Finally, the advance would create
free liquid surface of amount Lδx cos(θ ), so that the free liquid exerts a
force Lα cos(θ ) opposing the advance. The sum of these forces must be
zero; otherwise the fluid would advance or retreat spontaneously. Thus
This is Young’s law of partial wetting [6]. Young’s law determines the con-
tact angle; for a circular drop the contact angle in turn determines the
radius R. The free surface of a drop must have a constant spherical curvature
everywhere; otherwise, the Laplace pressure would be different in different
places. When the solid neither attracts nor repels the liquid, αs = αs0 . This
is the case of neutral wetting. Here the contact angle θ is 90◦ for any value
of α. In general αs is different from αs0 . Then if α is made smaller, the Fig. 6.8
contact angle decreases to zero to produce complete wetting or increases to Left: a partially wetting droplet on a
180◦ to produce complete nonwetting. horizontal surface, viewed from the side.
The contact angle θ is indicated. The liquid
The surface tension force of a liquid against its supporting structures can surface meets the solid surface at the
be mechanically significant. Such capillary forces can make fluid rise to contact line. As θ approaches 0 the droplet
macroscopic heights in a thin tube. Fluid transport in plants is accomplished approaches the state of complete wetting.
largely by capillary forces. In fine granular materials with large surface to As θ approaches π , the droplet approaches
the state of complete nonwetting (center).
Right: motion of the contact line by an
amount dx along the surface stretches the
s0 liquid surface by an amount dx cos θ. The
s vectors show the forces on a unit length of
the contact line, including the bare solid
dx cos
surface tension αs0 , the liquid-covered
solid surface tension αs , and the liquid–air
dx surface tension α.
158 Interfaces
volume ratios, these capillary forces can distort their supporting structures
or bind them together. Thus small amounts of water in the air can readily
form wetting layers on fine particles, turning a free-flowing powder into a
solid cake.
6.2. Capillary solid A person standing on a cube of wet sand exerts a stress of
104 N/m2 . This stress must be supported by sand grains whose size is a.
(a) Roughly how much force does a given grain support for a given size a?
On some of the grains this force is a tensile force trying to separate two
adjacent grains. Treat these two grains as cubes of length a and suppose
that the cube faces are separated by a thin water film. Separating the cubes
increases the thickness of this film and reduces its cross-sectional area to
maintain constant volume.
(b) For what size a is the interfacial tension equal to the applied force? Grains
smaller than this ought to support the weight of the person.
This equation readily tells how R varies with time t, if we neglect the slowly
varying log(R/b):
R
[(V 3 α/η)t]1/10 . (6.2)
This is the well-known Tanner law of wetting [8]. It says that the spreading
slows in a universal way with time in the marginal case where S0
0.
Remarkably, Tanner’s law is observed even when the flat spreading pres-
sure S0 is large and should dominate the full spreading pressure. Instead
of advancing rapidly in response to the full spreading pressure, the macro-
scopic contact line advances gradually as though there were no flat spreading
pressure. As anticipated above, the extra spreading pressure is opposed by
an extra dissipation mechanism, not considered above, that is independent
of the angle θ. As explained by de Gennes and Hervet [9], the fluid spreads
far ahead of the macroscopic contact line in a thin precursor film or foot. The
larger the flat spreading pressure is, the faster this film advances. The flat
spreading pressure is balanced by dissipation in the precursor film. Only the
residual spreading pressure 12 αθ 2 remains to be balanced by the dissipation
in the macroscopic drop.
We notice that if the interface has no effect on the solute, f (z) is a constant
and the expression for F2 reduces to our initial expression F1 − T log().
162 Interfaces
driven to the surface, thus lowering the surface energy. Making the solute
more immiscible in the liquid, amounts to increasing the free energy cost
F1 defined above. But in fact, making the solute immiscible is generally not
a good strategy for lowering surface tension. Immiscible solute molecules
generally experience an effective attraction for each other owing to their
mutual repulsion from the solvent. Such attractions are present even when
the solute molecules are driven to the surface. Thus instead of spreading
over the surface and producing surface pressure, the solute molecules may
phase separate on the surface, thus cutting the desired surface pressure.
A way around this dilemma is to use amphiphilic solute molecules. Such
molecules have a part that by itself would be insoluble and another part
that by itself would be highly soluble. As we shall see in the next chapter,
such amphiphilic molecules can segregate strongly to the interface without
phase separating there.
ξu
Fig. 6.10
Chains adsorbed on a phantom surface.
Top: a single chain adsorbs in a thickness
ξu as described in the text. Middle: a chain
from the bulk solution approaches a
saturated adsorbed layer. Bottom: blobs of
the bulk chain exchange with those of the
surface, thereby binding the bulk chain and
weakening the attachment of the
pre-adsorbed chain.
pressure, tending to push it away from the surface. Since the diameter of our
segment is roughly z, the force that it feels is roughly ((z/2) − (z))z2 .
Now, the blobs at height z have a size that is also roughly z according to
our picture. Thus (z)
T /z3 . The force is thus of order (T /z3 )z2 or
T /z. This is of the same order as the tensile force we supposed above. Thus
under our assumption, the forces acting on any typical segment are of the
same order. Any adjustment required in order to bring the segments to an
equilibrium state of balanced forces is a minor adjustment that does not
alter the scaling behavior we have deduced.
This picture of the height concentration profile, originally postulated by
de Gennes [15], has been supported by several calculational approaches.
It has also been confirmed by realistic computer simulations [16]. Direct
experimental confirmation is difficult because the profile involves only
a tiny fraction of the adsorbed monomers. Indeed, the above argument
says that φ(z) ∼ z−4/3 , so that the total surface excess ∼ dzφ(z)
is dominated by the smallest z’s—i.e., z
ξu . Thus
ξu φ(ξu ), as
argued above. Still, scattering investigations of the profile show behavior
consistent with the de Gennes picture [17].
proximal power law can lead to measurable effects on the surface excess
and the concentration profile. The reasons for the proximal power are dif-
ficult to account for in simple terms, and direct evidence for it has proved
difficult to establish experimentally [19]. There is in particular no simple
connection between the proximal power and the fractal dimension D of a
self-avoiding polymer.
Our discussion of the density profile assumed that the chains had indefin-
itely long length. When their length is finite, the power-law regime is limited
to distances z smaller than the unperturbed size of the chains. In addition, the
finite-length chains appear to show a second power-law regime that implies
a falloff of concentration slightly weaker than the z−4/3 found above for
sufficiently large z, owing to tails of chains adsorbed at only one end [20].
approached, yielding a uniform shear rate and a uniform viscous stress. But
flow through an adsorbed layer need not have a uniform viscous stress.
Instead, force can be transmitted from the fluid to the polymers, and thence
directly to the surface. We saw in Chapter 4 that fluid velocity decays expo-
nentially in a semidilute solution of blob size ξ , with a decay length of the
order of ξ . Similar decay occurs in an adsorbing layer. We may picture the
layer as having an outermost sublayer of thickness R/2, extending from
height R/2 to height R and having blobs of size R/2. Then the adjacent
denser sublayer has thickness R/4 and blob size R/4, followed by suc-
cessively denser sublayers. We suppose that the flow velocity at the outer
boundary of the layer is v0 . Then on traversing the outer sublayer, the
velocity must decay to a finite fraction of v0 , say v0 /χ. If this χ is larger
than 2, the velocity decreases faster than the height z. After k sublayers,
z = R × 2−k and v = v0 × χ −k . Solving for −k yields
Fig. 6.11
Top: sketch of a saturated adsorbed
polymer layer in equilibrium with a
solution. Layer consists of loops at various
heights such that the distance between
loops at height z is comparable to z.
Speed v
Bottom: illustration of the procedure in the
text for accounting for the velocity profile.
Beyond the adsorbed layer, the velocity
v0
grows linearly with distance from the
surface. Within the adsorbed layer, the
velocity decreases by more than a factor 2
for each halving of the distance, resulting v0/2
in a power-law increase of velocity whose v0/χ
exponent is larger than 1. The velocity at
large distances extrapolates to a height t, t
the effective hydrodynamic thickness. It is Height z
comparable to the layer thickness R. R/4 R/2 R
Stress
Fig. 6.12
Flow of an entangled polymer solution past Velocity
a wall, sketched at top. Upper curve shows
the (uniform) shear stress. Lower broken
curve shows the velocity profile required to
attain the uniform stress. The extrapolation Height z
length b is shown on the horizontal axis. b
Here the stress is conventional fluid stress, caused by fluid and polymer
atoms hitting wall atoms. It can be characterized by the viscosity of the
solvent or of disconnected monomers. Though the means of transmitting
stress in these two regions are very different, the amount of stress must be
the same. The stress must pass from the bulk of the fluid to the container
via the wall layer. To attain the same stress as the bulk, the wall layer must
have a much larger velocity gradient, as sketched in Fig. 6.12. The velocity
extrapolates to zero beyond the wall, at a distance b called the slip length.
For simple fluids the slip length is of the order of an atomic size a. But if the
References 171
polymer viscosity is larger than the solvent viscosity by a large factor, the
slip length exceeds a by an equally large factor. Such large slip lengths are
readily observed, and they may reach the scale of microns [23]. Sometimes
one wants to avoid this slipping behavior. One way to achieve this is to
fasten polymers to the wall. Then the fluid polymers must disentangle from
the wall polymers in order to flow. One may envisage all kinds of exotic
behavior when these flows are fast enough to perturb the fastened poly-
mers or to inhibit re-entanglement. The engineering literature has no lack
of examples of such exotic behavior. Perhaps the best known are the “shark-
skin” patterns caused by stick-slip flow past a wall under these nonlinear
conditions [24].
6.5 Conclusion
In this chapter we have accounted for the energy scales associated with
liquid surfaces and interfaces. We have seen how these energies create geo-
metric and kinetic effects e.g. in wetting phenomena. Polymers and colloidal
suspensions interact strongly with surfaces and this interaction can alter the
interaction between two surfaces. It can also alter flow past a surface. With
these interfacial phenomena in mind, we are ready to explore the distinctive
phenomena caused by surface-loving molecules, or surfactants.
References
1. Tantec contact angle micrometer Tantek Co. 630 Estes Avenue, Schaumburg
IL 60193 http://www.tantecusa.com/camplus2.html (2002).
2. Pamphlet “Spinning Drop Tensiometer Site 04” KRUESS USA, Instruments
for Surface Chemistry 9305 Monroe Road, Suite B, Charlotte, NC 28270 -
1488 http://www.kruss-usa.com/public_pdf/Site04e.pdf
3. SE400: Sentech Discrete Wavelength Ellipsometer http://www.sentech.
com/ellipsometer.htm SENTECH Instruments GmbH Carl-Scheele-
Strae 16, 12489 Berlin, Germany.
4. F. Rondelez, D. Ausserre, and H. Hervet, Ann. Rev. of Physical Chemistry 38
317 (1987).
5. J. Israelachvili, Intermolecular and Surface Forces (New York: Academic,
1985) Table XVII, p. 158.
6. J. Israelachvili, op. cit. Chapter 14.
7. P. G. de Gennes, Rev. Mod. Phys. 57 827 (1985).
8. L. Tanner, J. Phys. D. 12 1473 (1979).
9. H. Hervet and P. G. de Gennes, Comptes. Rendus Acad. Sci. 299II 499 (1984).
10. L. Leger and J. F. Joanny, Rep. Prog. Phys. 55 431 (1992).
11. C. Redon, F. Brochard-Wyart, and F. Rondelez, Phys. Rev. Lett. 66 715 (1991).
12. R. D. Deegan, O. Bakajin, T. F. Dupont, G. Huber, S. R. Nagel, and T. A. Witten,
Nature 389 827 (1997); Phys. Rev. E. 62, 756–765 (2000).
13. S. H. Davis, Ann. Rev. Fluid Mech. 19 403 (1987).
14. D. Jasnow, Rep. Prog. Phys. 47 1059 (1984).
15. P. G. de Gennes, Adv. Colloid Interface Sci. 27 (1987).
16. J. de Joannis, R. K. Ballamudi, C. W. Park, J. Thomatos, and I. A. Bitsantis,
Europhysics Lett. 56 200 (2001).
17. L. Auvray and J. P. Cotton, Macromolecules 20 202 (1987).
18. E. Eisenriegler, J. Chem. Phys. 79 1052 (1983).
19. J. M. DiMeglio and C. Taupin, Macromolecules 22 2388 (1989).
172 Interfaces
7.1 Introduction
In the previous chapter we saw that the boundary of a liquid is respons-
ible for significant, controllable forces—capillary forces. In this chapter
we encounter an enormous expansion of that control through the use of
“amphiphilic” molecules. An amphiphilic molecule contains a part that by
itself would be soluble in the liquid and another part that would be insol-
uble in it. The insoluble parts have a more favorable free energy when they
are away from the liquid; thus these molecules tend to concentrate at the
liquid boundary. Because they influence surfaces and other boundaries,
these molecules are called surfactants.
In a sense a surfactant incorporates a bit of interface within itself. The
two dissimilar parts of the molecule, having opposite affinities for the
liquid, generally have a repulsive interaction with each other. This repul-
sion amounts to a free energy that was stored in each surfactant molecule
when it was made. The boundary between the two parts of the molecule
often plays the role of an interface. To minimize the repulsive energy
inherent in the amphiphilic molecules, they often organize to create real
interfaces throughout the liquid, as shown in Fig. 1.2. The resulting struc-
tures can dramatically alter the forces within a liquid. In this way capillary
forces come to control the bulk properties of the liquid, not just its surface
properties.
A second feature of surfactants is that they can cause immiscible liquids
to mix. Here one considers surfactant amphiphiles one of whose parts is
soluble only in the first liquid and the other of whose parts is soluble only
in the second liquid. Such surfactants are at equilibrium at the interface
between the two liquids. Adding more surfactants creates more interface.
With a sufficient amount of surfactant, this interface may grow so much that
any point in either fluid is only a microscopic distance from the interface.
Then the immiscible fluids are in effect mixed.
This chapter explores the special properties that amphiphilic molecules
impart to a liquid. As in previous chapters, we focus on the simplest level
of quantitative, predictive analysis. We identify the characteristic spatial
scales, energy scales, and time scales that make these fluids act the way they
do. We begin with a section on mixing; we recall the empirical rules that
describe the mixing of different molecules. The next section describes the
most common surfactant molecules, noting their resemblances and dif-
ferences. The following section examines how strong amphiphiles should
interact in a liquid, and why they should aggregate into micelles when the
174 Surfactants
in the same way. By this approach one may justify five principles that are
generally borne out empirically. (1) Positivity: the work required to trans-
fer our B molecule is generally positive and increases with the molecular
size. (2) Additivity: the work W may be found by adding contributions for
each atomic constituent of the B molecule being mixed. (3) Ordering: the
work is greater when the two liquids differ more in their static dielectric
constant and their optical index of refraction. (4) Reciprocity: if A and B
are of comparable size and if A is miscible in B, then B is miscible in A.
(5) Transitivity: if A, B, and C are of comparable size and A and B are
miscible in C, then A is miscible in B. We discuss each of these principles
briefly below. These important aspects of physical chemistry are discussed
much more carefully in Israelachvili [1] and Denbigh [2].
7.2.1 Positivity
In the discussion of van der Waals forces in Chapter 5 we found that two
identical polarizable B atoms have an attractive interaction in vacuum or in
any polarizable medium. Now, the van der Waals interaction perturbs each
atom only slightly; thus the net interaction of many atoms is nearly the
sum of pairwise interactions. This means that a B atom near a boundary
between B and A liquids feels an attraction towards the B liquid—i.e.,
the work W is positive. Now if a few atoms are joined to form A and B
molecules, the atomic polarizability of each atom is altered by its bonding
environment. Still, each atom of B feels a net attraction for its counterpart
on other B atoms even in an environment of A. Thus two B molecules feel
a net attraction, and again the work W is positive. The removal of our B
molecule to an A environment also affects the packing of the neighboring
molecules around our B molecule. But for the small organic molecules we
encounter in practice, these differences in packing are usually not dominant.
Furthermore, the group of atoms is itself a polarizable entity that must
feel a net attraction towards similar entities in an environment of different
polarizability. Our small B molecule will in general have a nonspherical
distribution of charge, though it has no net charge. At large distances, this
nonspherical distribution leads to a dipole field. This dipole may be oriented
by electric fields; in other words it leads to a molecular polarizability. This
polarizability creates van der Waals attraction just as atomic polarizability
does. Thus it adds another positive contribution to the work W .
7.2.2 Additivity
Many molecules, such as hydrocarbon chains, are too large to be considered
as small, fixed clusters of atoms. How can we understand the work W for
carrying a hydrocarbon chain B from a liquid of similar chains into some
other liquid A? As we have seen, this work results from the altered interac-
tions between each hydrocarbon group and its immediate surroundings in
the A medium compared with its original B surroundings. The interactions
are local. Thus it is reasonable to approximate W as the sum of contrib-
utions W0 that each hydrocarbon group would experience on being carried
from the B liquid to the A liquid. Adopting this point of view leads us to an
176 Surfactants
7.2.4 Reciprocity
For atomic liquids of the same atomic size, only the contrast in polarizability
leads to attraction and immiscibility. If A and B atoms differ in polarizability,
this limits the miscibility of A atoms in B and B atoms in A to a comparable
degree. The same is expected when A and B are small molecules. This
reciprocity becomes less valid as the contrast between A and B increases.
Mixing principles 177
7.2.5 Transitivity
The principle of transitivity says that two molecules that are miscible in a
third should be miscible in each other. This is another consequence of the
like-dissolves-like principle. We have argued that small-molecule solvents
that are similar in dielectric constant and in index of refraction should be
miscible in each other. If these properties are similar for solvents A and
C, and for B and C, they must be so for A and B. Thus A and B should
be miscible. If the molecules differ greatly in size, the reasoning breaks
down, and the transitivity property need not hold. The most obvious counter
example is the case where A and B are large hydrocarbon polymers and C is
a small hydrocarbon molecule.
with the un-ionized state of the native B medium. Thus the ionization gives
a negative contribution to W and promotes mixing.
These benefits of ionization do not come without a cost: namely, the
electrostatic energy of separating each ion from its oppositely charged part-
ner(s) in the native B medium. The electrostatic cost must be smaller than
about T . Otherwise, there would be a net positive cost to ionization even
considering the entropy effect above. To avoid this cost, the mutual elec-
trostatic energy of the two opposite ions must be made small. This energy
U is reduced in proportion to the reduction of the electric field of, say, the
negative B ion; it is smaller than in vacuum by a factor . Thus large
means small U and small electrostatic cost of separating. A useful measure
of this electrostatic energy is the Bjerrum length , defined in Chapter 3. It
is defined so that the mutual energy of two singly charged point ions in the
liquid at large separation r is given by U (r) = T (/r). Evidently varies
as 1/. For water at room temperature
0.7 nm. Accordingly, the energy
cost to separate two ions from r = 0.7 nm to a large separation is just the
thermal energy T .
In water most small singly charged ions dissociate and mix readily. The
energy gained by bringing them to a distance of, say, 1 nm from infinity is
only a fraction of T , in accordance with the formula above. At this distance,
the organized layers of polarized water molecules around each ion begin to
interfere with each other. Bringing the ions closer disrupts these “solvation
shells” and costs free energy. The result is that little further net energy would
be gained by bringing the ions into contact. The total energy of separation is
of the order of T , and the ions may readily dissociate. For example, sodium
and chloride ions dissolve in water at volume fractions up to about 10%.
For volume fractions smaller than this, the gain in translational entropy is
more than enough to compensate for any loss in electrostatic and solvation
energy.
Evidently, the special properties of water are important in allowing ions
to dissociate. In most solvents dissociation is a minor effect. Even in
water dissociation is difficult for ions bearing more than one charge. Two
divalent ions must be further than 2.8 nm apart in water in order to have an
electrostatic energy smaller than T .
CH3
O
O
CH3
CH3 K+ O S O
– O CH3
O
O O
O O
H3C O O CH3
+
N P – O
H3C CH O O
3
H3C
H3C O O O OH
O O
Fig. 7.1
Five common surfactant molecules from [4]. Top row: SDS, the cationic surfactant cetyl trimethyl ammonium bromide (CTAB). Bottom,
the phospholipid 1-palmitoyl-2-oleoylphosphatidylcholine (POPC), sodium bis(2-ethylhexyl)sulfosuccinate (AOT), pentaethylene glycol
monododecyl ether (C12 E5).
sodium dodecyl sulfate or SDS, is much less miscible than its miscible part.
Figure 7.1 shows a diagram of it.
SDS, also called sodium lauryl sulfate, is the leading ingredient in house-
hold cleaning products, like soap, detergent, and shampoo. Many other
surfactants are variations on this theme. The hydrocarbon part, called the
tail, may be lengthened by four or more carbons. The polar sulfate head
group SO− −
4 can be replaced by carboxylate, CO3 , or by phosphate, PO4 .
−
Finally, some head groups may have two hydrocarbon tails attached. As we
shall see below, the significant differences in these molecules lie not only
in the solubility of the different parts but also in the relative bulk of the head
and tails and in the extensibility of the tails.
The polar head is most commonly a negative ion like the ones above.
Then the molecule is called an anionic surfactant. Sometimes one chooses
180 Surfactants
a cationic group such as ammonium, NH+ 4 as the polar head as in the CTAB
molecule shown in Fig. 7.1. Naturally, cationic and anionic surfactants
have a strong mutual electrostatic attraction. The dissociating ion is most
commonly a singly charged metal ion like sodium. But this ion can be an
organic group as well. Likewise, the ions may be weakly dissociating. An
important example is an OH− H+ pair whose dissociation may be controlled
by the ambient H+ concentration or pH. Amine groups NH− are a second
important weakly dissociating group.
Some surfactants, called nonionic, have polar heads without dissoci-
ating ions. One example is a zwitter ion, in which an organic cation is
attached to an organic anion by e.g., a short hydrocarbon chain. A second
important example is the ethylene glycol group, whose polar oxygens asso-
ciate strongly with surrounding water molecules and lead to a negative
contribution to the work of mixing W .
To summarize, surfactants vary in strength from weak to strong. Weak
surfactants have parts with little difference in miscibility. Strong surfact-
ants have parts that differ greatly. Beyond this difference, surfactants
differ in the relative bulk and deformability of their polar and nonpo-
lar parts. Further, the polar parts differ in being anionic, cationic, or
nonionic.
Fig. 7.2
Three types of micelles. Left to right,
spherical, cylindrical, and bilayer.
KB ←→ BK .
For large aggregation number K, [B] grows very slowly with C—[B] is
practically constant.
What determines the aggregation number K? Evidently it has something
to do with the shape and deformability of the surfactant molecules. For
example, an aggregate of SDS with K = 2 would not be very favorable,
182 Surfactants
since the two tails could not avoid contact with the solvent. Thus nearly
as much work would be required to put the aggregate into the solvent as
to put the dispersed molecules there. But a spherical micelle with tens of
molecules can readily shield the tails from the solvent. Beyond this point
there is little to be gained by increasing the aggregation number. In addition,
there is something to lose: the charged head groups repel each other, and
enlarging the aggregation number clearly increases this repulsion.
7.1. Law of mass action Molecules A, and B in a dilutesolution can form a com-
bined species AB. The work required to separate an AB into its constituents
A and B at two given positions far apart is W . Show that in equilibrium the
concentrations [A], [B], and [AB] are related by
[K + 1]/[K] = [1]/C0 .
In terms of the concentration of micelles Cp ≡ C/K , this says K → Our scaling law for the sums shows that
Cp /C0 . Thus the micelle size can be controlled and predicted over a wide Mmn
q(C)n−m . That is, all moments
range of concentration based on fundamental principles. behave as though they had dimensions of
This prediction of how concentration affects micelle size can readily be q(C)n−m . The averages discussed in the text
tested by scattering experiments. For example, are moments of this form: K = M0,1
static scattering
measures
and K w = M1,2 . They both grow with
the weight-averaged micelle size K w ≡ K 2 [K]/( K[K]). All such the same power of q(C), so that their ratio
averages scale with the same 12 power of overall concentration C † . Thus remains finite as either goes to infinity.
184 Surfactants
one may verify the square root law by measuring scattering intensity versus
concentration.
If the limiting monomer concentration C0 is small, we can create solu-
tions with very long worms. As these grow longer, they curve more and
more, eventually taking the form of random walks like the polymers stud-
ied in Chapters 3 and 4. In our study of these polymers we learned that
the interaction of the chain with itself affects its size in a qualitative way.
Repulsive interactions lead to the fractal properties of self-avoiding walks
rather than those of simple random walks. Repulsive interactions between
parts of a wormlike micelle must ultimately have the same effect. In par-
ticular, the work W needed to add a monomer chain now depends on the
chain size K. Expressed in terms of the equilibrium constant for K-mers
assembled from monomers, the K dependence has the form
[K] 1
= K−1 K γ −1 (7.4)
[1]K C̃0
The K γ −1 factor arises from the self repulsion, and the constant C̃0 plays
the role of C0 . Here γ is an universal exponent like the fractal dimension
D characteristic of any self-repelling polymer [6]. The value of γ is about
1.17 [6]. This factor alters the dependence of K on concentration, as we
may deduce by repeating the reasoning above. We first note that K ∼ 1/δ
as before. However, the overall concentration C has a different scaling. For
large K we may follow the scaling argument of the footnote to infer
C → constant δ −1−γ
molecules is great and the energy barrier is high. Thus the micelle relaxa-
tion time—denoted τ2 —is many times longer than the surfactant relaxation
time τ1 . For SDS in typical conditions it is on the scale of milliseconds.
It may be much longer in very pure solution where there are few impurity
molecules to aid the nucleation of micelles [5].
These two characteristic times tell only part of the story of the dynamics
of micelles. In the rheology section below we shall encounter further forms
of relaxation that become important for interacting micelles.
The micelles on the right are subject to the effects of repulsion and thus
obey Eq. (7.6)
[g + g]
C̃0−1 2γ −1 g 1−γ .
[g][g]
The repulsion suppresses the connection of the two pieces by a factor g 1−γ .
Micelle interaction 187
Here x is defined as the contents of the large (. . .). The repulsion has added
a new feature: the g-dependent prefactor. Again [K] falls off exponentially
with the controlling parameter x, as it did without repulsion. The only
change is to replace C0 by C̃0 g γ −1 . This [K] varies smoothly with x, so
it is approximately valid even when K is not an integer multiple of g.
Repeating the reasoning after Eq. (7.2), we find
The self-interaction has changed the growth exponent from 0.5 to 0.6. This
subtle change in how micelles grow with concentration will prove important
when we consider the viscosity of these solutions below.
7.2. Worm interpenetration As noted above, the growth of wormlike micelles with
concentration leads to interpenetration.
(a) How does the overlap ratio φ/φ ∗ increase with concentration, assuming
(i) self-avoidance effects are negligible so that D = 2, and (ii) self-
avoidance effects are important, so that D = 53 ?
(b) Suppose that the cylindrical micelles have a radius of 2 nm and that each
surfactant occupies a volume of 0.4 nm3 . Suppose that the persistence
length is 20 nm, so that the persistent segments are long and slender. For
significant entanglement effects to occur, one should have φ/φ ∗ > 10.
But if the volume fraction of worms exceeds 10–20%, their packing leads
to nematic order as treated in Problem 5.10. Then the resemblance to a
polymer liquid diminishes. If nematic order to be avoided, φ < 0.1, so
that φ ∗ must be less than 0.01. Find the condition on the characteristic
concentrations C0 or C̃0 that permits this much entanglement. What is the
corresponding volume fraction of free monomers?
Wormlike micelles often have substantial rigidity, as noted above: they
remain straight over a length p that is several times their width. Thus at
sufficiently high concentration, the micelles may align to form a nematic
liquid crystal of the kind discussed at the end of Chapter 5. To understand
what concentration is necessary, it is useful to think of the micelles as being
188 Surfactants
(a) (b) (c)
Fig. 7.3
Sketch of three periodic bicontinuous
micellar structures discovered in
concentrated surfactant solutions, courtesy
John Seddon [12] with permission from
Elsevier. (c) is known as the plumber’s
nightmare [11]. a a a
cut up into straight rods of length p . Then one may adopt the reasoning
used in Chapter 5 section on organized states.
Surfactants whose shape favors branch junctions can at sufficiently high
concentration form structures that are branched everywhere. Such struc-
tures resemble a regular lattice of connected cylinders. The simplest such
structure is a simple cubic lattice resembling a tangle of pipes, known as
the plumber’s nightmare, as shown in Fig. 7.3. For example, 60 volume
percent glycerol monoolein (a C18 hydrocarbon chain with polar OH and
C== O groups at one end) in water is believed to form this phase [10]. Many
similar periodic network structures have been discovered. These are called
Luzzati phases [11].
Here κ, κ̄, and c0 are material constants of the interface. The constant κ is
called the bending modulus or bending stiffness. The κ energy is minimal
when the mean curvature has the value c0 ; this c0 is known as the spon-
taneous curvature. The term in κ̄ is not relevant for our purposes. Thermal
fluctuations do not affect this energy. Thereason is a geometric fact called
the Gauss–Bonnet Theorem. It says that dsC1 C2 over any smooth sur- R⬘
face with a fixed boundary is unaffected by continuous deformation of the S
surface. That is, the integral is a topological invariant. If the deformation R
increases the Gaussian curvature in one place, it must decrease it in another
so as to compensate [12]. Thus compensation is illustrated in Fig. 7.5 and
discussed in the Appendix.
Knowing the energetic cost of bending, we can now analyze the fluc-
tuating shapes that occur in equilibrium. First we consider sections of the
surface that are small enough to be nearly flat. This must be true for small Fig. 7.5
enough sections; otherwise we would not be able to characterize our micelle Bump construction on a flat surface used
as a smooth, two-dimensional surface. For such a nearly flat surface, there to demonstrate the Gauss–Bonnet theorem.
is a convenient way to describe the fluctuating shapes. We represent the A flat surface was deformed to make the
circular bump shown, with positive
surface by giving its local height h above some reference plane (x, y), that Gaussian curvature at the middle. In order
makes only a small angle with the surface everywhere. The smallness of to rejoin the flat surface smoothly, the
these angles means that ∂h/∂x and ∂h/∂y are much smaller than unity. This surface must have negative Gaussian
h(x, y) scheme is known as the Monge representation. curvature at the rim. The Appendix shows
The energy of Eq. (7.9) can be readily expressed in the Monge rep- that the Gaussian curvature averages
to zero.
resentation. Indeed, the mean curvature is simply 12 ∇ 2 h, so that E =
1
2 2†
2 dx dy κ|∇ h| . Noting that the energy density is invariant under trans-
†This is the energy we encountered in
lation, we may simplify by expressing h as a sum of plane waves. We define Problem 5.4, with ψ corresponding to h.
h̃(q) via
h(x, y) = L−2 h̃(q)eiq·r . (7.10)
q
The wave vectors q are a complete set of waves needed to represent functions
on this square region. Then
1 −2 2 2
E= κL (q ) |h̃q |2 . (7.11)
2 q
190 Surfactants
1 2
T
h2 L = L−4 |h̃q |2 = L−2 .
q
κ q
q2
T −2 −2
μ2 = L−4 q 2 |h̃(q)|2 = L q .
q
κ q
2
T L qmax d 2q T qmax
μ = L−2
2
log . (7.13)
κ 2π qmin q2 κ qmin
We recall that qmin is of order L−1 . As for qmax , 1/qmax is the smal-
lest wavelength at which significant fluctuations occur. This smallest
wavelength is roughly the thickness b of the membrane. Waves with length
shorter than the thickness are not adequately treated by the thin-membrane
formula (7.9). Such waves cost more energy than this formula suggests,
and their effect can be neglected. We conclude that the typical local slope
of a patch of surface of size L is of order (T /κ log(L/b))1/2 —it increases
logarithmically with the patch size. If the patch size is small, the slope is
smaller than unity. (Otherwise there would be no scale on which the surface
was smooth.) As one enlarges the patch, there comes point, L ≡ p where
Micelle interaction 191
the typical slope becomes unity. For patches larger than this, our Monge
approximation breaks down. Still, we expect the slopes to continue to grow
with L. This means that for patches much larger than p , two arbitrarily
chosen points on the membrane will typically have completely independ-
ent orientations. Thus the orientations become uncorrelated over distances
longer than p . Thus p represents a persistence length like the persistence
length encountered above for wormlike micelles. The orientation at a given
point persists out to distances of order p . Since 1 = T /κ log(p /b), we
conclude
p
b exp(κ/T ). (7.14)
7.3. Size variability of micelles In our example we supposed that the solvent
containing micelles all had the same radius R. To examine this supposition,
consider the energetic cost of changing this radius. Assume that the micelle
radius R = 1/c0 . For a given bending stiffness κ, how much may R decrease
from 1/c0 before the bending energy Eκ increases by T ?
Our surfactant interface has reduced its curvature towards the spontan-
eous curvature c0 . This has reduced the energy of our system. Eκ (V ) =
2κ(/(3V ) − c0 )2 . As long as the micelles have a radius much smal-
ler than 1/c0 , this energy decreases as the amount V of included solvent
increases:
Eκ (V )
2κ 3 /(9V 2 ).
This makes a negative contribution to the work required. But there are other
contributions to consider. One such contribution is the interaction energy
of the transferred B molecules. Most of the B molecules transferred remain
surrounded by B, so that they feel no change in interaction energy. Only
the ones adjacent to the surfactant layer experience a change; the associated
work is simply proportional to the area , and is independent of V . We
shall return to this contribution below. There is a further form of energy
to be noted. Adding the B solvent reduces the number of micelles. This
number Nm is evidently /(4π R 2 ), or 3 /(36π V 2 ). This reduction in the † In more detail, each small increment of
number of micelles gives the solution fewer degrees of freedom. As we saw micelle number δNm costs a work δW =
T δNm log(/Nm ), where is the volume
in the section on mixing principles, this requires a work T Nm log(v)+const of the system in units of the micelle volume.
where v is the volume fraction per micelle† . In terms of the solvent volume To find the total work, one adds (integrates)
V , this amounts to a free energy cost ET : the work δW for each increment of micelle
number. The resulting integral simplifies to
T Nm log(v)+const in the regime where the
ET
(T 3 /36π V −2 ) (log v + const). volume fraction of micelles is small.
ways as the large colloidal particles discussed in Chapter 5. The larger they
become, the larger their attractive interaction. Beyond some given size the
attraction leads the micelles to precipitate into their own micelle-rich phase
rather than dispersing [15].
In our example we assumed that κ was large enough that the thermal
fluctuations of the surfactant interface would be small. This meant that the
micelles had to be smaller than the persistence length p . If R becomes
larger than p then the micelles fluctuate strongly in shape. If the distance
between micelles is comparable to R, the micelles can even merge to make
larger irregular objects. A limiting case of these fluctuations is the situ-
lρ ation where c0 = 0, the volumes of A and B solvent are equal, and the
radius R
p . Then the distinction between the inside and the outside
Fig. 7.8 of the droplets disappears. The microemulsion can then be viewed as a
Left: Schematic sketch of a symmetric, collection of cells of size p , each containing randomly A or B solvent,
bicontinuous microemulsion. Right: with a surfactant membrane between cells of different solvent, as shown
cellular model of the microemulsion in Fig. 7.8. Such bicontinuous microemulsions were first discussed in these
according to DeGennes and Taupin [13].
terms by DeGennes and Taupin [13]. Since that time, the so-called Winsor
The persistence length p is shown.
III microemulsions [16] known to chemical engineers were identified as
† Awell-known Winsor microemulsion con-
bicontinuous microemulsions [17]† (cf. Fig. 1.2).
sists of SDS, brine, and toluene, with a small The larger the typical size of the microemulsion droplets, the smaller
amount of butanol. The SDS has a concen-
the surfactant fraction becomes. One may ask how small the fraction of
tration of 6–10 g/l, the butanol has a weight
fraction of 3%, the brine contains 6.5 weight surfactant can become. Is there a limit to how little surfactant is needed
percent salt, and the toluene and water are in order to mix oil and water in an equilibrium dispersion? Many Winsor
in equal volumes within a factor 2. microemulsions achieve a surfactant volume fraction of as low as 10−4 or
10−5 . The typical size R̂ of the oil or water volumes appears in scattering
experiments as a peak at some q ∗
1/R̃, as shown in Fig. 7.10. Such
measurements indicate typical sizes of the order of 100 nm [15]. But it
is difficult in practice to increase R̂ much beyond this scale, or to reduce
the amount of surfactant below this level. To increase R̂ further requires
that the bending stiffness κ be increased so that the persistence length p
increases. But in addition, the spontaneous curvature c0 must decrease so
that c0 1/R̂. Thus 1/c0 must become much larger than the surfactant
layer thickness b. However, slight changes in the A or B solvents can make
small changes in c0 and thus reduce R̂ significantly.
There is a great incentive to create stable microemulsions in order to
mix immiscible fluids. Accordingly, the structure and phase stability of
microemulsions is extensively studied. The simple scheme of random
water and oil volumes sketched above has given way to much more
systematic and powerful analyses [19]. A major focus of modern the-
ories is to explain the occurrence and the coexistence of phases in the
oil–water–surfactant systems. Figure 7.9 shows a typical phase diagram,
indicating the multiple phases that can occur and the involuted boundaries
separating them.
All of this discussion presupposes that the solvent B enters the bare
micelles. In order to enter, the first B molecules must create a surfactant-B
interface and pay any associated interfacial energy α, as noted above. If
this energy price is higher than the gain in curvature energy discussed above,
the energy to absorb the next B is unfavorable. For a favorable balance,
Mixing immiscible liquids: microemulsions 195
105
= 6 wt%
104
I(q) [cm–1]
12
103
21
Fig. 7.10
Neutron scattering intensity for 102
bicontinuous microemulsions, reproduced
from [18, figure 17], with permission from
Elsevier. The microemulsion consists of
equal volumes of D2 O and n-octane. With 101
the indicated weight percent of C12 E5
surfactant, and closely resembles the
microemulsion of Fig. 1.2. The wave
vector q at the peak decreases as the 100
0 0.2 0.4 0.6 0.8
amount of surfactant decreases, indicating
larger droplets of oil and water. q [nm–1]
Fig. 7.11
Interfacial tension versus temperature in a
bicontinuous microemulsion, from [18, 32
figure 5b], with permission from Elsevier.
Solution contains three phases, with 30
overall composition of 58.5% water,
cb
28 Tu
36.5% n-decane, and 5% C8 E5 nonionic
T (°C)
is thus of order κ/R 2 . We recall that typical interfacial tensions are of the
order T per interfacial degree of freedom, or several T per square nanometer.
Our micelles have a diameter of many nanometers. Thus the interfacial ten-
sion can be orders of magnitude smaller than that of a simple-liquid interface
(cf. Fig. 7.11).
Mixing immiscible liquids: microemulsions 197
are dispersed in the liquid. Then when the pressure is released, the gas may escape
by expanding these preexisting bubbles: there is no barrier to overcome. Thus the
release of gas is much faster. This project investigates this effect.
A 1 L can of 10 cm diameter is shaken with an amplitude of 10 cm at 10 times
per second. You can assume that the center of the liquid remains stationary while
the can shakes.
(a) Estimate the shear rate γ̇ inside.
(b) What is the order of magnitude of the viscous stress in the liquid, assuming it
has the viscosity of water?
The shearing force pushes the air at the top of the can into the fluid and breaks it
up into small bubbles. The bubble size is such that the interfacial tension pushing
bubbles towards a spherical shape is comparable to the viscous force tending to
elongate the bubbles.
(c) For the given shear rate, what is the size of bubble for which these two forces
are comparable? You may assume that the interfacial tension is the same as that
between water and air.
(d) Supposing that the 5% of the original 1 L volume that was air is dispersed into
these small bubbles, what is the distance between these bubbles compared to
the average distance to the wall of the can?
(e) Suppose the bubbles rise as though they were Stokes spheres as discussed in
Chapter 4. (Bubbles actually rise somewhat faster than solid spheres of the
same size. Also a cloud of bubbles rises faster than an individual bubble.) How
long will they stay in the liquid? This is the time you must wait before opening
the can in order to avoid a mess.
the same in every way except that they are mutually immiscible. The chain
contains N monomers, with N /2 in each subchain or block. We label the
two blocks A and B. We first recall the criterion for being immiscible. We
denote the work to insert one isolated B monomer into a liquid of pure A
polymers (or the reverse) as χ T † . Most A and B monomers encountered in † In nonsymmetric polymers, χ T as con-
practice have χ parameters that are only a few times 10−2 or smaller. In our ventionally defined [26] is the average of
discussion below, we will make the simplifying assumption that χ 1. the work of inserting a B into A and the work
of inserting A into B. In our symmetric case
The χ parameter measures the immiscibility. As noted at the beginning these two works are the same.
of this chapter, two liquids remain miscible as long as the work to insert
remains below about T . The work to insert a random coil B block into pure
A is of order N χ T . Unless N is sufficiently large, this work is negligible,
and free A and B chains mix readily. However, as N increases, the work
increases; ultimately, phase separation occurs. The same is true when A and
B chains are attached to make copolymers. If χ N 1, the A and B blocks
must phase separate, so that virtually all the A monomers are surrounded by
A blocks and likewise for the B blocks. In the limit of large N any diblock
polymer liquid with a positive χ must become strongly phase separated in
this sense.
This strong phase separation imposes a serious constraint on the indi-
vidual chains. Each junction point between the two blocks is effectively
confined to the interface between A and B regions. To see why, we imagine
moving a junction point from the interface into an adjacent B region, for
example. This movement necessarily brings monomers from the A block
of that polymer into the B region. Each monomer thus pulled into the B
region requires a work of order χ T . Since χ N 1, even a small fraction
of the A block thus pulled requires a work exceeding T . This work must
decrease the probability that the junction point moves spontaneously into
the B region according to the Boltzmann Principle. The larger N is, the less
free the junction point is to explore the B region. Below, we shall treat the
junction points as being confined near the A–B interfaces.
Our liquid consists of regions of nearly pure A or pure B. Since this
is polymer melt, each point in an A region is uniformly occupied by A
monomers†† . The A monomers encountered at this point have a pecu- †† The cost of altering this uniform dens-
liar property not shared with ordinary polymer melts. Each such monomer ity, by e.g., mixing empty space into the A
region, can be viewed as a mixing process
belongs to a block that ends at the A–B interface. The same is true for the B
with its own χ parameter. This cost is pre-
monomers. This constraint limits the size of the A or B regions. If a region is sumed arbitrarily larger than the χ for AB
large, then blocks passing through the center must necessarily be elongated mixing. Thus we can neglect the effect of
in order to extend to the A–B interface as required. Such elongation costs A–B phase separation on the individual A
free energy and this cost makes large regions unfavorable. and B densities.
readily estimate how the free energy of the liquid depends on h for chains
of a given length. It is most convenient to measure chain length in terms of
the volume V that a chain displaces in the melt. Evidently V is proportional
to the number of monomers N . A free chain in the melt state is a random
walk, whose mean-squared end-to-end distance R 2 is also proportional to
V : R 2 = V /a. The proportionality constant a is a characteristic length
that depends on the chain structure and the packing of the monomers in the
melt state, but it is independent of the chain length. We call it the packing
length.
If we change the size h, one major effect is to alter the interfacial energy.
As with any liquid interface, the energy is proportional to the area of the
interface and is characterized by an interfacial energy α. The ultimate origin
of this energy is the immiscibility of the A and B blocks, measured by χ. For
our purposes, the connection between χ and α is unimportant. We need only
observe that α arises because of a thin mixing region that is independent of
† The mixing region consists of chain seg- chain length† . Thus α is independent of chain length for long enough chains.
ments of length k such that χ k
1. This The interfacial energy per unit volume is the interfacial area per unit volume
mixing region is the same as that at the
times α. The chains for a given area A of A–B interface occupy a volume
interface between A and B homopolymers.
2Ah. Thus the interfacial energy per unit volume is αA/(2Ah) = 12 α/h.
The interfacial energy per chain, denoted U , is thus given by U = 12 V α/h.
Evidently this U is smallest when h is largest. The interfacial energy favors
large phase-separated regions.
As noted above, there is an opposing cost to increase h. This is the cost of
the required chain elongation. If the slab thickness is h, chains must stretch
to lengths of order h. Indeed, any monomer at the center of a slab must
belong to a chain whose end-to-end distance is atleast h. Any monomer
picked at random lies at a distance of order h from the nearest interface.
Thus to increase h is to increase the end-to-end distances of the chains in
proportion to h. The free energy S required for this stretching is a basic
property of random-walk molecules. As discussed in Chapter 3, the work
needed to stretch a chain to an end-to-end distance h perpendicular to the
slab is given by
1 h2 ah2
S(h) = T 2
T .
2 R /3 V
Since most of the chains extend to a distance of order h, the chains have an
average stretching energy S of order S(h).
If the system is arranged into slabs of a given size h, the energy per chain
is S + U (h). If the system is free to adjust h, the relative probability of a
given h is e−(S +U (h))/T . The most likely value h∗ is that which minimizes
S + U (h). Evidently
∂
0= [S + U (h)]h∗ .
∂h
Since S varies as h2 while U varies as h−1 , we have
2S − U (h)
0= ∗. (7.16)
h h
Amphiphilic polymers 201
This equation tells us three things. First it gives the scaling of the preferred
size h∗ with chain size V
so that
h∗ ∼ V 2/3 (α/aT )1/3 . (7.17)
We note that h∗ becomes arbitrarily larger than the unperturbed size
R(∼ V 1/2 ) in the long-chain limit. This confirms that the main distortion
of the blocks resulting from the phase separation is indeed one of stretching,
as implicitly assumed above.
Second Eq. (7.16) tells us how the free energy per chain F ∗ grows with
the chain length V :
Since h∗ grows arbitrarily bigger than the resting size R, this energy grows
indefinitely relative to T . Thus, the stable structure contains significant
energy. Even when the stretching energy is little more than T per chain,
the total stretching energy is large, because there are many chains in the
structural unit. The number of chains in a cube of size h∗ is (h∗ )3 /V ∼
V α/(aT ). Both α and a have length dimensions of atomic scale. But the
chain volume V is typically thousands of times the volume of an atom. Thus
even when the immiscibility is weak, the number of chains in a region of
size h∗ is typically hundreds, and the energy associated with the structure
is many times T .
This large energy means that the region size h is well defined and varies
little from h∗ . If h departs slightly from h∗ , the free energy per chain can be
found by expanding S + U (h) around its minimum. The free energy in a
volume (h∗ )3 thus has the form (h∗ )3 V F ∗ [1+(const)(h−h∗ )/h∗ ]2 , where
(const) is of order unity. Since F ∗ is large, as argued above, even a small
departure of h from h∗ costs more than T and is thus improbable. The size
of the regions is very well defined and becomes more so with longer chain
length. Likewise, significant distortions of the shape of the A or B regions
creates stretching energy of order S and is thus strongly suppressed.
Finally the equilibrium condition for h∗ of Eq. (7.16) gives the ratio of
stretching energy S to interfacial energy U at equilibrium. Though we do
not know the coefficient relating S to h2 , we do know that the power is 2.
Thus we know rigorously that dS /dh = 2S /h. U also varies as a known
power of h, and its derivative is rigorously given by dU /dh = −U /h. Thus † Such “virial relations” appear whenever
Eq. (7.16) does not have unknown coefficients. When it is satisfied, we must the free energy of a system is the sum of
have 2S = U . The interfacial energy must be twice the stretching energy two parts, each of which varies as a power
at equilibrium. In other terms, the stretching energy is 13 of the total energy† . of some variational parameter.
the A–B interface as before. As solvent enters the B regions, their volume
increases, and the B blocks expand and stretch to dissolve in the solvent as
much as possible. The result is strong stretching of the B blocks and great
asymmetry between the A and B regions. The calculation of the region size
h is modified, because the volume per chain is no longer fixed and the
quality of the solvent must be considered. These changes, though, prove
minor so that the same mathematical methods used for diblock melts can
be modified to describe these solutions. Now in addition to the structures
above, one can form free spherical or wormlike micelles in the solution.
One can also attach chains to solid surfaces to form the brush structures
discussed in Chapter 5. That chapter gives a hint of how the concentration
profile of brushes or other micelles arises.
One may wonder how well amphiphilic polymers function as surfactants,
dispersing immiscible solvents and reducing interfacial tension. In practice,
these polymers are little used for these functions. Part of the reason lies in
the large energy of these stretched polymer blocks. Many of the functions
of surfactants require their micelles to be flexible and deformable. For
example, the Helfrich repulsion between surfactant bilayers relies on this
flexibility. But polymer micelles have little of this flexibility, because their
deformation energy is large, as we saw above. For example, the bending
modulus of a lamellar layer is of the order of the energy within a cube
of thickness h as thick as the layer. As shown in this section, this energy
grows rapidly with polymer molecular weight. Thus the main interest in
amphiphilic polymers lies in the great range of ordered and controllable
structures that they make.
Fig. 7.12
101
Elastic modulus for emulsions made of
spherical droplets of uniform size 100
from [28], © American Institute of
Physics, courtesy T. Mason. Droplets of 10–1
Scaled modulus
silicone oil coated with SDS were
dispersed in water. Droplets had radii r of 10–2
0.25 μm (solid circles) 0.37 μm
10–3
(triangles), 0.53 μm (squares), and 0.74
μm (diamonds). The horizontal axis is the 10–4
volume fraction of the droplets, including
the volume displaced by the polar heads. 10–5
The vertical axis is the measured elastic
modulus divided α/r. The large circles 10–6
0.5 0.6 0.7 0.8 0.9 1
report measurements of osmotic pressure
for r
0.48, scaled in the same way. The Effective volume fraction
modulus and osmotic pressure are very
similar in magnitude for all but the largest
volume fractions. They rise sharply as the The relaxation of stress in these packed micelle or droplet structures
volume fraction increases so that the
droplets push against each other.
can be very slow. When the droplets are concentrated enough to touch, an
applied stress can only relax if the droplets can move past each other. To
do this, a droplet must increase its area by a nonzero fraction. This motion
thus requires an energy of order αR 2 , which can be much larger than T .
Thus the spontaneous relaxation of stress requires activated motion as is
exponentially slow. For many purposes the system can be regarded as a solid.
or
dσ/dt
−σ (t)(τb /τrep )1/2 /τb
−σ (t)(τb τrep )−1/2 . (7.20)
This suggests that the stress decays exponentially with a relaxation time
τ
(τb τrep )1/2 . Remarkably, the relaxation time is both much shorter than
the reptation time and much longer than the breaking time.
It is interesting to see how τ increases with the mean length L̄, with
everything else held constant, using τb ∼ 1/L̄ and τrep ∼ L̄3 , we find
τ ∼ L̄. It increases much slower with length than the reptation time does.
This kind of breaking-dominated relaxation has been confirmed in rhe-
ological experiments [7]. Its most important effect is on the viscosity η of
206 Surfactants
the fluid. As for any fluid, this viscosity is roughly the product of the step-
strain modulus G0 and the stress relaxation time. The step-strain modulus,
like the osmotic pressure is roughly T /ξ 3 , where ξ is the distance between
chains in the semidilute solution. The relaxation time is the τ discussed
above. Knowing the dependence of L, and ξ on concentration φ, we can
readily predict how the viscosity depends on concentration. Self-repulsion
effects appear important for these solutions [7], so we must use the scaling
exponents D and γ introduced above. For wormlike micelles, L̄ ∼ φ y ,
where y is given by Eq. (7.8). We have deduced the dependence of τrep on
ξ in Eq. (7.17): τrep ∼ L̄3 ξ 3−3D . Thus for ordinary reptating polymers in
these conditions,
The micelles have virtually the same growth of viscosity with concentration
as a solution of ordinary polymers. The competing effects of increasing L̄
and the breaking kinetics nearly compensate.
As noted above, wormlike micelles can sometimes branch or join together
to form a network. The stress relaxation of branched polymers is a tricky
subject that depends on the state of branching, as noted in Chapter 4. The
crosslinks between micelles can inhibit stress relaxation altogether. We shall
not attempt to give a general account of how stress relaxes in these cases.
If the imposed shear rate is significant, the solution may be altered by the
shear. Some types of alteration were discussed in Chapter 4. The chief effect
is mutual alignment of polymers that reduces the degree of entanglement.
In a micellar solution, new forms of alteration appear [29]. The length
distribution of the worms can change, and can be different for worms with
different orientations. The breaking rates can also increase.
(b)
Fig. 7.13
The three principal orientations of a
(c) lamellar liquid under shear.
(a) Layer-sliding orientation,
(b) Layer-shearing orientation,
(c) Layer-stretching orientation. Shear
strain is in the plane of the paper in all
three pictures. The layer normal is
indicated by an arrow. In (a) a fluctuation
in two of the layers is shown to illustrate
the resulting compression under shear.
an energy cost even with no motion; the bending cost is controlled by the
bending stiffness κ; the compressional cost is controlled by the Helfrich
repulsion or by Coulomb repulsion. Here steady flow is impossible without
breaking or strongly distorting the layers.
If the shear is applied in a generic direction, one expects the lamellae
to reorient so as to allow easier flow. That is, the lamellae should adopt
orientation (a), and the viscosity should be little greater than that of the
solvent. This behavior is observed sometimes when the lamellae are care-
fully prealigned [30]. But typically a lamellar solution is disordered, with
different orientations in different places. Such solutions do not flow readily
at all. It is not merely that they show a high viscosity. They have a dynamic
modulus that varies roughly as the square root of frequency. This contrasts
with the behavior of a viscous liquid, where the dynamic modulus is pro-
portional to frequency. This square-root behavior may persist down to the
lowest frequencies measurable, implying a stress qualitatively stronger than
that of a liquid.
The reasons for this square-root behavior are ill understood. It is believed
that the lamellae are forced to store elastic energy even when the shear is
weak and slow. The amount of energy stored decreases as the frequency is
lowered because then the lamellae have a chance to relax over larger length
scales. These effects are seen most clearly when the surfactants used are
amphiphilic polymers. To date the justifications for this behavior remain
complicated and tentative [31].
the stress for a given strain rate, as anticipated above. Surprisingly, the
alignment direction is not always in orientation (a) above. Instead, orient-
ation (b) can occur. Sometimes orientation (a) occurs for the lowest shear
rates and orientation (b) occurs at higher shear rates. This strange behavior
can be understood if one looks further into the shape of the layers being
sheared.
Surfactant layers fluctuate in shape; it is these fluctuations that give
rise to the Helfrich repulsion discussed earlier in this chapter. These shape
fluctuations persist for a certain time. If the shear occurs faster than this time,
one may encounter a situation like Fig. 7.13(a) when one tries to shear in
orientation (a). Instead of minimizing stress, the shear drives the bumpy
lamellae into their neighbors and increases the stress. The orientation (b)
avoids this situation. The shear no longer causes rippled lamellae to run
into each other.
What governs the persistence time for these ripples? We first recall
that the fluctuations that control the local distance between bilayers have
wavelengths λ much larger than the layer spacing d. Thus to gauge the
persistence time, we imagine a thought experiment involving a large scale
distortion. We pinch the layers together over long distance λ so that in the
pinched region the spacing is reduced by a small amount A from its equi-
librium value d. We then release the pinching force and watch the spacing
return to its equilibrium value. In order to restore the spacing, solvent must
flow into the pinched region. The flow causes dissipation, which retards
the relaxation. We denote the typical shear rate in this region as γ̇ and the
pinched volume by . During this flow, the energy is dissipated at a rate
Ẇ
η(γ̇ )2 , where η is the solvent viscosity.
The return of the spacing to its equilibrium value releases compressional
free energy. As we have seen, these fluctuating bilayers repel each other
because compression constrains their random fluctuations. According to
the Helfrich formula of Eq. (7.15) above, the energy per unit area u can be
expressed in terms of the temperature T and the bending stiffness κ: u
(T /d 2 )(T /κ). Thus the Helfrich energy U in the volume is given by
U
(/d)(T /d 2 )(T /κ). Even if there is additional attraction between
the layers that holds them at an equilibrium spacing, the Helfrich energy is
a substantial fraction of the energy and thus the overall energy is proportional
to this U .
We now express these energies in terms of our perturbed spacing d − A.
If A increases at a rate Ȧ, there must be an inward flow to supply the needed
volume of fluid. The flow enters a given lamellar gap at the boundary of
the pinched region, whose height is roughly d and whose circumference
is roughly λ. Thus the inward current I is roughly λdv. This current must
match the increase of volume λ2 Ȧ, so that v
Ȧλ/d. Finally, the shear
rate γ̇ is the gradient of velocity, or γ̇
v/d
Ȧλ/d 2 . Thus Ẇ
η(Ȧλ/d 2 )2 . To express the perturbation in the free energy, δU , in terms
of A, we imagine that d is close to an equilibrium, so that δU ∼ A2 . In
view of the order of magnitude of U given above, we conclude
Now we can estimate the relaxation rate Ȧ by requiring that the release rate
of stored energy, balance the dissipation rate:
Thus A relaxes in a time τ 1/[. . .]. We can rearrange the factors to write
40 μm
Fig. 7.15
Light micrograph of myelin structures
emerging from a drop of concentrated
C12 E5 surfactant in water, from [34],
courtesy Mark Buchanan
References
1. J. N. Israelachvili, Intermolecular and Surface Forces, 2nd ed. (London: San
Diego: Academic, 1992).
2. K. G. Denbigh, The Principles of Chemical Equilibrium: With Applications in
Chemistry and Chemical Engineering, 4th ed. (Cambridge: Cambridge, 1981).
3. D. L. Ho, R. M. Briber, R. L. Jones, S. K. Kumar, and T. P. Russell,
Macromolecules 31 9247 (1998).
4. These drawings are provided by the ChemIDplus
(http://chem.sis.nlm.nih.gov/chemidplus/) resource using
molecular data from the National Library of Medicine, 8600 Rockville Pike,
Bethesda, MD 20894 (http://sis.nlm.nih.gov).
5. J. Lang and R. Zana, in Surfactant Solutions: New Methods of Investigation,
ed. R. Zana (New York: Marcel Dekker, 1987).
6. P.-G. deGennes, Scaling Concepts in Polymer Physics (Cornell: Ithaca, 1979),
Chapter 1.
7. M. E. Cates and S. J. Candau, J. Phys.: Condens. Matter. 2 6869 (1990).
8. T. Tlusty and S. A. Safran, J. Phys.: Condens. Matter. 12 A253 (2000).
9. S. A. Safran, Statistical Thermodynamics of Surfaces, Interfaces, and Mem-
branes (Reading MA, USA: Addison-Wesley, 1994).
10. S. T. Hyde, in Handbook of Applied Surface and Colloid Chemistry, ed., Krister
Holmberg (New York: Wiley, 2001), Chapter 16.
11. V. Luzzati, P. Mariani, and T. Gulik-Krzywicki, in Physics of Amphiphilic
Layers, eds. J. Meunier and D. Langevin (Berlin: Springer, 1987).
212 Surfactants
activated hopping 24, 25 Bjerrum length 73, 128, 129, 178 condensation 24, 45, 46, 113, 114
activation barrier 53, 134 blob 86, 87, 90, 99, 100, 104, 105, 106, conjugated 43
activation energy 25 121, 124, 164–9, 186 conjugated polymers 43
addition polymerization 45 block copolymer fluids 3, 10 constant of proportionality 28, 74, 104,
additivity 175–6, 178 block of a copolymer 199 147
adsorbates 36, 152, 161 Boltzmann constant 7, 16 contact angle 152, 159, 171
adsorbed polymer 125–7, 165, 167, 168, Boltzmann distribution 15, 26 interfaces 156–8
170 Boltzmann factor 121 cooperative diffusion constant 100, 101,
affine deformation 28 Boltzmann normalization 22 110
aggregates 4, 5, 7, 113, 116, 123, 134, Boltzmann principle 16, 22, 65, 91, 144, cooperative diffusion and permeability
142, 145, 147 199 101
colloidal 8, 10, 79, 139, 140, 143, 146, Boltzmann probabilities 17, 128 copolymers 44, 201–2
147, 163 Boltzmann–Poisson Equation 128, 132 block copolymers 2, 10, 53, 202
diffusion-limited 134, 144, 145, 147 Boltzmann weight 68, 69, 177 diblock copolymers 44, 198
elasticity 147 bond angles see local structure peptide copolymers 45
fractal 113, 134 breaking-dominated relaxation 205 triblock copolymers 202
properties 146–8 Brownian motion 83, 90, 91, 98–9, 104, copolymerizing 44
reaction-limited 134, 141, 144, 145, 136, 140, 145 corona 125
147 of aggregates 134 correlation functions 76, 84, 101
silica 10 of the polymers 100 correlation length 86, 100, 102, 111, 120,
air–solid interface 157 of a sphere 99 121, 123, 127, 169, 186
air–water interface 152, 153 brush regime 124 correlation time 37
Alexander–de Gennes approximation Coulomb energy 74, 75, 131
124 Coulomb repulsion 73, 127, 130, 207
Alexander, M. 124 capillary flows 160 counterions 43, 129, 130, 131–2, 136
alpha helix 43 capillary forces 157, 163, 173, Cox–Merz rule 109
amphiphilic molecule 11, 163, 173, 198 capillary solid 158 creaming 116
amphiphilic polymers 174, 198–202, 207 carbon–deuterium bonds 39 critical micelle concentration (CMC)
as surfactants 203 cartesian coordinate 14 181, 183, 185, 192
anisotropic interactions 134 cationic surfactant cetyl trimethyl critical mixture 162
annealed variables 24 ammonium bromide (CTAB) 9, crosslinking 43, 44
Argon 28 179, 180, 184 curvilinear diffusion constant 104, 108
Arrhenius plot 25 chain swelling 127 Cyrus, S. 36
“association” of mesoscopic structures charge renormalization 130, 131, 133
10, 43 charge separation 127
asymptotic dependence 6 chemical potential 22 Debye–Hückel approximation 129,
asymptotic equation 52 Cis–trans isomerization 44 130–1
atactic structures 44 configuration energy 20 Debye–Hückel equation: 128, 129
atomic collision times 31 coil-stretch transition 52, 53 Debye–Hückel potential 129
atomic force microscope (AFM) 35 collision 14, 24, 29, 31 Debye–Hückel screening 128, 129, 135
atomic polarizability 115, 175, 176 colloidal aggregate 7–8, 10, 79, 139–46, Debye screening length 128, 130, 131–3
attempt rate 25 163 deformability
average local density 54 a form of self-organization 7 of polymers 9, 96
colloidal crystals/crystallization 132, 136 of rubber 44
colloidal dispersions 113, 125 of surfactants 180, 181
Bakajin, O. 76 colloidal motion 135 deformation of random chains 8
ballistic aggregation 145 colloidal particles 4, 113, 120 De Gennes, PG 29, 64, 104, 124, 126,
Ball, R.C. 141, 147 colloidal stabilization 113–14 133, 158, 159, 167, 191, 194
baroclinic mode 209 colloidal suspensions see colloidal degrees of freedom 14, 15, 19, 24, 117,
beam direction 36, 90 dispersions 164, 184, 193
bending modulus 189, 195, 203 complete nonwetting 157 density correlation function 58
bending stiffness 189, 193, 203, 207, complete wetting on high-energy solid depletion attractions between colloidal
208, 209 surfaces 157 particles 121
bicontinuous microemulsions 194, 196 compressibility 92 Derjaguin approximation 116, 117, 123
bilayers 182, 203, 208, 209 concentration profile 99, 101, 124, 162, Des Cloizeaux J. 56, 64
biosynthesis 45 165–7, 168, 203 dewetting 160
214 Index
macroscopic responses: properties of Oseen tensor 93–4, 94–5, 95–6, 99 living polymerization 46
structured-fluid systems 31–4 origin 109–10 polymers 7, 32–3, 86, 96, 41, 66
magnetic resonance 38 osmotic pressure 33, 69, 70, 83–5, 87, flexibility 41
Mandelbrot, B.B. 55 101–3, 109–11, 121, 124, 127, 132, in good-solvent conditions 70
Marangoni flow 160 146, 165, 168, 186, 204, 206 in poor-solvent conditions 70
mass distribution 63, 79, 143, 144 of hard-sphere gas 84 randomness 41
mean curvature 188, 189, 191, 192, 203 variation 100 random-walk 41, 47, 53, 60, 64, 69,
mean-squared displacement 91, 103 osmotic repulsion 123, 125 90, 94, 188
mean-squared fluctuations 89 osmotic sensitivity 33 spatial structure 86
melt state 8, 90, 200 overcharging 131 spontaneous internal motions 98
solvent-free limit 8 overlap concentration 85, 87, 99, 148, theta condition 70
membrane fluctuations 191–2 185 polymer screening 81
Metropolis algorithm 27 polymer solution 89–90, 99
Metropolis dynamics 27 packing length of monomers 200 turbidity 89
micelle 180 “pair correlation function” 84 motion in 90–109
micelle interaction 186–92 parallel polymers 53 polypeptides 42, 45
micelles 9, 10, 11, 130, 173, 180–5 partial wetting 157 polypropylene 42
non-repelling micelles 187 particle–particle interaction 7 polysaccharides 45
spherical micelles 130, 187 partition function 20, 39 polystyrene sulfonate 43
two-dimensional micelles 185, 188–91 pentaethylene glycol monododecyl ether polystyrene: unsaturated hydrocarbon 42
wormlike micelles 182, 184, 187–8, (C12 E5) 179 poor solvent 70
202, 204, 206 permanent dipoles 177–8 positivity: miscibility in liquids 175
microemulsion 3, 10, 34, 174, 190–5, permeability 101, 110 potential of mean force 19, 24, 67–9
203 persistence length 48, 75, 184, 187, 191, power-law regime 168
or equilibrium dispersed state 193 194, 204 principal curvatures 188
miscibility in liquids: principles perturbation–attraction theorem 117–19, probability current 28, 91
controlling 174–8 120, 126–7, 137–9 probability and work 17–20
modulus 74 pervaded volume 7, 54, 62, 78, 80, 81, probes 34, 38
Mohanty, J. 114 147 of atomic environment 38
molecular-layer oscillations 34 phantom attraction 163 of spatial structure 34–7
molecular polarizability 175 phantom polymers 48 of structured fluids 31
moments 78, 81 phase stability 33, 194
momentum-absorbing blobs 99 phase transition theory 140
quantum tunneling 26
Monge representation 189, 191, 210 phospholipid 1-palmitoyl-2-
quenched variables 24
Mondello, M. 8 oleoylphosphatidylcholine (POPC)
monomer-centered spheres 78 179
Monte-Carlo simulations 26, 76 pinned contact line 159 radius of gyration 60, 64, 72, 85, 97,
multifractals 79 plateau region 108 121, 124
multilamellar vesicle morphology 209 plumber’s nightmare 188 randomly branched polymer 146
multivalent counterions 131 polarity 42 random thermal fluctuations 24, 101
myelin morphology 210 polarizability 114 random-walk 184, 200, 205
polyacetylene 43, 48 avoidance 145
nematic liquid crystal 133, 187 polyamide nylon 45 chain 62, 74, 87
“network fluid” 10 polyatomic solutes 163–71 fluctuations 204
neutral wetting 157 polybutadiene 42 history 144
neutrons 37, 38, 56, 57 1,2/1,4 polybutadiene 42, 46 molecules 200
neutron scattering intensity 196 polycondensation 46, 47 properties 7
neutron spin echo technique 38 polydimethyl siloxane 42 self-repelling 7, 164
Ninham, B.W. 114 polydispersity 142 reaction-limited aggregates 145, 147
non-scalar perturbations 120 polyelectrolytes 9, 43, 73–5, 128 reaction-limited regimes 144–6
normalized probability 21 polyesters 42 real system 19
normal stress 108 polyethylene 41 reciprocity: miscibility in liquids 175–7
nuclear magnetic resonance (NMR) 38 polyethylene glycol (PEG) 42 recombinant DNA technology 45
polyethylene oxide 41, 42 renormalized variables 50
onion morphology 209 polyethylene–propylene 44 reptation 104–5, 107, 205
Onsager, L. 133 polyisoprene 42 leads to stress relaxation 107
opaque structures 63, 65, 67, 84, 96–9, polymer deformability 96 model 103, 104, 105
135, 146 polymer melt 32, 85, 106, 107, 199 repulsive forces 113, 122–31, 164
open aggregation 182–5 polymeric amphiphiles 198, 202 requirement of homogeneity 51
optical fluorescence of small dye in solution 203 reservoir configurations 15
molecules 39 polymeric solvents and screening 79–81 Roe, R.J. 8
ordering: miscibility in liquids 176–7 polymer-induced interactions 121, 124 Rouse diffusion constant 104
order of magnitude 11, 27, 30, 31, 96, polymerization, types 41–7 Rouse model 96, 104
115, 161, 208 addition polymerization 45 Roux, D. 36
oscillating strain 32 condensation polymerization 45–6 rubber 4, 8, 28, 39, 42, 44, 106, 146, 151
216 Index