Tailoring Ca2Mn2O5 Based Perovskites for Improved Oxygen

Oct 29, 2018 - Advanced Quantum Theory: Molecules to Materials Group, School of ... evolution reactions (OERs) and oxygen reduction reactions (ORRs)...
0 downloads 0 Views 1MB Size
Subscriber access provided by The University of Texas at El Paso (UTEP)

Article 2

2

5

Tailoring CaMnO Based Perovskites for Improved Oxygen Evolution Reaction Neha Bothra, Sandhya Rai, and Swapan K. Pati ACS Appl. Energy Mater., Just Accepted Manuscript • DOI: 10.1021/acsaem.8b01301 • Publication Date (Web): 29 Oct 2018 Downloaded from http://pubs.acs.org on October 31, 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 24 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

Tailoring Ca2Mn2O5 Based Perovskites for Improved Oxygen Evolution Reaction Neha Bothra,†,‡ Sandhya Rai,†,‡ and Swapan K. Pati∗,† †Advanced Quantum Theory: Molecules to Materials Group, School of Advanced Materials (SAMat), Theoretical Science Unit, Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur, Bengaluru, 560064, India. ‡Authors contributed equally E-mail: *[email protected] Abstract The commercialization of numerous electrochemical energy systems is inhibited by the sluggish kinetics of oxygen evolution reactions (OER) and oxygen reduction reactions (ORR). OER catalysts should be stable, active and cheap for their application in electrochemical reactions. Recently, Mn based double perovskites have emerged as suitable catalysts for the OER. We present a density functional theory based comprehensive study on the pristine and doped Ca2 Mn2 O5 , where Ca is replaced with inexpensive alkaline earths, main group elements and lanthanides to elucidate how the electronic structure can aptly be tuned by doping, for the best catalytic performance. The catalytic performance has been assessed based on the theoretical overpotential calculated for each case. Doping assisted modification of the eg filling of the metal centre is observed which affects the catalytic activity. 30% Ce doped Ca2 Mn2 O5 shows the best activity in terms of the lowest overpotential value of 0.14 V. We have developed a simple predictive model using 4 electronic structure based descriptors. This study suggests Ce0.7 Ca1.3 Mn2 O5 as an efficient, cost effective catalyst for OER.

1

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

Abbreviations OER,DFT

Keywords Catalysis, Oxygen Evolution Reaction (OER), Density Functional Theory (DFT), overpotential, electrochemistry

Introduction Increasing energy demands and growing concerns about global warming have led to a surge in the discovery of economical and efficient catalysts that can be used for the electrochemical energy conversion processes, like oxygen evolution reaction (OER), oxygen reduction reaction (ORR) and hydrogen evolution reaction (HER). 1–6 However, the sluggish kinetics of the OER process, demands for high overpotential compared to the other two, inhibiting its extensive use in generating green energies. This is also responsible for making it one of the most discussed and explored research areas. Till date, rutile-type RuO2 is considered to be the most efficient electrocatalyst for OER used on industrial scale, followed by IrO2 in acidic medium. 7 But, low abundance and high cost of the metals along with the tendency to get further oxidized ceases its usage at industrial scale. 8–12 Another goal of catalyst design is to propose catalysts that are sufficiently economical from the synthesis point of view. In this respect, various perovskites have emerged as a promising substitute to the precious metals. 3,4,13 The catalytic activity of 3d transition metal based perovskites with varying doping percentage was identified by Bockris and Otagawa, whose work pioneered the way in understanding the reaction pathway and also the property-activity relationships. 14 These perovskites also serve with an extra benefit in terms of the compositional flexibility provided by them, giving a platform to tune the ac-

2

ACS Paragon Plus Environment

Page 2 of 24

Page 3 of 24 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

tivity via its electronic structure. In an ABO3 type perovskite, both ‘A’ and ‘B’ sites can be systematically varied with metals of different radii and valences. Partial doping is also a very common practice, bringing a systematic change in the physiochemical properties of the material. There is plethora of literature identifying different combinations of elements, 13,15–18 but it has been observed that the oxygen-deficiency enhances the catalytic activity incurred by the increased electrical conductivity, leading to the development of a series of double perovskites. 19 In addition to this, porosity in the nanostructured materials leads to defects and variable oxygen stoichiometry, resulting in enhanced electrocatalytic performance for ORR/OER reactions. 20–22 Plethora of literature is available on Co based perovskites that have shown extremely good catalytic activity for both ORR and OER processes. 1–4,23,24 A recent work by Ashok et al. on oxygen reduction reactions, reveals that in terms of stability, LaCoO3 lies above LaMnO3 , making the former a better catalyst. However, Mn based perovskite gives a higher current density than its Co counter-part, suggesting its potential to show better performance. 25 Hence, it is desirable to design Mn based perovskites with enhanced stability. It has been found that double perovskites have the potential in bridging this gap. An Mn based double perovskite was found to be stable and shows good catalytic activity towards OER/ORR. 22 The superior activity is a result of inherent surface defects and variable oxygen stoichiometry. 22 In addition to that, these double perovskites result in the 5 coordinated square pyramidal Mn centre, where as, a conventional perovskite has a 6 coordinated octahedral symmetry. This results in the further splitting of eg levels which were degenerate earlier. This eg filling is a well established descriptor for the activity of a catalyst. 4,26 The eg splitting causes a more controlled change over the occupancy which ideally is supposed to be ∼1.2 for a good OER catalyst. Though, Ca2 Mn2 O5 shows good catalytic activity, it still does not occupy the top position of the traditional volcano plot, which is widely used to scrutinize the efficiency of catalysts. With the idea of enhancing the catalytic activity of the pristine Ca2 Mn2 O5 , we have done a comprehensive theoretical modeling of Ca2 Mn2 O5 by replacing Ca with alkaline earth,

3

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

lanthanides and main group elements in accordance with maintaining the Goldschmidt tolerance factor (t) between 0.71 - 0.90. Goldschmidt tolerance factor is an indicator for the stability and distortion of crystal structures. It can be used to calculate the compatibility of an ion with the crystal structure. This factor helps in choosing the substituent dopants keeping the overall geometry of the perovskite unchanged. In addition to that, it is already known that for A2 B2 O5 , the valency of A can vary between +2 to +4. Thus, we have chosen the dopants of different valency in order to see its impact on the catalytic activity of the surface. The doping is sequentially done by varying the percentage of dopant from 10% to 40%, which eventually tunes the electronic structure of the catalyst. We have also developed a linear regression based model that can predict the activity of the catalyst from its pristine structure. We hope that this work will widen the horizons of rational catalyst designing.

Methods First principles calculations were carried out using Quantum Espresso 5.2.0 package. 27,28 Generalized gradient approximation was incorporated via the exchange correlation functional parameterized by Perdew-Burke-Ernzerhoff (PBE). 29,30 The core electrons were treated in the Rappe-Rabe-Kaxiras-Joannopoulos ultrasoft pseudopotential method. 31 Due to strong correlation of d–electrons in Mn, the Hubbard U correction 32 was applied to Mn with a specific value of 4.0 eV, in accordance with some recent calculations on Mn. 33 Spin polarization was applied on Mn with spin structure as ferromagnetic, as is reported about the spin state for experimentally synthesized Ca2 Mn2 O5 . 22 The 3×3×1 Monkhorst Pack 34 reciprocal grid was used to sample out the total energy numerically. In order to accelerate the computation near the Fermi level, a Methfessel-Paxton smearing of 0.01 was used. 35 The Broyden–Fletcher–Goldfarb–Shanno algorithm 36 was used for full relaxation of the atomic positions of both, slab and adsorbed intermediates. Conjugate gradient method was used for the electronic relaxation and simultaneous DOS, PDOS and Bader charge calculations. 37

4

ACS Paragon Plus Environment

Page 4 of 24

Page 5 of 24 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 wave function kinetic energy cutoff was 30 Ry, for energy accuracy. To avoid the asymmetric slab induced dipole moments, all the slabs were constructed symmetrically considering the (001) plane A–Mn–O (A=Ca/Dopants) as the terminating layer. As an initial guess, the adsorbates were placed over the Ca and Mn centres. But on structural relaxation, the system always preferred the Mn site for binding. There were two Mn sites on the surface, and we have chosen the one giving lower value of adsorption energy corresponding to OH binding as the catalytically active Mn to proceed with the further calculations in all the cases. All the adsorbates, ∗ O, ∗ OO, ∗ OH and ∗ OOH were bonded to the catalytically active Mn. This is reasonable assumption for modeling OER, which has been used in recent literature. 38 In order to calculate the energy of the ions in solvent, computational hydrogen electrode 39 (CHE) developed by Nørskov et al. was used as the reference electrode. Thus, all the electrochemical reaction energies were mapped to the value with reference to reversible hydrogen electrode (RHE). Solvation is expected to influence the adsorption of intermediates on the surface. However, computation of the solvation free energy corrections requires further validation of the models itself with the experiments, which is demanding. In addition to that, little data is available on the explicit water models on perovskites. Considering the results on Pt(111) surface, it was found that the *OH adsorption is stabilized by 0.50 eV and that of *OOH by 0.25 eV, on inclusion of solvation effects. 40 The net effect of the additional water layer on the overpotential will be 0.25 V. 40 Even if the solvation effect is applied, the qualitative interpretation of the results remains unchanged, in view of the fact that the rate determining step remains the same for all the catalysts considered. 41 The thermodynamic and entropic aspect was included in the 0K electronic energy of the slab+ adsorbate configuration from the previous works of Montoya et al. 41 It should be noted that for our calculations, all the reactions are considered at U=0 and pH=0. For all the calculations on (un)doped Ca2 Mn2 O5 , a 5–atomic layer slab with (001) plane A–O–Mn exposed layer was constructed with 12Å vacuum layer. The (001) surface with transition metal ions has been widely used as a catalytically active and stable surface for

5

ACS Paragon Plus Environment

ACS Applied Energy Materials

simple perovskites in DFT calculations. 42,43 We hypothesize this to be the case of Ca2 Mn2 O5 . The simulation slabs used for DFT calculations were fully oxidized and the structures were relaxed. This is inline with other calculations considering the fact that the OER process occurs in KOH where these slabs would expected to be oxidized.

Results and discussion Electronic structure of the catalyst Ca2 Mn2 O5 has an orthorhombic crystal existing in high spin state, where Mn shows a 5 coordinated square pyramidal geometry with O atoms as its neighbours. 22,44 This is a striking 1.5 Mn(dz2) Mn(dx2-y2)

1.0

O(px) O(pz)

0.5 PDOS

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 6 of 24

0.0

-0.5

-1.0

-1.5 -10.0

-5.0

0.0 E-EF (eV)

5.0

10.0

Figure 1: The partial density of states (PDOS) of the pristine Ca2 Mn2 O5 for the Mn(dx2 −y2 ), Mn(dz2 ), O(px ) and O(pz ), where Mn is the catalytically active site. difference from the orthorhombic CaMnO3 , where Mn centre has an octahedral symmetry and is surrounded by 6 O atoms. It is probably due to this symmetry breaking, that the eg levels become non-degenerate and give rise to more catalytically active Mn centres. The five coordinated Mn can exist in two geometries, i.e. trigonal bipyramidal and the square pyramidal and the d-orbital splitting is strikingly different in the two cases (Figure S1 of SI). The two geometries can interchange which may lead to a significant rearrangement of the 6

ACS Paragon Plus Environment

Page 7 of 24 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

d-orbitals. 45 This interchange is mainly governed by the electronic character of the ligands around, the d-shell configuration and the covalency of the metal-ligand bond. 45 The PDOS of Ca2 Mn2 O5 reveals the splitting of the eg levels (Figure 1). This was also observed in some of the earlier works on cobalt based perovskites. 24,46

Selecting the dopants The above discussion leads to the fact that tuning the eg level of the transition metal centre may lead to good catalytic activity. 3,4,13 In this work, eg levels are tuned by replacing Ca

Surface

Bulk

(a)

(b)

Figure 2: (a) Relaxed structure of Ca2 Mn2 O5 (001) surface with 5 layers. The bottom three layers were constrained to represent the bulk, and the top two layers were allowed to change their positions in the presence of adsorbate. The colour code: limegreen=Mn, Cyan=Ca, Red=O (b) The highlighted Ca are the positions that were considered for sequential doping. (Site ‘A’ in A2 B2 O5 ) on the top two layers of the (001) plane of pristine Ca2 Mn2 O5 with Sr, In, Bi, La, Ce and Eu (Figure 2). The doping percentage is sequentially varied from 10 to 40%. Positions of the dopant were varied for all the cases and the minimum energy confirmation was chosen for further study. Similarly, active site of the catalyst was analyzed on the basis of adsorption energies with respect to OH to find the minimum energy catalyst– adsorbate complex. The details of these calculations are provided in the methods section of 7

ACS Paragon Plus Environment

ACS Applied Energy Materials

10% Ce Doped 1.5

1.5

0.5 0.0

0.5 0.0

-0.5

-0.5

-1.0

-1.0

-1.5 -10.0

-5.0

0.0 E-EF (eV)

5.0

Mn(dz2) Mn(dx2-y2) O(px) O(pz)

1.0 PDOS

PDOS

20% Ce Doped

Mn(dz2) Mn(dx2-y2) O(px) O(pz)

1.0

-1.5 -10.0

10.0

-5.0

30% Ce Doped 1.5

0.5 0.0

-0.5

-0.5

-1.0

-1.0 -5.0

0.0 E-EF (eV)

10.0

5.0

Mn(dz2) Mn(dx2-y2) O(px) O(pz)

1.0

0.0

-1.5 -10.0

5.0

40% Ce Doped

PDOS

0.5

0.0 E-EF (eV)

1.5

Mn(dz2) Mn(dx2-y2) O(px) O(pz)

1.0 PDOS

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 8 of 24

-1.5 -10.0

10.0

-5.0

0.0 E-EF (eV)

5.0

10.0

Figure 3: The partial density of states (PDOS) of the Ce doped Ca2 Mn2 O5 of the Mn(dx2 −y2 ), Mn(dz2 ), O(px ) and O(pz ), where Mn is the catalytically active site. supporting information (SI). It is known that due to the difference in radii and electronegativity of site ‘A’, the d-filling of the corresponding ‘B’ site gets altered due to modification in the B–O bond. 46 Change in electronic environment due to doping is analyzed in terms of the projected density of states (PDOS) plots. As is evident from Figure 3, varying the doping percentage of ‘A’ modifies the O-p and metal d-band. The PDOS of pristine catalyst (Figure 1) shows the presence of dz 2 orbital crossing the Fermi level which is in close proximity to the occupied O-pz orbital. This observation is in line with an earlier reported result by Lu et al. 44 As a result of doping (Ca with 10% of another element), the density of states corresponding to Mn-dz2 are shifted and reduced and that of O-pz is increased. As doping percentage is increased, the Mn-dz2 , dxz and dyz tend to get more stabilized and shift away from the Fermi level. On 30% doping with Ce, Mn-dx2 −y2 and O-pz and O-px become 8

ACS Paragon Plus Environment

Page 9 of 24 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

degenerate and lie near the Fermi level. On further increase in doping, the density of states is shifted away from the Fermi level. It is important to note that, this reduction in catalytic activity on doping beyond 30% may occur at any composition between 30% and 40%. Our modelled catalysts are restricted by a resolution of 10% with respect to each doping concentration. This suggests that doping ‘A’ has a huge implication on the electronic structure around ‘B’ and hence can influence the catalytic activity considerably.

Free energy calculations As proposed by Goodenough et al., the OER process on the perovskite surface takes place in alkaline medium. 3 The reaction scheme proposed by him is given in Figure 4, and involves 4 electron transfer steps. As per recent studies, there are chances of the lattice oxygen to take part in the reaction, making the OH− and e− transfer non-concerted, as a result of which reaction becomes pH dependent. 24 This is thermodynamically favoured when the O 2p states at the perovskite Fermi level lie above the redox energy of the O2 /H2 O couple. As per our analysis of the PDOS plots, this does not happen in any of the cases, thus ruling out the possibility of a lattice oxygen participation in the reaction (See SI for details). In order to model the thermodynamics of the reaction, computational hydrogen electrode model proposed by Nørskov et al. was used. 39 The reactions were modeled considering the OH− and e− transfer as a concerted step. 47–49 The details are available in the SI. The linear free energy relationships (LFER) allow us to relate the thermodynamics to the kinetics of the reaction for an isoentropic process via Brønsted–Evans–Polyani relation. Thus, a thermodynamic picture is deemed accurate to evaluate the performance of the catalysts. 50,51 It is already known that the free energies scale up/down with the change in applied external potential (U). Hence, the potential determining step remains unchanged, making the theoretical overpotential independent of U and we represent this potential at 0 pH as Uo . Figure 5 shows the standard equilibrium potential (Uox , x=1,2,3,4) for all the four steps. For an ideal catalyst this value 9

ACS Paragon Plus Environment

ACS Applied Energy Materials

OH

2

O

-

e

OH -

+

4

-

O H2

OHn+ O-M-O 1

OO2(n+1)+ O-M-O

OER

O2(n+1)+ O-M-O

S

-

OH

3

RD + -

e

OOHn+ O-M-O

RD S

2

O H2

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 24

-

e

OH -

Figure 4: The reaction scheme as proposed by Goodenough et al. 3 for an OER process taking place on a perovskite surface in alkaline medium. should be the same for all the steps i.e. 1.23 V, at zero U. However, in real catalysts this does not happen. We find U2o , which involves the O–O bond formation, as the potential determining step for all cases. This further determines the overpotential of the reaction via Equation 1. 52 UOER = max[Uo1 , Uo2 , Uo3 , Uo4 ];

η OER = UOER − 1.23 V

(1)

It is interesting to note that the η predicted for Ca2 Mn2 O5 using the above formulation gives a value of 0.26 V which is close to the experimentally reported value of 0.27 V, providing a good confidence in our theoretical predictions. The rate determining step (RDS) is controlled by the eg filling of the underlying catalyst metal centre (Mn, in our case). It is either the O– O bond formation or the proton abstraction from the oxy-hydroxide group, which determine the rate OER. 3 The O–O bond formation is the RDS when eg > 1, which is consistent with our calculations. 3,53 It is evident from Figure 5 that increasing the doping concentration does not always reduce the Uo . Rather, in all the cases, Uo increases in going from 30% to 40%. In most of the cases, 30% doping concentration gives minimum overpotential value. This

10

ACS Paragon Plus Environment

Page 11 of 24

2.00

o

U1

Uo2 Uo3 Uo4

Equilibrium Potential (V)

1.50

1.00

0.50

0.00 %

0%

40

e4

C

La

%

0%

40

i4

B

In

%

0%

30

e3

C

La

%

0%

30

i3

B

In

%

%

20

30

Sr

Eu

%

0%

20

e2

C

La

%

0%

20

i2

B

In

%

%

10

20

Sr

Eu

%

e

0%

10

e1

C

La

%

0%

10

i1

B

In %

tin

10

is

Sr

Pr

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 5: A comparison of all the Uox (x=1,2,3,4) for all the steps of the OER process in all the studied (un)doped systems. The red line corresponds to the thermodynamic equilibrium potential of 1.23 V, for an ideally reversible process. suggests that electronic modification caused at this percentage is ideal for OER. It is known that the nature of ‘A’ can tune the concentration of oxygen vacancies and the oxidation state of the central metal, which in turn govern the overall OER activity. 46 As discussed earlier, PDOS for Ce (Figure 3) shows rearrangement of t2g and eg levels as the doping percentage is varied. At 30%, the dx2 −y2 overlaps with the p-band of lattice oxygen. This allows for more covalency in the metal-oxygen bonding, thereby modifying the eg filling of the metal. It can be hypothesized from here that the electronic modulation at 30% doping of Ce is perfectly ideal for OER catalysis, thereby resulting in lower overpotential values. The eg filling also depends on the oxidation state of the metal centre. It has been found earlier that the oxidation state of metal centre varies with the radius of ‘A’. 46 Grimaud et al. have reported an increase in the oxidation state of the metal, with respect to the increase in the ionic radius of ‘A’. Bader charge analysis shows a reduction in charge of Mn which is inversely proportional to the radii of dopant ‘A’. Interestingly, it is also observed that the electronegativity of ‘A’ also effects the oxidation state of ‘B’. Bader charge analysis suggests that as the doping percentage is increased, the oxidation state of the metal centre tends to 11

ACS Paragon Plus Environment

ACS Applied Energy Materials

decrease (details in the SI).This points towards an increase in the eg filling that exceeds more than 1 (ideal case for Mn3+ with weak ligands). From our estimation we conjecture that both, over oxidation or over reduction of Mn leads to poor overpotential and we identify 30% doping as ideal for these Mn based perovskites.

Scaling relations The adsorption energies of all the intermediates suggests that the bond strengths follow the order: E∗ OH > EO∗ > E∗ OOH > EOO∗ . It is evident from Figure 6 that all the adsorp4.00

∆GO* 3.50 ∆G*OOH 3.00 ∆GO*, ∆G*OOH (eV)

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 24

2.50 2.00 1.50 1.00 0.50 0.00 -0.50 -2.00

-1.50

-1.00 -0.50 ∆G*OH (eV)

0.00

0.50

Figure 6: Scaling relations between adsorption free energy values of ∗ O and ∗ OOH with that of ∗ OH for various (un)doped catalysts. ∆G∗ O = 1.70∆G∗ OH + 1.97, ∆G∗ OOH = 0.93∆G∗ OH + 3.18. tion free energies of intermediates scale with ∆G∗ OH which is in line with earlier reported results. 50–52,54,55 Figure 6 shows that ∆G∗ OH strongly correlates with ∆G∗ OOH and ∆G∗ O giving a Pearson correlation 56 coefficient of 0.92 and 0.90, respectively. These scaling relations consequently affect the catalyst design, suggesting that even the most active catalyst will have a non-zero overpotential. The binding energy of ∗ OH and ∗ OOH have 1:1 scaling with a constant difference of 3.18 eV which is in agreement with the earlier finding by Rossmeisl et al. where it is mentioned that this scaling is independent of the surface. 48 However, this indicates towards a non-zero overpotential because ideally this difference should be 2.46 12

ACS Paragon Plus Environment

Page 13 of 24 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

ron St

g

We ak bin din g

eg gl din n i b

leg

Figure 7: Activity volcano relation for (un)doped Mn based double perovskites. η OER = 0.069(∆G∗ O − ∆G∗ OH ) + 0.086 (2*1.23) V. The ∆G∗ O and ∆G∗ OH have a 2:1 scaling owing to the double bond character of the metal–∗ O bond. A constant intercept also implies that these intermediates prefer same type of binding site (Mn centre in this case).

Identifying Activity Descriptors Existence of a universal scaling relation between ∆G∗ OOH and ∆G∗ OH , eliminates the possibility of a variation in η OER with it. Earlier reports state that the adsorption of ∗ O on the catalyst surface determines the η OER of the reaction, suggesting step (2) or step (3) to be the rate determining step (RDS). Thus, ∆G∗ O –∆G∗ OH is identified as a universal descriptor for the OER process. 51,52 In agreement with these findings, Figure 7 shows a linear correlation of ∆G∗ O –∆G∗ OH with η OER with a Pearson correlation of 0.93. The difference in the adsorption energies of ∗ O and ∗ OH is thus determining the overpotential of the reaction, which for an ideal catalyst should be zero. The colinearity in the plot suggests a similar RDS for all the catalysts considered in this study. This is known as the universal volcano relation, where in our case, the adsorbate binds weakly to the catalyst. The higher the catalyst lies towards the peak, the better is its catalytic performance. From Figure 7, we identify 30% Ce doped Ca2 Mn2 O5 (Ce0.7 Ca1.3 Mn2 O5 ) for the class of catalysts considered in this study. η OER for this system is 0.14 V, which is less than that for pristine Ca2 Mn2 O5 (η OER =0.26 V). This difference in the η OER can be attributed to the stability of ∗ O on the Ce doped sur-

13

ACS Paragon Plus Environment

ACS Applied Energy Materials

Test set Train set

0.28 0.26 0.24 Calculated

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 24

0.22 0.20 0.18 0.16 0.14 0.14

0.16

0.18

0.20 0.22 Predicted

0.24

0.26

0.28

Figure 8: The predicted versus the calculated values of η OER using the linear regression based model. The R2 for this model is 0.83. face owing to the oxophilicity of the Ce atom. Presence of Ce in vicinity of the reactive site facilitates in the O-O bond formation via a non covalent interaction of Ce with O∗ resulting in a reduction of electron density over O, facilitating the bond formation with nucleophilic OH− . The overpotential estimated for Ce doped Ca2 Mn2 O5 is even lower than the experimentally reported overpotential for some Ni-based hydroxides hierarchical nanoarrays, which is 0.23 V, and is reportedly the best catalyst for OER till date, to the best of our knowledge. 57

Predictive model for OER activity Development of quantitative structure–activity relationship is central to the rational design of novel catalysts. These models enable to develop relationship between the desired catalytic activity and the structural/electronic properties of the material. This enables a fast screening of the probable candidates for catalysis and saves a lot of time, effort and expenses. Such models have been developed in the past which used around 14 structural descriptors on a data set for 51 perovskite catalysts. 58 Here, we use a set of 4 descriptors on a data-set of 21 catalysts where only the ‘A’ site is varied. The descriptors are based on the geometry, like the tolerance factor, percentage doping, and on the electronic structure, like, Fermi energy, pband centre for the surface lattice oxygen and metal-oxygen covalency. These descriptors are

14

ACS Paragon Plus Environment

Page 15 of 24 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

easy to estimate from the very basic theoretical calculations and are highly intuitive about the nature of the material under study. Figure 8 shows the predicted versus the calculated values of η OER . The corresponding model is given in Table 1. It is evident from here, that Table 1: Details of the multiple regression based predictive model. The intercept value for the model is 0.2000 eV. The importance value is estimated as the difference in the corresponding R2 values, when the corresponding variable is replaced with its average value. See SI for details. Variable Coefficient Importance Fermi Level (eV) -0.0075 0.10 p–band centre (eV) 0.0100 0.15 t×%doping 0.0393 (eV) 0.10 Covalency 0.1395 (eV) 0.17 this simple model is able to predict the η OER with a reasonably good accuracy (R2 =0.86). In order to develop the model, data are randomly splitted into training set (composing 70% of the total) and test set. It may be argued that with such a small data-set and low variation in the test systems, this model may be highly biased. However, we would like to emphasize on the fact that a relatively simple model is able to capture a wide range of elements of the periodic table and can stand useful for further generalization by testing it on a wider range of data. Moreover, the model is relatively simple and thus, much more intuitive for rational materials design. An importance analysis is done in order to evaluate the contribution of the variables towards the model. It is found that the covalency of the metal-oxygen bond is the most dominant factor in determining the OER activity, followed by p-band centre. The Fermi energy and the tolerence factor parameter have equal effect on the predictive ability of the model.

Conclusions In summary, a comprehensive study has been made on Mn based oxygen deficient double perovskite i.e. Ca2 Mn2 O5 by doping Ca with Sr, In, Bi, La, Ce and Eu. The dopant 15

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

percentage is varied from 10% to 40%, in order to see the effect of these different percentages on the catalytic activity of the perovskite. It is found that due to doping, the density of states near the Fermi level are affected, leading to a square pyramidal type d-orbital ordering. This leads to the tuning of the eg level and hence the catalytic activity of the system. We identify 30% doping as ideal, which tunes the electronic structure optimally for good OER activity. For all the cases, we identify the O–O bond formation as the rate determining step, and it correlates well with being the RDS for eg occupancy > 1. Of all the elements studied, 30% Ce doped catalyst shows the best results in terms of the lowest thermodynamic overpotential. We develop a linear regression based predictive model that uses simple, but intuitive descriptors, like metal-oxygen bond covalency, p–band centre, Fermi level and tolerance factor. These kind of models would make the catalyst design faster and efficient by considerably reducing the effort and computational cost that goes into the screening of the potential catalysts. With these results, we hope for the development of novel cost-effective perovskite based catalysts for OER.

Acknowledgement NB thanks CSIR for SRF. SR thanks DST-Nano mission for Post doctoral Fellowship. SKP thanks J. C. Bose National Fellowship and SERB/DST for funding.

Supporting Information Available The following files are available free of charge. • MSPerovskiteSI.pdf: Contains methodology in detail, MLR model, PDOS plots for pristine catalysts, standard equilibrium potential for all the catalysts, structures of the intermediates sitting on the catalyst surface.

16

ACS Paragon Plus Environment

Page 16 of 24

Page 17 of 24 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

References 1. Hwang, J.; Rao, R. R.; Giordano, L.; Katayama, Y.; Yu, Y.; Shao-Horn, Y. Perovskites in Catalysis and Electrocatalysis. Science 2017, 358, 751–756. 2. Han, B.; Stoerzinger, K. A.; Tileli, V.; Gamalski, A. D.; Stach, E. A.; Shao-Horn, Y. Nanoscale Structural Oscillations in Perovskite Oxides Induced by Oxygen Evolution. Nat Mater. 2017, 16, 121–127. 3. Suntivich, J.; May, K. J.; Gasteiger, H. A.; Goodenough, J. B.; Shao-Horn, Y. A Perovskite Oxide Optimized for Oxygen Evolution Catalysis from Molecular Orbital Principles. Science 2011, 334, 1383–1385. 4. Suntivich, J.; Gasteiger, H. A.; Yabuuchi, N.; Nakanishi, H.; Goodenough, J. B.; ShaoHorn, Y. Design Principles For Oxygen-Reduction Activity On Perovskite Oxide Catalysts For Fuel Cells And Metal–Air Batteries. Nat. Chem. 2011, 3, 546–550. 5. Muy, S.; Bachman, J. C.; Giordano, L.; Chang, H.-H.; Abernathy, D. L.; Bansal, D.; Delaire, O.; Hori, S.; Kanno, R.; Maglia, F. Tuning Mobility And Stability Of Lithium Ion Conductors Based On Lattice Dynamics. Energy Environ. Sci. 2018, 11, 850–859. 6. Ahmad, E. A.; Tileli, V.; Kramer, D.; Mallia, G.; Stoerzinger, K. A.; Shao-Horn, Y.; Kucernak, A. R.; Harrison, N. M. Optimizing Oxygen Reduction Catalyst Morphologies from First Principles. J. Phys. Chem C 2015, 119, 16804–16810. 7. Kim, N.-I.; Sa, Y. J.; Yoo, T. S.; Choi, S. R.; Afzal, R. A.; Choi, T.; Seo, Y.-S.; Lee, K.S.; Hwang, J. Y.; Choi, W. S. Oxygen-Deficient Triple Perovskites As Highly Active And Durable Bifunctional Electrocatalysts For Oxygen Electrode Reactions. Sci. Adv. 2018, 4, 9360–9367. 8. Wang, H.; Lee, H.-W.; Deng, Y.; Lu, Z.; Hsu, P.-C.; Liu, Y.; Lin, D.; Cui, Y. Bifunc-

17

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

tional Non-Noble Metal Oxide Nanoparticle Electrocatalysts Through Lithium-Induced Conversion For Overall Water Splitting. Nat. Commun. 2015, 6, 7261–7264. 9. Pfeifer, V.; Jones, T. E.; Wrabetz, S.; Massué, C.; Vélez, J. J. V.; Arrigo, R.; Scherzer, M.; Piccinin, S.; Hävecker, M.; Knop-Gericke, A. Reactive Oxygen Species In Iridium-Based Oer Catalysts. Chem. Sci. Trans. 2016, 7, 6791–6795. 10. Ng, J. W. D.; García-Melchor, M.; Bajdich, M.; Chakthranont, P.; Kirk, C.; Vojvodic, A.; Jaramillo, T. F. Gold-Supported Cerium-Doped NioX Catalysts For Water Oxidation. Nat. Energy. 2016, 1, 16053–16058. 11. Lu, Y.-C.; Xu, Z.; Gasteiger, H. A.; Chen, S.; Hamad-Schifferli, K.; Shao-Horn, Y. Platinum- Gold Nanoparticles:

A Highly Active Bifunctional Electrocatalyst For

Rechargeable Lithium- Air Batteries. J. Am. Chem. Soc. 2010, 132, 12170–12171. 12. Park, S.; Shao, Y.; Liu, J.; Wang, Y. Oxygen electrocatalysts for water electrolyzers and reversible fuel cells: status and perspective. Energy Environ. Sci. 2012, 5, 9331–9344. 13. Suntivich, J.; Gasteiger, H. A.; Yabuuchi, N.; Shao-Horn, Y. Electrocatalytic Measurement Methodology Of Oxide Catalysts Using A Thin-Film Rotating Disk Electrode. J. Electrochem. Soc. 2010, 157, B1263–B1268. 14. Bockris, J. O.; Otagawa, T. Mechanism Of Oxygen Evolution On Perovskites. J. Phys. Chem 1983, 87, 2960–2971. 15. Sunarso, J.; Hashim, S. S.; Zhu, N.; Zhou, W. Perovskite Oxides Applications In High Temperature Oxygen Separation, Solid Oxide Fuel Cell And Membrane Reactor: A Review. Prog. Energy Combust. Sci. 2017, 61, 57–77. 16. Sun, H.; Chen, G.; Sunarso, J.; Dai, J.; Zhou, W.; Shao, Z. A Molybdenum And Niobium Co-Doped B-Site Ordered Double Perovskite Catalyst For Efficient Oxygen Evolution Reaction. ACS Appl. Mater. Interfaces 2018, 10, 16939–16942. 18

ACS Paragon Plus Environment

Page 18 of 24

Page 19 of 24 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

17. Jiang, M.; Li, J.; Zhao, Y.; Pan, L.; Cao, Q.; Wang, D.; Du, Y. Double Perovskites as Model Bifunctional Catalysts toward Rational Design: The Correlation between Electrocatalytic Activity and Complex Spin Configuration. ACS Appl. Mater. Interfaces 2018, 10, 19746–19754. 18. De Angelis, F.; Meggiolaro, D.; Mosconi, E.; Petrozza, A.; Nazeeruddin, M. K.; Snaith, H. J. Trends in Perovskite Solar Cells and Optoelectronics: Status of Research and Applications from the PSCO Conference. ACS Energy Lett. 2017, 2, 857–861. 19. Du, J.; Zhang, T.; Cheng, F.; Chu, W.; Wu, Z.; Chen, J. Nonstoichiometric Perovskite Camno3- δ For Oxygen Electrocatalysis With High Activity. Inorg. Chem. 2014, 53, 9106–9114. 20. Oh, S. H.; Black, R.; Pomerantseva, E.; Lee, J.-H.; Nazar, L. F. Synthesis of a metallic mesoporous pyrochlore as a catalyst for lithium–O2 batteries. Nat. Chem. 2012, 4, 1004– 1009. 21. Oh, S. H.; Nazar, L. F. Oxide Catalysts for Rechargeable High-Capacity Li–O2 Batteries. Adv. Energy Mater. 2012, 2, 903–910. 22. Kim, J.; Yin, X.; Tsao, K.-C.; Fang, S.; Yang, H. Ca2 Mn2 O5 As Oxygen-Deficient Perovskite Electrocatalyst For Oxygen Evolution Reaction. J. Am. Chem. Soc. 2014, 136, 14646–14649. 23. Jung, J.-I.; Risch, M.; Park, S.; Kim, M. G.; Nam, G.; Jeong, H.-Y.; Shao-Horn, Y.; Cho, J. Optimizing Nanoparticle Perovskite For Bifunctional Oxygen Electrocatalysis. Energy Environ. Sci. 2016, 9, 176–183. 24. Grimaud, A.; Diaz-Morales, O.; Han, B.; Hong, W. T.; Lee, Y.-L.; Giordano, L.; Stoerzinger, K. A.; Koper, M. T.; Shao-Horn, Y. Activating Lattice Oxygen Redox Reactions In Metal Oxides To Catalyse Oxygen Evolution. Nat. Chem. 2017, 9, 457–465.

19

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

25. Ashok, A.; Kumar, A.; Bhosale, R. R.; Almomani, F.; Malik, S. S.; Suslov, S.; Tarlochan, F. Combustion Synthesis Of Bifunctional LaMO3 (M= Cr, Mn, Fe, Co, Ni) Perovskites For Oxygen Reduction And Oxygen Evolution Reaction In Alkaline Media. J. Electroanal. Chem. 2018, 809, 22–30. 26. Maitra, U.; Naidu, B.; Govindaraj, A.; Rao, C. Importance Of Trivalency And The Eg1 Configuration In The Photocatalytic Oxidation Of Water By Mn And Co Oxides. Proc. Natl. Acad. Sci. USA 2013, 110, 11704–11707. 27. Giannozzi, P.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Chiarotti, G. L.; Cococcioni, M.; Dabo, I. Quantum Espresso: A Modular And OpenSource Software Project For Quantum Simulations Of Materials. J. Phys.: Condens. Matter 2009, 21, 395502–395508. 28. Barnes, T. A.; Kurth, T.; Carrier, P.; Wichmann, N.; Prendergast, D.; Kent, P. R.; Deslippe, J. Improved Treatment Of Exact Exchange In Quantum Espresso. Comput. Phys. Commun. 2017, 214, 52–58. 29. Ernzerhof, M.; Scuseria, G. E. Assessment of the Perdew–Burke–Ernzerhof ExchangeCorrelation Functional. J. Chem. Phys. 1999, 110, 5029–5036. 30. Adamo, C.; Barone, V. Toward reliable density functional methods without adjustable parameters: The PBE0 model. J. Chem. Phys. 1999, 110, 6158–6170. 31. Rappe, A. M.; Rabe, K. M.; Kaxiras, E.; Joannopoulos, J. Optimized Pseudopotentials. Phys. Rev. B 1990, 41, 1227–1230. 32. Dudarev, S.; Botton, G.; Savrasov, S.; Humphreys, C.; Sutton, A. Electron-Energy-Loss Spectra and the Structural Stability of Nickel Oxide: An LSDA+ U Study. Phys. Rev. B 1998, 57, 1505–1509.

20

ACS Paragon Plus Environment

Page 20 of 24

Page 21 of 24 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

33. Lee, Y.-L.; Gadre, M. J.; Shao-Horn, Y.; Morgan, D. Ab Initio Gga+ U Study Of Oxygen Evolution And Oxygen Reduction Electrocatalysis On The (001) Surfaces Of Lanthanum Transition Metal Perovskites Labo 3 (B= Cr, Mn, Fe, Co And Ni). Phys. Chem. Chem. Phys. 2015, 17, 21643–21663. 34. Pack, J. D.; Monkhorst, H. J. Special Points For Brillouin-Zone Integrations - A Reply. Phys. Rev. B 1977, 16, 1748–1751. 35. Methfessel, M.; Paxton, A. High-Precision Sampling For Brillouin-Zone Integration In Metals. Phys. Rev. B 1989, 40, 3616–3621. 36. Fletcher, R. Practical Methods Of Optimization; John Wiley & Sons, 2013. 37. Arias, T. A.; Payne, M.; Joannopoulos, J. Ab Initio Molecular Dynamics: Analytically Continued Energy Functionals And Insights Into Iterative Solutions. Phys. Rev. Lett. 1992, 69, 1077–1080. 38. Zhao, B.; Zhang, L.; Zhen, D.; Yoo, S.; Ding, Y.; Chen, D.; Zhang, Q.; Doyle, B.; Xiong, X.; Liu, M. A Tailored Double Perovskite Nanofiber Catalyst Enables Ultrafast Oxygen Evolution. Nat. Commun. 2017, 8, 14586. 39. Nørskov, J. K.; Rossmeisl, J.; Logadottir, A.; Lindqvist, L.; Kitchin, J. R.; Bligaard, T.; Jonsson, H. Origin Of The Overpotential For Oxygen Reduction At A Fuel-Cell Cathode. J. Phys. Chem B 2004, 108, 17886–17892. 40. Tripković, V.; Skúlason, E.; Siahrostami, S.; Nørskov, J. K.; Rossmeisl, J. The oxygen reduction reaction mechanism on Pt (1 1 1) from density functional theory calculations. Electrochimica Acta 2010, 55, 7975–7981. 41. Montoya, J. H.; Garcia-Mota, M.; Nørskov, J. K.; Vojvodic, A. Theoretical Evaluation Of The Surface Electrochemistry Of Perovskites With Promising Photon Absorption Properties For Solar Water Splitting. Phys. Chem. Chem. Phys. 2015, 17, 2634–2640. 21

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

42. Kushima, A.; Yip, S.; Yildiz, B. Competing Strain Effects In Reactivity Of LaCoO3 With Oxygen. Phys. Rev. B 2010, 82, 115435–1–115435–6. 43. Wang, Y.; Cheng, H.-P. Oxygen Reduction Activity On Perovskite Oxide Surfaces: A Comparative First-Principles Study Of LaMnO3, LaFeO3, and LaCrO3. J. Phys. Chem C 2013, 117, 2106–2112. 44. Lu, Y.; Lu, F.; Yang, Z.; Wu, J.; Yu, H.; Xie, X.; Xu, J.; Cheng, F.; Chen, J.; Xiong, K. Elucidating Dz 2 Orbital Selective Catalytic Activity In Brownmillerite Ca2 Mn2 O5 . AIP Adv. 2016, 6, 095210–095216. 45. SAccoNI, L. Five Coordination In 3d Metal Complexes. Pure Appl. Chem. 1968, 17, 95–128. 46. Grimaud, A.; May, K. J.; Carlton, C. E.; Lee, Y.-L.; Risch, M.; Hong, W. T.; Zhou, J.; Shao-Horn, Y. Double Perovskites As A Family Of Highly Active Catalysts For Oxygen Evolution In Alkaline Solution. Nat. Commun. 2013, 4, 2439–2446. 47. Lee, J. G.; Hwang, J.; Hwang, H. J.; Jeon, O. S.; Jang, J.; Kwon, O.; Lee, Y.; Han, B.; Shul, Y.-G. A New Family Of Perovskite Catalysts For Oxygen-Evolution Reaction In Alkaline Media: BaNiO3 And BaNi0. 83O2. 5. J. Am. Chem. Soc. 2016, 138, 3541–3547. 48. Rossmeisl, J.; Logadottir, A.; Nørskov, J. K. Electrolysis Of Water On (oxidized) Metal Surfaces. Chem. Phys. 2005, 319, 178–184. 49. Han, B.; Viswanathan, V.; Pitsch, H. First-Principles Based Analysis Of The Electrocatalytic Activity Of The Unreconstructed Pt (100) Surface For Oxygen Reduction Reaction. J. Phys. Chem. C 2012, 116, 6174–6183. 50. Bothra, P.; Pati, S. K. Activity of Water Oxidation on Pure and (Fe, Ni, and Cu)Substituted Co3O4. ACS Energy Lett. 2016, 1, 858–862.

22

ACS Paragon Plus Environment

Page 22 of 24

Page 23 of 24 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

51. Hong, W. T.; Risch, M.; Stoerzinger, K. A.; Grimaud, A.; Suntivich, J.; Shao-Horn, Y. Toward The Rational Design Of Non-Precious Transition Metal Oxides For Oxygen Electrocatalysis. Energy Environ. Sci. 2015, 8, 1404–1427. 52. Man, I. C.; Su, H.-Y.; Calle-Vallejo, F.; Hansen, H. A.; Martínez, J. I.; Inoglu, N. G.; Kitchin, J.; Jaramillo, T. F.; Nørskov, J. K.; Rossmeisl, J. Universality in Oxygen Evolution Electrocatalysis on Oxide Surfaces. Chem. Cat. Chem. 2011, 3, 1159–1165. 53. Gupta, S.; Kellogg, W.; Xu, H.; Liu, X.; Cho, J.; Wu, G. Bifunctional Perovskite Oxide Catalysts For Oxygen Reduction And Evolution In Alkaline Media. Chem. Asian J. 2016, 11, 10–21. 54. Nørskov, J. K.; Bligaard, T.; Rossmeisl, J.; Christensen, C. H. Towards The Computational Design Of Solid Catalysts. Nat. Chem. 2009, 1, 37–46. 55. Greeley, J.; Stephens, I.; Bondarenko, A.; Johansson, T. P.; Hansen, H. A.; Jaramillo, T.; Rossmeisl, J.; Chorkendorff, I.; Nørskov, J. K. Alloys Of Platinum And Early Transition Metals As Oxygen Reduction Electrocatalysts. Nat. Chem. 2009, 1, 552–559. 56. Adler, J.; Parmryd, I. Quantifying Colocalization By Correlation: The Pearson Correlation Coefficient Is Superior To The Mander’s Overlap Coefficient. Cytometry A 2010, 77, 733–742. 57. Zhou, T.; Cao, Z.; Zhang, P.; Ma, H.; Gao, Z.; Wang, H.; Lu, Y.; He, J.; Zhao, Y. Transition Metal Ions Regulated Oxygen Evolution Reaction Performance Of Ni-Based Hydroxides Hierarchical Nanoarrays. Sci. Rep. 2017, 7, 46154–46159. 58. Hong, W. T.; Welsch, R. E.; Shao-Horn, Y. Descriptors Of Oxygen-Evolution Activity For Oxides: A Statistical Evaluation. J. Phys. Chem C 2015, 120, 78–86.

23

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

Graphical TOC Entry

24

ACS Paragon Plus Environment

Page 24 of 24