Stabilities, Electronic Structures, and Bonding Properties of Iron Complexes (E1E2)Fe(CO)2(CNArTripp2)2 (E1E2=BF, CO, N2, CN-, or NO+).

The coordination of 10-electron diatomic ligands (BF, CO N2) to iron complexes Fe(CO)2(CNArTripp2)2 [ArTripp2=2,6-(2,4,6-(iso-propyl)3C6H2)2C6H3] have been realized in experiments very recently (Science, 2019, 363, 1203-1205). Herein, the stability, electronic structures, and bonding properties of (E1E2)Fe-(CO)2(CNArTripp2)2 (E1E2=BF, CO, N2, CN-, NO+) were studied using density functional (DFT) calculations. The ground state of all those molecules is singlet and the calculated geometries are in excellent agreement with the experimental values. The natural bond orbital analysis revealed that Fe is negatively charged while E1 possesses positive charges. By employing the energy decomposition analysis, the bonding nature of the E2E1-Fe(CO)2(CNArTripp2)2 bond was disclosed to be the classic dative bond E2E1→Fe(CO)2(CNArTripp2)2 rather than the electron-sharing double bond. More interestingly, the bonding strength between BF and Fe(CO)2(CNArTripp2)2 is much stronger than that between CO (or N2) and Fe(CO)2(CNArTripp2)2, which is ascribed to the better σ-donation and π back-donations. However, the orbital interactions in CN-→Fe(CO)2(CNArTripp2)2 and NO+→Fe(CO)2(CNArTripp2)2 mainly come from σ-donation and π back-donation, respectively. The different contributions from σ donation and π donation for different ligands can be well explained by using the energy levels of E1E2 and Fe(CO)2(CNArTripp2)2 fragments.

Introduction

Carbon monoxide (CO), which is a diatomic molecule with 10 valence electrons, is among the most widely studied ligands in organometallic chemistry. Carbonyl complexes, especially the homoleptic carbonyl complexes, are prototypical examples of the Dewar–Chatt–Duncanson (DCD) model.[ , ] By using the DCD model, the bonding interaction between CO and transition metals could be mainly classified as two components, including σ‐donation from the lone pair of CO to an empty orbital on the metal and π back‐donation from occupied d‐orbitals of the metal to the π* orbitals of CO. Although the former usually plays as a primary bonding interaction, π back‐donations are important and sometimes are the dominant contribution in bonding interactions.[ , ] Several other 10‐electron diatomic ligands are isoelectronic and isolobal to CO, including neutral ligands N2 and BF and ionic ligands CN−, and NO+. Compared with CO, N2 has lower‐energy σ‐donor orbitals and higher‐energy π* orbitals, revealing its lower binding ability. Both CN– and NO+ are weak ligands because of less π‐acceptor ability and σ‐donor ability, respectively. The binding nature of N2, CN− and NO+ to transition metals has been characterized experimentally and theoretically.[ , ]

The boron monofluoride (BF) has a decreased HOMO‐LUMO gap (HOMO, highest occupied molecular orbital; LUMO, lowest unoccupied molecular orbital) compared to CO. On one hand, the resulted in better σ‐donor and π‐acceptor abilities lead to more favorably ligating properties to transition metals. On the other hand, the small HOMO‐LUMO gap gives rise to more reactivity and instability at room temperature. Lots of efforts have been devoted to fluoroborylene complexes since Timms postulated the first fluoroborylene complex [Fe(BF)(PF3)4] in 1973.[ , , ] Braunschweig et al. detected the bridged fluoroborylene complex [{(OC)5Mn}2(μ‐BF)] by using 11B NMR spectrum. A synthesis and structural characterization of fluoroborylene ruthenium complex Cp2Ru2(CO)4(μ‐BF), in which BF ligand is bridged between transition metal centers, were reported by Vidovic and Aldridge.[ , ] By using reactions of laser‐ablated metal atoms with BF3, Wang, Andrews, and co‐workers synthesized terminal fluoroborylene complexes FBMF2 (M=transition metal) in the gas phase.[ , , ] Fluoroborylene iron carbonyls and homoleptic fluoroborylenes complexes have also been proposed theoretically by Hoffmann and Schaefer and King.[ , , , , , , ]

Very recently, Figueroa et al. successfully realized the isolation of the iron complex (BF)Fe(CO)2(CNArTripp2)2 [ArTripp2=2,6‐(2,4,6‐(iso‐propyl)3C6H2)2C6H3] with a terminal BF ligand, along with the isoelectronic dinitrogen and CO complexes (N2)Fe(CO)2(CNArTripp2)2 and Fe(CO)3(CNArTripp2)2. Further single‐crystal x‐ray diffraction, spectroscopic, and electron‐density topology calculation studies demonstrated that the terminal BF ligand possesses particularly strong σ‐donor and π‐acceptor properties. However, the nature of the chemical bond could be altered from a dative bond to an electron‐sharing double bond if the π‐bonding interaction is strong enough. Therefore, the nature of B−Fe bond in (BF)Fe(CO)2(CNArTripp2)2 is unclear. Moreover, the σ‐bonding interaction and π‐bonding interaction strengths, which play as the dominant bonding interaction between Fe and BF, remain unknown.

In the present study, we report a theoretical study on geometries, electronic structures, and bonding properties of the 10‐electron diatomic ligand to iron complexes (E1E2)Fe(CO)2(CNArTripp2)2 (E1E2=BF, CO, N2, CN−, NO+) by using DFT calculations. The nature of the Fe–E1 bonds has been clarified. The bonding strengths of σ‐interaction and π‐interaction were analyzed and well explained by using the energy levels of 10‐electron diatomic ligand.

Computational Methods

Test calculations on (BF)Fe(CO)2(CNArTripp2)2 by using BP86,[ , ] B3LYP,[ , ] BP86‐D3(BJ), and ωB97XD functional with def2‐SVP basis set shows that both BP86 and ωB97XD could produce the experimental structure parameters well (please see ESI for details). Regarding the large calculation task in the present work, the BP86/def2‐SVP was used to optimize all structures without any symmetry restriction. Frequency calculations were performed at the same level of theory to verify that all the structures are local minima. The functional dependence of the optimized structures and the spin states was examined as shown in Tables S1 and S2. The above DFT calculations were conducted using Gaussian 09 software. The NBO partial charges and Wiberg bond orders were computed at BP86/def2‐TZVPP using NBO 3.1 as implemented in Gaussian 09.

The nature of the chemical bond E−Fe was investigated by means of an energy decomposition analysis (EDA) developed by Ziegler and Rauk. The EDA focuses on the instantaneous interaction energy ΔE int of a bond A–B between the fragments A and B in the particular electronic reference state at the frozen geometry of the molecule AB.[ , , , ] The interaction energy ΔE int is divided into three main components [eqn (1)]. In the present calculations, the dispersion correction ΔE disp term was also computed:

The term ΔE elstat corresponds to the quasi‐classical electrostatic interaction between the unperturbed charge distributions of the prepared atoms and is usually attractive. The Pauli repulsion ΔE Pauli arises as the energy change associated with the transformation from the superposition of the unperturbed electron densities ρ A + ρ B of the isolated fragments to the wavefunction which properly obeys the Pauli principle through explicit antisymmetrization (Â operator) and renormalization (N=constant) of the product wavefunction. ΔE Pauli comprises the destabilizing interactions between electrons of the same spin on either fragment. The orbital interaction ΔE orb accounts for charge transfer, polarization effects, and electron‐pair bonding. Finally, the dispersion interaction ΔE disp is involved with the pairwise correction (DFT−D3) from Grimme.

The EDA‐NOCV[ , ] method combines the EDA with the natural orbitals for chemical valence (NOCV)[ , , ] to decompose the orbital interaction term ΔE orb into pairwise contributions. The NOCVs are defined as the eigenvector of the valence operator, V given by eqn 2.

In the EDA‐NOCV scheme the orbital interaction term, ΔE orb is given by eqn 3,

in which and are diagonal transition state Kohn–Sham matrix elements corresponding to NOCVs with the eigenvalues and , respectively. The term for a particular type of bond is assigned by visual inspection of the shape of the deformation density Δρ. The latter term is a measure of the size of the charge deformation and it provides a visual notion of the charge flow that is associated with the pairwise orbital interaction. The EDA‐NOCV scheme thus provides both qualitative and quantitative information about the strength of orbital interactions in chemical bonds. The EDA‐NOCV has been proven very useful for clarifying bond nature and analysing the σ‐bonding and π‐bonding interactions of various main‐group compounds and transition metal complexes.[ , , ]

The EDA‐NOCV calculations were carried out with the program package ADF2019 with the zeroth‐order regular approximation (ZORA) using DFT functional BP86 and Grimme's Beck‐Johnson damping dispersion correction D3(BJ) with uncontracted Slater‐type orbitals (STOs) with TZ2P+ quality as basis functions on the BP86/def2‐SVP optimized geometries.

Results and Discussion

To determine the spin state of the (E1E2)Fe(CO)2(CNArTripp2)2 (E1E2=BF, CO, N2, CN−, NO+) compounds, the singlet and triplet states have been examined. As shown in Table 1, the singlet‐triplet splitting energy values for all five compounds are positive. Except for (N2)Fe(CO)2(CNArTripp2)2, the singlet state is more stable by over 24 kcal/mol than the triplet state. These results revealed that the ground state of five compounds is the singlet state. The calculations by using B3LYP and M06‐L also gave singlet as the ground state.

Table 1

The singlet‐triplet splitting electronic energy ES‐T and Gibbs free energy GS‐T for (E1E2)Fe(CO)2(CNArTripp2)2 (E1E2=BF, CO, N2, CN−, NO+) given in kcal/mol.

E1E2	BF	CO	N2	CN−	NO+
ES‐T	25.87	43.79	7.72[a]	32.99	28.81
GS‐T	24.27	40.87	2.93[a]	29.81	25.26

[a] The electronic energy and Gibbs free energy were calculated at the ωB97XD/def2‐SVP level because BP86/def2‐SVP met the serious convergence problem in calculating the triplet state.

Two selected optimized geometric structures of (BF)Fe‐(CO)2(CNArTripp2)2 and Fe(CO)3(CNArTripp2)2 are shown in Figure 1, and the calculated and experimental bond length and bond angle values are presented in Table 2. It can be seen from Table 2 that the calculated results are in excellent agreement with the experimental values. The optimized structures of (E1E2)Fe(CO)2(CNArTripp2)2 (E1E2=CO, N2, NO+) maintain the almost standard trigonal bipyramidal coordination geometries, as indicated by C1−Fe−C2 angles of 178.7°,177.9°and 171.8° for CO, N2, and NO+, respectively. However, the trigonal bipyramidal geometry in (BF)Fe(CO)2(CNArTripp2)2 has an obvious distortion and the C1‐Fe−C2 bond angle is only 157.5°, which may come from the better σ‐acceptor properties of BF.[ , ] The trigonal bipyramidal geometry in [(CN)Fe(CO)2(CNArTripp2)2]− has been obviously distorted, in which the bond angle of C1−Fe−C2 in is 144.3° while the CN(equatorial) −Fe−CO(equatorial) is as large as 167.6°. Besides, the E1−E2 bond in (E1E2)Fe(CO)2(CNArTripp2)2 is elongated moderately by about 0.02 Å with respect to the free E1E2 except for the distorted [(NO)Fe(CO)2(CNArTripp2)2]+.

Figure 1

Optimized geometries of (BF)Fe(CO)2(CNArTripp2)2 and Fe(CO)3(CNArTripp2)2 at the BP86/def2‐SVP level.

Table 2

Calculated and experimental results for selected bond length (L, in Å), bong angle (A, in degree), E1–E2 stretching wavenumbers (v, in cm−1) of (E1E2)Fe(CO)2(CNArTripp2)2 (E1E2=BF, CO, N2, CN−, NO+) and isolated E1E2 ligand. The experimental data are from ref. 25.

E1E2		L (Fe−E1)	L (E1−E2)	A (C1−Fe−C2)	V (E1−E2)
BF	Calc.	1.782	1.301	157.489	1401
	Expt.	1.770	1.277	160.378	1407
free BF	Calc.		1.277		1366
CO	Calc.	1.792	1.167	178.652	1972
	Expt.	1.809	1.144	178.606	1940
free CO	Calc.		1.142		2150
N2	Calc.	1.868	1.130	177.945	2180
	Expt.	1.885	1.105	179.038	2194
free N2	Calc.		1.112		2384
CN−	Calc.	1.946	1.184	144.295	2127
free CN−	Calc.		1.192		2081
NO+	Calc.	1.679	1.156	171.809	1934
free NO+	Calc.		1.077		2401

The calculated E1−E2 stretching wavenumbers are in conformity with the experimental data. Noteworthy, it has been proved that the BP86 could reproduce very well the C−O stretching mode in carbonyl complexes. Table 3 presents the energy levels and gaps of HOMO‐LUMO of five compounds. The HOMO‐LUMO gaps are in the range of 2.38–2.77 eV for (E1E2)Fe(CO)2(CNArTripp2)2 (E1E2=BF, CO, and N2), and these large values suggest their high kinetic stabilities. Meanwhile, these three compounds also exhibit similar HOMO and LUMO energy levels. In contrast, the HOMO and LUMO energy levels for [(CN)Fe(CO)2(CNArTripp2)2]− are close to zero while those for [(NO)Fe(CO)2(CNArTripp2)2]+ are much lower, which are attributed to their charged state. The small gapHOMO‐LUMO also suggests their low kinetic stability. Figure 2 depicts the HOMOs and LUMO of (BF)Fe(CO)2(CNArTripp2)2. The HOMO and HOMO‐1 mainly come from the Fe‐3d orbital and B‐2p orbital, whereas Fe‐3d orbital, B‐2p orbital, and N‐2p orbital, as well as important contributions from the Tripp2 ligand, constitute the LUMO.

Table 3

The energy levels (in eV) of HOMO and LUMO and HOMO‐LUMO gap (gapHOMO‐LUMO, in eV) of (E1E2)Fe(CO)2(CNArTripp2)2 (E1E2=BF, CO, N2, CN−, NO+).

E1E2	BF	CO	N2	CN−	NO+
LUMO	−2.49	−2.41	−2.42	+0.06	−6.63
HOMO	−5.22	−5.18	−4.80	−1.63	−7.98
gapHOMO‐LUMO	+2.73	+2.77	+2.38	+1.69	+1.35

Figure 2

Plot of the unoccupied and occupied molecular orbitals of (BF)Fe(CO)2(CNArTripp2)2 at the BP86‐D3(BJ)/TZ2P+ level.

Table 4 shows the calculated NBO partial charges, the Wiberg bond orders P of selected atoms, and Fe–E1 bond dissociation energy. (E1E2)Fe(CO)2(CNArTripp2)2 (E1E2=BF, CO, CN−) exhibit close bond dissociation energies ΔG BDE in the range from +44.4 to +76.1 kcal/mol while a smaller value is identified for (N2)Fe(CO)2(CNArTripp2)2. The (NO)+−Fe(CO)2(CNArTripp2)2 bond is much stronger with the ΔG BDE value of 147.5 kcal/mol. The ΔG BDE calculations on the five complexes also reveal that all the substitution reactions are exothermic. For all five compounds, the iron atom carries negative charges ranging from −1.04 to −1.60. The E1 is positively charged, whereas E2 possesses a negative charge except for N2. In particular, the boron atom in (BF)Fe(CO)2(CNArTripp2)2 has an obviously positive charge of +1.04, suggesting that the BF ligand has the strong ability of electron donation. The Wiberg bond orders of the Fe‐E1 bond range from 0.79 to 1.49, indicating that the Fe−E1 is a single bond. However, the Wiberg bond orders of the E1−E2 bond are quite different. For BF, it is suggested that B−F is a single bond whereas N2 and CN− shown the bond orders range from 2.70 to 2.78, revealing triple bond character. In the case of CO and NO+, the C−O and N−O bonds both are double bonds. To disclose more detailed information on the nature of the Fe–E1 bonds in those five compounds, we carried out the EDA‐NOCV calculations on (E1E2)Fe(CO)2(CNArTripp2)2 (E1E2=BF, CO, N2, CN−, NO+). Although five 10‐electron diatomic ligands usually form a dative bond with transition metals, two bonding models including the dative bond and electron‐sharing double bond have been analyzed, as shown in Scheme 1. The calculations for the dative bond use the neutral (or charged) singlet fragments E1E2 and the remaining neutral singlet fragment. In the case of the electron‐sharing double bond, the neutral (or charged) triplet fragments E1E2 and triplet fragment were employed.

Table 4

Calculated NBO partial charges q, Wiberg bond orders P, and Fe–E1 bond dissociation energy (ΔG BDE, in kcal/mol) in (E1E2)Fe(CO)2(CNArTripp2)2 (E1E2=BF, CO, N2, CN−, NO+).

E1E2	q (Fe)	q (E1)	q (E2)	P (Fe−E1)	P (E1−E2)	ΔG BDE
BF	−1.60	+1.04	−0.46	1.26	0.87	+66.5
CO	−1.40	+0.69	−0.45	1.25	2.07	+44.4
N2	−1.12	+0.06	+0.02	0.79	2.70	+19.9
CN−	−1.26	+0.25	−0.55	0.85	2.78	+76.1
NO+	−1.04	+0.40	−0.07	1.49	1.97	+147.5

Scheme 1

Schematic representation for the two bonding models of the E1–Fe bond in (E1E2)Fe(CO)2(CNArTripp2)2 (E1E2=BF, CO, N2, CN−, NO+).

Table 5 presents the numerical EDA‐NOCV results for all the five compounds (E1E2)Fe(CO)2(CNArTripp2)2 (E1E2=BF, CO, N2, CN−, NO+). The calculations revealed that the dative interactions E2E1→Fe are considered to be the best representation for the bonding situation in (E1E2)Fe(CO)2(CNArTripp2)2 because of the smaller absolute orbital interaction term ΔE orb for the dative bond with respect to the electron‐sharing double bond. The calculations suggest that the strength of the Fe‐ligand interaction ΔE int follows the order of NO+>CN−>BF>CO>N2.

Table 5

EDA‐NOCV results of the E2E1–Fe(CO)2(CNArTripp2)2 bond in (E1E2)Fe(CO)2(CNArTripp2)2 (E1E2=BF, CO, N2, CN−, NO+) at the BP86‐D3(BJ)/TZ2P+ level using different fragments. All values are in kcal/mol. The smallest ΔE orb values for each species are indicated in bold.

	BF	CO	N2	CN−	NO+
Dative bond with (singlet) and Fe(CO)2(CNArTripp2)2 (singlet)
ΔE int	−81.6	−58.3	−34.3	−93.3	−178.0
ΔE Pauli	+273.3	+165.1	+99.4	+153.8	+137.5
ΔE disp	−9.3	−8.2	−8.9	−10.9	−7.3
ΔE elstat [a]	−198.2 (57.4 %)	−119.0 (55.3 %)	−62.7 (50.3 %)	−132.1 (55.9 %)	−49.5 (16.1 %)
ΔE orb [a]	−147.4 (42.6 %)	−96.2 (44.7 %)	−62.1 (49.7 %)	−104.2 (44.1 %)	−258.8 (83.9 %)
ΔE orb(1) [b] E1→Fe σ donation	−71.1 (48.2 %)	−31.7 (33.0 %)	−21.9 (35.3 %)	−64.5 (61.9 %)	−14.8 (5.7 %)
ΔE orb(2) [b] E1←Fe π back‐donation	−37.1 (25.2 %)	−31.7 (33.0 %)	−19.7 (31.7 %)	−7.8 (7.5 %)	−131.3 (50.8 %)
ΔE orb(3) [b] E1←Fe π back‐donation	−22.0 (15.0 %)	−20.3 (21.1 %)	−12.4 (20.1 %)	−5.9 (5.7 %)	−73.0 (28.2 %)
ΔE orb(rest) [b]	−17.2 (11.6 %)	−12.5 (12.9 %)	−8.1 (12.9 %)	−26.0 (24.9 %)	−39.7 (15.3 %)

[a] The values in parentheses give the percentage contribution to the total attractive interactions ΔE elstat + ΔE orb. [b] The values in parentheses give the percentage contribution to the total orbital interactions ΔE orb.

Electron‐sharing double bond with E1E2 (triplet) and Fe(CO)2(CNArTripp2)2 (triplet)

−158.6

−209.1

−216.6

−227.0

−339.3

+181.1

+140.3

+128.1

+154.3

+169.4

−142.2 (43.0 %)

−104.1 (30.5 %)

−80.8 (24.1 %)

−108.0 (29.2 %)

−87.2 (17.4 %)

−188.2 (57.0 %)

−237.1 (69.5 %)

−255.0 (75.9 %)

−262.5 (70.8 %)

−414.2 (82.6 %)

ΔE orb(1) [b] E1−Fe σ bond

−87.1 (46.3 %)

−138.6 (58.4 %)

−177.3 (69.5 %)

−95.2 (36.3 %)

−261.2 (63.1 %)

ΔE orb(2) [b] E1−Fe π bond

−67.9 (36.1 %)

−69.8 (29.4 %)

−57.4 (22.5 %)

−139.6 (53.2 %)

−50.8 (12.3 %)

−33.2 (17.6 %)

−28.8 (12.2 %)

−20.3 (8.0 %)

−27.6 (10.5 %)

−102.2 (24.6 %)

For the three neutral ligands, it is found that the Fe‐BF interaction strength is −81.6 kcal/mol, which is larger than the cases with CO and N2. The covalent (orbital) interactions ΔE orb and the Coulomb term ΔE elstat have similar strengths for N2 while ΔE elstat is moderately stronger than ΔE orb for BF and CO. The breakdown of ΔE orb into the most important pairwise orbital interactions shows that the σ donation from BF to Fe fragment provides about 48.2 % to the total orbital interactions while two π back‐donations from the Fe fragment to BF π* contribute contributing about 40 %. The σ donation from CO (or N2) to Fe fragment is only about 33 %–35 %. However, the two π back‐donations are stronger, which provides about 52 %–54 % to the total orbital interaction. Further examination on the energy values of σ interaction and π back‐donations is stronger than that for CO and N2, which agrees with the fact that BF has better σ‐donation and π accepting properties than CO and N2. Therefore, it is reasonable to conclude that the σ donation is more important than the π back‐donations for BF−Fe interaction, whereas for CO−Fe and N2‐Fe interactions, showing the stronger π back‐donations.

Although the charged ligand NO+ and CN− have stronger interaction strengths, the contributions from the main attractive interactions, ΔE elstat and ΔE orb, to the interaction strengths are quite different. As for CN−, the electrostatic interaction ΔE elstat is moderately larger than the orbital interaction ΔE orb, in which the σ donation provides more than 60 %. However, ΔE orb plays a dominant role in the case of NO+ and the π back‐donations from Fe to NO+ contribute about 79 % to the total orbital interaction. These results agree well with the above Wiberg bond order results that the P (Fe−CN−) is only 0.85 while the P (Fe−NO+) is 1.49.

Figure 3 (a) shows the deformation density Δρ associated with the BF→Fe(CO)2(CNArTripp2)2 σ‐donation. It comes from the σ lone pair orbital of the BF(s) ligand to the formally empty d orbitals of Fe(CO)2(CNArTripp2)2 with a large eigenvalue of ν=0.98 and the stabilization energy of ΔE orb=−71.1 kcal/mol. Two BF (π*)←Fe(CO)2(CNArTripp2)2 π back‐donations are shown by the deformation density Δρ and Δρ in Figure 3 (b) and (c), both from the doubly occupied orbitals in Fe(CO)2(CNArTripp2)2 get into the formally empty π* orbitals of the BF ligand. The charge transfer along with π back‐donations is relatively smaller of 0.70 and 0.41. While the σ‐donation from the singly occupied orbital in CO(s) ligand into the formally empty d orbitals of Fe(CO)2(CNArTripp2)2, CO(σ)→Fe(CO)2(CNArTripp2)2, is −31.7 kcal/mol with the charge flow ν1 = 0.63, as depicted in Figure 3 (d). The Δρ and Δρ in Figure 3 (e) and (f) show that the CO π*(p)←Fe(CO)2(CNArTripp2)2 π back‐donations that come from the doubly occupied orbital in Fe(CO)2(CNArTripp2)2 into the formally empty π*(p) orbitals of the CO ligands. As shown in Figure S2 (g), (h) and (i), the shapes of deformation densities in (N2)Fe(CO)2(CNArTripp2)2 are analogous to Fe(CO)3(CNArTripp2)2. The deformation densities Δρ of (E1E2)Fe(CO)2(CNArTripp2)2 (E1E2=CN−, NO+) are depicted in Figure S2 (j)∼(o).

Figure 3

Plot of deformation densities Δρ (isovalue=0.004) of the pairwise orbital interactions in (E1E2)Fe(CO)2(CNArTripp2)2 (E1E2=BF, CO) at BP86‐D3(BJ)/TZ2P+ level. The charge flow of the electronic charge is red→blue. The associated orbital interaction energies ΔE orb are given in kcal/mol. The eigenvalues ν indicate the size of the charge flow.

To further explain the different contributions from σ donation and π back‐donations for different ligands, we calculated the energy levels of E1E2 and Fe(CO)2(CNArTripp2)2 fragments in their deformed configurations, as shown in Table 6. For the three neutral ligands, the HOMO energy levels of E1E2 fragment are quite different and follow the order of BF>CO>N2. However, the LUMO energy level of E1E2 fragment and HOMO and LUMO energy levels of Fe(CO)2(CNArTripp2)2 are close. Thus, the energy gap between and follows the order of BF < CO < N2, indicating that the strength of σ donation adopts the order of BF>CO>N2, in agreement with the above EDA‐NOCV results. Both HOMO and LUMO energy levels of the charged CN− fragment are higher than zero and the σ donation should be much stronger than π back‐donations, in conformity with the above results that the σ donation and π back‐donations provide about 62 % and 13.2 % to the total orbital interactions, respectively. In contrast, the lower HOMO and LUMO levels for the negatively charged NO+ suggest that the − interaction is important, agreeing with the strong π back‐donations from Fe(CO)2(CNArTripp2)2 fragment to NO+, as shown in Table 6.

Table 6

The energy levels (eV) of HOMO and LUMO of E1E2 and Fe(CO)2(CNArTripp2)2 fragments of (E1E2 = BF, CO, N2, CN−, NO+).

Ligand	E1E2 fragment			Fe(CO)2(CNArTripp2)2 fragment
	HOMO	LUMO		HOMO	LUMO
BF	−6.98	−2.35		−4.17	−3.91
CO	−9.23	−2.41		−4.18	−3.28
N2	−10.21	−2.30		−4.26	−3.22
CN−	+0.55	+6.94		−4.52	−4.27
NO+	−23.12	−15.62		−3.80	−3.04

Conclusion

The stability, electronic structures, and bonding properties of T(E1E2)Fe(CO)2(CNArTripp2)2 (E1E2=BF, CO, N2, CN−, NO+) have been studied using density functional calculations. All those molecules have a singlet ground state. NBO analysis revealed that Fe is negatively charged while E1 possesses positive charges. By employing the energy decomposition analysis, the binding nature of the E2E1−Fe(CO)2(CNArTripp2)2 bond was disclosed to be the classic dative bond E2E1→Fe(CO)2(CNArTripp2)2 rather than the electron‐sharing double bond. The bonding strength between BF and Fe(CO)2(CNArTripp2)2 is much stronger than that between CO (or N2) and Fe(CO)2(CNArTripp2)2, which comes from the better σ‐donation and π back‐donations. Thus, it can be sketched with the Lewis structure FB Fe(CO)2(CNArTripp2)2 because of the obvious π back‐donations. However, the orbital interactions in CN−→Fe(CO)2(CNArTripp2)2 and NO+→Fe(CO)2(CNArTripp2)2 mainly comes from σ‐donation and π back‐donations, respectively. The different contributions from σ donation and π donation for different ligands were well explained by using the energy levels of E1E2 and Fe(CO)2(CNArTripp2)2 fragments.

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.