The electronic structure and properties of Ni(C
6
H
4-n
F
n
)(CO)
2
(C
6
H
4
=benzyne, n=1-4) complexes have been investigated using hybrid density functional B3LYP theory. Both aromatic natures and nucleus independent chemical shift (NICS) of the benzyne rings have been analyzed. Among mono-, di-, and tri-fluorinated complexes,
3-F, 3,6-F,
and
4-H
are the most stable isomers, respectively. NICS values calculated at the several points above the ring centers are consistent with those based on the relative energies of the complexes. The atoms in molecules (AIM) analysis indicates that Ni-C bond distance is well correlated with the electron density of a ring critical point (ρ
rcp
) in all species.
INTRODUCTION
Benzyne has been the subject of many theoretical and experimental investigations
1
-
6
because it is found as an intermediate of many important organic or biochemical reactions.
3
,
7
-
8
Since the first attempt for making metal complexes using benzyne by by
Wittig
and Bickelhaupt in 1958,
9
many benzyne complexes have been successfully prepared,
10
-
13
for example, M. A. Bennett
et al
.
14
synthesized organometallic compounds, NiL
2
(C
6
H
4
)(L = PCy
3
, P
i
Pr
3
; Cy = cyclohexyl,
i
Pr = isopropyl). K. R. Deaton and M. S. Gin studied the reactions of nickel(0)-benzyne complexes with symmetrically substituted 1,3-diynes in the presence of triethylphosphine, which lead to the regioselective formation of 2,3-dialkynyl naphthalenes
15
. In the present study, the quantum chemical methods were used in order to gain a deeper insight into the structure and bonding of Ni(C
6
H
4-n
F
n
)(CO)
2
(C
6
H
4
= benzyne, n = 1 - 4) complexes and phenomena of the substituent effect in a benzyne ring.
- Computational Methods
All calculations were carried out with the Gaussian 03 suite of program.
16
Light atoms (C, H, O and F) were described by the standard 6-31+G(d) basis set.
17
Ni was described by the effective core potential (ECP) of Wadt and Hay pseudo-potential
21
with a double-ξ valance using the LANL2DZ basis set.
18
-
20
Geometry optimization was performed using Becke’s hybrid three-parameter exchange functional and the nonlocal correlation functional of Lee, Yang, and Parr (B3LYP).
22
A vibrational analysis was performed at each stationary point which corresponds to an energy minimum.
The nucleus-independent chemical shift (NICS)
23
,
24
has been defined as the absolute magnetic shielding computed at the center of a ring in a molecule. NICS(0.0), NICS(0.5), NICS(1.0), NICS(1.5) and NICS(2.0) were calculated at 0 (center), 0.5, 1.0, 1.5, and 2.0 Å above the ring, respectively.
The AIM2000 program
25
was used for the topological analysis of electron density, and the characteristics of ring critical points (RCPs) were taken into account: density at RCP (ρ(r
c
)), and its Laplacian (∇
2
ρ(r
c
)).
RESULT AND DISCUSSION
- Energy and geometry
The calculated energy values of the compounds in this work are listed along with the selective bond angles in
1.
When the compounds are classified with geometric isomers, 3-F, 3,6-F, and 4-H are the most stable isomer in each group.
Optimized structures of the compounds are displayed with bonds distances. Notice that the compounds can be classified into six groups according to their chemical natures, complexation or the number of F atoms: (benzyne), (H4), (3-F, 4-F), (5,6-F, 4,6-F, 3,6-F, 4,5-F), (3-H, 4-H), and (F4). In the multi-fluorinated compounds, 3-F, 3,6-F, and 4-H are the most stable ones which are used for the discussion of the their representative properties in this work.
Energies (Hartree), relative energies (kcal/mol), and selected bond angles (deg) for of benzyne and Ni(C6H4-nFn)(CO)2(C6H4= benzyne, n = 1-4) complexes. SeeFig. 1for their structures
Energies (Hartree), relative energies (kcal/mol), and selected bond angles (deg) for of benzyne and Ni(C6H4-nFn)(CO)2 (C6H4 = benzyne, n = 1-4) complexes. See Fig. 1 for their structures
The structures of the Ni(C
6
H
4-n
F
n
)(CO)
2
complexes were optimized by DFT calculations (
. 1
). The C-C bond lengths in the benyne-Ni complexes seem to be related with the changes in aromaticity. The C1-C6 bond distance of benzyne (1.386 Å) becomes elongated to 1.393 Å when benzyne forms a complex with Ni(CO)
2
in
H4
. In contrast, this bond tends to decrease as the number of F atoms increase in the complex if the most stable isomers of the fluorinated complexes are considered: 1.382 Å (
3-F
), 1.380 Å (
3,6-F
), 1.377 Å (
4-H
), and 1.378 Å (
F4
).
Charges of Ni, C and O atoms calculated by NBO analysis for benzyne and Ni(C6H4-nFn)(CO)2complexes.
Charges of Ni, C and O atoms calculated by NBO analysis for benzyne and Ni(C6H4-nFn)(CO)2 complexes.
The bond distances of Ni-C(benzyne) decrease, and the Ni-C(O) bonds increase, revealing that the fluorination of benzyne reduces the back-donation of electron densities of Ni to the carbonyls. This observation is supported by the increase of Ni charges in the fluorinated complexes (
2
): 0.420 (
H4
) to 0.421 (
3-F
), 0.432 (
3,6-F
), 0.432 (
4-H
), and 0.438 (
F4
). The bond angles around Ni atom also show a systematic trend that the C1-Ni-C2 bond angles increase while the C7-Ni-C8 ones decrease in the fluorinated complexes.
- Frontier orbitals
The HOMO and LUMO energies of the compounds in
3
show that all the Ni-complexes have higher HOMO and lower LUMO energies than those of benzyne, respectively. The HOMO—LUMO gap energies tend to decrease as the number of F atoms increases.
HOMO, LUMO and HOMO-LUMO gap energies for benzyne and Ni(C6H4-nFn)(CO)2(C6H4=benzyne, n=1-4) complexes. The vibration frequencies involving C≡C and carbonyls are also listed
HOMO, LUMO and HOMO-LUMO gap energies for benzyne and Ni(C6H4-nFn)(CO)2 (C6H4=benzyne, n=1-4) complexes. The vibration frequencies involving C≡C and carbonyls are also listed
- Vibration Analysis
The symmetric and asymmetric stretching modes of CO are displayed in
2
and the calculated frequencies of CO ligands, ν(CO), have been presented in
3
. These values increase in fluorinated complexes when the most stable isomers are considered.
The IR spectroscopy experimental results show that in NiL
2
(C
6
H
4
) (L=PCy
3
, Cy = cyclohexyl), the CI-C2 bond length of 1.332 Å is greater than the value of 1.29 Å typical of most ML
2
(alkyne) complexes, but the lengthening on coordination relative to the free ligand (Δr = 0.08-0.09 Å) is about the same.
11
Correspondingly, the ν(C≡C) value of ca. 1580 cm
-1
in ML
2
(C
6
H
4
) (M = Ni, Pt) complexes is less than the value of ca. 1700 cm
-1
found in ML
2
(alkyne) complexes, but the decrease relative to the free ligand is of the same order (ν(C≡C)= 400-500 cm
-1
) in both series. In the Ni(C
6
H
4-n
F
n
)(CO)
2
systems, the r(C1-C2) in complexes is greater than free benzyne (Δr
max
= 0.083 Å in F4 molceule). On the other hand, the ν(C≡C) values increase in complexes rather than free benzyne (
3
).
Stretching modes of carbonyl groups in the Ni-complexes: (a) asymmetric, (b) symmetric.
- 1H,13C-NMR
Along with the vibration analysis, we could compute
1
H-NMR chemical shifts for the set of compounds by GIAO method. The calculated chemical shifts of protons are listed in
4.
The H3 - H6 signals of
H4
are expected to appear at down-field compared with those in benzyne, indicating that the complexation makes the benzyne ring more electron-deficient. When the F atom is introduced to H4, some of protons shift to up-field although a systematic trend is not clearly detected. This is probably due to the induction effect of the electronegative F atom which attracts electrons again from Ni to the benzyne moiety.
Calculated proton chemical shifts (ppm) values for of benzyne and Ni(C6H4-nFn)(CO)2(C6H4=benzyne, n=1-4) complexes
Calculated proton chemical shifts (ppm) values for of benzyne and Ni(C6H4-nFn)(CO)2 (C6H4=benzyne, n=1-4) complexes
The benzyne
13
C chemical shifts provide some support for aromaticity. The experimental d values for C1 and C2 in NiL
2
(C
6
H
4
) (L=PCy
3
, Cy=cyclohexyl) are 145.2 ppm.
11
Also, the computed
13
C NMR chemical shifts of the Ni(C
6
H
4-n
F
n
)(CO)
2
are compatible with aromticity (
1
).
- Nucleus-independent chemical shift (NICS)
As an effort to discuss the use of NICS as a measure of aromaticity, we have calculated NICS values from the center of the ring to 2.0 Å above the plane along the z-axis of the benzyne ring. The shape of NICS profile with respect to the distance from the ring center falls into two categories. In addition, for all species, we have localized both the NICS maxima and minima, and determined the distances to the center of the ring at which they occur (
5
). For each benzyne, tri-, and tetra-fluorinated species and
H4
, the highest absolute value of NICS closes to the center of the ring. Both mono- and di-fluorinated species have a maximum about 0.5 Å to the ring center. It is possible that induced magnetic fields generated by the σ aromaticity are particularly large in the center of the ring. However, the molecular systems having π aromaticity have a minimum NICS at the certain distances from the center of the ring. There is a linear correlation between NICS(0.0) and NICS(0.5) values in all complexes: R
2
= 0.968, not shown.
NICS(0.0), NICS(0.5), NICS(1.0), NIS(1.5), and NICS(2.0) values for benzyne ring in benzyne and Ni(C6H4-nFn)(CO)2(C6H4=benzyne, n=1-4) complexes. The values in the parentheses show the distances in Å
NICS(0.0), NICS(0.5), NICS(1.0), NIS(1.5), and NICS(2.0) values for benzyne ring in benzyne and Ni(C6H4-nFn)(CO)2 (C6H4=benzyne, n=1-4) complexes. The values in the parentheses show the distances in Å
- AIM analysis
As it is difficult to separate the σ and π contributions to the electron density at the bond critical point, the ρ(r) values can be used to evaluate bond strength for different types of bonds (
6
). The different values of ρ(r) and ∇
2
ρ(r) for the Ni-C bonds evidently indicate the relative Ni-C bond strengths. This result is in agreement with the geometrical analysis, showing that the Ni-C bond of Ni(C
6
H
4
)(CO)
2
is longer than other species. On the other hand, the Ni-C bonds in all structures have positive values of ∇
2
ρ(r) which is indicative of the close shell interaction.
The value of electron density and its Laplacian estimated at bond critical point of NiC(benzyne) correlate very well with the strength of the bond, as well as with its length, since, as it is well known, both the strength and length of a bond are mutually dependent. A good relationship is present between Σρ(NiC
benzyne
) values and Σr(NiC
benzyne
) (R
2
=0.994). Similarity, Σ∇
2
ρ(NiC
benzyne
) values obeys a linear relationship (R
2
=0.941), too.
Selected AIM based parameters for (a) Ni-C(benzyne) and (b) NiC(O) bonds: Electron density (ρ), Laplacian of electron density (∇2ρ), kinetic electron energy density, G(ρ), the total electron energy density, H(ρ), potential electron energy density, V(ρ), and elipticity of the Ni(C6H4-nFn)(CO)2(C6H4=benzyne, n=1-4) complexes
Selected AIM based parameters for (a) Ni-C(benzyne) and (b) NiC(O) bonds: Electron density (ρ), Laplacian of electron density (∇2ρ), kinetic electron energy density, G(ρ), the total electron energy density, H(ρ), potential electron energy density, V(ρ), and elipticity of the Ni(C6H4-nFn)(CO)2 (C6H4=benzyne, n=1-4) complexes
The comparison of electron density in the bond critical points of Ni-C(benzyne) and Ni-C(O) shows that ρ(Ni-C
benzyne
) is smaller than ρ(Ni-C(O)). This trend is wellmatched with the results of the geometrical analysis.
The bond ellipticity is defined as ε=(λ
1
/λ
2
)-1, where |λ
1
|≥|λ
2
|. It provides a quantitative measurement of the π character of the bond. The plane of the π distribution is uniquely specifies by the direction of the axis associated with the curvature of smallest magnitude, λ
2
. The ε(NI-C) values show that the Ni-C bond in fluorinated rings has a smaller π-character in comparison with Ni(C
6
H
4
)(CO)
2
(
6
).
Further useful information on the chemical bond properties is obtainable from the total electron energy density (H(ρ)) and its components, a kinetic electron energy density (G(ρ), positive by definition) and a potential electron energy density (V(ρ), negative by definition). The following relation is known for H(ρ) and its components:
In the region of the bond critical point for the weak closed-shell inter-atomic interactions (Ni-C), the kinetic energy density dominates with the G(ρ) magnitude being slightly greater than the potential energy density |V(ρ)|, which implies that the total energy density H(ρ) is positive and closes to zero. For the strong covalent interactions (Ni-C), V(ρ) dominates over the kinetic energy density and H(λ) < 0. This usually accompany with ∇
2
ρ > 0 for H(λ) > 0, and ∇
2
ρ < 0 for H(ρ) < 0.
Electron density (ρ(3,+1)), Laplacian of electron density (∇2ρ(3,+1)) at the ring critical point (RCP) of a benzyne ring in each Ni(C6H4-nFn)(CO)2(C6H4=benzyne, n=1-4) compound
Electron density (ρ(3,+1)), Laplacian of electron density (∇2ρ(3,+1)) at the ring critical point (RCP) of a benzyne ring in each Ni(C6H4-nFn)(CO)2 (C6H4=benzyne, n=1-4) compound
At the ring critical point of benzyne, both the electron density ρ(3,+1) and Laplacian of electron density ∇
2
ρ(3,+1) have been calculated for all complexes (
7
). It is observed that these values become smaller in the fluorinated complexes when they are compared with those of benzyne and
H4
.
CONCLUSION
We investigated the structures and frontier orbitals of the Ni(C
6
H
4-n
F
n
)(CO)
2
(C
6
H
4
=benzyne, n=1-4) complexes. The results suggest that
3-F, 3,6-F, 4-H
isomers are most stable among the mono-, di-, tri-fluorinated complexes, respectively. Using the stable ones, both vibration and
1
H-NMR shifts were calculated and analyzed. The NICS calculations confirmed the aromaticity in the benzyne rings of the compounds. Using the analyses of both electron densities and energy densities, we could explain the characters of the Ni-C bonds in complexes.
Proft F. D.
,
Schleyer P. V. R.
,
Lenthe J. H. V.
,
Stahl F.
,
Geerlings P,
2002
Chem. Eur. J.
8
3402 -
Liu H.
,
Yang S.
,
Balteanu I.
,
Balaj O. P.
,
Fox-Beyer B. S.
,
Beyer M. K.
,
Bondybey V. E.
2004
Rapid Commun. Mass Spectrom
18
1479 -
DOI : 10.1002/rcm.1512
Andino J. G.
,
Kilgore U. J.
,
Pink M.
,
Ozarowski A.
,
Krzystek J.
,
Telser J.
,
Baik M.-H
,
Mindiola D. J.
2010
Chem. Sci.
Advance Article.
Hughes R. P.
,
Laritchev R. B.
,
Williamson A.
,
Incarvito C. D.
,
Zakharov L. N.
,
Rheingold A. L.
2002
Organometallics
21
4873 -
DOI : 10.1021/om0204787
Retbøll M.
,
Edwards A. J.
,
Rae A. D.
,
Willis A. C.
,
Bennett M. A.
,
Wenger E.
2002
J. Am. Chem. Soc.
124
8348 -
DOI : 10.1021/ja0264091
Deaton K. R.
,
Gin M. S.
2003
Organometallic
5
2477 -
Frisch M. J.
,
Trucks G. W.
,
Schlegel H. B.
,
Scuseria G. E.
,
Robb M. A.
,
Cheeseman J. R.
,
Montgomery J. A.
,
Vreven T.
,
Kudin K. N.
,
Burant J. C.
,
Millam J. M.
,
Iyengar S. S.
,
Tomasi J.
,
Barone V.
,
Mennucci B.
,
Cossi M.
,
Scalmani G.
,
Rega N.
,
Petersson G. A.
,
Nakatsuji H.
,
Hada M.
,
Ehara M.
,
Toyota K.
,
Fukuda R.
,
Hasegawa J.
,
Ishida M.
,
Nakajima T.
,
Honda Y.
,
Kitao O.
,
Nakai H.
,
Klene M.
,
Li X.
,
Knox J. E.
,
Hratchian H. P.
,
Cross J. B.
,
Adamo C.
,
Jaramillo J.
,
Gomperts R.
,
Stratmann R. E.
,
Yazyev O.
,
Austin A. J.
,
Cammi R.
,
Pomelli C.
,
Ochterski J. W.
,
Ayala P. Y.
,
Morokuma K.
,
Voth G. A.
,
Salvador P.
,
Dannenberg J. J.
,
Zakrzewski V. G.
,
Dapprich S.
,
Daniels A. D.
,
Strain M. C.
,
Farkas O.
,
Malick D. K.
,
Rabuck A. D.
,
Raghavachari K.
,
Foresman J. B.
,
Ortiz J. V.
,
Cui Q.
,
Baboul A. G.
,
Clifford S.
,
Cioslowski J.
,
Stefanov B. B.
,
Liu G.
,
Liashenko A.
,
Piskorz P.
,
Komaromi I.
,
Martin R. L.
,
Fox D. J.
,
Keith T.
,
Al-Laham M. A.
,
Peng C. Y.
,
Nanayakkara A.
,
Challacombe M.
,
Gill P. M. W.
,
Johnson B.
,
Chen W.
,
Wong M. W.
,
Gonzalez C.
,
Pople J. A.
2003
Revision B.03 ed.
Gaussian, Inc.
Pittsburgh PA
Chen Z.
,
Wannere C. S.
,
Corminboeuf C.
,
Puchta R.
,
Schleyer P. V. R.
2005
Chem. Rev.
105
3842 -
DOI : 10.1021/cr030088+
Schleyer P. V. R.
,
Maerker C.
,
Dransfeld A.
,
Jiao H.
,
Hommes N. J. R. v. E
1996
J. Am. Chem. Soc.
118
6317 -
DOI : 10.1021/ja960582d
Bader R. F. W
2000
ver 2.0, ed.
McMaster University
Hamilton