The Flory-Huggins Theory of Polymer Solutions. Purdue University Biomedical Engineering. Professors Kinam Park & Luis Solorio

Transcription

1 The Flory-Huggins Theory of Polymer Solutions Purdue University Biomedical Engineering Professors Kinam Park & Luis Solorio

2 Thermodynamics of Ideal Solutions G m G Mixture_ A&B G A G B G H TS G m H m TS m Ideal Solution: Mixture of moleules that are identical in. size. energy of like and unlike molecular interactions are equal Results in athermal mixing, or ΔH m =0 G m TS m Solutions are homogenous mixtures of two or more components on the molecular scale G 0 for spontaneous events, Gibbs free energy is related to enthalpy and entropy Driven towards creating the greatest chaotic dispersal Gasses don t spontaneously contract How does athermal mixing affect rotation, vibrational, and translation of molecules in the system?

3 Boltzmann s Equation S k ln S comb m k ln ideal _ mix solvent solute S comb m k ln ideal _ mix ln solvent ln solute solvent solute S m comb k ln ideal _ mix Solvent Lattice Solute From statistical mechanics S=entropy Ω=number of distinguishable equal energy states k=boltzmann constant Identical pure substances have only spatial arrangement

4 Ideal Liquids Sterling's approximation: ln!=ln- (for significantly large) S comb m k ln! kln!!!!! S comb m k ln ln Latice_ size Solvent Solute ideal _ mix!!! A ( ) and n i = i S comb m kn ln X n ln X k R A S m comb Rn ln X n ln X G m TdS RT n ln X n ln X

5 Ideal Gasses G H TS 0 dg dh SdT TdS dsdt H U PV 0 dh du PdV VdP dvdp dg du PdV VdP SdT TdS du TdS PdV dg TdS TdS PdV PdV VdP SdT dg VdP SdT At Constant Temperature dg=vdp Δ Δ dg V nrt P dg nrt VdP P P G nrt ln P P G molar RT ln P P dp P

6 Raoult s Law G m TdS RT n ln X n ln X Look at a single component system: G RTn ln X G n P,T G n G molar G molar RT ln X i i 0 RT ln X i G molar RT ln P P Raoult s Law X i P i P 0 Raoult s law is used to define an ideal solution It is used to identify how solutes affect factors such as changes in boiling temperature, freezing temperature, or to calculate the osmotic pressure It can also be used to determine how the solute is affecting the vapor pressure of the solvent above the solution Could these equations be used to accurately describe a polymer system?

7 Paul J. Flory, Ph.D. The obel Prize in Chemistry Born on 9 June, 90, in Sterling, Illinois. B.S. Manchester College in Indiana, 93. Ph.D. Physical Chemistry, Ohio State University, 934. I joined the Central Research Department of the DuPont Company. There it was my good fortune to be assigned to the small group headed by Dr. Wallace H. Carothers, inventor of nylon and neoprene, and a scientist of extraordinary breadth and originality. It was through the association with him that I first became interested in exploration of the fundamentals of polymerization and polymeric substances. His conviction that polymers are valid objects of scientific inquiry proved contagious. The time was propitious, for the hypothesis that polymers are in fact covalently linked macromolecules had been established by the works of Staudinger and of Carothers only a few years earlier. A year after the untimely death of Carothers, in 937, I joined the Basic Science Research Laboratory of the University of Cincinnati for a period of two years. With the outbreak of World War II and the urgency of research and development on synthetic rubber, supply of which was imperiled, I returned to industry, first at the Esso (now Exxon) Laboratories of the Standard Oil Development Company (940-43) and later at the Research Laboratory of the Goodyear Tire and Rubber Company (943-48). Provision of opportunities for continuation of basic research by these two industrial laboratories to the limit that the severe pressures of the times would allow, and their liberal policies on publication, permitted continuation of the beginnings of a scientific career which might otherwise have been stifled by the exigencies of those difficult years.

8 Paul J. Flory, Ph.D. The obel Prize in Chemistry 974 In the Spring of 948 it was my privilege to hold the George Fisher Baker on-resident Lectureship in Chemistry at Cornell University. The invitation on behalf of the Department of Chemistry had been tendered by the late Professor Peter J.W. Debye, then Chairman of that Department. The experience of this lectureship and the stimulating asociations with the Cornell faculty led me to accept, without hesitation, their offer of a professorship commencing in the Autumn of 948. There followed a most productive and satisfying period of research and teaching "Principles of Polymer Chemistry," published by the Cornell University Press in 953, was an outgrowth of the Baker Lectures It was during the Baker Lectureship that I perceived a way to treat the effect of excluded volume on the configuration of polymer chains. I had long suspected that the effect would be non-asymptotic with the length of the chain; that is, that the perturbation of the configuration by the exclusion of one segment of the chain from the space occupied by another would increase without limit as the chain is lengthened. The treatment of the effect by resort to a relatively simple "smoothed density" model confirmed this expectation and provided an expression relating the perturbation of the configuration to the chain length and the effective volume of a chain segment. It became apparent that the physical properties of dilute solutions of macromolecules could not be properly treated and comprehended without taking account of the perturbation of the macromolecule by these intramolecular interactions. Out of these developments came the formulation of the hydrodynamic constant called theta, and the recognition of the Theta point at which excluded volume interactions are neutralized. Criteria for experimental identification of the Theta point are easily applied. Ideal behavior of polymers, natural and synthetic, under Theta conditions has subsequently received abundant confirmation in many laboratories. These findings are most gratifying. More importantly, they provide the essential basis for rational interpretation of physical measurements on dilute polymer solutions, and hence for the quantitative characterization of macromolecules.

9 Polymer Solution. Flory-Huggins Mean Field Theory S k ln Polymer ideal _ mixture Solvent S comb m k ln ideal _ mixture Polymer Solvent S comb m k ln ideal _ mixture Polymer Solvent Lattice Solute Spatial arrangement of polymers is increased due to the different chain segment conformations 0 = + x 0 : number of lattice positions : number of solvent molecules : number of polymer molecules x : ratio of molecular volumes of solvent and polymer(~number of mers)

10 Polymer Solution. Flory-Huggins Mean Field Theory Mean field approximation: Effect of the other individuals on any single individual is approximated by a single average effect ideal _ mixture Continuous product of ν, i ν=possible conformations of the polymer! i in the lattice /! accounts for identical indistinguishable conformations (over counting) Solvent Lattice Each polymer molecule occupies x lattices, and i polymer molecules have already been placed in the lattice. 0 -x i : The number of vacant positions: f i : The fraction occupied by i molecules ( - f i ) : The fraction of remaining vacant positions z. ( - f i ) : The number of ways arranging the next segment Repeating unit (z-). ( - f i ) : The number of ways arranging the third (and all z neighboring lattice sites the remaining) segment (Z-)(-f i ) : The number of ways arranging the (i+)th polymer molecule ( o xi) Z( f )( )( )( )( ) i z fi Z fi i st segment nd segment 3rd segment 4th segment... f i 0 xi z(z ) x- (z-) x- 0

11 Polymer Solution. Flory-Huggins Mean Field Theory 0 0 x x i Z xi V The Stirling approximation 0 0 0!! x i Z i x xi V V i+ : The number of configurations of just one polymer molecule in the lattice : The number of ways to place the indistinguishable polymer molecules 0!! i i i i V V !!!!!! x x Z Z x There are many fewer configurations possible for the polymer solutions compared to small-molecule solutions. ideal _ mixture z (x) Polymer x e x Polymer Solvent solvent x Polymer Polymer Total conformations of the polymer molecules with complete freedom in the lattice Conformations when polymer molecules are competing for cells in a pure amorphous polymer The number of conformations of the polymer in a solvent Ω polymer is provided by the first two factors

12 Polymer Solution. Flory-Huggins Mean Field Theory The configurational entropy represents the entropy of mixing of the perfectly ordered pure solid polymer, for which S = 0, with pure solvent. S c k ln ln x x e Z x ln S c = S dis + S mix Disordering & dissolution of a polymer. Formation of an amorphous polymer (disorientation) S dis. Dissolution of the amorphous polymer in the solvent S mix Z Sdis lim Sc k ln x x ln Smix 0 e k ln x x ln x x

13 Polymer Solution. Flory-Huggins Mean Field Theory Assume x V V S mix R n ln n ln Polymer molecules S mix R n lnx n ln X Small molecules n : number of moles ϕ : volume fraction G=H-TS This provides the entropy and is sufficient for the cases when no enthalpy occurs However that is an unrealistic assumption with polymers

14 Polymer Solution. Flory-Huggins Mean Field Theory Mixing Enthalpy: Interaction energy between solvent molecules and solute segments,,, w, w, w w H mix RT n The energy change (in units of RT) that occurs when a mole of solvent molecules is removed from the pure solvent (where ϕ = 0) and is immersed in an infinite amount of pure polymer (where ϕ = 0). H mix < 0 : Exothermal H mix = 0 : Athermal H mix > 0 : Endothermal χ = z w /RT : Interaction parameter

15 Polymer Solution. Flory-Huggins Mean Field Theory G mix G mix RT 0 n n ln n ln n n ln n ln Polymer dissolves in solvent Polymer does not dissolve in solvent ϕ and ϕ are always lower than. Thus, lnϕ and lnϕ have always negative value. For complete miscibility over all concentrations, for the solute-solvent pair at the temperature (T) of interest must be less than 0.5. If > 0.5, then G mix > 0 and phase separation occurs If < 0.5, then G mix < 0 over the whole composition range. The temperature at which = 0.5 is the θ temperature. o G n o mix p, T, n o : Pure solvent activity : Solvent activity in solution RTln( ) x

16 Polymer Solution. Flory-Huggins Mean Field Theory ϕ G mix RT n n ln n ln and ϕ are always lower than. Thus, lnϕ and lnϕ have always negative value. Flory-Huggins equation makes it possible to account for the thermodynamic properties of polymer solutions Accounts for the large derivations form Raoult s Law It is great for predicting general trends, but theory still deviates from the experimental data Deviation occurs as a result of the necessary assumptions needed to create a solution -Does not account for volume change during mixing -Does not exclude self-intersections of the chain -Mean field approximations are more accurate for high polymer concentration Regardless, it was a huge step forward in understanding how polymers behave

17 Application: Polymer Solubility Polymer with higher molecular weight precipitates more readily than the polymer with lower molecular weight. Add non-solvent until solution becomes turbid Heat, cool slowly and separate precipitate Finite drop in temperature always renders finite range of molecular weight insoluble Some high molecular weight polymer will also remain soluble! Solubility of polymers generally decreases with increasing molecular weight. T p : Precipitation temperature at which detectable turbidity can be first observed. T p phases phase Mol. Wt Volume fraction polymer

18 The θ Conditions Described in Movies

19 Viscosity of Polymer Solutions

20 Dynamic Viscosity (Absolute or Shear Viscosity) Dynamic viscosity is a measure of the resistance of a fluid to gradual deformation by shear stress. The resistance is a result of friction caused by sliding layers of fluid. c y m / m s ewton s Law of Friction / m s s Kg m m m s : Shear stress in fluid in Pascal (/m ) : Dynamic viscosity of fluid (. s/m ) : Shear rate (s - ) = c/ y c: Unit velocity (m/s) y: Unit distance between layers (m) s m Pa s 0P 0 Pa s (Pascal second) P: Poise g cm s Pa s: Two plates are separated by m with a fluid in between, and one plate is moving parallel to the other plate with a shear stress of pascal to move m/s. If m/s, then the viscosity is ½ Pa s. The viscosity of water = cp ( = 0.0 g/cm. sec = mpa. sec = 0.00 Pa. sec = s/m ).

21 Viscosity of Polymer Solutions Viscosity of a polymer solution can be measured by various methods. Viscosity of a polymer solution is a function of temperature, solvent, polymer concentration, and polymer molecular weight. When a polymer molecular weight is determined Viscosity of a polymer solution provides insightful information on the configuration of the polymer molecules in solution. At high polymer concentrations, polymer molecules tend to entangle each other, and this affect the measurement of viscosity. Thus, it is necessary to estimate the effect of solvent viscosity and polymer concentration on the overall viscosity. Since the viscosity measurement procedures do not remove the effects of polymer-solvent interactions, the viscosity average molecular weight ( ) depends on some extent on the solvent used In a dilute solution, the disturbance of the flow pattern of the suspending medium by one particle does not overlap with that caused by another. Solvent trapped inside a polymer molecule has lower viscosity than the bulk solution. Thus, the polymer coil behaves as an impenetrable sphere.

22 Viscosity of Polymer Solutions Viscosity of spheres and random coil molecules. where A is the Avogadro s number, c is the polymer concentration (g/ml), M is the polymer molecular weight, and V h is the volume of a solvated molecule. Relative viscosity ( rel ): Specific viscosity ( sp ): The specific viscosity accounts for the effect of solvent viscosity. The specific viscosity ( sp ) is a suitable indicator for the changes in shape. Spherical polymer has low resistance to flow, and so, small sp. Stretched, rod-like shape has higher resistance to flow, leading to high sp.

23 Viscosity of Polymer Solutions Reduced viscosity ( red ): The concentration effect is adjusted and it is valid only for dilute solution. Intrinsic viscosity ( ) lim The contribution of the individual macromolecules to the viscosity increase will be independent and additive only when the polymer molecules are infinitely far from each other. A and M are constants, so V h can be calculated from the viscosity measurement. Huggins Equation: The reduced viscosity is related to [η] by a power series of the form. The constant k H is termed the Huggins constant and has values ranging from 0.3 in good solvents to 0.5 in poor solvents. It contains information about hydrodynamic and thermodynamic interactions between coils in solution.

24 Solubility Parameters Poly(lactide-co-glycolide) (PLGA)

25 Solubility Parameters Solubility parameters are frequently used to predict compatibility between polymers and solvents. The dissolution rate and extent of polymer dissolution in a solvent depend on the thermodynamic compatibility of the polymer with solvent and the T g of the polymer, the nature of solvent (solubility parameters), temperature, chemical composition and molecular structure of the polymer, and crystalline form of solid polymorphism. The solubility of a given polymer in different solvents largely depends on its chemical structure. Polymers dissolve in solvents whose solubility parameters are similar to their own, and in general, structural similarity favors solubility, known as the like dissolves like principle. Dissolution of a polymer in a solvent requires the energy (E vap ) to break all intermolecular physical links in a unit volume, which is known as the cohesive energy density (CED): where H vap is the enthalpy of vaporization, RT is the thermal energy, and V m is the molar volume of the polymer. Thus, CED is the energy of vaporization in calories (or Joules) per cubic centimeter. CED reflects the degree of van der Waals forces that hold the molecules together.

26 The Hildebrand Solubility Parameter (Total Solubility Parameter), t The Hildebrand solubility parameter (or total solubility parameter), t, is defined as: / / H vap of polymers cannot be measured directly, and thus, indirect methods are used to calculate the solubility parameters of polymers. The enthalpy of mixing a polymer with a solvent per unit volume, H mix /V mix, is given by:,, where V mix is the volume of mixing, and ϕ and ϕ are the volume fractions of solvent and polymer, respectively. As a general rule, (δ t, - δ t, ), the difference in solubility parameters of a solvent and a polymer, must be small for dissolution to occur over the entire volume fraction range. In fact, it is related to the Flory-Huggins parameter χ by the following relation:,, The Hildebrand solubility parameter does not account for any specific interaction, especially hydrogen bonding. Thus, if polar contributions are significant to the total CED of a solvent, the Hildebrand solubility parameter may not be able to explain the lack of interaction between the solvent and the polymer.

27 The Hansen s Partial solubility parameters Since slight differences in polar contributions cause significant differences in solubility behavior, different types of polar contributions were identified to more accurately reflect the solubility behavior of solvents. The Hildebrand solubility parameter can include the hydrogen bonding based on the assumption that solubility is highest if the two materials have similar polarities. In this approach, solvents are divided into three categories of hydrogen bonding, such as poor (or weak), moderate, and strong hydrogen bonding capabilities. Hansen also developed solubility parameters based on three specific interactions: dispersive (non-polar), polar cohesive (dipole-dipole), and hydrogen bonding interactions. Thus, it is called Hansen s three-component solubility parameters (or partial solubility parameters) in comparison with Hildebrand s total solubility parameter. In this approach, the cohesive energy (E) has dispersive (d), polar (p), and hydrogen bonding (h) components. / / Thus, For selecting good solvents for PLGA, each solvent needs to be examined by individual contributions, i.e., δ d, δ p, and δ h separately, instead of the total solubility parameter δ t. If a solvent has high polar contribution, the solvent molecules may be attracted to each other than to PLGA. Thus, to avoid any inconsistencies in the like dissolves like principle, all 4 solubility parameters (δ d, δ p, δ h, and δ t ) need be measured for each PLGA.

28 Solubility Parameters of Solvents Solubility parameters of solvents can be determined from experimentally measuring H vap, and the values of most solvents are available in the literature. The H vap value of a solvent measured at one temperature can be converted to the values at any other temperature using the following empirical equation:,,. where T c is the critical temperature, i.e., the theta (θ) temperature, at which a polymer of infinite size mixes with the solvent completely. The solubility parameters of the solvent blends can be calculated by averaging the solubility parameter values of the individual solvents by their volume fractions (ϕ i ): This equation is useful when different solvents need to be mixed to generate a solvent with particular solubility parameters, and/or solvents in need are not readily available or too toxic to use. It may sound strange but the mixture of two bad solvents can result in a good solvent, assuming that the two mix well.

29 Solubility Parameters of Polymers The solubility parameters of polymers are difficult to obtain, because there is no measurable value of H vap or boiling point for polymers. Thus, indirect methods are used to obtain the values. The classical method of determining Hansen solubility parameters of a polymer involves testing of polymer solubility in different solvents with known partial solubility parameters (Hansen s three solubility parameters). The partial solubility parameters of solvents that interact with the given solvent (i.e., either dissolves the polymer or the polymer swells in the solvent) are plotted in 3 planes of the δ d, δ p and δ h axes to form circles of radius R, which is the radius of interaction sphere. The coordinates of the center of the spheres represent Hansen solubility parameters of polymers. This classical method has a shortcoming that good solvents may not be included in the interaction sphere or poor solvents may be included in the sphere. This shortcoming was overcome by calculating the maximum separation distance and using it as the diameter of the interaction sphere to find the coordinates of its middle as the Hansen parameters. The limitation of this method is that the position of the middle is not well defined, if there exist several equal maximum separation distances. In addition, this method treats good and limiting solvents of a polymer with equal weight, i.e., there is no mechanism of distinguishing the good from the limiting solvents. Limiting solvents are those that have limited interactions, e.g., the polymer may not dissolve but can swell in them.

30 Solubility Parameters of Polymers The more accurate radius of interaction sphere can be obtained by assigning different weights to good and limiting solvents of a polymer. The weights of the solvents can be easily determined by the intrinsic viscosities [η] which can be readily measured by a triple-detector GPC. The measured [η] is normalized to unity as follows.,,,, where the subscript i indicates individual solvents. Solvents that do not dissolve the polymer are given an intrinsic viscosity of zero. The good solvents result in higher [η], while the limiting solvents results in much lower [η]. The accuracy of the experimentally measured δ d,i, δ p,i and δ h,i increases as more solvents are tested. Once the coordinates of δ d, δ p and δ h are plotted as the center, the separating distances from the center to points corresponding to good and limiting solvents are plotted to calculate R i values. The largest value of R i is considered the radius of interaction sphere, R. 4,,, This improved method based on weighted intrinsic viscosities of a polymer in various solvents results in the interaction sphere including all good and limiting solvents.

31 Solubility Parameters of PLGA

32 Solubility Parameters of PLGA + = Fig. 3 Fig. 3. Comparison of Hansen s solubility parameters for various solvents: contributions of dispersion forces ( f d ), polar interactions (f p ), and hydrogen bonding (f h ). ( ) Good solvent; ( ) intermediately good solvent; ( ) intermediately poor solvent; and (x) poor solvent. (Yeo et al., J. Controlled Rel. 93: 6-73, 003) Dried PLGA (L:G=50:50, [η] = 0.58 dl/g) (5 mg) was added to glass vials containing 5 ml of the test organic solvent. The vials were agitated overnight at room temperature. Solubility of PLGA was judged by visual examination. Solvents were classified into four groups: good solvents that formed clear polymer solutions; intermediately good solvents that formed turbid polymer solutions upon heating; intermediately poor solvents that were marginally able to swell the polymer; and poor solvents in which the polymer remained intact.

33 The Group Contribution Methods The solubility parameters of polymers can be measured by the weighted intrinsic viscosity ([η]) method, but measurement of [η] has not been practical until recently when GPC-TDS became available. Thus, the solubility parameter of polymers has been calculated using the group contribution methods, which consider the energy (E) as an additive property for low molecular weight materials. / / / where h i is the contribution of ith group to the molar heat of vaporization. The solubility parameter can also be calculated using Small s molar attraction constants (F). where ρ is the density. From the general structure of PLGA, the F value of each group found in the literature can be used to calculate the Hildebrand solubility parameter of PLGA.

34 The Group Contribution Methods Group CH 3 CH CH C=O O F F: (cal. cm 3 ) / cm 3 /mol at 5 C The Units of the solubility parameters: / / / Although the solubility parameters of polymers can be estimated using the group contribution methods, the estimated values are likely to be different from the experimentally measured values. This is because the group contribution methods only consider the chemical structure without considering the molecular weight of the polymer, and a polymer as a whole may have different interactions with solvents unlike small groups do. Even though the polymer has the same structure regardless of the molecular weight, the reality is that the same polymer often has different solubility in a given solvent depending on the molecular weight. Thus, it is important to obtain the partial solubility parameters based on experimental values, such as weighted intrinsic viscosity values.

35 The Mark-Houwink Plot The dependence of intrinsic viscosity [η] of a polymer on its molecular weight M is given by the Mark-Houwink(-Sakurada) equation: or where K and α are two parameters that depend on the solvent, polymer, and temperature. It is common to plot the Mark-Houwink equation in a log-log graph to calculate the K and α values from the intercept and the slope. The slope is related to the shape of the polymer molecules and the polymer-solvent interactions. α = 0.5 for a polymer under θ conditions (i.e., an unperturbed random coil). α = 0.8 for a polymer in a good solvent, while α = for rod-like polymers. The slope α is related to the solubility parameters of the polymer and the solvent. log[η] Good Solvent log M Poor Solvent Mark-Houwink plots of a polymer in different solvents. The molecular weight (M) remains the same but the V h of the polymer changes in different solvents. The intrinsic viscosity [η] increases as the V h increases in good solvents. Thus, the solvent quality for each PLGA can be characterized by using the K and α values.

36 The Intrinsic Viscosity and Solubility Parameters ACE: acetone, H O: water, DCM: dichloromethane Solvents in which the polymer is insoluble were given an intrinsic viscosity of zero. The calculated Hildebrand and Hansen solubility parameters of PLGA are indicated by broken lines. The shaded line indicates best cubic fit (R 0.95). The 90/0% ACE/H O blend was excluded from the fit. Intrinsic viscosities and solubility parameters of PLGA (L:G=50:50) in various solvents. C.G. Madsen, A. Skov, S. Baldursdottir, T. Rades, L. Jorgensen,.J. Medlicott. Simple measurements for prediction of drug release from polymer matrices Solubility parameters and intrinsic viscosity. Eur. J. Pharm. Biopharm., 9 (05) -7.