Molecular Modeling Investigations of Sorption and Diffusion of Small Molecules in Glassy Polymers

Polyimide Polymer Simulations through Coarse-Grained Modeling: Prediction of Structure, Physical Properties, and Gas Separation Properties

In this study, we introduce a set of coarse-grained (CG) force field parameters for simulating a series of 6FDA-based polyimides. Utilizing atomistic descriptors, we developed CG models that accurately predict the specific volume of the polymers under investigation. Our findings suggest that certain parameters, particularly those associated with specific diamines, can be employed to predict properties such as density using a multiple linear regression. Our study further explores the halogenation of diamines and proposes methods for estimating intermolecular interaction parameters. Our calculations refer to various structural properties, including the radius of gyration, end-to-end distance, glass transition temperature, and diffusion coefficients. Utilizing the newly developed CG force field parameters, we conducted gas separation simulations for 6FDA-DAM polyimide, particularly to predict both sorption- and diffusion-separation mechanisms within the polymer. These simulations provided excellent agreement with experimental data on solubility, diffusion, and permeability selectivity for CO2/CH4, O2/N2, and propylene/propane. The results contribute significantly to our understanding of polyimide behavior, and the parameters proposed here offer a promising tool for the development of new materials with tailored properties for targeted applications.

Safety challenges of mechanically recycled polyethylene terephthalate for food contact materials: A review

The global market for polyethylene terephthalate (PET) production and consumption is gradually expanding, with a large amount of plastic wasted. PET is the primary plastic used for direct food contact and the most commonly used mono-material. As it is an appreciable candidate for recycling, PET has emerged as a core concept in the global circular economy of plastics in terms of mechanical and chemical recycling. Chemical recycling is a depolymerization method involving gasification, pyrolysis, and depolymerization, whereas mechanical recycling is a physical process involving sorting, washing, drying, grinding, and remelting. Mechanical recycling, which has the benefits of lower carbon emissions and relative affordability, is more cost- and time-efficient than chemical recycling. Additionally, various safety assessment methods have been established for applying recycled PET (rPET) to food-contact materials (FCMs). Contaminants excised from mechanically recycled PET are reviewed according to different detection methods and types of rPET. The results suggest that various contaminants in rPET produced during each PET recycling process can be monitored using overall migration and specific migration methods. Additionally, the super-cleaning process with filtration is crucial in PET mechanical recycling in terms of the application of rPET as FCMs and the removal of residual migrants. Therefore, this review provides a comprehensive analysis of the safety evaluation of rPET from mechanical recycling, focusing on migration testing and detection methods. It highlights the reliability and limitations of the current regime and identifies gaps in mechanical PET recycling for future research, contributing to the advancement of rPET safety and sustainable applications.

Sorption and release of small molecules in PDMS and COC for Organs on chip

Accurate risk assessment in drug development is crucial, as conventional in vitro and in vivo models often fail to predict human-specific responses. Organs on chips offer a promising alternative, but widespread use of polydimethylsiloxane introduces challenges due to its sorption of small lipophilic molecules, distorting pharmacokinetic and pharmacodynamic data. Cyclic olefin copolymer, a chemically stable alternative with minimal sorption, has emerged as a potential solution. This study investigates the sorption behavior of seven pharmaceutically active compounds in microfluidic devices and washout of these compounds, using high-performance liquid chromatography-mass spectrometry to evaluate recovery of compounds. Lipophilic molecules exhibited substantial sorption in polydimethylsiloxane and lower retention in cyclic olefin copolymer. Imipramine (logP = 4.80) decreased from 100 µM to 0.0384 µM for polydimethylsiloxane and 31.5 µM for cyclic olefin copolymer after 24 h incubation. Sorption was governed by multiple factors - lipophilicity and rotatable bond count were critical for both materials, hydrogen bond acceptors and molecular weight played a larger role in cyclic olefin copolymer, whereas topological polar surface area was critical for polydimethylsiloxane. Washout studies revealed that polydimethylsiloxane retains lipophilic compounds through bulk absorption, causing slow release, while cyclic olefin copolymer facilitated easier desorption. The cumulative sum of the first 5 h washout of loperamide (logP = 5.13) is 37.8% for polydimethylsiloxane and 71.5% for cyclic olefin copolymer. These findings highlight the importance of material selection and molecular properties in minimizing sorption and ensuring reliable experimental outcomes especially in microfluidic systems with distinctly different surface to volume ratios to other models.

Measuring Oxygen Solubility in Amorphous and Semicrystalline Polyolefins Using Test Particle Insertion: A Comparative Study of Polyethylene and Isotactic Polypropylene

The test particle insertion method is used to study the solubility of oxygen in two commonly used polymers: polyethylene (PE) and isotactic polypropylene (iPP). Amorphous samples for both polymers were prepared by means of Monte Carlo and molecular dynamics simulations, and the oxygen solubility was measured across different temperatures. The solubility-temperature dependence for iPP proved to be nonmonotonic due to the interplay between binding and reorganizational enthalpy, while for PE, it appeared to be monotonic based on the available data in the studied temperature range. A broad comparison to experiments and simulations is included. Further oxygen insertions were also performed in semicrystalline PE and iPP samples at ambient temperature, and the obtained values were compared to a linear relationship which correlates the solubility in the purely amorphous phase with the solubility in the crystalline phase. The solubility of PE closely follows the linear relationship, while iPP exhibits some divergence. All the semicrystalline samples were previously annealed at elevated temperatures for long periods (a few μs), and a strong effect of annealing was observed on the structure and the solubility of iPP. A well-developed iPP lamellar structure emerged at longer annealing times, while PE develops that structure already in the early crystallization stages. The solubility of semicrystalline iPP samples with lamellar morphology exhibited better agreement with extrapolated solubility values of the amorphous state─the extrapolation was made using a linear relationship connecting solubility in the purely amorphous phase and solubility in mixed phases (amorphous and crystalline). Results on the correlation of the solubility with the local structural ordering are also present.

GMXPolymer: a generated polymerization algorithm based on GROMACS

CONTEXT

This work introduces a method for generating generalized structures of amorphous polymers using simulated polymerization and molecular dynamics equilibration, with a particular focus on amorphous polymers. The techniques and algorithms used in this method are described in the main text, and example input scripts are provided for the GMXPolymer code, which is based on the GROMACS molecular dynamics package. To demonstrate the efficacy of our method, we apply it to different glassy polymers exhibiting varying degrees of functionality, polarity, and rigidity. The reliability of the method is validated by comparing simulation results with experimental data in various structural and thermal properties, both of which show excellent agreement.

METHODS

This work implements the GMXPolymer simulated polymerization algorithm on the GROMACS program. GMXPolymer code controls the main polymerization loop. The energy minimizations and molecular dynamics simulations use the GROMACS program called by the GMXPolymer code. A new ITP file is generated when a new bond is formed, and the necessary additions to the ITP file are made to include new bonds, angles, and dihedrals. In preparing the ITP file of the monomer, the charge of the reactive atom must be modified before the code runs so that it is a correct value after bonding.

Study of oxygen transport in glassy polymers on a nanometer length scale utilizing the kinetic Monte Carlo simulations

Diffusion-controlled deactivation of excited phenanthrene and oxidation of triplet aryl-nitrene by molecular oxygen were used to determine the energetics of oxygen jump rates in the set of glassy polymers: poly(methyl methacrylate), poly(n-butyl methacrylate), polycarbonate, polystyrene, and polysulfone. To interpret experimental results, a simple model based on the transition state theory of diffusion jump has been used. The kinetic Monte Carlo simulations of phenanthrene deactivation and nitrene oxidation were carried out in a cubic lattice that modeled a polymer matrix. The bonds of the lattice were assigned to be activation barriers for the diffusion jumps of oxygen molecules from one site of the lattice to another. The standard deviation, σ‡, and spatial correlation length, rc, of the free energy of diffusion jump have been determined. It is shown that the spatial correlation of oxygen jump rates on a nanometer scale and the entropic nature of the dynamic heterogeneity are common features of all the studied polymers.

Comparative Study of Sorption Phenomena Between Three Medications and Syringes Made of Cyclic Olefin Copolymer or Polypropylene

INTRODUCTION

Medical syringes are widely used in hospitals to store and administer drugs, and the contact time between the drugs and these syringes can vary from a few minutes to several weeks like for pharmaceutical preparations. The aim of this comparative study was to evaluate the potential sorption phenomena occurring between three drugs (paracetamol, diazepam and insulin aspart) and polypropylene syringes (PP) or syringes made of Cyclic Olefin Copolymer (COC).

MATERIALS AND METHODS

50 mL 3-part syringes made of either COC with crosslinked silicone on the barrel inner surface (COC-CLS) and a bromobutyl plunger seal, or PP lubricated with silicone oil (PP-SOL) with a polyisoprene plunger seal were used.

RESULTS

COC-CLS syringes induced less sorption of diazepam and insulin than PP-SOL syringes and the plunger seal material seemed to be the main cause of these interactions. An alkalinization of the medications in contact with the PP-SOL syringes was observed. It could be caused by leachable compounds and should be investigated further.

CONCLUSION

This work shows once again that it is essential to consider content-container interactions to help improve the safe use of parenteral drugs.

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.

Anomalous Diffusion of Helium and Neon in Low-Density Silica Glass

The diffusion properties of low-density non-porous silica glasses (expanded silica glasses) were researched with the aim of searching for the molecular structure of membrane materials intended for the effective separation of helium-neon gas mixtures. It has been shown on a large number (84) of computer models of such glasses that there are molecular structures of silica in which various helium and neon diffusion mechanisms are simultaneously implemented: superdiffusion for helium and subdiffusion for neon. This makes it possible to significantly (by 3-5 orders of magnitude) increase the helium permeability of such glasses at room temperature and maintain a high selectivity for the separation of helium and neon (at the level of 104-105) at the same time.

Prediction of Adsorption and Diffusion of Shale Gas in Composite Pores Consisting of Kaolinite and Kerogen using Molecular Simulation

Natural gas production from shale formations is one of the most recent and fast growing developments in the oil and gas industry. The accurate prediction of the adsorption and transport of shale gas is essential for estimating shale gas production capacity and improving existing extractions. To realistically represent heterogeneous shale formations, a composite pore model was built from a kaolinite slit mesopore hosting a kerogen matrix. Moreover, empty slabs (2, 3, and 4 nm) were added between the kerogen matrix and siloxane surface of kaolinite. Using Grand-Canonical Monte Carlo (GCMC) and molecular dynamics (MD) simulations, the adsorption and diffusion of pure methane, pure ethane, and a shale gas mixture were computed at various high pressures (100, 150, and 250 atm) and temperature of 298.15 K. The addition of an inner slit pore was found to significantly increase the excess adsorption of methane, as a pure component and in the shale gas mixture. The saturation of the composite pore with methane was observed to be at a higher pressure compared to ethane. The excess adsorption of carbon dioxide was not largely affected by pressure, and the local number density profile showed its strong affinity to kerogen micropores and the hydroxylated gibbsite surface under all conditions and pore widths. Lateral diffusion coefficients were found to increase with increasing the width of the empty slab inside the composite pore. Statistical errors of diffusion coefficients were found to be large for the case of shale gas components present at low composition. A larger composite pore configuration was created to investigate the diffusion of methane in different regions of the composite pore. The calculated diffusion coefficients and mean residence times were found to be indicative of the different adsorption mechanisms occurring inside the pore.

Integrating Machine Learning in the Coarse-Grained Molecular Simulation of Polymers

Machine learning (ML) is having an increasing impact on the physical sciences, engineering, and technology and its integration into molecular simulation frameworks holds great potential to expand their scope of applicability to complex materials and facilitate fundamental knowledge and reliable property predictions, contributing to the development of efficient materials design routes. The application of ML in materials informatics in general, and polymer informatics in particular, has led to interesting results, however great untapped potential lies in the integration of ML techniques into the multiscale molecular simulation methods for the study of macromolecular systems, specifically in the context of Coarse Grained (CG) simulations. In this Perspective, we aim at presenting the pioneering recent research efforts in this direction and discussing how these new ML-based techniques can contribute to critical aspects of the development of multiscale molecular simulation methods for bulk complex chemical systems, especially polymers. Prerequisites for the implementation of such ML-integrated methods and open challenges that need to be met toward the development of general systematic ML-based coarse graining schemes for polymers are discussed.

How Segmental Dynamics and Mesh Confinement Determine the Selective Diffusivity of Molecules in Cross-Linked Dense Polymer Networks

The diffusion of molecules ("penetrants") of variable size, shape, and chemistry through dense cross-linked polymer networks is a fundamental scientific problem broadly relevant in materials, polymer, physical, and biological chemistry. Relevant applications include separation membranes, barrier materials, drug delivery, and nanofiltration. A major open question is the relationship between transport, thermodynamic state, and penetrant and polymer chemical structure. Here we combine experiment, simulation, and theory to unravel these competing effects on penetrant transport in rubbery and supercooled polymer permanent networks over a wide range of cross-link densities, size ratios, and temperatures. The crucial importance of the coupling of local penetrant hopping to polymer structural relaxation and the secondary importance of mesh confinement effects are established. Network cross-links strongly slow down nm-scale polymer relaxation, which greatly retards the activated penetrant diffusion. The demonstrated good agreement between experiment, simulation, and theory provides strong support for the size ratio (penetrant diameter to the polymer Kuhn length) as a key variable and the usefulness of coarse-grained simulation and theoretical models that average over Angstrom scale structure. The developed theory provides an understanding of the physical processes underlying the behaviors observed in experiment and simulation and suggests new strategies for enhancing selective polymer membrane design.

Novel Injectable Fluorescent Polymeric Nanocarriers for Intervertebral Disc Application

Damage to intervertebral discs (IVD) can lead to chronic pain and disability, and no current treatments can fully restore their function. Some non-surgical treatments have shown promise; however, these approaches are generally limited by burst release and poor localization of diverse molecules. In this proof-of-concept study, we developed a nanoparticle (NP) delivery system to efficiently deliver high- and low-solubility drug molecules. Nanoparticles of cellulose acetate and polycaprolactone-polyethylene glycol conjugated with 1-oxo-1H-pyrido [2,1-b][1,3]benzoxazole-3-carboxylic acid (PBC), a novel fluorescent dye, were prepared by the oil-in-water emulsion. Two drugs, a water insoluble indomethacin (IND) and a water soluble 4-aminopyridine (4-AP), were used to study their release patterns. Electron microscopy confirmed the spherical nature and rough surface of nanoparticles. The particle size analysis revealed a hydrodynamic radius ranging ~150-162 nm based on dynamic light scattering. Zeta potential increased with PBC conjugation implying their enhanced stability. IND encapsulation efficiency was almost 3-fold higher than 4-AP, with release lasting up to 4 days, signifying enhanced solubility, while the release of 4-AP continued for up to 7 days. Nanoparticles and their drug formulations did not show any apparent cytotoxicity and were taken up by human IVD nucleus pulposus cells. When injected into coccygeal mouse IVDs in vivo, the nanoparticles remained within the nucleus pulposus cells and the injection site of the nucleus pulposus and annulus fibrosus of the IVD. These fluorescent nano-formulations may serve as a platform technology to deliver therapeutic agents to IVDs and other tissues that require localized drug injections.

Adsorption and membrane separation for removal and recovery of volatile organic compounds

Volatile organic compounds (VOCs) are a crucial kind of pollutants in the environment due to their obvious features of severe toxicity, high volatility, and poor degradability. It is particularly urgent to control the emission of VOCs due to the persistent increase of concentration and the stringent regulations. In China, clear directions and requirements for reduction of VOCs have been given in the "national plan on environmental improvement for the 13th Five-Year Plan period". Therefore, the development of efficient technologies for removal and recovery of VOCs is of great significance. Recovery technologies are favored by researchers due to their advantages in both recycling VOCs and reducing carbon emissions. Among them, adsorption and membrane separation processes have been extensively studied due to their remarkable industrial prospects. This overview was to provide an up-to-date progress of adsorption and membrane separation for removal and recovery of VOCs. Firstly, adsorption and membrane separation were found to be the research hotspots through bibliometric analysis. Then, a comprehensive understanding of their mechanisms, factors, and current application statuses was discussed. Finally, the challenges and perspectives in this emerging field were briefly highlighted.

Upcoming Transformations in Integrated Energy/Chemicals Sectors: Some Challenges and Several Opportunities

The sociopolitical events over the past few years led to transformative changes in both the energy and chemical sectors. One of the most evident consequences of these events is the significant focus on sustainability. In fact, rather than an engaging discussion within elite social circles, the search for sustainability is now one of the hard requirements investors impose on companies. The concept of sustainability itself has developed since its inception, and now it encompasses environmental as well as socioeconomic aspects. The major players in the energy and chemical sectors seem to embrace these changes and the related challenges; in most cases, tangible ambitious goals have been proposed. For example, bp aims "to become a net zero company by 2050 or sooner, and to help the world get to net zero". Although tragic events such as the war in Ukraine directly affect global supply chains, leading to some reconsiderations in medium-term industrial and political strategies, trends and public demands seem determined to pursue ambitious sustainable goals, as tangible as the European Union's "Fit for 55" climate package, approved on May 12, 2022, which effectively bans internal combustion engines for new passenger cars and light commercial vehicles from 2035. These trends will likely lead to profound changes in both the chemical and energy sectors. While some predictions may miss the target, speculating about upcoming challenges and opportunities could help us prepare for the future. This is the purpose of this brief Perspective.

A systematic review of the molecular simulation of hybrid membranes for performance enhancements and contaminant removals

Number of research on molecular simulation and design has emerged recently but there is currently a lack of review to present these studies in an organized manner to highlight the advances and feasibility. This paper aims to review the development, structural, physical properties and separation performance of hybrid membranes using molecular simulation approach. The hybrid membranes under review include ionic liquid membrane, mixed matrix membrane, and functionalized hybrid membrane for understanding of the transport mechanism of molecules through the different structures. The understanding of molecular interactions, and alteration of pore sizes and transport channels at atomistic level post incorporation of different components in hybrid membranes posing impact to the selective transport of desired molecules are also covered. Incorporation of molecular simulation of hybrid membrane in related fields such as carbon dioxide (CO₂) removal, wastewater treatment, and desalination are also reviewed. Despite the limitations of current molecular simulation methodologies, i.e., not being able to simulate the membrane operation at the actual macroscale in processing plants, it is still able to demonstrate promising results in capturing molecule behaviours of penetrants and membranes at full atomic details with acceptable separation performance accuracy. From the review, it was found that the best performing ionic liquid membrane, mixed matrix membrane and functionalized hybrid membrane can enhance the performance of pristine membrane by 4 folds, 2.9 folds and 3.3 folds, respectively. The future prospects of molecular simulation in hybrid membranes are also presented. This review could provide understanding to the current advancement of molecular simulation approach in hybrid membranes separation. This could also provide a guideline to apply molecular simulation in the related sectors.

Modelling Sorption and Transport of Gases in Polymeric Membranes across Different Scales: A Review

Professor Giulio C. Sarti has provided outstanding contributions to the modelling of fluid sorption and transport in polymeric materials, with a special eye on industrial applications such as membrane separation, due to his Chemical Engineering background. He was the co-creator of innovative theories such as the Non-Equilibrium Theory for Glassy Polymers (NET-GP), a flexible tool to estimate the solubility of pure and mixed fluids in a wide range of polymers, and of the Standard Transport Model (STM) for estimating membrane permeability and selectivity. In this review, inspired by his rigorous and original approach to representing membrane fundamentals, we provide an overview of the most significant and up-to-date modeling tools available to estimate the main properties governing polymeric membranes in fluid separation, namely solubility and diffusivity. The paper is not meant to be comprehensive, but it focuses on those contributions that are most relevant or that show the potential to be relevant in the future. We do not restrict our view to the field of macroscopic modelling, which was the main playground of professor Sarti, but also devote our attention to Molecular and Multiscale Hierarchical Modeling. This work proposes a critical evaluation of the different approaches considered, along with their limitations and potentiality.

Accelerated Diffusion Following Deprotection in Chemically Amplified Resists

Polymeric chemically amplified resists (CARs) are critical materials for high-throughput lithographic processes. A photoactivated acid-anion catalyst changes the polymer's solubility via a deprotection reaction, which enables pattern development through selective dissolution. To capture observed reaction kinetics, reaction-diffusion models employ a catalyst diffusivity that is accelerated by reaction. However, the microscopic origin and factors contributing to this phenomena remain unclear. Herein, we employ detailed atomistic molecular dynamics simulations to examine the impact of protecting group removal and material relaxation on catalyst mobility. We report data on polymer density, catalyst dispersion, excess free volume, and segmental dynamics with increasing time/extent of deprotection. We then propose simple kinetic Monte Carlo algorithms that can describe both molecular dynamics simulations of deprotection reactions and experimental data.

Temperature Effects in Flexible Adsorption Processes for Amorphous Microporous Polymers

A collection of atomistic molecular simulations is reported that illustrate the impact of adsorption temperature on species uptake and adsorbate-induced structural rearrangement for amorphous polymers of intrinsic microporosity. Temperature-sensitive structural rearrangement is evaluated by contrasting two methods: standard grand canonical Monte Carlo simulations using a rigid framework approximation and a combined Monte Carlo/molecular dynamics approach that fully incorporates framework flexibility. We report single-component gas phase adsorption isotherms for CH₄, C₂H₄, C₂H₆, C₃H₆, C₃H₈, and CO₂ across a temperature range of 250-400 K for models of an archetypal polymer of intrinsic microporosity, PIM-1. A quadratic model is presented that captures two main mechanisms of temperature-dependent adsorption-induced deformation of PIM-1 up to a relative swelling of 1.15: thermal expansion and an increased propensity to swell as a function of species uptake. Two case studies are reported that highlight the critical role of operating temperature in industrial storage and separation applications. The first study focuses on methane storage and delivery applications using a pressure-temperature swing adsorption application (PTSA). We demonstrate that larger working capacities are accompanied by increased volumetric strain between adsorption-desorption steps. The second case study considers PIM-1 as an adsorbent to separate an exemplar ternary syngas mixture at operating temperatures ranging 300-550 K. A temperature threshold of ∼400 K is identified, beyond which adsorption-induced PIM-1 swelling is negligible and the solubility selectivity-loading curve transitions to exhibiting a nearly linear relationship.

Atomistic Modeling of Hydrogen and Oxygen Solubility in Semicrystalline PA-6 and HDPE Materials

Hydrogen is a clean and sustainable energy carrier which plays a major role in the transition of the global energy market to a less fossil fuel dependent future. Polymer-based materials are crucial in the production, storage, transportation, and energy extraction of hydrogen. More insights in the hydrogen-polymers interactions are required to guide material design and product development, especially for hydrogen solubility in polymers, which is crucial in many applications. The current study aims at rationalizing the determining factors of hydrogen solubility in two relevant polymers: polyamide-6 (PA-6) and high density polyethylene (HDPE). Based on atomistic molecular dynamics simulations and experimental data, we have reached several conclusions related to hydrogen and oxygen solubility in these two polymers. The crystal phases of PA-6 and HDPE are impenetrable to hydrogen and oxygen at elevated pressures, despite the small molecular size of hydrogen and oxygen. The practical implication for gas barrier applications is that polymer crystals act as impermeable obstacles and gas migration takes place primarily in the amorphous phase. Experimental hydrogen and oxygen solubilities in PA-6 and HDPE at elevated pressures can be predicted in a semiquantitative manner by molecular simulations. The discrepancies between experimental and predicted values could be attributed to neglect of the effect of crystal regions on the amorphous polymer domains. Although hydrogen is smaller than oxygen, it has been experimentally observed that hydrogen has a lower solubility in PA-6 and HDPE than oxygen. This observation has been confirmed by molecular simulations and attributed to the more favorable energetic interactions of oxygen with PA-6 and PE than of hydrogen. These interactions dominate the solubility behavior over the distribution of the accessible volume in the polymers.

Single-gas and mixed-gas permeation of N2/CH4 in thermally-rearranged TR-PBO membranes and their 6FDA-bisAPAF polyimide precursor studied by molecular dynamics simulations

High-performance polymers with polybenzoxazole (PBO) structures, formed via thermal rearrangement (TR) of aromatic polyimide precursors, have been developed for gas separation applications. The present work compares the transport of N₂ and CH₄ in a 6FDA-bisAPAF polyimide precursor and in its TR-PBO derivative using molecular dynamics (MD) simulations. The modelling closely mimicked the experimental approach by transforming a 6FDA-bisAPAF atomistic model into its corresponding TR-PBO structure via a specific algorithm. The densities and void spaces of both precursor and TR polymers were found to compare well to experimental data. An iterative technique was used to obtain the single-gas sorption isotherms of N₂ and CH₄ at 338.5 K in both polymers over a range of feed pressures up to and exceeding 65 bar. CH₄ was systematically found to be more soluble than N₂. Solubilities in both matrices were quite similar with those in TR-PBO being slightly higher due to its larger fraction of significant volume. Volume dilation analyses confirmed a higher resistance to plasticization for TR-PBO. Extended single-gas N₂ and CH₄ simulations and 2 : 1 binary CH₄/N₂ mixed-gas simulations were then conducted in both matrices at 338.5 K and at a pressure of ∼65 bar corresponding to natural gas processing conditions. Mixed-gas sorption was modelled using a modification of the aforementioned iterative method, which fixed the pressure and iterated to convergence the number of molecules of each type of penetrant. The gas diffusion coefficients were estimated using the Trajectory-Extending Kinetic Monte Carlo (TEKMC) procedure. As found experimentally, significantly higher diffusivities and permeabilities were observed in the TR polymer, which led to a slightly lower ideal N₂/CH₄ permselectivity for TR-PBO (∼2.6) when compared to its 6FDA-bisAPAF precursor (∼3.8). However, both models showed a reduced N₂/CH₄ separation efficiency under 2 : 1 binary CH₄/N₂ mixed-gas conditions bordering on the loss of selectivity. For 6FDA-bisAPAF, both permeabilities decreased in the mixed-gas case, but more for N₂ than for CH₄. For TR-PBO, the permeability of the faster N₂ decreased while the permeability of the slower CH₄ increased under mixed-gas conditions. This confirms that single-gas simulations are not sufficient for the prediction of the actual mixed-gas permselectivity behaviour in such polymers.

Predicting Solute Diffusivity in Polymers Using Time-Temperature Superposition

We demonstrate a novel application of the time-temperature superposition (TTS) principle to predict solute diffusivity D in glassy polymers using atomistic molecular dynamics simulations. Our TTS approach incorporates the Debye-Waller factor ⟨u²⟩, a measure of solute caging, along with concepts from thermodynamic scaling methods, allowing us to balance contributions to the dynamics from temperature and ⟨u²⟩ using adjustable parameters. Our approach rescales the solute mean-squared displacement curves at several temperatures into a master curve that approximates the diffusive dynamics at a reference temperature, effectively extending the simulation time scale from nanoseconds to seconds and beyond. With a set of "universal" parameters, this TTS approach predicts D with reasonable accuracy in a broad range of polymer/solute systems. Using TTS greatly reduces the computational cost compared to standard MD simulations. Thus, our method offers a means to rapidly and routinely provide order-of-magnitude estimates of D using simulations.

Molecular Characterization of Membrane Gas Separation under Very High Temperatures and Pressure: Single- and Mixed-Gas CO2/CH4 and CO2/N2 Permselectivities in Hybrid Networks

This work illustrates the potential of using atomistic molecular dynamics (MD) and grand-canonical Monte Carlo (GCMC) simulations prior to experiments in order to pre-screen candidate membrane structures for gas separation, under harsh conditions of temperature and pressure. It compares at 300 °C and 400 °C the CO₂/CH₄ and CO₂/N₂ sieving properties of a series of hybrid networks based on inorganic silsesquioxanes hyper-cross-linked with small organic PMDA or 6FDA imides. The inorganic precursors are the octa(aminopropyl)silsesquioxane (POSS), which degrades above 300 °C, and the octa(aminophenyl)silsesquioxane (OAPS), which has three possible meta, para or ortho isomers and is expected to resist well above 400 °C. As such, the polyPOSS-imide networks were tested at 300 °C only, while the polyOAPS-imide networks were tested at both 300 °C and 400 °C. The feed gas pressure was set to 60 bar in all the simulations. The morphologies and densities of the pure model networks at 300 °C and 400 °C are strongly dependent on their precursors, with the amount of significant free volume ranging from ~2% to ~20%. Since measurements at high temperatures and pressures are difficult to carry out in a laboratory, six isomer-specific polyOAPS-imides and two polyPOSS-imides were simulated in order to assess their N₂, CH₄ and CO₂ permselectivities under such harsh conditions. The models were first analyzed under single-gas conditions, but to be closer to the real processes, the networks that maintained CO₂/CH₄ and CO₂/N₂ ideal permselectivities above 2 were also tested with binary-gas 90%/10% CH₄/CO₂ and N₂/CO₂ feeds. At very high temperatures, the single-gas solubility coefficients vary in the same order as their critical temperatures, but the differences between the penetrants are attenuated and the plasticizing effect of CO₂ is strongly reduced. The single-gas diffusion coefficients correlate well with the amount of available free volume in the matrices. Some OAPS-based networks exhibit a nanoporous behavior, while the others are less permeable and show higher ideal permselectivities. Four of the networks were further tested under mixed-gas conditions. The solubility coefficient improved for CO₂, while the diffusion selectivity remained similar for the CO₂/CH₄ pair and disappeared for the CO₂/N₂ pair. The real separation factor is, thus, mostly governed by the solubility. Two polyOAPS-imide networks, i.e., the polyorthoOAPS-PMDA and the polymetaOAPS-6FDA, seem to be able to maintain their CO₂/CH₄ and CO₂/N₂ sieving abilities above 2 at 400 °C. These are outstanding performances for polymer-based membranes, and consequently, it is important to be able to produce isomer-specific polyOAPS-imides for use as gas separation membranes under harsh conditions.

In Silico Prediction of Food Properties: A Multiscale Perspective

Several open software packages have popularized modeling and simulation strategies at the food product scale. Food processing and key digestion steps can be described in 3D using the principles of continuum mechanics. However, compared to other branches of engineering, the necessary transport, mechanical, chemical, and thermodynamic properties have been insufficiently tabulated and documented. Natural variability, accented by food evolution during processing and deconstruction, requires considering composition and structure-dependent properties. This review presents practical approaches where the premises for modeling and simulation start at a so-called “microscopic” scale where constituents or phase properties are known. The concept of microscopic or ground scale is shown to be very flexible from atoms to cellular structures. Zooming in on spatial details tends to increase the overall cost of simulations and the integration over food regions or time scales. The independence of scales facilitates the reuse of calculations and makes multiscale modeling capable of meeting food manufacturing needs. On one hand, new image-modeling strategies without equations or meshes are emerging. On the other hand, complex notions such as compositional effects, multiphase organization, and non-equilibrium thermodynamics are naturally incorporated in models without linearization or simplifications. Multiscale method’s applicability to hierarchically predict food properties is discussed with comprehensive examples relevant to food science, engineering and packaging. Entropy-driven properties such as transport and sorption are emphasized to illustrate how microscopic details bring new degrees of freedom to explore food-specific concepts such as safety, bioavailability, shelf-life and food formulation. Routes for performing spatial and temporal homogenization with and without chemical details are developed. Creating a community sharing computational codes, force fields, and generic food structures is the next step and should be encouraged. This paper provides a framework for the transfer of results from other fields and the development of methods specific to the food domain.

Glassy Polymers—Diffusion, Sorption, Ageing and Applications

For the past few decades, researchers have been intrigued by glassy polymers, which have applications ranging from gas separations to corrosion protection to drug delivery systems. The techniques employed to examine the sorption and diffusion of small molecules in glassy polymers are the subject of this review. Diffusion models in glassy polymers are regulated by Fickian and non-Fickian diffusion, with non-Fickian diffusion being more prevalent. The characteristics of glassy polymers are determined by sorption isotherms, and different models have been proposed in the literature to explain sorption in glassy polymers over the last few years. This review also includes the applications of glassy polymers. Despite having many applications, current researchers still have difficulty in implementing coating challenges due to issues such as physical ageing, brittleness, etc., which are briefly discussed in the review.

A Molecular Simulation Study of Silica/Polysulfone Mixed Matrix Membrane for Mixed Gas Separation

Polysulfone-based mixed matrix membranes (MMMs) incorporated with silica nanoparticles are a new generation material under ongoing research and development for gas separation. However, the attributes of a better-performing MMM cannot be precisely studied under experimental conditions. Thus, it requires an atomistic scale study to elucidate the separation performance of silica/polysulfone MMMs. As most of the research work and empirical models for gas transport properties have been limited to pure gas, a computational framework for molecular simulation is required to study the mixed gas transport properties in silica/polysulfone MMMs to reflect real membrane separation. In this work, Monte Carlo (MC) and molecular dynamics (MD) simulations were employed to study the solubility and diffusivity of CO₂/CH₄ with varying gas concentrations (i.e., 30% CO₂/CH₄, 50% CO₂/CH₄, and 70% CO₂/CH₄) and silica content (i.e., 15–30 wt.%). The accuracy of the simulated structures was validated with published literature, followed by the study of the gas transport properties at 308.15 K and 1 atm. Simulation results concluded an increase in the free volume with an increasing weight percentage of silica. It was also found that pure gas consistently exhibited higher gas transport properties when compared to mixed gas conditions. The results also showed a competitive gas transport performance for mixed gases, which is more apparent when CO₂ increases. In this context, an increment in the permeation was observed for mixed gas with increasing gas concentrations (i.e., 70% CO₂/CH₄ > 50% CO₂/CH₄ > 30% CO₂/CH₄). The diffusivity, solubility, and permeability of the mixed gases were consistently increasing until 25 wt.%, followed by a decrease for 30 wt.% of silica. An empirical model based on a parallel resistance approach was developed by incorporating mathematical formulations for solubility and permeability. The model results were compared with simulation results to quantify the effect of mixed gas transport, which showed an 18% and 15% percentage error for the permeability and solubility, respectively, in comparison to the simulation data. This study provides a basis for future understanding of MMMs using molecular simulations and modeling techniques for mixed gas conditions that demonstrate real membrane separation.

Optical Analysis of the Internal Void Structure in Polymer Membranes for Gas Separation

Global warming by greenhouse gas emissions is one of the main threats of our modern society, and efficient CO2 capture processes are needed to solve this problem. Membrane separation processes have been identified among the most promising technologies for CO2 capture, and these require the development of highly efficient membrane materials which, in turn, requires detailed understanding of their operation mechanism. In the last decades, molecular modeling studies have become an extremely powerful tool to understand and anticipate the gas transport properties of polymeric membranes. This work presents a study on the correlation of the structural features of different membrane materials, analyzed by means of molecular dynamics simulation, and their gas diffusivity/selectivity. We propose a simplified method to determine the void size distribution via an automatic image recognition tool, along with a consolidated Connolly probe sensing of space, without the need of demanding computational procedures. Based on a picture of the void shape and width, automatic image recognition tests the dimensions of the void elements, reducing them to ellipses. Comparison of the minor axis of the obtained ellipses with the diameters of the gases yields a qualitative estimation of non-accessible paths in the geometrical arrangement of polymeric chains. A second tool, the Connolly probe sensing of space, gives more details on the complexity of voids. The combination of the two proposed tools can be used for a qualitative and rapid screening of material models and for an estimation of the trend in their diffusivity selectivity. The main differences in the structural features of three different classes of polymers are investigated in this work (glassy polymers, superglassy perfluoropolymers and high free volume polymers of intrinsic microporosity), and the results show how the proposed computationally less demanding analysis can be linked with their selectivities.

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.

Transport Properties of Thermoplastic R-BAPB Polyimide: Molecular Dynamics Simulations and Experiment

The present work evaluates the transport properties of thermoplastic R-BAPB polyimide based on 1,3-bis(3,3',4,4'-dicarboxyphenoxy)benzene (dianhydride R) and 4,4'-bis(4-aminophenoxy)biphenyl (diamine BAPB). Both experimental studies and molecular dynamics simulations were applied to estimate the diffusion coefficients and solubilities of various gases, such as helium (He), oxygen (O2), nitrogen (N2), and methane (CH4). The validity of the results obtained was confirmed by studying the correlation of the experimental solubilities and diffusion coefficients of He, O2, and N2 in R-BAPB, with their critical temperatures and the effective sizes of the gas molecules, respectively. The solubilities obtained in the molecular dynamics simulations are in good quantitative agreement with the experimental data. A good qualitative relationship between the simulation results and the experimental data is also observed when comparing the diffusion coefficients of the gases. Analysis of the Robeson plots shows that R-BAPB has high selectivity for He, N2, and CO2 separation from CH4, which makes it a promising polymer for developing gas-separation membranes. From this point of view, the simulation models developed and validated in the present work may be put to effective use for further investigations into the transport properties of R-BAPB polyimide and nanocomposites based on it.