Abstract:We present atomistic molecular dynamics simulations of two Polyethylene systems where all entanglements are trapped: a perfect network, and a melt with grafted chain ends. We examine microscopically at what level topological constraints can be considered as a collective entanglement effect, as in tube model theories, or as certain pairwise uncrossability interactions, as in slip-link models.A pairwise parameter, which varies between these limiting cases, shows that, for the systems studied, the character of the entanglement environment is more pairwise than collective. We employ a novel methodology, which analyzes entanglement constraints into a complete set of pairwise interactions, similar to slip links. Entanglement confinement is assembled by a plethora of links, with a spectrum of confinement strengths, from strong to weak. The strength of interactions is quantified through a link 'persistence', which is the fraction of time for which the links are active. By weighting links according to their strength, we show that confinement is imposed mainly by the strong ones, and that the weak, trapped, uncrossability interactions cannot contribute to the low frequency modulus of an elastomer, or the plateau modulus of a melt. A selfconsistent scheme for mapping topological constraints to specific, strong binary links, according to a given entanglement density, is proposed and validated. Our results demonstrate that slip links can be viewed as the strongest pairwise interactions of a collective entanglement environment. The methodology developed provides a basis for bridging the gap between atomistic simulations and mesoscopic slip link models.

Models of the mechanical response of glassy polymers commonly assume that entanglements inherited from the melt act like chemical cross-links. The predictions from these network models and the physical picture they are based on are tested by following the evolution of topological constraints in simulations of model polymer glasses. The same behavior is observed for polymers with entanglement lengths N e that vary by a factor of 3. A prediction for the craze extension ratio Λ based on the network model describes trends with N e , but polymers do not have the taut configurations it assumes. There is also no evidence of the predicted geometrically necessary entanglement loss. While the number of entanglements remains constant, the identity of the chains forming constraints changes. The same relation between the amount of entanglement exchange and nonaffine displacement of monomers is found for crazing and thermal diffusion in end-constrained melts. In both cases, about 1/ 3 of the constraints change when monomers move by the tube radius. The results show that chains in deformed glassy polymers are constrained by their rheological tubes rather than by entanglements that act like discrete cross-links.

Surfactants are amphiphilic molecules with multiple uses and industrial applications as detergents, wetting agents, emulsifiers, and so forth. They can be divided into three main categories: nonionic, ionic, and zwitterionic. The development of a universal computational framework able to predict key properties such as their critical micelle concentration (cmc) and the size of the micelles they form and to ultimately extract phase diagrams for their aqueous solutions, possibly in the presence of salts and oils, using their chemical constitution as input, would provide valuable information for the design and the production of these materials. In this work, we focus on ionic surfactants and investigate a possible route toward the development of such a framework based on coarse-grained simulations using the MARTINI forcefield in two versions: its implicit solvent version, called Dry MARTINI, and its explicit solvent version, called Wet MARTINI. The surfactants considered in our efforts are the anionic sodium dodecyl sulfate (SDS) and the three cationic cetyl, dodecyl, and octyl trimethyl ammonium bromide (CTAB, DTAB, and OTAB, respectively). First, we choose their mapping onto coarse-grained MARTINI beads. Next, we estimate their cmc’s, their peak aggregation numbers, N agg, and in the case of SDS, its small angle neutron scattering pattern at low concentrations, applying the Dry MARTINI forcefield. With a single modification to the Lennard-Jones energy parameter between hydrophobic beads and invoking Ewald summation with a physically meaningful dielectric constant for electrostatic interactions, our estimates are in very good agreement with experimental results. Furthermore, we predict the phase behavior of SDS/water and CTAB/water binary solutions using Wet MARTINI and find semiquantitative agreement with experimental phase diagrams. We conclude that the MARTINI forcefield, with careful treatment of electrostatic interactions and appropriate modification of parameters for some key functional groups, can be a powerful ally in the quest for a universal computational framework for the design of new surfactants with improved properties.

No abstract

Tacticity is one of the most important features of stereoregular polymers, with significant impact on morphology and on a variety of properties such as conformational, thermal, rheological, mechanical, etc. In this work we focus on tacticity effects on the conformations and more specifically on the unperturbed dimensions of single polypropylene (PP) homopolymer and both block and random poly(ethylene−propylene) copolymer chains. The equilibration of all chains is achieved by applying the single-chain Monte Carlo algorithm of Tzounis et al. [Macromolecules 2017, 50, 4575−4587]. In agreement with past studies, we find that tacticity has a significant effect on the stiffness of PP homopolymer. The characteristic ratio exhibits a nonmonotonic dependence on the fraction of meso dyads along the PP chains, which results from two competing mechanisms, revealed by analysis of the torsional states of the PP chain backbones. A simple theoretical model is developed to describe the dependence of the stiffness of poly(ethylene-block-propylene) copolymer on its propylene content and on the tacticity of its PP block. Finally, the effect of both tacticity and propylene content on the stiffness of poly(ethylene-randompropylene) chains, used to model ethylene−propylene monomer (EPM) materials, is found to be in very good qualitative agreement with available rotational isomeric state (RIS) model predictions.