We probe and rationalize the complex surface chemistry of wurtzite ZnO by employing interatomic potential calculations coupled with a Monte Carlo procedure that sampled over 0.5 million local minima. We analyze the structure and stability of the ( 0001) and (0001̅ ) ZnO surfaces, rationalizing previous patterns found in STM images and explaining the (1 × 1) periodicity reported by LEED analysis. The full range of Zn/O surface occupancies was covered for a (5 × 5) supercell, keeping |m Zn − m O |/N ≈ 0.24 where m and N are the numbers of occupied surface sites and total surface sites, respectively. Our calculations explain why the (5 × 5) reconstructions seen in some experiments and highlight the importance of completely canceling the inherent dipole of the unreconstructed polar surfaces. The experimentally observed rich reconstruction patterns can be traced from the lowest occupancy, showing the thermodynamically most stable configurations of both polar surfaces. Triangular and striped reconstructions are seen, inter alia, on both polar surfaces, and hexagonal patterns also appear on the O terminated surface. Our results explain the main experimental structures observed on these complex surfaces. Moreover, grand canonical simulations of ZnO polar surfaces reveal that disorder is favored and, thus, configurational entropic factors is the the cause of their stability.

Neutron scattering methods observed complete room temperature conversion of methanol to framework methoxy in a commercial sample of methanol-to-hydrocarbons (MTH) catalyst H-ZSM-5, evidenced by methanol immobility and vibrational spectra matched by ab initio calculations. No methoxylation was observed in a commercial HY sample, attributed to the dealumination involved in high silica HY synthesis.

The source of n-type conductivity in undoped transparent conducting oxides has been a topic of debate for several decades. The point defect of most interest in this respect is the oxygen vacancy, but there are many conflicting reports on the shallow versus deep nature of its related electronic states. Here, using a hybrid quantum mechanical/molecular mechanical embedded cluster approach, we have computed formation and ionization energies of oxygen vacancies in three representative transparent conducting oxides: In 2 O 3 , SnO 2 , and ZnO. We find that, in all three systems, oxygen vacancies form well-localized, compact donors. We demonstrate, however, that such compactness does not preclude the possibility of these states being shallow in nature, by considering the energetic balance between the vacancy binding electrons that are in localized orbitals or in effective-mass-like diffuse orbitals. Our results show that, thermodynamically, oxygen vacancies in bulk In 2 O 3 introduce states above the conduction band minimum that contribute significantly to the observed conductivity properties of undoped samples. For ZnO and SnO 2 , the states are deep, and our calculated ionization energies agree well with thermochemical and optical experiments. Our computed equilibrium defect and carrier concentrations, however, demonstrate that these deep states may nevertheless lead to significant intrinsic n-type conductivity under reducing conditions at elevated temperatures. Our study indicates the importance of oxygen vacancies in relation to intrinsic carrier concentrations not only in In 2 O 3 , but also in SnO 2 and ZnO.

The atomic structure of inorganic nanoclusters obtained via a search for low lying minima on energy landscapes, or hypersurfaces, is reported for inorganic binary compounds: zinc oxide (ZnO)n, magnesium oxide (MgO)n, cadmium selenide (CdSe)n, and potassium fluoride (KF)n, where n = 1-12 formula units. The computational cost of each search is dominated by the effort to evaluate each sample point on the energy landscape and the number of required sample points. The effect of changing the balance between these two factors on the success of the search is investigated. The choice of sample points will also affect the number of required data points and therefore the efficiency of the search. Monte Carlo based global optimisation routines (evolutionary and stochastic quenching algorithms) within a new software package, viz. Knowledge Led Master Code (KLMC), are employed to search both directly and after pre-screening on the DFT energy landscape. Pre-screening includes structural relaxation to minimise a cheaper energy function - based on interatomic potentials - and is found to improve significantly the search efficiency, and typically reduces the number of DFT calculations required to locate the local minima by more than an order of magnitude. Although the choice of functional form is important, the approach is robust to small changes to the interatomic potential parameters. The computational cost of initial DFT calculations of each structure is reduced by employing Gaussian smearing to the electronic energy levels. Larger (KF)n nanoclusters are predicted to form cuboid cuts from the rock-salt phase, but also share many structural motifs with (MgO)n for smaller clusters. The transition from 2D rings to 3D (bubble, or fullerene-like) structures occur at a larger cluster size for (ZnO)n and (CdSe)n. Differences between the HOMO and LUMO energies, for all the compounds apart from KF, are in the visible region of the optical spectrum (2-3 eV); KF lies deep in the UV region at 5 eV and shows little variation. Extrapolating the electron affinities found for the clusters with respect to size results in the qualitatively correct work functions for the respective bulk materials.