S. J. Sabounchei*a,
M. Ahmadianpoora,
A. Yousefia,
A. Hashemia,
M. Bayata,
A. Sedghia,
F. Akhlaghi Bagherjerib and
R. W. Gableb
aFaculty of Chemistry, Bu-Ali Sina University, Hamedan, 65174, Iran. E-mail: jsabounchei@yahoo.co.uk; Fax: +98 8138273231; Tel: +98 8138272072
bSchool of Chemistry, University of Melbourne, Victoria 3010, Australia
First published on 11th March 2016
The reaction of sulfur ylides (L) SMe2C(H)C(O)R (R = 4-nitrophenyl, phenyl, and 3-nitrophenyl) with the dichloro(1,5-cyclooctadiene)palladium(II) complex, [PdCl2(cod)], in a 2:
1 ratio gives the new Pd(II) complexes of type cis- and trans-[PdCl2(SMe2C(H)C(O)R)2] (R = phenyl (1), 3-nitrophenyl (2), and 4-nitrophenyl (3)). Characterization of the obtained complexes was performed by elemental analysis, IR, 1H, 13C NMR and mass spectroscopies. Also, the structure of complex trans-[PdCl2(SMe2C(H)C(O)C6H4-p-NO2)2] (3) was characterized by single crystal X-ray analysis. The X-ray crystallography results reveal that the structure of complex 3 contains two Cα-coordinated sulfur ylide ligands in trans geometry. These air/moisture stable complexes were also employed as efficient catalysts for the Mizoroki–Heck cross-coupling reaction of several aryl halides with olefins. The coupled products of these reactions were obtained in good to excellent yields and purity, short reaction times and low catalyst loading. Also a theoretical study on the structure of the complexes 1–3 have been investigated at the BP86/def2-SVP level of theory. The bonding situation between the [PdCl2] and L2 fragments [(Ylide)2] in the [L2PdCL2] complexes, were carried out by NBO and energy-decomposition analysis (EDA), as well as their natural orbitals for chemical valence variation (EDA-NOCV). The results confirmed that the contribution of the electrostatic interactions in the Pd–C bond in the complexes is about 50%.
Coordination modes of these compounds in several metal complexes have been studied and it was found that these ylides can coordinate to metal ions in three modes as depicted in Chart 1: mononuclear, dinuclear and polymeric structure modes. When a metal ion such as Pd(II) interacts only with the Cα atoms of two ylidic groups, these ligands can be placed in the cis or trans positions to form the mononuclear structure.8 Dinuclear structure was also observed in case of Hg(II) halides, in which two ylides were coordinate from the Cα atoms to dimeric Hg2X4 (X = Cl, Br, I) unit.9,10 Finally, the polymeric structure was seen when metal ions have bridging ligands. For example, it was reported a polymeric structure of sulfur ylide with AgNO3 by our group and the X-ray analysis shows a structure with a metallic skeleton based on the C-coordinated sulfur ylide and bridged nitrate group.11
It was reported that some of phosphorus ylide palladium complexes have useful applications in catalytic C–C coupling reactions.12–15 Although these complexes have some advantages such as high efficiency, low catalyst loading and moisture/air stability, but most of them are expensive and have environmental considerations. In the past years, phosphine free catalytic reactions were developed to eliminate or reduce costs, operational hazards and environmental pollution.16–19 Among these, there are some catalytic reactions that sulfur containing ligands and their Pd(II) complexes were employed as catalyst in the reaction.20–22
Heck reaction is one of the most widely adopted methods for the construction of carbon–carbon bonds in modern organic chemistry and it is commonly carried out in the presence of a palladium complex.23–28 For example, in the reported work by the group of Gruber et al. palladium catalysts containing modified sulfur ligands could activate the coupling reactions with less reactive aryl bromides and aryl chlorides as well as nonactivated alkenes.29 It appears that the release of Pd(0) species in Heck reactions is depended on the ligands around the palladium complex. Thus, organochalcogen ligands such as sulfur ylides, show rivalry with their respective phosphorus ylide analogues for releasing of Pd(0) species in the Heck C–C coupling reaction.
Although it has been published the synthesis of palladium complexes with sulfonium ylides by Bravo et al.,8 but only the trans structure of complexes was reported.
In this paper, we report the synthesis and theoretical study of the both cis- and trans-[Pd((Me)2SCHC(O)R)2] (R = phenyl, 3-nitrophenyl and 4-nitrophenyl) complexes (1–3). In addition, the X-ray crystal structure of complex 3 and further spectroscopic characterization of complexes was performed. Finally, Mizoroki–Heck reactions of various aryl halides with styrene and ethyl acrylate using these complexes as efficient catalysts were performed, giving good to excellent yields of products under aerobic conditions.
In the far IR spectra of complex 3, two Pd–Cl stretching band was observed in the 327 and 258 cm−1 that is in accord with both cis- and trans-dichlorobis(sulfur ylide) structure in the solid state. IR spectra of complexes 1–3 show a significantly frequency shift of ν(CO) than those of the related sulfur ylides. An increase in ν(CO), indicates the coordination of the ylide was occurred through carbon; whereas decrease of ν(CO) was expected when the ylide coordinates to metal from oxygen atom. The observed absorption bands about 1514–1540 cm−1 in ylides Y1–Y3, were shifted to 1629–1647 cm−1 in complexes 1–3, indicating the coordination of ylide through the Cα atoms. The ν(S+–C) band of the ylides is also sensitive to the coordination. It was seen that a lower frequency shift was occurred to about 836–860 cm−1 in comparison to the parent ylides (810–823 cm−1), that is in consistent with some removal of electron density in the S–C bonds. However, infrared spectroscopy could not differentiates between geometric isomers (cis and trans) of each complexes.
The 1H NMR signal for the SCH group of the complexes 1–3 shows a downfield shift compared to those of the free ylide. The coexistence of the cis and trans isomers is confirmed by the 1H NMR spectra. Complexes 1–3 in dimethyl sulfoxide displayed two doublets for SMe2 and two doublet of methine signals. The higher field pair of SMe2 and methine proton signals may be assigned to the trans isomer. Since methinic carbons are chiral, two (RR/SS or RS/SR) are possible for each cis- and trans-structure of these complexes. In addition, integrated intensity of the SMe2 signals revealed that the ratio of the cis to trans isomer was about 1:
2. The fact that PdCl2(ylide)2 has a larger portion of the trans isomer would be mainly attributed to a smaller steric repulsion between ylide ligands.
In the 13C NMR spectra, a similar downfield shift for characteristic groups SCH and carbonyl was observed. Coordination of ylide through the Cα causes more electron density on this atom, led to the appearance of related 13C signal in higher field. In addition, presence of two geometric isomers in these complexes was also confirmed by 13C NMR. Observation of four CH3 signals for the two SMe2 groups and two methinic carbon signals indicates the presence of such geometrical and optical isomers. However, since carbonyl carbons do not have a hydrogen attached to them, they appear less intense than other carbons in the spectra and the signals were not seen as two separate peaks.
Empirical formula | C20H22Cl2N2O6PdS2 |
Formula weight | 627.81 |
Temperature (K) | 130.01(10) |
Crystal system | Triclinic |
Space group | P![]() |
a (Å) | 7.1934(5) |
b (Å) | 7.5432(5) |
c (Å) | 11.4096(8) |
α (°) | 94.723(6) |
β (°) | 103.384(6) |
γ (°) | 91.985(6) |
V (Å3) | 599.32(7) |
Z | 1 |
Dcalc g cm−3 | 1.739 |
Absorption coefficient (mm−1) | 10.279 |
F(000) | 316.0 |
Crystal size (mm3) | 0.1272 × 0.0933 × 0.0452 |
θ range for data collection (°) | CuKα (λ = 1.54184) |
Reflections collected | 8 to 154.03 |
Index ranges | −9 ≤ h ≤ 7, −9 ≤ k ≤ 8, −11 ≤ l ≤ 14 |
Independent reflections | 5152 |
Data/restraints/parameters | 2488 [Rint = 0.0419, Rsigma = 0.0562] |
Goodness-of-fit on F2 | 2488/0/153 |
Final R indices [I > 2σ(I)] | 1.090 |
Final R indices (all data) | R1 = 0.0475, wR2 = 0.1268 |
Largest diff. peak and hole (e Å−3) | R1 = 0.0532, wR2 = 0.1331, 1.96/−1.09 |
The Pd atom is in a square planar arrangement, coordinated by two ligands and two chlorides, both ligands and both chlorides are in a trans configuration. The ylide aromatic ring is almost perpendicular to the PdCl2 C2 plane. The nitro group and the keto group are oriented slightly out of the aromatic ring, the dihedral angles being 13.7(3)° and 7.2(7)° respectively. There are π⋯π interactions between the aromatic rings of two adjacent molecules, the interaction being such that C5 is superimposed over C7 (C5–C7′: 3.550(8) Å, C4–C5–C7′: 86.3(3)°) and C6 is superimposed over the center of the other aromatic ring. The molecules are linked together by weak C–H⋯X interactions into a 3D network.
![]() | ||
Fig. 2 Optimized structures of trans-like (a, c, and e) and cis-like (b, d, and f) of complexes 1–3 as well as the energy of cis-like isomer relative to trans-like isomer at BP86/SVP level of theory. |
The optimized structures of both trans like and cis like complexes 1–3 at the BP86/def2-SVP level of theory are shown in Fig. 2, and trends for the variation of the corresponding bond lengths and bond angles extracted from optimizing geometries and are given in Table S3.† The calculated bond lengths and bond angles of complex 3 are in good agreement with the experimental values (Table S3†).
The computed Pd–C and Pd–Cl bond distances are about 0.1 Å longer than the X-ray values. Also the result showed that the R substituent has an insignificant effect on both the Pd–C and Pd–Cl bond distances and C–Pd–C, C–Pd–Cl and Cl–Pd–Cl bond angles in the complexes 1–3 (Table S3†).
For studying the nature of metal–ligand bonds and/or possible interactions in these complexes, the natural bond orbital (NBO) analysis was performed and used to estimate the delocalization of electron density between occupied Lewis-type orbitals and formally unoccupied non-Lewis NBOs (antibonding or Rydberg), which corresponds to stabilizing donor–acceptor interaction.45 The values of the partial charge on the Pd, C and Cl atoms involved in interaction between the PdCl2 and L2 in of the complexes 1–3 were investigated at BP86/def2-SVP level of theory (Table 2). The results show that the Pd and C atoms of the Pd–C bonds carry positive and negative partial charges and this confirm that there is σ donation between the lone pair of C atoms and empty orbitals of Pd in Pd–C bonds (Table 2).
Compound | WBIs | NPA | |||
---|---|---|---|---|---|
Pd–C | P–Cl | Pd | C | Cl | |
1 | 0.406 | 0.6453 | 0.067 | −0.67 | −0.397 |
2 | 0.410 | 0.6751 | 0.034 | −0.634 | −0.390 |
3 | 0.402 | 0.6488 | 0.067 | −0.671 | −0.390 |
Also the results show that the R substituent has an insignificant effect on the values of WBIs for Pd–C bonds (Table 2). A second-order perturbation theory analysis of the Fock matrix was carried out to evaluate the donor–acceptor interactions on the base of NBO analysis. The values of donor–acceptor interactions in the compounds are given in Table S4.† In these compounds, the important donor–acceptor interactions in the case of Pd–C bonds in complexes 1–3 were carry out from the σ(Pd–C) and σ(Pd–Cl) as donors to σ*(Pd–C), σ*(Pd–Cl), σ*(S–C) and σ*(O–C) as acceptors (Table S4†).
The calculated energies for HOMO and LUMO of 1, 2 and 3 complexes are also given in Table 3. The result showed that the calculated energy gap (ΔE) between the latter orbitals for the complex 3 containing NO2 as a substituent in para position is smaller than those for other complexes. Also, the latter complex is softer than the other complexes described here. The HOMO and LUMO for all three mononuclear complexes are illustrated in Fig. S1 (see ESI†).
Compound | LUMO (hartree) | HOMO (hartree) | η (eV) | ΔE (eV) |
---|---|---|---|---|
1 | −0.0989 | −0.1700 | 0.97 | 1.94 |
2 | −0.1261 | −0.1836 | 0.78 | 1.56 |
3 | −0.1364 | −0.1886 | 0.71 | 1.42 |
To better understand the nature of the Pd–C bonding between the L2 fragments [(Ylide)2] and [PdCl2] fragment (see Scheme 1) in 1–3 complexes, we carried out quantum chemical calculations in the terms of energy-decomposition analysis at the BP86/TZ2P(ZORA)//BP86/def2-SVP with C1 symmetry with the program package ADF2009.01. The ΔEint between two investigated fragments in complexes 1–3 (see Scheme 1) are about 149.4, 145.83 and 139.26 kcal mol−1 respectively (Table 4).
Compound | |||
---|---|---|---|
1 | 2 | 3 | |
ΔEint | −149.38 | −145.83 | −139.26 |
ΔEPauli | 384.16 | 378.20 | 367.60 |
ΔEelstat | −264.61(49.6%) | −261.42(49.9%) | −257.42(50.8%) |
ΔEorb | −246.98(46.3%) | −240.51(45.9%) | −227.31(44.8%) |
ΔEdisp | −21.94(4.1%) | −22.09(4.2%) | −22.13(4.3%) |
The results also show that there is a decreasing trend in the interaction energies in the complexes as following 1 > 2 > 3. Thus the interaction energies of complex 1 with H as substituent is greater than complex 2 with NO2 as substituent in meta position and both of them are greater than complex 3 with NO2 as substituent in para position. The breakdown of the ΔEint values into the Pauli repulsion ΔEPauli and the three attractive components shows that roughly 50–51% comes from the electrostatic attraction ΔEelstat while ∼45–46% comes from the orbital term ΔEorb and the remaining concern to ΔEdispersion. Thus the values of ΔEelstat the covalent bonding between the [PdCl2] fragment and L2 fragments [(Ylide)2] in 1–3 complexes becomes visible by the calculated deformation densities Δρ. Fig. 3 show the important deformation densities Δρ and the associated energy values which provide about 94% of the overall orbital interactions for the complex 1. Visual inspection of Fig. 3a indicates the Δρ1, which show the mainly contribution to ΔEorb and comes from the Cl2Pd←L2 σ donation. Fig. 3c displays the deformation density Δρ3 which come from charge transfer from L2 to [PdCl2] fragment. Also Fig. 3b, d and e display the deformation densities Δρ2, Δρ4 and Δρ5 which comes from the [PdCl2] fragment to Y2 ligands σ-back donation.
Using the optimized reaction conditions (Cs2CO3 (1.5 mmol), DMF (2 ml), catalyst 3 (0.005 mmol), reflux temperature) we explored the general applicability of Pd–sulfur ylide complex 3 with different olefins and aryl halides containing electron withdrawing to electron donating substituents. It is well known that the catalytic activity was extensively depends on the type of aryl halides substituents. Aryl halides with electron-withdrawing groups on the aryl ring show more reactivity, while electron donating groups make the aryl halides relatively inactive.46 As expected, the Mizoroki–Heck cross-coupling reaction of aryl iodides with olefins gave better yields of the desired product than that of aryl bromides and aryl chlorides (Table 5). Both electron-neutral and electron-rich aryl iodides were converted efficiently to the desired coupled products in excellent yields (Table 5, entries 1, 4 and 7). The high reactivity of aromatic iodides is attributed to their small C–I bond dissociation energy and easy oxidative addition to Pd(0) species.47 As shown in Table 5, reaction of activated aryl bromides such as 4-bromonitrobenzene, 4-bromobenzaldeyde and 4-bromoacetophenone with styrene and ethyl acrylate afford the corresponding coupled products in high yields (Table 5, entries 9, 10, 11, 14 and 17).
Entry | Arylhalide | Olefin | Product | Yieldc (%) | Time (h) |
---|---|---|---|---|---|
a Reaction conditions for Mizoroki–Heck cross-coupling reaction: arylhalide (0.5 mmol), olefin (0.75 mmol), Cs2CO3 (1.5 mmol), DMF (2 ml), catalyst 3 (0.005 mmol), in the air.b TBAB (1 mmol) was added.c All are isolated yields. | |||||
1 | Iodobenzene | Styrene | 1 g | 92 | 4 |
2 | Bromobenzene | Styrene | 1 g | 89 | 5 |
3 | Chlorobenzene | Styrene | 1 g | 60 | 6 |
86b | |||||
4 | 4-Iodotoluene | Styrene | 2 g | 88 | 5 |
5 | 4-Bromotoluene | Styrene | 2 g | 85 | 5.5 |
6 | 4-Chlorotoluene | Styrene | 2 g | 58 | 7 |
83b | |||||
7 | 4-Iodoanisole | Styrene | 3 g | 89 | 5 |
8 | 4-Bromoanisole | Styrene | 3 g | 88 | 5 |
9 | 4-Bromonitrobenzene | Styrene | 4 g | 90 | 4 |
10 | 4-Bromobenzaldehyde | Styrene | 5 g | 82 | 4.5 |
11 | 4-Bromoacetophenone | Styrene | 6 g | 86 | 5 |
12 | 1-Bromonaphthalene | Styrene | 7 g | 81 | 16 |
13 | Bromobenzene | Ethyl acrylate | 8 g | 92 | 3.5 |
14 | 4-Bromonitrobenzene | Ethyl acrylate | 9 g | 93 | 3 |
15 | 4-Bromotoluene | Ethyl acrylate | 10 g | 90 | 4 |
16 | 1-Bromonaphthalene | Ethyl acrylate | 11 g | 85 | 12 |
17 | 4-Bromobenzaldehyde | Ethyl acrylate | 12 g | 88 | 4 |
Whereas, deactivated aryl bromides such as 4-bromoanisole and 4-bromotoluene gave lower yields (Table 5, entries 5 and 8). Also, the electronically neutral bromobenzene produced good amounts of desired coupling products (Table 5, entry 2). Additionally, trying to apply 1-bromonaphthalene as efficient aryl bromide resulted eligible yield of coupling product (Table 5, entries 12 and 16). The usually unreactive aryl chlorides were further explored as the coupling partner with styrene using the complex 3 catalyst system. Although, a number of reported catalysts in literature need to be used in high loadings, and showed little or no activity with aryl chloride substrates,48 due to the economical aspect of aryl chlorides and the challenge to activate the C–Cl bond, there was much attention from the academic communities to couple aryl chlorides with olefins.49 Both electron rich and electron neutral 4-chlorotoluene and chlorobenzene show the low reactivity in coupling reactions, but we achieved more yield by using of TBAB as additive (Table 5, entries 3 and 6). TBAB is proposed to increase the solubility of the inorganic base Cs2CO3 in the organic reaction medium and strengthened the stabilizations of Pd(0) species.50
To extend the scope of our work, we next investigated the coupling reaction of various electron rich and deficient aryl halides with activated ethyl acrylate olefin. Palladium complex 3 exhibited a better catalytic activity for the Mizoroki–Heck reaction of activated ethyl acrylate than that of inactive styrene under same reaction conditions. As expected, electron-withdrawing substituent on olefin gives the desired coupling products in more purity and less reaction time (Table 5, entries 13–17).
To evaluate the presence of Pd(0) nanoparticles (NPs) in catalytic system, we performed the mercury drop test,51 was first introduced by Whitesides. Since mercury leads to amalgamation of the surface of a heterogeneous catalyst, and in contrast, Hg(0) is not expected to have a poisoning effect on homogeneous palladium complexes, the results of mentioned test help us to determine the nature of the active species. Addition of 100 equiv. Hg(0) to the reaction solution at t = 0, was not affected the conversion of the reaction, which suggests that the catalyst is homogeneous in nature.
Also, we analyzed an aliquot of the reaction mixture by SEM analysis. After completion of the reaction, 0.5 ml of the solution was separated and evaporated under reduced pressure. The crude residue was washed with water to remove base and then analyzed by SEM. As expected, the SEM images does not showed presence of any Pd nanoparticles in recovered catalyst (see ESI†).
Footnote |
† Electronic supplementary information (ESI) available: Physical measurements and selected 13C and 1H NMR spectra of some compounds. CCDC 1438729 contains the supplementary crystallographic data for the complex 3. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c6ra01390b |
This journal is © The Royal Society of Chemistry 2016 |