Yuki
Sakamoto
ab,
Yusuke
Noda
ac,
Kaoru
Ohno
ad,
Kayo
Koike
e,
Katsushi
Fujii
e,
Tomiko M.
Suzuki
f,
Takeshi
Morikawa
f and
Shinichiro
Nakamura
*a
aCluster for Science, Technology and Innovation, Nakamura Laboratory, RIKEN, 2-1, Hirosawa, Wako, Saitama 351-0198, Japan. E-mail: snakamura@riken.jp; Fax: +81-48-467-8503; Tel: +81-48-467-9477
bDepartment of Biological Information, Tokyo Institute of Technology, 4259 Nagatsuta, Midori-ku, Yokohama 226-8501, Japan
cDepartment of Materials Physics, Nagoya University, Furo, Chikusa, Nagoya, Aichi 464-8603, Japan
dDepartment of Physics, Graduate School of Engineering, Yokohama National University, 79-5, Tokiwadai, Hodogaya, Yokohama 240-8501, Japan
ePhotonics Control Technology Team, RIKEN Center for Advanced Photonics, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
fToyota Central R&D Labs., Inc., 41-1 Yokomichi, Nagakute 480-1192, Japan
First published on 15th May 2019
We report a theoretical study on iron oxyhydroxide (FeOOH). The FeOOH surface is expected to act as an efficient electrochemical catalyst for the oxygen evolution reaction (OER), because it is based on iron, an element of the fourth highest Clarke number. Experimentally, the OER activity of β-FeOOH is known to be higher than that of γ-FeOOH. However, the details of the OER mechanism and the surface reactivities of the FeOOH polymorphs have not yet been fully understood. We performed first-principles calculations of bulk and surfaces of β-FeOOH and γ-FeOOH using density functional theory, to investigate their electronic structures and catalytic activities. The calculations suggest that depending on the surface indices, several surfaces may be favored for catalytic activities.
Among these studies on 3d transition metal OER catalysts, Burke et al. reported remarkable results; iron doped nickel hydroxides showed the highest OER activity among 3d transition metal (oxy) hydroxides.5 To enhance the activity of nickel hydroxides, it is crucial to dope them with iron.6 Therefore, doping with a high concentration of iron is expected to provide a higher OER catalytic activity. However, subsequent studies have revealed that the reactivity increases only up to about 30% of iron content,7 which indicates that iron itself is the active site for the OER on these Ni-rich Ni–Fe–O systems. In contrast, when the iron content is greater than 30%, the catalytic activity decreases and accordingly iron oxyhydroxides themselves have been considered to possess a low OER activity.
Recently two of the authors, Suzuki and Morikawa, including others synthesized a highly crystalline β-FeOOH(Cl) nanorod, which is another polymorph of FeOOH8 with a much smaller size of 3 × 13 nanometers, and reported a remarkably high OER activity. According to the results of the OER activity in alkaline electrolytes, the β-FeOOH nanorod catalysts showed the highest OER activity among FeOOH polymorphs such as amorphous, α-, and γ-FeOOH.8,9 Furthermore, doping Ni into β-FeOOH highly enhanced its OER activity. Therefore, further understanding of electronic structures of bulk and surfaces of the FeOOH polymorphs will provide important clues for further improvement of the OER. The OER profiles of these catalysts must be compared in detail by various methods such as ECSA methods,10 in order to develop ideal OER electrocatalytic devices.
To the best of our knowledge, FeOOH polymorphs have been studied theoretically as well as experimentally by researchers mainly in the field of mineralogy.11–13 According to their reports, FeOOH consists of four polymorphs; α-FeOOH (goethite), β-FeOOH (akaganeite), γ-FeOOH (lepidocrocite), and ε-FeOOH. Among these polymorphs, ε-FeOOH is found under high pressure.12 These crystal structures are shown in Fig. 1. In these theoretical studies, several physical properties of FeOOH polymorphs such as magnetic configurations and electronic conduction mechanisms, have been revealed.12,14 On the other hand, theoretical investigations focusing on the catalyst, especially on their electrocatalytic activities and the origin of these differences, have been very limited.7,15
Previously, we reported a study on the band structures of various polymorphs of manganese oxides (MnO2) by density functional theory,16 as an approach to mimic natural photosynthesis, which follows up on our previous studies on the MnCaO5 cluster.17,18 We also reported a theoretical study on the electronic structure of iron doped vs. un-doped nickel hydroxides.19 In a previous study,19 we investigated the solid and surface electronic properties (obtained by periodic system calculations) in terms of molecular orbital interactions.20
As a natural extension of this study, in the present report we focus on the electronic structures of bulk and surfaces of the FeOOH catalyst, using first principles calculations. Here it is important to note that the catalytic activities are determined by a number of factors such as bulk and exposed structures, and orientation and electronic features. Amongst all, for the electrocatalysts of FeOOH polymorphs, the key factor determining its catalytic activities has not yet been understood. Therefore, the present study is aimed to obtain useful chemical insights to investigate electrocatalytic reactions. As will be explained in detail below, a series of calculations performed revealed that there are qualitative differences in the surface electronic structures between β-FeOOH and γ-FeOOH, which possess excellent and moderate OER activity among FeOOH polymorphs as described later. The result provides useful information on the relation of the OER catalytic activities with the crystal polymorphs and strategies for further improvement in the activities.
The crystal structures of the electrodes were assessed by X-ray diffraction (XRD: Rigaku, Ultima IV) using CuKα radiation at 40 kV and 40 mA. The electrochemical characteristics of the electrodes (1 × 1 cm) were investigated in 1 M KOH aqueous solution (pH 13.6) with a three-electrode configuration using an Ag/AgCl reference electrode and a Pt-wire counter electrode.
The geometries of the β and γ phases of FeOOH were prepared and optimized, owing to the references by Post et al.24 and Christensen et al.,25 respectively. For the sake of simplicity and to focus on the crystalline framework, we assumed that β-FeOOH contains neither chloride anions (Cl−) nor additional protons to neutralize the Cl− in the unit cell. As for the magnetic configuration, we assumed the antiferromagnetic configuration, following the previous studies by Alexandrov et al.14 and Guo et al.11 Notice that the 57Fe Mössbauer spectra revealed that pure β-FeOOH was paramagnetic at room temperature.8 To perform the band structure calculation, we prepared initial geometries with a magnetic primitive cell, which is the minimum cell to describe the specified magnetic configurations, based on these structures. For the k-point sampling in the geometrical optimization, we used a Monkhorst–Pack grid26 with 3 × 9 × 3 and 6 × 8 × 4 for β-FeOOH and γ-FeOOH, respectively.
To calculate the surface properties of β-FeOOH and γ-FeOOH, initial geometries were prepared using fully relaxed conventional unit cells. We take the unit cell so that the exposed surfaces were perpendicular to the z-axis (to be shown later). We used supercells including the vacuum region of at least 12 Å. The surface iron atoms were assumed to be six-coordinated and terminated by the OH− group, following previous studies.13 For structural relaxation, cell axes were fixed. Moreover, iron atoms located on the layers deeper than the first and second layers from the top surface were also fixed. For geometrical optimizations of these slab structures, we used 9 × 3 × 1 and 3 × 3 × 1 k-point grids for the (100) and the (010) surface exposed structure of β-FeOOH, whereas 3 × 4 × 1 and 8 × 2 × 1 k-point grids for the (010) and the (001) surface exposed structure of γ-FeOOH, respectively.
For the density of states (DOS) and partial DOS calculations, twice as many k-points along each axis were sampled. All symmetry k-points and symmetry k-lines in Brillouin zone were determined by the automatic flow program (AFLOW).27,28 The structures and surfaces were visualized with VESTA.29
The current–potential characteristics of the FeOOH polymorphs on CP were evaluated in 1 M KOH aqueous solutions. As shown in Fig. 2(a), the overpotential during the OER was significantly different among the FeOOH polymorphs. The overpotential of OER was the lowest on β-FeOOH, followed by γ-FeOOH and α-FeOOH. The time course of the reaction current measured at +1.60 V (vs. RHE) in Fig. 2(b) show that the OER current density of β-FeOOH was 1.3 mA cm−2, which was 7.2 times and 130 times higher than that of the γ-phase and α-phase, respectively. These results indicate that β-FeOOH exhibits the most excellent OER activity depending on its lowest overpotential among typical FeOOH polymorphs as previously reported.8
Fig. 2 (a) Current–potential characteristics and (b) time courses of currents (at +1.60 V vs. RHE) of β-FeOOH, γ-FeOOH, and α-FeOOH/CP electrodes, acquired in 1 M KOH aqueous solution. |
The calculated total and partial DOS and the band structure of β-FeOOH are shown in Fig. 4(a and b), respectively. According to the calculated band structure, there is almost no electronic interaction around the tunnel structure (plane normal to the tunnel axis; see Fig. 3(a)), since the band shows almost no large dispersion along the Γ–Y and Γ–Z lines in the reciprocal space as shown in Fig. 4(b). On the other hand, valence bands have larger dispersion along the Γ–X line (white blurred regions). The X point is at a distance half of the reciprocal vector b2 (Fig. 3(c)). These bands indicate that electronic interactions mainly exist along the tunnel structure. The bulk β-FeOOH has a band gap of 1.74 eV.
As for γ-FeOOH, there are two possibilities for its crystal system; Cmcm and Cmc21. Previously, the magnetic configuration of γ-FeOOH was reported by Alexandrov et al. using DFT calculations.14 In their study, they adopted the Cmc21 space group, which is non-centrosymmetric and has a lower symmetry than the Cmcm space group. Furthermore, they concluded that γ-FeOOH prefers the antiferromagnetic configuration, as shown in Fig. 3(b). Following their preceding study, we also assumed that γ-FeOOH has the Cmc21 space group and the antiferromagnetic configuration. The calculated total and partial DOS and band structure of γ-FeOOH are shown in Fig. 5(a and b), respectively. The valence band maximum (VBM) is located on the symmetry line connecting the E and M1 points, but not at the Γ point. The conduction band minimum (CBM) is located at the H point with the 1.18 eV band gap.
In order to check the influence of the Hubbard-U parameter for the results,30 we calculated the electronic structures of β-FeOOH and γ-FeOOH also with 3.0 eV of Hubbard-U parameter for iron. The results are shown in Fig. S2 and S3 (ESI†). The resulting electronic structures are almost the same as described above although their band gaps become slightly narrower. Therefore, we adopted the 5.0 eV for Hubbard-U value for iron in the later calculations.
The calculated total DOS and the Fe-3d and O-2p partial DOS of the (100) and (010) surface exposed structures are shown in Fig. 7(a and b), respectively. As seen from these calculated DOS, the (010) surface of β-FeOOH has a narrower band gap (0.6 eV) than the bulk structure, while the (100) surface exposed structure has almost the same band gap as that of the bulk. From these observations, the (010) surface may possibly have a higher reactivity than the (100) surface due to the narrower bandgap.
Fig. 7 Calculated total DOS (gray), the Fe-3d (blue) and O-2p (red) partial DOS of (a) for the (100) and (b) for the (010) surface exposed structures of β-FeOOH. |
In addition, we compared their surface stabilities from our calculations. Since the chemical composition ratio of these surface structures is different from the stoichiometric composition (Fe8O16H8 for the bulk structure) due to the OH surface terminations, we used the following expression to evaluate their surface stabilities ε.
The calculated surface stabilities of (100) and (010) surfaces of β-FeOOH are +0.018 and +0.197 eV Å−2, respectively. These results show that the stable (100) surface possesses a similar electronic structure to that of the bulk structure, on the other hand, the less stable (010) surface shows the electronic structure different from that of the bulk.
Fig. 9 Calculated total DOS (gray), Fe-3d (blue) and O-2p (red) partial DOS of (a) for the (010) and (b) for the (001) surface exposed structure of γ-FeOOH. |
The band gap of the (010) surface is almost identical to that of the bulk. On the other hand, in the (001) surface, the valence band crosses the Fermi energy. From these observations, the (001) surface may show a higher reactivity than the (010) surface.
We evaluated the surface stabilities in the same way as β-FeOOH surfaces. The calculated surface stabilities of (010) and (001) surfaces are +0.021 and +0.150 eV Å−2, respectively. As in the case of β-FeOOH, the less stable (001) surface showed the electronic structure different from that of the bulk structure, while the stable (010) surface showed an electronic structure similar to it.
Fig. 10 Index of iron atoms (a) for the (010) surface exposed structure of β-FeOOH and (c) for the (001) surface exposed structure of γ-FeOOH. The decomposed local DOS of Fe-3d orbitals are shown in pink for (b) the (010) surface exposed structure of β-FeOOH and (d) for the (001) surface exposed structure of γ-FeOOH. In (b and d), the numbers correspond to the layers in (a and c), respectively. In (b and d), the blue curves mean the sum of the partial DOS of Fe-3d orbitals. (They are the same as the Fe-3d DOS presented in Fig. 7(b) and 9(b).) |
According to Fig. 10(b), in the (010) surface exposed structure of β-FeOOH, the lowest conduction band is mainly assigned to the iron atoms located on layers 1 and 2, i.e., on the first and second layers from the top surface. This indicates that the (010) surface of β-FeOOH can accept electrons better than the other surfaces. On the other hand, in γ-FeOOH (Fig. 10(d)), the highest valence band, which crosses the Fermi energy, is mainly comprised of iron atoms located on the first layer, the exposed iron atoms. Thus, in both cases, the local DOS of iron atoms reveals a peculiar behavior of the surfaces, which is different from the bulk, even though all iron atoms are six-coordinated. According to the studies on the OER mechanisms under alkaline conditions, most of the proposed mechanisms contain the reactions to form intermediates, as shown below.2 On the other hand, as for the later steps, several O2 forming mechanisms are proposed.
M + OH− → MOH |
MOH + OH− → MO + H2O(l) |
Focusing on the (010) surface of β-FeOOH, we show the spatial distribution of the Kohn–Sham wave functions at the Γ point in real space. In Fig. 11, the up-spin pseudo wave functions of the highest occupied level and the lowest unoccupied level of the β-FeOOH(010) surface exposed structure are visualized analogous to the HOMO/LUMO visualization in molecular quantum chemistry. As shown in Fig. 11(a), the highest occupied orbital is widely distributed on the oxygen atoms. On the other hand, the lowest unoccupied orbital is localized on the up-spin polarized iron atoms located on the surface. The down-spin pseudo wave functions have a similar behavior due to the antiferromagnetic configuration. The distributions of these pseudo wave functions are consistent with the behavior described by the decomposed local DOS shown in Fig. 10(b), indicating again that the (010) of β-FeOOH is favorable for the OER activity. The DOS and pseudo wave function of the lowest conduction bands of the β-FeOOH(010) surface exposed structure suggest that the accepter site is localized on the atoms near the surface. On the other hand, the highest valence band distributes widely in the whole structure, not localized on specific location. In view of the fact that a chemical reaction is an event occurring at the localized sites, the (010) surface may prefer more electron accepting reactions than electron donating reactions.
In the previous DFT study, Otte et al.13 compared the most common surface of α-FeOOH, β-FeOOH and γ-FeOOH with various termination structures. They concluded that the electronic structures of surface iron atoms mainly depend on their coordination numbers and termination structures, but not on the crystalline phase. In the present study also, the calculated electronic structures of the same surface: the (100) surface of β-FeOOH and the (010) surface of γ-FeOOH, are consistent in their results qualitatively. On the other hand, our calculations show that some surfaces have distinct properties depending on the surface indices even in the same crystal structures. It is natural to infer that the catalytic activity depends on the surface properties. Therefore, the stability of the active surface is a critical problem for the catalysts. In general, the most exposed surface of crystals is relatively more stable than the others. On the other hand, the catalytic activity may not always belong to the most stable surface.
As seen above, the (010) surface of β-FeOOH have lower conduction bands than the bulk, while the (001) surface of γ-FeOOH has higher valence bands than the bulk. These features are contrasting with each other. Since the highest occupied states are localized near the (001) surface in γ-FeOOH, it is expected that the electrolytes existing nearby can easily access the surface (of the electronic structure) and may cause interaction. On the other hand, since the lowest unoccupied states are localized on the (010) surface of β-FeOOH, the surface electronic structure is expected to be suitable for the initial OER steps. This observation can be one explanation of the difference in the catalytic activities and stabilities of the β-FeOOH and γ-FeOOH.
Footnote |
† Electronic supplementary information (ESI) available. Relaxed coordinates of all the calculated structures are presented. See DOI: 10.1039/c9cp00157c |
This journal is © the Owner Societies 2019 |