Molecular dynamics simulation study of solvent effects on conformation and dynamics of polyethylene oxide and polypropylene oxide chains in water and in common organic solvents

New Insights on Solvent-Induced Changes in Refractivity and Specific Rotation of Poly(propylene oxide) Systems Extracted from Channeled Spectra

Investigation of chiroptical polymers in the solution phase is paramount for designing supramolecular architectures for photonic or biomedical devices. This work is devoted to the case study of poly(propylene oxide) (PPO) optical activity in several solvents: benzonitrile, carbon disulfide, chloroform, ethyl acetate, and p-dioxane. To attain information on the interactions in these systems, rheological testing was undertaken, showing distinct variations of the rheological parameters as a function of the solvent type. These aspects are also reflected in the refractive index dispersive behavior, from which linear and non-linear optical properties are extracted. To determine the circular birefringence and specific rotation of the PPO solutions, the alternative method of the channeled spectra was employed. The spectral data were correlated with the molecular modeling of the PPO structural unit in the selected solvents. Density functional theory (DFT) computational data indicated that the torsional potential energy-related to the O1-C2-C3-O4 dihedral angle from the polymer repeating unit-was hindered in solvation environments characterized by high polarity and the ability to interact via hydrogen bonding. This was in agreement with the optical characterization of the samples, which indicated a lower circular birefringence and specific rotation for the solutions of PPO in ethyl acetate and p-dioxane. Also, the shape of optical rotatory dispersion curves was slightly modified for PPO in these solvents compared with the other ones.

Unravelling the molecular interactions behind the formation of PEG/PPG aqueous two-phase systems

The understanding of molecular interactions that control phase separation in polymer/polymer aqueous two-phase systems (ATPS) has been a subject of debate up to this day. In light of this, we set out to investigate the molecular interactions occurring in ternary mixtures composed of polyethylene glycol (PEG600), polypropylene glycol (PPG400) and water. The ternary phase diagram was plotted at two temperatures (298 K and 323 K), revealing a transition from a type 0 to a type I diagram. Molecular dynamics (MD) simulations were performed to elucidate the polymer-polymer and polymer-water interactions occurring at different temperatures and water concentrations. COnductor-like Screening Model for Realistic Solvents (COSMO-RS) was used to assess the thermodynamic properties of the polymer-water binary mixtures and their correlation with ATPS formation. The MD simulations clearly demonstrate the effect of segregation/separation with increasing water content and temperature, highlighting a significant reduction in PPG-water interactions compared to PEG-water counterparts. Polymer-water interactions were identified as those controlling the phase separation mechanism, and the thermodynamic properties determined with COSMO-RS for the polymer-water binary systems further support this view.

Structure and Thermodynamics of Linear, Ring, and Catenane Polymers in Solutions and at Liquid-Liquid Interfaces

In recent decades, advances in the syntheses of mechanically interlocked macromolecules, such as catenanes, have led to much greater interest in the applications of these complexes, from molecular motors and actuators to nanoscale computational memory and nanoswitches. Much remains to be understood, however, regarding how catenated ring compounds behave as a result of the effects of different solvents as well as the effects of solvent/solvent interfaces. In this work, we have investigated, using molecular dynamics simulations, the effects of solvation of poly(ethylene oxide) chains of different topologies─linear, ring, and [2]catenane─in two solvents both considered favorable toward PEO (water, toluene) and at the water/toluene interface. Compared to ring and [2]catenane molecules, the linear PEO chain showed the largest increase in size at the water/toluene interface compared to bulk water or bulk toluene. Perhaps surprisingly, observations indicate that the tendency of all three topologies to extend at the water/toluene interface may have more to do with screening the interaction between the two solvents than with optimizing specific solvent-polymer contacts.

Methane Hydrate Behavior for Water–Oil Systems Containing CTAB and Synperonic PE/F127 Surfactants

Methane hydrates were studied in systems containing aqueous dissolved surfactants in oil emulsions with a volume ratio of 40/60. Two commercial surfactants, named synperonic PE/F127 and cetyltrimethylammonium bromide, were evaluated at 0, 350, 700 and 1500 ppm. Experiments were made by applying the cooling–heating path in an isochoric high-pressure cell at different initial pressures of 5.5, 8.0, 10.0 and 12.0 MPa. The obtained parameters were induction time, temperature onset, pressure drop, and dissociation conditions. The results revealed that the dissociation curve for methane in water-in-oil emulsions was not modified by the surfactants. The crystallization (onset) temperature was higher using synperonic PE/F127 in comparison with zero composition, while the opposite occurred with cetyltrimethylammonium bromide. Both surfactants induced a delaying effect on the induction time and a lesser pressure drop.

Hydrogen Bonding Sequence Directed Coil-Globule Transition in Water Soluble Thermoresponsive Polymers

The origin of the coil-globule transition for water-soluble thermoresponsive polymers frequently used in nanomaterials remains elusive. Using polypropylene oxide as an example we demonstrate by means of atomistic molecular dynamics simulations that temperature-induced increase in the sequence length of monomers that are not hydrogen bonded to water drives the coil-globule transition. Longer chains statistically exhibit longer sequences which serve as nucleation sites for hydrophobic cluster formation, facilitating chain collapse at lower temperature in agreement with experimental data.

Growth Mechanism of Polymer Membranes Obtained by H-Bonding Across Immiscible Liquid Interfaces

Complexation of polymers at liquid interfaces is an emerging technique to produce all-liquid printable and self-healing devices and membranes. It is crucial to control the assembly process, but the mechanisms at play remain unclear. Using two different reflectometric methods, we investigate the spontaneous growth of H-bonded PPO-PMAA (polypropylene oxide-polymetacrylic acid) membranes at a flat liquid-liquid interface. We find that the membrane thickness h grows with time t as h ∼ t^1/2, which is reminiscent of a diffusion-limited process. However, counterintuitively, we observe that this process is faster as the PPO molar mass increases. We are able to rationalize these results with a model which considers the diffusion of the PPO chains within the growing membrane. The architecture of the latter is described as a gel-like porous network, with a pore size much smaller than the radius of the diffusing PPO chains, thus inducing entropic barriers that hinder the diffusion process. From the comparison between the experimental data and the result of the model, we extract some key piece of information about the microscopic structure of the membrane. This study opens the route toward the rational design of self-assembled membranes and capsules with optimal properties.

The coarse-grained models of poly(ethylene oxide) and poly(propylene oxide) homopolymers and poloxamers in big multipole water (BMW) and MARTINI frameworks

Polyethylene oxide (PEO) and poly(propylene oxide) (PPO), especially their tri-block copolymers PEO-PPO-PEO (poloxamers), have a broad range of applications in biotechnology and medical science. Understanding their specific interactions with biomembranes is the key to unveil the unique features of poloxamers either as membrane-healing or membrane pore-forming agents. Based on the coarse-graining convention of the MARTINI force field and the big multipole water (BMW) model, which has a three charged site topology and can reproduce the correct dipole moment of four-water clusters, we generated coarse-grained (CG) models with analytical and numerical potentials for PEO and PPO homopolymers and poloxamers in dilute solution. The effective bonded interaction potentials between CG beads were determined from the probability distributions of bond lengths, angles and dihedrals that are determined from atomistic simulations. The nonbonded interaction parameters were fine-tuned to reproduce the conformational properties of atomistic PEO and PPO homopolymers and poloxamers via extensive CG simulations of PEO and PPO homopolymers and poloxamers in a BMW water environment. The reported CG models provide a promising framework for a comprehensive understanding of the microstructural, conformational, and dynamic properties of poloxamers and their delicate interactions with other species in an explicit water environment.

Spatial Distribution of PEO-PPO-PEO Block Copolymer and PEO Homopolymer in Lipid Bilayers

Maintaining the integrity of cell membranes is indispensable for cellular viability. Poloxamer 188 (P188), a poly(ethylene oxide)-b-poly(propylene oxide)-b-poly(ethylene oxide) (PEO-PPO-PEO) triblock copolymer with a number-average molecular weight of 8700 g/mol and containing 80% by mass PEO, protects cell membranes from various external injuries and has the potential to be used as a therapeutic agent in diverse applications. The membrane protection mechanism associated with P188 is intimately connected with how this block copolymer interacts with the lipid bilayer, the main component of a cell membrane. Here, we report the distribution of P188 in a model lipid bilayer comprising 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) using neutron reflectivity (NR) and atomic force microscopy (AFM). We also investigated the association of a PEO homopolymer (PEO8.4K; M_n = 8400 g/mol) that does not protect living cell membranes. These experiments were conducted following incubation of a 4.5 mmol/L polymer solution in a buffer that mimics physiological conditions with supported POPC bilayer membranes followed by washing with the aqueous medium. In contrast to previous reports, which dealt with P188 and PEO in salt-free solutions, both P188 and PEO8.4K penetrate into the inner portion of the lipid bilayer as revealed by NR, with approximately 30% by volume occupancy across the membrane without loss of bilayer structural integrity. These results indicate that PEO is the chemical moiety that principally drives P188 binding to bilayer membranes. No defects or phase-separated domains were observed in either P188- or PEO8.4K-incubated lipid bilayers when examined by AFM, indicating that polymer chains mingle homogeneously with lipid molecules in the bilayer. Remarkably, the breakthrough force required for penetration of the AFM tip through the bilayer membrane is unaffected by the presence of the large amount of P188 and PEO8.4K.

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.

What causes the anomalous aggregation in pluronic aqueous solutions?

Pluronic (PL) block copolymers have been widely used as delivery carriers, molecular templates for porous media, and process additives for affecting rheological behavior. Unlike most surfactant systems, where unimer transforms into micelle with increased surfactant concentration, anomalous large PL aggregates below the critical micelle concentration (CMC) were found throughout four types of PL (F108, F127, F88 and P84). We characterized their structures using dynamic light scattering and small-angle X-ray/neutron scattering. Molecular dynamics simulations suggest that the PPO segments, though weakly hydrophobic interaction (insufficient to form micelles), promote the formation of large aggregates. Addition of acid or base (e.g. citric acid, acetic acid, HCl and NaOH) in F108 solution significantly suppresses the aggregate formation for up to 20 days due to the repulsion force from the attached H3O+ molecules on the EO segment in both PEO and PL and the reduction of CMC through the salting out effect, respectively.

Molecular Dynamics Simulations of Strain-Induced Phase Transition of Poly(ethylene oxide) in Water

We study the dilute aqueous solutions of poly(ethylene oxide) (PEO) oligomers that are subject to an elongating force dipole acting on both chain ends using atomistic molecular dynamics. By increasing the force, liquid-liquid demixing can be observed at room temperature far below the lower critical solution temperature. For forces above 35 pN, fibrillar nanostructures are spontaneously formed related to a decrease in hydrogen bonding between PEO and water. Most notable is a rapid decrease in the bifurcated hydrogen bonds during stretching, which can also be observed for isolated single chains. The phase-segregated structures display signs of chain ordering, but a clear signature of the crystalline order is not obtained during the simulation time, indicating a liquid-liquid phase transition induced by chain stretching. Our results indicate that the solvent quality of the aqueous solution of PEO depends on the conformational state of the chains, which is most likely related to the specific hydrogen-bond-induced solvation of PEO in water. The strain-induced demixing of PEO opens the possibility to obtain polymer fibers with low energy costs because crystallization starts via the strain-induced demixing in the extended state only.

Charging and Release Mechanisms of Flexible Macromolecules in Droplets

We study systematically the charging and release mechanisms of a flexible macromolecule, modeled by poly(ethylene glycol) (PEG), in a droplet by using molecular dynamics simulations. We compare how PEG is solvated and charged by sodium Na⁺ ions in a droplet of water (H₂O), acetonitrile (MeCN), and their mixtures. Initially, we examine the location and the conformation of the macromolecule in a droplet bearing no net charge. It is revealed that the presence of charge carriers do not affect the location of PEG in aqueous and MeCN droplets compared with that in the neutral droplets, but the location of the macromolecule and the droplet size do affect the PEG conformation. PEG is charged on the surface of a sodiated aqueous droplet that is found close to the Rayleigh limit. Its charging is coupled to the extrusion mechanism, where PEG segments leave the droplet once they coordinate a Na⁺ ion or in a correlated motion with Na⁺ ions. In contrast, as PEG resides in the interior of a MeCN droplet, it is sodiated inside the droplet. The compact macro-ion transitions through partially unwound states to an extended conformation, a process occurring during the final stage of desolvation and in the presence of only a handful of MeCN molecules. For charged H₂O/MeCN droplets, the sodiation of PEG is determined by the H₂O component, reflecting its slower evaporation and preference over MeCN for solvating Na⁺ ions. We use the simulation data to construct an analytical model that suggests that the droplet surface electric field may play a role in the macro-ion-droplet interactions that lead to the extrusion of the macro-ion. This study provides the first evidence of the effect of the surface electric field by using atomistic simulations. Graphical Abstract ᅟ.

Entropy of aqueous surfaces. Application to polymeric Langmuir films

Measuring surface (excess) entropies provides a bounty of valuable structural information that is hard to obtain otherwise. In the paper these quantities are defined and procedures of measurements discussed. Mostly they involve measurements at different temperatures. A review is given for interfaces with aqueous solutions in the absence of polymers. This review illustrates how, sometimes unanticipated, pieces of information are obtained, for example with cloud seeding and a possible explanation of the Jones-Ray effect. As a novel extension the procedure is applied to deposited, or Langmuir, monolayers of poly(ethylene oxide)-poly(propylene oxide) block copolymers. It will be shown how the various phase transitions and associated configurations of these polymers can be recognized and monitored.

Effect of Chain Length on the Partition Properties of Poly(ethylene oxide): Comparison between MARTINI Coarse-Grained and Atomistic Models

The MARTINI coarse-grained beads are parameterized to match the partition coefficients of several organic molecules in different solvents. Here, we test the method when modeling the partitioning properties of poly(ethylene oxide) between solvents of different polarities. We show that, among the existing models, the latest model developed by Lee and co-workers [ Lee , H. ; Pastor , R. W. J. Phys. Chem. B 2011 , 115 , 7830 - 7837 ] is the one that most successfully reproduces the hydration free energy of short oligomers, although it predicts highly negative solvation free energies in octanol and hexane. We develop a new CG model matching the solvation free energy of the monomer in different solvents and propose a simple method to select the Lennard-Jones parameters that reproduce the desired partition coefficients. The model correctly reproduces water/hexane partition properties for oligomers up to 10 monomers but still suffers from a transferability problem for larger molecular weight.

The dynamics of solvation dictates the conformation of polyethylene oxide in aqueous, isobutyric acid and binary solutions

The dynamics of solvation dictates the conformation of polyethylene oxide in water and isobutyric acid causing a helix–coil transition in a mixed isobutyric acid/water solvent.

Native Silk Feedstock as a Model Biopolymer: A Rheological Perspective

Variability in silk's rheology is often regarded as an impediment to understanding or successfully copying the natural spinning process. We have previously reported such variability in unspun native silk extracted straight from the gland of the domesticated silkworm Bombyx mori and discounted classical explanations such as differences in molecular weight and concentration. We now report that variability in oscillatory measurements can be reduced onto a simple master-curve through normalizing with respect to the crossover. This remarkable result suggests that differences between silk feedstocks are rheologically simple and not as complex as originally thought. By comparison, solutions of poly(ethylene-oxide) and hydroxypropyl-methyl-cellulose showed similar normalization behavior; however, the resulting curves were broader than for silk, suggesting greater polydispersity in the (semi)synthetic materials. Thus, we conclude Nature may in fact produce polymer feedstocks that are more consistent than typical man-made counterparts as a model for future rheological investigations.

Unexpected Temperature Behavior of Polyethylene Glycol Spacers in Copolymer Dendrimers in Chloroform

We have studied copolymer dendrimer structure: carbosilane dendrimers with terminal phenylbenzoate mesogenic groups attached by poly(ethylene) glycol (PEG) spacers. In this system PEG spacers are additional tuning to usual copolymer structure: dendrimer with terminal mesogenic groups. The dendrimer macromolecules were investigated in a dilute chloroform solution by ¹H NMR methods (spectra and relaxations). It was found that the PEG layer in G = 5 generations dendrimer is “frozen” at high temperatures (above 260 K), but it unexpectedly becomes “unfrozen” at temperatures below 250 K (i.e., melting when cooling). The transition between these two states occurs within a small temperature range (~10 K). Such a behavior is not observed for smaller dendrimer generations (G = 1 and 3). This effect is likely related to the low critical solution temperature (LCST) of PEG and is caused by dendrimer conformations, in which the PEG group concentration in the layer increases with growing G. We suppose that the unusual behavior of PEG fragments in dendrimers will be interesting for practical applications such as nanocontainers or nanoreactors.

Tribological Behavior of Aqueous Copolymer Lubricant in Mixed Lubrication Regime

Although a number of experiments have been attempted to investigate the lubrication of aqueous copolymer lubricant, which is applied widely in metalworking operations, a comprehensive theoretical investigation at atomistic level is still lacking. This study addresses the influence of loading pressure and copolymer concentration on the structural properties and tribological performance of aqueous copolymer solution of poly(propylene oxide)-poly(ethylene oxide)-poly(propylene oxide) (PPO-PEO-PPO) at mixed lubrication using a molecular dynamic (MD) simulation. An effective interfacial potential, which has been derived from density functional theory (DFT) calculations, was employed for the interactions between the fluid's molecules and iron surface. The simulation results have indicated that the triblock copolymer is physisorption on iron surface. Under confinement by iron surfaces, the copolymer molecules form lamellar structure in aqueous solution and behave differently from its bulk state. The lubrication performance of aqueous copolymer lubricant increases with concentration, but the friction reduction is insignificant at high loading pressure. Additionally, the plastic deformation of asperity is dependent on both copolymer concentration and loading pressure, and the wear behavior shows a linear dependence of friction force on the number of transferred atoms between contacting asperities.

Multiscale Molecular Dynamics Simulations of Model Hydrophobically Modified Ethylene Oxide Urethane Micelles

The flower-like micelles of various aggregation numbers of a model hydrophobically modified ethylene oxide urethane (HEUR) molecule, C16E45C16, and their corresponding starlike micelles, containing the surfactants C16E22 and C16E23, were studied by atomistic and coarse-grained molecular dynamic (MD) simulations. We used free energies from umbrella sampling to calculate the size distribution of micelle sizes and the average time for escape of a hydrophobic group from the micelle. Using the coarse-grained MARTINI force field, the most probable size of the model HEUR molecule was thereby determined to be about 80 hydrophobes per micelle and the average hydrophobe escape time to be about 0.1 s, both of which are consistent with previous experimental studies. Atomistic simulations reveal that hydrogen bond formation and the mean lifetime of hydration waters of the poly(ethylene oxide) (or PEO) groups are location-dependent in the HEUR micelle, with PEO groups immediately adjacent to the C16 groups forming the fewest hydrogen bonds with water and having hydration waters with longer lifetimes than those of the PEO groups located further away from the C16 groups.

Surfactant or block copolymer micelles? Structural properties of a series of well-defined n-alkyl-PEO micelles in water studied by SANS

Here we present an extensive small-angle neutron scattering (SANS) structural characterization of micelles formed by poly(ethylene oxide)-mono-n-alkyl ethers (Cn-PEOx) in dilute aqueous solution. Chemically, Cn-PEOx can be considered as a hybrid between a low-molecular weight surfactant and an amphiphilic block copolymer. The present system, prepared through anionic polymerization techniques, is better defined than other commercially available polymers and allows a very precise and systematic testing of the theoretical predictions from thermodynamical models. The equilibrium micellar properties were elaborated by systematically varying the n-alkyl chain length (n) at constant PEO molecular weight or increasing the soluble block size (x), respectively. The structure was reminiscent of typical spherical star-like micelles i.e. a constant core density profile, ∼r(0), and a diffuse corona density profile, ∼r(-4/3). Through a careful quantitative analysis of the scattering data, it is found that the aggregation number, Nagg initially rapidly decreases with increasing PEO length until it becomes independent at higher PEO molecular weight as expected for star-like micelles. On the other hand, the dependency on the n-alkyl length is significantly stronger than that expected from the theories for star-like block copolymer micelles, Nagg ∼ n(2) similar to what is expected for surfactant micelles. Hence the observed aggregation behavior suggests that the Cn-PEOx micelles exhibit a behavior that can be considered as a hybrid between low-molecular weight surfactant micelles and diblock copolymer micelles.

Coarse-graining poly(ethylene oxide)-poly(propylene oxide)-poly(ethylene oxide) (PEO-PPO-PEO) block copolymers using the MARTINI force field

The MARTINI coarse-grain (CG) force field is extended for a class of triblock block copolymers known as Pluronics. Existing MARTINI bead types are used to model the non-bonded part of the potential while single chain properties for both homopolymers, poly(ethylene oxide) (PEO) and poly(propylene oxide) (PPO), are used to develop the bonded interactions. The new set of force field parameters reproduces structural and dynamical properties of high molecular weight homo- and copolymers. The CG model is moderately transferable in solvents of different polarity and concentration; however, the PEO homopolymer model presents a reduced thermodynamic transferability especially in water probably due to the lack of hydrogen bonds with the solvent. Our simulations of a monolayer of Pluronic L44 show polymer-brush-like characteristics for the PEO segments which protrude into the aqueous phase. Other membrane properties not easily accessible using experimental techniques such as its membrane thickness are also calculated.