Structured Ionomer Thin Films at Water Interface: Molecular Dynamics Simulation Insight.

Controlling the structure and dynamics of thin films of ionizable polymers at water interfaces is critical to their many applications. As the chemical diversity within one polymer is increased, controlling the structure and dynamics of the polymer, which is a key to their use, becomes a challenge. Here molecular dynamics simulations (MD) are used to obtain molecular insight into the structure and dynamics of thin films of one such macromolecule at the interface with water. The polymer consists of an ABCBA topology with randomly sulfonated polystyrene (C), tethered symmetrically to flexible poly(ethylene-r-propylene) blocks (B), and end-capped by a poly(t-butylstyrene) block (A). The compositions of the interfacial and bulk regions of thin films of the ABCBA polymers are followed as a function of exposure time to water. We find that interfacial rearrangements take place where buried ionic segments migrate toward the water interface. The hydrophobic blocks collapse and rearrange to minimize their exposure to water. The water that initially drives interfacial reengagements breaks the ionic clusters within the film, forming a dynamic hydrophilic internal network within the hydrophobic segments.

Coarse-Grained Modeling of Polyethylene Melts: Effect on Dynamics.

The distinctive viscoelastic behavior of polymers results from a coupled interplay of motion on multiple length and time scales. Capturing the broad time and length scales of polymer motion remains a challenge. Using polyethylene (PE) as a model macromolecule, we construct coarse-grained (CG) models of PE with three to six methyl groups per CG bead and probe two critical aspects of the technique: pressure corrections required after iterative Boltzmann inversion (IBI) to generate CG potentials that match the pressure of reference fully atomistic melt simulations and the transferability of CG potentials across temperatures. While IBI produces nonbonded pair potentials that give excellent agreement between the atomistic and CG pair correlation functions, the resulting pressure for the CG models is large compared with the pressure of the atomistic system. We find that correcting the potential to match the reference pressure leads to nonbonded interactions with much deeper minima and slightly smaller effective bead diameter. However, simulations with potentials generated by IBI and pressure-corrected IBI result in similar mean-square displacements (MSDs) and stress autocorrelation functions G(t) for PE melts. While the time rescaling factor required to match CG and atomistic models is the same for pressure- and non-pressure-corrected CG models, it strongly depends on temperature. Transferability was investigated by comparing the MSDs and stress autocorrelation functions for potentials developed at different temperatures.

Cluster Morphology-Polymer Dynamics Correlations in Sulfonated Polystyrene Melts: Computational Study.

Reaching exceptionally long times up to 500 ns in equilibrium and nonequilibrium molecular dynamics simulations studies, we have attained a fundamental molecular understanding of the correlation of ionomer clusters structure and multiscale dynamics, providing new insight into one critical, long-standing challenge in ionic polymer physics. The cluster structure in melts of sulfonated polystyrene with Na^{+} and Mg^{2+} counterions are resolved and correlated with the dynamics on multiple length and time scales extracted from measurements of the dynamic structure factor and shear rheology. We find that as the morphology of the ionic clusters changes from ladderlike for Na^{+} to disordered structures for Mg^{2+}, the dynamic structure factor is affected on the length scale corresponding to the ionic clusters. Rheology studies show that the viscosity for Mg^{2+} melts is higher than for Na^{+} ones for all shear rates, which is well correlated with the larger ionic clusters' size for the Mg^{2+} melts.

Resolving Dynamic Properties of Polymers through Coarse-Grained Computational Studies.

Coupled length and time scales determine the dynamic behavior of polymers and underlie their unique viscoelastic properties. To resolve the long-time dynamics it is imperative to determine which time and length scales must be correctly modeled. Here we probe the degree of coarse graining required to simultaneously retain significant atomistic details and access large length and time scales. The degree of coarse graining in turn sets the minimum length scale instrumental in defining polymer properties and dynamics. Using linear polyethylene as a model system, we probe how the coarse-graining scale affects the measured dynamics. Iterative Boltzmann inversion is used to derive coarse-grained potentials with 2-6 methylene groups per coarse-grained bead from a fully atomistic melt simulation. We show that atomistic detail is critical to capturing large-scale dynamics. Using these models we simulate polyethylene melts for times over 500 µs to study the viscoelastic properties of well-entangled polymer melts.

Clustering effects in ionic polymers: Molecular dynamics simulations.

Ionic clusters control the structure, dynamics, and transport in soft matter. Incorporating a small fraction of ionizable groups in polymers substantially reduces the mobility of the macromolecules in melts. These ionic groups often associate into random clusters in melts, where the distribution and morphology of the clusters impact the transport in these materials. Here, using molecular dynamic simulations we demonstrate a clear correlation between cluster size and morphology with the polymer mobility in melts of sulfonated polystyrene. We show that in low dielectric media ladderlike clusters that are lower in energy compared with spherical assemblies are formed. Reducing the electrostatic interactions by enhancing the dielectric constant leads to morphological transformation from ladderlike clusters to globular assemblies. Decrease in electrostatic interaction significantly enhances the mobility of the polymer.

Phases in the (1 - x)Na0.5Bi0.5TiO3-(x)CaTiO3 system.

The structures of (1 - x)Na(0.5)Bi(0.5)TiO(3)-(x)CaTiO(3) at room temperature have been investigated using neutron powder diffraction and dielectric studies. The system exhibits an orthorhombic (Pbnm) structure for x ≥ 0.15 and rhombohedral (R3c) for x ≤ 0.05. For x = 0.10, though the neutron diffraction pattern shows features of the orthorhombic (Pbnm) structure, Rietveld refinement using this structure shows a drastic reduction in the in-phase tilt angle (∼4°) as compared to the corresponding value (∼8°) for a neighbouring composition x = 0.15. The neutron diffraction pattern of x = 0.10 could be fitted equally well using a two-phase model (R3c + Pbnm) with orthorhombic as the minor phase (22%), without the need for a drastic decrease in the in-phase tilt angle. The dielectric studies of x = 0.10 revealed the presence of the polar R3c phase, thereby favouring the phase coexistence model, instead of a single-phase Pbnm structure, for this composition.