Abstract

Crystal structure prediction methods can enable the in silico design of functional molecular crystals, but solvent effects can have a major influence on relative lattice energies, sometimes thwarting predictions. This is particularly true for porous solids, where solvent included in the pores can have an important energetic contribution. We present a Monte Carlo solvent insertion procedure for predicting the solvent filling of porous structures from crystal structure prediction landscapes, tested using a highly solvatomorphic porous organic cage molecule, CC1. Using this method, we can understand why the predicted global energy minimum structure for CC1 is never observed from solvent crystallisation. We also explain the formation of three different solvatomorphs of CC1 from three structurally-similar chlorinated solvents. Calculated solvent stabilisation energies are found to correlate with experimental results from thermogravimetric analysis, suggesting a future computational framework for a priori materials design that factors in solvation effects.

For reproduction of material from NJC:
[Original citation] - Published by The Royal Society of Chemistry (RSC) on behalf of the
Centre National de la Recherche Scientifique (CNRS) and the RSC.

For reproduction of material from PCCP:
[Original citation] - Published by the PCCP Owner Societies.

For reproduction of material from PPS:
[Original citation] - Published by The Royal Society of Chemistry (RSC) on behalf of the
European Society for Photobiology, the European Photochemistry Association, and
RSC.

For reproduction of material from all other RSC journals:
[Original citation] - Published by The Royal Society of Chemistry.

Information about reproducing material from RSC articles with different licences
is available on our Permission Requests page.