Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Nature Communications volume 13, Article number: 4602 (2022 ) Cite this article Nickel Boride
Designing efficient catalysts and understanding the underlying mechanisms for anodic nucleophile electrooxidation are central to the advancement of electrochemically-driven technologies. Here, a heterostructure of nickel boride/nickel catalyst is developed to enable methanol electrooxidation into formate with a Faradaic efficiency of nearly 100%. Operando electrochemical impedance spectroscopy and in situ Raman spectroscopy are applied to understand the influence of methanol concentration in the methanol oxidation reaction. High concentrations of methanol inhibit the phase transition of the electrocatalyst to high-valent electro-oxidation products, and electrophilic oxygen species (O* or OH*) formed on the electrocatalyst are considered to be the catalytically active species. Additional mechanistic investigation with density functional theory calculations reveals that the potential-determining step, the formation of *CH2O, occurs most favorably on the nickel boride/nickel heterostructure rather than on nickel boride and nickel. These results are highly instructive for the study of other nucleophile-based approaches to electrooxidation reactions and organic electrosynthesis.
Wensheng Fang, Wei Guo, … Bao Yu Xia
Rongli Liu, Ziyang Wei, … Xiangfeng Duan
Martin Keilwerth, Weiqing Mao, …Karsten Meyer
Electrosynthesis powered by renewable electricity has emerged as a potentially sustainable and environmentally benign route for the production of chemicals and fuels1,2. In most cases, cathodic electrochemical reduction reactions, typically including hydrogen evolution3,4,5, O26,7 and CO2 electroreduction8,9 as well as organic electrocatalytic hydrogenation10, are paired with anodic oxygen evolution reaction (OER). The O2 product of OER holds little economic value11. Besides, thermodynamic analysis suggests that approximately 90% of the electricity driving these above electrosynthesis processes are consumed by the OER12. In other words, anodic OER is main efficiency bottleneck for these processes. In addition to accelerate the OER, electrooxidation reactions of organics such as methanol13, urea14, 5-hydroxymethylfurfural15,16, furfural17 and glucose18 have been proposed to replace sluggish OER. These organic electrooxidation reactions possess lower onset potentials than OER, thereby reducing the cost of electrolysis. Meanwhile, value-added products can also be obtained at the anode17; or the purpose of degrading wastewater can be achieved19.
These organic electrooxidation reactions are also called as nucleophile oxidation reactions (NOR)20. During the reactions, nucleophilic groups containing active hydrogen atoms in these organic substrates undergo dehydrogenation processes to form oxygen-rich products. Nickel-based electrocatalysts have been demonstrated as potential candidates for NOR21,22. However, much ambiguity is present regarding the mechanism of NOR. A typical work reported by Wang’s group proposed that NOR is a two-step, one-electron reaction including an electrogenerated catalyst dehydrogenation reaction from lattice hydroxyl to electron-deficient lattice oxygen and a spontaneous nucleophile dehydrogenation reaction on model Ni(OH)2 catalysts20. Compared to Ni(OH)2, some recent studies have shown that nickel non-hydroxides exhibit better performance23,24,25. Nevertheless, the surface species transformation of metallic catalysts often occurs during electrooxidation, leading to much controversy regarding the activity origin20. Some researches suggested that the oxyhydroxides in situ generated from the non-(hydro)oxide pre-catalysts are identified as a catalytically active phase26,27. Another statement is that the NOR proceeds via adsorbate oxidation mechanism pathway28. Designing efficient electrocatalyst and clarifying its NOR mechanism are essential for the development of nucleophile electrooxidation.
In this work, we focused on a typical nucleophile, methanol, in consideration of the advantages of its simple structure and high stability under alkaline condition. A heterostructure of nickel boride/nickel catalyst was successfully developed for methanol electrooxidation reaction (MOR), delivering nearly 100% formate Faradaic efficiency and long-term stability. Operando electrochemical impedance spectroscopy (EIS) was carried out to track the MOR and OER processes. It was revealed that the OER occurred on the oxide layer formed after catalyst electrooxidation, while the MOR process was not accompanied by observable catalyst electrooxidation when methanol concentration is high, which were also evidenced by in situ Raman spectroscopy. The electrooxidation of catalyst could completely be suppressed in high concentrations of methanol. The electrophilic OH* or O* was produced as the potential increases, and immediately the methanol molecule underwent a spontaneous nucleophile dehydrogenation. At low methanol concentrations, the high-valent products formed by catalyst electrooxidation are catalytically active species for MOR. The proposed mechanism well explained the effect of the presence of methanol on the selectivity towards the anodic reactions, catalyst passivation and catalyst electrooxidation. Furthermore, in situ attenuated total reflection surface enhanced infrared absorption spectroscopy (ATR-SEIRAS) measurement was performed to probe reaction intermediates and identify the methanol electrooxidation pathway. Density functional theory (DFT) calculations were also employed to understand the reasons for efficient methanol electrooxidation on the nickel boride/nickel heterostructure.
The nickel boride/nickel catalyst was prepared by reacting NaBH4 with nickel nitrate, followed by annealing in H2-Ar atmosphere (Fig. 1a, see details in methods)29. The NaBH4 acted as a reducing agent to reduce the nickel ion to form nickel boride (NiBx). The nickel nitrate solution was changed from green to black after the addition of NaBH4, indicating that Ni2+ was reduced. NiBx was used as a precursor, which was further transformed into nickel boride/nickel (NiBx/Ni) after annealing at different temperature in H2-Ar atmosphere. The discussion hereafter mainly focused on the sample annealed at 400 °C (denoted as NiB-400).
a Schematic diagram of the catalyst preparation. b SEM image of the prepared NiB-400 catalyst. c TEM image of the prepared NiB-400 catalyst. d STEM-EDX images of NiB-400 catalyst. e HRTEM image of the prepared NiB-400 catalyst. f–i Local enlargement of the four regions in e.
X-ray powder diffraction (XRD) was conducted to investigate the phase composition of the samples (Supplementary Fig. 1–4). The XRD pattern of NiB-400 is shown in Supplementary Fig. 1. The main peaks are assigned to orthorhombic Ni3B and face-centered cubic Ni. Weak peaks assigned to cubic B2O3 could also be observed. Unannealed precursor (NiBx) is amorphous, and the crystalline phases appear when annealing temperature is above 300 °C (Supplementary Fig. 2). Inductive Coupled Plasma Emission Spectrometer was utilized to determine the elemental content of NiBx, and the result shows that the content ratio of Ni to B is about 2.56:1 (Supplementary Table 1). The surface chemical state was investigated by X-ray photoelectron spectroscopy (XPS). As shown in Ni 2p spectrum of NiB-400 (Supplementary Fig. 5a), the peak at 852.8 eV is attributed to Ni0 in metal nickel and Ni3B. The peak at 856.9 eV is assigned to oxidized Ni. The peaks at 188.2 eV and 192.9 eV in the B 1 s spectrum were attributed to boride and oxidized boron, respectively (Supplementary Fig. 5b). XPS spectra of as-prepared Ni3B and Ni are shown in Supplementary Fig. 6, 7. Scanning electron microscopy (SEM) image (Fig. 1b) and transmission electron microscopy (TEM) images (Fig. 1c, Supplementary Fig. 8, 9) show that the as-prepared nickel boride/nickel, nickel boride and nickel exhibit nanoparticle morphology. Scanning TEM-energy dispersive X-ray (STEM-EDX) shows the homogeneous distribution of Ni element and B element (Fig. 1d). STEM-EDX and high-resolution transmission electron microscopy (HRTEM) reveal the existence of numerous interfaces in NiB-400 (Fig. 1d, e), which is consistent with previous report29. Figure 1f–i are enlarged views of the four regions in Fig. 1e. The interplanar spacings of the four regions are measured to be 0.203, 0.197, 0.176 and 0.186 nm, corresponding to the (111) crystal plane of fcc Ni, the (031) crystal plane of orthorhombic Ni3B, the (200) crystal plane of fcc Ni and the (221) plane of orthorhombic Ni3B (Supplementary Fig. 10).
Figure 2a shows the cyclic voltammetry (CV) curves obtained in 1 M KOH electrolyte and 1 M KOH + 1 M methanol electrolyte. When methanol is added to the KOH electrolyte, the onset potential of the oxidation reaction is significantly reduced, and the potential required for a current density of 10 mA cm−2 is reduced by about 235 mV. By comparing the potentials required to achieve the same current densities for samples annealed at different temperatures, it is found that the sample annealed at 400 °C exhibits the best performance (Fig. 2b). Compared with the recently reported Ni-based electrocatalysts, NiB-400 exhibited superior performance for MOR (Supplementary Table 2). The electrochemical double-layer capacitance (Cdl) of samples annealed at different temperatures shows a trend of decreasing with increasing annealing temperature, indicating the high intrinsic activity of NiB-400 (Supplementary Fig. 11). Combined with the XRD patterns (Supplementary Fig. 2), it can be seen that the improvement of performance is synchronized with the generation of Ni and Ni3B crystalline phases after annealing. With the further increase of the annealing temperature, the content of metal Ni increases due to the transformation of Ni3B into Ni29, and the performance decreases accordingly. Considering that the XRD pattern showing that the material contains a small amount of B2O3 (Supplementary Fig. 1), the methanol oxidation performance of pure B2O3 and clean glassy carbon electrode were tested. The pure B2O3 exhibits the same cubic phase as the B2O3 in NiB-400 (Supplementary Fig. 12). The performance of pure B2O3 and clean glassy carbon on methanol electrooxidation is extremely poor (Supplementary Fig. 13). Further, the performance of Ni3B and metal Ni on methanol oxidation was tested. The performance of Ni3B and metal Ni is far inferior to that of NiB-400 (Fig. 2b, Supplementary Fig. 14), which indicates that the catalytic activity may originate from the nickel boride/nickel heterostructure rather than pure Ni3B or Ni. Faradaic efficiency (FE) of formate obtained after 1 h of chronoamperometry (CA) test at 1.41–1.71 V in 1 M KOH + 1 M methanol electrolyte is very close to 100%, indicating high selectivity to formate product (Fig. 2c, Supplementary Fig. 15). During the continuous reaction at 1.61 V for 5 h, the Faradaic efficiency of formate is maintained around 100% (Fig. 2d). Chronopotentiometry (CP) measurement shows that an additional potential of only 0.07 V, above the initial 1.49 V, is required to maintain a current density of 100 mA cm−2 for 24 h (Fig. 2e). After compensating for the voltage drop, a large current density of 500 mA cm−2 can be achieved at 1.54 V using NiB-400 as the MOR electrocatalyst (Supplementary Fig. 16), which meets the requirement of high current density for potentially industrial applications. After MOR, NiB-400 maintains its original morphology and crystalline phase (Supplementary Fig. 17). The Ni 2p XPS spectra of NiB-400 after OER and MOR are shown in Supplementary Fig. 18.
a CV curves of NiB-400 in 1.0 M KOH electrolyte with and without 1 M methanol. b Comparison of the applied potentials required to achieve varied current densities in 1.0 M KOH with 1 M MeOH electrolyte for samples annealed at different temperatures. c Faradaic efficiency of formate obtained after 1 h of electrolysis at different potentials in 1.0 M KOH electrolyte with 1 M MeOH. d Faradaic efficiency of formate obtained after 5 h continuous electrolysis at 1.61 V in 1.0 M KOH electrolyte with 1 M MeOH. e CP profile at a current density of 100 mA cm−2. Source data are provided as a Source Data file.
Operando electrochemical impedance spectroscopy (EIS) recorded at various potentials was utilized to track the electrochemical processes of MOR and OER. In 1 M KOH electrolyte (Fig. 3a), the response appearing in the high frequency region after 1.36 V is associated with the oxidation of the electrode inner30. The response appearing in the low frequency region is related to the nonhomogeneous charge distribution caused by surface oxidized species30,31. It can be seen that the low-frequency interface appears with the electrooxidation of the catalyst, indicating that the low-frequency interface is between the diffuse double layer (DDL) and the oxidized layer formed by the electrooxidation of the catalyst. OER occurs on this oxidized layer formed by the electrooxidation of the electrocatalyst20. Correspondingly, the Nyquist plots show two semicircles after 1.36 V (Fig. 3b). The semicircles in the high and low frequency regions correspond to the electrooxidation of the catalyst and the OER, respectively. The equivalent circuit models for the operando EIS measured in 1 M KOH are shown in Fig. 2c. Before 1.36 V, the equivalent resistance R1 is large, indicating that the charge transfer is extremely weak. After reaching 1.36 V, R1 suddenly decreases significantly, indicating the onset of catalyst electrooxidation, while the equivalent resistance R2, reflecting the OER, appears and decreases with increasing potential. However, in 1 M KOH electrolyte containing 1 M methanol (Fig. 3e), only one response in the frequency region around 100 Hz is displayed in the Bode plot after 1.36 V. No phase angle peak in low frequency and higher frequency indicates the difference in the reaction interface and no electrooxidation of the catalyst. The MOR in 1 M KOH + 1 M methanol electrolyte does not occur on the surface of oxidized layer (NiOOH, NiOx) formed by the electrooxidation of electrocatalyst, but directly occurs between the DDL and the surface of the catalyst. In the Nyquist plot obtained in 1 M KOH containing 1 M methanol electrolyte (Fig. 3f), only a semicircle corresponding to the MOR reaction is shown after 1.36 V. Interestingly, Nyquist plot shows that the real part at low frequency contracted at 1.36 V, which is due to the adsorption of electroactive species. The equivalent circuit and equivalent resistance of the MOR process are shown in Fig. 3g–h. Both in 1 M KOH electrolyte and 1 M KOH + 1 M methanol electrolyte, the equivalent resistance begins to drop sharply at 1.36 V (Fig. 3d, h, Supplementary Table 3-4), suggesting that MOR and OER may involve some of the same intermediate species. However, the same intermediate species of the two reactions subsequently undergo different changes, resulting in differences in the reaction interfaces of MOR and OER.
a Bode plots and b Nyquist plots of NiB-400 for OER in different potentials. c Equivalent circuit models for OER. d Correlation of the equivalent resistances (R1 and R2) and potentials for NiB-400 during OER. e Bode plots. f Nyquist plots of NiB-400 for MOR in different potentials. g Equivalent circuit models for MOR. h Correlation of the equivalent resistance and potentials for NiB-400 during MOR. In situ Raman spectra of NiB-400 in OER i and in MOR j processes. Source data are provided as a Source Data file.
In situ Raman spectroscopy was used to investigate the changes in catalyst during the OER and MOR processes (Fig. 3i–j). The ex situ Raman spectrum of 1 M KOH + 1 M MeOH electrolyte is shown in Supplementary Fig. 19. In the OER process (1 M KOH electrolyte), two signal peaks appear at ~473 cm−1 and ~556 cm−1 at 1.35 V, corresponding to the bending vibration and stretching vibration of Ni3+-O of γ-NiOOH, respectively (Fig. 3i). At high potentials, the broad peak observed in the 800–1150 cm−1 wavenumber region is attributed to the formation of NiOO−32,33. However, no such features can be found in the MOR process (1 M KOH + 1 M MeOH electrolyte), which means no phase transformation to NiOOH occurs with the existence of high concentrations of methanol (Fig. 3j). The results of in situ Raman spectroscopy are consistent with the results of operando EIS. Therefore, it may be inaccurate to assume that the MOR is an process catalyzed by NiOOH in the presence of high concentrations of methanol.
Methanol concentration has an effect on the shape of the linear sweep voltammetry (LSV) curve. Many Ni-based materials exhibit the similar characteristic23,34,35, but the reason for this characteristic has not been further investigated. The change of the shape of LSV curve may reflect the change of reaction process or even reaction mechanism. For this reason, the LSV curves, Bode plots of operando EIS and in situ Raman spectra measured at different methanol concentrations are compared in Fig. 4a–d. The presence of the oxidation peak of NiII to NiIII in the LSV curve obtained in 1 M KOH + 0.1 M MeOH electrolyte indicates the occurrence of catalyst electrooxidation, which is consistent with that in 1 M KOH (Fig. 4a). Correspondingly, the characteristics of catalyst electrooxidation can be observed in the Bode plots of operando EIS (Fig. 4b) and in situ Raman spectra (Fig. 4c, Supplementary Fig. 20). The current then rises after 1.36 V, indicating the onset of methanol oxidation. However, unlike the case of 1 M MeOH, the response frequency of MOR in the Bode plots is in the low frequency region in the presence of 0.1 M MeOH (Fig. 4b middle), exhibiting similar characteristic to that of OER. This indicates that when the methanol concentration is low, the MOR occurs on the oxide layer formed by the electrooxidation of the catalyst. As shown in Fig. 4e and Supplementary Fig. 21, the response frequency of the reaction in the Bode plots shifts to the higher frequency with the increase of methanol concentration and stabilizes after the methanol concentration reaches 0.5 M, remaining at ~100 Hz. This indicates the change of the catalytic reaction interface from the low frequency interface to the high frequency interface. In addition, the current collapse starts to appear in the LSV curve as the potential increases. As can be seen from Fig. 4b, d, the passivation disappears with further increase of methanol concentration, and it can be found that the passivation phenomenon and the oxidation peak are disappearing simultaneously, which means that the passivation phenomenon is associated with the electrooxidation of the catalyst. The small decrease in current appears in LSV curves after the methanol concentration exceeds 0.3 M stems from the increase in solution resistance (Supplementary Fig. 22). In the Bode plots of EIS, the passivation phenomenon is usually manifested by the abnormal change of the peak position and peak height of the phase angle peak, obvious signal fluctuation or phase angle exceeding 90 degrees, etc. The reason for the passivation phenomenon is that species with poor MOR activity are generated20 due to the oxidation of the electrocatalyst and OER begins to compete with MOR36 as the potential increases. As shown in Fig. 4f, after applying a potential higher than 1.36 V, methanol was injected into 1 M KOH and kept the open circuit for 60 s. When a potential of 1.01 V is applied, the reduction current is significantly lower (red line) than that without methanol injection (black line). This indicates that methanol can react spontaneously and rapidly with the NiOOH without dependence on potential, leading to the consumption of NiOOH and thus the reduction current is decreased. For the case where methanol is present throughout (blue line), no reduction current is present when a potential of 1.01 V is applied, indicating the inhibitory effect of methanol on NiOOH formation.
a LSV curves (Scan rate was set as 20 mV/s), b Bode plots of operando EIS and (c) in situ Raman spectra measured in 1 M KOH, 1 M KOH + 0.1 M MeOH and 1 M KOH + 1 M MeOH. d LSV curves in electrolytes with different concentrations of methanol. (Scan rate was set as 20 mV/s). e The response frequency of the reaction at 1.41 V in the Bode plots as a function of methanol concentration. f Periodic electrochemical measurement over NiB-400. Source data are provided as a Source Data file.
All the above-mentioned evidences show that the high concentration of methanol can inhibit the catalyst electrooxidation to NiOOH, making the OER active site lacking, thus reducing the vicious competition between OER and MOR. Therefore, the Faradaic efficiency of formate is increased to nearly 100% with the increase of methanol concentration in the electrolyte (Supplementary Fig. 23).
Given the discussion above, the statement about MOR is catalyzed by the electrooxidation products of the catalyst such as NiOOH may be inaccurate, since the effect of methanol concentration on the structural changes of the catalyst is ignored. In alkaline electrolyte, the surface of the electrocatalyst adsorbs OH− to form OH*, making the surface appear divalent, which is consistent with the XPS result of the reacted electrocatalyst (Supplementary Fig. 18b). Strong OH* adsorption brings about a phase transition of Ni (oxy)hydroxide, but it becomes different when methanol molecules are present. From the LSV curves, the onset potential of MOR is at about 1.36 V. This potential is consistent with the oxidation potential of NiII to NiIII. For this oxidation process, two possible oxidation pathways are proposed here: (1) Ni atom lost electrons leading to electron deficiency of OH* adjacent to it; (2) Dehydrogenation of OH* adjacent to the Ni atom. These two possible pathways lead to the formation of electrophilic OH* and electrophilic O*, respectively. When the concentration of nucleophilic reagent methanol is high (≥ 0.5 M), the electrophilic oxygen species (OH* or O*) can rapidly capture hydrogen atoms from the methanol molecules, which leads to the oxidation of the methanol molecules while the catalyst eventually returns to its initial state (highlighted by dark colors) and the reaction is not accompanied by electrooxidation products of the catalyst such as NiOOH (Fig. 5a). When the methanol concentration is low (Fig. 5b), the electrophilic oxygen species cannot capture hydrogen atoms from methanol molecules in time, leading to the accumulation of O*. The accumulation of O* induces the occurrence of phase transition to form electrooxidation products of the catalyst such as NiOOH, which is consistent with the results of operando EIS (Fig. 4b) and in situ Raman spectroscopy (Fig. 4c, Supplementary Fig. 20). O* that cannot be consumed in time will combine with OH− in the electrolyte and then form NiOO− species that can be observed in the in situ Raman spectra (Fig. 4c, Supplementary Fig. 20). Thus, for the case of low concentrations of methanol, methanol oxidation can be catalyzed by both electrophilic oxygen species and products of catalyst electrooxidation (i.e., the widely proposed electrochemical-chemical process). Since the electrooxidation layer of the catalyst is also catalytically active for OER, OER competes with MOR when the onset potential of the OER is reached, thus leading to the current collapse (Fig. 4a, d). The presence of high concentrations of methanol is not accompanied by the formation of electrooxidation layer, so the OER is inhibited (Fig. 4a, d).
a Mechanisms of methanol oxidation reaction in high methanol concentration case (Upper: Ni atom loses electrons and the reaction involves electrophilic OH*. Lower: OH* on the catalyst surface loses hydrogen atoms and the reaction involves electrophilic O*). b Mechanisms of methanol oxidation reaction in low methanol concentration case ((1) Ni atom loses electrons. (2) OH* loses hydrogen atoms).
In situ ATR-SEIRAS performed in 1 M KOH containing 1 M MeOH was utilized to determine the reaction intermediates to identify the reaction path for assisting subsequent DFT calculations. The schematic diagram of the device used is shown in Fig. 6a. As shown in Fig. 6b, with the onset of the reaction after 1.36 V, some new bands start to appear in the spectra. The band at 1582 cm−1 is assigned to va(COO), the band at 1381 cm−1 to δ(CH) or ρr(COO), and the band at 1352 cm−1 to vs(COO) of formate37,38. The band at 1640 cm−1 can be assigned to water deformation37. No band of bridge-bonded CO(COB) or linearly-bonded CO(COL) can be observed, indicating that no CO intermediate is involved in the MOR process. Besides, the band of CO32− at around 1400 cm−1 cannot be found, indicating that further oxidation after the formation of HCOOH is unfavorable. Therefore, the product of MOR on NiB-400 is formate accompanied by water generation during the process. Based on the results of in situ ATR-SEIRAS, it is proposed that the methanol oxidation pathway on NiB-400 is CH3OH → *CH3OH → *CH3O → *CH2O → *CHO → *HCOOH → HCOOH. Unlike the previously reported pathway where the oxidation product is CO239, there is no further dehydrogenation of *CHO to generate *CO. According to the observed bands of COO and previous report13, it is proposed that the OH− in the electrolyte interacts with *CHO leading to the formation of *HCOOH.
a Schematic diagram of the ATR-SEIRAS set-up. b ATR-SEIRAS spectra during MOR in 1 M KOH containing 1 M methanol electrolyte. c Differential charge density plot of nickel boride/nickel heterostructure. d Gibbs free energy diagram of MOR occurring on Ni3B(001)/Ni(111) heterostructure (upper) and the change in Gibbs free energy between the steps on different surfaces (lower). e Partial density of states (PDOS) of d-projected DOS of Ni active sites for Ni3B(001)/Ni(111) heterostructure, Ni3B and Ni respectively. Source data are provided as a Source Data file.
The catalytically active species for MOR were discussed above, but the reason for the higher activity of nickel boride/nickel heterostructure over pure nickel boride and pure nickel is still poor understood. Therefore, DFT calculations were further performed. Based on the discussion of catalytically active species above, the catalyst surfaces with adsorbed OH were constructed. On the basis of the finite strain theory, the most matched Ni3B(001)/Ni(111) interface was constructed to demonstrate the heterostructure of nickel boride/nickel29,40. In addition, DFT calculations were also performed for other surfaces and DFT-optimized configurations of the favorable adsorbed intermediates for the oxidation of MeOH on different surfaces are depicted in Supplementary Fig. 24. Differential charge density plot (Fig. 6c) and Bader charge analysis (Supplementary Table 5) found that the electrons are transferred from Ni to Ni3B in the nickel boride/nickel interface, which could facilitate the subsequent oxidation reaction. The Gibbs free energy profiles for the MOR process on the various surfaces are illustrated in Fig. 6d and supplementary Fig. 25. It can be seen that the potential-determining step for MOR are all contributed by the process of *CH3O → *CH2O + H+ + e− on the four surfaces of Ni3B(001)/Ni(111) heterostructure, Ni3B(221)/Ni(111) heterostructure, Ni3B and Ni. Interestingly, Ni3B(001)/Ni(111) heterostructure exhibits the lowest Gibbs energy barrier of 0.57 eV in comparison to 0.74 and 0.70 eV for Ni3B and Ni, respectively, indicating that MOR can occur most favorably on the Ni3B/Ni surface. The ΔG of the HCOOH generation is −0.19 eV, which is much lower than that for the further oxidization to *COOH intermediate (−0.01 eV), suggesting that the final product of the MOR on Ni3B/Ni heterostructure is formic acid rather than carbon dioxide. This is also consistent with the in situ ATR-SEIRAS results. In addition, the produced formic acid can readily exist in the form of formate in alkaline medium, spontaneously making it difficult to be further oxidized. The partial density of states (PDOS) for the adsorption behavior of *CH2O intermediate on Ni3B(001)/Ni(111) heterostructure, Ni3B, and Ni structures were also calculated. As shown in Fig. 6e, the d-band centers of Ni3B/Ni heterostructure, Ni3B, and Ni are separetly −1.22, −1.26, and −1.34 eV. Obviously, the value of Ni3B/Ni heterostructure is most closely approaching to the Fermi energy level, which endows the strongest adsorption capacity for *CH2O. As a result, it can effectively promote the binding of key intermediate and lower the associated energy barrier. Therefore, the nickel boride/nickel heterostructure exhibits better activity for methanol oxidation than nickel boride and nickel (Supplementary Fig 14).
Taken together, the methanol oxidation process on nickel boride/nickel and the evolution of the catalyst are proposed (Fig. 7). When the methanol concentration is high, the electrooxidation of the catalyst is inhibited and the electrophilic oxygen species is the catalytically active species. In this case, MOR is most favorably occurring on the nickel boride/nickel heterostructure due to the fact that the heterstructure possesses the lowest ΔG for the potential-determining step of MOR and has the optimal adsorption for the key intermediates *CH2O. But when the methanol concentration is very low, the inhibitory effect of methanol on the electrooxidation of the catalyst becomes weak and the catalyst undergoes phase transition to generate high-valent products such as NiOOH. Methanol molecules are easily oxidized by the oxidative NiOOH. NiOOH is reactive to both OER and MOR, which can lead to a vicious competition. Therefore, increasing the methanol concentration can inhibit the OER.
Schematic representation of the methanol oxidation process on nickel boride/nickel and the evolution of the catalyst.
In summary, we have reported a heterostructure of nickel boride/nickel catalyst for methanol electrooxidation. Only an anodic potential of about 1.51 V has been required to achieve a current density of 100 mA cm−2 on the optimized nickel boride/nickel catalyst. Nearly 100% formate Faradaic efficiency and good long-term stability have also been delivered. Combining electrochemical measurements, operando electrochemical impedance spectroscopy and in situ Raman spectroscopy, the effect of methanol concentration on the mechanism of MOR is presented. The formation of high-valent oxidation products of the electrocatalyst is inhibited at high methanol concentrations, and the electrophilic oxygen species (OH* or O*) is the catalytically active species for MOR. DFT calculations indicate that the heterostructure can effectively promote the binding of key intermediate and lower the associated energy barrier, and thus has superior MOR performance over nickel boride and nickel. When the methanol concentration is low, the inhibitory effect of methanol on catalyst electrooxidation decreases and methanol is easily oxidized by the high-valent electrooxidation products of the electrocatalyst. The proposed mechanism can well explain the effect of the presence of methanol on the selectivity towards the anodic reactions, catalyst passivation and catalyst electrooxidation. The understanding can readily be generalized to other nucleophile electrooxidations and organic electrosynthesis.
Sodium borohydride (NaBH4, ≥ 98%, Weng Jiang Reagent), Nickel(II) nitrate hexahydrate ((Ni(NO3)2·6H2O, 99%, Acros Organics), Methanol (≥ 99.9%, Aladdin), Potassium hydroxide (KOH, 95%, Macklin), Nafion (5 wt. % in mixture of lower aliphatic alcohols and water, contains 45% water), Deuterium oxide(99.9% atom% D, contains 0.05 wt% TMSP, SIGMA-ALDRICH), Diboron trioxide (B2O3, ≥ 98%, SCR),Ethanol (≥99.7%, GENERAL-REAGENT), Tetraethylene glycol (99%, Shanghai Canspec), DI water (18.25 MΩ·cm−1), Nickel (20–100 nm, 99.9%, Macklin).
The morphology was characterized by Field-emission scanning electron microscopy (FESEM, Hitachi, S-4800) with an acceleration voltage of 15 kV. The microstructure observation was performed by Transmission electron spectroscopy (TEM, FEI F30). Phase compositions of the electrocatalysts were examined using a Rigaku Ultima IV operated under a voltage of 40 kV and a current of 40 mA. X-ray photoelectron spectroscopy (XPS) was analyzed by VG Escalab 250xi. All XPS patterns were corrected with C 1 s signal at 284.8 eV as the standard. Raman spectroscopic measurements were performed using a LabRAM HR800 confocal microscope (Horiba Jobin Yvon) (Excitation wavelength: 532 nm; Power: 25%; Acquisition time: 80 s). ICP was performed on Agilent 725. Attenuated Total Reflection Surface-Enhanced Infrared Absorption (ATR-SEIRA) measurements were performed on ThermoFisher Nicolet is50 (dector: MCT/A; number of scans:16; moving mirror speed: 0.4747; ink dosage: 800 μg in 80 μL). The optical crystal used here is a silicon prism, and a gold film is sputtered on its working surface using a sputter coater to enhance the signal (parameters of the sputter coater: current: 20 mA; duration: 250 s). Previous studies have shown that the effect of Au film on the MOR is negligible for the selected potential range23.
Ni(NO3)2·6H2O (1.175 g) was dissolved completely in 100 mL of DI water. NaBH4 solution (0.232 g, dissolved in 10 mL DI water) was added into the solution later and stirred for 2 h. It should be noted that the reaction should be carried out in a cold water bath at around 10 °C. The precipitate was then washed three times with cold DI water (~10 °C) and cold ethanol (~10 °C) to remove impurities. All the water and ethanol used in the synthesis and washing process are ventilated with nitrogen for 15 min to remove oxygen. The obtained black precipitate was vacuum dried at 60 °C. Completely dry precipitate was named NiBx. NiBx was divided into several parts and heated to 300, 350, 400, 450, 500 °C respectively with a heating rate of 5 °C/min and held for 3 h in H2/Ar (v/v = 5/95). The obtained black powders were named NiB-300, NiB-350, NiB-400, NiB-450 and NiB-500, respectively.
The synthetic process of Ni3B referred to the reported literature with some modifications41. In detail, 238 mg of nickel chloride and 30 mL of Tetraethylene glycol (TEG) were added into the 150 mL Round-bottom flask. The solution was sonicated for 60 min until most of the salts were dissolved and then allowed to stir and heat to 45 °C on an oil bath pan while purging with Ar all the time. Next, 30 mL of TEG with 1 g of NaBH4 was added to the metal salt solution within 1 min. The color of the solution quickly turned black. The temperature was slowly warmed up to 280 °C and maintained for 5 min with continuous stirring under an Ar purge. Lastly, after being naturally cooled down, the product was collected by centrifugation with ethanol for several times.
Ni(NO3)2·6H2O (1.175 g) was dissolved completely in 100 mL of DI water. NaBH4 solution (0.232 g, dissolved in 10 mL DI water) was added into the solution later and stirred for 2 hours. The reaction was carried out at around 28 °C. The precipitate was then washed three times with DI water and ethanol to remove impurities. The obtained gray precipitate was vacuum dried at 60 °C. The dried precipitate was then heated to 400 °C with a heating rate of 5 °C/min and held for 3 hours in H2/Ar (v/v = 5/95).
The electrochemical tests were performed in a three-electrode system (25 ± 1 °C) on CHI 760E electrochemical workstation. The electrolyte was 1 M KOH aqueous solution for OER and 1 M KOH (containing 1 M methanol) for MOR if no special instruction. To prepare the working electrode, 5 mg catalysts were mixed with 460 μL ethanol and 40 μL Nafion solution with ultrasonication. 5 μL of the obtained suspension was dropped onto a φ5 mm glassy carbon electrode to be used as the working electrode. Hg/HgO (1 M KOH) and graphite rod were used as a reference electrode and a counter electrode, respectively. Unless otherwise specified, the scanning rate of LSV and CV was set as 5 mV/s. Electrochemical impedance spectroscopy (EIS) was performed over a frequency range from 0.01 Hz to 100000 Hz and AC amplitude was set as 5 mV. Electrochemical data was not adjusted by iR compensation if no special instruction is given. All potentials were calibrated to the reversible hydrogen electrode (RHE) reference scale using the formulas as follows:
Electrochemical double-layer capacitance (Cdl) was measured by cyclic voltammetry. The potential window was selected in a range of OCP ± 50 mV. The scan speed was set to 5, 10, 15, 20, 25, 30, 35, 40, 45 mV/s. The Δj at the open circuit potential was plotted against the scan rate and the slope fitted was the estimated value of Cdl.
For the OER process, the equivalent circuit before 1.35 V consists of a solution resistor (Rs) in series with R1 | | CPE1, which reflects the electrooxidation of the catalyst. After 1.35 V, the equivalent circuit consists of a solution resistor, an R1 | | CPE1 and R2 | | CPE2 in series. A new R2 | | CPE2 appears to reflect the OER process, and R2 is the charge transfer resistance of the OER. For the MOR process (1 M KOH + 1 M methanol), since the electrooxidation of the catalyst is suppressed, the equivalent circuit is composed of a solution resistance (Rs) and R1 | | CPE1 in series, and R1 | | CPE1 reflects the MOR process.
Often a CPE is used in a model in place of a capacitor to compensate for non-homogeneity in the system. The CPE is defined by two values, CPE-T and CPE-P. If CPE-P equals 1, then CPE is equivalent to a capacitor. The contributions of the capacitive resistances of the EIS results for the OER and MOR processes were fitted according to the above circuits. The fitting results are detailed in Supplementary Table 3, 4.
Nuclear magnetic resonance (Ascend 600) was used to analyze the product content of the methanol oxidation reaction. First, a series of potassium formate (HCOOK) solutions of different concentrations (1 mM, 3 mM, 5 mM, 7 mM, 10 mM, 15 mM, 20 mM) were prepared. 500 μL of potassium formate solution and 100 μL of D2O (containing 0.05 wt% TMSP) were added into a NMR tube and mixed thoroughly for NMR test. A standard curve was drawn based on the proportional relationship between the integrated intensity of the formate ion signal peak and the concentration of formate ion. Typically, for the analysis of the product in the electrolyte, 500 μL of electrolyte and 100 μL of D2O (containing 0.05 wt% TMSP) were added in a NMR tube and mix them thoroughly. Then the nuclear magnetic resonance test was performed to obtain the integrated intensity of the formate ion signal. The content of formate produced was calculated according to the standard curve.
Geometry optimization and electronic structure computation have been carried out by using the Vienna abinitio simulation package (VASP)42. The projector-augmented-wave (PAW) method together with the generalized gradient approximation-Perdew-Burke-Ernzerhof (GGA-PBE) functional were employed in the simulation process43. An energy cutoff of 450 eV was used for the plane-wave expansion of the electronic wave function. All geometric structures were relaxed until the energy and forces were converged to 10−5 eV and −0.02 eV Å−1, respectively. A system of 2 × 2 slab with 4 layers was employed to model the Ni surfaces and a 2 × 2 slab with 1 layer was used to simulate the Ni3B surface. A vacuum layer of 15 Å along the z-direction was applied to separate the surfaces. The geometric model of Ni (111) and Ni3B (001) were chosen to demonstrate the pristine Ni and Ni3B materials. To simulate the Ni3B/Ni heterostructure, we built the layered model to investigate the interaction effect between Ni and Ni3B, as well as the complete mechanism of MOR on such hybrid catalyst. The Brillouin zone was sampled with 3 × 2 × 1 Monkhorst-Pack k-meshes for the geometry optimization, while 6 × 4 × 2 grid of K-points were used for the more accurate electronic states analysis.
The Gibbs free energy change (ΔG) of each reaction step was calculated according to the following formula:
where ΔEads is the adsorption energy of reaction intermediates; ΔEZPE and ΔS are the energy difference in zero point energy and entropy, respectively. Therein, a more negative Eads implies that the adsorption is thermodynamically more favorable.
The data that support the findings of this study are available from the corresponding author upon reasonable request. Source data are provided with this paper.
Blanco, D. & Modestino, M. Organic electrosynthesis for sustainable chemical manufacturing. Trends Chem. 1, 8–10 (2019).
Lam, K., Dobbs, A., Leech, M., Petti, A. & Garcia, A. Organic electrosynthesis: From academia to industry. React. Chem. Eng. 5, 977–990 (2020).
Sultan, S. et al. Single atoms and clusters based nanomaterials for hydrogen evolution, oxygen evolution reactions, and full water splitting. Adv. Energy Mater. 9, 1900624 (2019).
Zhu, Y., Tian, P., Jiang, H., Mu, J. & Li, C. Synergistic effect of platinum single atoms and nanoclusters boosting electrocatalytic hydrogen evolution. CCS Chem. 3, 2539–2547 (2021).
Liu, T. et al. CoP-Doped MOF-based electrocatalyst for pH-universal hydrogen evolution reaction. Angew. Chem. Int. Ed. 58, 4679–4684 (2019).
Zhang, T., Zhang, B. K., Peng, Q., Zhou, J. & Sun, Z. M. Mo2B2 MBene-supported single-atom catalysts as bifunctional HER/OER and OER/ORR electrocatalysts. J. Mater. Chem. A 9, 433–441 (2021).
Liang, Z. et al. Highly curved nanostructures coated Co, N-doped carbon materials for oxygen electrocatalysis. Angew. Chem. Int. Ed. 60, 12759–12764 (2021).
Article ADS CAS Google Scholar
Lee, W. H. et al. Carbon-Supported IrCoOx nanoparticles as an efficient and stable OER electrocatalyst for practicable CO2 electrolysis. Appl. Catal. B 269, 118820 (2020).
Pang, R. et al. Tracking structural evolution: Operando regenerative CeOx/Bi interface structure for high-performance CO2 electroreduction. Natl Sci. Rev. 8, nwaa187 (2021).
Article CAS PubMed Google Scholar
Akhade, S. A. et al. Electrocatalytic hydrogenation of biomass-derived organics: A Review. Chem. Rev. 120, 11370–11419 (2020).
Article CAS PubMed Google Scholar
Sherbo, R., Delima, R., Chiykowski, V., MacLeod, B. & Berlinguette, C. Complete electron economy by pairing electrolysis with hydrogenation. Nat. Catal. 1, 501–507 (2018).
Dong, C. et al. Dechlorination-facilitated deprotonation of CoFe (Oxy)hydroxide catalysts under electrochemical oxygen evolution. Chem. Eng. Sci. 252, 117270 (2022).
Hao, J. et al. In situ facile fabrication of Ni(OH)2 nanosheet arrays for electrocatalytic co-production of formate and hydrogen from methanol in alkaline solution. Appl. Catal. B 281, 119510 (2021).
Wang, X. et al. Hierarchical coral-like NiMoS nanohybrids as highly efficient bifunctional electrocatalysts for overall urea electrolysis. Nano Res. 11, 988–996 (2018).
Yang, G., Jiao, Y., Yan, H., Xie, Y. & Fu, H. Interfacial Engineering of MoO2-FeP Heterojunction for Highly Efficient Hydrogen Evolution Coupled with Biomass Electrooxidation. Adv. Mater. 32, 2000455 (2020).
Lu, Y. X. et al. Identifying the geometric site dependence of spinel oxides for the electrooxidation of 5-Hydroxymethylfurfural. Angew. Chem. Int. Ed. 59, 19215–19221 (2020).
Bharath, G. & Banat, F. High-grade biofuel synthesis from paired electrohydrogenation and electrooxidation of furfural using symmetric Ru/Reduced graphene oxide electrodes. ACS Appl. Mater. Inter. 13, 24643–24653 (2021).
Liu, W. J. et al. Efficient electrochemical production of glucaric acid and H2 via glucose electrolysis. Nat. Commun. 11, 265 (2020).
Article ADS CAS PubMed PubMed Central Google Scholar
Sha, L. et al.In situ grown 3D hierarchical MnCo2O4.5@Ni(OH)2 nanosheet arrays on Ni foam for efficient electrocatalytic urea oxidation.Chem.Eng.J. 381, 122603 (2020).
Chen, W. et al. Activity origins and design principles of nickel-based catalysts for nucleophile electrooxidation. Chemistry 6, 2974–2993 (2020).
Liang, Z. et al. Molecular Engineering to Tune the Ligand Environment of Atomically Dispersed Nickel for Efficient Alcohol Electrochemical Oxidation. Adv. Funct. Mater. 31, 2106349 (2021).
Gao, P. et al. Promoting effect of Bi in Ni–Bi oxide electrocatalysts for methanol oxidation reaction. J. Mater. Sci. Mater. Electron. 31, 13219–13228 (2020).
Li, J. et al. Selective methanol-to-formate electrocatalytic conversion on branched nickel carbide. Angew. Chem. Int. Ed. 59, 20826–20830 (2020).
Wu, T. H. & Hou, B. W. Superior catalytic activity of α-Ni(OH)2 for urea electrolysis. Catal. Sci. Technol. 11, 4294–4300 (2021).
Li, R. Q. et al. Hierarchical Ni3N/Ni0.2Mo0.8N heterostructure nanorods arrays as efficient electrocatalysts for overall water and urea electrolysis. Chem. Eng. J. 409, 128240 (2021).
Chen, S. et al. Methanol electro-oxidation in alkaline on a titanium electrode supporting Ni-P or Ni deposits. Recent Pat. Corros. Sci. 2, 42–47 (2012).
Chen, D. & Minteer, S. D. Mechanistic study of nickel based catalysts for oxygen evolution and methanol oxidation in alkaline medium. J. Power Sources 284, 27–37 (2015).
Article ADS CAS Google Scholar
Zhou, B. et al. Activity origin and alkalinity effect of electrocatalytic biomass oxidation on nickel nitride. J. Energy Chem. 61, 179–185 (2021).
Yang, F. et al. Inter-regulated d-band centers of the Ni3B/Ni heterostructure for boosting hydrogen electrooxidation in alkaline media. Chem. Sci. 11, 12118–12123 (2020).
Article CAS PubMed PubMed Central Google Scholar
Lu, Y. X. et al. Tailoring competitive adsorption sites by oxygen-vacancy on cobalt oxides to enhance the electrooxidation of biomass. Adv. Mater. 34, 2107185 (2022).
Wang, H. Y. et al. In operando identification of geometrical-site-dependent water oxidation activity of spinel Co3O4. J. Am. Chem. Soc. 138, 36–39 (2016).
Article CAS PubMed Google Scholar
Diaz-Morales, O., Ferrus-Suspedra, D. & Koper, M. T. M. The importance of nickel oxyhydroxide deprotonation on its activity towards electrochemical water oxidation. Chem. Sci. 7, 2639–2645 (2016).
Article CAS PubMed PubMed Central Google Scholar
Cole, K. M., Kirk, D. W. & Thorpe, S. J. In situ raman study of amorphous and crystalline Ni-Co alloys for the alkaline oxygen evolution reaction. J. Electrochem. Soc. 165, J3122–J3129 (2018).
Rezaee, S. & Shahrokhian, S. Facile synthesis of petal-like NiCo/NiO-CoO/nanoporous carbon composite based on mixed-metallic MOFs and their application for electrocatalytic oxidation of methanol. Appl. Catal. B 244, 802–813 (2019).
Li, J. et al. Superior methanol electrooxidation performance of (110)-faceted nickel polyhedral nanocrystals. J. Mater. Chem. A 7, 22036–22043 (2019).
Tao, H. B. et al. A general method to probe oxygen evolution intermediates at operating conditions. Joule 3, 1498–1509 (2019).
Kowal, A., Port, S. N. & Nichols, R. J. Nickel hydroxide electrocatalysts for alcohol oxidation reactions: An evaluation by infrared spectroscopy and electrochemical methods. Catal. Today 38, 483–492 (1997).
Christensen, P. A. & Linares-Moya, D. The role of adsorbed formate and oxygen in the oxidation of methanol at a polycrystalline Pt Electrode in 0.1 M KOH: An In Situ Fourier Transform Infrared Study. J. Phys. Chem. C. 114, 1094–1101 (2010).
Wang, X. et al. Materializing efficient methanol oxidation via electron delocalization in nickel hydroxide nanoribbon. Nat. Commun. 11, 4647 (2020).
Article ADS CAS PubMed PubMed Central Google Scholar
Liu, T. et al. Interfacial electron transfer of Ni2P–NiP2 polymorphs inducing enhanced electrochemical properties. Adv. Mater. 30, 1803590 (2018).
Schaefer, Z. L., Ke, X., Schiffer, P. & Schaak, R. E. Direct Solution Synthesis, Reaction Pathway Studies, and Structural Characterization of Crystalline Ni3B Nanoparticles. J. Phys. Chem. C. 112, 19846–19851 (2008).
Kresse, G. & Furthmiiller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. J. Comput. Mater. Sci. 6, 15–50 (1996).
Hammer, B., Hansen, L. B. & Nørskov, J. K. Improved adsorption energetics within density-functional theory using revised Perdew-Burke-Ernzerhof functionals. Phys. Rev. B 59, 7413–7421 (1999).
This work was supported by the National Natural Science Foundation of China (21838003, 91834301, 21978278, 21803063 and 21978087), the Shanghai Scientific and Technological Innovation Project (18JC1410500 and 19JC1410400), and the Fundamental Research Funds for the Central Universities (222201718002). The authors thank the Shanghai synchrotron Radiation Facility (14W1, SSRF), the Beijing Synchrotron Radiation Facility (1W1B and soft-X-ray endstation, BSRF), the Hefei Synchrotron Radiation Facility (Photoemission, Magnetic Circular Dichroism and Catalysis/Surface Science Endstations, Endstations at National Synchrotron Radiation Laboratory).
These authors contributed equally: Yanbin Qi, Yue Zhang.
Key Laboratory for Ultrafine Materials of Ministry of Education, School of Chemical Engineering, East China University of Science and Technology, Shanghai, 200237, China
Yanbin Qi, Hongliang Jiang & Chunzhong Li
Shanghai Engineering Research Center of Hierarchical Nanomaterials, School of Materials Science and Engineering, East China University of Science and Technology, Shanghai, 200237, China
Yanbin Qi, Yuhan Zhao, Yihua Zhu & Chunzhong Li
Institutes of Physical Science and Information Technology, Anhui University, Hefei, 230601, China
Yue Zhang & Li Yang
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
Y.B.Q., H.L.J., and C.Z.L. conceived the idea, designed the experiments and wrote the draft. Y.B.Q. performed main sample synthesis and electrochemical measurements. Y.Z. and L.Y. carried out DFT calculations and analyzed the calculation results. Yu.H.Z performed part of sample synthesis and analysis. Yi.H.Z contributed to the analysis of in situ Raman data.
Correspondence to Li Yang, Hongliang Jiang or Chunzhong Li.
The authors declare no competing interests.
Nature Communications thanks Guofeng Wang, Shuangyin Wang and Hongyan Liang for their contribution to the peer review of this work.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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 license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. 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/.
Qi, Y., Zhang, Y., Yang, L. et al. Insights into the activity of nickel boride/nickel heterostructures for efficient methanol electrooxidation. Nat Commun 13, 4602 (2022). https://doi.org/10.1038/s41467-022-32443-5
DOI: https://doi.org/10.1038/s41467-022-32443-5
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Nature Communications (Nat Commun) ISSN 2041-1723 (online)
Sand Making Machine Raw Material Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.