Clarifying the Roles of Oxygen Vacancy in W-Doped BiVO4 for Solar

Jun 28, 2018 - ACS Applied Energy Materials. Li, Wu, Zhang, and ... Enabling Solar Water Oxidation by BiVO4 Photoanodes in Basic Media. Chemistry of ...
0 downloads 0 Views 3MB Size
Subscriber access provided by The University of British Columbia Library

Article

Clarifying the roles of oxygen vacancy in W-doped BiVO4 for solar water splitting Xin Zhao, Jun Hu, Xin Yao, Shi Chen, and Zhong Chen ACS Appl. Energy Mater., Just Accepted Manuscript • DOI: 10.1021/acsaem.8b00559 • Publication Date (Web): 28 Jun 2018 Downloaded from http://pubs.acs.org on July 2, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Energy Materials

Clarifying the roles of oxygen vacancy in W-doped BiVO4 for solar water splitting §

Xin Zhao,‡, Jun Hu,‡*

§, #

§

Xin Yao, Shi Chen, ∥ and Zhong Chen*,§

§

School of Materials Science and Engineering, Nanyang Technological University, 50 Nanyang Avenue, 639798, Singapore, Singapore. School of Physical and Mathematical Sciences, Nanyang Technological University, 637371, Singapore, Singapore



#

School of Chemical Engineering, Northwest University, Xi’an, P. R. China 710069

E-mail: [email protected]; [email protected] ABSTRACT Most oxide semiconductor photoanode materials for water splitting are synthesized in ambient environment. Oxygen vacancy exists in these samples making them intrinsically n-type at the assynthesized state. Oxygen vacancy has been widely reported for enhancing the performance of a photoanode by improving the electron conductivity. Besides the effect on the bulk materials properties, oxygen vacancy also plays an important role in the interfacial charge transfer to electrolyte, on which much less attention has been paid in the past. Herein, we found that although air-annealed W-doped BiVO4 has a higher electron density, lower surface charge transfer resistance, and a slightly better light absorption than the O2-annealed sample, the latter displays a higher photocurrent density. Experimentally we found that the enhanced performance comes from a better charge separation efficiency, despite that the presence of oxygen vacancy does lead to a better charge transfer efficiency. Theoretical calculation finds that there is a localized state formed inside the bandgap in W-doped BiVO4 with oxygen vacancy, which serves as recombination center to reduce its charge separation efficiency. Oxygen vacancy on the V site activates two different kinds of V into reactive sites for improved surface catalysis. Oxygen vacancy also facilitates the adsorption of the OHads, Oads and OOHads involved in a water splitting process, which benefits the surface catalytic process. It is predicted from this study that better performance can be achieved by introducing oxygen vacancy on the surface of a doped BiVO4 and simultaneously avoiding oxygen vacancy in the bulk. The current study provides an important understanding of the roles played by oxygen vacancy in doped photoanode materials.

ACS Paragon Plus Environment

ACS Applied Energy Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

KEYWORDS: BiVO4; Oxygen vacancy; surface catalysis; DFT calculation; reactive sites; water splitting. 1. INTRODUCTION Photoelectrochemical (PEC) water splitting is a promising route to generate hydrogen fuel to solve energy and environmental issues by converting and storing intermittent and inexhaustible solar energy into chemical bonds.1-3 Compared with photocathode semiconductor for hydrogen evolution, the low efficiency in oxygen evolving of photoanode faces a greater challenge to an overall water splitting.4 As a result, much efforts have been made to enhance the performance of photoanode. Most oxide photoanode materials, such as TiO2, WO3, Fe2O3 and BiVO4, are synthesized in ambient environment and contain native oxygen vacancy in the as-synthesized state. The presence of oxygen vacancy is usually seen as an advantage since they endow these oxides with the electron conductivity (n-type) needed for the photoanodes. However, the electron conductivity usually remains at a lower level than expected for water splitting applications. Many researchers have employed elemental doping such as W or Mo in BiVO4, Ti, Sn, Si, W, Nb in Fe2O3 to improve the electron conductivity by increasing the electron density.2, 5-19 An alternative to elemental doping for better electron conductivity is to introduce oxygen vacancy in an oxide semiconductor. It has been reported oxygen vacancies generated by hydrogen treatment can improve the performance of TiO2, WO3 and BiVO4 for PEC water splitting.20-22 Similarly, oxygen vacancies can be created by annealing in a controlled oxygen-deficient atmosphere. This has been demonstrated in Fe2O3 photoanode, which has tremendously improved the photocurrent density to around 3 mA/cm2.23 The roles of oxygen vacancies can be viewed from two perspectives, the bulk and surface. For the bulk material, free electrons are generated with the creation of oxygen vacancies, which can be seen as electron donors.24 Thus, proper amount of oxygen vacancy can efficiently enhance the conductivity and improve the PEC performance. On the other hand, oxygen vacancy in the bulk also become the recombination center of the photo-generated electrons and holes, therefore too much oxygen vacancy would harm the photocatalytic performance.25 Hence, a good performance is a result of proper balance of the electron

ACS Paragon Plus Environment

Page 2 of 26

Page 3 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Energy Materials

carrier density and charge recombination at the oxygen vacancy site. For the surface, oxygen vacancy plays an important role in the interfacial charge transfer by modulating many surface properties such as water adsorption sites and active reaction sites.26-27 Recently, our theoretical work found that oxygen vacancy has two roles in enhancing the surface catalysis on pristine BiVO4. First, oxygen vacancy increases the active sites by turning the V site into active site. Second, stronger adsorption of H2Oads, OHads, and Oads are observed in the presence of oxygen vacancy, which indicates an easier charge transfer to the electrolyte and improved surface catalysis.26 In the Fe2O3 system, it was found from experiment that oxygen vacancy plays different roles on the surface catalysis for pristine and Ti-doped Fe2O3: it contributes positively to the activity of pristine Fe2O3 but negatively affects the surface activity of Tidoped Fe2O3.27 The reason lies with the fact that oxgen vacancy has introduced an acitve surface state in pristine Fe2O3, which has greatly enhanced the surface hole transfer. However, for Ti-doped sample, Ti doping generated two active surface states and oxygen vacancy had some negative effect on one of the states, which has resulted in a lower charge transfer efficiency of Ti-doped Fe2O3 with oxygen vacancy compared with that without oxygen vacancy.27 Hence, oxygen vacancy plays different roles in different materials. Based on previous studies, we believe it is of great interest and importance to fully understand the roles of oxygen vancancy in a doped BiVO4 photoanode materials. BiVO4 is one of the most promising photoanode materials for PEC water-splitting, which has a theoretical conversion efficiency of 9.1% with a bandgap of 2.4 eV.8, 10, 28 BiVO4 is an n-type semiconductor materials, oxygen vacancies exist when annealed in ambient condition. However, BiVO4 has a poor electron conductivity, thus n-type doping (W or Mo for most cases) is necessary to increase its PEC activity.5-8, 10 In this study, we employed both experimental and theoretical approaches to investigate the roles of oxygen vacancy in a W-doped BiVO4. The W-doped BiVO4 thin films were prepared by a drop-casting method reported by our previous work.29 Heat treatment in ambient condition and pure oxygen environment was conducted to generate different state in oxygen vacancy and to study its effect. Oxygen vacancies can be produced under oxygen-

ACS Paragon Plus Environment

ACS Applied Energy Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

deficient conditions or in a reductive environment such hydrogen.30-31 Hydrogen treatment method has the potential influence of hydrogen impurities.32 And annealing in N2 or Ar gas produced black BiVO4 based on our past experience. Thus, in this study, we use air and pure oxygen environment to generate and compare oxygen vacancy effect, considering the difference of O2 content in air (about 21%) and pure oxygen (99.99% used in this study). The effects on bulk properties like electronic and optical properties and surface catalysis of oxygen vacancy in W-doped BiVO4 were analysed. It was found from experiment that oxygen vacancy decreased the charge separation and improved the surface charge transfer activity. The mechanisms were then understood through theoretical computation. These results enrich our understanding of the roles that oxygen vacancy plays in W-doped BiVO4 photoanode on PEC performance.

2. METHODS AND CHARACTERIZATION 2.1 Preparation of W-doped BiVO4 thin films W-doped BiVO4 thin films were fabricated by a drop-casting method reported before.29 Bismuth trioxide, ammonium metavanadate, and ammonium tungstate hydrate were dissolved in DI water with proper amount of nitric acid to form 0.1 M Bi, V and W aqueous precursor solutions respectively. The solutions were mixed well according to stoichiometric ratio of Bi:V:W=100:97:3, which corresponded to a volume ratio of 500 µl : 485 µl : 15 µl) for 3% W-doped BiVO4. Citric acid (CA) was also added according to stoichiometric ratio of CA:M=1.5:1 (M is the total amount of cation). After that, the mixture was added with 1 ml water and 1 ml ethylene glycol (EG). After mixing, 60 µl of the precursor solution were dropped cast onto the exposed half side of a 1 cm × 2 cm FTO substrates with the other half of the substrate covered by a thermal tape. The samples were then dried at 120°C for 30 minutes to obtain gel films. After that, the tape was peeled off, and the sample was placed into a tube furnace for heat treatment at 500 °C for 2 h with a ramping rate of 10 °C/min. The annealing was carried out in air or pure O2 gas

ACS Paragon Plus Environment

Page 4 of 26

Page 5 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Energy Materials

environment. Before the annealing, the tube furnace was first vacuumed, after which the required gas was flowed through at a rate of 50 cc/min. 2.2 Characterization The morphologies and thickness of the W-doped BiVO4 films treated in air and O2 gas were observed using field emission scanning electron microscopy (FESEM, JEOL JSM-7600F). The light absorption efficiencies of the obtained samples were performed by measuring the reflectance and transmittance spectra using UV-Vis-NIR spectrophotometer with an integrating sphere (Lambda 950, Perkin-Elmer). Crystallinity of the W-doped BiVO4 films were identified by X-ray diffraction (XRD, Shimadzu 6000 Xray diffractometer) with Cu Kα radiation (λ = 0.154 nm). The chemical state of O was investigated by Xray photoelectron spectroscopy (XPS, Omicron EA125). The binding energy was calibrated by C1s (284.6 eV). PEC measurements of the W-doped BiVO4 films were measured using a three-electrode configuration (PCI4/300™ potentiostat with PHE200™ software, Gamry Electronic Instruments, Inc.), with the W-doped BiVO4 film as the working electrode, Pt mesh as the counter electrode, and Ag/AgCl as the reference electrode. A solar simulator (HAL-320, Asahi Spectra Co., Ltd.) with power intensity of 100 mW·cm-2 calibrated by a solar reference cell was used as light source for PEC measurement. The photocurrents of water oxidation were measured in 0.5 M Na2SO4 aqueous solution with a scan rate of 30 mV·s-1. The photocurrents tested with hole scavenger were measured in 0.5 M Na2SO4 aqueous solution mixed with 0.1 M Na2SO3. Na2SO3 was selected as the hole scavenger rather than others because it has a relatively large electrochemical window. Electrochemical impedance spectroscopy (EIS) measurements were carried out under AM 1.5G solar simulator illumination in 0.5 M Na2SO4 electrolyte at the applied potential of 1.23 V vs. RHE using an AUTOLAB potentiostat-galvanostat (AUTOLAB PGSTAT302 N). The Mott-Schottky plots were also carried out using the same AUTOLAB potentiostat-galvanostat at a fixed frequency of 1 kHz in 0.5 M Na2SO4 solution. 2.3 Calculation methods

ACS Paragon Plus Environment

ACS Applied Energy Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The CASTEP module of the Materials Studio software (Accelrys Inc.) was employed for the quantum chemistry calculations. During the calculations, a 2×1×2 BiVO4 (space group 15) supercell with crystallographic parameters of 14.602 Å×11.704 Å×10.290 Å was used for bulk property calculation. One V atom is replaced by W atom to represent the W doped ms-BiVO4. Two different kinds of oxygen vacancies were simulated. One is the one connects V and Bi (V-O-Bi), denoted as OV on V. The other one is the one connects W and Bi (W-O-Bi), denoted as OV on W. 2×2 supercells of W doped BiVO4 with and without oxygen vacancy on the (010) surface with a vacuum region of 15 Å were selected to investigate the surface property. The (010) surface was chosen since it is the most stable surface under the realistic conditions.33-35 Self-consistent periodic density functional theory (DFT) calculations were used. The generalized gradient approximation (GGA), in the form of the Perdew-Burke-Ernzerhof (PBE) approximation, was used to calculate the exchange-correlation energy. Broyden-Fletcher-GoldfarbShanno (BFGS) scheme was chosen as the minimization algorithm, and the DFT-D correction was used for dispersion corrections. Hubbard U-corrections to the d electrons of V (LDA+U, effective U(V)=2.5 eV) and spin-polarized were performed during the calculations. The energy cutoff is 380 eV and the SCF tolerance is 1.0×10-6 eV/atom. The optimization finishes when the energy, maximum force, maximum stress and maximum displacement are smaller than 5.0×10-6 eV/atom, 0.01 eV/ Å, 0.02 GPa and 5.0×10-4 Å, respectively. 2×2×3 k-points samplings were used for bulk and 2×1×1 for surface calculations. The Fermi level is simply defined as the valence band maximum (VBM) for semiconductors and insulators in the CASTEP code and some other codes.36-37

3. RESULTS 3.1 Experimental results

ACS Paragon Plus Environment

Page 6 of 26

Page 7 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Energy Materials

Figure 1. SEM images of W-doped BiVO4 (a) heated in air and (b) heated in pure O2. (c) and (d) are cross section images of W-doped BiVO4 heated in air (denoted as Air) and pure O2 (denoted as O2) respectively. (e) X-ray diffraction of W-doped BiVO4 heated in air and pure O2 respectively. (f) O 1s XPS spectroscopic spectra (intensity normalized for easy comparison) of the W-doped BiVO4 samples obtained by annealing in different atmosphere. Figure 1a-b illustrates the surface structure of W-doped BiVO4 with the heat treatment in air (denoted as sample Air) and pure O2 (denoted as sample O2). Both samples show similar porous structures. The size of pores in sample Air is slightly larger than that in sample O2, while the amount of pores in sample Air is slightly smaller than that in sample O2. Thus, it is difficult to directly tell which sample has a larger surface area. Considering this, electrochemically active surface areas were measured from the capacitive

ACS Paragon Plus Environment

ACS Applied Energy Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

region of cyclic voltammograms (Figure S1). It is found that the electrochemically active surface areas for both samples are almost the same (Figure S1). The cross section images (Figure 1c-d) further prove that the films are composed of porous structure with the thickness of around 400 nm, and adhere strongly to the FTO substrates. No obvious difference was observed for the W-doped BiVO4 heated in different gas atmosphere. X-ray diffraction (Figure 1e) shows all peaks of obtained W-doped BiVO4 thin films agree well with the ones of BiVO4 (PDF#14-0688). No peaks belonging to other phases are present except the ones from the substrate FTO (SnO2). This demonstrates different gas atmosphere heat treatment has no obvious changes on the morphology and crystallinity of W-doped BiVO4. The chemical state of O was investigated by XPS as shown in Figure 1f for the O1s. The low binding energy component located at 526.6 eV is assigned as lattice oxygen (Figure S2), which are almost the same for both samples.38-39 For the air treated sample, there is a peak at the higher binding energy, which corresponds to oxygen vacancy.38-39 Thus, it confirms that there are some oxygen vacancies in the air treated W-doped BiVO4. Compared with the air treated sample, pure O2 gas treated sample has no such peak, which shows that Wdoped BiVO4 heated in pure O2 gas are free of oxygen vacancy (or minimized at least).

Figure 2. Photocurrents of samples treated in air and O2. Figure 2 shows the photocurrents of samples treated in air and O2 after optimizing the annealing

time in air and O2 considering that annealing time might have great impact on the oxygen vacancy

ACS Paragon Plus Environment

Page 8 of 26

Page 9 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Energy Materials

concentration (Figure S3). The sample treated in air has a lower photocurrent (around 1.3 mA/cm2 at 1.23 V vs. RHE), while the sample treated in O2 has relatively better performance (around 1.6 mA/cm2 at 1.23 V vs. RHE). This corresponds to a ~23% increase. In order to confirm the observed enhancement, the experiment was repeated for 3 times with consistent trend observed (Figure S4). At 1.8 V vs. RHE, a ~33% increase is observed (current from 2 mA/cm2 to 2.66 mA/cm2). The observed enhancement in photocurrent density is believed to be attributed to the oxygen vacancy, in view of early evidence of almost identify film morphology, thickness, and crystallinity. Based on the mass action law on the equilibrium Equation (1), oxygen vacancy concentration, ݊௖ , in W-doped BiVO4 is dependent on oxygen −1

partial pressure , ‫݌‬ைమ , according to nc ∝ pO2 6 .24

OO ⇔

1 O2 + 2e '+ VO⋅⋅ 2

(1)

Under air treatment, oxygen vacancies will be created due to the lower oxygen partial pressure compared pure O2 gas. According to Equation (1), free electrons are created with oxygen vacancy generation. Thus, oxygen vacancy is a shallow donor which provides free electrons as charge carriers.23, 25, 40-41 Based on this analysis, the difference between the samples treated in air and O2 should be reflected in the carrier density.

Figure 3. (a) Mott-Schottky plots of sample Air and O2 measured at a fixed frequency of 1 kHz in 0.5M Na2SO4 aqueous solution. (b) Electrochemical impedance spectra (EIS) of sample Air and O2 measured at

ACS Paragon Plus Environment

ACS Applied Energy Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 26

the applied potential of 1.23 V vs. RHE in 0.5M Na2SO4 electrolyte under simulated solar illumination. The inset in Figure 3b is the equivalent circuit to simulate the EIS. Rbulk represents total resistance from external circuit; Q is the constant phase element; Rct is charge transfer resistance of holes from film surface to solution.

Mott-Schottky plots of sample Air and O2 are obtained to investigate the carrier density (Figure 3a). The majority carrier (i.e., electron) density can be estimated by the following Equation (2):

Nd =

−1 2 d (1/ C 2 ) / dV  A e0εε 0 2

(2)

where A is electrode area (1 cm2 in this study), e0 the electron charge (1.60×10-19 C), ε the dielectric constant of BiVO4 (68), ε0 the permittivity of vacuum (8.85×10-12 F·m-1), Nd the donor density and V the potential applied at the electrode. The slope of sample Air is lower than sample O2, which indicates a higher carrier density of the sample treated in air. The enhanced carrier density of sample Air is consistent with the results told from XPS, which proves that air treat sample have created oxygen vacancy as shallow donors and therefore increase the donor density. The calculated carrier density of sample Air and O2 are 6.1×1020 cm-3 and 4.8×1020 cm-3 respectively according to Equation (2). Generally higher carrier density improves the conductivity and a better charge separation efficiency,27 which leads to a higher photocurrent performance. However, in this study sample O2 with a lower carrier density has a better photocurrent than sample Air. Thus, surface properties have to be investigated to understand the charge transfer process next. Electrochemical impedance spectra (EIS) was then measured. In order to reflect the real situation of the surface reaction, EIS was performed under simulated solar illumination. Figure 3b shows the Nyquist plots, in which the radius of the semicircle for the sample O2 is larger than that of sample Air, indicating a larger surface resistance sample O2. According to previous work, an equivalent circuit is taken to analyze the surface charge transfer (inset in Figure 3b).42 In the equivalent circuit, Rs represents the resistances of

ACS Paragon Plus Environment

Page 11 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Energy Materials

the FTO film, the external electrical contacts and the liquid electrolyte, Rct (charge transfer resistance) and Q (constant phase element) represents the direct charge transfer at the semiconductor/electrolyte interface. Ideally, an equivalent circuit Rs (Cbulk Rct) is taken to simulate the EIS. However, non-typical semicircle is observed in this study and thus a constant phase element (Q), representing a non-intuitive circuit element Q, is utilized to simulate such type of spectra.43 The fitting results show that Rs values are independent on the film annealing environment, which are around 40 Ω. The difference lies with the surface charge transfer resistance Rct. Rct of sample Air is about 1410 Ω, while that of sample O2 is larger at around 2224 Ω. This result also indicates that sample O2 shows no advantage on the surface charge transfer over sample Air. Hence, both Mott-Schottky and EIS evidences indicate unfavorable trend for sample O2. In order to gain a clear understanding of this apparent parodox, we have to analyze in details the light absorption efficiency ηabs, charge separation efficiency in the bulk of the film ηsep, and interfacial charge transfer efficiency ηtran in a PEC process. The product of the three efficiencies will be able to explain the observed photocurrent densities.

ACS Paragon Plus Environment

ACS Applied Energy Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 26

Figure 4. (a) Light absorption efficiency (b) charge separation efficiency and (c) charge transfer efficiency of sample Air and O2.

The water oxidation photocurrent JH2O is expressed by:44 J H 2O = J 0 ×η abs ×η sep ×ηtran = J abs ×η sep ×ηtran

(3)

where J0 is the theoretical solar photocurrent density which assumes that all the solar energy corresponding to the band edge can be fully converted to fuel energy (7.3 mA/cm2 for BiVO4).45 Water oxidation has a slow oxidation kinetics. When an effective hole scavenger Na2SO3 is utilized, the charge transfer efficiency can be approximated to be 100% (ηinj=1) due to the fast oxidation kinetics of Na2SO3.44 The Na2SO3 oxidation photocurrent can be determined by: 44

J Na2 SO3 = J 0 ×ηabs ×ηsep = J abs ×ηsep

(4)

where JNa2SO3 is the oxidation photocurrent using Na2SO3. Combining Equations (3) and (4), we can obtain the charge separation efficiency ηsep=JNa2SO3/Jabs and the charge transfer efficiency ηtran=JH2O/JNa2SO3. Figure 4a shows the light absorption efficiency of sample Air and O2 gas by measuring the light transmittance and reflectance.45 The sample Air shows slightly enhanced light absorption compared with sample O2, which is probably due to the presence of oxygen vacancy. Jabs was calculated by integrating the light absorption efficiency over the solar spectrum (AM 1.5 G). The wavelength was integrated up to 510 nm, the light absorption edge of BiVO4. Jabs is estimated to be 4.56 mA/cm2 and 4.24 mA/cm2 for sample Air and O2, respectively. The calculated charge separation efficiency, ηsep, of sample O2 shows about 37% increase at 1.23 V vs. RHE compared with sample Air, while the charge transfer efficiency, ηtran, has a smaller decrease of about 12% at 1.23 V vs. RHE for sample O2 compared with sample Air. Thus, the good charge separation is the main reason for the observed enhanced photocurrent for the sample annealed in O2. In addition, it is also found that the oxygen vacancy can facilitate the surface charge transfer, which is consistent with predicted results from first principle calculation in our previous

ACS Paragon Plus Environment

Page 13 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Energy Materials

work in pristine BiVO4. The theoretical calculation indicates that the adsorption energies of some intermediates such as H2Oads, OHads, and Oads are higher on the surface with oxygen vacancy, which will enhance holes transfer from the pristine BiVO4 surface to the electrolyte.26 Considering that the air atmosphere will not necessarily result a lower amount of oxygen vacancy, the O2 annealed sample was treated in the N2 gas at 300 °C for 2 hours. The photocurrent dropped compared with that before N2 treatment (Figure S5). The charge injection efficiency of N2 treated sample was also calculated as shown in Figure S6, which is similar with that of Air annealed sample. This confirms that oxygen vacancy has positive effect on the surface catalysis. 3.2 Theoretical calculations

ACS Paragon Plus Environment

ACS Applied Energy Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 5. (a) Crystal structure of W-doped BiVO4, W-doped BiVO4 with oxygen vacancy (Ov) on V atom, and W-doped BiVO4 with oxygen vacancy (Ov) on W atom. Bi, V, O, W atoms are represented by purple, grey, red, and blue colours. (b) Density of States (DOS) of W-doped BiVO4, W-doped BiVO4 with oxygen vacancy (Ov) on V atom, and W-doped BiVO4 with oxygen vacancy (Ov) on W atom. (c) Optical absorption spectra of W-doped BiVO4, W-doped BiVO4 with oxygen vacancy (Ov) on V atom, and Wdoped BiVO4 with oxygen vacancy (Ov) on W atom on the (100) plane derived from frequency dependent imaginary part of dielectric function. The 2 × 1 × 2 supercell model of W-doped BiVO4, W-doped BiVO4 with oxygen vacancy (Ov) on V atom (the missing oxygen atom is the one connects V and Bi, V-O-Bi), and W-doped BiVO4 with oxygen vacancy (Ov) on W atom (the missing oxygen atom is the one connects W and Bi, W-O-Bi are shown in Figure 5a. Monoclinic BiVO4 is composed of a layered structure along the crystallographic c axis (Figure 5a), in which two polyhedra units, VO4 tetrahedra and BiO8 dodecahedra link to each other by edge sharing. The doping element, W, substitutes V shown in Figure 5a which is represented in blue color. As seen from Figure 5a, the generation of Ov causes some deformation around the Ov site, around which the well overlapped O atoms (red spheres) in the W-doped BiVO4 now are separated. The affected area is circled in the two Ov contained cells. Figure 5b shows the DOS of W-doped BiVO4, W-doped BiVO4 with Ov on V atom, and W-doped BiVO4 with Ov on W atom. The valence band (VB) and conduction band (CB) are mainly composed of O 2p states and V 3d states. The Bi 6p and 6s states also contribute to the composition of VB and CB. Compared with W-doped BiVO4, one feature is that the electron density near the bottom of the conduction band is enhanced for the Ov containing W-doped BiVO4. This finding is consistent with the enhanced carrier density measured by Mott-Schottky. The other obvious feature is the formation of localized, filled states inside the bandgap, which are composed of O 2p and V 3d states. Such kind of localized states usually serves as recombination centers. Figure 5c shows optical absorption spectra of W-doped BiVO4, W-doped BiVO4 with oxygen vacancy (Ov) on V atom, and W-doped BiVO4 with oxygen vacancy (Ov) on W atom. From Figure 5c, oxygen vacancy can enhance the light absorption

ACS Paragon Plus Environment

Page 14 of 26

Page 15 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Energy Materials

due to more free electrons provided by oxygen vacancy, which is consistent with the experimental measurement of light absorption. Such kind of light absorption change was also reported in H2-treated SrTaO2N.46 Considering the anisotropic optical property of BiVO4, the light absorption caused by oxygen vacancy on the (010) plane was also calculated (Figure S7), the same trend was observed.

ACS Paragon Plus Environment

ACS Applied Energy Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 6. (a) Side view, (b) top view, (c) HOMO and (d) LUMO of W-doped BiVO4 (010) facet and Wdoped BiVO4 (010) facet with oxygen vacancy. The optimized W-doped BiVO4 (010) facet with and without oxygen vacancy are shown in Figure 6. The loss of oxygen between Bi and V (denoted as Ov on V) has negligible influence to the surrounding atoms. The loss of oxygen between Bi and W (denoted as Ov on W) causes the change of the position of O(2) atom, which becomes top site from bridge site, but the bond length has no obvious change. While other atom positions are not affected by the Ov on W. According to the frontier molecular orbital theory of chemical reactivity by Fukui, most of the electron transition will happen at the frontier orbitals, named the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) of reacting species.26, 47 As shown in Figure 6, the HOMO region of W-doped BiVO4 (010) facet is mainly distributed on one V atom while the LUMO region of W-doped BiVO4 (010) facet is also distributed on the V atom but disperses to other V atoms (i.e., not localized on one V atom). Similar observation can be seen when introducing Ov on W. However, for Ov on V, LUMO has little change, but for HOMO one more V atom on the top atom layer also becomes a distribution sit. As well known, the energy of the HOMO orbital is directly related to the ionization potential and characterizes the susceptibility of the molecule to attacks by electrophiles.26 Thus, HOMO distribution site has electron acceptor ability from intermediates of H2O such as O, OH, OOH. From this view, W-doped BiVO4 (010) facet with Ov on V has more active sites for water oxidation. The PDOS of W-doped BiVO4 (010) with and without Ov (Figure S8) also shows that filled states inside the bandgap are observed, which is contributed by V 3d and O 2p orbitals. As seen in Figure S8, two different states in the bandgap emerge for the Ov on V atom, while only one state forms for W-doped BiVO4 without vacancy and W-doped BiVO4 with Ov on W atom. These calculation results indicate that oxygen vacancy on V site will activate two different kinds of V atoms into reactive sites for surface catalysis.

ACS Paragon Plus Environment

Page 16 of 26

Page 17 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Energy Materials

Figure 7. Adsorption structures and energies of OHads, Oads and OOHads involved in splitting processes on (a) W-doped BiVO4 (010) facet, (b) W-doped BiVO4 (010) facet with Ov on W site and (c) W-doped BiVO4 (010) facet with Ov on V site. The “-” and “@” signs stand for bond and adsorption state on the facet, respectively. Figure 7 displays the adsorption of several key species involved in the water oxidation on W-doped BiVO4 (010) facet with and without Ov. The adsorption energies of OHads, Oads and OOHads on Bi atoms of the (010) facet of oxygen vacancy containing W-doped BiVO4 are much smaller than those of W-

ACS Paragon Plus Environment

ACS Applied Energy Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

doped BiVO4 without oxygen vacancy, especially Oads, which indicates that presence of oxygen vacancy favors the adsorption of the surface species. For V atom sites, the adsorption energies of Oads and OOHads show a decrease in the presence of oxygen vacancy, which indicates a better adsorption of these species. For the OHads, the situation is complicated. The adsorption energy increases for the oxygen vacancy generated between W atom and Bi atom, which means a poorer adsorption under this condition. However, if the oxygen vacancy is generated between V and Bi, it facilitates the adsorption of OHads. Considering that the doping level is small (3%), most generated oxygen vacancies should be adjacent to V atoms. Thus, oxygen vacancy will favor the adsorption of all the surface species. For the W sites, the adsorption energies of OHads, Oads and OOHads on W atoms also decrease when oxygen vacancy is introduced on W site, indicating a better adsorption. The comparison of adsorption of surface species involved in the water oxidation apparently shows that oxygen vacancy favors the surface catalytic process.

4. DISCUSSION

Figure 8. Schematic description of the PEC water splitting process on a W-doped BiVO4 with or without oxygen vacancy. Based on the above analysis, the mechanism of how oxygen vacancy affects the PEC water splitting in a W-doped BiVO4 photoanode, as shown in Figure 8. For the bulk charge separation, the oxygen vacancy

ACS Paragon Plus Environment

Page 18 of 26

Page 19 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Energy Materials

plays a negative role, as it serves as charge recombination center. The generation of oxygen vacancy usually improves the electron conductivity that favors charge separation.27 However, its role in hindering charge separation through charge recombination should not be neglected when selecting dopant element and concentration. In the current study, W doping has already improved the electron conductivity in BiVO4 sufficiently; additional oxygen vacancy mainly promotes the charge recombination. If only oxygen vacancy is generated in undoped BiVO4, the photocurrent density improves with higher charge separation and injection efficiencies (Figure S9).48 Hence, the good PEC performance of oxygen vacancy containing sample is the balance of good conductivity and recombination through vacancies. Similar phenomenon exists in other metal oxides such as hematite. It was reported that the photocurrent of hematite first enhanced with increasing oxygen vacancy then showed decrease when the oxygen vacancy concentration is too much.25, 41 For the surface catalysis, two kinds of oxygen vacancies are studied: the one between Bi and V, and the other between Bi and W. The calculated results show that Ov on V can activate the surface V atom into a reactive site, which increase the surface reaction sites. This is likely to be the dominant type of oxygen vacancies, because the doping level of W is very small. Thus, air treated W-doped BiVO4 will favors the surface catalysis of water splitting. In addition, oxygen vacancy influences the adsorption of several key species on W-doped BiVO4 (010) facet, which plays a key role in enhancing the hole transfer for water oxidation. The calculation well explains the observed lower hole transfer resistance and consequent enhanced charge transfer efficiency in experiment. Moreover, we also compared the pristine BiVO4, which also showed that oxygen vacancy containing pristine BiVO4 had a higher charge injection efficiency (Table S1). Thus, for both doped and undoped BiVO4, oxygen vacancy is favorable in the surface catalysis. In summary, oxygen vacancy serves as recombination center in the bulk but facilitates the surface catalysis. Based on the above analysis, it can be predicted that in a doped BiVO4 photoanode better

ACS Paragon Plus Environment

ACS Applied Energy Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

performance can be achieved by introducing oxygen vacancy on the surface and simultaneously avoiding oxygen vacancy in the bulk.

5. CONCLUSION Porous W-doped BiVO4 thin films are prepared by drop casting, followed by heat treatment in air and in O2. The heat treatment does not change the morphology and crystallinity of the thin films, but the airannealed sample contains oxygen vacancy in this intrinsically n-type semiconductor. Air-annealed sample shows better electron density, lower surface charge transfer resistance, and a slightly higher light absorption than the O2-annealed sample. However, the O2-annealed sample displays better photocurrent density than the air-treated sample, as it has a much better charge separation efficiency. Theoretical calculation confirms that the better charge transfer in the air-annealed sample is attributed to the presence of oxygen vacancy for enhanced surface catalysis. Oxygen vacancy enhances the electron density near the bottom of the conduction band, and simultaneously improves light absorption. Oxygen vacancy also facilitates the adsorption of the OHads, Oads and OOHads involved in a water splitting process, which benefits the surface catalysis process. However, as oxygen vacancy acts as a recombination center for the photo-generated electrons and holes, care should be exercised when choosing this strategy to improve the photoanode materials for a PEC performance, particularly in a higher-valence doped sample where the electron density is high enough to ensure good charge separation efficiency. Current study predicts that in a doped BiVO4 photoanode better performance can be achieved by introducing oxygen vacancy on the surface and simultaneously avoiding oxygen vacancy in the bulk. ASSOCIATED CONTENT Supporting Information The Supporting Information is available free of charge on the ACS Publications website.

AUTHOR INFORMATION

ACS Paragon Plus Environment

Page 20 of 26

Page 21 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Energy Materials

‡These authors contributed equally. The authors declare no competing financial interest. ACKNOWLEDGEMENTS Financial supports from Ministry of Education (grant RG15/16), Nanyang Technological University in form of SUG, Singapore National Research Foundation through the Singapore-Berkeley Initiative for Sustainable Energy (SINBERISE) CREATE Programme, the National Natural Science Foundation of China (No. 21676216).

REFERENCES (1) Walter, M. G.; Warren, E. L.; McKone, J. R.; Boettcher, S. W.; Mi, Q.; Santori, E. A.; Lewis, N. S. Solar water splitting cells. Chem. Rev. 2010, 110, 6446-6473. (2) Li, Z.; Luo, W.; Zhang, M.; Feng, J.; Zou, Z. Photoelectrochemical cells for solar hydrogen production: current state of promising photoelectrodes, methods to improve their properties, and outlook. Energy Environ. Sci. 2013, 6, 347-370.

(3) Grätzel, M. Photoelectrochemical cells. Nature 2001, 414, 338-344. (4) Rao, P. M.; Cai, L.; Liu, C.; Cho, I. S.; Lee, C. H.; Weisse, J. M.; Yang, P.; Zheng, X. Simultaneously efficient light absorption and charge separation in WO3/BiVO4 Core/shell nanowire photoanode for photoelectrochemical water oxidation. Nano Lett. 2014, 14, 1099-1105. (5) Rettie, A. J.; Lee, H. C.; Marshall, L. G.; Lin, J.-F.; Capan, C.; Lindemuth, J.; McCloy, J. S.; Zhou, J.; Bard, A. J.; Mullins, C. B. Combined charge carrier transport and photoelectrochemical characterization of BiVO4 single crystals: intrinsic behavior of a complex metal oxide. J. Am. Chem. Soc. 2013, 135, 11389-11396. (6) Luo, W.; Li, Z.; Yu, T.; Zou, Z. Effects of surface electrochemical pretreatment on the photoelectrochemical performance of Mo-doped BiVO4. J. Phys. Chem. C 2012, 116, 5076-5081. (7) Luo, W.; Wang, J.; Zhao, X.; Zhao, Z.; Li, Z.; Zou, Z. Formation energy and photoelectrochemical properties of BiVO4 after doping at Bi3+ or V5+ sites with higher valence metal ions. Phys. Chem. Chem. Phys. 2013, 15, 1006-1013.

(8) Luo, W.; Yang, Z.; Li, Z.; Zhang, J.; Liu, J.; Zhao, Z.; Wang, Z.; Yan, S.; Yu, T.; Zou, Z. Solar hydrogen generation from seawater with a modified BiVO4 photoanode. Energy Environ. Sci. 2011, 4, 4046-4051.

ACS Paragon Plus Environment

ACS Applied Energy Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(9) Abdi, F. F.; Firet, N.; van de Krol, R. Efficient BiVO4 thin film photoanodes modified with Cobalt Phosphate catalyst and W-doping. ChemCatChem 2013, 5, 490-496. (10) Abdi, F. F.; Han, L.; Smets, A. H.; Zeman, M.; Dam, B.; van de Krol, R. Efficient solar water splitting by enhanced charge separation in a bismuth vanadate-silicon tandem photoelectrode. Nat. Commun. 2013, 4, 2195.

(11) Zhang, M.; Luo, W.; Li, Z.; Yu, T.; Zou, Z. Improved photoelectrochemical responses of Si and Ti codoped α-Fe2O3 photoanode films. Appl. Phys. Lett. 2010, 97, 042105. (12) Sivula, K.; Le Formal, F.; Grätzel, M. Solar water splitting: progress using hematite (α-Fe2O3) photoelectrodes. ChemSusChem 2011, 4, 432-449. (13) Bak, A.; Choi, S. K.; Park, H. Photoelectrochemical Performances of Hematite (α-Fe2O3) Films Doped with Various Metals. Bull. Korean Chem. Soc. 2015, 36, 1487-1494. (14) Cao, D.; Luo, W.; Feng, J.; Zhao, X.; Li, Z.; Zou, Z. Cathodic shift of onset potential for water oxidation on a Ti4+ doped Fe2O3 photoanode by suppressing the back reaction. Energy Environ. Sci. 2014, 7, 752-759.

(15) Deng, J.; Zhong, J.; Pu, A.; Zhang, D.; Li, M.; Sun, X.; Lee, S.-T. Ti-doped hematite nanostructures for solar water splitting with high efficiency. J. Appl. Phys. 2012, 112, 084312. (16) Glasscock, J. A.; Barnes, P. R.; Plumb, I. C.; Savvides, N. Enhancement of photoelectrochemical hydrogen production from hematite thin films by the introduction of Ti and Si. J. Phys. Chem. C 2007, 111, 16477-16488.

(17) Li, M.; Yang, Y.; Ling, Y.; Qiu, W.; Wang, F.; Liu, T.; Song, Y.; Liu, X.-X.; Fang, P.-P.; Tong, Y. Morphology and Doping Engineering of Sn-Doped Hematite Nanowire Photoanodes. Nano Lett. 2017, 17, 2490-2495. (18) Liu, J.; Liang, C.; Zhang, H.; Tian, Z.; Zhang, S. General strategy for doping impurities (Ge, Si, Mn, Sn, Ti) in hematite nanocrystals. J. Phys. Chem. C 2012, 116, 4986-4992. (19) Kay, A.; Cesar, I.; Grätzel, M. New benchmark for water photooxidation by nanostructured α-Fe2O3 films. J. Am. Chem. Soc. 2006, 128, 15714-15721. (20) Wang, G.; Wang, H.; Ling, Y.; Tang, Y.; Yang, X.; Fitzmorris, R. C.; Wang, C.; Zhang, J. Z.; Li, Y. Hydrogen-treated TiO2 nanowire arrays for photoelectrochemical water splitting. Nano lett. 2011, 11, 3026-3033. (21) Wang, G.; Ling, Y.; Wang, H.; Yang, X.; Wang, C.; Zhang, J. Z.; Li, Y. Hydrogen-treated WO3 nanoflakes show enhanced photostability. Energy Environ. Sci. 2012, 5, 6180-6187. (22) Cooper, J. K.; Scott, S. B.; Ling, Y.; Yang, J.; Hao, S.; Li, Y.; Toma, F. M.; Stutzmann, M.; Lakshmi, K.; Sharp, I. D. Role of Hydrogen in Defining the n-Type Character of BiVO4 Photoanodes. Chem. Mater. 2016, 28, 5761-5771.

ACS Paragon Plus Environment

Page 22 of 26

Page 23 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Energy Materials

(23) Ling, Y.; Wang, G.; Reddy, J.; Wang, C.; Zhang, J. Z.; Li, Y. The influence of oxygen content on the thermal activation of hematite nanowires. Angew. Chem. 2012, 124, 4150-4155. (24) Steier, L.; Herraiz-Cardona, I.; Gimenez, S.; Fabregat-Santiago, F.; Bisquert, J.; Tilley, S. D.; Grätzel, M. Understanding the Role of Underlayers and Overlayers in Thin Film Hematite Photoanodes. Adv. Funct. Mater. 2014, 24, 7681-7688.

(25) Yang, T.-Y.; Kang, H.-Y.; Sim, U.; Lee, Y.-J.; Lee, J.-H.; Koo, B.; Nam, K. T.; Joo, Y.-C. A new hematite photoanode doping strategy for solar water splitting: oxygen vacancy generation. Phys. Chem. Chem. Phys. 2013, 15, 2117-2124.

(26) Hu, J.; Zhao, X.; Chen, W.; Su, H.; Chen, Z. Theoretical Insight Into the Mechanism of Photoelectrochemical Oxygen Evolution Reaction on BiVO4 Anode with Oxygen Vacancy. J. Phys. Chem. C 2017, 121, 18702–18709.

(27) Yao, X.; Wang, D.; Zhao, X.; Ma, S.; Bassi, P. S.; Yang, G.; Chen, W.; Chen, Z.; Sritharan, T. ScaleUp of BiVO4 Photoanode for Water Splitting in a Photoelectrochemical Cell: Issues and Challenges. Energy Technol. 2018, 6, 100 –109.

(28) Kim, T. W.; Choi, K.-S. Nanoporous BiVO4 photoanodes with dual-layer oxygen evolution catalysts for solar water splitting. Science 2014, 343, 990-994. (29) Zhao, X.; Chen, Z. Enhanced photoelectrochemical water splitting performance using morphologycontrolled BiVO4 with W doping. Beilstein J. Nanotechnol. 2017, 8, 2640-2647. (30) Wang, G.; Ling, Y.; Li, Y., Oxygen-Deficient Metal Oxide Nanostructures for Photoelectrochemical Water Oxidation and Other Applications. Nanoscale 2012, 4, 6682-6691. (31) Kim, J. H.; Jo, Y.; Kim, J. H.; Jang, J. W.; Kang, H. J.; Lee, Y. H.; Kim, D. S.; Jun, Y.; Lee, J. S. Wireless solar water splitting device with robust cobalt-catalyzed, dual-doped BiVO4 photoanode and perovskite solar cell in tandem: a dual absorber artificial leaf. ACS Nano 2015, 9, 11820-11829. (32) Chen, X.; Zhang, Z.; Chi, L.; Nair, A. K.; Shangguan, W.; Jiang, Z., Recent Advances in VisibleLight-Driven Photoelectrochemical Water Splitting: Catalyst Nanostructures and Reaction Systems. Nano-Micro Lett. 2016, 8, 1-12.

(33) Li, R.; Zhang, F.; Wang, D.; Yang, J.; Li, M.; Zhu, J.; Zhou, X.; Han, H.; Li, C. Spatial separation of photogenerated electrons and holes among {010} and {110} crystal facets of BiVO4. Nat. Commun. 2013, 4, 1432.

(34) Yang, J.; Wang, D.; Zhou, X.; Li, C. A theoretical study on the mechanism of photocatalytic oxygen evolution on BiVO4 in aqueous solution. Chem. Eur. J. 2013, 19, 1320-1326. (35) Zhao, Z.; Li, Z.; Zou, Z. Structure and energetics of low-index stoichiometric monoclinic clinobisvanite BiVO4 surfaces. RSC Adv. 2011, 1, 874-883.

ACS Paragon Plus Environment

ACS Applied Energy Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(36) Liu, W.; Zheng, W.; Jiang, Q. First-principles study of the surface energy and work function of III-V semiconductor compounds. Phys. Rev. B 2007, 75, 235322. (37) Engels, B.; Richard, P.; Schroeder, K.; Blügel, S.; Ebert, P.; Urban, K. Comparison between ab initio theory and scanning tunneling microscopy for (110) surfaces of III-V semiconductors. Phys. Rev. B 1998, 58, 7799.

(38) Wang, Y.; Feng, C.; Zhang, M.; Yang, J.; Zhang, Z. Enhanced visible light photocatalytic activity of N-doped TiO2 in relation to single-electron-trapped oxygen vacancy and doped-nitrogen. Appl. Catal. B 2010, 100, 84-90. (39) Wang, J.; Wang, Z.; Huang, B.; Ma, Y.; Liu, Y.; Qin, X.; Zhang, X.; Dai, Y. Oxygen vacancy induced band-gap narrowing and enhanced visible light photocatalytic activity of ZnO. ACS Appl. Mater. Interfaces 2012, 4, 4024-4030.

(40) Le Formal, F.; Tétreault, N.; Cornuz, M.; Moehl, T.; Grätzel, M.; Sivula, K. Passivating surface states on water splitting hematite photoanodes with alumina overlayers. Chem. Sci. 2011, 2, 737-743. (41) Pu, A.; Deng, J.; Li, M.; Gao, J.; Zhang, H.; Hao, Y.; Zhong, J.; Sun, X. Coupling Ti-doping and oxygen vacancies in hematite nanostructures for solar water oxidation with high efficiency. J. Mater. Chem. A 2014, 2, 2491-2497.

(42) Hernández, S.; Thalluri, S. M.; Sacco, A.; Bensaid, S.; Saracco, G.; Russo, N. Photo-catalytic activity of BiVO4 thin-film electrodes for solar-driven water splitting. Appl. Catal. A 2015, 504, 266-271. (43) Macdonald, J. R. Impedance spectroscopy; Wiley: New York, 1987. (44) Dotan, H.; Sivula, K.; Grätzel, M.; Rothschild, A.; Warren, S. C. Probing the photoelectrochemical properties of hematite (α-Fe2O3) electrodes using hydrogen peroxide as a hole scavenger. Energy Environ. Sci. 2011, 4, 958-964.

(45) Zhao, X.; Luo, W.; Feng, J.; Li, M.; Li, Z.; Yu, T.; Zou, Z. Quantitative Analysis and Visualized Evidence for High Charge Separation Efficiency in a Solid-Liquid Bulk Heterojunction. Adv. Energy Mater. 2014, 4, 1301785.

(46) Zhong, Y.; Li, Z.; Zhao, X.; Fang, T.; Huang, H.; Qian, Q.; Chang, X.; Wang, P.; Yan, S.; Yu, Z. Enhanced Water-Splitting Performance of Perovskite SrTaO2N Photoanode Film through Ameliorating Interparticle Charge Transport. Adv. Funct. Mater. 2016, 26, 7156-7163. (47) Costa, L. A.; Breyer, H. S.; Rubim, J. C. Surface-enhanced Raman scattering (SERS) on copper electrodes in 1-n-butyl-3-methylimidazoliun tetrafluorbarate (BMI. BF4): The adsorption of benzotriazole (BTAH). Vib. Spectrosc. 2010, 54, 103-106. (48) Wang, S.; Chen, P.; Yun, J. H.; Hu, Y.; Wang, L., An Electrochemically Treated BiVO4 Photoanode for Efficient Photoelectrochemical Water Splitting. Angew. Chem. 2017, 129, 8620-8624.

ACS Paragon Plus Environment

Page 24 of 26

Page 25 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Applied Energy Materials

ACS Paragon Plus Environment

ACS Applied Energy Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

189x135mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 26 of 26