Open Access Article
This Open Access Article is licensed under a
Creative Commons Attribution 3.0 Unported Licence

Engineering hydrophobic–aerophilic interfaces to boost N2 diffusion and reduction through functionalization of fluorine in second coordination spheres

Sakshi Bhardwaj a, Sabuj Kanti Das a, Ashmita Biswas a, Samadhan Kapse b, Ranjit Thapa b and Ramendra Sundar Dey *a
aInstitute of Nano Science and Technology (INST), Sector-81, Mohali 140306, Punjab, India. E-mail: rsdey@inst.ac.in
bDepartment of Physics, SRM University, Andhra Pradesh 522240, India

Received 13th June 2023 , Accepted 31st July 2023

First published on 1st August 2023


Abstract

Ammonia is a crucial biochemical raw material for nitrogen containing fertilizers and a hydrogen energy carrier obtained from renewable energy sources. Electrocatalytic ammonia synthesis is a renewable and less-energy intensive way as compared to the conventional Haber–Bosch process. The electrochemical nitrogen reduction reaction (eNRR) is sluggish, primarily due to the deceleration by slow N2 diffusion, giving rise to competitive hydrogen evolution reaction (HER). Herein, we have engineered a catalyst to have hydrophobic and aerophilic nature via fluorinated copper phthalocyanine (F-CuPc) grafted with graphene to form a hybrid electrocatalyst, F-CuPc-G. The chemically functionalized fluorine moieties are present in the second coordination sphere, where it forms a three-phase interface. The hydrophobic layer of the catalyst fosters the diffusion of N2 molecules and the aerophilic characteristic helps N2 adsorption, which can effectively suppress the HER. The active metal center is present in the primary sphere available for the NRR with a viable amount of H+ to achieve a substantially high faradaic efficiency (FE) of 49.3% at −0.3 V vs. RHE. DFT calculations were performed to find out the rate determining step and to explore the full energy pathway. A DFT study indicates that the NRR process follows an alternating pathway, which was further supported by an in situ FTIR study by isolating the intermediates. This work provides insights into designing a catalyst with hydrophobic moieties in the second coordination sphere together with the aerophilic nature of the catalyst that helps to improve the overall FE of the NRR by eliminating the HER.


Introduction

Ammonia has been realized as an essential commodity in the current world scenario where there is an utmost urge towards sustainability.1 At the moment the research community is engrossed in scaling up the electrochemical method of ammonia synthesis to get into the competition of large-scale ammonia production with the energy-intensive Haber–Bosch process.2 To date, a lot of progress has been witnessed in the eNRR in terms of catalyst development, method optimization and identification of false positive responses. However, the catalyst development genre is under extensive investigation and still has enough room for progress. Mainly, an electrocatalyst for the NRR is developed keeping in mind a few bottlenecks of the NRR like (a) selectivity towards N2 adsorption, (b) competitive hydrogen evolution reaction (HER),3 (c) catalysts having Lewis acid sites that influence N2 bond polarization, etc.4 Among them, the HER is looked down upon as it causes acute hindrance for a smooth NRR on the catalyst surface making industrial-scale ammonia synthesis a distant dream. This is because the NRR is a sluggish process involving a multiphase reaction pathway with six electron–proton transfer processes.5 At the same time, the HER needs only two protons and two electrons and is relatively fast.6 Moreover, in the aqueous medium, materials favour H adsorption due to unrushed N2 molecules which occupy most of the active sites that lead to parasitic hydrogen evolution by using up the majority of available electrons and reduces the selectivity for the NRR.7

Therefore, HER suppression will resolve the bottleneck problem of the NRR and will alleviate its selectivity and conversion efficiency for its practical application. Different strategies have been adopted to suppress the HER: (1) by regulating the accessibility of electrons and protons;8 (2) by modulating the chemical equilibrium of the HER; (3) by designing suitable catalysts. Out of these strategies, although the catalyst designing is an emerging field to mitigate the adsorption of H atoms, it has been less explored. In this respect, early transition metals (Ti, Zr, Y, and Sc)9,10 are in use due to their intrinsic HER restriction. Recently, atomically dispersed catalysts with late transition metals are widely explored for the NRR due to their unsaturated co-ordination configuration and maximum utilization of atoms.3,11 However, this metal specificity limits exploring wide-ranging materials for the NRR. As the NRR is highly sensitive to microenvironments, it is urgently needed to construct catalysts with rationally engineered active sites.12 So, it could be implemented for any metal/non-metal irrespective of its inherent properties and induce its reactivity towards the NRR, bypassing the otherwise feasible HER process. In this realm, for HER suppression and maximizing the mass transfer efficiency for the NRR, catalyst designing with induced hydrophobicity and aerophilicity on its surface is fascinating.13,14 Hydrophobicity introduces a triple-phase interface with a gas diffusion layer between the electrolyte and catalyst, which entraps N2 and provides more accessible region for the NRR. In the absence of the triple-phase interface, the catalyst surface is entirely in contact with the electrolyte making proton adsorption more susceptible and enhancing the HER.15 Despite HER suppression, the poor solubility of N2 in aqueous solution16 hinders the catalyst performance in the NRR.17 It means that N2 fixation will occur effectively in the case of fast diffusion of N2 towards the catalyst and the N2 bubbles formed on the catalyst surface should adhere purposely, showing an aerophilic effect.18 All in all, the catalyst should be designed with a hydrophobic and aerophilic nature for selective NRR. Aiming at this, tetrafluoroethylene (PTFE) or alkanethiols are used in some catalysts to improve hydrophobicity.19,20 This catalyst architecture accelerates the mass transport process of gas but the addition of an extra insulating layer reduces the electron transport and slows down the electrochemical process. Therefore, the challenge is to improve gas diffusion without hampering the active catalytic sites. More importantly, besides this catalyst improvisation, it is crucial to realize the actual mechanism ongoing in these types of materials such that it could be further vitalized towards a more efficient NRR.

In general, most of the electrocatalytic processes take place in the transition metal centre, if present, where the metal–ligand moiety is present in the first coordination sphere (CS). Although the first CS is directly involved in the catalytic process, the role of the second CS is sometimes found to be exorbitant. Second coordination spheres, usually stabilize by H-bonding, π–π interactions, chemical bonding and even the electrostatic effect of the neighbouring atoms/molecules, can alter the reaction pathway in many ways. This leads to enhanced reaction kinetics via facilitating the diffusion of reactant molecules to the catalyst surface or the adsorption of the intermediates during the reaction's progress. In recent years, there have been a few reports, where the catalysts are designed based on the presence of second CSs to enhance the O2, N2 and CO2 reduction reactions.21–23 Unfortunately, most of the literature has not shown the detailed reaction pathway on how the protonation takes place after they got adsorbed on the catalyst surface, which is a major gap in previous studies. It is therefore indispensable to design an electrocatalyst for the NRR, where the second CS would help N2 diffusion and it gets adsorbed in the first CS, where the water concentration is reasonably better and sufficient for protonation during its reduction.

Therefore, to address this long-standing issue, for the first time, we have introduced the role of primary and secondary spheres in N2 adsorption and its subsequent reduction for the synthesis of ammonia. Herein, we have synthesized fluorinated copper phthalocyanine grafted on exfoliated graphene (F-CuPc-G) as a catalyst, which, as expected, is found to be highly hydrophobic and aerophilic in nature in aqueous electrolytes. The catalyst is engineered in such a way that fluorine moieties are available in the second CS (hydrophobic layer) and the active metal center is in the first CS. Therefore, there occurs fast N2 diffusion due to the formation of a triple-phase interface in the second CS and also active sites in the first CS available for the adsorption of N2 molecules and their subsequent reduction, leading to an efficient eNRR. Along with its hydrophobic nature, the catalyst is found to be aerophilic, which adheres gaseous N2 molecules on the catalyst surface and slows down the HER process, which is almost unavoidable in aqueous electrolytes and promotes the ammonia conversion.8 Moreover, a theoretical study demonstrates how hybridizing F-CuPc with graphene provides better electronic conductivity by altering the electronic structure of active sites and more density of states (DOS) near the Fermi level, facilitating better electronic transport from active sites to the electrode surface. Thus, this work highlights the importance of the co-existence of hydrophobic and aerophilic nature in the catalyst that assists in suppressing the HER and accelerating the NRR with high faradaic efficiency.

Results and discussion

The wettability of rough surfaces in nature remains an exciting subject for exploration and inspiration of many innovations. Surface roughness in combination with the chemistry of materials determines the wettability of the surface.24 For the conversion of a gaseous molecule in aqueous medium at the solid catalyst surface, the concept of hydrophobicity plays a crucial role in terms of gas entrapment at the solid–liquid–gas (three-phase) interface. Thus, getting inspiration from biological surfaces like rose petals (Fig. 1a) having nanostructures as seen in their field emission scanning electron microscopy (FESEM) image (Fig. 1b) responsible for their hydrophobicity, we added some functional groups to the existing catalyst (copper phthalocyanine: CuPc) specifically to improve its hydrophobicity and for swift N2 diffusion and entrapment on the catalyst active site. Here, we synthesized fluorinated copper phthalocyanine (F-CuPc) using a reported method25 as explained in the Experimental section and grafted it on electrophoretically exfoliated graphene as an underlying conductive support. The catalyst F-CuPc-G displayed a structural resemblance to that of rose petals as shown in FESEM (Fig. 1c). Hence, as expected, the catalyst was highly hydrophobic (Fig. 1d), as can be seen from the persistent water droplet over the catalyst surface similar to that on the rose petals (Fig. 1a). During the eNRR, the CuPc catalyst was all over flooded with electrolytes, as shown in Fig. 1e and hence, there occurred slow diffusion of N2 and was more prone to the HER. On the other hand, the as-synthesized catalyst (F-CuPc-G) manifested a better way for N2 diffusion owing to the hydrophobic layer, dominant at the three-phase interface (Fig. 1f) which was attributed to the presence of fluorine edge-functionalities that constitute the secondary sphere of the catalyst F-CuPc-G. Once the N2 molecules swept through the secondary sphere, their adsorption over the primary coordination sphere, followed by the subsequent protonation, got reinforced due to the decent availability of protons in the local hydrophilic environment of the first CS. This not only helped to mitigate the competitive HER over the catalyst active site successfully but also provided a route to justify the source of protons required for the protonation steps of the NRR despite the hydrophobic nature of the catalyst. The extent of hydrophobicity in the primary material as well as the control samples was analyzed by contact angle measurement. Fig. 2a–d show that with fluorination and graphitization the hydrophilic behaviour of CuPc turned to hydrophobic in F-CuPc-G as the contact angle changed from 71.6° to 137°.26 The high surface roughness and low surface energy attained due to the presence of graphene and fluorine, respectively, were synergistically responsible for the hydrophobicity of F-CuPc-G.20 Moreover, these properties imparted the characteristic of water adhesion to F-CuPc-G, similar to rose petal adhesion (Fig. S1). The adhesive ability of the material keeps the water in contact, which is helpful for the protonation step in the NRR process. Our hydrophobic and adhesive catalyst (F-CuPc-G) possessed another significant feature, that is the aerophilic property, which presented a strong N2 adsorption layer and interaction of the nitrogen bubbles with the catalyst as compared to CuPc as seen in Fig. 2e. The hydrophobic catalyst with aerophilic nature provided rapid N2 transport and kept it trapped even under the electrolyte flooding state (Fig. 2f).27 It was revealed that nitrogen bubbles were in the “pinning” state28i.e., aerophobic on CuPc and CuPc-G (Fig. 2g and h), while its angle decreased on the surface of F-CuPc and F-CuPc-G showing that the aerophilicity induced N2 gas layer on the catalyst slowed down the HER process on the active site, due to less availability of H+ ions (Fig. 2i and j).
image file: d3sc03002d-f1.tif
Fig. 1 (a) Optical image of rose petals with water droplets, FE-SEM image of (b) a rose petal, (c) F-CuPc-G, (d) optical image of F-CuPc-G with an electrolyte droplet, (e) scheme representing the slow diffusion of N2 in the electrolyte in the case of CuPc, and (f) fast diffusion of N2 in the electrolyte in the case of F-CuPc-G. The yellow, grey, blue, red, green and pink spheres represent the hydrogen, carbon, nitrogen, oxygen, fluorine and copper atoms.

image file: d3sc03002d-f2.tif
Fig. 2 The contact angles of electrolyte bubbles measured on the surface of (a) CuPc, (b) CuPc-G, (c) F-CuPc, and (d) F-CuPC-G; (e) scheme illustrating the aerophobic behaviour of CuPc and aerophilic behaviour of F-CuPc-G, (f) overall scheme showing the three phase interface of F-CuPc-G presenting hydrophobicity as well as aerophilicity, and the contact angles of N2 bubbles measured on the surface of (g) CuPc, (h) CuPc-G, (i) F-CuPc, and (j) F-CuPc-G.

The synthesized F-CuPc material was confirmed by Fourier-transform infrared spectroscopy (FTIR) and Raman analysis (Fig. S2 and S3).29 The FTIR spectrum of F-CuPc shows peaks at 603, 765, 965, and 1152 cm-1 due to hexa-decafluoro substituents, respectively, which are absent in the FTIR spectrum of CuPc, as shown in Fig. S2, indicating successful functionalization of fluorine in the CuPc matrix. However, owing to its very low conductivity in the order of 1–8 MΩ, F-CuPc receded in the NRR performance as discussed in the later sections. Therefore, to improve its conductivity F-CuPc was intercalated with graphene nanosheets using single-step electrophoretic exfoliation of graphite plates in an acidic medium to form a F-CuPc-graphene π-conjugated system as shown in Fig. S4 with three different contents (catalysts F-CuPc-G0.01, F-CuPc-G0.02 (F-CuPc-G), and F-CuPc-G0.03) and as expected, the conductivity of the hybrid F-CuPc-G catalyst improved (Fig. S5). The loading of F-CuPc into the graphene was measured by inductively coupled plasma mass spectroscopy (ICP-MS) analysis. From the calculation we found the presence of 120.96, 147 and 181.44 ppb of Cu in the catalysts F-CuPc-G0.01, F-CuPc-G0.02 (F-CuPc-G), and F-CuPc-G0.03 respectively in 0.07 mg catalyst used in the experiment (Table S1). Theoretical calculations were further carried out to determine the density of states (DOS) for the heterostructure of F-CuPc with graphene to understand the effect of underneath graphene (Fig. 3a and b). It was observed that a larger number of states occurred near the Fermi level for F-CuPc-G as compared to pristine F-CuPc, which was the reason for the higher conductivity of the F-CuPc-G system that was responsible for the better electron transfer from the catalyst to the electrode during the eNRR, whereas pristine graphene did not have the states at the Fermi level.30 Through charge density difference analysis, we visualized the charge re-distribution over F-CuPc molecules due to the underneath graphene layer (Fig. 3c).


image file: d3sc03002d-f3.tif
Fig. 3 (a) Optimised model structure of the F-CuPc-G heterostructure. (b) Plot of the density of states for the F-CuPc and F-CuPc-G structures, and (c) the charge density difference plots for the F-CuPc-G system. Yellow and blue lobes indicate electron accumulation and depletion layers, respectively (isosurface value = 0.00032 e Å−3).

The FE-SEM image in Fig. S6 presented the exfoliated graphene sheets. The granular structures of F-CuPc as shown in inset of Fig. 4a were also found to be dispersed on the intercalated graphene sheets depicting the presence of F-CuPc throughout the sheet (Fig. 4a). The fine graphene sheets with F-CuPc moieties stacked between them could be observed in the transmission electron microscopy (TEM) image (Fig. S7a and b). Fig. 4b distinctly shows the encircled F-CuPc dispersed over the graphene sheets. Furthermore, the images of elemental-mapping of carbon, copper, fluorine and nitrogen were displayed vividly, presenting their distribution (Fig. 4c).


image file: d3sc03002d-f4.tif
Fig. 4 (a) FE-SEM image of F-CuPc-G; inset: F-CuPc, (b) TEM image of F-CuPc-G; (c) TEM elemental mapping of F-CuPc-G, shows the distribution of C, F, N and Cu; (d) XRD pattern and (e) Raman spectra of F-CuPc-G, F-CuPc and exfoliated graphene. (f) The full survey XPS spectra of F-CuPc-G; high resolution XPS spectra of (g) Cu 2p, and (h) F 1s in F-CuPc-G.

The X-ray diffraction (XRD) analysis (Fig. 4d) showed the characteristic peak at 2θ = 26.1° denoting the (002) plane representing the sp2 hybridization of graphitic carbon-like exfoliated graphene.31 The XRD peaks of F-CuPc reflected well in the XRD pattern of F-CuPc-G and the little shift in the peaks may be due to the π–π interactions between F-CuPc and graphene sheets. In the Raman spectra of F-CuPc-G, two distinct peaks were obtained at 1347.4 cm−1 and 1577.2 cm−1 for D and G-bands reflecting the degree of graphitization in the material as in exfoliated graphene32 (Fig. 4e). The peaks found in F-CuPc overlaid some peaks in F-CuPc-G depicting the good incorporation of F-CuPc on the graphene sheets.33 Furthermore, to get more information about the surface properties of F-CuPc-G, an X-ray photoelectron spectroscopy (XPS) study was conducted.34 The XPS survey spectrum of F-CuPc-G showed distinct peaks for C, N, O, F and Cu confirming their presence in the hybrid material (Fig. 4f). The high-resolution XPS spectra of Cu (Fig. 4g) depicted peaks centered at 935.2, 955.1 and 943.5 eV for Cu 2p3/2, Cu 2p1/2 and Cu satellite peaks respectively.35 The F 1s peak at 688 eV represents the C–F bonds (Fig. 4h). The deconvoluted C 1s and N 1s spectra are shown in Fig. S8 and S9.36 All these characterisation methods confirm the successful incorporation of F-CuPc into the 2D layered graphene sheets to form F-CuPc-G. The N2 adsorption/desorption study is crucial to the catalytic efficiency of a porous catalyst. To analyse the surface area and pore size distribution (PSD) of the samples (Fig. S10a–c), the isotherms were obtained where the hybrid material possessed a surface area of 210 m2 g−1 with a PSD of 3.8 nm. Besides this, the exfoliated graphene and F-CuPc complex have a surface area of 38 m2 g−1 and 3 m2 g−1 along with a PSD of 1.5 nm, 3.8 nm and 2.7 nm, 7.4 nm, respectively. In the isotherm of F-CuPc, a very sharp increase of N2 adsorption/desorption in the higher pressure region was observed, indicating the presence of interparticle porosity (mesoporous) of the complex. In the isotherm of exfoliated graphene both micro- and meso-pores are present. The electrochemical exfoliation technique in the presence of pi electron containing F-CuPc facilitated the exfoliation of graphene sheets through pi–pi interaction. Due to incorporating the F-CuPc complex into the pi–pi stacked layered graphene material, the pore size distribution region shifted towards the mesoporous region and the surface area was also increased. The mesoporous hybrid structure increased the accessibility of the N2 containing electrolyte towards the catalytically active centre of the catalyst and consequently increased the catalytic efficiency. Also, the fast kinetics between the electrode and electrolyte was possible due to low internal resistance in F-CuPc-G as suggested by the Nyquist plot (Fig. S11). The electrochemical surface area (ECSA) was also found to be maximum for F-CuPc-G as it had the highest value of double layer capacitance (Cdl) as shown in Fig. S12.

The eNRR is mainly governed by a few parameters that primarily include nitrogen adsorption followed by its bond polarizability and protonation steps to yield NH3. The conductivity of an electrocatalyst and sufficiency of d-electrons enable a better bond polarizability of N2 by the facile electron donation and π-back donation strategy.37 This was the reason why we chose CuPc with the Cu active site coordinated with four N atoms and graphitized the material by means of electrophoretic exfoliation to obtain CuPc-G. Nevertheless, the main difficulty in aqueous electrolytes is that, although the proton sufficiency enables a better N2 protonation, there remains a predominant probability of proton adsorption over N2 adsorption. This, in turn, retards the yield rate and faradaic efficiency of NH3 synthesis. It is challenging to simultaneously play with both adsorption and protonation parameters to elevate the overall reaction kinetics. This was the primary focus of this work to enhance the N2 sufficiency around the active center without any detrimental effect on the subsequent NRR steps, which was successfully brought about by incorporating fluorine as a part of the active catalyst. The edge fluorine atoms induced local hydrophobicity and aerophilicity at the three-phase interface that made way for easy N2 diffusion across the secondary coordination sphere of the catalyst. Emphatically, the aerophilicity of the material leads to the formation of an N2 gas diffusion layer over the catalyst, which bypassed the tendency of water molecules to make their way to the catalyst before proton adsorption. So, here the role of the primary coordination sphere came into play, wherein the local hydrophilic environment provided protons to the Cu-adsorbed-N2 molecules to proceed further the protonation steps as schematically elaborated in Fig. 5a and b. The Gibbs free energy calculations of F-CuPc-G revealed a negative free energy and hence spontaneous N2 adsorption on the active site while disfavouring the high positive energy (2.11 eV) driven proton adsorption. In fact, the first protonation step (*NN to *NNH), which is the potential determining step for the NRR, involved 1.88 eV, better than the proton adsorption on the active site (Fig. 5c). This free energy value for the HER was highest for F-CuPc-G compared to CuPc and F-CuPc inherently restricting the HER (Fig. 5d). Thus, our catalyst designing played a dual role of promoting N2 adsorption and its first protonation despite the locally available protons in the primary coordination sphere of the catalyst F-CuPc-G.


image file: d3sc03002d-f5.tif
Fig. 5 (a) Scheme representing the limitations and improvements of different catalysts, (b) scheme showing the mechanistic pathway of N2 diffusion through primary and secondary spheres in F-CuPc-G, (c) free energy profile of N2 and NNH adsorption on F-CuPc-G, (d) free energy profile of the HER for F-CuPc-G, F-CuPc and CuPc systems. The bar plot diagrams illustrating the FE and yield rate (e) at different potentials for F-CuPc-G, and (d) for different catalysts.

This approach is expected to have a broader impact on the FE of ammonia production by the electrocatalyst during the NRR. In order to verify it, the electrochemical process was carried out in a three-electrode system in a 1.0 M Na2SO4 electrolyte. Before the electrochemical studies, the N2 gas was trapped in 0.5 M H2SO4 and 0.1 M KOH to free it from any kind of gaseous impurity. The trapped gas was hence checked with the gas chromatography technique, as shown in Fig. S13 and then purged into the electrolyte for the experiments. The yielded NH3 in each experiment was quantified by using a UV-vis spectrophotometer by the indophenol blue method using the calibration curve obtained with a known concentration of NH3 (Fig. S14 and S15). This was also done for acid and base trap solutions for the respective adventitious NH3 and NOx impurities to check the purity of N2 feed gas (Fig. S16–S21). For both the trap solutions, no detectable absorbance was obtained under the UV-visible spectrophotometer (Fig. S22 and S23). The side product hydrazine (N2H4) was detected by the Watt and Chrisp method38 (Fig. S24 and S25). The fast N2 diffusion by fluorination of CuPc and improvement in the conductivity by graphitization make F-CuPc-G a contender for an efficient NRR that can be seen from change in the current densities in linear sweep voltammetry (LSV) in Ar and N2 saturated atmospheres (Fig. S26). Over the potential range, 0.0 to −0.4 V vs. RHE chronoamperometry (CA) was run for an hour at each potential as shown in Fig. S27. The corresponding aliquot electrolyte was incubated for 2 hours for NH3 and 15 minutes for N2H4 detection using their respective coloring agents as explained in the ESI.39 Then, UV-vis spectra were taken for NH3 and N2H4 at different potentials (Fig. S28 and S29). The bar plot (Fig. 5e) presented the maximum NH3 yield of 80.58 μmol h−1 mgcat.−1 at −0.2 V and highest FE (49.3%) at −0.3 V vs. RHE using eqn (2) and (3) and the data have been compiled in Table S2. As expected, this FE value indicated a successful suppression of the HER (Fig. S30) at the active site, which is an essential requirement for improving the NRR kinetics in ambient aqueous electrolytes. There was no trace of any peak in the UV-visible spectra at 460 nm corresponding to the side product (N2H4) formation as determined by the Watt and Chrisp method,38 which forecasted the material selectivity towards NH3 synthesis over the entire potential range. Such high ammonia production could be attributed to the fast mass transfer of N2 due to the induced hydrophobic–aerophilic nature of the catalyst. To further confirm the role of F and underlying graphene support, various control samples such as CuPc, CuPc-G and F-CuPc were investigated (Fig. S31) for one-hour CA tests. Fig. 5f clearly reveals that F-CuPc-G wins the league with the highest NH3 yield and FE (Fig. S32 and Table S3) and also, the comparison with the literature is given in Table S4. F-CuPc acts as a n-type semiconductor40 that has low intrinsic conductivity. Therefore, the graphitization of F-CuPc improved the proton coupled electron transfer for nitrogen fixation that was anyway reflected in the electrocatalytic performance. Also, the amount of F-CuPc in graphene solution was varied (0.01 mmol, 0.02 mmol and 0.03 mmol) to optimize its electrocatalytic NRR response and the optimized F-CuPc-G (0.02 mmol) showed the highest FE and yield (Fig. S33 and S34).

Besides this, various blank CA experiments for F-CuPc-G with Ar and N2 saturation and at open circuit potential (OCP) in 1 M Na2SO4 were performed to check that the NH3 produced during the experiment was only from the feed N2 gas and not from any other nitrogenous impurity. Almost no NH3 was detected in either case, as can be seen from the spectra (Fig. S35) and the bar plot (Fig. S36). However, the active material F-CuPc-G contained N atoms coordinated to the metal active center. This adds to the benefit of the doubt whether the material N contributes to the NH3 so formed. This necessitates the verification of ammonia production using an isotope labeling experiment. In this, the distinguishable signals with different coupling constants were obtained in 1H-NMR for 14NH4+ and 15NH4+. Here, the triplet with 52 Hz and the doublet with 72 Hz were attributed to 14NH4+ and 15NH4+, respectively in 1H-NMR (Fig. 6a). Furthermore, quantitative analysis was done using maleic acid using eqn (4) in the ESI. The FE calculated from this method was almost similar in both cases of 14N2 and 15N2 as the feed gas.


image file: d3sc03002d-f6.tif
Fig. 6 (a) Isotope labelling experiment displaying 1H NMR spectra obtained after the NRR in 1 M Na2SO4 with 14N2 (pink) and 15N2 (blue) as feed gases, (b) Free energy profile of the NRR for F-CuPc-G, (c) scheme showing the reaction mechanism of F-CuPc-G during NRR follows alternating pathway. (d and e) In situ ATR-FTIR spectra obtained for F-CuPc-G.

Besides noteworthy NRR activity, it is very important that the catalyst be stable and maintain hydrophobic–aerophilic nature even after long run in the aqueous electrolyte. As fluorine functionalization of CuPc provides a hydrophobic–aerophilic interface, it is essential to check the elemental composition and contact angle after the stability test. It was found from the CA run that the catalyst was stable up to 16 h (Fig. S37) with almost no change in the curve and as expected almost the same NH3 yield rate was obtained (inset of Fig. S37). The elemental composition and structural conformation after the stability experiment were examined by XRD, Raman, XPS and elemental composition analysis (Fig. S38 and Table S5). All the characterization methods showed almost no change or insignificant change in the results obtained after the stability test compared to their respective data obtained before the stability test. This proved the robustness of our catalyst even after prolonged NRR activity. The hydrophobic–aerophilic interface was checked via contact angle measurements after a long-term stability test. It was found that the contact angle in the case of hydrophobicity changes from 137° to 132° (Fig. S39) and from 45° to 40° (Fig. S39) in the case of aerophilicity, illustrating persistent hydrophobic–aerophilic behaviour of the catalyst. This further proved the swift movement of N2 molecules by restricting the protons and the consistent N2 adsorption ability of the catalyst throughout the prolonged experimental hours. In addition, F-CuPc-G exhibited very good cyclability over five cycles of CA for one hour each under the same conditions (Fig. S40). We obtained almost similar UV-spectra for all (Fig. S41) and consequently, the NH3 yield rate and FE were repetitively reproducible (Fig. S42).

A theoretical study was performed to understand the effect of F ends and the origin of higher NRR activity of F-CuPc based systems. In this context, two model structures of CuPc with H ends and F-CuPc with F ends (Fig. S43a and b) were considered. To find the thermodynamic stability of F-CuPc, we performed molecular dynamics simulations at various temperatures (300 K, 500 K, 700 K, and 900 K).41,42 Here, we did not observe any bond breaking or distortion of the structure (Fig. S44a–d). However, we observed only a small change in the planarity of the F-CuPc molecule. Therefore, we confirmed the high thermodynamic stability of F-CuPc. Furthermore, DFT calculations were performed to elucidate the reaction mechanism of the NRR. The free energy profile of N2 and NNH indicates a small increment in the adsorption energies in F-CuPc as compared to CuPc (as shown in Fig. S45a). So, we plotted the full free energy profile of the NRR for F-CuPc (Fig. S45b) and F-CuPc-G (Fig. 6b) that gives the step of *N2 to *NNH as a potential determining step (PDS).43 It was found to be minimum for F-CuPc-G. Following theoretical studies44 from the literature and the free energy pathway (Fig. 6b), F-CuPc-G was found to obey an alternating pathway with the least Gibbs free energy. Therefore, Fig. 6c presents the full alternating pathway adopted on the catalyst, F-CuPc-G. It is more preferred over the distal pathway due to the higher adsorption energy of *NHNH than the *NNH2 intermediate.45 The synergetic roles of fluorine ends and the underneath graphene layer are influential in fast N2 diffusion and reducing the competitive hydrogen evolution reaction on the Cu site. Therefore, it is confirmed that F-CuPc-G is highly selective for nitrogen adsorption as compared to hydrogen, which leads to improve NRR activity over the HER. Thus, this work provided a strong impression that the hydrophobic–aerophilic catalyst is highly recommendable for gaseous electrocatalytic reactions.

To better understand the reaction intermediates, an in situ attenuated total reflectance-Fourier transform infrared (ATR-FTIR) spectroscopy study of F-CuPc-G during the eNRR was carried out. As shown in Fig. 6d, the interfacial FTIR spectrum in an argon medium does not show any characteristic peak. However, in a N2 saturated atmosphere, at different potentials from 0.1 to −0.5 V vs. RHE, after running chronoamperometry for 600 s, three peaks located at 1106, 1339 and 1456 cm−1 were obtained, which are assigned to N–N stretching, –NH2 wagging and H–N–H bending, respectively.46 The peak at 2326 cm−1 appeared due to the adsorption of the nitrogen molecule as N[triple bond, length as m-dash]N got polarised (Fig. 6e). Additionally, there was no peak for N2H4 observed, which is consistent with the UV-vis spectroscopy results.

Experimental section

Synthesis

Tetrafluorophthalonitrile and CuCl2 were purchased from Sigma Aldrich and were used as received. Potassium hydroxide (KOH) (analytical grade) obtained from Sigma Aldrich, hydrochloric acid (HCl), sulphuric acid (H2SO4) and potassium sulphate (Na2SO4) were purchased from Merck chemicals. All other reagents used in this study were of pure analytical grade and were used without any further purification. Graphite plates were purchased from Merck chemicals. All aqueous solution was prepared using Millipore water.

0.2 g of 1 mmol of tetrafluorophthalonitrile and 0.135 g of 1 mmol CuCl2 were added and heated at 200 °C for 3 h, giving rise to a dark blue/green solid. After it was allowed to cool down to an ambient temperature, it was vacuum filtrated and washed with water and acetone, respectively. The acetone filtrate was evaporated and F-CuPc was obtained.

The as-synthesised F-CuPc was further graphitized by a facile graphene exfoliation method. The experimental setup for this graphene synthesis is very simple equipped with two electrodes (graphite plates), an electrolyte and a DC source. Two graphite plates (1 × 1 cm2) having identical dimensions (one acts as a cathode and other as an anode) were dipped in an acidic solution containing 0.5 M H2SO4 and different amounts of F-CuPc (0.1 mmol, 0.2 mmol and 0.3 mmol named F-CuPc-G0.01, F-CuPc-G0.02, and F-CuPc-G0.03) at a static potential of 2.5 V for 30 minutes. There occurred exfoliation of graphene plate after a few minutes. The loose black flakes split up from the plate which would lead to a dark suspension of graphene flakes peeled off from the graphite plate and this plate started vanishing gradually. This black solution was then subjected to ultra-sonication for 30 to 120 minutes to get well dispersed graphene solution containing F-CuPc.

Conclusions

In summary, here we present the concept of hydrophobic–aerophilic behaviour of a F-CuP-G catalyst and the significance of the primary and secondary coordination spheres associated with it. This approach not only favoured N2 gas diffusion and adsorption across the three-phase interface but also clarified the source of protons required to carry forward the successive protonation steps despite the hydrophobic layer over the catalyst surface. This is a significant achievement of this work, as the hydrogenation step of NRR on the hydrophobic materials reported in the literature so far was unexplored. The catalyst significantly suppressed the parasitic HER process with a high NH3 yield rate of 80.53 μmol h−1 mgcat.−1 at −0.2 V and exceptionally high FE (49.3%) at −0.3 V vs. RHE. DFT and in situ FTIR studies were carried out to determine the reaction intermediates and unravel the energetically favoured associative alternating reaction pathway for the NRR on the catalyst, F-CuPc-G. This study thus opens up a new domain of research to explore the field of electrochemical NH3 synthesis with a strategy of fast N2 diffusion and HER suppression.

Data availability

All the experimental data and computational details related to this work are available either in the manuscript or in the ESI.

Author contributions

S. Bhardwaj, S. K. Das and R. S. Dey designed the work. S. K. Das synthesized the catalyst. S. Bhardwaj has done all the material charaterizations and electrochemical test. S. Bhardwaj has written the manuscript. S. K. Das and A, Biswas have assisted S. Bhardwaj in writing and mechanism study. S. Kapse and R. Thapa have carried out all the theoretical calculation and corresponding writing part. R. S. Dey supervised the work and corrected the manuscript. All the authors have checked the final version of the manuscript and approve it for submission.

Conflicts of interest

There are no conflicts to declare.

Acknowledgements

S. B., S. K. D. and A. B. acknowledge INST Mohali for providing fellowship and instrumental support. R. S. D. acknowledges DST SERB, India (Grant No. CRG/2020/005683) for financial support. RT thanks the Board of Research in Nuclear Sciences (BRNS), India, for financial support (Grant No. 37(2)/20/14/2018-BRNS/37144), Science and Engineering Research Board (SERB), India, for financial support (Grant No. CRG/2021/000620) and National Supercomputer Mission (NSM), India for financial support (Ref No. DST/NSM/R&D_HPC_Applications/2021/19). The authors thank the High-Performance Computing Center, SRM Institute of Science and Technology for providing computational facilities. We acknowledge Dr Uttam Manna and his student Ms Manideepa Dhar from IIT, Guwahati for the force adhesion experiment.

References

  1. B. Cui, J. Zhang, S. Liu, X. Liu, W. Xiang, L. Liu, H. Xin, M. J. Lefler and S. Licht, Green Chem., 2017, 19, 298–304 RSC.
  2. A. Biswas, S. Bhardwaj, T. Boruah and R. S. Dey, Mater. Adv., 2022, 3, 5207–5233 RSC.
  3. S. Murmu, S. Paul, A. Santra, M. Robert and U. K. Ghorai, Catal. Today DOI:10.1016/j.cattod.2022.10.020.
  4. C. Tang and S. Z. Qiao, Chem. Soc. Rev., 2019, 48, 3166–3180 RSC.
  5. L.-H. Zhang, F. Yu and N. R. Shiju, ACS Sustainable Chem. Eng., 2021, 9, 7687–7703 CrossRef CAS.
  6. H. Xu, K. Ithisuphalap, Y. Li, S. Mukherjee, J. Lattimer, G. Soloveichik and G. Wu, Nano Energy, 2020, 69, 104469 CrossRef CAS.
  7. L. Hu, Z. Xing and X. Feng, ACS Energy Lett., 2020, 2020, 430–436 CrossRef.
  8. P. Wang, T. Hayashi, Q. Meng, Q. Wang, H. Liu, K. Hashimoto and L. Jiang, Small, 2017, 13, 1601250 CrossRef PubMed.
  9. T. Wu, X. Zhu, Z. Xing, S. Mou, C. Li, Y. Qiao, Q. Liu, Y. Luo, X. Shi, Y. Zhang and X. Sun, Angew. Chem., Int. Ed. Engl., 2019, 58, 18449–18453 CrossRef CAS PubMed.
  10. Y. Luo, G. Chen, X. Chen, H. Wang, Y. Luo, G. Chen, L. Ding, X. Chen, L. Ding and H. Wang, Joule, 2019, 3, 279–289 CrossRef CAS.
  11. J. Li, S. Chen, F. Quan, G. Zhan, F. Jia, Z. Ai and L. Zhang, Chem, 2020, 6, 885–901 CAS.
  12. F. Zhang, J. Chen, G. G. Wallace and J. Yang, Prog. Mater. Sci., 2023, 133, 101069 CrossRef CAS.
  13. L. Xiao, S. Zhu, Y. Liang, Z. Li, S. Wu, S. Luo, C. Chang and Z. Cui, Appl. Catal., B, 2021, 286, 119895 CrossRef CAS.
  14. J. Zheng, Y. Lyu, M. Qiao, R. Wang, Y. Zhou, H. Li, C. Chen, Y. Li, H. Zhou, S. P. Jiang and S. Wang, Chem, 2019, 5, 617–633 CAS.
  15. D. Wakerley, S. Lamaison, F. Ozanam, N. Menguy, D. Mercier, P. Marcus, M. Fontecave and V. Mougel, Nat. Mater., 2019, 18, 1222–1227 CrossRef CAS PubMed.
  16. Y. Ren, C. Yu, X. Tan, H. Huang, Q. Wei and J. Qiu, Energy Environ. Sci., 2021, 14, 1176–1193 RSC.
  17. A. Biswas, S. Kapse, B. Ghosh, R. Thapa and R. S. Dey, Proc. Natl. Acad. Sci. U. S. A., 2022, 119, e2204638119 CrossRef CAS PubMed.
  18. J. Wang, Y. Zheng, F.-Q. Nie, J. Zhai and L. Jiang, Langmuir, 2009, 25, 14129–14134 CrossRef CAS PubMed.
  19. P. Wang, T. Hayashi, Q. Meng, Q. Wang, H. Liu, K. Hashimoto and L. Jiang, Small, 2017, 13, 1601250 CrossRef PubMed.
  20. C. Yu, P. Zhang, J. Wang and L. Jiang, Adv. Mater., 2017, 29, 1703053 CrossRef PubMed.
  21. C. Tang, L. Chen, H. Li, L. Li, Y. Jiao, Y. Zheng, H. Xu, K. Davey and S.-Z. Qiao, J. Am. Chem. Soc., 2021, 143, 7819–7827 CrossRef CAS PubMed.
  22. Y. Zhou, R. Lu, X. Tao, Z. Qiu, G. Chen, J. Yang, Y. Zhao, X. Feng and K. Müllen, J. Am. Chem. Soc., 2023, 145, 3647–3655 CrossRef CAS PubMed.
  23. D. Z. Zee, M. Nippe, A. E. King, C. J. Chang and J. R. Long, Inorg. Chem., 2020, 59, 5206–5217 CrossRef CAS PubMed.
  24. B. Bhushan and M. Nosonovsky, Philos. Trans. R. Soc., A, 2010, 368, 4713–4728 CrossRef CAS PubMed.
  25. E. Kol’tsov, T. Basova, P. Semyannikov and I. Igumenov, Mater. Chem. Phys., 2004, 86, 222–227 CrossRef.
  26. K.-Y. Law, J. Phys. Chem. Lett., 2014, 5, 686–688 CrossRef CAS PubMed.
  27. S. Liu, Z. Cao, Y. Meng, Y. Li, W. Yang, Z. Chang, W. Liu and X. Sun, ACS Appl. Mater. Interfaces, 2021, 13, 26853–26860 CrossRef CAS PubMed.
  28. Z. Lu, W. Xu, J. Ma, Y. Li, X. Sun and L. Jiang, Adv. Mater., 2016, 28, 7155–7161 CrossRef CAS PubMed.
  29. I. M. Denekamp, F. L. P. Veenstra, P. Jungbacker and G. Rothenberg, Appl. Organomet. Chem., 2019, 33, e4872 CrossRef.
  30. F. Shayeganfar, J. Phys.: Conf. Ser., 2015, 640, 012028 CrossRef.
  31. J. Zhao, J. Yang, L. Ji, H. Wang, H. Chen, Z. Niu, Q. Liu, T. Li, G. Cui and X. Sun, Chem. Commun., 2019, 55, 4266–4269 RSC.
  32. H. Zhang, C. Wang, H. Luo, J. Chen, M. Kuang and J. Yang, Angew. Chem., Int. Ed., 2023, 62, e202217071 CAS.
  33. T. Zou, J. Chang, Q. Chen, Z. Nie, L. Duan, T. Guo, Y. Song, W. Wu and H. Wang, ACS Omega, 2020, 5, 12067–12072 CrossRef CAS PubMed.
  34. M. Li, N. Song, W. Luo, J. Chen, W. Jiang and J. Yang, Adv. Sci., 2023, 10, 2204579 CrossRef CAS PubMed.
  35. D. G. De Oteyza, A. El-Sayed, J. M. Garcia-Lastra, E. Goiri, T. N. Krauss, A. Turak, E. Barrena, H. Dosch, J. Zegenhagen, A. Rubio, Y. Wakayama and J. E. Ortega, J. Chem. Phys., 2010, 133, 214703 CrossRef CAS PubMed.
  36. Y. Li, B. Wang, Z. Yu, X. Zhou, D. Kang, Y. Wu, Z. Chen, C. He and X. Zhou, RSC Adv., 2017, 7, 34215–34225 RSC.
  37. A. Biswas, S. Nandi, N. Kamboj, J. Pan, A. Bhowmik and R. S. Dey, ACS Nano, 2021, 15, 20364–20376 CrossRef CAS PubMed.
  38. G. W. Watt and J. D. Chrisp, Anal. Chem., 1952, 24, 2006–2008 CrossRef CAS.
  39. A. Biswas, B. Ghosh and R. S. Dey, Langmuir, 2023, 39, 3810–3820 CrossRef CAS PubMed.
  40. H. Jiang, J. Ye, P. Hu, F. Wei, K. Du, N. Wang, T. Ba, S. Feng and C. Kloc, Sci. Rep., 2014, 4, 1–6 Search PubMed.
  41. S. Kapse, S. Narasimhan and R. Thapa, Chem. Sci., 2022, 13, 10003–10010 RSC.
  42. Z. Lu, G. Xu, C. He, T. Wang, L. Yang, Z. Yang and D. Ma, Carbon, 2015, 84, 500–508 CrossRef CAS.
  43. S. Murmu, S. Paul, S. Kapse, R. Thapa, S. Chattopadhyay, A. Nagendra, S. N. Jha, D. Bhattacharyya and U. K. Ghorai, J. Mater. Chem. A, 2021, 9, 14477–14484 RSC.
  44. Y. Li, Y. Ji, Y. Zhao, J. Chen, S. Zheng, X. Sang, B. Yang, Z. Li, L. Lei, Z. Wen, X. Feng and Y. Hou, Adv. Mater., 2022, 34, 2202240 CrossRef CAS PubMed.
  45. U. K. Ghorai, S. Paul, B. Ghorai, A. Adalder, S. Kapse, R. Thapa, A. Nagendra and A. Gain, ACS Nano, 2021, 15, 5230–5239 CrossRef CAS PubMed.
  46. Y. Yao, S. Zhu, H. Wang, H. Li and M. Shao, J. Am. Chem. Soc., 2018, 140, 1496–1501 CrossRef CAS PubMed.

Footnotes

Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3sc03002d
Authors have equal contribution.

This journal is © The Royal Society of Chemistry 2023