Open Access Article
María
Romero-Angel‡
a,
Roumayssa
Amrine‡
b,
Beatriz
Ávila-Bolívar
b,
Neyvis
Almora-Barrios
a,
Carolina R.
Ganivet
a,
Natalia M.
Padial
a,
Vicente
Montiel
b,
José
Solla-Gullón
*b,
Sergio
Tatay
*a and
Carlos
Martí-Gastaldo
*a
aInstituto de Ciencia Molecular, Universidad de València, c/Catedrático José Beltrán, 2, 46980, Paterna, Spain. E-mail: sergio.tatay@uv.es; carlos.marti@uv.es
bInstitute of Electrochemistry, University of Alicante, Apdo. 99, E-03080 Alicante, Spain. E-mail: jose.solla@ua.es
First published on 25th March 2024
The combination of compositional versatility and topological diversity for the integration of electroactive species into high-porosity molecular architectures is perhaps one of the main appeals of metal–organic frameworks (MOFs) in the field of electrocatalysis. This premise has attracted much interest in recent years, and the results generated have also revealed one of the main limitations of molecular materials in this context: low stability under electrocatalytic conditions. Using zirconium MOFs as a starting point, in this work, we use this stability as a variable to discriminate between the most suitable electrocatalytic reaction and specific topologies within this family. Our results revealed that the PCN-224 family is particularly suitable for the electroreduction of molecular nitrogen for the formation of ammonia with faradaic efficiencies above 30% in the presence of Ni2+ sites, an activity that improves most of the catalysts described. We also introduce the fluorination of porphyrin at the meso position as a good alternative to improve both the activity and stability of this material under electrocatalytic conditions.
Provided the samples have sufficient chemical stability, their use in photo and electrocatalytic processes is of great importance to demonstrate the potential of reticular design in catalysis.7 MOFs have been used as catalysts with promising results in reactions as the hydrogen evolution reaction (HER),8 oxygen reduction/evolution reaction (ORR/OER),9 CO2 reduction reaction (CO2RR),10 and nitrogen reduction reaction (NRR).11 Porphyrin molecules are particularly appealing in this context. Besides their intrinsic electrocatalytic activity in solution,12–14 their rigid and robust backbones allow for the assembly of reticular frameworks with fine control over their intermolecular distance, their density, the chemical environment, and greater/lower accessibility to substrates depending on the net porosity of the material. This ease to direct framework assembly is supported by the huge number of porphyrin-MOFs reported with different metal nodes and secondary building units (SBUs).15–17 Among them, Zr/carboxylate-based SBUs, in particular, those based on the Zr6O4(OH)4 cluster (Zr6), seem to be one of the most synthesized, at least in terms of electrocatalytic efficiency, with at least one member of this subfamily standing out in the HER, OER, CO2RR or NRR. This is arguably due to the superior thermal and water/moisture stability of Zr-MOFs,18,19 that are expected to prevent their degradation under the conditions imposed by electrocatalytic reactions. The Zr6 SBU also combines a high versatility in its connectivity index, from 12 to 6 connection points, with a high persistence in dissolution, thus facilitating the synthesis of different topologies with variable porosities for a large pool of potential catalysts (Fig. 1).
![]() | ||
| Fig. 1 Some reported topologies for common Zr-oxo clusters and TCPP organic linkers. Representative MOF examples are given for each corresponding net connectivity and topology. | ||
Most porphyrin-based molecular frameworks have been built from tetratopic meso-tetra(carboxyphenyl)porphyrin (TCPP) linkers, or some other related porphyrins with reduced/increased topicity and/or extended backbones.20,21 Here we use several representative TCPP–MOFs to systematically study the electrocatalytic activity of such frameworks for ammonia generation as a function of their chemical stability, metal substitution (MTCPP), or the functionalization of the porphyrin ring in the metal position with fluorine atoms (TCPP–F) for tailorable pore environments.
In the search for a synthetically robust platform compatible with the iterative variation of metals and linkers, we chose to test in our laboratory some of the most representative members of the Zr–TCPP MOF family (Fig. 1) including: MOF-525,24 PCN-223,25 MOF-545
24 (also referred to as PCN-22226 and MMPF-6
27), PCN-224,28 and PCN-226.29 The reader should note that in this list we included not only the classical Zr6 clusters, but also the closely related Zr8 node and the rod-shape Zr-oxo SBU in PCN-226 for the shake of completeness. We attempted the synthesis of all these MOFs based on the original reported procedures,24–26,29 except for the case of PCN-224, for which we chose to follow the revisited synthetic protocol recently reported by Lotsch's group.22 According to our powder X-ray diffraction (PXRD) data (Fig. 2), we straightforwardly obtained phase pure crystalline phases of MOF-525, PCN-224 and PCN-226, while our attempts for the synthesis of MOF-545 and PCN-223 were unsuccessful and yielded amorphous phases. These results highlight the difficulty in controlling the synthesis of Zr–TCPP MOF, where different members of the family share very similar or even equal compositions. Our results are in line with recent studies that demonstrate the acute influence of water content on controlling phase formation in Zr–porphyrin-based MOFs30 and the problems for their synthesis reproducibility due to their sensitivity to slight changes in the experimental conditions.23,31,32
As a result of our initial screening, we selected MOF-525, PCN-224, whose crystal structure is represented in Fig. 2c, and PCN-226 as potential platforms to implement our meta-substituted TCPP roadmap.
Our chemical stability tests in conventional aqueous electrolytic solutions suggest that MOF-525, PCN-226 and PCN-224 do not display sufficient chemical stability to be used for the electrochemical CO2RR but they show comparatively higher stability in the medium used for the NRR. Among them, PCN-224 showed the highest chemical stability, displaying a minimum percentage of metal leaching after 24 hours (0.04%). Consequently, PCN-224 is the ideal candidate to investigate the electrochemical activity of the pristine framework or other derivatives belonging to the same family towards the production of ammonia from nitrogen.
:
4 M/Zr ratio consistent with the corresponding unit formula. Crystal size and morphology were studied by using scanning electron microscopy (SEM). Fig. 3b shows little variation in morphology across the series. All solids display cube-like crystals; however, two different size ranges are observed depending on the sample. While PCN-224–Cu shows crystals of about 1 μm very similar to those of PCN-224, the Co and Ni analogues are comparatively bigger. Leaving aside this small divergence in size, the properties of the PCN-224–M series agree well with the characteristics reported for this MOF and rule out the formation of contaminant oxide phases. Permanent porosity was also analysed with N2 isotherms at 77 K (Fig. S12†).28 All samples exhibit a reversible type-I isotherm. While PCN-224 shows the maximum N2 uptake with a value consistent with that originally reported,15 the metalated analogues display comparatively lower uptakes as previously reported by others.42,43
Next, the electrocatalytic activity of the materials toward the NRR was systematically investigated. For this purpose, electrodes of the different samples were prepared, by mixing the corresponding solids with the binder agent, and sprayed over a Toray carbon paper using an air brushing technique for a final electrode loading of about 0.7 mg cm−2. A platinum wire and a leakless AgCl/Ag (saturated KCl) were used as the counter and reference electrodes, respectively. All potentials were referred to the reversible hydrogen electrode (RHE). Cyclic voltammetry (CV) curves were recorded at 50 mV s−1 in Ar and N2 saturated 0.1 M Na2SO4 solution. Chronoamperometry (CA) measurements (electrolysis) were subsequently carried out at −0.04 V, −0.14 V and −0.34 V vs. RHE for 2 h. The experimental details are included in ESI Sections S2 and S7†. It is worth highlighting that these experiments were performed in the same Na2SO4 aqueous electrolyte for which we confirmed the stability of PCN-224 under open circuit potential conditions with negligible metal leaching after 24 hours. Faradaic efficiencies (FEs) and ammonia yields were determined using the amount of ammonia produced as determined by the indophenol blue method.42 Before each electrolysis, five voltammetry cycles between 0.56 V and −0.94 V vs. RHE were recorded to verify the status of the electrodes. Fig. 3c displays the voltammetric profiles corresponding to PCN-224–Co, PCN-224–Ni and PCN-224–Cu in the potential range from 0.2 V to −0.6 V vs. RHE. The voltammetric responses of the electrodes in Ar (gray line) and N2 (color line) are very similar, with a slight increase of the negative currents in the N2 saturated solution electrolyte which could be related to an enhanced activity for the electrochemical reduction of N2 to ammonia. To better evaluate the electrocatalytic properties of the samples towards N2 electroreduction, CAs were performed at different potentials of interest. Fig. 3d shows the variation of current density versus time at a constant potential of −0.04 V (the potential at which ammonia production reached its maximum) for all catalysts. The results demonstrate that the behaviour of the samples is strongly determined by the type of metal incorporated. In agreement with the voltametric curves (Fig. 3c), Ni and Cu samples exhibited lower current density compared to Co. In addition, all samples showed a relatively stable current during the electrolysis time. Regarding the FE (Table S3†), the results indicate that the FE towards ammonia ranges from 4 to 32% depending on the nature of the metal. Interestingly, in all cases, the highest FE value is achieved at −0.04 V; the Co, Ni, and Cu samples show FE values of about 11.4, 32.3 and 21.7%, respectively. In all cases, the FEs systematically decrease at more negative potentials due to the increased effect of the HER. Thus, at −0.34 V vs. RHE, the FE values are the lowest independent of the nature of the electrode. The FE values discussed previously clearly indicate that at a relatively positive potential value of −0.04 V vs. RHE, the incorporation of a metal into the PCN-224 framework significantly enhances the electrocatalytic activity of the electrocatalyst, with the Ni metalated being the most efficient one with an FE of about 32%. To rule out the possible contribution of the Toray paper, control experiments were performed with the bare Toray paper electrode. The results are shown in Fig. S27 and S28.† The results indicate the absence of detectable ammonia during these control experiments. To better understand our finding, Fig. 4 compares this 32% FE value with the FEs reported with other electrocatalysts, measured in a comparable range of potentials of ± 0.2 V and obtained in alkaline electrolyte which is expected to suppress the hydrogen evolution reaction for improved NRR activities (Table S4†). The FE of 32.29% registered for PCN-224–Ni is quite higher than the 16.2% reported for the related 2D microporous MOF Fe–TCPP, built from Zn(II) paddlewheel clusters and the metalated TCPP linker.44,45 Other electrocatalysts prepared by the calcination/pyrolysis of MOFs to generate the corresponding Bi nanoparticles,46 CeO2 nanorods47 Y-stabilized ZrO2,48 Co3O4 nanopolyhedra,49 iron doped nanocarbons,50 or NiO/Ni microtubes51 all supported on carbon also display comparatively lower activities below a FE of 13%. Higher activities have been reported for other MOF composites reliant on the use of MXenes as ZIF-67@Ti3C2,51 but they are comparatively lower with efficiencies near 20%. The activity associated with Zr6 clusters in defective UiO-66 frameworks has been also shown to be very high for this type of reaction with efficiencies up to 48%.52 This comparison clearly points out the relevance of our results as it implicates the MOF exclusively as a crystalline support of porphyrin complexes acting as electrocatalytic centers.51 Despite the excellent results in terms of FE, it is important to mention that the electrodes cannot be reused, and the FE values decrease significantly during their second use. This suggests a possible degradation of the electrodes under electrochemical working conditions, which is not obvious for the Cu and Ni derivatives (Fig. S26a†), but is more evident for the cobalt derivative (Fig. S26b†). In the last case, we observed a strong colouration of the solution that agrees well with the characteristic MTCPP UV-vis signals, which was concomitant with the loss of performance. In addition, ICP analysis of the isolated electrolyte solution after electrocatalytic experiments with PCN-224–Ni confirmed an increase in the nickel concentration in the electrolyte buffer of 29.2 ppb compared to the 0.1 ppb present before the reaction, which might be indicative of framework degradation. These results highlight the need to improve the stability of the samples under electrochemical working conditions.
![]() | ||
| Fig. 4 Overview of the FE of the PCN-224–M family compared to that of other NRR electrocatalysts measured in alkaline electrolytes and under comparable voltage conditions. See ESI Table S1† for more details. | ||
To evaluate the impact of the electron withdrawing fluorine groups on the electrocatalytic NRR we prepared the corresponding electrodes and tested the system under the same conditions used for the PCN-224–M family. The CV of the sample in Ar and N2-saturated solution was essentially similar to those shown by PCN-224–Ni in Fig. 3c, although for the PCN-224–Ni(F) sample (Fig. 5c), the voltammetric responses revealed an increase in current in the N2-saturated electrolyte which could be related to an enhanced activity for the electrochemical reduction of N2 to ammonia. With respect to the FE values (Table S4†), the incorporation of F also improved the FE at −0.04 V, from 32.3% to about 34.5%. As for the unfluorinated samples, for more negative potentials than −0.04 V, the FEs systematically decrease due to the increased effect of the HER.37 Regarding the ammonia yield rates, a maximum value of 0.9 μgNH3 h−1 mgcat−1 was obtained at −0.04 V vs. RHE. Remarkably, compared to PCN-224–Ni, PCN-224–Ni(F) showed constant FE values up to the third electrode use. This implies an increased stability of the samples under electrochemical working conditions. Hydrophobic fluorination has been used before to enhance MOF stability,53–55 and in the particular case of the NRR, to increase local concentration of N2 around the catalytic centre, thus limiting proton transportation and favouring the NRR over the HER.56–58 To confirm the improved stability associated with the incorporation of fluorinated substituents into the porphyrin ring, we analysed the electrolyte solutions after the first electroreduction reaction by ICP. Compared to the near 30 ppb leached by PCN-224–Ni, its fluorinated equivalent is only 4.3 ppb confirming a significant improvement in the stability of this MOF under electrocatalytic conditions. We also incubated both MOFs in the electrolyte used in the electrocatalytic experiment for 24 h. Fig. S30† shows that while PCN-224–Ni amorphizes considerably under these conditions, the crystallinity of the fluorinated derivative is retained without apparent amorphization. The water isotherms of PCN-224–Ni and PCN-224–Ni(F) at 298 K (Fig. S29†) provide additional evidence of the enhanced stability towards hydrolysis of the fluorinated framework, which is less prone to collapse and displays a much higher water uptake than PCN-224–Ni despite the increase in hydrophobicity.
Although all these experiments support the enhanced stability and reusability of the fluorinated framework, they do not provide a clear answer for the origin of the experimentally observed enhancement of FE. To clarify which was the most likely cause for the observed change in performance, we next carried out DFT calculations.
Based on these preliminary results we analyzed separately the two generally accepted mechanisms for the NRR. First, we considered an associative Heyrovsky pathway60 that relies on the adsorption of N2 molecules to the active site followed by hydrogenation, from the combination of protons in solution and electrons from the electrode. However, our simulations reveal that this mechanism is quite unlikely in our case. The binding of N2 with the MTCPP–X units is very weak or metastable with N2–M distances near 3.3 Å in all cases, compared to the 2.0 (Co), 2.5 (Cu), and 2.7 Å (Ni) that would be required for nitrogen adsorption and activation (Fig. S32†).
We next considered the Tafel–Volmer mechanism,61 which involves hydrogen dissociation (1), electron transfer to the catalyst (2), the adsorption of solvated protons from solution (3), and the adsorption of N2 for its reduction (4). We calculated the adsorption energy for all systems by combining steps (2)–(4) correlating the trends in catalytic activity with the transition metal and fluorination of the porphyrin ring (Fig. 6a). As anticipated by our experimental results, the N2 molecule was preferentially adsorbed in NiTCPP–H and NiTCPP–F, compared to the Co and Cu analogues. According to our calculations, the N2 molecule not only interacts preferably with the Ni site, but also forms a hydrogen bonding with the proton adsorbed by the porphyrin for enhancing adsorption and reducing the energy barrier in the rate determining step (Fig. 6b). This effect is even more acute for a fluorinated framework, for which the electron transfer and nitrogen adsorption are comparatively more favourable making it the best electrocatalyst of the family (Table S6†). The addition of the proton to the porphyrin core in CuTCPP–H is stabilized by 0.01 eV, although it elongates the Cu–N interactions with the pyrrole nitrogens by 0.5 Å, suggesting that the N2 is weakly bound as suggested before by the Heyrovsky mechanism. Finally, the electrons in CoTCPP–H are located at the lowest chemical potentials (2.65 eV) across the series (Table S6†). This suggests a poorer ability for accepting electrons in the first step of the NRR for the lowest electrocatalytic activity of the family.
Footnotes |
| † Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3ta07004b |
| ‡ These authors contributed to the work equally. |
| This journal is © The Royal Society of Chemistry 2024 |