Modeling micelle-templated mesoporous material SBA-15: atomistic model and gas adsorption studies

Unravelling the structure of CO2 in silica adsorbents: an NMR and computational perspective

This comprehensive review describes recent advancements in the use of solid-state NMR-assisted methods and computational modeling strategies to unravel gas adsorption mechanisms and CO2 speciation in porous CO2-adsorbent silica materials at the atomic scale. This work provides new perspectives for the innovative modifications of these materials rendering them more amenable to the use of advanced NMR methods.

Recent Advances in Computational Protocols Addressing Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDP) are abundant in the human genome and have recently emerged as major therapeutic targets for various diseases. Unlike traditional proteins that adopt a definitive structure, IDPs in free solution are disordered and exist as an ensemble of conformations. This enables the IDPs to signal through multiple signaling pathways and serve as scaffolds for multi-protein complexes. The challenge in studying IDPs experimentally stems from their disordered nature. Nuclear magnetic resonance (NMR), circular dichroism, small angle X-ray scattering, and single molecule Förster resonance energy transfer (FRET) can give the local structural information and overall dimension of IDPs, but seldom provide a unified picture of the whole protein. To understand the conformational dynamics of IDPs and how their structural ensembles recognize multiple binding partners and small molecule inhibitors, knowledge-based and physics-based sampling techniques are utilized in-silico, guided by experimental structural data. However, efficient sampling of the IDP conformational ensemble requires traversing the numerous degrees of freedom in the IDP energy landscape, as well as force-fields that accurately model the protein and solvent interactions. In this review, we have provided an overview of the current state of computational methods for studying IDP structure and dynamics and discussed the major challenges faced in this field.

Critical energy barrier for capillary condensation in mesopores: Hysteresis and reversibility

Capillary condensation in the regime of developing hysteresis occurs at a vapor pressure, Pcond, that is less than that of the vapor-like spinodal. This is because the energy barrier for the vapor-liquid transition from a metastable state at Pcond becomes equal to the energy fluctuation of the system; however, a detailed mechanism of the spontaneous transition has not been acquired even through extensive experimental and simulation studies. We therefore construct accurate atomistic silica mesopore models for MCM-41 and perform molecular simulations (gauge cell Monte Carlo and grand canonical Monte Carlo) for argon adsorption on the models at subcritical temperatures. A careful comparison between the simulation and experiment reveals that the energy barrier for the capillary condensation has a critical dimensionless value, Wc (*) = 0.175, which corresponds to the thermal fluctuation of the system and depends neither on the mesopore size nor on the temperature. We show that the critical energy barrier Wc (*) controls the capillary condensation pressure Pcond and also determines a boundary between the reversible condensation/evaporation regime and the developing hysteresis regime.

Effect of surface texture on freezing in nanopores: surface-induced versus homogeneous crystallization

Freezing of argon in ordered and disordered carbon pores of a similar diameter D ∼ 2.4 nm is investigated using extensive molecular simulations with large system sizes up to 10(4) atoms. While crystallization in the atomistically smooth pore consists in a surface-induced phase transition occurring at a temperature larger than the bulk, crystallization in the disordered pores, which is only partial as it is spatially restricted to the pore center, occurs through homogeneous crystallization. These results shed light on solidification in pores by showing that there is a crossover between surface-induced and homogeneous crystallization upon increasing the surface disorder of the host material. In the latter case, the Gibbs-Thomson equation, in which crystallization is assumed to occur when the crystal size equals the pore size corrected for the thickness of the unfreezable layer at the pore surface, is in reasonable agreement with the observed freezing temperature.

Ion-specific adsorption and electroosmosis in charged amorphous porous silica

Aqueous electrolyte solutions (NaCl, KCl, CsCl, and SrCl₂) confined in a negatively charged amorphous silica slit pore.

Surfactant behavior of ionic liquids involving a drug: from molecular interactions to self-assembly

Aggregates formed in an aqueous medium by three ionic liquids CnMImIbu made up of 1-alkyl-3-methyl-imidazolium cation (n = 4, 6, 8) and ibuprofenate anion are investigated. Dynamic light scattering (DLS), cryogenic transmission electron microscopy (cryo-TEM), (1)H nuclear magnetic resonance measurements, and atom-scale molecular dynamics simulations are used to shed light on the main interactions governing the formation of the aggregates and their composition. At high concentration, mixed micelles are formed with a composition that depends on the imidazolium alkyl chain length. For the shortest alkyl chain, micelles are mainly composed of ibuprofenate anions with some imidazolium cations intercalated between the anions. Upon increasing the alkyl chain length, the composition of the aggregates gets enriched in imidazolium cations and aggregates of stoichiometric composition are obtained. Attractive interactions between these aggregates led to the formation of larger aggregates. As suggested by molecular simulations, these larger aggregates might constitute the early stage of phase separation. Transitions from micelles to vesicles or ribbons are observed due to dilution effects and changes in the chemical composition of the aggregates. We also show that aggregation can be probed using simple microscopic quantities such as radial distribution functions and average solvation numbers.

Molecular simulation of adsorption and transport in hierarchical porous materials

Adsorption and transport in hierarchical porous solids with micro- (~1 nm) and mesoporosities (>2 nm) are investigated by molecular simulation. Two models of hierarchical solids are considered: microporous materials in which mesopores are carved out (model A) and mesoporous materials in which microporous nanoparticles are inserted (model B). Adsorption isotherms for model A can be described as a linear combination of the adsorption isotherms for pure mesoporous and microporous solids. In contrast, adsorption in model B departs from adsorption in pure microporous and mesoporous solids; the inserted microporous particles act as defects, which help nucleate the liquid phase within the mesopore and shift capillary condensation toward lower pressures. As far as transport under a pressure gradient is concerned, the flux in hierarchical materials consisting of microporous solids in which mesopores are carved out obeys the Navier-Stokes equation so that Darcy's law is verified within the mesopore. Moreover, the flow in such materials is larger than in a single mesopore, due to the transfer between micropores and mesopores. This nonzero velocity at the mesopore surface implies that transport in such hierarchical materials involves slippage at the mesopore surface, although the adsorbate has a strong affinity for the surface. In contrast to model A, flux in model B is smaller than in a single mesopore, as the nanoparticles act as constrictions that hinder transport. By a subtle effect arising from fast transport in the mesopores, the presence of mesopores increases the number of molecules in the microporosity in hierarchical materials and, hence, decreases the flow in the micropores (due to mass conservation). As a result, we do not observe faster diffusion in the micropores of hierarchical materials upon flow but slower diffusion, which increases the contact time between the adsorbate and the surface of the microporosity.

Atomistic model of micelle-templated mesoporous silicas: structural, morphological, and adsorption properties

The structural, morphological, and adsorption properties of MCM-41 porous silicas are investigated using a realistic numerical model obtained by means of ab initio calculations [Ugliengo, P.; et al. Adv. Mater.2008, 20, 1]. Simulated X-ray diffraction, small angle neutron scattering, and electronic microscopy for the atomistic model are in good agreement with experimental data. The morphological features are also assessed from chord length distributions and porous volume and specific geometrical surface calculations, etc. The N(2), CO(2), and H(2)O adsorption isotherms in the atomistic model of MCM-41 are also in reasonable agreement with their experimental counterpart. An important finding of the present work is that water forms a film adsorbed on specific hydrophilic regions of the surface while the rest of the surface is depleted in water molecules. This result suggests that the surface of MCM-41 materials is heterogeneous, as it is made up of both hydrophilic and hydrophobic patches. While adsorption and irreversible capillary condensation can be described using the thermodynamical approach by Derjaguin (also known as the Derjaguin-Broekhoff-De Boer model), the Freundlich equation fits nicely the data for reversible and continuous filling in small pores.

Capillary condensation of adsorbates in porous materials

Hysteresis in capillary condensation is important for the fundamental study and application of porous materials, and yet experiments on porous materials are sometimes difficult to interpret because of the many interactions and complex solid structures involved in the condensation and evaporation processes. Here we make an overview of the significant progress in understanding capillary condensation and hysteresis phenomena in mesopores that have followed from experiment and simulation applied to highly ordered mesoporous materials such as MCM-41 and SBA-15 over the last few decades.

Removal of patulin from aqueous solutions by propylthiol functionalized SBA-15

Propylthiol functionalized SBA-15 silica was investigated to detoxify aqueous solutions contaminated with the regulated mycotoxin patulin. Micelle templated silicas with a specific pore size were synthetically modified to possess propylthiol groups, a functional group known to form Michael reaction products with the conjugated double bond system of patulin. BET surface area analysis indicated the propylthiol functionalized SBA-15 possesses channels with the pore size of 5.4 nm and a surface area of 345 m(2)g(-1). Elemental analysis indicates the silicon/sulfur ratio to be 10:1, inferring one propylthiol substituent for every ten silica residues. The propylthiol modified SBA-15 was effective at significantly reducing high levels of patulin from aqueous solutions (pH 7.0) in batch sorption assays at room temperature. The material was less effective at lower pH; however heating low pH solutions and apple juice to 60 °C in the presence of propylthiol functionalized SBA-15 significantly reduced the levels of patulin in contaminated samples. Composite molecular models developed by semi-empirical PM3 and empirical force field methods support patulin permeation through the mesoporous channels of propylthiol functionalized SBA-15. Density functional study at the B3LYP/6-31G(d,p) level predicts the proposed patulin adducts formed by reaction with the thiol residues exhibit less electrophilic properties than patulin. It is demonstrated the use of propylthiol functionalized SBA-15 is a viable approach to reduce patulin levels in aqueous solutions, including contaminated apple juice.

Defocus image contrast in hexagonally-ordered mesoporous material

A transmission electron microscope was used to characterize a powder form of hexagonally-ordered mesoporous silica material. The structural symmetry built into this amorphous material allowed one to obtain three characteristic images, i.e. a hexagonal honey-comb structure and wide/narrow parallel lines. These images were found to originate primarily from phase contrast, which changed sensitively with defocusing. To further understand the contrast behaviour of these images, an analytical form of the defocus contrast theory was developed and applied to the simulation of the characteristic wide/narrow parallel line images. The simulation was found to be in good qualitative agreement with experiments, where changes in focus conditions and specimen thickness were predicted to alter the contrast in the resulting parallel-line type images.

Capillary condensation and evaporation in alumina nanopores with controlled modulations

Capillary condensation in nanoporous anodic aluminum oxide presenting not interconnected pores with controlled modulations is studied using adsorption experiments and molecular simulations. Both the experimental and simulation data show that capillary condensation and evaporation are driven by the smallest size of the nanopore (constriction). The adsorption isotherms for the open and closed pores are almost identical if constrictions are added to the system. The latter result implies that the type of pore ending does not matter in modulated pores. Thus, the presence of hysteresis loops observed in adsorption isotherms measured in straight nanopores with closed bottom ends can be explained in terms of geometrical inhomogeneities along the pore axis. More generally, these results provide a general picture of capillary condensation and evaporation in constricted or modulated pores that can be used for the interpretation of adsorption in disordered porous materials.

Molecular simulation of nitrogen adsorption in nanoporous silica

This article reports on a molecular simulation study of nitrogen adsorption and condensation at 77 K in atomistic silica cylindrical nanopores (MCM-41). Two models are considered for the nitrogen molecule and its interaction with the silica substrate. In the "pea" model, the nitrogen molecule is described as a single Lennard-Jones sphere and only Lennard-Jones interactions between the nitrogen molecule and the oxygens atoms of the silica substrate are taken into account. In the "bean" model (TraPPE force field), the nitrogen molecule is composed of two Lennard-Jones sites and a linear array of three charges on the atomic positions and at the center of the nitrogen-nitrogen bond. In the bean model, the interactions between the sites on the nitrogen molecule and the Si, O, and H atoms of the substrate are the sum of the Coulombic and dispersion interactions with a repulsive short-range contribution. The data obtained with the pea and bean models in silica nanopores conform to the typical behavior observed in the experiments for adsorption/condensation in cylindrical MCM-41 nanopores; the adsorbed amount increases continuously in the multilayer adsorption regime until an irreversible jump occurs because of capillary condensation and evaporation of the fluid within the pore. Our results suggest that the pea model can be used for characterization purposes where one is interested in capturing the global experimental behavior upon adsorption and desorption in silica nanopores. However, the bean model is more suitable to investigating the details of the interaction with the surface because this model, which accounts for the partial charges located on the nitrogen atoms of the molecule (quadrupole), allows a description of the specific interactions between this adsorbate and silica surfaces (silanol groups and siloxane bridges) or grafted silica surfaces. In particular, the bean model provides a more realistic picture of nitrogen adsorption in the vicinity of silica surfaces or confined in silica nanopores, where the isosteric heat of adsorption curves show that the nitrogen molecule in this model is sensitive to the surface heterogeneity.

Adsorption and structure of benzene on silica surfaces and in nanopores

Grand canonical Monte Carlo simulations are used to study the adsorption of benzene on atomistic silica surfaces and in cylindrical nanopores. The effect of temperature and surface chemistry is addressed by studying the adsorption of benzene at 293 and 323 K on both fully and partially hydroxylated silica surfaces or nanopores. We also consider the adsorption of benzene in a cylindrical nanopore of diameter D=3.6 nm and compare our results with those obtained for planar surfaces. The structure of benzene in the vicinity of the planar surface and confined in the cylindrical nanopore is determined by calculating orientational order parameters and examining positional pair correlation functions. The density profiles of adsorbed benzene reveal the strong layering of the adsorbate, which decays with the distance from the silica surface. At a given temperature and at low pressures, the film adsorbed at the fully hydroxylated silica surface is larger than that for the partially hydroxylated silica surface. This result is due to an increase in the density of silanol groups that induces an increase in the polarity of the silica surface, which becomes more attractive for the adsorbate. Our results also suggest that the benzene molecules prefer an orientation in which their ring is nearly perpendicular to the surface when fully hydroxylated surfaces are considered. When partially hydroxylated surfaces are considered, a second preferential orientation is observed where the benzene ring forms an angle of approximately 50 degrees with the silica surface. In this case, the average orientation of the benzene molecules appears disordered as in the bulk phase. These results suggest that determining the experimental orientation of benzene in the vicinity of a silica surface is a difficult task even when the surface chemistry is known. Capillary condensation in the nanopores involves a transition from a partially filled pore (a thin film adsorbed at the pore surface) to a completely filled pore configuration where the confined liquid coexists at equilibrium with the external gas phase. The disordered orientation of the adsorbed benzene molecules in the case of the partially hydroxylated surface favors the condensation of benzene molecules (the condensation pressure for this substrate is lower than that for the fully hydroxylated surface). Finally, these results are consistent with the structural analysis, showing that (1) benzene tends to relax its liquid structure a little in order to optimize its molecular arrangement near the pore wall and (2) the disordering of the liquid structure induced by the surface becomes stronger as the interaction with the pore wall increases.