Bronsted-Evans-Polanyi relation of multistep reactions and volcano curve in heterogeneous catalysis

Multistep surface processes involving a number of association reactions and desorption processes may be considered as hypothetical one-step desorption processes. Thus, heterogeneous catalytic reactions can be treated kinetically as consisting of two steps: adsorption and desorption. It is also illustrated that the hypothetical one-step desorption process follows the BEP relation. A volcano curve can be obtained from kinetic analysis by including both adsorption and desorption processes.

Possibility of designing catalysts beyond the traditional volcano curve: A theoretical framework for multi-phase surfaces

The current theory of catalyst activity in heterogeneous catalysis is mainly obtained from the study of catalysts with mono-phases, while most catalysts in real systems consist of multi-phases, the understanding of which is far short of chemists' expectation. Density functional theory (DFT) and micro-kinetics simulations are used to investigate the activities of six mono-phase and nine bi-phase catalysts, using CO hydrogenation that is arguably the most typical reaction in heterogeneous catalysis. Excellent activities that are beyond the activity peak of traditional mono-phase volcano curves are found on some bi-phase surfaces. By analyzing these results, a new framework to understand the unexpected activities of bi-phase surfaces is proposed. Based on the framework, several principles for the design of multi-phase catalysts are suggested. The theoretical framework extends the traditional catalysis theory to understand more complex systems.

Catalyst screening: Refinement of the origin of the volcano curve and its implication in heterogeneous catalysis

Understanding the overall catalytic activity trend for rational catalyst design is one of the core goals in heterogeneous catalysis. In the past two decades, the development of density functional theory (DFT) and surface kinetics make it feasible to theoretically evaluate and predict the catalytic activity variation of catalysts within a descriptor-based framework. Thereinto, the concept of the volcano curve, which reveals the general activity trend, usually constitutes the basic foundation of catalyst screening. However, although it is a widely accepted concept in heterogeneous catalysis, its origin lacks a clear physical picture and definite interpretation. Herein, starting with a brief review of the development of the catalyst screening framework, we use a two-step kinetic model to refine and clarify the origin of the volcano curve with a full analytical analysis by integrating the surface kinetics and the results of first-principles calculations. It is mathematically demonstrated that the volcano curve is an essential property in catalysis, which results from the self-poisoning effect accompanying the catalytic adsorption process. Specifically, when adsorption is strong, it is the rapid decrease of surface free sites rather than the augmentation of energy barriers that inhibits the overall reaction rate and results in the volcano curve. Some interesting points and implications in assisting catalyst screening are also discussed based on the kinetic derivation. Moreover, recent applications of the volcano curve for catalyst design in two important photoelectrocatalytic processes (the hydrogen evolution reaction and dye-sensitized solar cells) are also briefly discussed.

First principles simulation of amorphous silicon bulk, interfaces, and nanowires for photovoltaics

Amorphous silicon has become the material of choice for many technologies, with major applications in large area electronics: displays, image sensing and thin film photovoltaic cells. This technology development has occurred because amorphous silicon is a thin film semiconductor that can be deposited on large, low cost substrates using low temperature. In this thesis, classical molecular dynamics and first principles DFT calculations have been performed to generate structural models of amorphous and hydrogenated amorphous silicon and interfaces of amorphous and crystalline silicon, with the ultimate aim of understanding the photovoltaic properties of core-shell crystalline amorphous Si nanowire structures. We have shown, unexpectedly, from the simulations, that our understanding of hydrogenated bulk a-Si needs to be revisited, with our robust finding that when fully saturated with hydrogen, bulk a-Si exhibits a constant optical energy gap, irrespective of the hydrogen concentration in the sample. Unsaturated a-Si:H, with a lower than optimum hydrogen content, shows a smaller optical gap, that increases with hydrogen content until saturation is reached. The mobility gaps obtained from an analysis of the electronic states show similar behavior. We also obtained that the optical and mobility gaps show a volcano curve as the H content is varied from 7% (undersaturation) to 18% (mild oversaturation). In the case of mild over saturation, the mid-gap states arise exclusively from an increase in the density of strained Si-Si bonds. Analysis of our structures shows the extra H atoms in this case form a bridge between neighboring silicon atoms which increases the corresponding Si-Si distance and promotes bond length disorder in the sample. That has the potential to enhance the Staebler-Wronski effect. Planar interface models of amorphous-crystalline silicon have been generated in Si (100), (110) and (111) surfaces. The interface models are characterized by structure, RDF, electronic density of states and optical absorption spectrum. We find that the least stable (100) surface will result in the formation of the thickest amorphous silicon layer, while the most stable (110) surface forms the smallest amorphous region. We calculated for the first time band offsets of a-Si:H/c-Si heterojunctions from first principles and examined the influence of different surface orientations and amorphous layer thickness on the offsets and implications for device performance. The band offsets depend on the amorphous layer thickness and increase with thickness. By controlling the amorphous layer thickness we can potentially optimise the solar cell parameters. Finally, we have successfully generated different amorphous layer thickness of the a-Si/c-Si and a-Si:H/c-Si 5 nm nanowires from heat and quench. We perform structural analysis of the a-Si-/c-Si nanowires. The RDF, Si-Si bond length distributions, and the coordination number distributions of amorphous regions of the nanowires reproduce similar behaviour compared to bulk amorphous silicon. In the final part of this thesis we examine different surface terminating chemical groups, -H, - OH and –NH2 in (001) GeNW. Our work shows that the diameter of Ge nanowires and the nature of surface terminating groups both play a significant role in both the magnitude and the nature of the nanowire band gaps, allowing tuning of the band gap by up to 1.1 eV. We also show for the first time how the nanowire diameter and surface termination shifts the absorption edge in the Ge nanowires to longer wavelengths. Thus, the combination of nanowire diameter and surface chemistry can be effectively utilised to tune the band gaps and thus light absorption properties of small diameter Ge nanowires.

Understanding the Optimal Adsorption Energies for Catalyst Screening in Heterogeneous Catalysis

The fundamental understanding of the activity in heterogeneous catalysis has long been the major subject in chemistry. This paper shows the development of a two-step model to understand this activity. Using the theory of chemical potential kinetics with Bronsted-Evans-Polanyi relations, the general adsorption energy window is determined from volcano curves, using which the best catalysts can be searched. Significant insights into the reasons for catalytic activity are obtained.

Some Understanding of Fischer-Tropsch Synthesis from Density Functional Theory Calculations

The combination of density functional theory (DFT) calculations and kinetic analyses is a very useful approach to study surface reactions in heterogeneous catalysis. The present paper reviews some recent work applying this approach to Fischer-Tropsch (FT) synthesis. Emphasis is placed on the following fundamental issues in FT synthesis: (i) reactive sites for both hydrogenation and C-C coupling reactions; (ii) reaction mechanisms including carbene mechanism, CO-insertion mechanism and hydroxyl-carbene mechanism; (iii) selectivity with a focus on CH(4) selectivity, alpha-olefin selectivity and chain growth probability; and (iv) activity.

Density Functional Theory Study of Iron and Cobalt Carbides for Fischer-Tropsch Synthesis

Carbides are important phases in heterogeneous catalysis. However, the understanding of carbide phases is inadequate: Fe and Co are the two commercial catalysts for Fischer-Tropsch (FT) synthesis, and experimental work showed that Fe carbide is the active phase in FT synthesis, whereas the appearance of Co carbide is considered as a possible deactivation cause, TO understand very different catalytic roles of carbides, all the key elementary steps in FT synthesis, that is, CO dissociation, C(1) hydrogenation, and C(1)+C(1) coupling, are extensively investigated on both carbide surfaces using first principles calculations. In particular, the most important issues in FT synthesis, the activity and methane selectivity, on the carbide surfaces are quantitatively determined and analyzed. They are also discussed together with metallic Fe and Co surfaces. It is found that (i) Fe carbide is more active than metallic Fe and has similar methane selectivity to Fe, being consistent with the experiments; and (ii) Co carbide is less active than Co and has higher methane selectivity, providing evidence on the molecular level to support the suggestion that the formation of Co carbide is a cause of relatively high methane selectivity and deactivation on Co catalysts.

An understanding and implications of the coverage of surface free sites in heterogeneous catalysis

The crucial roles of the coverage of surface free sites in determining catalytic activity trend are quantitatively addressed with the help of density functional theory and microkinetics. First, by analyzing activity trends of NO oxidation catalyzed by Ru, Rh, Pd, Os, Ir, and Pt surfaces with full kinetic considerations, we identify that the activity trend is in general determined by the competition between the reaction barrier and the coverage of surface free sites. Second, since the dissociation of many important molecules, such as the dissociation of N(2), O(2), and CO, follows the same Bronsted-Evans-Polanyi relationship, the coverage of surface free sites is usually a decisive term that affects the overall activity. Third, an equation is derived for the coverage of surface free sites and it is found that the coverage of surface free sites contains not only all the key thermodynamic parameters but also all the kinetic properties in the catalytic system. (C) 2009 American Institute of Physics. [DOI: 10.1063/1.3140202]

A first principles analysis of the effect of hydrogen concentration in hydrogenated amorphous silicon on the formation of strained Si-Si bonds and the optical and mobility gaps

In this paper, we use a model of hydrogenated amorphous silicon generated from molecular dynamics with density functional theory calculations to examine how the atomic geometry and the optical and mobility gaps are influenced by mild hydrogen oversaturation. The optical and mobility gaps show a volcano curve as the hydrogen content varies from undersaturation to mild oversaturation, with largest gaps obtained at the saturation hydrogen concentration. At the same time, mid-gap states associated with dangling bonds and strained Si-Si bonds disappear at saturation but reappear at mild oversaturation, which is consistent with the evolution of optical gap. The distribution of Si-Si bond distances provides the key to the change in electronic properties. In the undersaturation regime, the new electronic states in the gap arise from the presence of dangling bonds and strained Si-Si bonds, which are longer than the equilibrium Si-Si distance. Increasing hydrogen concentration up to saturation reduces the strained bonds and removes dangling bonds. In the case of mild oversaturation, the mid-gap states arise exclusively from an increase in the density of strained Si-Si bonds. Analysis of our structure shows that the extra hydrogen atoms form a bridge between neighbouring silicon atoms, thus increasing the Si-Si distance and increasing disorder in the sample.

Elastic Propagation in random media: applications to the imaging of volcano structures

High-frequency seismograms contain features that reflect the random inhomogeneities of the earth. In this work I use an imaging method to locate the high contrast small- scale heterogeneity respect to the background earth medium. This method was first introduced by Nishigami (1991) and than applied to different volcanic and tectonically active areas (Nishigami, 1997, Nishigami, 2000, Nishigami, 2006). The scattering imaging method is applied to two volcanic areas: Campi Flegrei and Mt. Vesuvius. Volcanic and seismological active areas are often characterized by complex velocity structures, due to the presence of rocks with different elastic properties. I introduce some modifications to the original method in order to make it suitable for small and highly complex media. In particular, for very complex media the single scattering approximation assumed by Nishigami (1991) is not applicable as the mean free path becomes short. The multiple scattering or diffusive approximation become closer to the reality. In this thesis, differently from the ordinary Nishigami’s method (Nishigami, 1991), I use the mean of the recorded coda envelope as reference curve and calculate the variations from this average envelope. In this way I implicitly do not assume any particular scattering regime for the "average" scattered radiation, whereas I consider the variations as due to waves that are singularly scattered from the strongest heterogeneities. The imaging method is applied to a relatively small area (20 x 20 km), this choice being justified by the small length of the analyzed codas of the low magnitude earthquakes. I apply the unmodified Nishigami’s method to the volcanic area of Campi Flegrei and compare the results with the other tomographies done in the same area. The scattering images, obtained with frequency waves around 18 Hz, show the presence of high scatterers in correspondence with the submerged caldera rim in the southern part of the Pozzuoli bay. Strong scattering is also found below the Solfatara crater, characterized by the presence of densely fractured, fluid-filled rocks and by a strong thermal anomaly. The modified Nishigami’s technique is applied to the Mt. Vesuvius area. Results show a low scattering area just below the central cone and a high scattering area around it. The high scattering zone seems to be due to the contrast between the high rigidity body located beneath the crater and the low rigidity materials located around it. The central low scattering area overlaps the hydrothermal reservoirs located below the central cone. An interpretation of the results in terms of geological properties of the medium is also supplied, aiming to find a correspondence of the scattering properties and the geological nature of the material. A complementary result reported in this thesis is that the strong heterogeneity of the volcanic medium create a phenomenon called "coda localization". It has been verified that the shape of the seismograms recorded from the stations located at the top of the volcanic edifice of Mt. Vesuvius is different from the shape of the seismograms recorded at the bottom. This behavior is justified by the consideration that the coda energy is not uniformly distributed within a region surrounding the source for great lapse time.