# Investigation of polymer-derived Si-(B)-C-N ceramic/reduced graphene oxide composite systems as active catalysts towards the hydrogen evolution reaction

Quentin Hanniet, Moustapha Boussmen, Jonathan Barés, Vincent Huon, Igor Iatsunskyi, Emerson Coy, Mikhael Bechelany, Christel Gervais, Damien Voiry, Philippe Miele, et al.

► To cite this version:

Quentin Hanniet, Moustapha Boussmen, Jonathan Barés, Vincent Huon, Igor Iatsunskyi, et al. Investigation of polymer-derived Si-(B)-C-N ceramic/reduced graphene oxide composite systems as active catalysts towards the hydrogen evolution reaction. Scientific Reports, 2020, 10, pp.22003, 1-15. 10.1038/s41598-020-78558-x. hal-03145377

HAL Id: hal-03145377

https://hal.sorbonne-universite.fr/hal-03145377

Submitted on 18 Feb 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

# Investigation of polymer-derived Si-(B)-C-N ceramic/reduced graphene oxide composite systems as active catalysts towards the hydrogen evolution reaction

Quentin Hanniet<fnref n="1" />, Moustapha Boussmen<fnref n="1" />, Jonathan Barés², Vincent Huon², Igor Iatsunskyi³, Emerson Coy³, Mikhael Bechelany<fnref n="1" />, Christel Gervais⁴,⁵, Damien Voiry<fnref n="1" />, Philippe Miele<fnref n="1" />,⁵✉ & Chrystelle Salameh<fnref n="1" />✉

Hydrogen Evolution Reaction (HER) is an attractive technology for chemical conversion of energy. Replacement of platinum with inexpensive and stable electrocatalysts remains a major bottleneck hampering large-scale hydrogen production by using clean and renewable energy sources. Here, we report electrocatalytically active and ultra-stable Polymer-Derived Ceramics towards HER. We successfully prepared ultrathin silicon and carbon (Si-C) based ceramic systems supported on electrically conducting 2D reduced graphene oxide (rGO) nanosheets with promising HER activity by varying the nature and the composition of the ceramic with the inclusion of nitrogen, boron and oxygen. Our results suggest that oxygen-enriched Si-B-C-N/rGO composites (O-SiBCN/rGO) display the strongest catalytic activity leading to an onset potential and a Tafel slope of ~340 mV and ~120 mV dec⁻<fnref n="1" /> respectively. O-SiBCN/rGO electrodes display stability over 170 h with minimal increase of 14% of the overpotential compared to ~1700% for commercial platinum nanoparticles. Our study provides new insights on the performance of ceramics as affordable and robust HER catalysts calling for further exploration of the electrocatalytic activity of such unconventional materials.

Thanks to its high energy density and environment-friendly nature, hydrogen is considered to be one of the most promising alternative energy carrier for the future<fnref n="1" />. Producing hydrogen in an efficient, sustainable and cost-effective way is prerequisite for the hydrogen economy. The most common method for producing hydrogen consists in natural gas reforming but such a process is responsible for the emission of greenhouse gases such as CO₂. Alternatively, considerable efforts have been dedicated towards the development of electrocatalytic water splitting². Platinum (Pt) remains the most active electrocatalyst for the Hydrogen Evolution Reaction (HER: 2H⁺ + 2e⁻ → H₂) and is used in polymer electrolyte membrane electrolyzers. However, Pt is scarce and suffers from high-cost preventing its large-scale utilization³. Platinum is also unstable with time due to poisoning and leaching. Therefore, many electrocatalysts with a small onset potential and a low Tafel slope are currently being explored as alternatives to Pt⁴⁻¹². In particular, low-dimensional materials have been identified as potential catalyst candidates for the electrocatalytic hydrogen evolution. 0D materials such as metal nanoparticles and 1D materials in the form of nanowires and urchin-like structures have been investigated for the HER¹³. As for 2D materials, transition metal dichalcogenides (TMDs) such as MoS₂ have been the most examined¹⁴. Other 2D catalysts include doped graphene¹⁵, graphene-based materials¹⁶⁻¹⁸ and double-layer hydroxides¹⁹. The use

<footnote n="1">¹Institut Européen Des Membranes, IEM, UMR 5635, Univ Montpellier, CNRS, ENSCM, 34095 Montpellier Cedex 5, France. ²Laboratoire de Mécanique Et Génie Civil, UMR 5508, CNRS, Université de Montpellier, 34090 Montpellier, France. ³NanoBioMedical Centre, Adam Mickiewicz University, Wszechnicy Piastowksiej str. 3, 61-614 Poznan, Poland. ⁴Sorbonne Université, UMR 7574 CNRS, Laboratoire de Chimie de la Matière Condensée de Paris, LCMCP, 75005 Paris, France. ⁵Institut Universitaire de France (IUF), MENESR, 1 rue Descartes, 75231 Paris Cedex 05, France. ✉email: philippe.miele@umontpellier.fr; chrystelle.salameh@enscm.fr</footnote>
of HER catalysts is however largely limited by their poor stability in acidic medium due to parasite reactions or leaching, while the recycling of metallic nanocatalysts remains difficult and costly.

Polymer-Derived Ceramics (PDCs) are a class of ceramics synthesized by the thermal decomposition of polymeric precursors that are best known for their excellent high-temperature stability, their remarkable thermo-mechanical properties and their chemical stability; attributed to their unique composition and microstructures<sup>20,21</sup>. Moreover, the final properties of PDCs can be modified in a controlled manner by tailoring the precursor's composition and the processing parameters. The first generation of PDCs was the SiC fibers developed by Yajima in 1975 for high-temperature applications<sup>22</sup>. Since then, many other ternary and quaternary carbides and nitrides have been extensively investigated such as silicon oxycarbide (SiOC)<sup>23</sup>, boron modified silicon carbide (SiBC)<sup>24,25</sup>, silicon carbonitride (SiCN)<sup>26</sup>, boron and aluminum modified silicon carbonitride (SiBCN<sup>27-29</sup>, SiAlCN<sup>21,30</sup>). PDCs display semi-conducting behavior allowing them to be used in Micro-Electro-Mechanical Systems (MEMS)<sup>31</sup>, p-n junction devices<sup>32</sup>, micro glow plugs<sup>33</sup> and pressure sensors<sup>34</sup>. The electronic structure of PDCs has however received little attention because of the amorphous nature of the materials complicating the theoretical calculations of their structure. The fact that PDCs have proved themselves eligible for electrochemical applications e.g. electrodes for lithium-ions batteries has motivated the exploration of PDC-based materials as potential electrocatalysts<sup>35-39</sup>.

Besides, several experimental and theoretical works have demonstrated that doping graphene or other nano-carbons with non-metallic heteroatoms e.g. oxygen, boron and nitrogen strongly enhances the HER<sup>17,40</sup>. Heteroatom doping can indeed tune the hydrogen adsorption strength on graphene in order to comply with the Sabatier principle which suggests that the free energy of hydrogen adsorption (ΔG<sub>H</sub>) should be thermo-neutral<sup>40</sup>. The rational design of inorganic polymer precursors allows the preparation of PDCs with the inclusion of oxygen, nitrogen and boron atoms and therefore holds promises for the optimization of the HER activity of the ceramics.

Here we report the preparation of PDCs and graphene based composites as porous and ultra-stable electro-catalysts towards the production of hydrogen. The PDC/graphene composites were prepared by impregnating 3D reduced graphene oxide (rGO) macroporous foams with different ratios of SiCN and SiBCN preceramic polymer solutions in order to vary their porosity as well as their mechanical properties. The electrochemical activity towards HER was systematically investigated in acidic medium varying the PDCs content and composition. Our results show that PDCs are clearly HER active when supported on conducting graphene and the best HER performance is obtained in the case of O-SiBCN ceramics with an overpotential of ~580 mV vs. reversible hydrogen electrode (RHE) for a current density of 10 mA cm⁻² and a Tafel slope of ~120 mV dec⁻¹. Importantly, the PDC/graphene composites display stable HER activity over almost 170 h with a minimal increase of the potential at 10 mA cm⁻² of only 14%. When testing under the same conditions, supported commercial platinum nanoparticles exhibit a tenfold increase of the overpotential highlighting the clear instability of the Pt catalyst. By taking advantage of the mechanical properties of the PDC/graphene composite, we demonstrate a free-standing electrode for the production of hydrogen that can sustain a current of 100 mA at an overpotential of 660 mV versus RHE. Together with the high stability and durability assessment, the electrochemical investigation was carefully conducted in order to correlate the catalytic activity to the nature, the porosity, the structure and the electrical conductivity of the composites. Our work sheds light on the electrocatalytic activity of inorganic ceramics as affordable and ultra-stable electocatalysts.

## Results and discussion

### Characterization of the PDC/rGO composites.
The chemical structure and the composition were first investigated for two preceramic polymers: PVZ (polyvinylsilazane) and PBVZ (boron-modified polyvinylsilazane) used as precursors of SiCN and SiBCN respectively. PVZ, a yellowish transparent liquid, is a commercial precursor of SiCN. We prepared the PBVZ polymer by treating PVZ with BDMS (borane dimethylsulfide complex) as a hydroboration agent<sup>41</sup>. We chose a controlled Si:B molar ratio equal to 3 to induce a complete hydroboration of the vinyl groups in the PVZ. The reaction between PVZ and BDMS also causes dehydrocoupling reactions between B-H from the BDMS and the N-H units present in the PVZ according to the general chemical formula of PVZ. We first used FT-IR spectroscopy to confirm that the reaction between PVZ and BDMS results in the formation of B-N and B-C bonds (Fig. 1a). The PVZ polymer displays the conventional absorption bands of the different chemical functions with notably the N-H bonds at 3382 cm⁻¹ coupled to the vibration of Si-N at 1170 cm⁻¹<sup>42,43</sup>. In the case of PBVZ, the intensities of the signals attributed to the vinyl groups decrease relative to those of PVZ suggesting thus that the hydroboration reaction occurred. In addition, the increase of the signal around 1325 cm⁻¹ is associated to the emergence of B-N bonds via dehydrocoupling reactions between N-H of PVZ and B-H of the borane complex. The intensity of the peak of N-H bonds at 3382 cm⁻¹ also decreased significantly after the incorporation of boron.

To get additional evidence of the formation of PBVZ, we performed solid-state NMR for identifying the local environment around boron, carbon and silicon atoms. The solid state NMR spectra of <sup>11</sup>B, <sup>13</sup>C and <sup>29</sup>Si recorded for PBVZ are shown in Fig. 1b-d respectively. The <sup>11</sup>B MAS NMR spectra are rather complex due to the heterogeneity in the local chemical environment of the boron nuclei and the presence of quadrupolar interactions. The presence of tri-coordinated boron atoms is however clearly observed. The signals can be fitted with components corresponding to $BC_2N$, $BCN_2$ and $BN_3$ (or also $BN_2H$ since B-H unit is observed by FTIR) environments appearing at $\delta_{iso} \sim 55$ ppm, 40 ppm and 30 ppm respectively in agreement with the literature<sup>41</sup>. These results thus confirm a dehydrocoupling reaction between B-H and N-H resulting in B-N. We note the presence of negligible amount of oxygen detected by the appearance of $BO_3$ and $BO_4$ environments at $\delta_{iso} \sim 15$ ppm, and 0 ppm respectively. The <sup>13</sup>C CP MAS NMR spectrum shown in Fig. 1c displays a signal centered around 10 ppm in the aliphatic region that can be attributed to the superposition of $SiCH_3$ (~2 ppm), $SiCH_2CH_2-B$ (~6 ppm), $SiCHCH_3-B$ (~11 ppm) and $SiCH_2CH_2-B$ (~16 ppm). The presence of these signals strongly supports
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 1.** (a) FTIR spectra of PVZ and PBVZ preceramic polymers. (b) <sup>11</sup>B MAS NMR spectra of PBVZ, O-PBVZ polymers and O-PBVZ/rGO composite. (c) <sup>13</sup>C CP MAS NMR spectrum of PBVZ polymer and (d) <sup>29</sup>Si MAS NMR spectrum of PBVZ polymer.</figcaption></figure>

the hydroboration reaction in perfect agreement with the FT-IR results. The <sup>29</sup>Si MAS spectrum of the polymer in Fig. 1d displays a relatively complex distribution of SiNₓC₄⁻ₓ environments with the obvious presence of H-SiN₂-CH₃ units centered around –20 ppm<sup>44</sup>. The two other peaks observed around –5 ppm and –30 ppm are attributed to SiC₂N₂ and SiCN₃ respectively.

The beneficial role of oxygen as heteroatom in electrocatalysis was recently identified<sup>45,46</sup>. We prepared oxygen-rich PBVZ by exposing the preceramic polymer in air for 48 h to allow the diffusion of oxygen within the polymer and its reaction with the excess of B–H units. Similarly oxygen-rich preceramic polymers were prepared in order to produce O–SiC and O–SiCN. The modification of the C, Si and B environments was explored using NMR. The partial oxidation of PBVZ is confirmed by the increase of the BO₃ and BO₄ signals in the O-PBVZ compared to those in PBVZ (Fig. 1b). The oxidation of the SiNₓ bonds into SiNₓCᵧOz is also identified at –80 ppm as depicted in the <sup>29</sup>Si MAS NMR spectrum (Fig. S1a). The <sup>13</sup>C CP MAS NMR spectrum of O-PBVZ (Fig. S1b) is very similar to that of PBVZ. We also performed <sup>11</sup>B MAS NMR on the O-SiBCN/rGO composite (Fig. 1b). Noticeably, the poor signal-to-noise ratio is explained by the conductive nature of rGO or the possible presence of free radicals preventing the correct tuning of the probe. Importantly, careful examination of the <sup>11</sup>B spectrum reveals that signatures of the boron environment in the composite match those from the PBVZ preceramic polymer. O-SiBCN/rGO composite exhibits the same environments with nonetheless a higher proportion of BO₃
<figure><img src="image_9.png" /><figcaption>**Figure 2.** Low-magnification SEM image of O-SiBCN/rGO (top left panel) and the corresponding elemental maps.</figcaption></figure>

and $BO_4$ environments compared to $BC_2N$, $BCN_2$ and $BN_3$ confirming that mixing rGO and PDCs does not alter the chemical structure of the ceramic in the composite. Importantly, the absence of changes from the chemical environment of boron points out that the formation of the ceramics upon annealing does not induce any doping of the graphene support. Thermogravimetric analyses were performed on the different PDC/rGO composites in order to study their intrinsic thermal behavior during the polymer-to-ceramic conversion. The TGA curves display the weight loss of the PVZ, PBVZ and O-PBVZ when increasing the temperature up to 1000 °C under nitrogen (Fig. S2). Our results show that the addition of boron significantly decreases the weight losses of the PVZ improving the ceramic yield from 60 to 78% for PVZ and PBVZ respectively. We note that a high ceramic yield is strongly desired for the retention of the porosity, shape and structural properties of the impregnated rGO foam. Remarkably, the addition of oxygen does not drastically affect the weight loss of the PBVZ and the ceramic yield is maintained as high as 77%.

Next, we characterized the structure of the rGO and O-SiBCN/rGO composites using SEM and TEM. Figure 2 shows the low-magnification SEM images and the EDX elemental mapping of the O-SiBCN/rGO composite with a 1:1 ratio. The composite exhibits clear porosity, while the presence of B, N and Si suggests the successful infiltration of the ceramic inside the porosity of the rGO foam. High-magnification SEM images of rGO and O-SiBCN/rGO (Fig. 3a) confirm that the impregnation of rGO with the polymer occurred without altering the porous structure of the foam. This supports the fact that the pores are uniformly interconnected in the 3D composite and the foams retain their high porosity. Because of the macroscopic nature of our composite, we performed mercury intrusion porosimetry and we found a porosity of 93% for O-SiBCN/rGO vs 97% for rGO with an average pore diameter around 5 µm. The increase of the ceramic ratio in the composite leads to the progressive filling of the pores of the rGO scaffold (Figs. 3a and S3a). In order to further confirm the formation of ceramics on the rGO nanosheets, TEM and EELS observations were carried out on the composite (Fig. 3b,c). The two-dimensional nature is identified from the TEM images and attributed to rGO nanosheets. High-resolution TEM (HRTEM) images reveal that the wall of the foam consists of 4-8 rGO layers with sizes ranging from a 50 to few hundred nanometers (Fig. 3b). The interplanar distance ($d_{002}$) for rGO is estimated to ~3.85 Å—close to interlayer distance of stacked graphene nanosheets of 3.43 Å. The presence of O-SiBCN on the top of rGO can be observed in Fig. 3c, while the EELS elemental mapping clearly confirms that all the components of the ceramic i.e. Si, B, C and N are present and homogenously distributed in the sample. We note that the diffraction pattern suggests a rather disordered structure of the composites as expected from the turbostatic stacking of the graphene nanosheets and the amorphous nature of the ceramic (Inset Fig. 3b). O-SiCN/rGO and O-SiC/rGO composites with a 1:1 ratio were also observed (Fig. S3b,c). Similarly to O-SiBCN/rGO, the porosity of the composite decreases with
<figure><img src="image_10.png" /><figcaption>O-SiBCN/rGO R=1</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>O-SiBCN/rGO R=2</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /></figure>

<figure><img src="image_18.png" /></figure>

<figure><img src="image_19.png" /></figure>

<figure><img src="image_20.png" /></figure>

<figure><img src="image_21.png" /><figcaption>**Figure 3.** (a) High-magnification SEM images of rGO foam, O-SiBCN/rGO and O-SiBCN/rGO composite with a PDC/rGO ratio of 1 and 2 respectively. (b) TEM images of O-SiBCN/rGO (R<sub>O-SiBCN/rGO</sub> = 1). (c) TEM images and the corresponding EELS chemical maps of O-SiBCN/rGO (R<sub>O-SiBCN/rGO</sub> = 1) (Silicon (blue), Boron (green), Carbon (yellow), Nitrogen (red)).</figcaption></figure>
<figure><img src="image_22.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_23.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_24.png" /></figure>

**Figure 4.** (a) XPS spectra of the N1s, C1s, B1s, Si2p regions of O-SiBCN/rGO (top) and O-SiCN/rGO composites (down). (b) TGA of the rGO foam, O-SiBCN ceramic, O-SiBCN/rGO composite under air up to 1000 °C and (c) electrical conductivity of PDCs/rGO composite systems.

the increase of the PDC ratio (Fig. S3a). The presence of silicon carbide (from Si edge) and graphitic carbon (C edge) with small contents of amorphous carbon are detected from the EELS spectra.

X-ray photoelectron spectroscopy was then conducted to investigate the chemical composition of O-SiBCN and O-SiCN based composites pyrolyzed at 1000 °C (Figs. 4a and S4). The survey spectra of the composites confirm the presence of all the elements from the ceramics: Si, B, C, N and O, while no metals impurities can be detected (Fig. S4). High-resolution spectra of O-SiBCN/rGO from the core of the composites are presented in Fig. 4a. The intense peak at 284.4 eV is attributed to the C=C bonds from sp² carbons in rGO. The peak at 285.2 eV suggests the presence of amorphous carbon while the existence of remaining C–O and C=O bonds comes from the incomplete reduction of the GO foams during the reduction step⁴⁷. The Si–O, Si–N, N–SiOₓ and N–Si peaks at 103.4 eV, 102.2 eV, 398.4 eV and 397.8 eV respectively confirm the existence of the SiNₓOᵧ
phase in the microstructure of the composite. The B1s and N1s regions can be deconvoluted with peaks at 190.7 eV, 192.3 eV and 193.4 eV from B–C, B–N and B–O respectively and N–B at 399.1 eV in agreement with the formation of the BNₓOᵧ phase. Remarkably, the value of the Si:B ratio is estimated to 3—fully consistent with the expected Si:B ratio for O-SiBCN. This hints that no reaction occurred between the polymer and the rGO during the impregnation and the polymer-to-ceramic conversion process, strongly supporting the observations of NMR. The formation of SiCN ceramic in O-SiCN/rGO composites was also confirmed from the XPS analyses (Fig. 4a). Si–C, Si–N, and Si–O bonds are detected at 101.4, 102.4 and 103.7 eV respectively from the Si2p region, while the peaks of N–SiOₓ, N–Si, Si–N and Si–C are identified at 399.6, 398, 102.3 and 99.3 eV. XPS spectra of deconvoluted O1s are shown in Fig. S5.

To further confirm the presence of the ceramic in the composite, we performed TGA under air atmosphere (Fig. 4b). The weight loss at 1000 °C reaches 63% for O-SiBCN/rGO composites compared to 100% and 0.5% for rGO and O-SiBCN respectively. The weight loss for the composite is attributed to the combustion of rGO leaving the O-SiBCN ceramic as confirmed by our EDX analyses (Fig. S6). Overall, TGA and XPS confirm the presence of the ceramics within the different PDC/rGO composites. The electrical conductivities of the catalytic systems were examined for different ceramic ratios (Fig. 4c). We found that the conductivity rapidly decreases when increasing the ceramic ratio in the composite. For a ratio PDC:rGO equal to 1, the electrical conductivity is estimated to 713 S m⁻¹ close to the conductivity of pristine rGO foam at ~1000 S m⁻¹. This, combined with the high porosity of the PDC/rGO composites, is important for the fabrication of electrodes with high specific surface area.

**Electrocatalytic measurements.** The electrocatalytic performance of the pristine ceramics and the rGO-based composites was investigated towards the evolution of hydrogen in 0.5 mol L⁻¹ H₂SO₄. The HER activity from the glassy carbon was also measured to confirm that no activity arises from the electrode support. Figure 5 shows the electrochemical behavior of the different composite electrodes for a PDC:rGO ratio of 1—larger ratios of ceramics were also tested but do not lead to higher HER activity (Fig. S7). First, we estimated the electrochemical active surface area (ECSA) by measuring the double layer capacitance (C₉) of the composites (Fig. 5). The ECSA was calculated assuming an ideal capacitance of 0.015 mF cm⁻² for a flat surface of graphene—note that the contribution of the ceramic can be neglected due to its insulating nature (Fig. S8)⁴⁸. All the composites exhibit larger C₉ compared to the glassy carbon support revealing the porosity of the electrode in agreement with our SEM observations (Fig. 5b). The porous network combined with the high conductivity of the rGO scaffold provides (i) accessibility to more active sites and (ii) efficient transport pathways of the electrons to the active sites. Interestingly the composite electrodes display larger C₉ than pristine rGO (420 vs. 316 µF cm⁻²). This is attributed to the mechanical reinforcement of rGO by the ceramics that prevents the porous structure from collapsing. Of all the different PDC/rGO composites, O-SiBCN/rGO exhibits the largest C₉ and therefore the largest ECSA of 28 cm⁻² ECSA.

The polarization curves from the different composite catalysts i.e. SiC/rGO, O-SiC/rGO, SiCN/rGO, O-SiCN/rGO, SiBCN/rGO and O-SiBCN/rGO are presented in Fig. 5c. The HER activity is clearly increased in the case of the PDC/rGO compared to glassy carbon and rGO alone—ruling out any possible activity from impurities present in rGO or the conducting support. To examine the role of rGO on the catalytic performance, we assess the electrochemical behavior of the ceramics alone (Fig. S8a,b). We can see that the ceramics display poor HER activity demonstrating that rGO allows activating the PDCs via a synergistic effect as previously described for supported semiconducting materials⁴⁹.

Remarkably O-SiBCN/rGO demonstrated the best performance among the different PDC/rGO electrodes with an onsetpotential of -340 mV vs. RHE and a current density as high as 60 mA cm⁻² at an overpotential of η=660 mV. O-SiBCN/rGO also exhibits the lowest Tafel slope at 127 mV dec⁻¹ suggesting faster kinetics and the overpotential at 10 mA cm⁻² reaches 590 mV (Fig. 5c). For comparison, overpotentials of η=51 mV and η=20 mV for driving a current of 10 mA cm⁻² have been measured using respectively Pt-based alloyed catalysts such as Pt–Pd supported on graphene oxide and Pt–Ni–Co⁵⁰,⁵¹. As for noble metal-free catalysts such as phosphides, overpotentials of 130–200 mV are requested to produce a current of 10 mA cm⁻²⁵²; comparable to that of MoS₂ deposited on graphene⁴⁹. The decrease in Tafel values for the composites compared to their ceramic counterparts emphasizes the role of rGO as a conductive support. It has been shown that a Tafel slope value of 120 mV dec⁻¹ corresponds to a kinetically limiting stage attributed to the hydrogen adsorption: Volmer step, H⁺ + e⁻ → H* ⁵³. In the case of O-SiBCN/rGO, the value of 127 mV dec⁻¹ suggests that the reaction on the catalyst surface proceeds via the Volmer-Heyrovsky mechanism with the Volmer step being the rate determining step (RDS). To obtain a fair comparison of the electrodes and avoid any contribution from the porosity of the different composites, we normalized the geometrical current density to the ECSA (Fig. 5e). We estimated the ECSA-normalized current density at an overpotential of η=650 mV to compare the different composites and the performances follow the trend: rGO < O-SiC/rGO ≈ SiC/rGO < SiBCN/rGO ≈ SiCN/rGO < O-SiCN/rGO < O-SiBCN/rGO (Fig. 5f).

Based on our electrochemical measurements and the physical characterization of the composite, we speculate that the active sites in the composite O-SiBCN/rGO are likely to be located at BN₃, BCₓNᵧ sites of the ceramic. This calls for additional investigation notably using first principle calculations that are however beyond the scope of this paper.

To get further understanding of the nature of the active sites of the O-SiBCN ceramics, we investigated the effect of the crystallinity of the ceramics on the electrocatalytic activity of O-SiBCN/rGO and O-SiCN/rGO by annealing the composites at increasing temperatures (Fig. 6a,b). Other amorphous catalysts such as metal chalcogenides, metal alloys, metal phosphides have already demonstrated larger catalytic activity relative to their crystalline counterparts⁵⁴⁻⁵⁶. Both O-SiBCN and O-SiCN catalysts remain active after annealing at temperatures
<figure><img src="image_25.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_26.png" /><figcaption>e</figcaption></figure>

<figure><img src="image_27.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_28.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_29.png" /><figcaption>f</figcaption></figure>

<figure><img src="image_30.png" /></figure>

**Figure 5.** (a) Cyclic voltammetry curves for O-SiBCN/rGO composite. (b) Estimation of the double layer capacitance (C<sub>d</sub>) for the different electrodes. (c,d) Polarization curves and Tafel slopes of rGO, SiC/rGO, O-SiC/rGO, SiCN/rGO, O-SiCN/rGO, SiBCN/rGO and O-SiBCN/rGO catalysts. (e) Corresponding ECSA-normalized polarization curves and (f) Histograms displaying the ECSA-normalized current density (J) at η = 650 mV for the different catalysts.

up to 1200 °C but then lose most of the activity at 1400 °C and above. To elucidate the origin of the activity loss,
the PDCs were characterized using X-ray diffraction and the PDC/rGO composites using Raman spectroscopy.
Both O-SiBCN and O-SiCN ceramics progressively crystallize at temperature higher than 1200 °C as shown by
the rise of the peaks at 35.7° (111), 60.3° (220) and 71.25° (311) corresponding to the formation of nanosized
β-SiC crystals (Fig. 6c,d). These peaks are more obvious in the case of O-SiCN (Fig. 6d). This is attributed to the
presence of boron within the ceramic that delays the crystallization of Si–B–C–N systems compared to Si–C–N
ones. The incapacity of Si₃N₄ to react with free carbon via the carboreduction reaction: Si₃N₄ + 3C → SiC + 2N₂
primarily induces the crystallization into β-Si₃N₄ as confirmed by the peaks detected at 26.5, 33.7° and 36.1°.
Of note, the presence of oxygen has no deleterious effect on the amorphous behavior, which may suggest that
<figure><img src="image_31.png" /><figcaption>**Figure 6.** Polarization curves of O-SiBCN/rGO (a) and O-SiCN/rGO (b) composites at different temperatures. (c,d) X-Ray diffractograms of O-SiBCN ceramics and O-SiCN ceramics at different temperatures. (e,f) Corresponding Raman spectra of O-SiBCN/rGO and O-SiCN composites.</figcaption></figure>

the selected ceramics SiBCN and O-SiBCN share similar microstructural framework. The increase of the composites crystallinity is also confirmed by the strengthening of Raman signatures from rGO graphitic domains (Fig. 6e,f). The position of the G band of O-SiBCN/rGO pyrolyzed between 1000 °C and 1600 °C continuously blue-shifts from 1587.5 to 1579.7 cm⁻¹ due to the transition from amorphous disordered carbon to graphitic domains in agreement with the increase of the 2D band⁵⁸. Overall our results point out the importance of the amorphous ceramics—notably O-SiBCN—on the electrocatalytic activity towards the production of hydrogen.

To further explore the potential of the O-SiBCN/rGO composites, we studied the stability of the HER activity in acidic medium. The electrode was continuously subjected to various current densities for a total of 168 h. During the stability test, the applied current density was set to 2.5 mA cm⁻², 5 mA cm⁻² and 7.5 mA cm⁻² for 48,
<figure><img src="image_32.png" /><figcaption>**Figure 7.** (a) Chronopotentiometry measurements for supported Pt/C electrodes and O-SiBCN/rGO composite over 168 h. (b) Change in overpotential (%) for O-SiBCN/rGO and other electrocatalysts in the literature: (Ni-FeP/C-hnr⁶² (hnr=hallow nanorod), Ni₂P/CNFs⁶³ (CNFs=carbon nanofibers), NiP₂/CNFs⁶³, Ni foam⁶⁴, NiCo-OH/NF⁶⁴ (NF=nickel foam), Ni₂P/NF⁶⁴, NiCoP/NF⁶⁴, Ni-C-N NSs⁶⁵ (NSs=nanosheets), Ni-Mo⁶⁶, Fe-Ni-P⁶⁷, CM/Co/SiO₂⁶⁸ (CM=carbon material), SAS/Co/HOPNC⁶⁹ (SAS=single-atomic Co sites, HOPNC=embedded in hierarchically ordered porous N-doped carbon), Ir/C₃N₄/NG⁷⁰, Pt/Ru ALD⁷¹ (ALD=atomic layer deposition), Pt/SA⁷¹ (SA=single atomic), Pt/C⁷¹, Ru-MoS₂⁷², MoS₂ 1T⁷³, Et₂NPh fct MoS₂⁷³, TiO₂/MoS₂⁷⁴, O-WS2-1T⁷⁵, O-WS2-2H⁷⁵, WSe₂ NFlw⁷⁶ (NFlw=nano flowers). (c) Compression curves of the rGO foam and the O-SiBCN/rGO composite. The inset shows that the composite remains stable after compression. (d) Comparison of the polarization curves from the O-SiBCN/rGO composite catalyst deposited on the glassy carbon electrode and free-standing O-SiBCN/rGO. Note that the potentials of the free-standing electrode were not iR loss corrected.</figcaption></figure>

24 and 96 h, respectively. The modulation of the requested current density aimed at simulating the variation of the supply power similarly to the case of an electrolyzer connected to a photovoltaic cell. For a fair comparison commercial 5%Pt/C electrode was subjected to the same procedure (Fig. 7a). As expected commercial Pt nano-particles initially show better HER performance as shown by the low overpotential of 51 mV at 2.5 mA cm⁻²—smaller than in the case of O-SiBCN/rGO of 454 mV. Remarkably, the overpotential for Pt rapidly increases by 493% during the first 12 h suggesting a rapid loss of activity. This is likely due to surface poisoning of the Pt surface by impurities from the electrolyte solution as well as possible leaching of the Pt nanoparticles in the acid electrolyte. Conversely, O-SiBCN/rGO exhibits significantly better retention of the HER performance over ≈170 h with an increase of the overpotential of 18.4%, 0.3% and ≈0% after 12, 24 and 48 h at 2.5 mA cm⁻², 5 mA cm⁻² and 7.5 mA cm⁻², respectively. After 168 h, the difference of the overpotential between O-SiBCN/rGO and 5%Pt/C sharply decreases from 424 mV to only 85 mV. These results reveal the interest of O-SiBCN/rGO composite as an ultra-stable electrocatalyst for the HER in acidic medium. We also compared the stability of O-SiBCN/rGO electrode with that of several electrocatalysts reported previously in the literature. Our results clearly show outstanding stability of O-SiBCN/rGO over a long period of time demonstrating thus the potential of such materials towards realistic operation (Fig. 7b).
To demonstrate the advantage of the PDC/rGO composites for electrocatalysis, the mechanical properties of O-SiBCN/rGO and pristine rGO foams were evaluated by compression tests. The Young's moduli for both pristine and PDC/rGO composites were deduced from the first linear part of stress-strain curves (Fig. 7c). Our results show that the O-SiBCN/rGO and rGO samples have a plastic behavior with a plasticity cross-over around 13% and 11% strain respectively. The Young's modulus of O-SiBCN/rGO is estimated to almost 50 kPa—approximately fourfold larger than that of rGO. This considerable improvement of the stiffness of the composites is attributed to the presence of the ceramic on the rGO nanosheets. The plasticity modulus, extrapolated from the second linear part of the strain-stress curves, also decreases from 50 kPa for rGO to 25 kPa for O-SiBCN. Our mechanical investigations thus reveal that the controlled coating of the rGO foam with the ceramics strongly benefits the mechanical properties of the composite and allows the preparation of mechanically-robust electrodes. To validate our concept, we fabricated a free-standing electrode based on O-SiBCN/rGO. Thanks to the mechanical properties of the ceramics reinforcement and the electrical conductivity of the rGO scaffold, the electrode can directly be used without any conducting support. Such rational design of the composite electrode demonstrates a high porosity thanks to the rGO nanosheets, while the HER activity originating from the O-SiBCN ceramic allows reaching a current of almost 130 mA (Fig. 7d).

## Conclusion

Ceramics have been largely disregarded in electrocatalysis because of their insulating behavior. Here, we demonstrated an example of engineered PDC/rGO based composites as electrocatalytic active materials. Within the composite, rGO acts as a porous conductive network, with a high porosity of 93% and a high electrical conductivity allowing fast transport of the charge carriers and consequently activating the supported ceramic. In-depth investigations revealed that the HER activity of the ceramics can be modulated by the presence of O, N and B in the Si–C systems. Our careful measurements based on the ECSA and linear voltammetry allowed identifying O-SiBCN ceramic as stable HER catalyst. O-SiBCN/rGO exhibits the best performance with an onset potential of –340 mV vs. RHE, while exhibiting an outstanding stability over almost 170 h with a minimal increase of the overpotential. Importantly, our results demonstrate the potential of ceramics for electrocatalysis when prepared in the form of an ultrathin film deposited on a conducting substrate. Considering the great stability of ceramics, we believe our work opens perspectives in the rational design of ceramic/graphene composites with controlled composition and rational structure and sheds light on novel applications of inorganic ceramics in electrocatalysis.

## Methods

**Reagents and materials.** Polyvinylsilazane (PVZ, also called durazane)—a commercial preceramic polymer for SiCN—was purchased from Merck. AllylHydridoPolyCarboSilane (SMP-10, abbreviated AHPCS)—a commercial preceramic precursor for SiC—was purchased from Starfire. Toluene (anhydrous 99.8%, Sigma Aldrich), borane dimethylsulfide complex (BDMS, 2 M in toluene under argon, Alfa Aesar), Nafion (D-520 dispersion, 5%w/w in water and 1-propanol, Sigma Aldrich), isopropanol (100%, VWR chemical), potassium permanganate (KMNO₄, 99.0%min., Alfa Aesar), hydrogen peroxide (H₂O₂, ACS, 29–32% aqueous solution, Alfa Aesar) and sulfuric acid (H₂SO₄, ACS reagent 95–98%, Sigma Aldrich) were used as received without further purification.

**Materials synthesis.** Synthesis of the preceramic polymer precursor of SiBCN. All chemicals and synthesis steps were handled in an argon filled glovebox (MBraun, MB-Labstar, O₂ and H₂O < 0.5 ppm) and using Schlenk techniques on a vacuum/argon line (with argon purified through a Siccapent phosphorous pentoxide column and BTS catalyst). PVZ and BDMS reacted through hydroboration and dehydrocoupling reactions in a quantitative yield leading to a boron-modified polyvinylsilazane abbreviated as PBVZ. BDMS was added dropwise at 0 °C into a three-neck round-bottom flask containing PVZ and toluene under stirring. When the addition was complete, the mixture was left under stirring at RT for 3 days. The solvent was then removed by distillation at reduced pressure. A moisture and air-sensitive glass-like white polymer was obtained with a yield of 81 wt%. The same procedure was applied to synthesize the precursor of SiBC by reacting AHPCS with the borane complex yielding also a B-AHPCS glass-like white polymer.

*Synthesis of rGO foams.* In a first step, GO was synthesized according to a modified Hummers method starting from graphite powder⁵⁹. In a typical synthesis, the graphite was oxidized by concentrated sulfuric acid and potassium permanganate. Deionized water was slowly added, followed by hydrogen peroxide. GO was then washed with deionized water and subsequently centrifuged several times until a concentrated solution with a neutral pH was obtained. GO foams were synthesized by hydrothermal synthesis; where aqueous solutions of GO (1.5 mg mL⁻¹) were introduced in an autoclave (Parr instrument company) and heated at 180 °C for 12 h. The as-formed wet gels were then fully frozen and consequently freeze-dried (Labonco Freezone 4.5; 0.05 mBar) yielding 3D GO foams. The reduction of the GO into reduced graphene oxide (rGO) foams was carried out by thermal treatment at 900 °C for 2 h under reducing atmosphere (Varigon, 5%H₂/Ar) in a Vecstar tubular furnace.

*Preparation of PDC/rGO composites.* PDC/rGO composites were prepared via a vacuum-assisted impregnation process in a Schlenk-type flask. The preparation of SiCN/rGO is described hereafter as an example. The preceramic precursor i.e. PVZ was dissolved in toluene under stirring and added to a rGO foam. rGO was previously dehydrated at 120 °C under vacuum for 12 h to remove all the trapped air and moisture between the pores of the foam. The mixture was allowed to infiltrate the pores of the foam by capillarity under static vacuum for 24 h. Different ratios between the preceramic precursors and the foams were prepared varying from 1 to 10 in weight. After washing with toluene (4×10 mL) and drying, the PVZ/rGO composite was then transferred into a quartz
tube in a tubular furnace (Nabertherm type RS 80/500/11 Germany) and subjected to a pyrolysis under nitrogen at 1000 °C for 2 h with a heating rate of 1 °C min⁻¹ and a cooling rate of 2 °C min⁻¹. The slow heating rate generates nanoporosity within the ceramic enhancing thus the surface area as previously described⁶⁰ while avoiding the collapse of the composite structure. This thermal treatment leads to the ceramic conversion (PVZ → SiCN) and further reduces the rGO in order to have a chemical structure as close as possible to graphene. Other PDC-based composites were prepared accordingly with the following compositions: SiBCN/rGO, SiC/rGO and SiBC/rGO by using PBVZ, AHPCS and B-AHPCS as preceramic polymers respectively.

**Crystallization of PDC/rGO composites.** The PDC/rGO composites were subsequently annealed at high temperatures (1200–1600 °C) for 2 h under nitrogen in a graphitic furnace (Gero, Model HTK8) in order to induce crystallization of the ceramics. During all the thermal treatments, the materials were loaded in the furnace, after what the furnace chamber was evacuated for 1 h then filled with nitrogen with a 200 mL.min⁻¹ flow rate during the treatment.

**Materials characterization.** The chemical structure of the preceramic polymers was characterized by Fourier Transform InfraRed (FT-IR) spectroscopy with a Nicolet iS50 Thermo scientific spectrometer using KBr pellets (2 wt% of polymer mixed with previously dried KBr powder followed by compaction into a dense pellet). It was also characterized by solid-state NMR spectroscopy. Solid-state ¹³C CP MAS and ²⁹Si MAS NMR were recorded on a Bruker Avance 300 spectrometer (B₀ = 7.0 T, v₀ (¹H) = 300.29 MHz, v₀ (¹³C) = 75.51 MHz, v₀ (²⁹Si) = 59.66 MHz) using a 4 mm Bruker probe and spinning frequency of 10 kHz. While ¹¹B MAS spectra were recorded on a Bruker Avance 700 spectrometer (B₀ = 716.4 T, v₀ (¹¹B) = 224.6 MHz) using a 3.2 mm Bruker probe and spinning frequency of 20 kHz. ¹¹B spectra were acquired using a spin-echo θ-τ-2θ pulse sequence with θ = 90°, to overcome problems related to the probe signal. The τ delay (50 µs) was synchronized with the spinning frequency. The chemical shift values were referenced to TMS for ¹³C and ²⁹Si, and BF₃.OEt₂ for ¹¹B. The NMR spectra were simulated with the DMFIT program⁶¹.

Thermogravimetric analysis (TGA) was carried out on a TGA-STD Q600 thermal analysis device under atmospheric pressure from RT up to 1000 °C with a heating rate of 5 °C.min⁻¹. TGA was performed under nitrogen in the case of the preceramic polymers to investigate their ceramic conversion and under air for the PDC/rGO composites to study the effect of the thermal annealing in a reactive atmosphere. Scanning Electron Microscopy (SEM) (Hitashi S4800, operating with an acceleration voltage between 0.1 kV and 30 kV) and Transmission Electron Microscopy (TEM) measurements (JEOL ARM 200F operating at 200 eV, equipped with an energy dispersive X-Ray spectroscopy—EDX and an energy electron loss spectroscopy—EELS) were used to investigate the morphology and microstructure of the composites. Energy dispersive X-ray spectroscopy (Detector: Oxford Instruments X-Max N SDD; Microscopy: Zeiss EVO HD15) and X-ray photoelectron spectroscopy (XPS (ESCALAB 250, ThermoElectron equipped with an Al Kα line 1486.6 eV as monochromatic excitation source with a 400-µm spot size) were used to analyze the chemical composition of the ceramics and composites. The samples for XPS were prepared on a silicon wafer coated with gold. Raman spectroscopy (LABRAM 1B) and X-ray diffraction (Philips X'pert PRO X-ray diffraction system operating at 20 mA and 40 kV from 10 to 90° with a step size of 0.0167, using a Kα1 of copper as source) were used to study the evolution of the structure of the PDCs with and without rGO with the temperature. The electrical conductivity was measured by using a two-probe method with silver as contact metal. Compression measurements were carried out on the rGO foam and on the O-SiBCN/rGO composite using a LLOYD 01/LF compression device equipped with a 50 N capacity force sensor. The porosity of the foams was characterized by mercury intrusion porosimetry using a Micromeritics Autopore IV 9500 penetrometer.

**Electrochemical measurements.** Electrochemical characterizations were performed in a traditional three-electrode setup using a VSP500 potentiostat from Biologic in a 0.5 mol L⁻¹ H₂SO₄ electrolyte solution. A Hg₂Cl₂/Hg electrode and graphite rod (Alfa Aesar) were used as the reference and the counter electrodes respectively. The working electrodes were prepared by dispersing 5 mg of the catalyst in 1.0 mL of deionized water, isopropanol and 5 wt% Nafion solution (at 62.5; 25 and 12.5 vol%, respectively). 10 µL of the ink was loaded onto a glassy carbon (Sigradur) electrode with 0.2 cm² active area corresponding to a weight density of 1.25 mg cm⁻². All potentials were corrected to the reversible hydrogen electrode (RHE). The electrochemical impedance spectroscopy has been recorded at 550 mV vs. RHE with an AC voltage of 5 mV from 10⁶ to 0.01 Hz. The double-layer charge capacitance (C₉₉) was measured by cycling the electrodes at increasing steps from 5 mV s⁻¹ up to 250 mV s⁻¹ between 241 and 91 mV vs. RHE. The C₉₉ for the different electrodes was estimated from the slope of the anodic and cathodic current densities at 166 mV vs. RHE. All potentials presented were corrected for iR losses unless specifically noted.

Received: 7 October 2020; Accepted: 23 November 2020

Published online: 15 December 2020

**References**

1. Dresselhaus, M. S. & Thomas, I. L. Alternative energy technologies. *Nature* **414**, 332 (2001).

2. Voiry, D. *et al*. The role of electronic coupling between substrate and 2D MoS 2 nanosheets in electrocatalytic production of hydrogen. *Nat. Mater.* **15**, 1003 (2016).

3. Mahmoud, M. H. H. Leaching platinum-group metals in a sulfuric acid/chloride solution. *JOM* **55**, 37–40 (2003).
4. Le Goff, A. et al. From hydrogenases to noble metal-free catalytic nanomaterials for H2 production and uptake. *Science* **326**, 1384–1387 (2009).

5. Subbaraman, R. et al. Enhancing hydrogen evolution activity in water splitting by tailoring Li+-Ni (OH) 2-Pt interfaces. *Science* **334**, 1256–1260 (2011).

6. Kibsgaard, J., Chen, Z., Reinecke, B. N. & Jaramillo, T. F. Engineering the surface structure of MoS 2 to preferentially expose active edge sites for electrocatalysis. *Nat. Mater.* **11**, 963 (2012).

7. Voiry, D. et al. Enhanced catalytic activity in strained chemically exfoliated WS 2 nanosheets for hydrogen evolution. *Nat. Mater.* **12**, 850 (2013).

8. Voiry, D., Yang, J. & Chhowalla, M. Recent strategies for improving the catalytic activity of 2D TMD nanosheets toward the hydrogen evolution reaction. *Adv. Mater.* **28**, 6197–6206 (2016).

9. Yang, J. et al. Two-dimensional hybrid nanosheets of tungsten disulfide and reduced graphene oxide as catalysts for enhanced hydrogen evolution. *Angew. Chemie Int. Ed.* **52**, 13751–13754 (2013).

10. Liu, Q. et al. Electron-doped 1T-MoS2 via interface engineering for enhanced electrocatalytic hydrogen evolution. *Chem. Mater.* **29**, 4738–4744 (2017).

11. Li, L. et al. Role of sulfur vacancies and undercoordinated MO regions in MoS2 nanosheets towards the evolution of hydrogen. *ACS Nano* **13**, 6824–6834 (2019).

12. Coy, E. et al. High electrocatalytic response of a mechanically enhanced NbC nanocomposite electrode toward hydrogen evolution reaction. *ACS Appl. Mater. Interfaces* **9**, 30872–30879 (2017).

13. Yang, H., Zhang, Y., Hu, F. & Wang, Q. Urchin-like CoP nanocrystals as hydrogen evolution reaction and oxygen reduction reaction dual-electrocatalyst with superior stability. *Nano Lett.* **15**, 7616–7620 (2015).

14. Yang, J. & Shin, H. S. Recent advances in layered transition metal dichalcogenides for hydrogen evolution reaction. *J. Mater. Chem. A* **2**, 5979–5985 (2014).

15. Ferrari, A.G.-M., Brownson, D. A. C. & Banks, C. E. Investigating the integrity of Graphene towards the electrochemical Hydrogen evolution Reaction (HeR). *Sci. Rep.* **9**, 1–10 (2019).

16. Huang, H. et al. Graphene nanoarchitectonics: recent advances in graphene-based electrocatalysts for hydrogen evolution reaction. *Adv. Mater.* **31** 1903415 (2019).

17. Ito, Y., Cong, W., Fujita, T., Tang, Z. & Chen, M. High catalytic activity of nitrogen and sulfur co-doped nanoporous graphene in the hydrogen evolution reaction. *Angew. Chemie Int. Ed.* **54**, 2131–2136 (2015).

18. Meng, X. et al. Three-dimensionally hierarchical MoS2/graphene architecture for high-performance hydrogen evolution reaction. *Nano Energy* **61**, 611–616 (2019).

19. Luo, J. et al. Water photolysis at 12.3% efficiency via perovskite photovoltaics and Earth-abundant catalysts. *Science* **345**, 1593–1596 (2014).

20. Colombo, P., Mera, G., Riedel, R. & Soraru, G. D. Polymer-derived ceramics: 40 years of research and innovation in advanced ceramics. *J. Am. Ceram. Soc.* **93**, 1805–1837 (2010).

21. Salameh, C. et al. Chemistry of a series of aluminum-modified polysilazanes: synthesis, pyrolysis behaviour and microstructural evolution. *J. Eur. Ceram. Soc.* **39**, 183–194 (2019).

22. Yajima, S., Hayashi, J. & Omori, M. Continuous silicon carbide fiber of high tensile strength. *Chem. Lett.* **4**, 931–934 (1975).

23. Soraru, G. D., D'andrea, G., Campostrini, R., Babonneau, F. & Mariotto, G. Structural characterization and high-temperature behavior of silicon oxycarbide glasses prepared from sol–gel precursors containing Si–H bonds. *J. Am. Ceram. Soc.* **78**, 379–387 (1995).

24. Kaur, S., Mera, G., Riedel, R. & Ionescu, E. Effect of boron incorporation on the phase composition and high-temperature behavior of polymer-derived silicon carbide. *J. Eur. Ceram. Soc.* **36**, 967–977 (2016).

25. Schmidt, M. et al. Molecular-level processing of Si-(B)-C materials with tailored nano/microstructures. *Chem. A Eur. J.* **23**, 17103–17117 (2017).

26. Mera, G., Navrotsky, A., Sen, S., Kleebe, H.-J. & Riedel, R. Polymer-derived SiCN and SiOC ceramics–structure and energetics at the nanoscale. *J. Mater. Chem. A* **1**, 3826–3836 (2013).

27. Viard, A. et al. Molecular design of melt-spinnable co-polymers as Si–B–C–N fiber precursors. *Dalt. Trans.* **46**, 13510–13523 (2017).

28. Bechelany, M. C. et al. Preparation of polymer-derived Si-B-C-N monoliths by spark plasma sintering technique. *J. Eur. Ceram. Soc.* **35**, 1361–1374 (2015).

29. Viard, A. et al. Molecular chemistry and engineering of Boron-modified polyorganosilazanes as new processable and functional SiBCN precursors. *Chem. A Eur. J.* **23**, 9076–9090 (2017).

30. Singh, G. & David, L. Aluminum-modified polysilazanes for polymer-derived ceramic nanocomposites (2018).

31. Liew, L.-A. et al. Fabrication of SiCN MEMS by photopolymerization of pre-ceramic polymer. *Sens. Actuat. A Phys.* **95**, 120–134 (2002).

32. Chu, C.-J., Soraru, G. D., Babonneau, F. & Mackenzie, J. D. Preparation and characterization of amorphous SiC film by a liquid route. In *Amorphous and Crystalline Silicon Carbide II* 66–71 (Springer, 1989).

33. Liew, L.-A., Bright, V. M. & Raj, R. A novel micro glow plug fabricated from polymer-derived ceramics: in situ measurement of high-temperature properties and application to ultrahigh-temperature ignition. *Sens. Actuators A Phys.* **104**, 246–262 (2003).

34. Riedel, R. et al. Piezoresistive effect in SiOC ceramics for integrated pressure sensors. *J. Am. Ceram. Soc.* **93**, 920–924 (2010).

35. Mukherjee, S., Ren, Z. & Singh, G. Molecular polymer-derived ceramics for applications in electrochemical energy storage devices. *J. Phys. D. Appl. Phys.* **51**, 463001 (2018).

36. David, L., Bhandavat, R., Barrera, U. & Singh, G. Silicon oxycarbide glass-graphene composite paper electrode for long-cycle lithium-ion batteries. *Nat. Commun.* **7**, 10998 (2016).

37. Bhandavat, R. & Singh, G. Improved electrochemical capacity of precursor-derived Si (B) CN-carbon nanotube composite as Li-ion battery anode. *ACS Appl. Mater. Interfaces* **4**, 5092–5097 (2012).

38. Kaspar, J., Mera, G., Nowak, A. P., Graczyk-Zajac, M. & Riedel, R. Electrochemical study of lithium insertion into carbon-rich polymer-derived silicon carbonitride ceramics. *Electrochim. Acta* **56**, 174–182 (2010).

39. Feng, Y. et al. Nowotny phase Mo3+ 2xSi3C0.6 dispersed in a porous SiC/C matrix: a novel catalyst for hydrogen evolution reaction. *J. Am. Ceram. Soc.* **103**, 508–519 (2019).

40. Jiao, Y., Zheng, Y., Davey, K. & Qiao, S.-Z. Activity origin and catalyst design principles for electrocatalytic hydrogen evolution on heteroatom-doped graphene. *Nat. Energy* **1**, 16130 (2016).

41. Viard, A. et al. Molecular chemistry and engineering of Boron-modified polyorganosilazanes as new processable and functional SiBCN precursors. *Chem. Eur. J.* **23**, 9076–9090 (2017).

42. Lücke, J., Hacker, J., Suttor, D. & Ziegler, G. Synthesis and characterization of silazane- based polymers as precursors for ceramic. *Appl. Organomet. Chem.* **11**, 181–194 (1997).

43. Seifollahi Bazarjani, M. et al. Nanoporous silicon oxycarbonitride ceramics derived from polysilazanes in situ modified with nickel nanoparticles. *Chem. Mater.* **23**, 4112–4123 (2011).

44. Li, Y. et al. Thermal cross-linking and pyrolytic conversion of poly(ureamethylvinyl)silazanes to silicon-based ceramics. *Appl. Organomet. Chem.* **15**, 820–832 (2001).
45. Silva, R., Voiry, D., Chhowalla, M. & Asefa, T. Efficient metal-free electrocatalysts for oxygen reduction: polyaniline-derived N-and O-doped mesoporous carbons. *J. Am. Chem. Soc.* **135**, 7823–7826 (2013).

46. Lu, Z. *et al.* High-efficiency oxygen reduction to hydrogen peroxide catalysed by oxidized carbon materials. *Nat. Catal.* **1**, 156–162 (2018).

47. Dave, K., Park, K. H. & Dhayal, M. Two-step process for programmable removal of oxygen functionalities of graphene oxide: functional, structural and electrical characteristics. *RSC Adv.* **5**, 95657–95665 (2015).

48. Khalid, M. & Varela, H. A general potentiodynamic approach for red phosphorus and sulfur nanodot incorporation on reduced graphene oxide sheets: metal-free and binder-free electrodes for supercapacitor and hydrogen evolution activities. *J. Mater. Chem. A* **6**, 3141–3150 (2018).

49. Li, Y. *et al.* MoS2 nanoparticles grown on graphene: an advanced catalyst for the hydrogen evolution reaction. *J. Am. Chem. Soc.* **133**, 7296–7299 (2011).

50. Shi, Y.-C. *et al.* One-step hydrothermal synthesis of three-dimensional nitrogen-doped reduced graphene oxide hydrogels anchored PtPd alloyed nanoparticles for ethylene glycol oxidation and hydrogen evolution reactions. *Electrochim. Acta* **293**, 504–513 (2019).

51. Chen, H.-Y. *et al.* Simple fabrication of trimetallic platinum-nickel-cobalt hollow alloyed 3D multipods for highly boosted hydrogen evolution reaction. *J. Colloid Interface Sci.* **570**, 205–211 (2020).

52. Popczun, E. J. *et al.* Nanostructured nickel phosphide as an electrocatalyst for the hydrogen evolution reaction. *J. Am. Chem. Soc.* **135**, 9267–9270 (2013).

53. Voiry, D., Shin, H. S., Loh, K. P. & Chhowalla, M. Low-dimensional catalysts for hydrogen evolution and CO 2 reduction. *Nat. Rev. Chem.* **2**, 0105 (2018).

54. Lian, K., Thorpe, S. J. & Kirk, D. W. The electrocatalytic activity of amorphous and crystalline Ni-Co alloys on the oxygen evolution reaction. *Electrochim. Acta* **37**, 169–175 (1992).

55. Morales-Guio, C. G. & Hu, X. Amorphous molybdenum sulfides as hydrogen evolution catalysts. *Acc. Chem. Res.* **47**, 2671–2681 (2014).

56. McEnaney, J. M. *et al.* Electrocatalytic hydrogen evolution using amorphous tungsten phosphide nanoparticles. *Chem. Commun.* **50**, 11026–11028 (2014).

57. Riedel, R., Mera, G., Hauser, R. & Klonczynski, A. Silicon-based polymer-derived ceramics: synthesis properties and applications-a review dedicated to Prof. Dr. Fritz Aldinger on the occasion of his 65th birthday. *J. Ceram. Soc. Jpn.* **114**, 425–444 (2006).

58. Ferrari, A. C. *et al.* Raman spectrum of graphene and graphene layers. *Phys. Rev. Lett.* **97**, 187401–187404 (2006).

59. Hummers, W. S. Jr & Offeman, R. E. Preparation of graphitic oxide. *J. Am. Chem. Soc.* **80**, 1339 (1958).

60. Bazarjani, M. S. *et al.* In situ formation of tungsten oxycarbide, tungsten carbide and tungsten nitride nanoparticles in micro-and mesoporous polymer-derived ceramics. *J. Mater. Chem. A* **2**, 10454–10464 (2014).

61. Massiot, D. *et al.* Modelling one- and two-dimensional solid-state NMR spectra. *Magn. Reson. Chem.* **40**, 70–76 (2002).

62. Lu, X. F., Yu, L. & Lou, X. W. D. Highly crystalline Ni-doped FeP/carbon hollow nanorods as all-pH efficient and durable hydrogen evolving electrocatalysts. *Sci. Adv.* **5**, eaav6009 (2019).

63. Kim, J.-Y. *et al.* Predictive fabrication of Ni phosphide embedded in carbon nanofibers as active and stable electrocatalysts. *J. Mater. Chem. A* **7**, 7451–7458 (2019).

64. Liang, H. *et al.* Plasma-assisted synthesis of NiCoP for efficient overall water splitting. *Nano Lett.* **16**, 7718–7725 (2016).

65. Yin, J. *et al.* Ni–C–N nanosheets as catalyst for hydrogen evolution reaction. *J. Am. Chem. Soc.* **138**, 14546–14549 (2016).

66. Gong, M., Wang, D.-Y., Chen, C.-C., Hwang, B.-J. & Dai, H. A mini review on nickel-based electrocatalysts for alkaline hydrogen evolution reaction. *Nano Res.* **9**, 28–46 (2016).

67. Costa, J. D. *et al.* Electrocatalytic performance and stability of nanostructured Fe–Ni pyrite-type diphosphide catalyst supported on carbon paper. *J. Phys. Chem. C* **120**, 16537–16544 (2016).

68. Wu, Z.-Y. *et al.* Transition metal-assisted carbonization of small organic molecules toward functional carbon materials. *Sci. Adv.* **4**, eaat0788 (2018).

69. Sun, T. *et al.* Single-atomic cobalt sites embedded in hierarchically ordered porous nitrogen-doped carbon as a superior bifunctional electrocatalyst. *Proc. Natl. Acad. Sci.* **115**, 12692–12697 (2018).

70. Jiang, B. *et al.* Ir/g-C3N4/nitrogen-doped graphene nanocomposites as bifunctional electrocatalysts for overall water splitting in acidic electrolytes. *ACS Appl. Mater. Interfaces* **10**, 39161–39167 (2018).

71. Zhang, L. *et al.* Atomic layer deposited Pt-Ru dual-metal dimers and identifying their active sites for hydrogen evolution reaction. *Nat. Commun.* **10**, 1–11 (2019).

72. Luo, Z. *et al.* Reactant friendly hydrogen evolution interface based on di-anionic MoS2 surface. *Nat. Commun.* **11**, 1–9 (2020).

73. Benson, E. E. *et al.* Balancing the hydrogen evolution reaction, surface energetics, and stability of metallic MoS2 nanosheets via covalent functionalization. *J. Am. Chem. Soc.* **140**, 441–450 (2018).

74. Tahira, A. *et al.* Advanced Electrocatalysts for Hydrogen Evolution Reaction Based on Core-Shell MoS2/TiO2 Nanostructures in Acidic and Alkaline Media. *ACS Appl. Energy Mater.* **2**, 2053–2062 (2019).

75. Sarma, P. V, Tiwary, C. S., Radhakrishnan, S., Ajayan, P. M. & Shaijumon, M. M. Oxygen incorporated WS 2 nanoclusters with superior electrocatalytic properties for hydrogen evolution reaction. *Nanoscale* **10**, 9516–9524 (2018).

76. Henckel, D. A., Lenz, O. M., Krishnan, K. M. & Cossairt, B. M. Improved HER catalysis through facile, aqueous electrochemical activation of nanoscale WSe2. *Nano Lett.* **18**, 2329–2335 (2018).

## Acknowledgements

C.S. acknowledges financial supports from CNRS Cellule Energie exploratory project “CeraMicroPac” and Axe Transverse Membranes du Futur project “CeraMesoMem”. D.V. acknowledges the ERC under the European Union’s Horizon 2020 research and innovation programme (Grant No. 804320). II and EC acknowledge Can-BioSe (Grant No. 778157) and the National Science Center of Poland (UMO-2019/35 /B/ST5/00248). C.G. acknowledges The French Région Ile de France—SESAME program for support for the NMR measurements (700 MHz NMR spectrometer). C.S. and P.M. acknowledge the Ecole Nationale Supérieure de Chimie de Montpellier (ENSCM) through the Doctoral School “Ecole doctorale des Sciences Chimiques Balard, ED 459” for the PhD scholarship of Q.H. Didier Cot and Bertrand Rebière are acknowledged for support with SEM and EDX. Valérie Flaud is acknowledged for the XPS measurements.

## Author contributions

C.S. conceived the project, designed the research, supervised the experiments, analyzed the data and wrote the manuscript. Q.H. performed the experiments and analyzed the data. M.Bo. performed the electrical conductivity measurements. D.V. supervised the electrochemical measurements and analyzed the related data. J.B. and V.H. performed the mechanical measurements. I.I. and E.C. performed the TEM measurements. C.G. carried out the
NMR spectroscopy measurements. P.M. supervised the work. M.Be. D.V. and P.M. discussed the results with C.S and Q.H. All of the authors edited the manuscript before submission.

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-020-78558-x.

**Correspondence** and requests for materials should be addressed to P.M. or C.S.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_33.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits 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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020