Computer simulations of amorphous polymers: From quantum mechanical calculations to mesoscopic models

Doped photo-crosslinked polyesteramide hydrogels as solid electrolytes for supercapacitors

High-performance hydrogels play a crucial role as solid electrolytes for flexible electrochemical supercapacitors (ESCs). More specifically, all solid-state ESCs based on renewable, biodegradable and/or biocompatible hydrogels doped with inorganic salts as electrolytes are attractive not only because of their contribution to reducing resource consumption and/or the generation of electronic garbage, but also due to their potential applicability in the biomedical field. Here, computer simulations have been combined with experimental measurements to probe the outstanding capability as solid electrolytes of photo-crosslinked unsaturated polyesteramide hydrogels containing phenylalanine, butenediol and fumarate, and doped with NaCl (UPEA-Phe/NaCl). Atomistic molecular dynamics simulations have shown the influence of the hydrogel pore structure in the migration of Na⁺ and Cl^- ions, suggesting that UPEA-Phe/NaCl hydrogels prepared without completing the photo-crosslinking reaction will exhibit better behavior as solid electrolytes. Theoretical predictions have been confirmed by potentiodynamic and galvanostatic studies on ESCs fabricated using poly(3,4-ethylenedioxythiophene) electrodes and UPEA-Phe/NaCl hydrogels, which were obtained using different times of exposure to UV radiation (i.e. 4 and 8 h for incomplete and complete photo-crosslinking reaction). Moreover, the behavior as a solid electrolyte of the UPEA-Phe/NaCl hydrogel prepared using a photo-polymerization time of 4 h has been found to be significantly superior to those exhibited by different polypeptide and polysaccharide hydrogels, which were analyzed using ESCs with identical electrodes and experimental conditions.

Atomistic simulations of the structure of highly crosslinked sulfonated poly(styrene-co-divinylbenzene) ion exchange resins

The microscopic structures of highly crosslinked sulfonated poly(styrene-co-divinylbenzene) resins have been modeled by generating atomistic microstructures using stochastic-like algorithms, which are subsequently relaxed using molecular dynamics. Two different generation algorithms have been tested. The relaxation of the microstructures generated by the first algorithm, which is based on a homogeneous construction of the resin, leads to a significant overestimation of the experimental density as well as to an unsatisfactory description of the porosity. In contrast, the generation approach that combines algorithms for the heterogeneous growing and branching of the chains enables the formation of crosslinks with different topologies. In particular, the intrinsic heterogeneity observed in these resins is efficiently reproduced when the topological loops, which are defined by two or more crosslinks closing a cycle, are present in their microscopic description. Thus, the apparent density, porosity and pore volume estimated using microstructures with these topological loops, called super-crosslinks, are in very good agreement with the experimental results. Although the backbone dihedral angle distribution of the generated and relaxed models is not influenced by the topology, the number and type of crosslinks affect the medium- and long-range atomic disposition of the backbone atoms and the distribution of sulfonic groups. An analysis of the distribution of the local density indicates that super-crosslinks are responsible for the heterogeneous homogenization observed during the MD relaxation. Finally the π-π stacking interactions have been analyzed. Results indicate that those in which the two rings adopt a T-shaped disposition are considerably more abundant as compared to those with the co-facially oriented rings, independently of the resin topology.

Spontaneous crystallization in athermal polymer packings

We review recent results from extensive simulations of the crystallization of athermal polymer packings. It is shown that above a certain packing density, and for sufficiently long simulations, all random assemblies of freely-jointed chains of tangent hard spheres of uniform size show a spontaneous transition into a crystalline phase. These polymer crystals adopt predominantly random hexagonal close packed morphologies. An analysis of the local environment around monomers based on the shape and size of the Voronoi polyhedra clearly shows that Voronoi cells become more spherical and more symmetric as the system transits to the ordered state. The change in the local environment leads to an increase in the monomer translational contribution to the entropy of the system, which acts as the driving force for the phase transition. A comparison of the crystallization of hard-sphere polymers and monomers highlights similarities and differences resulting from the constraints imposed by chain connectivity.

Systematic coarse-graining of the dynamics of entangled polymer melts: the road from chemistry to rheology

For optimal processing and design of entangled polymeric materials it is important to establish a rigorous link between the detailed molecular composition of the polymer and the viscoelastic properties of the macroscopic melt. We review current and past computer simulation techniques and critically assess their ability to provide such a link between chemistry and rheology. We distinguish between two classes of coarse-graining levels, which we term coarse-grained molecular dynamics (CGMD) and coarse-grained stochastic dynamics (CGSD). In CGMD the coarse-grained beads are still relatively hard, thus automatically preventing bond crossing. This also implies an upper limit on the number of atoms that can be lumped together (up to five backbone carbon atoms) and therefore on the longest chain lengths that can be studied. To reach a higher degree of coarse-graining, in CGSD many more atoms are lumped together (more than ten backbone carbon atoms), leading to relatively soft beads. In that case friction and stochastic forces dominate the interactions, and action must be undertaken to prevent bond crossing. We also review alternative methods that make use of the tube model of polymer dynamics, by obtaining the entanglement characteristics through a primitive path analysis and by simulation of a primitive chain network. We finally review super-coarse-grained methods in which an entire polymer is represented by a single particle, and comment on ways to include memory effects and transient forces.

Combined molecular algorithms for the generation, equilibration and topological analysis of entangled polymers: methodology and performance

We review the methodology, algorithmic implementation and performance characteristics of a hierarchical modeling scheme for the generation, equilibration and topological analysis of polymer systems at various levels of molecular description: from atomistic polyethylene samples to random packings of freely-jointed chains of tangent hard spheres of uniform size. Our analysis focuses on hitherto less discussed algorithmic details of the implementation of both, the Monte Carlo (MC) procedure for the system generation and equilibration, and a postprocessing step, where we identify the underlying topological structure of the simulated systems in the form of primitive paths. In order to demonstrate our arguments, we study how molecular length and packing density (volume fraction) affect the performance of the MC scheme built around chain-connectivity altering moves. In parallel, we quantify the effect of finite system size, of polydispersity, and of the definition of the number of entanglements (and related entanglement molecular weight) on the results about the primitive path network. Along these lines we approve main concepts which had been previously proposed in the literature.