Evidence for a chemical clock in oscillatory formation of UiO-66.

Chemical clocks are often used as exciting classroom experiments, where an induction time is followed by rapidly changing colours that expose oscillating concentration patterns. This type of reaction belongs to a class of nonlinear chemical kinetics also linked to chaos, wave propagation and Turing patterns. Despite its vastness in occurrence and applicability, the clock reaction is only well understood for liquid-state processes. Here we report a chemical clock reaction, in which a solidifying entity, metal-organic framework UiO-66, displays oscillations in crystal dimension and number, as shown by X-ray scattering. In rationalizing this result, we introduce a computational approach, the metal-organic molecular orbital methodology, to pinpoint interaction between the tectonic building blocks that construct the metal-organic framework material. In this way, we show that hydrochloric acid plays the role of autocatalyst, bridging separate processes of condensation and crystallization.

Achemical clock is a particular type of chemical reaction that is frequently applied in chemistry education. The precise definition of a chemical clock is cause for some debate1, but it is generally accepted that the reaction possesses an induction time after which a sudden and marked change in concentration of the reactant and product species ensues2. A reaction based on hydrogen peroxide, iodide and thiosulfate is popular3:

The first reaction is slow, the second one fast. As a result of the difference in reaction rates and the fact concentrations of the autocatalytic redox-pair iodide/iodine are coupled in nonlinear kinetics, the iodine concentration is seen to shoot up after an initial lag time after which it oscillates to equilibrium. Iodide is not consumed in the overall reaction and serves as autocatalyst. If starch is added, the oscillations can be followed by corresponding changes in the solution colour. This makes it a popular classroom reaction. In broad context, clock reactions belong to a class of nonlinear chemical phenomena linked to several natural processes and rhythms. Great interest exists in coupling oscillatory clock reactions to processes occurring in the solid state4. A particularly exciting development is the use of chemical actuation to fabricate a synthetic muscle. It is possible to use a clock reaction to produce oscillatory pH or temperature trajectories that drive a muscle-like mechanical response from a polymeric material5. For example, Howse et al.6 fabricated in this way a device capable of producing 20 mW kg−1 of power. These are solution-driven oscillations in the solid state, but ‘direct' oscillatory behaviour exists in conversion and temperature during certain polymerization reactions7. Clock reactions involving phase transitions, in which for instance precipitating or crystallizing particle dimensions oscillate during the course of reaction, have never been shown to exist, despite nonlinear dynamics suspected to govern several important condensation processes, such as the hardening of cement8.

In crystal engineering studies—the design and synthesis of functional crystalline structures based on a bottom-up approach from smaller building blocks—scientists have shown remarkable chemical mechanisms to prelude and govern synthesis of metal–organic framework (MOF) crystals910111213141516. MOFs are porous materials consisting of metal clusters that are linked to each other by organic bridging ligands. MOFs carry, arguably more than any other type of material, the potential to fulfil crystal engineering, for the cluster and ligand connectivity, and properties, can be to certain extent controlled17. This control is limited, much because of the unpredictability of coordination chemistry around the tectonic building blocks of MOFs1418. One example is a report by Farha et al. on the synthesis of the ultra-stable and highly popular MOF UiO-66, which is based on hexanuclear zirconium/hafnium oxoclusters, and terephthalate. The authors showed that the addition of hydrochloric acid (HCl) considerably speeds up the synthesis19. This is surprising, because HCl is formed during the synthesis, which (based on naive grounds) would suggest that its neutralization, rather than addition would lead to improvement in synthesis by favourably shifting the equilibrium towards the formation of the MOF. The only scenario, in which HCl is formed in situ, while being kinetically beneficial matters the scenario of HCl being autocatalytic, fulfilling a similar role to iodide in reactions 1 and 2, and possibly sparking spectacular phenomena during the formation of UiO-66.

Using a combination of in situ synchrotron X-ray scattering and a quantum chemical methodology, we herein show that HCl is indeed autocatalytic; as a consequence, the formation of UiO-66 involves a clock reaction with a lag time followed by an oscillatory wave front in crystal dimensions.

Results

X-ray scattering

In situ small- and wide-angle X-ray scattering (SAXS/WAXS) is well suited to follow crystallization and condensation processes. If synchrotron radiation is used, time-resolved data of high resolution can be obtained on multiple length scales of self-assembly. The small-angle camera follows events on the colloidal (2–50 nm) length scale, permitting time-resolved analysis of the solidifying entity, while the wide-angle camera follows the evolution of Bragg peaks in case of evolution of a crystalline compound16. Important, and of great use in this work, is the possibility to deconvolute condensation from crystallization by respectively analysing SAXS versus WAXS data.

The SAXS/WAXS measurements were carried out at the X9 beamline at the National Synchtrotron Light Source (NSLS), Brookhaven, NY, USA. SAXS data were processed to obtain temporal information on the following integral parameters:

·Porod volume, Vp(nm3): volume of the particle.

·Radius of gyration Rg (nm): root mean square average of all distances from the particle boundary to its centre of mass.

·Number densiy N/V (−): the number of particles per unit volume, in arbitrary units. Only the trend is followed.

·Porod surface-to-volume S/V (nm−1): specific surface area of the precipitating particle.

·During the experiments, WAXS was used to follow:

·Extent of crystallization α (−): extent parameter obtained via Bragg peak evolution. In this case focused on the most dominant reflection of UiO-66 at q=0.5 nm−1.

Figure 1 displays the time derivative of the calculated Porod volume, Vp′ (or dVp/dt), versus time. Vp′ is used as a parameter for the particle growth rate. With Zr in blue and Hf in red, the evolution is displayed for the slower syntheses; ½ eq. HCl at 100 °C and 1 eq. HCl at 80 °C (eq. concerns a defined equivalent of HCl specified in the Methods section). We see that after a certain induction time and rapid condensation thereafter, the particle size oscillates to final values, that is, the derivative in Fig. 1 oscillates to zero. Oscillations were, in fact, found for the trends of all integral parameters (Supplementary Figs 10–14), and fully reproducible over several numbers of experiments. There is a difference between Zr and Hf in terms of dependence of maximum growth rate on temperature and HCl concentration. For Zr, decreasing the temperature by 20 °C and doubling the amount of HCl leads to a decrease of this maximum, whereas for Hf, there is only a small change and it follows the opposite direction. In further analysing the patterns, we found the evolution of the Porod volume to mimic the response of an underdamped harmonic oscillator. All trends, with exception of the fastest syntheses at 120 °C or 2 eq. of HCl, could be well fitted to V(t)=exp[−γt]·a·cos(ωt−α), where γ (min−1) is the damping coefficient, ω (min−1) the frequency and a (−) and α (min−1) amplitude and offset constants, respectively (see Supplementary Methods and Supplementary Fig. 1 for a more detailed explanation). What follows directly from our mathematical analysis is that syntheses involving higher temperatures and/or higher concentrations of HCl have their oscillations more effectively damped, and induction times are decreased. This is expected for a chemical clock: increasing the concentration of the autocatalytic species, generally, carries the same effect as increasing the temperature: the rate of the overall reaction goes up. It is interesting to check the response of oscillation frequency ω to changes in temperature and HCl concentration for both metals, by comparing it to the natural frequency ω0, which is the frequency of the undamped system. In doing so, Supplementary Fig. 9, we see that Zr-based UiO-66 reacts more strongly to an increase of temperature and HCl concentration with respect to Hf, and its oscillations are more easily damped. So, there are differences between Zr and Hf. This, and certainly Hf's sluggish behaviour with respect to Zr, is not new in coordination chemistry and has been observed before20.

Figure 1

Particle growth rates.

Time-resolved plot of particle growth rate, as monitored by Vp′ (nm3 min−1) for zirconium (blue) and hafnium (red) based UiO-66, at the conditions of ½ eq. HCl and 100 °C, and 1 eq. HCl and 80 °C.

Let us go back to the oscillations of all integral parameters; we remark that the final, equilibrium values for particle dimensions (Vp, Rg, S/V) in synthesis were found to be independent of temperature and or HCl concentration, consistent with scanning electron microscope analysis on the UiO-66 crystals by Farha and colleagues19. Notably, no oscillatory behaviour was observed for the extent of crystallization α obtained from the WAXS data, which is in line with observations by Ragon et al.21 in a study on formation of UiO-66 with in situ synchrotron X-ray diffraction. Moreover, the start of bulk crystallization was in all cases seen to carry a small delay, while proceeding gradually with respect to the process of condensation. This indicates that the condensed coordination polymer is briefly amorphous before rearranging towards a crystalline structure.

It is now important to explain the induction time and ensuing oscillatory behaviour, which is something expected for clock reactions, such as the clock reaction discussed in the introduction. In that reaction, iodide served as autocatalytic species serving an overall reaction without being consumed. The only species in the synthesis of UiO-66 that could carry such a role is HCl. A starting point in explaining the influence of HCl lies in the observation that as-synthesized UiO-66 contains four μ3-OH ligands, next to four μ3-O ligands. This is striking, because the μ3-OH ligand is known to be highly Brønsted acidic. Indeed, in the case of UiO-66, it has been used to play either an active role in acid catalysis22, or be prone to exchange for other monovalent moieties in targeted functionalization of the inorganic node23. These ligands are thus expected to materialize in rather acidic local environments, which implies that HCl might be actively protonating the inorganic cluster during the synthesis of UiO-66. This by itself does not explain the clock-like phenomenon, for which HCl needs to be released by the crystallizing solid during the oscillations. But since we know that only four out of eight bridging oxygen ligands are protonated in the final framework for charge-balancing reasons, this would imply that in the short period during which the solid is amorphous, the number of protonated oxygen ligands could exceed the number of four. Thus, the question is whether protonated inorganic clusters may be vital during UiO-66 synthesis. In fact, this would not be a major surprise, for the decoration of building blocks by small moieties to promote MOF topologies is commonplace. For instance, we have shown by quantum chemical calculations before that the addition of an amine group to terephthalate is important for the ligand's ability to coordinate to the μ3-O-centred Al3 tectonic unit that builds NH2-Material Institute Lavoisier (MIL)-101(Al)15. Despite confirmation of experiment by computation, this, and analogous cases, have remained unexplained from the perspective of chemical theory, for calculations have merely focused on quantitative determination of energies instead of qualitative understanding14.

As Kohn–Sham density functional theory (DFT) is a molecular orbital (MO) method, it is insightful to picture the metal–organic bond as an MO interaction between an inorganic and an organic fragment, the nature of which we will try to unravel physically. We introduce the metal–organic MO (MOMO) analysis, which we believe will be of general use in providing chemically intuitive interpretation of DFT computation on MOF self-assembly. Herein, it is employed to investigate metal–organic interaction as function of degree of protonation of the inorganic unit.

Density functional calculations

All calculations were carried out with Amsterdam density functional (ADF) at the ZORA-Becke-Lee-Yang-Parr (BLYP)-D3(BJ)/TZ2P level of the generalized gradient approximation (GGA) of Kohn–Sham DFT242526, where ZORA stands for zeroth-order regular approximation, used to take relativistic effects into account, expected to be particularly important for calculations involving hafnium2728. More details are given in the Methods section. In this first use of the MOMO method, a fragment of 12 geometrically optimized formate molecules is deformed to its final, crystallographic state in the UiO-66 framework using the activation-strain model29. The final crystallographic state is derived from Rietveld refinements provided by Lamberti et al.3031 (Supplementary Table 3). Formate, rather than terephthalate or benzoate, was chosen to make the high-level calculations less expensive. The specified [formato]12 fragment is allowed to interact with an inorganic fragment of the form M6O8Hx, with M=Zr or Hf, and x=0–8, rendering a total of 2 × 9 isoelectronic (inorganic) fragments, and 18 metal–organic interactions with the degree of protonation as main parameter in the construction of Zr- and Hf-based UiO-66. As explained above, it is our target to obtain a physical picture of the bonding situation.

We used an energy decomposition analysis (EDA) to analyse the metal–organic interaction. This allows us to split the energy associated with interaction into electrostatic, orbital and Pauli terms: ΔEint=ΔVelstat+ΔEoi+ΔEPauli (Fig. 2a)32.

Figure 2

Activation strain model and degree of covalency.

(a) Implementation of the activation-strain model and energy decomposition analysis (EDA), in which the [formato]12 fragment reacts with 18 isoelectronic M6O8Hx (x=0..8) fragments. ζ represents the reaction coordinate. The fragments deform towards their final state in the framework, and this metal–organic interaction is decomposed into energetic terms for electrostatic interaction, orbital interaction and Pauli repulsion. (b) Degree of covalency defined as ΔEoi/ΔVelstat as calculated for the 18 fragment interactions (nine for zirconium and nine for hafnium based UiO-66).

As expected, the bonding becomes stronger upon protonation of the inorganic fragment. This is anticipated, because it involves increase of charge difference between the fragments. But it is here that EDA manifests itself as a valuable tool by providing the possibility of calculating a measure of covalency, herein defined as ΔEoi/ΔVelstat. This term is seen to increase linearly and significantly upon increasing the degree of protonation of the inorganic fragment (Fig. 2b). Clearly, highly protonated clusters engage in stronger orbital interactions and precipitate more irreversibly. We relate this to the additional strengthening of the bond, and arguably also the reduction in polarity that prevents solvent attack on the bond, and thus also dissociation and solvation. Mulliken populational analysis (Supplementary Table 2) also indicates that electron transfer from the [formato]12 fragment increases upon increasing the number of protons on the (electron-accepting) inorganic fragment.

The MOMO approach provides insight in disclosing the nature of this increase in metal–organic donor–acceptor interaction. We first recognize that the bonding between the two fragments involves 36 electron pairs (24σ and 12π, in line of expectation for 12 μ2-formato ligands). We then note that in many cases the orbital interactions between the two fragments involve more than two fragment orbitals, which makes it rather challenging to construct a simple MO picture that provides a good qualitative picture of the influence of protonation on the metal–organic bond formation. However, if one only considers the interactions involving the inorganic fragments with zero and eight protons, one does not only consider the extremes that span the range of study, but also the interactions in which both interacting fragments carry O symmetry. This considerably facilitates the analysis: the use of group theory indicates that solely considering the high symmetries A1g, A2g, A2u, thus three MO diagrams per model complex, a reasonably complete picture can be obtained. Only the four MOs of totally symmetric A1g are displayed in Fig. 3, since they represent the most dominant interactions (MO diagrams for A2g and A2u can be found in Supplementary Figs 15–23 and Supplementary Table 1). From the MO coefficients displayed in Fig. 3, it is clear that for both Zr- and Hf-based UiO-66, metal–organic bonding becomes decreasingly polar as the degree of protonation increases. Note that the stronger bonding is not a result of the metal–organic orbital overlap, which remains similar if not slightly decreasing for the eight proton inorganic clusters with respect to their zero proton analogues.

Figure 3

MOMO diagrams.

MOMO diagrams showing the metal–organic interaction between the totally symmetric (A1g) fragment orbitals, with the inorganic fragment either carrying zero or eight protons, in O symmetry. The overlap matrices are depicted in the top left corner for each diagram.

The analysis of the metal–organic bonding thus indicates that while bonds are spatially similar, covalency—as we measured it—increases. In a simple description of this phenomenon, we consider the postulate of the bond-order conservation principle, without the Morse-potential approach as frequently applied by Shustorovich and others33, and observe the inorganic cluster that undergoes protonation at μ3-O. According to Shustorovich' theory this weakens the metal–oxygen bond and therefore increases the acidity of the metal ion, thus enhancing the metal–organic interaction in this case.

Discussion

Figure 4 displays a schematic pathway based on the observations in SAXS/WAXS and DFT analysis. Step 1 in the scheme occurs in the preparation of the precursor solution and concerns the water-promoted hydrolysis of ZrCl4 and formation of zirconyl chloride species, which are known to form multinuclear species in solution, such as [Zr4(μ2-OH)8]8+34. The rearrangement of such clusters into the hexanuclear UiO-66 building block follows an unknown path, but what is certain is that the pathway to condensation must be both water and acid catalysed, since in the appearance of μ3-O(H) bridges, O's remain largely protonated. This is indeed observed by the groups of Serre and Farha1921. The highly protonated clusters then engage in metal–organic bond formation in step 2 and fast condensation of the solid coordination polymer ensues. The solid is not fully crystalline at this brief stage, and the situation is metastable. Reorganization however, evidenced by SAXS versus WAXS comparison, is swift: the crystalline UiO-66 lattice is formed and the oxido bridges deprotonate, leaving the scaffold charge-neutralized (step 3). This process of bulk crystallization within the solid state deserves some notion. By virtue of established laws in solid-state chemistry, such rapid transition requires a large degree of chemical inheritance between parent and resulting solid, that is, involving no substantial bond breakage and reformation35. Thus, in line with works by Cheetham36, we expect the solid to exist as amorphous MOF, chemically cognate to UiO-66, but with protonated oxido bridges countered by chloride ions in solution.

Figure 4

Schematic view of the formation of UiO-66 and the autocatalytic role of H+.

Step 1: formation of multinuclear clusters; step 2: consumption of H+ and fast condensation; and step 3: development of crystallinity upon release of H+.

The drive behind the process of crystallization is in all likelihood caused by lattice enthalpy, similar to what we observed before in the crystal-to-crystal transformation of NH2-MOF-235(Al) to NH2-MIL-101(Al)15, and also seen in the formation of zeolites, where charged templates are squeezed out of the network during the process of crystallization37. The release of protons in step 3 is vital, since it accelerates step 2, which is prohibited from moving forward, and even reversed by the low concentration of protons after their consumption. This is witnessed by the negative values for V′ in Fig. 1. The oscillating proton concentration, caused by the coupling of step 2 and step 3, thus leads to oscillating temporary equilibria in step 2, and dissolution recrystallization near the surface of the particle governs commensurate oscillations for all integral parameters.

Despite this being the first report on a clock reaction that generates oscillations in the condensed state, we do not rule out the possibility that similar processes occur in the formation of other MOFs; on the contrary, we expect them to be present in the coordination of other high-coordinate oxoclusters.

Again, this work demonstrates that in resolving molecular pathways of MOF crystallization, a multi-scale approach is required: in this work, the use of both SAXS and WAXS was required to separate the process of the condensation from crystallization and show them to exist as independent processes in MOF synthesis. At the molecular scale, the MOMO Kohn–Sham DFT approach extends the information obtained with X-ray scattering to provide necessary insight from the perspective of coordination chemistry. We believe that the MOMO approach should be encouraged in theoretical studies on MOF formation, since interaction between the organic struts and inorganic nodes can be highlighted, essential in the understanding of metal–organic network formation. It is herein demonstrated that this is a viable new theoretical approach, while overall, the result adds to the exciting array of phenomena associated with nonlinear dynamics.

Methods

Synthesis

An in-house developed SAXS/WAXS cell (Supplementary Figs 2 and 3), further described elsewhere16, was used for the experiments. For a typical solution, 0.75 mmol terephthalic acid (134 mg), 0.54 mmol anhydrous ZrCl4/HfCl4 (125/172 mg), 1 ml HCl (37%) and 15 ml dried dimethylformamide (DMF) were mixed and injected in the synchrotron cell. A measure of 1 ml stands for 1 eq. of HCl in the text, added to the starting solution.

SAXS/WAXS measurements

Time-resolved SAXS/WAXS experiments were performed at the X9 beamline at the NSLS at Brookhaven National Laboratory (New York, USA). The undulator-based X9 beamline at the NSLS applied a marCCD SAXS detector and photonic science WAXS detector. The wavelength was in all cases 0.879 Å, and the sample-to-detector distance was 3.099 and 0.4500, m for the SAXS and WAXS camera, respectively. In the cell, the synthesis solution is loaded between two mica windows separated by polytetrafluoroethylene (PTEF) spacers (1.5 mm thickness). The depth of the solution probed by X-rays in measurement was ca. 75 mm. The beam-spot size was 20 by 20 μm. Heating was provided by four cartridge electrical resistances and temperature was controlled at the external wall of the PTFE inserts with a temperature controller.

In calculation of the integral parameters, the Guinier regime for the calculation of Rg was for all data sets based on the scattering vector range of q=(0.07–0.15), the Porod regime was in all cases q=(3.2…4). The number density is based on both the Guinier and Porod regime.

All data sets were background subtracted, azimuthally integrated and normalized for beam intensity to yield time-resolved one-dimensional data that was processed to obtain temporal information on the following integral parameters: radius of gyration Rg (nm),

Porod volume Vp (nm3)

with Porod invariant Q

Porod surface-to-volume S/V (nm−1)

and number densiy N/V (−)

Where forward scattering, I(q=0) was calculated using (3) to extrapolate q to 0.

Extent of crystallization α was derived from the most intense Bragg peak in UiO-66, at q=0.5 nm−1. Supplementary Figures 3–8 show the in situ collected diffraction data.

DFT calculations

All standard calculations were performed using the ADF program27283738, using the BLYP functional. The numerical integration was performed using the procedure developed by Boerrigter, te Velde and Baerends3940. The MOs were expanded in a large uncontracted set of Slater type orbitals containing diffuse functions, which is of triple-ζ quality for all atoms and was augmented with two sets of polarization functions (TZ2P). The core shells of all atoms were treated by the frozen-core approximation41. All energies and geometries were calculated using the GGA of DFT at the BLYP level4243. GGA proceeds from the local density approximation, where exchange is described by Slater's Xα potential and correlation is treated in the Vosko–Wilk–Nusair parametrization44, which is augmented with nonlocal corrections to exchange due to Becke45, and correlation due to Perdew46 added self-consistently47. For all calculations, relativistic effects were accounted for by ZORA48 and dispersion by Grimme's correction, D345. The (small) contribution from dispersion, obtained via EDA can be found in the Supplementary Information.

The MOMO approach is based on the activation-strain model, elaborated on elsewhere32495051. Activation-strain model was only applied to the organic fragment, since the unprotonated inorganic fragment is not stable upon geometry optimization. To verify that formate ligands are representative in calculation, an additional (expensive) calculation with four protons and benzoate fragments was carried out at the same computational level.

Data availability

The authors declare that all other relevant data not included in the Supplementary Information and supporting the findings of this study are available on request.

Additional information

How to cite this article: Goesten, M. G. et al. Evidence for a chemical clock in oscillatory formation of UiO-66. Nat. Commun. 7:11832 doi: 10.1038/ncomms11832 (2016).