Polymer Viscosity

Quick Start

We would love to know how the viscosity of our polymer solution depends on concentration and on MWt. Unfortunately it's complicated, as can be visualised in the Polymers in Solvents app. Here we do our best. It's probably best to play with any parameter you recognise (if you've tried Polymers in Solvents then you already know a lot) then you will have to read things in details. It's complicated and there seem to be no short cuts.

Polymer Viscosity

Rg @0.5+40K nm
Mmono
Polymer MWt
Polymer Mc
Flory-Huggin χ
Cmax g/l
Rg @χ+MWt nm
Ne
C* g/dl
Ce g/dl
|η| dl/g
N.bonds
Mbond
Bond length nm
MWt
Rg nm

Here's a simple problem. You have a polymer and a solvent and you want some idea of how the viscosity of that solution depends on the polymer's MWt and how "happy" it is in the solvent.

It turns out that the answer is far from simple and I am not aware of any really good ways to get an answer. So I have created this plausible app based on a formula1 that makes plenty of sense but which does not seem widly known or used.

The viscosity we want is the Relative Viscosity, ηr which is that of the polymer solution ηp divided by that of the pure solvent ηs. We want to know how it depends on concentration, C, (in units of g/dl, i.e. g/100ml which is ~ weight %) up to a maximum value, Cmax. The first of many equations tells us that:

ηrps

First start with the default inputs and vary MWt and then the χ parameter (0=perfect solvent up to 0.5=theta solvent) get a feel for what happens, then read more details below.

Now we can take things step by step. The key facts of polymer life is that at a concentration C* the polymer chains start touching each other, significantly increasing viscosity, transitioning from the "dilute" to "semi-dilute" regime, then at Ce they enter the entangled domain and viscosity really starts to increase. It is this concentration zone where most formulators find themselves. Incidentally, most polymer scientists are much happier discussing things in the entirely uninteresting dilute regime below C*.

We can calculate C* by knowing the MWt, Avogadro's number Av and the Radius of Gyration, Rg.

C* = MWt/(Av.Rg³)

So now we need to know Rg. There is a well-known formula based on the number N of pseudo-monomers of length b, and in the bottom row you can try out this formula. What is b? Typically it's a few real monomer units, but for your specific polymer, who knows. Use your instincts. If you set things so that the MWt of this pseudo-polymer is 40K then you can use the estimated Rg in the real calculation. Alternatively Rg can be measured via scattering experiments for those who are keen.

But Rg depends on MWt and on how "happy" the polymer is in the solvent. To simplify the app as much as possible, you input Rg @0.5+40K which means the value in a theta solvent (χ=0.5) for a 40K MWt polymer, a value I chose for no special reason. As you change the MWt and χ, a new Rg@χ+MWt is calculated and shown in the row below.

To calculate the viscosity in the domain above C* we use the Rouseian formula:

ηr = (|η|C)1.3

Which means that we need to calculate the Intrinsic Viscosity |η| which is a limiting value of (ηps)/(C.ηs) in units of 1/concentration. We get it from:

|η| = 5AvRg³/MWt

Finally, we can calculate our relative viscosity by multiplying the Rouseian viscosity with a term dependent on C/Ce:

ηr = (|η|C)1.345(C/Ce)³

Ah, but for this we need Ce which comes from:

Ce = Ne3v-1C*/6.6

And for this we need v (which varies from 0.5 at theta to 0.585 for χ=0) and Ne, which is the number of monomer units between entanglements which depends on the two inputs we've not mentioned so far, Mc, the critical entanglement MWt above which the polymer tangles, and Mmono which is the monomer MWt.

Ne = Mc/3/Mmono

Should we bother?

Thank goodness for apps (and for some VB code I wrote 10 years ago when I first tried to understand all this). At least the tedious calculations are taken care of in a plausible manner. But how does it apply to your polymer in your solvent system? If you don't know Rg and Mc (and most of us don't) then it might all seem a waste of time. But hold on, we're not completely lost. If you happen to have just a single ηr measurement and you know your MWt and can guess whether the solvent is really good (χ=0), neutral (χ=0.5) or in between, then you can make some plausible guesses for Mc and Rg, so you can at least get a feel of what might happen if you changed concentration or MWt. With two datapoints then you will have a quite reasonable idea of the whole system. And, finally:

If you happen to be a polymer physicist with some much better ideas of how to make a usable app for the formulation community, and/or if you see issues with the current version, please let me know.

1The formula, and the calculations of the intermediate values come from: Youngsuk Heo and Ronald G. Larson, The scaling of zero-shear viscosities of semidilute polymer solutions with concentration, J. Rheol. 49, 1117-1128, 2005. Heo's excellent PhD thesis with more details can readily be found online.