Fe-containing transition-metal (oxy)hydroxides are highly active oxygen-evolution reaction (OER) electrocatalysts in alkaline media and ubiquitously form across many materials systems. The complexity and dynamics of the Fe sites within the (oxy)hydroxide have slowed understanding of how and where the Fe-based active sites form-information critical for designing catalysts and electrolytes with higher activity and stability. We show that where/how Fe species in the electrolyte incorporate into host Ni or Co (oxy)hydroxides depends on the electrochemical history and structural properties of the host material. Substantially less Fe is incorporated from Fe-spiked electrolyte into Ni (oxy)hydroxide at anodic potentials, past the nominally Ni2+/3+ redox wave, compared to during potential cycling. The Fe adsorbed under constant anodic potentials leads to impressively high per-Fe OER turn-over frequency (TOFFe) of ~40 s-1 at 350 mV overpotential which we attribute to under-coordinated "surface" Fe. By systematically controlling the concentration of surface Fe, we find TOFFe increases linearly with the Fe concentration. This suggests a changing OER mechanism with increased Fe concentration, consistent with a mechanism involving cooperative Fe sites in FeOx clusters.
Many key features in photovoltaic perovskites occur in relatively long time scales and involve mixed compositions. This requires realistic but also numerically simple models. In this work we present a transferable classical force field to describe the mixed hybrid perovskite MAxFA1-xPb(BryI1-y)3 for variable composition (∀x, y ∈ [0, 1]). The model includes Lennard-Jones and Buckingham potentials to describe the interactions between the atoms of the inorganic lattice and the organic molecule, and the AMBER model to describe intramolecular atomic interactions. Most of the parameters of the force field have been obtained by means of a genetic algorithm previously developed to parametrize the CsPb(BrxI1-x)3 perovskite (Balestra et al. J. Mater. Chem. A. 2020, DOI: 10.1039/d0ta03200j). The algorithm finds the best parameter set that simultaneously fits the DFT energies obtained for several crystalline structures with moderate degrees of distortion with respect to the equilibrium configuration. The resulting model reproduces correctly the XRD patterns, the expansion of the lattice upon I/Br substitution, and the thermal expansion coefficients. We use the model to run classical molecular dynamics simulations with up to 8600 atoms and simulation times of up to 40 ns. From the simulations we have extracted the ion diffusion coefficient of the pure and mixed perovskites, presenting for the first time these values obtained by a fully dynamical method using a transferable model fitted to first-principles calculations. The values here reported can be considered as the theoretical upper limit, that is, without grain boundaries or other defects, for ion migration dynamics induced by halide vacancies in photovoltaic perovskite devices under operational conditions.
Metal-organic frameworks (MOFs) are commended as photocatalysts for H2 evolution and CO2 reduction as they combine light-harvesting and catalytic functions with excellent reactant adsorption capabilities. For dynamic processes in liquid phase, the accessibility of active sites becomes a critical parameter as reactant diffusion is limited by the inherently small micropores. Our strategy is to introduce additional mesopores by selectively removing one ligand in mixed-ligand MOFs via thermolysis. Here we report photoactive MOFs of the MIL-125-Ti family with two distinct mesopore architectures resembling either large cavities or branching fractures. The ligand removal is highly selective and follows a 2-step process tunable by temperature and time. The introduction of mesopores and the associated formation of new active sites have improved the HER rates of the MOFs by up to 500%. We envision that this strategy will allow the purposeful engineering of hierarchical MOFs and advance their applicability in environmental and energy technologies.
Hematite is a possible photoanode for photoelectrochemical cells (PECs) that is widely studied using density functional theory (DFT). In this paper we perform more accurate calculations of the absorption spectrum of hematite using the one-shot Green's function (G0W0) and Bethe-Salpeter equation (BSE) methods, which take excited states into account and compare the spectrum to experimental data. We found a match between our calculations and the observed absorption spectra in peak locations. Furthermore, there is anisotropy of the absorption spectra that is concurrent with the crystal structure. We also calculated the absorption spectrum of hematite intermediates during catalysis of the oxygen evolution reaction to better understand which intermediate is dominant during the reaction and the contribution of excited states to catalysis. The *O intermediate was found to be the most optically and chemically dominant species during catalysis.
Iron cations are essential for the high activity of nickel and cobalt-based (oxy)hydroxides for the oxygen evolution reaction, but the role of iron in the catalytic mechanism remains under active investigation. Operando X-ray absorption spectroscopy and density functional theory calculations are used to demonstrate partial Fe oxidation and a shortening of the Fe-O bond length during oxygen evolution on Co(Fe)Ox Hy . Cobalt oxidation during oxygen evolution is only observed in the absence of iron. These results demonstrate a different mechanism for water oxidation in the presence and absence of iron and support the hypothesis that oxidized iron species are involved in water-oxidation catalysis on Co(Fe)Ox Hy .
In recent years, hematite's potential as a photoanode material for solar hydrogen production has ignited a renewed interest in its physical and interfacial properties, which continues to be an active field of research. Research on hematite photoanodes provides new insights on the correlations between electronic structure, transport properties, excited state dynamics, and charge transfer phenomena, and expands our knowledge on solar cell materials into correlated electron systems. This research news article presents a snapshot of selected theoretical and experimental developments linking the electronic structure to the photoelectrochemical performance, with particular focus on optoelectronic properties and charge carrier dynamics.
Due to the high oxygen evolution reaction (OER) activity, stability, and abundance of NiOx materials, they are found to be promising catalysts, competitive with expensive metal oxides such as IrO2 and RuO2. From a theoretical point of view, studies reported in the literature so far are mostly based on density functional theory using periodic slab models for the bulk and surface of ß-NiOOH, one of the active NiOx phases. However, cluster models are a valid method to investigate many aspects about structure, charge carrier transport properties, and OER activity of ß-NiOOH. Hence, here we present new cluster models for the surface of ß-NiOOH, where the oxygen atoms are bonded to Mg effective core potentials (ECPs) mimicking neighboring atom cores. This cluster embedding procedure is superior to saturating the cluster with hydrogen atoms, and to using other atomic ECPs for ß-NiOOH. We find that layered materials such as ß-NiOOH are more vulnerable to geometrical rupture and therefore a cluster approach requires additional care in choosing the embedding approach. We evaluated the models by using them to calculate the energy required for water adsorption and deprotonation, which are essential ingredients for OER. Specifically, our results agree with previous slab models that the first deprotonation reaction step requires a large amount of energy. In addition, we find that water and hydroxyl groups have high adsorption energy and therefore the first deprotonation step is limiting the reaction efficiency.
Discovering better catalysts for water splitting is the holy grail of the renewable energy field. One of the most successful water oxidation catalysts is nickel oxyhydroxide (NiOOH), which is chemically active only as a result of doping with Fe. In order to shed light on how Fe improves efficiency, we perform Density Functional Theory +U (DFT+U) calculations of water oxidation reaction intermediates of Fe substitutional doped NiOOH. The results are analyzed while considering the presence of vacancies that we use as probes to test the effect of adding charge to the surface. We find that the smaller electronegativity of the Fe dopant relative to Ni allows the dopant to have several possible oxidation states with less energy penalty. As a result, the presence of vacancies which alters local oxidation states does not affect the low overpotential of Fe-doped NiOOH. We conclude that the secret to the success of doping NiOOH with iron is the ability of iron to easily change oxidation states, which is critical during the chemical reaction of water oxidation.
Improving water splitting efficiency has been the holy grail of hydrogen fuel production. Major efforts have been invested in an attempt to enhance efficiency of a common water oxidation catalyst, α-Fe2O3, through doping and alloying. Recent experiments show that higher efficiency is achieved when niobium (Nb) is added beyond the solubility limit to generate a mixture of two phases: Nb-doped and Nb-alloyed α-Fe2O3. In order to understand why adding high concentrations of Nb is beneficial, we provide a thorough first principles study of the bulk and the surface of pure, Nb-doped, and Nb-alloyed α-Fe2O3 with several surface facets and terminations. We find that the addition of Nb changes the band edge and Fermi level positions. Therefore, we propose a mechanism by which having different Nb doping levels within and above the solubility limit has an advantage: electrons and holes could separate better between doped and alloyed regions that have different band edge positions or between regions with different doping concentrations. Furthermore, the holes' driving force to oxidize water can be increased by placing on the surface the undoped or alloyed phases, since they have a lower valence band maximum. We suggest that obtaining two material phases or gradual doping can be used as a design strategy for next generation catalysts.
The computational design of solid catalysts has become a very "hot" field during the last decades, especially with the recent increase in computational tool performance. However, theoretical techniques are still very time demanding because they require the consideration of many adsorption configurations of the reaction intermediates on the surface. Herein, we propose to use the metal-oxygen (M-O) bond ionicity as a descriptor for the photocatalytic activity of one of the best catalysts for the oxygen evolution reaction (OER). Ionicity is a bulk property and thus carries the advantage of being easily obtainable from a simple Bader charge analysis by using density functional theory (DFT). We will show that this criterion can be used successfully to design efficient dopants for NiOOH material. This catalyst is known to exhibit interesting photoelectrochemical properties for OER if it is doped with specific transition metals. Finally, we demonstrate that other electronic properties that relate to bulk calculation, such as oxidation states and density of states, are not alone sufficient to explain the photocatalytic activity of the material. Thus, M-O bond ionicity attracts significant interest compared with other bulk observables obtained by using DFT computations.
Metal/oxide interfaces have long been studied for their fundamental importance in material microstructure as well as their broad applicability in electronic devices. However, the challenge involved in characterizing the relation between structure and electron transport of a large number of metal/oxide combinations inhibits the search for interfaces with improved functionality. Therefore, we develop a novel high-throughput screening approach that combines computational and theoretical techniques. We use a Density Functional Theory + U (DFT+U) quantum mechanical formalism to produce effective Schrödinger equations, which are solved by wave packet propagation to simulate charge transport across the metal/oxide interface. We demonstrate this method on α-Fe2O3/Mt interfaces, for Mt = Ag, Al, Au, Ir, Pd, or Pt metals. We use this novel method to screen for binary alloys of these metals at the α-Fe2O3/Mt interface and perform a successful validation test of the methodology. Finally, we correlate the interface potential energy and the charge transport permeability through the interface. Counterintuitively, among the interfaces studied, we find that higher mismatch interfaces have better charge transport permeability. We anticipate that this method will be useful as a computationally tractable strategy to perform high-throughput screening for new metal/oxide interfaces.
The challenge of improving the efficiency of photo-electrochemical devices is often addressed through doping. However, this strategy could harm performance. Specifically, as demonstrated in a recent experiment, doping one of the most widely used materials for water splitting, iron (III) oxide (Fe2O3), with niobium (Nb) can still result in limited efficiency. In order to better understand the hazardous effect of doping, we use Density Functional Theory (DFT)+U for the case of Nb-doped Fe2O3. We find a direct correlation between the charge of the dopant, the charge on surface of the Fe2O3 material, and the overpotential required for water oxidation reaction. We believe that this work contributes to advancing our understanding of how to select effective dopants for materials.
Electrochemistry , Ferric Compounds/chemistry , Models, Theoretical , Niobium/chemistry , Photochemistry
Metal/metal-oxide interfaces appear in a wide variety of disciplines including electronics, corrosion, electrochemistry, and catalysis. Specifically, covering a metal-oxide with a metal is often thought to enhance solar energy absorption and to improve photocatalytic activity. For example, the platinum/hematite (Pt/α-Fe2O3) interface has demonstrated improved functionality. In order to advance our understanding of how metal coverage over an oxide helps performance, we characterize the geometry and electronic structure of the Pt/α-Fe2O3 interface. We investigate the interface using density functional theory +U, and find a stable crystallographic orientation relationship: Fe2O3(0001)[11[combining macron]00]â¥Pt(111)[101[combining macron]] that agrees with experiment. Furthermore, there are significant changes in the electronic structure of α-Fe2O3 as a result of Pt coverage. We therefore suggest the concept of "judging" the electronic properties of an oxide only with its cover. Specifically, covering Fe2O3 with Pt reduces carrier effective mass and creates a continuum of states in the band gap. The former could be beneficial for catalytic activity, while the latter may cause surface recombination. In order to circumvent this problem, we suggest putting metal coverage behind the oxide and far from the electrolyte in a photoelectrochemical device in order to quickly collect electron carriers and avoid recombination with vulnerable holes accumulating as a result of catalysis at the surface.
We use two different ab initio quantum mechanics methods, complete active space self-consistent field theory applied to electrostatically embedded clusters and periodic many-body G0W0 calculations, to reanalyze the states formed in nickel(II) oxide upon electron addition and ionization. In agreement with interpretations of earlier measurements, we find that the valence and conduction band edges consist of oxygen and nickel states, respectively. However, contrary to conventional wisdom, we find that the oxygen states of the valence band edge are localized whereas the nickel states at the conduction band edge are delocalized. We argue that these characteristics may lead to low electron-hole recombination and relatively efficient electron transport, which, coupled with band gap engineering, could produce higher solar energy conversion efficiency compared to that of other transition-metal oxides. Both methods find a photoemission/inverse-photoemission gap of 3.6-3.9 eV, in good agreement with the experimental range, lending credence to our analysis of the electronic structure of NiO.
Long-term sustainable solar energy conversion relies on identifying economical and versatile semiconductor materials with appropriate band structures for photovoltaic and photocatalytic applications (e.g., band gaps of â¼ 1.5-2.0 eV). Nickel oxide (NiO) is an inexpensive yet highly promising candidate. Its charge-transfer character may lead to longer carrier lifetimes needed for higher efficiencies, and its conduction band edge is suitable for driving hydrogen evolution via water-splitting. However, NiO's large band gap (â¼ 4 eV) severely limits its use in practical applications. Our first-principles quantum mechanics calculations show band gaps dramatically decrease to â¼ 2.0 eV when NiO is alloyed with Li2O. We show that Lix Ni1-x O alloys (with x=0.125 and 0.25) are p-type semiconductors, contain states with no impurity levels in the gap and maintain NiO's desirable charge-transfer character. Lastly, we show that the alloys have potential for photoelectrochemical applications, with band edges well-placed for photocatalytic hydrogen production and CO2 reduction, as well as in tandem dye-sensitized solar cells as a photocathode.
Lithium/chemistry , Nickel/chemistry , Alloys/chemistry , Semiconductors , Solar Energy
The positions of electronic band edges are one important metric for determining a material's capability to function in a solar energy conversion device that produces fuels from sunlight. In particular, the position of the valence band maximum (conduction band minimum) must lie lower (higher) in energy than the oxidation (reduction) reaction free energy in order for these reactions to be thermodynamically favorable. We present first principles quantum mechanics calculations of the band edge positions in five transition metal oxides and discuss the feasibility of using these materials in photoelectrochemical cells that produce fuels, including hydrogen, methane, methanol, and formic acid. The band gap center is determined within the framework of DFT+U theory. The valence band maximum (conduction band minimum) is found by subtracting (adding) half of the quasiparticle gap obtained from a non-self-consistent GW calculation. The calculations are validated against experimental data where possible; results for several materials including manganese(ii) oxide, iron(ii) oxide, iron(iii) oxide, copper(i) oxide and nickel(ii) oxide are presented.
Hematite (α-Fe(2)O(3)) is a promising candidate for photoelectrochemical splitting of water. However, its intrinsically poor conductivity is a major drawback. Doping hematite to make it either p-type or n-type enhances its measured conductivity. We use quantum mechanics to understand how titanium, zirconium, silicon, or germanium n-type doping affects the electron transport mechanism in hematite. Our results suggest that zirconium, silicon, or germanium doping is superior to titanium doping because the former dopants do not act as electron trapping sites due to the higher instability of Zr(III) compared to Ti(III) and the more covalent interactions between silicon (germanium) and oxygen. This suggests that use of n-type dopants that easily ionize completely or promote covalent bonds to oxygen can provide more charge carriers while not inhibiting transport.
Ferric Compounds/chemistry , Models, Chemical , Computer Simulation , Electron Transport , Quantum Theory
The ability to control electronic tunneling in complex molecular networks of multiple donor/acceptor sites is studied theoretically. Our past analysis, demonstrating the phenomenon of site-directed transport, was limited to the coherent tunneling regime. In this work we consider electronic coupling to a dissipative molecular environment including the effect of decoherence. The nuclear modes are classified into two categories. The first kind corresponds to the internal molecular modes, which are coupled to the electronic propagation along the molecular bridges. The second kind corresponds to the external solvent modes, which are coupled to the electronic transport between different segments of the molecular network. The electronic dynamics is simulated within the effective single electron picture in the framework of the tight binding approximation. The nuclear degrees of freedom are represented as harmonic modes and the electronic-nuclear coupling is treated within the time-dependent Redfield approximation. Our results demonstrate that site-directed tunneling prevails in the presence of dissipation, provided that the decoherence time is longer than the time period for tunneling oscillations (e.g., at low temperatures). Moreover, it is demonstrated that the strength of electronic coupling to the external nuclear modes (the solvent reorganization energy) controls the coherent intramolecular tunneling dynamics at short times and may be utilized for the experimental control of site-directed tunneling in a complex network.
Electrons , Models, Molecular , Computer Simulation , Electron Transport , Feasibility Studies , Temperature , Time Factors
