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

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 identication 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 N 2 adsorption, (b) competitive hydrogen evolution reaction (HER), 3 (c) catalysts having Lewis acid sites that inuence N 2 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 N 2 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 eld 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 conguration and maximum utilization of atoms. 3,11 However, this metal speci-city 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 N 2 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 N 2 in aqueous solution 16 hinders the catalyst performance in the NRR. 17 It means that N 2 xation will occur effectively in the case of fast diffusion of N 2 towards the catalyst and the N 2 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, tetrauoroethylene (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 metalligand moiety is present in the rst coordination sphere (CS). Although the rst 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, p-p 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 O 2 , N 2 and CO 2 reduction reactions. [21][22][23] Unfortunately, most of the literature has not shown the detailed reaction pathway on how the protonation takes place aer 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 N 2 diffusion and it gets adsorbed in the rst CS, where the water concentration is reasonably better and sufficient for protonation during its reduction.
Therefore, to address this long-standing issue, for the rst time, we have introduced the role of primary and secondary spheres in N 2 adsorption and its subsequent reduction for the synthesis of ammonia. Herein, we have synthesized uorinated copper phthalocyanine graed 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 uorine moieties are available in the second CS (hydrophobic layer) and the active metal center is in the rst CS. Therefore, there occurs fast N 2 diffusion due to the formation of a triple-phase interface in the second CS and also active sites in the rst CS available for the adsorption of N 2 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 N 2 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 eld emission scanning electron microscopy (FESEM) image ( Fig. 1b) responsible for their hydrophobicity, we added some functional groups to the existing catalyst (copper phthalocyanine: CuPc) specically to improve its hydrophobicity and for swi N 2 diffusion and entrapment on the catalyst active site. Here, we synthesized uorinated copper phthalocyanine (F-CuPc) using a reported method 25 as explained in the Experimental section and graed 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 ooded with electrolytes, as shown in Fig. 1e and hence, there occurred slow diffusion of N 2 and was more prone to the HER. On the other hand, the assynthesized catalyst (F-CuPc-G) manifested a better way for N 2 diffusion owing to the hydrophobic layer, dominant at the three-phase interface (Fig. 1f) which was attributed to the presence of uorine edge-functionalities that constitute the secondary sphere of the catalyst F-CuPc-G. Once the N 2 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 rst 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 uorination 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 uorine, 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 signicant feature, that is the aerophilic property, which presented a strong N 2 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 N 2 transport and kept it trapped even under the electrolyte ooding state (Fig. 2f). 27 It was revealed that nitrogen bubbles were in the "pinning" state 28 i.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 N 2 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).
The synthesized F-CuPc material was conrmed by Fouriertransform 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-decauoro substituents, respectively, which are absent in the FTIR spectrum of CuPc, as shown in Fig. S2, † indicating successful functionalization of uorine in the CuPc matrix. However, owing to its very low conductivity in the order of 1-8 MU, 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-CuPcgraphene p-conjugated system as shown in Fig. S4 † with three different contents (catalysts F-CuPc-G 0.01 , F-CuPc-G 0.02 (F-CuPc-G), and F-CuPc-G 0.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-G 0.01 , F-CuPc-G 0.02 (F-CuPc-G), and F-CuPc-G 0.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).
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 ne 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 elementalmapping of carbon, copper, uorine and nitrogen were displayed vividly, presenting their distribution (Fig. 4c).
The X-ray diffraction (XRD) analysis (Fig. 4d) showed the characteristic peak at 2q = 26.1°denoting the (002) plane representing the sp 2 hybridization of graphitic carbon-like exfoliated graphene. 31 The XRD peaks of F-CuPc reected well in the XRD pattern of F-CuPc-G and the little shi in the peaks may be due to the p-p 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 Gbands reecting the degree of graphitization in the material as in exfoliated graphene 32 (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 conrming 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 2p 3/ 2 , Cu 2p 1/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 conrm the successful incorporation of F-CuPc into the 2D layered graphene sheets to form F-CuPc-G. The N 2 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 m 2 g −1 with a PSD of 3.8 nm. Besides this, the exfoliated graphene and F-CuPc complex have a surface area of 38 m 2 g −1 and 3 m 2 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 N 2 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 shied towards the mesoporous region and the surface area was also increased. The mesoporous hybrid structure increased the accessibility of the N 2 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 (C dl ) 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 NH 3 . The conductivity of an electrocatalyst and sufficiency of d-electrons enable a better bond polarizability of N 2 by the facile electron donation and p-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 N 2 protonation, there remains a predominant probability of proton adsorption over N 2 adsorption. This, in turn, retards the yield rate and faradaic efficiency of NH 3 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 N 2 sufficiency around the active center without any detrimental effect on the subsequent NRR steps, which was successfully brought about by incorporating uorine as a part of the active catalyst. The edge uorine atoms induced local hydrophobicity and aerophilicity at the three-phase interface that made way for easy N 2 diffusion across the secondary coordination sphere of the catalyst. Emphatically, the aerophilicity of the material leads to the formation of an N 2 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-N 2 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 N 2 adsorption on the active site while disfavouring the high positive energy (2.11 eV) driven proton adsorption. In fact, the rst 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 N 2 adsorption and its rst protonation despite the locally available protons in the primary coordination sphere of the catalyst F-CuPc-G.
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 Na 2 SO 4 electrolyte. Before the electrochemical studies, the N 2 gas was trapped in 0.5 M H 2 SO 4 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 NH 3 in each experiment was quantied by using a UV-vis spectrophotometer by the indophenol blue method using the calibration curve obtained with a known concentration of NH 3 ( Fig.   S14 and S15 †). This was also done for acid and base trap solutions for the respective adventitious NH 3 and NO x impurities to check the purity of N 2 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 (N 2 H 4 ) was detected by the Watt and Chrisp method 38 (Fig. S24 and S25 †). The fast N 2 diffusion by uorination 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 N 2 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 NH 3 and 15 minutes for N 2 H 4 detection using their respective coloring agents as explained in the ESI. † 39 Then, UV-vis spectra were taken for NH 3 and N 2 H 4 at different potentials ( Fig. S28 and S29 †). The bar plot (Fig. 5e (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 UVvisible spectra at 460 nm corresponding to the side product (N 2 H 4 ) formation as determined by the Watt and Chrisp method, 38 which forecasted the material selectivity towards NH 3 synthesis over the entire potential range. Such high ammonia production could be attributed to the fast mass transfer of N 2 due to the induced hydrophobic-aerophilic nature of the catalyst. To further conrm 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 NH 3 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 semiconductor 40 that has low intrinsic conductivity. Therefore, the graphitization of F-CuPc improved the proton coupled electron transfer for nitrogen xation that was anyway reected 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 N 2 saturation and at open circuit potential (OCP) in 1 M Na 2 SO 4 were performed to check that the NH 3 produced during the experiment was only from the feed N 2 gas and not from any other nitrogenous impurity. Almost no NH 3 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 benet of the doubt whether the material N contributes to the NH 3 so formed. This necessitates the verication of ammonia production using an isotope labeling experiment. In this, the distinguishable signals with different coupling constants were obtained in 1 H-NMR for 14  , respectively in 1 H-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 14 N 2 and 15 N 2 as the feed gas.
Besides noteworthy NRR activity, it is very important that the catalyst be stable and maintain hydrophobic-aerophilic nature even aer long run in the aqueous electrolyte. As uorine functionalization of CuPc provides a hydrophobic-aerophilic interface, it is essential to check the elemental composition and contact angle aer 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 Fig. 6 (a) Isotope labelling experiment displaying 1 H NMR spectra obtained after the NRR in 1 M Na 2 SO 4 with 14 N 2 (pink) and 15 N 2 (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. NH 3 yield rate was obtained (inset of Fig. S37 †). The elemental composition and structural conformation aer 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 insignicant change in the results obtained aer the stability test compared to their respective data obtained before the stability test. This proved the robustness of our catalyst even aer prolonged NRR activity. The hydrophobic-aerophilic interface was checked via contact angle measurements aer 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 swi movement of N 2 molecules by restricting the protons and the consistent N 2 adsorption ability of the catalyst throughout the prolonged experimental hours. In addition, F-CuPc-G exhibited very good cyclability over ve 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 NH 3 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 nd 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 conrmed the high thermodynamic stability of F-CuPc. Furthermore, DFT calculations were performed to elucidate the reaction mechanism of the NRR. The free energy prole of N 2 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 prole of the NRR for F-CuPc (Fig. S45b †) and F-CuPc-G (Fig. 6b) that gives the step of *N 2 to *NNH as a potential determining step (PDS). 43 It was found to be minimum for F-CuPc-G. Following theoretical studies 44 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 *NNH 2 intermediate. 45 The synergetic roles of uorine ends and the underneath graphene layer are inuential in fast N 2 diffusion and reducing the competitive hydrogen evolution reaction on the Cu site. Therefore, it is conrmed 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 reectance-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 N 2 saturated atmosphere, at different potentials from 0.1 to −0.5 V vs. RHE, aer running chronoamperometry for 600 s, three peaks located at 1106, 1339 and 1456 cm −1 were obtained, which are assigned to N-N stretching, -NH 2 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^N got polarised (Fig. 6e). Additionally, there was no peak for N 2 H 4 observed, which is consistent with the UV-vis spectroscopy results.

Synthesis
Tetrauorophthalonitrile and CuCl 2 were purchased from Sigma Aldrich and were used as received. Potassium hydroxide (KOH) (analytical grade) obtained from Sigma Aldrich, hydrochloric acid (HCl), sulphuric acid (H 2 SO 4 ) and potassium sulphate (Na 2 SO 4 ) were purchased from Merck chemicals. All other reagents used in this study were of pure analytical grade and were used without any further purication. Graphite plates were purchased from Merck chemicals. All aqueous solution was prepared using Millipore water. 0.2 g of 1 mmol of tetrauorophthalonitrile and 0.135 g of 1 mmol CuCl 2 were added and heated at 200°C for 3 h, giving rise to a dark blue/green solid. Aer it was allowed to cool down to an ambient temperature, it was vacuum ltrated and washed with water and acetone, respectively. The acetone ltrate 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 cm 2 ) having identical dimensions (one acts as a cathode and other as an anode) were dipped in an acidic solution containing 0.5 M H 2 SO 4 and different amounts of F-CuPc (0.1 mmol, 0.2 mmol and 0.3 mmol named F-CuPc-G 0.01, F-CuPc-G 0.02, and F-CuPc-G 0.03 ) at a static potential of 2.5 V for 30 minutes. There occurred exfoliation of graphene plate aer a few minutes. The loose black akes split up from the plate which would lead to a dark suspension of graphene akes 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 signicance of the primary and secondary coordination spheres associated with it. This approach not only favoured N 2 gas diffusion and adsorption across the three-phase interface but also claried the source of protons required to carry forward the successive protonation steps despite the hydrophobic layer over the catalyst surface. This is a signicant achievement of this work, as the hydrogenation step of NRR on the hydrophobic materials reported in the literature so far was unexplored. The catalyst signicantly suppressed the parasitic HER process with a high NH 3 yield rate of 80.53 mmol h −1 mg cat. −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 eld of electrochemical NH 3 synthesis with a strategy of fast N 2 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. †

Conflicts of interest
There are no conicts to declare.