DFT based study of transition metal nano-clusters for electrochemical NH3 production.

Theoretical studies of the possibility of producing ammonia electrochemically at ambient temperature and pressure without direct N2 dissociation are presented. Density functional theory calculations were used in combination with the computational standard hydrogen electrode to calculate the free energy profile for the reduction of N2 admolecules and N adatoms on transition metal nanoclusters in contact with an acidic electrolyte. This work has established linear scaling relations for the dissociative reaction intermediates NH, NH2, and NH3. In addition, linear scaling relations for the associative reaction intermediates N2H, N2H2, and N2H3 have been determined. Furthermore, correlations between the adsorption energies of N, N2, and H have been established. These scaling relations and the free energy corrections are used to establish volcanoes describing the onset potential for electrochemical ammonia production and hence describe the potential determining steps for the electrochemical ammonia production. The competing hydrogen evolution reaction has also been analyzed for comparison.

Universal transition state scaling relations for (de)hydrogenation over transition metals.

We analyse the transition state energies for 249 hydrogenation/dehydrogenation reactions of atoms and simple molecules over close-packed and stepped surfaces and nanoparticles of transition metals using Density Functional Theory. Linear energy scaling relations are observed for the transition state structures leading to transition state scaling relations for all the investigated reactions. With a suitable choice of reference systems the transition state scaling relations form a universality class that can be approximated with one single linear relation describing the entire range of reactions over all types of surfaces and nanoclusters.

On the behavior of Brønsted-Evans-Polanyi relations for transition metal oxides.

Versatile Brønsted-Evans-Polanyi (BEP) relations are found from density functional theory for a wide range of transition metal oxides including rutiles and perovskites. For oxides, the relation depends on the type of oxide, the active site, and the dissociating molecule. The slope of the BEP relation is strongly coupled to the adsorbate geometry in the transition state. If it is final state-like the dissociative chemisorption energy can be considered as a descriptor for the dissociation. If it is initial state-like, on the other hand, the dissociative chemisorption energy is not suitable as descriptor for the dissociation. Dissociation of molecules with strong intramolecular bonds belong to the former and molecules with weak intramolecular bonds to the latter group. We show, for the prototype system La-perovskites, that there is a "cyclic" behavior in the transition state characteristics upon change of the active transition metal of the oxide.

Using scaling relations to understand trends in the catalytic activity of transition metals.

A method is developed to estimate the potential energy diagram for a full catalytic reaction for a range of late transition metals on the basis of a calculation (or an experimental determination) for a single metal. The method, which employs scaling relations between adsorption energies, is illustrated by calculating the potential energy diagram for the methanation reaction and ammonia synthesis for 11 different metals on the basis of results calculated for Ru. It is also shown that considering the free energy diagram for the reactions, under typical industrial conditions, provides additional insight into reactivity trends.

Predicting catalysis: understanding ammonia synthesis from first-principles calculations.

Here, we give a full account of a large collaborative effort toward an atomic-scale understanding of modern industrial ammonia production over ruthenium catalysts. We show that overall rates of ammonia production can be determined by applying various levels of theory (including transition state theory with or without tunneling corrections, and quantum dynamics) to a range of relevant elementary reaction steps, such as N(2) dissociation, H(2) dissociation, and hydrogenation of the intermediate reactants. A complete kinetic model based on the most relevant elementary steps can be established for any given point along an industrial reactor, and the kinetic results can be integrated over the catalyst bed to determine the industrial reactor yield. We find that, given the present uncertainties, the rate of ammonia production is well-determined directly from our atomic-scale calculations. Furthermore, our studies provide new insight into several related fields, for instance, gas-phase and electrochemical ammonia synthesis. The success of predicting the outcome of a catalytic reaction from first-principles calculations supports our point of view that, in the future, theory will be a fully integrated tool in the search for the next generation of catalysts.

Towards the computational design of solid catalysts.

Over the past decade the theoretical description of surface reactions has undergone a radical development. Advances in density functional theory mean it is now possible to describe catalytic reactions at surfaces with the detail and accuracy required for computational results to compare favourably with experiments. Theoretical methods can be used to describe surface chemical reactions in detail and to understand variations in catalytic activity from one catalyst to another. Here, we review the first steps towards using computational methods to design new catalysts. Examples include screening for catalysts with increased activity and catalysts with improved selectivity. We discuss how, in the future, such methods may be used to engineer the electronic structure of the active surface by changing its composition and structure.

Cyclic voltammograms for H on Pt(111) and Pt(100) from first principles.

Cyclic voltammetry is a fundamental experimental method for characterizing electrochemical surfaces. Despite its wide use, a way to quantitatively and directly relate cyclic voltammetry to ab initio calculations has been lacking. We derive the cyclic voltammogram for H on Pt(111) and Pt(100), based solely on density functional theory calculations and standard molecular tables. By relating the gas phase adsorption energy to the electrochemical electrode potential, we provide a direct link between surface science and electrochemistry.

Scaling properties of adsorption energies for hydrogen-containing molecules on transition-metal surfaces.

Density functional theory calculations are presented for CHx, x=0,1,2,3, NHx, x=0,1,2, OHx, x=0,1, and SHx, x=0,1 adsorption on a range of close-packed and stepped transition-metal surfaces. We find that the adsorption energy of any of the molecules considered scales approximately with the adsorption energy of the central, C, N, O, or S atom, the scaling constant depending only on x. A model is proposed to understand this behavior. The scaling model is developed into a general framework for estimating the reaction energies for hydrogenation and dehydrogenation reactions.

Combined electronic structure and evolutionary search approach to materials design.

We show that density functional theory calculations have reached an accuracy and speed making it possible to use them in conjunction with an evolutionary algorithm to search for materials with specific properties. The approach is illustrated by finding the most stable four component alloys out of the 192 016 possible fcc and bcc alloys that can be constructed out of 32 different metals. A number of well known and new "super alloys" are identified in this way.