Next Article in Journal
Locating the Sets of Exceptional Points in Dissipative Systems and the Self-Stability of Bicycles
Next Article in Special Issue
Use of the Complex Zeros of the Partition Function to Investigate the Critical Behavior of the Generalized Interacting Self-Avoiding Trail Model
Previous Article in Journal
On the Significance of the Quantum Mechanical Covariance Matrix
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Entropy Driven Phase Transition in Polymer Gels: Mean Field Theory

Physics Department, Cleveland State University, Cleveland, OH 44115, USA
Entropy 2018, 20(7), 501; https://doi.org/10.3390/e20070501
Submission received: 16 April 2018 / Revised: 21 June 2018 / Accepted: 27 June 2018 / Published: 30 June 2018

Abstract

:
We present a mean field model of a gel consisting of P polymers, each of length L and Nz polyfunctional monomers. Each polyfunctional monomer forms z covalent bonds with the 2P bifunctional monomers at the ends of the linear polymers. We find that the entropy dependence on the number of polyfunctional monomers exhibits an abrupt change at Nz = 2P/z due to the saturation of possible crosslinks. This non-analytical dependence of entropy on the number of polyfunctionals generates a first-order phase transition between two gel phases: one poor and the other rich in poly-functional molecules.

1. Motivation

A polymer gel like polyacrylamide changes the volumes by a large factor of ~1000 when a small quantity of solvent like acetone is added to the solution or when the temperature is varied slightly. Central to the understanding of this phase transition are the covalent crosslinks [1,2,3,4,5] between the linear chains. Polymer chains, such as in hydroxypropylcellulose (HPC) immersed in a water solution, form hydrogen bonds with the water. As the temperature, the pH, or some other external condition is varied, there is a change in the strength of the hydrogen bonds that results in the formation or destruction of aggregates of crosslinked polymer chains. Light scattering experiments [6,7] are used to study the influence of the amount of crosslinkers on the properties of the gel phases. However, the theoretical work to date [1,2,3,4,5] explained the swelling transition as the sole result of the competition between the elastic energy of the gel network and the Flory-Huggins entropy and energy of linear polymers in solution, without accounting for the crosslinkers. Here, we study the mean-field theory of the thermodynamics of such a system that explicitly includes polyfunctional molecules that crosslink with bifunctional molecules at the ends of linear polymers.
The main goal of this work is to highlight the role of the crosslinking entropy and of the saturation effect that occurs when the number of polyfunctional molecules increases to become comparable to the number of polymers. With this in mind, and attempting not to clutter the issue, we present here numerical results for the special case of dry gels, no solvent, and we ignore the network elastic energy [1,2,3,4,5]. We find that the saturation of crosslinks induces a non-analytical dependence of entropy on the number of polyfunctional molecules, and consequently a strong first-order transition occurs even for zero energy. In other words, the crosslinking entropy drives a first-order transition.
In future work we will include all energy contributions (Flory-Huggins, elastic) and study the interplay of energy and entropy to achieve a more complete understanding of the thermodynamics of polymer gels. A mean field theory ignores the range of interactions, and this important aspect needs to be addressed by other techniques, such as Monte Carlo simulations or renormalization group calculations.

2. Model

The fundamental quantity of statistical mechanics that determines the thermodynamics is the partition function Ξ , which is the sum of the un-normalized Boltzmann probability function over all of the possible configurations. The logarithm of the partition function yields the thermodynamic potential Ω. Thus, the fundamental equation that gives the thermodynamic potential as a function of temperature, chemical potentials, etc. is determined.
We consider a lattice of N sites and volume Nv that is occupied by Nz polyfunctionals, N2 bifunctionals supporting P linear chains each of length L, N2 = P(L + 1), and N0 solvent molecules. The total number of molecules is: N = Nz + N2 + N0. The mean field [8,9] grand canonical partition function is:
Ξ = N z , N 0 Γ ( N z , N 0 ) ( 1 N ) P L + B exp ( E ( N z , N 0 ) / T ) exp ( μ z N z T ) exp ( μ 0 N 0 T )
where Г counts the number of configurations for a given Nz and N0; B is the number of crosslinks formed between the Nz polyfunctionals and the 2P bifunctionals at the ends of the chains; μz is the chemical potential conjugated to the number of polyfunctionals, μ0 is the chemical potential conjugated to the number of solvent molecules, and we use units so that kB = 1. The mean field theory is exact for the equivalent neighbor lattice (also known as the Renyi-Erdos network). The equivalent-neighbor lattice has N vertices and N(N − 1)/2 edges connecting each pair of vertices. There is no distance dependence of the interactions and, as a result, the terms in the partition function of Equation (1) depend only on the number of molecules.
We assume that it is energetically favorable to establish as many crosslinks as possible. Hence: B = min(zNz, 2P). Each end of the P polymers can accept one crosslink, and the Nz polyfunctionals can each accept z bonds. The energy is expressed as a quadratic function of the three numbers of molecules: N0, Nz, and N2:
  E T = χ 0 N 0 ( 1 N 0 / N ) + χ z N z ( 1 N Z   N ) + χ 2 N 2 ( 1 N 2 / N ) +   χ 0 z N 0 N z   / N + χ 2 z N 2 N z   / N + χ 02 N 2 N 0   / N
An important contribution to the energy that needs to be included in a complete model of gels is the Flory-Rehner elastic energy [10] of the gel network. One can also include the concentration dependence of the various χ coefficients (see Hirotsu in [1,2,3,4,5]). In Equation (2) above, we only listed terms that are quadratic in densities.
The combinatorical factor Γ = exp(S) and S is the entropy. There are three contributions to the entropy S: (1) the entropy of mixing the Nz polyfunctionals with the N2 byfunctionals and with the N0 solvent molecules; (2) the Flory-Huggins entropy [8,9,10] associated with arranging P linear chains on the byfunctionals, N2 = P(L + 1); (3) the crosslinking entropy associated with establishing covariant bonds between the 2P ends of the chains and the Nz polyfunctionals. The latter entropy exhibits non-analytical behavior as a function of Nz at 2P/z, which triggers a phase transition.
In the thermodynamic limit (large N), the partition function = e N Ψ is approximated by the largest term [8,9] in the sum. Then, the grand-canonical thermodynamic potential Ω is:
Ω T = ln ( ) = max [ S ( P L + B ) ln ( N ) + E ( N z , N 0 ) + μ z N z + μ 0 N 0 T ]
Using the Euler equation of thermodynamics, E = TSpV + μzNz + μ0N0, one can relate the thermodynamic potential, Ω = ETSμzNzμ0N0, to the osmotic pressure:
Ω =   p N v
where v is the volume per molecule. The Euler equation is the direct consequence of the extensiveness of energy and entropy. Equation (3) implies:
p v T =   max [ Ψ ]
where
Ψ = s ^ + e ( n z , n 0 ) + μ z n z + μ 0 n 0 T , s = S / N , p = P / N , b = B / N , n z = N z / N , n 0 = N 0 / N , and s ^ = s ( p L + b ) ln ( N )

3. Entropy

The entropy is the sum of three physically distinct contributions: S = S1 + S2 + S3.
The entropy of mixing the Nz polyfunctionals, the N2 bifunctionals, and the solvent molecules N0:
e S 1 =   N ! N 0 ! N 2 ! N z !
Here, N = N0 + N2 + Nz.
The entropy that is associated with arranging P linear chains, each of length L, on the set of N2 = P(L + 1) bifunctionals is [8,10]:
e S 2 =   N 2 ! P ! 2 P
Previous models [1,2,3,4,5] have accounted only for the above two contributions to the entropy: Flory-Huggins entropy. The main contribution of this paper is the calculation of the entropy S3 that is associated with establishing covariant bonds between the 2P byfunctionals at the ends of the chains and the Nz polyfunctionals. One gets different formulas depending on relative sizes of 2P and zNz.
We assume that energetically it is favorable to generate as many covalent bonds as possible between the Nz crosslinkers (each accepts z bonds) and the 2P bifunctionals at the ends of the linear polymers (each accepts one bond). The number of configurations is obtained by using combinations C(N, K) = N!/K!(NK)! of N objects that are taken in groups of K objects.
If 2P > zNz, the number of configurations of crosslinking bonds is:
e S 3 =   C ( 2 P , z ) C ( 2 P z , z ) C ( 2 P ( N z 1 ) z , z )   = ( 2 P ) ! ( 2 P z N z ) ! ( z ! ) N z
Here, C(2P, z) counts the number of configurations of z bonds coming out of one polyfunctional and each bond connecting to an end of a linear polymer. The second term counts the same type of bond configurations from a second polyfunctional and bonds ending on z of the 2Pz free ends of polymers. This process continues until all polyfunctionals are saturated: there are Nz terms in the product.
If 2P < zNz, the number of configurations of crosslinking bonds is
e S 3 = C ( 2 P , z ) C ( 2 P z , z ) C ( 2 P ( 2 P z 1 ) z , z ) C ( N z , 2 P z ) =       = ( 2 P ) ! ( N z ) ! ( z ! ) 2 P z ( 2 P z ) ! ( N z     2 P z ) !
Here, the process of counting is similar. There are 2P/z factors counting the crosslinks configurations: from each polyfunctional, z bonds connect to 2P ends of linear polymers. The last factor in Equation (9) counts the number of groups of 2P/z polyfunctionals that can be formed out of the available Nz molecules.
The entropy per molecule s = S/N is obtained by using the equations above and the Stirling approximation. We find:
s ^ = n 0 ln ( n 0 ) n z ln ( n z ) ( 2 p z n z ) ln ( 2 p z n z ) + p ln ( p ) ( L ln ( 2 ) ) p ( z + ln ( z ! ) ) n z                                                                                             if     z n z < 2 p
s ^ = n 0 ln ( n 0 ) ( n z 2 p z ) ln ( n z 2 p z ) + z 2 z p ln ( p ) ( L ln ( 2 ) + 2 + 2 z ln ( 2 ( z 1 ) ! ) ) p                                                 if     z n z > 2 p

4. Results

We solved numerically the maximization Equation (4) for the special case of no solvent molecules n 0 = 0 and for χ 2 = χ 2 z = 0 . These parameter values are chosen to highlight the role of the crosslinking entropy in triggering the first-order transition, even for zero energy. This is a novel entropic effect that is not included in previous studies [1,2,3,4,5]. This model exhibits strong first-order transitions originating in the non-analytical behavior of the entropy as a function of densities at 2P = zNz. Its physical origin is the transition from satisfying all the byfunctionals if 2P < zNz to satisfying all of the polyfunctionals if zNz < 2P. The crossover condition zNz = 2P corresponds to a density of polyfunctionals of n z = 1 1 + z 2 ( L + 1 ) . Hence for large z, z >> 1, the crossover occurs for very low volume densities of polyfunctionals, while for smaller z~1 (such as for water), the crossover occurs at larger volume densities of crosslinkers. The numerical results below are for z = 3 and L= 3.
To demonstrate this, we first consider the model with zero energy: χ z = 0 , i.e., at very high temperatures. In Figure 1, we show the density of polyfunctionals nz as a function of the chemical potential μz/T. Phase transitions occur around μz/T~−0.5. Two first-order transitions are exhibited: from nz~0 to nz~0.9 to nz~1.
The latter transition occurs at μz/T ~ −0.52. To understand it, we show in Figure 2, Ψ vs. nz for μz/T = −0.53 (lower curve), −0.52 (middle curve), and −0.51 (upper curve). The maximum in Ψ occurs at nz ~ 0.89 and at 1 for μz/T ~ −0.52.
In Figure 3, we show the thermodynamic potential, i.e., max(Ψ), as a function of the chemical potential. The graph also includes the branches corresponding to the two formulas for the entropy (Equations (10) and (11)). This shows that the phase transitions are induced by the changeover in the entropy dependence on the polyfunctionals’ density.
Next, we consider χ z = 1 . We show in Figure 4 and Figure 5 the isotherms for T = 1. There is a phase transition from a phase with very few polyfunctionals to a phase with nz ~ 1 at μz ~ −0.575.
Figure 5 shows the thermodynamic potential vs. the chemical potential. It confirms that the first-order transition is triggered by the entropic crossover from nz < 2P/z to nz > 2P/z. The presence of the energy term χznz(1 − nz) strengthen this transition, as this energy is minimized at nz = 0 and 1. As a result, the intermediate phase that is apparent at zero energy (see Figure 1) does not show up here.

5. Conclusions

The importance of crosslinks that are formed between polymers and other chemicals in solutions was recognized [1,2,3,4,5] as central to understanding phase transitions (such as swelling) in polymeric systems. In this work, we have explicitly calculated the crosslinking entropy. We have shown that the entropy that is associated with the formation of crosslinking bonds between the ends of linear polymers and the polyfunctional molecules crosses over from one dependence for 2P < zNz (when 2P bonds are formed) to another one for 2P > zNz (when zNz bonds are formed). The crossover results in strong first-order transitions between a phase poor in crosslinkers and a phase with a high concentration of crosslinkers. The important result demonstrated in this work is that the crosslinking entropy that measures the number of possibilities of establishing bonds between the ends of the linear polymers and the polyfunctional molecules, drives a strong first order transition.
This model can be expanded to include the elasticity energy of the linear polymers, which was considered by Tanaka and Hirotsu [1,2,3,4,5]. We will also analyze the phase diagram in the five-dimensional space of interactions, three couplings, and two chemical potentials, which is quite similar to the three-component model [11]. Furthermore, there is experimental evidence that the polymers are polydisperse as opposed to monodisperse (as was assumed here). The influence of the size distribution of polymers on the phase transitions is expected [8,10] to be quite interesting and will be considered. The role of the density of crosslinks on the thermodynamic properties can be analyzed by varying the functionality parameter z. This is relevant to slightly crosslinked gels, z >> 1, that are known to be superabsorbent [12].
The mean-field theory is exact for the equivalent-neighbor lattice [8,9] where the interactions are infinitely ranged: they are independent of distance. The success of the Flory-Huggins [10] theory of polymer solutions may be due to its realizability (in Berker’s sense [13,14]), i.e., it is exact for the equivalent neighbor lattice model, though it is approximate for the three dimensional system. This ensures that the thermodynamics laws are satisfied [13,14]. To introduce a distance in the model, one considers energy terms such as Jijninj, where ni is one if a particle is present at location i or zero if it is absent, and the coupling J depends on the distance between i and j. The mean-field theory is generally reliable in predicting qualitatively correct phase diagrams. However, for realistic finite-range interactions, the critical exponents [15] are generally modified, and the type of transition [16]—first- or second-order—may differ. The role of the range of interactions can be studied by means of Monte-Carlo simulations [17,18] and renormalization group calculations [19].

Funding

This work was partially supported by the Institute of Advanced Studies of the University of Cergy Pontoise under the Paris Seine Initiative for Excellence (“Investissements d’Avenir” ANR-16-IDEX-0008).

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Tanaka, T. Collapse of Gels and the Critical Endpoint. Phys. Rev. Lett. 1978, 40, 820. [Google Scholar] [CrossRef]
  2. Tanaka, T.; Fillmore, D.; Sun, S.-T.; Nishio, I.; Swislow, G.; Shah, A. Phase Transitions in Ionic Gels. Phys. Rev. Lett. 1980, 45, 1636–1639. [Google Scholar] [CrossRef]
  3. Amyia, T.; Tanaka, T. Phase transitions in crosslinked gels of natural polymers. Macromolecules 1987, 20, 1162–1164. [Google Scholar] [CrossRef]
  4. Hirotsu, S. Critical points of the volume phase transition in N-isopropylacrylamide gels. J. Chem. Phys. 1988, 88, 427. [Google Scholar] [CrossRef]
  5. Hirotsu, S. Static and time-dependent properties of polymer gels around the volume phase transition. Phase Trans. 1994, 47, 183–240. [Google Scholar] [CrossRef]
  6. Streletzky, K.A.; McKenna, J.T.; Mohieddine, R. Spectral time moment analysis of microgel structure and dynamics. J. Polym. Sci. Part B Polym. Phys. 2008, 46, 771–781. [Google Scholar] [CrossRef]
  7. Streletzky, K.A.; McKenna, J.T. Spectral time moment analysis of microgel deswelling. Effect of the heating rate. J. Polym. Sci. Part B Polym. Phys. 2008, 46, 2792–2802. [Google Scholar] [CrossRef] [Green Version]
  8. Kaufman, M. Equilibrium polymerization on the equivalent-neighbor lattice. Phys. Rev. B 1989, 39, 6898–6906. [Google Scholar] [CrossRef] [Green Version]
  9. Kaufman, M.; Kardar, M. Short-range and infinite-range bond percolation. Phys. Rev. B 1984, 29, 5053–5059. [Google Scholar] [CrossRef]
  10. Flory, P.J. Principles of Polymer Chemistry; Cornell University Press: Ithaca, NY, USA, 1986. [Google Scholar]
  11. Kaufman, M.; Griffiths, R.B.; Yeomans, J.M.; Fisher, M.E. Three-component model and tricritical points: A renormalization-group study. Two dimensions. Phys. Rev. B 1981, 23, 3448–3459. [Google Scholar] [CrossRef]
  12. Philippova, O.E. Responsive Polymer Gels. Polym. Sci. Ser. C 2000, 42, 208–228. [Google Scholar]
  13. Berker, A.N.; Ostlund, S. Renormalisation-group calculations of finite systems: order parameter and specific heat for epitaxial ordering. J. Phys. C 1979, 12, 4961. [Google Scholar] [CrossRef]
  14. Kaufman, M.; Griffiths, R.B. Convexity of the free energy in some real-space renormalization-group approximations. Phys. Rev. B 1983, 28, 3864–3865. [Google Scholar] [CrossRef]
  15. Fisher, M.E. The renormalization group in the theory of critical behavior. Rev. Mod. Phys. 1974, 46, 597–616. [Google Scholar] [CrossRef]
  16. Wu, F.Y. The Potts model. Rev. Mod. Phys. 1982, 54, 235–268. [Google Scholar] [CrossRef]
  17. Kenkare, N.R.; Hall, C.K.; Khan, S.A. Theory and simulation of the swelling of polymer gels. J. Chem. Phys. 2000, 113, 404–418. [Google Scholar] [CrossRef]
  18. Diep, H.T. Statistical Physics—Fundamentals and Application to Condensed Matter; Word Scientific: Singapore, 2015. [Google Scholar]
  19. Kaufman, M.; Berger, T.; Gujrati, P.D.; Bowman, D. Polymerization on the diamond hierarchical lattice: The Migdal-Kadanoff renormalization-group scheme. Phys. Rev. A 1990, 41, 4371–4378. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. The density of polyfunctionals nz versus the chemical potential μz/T for zero energy.
Figure 1. The density of polyfunctionals nz versus the chemical potential μz/T for zero energy.
Entropy 20 00501 g001
Figure 2. Ψ vs. nz. There is a phase transition between nz~0.89 and nz = 1 for μz/T = −0.52 (middle curve). For μz/T = −0.53 (lower curve), the maximum of Ψ is at nz = 0.89. For μz/T = −0.51 (upper curve), the maximum in Ψ occurs at nz = 1.
Figure 2. Ψ vs. nz. There is a phase transition between nz~0.89 and nz = 1 for μz/T = −0.52 (middle curve). For μz/T = −0.53 (lower curve), the maximum of Ψ is at nz = 0.89. For μz/T = −0.51 (upper curve), the maximum in Ψ occurs at nz = 1.
Entropy 20 00501 g002
Figure 3. pv/T vs. μz/T for zero energy; the entropic crossover between znz < 2p and znz > 2p is responsible for the transition. The two lines correspond to the entropy formula, Equations (10) and (11). In the thermodynamic limit, the maximization procedure, Equation (4), choses for a given chemical potential the higher osmotic pressure, heavy curve.
Figure 3. pv/T vs. μz/T for zero energy; the entropic crossover between znz < 2p and znz > 2p is responsible for the transition. The two lines correspond to the entropy formula, Equations (10) and (11). In the thermodynamic limit, the maximization procedure, Equation (4), choses for a given chemical potential the higher osmotic pressure, heavy curve.
Entropy 20 00501 g003
Figure 4. Density of polyfunctionals nz versus chemical potential μz for χz = 1 and T = 1.
Figure 4. Density of polyfunctionals nz versus chemical potential μz for χz = 1 and T = 1.
Entropy 20 00501 g004
Figure 5. pv/T vs. μz for χz = 1 and T = 1. The two lines correspond to the entropy formula, Equations (10) and (11). In the thermodynamic limit, the maximization procedure, Equation (4), choses for a given chemical potential the higher osmotic pressure, heavy curve.
Figure 5. pv/T vs. μz for χz = 1 and T = 1. The two lines correspond to the entropy formula, Equations (10) and (11). In the thermodynamic limit, the maximization procedure, Equation (4), choses for a given chemical potential the higher osmotic pressure, heavy curve.
Entropy 20 00501 g005

Share and Cite

MDPI and ACS Style

Kaufman, M. Entropy Driven Phase Transition in Polymer Gels: Mean Field Theory. Entropy 2018, 20, 501. https://doi.org/10.3390/e20070501

AMA Style

Kaufman M. Entropy Driven Phase Transition in Polymer Gels: Mean Field Theory. Entropy. 2018; 20(7):501. https://doi.org/10.3390/e20070501

Chicago/Turabian Style

Kaufman, Miron. 2018. "Entropy Driven Phase Transition in Polymer Gels: Mean Field Theory" Entropy 20, no. 7: 501. https://doi.org/10.3390/e20070501

APA Style

Kaufman, M. (2018). Entropy Driven Phase Transition in Polymer Gels: Mean Field Theory. Entropy, 20(7), 501. https://doi.org/10.3390/e20070501

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop