RESUMO
The physicochemical properties of molecular crystals, such as solubility, stability, compactability, melting behaviour and bioavailability, depend on their crystal form1. In silico crystal form selection has recently come much closer to realization because of the development of accurate and affordable free-energy calculations2-4. Here we redefine the state of the art, primarily by improving the accuracy of free-energy calculations, constructing a reliable experimental benchmark for solid-solid free-energy differences, quantifying statistical errors for the computed free energies and placing both hydrate crystal structures of different stoichiometries and anhydrate crystal structures on the same energy landscape, with defined error bars, as a function of temperature and relative humidity. The calculated free energies have standard errors of 1-2 kJ mol-1 for industrially relevant compounds, and the method to place crystal structures with different hydrate stoichiometries on the same energy landscape can be extended to other multi-component systems, including solvates. These contributions reduce the gap between the needs of the experimentalist and the capabilities of modern computational tools, transforming crystal structure prediction into a more reliable and actionable procedure that can be used in combination with experimental evidence to direct crystal form selection and establish control5.
RESUMO
We exploit the possible link between structural surface roughness and difficulty of crystallisation. Polymorphs with smooth surfaces may nucleate and crystallise more readily than polymorphs with rough surfaces. The concept is applied to crystal structure prediction landscapes and reveals a promising complementary way of ranking putative crystal structures.
RESUMO
Loratadine, an over-the-counter antihistamine medication, has two known monotropically related polymorphs, both of which feature disorder. A combined experimental and computational approach using variable temperature single crystal X-ray diffraction (VT-SCXRD) analysis and dispersion corrected density functional theory (DFT-D) reveals that the nature of the disorder in each form is markedly different and cannot be described by a simple isolated-site model with thermally populated conformations in either of the two cases. In Form I, the ethyl carbamate functionality adopts two different configurations, with adjacent moieties interacting along one-dimensional chains. The most stable arrangement features alternating configurations, but because of the low energetic cost of stacking faults, the domain sizes are short and an average crystal structure is observed experimentally. The configurational free energy of the disordered structure is lower than the energy of the two corresponding ordered crystal structures, but the energy decrease is dominated by the lower lattice energy of the alternating arrangement with a small entropic contribution. In Form II, the flexible cycloheptane bridge adopts two different configurations. The disorder is not an equilibrium property but is instead frozen-in during the crystallisation process. The configurational free energy of the disordered structure falls in between the lattice energies of the two corresponding ordered structures. The two ordered components of each disordered structure are all found in a crystal structure prediction (CSP) study with the GRACE programme. However, the experimentally observed stability relationship is only reproduced when the energy contribution of disorder is taken into account. The disordered model of Form I is found to be lower in energy than all the other predicted structures and there is no indication of a missing, thermodynamically more stable, form. The case of loratadine demonstrates that experimentally observed disorder close to 50/50 does not necessarily correspond to a free energy decrease by kT ln 2.