Artem S.
Poryvaev
a,
Aleksandr A.
Efremov
ab,
Dmitry V.
Alimov
ab,
Kristina A.
Smirnova
ab,
Daniil M.
Polyukhov
a,
Renad Z.
Sagdeev
a,
Samuel
Jacoutot
c,
Sylvain R. A.
Marque
*c and
Matvey V.
Fedin
*ab
aInternational Tomography Center SB RAS, Novosibirsk, 630090, Russia. E-mail: mfedin@tomo.nsc.ru
bNovosibirsk State University, Novosibirsk, 630090, Russia
cAix Marseille University, CNRS, UMR, Avenue Escadrille Normandie-Niemen, 7273 Marseille, 13397 CEDEX 20, France. E-mail: sylvain.marque@univ-amu.fr
First published on 1st March 2024
Metal–organic frameworks (MOFs) draw increasing attention as nanoenvironments for chemical reactions, especially in the field of catalysis. Knowing the specifics of MOF cavities is decisive in many of these cases; yet, obtaining them in situ remains very challenging. We report the first direct assessment of the apparent polarity and solvent organization inside MOF cavities using a dedicated structurally flexible spin probe. A stable β-phosphorylated nitroxide radical was incorporated into the cavities of a prospective MOF ZIF-8 in trace amounts. The spectroscopic properties of this probe depend on local polarity, structuredness, stiffness and cohesive pressure and can be precisely monitored by Electron Paramagnetic Resonance (EPR) spectroscopy. Using this approach, we have demonstrated experimentally that the cavities of bare ZIF-8 are sensed by guest molecules as highly non-polar inside. When various alcohols fill the cavities, remarkable self-organization of solvent molecules is observed leading to a higher apparent polarity in MOFs compared to the corresponding bulk alcohols. Accounting for such nanoorganization phenomena can be crucial for optimization of chemical reactions in MOFs, and the proposed methodology provides unique routes to study MOF cavities inside in situ, thus aiding in their various applications.
Electron Paramagnetic Resonance (EPR) is a powerful method for studying coordination compounds. Its applications to MOFs are currently expanding, mostly focusing on the properties of MOFs such as sorption,28–31 structural flexibility32–36 and catalysis.37–39 Since most of the MOFs are diamagnetic, the application of EPR requires incorporation of some paramagnetic species into the MOFs. Owing to the high sensitivity of EPR, this can be done in trace amounts (1 paramagnetic center per 1000 cavities or less), not disturbing the original structure of the framework. One of the ways to incorporate paramagnetic reporters into MOFs relies on postsynthetic sorption of paramagnetic probes such as 2,2,6,6-tetramethyl(piperidin-1-yl)oxyl (TEMPO),35,36,40,41 nitrogen monoxide32 or nitrogen dioxide.42 This approach is commonly called a method of spin probes. Previously, this method with a TEMPO probe was applied for monitoring gas adsorption on a ZIF-8 surface.41 Later on, a very fruitful modification has been proposed, called the method of encapsulated spin probes, where nitroxide (TEMPO) radicals were permanently entrapped in the cavities of a MOF (TEMPO@MOF) and used as reporters located exclusively in the cavities (not on the outer surface) of the MOF. This method allows numerous applications, including precise measurements of the window sizes vs. temperature,43 adjusting conditions for various separations.44,45 assessing pressure effects and aiding in the formation of MOF pellets and composites.46–50 All information obtained using TEMPO@MOF systems relied on the sensitivity of the EPR spectrum to the presence of oxygen or solvent molecules in the cavity, as well as to the integrity of MOF cavities. However, this TEMPO-based approach provides only basic information on the properties of MOF cavities inside and their changes upon adsorption of various guests.
In a recent series of studies, β-phosphorylated nitroxide radicals were successfully employed to study the local polarity and composition of various solutions.51–55 The hyperfine interaction (HFI) constants on nitrogen and phosphorus nuclei were found to be extremely sensitive to local polarity, structuredness, stiffness and cohesive pressure, and could be deduced from continuous wave (CW) EPR experiments in solution at ambient temperatures.
In this work, we combined these two independent research studies and incorporated a β-phosphorylated nitroxide into a ZIF-8 cavity as a probe of local polarity. The application of EPR to this system allowed the first direct assessment of the polarity inside the ZIF-8 cavity and its guest-induced changes. Moreover, when ZIF-8 was impregnated with alcohols, the polarity sensed by the probe was found anomalous, in some cases exceeding the bulk polarity of the same solvent, to be assigned to specific guest–host interactions of alcohols with the ZIF-8 inner surface. Below we describe the developed methodology, obtained results and important conclusions obtained from these studies.
The conformation of this radical and, correspondingly, observed HFI splitting strongly depend on the polarity of surrounding media and characteristics such as structuredness, stiffness and cohesive pressure of the surrounding solvent.51–55
First we have ensured that the radicals in the synthesized S1@ZIF-8 material are encapsulated (permanently entrapped) inside the cavities of the MOF. For this sake we confirmed the formation of the ZIF-8 material using PXRD (Fig. S1, ESI†) and studied radical behavior by EPR (Fig. 2). The radicals in ZIF-8 demonstrate fast motion EPR spectra characteristics of unrestricted rotation in the cavity (Fig. 2). These spectra remain unchanged upon evacuation with heating, as well as upon flushing with ethanol. Consequently, once these radicals entered the cavities during synthesis, they cannot exit outside due to their large size compared to inter-cavity windows. The room-temperature EPR spectrum of S1@ZIF-8 features a single fraction with narrow lines, meaning that the radical is evenly distributed over the volume of the MOF.
Fig. 2 compares the CW EPR spectra of S1 dissolved in water and S1@ZIF-8. The simulation of experimental data shown in Fig. 2 allows obtaining the following set of isotropic HFI constants on nitrogen (AN) and β-phosphorus (AP) nuclei (Table 1). The HFIs of this radical were obtained previously in a large set of solvents.51 Our results in water (AN = 1.612 mT and AP = 2.231 mT) agree with previous data within an accuracy of ±0.04 and ±0.05 mT for AN and AP, respectively (these values can serve as estimates of measurement error). The spectra in water and in ZIF-8 are markedly different, owing to the different values of AN and AP. Both HFIs of these radicals are sensitive to the polarity of the media,51–55 and the changes of their values between water and the ZIF-8 interior indicate stronger hydrophobicity of the latter. Previous studies have also demonstrated that AN and AP show opposite trends vs. normalized Reichard polarity constants (ENT): AN moderately grows with ENT, whereas AP strongly decreases upon polarity increase. The obtained HFIs for S1@ZIF-8 (AN = 1.420 mT and AP = 4.033 mT) are very close to the values obtained for S1 in n-pentane (AN = 1.430 mT and AP = 4.100 mT), meaning that the polarity sensed by a probe inside bare ZIF-8 is approximately as small as in n-pentane. Interestingly, however, that both AN and AP values are slightly lower in ZIF-8 than in n-pentane, while one would expect differences of the opposite sign. This might indicate a subtle confinement effect on the conformation of S1 in ZIF-8 that ultimately influences the observed HFIs. However, the magnitude of the difference observed between S1@ZIF and S1 in n-pentane is comparable to the accuracy of measurements. Thus, we provide the first experimental evidence that the internal surface of the ZIF-8 cavity sensed by guest molecules is strongly non-polar, in perfect agreement with theoretical expectations for this structure and with the analysis of water vapor sorption isotherms.23,27
Environment | A N/mT ±0.040 mT | A P/mT ±0.050 mT | E N T | θ | V X | V M | c | No. |
---|---|---|---|---|---|---|---|---|
Water | 1.612 | 2.231 | 16.7 | |||||
1.650 (ref. 52) | 2.280 (ref. 52) | 1.000 (ref. 51) | 54 | 18 | 2294 | |||
ZIF-8 | 1.420 | 4.033 | 34 | 1 | ||||
n-Pentane | 1.406 (ref. 51) | 4.135 (ref. 51) | 0.009 (ref. 51) | 34 | 81.3 | 114.52 | 205 | |
ZIF-8 + isopropanol | 1.514 | 2.271 | 53 | 2 | ||||
Isopropanol | 1.520 (ref. 61) | 2.760 (ref. 61) | 0.546 (ref. 51) | 50 | 59 | 76.51 | 558 | |
ZIF-8 + propanol-1 | 1.520 | 2.038 | 55 | 3 | ||||
Propanol-1 | 1.525 | 2.453 | 0.617 (ref. 51) | 51 | 59 | 74.8 | 600 | |
ZIF-8 + ethanol | 1.535 | 1.840 | 58 | 4 | ||||
Ethanol | 1.560 (ref. 61) | 2.400 (ref. 61) | 0.654 (ref. 51) | 52 | 44.9 | 58.41 | 676 | |
ZIF-8 + methanol | 1.543 | 1.543 | 61 | 5 | ||||
Methanol | 1.570 (ref. 61) | 2.190 (ref. 61) | 0.762 (ref. 51) | 55 | 30.8 | 40.43 | 858 | |
ZIF-8 + ethylene glycol | 1.565 | 1.470 | 62 | 6 | ||||
Ethylene glycol | 1.540 (ref. 51) | 2.159 (ref. 51) | 0.790 (ref. 51) | 54 | 50.8 | 56.01 | 1050 | |
ZIF-8 + TFE | 1.560 | 1.430 | 62 | 7 | ||||
TFE | 1.590 (ref. 51) | 2.042 (ref. 51) | 0.898 (ref. 51) | 56 | 41.5 | 72.4 | 573 |
As an example, Fig. 3 shows CW EPR of S1@ZIF-8 impregnated with different alcohols. The empty cavity of ZIF-8 can accommodate 8–10 molecules of alcohols; however, in the presence of nitroxide in the same cavity, this number can be reduced to ∼4–5. The spectra are markedly different, and the obtained HFI constants are listed in Table 1.
However, the most astonishing observation is that apparent polarity sensed by the spin probe inside ZIF-8 impregnated with alcohols differs from that in the corresponding bulk alcohols in a peculiar way. Fig. 4a shows the AN and AP values vs. Reichardt polarity ENT for ZIF-8 impregnated with a set of solvents, as well as for the same bulk solvents. The AN values change very insignificantly with ENT, and there is also no solid trend between AN in the impregnated MOF and in bulk solvent. At the same time, the AP values change noticeably vs. ENT, and they are markedly different in the MOF vs. bulk. Note that the other non-hydrogen bonding solvents show the same trend (see more data in the ESI†).
![]() | ||
Fig. 4 (a) AN and AP values vs. Reichardt polarity ENT in the bulk and in ZIF-8. (b) AP/AN ratio vs. ENT in the bulk (red) and in ZIF-8 (blue). The lines guide the eye. |
Fig. 4a can be understood in a way that the polarity sensed via AP is much higher in impregnated ZIF-8 compared to the same bulk solvent. At the same time, AN shows the same trend only for some solvents, and the results based on AN are inconclusive. In fact, they fall within the uncertainty introduced by a minor confinement effect discussed above on the basis of comparison with n-pentane (∼0.06 mT). Thus, we should hold to the more reliable data delivered by AP, whose dependence on solvent polarity spreads over a much broader range.
The electron density on the phosphorus atom originates from the hyperconjugation mechanism, which depends on the angle θ shown in Fig. 1cvia the Heller–McConnell relationship:
AP = B0 + B1ρπN![]() ![]() ![]() ![]() | (1) |
Fig. 4b shows the ratio AP/ANvs. ENT, which can be understood using eqn (1) as indirect dependence of θ on ENT. As is shown in Fig. 4b, the dependence of AN and AP in the MOF is more pronounced than in the bulk, implying that angle θ changes more vs. solvent polarity for radical encapsulated in the MOF.
One might think that the confinement of radicals in the MOF filled with the solvent imposes steric hindrance leading to specific distortions of radical geometry and amplification of angle θ in the MOF compared to the corresponding bulk solvent. However, if pronounced steric hindrances take place, one would simultaneously expect more restricted and potentially anisotropic mobility of the radical. However, all EPR spectra shown in Fig. 3 have the lineshapes close to the isotropic limit, which implies that the rotation of the radical is rather fast and thus not strongly constrained. Recently, we have shown that TEMPO nitroxide exhibits slightly anisotropic mobility at room temperatures being confined in ZIF-8 impregnated with alcohols.62 However, despite some anisotropy of rotation imposed by the cavity partly occupied by solvent molecules, the observed rotational correlation times were rather short (<130 ps). Thus, as long as the single set of HFI couplings is observed, one can assume that the values of these couplings correctly reflect the apparent polarity of the MOF cavity.
The most evident difference between bulk and in-MOF AP values is found for methanol, where we observe AN = AP = 1.543 mT for S1@ZIF-8 and [AN, AP] = [1.570, 2.190] mT for S1 in bulk methanol. The nitrogen HFI AN is nearly the same (within given accuracy) in both cases, whereas phosphorus HFI AP is dramatically smaller in the MOF. In fact, such a small value AP = 1.543 mT was not obtained in any bulk organic solvent, even in water.63,64
The decrease in Ap (Fig. 4a) is dependent on the dihedral angle θ (eqn (1)) and denotes that conformer B is favored over conformer A in a polar environment due to its N+˙—O–⋯P+—O– dipole–dipole interaction (see Fig. 1c). We suggest there are two plausible mechanisms that can explain larger apparent polarity (larger θ) obtained in ZIF-8 impregnated with solvents compared to the same bulk solvents. Both mechanisms assume that the ZIF-8 cavity introduces nanoordering of solvent molecules, which preferentially cover the walls of the cavities (the latter is well justified by rather fast and isotropic rotation of the radical, as was discussed above).
First the ‘structuredness-driven’ mechanism involves interplay between the polarity and factors such as solvent organization given by the molar volume VM. The reorganization of solvent forced by the ZIF cavity leads to a striking increase of the organization of the cybotactic layer depending on the size of the solvent. Therefore, solvent molecules covering the cavity walls (higher structuredness) cause the release of constraints around the C–P bond of a radical compared to the bulk solvent, thus affording easier rotation (lower stiffness) and favoring conformer B (larger θ). This loss of organization is balanced with the increase of the size, which reduces the space in the cavity (as given by the negative value of the parameter, see below). Once there are polar solvent molecules in the cavity of ZIF-8, they impose electric fields on a radical (favoring conformer B), but at the same time do not provide as strong steric constraints as the same molecules in the bulk, resulting in larger θ values compared to those in the same bulk solvent.
The second ‘hyperpolarity’ mechanism relies on the fact that alcohols, for example methanol, are composed of a hydrophobic part (alkyl fragment) and a hydrophilic part (OH function). Therefore, it is plausible that guest molecules inside the cavity of a MOF are nanoordered in a way that non-polar (hydrophobic) moieties are directed towards the cavity walls, while polar moieties are directed towards the cavity center (sketched in Fig. 5). In this case, the electric field imposed onto a radical in the MOF impregnated with solvent (e.g. methanol) can be higher than that in the bulk, where all solvent molecules around the radical are disordered, thus leading to a higher apparent polarity (‘hyperpolarity’) sensed by the radical. We attempted to confirm this mechanism using pulse EPR (ESEEM) and partly deuterated methanol; however, the observed trend was too weak to provide ultimate experimental proof (see the ESI†).
![]() | ||
Fig. 5 Proposed methanol localization in the cavity of S1@ZIF-8. Hydrophilic environment around the probe is shaded in red, whereas the hydrophobic part of the cavity is marked in gray. |
In fact, both mechanisms can work simultaneously; however, let us briefly discuss the pros and cons of each one in view of available experimental data.
When other alcohols are used (ethanol, isopropanol, ethylene glycol, and TFE), the apparent polarity inside the cavity sensed by a phosphorus HFI (AP) progressively decreases (AP increases) (Table 1). This trend follows the one for bulk solvents; but still the apparent polarity in the cavity is systematically larger than that in the bulk solvent. The decrease of the apparent polarity in the order methanol > ethanol > isopropanol agrees well with the qualitative expectations that the ratio of polar/non-polar moieties in these molecules decreases. In this case, even though ordering of solvent molecules in the MOF is still effective (leading to a higher polarity vs. the corresponding bulk solvent), the apparent polarity can be comparable to that in some common organic solvents (though different from the impregnating one). At the same time, the volume of solvent molecules grows in the order isopropanol > ethanol > methanol, implying possible increase of steric hindrances for the rotation around the C–N bond and inhibiting large θ-values. Therefore, the trends observed for a series of alcohols are consistent with both structuredness-driven and hyperpolarity mechanisms.
Let us consider in more detail the absolute values of the observed spectroscopic parameters. Assuming B1ρπN = 5.9 mT in non-polar solvent such as pentane,51–55 the corresponding θ value can be estimated as 34° based on eqn (1), both in bulk pentane and in solvent-free ZIF-8. To determine θ in other solvents, changes in ρπN due to the change in the polarity of solvent are accounted by assuming ρπN ≈ AN. Thus, θ is given by eqn (2) (AP,pentane and AN,pentane refer to the 31P and 14N HFI constants in bulk n-pentane, respectively, whereas AP,n and AN,n stand for similar values in the arbitrary case; θpentane is given by eqn (1)).
![]() | (2) |
Previously, it has been shown that AN and AP values can be described using multiparameter Kamlet–Abboud–Taft and Koppel–Palm equations, which take into account polarity, H-bonding, bulk/cavity terms given by the structuredness/stiffness, and volume of solvent (c, VX or VM, respectively).52–55 Due to the small number of solvents investigated (see Table 1), the KAT relationship is preferred in the form where polarity and H-bonding are accounted in Reichardt constant ETN, reducing the number of parameters. It is combined, in eqn (3), with the molar volume VM, which provides the best correlation (see the ESI†) describing the bulk/cavity effect (Fig. 6 and eqn (4)): ETN (N = 11, t-test at 99.99% for θ0, at 99.88% for ETN, and at 99.52% for VM, R2 = 0.976, F = 165, and F-test at 99.99%).
θ = θ0 + a1ETN + a4VM | (3) |
θ = 53.1(±1.7) + 19.4(±2.8)·ETN − 0.13(±0.03)·VM | (4) |
![]() | ||
Fig. 6 Plot of dihedral angle θ vs. f(ENT,VM) given by eqn (4). |
Fig. 6 shows that the combination of the Reichardt polarity constant ETN and cohesive pressure VM allows for good linear correlation with the observed θ values. This means that the observed anomalous HFI constants in ZIF-8 impregnated with solvents stem from both structuredness-driven and hyperpolarity mechanisms. Their relative contributions are different for different solvents. Although the detailed analysis is still to be done in the future, based on the values reported in Table 1 we conclude that the second mechanism dominates for the studied set of solvents. It is noteworthy that both mechanisms rely on the solvent ordering in the cavities, which favors the B conformer much stronger in-MOF compared to the same bulk solvents.
To confirm that the rotation around the C–N bond (Fig. 1c) is the major factor influencing AP values, we performed a series of DFT calculations (see the ESI†). We calculated the dependence of AP on cos2θ upon variations of other geometric parameters (bond lengths, angles, and dihedrals). In all cases, the correlation between AP and cos2
θ was maintained, and the deviations were smaller than the changes of APvs. θ.
Remarkably, the observed solvent ordering is an important effect for carrying out the reactions in ZIF-8, since the structure and properties (such as polarity) of the active site environment might be decisive for the enhancement/inhibition of particular reaction routes. Moreover, polarity-dependent stabilization of certain transition states has high impact on the efficiency of enzymatic reactions and selective catalysts.65–69 Therefore, we believe that this finding is extremely important for understanding chemical reactions occurring in MOFs and must be taken into consideration when such particular applications are being designed.
In both cases we have obtained very close values for minor components (n-pentane: AN = 1.431 mT and AP = 3.650 mT; n-hexane: AN = 1.431 mT and AP = 3.798 mT). Moreover these values of AN and AP are close to those observed in the corresponding non-polar bulk solvents. At the same time, the sets of parameters for the major components have noticeably different values of AP in n-hexane and n-pentane, while AN values are closely the same. Therefore, we speculate that for each of two above solvents there are two local surroundings of the radical in the ZIF-8 cavity: one that does not impose steric hindrances leading to enforced conformation and another one that does impose such hindrances. Most plausibly, the latter situation refers to the cavities with a larger number of solvent molecules inside compared to the former one. In addition, when we deal with a ‘tighter’ radical surrounding, the resulting enforced conformation should naturally depend on the structure of solvent. This explains why the second (major, enforced) components of the spectra are described by different HFIs for n-pentane and n-hexane, while the first ones (minor, unperturbed) are rather close.
In summary, the above consideration shows that the minor components of these complex two-component EPR spectra reflect the polarity and structuredness of the radical environment, whereas the major components refer to ‘unnatural’ radical conformations enforced by severe steric hindrances. In principle, in a systematic study the component of interest can be safely identified and interpreted.
All these benefits were demonstrated using one of the most appealing MOFs nowadays – ZIF-8. Exploiting the radical sensitive to a microenvironment we have experimentally evidenced that the ZIF-8 cavity is non-polar. Also, we determined that impregnation of ZIF-8 with different alcohols leads to a drastic change in the apparent polarity inside the cavity (i.e. a polarity sensed by solutes) from non-polar to a highly polar one. Remarkably, in the cases of methanol, ethanol and isopropanol, the apparent polarity inside ZIF-8 was found to be higher than that in the bulk solvent. In the cases of methanol and ethanol, it is even higher than that in any bulk organic solvent. This was rationalized by nanoordering of alcohol molecules inside the cavity of ZIF-8 driven by hydrophobic interactions.
Many applications of MOFs require their impregnation with a solution of reactants, which is especially true for the MOF catalysis.70–74 Moreover, the design of chemical processes in MOFs needs understanding of local properties, in particular the polarities of MOF internal surfaces sensed by solutes and carrier solvents, as well as the structuredness of the corresponding cybotactic layers. As is shown in the present work, the values of apparent polarity inside MOFs and local solvent organization can be drastically different from expectations based on the bulk properties of the same molecules, and this needs to be carefully accounted for in the future. The methodology developed by us has a broad scope of applications, and it is quite easy to implement and provides a lot of opportunities for understanding and rational optimization of chemical reactions in MOFs. In this work it was exemplified using ZIF-8; however, the properties of many other MOFs (for instance, other MOFs of the ZIF family and recent example of UiO-66 (ref. 50)) can be accessed in the same way. Therefore, we believe that the methodology and conclusions of this paper may find broad implementation in MOF science in the future.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3sc05724k |
This journal is © The Royal Society of Chemistry 2024 |