Elena
Sanna
a,
Eduardo C.
Escudero-Adán
bc,
Antonio
Bauzá
a,
Pablo
Ballester
bc,
Antonio
Frontera
a,
Carmen
Rotger
a and
Antonio
Costa
*a
aDepartament de Química, Universitat de les Illes Balears, Palma, 07122, Spain. E-mail: antoni.costa@uib.es; Fax: +34 971 172436; Tel: +34 971 173266
bInstitut of Chemical Research of Catalonia (ICIQ), Avda. Països Catalans 16, 43007 Tarragona, Spain
cCatalan Institution for Research and Advanced Studies (ICREA), Passeig Lluis Companys 23, 08010, Barcelona, Spain
First published on 14th July 2015
A crystalline porous material showing one-dimensional (1-D) rectangular micropores (12 × 9 Å2) has been assembled from a semirigid macrocyclic tetraimine and EtOAc as the templating agent. The 1-D nature of the material is intrinsic to the conformationally rigid structure of a macrocyclic sub-unit bearing four cyclohexylidene residues. The multiple dispersive forces established among the aliphatic residues glutted the 1-D channels and provided thermal stability to the material at temperatures below 160 °C. Upon removal of the template, the structure of the empty solid exhibited permanent microporosity (SBET = 342 m2 g−1). Being a true molecular sponge, the channel framework of this material allowed the inclusion of a variety of molecular sample guests without compromising its crystalline nature. Remarkably, this crystalline material enabled the structure determination by X-ray diffraction of the included molecules. Theoretical studies demonstrated the vital role played by the dispersive forces in the overall stabilization of the crystal packing.
In this work, we present a new organic material featuring permanent porosity based on the macrocyclic tetraimine 3 (Scheme 1). A key aspect of our design is the use of multiple weak dispersive interactions as the glue element to organize and maintain together the components of the stable porous molecular material. Our report constitutes a paradigmatic example of the so-called Gulliver's principle10 applied to the design of porous organic materials.
Previously, we have described the synthesis of the 1,3-phenylene-bis-propargylic diamine 1 used in this work as rigid spacer to avoid the collapse of the macrocycle 3.11 The hierarchical self-assembly of the macrocycle 3 into a porous material arises from the combination of multiple van der Waals and aromatic interactions.12 These interactions are established between the phenyl and cyclohexylidene residues, the “sticky groups” of 3. The resulting porous material displays 1-D channels that are large enough for liquid–solid sorption applications.13
More specifically, we were interested in knowing if a crystalline porous material based only on dispersive forces, could be used as the crystalline solid support in X-ray diffraction studies by applying the “crystalline sponge method” developed by Fujita.14 In this method, accurate molecular structures of samples can be obtained by analyzing the diffraction patterns of single crystals without requiring the crystallization of the samples. However, the crystalline sponges reported so far are based on metal–organic frameworks (MOFs) that include electron-rich atoms, such as iodine, bromine, chlorine as well as transition metals in their structures.15 These atoms contribute heavily to the overall diffraction pattern masking the information available for the guest. From this point of view, our system based only in light-atoms, could offer an alternative to MOF-based crystalline sponges.
The single crystal X-ray structure showed the macrocycle 3 molecules stacked one on top of another forming a columnar assembly. The solid material displayed one-dimensional channels with distorted rectangular morphology along the b axis (Fig. 1). The total solvent accessible void volume within the channels calculated from PLATON/Squeeze, was 1096 Å3 (274 Å3 per each macrocyclic imine) representing 19.4% of the unit cell volume.17 We observed that the channels were filled with partially disordered EtOAc molecules (Fig. 1a). In the crystal packing, the four aromatic rings of 3 were tilted ca. 60° with regard to the plane of the macrocycle. Furthermore, the two phenyl rings located at the wider sides of 3 oscillated between two positions (75:25 occupancy) and acted as revolving doors (Fig. 1b).
Fig. 1 (a) Representative X-ray perspective side-view in b direction of 3@EtOAc crystals displaying disordered EtOAc molecules inside the channel. (b) Partial view showing the relative orientation of the 1,4-phenylene rings (see text). (c) Top view of a molecular cell unit of the monoclinic crystals (Vc = 5642 Å3, Z = 4) of 3@EtOAc (DC = 1.121 g cm−3; L ∼ 11.98 Å, W1 ∼ 9.13 Å; W2 ∼ 6.11–6.90 Å) (CCDC no. 1012389). Disorder and EtOAc solvent molecules are removed for clarity. The inset shows multiple van der Waals contacts between cyclohexylidene residues. |
From this structure, it became clear that the combined action of intermolecular displaced π–π aromatic stacking, as well as multiple CH (sp3)–CH (sp3) van der Waals interactions extending throughout the channel, provided the necessary stabilization of the assembly.
The results obtained in 13C CP-MAS NMR experiments confirmed the role of EtOAc as guest template. The carbon signals of the relatively mobile molecules of EtOAc included in the solid, but not those of the host tetraimine 3 that formed the porous framework were observed in the standard mode 13C MAS spectra (Fig. S7 in ESI†).
Fig. 2 X-ray perspective top-view of the apohost 3(I) (Vc = 5639 Å3, Z = 4) (DC = 0.9859 g cm−3; L ∼ 11.98 Å, W1 ∼ 9.20 Å; W2 ∼ 6.90 Å) (CCDC no. 1012398). |
More important, the crystal packing of 3(I) was almost identical to that of the 3@EtOAc solvate (see Table S1 in ESI†), indicating that the porous structure remained intact after vacuum treatment.
The relatively easy removal of the included EtOAc molecules at room temperature implied that they were weakly retained in the channels.
Remarkably, the 3@EtOAc solvate could be regenerated by simply soaking crystals of the empty material 3(I) in EtOAc. These findings were consistent with the absence of strong dipole–dipole or hydrogen bonding interactions between the cycloimine 3 and the included EtOAc molecules.
The dynamic nature of the guest exchange constituted an essential feature of this material. The kinetics of the guest exchange process was monitored by using 13C NMR. Crystals of 3@13C-[1.2]-EtOAc obtained by treating 3(I) with 13C labelled EtOAc were suspended in unlabelled EtOAc.
We observed that the combined integral values of the isotopic 13C doublet grew steadily with time and reached a plateau after 8–10 h (Fig. 3). This experiment demonstrated that the almost complete exchange of the included 13C labeled EtOAc was slow in the human time-scale. Desolvation of EtOAc was also monitored using thermogravimetry. The TGA trace of 3@EtOAc revealed an 11.8% weight loss in a temperature range from room temperature to 175 °C. This weight loss corresponded to the vaporization of 1.3 molecules of EtOAc per macrocycle, therefore indicating the 3:4 stoichiometry found for this complex (see Fig. S9 in ESI†). In complete agreement, the differential scanning calorimetry (DSC) thermogram of 3@EtOAc exhibited an endothermic peak at around 160 °C due to vaporization of the entrapped guest solvent (Fig. 4b).
On the other hand, the porosity of 3(I) was also confirmed by gas absorption of nitrogen at 77 K (Fig. 4a). The empty solid 3(I) exhibited a typical type I adsorption isotherm indicative of a microporous (<2 nm) structure with significant specific surface area (SBET = 342 m2 g−1).18 In addition, the DSC trace of the apohost 3(I) displayed an exothermic and irreversible process at around 175 °C (Fig. 4b and S10†), which we assigned to a monotropic solid-phase transition of 3(I). In both cases, the X-ray structure analysis of the material obtained after thermal treatment above 160 °C, revealed the formation of a new non-porous polymorphic form 3(II).
The X-ray structure of 3(II) (see Fig. S11 in ESI†) showed densely packed layers of displaced macrocyclic imine units 3 (DC = 1.127 g cm−3) and the concomitant loss of the columnar assembly characteristic of the 3(I) polymorphic form.
To test this functionality, we soaked single crystals of 3@EtOAc and 3(I) in different neat liquid samples with molecular volumes ranging from 55–60 Å3i.e. nitromethane and ethylene glycol, to >200 Å3 of Z-stilbene (see, Table S2 in ESI†). In favourable cases, i.e. less than 5 min for nitromethane, the initially floating crystals of 3@EtOAc gradually sank to the bottom of the vial providing a visual indication of the exchange status. During the process, we did not observe any apparent degradation of the crystal integrity. The exchange of the guest samples was confirmed by 1H NMR of CDCl3 solutions of the solvates. Also, 13C CP-MAS NMR experiments confirmed the inclusion of the liquid guests within the porous framework. The carbon signals of non-volatile guests (i.e. diethyl squarate) appeared slightly upfield shifted (0.5–2 ppm) with respect to those of diethyl squarate in the bulk (Fig. S8 in ESI†).
Owing to the poor solubility of the solvates in CDCl3, the molar ratios were better evaluated by digesting the solvates with DCl/D2O (50 μL) before registering the 1H NMR spectra in a mixture of DMSO-d6:CDCl3 (2.5:1 v/v, 600 μL) (see, ESI†).15c In general, the molar ratios of 3:sample evaluated by 1H NMR after a common period of 72 h, agree well with the crystallographic stoichiometries (Table S2 in ESI†). In one case, the average experimental molar ratio calculated for Z-stilbene of 0.8 was consistently inferior to that expected for a 1:1 3@Z-stilbene solvate. This result is likely due to the incomplete inclusion of Z-stilbene within the channels of the 3(I). Therefore, the molecular volume of Z-stilbene of 216 Å3, approaching the maximum available void volume of 274 Å3 per macrocycle calculated using PLATON/Squeeze,17 suggested us a practical upper limit for the effective inclusion of samples in 3(I).
Fig. 5a–h, depict ORTEP views of representative X-ray structures obtained by guest exchange from single crystals of 3@EtOAc or alternatively, by soaking a single crystal of the empty apohost 3(I) with ca. 1 mL of the liquid sample. Both methods produced identical results in all tested cases (see ESI for X-ray structural data†). Both, the empty material 3(I) and the solvate 3@EtOAc allowed the formation of inclusion solid complexes (solvates) with aromatic, unsaturated and aliphatic neutral substrates bearing a variety of functional groups such as hydroxy, aldehyde, ester and nitro groups, among others.
Fig. 5 (a to h) ORTEP views of different liquid molecules diffused into the channel structure of 3(I) with thermal ellipsoids set at 50% probability. Disorder is removed for clarity. CCDC codes 1012389 to 1012397 and 1063714 contain the structural X-ray data (see ESI for RX structural details†). |
Since the included guests could, eventually, establish hydrogen bonds or dipole–dipole interactions with 3 that would lead to degradation of the porous structure of the solids, the obtained results in the guest exchange/inclusion experiments confirmed that the interactions established between 3 and the incoming sample guests must be feeble.
Remarkably, the inclusion of the different guests into the crystal lattice did not disrupt the pore size significantly (ΔL < 2%; ΔW1 < 2.6%). However, to obtain a better fit between the size and shape of the pore and those of the included molecules (induce fit), the two opposite phenyl rings located on the wider sides of the macrocycle 3 oscillated with a mechanism that resembles the functioning of the revolving doors. The adaptation process experienced by the pore was revealed by careful analysis of the X-ray structures of 3@EtOAc, 3@p-anisaldehyde, and 3(I)@S-(−)-nicotine (ΔW2 ∼ 20%).
The induced fit process, in which the width of the channels responded to the size and shape of the included sample guests took place without any loss of crystal integrity (Table S1 in ESI†). In contrast, the X-ray structure of the non-porous polymorphic form 3(II) displayed a narrowed entrance of the pore that can be referred as “closed” doors. Most likely, this structural change was responsible for inhibiting the guest sorption of liquid guests (W1 ∼ 9.91 Å; W2 ∼ 5.54 Å).
In addition, the interaction energy computed for the dimer using the crystallographic coordinates (without optimization) gave a very similar value (−32.9 kcal mol−1). Removal of the cyclohexylidene residues led to a significant reduction of the interaction energy (−14.5 kcal mol−1). This finding supported the strong influence exerted by the C–H⋯H–C interactions to the stability of the porous material.19
The energy difference calculated for the two energy-minimized dimers (18.4 kcal mol−1) represented an estimation of the contribution of these interactions to the overall dimerization process. Since four cyclohexylidene residues are present in the macrocycle 3, each intermolecular cyclohexylidene⋯cyclohexylidene interaction contributes in 4.6 kcal mol−1. This value was in good agreement with previously reported high level ab initio results for dodecahedrane.19
To further analyze the non-covalent interactions in the dimer 32 we used the NCI (Non Covalent Interaction) index.20 This visualization index enables the identification and visualization of weak non-covalent interactions efficiently.
The isosurfaces corresponded to both favourable and unfavourable interactions, as differentiated by the sign of the second density Hessian eigenvalue and defined by the isosurface color. NCI analysis allowed an assessment of host–guest complementarity and the extent to which weak interactions stabilize a complex. Fig. 7, shows the representation of the NCI plot computed for the dimer 32. In addition to a panoply of C–H⋯H–C interactions, two weak equivalent and oppositely located parallel-displaced π–π interactions and four C H/π interactions involving the π-system of the triple bond were also observed. The largest isosurface was located around the cyclohexylidene rings that established a multitude of C–H⋯H–C interactions with another cyclohexylidene and the aromatic ring.
The process was favourable by 4.1 kcal mol−1 (Fig. 8). Subsequently, we assessed the energy required to exchange the expelled EtOAc by a new SQA molecule. This second process was also favourable by 3 kcal mol−1. This sequence was repeated to evaluate the energy of the exchange processes of the second and third molecules of EtOAc by SQA molecules (see Fig. S15–S17 in ESI†).
Fig. 8 Favourable energy change due to the substitution of one EtOAc molecule by one diethyl squarate molecule in a minimalist representative channel of 3@EtOAc. |
These processes were also energetically favourable and showed similar magnitudes for the exchange energies. Hence, the results of this theoretical study agreed well with the aforementioned experimental results and supported the viability of the liquid–solid guest exchange process.
Footnotes |
† Electronic supplementary information (ESI) available: Synthesis of macrocycle 3 together with characterization data. TGA and DSC traces. CIF files of all structures. Computational details and additional calculations. CCDC 1012389–1012397 and 1063714. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c5sc01838b |
‡ Crystal data for 3@EtOAc (CCDC no. 1012389): C65.36H70.72N4O2.68, M = 955.18, monoclinic, a = 34.794(16) Å, b = 5.836(3) Å, c = 28.131(14) Å, α = 90.00°, β = 97.862(11)°, γ = 90.00°, V = 5659(5) Å3, T = 100(2) K, space group C2/c, Z = 4, 35400 reflections measured, 10696 independent reflections (Rint = 0.0434). The final R1 values were 0.0597 (I > 2σ(I)). The final wR(F2) values were 0.1624 (I > 2σ(I)). The final R1 values were 0.0931 (all data). The final wR(F2) values were 0.1931 (all data). Crystal data for 3(I), (CCDC no. 1012398): C60H60N4, M = 837.12, monoclinic, a = 34.8696(19) Å, b = 5.8099(3) Å, c = 28.1212(16) Å, α = 90°, β = 98.143(2)°, γ = 90°, V = 5639.6(5) Å3, T = 100(2) K, space group C2/c, Z = 4, 26236 reflections measured, 8377 independent reflections (Rint = 0.0577). The final R1 values were 0.0602 (I > 2σ(I)). The final wR(F2) values were 0.1686 (I > 2σ(I)). The final R1 values were 0.0825 (all data). The final wR(F2) values were 0.1883 (all data). Crystal data for 3(II), (CCDC no. 1012399): C30H30N2, M = 418.56, triclinic, a = 6.5695(4) Å, b = 13.4759(8) Å, c = 14.1239(9) Å, α = 91.0310(10)°, β = 99.2650(10)°, γ = 91.6850(10)°, V = 1233.21(13) Å3, T = 100(2) K, space group P, Z = 2, 13275 reflections measured, 13275 independent reflections (Rint = 0.0000). The final R1 values were 0.0529 (I > 2σ(I)). The final wR(F2) values were 0.1334 (I > 2σ(I)). The final R1 values were 0.0860 (all data). The final wR(F2) values were 0.1550 (all data). |
This journal is © The Royal Society of Chemistry 2015 |