Modeling of aqueous poly(oxyethylene) solutions: 1. Atomistic simulations

Computing Accurate True Self-Diffusion Coefficients and Shear Viscosities Using the OrthoBoXY Approach

In a recent paper [Busch, J.; Paschek, D. J. Phys. Chem. B 2023, 127, 7983-7987], we have shown that for molecular dynamics (MD) simulations using orthorhombic periodic boundary conditions with "magic" box length ratios of Lz/Lx = Lz/Ly = 2.7933596497, the self-diffusion coefficients Dx and Dy in x- and y-directions are independent of the system size. They both represent the true self-diffusion coefficient D0 = (Dx + Dy)/2, while the shear viscosity can be calculated from diffusion coefficients in x-, y-, and z-directions, using η = kBT·8.1711245653/[3πLz(Dx + Dy - 2Dz)]. In this contribution, we test this "OrthoBoXY" approach by its application to a variety of different systems: liquid water, dimethyl ether, methanol, triglyme, water/methanol mixtures, water/triglyme mixtures, and imidazolium-based ionic liquids. The chosen systems range from small-sized molecular liquids to complex mixtures and ionic liquids, while spanning a viscosity range of almost 3 orders of magnitude. We assess the efficiency of the method for computing true self-diffusion and viscosity data and provide simple formulas for estimating the required MD simulation lengths and sizes for delivering reliable data with targeted uncertainty levels. Our analysis of the system size dependence of statistical uncertainties for both the viscosity and the self-diffusion coefficient leads us to the conclusion that it is preferable to extend the simulation length instead of increasing the system size. MD simulations consisting of 768 molecules or ion pairs seem to be perfectly adequate.

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.

Simulational Tests of the Rouse Model

An extensive review of literature simulations of quiescent polymer melts is given, considering results that test aspects of the Rouse model in the melt. We focus on Rouse model predictions for the mean-square amplitudes ⟨(Xp(0))2⟩ and time correlation functions ⟨Xp(0)Xp(t)⟩ of the Rouse mode Xp(t). The simulations conclusively demonstrate that the Rouse model is invalid in polymer melts. In particular, and contrary to the Rouse model, (i) mean-square Rouse mode amplitudes ⟨(Xp(0))2⟩ do not scale as sin-2(pπ/2N), N being the number of beads in the polymer. For small p (say, p≤3) ⟨(Xp(0))2⟩ scales with p as p-2; for larger p, it scales as p-3. (ii) Rouse mode time correlation functions ⟨Xp(t)Xp(0)⟩ do not decay with time as exponentials; they instead decay as stretched exponentials exp(-αtβ). β depends on p, typically with a minimum near N/2 or N/4. (iii) Polymer bead displacements are not described by independent Gaussian random processes. (iv) For p≠q, ⟨Xp(t)Xq(0)⟩ is sometimes non-zero. (v) The response of a polymer coil to a shear flow is a rotation, not the affine deformation predicted by Rouse. We also briefly consider the Kirkwood-Riseman polymer model.

Extension of the TraPPE Force Field for Battery Electrolyte Solvents

Optimizing electrolyte formulations is key to improving performance of Li-/Na-ion batteries, where transport properties (diffusion coefficient, viscosity) and permittivity need to be predicted as functions of temperature, salt concentration and solvent composition. More efficient and reliable simulation models are urgently needed, owing to the high cost of experimental methods and the lack of united-atom molecular dynamics force fields validated for electrolyte solvents. Here the computationally efficient TraPPE united-atom force field is extended to be compatible with carbonate solvents, optimizing the charges and dihedral potential. Computing the properties of electrolyte solvents, ethylene carbonate (EC), propylene carbonate (PC), dimethyl carbonate (DMC), diethyl carbonate (DEC), and dimethoxyethane (DME), we observe that the average absolute errors in the density, self-diffusion coefficient, permittivity, viscosity, and surface tension are approximately 15% of the corresponding experimental values. Results compare favorably to all-atom CHARMM and OPLS-AA force fields, offering computational performance improvement of at least 80%. We further use TraPPE to predict the structure and properties of LiPF₆ salt in these solvents and their mixtures. EC and PC form complete solvation shells around Li⁺ ions, while the salt in DMC forms chain-like structures. In the poorest solvent, DME, LiPF₆ forms globular clusters despite DME's higher permittivity than DMC.

Strongly Solvating Ether Electrolytes for High-Voltage Lithium Metal Batteries

Ethers are promising electrolytes for lithium (Li) metal batteries (LMBs) because of their unique stability with Li metal. Although intensive research on designing anion-enriched electrolyte solvation structures has greatly improved their electrochemical stabilities, ether electrolytes are approaching an anodic bottleneck. Herein, we reveal the strong correlation between electrolyte solvation structure and oxidation stability. In contrast to previous designs of weakly solvating solvents for enhanced anion reactivities, the triglyme (G3)-based electrolyte with the largest Li⁺ solvation energy among different linear ethers demonstrates greatly improved stability on Ni-rich cathodes under an ultrahigh voltage of 4.7 V (93% capacity retention after 100 cycles). Ether electrolytes with a stronger Li⁺ solvating ability could greatly suppress deleterious oxidation side reactions by decreasing the lifetime of free labile ether molecules. This study provides critical insights into the dynamics of the solvation structure and its significant influence on the interfacial stability for future development of high-efficiency electrolytes for high-energy-density LMBs.

Phase Behavior of Alkyl Ethoxylate Surfactants in a Dissipative Particle Dynamics Model

We present a dissipative particle dynamics (DPD) model capable of capturing the liquid state phase behavior of nonionic surfactants from the alkyl ethoxylate (C_nE_m) family. The model is based upon our recent work [Anderson et al. J. Chem. Phys. 2017, 147, 094503] but adopts tighter control of the molecular structure by setting the bond angles with guidance from molecular dynamics simulations. Changes to the geometry of the surfactants were shown to have little effect on the predicted micelle properties of sampled surfactants, or the water-octanol partition coefficients of small molecules, when compared to the original work. With these modifications the model is capable of reproducing the binary water-surfactant phase behavior of nine surfactants (C₈E₄, C₈E₅, C₈E₆, C₁₀E₄, C₁₀E₆, C₁₀E₈, C₁₂E₆, C₁₂E₈, and C₁₂E₁₂) with a good degree of accuracy.

Dynamics of Polymer Chains in Poly(ethylene oxide)/Silica Nanocomposites via a Combined Computational and Experimental Approach

The dynamics of polymer chains in poly(ethylene oxide)/silica (PEO/SiO₂) nanoparticle nanohybrids have been investigated via a combined computational and experimental approach involving atomistic molecular dynamics simulations and dielectric relaxation spectroscopy (DRS) measurements. The complementarity of the approaches allows us to study systems with different polymer molecular weights, nanoparticle radii, and compositions across a broad range of temperatures. We study the effects of spatial confinement, which is induced by the nanoparticles, and chain adsorption on the polymer's structure and dynamics. The investigation of the static properties of the nanocomposites via detailed atomistic simulations revealed a heterogeneous polymer density layer at the vicinity of the PEO/SiO₂ interface that exhibited an intense maximum close to the inorganic surface, whereas the bulk density was reached for distances ∼1-1.2 nm away from the nanoparticle. For small volume fractions of nanoparticles, the polymer dynamics, probed by the atomistic simulations of low-molecular-weight chains at high temperatures, are consistent with the presence of a thin adsorbed layer that exhibits slow dynamics, with the dynamics far away from the nanoparticle being similar to those in the bulk. However, for high volume fractions of nanoparticles (strong confinement), the dynamics of all polymer chains were predicted slower than that in the bulk. On the other hand, similar dynamics were found experimentally for both the local β-process and the segmental dynamics for high-molecular-weight systems measured at temperatures below the melting temperature of the polymer, which were probed by DRS. These differences can be attributed to various parameters, including systems of different molecular weights and nanoparticle states of dispersion, the different temperature range studied by the different methods, the potential presence of a reduced-mobility PEO/SiO₂ interfacial layer that does not contribute to the dielectric spectrum, and the presence of amorphous-crystalline interfaces in the experimental samples that may lead to a different dynamical behaviors of the PEO chains.

Solutions and Condensed Phases of PEG2000 from All-Atom Molecular Dynamics

Extensive all-atom molecular dynamics studies of polyethylene glycol (PEG₂₀₀₀) when solvated and in the polymer bulk condensed phases were performed across a wide temperature range. We proposed two modified all-atom force field and observed the fate of the PEG₂₀₀₀ macromolecule when solvated in water, water with 4% ethanol, and ethyl acetate. In aqueous solutions, the macromolecule collapsed into a prolate spheroidal ball-like structure while adopting a rather elongated coiled structure in ethyl acetate. Inspection of the polymer-condensed phases across the 150-340 K temperature range enabled the atomistic view of the solid glass below the glass transition temperature of 230 K < T_g < 250 K and the rubber behavior above T_g. Predicted properties include the enthalpy, density, and cohesive energy temperature behavior, the specific heat, thermal expansivity, thermal compressibility, bulk modulus, and Hildebrand solubility parameter both below and above T_g. Within the polymer matrix, the PEG₂₀₀₀ macromolecules were entangled displaying a wide distribution of sizes that persisted when transitioning from the glass to the rubbery phases. Calculated properties agree very well with experiments when available or stand as crucial predictions while awaiting experimental measurement. Understanding the thermodynamics and structure of this useful polymer enables the efficient prediction of its behavior when building novel composite materials for nanomedicine and nanotherapeutics.

Analysis of folded structure and folding thermodynamics in heterogeneous-backbone proteomimetics

Recent years have seen a growing number of examples of designed oligomeric molecules with artificial backbone connectivity that are capable of adopting complex folded tertiary structures analogous to those seen in natural proteins. A range of experimental techniques from structural biology and biophysics have been brought to bear in the study of these proteomimetic agents. Here, we discuss some considerations encountered in the characterization of high-resolution folded structure as well as folding thermodynamics of protein-like artificial backbones. We provide an overview of the use of X-ray crystallography and NMR spectroscopy in such systems and review example applications of these methods in the primary literature. Further, we provide detailed protocols for two experiments that have proved useful in our prior and ongoing efforts to compare folding thermodynamics between natural protein domains and heterogeneous-backbone counterparts.

Mikto-Arm Stars as Soft-Patchy Particles: From Building Blocks to Mesoscopic Structures

We present an atomistic molecular dynamics study of self-assembled mikto-arm stars, which resemble patchy-like particles. By increasing the number of stars in the system, we propose a systematic way of examining the mutual orientation of these fully penetrable patchy-like objects. The individual stars maintain their patchy-like morphology when creating a mesoscopic (macromolecular) self-assembled object of more than three stars. The self-assembly of mikto-arm stars does not lead to a deformation of the stars, and their shape remains spherical. We identified characteristic sub-units in the self-assembled structure, differing by the mutual orientation of the nearest neighbor stars. The current work aims to elucidate the possible arrangements of the realistic, fully penetrable patchy particles in polymer matrix and to serve as a model system for further studies of nanostructured materials or all-polymer nanocomposites using the mikto-arm stars as building blocks.

Insights into the Thermal Response of a Poly(ethylene oxide)-poly(propylene oxide)-poly(ethylene oxide) Triblock Polymer in Water

A thermal responsive block copolymer made up of ethylene oxide (EO) and propylene oxide (PO) blocks was simulated with optimized atomistic potentials and enhanced sampling methods over a range of temperatures. The results for the L42 pluronic polymer (EO)₄(PO)₂₂(EO)₄, which is known to undergo a transition in this temperature range, and the similarly sized (EO)₃₀ polymer, which does not, are compared. The thermal responsive L42 polymers in a dilute solution tend to aggregate, and this tendency gets stronger as temperature increases. The poly(ethylene oxide) polymer shows no such tendency. The aggregation is stabilized by the hydrophobic contact of the propylene oxide methyl groups, which outweighs a small loss in hydrogen bonds between the ether oxygens and water.

Effectiveness of coarse graining degree and speedup on the dynamic properties of homopolymer

Evaluation of effective coarse graining (CG) degree and reasonable speedup relative to all-atomistic (AA) model was conducted to provide a basis for building appropriate larger-scale model. The reproducibility of atomistic conformation and temperature transferability both act as the analysis criteria to resolve the maximum acceptable CG degree. Taking short- and long time spans into account simultaneously in the estimation of computational speedup, a dynamic scaling factor is accessible in fitting mean squared displacement ratio of CG to AA as an exponential function. Computing loss in parallel running is an indispensable component in acceleration, which was also added in the evaluation. Subsequently, a quantified prediction of CG speedup arises as a multiplication of dynamic scaling factor, computing loss, time step, and the square of reduction in the number of degrees of freedom. Polyethylene oxide was adopted as a reference system to execute the direct Boltzmann inversion and iterative Boltzmann inversion. Bonded and non-bonded potentials were calculated in CG models with 1~4 monomers per bead. The effective CG degree was determined as two at the most with a speedup of four orders magnitude over AA in this study. Determination of effectiveness CG degree and the corresponding speedup prediction provide available tools in larger spatiotemporal-scale calculations.

How Does the Number of Arms Affect the Properties of Mikto-Arm Stars in a Selective Oligomeric Matrix? Insights from Atomistic Simulations

We present a simulation study of amphiphilic mikto-arm star copolymers in a selective polymer host. By means of atomistic molecular dynamics simulations, we examine the structural and dynamical properties of mikto-arm stars with varying number, n, of poly(ethylene oxide) (PEO) and polystyrene (PS) arms, (PEO) n (PS) n in a 33% wt blend with an oligomeric PEO host (o-PEO). As the number of arms increases, the stars resemble more spherical particles with less separated PEO and PS intramolecular domains. As a result of their internal morphology and associated geometrical constraints, the mikto-arm stars self-assemble either into cylindrical-like objects or a percolated network with increasing n, within the o-PEO matrix. The segmental dynamics is mostly governed by the star architecture and the heterogeneous local environment, formed by the intra- and intermolecular nanosegregation. We discuss the role of each factor and compare the results with previously published studies on mikto-arm stars.

Analysis of Atomistic Potentials for Poly(ethylene glycol) Ethers

Two different potentials, the modified TraPPE-UA model of Fischer, J. [ J. Phys. Chem. B 2008, 112, 2388-2398] and the modified general AMBER force-field (GAFF) model of Barbosa, N. S. V. [ J. Mol. Model. 2017, 23, 194]are tested for a variety of temperature-dependent properties for neat and aqueous solutions for small poly(oxyethylene) (PEO) oligomers and larger polymers. A set of charges for PEO of arbitrary size is adapted for the modified GAFF model. Both models accurately reproduce experimental properties, but we find that the modified TraPPE-UA model is more accurate.

Water dynamics affects thermal transport at the surface of hydrophobic and hydrophilic irradiated nanoparticles

Plasmonic nanoparticles, such as Au nanoparticles (NPs) coated with bio-compatible ligands, are largely studied and tested in nanomedicine for photothermal therapies. Nevertheless, no clear physical interpretation is currently available to explain thermal transport at the nanoparticle surface, where a solid-liquid (core-ligand) interface is coupled to a liquid-liquid (ligand-solvent) interface. This lack of understanding makes it difficult to control the temperature increase imposed by the irradiated NPs to the surrounding biological environment, and it has so far hindered the rational design of the NP surface chemistry. Here, atomistic molecular dynamics simulations are used to show that thermal transport at the nanoparticle surface depends dramatically on solvent diffusivity at the ligand-solvent interface. Furthermore, using physical indicators of water confinement around hydrophobic and hydrophilic ligands, a predictive model is developed to allow the engineering of NP coatings with the desired thermal conductivities at the nanoscale.

Distance Distribution between Two Iodine Atoms Derived from Small-Angle X-ray Scattering Interferometry for Analyzing a Conformational Ensemble of Heavy Atom-Labeled Small Molecules

To obtain unbiased information about the dynamic conformational ensemble of a molecule in solution, one promising approach is small-angle X-ray scattering (SAXS). Conventionally, SAXS data are converted to a pair distribution function, which describes the distance distribution between all pairs of atoms within a molecule. If two strong X-ray scatterers are introduced and the background contributions from the other atoms are suppressed, then the distance distribution between the two scatterers provides spatial information about a flexible molecule. Gold nanocrystals can provide such information for distances of >50 Å. Here, we synthesized a chemical compound containing two iodine atoms attached to the ends of a flexible polyethylene glycol chain and used the relevant singly labeled and unlabeled compounds to suppress the background contribution. This is a feasibility demonstration to prove that the distance distribution in the range of 10-30 Å can be experimentally accessed by SAXS.

Spatio-temporal heterogeneities in nanosegregated single-molecule polymeric nanoparticles

The study of the coupling between structural and dynamical heterogeneities in nanostructured systems is essential for the design of hybrid materials with the desired properties. Here, we use atomistic molecular dynamics simulations to closely examine the dynamical heterogeneities in nanostructured single-molecule nanoparticles consisting of mikto-arm star copolymers with poly(ethylene oxide), PEO, and polystyrene, PS, arms. The particles exhibit an internally nanostructured morphology, resembling either "Janus-like" or "patchy-like" morphology when the functionality of the stars varies. The differences in the local environment result in strong intramolecular dynamical heterogeneities. In the proximity of the star core, geometric constraints promote unfavorable PEO:PS contacts that lead to a behavior similar to dynamically asymmetric miscible polymer blends or disordered copolymers. In contrast, further away from the core, the nanosegregation induces segmental dynamics very similar to the one found in the homopolymer star analogues.

Dynamic Heterogeneity in Ring-Linear Polymer Blends

We present results from a direct statistical analysis of long molecular dynamics (MD) trajectories for the orientational relaxation of individual ring molecules in blends with equivalent linear chains. Our analysis reveals a very broad distribution of ring relaxation times whose width increases with increasing ring/linear molecular length and increasing concentration of the blend in linear chains. Dynamic heterogeneity is also observed in the pure ring melts but to a lesser extent. The enhanced degree of dynamic heterogeneity in the blends arises from the substantial increase in the intrinsic timescales of a large subpopulation of ring molecules due to their involvement in strong threading events with a certain population of the linear chains present in the blend. Our analysis suggests that the relaxation dynamics of the rings are controlled by the different states of their threading by linear chains. Unthreaded or singly-threaded rings exhibit terminal relaxation very similar to that in their own melt, but multiply-threaded rings relax much slower due to the long lifetimes of the corresponding topological interactions. By further analyzing the MD data for ring molecule terminal relaxation in terms of the sum of simple exponential functions we have been able to quantify the characteristic relaxation times of the corresponding mechanisms contributing to ring relaxation both in their pure melts and in the blends, and their relative importance. The extra contribution due to ring-linear threadings in the blends becomes immediately apparent through such an analysis.

Investigation of the properties of nanographene in polymer nanocomposites through molecular simulations: dynamics and anisotropic Brownian motion

The dynamical behavior of nanographene sheets dispersed in polymer matrices is investigated through united-atom molecular dynamics simulations. The Brownian motion of the sheet and the anisotropy in its translational and orientational diffusion are the topics of the current study. Different polymer matrices and pristine and functionalized graphene constitute various nanocomposite systems. Interactions between the nanographene flake and the matrix determine the dynamics of the systems. The dynamics is reduced in polyethylene oxide compared to polyethylene matrix, whereas carboxylated sheets move considerably slower than the pristine nanographene in any matrix. Diffusion is anisotropic for short times, while it becomes isotropic in the long time limit. The in-plane motion of the nanographene sheet is faster than the out-of-plane component, in agreement with the diffusion of perfectly oblate ellipsoids. In functionalized graphene, the anisotropy is suppressed. By exploring the temperature effect on both the nanographene sheet and polymer close to the surface, indications for coupling in the motion of the two components are revealed. The strong effect of edge functional groups on the dynamics can be used as a way to control the Brownian motion of nanographene sheets in polymer nanocomposites and consequently tailor the properties of the materials.

Shear Rheology of Unentangled and Marginally Entangled Ring Polymer Melts from Large-Scale Nonequilibrium Molecular Dynamics Simulations

We present results for the steady state shear rheology of non-concatenated, unentangled and marginally entangled ring poly(ethylene oxide) (PEO) melts from detailed, atomistic nonequilibrium molecular dynamics (NEMD) simulations, and compare them to the behavior of the corresponding linear melts. The applied flow field spans a wide range of shear rates, from the linear (Newtonian) to the highly non-linear (described by a power law) regime. For all melts studied, rings are found to exhibit shear thinning but to a lesser degree compared to linear counterparts, mostly due to their reduced deformability and stronger resistance to alignment in the direction of flow. These features are attributed to the more compact structure of ring molecules compared to linear chains; the latter are capable of adopting wider and more open conformations even under shear due to the freedom provided by the free ends. Similar to linear melts, rings also exhibit a first and a second normal stress coefficient; the latter is negative. The ratio of the magnitude of the two coefficients remains practically constant with shear rate and is systematically higher than the corresponding one for linear melts. Emphasis was also given to the statistics of terminal (re-orientational) relaxation times which we computed by analyzing all chains in the simulated systems one by one; it was demonstrated that long time dynamics are strongly heterogeneous both for rings and (especially) linears. Repeating the analysis under flow conditions, and as expected, we found that the applied flow field significantly suppresses dynamic heterogeneity, especially for high shear rates well beyond the Newtonian plateau. Finally, a detailed geometrical analysis revealed that the average population of ring-ring threading events in the longest melt studied here (the PEO-5k ring) remains practically unaffected by the imposed flow rate even at strong shear rates, except for multi-threadings which disappear. To further analyze this peculiar and rather unexpected effect, we computed the corresponding survival times and penetration lengths, and found that the overwhelming majority of threadings under shear are extremely weak constraints, as they are characterized by very small penetration lengths, thus also by short survival times. They are expected therefore to play only a minor (if any) role on chain dynamics.

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.

Nanostructuring Single-Molecule Polymeric Nanoparticles via Macromolecular Architecture

Heterogeneous polymer-based nanoparticles comprise a very promising family of materials for a broad range of applications. Here, we present a detailed study of structural heterogeneities in nanostructured single-molecule nanoparticles in various environments by means of atomistic molecular dynamics simulations. The nanoparticles consist of mikto-arm star copolymers with two types of chemically incompatible arms, namely poly(ethylene oxide) (PEO) and polystyrene (PS), (PS) _n,(PEO) _n, where n is the number of arms. The immiscibility between the two components gives rise to intramolecularly nanostructured particles. The nanostructured objects resemble either "Janus-like" or "patchy-like" particles, depending on the number or the length of the arms (or both) as well as the interaction with the surrounding medium. The degree of intramolecular heterogeneity increases with increasing number of arms and with decreasing affinity of star components to the polymer host. We provide a detailed analysis of the internal structure of the star-shaped particles, focusing on the intramolecular packing and the spatial arrangement of the arms. The results of our study can be used to design heterogeneous, internally nanostructured particles with two phases of distinct static properties for challenging specific applications of next-generation materials.

Effect of macromolecular architecture on the self-assembly behavior of copolymers in a selective polymer host

We present a detailed simulation study of the structural and dynamical behavior of star-shaped mikto-arm (polystyrene)8(poly(ethylene oxide))8, (PS)8(PEO)8, copolymers with eight arms of each type, versus that of a linear polystyrene-block-poly(ethylene oxide), PS-b-PEO, diblock, in a selective homopolymer host. Both copolymers are blended at the same weight fraction 33% with an oligomeric PEO host. We use atomistic molecular dynamics simulations to account for the molecular interactions present in the blends and to study quantitatively the dynamical and structural properties of these systems. The presence of the selective oligomeric PEO host leads to the formation of complex self-assembled structures. While cylindrical structures are formed in the case of linear diblock copolymers, mikto-arm star copolymers form percolated interconnected assemblies within the PEO host. The cylindrical objects formed by the linear diblock copolymers exhibit a higher degree of compactness and a weaker temperature dependence than the percolated network formed by their star-shaped analogues. The dynamics is governed primarily by the local structural heterogeneity, i.e., the environment around a segment, which is determined by the interaction between the different components, the macromolecular architecture of the copolymer as well as the associated geometrical constrains. Our data further stress the fact that the structural and dynamical properties in these blends may be controlled/tuned by the macromolecular architecture of the copolymer and/or by adjusting the temperature.

Scaling Laws for the Conformation and Viscosity of Ring Polymers in the Crossover Region around Me from Detailed Molecular Dynamics Simulations

We present results from detailed, atomistic molecular dynamics (MD) simulations of pure, strictly monodisperse linear and ring poly(ethylene oxide) (PEO) melts under equilibrium and nonequilibrium (shear flow) conditions. The systems examined span the regime of molecular weights (M_w) from sub-Rouse (M_w < M_e) to reptation (M_w ∼ 10 M_e), where M_e denotes the characteristic entanglement molecular weight of linear PEO. For both PEO architectures (ring and linear), the predicted chain center-of-mass self-diffusion coefficients D_G as a function of PEO M_w are in remarkable agreement with experimental data. From the flow simulations under shear, we have extracted and analyzed the zero-shear viscosity of ring and linear PEO melts as a function of M_w.

Coarse-Grained Simulations of Aqueous Thermoresponsive Polyethers

Thermoresponsive polymers can change structure or solubility as a function of temperature. Block co-polymers of polyethers have a response that depends on polymer molecular weight and co-polymer composition. A coarse-grained model for aqueous polyethers is developed and applied to polyethylene oxide and polyethylene oxide-polypropylene oxide-polyethylene oxide triblock co-polymers. In this model, no interaction sites on hydrogen atoms are included, no Coulombic interactions are present, and all interactions are short-ranged, treated with a combination of two- and three-body terms. Our simulations find that The triblock co-polymers tend to associate at temperatures above 350 K. The aggregation is stabilized by contact between The hydrophobic methyl groups on The propylene oxide monomers and involves a large, favorable change in entropy.

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.

Molecular dynamics and a spectroscopic study of sulfur dioxide absorption by an ionic liquid and its mixtures with PEO

An investigation comprising experimental techniques (absorption capacity of SO₂ and vibrational spectroscopy) and molecular simulations (thermodynamics, structure, and dynamics) has been performed for the polymer poly(ethylene oxide) (PEO), the ionic liquid butyltrimethylammonium bis(trifluoromethylsulfonyl)imide ([N₄₁₁₁][Tf₂N]) and their mixtures as sulfur dioxide (SO₂) absorbing materials. The polymer PEO has higher capacity to absorb SO₂ than the neat ionic liquid, whereas the mixtures presented intermediary absorption capacities. The band assigned to the symmetric stretching band of SO₂ at ca. 1140 cm^-1, which is considered a spectroscopic probe for the strength of SO₂ interactions with its neighborhood, shifts to lower wavenumbers as more negative total interaction energy values of SO₂ were evaluated from the simulations. The solvation free energy of SO₂, ΔG_sol, correlates linearly with the absorption capacity of SO₂. The negative values of ΔG_sol are due to negative and positive values of enthalpy and entropy, respectively. In the ionic liquid, SO₂ weakens the cation-anion interactions, whereas in the mixture with a high content of PEO these interactions are slightly increased. Such effects were correlated with the relative population of cisoid and transoid conformers of Tf₂N anions as revealed by Raman spectroscopy. Moreover, the presence of SO₂ in the systems provokes the increase of diffusion coefficients of the absorbing species in comparison with the systems without the gas. Proper to the slow dynamics of the polymer, the diffusion coefficient of ions and SO₂ diminishes with the increase of the PEO content.

Microscopic Dynamics and Topology of Polymer Rings Immersed in a Host Matrix of Longer Linear Polymers: Results from a Detailed Molecular Dynamics Simulation Study and Comparison with Experimental Data

We have performed molecular dynamics (MD) simulations of melt systems consisting of a small number of long ring poly(ethylene oxide) (PEO) probes immersed in a host matrix of linear PEO chains and have studied their microscopic dynamics and topology as a function of the molecular length of the host linear chains. Consistent with a recent neutron spin echo spectroscopy study (Goossen et al., Phys. Rev. Lett. 2015, 115, 148302), we have observed that the segmental dynamics of the probe ring molecules is controlled by the length of the host linear chains. In matrices of short, unentangled linear chains, the ring probes exhibit a Rouse-like dynamics, and the spectra of their dynamic structure factor resemble those in their own melt. In striking contrast, in matrices of long, entangled linear chains, their dynamics is drastically altered. The corresponding dynamic structure factor spectra exhibit a steep initial decay up to times on the order of the entanglement time τe of linear PEO at the same temperature but then they become practically time-independent approaching plateau values. The plateau values are different for different wavevectors; they also depend on the length of the host linear chains. Our results are supported by a geometric analysis of topological interactions, which reveals significant threading of all ring molecules by the linear chains. In most cases, each ring is simultaneously threaded by several linear chains. As a result, its dynamics at times longer than a few τe should be completely dictated by the release of the topological restrictions imposed by these threadings (interpenetrations). Our topological analysis did not indicate any effect of the few ring probes on the statistical properties of the network of primitive paths of the host linear chains.

The Differential Response of Proteins to Macromolecular Crowding

The habitat in which proteins exert their function contains up to 400 g/L of macromolecules, most of which are proteins. The repercussions of this dense environment on protein behavior are often overlooked or addressed using synthetic agents such as poly(ethylene glycol), whose ability to mimic protein crowders has not been demonstrated. Here we performed a comprehensive atomistic molecular dynamic analysis of the effect of protein crowders on the structure and dynamics of three proteins, namely an intrinsically disordered protein (ACTR), a molten globule conformation (NCBD), and a one-fold structure (IRF-3) protein. We found that crowding does not stabilize the native compact structure, and, in fact, often prevents structural collapse. Poly(ethylene glycol) PEG500 failed to reproduce many aspects of the physiologically-relevant protein crowders, thus indicating its unsuitability to mimic the cell interior. Instead, the impact of protein crowding on the structure and dynamics of a protein depends on its degree of disorder and results from two competing effects: the excluded volume, which favors compact states, and quinary interactions, which favor extended conformers. Such a viscous environment slows down protein flexibility and restricts the conformational landscape, often biasing it towards bioactive conformations but hindering biologically relevant protein-protein contacts. Overall, the protein crowders used here act as unspecific chaperons that modulate the protein conformational space, thus having relevant consequences for disordered proteins.

Most in vitro and in silico biophysical experiments generally study proteins in an isolated environment, overlooking that their natural environment—the cell cytoplasm—is a solution that is highly populated by proteins. To address this knowledge gap, here we explored how a crowded environment alters the conformational sampling of three proteins, each with a different degree of disorder and flexibility. We simulated a crowded system composed by the three proteins and reaching a cell-like concentration and compared the protein behavior observed with that induced by PEG500, a synthetic crowding agent. Despite some similarities between the environments, protein crowders showed a number of characteristics that raise concerns about the use of diluted solutions or synthetic agents when studying protein behavior.

Mid-infrared spectroscopic investigation of the perfect vitrification of poly(ethylene glycol) aqueous solutions

Crystallization/recrystallization behaviors of poly(ethylene glycol) (PEG) aqueous solutions with water contents (WC's) of ∼36-51 wt % were investigated by temperature-variable mid-infrared spectroscopy. At a WC of 43.2 wt %, crystallization and recrystallization of water and PEG were not observed. At this specific WC value (WCPV), perfect vitrification occurred. Below and above the WCPV value, crystallization/recrystallization behaviors changed drastically. The crystallization temperature below WCPV (237 K) was ∼10 K greater than that above WCPV (226 K). Recrystallization above and below WCPV occurred in one (213 K) and two (198 and 210 K) steps, respectively. These findings resulted from the difference in the (re)crystallization behaviors of water molecules associated with PEG chains with helical and random-coil conformations. These two types of water molecules might have limiting concentrations for their (re)crystallization, indicating that perfect vitrification might have occurred when the concentrations of the two types of water molecules were less than the limiting concentrations of their (re)crystallization.

Atomistic potentials for trisiloxane, alkyl ethoxylate, and perfluoroalkane-based surfactants with TIP4P/2005 and application to simulations at the air-water interface

The mechanism of superspreading, the greatly enhanced spreading of water droplets facilitated by trisiloxane surfactants, is still under debate, largely because the role and behavior of the surfactants cannot be sufficiently resolved by experiments or continuum simulations. Previous molecular dynamics studies have been performed with simple model molecules or inaccurate models, strongly limiting their explanatory power. Here we present a force field dedicated to superspreading, extending existing quantum-chemistry-based models for the surfactant and the TIP4P/2005 water model ( Abascal et al. J. Chem. Phys. , 2005 , 123 , 234505 ). We apply the model to superspreading trisiloxane surfactants and nonsuperspreading alkyl ethoxylate and perfluoroalkane surfactants at various concentrations at the air-water interface. We show that the developed model accurately predicts surface tensions, which are typically assumed important for superspreading. Significant differences between superspreading and traditional surfactants are presented and their possible relation to superspreading discussed. Although the force field has been developed for superspreading problems, it should also perform well for other simulations involving polymers or copolymers with water.