1. Introduction
Metal–organic frameworks, or MOFs, have been widely studied over the past 20 years for possible applications in gas storage, separation, and sensors [
1]. The modular nature of the synthesis has produced more than 50,000 structures, since a wide range of inorganic cations (“nodes”) and organic ligands can be combined together [
2]. The appropriate choice of the chemical building blocks affords a porous three-dimensional solid that exhibits extremely high surface area. However, the pores are often filled by the solvent after the synthesis and must be removed for the applications mentioned above. In addition, the mechanical stability is often limited after the solvent is removed from the porous solid. In 2008 a significant advance in MOF synthesis was reported by Lillerud and co-workers at the University of Oslo [
3]. Cavka et al. reported the first synthesis of UiO-66-H using 1,4-benzenedicarboxylic acid as the organic linker and Zr clusters as the inorganic “nodes” [
3]. The material is stable over a large pH range in solution and showed excellent thermal stability [
4,
5,
6]. The ideal structure of a UiO-66 unit-cell consists of a zirconium-oxo cluster, Zr
6O
4(OH)
4, to form a node and six terephthalic acid linkers. In the ideal UiO-66 framework, each Zr ion takes the vertex in an octahedron with the μ
3-O (triply bridging oxygen) and μ
3-OH groups forming the eight faces [
7,
8]. The external coordination spheres of the zirconium ions are formed by carboxylate donors forming the twelve edges of the octahedron. In this way, the overall metal–organic framework forms as twelve clusters connected by terephthalic acid linkers to each other into a face-centered cubic array (fcu). As the nodes occupy corner and face sites of the fcu-type unit cell, the linkers will divide the cell into two types of pores, which are common in fcu-type frameworks. A single large central octahedral pore (~11 Å) in the center of the unit cell is formed by the space between the six facial nodes. Smaller tetrahedral pores (~8 Å) are then formed around the central pore by the linkers between three facial and one corner node, one for each corner of the cube. There are smaller triangular windows (~6 Å) that connect the two pores. The effective size of the triangular window is temperature dependent because the linker can rotate [
9]. Otherwise, molecules with kinetic diameters larger than the window would not be able to pass through the framework.
The fundamental adsorption, desorption and diffusion processes of hydrocarbons into MOFs has been studied because of their industrial importance [
10,
11,
12,
13,
14,
15,
16,
17,
18,
19,
20]. In particular, the separation of the BTEX family of compounds is very challenging due to the similar physical properties. Industrially the separation is accomplished with a simulated moving bed using zeolites. An interesting observation was discovered by Barcia et al. when they performed breakthrough measurements with ethylbenzene and xylene isomers on UiO-66-H pellets [
21]. At 398 K and 423 K,
o-xylene was retained longer than
p-xylene from a mixture of vapors containing equal concentration. This result revealed a “reverse shape selectivity” phenomena, which states that the more bulky molecule is retained in favor the least bulky molecule. After this initial finding, others observed the reverse-shape selectivity phenomena for xylene isomers in gas and liquid chromatography experiments. For example, Chang and Yan used UiO-66-H as a stationary phase in capillary columns to separate a mixture xylene isomers and ethylbenzene [
22]. They found that
p-xylene eluted first, followed by
m-xylene, ethylbenzene, and
o-xylene, and calculated the adsorption enthalpy and entropy from the chromatograms. Subsequently, Duerinck et al. used inverse pulse gas chromatography to calculate the Henrys constants for BTEX compounds on a packed gas chromatography column of UiO-66-X where X was –H, –NO
2, and –CH
3 groups [
14]. In general, they found that decreasing the pore diameter had the largest influence on adsorption properties, rather than the specific nature of the functional group.
In addition to the experimental work, some theoretical studies have been done to better understand the reverse shape selectivity. The first attempt to theoretically study the adsorption of xylene isomers in UiO-66-H was performed by Granato et al. using grand-canonical Monte Carlo simulations [
23]. They used four different force fields (all models assumed a rigid UiO-66-H structure) to predict isotherms at 423 K for ethylbenzene and the xylene isomer. In general, the simulations confirmed the experimental ortho isomer selectivity. Subsequently, molecular dynamics simulations by Lennox and Düren were aimed at better understanding the diffusion mechanism of xylene isomers in UiO-66-H [
9]. They found that rotation of the benzene dicarboxylate linker was crucial to allow passage of xylene isomers from pore to pore.
In this work, static vapor-phase adsorption of BTEX compounds into UiO-66-X thin films was studied to better understand the effect of functional groups on the adsorption of aromatic compounds. Our experimental approach is different than the previous breakthrough measurements and chromatographic studies in two ways. First, we measure the adsorption of a single compound to the MOF and second, there is no competitive adsorption from a non-adsorbing solvent. In addition, this is the first time UiO-66-NH2 has been evaluated as a sorbent for BTEX. Lastly, by performing measurements near room temperature we are able to reach saturation adsorption capacity, which reveals insight about the role of missing linker defects.
2. Results and Discussion
X-ray diffraction was used to assess the crystallinity of the UiO-66 thin films.
Figure 1 shows the low angle Bragg peaks for the thin films. The most intense peak near 7.5° corresponds to the (111) plane. A weaker reflection near 8.5° is due to the (200) plane. In the UiO-66-H film, the (111) peak is the most intense and narrow peak, which indicates better crystallinity. In contrast, in the UiO-66-NO
2 film the (111) peak intensity is weaker and broader, which indicates less crystallinity. Previous work has shown the ratio of the intensity of the (111)/(200) peaks is a measure of missing linker defects [
24]. The higher the intensity of the ratio of the (111)/(200) peaks, the fewer missing linkers. Consequently, the data in
Figure 1 suggest that the UiO-66-H films have the fewest missing linkers and the UiO-66-NO
2 films have more missing linkers. Scanning electron microscopy was used to investigate the morphology and thickness of the films.
Figure 2 shows SEM images of the top view of the UiO-66-X films. In general, each film is continuous with evidence of grain boundaries. There are some regions (not shown) where microscopic cracks are observed, which are attributed to stress while drying. Cross-section SEM images reveal the films are about 1–2 μm thick.
Figure 3 shows IR spectra in the fingerprint region of the synthesized UiO-66-H, UiO-66-NH
2, and UiO-66-NO
2 films. The spectra are consistent with reports in literature for bulk powders [
4,
19]. There is a small peak around 1705–1655 cm
−1 that can be assigned to the carbonyl group of residual DMF solvent in the pores (which is removed upon activation before the adsorption studies). In all three RAIRS spectra, there are a set of strong peaks between 1580 and 1395 cm
−1 due to the symmetric and asymmetric carboxylate modes. There is also a strong C–H bending mode near 770 cm
−1 in all three films. There are several peaks unique to the amino and nitro functional groups. For example, in the UiO-66-NH
2 spectrum, the peaks at 1627 cm
−1 and 1260 cm
−1 are associated with the –NH
2 bending mode and C–N stretching of the amino group. In the case of the UiO-66-NO
2 spectrum, the peaks at 1540 cm
−1 and 1346 cm
−1 can be assigned to the asymmetric and symmetric stretching modes of the NO
2 group.
Figure 4 shows the Raman spectra from 650 to 1800 cm
−1. In the spectrum for UiO-66-H, there are several strong peaks: (i) the peak at 1615 cm
−1 can be assigned to the carbon–carbon bond stretching mode of the aromatic ring of the linker, (ii) at peaks at 1450 cm
−1 and 1430 cm
−1 are two modes associated with the in-phase carboxylate stretching normal mode, (iii) the peak at 1145 cm
−1 can be assigned to the breathing mode of the terephthalate ring and (iv) and a peak at 860 cm
−1 due to a C–H in-plane bending mode of the linker. In the UiO-66-NO
2 Raman spectrum, there is an additional peak at 1351 cm
−1 assigned to the symmetric stretching mode of the –NO
2 group, and, in the UiO-66-NH
2 film, there is a new peak at 1265 cm
−1 that can be assigned to the –NH
2 bending mode.
Figure 5 shows survey X-ray photoelectron spectra. The main peaks for the MOF nodes are the O(1s), Zr(3d) and Zr(3p) peaks as seen near 520 eV, 183 eV and 330 eV, respectively. The main peak for the linker is the C(1s) at 285 eV. The N(1s) binding energy for the NH
2 and NO
2 groups are observed at 404 and 400 eV, respectively. There is a small Au(4f) peak seen in the UiO-66-NO
2 spectrum, which is from the substrate. The absence of Cl suggests there is no unreacted starting material.
Figure 6 shows a representative set of QCM data on a UiO-66-NH
2 film. Similar data sets for the UiO-66-H and UiO-66-NO
2 films are shown in the
Supporting Information. In each experiment, the crystal frequency was recorded at each temperature until it stabilized. Multiple runs were measured after baking the chamber. The results of the multiple runs were averaged and used to calculate the Henry’s constant at each temperature (
Table S1 in Supplementary Materials). There are two global trends worth noting. First, the Henry’s constants on the UiO-66-H are an order of magnitude smaller than the other two MOF films. Second, the Henry’s constants for the UiO-66-NO
2 are greater than UiO-66-NH
2.
Our results indicate that larger substituents on the organic linker decrease the pore size and increase the interaction between the molecules and MOF framework. Previous nitrogen adsorption isotherms on UiO-66-X have shown a direct correlation between the functional group and the pore volume and surface area [
4]. For example, in UiO-66-H the surface area has been measured to be near 1300 m
2/g. In UiO-66-NH
2, the surface area decreases to 1250 m
2/g. When the amine is replaced by a nitro group in UiO-66-NO
2, the surface area decreases to 856 m
2/g. The porosity of the framework is sensitive to the synthesis conditions. In fact, the porosity can be tuned by the addition of monocarboxylic acids that act as modulators of the crystallization process [
25,
26,
27]. The monocarboxylic acids compete with binding to the nodes, which slows the crystallization kinetics and can increase the crystallite size. In addition, the monocarboxylic acid promotes the formation of missing linker defects. It is well-established that using a stronger monocarboxylic acid in the synthesis increases the surface area [
25]. The increase in surface area can be attributed to either missing linkers or missing nodes. It is very challenging to provide direct evidence for the missing linker defects. Initial work using thermal gravimetric analysis inferred the missing linkers based on stoichiometry arguments.
In 2013, Wu et al. used high-resolution powder neutron diffraction to provide the first direct structural evidence for missing linkers in the UiO-66 framework [
28]. Ideally, a perfect UiO-66-H framework should have pore volume and surface area of 0.426 cm
3/g and 954 m
2/g, respectively. If one of twelve linkers are removed, then the pore volume and surface area increase to 0.502 cm
3/g and 1433 m
2/g, respectively. Several other groups have established the importance of missing linkers in gas adsorption studies in UiO-66 [
24,
25,
26,
27,
28,
29,
30,
31,
32,
33,
34]. We have calculated the theoretical adsorption capacity (as weight percent) for each MOF assuming one BTEX molecule per pore with no missing linkers. In the ideal UiO-66-X unit cell there are four octahedral pores and eight tetrahedral pores. The Henry’s constants can be converted to a weight percent to compare with the theoretical adsorption capacity assuming no defects. From the experimental measurements at 30 °C we find that, on average, there is about one BTEX molecule per unit cell for the UiO-66-H MOF. In the UiO-66-NH
2 MOF, there are about 4–6 BTEX molecules per unit cell, and, within the UiO-66-NO
2 MOFs, there are about 9–11 BTEX molecules per unit cell.
Table 1 shows the measured adsorption capacities in UiO-66-NH
2 and UiO-66-NO
2 films exceed the theoretical adsorption capacity, which suggests that missing linkers are present. The adsorption capacity is calculated for each MOF assuming a defect free and missing-linker structure(s). In the case of the defect free unit cell, we divide the mass of 12 sorbent molecules by the mass of the perfect unit cell. In the case of one missing linker, we divide the mass of 13 sorbent molecules by the mass of the unit cell less one linker, etc. In the case of two missing linkers, we divide the mass of 14 sorbent molecules by the mass of the unit cell less one linker, etc. It can be seen the adsorption capacity for UiO-66-NH
2 and UiO-66-NO
2 is consistent with 1–2 and 3–4 missing linkers, respectively.
In general, the role of missing linkers in UiO-66 thin films has not been studied in detail. In fact, there are very few reports of UiO-66 films grown directly onto surfaces [
35,
36,
37,
38]. This is due, in part, to the fact that the common methods for quantifying missing linkers e.g., thermal gravimetric analysis, are not amenable to thin films. To the best of our knowledge, there are only two reports in the literature of UiO-66 films grown on QCM crystals. Virmani et al. measured the adsorption of ethanol at 25 °C into UiO-66-NH
2 films [
35]. Semrau et al. measured the adsorption of methanol into UiO-66-H at 40 °C [
36]. In the latter study, the authors used X-ray photoelectron spectroscopy to measure the elemental composition of films made using different amounts of a modulator. A comparison of the ideal, defect-free stoichiometry to the measured Zr, C, and O atomic composition can estimate the extent of defects in the films. Semrau et al. reported that films with more defects had an increased adsorption capacity of ethanol.
The role of the functional group is much more difficult to assess experimentally because it is nearly impossible to synthesize ideal materials. Several studies have reported higher adsorption capacity in UiO-66-NH
2 and UiO-66-NO
2 than UiO-66-H [
39,
40,
41,
42,
43,
44,
45]. In most cases, the smaller pore size of UiO-66-NO
2 and/or UiO-66-NH
2 compared to UiO-66-H is thought to contribute to stronger guest–host interactions. Demir et al. computationally screened 36 functional groups in the UiO-66 framework to assess the relative binding energies for NH
3, H
2O and CO
2. The calculations revealed that the dispersion energy at the UiO-66 pore center was similar for the NH
2 and NO
2 functional groups. In our study, it is most likely that pi–pi interactions between the aromatic molecules and organic linkers dominate the strength of adsorption, as seen in a recent study by Akpinar et al. [
46]. Overall, we conclude that the missing linkers play a more important role in our films than the nature of the functional group.