Angewandte
Chemie
The description of N2(PPh3)2 with two lone-pairs at each
nitrogen atom suggests that the compound should be a strong
proton acceptor and that even the second proton affinity,
where the same nitrogen atom is protonated, should be rather
high. This idea is verified by the calculated proton affinities
(PAs). The first PA of N2(PPh3)2 is 245.2 kcalmolꢀ1 which is
larger than the PA of most amines.[16] The second PA yielding
Ph3PN(H+)2NPPh3 is 170.5 kcalmolꢀ1 which is only a bit
smaller than the second PA yielding Ph3PN(H+)N(H+)PPh3
(195.6 kcalmolꢀ1). It has been shown that the second PA is
a sensitive probe for the occurrence of two lone pairs in
carbon bases.[17] The calculated second PA for one nitrogen
atom is a further evidence for the existence of two lone-pair
orbitals at nitrogen.[18]
a RASSCF level which is not possible for us. However, our
calculations do provide information about the dissociation
pathway. We calculated N2(PPh3)2 at RI-BP86/def2-TZVPP
ꢀ
where one P N separation was fixed at longer distances than
ꢀ
the equilibrium values with intervals of 0.1 ꢀ. The second P
ꢀ
N bond becomes longer while the N N bond becomes clearly
shorter up to a value of 2.0 ꢀ for the frozen P N bond. At this
ꢀ
ꢀ
point, the second P N bond length was 1.606 ꢀ, the N–N
separation was 1.379 ꢀ, and the energy was 25.1 kcalmolꢀ1
higher than the optimized structure. When the P–N frozen
ꢀ
separation was further elongated to 2.1 ꢀ, the second P N
bond broke and the associated PPh3 ligand dissociated. A
geometry optimization of a possible intermediate N2(PPh3)
did not give an equilibrium structure. The calculations suggest
that the reaction N2(PPh3)2!N2 + 2PPh3 which proceeds as
Why is triphenylphosphinazine an amazingly stable com-
pound although the dissociation reaction into N2 + 2PPh3 is
strongly exergonic? As shown above, the intrinsic interaction
energy DEint between the fragments N2 in the (1)1Gg reference
state and PPh3 in the frozen geometry of N2(PPh3)2 is very
large (ꢀ300.1 kcalmolꢀ1). Thus, dinitrogen in the (1)1Gg state
is a very powerful Lewis acid which is due to two reasons. One
reason is the rather high electronegativity of nitrogen while
the other reason is that the nitrogen atoms in jN ꢀ Nj have
only an electron sextet. However, the large interaction energy
does not compensate for the relaxation of the fragments into
the equilibrium geometries and electronic ground states
which are calculated to be 8.3 kcalmolꢀ1 for both phosphine
ligands and 341.0 kcalmolꢀ1 for the process N2 (11Gg)fr!
ꢀ
concerted but not necessarily synchronous rupture of the P N
bonds with a rather high barrier. The DFT calculations
suggest that the activation barrier for the concerted but not
necessarily synchronous reaction is higher than 25.1 kcal
molꢀ1 but lower than 67.6 kcalmolꢀ1 which is the bond
ꢀ
dissociation energy for the N N bond of N2(PPh3)2 yielding
NPPh3 in the electronic doublet ground state [Eq. (2)].
ð2Þ
N2ðPPh3Þ2 ! 2 NPPh3
þ 67:6 kcal molꢀ1
The above results suggest that donor–acceptor interac-
tions[20] in main-group compounds of atoms of the first octal-
row may be more important than hitherto realized. In 2006 it
was recognized by one of us[21] that carbodiphosphoranes
C(PR3)2 which were synthesized as early as 1961[22] are
examples for the compound class of carbones CL2 which have
N2(X1Sg ). High-level ab initio calculations predict that the
+
(1)1Gg state of N2 has an equilibrium separation of 1.608 ꢀ and
is 294.3 kcalmolꢀ1 above the X1Sg ground state.[19] The RI-
+
BP86/def2-TZVPP calculations for N2 (1)1Gg give a N–N
separation of 1.661 ꢀ and an excitation energy of 328.6 kcal
molꢀ1. This is a remarkably good agreement with the ab initio
result in light of the approximations of the DFT method. The
N–N separation becomes shorter in N2(PPh3)2 because the
!
donor–acceptor bonds L!C L to a naked carbon atom in
the singlet 1D state that has two lone pairs, and is thus
a double Lewis base. This was the starting point for intensive
studies which led to the synthesis of new carbones and
carbone complexes.[23] Herein we show that that triphenyl-
phosphinazine which was isolated in 1964[10] is a donor–
ꢀ
donation into the N N bonding orbital 1pu’ is likely to affect
the N–N distance more strongly than the donation into the
antibonding 1pg’ acceptor orbital (See Figure 1d and
Figure 3). The 1pg’ orbital mixes with the 2sg orbital of N2
!
acceptor complex of dinitrogen Ph3P!N2 PPh3 which can
be considered as an example of dinitrogen complexes L!
ꢀ
(Figure 1) which decreases the N N antibonding character.
!
N2 L. We think that triphenylphosphinazine is an excep-
The calculations thus suggest that the very strong Ph3P!
tional compound, because it melts at 1848 and it decomposes
into N2 and PPh3 only above 2158 although the associated
bond breaking reaction is calculated to be exergonic by 74.5–
87.7 kcalmolꢀ1. Further studies might lead to the synthesis of
complexes N2(L)2 with other ligands L and to transition-metal
complexes in which the nitrogen atoms of N2(L)2 serve as
double donors to the metal. Also, it has not escaped our
attention that the direct synthesis of N2(PPh3)2 from N2 and
PPh3 would be a significant contribution to the topics of
nitrogen activation and chemical energy storage. The results
which are reported herein are a challenge for experiment.
!
N2(1)1Gg PPh3 attraction makes triphenylphosphinazine
a
kinetically stable compound, because the activation
energy for breaking the donor–acceptor bonds is very high.
We carried out preliminary calculations to estimate the
activation barrier for dissociation of the phosphine ligands
but we experienced severe convergence problems. The
transition state for the dissociation reaction N2(PPh3)2!
N2 + 2PPh3 which involves a change in the configuration of
+
the N2 moiety from the (1)1Gg excited state to the X1Sg
ground state can be optimized only at a multi-reference
level. This is a formidable task which requires substantial
computational resources, because a size reduction of the
molecule to a model system such as N2(PH3)2 is not feasible
because hydrogen migration may take place during the
geometry optimization. A calculation of the transition state
for the molecule which has C1 symmetry involves the
optimization of 3Nꢀ6 = 214 geometry variables at
Received: August 6, 2012
Revised: October 16, 2012
Published online: February 4, 2013
Keywords: bonding analysis · DFT calculations ·
dinitrogen complex · donor–acceptor interactions
.
ꢀ 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim
Angew. Chem. Int. Ed. 2013, 52, 3004 –3008