Quantum effect induced reverse kinetic molecular sieving in microporous materials

We report kinetic molecular sieving of hydrogen and deuterium in zeolite rho at low temperatures, using atomistic molecular dynamics simulations incorporating quantum effects via the Feynman-Hibbs approach. We find that diffusivities of confined molecules decrease when quantum effects are considered, in contrast with bulk fluids which show an increase. Indeed, at low temperatures, a reverse kinetic sieving effect is demonstrated in which the heavier isotope, deuterium, diffuses faster than hydrogen. At 65 K, the flux selectivity is as high as 46, indicating a good potential for isotope separation.

Quantum effects on adsorption and diffusion of hydrogen and deuterium in microporous materials

Monte Carlo and molecular dynamics simulations and neutron scattering experiments are used to study the adsorption and diffusion of hydrogen and deuterium in zeolite Rho in the temperature range of 30-150 K. In the molecular simulations, quantum effects are incorporated via the Feynman-Hibbs variational approach. We suggest a new set of potential parameters for hydrogen, which can be used when Feynman-Hibbs variational approach is used for quantum corrections. The dynamic properties obtained from molecular dynamics simulations are in excellent agreement with the experimental results and show significant quantum effects on the transport at very low temperature. The molecular dynamics simulation results show that the quantum effect is very sensitive to pore dimensions and under suitable conditions can lead to a reverse kinetic molecular sieving with deuterium diffusing faster than hydrogen.

Structural evolution in a hydrothermal reaction between Nb2O5 and NaOH solution: From Nb2O5 grains to microporous Na2Nb2O6 center dot(2)/3H2O fibers and NaNbO3 cubes

Niobium pentoxide reacts actively with concentrate NaOH solution under hydrothermal conditions at as low as 120 degrees C. The reaction ruptures the corner-sharing of NbO7 decahedra and NbO6 octahedra in the reactant Nb2O5, yielding various niobates, and the structure and composition of the niobates depend on the reaction temperature and time. The morphological evolution of the solid products in the reaction at 180 degrees C is monitored via SEM: the fine Nb2O5 powder aggregates first to irregular bars, and then niobate fibers with an aspect ratio of hundreds form. The fibers are microporous molecular sieve with a monoclinic lattice, Na2Nb2O6 center dot(2)/3H2O. The fibers are a metastable intermediate of this reaction, and they completely convert to the final product NaNbO3 Cubes in the prolonged reaction of 1 h. This study demonstrates that by carefully optimizing the reaction condition, we can selectively fabricate niobate structures of high purity, including the delicate microporous fibers, through a direct reaction between concentrated NaOH solution and Nb2O5. This synthesis route is simple and suitable for the large-scale production of the fibers. The reaction first yields poorly crystallized niobates consisting of edge-sharing NbO6 octahedra, and then the microporous fibers crystallize and grow by assembling NbO6 octahedra or clusters of NbO6 octahedra and NaO6 units. Thus, the selection of the fibril or cubic product is achieved by control of reaction kinetics. Finally, niobates with different structures exhibit remarkable differences in light absorption and photoluminescence properties. Therefore, this study is of importance for developing new functional materials by the wet-chemistry process.

Systematic selection of solvents for the fabrication of 3D combined macro- and microporous polymeric scaffolds for soft tissue engineering

In this study, we investigate the fabrication of 3D porous poly(lactic-co-glycolic acid) (PLGA) scaffolds using the thermally-induced phase separation technique. The current study focuses on the selection of alternative solvents for this process using a number of criteria, including predicted solubility. toxicity, removability and processability. Solvents were removed via either vacuum freeze-drying or leaching, depending on their physical properties. The residual solvent was tested using gas chromatography-mass spectrometry. A large range of porous, highly interconnected scaffold architectures with tunable pore size and alignment was obtained, including combined macro- and microporous structures and an entirely novel 'porous-fibre' structure. The morphological features of the most promising poly(lactic-co-glycolic acid) scaffolds were analysed via scanning electron microscopy and X-ray micro-computed tomography in both two and three dimensions. The Young's moduli of the scaffolds under conditions of temperature, pH and ionic strength similar to those found in the body were tested and were found to be highly dependent on the architectures.

High-pressure adsorption of supercritical gases on activated carbons: An improved approach based on the density functional theory and the bender equation of state

Adsorption of nitrogen, argon, methane, and carbon dioxide on activated carbon Norit R1 over a wide range of pressure (up to 50 MPa) at temperatures from 298 to 343 K (supercritical conditions) is analyzed by means of the density functional theory modified by incorporating the Bender equation of state, which describes the bulk phase properties with very high accuracy. It has allowed us to precisely describe the experimental data of carbon dioxide adsorption slightly above and below its critical temperatures. The pore size distribution (PSD) obtained with supercritical gases at ambient temperatures compares reasonably well with the PSD obtained with subcritical nitrogen at 77 K. Our approach does not require the skeletal density of activated carbon from helium adsorption measurements to calculate excess adsorption. Instead, this density is treated as a fitting parameter, and in all cases its values are found to fall into a very narrow range close to 2000 kg/m(3). It was shown that in the case of high-pressure adsorption of supercritical gases the PSD could be reliably obtained for the range of pore width between 0.6 and 3 run. All wider pores can be reliably characterized only in terms of surface area as their corresponding excess local isotherms are the same over a practical range of pressure.

Modeling molecular transport in slit pores

We examine the transport of methane in microporous carbon by performing equilibrium and nonequilibrium molecular dynamics simulations over a range of pore sizes, densities, and temperatures. We interpret these simulation results using two models of the transport process. At low densities, we consider a molecular flow model, in which intermolecular interactions are neglected, and find excellent agreement between transport diffusion coefficients determined from simulation, and those predicted by the model. Simulation results indicate that the model can be applied up to fluid densities of the order to 0.1-1 nm(-3). Above these densities, we consider a slip flow model, combining hydrodynamic theory with a slip condition at the solid-fluid interface. As the diffusion coefficient at low densities can be accurately determined by the molecular flow model, we also consider a model where the slip condition is supplied by the molecular flow model. We find that both density-dependent models provide a useful means of estimating the transport coefficient that compares well with simulation. (C) 2004 American Institute of Physics.

Permeation of binary gas mixtures in ultramicroporous membranes

High-quality nanometer thick ultramicroporous membranes were prepared from silica sol-gel processes and tested for the permeation of binary gas mixtures of He, H-2, CO2, and CH4 across different temperature and partial pressure regimens. Pore size distribution by molecular probing showed that the majority of pore sizes had dimensions below 2.9 Angstrom. In 50:50 binary mixtures, the fluxes of gases increased as a function of temperature, indicating an activated transport mechanism. The ultramicroporous membranes showed high selectivities at 150 degreesC for He/CO2 (30), He/CH4 (93), H-2/CO2 (10), and H-2/CH4 (9) with lower selectivities for CO2/CH4 (5). High activation energies (E-a) were observed for the permeance of 50:50 binary mixtures containing He and H-2 of 22.1-27.5 and 17.6-23.1 kJ.mol(-1), respectively. The E-a for the permeance of the total mixture approached the E-a for the permeance of the molecule with the smaller kinetic diameter (He or H-2).

Non-additivity of attractive potentials in modeling of N-2 and Ar adsorption isotherms on graphitized carbon black and porous carbon by means of density functional theory

We present a new approach accounting for the nonadditivity of attractive parts of solid-fluid and fluidfluid potentials to improve the quality of the description of nitrogen and argon adsorption isotherms on graphitized carbon black in the framework of non-local density functional theory. We show that the strong solid-fluid interaction in the first monolayer decreases the fluid-fluid interaction, which prevents the twodimensional phase transition to occur. This results in smoother isotherm, which agrees much better with experimental data. In the region of multi-layer coverage the conventional non-local density functional theory and grand canonical Monte Carlo simulations are known to over-predict the amount adsorbed against experimental isotherms. Accounting for the non-additivity factor decreases the solid-fluid interaction with the increase of intermolecular interactions in the dense adsorbed fluid, preventing the over-prediction of loading in the region of multi-layer adsorption. Such an improvement of the non-local density functional theory allows us to describe experimental nitrogen and argon isotherms on carbon black quite accurately with mean error of 2.5 to 5.8% instead of 17 to 26% in the conventional technique. With this approach, the local isotherms of model pores can be derived, and consequently a more reliab * le pore size distribution can be obtained. We illustrate this by applying our theory against nitrogen and argon isotherms on a number of activated carbons. The fitting between our model and the data is much better than the conventional NLDFT, suggesting the more reliable PSD obtained with our approach.

Application of density functional theory to analysis of energetic heterogeneity and pore size distribution of activated carbons

A new approach based on the nonlocal density functional theory to determine pore size distribution (PSD) of activated carbons and energetic heterogeneity of the pore wall is proposed. The energetic heterogeneity is modeled with an energy distribution function (EDF), describing the distribution of solid-fluid potential well depth (this distribution is a Dirac delta function for an energetic homogeneous surface). The approach allows simultaneous determining of the PSD (assuming slit shape) and EDF from nitrogen or argon isotherms at their respective boiling points by using a set of local isotherms calculated for a range of pore widths and solid-fluid potential well depths. It is found that the structure of the pore wall surface significantly differs from that of graphitized carbon black. This could be attributed to defects in the crystalline structure of the surface, active oxide centers, finite size of the pore walls (in either wall thickness or pore length), and so forth. Those factors depend on the precursor and the process of carbonization and activation and hence provide a fingerprint for each adsorbent. The approach allows very accurate correlation of the experimental adsorption isotherm and leads to PSDs that are simpler and more realistic than those obtained with the original nonlocal density functional theory.

Multicomponent adsorption on activated carbons under supercritical conditions

Adsorption of binary mixtures onto activated carbon Norit R1 for the system nitrogen-methane-carbon dioxide was investigated over the pressure range up to 15 MPa. A new model is proposed to describe the experimental data. It is based on the assumption that an activated carbon can be characterized by the distribution function of elements of adsorption volume (EAV) over the solid-fluid potential. This function may be evaluated from pure component isotherms using the equality of the chemical potentials in the adsorbed phase and in the bulk phase for each EAV. In the case of mixture adsorption a simple combining rule is proposed, which allows determining the adsorbed phase density and its composition in the EAV at given pressure and compositions of the bulk phase. The adsorbed concentration of each adsorbate is the integral of its density over the set of EAV. The comparison with experimental data on binary mixtures has shown that the approach works reasonably well. In the case of high-pressure binary mixture adsorption, when only total amount adsorbed was measured, the proposed model allows reliably determining partial amounts of the adsorbed components. (C) 2004 Elsevier Inc. All rights reserved.

Prediction of high-pressure adsorption equilibrium of supercritical gases using density functional theory

In this paper, we present the results of the prediction of the high-pressure adsorption equilibrium of supercritical. gases (Ar, N-2, CH4, and CO2) on various activated carbons (BPL, PCB, and Norit R1 extra) at various temperatures using a density-functional-theory-based finite wall thickness (FWT) model. Pore size distribution results of the carbons are taken from our recent previous work 1,2 using this approach for characterization. To validate the model, isotherms calculated from the density functional theory (DFT) approach are comprehensively verified against those determined by grand canonical Monte Carlo (GCMC) simulation, before the theoretical adsorption isotherms of these investigated carbons calculated by the model are compared with the experimental adsorption measurements of the carbons. We illustrate the accuracy and consistency of the FWT model for the prediction of adsorption isotherms of the all investigated gases. The pore network connectivity problem occurring in the examined carbons is also discussed, and on the basis of the success of the predictions assuming a similar pore size distribution for accessible and inaccessible regions, it is suggested that this is largely related to the disordered nature of the carbon.

Effects of the juxtaposition of carbonaceous slit pores on the overall transport behavior of adsorbed fluids

It is a common approximation in the modeling of adsorption in microporous carbons to treat the pores as slit pores, whose walls are considered to consist of an infinite number of graphitic layers. In practice, such an approximation is appropriate as long as the number of graphitic layers in the wall is greater than three. However, it is understood that pore walls in microporous carbons commonly consist of three or fewer layers. As well as affecting the solid-fluid interaction within a pore, such narrow walls permit the interaction of fluid molecules through the wall, with consequences for the adsorption characteristics. We consider the effect that a distributed pore-wall thickness model can have on transport properties. At low density we find that the only significant deviation in the transport properties from the infinite pore-wall thickness model occurs in pores with single-layer walls. For a model of activated carbons with a distribution of pore widths and pore-wall thicknesses, the transport properties are generally insensitive to the effects of finite walls, in terms of both the solid-fluid interaction within a pore and fluid-fluid interaction through the pore walls.

Grand canonical Monte Carlo simulation study of methane adsorption at an open graphite surface and in slitlike carbon pores at 273 K

Grand canonical Monte Carlo (GCMC) simulation was used for the systematic investigation of the supercritical methane adsorption at 273 K on an open graphite surface and in slitlike micropores of different sizes. For both considered adsorption systems the calculated excess adsorption isotherms exhibit a maximum. The effect of the pore size on the maximum surface excess and isosteric enthalpy of adsorption for methane storage at 273 K is discussed. The microscopic detailed picture of methane densification near the homogeneous graphite wall and in slitlike pores at 273 K is presented with selected local density profiles and snapshots. Finally, the reliable pore size distributions, obtained in the range of the microporosity, for two pitch-based microporous activated carbon fibers are calculated from the local excess adsorption isotherms obtained via the GCMC simulation. The current systematic study of supercritical methane adsorption both on an open graphite surface and in slitlike micropores performed by the GCMC summarizes recent investigations performed at slightly different temperatures and usually a lower pressure range by advanced methods based on the statistical thermodynamics.

Adsorption of pure fluids in activated carbon with a modified lattice gas model

In this article we study the effects of adsorbed phase compression, lattice structure, and pore size distribution on the analysis of adsorption in microporous activated carbon. The lattice gas approach of Ono-Kondo is modified to account for the above effects. Data of nitrogen adsorption at 77 K onto a number of activated carbon samples are analyzed to investigate the pore filling pressure versus pore width, the packing effect, and the compression of the adsorbed phase. It is found that the PSDs obtained from this analysis are comparable to those obtained by the DFT method. The discrete nature of the PSDs derived from the modified lattice gas theory is due to the inherent assumption of discrete layers of molecules. Nevertheless, it does provide interesting information on the evolution of micropores during the activation process.

Adsorption in slit pores and pore-size distribution: A molecular layer structure theory

A new approach is developed to analyze the thermodynamic properties of a sub-critical fluid adsorbed in a slit pore of activated carbon. The approach is based on a representation that an adsorbed fluid forms an ordered structure close to a smoothed solid surface. This ordered structure is modelled as a collection of parallel molecular layers. Such a structure allows us to express the Helmholtz free energy of a molecular layer as the sum of the intrinsic Helmholtz free energy specific to that layer and the potential energy of interaction of that layer with all other layers and the solid surface. The intrinsic Helmholtz free energy of a molecular layer is a function (at given temperature) of its two-dimensional density and it can be readily obtained from bulk-phase properties, while the interlayer potential energy interaction is determined by using the 10-4 Lennard-Jones potential. The positions of all layers close to the graphite surface or in a slit pore are considered to correspond to the minimum of the potential energy of the system. This model has led to accurate predictions of nitrogen and argon adsorption on carbon black at their normal boiling points. In the case of adsorption in slit pores, local isotherms are determined from the minimization of the grand potential. The model provides a reasonable description of the 0-1 monolayer transition, phase transition and packing effect. The adsorption of nitrogen at 77.35 K and argon at 87.29 K on activated carbons is analyzed to illustrate the potential of this theory, and the derived pore-size distribution is compared favourably with that obtained by the Density Functional Theory (DFT). The model is less time-consuming than methods such as the DFT and Monte-Carlo simulation, and most importantly it can be readily extended to the adsorption of mixtures and capillary condensation phenomena.