Flory–Huggins solution theory

<h2 id="theory">Theory</h2>
The <a href="/facts/Thermodynamic_potentials/rBnf2hB6">thermodynamic equation</a> for the <a href="/facts/Gibbs_energy/bvenDuzh">Gibbs energy</a> change accompanying mixing at constant <a href="/facts/Temperature/QyGe4gmj">temperature</a> and (external) <a href="/facts/Pressure/eyc37GRG">pressure</a> is

Δ
        
          G
          
            
              m
              i
              x
            
          
        
        =
        Δ
        
          H
          
            
              m
              i
              x
            
          
        
        −
        T
        Δ
        
          S
          
            
              m
              i
              x
            
          
        
      
    
    {\displaystyle \Delta G_{\rm {mix}}=\Delta H_{\rm {mix}}-T\Delta S_{\rm {mix}}}

A change, denoted by 
 
 
 
 Δ
 
 
 {\displaystyle \Delta }
 
, is the <a href="/facts/Number/auQSAE2m">value</a> of a <a href="/facts/Variable_(mathematics)/VbQlrvKz">variable</a> for a <a href="/facts/Solution_(chemistry)/MGchhMQV">solution</a> or <a href="/facts/Mixture/tms8rbmn">mixture</a> minus the values for the pure <a href="/facts/Component_(thermodynamics)/XaLQE35f">components</a> considered separately. The objective is to find explicit <a href="/facts/Formula/BosTcKnv">formulas</a> for 
 
 
 
 Δ
 
 H
 
 
 m
 i
 x
 
 
 
 
 
 {\displaystyle \Delta H_{\rm {mix}}}
 
 and 
 
 
 
 Δ
 
 S
 
 
 m
 i
 x
 
 
 
 
 
 {\displaystyle \Delta S_{\rm {mix}}}
 
, the <a href="/facts/Enthalpy/bBou8lCE">enthalpy</a> and <a href="/facts/Entropy/s52IXeFz">entropy</a> increments associated with the mixing <a href="/facts/Process_(science)/182j8RWE">process</a>.
The result obtained by <a href="/facts/Paul_Flory/p5n6bFca">Flory</a><a href="https://en.wikipedia.org/wiki/Flory%E2%80%93Huggins_solution_theory#endnote_1">[1]</a> and <a href="/facts/Maurice_Loyal_Huggins/KDjfm6J9">Huggins</a><a href="https://en.wikipedia.org/wiki/Flory%E2%80%93Huggins_solution_theory#endnote_2">[2]</a> is

Δ
        
          G
          
            
              m
              i
              x
            
          
        
        =
        R
        T
        [
        
        
          n
          
            1
          
        
        ln
        ⁡
        
          ϕ
          
            1
          
        
        +
        
          n
          
            2
          
        
        ln
        ⁡
        
          ϕ
          
            2
          
        
        +
        
          n
          
            1
          
        
        
          ϕ
          
            2
          
        
        
          χ
          
            12
          
        
        
        ]
      
    
    {\displaystyle \Delta G_{\rm {mix}}=RT[\,n_{1}\ln \phi _{1}+n_{2}\ln \phi _{2}+n_{1}\phi _{2}\chi _{12}\,]}

The right-hand side is a <a href="/facts/Function_(mathematics)/CdReEKCh">function</a> of the number of <a href="/facts/Mole_(unit)/kCMFwCc0">moles</a> 
 
 
 
 
 n
 
 1
 
 
 
 
 {\displaystyle n_{1}}
 
 and volume fraction 
 
 
 
 
 ϕ
 
 1
 
 
 
 
 {\displaystyle \phi _{1}}
 
 of <a href="/facts/Solvent/LoMPMz8S">solvent</a> (<a href="/facts/Component_(thermodynamics)/XaLQE35f">component</a> 
 
 
 
 1
 
 
 {\displaystyle 1}
 
), the number of moles 
 
 
 
 
 n
 
 2
 
 
 
 
 {\displaystyle n_{2}}
 
 and volume fraction 
 
 
 
 
 ϕ
 
 2
 
 
 
 
 {\displaystyle \phi _{2}}
 
 of polymer (component 
 
 
 
 2
 
 
 {\displaystyle 2}
 
), with the introduction of a parameter 
 
 
 
 χ
 
 
 {\displaystyle \chi }
 
 to take account of the <a href="/facts/Energy/wd8PrQM7">energy</a> of interdispersing polymer and solvent molecules. 
 
 
 
 R
 
 
 {\displaystyle R}
 
 is the <a href="/facts/Gas_constant/RUXszrMH">gas constant</a> and 
 
 
 
 T
 
 
 {\displaystyle T}
 
 is the <a href="/facts/Thermodynamic_temperature/JXeJJPa6">absolute temperature</a>. The volume fraction is analogous to the <a href="/facts/Mole_fraction/nW4EepmJ">mole fraction</a>, but is weighted to take account of the relative sizes of the molecules. For a small solute, the mole fractions would appear instead, and this modification is the innovation due to Flory and Huggins. In the most general case the mixing parameter, 
 
 
 
 χ
 
 
 {\displaystyle \chi }
 
, is a free energy parameter, thus including an entropic component.<a class="footnote-ref" id="fnref:1" href="#fn:1">1</a><a class="footnote-ref" id="fnref:2" href="#fn:2">2</a>

<h2 id="derivation">Derivation</h2>
We first calculate the <a href="/facts/Entropy_of_mixing/zDdJNHTe">entropy of mixing</a>, the increase in the <a href="/facts/Entropy/s52IXeFz">uncertainty</a> about the locations of the molecules when they are interspersed. In the pure condensed <a href="/facts/Phase_(matter)/NkHXDdN3">phases</a> – <a href="/facts/Solvent/LoMPMz8S">solvent</a> and polymer – everywhere we look we find a molecule.<a href="https://en.wikipedia.org/wiki/Flory%E2%80%93Huggins_solution_theory#endnote_3">[3]</a> Of course, any notion of "finding" a molecule in a given location is a <a href="/facts/Thought_experiment/GXlr0lz9">thought experiment</a> since we can't actually examine <a href="/facts/Space/tq6LUM3C">spatial</a> locations the size of molecules. The <a href="/facts/Expression_(mathematics)/MPrMlYbE">expression</a> for the <a href="/facts/Entropy_of_mixing/zDdJNHTe">entropy of mixing</a> of small molecules in terms of <a href="/facts/Mole_fraction/nW4EepmJ">mole fractions</a> is no longer reasonable when the <a href="/facts/Solute/MGchhMQV">solute</a> is a <a href="/facts/Macromolecule/039ULSLg">macromolecular</a> <a href="/facts/Ideal_chain/JXqaeaPC">chain</a>. We take account of this dis<a href="/facts/Symmetry/hpu7DK8p">symmetry</a> in molecular sizes by assuming that individual polymer segments and individual solvent molecules occupy sites on a <a href="/facts/Lattice_(group)/2ls9Iv9N">lattice</a>. Each site is occupied by exactly one molecule of the solvent or by one <a href="/facts/Monomer/Gw3CltMb">monomer</a> of the polymer chain, so the total number of sites is

N
        =
        
          N
          
            1
          
        
        +
        x
        
          N
          
            2
          
        
      
    
    {\displaystyle N=N_{1}+xN_{2}}

where 
 
 
 
 
 N
 
 1
 
 
 
 
 {\displaystyle N_{1}}
 
 is the number of solvent molecules and 
 
 
 
 
 N
 
 2
 
 
 
 
 {\displaystyle N_{2}}
 
 is the number of polymer molecules, each of which has 
 
 
 
 x
 
 
 {\displaystyle x}
 
 segments.<a href="https://en.wikipedia.org/wiki/Flory%E2%80%93Huggins_solution_theory#endnote_4">[4]</a>
For a <a href="/facts/Random_walk/08v0jmfv">random walk</a> on a lattice<a class="footnote-ref" id="fnref:3" href="#fn:3">3</a> we can calculate the <a href="/facts/Entropy/s52IXeFz">entropy</a> change (the increase in <a href="/facts/Space/tq6LUM3C">spatial</a> <a href="/facts/Information_entropy/NLg4NLvt">uncertainty</a>) as a result of mixing solute and solvent.

Δ
        
          S
          
            
              m
              i
              x
            
          
        
        =
        −
        
          k
          
            
              B
            
          
        
        
          [
          
            
              N
              
                1
              
            
            ln
            ⁡
            
              
                
                  
                    N
                    
                      1
                    
                  
                  N
                
              
            
            +
            
              N
              
                2
              
            
            ln
            ⁡
            
              
                
                  
                    x
                    
                      N
                      
                        2
                      
                    
                  
                  N
                
              
            
          
          ]
        
      
    
    {\displaystyle \Delta S_{\rm {mix}}=-k_{\rm {B}}\left[N_{1}\ln {\tfrac {N_{1}}{N}}+N_{2}\ln {\tfrac {xN_{2}}{N}}\right]}

where 
 
 
 
 
 k
 
 
 B
 
 
 
 
 
 {\displaystyle k_{\rm {B}}}
 
 is the <a href="/facts/Boltzmann_constant/5DIejyR3">Boltzmann constant</a>. Define the lattice volume fractions 
 
 
 
 
 ϕ
 
 1
 
 
 
 
 {\displaystyle \phi _{1}}
 
 and 
 
 
 
 
 ϕ
 
 2
 
 
 
 
 {\displaystyle \phi _{2}}

ϕ
          
            1
          
        
        =
        
          
            
              N
              
                1
              
            
            N
          
        
        ,
        
        
          ϕ
          
            2
          
        
        =
        
          
            
              x
              
                N
                
                  2
                
              
            
            N
          
        
      
    
    {\displaystyle \phi _{1}={\frac {N_{1}}{N}},\quad \phi _{2}={\frac {xN_{2}}{N}}}

These are also the probabilities that a given lattice site, chosen at <a href="/facts/Randomness/QJQBoEX5">random</a>, is occupied by a solvent molecule or a polymer segment, respectively. Thus

Δ
        
          S
          
            
              m
              i
              x
            
          
        
        =
        −
        
          k
          
            
              B
            
          
        
        [
        
        
          N
          
            1
          
        
        ln
        ⁡
        
          ϕ
          
            1
          
        
        +
        
          N
          
            2
          
        
        ln
        ⁡
        
          ϕ
          
            2
          
        
        
        ]
      
    
    {\displaystyle \Delta S_{\rm {mix}}=-k_{\rm {B}}[\,N_{1}\ln \phi _{1}+N_{2}\ln \phi _{2}\,]}

For a small solute whose molecules occupy just one lattice site, 
 
 
 
 x
 
 
 {\displaystyle x}
 
 equals one, the volume fractions reduce to <a href="/facts/Mole_fraction/nW4EepmJ">molecular or mole fractions</a>, and we recover the usual <a href="/facts/Entropy_of_mixing/zDdJNHTe">entropy of mixing</a>.
In addition to the entropic effect, we can expect an enthalpy change.<a href="https://en.wikipedia.org/wiki/Flory%E2%80%93Huggins_solution_theory#endnote_5">[5]</a> There are three molecular interactions to consider: solvent-solvent 
 
 
 
 
 w
 
 11
 
 
 
 
 {\displaystyle w_{11}}
 
, monomer-monomer 
 
 
 
 
 w
 
 22
 
 
 
 
 {\displaystyle w_{22}}
 
 (not the <a href="/facts/Covalent_bond/EzMLf3Pz">covalent bonding</a>, but between different chain sections), and monomer-solvent 
 
 
 
 
 w
 
 12
 
 
 
 
 {\displaystyle w_{12}}
 
. Each of the last occurs at the expense of the average of the other two, so the energy increment per monomer-solvent contact is

Δ
        w
        =
        
          w
          
            12
          
        
        −
        
          
            
              1
              2
            
          
        
        (
        
          w
          
            22
          
        
        +
        
          w
          
            11
          
        
        )
      
    
    {\displaystyle \Delta w=w_{12}-{\tfrac {1}{2}}(w_{22}+w_{11})}

The total number of such contacts is

x
        
          N
          
            2
          
        
        z
        
          ϕ
          
            1
          
        
        =
        
          N
          
            1
          
        
        
          ϕ
          
            2
          
        
        z
      
    
    {\displaystyle xN_{2}z\phi _{1}=N_{1}\phi _{2}z}

where 
 
 
 
 z
 
 
 {\displaystyle z}
 
 is the coordination number, the number of nearest neighbors for a lattice site, each one occupied either by one chain segment or a solvent molecule. That is, 
 
 
 
 x
 
 N
 
 2
 
 
 
 
 {\displaystyle xN_{2}}
 
 is the total number of polymer segments (monomers) in the solution, so 
 
 
 
 x
 
 N
 
 2
 
 
 z
 
 
 {\displaystyle xN_{2}z}
 
 is the number of nearest-neighbor sites to all the polymer segments. Multiplying by the probability 
 
 
 
 
 ϕ
 
 1
 
 
 
 
 {\displaystyle \phi _{1}}
 
 that any such site is occupied by a solvent molecule,<a href="https://en.wikipedia.org/wiki/Flory%E2%80%93Huggins_solution_theory#endnote_6">[6]</a> we obtain the total number of polymer-solvent molecular interactions. An approximation following <a href="/facts/Mean_field_theory/tJrwSURE">mean field theory</a> is made by following this procedure, thereby reducing the complex problem of many interactions to a simpler problem of one interaction.
The enthalpy change is equal to the energy change per polymer monomer-solvent interaction multiplied by the number of such interactions

Δ
        
          H
          
            
              m
              i
              x
            
          
        
        =
        
          N
          
            1
          
        
        
          ϕ
          
            2
          
        
        z
        Δ
        w
      
    
    {\displaystyle \Delta H_{\rm {mix}}=N_{1}\phi _{2}z\Delta w}

The polymer-solvent interaction parameter <a href="/facts/Chi_(letter)/pj4UT8Pf">chi</a> is defined as

χ
          
            12
          
        
        =
        
          
            
              z
              Δ
              w
            
            
              
                k
                
                  
                    B
                  
                
              
              T
            
          
        
      
    
    {\displaystyle \chi _{12}={\frac {z\Delta w}{k_{\rm {B}}T}}}

It depends on the nature of both the solvent and the solute, and is the only material-specific parameter in the model. The enthalpy change becomes

Δ
        
          H
          
            
              m
              i
              x
            
          
        
        =
        
          k
          
            
              B
            
          
        
        T
        
          N
          
            1
          
        
        
          ϕ
          
            2
          
        
        
          χ
          
            12
          
        
      
    
    {\displaystyle \Delta H_{\rm {mix}}=k_{\rm {B}}TN_{1}\phi _{2}\chi _{12}}

Assembling terms, the total free energy change is

where we have converted the expression from molecules 
 
 
 
 
 N
 
 1
 
 
 
 
 {\displaystyle N_{1}}
 
 and 
 
 
 
 
 N
 
 2
 
 
 
 
 {\displaystyle N_{2}}
 
 to moles 
 
 
 
 
 n
 
 1
 
 
 
 
 {\displaystyle n_{1}}
 
 and 
 
 
 
 
 n
 
 2
 
 
 
 
 {\displaystyle n_{2}}
 
 by transferring the <a href="/facts/Avogadro_constant/SwajCq22">Avogadro constant</a> 
 
 
 
 
 N
 
 A
 
 
 
 
 {\displaystyle N_{\text{A}}}
 
 to the <a href="/facts/Gas_constant/RUXszrMH">gas constant</a> 
 
 
 
 R
 =
 
 k
 
 
 B
 
 
 
 
 N
 
 A
 
 
 
 
 {\displaystyle R=k_{\rm {B}}N_{\text{A}}}
 
.
The value of the interaction parameter can be estimated from the <a href="/facts/Hildebrand_solubility_parameter/U8jmfFGH">Hildebrand solubility parameters</a> 
 
 
 
 
 δ
 
 a
 
 
 
 
 {\displaystyle \delta _{a}}
 
 and 
 
 
 
 
 δ
 
 b
 
 
 
 
 {\displaystyle \delta _{b}}

χ
          
            12
          
        
        =
        
          
            
              
                V
                
                  
                    s
                    e
                    g
                  
                
              
              (
              
                δ
                
                  a
                
              
              −
              
                δ
                
                  b
                
              
              
                )
                
                  2
                
              
            
            
              R
              T
            
          
        
      
    
    {\displaystyle \chi _{12}={\frac {V_{\rm {seg}}(\delta _{a}-\delta _{b})^{2}}{RT}}}

where 
 
 
 
 
 V
 
 
 s
 e
 g
 
 
 
 
 
 {\displaystyle V_{\rm {seg}}}
 
 is the actual volume of a polymer segment.
In the most general case the interaction 
 
 
 
 Δ
 w
 
 
 {\displaystyle \Delta w}
 
 and the ensuing mixing parameter, 
 
 
 
 χ
 
 
 {\displaystyle \chi }
 
, is a free energy parameter, thus including an entropic component.<a class="footnote-ref" id="fnref:4" href="#fn:4">4</a><a class="footnote-ref" id="fnref:5" href="#fn:5">5</a> This means that aside to the regular mixing entropy there is another entropic contribution from the interaction between solvent and monomer. This contribution is sometimes very important in order to make quantitative predictions of thermodynamic properties.
More advanced solution theories exist, such as the Flory–Krigbaum theory.

<h2 id="liquid-liquid-phase-separation">Liquid-liquid phase separation</h2>

Polymers can separate out from the solvent, and do so in a characteristic way.<a class="footnote-ref" id="fnref:6" href="#fn:6">6</a> The Flory–Huggins free energy per unit volume, for a polymer with 
 
 
 
 N
 
 
 {\displaystyle N}
 
 monomers, can be written in a simple dimensionless form

f
        =
        
          
            ϕ
            N
          
        
        ln
        ⁡
        ϕ
        +
        (
        1
        −
        ϕ
        )
        ln
        ⁡
        (
        1
        −
        ϕ
        )
        +
        χ
        ϕ
        (
        1
        −
        ϕ
        )
      
    
    {\displaystyle f={\frac {\phi }{N}}\ln \phi +(1-\phi )\ln(1-\phi )+\chi \phi (1-\phi )}

for 
 
 
 
 ϕ
 
 
 {\displaystyle \phi }
 
 the volume fraction of monomers, and 
 
 
 
 N
 ≫
 1
 
 
 {\displaystyle N\gg 1}
 
. The <a href="/facts/Osmotic_pressure/6oWKTsXV">osmotic pressure</a> (in reduced units) is

Π
 =
 
 
 ϕ
 N
 
 
 −
 ln
 ⁡
 (
 1
 −
 ϕ
 )
 −
 ϕ
 −
 χ
 
 ϕ
 
 2
 
 
 
 
 {\displaystyle \Pi ={\frac {\phi }{N}}-\ln(1-\phi )-\phi -\chi \phi ^{2}}
 
.
The polymer solution is stable with respect to small fluctuations when the second derivative of this free energy is positive. This second derivative is

f
          ″
        
        =
        
          
            1
            
              N
              ϕ
            
          
        
        +
        
          
            1
            
              1
              −
              ϕ
            
          
        
        −
        2
        χ
      
    
    {\displaystyle f''={\frac {1}{N\phi }}+{\frac {1}{1-\phi }}-2\chi }

and the solution first becomes unstable when this and the third derivative

f
          ‴
        
        =
        −
        
          
            1
            
              N
              
                ϕ
                
                  2
                
              
            
          
        
        +
        
          
            1
            
              (
              1
              −
              ϕ
              
                )
                
                  2
                
              
            
          
        
      
    
    {\displaystyle f'''=-{\frac {1}{N\phi ^{2}}}+{\frac {1}{(1-\phi )^{2}}}}

are both equal to zero. A little algebra then shows that the polymer solution first becomes unstable at a critical point at

χ
          
            cp
          
        
        ≃
        1
        
          /
        
        2
        +
        
          N
          
            −
            1
            
              /
            
            2
          
        
        +
        ⋯
        
        
          ϕ
          
            cp
          
        
        ≃
        
          N
          
            −
            1
            
              /
            
            2
          
        
        −
        
          N
          
            −
            1
          
        
        +
        ⋯
      
    
    {\displaystyle \chi _{\text{cp}}\simeq 1/2+N^{-1/2}+\cdots \qquad \phi _{\text{cp}}\simeq N^{-1/2}-N^{-1}+\cdots }

This means that for all values of 
 
 
 
 0
 <
 χ
 ≲
 1
 
 /
 
 2
 
 
 {\displaystyle 0<\chi \lesssim 1/2}
 
 the monomer-solvent effective interaction is weakly repulsive, but this is too weak to cause liquid/liquid separation. However, when 
 
 
 
 χ
 >
 1
 
 /
 
 2
 
 
 {\displaystyle \chi >1/2}
 
, there is separation into two coexisting phases, one richer in polymer but poorer in solvent, than the other.
The unusual feature of the liquid/liquid phase separation is that it is highly asymmetric: the volume fraction of monomers at the critical point is approximately 
 
 
 
 
 N
 
 −
 1
 
 /
 
 2
 
 
 
 
 {\displaystyle N^{-1/2}}
 
, which is very small for large polymers. The amount of polymer in the solvent-rich/polymer-poor coexisting phase is extremely small for long polymers. The solvent-rich phase is close to pure solvent. This is peculiar to polymers, a mixture of small molecules can be approximated using the Flory–Huggins expression with 
 
 
 
 N
 =
 1
 
 
 {\displaystyle N=1}
 
, and then 
 
 
 
 
 ϕ
 
 cp
 
 
 =
 1
 
 /
 
 2
 
 
 {\displaystyle \phi _{\text{cp}}=1/2}
 
 and both coexisting phases are far from pure.

<h2 id="polymer-blends">Polymer blends</h2>
Synthetic polymers rarely consist of chains of uniform length in solvent. The Flory–Huggins free energy density can be generalized<a class="footnote-ref" id="fnref:7" href="#fn:7">7</a> to an N-component mixture of polymers with lengths 
 
 
 
 
 r
 
 i
 
 
 
 
 {\displaystyle r_{i}}
 
 by

f
        
          
            (
          
        
        {
        
          ϕ
          
            i
          
        
        ,
        
          r
          
            i
          
        
        }
        
          
            )
          
        
        =
        
          ∑
          
            i
            =
            1
          
          
            N
          
        
        
          
            
              ϕ
              
                i
              
            
            
              r
              
                i
              
            
          
        
        ln
        ⁡
        
          ϕ
          
            i
          
        
        +
        
          
            1
            2
          
        
        
          ∑
          
            i
            ,
            j
            =
            1
          
          
            N
          
        
        
          ϕ
          
            i
          
        
        
          ϕ
          
            j
          
        
        
          χ
          
            i
            j
          
        
      
    
    {\displaystyle f{\Bigl (}\{\phi _{i},r_{i}\}{\Bigr )}=\sum _{i=1}^{N}{\frac {\phi _{i}}{r_{i}}}\ln \phi _{i}+{\frac {1}{2}}\sum _{i,j=1}^{N}\phi _{i}\phi _{j}\chi _{ij}}

For a binary <a href="/facts/Polymer_blend/yoNz73uk">polymer blend</a>, where one species consists of 
 
 
 
 
 N
 
 A
 
 
 
 
 {\displaystyle N_{A}}
 
 monomers and the other 
 
 
 
 
 N
 
 B
 
 
 
 
 {\displaystyle N_{B}}
 
 monomers this simplifies to

f
        (
        ϕ
        )
        =
        
          
            ϕ
            
              N
              
                A
              
            
          
        
        ln
        ⁡
        ϕ
        +
        
          
            
              1
              −
              ϕ
            
            
              N
              
                B
              
            
          
        
        ln
        ⁡
        (
        1
        −
        ϕ
        )
        +
        χ
        ϕ
        (
        1
        −
        ϕ
        )
      
    
    {\displaystyle f(\phi )={\frac {\phi }{N_{A}}}\ln \phi +{\frac {1-\phi }{N_{B}}}\ln(1-\phi )+\chi \phi (1-\phi )}

As in the case for dilute polymer solutions, the first two terms on the right-hand side represent the entropy of mixing. For large polymers of 
 
 
 
 
 N
 
 A
 
 
 ≫
 1
 
 
 {\displaystyle N_{A}\gg 1}
 
 and 
 
 
 
 
 N
 
 B
 
 
 ≫
 1
 
 
 {\displaystyle N_{B}\gg 1}
 
 these terms are negligibly small. This implies that for a stable mixture to exist 
 
 
 
 χ
 <
 0
 
 
 {\displaystyle \chi <0}
 
, so for polymers A and B to blend their segments must attract one another.<a class="footnote-ref" id="fnref:8" href="#fn:8">8</a>

<h2 id="limitations">Limitations</h2>
Flory–Huggins theory tends to agree well with experiments in the semi-dilute concentration regime and can be used to fit data for even more complicated blends with higher concentrations. The theory qualitatively predicts phase separation, the tendency for high molecular weight species to be immiscible, the 
 
 
 
 χ
 ∝
 
 T
 
 −
 1
 
 
 
 
 {\displaystyle \chi \propto T^{-1}}
 
 interaction-temperature dependence and other features commonly observed in polymer mixtures. However, unmodified Flory–Huggins theory fails to predict the <a href="/facts/Lower_critical_solution_temperature/DJQmzsIq">lower critical solution temperature</a> observed in some polymer blends and the lack of dependence of the critical temperature 
 
 
 
 
 T
 
 c
 
 
 
 
 {\displaystyle T_{\text{c}}}
 
 on chain length 
 
 
 
 
 r
 
 i
 
 
 
 
 {\displaystyle r_{i}}
 
.<a class="footnote-ref" id="fnref:9" href="#fn:9">9</a> Additionally, it can be shown that for a binary blend of polymer species with equal chain lengths 
 
 
 
 (
 
 N
 
 A
 
 
 =
 
 N
 
 B
 
 
 )
 
 
 {\displaystyle (N_{A}=N_{B})}
 
 the critical concentration should be 
 
 
 
 
 ψ
 
 c
 
 
 =
 1
 
 /
 
 2
 
 
 {\displaystyle \psi _{\text{c}}=1/2}
 
; however, polymers blends have been observed where this parameter is highly asymmetric. In certain blends, mixing entropy can dominate over monomer interaction. By adopting the mean-field approximation, 
 
 
 
 χ
 
 
 {\displaystyle \chi }
 
 parameter complex dependence on <a href="/facts/Temperature-responsive_polymer/3FUmGRYW">temperature</a>, blend composition, and chain length was discarded. Specifically, interactions beyond the nearest neighbor may be highly relevant to the behavior of the blend and the distribution of polymer segments is not necessarily uniform, so certain lattice sites may experience interaction energies disparate from that approximated by the mean-field theory.
One well-studied<a class="footnote-ref" id="fnref:10" href="#fn:10">10</a><a class="footnote-ref" id="fnref:11" href="#fn:11">11</a> effect on interaction energies neglected by unmodified Flory–Huggins theory is chain correlation. In dilute polymer mixtures, where chains are well separated, intramolecular forces between monomers of the polymer chain dominate and drive demixing leading to regions where polymer concentration is high. As the polymer concentration increases, chains tend to overlap and the effect becomes less important. In fact, the demarcation between dilute and semi-dilute solutions is commonly defined by the concentration where polymers begin to overlap 
 
 
 
 
 c
 
 ∗
 
 
 
 
 {\displaystyle c^{*}}
 
 which can be estimated as

c
          
            ∗
          
        
        =
        
          
            m
            
              
                
                  4
                  3
                
              
              π
              
                R
                
                  g
                
                
                  3
                
              
            
          
        
      
    
    {\displaystyle c^{*}={\frac {m}{{\frac {4}{3}}\pi R_{\text{g}}^{3}}}}

Here, m is the mass of a single polymer chain, and 
 
 
 
 
 R
 
 g
 
 
 
 
 {\displaystyle R_{\text{g}}}
 
 is the chain's <a href="/facts/Radius_of_gyration/VWdacrjX">radius of gyration</a>.

<h2 id="footnotes">Footnotes</h2>

<ol><li>^ "<a href="/facts/Thermodynamics/q4hIx3yP">Thermodynamics</a> of High <a href="/facts/Polymer/eF32E50K">Polymer</a> <a href="/facts/Solution_(chemistry)/MGchhMQV">Solutions</a>", <a href="/facts/Paul_Flory/p5n6bFca">Paul J. Flory</a> Journal of Chemical Physics, August 1941, Volume 9, Issue 8, p. 660 <a href="https://aip.scitation.org/doi/10.1063/1.1750971">Abstract</a>. Flory suggested that Huggins' name ought to be first since he had published several months earlier: Flory, P.J., "Thermodynamics of high polymer solutions", J. Chem. Phys. 10:51-61 (1942) <a href="http://www.garfield.library.upenn.edu/classics1985/A1985AFW3100001.pdf">Citation Classic No. 18, May 6, 1985</a></li>
<li>^ "Solutions of Long Chain <a href="/facts/Chemical_compound/37FXqqtv">Compounds</a>", <a href="/facts/Maurice_Loyal_Huggins/KDjfm6J9">Maurice L. Huggins</a> Journal of Chemical Physics, May 1941 Volume 9, Issue 5, p. 440 <a href="https://archive.today/20130223092410/http://link.aip.org/getabs/servlet/GetabsServlet?prog=normal&id=JCPSA6000009000005000440000001&idtype=cvips&gifs=yes">Abstract</a></li>
<li>^ We are ignoring the free volume due to molecular disorder in liquids and amorphous solids as compared to <a href="/facts/Crystal/e2R2mk3E">crystals</a>. This, and the assumption that <a href="/facts/Monomer/Gw3CltMb">monomers</a> and solute molecules are really the same size, are the main geometric approximations in this model.</li>
<li>^ For a real <a href="/facts/Chemical_synthesis/3KXsVVwC">synthetic</a> polymer, there is a <a href="/facts/Statistics/DNPsKGYU">statistical</a> <a href="/facts/Random_variable/TwTBXnLT">distribution</a> of <a href="/facts/Ideal_chain/JXqaeaPC">chain</a> lengths, so 
 
 
 
 x
 
 
 {\displaystyle x}
 
 would be an <a href="/facts/Average/etawWjw4">average</a>.</li>
<li>^ The <a href="/facts/Enthalpy/bBou8lCE">enthalpy</a> is the <a href="/facts/Internal_energy/IYdsoxhR">internal energy</a> corrected for any <a href="/facts/Pressure/eyc37GRG">pressure</a>-<a href="/facts/Volume/aeAXCBUd">volume</a> <a href="/facts/Mechanical_work/mbclFPBs">work</a> at constant (external) 
 
 
 
 P
 
 
 {\displaystyle P}
 
. We are not making any distinction here. This allows the approximation of <a href="/facts/Helmholtz_free_energy/2AK0HpPe">Helmholtz free energy</a>, which is the natural form of free energy from the Flory–Huggins lattice theory, to Gibbs free energy.</li>
<li>^ In fact, two of the sites adjacent to a polymer segment are occupied by other polymer segments since it is part of a <a href="/facts/Ideal_chain/JXqaeaPC">chain</a>; and one more, making three, for <a href="/facts/Branching_(polymer_chemistry)/KYF9YA1j">branching</a> sites, but only one for <a href="/facts/Polymer/eF32E50K">terminals</a>.</li></ol>

<h2 id="external-links">External links</h2>
<ul><li><a href="http://www.informit.com/content/images/chap3_0130181684/elementLinks/chap3_0130181684.pdf">"Conformations, Solutions and Molecular Weight" (book chapter)</a>, Chapter 3 of Book Title: Polymer Science and Technology; by Joel R. Fried; 2nd Edition, 2003</li></ul>

<h2 id="references">References</h2>

<ol>
<li id="fn:1">Burchard, W (1983). "Solution Thermodyanmics of Non-Ionic Water Soluble Polymers.". In Finch, C. (ed.). Chemistry and Technology of Water-Soluble Polymers. Springer. pp. 125–142. ISBN 978-1-4757-9661-2. <a href="978-1-4757-9661-2" target="_blank">978-1-4757-9661-2</a> <a href="#fnref:1" class="footnote-back-ref">↩</a></li>
<li id="fn:2">Franks, F (1983). "Water Solubility and Sensitivity-Hydration Effects.". In Finch, C. (ed.). Chemistry and Technology of Water-Soluble Polymers. Springer. pp. 157–178. ISBN 978-1-4757-9661-2. <a href="978-1-4757-9661-2" target="_blank">978-1-4757-9661-2</a> <a href="#fnref:2" class="footnote-back-ref">↩</a></li>
<li id="fn:3">Dijk, Menno A. van; Wakker, Andre (1998-01-14). Concepts in Polymer Thermodynamics. CRC Press. pp. 61–65. ISBN 978-1-56676-623-4. <a href="978-1-56676-623-4" target="_blank">978-1-56676-623-4</a> <a href="#fnref:3" class="footnote-back-ref">↩</a></li>
<li id="fn:4">Burchard, W (1983). "Solution Thermodyanmics of Non-Ionic Water Soluble Polymers.". In Finch, C. (ed.). Chemistry and Technology of Water-Soluble Polymers. Springer. pp. 125–142. ISBN 978-1-4757-9661-2. <a href="978-1-4757-9661-2" target="_blank">978-1-4757-9661-2</a> <a href="#fnref:4" class="footnote-back-ref">↩</a></li>
<li id="fn:5">Franks, F (1983). "Water Solubility and Sensitivity-Hydration Effects.". In Finch, C. (ed.). Chemistry and Technology of Water-Soluble Polymers. Springer. pp. 157–178. ISBN 978-1-4757-9661-2. <a href="978-1-4757-9661-2" target="_blank">978-1-4757-9661-2</a> <a href="#fnref:5" class="footnote-back-ref">↩</a></li>
<li id="fn:6">de Gennes, Pierre-Gilles (1979). Scaling concepts in polymer physics. Ithaca, N.Y.: Cornell University Press. ISBN 080141203X. OCLC 4494721. <a href="080141203X" target="_blank">080141203X</a> <a href="#fnref:6" class="footnote-back-ref">↩</a></li>
<li id="fn:7">Berry, J; et al. (2018). "Physical principles of intracellular organization via active and passive phase transitions". Reports on Progress in Physics. 81 (46601): 046601. Bibcode:2018RPPh...81d6601B. doi:10.1088/1361-6633/aaa61e. PMID 29313527. S2CID 4039711. <a href="/wiki/Bibcode_(identifier)" target="_blank">/wiki/Bibcode_(identifier)</a> <a href="#fnref:7" class="footnote-back-ref">↩</a></li>
<li id="fn:8">Doi, Masao (2013). Soft Matter Physics. Great Clarendon Street, Oxford, UK: Oxford University Press. ISBN 9780199652952. <a href="9780199652952" target="_blank">9780199652952</a> <a href="#fnref:8" class="footnote-back-ref">↩</a></li>
<li id="fn:9">Schmid, Friederike (2010). "Theory and Simulation of Multiphase Polymer Systems". arXiv:1001.1265 [cond-mat.soft]. <a href="/wiki/Friederike_Schmid" target="_blank">/wiki/Friederike_Schmid</a> <a href="#fnref:9" class="footnote-back-ref">↩</a></li>
<li id="fn:10">de Gennes, Pierre-Gilles (1979). Scaling concepts in polymer physics. Ithaca, N.Y.: Cornell University Press. ISBN 080141203X. OCLC 4494721. <a href="080141203X" target="_blank">080141203X</a> <a href="#fnref:10" class="footnote-back-ref">↩</a></li>
<li id="fn:11">Doi, Masao (2013). Soft Matter Physics. Great Clarendon Street, Oxford, UK: Oxford University Press. ISBN 9780199652952. <a href="9780199652952" target="_blank">9780199652952</a> <a href="#fnref:11" class="footnote-back-ref">↩</a></li>
</ol>

Flory–Huggins solution theory open-in-new

Flory–Huggins solution theory