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.

ChemShell is a scriptable computational chemistry environment with an emphasis on multiscale simulation of complex systems using combined quantum mechanical and molecular mechanical (QM/MM) methods. Motivated by a scientific need to efficiently and accurately model chemical reactions on surfaces and within microporous solids on massively parallel computing systems, we present a major redevelopment of the ChemShell code, which provides a modern platform for advanced QM/MM embedding models. The new version of ChemShell has been re-engineered from the ground up with a new QM/MM driver module, an improved parallelization framework, new interfaces to high performance QM and MM programs, and a user interface written in the Python programming language. The redeveloped package is capable of performing QM/MM calculations on systems of significantly increased size, which we illustrate with benchmarks on zirconium dioxide nanoparticles of over 160,000 atoms.

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.

Atomic structure prediction, using KLMC (Lamarckian evolutionary algorithm search), and properties comparison of (KF) _n , (MgO) _n , (ZnO) _n and (CdSe) _n nanoclusters.

We integrate the all-electron electronic structure code FHI-aims into the general ChemShell package for solid-state embedding quantum and molecular mechanical (QM/MM) calculations. A major undertaking in this integration is the implementation of pseudopotential functionality into FHI-aims to describe cations at the QM/MM boundary through effective core potentials and therewith prevent spurious overpolarization of the electronic density. Based on numeric atomic orbital basis sets, FHI-aims offers particularly efficient access to exact exchange and second order perturbation theory, rendering the established QM/MM setup an ideal tool for hybrid and double-hybrid level density functional theory calculations of solid systems. We illustrate this capability by calculating the reduction potential of Fe in the Fe-substituted ZSM-5 zeolitic framework and the reaction energy profile for (photo-)catalytic water oxidation at TiO₂(110).