COMPUTATIONAL SIMULATION OF ELASTOMER NANOCOMPOSITES: CURRENT PROGRESS AND FUTURE CHALLENGES

Molecular dynamics simulation of the impact of the surface topology of carbon black on the mechanical properties of elastomer nanocomposites

Carbon black has always played a pivotal role in reinforcing elastomers because it remarkably improves the mechanical properties. The reinforcing effect of carbon black is influenced by its grades, which mainly depend on the difference in the structure of the carbon black particles. Despite many traditional experiments on the performance of carbon black composites, there has been less emphasis on reinforcement mechanisms due to the challenges associated with unraveling the intermolecular interactions. In this paper, a coarse grained molecular dynamics simulation was employed to examine the relationship between the morphology of the carbon black particles and the mechanical properties of the elastomer nanocomposites. Specifically, three different morphological carbon black nanoparticle models, including the smooth particle model, rough particle model, and the rough ellipsoid model, were constructed first. We then focused on investigating the changes of the mechanical properties by systematically varying the filling fraction of the carbon black particles, and the strength of the interfacial interaction between the filler and the rubber. The results indicated that the surface roughness and the filler's shape had a significant impact on the mechanical properties of the filled rubber models. The mechanical enhancement effect of the rough ellipsoidal carbon black is around 50-400% higher than that of the smooth carbon black, and the stronger the interfacial interactions, the more pronounced the enhancement. In addition, the rough ellipsoid filled system has low hysteresis, low permanent deformation, and high fatigue resistance. In general, this work explores the strengthening mechanism of carbon black on the elastomer at the molecular level and generates new insight into the design and fabrication of novel reinforcing fillers.

Fracture in Silica/Butadiene Rubber: A Molecular Dynamics View of Design-Property Relationships

Despite intense investigation, the mechanisms governing the mechanical reinforcement of polymers by dispersed nanoparticles have only been partially clarified. This is especially true for the ultimate properties of the nanocomposites, which depend on their resistance to fracture at large deformations. In this work, we adopt molecular dynamics simulations to investigate the mechanical properties of silica/polybutadiene rubber, using a quasi-atomistic model that allows a meaningful description of bond breaking and fracture over relatively large length scales. The behavior of large nanocomposite models is explored systematically by tuning the cross-linking, grafting densities, and nanoparticle concentration. The simulated stress-strain curves are interpreted by monitoring the breaking of chemical bonds and the formation of voids, up to complete rupture of the systems. We find that some chemical bonds, and particularly the S-S linkages at the rubber-nanoparticle interface, start breaking well before the appearance of macroscopic features of fracture and yield.

Molecular Dynamics Study on Mechanical Properties of Interface between Urea-Formaldehyde Resin and Calcium-Silicate-Hydrates

Microcapsule based self-healing concrete can automatically repair damage and improve the durability of concrete structures, the performance of which depends on the binding behavior between the microcapsule wall and cement matrix. However, conventional experimental methods could not provide detailed information on a microscopic level. In this paper, through molecular dynamics simulation, three composite models of Tobermorite (Tobermorite 9 Å, Tobermorite 11 Å, Tobermorite 14 Å), a mineral similar to Calcium-Silicate-Hydrate (C-S-H) gel, with the linear urea-formaldehyde (UF), the shell of the microcapsule, were established to investigate the mechanical properties and interface binding behaviour of the Tobermorite/UF composite. The results showed that the Young's modulus, shear modulus and bulk modulus of Tobermorite/UF were lower than that of 'pure' Tobermorite, whereas the tensile strength and failure strain of Tobermorite/UF were higher than that of 'pure' Tobermorite. Moreover, through radial distribution function (RDF) analysis, the connection between Tobermorite and UF found a strong interaction between Ca, N, and O, whereas Si from Tobermorite and N from UF did not contribute to the interface binding strength. Finally, high binding energy between the Tobermorite and UF was observed. The research results should provide insights into the interface behavior between the microcapsule wall and the cement matrix.

Tailoring the mechanical properties by molecular integration of flexible and stiff polymer networks

Designing a multiple-network structure at the molecular level to tailor the mechanical properties of polymeric materials is of great scientific and technological importance. Through the coarse-grained molecular dynamics simulation, we successfully construct an interpenetrating polymer network (IPN) composed of a flexible polymer network and a stiff polymer network. First, we find that there is an optimal chain stiffness for a single network (SN) to achieve the best stress-strain behavior. Then we turn to study the mechanical behaviors of IPNs. The result shows that the stress-strain behaviors of the IPNs appreciably exceed the sum of that of the corresponding single flexible and stiff network, which highlights the advantage of the IPN structure. By systematically varying the stiffness of the stiff polymer network of the IPNs, optimal stiffness also exists to achieve the best performance. We attribute this to a much larger contribution of the non-bonded interaction energy. Last, the effect of the component concentration ratio is probed. With the increase of the concentration of the flexible network, the stress-strain behavior of the IPNs is gradually enhanced, while an optimized concentration (around 60% molar ration) of the stiff network occurs, which could result from the dominant role of the enthalpy rather than the entropy. In general, our work is expected to provide some guidelines to better tailor the mechanical properties of the IPNs made of a flexible network and a stiff network, by manipulating the stiffness of the stiff polymer network and the component concentration ratio.

Mechanisms of reinforcement in polymer nanocomposites

Qualitatively different stress–strain responses of polymer nanocomposites are shown to result from the dynamical evolution of three principal molecular structural motifs in the polymer–filler network.

Self-Assembly of Block Copolymer Chains To Promote the Dispersion of Nanoparticles in Polymer Nanocomposites

In this paper we adopt molecular dynamics simulations to study the amphiphilic AB block copolymer (BCP) mediated nanoparticle (NP) dispersion in polymer nanocomposites (PNCs), with the A-block being compatible with the NPs and the B-block being miscible with the polymer matrix. The effects of the number and components of BCP, as well as the interaction strength between A-block and NPs on the spatial organization of NPs, are explored. We find that the increase of the fraction of the A-block brings different dispersion effect to NPs than that of B-block. We also find that the best dispersion state of the NPs occurs in the case of a moderate interaction strength between the A-block and the NPs. Meanwhile, the stress-strain behavior is probed. Our simulation results verify that adopting BCP is an effective way to adjust the dispersion of NPs in the polymer matrix, further to manipulate the mechanical properties.

Reconsideration of the Rheology of Silica Filled Natural Rubber Compounds

There is substantial progress along with giant debate in reinforcement mechanisms in relation to structured filler network and heterogeneously retarded polymer dynamics, while the dissipation behaviors have never been clarified for nanoparticle filled polymers. Herein dynamic rheological behaviors of silica filled natural rubber were investigated. Master curves of linear rheology in the hydrodynamic regime and those of the nonlinear Payne effect at a predetermined frequency were created, disclosing a leading role of dynamically retarded bulk rubbery phase to the hydrodynamic regime and a leading role of molecular disentanglement in the bulk phase to the Payne effect. The methodology is able to account for both reinforcement and dissipation of the compounds as a function of filler content. Furthermore, a frequency-dependent hydrodynamic to non-hydrodynamic transition is revealed, revealing the importance of the relaxation of chains in the bulk phase to both reinforcement and dissipation of the compounds. It is suggested that the dynamics of the bulk phase play a critical role for the rheology in the hydrodynamic regime while the fractal filler aggregates become dominative only in the terminal non-hydrodynamic regime where the bulk phase relaxes sufficiently.

Designing polymer nanocomposites with a semi-interpenetrating or interpenetrating network structure: toward enhanced mechanical properties

Semi-interpenetrating and interpenetrating network structures for the uniform dispersion of NPs and the reinforced mechanical properties of polymer nanocomposites.

Molecular insights into the effect of graphene packing on mechanical behaviors of graphene reinforced cis-1,4-polybutadiene polymer nanocomposites

We adopt molecular dynamics simulation to study the graphene packing patterns on chain structure, dynamics, uniaxial tension and visco-elastic behaviors.

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.

Filler reinforcement in cross-linked elastomer nanocomposites: insights from fully atomistic molecular dynamics simulation

Using a fully atomistic model, we perform large-scale molecular dynamics simulations of sulfur-cured polybutadiene (PB) and nanosilica-filled PB composites. A well-integrated network without sol fraction is built dynamically by cross-linking the coarse-grained precursor chains in the presence of embedded silica nanoparticles. Initial configurations for subsequent atomistic simulations are obtained by reverse mapping of the well-equilibrated coarse-grained systems. Based on the concept of "maximally inflated knot" introduced by Grosberg et al., we show that the networks simulated in this study behave as mechanically isotropic systems. Analysis of the network topology in terms of graph theory reveals that mechanically inactive tree-like structures are the dominant structural components of the weakly cross-linked elastomer, while cycles are mainly responsible for the transmission of mechanical forces through the network. We demonstrate that quantities such as the system density, thermal expansion coefficient, glass transition temperature and initial Young's modulus can be predicted in qualitative and sometimes even in quantitative agreement with experiments. The nano-filled system demonstrates a notable increase in the glass transition temperature and an approximately two-fold increase in the nearly equilibrium value of elastic modulus relative to the unfilled elastomer even at relatively small amounts of filler particles. We also examine the structural rearrangement of the nanocomposite subjected to tensile deformation. Under high strain-rate loading, the formation of structural defects (microcavities) within the polymer bulk is observed. The nucleation and growth of cavities in the post-yielding strain hardening regime mainly take place at the elastomer/nanoparticle interfaces. As a result, the cavities are concentrated just near the embedded nanoparticles. Therefore, while the silica nanofiller increases the elastic modulus of the elastomer, it also creates a more defective structure of higher energy in comparison with the unfilled network.

Effect of chemical structure of elastomer on filler dispersion and interactions in silica/solution-polymerized styrene butadiene rubber composites through molecular dynamics simulation

The dynamic properties, filler–rubber interactions, and filler dispersion in silica/SSBR composites with various chemical structures of SSBR were studied using MD. Competing effects led to the existence of an optimum modifier content of 14.2 wt%.

Impact of Interaction Strength and Surface Heterogeneity on the Dynamics of Adsorbed Polymers

We present molecular dynamics simulations of bead-and-spring polymer chains on chemically heterogeneous, energetically disordered surfaces at near-monolayer coverages. The surfaces consist of random mixtures of weakly (W) and strongly (S) attractive sites. We explore systematically the effect of surface composition on the diffusive dynamics of the chains. The polymer diffusion coefficients have a near-Arrhenius temperature dependence, with activation energies which have a nonmonotonic dependence on the fraction of S sites. In other words, we see a nonmonotonic dependence of the interfacial polymer dynamics on its affinity with the surface, when the latter involves some heterogeneity. The maximum activation energy belongs to the surface containing 75% S and 25% W sites, which combines near-maximum average polymer-surface interactions with near-maximum spread or disorder in these interactions. Our results have interesting implications for polymer adhesion and friction and structure-property relationships in polymer nanocomposites.

Elucidating and tuning the strain-induced non-linear behavior of polymer nanocomposites: a detailed molecular dynamics simulation study

By setting up a coarse-grained model of polymer nanocomposites, we monitored the change in the elastic modulus as a function of the strain, derived from the stress-strain behavior by determining uniaxial tension and simple shear of two typical spatial distribution states (aggregation and dispersion) of nanoparticles (NPs). In both these cases, we observed that the elastic modulus decreases non-linearly with the increase of strain and reaches a low plateau at larger strains. This phenomenon is similar to the so-called "Payne effect" for elastomer nanocomposites. Particularly, the modulus of the aggregation case is more sensitive to the imposed strain. By examining the structural parameters, such as the number of neighboring NPs, coordination number of NPs, root-mean-squared average force exerted on the NPs, local strain, chain conformations (bridge, dangle, loop, interface bead and connection bead), and the total interaction energy of NP-polymer and NP-NP, we inferred that the underlying mechanism of the aggregation case is the disintegration of the NP network or clusters formed through direct contact; however, for the dispersion case, the non-linear behavior is attributed to the destruction of the NP network or clusters formed through the bridging of adsorbed polymer segments among the NPs. The former physical network is influenced by NP-NP interaction and NP volume fraction, while the latter is influenced by NP-polymer interaction and NP volume fraction. Lastly, we found that for the dispersion case, further increasing the inter-particle distance or grafting NPs with polymer chains can effectively reduce the non-linear behavior due to the decrease of the physical network density. In general, this simulation work, for the first time, establishes the correlation between the micro-structural evolution and the strain-induced non-linear behavior of polymer nanocomposites, and sheds some light on how to reduce the "Payne effect".