Remarkable catalytic activity of dinitrogen-bridged dimolybdenum complexes bearing NHC-based PCP-pincer ligands toward nitrogen fixation

Abstract

Intensive efforts for the transformation of dinitrogen using transition metal-dinitrogen complexes as catalysts under mild reaction conditions have been made. However, limited systems have succeeded in the catalytic formation of ammonia. Here we show that newly designed and prepared dinitrogen-bridged dimolybdenum complexes bearing N-heterocyclic carbene- and phosphine-based PCP-pincer ligands [{Mo(N₂ )₂ (PCP)}₂ (μ-N₂ )] (1) work as so far the most effective catalysts towards the formation of ammonia from dinitrogen under ambient reaction conditions, where up to 230 equiv. of ammonia are produced based on the catalyst. DFT calculations on 1 reveal that the PCP-pincer ligand serves as not only a strong σ-donor but also a π-acceptor. These electronic properties are responsible for a solid connection between the molybdenum centre and the pincer ligand, leading to the enhanced catalytic activity for nitrogen fixation.

ARTICLE

Received 24 Jun 2016 | Accepted 10 Feb 2017 | Published 4 Apr 2017

Aya Eizawa1, Kazuya Arashiba1, Hiromasa Tanaka2, Shogo Kuriyama1, Yuki Matsuo2, Kazunari Nakajima1, Kazunari Yoshizawa2,3 & Yoshiaki Nishibayashi1

Intensive efforts for the transformation of dinitrogen using transition metaldinitrogen complexes as catalysts under mild reaction conditions have been made. However, limited systems have succeeded in the catalytic formation of ammonia. Here we show that newly designed and prepared dinitrogen-bridged dimolybdenum complexes bearing N-heterocyclic carbene- and phosphine-based PCP-pincer ligands [{Mo(N2)2(PCP)}2(m-N2)] (1) work as so far the most effective catalysts towards the formation of ammonia from dinitrogen under ambient reaction conditions, where up to 230 equiv. of ammonia are produced based on the catalyst. DFT calculations on 1 reveal that the PCP-pincer ligand serves as not only a strong

s-donor but also a p-acceptor. These electronic properties are responsible for a solid connection between the molybdenum centre and the pincer ligand, leading to the enhanced catalytic activity for nitrogen xation.

DOI: 10.1038/ncomms14874 OPEN

Remarkable catalytic activity of dinitrogen-bridged dimolybdenum complexes bearing NHC-based PCP-pincer ligands toward nitrogen xation

1 Department of Systems Innovation, School of Engineering, The University of Tokyo, Bunkyo-ku, Tokyo 113-8656, Japan. 2 Institute for Materials Chemistry and Engineering, Kyushu University, Nishi-ku, Fukuoka 819-0395, Japan. 3 Elements Strategy Initiative for Catalysts and Batteries, Kyoto University, Nishikyo-ku, Kyoto 615-8520, Japan. Correspondence and requests for materials should be addressed to K.Y. (email: mailto:[email protected]

Web End [email protected] )or to Y.N. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 8:14874 | DOI: 10.1038/ncomms14874 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14874

Nitrogen is an essential element for all living things on earth. Since most of the nitrogen atoms on earth exist as the form of inert dinitrogen gas, the xing of molecular

dinitrogen is necessary to be utilized. The industrial dinitrogen xation system, called the HaberBosch process, plays an important role in producing ammonia from dinitrogen gas today1. The operation of the process, however, requires high temperature and high pressure, resulting in large consumption of fossil fuels1. On the other hand, nitrogenases transform dinitrogen gas into ammonia under ambient reaction conditions, where the active sites of nitrogenases include iron, molybdenum and vanadium as essential transition metals25. Studies on the active sites of nitrogenases are considered to be important to elaborate an efcient articial system for the catalytic ammonia formation from dinitrogen gas68.

Despite intensive efforts for the transformation of dinitrogen gas using transition metaldinitrogen complexes as catalysts under mild reaction conditions916, only a few systems have succeeded in the catalytic formation of ammonia from dinitrogen gas1728. In 2003, Yandulov and Schrock2932 reported the rst successful example of the catalytic conversion of dinitrogen

gas into ammonia using a molybdenumdinitrogen complex as a catalyst and in 2013 Peters and co-workers3339 reported the iron-catalysed transformation using an irondinitrogen complex as a catalyst. We also found that several dinitrogen-bridged dimolybdenum complexes such as [{Mo(N2)2(PNP)}2(m-N2)]

(2; PNP 2,6-bis(di-tert-butylphosphinomethyl)pyridine)4042

and molybdenumnitride complexes bearing PNP-type pincer ligands43 or mer-tridentate triphosphine44 worked as more effective catalysts towards ammonia formation under ambient reaction conditions, where up to 63 equiv. of ammonia were produced based on the molybdenum atom of the catalyst.

During our continuous study, we have realized three promising clues to develop more effective catalysts. The rst clue is the introduction of an electron-donating group to the pincer ligands to increase the backdonating ability of the molybdenum atom to the dinitrogen ligand. In fact, dinitrogen-bridged dimolybdenum complexes bearing the electron-donating group-substituted PNP-pincer ligands worked as more effective catalysts in our previous reaction system41. The second clue is the inhibition of the dissociation of the pincer ligand from the molybdenum atom to increase the stability of the molybdenum complex. We generally

a

PCP-type pincer ligands (PCP, phosphoruscarbenephosphorus)

P

P

cf.

N

N

N

N

P P

P

P

P

N

N

N

N

P

P

Bim-PCP[1] Im-PCP[2] Bim-PCP[2] Im-PCP[1] P = PtBu2 P = PtBu2 P = PtBu2 P = PtBu2

P = PtBu2

P = PtBu2

b

N N

P

R R

R

N

N

N2 (1 atm) Na-Hg

THF rt, 17 h

N2 (1 atm) Na-Hg

THF rt, 17 h

P

N

N

N

N

N N

N

Mo N

N N

N N

R

Mo

Cl

Cl

Cl

N

Mo

P R

R

P

P

R = H [MoCl3(Bim-PCP[1])] (3a) [{Mo(N2)2(Bim-PCP[1])}2(-N2)] (1a), 37%[{Mo(N2)2(Me-Bim-PCP[1])}2(-N2)] (1c), 46%

[{Mo(N2)2(lm-PCP[2])}2(-N2)] (1b), 53%

R = H

R = Me

R = Me

[MoCl3(Me-Bim-PCP[1])] (3c)

[MoCl3(lm-PCP[2])] (3b)

N(7)*

N(7) N(7)

N(3)

N(4)

N N

Cl

Cl

N N

P P

P P

N N

N

P

N N

N

N N

N N

N

Mo

Cl

Mo

Mo

P

c

P(2) P(2)

N(6)

Mo(1)

P(1)

Mo(1)

N(6)

N(5)

N(2)

N(1)

Mo(1)*

N(5)

N(7)*

N(3)

N(4)

N(2)

N(1)

C(1)

C(1)

P(1)

1a 1b

Figure 1 | Design and synthesis of dinitrogen-bridged dimolybdenum complexes. (a) Newly designed PCP-type pincer ligands (PCP, phosphorus carbenephosphorus). (b) Synthesis of dinitrogen-bridged dimolybdenum complexes 1a1c. (c) ORTEP drawings of 1a (left) and 1b (right). Thermal ellipsoids are shown at the 50% probability level. Hydrogen atoms and solvated molecules are omitted for clarity.

2 NATURE COMMUNICATIONS | 8:14874 | DOI: 10.1038/ncomms14874 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14874 ARTICLE

observed the dissociation of the PNP-pincer ligand after the ammonia formation in the catalytic reaction40. The third clue is the preservation of the dinitrogen-bridged dimolybdenum core to promote the catalytic ammonia formation from the coordinated dinitrogen43. Density functional theory (DFT) calculations demonstrated that the dinitrogen-bridged dimolybdenum structure plays a vital role in the protonation of a dinitrogen ligand, where one molybdenum moiety of the dinuclear molybdenumdinitrogen complex works as a mobile ligand to the other molybdenum moiety as an active site43.

Taking account of these clues, we have now planned to design an N-heterocyclic carbene- (NHC-)45,46 and phosphine-based PCP-type pincer ligand (a PCP-type pincer ligand composed of NHC and two phosphines) as a tridentate ligand in place of the so far employed PNP-type pincer ligand for preparing a new molybdenumdinitrogen complex. It is known that NHC works as a stronger electron-donating ligand than pyridine and binds to a transition metal centre more strongly than pyridine47,48. In this article, we demonstrate that dinitrogen-bridged dimolybdenum complexes bearing PCP-pincer ligands [{Mo(N2)2(PCP)}2 (m-N2)] (1) worked as effective catalysts towards ammonia formation under ambient reaction conditions, where up to 230 equiv. of ammonia were produced based on the catalyst (115 equiv. of ammonia based on the molybdenum atom of the catalyst). This is so far the most effective catalytic reduction of dinitrogen gas into ammonia under ambient reaction conditions using transition metaldinitrogen complexes as catalysts.

ResultsPreparation and characterization of 1. On the basis of our proposal, we designed two types of dinitrogen-bridged dimolybdenum complexes bearing PCP-type pincer ligands with two tert-butyl groups on each phosphorus atom (Fig. 1a). One is the complex bearing methylene linkers between the NHC skeleton and the phosphorus atom, where a similar PCP-pincer ligand bearing two phenyl groups on each phosphorus atom has quite recently been reported by Rieger and co-workers49. The other is the complex bearing ethylene linkers, where similar PCP-pincer ligands bearing two phenyl groups on each phosphorus atom have already been reported by some research groups5055.

According to our previous procedure4042, we newly prepared three dinitrogen-bridged dimolybdenum complexes bearing the PCP-type pincer ligands. Treatment of [MoCl3(PCP)] (3a,

PCP 1,3-bis((di-tert-butylphosphino)methyl)benzimidazol-2-

ylidene (Bim-PCP[1]); 3b, PCP 1,3-bis(2-(di-tert-butylphosphino)

ethyl)imidazol-2-ylidene (Im-PCP[2]); 3c, PCP 5,6-dimethyl-1,

3-bis((di-tert-butylphosphino)methyl)benzimidazol-2-ylidene (Me-Bim-PCP[1])) with 6 equiv. of NaHg in tetrahydrofuran (THF) at room temperature for 17 h under an atmospheric pressure of molecular dinitrogen gave the corresponding dinitrogen-bridged dimolybdenum complexes [{Mo(N2)2(PCP)}2 (m-N2)] (1a; PCP Bim-PCP[1], 1b; PCP Im-PCP[2], 1c;

PCP Me-Bim-PCP[1]) in 37%, 53% and 46% yields,

respectively (Fig. 1b). Detailed synthetic procedures for ligand precursors are included in Supplementary Methods and synthetic procedures for metal precursors 3a3c are included in Methods section. These dinitrogen-bridged dimolybdenum complexes were characterized by 1H and 31P{1H} NMR. Detailed molecular structures of these complexes 1a1c were determined by X-ray crystallography (Fig 1c for 1a and 1b; Supplementary Fig. 1 for 1c), which were similar to those of the dinitrogen-bridged dimolybdenum complexes bearing PNP-pincer ligands40

42. However, the bond lengths, the bond angles and dihedral angles were signicantly different between 1a and 1b according to the nature of the linkers in the PCP-pincer ligands

(Supplementary Tables 3 and 4). The bond lengths dened by Mo(1)-C(1) of 1a and 1b were 2.064(2) and 2.153(4) , respectively, and the bond angles dened by P(1)-Mo(1)-P(2) of 1a and 1b were 153.95(3) and 163.18(8), respectively. On the other hand, the dihedral angles dened by N(1)-C(1)-Mo(1)-N(5) of 1a and 1b were 82.35(5) and 43.74(7), respectively. The shortened bond length of Mo(1)-C(1) of 1a suggests the stronger p-backdonation from the molybdenum centre to the NHC unit due to the almost-perpendicular orientation of the NHC unit. On the other hand, the longer CH2 linker of Im-PCP[2] forces the twisted coordination of the NHC and is likely to weaken the p-backdonation in 1b. Further information on this topic is discussed based on DFT calculations (vide infra).

Infrared spectra of 1a1c in the solid state showed a strong absorption peak assignable to terminal dinitrogen ligands at 1,978, 1,911 and 1,969 cm 1 respectively. The single peak of each complex corresponds to the dinitrogen-bridged dimolybdenum structure, as determined by X-ray crystallography. Compared with the infrared spectrum of 1b, those of 1a and 1c showed the peak at much higher frequency due to the strong p-backdonation from the molybdenum centre to the NHC. The infrared spectra of 1a and 1c in THF solution showed one strong absorption peak assignable to terminal dinitrogen ligands at 1,979 and 1,973 cm 1, respectively. Comparison of the infrared spectra of 1a and 1c in the solid state with that in the solution state indicates that the dinitrogen-bridged dinuclear structures of 1a and 1c are preserved even in solution. Furthermore, the 15N{1H} NMR spectrum of 15N2-labelled 1a in C6D6 under 15N2 showed two singlet and one doublet signals; d 7.2 (s, MoN NMo), 13.0

(d, 1JNN 5.4 Hz, MoN N), 32.0 (br s, MoN N), which

are consistent with the dinuclear structure (Fig. 2a)40,56. In contrast, the infrared spectrum of 1b in THF solution showed two peaks assignable to terminal dinitrogen ligands, suggesting that the structure of 1b in THF is no longer the same with that in the solid state. To obtain more information on real species of 1b in the THF solution, the 15N{1H} NMR spectrum of 1b was measured in THF-d8 solution under an atmospheric pressure of 15N2 gas. The spectrum showed two doublet and two double triplet signals; d 1.2 (d, 1JNN 6.0 Hz,

MoN N(equatorial)), 20.8 (d, 1JNN 6.1Hz, MoN N(axial)),

27.8 (dt, 1JNN 6.0 Hz and

2JNP 1.6 Hz, MoN N

(equatorial)), 30.5 (dt,

1JNN 6.1 Hz and

2JNP 2.3 Hz,

MoN N(axial); Fig. 2b), demonstrating the formation of the

corresponding mononuclear dinitrogen complex [Mo(N2)3

(Im-PCP[2])] (1b0) in THF. In fact, these spectroscopic features of 1b0 in THF are consistent with those of similar mononuclear dinitrogen complexes such as mer-[Mo(N2)3L3] structures57,58. The instability of 1b in solution may be derived from the steric repulsion between the two molybdenum moieties bearing Im-PCP[2].

As shown in Fig. 1, we prepared two types of the PCP-pincer ligands, where Bim-PCP[1] has a benzimidazol-2-ylidene skeleton and Im-PCP[2] has an imidazol-2-ylidene skeleton. Unfortunately, we were unable to synthesize the corresponding two PCP-pincer ligands based on the same NHC skeleton bearing linkers of the different lengths between the NHC skeleton and each phosphorus atom (Bim-PCP[2] and Im-PCP[1] in Fig. 1a). However, we consider that the presence or absence of benzene ring has little inuence on either the thermodynamic stability of the MoN NMo structures or the electron-donating ability

of the pincer ligands for the following reasons. To assess the inuence of the benzene moiety in Bim-PCP[1] on the thermodynamic stability of the dinitrogen-bridged dimolybdenum structure, we calculated a model complex [{Mo(N2)2(Im-

PCP[1])}2(m-N2)], where the benzimidazol-2-ylidene skeleton in Bim-PCP[1] is replaced by the imidazol-2-ylidene skeleton. The optimized distance of the MoN2(bridging) bond and its bond

NATURE COMMUNICATIONS | 8:14874 | DOI: 10.1038/ncomms14874 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14874

a

N NP

Mo Mo P P

N

N

-N

P N

N

N

N

N = 15N

15N2-1a

N N

N N

N N

N N

P = PtBu2

-N

N N

N N

10

0

10 20 30 p.p.m.

b

N N

N

NNMo Mo

N

N

P

N

N N

N

N

P

P

N

15N2 (1 atm)

THF-d8, rt

N

P

Mo

N

2

N

P

N

N N N

N N'

P

N N

N'

N = 15N

P = PtBu2

N

1b

[Mo(15N2)3(Im-PCP[2])] (15N2-1b)

N'

N'

10 10 20 30 p.p.m.

0

Figure 2 | 15N{1H} NMR spectra of 1a and 1b. (a) 15N{1H} NMR spectrum of 15N2-1a in C6D6 under 15N2. (b) 15N{1H} NMR spectrum of 15N2-1b0 in THF-d8 under 15N2.

dissociation energy (BDE) are 2.105 and 18.4 kcal mol 1, respectively, both of which are almost identical to those calculated for 1a (2.108 and 18.8 kcal mol 1). The denition of BDE is described in the Methods section. On the other hand, Tuczek and co-workers52,53 have previously prepared two PCP ligands based on benzimidazol-2-ylidene and imidazol-2-ylidene, where these ligands have similar s-donating ability. Gusev59,60 has previously estimated the donor ability of various NHC ligands based on the computational evaluation of nCO (A1) of [Ni(CO)3(NHC)]. The author showed that 1,3-dimethylbenzimidazol-2-ylidene and 1,3-dimethylimidazol-2-ylidene have almost the same values of nCO (2,057 and 2,054 cm 1, respectively), suggesting that the electron-donating ability of the NHC ligands was scarcely inuenced by the difference between benzimidazol-2-ylidene and imidazol-2-ylidene. We therefore expect that the introduction of a benzene ring to the NHC skeleton little inuences either the thermodynamic stability of the MoN N

Mo structure or the electron-donating ability of the pincer ligand.

We then performed DFT calculations to elucidate how the length of linkers connecting the NHC skeleton with the PtBu2 groups in Bim-PCP[1] and Im-PCP[2] inuences the thermodynamic stability of the dinitrogen-bridged dimolybdenum structures in 1a and 1b. Figure 3 shows optimized structures of dimolybdenum complexes 1a and 1b. The N N stretching

frequencies of terminal dinitrogen ligands calculated for 1a (2,012 cm 1) and 1b (1,969 cm 1) reproduced the experimental trend. The MoN and NN distances of a terminal dinitrogen ligand are calculated to be 2.032 and 1.137 for 1a and 2.016 and1.142 for 1b, respectively, indicating that the terminal dinitrogen ligand in 1b is more activated than that in 1a. As a

result, the BDE of a MoN2(terminal) of 1a (11.9 kcal mol 1) is considerably lower than that of 1b (16.7 kcal mol 1). On the other hand, the MoN distance of the bridging dinitrogen ligand of 1a (2.108 ) is shorter than that of 1b (2.133 ). The BDE of the MoN2(bridging) bond of 1a is 18.8 kcal mol 1, which is more than twice as high as that of 1b, 9.0 kcal mol 1. The very low BDE of the MoN2(bridging) bond of 1b can be associated with the experimental observation that the dinuclear complex 1b is labile to be separated into two mononuclear complexes in solution.

Differences in thermodynamic stability of the MoN NMo

structure between 1a and 1b can be rationalized by optimized structures of the corresponding mononuclear dinitrogen complexes. Figure 4 presents space-lling models of [Mo(N2)3(PCP)]

(1a0; PCP Bim-PCP[1], 1b0; PCP Im-PCP[2]). Comparison of

the MoN distance for the equatorial dinitrogen ligand of 1a0(2.084 ) with that of 1b0 (2.041 ) suggests that 1b0 bearing Im-PCP[2] strongly binds dinitrogen at the equatorial position. Contrary to the BDEs of the MoN2(bridging) bond calculated for dinuclear complexes 1a and 1b, the BDE of the MoN2(equatorial) bond of 1b0 (21.5 kcal mol 1) is almost the same with that of 1a0 (21.2 kcal mol 1). The dramatic decrease in the MoN2(equatorial) BDE of 1b0 in the formation of the dinuclear structure can be ascribed to steric hindrance caused by the tert-butyl groups on each phosphorus atom in Im-PCP[2]. The optimized structure of 1a0 bearing Bim-PCP[1] with the methylene linkers has the PMoP bond angle of 151.2, while that of 1b0 bearing Im-PCP[2] with the ethylene linkers has the bond angle of 164.3. The extension of the CH2 linkers in 1b0 forces the tert-butyl groups on each phosphorus atom in

4 NATURE COMMUNICATIONS | 8:14874 | DOI: 10.1038/ncomms14874 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14874 ARTICLE

2.032

1.137

1.148

2.093

2.108

[afii9840]NN = 2,012 cm1

[afii9840]NN = 1,969 cm1

1a

1b

2

N

N P P

P

P

tBu2 tBu2

tBu2

tBu2

N

N

N

N

N

N

N

N

Mo Mo

PMoP = 151.2 PMoP = 164.3 1a' 1b'

Figure 4 | Space-lling models of mononuclear MoN2 complexes 1a0 and 1b0. The dashed lines represent the projection of tert-butyl groups on phosphine atoms.

N

N N N

N

2.016 1.142

1.151

N

2.172

2.133

2.052

1.139

2.022

2.235

1.154

[afii9840]NN = 1,993 cm1

Figure 3 | Optimized structures of dinuclear complexes. Bond distances are shown in . The values of nNN present the N N stretching frequencies

of terminal dinitrogen ligands. Hydrogen atoms are omitted for clarity.

Im-PCP[2] to project towards the space surrounding the equatorial dinitrogen ligand (Fig. 4). As a result, the formation of a thermodynamically stable dimolybdenum complex bearing Im-PCP[2] is encumbered by steric repulsions between tert-butyl groups in two mononuclear molybdenum units facing each other.

Catalytic nitrogen xation using 1 as catalysts. The catalytic reduction of molecular dinitrogen into ammonia using 1 as catalysts was carried out according to the following procedure of the previous method4042. To a mixture of 1 and 2,6-lutidinium triuoromethanesulfonate (96 equiv. to 1; [LutH]OTf) as a proton source in toluene was added a solution of metallocene (72 equiv. to 1) as a reductant in toluene via a syringe pump at room temperature over a period of 1 h, followed by stirring at room temperature for another 19 h under an atmospheric pressure of dinitrogen. After the reaction, the amounts of ammonia and molecular dihydrogen were determined by indophenol method61 and gas chromatography (GC), respectively. The yields of ammonia and molecular dihydrogen were calculated based on the metallocene. Typical results are shown in Table 1. In all cases, no formation of other products such as hydrazine was observed at all.

First, we carried out the catalytic reaction in the presence of 1a as a catalyst using either cobaltocene (CoCp2; Cp Z5C5H5),

decamethylchromocene (CrCp*2; Cp* Z5C5Me5), and deca

methylcobaltocene (CoCp*2) as reductants, to produce 5.7, 17.6 and 11.8 equiv. of ammonia based on the catalyst, respectively (Table 1, runs 13). In the absence of a reductant, only 0.2 equiv. of ammonia were produced based on 1a (Table 1, run 4). We have previously obtained the result that 12.2 equiv. of ammonia were produced based on the catalyst from the reaction with CrCp*2 as a reductant in the presence of [{Mo(N2)2(PNP)}2 (m-N2)] 2 as a catalyst (Table 1, run 9)40. This means that 1a promoted the catalytic nitrogen xation more effectively than 2. Separately, we conrmed the direct conversion of molecular dinitrogen into ammonia by using 15N2 gas in place of normal 14N2 gas (see Supplementary Methods for detailed procedure).

In stark contrast to the catalytic activity of 1a, 1b did not work so effectively towards the formation of ammonia under the same reaction conditions. When CoCp2, CrCp*2 and CoCp*2 were used as reductants, only 1.4, 3.2 and 2.9 equiv. of ammonia were produced based on the catalyst, respectively (Table 1, runs 57). In the absence of a reductant, 1.5 equiv. of ammonia were produced based on 1b (Table 1, run 8).

Next, we investigated the inuence of a proton source in the catalytic nitrogen xation using 1a as a catalyst. Typical results are shown in Table 2, where larger amounts of both reductant and proton source were employed in order to sharpen the difference among the results (see Supplementary Methods for the detailed procedure). The catalytic reaction using larger amounts of reductant CrCp*2 (360 equiv. to 1a) and proton source [LutH]OTf (480 equiv. to 1a) afforded 79 equiv. of ammonia based on the catalyst (Table 2, run 1). When 2-picolinium triuoromethanesulfonate ([PicH]OTf; Pic 2-picoline) was

used in place of [LutH]OTf, only a small amount of ammonia was produced based on the catalyst (Table 2, run 2). On the other hand, 2,4,6-collidinium triuoromethanesulfonate ([ColH]OTf; Col 2,4,6-collidine) worked rather effectively, where 61 equiv. of

ammonia were produced (Table 2, run 3). When a non-coordinating anion BArF4 (ArF 3,5-(CF3)2C6H3) was used in

place of OTf in [LutH]OTf, only a small amount of ammonia was produced (Table 2, run 4). These results indicate that the use of [LutH]OTf as a proton source is an essential factor to achieve the high performance of 1a as a catalyst.

NATURE COMMUNICATIONS | 8:14874 | DOI: 10.1038/ncomms14874 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14874

Table 1 | Catalytic formation of ammonia from dinitrogen gas employing 1a or 1b as a catalyst.

N2 +

1 atm

cat.(10.0 mol)

6 reductant

72 equiv.* 96 equiv.*

Run Catalyst Reductant E1/2 (V)w NH3 (equiv.)* NH3 (%)z H2 (equiv.)* H2 (%)z 1 1a CoCp2 1.15 5.70.6 24 15.21.5 42

2 1a CrCp*2 1.35 17.60.7 73 3.90.4 11

3 1a CoCp*2 1.78 11.81.0 49 3.30.8 9

4 1a none 0.2 0.8 0.1 0.3

5 1b CoCp2 1.15 1.40.1 6 16.71.9 47

6 1b CrCp*2 1.35 3.20.3 13 15.10.7 42

7 1b CoCp*2 1.78 2.91.0 12 3.01.0 8

8 1b none 1.5 6 0.1 0.3

9y 2 CrCp*2 1.35 12.2 51 4.2 12

A solution of a reductant in 4 ml of toluene was added to a mixture of the catalyst and [LutH]OTf in 1 ml of toluene at room temperature over a period of 1 h, followed by stirring at room temperature for another 19 h under 1 atm of dinitrogen gas.*Equivs based on catalystwE values of reductant in MeCN versus Ag/Ag (ref. 69).

zYields based on reductant. yref. 40.

+ 6 N

OTf

2 NH3

H

toluene rt, 20 h

Table 2 | Catalytic formation of ammonia from dinitrogen gas using 1a as a catalyst.

N2

cat. 1a(2.0 mol)

+ 6 Cr + 6 proton

source 2 NH3

toluene rt, 20 h

1 atm 360 equiv.* 480 equiv.*

Run Proton source pKaw NH3 (equiv.)* NH3 (%)z H2 (equiv.)* H2 (%)z 1 [LutH]OTf 6.77 794 66 122 62 [PicH]OTf 5.97 191 16 5811 283 [ColH]OTf 7.48 611 51 288 124 [LutH]BArF4 6.77 155 13 4215 16

A solution of CrCp* in 4 ml of toluene was added to a mixture of 1a and a proton source in 1 ml of toluene at room temperature over a period of 1 h, followed by stirring at room temperature for another 19 h under 1 atm of dinitrogen gas.*Equivs based on catalyst.

wpK value of proton source in H O (ref. 70).

zYield based on CrCp* .

N2

1 atm

+ +

Cr

OTf cat.

(1.0 mol)

toluene rt, 20 h

2 NH3 + H2

N H

1,440 equiv. 1,920 equiv.

cat.

N NP

P

N

N N

N

N

NN

N

P

N

N

N

N

N

N P

Mo Mo

N N

N

N

N

N

N

N

N N

1a 1c

NH3 H2

P

P

Mo Mo

P = PtBu2

N N

P

P

P = PtBu2

(2.00.2)102 equiv. (42%)(1.00.4)102 equiv. (14%)

NH3 H2

2.3102 equiv. (48%)1.2102 equiv. (16%)

Figure 5 | Catalytic formation of ammonia using larger amounts of a reductant and a proton source in the presence of 1a or 1c as a catalyst. A solution of CrCp*2 in 5 ml of toluene was added to a mixture of the catalyst and [LutH]OTf in 1 ml of toluene at room temperature over a period of 1 h (for 1a) or 5 h (for 1c), followed by stirring at room temperature for another 19 h (for 1a) or 15 h (for 1c) under 1 atm of dinitrogen gas. The amounts of ammonia and hydrogen (equiv.) are based on the catalyst. Yields are based on CrCp*2.

6 NATURE COMMUNICATIONS | 8:14874 | DOI: 10.1038/ncomms14874 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14874 ARTICLE

N2 + 6 e + 6 H+ cat. toluene, rt

288 equiv.

216 equiv.

1 atm

60

1c

1a

2

On the basis of the results shown in Tables 1 and 2, we carried out the catalytic reaction using much larger amounts of CrCp*2 and [LutH]OTf as a reductant and a proton source to the catalyst, respectively (see Supplementary Methods for the detailed procedure). The reaction using 1,440 equiv. of CrCp*2 as a reductant and 1,920 equiv. of [LutH]OTf as a proton source in the presence of 1a as a catalyst under ambient reaction conditions gave 2.0 102 equiv. of ammonia based on the catalyst (Fig. 5).

The catalytic activity of 1a towards the formation of ammonia is an order of magnitude greater than that of [{Mo(N2)2(PNP)}2 (m-N2)] 2 (up to 23 equiv. of ammonia based on 2). Ammonia was obtained in 42% yield together with molecular dihydrogen (14% yield). Furthermore, a higher catalytic activity has been achieved when 1c was used as a catalyst, where up to 2.3

102 equiv. of ammonia based on the catalyst were produced under similar reaction conditions (Fig. 5). We have not yet obtained the exact reason why 1c worked as the more effective catalyst than 1a, but we consider that the introduction of two methyl groups to the benzimidazol-2-ylidene skeleton in the PCP-pincer ligand may increase the backdonating ability of molybdenum centres to the coordinated dinitrogen ligand and activated the terminal dinitrogen ligands more strongly than 1a. The infrared spectrum of 1c in the solid state showed a strong absorption peak assignable to terminal dinitrogen ligands at 1,969 cm 1, which is lower than that of 1a (vide supra). Previously, we reported that the introduction of electron-donating groups such as methyl and methoxy groups to the pyridine ring of the PNP-pincer ligand in 2 markedly enhanced the catalytic activity under the same reaction conditions41.

The time prole of the catalytic reactions using 1a and 1c as catalysts was monitored (see Supplementary Methods for the detailed procedure). Typical results are shown in Fig. 6 together with the time prole using 2 as a catalyst41,42. The turnover frequency (TOF) for ammonia formation, which was determined as mols of ammonia (based on the catalyst) produced in initial 1 h, was 42 h 1 for 1a and 53 h 1 for 1c. The TOFs for ammonia using 1a and 1c are ca. 2.5 and 3.1 times larger than that using 2 (17 h 1), respectively. This result indicates that the dinitrogen-bridged dimolybdenum complexes bearing PCP[1]-type pincer ligands such as 1a and 1c have the effective performance not only on the catalytic activity but also on the rate for ammonia formation.

For comparison of the stability of the dinitrogen-bridged dimolybdenum complexes bearing PCP-pincer ligands with that of PNP-pincer ligands, we carried out the following catalytic reactions using 1a and 2 as catalysts. After the formation of 47 equiv. of ammonia based on the catalyst from molecular dinitrogen following the same procedure of the time prole experiment using 1a, the same amounts of [LutH]OTf and CrCp*2 were further added at room temperature, and the mixture was stirred for another 2 h to afford a total of 69 equiv. of ammonia based on the catalyst (Fig. 7; see Supplementary Methods for the detailed procedure). In this reaction system, 22 equiv. of ammonia were produced from a further reaction of molecular dinitrogen with excess amounts of a proton source and a reductant. This experimental result indicates that the active species derived from 1a remain even after the catalytic reaction. In fact, no free PCP-pincer ligand was observed from the reaction mixture after the catalytic reaction using 1a as a catalyst, suggesting that the active species derived from 1a were still active. In sharp contrast, no additional ammonia was produced from similar treatment using 2 as a catalyst (Fig. 7), where free PNP-pincer ligand was observed from the reaction mixture after the catalytic reaction using 2 as a catalyst40. These results indicate that the stability of 1a is much improved compared with that of 2.

Comparison of PCP and PNP ligands. In this section, we compare the electronic properties and reactivity of 1a and 2. We have previously reported that the catalytic activity of 2 was improved by the introduction of electron-donating groups to the 4-position of the pyridine ring in PNP41. In this report, DFT calculations demonstrated that the introduction of electron-donating groups to PNP enhances the backdonating ability of molybdenum centres and thereby leads to activation of dinitrogen ligands. As described in the Introduction, the NHC-based PCP ligand was expected to work as a strong electron donor to activate dinitrogen ligands coordinated to the molybdenum centre. For understanding the geometric and electronic structures of 1a and 2 upon the coordination of the pincer ligands, mononuclear molybdenum complexes 1a0 and 20 were investigated in detail.

Figure 8a compares the electron-donating ability of Bim-PCP[1] and PNP in terms of differences in atomic charge (Dq)

between the dinitrogen complexes (1a0 and 20) and the free ligands (Bim-PCP[1] and PNP). The charges of the Mo atom and three N2 ligands obtained with the natural population analysis (NPA)62 were set to zero for the free ligands, and hence the sum of the charges is identical to the Dq value of the Mo(N2)3 moiety.

The gross NPA charges of the Mo(N2)3 moiety in 1a0 and 20 can be regarded as the amount of electron donated from the pincer ligands during complexation. To evaluate the electron-donating ability of Bim-PCP[1] and PNP, the Dq values of the PtBu2 groups and the carbene or pyridine moiety containing the methylene linkers are separately given in Fig. 8a. The Dq values of the PtBu2 groups are identical in both Bim-PCP[1] and PNP ( 0.29), indicating that the electron-donating ability of Bim-

PCP[1] and PNP is controlled by the carbene and pyridine moieties. Since the Dq values of the carbene and pyridine moieties are 0.23 and 0.12, respectively, the pincer ligands donate

0.81e (Bim-PNP[1]) and 0.70e (PNP) to the Mo(N2)3 moiety during complexation. As we expected, the NHC-based pincer ligand exhibits a stronger electron-donating ability than the pyridine-based one from a viewpoint of atomic charge.

Optimized structures and BDEs between the molybdenum centre and dinitrogen ligands also reect the strength of the

2 NH3

50

40

30

20

10

0

0 5 10 15 20 Time (h)

1

NH 3/equiv. cat.

Figure 6 | Time proles of the formation of ammonia from dinitrogen gas. A solution of CrCp*2 (0.72 mmol) in toluene was added to a mixture of 1a or 1c (0.0033 mmol) and [LutH]OTf (0.96 mmol) at room temperature over a period of 1 h under 1 atm of dinitrogen gas, followed by stirring for the indicated time (0.33, 0.67, 1, 2 and 20 h). A solution of CoCp2 (2.16 mmol)

in toluene was added to a mixture of 2 (0.010 mmol) and [LutH]OTf(2.88 mmol) at room temperature over a period of 1 h under 1 atm of dinitrogen gas, followed by stirring for the indicated time (0.33, 0.67, 1, 2 and 20 h). The amount of ammonia (equiv.) is based on the catalyst.

NATURE COMMUNICATIONS | 8:14874 | DOI: 10.1038/ncomms14874 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14874

e (216 equiv.) H+ (288 equiv.)

toluene, rt, 20 h

e (216 equiv.) H+ (288 equiv.)

toluene, rt, 2 h+ + +

cat. NH3 NH3

NH3 first NH3 second

cat. 1a

cat. 2

N2 N2 N2

47 equiv. 22 equiv.

0 equiv.

23 equiv.

Figure 7 | Reactions of further addition of proton source and reductant. NH3 rst and NH3 second were collected in separated runs. Each NH3 rst is the same value as the time prole experiment. Each NH3 second is collected by the following procedure. A solution of a reductant (CrCp*2 for 1a and CoCp2 for 2; 216 equiv) in toluene (4 ml) was added to a mixture of 1a or 2 (0.0033 mmol for 1a and 0.010 mmol for 2) and [LutH]OTf (288 equiv) at room temperature over a period of 1 h under 1 atm of dinitrogen gas, followed by stirring for 20 h. Then, [LutH]OTf (288 equiv) was added in one portion and another solution of the same reductant (216 equiv) in toluene (4 ml) was added over a period of 1 h, followed by stirring at room temperature for another 1 h under 1 atm of dinitrogen gas. The difference between the amount of ammonia obtained in this experiment and the NH3 rst is the NH3 second. The amount of ammonia (equiv.) is based on the catalyst.

electron-donating ability of Bim-PCP[1] and PNP. Table 3 summarizes geometric parameters around the molybdenum centre in mononuclear molybdenumdinitrogen complexes 1a0 and 20, together with the BDEs of the MoN2(axial) and

MoN2(equatorial) bonds. The Mo-C(carbene) bond distance(2.099 ) in 1a0 is signicantly shorter than the MoN(pyridine) distance (2.240 ) in 20. The Mayer bond order (b.o.)63 of the Mo-C(carbene) bond is calculated to be 0.91, which is much larger than that of the MoN(pyridine) bond (0.39). The presence of a strong bonding interaction between the molybdenum centre and Bim-PCP[1] is consistent with the experimental fact that 1a bearing Bim-PCP[1] works as a long-lived catalyst for the catalytic nitrogen xation with a high turnover number compared with 2 bearing PNP (vide supra). On the other hand, the strong trans inuence of the carbene ligand weakens the Mo N2(equatorial) bond in 1a0. The MoN2(equatorial) bond distance of 2.084 (b.o. 0.50) in 1a0 is much longer than

that of 2.018 (b.o. 0.62) in 20, and the BDEs of the

MoN2(equatorial) bond are 21.2 kcal mol 1 for 1a0 and30.1 kcal mol 1 for 20. Interestingly, the coordination of Bim-PCP[1] to molybdenum would inuence all dinitrogen ligands at both trans- and cis-positions so as to weaken all the MoN2 bonds. The MoN2(axial) bond distances (b.o.) are calculated to be 2.034 (0.53) for 1a0 and 2.024 (0.54) for 20. The BDE of the MoN2(axial) bond of 1a0 (12.5 kcal mol 1) is also lower than that of 20 (14.0 kcal mol 1). A similar trend was observed for the MoN2(bridging) and MoN2(terminal) bonds in dimolybdenumdinitrogen complexes 1a and 2. The bond dissociation energies are 18.8 and 11.9 kcal mol 1 for 1a, both of which are smaller than those obtained for 2 (24.9 and14.4 kcal mol 1).

The origin of the weaker MoN2 bonds in 1a0 and 1a is

understood by looking at frontier orbitals responsible for the bonding between the molybdenum centre and the carbene C atom of Bim-PCP[1]52,53,64,65. As depicted in Fig. 8b, the HOMO-6 (1a0) and HOMO-5 (20) contribute to a s-bond between the Mo atom and the carbene C atom (or the pyridine N atom). The large size of the lobe between the Mo and C atoms indicates that Bim-PCP[1] works as a strong s-donor compared to PNP. The HOMO-1 in Fig. 8c mainly contributes to p-backdonation from an out-of-plane d orbital of Mo to a p* orbital of dinitrogen ligands. The backdonation from metal to dinitrogen is essential for the activation of dinitrogen upon the

formation of metaldinitrogen complexes. Occupation of the HOMO-1 strengthens all of the MoN2 bonds in 1a0 and 20 because of their symmetrical structures. By comparing the HOMO-1 of 1a0 and 20, one can nd a bonding interaction between the Mo atom and the carbene C atom through p-backdonation from the d orbital of Mo to the vacant p orbital of C perpendicular to the carbene ring in 1a0. This backdonation decreases the amount of electron transferred to both the equatorial and axial dinitrogen ligands, leading to the lower BDEs of the MoN2 bonds in 1a0 (1a). As presented in

Fig. 8a, the Dq value of the dinitrogen ligands in 1a0 ( 0.26) is

smaller than that in 20 ( 0.31) in spite of the electron-donating

ability of Bim-PCP[1] superior to PNP. The backdonation from the Mo atom to the carbene C atom also contributes to the strong binding of Bim-PCP[1] to Mo.

On the other hand, the Mo-C bond distance (b.o.) in 1b0 bearing Im-PCP[2] (2.178 (0.79)) indicates that the Mo-C bond in 1b0 is weaker than that in 1a0, as summarized in Table 3.

Owing to the longer CH2 linkers, the coordination of the carbene moiety to the molybdenum centre in 1b0 is highly twisted compared to 1a0; the dihedral angles of N(1)-C(1)-Mo(1)-N(3) are 69.6 for 1a0 and 43.8 for 1b0 (Supplementary Fig. 21). The twisted coordination of the carbene moiety in 1b0 reduces the overlap between the out-of-plane d orbital of the Mo atom and the vacant p orbital of the carbene C atom. As a result, Im-PCP[2] works only as a very strong s-donor (0.90e donation to Mo;

Fig. 8a). The gross NPA charge on the dinitrogen ligands ( 0.36) as well as the large BDE of the MoN2(axial) bond of 1b0

(14.3 kcal mol 1) implies that the coordination of Im-PCP[2] to the molybdenum centre effectively activates the coordinated dinitrogen ligands. However, we theoretically conrmed that the mononuclear molybdenumdinitrogen complexes such as 1b0 and [cis-Mo(N2)2(Im-PCP[2])] cannot be protonated by

LutH , similar to 1a0 and [Mo(N2)3(PNP)]43. All attempts to optimize a product complex comprise the protonated 1a0 (1b0), and Lut resulted in formation of a reactant complex comprising 1a0 (1b0) and LutH , even though the optimization started from a structure with the N2yH distance of 5 . Thus, the lack of the catalytic activity of 1b for nitrogen xation can be attributed to the thermodynamic instability of the dinitrogen-bridged dimolybdenum structure, as mentioned in the former section.

On the basis of the catalytic reaction pathway previously proposed for nitrogen xation using 2 (ref. 43), we have

8 NATURE COMMUNICATIONS | 8:14874 | DOI: 10.1038/ncomms14874 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14874 ARTICLE

a

+0.29

+0.29

+0.29 +0.30

+0.30

H2 H2 H4

H4

C2

C2

+0.23 +0.12 +0.30

PtBu2

PtBu2 P

C C

C

0.26

(N2)3

N

0.31 0.36

(N2)3 (N2)3

0.55 0.39 0.54

Mo Mo Mo

1a

N

N NN

P

N

N

tBu2

PtBu2

N

N

C

H2

PtBu2 +0.29

PtBu2

H2

2 1b

b

N

N N

P

P

Mo N N

P = PtBu2

[Mo(N2)3(Bim-PCP[1]) (1a)

HOMO-6 (1a) HOMO-5 (2)

C Mo

N2 (axial)

N2 (axial)

N2 (equatorial)

P

Mo N N

N

N

N

N

P = PtBu2

[Mo(N2)3(PNP)] (2)

d

c

Mo

C

HOMO-1 (1a) HOMO-1 (2)

Figure 8 | Electronic properties of mononuclear molybdenumdinitrogen complexes. (a) Changes in the NPA atomic charge (Dq) in the coordination of the pincer ligands to the Mo(N2)3 moiety. The values of Dq are obtained as differences between mononuclear MoN2 complexes (1a0, 20 and 1b0) and free ligands (Bim-PCP[1] for 1a0, PNP for 20 and Im-PCP[2] for 1b0). (b) Spatial distribution of frontier orbitals of 1a0 and 20 that contribute to s donation from the pincer ligand to Mo. (c) Spatial distribution of frontier orbitals of 1a0 and 20 that contribute to p back donation from Mo to both equatorial and axial dinitrogen ligands. The molecular structures are rotated by 90 along the MoN2(equatorial) bond from those in Fig. 8b. (d) A schematic drawing of the bonding interactions between the Mo atom and the carbene C atom in 1a0.

theoretically investigated possible reaction pathways catalysed by 1a. In the present article, we particularly focus on the rst protonation process shown in Fig. 9, since the protonation of a terminal dinitrogen ligand in 2 by [LutH]OTf is energetically the most unfavourable process in the catalytic cycle43. In the calculated reaction pathway, a terminal dinitrogen ligand in 1a is rst protonated by LutH (1a-A-PCP), and then the dinitrogen ligand trans to the generated NNH group is eliminated (A-PCP-B-PCP). The protonation of 1a yielding A-PCP is endothermic by 8.1 kcal mol 1 with an activation energy of8.3 kcal mol 1. This energy prole indicates that proton detachment from A-PCP can easily occur like the PNP system.

On the other hand, the following N2 elimination yielding B-PCP is exothermic by 5.2 kcal mol 1 with a low activation energy of 4.0 kcal mol 1. The coordination of OTf to B-PCP is highly

exothermic by 20.7 kcal mol 1, and thus the whole reaction pathway leading to C-PCP is energetically downhill. Comparison of the energy proles of the PCP and PNP systems suggests that the reactivity of the dinitrogen complexes 1a and 2 with [LutH]OTf would not be a major factor for rationalizing the high catalytic activity of 1a.

DiscussionOn the basis of our previous ndings on the catalytic nitrogen xation, we have newly designed and prepared novel dinitrogen-bridged dimolybdenum complexes bearing NHC and phosphine-based PCP-pincer ligands, Bim-PCP[1] and Im-PCP[2]. The dimolybdenumdinitrogen complexes bearing Bim-PCP[1] as PCP-pincer ligands have been found to work as so far the most effective catalysts towards the ammonia formation from

NATURE COMMUNICATIONS | 8:14874 | DOI: 10.1038/ncomms14874 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14874

Table 3 | Selected bond distances in and MoN2 BDEs in kcal mol 1 of mononuclear molybdenumdinitrogen complexes.

1a0 20 1b0

Mo-C(carbene)/N(pyridine) 2.099 (0.91) 2.240 (0.39) 2.178 (0.79) MoN(equatorial) 2.084 (0.50) 2.018 (0.62) 2.041 (0.55) MoN(axial) 2.034 (0.53) 2.024 (0.54) 2.026 (0.53) BDE(equatorial N2) 21.2 30.1 21.5

BDE(axial N2) 12.5 14.0 14.3

N NP

PN P = PtBu2 P = PtBu2 P = PtBu2

N

[Mo(N2)3(Bim-PCP[1])] (1a') [Mo(N2)3(PNP)] (2') [Mo(N2)3(Im-PCP[2])] (1b')

NN Mo N N

NN P

P N N

N Mo N N

N N P

P N N

NN Mo N N

BDE, bond dissociation energy.

The Mayer bond orders are presented in parenthesis.

1a

or

2

H

+

H

+

H

N N

N N

N

N

+ H+

8.1 (8.3)

N N

N N

N2

5.2 (4.0) 20.7

15.5

1.9 (4.4)

N

N

+ OTf

N

N

N N

N

N

Mo

L3

L3 L3 L3 L3

Mo

N

N N

N

Mo Mo Mo

OTf

N N

6.3 (8.2)

L3 Mo

N N

A B C

N N

N N

N

N

N N

Mo Mo

N N

Mo Mo Mo Mo

N N

N N

N

N N N

N N

N N

N

N

L3 L3

N N

NN P

P P

P

N N

N N

P P

P

P

P = PtBu2 P = PtBu2

N N

1a: L3 = Bim-PCP[1]

2: L3 = PNP

NN[{Mo(N2)2(Bim-PCP[1])}2(-N2)] (1a) [{Mo(N2)2(PNP)}2(-N2)] (2)

Figure 9 | A possible reaction pathway and energy proles of the rst protonation process on a terminal dinitrogen ligand in 1a (highlighted in yellow) and 2. Energy changes (activation energies in parentheses) are presented in kcal mol 1.

molecular dinitrogen under ambient reaction conditions, where up to 230 equiv. of ammonia were produced based on the catalyst (115 equiv. of ammonia based on the molybdenum atom of the catalyst). The superior activity of dimolybdenumdinitrogen complexes bearing Bim-PCP[1] included the high TOF for ammonia formation and the catalyst stability. DFT calculations on 1 reveal that Bim-PCP[1] as a PCP-pincer ligand serves as not only a strong s-donor but also a p-acceptor. These electronic properties are responsible for a solid connection between the molybdenum centre and the pincer ligand, leading to the enhanced catalytic activity for nitrogen xation.

Methods

General information. Detailed experimental procedures, characterization of compounds and the computational details can be found in the Supplementary Figs 121, Supplementary Tables 110 and Supplementary Methods. Cartesian coordinates are available in Supplementary Data 1.

Synthesis of [MoCl3(PCP)] (3a3c). A typical procedure for the preparation of 3a is described below. To a mixed solid of 1,3-bis((di-tert-butylphosphino)methyl)-1H-benzo[d]imidazol-3-ium hexauorophosphate (7a, 1.16 g, 2.00 mmol) and KN(SiMe3)2 (559 mg, 2.80 mmol) was added toluene (40 ml), and the resulting suspension was stirred for 20 min at room temperature. [MoCl3(thf)3] (733 mg,1.75 mmol) and toluene (15 ml) were added to the suspension and stirred at 80 C for 19 h. The solvent was removed under vacuum, and the residue was washed with hexane (5 ml 2), toluene (10 ml) and hexane (5 ml 2). The solid was dried

under vacuum. The solid was extracted with CH2Cl2 (10 ml 1, 5 ml 7),

recrystallized from CH2Cl2hexane and dried under vacuum to afford3a 0.5CH2Cl2 (712 mg, 1.05 mmol, 60%). Anal. Calcd. for C25.5H45Cl4MoN2P2

(3a 0.5CH2Cl2): C, 45.08; H, 6.68; N, 4.12. Found: C, 45.41; H, 6.61; N, 4.44.

Crystals suitable for preliminary X-ray analysis were prepared by recrystallizing from CH2Cl2hexane to give 3a CH2Cl2 as orange-brown crystals.

3b: Recrystallization from CH2Cl2hexane gave 3b 0.5CH2Cl2 as orange

crystals. 33% yield. Anal. Calcd. for C23.5H47Cl4MoN2P2 (3b 0.5CH2Cl2): C, 42.94;

H, 7.21; N, 4.26. Found: C, 43.04; H, 7.44; N, 4.19. Crystals suitable for X-ray analysis were prepared by recrystallizing from 1,2-dichloroethanehexane to afford 3b 1/3C6H14. The structure is included in Supplementary Fig. 2 and selected bond

lengths and angles in 3b are included in Supplementary Table 6.3c: Recrystallization from CH2Cl2hexane afforded 3c 0.5CH2Cl2 as orange

crystals. 48% yield. Anal. Calcd. for C27.5H49Cl4MoN2P2 (3c 0.5CH2Cl2): C, 46.69;

H, 6.98; N, 3.96. Found: C, 46.58; H, 6.79; N, 4.07. The structure is included in Supplementary Fig. 3, and selected bond lengths and angles in 3c are included in Supplementary Table 7.

Synthesis of [{Mo(N2)2(PCP)}2(l-N2)] (1a1c). A typical procedure for the preparation of 1a is described below. To a suspension of NaHg (0.5 wt% Na, 13.8g,3.00 mmol) in THF was added [MoCl3(Bim-PCP[1])] 0.5CH2Cl2 (341 mg,

0.501 mmol), and the resulting suspension was stirred under atmospheric pressure of N2 at room temperature for 17 h. The supernatant suspension was ltered through

Celite, and the solvent was removed under vacuum. The resulting solid was extracted with benzene (5 ml) and ltered through Celite. The lter cake was washed with benzene (2 ml 9), and the solvent of the combined solution was removed under

vacuum. The resulting solid was washed with pentane (2 ml 20) to give

1a 1.3C4H8O 0.4C5H12 as a dark purple solid, where 1.3equiv. of THF and 0.4 equiv.

of hexane were determined by 1H NMR (141 mg, 0.102 mmol, 37%). Analytically pure sample was prepared by recrystallization from THF at 18 C to afford 1a.

1a: 1H NMR (C6D6): d 7.016.98 (m, 4H, ArH), 6.756.71 (m, 4H, ArH), 3.84 (s, 8H, NCH2P), 1.38 (pseudo t, 3JPH 5.7 Hz, 72H, PtBu2). 31P{1H} NMR (C6D6):

105.6 (s, PtBu2). Infrared (KBr, cm 1): 1,978 (s, nNN). Infrared (THF, cm 1): 1,979 (s, nNN). Anal. Calcd. for C50H88Mo2N14P4: C, 50.00; H, 7.38;

N, 16.33. Found: C, 49.66; H, 6.94; N, 14.18. The lower content of nitrogen is

10 NATURE COMMUNICATIONS | 8:14874 | DOI: 10.1038/ncomms14874 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14874 ARTICLE

considered to be due to the labile property of the coordinated dinitrogen ligand in 1a under the analytical conditions.

1b: Recrystallization from benzenehexane afforded 1b 2/3C6H14 as dark-brown

crystals. 53% yield. Anal. Calcd. for C50H101.33Mo2N14P4 (1b 2/3C6H14): C, 49.44; H,

8.41; N, 16.14. Found: C, 49.82; H, 8.40; N, 15.35. 1H and 31P{1H} NMR were measured in THF-d8 as a mixture of 1b and 1b0. 1H NMR (THF-d8): 1b, d 6.76 (s, 2H,

NCHCHN), 1.44 (pseudo t, 3JPH 11.2 Hz, PtBu2), 0.95 (pseudo t, 3JPH 11.2 Hz,

PtBu2): 1b0, d 6.85 (s, 2H, NCHCHN), 4.294.19 (br m, 4H, NCH2), 1.981.91 (br m, 4H, CH2P), 1.18 (br s, 36H, PtBu2). 31P{1H} NMR (THF-d8): 1b, d 71.1 (s, tBu2P): 1b0, d 69.2 (s, tBu2P). Infrared (KBr, cm 1): 1,911 (s, nNN for 1b). Infrared (THF under N2, cm 1): 2,041 (m, nNN for 1b0), 1,945 (s, nNN for 1b0).

1c: Reprecipitation from THFhexane afforded 1c 0.5C6H14 as a brown solid.

46% yield. Crystals suitable for X-ray analysis were prepared by recrystallization from THF at 18 C to afford 1c. The structure is included in Supplementary

Fig. 1, and selected bond lengths and angles are included in Supplementary Table 5.

1H NMR (C6D6): d 6.68 (s, 4H, ArH), 3.89 (s, 8H, NCH2P), 2.26 (s, 12H, ArCH3),1.41 (pseudo t, 3JPH 5.5 Hz, 72H, PtBu2). 31P{1H} NMR (C6D6): d 105.8

(s, PtBu2). Infrared (KBr, cm 1): 1,969 (s, nNN). Infrared (THF, cm 1): 1,973 (s, nNN). Anal. Calcd. for C57H103Mo2N14P4 (1c 0.5C6H14): C, 52.65; H, 7.98;

N, 15.08. Found: C, 52.80; H, 7.70; N, 13.37.

Catalytic reduction of dinitrogen to ammonia under N2. In a 50 ml Schlenk ask were placed 1a 1.3C4H8O 0.4C5H12 (12.9 mg, 0.00970 mmol) and 2,6-lutidinium

triuoromethanesulfonate [LutH]OTf (247 mg, 0.960 mmol). Toluene (1.0 ml) was added under N2 (1 atm), and then a solution of CrCp*2 (232 mg, 0.719 mmol) in toluene (4.0 ml) was added to the stirred suspension in the Schlenk ask with a syringe pump at a rate of 4.0 ml h 1. After the addition of CrCp*2, the mixture was further stirred at room temperature for 19 h. The reaction mixture was evaporated under reduced pressure, and the distillate was trapped in dilute H2SO4 solution(0.5 M, 10.00 ml). Aqueous solution of potassium hydroxide (30 wt%; 5 ml) was added to the residue to fully liberate ammonia, and the mixture was distilled into another dilute H2SO4 solution (0.5 M, 10.00 ml). The amount of NH3 present in each of the H2SO4 solutions was determined by the indophenol method.61

Computational method. DFT calculations were performed to search all intermediates and transition structures on potential energy surfaces using the Gaussian 09 programme66. Similar to the previous study43, we adopted the B3LYP* functional, which is a reparametrized version of the B3LYP hybrid functional developed by Reiher et al.67. For optimization, the StuttgartDresden pseudopotentials and 631G(d) basis sets were chosen for the Mo atom and the other atoms, respectively. To determine the energy prole of the rst protonation process, we performed single-point energy calculations at the optimized geometries using the 6311 G(d,p) basis sets in place of the 631G(d) basis sets. Zero-point

energy corrections were applied for energy changes (DE) and activation energies (Ea) calculated for each reaction step. Solvation effects (toluene) were taken into account by using the polarizable continuum model in the single-point energy calculations68. More details are described in Supplementary Methods. Throughout the paper, the BDE of an MoN2 (terminal, axial or equatorial) bond is dened as the energy change for dissociation of the corresponding dative N2 ligand, for example, [{Mo(N2)2(PNP)}2(m-N2)]-[{Mo(N2)(PNP)}-NN-{Mo(N2)2(PNP)}] N2.

The BDE of the MoN2 (bridging) bond is dened as the energy change for separation of a dimolybdenum complex into two mononuclear MoN2 complexes, such as [{Mo(N2)2(PNP)}2(m-N2)]-[Mo(N2)3(PNP)] cis-[Mo(N2)2(PNP)].

Data availability. The X-ray crystallographic coordinates for structures reported in this article have been deposited at the Cambridge Crystallographic Data Centre (CCDC), under deposition number CCDC 1482254 (1a), 1482255 (1b), 1482256 (1c), 1482258 (3b) and 1482259 (3c). These data can be obtained free of charge from The Cambridge Crystallographic Data Centre via http://www.ccdc.cam.ac.uk/data_request/cif

Web End =www.ccdc.cam.ac.uk/ http://www.ccdc.cam.ac.uk/data_request/cif

Web End =data_request/cif . All other data are available from the authors upon reasonable request.

References

1. Liu, H. Ammonia Synthesis Catalysts: Innovation and Practice (World Scientic, 2013).

2. Eady, R. R. Current status of structure function relationships of vanadium nitrogenase. Coord. Chem. Rev. 237, 2330 (2003).

3. Spatzal, T. et al. Evidence for interstitial carbon in nitrogenase FeMo cofactor. Science 334, 940 (2011).

4. Lancaster, K. M. et al. X-ray emission spectroscopy evidences a central carbon in the nitrogenase iron-molybdenum cofactor. Science 334, 974977 (2011).

5. Wiig, J. A., Hu, Y., Lee, C. C. & Ribbe, M. W. Radical SAM-dependent carbon insertion into the nitrogenase M-cluster. Science 337, 16721675 (2012).

6. Hoffman, B. M., Lukoyanov, D., Yang, Z.-Y., Dean, D. R. & Seefeldt, L. C. Mechanism of nitrogen xation by nitrogenase: the next stage. Chem. Rev. 114, 40414062 (2014).

7. Hidai, M. & Mizobe, Y. Recent advances in the chemistry of dinitrogen complexes. Chem. Rev. 95, 11151133 (1995).

8. MacKay, B. A. & Fryzuk, M. D. Dinitrogen coordination chemistry: on the biomimetic borderlands. Chem. Rev. 104, 385402 (2004).

9. Khoenkhoen, N., de Bruin, B., Reek, J. N. H. & Dzik, W. I. Reactivity of dinitrogen bound to mid- and late-transition-metal centers. Eur. J. Inorg. Chem. 2015, 567598 (2015).

10. Nishibayashi, Y. Recent progress in transition-metal-catalyzed reduction of molecular dinitrogen under ambient reaction conditions. Inorg. Chem. 54, 92349247 (2015).

11. Tanaka, H., Nishibayashi, Y. & Yoshizawa, K. Interplay between theory and experiment for ammonia synthesis catalyzed by transition metal complexes. Acc. Chem. Res. 49, 987995 (2016).

12.ori, I., Mercado, B. Q., Bill, E., Vinyard, D. J. & Holland, P. L. Binding of dinitrogen to an iron-sulfur-carbon site. Nature 526, 9699 (2015).

13. McSkimming, A. & Harman, W. H. A terminal N2 complex of high-spin iron(I) in a weak, trigonal ligand eld. J. Am. Chem. Soc. 137, 89408943 (2015).

14. Lee, Y. et al. Dinitrogen activation upon reduction of a triiron(II) complex. Angew. Chem. Int. Ed. 54, 14991503 (2015).

15. Keane, A. J., Farrell, W. S., Yonke, B. L., Zavalij, P. Y. & Sita, L. R. Metal-mediated production of isocyanates, R3EN C O from dinitrogen, carbon

dioxide, and R3ECl. Angew. Chem. Int. Ed. 54, 1022010224 (2015).16. Klopsch, I., Kinauer, M., Finger, M., Wrtele, C. & Schneider, S. Conversion of dinitrogen into acetonitrile under ambient conditions. Angew. Chem. Int. Ed. 55, 47864789 (2016).

17. Bazhenova, T. A. & Shilov, A. E. Nitrogen xation in solution. Coord. Chem. Rev. 144, 69145 (1995).

18. Shilov, A. E. Catalytic reduction of molecular nitrogen in solutions. Russ. Chem. Bull. 52, 25552562 (2003).

19. Shiina, K. Reductive silylation of molecular nitrogen via xation to tris(trialkylsilyl)amine. J. Am. Chem. Soc. 94, 92669267 (1972).

20. Komori, K., Oshita, H., Mizobe, Y. & Hidai, M. Preparation and properties of molybdenum and tungsten dinitrogen complexes. 25. Catalytic conversion of molecular nitrogen into silylamines using molybdenum and tungsten dinitrogen complexes. J. Am. Chem. Soc. 111, 19391940 (1989).

21. Mori, M. Activation of nitrogen for organic synthesis. J. Organomet. Chem. 689, 42104227 (2004).

22. Tanaka, H. et al. Molybdenum-catalyzed transformation of molecular dinitrogen into silylamine: experimental and DFT study on the remarkable role of ferrocenyldiphosphine ligands. J. Am. Chem. Soc. 133, 34983506 (2011).

23. Liao, Q., Saffon-Merceron, N. & Mzailles, N. Catalytic dinitrogen reduction at the molybdenum center promoted by a bulky tetradentate phosphine ligand. Angew. Chem. Int. Ed. 53, 1420614210 (2014).

24. Siedschlag, R. B. et al. Catalytic silylation of dinitrogen with a dicobalt complex.J. Am. Chem. Soc. 137, 46384641 (2015).25. Imayoshi, R. et al. Cobalt-catalyzed transformation of molecular dinitrogen into silylamine under ambient reaction conditions. Chem. Eur. J. 21, 89058909 (2015).

26. Liao, Q., Saffon-Merceron, N. & Mzailles, N. N2 reduction into silylamine at tridentate phosphine/Mo center: catalysis and mechanistic study. ACS Catal 5, 69026906 (2015).

27. Kuriyama, S. et al. Catalytic transformation of dinitrogen into ammonia and hydrazine by iron-dinitrogen complexes bearing pincer ligand. Nat. Commun. 7, 12181 (2016).

28. Kuriyama, S. et al. Direct transformation of molecular dinitrogen into ammonia catalyzed by cobalt dinitrogen complexes bearing anionic PNP pincer ligands. Angew. Chem. Int. Ed. 55, 1429114295 (2016).

29. Yandulov, D. V. & Schrock, R. R. Catalytic reduction of dinitrogen to ammonia at a single molybdenum center. Science 301, 7678 (2003).

30. Ritleng, V. et al. Molybdenum triamidoamine complexes that contain hexa-tert-butylterphenyl, hexamethylterphenyl, or p-bromohexaisopropylterphenyl substituents. An examination of some catalyst variations for the catalytic reduction of dinitrogen. J. Am. Chem. Soc. 126, 61506163 (2004).

31. Schrock, R. R. Catalytic reduction of dinitrogen to ammonia at a single molybdenum center. Acc. Chem. Res. 38, 955962 (2005).

32. Schrock, R. R. Catalytic reduction of dinitrogen to ammonia by molybdenum: theory versus experiment. Angew. Chem. Int. Ed. 47, 55125522 (2008).33. Anderson, J. S., Rittle, J. & Peters, J. C. Catalytic conversion of nitrogen to ammonia by an iron model complex. Nature 501, 8487 (2013).

34. Creutz, S. E. & Peters, J. C. Catalytic reduction of N2 to NH3 by an Fe-N2 complex featuring a C-atom anchor. J. Am. Chem. Soc. 136, 11051115 (2014).

35. Ung, G. & Peters, J. C. Low-temperature N2 binding to two-coordinate L2Fe0 enables reductive trapping of L2FeN2 and NH3 generation. Angew. Chem. Int.

Ed. 54, 532535 (2015).36. Rittle, J. & Peters, J. C. Fe-N2/CO complexes that model a possible role for the interstitial C atom of FeMo-cofactor (FeMoco). Proc. Natl Acad. Sci. USA 110, 1589815903 (2013).

37. Rittle, J., McCrory, C. C. L. & Peters, J. C. A 106-fold enhancement in N2-binding afnity of an Fe2(m-H)2 core upon reduction to a mixed-valence FeIIFeI state. J. Am. Chem. Soc. 136, 1385313862 (2014).

NATURE COMMUNICATIONS | 8:14874 | DOI: 10.1038/ncomms14874 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 11

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms14874

38. Rittle, J. & Peters, J. C. An Fe-N2 complex that generates hydrazine and ammonia via Fe NNH2: demonstrating a hybrid distal-to-alternating pathway

for N2 reduction. J. Am. Chem. Soc. 138, 42434248 (2016).39. Del Castillo, T. J., Thompson, N. B. & Peters, J. C. A synthetic single-site Fe nitrogenase: high turnover, freeze-quench 57Fe Mssbauer data, and a hydride resting state. J. Am. Chem. Soc. 138, 53415350 (2016).

40. Arashiba, K., Miyake, Y. & Nishibayashi, Y. A molybdenum complex bearing PNP-type pincer ligands leads to the catalytic reduction of dinitrogen into ammonia. Nat. Chem. 3, 120125 (2011).

41. Kuriyama, S. et al. Catalytic formation of ammonia from molecular dinitrogen by use of dinitrogen-bridged dimolybdenum-dinitrogen complexes bearing PNP-pincer ligands: remarkable effect of substituent at PNP-pincer ligand. J. Am. Chem. Soc. 136, 97199731 (2014).

42. Kuriyama, S. et al. Nitrogen xation catalyzed by ferrocene-substituted dinitrogen-bridged dimolybdenum-dinitrogen complexes: unique behavior of ferrocene moiety as redox active site. Chem. Sci. 6, 39403951 (2015).

43. Tanaka, H. et al. Unique behaviour of dinitrogen-bridged dimolybdenum complexes bearing pincer ligand towards catalytic formation of ammonia. Nat. Commun. 5, 3737 (2014).

44. Arashiba, K. et al. Catalytic reduction of dinitrogen to ammonia by use of molybdenum-nitride complexes bearing a tridentate triphosphine as catalysts.J. Am. Chem. Soc. 137, 56665669 (2015).45. Hopkinson, M. N., Richter, C., Schedler, M. & Glorius, F. An overview of N-heterocyclic carbenes. Nature 510, 485496 (2014).

46. Nolan, S. P. (ed.) N-Heterocyclic Carbenes (Wiley-VCH, 2014).47. Trnka, T. M. & Grubbs, R. H. The development of L2X2Ru CHR olen

metathesis catalysts: an organometallic success story. Acc. Chem. Res. 34, 1829 (2001).48. Ohki, Y. & Seino, H. N-heterocyclic carbenes as supporting ligands in transition metal complexes of N2. Dalton Trans. 45, 874880 (2016).

49. Plikhta, A., Pthig, A., Herdtweck, E. & Rieger, B. Toward new organometallic architectures: synthesis of carbene-centered rhodium and palladium bisphosphine complexes. Stability and reactivity of [PCBImPRh(L)][PF6]

pincers. Inorg. Chem. 54, 95179528 (2015).50. Lee, H. M., Zeng, J. Y., Hu, C.-H. & Lee, M.-T. A new tridentate pincer phosphine/N-heterocyclic carbene ligand: palladium complexes, their structures, and catalytic activities. Inorg. Chem. 43, 68226829 (2004).

51. Chiu, P. L. & Lee, H. M. Chemistry of the PCNHCP ligand: silver and ruthenium complexes, facial/meridional coordination, and catalytic transfer hydrogenation. Organometallics 24, 16921702 (2005).

52. Gradert, C., Krahmer, J., Snnichsen, F. D., Nather, C. & Tuczek, F. Small-molecule activation with molybdenum(0) complexes supported by mixed inidazol-2-ylidene/phosphanyl hybrid ligandselectronic and structural consequences of substituting a phosphane by a carbene group. Eur. J. Inorg. Chem. 2013, 39433955 (2013).

53. Gradert, C., Krahmer, J., Snnichsen, F. D., Nather, C. & Tuczek, F. Molybdenum(0)-carbonyl complexes supported by mixed benzimidazol-2-ylidene/phosphine ligands: inuence of benzannulation on the donor properties of the NHC groups. J. Organomet. Chem. 770, 6168 (2014).

54. Bestgen, S., Gamer, M. T., Lebedkin, S., Kappes, M. M. & Roesky, P. W. Di- and trinuclear gold complexes of diphenylphosphinoethyl-functionalised imidazolium salts and their N-heterocyclic carbenes: synthesis and photophysical properties. Chem. Eur. J. 21, 601614 (2015).

55. Gradert, C., Stucke, N., Krahmer, J., Nather, C. & Tuczek, F. Molybdenum complexes supported by mixed NHC/phosphine ligands: activation of N2 and reaction with P(OMe)3 to the rst meta-phosphite complex. Chem. Eur. J. 21, 11301137 (2015).

56. Anderson, S. N. & Richards, R. L. Dinitrogen complexes of tungsten with bulky phosphine co-ligands. Preparation and crystal structures of [{W(N2)2(PEt2Ph)3}2(m-N2)], trans-[W(N2)2(PEt2Ph)4] C4H8O, and

[W(Z6-C6H5PPrn2)(N2)(PPrn2Ph)2] and their reactions to give hydrazine or ammonia. J. Chem. Soc. Dalton Trans. 245252 (1986).57. Anderson, S. N., Hughes, D. L. & Richards, R. L. Preparation of mer-[Mo(N2)3(PPrn2Ph)3] from trans-[Mo(N2)2(PPrn2Ph)4]; its X-ray crystal structure, spectroscopic properties, and protonation to give ammonia. J. Chem. Soc. Dalton Trans. 15911594 (1986).

58. Ohki, Y., Aoyagi, K. & Seino, H. Synthesis and protonation of N-heterocycliccarbene-supported dinitrogen complexes of molybdenum(0). Organometallics 34, 34143420 (2015).

59. Gusev, D. G. Donor properties of a series of two-electron ligands. Organometallics 28, 763770 (2009).

60. Gusev, D. G. Electronic and steric parameters of 76 N-heterocyclic carbenes in Ni(CO)3(NHC). Organometallics 28, 64586461 (2009).

61. Weatherburn, M. W. Phenol-hypochlorite reaction for determination of ammonia. Anal. Chem. 39, 971974 (1967).

62. Glendening, E. D. et al. NBO 5.9 (Theoretical Chemistry Institute, University of Wisconsin, 2009).

63. Mayer, I. Charge, bond order and valence in the AB initio SCF theory. Chem. Phys. Lett. 97, 270274 (1983).

64. Comas-Vives, A. & Harvey, J. N. How important is backbonding in metal complexes containing N-heterocyclic carbenes? Structural and NBO analysis. Eur. J. Inorg. Chem. 2011, 50255035 (2011).

65. Nelson, D. J. & Nolan, S. P. Quantifying and understanding the electronic properties of N-heterocyclic carbenes. Chem. Soc. Rev. 42, 67236753 (2013).

66. Frisch, M. J. et al. Gaussian 09, Revision C.01 (Gaussian, Inc., 2011).67. Reiher, M., Salomon, O. & Hess, B. A. Reparametrization of hybrid functionals based on energy differences of states of different multiplicity. Theor. Chem. Acc. 107, 4855 (2001).

68. Tomasi, J., Mennucci, B. & Cammi, R. Quantum mechanical continuum solvation models. Chem. Rev. 105, 29993094 (2005).

69. Castellani, M. P., Geib, S. J., Rheingold, A. L. & Trogler, W. C. Synthesis, reactivities, molecular structures, and physical properties of paramagnetic bis(tetraphenylcyclopentadienyl) complexes of vanadium, chromium, cobalt, and nickel. Organometallics 6, 17031712 (1987).

70. Clarke, K. & Rothwell, K. A kinetic study of the effect of substituents on the rate of formation of alkylpyridinium halides in nitromethane solution. J. Chem. Soc. 1960, 18851895 (1960).

Acknowledgements

The present project is supported by CREST, JST. We thank KAKENHI (Nos JP26288044, JP26105708, JP15K13687, JP15H05798 to Y.N., No. JP24109014 to K.Y. and No. JP26888008 to H.T.) from JSPS and MEXT. A.E. and S.K. are recipients of the JSPS Predoctoral Fellowships for Young Scientists. We also thank the Research Hub for Advanced Nano Characterization at The University of Tokyo for X-ray analysis.

Author contributions

K.Y. and Y.N. directed and conceived this project. A.E., K.A., S.K. and K.N. conducted the experimental work. H.T. and Y.M. conducted the computational work. All authors discussed the results and wrote the manuscript.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Eizawa, A. et al. Remarkable catalytic activity of dinitrogen-bridged dimolybdenum complexes bearing NHC-based PCP-pincer ligands toward nitrogen xation. Nat. Commun. 8, 14874 doi: 10.1038/ncomms14874 (2017).

Publishers note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional afliations.

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

r The Author(s) 2017

12 NATURE COMMUNICATIONS | 8:14874 | DOI: 10.1038/ncomms14874 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications