Challenges of modeling nanostructured materials for photocatalytic water splitting.

Understanding the water splitting mechanism in photocatalysis is a rewarding goal as it will allow producing clean fuel for a sustainable life in the future. However, identifying the photocatalytic mechanisms by modeling photoactive nanoparticles requires sophisticated computational techniques based on multiscale modeling. In this review, we will survey the strengths and drawbacks of currently available theoretical methods at different length and accuracy scales. Understanding the surface-active site through Density Functional Theory (DFT) using new, more accurate exchange-correlation functionals plays a key role for surface engineering. Larger scale dynamics of the catalyst/electrolyte interface can be treated with Molecular Dynamics albeit there is a need for more generalizations of force fields. Monte Carlo and Continuum Modeling techniques are so far not the prominent path for modeling water splitting but interest is growing due to the lower computational cost and the feasibility to compare the modeling outcome directly to experimental data. The future challenges in modeling complex nano-photocatalysts involve combining different methods in a hierarchical way so that resources are spent wisely at each length scale, as well as accounting for excited states chemistry that is important for photocatalysis, a path that will bring devices closer to the theoretical limit of photocatalytic efficiency.

Simulations to Cover the Waterfront for Iron Oxide Catalysis.

Hematite has been widely studied for catalytic water splitting, but the role of the interactions between catalytic sites is unknown. In this paper, we calculate the oxygen evolution reaction free energies and the surface adsorption distribution using a combination of density functional theory and Monte Carlo simulations to "cover the waterfront," or cover a wide range of properties with a simulation of the hematite surface under working conditions. First, we show that modeling noninteracting catalytic sites provides a poor explanation of hematite's slow reaction kinetics. The interactions between the catalytic site may hinder catalysis through the strong interactions of *OH2 and *OOH intermediates, which cause the reaction to revert back to the *O intermediate. Hence, neighboring interactions may be a possible reason for the abundant, experimentally observed *O intermediate on the surface. This study demonstrates how neighboring sites impact the energy required for catalytic steps, thus providing new avenues to improve catalysis by controlling neighboring site interactions.

The effect of interlayer stacking arrangements in two dimensional NiOOH on water oxidation catalysis.

Electrolysis of water to produce green and renewable hydrogen fuel is of great interest in the clean energy field. Water molecules can be decomposed to hydrogen and oxygen through catalysis. Catalytic materials under electrochemical operation are subject to harsh chemical environments, and as a result mechanical changes may appear in the material. In two dimensional materials, the weak van der Waals (vdW) forces holding the layers together may cause a change in the stacking order of the material. The big challenge is to understand the effect of the interlayer arrangements of two dimensional materials on their catalytic performance. In this research we use Density Functional Theory in order to explore the catalytic performance of ß-NiOOH, a two dimensional material that is one of the best known catalysts for the oxygen evolution reaction (OER), under different displacements. Our results indicate that changes in the structural stacking of NiOOH could affect the catalytic properties of the system. Particularly, we find that small shifts between the layers enhance the OER activity by reducing the overpotential down to 240 [mV] due to the formation of an unstable state and the formation of new vdW bonds between the layers. The potential ability to lower the overpotential of NiOOH could give exceptional results in increasing the efficiency of the OER.

Enhanced Li-ion diffusion and electrochemical performance in strained-manganese-iron oxide core-shell nanoparticles.

Efforts to improve energy storage depend greatly on the development of efficient electrode materials. Recently, strain has been employed as an alternate approach to improve ion mobility. While lattice strain has been well-researched in catalytic applications, its effects on electrochemical energy storage are largely limited to computational studies due to complexities associated with strain control in nanomaterials as well as loss of strain due to the phase change of the active material during charging-discharging. In this work, we overcome these challenges and investigate the effects of strain on supercapacitor performance in Li-ion-based energy devices. We synthesize epitaxial Fe3O4@MnFe2O4 (core@shell) nanoparticles with varying shell thickness to control the lattice strain. A narrow voltage window for electrochemical testing is used to limit the storage mechanism to lithiation-delithiation, preventing a phase change and maintaining structural strain. Cyclic voltammetry reveals a pseudocapacitive behavior and similar levels of surface charge storage in both strained- and unstrained-MnFe2O4 samples; however, diffusive charge storage in the strained sample is twice as high as the unstrained sample. The strained-MnFe2O4 electrode exceeds the performance of the unstrained-MnFe2O4 electrode in energy density by ∼33%, power density by ∼28%, and specific capacitance by ∼48%. Density functional theory shows lower formation energies for Li-intercalation and lower activation barrier for Li-diffusion in strained-MnFe2O4, corresponding to a threefold increase in the diffusion coefficient. The enhanced Li-ion diffusion rate in the strained-electrodes is further confirmed using the galvanostatic intermittent titration technique. This work provides a starting point to using strain engineering as a novel approach for designing high performance energy storage devices.

Pathways for charge transport through material interfaces.

Modeling charge transport across material interfaces is important for understanding the limitations of electronic devices such as transistors, electrochemical cells, sensors, and batteries. However, modeling the entire structure and full dimensionality of an interface can be computationally demanding. In this study, we investigate the validity of an efficient reduced one-dimensional Hamiltonian for calculating charge transport along interfaces by comparing to a two-dimensional model that accounts for additional charge transport pathways. We find that the one-dimensional model successfully predicts the qualitative trend of charge transmission probability among Pt/Fe2O3 and Ag/Fe2O3 interfaces. However, the two-dimensional model provides additional information on possible pathways that are not perpendicular to the interface direction. These charge transport pathways are directed along the lowest potential energy profile of the interface that correlates with the crystal structure of the constituting materials. However, the two-dimensional paths are longer and take more scattering time. Therefore, the one-dimensional model may hold sufficient information for qualitative estimation of charge transport through some material interfaces.

Hydrogen transfer through different crystal phases of nickel oxy/hydroxide.

Nickel hydroxide phases are common in several energy conversion devices including battery electrochemical cell electrodes. These materials have a unique layered structure that may facilitate hydrogen transfer between oxygen sites and allow using this material also for proton conducting fuel cells. In order to assess this functionality, we use Density Functional Theory+U together with the Nudged Elastic Band method to calculate minimum energy diffusion paths and hydrogen vacancy formation energies of different crystal phases of NiOOH including ß-NiOOH, ß-Ni(OH)2 and α-Ni(OH)2. We follow several diffusion paths and mechanisms in several phases, both across layers and through them. We pin down the reason for efficient diffusion laterally on a layer and explain why diffusion through a layer is impossible. Our results suggest that hydrogen transfer may be possible for the ß-NiOOH phase with hydrogen added interstitially and transferred along the layers. This study significantly advances our understanding of diffusion in an uncommonly structured material.

Communication: Nickel hydroxide as an exceptional deviation from the quantum size effect.

The quantum size effect is a well-known fundamental scientific phenomenon. Due to quantum confinement, downscaling a system to small sizes should increase the bandgap of a solid state material. However, in this work, we present an exception: monolayers of nickel hydroxide have smaller bandgaps than their bulk analogues, due to the surface states appearing at energies within the bandgap region. Our findings are obtained by several state-of-the-art first principles calculations.

Identifying the bottleneck of water oxidation by ab initio analysis of in situ optical absorbance spectrum.

Hematite's (α-Fe2O3) major limitation to efficiently splitting water using sunlight is the low rate of the oxygen evolution reaction (OER). Thus, identifying the OER rate limiting step is a cornerstone to enhancing the current under low applied potential. Different measurement techniques showed similar absorption difference spectra during a change in applied potential on the hematite anode below and above the onset of the OER in the dark and under light. This absorption change was shown to result from surface modification during the OER, but the specific surface species could not be resolved. On the basis of ab initio calculations, we analyze the calculated absorption spectra in relation to previous measurements. We provide for the first time solid evidence to specify H2O + *O → *OOH + H+ + e- as the rate limiting step and *O as the bottleneck intermediate of the hematite OER.

Play the heavy: An effective mass study for α-Fe2O3 and corundum oxides.

Iron(iii) oxide (α-Fe2O3) is a known water splitting catalyst commonly used in photoelectrochemical cells. These cells are severely impaired by poor conductivity in α-Fe2O3, and resolving the conductivity issue is therefore crucial. One of the most intrinsic properties of matter, which governs conductivity, is the carrier effective masses. In this work, we investigate the carrier effective masses in α-Fe2O3 and other corundum oxides, including Al2O3, Cr2O3, Ga2O3, and In2O3 with different theoretical constructs: density functional theory (DFT), DFT+U, hybrid DFT, and G0W0. We find DFT sufficiently describes the carrier masses and a quasi-particle theory is only required for accuracies better than 30% for the conduction band effective mass. Additionally, we compare the density of states (DOS) and band effective mass approximations and conclude the DOS effective mass provides poor results whenever the band structure is anisotropic. We find that the charge carriers in Fe2O3 "play the heavy" since they have large effective masses that reduce conductivity and device efficiency. Finally, we conclude that the less heavy electron effective masses of other corundum oxides studied relative to Fe2O3 could contribute to efficiency improvements in Fe2O3 upon Al2O3, Ga2O3, and In2O3 coverage.

From density functional theory to machine learning predictive models for electrical properties of spinel oxides.

This work focuses on predicting and characterizing the electronic conductivity of spinel oxides, which are promising materials for energy storage devices and for the oxygen evolution and oxygen reduction reactions due to their attractive properties and abundance of transition metals that can act as active sites for catalysis. To this end, a new database was developed from first principles, including band structure and conductivity properties of spinel oxides, and machine learning algorithms were trained on this database to predict electronic conductivity and band gaps based solely on the compositions. The models developed in this study are scaled from the quantum level up to a continuum conductivity model. The relatively small database used in this study allowed for accurate predictions of band gap and conductivity. By altering the composition of spinel oxides, the model was able to predict high conductivity for spinels with high nickel content and to match experimental trends for manganese cobalt spinels. The ability to predict material properties is especially important in energy conversion devices such as batteries and supercapacitors where redox reactions take place.

Enhancing Photocatalysis: Understanding the Mechanistic Diversity in Photocatalysts Modified with Single-Atom Catalytic Sites.

Surface modification of heterogeneous photocatalysts with single-atom catalysts (SACs) is an attractive approach for achieving enhanced photocatalytic performance. However, there is limited knowledge of the mechanism of photocatalytic enhancement in SAC-modified photocatalysts, which makes the rational design of high-performance SAC-based photocatalysts challenging. Herein, a series of photocatalysts for the aerobic degradation of pollutants based on anatase TiO2 modified with various low-cost, non-noble SACs (vanadate, Cu, and Fe ions) is reported. The most active SAC-modified photocatalysts outperform TiO2 modified with the corresponding metal oxide nanoparticles and state-of-the-art benchmark photocatalysts such as platinized TiO2 and commercial P25 powders. A combination of in situ electron paramagnetic resonance spectroscopy and theoretical calculations reveal that the best-performing photocatalysts modified with Cu(II) and vanadate SACs exhibit significant differences in the mechanism of activity enhancement, particularly with respect to the rate of oxygen reduction. The superior performance of vanadate SAC-modified TiO2 is found to be related to the shallow character of the SAC-induced intragap states, which allows for both the effective extraction of photogenerated electrons and fast catalytic turnover in the reduction of dioxygen, which translates directly into diminished recombination. These results provide essential guidelines for developing efficient SAC-based photocatalysts.

On the Interplay Between Oxygen Vacancies and Small Polarons in Manganese Iron Spinel Oxides.

Ternary spinel oxides are promising materials due to their potentially versatile properties resulting from the disorder inherent in their crystal structure. To fully unlock the potential of these materials, a deeper understanding of their electronic structures, both as pristine and defective crystals, is required. In the present work, we investigate the effects of oxygen vacancies on the electronic structure and charge transport properties of the ternary spinel oxide Mn x Fe3-x O4, modeled on epitaxial thin films of the material, using density functional theory + U (DFT + U). The formation energy of a single oxygen vacancy in the spinel cell is found to be large and unaffected by changes in stoichiometry, in agreement with experimental results. We find that the immediate vicinity of the vacancy has a marked impact on the formation energy. In particular, Mn cations are found to be preferred over Fe as sites for charge localization around the vacancy. Finally, we examine the charge transport in the defective cell using the formalism of Marcus theory and find that the activation barrier for electron small-polaron hopping between sites not adjacent to the vacancy is significantly increased, with a large driving force toward sites that reside on the same (001) plane as the vacancy. Hence, vacancies delay charge transport by increasing the activation barrier, attributed to a rearrangement of vacancy-released charge on the cations immediately neighboring the vacancy site. These results highlight the impact of oxygen vacancies on charge transport in spinel oxides.

NiFe-mixed metal porphyrin aerogels as oxygen evolution reaction catalysts in alkaline electrolysers.

Aerogels are a very interesting group of materials owing to their unique physical and chemical properties. In the context of electrocatalysis, the focus has been on their physical properties, and they have been used primarily as catalyst supports so far. In this work, we synthesized porphyrin aerogels containing Ni and NiFe mixed metal materials and studied them as catalysts for the oxygen evolution reaction (OER). Different Ni : Fe ratios were synthesized and studied in electrochemical cells, and DFT calculations were conducted in order to gain insight into their behavior. The activity trends were dependent on the metal ratios and differ from known NiFeOOH materials due to the change in the oxidation states of the metals to higher numbers. Herein, we show that Ni and Fe have a synergistic effect on the OER, despite being structurally separated. They are connected electronically, though, through a large organic aromatic system that facilitates electron sharing between them. Among the mixed metal porphyrin aerogels, the best ratio was found to be Ni : Fe = 35 : 65, in contrast to oxide/oxyhydroxide materials in which a ratio of 80 : 20 was found to be ideal.

Dipole-Induced Raman Enhancement Using Noncovalent Azobenzene-Functionalized Self-Assembled Monolayers on Graphene Terraces.

Graphene is a promising material in the field of interface science, especially for noncovalent functionalization, sensing, and for applications in catalysis and nanoelectronics. The noncovalent self-assembly of aromatic molecules on graphene promotes electronic coupling through π-π interactions that allows for quenching of the fluorescence of adsorbent molecules and the enhancement of their Raman spectra via graphene-enhanced Raman spectroscopy (GERS). Although recent work has explored the Raman enhancement on mono- and bilayer graphene, the layer dependence of both electronic phenomena (i.e., fluorescence quenching and Raman enhancement) has largely remained underexplored. Similarly, the effect of near-surface molecular dipoles on GERS has sparsely been examined. In this work, we employ self-assembled monolayers of azobenzene-decorated triazatriangulene molecules (AzoTATA) on graphene terraces to examine the effect of switchable molecular dipoles on the GERS effect, which occurs as a function of azobenzene photoisomerization. Furthermore, using empirical and computational methods, we present a systematic study for deriving the mechanism of GERS enhancement and fluorescence quenching on graphene terraces.

Bifunctional PGM-free metal organic framework-based electrocatalysts for alkaline electrolyzers: trends in the activity with different metal centers.

In order to solely rely on renewable and efficient energy sources, reliable energy storage and production systems are required. Hydrogen is considered an ideal solution as it can be produced electrochemically by water electrolysis and renewably while no pollutants are released when consumed. The most common catalysts in electrolyzers are composed of rare and expensive precious group metals. Replacing these materials with Earth-abundant materials is important to make these devices economically viable. Metal organic frameworks are one possible solution. Herein we demonstrate the synthesis and characterization studies of metal benzene-tri-carboxylic acid-based metal-organic frameworks embedded in activated carbon. The conductive composite material was found to be electrocatalytically active for both the oxygen evolution reaction and the hydrogen evolution reaction. Furthermore, several metal organic frameworks sharing the same ligand but with different first-row transition metals (M = Co, Cu, Fe, Mn) were compared, and the trend of their activity is discussed. Cobalt was found to have the highest activity among the studied metal centers, and therefore has the best potential to serve as a bifunctional catalyst for alkaline electrolyzers.

Optimization of Ni-Co-Fe-Based Catalysts for Oxygen Evolution Reaction by Surface and Relaxation Phenomena Analysis.

Trimetallic double hydroxide NiFeCo-OH is prepared by coprecipitation, from which three different catalysts are fabricated by different heat treatments, all at 350 °C maximum temperature. Among the prepared catalysts, the one prepared at a heating and cooling rate of 2 °C min-1 in N2 atmosphere (designated NiFeCo-N2 -2 °C) displays the best catalytic properties after stability testing, exhibiting a high current density (9.06 mA cm-2 at 320 mV), low Tafel slope (72.9 mV dec-1 ), good stability (over 20 h), high turnover frequency (0.304 s-1 ), and high mass activity (46.52 A g-1 at 320 mV). Stability tests reveal that the hydroxide phase is less suitable for long-term use than catalysts with an oxide phase. Two causes are identified for the loss of stability in the hydroxide phase: a) Modeling of the distribution function of relaxation times (DFRT) reveals the increase in resistance contributed by various relaxation processes; b) density functional theory (DFT) surface energy calculations reveal that the higher surface energy of the hydroxide-phase catalyst impairs the stability. These findings represent a new strategy to optimize catalysts for water splitting.

Breakdown of the Small-Polaron Hopping Model in Higher-Order Spinels.

The small-polaron hopping model has been used for six decades to rationalize electronic charge transport in oxides. The model was developed for binary oxides, and, despite its significance, its accuracy has not been rigorously tested for higher-order oxides. Here, the small-polaron transport model is tested by using a spinel system with mixed cation oxidation states (Mnx Fe3- x O4 ). Using molecular-beam epitaxy (MBE), a series of single crystal Mnx Fe3- x O4 thin films with controlled stoichiometry, 0 ≤ x ≤ 2.3, and lattice strain are grown, and the cation site-occupation is determined through X-ray emission spectroscopy (XES). Density functional theory + U analysis shows that charge transport occurs only between like-cations (Fe/Fe or Mn/Mn). The site-occupation data and percolation models show that there are limited stoichiometric ranges for transport along Fe and Mn pathways. Furthermore, due to asymmetric hopping barriers and formation energies, the Mn O h 2 + polaron is energetically preferred to the Fe O h 2 + polaron, resulting in an asymmetric contribution of Mn/Mn pathways. All of these findings are not contained in the conventional small-polaron hopping model, highlighting its inadequacy. To correct the model, new parameters in the nearest-neighbor hopping equation are introduced to account for percolation, cross-hopping, and polaron-distribution, and it is found that a near-perfect correlation can be made between experiment and theory for the electronic conductivity.

Water Oxidation Catalysis for NiOOH by a Metropolis Monte Carlo Algorithm.

Understanding catalytic mechanisms is important for discovering better catalysts, particularly for water splitting reactions that are of great interest to the renewable energy field. One of the best performing catalysts for water oxidation is nickel oxyhydroxide (NiOOH). However, only one mechanism has been adopted so far for modeling catalysis of the active plane: ß-NiOOH(01̅5). In order to understand how a second reaction mechanism affects catalysis, we perform Density Functional Theory + U (DFT+U) calculations of a second mechanism for water oxidation reaction of NiOOH. Then, we use a Metropolis Monte Carlo algorithm to calculate how many catalytic cycles are completed when two reaction mechanisms are competing. We find that within the Metropolis algorithm, the second mechanism has a higher overpotential and is therefore not active even for large applied biases.

Benchmarking Density Functional Theory Based Methods To Model NiOOH Material Properties: Hubbard and van der Waals Corrections vs Hybrid Functionals.

NiOOH has recently been used to catalyze water oxidation by way of electrochemical water splitting. Few experimental data are available to rationalize the successful catalytic capability of NiOOH. Thus, theory has a distinctive role for studying its properties. However, the unique layered structure of NiOOH is associated with the presence of essential dispersion forces within the lattice. Hence, the choice of an appropriate exchange-correlation functional within Density Functional Theory (DFT) is not straightforward. In this work, we will show that standard DFT is sufficient to evaluate the geometry, but DFT+U and hybrid functionals are required to calculate the oxidation states. Notably, the benefit of DFT with van der Waals correction is marginal. Furthermore, only hybrid functionals succeed in opening a bandgap, and such methods are necessary to study NiOOH electronic structure. In this work, we expect to give guidelines to theoreticians dealing with this material and to present a rational approach in the choice of the DFT method of calculation.

Electronic transport through molecular junctions with nonrigid molecule-leads coupling.

The Landauer-type formulation of current through a molecular junction with electronic-nuclear coupling introduced by Troisi et al. [J. Chem. Phys. 118, 6072 (2003)] is generalized to account for the dependence of the molecule-leads coupling terms on the nuclear coordinates. Although this electronic-nuclear coupling is external to the molecule there is no need to extend the molecular subspace when projection operators are employed for calculations of the current through the junction. A test case of a conductor with vibrating contacts to the leads is studied numerically. It is demonstrated that contact vibrations lead to inelastic contributions to the current and to characteristic features in the I-V curve and its derivatives, similar to the ones observed for internal (molecular) electronic-nuclear coupling.