A computational study of direct CO₂ 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

Mechanistic and kinetic investigation on maximizing the formation of levoglucosan from cellulose during biomass pyrolysis

Attempts to understand and control the formation of levoglucosan, as a key bio-oil compound and an important intermediate chemical, during biomass pyrolysis have recently generated a large number of experimental and computational studies. Whilst promising mechanisms have been put forward to explain levoglucosan formation, there has yet to be clarity on the factors that inhibit the stoichiometric yields from cellulose, and, current available mechanisms adopt model compounds of short chain oligosaccharide structures rather than cellulose polymer. In this work, kinetic models describing cellulose (with Degree of Polymerisation DP of 2048) decomposition and levoglucosan formation have been constructed. It was found that temperature and residence time co-ordinately influenced the yield of levoglucosan, i.e. maximum yield can be reached within 80 s at 600 °C, however, around 10 min is required at 550 °C. Anhydro-oligosaccharides greater than DP8 will be generated and then consumed in 60 s with pyrolysis temperatures of 500 °C and over. The results also showed that the current available concerted initiation and depropagation reactions are insufficient to explain the experimentally observed high rates of cellulose depolymerisation and levoglucosan formation, other key factors such as side reactions have to be considered. Additionally, it was found that the rate of levoglucosan generation is independent of the initial cellulose chain length, which is a finding that will facilitate future studies in the area of bioenergy and biomass decomposition

Investigation of the Pd(1−x)Znx alloy phase diagram using ab initio modelling approaches

The identification of the stable phases in alloy materials is challenging because composition affects the structural stability of different intermediate phases. Computational simulation, via multiscale modelling approaches, can significantly accelerate the exploration of phase space and help to identify stable phases. Here, we apply such new approaches to understand the complex phase diagram of binary alloys of PdZn, with the relative stability of structural polymorphs considered through application of density functional theory coupled with cluster expansion (CE). The experimental phase diagram has several competing crystal structures, and we focus on three different closed-packed phases that are commonly observed for PdZn, namely the face-centred cubic (FCC), body-centred tetragonal (BCT) and hexagonal close packed (HCP), to identify their respective stability ranges. Our multiscale approach confirms a narrow range of stability for the BCT mixed alloy, within the Zn concentration range from 43.75% to 50%, which aligns with experimental observations. We subsequently use CE to show that the phases are competitive across all concentrations, but with the FCC alloy phase favoured for Zn concentrations below 43.75%, and that the HCP structure favoured for Zn-rich concentrations. Our methodology and results provide a platform for future investigations of PdZn and other close-packed alloy systems with multiscale modelling techniques

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

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

Investigation of the Pd (1− x ) Zn x alloy phase diagram using ab initio modelling approaches

The identification of the stable phases in alloy materials is challenging because composition affects the structural stability of different intermediate phases. Computational simulation, via multiscale modelling approaches, can significantly accelerate the exploration of phase space and help to identify stable phases. Here, we apply such new approaches to understand the complex phase diagram of binary alloys of PdZn, with the relative stability of structural polymorphs considered through application of density functional theory coupled with cluster expansion (CE). The experimental phase diagram has several competing crystal structures, and we focus on three different closed-packed phases that are commonly observed for PdZn, namely the face-centred cubic (FCC), body-centred tetragonal (BCT) and hexagonal close packed (HCP), to identify their respective stability ranges. Our multiscale approach confirms a narrow range of stability for the BCT mixed alloy, within the Zn concentration range from 43.75% to 50%, which aligns with experimental observations. We subsequently use CE to show that the phases are competitive across all concentrations, but with the FCC alloy phase favoured for Zn concentrations below 43.75%, and that the HCP structure favoured for Zn-rich concentrations. Our methodology and results provide a platform for future investigations of PdZn and other close-packed alloy systems with multiscale modelling techniques

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

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

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 shows 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, while the Pd2Ga alloy formed for the Pd/Ga2O3 sample. The heats of mixing were calculated for a variety of alloy compositions with high heats of mixing calculated for both PdZn and Pd2Ga alloys, with values of ca. -0.6 eV/atom and ca. -0.8 eV/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

The critical role of βPdZn alloy in Pd/ZnO catalysts for the hydrogenation of carbon dioxide to methanol

The rise in atmospheric CO2 concentration and the concomitant rise in global surface temperature have prompted massive research effort in designing catalytic routes to utilize CO2 as a feedstock. Prime among these is the hydrogenation of CO2 to make methanol, which is a key commodity chemical intermediate, a hydrogen storage molecule, and a possible future fuel for transport sectors that cannot be electrified. Pd/ZnO has been identified as an effective candidate as a catalyst for this reaction, yet there has been no attempt to gain a fundamental understanding of how this catalyst works and more importantly to establish specific design criteria for CO2 hydrogenation catalysts. Here, we show that Pd/ZnO catalysts have the same metal particle composition, irrespective of the different synthesis procedures and types of ZnO used here. We demonstrate that all of these Pd/ZnO catalysts exhibit the same activity trend. In all cases, the β-PdZn 1:1 alloy is produced and dictates the catalysis. This conclusion is further supported by the relationship between conversion and selectivity and their small variation with ZnO surface area in the range 6–80 m2g–1. Without alloying with Zn, Pd is a reverse water-gas shift catalyst and when supported on alumina and silica is much less active for CO2 conversion to methanol than on ZnO. Our approach is applicable to the discovery and design of improved catalysts for CO2 hydrogenation and will aid future catalyst discovery