Molecular dynamics simulation of diffusion and permeation of gases in polystyrene

Molecular Simulation Study on the Hydrogen Permeation Behavior and Mechanism of Common Polymers

This research aimed to provide an understanding of the selection and safe application of pipeline liner materials for hydrogen transport by examining the permeation properties and mechanisms of hydrogen within polymers commonly used for this purpose, such as high-density polyethylene (HDPE) and ethylene-vinyl alcohol copolymer (EVOH), through molecular simulation. The study was carried out within defined operational parameters of temperature (ranging from room temperature to 80 °C) and pressure (from 2.5 to 10 MPa) that are pertinent to hydrogen pipeline infrastructures. The results reveal that with an increase in temperature from 30 °C to 80 °C, the solubility, diffusion, and permeability coefficients of hydrogen in HDPE increase by 18.7%, 92.9%, and 129.0%, respectively. Similarly, in EVOH, these coefficients experience increments of 15.9%, 81.6%, and 112.7%. Conversely, pressure variations have a negligible effect on permeability in both polymers. HDPE exhibits significantly higher hydrogen permeability compared to EVOH. The unique chain segment configuration of EVOH leads to the formation of robust hydrogen bonds among the hydroxyl groups, thereby impeding the permeation of hydrogen. The process by which hydrogen is adsorbed in polymers involves aggregation at low potential energy levels. During diffusion, the hydrogen molecule primarily vibrates within a limited range, with intermittent occurrences of significant hole-to-hole transitions over larger distances. Hydrogen exhibits a stronger interaction with HDPE compared to EVOH, leading to a higher number of adsorption sites and increased hydrogen adsorption capacity in HDPE. Hydrogen molecules move more actively in HDPE than in EVOH, exhibiting greater hole amplitude and more holes in transition during the diffusion process.

Deep eutectic solvent-based blended membranes for ultra-super selective separation of SO2

SO2 is a major atmospheric pollutant leading to acid rain and smog. As a new generation of green solvents, deep eutectic solvents (DESs) have been widely investigated for gas capture. Nevertheless, studies on DES-based membranes for SO2 separation are yet minimal. Herein, we devised polymer/DES blended membranes comprising 1-butyl-3-methyl-imidazolium bromide ([Bmim]Br)/diethylene glycol (DEG) DES and poly (vinylidene fluoride) (PVDF), and these membranes were firstly used for selective separation of SO2 from N2 and CO2. The permeability of SO2 reaches up to 17480 Barrer (0.20 bar, 40 ºC) in PVDF/DES blended membrane containing 50 wt% of [Bmim]Br/DEG (2:1), with ultrahigh SO2/N2 and SO2/CO2 selectivity of 3690 and 211 obtained, respectively, far exceeding those in the state-of-the-art membranes reported in literature. The highly-reversible multi-site interaction between SO2 and [Bmim]Br/DEG DES was revealed by spectroscopic analysis. Furthermore, the PVDF/DES blended membrane was also able to efficiently and stably separate SO2/CO2/N2 (2.5/15/82.5%) mixed gas for at least 100 h. This work demonstrates for the first time that [Bmim]Br-based DESs are very efficient media for membrane separation of SO2. The easy preparation, low cost and high performance enable polymer/DES blended membranes to be promising candidates for flue gas desulfurization.

A theoretical probe into the separation of CO2/CH4/N2 mixtures with polysulfone/polydimethylsiloxane-nano zinc oxide MMM

In the current investigation, molecular dynamics (MD) and Grand Canonical Monte Carlo (GCMC) simulation as remarkable and competent approaches have been employed for understanding structural and transport properties of MMMs in the realm of gas separation. The two commonly used polymers i.e. polysulfone (Psf) and polydimethylsiloxane (PDMS) as well as zinc oxide (ZnO) nanoparticle (NP) were used to carefully examine the transport properties of three light gasses (CO₂, N₂ and CH₄) through simple Psf, Psf/PDMS composite loaded by different amounts of ZnO NP. Also, the fractional free volume (FFV), X-ray diffraction (XRD), glass transition temperature (T_g), and Equilibrium density were calculated to scrutinize the structural characterizations of the membranes. Moreover, the effect of feed pressure (4-16 bar) on gas separation performance of simulated MMMs was investigated. Results obtained in different experiments showed a clear improvement in the performance of simulated membranes by adding PDMS to PSf matrix. The selectivity of studied MMMs was in the range from 50.91 to 63.05 at pressures varying from 4 to 16 bar for the CO₂/N₂ gas pair, whereas the corresponding value for CO₂/CH₄ system was found to be in the range 27.27-46.24. For 6 wt% ZnO in 80%PSf + 20%PDMS membrane, high permeabilities of 78.02, 2.86 and 1.33 barrers were observed for CO₂, CH₄ and N₂ gases, respectively. The 90%PSf + 10%PDMS membrane with 2% ZnO had a highest CO₂/N₂ selectivity value of 63.05 and its CO₂ permeability at 8 bar was 57 barrer.

Molecular Structure and Dynamics in Wet Gecko β-Keratin

Molecular dynamics simulations are performed to investigate the molecular picture of water sorption in gecko keratin and the influence of relative humidity (RH) on the local structure and dynamics in water-swollen keratin. At low RHs, water sorption occurs through hydrogen bonding of water with the hydrophilic groups of keratin. At high RHs (>80%), additional water molecules connect to the first "layer" of amide-connected water molecules (multimolecular sorption) through hydrogen bonds, giving rise to a sigmoidal shape of the sorption isotherm. This causes the formation of large chain-like clusters surrounding the hydrophilic groups of keratin, which upon a further increase of the RH form a percolating water network. An examination of the dynamics of water molecules sorbed in keratin demonstrates that there are two states, bound and free, for water. The dynamics of water in these states depends on the RH. At low RHs, large-scale translational motions of tightly bound water molecules to keratin are needed to remake the entire hydration shell of the keratin. At high RHs (>80%), the water molecules more quickly exchange between the two states. The center-of-mass mean-square displacement of water molecules indicates a hopping motion of water molecules in the keratin solvation shell. The hopping mechanism is more pronounced at RHs < 80%. At higher RHs, water translation through water clusters (water network) dominates. We have observed two regimes for the dependence of dynamical properties on the RH: a regime of gradual increase of the dynamics over 10% < RH < 80% and a regime of drastic dynamic acceleration at RH > 80%. The latter regime begins exactly where the water uptake and the volume swelling also increase much more and where a drastic change in the elastic properties of gecko keratin has been observed. A nearly linear relation between the relaxation times for all dynamical processes and the water content of gecko keratin is observed.

Diffusion Mechanisms of Dissolved Gases in Transformer Oil Influenced with Moisture Based on Molecular Dynamics Simulation

Dissolved Gas Analysis (DGA) of insulating oil is widely used for diagnosing transformer incipient faults. Moisture is a major contaminant and degradation byproduct of transformer insulating oil. In this paper, molecular dynamics simulation was used to study the influence of moisture on the diffusion movement of dissolved gases in the insulating oil. Cycloalkanes (C₂₀H₄₂), alkanes (C₂₀H₃₈), and aromatic hydrocarbons (C₂₀H₂₆) are selected as the basic structural units in the constructed transformer oil simulation system. 0%, 1%, 3%, and 5% moisture molecules are added to insulating oil, respectively, and the insulating oil generates seven kinds of gas molecules through cracking. With an anhydrous model used as a benchmark, we researched the diffusion trajectory, the diffusion coefficient (D), free volume (V _F), and the moisture-gas interaction energy of each gas molecule as a function of moisture content. Through this study, we found that the increase of moisture content enlarges the V _F value of dissolved gas in insulating oil, which makes the gas in oil easier to diffuse. Besides, the moisture can also alter the mean square displacement (MSD) of dissolved gases; the mutual energy of molecules is mainly affected by the electrostatic interaction energy. This study can contribute to a better understanding of the influence of different moisture contents on the diffusion movement of dissolved gas in transformer oil from the micro level.

Permeation Characteristics of CH4 in PVDF with Crude Oil-Containing

The liner of reinforced thermoplastic composite pipes (RTPs) used for oil and gas gathering and transportation experienced blister failure due to gas permeation. Few reports have appeared on the problem of gas permeation in thermoplastics with absorbed crude oil. Accordingly, the permeability of CH₄ in polyvinylidene fluoride (PVDF) containing crude oil was studied at the normal service conditions by molecular simulations. The results showed that the solubility coefficients of CH₄ in PVDF containing crude oil were much lower than those in pure PVDF. It can be concluded that the crude oil molecules absorbed into PVDF occupied certain adsorption sites, resulting in a decrease in the adsorption capacity of CH₄ molecules in PVDF. The diffusion coefficients of CH₄ in oil-containing PVDF were significantly greater than in PVDF. This is because the absorption of oil molecules leads to the volume swelling of PVDF and then increases the free volume for diffusion. The permeation process showed that CH₄ molecules were selective-aggregate adsorbed in the region with low potential energy in oil-containing PVDF firstly, and then they vibrated within the holes of PVDF containing oil in most cases and jumped into the neighboring holes at high temperatures and pressures.

Molecular Simulation on Permeation Behavior of CH4/CO2/H2S Mixture Gas in PVDF at Service Conditions

Reinforced thermoplastic composite pipes (RTPs) have been widely used for oil and gas gathering and transportation. Polyvinylidene fluoride (PVDF) has the greatest potential as a thermoplastic liner of RTPs due to its excellent thermal and mechanical properties. However, permeation of gases is inevitable in the thermoplastic liner, which may lead to blister failure of the liner and damage the safe operation of the RTPs. In order to clarify the permeation behavior and obtain the permeation mechanism of the mixture gas (CH4/CO2/H2S) in PVDF at the normal service conditions, molecular simulations were carried out by combining the Grand Canonical Monte Carlo (GCMC) method and the Molecular Dynamics (MD) method. The simulated results showed that the solubility coefficients of gases increased with the decrease in temperature and the increase in pressure. The adsorption isotherms of all gases were consistent with the Langmuir model. The order of the adsorption concentration for different gases was H2S > CO2> CH4. The isosteric heats of gases at all the actual service conditions were much less than 42 kJ/mol, which indicated that the adsorption for all the gases belonged to the physical adsorption. Both of the diffusion and permeation coefficients increased with the increase in temperature and pressure. The diffusion belonged to Einstein diffusion and the diffusion coefficients of each gas followed the order of CH4 > CO2 > H2S. During the permeation process, the adsorption of gas molecules in PVDF exhibited selective aggregation, and most of them were adsorbed in the low potential energy region of PVDF cell. The mixed-gas molecules vibrated within the hole of PVDF at relatively low temperature and pressure. As the temperature and pressure increase, the gas molecules jumped into the neighboring holes occasionally and then dwelled in the holes, moving around their equilibrium positions.

Water Dissolved in a Variety of Polymers Studied by Molecular Dynamics Simulation and a Theory of Solutions

The performance of a polymer medium as a separation membrane is determined by the dissolution free energy ΔG and diffusion coefficient D of the permeant. In this work, ΔG and D of water are investigated with all-atom molecular dynamics simulation in a wide variety of polymer species in the amorphous state. The computed ΔG is shown to agree well with the experimental value for linear homopolymers, and the degrees of polymerization of the homopolymers do not affect ΔG when they are beyond ∼10. The copolymers of ethylene-vinylidene difluoride, ethylene-vinyl acetate, and ethylene-acrylamide are then examined by changing the repeating patterns of the constituent monomers in both the periodic and graft forms. It is found that ΔG is determined primarily by the overall compositions of the monomers and is not affected by the copolymerization topology (periodic or graft). The hydrophobicity of the copolymer is enhanced, furthermore, when the hydrophobicity and hydrophilicity of the ethylene and non-ethylene parts are well contrasted and those parts are fragmented along the polymer chain. According to the computed D, the diffusivity of water tends to be larger when the (co)polymer is more hydrophobic and ΔG is more positive. D is actually seen to vary by orders of magnitude with the polymer structures, while the effect of the polymer species on the water permeation is stronger for ΔG than for D.

Relations Between Dynamic Localization and Solute Diffusion in Polymers

Various public health concerns can arise from the unintended leaching of additives and impurities from polymeric medical devices or food packaging, which is directly related to each solute's diffusivity D. Both experimental and simulation methods can be used to quantify D, but slow diffusion at physiologic temperature in glassy polymers can render these approaches impractical. Here, we investigate a simulation approach with the potential to more rapidly calculate D. Specifically, we examine links between dynamic localization, characterized by the Debye-Waller factor, ⟨u²⟩, and D in a variety of polymer/solute systems using atomistic molecular dynamics (MD) simulations. Using short, high-temperature MD simulations to estimate D at physiologic temperature, we find that the relation ln D ∝ 1/⟨u²⟩ quantitatively predicts D for small solutes and produces an upper-bound estimate of D for larger solutes. Upper-bound estimates are useful in certain contexts, and we compare our results with another approach for determining upper bounds, the Piringer model, to show where each method may be useful. Then, we examine a modified relation where the Debye-Waller factor is rescaled by the mode coupling temperature T_c, which can produce better estimates of D if T_c is carefully chosen. Last, we compare our approach with several other models that relate temperature or localized dynamics with diffusivity. Although each of these approaches can be used to model D across wide temperature ranges using one or more adjustable parameters, none of them are truly predictive in glassy polymers. Further developments are needed to predict the optimal values of the adjustable parameters a priori.

The structure of CO2 and CH4 at the interface of a poly(urethane urea) oligomer model from the microscopic point of view

The world desperately needs new technologies and solutions for gas capture and separation. To make this possible, molecular modeling is applied here to investigate the structural, thermodynamic, and dynamical properties of a model for the poly(urethane urea) (PUU) oligomer model to selectively capture CO₂ in the presence of CH₄. In this work, we applied a well-known approach to derive atomic partial charges for atoms in a polymer chain based on self-consistent sampling using quantum chemistry and stochastic dynamics. The interactions of the gases with the PUU model were studied in a pure gas based system as well as in a gas mixture. A detailed structure characterization revealed high interaction of CO₂ molecules with the hard segments of the PUU. Therefore, the structural and energy properties explain the reasons for the greater CO₂ sorption than CH₄. We find that the CO₂ sorption is higher than the CH₄ with a selectivity of 7.5 at 298 K for the gas mixture. We characterized the Gibbs dividing surface for each system, and the CO₂ is confined for a long time at the gas-oligomer model interface. The simulated oligomer model showed performance above the 2008 Robeson's upper bound and may be a potential material for CO₂/CH₄ separation. Further computational and experimental studies are needed to evaluate the material.

Molecular dynamics simulation for membrane separation and porous materials: A current state of art review

Computational frameworks have been under specific attention within the last two decades. Molecular Dynamics (MD) simulations, identical to the other computational approaches, try to address the unknown question, lighten the dark areas of unanswered questions, to achieve probable explanations and solutions. Owing to their complex microporous structure on one side and the intricate biochemical nature of various materials used in the structure, separative membrane materials possess peculiar degrees of complications. More notably, as nanocomposite materials are often integrated into separative membranes, thin-film nanocomposites and porous separative nanocomposite materials could possess an additional level of complexity with regard to the nanoscale interactions brought to the structure. This critical review intends to cover the recent methods used to assess membranes and membrane materials. Incorporation of MD in membrane technology-related fields such as desalination, fuel cell-based energy production, blood purification through hemodialysis, etc., were briefly covered. Accordingly, this review could be used to understand the current extent of MD applications for separative membranes. The review could also be used as a guideline to use the proper MD implementation within the related fields.

Impact of Backbone Amide Substitution at the Meta- and Para-Positions on the Gas Barrier Properties of Polyimide

This study designed and synthesised a meta-amide-substituted dianiline monomer (m-DABA) as a stereoisomer of DABA, a previously investigated para-amide-substituted dianiline monomer. This new monomer was polymerised with pyromellitic dianhydride (PMDA) to prepare a polyimide film (m-DABPI) in a process similar to that employed in a previous study. The relationship between the substitution positions on the monomer and the gas barrier properties of the polyimide film was investigated via molecular simulation, wide-angle X-ray diffraction (WXRD), and positron annihilation lifetime spectroscopy (PALS) to gain deeper insights into the gas barrier mechanism. The results showed that compared with the para-substituted DABPI, the m-DABPI exhibited better gas barrier properties, with a water vapour transmission rate (WVTR) and an oxygen transmission rate (OTR) as low as 2.8 g·m⁻²·d⁻¹ and 3.3 cm³·m⁻²·d⁻¹, respectively. This was because the meta-linked polyimide molecular chains were more tightly packed, leading to a smaller free volume and lower molecular chain mobility. These properties are not conducive to the permeation of small molecules into the film; thus, the gas barrier properties were improved. The findings have significant implications for the structural design of high-barrier materials and could promote the development of flexible display technology.

Structure and Gas Barrier Properties of Polyimide Containing a Rigid Planar Fluorene Moiety and an Amide Group: Insights from Molecular Simulations

A novel diamine (FAPDA) bearing rigid planar fluorene and amide groups was successfully synthesized. Using such diamine and pyromellitic dianhydride (PMDA), a high-barrier polyimide (FAPPI) was obtained. FAPPI exhibits an outstanding gas barrier. Its water vapor transmission rate (WVTR) and oxygen transmission rate (OTR) are as low as 0.51 g·m^-2·day^-1 and 0.43 cm³·m^-2·day^-1, respectively. Additionally, FAPPI shows excellent thermal stability with a coefficient of thermal expansion (CTE) of 5.8 ppm·K^-1 and a glass transition temperature (T _g) of 416 °C. Molecular simulations, positron annihilation, and X-ray diffraction were utilized to gain insight on the microstructures for the enhanced barrier properties. Introducing fluorene moieties and amide groups improves the regularity and rigidity of molecular chains and increases interchain interaction of PI, resulting in low free volumes and decreased movement capacity of the chain. The low free volumes of FAPPI restrain the gas diffusivity and solubility. Meanwhile, the decreased chain movement reduces the diffusivity of gases. Consequently, barrier performances of FAPPI are improved. The polyimide possesses widespread application in the microelectronics packaging fields.

High-Barrier Polyimide Containing Carbazole Moiety: Synthesis, Gas Barrier Properties, and Molecular Simulations

A high-barrier polyimide (2,7-CPI) was synthesized through the polymerization of pyromellitic dianhydride (PMDA) and a novel diamine (2,7-CDA) containing carbazole moiety. The synthesized diamine and polyimide were fully characterized by elemental analyses, FTIR and NMR. The 2,7-CPI displays very attractive barrier performances, with oxygen transmission rate (OTR) and water vapor transmission rate (WVTR) low to 0.14 cm³·m⁻²·day⁻¹ and 0.05 g·m⁻²·day⁻¹, respectively. Meanwhile, 2,7-CPI also exhibits exceptional thermal stability with a glass transition temperature (T_g) of 467 °C, 5% weight-loss temperature (T_d5%) of 550 °C under N₂ and coefficient of thermal expansion (CTE) of 3.4 ppm/K. The barrier performances of 2,7-CPI are compared with those of a structural analogue (2,7-CPPI) and a typical polyimide (Kapton). Their barrier performances with respect to microstructure were investigated by molecular simulations, wide angle X-ray diffraction (WAXD), and positron annihilation lifetime spectroscopy (PALS). The results show that 2,7-CPI possesses better coplanar structure and more number of intermolecular hydrogen bonds among the three PIs, which result in tight chain packing and thereby high crystallinity, low free volume, and decreased chains mobility. That is, the high crystallinity and low free volume of 2,7-CPI reduce the diffusion and solubility of gases. Meanwhile, the poor chains mobility further decreases the gases diffusion. The reduced diffusion and solubility of gases consequently promote the improvement of barrier properties for 2,7-CPI. The polyimide has a wide application prospect in the field of flexible electronic packaging industries.

Molecular Simulations and Mechanistic Analysis of the Effect of CO2 Sorption on Thermodynamics, Structure, and Local Dynamics of Molten Atactic Polystyrene

A simulation strategy encompassing different scales was applied to the systematic study of the effects of CO₂ uptake on the properties of atactic polystyrene (aPS) melts. The analysis accounted for the influence of temperature between 450 and 550 K, polymer molecular weights (M _w) between 2100 and 31000 g/mol, and CO₂ pressures up to 20 MPa on the volumetric, swelling, structural, and dynamic properties of the polymer as well as on the CO₂ solubility and diffusivity by performing molecular dynamics (MD) simulations of the system in a fully atomistic representation. A hierarchical scheme was used for the generation of the higher M _w polymer systems, which consisted of equilibration at a coarse-grained level of representation through efficient connectivity-altering Monte Carlo simulations, and reverse-mapping back to the atomistic representation, obtaining the configurations used for subsequent MD simulations. Sorption isotherms and associated swelling effects were determined by using an iterative procedure that incorporated a series of MD simulations in the NPT ensemble and the Widom test particle insertion method, while CO₂ diffusion coefficients were extracted from long MD runs in the NVE ensemble. Solubility and diffusivity compared favorably with experimental results and with predictions of the Sanchez-Lacombe equation of state, which was reparametrized to capture the M _w dependence of polymer properties with greater accuracy. Structural features of the polymer matrix were correctly reproduced by the simulations, and the effects of gas concentration and M _w on structure and local dynamics were thoroughly investigated. In the presence of CO₂, a significant acceleration of the segmental dynamics of the polymer occurred, more pronouncedly at low M _w. The speed-up effect caused by the swelling agent was not limited to the chain ends but affected the whole chain in a similar fashion.

Molecular dynamics simulation of polystyrene copolymer with octyl short-chain branches in toluene

In this study, dimensional, conformational and dynamic behaviors of a short-chain branched styrene/1-octene copolymer chain with different 1-octene percentages, i.e., 0, 2, 4 and 6%, in toluene are investigated at the temperature of 298.15 K via molecular dynamics simulation. The chain dimensions and flexibility in the solvent are evaluated by calculating the radius of gyration (R_g), end-to-end distance (<r₀>), surface area (A_ch), and volume (V_ch) of the copolymer chain. The mean square displacement (MSD) and diffusivity coefficient for each copolymer chain are measured to determine its dynamic behavior and mobility in aromatic media. To consider the effect of increasing the 1-octene co-monomer percentage on the copolymer chain affinity to the solvent molecules, the interaction energy (E_int) and Flory-Huggins (FH) interaction parameter are calculated for each equilibrated solution model. The simulation results indicate that the co-monomer level increment in the copolymer structure reduces the chain R_g amount and its interaction with the solvent. The <r₀> of the chain increases up to 4% co-monomer content, while further co-monomer content decreases the <r₀> value. Additionally, the viscosity of the equilibrated dilute solutions is calculated via non-equilibrium molecular dynamics simulation (NEMD). Moreover, the steric hindrance of the copolymers and the solvent molecules capturing in the dilute solution is determined via radial distribution function (RDF) analysis. Helmholtz free energy and the system entropy changes are calculated to evaluate the tendency of the copolymer to the solvent molecules and its dilute solution irregularity, respectively. Graphical abstract The figure shows the variations trend of the poly(styrene-co-1-octene) chain dimensions in toluene aromatic solvent by increasing the 1-octene content (x), after the equilibration state. Red and blue colors represent the carbon atoms of the copolymer chain backbone and 1-octene side chains, respectively. The styrene rings and the hydrogen atoms of the chains were removed for better view.

Molecular dynamics simulation of three plastic additives' diffusion in polyethylene terephthalate

Accurate diffusion coefficient data of additives in a polymer are of paramount importance for estimating the migration of the additives over time. This paper shows how this diffusion coefficient can be estimated for three plastic additives [2-(2'-hydroxy-5'-methylphenyl) (UV-P), 2,6-di-tert-butyl-4-methylphenol (BHT) and di-(2-ethylhexyl) phthalate (DEHP)] in polyethylene terephthalate (PET) using the molecular dynamics (MD) simulation method. MD simulations were performed at temperatures of 293-433 K. The diffusion coefficient was calculated through the Einstein relationship connecting the data of mean-square displacement at different times. Comparison of the diffusion coefficients simulated by the MD simulation technique, predicted by the Piringer model and experiments, showed that, except for a few samples, the MD-simulated values were in agreement with the experimental values within one order of magnitude. Furthermore, the diffusion process for additives is discussed in detail, and four factors - the interaction energy between additive molecules and PET, fractional free volume, molecular shape and size, and self-diffusion of the polymer - are proposed to illustrate the microscopic diffusion mechanism. The movement trajectories of additives in PET cell models suggested that the additive molecules oscillate slowly rather than hopping for a long time. Occasionally, when a sufficiently large hole was created adjacently, the molecule could undergo spatial motion by jumping into the free-volume hole and consequently start a continuous oscillation and hop. The results indicate that MD simulation is a useful approach for predicting the microstructure and diffusion coefficient of plastic additives, and help to estimate the migration level of additives from PET packaging.

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.

Experimental and simulation-based investigation of He, Ne and Ar irradiation of polymers for ion microscopy

Secondary ion mass spectrometry (SIMS) on the helium ion microscope (HIM) promises higher lateral resolution than on classical SIMS instruments. However, full advantage of this new technique can only be obtained when the interaction of He(+) or Ne(+) primary ions with the sample is fully controlled. In this work we investigate how He(+) and Ne(+) bombardment influences roughness formation and preferential sputtering for polymer samples and how they compare to Ar(+) primary ions used in classical SIMS by combining experimental techniques with Molecular Dynamics (MD) simulations and SD_TRIM_SP modelling. The results show that diffusion coefficients for He, Ne and Ar in polymers are sufficiently high to prevent any accumulation of rare gas atoms in the polymers which could lead to some swelling and bubble formation. Roughness formation was also not observed. Preferential sputtering is more of a problem, with enrichment of carbon up to surface concentrations above 80%. In general, the preferential sputtering is largely depending on the primary ion species and the impact energies. For He(+) bombardment, it is more of an issue for low keV impact energies and for the heavier primary ion species the preferential sputtering is sample dependent. For He(+) steady state conditions are reached for fluences much higher than 10(18) ions/cm(2). For Ne(+) and Ar(+), the transient regime extends up to fluences of 10(17)-10(18) ions/cm(2). Hence, preferential sputtering needs to be taken into account when interpreting images recorded under He(+) or Ne(+) bombardment on the HIM.

Prediction and validation of diffusion coefficients in a model drug delivery system using microsecond atomistic molecular dynamics simulation and vapour sorption analysis

Diffusion of small to medium sized molecules in polymeric medical device materials underlies a broad range of public health concerns related to unintended leaching from or uptake into implantable medical devices. However, obtaining accurate diffusion coefficients for such systems at physiological temperature represents a formidable challenge, both experimentally and computationally. While molecular dynamics simulation has been used to accurately predict the diffusion coefficients, D, of a handful of gases in various polymers, this success has not been extended to molecules larger than gases, e.g., condensable vapours, liquids, and drugs. We present atomistic molecular dynamics simulation predictions of diffusion in a model drug eluting system that represent a dramatic improvement in accuracy compared to previous simulation predictions for comparable systems. We find that, for simulations of insufficient duration, sub-diffusive dynamics can lead to dramatic over-prediction of D. We present useful metrics for monitoring the extent of sub-diffusive dynamics and explore how these metrics correlate to error in D. We also identify a relationship between diffusion and fast dynamics in our system, which may serve as a means to more rapidly predict diffusion in slowly diffusing systems. Our work provides important precedent and essential insights for utilizing atomistic molecular dynamics simulations to predict diffusion coefficients of small to medium sized molecules in condensed soft matter systems.

Sorption and diffusion of carbon dioxide and nitrogen in poly(methyl methacrylate)

Molecular dynamics simulations are performed to determine the solubility and diffusion coefficient of carbon dioxide and nitrogen in poly(methyl methacrylate) (PMMA). The solubilities of CO2 in the polymer are calculated employing our grand canonical ensemble simulation method, fixing the target excess chemical potential of CO2 in the polymer and varying the number of CO2 molecules in the polymer matrix till establishing equilibrium. It is shown that the calculated sorption isotherms of CO2 in PMMA, employing this method well agrees with experiment. Our results on the diffusion coefficients of CO2 and N2 in PMMA are shown to obey a common hopping mechanism. It is shown that the higher solubility of CO2 than that of N2 is a consequence of more attractive interactions between the carbonyl group of polymer and the sorbent. While the residence time of CO2 beside the carbonyl group of polymer is about three times higher than that of N2, the diffusion coefficient of CO2 in PMMA is higher than that of N2. The higher diffusion coefficient of CO2, compared to N2, in PMMA is shown to be due to the higher (≈3 times) swelling of polymer upon CO2 uptake.