Open-Source, Python-Based Redevelopment of the ChemShell Multiscale QM/MM Environment

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 160000 atoms.

Thermodynamically accessible titanium clusters TiN, N = 2–32

We have performed a genetic algorithm search on the tight-binding interatomic potential energy surface (PES) for small Ti_N (N = 2–32) clusters.

The adsorption of Cu on the CeO2(110) surface

A detailed density functional theory (DFT) study coupled with extended X-ray absorption fine structure (EXAFS) experiments on the geometrical and electronic properties of copper species on CeO₂ surface demonstrating the effects of oxidation state and solvent environment.

Room temperature methoxylation in zeolites: insight into a key step of the methanol-to-hydrocarbons process

Neutron scattering methods studying mobility and vibrational spectra observed complete room temperature methoxylation in a commercial sample of methanol-to-hydrocarbons (MTH) catalyst H-ZSM-5.

Embedded-cluster calculations in a numeric atomic orbital density-functional theory framework

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 TiO2(110).

Double bubbles: a new structural motif for enhanced electron-hole separation in solids

Electron-hole separation for novel composite systems comprised of secondary building units formed from different compounds is investigated with the aim of finding suitable materials for photocatalysis. Pure and mixed SOD and LTA superlattices of (ZnO)12 and (GaN)12, single-shell bubbles are constructed as well as core@shell single component frameworks composed of larger (ZnO)48 and (GaN)48 bubbles with each containing one smaller bubble. Enthalpies of formation for all systems are comparable with fullerenes. Hole and electron separation is achieved most efficiently by the edge sharing framework composed of (GaN)12@(ZnO)48 double bubbles, with the hole localised on the nitrogen within the smaller bubbles and the excited electron on zinc within the larger cages.