Emily G.
Meekel
and
Andrew L.
Goodwin
*
Department of Chemistry, Inorganic Chemistry Laboratory, University of Oxford, South Parks Road, Oxford OX1 3QR, UK. E-mail: andrew.goodwin@chem.ox.ac.uk; Tel: +44 (0)1865 272137
First published on 1st March 2021
We survey the importance of correlated—as opposed to random—disorder in metal–organic frameworks (MOFs), focussing on correlations in vacancy arrangements, composition, and linker orientations in some canonical systems. We explore the link to physical and chemical properties, and highlight the various experimental and computational tools available for characterising and interpreting correlated disorder in MOFs.
Disorder is sometimes random; more often it is correlated.12 A vacancy at one linker site in an otherwise crystalline MOF, for example, may affect the likelihood for neighbouring linkers also to be vacant. Such vacancy correlations can propagate throughout a structure, collectively determining a variety of bulk properties—including porosity.2,13 Likewise, mixed-metal and/or mixed-linker MOFs can have non-statistical distributions that determine their corresponding chemistry; such is the basis, after all, of the ‘multivariate’ MOF approach which exploits compositional heterogeneity to improve functionality.14 And in orientationally-disordered MOFs, correlations in linker orientations can also affect mechanical and host/guest responses.15 These various distinctions are illustrated schematically in Fig. 1. It is no accident then that controlling the correlations within disordered MOFs—whether involving vacancy, composition, or linker orientation distributions—is an essential aspect of defect engineering.2,16–18
Distinguishing correlated disorder from random disorder is often a nontrivial experimental and computational challenge. Conventional crystallography, which has played such an important historical role in MOF science, probes only the average structure of materials: it is demonstrably insensitive to correlations within disordered states.19 Yet a number of developments have now converged to equip us, as a community, with the tools needed to probe correlated disorder in MOFs. Experimentally, recent advances in solid-state NMR, X-ray diffuse scattering, pair distribution function (PDF), and electron microscopy measurements are providing valuable empirical insight into MOF defect structures.13,20–22 At the same time, ab initio density functional theory (DFT) calculations of a sufficient scale to interrogate non-crystalline MOF architectures are becoming increasingly tractable.6,23 The development of efficient high-quality MOF force-field methodologies extend these approaches beyond the nanometre scale.24,25
With these recent developments in mind, we survey in this Highlight article the phenomenology of correlated disorder in MOFs. Our focus is on disorder of three primary types—correlated vacancies, compositional disorder, and orientational disorder—and we place our discussion in the context of the behaviour of a small selection of canonical systems. Our emphasis is also threefold: (i) the nature of correlated disorder in these MOFs; (ii) the use of experimental and computational tools to characterise correlations; and (iii) the implications for physical and chemical properties. Our article concludes with a forward-looking discussion of the future challenges and opportunities in the field.
Fig. 2 Correlated vacancies in MOFs. (a) Nanodomain model for UiO-66(Hf) developed in ref. 27: cluster vacancies are correlated to give domains of reo topology embedded within the parent fcu framework. Reproduced from ref. 27 with permission from Springer Nature. (b) High-resolution transmission electron microscopy images of missing-cluster/missing-linker nanodomains in UiO-66.33 The scale bar measures 5 nm. The lower panels denote (left–right) atomistic model, calculated potential map, and observed potential map. Adapted from ref. 33 with permission from Springer Nature. (c) Comparison of sql–kgm transformation mechanisms (left–right) in non-defective AlTz-53-DMF (top) and defective AlTz-53-DEF (bottom); the coordinating strength of included solvent plays a key role in facilitating network rearrangement.34 The green nodes denote aluminium hydroxide/carboxylate columns (illustrated on left). Adapted from ref. 34 with permission from Springer Nature. (d) The single-crystal X-ray diffuse scattering pattern of Zn3[Co(CN)6]2 PBA is characteristic of non-random correlations in the underlying pore network structure (green tubes) that permeates the cubic transition-metal cyanide framework. As disordered as the pore network is, it is not random.13 Adapted from ref. 13 with permission from Springer Nature. |
The driving force for correlated disorder in the UiO-66 family is thought to be the mechanical stability of reo-type connectivity within the fcu net.35 Other, lower symmetry, vacancy distributions are mechanically unstable with respect to shear and so do not survive during crystal growth; this biases the system away from purely statistical vacancy arrangements. An additional consideration is the structure-directing interaction between the UiO-66 framework and both counterions and solvent in the reaction mixture.36 The picture that emerges is one in which vacancy incorporation in UiO-66 is a kinetic effect—and hence why defect density can be controlled by e.g. synthesis temperature,37 solvent,38 or modulator concentration or pKa.16
From a functional perspective, there are a number of implications of random vs. correlated vacancy defects is this family. Perhaps the most obvious is the effect on sorption behaviour, because this depends on the nature of the pore network generated by connected vacancies.16 Related is the dramatic increase in proton conductivity observed for UiO-66 samples with high linker-vacancy concentrations—although the specific link between missing-linker arrangements and efficiency of proton transport in this system remains to be explored in detail.39 Correlated cluster absences affect mechanical properties, such as thermal expansivity40 and compressibility.10 And DFT calculations show the distribution of both linker and cluster vacancies control the spatial localisation of frontier orbitals, in turn affecting both electronic and catalytic behaviour.36,41
The importance of vacancy correlations in other MOF systems is also increasingly apparent. In the aluminium carboxylate MOF AlTz-53-DEF, for example, correlated linker vacancies facilitate a remarkable reversible switching between sql and kgm network topologies during desolvation/resolvation cycles [Fig. 2(c)].34 Even for the well-known MOF-5, adsorption isotherms are sensitive to the distribution of linker vacancies—a result demonstrated using grand canonical Monte Carlo (GCMC) simulations.42 A similar link between linker vacancy arrangements and water adsorption behaviour has been made particularly clear in the case of MOF-801—a fumarate-containing relative of UiO-66.43 And in DUT-49, a MOF famous for its negative gas adsorption (NGA) behaviour,44 cluster vacancies in its fcu network quench the NGA effect, with vacancy correlations thought to play a key role.45
Dense coordination networks can harbour disordered correlated vacancy distributions, the form and origin of which we expect may prove relevant in due course to the behaviour of MOFs. For example, Prussian blue analogues (PBAs) are heavily defective cyanide frameworks with the pcu topology.46,47 Single-crystal X-ray diffraction measurements of PBAs reveal highly structured diffuse scattering,13,48,49 which has recently been interpreted using the three-dimensional difference pair distribution function (3D-ΔPDF) approach [Fig. 2(d)].13 The scattering arises from correlated hexacyanometallate vacancies, which avoid occupying neighbouring sites. A key point is that the precise correlations—which control the sorption and transport characteristics of the resulting pore network—can be tuned by variation in PBA composition and synthesis route.13 Conceptually similar vacancy patterns are found in the thiocyanate analogues, albeit with long-range order.50 And in formate perovskites, cation vacancies again interact by avoiding neighbouring sites; this leads to a disorder/order transition in [GUA] Mn1−x2+Fe2x/33+□x/3(HCOO)3 (GUA+ = guanidinium cation) at xord ≃ 0.6.51
The sense one gets is that, other than for very dilute mixtures, truly-random compositional distributions are conspicuously rare in MOFs. Clearly, kinetics of crystallisation play an important role here; for example, the different nucleation rates of Zn- and Co-containing imidazolate frameworks allow the one-port formation of mixed-metal frameworks with core–shell architectures.57 Likewise, precursors with different kinetic stabilities can be chosen cleverly to control the degree of mixing during synthesis of mixed-metal coordination polymers.58 But there are thermodynamic considerations at play too: a comprehensive DFT study has highlighted the dual roles of charge distribution and size-mismatch in determining the balance between mixing and segregation.59,60 Such considerations may rationalise the preferential formation of CeZr5 clusters in mixed-Ce/Zr UiO-66 samples, as evidenced by X-ray absorption fine structure (EXAFS) measurements and as implicated in their catalytic behaviour.61
In addition to DFT calculations and EXAFS measurements, a variety of techniques offer sensitivity to the spatial distribution of components in compositionally disordered MOFs. Electron microscopy allows direct visualisation of compositional gradients on the micron scale [Fig. 3(a)]; so too does fluorescence imaging microscopy.64 Even peak asymmetry in synchrotron powder X-ray diffraction measurements has some indirect sensitivity to this microstructure.65 Rotational-echo double-resonance (REDOR) NMR measurements are sensitive to the through-space distance between nuclear spins; by selective labelling of different linkers it is possible to differentiate random, clustered and anti-clustered mixed-linker arrangements on the unit-cell lengthscale. A particularly elegant REDOR-NMR study of a multivariate MOF-5 family showed how different linker mixtures favoured different substitutional correlation types [Fig. 3(b)].54 Very recently, atom probe tomography (APT) has even been used to ‘sequence’ metals in mixed-metal MOF-74 rods [Fig. 3(c)], demonstrating that both random and correlated sequence types emerge naturally from different metal combinations and different synthesis temperatures.62 While both REDOR-NMR and APT provide very detailed information, they are nonetheless demanding techniques. By contrast, we have found that humble infrared spectroscopy can offer a surprisingly straightforward measure of cation distributions in mixed-metal MOFs and coordination polymers [Fig. 3(d)].63
Fig. 3 Determining compositional correlations in compositionally disordered MOFs. (a) STEM-EDX images of transition-metal segregation in the mixed-metal Zn/Ni/MnCID-5 ⊃ G MOF.58 Adapted from ref. 58 with permission from the American Chemical Society. (b) REDOR-NMR decay curves are sensitive to the degree of clustering/anticlustering of individual linkers in mixed-linker MOFs (LC = large clusters, SC = small clusters, Alt = alternating).54 Adapted from ref. 54 with permission from the American Association for the Advancement of Science. (c) Atom probe tomography allows the sequencing of metal-oxide rods in multivariate MOF-74.62 Here the method is shown to differentiate ‘duplicate’ (D), random (R), and ‘insertion’ (I) sequences. Adapted from ref. 62 with permission from the American Association for the Advancement of Science. (d) The in-plane methylimidazolate (mIm) stretching region of the IR spectra of mixed-metal (Zn/Cd) ZIF-8 derivatives is sensitive to the relative populations of Zn–mIm–Zn, Zn–mIm–Cd and Cd–mIm–Cd linkages. These populations can be extracted using non-negative matrix factorisation and used to derive models of the Zn/Cd distribution, which is found to involve clustering.63 Note the greater-than-statistical population of Zn4/Cd4 four-rings (experimental ZnnCd4−n populations as solid bars; statistical populations as open bars). Adapted from ref. 63 with permission from the Royal Society of Chemistry. |
Just as enzymes derive their function from the precise spatial arrangement of different components—e.g. the amino-acids of their protein backbone—the vision with compositionally complex MOFs has always been one of controlling heterogeneity to optimise a particular chemical or physical property.11 A simple but effective example is the recent demonstration of correlation-dependent photocatalysis in a mixed-linker UiO-68 derivative.66 But this vision is not unique to the MOF field: indeed, conventional solid-state chemistry has long exploited the use of solid solutions and inhomogeneity in this general way. A good example is that of the relaxor ferroelectrics, such as PbMg1/3Nb2/3O3, where the structural complexity responsible for their useful dielectric properties emerges naturally from an interplay between composition and lattice geometry.67 We anticipate the design rules developed to control complexity in these inorganic materials may yet have a role to play in MOF science.12
Fig. 4 Orientational disorder in MOFs. (a) The α- and β-polymorphs of the substituted MOF-5 derivative IRMOF-3-AMPh have different linker orientation correlations and H2 (de)sorption profiles. Using the fully substituted AMPh-terephthalate linker during MOF synthesis gives the ordered α phase, with reduced and hysteretic H2 uptake is sorption profile. By contrast, post-synthetic modification of an aminoterephthalate MOF-5 derivative gives the disordered β phase, which has a greater accessible pore volume.71 Adapted from ref. 71 with permission from the Royal Society of Chemistry. (b) In the pillared MOF DUT-8(Ni), naphthalenedicarboxylate linkers introduce a ‘step-up’/‘step-down’ between neighbouring columns (exaggerated for effect here). There are six feasible NDC arrangements—i.e. two steps up and two steps down—around the lozenge-shaped pores. These arrangements give pores with either D2d or C2v point symmetry (grayscale and coloured tiles, respectively) that propagate according to a set of simple rules (jigsaw pairings shown here). Different guest molecules have different affinities for these two pore geometries, which allows a guest-induced transition between complex phases; for example, toluene and dichloromethane (DCM) favour configurations with C2v- and D2d-rich pore geometries, respectively.15 |
Just as linker orientations influence pore chemistry, so is it that guest adsorption can affect the correlations in linker orientations. Perhaps the clearest example in a disordered MOF is found in DUT-8(Ni), a pillared paddle-wheel MOF with 2,6-naphthalenedicarboxylate (NDC) linkers.72 The origin of disorder in this system is the existence of two degenerate NDC orientations.73 Each orientation introduces an offset—either up or down—in the positions of the neighbouring paddlewheel pillars. Since the lozenge-shaped pores of DUT-8(Ni) are flanked by four NDCs, there is a natural constraint that two must step up and two step down—after all, the total offset on traversing a loop around the pore must vanish to zero.68 The six possible NDC arrangements satisfying this rule partition into a set of four with a C2v-symmetric pore structure, and a set of two with a D2d pores. Neighbouring pores share edges, such that a large and complex configurational landscape emerges, characterised by different spatial arrangements of pore geometries.15 Remarkably, the adsorption and desorption of different solvent molecules drives the system reversibly across this landscape—the strength of host/guest interactions overcoming the ∼20 kJ mol−1 barrier to NDC ‘flips’ [Fig. 4(b)]. The recent development of efficient and accurate MOF force fields24,25 may allow computational insight into the mechanism for this unusual disorder–disorder switching.
The extent to which orientational disorder is static or dynamic depends on the reorientation barrier height, and if this barrier is commensurate with the available thermal energy then one expects to observe an order/disorder transition.12 There are plenty of examples in the literature, but we focus here on the case of the hybrid perovskite [NPr4][Mn(N(CN)2)3]—not technically a MOF but relevant nonetheless because its MOF-like framework structure exhibits a particularly useful barocaloric effect exploitable in solid-state cooling.74 The conformations of the flexible [NPr4]+ extra-framework cations and orientations of the dicyanamide linkers couple: they are disordered at temperatures above 331 K and ordered below. The transition between the two states is accompanied by a large change in both entropy and volume. The latter confers sensitivity to pressure—the disordered state can be re-ordered by squeezing, which is the basis for a barocaloric heating/cooling cycle [Fig. 5].75 Correlations within the disordered state, although still poorly understood, influence the barocaloric entropy change and reduce the critical pressure needed for each cooling cycle.12,74 In the parallel field of magnetocalorics—based on the order and disorder of magnetic spins—correlated disorder and the geometric frustration from which it arises are well-known to increase cooling efficiency.76–78 We therefore anticipate that similar design principles will prove valuable in exploiting orientational disorder in hybrid frameworks and MOFs alike.
Fig. 5 A simple barocaloric cooling strategy that exploits the pressure-dependence of the order/disorder transition in [NPr4][Mn(N(CN)2)3]. Starting at ambient conditions where the system is orientationally disordered (top left), ordering is induced by (1) application of pressure (note the change in unit cell dimensions is exaggerated here for effect). The loss of configurational entropy results in heating, and this excess heat (2) is dumped to the environment. (3) Removing the applied pressure forces the system to disorder; the concomitant increase in entropy now requires the system to cool. The cold solid (4) extracts heat from its environment.74,75 Adapted from ref. 74 with permission from Springer Nature. |
An historical barrier to studying and exploiting correlated disorder in MOFs has been the ostensible insensitivity of experimental techniques to its presence and nature: conventional crystallography is the obvious culprit here. But there are now many different experimental and computational tools available with the relevant capabilities. The hurdle that unquestionably remains is that of developing rigorous strategies for exercising synthetic control over correlated disorder. To this end we have flagged the interplay between composition and network geometry, and also the role of local-symmetry lowering as potential mechanisms for driving very specific types of correlated disordered states. The balance of kinetic and thermodynamic control during synthesis is clearly also to be decisive in many cases.
Looking forward, the field shares much of its aims with those of multivariate MOF chemistry: ultimately, correlated disorder is a form of structural complexity that might be exploited in e.g. catalysis or information storage or selective host/guest chemistry.11,68 An intriguing recent development along similar lines is the concept of using correlated disordered states as intermediates in materials synthesis to target low-dimensional phases.17 But we note in conclusion the possible additional importance of different types of correlated disorder in MOFs to those outlined here—involving e.g. frustrated magnetism,78 local charge ordering,81 or cooperative spin-state transitions.82 We see enormous potential in combining the unconventional physics of correlated disordered states with the scope for exquisite control over geometry for which MOFs are renowned.
This journal is © The Royal Society of Chemistry 2021 |