Bead–Spring Simulation of Ionomer Melts—Studying the Effects of Chain-Length and Associating Group Fraction on Equilibrium Structure and Extensional Flow Behavior

Ionomers are associative polymers with diverse applications ranging from selective membranes and high-performance adhesives to abrasion- and chemical-resistant coatings, insulation layers, vacuum packaging, and foamed sheets. Within equilibrium melt, the ionic or associating groups are known to form...

Full description

Bibliographic Details
Main Authors: Supun S. Mohottalalage, Andrew P. Saab, Amitesh Maiti
Format: Article
Language:English
Published: MDPI AG 2023-11-01
Series:Polymers
Subjects:
Online Access:https://www.mdpi.com/2073-4360/15/23/4560
Description
Summary:Ionomers are associative polymers with diverse applications ranging from selective membranes and high-performance adhesives to abrasion- and chemical-resistant coatings, insulation layers, vacuum packaging, and foamed sheets. Within equilibrium melt, the ionic or associating groups are known to form thermally reversible, associative clusters whose presence can significantly affect the system’s mechanical, viscoelastic, and transport properties. It is, thus, of great interest to understand how to control such clusters’ size distribution, shape, and stability through the designed choice of polymer architecture and the ionic groups’ fraction, arrangement, and interaction strength. In this work, we represent linear associating polymers using a Kremer–Grest type bead–spring model and perform large-scale MD simulations to explore the effect of polymer chain-length (<inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>l</mi></mrow></semantics></math></inline-formula>) and fraction (<inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><msub><mrow><mi>f</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow></semantics></math></inline-formula>) of randomly placed associating groups on the size distribution and stability of formed clusters. We consider different chain-lengths (below and above entanglement), varying fractions of associating groups (represented by ‘sticky’ beads) between 5 and 20%, and a fixed sticky–sticky nonbond interaction strength of four times that between regular non-associating beads. For all melts containing associating groups the equilibrium structure factor <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>S</mi><mo>(</mo><mi>q</mi><mo>)</mo></mrow></semantics></math></inline-formula> displays a signature ionomer peak at low wave vector <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>q</mi></mrow></semantics></math></inline-formula> whose intensity increases with increasing <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><msub><mrow><mi>f</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow></semantics></math></inline-formula> and <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>l</mi></mrow></semantics></math></inline-formula>. The average cluster size <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi>c</mi></mrow></msub></mrow></semantics></math></inline-formula> increases with <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><msub><mrow><mi>f</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow></semantics></math></inline-formula>. However, the effect of chain-length on <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi>c</mi></mrow></msub></mrow></semantics></math></inline-formula> appears to be pronounced only at higher values of <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><msub><mrow><mi>f</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow></semantics></math></inline-formula>. Under extensional flows, the computed stress (and viscosity) is higher at higher <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><msub><mrow><mi>f</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow></semantics></math></inline-formula> and <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><mi>l</mi></mrow></semantics></math></inline-formula> regardless of strain rate. Beyond a critical strain rate, we observe fragmentation of the associative clusters, which has interesting effects on the stress/viscous response.
ISSN:2073-4360