Witold M.
Bloch
*a,
Ravichandar
Babarao
bc and
Matthew L.
Schneider
a
aDepartment of Chemistry, The University of Adelaide, Adelaide, Australia. E-mail: witold.bloch@adelaide.edu.au; Tel: +61 8 8313 5039
bSchool of Science, RMIT University, Melbourne, Victoria 3001, Australia
cCSIRO, Normanby Road, Clayton 3168, Victoria, Australia
First published on 5th March 2020
Synthetic porous materials composed of metal–organic polyhedra (MOPs) have found application in topical areas such as gas storage, separation and catalysis. Control over their physical properties (e.g. porosity) has typically been achieved through ligand design or judicious choice of metal ions. Here, we demonstrate pore-size control and on/off porosity in Cu4L4 MOPs by exploiting their structural non-rigidity. We report an aldehyde-functionalised MOP (1) that can be isolated in five distinct solvatomorphs, each exhibiting different structural flexibility. When soaked in MeOH, two of these solvatomorphs undergo a rapid transformation to a thermodynamically favoured phase, whilst in acetone they template the crystallisation of an entirely new crystal packing. We support these findings by single and powder X-ray diffraction and rationalise the observed phase transformations by lattice energy calculations. Of the five solvatomorphs, three can be obtained as solvent-exchanged pseudo-polymorphs with distinct porosities in their activated form (SABET = 35–455 m2 g−1). Further control over the crystal packing of MOPs is achieved through covalent post-assembly modifications, which promote the crystallisation of isoreticular 2-D sheet-like structures.
Metal–organic polyhedra (MOPs) are a class of discrete nanocages commonly formed by the self-assembly of dicarboxylate linkers and bimetallic paddle-wheels nodes (e.g. Cu2, Rh2 and Cr2).2,12,15 Solvent molecules typically occupy the apical coordination sites of the paddle-wheel nodes, and thus, desolvation allows access to unsaturated metal centres. This feature has motivated extensive studies examining MOPs for selective gas adsorption16 and heterogeneous catalysis.14,17,18 However, the crystal engineering of MOP solids is an area dominated more by serendipity than design; the weak intermolecular interactions that direct the crystal packing of MOPs can depend on the nature of the paddlewheel-bound solvent, making it possible to crystallise different solvatomorphs of the same cage compound.19 To add to the complexity, solvent exchange processes employed en-route to activation can lead to unknown crystalline or amorphous phases,12 making it difficult to control the pore-size and porosity of MOPs. Therefore, gaining control over the crystal packing of MOPs, and hence their physicochemical properties, remains a challenge.
Much of the efforts aiming to design MOP-based solids have focused on immobilising MOP particles in mesoporous supports,20,21 or on approaches that increase intermolecular interactions in the solid-state. For example, introducing ligand-based functionalities such as alkyl substituents,22 bulky hydrophobic groups,23 or coordinative O-donors24 has led to more robust and structurally rigid crystal packings. Owing to the solubility of MOPs, covalent post-assembly modifications (PAM)25 have also been carried out to modify the interior26 or exterior surface functionality.27 However, the development of covalent PAM on paddlewheel-based MOPs is still an emerging field of study, with only a few examples reported.
Herein we show that the crystal-packing and porosity of MOPs can be controlled in the same way as their purely organic analogues; by taking advantage of their rich phase behaviour and amenability to PAM. Thus, we report the synthesis and characterisation of an aldehyde-functionalised Cu4L4 MOP (1, Fig. 1a). By varying the crystallisation conditions, 1 can be obtained in five distinct solvatomorphs with different structural flexibilities. This is demonstrated by phase-dependant transformations which we characterise by X-ray diffraction and rationalise by lattice energy calculations. We show that the activation and porosity of 1 is phase dependant, and non-porous forms can be switched to porous forms by solvent treatment (Fig. 1b). Owing to the solubility of 1, PAM carried out at the cage exterior facilitate control over the crystal packing by favouring the formation of 2-D sheet-like motifs (Fig. 1c).
Fig. 2 (a) A scheme showing the self-assembly of 1 in solution from LH2 and Cu(OAc)2 (S = solvent); (b) simulated and experimental PXRD patterns of phase 1a; (c) 1H NMR spectrum (500 MHz/DMF-d7, 25 °C): (i) LH2, (ii) LH2 + Cu(OAc)2 (* indicates solvent peak).‡ Above (b) and (c) are photos of crystals of 1a (left) and a DMF solution of 1 (right). |
Self-assembly of LH2 with 1.1 equivalents of Cu(OAc)2·2H2O in DMF resulted in the formation of a soluble Cu4L4 lantern cage (1) as indicated by 1H NMR spectroscopy and UV-Vis spectroscopy (Fig. S6†).26,28 Owing to the anomalous magnetic moment of the Cu2 paddlewheel,29 proton resonances not immediately adjacent to the metal nodes could be resolved in the 1H NMR spectrum of 1.‡ For example, the aldehyde resonance of 1 appears upfield shifted relative to the free ligand (ca. −0.2 ppm, Fig. 2c). Slow-vapour diffusion of MeOH into a DMF solution of 1 over a period of 5 days yielded rhombohedral crystals of 1a in 87% yield (based on CHN analysis). Single-crystal X-ray analysis (SCXRD) revealed that 1a crystallises in the triclinic space group P with half of two crystallographically independent cages in the asymmetric unit. Closer inspection revealed that the ligand backbones in one cage molecule adopt a flattened conformation when viewing the cage down the axis of the Cu2 paddlewheels. These nodes are coordinated by four carboxylate ligands and two axial DMF ligands. The second crystallographically unique cage adopts a pseudo D4h-symmetric conformation where the ligand backbones arrange in a mutually perpendicular fashion.19 Here, DMF and H2O ligands coordinate to the exterior and interior sites of the paddlewheels, respectively. The two chemically unique cage molecules present in the crystal packing of 1a give an overall formula of [Cu4L4(DMF)4]·[Cu4L4(DMF)2(H2O)2] (Fig. 3a and S33†). The extended packing of 1a is stabilised by multiple intermolecular π–π interactions that exist between the benzaldehyde cores of 1 (closest contact = 3.33 Å). This motif produces ordered channels along the a axis (Fig. 3a) with aperture dimensions measuring 3.8 × 6.1 Å (when omitting coordinating solvents).
Fig. 3 (a) Single crystal X-ray structures of solvatomorphs 1a–e in a scheme portraying their structural transformations. The crystal packing of 1a–e is represented as space-filling (coordinating solvents removed for clarity) and alongside are depictions of the orientations of cage molecules within each phase. Each phase can be regenerated through recrystallisation (*co-solvent = H2O); (b) PXRD patterns highlighting the solid-state transformations of 1c and 1d → 1b upon exposure to MeOH. Full details of the transformations shown in (a) can be found in the ESI (Fig. S11–S16†). |
Needle-shaped crystals of 1b were obtained when a mixture of LH2 and Cu(OAc)2·2H2O was heated in DMSO at 85 °C for 16 h. The same crystal packing could also be obtained by SVD of MeOH into a DMA solution of 1 (1b-DMA, Fig. S13c†), although only the former conditions produced crystals suitable for SCXRD. 1b-DMSO crystallises in the monoclinic space group C2/c with half of the lantern cage in the asymmetric unit. Both interior and exterior paddlewheel sites are occupied by disordered DMSO ligands to give an overall formula of [Cu4L4(DMSO)4]. The crystal packing of 1b is stabilised by multiple intermolecular π–π interactions between the phenyl carboxylate and benzaldehyde moieties of 1 (shortest distance = 3.32 Å). This bilateral π-stacking motif results in molecules of 1 arranging in an edge-to-edge fashion. When viewing the extended structure along hkl = 111 plane, open channels intrinsic to the cage structure are evident and measure approximately 4.4 × 4.6 Å in diameter (Fig. 3a).
Slow diffusion of atmospheric moisture into a DMF solution of 1 resulted in the overnight formation of large block-shaped crystals of a new phase (1c), as indicated by SCXRD and PXRD analysis. We postulate that under these conditions, 1c forms in favour of 1a due to the rapid rate of crystallisation (24 h) compared to DMF/MeOH (5 days).91c crystallises in the monoclinic space group I2/m, with a quarter of two crystallographically independent cages in the asymmetric unit. In contrast to 1a, both cage molecules have identical paddlewheel compositions; DMF and H2O ligands occupy the exterior and interior apical sites, respectively. The overall formula of 1c is therefore [Cu4L4(DMF)2(H2O)2] (Fig. 3a and S33†). Interestingly, the extended structure of 1c comprises of a mutually perpendicular face-to-face cage packing. This arrangement creates open channels along the a and c axis, with the largest pore dimensions measuring 5.4 × 6.0 Å. However, this packing motif lacks defined π–π interactions and is stabilised by weak hydrogen bonding between the coordinated DMF ligands and the aldehyde group of 1 (C–H⋯O = 2.31 Å, 119°).
Another distinct crystalline phase could be obtained by slow diffusion of atmospheric moisture into a DMA solution of 1. 1d crystallises in the monoclinic space group P21/c, with half of the lantern cage in the asymmetric unit. In this solvatomorph, both apical coordination sites of the Cu2 paddlewheel are occupied by DMA ligands (Fig. S33†). In the crystal packing of 1d, cage molecules arrange in an edge-to-edge fashion and exhibit intermolecular π–π interactions between the benzaldehyde and phenyl carboxylate moieties. Despite possessing a similar crystal density to that of 1b (Table S2†) the arrangement of cage molecules in 1d does not result in defined solvent-accessible channels (Fig. 3a).
Owing to the solubility of 1 in DMF and DMA, crystals of each phase could be recrystallised to recycle between 1a–d (Fig. 3a).
The crystallisation of 1a–d indicated to us that the solvent and precipitant both contribute in determining the extended packing of the MOP structure. As such, we posited that in the four structures of 1, the apical paddlewheel solvent stabilises the crystal packing to different degrees. In order to examine this hypothesis, we subjected solvatomorphs 1a–d to solvent exchange with MeOH and acetone, as these solvents are commonly used en-route to activation.32,33 With the X-ray structures of 1a–d in hand, we were able to use powder diffraction to map out the solvent-induced transformations in bulk samples of 1 (Fig. 3b and S11–S16†).
PXRD analysis revealed that both 1a and 1b retain their overall crystal packing after solvent exchange with MeOH (Fig. 3b, S11 and S13†). Whilst the methanol exchanged form of 1b was stable indefinitely, a time-dependant PXRD experiment of 1a revealed that after 9 h, a slow transformation to 1b begins to occur and is complete within ∼24 h (Fig. 3a and S12†). In contrast, washing 1c and 1d in MeOH resulted in a rapid (<5 minutes) and quantitative transformation to 1b as indicated by PXRD (Fig. 3b, S14 and S15†). Notably, these conversions proceed entirely in the solid-state as 1 is insoluble in MeOH and DMF/MeOH mixtures (Fig. S8†).§ Although the mechanism of these types of transformations is difficult to elucidate,8 we speculate that the relative cage re-orientations that occur during this process are different for each phase (Fig. S34†).
Interestingly, we observed a different transformation when solvatomorphs of 1 were immersed in acetone. Whilst 1a and 1b retained their crystal packing even after 3 months (Fig. S11 and S13†), 1c and 1d both transformed to a new phase (1e) over a period of 3 weeks (Fig. S14 and S15†). The crystals of the former sample were of sufficient quality for structural determination by synchrotron SCXRD. 1e crystallises in the monoclinic space group C2/m and contains a quarter of the cage in the asymmetric unit. Here, the apical sites of the Cu2 paddlewheel nodes are occupied by exterior acetone (DMK) and interior water ligands to give the formula [Cu4L4(DMK)2(H2O)2]. The lantern cage adopts an extremely flattened conformation that facilitates close intermolecular π-stacking between the benzaldehyde and phenyl carboxylate moieties (shortest distance = 3.29 Å). The eclipsed face-to-face cage packing results in open channels along the c axis with apertures of 6.0 × 8.9 Å (Fig. 3a). Unsurprisingly, we found that the PXRD pattern of 1e remains unchanged upon solvent exchange with MeOH (Fig. S16†).
To further investigate the formation of 1e, we soaked the as-synthesised crystals of 1c and 1d in acetone-d6. Indeed, we observed trace quantities of 1 in the 1H NMR spectrum of the suspension, albeit only in the presence of DMA or DMF; washing the samples further with acetone-d6 (in order to exchange out the remaining solvent) led to the disappearance of resonances originating from 1 (Fig. S8†). As the crystals of 1c and 1d do not completely dissolve during the transformation, we suggest that the conversion to 1e is solvent assisted,34,35 although a traditional recrystallisation cannot be ruled out. It is noteworthy that 1e could not be obtained through conventional crystallisation methods or by the conversion of other solvatomorphs (Fig. S17, Table S1†), suggesting that 1c and 1d act as crystalline templates for this transformation.
To rationalise the observed transformations we performed lattice energy calculations on the set of five structures based on periodic density functional theory (DFT) calculations as implemented in the software package VASP 5.4.4.36 We note that lattice energy comparisons are commonly made for purely organic polymorphs (where solvent does not play a structural role).37,38 Nevertheless, we expected that the relative lattice energies of the desolvated forms of 1a–e should give insight into their transformations despite their experimental accessibility. In our DFT calculations, electron exchange and correlation were described using the generalised gradient approximation Perdew, Burke, and Ernzerhof (PBE) form39 and the projector-augmented wave potentials were used to treat core and valence electrons.40 To account for the long-range dispersion interactions, we included the dispersion corrections using DFT-D3 method.41 A full optimisation of the structures was performed allowing both the atomic and cell parameters to relax in the absence of the apical paddlewheel solvents. We observed no substantial change in the structural conformation and the lattice parameters for each of the five structures in the DFT optimisation (full details are given in the ESI†). The results of the DFT calculations were consistent with the expected correlation of higher density structures (e.g.1a, 1b and 1e) being lower in energy compared to structures with lower density (e.g.1c and 1d, Table S6†).37 As observed experimentally, 1b is indeed the lowest energy structure within the set of five phases (Table S4†), closely followed by 1e (+15.3 kJ mol−1, Fig. 4). The fact that these two phases are energetically similar yet do not interconvert through solvent exchange suggests that their crystal packing is stabilised predominantly by intermolecular cage interactions rather than the apical paddlewheel solvent. Relative to 1b, the lattice energy of 1a differs by only +22.3 kJ mol−1 which is in agreement with its metastability and slow transformation to 1b in MeOH. On the other hand, the lattices energies of 1c and 1d are comparatively higher (+106 kJ mol−1 and +69 kJ mol−1, respectively) which also supports their transformation to the low energy structures 1b (in MeOH) or 1e (in acetone). Therefore, compared to the other phases, 1c and 1d are largely stabilised by intermolecular interactions involving the coordinated paddlewheel and pore-bound solvents.
The structure of the activated phases was investigated by deriving pore-size distributions (PSD) from the low-pressure region of the isotherm data. PSD analysis revealed that indeed, 1a, 1b and 1e display distinct pore-sizes in the microporous region, despite possessing the same chemical structure. 1a shows a broad pore size centered around 10 Å, whilst 1b and 1e display much sharper distributions. 1b has pore volume maxima at 8.0 Å and 11.8 Å, whilst 1e possesses a sharp maximum at 5.9 Å (Fig. 5c).
Due to the weak, non-covalent interactions that stabilise their crystal packings, the PXRD patterns of the desolvated and activated samples showed shifting and broadening of the low-angle peaks with a reduction of long-range order compared to their parent MeOH-exchanged phases (Fig. 5b and S21†). This is common for carboxylate paddlewheel MOPs and general strategies to totally prevent lattice changes upon desolvation remain scarce.42,44 Nevertheless, the low angle diffraction peaks of these samples suggest that the activated forms of 1a, 1b and 1e possess unique crystal packings, which is also supported by their pore-size distributions. Therefore, we infer that the structural changes accompanying activation are different for each pseudo-polymorph,32 especially since different adsorbing forms of 1 are the end result. Upon re-solvating the activated forms of 1b and 1e in MeOH, we observed the appearance of sharp PXRD peaks consistent with 1b-MeOH. Resolvation of the activated samples of 1a, 1b, and 1e in their mother liquor, however, resulted in the regeneration of their crystalline, as-synthesised forms (Fig S20†). Thus, the resolvation experiments indicate that the activated pseudo-polymorphs retain structural flexibility and further support 1b as the lowest energy structure in the group of five solvatomorphs.
Indeed, the deliberate instalment of the exterior aldehyde groups in 1 enabled us to carry out Lewis-acid catalysed imine condensation. Heating a DMF-d7 solution of 1 at 45 °C in the presence of 0.1 equivalents of Sc(OTf)3 (ref. 45) and 16 equivalents of o-toluidine resulted in the formation of 2 after 16 h. The 1H NMR spectrum of the reaction mixture revealed a characteristic imine resonance at 8.6 ppm and an absence of the aldehyde resonance normally observed at 9.9 ppm (Fig. 6b). Subjecting the reaction mixture to slow-vapour diffusion with diethyl ether yielded rod-like crystals of 2 after three days. SCXRD analysis confirmed that the complete conversion of all four aldehyde substituents had occurred without interfering with the Cu4L4 connectivity (Fig. 6c). 2 crystallises in the triclinic space group P with half of the formula in the asymmetric unit. In the extended structure, the newly formed cage packs in an interdigitated edge-to-edge fashion. This arrangement is stabilised by numerous intermolecular π–π interactions and results in an overall 2D sheet-like packing (Fig. 6d). Within this motif, the average Cu⋯Cu distance between cage units measures 18.7 Å, which is an increase of 3.7 Å compared to the furthest Cu⋯Cu distance observed in the crystal packing of 1 (1c = 15.0 Å).
Fig. 6 (a) Scheme illustrating the covalent PAM of 1 to give 2 or 3; (b) 1H NMR (500 MHz/DMF-d7) of (ii) the reaction mixture of 2 after heating and (i) before heating.‡ Excess aniline shown in grey, *DMF solvent peak; (c) X-ray structure of 2; (d) perspective view of the crystal packing of 2; (e) perspective view of the crystal packing of 3. |
We also performed covalent PAM of 1 with 9-ethyl-9H-carbazol-3-amine due to its solubility and extended aromatic structure. This allowed us to thus further examine extensions to the cage architecture and their subsequent effect on the crystal packing. Once again, 1H NMR spectroscopy revealed a >98% conversion of 1 to 3 under the similar reaction conditions to that of 2 (Fig. S9a†). Slow vapour diffusion of diethyl ether into the reaction mixture yielded X-ray quality crystals of a monoclinic phase with the space group P21/n. Interestingly, in the extended structure of 3, molecules of interdigitate in an identical manner to that of 2. The bulkier carbazole moiety pushes the cages further apart within the 2D layers to give an average inter-cage Cu⋯Cu distance of 21.1 Å (Fig. 6e).
We confirmed the bulk purity of 2 and 3 by IR spectroscopy and PXRD (Fig. S9b and S18†). The former revealed new CN stretches at 1619 cm−1 and 1630 cm−1 for 2 and 3 respectively, as well as the absence of the aldehyde stretch (1698 cm−1). Attempting other crystallisation conditions to obtain solvatomorphs of 2 and 3 led to the same crystal packing or the isolation of poorly crystalline precipitates that were not further investigated. By providing more opportunities for π–π interactions, it is evident that the range of solvatomorphs may be limited compared to 1. However, this may be advantageous for the crystal engineering of MOP solids; other functionalities (e.g. H-bonding) could be introduced through covalent PAM to control the structural rigidity, thus limiting the effect that coordinating solvent plays in determining the crystal packing.
Finally, we demonstrated that 1 can be covalently modified in solution via imine condensation. In the solid-state, the covalent extensions of 1 led to an interdigitated 2D sheet-like motif, which was observed in both of the modified samples. We think that this approach may open up new avenues for the design and crystal engineering of porous MOP materials. Indeed, the assembly of hierarchical materials is an emerging field of research,46–48 and the exterior modifications of 1 show promise for assembling COF-like materials from MOP building blocks. Efforts towards this goal are currently underway in our laboratory.
Footnotes |
† Electronic supplementary information (ESI) available: Full experimental details and supporting analysis. CCDC 1969619–1969625. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/d0sc00070a |
‡ Owing to the paramagnetic Cu2 paddlewheel node of 1–3, proton resonances in close proximity to the metal node are significantly broadened.26 |
§ During the transformation of 1c/d to 1b, the crystals of the parent phase fragment rapidly. Despite the insolubility of 1 in MeOH and even mixtures of DMF/MeOH, we cannot completely rule out the possibility of an extremely rapid recrystallisation. |
This journal is © The Royal Society of Chemistry 2020 |