Modelling metal centres, acid sites and reaction mechanisms in microporous catalysts

We discuss the role of QM/MM (embedded cluster) computational techniques in catalytic science, in particular their application to microporous catalysis. We describe the methodologies employed and illustrate their utility by briefly summarising work on metal centres in zeolites. We then report a detailed investigation into the behaviour of methanol at acidic sites in zeolites H-ZSM-5 and H-Y in the context of the methanol-to-hydrocarbons/olefins process. Studying key initial steps of the reaction (the adsorption and subsequent methoxylation), we probe the effect of framework topology and Brønsted acid site location on the energetics of these initial processes. We find that although methoxylation is endothermic with respect to the adsorbed system (by 17-56 kJ mol(-1) depending on the location), there are intriguing correlations between the adsorption/reaction energies and the geometries of the adsorbed species, of particular significance being the coordination of methyl hydrogens. These observations emphasise the importance of adsorbate coordination with the framework in zeolite catalysed conversions, and how this may vary with framework topology and site location, particularly suited to investigation by QM/MM techniques

Segregation effects on the properties of (AuAg)₁₄₇.

AuAg nanoclusters are promising supported co-catalysts for photocatalytic hydrogen reduction. However, beyond the quantum regime (N > 100) little is known about how the electronic properties of these nanoparticles are affected by chemical ordering. We investigate the effects of chemical ordering on the properties of 147-atom cuboctahedral AuAg nanoclusters, using empirical potentials coupled with an atomic-swap basin-hopping search to optimise the elemental distribution, with the lowest energy arrangements then reminimised using Density Functional Theory (DFT). Force-field calculations show Au atoms preferentially occupy sub-surface positions in the bimetallic structures, which results in the formation of a pseudo-onion structure for Ag-rich compositions. At the DFT-level, however, an Ag core surrounded by an Au shell (Ag@Au) is energetically favoured, as electron density can be drawn more readily when Au atoms are positioned on the nanocluster surface, thus resulting in a partial negative charge. Core@shell configurations are analogous to structures that can be chemically synthesised, and further detailed electronic analysis is discussed in the context of nanocluster applications to co-catalysed photocatalysis

Phase control during the synthesis of nickel sulfide nanoparticles from dithiocarbamate precursors

Square-planar nickel bis(dithiocarbamate) complexes, [Ni(S2CNR2)2], have been prepared and utilised as single source precursors to nanoparticulate nickel sulfides. While they are stable in the solid-state to around 300 °C, heating in oleylamine at 230 °C, 5 mM solutions afford pure α-NiS, where the outcome is independent of the substituents. DFT calculations show an electronic effect rather than steric hindrance influences the resulting particle size. Decomposition of the iso-butyl derivative, [Ni(S2CN(i)Bu2)2], has been studied in detail. There is a temperature-dependence of the phase of the nickel sulfide formed. At low temperatures (150 °C), pure α-NiS is formed. Upon raising the temperature, increasing amounts of β-NiS are produced and at 280 °C this is formed in pure form. A range of concentrations (from 5-50 mM) was also investigated at 180 °C and while in all cases pure α-NiS was formed, particle sizes varied significantly. Thus at low concentrations average particle sizes were ca. 100 nm, but at higher concentrations they increased to ca. 150 nm. The addition of two equivalents of tetra-iso-butyl thiuram disulfide, ((i)Bu2NCS2)2, to the decomposition mixture was found to influence the material formed. At 230 °C and above, α-NiS was generated, in contrast to the results found without added thiuram disulfide, suggesting that addition of ((i)Bu2NCS2)2 stabilises the metastable α-NiS phase. At low temperatures (150-180 °C) and concentrations (5 mM), mixtures of α-NiS and Ni3S4, result. A growing proportion of Ni3S4 is noted upon increasing precursor concentration to 10 mM. At 20 mM a metastable phase of nickel sulfide, NiS2 is formed and as the concentration is increased, α-NiS appears alongside NiS2. Reasons for these variations are discussed

Strain and orientation modulated bandgaps and effective masses of phosphorene nanoribbons

Passivated phosphorene nanoribbons, armchair (a-PNR), diagonal (d-PNR), and zigzag (z-PNR), were investigated using density functional theory. Z-PNRs demonstrate the greatest quantum size effect, tuning the bandgap from 1.4 to 2.6 eV when the width is reduced from 26 to 6 Å. Strain effectively tunes charge carrier transport, leading to a sudden increase in electron effective mass at +8% strain for a-PNRs or hole effective mass at +3% strain for z-PNRs, differentiating the (mh*/me*) ratio by an order of magnitude in each case. Straining of d-PNRs results in a direct to indirect band gap transition at either -7% or +5% strain and therein creates degenerate energy valleys with potential applications for valleytronics and/or photocatalysis

Interlayer cation exchange stabilizes polar perovskite surfaces.

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

Effect of nickel monolayer deposition on the structural and electronic properties of the low miller indices of (bcc) iron: A DFT study

Metal clusters of both iron (Fe) and nickel (Ni) have been found in nature as active electro-catalytic sites, for example in the enzyme carbon mono-oxide dehydrogenase found in autotrophic organisms. Thus, surface modification of iron with nickel could improve the surface work function to enhance catalytic applications. The effects of surface modifications of iron by nickel on the structural and electronic properties have been studied using spin-polarised density functional theory calculations within the generalised gradient approximation. The thermodynamically preferred sites for Ni adsorption on the Fe (100), (110) and (111) surfaces have been studied at varying monolayer coverages (including 0.25 ML and 1 ML). The work function of the bare Fe surfaces is found to be of the order (100) ∼ (111) (110) > (100), i.e. 6.78% > 5.76% > 1.99%. The extent of relaxation is magnified on the stepped (111) surface (by 1.09% to 30.88%), where the Ni coordination number is highest at 7 compared to 5 on the (100) facet and 4 on the (110) facet. The Ni deposition changes the work functions of the various surfaces due to charge reordering illustrated by charge density plots, where the work function is reduced only on the (110) surface by 0.04 eV, 0.16 eV and 0.17 eV at 1 ML, 0.5 ML and 0.25 ML respectively, with a concomitant increase in the surface dipole (polarity). This result implies enhanced electron activity and electrochemical reactivity on the most stable and therefore frequently occurring Ni-doped (110) facet compared to the clean (110) facet, which has implications for the development of improved Fe electro-catalysts (for example for CO2 activation and reduction). These findings improve our understanding of the role of surface topology and stability on the extent of Ni interactions with Fe surfaces and the extent to which the Fe surface structures and properties are altered by the Ni deposition

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

Surface energies control the self-organization of oriented In2O3 nanostructures on cubic zirconia.

Highly aligned one-dimensional (1D) nanorods of the transparent conducting oxide In(2)O(3) have been grown on (110)-oriented Y-stabilized ZrO(2) substrates, whereas growth on (100) and (111) substrates leads respectively to blocklike 3D islands and continuous films. It is shown that the striking influence of substrate orientation on the growth morphology is controlled by differences in energies between the low index surfaces of In(2)O(3) and that spontaneous self-organization is driven by minimization of surface energies

CO2 activation and dissociation on the low miller index surfaces of pure and Ni-coated iron metal: a DFT study

We have used spin polarized density functional theory calculations to perform extensive mechanistic studies of CO2 dissociation into CO and O on the clean Fe(100), (110) and (111) surfaces and on the same surfaces coated by a monolayer of nickel. CO2 chemisorbs on all three bare facets and binds more strongly to the stepped (111) surface than on the open flat (100) and close-packed (110) surfaces, with adsorption energies of −88.7 kJ mol−1, −70.8 kJ mol−1 and −116.8 kJ mol−1 on the (100), (110) and (111) facets, respectively. Compared to the bare Fe surfaces, we found weaker binding of the CO2 molecules on the Ni-deposited surfaces, where the adsorption energies are calculated at +47.2 kJ mol−1, −29.5 kJ mol−1 and −65.0 kJ mol−1 on the Ni-deposited (100), (110) and (111) facets respectively. We have also investigated the thermodynamics and activation energies for CO2 dissociation into CO and O on the bare and Ni-deposited surfaces. Generally, we found that the dissociative adsorption states are thermodynamically preferred over molecular adsorption, with the dissociation most favoured thermodynamically on the close-packed (110) facet. The trends in activation energy barriers were observed to follow that of the trends in surface work functions; consequently, the increased surface work functions observed on the Ni-deposited surfaces resulted in increased dissociation barriers and vice versa. These results suggest that measures to lower the surface work function will kinetically promote the dissociation of CO2 into CO and O, although the instability of the activated CO2 on the Ni-covered surfaces will probably result in CO2 desorption from the nickel-doped iron surfaces, as is also seen on the Fe(110) surface