Next Article in Journal
Bacterial Nanocellulose Hydrogel: A Promising Alternative Material for the Fabrication of Engineered Vascular Grafts
Next Article in Special Issue
Self-Assembly of Symmetric Copolymers in Slits with Inert and Attractive Walls
Previous Article in Journal
Recyclable Superhydrophobic Surface Prepared via Electrospinning and Electrospraying Using Waste Polyethylene Terephthalate for Self-Cleaning Applications
Previous Article in Special Issue
Molecular Weight Segregation and Thermal Conductivity of Polydisperse Wax–Graphene Nanocomposites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Molecular Dynamics Study of Single-Gas and Mixed-Gas N2 and CH4 Transport in Triptycene-Based Polyimide Membranes

1
Univ. Savoie Mont Blanc, Univ. Grenoble Alpes, CNRS, Grenoble INP, LEPMI, 38000 Grenoble, France
2
Saudi Aramco, Research & Development Center, P.O. Box 62, Dhahran 31311, Saudi Arabia
*
Author to whom correspondence should be addressed.
Current address: Siemens Digital Industries Software, 2595 BN The Hague, The Netherlands.
Polymers 2023, 15(18), 3811; https://doi.org/10.3390/polym15183811
Submission received: 27 July 2023 / Revised: 4 September 2023 / Accepted: 13 September 2023 / Published: 18 September 2023

Abstract

:
Fluorinated polyimides incorporated with triptycene units have gained growing attention over the last decade since they present potentially interesting selectivities and a higher free volume with respect to their triptycene-free counterparts. This work examines the transport of single-gas and mixed-gas N2 and CH4 in the triptycene-based 6FDA-BAPT homopolyimide and in a block 15,000 g mol−1/15,000 g mol−1 6FDA-mPDA/BAPT copolyimide by using molecular dynamics (MD) simulations. The void-space analyses reveal that, while the free volume consists of small-to-medium holes in the 6FDA-BAPT homopolyimide, there are more medium-to-large holes in the 6FDA-mPDA/BAPT copolyimide. The single-gas sorption isotherms for N2 and CH4 over the 0–70 bar range at 338.5 K show that both gases are more soluble in the block copolyimide, with a higher affinity for methane. CH4 favours sites with the most favourable energetic interactions, while N2 probes more sites in the matrices. The volume swellings remain limited since neither N2 nor CH4 plasticise penetrants. The transport of a binary-gas 2:1 CH4/N2 mixture is also examined in both polyimides under operating conditions similar to those used in current natural gas processing, i.e., at 65.5 bar and 338.5 K. In the mixed-gas simulations, the solubility selectivities in favour of CH4 are enhanced similarly in both matrices. Although diffusion is higher in 6FDA-BAPT/6FDA-mPDA, the diffusion selectivities are also close. Both triptycene-based polyimides under study favour, to a similar extent, the transport of methane over that of nitrogen under the conditions studied.

1. Introduction

The separation of small gas molecules is of major importance for the recovery of valuable gases and fuel production [1,2,3], as well as in other industrial processes, such as air separation, food packaging, and pollution control [2,4,5,6]. Polymeric-membrane-mediated gas separation appears as one of the preferred choices for the gas separation industry over traditional processes, such as cryogenic distillation, filtration, ion exchange, and thermal and pressure-swing adsorption, due to its simpler fabrication, lower costs, energy efficiency, and reduced environmental impact [7,8,9,10,11,12,13,14,15]. Recent studies report a class of hydrocarbon polymers that can achieve both high selectivity and high gas permeability in membrane separations for many gas mixtures [16], whereas other current works highlight the contribution of machine learning models to the optimisation of membrane fabrication and design [17,18,19,20,21,22,23].
Applications of polymer membranes in natural gas processing often involve denitrogenation, i.e., the removal of excess nitrogen from natural gas reservoirs, as the former lowers the thermal efficiency (heating value) of natural gas and increases transport volumes. In addition, most pipeline standards require that natural gas contains less than 4% N2 [24,25,26]. In general, nitrogen is difficult to separate from methane due to their similar kinetic diameters (3.64 and 3.8 Å, respectively) [27] and physical properties. However, membrane-mediated separation is accomplished via differences in both gas solubilities and diffusivities [28,29,30].
The transport of gases through dense membranes is a solution–diffusion process, where gas molecules are initially sorbed into the membrane surface, diffuse across the membrane via a series of jumps through inter-connected cavities, and, finally, desorb at the downstream compartment of the membrane [31,32]. The ability of a membrane to separate a given gas pair is usually assessed by its permeability and its selectivity or separation factor. In the case where the downstream pressure is negligible compared to the upstream pressure of the gas feed, the permeability, P, is usually expressed as the product of the solubility and diffusion coefficients (Equation (1)):
P i = S c i D i
The diffusion coefficient, D i , is the kinetic contribution to the permeability, whereas the solubility coefficient, S c i , characterises the partitioning of the gas penetrants into the polymer and is given by the ratio of the concentration of the gas in the polymer, C i , to the partial pressure, p i , of the gas to which it is exposed (Equation (2)):
S c i = C i p i
Assuming that a gas, A, is more permeable than a gas, B, the separation capability of a membrane for a given gas pair, αA/B, is given by the ratio of the permeabilities (Equation (3)):
α A / B = P A P B = D A D B S A S B
Amongst the various polymer families used in gas separations, aromatic polyimides (PIs) have attracted special attention due to their satisfactory chemical and thermal stabilities, as well as their resistance to mechanical stresses [7,33,34,35]. Moreover, hexafluorodianhydride (6FDA)-based PIs appear as particularly appealing candidates in separation processes as they provide considerably higher selectivities than other dianhydrides [36,37,38].
Despite the exceptional properties of aromatic PIs, there still exist several challenges for long-term industrial applications, which are primarily related to the permeability/selectivity tradeoff, physical ageing, and membrane plasticisation in the presence of highly condensable gases [39]. These issues can partly be circumvented at the macromolecular level by integrating bulky moieties and/or rigid tortuous elements in the polymer backbone in order to inhibit polymer chain packing and increase the fractional free volume (FFV) of the material [39,40].
Recent studies have highlighted the usage of iptycene molecules as building elements to address some of the challenges mentioned earlier [7,41,42,43,44]. Triptycene is a rigid three-dimensional molecular unit with three blades, each composed of a benzene ring [45,46]. Several groups have synthesised novel triptycene-based PIs and examined structure–property relationships that can guide the rational design of improved-performance membranes targeted for gas separations. One of the first-synthesised triptycene-based PIs bears the commercial name 6FDA-DATRI [47]. This PI was found to have a higher FFV than other 6FDA-based PIs and a CO2/CH4 selectivity that outperforms any commercial PI. In addition, 6FDA-DATRI presents a satisfactory resistance to plasticisation, which was attributed to the π–π interactions of the triptycene benzene blades perpendicular to the polymer backbone [47]. Wiegand et al. synthesised a 6FDA-1,4-triptycene polyimide, also known as 6FDA-BAPT, and compared its gas separation performance with other polyimides [40]. They measured an ideal selectivity of 1.78 for the N2/CH4 pair at T = 35° and p = 9 bar and observed a favourable free-volume size distribution induced by the symmetric nature of the triptycene unit [40]. Weidman and co-workers examined the physical ageing properties of 6FDA-BABT and found that the addition of substituent groups resulted in an improved permeability over time due to steric hindrance effects [48]. Other groups examined the permselectivity properties of thermally rearranged triptycene-containing polybenzoxazoles (PBO) derived from polyimide precursors with triptycene units [49,50]. Specifically, triptycene-containing TR-PBO membranes synthesised by Luo et al. were found to have an enhanced N2/CH4 ideal selectivity greater than 2 at p = 11 atm and T = 35 °C [49].
As a complement to experiments, molecular dynamics (MD) simulations are a useful tool to study structure–property relationships and the gas separation performance of triptycene-based PIs. In particular, Cho et al. [47] examined the free-volume distribution and gas transport in 6FDA-DATRI membranes in a combined experimental and simulation study. Their calculations yielded an FFV higher than other polyimide membranes, whereas a selectivity of 1.3 for N2/CH4 at p = 1 bar and T = 35 °C was reported [47]. A combination of MD, Monte Carlo (MC), and ab initio simulations was used by Chen et al. [51] to examine the structural properties and transport behaviour of a triptycene-based polymer with intrinsic microporosity (PIM) membranes. They showed that the incorporation of the triptycene moiety rigid structure disrupts the polymer chain packing, leading to enhanced diffusion selectivities for the O2/N2 pair [51]. In an analogous study of Ghasemnejad-Afshar et al., it was shown that PIM membranes with larger side branch groups in their polymeric chain structure are more rigid, which is due to restriction in chain packing and cavity formation between polymer chains [52]. More recently, Balcik et al. [53] investigated the pure- and mixed-gas CO2/CH4 separation properties of a triptycene–Tröger’s base ladder PIM-Trip-TB polymer. They found that the CO2/CH4 solubility selectivity is enhanced under mixed-gas conditions, but despite the intra-chain rigidity, an increase in FFV upon exposure to plasticising CO2 leads to a loss in diffusion solubility [53].
The aim of this contribution is to examine the transport of the less plasticising N2 and CH4 in triptycene-based homo- and block copolyimides using MD simulations. In particular, we assess the free-volume distribution and transport behaviour of 6FDA-1,4-triptycene, i.e., 6FDA-BAPT, for both single-gas and mixed-gas feeds. Given its promising transport properties [40], it appears interesting to combine it with the 6FDA-mPDA polyimide, which exhibits better packing on its own and is selective over nitrogen [54,55]. To this end, we also examine the performance of a block copolyimide made of the aforementioned base homopolyimides. The chemical formulae of the polyimides under study are shown in Figure 1.
As shown in previous work, an estimation of the permselectivity behaviour of the polyimides exclusively from single-gas conditions can lead to misleading findings [56,57,58,59]. Therefore, permeation studies were carried out for pure N2 and CH4 loadings and also for binary mixtures at a ratio equivalent to those found in some natural gas reservoirs, i.e., as used before [59], a 2:1 mixture of CH4/N2 at 65.5 bar and at 338.5 K. These are the processing conditions that were specified by the industrial partner. The preparation of the amorphous polyimide models, along with the examination of their void space, is described in Section 2. Section 3 presents the simulation protocols and the methods to measure model gas permeation in the polyimide matrices for both pure- and mixed-gas feeds at 338.5 K. The compared single-gas and binary-gas analyses of N2 and CH4 transport in the 6FDA/BAPT polyimide and its 6FDA-mPDA/BAPT copolyimide are presented in Section 4.

2. Modelling the Pure Triptycene-Based 6FDA-Polyimides

2.1. Preparation of the Amorphous Polyimide and Copolyimide Models

All simulations were carried out with the aid of the scalar and parallel versions of the general-purpose gmq software version VI [60]. As the transport characteristics of a membrane are governed by the core rather than by the membrane–gas interface, the triptycene-based polyimide representations were bulk models in three dimensions with periodic boundary conditions [61]. The fully atomistic amorphous models of 6FDA-BAPT and 6FDA-mPDA/BAPT contained each 6 chains and were prepared via the hybrid pivot Monte Carlo–molecular dynamics (PMC-MD) method [62,63]. This allows for the generation of a polymer sample at the desired density with the polymer chains bearing conformations that correspond to the equilibrium melt at the required temperature. The underlying theory stems from Flory’s hypothesis [64], according to which the configurations of polymer chains in the pure melt are determined by only a certain number of near-neighbour intramolecular interactions. The method has been thoroughly described elsewhere and validated against various polymer families [62,63]. Earlier studies on polyimide systems reported that the chain characteristic ratio, R 2 / n m o n o m e r s , i.e., the ratio of the mean square end-to-end distance of the chains <R2> over the number of monomers, becomes number-independent from nmonomers ~30 [38,39]. On these grounds, the mean molecular weight of the 6FDA-BAPT chains was ~40,000 g mol−1, with each chain comprising l = 45 monomers (Figure 1), which is consistent with experimental studies [40]. Motivated by a previous study demonstrating that block copolyimides of a similar block size show promising gas transport properties [54], we considered a system comprising two types of chains. One chain of type ABA bearing the simplified sequence (mPDA-6FDA)30-(BAPT-6FDA)17-(mPDA-6FDA)29 resulting in 4323 atoms per chain and another chain of type BAB with the sequence (BAPT-6FDA)18-(mPDA-6FDA)29-(BAPT-6FDA)17 corresponding to a total number of 4551 atoms. The number of monomers of each block is chosen as such in order to obtain a mean molar mass ratio of 15,000/15,000/150,000. A schematic representation of a chain sequence comprising l 6FDA-mPDA blocks and m 6FDA-BAPT blocks is presented in Figure 1. For both systems, the required number of uncorrelated chains was generated via PMC-MD at a temperature just above their glass transition temperature. Specifically, as the experimental measurements of the 6FDA-BAPT and 6FDA-mPDA Tg yield values close to 600 K [40,54,55], the generation temperature was set to 700 K for both polyimide systems. The melts generated by PMC-MD were, thus, obtained at a temperature where the initial configurations should be close to those of the lower-temperature glass [65].
All bond lengths were rigidly constrained to standard values via the SHAKE algorithm in order to ensure an equipartition of the kinetic energy [66]. The simulation time step was set to 1 fs. Polymer chain intra- and intermolecular interactions were parametrised via the TRIPOS 5.2 force field [67], and Lorentz–Berthelot combining rules were applied to all unlike-atom Lennard-Jones (LJ) interactions. Partial charges on the atoms (provided in the Supplementary Materials) were calculated using the Gaussian code [68] at the B3LYP/6-31G** level on representative three-fragment structures of the polyimides under study. Electrostatic interactions were evaluated via an optimised Ewald summation [69]. Dispersive interactions were evaluated for pairs of atoms at distances up to 10 Å for all systems, whereas the energy and pressure beyond the cut-off radius were evaluated by applying standard long-range corrections. The systems generated at 700 K were run for 500 ps with constant-volume constant-temperature dynamics (NVT conditions) in order to allow for local relaxations. A gradual cooling to room temperature (298 K) was performed afterwards, at a rate of −0.5 K/ps. A short relaxation run at constant volume and constant temperature (NVT) followed, and afterwards, the systems were switched to isobaric–isothermal (NPT) conditions. The temperature was fixed at 298 K by loose coupling to a thermostat with a coupling constant of 0.1 ps [70]. The simulation cells were allowed to relax at 1 bar by loose coupling to a barostat using a coupling constant of 5 ps [71]. Production runs were conducted for 4 ns and configurations were saved every 10 ps.

2.2. Densities and Void-Space Analyses of the Pure Polyimide Models

Model densities were evaluated by analysing the last 2 ns of the production trajectories at 298 K. The average densities were found to be 1.306 g cm−3 for 6FDA-BAPT and 1.378 g cm−3 for 6FDA-mPDA/BAPT with error bars being of the order of O(10−3). Wiegand et al. [40] measured an average density of 1.323 g cm−3 for 6FDA-BAPT chains of Mw = 41,800 g mol−1 at 308 K, which is in good agreement with our calculations. No experimental density measurements are available for 6FDA-mPDA/BAPT, but considering the densities of the base polyimides (~1.42 g cm−3 for 6FDA-mPDA [72] and ~1.32 g cm−3 for 6FDA-BAPT [73]) and the fact that the block copolyimide chains are of a 50:50 molar mass ratio, an average model density of ~1.38 g cm−3 is reasonable.
The polyimide void space was first analysed with the aid of a geometric technique similar to the “phantom sphere approach” that allows the determination of the volume accessible to a ghost spherical probe molecule of a preset radius inserted repeatedly into the stored polymer configurations [74,75]. We refer to our approach as the “probe-accessible volume” (PAV) [60]. The insertions of the spherical probes are performed independently, i.e., the current ghost probe does not interact with the other probes. An insertion is considered as “accepted” when the centre of the probe does not overlap with the polymer atoms. The latter is represented by hard spheres with standard van der Waals radii [76,77]. An illustration of the percentage of probe-accessible volume, %PAV, as a function of the probe size, is given in Figure 2 for both pure polymer models under study.
As expected, the %PAV drops systematically with the probe size (Figure 2). For small probe radii (<0.4 Å), the 6FDA-BAPT appears to contain higher %PAV, suggesting that it has more small holes than the copolyimide. The trend changes for larger probe radii, as 6FDA-mPDA/BAPT appears to have more accessible space than the homopolyimide. Its 6FDA-mPDA block with a molecular weight of 15,000 g mol−1 is probably not long enough to induce the more selective behaviour of its homopolyimide, and indeed, it further disrupts the 6FDA-BAPT block. As such, adding the 6FDA-mPDA block leads to larger-size holes, which are likely to induce even more free volume for gas transport in the copolyimide.
To obtain a more visual perception, an illustration of the free-volume distribution for each system is provided in Figure 3. The probe size radius is 1.9 Å, i.e., its diameter is close to the kinetic diameter of both N2 and CH4. This reveals that the 6FDA-BAPT homopolyimide (left picture) is indeed characterised by a rather heterogeneous distribution of accessible volume, the major part of which manifests itself in the form of small-to-medium hole sizes. The coexistence of the mPDA and BAPT diamines in the 6FDA-mPDA/BAPT copolyimide (right picture) increases the accessible space, which, in this case, is dominated by medium-to-large holes. This is consistent with the trends detected in Figure 2 for probe radii greater than 0.4 Å.

3. Modelling Gas Sorption and Diffusion in the Polyimide Matrices

3.1. MD Simulations of N2 and CH4 in the Pure-Gas Phase

Simulations of gas transport were all carried out at 338.5 K since this is a typical temperature for gas wells. As for the chemical bonds in the polymer chains, gas molecules were modelled as rigid to ensure the equipartition of kinetic energy. The parameters for the N2 diatomic molecule were obtained from the symmetrical two-centre Lennard-Jones plus point quadrupole pair potential (2CLJQ) from Vrabec et al. [78]. For CH4, a rigid 5-site model was used, with the potential parameters and partial charges being taken from the model of Yin and Mackerell [79], which has successfully been used in analogous studies [31,58]. Both the N2 and the CH4 force fields used here satisfactorily describe the experimental boiling points and self-diffusion coefficients [56].
To extract the penetrant sorption isotherms in the polymer matrices, single-gas simulations need to be conducted in the pure phase at the temperature and pressures of interest, i.e., a 0–80 bar range. For this reason, a series of NVT MD simulations were performed for 4 ns on systems comprising 512 CH4 and 1000 N2 molecules respectively, using box sizes calculated from experimental data in the NIST database [80].
The last 2 ns of these simulations were used for post-analyses with the configurations stored at intervals of 10 ps. After verifying that the gas systems had reached equilibrium, the gas average pressures, p, were extracted directly from the MD simulations. Solubilities were calculated from the standard Widom test particle insertion (TPI) method [81,82]. In order to obtain the concentrations Cgas(p) and the corresponding solubilities in the gas phase Sgas(p) at any pressure up to 80 bar, the data were fitted to analytical functions having the correct limiting ideal gas behaviour at low pressures. For the binary-gas mixture, the same analyses were carried out on a 2:1 CH4/N2 system with 512 CH4 and 256 N2 molecules, albeit for one set of temperature and pressure only (T = 338.5 K and p = 65.5 bar).

3.2. Single-Gas Uptakes in the Polyimides

In permeation experiments, the number of gas molecules sorbed in the polymer matrix depends on the external pressure, p, of the gas reservoir in contact with the matrix and on the temperature, T. In this work, an approach that iterates the pressure and accounts for the matrix changes during sorption was used to calculate the single-gas sorption isotherms for N2 and CH4 in both 6FDA-BAPT and 6FDA-mPDA/BAPT matrices [61]. The details are given elsewhere [35,62], so only the basic concepts will be presented here.
When gas molecules in the pure-gas phase are in equilibrium with the gas sorbed into the matrix, the chemical potential for the gas in the gas phase,   μ gas , needs to be equal to the chemical potential for the gas in the polymer phase, μ pol . For rigid molecules, this equilibrium condition implies that the difference between the excess chemical potential, μ ex , for the gas in its gas phase, μ ex gas , and in its polymer phase, μ ex pol , is related to the ratio of their concentrations in both phases, with kB being Boltzmann’s constant (Equation (4)) [83]:
Δ μ ex = μ ex pol μ ex gas = k B T   ln C gas C pol
In statistical mechanics, μ e x is given by Equation (5) [60]:
  μ ex = k B T   ln V exp Δ Φ k B T V
where Δ Φ is the change in potential energy estimated using the TPI test particle insertion method, in which a penetrant molecule is virtually and repeatedly introduced into stored configurations of the corresponding model. In dense polymer systems, many millions of particle insertions have to be made in order to reach converged average values for   μ ex . From Equation (5), an expression for the (dimensionless) solubility is obtained (Equation (6)):
S = exp μ ex k B T = V exp Δ Φ k B T V exp Δ Φ k B T
The approximation in Equation (6) holds for fairly incompressible materials such as glassy polymers where volume fluctuations are minor.
As TPI methods are inefficient at high densities due to the low probability of favourable insertions, it was combined here with the excluded-volume map sampling (EVMS) approach [84,85,86], which accelerates the calculations by pre-eliminating the regions of very low insertion probabilities. This is referred to as the EVMS-TPI method. For the solubilities of the gas in the gas phases (Section 3.1), the conventional TPI is sufficient.
Equations (4)–(6) can be combined to give a convenient relationship between the ratios of the concentrations and solubilities in the two phases (Equation (7)):
C gas C pol = exp Δ μ ex k B T = exp μ ex gas k B T exp μ ex pol k B T = S gas S pol     o r       C pol C gas = S pol S gas
Equation (7) is used to determine the single-gas sorption isotherms for both N2 and CH4. For the gas phases, the gas concentration, C gas(p), and solubility, Sgas(p), are obtained for pressures, p, up to ~80 bar at 338.5 K, as described in Section 3.1. For the polymer phases, both models were initially heated to 338.5 K at a rate of 0.5 K/ps. Then, based on the infinite dilution solubility coupled with Cgas(p) and Sgas(p) over a small pressure interval, an initial guess for the number of gas molecules to insert into a given polymer at an initial guess pressure, p1, was made. This was found to be ~50 molecules for N2 and ~20 molecules for CH4. Once inserted with EVMS-TPI, a short NVT 100 ps equilibration run then follows before switching to NPT conditions for a production run of at least 2 ns, the last 1 ns of which is used for post-processing. The average number density, Cpol(p1), and solubility, Spol(p1), of the gas in the polymer phase are determined at pressure p1. Following this, the quantities Cpol(p1)/Cgas(p) and Spol(p1)/Sgas(p) are plotted separately as a function of the pressure, p, of the gas phase. The point of intersection of the two curves gives the second approximation, p2. Then, another simulation of the gas in the polymer phase is carried out at p2, which provides a third estimate, and so on. The maximum number of iterations to reach convergence at 338.5 K was three for both N2 and CH4, as dense polymer matrices are not subject to large changes in volume for moderate pressure increases.

3.3. Binary-Gas Uptakes in the Polyimides

In order to study the sorption of gases under conditions closer to industrial gas separations, the exposure of the polymer membranes to a mixed-gas reservoir was also examined [33]. In this case, the equilibrium condition between the gas phase mixture of fixed composition and the polymer + penetrants system is the same as before, but this now applies to each penetrant. For rigid gas molecules, this implies that the equality of the concentration and solubility ratios, i.e., Cpol/Cgas = Spol/Sgas (Equation (7)), must hold for each type of penetrant.
The iterative technique for single-gas sorption, as described in Section 3.2., has previously been adapted to estimate a binary-gas uptake [56]. A detailed description has been given for the sorption of a binary CH4/N2 mixture in the 6FDA-mPDA, 6FDA-durene, and TR polyimides [59]. The main difference with the single-gas procedure is that, for a binary-gas system, the applied pressure is fixed and the numbers of each type of molecule inserted into the polymer matrix are iterated to convergence.
For the polymer membranes exposed to the binary mixture reservoir, the number of penetrants of both gas types to be initially inserted into the matrices was estimated from the single-gas uptake results. Each polyimide + N2 + CH4 system was relaxed under NVT conditions for 100 ps, and afterwards, NPT runs of at least 4 ns were conducted to allow the systems to reach their equilibrium size and density. The last 2 ns of these runs were used for subsequent analyses. To determine whether penetrants of each gas species had to be added or removed, a comparison of the solubility, Sgas/Spol, and concentration, Cgas/Cpol, ratios for each gas was carried out. The absolute value of the difference between the two aforementioned ratios was iteratively minimised in order to reach convergence. In the systems under study, the difference between the concentration and solubility ratios was minimised within a maximum of five iterations in all models, taking into account errors in calculations. Snapshots of the 6FDA-BAPT homopolyimide both in the pure state as well as loaded with N2 and CH4 at 338.5 K and 65.5 bar are provided in Figure 4. The corresponding snapshots for 6FDA-mPDA/BAPT are provided in the Supplementary Materials.

3.4. Diffusion Coefficients in the Polyimides

As noted above, the permeation studies were carried out at 65.5 bar and 338.5 K for pure N2, pure CH4, and a binary 2:1 mixture of CH4/N2 in both matrices. To estimate their respective diffusion coefficients, the six simulations at 65.5 bar and 338.5 K were extended up to 15–20 ns with NPT MD, which allowed the Fickian regime for diffusion to be attained (see later). On these grounds, the self-diffusion coefficients, D, could be extracted using Einstein’s equation (Equation (8)):
D = lim t r i t + t 0 r i t 0 2 6 t
using the mean square displacements, MSD = <(ri (t + t0) – ri(t0))2>, of a given type of penetrant averaged over all time origins, t0, and over all penetrants of the same type.

4. Results and Discussion

4.1. Single-Gas Sorption of N2 and CH4 in the Triptycene-Based Polyimides

In this section, we present the results obtained from the pure-gas simulations at 338.5 K in which N2 and CH4 were gradually loaded separately in the polyimides’ matrices and the sorption isotherms that were obtained over the pressure range of 0–80 bar through the iterative method described in Section 3.2. In order to be consistent with experiments, the gas concentrations in the polymer at a gas reservoir pressure p, Cpol(p), are given in terms of the volume that the number of gas molecules in the matrix, npol(p), would occupy if they were an ideal gas at the standard temperature and pressure (STP) conditions, divided by the volume of the pure polymer volume, V0. These nominal concentrations, C0(p), do not take into account the limited volume swelling, which is very difficult to measure experimentally [88,89] (Equation (9)):
C 0 p = n pol p   k B   T STP V 0   p STP
where pSTP = 1.01325 bar and TSTP = 273.15 K are the standard temperature and pressure, respectively. The nominal solubility coefficient is then defined as (Equation (10)):
S c 0 p = C 0 p p
The nominal solubility coefficients are plotted with symbols for both gases and both matrices in Figure 5 as a function of p. The limiting values of the solubility coefficient for the penetrant molecule in the matrices can be obtained from (Equation (11)) [90]:
lim p 0 S c 0 p = S pol   T STP T   p STP
In addition, the added smooth curves in Figure 5 are nonlinear least-squares regression fits to the dual-mode sorption (DMS) model (Equation (12)) [58]:
S c 0 p = k D + C H   b 1 + b p
according to which gas sorption in glassy polymers is manifested both by a “hole-filling” mechanism described by the Langmuir sorption capacity, C H , and the hole affinity, b, as well as by the dissolution into the dense regions obeying Henry’s law characterised by the constant kD. A drawback of the DMS model is that it is difficult to calculate its parameters from the known properties of the polymers and the dissolved gases [91]. However, the functional form of Equation (12) is often used to fit experimental single-gas uptake curves. Figure 5 clearly shows that, as in the experiment, the model results can be fitted by DMS too.
In accordance with other studies [56,59], methane is more soluble than nitrogen in both matrices. The block copolyimide exhibits consistently higher solubilities over the whole pressure range for both gases, which is consistent with the larger holes in 6FDA-mPDA/BAPT compared to 6FDA-BAPT (Figure 2 and Figure 3). Nevertheless, the differences between both matrices are rather small. This is expected since it is well known that polyimides usually have fairly close penetrant solubilities and that the separations of gases in these polymers are often mainly due to the diffusion selectivities [92].
The TPI analyses provide the weighted mean test particle insertion energy (Equation (13)):
Δ Φ = Δ Φ     V e x p Δ Φ k B T V e x p Δ Φ k B T
To ensure that the averages in Equation (13) represent the typical mean interaction energies that the probe molecule would have in the matrices, they can be weighted by the Boltzmann’s factor [59]. In general, the absolute values of the mean insertion energies for both systems follow the same hierarchy as the solubility coefficients, and indeed, that for   nitrogen   (≈15 kT is lower than that for methane (≈20 kT. They differ marginally between the two models.
From the TPI analyses, it is also possible to measure the fraction of the “significant accessible volume” (SAV) (also sometimes called “fraction of significant volume” (FSV) [89]), which is actually available for the insertion of a probe molecule. Unlike the geometric criterion used in Section 2.2, this approach is based on an energetic criterion, i.e., it estimates the available volume that is responsible for 99.9% of the solubility [60]. In the case of incompressible systems, the probability densities of the insertion energies, ρ Δ Φ , can be related to the solubility via (Equation (14)) [60]:
+ ρ Δ ϕ e x p Δ ϕ k B T d Δ ϕ
The associated integrals of these functions reach a plateau relatively quickly as the Boltzmann’s weighting factor rapidly diminishes with the increase in insertion energy. The critical upper limit, ΔΦc, for the solubility (by which the integrals are converged to more than 99.9%) is first determined. Integrating the probability density distribution functions up to this critical energy will then yield the fraction of SAV for the molecule being inserted (Equation (15)) [60]:
S A V = Δ Φ c ρ Δ Φ d Δ Φ
Given that this approach implicitly takes into account the interactions of each probe with its matrix, it allows for a better estimation of the volume really used by the probe molecule in question. Figure 6 presents the %SAV in both the pure matrices (C0(p) = 0) and in the polymer + penetrants systems as a function of the penetrant nominal concentrations over the whole range of p studied for the single-gas simulations.
The relative order of the %SAV between both polyimides is in accordance with the predictions from the phantom sphere approach for large probe radii in pure polymers (Figure 2 and Figure 3). This shows, again, that the incorporation of the mPDA diamine into the block copolyimide results in a significant gain in accessible volume with respect to 6FDA-BAPT. The %SAV is higher for the less soluble N2 penetrant, which means that it needs to reach more sites in the matrix than CH4 to account for its solubility. On the other hand, the latter favours those sites with the most favourable energetic interactions. The %SAV drops with penetrant concentrations as sites are progressively filled up.

4.2. Volume Dilations upon Single-Gas Sorption

Gas-sorption-induced volume dilation effects are difficult to capture experimentally [67,68]. In contrast, volume swelling can be extracted easily in MD simulations in 3D periodic boundary conditions by calculating the volume change in the simulation box with respect to the pure polymer volume, V0 [56,93].
To our knowledge, there are no available volume dilation data for N2 or CH4 in triptycene-based polyimides. However, a few studies examined volume swelling in glassy polymers due to these penetrants [88,94,95]. They report CH4-induced volume swelling up to 2% at pressures from 60 to 70 bar [88,94], whereas the volume changes induced by N2 in a glassy copolymer reach 0.6% at 30 bar [95]. Methane can, thus, be considered as moderately plasticizing, whereas nitrogen is non-plasticizing.
In a previous simulation of the N2/CH4 gas pair in fluorinated polyimides [56], CH4-induced volume dilation reached ~5% at pressures approaching 70 bar, whereas there was limited swelling upon N2 sorption. The dilation effects in a 6FDA-6FpDA polyimide were found to be ~2.6% for CH4 and ~1.4% for N2 at 60 bar [96]. Another work addressing the sorption of N2 and CH4 in thermally rearranged (TR) polymers reported a volume swelling for the 6FDA-bisAPAF polyimide precursor of ~3.6% for CH4 and ~1.8% for N2 at pressures of ~63 and ~73 bar, respectively, whereas both gases induced negligible dilations for the TR structure [59]. In the present case, N2 and CH4 induce volume changes in 6FDA-BAPT of about 1.3% and 1.8%, respectively, at the highest gas loadings, i.e., at feed pressures approaching 60–70 bar. These findings support experimental work demonstrating the high plasticisation resistance for triptycene-based polyimides [39,40] and indicate that 6FDA-BAPT is less prone to volume dilation than other non-thermally rearranged 6FDA-based polyimides [55]. For the 6FDA-mPDA/BAPT copolyimide, N2 sorption induces negligible swelling, i.e., less than 0.4%, which is probably due to its larger holes (Figure 2) but CH4 swells the matrix up to 3% at a pressure of 64 bar. The latter can be attributed to the mPDA diamine structure, which is known to be less resistant to plasticisation than the more rigid and bulky BAPT counterpart [62].

4.3. Binary-Gas Sorption of a 2:1 CH4/N2 Mixture in the Triptycene-Based Polyimides

When in equilibrium with a 2:1 binary-gas CH4/N2 mixture at 65.5 bar and at 338.5 K, the mixed-gas compositions in the polyimide models were obtained via the implementation of the method described in Section 3.3. The nominal concentrations, C0, issued from the MD simulations are presented in Table 1, along with the C0 predictions obtained from the extension of the DMS model to binary mixtures [97,98] (Equation (16)):
C 0 A p = K D A p A + C H A b A p A     1 + b A p A + b B p B    
where A and B are the two different gases, and pA and pB are the partial pressures of each gas. The constants for each gas are estimated from the DMS fits to the single-gas sorption curves (Figure 5).
Slight differences in mixed-gas compositions are observed between 6FDA-BAPT and 6FDA-mPDA/BAPT. Both gases exhibit higher concentrations in the block copolyimide, in accordance with the single-gas uptake curves presented in Figure 5. The molar gas ratios in the polymer phases are higher than 4:1, deviating both from the ratios of CH4 and N2 molecules based on the pure-gas solubilities (Figure 5), as well as from the 2:1 molar ratio of the pure binary-gas phase. A similar trend is observed for the nominal C0 concentrations obtained from the DMS theory extended to binary mixtures (Equation (16)), which are in reasonable agreement with the simulated C0 (Table 1). The higher sorption of methane over nitrogen, which is expected from the higher solubility of the former and from its higher molar fraction in the binary-gas phase, leads to a similar volume swelling of ~2% for both matrices, with the degree of volume swelling being intermediate between the single-gas swellings (Section 4.2).
These results lie in accordance with analogous studies on N2/CH4 mixed-gas sorption in thermally rearranged polymer membranes [59]. The data in Table 1 also demonstrate the robustness of the DMS theory for binary mixtures as it allows for a fairly satisfactory prediction of the penetrant concentrations in the polymer phase without the need to conduct mixed-gas simulations. However, the adequacy of the simulated C0 is both matrix- and penetrant-dependent, and in the case of the block copolyimide, the DMS theory does underestimate the ratio of the equilibrium gas concentrations by ~10–20%.
To examine possible preferential interactions between the polymer and one of both gases, intermolecular radial distribution functions, ginter(r), were calculated between the nitrogen centre of mass, the methane carbon and the different types of polymer atoms, with averages being taken over all frames of the production runs. These functions were, thus, analysed in the case of all the mixed-gas systems; examples are given for the 6FDA-mPDA/BAPT + N2 + CH4 system comprising the gas concentrations listed in Table 1.
As can be inferred from Figure 7, methane shows a higher probability of interaction with the polymer aromatic carbon than nitrogen, whereas the inverse tendency is observed when the gases interact with the fluorocarbon group. This implies that not all polymer sites interact with the gas components in the same manner.

4.4. Single-Gas vs. 2:1 Binary-Gas CH4/N2 Gas Diffusivities, Permeabilities, and Permselectivities in the Triptycene-Based Polyimides

This section compares the diffusion, permeation, and separation abilities for both mixed-gas systems at 65.5 bar and 338.5 K with the pure-gas systems. As the simulations for the latter were performed at a fixed number of penetrants, the pure-gas systems closest in pressure to 65.5 bar were used for the comparison.
As explained in Section 3.4, the gas diffusion coefficients were extracted from the MSD of the gas centres of mass, which were calculated from MD trajectories extrapolated up to at least 15 ns. Figure 8 presents the MSD/6t vs. simulation time, t, for the gases in the 6FDA-BAPT matrix and shows that each MSD/6t reaches a plateau value, which can provide a reasonable estimation of D (Equation (8)). This is also the case for 6FDA-mPDA/BAPT, although the plateau values were attained a bit faster, as the penetrants were more mobile in the copolyimide.
The gas permeabilities, P, were estimated using Equation (1) from the product of the Sc solubility coefficients extracted from Equation (2) and the diffusion coefficients, D, obtained from Equation (8). The separation factors or permselectivities, α, were obtained from Equation (3). Table 2 lists the gas solubility, diffusion, and permeability in the mixed-gas systems at 65.5 bar together with the corresponding data of the single-gas systems at the feed pressures closest to this value. The mixed N2/CH4 separation factors are obtained from the binary-gas runs, and the pure N2/CH4 separation factors are the ideal permselectivities from the single-gas runs. We note here that the Sc are reported in units of cm3(STP)/cm3(cm Hg) as this is a straightforward way to calculate the P in Barrer afterwards (1 Barrer = 10−10 cm3(STP) cm2/(cm3 s (cm Hg)). D are provided in units of cm2 s−1. The maximum errors on Sc are 0.0001 cm3(STP)/cm3(cm Hg), while those on D are larger due to the plateau values in Figure 8. They are all typically of the order of 10−7 cm2 s−1. As such, only the specific errors on P and on the separation factors are reported in Table 2.
Upon mixing, there is a slight decrease in the solubility coefficients, Sc, for N2 in both matrices, which is to be expected because of the competition with the more soluble CH4. Indeed, the Sc of CH4 increases simultaneously from ~1.5 in the single gas to ~2.5 larger than for N2 in the mixed-gas systems. The CH4/N2 solubility selectivity is, thus, enhanced under mixed-gas conditions. However, as observed in the sorption isotherms (Figure 5), the Sc for both gas components exhibits only small differences between both 6FDA-BAPT and 6FDA-BAPT/mPDA. Since the permeabilities, P, are quite different, this suggests that the kinetic aspect also plays a role [37], which is similar to analogous studies on fluorinated polyimides [56].
N2 diffuses faster than CH4, as expected from its slightly smaller kinetic diameter [27]. In absolute terms, the fastest diffusivities, D, and, hence, the highest gas permeabilities, P, are observed in the 6FDA-BAPT/mPDA copolyimide. Indeed, the incorporation of the mPDA diamine induces an increase in the significant accessible volume for both gas penetrants, thus enhancing their overall mobilities by a factor of 2–3 with respect to 6FDA-BAPT. However, in both matrices, the N2/CH4 diffusion selectivities are all less than ~2.
The ideal permselectivity of 6FDA-BAPT for the N2/CH4 pair has been reported experimentally to be ~1.8 at T = 35° and p = 9 bar. Although the conditions are different, this is quite coherent with the ideal permselectivity of our model [40] (Table 2). Due to its higher available volume, 6FDA-mPDA/BAPT seems to lose its permselectivity under single-gas conditions.
However, both matrices behave in different ways in the mixed-gas systems. The enhancements in solubility selectivities lead to larger CH4 and lower N2 permeabilities. This behaviour has also been detected in earlier experimental studies [97,99], as well as in simulations of fluorinated and TR polyimides [56,59]. Both triptycene-based polyimides, thus, favour the transport of methane over that of nitrogen under the mixed-gas conditions studied (2:1 CH4:N2 binary mixtures at 65.5 bar and 338.5 K). It is, thus, important to test gas mixtures in addition to single-gas studies [33].

5. Conclusions

MD simulations were performed to compare the transport of N2 and CH4 in two triptycene-based polyimides. The first one was the 6FDA-BAPT homopolyimide, and the second one was a block 6FDA-mPDA/BAPT copolyimide, combining the high FFV of the triptycene diamine with the mPDA diamine, with the latter expected to be more selective. The densities of the pure polyimide models were found to be in agreement with experiment. However, the void-space analyses revealed that, while the free volume consists of small-to-medium holes in 6FDA-BAPT, there are more medium-to-large holes in 6FDA-mPDA/BAPT. Indeed, the 6FDA-mPDA block with a molecular weight of 15,000 g mol−1 is probably not long enough to induce a selective behaviour, and it further disrupts the 6FDA-BAPT block, hence leading to a copolyimide with more free volume.
Following the generation of the polymer matrices, a single-gas sorption technique iterating the pressure for a constant number of penetrants was applied to obtain the sorption isotherms for N2 and CH4 over a 0–70 bar range at 338.5 K. Both gases were found to be more soluble in the block copolyimide with a higher affinity for methane. This was confirmed by the fraction of significant accessible volume (%SAV), which showed that CH4 favours sites with the most favourable energetic interactions, unlike the lower-solubility N2, which probes more sites in the matrices. 6FDA-BAPT swelled less than 2% at the highest gas loadings, while 6FDA-mPDA/BAPT swelled up to 3% due to its less resistant mPDA diamine [39].
An adapted mixed-gas sorption technique iterating the number of gas penetrants for a constant external gas pressure was used to estimate the mixed-gas equilibrium concentrations of a 2:1 CH4/N2 mixture at 338.5 K and 65.5 bar in both matrices. The molar CH4/N2 ratios in the polymer phases were found to be higher than 4:1, i.e., with an enhancement in the solubility selectivity in favour of CH4. N2 and CH4 diffusivities were probed via the mean square displacements. Although diffusion was faster for N2 and generally higher in 6FDA-BAPT/mPDA, the diffusion selectivities were not significantly different from each other. As such, the higher solubility selectivities meant that the triptycene-based polyimides both favoured, to a similar extent, the transport of methane over that of nitrogen under the conditions studied.
Concerning the permselectivity performance of the triptycene-based PIs examined in this study as compared to other 6FDA-based PI membranes previously examined by our group, it appears that all polymer matrices behave in different ways in the mixed-gas systems compared to the single-gas case. 6FDA-mPDA/BAPT is highly permeable to methane in the mixed-gas case, which is similar to the behaviour observed in a thermally rearranged TR-PBO polyimide [59]. Since the mPDA diamine is methane-selective, when it is combined with nitrogen-selective diamines, such as the durene diamine, it induces a reduction in the nitrogen-selectivity of the latter [56]. Finally, the BAPT polyimide appears reasonably more resistant to plasticisation compared to the mPDA and durene diamines, with TR-PBO being the only polymer matrix that undergoes lower volume dilation under the same operating conditions.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/polym15183811/s1.

Author Contributions

Conceptualization, D.B.; Methodology, I.T. and D.B.; Software, D.B.; Validation, I.T.; Formal analysis, I.T.; Investigation, I.T.; Resources, D.B.; Data curation, I.T., D.B. and S.N.; Writing—original draft, I.T.; Writing—review & editing, D.B. and S.N.; Visualization, I.T. and D.B.; Supervision, D.B. and S.N.; Project administration, D.B., J.-P.B., M.V., S.D. and A.B.; Funding acquisition, D.B. and S.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Saudi Aramco. Contract number: 6600025762.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to restrictions imposed by the contract number 6600025762.

Acknowledgments

We gratefully thank Saudi Aramco for providing a 3-year postdoctoral grant to I.T., as well as for the generous financial assistance for this project. This work was granted access to the HPC resources of CCRT/CINES/IDRIS made by GENCI (Grand Equipment National de Calcul Intensif), France. The MUST cluster at the University Savoie Mont Blanc (France) is also acknowledged for the allocation of computer time. Nicolas Charvin is thanked for his help with the installation and maintenance of the laboratory local servers.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, Z.; Fuoco, A.; He, G. Membranes for Gas Separation. Membranes 2021, 11, 755. [Google Scholar] [CrossRef] [PubMed]
  2. Sridhar, S.; Smitha, B.; Aminabhavi, T.M. Separation of Carbon Dioxide from Natural Gas Mixtures through Polymeric Membranes—A Review. Sep. Purif. Rev. 2007, 36, 113–174. [Google Scholar] [CrossRef]
  3. Baker, R.W. Future Directions of Membrane Gas Separation Technology. Ind. Eng. Chem. Res. 2002, 41, 1393–1411. [Google Scholar] [CrossRef]
  4. Sholl, D.S.; Lively, R.P. Seven chemical separations to change the world. Nature 2016, 532, 435–437. [Google Scholar] [CrossRef] [PubMed]
  5. Figoli, A.; Simone, S. Membranes for sustainable food packaging. Membr. Process. Sustain. Growth 2013, 101–119. [Google Scholar]
  6. Korpyś, M.; Wójcik, J.; Synowiec, P. Methods for sweetening natural and shale gas. Chem. Sci. Tech. Mark. 2014, 68, 213–215. [Google Scholar]
  7. Sanaeepur, H.; Ebadi Amooghin, A.; Bandehali, S.; Moghadassi, A.; Matsuura, T.; Van der Bruggen, B. Polyimides in membrane gas separation: Monomer’s molecular design and structural engineering. Prog. Polym. Sci. 2019, 91, 80–125. [Google Scholar] [CrossRef]
  8. Maier, G. Gas Separation by Polymer Membranes: Beyond the Border. Angew. Chem. Int. Ed. 2013, 52, 4982–4984. [Google Scholar] [CrossRef]
  9. Siagian, U.W.R.; Raksajati, A.; Himma, N.F.; Khoiruddin, K.; Wenten, I.G. Membrane-based carbon capture technologies: Membrane gas separation vs. membrane contactor. J. Nat. Gas Sci. Eng. 2019, 67, 172–195. [Google Scholar] [CrossRef]
  10. Castel, C.; Bounaceur, R.; Favre, E. Engineering of Membrane Gas Separation Processes: State of The Art and Prospects. J. Membr. Sci. Res. 2020, 6, 295–303. [Google Scholar] [CrossRef]
  11. Sidhikku Kandath Valappil, R.; Ghasem, N.; Al-Marzouqi, M. Current and future trends in polymer membrane-based gas separation technology: A comprehensive review. J. Ind. Eng. Chem. 2021, 98, 103–129. [Google Scholar] [CrossRef]
  12. Iulianelli, A.; Drioli, E. Membrane engineering: Latest advancements in gas separation and pre-treatment processes, petrochemical industry and refinery, and future perspectives in emerging applications. Fuel Process. Technol. 2020, 206, 106464. [Google Scholar] [CrossRef]
  13. Chen, X.Y.; Vinh-Thang, H.; Ramirez, A.A.; Rodrigue, D.; Kaliaguine, S. Membrane gas separation technologies for biogas upgrading. RSC Adv. 2015, 5, 24399–24448. [Google Scholar] [CrossRef]
  14. Lasseuguette, E.; Comesaña-Gándara, B. Polymer Membranes for Gas Separation. Membranes 2022, 12, 207. [Google Scholar] [CrossRef] [PubMed]
  15. Liang, C.Z.; Chung, T.-S.; Lai, J.-Y. A review of polymeric composite membranes for gas separation and energy production. Prog. Polym. Sci. 2019, 97, 101141. [Google Scholar] [CrossRef]
  16. Lai, H.; Benedetti, F.M.; Ahn, J.; Robinson, A.; Wang, Y.; Pinnau, I.; Smith, Z.; Xia, Y. Hydrocarbon ladder polymers with ultrahigh permselectivity for membrane gas separations. Science 2022, 375, 1390–1392. [Google Scholar] [CrossRef] [PubMed]
  17. Yang, J.; Tao, L.; He, J.; McCutcheon, J.R.; Li, Y. Machine learning enables interpretable discovery of innovative polymers for gas separation membranes. Sci. Adv. 2022, 8, eabn9545. [Google Scholar] [CrossRef]
  18. Chen, G.; Shen, Z.; Iyer, A.; Ghumman, U.F.; Tang, S.; Bi, J.; Chen, W.; Li, Y. Machine-Learning-Assisted De Novo Design of Organic Molecules and Polymers: Opportunities and Challenges. Polymers 2020, 12, 163. [Google Scholar] [CrossRef]
  19. Audus, D.J.; de Pablo, J.J. Polymer Informatics: Opportunities and Challenges. ACS Macro Lett. 2017, 6, 1078–1082. [Google Scholar] [CrossRef]
  20. Wang, J.; Tian, K.; Li, D.; Chen, M.; Feng, X.; Zhang, Y.; Wang, Y.; Van der Bruggen, B. Machine learning in gas separation membrane developing: Ready for prime time. Sep. Purif. Technol. 2023, 313, 123493. [Google Scholar] [CrossRef]
  21. Tayyebi, A.; Alshami, A.S.; Yu, X.; Kolodka, E. Can machine learning methods guide gas separation membranes fabrication? J. Membr. Sci. Lett. 2022, 2, 100033. [Google Scholar] [CrossRef]
  22. Madaeni, S.S.; Hasankiadeh, N.T.; Kurdian, A.R.; Rahimpour, A. Modeling and optimization of membrane fabrication using artificial neural network and genetic algorithm. Sep. Purif. Technol. 2010, 76, 33–43. [Google Scholar] [CrossRef]
  23. Rall, D.; Schweidtmann, A.M.; Aumeier, B.M.; Kamp, J.; Karwe, J.; Ostendorf, K.; Mitsos, A.; Wessling, M. Simultaneous rational design of ion separation membranes and processes. J. Membr. Sci. 2020, 600, 117860. [Google Scholar] [CrossRef]
  24. Rufford, T.E.; Smart, S.; Watson, G.C.Y.; Graham, B.F.; Boxall, J.; Diniz da Costa, J.C.; May, E.F. The removal of CO2 and N2 from natural gas: A review of conventional and emerging process technologies. J. Pet. Sci. Eng. 2012, 94–95, 123–154. [Google Scholar] [CrossRef]
  25. Baker, R.W.; Lokhandwala, K. Natural Gas Processing with Membranes:  An Overview. Ind. Eng. Chem. Res. 2008, 47, 2109–2121. [Google Scholar] [CrossRef]
  26. Butwell, K.F.; Dolan, W.B.D.; Kuznicki, S.M. Selective Removal of Nitrogen from Natural Gas by Pressure Swing Adsorption. US Patent US6444012B1, 3 September 2002. [Google Scholar]
  27. Robeson, L.M.; Freeman, B.D.; Paul, D.R.; Rowe, B.W. An empirical correlation of gas permeability and permselectivity in polymers and its theoretical basis. J. Membr. Sci. 2009, 341, 178–185. [Google Scholar] [CrossRef]
  28. Robeson, L.M. The upper bound revisited. J. Membr. Sci. 2008, 320, 390–400. [Google Scholar] [CrossRef]
  29. Murali, R.S.; Sankarshana, T.; Sridhar, S. Air Separation by Polymer-based Membrane Technology. Sep. Purif. Rev. 2013, 42, 130–186. [Google Scholar] [CrossRef]
  30. Edens, S.J.; McGrath, M.J.; Guo, S.; Du, Z.; Zhou, H.; Zhong, L.; Shi, Z.; Wan, J.; Bennett, T.D.; Qiao, A.; et al. An Upper Bound Visualization of Design Trade-Offs in Adsorbent Materials for Gas Separations: CO2, N2, CH4, H2, O2, Xe, Kr, and Ar Adsorbents. Adv. Sci. 2023, 10, 2206437. [Google Scholar] [CrossRef]
  31. Wijmans, J.G.; Baker, R.W. The solution-diffusion model: A review. J. Membr. Sci. 1995, 107, 1–21. [Google Scholar] [CrossRef]
  32. Wijmans, J.G.; Baker, R.W. The Solution-Diffusion Model: A Unified Approach to Membrane Permeation. In Materials Science of Membranes for Gas and Vapor Separation; Yampolskii, Y., Pinnau, I., Freeman, B., Eds.; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2006; pp. 159–189. [Google Scholar]
  33. Sanders, D.F.; Smith, Z.P.; Guo, R.; Robeson, L.M.; McGrath, J.E.; Paul, D.R.; Freeman, B.D. Energy-efficient polymeric gas separation membranes for a sustainable future: A review. Polymer 2013, 54, 4729–4761. [Google Scholar] [CrossRef]
  34. Xu, S.; Ren, X.; Zhao, N.; Wu, L.; Zhang, Z.; Fan, Y.; Li, N. Self-crosslinking of bromomethylated 6FDA-DAM polyimide for gas separations. J. Membr. Sci. 2021, 636, 119534. [Google Scholar] [CrossRef]
  35. Xu, Z.; Croft, Z.L.; Guo, D.; Cao, K.; Liu, G. Recent development of polyimides: Synthesis, processing, and application in gas separation. J. Polym. Sci. 2021, 59, 943–962. [Google Scholar] [CrossRef]
  36. Ghosh, M.K.; Mittal, K.L. Polyimides: Fundamentals and Applications; Marcel Dekker, Inc.: New York, NY, USA, 1996. [Google Scholar]
  37. Kim, T.H.; Koros, W.J.; Husk, G.R.; O’Brien, K.C. Relationship between gas separation properties and chemical structure in a series of aromatic polyimides. J. Membr. Sci. 1988, 37, 45–62. [Google Scholar] [CrossRef]
  38. White, L.S.; Blinka, T.A.; Kloczewski, H.A.; Wang, I.f. Properties of a polyimide gas separation membrane in natural gas streams. J. Membr. Sci. 1995, 103, 73–82. [Google Scholar] [CrossRef]
  39. Weidman, J.R. Triptycene as an Architectural Motif in the Macromolecular Design of Polyimides for Gas Separation Membranes. Ph.D. Thesis, University of Notre Dame, Notre Dame, IN, USA, 2018. [Google Scholar]
  40. Wiegand, J.R.; Smith, Z.P.; Liu, Q.; Patterson, C.T.; Freeman, B.D.; Guo, R. Synthesis and characterization of triptycene-based polyimides with tunable high fractional free volume for gas separation membranes. J. Mater. Chem. A 2014, 2, 13309–13320. [Google Scholar] [CrossRef]
  41. Sydlik, S.A.; Chen, Z.; Swager, T.M. Triptycene Polyimides: Soluble Polymers with High Thermal Stability and Low Refractive Indices. Macromolecules 2011, 44, 976–980. [Google Scholar] [CrossRef]
  42. Swaidan, R.J.; Ghanem, B.; Swaidan, R.; Litwiller, E.; Pinnau, I. Pure- and mixed-gas propylene/propane permeation properties of spiro- and triptycene-based microporous polyimides. J. Membr. Sci. 2015, 492, 116–122. [Google Scholar] [CrossRef]
  43. Ghanem, B.S.; Swaidan, R.; Litwiller, E.; Pinnau, I. Ultra-Microporous Triptycene-based Polyimide Membranes for High-Performance Gas Separation. Adv. Mater. 2014, 26, 3688–3692. [Google Scholar] [CrossRef]
  44. Luo, S.; Wiegand, J.R.; Kazanowska, B.; Doherty, C.M.; Konstas, K.; Hill, A.J.; Guo, R. Finely Tuning the Free Volume Architecture in Iptycene-Containing Polyimides for Highly Selective and Fast Hydrogen Transport. Macromolecules 2016, 49, 3395–3405. [Google Scholar] [CrossRef]
  45. Hsiao, S.-H.; Wang, H.-M.; Chen, W.-J.; Lee, T.-M.; Leu, C.-M. Synthesis and properties of novel triptycene-based polyimides. J. Polym. Sci. Part A Polym. Chem. 2011, 49, 3109–3120. [Google Scholar] [CrossRef]
  46. Hsiao, S.-H.; Guo, W.; Tsai, T.-H.; Chiu, Y.-T. Synthesis of soluble and thermally stable triptycene-based poly(amide-imide)s. J. Polym. Res. 2014, 21, 391. [Google Scholar] [CrossRef]
  47. Cho, Y.J.; Park, H.B. High Performance Polyimide with High Internal Free Volume Elements. Macromol. Rapid Commun. 2011, 32, 579–586. [Google Scholar] [CrossRef] [PubMed]
  48. Weidman, J.R.; Luo, S.; Doherty, C.M.; Hill, A.J.; Gao, P.; Guo, R. Analysis of governing factors controlling gas transport through fresh and aged triptycene-based polyimide films. J. Membr. Sci. 2017, 522, 12–22. [Google Scholar] [CrossRef]
  49. Luo, S.; Liu, J.; Lin, H.; Kazanowska, B.A.; Hunckler, M.D.; Roeder, R.K.; Guo, R. Preparation and gas transport properties of triptycene-containing polybenzoxazole (PBO)-based polymers derived from thermal rearrangement (TR) and thermal cyclodehydration (TC) processes. J. Mater. Chem. A 2016, 4, 17050–17062. [Google Scholar] [CrossRef]
  50. Loianno, V.; Zhang, Q.; Luo, S.; Guo, R.; Galizia, M. Modeling Gas and Vapor Sorption and Swelling in Triptycene-Based Polybenzoxazole: Evidence for Entropy-Driven Sorption Behavior. Macromolecules 2019, 52, 4385–4395. [Google Scholar] [CrossRef]
  51. Chen, Y.-R.; Chen, L.-H.; Chang, K.-S.; Chen, T.-H.; Lin, Y.-F.; Tung, K.-L. Structural characteristics and transport behavior of triptycene-based PIMs membranes: A combination study using ab initio calculation and molecular simulations. J. Membr. Sci. 2016, 514, 114–124. [Google Scholar] [CrossRef]
  52. Ghasemnejad-Afshar, E.; Amjad-Iranagh, S.; Zarif, M.; Modarress, H. Effect of side branch on gas separation performance of triptycene based PIM membrane: A molecular simulation study. Polym. Test. 2020, 83, 106339. [Google Scholar] [CrossRef]
  53. Balcik, M.; Wang, Y.; Pinnau, I. Exploring the effect of intra-chain rigidity on mixed-gas separation performance of a Triptycene-Tröger’s base ladder polymer (PIM-Trip-TB) by atomistic simulations. J. Membr. Sci. 2023, 677, 121614. [Google Scholar] [CrossRef]
  54. Heck, R.; Qahtani, M.S.; Yahaya, G.O.; Tanis, I.; Brown, D.; Bahamdan, A.A.; Ameen, A.W.; Vaidya, M.M.; Ballaguet, J.P.R.; Alhajry, R.H.; et al. Block copolyimide membranes for pure- and mixed-gas separation. Sep. Purif. Technol. 2017, 173, 183–192. [Google Scholar] [CrossRef]
  55. Wang, L.; Cao, Y.; Zhou, M.; Ding, X.; Liu, Q.; Yuan, Q. The gas permeation properties of 6FDA-2, 4, 6-trimethyl-1, 3-phenylenediamine (TMPDA)/1, 3-phenylenediamine (mPDA) copolyimides. Polym. Bull. 2008, 60, 137–147. [Google Scholar] [CrossRef]
  56. Tanis, I.; Brown, D.; Neyertz, S.; Heck, R.; Mercier, R.; Vaidya, M.; Ballaguet, J.-P. A comparison of pure and mixed-gas permeation of nitrogen and methane in 6FDA-based polyimides as studied by molecular dynamics simulations. Comput. Mater. Sci. 2018, 141, 243–253. [Google Scholar] [CrossRef]
  57. O’Brien, K.C.; Koros, W.J.; Barbari, T.A.; Sanders, E.S. A new technique for the measurement of multicomponent gas transport through polymeric films. J. Membr. Sci. 1986, 29, 229–238. [Google Scholar] [CrossRef]
  58. Yampolskii, Y.P.; Pinnau, I.; Freeman, B.D. Materials Science of Membranes for Gas and Vapour Separation; John Wiley & Sons Ltd.: Hoboken, NJ, USA, 2006. [Google Scholar]
  59. Tanis, I.; Brown, D.; Neyertz, S.; Vaidya, M.; Ballaguet, J.-P.; Duval, S.; Bahamdan, A. 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. Phys. Chem. Chem. Phys. 2022, 24, 18667–18683. [Google Scholar] [CrossRef] [PubMed]
  60. Brown, D. The gmq User Manual Version 6. 2021. Available online: http://www.lmops.univ-savoie.fr/brown/gmq.html (accessed on 10 April 2023).
  61. Neyertz, S.; Douanne, A.; Brown, D. Effect of Interfacial Structure on Permeation Properties of Glassy Polymers. Macromolecules 2005, 38, 10286–10298. [Google Scholar] [CrossRef]
  62. Tanis, I.; Brown, D.; Neyertz, S.J.; Heck, R.; Mercier, R. A comparison of homopolymer and block copolymer structure in 6FDA-based polyimides. Phys. Chem. Chem. Phys. 2014, 16, 23044–23055. [Google Scholar] [CrossRef] [PubMed]
  63. Neyertz, S.; Brown, D. Preparation of bulk melt chain configurations of polycyclic polymers. J. Chem. Phys. 2001, 115, 708–717. [Google Scholar] [CrossRef]
  64. Neyertz, S.; Brown, D.; Clarke, J.H.R. The local energy approximation and the predictability of chain configurations in polymer melts. J. Chem. Phys. 1996, 105, 2076–2088. [Google Scholar] [CrossRef]
  65. Neyertz, S. Tutorial: Molecular Dynamics Simulations of Microstructure and Transport Phenomena in Glassy Polymers. Soft Mater. 2007, 4, 15–83. [Google Scholar] [CrossRef]
  66. Ciccotti, G.; Ferrario, M.; Ryckaert, J.P. Molecular dynamics of rigid systems in cartesian coordinates A general formulation. Mol. Phys. 1982, 47, 1253–1264. [Google Scholar] [CrossRef]
  67. Clark, M.; Cramer, R.D.; Van Opdenbosch, N. Validation of the general purpose tripos 5.2 force field. J. Comput. Chem. 1989, 10, 982–1012. [Google Scholar] [CrossRef]
  68. Frisch, M.J.T.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Montgomery, J.A., Jr.; Vreven, T.; Kudin, K.N.; Burant, J.C.; et al. Gaussian, Inc. Available online: www.gaussian.com (accessed on 5 April 2023).
  69. Ewald, P.P. Die Berechnung optischer und elektrostatischer Gitterpotentiale. Ann. Phys. 1921, 369, 253–287. [Google Scholar] [CrossRef]
  70. Berendsen, H.J.C.; Postma, J.P.M.; van Gunsteren, W.F.; DiNola, A.; Haak, J.R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1984, 81, 3684–3690. [Google Scholar] [CrossRef]
  71. Brown, D.; Neyertz, S. A general pressure tensor calculation for molecular dynamics simulations. Mol. Phys. 1995, 84, 577–595. [Google Scholar] [CrossRef]
  72. Cheng, S.-X.; Chung, T.-S.; Wang, R.; Vora, R.H. Gas-sorption properties of 6FDA–durene/1,4-phenylenediamine (pPDA) and 6FDA–durene/1,3-phenylenediamine (mPDA) copolyimides. J. Appl. Polym. Sci. 2003, 90, 2187–2193. [Google Scholar] [CrossRef]
  73. Han, S.H.; Misdan, N.; Kim, S.; Doherty, C.M.; Hill, A.J.; Lee, Y.M. Thermally Rearranged (TR) Polybenzoxazole: Effects of Diverse Imidization Routes on Physical Properties and Gas Transport Behaviors. Macromolecules 2010, 43, 7657–7667. [Google Scholar] [CrossRef]
  74. Boyd, R.H.; Pant, P.V.K. Simulation of glassy polymethylene starting from the equilibrated liquid. Macromolecules 1991, 24, 4078–4083. [Google Scholar] [CrossRef]
  75. Pant, P.V.K.; Boyd, R.H. Molecular-dynamics simulation of diffusion of small penetrants in polymers. Macromolecules 1993, 26, 679–686. [Google Scholar] [CrossRef]
  76. Van Krevelen, D.W. Properties of Polymers; Elsevier Science: Amsterdam, The Netherlands, 1997. [Google Scholar]
  77. Bondi, A. van der Waals Volumes and Radii. J. Phys. Chem. 1964, 68, 441–451. [Google Scholar] [CrossRef]
  78. Vrabec, J.; Stoll, J.; Hasse, H. A Set of Molecular Models for Symmetric Quadrupolar Fluids. J. Phys. Chem. B 2001, 105, 12126–12133. [Google Scholar] [CrossRef]
  79. Yin, D.; MacKerell, A.D. Combined ab initio/empirical approach for optimization of Lennard–Jones parameters. J. Comput. Chem. 1998, 19, 334–348. [Google Scholar] [CrossRef]
  80. Lemmon, E.W.; McLinden, M.O.; Friend, D.G. Thermophysical Properties of Fluid Systems. In NIST Chemistry WebBook, NIST Standard Reference Database Number 69; Linstrom, P.J., Mallard, W.G., Eds.; U.S Department of Commerce: Washington, DC, USA, 2021. [Google Scholar]
  81. Frenkel, D.; Smit, B. Understanding Molecular Simulation, Second Edition: From Algorithms to Applications (Computational Science); Academic Press: Cambridge, MA, USA, 2001. [Google Scholar]
  82. Van der Vegt, N.F.A.; Briels, W.J.; Wessling, M.; Strathmann, H. The sorption induced glass transition in amorphous glassy polymers. J. Chem. Phys. 1999, 110, 11061–11069. [Google Scholar] [CrossRef]
  83. Ben-Naim, A.; Marcus, Y. Solvation thermodynamics of nonionic solutes. J. Chem. Phys. 1984, 81, 2016–2027. [Google Scholar] [CrossRef]
  84. Cuthbert, T.R.; Wagner, N.J.; Paulaitis, M.E. Molecular Simulation of Glassy Polystyrene:  Size Effects on Gas Solubilities. Macromolecules 1997, 30, 3058–3065. [Google Scholar] [CrossRef]
  85. Dömötör, G.; Hentschke, R. Atomistically Modeling the Chemical Potential of Small Molecules in Dense Systems. J. Phys. Chem. B 2004, 108, 2413–2417. [Google Scholar] [CrossRef]
  86. Tamai, Y.; Tanaka, H.; Nakanishi, K. Molecular Simulation of Permeation of Small Penetrants through Membranes. 2. Solubilities. Macromolecules 1995, 28, 2544–2554. [Google Scholar] [CrossRef]
  87. Humphrey, W.; Dalke, A.; Schulten, K. VMD—Visual Molecular Dynamics. J. Molec. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  88. Böhning, M.; Springer, J. Sorptive dilation and relaxational processes in glassy polymer/gas systems—I. Poly(sulfone) and poly(ether sulfone). Polymer 1998, 39, 5183–5195. [Google Scholar] [CrossRef]
  89. Ismail, A.F.; Lorna, W. Penetrant-induced plasticization phenomenon in glassy polymers for gas separation membrane. Sep. Purif. Technol. 2002, 27, 173–194. [Google Scholar] [CrossRef]
  90. Pandiyan, S.; Brown, D.; Neyertz, S.; van der Vegt, N.F.A. Carbon Dioxide Solubility in Three Fluorinated Polyimides Studied by Molecular Dynamics Simulations. Macromolecules 2010, 43, 2605–2621. [Google Scholar] [CrossRef]
  91. Kirchheim, R. Sorption and partial molar volume of small molecules in glassy polymers. Macromolecules 1992, 25, 6952–6960. [Google Scholar] [CrossRef]
  92. Neyertz, S.; Brown, D. The effect of structural isomerism on carbon dioxide sorption and plasticization at the interface of a glassy polymer membrane. J. Membr. Sci. 2014, 460, 213–228. [Google Scholar] [CrossRef]
  93. Neyertz, S.; Brown, D.; Pandiyan, S.; van der Vegt, N.F.A. Carbon Dioxide Diffusion and Plasticization in Fluorinated Polyimides. Macromolecules 2010, 43, 7813–7827. [Google Scholar] [CrossRef]
  94. Pope, D.S.; Fleming, G.K.; Koros, W.J. Effect of various exposure histories on sorption and dilation in a family of polycarbonates. Macromolecules 1990, 23, 2988–2994. [Google Scholar] [CrossRef]
  95. De Angelis, M.G.; Merkel, T.C.; Bondar, V.I.; Freeman, B.D.; Doghieri, F.; Sarti, G.C. Gas Sorption and Dilation in Poly(2,2-bistrifluoromethyl-4,5-difluoro-1,3-dioxole-co-tetrafluoroethylene):  Comparison of Experimental Data with Predictions of the Nonequilibrium Lattice Fluid Model. Macromolecules 2002, 35, 1276–1288. [Google Scholar] [CrossRef]
  96. Neyertz, S.; Brown, D. Single- and mixed-gas sorption in large-scale molecular models of glassy bulk polymers. Competitive sorption of a binary CH4/N2 and a ternary CH4/N2/CO2 mixture in a polyimide membrane. J. Membr. Sci. 2020, 614, 118478. [Google Scholar] [CrossRef]
  97. Koros, W.J.; Chern, R.T.; Stannett, V.; Hopfenberg, H.B. A model for permeation of mixed gases and vapors in glassy polymers. J. Polym. Sci. Polym. Phys. Ed. 1981, 19, 1513–1530. [Google Scholar] [CrossRef]
  98. Koros, W.J. Model for sorption of mixed gases in glassy polymers. J. Polym. Sci. Polym. Phys. Ed. 1980, 18, 981–992. [Google Scholar] [CrossRef]
  99. Antonson, C.R.; Gardner, R.J.; King, C.F.; Ko, D.Y. Analysis of Gas Separation by Permeation in Hollow Fibers. Ind. Eng. Chem. Process. 1977, 16, 463–469. [Google Scholar] [CrossRef]
Figure 1. Schematic representation of the chemical structures of (a) the 6FDA-BAPT homopolyimide with l = 45 monomers and (b) the block 6FDA-mPDA/BAPT copolyimide with l = 29 monomers of 6FDA-mPDA followed by m = 17 monomers of 6FDA-BAPT.
Figure 1. Schematic representation of the chemical structures of (a) the 6FDA-BAPT homopolyimide with l = 45 monomers and (b) the block 6FDA-mPDA/BAPT copolyimide with l = 29 monomers of 6FDA-mPDA followed by m = 17 monomers of 6FDA-BAPT.
Polymers 15 03811 g001
Figure 2. The mean percentage of probe-accessible volume (%PAV) as a function of the probe radius used in the trial insertions.
Figure 2. The mean percentage of probe-accessible volume (%PAV) as a function of the probe radius used in the trial insertions.
Polymers 15 03811 g002
Figure 3. Visual characterisation of the free-volume distribution in the 6FDA-BAPT homopolyimide (left) and in the 6FDA-mPDA/BAPT copolyimide (right) obtained after performing trial insertions with a probe of radius 1.9 Å.
Figure 3. Visual characterisation of the free-volume distribution in the 6FDA-BAPT homopolyimide (left) and in the 6FDA-mPDA/BAPT copolyimide (right) obtained after performing trial insertions with a probe of radius 1.9 Å.
Polymers 15 03811 g003
Figure 4. Close-up view of the pure 6FDA-BAPT at 338.5 K (left), and the 6FDA-BAPT mixed-gas CH4/N2 system at 338.5 K and p = 65.5 bar 6FDA-BAPT system at 338.5 K (right). For clarity, all hydrogens in the polymer are omitted. The polymer is presented using bonds with the following colour codes: C—cyan, F—lime green, O—red, and N—blue. The penetrant molecules are presented in VDW format. For methane: C—purple and H—white. For nitrogen: N—yellow. Image created using version 1.9.4 of VMD [87].
Figure 4. Close-up view of the pure 6FDA-BAPT at 338.5 K (left), and the 6FDA-BAPT mixed-gas CH4/N2 system at 338.5 K and p = 65.5 bar 6FDA-BAPT system at 338.5 K (right). For clarity, all hydrogens in the polymer are omitted. The polymer is presented using bonds with the following colour codes: C—cyan, F—lime green, O—red, and N—blue. The penetrant molecules are presented in VDW format. For methane: C—purple and H—white. For nitrogen: N—yellow. Image created using version 1.9.4 of VMD [87].
Polymers 15 03811 g004
Figure 5. The nominal solubility coefficients for N2 (coloured symbols) and CH4 (white symbols) in 6FDA-BAPT (blue) and 6FDA-mPDA/BAPT (red) at 338.5 K as a function of the reservoir pressure, p. The smooth curves are fit to the dual-mode sorption model.
Figure 5. The nominal solubility coefficients for N2 (coloured symbols) and CH4 (white symbols) in 6FDA-BAPT (blue) and 6FDA-mPDA/BAPT (red) at 338.5 K as a function of the reservoir pressure, p. The smooth curves are fit to the dual-mode sorption model.
Polymers 15 03811 g005
Figure 6. Percentage of significant accessible volume (%SAV) for the insertion of N2 (coloured symbols) and CH4 (white symbols) in 6FDA-BAPT (blue) and 6FDA-mPDA/BAPT (red) as a function of the penetrant nominal concentration, C0(p).
Figure 6. Percentage of significant accessible volume (%SAV) for the insertion of N2 (coloured symbols) and CH4 (white symbols) in 6FDA-BAPT (blue) and 6FDA-mPDA/BAPT (red) as a function of the penetrant nominal concentration, C0(p).
Polymers 15 03811 g006
Figure 7. Intermolecular radial pair distribution functions between the centre of mass of a nitrogen molecule, N, or the methane carbon, C_met, and either (left) the dianhydride aromatic carbon or the CF3 group carbon (right).
Figure 7. Intermolecular radial pair distribution functions between the centre of mass of a nitrogen molecule, N, or the methane carbon, C_met, and either (left) the dianhydride aromatic carbon or the CF3 group carbon (right).
Polymers 15 03811 g007
Figure 8. The MSD/6t curves as a function of simulation time, t, for single-gas N2, single-gas CH4, and mixed-gas N2 and CH4 at 65.5 bar and 338.5 K. The plateau values allow for the estimation of the diffusion coefficients.
Figure 8. The MSD/6t curves as a function of simulation time, t, for single-gas N2, single-gas CH4, and mixed-gas N2 and CH4 at 65.5 bar and 338.5 K. The plateau values allow for the estimation of the diffusion coefficients.
Polymers 15 03811 g008
Table 1. Optimal numbers of each type of penetrant in the polyimide matrices when in equilibrium with 2:1 binary-gas CH4/N2 mixture at 65.5 bar and 338.5 K. Theoretical estimates of the nominal concentrations in units of cm3(STP) cm−3 are also shown together with the average volume dilations.
Table 1. Optimal numbers of each type of penetrant in the polyimide matrices when in equilibrium with 2:1 binary-gas CH4/N2 mixture at 65.5 bar and 338.5 K. Theoretical estimates of the nominal concentrations in units of cm3(STP) cm−3 are also shown together with the average volume dilations.
PolyimideGasNo. of Gas MoleculesC0
(Simulation)
C0
(Equation (16))
%Volume Swelling
6FDA-BAPTCH433139.756 ± 0.01037.691.99 ± 0.04
N2718.527 ± 0.0028.54
6FDA-mPDA/BAPTCH440045.926 ± 0.00242.292.10 ± 0.07
N2869.680 ± 0.00111.81
Table 2. Average permeabilities, P; solubility coefficients, Sc; diffusion coefficients, D; and selectivities, α, for N2 and CH4 in the single- and mixed-gas simulations conducted at 338.5 K. The average feed pressure of the single-gas systems is also given. For the mixed-gas systems, the pressure is 65.5 bar. P are expressed in Barrer, Sc in cm3(STP)/cm3(cm Hg), and D in cm2 s−1.
Table 2. Average permeabilities, P; solubility coefficients, Sc; diffusion coefficients, D; and selectivities, α, for N2 and CH4 in the single- and mixed-gas simulations conducted at 338.5 K. The average feed pressure of the single-gas systems is also given. For the mixed-gas systems, the pressure is 65.5 bar. P are expressed in Barrer, Sc in cm3(STP)/cm3(cm Hg), and D in cm2 s−1.
Polyimide6FDA-BAPT6FDA-mPDA/BAPT
p/bar66.063.3
N2Single gasP59 ± 7122 ± 10
S c 0.00620.0081
D9 × 10−715 × 10−7
Mixed gasP50 ± 6107 ± 7
S c 0.00500.0057
D10 × 10−719 × 10−7
p/bar65.065.0
CH4Single gasP44 ± 11140 ± 14
S c 0.01080.0120
D4 × 10−712 × 10−7
Mixed gasP74 ± 13146 ± 15
S c 0.01260.0134
D6 × 10−711 × 10−7
α(N2/CH4)Single gas 1.3 ± 0.5 0.9 ± 0.2
Mixed gas 0.7 ± 0.20.7 ± 0.1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tanis, I.; Brown, D.; Neyertz, S.; Vaidya, M.; Ballaguet, J.-P.; Duval, S.; Bahamdan, A. A Molecular Dynamics Study of Single-Gas and Mixed-Gas N2 and CH4 Transport in Triptycene-Based Polyimide Membranes. Polymers 2023, 15, 3811. https://doi.org/10.3390/polym15183811

AMA Style

Tanis I, Brown D, Neyertz S, Vaidya M, Ballaguet J-P, Duval S, Bahamdan A. A Molecular Dynamics Study of Single-Gas and Mixed-Gas N2 and CH4 Transport in Triptycene-Based Polyimide Membranes. Polymers. 2023; 15(18):3811. https://doi.org/10.3390/polym15183811

Chicago/Turabian Style

Tanis, Ioannis, David Brown, Sylvie Neyertz, Milind Vaidya, Jean-Pierre Ballaguet, Sebastien Duval, and Ahmad Bahamdan. 2023. "A Molecular Dynamics Study of Single-Gas and Mixed-Gas N2 and CH4 Transport in Triptycene-Based Polyimide Membranes" Polymers 15, no. 18: 3811. https://doi.org/10.3390/polym15183811

APA Style

Tanis, I., Brown, D., Neyertz, S., Vaidya, M., Ballaguet, J. -P., Duval, S., & Bahamdan, A. (2023). A Molecular Dynamics Study of Single-Gas and Mixed-Gas N2 and CH4 Transport in Triptycene-Based Polyimide Membranes. Polymers, 15(18), 3811. https://doi.org/10.3390/polym15183811

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop