Stabilization of Linear C3 by Two Donor Ligands: A Theoretical Study of L-C3 -L (L=PPh3 , NHCMe , cAACMe )*.

Quantum chemical studies using density functional theory and ab initio methods have been carried out for the molecules L-C3 -L with L=PPh3 (1), NHCMe (2, NHC=N-heterocyclic carbene), and cAACMe (3, cAAC=cyclic (alkyl)(amino) carbene). The calculations predict that 1 and 2 have equilibrium geometries where the ligands are bonded with rather acute bonding angles at the linear C3 moiety. The phosphine adduct 1 has a synclinal (gauche) conformation whereas 2 exhibits a trans conformation of the ligands. In contrast, the compound 3 possesses a nearly linear arrangement of the carbene ligands at the C3 fragment. The bond dissociation energies of the ligands have the order 1<2<3. The bonding analysis using charge and energy decomposition methods suggests that 3 is best described as a cumulene with electron-sharing double bonds between neutral fragments (cAACMe )2 and C3 in the respective electronic quintet state yielding (cAACMe )=C3 =(cAACMe ). In contrast, 1 and 2 possess electron-sharing and dative bonds between positively charged ligands [(PPh3 )2 ]+ or [(NHCMe )2 ]+ and negatively charged [C3 ]- fragments in the respective doublet state.

Introduction

The stabilization of transient species, small clusters, and bare atoms through coordination with Lewis bases has become a highly prolific discipline in modern chemistry. Stable singlet carbenes like cyclic (alkyl)(amino) carbene (cAAC) and N‐heterocyclic carbene (NHC) have turned out to be excellent ligands for this purpose besides the well‐known phosphine groups. Focusing on the carbon atoms as the species to be stabilized, the first divalent carbon(0) complex, carbodiphosphorane C(PPh3)2 was synthesized by Ramirez et al. in 1961, followed by a proper structural characterization through X‐ray analysis in 1978 by Hardy et al. Later, theoretical studies on the bonding situation revealed that this complex can be best described in terms of donor–acceptor interactions Ph3P→C(0)←PPh3 between the phosphine ligands and a bare carbon atom in the excited 1D state. This finding led to the theoretical prediction of stable carbodicarbene C(NHC)2 in 2007, a benzoannulated derivative C(NHCBz)2 that was shortly afterwards synthesized by the groups of Bertrand and Fürstner. Then, gradually the bonding concept of carbones L→C(0)←L was extended for other donor ligands L. The notion of dative bonding in carbones was not undisputed, but it has been generally accepted now as a valid description of main‐group compounds. It should be noted that the use of arrows for dative bonds was originally introduced by Sidgwick in the 1920s and that the suggestion of divalent carbon(0) compounds L→C←L was suggested already by Varshavskii in 1980.

The bonding analysis also showed that molecules like (R2N)2C=C=C(NR2)2 and (cAAC)=C=(cAAC) can be better represented with classical electron‐sharing double bonds rather than dative bonds.[ , ] It should be noted that the donor–acceptor interactions comprise σ donation L→C(0)←L as well as π back‐donation L←C(0)→L, where the extent of back‐donation depends on the availability of low‐lying π acceptor orbitals of the ligand L. For example, a strong π acceptor such as CO induces strong π back‐donation, which means that the dative bonds in the bent equilibrium structure of carbon suboxide C3O2 should be described with the formula OC⇄C(0)⇄CO.[ , ]

Regarding its next larger carbon homolog, Bestmann and co‐workers prepared in 1989 a push–pull molecule, R3P‐C2‐BR3, with a dicarbon core. In retrospect, it is interesting that the authors suggested in that paper that carbodiphosphorane C(PPh3)2 can be formally regarded as a “complex” consisting of two donor molecules and an electron‐rich, excited carbon atom citing earlier work with similar suggestions, which were then thought to be a purely hypothetical formulation. The authors postulated that carbon chains with an even number n of carbon atoms would be stabilized by a push–pull mechanism D→C→A with a donor ligand D and acceptor ligand A at either end, whereas carbon chains with an odd number of carbon atoms are stabilized by two donor ligands D→C←D. The hypothesis was later confirmed by the synthesis of the C3 complex Ph3P‐C3‐PPh3 and the C4 adduct Ph2(Me)P‐C4‐B(CH2Ph)3. The latter species could be structurally characterized by X‐ray analysis whereas the former C3 complex could only be identified by NMR and IR spectroscopy, which were interpreted in terms of a trans‐bent geometry.

The possibility of stabilizing C2 by two donor NHC ligands was theoretically explored in 2012 by Dutton and Wilson. They also considered the PMe3 ligand in their study and reported that, although in (NHC)2C2 the central C4 unit is almost linear, in (PMe3)2C2, P‐C‐C‐P has a trans‐bent geometry. More recently, successful attempts to synthesize cAAC‐stabilized C2 as neutral species (cAAC)C2 were reported independently by the Roesky and Bertrand groups. A similar effort to synthesize NHCBz‐stabilized C2, by the deprotonation of a doubly protonated precursor, [(NHCBz)2C2H2]2+, failed presumably because of the small HOMO–LUMO energy gap of (NHCBz)2C2. Similar to its calculated NHC analog, the C4 core in (cAAC)C2 only slightly deviates from linearity. Although a comparative bonding analysis for donor–acceptor versus electron‐sharing is still missing, the linear C4 core moiety suggests an electron‐sharing C4 unit with electron‐sharing double bonds in (NHC)=C=C=(NHC) and (cAAC)=C=C=(cAAC) rather than dative bonds (NHC)→C2←(NHC) and (cAAC)→C2←(cAAC). Later, Wu et al. isolated the 4‐pyridylidene‐supported C2 species, which exhibits a more significant deviation from linearity by 24.2° than that observed in NHC or cAAC analogs.

The stabilization of longer carbon chains C where n>2 in compounds R‐C‐R′ with various groups R has been the topic of numerous studies. The synthetic efforts to isolate polyynes that are stabilized by metal or non‐metal end groups R has been summarized in a seminal review by Szafert and Gladys in 2003. Research about longer cumulenes has been reviewed more recently by Tykwinski and co‐workers. Very little is known about ligand‐stabilized C species with donor ligands beyond the C2 unit. As mentioned above, the complex Ph3P‐C3‐PPh3 has been synthesized in the equilibrium geometry and the bonding situation is not known so far. Free C3 is known as an interstellar species having 1Σg + ground state, which was spectroscopically studied in the gas phase.

In the present study, we carried out a density functional theory (DFT) study on (L)2C3 complexes where L=PPh3, NHCMe, and cAACMe. The bonding situation is analyzed through natural bond orbital (NBO) and charge and energy decomposition analyses (EDA‐NOCV). The present results suggest that the phosphine complex Ph3P‐C3‐PPh3 unexpectedly has a synclinal (gauche) conformation and not a trans conformation as originally assumed. The calculations also suggest that the title molecules L‐C3‐L with L=NHCMe and cAACMe are also viable to be realized experimentally at ambient conditions. The complex (NHCMe)2C3 has a trans‐bent equilibrium geometry whereas the adduct (cAACMe)2C3 possesses a nearly linear arrangement of the central C4 core. The bonding situation of the latter complex is distinctively different from the other two adducts.

Computational Details

The geometry optimization followed by the vibrational frequencies calculations of (L)2C3 molecules where L=NHCMe, cAACMe, and PPh3 in their respective singlet and triplet electronic states were performed at the BP86‐D3(BJ)/def2‐TZVPP level. All the molecules are minima on the potential energy surfaces as ensured by the absence of any imaginary frequencies. For some molecules, additional calculations were carried out with the M06‐2X functional developed by Truhlar and co‐workers. These calculations were carried out with the Gaussian 16 program package. The NBO analysis was performed to evaluate the partial charge, Wiberg bond indices (WBI), and natural bond orbitals by using the NBO 6.0 program. Unless otherwise noted, we discuss geometries and energies at the BP86‐D3(BJ)/def2‐TZVPP level. Some ab initio calculations were carried out with the DLPNO‐CCSD(T) approach in conjunction with the cc‐pVTZ basis set by using the program package ORCA.

The bonding situations were analyzed by means of an energy decomposition analysis (EDA) together with the natural orbitals for chemical valence (NOCV) method by using the ADF 2018.105 program package. The EDA‐NOCV calculations were carried out at the BP86‐D3(BJ)/TZ2P level by using the BP86‐D3(BJ)/def2‐TZVPP optimized geometries. In this analysis, the intrinsic interaction energy (ΔΕ int) between two fragments can be divided into four energy components as follows [Eq. 1]:

Although the electrostatic ΔE elstat term represents the quasi‐classical electrostatic interaction between the unperturbed charge distributions of the prepared fragments, the Pauli repulsion ΔE Pauli corresponds to the energy change associated with the transformation from the superposition of the unperturbed electron densities of the isolated fragments to the wavefunction, which properly obeys the Pauli principle through explicit antisymmetrization and renormalization of the production wavefunction. As we included D3(BJ), it provides us with the dispersion interaction energy between two interacting fragments ΔE disp. The orbital term ΔE orb originates from the mixing of orbitals, charge transfer, and polarization between the isolated fragments, which can be further decomposed into contributions from each irreducible representation of the point group of the interacting system as follows [Eq. 2]:

The combination of the EDA with NOCV enables the partition of the total orbital interactions into pairwise contributions of the orbital interactions, which is vital to get a complete picture of the bonding. The charge deformation Δρ(r), resulting from the mixing of the orbital pairs ψ k(r) and ψ −k(r) of the interacting fragments presents the amount and the shape of the charge flow owing to the orbital interactions [Eq. 3], and the associated energy term ΔE orb provides the size of stabilizing orbital energy originating from such interactions [Eq. 4].

More details about the EDA‐NOCV method and its application are given in recent review articles.

Results and Discussion

Figure 1 shows the optimized geometries of the three (L)2C3 molecules with L=PPh3 (1), NHCMe (2), cAACMe (3), and the three free ligands L and the C3 moiety in the electronic singlet ground states. The three compounds show distinct differences in their equilibrium geometries. The molecules 1 and 2 have rather acute bonding angles L‐C1‐C2 of 124° and 139° whereas 3 exhibits a nearly linear arrangement with a bending angle L‐C1‐C2 of 171°. Another difference between the carbene complexes concerns the planes of the ligands in 2 and 3. The two NHCMe ligands in 2 are in the same plane as C3, whereas the two cAACMe ligands in 3 are perpendicular to each other, as can be seen from the torsional angle, t(N′‐C′‐C′′‐N′′) of 178.7° in the first molecule and 99.9° in the second complex.

Figure 1

Optimized geometries of (L)2C3 (L=PPh3, NHCMe, and cAACMe) and fragments at the BP86‐D3(BJ)/def2‐TZVPP level. Bond lengths are given in Å, angles in degrees. Calculated bond dissociation energies D e and free energies ΔG 298 for the reactions (L)2C3→2 L+C3.

The phosphine ligands L around the central C3 component in the energetically lowest lying form of 1 unexpectedly have a synclinal (gauche) conformation and not a trans conformation as originally suggested. The calculations at the BP86‐D3(BJ)/def2‐TZVPP level predict that the trans conformation 1 a is 8.9 kcal mol−1 higher in energy than the gauche form 1. The optimization of both isomers without the dispersion term at the BP86/def2‐TZVPP level suggests that 1 is only 0.4 kcal mol−1 lower in energy than 1 a, which suggests that dispersion interactions significantly contribute to the stability of 1. We further optimized both conformations by using the meta‐hybrid functional M06‐2X using a def2‐TZVPP basis set with and without explicit consideration of dispersion interactions using the D3 term of Grimme. The calculations at the M06‐2X‐D3/def2‐TZVPP level also predict that 1 is 7.0 kcal mol−1 lower in energy than 1 a whereas M06‐2X/def2‐TZVPP favors 1 over 1 a by 6.0 kcal mol−1. According to the DFT calculations, (PPh3)2C3 (1) has a gauche equilibrium geometry whereas (NHCMe)2C3 (2) has a trans‐bent conformation as the energetically lowest‐lying form. The same energy ordering is predicted by ab initio calculations at the DLPNO‐CCSD(T)/cc‐pVTZ level using the BP86‐D3(BJ)/def2‐TZVPP optimized geometries, which give the gauche isomer 1 as being 7.7 kcal mol−1 lower in energy than the trans isomer 1 a. For all three molecules, the corresponding triplet state is considerably higher in energy by 30.8–41.7 kcal mol−1 than the singlet state (see Figure S1 in the Supporting Information)

The experimental work on (PPh3)2C3 (1) reported the 31P NMR spectrum of the complex, which has a signal at δ=−10.4 ppm. We calculated the 31P NMR chemical shift values of 1 and 1 a by using the GIAO approach at the BP86/def2‐TZVPP level. The calculated chemical shifts for the 31P NMR spectrum are δ=+11.3 and +17.5 ppm for the non‐degenerate P atoms of 1 and δ=−5.6 ppm for the degenerate P of 1 a. The latter value is in much better agreement with the experimental result, which was recorded in C6D6. It is possible that in solution 1 a prevails over 1 or that a fast equilibrium between the gauche conformations leads to a 31P NMR spectrum that mimics a trans form. The calculation of 1 and 1 a with inclusion of the solvent effect of C6D6 by using the CPCM approach at the BP86‐D3(BJ)/def2‐TZVPP level even showed a slightly higher preference of 1 over 1 a by 9.2 kcal mol−1.

The comparison of the equilibrium geometries of 1–3 with the free ligands shows that the C−C bonds of the central C3 moiety in 1 (1.299/1.305 Å) and 2 (1.299 Å) have essentially the same length as in free C3 (1.301 Å) whereas 3 possesses somewhat shorter bonds (1.282 Å). The geometries suggest that the l−C3 bonds in 1 and 2 may have a different character than in 3. This is supported by the calculated bond strength of the three different ligands to the C3 fragment, which is also given in Figure 1. The cAACMe ligands in 3 are much more strongly bonded (D e=263.9 kcal mol−1; ΔG 298=233.5 kcal mol−1) than the NHCMe ligands in 2 (D e=166.9 kcal mol−1; ΔG 298=137.6 kcal mol−1) and the PPh3 ligands in 1 (D e=126.5 kcal mol−1; ΔG 298=92.7 kcal mol−1). The calculated bond dissociation energies suggest that 2 and particularly 3 are thermodynamically more stable than 1.

We analyzed the electronic structure and bonding situation of the three molecules with charge and energy decomposition methods. Table 1 shows the NBO results of the adducts. The NBO method finds a single P−C1 σ‐bond between the phosphine ligands and C3 in 1 whereas double bonds CL=C1 are calculated for the ligands in 2 and 3. The P−C1 σ‐bond in 1 is polarized toward the carbon end whereas the CL=C1 σ‐bonds in 2 and 3 are slightly polarized toward the ligand atom. The CL=C1 π‐bonds in 2 and 3 exhibit a small polarization in the opposite direction towards C3. The charge distribution suggests that the central carbon atom C2 of the central C3 moiety is essentially neutral whereas the terminal carbon atoms carry negative charges, which indicates the total charge donation L→C3←L in the order 1>2>3.

Table 1

NBO results of the compounds (L)2C3 at the BP86/def2‐TZVPP level of theory. Occupation number ON, polarization, and hybridization of the L−C1 bonds and partial charges q.

Complex	Bond	ON	l−C1 polarization		q
			and hybridization [%]		C1	C2	C3
1	P−C1 σ	1.96	P: 40.6 s(30.0), p(69.4)	C1: 59.4 s(49.0), p(51.0)	−0.80	0.09	−0.79
2	CL−C1 σ	1.97	CL: 52.6 s(42.1), p(57.9)	C1: 47.4 s(47.6), p(52.4)	−0.42	0.05	−0.42
	CL−C1 π	1.65	CL: 46.7 s(0.1), p(99.9)	C1: 53.3 s(0.1), p(99.9)
3	CL−C1 σ	1.97	CL: 51.2 s(38.6), p(61.4)	C1: 48.8 s(50.1), p(49.9)	−0.19	−0.01	−0.19
	CL−C1 π	1.80	CL: 46.1 s(0.1), p(99.9)	C1: 53.9 s(0.1), p(99.9)

Detailed insight into the nature of the ligand–C3 interactions is provided by the EDA‐NOCV approach. The versatility of the method makes it possible to distinguish (a) between the total interactions of the isolated fragments in their original electronic states, charges, and equilibrium geometries, and (b) between the interactions of the fragments in the electronic states, charges, and geometries that are present in the finally formed molecule. For example, the electronic structure of diatomic LiF can be analyzed either by using the neutral atoms Li and F as interacting fragments in the EDA‐NOCV calculations, taking into account the charge migration Li→F during bond formation. Alternatively, one can use the finally formed ions Li+ and F− as fragments, which are better suited to calculate the interactions in the eventually formed bond. Those fragments that give the smallest change in the orbital interactions ΔE orb are best suited to describe the final bond. This approach has proven to be very useful in cases where the electronic state or charge of the fragments is not directly apparent.

In the present case of the molecules l‐C3‐L, we focus on the best description of the l−C bonds in 1–3, which can be sketched with dative or electron‐sharing bonds between neutral or charged fragments. We considered four different pairs of interacting species L2 q and C3 q in different electronic states, which may be considered as chemically reasonable species, which lead to the molecules l‐C3‐L. These are (a) the neutral fragments L2 and C3 in the electronic singlet (S) state forming dative bonds; (b) the neutral fragments L2 and C3 in the electronic quintet (Q) state forming electron‐sharing double bonds; (c) the singly charged fragments L2 + and C3 − in the electronic doublet (D) state forming electron‐sharing and dative bonds; (d) the doubly charged fragments L2 2+ and C3 2− in the electronic triplet (T) state forming electron‐sharing and dative bonds. Table 2 gives the numerical results of the EDA‐NOCV calculations of 1–3 using the four fragments pairs (a)–(d) as interacting species.

Table 2

EDA‐NOCV results of (L)2C3 molecules using four different sets of fragments with different charges and electronic states (S=singlet, D=doublet, T=triplet, Q=quintet) and associated bond types at the BP86‐D3(BJ)/TZ2P level. Energies are in kcal mol−1. The most favorable fragmentation scheme and bond type is given by the smallest ΔE orb value written in bold.

Molecule	Bond type[a]	Fragments	ΔE int	ΔE Pauli	ΔE elstat	ΔE disp	ΔE orb
(PPh3)2C3	D	(PPh3)2(S)+C3 (S)	−488.5	674.6	−428.1	−16.3	−718.6
	E	(PPh3)2(Q)+C3 (Q)	−461.1	630.4	−371.4	−16.3	−703.7
	D+E	[(PPh3)2]+(D)+[C3]− (D)	−493.8	766.5	−565.5	−16.3	−678.5
	D+E	[(PPh3)2]2+(T)+[C3]2− (T)	−824.4	967.8	−918.8	−16.3	−857.4
(NHCMe)2C3	D	(NHCMe)2(S)+C3 (S)	−541.0	634.7	−466.2	−9.4	−700.1
	E	(NHCMe)2(Q)+C3 (Q)	−544.8	506.4	−373.8	−9.4	−667.9
	D+E	[(NHCMe)2]+(D)+[C3]− (D)	−551.9	621.5	−517.2	−9.4	−646.7
	D+E	[(NHCMe)2]2+(T)+[C3]2− (T)	−875.9	659.2	−798.2	−9.4	−727.5
(cAACMe)2C3	D	(cAACMe)2(S)+C3 (S)	−594.2	419.5	−359.6	−10.6	−643.5
	E	(cAACMe)2(Q)+C3 (Q)	−407.7	519.1	−348.5	−10.6	−567.6
	D+E	[(cAACMe)2]+(D)+[C3]− (D)	−586.8	555.8	−445.4	−10.6	−686.5
	D+E	[(cAACMe)2]2+(T)+[C3]2− (T)	−930.1	873.0	−819.7	−10.6	−972.7

[a] D=Dative bond; E=electron‐sharing bond.

The results in Table 2 suggest that the best description of the l‐C3‐L bonds for L=PPh3 (1) and L=NHCMe (2) comes from the interactions of the singly charged fragments L2 + and C3 − in the electronic doublet (D) state forming electron‐sharing and dative bonds. In contrast, the smallest ΔE orb value and thus the best description of the l‐C3‐L bonds for L=cAACMe (3) is provided by the neutral fragments L2 and C3 in the electronic quintet (Q) state forming electron‐sharing double bonds. This agrees with the NBO charges in Table 1, which suggest that the partial charges of the interacting fragments in 1 and 2 are clearly higher than in 3. We further analyzed the orbital interactions ΔE orb of 1–3 by using the most favorable fragments and studied the main pairwise orbital contributions ΔE orb(. The numerical results are given in Table 3.

Table 3

The EDA‐NOCV results at the BP86‐D3(BJ)/TZ2P level of L2C3 molecules 1–3 using (L2)+ and (C3)− in the electronic doublet (D) states as interacting fragments for 1 and 2 and neutral L2 and C3 in the quintet (Q) states as interacting fragments for 3. Energies are in kcal mol−1.

Energy	Interaction[c]	[(PPh3)2]+(D) +[C3]−(D)	[(NHCMe)2]+ (D) +[C3]−(D)	(cAACMe)2(Q) +C3 (Q)
ΔE int		−493.8	−551.9	−407.7
ΔE Pauli		766.5	621.5	519.1
ΔE disp [a]		−16.3 (1.2 %)	−9.4 (0.8 %)	−10.6 (1.1 %)
ΔE elstat [a]		−565.5 (44.8 %)	−517.2 (44.0 %)	−348.5 (37.6 %)
ΔE orb [a]		−678.5 (53.8 %)	−646.7 (55.1 %)	−567.6 (61.2 %)
ΔE orb(1) [b]	l‐C3‐L σ e− sharing (+,+)			−187.6 (33.1 %)
	L→C3←L σ donation (+,+)	−286.3 (42.1 %)	−235.1 (36.3 %)	–
ΔE orb(2) [b]	l‐C3‐L σ e− sharing (+,−)	−231.3 (34.1 %)	−208.7 (32.2 %)	−180.3 (31.8 %)
	L→C3←L σ donation (+,−)			–
ΔE orb(3) [b]	l‐C3‐L π/σ e− sharing			−77.7 (13.7 %)
	L←C3→L σ back donation	−51.8 (7.6 %)	−88.3 (13.6 %)	–
ΔE orb(4) [b]	l‐C3‐L π e− sharing			−76.7 (13.5 %)
	L←C3→L π back‐donation	−40.4 (5.9 %)	−52.9 (8.1 %)	–
ΔE orb(5) [b]	L←C3→L π back‐donation	−25.0 (3.6 %)	−18.4 (2.8 %)	–
ΔE orb(rest) [b]		−43.6 (6.4 %)	−43.3 (6.6 %)	−45.3 (8.0 %)

[a] The values in the parentheses show the contribution to the total attractive interaction ΔE elstat+ΔE orb+ΔE disp. [b] The values in parentheses show the contribution to the total orbital interaction ΔE orb.

–

–

–

–

–

There are five relevant orbital interactions ΔE orb(1)–ΔE orb(5) for molecules 1 and 2 and four terms ΔE orb(1)–ΔE orb(4) for complex 3. The nature of the orbital terms and the involved fragments orbitals can be identified with the help of the associate deformation densities Δρ and connected MOs that are shown in Figure 2, Figure 3, and Figure 4. The strongest orbital interaction ΔE orb(1) in 1 and 2 is due to the σ donation of the HOMO of the ligands L2 + into the LUMO of C3 −. It is remarkable and counter‐intuitive that the cation serves as donor and the anion as acceptor. However, this must be seen as only one part of the total orbital interactions. The slightly weaker orbital term ΔE orb(2) in 1 and 2 comes from the electron‐sharing interactions of the unpaired electrons in the SOMO of the fragments. The remaining three orbital terms ΔE orb(3)–ΔE orb(5) in the two molecules are due the π back‐donation of the degenerate HOMO and the HOMO−1 of C3 − to vacant orbitals of L2 +. The strength of the π back‐donation is much less than that of the σ donation ΔE orb(1) but the amount of charge migration, which is given by the eigenvalues υ, shows that π back‐donation has a similar magnitude as the σ donation, which is consistent with the NBO charges (Table 1). The amount of charge migration given by υ does not correlate with the associated stabilization energy ΔE orb(.

Figure 2

The shape of the deformation densities Δρ (1)–(5) that correspond to ΔE orb(1)–(5), and the associated MOs of (PPh3)2C3 and the fragments orbitals of [(PPh3)2]+ and [C3]− in the doublet state at the BP86‐D3(BJ)/TZ2P level. Isosurface values are 0.003 au. The eigenvalues |ν| give the size of the charge migration in e. The direction of the charge flow of the deformation densities is red→blue.

Figure 3

The shape of the deformation densities Δρ (1)–(5) that correspond to ΔE orb(1)–(5), and the associated MOs of (NHCMe)2C3 and the fragments orbitals of [(NHCMe)2]+ and [C3]− in the doublet state at the BP86‐D3(BJ)/TZ2P level. Isosurface values are 0.003 au. The eigenvalues |ν| give the size of the charge migration in e. The direction of the charge flow of the deformation densities is red→blue.

Figure 4

The shape of the deformation densities Δρ (1)–(5) that correspond to ΔE orb(1)–(5), and the associated MOs of (cAACMe)2C3 and the fragments orbitals of (cAACMe) and C3 in the quintet state at the BP86‐D3(BJ)/TZ2P level. Isosurface values are 0.003 au. The eigenvalues |ν| give the size of the charge migration in e. The direction of the charge flow of the deformation densities is red→blue.

The four significant orbital terms ΔE orb(1)–ΔE orb(4) between the neutral fragments in the quintet states in complex 3 are easily recognized as the strong in‐phase (+,+) and out‐of‐phase (+,−) electron‐sharing interactions ΔE orb(1) and ΔE orb(2) and the weaker but still rather strong orthogonal π interactions ΔE orb(3) and ΔE orb(4), which are slightly different, because the molecular planes of the cAACMe ligands are not perfectly orthogonal to each other.

Summary and Conclusion

The present study may be summarized as follows. Quantum chemical calculations suggest that the molecules l‐C3‐L with L=PPh3 (1) and NHCMe (2) have equilibrium geometries where the ligands are bonded with rather acute bonding angles at the linear C3 moiety. The phosphine adduct 1 has a synclinal (gauche) conformation whereas 2 exhibits a trans conformation of the ligands. In contrast, the complex with L=cAACMe (3) possesses a nearly linear arrangement of the carbene ligands at the C3 fragment. The bond dissociation energies of the ligands have the order 1<2<3. The bonding analysis using charge and energy decomposition methods suggests that 3 is best described as a cumulene with electron‐sharing double bonds between neutral fragments (cAACMe)2 and C3 in the respective electronic quintet state yielding (cAACMe)=C3=(cAACMe). In contrast, 1 and 2 possess electron‐sharing and dative bonds between positively charged ligands [(PPh3)2]+ or [(NHCMe)2]+ and negatively charged [C3]− fragments in the respective doublet state.

Conflict of interest

The authors declare no conflict of interest.

As a service to our authors and readers, this journal provides supporting information supplied by the authors. Such materials are peer reviewed and may be re‐organized for online delivery, but are not copy‐edited or typeset. Technical support issues arising from supporting information (other than missing files) should be addressed to the authors.

Supplementary

Click here for additional data file.