Michal Bajdich

Identification of highly active Fe sites in (Ni,Fe)OOH for electrocatalytic water splitting

Highly active catalysts for the oxygen evolution reaction (OER) are required for the development of photoelectrochemical devices that generate hydrogen efficiently from water using solar energy. Here, we identify the origin of a 500-fold OER activity enhancement that can be achieved with mixed (Ni,Fe)oxyhydroxides (Ni(1-x)Fe(x)OOH) over their pure Ni and Fe parent compounds, resulting in one of the most active currently known OER catalysts in alkaline electrolyte. Operando X-ray absorption spectroscopy (XAS) using high energy resolution fluorescence detection (HERFD) reveals that Fe(3+) in Ni(1-x)Fe(x)OOH occupies octahedral sites with unusually short Fe-O bond distances, induced by edge-sharing with surrounding [NiO6] octahedra. Using computational methods, we establish that this structural motif results in near optimal adsorption energies of OER intermediates and low overpotentials at Fe sites. By contrast, Ni sites in Ni(1-x)Fe(x)OOH are not active sites for the oxidation of water.

Homogeneously dispersed, multimetal oxygen-evolving catalysts

Modulating metal oxides The more difficult step in fuel cells and water electrolysis is the oxygen evolution reaction. The search for earth-abundant materials to replace noble metals for this reaction often turns to oxides of three-dimensional metals such as iron. Zhang et al. show that the applied voltages needed to drive this reaction are reduced for iron-cobalt oxides by the addition of tungsten. The addition of tungsten favorably modulates the electronic structure of the oxyhydroxide. A key development is to keep the metals well mixed and avoid the formation of separate phases. Science, this issue p. 333 The addition of tungsten to iron cobalt oxides lowers the overpotential required for the evolution of oxygen from water. Earth-abundant first-row (3d) transition metal–based catalysts have been developed for the oxygen-evolution reaction (OER); however, they operate at overpotentials substantially above thermodynamic requirements. Density functional theory suggested that non-3d high-valency metals such as tungsten can modulate 3d metal oxides, providing near-optimal adsorption energies for OER intermediates. We developed a room-temperature synthesis to produce gelled oxyhydroxides materials with an atomically homogeneous metal distribution. These gelled FeCoW oxyhydroxides exhibit the lowest overpotential (191 millivolts) reported at 10 milliamperes per square centimeter in alkaline electrolyte. The catalyst shows no evidence of degradation after more than 500 hours of operation. X-ray absorption and computational studies reveal a synergistic interplay between tungsten, iron, and cobalt in producing a favorable local coordination environment and electronic structure that enhance the energetics for OER.

Theoretical investigation of the activity of cobalt oxides for the electrochemical oxidation of water

The presence of layered cobalt oxides has been identified experimentally in Co-based anodes under oxygen-evolving conditions. In this work, we report the results of theoretical investigations of the relative stability of layered and spinel bulk phases of Co oxides, as well as the stability of selected surfaces as a function of applied potential and pH. We then study the oxygen evolution reaction (OER) on these surfaces and obtain activity trends at experimentally relevant electro-chemical conditions. Our calculated volume Pourbaix diagram shows that β-CoOOH is the active phase where the OER occurs in alkaline media. We calculate relative surface stabilities and adsorbate coverages of the most stable low-index surfaces of β-CoOOH: (0001), (0112), and (1014). We find that at low applied potentials, the (1014) surface is the most stable, while the (0112) surface is the more stable at higher potentials. Next, we compare the theoretical overpotentials for all three surfaces and find that the (1014) surface is the most active one as characterized by an overpotential of η = 0.48 V. The high activity of the (1014) surface can be attributed to the observation that the resting state of Co in the active site is Co(3+) during the OER, whereas Co is in the Co(4+) state in the less active surfaces. Lastly, we demonstrate that the overpotential of the (1014) surface can be lowered further by surface substitution of Co by Ni. This finding could explain the experimentally observed enhancement in the OER activity of Ni(y)Co(1-y)O(x) thin films with increasing Ni content. All energetics in this work were obtained from density functional theory using the Hubbard-U correction.

QWalk: A quantum Monte Carlo program for electronic structure

We describe QWalk, a new computational package capable of performing quantum Monte Carlo electronic structure calculations for molecules and solids with many electrons. We describe the structure of the program and its implementation of quantum Monte Carlo methods. It is open-source, licensed under the GPL, and available at the web site http://www.qwalk.org.

Pfaffian pairing wave functions in electronic-structure quantum Monte Carlo simulations.

We investigate the accuracy of trial wave functions for quantum Monte Carlo based on Pfaffian functional form with singlet and triplet pairing. Using a set of first row atoms and molecules we find that these wave functions provide very consistent and systematic behavior in recovering the correlation energies on the level of 95%. In order to get beyond this limit we explore the possibilities of multi-Pfaffian pairing wave functions. We show that a small number of Pfaffians recovers another large fraction of the missing correlation energy comparable to the larger-scale configuration interaction wave functions. We also find that Pfaffians lead to substantial improvements in fermion nodes when compared to Hartree-Fock wave functions.

On the chemical state of Co oxide electrocatalysts during alkaline water splitting

Resonant inelastic X-ray scattering and high-resolution X-ray absorption spectroscopy were used to identify the chemical state of a Co electrocatalyst in situ during the oxygen evolution reaction. After anodic electrodeposition onto Au(111) from a Co(2+)-containing electrolyte, the chemical environment of Co can be identified to be almost identical to CoOOH. With increasing potentials, a subtle increase of the Co oxidation state is observed, indicating a non-stoichiometric composition of the working OER catalyst containing a small fraction of Co(4+) sites. In order to confirm this interpretation, we used density functional theory with a Hubbard-U correction approach to compute X-ray absorption spectra of model compounds, which agree well with the experimental spectra. In situ monitoring of catalyst local structure and bonding is essential in the development of structure-activity relationships that can guide the discovery of efficient and earth abundant water splitting catalysts.

Pfaffian pairing and backflow wavefunctions for electronic structure quantum Monte Carlo methods

We investigate pfaffian trial wavefunctions with singlet and triplet pair orbitals by quantum Monte Carlo methods. We present mathematical identities and the key algebraic properties necessary for efficient evaluation of pfaffians. Following upon our previous study [Bajdich et al., Phys. Rev. Lett. 96, 130201 (2006)], we explore the possibilities of expanding the wavefunction in linear combinations of pfaffians. We observe that molecular systems require much larger expansions than atomic systems and linear combinations of a few pfaffians lead to rather small gains in correlation energy. We also test the wavefunction based on fully antisymmetrized product of independent pair orbitals. Despite its seemingly large variational potential, we do not observe additional gains in correlation energy. We find that pfaffians lead to substantial improvements in fermion nodes when compared to Hartree-Fock wavefunctions and exhibit the minimal number of two nodal domains in agreement with recent results on fermion nodes topology. We analyze the nodal structure differences of Hartree-Fock, pfaffian, and essentially exact large-scale configuration interaction wavefunctions. Finally, we combine the recently proposed form of backflow correlations [Drummond et al., J. Phys. Chem. 124, 22401 (2006); Rios et al., Phys. Rev. E. 74, 066701 (2006)] with both determinantal and pfaffian based wavefunctions.

Interface Controlled Oxidation States in Layered Cobalt Oxide Nanoislands on Gold

Layered cobalt oxides have been shown to be highly active catalysts for the oxygen evolution reaction (OER; half of the catalytic water splitting reaction), particularly when promoted with gold. However, the surface chemistry of cobalt oxides and in particular the nature of the synergistic effect of gold contact are only understood on a rudimentary level, which at present prevents further exploration. We have synthesized a model system of flat, layered cobalt oxide nanoislands supported on a single crystal gold (111) substrate. By using a combination of atom-resolved scanning tunneling microscopy, X-ray photoelectron and absorption spectroscopies and density functional theory calculations, we provide a detailed analysis of the relationship between the atomic-scale structure of the nanoislands, Co oxidation states and substrate induced charge transfer effects in response to the synthesis oxygen pressure. We reveal that conversion from Co(2+) to Co(3+) can occur by a facile incorporation of oxygen at the interface between the nanoisland and gold, changing the islands from a Co-O bilayer to an O-Co-O trilayer. The O-Co-O trilayer islands have the structure of a single layer of β-CoOOH, proposed to be the active phase for the OER, making this system a valuable model in understanding of the active sites for OER. The Co oxides adopt related island morphologies without significant structural reorganization, and our results directly demonstrate that nanosized Co oxide islands have a much higher structural flexibility than could be predicted from bulk properties. Furthermore, it is clear that the gold/nanoparticle interface has a profound effect on the structure of the nanoislands, suggesting a possible promotion mechanism.

Approximate and exact nodes of fermionic wavefunctions : Coordinate transformations and topologies

A study of fermion nodes for spin-polarized states of a few-electron ions and molecules with s,p,d one-particle orbitals is presented. We find exact nodes for some cases of two-electron atomic and molecular states and also the first exact node for the three-electron atomic system in {sup 4}S(p{sup 3}) state using appropriate coordinate maps and wave function symmetries. We analyze the cases of nodes for larger number of electrons in the Hartree-Fock approximation and for some cases we find transformations for projecting the high-dimensional node manifolds into three-dimensional space. The node topologies and other properties are studied using these projections. We also propose a general coordinate transformation as an extension of Feynman-Cohen backflow coordinates to both simplify the nodal description and as a new variational freedom for quantum Monte Carlo trial wave functions.

Subnanometer-sized Pt/Sn alloy cluster catalysts for the dehydrogenation of linear alkanes.

The reaction pathways for the dehydrogenation of ethane, propane, and butane, over Pt are analyzed using density functional theory (DFT). Pt nanoparticles are represented by a tetrahedral Pt4 cluster. The objectives of this work were to establish which step is rate limiting and which one controls the selectivity for forming alkenes as opposed to causing further dehydrogenation of adsorbed alkenes to produce precursors responsible for catalyst deactivation due to coking. Further objectives of this work are to identify the role of adsorbed hydrogen, derived from H2 fed together with the alkane, on the reaction pathway, and the role of replacing one of the four Pt atoms by a Sn atom. A comparison of Gibbs free energies shows that in all cases the rate-determining step is cleavage of a C-H bond upon alkane adsorption. The selectivity to alkene formation versus precursors to coking is dictated by the relative magnitudes of the activation energies for alkene desorption and dehydrogenation of the adsorbed alkene. The presence of an adsorbed H atom on the cluster facilitates alkene desorption relative to dehydrogenation of the adsorbed alkene. Substitution of a Sn atom in the cluster to produce a Pt3Sn cluster leads to a downward shift of the potential energy surface for the reaction and causes an increase of the activity of the catalyst as suggested by recent experiments due to the lower net activation barrier for the rate limiting step. However, the introduction of Sn does not alter the relative activation barriers for gas-phase alkene formation versus loss of hydrogen from the adsorbed alkene, the process leading to the formation of coke precursors.