Full text

Turn on search term navigation

ARTICLE

Received 28 Apr 2015 | Accepted 31 Jul 2015 | Published 14 Sep 2015

Hyunah Kim1,*, Jimin Park1,2,*, Inchul Park1, Kyoungsuk Jin1, Sung Eun Jerng3, Sun Hee Kim3, Ki Tae Nam1

& Kisuk Kang1,4

The development of efcient and stable water oxidation catalysts is necessary for the realization of practically viable water-splitting systems. Although extensive studies have focused on the metal-oxide catalysts, the effect of metal coordination on the catalytic ability remains still elusive. Here we select four cobalt-based phosphate catalysts with various cobalt- and phosphate-group coordination as a platform to better understand the catalytic activity of cobalt-based materials. Although they exhibit various catalytic activities and stabilities during water oxidation, Na2CoP2O7 with distorted cobalt tetrahedral geometry shows high activity comparable to that of amorphous cobalt phosphate under neutral conditions, along with high structural stability. First-principles calculations suggest that the surface reorganization by the pyrophosphate ligand induces a highly distorted tetrahedral geometry, where water molecules can favourably bind, resulting in a low overpotential (B0.42 eV). Our ndings emphasize the importance of local cobalt coordination in the catalysis and suggest the possible effect of polyanions on the water oxidation chemistry.

DOI: 10.1038/ncomms9253 OPEN

Coordination tuning of cobalt phosphates towards efcient water oxidation catalyst

1 Department of Materials Science and Engineering, Research Institute of Advanced Materials (RIAM), Seoul National University, 1 Gwanak-ro, Gwanak-gu, Seoul 151-742, Republic of Korea. 2 Center for Biomaterials, Korea Institute of Science and Technology, 5, Hwarang-ro 14-gil, Seongbuk-gu, Seoul 136-791, Republic of Korea. 3 Western Seoul Center, Korea Basic Science Institute (KBSI), 150, Bukahyeon ro, Seodaemun-gu, Seoul 120-140, Korea. 4 Center for Nanoparticle Research, Institute for Basic Science (IBS), Seoul National University, 1 Gwanak-ro, Gwanak-gu, Seoul 151-742, Republic of Korea.* These authors contributed equally to this work. Correspondence and requests for materials should be addressed to K.T.N. (email: mailto:[email protected]

Web End [email protected] )or to K.K. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 6:8253 | DOI: 10.1038/ncomms9253 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9253

Splitting water into hydrogen and oxygen molecules via solar energy has been considered as one of the most environment-friendly ways to effectively use renewable

energy resources from harvest to redistribution16. Although it has been studied for more than a half century, the inefciency of the oxygen evolution reaction (OER) has not yet been resolved and is still regarded as a bottleneck in the integration of an overall water splitting system18. This inefciency is primarily due to the rigid OO double bond formation and the sluggish proton-coupled electron-transfer reactions79. Noble-metal catalysts, such as Pt, IrOx and RuOx, exhibit outstanding OER catalytic activity; however, the high costs of these catalysts prohibit their practical use1014. In this regard, the development of cost-effective, abundant-element-based and efcient OER catalysts is highly demanded.

Inspired from the Mn-based water oxidation complex in photosystem II that generates oxygen from water with an exceptionally high efciency1517, numerous OER catalysts have been synthesized and evaluated using relatively cheaper 3d transition metals (Mn, Fe, Co and Ni)8,1836. In particular, cobalt-containing catalysts have been studied as promising OER catalysts because of their high catalytic activity for water oxidation, which is comparable to those of precious-metal-based catalysts8,2336. Various cobalt-based oxides, perovskites and layered double hydroxides have been developed and examined via experimental and theoretical approaches2336. Among these materials, cobalt phosphate (Co-Pi) prepared using electrodeposition exhibits one of the highest performances to date with efcient and stable catalytic performance8. The extraordinary performance of Co-Pi has cast fundamental questions on the correlation between cobalt chemistry and the OER mechanism, and has further promoted the exploration of other cobalt-based OER catalysts.

Although various cobalt-based materials have been investigated as OER catalysts, researches were carried out mostly on materials with specic cobalt local environments such as octahedral (Oh)

cobalt8,2628,3234. Even amorphous Co-Pi has been recently shown to adopt the local Oh coordination of cobalt26. Although the coordination and corresponding electronic structure of the transition metal in the catalysts is believed to greatly affect the OER activity1923,3134, there is a lack of systematic information relating the local cobalt coordination to the OER catalysis.

Moreover, in practice, most of the cobalt-based materials reported thus far have been prone to transform into amorphous Co-Pi or cobalt-oxide-like phases from the surface under the typical pH and voltage regions where the OER proceeds2731. This transformation makes it more difcult to apprehend the correlation between the structural cobalt geometry and the OER mechanisms, which has in turn retarded the development of new cobalt-based OER catalysts.

Transition metal phosphates have positive attributes in structural stability, especially in oxidative environments over oxides, and can adopt versatile structures with various orientations of the phosphate group present in the crystal3739. As shown in Supplementary Table 1, whereas conventional oxide materials have cobalt Oh geometry with higher-symmetry crystal structures, phosphate/pyrophosphate-containing compounds have various cobalt geometries including tetrahedral (Td), trigonal bipyramidal (TBP) and Oh, with lower-symmetry crystal structures due to the bulky phosphate/pyrophosphate ligands. Moreover, the exible coordination of phosphate groups can stabilize the intermediate state of the transition metal by changing their local positions with ease, ensuring an effective redox change in the transition metal20,40. Indeed, our recent study using Mn3(PO4)2 3H2O as a OER catalyst demonstrated that the

inherently unstable Mn(III) state is stabilized by the structural exibility of the phosphate group20. Along with the case of Co-Pi30, it implies that the interplay between the transition metal and phosphate groups in the structure can also affect the catalytic activity, which, if clearly understood, would aid in the discovery of new efcient OER catalysts.

Herein, we explore the OER catalytic capabilities of four cobalt-based phosphate materials with various cobalt environments, Na2CoP2O7, NaCoPO4, Li2CoP2O7 and LiCoPO4, which are well known as cobalt redox-active electrode materials in rechargeable lithium or sodium batteries4145. Except for LiCoPO4, other three materials have not been evaluated as OER catalysts. All of these materials have crystal structures in which cobalt polyhedrons are cross-linked by phosphate groups. The uniqueness of this selection of materials lies in the various local cobalt coordination, including Oh, Td and TBP congurations, with different networks between the cobalt polyhedron. These different cobalt environments provide a platform to better understand the

Na2CoP2O7 Li2CoP2O7 Edge sharing (Oh, TBP)

Isolated (Td)

a c

c Coordinated4 Co2+

5 Coordinated Co2+

NaCoPO4 Edge sharing (Oh) Corner sharing (Oh)

6 6 CoordinatedCo2+

Coordinated

6 Co2+

LiCoPO4

Coordinated Co2+

b c

Figure 1 | Crystal structure and cobalt crystal eld of cobalt phosphate catalysts. Crystal structures of (a) Na2CoP2O7, (b) Li2CoP2O7, (c) NaCoPO4 and (d) LiCoPO4. The inset shows the local environment around the cobalt subunit (blue) and the spin state of the cobalt atom. The phosphate units, Na atoms and Li atoms are depicted in grey, yellow and green, respectively.

2 NATURE COMMUNICATIONS | 6:8253 | DOI: 10.1038/ncomms9253 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9253 ARTICLE

catalytic activity of cobalt-based materials in the OER. Interestingly, among the four candidates, we nd that Na2CoP2O7 with Td geometry has excellent catalytic activity and high stability comparable to those of the well-known amorphous Co-Pi. Cyclic voltammetry (CV), high-resolution transmission electron microscopy (HR-TEM), X-ray photo-electron spectroscopy (XPS) and in-situ X-ray absorption near edge spectroscopy (XANES), scanning transmission electron microscopy energy-dispersive X-ray spectroscopy (STEM-EDX) and electron paramagnetic resonance (EPR) conrms the catalytic and structural stability of Na2CoP2O7. Theoretical calculations reveal that a unique OER pathway is present for Na2CoP2O7, which uses four coordinated cobalt atoms with low-energy barriers for the water-splitting reaction.

ResultsCobalt phosphates as a new group of OER catalysts. The Na2CoP2O7, Li2CoP2O7, NaCoPO4 and LiCoPO4 compounds were synthesized using conventional solid-state methods, as described previously4145, and the X-ray diffraction (XRD) patterns of the synthesized compounds well matched those of previous reports without impurites4145 (Supplementary Fig. 1). The local crystal structure and overall networks of the polyhedrons in these four materials are schematically shown in Fig. 1 based on Rietveld analysis of each material along with the lattice parameters (Supplementary Fig. 1 and Supplementary Table 2).

The phosphate and pyrophosphate groups in these four materials induced various cobalt coordination and local polyhedron networks. Specically, Na2CoP2O7 and Li2CoP2O7 have

Td and mixed Oh/TBP cobalt geometries, respectively, whereas the cobalt in NaCoPO4 and LiCoPO4 was only in the Oh. Td polyhedron with higher point symmetry in Na2CoP2O7 is not usual considering the fact that Co(II) is typically stabilized in Oh environments2729. Moreover, the Td polyhedrons are isolated by pyrophosphate groups, and with respect to the long-range ordering of the polyhedron the Td polyhedra in Na2CoP2O7 form a two-dimensional layer via pyrophosphate groups, which is alternately stacked with sodium layers. Similarly, NaCoPO4 exhibited a two-dimensional nature; however, each edge-sharing Oh layer was cross-linked by a phosphate ligand. LiCoPO4 exhibited similar structures as NaCoPO4 in that two-dimensional chains of corner-sharing CoO6 were linked by a phosphate ligand. However, in Li2CoP2O7, Co2O9 subunits were interconnected through P2O7 groups to build a 3D framework structure. These various structural features differ from the well-studied edge-shared Oh CoO6 observed in most cobalt-based catalysts8,2628,3234.

OER catalytic behaviour of the cobalt phosphate catalysts. The water oxidation catalytic properties of Na2CoP2O7, Li2CoP2O7,

NaCoPO4 and LiCoPO4 were evaluated by CV in 0.5 M sodium phosphate buffer at pH 7.0 using a previously reported method21,23. One hundred cycles of potential sweep from 0.7 to1.5 V versus normal hydrogen electrode (NHE) were performed for each substrate (Supplementary Fig. 2). The 1st and 100th CV curves of Na2CoP2O7, Li2CoP2O7, NaCoPO4 and LiCoPO4 are shown in Fig. 2a. We observed that the catalytic activities increased in the series of LiCoPO4 oLi2CoP2O7 oNaCoPO4 oNa2CoP2O7. Na2CoP2O7 exhibited a current density of2.62 mA cm 2 at 1.44 V versus NHE, which is signicantly higher than those of NaCoPO4 (0.90 mA cm 2), Li2CoP2O7(0.14 mA cm 2) and LiCoPO4 (0.08 mA cm 2) under identical conditions. The superior catalytic activity of Na2CoP2O7 compared with the other compounds was maintained even after 100 cycles. The activities of both Na2CoP2O7 and NaCoPO4

remained nearly unchanged during 100 cycles, whereas those of Li2CoP2O7 and LiCoPO4 slightly increased. The current densities of Li2CoP2O7 and LiCoPO4 changed from 0.14 to 0.29 mA cm 2 and from 0.08 to 0.19 mA cm 2 during 100 cycles at 1.44 V versus NHE, respectively. A similar current increase was previously observed in LiCoPO4 (ref. 27), which was explained by the phase transformation on the catalyst surface, as will be discussed in detail in the later section.

The observed catalytic trend was also conrmed for the Tafel and potentiostatic electrolysis (Fig. 2b,c and Supplementary Fig. 3). Figure 2b shows that the exchange current increases in the series of LiCoPO4 oLi2CoP2O7 oNaCoPO4 oNa2CoP2O7, where the exchange currents are high for NaCoPO4 and Na2CoP2O7 even after 100 cycles of CV scans. The Tafel slopes of Na2CoP2O7, NaCoPO4, Li2CoP2O7 and LiCoPO4 at the 100th cycle were calculated as 82, 85, 84 and 82 mV dec 1, respectively.

Figure 2c indicates that the current densities of Na2CoP2O7 and NaCoPO4 in the potentiostatic electrolysis were stably maintained at a high value under constant potential 1.4 V versus NHE; in contrast, those of Li2CoP2O7 and LiCoPO4 were almost an order lower even though they slightly increased at the beginning of the electrolysis. It was found that the catalytic performance of Na2CoP2O7 was almost comparable to that of the electro-deposited Co-Pi in Supplementary Fig. 4. Despite the lower non-catalytic current, which is related to the number of redox-active cobalt ions or electrochemically active surface area28,46, the OER current of Na2CoP2O7 was higher than that of the electro-deposited Co-Pi at 1.4 V versus NHE (Supplementary Fig. 4). Comparison with b-CoOOH whose local structure is similar to that of Co-Pi exhibited similar trends (see details in the later part of the manuscript). However, it should be noted that, as Na2CoP2O7 was synthesized using a conventional solid-state reaction and not by electrodeposition, a further study for comparing the catalytic abilities of the two samples using a similar preparation method is needed.

In Fig. 2d, we summarized the electrochemical performance of four materials with respect to their catalytic activities and stabilities, that is, OER current after 100 cycles and the ratio of minimum to maximum current during the 100 cycles. A higher OER current along the x axis indicates a more active nature of the catalyst. Likewise, a higher current ratio close to unity along the y axis indicates better stability. As shown in Fig. 2d, Na2CoP2O7 exhibits the best catalytic activity and stability among the four materials. Moreover, further long-term stability tests for Na2CoP2O7 revealed that the catalytic current can be stably maintained during 500 cycles of CV scans and 12 h of bulk electrolysis (Supplementary Fig. 5). In the following, a more detailed investigation is focused primarily on Na2CoP2O7.

For deeper investigation of the OER performance of Na2CoP2O7, CV curves of Na2CoP2O7 were recorded under various pH conditions from pH 5 to pH 13. As shown in Fig. 3a, the onset potential for Na2CoP2O7 shifted towards a more negative potential with an increase in pH. The catalytic current was stably maintained during the OER under basic condition, whereas degradation in the catalytic current was observed during the OER tests at pH 5 (Supplementary Fig. 6).

The OER kinetics of Na2CoP2O7 was examined by varying the pH and the phosphate concentration of the electrolytes (Fig. 3a). First, we investigated the pH-dependent potential change at a constant current density for understanding the dependence of the reaction rate on proton activity4749. A plot between the pH and applied potential at constant current density of 1 mA cm 2 exhibited linearity from pH 6.0 to 8.0 and the slope was calculated as approximately 67 mV per pH (Fig. 3b). The Tafel

slopes for pH 6, 6.5 and 7.5 were measured as B80 mV dec 1, similar to the value at pH 7 (Fig. 3c). Then, the dependence of the

NATURE COMMUNICATIONS | 6:8253 | DOI: 10.1038/ncomms9253 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9253

3.5

3.0

2.5

2.0

1.5

1.0

0.5

0.8 1.0 1.2 1.4 1.6

E (V vs NHE)

0.7

0.65

0.6

0.55

2.5 2.0 1.5 1.0 0.5 0 0.5

Log I (mA cm2catalyst)

100

Na2CoP2O7 NaCoPO4

LiCoPO4

Li2CoP2O7

1 1 100

1 100

100

I(mA cm2 catalyast)

Overpotential (V)

100

E(H2O/O2)

lnactive Active

2.0

1.5

1.0

0.5

0 0 2,500 5,000 7,500 10,000 Time (s)

0.1

1.4 V electrolysis

1.43 V 1.44 V

Na2CoP2O7

StableUnstable

I(mA cm2 catalyast)

Increasing

Na2CoP2O7

stable NaCoPO4 stable

Li2CoP2O7 LiCoPO4

Minimum current

Maximum current

1.0

0.75

0.5

0.25

0 0 0.5 1.0 1.5 2.0 2.5 OER activity at selected voltage I (mA cm2catalyst)

1.43 V 1.44 V1.43 V1.44 V

Li2CoP2O7

00 2,500

NaCoPO4

1.44 V

LiCoPO4

1.43 V

Figure 2 | Oxygen evolution performances of cobalt phosphate catalysts. (a) CV curves and (b) Tafel plots of Na2CoP2O7, Li2CoP2O7, NaCoPO4 and LiCoPO4 for the 1st and 100th cycle from 0.7 to 1.5 V versus NHE. The thermodynamic potential for water oxidation is marked at 0.816 V versus NHE (pH 7.0). (c) Bulk electrolysis was performed under a constant potential of 1.4 V versus NHE. The inset displays the enlarged bulk electrolysis proles of Li2CoP2O7 and LiCoPO4. (d) The electrochemical performances of Na2CoP2O7, Li2CoP2O7, NaCoPO4 and LiCoPO4 with respect to their catalytic activities and stabilities. OER current of 100th cycle at the selected voltage and the ratio of the minimum current to maximum current at the selected voltage during 100 cycles are depicted along the x axis and y axis, respectively. Error bars represent s.d. (n 5).

reaction rate on the proton activity was derived by following equation4749:

@E @pH

@E @logj

The Faradaic efciency of Na2CoP2O7 during OER was quantied using oxygen sensor system. About 60 mmol of O2 was produced during the electrolysis at 1.4 V versus NHE at a period of 5,000 s with a Faradaic efciency of around B100%, indicating the measured current by CV is originated from the OER (Fig. 4a). With its proven ability for OER, facile mass production of Na2CoP2O7 powder using conventional synthesis methods and electrode-fabrication techniques enable us to produce large-scale oxygen evolution systems. The catalytic currents can be increased by integrating stacks of cells or by enlarging the size of the conducting substrate, and 180 mA of OER current can easily be acquired by stacking two cells with a geometric area of 64 cm 2 (Supplementary Fig. 7).

Structural and chemical stability of Na2CoP2O7 during OER. To understand the relative stability of Na2CoP2O7 and NaCoPO4 compared with the other two compounds during OER cycles, XRD analysis was performed. The XRD patterns in Supplementary Fig. 8 reveal no appreciable peak changes after the electrolysis. For detecting the possible structural changes at the surface, we performed surface-sensitive HR-TEM and XPS analysis on the electrolysed samples. The HR-TEM image and Fast Fourier Transform of Na2CoP2O7 reveal the high crystallinity at the surface after 100 cycles (Fig. 4b). In addition, the NaCoPO4 exhibited high crystallinity, as shown in Supplementary Fig. 9. Point EDX analysis for bulk-electrolysed Na2CoP2O7 in Cs-corrected STEM mode showed no change in the Co/P ratio from the surface to the bulk region, further supporting their structural stabilities.

In contrast, an amorphous phase was formed on the surfaces of Li2CoP2O7 and LiCoPO4 after 100 potential cycles, whereas high crystallinity was found in pristine samples (Supplementary

Substituting @E@pHj 67 mV pH 1 and @E@logjpH

80 mV dec 1 into the equation (1) leads to @logj@pHE 0:84. This

result indicated that the reaction rate follows an approximately inverse rst-order dependence on the proton activity4749. In addition, the dependence of the reaction rate on the phosphate-buffer strength was also examined by varying the buffer strength from 100 mM to 1 M. No change in the catalytic current depending on the buffer strength was found, indicating zero-order dependence on the phosphate-buffer strength (Fig. 3d)47,48.

Although the measured Tafel slopes were somewhat far from2.3 RT F 1, typically 60B80 mV dec 1 of Tafel slopes indicate that a one electron-transfer step is involved before rate determinant step (RDS)48,49. The exact reason for the higher Tafel values of our materials than 2.3 RT F 1 is unclear, but one possible explanation is that limitation in ion or electron transport might inate the Tafel slope, as proposed in previous reports for LiCoPO4 and thick Co-Pi lms27.

Although there is some deviation between the theoretical values and our experimental values, based on the Tafel and pH-dependent analysis, we can postulate that a one-proton/ one-electron-involved equilibrium step may exist before the RDS in the OER mechanism of Na2CoP2O7 (refs 4749). Considering the calculated RDS of Na2CoP2O7 was Co(IV) O formation, as

described in a later section, we could expect that equilibrium between Co(II) and Co(III) is the prior step before the RDS in the OER.

@logj

@pH

4 NATURE COMMUNICATIONS | 6:8253 | DOI: 10.1038/ncomms9253 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9253 ARTICLE

Potential(V) @ 1mA cm2

1.5

1.4

1.3

0.45

0.40

0.35

pH 13

pH 12

6.0 6.5 7.0 7.5 8.0 pH

Current (a.u.)

Overpotential (V)

0.30

1.0

1.5

2.0

2.5

3.0

100 300

pH 11

pH 8pH 7.5 pH 7 pH 6.5 pH 6

pH 5

5 mA cm2catalyst

3.5 3.0 2.5

pH 6 pH 6.5 pH 7 pH 7.5 2.0 1.5 1.0 0.5

Log i (mA cm2catalyst)

0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 E (V vs NHE)

Log i(mA cm2 catalyst)

1.1 V 1.15 V 1.2 V 1.25 V

500 1,000 Phosphate buffer concentration (mM)

Figure 3 | OER catalytic performances of Na2CoP2O7 at various pH and phosphate buffer concentrations. (a) CV curves for Na2CoP2O7 from pH 5 to pH 13. (b) pH-dependent potential change at a constant current density of 1 mA cm 2. (c) Tafel plots of Na2CoP2O7 from pH 6 to pH 7.5.

(d) Phosphate buffer concentration dependency of the OER current at a constant applied potential (1.10, 1.15, 1.2 and 1.25 V versus NHE) at pH 7. Error bars represent s.d. (n 5).

Fig. 10). This observation is in good agreement with a previous study on LiCoPO4, which gradually amorphized from the surface, resembling the local structure of Co-Pi (ref. 27). Moreover, the point EDX spectra for the bulk-electrolysed Li2CoP2O7 and

LiCoPO4 showed that the Co/P ratio of the amorphized surface region is B2, which is well accordant to that of Co-Pi (refs 8,27). In contrast, the bulk crystalline regions of Li2CoP2O7 and

LiCoPO4 have Co/P ratios of 0.5 and 1, respectively, which were matched to their bulk compositional ratios. In addition, the CV curves for Li2CoP2O7, LiCoPO4 and Co-Pi showed similar redox features before the catalytic waves and onset-potential voltage (Supplementary Fig. 11). From the amorphous character and electrochemically similar behaviours of Li2CoP2O7, LiCoPO4 and

Co-Pi, it can be deduced that the surface structure undergo phase transformation into Co-Pi or Co-Pi analogues. Previous reports showed similar phase transformation into Co-Pi-like structures in LiCoPO4 and BaxSr1xCo0.8Fe0.2O3d during OER tests27,28. The high crystallinity and stability both at the surface and in the bulk of Na2CoP2O7 and NaCoPO4 are quite inspiring considering that most cobalt-containing materials undergo a transformation into an amorphous Co-Pi-like phase from the surfaces27,28,30.

In-situ XANES and XPS analysis of Na2CoP2O7 were performed to probe any irreversible variation of the cobalt valence states during the OER using an in-situ electrochemical cell20,50. The Co K-edge spectra of Na2CoP2O7 were measured from 2 to 3 h after initiation of electrolysis at 1.5 V versus NHE. The Co K-edge spectra remained nearly constant during and after bulk electrolysis at 1.5 V versus NHE (Fig. 4c), indicating that no major change in the cobalt oxidation state occurred. In contrast to

Na2CoP2O7, the Co K-edge spectra of the electrolysed Li2CoP2O7 and LiCoPO4 shifted towards the higher-energy region compared with that of the as-prepared sample, implying that cobalt oxidation occurred during the OER (Supplementary Fig. 12). Moreover, the Co 2p XPS peaks from the surface to a depth of 10 nm of the bulk-electrolysed Na2CoP2O7 and NaCoPO4 showed negligible variation (Fig. 4d and Supplementary Fig. 13). In contrast, the XPS depth proling from the surface to 12-nm depths of cycled Li2CoP2O7 and LiCoPO4 showed that the surface

Co 2p peaks shifted towards higher energy compared with those of the bulk, further supporting their structural change on electrolysis (Supplementary Fig. 13).

EPR analysis clearly showed the different structural stabilities of Na2CoP2O7, NaCoPO4, Li2CoP2O7 and LiCoPO4 during the

OER tests by direct observations of Co(II) and Co(IV) species. The continuous-wave X-band EPR spectrum for Na2CoP2O7 only exhibited a characteristic S 3/2, Co(II) signal at geffE5 after

bulk electrolysis at various potentials and times (Fig. 5a). Even 300 min of bulk electrolysis at 1.4 V versus NHE did not trigger a change in the Co EPR signal. Here, the signal at g 2 resulted

from the carbon support (Supplementary Fig. 14). Invariance of the EPR signals on electrolysis supported the high structural stability of Na2CoP2O7.

In contrast, the distinct EPR signal around geffE2.3, which can

be assigned to the low spin Co(IV) with its total spin S 1/2,

appeared on bulk electrolysis of Li2CoP2O7 at an applied voltage of above 1.2 V versus NHE. With the decrease in the Co(II) signal, the Co(IV) signal increased when we raised the applied potential or electrolysis time, suggesting that signicant oxidation

NATURE COMMUNICATIONS | 6:8253 | DOI: 10.1038/ncomms9253 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9253

[110]

(001)

(110)

Na2CoP2O7

(prisitine)

Na2CoP2O7

(after 100 cycles)

Faradic efficiency

Experimental Theoretical

After electrolysis

Stable

[111]

Oxygen (M)

(011)

(101)

(110)

Na CoP O (after 100 cycles)

+ 2

Point EDX - 1 Co : P ~ 1: 1.88

Point EDX - 2 Co : P ~ 1: 1.89

Point EDX - 3 Co : P ~ 1: 1.93

Co Co

1.4 V electrolysis

4,000 5,000

0 1,000 2,000 3,000

Electrolysis time (s)

0 2 4 6 8 10 0

2 4 6 8 10

Energy (keV)

As prepared

Normalized intensity (a.u.)

1.5

0.5

1.4 V electrolysis

During electrolysis

Intensity (a.u.)

10 nm

8 nm

6 nm

4 nm

2 nm

7,700 7,710 7,720

0.5

0 7,714 7,716 7,718

In-situ XANES

Surface

7,730

7,740 810 800 790 780 770

Energy (eV)

Binding energy (eV)

Figure 4 | Faradaic efciency and structural stability of Na2CoP2O7 during OER. (a) Faradic efciency of Na2CoP2O7 measured by an oxygen sensor system. Bulk electrolysis was performed in a gas-tight electrochemical cell with the O2 sensor placed in the headspace. After initiating electrolysisat 1.4 V versus NHE, the percentage of O2 measured in the sensor on the FOXY probe was converted into the partial pressure of O2 in the headspace. The amount of evolved O2 molecules detected during bulk electrolysis and the theoretical amount of evolved O2 assuming a Faradaic efciency of 100% are depicted by the blue line and red line, respectively. (b) HR-TEM images and FFTs at the surface region of Na2CoP2O7 before and after 100 continuous potential cycles from 0.7 to 1.5 V versus NHE (scale bar, 5 nm). STEM image of Na2CoP2O7 after potential cycles and point EDX spectrum from the surfaces to the inside region, showing similar Co/P ratios (scale bar, 20 nm). (c) In-situ XANES spectra corresponding to the Co K-edge of Na2CoP2O7 before, during and after bulk electrolysis at an applied potential of 1.5 V versus NHE for 2 h. The inset displays the enlarged XANES edge region. (d) XPS spectra of the Co 2p region of the surface of the bulk-electrolysed sample at an applied potential of 1.4 V versus NHE for 1 h.

of Co(II) into Co(IV) in Li2CoP2O7 occurred during the OER (Fig. 5b). Similar EPR trends were observed in previous reports using Co-Pi and they found that its catalytic resting state comprises Co(IV) species30,51. The increase in the Co(IV) signal with the concurrent decrease in Co(II) presented another evidence for structural change of Li2CoP2O7 into Co(IV)-containing Co-Pi-like structures during OER.

Similar with the case of Na2CoP2O7 and Li2CoP2O7, the EPR spectrum for bulk-electrolysed NaCoPO4 only showed a Co(II) signal at geffE5, whereas a distinct Co(IV) signal at geffE2.3 was found in bulk-electrolysed LiCoPO4 (Supplementary Fig. 15). In our studies the Co(III) signal did not appear, because the low-spin Co(III) with its total spin S 0 is EPR silent and high-spin

Co(III) are rare and only existed in weak-eld ligands, as described previously51. Although it is uncertain whether a small portion of Co(II) in Na2CoP2O7 or NaCoPO4 transformed into

Co(III) during the OER, based on negligible change in the Co(II) signal on bulk electrolysis and no appearance of a Co(IV) signal, we believe that their dominant Co species during OER may be Co(II) or Co(II)/Co(III) mixed states. In addition, these EPR results supported phase transformation of Li2CoP2O7 and

LiCoPO4 during the OER into Co(IV)-containing Co-Pi-like structures.

Calculated surface structures of Na2CoP2O7. To understand the catalytic activity of Na2CoP2O7 with respect to its local structure on an atomistic level, we performed a series of rst-principles calculations. The possible surface structures of Na2CoP2O7 were

rst calculated to estimate the equilibrium morphology at pH 7.0 (Supplementary Table 3). Based on the relative surface energies, the equilibrium crystal shape of Na2CoP2O7 was projected using

Wulff constructions, which revealed that Na2CoP2O7 is mainly composed of a (101) surface with slight (001) surface exposure (Fig. 6a). It is worth noting that in the main (101) surface of Na2CoP2O7, the pyrophosphate ligands rotated around the cobalt atoms and generated a new CoO bond, as illustrated in Fig. 6b. As shown in Supplementary Fig. 16, a general relaxation of surface atoms occurred and then new CoO bonds were formed by rotation of the pyrophosphate ligands. After new CoO bonds were formed, further relaxation occurred, creating the nal surface structures of Na2CoP2O7 (101) plane. It should be noted that the pyrophosphate-group rotation in surface spontaneously occurs following an energetically favourable pathway with the decrease of Gibbs free energy. The structural exibility of the P2O7 group could stabilize the under-coordinated transition metals on the surface40. As a result of the passivation by the pyrophosphate group, three-coordinated surface cobalt atoms initially formed after surface-cleaving step irreversibly transformed into a highly distorted tetrahedral geometry similar to the CoO5 TBP structure with a deciency of one oxygen.

OER catalytic mechanism of Na2CoP2O7. Given the well-established surface structures, we attempted to estimate the overpotential by sequentially calculating the Gibbs free energies for the four elementary steps occurring on the surface during the

6 NATURE COMMUNICATIONS | 6:8253 | DOI: 10.1038/ncomms9253 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9253 ARTICLE

Co(II) Na2CoP2O7

Li2CoP2O7

Intensity (a.u.) Intensity (a.u.)

No Co(IV)

0.8 V 30 min

0 1,000

Co(II)

1.2 V 30 min

1.4 V 300 min

0.8 V 30 min1.2 V 30 min

1.4 V 300 min

2,000 3,000 4,000 5,000 6,000Field (G)

0 1,000 2,000 3,000 4,000 5,000 6,000Field (G)

Co(IV)

1.4 V 30 min

Figure 5 | Continuous wave X-band EPR spectra for Na2CoP2O7 and Li2CoP2O7 on bulk electrolysis. (a) EPR spectra for Na2CoP2O7 at 0.8 V for 30 min (black), 1.2 V for 30 min (green) and 1.4 V for 300 min (blue). Only Co(II) signals at geffE5 appeared on bulk electrolysis. Here, the signal at g 2 is originated from carbon. (b) EPR spectra for Li2CoP2O7 at 0.8 V for

30 min (black), 1.2 V for 30 min (green), 1.4 V for 30 min (orange) and 1.4 V for 300 min (blue). Contrary to Na2CoP2O7, a decrease in Co(II) signals with a concurrent increase in Co(IV) signals (geffE2.3) was observed with

the increase in applied voltage and time.

OER based on acidbase mechanism18,22,31. These four steps are used for modelling the thermochemistry of the OER18,22,31 and are described as follows:

H2O ! OH e H 2

OH ! O e H 3

O H2O ! OOH e H 4

OOH ! O2 e H 5 (*represents the active surface site where OH and O species can adsorb).

During the process, water molecules rst adsorb on the active site on the surface and dissociate into a proton and OH* (step 1), which again splits into a proton and O* (step 2). The oxygen atom recombines with a water molecule to form an OOH* bond (step 3) and nally evolves as O2 gas (step 4). Our calculation model has limitations in neglecting the inuence of the electric eld on the double layer and the possible differences resulting from the types of proton donor (H3O or H2O). According to previous reports18,52,53, with this assumption, the Gibbs free energy of each step during the OER can be calculated as identical, regardless of the pH value. In this regard, we simply calculated

the Gibbs free energy of each step in Na2CoP2O7 using the four-step model in acidic conditions.

Moreover, it should be noted that an alternative OER mechanism exists, the so-called direct coupling mechanism54,55. Recently, Wang and Van Voorhis54 revealed the new OER mechanism in cobalt oxide, demonstrating direct-coupling OO bond formation can occur with a very low activation barrier, and steps other than OO bond formation act as the RDS in cobalt oxides. In this stage, based on the relatively larger inter-cobalt distance in Na2CoP2O7 (5.15 ) compared with that in cobalt oxide (B2.85 ) revealed by in-situ XAS analysis26, we focused on the acidbase mechanism for exploring the OER pathway in Na2CoP2O7.

During these four steps, the cobalt ion is sequentially oxidized from 2 to 3 (step 1) and 3 to 4 (step 2), and is reduced from 4 to 3 (step 3) and 3 to 2 (step 4), indicating the active participation of the cobalt redox reaction at every step of the OER. The detailed energy changes during the four elementary steps are shown for the two dominant surfaces, (101) and (001) surfaces, in Fig. 6c and Supplementary Fig. 17, respectively. Here, the dotted line represents the calculated Gibbs free energies at ideal charge transfer steps and the blue line represents those of Na2CoP2O7 when U 0, both at pH 7.0. The red line represents the net Gibbs

free energy changes of Na2CoP2O7 when U 0.817, which

represents the additionally required energy for the reaction to proceed (equivalent to the overpotential) in each step during the OER. Probing the Gibbs free energies for these steps reveals that the (101) surface exhibits a signicantly lower overpotential of0.417 eV than that of the (001) surface of 1.117 eV (Fig. 6c and Supplementary Fig. 17). It is noteworthy that the magnied CV curves in Supplementary Fig. 4 show that the current increased dramatically near the overpotential of 0.42 V, supporting the validity of our calculations, although experimental overpotential may slightly vary with the current density.

The individual steps of the water oxidation are closely investigated in the following. As a rst step, we monitored the Gibbs free energy change for OH* formation (step 1) on the (101) surface in our system. As the absorption of water molecules on the surface is the initial step of the OER, it is generally perceived that the absorption energy of water molecules highly affects the OER activity5260. Norskov and colleagues52,53,58 reported the importance of the adsorption energies of adsorbents, such as *OH and *OOH species, on the OER activities on metal-oxide surfaces and a universal scaling relation between *OH and *OOH adsorption energy was found. Recently, Liu et al.59 and Kim et al.60 showed that oxygen or cobalt vacancies with high binding energies with water molecules lead to a signicant enhancement of OER catalytic capabilities. Interestingly, the Gibbs free energy for step 1 in our system only requires 0.713 eV, which is lower than the energy needed (0.817 eV) to initiate water oxidation ideally at pH 7.0, implying a facile water adsorption on surface cobalt atoms in the (101) plane. This facile adsorption originates from the unique structural characteristic of the reorganized surface cobalt atoms; pyrophosphate ligands induced highly distorted Td coordination, which has an open site, where oxygen adsorbants can easily bind. Indeed, previous reports showed that similar distorted metal geometries with open coordination sites have an advantage in the uptake of oxygen adsorbants5962. This nding strongly suggests that the distorted local cobalt geometry induced by the reorganization of pyrophosphate ligands contributed to a favourable water adsorption and further water oxidation catalysis.

Compared with facile OH* formation, O* formation (step 2) requires the highest barrier (0.417 eV) among the four steps, indicating that this step is the rate-limiting process of the water oxidation catalysis of Na2CoP2O7. This theoretical value is

NATURE COMMUNICATIONS | 6:8253 | DOI: 10.1038/ncomms9253 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9253

3-Coordinated cobalt atom

Before P2O7 rotation

4-Coordinated cobalt atom

P2O7 P2O7

After P2O7 rotation

Co(III)

(101) Surface

Active site

Co(II) 4 Schematic of OER on (101) surface pH 7

3 OOH*

2H2O +*

O* + H2O

G (eV)

OH* + H2O

2H2O +*

U = 0.0

= 0.417 U = 0.817

G = 0.713 eVG = 1.234 eVG = 0.929 eVG = 0.392 eV

1 Reaction coordinate

Co(III)

Co(IV)

Figure 6 | Density functional theory calculation for OER mechanism of Na2CoP2O7. (a) Wulff structure of Na2CoP2O7, showing the equilibrium crystal shape based on the relative surface energies and the enlarged (101) surface with active sites. (b) Local environment around the active cobalt atom located on the (101) surface before and after P2O7 rotation; the active cobalt atom is depicted in pink. (c) Schematic of the Gibbs free energy changes for the four elementary steps during the OER based on DFTcalculations. The local structure and valence change of the cobalt atoms in the active sites are shown at the corners. The inset shows the free-energy landscape on the (101) surface of Na2CoP2O7 compared with an ideal catalyst at pH 7.0. An ideal catalyst,

Na2CoP2O7 for U 0 and Na2CoP2O7 for U 0.817 are depicted by the dashed line, blue line and red line, respectively.

comparable in magnitude to that of the active cobalt-oxide b-CoOOH (0.48 eV), which was selected as a representative for amorphous Co-Pi (ref. 31). Considering that RDS in conventional cobalt oxides including b-CoOOH is generally known as OOH* formation30,31, the limiting step of O* formation in Na2CoP2O7 is worthwhile to be noticed. However, it should be noted that steps other than OO bond formation can be the RDS by direct OO coupling54. We speculate that the highest barrier residing in step 2 of Na2CoP2O7 is related to the inherent nature of the high cobalt oxidation potential in the phosphate framework42,43,45. Much higher redox potentials for cobalt (Co(II)/Co(III) or Co(III)/Co(IV)) are generally observed in polyanion frameworks such as phosphates and pyrophosphates compared with oxides because of the inductive effect42,43,45. As step 2 requires the oxidation of the Co(III) to Co(IV) in pyrophosphates, a substantially higher energy will be required to oxidize the Co(III) ion than in oxides42,43,45. It is noteworthy that the redox potentials of cobalt ions in oxides typically shift up by B1 V in phosphates42,45. The higher energy for oxidizing

Co(III) in Na2CoP2O7 may make O* formation the RDS.

Another interesting characteristic of Na2CoP2O7 was found in the adaptable local cobalt structure during the OER (corners of Fig. 6c). Although the OER proceeds in ve- and six-coordinated polyhedron in conventional cobalt oxides26,30,31, Na2CoP2O7 used four- and ve-coordination geometries accompanying the change of cobalt valence. Specically, with the valence change from 2 to 3 in step 1, the oxygen coordination number of cobalt increased from 4 to 5 by adsorbing the OH species. Then, further oxidation to Co(IV) occurs by detaching the proton from the adsorbed OH and the remaining oxygen bound to Co(IV) recombines with a water molecule to make OOH, reducing the

cobalt to III. During these two redox steps (step 2 and 3), the ve-coordination geometry was retained without a signicant change in the atom positions. As under-coordinated transition metals are generally unstable and are prone to induce surface reconstruction63,64, the stable performance of the four/ve coordination of cobalt is remarkable. This performance was possible mainly because of the rotation of the exible pyrophosphate group, which readily forms additional CoO bonds at the surface, indicating the important interplay between the cobalt ion and phosphate polyanion.

We performed additional computations to explore the theoretical OER mechanism for NaCoPO4 and b-CoOOH (Supplementary Fig. 18). b-CoOOH was introduced for representing Co-Pi-like structures at the surface of Li2CoP2O7 and

LiCoPO4 during the OER. The four OER steps of the dominant NaCoPO4 surface, the (010) plane, showed that theoretical overpotential of NaCoPO4 was 0.60 V with RDS of O* formation. In case of b-CoOOH, the 10 14

plane was selected as the most

catalytically active surface, on the basis of a previous study based on the acidbase mechanism. The theoretical overpotential of b-CoOOH was calculated as 0.44 V with RDS of *OOH formation, which was in line with a previous report31.

Finally, we compared the similar and different characteristics of Na2CoP2O7 and Co-Pi (b-CoOOH) in Table 1. Whereas these two types of catalyst exhibited similar electrochemical performances, their underlying catalytic mechanisms showed strikingly different properties. In terms of structure, Na2CoP2O7 has four/ve-coordination geometry inside the crystalline structure during OER, whereas Co-Pi has Oh-based ve/six coordination in an amorphous structure8,30,31. The dominant Co species of Na2CoP2O7 during OER is Co(II) or Co(II)/Co(III), whereas

8 NATURE COMMUNICATIONS | 6:8253 | DOI: 10.1038/ncomms9253 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9253 ARTICLE

Table 1 | Comparison between Na2CoP2O7 and CoOx.

Catalyst Na2CoP2O7 CoOx (Co-Pi, b-CoOOH) Refs Electrochemical performance

Catalyst type Powder (crystalline) Electrodeposited thin lm (amorphous) 8

Experimental overpotential value 560 mV (1 mA cm 2) 410 mV (1 mA cm 2) 8

pH (electrolyte composition ) pH 7 (0.5 M phosphate) pH 7 (0.1 M phosphate) 8

Tafel slope Approximately 80B90 mVdec 1 60 mVdec 1 8, 30

pH dependence 67 mV 64 mV 47

Buffer concentration dependence Zeroth order dependence Zeroth order dependence 47

Structural/chemical property during OEREPR 2 or 2 /3 mixture 3 /4 mixture

In-situ XANES 2 Z3

HR-TEM at surface (EDX) High crystallinity (Co: P 1:2, TEM) Amorphous phase (Co: P 2B3:1, SEM)

OER mechanism

Cobalt coordination (Pristine state) 4 (Tetrahedral) 6 (Octahedral) 26, 30

Cobalt coordination (During OER) 425 526 31

OER active surface (101) 10 14

Rate determining Step OH*-O* O*-OOH* (acidbase mechanism) 31

Theoretical overpotential value 417 mV 480 mV (ref) 437 mV (our work) 31

EDX, scanning transmission electron microscopy; EPR, electron paramagnetic resonance; HR-TEM, high-resolution transmission electron microscopy; OER, oxygen evolution reaction; XANES, X-ray

absorption near edge spectroscopy.

that of Co-Pi is known as Co(III)/Co(IV)26,30. Finally, the theoretically calculated RDE of Na2CoP2O7 is *O formation, whereas that of Co-Pi is *OOH formation (acidbase mechanism)30,31. Based on these comprehensive comparisons, we think that Na2CoP2O7 can be a good model system to explore new OER chemistry beyond conventional Oh-based cobalt-oxide materials.

DiscussionIn summary, we explored the abilities of four phosphate/ pyrophosphate-based, well-known cobalt redox-active materials as water oxidation catalysts and found a highly active and stable material, Na2CoP2O7. Whereas previously reported cobalt-containing materials underwent signicant phase transformation during the OER, Na2CoP2O7 exhibited excellent catalytic activity and phase stability. Moreover, mechanistic studies for the OER of Na2CoP2O7 were performed by varying the pH and phosphate buffer strength of the electrolyte. Density functional theory calculations revealed that Na2CoP2O7 have lower theoretical overpotential (B0.42 eV) compared with the most active cobalt-oxide phase, CoOOH (B0.48 eV), theoretically verifying its high catalytic activity. In addition, we found the possible OER pathway in Na2CoP2O7 that used four- and ve-coordinated cobalt atoms stabilized by actively rotating pyrophosphate groups. Through the experimental and theoretical comparison, we found that the underlying catalytic mechanism of Na2CoP2O7 was notably different to that of the conventional cobalt oxides, in terms of RDS, local cobalt structure and cobalt valence state during the OER. We believe that using polyanion-based cobalt compounds can be an interesting new direction for the development of new efcient and stable catalysts, and the discovery of new catalytic pathways in various transition metal geometries. We expect that this result can broaden the current oxide-based picture of the water oxidation complex and give hints at a possible effect of the polyanion on OER catalytic abilities.

Methods

Material synthesis. Na2CoP2O7 and NaCoPO4 powders were synthesized using the conventional solid-state method. Stoichiometric amounts of Na2CO3 (ACS reagent, Z99%, Aldrich), CoC2O4 (Z99%, Kojundo) and (NH4)2HPO4 (ACS

reagent 98%, Aldrich) were wet-ball milled in acetone media and heated at 300 C under Ar media for 6 h. After grinding, the pelletized samples were calcined at 600 C and 750 C under Ar media for 6 h for Na2CoP2O7 and NaCoPO4, respectively. Li2CoP2O7 and LiCoPO4 powders were synthesized using the stoichiometric amounts of Li2CO3 (ACS reagent, Z99%, Aldrich), CoC2O4 (Z99%, Kojundo) and (NH4)2HPO4 (ACS reagent 98%, Aldrich). The precursors were mixed by ball milling, heated at 300 C under Ar media for 6 h, pelletized and calcined at 600 C for Li2CoP2O7 and 800 C for LiCoPO4 under

Ar media for 6 h. Amorphous Co-Pi was synthesized according to the previous report8.

Material characterization. Powder XRD patterns were measured using a Bruker D-8 Advance X-ray diffractometer using Cu Ka radiation (l 1.54056 ) from 10

to 70 with a step size of 0.01. Structural analysis of materials was performed by Rietveld renement method using FullProf software. TEM images were obtained using a eld-emission TEM (Philips Tecnai F20) operating at 200 kV. HR-TEM images and FFTs of the materials were used to analyse surface morphology. Point EDX analysis was performed under the Cs-corrected STEM mode and the spectra from three different regions were acquired with a Quanta 3D FEG (FEI,The Netherlands). The XAS experiments were performed at beamline 8C atthe Pohang Accelerator Laboratory, Republic of Korea. The oxidation state corresponding to the Co K-edge in each sample was investigated using XANES analysis. In-situ XANES spectra were collected in uorescence mode with an electron energy of 2.5 GeV and stored in the current of 250 mA top-up mode. XPS spectra were performed to examine the oxidation state of the Co 2p peaks from surface to depth of 10 nm of the catalyst after bulk electrolysis. XPS analysis was performed ve times at an increment of 2 nm etching from the surface. The data were attained with passing energy 23.5 eV with a step size of 0.05 eV by electron spectroscopy (Sigma Probe, Thermo VG Scientic, UK). BrunauerEmmettTeller analysis was performed to analyse the surface area of the catalyst. The sample was loaded in the BrunauerEmmettTeller analyser (Tristar II 3020 Micromeritics, USA) under an N2 adsorption environment.

Electrochemistry. The electrochemical measurements were performed in a three-electrode system using a potentiostat (CHI 760C, CH Instruments, Inc.).

Pt foil (2 cm 2 cm 0.1 mm, 99.997% purity, Alfa Aesar) was used as the counter

electrode and Ag/AgCl electrode (BASi, 3 M NaCl) was used as a reference electrode. The working electrode was prepared by a previously reported spin-coating method21,23. The catalyst with neutralized Naon solution was dispersed in puried deionized water using sonication and was spin-coated onto an uorine-doped tin oxide substrate. After the electrode was dried at 80 C for 5 min, the loaded catalyst weight was carefully measured using a micro weighing electronic scale (Sartorius Micro Balance). A pH 7.0 sodium phosphate buffer with buffer strength 0.5 M was used as the electrolyte, and the electrolyte was degassed with high-purity N2 before every electrochemical measurement. The OER curves were corrected using the averaged value of the forward and reverse CV currents for polarization compensation and further iR-compensation using V V

applied

was performed. After polarization and iR compensation, each OER value was

NATURE COMMUNICATIONS | 6:8253 | DOI: 10.1038/ncomms9253 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9253

following equations:

DG1 G OH

G mOH

E OH

E E H2O

1=2E H2

eU DGH pH

DZPE TDS

normalized to the total surface area of the catalysts on the substrate. The electrochemical potentials were converted to the NHE scale by following equation: E(NHE) E(Ag/AgCl) 0.197 V.

To prepare electrolytes with different pH values from 5.0 to 8.8, 0.5 M phosphate buffer was used with varying ratios of NaH2PO4 and Na2HPO4. Over pH 9, the pH was adjusted by addition of NaOH into the phosphate buffer (pH 8.8) in the customary manner. To determine the effect of the phosphate buffer concentration on the activity, different phosphate buffer concentrations from 0.1 to 1 M were prepared and different amounts of NaClO4 were added to maintain the same solution conductivity regardless of the phosphate buffer concentrations.

EPR measurement. X-band continuous-wave EPR spectra were obtained using a Bruker EMX/Plus spectrometer equipped with a dual-mode cavity (ER 4116DM). A liquid He quartz cryostat (Oxford Instruments ESR 900) combined with gas ow controller (Oxford Instruments ITC 503) was introduced to control the low temperatures (B5.7 K). For EPR measurement, 15 mg of Na2CoP2O7, NaCoPO4,

Li2CoP2O7 and LiCoPO4 catalyst electrodes were prepared as carbon-paste type. After bulk electrolysis at the specic voltages and times, the electrodes were removed from electrolyte, followed by quick drying using an Ar (99.999 %) gun. Then, the solid-state catalysts were carefully collected by a razor blade and were transferred to an EPR tube under Ar (99.999 %) atmosphere. The EPR tubewas immediately frozen and stored at 77 K using liquid nitrogen before EPR measurement. The elapsed time from the end of electrolysis to freezing of the EPR tube was B5 min. The following parameters were used: microwave frequency 9.64 GHz, modulation amplitude 10 G, modulation

frequency 100 kHz, microwave power 0.94 mW and temperature 5.7 K. T

en scans were added for obtaining each spectrum.

First-principles calculations. All the free energies in this work were calculated using the generalized gradient approximation (GGA) and the PerdewBurke Ernzerhof functional for the exchange correlation to the density functional theory65. The projected augmented wave method was used, as implemented in the Vienna ab initio simulation package66. The Hubbard U parameter (GGA U) with

U 3.5 eV was used to accurately calculate the electron correlation within the d

states in a cobalt ion6770. A 3D slab model with periodic boundary conditions was used with an energy cutoff of 500 eV and an appropriate g-point, only k-points mesh, was selected to ensure that the total energies converged within 5 meV per formula unit of Na2CoP2O7.

Wulff construction and surface properties. To calculate the surface energies and Wulff shape of Na2CoP2O7, the initial unrelaxed surface structures were carved out from the fully relaxed bulk structure of Na2CoP2O7 with P42mnm symmetry, which is a less distorted structure with higher symmetry than the experimental P21cn structure. The selection of surfaces to investigate was limited to low-index surfaces and guided by favouring surfaces with preservation of the PO4 tetrahedron, which led to the following set of surfaces: (100), (001), (110), (210), (101), (111) and 3 10

All the surface structures were calculated within the GGA U scheme as

dened above, using a supercell with atomic layers over 15 and vacuum slabs over 20 along the z direction using a convergence test to obtain the adsorption energies and to recover the bulk magnetic ordering of the centre layers (convergence within 1 meV 2 for the surface energies). The surface energies were calculated from equation (6) on either stoichiometric or non-stoichiometric slabs for which the top- and bottom-most cobalt pyrophosphate layers were relaxed below the force threshold of 0.05 eV 1.

Surface free energy : g Gslab NCoGbulk X

DG2 G O

mH G OH

E O

1=2E H2

E OH

eU DGH pH

DZPE TDS

DG3 G OOH

G O

mOH

E OOH

E O

E H2O

1=2E H2

eU DGH pH

DZPE TDS

DG4 4 1:23 eV eU DGH pH

DG1 DG2 DG3 13

where the density functional theory energies, zero point energy and entropy correction values are listed in the Supplementary Table 4.

References

1. Gratzel, M. Photoelectrochemical cells. Nature 414, 338344 (2001).2. Lewis, N. S. & Nocera, D. G. Powering the planet: chemical challenges in solar energy utilization. Proc. Natl Acad. Sci. USA 103, 1572915735 (2006).

3. Barber, J. Photosynthetic energy conversion: natural and articial. Chem. Soc. Rev. 38, 185196 (2009).

4. Gust, D., Moore, T. A. & Moore, A. L. Solar fuels via articial photosynthesis. Acc. Chem. Res. 42, 18901898 (2009).

5. Faunce, T. et al. Articial photosynthesis as a frontier technology for energy sustainability. Energy Environ. Sci. 6, 10741076 (2013).

6. Swierk, J. R. & Mallouk, T. E. Design and development of photoanodes for water-splitting dye-sensitized photoelectrochemical cells. Chem. Soc. Rev. 42, 23572387 (2013).

7. Betley, T. A., Wu, Q., Van Voorhis, T. & Nocera, D. G. Electronic design criteria for O O bond formation via metal oxo complexes. Inorg. Chem.

47, 18491861 (2008).8. Kanan, M. W. & Nocera, D. G. In situ formation of an oxygen-evolving catalyst in neutral water containing phosphate and Co2. Science 321, 10721075 (2008).

9. Zhao, Y. et al. Improving the efciency of water splitting in dye-sensitized solar cells by using a biomimetic electron transfer mediator. Proc. Natl Acad. Sci. USA 109, 1561215616 (2012).

10. Sartorel, A. et al. Polyoxometalate embedding of a tetraruthenium(IV)-oxo-core by template-directed metalation of [g-SiW10O36] 8 : a totally

inorganic oxygen-evolving catalyst. J. Am. Chem. Soc. 130, 50065007 (2008).11. Youngblood, W. J. et al. Photoassisted overall water splitting in a visible light-absorbing dye-sensitized photoelectrochemical cell. J. Am. Chem. Soc. 131, 926927 (2009).

12. Tong, L., Duan, L., Xu, Y., Privalov, T. & Sun, L. Structural modications of mononuclear ruthenium complexes: a combined experimental and theoretical study on the kinetics of ruthenium-catalyzed water oxidation. Angew. Chem. Int. Ed. 50, 445449 (2011).

13. Duan, L. et al. A molecular ruthenium catalyst with water-oxidation activity comparable to that of photosystem II. Nat. Chem. 4, 418423 (2012).

14. Joya, K. S., Subbaiyan, N. K., DSouza, F. & de Groot, H. J. Surface-immobilized single-site iridium complexes for electrocatalytic water splitting. Angew. Chem. Int. Ed. 51, 96019605 (2012).

15. Loll, B., Kern, J., Saenger, W., Zouni, A. & Biesiadka, J. Towards complete cofactor arrangement in the 3.0 resolution structure of photosystem II. Nature 438, 10401044 (2005).

16. Yano, J. et al. Where water is oxidized to dioxygen: structure of the photosynthetic Mn4Ca cluster. Science 314, 821825 (2006).

17. Umena, Y., Kawakami, K., Shen, J.-R. & Kamiya, N. Crystal structure of oxygen-evolving photosystem II at a resolution of 1.9 . Nature 473, 5560 (2011).

18. Su, H. Y. et al. Identifying active surface phases for metal oxide electrocatalysts: a study of manganese oxide bi-functional catalysts for oxygen reduction and water oxidation catalysis. Phys. Chem. Chem. Phys. 14, 1401014022 (2012).

19. Bergmann, A., Zaharieva, I., Dau, H. & Strasser, P. Electrochemical water splitting by layered and 3D cross-linked manganese oxides: correlating structural motifs and catalytic activity. Energy Environ. Sci. 6, 27452755 (2013).

20. Jin, K. et al. Hydrated manganese(II) phosphate (Mn(3)(PO(4))(2).3H(2)O) as a water oxidation catalyst. J. Am. Chem. Soc. 136, 74357443 (2014).

21. Park, J. et al. A new water oxidation catalyst: lithium manganese pyrophosphate with tunable Mn valency. J. Am. Chem. Soc. 136, 42014211 (2014).

22. Friebel, D. et al. Identication of highly active Fe sites in (Ni,Fe)OOH for electrocatalytic water splitting. J. Am. Chem. Soc. 137, 13051313 (2015).

23. Gardner, G. P. et al. Structural requirements in lithium cobalt oxides for the catalytic oxidation of water. Angew. Chem. Int. Ed. 51, 16161619 (2012).

Ni xiNCo

mi =2A; 6

with surface area 2A, free energy Gslab, chemical potentials mi, the number of excess adsorbed species over bulk Ni xiNCo (where xi is the number of atoms per bulk

formula).

The chemical potentials mi of the adsorbents were taken relative to the free energy of liquid water and hydrogen gas, written explicitly as functions of pH and applied potential U:

mO G H2O

G H2

2 eU DGH pH

mH G H2

=2 eU DGH pH

mNa G Na

eU DGH pH

9 where U is the potential measured against NHE at standard conditions and

DGH(pH) is kBT log(pH).

Gibbs free energy changes in water oxidation steps. We dened the OER steps similar to equations (25) in the manuscript. The Gibbs free energy changes for the water oxidation steps (equations (25) in the manuscript) were calculated using the

10 NATURE COMMUNICATIONS | 6:8253 | DOI: 10.1038/ncomms9253 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms9253 ARTICLE

24. Lu, Z. et al. Electrochemical tuning of layered lithium transition metal oxides for improvement of oxygen evolution reaction. Nat. Commun. 5, 4345 (2014).

25. Trotochaud, L., Ranney, J. K., Williams, K. N. & Boettcher, S. W. Solution-cast metal oxide thin lm electrocatalysts for oxygen evolution. J. Am. Chem. Soc. 134, 1725317261 (2012).

26. Kanan, M. W. et al. Structure and valency of a cobalt- phosphate water oxidation catalyst determined by in situ X-ray spectroscopy. J. Am. Chem. Soc. 132, 1369213701 (2010).

27. Lee, S. W. et al. The nature of lithium battery materials under oxygen evolution reaction conditions. J. Am. Chem. Soc. 134, 1695916962 (2012).

28. Risch, M. et al. Structural changes of cobalt-based perovskites upon water oxidation investigated by EXAFS. J. Phys. Chem. C 117, 86288635 (2013).

29. May, K. J. et al. Inuence of oxygen evolution during water oxidation on the surface of perovskite oxide catalysts. J. Phys. Chem. Lett. 3, 32643270 (2012).

30. Gerken, J. B. et al. Electrochemical water oxidation with cobalt-based electrocatalysts from pH 0-14: the thermodynamic basis for catalyst structure, stability, and activity. J. Am. Chem. Soc. 133, 1443114442 (2011).

31. Bajdich, M., Garcia-Mota, M., Vojvodic, A., Norskov, J. K. & Bell, A. T. Theoretical investigation of the activity of cobalt oxides for the electrochemical oxidation of water. J. Am. Chem. Soc. 135, 1352113530 (2013).

32. Grimaud, A. et al. Oxygen evolution activity and stability of Ba6Mn5O16, Sr4Mn2CoO9, and Sr6Co5O15: the inuence of transition metal coordination.J. Phys. Chem. C 117, 2592625932 (2013).33. Grimaud, A. et al. Double perovskites as a family of highly active catalysts for oxygen evolution in alkaline solution. Nat. Commun. 4, 2439 (2013).

34. Song, F. & Hu, X. Ultrathin cobalt-manganese layered double hydroxide is an efcient oxygen evolution catalyst. J. Am. Chem. Soc. 136, 1648116484 (2014).

35. Liu, Y. et al. Electrochemical tuning of olivine-type lithium transition-metal phosphates as efcient water oxidation catalysts. Energy Environ. Sci. 8, 17191724 (2015).

36. Wang, H. et al. Bifunctional non-noble metal oxide nanoparticle electrocatalysts through lithium-induced conversion for overall water splitting. Nat. Commun. 6, 7261 (2015).

37. Delacourt, C., Poizot, P., Tarascon, J.-M. & Masquelier, C. The existence of a temperature-driven solid solution in LixFePO4 for 0 ox o1. Nat. Mater. 4, 254260 (2005).

38. Ong, S. P., Jain, A., Hautier, G., Kang, B. & Ceder, G. Thermal stabilities of delithiated olivine MPO4 (M Fe, Mn) cathodes investigated using rst

principles calculations. Electrochem. Commun. 12, 427430 (2010).39. Masquelier, C. & Croguennec, L. Polyanionic (phosphates, silicates, sulfates) frameworks as electrode materials for rechargeable Li (or Na) batteries. Chem. Rev. 113, 65526591 (2013).

40. Kim, H. et al. Understanding the electrochemical mechanism of the new iron-based mixed-phosphate Na4Fe3(PO4)2(P2O7) in a Na rechargeable battery.

Chem. Mater. 25, 36143622 (2013).41. Erragh, F., Boukhari, A., Elouadi, B. & Holt, E. M. Crystal structures of two allotropic forms of Na2CoP2O7. J. Crystallogr. Spectrosc. Res. 21, 321326 (1991).

42. Padhi, A. K., Nanjundaswamy, K. & Goodenough, J. Phospho-olivines as positive-electrode materials for rechargeable lithium batteries. J. Electrochem. Soc. 144, 11881194 (1997).

43. Kim, H. et al. Neutron and X-ray diffraction study of pyrophosphate-based Li2x MP2O7 (M Fe, Co) for lithium rechargeable battery electrodes. Chem.

Mater. 23, 39303937 (2011).44. Shakoor, R. et al. Synthesis of NaFePO4/NaCoPO4 and their application to sodium batteries. J. Korean Battery Soc. 3, 0000 (2010).

45. Amine, K., Yasuda, H. & Yamachi, M. Olivine LiCoPO4 as 4.8 V electrode material for lithium batteries. Electrochem. Solid State Lett. 3, 178179 (2000).

46. Klingan, K. et al. Water oxidation by amorphous cobalt-based oxides: volume activity and proton transfer to electrolyte bases. ChemSusChem 7, 13011310 (2014).

47. Surendranath, Y., Kanan, M. W. & Nocera, D. G. Mechanistic studies of the oxygen evolution reaction by a cobalt-phosphate catalyst at neutral pH. J. Am. Chem. Soc. 132, 1650116509 (2010).

48. Gileadi, E. Electrode Kinetics for Chemists, Chemical Engineers, and Materials Scientists (Capstone, 1993).

49. Jin, K. et al. Partially oxidized sub-10 nm MnO nanocrystals with high activity for water oxidation catalysis. Sci. Rep. 5, 10279 (2015).

50. Gorlin, Y. et al. In situ X-ray absorption spectroscopy investigation of a bifunctional manganese oxide catalyst with high activity for electrochemical water oxidation and oxygen reduction. J. Am. Chem. Soc. 135, 85258534 (2013).

51. McAlpin, J. G. et al. EPR evidence for Co (IV) species produced during water oxidation at neutral pH. J. Am. Chem. Soc. 132, 68826883 (2010).

52. Rossmeisl, J., Qu, Z. -W., Zhu, H., Kroes, G. -J. & Nrskov, J. K. Electrolysis of water on oxide surfaces. J. Electroanal. Chem. 607, 8389 (2007).

53. Rossmeisl, J., Logadottir, A. & Nrskov, J. K. Electrolysis of water on (oxidized) metal surfaces. Chem. Phys. 319, 178184 (2005).

54. Wang, L. -P. & Van Voorhis, T. Direct-coupling O2 bond forming a pathway in cobalt oxide water oxidation catalysts. J. Phys. Chem. Lett. 2, 22002204 (2011).

55. Dogutan, D. K., McGuire, Jr. R. & Nocera, D. G. Electocatalytic water oxidation by cobalt (III) hangman b-octauoro corroles. J. Am. Chem. Soc. 133, 91789180 (2011).

56. Wang, L. -P., Wu, Q. & Van Voorhis, T. Acid-base mechanism for ruthenium water oxidation catalysts. Inorg. Chem. 49, 45434553 (2010).

57. Mavros, M. G. et al. What can density functional theory tell us about articial catalytic water splitting? Inorg. Chem. 53, 63866397 (2014).

58. Man, I. C. et al. Universality in oxygen evolution electrocatalysis on oxide surfaces. ChemCatChem 3, 11591165 (2011).

59. Liu, Y. et al. Low overpotential in vacancy-rich ultrathin CoSe2 nanosheets for water oxidation. J. Am. Chem. Soc. 136, 1567015675 (2014).

60. Kim, J., Yin, X., Tsao, K. -C., Fang, S. & Yang, H. Ca2Mn2O5 as oxygen-decient perovskite electrocatalyst for oxygen evolution reaction. J. Am. Chem. Soc. 136, 1464614649 (2014).

61. Yamazaki, Y., Babilo, P. & Haile, S. M. Defect chemistry of yttrium-doped barium zirconate: a thermodynamic analysis of water uptake. Chem. Mater. 20, 63526357 (2008).

62. Zhang, G. & Smyth, D. Protonic conduction in Ba2In2O5. Solid State Ionics 82, 153160 (1995).

63. Dunitz, J. T. & Orgel, L. Electronic properties of transition-metal oxides-II: cation distribution amongst octahedral and tetrahedral sites. J. Phys. Chem. Solids 3, 318323 (1957).

64. Kung, H. H. Transition Metal Oxides: Surface Chemistry and Catalysis (Elsevier, 1989).

65. Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 38653868 (1996).

66. Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B 47, 558561 (1993).

67. Wang, L., Maxisch, T. & Ceder, G. Oxidation energies of transition metal oxides within the GGA U framework. Phys. Rev. B. 73, 195107 (2006).

68. Khachatryan, V. et al. Search for supersymmetry in pp collisions at 7 TeV in events with jets and missing transverse energy. Phys. Lett. B 698, 196218 (2011).

69. Viswanathan, V. et al. Simulating linear sweep voltammetry from rst-principles: application to electrochemical oxidation of water on Pt (111) and Pt3Ni (111). J. Phys. Chem. C. 116, 46984704 (2012).

70. Anisimov, V. V., Zaanen, J. & Andersen, O. K. Band theory and Mott insulators: Hubbard U instead of Stoner I. Phys. Rev. B Condens. Matter 44, 943954 (1991).

Acknowledgements

This research was supported by Project Code (IBS-R006-G1) and by the Supercomputing Center/Korea Institute of Science and Technology Information with supercomputing resources including technical support (KSC-2013-C3-038), by Global Frontier R&D Program on Center for Multiscale Energy System funded by the National Research Foundation (NRF) under the Ministry of Science, ICT and Future, Korea (0420-20130104), and by Fusion Research Program for Green Technologies through the NRF funded by the Ministry of Science, ICT and Future (2012M3C1A1048863), and the National Reserch Council of Science and Technology through the Degree and Research Center Program (DRC-14-3-KBSI).

Author contributions

H.K., J.P., K.T.N. and K.K. designed the research and wrote the manuscript. H.K. and J.P. designed and performed the experiments, and anylsed the data. H.K. and J.P. contributed equally to this work. I.P. conducted the density functional theory calculations and K.J. synthesized Co-Pi catalysts. H.K., J.P., S.E.J. and S.H.K conducted EPR measurement. All authors discussed the experiments and contributed to writing the manuscript. K.T.N. and K.K. guided all aspects of the work.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Kim, H. et al. Coordination tuning of cobalt phosphates towards efcient water oxidation catalyst. Nat. Commun. 6:8253 doi: 10.1038/ncomms9253 (2015).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

NATURE COMMUNICATIONS | 6:8253 | DOI: 10.1038/ncomms9253 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 11

Word count: 10635

Show less

Abstract

Translate

The development of efficient and stable water oxidation catalysts is necessary for the realization of practically viable water-splitting systems. Although extensive studies have focused on the metal-oxide catalysts, the effect of metal coordination on the catalytic ability remains still elusive. Here we select four cobalt-based phosphate catalysts with various cobalt- and phosphate-group coordination as a platform to better understand the catalytic activity of cobalt-based materials. Although they exhibit various catalytic activities and stabilities during water oxidation, Na₂ CoP₂ O₇ with distorted cobalt tetrahedral geometry shows high activity comparable to that of amorphous cobalt phosphate under neutral conditions, along with high structural stability. First-principles calculations suggest that the surface reorganization by the pyrophosphate ligand induces a highly distorted tetrahedral geometry, where water molecules can favourably bind, resulting in a low overpotential (∼0.42 eV). Our findings emphasize the importance of local cobalt coordination in the catalysis and suggest the possible effect of polyanions on the water oxidation chemistry.

Details

Title

Coordination tuning of cobalt phosphates towards efficient water oxidation catalyst

Author

Kim, Hyunah; Park, Jimin; Park, Inchul; Jin, Kyoungsuk; Jerng, Sung Eun; Kim, Sun Hee; Nam, Ki Tae; Kang, Kisuk

Pages

8253

Publication year

2015

Publication date

Sep 2015

Publisher

Nature Publishing Group

e-ISSN

20411723

Source type

Scholarly Journal

Language of publication

English

DOI

https://doi.org/10.1038/ncomms9253

ProQuest document ID

1712333354

Coordination tuning of cobalt phosphates towards efficient water oxidation catalyst

Jump to:

Full text

Abstract

Details

Suggested sources