Computational Modeling of the Oxygen Evolution Reaction at Semiconductor-Water Interfaces: A Path Towards Breaking Linear Scaling Relationships

Photoelectrochemical water splitting is a promising source of clean, renewable fuel in the form of hydrogen. Despite extensive research endeavors, the widespread adoption of this technology is impeded due to suboptimal catalysts for the oxygen evolution reaction (OER) taking place on the anode. Studies have shown that linear scaling relationships exist between the binding energies of the OER intermediates. These relationships are responsible for the emergence of an overpotential which limits the catalytic efficiency of electrochemical cells. This work is devoted to the study of the OER, the emergence of linear scaling relationships, and possible avenues for overcoming them. We employ a selection of ab-initio computational free energy methods based on density functional theory in order to investigate the OER at atomistic model semiconductor interfaces.

An approach based on thermodynamic integration and \textit{ab-initio} molecular dynamics simulations is used to isolate the effect of the solvent on the OER free energy steps. The observed effect is twofold. First, the presence of explicit water affects the equilibrium configurations of the reaction intermediates. Second, explicit water modifies the OER free energy steps by up to 0.5 eV. This effect is rationalized in terms of electrostatic interactions at the interface. Next, we verify the existence of the linear scaling relationships across a set of semiconductor and metal interfaces. We establish the robustness of the scaling relationships with respect to the adopted energy functional. Two potential approaches of breaking the linear scaling are considered.

First, a bifunctional mechanism that involves two close but functionally different active sites is investigated. The performance of the proposed mechanism is studied on a set of model interfaces, and the results indicate that specific combinations of catalysts may overcome the limitations imposed by the linear scaling relationship. We were able to establish general trends governing the efficiency of the bifunctional mechanism. In particular, a correlation between the valence band maximum and the hydrogen binding energy may be used as a descriptor in future searches for suitable hydrogen acceptors within the bifunctional scheme. Next, a detailed study of NiOOH/FeOOH catalysts identifies features most favorable to the bifunctional mechanism. Overpotentials as low as 0.13 eV are found for specific configurations. These results support the bifunctional mechanism as a means to break the linear scaling relationships.

Second, hole polarons at the TiO $_2$ surface induced by electronegative adsorbates in the OER intermediates are studied. The computational hydrogen electrode method in conjunction with a hybrid density functional is employed to quantify the effect of the surface polarons on the OER free energy steps. We find that the hole bipolarons reduce the overpotential of the reaction-determining step leading to good agreement with experiment. The stability of the polarons is further confirmed at the hydrated surface through a free energy study involving the explicit treatment of the solvent. Finally, the polaron energy level is aligned with respect to the redox levels in water to better understand the role of the polarons in the OER. Since the occurrence of surface hole polarons is unrelated to the scaling relationships, it offers an additional handle in the search for improved catalysts.

Alignment of energy levels in amorphous oxides and at semiconductor-water interfaces

Zhendong Guo

We conduct a detailed investigation of defects in two representative amorphous oxides: amorphous Al2O3 (am-Al2O3) and TiO2 (am-TiO2), by combining ab initio molecular dynamics (MD) simulations and hybrid functional calculations. Our results indicate that oxygen vacancies and interstitials occur neither in am-Al2O3 nor in am-TiO2 due to structural rearrangements which assimilate the defect structure, and that the injection of extra holes can lead to the formation of O-O peroxy linkages. In am-Al2O3, hydrogen is found to be amphoteric. Based on localized Wannier functions, we identify the nominal charge state and the composition of the defect core units related to C and N impurities in am-Al2O3, which are found to depend on the total charge set in the simulation cell. Through the adopted electron counting rule, we assess that carbon and nitrogen impurities are only found in neutral and in singly positive charge states, respectively, indicating that none of them gives charge transition levels in am-Al2O3. In addition, those defect core units are shown to incorporate a varying number of oxygen atoms, by which their formation energy is dependent on the oxygen chemical potential. In addition, we propose an exchange mechanism for hole transport in am-TiO2 that relies on the simultaneous breaking and forming of O-O peroxy linkages that share one O atom. Through the use of nudged-elastic-band calculations, a hopping path as long as 1.2 nm with barriers of 0.3-0.5 eV is demonstrated, suggesting that hole diffusion through O-O peroxy linkages is viable in am-TiO2. In this work, we also determine the band alignment between various semiconductors and liquid water by combining MD simulations of atomistic interface models, electronic-structure calculations at the hybrid-functional and GW levels, and a computational standard hydrogen electrode (SHE). Our study comprises GaAs, GaP, GaN, CdS, ZnO, SnO2, rutile and anatase TiO2. For each semiconductor, we generate atomistic interface models with liquid water at the pH corresponding to the point of zero charge. The MD simulations are started from initial configurations, in which the water molecules are either molecularly (m) or dissociatively (d) adsorbed on the semiconductor surface. The calculated band offsets are found to be strongly influenced by the adsorption mode at the semiconductor-water interface, leading to differences larger than 1 eV between m and d models of the same semiconductor. We then assess the accuracy of various ab initio electronic-structure schemes in determining the band alignment. In the last part, we try to evaluate photocatalysts for water-splitting by considering the surface coverage and the energy alignment. We determine surface concentrations of water molecules, protons, and hydroxyl ions adsorbed at the semiconductor-water interfaces mentioned above as a function of pH. This is achieved through the calculation of the acidity constants at the surface sites, which are derived from ab initio MD simulations and a grand-canonical formulation of adsorbates. It is finally shown how the potential of a semiconductor material as photocatalysts for water splitting can be inferred by combining the nature of the surface coverage and the alignment of the band edges to the relevant redox levels.

EPFL2019

Coupling of electrocatalysts to photoabsorbers for solar fuels production

Carlos Gilberto Morales Guio

The direct transformation of sunlight into fuels and chemicals has attracted considerable attention for the potential impact on the storage of solar energies at large scales and the reduction of CO2 emissions. More than 80% of our current global energy consumption comes from the burning of non-renewable fossil fuels which produces enormous quantities of CO2 and contributes to global warming. We believe today that the key to reduce our dependence on exhaustible natural resources and limit the emission of greenhouse gases is the transition to CO2-free renewable energy sources. However, harvesting of energy from renewable sources (i.e. solar, wind, tidal) is intermittent and the deployment of devices for their transformation into useful forms of energy (i.e. electricity, chemicals, fuels, biomass) in a global scale is conditioned upon the improvement in efficiency of the energy storage mechanisms. A promising approach to store energy is the use of sunlight to drive the production of hydrogen fuel from water. Hydrogen produced from water and sunlight can be transformed back to electricity on demand or be used to generate mechanical energy in cars to power the transportation industry. The oxidation of hydrogen generates useable energy and at the same time exhaust clean, carbon-free water as the only product. Therefore, the production of hydrogen using a renewable energy source in an efficient and inexpensive device has the potential to discourage the continued use of fossil fuels and improve our environment. The viability of any mechanism for the storage of sunlight as fuel will depend on the reduction of energy penalties during the transformation process. Electrocatalysts are advanced materials that bring into one active site electrons and molecules in the appropriate orientation to lower activation energy barriers, accelerate rate of transformations and improve overall efficiencies for chemical reactions. The best catalysts for water splitting are precious metals such as Pt, Ir and Ru, which show superior rates of reaction at low overpotentials. However, these noble metals are too expensive and scarce for large scale applications. This thesis summarizes our recent research in the preparation, characterization, and application of Earth-abundant electrocatalysts for solar water splitting. Although in the last decade a great number of efficient and inexpensive electrocatalysts have been reported, methods for the deposition of these materials on the surface of photoabsorbers are rarely reported. This work shows various simple and scalable techniques for the deposition of hydrogen and oxygen evolving electrocatalysts on the surface of photocathodes and photoanodes, respectively. Surface-protected photocahodes were activated for solar hydrogen production in alkaline conditions with MoS2+x, Ni-Mo and Mo2C. A novel method for the deposition of an optically transparent amorphous iron nickel oxide oxygen evolution electrocatalyst is also reported. The catalyst was deposited on both thin film and high-aspect ratio nanostructured hematite photoanodes. The low catalyst loading combined with its high activity at low overpotential results in significant improvement on the onset potential for photoelectrochemical water oxidation. This transparent catalyst further enables the preparation of a stable hematite/perovskite solar cell tandem device, which performs unassisted water splitting.

EPFL2016

Solar photoelectrolysis of water with translucent nano-structured hematite photoanodes

Ilkay Cesar

In the context of the diminishing fossil energy resources and global warming, much research is focused on sustainable energy sources. The research that is presented in this thesis falls in this category as it contributes to the development of a device that stores solar energy into a chemical fuel. This device is based on the photoelectrolysis of water in a tandem-cell that produces hydrogen and oxygen under illumination of sunlight. This can be achieved by connecting a larger band gap semiconductor photoanode for oxygen evolution (such as WO3 with Eg = 2.6 eV or Fe2O3 with Eg = 2.0 eV) in series with two dye sensitized solar cells and a hydrogen evolving cathode. Red and green light transmitted by the photoanode is absorbed by the dye, so that a large part of the solar spectrum can be utilized. Fe2O3 absorbs a larger part of the visible solar spectrum as compared to WO3, which results in a 3 times higher theoretical conversion efficiency for the tandem-cell. However, the small hole diffusion length in Fe2O3 of 2 to 20nm is seen as an important reason for the lower practical conversion efficiencies reported in literature. The aim of this thesis project is to improve the photooxidation activity of iron oxide photoanodes in order to improve the practical conversion efficiency of the tandem-cell. For this purpose we have prepared thin films of nano-structured iron oxide with various deposition methods and characterized their performance on the photooxidation of water. The photoanodes are measured in aqueous solution of 1 M NaOH (pH=13.6) under illumination of simulated sunlight (AM 1.5 global of 1000 W/m2). We report the photocurrent at the reversible water oxidation potential of 1.23 V vs. RHE. The first chapter briefly discusses the potential of hydrogen as an energy carrier in a global hydrogen economy and places the tandem-cell in this context. In chapter two the current status of photoelectrolysis of water with iron oxide is briefly reviewed. The experimental set-up for the measurement of the photocurrent under simulated sun light and for the acquisition of photocurrent action spectra are described. This chapter also deals with the spectral mismatch error that is associated with the solar simulator calibration procedure. Chapter three introduces a new solution strategy for efficient photoanodes based on a nanocomposite structure with a hematite film thickness that is commensurate to the hole diffusion length alleviating the problem of poor charge transport. For this purpose nano-porous films of doped and undoped tinoxide and arrays of perpendicularly oriented nanorods of the same material are conformally coated with a thin film of hematite (2-20nm). These electrodes are characterized by SEM, TEM, XEDS elemental analysis, Raman spectroscopy in addition to the photoelectrochemical response. Although, evidence is presented for the successful preparation of the nano-composite structures with a sufficient optical density, the photocurrents obtained from these electrodes are very low. In chapter four we report on thin silicon-doped nanocrystalline α-Fe2O3 films (thickness of 250-500nm) that have been deposited on F-doped SnO2 substrates by ultrasonic spray pyrolysis (USP) and chemical vapor deposition at atmospheric pressure (APCVD). The photoanodes prepared by USP and APCVD gave a photocurrent of 1.17 and 1.45 mA/cm2 respectively. The morphology of the α-Fe2O3 was strongly influenced by the silicon doping, decreasing the feature size of the mesoscopic film. The silicon-doped α-Fe2O3 nano-leaflets show a preferred orientation with the (001) basal plane normal to the substrate. Chapter five treats the preparation and characterization of hematite photoanodes prepared by an improved APCVD setup. Under illumination water is oxidized at the Fe2O3 electrode with higher efficiency (IPCE = 42% at 370 nm and 2.2 mA/cm2) than at the best reported single crystalline Fe2O3 electrodes. This unprecedented efficiency is in part attributed to the dendritic nanostructure which minimizes the distance photo generated holes have to diffuse to reach the Fe2O3/electrolyte interface while still allowing efficient light absorption. Part of the gain in efficiency is obtained by depositing a thin insulating SiO2 interfacial layer between the SnO2 substrate and the Fe2O3 film, and a catalytic cobalt monolayer on the Fe2O3 surface. A mechanistic model for water photooxidation is presented, involving stepwise accumulation of four holes by two vicinal iron or cobalt surface sites. In chapter six we present a detailed investigation of the dependency of the photocurrent response on film thickness and feature size of doped and undoped hematite photoanodes prepared by the APCVD setup discussed in chapter five, in absence of a redox catalyst. Evidence is presented of an unusual high donor density of 1020 cm-3 in the silicon doped photoanodes which would allow a space charge layer to be formed inside the nano-sized features of the polycrystalline film. This suggests that electric field induced charge separation could be an important reason for the photocurrent improvement reported in chapter five. In addition, a strong correlation between feature size and the intensity of a Raman band at 660cm-1 is presented and discussed in the context of disorder in the crystal structure. Chapter seven presents a study of the influence that the deposition parameters have on the photoresponse, morphology, and crystal structure of films prepared by ultrasonic spray pyrolysis. Films are characterized by SEM, XEDS elemental analysis, TEM, XRD, Mössbauer and Raman spectroscopy. A strong correlation is found between the formation of β-Fe2O3 (Bixbyite) and low photocurrents. Silicon doped films enable the electron transport across a film thickness of 1300nm. The absorbed-photon-to-current-conversion-efficiency (APCE) increases with film thickness below 150nm, which indicates an increased recombination near the substrate interface. This could be an important reason for the low photoresponse of nanocomposite electrodes studied in chapter three. In chapter eight we draw some general conclusions with regards to the future development of these photoanodes.

EPFL2007