Monodisperse Co9S8 nanoparticles in situ embedded within N, S-codoped honeycomb-structured porous carbon for bifunctional oxygen electrocatalyst in a rechargeable Znâ€“air battery

The exploitation of cost-effective, highly active, and robust non-precious metal bifunctional oxygen electrocatalysts for both the oxygen reduction reaction (ORR) and oxygen revolution reaction (OER) is the key to promoting the application of regenerative fuel cells and metal–air batteries. Co9S8 is considered a promising non-precious metal bifunctional oxygen electrocatalyst. However, the electrocatalytic performance of cobalt sulfide-based nanocatalysts is still far from satisfactory because of their poor conductivity, insufficient exposed active sites, and aggregation-prone during continuous work. Herein, based on the inspiration from honeycombs in nature, we synthesized Co9S8/N, Scodoped honeycomb-structured porous carbon (Co9S8/NSC) in situ composites via a simple method. Benefiting from the unique honeycomb composite structure composed of monodisperse Co9S8 nanoparticles in situ embedded within the three-dimensional interconnected network carbon matrix with high conductivity, which facilitates not only the electron transport and charge transfer across the interface but also the exposure of active sites and rapid transport of ORR/OER-related species, the obtained Co9S8/NSC in situ composites exhibit high stability and activity in both the ORR in terms of more a positive half-wave potential (0.88 V vs. RHE) than that (0.86 V vs. RHE) of commercial 20% Pt/C and the OER in terms of a small overpotential (0.41 V vs. RHE) approaching that of commercial IrO2 (0.39 V vs. RHE) in alkaline electrolytes at 10 mA cm. Thus, as expected, the assembled rechargeable Zn–air batteries based on the bifunctional electrocatalysts exhibit a small discharge/charge overpotential (0.96 V) with a high voltaic efficiency of 55.1% at 10 mA cm and a long-term cycling stability of over 4000 min. Introduction With the gradual exhaustion of non-renewable resources, such as fossil fuels, and the increasing seriousness of environmental pollution, finding alternative energy sources has been a global concern. To solve this urgent problem, many attempts and painstaking experiments have been done to explore variously alternative energy storage and conversion technologies, such as solar power batteries, lithium-ion batteries, supercapacitors, fuel cells (FCs), metal–air batteries, etc. Among them, FCs and metal–air batteries, by virtue of the obvious advantages, such as low cost, resource abundant, eco-friendly, good safety, and high-energy density, are considered a promising solution for clean and sustainable energy generation. However, the indispensable catalysts, especially those with bifunctional catalytic activity for both the oxygen reduction reaction (ORR) and the oxygen © The Author(s) 2018 OpenAccessThis article is licensedunder aCreativeCommonsAttribution 4.0 International License,whichpermits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changesweremade. The images or other third partymaterial in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to thematerial. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. Correspondence: Ming-zai Wu (wumz@ustc.edu.cn) or Hai-bo Hu (haibohu@ahu.edu.cn) or Chun-yi Zhi (cy.zhi@cityu.edu.hk) School of Physics and Materials Science, Anhui University, Hefei, China Department of Physics and Materials Science, City University of Hong Kong, 83 Dachi Road, Kowloon, Hong Kong SAR, China 12 34 56 78 90 () :,; 12 34 56 78 90 () :,; 1 2 3 4 5 6 7 8 9 0 () :,; 12 34 56 78 90 () :,; revolution reaction (OER) in regenerative FCs and metalair batteries, have remained a technological bottleneck for their large-scale commercialization. Traditionally, platinum (Pt) and Pt-based materials are considered the most effective electrocatalysts for the ORR. However, they usually suffer issues of high cost, scarcity, and poor stability. On the other hand, ruthenium or iridium oxides (RuO2 or IrO2) exhibit excellent OER catalytic behavior but relatively poor ORR activity. Therefore, it is still a great challenge to prepare low-cost, resource-rich, and robust bifunctional oxygen electrocatalysts with both high OER and ORR catalytic activity for applications in highly efficient regenerative FCs and metal–air batteries to promote their practical applications. Recently, noticeable progress has been made in the development of low-cost, resource-rich, and highly efficient oxygen catalysts for the ORR and OER to replace traditional precious metals. Notably, both theoretical calculations and experimental results indicate that monodisperse Co9S8 nanocatalysts can serve as highly efficient bifunctional oxygen electrocatalysts for both the ORR and OER. However, the electrocatalytic performance of cobalt sulfide-based nanocatalysts is still far from satisfactory because they show poor conductivity, have insufficient exposed active sites and are prone to aggregation during continuous work. Therefore, determining how to effectively increase the catalytic active sites, enhance the electrical conductivity of electrocatalysts, and inhibit the aggregation of monodisperse Co9S8 nanocatalysts during continuous work is critical for the exploitation of their potential catalytic activities. Carbon materials, by virtue of their excellent conductivity and high specific surface area, are considered to be an ideal carrier for monodisperse nanocatalysts, such as Co9S8 to improve their poor electrical conductivity and inhibit aggregation. More interestingly, doping of carbon materials with highly electronegative heteroatoms, such as nitrogen can further significantly increase their conductivity by providing an advantageous channel for electron mobility. In addition, doping with nitrogen heteroatoms will provide a strong polarization in the carbon matrix, resulting in an improved binding capacity of adjacent carbon to molecular oxygen and then enhanced ORR and OER activities. In addition, the electronegativity of sulfur, close to carbon, is considered to be able to increase the exposed edge positions and disordered sites in carbon materials, such as graphene, which can enhance the electrochemical activity, stability, and alcohol resistance of carbon materials in alkaline medium. Therefore, to fully explore the potential catalytic activities of Co9S8 electrocatalysts in both the ORR and OER, a combination of monodisperse Co9S8 nanocatalysts with heteroatom-doped carbon carrier becomes appealing to improve the poor conductivity, expose more active sites, and obtain an enhanced synergistic catalytic effect. In this work, based on the inspiration from honeycombs in nature, we fabricated N, S-codoped honeycombstructured porous carbon with in situ embedded monodisperse Co9S8 nanoparticles (2–20 nm) for highly efficient bifunctional supported electrocatalysts in Zn–air batteries. The unique honeycomb-structured porous composite structure design has multiple advantages. First, the formation of a highly conductive N, S-codoped carbon matrix synchronously accompanied by the generation of monodisperse in situ embedded Co9S8 nanoparticles gives the obtained Co9S8 nanocatalysts improved conductivity and good stability during continuous work. Second, the three-dimensional (3D) interconnected porous network with high specific surface area is propitious to loading monodisperse Co9S8 nanoparticles, providing more active sites accessible to the electrolyte. Last but not the least, the 3D interconnected porous network is also beneficial for the mass transport of electrolyte ions and ORR/OERrelated species during the catalytic process, allowing a rapid and sustained reaction. As expected, benefiting from the reasonable structural design and appropriate materials selection, the as-fabricated Co9S8/N, S-codoped honeycomb-structured porous carbon (Co9S8/NSC) in situ composites exhibit not only excellent ORR catalytic activity, of which the half-wave (0.88 V vs. RHE) potentials clearly outperformed that of the commercial 20% Pt/C (0.86 V vs. RHE) but also good catalytic activity and stability for the OER in terms of a small overpotential (0.41 V vs. RHE at 10mA cm) approaching that of commercial IrO2 (0.39 V vs. RHE) in alkaline electrolytes. Thus, as a proof-of-concept application, the further assembled rechargeable Zn–air batteries based on the bifunctional electrocatalysts exhibit a small discharge/charge overpotential (0.96 V) with a high voltaic efficiency of 55.1% at 10mA cm and a long-term cycling stability over 4000 min. Experimental procedures Synthesis of SiO2 nanospheres SiO2 nanoparticles were prepared by the Stöber method. In a typical procedure, 4.5 mL of tetraethyl orthosilicate was added to a mixed solution containing 65.0 mL of ethanol and 25.0 mL of deionized (DI) water under vigorous stirring. Then, 9 mL of aqueous ammonia solution was added dropwise to the above mixture. The mixture was stirred for 8 h at room temperature. Subsequently, the product was centrifuged and washed several times with ethanol and DI water. Finally, the white SiO2 nanospheres with an average size of ∼200 nm were obtained after drying overnight at 60 °C. Cao et al. NPG Asia Materials (2018) 10: 670-684 671


Introduction
With the gradual exhaustion of non-renewable resources, such as fossil fuels, and the increasing seriousness of environmental pollution, finding alternative energy sources has been a global concern. To solve this urgent problem, many attempts and painstaking experiments have been done to explore variously alternative energy storage and conversion technologies, such as solar power batteries, lithium-ion batteries, supercapacitors, fuel cells (FCs), metal-air batteries, etc. Among them, FCs and metal-air batteries, by virtue of the obvious advantages, such as low cost, resource abundant, eco-friendly, good safety, and high-energy density, are considered a promising solution for clean and sustainable energy generation [1][2][3][4] . However, the indispensable catalysts, especially those with bifunctional catalytic activity for both the oxygen reduction reaction (ORR) and the oxygen revolution reaction (OER) in regenerative FCs and metalair batteries, have remained a technological bottleneck for their large-scale commercialization [5][6][7] . Traditionally, platinum (Pt) and Pt-based materials are considered the most effective electrocatalysts for the ORR. However, they usually suffer issues of high cost, scarcity, and poor stability [8][9][10][11] . On the other hand, ruthenium or iridium oxides (RuO 2 or IrO 2 ) exhibit excellent OER catalytic behavior but relatively poor ORR activity 12,13 . Therefore, it is still a great challenge to prepare low-cost, resource-rich, and robust bifunctional oxygen electrocatalysts with both high OER and ORR catalytic activity for applications in highly efficient regenerative FCs and metal-air batteries to promote their practical applications.
Recently, noticeable progress has been made in the development of low-cost, resource-rich, and highly efficient oxygen catalysts for the ORR and OER to replace traditional precious metals [14][15][16][17][18][19][20] . Notably, both theoretical calculations and experimental results indicate that monodisperse Co 9 S 8 nanocatalysts can serve as highly efficient bifunctional oxygen electrocatalysts for both the ORR and OER [21][22][23] . However, the electrocatalytic performance of cobalt sulfide-based nanocatalysts is still far from satisfactory because they show poor conductivity, have insufficient exposed active sites and are prone to aggregation during continuous work [24][25][26] . Therefore, determining how to effectively increase the catalytic active sites, enhance the electrical conductivity of electrocatalysts, and inhibit the aggregation of monodisperse Co 9 S 8 nanocatalysts during continuous work is critical for the exploitation of their potential catalytic activities. Carbon materials, by virtue of their excellent conductivity and high specific surface area, are considered to be an ideal carrier for monodisperse nanocatalysts, such as Co 9 S 8 to improve their poor electrical conductivity and inhibit aggregation. More interestingly, doping of carbon materials with highly electronegative heteroatoms, such as nitrogen can further significantly increase their conductivity by providing an advantageous channel for electron mobility [27][28][29][30] . In addition, doping with nitrogen heteroatoms will provide a strong polarization in the carbon matrix, resulting in an improved binding capacity of adjacent carbon to molecular oxygen and then enhanced ORR and OER activities [31][32][33] . In addition, the electronegativity of sulfur, close to carbon, is considered to be able to increase the exposed edge positions and disordered sites in carbon materials, such as graphene, which can enhance the electrochemical activity, stability, and alcohol resistance of carbon materials in alkaline medium [34][35][36][37][38] . Therefore, to fully explore the potential catalytic activities of Co 9 S 8 electrocatalysts in both the ORR and OER, a combination of monodisperse Co 9 S 8 nanocatalysts with heteroatom-doped carbon carrier becomes appealing to improve the poor conductivity, expose more active sites, and obtain an enhanced synergistic catalytic effect.
In this work, based on the inspiration from honeycombs in nature, we fabricated N, S-codoped honeycombstructured porous carbon with in situ embedded monodisperse Co 9 S 8 nanoparticles (2-20 nm) for highly efficient bifunctional supported electrocatalysts in Zn-air batteries. The unique honeycomb-structured porous composite structure design has multiple advantages. First, the formation of a highly conductive N, S-codoped carbon matrix synchronously accompanied by the generation of monodisperse in situ embedded Co 9 S 8 nanoparticles gives the obtained Co 9 S 8 nanocatalysts improved conductivity and good stability during continuous work. Second, the three-dimensional (3D) interconnected porous network with high specific surface area is propitious to loading monodisperse Co 9 S 8 nanoparticles, providing more active sites accessible to the electrolyte. Last but not the least, the 3D interconnected porous network is also beneficial for the mass transport of electrolyte ions and ORR/OERrelated species during the catalytic process, allowing a rapid and sustained reaction. As expected, benefiting from the reasonable structural design and appropriate materials selection, the as-fabricated Co 9 S 8 /N, S-codoped honeycomb-structured porous carbon (Co 9 S 8 /NSC) in situ composites exhibit not only excellent ORR catalytic activity, of which the half-wave (0.88 V vs. RHE) potentials clearly outperformed that of the commercial 20% Pt/C (0.86 V vs. RHE) but also good catalytic activity and stability for the OER in terms of a small overpotential (0.41 V vs. RHE at 10 mA cm −2 ) approaching that of commercial IrO 2 (0.39 V vs. RHE) in alkaline electrolytes. Thus, as a proof-of-concept application, the further assembled rechargeable Zn-air batteries based on the bifunctional electrocatalysts exhibit a small discharge/charge overpotential (0.96 V) with a high voltaic efficiency of 55.1% at 10 mA cm −2 and a long-term cycling stability over 4000 min.

Synthesis of SiO 2 nanospheres
SiO 2 nanoparticles were prepared by the Stöber method 39 . In a typical procedure, 4.5 mL of tetraethyl orthosilicate was added to a mixed solution containing 65.0 mL of ethanol and 25.0 mL of deionized (DI) water under vigorous stirring. Then, 9 mL of aqueous ammonia solution was added dropwise to the above mixture. The mixture was stirred for 8 h at room temperature. Subsequently, the product was centrifuged and washed several times with ethanol and DI water. Finally, the white SiO 2 nanospheres with an average size of ∼200 nm were obtained after drying overnight at 60°C.
Preparation of honeycomb-structured porous Co 9 S 8 /NSC in situ composite electrocatalyst First, 0.1 g of Co(NO 3 ) 2 ·6H 2 O was dissolved in 50 mL of alcohol solution (37.5% ethanol) to form a red solution. Then, 1.0 g of SiO 2 nanospheres, 1.0 g of thiourea (CS (NH 2 ) 2 ), and 0.5 g of glucose were successively added to the above solution under continuous stirring at room temperature for 2 h. Then, the as-obtained mud-like precursor was placed in tube furnace and heated to 900°C for 2 h in a N 2 atmosphere. After further natural cooling, the sample was immersed in 2 M NaOH solution at 50°C for 8 h to remove the SiO 2 template, and then centrifuged and washed with DI water to remove residual NaOH until reaching pH = 7. Finally, the honeycomb-structured porous in situ composites, i.e., Co 9 S 8 /NSC-900 porous composites, were obtained after drying overnight at 60°C. For comparison, a series of catalysts labeled C-900, NSC-900, Co 9 S 8 /NSC-800, and Co 9 S 8 /NSC-1000 were also prepared using similar procedures, but without (CS(NH 2 ) 2 ), without Co(NO 3 ) 2 ·6H 2 O, and at 800 and 1000°C, respectively.

Characterization
The microstructure and phase composition of the as obtained samples were characterized by X-ray powder diffraction (XRD Bruker D8-ADVANCE) with an 18 kW advanced X-ray diffractometer with Cu K α radiation (λ = 1.54056 Å), Raman spectroscopy (inVia-Reflex, Renishaw, UK), an elemental analyzer (EA, Perkin-Elmer 2400 CHN), field-emission scanning electron microscopy (FE-SEM, S-4800, Hitachi, Japan), transmission electron microscopy (TEM, JEM-2100, JEOL, Japan), and Fourier transform infrared spectroscopy (VERTEX 801HYPER-ION2000, Bruker Optics, Germany). The surface area was determined by a nitrogen adsorption-desorption method at 77 K (3H-2000PS2, BeiShiDe Instrument, China). X-ray photoelectron spectroscopy (XPS) was conducted with a Mg Kα achromatic X-ray source. Thermo-gravimetric analysis (TGA) measurements were performed on Diamond TG 6300 at a heating rate of 10 C min −1 in an air atmosphere in the range of 40-900°C. The surface areas and pore sizes of the as-prepared samples were analyzed by the nitrogen adsorption-desorption isotherms (BET, Quantachrome, NOVAe). Inductively coupled plasmaatomic emission spectrometry (ICP-AES) measurements were conducted to determine the cobalt and sulfur concentrations of the samples with an Optima 7400 DV instrument.

Electrochemical measurements
All the ORR tests were performed on the electrochemical work station (CHI 660E, Chenhua, Shang-hai) in a conventional three-electrode cell equipped with a Pt wire counter electrode, an Ag/AgCl reference electrode and a glassy carbon (GC) disk electrode. The measured potentials are referred to the RHE, E RHE = E Ag/AgCl + 0.059 pH + E θ Ag/AgCl . The rotating disk electrode (RDE) was used as the working electrode. For each sample, 5 mg of catalysts were ultrasonically dispersed in a suspension containing 700 μL of DI water, 300 μL of isopropanol and 50 μL of 5 wt% Nafion solution. Subsequently, the obtained homogeneous catalyst ink with a volume of 8 μL was dripped on a GC RDE (5 mm in diameter) and dried thoroughly at room temperature. The Pt/C catalyst (20 wt %) ink was prepared as described above for comparison. In alkaline media, the RDE test was carried out at room temperature in O 2 -saturated 0.1 M KOH solution. The cyclic voltammetry (CV) measurements were scanned from −0.9 to 0.1 V (vs. Ag/AgCl) at a scan rate of 20 mV s −1 . The RDE test was carried out by linear sweep voltammetry (LSV) at a potential range of −0.9 to 0.2 V (vs. Ag/AgCl) at a scanning rate of 5 mV s −1 at different rotational rates of O 2 -saturation at 400-2500 rpm alkaline electrolyte. Under the same conditions, OER measurements were carried out by LSV at a potential range from 0 to 1.0 V (vs. Ag/AgCl) at a scanning rate of 5 mV s −1 under 1600 rpm. The IrO 2 catalyst ink was prepared as described above for comparison. Koutecky-Levich (K-L) plots were obtained from the ORR polarization curves at different potentials. The electron transfer number (n) per O 2 molecule in the ORR was determined by the slopes of the linear fit lines on the basis of the following K-L equation: where j is the measured current density, j k is the kinetic current density, ω is the electrode rotation rate, F is the Faradic constant (96,485 C mol −1 ), C O2 is the saturated oxygen concentration in 0.1 M KOH aqueous solution (1.2 × 10 −6 mol cm −3 ), D O2 is the oxygen diffusion coefficient (1.73 × 10 −5 cm 2 s −1 ), v is the kinematic viscosity of the solution (0.01 cm 2 s −1 ) and n is the electron transfer number in the ORR. Rotating ring-disk electrode (RRDE) voltammogram measurements for the calculation of the ORR electron transfer number were conducted on an RRDE configuration with a 320 µm gap Pt ring electrode. The electron transfer numbers (n) and peroxide yield (HO 2 − %) can be obtained from the RRDE curve using Eqs. (3) and (4), respectively. 40 where I d is the disk current, I r is the ring current and N is the current collection efficiency of the Pt ring which was determined to be 0.37.

Fabrication and test of Zn-air batteries
The Zn-air battery test was conducted in a homemade battery, with a Zn plate serving as the anode and catalystcoated carbon cloth (1 cm 2 ) serving as the air electrode. Typically, 5 mg of catalysts were ultrasonically dispersed into 1 mL of ethanol and 50 μL of Nafion solution by sonication for 30 min. Then, some amount of the catalyst ink was dropped on the carbon paper to obtain a catalyst loading of 1 mg cm −2 . The electrolyte was a 6 M KOH solution containing 0.2 M Zn(Ac) 2 . Battery tests were carried out an electrochemical workstation at room temperature.

Characterizations of Co 9 S 8 /NSC in situ composites
The detailed preparation process for the Co 9 S 8 /NSC in situ composites is schematically illustrated in Fig. 1. The easily manipulated synthesis protocol first began with mixing SiO 2 nanospheres, thiourea, and glucose with Co (NO 3 ) 2 ·6H 2 O solution. Then, the obtained precursor was thermally treated, during which glucose was pyrolyzed into the N, S-codoped carbon matrix with the in situ formed monodisperse Co 9 S 8 nanoparticles (2-20 nm) and encapsulated SiO 2 nanospheres. Finally, the encapsulated SiO 2 nanospheres templates were removed by being dissolved in 2 M NaOH solution and washed away by DI water. Thus, the Co 9 S 8 /NSC in situ composites were obtained. It should be noted that the easily manipulated synthesis method does not involve tedious operations and any expensive equipment or reagents, which is beneficial for the scale-production of the low-cost catalysts and their commercial application.
X-ray diffraction (XRD) was first used to determine the crystalline phase of the synthesized samples. Figure 2a shows the XRD patterns of the Co 9 S 8 /NSC-900 product, and the contrast samples C-900 and NSC-900. As shown, the diffraction peaks at 2θ = 15.8°, 29.8°, 31.2°, 39.5°, 47.6°, and 51.9°are ascribable to the (111), (311), (222), (331), (511), and (440) planes, respectively, of Co 9 S 8 (JCPDS card 65-6801) 23,41,42 , directly indicating the successful of formation of Co 9 S 8 within the porous carbon skeleton. In addition, the energy-dispersive spectrum (EDS) (Figure S1a) also shows the contained Co, N, S, C, and O elements in the Co 9 S 8 /NSC-900 product. A representative Raman spectrum of the as-obtained Co 9 S 8 /NSC-900 is further shown in Fig. 2b, which can be divided into two parts: <800 and >800 cm −1 . The Raman bands below 800 cm −1 can be indexed to Co 9 S 8 , and those above 800 cm −1 are related to amorphous carbon in the Co 9 S 8 /NSC-900 composite. Typical D and G bands of amorphous carbon located at 1369 and 1583 cm −1 were clearly observed in this spectrum, indicating the formation of carbonaceous composites after pyrolysis. The relative intensity ratio of the D band to the G band (I D /I G ) has been widely adopted for evaluating the graphitization degree of carbon 43,44 . The I D /I G ratios of NSC-900 (0.95) and Co 9 S 8 /NSC-900 (0.94) are lower than that of C-900 (1.03), which can be attributed to the N and S heteroatom doping. The results show that N and S heteroatom doping causes many defects in the internal portion of honeycomb-like porous carbon skeleton, which can further enhance the degree of graphite, resulted in improved catalyst conductivity and catalytic activity.
To investigate the porosity of C-900, NSC-900, and Co 9 S 8 /NSC-900 samples, nitrogen adsorption-desorption isotherms were further measured. As shown in Fig. 2c, all the samples display a type-IV curve with a distinct hysteresis loop, indicating the presence of many mesopores and macropores in these samples 45 . The specific surface area of the Co 9 S 8 /NSC-900 is determined to be 556 m 2 g −1 , which is clearly larger than that of the C-900 (417 m 2 g −1 ) and NSC-900 (490 m 2 g −1 ). According to the Barrett −Joyner−Halenda (BJH) method, these samples mainly contain mesopores with diameters between 2 and 20 nm (Fig. 2d). The nonlocal density functional theory (NLDFT) method was further employed to analyze the pore size distribution of Co 9 S 8 /NSC-900 in situ composites more rigorously. As clearly shown in Figure S1b, in addition to the mesopores with diameters between 2 and 20 nm, the Co 9 S 8 /NSC-900 in situ composites also contain a large number of micropores with diameters between 0.89 and 2.0 nm. We think that these micropores and mesopores are mainly derived from the combined action of the Fig. 1 Flowchart of the easily-manipulated synthesis protocol for the Co 9 S 8 /N, S-codoped honeycomb-structured porous carbon in-situ composites released gases in the annealing process and following alkaline etching process, resulting in a large specific area of the as-prepared the Co 9 S 8 /NSC-900 in situ composites with a honeycomb-like porous structure. This is of great significance for catalytic applications, as it can effectively facilitate the exposure of active sites and rapid transport of ORR/OER-related species, allowing a rapid and sustained reaction. Figure 3a shows a typical SEM image of the SiO 2 nanospheres used in the synthesis with a uniform average diameter of ∼200 nm. By taking advantage of the SiO 2 nanosphere template, the as-obtained NSC-900, Co 9 S 8 / NSC-800, Co 9 S 8 /NSC-900, and Co 9 S 8 /NSC-1000 products have a similar honeycomb-like porous structure, as clearly demonstrated in the SEM image (Fig. 3b-d and Figure S1c-d). The typical optical photo of a natural honeycomb in the inset of Fig. 3b clearly shows a high degree of structural similarity to the as-obtained samples. As further shown in the enlarged SEM image (Fig. 3e), the prepared honeycomb-structured porous Co 9 S 8 /NSC-900 composites consist of a number of interconnected macropores composed of carbon film with high conductivity, whose surface is uniformly covered with in situ generated and embedded monodisperse Co 9 S 8 nanoparticles with diameters from 2 to 20 nm. The joint action of the dissolution of SiO 2 , the surface tension, as well as the mild chemical reaction promotes the formation of the honeycomb-like structure. We further carried out EDSmapping analysis of the Co 9 S 8 /NSC-900 sample. As shown in Fig. 3f-i, the uniform dispersion of C, N, S, and Co elements across the entire image directly indicates the formation of Co 9 S 8 /NSC in situ composites.
The morphology and microstructure of the Co 9 S 8 /NSC-900 in situ composites were further characterized with TEM and high-resolution TEM (HRTEM). As clearly shown in Fig. 4a, b, the Co 9 S 8 /NSC-900 in situ composites have a unique honeycomb-like porous hybrid structure composed of a large number of in situ generated monodisperse Co 9 S 8 nanoparticles uniformly dispersed on the carbon skeleton with high conductivity. Figure 4c, d are the HRTEM images of the monodisperse Co 9 S 8 nanoparticles with obvious lattice edges. The measured lattice distance of 0.29 nm corresponds to the (311) plane of Co 9 S 8 . Notably, as shown in Fig. 4c, d, the in situ generated Co 9 S 8 nanoparticles (2-20 nm) were tightly encapsulated by the carbon film. The unique in situ embedding Fig. 2 Characterization of composition and pore size distribution of the samples, a XRD patterns of C-900, NSC-900 and Co 9 S 8 /NSC-900; b Raman spectra of C-900, NSC-900 and Co 9 S 8 /NSC-900; c N2 adsorption-desorption isotherms, and d BJH pore-size distribution of C-900, NSC-900 and Co 9 S 8 / NSC-900 configuration would greatly reduce the tendency of agglomeration of monodisperse Co 9 S 8 nanoparticles and simultaneously improve their poor conductivity, resulting in enhanced ORR/OER catalytic activity. The electron diffraction rings shown in the inset of Fig. 4d further indicate the porous honeycomb structure of the Co 9 S 8 / NSC-900. Two typical partial diffraction rings were indexed in the SAED pattern and can be assigned to the (311) and (440) planes of the Co 9 S 8 phase. All the above results fully proved the successful fabrication of honeycomb-structured porous Co 9 S 8 /NSC-900 in situ composites, which could be expected to exhibit excellent ORR/OER catalytic activity depending on the produced graphitic carbon skeleton with high conductivity, codoped pyridinic-N and thiophene-S heteroatoms, and the in situ embedded monodisperse Co 9 S 8 nanoparticles.
To shed light on the chemical compositions and elemental valence states of the Co 9 S 8 /NSC-900 in situ composites, XPS, elemental analysis (EA), thermogravimetric analysis (TGA), and ICP-AES measurements were further carried out. Figure 5a is the XPS survey spectrum, which indicates the presence of C, N, S, and Co. The resulting estimated atomic weight percentages are 84.29% for C, 8.27% for N, 3.59% for S, and 3.85% for Co, respectively, which reveals that the atomic ratio of S:Co is approximately 1:1.07, consistent with the atomic ratio of S:Co in the chemical molecular formula of Co 9 S 8 . In addition, the atomic ratio is also congruent with that calculated from the EDS-mapping analysis (Table S1). These results forcefully prove the formation of cobalt sulfides in the composite. A typical TGA curve ( Figure S2a) illustrates the weight loss of the as-prepared Co 9 S 8 /NSC-900 sample with the temperature increase from 40 to 900°C in air. The sample has an initial loss of weight (∼18 wt%) below 400°C, which can be ascribed to the removal of absorbed water and air because of the large specific surface area of the as-prepared Co 9 S 8 /NSC-900 sample. However, there is a significant weight drop in the temperature range from 400 to 500°C, mainly attributed to the loss of carbon. The final remnant product (17.08 wt%) is mainly Co 3 O 4 transformed from Co 9 S 8 during the thermogravimetric process in the air, which can be proven from the obtained XRD pattern of the final remnant product, as shown in Figure S2b. It is worth noting that no diffraction peaks belonging to silicon dioxide are observed in the XRD pattern, which indicates that the as-prepared Co 9 S 8 /NSC-900 sample has a very low level of unremoved silicon dioxide. According to the EA and ICP-AES results for Co 9 S 8 /NSC-900, the estimated weight percentages of C, N, S, and Co, are approximately 53.87, 5.74, 7.44, and 15.03 wt%, respectively (Table S1). The weight percentage of the missing element is also consistent with the results from TGA. Detailed information about the bonding configuration of each element was further given by high-resolution XPS spectra; the schematic of different types of nitrogen and sulfur is shown in Fig. 5b. The high-resolution spectrum of C 1s (Fig. 5c) can be deconvoluted into several peaks, corresponding to C = C-C (284.5 eV), C-S-C (285.4 eV), C-N-C (286.5 eV), and C = O (288.4 eV), indicating that N and S were successfully doped into the carbon skeleton 29,46,47 . For the high-resolution N 1s XPS spectrum (Fig. 5d), the N 1s peak can also be deconvoluted into three peaks, corresponding to pyridinic N (398.6 eV), pyrrolic N (400.2 eV), and graphitic N (401.4 eV), indicating N doping in the graphitic carbon structure [48][49][50] . Figure 5e shows the high-resolution S 2p XPS spectrum for the Co 9 S 8 /NSC-900 in situ composites. In the S 2p XPS spectrum, the raw curve can be peak-fitted into four curves attributed to three different types of sulfur species. The peaks at 161.6 and 163.5 eV can be assigned to the S 2p 3/2 and S 2p 1/2 binding energies of S 2− in Co 9 S 8 , respectively, and the peaks at 164.6 and 168.1 eV are attributed to C-S and SO 4 2− , respectively. The presence of C-S may originate from the interfaces between Co 9 S 8 and the carbon matrix, whereas SO 4 2− arose from partly oxidized sulfur species or residual sulfate groups on the surface of the material. In the Co 2p spectrum of the Co 9 S 8 /NSC-900 (Fig. 5f) in situ composites, the peaks at 796.68 and 798.54 eV are, respectively, assigned to the 2p 1/2 of Co 2+ and Co 3+ , and the peaks at 781.03 and 783.18 eV, respectively, belong to the 2p 3/2 of Co 2+ and Co 3+51,52 , as well as the corresponding two shakeup satellites peaks (abbreviated as "Sat."), also suggesting the formation of cobalt sulfides in the carbon structure.
As is well-known, temperature is an important parameter that affects the structure and composition of the prepared catalyst, thus determining the ORR activities. Therefore, the Co 9 S 8 /NSC catalysts prepared at lower temperature (800°C) and higher temperature (1000°C) were further studied to determine the best heat treatment conditions. The crystalline structures of the as-obtained Co 9 S 8 /NSC-800 and Co 9 S 8 /NSC-1000 were first characterized by X-ray diffraction (XRD). As Figure S2c shows, all the corresponding XRD diffraction patterns confirm the successful formation of Co 9 S 8 in the 3D porous carbon skeleton at both lower or higher pyrolysis temperatures. However, it is worth noting that, as the pyrolysis temperature increases, the peak intensity of N 1s and S 2p gradually reduced, indicating that the contents of doped N and S decreased, as shown in Figure S2d. Even so, as further clearly shown in the high-resolution N 1s and S 2p spectra ( Figure S3a-h), the relative percentages of pyrrolopyridine-N and pyridine-N decreased with the increase in pyrolysis temperature, while the relative percentage of graphite-N and thiophene-S increased and the intensity of Co-S bonds became gradually higher. These results indicate that higher pyrolysis temperature favors the conversion of doped-N from pyrrolopyridine-N/pyridinic-N to graphitic-N and doped-S from SO x to thiophene-S. The graphitic-N groups have been demonstrated to be the active species for the ORR/OER [53][54][55][56][57] . In addition, thiophene-S is considered to be able to increase the exposed edge positions and disordered sites in carbon materials, such as graphene, which can enhance the electrochemical activity, stability, and alcohol resistance of carbon materials in alkaline media [34][35][36][37][38] . Electrochemical impedance spectrum (EIS) measurements were further carried out to investigate the potential impact of pyrolysis temperature on the catalytic performance of the ; c-f high-resolution spectra of C 1s, N 1s S 2p and Co 2p for Co 9 S 8 /NSC-900 in-situ composites, respectively as-obtained samples. As shown in Figure S3i, with the increase in the pyrolysis temperature, the as-obtained samples tend to show a lower charge-transfer resistance (Rct, 22.3 Ω for Co 9 S 8 /NSC-800, 8.4 for Co 9 S 8 /NSC-900, and 14.2 Ω for Co 9 S 8 /NSC-1000), which can be attributed to the enhanced conductivity of the 3D porous carbon matrix with increased graphitic-N and thiophene-S content, facilitating the electron transport across the interfaces between the Co 9 S 8 , carbon skeleton, and electrolyte. This is also very beneficial for the improvement of the catalytic activity of the prepared catalyst. Thus, based on the above characterization results, under the combined action of the higher proportion of graphite-N, thiophene-S, the in situ embedded active Co 9 S 8 nanocatalysts, as well as the unique honeycomb porous carbon skeleton with enhanced conductivity, the Co 9 S 8 /NSC-900 composite catalysts could be expected to exhibit better catalytic activity and electrochemical kinetics compared to those of the other two samples.

Electrochemical performance for ORR and OER
The electrochemical ORR activity properties of the Co 9 S 8 /NSC-900 in situ composites and all the additional samples in O 2 or N 2 -saturated 0.1 M KOH were first evaluated through CV measurements (Fig. 6a). Generally, higher peak potential and larger peak current mean higher catalytic activity. It is obvious that both the peak potential and peak current of honeycomb-like Co 9 S 8 /NSC-900 in situ composites are larger than those of C-900 (peak potential at 0.78 V vs. RHE) and the NSC-900 (peak potential at 0.84 V vs. RHE), indicating better ORR activity with a higher peak potential at 0.88 V (vs. RHE). To gain better insight into the kinetics of the ORR catalyzed by Co 9 S 8 /NSC-900, linear sweep voltammograms (LSVs) were also collected in O 2 -saturated 0.1 M KOH at an electrode rotation rate of 1600 rpm. Figure 6b and Figure S4 show that the Co 9 S 8 /NSC-900 in situ composites show an onset potential of 0.96 V (vs. RHE), approaching the value of the commercial 20% Pt/C (1.0 V vs. RHE). Meanwhile, the Co 9 S 8 /NSC-900 in situ composites exhibit a more positive half-wave potential (0.88 V vs. RHE) than that of C-900 (0.78 V vs. RHE), NSC-900 (0.84 V vs. RHE), Co 9 S 8 /NSC-800 (0.85 V vs. RHE), Co 9 S 8 /NSC-1000 (0.85 V vs. RHE), and the commercial 20% Pt/C (0.86 V vs. RHE), respectively. The current density (at 0.2 V vs. RHE) of the Co 9 S 8 /NSC-900 in situ composites (5.56 mA cm −2 ) was also higher than those found for C-900 (3.45 mA cm −2 ), NSC-900 (3.84 mA cm −2 ), Co 9 S 8 /NSC-800 (4.75 mA cm −2 ), Co 9 S 8 /NSC-1000 (5.18 mA cm −2 ), and the commercial 20% Pt/C (5.51 mA cm −2 V). To further evaluate the structural and compositional advantages of Co 9 S 8 /NSC-900, the Tafel plots of the Co 9 S 8 /NSC-900 and contrastive samples, derived from the ORR polarization curves, are shown in Fig. 6c. The corresponding Tafel slopes of Co 9 S 8 /NSC-900, Co 9 S 8 /NSC-1000, Co 9 S 8 /NSC-800, and 20% Pt/C were calculated to be 64, 66, 69, and 72 mV dec −1 , respectively. The smaller Tafel slope indicates a superior ORR activity of Co 9 S 8 /NSC-900. All of the above results fully indicate the better ORR activity of honeycomb-like Co 9 S 8 /NSC-900 in situ composites compared to that of commercial 20% Pt/C and the additional samples in the study.
LSVs were recorded at different rotation rates (Fig. 6d) to further investigate the ORR kinetics of the Co 9 S 8 /NSC-900 in situ composites. The K-L plots (Fig. 6d, inset) show a good linear relationship for the Co 9 S 8 /NSC-900 in situ composites and the measured electron transfer number is close to 4, suggesting that the ORR on the Co 9 S 8 /NSC-900 catalyst mainly follows the effective fourelectron pathway, in which the oxygen molecule (O 2 ) is reduced to OH − directly. In addition, to quantify the yield of the peroxide intermediate, as well as the electron transfer number, RRDE measurements with superior precision were further carried out (Fig. 6e). Based on the ring current and disk current, the average electrontransfer number (n) of Co 9 S 8 /NSC-900 is calculated to be 3.98 in the potential range from 0.8 to 0.2 V with a low H 2 O 2 yield of 3.4% (Fig. 6e, inset), which is close to those of 20% Pt/C (H 2 O 2 yield of 2.96%, n of 3.99, Figure S5a and b). This effectively indicates that the ORR process under the action of the Co 9 S 8 /NSC-900 catalyst was almost the four-electron pathway in the alkaline electrolyte, which can be depicted as follows: The durability and tolerance to the methanol poisoning effect of the Co 9 S 8 /NSC-900 in situ composites and commercial 20% Pt/C in O 2 -saturated 0.1 M KOH solution ( Figure S5c-f) are also evaluated through the LSV curves, respectively. Notably, the LSV curve of the Co 9 S 8 / NSC-900 in situ composites did not change significantly ( Figure S5c) when 0.5 M CH 3 OH was added to the O 2saturated 0.1 M KOH solution. In sharp contrast, the commercial 20% Pt/C catalyst underwent tremendous changes ( Figure S5d). The results suggested that the Co 9 S 8 /NSC-900 catalyst has superior resistance to methanol poisoning than the commercial 20% Pt/C catalyst. In addition, the respective cycling stability of the Co 9 S 8 /NSC-900 in situ composites and the commercial 20% Pt/C was also evaluated by cycling the catalyst between 0.2 and 1.1 V (vs. RHE) at a high scan rate of 300 mV s −1 in O 2 -saturated 0.1 M KOH. The current density observed in the LSV curve of the Co 9 S 8 /NSC-900 recorded after the 10,000 cycles decreased a little compared to that observed in the initial LSV curve ( Figure S5e). In comparison, the commercial 20% Pt/C exhibited an obvious current density decay only after 2000 cycles  Figure S5f). These results directly indicated the superior stability of the Co 9 S 8 /NSC-900 catalyst in comparison with the commercial 20% Pt/C.
The OER catalytic activity is another very important performance parameter for the bifunctional oxygen electrocatalysts used in rechargeable Zn-air batteries. Therefore, the LSVs for the OER of the Co 9 S 8 /NSC-900 in situ composite catalysts were also measured in O 2saturated 0.1 M KOH aqueous solution at 1600 rpm. Figure 6f shows that Co 9 S 8 /NSC-900 in situ composites exhibit a more positive onset potential (1.51 V vs. RHE) than that of C-900 (1.65 V vs. RHE), NSC-900 (1.64 V vs. RHE), Co 9 S 8 /NSC-800 (1.56 V vs. RHE), and Co 9 S 8 /NSC-1000 (1.54 V vs. RHE). Meanwhile, the Co 9 S 8 /NSC-900 in situ composites show a smaller overpotential (0.41 V vs. RHE) compared to that of Co 9 S 8 /NSC-800 (0.53 V vs. RHE), and Co 9 S 8 /NSC-1000 (0.47 V vs. RHE) at a current density of 10 mA cm −2 , approaching the value of commercial IrO 2 (0. 39 V vs. RHE) and indicated the good OER activity of Co 9 S 8 /NSC-900. The Tafel plots derived from LSV polarization curves further verify the good OER performance of Co 9 S 8 /NSC-900, as shown in Fig. 6g. The corresponding Tafel slope values for Co 9 S 8 /NSC-900, Co 9 S 8 /NSC-800, Co 9 S 8 /NSC-1000, and IrO 2 were calculated to be 99, 103, 105, and 94 mV dec −1 , respectively. Thus, apparently, Co 9 S 8 /NSC-900 also exhibits a close Tafel slope compared with that of the IrO 2 catalyst, also indicating the good OER activity of Co 9 S 8 /NSC-900. The OER stability of the Co 9 S 8 /NSC-900 in situ composites was further evaluated through 10,000 CV cycles at a high scan rate of 300 mV s −1 in O 2 -saturated 0.1 M KOH between 1.1 and 1.8 V (vs. RHE). The observed overpotential in the polarization curve recorded after 10,000 cycles only shows a small decline of 27 mV ( Figure S5g). In comparison, the commercial IrO 2 catalyst exhibited a larger attenuation of 41 mV after 10,000 cycles ( Figure S5h). These results directly indicate the superior stability of the Co 9 S 8 /NSC-900 catalyst in the OER in comparison with the commercial IrO 2 catalyst.
(see figure on previous page) Fig. 6 Electrochemical performance characterization for ORR and OER, a CVs of C-900, NSC-900 and Co 9 S 8 /NSC-900 in-situ composites in N 2 and O 2saturated 0.1 M KOH. Scan rate: 20 mV s-1; b Steadystate LSVs on different catalysts at a rotating rate of 1600 rpm in O 2 -saturated 0.1 M KOH; c Tafel plot of Co 9 S 8 /NSC-800, Co 9 S 8 /NSC-900, Co 9 S 8 /NSC-1000, and 20% Pt/C; d LSVs for the ORR/K-L plots of Co 9 S 8 /NSC-900 in-situ composites at different rotating rates; e peroxide yield and electron transfer numbers for Co 9 S 8 /NSC-900 obtained from RRDE voltammograms at a rotating rate of 1600 rpm in O 2 -saturated 0.1 M KOH; f LSVs for the OER of Co 9 S 8 /NSC-900 in-situ composites in 0.1 M KOH at 1600 rpm for different catalysts (scan rate, 5 mV s-1; electrolyte in 0.1 M KOH; rotation speed, 1600 rpm); g Tafel plot of Co 9 S 8 /NSC-800, Co 9 S 8 /NSC-900, Co 9 S 8 /NSC-1000, and IrO 2 ; h the overall LSV curves of all samples in the whole ORR and OER region in 0.1 M KOH The bold value means the best performance in this item compared to the previously reported works list in the table.
To intuitively grasp the catalytic performance of bifunctional honeycomb-like Co 9 S 8 /NSC-900 in situ composites, the key parameters for evaluating the ORR/ OER catalytic activity were compared to those of other currently reported state-of-the-art catalysts ( Table 1). As shown Fig. 6h, it is obvious that the Co 9 S 8 /NSC-900 in situ composites exhibit excellent ORR/OER bifunctional catalytic activity with a small potential difference (ΔE = 0.76 V) between the half-wave potential of the ORR and the potential of the OER at 10 mA cm −2 , which is usually used to estimate the performance of a bifunctional oxygen electrode catalyst [58][59][60] . We believe that the main reason for the excellent ORR/OER catalytic activity might be attributed to following reasons, First, the doping with the highly electronegative heteroatoms N and S can significantly increase the conductivity of the honeycomb-like carbon porous skeleton, which can further improve the poor conductivity of Co 9 S 8 nanocatalysts through the strong interface coupling between the in situ embedded monodisperse Co 9 S 8 nanoparticles and the outer carbon matrix, providing an advantageous channel for electron mobility. Second, the tendency of the monodisperse Co 9 S 8 nanoparticles to aggregate during continuous work can be greatly relieved through the in situ embedding within the porous carbon skeleton. Moreover, the 3D interconnected porous network with high specific surface area can guarantee a large nanocatalyst loading amount, thus providing more active sites generated from the in situ embedded Co 9 S 8 nanoparticles accessible to the electrolyte. Last but not the least, the 3D interconnected porous network is also beneficial for the mass transport of electrolyte ions and ORR/OER-related species during the catalytic process, allowing a rapid and sustained reaction.

Performance of rechargeable Zn-air battery
Finally, based on the excellent bifunctional oxygen catalytic activity (ORR/OER) of the as-prepared Co 9 S 8 / NSC-900 electrocatalyst, a two-electrode rechargeable zinc-air battery was fabricated. The custom-made zinc-air battery was assembled with a Zn plate serving as an anode, Co 9 S 8 /NSC-900 serving as the air-cathode catalyst, and a 6.0 M KOH solution containing 0.2 M Zn (Ac) 2 serving as an electrolyte. It is known that a minor overpotential region for the charge and discharge Fig. 7 Performance characterization of rechargeable Zn-air battery, a Charge and discharge polarization curves of the fabricated Zn-air battery; b Photograph of a electroluminescence panel powered by two Zn-air batteries with the Co 9 S 8 /NSC-900 based air electrode connected in series; c Galvanostatic discharge-charge cycling curves at 10 mA cm-2 of rechargeable Zn-air batteries with the Co 9 S 8 /NSC-900 and commercial Pt/C-IrO 2 catalyst based air electrode, respectively potential is beneficial for the application of an electrocatalyst in rechargeable zinc-air batteries with highperformance. Figure 7a shows that the discharge/charge voltages of the Co 9 S 8 /NSC-900 bifunctional oxygen electrocatalyst at the current density of 50 mA cm −2 are 1.10/2.25 V, superior to the discharge/charge voltage of commercial Pt/C (1.09/2.47 V) at the same current density, indicating the superior performance of the fabricated rechargeable zinc-air battery based on Co 9 S 8 /NSC-900 eletrocatalyst. Notably, Co 9 S 8 /NSC-900 discharged better at the larger densities. In addition, the electroluminescence panel was successfully powered by two tandem Zn-air batteries based on the Co 9 S 8 /NSC-900 eletrocatalyst (Fig. 7b), directly indicating the great potential of our designed and prepared eletrocatalyst to be used for actual energy storage and conversion devices. The battery cycling test was further conducted to confirm the recharging ability of the fabricated Zn-air battery based on the Co 9 S 8 /NSC-900 eletrocatalyst. When cycled at a current density of 10 mA cm −2 , the battery fabricated on the Co 9 S 8 /NSC-900 eletrocatalyst exhibited considerably stable charge/discharge potentials at 2.14/1.18 V, respectively (Fig. 7c). The whole charging and discharging overpotential for Co 9 S 8 /NSC-900 was ∼0.96 V, approaching that of the commercial Pt/C-IrO 2 composite catalyst (2.0/1.15 V, 0.85 V), corresponding to a high voltaic efficiency of 55.1%, comparable to that of the commercial Pt/C-IrO 2 composite catalyst (57.5%). However, after 4000 min, the polarization potential and current efficiency remained almost unchanged, indicating a better long-term cycling stability than that based on the commercial Pt/C-IrO 2 composite catalyst, for which the cycling stability tended to deteriorate after~1500 min, as shown in Fig. 7c. All of the results indicated the superior competitive power of the honeycomb-structured porous Co 9 S 8 /NSC-900 in situ composites in terms of bifunctional catalytic activities, stabilities, and efficiency for applications in energy storage and conversion devices.

Conclusion
In this work, inspired by honeycombs in nature, a unique honeycomb-structured porous Co 9 S 8 /NSC hybrid bifunctional oxygen electrocatalyst was designed and fabricated with N, S-codoped honeycomb-structured carbon with high conductivity serving as the porous skeleton, within which were homogeneously in situ-embedded monodisperse Co 9 S 8 nanoparticles (2-20 nm) serving as highly active catalysts providing more active sites. Thanks to the rational design of honeycomb-like porous hybrid structure and appropriate material selection, the as-fabricated honeycomb-structured porous Co 9 S 8 /NSC in situ composites not only exhibited higher catalytic activity but also better durability and methanol tolerance compared to that of the commercial Pt/C catalyst. In particular, the half-wave potential (0.88 V vs. RHE) clearly outperformed that of the commercial 20% Pt/C (0.86 V vs. RHE) catalyst in alkaline electrolyte. Furthermore, the Co 9 S 8 /NSC in situ composites also present good catalytic activity and stability for the OER in terms of a small overpotential (0.41 V vs. RHE at 10 mA cm −2 ). Thus, a rechargeable Zn-air battery was further assembled based on the honeycomb-structured porous Co 9 S 8 /NSC in situ composites served as the air electrode with dual electrocatalytic activity, which exhibited a narrow discharge-charge voltage gap (0.96 V) with a high voltaic efficiency of 55.1% at 10 mA cm −2 , comparable to those of an assembled battery based on the commercial Pt/ C-IrO 2 composite catalyst (0.85 V, 57.5%), and a better cycling stability over 4000 min than that of the commercial Pt/C-IrO 2 composite catalyst (~1500 min). Other distinctive features include scalable throughput and a simple, as well as cost-effective fabrication process, which jointly make our designed honeycomb-structured porous Co 9 S 8 / NSC-900 bifunctional eletrocatalyst with superior competitive power appropriate for applications in energy storage and conversion devices.