Molecular simulation of the effect of temperature and architecture on polyethylene barrier properties

Hop, Skip, and Jump: Hydrogen Molecular Transport through Amorphous Polyethylene Matrices Studied via Molecular Dynamics Simulations

In the pursuit of advancing and diversifying energy technologies for a more sustainable future, the possibilities of hydrogen (H2) usage will broaden, as will our understanding of its containment materials. Polyethylene (PE) has a vast assortment of uses and applications, which are growing with demands for alternative energy possibilities. One use of PE liner is as a prime candidate for nonmetallic piping and pressurized type IV storage devices. Such applications require PE to effectively prevent H2 transport through containment systems. To study the molecular transport mechanism of hydrogen through polymeric barriers, a system containing hydrogen molecules absorbed within amorphous PE is modeled here using molecular dynamics simulations. The simulations are conducted within a range of temperatures that span the glass transition temperature of amorphous PE. The simulated PE displays bulk density, radius of gyration, and self-diffusion coefficient that are consistent with experimental data. The simulated trajectories are interrogated to study the movement of the guest gas molecules. The results show that the diffusion coefficients increase with temperature, as expected. However, the mobility of the PE chains is found to affect the mobility of absorbed H2 molecules to a much lower extent than it affects that of CH4 molecules because of the much smaller size of the former than of the latter guest. From a molecular perspective, a "hopping" mechanism is observed, according to which H2 molecules hop between one vacant free volume space to another within the polymer matrix, in combination with longer, straight, undisturbed "jumps" or "skips" along directions aligned with regions of ordered PE chains. This suggests that the orientation of the crystallites within the semicrystalline PE matrix affects the H2 containment. Implications of these findings toward PE usage as containment material are discussed.

Sorption and Diffusion of Methane and Carbon Dioxide in Amorphous Poly(alkyl acrylates): A Molecular Simulation Study

Molecular simulations were carried out to understand the structural features and the sorption and diffusion behavior of methane and carbon dioxide in amorphous poly(alkyl acrylates) in the temperature range of 300-600 K. The hybrid Monte Carlo/molecular dynamics approach was employed to address the effects of polymer swelling and framework flexibility on the gas sorption. Simulations show that the glass-transition temperature decreases with the side-chain length of poly(alkyl acrylate), consistent with experiments. This is due to the fact that the shielding of the polar ester groups increases with the side-chain length. The simulated sorption isotherms for methane and carbon dioxide were in agreement with the experimental data. The polymer swelling becomes more pronounced, especially in the case of sorption of carbon dioxide. A significant swelling occurs, possibly because of the greater interaction between carbon dioxide and the polar ester groups in the polymers. The uptake of methane and carbon dioxide by poly(alkyl acrylates) generally increases with the side-chain length. Our simulations confirm the experimental findings that the diffusion coefficients of methane and carbon dioxide in poly(alkyl acrylates) increase with the side-chain length. Interestingly, the activation energies of gas diffusion decrease with the side-chain length. The diffusion coefficients of the penetrants have an exponential relationship with the void fraction, which is in agreement with the free volume theory.

Transport Diffusion of Light Gases in Polyethylene Using Atomistic Simulations

We explore the temperature dependence of the self-, corrected-, and transport-diffusivities of CO₂, CH₄, and N₂ in a polyethylene (PE) polymer membrane through equilibrium molecular dynamics simulations. We also investigate the morphology of the polymer membrane based on the intermolecular radial distribution function, free volume, and pore size distribution analysis. The results indicate the existence of 1.5-3 Å diameter pores in the PE membrane, and with the increase in the temperature, the polymer swells linearly with changing slope at 450 K in the absence of gas and exponentially in the presence of gas. The gas adsorption isotherms extracted via a two-step methodology, considering the dynamics and structural transitions in the polymer matrix upon gas adsorption, were fitted using a "two-mode sorption" model. Our results suggest that CO₂ adsorbs strongly, whereas N₂ shows weak adsorption in PE. The results demonstrate that CO₂ is more soluble, whereas N₂ is least soluble. Further, it is found that an increase in the temperature negatively impacts the solubility of CO₂ and CH₄ but positively for N₂; this reverse solubility behavior is due to increased availability of pores accessible to N₂, which are kinetically closed at the lowest temperatures. The reported self-diffusivities of the gases from our simulations are on the order of 10^-6 cm²/s, consistent with the experimental evidence, whereas transport-diffusivities are 2 orders of magnitude higher than self-diffusivities. Furthermore, the temperature dependence of the self-diffusivity follows Arrhenius behavior, whereas the transport-diffusivity follows non-Arrhenius behavior having different activation energies in low and high temperature regions. Also, it is seen that loading has little effect on the self- and corrected-diffusion coefficients of all gases in the PE membrane.

Spontaneous crystallization in athermal polymer packings

We review recent results from extensive simulations of the crystallization of athermal polymer packings. It is shown that above a certain packing density, and for sufficiently long simulations, all random assemblies of freely-jointed chains of tangent hard spheres of uniform size show a spontaneous transition into a crystalline phase. These polymer crystals adopt predominantly random hexagonal close packed morphologies. An analysis of the local environment around monomers based on the shape and size of the Voronoi polyhedra clearly shows that Voronoi cells become more spherical and more symmetric as the system transits to the ordered state. The change in the local environment leads to an increase in the monomer translational contribution to the entropy of the system, which acts as the driving force for the phase transition. A comparison of the crystallization of hard-sphere polymers and monomers highlights similarities and differences resulting from the constraints imposed by chain connectivity.

How to predict diffusion of medium-sized molecules in polymer matrices. From atomistic to coarse grain simulations

The normal diffusion regime of many small and medium-sized molecules occurs on a time scale that is too long to be studied by atomistic simulations. Coarse-grained (CG) molecular simulations allow to investigate length and time scales that are orders of magnitude larger compared to classical molecular dynamics simulations, hence providing a valuable approach to span time and length scales where normal diffusion occurs. Here we develop a novel multi-scale method for the prediction of diffusivity in polymer matrices which combines classical and CG molecular simulations. We applied an atomistic-based method in order to parameterize the CG MARTINI force field, providing an extension for the study of diffusion behavior of penetrant molecules in polymer matrices. As a case study, we found the parameters for benzene (as medium sized penetrant molecule whose diffusivity cannot be determined through atomistic models) and Poly (vinyl alcohol) (PVA) as polymer matrix. We validated our extended MARTINI force field determining the self diffusion coefficient of benzene (2.27·10⁻⁹m² s⁻¹) and the diffusion coefficient of benzene in PVA (0.263·10⁻¹² m² s⁻¹). The obtained diffusion coefficients are in remarkable agreement with experimental data (2.20·10⁻⁹m² s⁻¹ and 0.25·10⁻¹² m² s⁻¹, respectively). We believe that this method can extend the application range of computational modeling, providing modeling tools to study the diffusion of larger molecules and complex polymeric materials.

Combined molecular algorithms for the generation, equilibration and topological analysis of entangled polymers: methodology and performance

We review the methodology, algorithmic implementation and performance characteristics of a hierarchical modeling scheme for the generation, equilibration and topological analysis of polymer systems at various levels of molecular description: from atomistic polyethylene samples to random packings of freely-jointed chains of tangent hard spheres of uniform size. Our analysis focuses on hitherto less discussed algorithmic details of the implementation of both, the Monte Carlo (MC) procedure for the system generation and equilibration, and a postprocessing step, where we identify the underlying topological structure of the simulated systems in the form of primitive paths. In order to demonstrate our arguments, we study how molecular length and packing density (volume fraction) affect the performance of the MC scheme built around chain-connectivity altering moves. In parallel, we quantify the effect of finite system size, of polydispersity, and of the definition of the number of entanglements (and related entanglement molecular weight) on the results about the primitive path network. Along these lines we approve main concepts which had been previously proposed in the literature.