Molecular Dynamics Simulations of Oxygen Transport through a Fully Atomistic Polyimide Membrane

Multiscale equilibration of highly entangled isotropic model polymer melts

We present a computationally efficient multiscale method for preparing equilibrated, isotropic long-chain model polymer melts. As an application, we generate Kremer-Grest melts of 1000 chains with 200 entanglements and 25 000-2000 beads/chain, which cover the experimentally relevant bending rigidities up to and beyond the limit of the isotropic-nematic transition. In the first step, we employ Monte Carlo simulations of a lattice model to equilibrate the large-scale chain structure above the tube scale while ensuring a spatially homogeneous density distribution. We then use theoretical insight from a constrained mode tube model to introduce the bead degrees of freedom together with random walk conformational statistics all the way down to the Kuhn scale of the chains. This is followed by a sequence of simulations with carefully parameterized force-capped bead-spring models, which slowly introduce the local bead packing while reproducing the larger-scale chain statistics of the target Kremer-Grest system at all levels of force-capping. Finally, we can switch to the full Kremer-Grest model without perturbing the structure. The resulting chain statistics is in excellent agreement with literature results on all length scales accessible in brute-force simulations of shorter chains.

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.

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.

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.

Thermal properties of bulk polyimides: insights from computer modeling versus experiment

Due to the great importance for many industrial applications it is crucial from the point of view of theoretical description to reproduce thermal properties of thermoplastic polyimides as accurate as possible in order to establish "chemical structure-physical properties" relationships of new materials. In this paper we employ differential scanning calorimetry, dilatometry, and atomistic molecular dynamics (MD) simulations to explore whether the state-of-the-art computer modeling can serve as a precise tool for probing thermal properties of polyimides with highly polar groups. For this purpose the polyimide R-BAPS based on dianhydride 1,3-bis(3',4-dicarboxyphenoxy)benzene (dianhydride R) and diamine 4,4'-bis(4''-aminophenoxy)biphenyl sulphone) (diamine BAPS) was synthesized and extensively studied. Overall, our findings show that the widely used glass-transition temperature Tg evaluated from MD simulations should be employed with great caution for verification of the polyimide computational models against experimental data: in addition to the well-known impact of the cooling rate on the glass-transition temperature, correct definition of Tg requires cooling that starts from very high temperatures (no less than 800 K for considered polyimides) and accurate evaluation of the appropriate cooling rate, otherwise the errors in the measured values of Tg become undefined. In contrast to the glass-transition temperature, the volumetric coefficient of thermal expansion (CTE) does not depend on the cooling rate in the low-temperature domain (T < Tg) so that comparison of computational and experimental values of CTE provides a much safer way for proper validation of the theoretical model when electrostatic interactions are taken into account explicitly. Remarkably, this conclusion is most likely of generic nature: we show that it also holds for the commercial polyimide EXTEM, another polyimide with a similar chemical structure.