Improved Temperature Behavior of PNIPAM in Water with a Modified OPLS Model

We test the OPLS/AA force field for a single PNIPAM 40-mer in aqueous solution using replica exchange molecular dynamics simulations and find that the force field fails to reproduce the experimental temperature behavior. To resolve this issue, we apply a modification on the partial charges previously suggested to reproduce the liquid-liquid phase separation of NIPAM aqueous solutions. The modified force field features stronger amide-water electrostatic interactions than the original OPLS model, predicts a weaker water-mediated monomer-monomer attraction, and reproduces the experimental coil-globule collapse enthalpy of PNIPAM in water. We revisit the cononsolvency problem of PNIPAM in methanol/water mixtures with the modified model and show that the dependence of the coil-globule collapse enthalpy on methanol concentration follows the experimental trend of the lower critical solution temperature. The calculations with the modified force field confirm that polymer dehydration is the determining factor for chain collapse in the cononsolvency regime.