Hybrid quantum mechanical/molecular mechanical (QM/MM) methods are a powerful computational tool for the investigation of all forms of catalysis, as they allow for an accurate description of reactions occurring at catalytic sites in the context of a complicated electrostatic environment. The scriptable computational chemistry environment ChemShell is a leading software package for QM/MM calculations, providing a flexible, high performance framework for modelling both biomolecular and materials catalysis. We present an overview of recent applications of ChemShell to problems in catalysis and review new functionality introduced into the redeveloped Python-based version of ChemShell to support catalytic modelling. These include a fully guided workflow for biomolecular QM/MM modelling, starting from an experimental structure, a periodic QM/MM embedding scheme to support modelling of metallic materials, and a comprehensive set of tutorials for biomolecular and materials modelling.
We have studied the clusters involved in the initial stages of nucleation of Zeolitic Imidazolate Frameworks, employing a wide range of computational techniques. In the pre-nucleating solution, the prevalent cluster is the ZnIm4 cluster (formed by a zinc cation, Zn2+, and four imidazolate anions, Im-), although clusters such as ZnIm3, Zn2Im7, Zn2Im7, Zn3Im9, Zn3Im10, or Zn4Im12 have energies that are not much higher, so they would also be present in solution at appreciable quantities. All these species, except ZnIm3, have a tetrahedrally coordinated Zn2+ cation. Small ZnxImy clusters are less stable than the ZnIm4 cluster. The first cluster that is found to be more stable than ZnIm4 is the Zn41Im88 cluster, which is a disordered cluster with glassy structure. Bulk-like clusters do not begin to be more stable than glassy clusters until much larger sizes, since the larger cluster we have studied (Zn144Im288) is still less stable than the glassy Zn41Im88 cluster, suggesting that Ostwald's rule (the less stable polymorph crystallizes first) could be fulfilled, not for kinetic, but for thermodynamic reasons. Our results suggest that the first clusters formed in the nucleation process would be glassy clusters, which then undergo transformation to any of the various crystal structures possible, depending on the kinetic routes provided by the synthesis conditions. Our study helps elucidate the way in which the various species present in solution interact, leading to nucleation and crystal growth.
Cu impurities are reported to have significant effects on the electrical and optical properties of bulk ZnO. In this work, we study the defect properties of Cu in ZnO using hybrid quantum mechanical/molecular mechanical (QM/MM)-embedded cluster calculations based on a multi-region approach that allows us to model defects at the true dilute limit, with polarization effects described in an accurate and consistent manner. We compute the electronic structure, energetics, and geometries of Cu impurities, including substitutional and interstitial configurations, and analyze their effects on the electronic structure. Under ambient conditions, CuZn is the dominant defect in the d9 state and remains electronically passive. We find that, however, as we approach typical vacuum conditions, the interstitial Cu defect becomes significant and can act as an electron trap.
We review the current techniques used in the prediction of crystal structures and their surfaces and of the structures of nanoparticles. The main classes of search algorithm and energy function are summarized, and we discuss the growing role of methods based on machine learning. We illustrate the current status of the field with examples taken from metallic, inorganic and organic systems. This article is part of a discussion meeting issue 'Dynamic in situ microscopy relating structure and function'.
We report a new class of carbon nanostructures at a lower sub-nano end of the size scale with a surprising stability, as compared to the well-known carbon fullerenes. The octahedral carbon clusters contain tetragonal rings, which, in spite of a common belief, prove to be an energy efficient means of plying graphene sheets to make three-dimensional spheroid shapes, similar to fullerenes. The two families of structures are shown to be competitive at small sizes (â¼20 atoms) at room temperature, and for higher temperatures, at both small and large sizes (>200 atoms). Our calculations demonstrate that both vibrational and electronic spectra of these cluster families are similar, which thus might cloud their experimental identification. However, there is a sufficiently strong shift in vibrational frequencies below 160 and in the range of 600-800 cm-1, which should help to identify different types of carbon clusters experimentally. We propose octahedral clusters and other structures containing tetragonal rings as viable structural elements and building units in inorganic chemistry and materials science of carbon along with fullerenes.
Magnesium-rich silicates are ubiquitous both terrestrially and astronomically, where they are often present as small particles. Nanosized Mg-rich silicate particles are likely to be particularly important for understanding the formation, processing, and properties of cosmic dust grains. Although astronomical observations and laboratory studies have revealed much about such silicate dust, our knowledge of this hugely important class of nanosolids largely rests on top-down comparisons with the properties of bulk silicates. Herein, we provide a foundational bottom-up study of the structure and properties of Mg-rich nanosilicates based on carefully procured atomistic models. Specifically, we employ state-of-the-art global optimization methods to search for the most stable structures of silicate nanoclusters with olivine (Mg2SiO4) N and pyroxene (MgSiO3) N compositions with N = 1-10. To ensure the reliability of our searches, we develop a new interatomic potential that has been especially tuned for nanosilicates. Subsequently, we refine these searches and calculate a range of physicochemical properties of the most stable nanoclusters using accurate density functional theory based electronic structure calculations. We report a detailed analysis of structural and energy properties, charge distributions, and infrared vibrational spectra, where in all cases we compare our finding for nanosilicates with those of the corresponding bulk silicate crystals. For most properties considered, we find large differences with respect to the bulk limit, underlining the limitations of a top-down approach for describing these species. Overall, our work provides a new platform for an accurate and detailed understanding of nanoscale silicates.
To address the question posed in the title, we have created, and now report details of, an open-access database of cluster structures with a web-assisted interface and toolkit as part of the WASP@N project. The database establishes a map of connectivities within each structure, the information about which is coded and kept as individual labels, called hashkeys, for the nanoclusters. These hashkeys are the basis for structure comparison within the database, and for establishing a map of connectivities between similar structures (topologies). The database is successfully used as a key element in a data-mining study of (MX)12 clusters of three binary compounds (LiI, SrO and GaAs) of which the database has no prior knowledge. The structures are assessed on the energy landscapes determined by the corresponding bulk interatomic potentials. Global optimisation, using a Lamarckian genetic algorithm, is used to search for low lying minima on the same energy landscape to confirm that the data-mined structures form a representative sample of the landscapes, with only very few structures missing from the close energy neighbourhood of the respective global minima.
We have performed a genetic algorithm search on the tight-binding interatomic potential energy surface (PES) for small TiN (N = 2-32) clusters. The low energy candidate clusters were further refined using density functional theory (DFT) calculations with the PBEsol exchange-correlation functional and evaluated with the PBEsol0 hybrid functional. The resulting clusters were analysed in terms of their structural features, growth mechanism and surface area. The results suggest a growth mechanism that is based on forming coordination centres by interpenetrating icosahedra, icositetrahedra and Frank-Kasper polyhedra. We identify centres of coordination, which act as centres of bulk nucleation in medium sized clusters and determine the morphological features of the cluster.
We have developed and implemented a new global optimization technique based on a Lamarckian genetic algorithm with the focus on structure diversity. The key process in the efficient search on a given complex energy landscape proves to be the removal of duplicates that is achieved using a topological analysis of candidate structures. The careful geometrical prescreening of newly formed structures and the introduction of new mutation move classes improve the rate of success further. The power of the developed technique, implemented in the Knowledge Led Master Code, or KLMC, is demonstrated by its ability to locate and explore a challenging double funnel landscape of a Lennard-Jones 38 atom system (LJ38). We apply the redeveloped KLMC to investigate three chemically different systems: ionic semiconductor (ZnO)1-32, metallic Ni13 and covalently bonded C60. All four systems have been systematically explored on the energy landscape defined using interatomic potentials. The new developments allowed us to successfully locate the double funnels of LJ38, find new local and global minima for ZnO clusters, extensively explore the Ni13 and C60 (the buckminsterfullerene, or buckyball) potential energy surfaces.
We employ global optimisation to investigate how oxide nanoclusters of increasing size can best adapt their structure to lower the system energy when interacting with a realistic extended metal support. Specifically, we focus on the (ZnO)@Ag(111) system where experiment has shown that the infinite Ag(111)-supported ZnO monolayer limit corresponds to an epitaxially 7 : 8 matched graphene-like (Zn(3)O(3))-based hexagonal sheet. Using a two-stage search method based on classical interatomic potentials and then on more accurate density functional theory, we report global minina candidate structures for Ag-supported (ZnO)n cluster with sizes ranging from n = 1-24. Comparison with the respective global minina structure of free space (ZnO)n clusters reveals that the surface interaction plays a decisive role in determining the lowest energy Ag-supported (ZnO)n cluster structures. Whereas free space (ZnO)n clusters tend to adopt cage-like bubble structures as they grow larger, Ag-supported (ZnO)n clusters of increasing size become progressively more like planar cuts from the infinite graphene-like ZnO single monolayer. This energetic favourability for planar hexagonal Ag-supported clusters over their 3D counterparts can be partly rationalised by the ZnO-Ag(111) epitaxial matching and the increased number of close interactions with the Ag surface. Detailed analysis shows that this tendency can also be attributed to the capacity of 2D clusters to distort to improve their interaction with the Ag surface relative to more rigid 3D bubble cluster isomers. For the larger sized clusters we find that the adsorption energies and most stable structural types appear to be rather converged confirming that our study makes a bridge between the Ag-supported ZnO monomer and the infinite Ag-supported ZnO monolayer.
Global optimization is used to study the structure of the polar KTaO3 (001) surface. It is found that cation exchange near the surface leads to the most stable structure. This mechanism is likely to be general to metal oxides containing cations of differing charge.
We have investigated the suitability of Time-Dependent Density Functional Theory (TD-DFT) to describe vertical low-energy excitations in naked and hydrated titanium dioxide nanoparticles. Specifically, we compared TD-DFT results obtained using different exchange-correlation (XC) potentials with those calculated using Equation-of-Motion Coupled Cluster (EOM-CC) quantum chemistry methods. We demonstrate that TD-DFT calculations with commonly used XC potentials (e.g., B3LYP) and EOM-CC methods give qualitatively similar results for most TiO2 nanoparticles investigated. More importantly, however, we also show that, for a significant subset of structures, TD-DFT gives qualitatively different results depending upon the XC potential used and that only TD-CAM-B3LYP and TD-BHLYP calculations yield results that are consistent with those obtained using EOM-CC theory. Moreover, we demonstrate that the discrepancies for such structures originate from a particular combination of defects that give rise to charge-transfer excitations, which are poorly described by XC potentials that do not contain sufficient Hartree-Fock like exchange. Finally, we consider that such defects are readily healed in the presence of ubiquitously present water and that, as a result, the description of vertical low-energy excitations for hydrated TiO2 nanoparticles is nonproblematic.
We have investigated the description of excited state relaxation in naked and hydrated TiO2 nanoparticles using Time-Dependent Density Functional Theory (TD-DFT) with three common hybrid exchange-correlation (XC) potentials: B3LYP, CAM-B3LYP and BHLYP. Use of TD-CAM-B3LYP and TD-BHLYP yields qualitatively similar results for all structures, which are also consistent with predictions of coupled-cluster theory for small particles. TD-B3LYP, in contrast, is found to make rather different predictions; including apparent conical intersections for certain particles that are not observed with TD-CAM-B3LYP nor with TD-BHLYP. In line with our previous observations for vertical excitations, the issue with TD-B3LYP appears to be the inherent tendency of TD-B3LYP, and other XC potentials with no or a low percentage of Hartree-Fock like exchange, to spuriously stabilize the energy of charge-transfer (CT) states. Even in the case of hydrated particles, for which vertical excitations are generally well described with all XC potentials, the use of TD-B3LYP appears to result in CT problems during excited state relaxation for certain particles. We hypothesize that the spurious stabilization of CT states by TD-B3LYP even may drive the excited state optimizations to different excited state geometries from those obtained using TD-CAM-B3LYP or TD-BHLYP. Finally, focusing on the TD-CAM-B3LYP and TD-BHLYP results, excited state relaxation in small naked and hydrated TiO2 nanoparticles is predicted to be associated with a large Stokes' shift.
The most widely used oxide for photocatalytic applications owing to its low cost and high activity is TiO2. The discovery of the photolysis of water on the surface of TiO2 in 1972 launched four decades of intensive research into the underlying chemical and physical processes involved. Despite much collected evidence, a thoroughly convincing explanation of why mixed-phase samples of anatase and rutile outperform the individual polymorphs has remained elusive. One long-standing controversy is the energetic alignment of the band edges of the rutile and anatase polymorphs of TiO2 (ref. ). We demonstrate, through a combination of state-of-the-art materials simulation techniques and X-ray photoemission experiments, that a type-II, staggered, band alignment of ~ 0.4 eV exists between anatase and rutile with anatase possessing the higher electron affinity, or work function. Our results help to explain the robust separation of photoexcited charge carriers between the two phases and highlight a route to improved photocatalysts.
Titanium/chemistry , Catalysis , Models, Chemical , Photoelectron Spectroscopy
Al2O3 changes from transparent to a range of intense colours depending on the chemical impurities present. In blue sapphire, Fe and Ti are incorporated; however, the chemical process that gives rise to the colour has long been debated. Atomistic modelling identifies charge transfer from Ti(III) to Fe(III) as being responsible for the characteristic blue appearance.
