Modeling of Aqueous Poly(oxyethylene) Solutions: 1. Atomistic Simulations

Citations Over TimeTop 10% of 2008 papers

Abstract

The performance of different recently proposed force fields in combination with TIP4P-Ewald (TIP4P-Ew) water in reproducing experimental data of liquid 1,2-dimethoxyethane (DME) and its aqueous solutions for conformer populations, densities of solutions, and self-diffusion coefficients was explored. A modified version of the OPLS force field ("engineered") showed best performance in describing the conformer equilibria, but extremely high interconformational barriers reduce its applicability in dynamical simulations. The TraPPE-united atom force field (TraPPE-UA) by Siepmann et al. (J. Phys. Chem. B 2004, 108, 17596) was found to perform best in reproducing thermodynamic properties, but it showed some deficiency in describing the conformer equilibria. We reparameterized the dihedral potentials to match recent ab initio data by Anderson and Wilson (Mol. Phys. 2005, 103, 89) and could improve significantly the performance of description of conformer populations of DME in water. Subsequently, this modified TraPPE-UA was used in extensive simulations of poly(oxyethylene) oligomers H(CH2OCH2)nH (POEn) with n = 3, 5, 10, 12, 20, 30 repeat units at mass fractions between 3% and 80% at 298 K. Density, radii of gyration, and diffusion coefficients are in very good agreement with available experimental data. We conclude that this force field in combination with the TIP4P-Ew water model is very suitable for simulations of poly(oxyethylene) oligomers in aqueous solution. The application to real polymeric systems on the atomistic level is however hindered by very slow decorrelation of large-scale features and by slow diffusion.

Related Papers

• → Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties in SPC/E, TIP3P-fb, TIP4P/2005, TIP4P-Ew, and TIP4P-D(2021)44 cited
• → Molecular Dynamics Simulations of Proteins: Can the Explicit Water Model Be Varied?(2007)71 cited
• → Rotamer libraries and probabilities of transition between rotamers for the side chains in protein–protein binding(2012)26 cited
• → Solid and liquid phase Raman and infrared spectra of diallyl sulfide and diallyl disulfide(1989)16 cited
• → Explicit and GB/SA solvents: Each with two different force fields in multicanonical conformational sampling of a 25-residue polypeptide(2006)5 cited