Molecular Dynamics Study of Carbon Dioxide Sorption and Plasticization at the Interface of a Glassy Polymer Membrane

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.

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

With a wide range of applications, from energy and environmental engineering, such as in gas separations and water purification, to biomedical engineering and packaging, glassy polymeric materials remain in the core of novel membrane and state-of the art barrier technologies. This review focuses on molecular simulation methodologies implemented for the study of sorption and diffusion of small molecules in dense glassy polymeric systems. Basic concepts are introduced and systematic methods for the generation of realistic polymer configurations are briefly presented. Challenges related to the long length and time scale phenomena that govern the permeation process in the glassy polymer matrix are described and molecular simulation approaches developed to address the multiscale problem at hand are discussed.

Relaxation Dynamics of Thin Matrimid 5218 Films in Organic Solvents

Polyimides are interesting polymer materials for organic solvent nanofiltration (OSN) applications because of their high excess free volume and high chemical and temperature resistance. However, an open challenge that remains for glassy polymer materials (i.e., polyimides) is their tendency to swell in organic solvents which can lead to a loss of performance. An understanding on how swelling influences the polymer properties and performance is then of crucial importance for assessing polyimide suitability in OSN applications. Here, the combination of in situ spectroscopic ellipsometry (iSE), broadband dielectric spectroscopy (BDS), and diffuse reflectance Fourier transform infrared spectroscopy (DRIFT-FTIR) is applied to study the molecular interaction of two organic penetrants, toluene and n-hexane, with Matrimid 5218 in detail. iSE shows that slightly cross-linked Matrimid 5218 swells approximately seven times more in toluene (swelling degree ≈ 28%) compared to in n-hexane (swelling degree ≈ 4%). Combined BDS and DRIFT-FTIR results indicate that toluene interacts with the benzene ring present in the diamine via π–π interactions, while n-hexane likely fills up the excess free volume and interacts via local van der Waals interactions. This work highlights the insights into the exact nature of the molecular interactions between the penetrant and polymer that can be gained from a combination of BDS and other techniques and how these insights can be used to estimate or understand solvent-induced swelling of polymers used in OSN applications.

Evolution of Filament-Shaped Porous Structure in Polycarbonate by Stretching under Carbon Dioxide

We found that a filament-shaped porous structure with periodic distance was obtained in polycarbonate for optical disk grade (OD-PC) film by stretching under compressed carbon dioxide (CO₂). The evolution of the characteristic porous structure was investigated by in situ observation during the stretching under compressed CO₂ and the optical microscopic observation of the stretched specimen. The voids were obtained under high CO₂ pressure as in the case of elevated temperature, suggesting that the evolution of the voids was caused by crazing due to chain disentanglement by accelerated molecular motion owing to the plasticization effect of CO₂. The filament-shaped voids were initiated at around the yielding point and increased continuously by nucleation in the matrix around the surface of the pre-existing voids. The shape of the voids did not change to an ellipsoidal one during stretching due to suppression of the craze opening by the hydrostatic pressure effect. The stretching of the CO₂-absorbed depressurized OD-PC revealed that the initiation of the voids was not only caused by the plasticization effect, but the hydrostatic pressure effect was also required.

A Review of Multiscale Computational Methods in Polymeric Materials

Polymeric materials display distinguished characteristics which stem from the interplay of phenomena at various length and time scales. Further development of polymer systems critically relies on a comprehensive understanding of the fundamentals of their hierarchical structure and behaviors. As such, the inherent multiscale nature of polymer systems is only reflected by a multiscale analysis which accounts for all important mechanisms. Since multiscale modelling is a rapidly growing multidisciplinary field, the emerging possibilities and challenges can be of a truly diverse nature. The present review attempts to provide a rather comprehensive overview of the recent developments in the field of multiscale modelling and simulation of polymeric materials. In order to understand the characteristics of the building blocks of multiscale methods, first a brief review of some significant computational methods at individual length and time scales is provided. These methods cover quantum mechanical scale, atomistic domain (Monte Carlo and molecular dynamics), mesoscopic scale (Brownian dynamics, dissipative particle dynamics, and lattice Boltzmann method), and finally macroscopic realm (finite element and volume methods). Afterwards, different prescriptions to envelope these methods in a multiscale strategy are discussed in details. Sequential, concurrent, and adaptive resolution schemes are presented along with the latest updates and ongoing challenges in research. In sequential methods, various systematic coarse-graining and backmapping approaches are addressed. For the concurrent strategy, we aimed to introduce the fundamentals and significant methods including the handshaking concept, energy-based, and force-based coupling approaches. Although such methods are very popular in metals and carbon nanomaterials, their use in polymeric materials is still limited. We have illustrated their applications in polymer science by several examples hoping for raising attention towards the existing possibilities. The relatively new adaptive resolution schemes are then covered including their advantages and shortcomings. Finally, some novel ideas in order to extend the reaches of atomistic techniques are reviewed. We conclude the review by outlining the existing challenges and possibilities for future research.

The influence of the dianhydride precursor in hyper-cross-linked hybrid polyPOSS-imide networks

Hybrid organic/inorganic hyper-cross-linked membranes based on imides covalently bonded with polyhedral oligomeric silsesquioxanes (POSS) have recently been developed for gas-separation applications under high pressure and/or temperature conditions. Their molecular sieving capabilities have been shown to depend on the nature of the organic dianhydride precursor. In the present work, realistic molecular models of such polyPOSS-imide films based on the flexible 6FDA dianhydride are compared to those based on the shorter and more rigid PMDA dianhydride. The models creation procedure closely mimicks the mixing, polycondensation and imidization steps of the experimental scheme. The resulting networks are found to be highly heterogeneous in terms of both the number of links (from zero to the maximum possible of eight per POSS cage with an average of four) and their structure (interPOSS, intraPOSS, single-links, double-links) because of the eight-equivalent-arms nature of the POSS precursor. For both dianhydride precursors, crosslinking with POSS and the subsequent imidization step decrease the density, create additional void-space and increase the solubility of the resulting membranes. However, when compared to PMDA, the added flexibility of the central 6FDA bridge leads to a larger thermally-induced dilation of the networks and a larger volume loss per H₂O over the imidization step. With their better ability to redensify and to adapt to the added constraints, the cagecage distances and cage(organic bridge)cage angles in the 6FDA polyPOSS-imides span a larger range than in their PMDA counterparts. In addition, the stiffness of the PMDA moiety results in more unrelaxed free volume remaining trapped in the PMDA polyPOSS-imides upon imidization, and as such, to significantly more open structures with less favourable interactions. As expected from their enhanced flexibility, the thermomechanical properties of the 6FDA networks are slightly lower than those based on PMDA. However, the better mechanical resistance of PMDA over 6FDA does not really become significant before very large volume dilations.

In Silico Determination of Gas Permeabilities by Non-Equilibrium Molecular Dynamics: CO2 and He through PIM-1

We study the permeation dynamics of helium and carbon dioxide through an atomistically detailed model of a polymer of intrinsic microporosity, PIM-1, via non-equilibrium molecular dynamics (NEMD) simulations. This work presents the first explicit molecular modeling of gas permeation through a high free-volume polymer sample, and it demonstrates how permeability and solubility can be obtained coherently from a single simulation. Solubilities in particular can be obtained to a very high degree of confidence and within experimental inaccuracies. Furthermore, the simulations make it possible to obtain very specific information on the diffusion dynamics of penetrant molecules and yield detailed maps of gas occupancy, which are akin to a digital tomographic scan of the polymer network. In addition to determining permeability and solubility directly from NEMD simulations, the results shed light on the permeation mechanism of the penetrant gases, suggesting that the relative openness of the microporous topology promotes the anomalous diffusion of penetrant gases, which entails a deviation from the pore hopping mechanism usually observed in gas diffusion in polymers.

How does low-molecular-weight polystyrene dissolve: osmotic swelling vs. surface dissolution

By means of multiscale hierarchical modeling we study the real time evolution of low-molecular-weight polystyrene, below the glass transition temperature, in contact with its solvent, toluene. We observe two concurrent phenomena taking place: (1) the solvent diffuses into the polymer by a Case II mechanism, leading to osmotic driven swelling and progressive chain dilution (inside-out mechanism); (2) polymer chains are solvated, detach from the interface and move into the solvent before the film is completely swollen (outside-in mechanism). From our simulations we conclude that, below the entanglement length, a thin swollen layer, also observed in previous experiments, forms almost instantaneously, which allows for the outside-in mechanism to start a few tens of nanoseconds after the polymer-solvent initial contact. After this initial transient time the two mechanisms are concurrent. We furthermore observe that the presence of the solvent significantly enhances the mobility of the polymer chains of the surface layer, but only in the direction parallel to the interface.

A comparison of homopolymer and block copolymer structure in 6FDA-based polyimides

Two homopolyimides and the corresponding block copolyimide, all based on the 4,4'-(hexafluoroisopropylidene)diphthalic dianhydride (6FDA), have been synthesized and fully atomistic models have been studied using molecular dynamics (MD) simulation. The respective diamines were 1,3-phenylenediamine (mPDA) and 2,3,5,6-tetramethyl-1,4-phenylenediamine (durene). These polyimides are potential candidates for gas separation applications. The synthesized polymers were processed as dense flat membranes. The effects of diamine structure were investigated at the molecular level and an attempt to compare the structural features of homo- and block copolyimides was made. Amorphous models were generated using a hybrid pivot Monte Carlo-MD sampling preparation technique. Average model densities were validated against experimental measurements on the dense films. Cohesive energies, Hildebrand solubility parameters, conformational characteristics, intermolecular interactions and available void spaces were analysed for each system. The durene diamine was found to hinder stacking and increase the available space. This is associated with the steric effect of the methyl substituents. In general, 6FDA-mPDA/durene exhibits an intermediate behaviour with respect to its base polyimides. For most of the examined properties, the differences between different size simulated systems were minor with the exception of the free volume distribution.