Hydrogen-Bonded Water-Aminium Assemblies for Synthesis of Zeotypes with Ordered Heteroatoms.

Water plays a central role in the crystallization of a variety of organic, inorganic, biological, and hybrid materials. This is also true for zeolites and zeolite-like materials, an important class of industrial catalysts and adsorbents. Water is always present during their hydrothermal synthesis, either with or without organic species as structure-directing agents. Apart from its role as a solvent or a catalyst, structure direction by water in zeolite synthesis has never been clearly elucidated. Here, we report the crystallization of phosphate-based molecular sieves using rationally designed, hydrogen-bonded water-aminium assemblies, resulting in molecular sieves exhibiting the crystallographic ordering of heteroatoms. We demonstrate that a 1:1 assembly of water and diprotonated N,N-dimethyl-1,2-ethanediamine acts as a structure-directing agent in the synthesis of a silicoaluminophosphate material with phillipsite (PHI) topology, using SMARTER crystallography, which combines single-crystal X-ray diffraction and nuclear magnetic resonance spectroscopy, as well as ab initio molecular dynamics simulations. The molecular arrangement of the hydrogen-bonded assembly matches well with the shape and size of subunits in the PHI structure, and their charge distributions result in the strict ordering of framework tetrahedral atoms. This concept of structure direction by water-containing supramolecular assemblies should be applicable to the synthesis of many classes of porous materials.

3D-3D topotactic transformation in aluminophosphate molecular sieves and its implication in new zeolite structure generation.

Zeolites have unique pore structures of molecular dimensions and tunable compositions, making them ideal for shape selective catalysis and separation. However, targeted synthesis of zeolites with new pore structures and compositions remains a key challenge. Here, we propose an approach based on a unique 3D-3D topotactic transformation, which takes advantage of weak bonding in zeolites. This is inspired by the structure transformation of PST-5, a new aluminophosphate molecular sieve, to PST-6 by calcination. The structure of nano-sized PST-5 crystals is determined by 3D electron diffraction. We find that the 3D-3D topotactic transformation involves two types of building units where penta- or hexa-coordinated Al is present. We apply this approach to several other zeolite systems and predict a series of new zeolite structures that would be synthetically feasible. This method provides a concept for the synthesis of targeted zeolites, especially those which may not be feasible by conventional methods.

Molecular dynamics study of tridymite.

Structural changes in tridymite have been investigated by molecular dynamics simulation. Two thermal processes were carried out, one cooling from the high-temperature hexagonal structure of tridymite (HP-tridymite) and the other heating from the low-temperature monoclinic structure of tridymite (MX1-tridymite). The former process showed that HP, LHP (low-temperature hexagonal structure), OC (orthorhombic structure with C2221 symmetry) and OP (orthorhombic structure with P212121 symmetry)-like structures appeared in sequence. In contrast, the latter process showed that MX1, OP, OC, LHP and HP-like structures appeared in sequence. Detailed analysis of the calculated structures showed that the configuration underwent stepwise changes associated with several characteristic modes. First, the structure of HP-tridymite determined from diffraction experiments was identified as a time-averaged structure in a similar manner to ß-cristobalite, thus indicating the important role of floppy modes of oxygen atoms at high temperature - one of the common features observed in silica crystals and glass. Secondly, the main structural changes were ascribed to a combination of distortion of the six-membered rings in the layers and misalignment between layers. We suggest that the slowing down of floppy oxygen movement invokes the multistage emergence of structures with lower symmetry on cooling. This study therefore not only reproduces the sequence of the main polymorphic transitions in tridymite, except for the appearance of the monoclinic phase, but also explains the microscopic dynamic structural changes in detail.

Mechanisms of CO2 capture in ionic liquids: a computational perspective.

We present computational studies of CO2 sorption in two different classes of ionic liquid. The addition of carbon dioxide to four superbase ionic liquids, [P3333][Benzim], [P3333][124Triz], [P3333][123Triz] and [P3333][Bentriz], was studied using the DFT approach and considering anions alone and individual ion pairs. The addition of CO2 to the anion alone clearly resulted in the formation of a covalently-bound carbamate function with the strength of binding correlated to experimental capacity. In the ion pair however the cation significantly alters the nature of the bonding such that the overall cohesive energy is reduced. Formation of a strong carbamate function occurs at the expense of weakening the interaction between anion and cation. In [N1111][l-ALA], a representative amino acid ionic liquid, evidence was found for a low-energy monomolecular mechanism for carbamate formation, explaining the 1 : 1 molar uptake ratio observed in some amino acid ionic liquids. The mechanism involves proton transfer to the carboxylate group of the aminate anion.

The addition of CO2 to four superbase ionic liquids: a DFT study.

The addition of carbon dioxide to four superbase ionic liquids, [P3333][Benzim], [P3333][124Triz], [P3333][123Triz] and [P3333][Bentriz] was studied using a molecular DFT approach involving anions alone and individual ion pairs. Intermolecular bonding within the individual ion pairs is characterised by a number of weak hydrogen bonds, with the superbase anion geometrically arranged so as to maximize interactions between the heterocyclic N atoms and the cation. The pairing energies show no correlation to the observed CO2 adsorption capacity. Addition of CO2 to the anion alone clearly resulted in the formation of a covalently-bound carbamate function with the strength of binding correlated to experimental capacity. In the ion pair however the cation significantly alters the nature of the bonding such that the overall cohesive energy is reduced. Formation of a strong carbamate function occurs at the expense of weakening the interaction between anion and cation. In the more weakly absorbing ion pairs which contain [123Triz](-) and [Bentriz](-), the carbamate-functionalised systems are very close in energy to adducts in which CO2 is more weakly bound, suggesting an equilibrium between the chemi- and physisorbed CO2.

Advances in theory and their application within the field of zeolite chemistry.

Zeolites are versatile and fascinating materials which are vital for a wide range of industries, due to their unique structural and chemical properties, which are the basis of applications in gas separation, ion exchange and catalysis. Given their economic impact, there is a powerful incentive for smart design of new materials with enhanced functionalities to obtain the best material for a given application. Over the last decades, theoretical modeling has matured to a level that model guided design has become within reach. Major hurdles have been overcome to reach this point and almost all contemporary methods in computational materials chemistry are actively used in the field of modeling zeolite chemistry and applications. Integration of complementary modeling approaches is necessary to obtain reliable predictions and rationalizations from theory. A close synergy between experimentalists and theoreticians has led to a deep understanding of the complexity of the system at hand, but also allowed the identification of shortcomings in current theoretical approaches. Inspired by the importance of zeolite characterization which can now be performed at the single atom and single molecule level from experiment, computational spectroscopy has grown in importance in the last decade. In this review most of the currently available modeling tools are introduced and illustrated on the most challenging problems in zeolite science. Directions for future model developments will be given.

Cation-exchanged SAPO-34 for adsorption-based hydrocarbon separations: predictions from dispersion-corrected DFT calculations.

The influence of the nature of the cation on the interaction of the silicoaluminophosphate SAPO-34 with small hydrocarbons (ethane, ethylene, acetylene, propane, propylene) is investigated using periodic density-functional theory calculations including a semi-empirical dispersion correction (DFT-D). Initial calculations are used to evaluate which of the guest-accessible cation sites in the chabazite-type structure is energetically preferred for a set of ten cations, which comprises four alkali metals (Li(+), Na(+), K(+), Rb(+)), three alkaline earth metals (Mg(2+), Ca(2+), Sr(2+)), and three transition metals (Cu(+), Ag(+), Fe(2+)). All eight cations that are likely to be found at the SII site (centre of a six-ring) are then included in the following investigation, which studies the interaction with the hydrocarbon guest molecules. In addition to the interaction energies, some trends and peculiarities regarding the adsorption geometries are analysed, and electron density difference plots obtained from the calculations are used to gain insights into the dominant interaction types. In addition to dispersion interactions, electrostatic and polarisation effects dominate for the main group cations, whereas significant orbital interactions are observed for unsaturated hydrocarbons interacting with transition metal (TM) cations. The differences between the interaction energies obtained for pairs of hydrocarbons of interest (such as ethylene-ethane and propylene-propane) deliver some qualitative insights: if this energy difference is large, it can be expected that the material will exhibit a high selectivity in the adsorption-based separation of alkene-alkane mixtures, which constitutes a problem of considerable industrial relevance. While the calculations show that TM-exchanged SAPO-34 materials are likely to exhibit a very high preference for alkenes over alkanes, the strong interaction may render an application in industrial processes impractical due to the large amount of energy required for regeneration. In this respect, SAPOs exchanged with alkaline earth cations could provide a better balance between selectivity and energy cost of regeneration.

Thermal rearrangement mechanisms in icosahedral carboranes and metallocarboranes.

Ab initio MD and potential energy surface sampling has been used to study the rearrangement processes in carboranes and their derivatives. A new mechanism is found, in addition to those previously proposed. The fact that theoretical activation energies are lower than those observed experimentally, and the differing activity of technetium and rhenium complexes, are rationalised by orbital symmetry constraints.

Isomorphism between ice and silica.

Both ice and silica crystallize into solid-state structures composed of tetrahedral building units that are joined together to form an infinite four-connected net. Mathematical considerations suggest that there is a vast number of such nets and thus potential crystal structures. It is therefore perhaps surprising to discover that, despite the differences in the nature of interatomic interactions in these materials, a fair number of commonly observed ice and silica phases are based on common nets. Here we use computer simulation to investigate the origin of this symmetry between the structures formed for ice and silica and to attempt to understand why it is not complete. We start from a comparison of the dense phases and then move to the relationship between the different open (zeolitic and clathratic) structures formed for both materials. We show that there is a remarkably strong correlation between the energetics of isomorphic silica and water ice structures and that this correlation arises because of the strong link between the total energy of a material and its local geometric features. Finally, we discuss a number of as yet unsynthesized low-energy structures which include a phase of ice based on quartz, a silica based on the structure of ice VI, and an ice clathrate that is isomorphic to the silicate structure nonasil.

Impact of ligands on CO(2) adsorption in metal-organic frameworks: First principles study of the interaction of CO(2) with functionalized benzenes. II. Effect of polar and acidic substituents.

Intermolecular interactions between the CO(2) molecule and a range of functionalized aromatic molecules have been investigated using density functional theory. The work is directed toward the design of linker molecules which could form part of new metal-organic framework materials with enhanced affinity for CO(2) adsorption at low pressure. Here, the focus was on the effect of introducing polar side groups, and therefore functionalized benzenes containing -NO(2), -NH(2), -OH, -SO(3)H, and -COOH substituents were considered. The strongest types of intermolecular interactions were found to be: (i) between lone pair donating atoms (N,O) of the side groups and the C of CO(2) (enhancement in binding energy of up to 8 kJ mol(-1) compared to benzene); and (ii) hydrogen bond interactions between acidic protons (of COOH and SO(3)H groups) and CO(2) oxygen (enhancement of 3-4 kJ mol(-1)). Both of these types of interaction have the effect of polarizing the CO(2) molecule. Weaker types of binding include hydrogen-bond-like interactions with aromatic H and pi-quadrupole interactions. The strongest binding is found when more than one interaction occurs simultaneously, as in C(6)H(5)SO(3)H and C(6)H(5)COOH, where simultaneous lone pair donation and H-bonding result in binding energy enhancements of 10 and 11 kJ mol(-1), respectively.

Tracking the formation of cobalt substituted ALPO-5 using simultaneous in situ X-ray diffraction and X-ray absorption spectroscopy techniques.

We followed the formation of metal ion substituted aluminophosphate, a redox molecular sieve, using a newly developed in situ simultaneous X-ray diffraction and X-ray absorption spectroscopic technique. The study showed that a cobalt-phosphorous network forms prior to the crystallisation.

Nanoporous oxidic solids: the confluence of heterogeneous and homogeneous catalysis.

The several factors that render certain kinds of nanoporous oxidic solids valuable for the design of a wide range of new heterogeneous catalysts are outlined and exemplified. These factors include: (i), their relative ease of preparation, when both mesoporous siliceous frameworks (ca. 20 to 250 A diameter pores) and microporous framework-substituted aluminophosphates (ca. 4 to 14 A diameter pores) can be tailored to suit particular catalytic needs according to whether regiospecific or enantio- or shape-selective conversions are the goal; (ii), the enormous internal (three-dimensional) areas that these nanoporous solids possess (typically 10(3) m(2) g(-1)) and the consequential ease of access of reactants through the internal pores of the solids; (iii), the ability, by judicious solid-state preparative methods to assemble spatially isolated, single-site active centres at the internal surfaces of these open-structure solids, thereby making the heterogeneous catalyst simulate the characteristic features of homogenous and enzymatic catalysts; (iv), the wide variety of in situ, time-resolved and ex situ experimental techniques, coupled with computational methods, that can pin-point the precise structure of the active site under operating conditions and facilitate the formulation of reaction intermediates and mechanisms. Varieties of catalysts are described for the synthesis (often under environmentally benign and solvent-free conditions) of a wide range of organic materials including commodity chemicals (such as adipic and terephthalic acid), fine and pharmaceutical chemicals (e.g. vitamin B(3)), alkenes, epoxides, and for the photocatalytic preferential destruction of carbon monoxide in the presence of hydrogen. Nanoporous oxidic solids are ideal materials to investigate the phenomenology of catalysis because, in many of them, little distinction exists between a model and a real catalyst.

Impact of ligands on CO2 adsorption in metal-organic frameworks: first principles study of the interaction of CO2 with functionalized benzenes. I. Inductive effects on the aromatic ring.

Intermolecular interactions between the CO(2) molecule and a range of functionalized aromatic molecules have been investigated using density functional theory. The work is directed toward the design of linker molecules which could form part of new metal-organic framework materials with enhanced affinity for CO(2) adsorption at low pressure. Two classes of substituted benzene molecules were considered: (i) with halogen substituents (tetrafluoro-, chloro-, bromo-, and dibromobenzene) and (ii) with methyl substituents (mono-, di-, and tetramethylbenzene). In the benzene-CO(2) complex, the main interaction is between the delocalized pi aromatic system and the molecular quadrupole of CO(2). Halogen substituents have an electron-withdrawing effect on the ring which destabilizes the pi-quadrupole interaction. Weak "halogen-bond" and hydrogen bondlike interactions partially compensate for this, but not to the extent that any significant enhancement of the intermolecular binding energy is observed. Methyl groups, on the other hand, have a positive inductive effect which strengthens the CO(2)-aromatic interaction by up to 3 kJ mol(-1) in the case of tetramethylbenzene. Weak hydrogen bondlike interactions with methyl H also contribute to the stability of the complexes.

Isomorphism of anhydrous tetrahedral halides and silicon chalcogenides: energy landscape of crystalline BeF2, BeCl2, SiO2, and SiS2.

We employ periodic density functional theory calculations to compare the structural chemistry of silicon chalcogenides (silica, silicon sulfide) and anhydrous tetrahedral halides (beryllium fluoride, beryllium chloride). Despite the different formal oxidation states of the elements involved, the divalent halides are known experimentally to form crystal structures similar to known SiX2 frameworks; the rich polymorphic chemistry of SiO2 is however not matched by divalent halides, for which a very limited number of polymorphs are currently known. The calculated energy landscapes yield a quantitative match between the relative polymorphic stability in the SiO2/BeF2 pair, and a semiquantitative match for the SiS2/BeCl2 pair. The experimentally observed polymorphs are found to lie lowest in energy for each composition studied. For the two BeX2 compounds studied, polymorphs not yet synthesized are predicted to lie very low in energy, either slightly above or even in between the energy of the experimentally observed polymorphs. The experimental lack of polymorphism for tetrahedral halide materials thus does not appear to stem from a lack of low-energy polymorphs but more likely is the result of a lack of experimental exploration. Our calculations further indicate that the rich polymorphic chemistry of SiO2 can be potentially matched, if not extended, by BeF2, provided that milder synthetic conditions similar to those employed in zeolite synthesis are developed for BeF2. Finally, our work demonstrates that both classes of materials show the same behavior upon replacement of the 2p anion with the heavier 3p anion from the same group; the thermodynamic preference shifts from structures with large rings to structures with larger fractions of small two and three membered rings.

A zeolite family with chiral and achiral structures built from the same building layer.

Porosity and chirality are two of the most important properties for materials in the chemical and pharmaceutical industry. Inorganic microporous materials such as zeolites have been widely used in ion-exchange, selective sorption/separation and catalytic processes. The pore size and shape in zeolites play important roles for specific applications. Chiral inorganic microporous materials are particularly desirable with respect to their possible use in enantioselective sorption, separation and catalysis. At present, among the 179 zeolite framework types reported, only three exhibit chiral frameworks. Synthesizing enantiopure, porous tetrahedral framework structures represents a great challenge for chemists. Here, we report the silicogermanates SU-32 (polymorph A), SU-15 (polymorph B) (SU, Stockholm University) and a hypothetical polymorph C, all built by different stacking of a novel building layer. Whereas polymorphs B and C are achiral, each crystal of polymorph A exhibits only one hand and has an intrinsically chiral zeolite structure. SU-15 and SU-32 are thermally stable on calcination.

On the performance of DFT and interatomic potentials in predicting the energetics of (three-membered ring-containing) siliceous materials.

We compare the enthalpies of transition for a range of SiO2 phases, including siliceous zeolites and dense polymorphs, calculated using three different interatomic potentials (Sanders-Leslie-Catlow (SLC), Sastre-Gale (SG), van Beest-Kramers-van Santen (BKS)), and from B3LYP periodic DFT calculations, with the experimentally measured values. It is found that the calculated results show a linear correlation with the measured values but that they often either considerably underestimate or overestimate enthalpy differences compared to experiment. Care should thus be taken when comparing experimental and calculated results. A linear rescaling of the calculated enthalpies to put the data on the same energy scale is proposed. Furthermore, it is found that when comparing enthalpies of transitions for materials containing three-membered rings, for which there is no experimental data available, the values, rescaled to the experimental energy scale, are very similar for both DFT and interatomic potentials (except for the BKS potential). The latter result suggests that the energetics of three-membered ring containing materials is well described using both approaches. Finally, we discuss the transition enthalpies of four three-membered ring containing siliceous materials and demonstrate that three-membered ring containing materials are not necessarily energetically disadvantageous but do become so progressively with increasing number of three-membered rings.

Diffusion of methanol in zeolite NaY: a molecular dynamics study.

Molecular dynamics simulations were performed in order to obtain a detailed understanding of the self-diffusion mechanisms of methanol in the zeolite NaY system. We derived a new force-field term to describe the interactions between the methanol molecules and the extraframework cations. From the simulations, we show that diffusive behavior in the high-temperature range consists of a combination of both short- and long-range motions at low and intermediate loadings. This type of motion is characterized by an activation energy that decreases as the loading increases. At low loadings, we also observe short-range diffusive behavior based on a surface-mediated mechanism. The short-range behavior corresponds to motion only on the length scale of an FAU supercage, whereas the long-range behavior involves intercage diffusion. For the saturation loading corresponding to 96 methanol molecules per unit cell, only short-range motions within the same supercage predominate. Finally, the preferential arrangement of the adsorbate molecules around the extraframework cations are examined and compared with those previously deduced from experimental data.

Cation migration upon adsorption of methanol in NaY and NaX faujasite systems: a molecular dynamics approach.

Molecular dynamics simulations have been carried out to address the question of cation migration upon adsorption of methanol in NaY and NaX faujasite systems as a function of the loading. For NaY, it has been shown that, at low and intermediate loadings, SII cations can migrate toward the center of the supercage due to strong interactions with the adsorbates, followed by a hopping of SI' from the sodalite cage into the supercage to fill the vacant SII site. A SI' cation can also migrate across the double six ring and takes a SI' vacant position. SI cations mainly remain trapped in their initial sites whatever the loading. At high loading, only limited motions are observed for SII cations due to steric effects induced by the presence of adsorbates within the supercage. For NaX, the SIII' cations which occupy the most accessible adsorption sites are significantly moving upon coordination to the methanol molecules; the extent of this mobility exhibits a maximum for 48 methanol molecules per unit cell before decreasing at higher loadings due to steric hindrance. In addition, the SI' and SII cations remain almost trapped in their initial sites whatever the loading. Indeed, the most probable migration mechanism involves SIII' cation displacements into nearby SIII' sites.