Important Points from Last Lecture:

1. Important Points from Last Lecture: • For mixing to occur, the free energy (F) of the system must decrease; DFmix < 0. • At a constant T, the change in free energy upon mixing is determined by changes in internal energy (U) and entropy (S): DFmix = DU - TDS. • The c interaction parameter is a unitless parameter to compare the interaction energy between dissimilar molecules and their self-interaction energy. • The change of DFmix with f, c (and T) leads to stable, metastable, and unstable regions of the phase diagram. • For simple liquids, with molecules of the same size, assuming non-compressibility, the critical point occurs when c = 2. • At the critical point, interfacial energy, g = 0.

2. Constructing a Phase Diagram Spinodal where: Co-existence where: c >2 c =2 T1 T2 T3 T4 T5 fG T1<T2<T3….

3. Phase Diagram for Two Liquids Described by the Regular Solution Model Low T Immiscible Miscible High T Spinodal and co-existence lines meet at the critical point. fG

4. PH3-SM (PHY3032) Soft Matter Lecture 6 Polymer Phase Separation, Morphologies and Interfaces ● An Introduction to Colloids 8 November, 2011 See Jones’ Soft Condensed Matter, Chapt. 3, 9 and 4

5. Free Energy of Mixing for Polymers Polymers consist of N repeat units (or “mers”) each of length a. a The thermodynamic arguments applied to deriving DFmix for simple liquids can likewise be applied to polymers. The derivation must consider the connectivity of units when putting them on a lattice. N units are mixed all at once rather than individually.

6. The free energy change per polymer molecule is therefore: Free Energy of Mixing for Polymer Blends, DFmixpol We start with our expression for free energy of mixing per liquid molecule, DFmix, for simple liquids: When arranging the repeat units on the lattice, the probability is determined by the volume fractions,f, of the two polymers (assuming equal-sized units). N has no influence on DSmix per polymer molecule. But the change in U upon mixing polymers must be a function of N times the DUmix for each of the repeat units.

7. In general terms, for an asymmetric blend of chains with two different degrees of polymerisation (NA and NB), we can write the Flory-Huggins expression, where f is the volume fraction of one of the polymers: The points on the coexistence curve are found from the tangent found when the differential is set = 0: Polymer Phase Separation As the polymer consists of N repeat units (mers), we can find the free energy of mixing per mer by dividing through by N:

8. Coexistence Curve in a Symmetric Blend In a symmetric blend, where the chains have an equal number of repeat units: NA = NB = N The values ofcthat satisfy this equation define the co-existence curve, which is the line that separates the one-phase and two-phase regions of the phase diagram, are: This equation can be used to plot cN as a function of f on a phase diagram to define the coexistence curve.

9. Spinodal Line for a Polymer Blend The spinodal line is found by equating the second differential to 0: The curve corresponding to the inflection points defines the boundary between the unstable and metastable regions, called the spinodal line. Solving for c gives us: For given values of NA and NB, the spinodal line will follow the dependence on f given by this equation. For a symmetric blend (NA = NB), we see that the spinodal passes through cN = 2 when f = ½.

10. Spinodal Line for a Polymer Blend Then, Phase separation Single phase (blend) is stable The critical point can be found from For a symmetric blend, the phase diagram is symmetric, and the critical concentration is atfc= 1/2 or ccN = 2 at the critical point For polymers, cN is the key parameter - rather than c as for simple liquids.

11. Polymer Phase Separation in a Symmetric Blend NA = NB = N cN fA As c ~ 1/T, we can either increaseN or decreaseT to cause miscible polymers to phase separate.

12. Polymer Immiscibility N typically has a value of 1000 or more, so that cc = 2/N is very small. Entropic contributions in polymers that encourage mixing cannot easily compensate for unfavourable energies of mixing. (Remember that in simple liquids mixing will occur up to when c = 2 as a result of entropy.) Therefore, polymer mixing (miscibility) is only favoured when c is negative or exceedingly small or when N is very small (such that cN<2). This fact explains why most polymers are immiscible - making them difficult to recycle, unless they are blended with very similar molecules (low c) or else have a low N. Polymeric interfacial structure and phase separation are often studied by neutron scattering and reflectivity.

13. When moving from the one-phase to the unstable two-phase region of the phase diagram, ALL concentration fluctuations are stable. F . Fo fo 1 0 f1 f2 Structures Resulting from Phase Separation in the Unstable Region (cNif polymers!) Spinodal points define the unstable region boundaries. c Leads to “spinodal decomposition” f

14. Two-Phase Structure Obtained from Spinodal Decomposition The two phases have a characteristic size scale defined by a compromise. Fourier transform of image If the sizes of the phases are too small: energy cost of extra interfaces is too high. If the phases are quite large, it takes too long for the molecules to travel the distances required for phase separation. Poly(styrene) and poly(butadiene) undergoing spinodal decomposition. Simulation of spinodal decomposition: http://youtu.be/NSpOX9mfX3g

15. F Ff1 . f2* Fo Small fluctuations in composition are not stable. f1 fo f2 Structures Resulting from Phase Separation in the Metastable Region (cNfor polymers!) Free energy change (per unit volume) on de-mixing: c DFv = Fo - Ff1 f Onlyf1 and f2* are stable phases! The f2* composition must be nucleated and then it will grow.

16. Phases grow in size to reduce their interfacial area in a process called “coarsening”. Structures Obtained from Two Immiscible Polymers via Nucleation Poly(ethylene) and poly(styrene) blend AFM image 10 mm x 10 mm

17. Comparison of “Nucleated” and “Spinodal” Structures Nucleated structure: islands of one phase in another Spinodal structure: co-continuous phases From G. Strobl, Polymer Physics, Springer

18. Energy reduction through phase separation with growth of the nucleus with volume (4/3)pr3 Energy “cost” of creating a new interface with an area of4pr2 Nucleation of a Second Phase in the Metastable Region Growth of the second phase occurs only when a stable nucleus with radius r has been formed. f2* r f1 gis the interfacial energy between the two phases.

19. F* r * The free energy change in nucleating a phase, DFnucl, is maximum for a nucleus of a critical size, r *. Critical Size for a Stable Nuclei + DFnucl r - If r < r *, further growth of the nucleus will raise the free energy. The nucleus is unstable. If r > r *, the nucleus is stable, and its further growth will lower the free energy of the phase-separating system.

20. We can find the maximum of DFnucl from: Solving for r, we see: Substituting in our value of r *, we can find the energy barrier to nucleation: Simplifying, we see: Calculating the Size of the Critical Nucleus, r *

21. This probability is given by a Boltzmann factor: Estimating the Rate of Nucleation during Phase Separation Nucleation occurs when a fluctuation in DF during the formation of a nucleus is > DF *. The rate of nucleation is determined by the frequency of the fluctuations and their probability of exceeding DF*. The temperature dependence is complicated by the fact that DF* is a function of g and DFv, which are both temperature dependent. The barrier DF * can be lowered by the presence of a “nucleant” (a surface on which the phase can grow) in heterogeneous nucleation.

22. The system will separate into two “bulk” phases; droplets of any size are not favoured. Significance of Surface Tension,g droplet If g >0, then the system can lower its free energy by reducing the interfacial area: DF = gDA But if g = 0, then mixing of droplets - or molecules - does not “cost” any energy (no work required!). Thus, mixing is favoured at the critical point.

23. A B Polymer Interfacial Width, w w w The interface between immiscible polymers is never atomistically sharp. If the molecules are forbidden from crossing the boundary, their number of conformations (i.e. arrangements) is reduced. The entropy will then decrease.Entropy increases with interfacial broadening. Therefore, an interfacial width, w, can be defined for any polymer interface.

24. Critical angle Neutron Reflectivity from a Single Interface q q Usually a distribution of l are used to obtain a broad range of Q values.

25. Inversely related to film thickness Sensitivity of Neutron Reflectivity to Interfacial Roughness Polymer film on a Si substrate with increasing surface roughness,s.

26. Scattering density profile Reflectivity from a Polymer Multi-Layer Deuterium is often used to label polymers to increase their “contrast” with hydrogenated polymers when scattering neutrons.

27. Comparison of Polymers with DifferentcParameters and Interfacial Widths w Scattering density profile

28. Width between Two Polymer Phases when Approaching the Critical Point Experiments on immiscible polymers confirm that the interface broadens as the critical point is approached. Also as cN decreases toward 2, g approaches 0. Data from C. Carelli et al. Phys. Rev. E72, 031807 (2005)

29. Colloids Particles are much larger than the size of molecules. 1 mm Because the size of colloidal particles is on the order of the wavelength of light, they offer some interesting optical characteristics.

30. Optical Characteristics of Colloidal Films Colloidal particles are dispersed in a continuous matrix. l Diffraction condition: q q d

31. The effect is a result of the opal structure, which consists of silica spherical particles (typically 250 - 400 nm in diameter) about 1/2 the wavelength of light, leading to diffraction of the light by the regular spacing. Colloids in Nature: Opals Natural opal reflects various colours of light depending on the viewing angle.

32. The particles are then dissolved to leave a network of air voids. The space between the particles is filled with a solid through infiltration or deposition from the vapour phase. Using Colloids to Create “Inverse Opal” Structures Colloidal particles are packed into an ordered array. •Useful optical and magnetic properties. • Inverse opals have “optical band gaps”

33. Phase Equilibria in Colloidal Dispersions Equilibrium: Non-equilibrium (glassy): (Volume %) RCP = random close-packing; HCP = hexagonal close-packing Mono-sized particles can become ordered into crystals at f = 0.54 while still in the “wet” state. V. Prasad, D. Semwogerere and Eric R. Weeks, J. Phys.: Condens. Matter19 (2007) 113102 (25pp)

34. Forces Acting on Colloidal Particles • Drag force from moving through a viscous medium • Gravity: leads to sedimentation or creaming • Random, “thermal” forces from molecules: lead to Brownian motion • Coulombic: can be attractive or repulsive; screened by the intervening medium • van der Waals’: attractive for like molecules • Steric: caused by intervening molecules that prevent close approach

35. FS The Viscous Drag Force on Colloids The Reynolds’ number, Re, compares the effects of inertia (momentum) to viscous drag: V = velocity a r= density of the continuous medium h= viscosity of the continuous medium • When Re is low, viscosity dominates over inertia. There is no “coasting”! • At low Reynolds’ number, when rva << h, the viscosity of the liquid imposes a significant drag force on the particle’s movement. • The Stokes’ equation gives this force (Stokes’ or drag force) as: • Fs = 6phav • Observe that Fs applies when h is large in comparison to a and v.

36. Effect of Gravity on Particle Velocity FS • If the density of a particle is different than that of the surrounding fluid, it will be subject to a gravitational force, Fg, leading to settling (or rising). • If the difference in density is Dr(+ or -), then Fg = (4/3)pa3Drg, where g is the acceleration due to gravity. • At equilibrium the forces balance: Fs = Fg. • So, 6phav = (4/3)pa3Drg • The velocity at equilibrium, i.e. the terminal velocity, vt, is then found to be (2a2Drg)/9h. • Larger particles will settle out much faster than smaller particles - giving us a means to separate particles by size. The same principle applies for separation by size using centrifugation. a Fg

37. Experimental Observation of Brownian Motion Phenomenon was first reported by a Scottish botanist named Robert Brown (19th cent.) Brown observed the motion of pollen grains but realised that they were not living. Brownian motion Computer Simulation: http://youtu.be/6VdMp46ZIL8

38. Effect of Molecular Momentum Transfer: Random Brownian Paths Self-similarity: paths appear the same on different size scales 2-D representations of 3-D particle trajectories

39. Mean-Square Displacement Displacements of a large number of particles Mean displacement: Mean-square displacement: Root mean-square displacement:

40. Start Finish n 3 2 1 Then when observed over n steps at equal time units, the average particle displacement for several “walks” will be 0, but the mean square displacement is non-zero: Distance Travelled by Particles If in every unit of time, a particle takes a step of average distance, , in a random direction... Random walk Thus the mean-square displacement is proportional to the number of steps and hence to the time.

41. Einstein and Smoluchowski wrote an equation for the equation of motion for a Brownian particle in which the net random force exerted by the fluid molecules, Frand, balances the forces of the particle: where x is a drag coefficient equal to 6pha for an isolated, spherical particle in a viscous fluid. Writing and in terms of we see: Equation of Motion for Brownian Particles Einstein was unaware of Brown’s observation, but he predicted random particle motion in his work on molecular theory.

42. The Mean-Square Displacement And we see that Then, multiplying through by x: Substituting in for the first term, we find: If random, the mean displacements in the x, y and z directions must be equal, so But we recognise that:

43. The Mean-Square Displacement After substituting an identity and taking the average of each term: 0 kT This leaves us with: Because Frand, x and v are uncorrelated, the first two terms on the r.h.s. average to zero. Finally, the equipartition of energy says that for each d.o.f., (1/2)mv2 = (1/2)kT in thermal energy.

44. Integrating and multiplying by three, E and S thus showed that the mean squared displacement of a Brownian particle observed for a time, t, is A diffusion coefficient, D, which relates the distance (in 3D) to the time of travel, is defined as So it is apparent that Recall Stokes’ equation, x = 6pha for a spherical particle in a fluid. The Stokes-Einstein diffusion coefficient is thus: The Stokes-Einstein Diffusion Coefficient

45. Applications of the Stokes-Einstein Equation • Observe that the distance travelled, R (root-mean-square displacement, <R2>1/2) varies as the square root of time, t1/2. • Early work assumed that the distance should be directly proportional to time and made data interpretation impossible. • Experiments, in which the displacement of colloidal particles with a known size was measured, were used by Perrin to determine the first experimental value of k. • Brownian diffusion measurements can be used to determine unknown particle sizes. • The technique of light scattering from colloidal liquids is used to find particle size through a diffusion measurement.