Open Access Article
Michelle
Ernst
*a and
Tomasz
Poręba
*b
aInstitute of Geological Sciences, University of Bern, 3012 Bern, Switzerland. E-mail: michelle.ernst@unibe.ch
bLaboratory for Quantum Magnetism, Institute of Physics, École Polytechnique Fédérale de Lausanne, Lausanne CH-1015, Switzerland. E-mail: tomasz.poreba@epfl.ch
First published on 2nd October 2025
Cadmium–porphyrin metal–organic frameworks combine the structural rigidity of the porphyrin core with the coordination flexibility of closed-shell Cd2+ centres, enabling diverse framework topologies and versatile applications. Here, we report a single-crystal-to-single-crystal transformation in a CdCl2–meso-tetra(4-pyridyl)porphyrin (H2TPyP) MOF induced by mild thermal desolvation (120 °C in air). The solvated MOF (UB-MOF-1) adopts a non-interpenetrated CdSO4-type topology stabilized by 1,1,2,2-tetrachloroethane and 3-chloroaniline guests. Upon solvent removal, the framework undergoes irreversible drastic rearrangement to a two-fold interpenetrated cds-type network (ESRF-MOF-1), accompanied by more than a twofold reduction in accessible pore volume. Periodic DFT calculations show that, while guest inclusion substantially stabilizes the open framework UB-MOF-1, the interpenetrated, solvent-free ESRF-MOF-1 is thermodynamically favoured over the hypothetical empty UB-MOF-1 as interpenetration provides even greater stabilization. The transformation, driven by framework densification, is irreversible due to geometric pore constraints. This study demonstrates how controlled solvent manipulation can trigger topology changes in porphyrin-based MOFs, offering a route to design dense, structurally locked architectures with tailored accessibility.
Metal–organic frameworks (MOFs) represent a broader class of coordination polymers characterized by their permanent porosity and crystalline order. They have received ongoing attention due to their structural richness and potential for tuneable porosity and functionality. In this study, we explore solvent-assisted structural changes in two related Cd–porphyrin-based MOFs.
Previous work has highlighted the structural adaptability and solvent-responsive behaviour of Cd–TPyP MOFs. Subtle changes, such as metal type, synthetic conditions, recrystallization method, and guest inclusion can significantly influence the dimensionality, degree of interpenetration, and porosity of these porphyrinic coordination networks.6 Several 3D Cd–TPyP MOFs have been reported, featuring interpenetrated networks and remarkable coordination diversity: individual porphyrin ligands were observed to bind one, two, four, or five (through porphyrin metalation) cadmium centres, while Cd2+ nodes exhibited five distinct coordination geometries with various supramolecular architectures.6,7 For example, combining freebase or partially metalated TPyP with CdI2 subjected to a solvent-diffusion crystallization process yields layered, solvated coordination polymers that accommodated multiple tetrachloroethane molecules between porphyrin layers.8 This suggests that solvent inclusion during crystallization contributes to the thermodynamic stabilization of these layered structures in the solid state. In contrast, HMOF-1, a non-interpenetrated, solvent-free CdI2–TPyP framework (crystallized under similar conditions) displays a 3D hinged “lattice-fence” topology with extraordinary anisotropic reversible thermal expansion.5 This phenomenon was attributed to the large hinge-like motion range around the Cd centres within the framework.†
By contrast, despite having a similar van der Waals ligand radius, Cd(SH)2(H2TPyP) crystallized under hydrothermal conditions forms an interpenetrated, solvent-free, cds-type MOF and not flexible hinged structure.7 These examples highlight the remarkable structural diversity of Cd–TPyP MOFs, even when built from the same ligands, as a result of different synthetic and crystallization pathways. In this work, we present an alternative post-synthetic route to modify the framework packing: thermal activation removes the solvent molecules and induces an increase in the degree of interpenetration.
Changes in the degree of interpenetration upon solvent removal have been widely observed across a broad range of MOF systems.9 Transformations from non-interpenetrated to interpenetrated structures have been reported in multiple studies, demonstrating the role of solvent molecules in stabilizing non-interpenetrated frameworks and the tendency of the framework to collapse into a denser interpenetrated structure once these molecules are removed.10–15
In many cases, increases in the degree of interpenetration have been observed, such as two-fold to three-fold16,17 and two-fold to four-fold,18 as well as more complex transitions including four-fold to five-fold,19 five-fold to six-fold,20 and even seven-fold to eight-fold.21 This highlights the general tendency of interpenetration to increase as solvent molecules are removed and structural relaxation occurs. In contrast, transformations to a lower degree of interpenetration are exceedingly rare, and, we are not aware of any examples of transformations that decrease the degree of interpenetration upon simple solvent removal; instead, reported cases involve reactions leading to deinterpenetration,22,23 interpenetration increases induced by guest addition,24 or even reversible interpenetration changes triggered by solvent exchange.25
In this context, understanding why and how guest removal can trigger such radical topological rearrangements is essential for the rational design of adaptive and responsive materials. The ability to rationalize and predict interpenetration transitions remains a central goal in crystal engineering, with direct implications for designing MOFs that combine functionality with stability, porosity, and dynamic behaviour in targeted ways.
| Empirical formula | C40H26CdCl2N8 |
|---|---|
| Formula weight | 801.99 |
| Temperature/K | 298(1) |
| Crystal system | Orthorhombic |
| Space group | Pnna |
| a/Å | 12.8121(2) |
| b/Å | 15.6630(2) |
| c/Å | 20.3134(4) |
| α/° | 90 |
| β/° | 90 |
| γ/° | 90 |
| Volume/Å3 | 4076.41(11) |
| Z | 4 |
| ρ calc g cm−3 | 1.307 |
| μ/mm−1 | 0.729 |
| F(000) | 1616.0 |
| Crystal size/mm3 | 0.1 × 0.05 × 0.04 |
| Radiation | Synchrotron (λ = 0.40988) |
| 2Θ range for data collection/° | 3.31 to 29.686 |
| Index ranges | −15 ≤ h ≤ 16, −19 ≤ k ≤ 19, −21 ≤ l ≤ 25 |
| Reflections collected | 15 493 |
| Independent reflections | 4070 [Rint = 0.0463, Rsigma = 0.0680] |
| Data/restraints/parameters | 4070/0/234 |
| Goodness-of-fit on F2 | 1.129 |
| Final R indexes [I ≧ 2σ(I)] | R 1 = 0.0675, wR2 = 0.2276 |
| Final R indexes [all data] | R 1 = 0.0739, wR2 = 0.2327 |
| Largest diff. peak/hole/e Å−3 | 1.31/−0.54 |
The binding energy Ebind of UB-MOF-1 corresponds to the sum of the electronic interaction energy of the framework with the solvent molecules plus the reorganization energy (cf.eqn (1)).
| Ebind = Efilled − (Eempty + Esolvent) | (1) |
:
1 TCE–MeOH solution of H24TPyP layered with metal salts dissolved in MeOH. The layers of the resulting MOF are intercalated with TCE molecules in a 1
:
2 molar ratio. A similar synthesis was reported by DeVries in which H24TPyP was likewise dissolved in 3
:
1 TCE–MeOH, layered with 3-chloroaniline, and finally topped with a N,N′-dimethylacetamide solution of CdI2. In this case, a cds-type solvent-free MOF was obtained. Although we carefully followed the synthetic procedure of DeVries5 that yielded the cds-type solvent-free MOF (CCDC 854909), we obtained a solvated MOF (UB-MOF-1) containing both TCE and 3-CA. Moreover, the CdI2 clusters in the framework were replaced by CdCl2 ones. The non-incorporation of iodine in this structure has been previously attributed to competing Finkelstein-type halide exchange reaction with excess of the TCE, as discussed in our earlier work.26
Heating of UB-MOF-1 in air at 120 °C for 4 hours resulted in a fully emptied MOF (ESRF-MOF-1). Remarkably, the mild conditions for the desorption process resulted in single-crystal to single-crystal (SC–SC) transition, which allowed for subsequent SC-XRD structural determination. The integrity of the crystal is retained despite the drastic structural changes upon solvent removal (Fig. 1 and S2). Both structures, UB-MOF-1 and ESRF-MOF-1, are built of CdCl2trans-octahedral clusters bound together by tetradentate meso-tetra(4-pyridyl)porphyrin (H24TPyP) ligands.
However, topological analysis revealed significant differences between the two extended structures. In both cases, the Cd nodes and porphyrin cores function as square-planar four-connected units. The solvated UB-MOF-1 is an example of a CdSO4-type net, where guest 1,1,2,2-tetrachloroethane is sandwiched between the porphyrin rings, and 3-chloroaniline occupies the interstitial voids (Fig. 1a and b). On the other hand, the emptied MOF, ESRF-MOF-1, forms a much more compact structure made of two interpenetrating identical cds nets (Fig. 1c and d).
The nodes of the net are shifted by 1/2a, 1/2b and 0.150c with respect to each other. As a result of the tight interpenetration, the structure of ESRF-MOF-1 is much denser than the theoretical solvent-free UB-MOF-1 (1.307 and 0.965 g cm−3, respectively) and shows nearly half of the network-accessible helium volume (29.8% versus 55.1%).
The closed-shell electronic structure of Cd allows for certain structural flexibility around the metal centres, which can be noted as a “wavy” pattern (Fig. 1a) formed by the cadmium clusters bonded to the pyridin-4-yl ligand at 9.3(2)°. This feature combined with the rigidity of the porphyrin linkers likely allows for the controlled release of both solvents (TCE and 3-CA) before the structure collapse into the interpenetrated net. Noteworthy, at the applied annealing conditions (120 °C, 1 atm) both solvents reach relatively high vapour pressure, which fosters their full desorption (∼370 mmHg and ∼25 mmHg for TCE and 3-CA respectively, estimated from Antoine coefficients).33,34
On the other hand, even though the maximum pore diameter differs only slightly between the two structures (4.68 Å in UB-MOF-1 and 4.37 Å in ESRF-MOF-1), the latter could not be reversed back to the solvated form. Attempts to restore the solvated phase by immersing ESRF-MOF-1 in a 1% w/w 3-chloroaniline solution in MeOH for 3 hours at room temperature were unsuccessful. Similarly, incorporation of EtOH, MeOH and H2O by the same technique were unsuccessful, as evidenced by the subsequently collected XRD data (Table S2). The reason for this behaviour likely lies in the small pore limiting diameter, the narrowest constriction along the pore channel, of ESRF-MOF-1 of 1.34 Å (3.1 Å in UB-MOF-1), which does not allow for molecules larger than He to penetrate the crystal interior. Analysis of the residual electron density in the structure, performed using the PLATON SQUEEZE module,35 indicates 100 residual electrons in the voids per unit cell. Since there are four voids, each contains 25 electrons, which corresponds to ∼0.4 3-chloroaniline molecules per void. Modelling of this residual density does not reveal any connected network that would resemble (partially) ordered 3-CA or TCE molecules. Nonetheless, we cannot exclude the possibility that the MOF is not fully desolvated, with some remaining guests disordered beyond the resolution limit. The residual density is likely not due to TCE, given its much lower evaporation point, as discussed above.
For comparison, cluster-based interaction energies between UB-MOF-1 and the guest solvent molecules yielded values of −114.5 kJ mol−1 for TCE and −124.6 kJ mol−1 for 3-CA.26 This corresponds to an interaction energy of −727.4 kJ mol−1 per unit cell (CdCl2(4TPyP)2·C2H2Cl4·2C6H4NH2Cl, Z′ = 0.5, Z = 2). The close agreement between the cluster-based values and our periodic calculations is reassuring.
Our periodic DFT calculations underline the thermodynamic driving force of this transformation. The interpenetrated ESRF-MOF-1 framework is calculated to be 112.8 kJ mol−1 per unit cell more stable than the hypothetical empty non-interpenetrated UB-MOF-1, consistent with the experimentally observed transformation. While the solvated UB-MOF-1 is stabilized by a strong interaction energy between the framework and the guest solvent molecules (765.6 kJ mol−1 per unit cell), this interaction is lost upon desolvation, prompting the framework to irreversibly reorganize into a more energetically favourable and compact topology.
It is difficult to generally estimate the influence of the framework–solvent interaction strength on the formation of interpenetrated networks due to different ways XRD structural data is treated before publishing.42 CSD v.5.44 contains 108
881 entries which can be classified as metal–organic frameworks and/or coordination polymers (search criteria: contains transition metal, contains “catena” keyword, 3D structure determined), whereas in nearly 11% of them a disordered part of a molecule and/or solvent have been removed using the PLATON SQUEEZE routine35 (11
971 entries, keyword “SQUEEZE”). In some cases, by omitting the solvent molecules one may obscure the underlying activation phenomena. For example, in the isostructural iodine analogue of ESRF-MOF-1, HMOF-1,5 the bending around the Cd node has been identified to be responsible for the highly anisotropic thermal expansion (without increase of the degree of interpenetration). However, the reported results of the elemental analysis (in the SI of the cited work) indicate the presence of impurities and/or residual solvent, which has been eliminated from the reported crystal structure using the SQUEEZE routine. In this case, the source of such high structural mobility without the increase of degree of interpenetration might lie in the phase transition associated with the solvent removal. In fact, the collected X-ray powder diffractograms for this compound show splitting of certain Bragg reflections above 300 K.
Attempts to reverse the transformation by re-soaking the ESRF-MOF-1 crystals in solutions of 3-chloroaniline, ethanol, methanol (also at 50 °C, see Fig. S3), and water failed, highlighting the irreversibility of the transition under mild conditions (Table S2). Once the interpenetrated structure forms, its narrow channels prevent diffusion of new guest molecules, thus trapping the framework in a dense, closed-pore state. In this sense, the transformation is not simply temperature-dependent, but accessibility-limited, and ESRF-MOF-1 behaves more like an “exhale” MOF than a breathing one.37
The structural data supporting the findings of this study are available in the Cambridge Crystallographic Data Centre (CCDC): CCDC 2480007 contains the supplementary crystallographic data for this paper.43
All computation files are openly available in the Materials Cloud Archive under https://doi.org/10.24435/materialscloud:ja-j4.
Footnote |
| † Even though, the unit cell dimensions are very close to the ones found for the MOF reported in this study (CdCl2–TPyP), the reported structure crystallizes in an orthorhombic space group with β = 90° and a “hinge angle” between the vector connecting the axial ligands of Cd and a plane of TPyP around 26.7°, while our structure adopts a monoclinic cell with β = 93.5°, with a hinge angle of 43.5°. |
| This journal is © The Royal Society of Chemistry 2025 |