Predicting Low-Pressure O2 Adsorption in Nanoporous Framework Materials for Sensing Applications

A set of 98 nanoporous framework material (NFM) structures was investigated by classical Grand canonical Monte Carlo simulations for low-pressure O₂ adsorption properties (Henry’s constant and isosteric heat of adsorption). The set of materials includes those that have shown high O₂ uptake experimentally as well as a subset of more than 2000 structures previously screened for noble-gas uptake. While use of the general force field UFF is fruitful for noble-gas adsorption studies, its use is shown to be limited for the case of O₂ adsorption—one distinct limitation is a lack of sufficient O₂–metal interactions to be able to describe O₂ interaction with open metal sites. Nonetheless, those structures without open metal sites that have very small pores (<2.5 Å) show increased O₂/N₂ selectivity. Additionally, O₂/N₂ mixture simulations show that in some cases, H₂O or N₂ can hinder O₂ uptake for NFMs with small pores due to competitive adsorption.     

DFT study of Pt adsorption on low index SrTiO3 surfaces: SrTiO3(1 0 0), SrTiO3(1 1 1) and SrTiO3(1 1 0)

Density Functional Theory has been used to determine the energetically preferred structures of submonolayer, monolayer, and multilayer Pt films on both ideal terminations of SrTiO₃(1 0 0), SrTiO₃(1 1 1), and SrTiO₃(1 1 0). The strength of the resulting metal/metal oxide interfaces was characterized by the adsorption energy of the film and the film’s work of separation. The two polar surfaces, SrTiO₃(1 1 1) and SrTiO₃(1 1 0), form significantly stronger interfaces than the non-polar SrTiO₃(1 0 0) surface. Approximate criteria were applied to predict the growth mode of Pt on each surface.     

Rapid diffusion of CH4/H2 mixtures in single-walled carbon nanotubes

Equilibrium molecular dynamics (EMD) are used to examine the self-diffusion and macroscopic diffusion of CH4/H2 mixtures adsorbed inside a (10,10) single-walled carbon nanotube. EMD can be used to determine the macroscopic diffusion coefficients of adsorbed mixtures by evaluating the matrix of Onsager transport coefficients. Earlier studies have indicated the diffusion of light gases adsorbed as single components in carbon nanotubes is extremely rapid compared to that in other known nanoporous materials. The results presented here indicate that extremely rapid diffusion can also occur for mixtures of adsorbed molecules. The rapid diffusion of adsorbed molecules and the strong coupling between the fluxes of the adsorbed species in a mixture have interesting implications for uses of carbon nanotubes in membrane-based applications.     

Enantioselective separation on a naturally chiral surface

Kinked-stepped, high Miller index surfaces of metal crystals are chiral and, therefore, exhibit enantiospecific properties. Previous temperature-programmed desorption (TPD) spectra have shown that the desorption energies of R-3-methylcyclohexanone (R-3-MCHO) on the chiral Cu(643)(R) and Cu(643)(S) surfaces are enantiospecific (J. Am. Chem. Soc. 2002, 124, 2384). Here, a comparison of the TPD spectra from Cu(111), Cu(221), Cu(533), Cu(653)(R&S), and Cu(643)(R&S) surfaces reveals that the enantiospecific desorption occurs from the chiral kink sites on the Cu(643) surfaces. Titration of the chiral kink sites with I atoms confirms this assignment of desorption features in the TPD spectra. Finally, the enantiospecific difference in the desorption energies of R- and S-3-MCHO has been used as the basis for demonstration of an enantioselective, kinetic separation of racemic 3-MCHO into its purified components during adsorption and desorption on the Cu(643)(R&S) surfaces.     

Testing predictions of macroscopic binary diffusion coefficients using lattice models with site heterogeneity

Quantitatively predicting mass transport rates for chemical mixtures in porous materials is important in applications of materials such as adsorbents, membranes, and catalysts. Because directly assessing mixture transport experimentally is challenging, theoretical models that can predict mixture diffusion coefficients using only single-component information would have many uses. One such model was proposed by Skoulidas, Sholl, and Krishna (Langmuir, 2003, 19, 7977), and applications of this model to a variety of chemical mixtures in nanoporous materials have yielded promising results. In this paper, the accuracy of this model for predicting mixture diffusion coefficients in materials that exhibit a heterogeneous distribution of local binding energies is examined. To examine this issue, single-component and binary mixture diffusion coefficients are computed using kinetic Monte Carlo for a two-dimensional lattice model over a wide range of lattice occupancies and compositions. The approach suggested by Skoulidas, Sholl, and Krishna is found to be accurate in situations where the spatial distribution of binding site energies is relatively homogeneous, but is considerably less accurate for strongly heterogeneous energy distributions.     

Examining the robustness of first-principles calculations for metal hydride reaction thermodynamics by detection of metastable reaction pathways

First principles calculations have played a useful role in screening mixtures of complex metal hydrides to find systems suitable for H(2) storage applications. Standard methods for this task efficiently identify the lowest energy reaction mechanisms among all possible reactions involving collections of materials for which DFT calculations have been performed. The resulting mechanism can potentially differ from physical reality due to inaccuracies in the DFT functionals used, or due to other approximations made in estimating reaction free energies. We introduce an efficient method to probe the robustness of DFT-based predictions that relies on identifying reactions that are metastable relative to the lowest energy reaction path predicted with DFT. An important conclusion of our calculations is that in many examples DFT cannot unambiguously predict a single reaction mechanism for a well defined metal hydride mixture because two or more mechanisms have reaction energies that differ by a small amount. Our approach is illustrated by analyzing a series of single step reactions identified in our recent work that examined reactions with a large database of solids [Kim et al., Phys. Chem. Chem. Phys. 2011, 13, 7218].     

Density functional theory study of beta-hydride elimination of ethyl on flat and stepped Cu surfaces

Plane wave density functional theory calculations have been used to characterize the transition states for beta-hydride elimination of ethyl on Cu(100), Cu(110), Cu(111), and Cu(221). The reaction rates predicted by these calculations have been compared to experiments by including tunneling corrections within harmonic transition state theory. Tunneling corrections are found to be important in describing the peak temperatures observed using temperature programmed desorption experiments on Cu(110), Cu(111), and Cu(221). Once these corrections are included, the effective activation energies obtained from our calculations are in good agreement with previous experimental studies of this reaction on these four Cu surfaces. The transition states determined in our calculations are used to examine two general hypotheses that have been suggested to describe structure sensitivity in metal-catalyzed surface reactions.     

Effect of Pr-Ca substitution on the transport and magnetic behavior of LaMnO3 perovskite

The effect of simultaneous substitution of a fluctuating cation and a divalent cation in LaMnO₃ perovskite modifies the properties of the material to exhibit large valence colossal magnetoresistance (CMR) effect. A good example of these properties is (La_1−2x Pr _x Ca _x )MnO₃ (LPCMO) type CMR material. In this communication it is reported that, with the increase in x (for x=0.1, 0.15, 0.2), the T _c varies between 100 and 120 K with improvisation in metal-insulator transition. Interestingly, resistance increases with x from few hundred ohms to few kilo ohms with corresponding decrease in the unit cell volume. The results of the studies using X-ray diffraction (XRD), electrical resistivity, magnetoresistance and ac susceptibility measurements on LPCMO samples for understanding the structural, transport and magnetic properties are discussed in detail.     

Characterization of enantiospecific chemisorption on chiral Cu surfaces vicinal to Cu(111) and Cu(100) using density functional theory

Surfaces of simple fcc metals such as Cu with nonzero and unequal Miller indices are intrinsically chiral. Density functional theory (DFT) calculations are a useful way to study the enantiospecific adsorption of small chiral molecules on these chiral metal surfaces. We report DFT calculations of seven chiral molecules on several structurally distinct chiral Cu surfaces. These surfaces include two surfaces with (111)-oriented terraces and one with (100)-oriented terraces. Calculations are also described on a surface that was modified to mimic the surface structures that typically appear on real metal surfaces following thermally driven fluctuations in step edges. Our results provide initial information on how variation in the surface structure of intrinsically chiral metal surfaces can affect the enantiospecific adsorption of small molecules on these surfaces.     

Efficient simulation of binary adsorption isotherms using transition matrix Monte Carlo

Molecular simulations of binary adsorption in porous materials are a useful complement to experimental studies of mixture adsorption. Most molecular simulations of binary adsorption are performed using grand canonical Monte Carlo (GCMC) to independently examine a range of state points of interest. A disadvantage of this approach is that it only yields information at a discrete set of state points; therefore, if a complete isotherm is required for arbitrary conditions, some type of data fitting or interpolation must be used in combination with the GCMC data. We show that the transition matrix Monte Carlo (TMMC) method of Shen and Errington (Shen, V. K.; Errington, J. R. J. Chem.Phys. 2005, 122, 064508) is well-suited to simulation of binary adsorption in porous materials. At the completion of a TMMC simulation, the adsorption isotherm for all possible bulk phase compositions and pressures is available without data fitting or interpolation. It is also straightforward to use results from TMMC to compute derivatives of the isotherm such as the mixture thermodynamic correction factors, partial differential ln f(i)/partial differential ln c(j), again without data fitting or interpolation. This approach should be useful in contexts where information on the full adsorption isotherm is needed, such as the design of adsorption- or membrane-based separations.