Molecular dynamics simulations of gas separations using faujasite-type zeolite membranes

Gas separations with faujasite zeolite membranes have been examined using the method of molecular dynamics. Two binary mixtures are investigated, oxygen/nitrogen and nitrogen/carbon dioxide. These mixtures have been found experimentally to exhibit contrasting behavior. In O(2)/N(2) mixtures the ideal selectivity (pure systems) is higher than the mixture selectivity, while in N(2)/CO(2) the mixture selectivity is higher than the ideal selectivity. One of the key goals of this work was to seek a fundamental molecular level understanding of such divergent behavior. Our simulation results (using previously developed intermolecular models for both the gases and zeolites investigated) were found to replicate this experimental behavior. By examining the loading of the membranes and the diffusion rates inside the zeolites, we have been able to explain such contrasting behavior of O(2)/N(2) and N(2)/CO(2) mixtures. In the case of O(2)/N(2) mixtures, the adsorption and loading of both O(2) and N(2) in the membrane are quite competitive, and thus the drop in the selectivity in the mixture is primarily the result of oxygen slowing the diffusion of nitrogen and nitrogen somewhat increasing the diffusion of oxygen when they pass through the zeolite pores. In N(2)/CO(2) systems, CO(2) is rather selectively adsorbed and loaded in the zeolite, leaving very little room for N(2) adsorption. Thus although N(2) continues to have a higher diffusion rate than CO(2) even in the mixture, there are so few N(2) molecules in the zeolite in mixtures that the selectivity of the mixture increases significantly compared to the ideal (pure system) values. We have also compared simulation results with hydrodynamic theories that classify the permeance of membranes to be either due to surface diffusion, viscous flow, or Knudsen diffusion. Our results show surface diffusion to be the dominant mode, except in the case of N(2)/CO(2) binary mixtures where Knudsen diffusion also makes a contribution to N(2) transport.     

In silico screening of metal-organic frameworks in separation applications

Porous materials such as metal-organic frameworks (MOFs) and zeolitic imidazolate frameworks (ZIFs) offer considerable potential for separating a variety of mixtures such as those relevant for CO(2) capture (CO(2)/H(2), CO(2)/CH(4), CO(2)/N(2)), CH(4)/H(2), alkanes/alkenes, and hydrocarbon isomers. There are basically two different separation technologies that can be employed: (1) a pressure swing adsorption (PSA) unit with a fixed bed of adsorbent particles, and (2) a membrane device, wherein the mixture is allowed to permeate through a micro-porous crystalline layer. In view of the vast number of MOFs, and ZIFs that have been synthesized there is a need for a systematic screening of potential candidates for any given separation task. Also of importance is to investigate how MOFs and ZIFs stack up against the more traditional zeolites such as NaX and NaY with regard to their separation characteristics. This perspective highlights the potency of molecular simulations in determining the choice of the best MOF or ZIF for a given separation task. A variety of metrics that quantify the separation performance, such as adsorption selectivity, working capacity, diffusion selectivity, and membrane permeability, are determined from a combination of Configurational-Bias Monte Carlo (CBMC) and Molecular Dynamics (MD) simulations. The practical utility of the suggested screening methodology is demonstrated by comparison with available experimental data.     

Structure-transport correlation for the diffusive tortuosity of bulk, monodisperse, random sphere packings

The mass transport properties of bulk random sphere packings depend primarily on the bed (external) porosity ε, but also on the packing microstructure. We investigate the influence of the packing microstructure on the diffusive tortuosity τ=D(m)/D(eff), which relates the bulk diffusion coefficient (D(m)) to the effective (asymptotic) diffusion coefficient in a porous medium (D(eff)), by numerical simulations of diffusion in a set of computer-generated, monodisperse, hard-sphere packings. Variation of packing generation algorithm and protocol yielded four Jodrey-Tory and two Monte Carlo packing types with systematically varied degrees of microstructural heterogeneity in the range between the random-close and the random-loose packing limit (ε=0.366-0.46). The distinctive tortuosity-porosity scaling of the packing types is influenced by the extent to which the structural environment of individual pores varies in a packing, and to quantify this influence we propose a measure based on Delaunay tessellation. We demonstrate that the ratio of the minimum to the maximum void face area of a Delaunay tetrahedron around a pore between four adjacent spheres, (A(min)/A(max))(D), is a measure for the structural heterogeneity in the direct environment of this pore, and that the standard deviation σ of the (A(min)/A(max))(D)-distribution considering all pores in a packing mimics the tortuosity-porosity scaling of the generated packing types. Thus, σ(A(min)/A(max))(D) provides a structure-transport correlation for diffusion in bulk, monodisperse, random sphere packings.     

Highly hydrothermally stable microporous silica membranes for hydrogen separation

Fluorocarbon-modified silica membranes were deposited on gamma-Al2O3/alpha-Al2O3 supports by the sol-gel technique for hydrogen separation. The hydrophobic property, pore structure, gas transport and separation performance, and hydrothermal stability of the modified membranes were investigated. It is observed that the water contact angle increases from 27.2+/-1.5 degrees for the pure silica membranes to 115.0+/-1.2 degrees for the modified ones with a (trifluoropropyl)triethoxysilane (TFPTES)/tetraethyl orthosilicate (TEOS) molar ratio of 0.6. The modified membranes preserve a microporous structure with a micropore volume of 0.14 cm3/g and a pore size of approximately 0.5 nm. A single gas permeation of H2 and CO2 through the modified membranes presents small positive apparent thermal activation energies, indicating a dominant microporous membrane transport. At 200 degrees C, a single H2 permeance of 3.1x10(-6) mol m(-2) s(-1) Pa(-1) and a H2/CO2 permselectivity of 15.2 were obtained after proper correction for the support resistance and the contribution from the defects. In the gas mixture measurement, the H2 permeance and the H2/CO2 separation factor almost remain constant at 200 degrees C with a water vapor pressure of 1.2x10(4) Pa for at least 220 h, indicating that the modified membranes are hydrothermally stable, benefiting from the integrity of the microporous structure due to the fluorocarbon modification.     

The permeation of binary gas mixtures through support structures of composite membranes

The dusty gas model (DGM) is used to describe transport of binary gas mixtures through porous membrane supports to quantify the resistance towards permeation. The model equations account for three different transport mechanisms for the permeating components: conventional viscous pore flow, Knudsen diffusion, and binary diffusion. Experimental data obtained with the uncoated membrane supports are used to determine the morphological parameters needed in the DGM equations. Flat sheet and hollow fiber membrane supports are characterized by the permeation of a TCE/nitrogen vapor. The DGM shows an excellent fit to experimental data when the asymmetric structure of the membrane supports is taken into account, but the morphological parameters cannot necessarily be related to precise physical structure parameters such as pore size, porosity, and tortuosity. The DGM works well even when the membrane supports are modeled as a single homogenous structure. The membrane supports exhibit different resistances towards the various transport mechanisms that occur within the porous support and the resistances vary with process conditions so that support optimization is not straightforward. With the analysis presented in this paper and transport equations specific to the dense coating and module geometries, the influence of the support layer on gas or vapor separation can be quantified.     

Combination of MOFs and zeolites for mixed-matrix membranes

Mixed-matrix membranes (MMMs) were prepared by combinations of two different kinds of porous fillers [metal-organic frameworks (MOFs) HKUST-1 and ZIF-8, and zeolite silicalite-1] and polysulfone. In the search for filler synergy, the MMMs were applied to the separation of CO(2)/N(2), CO(2)/CH(4), O(2)/N(2), and H(2)/CH(4) mixtures and we found important selectivity improvements with the HKUST-1-silicalite-1 system (CO(2)/CH(4) and CO(2)/N(2) separation factors of 22.4 and 38.0 with CO(2) permeabilities of 8.9 and 8.4 Barrer, respectively).     

Influence of MgO template on carbon dioxide adsorption of cation exchange resin-based nanoporous carbon

In this work, porous carbons with well-developed pore structures were directly prepared from a weak acid cation exchange resin (CER) by the carbonization of a mixture with Mg acetate in different ratios. The effect of the Mg acetate-to-CER ratio on the pore structure and CO(2) adsorption capacities of the obtained porous carbons was studied. The textural properties and morphologies of the porous carbons were analyzed via N(2)/77K adsorption/desorption isotherms, SEM, and TEM, respectively. The CO(2) adsorption capacities of the prepared porous carbons were measured at 298 K and 1 bar and 30 bar. By dissolving the MgO template, the porous carbons exhibited high specific surface areas (326-1276 m(2)/g) and high pore volumes (0.258-0.687 cm(3)/g). The CO(2) adsorption capacities of the porous carbons were enhanced to 164.4 mg/g at 1 bar and 1045 mg/g at 30 bar by increasing the Mg acetate-to-CER ratio. This result indicates that CER was one of the carbon precursors to producing the porous structure, as well as for improving the CO(2) adsorption capacities of the carbon species.     

Multicomponent ballistic transport in narrow single wall carbon nanotubes: analytic model and molecular dynamics simulations

The transport of gas mixtures through molecular-sieve membranes such as narrow nanotubes has many potential applications, but there remain open questions and a paucity of quantitative predictions. Our model, based on extensive molecular dynamics simulations, proposes that ballistic motion, hindered by counter diffusion, is the dominant mechanism. Our simulations of transport of mixtures of molecules between control volumes at both ends of nanotubes give quantitative support to the model's predictions. The combination of simulation and model enable extrapolation to longer tubes and pore networks.     

Computational Study of Adsorption and Separation of CO(2), CH(4), and N(2) by an rht-Type Metal-Organic Framework

In this work, a computational study is performed to evaluate the adsorption-based separation of CO(2) from flue gas (mixtures of CO(2) and N(2)) and natural gas (mixtures of CO(2) and CH(4)) using microporous metal organic framework Cu-TDPAT as a sorbent material. The results show that electrostatic interactions can greatly enhance the separation efficiency of this MOF for gas mixtures of different components. Furthermore, the study also suggests that Cu-TDPAT is a promising material for the separation of CO(2) from N(2) and CH(4), and its macroscopic separation behavior can be elucidated on a molecular level to give insight into the underlying mechanisms. On the basis of the single-component CO(2), N(2), and CH(4) isotherms, binary mixture adsorption (CO(2)/N(2) and CO(2)/CH(4)) and ternary mixture adsorption (CO(2)/N(2)/CH(4)) were predicted using the ideal adsorbed solution theory (IAST). The effect of H(2)O vapor on the CO(2) adsorption selectivity and capacity was also examined. The applicability of IAST to this system was validated by performing GCMC simulations for both single-component and mixture adsorption processes.     

Dual control cell reaction ensemble molecular dynamics: a method for simulations of reactions and adsorption in porous materials

We present a simulation tool to study fluid mixtures that are simultaneously chemically reacting and adsorbing in a porous material. The method is a combination of the reaction ensemble Monte Carlo method and the dual control volume grand canonical molecular dynamics technique. The method, termed the dual control cell reaction ensemble molecular dynamics method, allows for the calculation of both equilibrium and nonequilibrium transport properties in porous materials such as diffusion coefficients, permeability, and mass flux. Control cells, which are in direct physical contact with the porous solid, are used to maintain the desired reaction and flow conditions for the system. The simulation setup closely mimics an actual experimental system in which the thermodynamic and flow parameters are precisely controlled. We present an application of the method to the dry reforming of methane reaction within a nanoscale reactor model in the presence of a semipermeable membrane that was modeled as a porous material similar to silicalite. We studied the effects of the membrane structure and porosity on the reaction species permeability by considering three different membrane models. We also studied the effects of an imposed pressure gradient across the membrane on the mass flux of the reaction species. Conversion of syngas (H2/CO) increased significantly in all the nanoscale membrane reactor models considered. A brief discussion of further potential applications is also presented.     

Theoretical and experimental correlations of gas dissolution, diffusion, and thermodynamic properties in determination of gas permeability and selectivity in supported ionic liquid membranes

Supported ionic liquid membranes (SILMs) has the potential to be a new technological platform for gas/organic vapour separation because of the unique non-volatile nature and discriminating gas dissolution properties of room temperature ionic liquids (ILs). This work starts with an examination of gas dissolution and transport properties in bulk imidazulium cation based ionic liquids [C(n)mim][NTf2] (n=2.4, 6, 8.10) from simple gas H(2), N(2), to polar CO(2), and C(2)H(6), leading to a further analysis of how gas dissolution and diffusion are influenced by molecular specific gas-SILMs interactions, reflected by differences in gas dissolution enthalpy and entropy. These effects were elucidated again during gas permeation studies by examining how changes in these properties and molecular specific interactions work together to cause deviations from conventional solution-diffusion theory and their impact on some remarkably contrasting gas perm-selectivity performance. The experimental perm-selectivity for all tested gases showed varied and contrasting deviation from the solution-diffusion, depending on specific gas-IL combinations. It transpires permeation for simpler non-polar gases (H(2), N(2)) is diffusion controlled, but strong molecular specific gas-ILs interactions led to a different permeation and selectivity performance for C(2)H(6) and CO(2). With exothermic dissolution enthalpy and large order disruptive entropy, C(2)H(6) displayed the fastest permeation rate at increased gas phase pressure in spite of its smallest diffusivity among the tested gases. The C(2)H(6) gas molecules "peg" on the side alkyl chain on the imidazulium cation at low concentration, and are well dispersed in the ionic liquids phase at high concentration. On the other hand strong CO(2)-ILs affinity resulted in a more prolonged "residence time" for the gas molecule, typified by reversed CO(2)/N(2) selectivity and slowest CO(2) transport despite CO(2) possess the highest solubility and comparable diffusivity in the ionic liquids. The unique transport and dissolution behaviour of CO(2) are further exploited by examining the residing state of CO(2) molecules in the ionic liquid phase, which leads to a hypothesis of a condensing and holding capacity of ILs towards CO(2), which provide an explanation to slower CO(2) transport through the SILMs. The pressure related exponential increase in permeations rate is also analysed which suggests a typical concentration dependent diffusion rate at high gas concentration under increased gas feed pressure. Finally the strong influence of discriminating and molecular specific gas-ILs interactions on gas perm-selectivity performance points to future specific design of ionic liquids for targeted gas separations.     

Molecular dynamics simulations of transport and separation of carbon dioxide-alkane mixtures in carbon nanopores

The configurational-bias Monte Carlo method, which is used for efficient generation of molecular models of n-alkane chains, is combined for the first time with the dual control-volume grand-canonical molecular-dynamics simulation, which has been developed for studying transport of molecules in pores under an external potential gradient, to investigate transport and separation of binary mixtures of n-alkanes, as well as mixtures of CO2 and n-alkanes, in carbon nanopores. The effect of various factors, such as the temperature of the system, the composition of the mixture, and the pore size, on the separation of the mixtures is investigated. We also report the preliminary results of an experimental study of transport and separation of some of the same mixtures in a carbon molecular-sieve membrane with comparable pore sizes. The results indicate that, for the mixtures considered in this paper, even in very small carbon nanopores the energetic effects still play a dominant role in the transport and separation properties of the mixtures, whereas in a real membrane they are dominated by the membrane's morphological characteristics. As a result, for the mixtures considered, a single pore may be a grossly inadequate model of a real membrane, and hence one must resort to three-dimensional molecular pore network models of the membrane.     

Combinatorial synthesis and characterization of metal-open frameworks in mild and friendly conditions: application to CO₂ adsorption

Combinatorial screening using precipitation methods at room temperature can lead to a great diversity of carboxylate based Metal Organic Frameworks (MOFs) including already known or original porous solids. The investigation of the synthesis of MOFs in different solvent and solvent mixtures includes the use of solvents such as alcohols and tetrahydrofuran (THF) which would greatly facilitate large scale production. We also show the application of Principal Component Analysis (PCA) and clustering techniques on large libraries of XRD diffraction files in order to identify classes of similar phases and peculiar phases. The combinatorial screening of 105 samples in the La/btc system has led to the identification of two phases which are solvent depending. On the La(btc) compound, the CO? adsorption measurements reveal a guest-host interactions as supported by XRD phase transformation upon thermal treatment. The mass transport can be assigned to a "single file diffusion" regime due to the one dimensional channel porous structure associated to small pore size.     

Modeling growth of organized nanoporous structures by anodic oxidation

Nanostructured porous oxides are produced by anodic dissolution of several metals. A scaling approach is introduced to explain pattern nucleation in an oxide layer, and a related microscopic model shows oxide growth with long nanopores. The scaling approach matches the time of ion transport across the thin oxide layer, which is related to metal corrosion, and the time of diffusion along the oxide/solution (OS) interface, which represents the extension of oxide dissolution. The selected pattern size is of order (dD(S)/v(O))(1/2), where d is the oxide thickness, v(O) is the migration velocity of oxygen ions across the oxide, and D(s) is the diffusion coefficient of H(+) ions along the oxide/solution interface. This result is consistent with available experimental data for those quantities, predicts the increase of pore size with the external voltage, and suggests the independence of pore size with the solution pH. Subsequently, we propose a microscopic model that expresses the main physicochemical processes as a set of characteristic lengths for diffusion and surface relaxation. It shows a randomly perturbed OS interface at short times, its evolution to pore nucleation and to stable growth of very long pores, in agreement with the mechanistic scenario suggested by two experimental groups. The decrease of the size of the walls between the pores with the interface tension is consistent with arguments for formation of titania nanotube arrays instead of nanopores. These models show that pattern nucleation and growth depend on matching a small number of physicochemical parameters, which is probably the reason for the production of nanostructured porous oxides from various materials under suitable electrochemical conditions.     

The effect of obstacle conductivity and electric field on effective mobility and dispersion in electrophoretic transport: a volume averaging approach

The method of volume averaging has been used to determine the effective electrophoretic mobility and dispersion coefficients for molecular transport of point-like solutes in a two-phase porous medium where the electrical conductivity and the diffusion and mobility coefficients may vary in both phases. The formal theory, derived in previous work, is numerically evaluated for cases where the obstacle phase has a large or small conductivity relative to the fluid phase and where the diffusion coefficient of the solute in the obstacle phase can be large or small relative to that in the fluid phase. In agreement with previous Monte Carlo methods, the effective electrophoretic mobility is not a function of media conductivity or electric field when the obstacles are impermeable to solute transport or have small diffusion solute diffusion coefficients. However, the dispersion coefficient is a strong function of electric field and varies with obstacle conductivity when diffusive transport is small in the obstacles relative to the fluid. In contrast, the effective electrophoretic mobility is a function of electric field when conductivity of the obstacles is much larger than the fluid and when the obstacles are very permeable to solute but have low electrical conductivity.     

Sulfonate-grafted porous polymer networks for preferential CO2 adsorption at low pressure

A porous polymer network (PPN) grafted with sulfonic acid (PPN-6-SO(3)H) and its lithium salt (PPN-6-SO(3)Li) exhibit significant increases in isosteric heats of CO(2) adsorption and CO(2)-uptake capacities. IAST calculations using single-component-isotherm data and a 15/85 CO(2)/N(2) ratio at 295 K and 1 bar revealed that the sulfonate-grafted PPN-6 networks show exceptionally high adsorption selectivity for CO(2) over N(2) (155 and 414 for PPN-6-SO(3)H and PPN-6-SO(3)Li, respectively). Since these PPNs also possess ultrahigh physicochemical stability, practical applications in postcombustion capture of CO(2) lie well within the realm of possibility.     

Novel general expressions that describe the behavior of the height equivalent of a theoretical plate in chromatographic systems involving electrically-driven and pressure-driven flows

Novel general expressions are constructed and presented that describe the behavior of the height equivalent of a theoretical plate (plate height), H, as a function of the linear velocity, Vx, along the axis, x, of the column and the kinetic parameters that characterize the mass transfer and adsorption mechanisms in chromatographic columns. Open tube capillaries as well as columns packed with either non-porous or porous particles are studied. The porous particles could have unimodal or bimodal pore-size distributions and intraparticle convective fluid flow and pore diffusion are considered. The expressions for the plate height, H, presented in this work could be applicable to high-performance liquid chromatography (HPLC) and capillary electrochromatography (CEC) systems, and could be used together with experimental plate height, H, versus linear velocity, Vx, data to determine the values of the parameters that characterize intraparticle convective fluid flow and pore diffusion. Furthermore, chromatographic systems under unretained as well as under retained conditions are examined. The experimental values of the plate height, H, versus the linear velocity, Vx, for a CEC system involving charged porous silica C8 particles and an uncharged analyte are compared with the theoretical results for the plate height, H, obtained from the expressions presented in this work. The agreement between theory and experiment is good, and the results indicate that the magnitude of the intraparticle electroosmotic flow (EOF) in the pores of the particles is substantial while the pore diffusion coefficient was of small magnitude. But the overall intraparticle mass transfer resistance in these particles was low because of the significant contribution of the intraparticle EOF. Simulation results are also presented (i) for a hybrid HPLC-CEC system, and (ii) for different CEC systems involving open capillaries as well as packed columns having non-porous or porous particles. The analysis of the results indicates (a) the reasons for the superior performance exhibited by the hybrid HPLC-CEC system over the performance obtained when the system is operated only in the HPLC mode, and (b) the operational configuration and the properties that the structure of the porous particles would have to have in CEC systems involving uncharged or charged analytes under unretained or retained conditions in order to obtain high CEC efficiency (low values of the plate height, H).