Zeolites, microporous aluminosilicates, are amongst the most widely used catalysts in the petrochemical industry. Zeolite catalytic functionality is influenced by the location of tetrahedral alumina and associated counter-cations in the aluminosilicate framework, yet little is definitively known about the factors that govern the framework aluminium arrangement. It is generally accepted that all zeolites obey Löwenstein's rule of “aluminium avoidance”, and that –Al–O–Al– bond formation is forbidden. Here, we describe an unprecedented screening of aluminium distribution in catalytically active zeolite SSZ-13 (CHA) in both its protonated and sodium-containing forms, H-SSZ-13 and Na-SSZ-13, using density functional theory (DFT). We predict violations of Löwenstein's rule in high and low silica H-SSZ-13 and other protonated frameworks considered in this investigation: H-LTA, H-RHO, H-ABW and H-MOR. The synthetic realisation of these zeolites could spur the development of new catalytic routes and materials, and the optimisation of existing zeolite catalysts.

The metal–organic framework (MOF) material UiO-66 has emerged as one of the most promising MOF materials due to its thermal and chemical stability and its potential for catalytic applications. Typically, as-synthesised UiO-66 has a relatively high concentration of missing linker defects. The presence of these defects has been correlated with catalytic activity but characterisation of defect structure has proved elusive. We refine a recent experimental determination of defect structure using static and dynamic first principles approaches, which reveals a dynamic and labile acid centre that could be tailored for functional applications in catalysis.

Many of the potential applications for metal–organic frameworks (MOFs) focus on exploiting their porosity for molecular storage, release, and separation, where the functional behavior is controlled by a subtle balance of host–guest interactions. Typically, the host structure is relatively unperturbed by the presence of guests; however, a subset of MOFs exhibit dramatic phase-change behavior triggered by the adsorption of guests or other stimuli, for which the MIL-53 material is an archetype. In this work, we use density functional approaches to examine the electronic structure changes associated with changes of phase and density and find the associated change in band gaps can be larger than 1 eV for known MIL-53-type materials and hypothecated structures. Moreover, we show that internal pressure (via guest molecules) and external pressure can exert a major influence on the band gap size and gap states. The large response in electronic properties to breathing transitions in MOFs could be exploitable in future applications in resistive switching, phase-change memory, piezoresistor, gas-sensing and thermochromic materials.

Ordinary water ice forms under ambient conditions and has two polytypes, hexagonal ice (Ih) and cubic ice (Ic). From a careful comparison of proton ordering arrangements in Ih and Ic using periodic density functional theory (DFT) and diffusion Monte Carlo (DMC) approaches, we find that the most stable arrangement of water molecules in cubic ice is isoenergetic with that of the proton ordered form of hexagonal ice (known as ice XI). We denote this potential new polytype of ice XI as XIc and discuss a possible route for preparing ice XIc.

Using density functional theory we investigate whether intrinsic defects in ice surface segregate. We predict that hydronium, hydroxide, and the Bjerrum L- and D-defects are all more stable at the surface. However, the energetic cost to create a D-defect at the surface and migrate it into the bulk crystal is smaller than its bulk formation energy. Absolute and relative segregation energies are sensitive to the surface structure of ice, especially the spatial distribution of protons associated with dangling hydrogen bonds. It is found that the basal plane surface of hexagonal ice increases the bulk concentration of Bjerrum defects, strongly favoring D-defects over L-defects. Dangling protons associated with undercoordinated water molecules are preferentially injected into the crystal bulk as Bjerrum D-defects, leading to a surface dipole that attracts hydronium ions. Aside from the disparity in segregation energies for the Bjerrum defects, we find the interactions between defect species to be very finely balanced; surface segregation energies for hydronium and hydroxide species and trapping energies of these ionic species with Bjerrum defects are equal within the accuracy of our calculations. The mobility of the ionic hydronium and hydroxide species is greatly reduced at the surface in comparison to the bulk due to surface sites with high trapping affinities. We suggest that, in pure ice samples, the surface of ice will have an acidic character due to the presence of hydronium ions. This may be important in understanding the reactivity of ice particulates in the upper atmosphere and at the boundary layer.

Ordinary water ice forms under ambient conditions and has two polytypes, hexagonal ice (Ih) and cubic ice (Ic). From a careful comparison of proton ordering arrangements in Ih and Ic using periodic density functional theory (DFT) and diffusion Monte Carlo (DMC) approaches, we find that the most stable arrangement of water molecules in cubic ice is isoenergetic with that of the proton ordered form of hexagonal ice (known as ice XI). We denote this potential new polytype of ice XI as XIc and discuss a possible route for preparing ice XIc.