Performance of Coarse Graining in Estimating Polymer Properties: Comparison with the Atomistic Model

Temperature dependence of aging dynamics in highly non-equilibrium model polymer glasses

A central feature of the non-equilibrium glassy “state” is its tendency to age toward equilibrium, obeying signatures identified by Kovacs over 50 years ago. The origin of these signatures, their fate far from equilibrium and at high temperatures, and the underlying nature of the glassy “state” far from equilibrium remain unsettled. Here, we simulate physical aging of polymeric glasses, driven much farther from equilibrium and at much higher temperatures than possible in experimental melt-quenched glasses. While these glasses exhibit Kovacs’ signatures of glassy aging at sufficiently low temperatures, these signatures disappear above the onset TA of non-Arrhenius equilibrium dynamics, suggesting that TA demarcates an upper bound to genuinely glassy states. Aging times in glasses after temperature up-jumps are found to obey an Arrhenius law interpolating between equilibrium dynamics at TA and at the start of the temperature up-jump, providing a zero-parameter rule predicting their aging behavior and identifying another unrecognized centrality of TA to aging behavior. This differs qualitatively from behavior of our glasses produced by temperature down-jumps, which exhibit a fractional power law decoupling relation with equilibrium dynamics. While the Tool–Narayanaswamy–Moynihan model can predict the qualitative single-temperature behavior of these systems, we find that it fails to predict the disappearance of Kovacs signatures above TA and the temperature dependence of aging after large temperature up-jumps. These findings highlight a need for new theoretical insights into the aging behavior of glasses at ultra-high fictive temperatures and far from equilibrium.

Ionic Polymer Nanocomposites Subjected to Uniaxial Extension: A Nonequilibrium Molecular Dynamics Study

We explore the behavior of coarse-grained ionic polymer nanocomposites (IPNCs) under uniaxial extension up to 800% strain by means of nonequilibrium molecular dynamics simulations. We observe a simultaneous increase of stiffness and toughness of the IPNCs upon increasing the engineering strain rate, in agreement with experimental observations. We reveal that the excellent toughness of the IPNCs originates from the electrostatic interaction between polymers and nanoparticles, and that it is not due to the mobility of the nanoparticles or the presence of polymer-polymer entanglements. During the extension, and depending on the nanoparticle volume fraction, polymer-nanoparticle ionic crosslinks are suppressed with the increase of strain rate and electrostatic strength, while the mean pore radius increases with strain rate and is altered by the nanoparticle volume fraction and electrostatic strength. At relatively low strain rates, IPNCs containing an entangled matrix exhibit self-strengthening behavior. We provide microscopic insight into the structural, conformational properties and crosslinks of IPNCs, also referred to as polymer nanocomposite electrolytes, accompanying their unusual mechanical behavior.

Development of Coarse-Grained Models for Poly(4-vinylphenol) and Poly(2-vinylpyridine): Polymer Chemistries with Hydrogen Bonding

In this paper, we identify the modifications needed in a recently developed generic coarse-grained (CG) model that captured directional interactions in polymers to specifically represent two exemplary hydrogen bonding polymer chemistries-poly(4-vinylphenol) and poly(2-vinylpyridine). We use atomistically observed monomer-level structures (e.g., bond, angle and torsion distribution) and chain structures (e.g., end-to-end distance distribution and persistence length) of poly(4-vinylphenol) and poly(2-vinylpyridine) in an explicitly represented good solvent (tetrahydrofuran) to identify the appropriate modifications in the generic CG model in implicit solvent. For both chemistries, the modified CG model is developed based on atomistic simulations of a single 24-mer chain. This modified CG model is then used to simulate longer (36-mer) and shorter (18-mer and 12-mer) chain lengths and compared against the corresponding atomistic simulation results. We find that with one to two simple modifications (e.g., incorporating intra-chain attraction, torsional constraint) to the generic CG model, we are able to reproduce atomistically observed bond, angle and torsion distributions, persistence length, and end-to-end distance distribution for chain lengths ranging from 12 to 36 monomers. We also show that this modified CG model, meant to reproduce atomistic structure, does not reproduce atomistically observed chain relaxation and hydrogen bond dynamics, as expected. Simulations with the modified CG model have significantly faster chain relaxation than atomistic simulations and slower decorrelation of formed hydrogen bonds than in atomistic simulations, with no apparent dependence on chain length.

Polymer Conformations, Entanglements and Dynamics in Ionic Nanocomposites: A Molecular Dynamics Study

We investigate nanoparticle (NP) dispersion, polymer conformations, entanglements and dynamics in ionic nanocomposites. To this end, we study nanocomposite systems with various spherical NP loadings, three different molecular weights, two different Bjerrum lengths, and two types of charge-sequenced polymers by means of molecular dynamics simulations. NP dispersion can be achieved in either oligomeric or entangled polymeric matrices due to the presence of electrostatic interactions. We show that the overall conformations of ionic oligomer chains, as characterized by their radii of gyration, are affected by the presence and the amount of charged NPs, while the dimensions of charged entangled polymers remain unperturbed. Both the dynamical behavior of polymers and NPs, and the lifetime and amount of temporary crosslinks, are found to depend on the ratio between the Bjerrum length and characteristic distance between charged monomers. Polymer-polymer entanglements start to decrease beyond a certain NP loading. The dynamics of ionic NPs and polymers is very different compared with their non-ionic counterparts. Specifically, ionic NP dynamics is getting enhanced in entangled matrices and also accelerates with the increase of NP loading.