Interactions of Common Synthesis Solvents with MOFs Studied via Free Energies of Solvation: Implications on Stability and Polymorph Selection

Solvents are central to metal–organic framework (MOF) solvothermal synthesis. However, how solvent–MOF interplay impacts MOF stabilization individually and across polymorphs is not well understood. To address this knowledge gap, here we perform data-driven analysis of 20,532 heats of adsorption at dilute conditions (Δ H ads ) and 447 free energies of solvation (Δ F sol ) for four solvents, namely, dimethylformamide (DMF), water (H 2 O), methanol (MeOH), and n -hexane (C6) (which was used as a control). To accelerate data collection, we developed a protocol to extrapolate Δ F sol from calculations with the MOFs only partially solvated. Free energies were obtained via thermodynamic integration. We found Δ F sol and Δ H ads to be only moderately correlated due to solvent–solvent interactions coming into play when the MOF is solvated. In any case, trends in Δ F sol were ultimately explained on the basis of solvent kinetic diameter and polarity, as well as MOF void fraction (V f ) and functionalization polarity. For instance, the correlation between Δ F sol and V f was one of the strongest correlations presented in this study (more so as the solvent size increases), indicating that small-pore MOFs are more easily stabilized by solvation than large-pore MOFs. We also found that solvation-induced MOF stabilization became more pronounced as solvent kinetic diameter (polarity) decreased (increased). We found differences in this solvation-induced stabilization between polymorphs capable of overcoming inherent (i.e., in vacuum) differences in polymorph stability, causing the most stable polymorph to “switch.” We found the probability to cause “switches” to increase as solvent kinetic diameter (polarity) decreased (increased). Inspection of multivariate linear regression coefficients suggested that differences in solvation-induced stabilization in polymorphs can be primarily explained by their differences in density, V f, and, to a lesser extent, volumetric surface area.