Open Access Article
Loris Lombardo
*ab,
Ellan K. Berdichevsky
c,
Farzaneh Fadaei-Tirani
d,
Shin-ichi Orimo
ae,
Andreas Züttel
b and
Satoshi Horike
*fgh
aInstitute for Materials Research, Tohoku University, Sendai, Miyagi 980-8577, Japan. E-mail: lombardo.loris.giovanni.a5@tohoku.ac.jp
bInstitute of Chemical Sciences and Engineering, Basic Science Faculty, École Polytechnique Fédérale de Lausanne (EPFL) Valais/Wallis, Energypolis, Rue de l’Industrie 17, 1951 Sion, Switzerland
cDepartment of Synthetic Chemistry and Biological Chemistry, Graduate School of Engineering, Kyoto University, Katsura, Nishikyo-ku, Kyoto 615-8510, Japan
dInstitut des Sciences et Ingénierie Chimiques, École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
eAdvanced Institute for Materials Research (WPI-AIMR), Tohoku University, Sendai, Miyagi 980-8577, Japan
fDepartment of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa-Oiwakecho, Sakyo-ku, Kyoto 606-8502, Japan. E-mail: horike.satoshi.3r@kyoto-u.ac.jp
gInstitute for Integrated Cell-Material Sciences, Institute for Advanced Study, Kyoto University, Yoshida-Honmachi, Sakyo-ku, Kyoto 606-8501, Japan
hDepartment of Materials Science and Engineering, School of Molecular Science and Engineering, Vidyasirimedhi Institute of Science and Technology, Rayong 21210, Thailand
First published on 6th April 2026
Zirconium metal–organic cages (Zr-MOCs) exhibit promising applications in adsorption and molecular separation, thanks to their high porosity, stability, and solution processability. However, their synthesis relies on traditional solvothermal methods. In this work, we developed a mechanochemical synthetic route for Zr-MOCs by milling pre-formed Zr-clusters with terephthalic acid. The Zr3-cluster is easily prepared by combining zirconocene dichloride with carboxylic acids such as acetic acid and benzoic acid. Typical yields of 70% were achieved in only 30 minutes, using a minimal amount of DMF (η = 0.5 µL g−1). This kinetically controlled strategy produces phase-pure tetrahedral cages with high surface area (>450 m2 g−1), comparable to those produced by conventional solvothermal processes. The method enables tunable cage architectures by varying ligand identity. The successful synthesis of Zr-MOCs with mono- and bi-functional bidentate ligands, as well as tri-dentate ligands, was achieved, demonstrating the versatility of the mechanochemical approach. This work opens new possibilities for the preparation of Zr-MOCs with insoluble ligands and facilitates scalable production.
Zr-MOCs commonly contain the [Cp3Zr3O(OH)3(COO)3]− (Cp = η5-C5H5) secondary building unit (SBU), and can exhibit different cage geometries when built with bidentate ligands, such as tetrahedrons (four SBUs per cage) or “cigars” (two SBUs per cage). Recently, several approaches have been explored to control the geometry and phase-purity of Zr-MOCs. Ligand-based phase control limits the choice of ligands to bulky ones for the formation of pure tetrahedral cages, or rigid, extended ligands for pure “cigar” cages.13 The type of solvent can also direct the formation of the tetrahedral or “cigar” structure, but it was demonstrated only for selected ligands.14 The separation of mixtures containing both Zr-MOC geometries has been demonstrated, although the process is laborious and has poor atom economy, requiring large amounts of solvent and AgOTf to solubilize the cages.15 These limitations underscore the necessity for innovative approaches that allow for the controlled and scalable synthesis of Zr-MOCs with well-defined geometries.
Mechanochemistry has emerged as a powerful green synthetic method for MOFs and other porous materials.16–22 This technique is environmentally benign (minimal amounts of solvent), allows short reaction times, and is scalable. Mechanochemical synthesis has proven to be effective for preparing phase-pure MOFs by limiting side reactions.23 It has been extensively studied for the preparation of UiO-type MOFs.24,25 Most of the reports first prepared the Zr-SBU before milling it with the target ligand.26–29 In the case of Zr-MOC, the SBU [Cp3Zr3O(OH)3(COO)3]− can be prepared using different carboxylic acids such as benzoic acid or acetic acid.12,30,31 Once the Zr-SBU is ready, ligand substitution in solution was demonstrated as a pathway to get Zr-MOC.32
Inspired by the work on UiO MOF synthesis using Zr-SBUs, we report a kinetically guided strategy that combines pre-formed Zr-cluster with bi-dentate and tri-dentate ligands for the preparation of tetrahedral Zr-cages in only 30 minutes. Mechanochemical synthesis can be conducted with or without DMF and exhibits broad tolerance toward diverse functional groups. Yields exceeding 70% have been achieved, and all resulting cages display high porosity. This synthetic approach may be applicable to insoluble ligands and is ready for scale-up.
space group. The structure is composed of trinuclear Zr-clusters with three BA coordinated, as expected. The clusters interact with each other via hydrogen bonding and weak interactions, forming a cage-like structure composed of two clusters (Fig. S3).
The crystallographic data is presented in Table S1. Both clusters exhibit thermal stability up to 225 °C and 245 °C for ZrCp-BA-cluster and ZrCp-OAc-cluster, respectively (Fig. S4).
Terephthalic acid (H2BDC) was selected as the ligand for mechanochemical optimization using ZrCp-OAc-cluster, as the corresponding cage has been reported to form cigar- or tetrahedral-type geometries.13 Milling was carried out at 500 rpm with DMF as the LAG solvent (η = 0.5 µL g−1), chosen for its ability to dissolve both precursors and its common use in solvothermal synthesis. Reaction times of 15–60 min were tested with 5 minute pauses every 15 min. The obtained Zr-MOC-BDC powders were washed and activated before measuring N2 isotherms. Increasing the time from 15 to 30 min improved both yield and surface area of Zr-MOC-BDC, giving >70% yield and a BET surface area of 509 m2 g−1 (Fig. 2A and S5), comparable to phase-pure Zr-BDC-tetrahedral cages prepared solvothermally (462–685 m2 g−1).32–34 In contrast, Zr-BDC-cigar cages typically show surface areas <100 m2 g−1.13 Mechanochemical synthesis can effectively provide the kinetic product (tetrahedral) thanks to the short reaction time, mild energy input, and low amount of solvent.35,36 Prolonged milling (60 min) lowered the yield by ∼20% and surface area by ∼40%. FT-IR spectra confirmed Cp* presence with sharp peaks at ∼809 and ∼1016 cm−1 from the Zr-SBU in all samples (Fig. 2B and S5). However, the 60 min sample displayed a stronger band at 1655 cm−1, consistent with the presence of coordinated DMF.37 Prolonged milling thus compromises cage integrity through the formation of defects. The rapid ligand exchange between OAc and BDC was confirmed by examining the PXRD patterns after 5 minutes of milling (Fig. S6 and S7). While BDC retains its crystalline structure when milled alone, the product becomes weakly crystalline when BDC is milled with ZrCp-OAc, confirming the occurrence of spontaneous ligand exchange. This behavior is readily understood, as the cage is thermodynamically more stable than the cluster.12,32 In addition, milling the cluster alone can re-arrange the acetate binding from bidentate to monodentate, promoting the exchange with BDC, consistent with the ligand exchange mechanism of Zr12 oxo-cluster (Fig. S8).38 With the milling time fixed at 30 min, DMF content was varied. Across all conditions, the BET surface area remained >450 m2 g−1 (Fig. 2C and S9). Notably, cage formation under neat grinding is uncommon in ball-milled MOFs, suggesting that acetic acid released during milling may act as an intrinsic LAG agent. The yield increased at η = 0.5 µL g−1 DMF but declined at higher loadings due to slurry formation, which reduces mechanical energy transfer and local heat.39 The optimized conditions of 30 min milling with η = 0.5 µL g−1 DMF were applied for subsequent screenings.
When ZrCp-BA-cluster was employed instead of ZrCp-OAc-cluster under the optimal conditions, a drastic drop in yield to 25% was observed (Fig. 2D). The surface area and FT-IR of the formed Zr-MOC are consistent with a tetrahedral cage (Fig. S10). We hypothesize two main reasons for the low yield: the release of solid benzoic acid upon milling, which modifies the LAG parameter, and the stronger intermolecular interactions between the ZrCp-BA-clusters, hindering ligand exchange with BDC, as indicated by single-crystal data (Fig. S3). Milling parameters' optimization is certainly possible with ZrCp-BA-cluster, but we focused on ZrCp-OAc-cluster for its higher atom economy. Nuclear magnetic resonance (NMR) and electron-spray ionization mass spectroscopy (ESI-MS) of the optimized Zr-MOC-BDC confirmed the successful preparation of the cage. The 1H NMR displays only 3 main peaks corresponding to the Cp*, ligand, and OH-bridging of Zr-clusters (Fig. S11). ESI-MS peaks at 781.5, 1043.1, and 1564.7 m/z are consistent with [Zr-MOC – 4 Cl−]4+, [Zr-MOC – 4 Cl− – H]3+, and [Zr-MOC – 4 Cl− – 2H]2+ m/z values for the tetrahedron phase (Fig. S12).33 With the optimized conditions, the mechanochemical synthesis of Zr-MOC-BDC is 16 times faster, and consumes around 7 times less energy and 115 times less solvent (excluding the cluster synthesis), than the conventional solvothermal protocol (Table S2).
To assess the versatility of our newly developed mechanochemical preparation method, we screened various bi-dentate and tri-dentate ligands known to form MOCs by solvothermal synthesis. For the bi-dentate case, 2-amino-terepthalic acid (H2BDC-NH2) and 2,6-napthalene dicarboxylic acid (H2ndc) were selected to examine the influence of functional groups and extended ligand length, respectively. The BM Zr-MOC-BDC-NH2 and Zr-MOC-ndc exhibit BET surface areas in the same range as previous reports (Fig. 3A and S19, S20).34,40 Notably, BM Zr-MOC-ndc shows the highest reported surface area (584 m2 g−1), confirming the phase purity achievable by mechanochemistry.13,14 The 1H NMR data of both cages are consistent with previous reports, with peaks from the Cp*, Zr-cluster, and ligand (Fig. S14 and S15). FT-IR of Zr-MOC-BDC-NH2 confirms the presence of primary amine moieties with a peak around 1335 cm−1 (Fig. S13).40 Moving from bi-dentate to tri-dentate ligands can be challenging due to the higher connectivity, which often leads to defect formation or phase diversity. The mechanochemically prepared cage with trimesic acid (H3BTC) displays a moderate surface area of 178 m2 g−1, approximately two times lower than that obtained from solvothermal synthesis (Fig. S18).33,41 When looking at the 1H NMR spectrum of the cage, we observed multiple peaks for the Cp* and BTC ligand (Fig. S14). The peak splitting of BTC was also observed in solvothermal synthesis and attributed to the presence or absence of DMF inside the cage.33,41 The pore aperture of Zr-MOC-BTC (4.54 Å) is smaller than the kinetic diameter of DMF (5.95 Å), blocking it from moving out. As mechanochemical synthesis is kinetically controlled in our case (short reaction time and low amount of solvent), the cage cannot rearrange to form guest-free products, which explains the reduced surface area. In the case of Cp*, structural defects could account for variations in the local environment, resulting in multiple signals. Unlike with BDC, milling without solvent didn't produce a porous product, highlighting the challenge of using tri-dentate ligands.
![]() | ||
| Fig. 3 Versatility of mechanochemical Zr-MOC synthesis. (A) Summary of ligands used and surface area. Reported Zr-cage surface areas: BDC-NH2,34,40 2,6-ndc,13,14 BTC,33,41 BTB.33 (B) Zr-MOC yields from mechanochemical preparation. | ||
Using the larger 1,3,5-tris(4-carboxyphenyl)benzene (H3BTB) gave a cage with a surface area close to 500 m2 g−1, on par with its solvothermal counterpart (Fig. S22).33 The wider pore width permits the proper activation of Zr-MOC-BTB. 1H NMR spectrum did not show peak splitting of BTB, confirming the formation of a guest-free cage (Fig. S18). However, the Cp* signal has a small peak shoulder, which can originate from the presence of a minor number of defects. The mechanochemical cage formation is promoted by using BTB instead of BTC, thanks to the size of BTB, which avoids steric hindrance between the Zr-clusters.
The tetrahedral geometry of the Zr-MOCs was further confirmed by ESI-MS analysis conducted in MeOH. The expected m/z peaks corresponding to +4, +3, and +2 tetrahedral ions were observed for Zr-MOC-BDC-NH2 and Zr-MOC-ndc (Fig. S23 and S24).13,34 For Zr-MOC-BTC, the main peaks were followed by smaller ones with m/z differences of approximately 32, 73, and 85, attributed to MeOH, DMF, and DCM molecules trapped within or bound to the cage, respectively (Fig. S25).33 This is in good agreement with the peak splitting observed by 1H NMR and explains the low surface area. ESI-MS of Zr-MOC-BTB was unsuccessful as it is not soluble in MeOH, but the 1H NMR and FT-IR data support the formation of the cage (Fig. S18 and S26). PXRD patterns of the as-synthesized BM cages show broad signals that partially match the simulated ones obtained from reported single-crystal data (Fig. S27–S30). The MOC crystals are stabilized by weak intermolecular forces among the cages, with solvent molecules playing a crucial role in crystal formation. The activated cages are largely amorphous and exhibit only broad signals, emphasizing their weak intermolecular interactions. The limited amount of DMF combined with the short reaction time likely prevents the development of long-range order, explaining the low crystallinity and deviation from single-crystal data. After washing, only broad signals remain.
As a final proof of concept for our mechanochemical approach, we synthesized a new tetrahedral cage with 2,5-dihydroxyterephthalic acid (BDC-(OH)2) and compared the data with the ball-milled sample. This ligand was selected to test the effect of two functional groups on the mechanochemical synthesis. Yellow cubic crystals were obtained after letting Cp2ZrCl2 and BDC-(OH)2 undisturbed in dimethylacetamide for two days, with a yield of around 30%. Single-crystal X-ray diffraction analysis revealed that the new MOC crystallized in the cubic space group Fm
m, similar to other Zr-MOC structures (Table S3).34 As expected, Zr-MOC-BDC-(OH)2 is a tetrahedral cage with four Cp3Zr3O(OH)3 extremities and six ligands (Fig. 4A). The charge is balanced by chloride anions.
Each cubic unit cell contains eight MOCs interconnected by hydrogen bonds between µ2-OH of the Zr-SBU and Cl− (Fig. 4B). The PXRD pattern of the product matched the simulated one (Fig. 4C). The solvothermal cage shows a BET surface area of 266 m2 g−1, which is lower than Zr-MOC-BDC and Zr-MOC-BDC-NH2 due to the smaller pore size induced by the hydroxyl groups. Preparing the cage by ball milling with the optimized conditions yields a yellow powder with a PXRD pattern that exhibits broad peaks (Fig. S32). The yield reaches 69% in only 30 min, more than double that of solvothermal synthesis, confirming the advantage of mechanochemistry. The porosity of the ball-milled Zr-MOC-BDC-(OH)2 is almost identical to the solvothermal sample, with a BET surface area of 256 m2 g−1 (Fig. 4D). The adsorption curve displays a non-well-defined plateau pressure attributed to the small particle size. 1H NMR signals and FT-IR also confirm the successful formation of the mechanochemical cage with signals from the ligand, cluster, and Cp* (Fig. S33 and S34). This demonstrates the functional group tolerance of the mechanochemical synthesis.
The data supporting this article have been included as part of the supplementary information (SI). Supplementary information: Fig. S1–S31, Tables S1, S2, and experimental details. See DOI: https://doi.org/10.1039/d6ta00338a.
| This journal is © The Royal Society of Chemistry 2026 |