-
Efficient Pourbaix diagrams of many-element compounds
Authors:
Anjli Patel,
Jens K. Nørskov,
Kristin A. Persson,
Joseph H. Montoya
Abstract:
Pourbaix diagrams have long been an essential tool for determining the phase stability of solids and their associated ionic species under electrochemical conditions. In recent years, Pourbaix diagrams have been used for applications ranging from corrosion-resistance alloy design to electrocatalysis, and the data from which they are generated has been enhanced by the availability of materials data…
▽ More
Pourbaix diagrams have long been an essential tool for determining the phase stability of solids and their associated ionic species under electrochemical conditions. In recent years, Pourbaix diagrams have been used for applications ranging from corrosion-resistance alloy design to electrocatalysis, and the data from which they are generated has been enhanced by the availability of materials data in various online databases. However, generation of multi-element Pourbaix diagrams has a critical bottleneck which makes 3-element systems difficult to analyze quickly, and 4 and 5 element systems intractable. In this work, we present a method for constructing Pourbaix diagrams which uses a pre-processing step to circumvent the most egregious computational bottleneck and make many-element Pourbaix diagrams computationally efficient.
△ Less
Submitted 30 August, 2019;
originally announced September 2019.
-
Selective electrochemical generation of hydrogen peroxide from water oxidation
Authors:
Venkatasubramanian Viswanathan,
Heine A. Hansen,
Jens K. Nørskov
Abstract:
Water is a life-giving source, fundamental to human existence, yet, over a billion people lack access to clean drinking water. Present techniques for water treatment such as piped, treated water rely on time and resource intensive centralized solutions. In this work, we propose a decentralized device concept that can utilize sunlight to split water into hydrogen and hydrogen peroxide. The hydrogen…
▽ More
Water is a life-giving source, fundamental to human existence, yet, over a billion people lack access to clean drinking water. Present techniques for water treatment such as piped, treated water rely on time and resource intensive centralized solutions. In this work, we propose a decentralized device concept that can utilize sunlight to split water into hydrogen and hydrogen peroxide. The hydrogen peroxide can oxidize organics while the hydrogen bubbles out. In enabling this device, we require an electrocatalyst that can oxidize water while suppressing the thermodynamically favored oxygen evolution and form hydrogen peroxide. Using density functional theory calculations, we show that the free energy of adsorbed OH$^*$ can be used as a descriptor to screen for selectivity trends between the 2e$^-$ water oxidation to H$_2$O$_2$ and the 4e$^-$ oxidation to O$_2$. We show that materials that bind oxygen intermediates sufficiently weakly, such as SnO$_2$, can activate hydrogen peroxide evolution. We present a rational design principle for the selectivity in electrochemical water oxidation and identify new material candidates that could perform H$_2$O$_2$ evolution selectively.
△ Less
Submitted 24 September, 2015;
originally announced September 2015.
-
Surface Energetics of Alkaline-Earth Metal Oxides: Trends in Stability and Adsorption of Small Molecules
Authors:
Michal Bajdich,
Jens K. Nørskov,
Aleksandra Vojvodic
Abstract:
We present a systematic theoretical investigation of the surface properties, stability and reactivity, of rock-salt type alkaline-earth metal oxides including MgO, CaO, SrO, and BaO. The accuracy of commonly used exchange-correlation density functionals (LDA, PBE, RPBE, PBEsol, BEEF-vdW and hybrid HSE) and random-phase approximation (RPA) is evaluated and compared to existing experimental values.…
▽ More
We present a systematic theoretical investigation of the surface properties, stability and reactivity, of rock-salt type alkaline-earth metal oxides including MgO, CaO, SrO, and BaO. The accuracy of commonly used exchange-correlation density functionals (LDA, PBE, RPBE, PBEsol, BEEF-vdW and hybrid HSE) and random-phase approximation (RPA) is evaluated and compared to existing experimental values. Calculated surface energies of the four most stable surface facets under vacuum conditions: the (100) surface, the metal and oxygen terminated octopolar (111), and the (110) surfaces exhibit a monotonic increase in stability from MgO to BaO. On the MgO(100) surface, adsorption of CO, NO, CH4 is characterized by physisorption while H2O chemisorbs, which is in agreement with experimental findings. We further use the on-top metal adsorption of CO and NO molecules to map out the surface energetics of each alkaline-earth metal oxide surface. The considered functionals all qualitatively predict similar adsorption energy trends. The ordering between the adsorption energies on different surface facets can be attributed to differences in the local geometrical surface structure and the electronic structure of the metal constituent of the alkaline-earth metal oxide. The striking observation that CO adsorption strength is weaker than NO adsorption on the (100) terraces as the period of the alkaline-earth metal in the oxide increases, is analyzed in detail in terms of charge redistribution within the σ and π channels of adsorbates. Finally, we also present oxygen adsorption and oxygen vacancy formation energies in these oxide systems.
△ Less
Submitted 21 March, 2015; v1 submitted 19 December, 2014;
originally announced December 2014.
-
Insights into carbon nanotube nucleation: Cap formation governed by catalyst interfacial step flow
Authors:
Rahul Rao,
Renu Sharma,
Frank Abild-Pedersen,
Jens K. Nørskov,
Avetik R. Harutyunyan
Abstract:
In order to accommodate an increasing demand for carbon nanotubes (CNTs) with desirable characteristics one has to understand the origin of helicity of their structures. Here, through in situ microscopy we demonstrate that the nucleation of a carbon nanotube is initiated by the formation of the carbon cap. Nucleation begins with the formation of a graphene embryo that is bound between opposite ste…
▽ More
In order to accommodate an increasing demand for carbon nanotubes (CNTs) with desirable characteristics one has to understand the origin of helicity of their structures. Here, through in situ microscopy we demonstrate that the nucleation of a carbon nanotube is initiated by the formation of the carbon cap. Nucleation begins with the formation of a graphene embryo that is bound between opposite step-edges on the nickel catalyst surface. The embryo grows larger as the step-edges migrate along the surface, leading to the formation of a curved carbon cap when the steps flow across the edges of adjacent facets. Further motion of the steps away from the catalyst tip with attached rims of the carbon cap generates the wall of the nanotube. Density Functional Theory calculations bring further insight into the process, showing that step flow occurs by surface self diffusion of the nickel atoms via a step-edge attachment-detachment mechanism. Since the fact that cap forms first in the sequence of stages involved in nanotube growth, we suggest that it originates the helicity of the nanotube. Therefore, the angular distribution of catalyst facets could be exploited as a new parameter for controlling the curvature of the cap and, presumably, the helicity of the nanotube.
△ Less
Submitted 8 August, 2014; v1 submitted 7 May, 2014;
originally announced May 2014.
-
Ab initio van der Waals interactions in simulations of water alter structure from mainly tetrahedral to high-density-like
Authors:
Andreas Møgelhøj,
André Kelkkanen,
K. Thor Wikfeldt,
Jakob Schiøtz,
Jens Jørgen Mortensen,
Lars G. M. Pettersson,
Bengt I. Lundqvist,
Karsten W. Jacobsen,
Anders Nilsson,
Jens K. Nørskov
Abstract:
The structure of liquid water at ambient conditions is studied in ab initio molecular dynamics simulations using van der Waals (vdW) density-functional theory, i.e. using the new exchange-correlation functionals optPBE-vdW and vdW-DF2. Inclusion of the more isotropic vdW interactions counteracts highly directional hydrogen-bonds, which are enhanced by standard functionals. This brings about a soft…
▽ More
The structure of liquid water at ambient conditions is studied in ab initio molecular dynamics simulations using van der Waals (vdW) density-functional theory, i.e. using the new exchange-correlation functionals optPBE-vdW and vdW-DF2. Inclusion of the more isotropic vdW interactions counteracts highly directional hydrogen-bonds, which are enhanced by standard functionals. This brings about a softening of the microscopic structure of water, as seen from the broadening of angular distribution functions and, in particular, from the much lower and broader first peak in the oxygen-oxygen pair-correlation function (PCF), indicating loss of structure in the outer solvation shells. In combination with softer non-local correlation terms, as in the new parameterization of vdW-DF, inclusion of vdW interactions is shown to shift the balance of resulting structures from open tetrahedral to more close-packed. The resulting O-O PCF shows some resemblance with experiment for high-density water (A. K. Soper and M. A. Ricci, Phys. Rev. Lett., 84:2881, 2000), but not directly with experiment for ambient water. However, an O-O PCF consisting of a linear combination of 70% from vdW-DF2 and 30% from experiment on low-density liquid water reproduces near-quantitatively the experimental O-O PCF for ambient water, indicating consistency with a two-liquid model with fluctuations between high- and low-density regions.
△ Less
Submitted 3 May, 2011; v1 submitted 28 January, 2011;
originally announced January 2011.
-
Trends in Metal Oxide Stability for Nanorods, Nanotubes, and Surfaces
Authors:
D. J. Mowbray,
J. I. Martinez,
F. Calle-Vallejo,
J. Rossmeisl,
K. S. Thygesen,
K. W. Jacobsen,
J. K. Norskov
Abstract:
The formation energies of nanostructures play an important role in determining their properties, including the catalytic activity. For the case of 15 different rutile and 8 different perovskite metal oxides, we find that the density functional theory (DFT) calculated formation energies of (2,2) nanorods, (3,3) nanotubes, and the (110) and (100) surfaces may be described semi-quantitatively by th…
▽ More
The formation energies of nanostructures play an important role in determining their properties, including the catalytic activity. For the case of 15 different rutile and 8 different perovskite metal oxides, we find that the density functional theory (DFT) calculated formation energies of (2,2) nanorods, (3,3) nanotubes, and the (110) and (100) surfaces may be described semi-quantitatively by the fraction of metal--oxygen bonds broken and the bonding band centers in the bulk metal oxide.
△ Less
Submitted 25 February, 2010;
originally announced February 2010.
-
Bayesian Error Estimation in Density Functional Theory
Authors:
J. J. Mortensen,
K. Kaasbjerg,
S. L. Frederiksen,
J. K. Norskov,
J. P. Sethna,
K. W. Jacobsen
Abstract:
We present a practical scheme for performing error estimates for Density Functional Theory calculations. The approach which is based on ideas from Bayesian statistics involves creating an ensemble of exchange-correlation functionals by comparing with an experimental database of binding energies for molecules and solids. Fluctuations within the ensemble can then be used to estimate errors relativ…
▽ More
We present a practical scheme for performing error estimates for Density Functional Theory calculations. The approach which is based on ideas from Bayesian statistics involves creating an ensemble of exchange-correlation functionals by comparing with an experimental database of binding energies for molecules and solids. Fluctuations within the ensemble can then be used to estimate errors relative to experiment on calculated quantities like binding energies, bond lengths, and vibrational frequencies. It is demonstrated that the error bars on energy differences may vary by orders of magnitude for different systems in good agreement with existing experience.
△ Less
Submitted 9 June, 2005;
originally announced June 2005.
-
Pareto-optimal alloys
Authors:
Thomas Bligaard,
Gisli H. Johannesson,
Andrei V. Ruban,
Hans L. Skriver,
Karsten W. Jacobsen,
Jens K. Norskov
Abstract:
Large databases that can be used in the search for new materials with specific properties remain an elusive goal in materials science. The search problem is complicated by the fact that the optimal material for a given application is usually a compromise between a number of materials properties and the price. In this letter we present a database consisting of the lattice parameters, bulk moduli,…
▽ More
Large databases that can be used in the search for new materials with specific properties remain an elusive goal in materials science. The search problem is complicated by the fact that the optimal material for a given application is usually a compromise between a number of materials properties and the price. In this letter we present a database consisting of the lattice parameters, bulk moduli, and heats of formation for over 64,000 ordered metallic alloys, which has been established by direct first-principles density-functional-theory calculations. Furthermore, we use a concept from economic theory, the Pareto-optimal set, to determine optimal alloy solutions for the compromise between low compressibility, high stability and price.
△ Less
Submitted 29 October, 2003;
originally announced October 2003.
-
Incomplete melting of the Si(100) surface from molecular-dynamics simulations using the Effective-Medium Tight-Binding model
Authors:
K. Stokbro K. W. Jacobsen,
J. K. Nørskov,
D. M. Deaven,
C. Z. Wang,
K. M. Ho
Abstract:
We present molecular-dynamics simulations of the Si(100) surface in the temperature range 1100-1750K. To describe the total energy and forces we use the Effective-Medium Tight-Binding model. The defect-free surface is found to melt at the bulk melting point, which we determine to be 1650 K, but for a surface with dimer vacancies we find a pre-melting of the first two layers 100 K below the melti…
▽ More
We present molecular-dynamics simulations of the Si(100) surface in the temperature range 1100-1750K. To describe the total energy and forces we use the Effective-Medium Tight-Binding model. The defect-free surface is found to melt at the bulk melting point, which we determine to be 1650 K, but for a surface with dimer vacancies we find a pre-melting of the first two layers 100 K below the melting point. We show that these findings can rationalize recent experimental studies of the high temperature Si(100) surface.
△ Less
Submitted 2 April, 1996;
originally announced April 1996.
-
A direct pathway for sticking/desorption of H$_2$ on Si(100)
Authors:
P. Kratzer,
B. Hammer,
J. K. Norskov
Abstract:
The energetics of H$_2$ interacting with the Si(100) surface is studied by means of {\em ab initio} total energy calculations within the framework of density functional theory. We find a direct desorption pathway from the mono-hydride phase which is compatible with experimental activation energies and demonstrate the importance of substrate relaxation for this process. Both the transition state…
▽ More
The energetics of H$_2$ interacting with the Si(100) surface is studied by means of {\em ab initio} total energy calculations within the framework of density functional theory. We find a direct desorption pathway from the mono-hydride phase which is compatible with experimental activation energies and demonstrate the importance of substrate relaxation for this process. Both the transition state configuration and barrier height depend crucially on the degree of buckling of the Si dimers on the Si(100) surface. The adsorption barrier height on the clean surface is governed by the buckling via its influence on the surface electronic structure. We discuss the consequences of this coupling for adsorption experiments and the relation between adsorption and desorption.
△ Less
Submitted 29 March, 1995; v1 submitted 27 March, 1995;
originally announced March 1995.
-
An Effective-Medium Tight-Binding Model for Silicon
Authors:
K. Stokbro,
N. Chetty,
K. W. Jacobsen,
J. K. Nørskov
Abstract:
A new method for calculating the total energy of Si systems is presented. The method is based on the effective-medium theory concept of a reference system. Instead of calculating the energy of an atom in the system of interest a reference system is introduced where the local surroundings are similar. The energy of the reference system can be calculated selfconsistently once and for all while the…
▽ More
A new method for calculating the total energy of Si systems is presented. The method is based on the effective-medium theory concept of a reference system. Instead of calculating the energy of an atom in the system of interest a reference system is introduced where the local surroundings are similar. The energy of the reference system can be calculated selfconsistently once and for all while the energy difference to the reference system can be obtained approximately. We propose to calculate it using the tight-binding LMTO scheme with the Atomic-Sphere Approximation(ASA) for the potential, and by using the ASA with charge-conserving spheres we are able to treat open system without introducing empty spheres. All steps in the calculational method is {\em ab initio} in the sense that all quantities entering are calculated from first principles without any fitting to experiment. A complete and detailed description of the method is given together with test calculations of the energies of phonons, elastic constants, different structures, surfaces and surface reconstructions. We compare the results to calculations using an empirical tight-binding scheme.
△ Less
Submitted 9 May, 1994;
originally announced May 1994.
-
Construction of transferable spherically-averaged electron potentials
Authors:
K. Stokbro,
N. Chetty,
K. W. Jacobsen,
J. K. Nørskov
Abstract:
A new scheme for constructing approximate effective electron potentials within density-functional theory is proposed. The scheme consists of calculating the effective potential for a series of reference systems, and then using these potentials to construct the potential of a general system. To make contact to the reference system the neutral-sphere radius of each atom is used. The scheme can sim…
▽ More
A new scheme for constructing approximate effective electron potentials within density-functional theory is proposed. The scheme consists of calculating the effective potential for a series of reference systems, and then using these potentials to construct the potential of a general system. To make contact to the reference system the neutral-sphere radius of each atom is used. The scheme can simplify calculations with partial wave methods in the atomic-sphere or muffin-tin approximation, since potential parameters can be precalculated and then for a general system obtained through simple interpolation formulas. We have applied the scheme to construct electron potentials of phonons, surfaces, and different crystal structures of silicon and aluminum atoms, and found excellent agreement with the self-consistent effective potential. By using an approximate total electron density obtained from a superposition of atom-based densities, the energy zero of the corresponding effective potential can be found and the energy shifts in the mean potential between inequivalent atoms can therefore be directly estimated. This approach is shown to work well for surfaces and phonons of silicon.
△ Less
Submitted 9 May, 1994;
originally announced May 1994.
-
Melting a copper cluster: Critical droplet theory
Authors:
Ole H. Nielsen,
James P. Sethna,
Per Stoltze,
Karsten W. Jacobsen,
Jens K. Norskov
Abstract:
We simulate the melting of a 71 A diameter cluster of Cu. At low temperatures the crystal exhibits facets. With increasing temperatures the open facets pre-melt, the melted regions coalesce into a liquid envelope containing a crystalline nucleus, and the nucleus finally goes unstable to the supercooled liquid. Using critical droplet theory and experimental data for Cu, we explain the thermodynam…
▽ More
We simulate the melting of a 71 A diameter cluster of Cu. At low temperatures the crystal exhibits facets. With increasing temperatures the open facets pre-melt, the melted regions coalesce into a liquid envelope containing a crystalline nucleus, and the nucleus finally goes unstable to the supercooled liquid. Using critical droplet theory and experimental data for Cu, we explain the thermodynamics of the coexistence region. The width of the transition scales as (Number of particles) to the power (-1/4).
△ Less
Submitted 20 July, 1992;
originally announced July 1992.