Josh
Abbenseth
a,
Markus
Finger
a,
Christian
Würtele
a,
Müge
Kasanmascheff
b and
Sven
Schneider
*a
aInstitut für Anorganische Chemie, Georg-August-Universität, Tammannstraße 4, 37077 Göttingen, Germany. E-mail: sven.schneider@chemie.uni-goettingen.de
bMax Planck Institute for Biophysical Chemistry, 37077 Göttingen, Germany
First published on 15th January 2016
The oxidative coupling of nitride ligands (N3−) to dinitrogen and its microscopic reverse, N2-splitting to nitrides, are important elementary steps in chemical transformations, such as selective ammonia oxidation or nitrogen fixation. Here an experimental and computational evaluation is provided for the homo- and heterocoupling of our previously reported iridium(IV) and iridium(V) nitrides [IrN(PNP)]n (n = 0, +1; PNP = N(CHCHPtBu2)2). All three formal coupling products [(PNP)IrN2Ir(PNP)]n (n = 0–+2) were structurally characterized. While the three coupling reactions are all thermodynamically feasible, homocoupling of [IrN(PNP)]+ is kinetically hindered. The contributing parameters to relative coupling rates are discussed providing qualitative guidelines for the stability of electron rich transition metal nitrides.
Sven Schneider |
Fig. 1 Qualitative molecular orbitals relevant for N2 splitting with Cummins’ complex (a) and schematic transition states that were proposed for N2-splitting/coupling (b–d). |
A similar picture arises for the reverse reaction, i.e. the coupling of terminal nitrides. Several authors reported the formation of dinitrogen bridged complexes upon reductive coupling of octahedral MV nitrides (M = Fe, Ru, Os).6 These transient nitrides are typically formed in situ, e.g. by photolysis of MIII-azides, reduction of stable MVI nitrides or oxidation of low-valent ammine complexes to give divalent ({MIIN2} or {MIIN2MII}) or mixed valent {MIIN2MIII} dinitrogen complexes. Except for a few cases, it is not fully clear whether the reaction proceeds via MVN homocoupling or coupling of MVN with parent MVIN. Importantly, within the simplified MO-picture (Fig. 1a) the MII/II and MII/III dinitrogen bridged complexes represent 12 and 11 π-electron {MNNM}-cores, respectively. However, some closed-shell nitrides also decay by nitride coupling underlining that the electronic configuration is not the only predictor for the thermochemistry of N2-splitting/coupling.7–9
In comparison, the relationship of electronic structure and kinetics received less attention. Ware and Taube examined the thermal decomposition of [OsVINCl3(py)2] by nitride coupling to [OsIIICl3(py)3].7 They pointed out that donation from a nitride lone-pair into an empty MN π*-orbital should be favored over a collinear approach of the two monomers and therefore proposed a semi-bent transition state (Fig. 1b). In line with such a polar transition state (TS), Seymore and Brown showed that heterocoupling of electrophilic OsVI- and nucleophilic MoVI-nitrides is considerably faster than the respective homocoupling reactions.10 However, Burger and coworkers emphasized that nitride homocoupling should favor a symmetrical “non least-motion pathway” or in-plane zig-zag TS, which features mutual σ/π donor–acceptor interactions of the two monomer fragments (Fig. 1c).11
Recently, we reported that photolysis of the iridium(II) azide [Ir(N3)(PNP)] (1, PNP = N(CHCHPtBu2)2) results in the formation of iridium(IV) nitride [IrN(PNP)] (2, Scheme 1).12 Spectroscopic and computational analysis suggested that 2 and the analogous rhodium compound are π-radicals with strong delocalization of the unpaired electron over the MN core, i.e. considerable “nitridyl” character.12,13 In solution, these transient nitrides undergo clean coupling, e.g. to 12-π-electron N2-complex [N2{Ir(PNP)}2] (3) in case of 2 (Scheme 1). The computed TS for a PMe2-truncated model featured a mutual π/π MO-interaction as the most favorable pathway with an out of plane zig-zag arrangement of the two monomer radicals (Fig. 1d).12 Alternatively, oxidation of 1 results in formation of isolable iridium(V) nitride [IrN(PNP)]+ (4), which in turn also gives 2 upon reduction (Scheme 1). In this context, we were interested whether the enhanced thermal stability of 4 can be attributed to thermodynamic or kinetic arguments. Here, we present a systematic experimental and computational study that addresses the stability of the iridium nitrides with respect to oxidative coupling.
Scheme 1 Syntheses of iridium(IV) and iridium(V) PNP nitrides 2 and 4 and coupling to iridium(I)/iridium(I) N2-complex 4.12 |
Chemical oxidation of 3 with one equivalent of [FeCp2]PF6 in dichloromethane at −35 °C affords the isolation of IrI/II N2-complex [N2{Ir(PNP)}2]PF6 (5) as a temperature, light, and air sensitive dark green solid (Scheme 2). The X-band EPR spectrum of 5 in frozen solution (4 K) reveals a rhombic signal corresponding to an S = 1/2 system (ESI†). The large anisotropy of the g-tensor (gx = 2.66, gy = 2.12, gz = 1.62) without resolved hyperfine coupling was similarly observed for several IrII(PNP) complexes with square-planar geometry.5,14 From the CV data a comproportionation constant KC is estimated to be around 107. A large KC like this is frequently stressed as a simple predictor for electronic delocalization in the mixed valent complex indicating stabilization due to resonance exchange.15 Note that KC data should be treated with caution and is easily over-interpreted due to the many other parameters possibly contributing to electrochemical potentials.16 In fact, the bond parameters in the solid state indicate valence localization in the crystal (see below) and the DFT model of 5 also exhibits localization of the spin density on one of the two Ir(PNP) moieties (ESI†). In contrast, the 1H NMR spectrum of 5 exclusively shows one set of paramagnetically shifted and broadened pincer ligand peaks (ESI†), suggesting rapid electron transfer between the N2-bridged moieties on the NMR timescale (1–10−4 s). Unfortunately, detailed analysis of the NIR spectrum was hampered by the thermal instability of 5 and dication 6 (see below) and did not allow for a reliable assignment of intervalence charge transfer vs. interconfigurational electronic transitions. However, IR spectroscopy reveals the presence of a relatively intense band assignable to the N2 stretching vibration both in solution (CH2Cl2: 1959 cm−1) and in the solid state (nujol: 1960 cm−1). The assignment was confirmed by DFT computations (ESI†) scaled to monomeric complex [Ir(N2)(PNP)] (2077 cm−1):12 The experimental shift ΔνN2 = 117 cm−1 is well reproduced by DFT (ΔνN2 = 115 cm−1). This comparison also excludes the formation of monomeric [Ir(N2)(PNP)]+, which should be blue shifted with respect to [Ir(N2)(PNP)]. Importantly, the observation of an N2-stretch confirms the presence of a permanent dipole, hence valence localization on the IR-timescale (∼10−14 s) also in solution providing an upper limit for the charge transfer rate. Similar properties between localized and delocalized mixed-valence (Class II–III) were found for N2-bridged OsII/III complexes.17 Unequivocal assignment of complex 5 requires further investigations.
Confirmation of the molecular structure of 5 is provided by single-crystal X-ray diffraction (Fig. 2 and Table 1). Both Ir ions are in distorted square-planar coordination geometries and bridged by a linear, end-on bound N2, as in parent 3. The structure features two crystallographically distinctly different Ir(PNP) fragments. For example, the Ir–NPNP and Ir–NN2 distances differ by 0.06 Å and 0.11 Å, respectively. Ir2 exhibits a longer bond to the PNP π-donor ligand and a shorter bond to the π-acceptor N2 than Ir1. In that respect, Ir2 behaves similar to Ir in 3. Hence, the structural parameters suggest trapped oxidation states in the solid state with Ir1 and Ir2 being associated with the +2 and +1 valence states, respectively. Importantly, the N–N distance (1.136(6) Å) is identical with that in 3 (1.135(4) Å) and in [N2{Ir(PCP)}2] (1.134(2) Å; PCP = C6H3-2,6-(CH2PtBu2)2)18 and close to free N2 (1.097 Å) revealing a weak degree of N2-activation upon oxidation.
Fig. 2 Molecular structure of 5 in the solid state (ellipsoids at 50% probability level, hydrogen atoms and solvent molecules omitted for clarity). |
3 | 5 | 6 | ||||
---|---|---|---|---|---|---|
Exp | DFT | Exp | DFT | Exp | DFT | |
a Computed spin ground–states: 13, 25, 36. | ||||||
Bond lengths (Å) | ||||||
Ir–N2 | 1.937(3)/1.933(3) | 1.922 | 1.984(5)/1.878(5) | 1.988/1.855 | 1.954(3) | 1.928 |
Ir–NPNP | 2.041(3)/2.035(3) | 2.055 | 1.973(4)/2.031(4) | 1.969/2.059 | 1.983(3) | 1.995 |
N–N | 1.135(4) | 1.136 | 1.136(6) | 1.141 | 1.138(6) | 1.140 |
Bond angles (°) | ||||||
Ir–N–N | 174.6(2)/172.3(3) | 172.8 | 175.3(4)/174.6(4) | 172.0/173.2 | 172.2(3) | 173.2 |
N–Ir–N | 173.6(2)/171.0(1) | 172.0 | 172.8(2)/173.3(2) | 169.9/171.6 | 168.0(1) | 169.0 |
The reaction of 3 with two equivalents of AgSbF6 at −70 °C in dichloromethane under the exclusion of light results in the immediate precipitation of metallic silver and formation of a deep red solution (Scheme 2). Monitoring by 1H NMR spectroscopy reveals the formation of a new paramagnetic compound with one set of signals assignable to a PNP pincer ligand, yet not as clean as in case of 5 (ESI†). No signals were found by X-band EPR spectroscopy as expected for a non-Kramers system. Decomposition in solution above −50 °C into several unidentified paramagnetic compounds is indicated by new 1H NMR signals upon warming. Importantly, the nitride 4 was not found as decomposition product.
The high lability of the product from double oxidation with respect to ambient temperatures and light hampered a more detailed spectroscopic characterization, so far. However, crystallization at low temperatures afforded red crystals of the dinitrogen bridged IrII/II complex [(N2){Ir(PNP)}2](SbF6)2 (6) that were suitable for X-ray diffraction. In contrast to mixed-valent complex 5, the molecular structure of 6 (Fig. 3 and Table 1) features two crystallographically equivalent Ir(PNP) fragments. The Ir–NPNP (1.983(3) Å) and Ir–NN2 (1.954(3) Å) bond lengths are very close to those of Ir1 in 5, supporting the valence state assignments for the structure of 5 and full oxidation of 6 to IrII/II. The N2 bond length (1.138(6) Å) is identical within error compared to the parent IrI/I (1.135(4) Å) and IrI/II (1.136(6) Å) complexes, indicating that N2 activation is negligible within the [N2{Ir(PNP)}2]n+ (n = 0–2) redox series.
Fig. 3 Molecular structure of 6 in the solid state (ellipsoids at 50% probability level, hydrogen atoms and solvent molecules omitted for clarity). |
Coupling of the iridium(V) and iridium(IV) nitrides to mixed-valent 5 was also examined (Scheme 3, B). Spectro-electrochemical reduction of 4 (Ered = −0.86 V) at room temperature results in clean conversion to IrI/I complex 3 (Eox = −0.27 V). The observation of an isosbestic point at 560 nm indicates the absence of long-lived intermediates on the experimental time-scale (ESI†). Moreover, the equimolar addition of 2 to 4 at −70 °C, where both pure nitrides are kinetically stable, results in immediate coupling to 5 (ESI†) confirming much more rapid coupling to 5 than to 3 or 6. Even the addition of a substoichiometric amount of 4 (10 mol%) to 2 at −70 °C results in the immediate formation with dimer 3 as the main product (ESI†). This observation can be rationalized with a redox-catalytic cycle (Scheme 4): Nitrides 2 and 4 undergo rapid coupling. The resulting dimer 5 is then reduced by starting material 2 to restore catalyst 4, which is in line with the potentials of redox couples 2/4 and 3/5, respectively. This mechanism also explains a previously reported observation for the two routes to iridium(IV) nitride 2 (Scheme 1): Unlike the photolysis route, the reduction route also yields considerable amounts of coupling product 3, even at low temperatures where coupling of pure 2 is not observed.
Computational evaluation of reactions A–C (Scheme 3) fully reproduces our experimental results (Table 2 and ESI†). Nitride coupling is thermodynamically strongly favored for all three reactions. However, the driving force within the redox series 3 (ΔG298 = −101.5 kcal mol−1), 5 (ΔG298 = −97.1 kcal mol−1), and 6 (ΔG298 = −73.9 kcal mol−1) decreases, as expected from the destabilizing population of an Ir–N π*-antibonding orbital in 2. For reaction A, a kinetic barrier was computed (ΔG‡298 = 21.0 kcal mol−1) that is in excellent agreement with experiment (≈ 22 kcal mol−1).12 It mainly arises from the entropic contribution with only a small electronic fraction. The transition state, BS(1,1)TSA (Fig. 4A), was located on the open shell singlet surface at a rather long N–N distance (2.06 Å). The spin is eventually quenched at significantly smaller N–N distances, leading to the closed shell ground state of dimer 3. The TS structure features the two Ir(PNP) monomer fragments being essentially perpendicular to each other. The two Ir–N groups are offset from collinearity along both relative orientations of the Ir(PNP) planes resulting in an out-of-plane zig-zag transition state. The orthogonalized sets of natural frontier orbitals represent hybrid orbitals from SOMO/LUMO mixing of the nitride monomers. π/π-Orbital overlap between the two monomers is further increased by slightly tilting the two monomer fragments towards each other (Ir–N–N–Ir dihedral angle: 152°). The previously computed PMe2-truncated model of reaction A also gave a good agreement for the driving force (ΔG298 = −108.2 kcal mol−1) but underestimation of the reaction barrier (ΔG‡298 = 13.7 kcal mol−1).12 In fact, the smaller model featured a distinctly different conformation for the TS with almost coplanar relative orientation of the two monomers. However, the minute electronic contributions to biradical coupling is almost identical for both models, while the entropic contribution seems underestimated in case of the small model presumably due to steric constraints for vibrational and rotational degrees of freedom of the tBu groups in the TS.
Fig. 4 Transition state structures of BS(1,1)TSA (A), 2TSB (B) and 1TSC (C) and schematic representations of the NBO interactions. |
Coupling of 22 and 14 to 25, was computed to be essentially electronically barrierless with a free energy of activation ΔG‡298 = 10.1 kcal mol−1, confirming the experimentally observed rapid coupling reaction B. The transition state 2TSB (Fig. 4B) exhibits strong structural similarities with BS(1,1)TSA concerning the N–N distance (DN–N = 2.00 Å) and the out-of-plane zig-zag relative orientation of the two nitride fragments (Ir–N–N–Ir dihedral angle: 142°). Accordingly, 2TSB also features a π/π-interaction between the two monomers, here arising from orbital overlap of the SOMO (2) and LUMO (4) monomer orbitals, respectively. The spin density is equally distributed over the mixed-valent dimer indicating strong electronic coupling and charge delocalization in the TS.
In contrast, reaction C is associated with a considerably higher computed kinetic barrier (ΔG‡298 = 29 kcal mol−1), which is in agreement with the enhanced thermal stability of 4. The transition state 1TSC (Fig. 4C) is located on the singlet surface. At a later stage there is a change to an open shell singlet and eventually the triplet ground state of 36 (ESI†). 1TSC exhibits significant differences compared with BS(1,1)TSA and 2TSB: The two Ir–N bond vectors are almost parallel (Ir–N–N–Ir dihedral angle: 178°) and the offset from collinearity only along one of the Ir(PNP) planes results in an in-plane zig-zag transition state. The closed-shell σ/π-orbital interaction between the two monomers resembles the TS previously considered by Burger.11 In the present case though different Ir–N π*-orbitals are involved on the two monomers, respectively, due to the perpendicular orientation of the two monomer fragments. Importantly, the N–N distance in 1TSC (1.84 Å) is considerably shorter compared with those in BS(1,1)TSA and 2TSB (both >2.00 Å) indicating that a stronger orbital interaction is necessary to overcome the kinetic barrier.
Hence, the computational evaluation confirms the order of relative coupling rates that was found experimentally (kB > kA > kC). Albeit partitioning of TS energies into the contributing terms was not done within this study, the TS structures allow for a qualitative rationalization that is in line with the observed order: Coulombic repulsion of the two positively charged nitrides is expected to contribute to the higher barrier of reaction Cvs.A. In turn, the TS of reaction B exhibits strong charge delocalization over both monomer moieties, tantamount a stabilizing Coulombic effect. Besides these charge considerations, barrierless coupling of 2 and 4 involves direct overlap of the respective SOMO and LUMO orbitals while diradical coupling reaction A necessitates some electronic rearrangement within BS(1,1)TSA with respect to 22. However, in all three TSs relatively large N–N distances were found and NBO analysis of 1TSC reflects only a donor–acceptor interaction between electronically weakly perturbed monomers. Hence, these results suggest for the Coulombic interaction to be the main contributor to the computed order of reaction barriers.
Irradiation was carried out with a Lot Oriel Xe-Short arc lamp. Cyclic voltammograms were recorded with a Metrohm Autolab PGSTAT101 using Ag/Ag+ reference-, glassy-carbon working- and Pt-wire counter-electrodes. Pt-net working-electrodes and a spectro-electrochemistry cuvette (0.1 cm optical pathway) were used for UV/VIS-spectroelectrochemical measurements in combination with an Avantes AvaSpec-2048×14 UV/VIS-spectrometer. Otherwise UV/VIS-spectra were recorded with an Agilent Cary 300 UV/VIS-spectrometer at room temperature with a 0.1 mM concentration in quartz cuvettes (1 cm optical pathway). NMR-spectra were recorded on a Bruker Avance 300/400 spectrometer and calibrated to the residual proton resonance of the solvent (CD2Cl2: δH = 5.32 ppm; δC = 54.00 ppm; Toluene-d8: δH = 7.09, δC = 137.86). 31P-NMR chemical shifts are reported relative to external phosphoric acid (δ = 0.0 ppm). Signal multiplicities are abbreviated as: s (singlet), h (heptet), br (broad). Experimental X-band EPR spectra were recorded on a Bruker ELEXSYS-II E500 CW-EPR. The spectra were simulated by iteration of the anisotropic g-values, (super)hyperfine coupling constants, and line widths using the EPR-simulation program W95EPR developed by Prof. Dr Frank Neese. ESI-MS-experiments employed a micrOTOF-Q II instrument (Bruker Daltonik). IR spectra were recorded with a Thermo Scientific Nicolet iZ10 FT/IR spectrometer at r.t. IR spectra in solution were recorded with a thin layer IR cell.
NBO analyses at the D3BJ-PBE0/def2-SV(P) level were performed with Gaussian 09 RevD.0127 coupled to NBO6.028 using the structures optimized as described above.
Footnote |
† Electronic supplementary information (ESI) available: Spectroscopic, electrochemical, kinetic, and mass spectrometric data and crystallographic and computational details. CCDC 1437698 and 1437699. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c5qi00267b |
This journal is © the Partner Organisations 2016 |