Sustainable catalysis: rational Pd loading on MIL-101Cr-NH2 for more efficient and recyclable Suzuki-Miyaura reactions.

Palladium nanoparticles have been immobilized into an amino-functionalized metal-organic framework (MOF), MIL-101Cr-NH2, to form Pd@MIL-101Cr-NH2. Four materials with different loadings of palladium have been prepared (denoted as 4-, 8-, 12-, and 16 wt%Pd@MIL-101Cr-NH2). The effects of catalyst loading and the size and distribution of the Pd nanoparticles on the catalytic performance have been studied. The catalysts were characterized by using scanning electron microscopy (SEM), transmission electron microscopy (TEM), Fourier-transform infrared (FTIR) spectroscopy, powder X-ray diffraction (PXRD), N2-sorption isotherms, elemental analysis, and thermogravimetric analysis (TGA). To better characterize the palladium nanoparticles and their distribution in MIL-101Cr-NH2, electron tomography was employed to reconstruct the 3D volume of 8 wt%Pd@MIL-101Cr-NH2 particles. The pair distribution functions (PDFs) of the samples were extracted from total scattering experiments using high-energy X-rays (60 keV). The catalytic activity of the four MOF materials with different loadings of palladium nanoparticles was studied in the Suzuki-Miyaura cross-coupling reaction. The best catalytic performance was obtained with the MOF that contained 8 wt% palladium nanoparticles. The metallic palladium nanoparticles were homogeneously distributed, with an average size of 2.6 nm. Excellent yields were obtained for a wide scope of substrates under remarkably mild conditions (water, aerobic conditions, room temperature, catalyst loading as low as 0.15 mol%). The material can be recycled at least 10 times without alteration of its catalytic properties.

Introduction

Because of increasing concern regarding environmental issues, there is a great demand for the development of efficient sustainable technologies. In the area of organic synthesis, the development of efficient heterogeneous catalytic procedures[1] has become of enormous importance, as shown by the number of reports in this field in recent years.[2] For the industrial syntheses of building blocks and pharmaceutically relevant compounds, the elaboration of environmentally benign and cost-effective methods for carbon–carbon or carbon–heteroatom bond-formating reactions is of paramount importance. Within this family of reactions, the Suzuki–Miyaura cross-coupling reaction is one of the most robust and versatile methods for the formation of the carbon backbone of a target molecule, which is of obvious relevance for further progress in various branches of this field. A number of challenges remain to be overcome to achieve high efficiency under mild reaction conditions (e.g., ambient temperature) in environmentally friendly solvents. In addition, the use of cheaper and less-reactive substrates (e.g., aryl chlorides) also remains a challenge. Typical catalysts for these transformations can sometimes become very expensive, owing to the necessity of using transition-metal complexes in combination with different types of ligands. Thus, their high cost and environmental impact are major drawbacks for their large-scale application.[3]

Metal–organic frameworks (MOFs), also known as coordination polymers, have gained wide interest over the last decade. They have large surface area and are highly functionalizable.[4] MOFs have shown great potential in gas uptake and separation,[5] as well as in various biomedical applications,[6] such as drug encapsulation and delivery. In addition, MOFs can be used as catalyst supports that allow for the recycling of metallic species and also prevent the contamination of the final products.[4] However, for a MOF to be a suitable catalytic support it must be stable under the reaction conditions and allow no leaching of the active species, whilst still allowing full accessibility of the starting materials to the reaction site. However, this requirement is not always a simple task. Gratifyingly, by a thoughtful design of the linkers[7] or subsequent postsynthetic modifications,[8] it is possible to carefully tailor the inside of the pores so that they can better accommodate the catalyst and prevent it from leaching. One of the materials that meets all of the criteria mentioned above is MIL-101Cr [(Cr3F(H2O)2O[(O2C)-C6H4-(CO2)]3⋅nH2O)], which was first reported in 2005 by Férey and co-workers.[9] MIL-101Cr has a very high specific surface area (BET, about 4100 m2 g−1) and contains two types of cages with diameters of 29 and 34 Å, which have pore apertures of 12 and 16 Å, respectively. MIL-101Cr is stable in water,[10] even under very acidic conditions, and can tolerate temperatures as high as 275 °C in air. These unique features of MIL-101Cr have attracted enormous attention from researchers worldwide to explore its potential in various applications, including catalysis.[11]

Supporting metallic nanoparticles inside a porous matrix is one of the three possible ways of exploiting MOFs in catalysis, along with supporting single-site metal complexes[12] or by making use of the structural metal atoms in the framework as catalytically active sites.[13, 14] The main advantage of MOFs is that, by controlling the pore size, one can limit the growth of the nanoparticles, thereby obtaining small and well-distributed nanoparticles and avoiding agglomeration. This result is achieved by the infiltration of a metallic precursor into the framework, followed by subsequent reduction, thus leading to the direct formation of active nanoparticles in a controlled manner inside the pores. Recent reviews by Dhakshinamoorthy and Garcia summarized the use of metal nanoparticles that were immobilized inside MOFs as heterogeneous catalysts for various reactions, including condensation, hydrogenation, carbon–carbon coupling, and alcohol-oxidation reactions, as well as the synthesis of MeOH.[15] Palladium nanoparticles have shown exceptional ability to catalyze a rich variety of transformations.[16] As such, great efforts have been made to immobilize palladium nanoparticles inside the robust MIL-101Cr to form Pd@MIL-101Cr as a heterogeneous catalyst, which has been used for Suzuki coupling reactions[17, 18] or other palladium-catalyzed reactions, such as the synthesis of indoles through Sonogashira cross-coupling followed by a cycloaddition reaction,[19] the aqueous hydrogenation of phenols into cyclohexanones,[20] the synthesis of methyl isobutyl ketone from acetone,[21] the hydrogenation of ketones,[22] and even tandem processes, as cleverly illustrated by Corma and co-workers.[23]

As a proof of concept, in 2010, Jiang and co-workers published the first report[17] of a Suzuki coupling reaction of aryl chlorides in water by using Pd@MIL-101Cr at 80 °C with NaOMe as a base (1.5 equiv) and tetrabutylammonium bromide (0.6 equiv) as a co-catalyst under an inert atmosphere. From this starting point, other groups have been working on the elaboration of this method to maintain high conversions whilst trying to perform the reactions under aerobic conditions at lower temperatures.[18] Significant progress in the field of functionalized MOFs was reported in 2011 by Stock and co-workers, who described, for the first time, a rapid and facile method[24] for the postsynthetic functionalization of MIL-101Cr by covalently grafting amino groups onto the terephthalate linkers. Bearing in mind the affinity of palladium for nitrogen atoms, this modification constituted a very important breakthrough, because it minimized the leaching of metallic species. Indeed, Cao and co-workers followed this idea and presented evidence of an increased stabilization effect of the nanoparticles in amino-functionalized MIL-53Al-NH2 compared to MIL-53Al.[18, 25]

Over the past few years, our groups have undertaken on-going research into the development of new and efficient transition-metal-catalyzed processes for the formation of carbon–carbon (C–C) and carbon–heteroatom (C–X) bonds[26] and the use of MOFs in heterogeneous catalysis.[12, 14] During our investigations, we became interested in understanding the parameters that affect the efficiency of Suzuki–Miyaura cross-coupling reactions using Pd nanoparticles that are immobilized inside MIL-101Cr-NH2. In particular, we were intrigued by the fact that the reported palladium loadings in the MOFs were often low, about 1–3 wt %,[16, 19, 21, 23a, 25] with only a few exceptions with higher Pd loadings (5–6 wt %).[20] Herein, we report our investigations into the potential and limitations of MIL-101Cr-NH2 as a catalyst support of Pd nanoparticles and its application in Suzuki–Miyaura cross-coupling reactions.[27] We investigated how the loading of Pd into MIL-101Cr-NH2 affected the stability of the material and its performance as a catalyst and we determined the optimal Pd loading into the framework. Moreover, the leaching of palladium under different conditions has been investigated. TEM and total scattering techniques were employed to study the Pd nanoparticles in the MIL-101Cr-NH2 crystals. The catalyst with the optimized amount of Pd (8 wt % loaded) in MIL-101Cr-NH2 was subjected to a systematic screening of the reaction parameters, thereby leading to efficient Suzuki–Miyaura coupling reactions, with excellent yields for a large substrate scope, including challenging heteroaromatic halides and bulky substrates, under environmentally friendly conditions (water or water/EtOH, ambient temperature, K2CO3). Thus, herein, we show the importance of tuning the loading of the nanoparticles in the material and of optimization of the reaction conditions for efficient and recyclable catalysis by using MOFs, which demonstrates their great potential as supports of metal nanoparticles for practical use and large-scale applications.

Results and Discussion

Synthesis and characterization of Pd@MIL-101Cr-NH: MIL-101Cr was synthesized according to a literature procedure.[9] The amino functionality was introduced into MIL-101Cr in two steps, as reported by Stock and co-workers.[24] First, nitration of the phenyl rings was achieved by the treatment of MIL-101Cr with HNO3/H2SO4. In the second step, reduction of the nitro groups by using SnCl2⋅2 H2O afforded Pd@MIL-101Cr-NH2. The quantitative functionalization of MIL-101Cr was confirmed by both elemental analysis and by 1H NMR spectroscopy after the digestion of MIL-101Cr-NO2 and MIL-101Cr-NH2 (see the Supporting Information, Figure S1, and the Experimental Section). Palladium loadings of 4, 8, 12, and 16 wt % were prepared by the impregnation of MIL-101Cr-NH2 with the corresponding amount of [PdCl2(MeCN)2] in dry CH2Cl2 overnight, thus yielding PdII@MIL-101Cr-NH2, in which PdII was coordinated to the amino groups. Subsequent reduction into Pd0 by using NaBH4 afforded nanosized Pd particles inside the MIL-101Cr-NH2 crystals, Pd0@MIL-101Cr-NH2.

The synthesized material was washed with water and EtOH to eliminate all of the residues. The purity of the four final materials (4-, 8-, 12-, and 16 wt %Pd@MIL-101Cr-NH2) was checked by energy-dispersive spectroscopy (EDS) that was connected to a SEM, elemental analysis, and by inductively coupled plasma-optical emission spectrometry (ICP-OES; Supporting Information, Table S2). The actual Pd contents, as determined by ICP-OES, for the samples of 4-, 8-, 12-, and 16 wt %Pd@MIL-101Cr-NH2 were 4.38, 8.06, 11.29, and 15.07 wt %, respectively. Importantly, no traces of either Sn from SnCl2 or Na from the reducing reagent were found in any of the Pd@MIL-101Cr-NH2 samples.

The functionalization of the terephthalate linker first with nitro groups and then with amino groups was also followed by FTIR spectroscopy. The FTIR spectra and vibration frequencies of our synthesized samples, as shown in the Supporting Information, Figures S2 and S3 and Table S3, corresponded well with literature values.[24] Strong bands at 1540 and 1357 cm−1, which corresponded to N–O stretches, were observed in the spectrum of MIL-101Cr-NO2 (Supporting Information, Figure S3). These bands disappeared after reduction of the nitro groups and a clear N–H bending frequency band at 1593 cm−1 appeared instead. The C–N stretch at 1258 cm−1 in MIL-101Cr-NH2 was clearly seen in the spectrum. All of the frequency bands of the MIL-101Cr framework remained unchanged after the incorporation and reduction of palladium.

The structure of MIL-101Cr was verified by powder X-ray diffraction. The PXRD patterns of our synthesized samples of MIL-101Cr matched with those reported in the literature.[9] The structure of MIL-101Cr remained intact after all of the postsynthetic modification steps and no obvious changes in the PXRD patterns were detected, which confirmed that MIL-101Cr-NH2 was chemically and hydrothermally stable and that it could withstand harsh conditions, such as those applied during the nitration step (with mixtures of concentrated acids) and during the reduction of PdII by using excess NaBH4 (Figure 1). The systematic changes in the relative peak intensities in the PXRD patterns with different Pd loadings, for example, those within the regions 2θ=5–6° and 8–9°, indicated that Pd had been incorporated into the pores of MIL-101Cr.

Figure 1

PXRD patterns of MIL-101Cr, MIL-101Cr-NO2, MIL-101Cr-NH2, and 4-, 8-, 12-, and 16 wt %Pd@MIL-101Cr-NH2. The palladium in all four Pd-loaded samples is Pd0.

PXRD patterns of MIL-101Cr, n class="Chemical">MIL-101Cr-NO2, MIL-101Cr-NH2, and 4-, 8-, 12-, and 16 wt %Pd@MIL-101Cr-NH2. The palladium in all four Pd-loaded samples is Pd0.

The weight losses of the synthesized samples upon heating in air were studied by using TGA (Supporting Information, Figure S4). MIL-101Cr was stable up to 300 °C. The thermal stability of MIL-101Cr seemed to be slightly weakened after amino-functionalization and even more so after immobilization of the Pd nanoparticles. The amount of water and solvent molecules in the pores differed between the samples, as shown by the different weight losses during heating.

N2-sorption measurements were performed for as-synthesized MIL-101Cr, amino-functionalized MIL-101Cr-NH2, and 4-, 8-, 12-, and 16 wt %Pd@MIL-101Cr-NH2 samples. The N2-adsorption isotherms, the BET surface areas, and the pore volumes are shown in Figure 2 and Table 1. The BET surface areas of MIL-101Cr and MIL-101Cr-NH2 are 2966 and 2869 m2g−1, respectively, which are in agreement with the reported values.[24] Compared with the parent MIL-101Cr-NH2, the BET surface area and pore volume decreased with higher Pd loading (Table 1). The adsorption isotherms of all samples show stepwise N2 uptake, at p/p0 values up to 0.25, thus indicating the existence of different types of cage.

Figure 2

Nitrogen-adsorption isotherms of MIL-101Cr, MIL-101Cr-NH2, and 4-, 8-, 12-, and 16 wt %Pd@MIL-101Cr-NH2. The measurements were performed at 77 K.

Table 1

Specific surface areas and pore volumes of MIL-101Cr and the functionalized MIL-101Cr materials

Material	BET surface area [m2 g−1]	Total pore volume [m3 g−1]
MIL-101Cr	2966	1.403
MIL-101Cr-NH2	2869	1.389
4 wt %Pd@MIL-101Cr-NH2	1831	0.870
8 wt %Pd@MIL-101Cr-NH2	1623	0.767
12 wt %Pd@MIL-101Cr-NH2	1396	0.678
16 wt %Pd@MIL-101Cr-NH2	1321	0.625

Nitrogen-adsorption isotherms of n class="Chemical">MIL-101Cr, MIL-101Cr-NH2, and 4-, 8-, 12-, and 16 wt %Pd@MIL-101Cr-NH2. The measurements were performed at 77 K.

Specific surface areas and pore volumes of MIL-101Cr and the functionalized n class="Chemical">MIL-101Cr materials

SEM analysis showed that the MIL-101Cr crystals were formed as octahedra, with sizes between 0.3 and 1.0 μm (Supporting Information, Figure S5). TEM analysis indicated that the octahedral morphology of MIL-101Cr was preserved throughout all post-modification steps, which include nitration, reduction of the nitro groups, introduction of PdII, and reduction of PdII into Pd0. The sizes of the Pd particles ranged from 20 to 30 Å, which were comparable to the pore diameters of unfunctionalized MIL-101Cr (29 and 34 Å, Figure 3).[9] However, on increasing the Pd loading in MIL-101Cr-NH2 from 11.3 to 15.1 wt %, the size of the Pd nanoparticles increased up to about 8 nm (Figure 3). These larger Pd nanoparticles were found on the surfaces of the 16 wt %Pd@MIL-101Cr-NH2 crystals, as shown by TEM. Electron tomography was performed on crystals of 8 wt %Pd@MIL-101Cr-NH2 to study the distribution of the Pd nanoparticles. The reconstructed 3D tomographic images show that the nanoparticles are embedded inside the crystals and are homogeneously distributed throughout the crystals (Figure 4 and the Supporting Information, Movies S1 and S2). Electron diffraction and high-resolution transmission electron microscopy (HRTEM) images show that the Pd nanoparticles are crystalline (Supporting Information, Figure S6).[28]

Figure 3

TEM images of samples of 4-, 8-, 12-, and 16 wt %Pd@MIL-101Cr-NH2 with different Pd loadings (4.4, 8.1, 11.3, and 15.1 wt % Pd, respectively). Insets show the size distributions of the Pd particles as histograms, each of which was formed from about 100 Pd particles in several crystals in the TEM images. The x axis in the histograms is the size in nm and the y axis is the number of particles.

Figure 4

2D section of the 3D tomographic image of two Pd@MIL-101Cr-NH2 crystals with 8 wt % Pd loading, as reconstructed by the electron tomography from a tilt series of 83 TEM images with a tilt range of 41° and a tilt step of 0.5°. For the tilt series and the full tomogram, see the Supporting Information, Movies S1 and S2, respectively.

TEM images of samples of 4-, 8-, 12-, and 16 wt %Pd@n class="Chemical">MIL-101Cr-NH2 with different Pd loadings (4.4, 8.1, 11.3, and 15.1 wt % Pd, respectively). Insets show the size distributions of the Pd particles as histograms, each of which was formed from about 100 Pd particles in several crystals in the TEM images. The x axis in the histograms is the size in nm and the y axis is the number of particles.

2D section of the 3D tomographic image of two Pd@n class="Chemical">MIL-101Cr-NH2 crystals with 8 wt % Pd loading, as reconstructed by the electron tomography from a tilt series of 83 TEM images with a tilt range of 41° and a tilt step of 0.5°. For the tilt series and the full tomogram, see the Supporting Information, Movies S1 and S2, respectively.

A better technique to study the size and distribution of the Pd nanoparticles is X-ray total scattering, because TEM only gives some local information and does not necessarily represent the nature of the bulk samples. Owing to their small particle sizes and the low concentrations of the Pd nanoparticles, no Bragg peaks for the Pd nanoparticles could be clearly seen in the PXRD patterns that were recorded on an in-house powder X-ray diffractometer. Therefore, we recorded the total scattering intensities of the 4-, 8-, 12-, and 16 wt %Pd@MIL-101Cr-NH2 samples at the beamline P02.1, Petra III, DESY, by using a high-energy synchrotron source (60 keV). The pair distribution function (PDF), G(r), was obtained by performing a Fourier transform of the total scattering intensities, including both Bragg peaks and diffuse scattering.[29] This function describes the statistical distribution of all of the interatomic distances in a sample and, thus, is particularly suitable for investigating nanoparticles,[30] which usually give very broadened or even no Bragg peaks by using conventional PXRD. The PDFs of the 4-, 8-, 12-, and 16 wt %Pd@MIL-101Cr-NH2 samples were extracted from the total X-ray scattering signal truncated at a maximum Q value (Qmax) of 19 Å−1 (Figure 5; Supporting Information, Figure S7). For comparison, the PDFs that were calculated from the structural model of MIL-101Cr and from the structure of an ideal Pd nanoparticle (space group Fmm) with a diameter of 25 Å are also shown in Figure 5 (for the detailed procedure see the Supporting Information). The peaks in the observed PDFs of the 4-, 8-, 12-, and 16 wt %Pd@MIL-101Cr-NH2 samples correspond to the interatomic distances in both MIL-101Cr and the Pd nanoparticle. For example, the peak at 2.0 Å corresponds to the nearest mean Cr–O distance, which is not present in the calculated PDF of the pure Pd nanoparticle (Figure 5, top). The peak at 2.9 Å corresponds to the shortest distance between two Pd atoms in the Pd nanoparticle. The “ripples” on either side of the first peak in the calculated PDF of the Pd nanoparticle originate from termination effects that are due to the finite Qmax of 19 Å−1 that was used for the calculations and can be neglected. Whereas the peaks in the calculated PDF of the Pd nanoparticle are relatively sharp, those above 2.5 Å in the calculated PDF of MIL-101Cr are rather broad, owing to multiple overlapping of different atom–atom correlations of similar distances. The peak at 1.3 Å results from overlapping of the nearest C–C and C–O distances in MIL-101Cr. The integrated intensity of a well-defined PDF peak yields the coordination number of that specific atomic pair.[28] For example, the integrated intensity of the first peak at 2.9 Å in the calculated PDF of the pure Pd nanoparticle will yield 12, which corresponds to the coordination number of Pd. In case of multielement samples, such as MIL-101Cr and 4-, 8-, 12-, and 16 wt %Pd@MIL-101Cr-NH2, the peak intensities in the PDF are multiplied by a weighting factor to take into account the different scattering factors of the chemical species. Thus, in the observed PDFs of 4-, 8-, 12-, and 16 wt %Pd@MIL-101Cr-NH2, the peak intensities can be related to the multiplicity and abundance of the atomic pairs, as well as to atomic scattering factors. By comparing the PDFs of the 4-, 8-, 12-, and 16 wt %Pd@MIL-101Cr-NH2 samples (Figure 5), a change in the relative peak intensities can be observed; in particular, with higher Pd loading, the intensity of the peak at 2.9 Å, which corresponds to the Pd–Pd distance, increases in comparison with that of the peak at 2.0 Å, which is assigned to the Cr–O distance. This result is reasonable because, with higher Pd loading, the number of Pd atoms and, hence, the number of Pd–Pd distances of 2.9 Å is expected to increase, whereas the number of Cr–O distances of 2.0 Å in the MIL-101Cr-NH2 should remain constant. When the PDFs of the 12- and 16 wt %Pd@MIL-101Cr-NH2 samples are compared to the PDF of the 25 Å Pd nanoparticles (calculated), one can see that the peaks at r>4.5 Å mainly originate from the Pd nanoparticles. In fact, the PDF of 16 wt %Pd@MIL-101Cr-NH2 can be almost fully described by the structure of the Pd nanoparticles with an average size of 25 Å. The PDF of 12 wt %Pd@MIL-101Cr-NH2 is similar to that of 16 wt %Pd@MIL-101Cr-NH2, except that the former attenuated at a shorter distance (r). This result indicates that the size of the Pd nanoparticles in 16 wt %Pd@MIL-101Cr-NH2 is slightly larger than that in 12 wt %Pd@MIL-101Cr-NH2. On the other hand, the PDFs of the 4- and 8 wt %Pd@MIL-101Cr-NH2 samples are very similar at r>4.5 Å and attenuate at much lower r values compared to the PDFs of 12- and 16 wt %Pd@MIL-101Cr-NH2. This result indicates that the sizes of the Pd particles are similar in 4- and 8 wt %Pd@MIL-101Cr-NH2 and are smaller than those in 12- and 16 wt %Pd@MIL-101Cr-NH2. It is worth mentioning that the PDF can also tell if there is agglomeration of the nanoparticles in the sample. Because no peaks are observed at large r values (>20 Å), it can be concluded that there is no significant agglomeration of Pd nanoparticles in the four Pd@MIL-101Cr-NH2 samples.

Figure 5

PDFs of samples of 4-, 8-, 12-, and 16 wt %Pd@MIL-101Cr-NH2, compared with the calculated PDFs of MIL-101Cr (bottom) and of a Pd nanoparticle (diameter: 25 Å, top). For the sake of clarity, the PDFs are vertically stacked by 10 Å−2.

PDFs of samples of 4-, 8-, 12-, and 16 wt %n class="Chemical">Pd@MIL-101Cr-NH2, compared with the calculated PDFs of MIL-101Cr (bottom) and of a Pd nanoparticle (diameter: 25 Å, top). For the sake of clarity, the PDFs are vertically stacked by 10 Å−2.

Suzuki–Miyaura cross-coupling reaction with Pd@MIL-101Cr-NH: In an initial series of reactions, we screened the reactivity of p-bromoanisol (1 a) with different organoboron reagents and bases (Table 2). A similar screening of the reactivity of p-bromotoluene is shown in the Supporting Information (Table S3). Unless otherwise specified, water was used as the solvent and the reactions were performed at room temperature. Based on the characterization of the different catalysts, screening of the reaction conditions was performed with 8 wt %Pd@MIL-101Cr-NH2 (3 mol % Pd loading). This material had small Pd nanoparticles, which ranged from 20 to 30 Å. In contrast to 12- and 16 wt %Pd@MIL-101Cr-NH2, the 8 wt % catalyst did not show agglomeration of palladium on its surface. Moreover, 8 wt %Pd@MIL-101Cr-NH2 contained a more favorable wt % Pd per gram of material than 4 wt %Pd@MIL-101Cr-NH2, which is desirable from an atom-economy viewpoint. Nevertheless, the catalytic activity of all of the catalysts has been compared (see below).

Table 2

Suzuki–Miyaura coupling reactions catalyzed by Pd@MIL-101Cr-NH2; optimization of the reaction conditions[a]

Entry	ArBr	Boron (equiv)	Base	Solvent	t [h]	Yield [%][b]
1	1 a	PhB(OH)2 (2 a, 1.5)	Cs2CO3	water	6	>99
2	1 a	2 a (1.5)	NaOH	water	6	>99
3	1 a	2 a (1.5)	Na2CO3	water	6	71
4	1 a	2 a (1.5)	K2CO3	water	6	82
5	1 a	2 a (1.5)	BaCO3	water	6	94
6	1 a	2 a (1.5)	CsF	water	6	71
7	1 a	2 a (1.5)	NaOAc	water	6	–[c]
8	1 a	2 a (1.5)	iPr2NH	water	6	89
9	1 a	2 a (1.5)	iPr2EtN	water	6	86
10	1 a	2 a (1.5)	Et3N	water	6	83
11	1 a	2 a (1.5)	C2H4(NH2)2	water	6	–
12	1 a	2 a (1.5)	No base	water	6	3
13	1 b	2 a (1.5)	Cs2CO3	water	6	>99
14	1 c	2 a (1.5)	Cs2CO3	water	6	84
15	1 c	2 a (1.5)	Cs2CO3	water/EtOH	0.5	>99
16	1 a	2 a (1.2)	Cs2CO3	water	4	68
17	1 a	PhBpin (2 b, 1.2)	Cs2CO3	water	4	99
18	1 a	PhBF3K (2 c, 1.2)	Cs2CO3	water	4	83
19	1 a	2 b (1.2)	K2CO3	water	6	>99
20	1 a	2 b (1.2)	K2CO3	water	3	58
21[d]	1 a	2 b (1.2)	K2CO3	water	3	56
22[e]	1 a	2 b (1.2)	K2CO3	water	3	33
23[f]	1 a	2 b (1.2)	K2CO3	water	3	36

[a] Unless otherwise stated, the organoboron reagent (0.12 or 0.15 mmol), aryl halide (0.1 mmol), base (0.2 mmol), and 8 wt %Pd@MIL-101Cr-NH2 (4 mg, 0.003 mmol, 3 mol % Pd loading) were suspended in deionized water (2 mL). The mixture was stirred vigorously at RT for 6 h in a glass vial in air (see the Supporting Information); [b] yields determined by 1H NMR spectroscopy with 1,2,4,5-tetrachloro-3-nitrobenzene as an internal standard; [c]>95 % conversion into anisole (i.e., dehalogenation); [d] 4 wt %Pd@MIL-101Cr-NH2 (3 mol % Pd loading); [e] 12 wt %Pd@MIL-101Cr-NH2 (3 mol % Pd loading); [f] 16 wt %Pd@MIL-101Cr-NH2 (3 mol % Pd loading).

Suzuki–Miyaura coupling reactions catalyzed by Pd@n class="Chemical">MIL-101Cr-NH2; optimization of the reaction conditions[a]

Influence of the base: The nature of the base in a Suzuki–Miyaura coupling reaction greatly affects the efficiency of the reaction.[31] In general, the reactivities of different bases in the coupling reactions of aryl iodides and bromides decreased in the order: OH−>CO32−≫AcO−. At high concentrations, bases can lead to the formation of unreactive boronates (ArB(OH)3−) and, thus, the [OH−]/[ArB(OH)2] ratio needs to be controlled. Moreover, the countercations can have an unexpected deceleration effect on the transmetalation step, by coordinating to the hydroxo ligand of [ArPd(OH)L2] intermediates.[31] With these arguments in mind, we screened a wide selection of bases and observed the stability of the 8 wt %Pd@MIL-101Cr-NH2 catalyst, as shown in Table 2. When PhB(OH)2 (2 a) and Cs2CO3 were used, 4-methoxybiphenyl (3 a) was obtained in quantitative yield within only 6 h (Table 2, entry 1). We were particularly interested in finding a less-expensive base than Cs2CO3 that could afford comparably high yields without harming the material. Thus, a range of other bases was tested (Table 2, entries 2–12). Sodium hydroxide also afforded excellent results (Table 2, entry 2); however, MIL-101Cr-NH2 was not stable under these strongly basic conditions (Supporting Information, Figure S8).[24] For this reason, carbonates were preferred because the release of hydroxy anions into the aqueous medium was slower. This property prevented the decomposition of the framework during the catalysis and also allowed for better control over the [OH−]/[ArB(OH)2] ratio. Other carbonates, such as Na2CO3, K2CO3, and BaCO3 afforded compound 3 a in 71 %, 82 %, and 94 % yield, respectively (Table 2, entries 3–5). CsF gave a comparable yield to those that were obtained with carbonates (Table 2, entry 6). Acetates have been shown to be a poor choice of base, because they lead to a slow transmetalation step.[31] Although various amounts of the coupling product (3 a) were detected in several experiments with NaOAc as the base, most of the starting material underwent dehalogenation to produce anisole, which occurred faster than the coupling reaction (Table 2, entry 7). We also tested a number of organic bases (Table 2, entries 8–11) because their liquid nature may facilitate the large-scale application of this reaction. Potential problems owing to pore blocking in the case of insufficient solubility of the base would also be avoided. We found that several bases, such as iPr2NH, iPr2EtN, and Et3N, afforded the cross-coupled product (3 a) in very high yields (Table 2, entries 8–10). However, the stability of the nanoparticles with these bases was compromised. A comparative leaching study was performed by using one inorganic base (Cs2CO3) and one organic base (iPr2NH) in the coupling reaction of compound 1 a with compound 2 b in water at room temperature. After completion of the reaction (6 h), the mixture was filtered through a short pad of Celite and the Pd content in the filtrate was analyzed by ICP-OES analysis. For the reaction with Cs2CO3, a remarkable stability of the catalyst was observed and only 0.17 ppm Pd was detected in the filtrate. Higher Pd content was detected in the filtrate when iPr2NH was used (5.64 ppm), which, although still relatively low, was not negligible. This higher leaching was a consequence of the high affinity of palladium for nitrogenated ligands. Furthermore, if a chelating base (ethylenediamine) was used (Table 2, entry 11), all of the palladium was extracted from the support and the material lost its catalytic activity. Experiments that were performed in the absence of base showed conversions from 3 % to 21 % (Table 2, entry 12). This result may be attributed to the presence of free amino groups in the framework, which can act as basic functionalities.

Solvent and substrate polarity: Under the conditions shown in Table 2, entry 1, we compared the reactivity of compound 1 a with those of p-bromotoluene (1 b) and electron-poor bromide 1 c (Table 2, entries 13 and 14). Surprisingly, in repeated experiments, the neutral (1 b) and electron-rich substrates (1 a) reached full conversion within 6 h, whereas the electron-poor compound 1 c was converted into the corresponding biphenyl in only 84 % yield. This unexpected reactivity, which was opposite to that in organic solvents, owing to the electronic nature of these substrates,[32] may be attributed to the different solubilities of compounds 1 a–1 c in water. Indeed, when water was replaced by a mixture of EtOH and water (1:1, in which differences in solubility were minimized), bromide 1 c afforded the corresponding coupling product in excellent yield within only 30 min (Table 2, entry 15).

Organoboron reagents: Next, we evaluated other organoboron derivatives. To clearly observe any differences in their reactivities, the reactions were stopped after only 4 h and only 1.2 equivalents of the organoboron reagents were used (Table 2, entries 16–18). Pinacol phenylboronate (2 b) afforded the highest yield, followed by phenyltrifluoroborate potassium salt (2 c) and phenylboronic acid (2 a). The lower reactivity of compound 2 a can be explained by its preference to form boroxine trimers, which are inactive under these mild reaction conditions.[33] Thus, we decided to continue our investigation by using pinacolate ester 2 b. Furthermore, the use of compound 2 b compensated for the loss of reactivity that was observed when Cs2CO3 was replaced by K2CO3 and quantitative yields could be obtained with the latter base by using only 1.2 equivalents of PhBpin (Table 2, entry 19).

Comparison of Pd@MIL-101Cr-NH: With the optimized conditions in hand, we compared the activities of the 4-, 8-, 12-, and 16 wt %Pd@MIL-101Cr-NH2 catalysts (Table 2, entries 20–23). The reactions were analyzed before full conversion had occurred (3 h) by using a Pd loading of 3 mol %. Whereas the reactions with 4- and 8 wt % Pd@MIL-101Cr-NH2 afforded very similar and satisfactory yields (Table 2, entries 20 and 21), the reactions with 12- and 16 wt %Pd@MIL-101Cr-NH2 afforded significantly lower yields (Table 2, entries 22 and 23). This behavior was in agreement with the TEM and PDF analyses, which showed that the nanoparticles in the 12- and 16 wt %Pd@MIL-101Cr-NH2 catalysts were larger with broader size distributions. The poorer performance of the MOFs with higher Pd loadings may be explained by an increased difficulty in the diffusion of the substrates and the products through these materials, owing to pore blocking.

Leaching tests: The stability and performance of the immobilized nanoparticles were evaluated according to the approach reported by Sheldon et al.[34] A hot-filtration test was performed by stirring the 8 wt %Pd@MIL-101Cr-NH2 catalyst in boiling water for 15 min, followed by filtration through a pad of Celite. The filtrate was subsequently used as the solvent for the reaction of compound 1 a with compound 2 a under similar reaction conditions to those in Table 2, entry 1, but in the absence of the palladium catalyst. After 6 h, only 3 % conversion was detected by 1H NMR spectroscopy. The Pd leached into solution after catalysis when carbonates were used as the base was only 0.17 ppm (see above). Together with the high recyclability of the catalyst (see below), these results clearly demonstrate the stability of the Pd nanoparticles in MIL-101Cr-NH2.

Substrate scope: Based on this series of observations (Table 2), we screened a variety of substrates with 8 wt %Pd@MIL-101Cr-NH2 as the catalyst (3 mol % Pd loading) and K2CO3 as the base, because K2CO3 is environmentally benign and about 20-times cheaper than Cs2CO3. Moreover, potassium has been shown to be the countercation that infers the least-intense decelerating effect in the transmetalation step.[31] These results are shown in Table 3. The reactions were either performed in water or in water/EtOH mixtures (1:1 v/v). Aryl bromides 1 a–1 c were easily coupled within short reaction times (6 h in water, 30 min in water/EtOH) at room temperature (Table 3, entries 1–3). 2-Bromonaphthalene (1 d), which showed very low solubility in water at ambient temperatures, afforded compound 3 d in only 22 % yield. However, 97 % yield was obtained after only 2 h in water/EtOH, in which the solubility of the substrate was increased (Table 3, entry 4). Interestingly, the reaction of 9-bromoanthracene (1 e) proceeded very slowly, even in water/EtOH after a prolonged reaction time (Table 3, entry 5). This result may be due to the large size of this substrate, which results in limited access to 8 wt %Pd@MIL-101Cr-NH2, and the catalysis mainly occurred at the Pd nanoparticles on the external surfaces of the MOF catalyst particles. Importantly, heteroaromatic substrates 1 f and 1 g could also be coupled in high yields (Table 3, entries 6–7). However, a lower yield was obtained for 2-bromopyridine (1 h; Table 3, entry 8). Slightly higher leaching of palladium was observed for substrates that contained heteroatoms: 3.10 ppm Pd was found in the filtrate after the reaction of bromopyridine (1 g), whereas only 0.17 ppm was found after the coupling of compound 1 a (see above). The coupling reactions of aryl chlorides required more elevated temperatures, which, disappointingly, were more favorable for the formation of byproducts (that is, dehalogenation (ArH) or homocoupling (Ar–Ar) products). Activated aryl chlorides 1 i and 1 j successfully reacted in water and they returned good-to-excellent yields with little or no formation of undesired byproducts (Table 3, entries 9 and 10). However, if EtOH was used as a co-solvent, both activated and inactivated substrates 1 j–1 m returned very good conversions but only fair yields, with the concomitant formation of significant amounts of homocoupling and dehalogenation products (Table 3, entries 10–13).

Table 3

Suzuki coupling reactions of various aryl bromides and chlorides catalyzed by 8wt %Pd@MIL-101Cr-NH2[a]

Entry	Substrate	Product	Yield (water, 6 h) [%][b]	Yield (water/EtOH, 0.5 h) [%][b]
1			>99 (94)	95
2			>99 (96)	>99
3			84	99 (96)
4			22	97 (94)[c]
5			0	19[d]
6			99 (89)	–
7			99 (92)	–
8			47	95 (90)
9			99 (95)[e]	–
10			70 (62)[f,e]	54[g,e]
11			52 (47)[h,e]	50[i,e]
12			4[j,e]	44[k,e]
13			3[l,e]	36[m,e]

[a] Unless otherwise stated, the pinacol phenylboronate (0.12 mmol), aryl halide (0.1 mmol), K2CO3 (0.2 mmol), and 8 wt %Pd@MIL-101Cr-NH2 (4 mg, 0.003 mmol, 3 mol % Pd loading) were suspended in either deionized water (2 mL) or water/EtOH (1:1). Then, the mixture was stirred vigorously at RT in air for 6 h (water) or 30 min (water/EtOH) in a glass vial (see the Supporting Information); [b] yield determined by 1H NMR spectroscopy with 1,2,4,5-tetrachloro-3-nitrobenzene as an internal standard; yield of the isolated product is given in parentheses; [c] 2 h; [d] 16 h; [e] 80 °C; [f] 96 % conversion, 13 % homocoupling; [g] 89 % conversion, 9 % homocoupling, 16 % dehalogenation; [h] 68 % conversion, 13 % homocoupling; [i] 89 % conversion, 12 % homocoupling, 15 % dehalogenation; [j] 98 % conversion; toluene was the main byproduct; [k] 100 % conversion, 12 % homocoupling, 32 % dehalogenation; [l] 44 % conversion; anisole was the main byproduct; [m] 79 % conversion, 10 % homocoupling, 23 % dehalogenation.

Suzuki coupling reactions of various aryl bromides and n class="Chemical">chlorides catalyzed by 8wt %Pd@MIL-101Cr-NH2[a]

Recycling tests: The recyclability of 8 wt %Pd@MIL-101Cr-NH2 was investigated in the reaction of p-bromotoluene. The experiments were performed on a larger scale (1 mmol) than in previous investigations to facilitate separation and minimize handling errors. After each run, the solid catalyst was separated by centrifugation, washed with water and EtOH (twice each), dried under vacuum, and reused. As shown in Figure 6, the catalyst remained highly active, even after 10 cycles.[35] Interestingly, the PXRD pattern after the tenth cycle was different from that of MIL-101Cr-NH2, but similar to that of another phase, presumably MIL-88B-NH2 (Supporting Information, Figure S9).[36] This phase-transformation phenomenon is currently under investigation in our laboratories. TEM analysis showed the presence of some limited aggregation of Pd nanoparticles on the surfaces of the MOF catalysts after the 10th cycle (Supporting Information, Figure S10). Because the size of the nanoparticles remained unchanged, the aggregation did not seem to significantly affect the efficiency of the catalyst. Elemental analysis of the recycled 8 wt %Pd@MIL-101Cr-NH2 material showed an overall loss of 0.96 wt % Pd after 10 runs. Yields of between 97 % and >99 % were obtained in the recycling runs. Based on these excellent results, we decreased the Pd loading by 20 times to 0.15 mol % and performed the reactions on the same scale and under similar conditions to those in Table 2, entry 19, and the reaction proceeded to completion within 4 h. This result indicates that Pd@MIL-101Cr-NH2 is a highly efficient catalyst.

Figure 6

Recycling experiments on a 1 mmol scale; reaction time of each run: 0.5 h.

Conclusion

In summary, we have provided evidence for a strong relationship between the amount of palladium that is loaded into MIL-101Cr-NH2 and the catalytic properties of the corresponding system. The MIL-101Cr-NH2 material was found to be an excellent support for palladium nanoparticles. The size, distribution, and reactivity of palladium nanoparticles were shown to be influenced by the amount of palladium metal that was impregnated into the framework. The nanoparticles were investigated by using both conventional and advanced methods, such as electron tomography and PDF analysis. The best performance was observed for MIL-101Cr-NH2 that contained about 8 wt % palladium nanoparticles. This result indicates that an optimal ratio was achieved between the amount of loaded catalyst and the number of pores that remained accessible for the diffusion of the substrates and the products. The catalyst 8 wt %Pd@MIL-101Cr-NH2 displayed remarkable activity in the Suzuki–Miyaura cross-coupling reaction. Aryl and heteroaryl bromides and activated aryl chlorides were successfully coupled within short reaction times under mild conditions (ambient temperature and 80 °C, respectively) in benign solvents (water and EtOH). Our investigations have also shown that the solubility and hydrophobicity of the substrates significantly affect their ability to access the catalyst pores in water, but this limitation is minimized by employing EtOH as a co-solvent. The excellent recyclability of the 8 wt %Pd@MIL-101Cr-NH2 catalyst makes it a promising candidate for large-scale Suzuki–Miyaura cross-coupling reactions. This thorough study of the structural and catalytic properties of metallic nanoparticles that are supported on a MOF material can stand as a model for detailed investigations of other related heterogeneous systems, in search for a better understanding of the underlying mechanism of heterogeneous catalysis with MOFs.

Experimental Section

General: All chemicals and solvents were used as received without further purification. Powder X-ray diffraction (PXRD) was performed in the Bragg–Brentano geometry on a PANalytical X’Pert PRO diffractometer that was equipped with a Pixel detector and CuKα1 radiation (λ=1.5406 Å). The samples were dispersed over zero-background Si plates. SEM and energy-dispersive spectroscopy (EDS) were performed on a JEOL-7000F field-emission scanning electron microscope operating at 12 kV. The samples were precoated with a thin carbon film. The presence of possible byproducts was investigated by using EDS. TEM observations were performed on a JEOL TEM (JEM-2100F) operating at 200 kV. The MIL-101Cr samples were dispersed in absolute EtOH and then treated by ultrasonication for 2 min. A droplet of the suspension was transferred onto a copper grid and dried in air. A tilt series of TEM images for electron tomography was recorded with a tilt step of 0.5° and an angular range of 41° (83 images in total). TomoJ[37, 38] was employed for the tomographic reconstruction. The images were aligned by using the Pd particles themselves as markers. Ten iterations of the ART[39] algorithm were performed in TomoJ with a relaxation coefficient of 0.1. ImageJ[40] was used to preprocess the tilt series of TEM images. Thermogravimetric analysis (TGA) was performed under a flow of air between 20 and 500 °C with a heating rate of 2 °C min−1 on a thermogravimetric analyzer (PERKIN ELMER TGA 7) that was equipped with a platinum pan. N2-adsorption isotherms were recorded at 77 K on a Micromeritics ASAP2020 analyzer. The samples for N2 adsorption were first exchanged with absolute EtOH before activation. The samples were degassed by heating at 120 °C overnight prior to the sorption measurements. Specific surface areas were calculated from the data in the adsorption branch at p/p0=0.05–0.30. The total pore volume was calculated from the uptake at p/p0=0.99. FTIR spectroscopy was performed on a Varian 670-IR spectrometer. Elemental analysis was performed on a Carlo Erba Flash 1112 elemental analyzer and the Cr and Pd contents were determined by inductively coupled plasma-optical emission spectrometry (ICP-OES) on a Varian Vista MPX ICP-OES at Medac Ltd, Chobham, UK. High-energy X-ray diffraction experiments were carried out at a fixed energy of approximately 60 keV at a wavelength of 0.2073 Å (determined by using a CeO2 standard) at the beamline P02.1 of Petra III, Desy, Hamburg. The 2D diffraction images were integrated into a linear scattering signal by using the Fit2D software.[41] The pair distribution function (PDF), G(r), was obtained by sine Fourier transform of the normalized scattering intensity to a maximum Q value, Qmax4=19 Å−1, by using the PDFgetX2 software.[42] The PDF of the Pd nanoparticles was calculated with the PDFgui software[43] (for a detailed procedure, see the Supporting Information). 1H NMR spectra were recorded at 400 MHz; 13C NMR spectra were recorded at 100 MHz on a Bruker Advance spectrometer. 1H and 13C NMR chemical shifts (δ) are reported in ppm relative to tetramethylsilane, with the solvent resonance as an internal standard (CDCl3: δH=7.26 ppm, δC=77.16 ppm; [D6]DMSO: δH=2.50 ppm, δC=39.5 ppm). Coupling constants (J) are given in Hz. High-resolution mass spectra (HRMS) were recorded on a Bruker microTOF ESI-TOF mass spectrometer.

Synthesis of Pd@n class="Chemical">MIL-101Cr-NH: Samples of MIL-101Cr (Cr3F(H2O)2O[(CO2)-C6H4-(CO2)]3⋅nH2O),[9, 44] MIL-101Cr-NO2 (Cr3(H2O)2O[(CO2)-C6H3NO2-(CO2)]3⋅nH2O),[24] and MIL-101Cr-NH2 (Cr3(H2O)2O[(CO2)-C6H3NH2-(CO2)]3⋅nH2O)[24] were prepared according to literature procedures. Samples of Pd@MIL-101Cr-NH2 with four different loadings of Pd nanoparticles that were intended to be 4, 8, 12, and 16 wt % Pd were prepared by using analogous procedures, as given below.

Synthesis of MIL-101Cr-NH: For MIL-101Cr-NH2-PdCl2 affording 8 wt %Pd@MIL-101Cr-NH2, the following procedure was used: A mixture of MIL-101Cr-NH2 (0.240 g, 3.15×10−4 mol), [PdCl2(MeCN)2] (0.055 g, 2.11×10−4 mmol), and anhydrous CH2Cl2 (30 mL) was stirred overnight at RT (18 h). The solid was obtained by centrifugation, washed with anhydrous CH2Cl2 (2×10 mL), and dried at 40 °C for 2 h under vacuum. To obtain samples with 4, 12, and 16 wt % Pd loading, the amount of [PdCl2(MeCN)2] that was used to afford 8 wt % Pd loading was multiplied by 0.5, 2, and 3, respectively; the amount of solvent was kept the same.

Synthesis of 8 wt %Pd@n class="Chemical">MIL-101Cr-NH: NaBH4 (0.076 g, 2 mmol) was added portion-wise to a mixture of MIL-101Cr-NH2-PdCl2 (0.140 g) and EtOH (30 mL) at 0 °C. The green mixture became darker after about 2 min. The reaction was stirred for an additional 2 h at 0 °C to complete the reduction. The solid was obtained by centrifugation, washed twice with water (45 mL) and once with EtOH (45 mL), and the dark-green solid was dried at 85 °C overnight (18 h) in air.

Digestion of MIL-101Cr: In a typical experiment, the material (10 mg) was suspended in an aqueous solution of NaOH (2.0 m, 1 mL) and ultrasonicated for 30 min. The water was evaporated under vacuum and the residue was analyzed by 1H NMR spectroscopy in deuterated DMSO (Supporting Information, Figure S1).

General procedure for the Suzuki–Miyaura cross-coupling reactions: Unless otherwise stated, pinacol phenylboronate (24.5 mg, 0.12 mmol), the aryl bromide (0.1 mmol), K2CO3 (27.6 mg, 0.2 mmol), and 8 wt %Pd@MIL-101Cr-NH2 (4 mg, 0.003 mmol) were suspended in deionized water (2 mL) and the mixture was stirred vigorously at RT for 6 h. Then, the mixture was extracted into CH2Cl2 (4 mL), the organic phase was separated, the volatile compounds were evaporated, and the residue was analyzed by 1H NMR spectroscopy. To isolate the final product, the reaction was repeated on a larger scale (×4) and the crude reaction mixture was purified by column chromatography on silica gel (pentane or pentane/EtOAc mixtures).