For the operando microscopy experiments, we prepared well-defined Cu2O cubes on the carbon working electrode of the EC-TEM chips via electro-deposition28,29 as shown schematically in Fig. 1a. The as-prepared Cu2O cubes have an average size of 250 nm and consist of six {100} facets without the exposure of other minor facets such as {110} or {111} (Suppl. Figure 1). All voltages indicated in this paper are referenced against a Ag/AgCl electrode and then converted to the reversible hydrogen electrode scale (RHE) using the Nerst equation and the bulk pH of the electrolyte. Intriguingly, the image sequences show that the cubes do not undergo significant change in the commonly employed 0.1 M Na2SO4 + 8 mM NaNO3 electrolyte for NO3RR (Fig. 1b), during the initial sweep towards cathodic potentials. Cu2O should reduce directly to metallic Cu at the higher overpotentials of the sweep, where according to the Pourbaix diagram, the redox potential for the transformation of Cu2O to metallic Cu is 0 VRHE in a solution with pH 720,21, and metallic Cu is the stable phase below − 0.2 VRHE onwards. The stability of the Cu2O cubes is remarkable because these redox transformations usually lead to morphological changes. For comparison, Fig. 1c depicts a Cu2O cube under CO2RR conditions in CO2-saturated 0.1 M KHCO3 at a potential similar to that applied in the NO3RR experiment. As we reported previously29, the latter cubes undergo fragmentation together with the re-deposition of small particles, a behavior that differs from the morphologically much more stable NO3RR samples at the same applied potentials. The linear sweep voltammograms acquired during these two experiments are available as Suppl. Figure 2.
Next, we studied these Cu2O cubes systematically at different sustained potentials from − 0.2 VRHE to -0.6 VRHE (Fig. 2a-e) to probe further their morphological stability during NO3RR. For these extended experiments, we adopted an intermittent imaging protocol (images captured at 15-minute intervals with the electron beam blanked the rest of the time) to minimize beam-induced dissolution of the Cu2O cubes (see discussion in Suppl. Note 1) and ensure that the catalyst re-structuring kinetics we extract from the collected data are as accurate as possible. The electrochemical current profiles over time at each potential measured in these EC-TEM experiments are provided as Suppl. Figure 3. At -0.2 VRHE (Fig. 2a), the cubes were stable during our entire observation, with no significant restructuring observed. From − 0.2 to -0.5VRHE, dissolution/re-deposition is the main restructuring pathway. At -0.3 VRHE, the cubic form persisted for almost 135 min (Fig. 2b), while the cube completely dissolved after 140 min at -0.4 VRHE (Fig. 2c) and after 90 min at -0.5 VRHE (Fig. 2d). The lighter contrast of the cube exterior in Fig. 2c compared to the middle of the cube at 60 and 80 min is explained by the cube corners and edges being etched first. The weaker contrast of the dissolving Cu2O cubes compared to that of the growing Cu NPs also suggests that the dissolving cubes were still in oxide form.
Two Cu2O cubes were captured in the images acquired at -0.6 VRHE (Fig. 2e). One cube shrank and restructured into a smaller cube with a void in the center and then became rougher due to small NPs attaching to its surface, while another completely dissolved within the same time frame. We also highlight that the intensity of the cubes in the TEM images obtained from − 0.3 to -0.5 VRHE gradually decrease, whereas the intensity of the cube at -0.6 VRHE is brighter, implying that cube-like frame at -0.6 VRHE is metallic. Moreover, the interplay of dissolution/re-deposition and direct reduction at the more cathodic potentials means that the terminal catalysts morphologies of oxide pre-catalysts vary depending on the applied potential.
Next, we repeated the NO3RR experiments in a H-type cell with Cu2O cubes electrodeposited on carbon paper to check for the consistency of the EC-TEM results with standard reaction geometries. Figures 2f and 2g show lower magnification images of samples from the EC-TEM experiments with scanning electron microscopy (SEM) images of samples extracted from H-type cell experiments after 2 hours of reaction at three different applied potentials, -0.2, -0.4 and − 0.6 VRHE. Electron diffraction patterns taken from samples extracted after reaction show that the cubes did not undergo significant restructuring at -0.2 VRHE and remain Cu2O whereas samples reacted at -0.6 VRHE were largely metallic (Suppl. Figure 4). Conversely, samples reacted at -0.4 VRHE show a mixture of residual Cu2O and metallic Cu structures (Suppl. Figure 4–6). The morphological differences between the sample after reaction in the H-type cell at -0.4 VRHE and − 0.6 VRHE further support that the catalysts indeed restructure through different pathways as described by our EC-TEM experiments. Inductively coupled plasma mass spectrometry measurements of the electrode and the electrolyte in the H-type cell after reaction further show that Cu dissolution happens at all the applied potentials (Suppl. Figure 7). Therefore, these experiments indicate that the Cu2O cubes undergo a gradual dissolution under NO3RR conditions, which in turn leads to the re-deposition of metallic particles elsewhere on the working electrode with shapes and sizes that are modulated by the applied potential.
To obtain unambiguously the oxidation state of the catalyst species present during reaction and rule out the possibility that the ex situ identified Cu2O phase are the result of re-oxidation during the return to open circuit potential30,31 (OCP), we performed operando EC-TXM measurements on the Cu2O cubes by transferring our EC-TEM holder into a TXM at the BESSY II synchrotron facility as illustrated in Suppl. Figure 8. This unique arrangement maintains the same reaction environment between the two experiments, while enabling time-resolved operando measurements of Cu absorption edges under applied potential without compromising the sustained electrolyte flow as X-rays are attenuated less by the electrolyte and enclosing membranes. Thus, the evolution of the electrode’s composition can be tracked during NO3RR. Figures 3a-d show the time-resolved evolution of the Cu2O catalysts during NO3RR at -0.4 VRHE as observed by EC-TXM in the form of the colored maps that were reconstructed from a XAS image stack using linear combination fitting (LCF)32. Cu2O is depicted in red and metallic Cu in yellow. The maps show that Cu2O and Cu are the dominant phases present over the duration of the NO3RR and that the oxide and metallic phases coexist under specific reaction conditions, but they are spatially separated. CuO is also detected (blue color) but it is present only in small quantities and is not clearly visible from the maps. The corresponding decomposed spectra are shown in Fig. 3e-p, where the yellow, red and blue colored lines represent the respective Cu species. The total contribution of the individual spectra in Fig. 3e-p represents the amount (thickness) of each species, which indicates that the content of metallic Cu species increases (Fig. 3e-h) during reaction, whereas Cu2O decreases (Fig. 3i-p). The details of the data acquisition and processing are discussed in Suppl. Note 3. Most importantly, these results confirm the sluggish reduction kinetics of the large Cu2O cubes in the Na2SO4 + NaNO3 electrolyte, and that the metallic phase forms when the dissolved Cu species re-deposit on the working electrode due to the reductive potential employed.
We further verified that the slow reduction of the Cu oxide cubes extends to larger reaction volumes with operando hard XAS measurements of samples electro-deposited on carbon paper in our home-built electrochemical XAS cell33. In Fig. 3q, we plot the changes in the Cu K-edge valence states (from 8950 to 9105 eV) that were obtained from operando XAS. The weight of the Cu valence state is extracted by LCF of the X-ray absorption near-edge structure (XANES) of the oxide-derived Cu catalyst collected at a constant potential of -0.4 VRHE in 0.1 M Na2SO4 + 8 mM NaNO3 electrolyte. As seen in Fig. 3r, the fraction of Cu2O decreased but not completely after more than 2 hours of electrolysis, while the fraction of metallic Cu increased correspondingly, eventually to almost a 1:1 ratio of Cu2O:Cu. The XANES results agree with the persistence of Cu2O and the continual evolution of the Cu species seen in the EC-TEM (Fig. 2c) and EC-TXM (Cu-L3 edges in Fig. 3r) results at -0.4 VRHE. Minute amounts of the CuO species were also detected during the experiment. The change in the weights of the three species over time exhibit similar trends in both TXM and XAS, confirming that the results we obtain in the EC-TEM cells indeed extrapolate to a larger ensemble of catalyst particles.
This overall agreement between different methods and experimental geometries means that we can use the in situ TEM image sequences to quantify the potential-dependent dissolution and re-deposition rates. Our method for fraction extraction from the EC-TEM images and additional analysis of the re-deposited particles is described in the Suppl. Note 4. As shown in Fig. 4a, the cube fraction decreases over time at an increasing rate as the potential decreases from − 0.2 V to -0.5 VRHE. The sample at -0.6 VRHE deviates from this trend (dark purple line) because of the direct reduction of Cu2O to metallic Cu. In Fig. 4b, we use the cube dissolution rate to estimate the Cu2O to Cu ratio at a certain potential and use it to visualize the majority phase (greater than 50%) at different times. We further compare the Cu2O/Cu ratio with NH3 conversion activity (current density) and selectivity (Faradaic Efficiency, FE) obtained from our bench top electrochemistry measurements. In Fig. 4c, we plot the linear sweep voltammogram of the Cu2O cubes prepared on carbon paper and in Fig. 4d-e, their product distribution as a function of the applied potentials. The measured yield rate and the FE towards NH3 was much higher at -0.6 VRHE as compared to -0.2 VRHE and − 0.4 VRHE, implying that the change in catalytic selectivity is related to the faster rate of oxide to metal conversion at -0.6 VRHE (Fig. 4b).
Next, we performed EC-TEM studies in various electrolyte compositions to elucidate the mechanism behind Cu2O stabilization. Suppl. Figure 9a and 9b describe experiments using pure 0.1 M Na2SO4 and 0.1 M Na2SO4 + 8 mM NaNO2 (i.e. nitrite reduction) respectively. In both cases, the cubes behaved similarly to their behavior with NO3RR where they gradually reduced in size until they fragment/reduce at longer reaction duration, which means the Cu2O stability is related to the Na2SO4 supporting electrolyte and not the reactant. Re-deposition was, however, much less in Na2SO4 compared to its NaNO3/NaNO2 containing counterparts. We attribute this difference to how the local pH during electro-reduction differs in the presence and absence of NOx species. Under applied cathodic potentials, the pH at the electrocatalyst surface increases as hydroxyl ions form34,35 due to the reduction of e.g. H2O, O2, NO2− and NO3−. In particular, NOxRR results in higher currents and consequently a steeper rise in the local pH as compared to when only hydrogen reduction takes place. This process can bring the pH of a neutral electrolyte above 1234, triggering the formation of soluble Cu hydroxides. To look at the effect of electrolyte pH, we further performed experiments in 0.1 M Na2SO4 where the pH was increased to 10 by adding NaOH. As shown in Suppl. Figure 9c, it altered the amount of re-deposition observed. Lastly, to probe the influence of ammonium ions, we deliberately added NH4OH into the 0.1 M Na2SO4 carrier electrolyte, which led to rapid restructuring of the cubes as shown in Suppl. Figure 9d.
To explain these results, we consider the phase stability of Cu as a function of pH and in the presence of NH3. A complex series of reactions encompassing different acid–base chemistries, Cu(OH)2 precipitation, and complex ion formation are known for the Cu-NH3 system36,37 (see Suppl. Note 5). Specifically, the equilibrium between solid Cu(OH)2 and the Cu(NH3)42+ complex depends on the NH3 concentration, with Cu(OH)2 precipitation being favored at low concentrations due to the poor solubility of Cu(OH)236. We hypothesize that the sluggish reduction observed may be the result of transient surface Cu(OH)2 formation induced by the interfacial pH rise in the course of the NO3RR. Cu(OH)2 formation may also be more favorable in Na2SO4, due to the electrolyte’s inability to buffer the local pH increase from electro-reduction34, as compared to the KHCO3 electrolyte used in CO2RR, thereby leading to the differences in re-structuring behaviors. To validate this hypothesis, we studied the chemical changes taking place on the surface of the Cu2O cubes with operando Raman spectroscopy measurements.
Figure 5a shows the results of experiments at constant applied potentials performed with cubes electrodeposited on glassy carbon plates. At OCP, the Raman spectrum shows three features of bands centered at 415, 520, and 630 cm− 1 respectively, which are in good agreement with the reported values of Cu2O38,39. When − 0.2 VRHE was applied, the band intensity at 520 and 630 cm− 1 decreased over time but continued to persist, which is consistent with the gradual dissolution of Cu2O. At -0.4 VRHE, a new peak at 475 cm− 1 started to emerge, while the peak at 630 cm− 1 flattened, indicating oxide to metal transition. At -0.6 VRHE, the characteristic bands of Cu2O were less pronounced and additional weak bands emerged at 450, 475 and 590 cm− 1. The band around 475 cm− 1 can be assigned to the Cu-O-H vibration39–41 or Cu(OH)240, while the peak at 590 cm− 1 is often assigned to CuO38 or adsorbed O-species on Cu41. These results suggest the formation of a transient intermediate oxide or hydroxide phase during electrolysis, which is also supported by the extended presence of Cu(I) signatures and the weak but persistent Cu(II) signatures in the EC-TXM and operando XAS measurements in Fig. 3. It was, however, difficult to obtain the Raman signatures of absorbates with these samples due to the relatively low loading of the cubes that we were able to electrodeposit and the overlap of the D/G bands of the glassy carbon support with surface absorbate bands, which limits the signal-to-noise ratios at those bands. Thus, we repeated the Raman measurements with cubes electrodeposited at higher loading on carbon paper to improve the signal-to-noise ratios to identify surface adsorbed species or intermediate species in the vicinity of electrode. These results are discussed in Supplementary Note 6.
Hence, we arrive at a restructuring mechanism illustrated in Fig. 5b where surface hydroxides first form on the Cu2O cubes due to the pH increase induced by electrolysis, which delays the oxide reduction. The subsequent NH3 production and added pH rise from continued NO3RR then destabilizes this hydroxide layer and forms soluble Cu complexes, thereby initiating catalyst evolution via re-deposition from Cu complex reduction or the aggregation of migrating NPs. We also performed more EC-TEM experiments to check that the delayed restructuring kinetics extend to other pre-catalyst geometries. Suppl. Figure 10 describes the evolution during NO3RR of Cu2O truncated octahedrons and metallic frames created by pre-reducing the Cu2O cubes. Both samples are stable during early-stage reaction. With more time, re-deposition similar to that seen in the cubes was observed in the Cu2O octahedrons whereas little re-deposition was noticeable with the metallic frames, which was likely due to less dissolution occurring when we start from metallic pre-catalysts.
The correlated microscopy and spectroscopy experiments presented here therefore indicate that the morphology of the Cu catalysts during NO3RR at a given pH is governed by a complex, time- and potential-dependent interplay of three processes: (i) dissolution of oxide and hydroxide species, (ii) metal re-deposition, and (iii) oxide catalyst reduction. According to the Pourbaix diagram20,21, metallic Cu is the only stable species under the specific applied potentials and pH of our experiments, but we have shown that oxidic and metallic phases can co-exist over extended durations and over a broad range of applied potentials, which has significant implications in terms of determining the active species for producing NH3. It has been suggested previously19, based on operando spectroscopy measurements of CuO pre-catalysts, that Cu/Cu2O interfaces are responsible for NH3 production19, but these methods cannot differentiate the distribution of these species on the nanoscale. As we have shown in this work, the presence of both spectroscopic signatures in ensemble-averaging measurements does not necessarily mean that the two phases are spatially connected. Furthermore, we have demonstrated that the de-coupled metallic Cu and Cu oxide phases can persist for the extended durations at mild cathodic potentials (less than − 0.5 VRHE), and that a high residual abundance of Cu2O in the operando XAS measurements corresponded to a low NH3 production efficiency in our electrolysis data of equivalent samples. The improvement of NH3 selectivity with increasing overall metallic character of the samples, therefore, advocates that metallic Cu is the active phase for producing NH3 compared to Cu2O, in agreement with recent work on the topic15–18. In this case, the strong stability of the Cu2O cubes and their sluggish reduction kinetics in the often-used Na2SO4 carrier electrolyte, are detrimental for the NH3 production by delaying the onset of selective NH3 formation.
By showing the diverse behaviors that can be derived in different electrolytes and reaction conditions, our work also illustrates the critical need to pay attention to how the electrolyte can influence the restructuring of catalysts and the stability of oxide, hydroxide, and metallic phases before we attempt to generalize results across different studies and reactions. So far, the description of electrolyte effects in electrolysis has been largely confined to cation adsorption effects42–45 and re-structuring induced by aggressive halide anions44–46, whereas studies of pH had focused on its impact on reaction mechanisms and NH3 selectivity8,10, and not catalyst phase stability. Significant additional microscopy work, such as the one presented here will be required to separate the impact of electrolyte-driven morphological transformation from the much better understood associated electronic and chemical changes. Furthermore, current computational models still cannot rationalize the impact of an explicit complex electrolyte on the catalyst restructuring and its associated influence on the creation of active sites. Efforts to improve these models and advance the theory describing electrocatalytic processes, will undoubtedly require more accurate representations of dynamic catalyst surfaces. The challenge here is serious, since theoretical mechanistic insight must consider two simultaneously occurring dynamic processes, namely the one that the catalyst material experiences, in parallel to that undergone by the reactants, both of which being coupled and driven by the local chemical potential.47 Our results revealing phase co-existence also opens the possibility that different species may be responsible for activating specific steps of the conversion reaction. Hence, we expect operando approaches that incorporate chemically resolved microscopy within multi-modal spectroscopic investigations, as we demonstrated here, to play a vital role moving forward in the understanding of electrocatalytic processes by providing a path towards mapping such complexity.
In summary, operando EC-TEM and EC-TXM measurements have revealed that the morphologies of Cu2O pre-catalysts during NO3RR and their evolutionary pathways are sensitive to the reaction time, applied potential and the nature of the electrolyte. As expected, the rate of oxide reduction accelerates with increasing negative applied potentials, but spatially-separated oxide and metallic phases can co-exist over extended reaction times under moderately reductive potentials. More importantly, the kinetics of the different re-structuring processes, which have been unveiled here, determine that final morphology of the catalysts. Our results also indicate that the nature of the electrolyte can introduce time-dependent selectivity changes in the early stages of the catalyst restructuring, which will help resolve on-going controversies regarding the active state of Cu for selective NH3 production. Finally, this work impacts our understanding of how electrocatalysts evolve under reaction condition through the discovery of Cu oxide and hydroxide stability. In addition, we unveiled local structural and chemical heterogeneities that develop under electrochemical working conditions, even on a pre-catalyst sample initially characterized by a narrow size, shape and compositional distribution. Thus, our findings emphasize the need of operando characterization methods to establish connections between materials’ structural and compositional characteristics under specific reaction environments and external stimuli and their electrocatalytic performance.