Nilmadhab
Mukherjee
a,
Ashadul
Adalder
a,
Narad
Barman
b,
Ranjit
Thapa
b,
Rajashri
Urkude
c,
Biplab
Ghosh
c and
Uttam Kumar
Ghorai
*a
aDepartment of Industrial Chemistry & Applied Chemistry, Swami Vivekananda Research Centre, Ramakrishna Mission Vidyamandira, Belur Math, Howrah - 711202, India. E-mail: uttam.indchem@vidyamandira.ac.in
bDepartment of Physics and Centre for Computational and Integrative Sciences, SRM University – AP, Amaravati 522240, Andhra Pradesh, India
cBeamline Development & Application Section, Bhabha Atomic Research Center, Trombay, Mumbai 400 085, India
First published on 2nd December 2023
The electrochemical reduction of nitrate to ammonia (NO3RR) catalyzed by metal organic frameworks (MOFs) is a promising and efficient method for reducing nitrate pollution in water while simultaneously producing a valuable product, ammonia. Herein, we report the 3D nanoarray architecture of the metal organic complex Fe-(tetracyanoquinodimethane)2 Fe(TCNQ)2 as an efficient electrocatalyst that exhibits a high ammonia yield rate of 11351.6 μg h−1 cm−2 and faradaic efficiency (FE) of 85.2% at −1.1 V vs. RHE and excellent catalytic stability up to 2 days. The excellent catalytic performance is evaluated by ATR-FTIR spectroscopy and a series of control experiments. Density functional-based theoretical calculations are carried out to identify Fe–N4 active sites in metal–organic network structures. This study showcases the advancement of transition metal-based organic frameworks as very effective electrocatalysts for the reduction of nitrate to ammonia (NH3).
In this work, a Fe(TCNQ)2 nanorod array has been fabricated on a copper foam (CF) substrate, which is a binder-free electrocatalyst-electrode assembly for NO3RR studies. FESEM observations confirm the distinct 3-dimensional (3D) nanoarray architecture. Extended X-ray absorption fine structure (EXAFS) spectra identify the distinctive Fe–N4 coordination environment which is the active center for the NO3RR. Hence, the 3D Fe(TCNQ)2/CF system delivered an ammonia yield rate 11351.6 μg h−1 cm−2 and FE 85.2% at −1.1 V vs. RHE along with electrocatalytic stability of more than 2 days. Attenuated total reflection–Fourier transform infrared (ATR-FTIR) spectroscopy confirms the various intermediate species, specifically *NO, *NH, and *NH2 formed during NH3 formation. Furthermore, density functional based theoretical calculations demonstrated the mechanism of the eight electron transfer associated NO3RR. Additionally, several control assays and isotopic labelling experiments confirm the source of nitrogen in the produced NH3, as well.
Fig. 1 (a) Growth mechanism of Fe(TCNQ)2/CF, (b, c) FESEM images of CuTCNQ/CF, (d) FESEM image of Fe(TCNQ)2/CF, and (e) TEM image of Fe(TCNQ)2/CF. |
Prior to this, a nanorod array of Cu(TCNQ) had been successfully synthesized on a copper foam substrate defined as Cu(TCNQ)/CF. The transformation of Cu(TCNQ)/CF into Fe(TCNQ)2/CF involved a controlled cationic exchange process, where Fe2+ ions effectively exchange Cu+ ions.42 This innovative Fe(TCNQ)2/CF nanorod array configuration served as a conductive electrode in the realm of electrochemical nitrate reduction, yielding notably enhanced performance in the context of the nitrate reduction reaction. The influence of the morphology of the resultant metal-TCNQ composites was systematically examined through scanning electron microscopy (SEM) and transmission electron microscopy (TEM). The FESEM depiction of the Cu(TCNQ) specimen unveils developed growth of Cu(TCNQ) nanorods achieved through the surface coverage of Cu foam, as demonstrated in Fig. 1b and c. Upon investigating the Fe(TCNQ)2 sample subsequent to the cationic exchange process, its distinct nanoarray architecture is observed to persist, as evident in Fig. 1d. And Fig. 1e shows a TEM image of Fe(TCNQ)2 where the nanorods are ∼150–160 nm in width.
To gain insight into the structural properties of the catalyst, we conducted X-ray diffraction (XRD) analysis. Fig. 2a shows the XRD patterns of both Cu(TCNQ)/CF and Fe(TCNQ)2/CF. The distinct peaks in the XRD plot correspond to the X-ray diffraction patterns of Cu(TCNQ) and Fe(TCNQ)2, thereby revealing important information about their crystal structures which is well matched by the previous data.42,43 The X-ray photoelectron spectroscopy (XPS) spectrum of Fe (TCNQ)2/CF is depicted in Fig. S1† and shows that the catalyst is composed of the elements Fe, C, N and Cu. The oxygen (O) element may appear due to material contamination or oxidation on the product surface.42 At about 398.6 eV, an XPS N 1s spectrum peak (Fig. 2b) is observed due to Fe–NC; also a shake-up peak comes at 400.1 eV in the Fe(TCNQ)2/CF system. The peaks with binding energy at 284.58 and 285.94 eV are produced by carbon bonds (C–H) and nitrile carbon bonds (C–N) in the TCNQ system, respectively (as shown in Fig. 2c). In addition, the unsaturated carbon bond produces a satellite at 288.1 eV. The two main binding energies for Fe 2p3/2 and Fe 2p1/2 are shown in Fig. 2d at 712.7 and 724.1 eV, respectively.41,42,44 To elucidate the distinctive Fe–N4 coordination and to unravel the intricate structural attributes of the catalytically active sites within the developed Fe(TCNQ)2/CF catalyst, a comprehensive investigation employing X-ray absorption spectroscopy (XAS), encompassing both X-ray absorption near edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) spectra, was carried out. Fig. 2e illustrates the XANES absorption edge of Fe in the Fe(TCNQ)2/CF catalyst, positioned intermediate to the reference standards of elemental Fe foil and Fe2O3. This positioning of the absorption edge signifies that the average valence state of Fe within the Fe(TCNQ)2/CF catalyst falls within the range of 0 to +3. The detailed structural insight derived from Fig. 2f showcases the k3-weighted Fourier transform (FT) of the Fe K-edge EXAFS spectra for Fe(TCNQ)2/CF, Fe foil, and Fe2O3. The distinct peak observed near 1.42 Å in the Fe(TCNQ)2/CF sample corresponds to Fe–N coordination, a significant indicator of the active sites for the NO3RR process. Contrastingly, the negligible presence of peaks at 2.19 Å and 2.55 Å, which signify Fe–Fe metal pair interaction, establishes the uniqueness of the Fe(TCNQ)2/CF catalyst.45,46 This analysis collectively points towards a practical coordination number of around 4 for Fe–N4 in the Fe(TCNQ)2/CF catalyst (Table S1†), strongly implicating its role as the dominant active site for the NO3RR.45,46 Additionally, a wavelet transform (WT) analysis was undertaken on the K3-weighted EXAFS signal to offer further clarity on atom backscattering, as depicted in Fig. 2g–i. In this context, the distinctive Fe foil exhibited a WT maximum at 7.3 Å−1, aligned with Fe–Fe interactions. Conversely, the reference Fe2O3 displayed two prominent WT maxima at 7.4 and 3.2 Å−1, corresponding to Fe–Fe and Fe–O interactions, respectively. Also, the Fe(TCNQ)2/CF catalyst showcased a singular WT maximum situated at 2.9 Å−1, unequivocally attributed to Fe–N4 interactions, thereby corroborating its structural significance in the catalytic process.45,46
Evaluation of the catalytic efficiency of metal-TCNQ/CF electrodes regarding nitrate reduction was performed using a conventional three-electrode system (H-type cell), as detailed in the Experimental methods section. Validation of the electrocatalytic performance for NO3− reduction was established by observing the amplified current density within the linear sweep voltammogram (LSV) profile following the introduction of 0.2 M NaNO3 in 0.1 M Na2SO4 (illustrated in Fig. 3a). Furthermore, electrochemically active surface area (ECSA), obtained via double-layer capacitance (Cdl) assessment (Fig. 3b and S20), revealed significant insights. In particular, the Fe(TCNQ)2/CF configuration exhibited an Cdl and corresponding ECSA of 11.2 mF cm−2 and 252 cm2, which exceeded both CuTCNQ/CF (4.2 mF cm−2 and 94.5 cm2) and CF (0.43 mF cm−2 and 9.7 cm2), thereby confirming its higher internal factoring intrinsic reduction rate. To comprehensively assess the catalyst's catalytic proficiency, electrochemical impedance spectroscopy (EIS) plots were harnessed to elucidate its charge transfer (CT) kinetics. The parameters extracted from these measurements were effectively modeled using a Randles circuit, incorporating solution resistance (Rs) and charge transfer resistance (Rct), yielding an accurate representation of the Nyquist plots' complexities (Fig. 3c). Impressively, Fe(TCNQ)2/CF exhibited notably diminished Rct values in stark comparison to CF and CuTCNQ, indicative of the heightened conductivity and enhanced proton transfer facilitation intrinsic to Fe(TCNQ)2/CF, ultimately underscoring its superior nitrate reduction efficacy. For a comprehensive exploration of the ammonia production efficiency within the catalytic system, chronoamperometry (CA) experiments were performed over 1 h reduction time, as exemplified in Fig. S13.† After 1 h of testing, all electrolyte solutions were used for indophenol blue method for ammonia detection by UV-vis spectroscopy (details in the ESI†). The UV-vis spectra of the electrolyte solutions (after electrolysis) exhibited increasing peak intensities at different applied potentials (Fig. 3d–f). Fe(TCNQ)2/CF manifested an exceptional electrocatalytic prowess in catalyzing the NO3RR, achieving an impressive maximum faradaic efficiency (FE) of 85.2% at −1.1 V vs. RHE (Fig. 3g), accompanied by a peak ammonia yield rate of up to 11351.6 μg h−1 cm−2 at −1.1 V vs. RHE (Fig. 3h). Similarly, both the FE and ammonia yield rate of CuTCNQ/CF (84.8% & 10562 μg h−1 cm−2) and CF (67.3% & 5661 μg h−1 cm−2) outperformed expectations for the NO3RR.
These collective findings serve to underline the remarkable potential of Fe(TCNQ)2/CF as a superior electrocatalytic agent for nitrate reduction, signifying significant strides in the field of electrochemical nitrogen fixation. To quantitatively gauge catalytic activity, the turnover frequency (TOF, s−1) was meticulously calculated based on ammonia production rates at various potentials. The distinctive histograms of TOF plotted against potential in Fig. 3i vividly depict the ammonia generation rates at Fe sites within the Fe (TCNQ)2/CF systems. Specifically, at a potential of −1.1 V, Fe (TCNQ)2/CF exhibited a noteworthy peak TOF of 2.6 s−1. A series of control experiments were systematically conducted to unequivocally confirm the generation of ammonia during the nitrate reduction reaction (NO3RR) (Fig. 4a). Initially, no signal was observed for nitrate free electrolyte solution also and a moderate signal was discerned in a solution composed of 0.2 M NaNO3 and 0.1 M Na2SO4 in the absence of metal-TCNQ catalyst growth copper foam (CF). Subsequently, the CuTCNQ/CF system produced ammonia during the NO3RR. In parallel, the Fe (TCNQ)2/CF system evaluation within a solution composed of 0.1 M Na2SO4 and 0.2 M NaNO3 unambiguously established the occurrence of ammonia production during the NO3RR. To unequivocally establish the provenance of ammonia from NO3− ions, we conducted 1H–nuclear magnetic resonance (NMR) analysis employing a 14N and 15N-labeled NO3− electrolyte solution, showcased in Fig. 4b.16 The resultant 1H-NMR spectra depicted a distinctive doublet peak corresponding to 15NH4+, providing unambiguous confirmation of ammonia's genesis via the electrochemical reduction of nitrate ions (15N-labeled), eliminating any confounding influences of extraneous factors. Equally noteworthy, in experiments using a 14NO3− containing electrolyte, a distinct triplet peak corresponding to 14NH4+ further reaffirmed the specific nitrate-driven ammonia formation. Crucially, in scenarios where nitrate was absent within the electrolyte during the electroreduction process, the absence of detectable ammonia served as a definitive confirmation of the dependency on nitrate for ammonia production. Subsequent to the electrolytic process, we conducted attenuated total reflection-Fourier transform infrared (ATR-FTIR) spectroscopy (Fig. 4c). Evident within the FTIR spectrum was a deformation vibration signal at 1637 cm−1, which is characteristic of H2O. Remarkably, discernible peaks at 1795 cm−1, 1395 cm−1, and 1101 cm−1 were attributed to various intermediate species, specifically *NO, *NH, and *NH2. The primary formation of NH3 was conclusively ascertained through the prominent peak at 1530 cm−1.47–49 Intriguingly, spectral analysis revealed a direct correlation between applied potentials and ammonia intensity, substantiating the influence of potential on product yield. As relevant to catalyst evaluation, stability analysis is paramount. Obviously, a prolonged electrolysis test at −1.1 V was performed to measure the long-term stability of Fe(TCNQ)2/CF (Fig. 4d). Impressively, as depicted in Fig. 4d, a marginal decrease in current density was observed only after 50 h of continuous operation, which underscores the exceptional stability of Fe(TCNQ)2/CF under practical working conditions. Addressing reusability, an important aspect for practical application, the nitrate reduction protocol was subjected to ten consecutive cycles, each lasting one hour (illustrated in Fig. 4e). Also, the catalyst consistently maintained high faradaic efficiency (FE) and yield rates throughout these cycles, firmly establishing its recyclability for the NO3RR. The XRD pattern (Fig. S21†) and FESEM image (Fig. S22†) observed for post-electrolysis Fe(TCNQ)2/CF specimens (as depicted in respective figures) provide further compelling evidence of the extraordinary durability of our catalyst under operational conditions. This comprehensive analysis substantiates the absence of deleterious species and further underscores the catalyst's impeccable stability.
Theoretical studies of the NO3RR, on the catalyst Fe-TCNQ, have been carried out by introducing an optimized model structure of Fe-TCNQ (Fig. S23†), where the lattice parameters of the system are 6.88 Å, 11.27 Å and 20 Å along a, b and c axes respectively. Along the c-axis, a length of 20 Å is incorporated to nullify the false interactions due to the virtual image. The theoretical study comprises three sections; the first one is about checking for improved NO3 adsorption on the metal site of the catalyst over hydrogen adsorption and in the next part we have illustrated the reaction mechanism and corresponding free energy diagram with a possible pathway for the NO3RR. In the last section, we have discussed the catalytic activity of the given catalyst.
Since the HER is taking part as a major competitor, hampering the reaction phenomena parallelly by blocking active sites, it plays a crucial role while forming NH3 from NO3 reduction. We have found that the adsorption process of NO3, on the Fe-TCNQ catalyst, is more exothermic in nature than the same for hydrogen with free energy changes of −1.04 eV and 0.13 eV respectively as shown in Fig. 5a, indicating that reaction selectivity moves towards the NO3RR.50,51 Now, when NO3 comes and sits on the active metal site, it gets reduced by taking charges from valence orbitals of the metal atom. We have verified this reduction phenomenon quantitatively and qualitatively by using Bader charge analysis and a known parameter, called the d-band centre35 respectively. The charges associated with the metal atom, before and after adsorption of NO3−, are estimated to be 6.82e and 6.686e respectively. The shifting of the d-band centre from 2.73 eV to 2.35 eV, after adsorption of NO3− towards the Fermi level, points out that charge transfer occurred. To visualise the difference in charge density i.e., loss and gain of charges by respective atoms after NO3− adsorption, we have presented Fig. 5b and shifting of the d-band centre is depicted in Fig. 5c and d.
To understand the reaction mechanism of the NO3RR on the given catalyst, we have considered three known possible pathways (Fig. 6a), reported in the literature52 for metal surfaces and single atom catalysts etc., where protonation takes place at the 2nd step of the possible reaction pathway i.e., NO3− → *NO3 → NO3H and all the free energy calculations of the respective intermediates are improved by adding the correction term (EZPE-TS). In Fig. 6b–d, we have demonstrated the reaction mechanism, along the most favourable pathway i.e., Path I (chosen by the lowest ΔG of the potential determining step), and corresponding full free energy profile for the NO3RR respectively. From the free energy diagram, it is observed that NO3− comes and easily gets adsorbed on the active metal site and forms the intermediate *NO, by crossing a series of exothermic steps with a free energy change of −1.04 eV (step 1), −0.26 eV (step 2), −1.62 eV (step 3), −0.35 eV (step 4) and −1.63 eV (step 5) respectively. During this five step process, two H2O molecules are liberated into the electrolyte in the following way: *NO3H + 3H+ + 3e− → *NO + 2H2O. In the next step, protonation takes place in an endothermic way by the (H+ + e−) pair and forms the intermediate *NOH, which is recognized as the potential determining step (PDS),53 with a free energy change of 0.30 eV. After this, the (H+ + e−) pair attacks serially and causes the generation of a series of intermediates that follow Path I (*NOH → *NHOH → *NH → *NH2 → *NH3), where the first step is endothermic and the rest of the intermediates are formed spontaneously with a free energy change of 0.22 eV, −0.63 eV, −1.53 eV, and −0.79 eV respectively. At last, desorption of NH3 also, from the surface of the catalyst, takes place spontaneously. Now, estimation of the highest energy barrier known as the PDS for two pathways, along Path I of *NOH → *NHOH and Path III of *NOH → *N, tells us that Path I is more favourable over Path III.
On the other side, the thermodynamic step *NHOH → *NH is more spontaneous than the step *NHOH → *NH2OH; hence there would be very little chance for Path II to occur. Therefore, based on the above discussion we can conclude that reduction of NO3− to NH3 on the Fe-TCNQ catalyst takes place smoothly through Path I over Path II and Path III. To find the catalytic activity of the proposed catalyst, there are so many factors to be investigated, e.g., limiting potential or the PDS, HER selectivity, change in free energy during desorption of NH3 from the catalyst surface etc. In this DFT study, we have obtained a lower value of limiting potential 0.30 V (Fig. 6d) and the given catalyst is observed to suppress the side reaction, HER phenomena that we already mentioned in the above section, resulting in higher selectivity towards the NO3RR. Besides, the change in free energy in the case of the desorption process of NH3 from the surface of the host catalyst is negative, −0.09 eV, which signifies that Fe-TCNQ is an efficient catalyst making the reduction process more feasible.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3ta05300h |
This journal is © The Royal Society of Chemistry 2024 |