Environment-Driven Variability in Absolute Band Edge Positions and Work Functions of Reduced Ceria.

The absolute band edge positions and work function (Φ) are the key electronic properties of metal oxides that determine their performance in electronic devices and photocatalysis. However, experimental measurements of these properties often show notable variations, and the mechanisms underlying these discrepancies remain inadequately understood. In this work, we focus on ceria (CeO2), a material renowned for its outstanding oxygen storage capacity, and combine theoretical and experimental techniques to demonstrate environmental modifications of its ionization potential (IP) and Φ. Under O-deficient conditions, reduced ceria exhibits a decreased IP and Φ with significant sensitivity to defect distributions. In contrast, the IP and Φ are elevated in O-rich conditions due to the formation of surface peroxide species. Surface adsorbates and impurities can further augment these variabilities under realistic conditions. We rationalize the shifts in energy levels by separating the individual contributions from bulk and surface factors, using hybrid quantum mechanical/molecular mechanical (QM/MM) embedded-cluster and periodic density functional theory (DFT) calculations supported by interatomic-potential-based electrostatic analyses. Our results highlight the critical role of on-site electrostatic potentials in determining the absolute energy levels in metal oxides, implying a dynamic evolution of band edges under catalytic conditions.

Heterogeneous Trimetallic Nanoparticles as Catalysts.

The development and application of trimetallic nanoparticles continues to accelerate rapidly as a result of advances in materials design, synthetic control, and reaction characterization. Following the technological successes of multicomponent materials in automotive exhausts and photovoltaics, synergistic effects are now accessible through the careful preparation of multielement particles, presenting exciting opportunities in the field of catalysis. In this review, we explore the methods currently used in the design, synthesis, analysis, and application of trimetallic nanoparticles across both the experimental and computational realms and provide a critical perspective on the emergent field of trimetallic nanocatalysts. Trimetallic nanoparticles are typically supported on high-surface-area metal oxides for catalytic applications, synthesized via preparative conditions that are comparable to those applied for mono- and bimetallic nanoparticles. However, controlled elemental segregation and subsequent characterization remain challenging because of the heterogeneous nature of the systems. The multielement composition exhibits beneficial synergy for important oxidation, dehydrogenation, and hydrogenation reactions; in some cases, this is realized through higher selectivity, while activity improvements are also observed. However, challenges related to identifying and harnessing influential characteristics for maximum productivity remain. Computation provides support for the experimental endeavors, for example in electrocatalysis, and a clear need is identified for the marriage of simulation, with respect to both combinatorial element screening and optimal reaction design, to experiment in order to maximize productivity from this nascent field. Clear challenges remain with respect to identifying, making, and applying trimetallic catalysts efficiently, but the foundations are now visible, and the outlook is strong for this exciting chemical field.

Influence of Solvent on Selective Catalytic Reduction of Nitrogen Oxides with Ammonia over Cu-CHA Zeolite.

The copper-exchanged zeolite Cu-CHA has received considerable attention in recent years, owing to its application in the selective catalytic reduction (SCR) of NOx species. Here, we study the NH3-SCR reaction mechanism on Cu-CHA using the hybrid quantum mechanical/molecular mechanical (QM/MM) technique and investigate the effects of solvent on the reactivity of active Cu species. To this end, a comparison is made between water- and ammonia-solvated and bare Cu species. The results show the promoting effect of solvent on the oxidation component of the NH3-SCR cycle since the formation of important nitrate species is found to be energetically more favorable on the solvated Cu sites than in the absence of solvent molecules. Conversely, both solvent molecules are predicted to inhibit the reduction component of the NH3-SCR cycle. Diffuse reflectance infrared fourier-transform spectroscopy (DRIFTS) experiments exploiting (concentration) modulation excitation spectroscopy (MES) and phase-sensitive detection (PSD) identified spectroscopic signatures of Cu-nitrate and Cu-nitrosamine (H2NNO), important species which had not been previously observed experimentally. This is further supported by the QM/MM-calculated harmonic vibrational analysis. Additional insights are provided into the reactivity of solvated active sites and the formation of key intermediates including their formation energies and vibrational spectroscopic signatures, allowing the development of a detailed understanding of the reaction mechanism. We demonstrate the role of solvated active sites and their influence on the energetics of important species that must be explicitly considered for an accurate understanding of NH3-SCR kinetics.

Experimental and theoretical investigations on the anti-perovskite nitrides Co3CuN, Ni3CuN and Co3MoN for ammonia synthesis.

The ammonia synthesis activities of the anti-perovskite nitrides Co3CuN and Ni3CuN have been compared to investigate the possible metal composition-activity relationship. Post-reaction elemental analysis showed that the activity for both nitrides was due to loss of lattice nitrogen rather than a catalytic process. Co3CuN was observed to convert a higher percentage of lattice nitrogen to ammonia than Ni3CuN and was active at a lower temperature. The loss of lattice nitrogen was revealed to be topotactic and Co3Cu and Ni3Cu were formed during the reaction. Therefore, the anti-perovskite nitrides may be of interest as reagents for the formation of ammonia through chemical looping. The regeneration of the nitrides was achieved by ammonolysis of the corresponding metal alloys. However, regeneration using N2 was shown to be challenging. In order to understand the difference in reactivity between the two nitrides, DFT techniques were applied to investigate the thermodynamics of the processes involved in the evolution of lattice nitrogen to the gas phase via conversion to N2 or NH3, revealing key differences in the energetics of bulk conversion of the anti-perovskite to the alloy phase, and in loss of surface N from the stable low-index N-terminated (111) and (100) facets. Computational modelling of the density of states (DOS) at the Fermi level was performed. It was shown that the Ni and Co d states contributed to the density of states and that the Cu d states only contributed to the DOS for Co3CuN. The anti-perovskite Co3MoN has been investigated as comparisons with Co3Mo3N may give an insight into the role structure type plays in the ammonia synthesis activity. The XRD pattern and elemental analysis for the synthesised material revealed that an amorphous phase was present that contained nitrogen. In contrast to Co3CuN and Ni3CuN, the material was shown to have steady state activity at 400 °C with a rate of 92 ± 15 µmol h-1 g-1. Therefore, it appears that metal composition has an influence on the stability and activity of the anti-perovskite nitrides.

Mechanism of ammonia synthesis on Fe3Mo3N.

Ammonia (NH3) synthesis is an essential yet energy-demanding industrial process. Hence, there is a need to develop NH3 synthesis catalysts that are highly active under milder conditions. Metal nitrides are promising candidates, with the η-carbide Co3Mo3N having been found to be more active than the industrial Fe-based catalyst. The isostructural Fe3Mo3N catalyst has also been identified as highly active for NH3 synthesis. In the present work, we investigate the catalytic ammonia synthesis mechanisms in Fe3Mo3N, which we compare and contrast with the previously studied Co3Mo3N. We apply plane-wave density functional theory (DFT) to investigate surface N vacancy formation in Fe3Mo3N, and two distinct ammonia synthesis mechanisms. The calculations reveal that whilst N vacancy formation on Fe3Mo3N is more thermodynamically demanding than for Co3Mo3N, the formation energies are comparable, suggesting that surface lattice N vacancies in Fe3Mo3N could facilitate NH3 synthesis. N2 activation was found to be enhanced on Fe3Mo3N compared to Co3Mo3N, for adsorption both at and adjacent to the vacancy. The calculated activation barriers suggest that, as for Co3Mo3N, the associative Mars van Krevelen mechanism affords a much less energy-demanding pathway for ammonia synthesis, especially for initial hydrogenation processes.

Methanol synthesis from CO2 and H2 using supported Pd alloy catalysts.

A number of Pd based materials have been synthesised and evaluated as catalysts for the conversion of carbon dioxide and hydrogen to methanol, a useful platform chemical and hydrogen storage molecule. Monometallic Pd catalysts show poor methanol selectivity, but this is improved through the formation of Pd alloys, with both PdZn and PdGa alloys showing greatly enhanced methanol productivity compared with monometallic Pd/Al2O3 and Pd/TiO2 catalysts. Catalyst characterisation shows that the 1 : 1 ß-PdZn alloy is present in all Zn containing post-reaction samples, including PdZn/Ga2O3, with the Pd2Ga alloy formed for the Pd/Ga2O3 sample. The heat of mixing was calculated for a variety of alloy compositions with high values determined for both PdZn and Pd2Ga alloys, at ca. -0.6 eV per atom and ca. -0.8 eV per atom, respectively. However, ZnO is more readily reduced than Ga2O3, providing a possible explanation for the preferential formation of the PdZn alloy, rather than PdGa, when in the presence of Ga2O3.

Exploring the length scales, timescales and chemistry of challenging materials (Part 1).

This themed issue explores the different length scales and timescales that determine the physics and chemistry of a variety of key materials, explored from the perspective of a wide range of disciplines, including physics, chemistry, materials science, Earth science and biochemistry. The topics discussed include catalysis, chemistry under extreme conditions, energy materials, amorphous and liquid structure, hybrid organic materials and biological materials. The issue is in two parts, with the present part exploring glassy and amorphous systems and materials at high pressure. This article is part of the theme issue 'Exploring the length scales, timescales and chemistry of challenging materials (Part 1)'.

Methanol diffusion and dynamics in zeolite H-ZSM-5 probed by quasi-elastic neutron scattering and classical molecular dynamics simulations.

Zeolite ZSM-5 is a key catalyst in commercially relevant processes including the widely studied methanol to hydrocarbon reaction, and molecular diffusion in zeolite pores is known to be a crucial factor in controlling catalytic reactions. Here, we present critical analyses of recent quasi-elastic neutron scattering (QENS) data and complementary molecular dynamics (MD) simulations. The QENS experiments show that the nature of methanol diffusion dynamics in ZSM-5 pores is dependent both on the Si/Al ratio (11, 25, 36, 40 and 140), which determines the Brønsted acid site density of the zeolite, and that the nature of the type of motion observed may vary qualitatively over a relatively small temperature range. At 373 K, on increasing the ratio from 11 to 140, the observed mobile methanol fraction increases and the nature of methanol dynamics changes from rotational (in ZSM-5 with Si/Al of 11) to translational diffusion. The latter is either confined localized diffusion within a pore in zeolites with ratios up to 40 or non-localized, longer-range diffusion in zeolite samples with the ratio of 140. The complementary MD simulations conducted over long time scales (1 ns), which are longer than those measured in the present study by QENS (≈1-440 ps), at 373 K predict the occurrence of long-range translational diffusion of methanol in ZSM-5, independent of the Si/Al ratios (15, 47, 95, 191 and siliceous MFI). The rate of diffusion increases slightly by increasing the ratio from 15 to 95 and thereafter does not depend on zeolite composition. Discrepancies in the observed mobile methanol fraction between the MD simulations (100% methanol mobility in ZSM-5 pores across all Si/Al ratios) and QENS experiments (for example, ≈80% immobile methanol in ZSM-5 with Si/Al of 11) are attributed to the differences in time resolutions of the techniques. This perspective provides comprehensive information on the effect of acid site density on methanol dynamics in ZSM-5 pores and highlights the complementarity of QENS and MD, and their advantages and limitations. This article is part of the theme issue 'Exploring the length scales, timescales and chemistry of challenging materials (Part 2)'.

Computational infrared and Raman spectra by hybrid QM/MM techniques: a study on molecular and catalytic material systems.

Vibrational spectroscopy is one of the most well-established and important techniques for characterizing chemical systems. To aid the interpretation of experimental infrared and Raman spectra, we report on recent theoretical developments in the ChemShell computational chemistry environment for modelling vibrational signatures. The hybrid quantum mechanical and molecular mechanical approach is employed, using density functional theory for the electronic structure calculations and classical forcefields for the environment. Computational vibrational intensities at chemical active sites are reported using electrostatic and fully polarizable embedding environments to achieve more realistic vibrational signatures for materials and molecular systems, including solvated molecules, proteins, zeolites and metal oxide surfaces, providing useful insight into the effect of the chemical environment on the signatures obtained from experiment. This work has been enabled by the efficient task-farming parallelism implemented in ChemShell for high-performance computing platforms.  This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.

The effect of dissolved chlorides on the photocatalytic degradation properties of titania in wastewater treatment.

We investigate the effect of chlorides on the photocatalytic degradation of phenol by titania polymorphs (anatase and rutile). We demonstrate how solubilised chlorides can affect the hydroxyl radical formation on both polymorphs with an overall effect on their photodegradative activity. Initially, the photocatalytic activity of anatase and rutile for phenol degradation is investigated in both standard water and brines. With anatase, a significant reduction of the phenol conversion rate is observed (from a pseudo-first-order rate constant k = 5.3 × 10-3 min-1 to k = 3.5 × 10-3 min-1). In contrast, the presence of solubilised chlorides results in enhancement of rutile activity under the same reaction conditions (from 2.3 × 10-3 min-1 to 4.8 × 10-3 min-1). Periodic DFT methods are extensively employed and we show that after the generation of charge separation in the modelled titania systems, adsorbed chlorides are the preferential site for partial hole localisation, although small energy differences are computed between partially localised hole densities over adsorbed chloride or hydroxyl. Moreover, chlorides can reduce or inhibit the ability of r-TiO2 (110) and a-TiO2 (101) systems to localise polarons in the slab structure. These results indicate that both mechanisms - hole scavenging and the inhibition of hole localisation - can be the origin of the effect of chlorides on photocatalytic activity of both titania polymorphs. These results provide fundamental insight into the photocatalytic properties of titania polymorphs and elucidate the effect of adsorbed anions over radical formation and oxidative decomposition of organic pollutants.

Multiscale QM/MM modelling of catalytic systems with ChemShell.

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.

Cd/Pt Precursor Solution for Solar H2 Production and in situ Photochemical Synthesis of Pt Single-atom Decorated CdS Nanoparticles.

Despite extensive efforts to develop high-performance H2 evolution catalysts, this remains a major challenge. Here, we demonstrate the use of Cd/Pt precursor solutions for significant photocatalytic H2 production (154.7 mmol g-1 h-1 ), removing the need for a pre-synthesized photocatalyst. In addition, we also report simultaneous in situ synthesis of Pt single-atoms anchored CdS nanoparticles (PtSA -CdSIS ) during photoirradiation. The highly dispersed in situ incorporation of extensive Pt single atoms on CdSIS enables the enhancement of active sites and suppresses charge recombination, which results in exceptionally high solar-to-hydrogen conversion efficiency of ≈1 % and an apparent quantum yield of over 91 % (365 nm) for H2 production. Our work not only provides a promising strategy for maximising H2 production efficiency but also provides a green process for H2 production and the synthesis of highly photoactive PtSA -CdSIS nanoparticles.

Bulk and Surface Contributions to Ionisation Potentials of Metal Oxides.

Determining the absolute band edge positions in solid materials is crucial for optimising their performance in wide-ranging applications including photocatalysis and electronic devices. However, obtaining absolute energies is challenging, as seen in CeO2 , where experimental measurements show substantial discrepancies in the ionisation potential (IP). Here, we have combined several theoretical approaches, from classical electrostatics to quantum mechanics, to elucidate the bulk and surface contributions to the IP of metal oxides. We have determined a theoretical bulk contribution to the IP of stoichiometric CeO2 of only 5.38 eV, while surface orientation results in intrinsic IP variations ranging from 4.2 eV to 8.2 eV. Highly tuneable IPs were also found in TiO2 , ZrO2 , and HfO2 , in which surface polarisation plays a pivotal role in long-range energy level shifting. Our analysis, in addition to rationalising the observed range of experimental results, provides a firm basis for future interpretations of experimental and computational studies of oxide band structures.

A computational study of direct CO2 hydrogenation to methanol on Pd surfaces.

The reaction mechanism of direct CO2 hydrogenation to methanol is investigated in detail on Pd (111), (100) and (110) surfaces using density functional theory (DFT), supporting investigations into emergent Pd-based catalysts. Hydrogen adsorption and surface mobility are firstly considered, with high-coordination surface sites having the largest adsorption energy and being connected by diffusion channels with low energy barriers. Surface chemisorption of CO2, forming a partially charged CO2δ-, is weakly endothermic on a Pd (111) whilst slightly exothermic on Pd (100) and (110), with adsorption enthalpies of 0.09, -0.09 and -0.19 eV, respectively; the low stability of CO2δ- on the Pd (111) surface is attributed to negative charge accumulating on the surface Pd atoms that interact directly with the CO2δ- adsorbate. Detailed consideration for sequential hydrogenation of the CO2 shows that HCOOH hydrogenation to H2COOH would be the rate determining step in the conversion to methanol, for all surfaces, with activation barriers of 1.41, 1.51, and 0.84 eV on Pd (111), (100) and (110) facets, respectively. The Pd (110) surface exhibits overall lower activation energies than the most studied Pd (111) and (100) surfaces, and therefore should be considered in more detail in future Pd catalytic studies.

Concluding remarks: Reaction mechanisms in catalysis: perspectives and prospects.

We consider the current status of our understanding of reaction mechanisms in catalysis in the light of the papers presented in this Discussion. We identify some of the challenges in both theoretical and experimental studies, which we illustrate by considering three key reactions.

How bulk and surface properties of Ti4SiC3, V4SiC3, Nb4SiC3 and Zr4SiC3 tune reactivity: a computational study.

We present several in silico insights into the MAX-phase of early transition metal silicon carbides and explore how these affect carbon dioxide hydrogenation. Periodic density functional methodology is applied to models of Ti4SiC3, V4SiC3, Nb4SiC3 and Zr4SiC3. We find that silicon and carbon terminations are unstable, with sintering occurring in vacuum and significant reconstruction taking place under an oxidising environment. In contrast, the metal terminated surfaces are highly stable and very active towards CO2 reduction. However, we predict that under reaction conditions these surfaces are likely to be oxidised. These results are compared to studies on comparable materials and we predict optimal values for hydrogen evolution and CO2 reduction.

Combination of theoretical and in situ experimental investigations of the role of lithium dopant in manganese nitride: a two-stage reagent for ammonia synthesis.

Manganese nitride related materials are of interest as two-stage reagents for ammonia synthesis via nitrogen chemical looping. However, unless they are doped with a co-cation, manganese nitrides are thermochemically stable and a high temperature is required to produce ammonia under reducing conditions, thereby hindering their use as nitrogen transfer materials. Nevertheless, when lithium is used as dopant, ammonia generation can be observed at a reaction temperature as low as 300 °C. In order to develop strategies for the improvement of the reactivity of nitride materials in the context of two-stage reagents, it is necessary to understand the intrinsic role of the dopant in the mechanism of ammonia synthesis. To this end, we have investigated the role of lithium in increasing the manganese nitride reactivity by in situ neutron diffraction studies and N2 and H2 isotopic exchange reactions supplemented by DFT calculations.

A combined periodic DFT and QM/MM approach to understand the radical mechanism of the catalytic production of methanol from glycerol.

The production of methanol from glycerol over a basic oxide, such as MgO, using high reaction temperatures (320 °C) is a promising new approach to improving atom efficiency in the production of biofuels. The mechanism of this reaction involves the homolytic cleavage of the C3 feedstock, or its dehydration product hydroxyacetone, to produce a hydroxymethyl radical species which can then abstract an H atom from other species. Obtaining a detailed reaction mechanism for this type of chemistry is difficult due to the large number of products present when the system is operated at high conversions. In this contribution we show how DFT based modelling studies can provide new insights into likely reaction pathways, in particular the source of H atoms for the final step of converting hydroxymethyl radicals to methanol. We show that water is unlikely to be important in this stage of the process, C-H bonds of C2 and C3 species can give an energetically favourable pathway and that the disproportionation of hydroxymethyl radicals to methanol and formaldehyde produces a very favourable route. Experimental analysis of reaction products confirms the presence of formaldehyde. The calculations presented in this work also provide new insight into the role of the catalyst surface in the reaction showing that the base sites of the MgO(100) are able to deprotonate hydroxymethyl radicals but not methanol itself. In carrying out the calculations we also show how periodic DFT and QM/MM approaches can be used together to obtain a rounded picture of molecular adsorption to surfaces and homolytic bond cleavage which are both central to the reactions studied.

QM/MM study of the stability of dimethyl ether in zeolites H-ZSM-5 and H-Y.

The methanol-to-hydrocarbons (MTH) process transforms C1 carbon sources to higher hydrocarbons, but details of the mechanism that leads to the formation of the first carbon-carbon bond remain unclear. Here, we present a computational investigation of how a crucial intermediate, dimethyl ether (DME), interacts with different zeolite catalysts (H-ZSM-5, H-Y) to gain insight into the initial stages in the MTH process. We use QM/MM computational simulations to model the conversion of methanol to DME in H-ZSM-5, which is a well characterised and important reaction intermediate. We analyse and compare the stability of DME on several acid sites in H-ZSM-5 and H-Y, and show that the more acidic and open "intersection sites" in the H-ZSM-5 framework are able to bond strongest with DME, with complete deprotonation of the acid site occurring. The conversion of methanol to DME in H-ZSM-5 is calculated as requiring a higher activation energy than framework methoxylation, which indicates that a stepwise (indirect) mechanism, through a methoxy intermediate, is the most likely route to DME formation during the initiation of the MTH process.

A computational study of the properties of low- and high-index Pd, Cu and Zn surfaces.

We report a detailed Density Functional Theory (DFT) based investigation of the structure and stability of bulk and surface structures for the Group 10-12 elements Pd, Cu and Zn, considering the effect of the choice of exchange-correlation density functional and computation parameters. For the initial bulk structures, the lattice parameter and cohesive energy are calculated, which are then augmented by calculation of surface energies and work functions for the lower-index surfaces. Of the 22 density functionals considered, we highlight the mBEEF density functional as providing the best overall agreement with experimental data. The optimal density functional choice is applied to the study of higher index surfaces for the three metals, and Wulff constructions performed for nanoparticles with a radius of 11 nm, commensurate with nanoparticle sizes commonly employed in catalytic chemistry. For Pd and Cu, the low-index (111) facet is dominant in the constructed nanoparticles, covering ∼50% of the surface, with (100) facets covering a further 10 to 25%; however, non-negligible coverage from higher index (332), (332) and (210) facets is also observed for Pd, and (322), (221) and (210) surfaces are observed for Cu. In contrast, only the (0001) and (10-10) facets are observed for Zn. Overall, our results highlight the need for careful validation of computational settings before performing extensive density functional theory investigations of surface properties and nanoparticle structures of metals.