Mechanical testing of glassy and rubbery polymers in numerical simulations: Role of boundary conditions in tensile stress experiments

2D MOFs and Zeolites for Composite Membrane and Gas Separation Applications: A Brief Review

Commercial membranes have predominantly been fabricated from polymers due to their economic viability and processability. This choice offers significant advantages in energy efficiency, cost-effectiveness, and operational simplicity compared to conventional separation techniques like distillation. However, polymeric membranes inherently exhibit a trade-off between their permeability and selectivity, which is summarized in the Robeson upper bound. To potentially surpass these limitations, mixed-matrix membranes (MMMs) can be an alternative solution, which can be constructed by combining polymers with inorganic additives such as metal-organic frameworks (MOFs) and zeolites. Incorporating high-aspect-ratio fillers like MOF nanosheets and zeolite nanosheets is of significant importance. This incorporation not only enhances the efficiency of separation processes but also reinforces the mechanical robustness of the membranes. We outline synthesis techniques for producing two-dimensional (2D) crystals (including nanocrystals with high aspect ratio) and provide examples of their integration into membranes to customize separation performances. Moreover, we propose a potential trajectory for research in the area of high-aspect-ratio materials-based MMMs, supported by a mathematical-model-based performance prediction.

Effect of Polymer Composition and Morphology on Mechanochemical Activation in Nanostructured Triblock Copolymers

The effect of composition and morphology on mechanochemical activation in nanostructured block copolymers was investigated in a series of poly(methyl methacrylate)-block-poly(n-butyl acrylate)-block-poly(methyl methacrylate) (PMMA-b-PnBA-b-PMMA) triblock copolymers containing a force-responsive spiropyran unit in the center of the rubbery PnBA midblock. Triblock copolymers with identical PnBA midblocks and varying lengths of PMMA end-blocks were synthesized from a spiropyran-containing macroinitiatior via atom transfer radical polymerization, yielding polymers with volume fractions of PMMA ranging from 0.21 to 0.50. Characterization by transmission electron microscopy revealed that the polymers self-assembled into spherical and cylindrical nanostructures. Simultaneous tensile tests and optical measurements revealed that mechanochemical activation is strongly correlated to the chemical composition and morphologies of the triblock copolymers. As the glassy (PMMA) block content is increased, the overall activation increases, and the onset of activation occurs at lower strain but higher stress, which agrees with predictions from our previous computational work. These results suggest that the self-assembly of nanostructured morphologies can play an important role in controlling mechanochemical activation in polymeric materials and provide insights into how polymer composition and morphology impact molecular-scale force distributions.

Computational Study of Mechanochemical Activation in Nanostructured Triblock Copolymers

Force-driven chemical reactions have emerged as an attractive platform for diverse applications in polymeric materials. However, the microscopic chain conformations and topologies necessary for efficiently transducing macroscopic forces to the molecular scale are not well-understood. In this work, we use a coarse-grained model to investigate the impact of network-like topologies on mechanochemical activation in self-assembled triblock copolymers. We find that mechanochemical activation during tensile deformation depends strongly on both the polymer composition and chain conformation in these materials. Activation primarily occurs in the tie chains connecting different glassy domains and in loop chains that are hooked onto each other by physical entanglements. Activation also requires a higher stress in materials having a higher glassy block content. Overall, the lamellar samples show the highest percent activation at high stress. In contrast, at low stress, the spherical morphology, which has the lowest glassy fraction, shows the highest activation. Additionally, we observe a spatial pattern of activation, which appears to be tied to distortion of the self-assembled morphology. Higher activation is observed in the tips of the chevrons formed during deformation of lamellar samples as well as in the centers between the cylinders in the cylindrical morphology. Our work shows that changes in the network-like topology in different morphologies significantly impact mechanochemical activation efficiencies in these materials, suggesting that this area will be a fruitful avenue for further experimental research.

Effects of Mechanical Stress on Insulation Structure and Performance of HV Cable

Mechanical stresses generated during manufacturing and laying process of high voltage cables can result in degradation of insulation properties, affecting the stable operation of the transmission system. Traditional test methods for testing the effect of mechanical stress on the insulation properties of polyethylene still have some shortcomings to be explored and it is able to explain the changes of the insulation properties of polyethylene under mechanical stress from a microscopic perspective. In order to further study the effect of stress on the insulation properties of polyethylene, microstructural changes, the breakdown field strength, conductivity and charge distribution of polyethylene at different elongation rates are investigated by a combination of experimental and molecular dynamics simulations. The results show that the increase in stress leads to a decrease in crystallinity and microcrystalline size of the material decrease. The untwisting and orientation of the polyethylene molecular chains during the stretching process can create cavities, resulting in an uneven sample distribution and thickness reduction, leading to a reduction in the breakdown field strength. Meanwhile, some crystal regions are transformed into amorphous regions. The loose amorphous regions facilitate the directional migration of carriers, resulting in the increase of conductivity. When the elongation ratio is smaller, the distance between the molecular chains increases and the trap depth of the specimen becomes shallower. This facilitates the migration of ions and electrons and increases the rate of decay of the surface potential. When the stretch is further increased, new traps are created by broken molecular chains to limit the movement of charges, decreasing the decay rate of the surface potential and reducing the insulation properties of the polyethylene. Meanwhile, the molecular dynamics model of semi-crystalline polyethylene was developed to observe the microstructure and energy changes during the stretching process. The conclusions in terms of tensile tests were verified from a microscopic perspective.

Coarse-grained molecular dynamics simulations of void generation and growth processes in the fracture of the lamellar structure of polyethylene

We investigate the void generation and growth process in semicrystalline polymers by large-scale coarse-grained molecular dynamics simulations. Voids are generated in the amorphous layers and grow spherically and then cylindrically, consistent with the results of previous experiments. Interestingly, the fusion of voids is indicated to enlarge the voids in the direction perpendicular to the stretching direction, but not beyond the crystalline layers. The orientational order along the stretching direction increased sharply before void generation, but the increase leveled off afterward. The simulations also clearly reveal that the crystalline layers break but do not bend in the fragmentation process. The dependence of the void growth process on stretching velocity and the stability levels of voids at constant strain are also discussed.

Crystallization of finite-extensible nonlinear elastic Lennard-Jones coarse-grained polymers

The ability of a simple coarse-grained finite-extensible nonlinear elastic (FENE) Lennard-Jones (LJ) polymer model to be crystallized is investigated by molecular dynamics simulations. The optimal FENE Lennard-Jones parameter combinations (ratio between FENE and LJ equilibrium distances) and the optimal lattice parameters are calculated for five different perfect crystallite structures: simple tetragonal, body-centered tetragonal, body-centered orthorhombic, hexagonal primitive, and hexagonal close packed. It was found that the most energetically favorable structure is the body-centered orthorhombic. Starting with an equilibrated polymer liquid and with the optimal parameters found for the body-centered orthorhombic, an isothermal treatment led to the formation of large lamellar crystallites with a typical chain topology: folded, loop, and tie chains, and with a crystallinity of about 60%-70%, similar to real semicrystalline polymers. This simple coarse-grained Lennard-Jones model provides a qualitative tool to study semicrystalline microstructures for polymers.

Analysis of local properties during a scratch test on a polymeric surface using molecular dynamics simulations

This work demonstrates a possible route to connect a particle (chain) based understanding with continuum mechanical questions about contact mechanics. The bond orientation, chain conformation and stress field of a polymer film were analyzed during scratch tests (tangential contact) using a molecular dynamics (MD) simulation approach. Scratch tests with a conical tip at constant scratching velocity were simulated on linear amorphous polymer surfaces at various temperatures and roughnesses of the tip and for various interactions between the tip and the particles of the polymer chains. The second Legendre polynomial (computed for small domains around the tip) gave the bond orientation inside the polymer film during sliding of the tip. The gyration tensor (layer-resolved in the direction of the polymer film thickness) provided information about the conformation of the polymer chains. These results allowed us to argue in favor of Briscoe's hypothesis (thin film sheared vs. "bulk" compressive behavior) concerning the friction properties of the polymer surfaces. Finally, the first stress measurements of the virial stress tensor (in sub-boxes placed in the MD cell) revealed a complex combination between compressive hydrostatic pressure and shear stress, which may be interpreted as a complex sheared domain at the interface.

Molecular dynamics simulation of rupture in glassy polymer bridges within filler aggregates

We present a series of nonequilibrium molecular dynamics simulations, investigating the rupture mechanisms in glassy polymer films confined between two solid surfaces. Such systems provide a useful model for the strong nonlinear reinforcement of rubber by colloidal filler particles. Depending on the degree of confinement three qualitatively different rupture modes have been found, which originate from the interplay of internal (polymer-polymer) and external (polymer-wall) interactions. In very thin films we observe the formation and stretching of many single-chain bridges between the confining surfaces. Progressing to thicker samples we observe fewer bridges, consisting of bundled polymer chains, and eventually just one large bridge in thick specimens. The yield stress and the elongational modulus of the polymer films have been calculated from the stress-strain curves at various temperatures and confinements and their behavior has been analyzed in terms of polymer-polymer and polymer-surface interaction energies. The thinnest films (5 monomer diameters) are always glassy in our simulations, while the others display a glass transition temperature around 0.50-0.55 (in units ε(0)/k(B) of the Lennard-Jones interaction energy), depending on their thickness. This range of values, which has been determined using both the nonequilibrium tensile simulations and equilibrium diffusion data, agrees with the transition temperature previously found by shear simulations [Baljon and Robbins, Science 271, 482 (1996)].

Mechanical behavior of linear amorphous polymers: comparison between molecular dynamics and finite-element simulations

This paper studies the rheology of weakly entangled polymer melts and films in the glassy domain and near the rubbery domain using two different methods: molecular dynamics (MD) and finite element (FE) simulations. In a first step, the uniaxial mechanical behavior of a bulk polymer sample is studied by means of particle-based MD simulations. The results are in good agreement with experimental data, and mechanical properties may be computed from the simulations. This uniaxial mechanical behavior is then implemented in FE simulations using an elasto-viscoelasto-viscoplastic constitutive law in a continuum mechanics (CM) approach. In a second step, the mechanical response of a polymer film during an indentation test is modeled with the MD method and with the FE simulations using the same constitutive law. Good agreement is found between the MD and CM results. This work provides evidence in favor of using MD simulations to investigate the local physics of contact mechanics, since the volume elements studied are representative and thus contain enough information about the microstructure of the polymer model, while surface phenomena (adhesion and surface tension) are naturally included in the MD approach.

Mechanical properties of the beetle elytron, a biological composite material

We determined the relationship between composition and mechanical properties of elytra (modified forewings that are composed primarily of highly sclerotized dorsal and less sclerotized ventral cuticles) from the beetles Tribolium castaneum (red flour beetle) and Tenebrio molitor (yellow mealworm). Elytra of both species have similar mechanical properties at comparable stages of maturation (tanning). Shortly after adult eclosion, the elytron of Tenebrio is ductile and soft with a Young's modulus (E) of 44 ± 8 MPa, but it becomes brittle and stiff with an E of 2400 ± 1100 MPa when fully tanned. With increasing tanning, dynamic elastic moduli (E') increase nearly 20-fold, whereas the frequency dependence of E' diminishes. These results support the hypothesis that cuticle tanning involves cross-linking of components, while drying to minimize plasticization has a lesser impact on cuticular stiffening and frequency dependence. Suppression of the tanning enzymes laccase-2 (TcLac2) or aspartate 1-decarboxylase (TcADC) in Tribolium altered mechanical characteristics consistent with hypotheses that (1) ADC suppression favors formation of melanic pigment with a decrease in protein cross-linking and (2) Lac2 suppression reduces both cuticular pigmentation and protein cross-linking.

Evolution of entanglements during the response to a uniaxial deformation of lamellar triblock copolymers and polymer glasses

Using coarse-grained molecular-dynamics simulations, a generic styrene-block-butadiene-block-styrene triblock copolymer under lamellar conformation is used in order to investigate the mutual entanglement evolution when a structure of alternating glassy (S)/rubbery (B) layers is submitted to an imposed deformation. By varying the amount of loop chains between each phase, i.e., noncrossing chains, it is possible to generate different types of S/B interface definitions. A specific boundary driven tensile strain protocol has been developed in order to mimic "real" experiments and measure the stress-strain curve. The same protocol is also applied to a reference state consisting in a directed glassy homopolymers, as well as to an isotropic glassy polymer. The evolution of initial mutual entanglements from the undeformed samples during the whole deformation process is monitored. It is shown for all considered systems that initial entanglements mostly participate to the preyield regime of the stress-strain curve and that this network is debonded during the strain-hardening regime. For triblocks with a non-null amount of crossing chains, the lower the amount is, the longer the memory effect of the initial entanglement network in the postyield regime is. On the fly distributions of entanglements, which depart from the postyield regime, depict memory effects and long-time correlations during the strain-hardening regime. For triblocks, loop chains reinforce these effects.

Fracture in glassy polymers: a molecular modeling perspective

Over the past 25 years, molecular modeling and simulations have provided important insights into the physics of deformation and fracture of glassy polymers. This review presents an overview of key results discussed in the context of experimentally observed polymer behavior. Both atomistic and coarse-grained polymer models have been used in different deformation protocols to study elastic properties, shear yielding, creep, physical aging, strain hardening and crazing. Simulations reproduce most of the macroscopic features of plasticity in polymer glasses such as stress-strain relations and creep response, and reveal information about the underlying atomistic processes. Trends of the shear yield stress with loading conditions, temperature and strain rate, and the atomistic dynamics under load have been systematically explored. Most polymers undergo physical aging, which leads to a history-dependent mechanical response. Simulations of strain hardening and crazing demonstrate the nature of polymer entanglements in the glassy state and the role of local plasticity and provide insight into the origin of fracture toughness of amorphous polymers.