Molecular simulations of gas transport in nitrile rubber and styrene butadiene rubber

Diffusion Behavior of VOC Molecules in Polyvinyl Chloride Investigated by Molecular Dynamics Simulation

Due to the threats posed by many volatile organic compounds (VOCs) to human health in indoor spaces via air, the mass transfer characteristics of VOCs are of critical importance to the study of their mechanism and control. As a significant part of the mass transfer process, diffusion widely exists in emissions from floors (e.g., PVC floors) and in sorption in porous materials. Molecular simulation studies by can provide unparalleled insights into the molecular mechanisms of VOCs. We construct the detailed atomistic structures of PVC blend membranes to investigate the diffusion behavior of VOC molecules (n-hexane) in PVC by molecular dynamics (MD). The variation in the diffusion coefficient of n-hexane in PVC with respect to temperature is in line with Arrhenius' law. The effect of temperature on the diffusion mechanism was investigated from the perspectives of free volume, cavity distribution and polymer chain mobility. It was found that the relationships between the diffusion coefficients of n-hexane in the polymer and the inverse fractional free volume are exponential and agree well with the free volume theory. Hopefully, this study will offer quantitative insights into the mass transport phenomena of VOCs within polymeric materials.

Molecular Dynamics Simulation of Oxidative Aging Effect on Diffusion Behaviors of Oxygen and Cyclohexane in NBR

The influences of thermal-oxidative aging on the diffusion behaviors of oxygen and cyclohexane in nitrile-butadiene rubber (NBR) at the micro-scale were investigated by molecular dynamics (MD) simulation. The two types of aged rubber models were established on the basis of rubber oxidative chains modified by the introduction of hydroxyl groups and carbonyl groups in rubber chains. The diffusion behaviors of oxygen and cyclohexane in NBR under different conditions were characterized by the fractional free volume (FFV), mean square displacement (MSD), diffusion coefficients, and diffusion trajectory. It turns out that the elevated temperature contributed to the increase in the free volume and diffusion range of oxygen and cyclohexane, while the compressive stress showed the reverse influence. Additionally, the introduction of oxidative polar functional groups (hydroxyl groups and carbonyl groups) in rubber chains lowered the flexibility of the rubber chains and promoted the formation of strong polar interaction, which further inhibits the diffusion of oxygen and cyclohexane.

M SB, George JJ, Rasulev B, Raghavan P. Experimental and Simulation Studies on Nonwoven Polypropylene-Nitrile Rubber Blend: Recycling of Medical Face Masks to an Engineering Product

The battle against the COVID-19 pandemic counters the waste management system, as billions of single-use face masks are used per day all over the world. Proper disposal of used face masks without jeopardizing the health and the environment is a challenge. Herein, a novel method for recycling of medical face masks has been studied. This method incorporates the nonwoven polypropylene (PP) fiber, which is taken off from the mask after disinfecting it, with acrylonitrile butadiene rubber (NBR) using maleic anhydride as the compatibilizer, which results in a PP-NBR blend with a high percentage economy. The PP-NBR blends show enhanced thermomechanical properties among which, 70 wt % PP content shows superior properties compared to other composites with 40, 50, and 60 wt % of PP. The fully Atomistic simulation of PP-NBR blend with compatibilizer shows an improved tensile and barrier properties, which is in good agreement with the experimental studies. The molecular dynamics simulation confirms that the compatibility between non-polar PP and polar NBR phases are vitally important for increasing the interfacial adhesion and impeding the phase separation.

Molecule Dynamics Simulation of the Effect of Oxidative Aging on Properties of Nitrile Rubber

The effects of oxidative aging on the static and dynamic properties of nitrile rubber at the molecular scale were investigated by molecular dynamics simulation. The aged nitrile rubber models were constructed by introducing hydroxyl groups and carbonyl groups into rubber molecular chains to mimic oxidative aging. The static and dynamic properties of the unaged and aged nitrile rubber under different conditions were evaluated by mean square displacement, self-diffusion coefficients, hydrogen bond, fractional free volume, radial distribution function, cohesive energy density and solubility parameter. The results show that the elevated temperature intensified significantly the mobility of rubber molecular chains and fractional free volume, while the compressive strain displayed the opposite effect resulting in packing and rearrangement of rubber chains. The introduction of hydroxyl groups and carbonyl groups enhanced the polarity, intermolecular interactions, the volume and rigidity of molecular chains, implying weaker mobility of molecular chains as compared to unaged models. The compressive strain and oxidative aging both decreased the fractional free volume, which inhibited gaseous and liquid diffusion into the rubber materials, and slowed down the oxidative aging rate. This study provides insights to better understand the effect of molecular changes due to oxidative aging on the structural and dynamic properties of rubber materials at the molecular level.

Exploring secondary interactions and the role of temperature in moisture-contaminated polymer networks through molecular simulations

Leveraging the state of absorbed moisture within a polymer network to identify physical and chemical features of the host material is predicated upon a clear understanding of the interaction between the polymer and a penetrant water molecule; an understanding that has remained elusive. Recent work has revealed that a novel damage detection method that exploits the very low baseline levels of water typically found in polymer matrix composites (PMC) may be a valuable tool in the composite NDE arsenal, provided that a clear understanding of polymer-water interaction can be obtained. Precise detection, location, and possible quantification of the extent of damage can be performed by characterizing the physical and chemical states of moisture present in an in-service PMC. Composite structures have a locally elevated dielectric constant near the damage sites due to a higher fraction of bulk ("free") water, which has a higher dielectric constant when compared to water molecules bound to the polymer network through secondary bonding interactions. In this study, we aim to get a clear atomistic scale picture of the interactions which drive the dielectric signature variations necessary for tracking damage. Molecular Dynamics (MD) simulations were used to explore the effect of temperature on the state of moisture in two epoxy matrices with identical chemical constituents but different morphologies. The motivation was to understand whether higher polarity binds a greater fraction of moisture even at higher temperatures, leading to suppressed dielectric activity. Consequently, the influence of secondary bonding interactions was investigated to understand the impact of temperature on the absorbed water molecules in a composite epoxy matrix.

Investigation of nanoparticle-polymer interaction in bio-based nanosilica-filled PLA/NR nanocomposites: molecular dynamics simulation

Molecular dynamics (MD) simulation, by employing the COMPASS force field, was utilized to investigate structural and thermal characteristics as well as interfacial interactions between components of nanocomposite consisting of poly(lactic acid) (PLA)/natural rubber (NR)/nanosilica, abbreviated as PSxN, where 1 ≤ x ≤ 7 and it represents the parts of SiO₂ nanoparticles added to the PLA/NR (PN) blend. Analysis of the obtained results including density (ρ), fractional free volume (FFV), glass transition temperature (T_g), interaction energy (E_interaction), and radial distribution function (RDF) of these nanocomposites was performed. Comparing E_interaction of nanocomposites with that of the PN blend showed that the interactions between the chains of the two polymers are highly dependent on the added amounts of silica nanoparticles, so that by adding silica to the PN blend to obtain PS1N and PS3N nanocomposites, the amount of E_interaction was reduced to a smaller positive value, which indicates the tendency of the nanocomposite's components to interact with each other. By further addition of silica nanoparticles to have PS5N and PS7N nanocomposites and then by analysis of the RDF results, it was found that the nanoparticles were not well dispersed in these two nanocomposites and they were accumulated in the NR rubbery phase. Therefore, the percolation threshold for silica loading on the PN blend is at most 3 parts (x = 3). These results as well as the other obtained simulation results were compared with the available experimental data, and the agreement observed between them approved the simulation procedure and validated the obtained results.

Molecular simulations of gas transport in hydrogenated nitrile butadiene rubber and ethylene-propylene-diene rubber

Diffusion and sorption of five gases (H₂, N₂, O₂, CO₂, CH₄) in hydrogenated nitrile butadiene rubber (HNBR) and ethylene–propylene–diene rubber (EPDM) have been investigated by molecular dynamics (MD) and grand canonical Monte Carlo (GCMC) simulations. The diffusion coefficients of gas molecules in HNBR and EPDM are well correlated with the effective penetrant diameter except for CO₂. CO₂ shows a lower diffusion coefficient due to its linear shape. Additionally, the favorable interaction between CO₂ and HNBR is another factor for its lower diffusion coefficient in HNBR. HNBR shows lower diffusion coefficients than EPDM. This is because the polar –CN groups in HNBR chains increase interchain cohesion and result in tight intermolecular packing, low free volume and poor chain mobility, which decreases the diffusion coefficients of HNBR. The solubility coefficients of CH₄, O₂, N₂ and H₂ in HNBR are lower than those in EPDM, which is a result of the weak HNBR–penetrant interactions and low free volume of HNBR. However, the solubility coefficient of CO₂ in HNBR is higher than in EPDM. This is attributed to the strong interaction between CO₂ and HNBR. H₂, O₂, N₂ and CH₄ show lower permeability coefficients in HNBR than in EPDM, while CO₂ has higher permeability coefficients in HNBR. These molecular details provide critical information for the understanding of structures and gas transport between HNBR and EPDM.

Diffusion and sorption of five gases (H₂, N₂, O₂, CO₂, CH₄) in HNBR and EPDM were explored by MD and GCMC simulations.

Atomistic QM/MM simulations of the strength of covalent interfaces in carbon nanotube–polymer composites

We use a QM/MM approach to simulate CNT pull-out from crosslinked polyethylene studying how interfacial strength depends on its chemistry.

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.

Molecular Simulation of Gas Solubility in Nitrile Butadiene Rubber

Molecular simulation is used to compute the solubility of small gases in nitrile butadiene rubber (NBR) with a Widom particle-insertion technique biased by local free volume. The convergence of the method is examined as a function of the number of snapshots upon which the insertions are performed and the number of insertions per snapshot and is compared to the convergence of the unbiased Widom insertion technique. The effect of varying the definition of local free volume is also investigated. The acrylonitrile content of the polymer is altered to examine its influence on the solubility of helium, CO₂, and H₂O, and the solubilities of polar gases are found to be enhanced relative to those of nonpolar gases, in qualitative agreement with experiment. To probe this phenomenon further, the solubilities are decomposed into contributions from the neighborhoods of different atoms, using a Voronoi cell construction, and a strong bias is found for CO₂ and H₂O in particular to be situated near nitrogen sites in the elastomer. Temperature is shown to suppress the solubility of CO₂ and H₂O but to increase that of helium. Increasing pressure is found to suppress the solubility of all gases but at different rates, according to a balance between their molecular sizes and electrostatic interactions with the polymer. These results are relevant to the use of NBR seals at elevated temperatures and pressures, such as in oil and gas wells.

A study on the relationship between polycarbonate microstructure and performance as determined by a combined experimental and molecular dynamics simulation method

The influence of the microstructure of polycarbonate (PC) on performance was systematically investigated by both experimental method and molecular simulation. Yield stress, impact strength, molecular weight, and transmittance were used to distinguish the degradation processes between different PCs, and thermal degradation kinetics was studied to obtain the activation energy. At the molecular level, through 13C nuclear magnetic resonance (NMR) spectroscopy, it was observed that PCs have a more polar group of benzene rings, resulting in the high density, dielectric constant, and tensile modulus. Meanwhile, molecular dynamics (MD) simulation was employed under a polymer consistent force field force field. Specific volume and mechanical property were analyzed to investigate the thermodynamic property. The molecular dynamics simulation and experimental results on half decomposition temperature (T1/2), refraction index, flow activation energy, average density, cohesive energy density, glass transition temperature (Tg), and elastic modulus had good agreement. Therefore, it was indicated that the molecular simulation could successfully study the characteristics and properties. The fundamental studies would be expected to supply useful information for designing materials and optimizing processing technology.

Pure and Modified Co-Poly(amide-12-b-ethylene oxide) Membranes for Gas Separation Studied by Molecular Investigations

This paper deals with a theoretical investigation of gas transport properties in a pure and modified PEBAX block copolymer membrane with N-ethyl-o/p-toluene sulfonamide (KET) as additive molecules. Molecular dynamics simulations using COMPASS force field, Gusev-Suter Transition State Theory (TST) and Monte Carlo methods were used. Bulk models of PEBAX and PEBAX/KET in different copolymer/additive compositions were assembled and analyzed to evaluate gas permeability and morphology to characterize structure-performance relationships.