Primitive chain network study on uncrosslinked and crosslinked cis-polyisoprene polymers

Phase separation dynamics in wetting ridges of polymer surfaces swollen with oils of different viscosities

When drops are placed on a sufficiently soft surface, the drop surface tension drives an out of plane deformation around the contact line (i.e., a wetting ridge). For soft elastomeric surfaces that are swollen with a liquid, capillarity from a drop can induce a phase separation in the wetting ridge. Using confocal microscopy, we study the dynamics of phase separation at the wetting ridge of glycerol drops on silicone elastomers, which are swollen with silicone oils of varying viscosity (i.e., molecular weight). We show that the viscosity of the swelling oil plays a large role in the oil separation size and separation rate. For networks swollen to near their maximum swelling (i.e., saturated), lower viscosity oil separates more and separates faster at early times compared to larger viscosity oil. During late-stage wetting, the growth rate of the separation is a function of viscosity and swelling ratio, which can be described by a simple diffusive model and a defined wetting ridge geometry. In this late-stage wetting, the higher viscosity oil evidently grows faster, likely because it is further from reaching equilibrium. Interestingly, the separated oil phase region grows with a nearly constant, geometrically similar shape. Understanding how phase separation occurs on swollen substrates should provide information on how to control drop spreading, sliding, adhesion, or friction on such surfaces.

Exploring the Role of Hydroxy- and Phosphate-Terminated cis-1,4-Polyisoprene Chains in the Formation of Physical Junction Points in Natural Rubber: Insights from Molecular Dynamics Simulations

This study elucidates the pivotal role of terminal structures in cis-1,4-polyisoprene (PI) chains, contributing to the exceptional mechanical properties of Hevea natural rubber (NR). NR's unique networking structure, crucial for crack resistance, elasticity, and strain-induced crystallization, involves two terminal groups, ω and α. The proposed ω terminal structure is dimethyl allyl-(trans-1,4-isoprene)2, and α terminals exist in various forms, including hydroxy, ester, and phosphate groups. Among others, we investigated three types of cis-1,4-PI with different terminal combinations: HPIH (pure PI with H terminal), ωPIα6 (PI with ω and α6 terminals), and ωPIPO4 (PI with ω and PO4 terminals) and revealed significant dynamics variations. Hydrogen bonds between α6 and α6 and PO4 and PO4 residues in ωPIα6 and ωPIPO4 systems induce slower dynamics of hydroxy- and phosphate-terminated PI chains. Associations between α6 and α6 and PO4 and PO4 terminals are markedly stronger than ω and ω, and hydrogen terminals in HPIH and ω PIα6,PO4 systems. Phosphate terminals exhibit a stronger mutual association than hydroxy terminals. Potentials of mean force analysis and cluster-formation-fraction computations reveal stable clusters in ωPIα6 and ωPIPO4 , supporting the formation of polar aggregates (physical junction points). Notably, phosphate terminal groups facilitate large and highly stable phosphate polar aggregates, crucial for the natural networking structure responsible for NR's outstanding mechanical properties compared to synthetic PI rubber. This comprehensive investigation provides valuable insights into the role of terminal groups in cis-1,4-PI melt systems and their profound impact on the mechanical properties of NR.

Role of Terminal Groups of cis-1,4-Polyisoprene Chains in the Formation of Physical Junction Points in Natural Rubber

The terminal structures of cis-1,4-polyisoprene (PI) chains play a vital role in the excellent comprehensive performance of Hevea natural rubber (NR) with properties such as high toughness, tear-resistance, and wet skid resistance. The cis-1,4-polyisoprene chain constituting NR exhibits a distinct composition of terminal groups comprising two distinct types, namely, the ω and α terminal groups. The structures of the ω terminal [dimethyl allyl (DMA)-(trans-1,4-isoprene)2] and six kinds of α end groups of the polymer chain of NR have been explored by utilizing a newly developed 2D NMR method. In the present work, we examine different kinds of PI melt systems, and we choose various combinations of terminal groups: Hydrogen, one DMA unit with two trans isoprene units as ω end groups and ester-terminated isopentene (α1), hydroxy-terminated isopentene (α2), ester-terminated isobutane (α3), hydroxy-terminated isobutane (α4), ester-terminated 1,4-cis-isoprene (α5), and hydroxy-terminated 1,4-cis-isoprene (α6), i.e., HPIH (PI0)-pure PI (Hydrogen terminal), ωPIα1 (PII), ωPIα2 (PIII), ωPIα3 (PIIII), ωPIα4 (PIIV), ωPIα5 (PIV), and ωPIα6 (PIVI). We evaluated dynamic and static properties of PI chains such as the end-to-end vector autocorrelation function (C(t)), its average relaxation time (τ), end-to-end distance (Ree), and radius of gyration (Rg). We also estimated the diffusion coefficients of polyisoprene chains and pair correlation functions [radial distribution functions (RDFs)], potentials of mean force (PMFs) in between end residues, and survival probability (P(τ)) of end groups around the end group by analyzing the equilibrated trajectories of full-atom MD simulations. As per the examination of C(t), rotational relaxation time τ, and RDFs, we discovered that the existence of a strong hydrogen bond in α2-α2, α4-α4, and α6-α6 residues makes the dynamics of hydroxy-terminated polyisoprene chains in ωPIα2,α4,α6 melt systems slower. From the analyses of RDFs and PMFs (W(r)), the association between [α2]-[α2], [α4]-[α4], and [α6]-[α6] terminals in ωPIα2,α4,α6 melt systems is significantly stronger than in [ISO]-[ISO] [Hydrogen terminated 1,4-cis-isoprene:(ISO)] in HPIH and ω-ω, [α1]-[α1], [α3]-[α3], and [α5]-[α5] in ωPIα1,α3,α5 systems. We quantified the fraction of cluster formation of terminal groups of a given size in the seven PI melt systems by employing the criteria of PMFs. It is revealed that no stable cluster exists in the HPIH, ωPIα1, ωPIα3, and ωPIα5 melt systems. Conversely, in the ωPIα2, ωPIα4, and ωPIα6 systems, we perceived stable clusters of [(α2)p] [(α4)p] and [(α6)p] end groups where p (2 ≤ x ≤ 6). These stable clusters validate the presence of physical junction points in between hydroxy-terminated polyisoprene chains through their α2, α4, and α6 terminals. These physical junction points might be crucial for superior properties of NR such as high toughness, crack growth resistance, and strain-induced crystallization.

Multiscale equilibration of highly entangled isotropic model polymer melts

We present a computationally efficient multiscale method for preparing equilibrated, isotropic long-chain model polymer melts. As an application, we generate Kremer-Grest melts of 1000 chains with 200 entanglements and 25 000-2000 beads/chain, which cover the experimentally relevant bending rigidities up to and beyond the limit of the isotropic-nematic transition. In the first step, we employ Monte Carlo simulations of a lattice model to equilibrate the large-scale chain structure above the tube scale while ensuring a spatially homogeneous density distribution. We then use theoretical insight from a constrained mode tube model to introduce the bead degrees of freedom together with random walk conformational statistics all the way down to the Kuhn scale of the chains. This is followed by a sequence of simulations with carefully parameterized force-capped bead-spring models, which slowly introduce the local bead packing while reproducing the larger-scale chain statistics of the target Kremer-Grest system at all levels of force-capping. Finally, we can switch to the full Kremer-Grest model without perturbing the structure. The resulting chain statistics is in excellent agreement with literature results on all length scales accessible in brute-force simulations of shorter chains.

Determination of elastic constants of functionalized graphene-based epoxy nanocomposites: a molecular modeling and MD simulation study

Functionalization of graphene is the best way to create a high degree of dispersion and bonding to polymer matrix in order to obtain high performance composites. The effects of carboxyl (-COOH) functionalized graphene (FG) on the mechanical properties of its epoxy-based nanocomposites have been examined by molecular dynamics (MD) simulations. Simulations cells of nanocomposites with varying wt% of FG (1, 2, and 3 wt%) were constructed using Material Studio 6.0. The MD simulation findings of nanocomposites reveal that they have better mechanical properties such as elastic modulus, bulk modulus, shear modulus, and the Poisson's ratio than pure epoxy. Furthermore, the computational results of nanocomposites have been effectively confirmed with available experimental data. Therefore, the current MD simulation shows a decent computational sign for the existing experimental and simulation outcomes on mechanical properties of FG/epoxy nanocomposites.

Multiscale Modeling and Simulation of Polymer Blends in Injection Molding: A Review

Modeling and simulation of the morphology evolution of immiscible polymer blends during injection molding is crucial for predicting and tailoring the products’ performance. This paper reviews the state-of-the-art progress in the multiscale modeling and simulation of injection molding of polymer blends. Technological development of the injection molding simulation on a macroscale was surveyed in detail. The aspects of various models for morphology evolution on a mesoscale during injection molding were discussed. The current scale-bridging strategies between macroscopic mold-filling flow and mesoscopic morphology evolution, as well as the pros and cons of the solutions, were analyzed and compared. Finally, a comprehensive summary of the above models is presented, along with the outlook for future research in this field.

Smart material based on boron crosslinked polymers with potential applications in cancer radiation therapy

Organoboron compounds have been playing an increasingly important role in analytical chemistry, material science, health applications, and particularly as functional polymers like boron carriers for cancer therapy. There are two main applications of boron isotopes in radiation cancer therapy, Boron Neutron Capture Therapy and Proton Boron Fusion Therapy. In this study, a novel and original material consisting of a three-dimensional polymer network crosslinked with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{10}$$\end{document}10B enriched boric acid molecules is proposed and synthesized. The effects of the exposition to thermal neutrons were studied analyzing changes in the mechanical properties of the proposed material. Dedicated Monte Carlo simulations, based on MCNP and FLUKA main codes, were performed to characterize interactions of the proposed material with neutrons, photons, and charged particles typically present in mixed fields in nuclear reactor irradiations. Experimental results and Monte Carlo simulations were in agreement, thus justifying further studies of this promising material.

Effectiveness of coarse graining degree and speedup on the dynamic properties of homopolymer

Evaluation of effective coarse graining (CG) degree and reasonable speedup relative to all-atomistic (AA) model was conducted to provide a basis for building appropriate larger-scale model. The reproducibility of atomistic conformation and temperature transferability both act as the analysis criteria to resolve the maximum acceptable CG degree. Taking short- and long time spans into account simultaneously in the estimation of computational speedup, a dynamic scaling factor is accessible in fitting mean squared displacement ratio of CG to AA as an exponential function. Computing loss in parallel running is an indispensable component in acceleration, which was also added in the evaluation. Subsequently, a quantified prediction of CG speedup arises as a multiplication of dynamic scaling factor, computing loss, time step, and the square of reduction in the number of degrees of freedom. Polyethylene oxide was adopted as a reference system to execute the direct Boltzmann inversion and iterative Boltzmann inversion. Bonded and non-bonded potentials were calculated in CG models with 1~4 monomers per bead. The effective CG degree was determined as two at the most with a speedup of four orders magnitude over AA in this study. Determination of effectiveness CG degree and the corresponding speedup prediction provide available tools in larger spatiotemporal-scale calculations.

Temperature-Dependent Gas Transport Behavior in Cross-Linked Liquid Crystalline Polyacrylate Membranes

Stable, cross-linked, liquid crystalline polymer (LCP) films for membrane separation applications have been fabricated from the mesogenic monomer 11-(4-cyanobiphenyl-4′-yloxy) undecyl methacrylate (CNBPh), non-mesogenic monomer 2-ethylhexyl acrylate (2-EHA), and cross-linker ethylene glycol dimethacrylate (EGDMA) using an in-situ free radical polymerization technique with UV initiation. The phase behavior of the LCP membranes was characterized using differential scanning calorimetry (DSC) and X-ray scattering, and indicated the formation of a nematic liquid crystalline (LC) phase above the glass transition temperature. The single gas transport behavior of CO₂, CH₄, propane, and propylene in the cross-linked LCP membranes was investigated for a range of temperatures in the LC mesophase and the isotropic phase. Solubility of the gases was dependent not only on the condensability in the LC mesophase, but also on favorable molecular interactions of penetrant gas molecules exhibiting a charge separation, such as CO₂ and propylene, with the ordered polar mesogenic side chains of the LCP. Selectivities for various gas pairs generally decreased with increasing temperature and were discontinuous across the nematic–sotropic transition. Sorption behavior of CO₂ and propylene exhibited a significant change due to a decrease in favorable intermolecular interactions in the disordered isotropic phase. Higher cross-link densities in the membrane generally led to decreased selectivity at low temperatures when the main chain motion was limited by the lack of mesogen mobility in the ordered nematic phase. However, at higher temperatures, increasing the cross-link density increased selectivity as the cross-links acted to limit chain mobility. Mixed gas permeation measurements for propylene and propane showed close agreement with the results of the single gas permeation experiments.

Slip Spring-Based Mesoscopic Simulations of Polymer Networks: Methodology and the Corresponding Computational Code

In previous work by the authors, a new methodology was developed for Brownian dynamics/kinetic Monte Carlo (BD/kMC) simulations of polymer melts. In this study, this methodology is extended for dynamical simulations of crosslinked polymer networks in a coarse-grained representation, wherein chains are modeled as sequences of beads, each bead encompassing a few Kuhn segments. In addition, the C++ code embodying these simulations, entitled Engine for Mesoscopic Simulations for Polymer Networks (EMSIPON) is described in detail. A crosslinked network of cis-1,4-polyisoprene is chosen as a test system. From the thermodynamic point of view, the system is fully described by a Helmholtz energy consisting of three explicit contributions: entropic springs, slip springs and non-bonded interactions. Entanglements between subchains in the network are represented by slip springs. The ends of the slip springs undergo thermally activated hops between adjacent beads along the chain backbones, which are tracked by kinetic Monte Carlo simulation. In addition, creation/destruction processes are included for the slip springs at dangling subchain ends. The Helmholtz energy of non-bonded interactions is derived from the Sanchez⁻Lacombe equation of state. The isothermal compressibility of the polymer network is predicted from equilibrium density fluctuations in very good agreement with the underlying equation of state and with experiment. Moreover, the methodology and the corresponding C++ code are applied to simulate elongational deformations of polymer rubbers. The shear stress relaxation modulus is predicted from equilibrium simulations of several microseconds of physical time in the undeformed state, as well as from stress-strain curves of the crosslinked polymer networks under deformation.

Prediction of the glass transition temperature and design of phase diagrams of butadiene rubber and styrene-butadiene rubber via molecular dynamics simulations

To prevent car accidents, it is important to evaluate the thermal stability of tire rubbers, such as natural rubber (NR), butadiene rubber (BR), and styrene-butadiene rubber (SBR). Controlling the glass transition temperature (T_g) is the main factor for obtaining desirable thermal stability. Here, we developed an optimized equation for the prediction of the T_g of the various rubber systems using molecular dynamics (MD) simulations. We modeled a random copolymer system, blended monomers, and calculated the T_g of butadiene isomers in each composition. From these results, we designed the T_g contour of ternary cis-trans-vinyl butadiene and derived an equation of T_g for the ternary system. Moreover, we developed an equation to evaluate the pseudo-ternary T_g of quaternary SBR and plotted it. Our results present a novel way of predicting the T_g of ternary BR and quaternary SBR, which is critical for rational tire design with optimized thermal and mechanical stability.

Importance of Compact Random Walks for the Rheology of Transient Networks

Controlling the mechanical behavior of novel supramolecular materials is of the utmost importance and requires a fundamental understanding of the underlying physical processes. We present a multimethods approach to the dynamics of entangled transient polyisoprene networks. Small-angle neutron scattering (SANS) on randomly functionalized chains shows homogeneous supramolecular melts with Gaussian chain conformations. The H-bond lifetimes (dielectric α*-process) and the rheological response in terms of the loss modulus G″ differ by 2 orders of magnitude in time. Within the concept of a compact random walk (RW), where the random walker (urazole group acting as a sticker) undergoes multiple returns to its starting point and following the concept of theoretical proposed renormalized sticky bond lifetimes, we quantitatively solve this longstanding and unexplained large discrepancy: While the bond opening gives rise to the dielectric response, for rheological relaxation the association with a new partner is relevant. This takes place only after multiple returns to the original binding partner.

Tailoring the dispersion of nanoparticles and the mechanical behavior of polymer nanocomposites by designing the chain architecture

The dispersion of nanoparticles with different polymer–nanoparticle interaction strengths and chain architectures.

Reversible wrinkles of monolayer graphene on a polymer substrate: toward stretchable and flexible electronics

The surface instability of monolayer graphene supported by a soft (polymer) substrate under equal-biaxial compression has been explored through large scale coarse-grained molecular simulations. Regardless of the interfacial adhesion strength between the graphene and the substrate, herringbone wrinkles have always been observed due to their lowest energy status, compared with the checkerboard, hexagonal, triangular and one dimensional sinusoidal modes. Moreover, the graphene-polymer substrate interaction energy has a negligible effect on the critical strain for the onset of these wrinkles. Yet, if the graphene is bonded to a rigid (non-deformable) substrate, the critical strain increases with increasing graphene-substrate interfacial strength. The surface wrinkles of graphene are delayed and suppressed by the strong bonding of graphene to the rigid substrate. Besides, only localized folds and crumples have been observed on the surface of graphene, when graphene-substrate interaction energy is strong enough. All these observations signal that the deformability (stiffness) of the substrate plays an essential role in determining the morphology of supported graphene under compression. In addition, when a flat graphene is attached on a highly pre-strained (50%) polymer substrate, wrinkles will be formed on its surface during the relaxation of pre-strain within the polymer substrate. The wrinkled graphene could be stretched up to 50% without fracture, accompanied by the diminishing of surface wrinkles. Therefore, it opens a new avenue to enhance the stretchability of graphene materials, and enables the future applications of graphene and other 2D materials in stretchable and flexible electronics.

USNCTAM perspectives on mechanics in medicine

Over decades, the theoretical and applied mechanics community has developed sophisticated approaches for analysing the behaviour of complex engineering systems. Most of these approaches have targeted systems in the transportation, materials, defence and energy industries. Applying and further developing engineering approaches for understanding, predicting and modulating the response of complicated biomedical processes not only holds great promise in meeting societal needs, but also poses serious challenges. This report, prepared for the US National Committee on Theoretical and Applied Mechanics, aims to identify the most pressing challenges in biological sciences and medicine that can be tackled within the broad field of mechanics. This echoes and complements a number of national and international initiatives aiming at fostering interdisciplinary biomedical research. This report also comments on cultural/educational challenges. Specifically, this report focuses on three major thrusts in which we believe mechanics has and will continue to have a substantial impact. (i) Rationally engineering injectable nano/microdevices for imaging and therapy of disease. Within this context, we discuss nanoparticle carrier design, vascular transport and adhesion, endocytosis and tumour growth in response to therapy, as well as uncertainty quantification techniques to better connect models and experiments. (ii) Design of biomedical devices, including point-of-care diagnostic systems, model organ and multi-organ microdevices, and pulsatile ventricular assistant devices. (iii) Mechanics of cellular processes, including mechanosensing and mechanotransduction, improved characterization of cellular constitutive behaviour, and microfluidic systems for single-cell studies.

Direct observation of nonaffine tube deformation in strained polymer networks

We present a one-to-one comparison of polymer segmental fluctuations as measured by small angle neutron scattering in a network under deformation with those obtained by neutron spin echo spectroscopy. This allows an independent proof of the strain dependence of the chain entanglement length. The experimentally observed nonaffine square-root dependence of the tube channel on strain is in excellent agreement with theoretical predictions and permits us to exclude an often invoked nondeformed as well as affinely deformed tube.

COMPUTATIONAL SIMULATION OF ELASTOMER NANOCOMPOSITES: CURRENT PROGRESS AND FUTURE CHALLENGES

In the field of elastomer nanocomposites (ENCs), computational simulation technique is becoming more and more essential, as a result of its ability to provide important and clear information at the molecular level, which is always difficult to obtain or not accessible through experimental investigations. We focus on summarizing the progress achieved in the simulation research of three critical topics of ENCs, namely, (i) the dispersion mechanism (particularly polymer-mediated interparticle interaction, the “many-body” effect at high filler loading), (ii) the characterization of the nanoscale/microscale structure and dynamics [the modified chain configuration in the presence of nanoparticles (NPs), the interfacial binding strength determining the efficiency of the stress transfer, the possibly altered interfacial chain structure, interfacial segmental dynamics leading to the shift of the glass transition temperature Tg, the formation of the filler network and its structure, the chemical cross-linking process], and (iii) the macroscopic viscoelasticity (the Payne effect), mechanical reinforcement, and physical property (thermal conductivity). Since recently only limited simulation work has been carried out pertaining to ENCs, we discuss these three topics in light of the simulation and theoretical achievements of polymer nanocomposites (mainly polymer melts filled with NPs). Meanwhile, some relevant experimental studies are also included for better illustration. Furthermore, for each topic, three typically different reinforcing fillers, such as three-dimensional spherical, two-dimensional sheet, and one-dimensional rod NPs, separately corresponding to carbon black or silica, clay sheets, and carbon nanotubes intensively used in the practical applications of ENCs, are illustrated in order. In order to realize a comprehensive understanding of the structure–property relation and in the meantime to provide more practical guidelines for the engineering applications of ENCs, we investigate future simulation opportunities and difficulties.