Helge
Reinsch
*a,
Renjith S.
Pillai
b,
Renée
Siegel
c,
Jürgen
Senker
c,
Alexandra
Lieb
d,
Guillaume
Maurin
b and
Norbert
Stock
*e
aDepartment of Chemistry, University of Oslo, P.O. Box 1033 Blindern, 0315 Oslo, Norway. E-mail: helge.reinsch@smn.uio.no
bInstitut Charles Gerhard Montpellier, UMR-5253 Université Montpellier CNRS ENSCM, Place E. Bataillon 34095, Montpellier cedex 05, France. E-mail: guillaume.maurin@univ-montp2.fr
cInorganic Chemistry III, University Bayreuth, Universitätsstaße. 30, 95447 Bayreuth, Germany
dInstitute of Chemistry, Universitätsplatz 2, 39106, Magdeburg, Germany
eInstitute of Inorganic Chemistry, Christian-Albrechts-Universität Kiel, Max-Eyth-Straße 2, 24118 Kiel, Germany. E-mail: stock@ac.uni-kiel.de
First published on 16th October 2015
The new aluminium based metal–organic framework [Al(OH)(O2C–C4H8–CO2)]·H2O denoted as Al-MIL-53-ADP-lp (lp stands for large pore) was synthesised under solvothermal conditions. This solid is an analogue of the archetypical aluminium terephthalate Al-MIL-53 based on the aliphatic single-chain linker molecule adipic acid (H2ADP, hexanedioic acid). In contrast to its aromatic counterparts, Al-MIL-53-ADP exhibits a structural breathing behaviour solely upon dehydration/rehydration. The crystal structure of the anhydrous compound denoted as Al-MIL-53-ADP-np (np stands for narrow pore) was determined by a combination of forcefield-based computations and Rietveld refinement of the powder X-ray diffraction data while the structure of the hydrated form Al-MIL-53-ADP-lp was derived computationally by a combination of force field based methods and Density Functional Theory calculations. Both structures were further supported by 1H, 13C and 27Al high-resolution NMR MAS 1D data coupled again with simulations. Al-MIL-53-ADP was further characterised by means of vibrational spectroscopy, elemental analysis, thermogravimetry and water vapour sorption.
Herein, we report the aluminium based metal–organic framework [Al(OH)(O2C–C4H8–CO2)] or Al-MIL-53-ADP (ADP stands for adipate) and its breathing behaviour, which is the first example of a MIL-53 framework based on a single-chain saturated aliphatic linker molecule.
All solid-state NMR experiments were acquired on Bruker Avance-III HD spectrometers operating at a B0 field of 14.1 T (1H and 27Al) and 9.4 T (13C). 1H (600.15 MHz) high-resolution spectra were acquired after a 90° pulse of 1.4 μs with a spinning speed of 62.5 kHz using a commercial 1.3 mm MAS double resonance probe (Bruker). 13C (100.6 MHz) MAS spectra were obtained with a ramped cross-polarization (CP) experiment where the nutation frequency νnut on the proton channel was varied linearly by 50%. The samples were spun at 12.5 kHz in a 4 mm MAS double resonance probe (Bruker). The corresponding νnut on the 13C channel and the contact time were adjusted to 70 kHz and 3.0 ms, respectively. Proton broadband decoupling with spinal-64 was applied during acquisition. Both the 1H and 13C spectra are referenced with respect to TMS (tetramethylsilane) using the secondary standard adamantane. Finally, the 27Al (156.4 MHz) 1D MAS experiments were acquired using a rotor-synchronized spin-echo with central transition selective 90° and 180° pulses of 12.5 and 25 μs (to remove the resonances from the satellite transitions) at a spinning speed of 20 kHz (3.2 mm MAS triple resonance probe, Bruker). The reference was set using a solution of AlCl3 adjusted to a pH of one with hydrochloric acid. The 27Al spectra were fitted using the simulation package SOLA included in Bruker's Topsin3.2 software.
The force field based-geometry optimizations were conducted using the Forcite software available in the Materials Studio package and the generic Universal Force Field Lennard–Jones parameters for all atoms of the MOF framework.29
The DFT geometry optimizations were performed with the QUICKSTEP module, available as part of the CP2K code using the PBE functional and the triple zeta basis set (TZVP-MOLOPT) for all atoms, except for the Al centres, where double zeta functions (DZVP-MOLOPT) were considered. Semi-empirical dispersion corrections as implemented in the DFT-D3 method were taken into account.30 These calculations considered simulation boxes consisting of (1 × 2 × 2) and (2 × 2 × 2) unit cells for the anhydrous and hydrated phases respectively. The first-principles calculations of the NMR parameters were performed from the structural models for both the large and the narrow pore forms with CASTEP31,32 using the PBE functional33 and ultra-soft pseudo-potentials generated ‘on the fly’.34 The wave functions were expanded on a plane wave basis set with a kinetic energy cut-off of 610 eV. The projector augmented waves (PAW)35 and gauge included projector augmented waves (GIPAW)36 algorithms for the EFGs and NMR chemical shifts, respectively.
Ideal formula | [Al(OH)(O2C–C4H8–CO2)] | [Al(OH)(O2C–C4H8–CO2)]·H2O |
Crystal system | Monoclinic | Triclinic |
Space group | P21 | P1 |
a/Å | 15.8741(8) | 6.6289 |
b/Å | 7.9646(3) | 8.6458 |
c/Å | 6.6187(3) | 9.0811 |
α/° | 90 | 62.429 |
β/° | 93.669(5) | 87.743 |
γ/° | 90 | 88.372 |
V/Å3 | 835.1(3) | 461.0 |
Method | Rietveld refinement | DFT calculation |
R wp/% | 8.1 | — |
R Bragg/% | 1.5 | — |
GoF | 2.2 | — |
For the hydrated form Al-MIL-53-ADP-lp we started with the indexed triclinic cell parameters (a = 6.623(2) Å, b = 8.624(2) Å, c = 9.058(2) Å, α = 117.38(2)°, β = 92.54(1)°, γ = 88.27(2)°) which were obtained after a structure-less cell parameter refinement.
These lattice parameters are similar to the ones reported for CAU-13,25 which contains cyclohexanedicarboxylate ions and also exhibits the MIL-53-type structure. Hence, a ligand replacement strategy was applied to the structure of CAU-13 while the space group was converted from P to P1. This structural model without water molecules was then geometry optimized at the force field level by imposing the cell parameters to the indexed values. The resulting structural model could not be refined by Rietveld methods due to a lower crystallinity of this form compared to the anhydrous solid. Hence, this structural model was completed by including one water molecule per formula unit as determined from TGA and water sorption experiments and subsequently optimised by DFT calculations, retaining the indexed cell parameters and taking also the hydrogen atoms into account in order to get a more reliable model. Such an approach has been proven highly valuable for the structure elucidation of complex and weakly ordered materials.22,39 The final optimised model was converted to conventional cell parameter settings (Table 1) and led to a simulated PXRD pattern in reasonable agreement with that observed (Fig. S1†). However, the intensities especially at low 2θ values are much lower than expected for a highly crystalline structure. We attribute this to severe disorder of the H2O guest molecules as well as to the relatively easy conformational changes of the aliphatic linker molecules induced by small rotational movements which require small activation barriers to be overcome. Accounting for this by using the spherical harmonics corrections for an initial fit without any structural refinement gives a reasonable agreement of the calculated and simulated PXRD pattern (see ESI†). The spherical harmonics correction is usually applied for modelling of deviating intensities in case of strong textural effects, however, in this case we attribute these deviations to the pronounced disorder in Al-MIL-53-ADP-lp.
For compounds based on aromatic linker molecules, the channels exhibit rectangular shape and opposite channel walls that are strictly parallel to each other. The shape of the channels in Al-MIL-53-ADP-np rather resembles a “Z” and thus the observed packing does not show parallel opposite channel walls (Fig. 3).
Fig. 3 The structure of Al-MIL-53-ADP-np as seen along the c axis, showing the packing of the “Z”-shaped channels. AlO polyhedra in grey and carbon atoms in black. |
The crystal structure of Al-MIL-53-ADP-lp does not show any change in connectivity but represents a hydrated analogue of the narrow pore form. Due to the insertion of water molecules the channels open up to a size suitable to accommodate water molecules. Thus the channel diameter as determined by pore size distribution analysis carried out on the structural model is ca. 3.2 Å (Fig. 4 and S3†). The shape of the channels rather resembles a kinked variant of the normal rhombic form.
Fig. 4 One channel in the structure of Al-MIL-53-ADP-lp as seen along the a axis. Hydrogen atoms are omitted for clarity. Aluminium atoms in grey, oxygen atoms in red and carbon atoms in black. |
The hydrated structure reveals that the water molecules are arranged in such a way to form hydrogen bonds with both the μ2-(OH) groups and one oxygen atom of one carboxylate group with characteristic distances between two water/host oxygen atoms ranging from 2.88 to 3.04 Å. In addition to this, the confined water molecules form a relatively strong hydrogen bond network with corresponding Ow_Ow distances (2.92 Å). This spatial distribution and the relevant interactions are very similar to that previously pointed in the MIL-53 series.22,40 Remarkably, the adsorption of water induces a rotation of the AlO-polyhedra chains and a framework expansion of ≈10% (vs. 3–4% and 33% for Al-MIL-53-FUM22 and Al-MIL-53-BDC13 respectively) while the length of the linker molecules is decreased due their conformational changes. The distances between carboxylate groups in the independent linker molecules in Al-MIL-53-ADP-np are 5.79(2) and 5.90(2) Å while they are 5.43 and 5.7 Å in the hydrated form Al-MIL-53-ADP-lp. This “shortening” of the linker molecules could also allow for an easier rotation around the C–C bonds and therefore might be the reason for the observed disordering in the crystal structure. A comparison with the well known highly crystalline layered compound zinc adipate41 pronounces the conformationally compressed character of the linker molecules. Zinc adipate incorporates adipate ions in their longest possible conformation and therefore the average distance between carboxylate groups amounts to 6.432(4) Å (Fig. 5).
Fig. 5 Conformation and lengths (carboxylate to carboxylate) of the adipate molecules in Al-MIL-53-ADP and for comparison also in zinc adipate.41 |
This change in conformation also seems to be the dominant mechanism controlling the breathing behaviour. In the aromatic parent compound Al-MIL-5313 the knee-cap mechanism results in profound changes in the angle between adjacently coordinated carboxylate groups. These values amount to 100° and 128° in the large pore and the narrow pore forms of the Al-MIL53-BDC solids, respectively. Here the scenario is different: for the hydrated Al-MIL-53-ADP-lp, the angles between crystallographically independent groups are 97 and 103° and upon dehydration to the narrow pore form they only slightly change to 97 and 107° (Fig. S4†).
It is worth mentioning that we also tried to intercalate several other solvent molecules like ethanol, acetone, toluene or dimethylformamide by dispersion of the dried MOF at elevated temperature, but all attempts proved unsuccessful.
The thermogravimetric curve (Fig. 7) was measured for a sample prepared directly after filtration from an aqueous dispersion and shows two distinct weight losses. The first step is due to the dehydration of the structure (10 wt%) while the second step of 63.4% above 360 °C is attributed to the decomposition of the framework. The first weight loss of 10 wt% is slightly larger than the expected value (9 wt%) while the second step (63.4%) is slightly smaller than the expected loss (66.5%) based on the deduced composition [Al(OH)(O2C–C4H8–CO2)]·H2O. Taking into account that there are still sulphur containing impurities present, both values are in reasonable agreement with the expected ones.
In order to investigate the sorption properties of Al-MIL-53-ADP, several sorption isotherms were measured. However, no uptake was observed for N2 at −196 °C or for CO2 at 25 °C, regardless of activation conditions (0.1 mbar, between 100 °C and 200 °C). Only for water vapour we observed a considerable adsorption. The isotherm is shown in Fig. 8 and exhibits a strong hysteresis loop. By measuring PXRD patterns depending on relative humidity we were able to identify the respective framework conformation during the sorption process (Fig. 8).
The adsorption of water at 25 °C starts at a relative humidity of ≈55% and proceeds in one step. The adsorbed amount of 1.1 H2O per aluminium atom is in good agreement with the expected amount of 1 molecule. Compared to the adsorption step at ≈55% relative humidity, the desorption starts at a much lower value of ≈25%. Such hysteretic water sorption isotherms are often observed in breathing MOFs, for example for Al-MIL-53-(OH)217 or Al-MIL-53-(1,4-NDC).43 The PXRD data indicate a rather slow equilibration during the uptake and release of H2O molecules. Diffractograms measured at 60% relative humidity during uptake show in general both forms, Al-MIL-53-ADP-np and Al-MIL-53-ADP-lp.
Nearly all of the dry form is converted after 20 minutes of equilibration at 70% relative humidity. The same mixture of framework conformers is observed for the desorption curve at 20% relative humidity. Moreover, even after extended equilibration under very moist (relative humidity >70%) or very dry (relative humidity <20%) conditions, there are still traces of the respective counterpart conformation present. A possible explanation could be that the diffusion into and out of the very narrow pores in Al-MIL-53-ADP is very slow and thus impedes fast equilibration times below 30 minutes. In addition the MOF crystals show a strong tendency for intergrowth, indicated by the formation of a thick layer of product at the glass reactor walls during the synthesis. This might also slow down the interparticle diffusion between intergrown crystals and thus delay equilibration.
The 13C spectra (Fig. 10) show three well-separated peak regions characteristic of the three different carbon types in the structures. While the carboxyl units are observed between 179 and 182 ppm, the outer and inner pairs of CH2 groups resonate in the ranges 37.0–38.2 ppm and 23.7–25.8 ppm, respectively. In the case of Al-MIL-53-ADP-np, six resonances – each two for the inner and outer CH2 groups (23.7/25.8 ppm and 37.1/38.1 ppm) as well as for the CO2 units (179.7/179.8 ppm)) – were observed. These observations are again in very good agreement with our DFT predicted 13C NMR parameters for both CH2 groups (22.5/25.2 ppm and 36.4/37.2 ppm) and CO2 units (183.4 ppm), this accounts for half of the possible resonances with respect to the structural model. Since the two adipate ions in the asymmetric unit are almost inversion symmetric we expect that the intramolecular splitting between resonances of similar chemical units of the same adipate molecule could not be resolved.
Fig. 10 13C CP MAS NMR spectra of Al-MIL-53-ADP. The resonances marked with an asterisk represent a small impurity of the np form in the lp phase. |
Due to the pronounced dynamical disorder of the CH2 units of adipate molecules and the water molecules the shift differences between similar inner and outer CH2 units, respectively, are averaged for Al-MIL-53-ADP-lp. This leads to only two resonances at 24.9 and 37.6 ppm in this region. In contrast, for the carboxylic groups of Al-MIL-53-ADP-lp two well resolved resonances at 180.8 and 181.4 ppm with a splitting of 0.6 ppm are visible. These experimental data are consistent with the theoretical findings reported in Table S2.†
The 27Al NMR spectra of both phases are displayed in Fig. 11 along with their simulated spectra including chemical shift and second order quadrupolar interactions (Table 2). For both the np and the lp forms two individual resonances were observed which is in line with the structure solutions of the PXRD data. All resonances exhibit isotropic chemical shifts close to zero reflecting the octahedral coordination of the Al3+ ions. The significantly different quadrupolar coupling constants CQ and asymmetry parameters η for the np phase demonstrate that the local environment of both Al octahedra differs markedly (Table 2). The CQ (9.06 MHz, 7.96 MHz) obtained for the anhydrous form are consistent with the experimental findings previously reported for the dry Al-MIL-53-BDC (CQ = 8.3 MHz) and Al-MIL-53-FUM (8.7 MHz) materials.22 And they also match with the corresponding first-principle calculated NMR parameters (8.60 MHz, 7.90 MHz).22 The quadrupolar interactions of the Al nuclei in the lp phase turn out to be very similar (8.32 and 8.37 MHz see Table 2) and could almost be assimilated into one single resonance. This is in contrast to the first-principles calculations where two types of Al3+, with CQ values of 8.3 MHz and 12.5 MHz respectively have been found (see Table S2†). This deviation probably results from a dynamical disorder, which was experimentally observed for this phase.
Fig. 11 27Al MAS NMR spectra of Al-MIL-53-ADP along with their respective fitted spectra (— individual resonances, -- sum of all resonances). |
[Al(OH)(O2C–C4H8–CO2)] | [Al(OH)(O2C–C4H8–CO2)]·H2O | |
---|---|---|
δ iso/ppm | 0.2/−0.5 | 0.2/0.3 |
C Q/MHz | 9.06/7.96 | 8.32/8.37 |
η | 0.0/0.19 | 0.0/0.13 |
Overall, the NMR experimental findings are in agreement with the corresponding first principles calculated parameters determined from the resolved crystal structures and support the reliability of the proposed structural models.
Footnote |
† Electronic supplementary information (ESI) available: Additional XRPD-data, simulated PSD, crystallographic information files, simulated NMR data. See DOI: 10.1039/c5dt03510d |
This journal is © The Royal Society of Chemistry 2016 |