The structure of poly(ethylene oxide) liquids: comparison of integral equation theory with molecular dynamics simulations and neutron scattering

Effects of intramolecular chain conformation on the hydration and miscibility of polyethylene glycol in water studied by means of polymer reference interaction site model theory

To examine the conventional idea that the gauche conformation of the OCCO dihedral angle promotes the dissolution of polyethylene glycol (PEG) in water through strong hydration, the thermodynamic properties of liquid mixtures of PEG and water were studied by means of polymer reference interaction site model (PRISM) theory. The intramolecular correlation functions required as input for PRISM theory were calculated by the generator matrix method, accompanied by changes in the distribution of dihedral angles. In the infinite dilution limit, the increased probability of gauche conformation of the OCCO dihedral angles stabilizes the hydration of PEG through enhanced hydrogen bonding between the ether oxygen of PEG and water. The mixing Gibbs energies of the liquid mixtures were also calculated in the whole concentration range based on the Gibbs-Duhem equation, as per our recent proposal. A liquid-liquid phase separation was observed when all the dihedral angles of PEG were in the trans conformation; for the liquid mixture to be miscible in the whole concentration range, the introduction of the OCCO gauche conformation was found to be indispensable. The above theoretical results support the conventional idea that the OCCO gauche conformation is important for the high miscibility of PEG and water.

Hydrogen Bond Network Disruption by Hydration Layers in Water Solutions with Salt and Hydrogen-Bonding Polymers (PEO)

A mean field theory model describing the interaction of ion hydration layers with the network of hydrogen bonds of both water and the nonionic polymer poly(ethylene oxide) (PEO) is presented. The predictions of the model for types and statistics of hydrogen bonds, the number of water molecules bound to PEO, or their dependence on temperature are successfully verified from all-atom simulations at different NaCl and PEO concentrations. Furthermore, our simulations show that the binding of cations to PEO increases monotonically with salt concentration, in agreement with recent experimental results, through a mechanism in which the sum of the number of bound water and cations is independent of salt concentration. The model introduced is general and can describe any salt or hydrogen-bond-forming polymer.

On the origin of the extremely different solubilities of polyethers in water

The solubilities of polyethers are surprisingly counter-intuitive. The best-known example is the difference between polyethylene glycol ([-CH₂-CH₂-O-]_n) which is infinitely soluble, and polyoxymethylene ([-CH₂-O-]_n) which is completely insoluble in water, exactly the opposite of what one expects from the C/O ratios of these molecules. Similar anomalies exist for oligomeric and cyclic polyethers. To solve this apparent mystery, we use femtosecond vibrational and GHz dielectric spectroscopy with complementary ab initio calculations and molecular dynamics simulations. We find that the dynamics of water molecules solvating polyethers is fundamentally different depending on their C/O composition. The ab initio calculations and simulations show that this is not because of steric effects (as is commonly believed), but because the partial charge on the O atoms depends on the number of C atoms by which they are separated. Our results thus show that inductive effects can have a major impact on aqueous solubilities.

Integral equation theory for atactic polystyrene nanocomposite melts with a multi-site model

In this work, a multi-site chain model was incorporated into the polymer reference interaction site model to investigate the structure and properties of atactic polystyrene (aPS) melt and the structural correlations of dilute spherical nanoparticles dissolved in aPS melt. The theoretically calculated X-ray scattering intensities, solubility parameters and intermolecular correlation functions of aPS and its nanocomposites are found to be in agreement with the corresponding molecular simulation and experimental data. The theory was further employed to investigate the distribution functions of different size effects of aPS-nanoparticle system with consideration of the potential of mean force and depletion force. The aggregation of large nanoparticles increases with the increase of the nanoparticle-site size ratio in the infinitely dilute limit. The results show that the present theory can be used to investigate the structure of aPS melt and its nanocomposite, and give a further understanding of the filler dispersion and aggregation. All the observations indicate molecular-level details of the underlying mechanisms, providing useful information for the future design control of new aPS-nanocomposite materials with tailored properties.

Retention in gas–liquid chromatography with a polyethylene oxide stationary phase: Molecular simulation and experiment

Configurational-bias Monte Carlo simulations in the isobaric-isothermal Gibbs ensemble were carried out to investigate the partitioning of normal alkanes, primary and secondary alcohols, symmetric alkyl ethers and arenes between a helium vapor phase and a polyethylene oxide stationary phase (M(W)=382 g mol(-1)). The united-atom version of the transferable potentials for phase equilibria force field was used to model all solutes, polyethylene oxide and helium. The Gibbs free energies of transfer and Kovats retention indices of the solutes were calculated directly from the partition constants at two different temperatures, 353 and 393 K. Chromatographic experiments on a Carbowax 20M retentive phase were performed for the same set of solutes and temperatures ranging from 333 to 413 K. The predicted retention indices for alcohols, ethers and arenes are overestimated by about 120, 70 and 20 retention index units, respectively, pointing to an overestimation of the first-order electrostatic interactions in the model system. Molecular-level analysis shows that hydrogen-bonding and dipole-dipole interactions lead to orientational ordering for the alcohol and ether analytes, whereas the weaker dipole-quadrupole interactions for the arene solutes are not sufficient to induce orientational ordering. The retention indices of alcohols and ethers decrease with increasing temperature because of the large entropic cost of hydrogen-bonding and orientational ordering. In contrast, the retention indices for arenes increase with increasing temperature because the entropic cost of cavity formation is smaller for arenes than for comparable alkanes.