Molecular modeling of poly(benzoxazole-co-imide) membranes: A structure characterization and performance investigation

Molecular Dynamics Simulations of Gas Transport Properties in Cross-Linked Polyamide Membranes: Tracing the Morphology and Addition of Silicate Nanotubes

This study employs molecular dynamics (MD) simulations to fundamentally provide insight into the role of cross-link density in the CO2 separation properties of interfacially polymerized polyamide (PA) membranes. For this purpose, two atomistic models of pure polyamide membranes with different cross-link densities are constructed by MD simulations to conceptually determine how the fractional free volume of polyamide affects the gas separation performance of the membrane. The PA membrane with a lower cross-link density (LCPA) shows a higher gas diffusion coefficient, a lower gas solubility coefficient, and a higher gas permeability than the PA membrane with a higher cross-link density (HCPA). Moreover, the pristine and modified silicate nanotubes (SNTs) as the fast gas transport channels are incorporated into the polyamide membranes to assess the effect of the SNT/PA interface chemistry on the CO2 separation properties of the membranes. SNTs are systematically modified by three modifying agents with different CO2-philic groups and different interfacial interaction energies with the polyamide matrix. The results of MD simulations demonstrate that the incorporation of silicate nanotubes into the PA matrix increases the gas diffusivity and permeability and decreases the CO2/gas selectivity. Moreover, the membranes containing modified SNTs possessing high CO2-philicity and high SNTs/PA interfacial interactions show a high CO2 separation performance.

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.

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.

High-Drug-Loading Amorphous Solid Dispersions via In Situ Thermal Cross-Linking: Unraveling the Mechanisms of Stabilization

This article takes a step forward in understanding the mechanisms involved during the preparation and performance of cross-linked high-drug-loading (HDL) amorphous solid dispersions (ASDs). Specifically, ASDs, having 90 wt % poorly water-soluble drug indomethacin (IND), were prepared via in situ thermal cross-linking of poly(acrylic acid) (PAA) and poly(vinyl alcohol) (PVA) and thoroughly evaluated in terms of physical stability and in vitro supersaturation. Results showed that HDL ASDs having excellent active pharmaceutical ingredient (API) amorphous stability and prolonged in vitro supersaturation were prepared by fine tuning the cross-linking procedure. Unraveling of the processes involved during ASD's formation shed light on the significant role of the cross-linking conditions (i.e., temperature and time), the physicochemical properties of the API, and the hydrolysis level of the cross-linker as key factors in modulating ASD's stability. In-depth analysis of the prepared systems revealed the (1) reduction of API's molecular motions within the cross-linked polymeric networks (through API's strong spatial confinement), (2) the structural changes in the prepared cross-linked matrices (induced by the high API drug loading), and (3) the tuning of the cross-linking density via utilization of low-hydrolyzed PVA as the major mechanisms responsible for ASD's exceptional performance. Complementary analysis by means of molecular dynamics simulations also highlighted the vital role of strong drug-polymer intermolecular interactions evolving among the ASD components. Overall, the impression of the complexity of in situ cross-linked ASDs has been reinforced with the excessive variation of parameters investigated in the current study, offering thus insights up to the submolecular level to lay the groundwork and foundations for the comprehensive assessment of a new emerging class of HDL amorphous API formulations.

Revealing the High-Modulus Mechanism of Polyimide Films Prepared with 3,4'-ODA

To prepare PIs (polyimides) with desirable thermal and mechanical properties is highly demanded due to their widespread applications in flexible optoelectronic devices and printed circuit boards. Here, the PI films of BPDA/4,4′-ODA, BPDA/3,4′-ODA, PMDA/4,4′-ODA, PMDA/3,4′-ODA systems were prepared, and it was found that the PIs with 3,4′-ODA always exhibit a high modulus compared with the PIs with 4,4′-ODA. To disclose the mechanism of high-modulus PI films with 3,4′-ODA, amorphous PI models and uniaxial drawing PI models were established and calculated based on MD simulation. The PI structural deformations at different length scales, i.e., molecular chain cluster scale and repeat unit scale, under the same stress were detailed and analyzed, including the variation of chain conformation, bond length, bond angle, internal rotation energy, and torsion angle. The results indicate that PIs with 3,4-ODA have higher internal rotation energy and smaller deformation with the same stress, consistent with the high modulus.

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.

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.

Molecular Simulation Study on Adsorption and Diffusion Behaviors of CO2/N2 in Lignite

Understanding the adsorption and diffusion of CO₂ and N₂ in lignite at high temperature is of great significance for fire prevention and control. Considering the influence of temperature on coal structure, molecular structure models of lignite at 298.15, 323.15, 423.15, and 523.15 K were constructed by molecular mechanics and dynamics, and grand canonical Monte Carlo molecular simulation was conducted for single-component and two-component systems under different temperatures, pressures, and gas ratios. The adsorption capacity was positively correlated with the pressure and molar ratio of CO₂ but negatively correlated with the temperature. The adsorption amount of CO₂ (1.060 mmol/g) was generally larger than that of N₂ (0.069 mmol/g), showing a greater selectivity. However, CO₂ was more sensitive to temperature, and the adsorption amount decreased faster with the increase in temperature. At high temperature, the adsorption amount of CO₂ and N₂ is basically equal, both of which are at a low level. The CO₂ isosteric heat of adsorption (7.46-8.84 kcal/mol) varies significantly with temperature. The interaction energy is consistent with the change trend of adsorption quantity, and van der Waals energy plays a dominant role in adsorption. Injecting CO₂ and N₂ at the high temperature stage has a poor extinguishing effect, which can only dilute oxygen content and exchange heat, and the advantage of CO₂ will be lost. It may be more effective using liquid nitrogen which has the properties of low temperature and high specific heat capacity. The results are of great significance to improve the efficiency of fire prevention and suppression in underground coal mines.

Effect of introduction of fluoromonomer copolymerization on properties of polyimide hollow fibers

The novel copolymerization polyimide (PI) hollow fibers (HFs) of 4,4′-(hexafluoroisopropylidene) diphthalic anhydride (6FDA) containing –CF3 groups were prepared and investigated through both simulation and experiment. To demonstrate the alteration attributable to the introduction of fluoromonomers, the condensed states of pyromellitic dianhydride (PMDA)/4,4′-oxybisbenzenamine (ODA), PMDA/6FDA/ODA, and 6FDA/ODA PI were constructed by Material Studio, and we simulated the mobility of molecular chain, free volume fraction, and O2/N2 dissolution–diffusion process. The molecular dynamics simulation results demonstrated that the properties of the copolymerized PI system with 6FDA were significantly improved, while the selectivity remained almost unchanged. Then, the films of copolymerized PI and HFs were prepared by the two-step method, and O2/N2 permeability of the PI copolymer films was characterized, indicating that although the gas permeation performance was greatly improved, the selectivity was not so satisfactory. However, the selection factor increased heavily after polydimethylsiloxane coating.

Molecular Model Construction of the Dense Medium Component Scaffold in Coal for Molecular Aggregate Simulation

Coal as an important fossil energy has been comprehensively studied in terms of its structure, reactivity, and application. However, there are few publications reported about the formation mechanism of coal. In order to explore the molecular mechanism of the formation of the dense medium component (DMC) aggregate, which is extracted from coal, the molecular model of the DMC scaffold (DMC-S) was constructed based on a number of X-ray photoelectron spectroscopy, ¹³C NMR, and ultimate analysis. Then, DMC-S was further optimized, and the periodic boundary condition was added for molecular mechanics and molecular dynamics simulation. The DMC-S molecule model with a density of 1.05 g/cm³ and a different number of unit cells was obtained after the aforementioned experiments and simulations. When the unit cell contained 12 DMC-S molecules, the absolute value of electrostatic energy significantly increased and the peripheral branch chains in DMC-S interlaced with each other, forming a compact aggregate. The density and macrosize calculated values are all slightly lower than the true relative values because the presence of minerals or small molecules was not included in the model construction. Despite some unavoidable defects, the comparison between the simulated and experimental results validates the DMC-S aggregate model and lays a solid foundation for an in-depth study of DMC and its reactivity.

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.

Polyimides Containing Phosphaphenanthrene Skeleton: Gas-Transport Properties and Molecular Dynamics Simulations

A series of new semifluorinated polyimide (PI) films with phosphaphenanthrene skeleton were prepared by thermal imidization of poly(amic acid)s derived from a diamine monomer: 1,1-bis[2'-trifluoromethyl-4'-(4″-aminophenyl)phenoxy]-1-(6-oxido-6H-dibenz⟨c,e⟩⟨1,2⟩oxaphosphorin-6-yl)ethane on reaction with four structurally different aromatic dianhydrides. The chemical structures of the polymers were established by Fourier transform infrared and 1H NMR spectroscopy techniques. The polymers showed a good combination of thermal and mechanical properties (T d10 up to 416 °C under synthetic air and tensile strength up to 91 MPa), low dielectric constant (2.10-2.55 at 1 MHz), and T g values as high as 261 °C. Gas permeabilities of these films were investigated for four different gases CO2, O2, N2, and CH4. The PI films showed high gas permeability (P CO2 up to 175 and P O2 up to 64 barrer) with high permselectivity (P CO2 /P CH4 up to 51 and P O2 /P N2 up to 7.1), and the values are better than those of many other similar polymers reported earlier. For the O2/N2 gas pair, the PIs (PI A) surpassed the present upper boundary limit drawn by Robeson. A detailed molecular dynamics (MD) simulation study has been conducted to understand better the gas-transport properties. The effect of phosphaphenanthrene skeleton, its spatial arrangement, and size distribution function of the free volume were studied using molecular dynamics (MD) simulation and the results are correlated with the experimental data.

Soluble, optically transparent polyamides with a phosphaphenanthrene skeleton: synthesis, characterization, gas permeation and molecular dynamics simulations

Soluble, optically transparent polyamides with a phosphaphenanthrene skeleton: synthesis, characterization, gas permeation and molecular dynamics simulations.