Catalytic production of ammonia from dinitrogen employing molybdenum complexes bearing N-heterocyclic carbene-based PCP-type pincer ligands

Mechanistic insight into the catalytic production of ammonia from dinitrogen is needed to improve the synthesis of this vital molecule. Here we study the use of samarium diiodide (SmI2) and water in the presence of molybdenum complexes that bear PCP-type pincer ligands to synthesize ammonia. The proton-coupled electron transfer during the formation of a N–H bond on the molybdenum imide complex was found to be the rate-determining step at high catalyst concentrations. Additionally, the dimerization step of the catalyst became the rate-determining step at low catalyst concentrations. We designed PCP-type pincer ligands with various substituents at the 5- and 6-positions and observed that electron-withdrawing groups promoted the reaction rate, as predicted by density functional theory calculations. A molybdenum trichloride complex that bears a trifluoromethyl group functioned as the most effective catalyst and produced up to 60,000 equiv. ammonia based on the molybdenum atom of the catalyst, with a molybdenum turnover frequency of up to 800 equiv. min−1. The findings reported here can contribute to the development of an environmentally friendly next-generation nitrogen-fixation system. The synthesis of ammonia from dinitrogen is a vital reaction. Now, ligands prepared based on density functional theory calculations are used with a molybdenum trichloride complex for the production of up to 60,000 equiv. ammonia based on the catalyst, with a molybdenum turnover frequency of 800 equiv. min−1.

Mechanistic insight into the catalytic production of ammonia from dinitrogen is needed to improve the synthesis of this vital molecule. Here we study the use of samarium diiodide (SmI 2 ) and water in the presence of molybdenum complexes that bear PCP-type pincer ligands to synthesize ammonia. The proton-coupled electron transfer during the formation of a N-H bond on the molybdenum imide complex was found to be the rate-determining step at high catalyst concentrations. Additionally, the dimerization step of the catalyst became the rate-determining step at low catalyst concentrations. We designed PCP-type pincer ligands with various substituents at the 5-and 6-positions and observed that electron-withdrawing groups promoted the reaction rate, as predicted by density functional theory calculations. A molybdenum trichloride complex that bears a trifluoromethyl group functioned as the most effective catalyst and produced up to 60,000 equiv. ammonia based on the molybdenum atom of the catalyst, with a molybdenum turnover frequency of up to 800 equiv. min −1 . The findings reported here can contribute to the development of an environmentally friendly next-generation nitrogen-fixation system. Ammonia plays an essential role globally as a raw material to produce fertilisers and nitrogen-containing materials, which include pharmaceuticals, plastics, textiles and explosives. Presently, the Haber-Bosch process, with which ~144 metric tonnes of nitrogen are converted into ammonia annually, accounts for the main method to synthesize ammonia at the industrial level 1 . However, the conversions of dinitrogen and dihydrogen gases into ammonia via the Haber-Bosch process require a high temperature and pressure. Additionally, the preparation of dihydrogen feedstocks consumes a substantial amount of fossil fuels and is accompanied by considerable carbon dioxide (CO 2 ) emission 2,3 .
Ammonia recently attracted attention as a candidate transporter of sustainable energy because of its ease of liquefication, which is suitable for storage and transportation; moreover, only water and dinitrogen gas are emitted via combustion [4][5][6] . Thus, the past decade witnessed a demand for an environmentally friendly and tractable production of ammonia 7,8 . Notably, extensive studies were conducted on nitrogen fixation catalysed by heterogeneous catalysts to achieve the production of ammonia at a low temperature and pressure 9-12 . Since the breakthrough reported by Yandulov and Schrock, the development of nitrogen fixation under mild reaction conditions in the In our previous study, which was based on the stoichiometric and catalytic reactions of molybdenum(III) trihalide complexes that bear a pyridine-based PNP-type pincer ligand, [MoX 3 (PNP)] (X = Cl, Br and I)in which the corresponding molybdenum(IV)-nitride complex [Mo(≡N) I(PNP)] was a key reactive intermediate in the catalytic formation of ammonia under ambient reaction conditions-the molybdenum(III) trichloride complex that bears the NHC-based PCP-type pincer ligand [MoCl 3 (PCP)] (1a; PCP = 1,3-bis((di-tert-butylphosphino)methyl) benzimidazol-2-ylidene) exhibited an excellent catalytic activity 30 . However, the corresponding molybdenum(IV)-nitride complex that bears the NHC-based PCP-type pincer ligand [Mo(≡N)I(PCP)] (2a) was not prepared.
Therefore, we first prepared 2a in this study following the same procedure as that to prepare [Mo(≡N)I(PNP)]. The reaction of 1a with 5 equiv. SmI 2 in atmospheric dinitrogen and THF for five minutes at room temperature (25 °C) afforded 2a in a 72% NMR yield (Fig. 1b). Next, via the reaction of [MoI 3 (PCP)] with KC 8 as a reductant, we isolated the nitride complex (2a) as a pure form. The detailed molecular structure of 2a was confirmed via X-ray analysis ( Supplementary Fig. 7). This experimental result, which is illustrated in Fig. 1b 13 . Employing the homogenous catalytic reaction system, dinitrogen gas can be converted into ammonia by reactions with reductants and proton sources (as the chemical reagents) in the presence of a catalytic amount of transition metal complexes as the catalysts at atmospheric pressure or room temperature or lower, for example, −78 °C (refs. [14][15][16][17]. At the early stage of these reactions, the system often exhibited a limited catalytic activity owing to the deactivation of the catalysts and/or formation of dihydrogen as a side product from the high reactivity of the reductants, as well as the proton sources 13,[18][19][20][21][22][23][24][25][26][27][28] . Very recently, we reported the catalytic formation of ammonia from dinitrogen gas by combining samarium diiodide (SmI 2 ) with water (as the one-electron reductant and the proton source, respectively) in the presence of molybdenum complexes that bear a N-heterocyclic carbene (NHC)-based PCP-type pincer ligand under ambient reaction conditions 29 (Fig. 1a) 30 . Our reaction system exhibited a high catalytic activity, and demonstrated a high production of ammonia; the amount of ammonia produced reached 4,350 equiv. based on the molybdenum atom of the catalysts (the turnover frequency (TOF) was ~120 min -1 ).
Here we designed molybdenum complexes that bear different substituted PCP-type pincer ligands based on experiments, which included the isolation of nitride complexes as key reactive intermediates, as well as kinetic and theoretical studies of the catalytic reaction. Thereafter, we investigated the catalytic activity of the molybdenum complexes in the production of ammonia under ambient reaction conditions with SmI 2 and water as the reductant and the proton source, respectively. Ultimately, our system dramatically increased the quantity of ammonia produced, as well as the reaction rate, compared with the results of our previous study 30    into [Mo(≡N)I(PNP)] proceeded via a dinuclear complex, such as the dinitrogen-bridged dimolybdenum complex as a key reactive intermediate 31 . In addition, we more recently investigated cycling between molybdenum(I)-dinitrogen and molybdenum(IV)-nitride complexes to support our proposed reaction pathway for the catalytic formation of ammonia from dinitrogen, together with density functional theory (DFT) calculations 32 . Based on our previous results, here we quantitatively confirmed the conversion of the nitride ligand in 2a into ammonia via the reaction of 2a with 5 equiv. SmI 2 and H 2 O in an argon atmosphere and THF for 15 minutes at 25 °C (Fig. 1c). These experimental results indicated that dinitrogen can be transformed into ammonia employing 2a as a reactive intermediate.
Next, we compared the catalytic activities of 1a and 2a in ammonia formation under the optimized reaction conditions (Fig. 1d). The reaction of dinitrogen (1 atm) with 28,800 equiv. SmI 2 (the reductant) and 28,800 equiv. H 2 O (the proton source) in the presence of a catalytic amount (25 nmol) of 1a in THF for four hours at 25 °C yielded 8,410 equiv. ammonia based on the molybdenum atom of the catalyst (88% yield based on SmI 2 ). The catalytic reaction that employed 2a rather than 1a as the catalyst exhibited almost the same catalytic activity. Based on these stoichiometric and catalytic reactions, we expected that the molybdenum(IV)-nitride complex (2a) would function as a key reactive intermediate to effectively promote the catalytic nitrogen.
Based on our previous [30][31][32] and present findings, we proposed the following reaction pathway for the production of ammonia employing 1a as the precatalyst and 2a as key reactive intermediate ( Fig. 1e): first, the two-electron reduction of 1a with 2 equiv. SmI 2 , as well as the ligand exchange of chloride with the iodide obtained from SmI 2 , yields the corresponding dinitrogen-bridged dimolybdenum(I) complex (A). Afterwards, the cleavage of the bridged dinitrogen ligand in A affords the corresponding molybdenum(IV)-nitride complex (2a). Subsequently, 2a is converted into the corresponding molybdenum(I)ammonia complex (B) via three reduction and protonation steps that employed imide and amide complexes. The proton-coupled electron transfer (PCET) process [33][34][35] in the reaction of dinitrogen with SmI 2 and H 2 O is generally proposed as the key steps [36][37][38] . After the formation of B, its subsequent dimerization facilitates the formation of a six-coordinated A (C). Finally, the elimination of the ammonia ligand from C regenerates A, and thereby completes the catalytic cycle. As pointed out in the previous paragraph, this proposed reaction pathway was supported by DFT calculations in the reaction that employed the corresponding molybdenum complexes that bear the PNP-type pincer ligand 32 . We previously investigated the stoichiometric reactions of the molybdenum-nitride complex that bear the PNP-type pincer ligand, [Mo(≡N)I(PNP)], as well as their kinetic isotope effect (KIE) 30 . These experimental results revealed that the transformation of 2a into B also proceeded via PCET.

Kinetic study of the catalytic reaction
To gain more insight into the reaction pathway of the formation of ammonia (Fig. 1e), we determined the kinetic parameters of the catalytic reduction of dinitrogen via the reaction with SmI 2 and H 2 O in the presence of 1a in THF at room temperature. A THF solution that contained SmI 2 (0.12 M) and H 2 O (0.12 M) was mixed in atmospheric dinitrogen (1 atm) in the presence of 6-0.5 μM 1a at room temperature. The initial rate of producing ammonia (v NH 3 ) was determined, following the ammonia yield from the catalytic reaction that was quenched at the initial stage. The reaction order with respect to 1a was obtained from the slope of the plot of log(v NH 3 ) versus log([1a]), where [1a] is the concentration of 1a (Fig. 2a). The reaction order exhibited two different kinetic regions based on the concentration of the catalyst. One was the first order, which employed the catalyst at high concentrations (region A, represented by the green line), and the other was the second order, which employed the catalyst at low concentrations (region B, represented by the orange line). Similar phenomena were observed in the ruthenium-catalysed oxidation of water in which the dimerization of the ruthenium complexes was the key step 39,40 .
To investigate the rate-determining steps of both regions, we determined the rate orders in the presence of SmI 2 and H 2 O at high ([1a] = 4.17 μM) and low ([1a] = 0.7 μM) concentrations of 1a (the catalyst). Thus, employing a high concentration of the catalyst (region A), we estimated that SmI 2 and H 2 O exhibited first-order kinetics (0.95 and 0.76 for SmI 2 and H 2 O, respectively), as estimated for 1a (Fig. 2b,c, green line). These experimental results indicated that the reaction involving the three molecules, 1a, SmI 2 and H 2 O, was the rate-determining step in region A. Such reaction steps were attributable to the reduction and protonation steps in PCET (Fig. 1e). Conversely, the estimated small rate orders in SmI 2 and H 2 O employing a low concentration of the catalyst were not really zeroth order (0.34 and 0.39 for SmI 2 and H 2 O, respectively) ( Fig. 2b,c, orange line). These small orders may be due to the result of a background reduction of water into dihydrogen 41 . Based on these results, we deduced that only the concentration of the y = 1.93x + 9.89 R 2 = 1.00 y = 1.28x + 6.10 R 2 = 0.99 Article https://doi.org/10.1038/s44160-023-00292-9 molybdenum catalyst affected the second-order reaction rate at a low concentration of the catalyst. Thus, at a low catalyst concentration, the dimerization of the molybdenum complexes accounts for the rate-determining step (Fig. 1e). Further, we measured the kinetic isotope effect of the catalytic reaction with H 2 O or D 2 O in both concentration regions. A k H /k D value of 2.61 was observed in region A; this value confirmed that the PCET process was the rate-determining step [36][37][38] . Conversely, a k H /k D value of 1.42 was observed in region B, although no kinetic isotope effect was predicted by the kinetic study because of the lack of a relationship between the proton sources and the dimerization step of the catalyst as the rate-determining step. Although the details are unclear, the k H /k D differences between the two regions are consistent with a change in the rate-determining step at low and high catalyst concentrations.

Density functional theory calculations of the bond dissociation free energies
Next, we employed two strategies, based on the kinetic study (Fig. 2), to improve the catalytic activity of the system by tuning the pincer ligand of the molybdenum catalysts. The first strategy involved accelerating the PCET process in the presence of a high concentration of the catalyst, and the other involved accelerating the dimerization process in the presence of a low concentration of the catalyst. Here we focused on adding hydrogen atoms to the nitride complex (2a) to produce an ammine complex (IIIa) via PCET, after which we computationally predicted the influence of introducing the substituents to the pincer ligand of 2a on the thermodynamic stability of the hydrogenated intermediates. We previously reported that the nitrogen-fixing activity of the molybdenum-dinitrogen complexes that bear PNP-type pincer ligands could be tuned by introducing electron-donating and withdrawing groups to the pincer ligands 42,43 .
First, we performed the DFT calculations at the B3LYP-D3 theory level [44][45][46][47][48] to evaluate the bond dissociation free energies (BDFEs) of the N-H bonds in three possible reaction intermediates, beginning with 2a, namely the imide (Ia), amide (IIa) and ammine (IIIa) complexes. BDFEs were calculated according to the reaction [MoI lost a proton and an electron simultaneously. Chirik and co-workers utilized the BDFEs of N-H bonds to discuss the photoinduced PCET-driven reductions of different transition-metal-nitrogen complexes, such as manganese-nitride 49,50 , titanium-amide 51 and cobalt-imide complexes 51 . In our previous study on the molybdenum-catalysed nitrogen fixation reaction that employed [Mo(N 2 )(PMePh 2 ) 4 ], the BDFE of the N-H bond of a molybdenum-diazenide (Mo-NNH) complex was employed to evaluate the hydrogen atom affinity 52 . Table 1   Both molecular orbital energy levels correlated with the Hammett substituent constant, and a strong electron-withdrawing substituent, such as -CF 3 , lowered (or stabilized) the LUMO energy of the nitride complex more effectively than it did the HOMO energy. The lowering of the LUMO energy of 2a is expected to improve the PCET reaction to yield the corresponding imide complex, as it will enhance the electron affinity of 2a.
NHCs coordinated to a transition-metal centre are known to serve not only as a strong σ-donor but also as a π acceptor 29,55,56 . As shown in  Fig. 3d, the α-HOMO of 1a in the doublet spin state, which corresponds to the LUMO of 2a, is responsible for the π backdonation from the singly occupied d zx orbital of Mo to the vacant p z orbital of the carbene C atom. We previously reported that this unique π-accepting ability of the NHC-based PCP ligand contributed to the high thermodynamic stability of [{Mo(0)(N 2 ) 2 (PCP)} 2 (μ-N 2 )] as an N 2 -fixing catalyst 29,58 , Therefore, we can expect that the use of the substituted PCP ligands improves the catalytic activity of the present Mo-PCP system for nitrogen fixation.
In summary, DFT calculations predicted that the introduction of strong electron-withdrawing substituents, such as -CF 3 , to the PCP ligand would effectively lower the LUMO energy level. This increases the electron affinity of 2a, which is essential to accelerate the rate-determining N-H bond formation reactions via PCET. In addition, the introduction of the electron-withdrawing substituents enhances the π-accepting ability of the PCP ligand, which leads to a solid connection between the Mo centre and the PCP ligand through π backdonation to give an improvement in the catalytic activity for nitrogen fixation.

Properties of the molybdenum complexes that bear PCP ligands
Based on the DFT results, we designed PCP-type pincer ligands that bear different substituents at positions 5 and/or 6 of the benzimidazole ring (3a-3f), as well as their corresponding molybdenum trichloride complexes [MoCl 3 (R-PCP)] (1a-1f) (Fig. 4a,b, respectively). The electronic properties of these ligands and molybdenum complexes were estimated via NMR spectroscopy of the corresponding selenocarbene compounds of the PCP-type pincer ligands (4a-4f), the infrared spectroscopy of the corresponding molybdenum(0) tricarbonyl complexes that bear the same PCP-type pincer ligands, [Mo(CO) 3 (R-PCP)] (6a-6f) 59 2a and 2e).
The chemical shift of the 77 Se NMR spectra of the selenocarbene compounds generally reveals the π-backdonation ability of the corresponding carbene ligands 54,61 ; the CO-stretching frequency of the carbonyl complexes measured via infrared spectroscopy reveals all the electronic properties of the pincer ligands, which includes those of carbene and two phosphine moieties. The redox potentials measured by cyclic voltammetry of the molybdenum(IV)-nitride complexes provide direct information on the electronic properties with regard to reactive intermediates. First, we prepared the selenocarbene compounds (4a-4f) via the reactions of the corresponding carbene ligands with elemental selenium powder (see Supplementary Section 2 for the details) and performed 77 Se NMR spectroscopy of the corresponding selenocarbene compounds in THF-d 8 (Fig. 4a). The chemical shifts of the selenocarbene compounds, which were identified as selenium atoms that were bonded to carbon atoms, are listed in Table 2. The 77 Se NMR spectroscopy of the selenocarbene compounds revealed that the enhanced π backdonation contributes to the paramagnetic shielding term with a substantial downfield 77 Se-NMR signal 54,61 . The shift of the selenocarbene compounds that introduced two methyl groups as the electron-donating groups (4b, 135.2 ppm) was more upfield than that for the non-substituted selenocarbene compound (4a, 147.1 ppm). Contrarily, the introductions of the fluoro-(4c, 170.2 ppm), chloro-(4d, 177.3 ppm) and trifluoromethyl (4e and 4f with one and two CF 3 groups, respectively, 174.3 and 196.6 ppm) groups as the electron-withdrawing groups facilitated downfield shifts. These results demonstrated that the π-backdonation ability of the carbene moiety was effectively tuned via the introduction of substituent groups in positions 5 and/or 6 of the benzimidazole ring, as predicted by DFT. Next, the molybdenum tricarbonyl complexes that bear the PCP-type pincer ligands, [Mo(CO) 3 (R-PCP)] (6a-6f) were synthesized (Fig. 4b). The treatment of the pincer ligands (3a-3f), which were generated in situ via the reactions of the precursors of the PCP-pincer ligands (5a-5f) with 1.4 equiv. potassium hexamethyldisilazide (KN(SiMe 3 ) 2 ) in toluene for one hour at room temperature and 0.9 equiv. [Mo(CO) 3 (η 6 -C 6 H 5 CH 3 )] in toluene for 18 hours at 80 °C yielded the corresponding molybdenum tricarbonyl complexes that bear the PCP-type pincer ligands, [Mo(CO) 3 (R-PCP)] (6a-6f), in yields of 8-26%. X-ray analysis confirmed that the molecular structures of these complexes were almost the same ( Supplementary  Figs. 9-14). We measured the CO-stretching frequencies of 6a-6f in a THF solution. Two signals, which corresponded to the CO-stretching frequency, were observed in all the complexes. A slightly smaller wavenumber of the CO-stretching frequency of 6b, which introduced two methyl groups (1,831 and 1,935 cm −1 ) as the electron-donating groups, than that of the non-substituted complex (6a, 1,833 and 1,936 cm −1 ) was obtained. Conversely, the complexes that contained the fluoro-, chloro-and trifluoromethyl groups as electron-withdrawing groups exhibited higher wavenumbers (6c, 1,836 and 1,939 cm −1 ; 6d, 1,838 and 1,940 cm −1 ; 6e, 1,837 and 1,940 cm −1 and 6f, 1,842 and 1,942 cm −1 ). These infrared results indicated that the tendency of the effect of the substituent groups was consistent with that observed by NMR spectroscopy for the selenocarbene compounds (4), although we observed only a small difference in the CO-stretching frequencies of the molybdenum(0) tricarbonyl complexes (6) by the introduction of the substituent groups to positions 5 and/or 6 of the benzimidazole ring.
To obtain more detailed and direct information on the electronic property with regard to reactive intermediates, we measured the cyclic voltammetry of the molybdenum(IV)-nitride complex [Mo(N)I(PCP)] (2a). In this case, an irreversible reduction wave was observed at −2.85 V versus FeCp 2 +/0 , which was estimated by differential pulse voltammetry measurements. However, in the case of the molybdenum(IV)-nitride complex that bear a CF 3 -substituted PCP-type pincer ligand, [Mo(N) I(CF 3 -PCP)] (2e), the comparable reduction potential shifted to −2.63 V versus FeCp 2 +/0 ( Supplementary Fig. 1). These experimental results indicated that the introduction of electron-withdrawing groups, such as a CF 3 group at the 5-position in the PCP-type pincer ligand, substantially contributed to the electron-deficient molybdenum centre at the molybdenum(IV)-nitride complex. Thus, the electronic properties of the substituent groups that were introduced to positions 5 and/or 6 of the benzimidazole ring strongly affected the overall electronic properties of the pincer ligands.

Catalytic activities of the molybdenum PCP complexes
We examined the catalytic activities of 1a-1f in the catalytic reduction of dinitrogen into ammonia by investigating the degree of progress of the catalytic reaction. The amount of ammonia produced from the reaction of dinitrogen with 28,800 equiv. SmI 2 as the reductant and 28,800 equiv. H 2 O as the proton source in the presence of a catalytic amount (25 nmol) each of 1a-1f in THF at 25 °C was measured in the reaction time range between five minutes and four hours (Fig. 5a). The time profiles of the amounts of ammonia in the catalytic reaction are shown in Fig. 5b,c. The amount of ammonia decreased significantly when a complex that bear two methyl groups, 1b, was employed as the electron-donating group. Conversely, the utilization of complexes that bear the electron-withdrawing groups (1c-1f), as well as 1a, accelerated the reaction in the early stage. The turnover number (TON) values, as calculated by the amount of ammonia based on the Mo atom in the reaction after 60 minutes, are presented in Table 2. Notably, the TON tended to increase as the electron-withdrawing strength of the substituents at positions 5 and/or 6 of the benzimidazole ring increased; maximum values were obtained in the presence of 1e and 1f as the catalysts. However, the TOFs, as calculated by the amount of ammonia based on the Mo atom in the reaction within five minutes (Table 2), slightly decreased when the complex (1f), which exhibited a stronger electron-withdrawing property than that of 1e was utilized. Conversely, the amount of hydrogen gas produced as a by-product in this catalytic reaction system exhibited an inverse trend against the produced amount of ammonia when employing the catalysts (Fig. 5d). Thus, the amount of hydrogen gas increased significantly when the complex that bear two methyl groups was employed as the electron-donating groups (1b). The experimental results (Fig. 5) indicated that the DFT-based predictions of the designed molybdenum complexes as to the LUMO energies of the corresponding molybdenum-nitride complexes availed an efficient method to develop very effective molybdenum complexes as catalysts for the production of ammonia under ambient reaction conditions. In the present reaction system, the rate of ammonia formation decreased over time as reactants, such as SmI 2 and H 2 O, were consumed after the start of the reaction. As a result, we consider that the reason for the decrease in the rate of ammonia formation over time is due to the consumption of reactants such as SmI 2 and H 2 O. Thus, we do not believe that the reason is due to the effect of the produced ammonia and dihydrogen or the deactivation of the catalyst. Finally, we conducted the catalytic reaction in the presence of 1e, which exhibited the highest TON and TOF values among the prepared catalysts under the present reaction conditions higher amounts of the reductant and a proton source were employed (Fig. 5e)  Article https://doi.org/10.1038/s44160-023-00292-9 larger than our previously reported one (4,350 equiv. ammonia) 30 , and the highest achieved among the catalytic activities of reported reactions that employed transition metal complexes as the catalysts.

Conclusion
In summary, we confirmed that molybdenum trichloride complexes bearing a trifluoromethyl-substituted PCP-type pincer ligand functioned as the most effective catalyst for the catalytic production of ammonia from dinitrogen under ambient reaction conditions based on the mechanistic insight afforded by the experimental results, as well as on the prediction of the DFT calculations on the reactive intermediates. In this novel reaction system, up to 60,000 equiv. ammonia were produced based on the Mo atom of the catalyst, demonstrating a TOF of up to 800 equiv. Mo min −1 . The amount of the ammonia produced in this study, as well as the production rate, were approximately one order of magnitude higher than those observed under previous reaction conditions 30 .
As described in here, we consider that the use of the BDFE(N-H) values of key reactive intermediates, such as molybdenum-imide, molybdenum-amide and molybdenum-ammine complexes estimated by DFT calculations provides a suitable and reliable predicting method to develop more effective catalysts under ambient reaction conditions. The catalysts successfully developed will be applied not only to the catalytic nitrogen fixation driven by visible light 62 and electrochemical energy, but also to the catalytic formation of nitrogen-containing organic compounds directly from nitrogen gas under mild reaction conditions 63 . We believe that these findings can contribute to the development of an environmentally friendly next-generation nitrogen-fixation system in the near future.

Methods
In a 50 ml Schlenk flask was placed a CH 2 Cl 2 solution of 1e (0.05 mM, 500 μl, 25 nmol) and the solvent was removed under reduced pressure. To the flask were added SmI 2 (THF) 2 (1.05 g, 1.92 mmol) and THF (6 ml) under N 2 . Then a THF solution (1 ml) that contained H 2 O (1.92 mmol) was added to the stirred solution in the Schlenk flask in one portion. After the addition of the water, the mixture was further stirred at 25 °C for 24 h. This procedure was then repeated twice (in total, SmI 2 (THF) 2 (1.92 mmol × 3), H 2 O (1.92 mmol × 3), 24 h × 3). Aqueous potassium hydroxide solution (30 wt%, 5 ml) was added to the reaction mixture. The mixture was evaporated under reduced pressure, and the distillate trapped in a dilute H 2 SO 4 solution (0.5 M, 10 ml). The amount of ammonia (1.42 mmol, 56,800 equiv. based on the molybdenum atom, a 74% yield based on SmI 2 (THF) 2 ) present in the H 2 SO 4 solution was determined by the indophenol method.