A coarse-grained model for polyethylene oxide and polyethylene glycol: conformation and hydrodynamics

A coarse-grained (CG) model for polyethylene oxide (PEO) and polyethylene glycol (PEG) developed within the framework of the MARTINI CG force field (FF) using the distributions of bonds, angles, and dihedrals from the CHARMM all-atom FF is presented. Densities of neat low molecular weight PEO agree with experiment, and the radius of gyration R(g) = 19.1 A +/- 0.7 for 76-mers of PEO (M(w) approximately 3400), in excellent agreement with neutron scattering results for an equal sized PEG. Simulations of 9, 18, 27, 36, 44, 67, 76, 90, 112, 135, and 158-mers of the CG PEO (442 < M(w) < 6998) at low concentration in water show the experimentally observed transition from ideal chain to real chain behavior at 1600 < M(w) < 2000, in excellent agreement with the dependence of experimentally observed hydrodynamic radii of PEG. Hydrodynamic radii of PEO calculated from diffusion coefficients of the higher M(w) PEO also agree well with experiment. R(g) calculated from both all-atom and CG simulations of PEO76 at 21 and 148 mg/cm(3) are found to be nearly equal. This lack of concentration dependence implies that apparent R(g) from scattering experiments at high concentration should not be taken to be the chain dimension. Simulations of PEO grafted to a nonadsorbing surface yield a mushroom to brush transition that is well described by the Alexander-de Gennes formalism.