A DFT+U Study - ACS Publications - American Chemical Society

Feb 27, 2015 - In order to understand how platinum increases efficiency, we use the density functional theory + U (DFT+U) method to model the bulk and...
6 downloads 9 Views 4MB Size
Subscriber access provided by UNIV PRINCE EDWARD ISLAND

Article 2

3

Pt-Doped #-FeO for Enhanced Water Splitting Efficiency: A DFT+U Study Ofer Neufeld, and Maytal Caspary Toroker J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/jp512002f • Publication Date (Web): 27 Feb 2015 Downloaded from http://pubs.acs.org on March 3, 2015

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 free 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 accessible to all readers and 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.

The Journal of Physical Chemistry C 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 44

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 Journal of Physical Chemistry

82x38mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

82x27mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 2 of 44

Page 3 of 44

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 Journal of Physical Chemistry

174x244mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

82x22mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 4 of 44

Page 5 of 44

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 Journal of Physical Chemistry

82x31mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

177x231mm (156 x 156 DPI)

ACS Paragon Plus Environment

Page 6 of 44

Page 7 of 44

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 Journal of Physical Chemistry

82x36mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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

177x125mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 8 of 44

Page 9 of 44

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 Journal of Physical Chemistry

191x190mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 44

Pt-Doped α-Fe2O3 for Enhanced Water Splitting Efficiency: A DFT+U Study Ofer Neufeld1 and Maytal Caspary Toroker2,* 1

The Nancy and Stephen Grand Technion Energy Program and 2Department of Materials Science and Engineering, Technion - Israel Institute of Technology, Haifa 32000, Israel Keywords: Water splitting, Density Functional Theory, DFT+U, Iron oxides, Pt-doping. Abstract Hematite (α-Fe2O3) is commonly considered for converting solar energy into hydrogen fuel through water splitting. Recent experiments performed in 2013 reached a maximum efficiency in Fe2O3 photoelectrochemical cells, while using platinum-doped Fe2O3. In order to understand how platinum increases efficiency we use density functional theory +U (DFT+U) method to model the bulk and the (0001) surface of platinum-doped Fe2O3. We also give a unique ligand field theory combined with Bader charge analysis to explain changes resulting from symmetry breaking by the dopant. First, we find that although platinum has a lower oxidation state than usual n-type dopants, platinum donates electrons. We find a theoretical ideal doping range of 0.64-2.96%at for enhanced electron conductivity, which is within the optimal range obtained by previous experiments. Second, we find that the energy gap decreases upon doping, improving solar energy absorption. Third, in agreement with previous experiments, we calculate an unfavorable increase in overpotential for oxidizing water upon platinum doping. Since platinum has both good and bad affects, we recommend bypassing this duality by platinum doping with a gradient-based strategy: high doping in the bulk for enhanced conductivity, and low doping at the surface to not interfere with catalysis. We anticipate that experimentally testing our proposed strategy will advance the development of better electrodes for photoelectrochemistry. *Corresponding author: E-mail: [email protected] , Tel.: +972 4 8294298. 1 ACS Paragon Plus Environment

Page 11 of 44

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

I.

The Journal of Physical Chemistry

Introduction

A promising route in the search for renewable energy sources is using solar energy to split water into oxygen and hydrogen.1-3 This process is environmentally “clean”, i.e., does not produce greenhouse gases. Other advantages include earth’s large water reservoir and the high energy density of hydrogen fuel.

Solar energy conversion into hydrogen can be accomplished in a photoelectrochemical (PEC) cell (see Figure 1), where solar energy is absorbed at the semiconductor photoanode. As a result of absorbing photons, excited electrons transport through the circuit and arrive at the opposing cathode where they take part in the half-cell reaction of hydrogen reduction. Holes transport in the opposite direction toward the anode surface and complete the other half-cell reaction of water oxidation. Overall, this process splits water to produce hydrogen and oxygen gas 4.

Figure 1 - Diagram for photoelectrochemical water splitting with a Fe 2O3 photoanode performing oxygen evolution (oxidation reaction) and a cathode performing hydrogen evolution (reduction reaction).𝚽 is the externally applied bias,

is the solar photon frequency, Wd is the width of the depletion layer, and Eg is the

band gap.

2 ACS Paragon Plus Environment

The Journal of Physical Chemistry

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 44

Conversion efficiencies of PEC cells are limited by the absorbing electrode's intrinsic qualities. The electrode should be a good conductor, a good solar light absorber, a good catalyst, and have an appropriate band edge alignment for water splitting. Furthermore, the electrode should be relatively inexpensive, abundant, nontoxic, and maintain stability in operating conditions. Most of these qualities can be found in hematite (α-Fe2O3, α will be dropped from this point on),4-5 and therefore this material has been widely studied for PEC cells.6-9 However, Fe2O3 also has disadvantages, including low charge mobility,10-11 high electron-hole recombination rate, and a large overpotential of 0.5-0.6[V] that is required for water oxidation.4 The efficiency of a PEC cell containing Fe2O3 can be amplified through doping.12-16 In particular, measurements show a significant increase in efficiency and photo-current in platinum (Pt)-doped cells.16-21 In these measurements, Pt increases electron conductivity and therefore was regarded as an n-type dopant. Pt also changes the electrode's morphology, causing smaller grain size, larger surface area, and a more uniform and dense Fe2O3 film, which are thought to aid in charge transport throughout the electrode.19-21 In fact, Jae Young Kim et al.20, using Pt-doped Fe2O3 with a single-crystalline “wormlike” morphology and a cobalt phosphate co-catalyst manufactured the world's highest record for a Fe2O3 PEC cell current density in 2013. All of these studies report an optimum in Pt-doping in the 0.1-4%at range, yet no study has fully explained why this is the optimal range. In addition, no study has given a complete explanation to why Pt is a successful dopant.

Therefore, the goal of this research is to use electronic structure theory to shed light on the influences Pt-doping has on the properties of Fe2O3 which are important for solar energy conversion. We go beyond previous work on first principles calculations of bulk Pt-doped

3 ACS Paragon Plus Environment

Page 13 of 44

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 Journal of Physical Chemistry

Fe2O316 by extending the analysis to account for atomic Bader charges. We give a unique assignment to atoms according to their site symmetry relative to the Pt dopant position and provide a novel explanation for routes of electron transport and for the observed ideal doping range. Furthermore, we calculate the free energies for water oxidation on Pt-doped Fe2O3(0001) and find that Pt reduces the catalytic activity of the surface.

The paper includes a computational details section with details on a Density Functional Theory (DFT) method that includes an exchange term (DFT+U) needed for calculating the electronic structure of Fe2O3. This section also includes details about calculations of the free energies required for intermediate steps in the water oxidation reaction. The results and discussion section covers our electronic structure analysis of the bulk and the (0001) surface of Pt-doped Fe2O3 and our explanation for the ideal doping concentration range. Finally, the conclusion section contains our proposed gradient doping scheme.

II.

Computational Methods

All calculations were performed using the Vienna Ab-initio Simulation Package (VASP).22-24 We chose spin polarized DFT+U formalism of Dudarev et al.25, since Fe2O3 is a first row transition metal oxide containing highly correlated 3d electrons, whose electronelectron interactions are ill described by regular exchange-correlation (XC) approximations of DFT.26-28 A U value of 4.3[eV] that was derived ab-initio was chosen for iron (Fe) atoms in order to best describe the ground state properties of Fe2O3.29-31 A Perdue-Burke-Ernzerhof (PBE)32-33 XC functional was chosen due to previous results correctly describing similar systems containing Fe2O3 and Pt metal.12,

4 ACS Paragon Plus Environment

16, 34-36

Projected

The Journal of Physical Chemistry

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 14 of 44

augmented wave (PAW) potentials37-38 were used to represent the frozen core electrons and nuclei of each atom. For all bulk calculations the 6s, 5p and the 5d orbitals of Pt, 4s and 3d of Fe and 2s and 2p of O were used as valence shells and all other inner electrons replaced by appropriate PAW potentials. In surface slab calculations the 3p Fe electrons were also used in the Fe valence shell since surface properties required the addition of Fe p orbitals and since we wanted to compare to previous work on doped Fe2O3(0001) that included explicit Fe p orbitals.12

The Kohn-Sham equations were solved with a plane wave basis set in three dimensional periodic boundary conditions to self-consistency with a tolerance of 10-5[eV] in total energy. Symmetry was not imposed to allow an accurate description of distortions within the lattice. Kspace integration was performed using the tetrahedron method with Blöchl corrections.39-40 Ionic charges were calculated using the Bader charge scheme.41 The Bader Fast Fourier Transform (FFT) grids were converged for ionic charges tolerances of 0.05e.

Our calculations involved three main structures, including pure Fe2O3 in the rhombohedral primitive cell, Pt-doped Fe2O3 in the rhombohedral 2x2x2 supercell, and a (0001) (1x1) surface slab for both pure and Pt-doped Fe2O3. For the slab structure, we calculated the free energies required for intermediate reaction steps of water oxidation. Unit cells can be found in supplementary information. The various details of each computation are presented below.

1. Pure Fe2O3: The pure Fe2O3 structure was used for comparison to Pt-doped Fe2O3. Pure Fe2O3 has a 10-atom primitive rhombohedral cell (Figure 2a). The full long range antiferromagnetic ordering of Fe2O3 was taken into account.42-43 A 7x7x7 gamma-centered k-space grid and a plane wave energy cutoff of 700[eV] converged the cell energy to a tolerance of