Frontiers of molecular crystal structure prediction for pharmaceuticals and functional organic materials

The reliability of organic molecular crystal structure prediction has improved tremendously in recent years. Crystal structure predictions for small, mostly rigid molecules are quickly becoming routine. Structure predictions for larger, highly flexible molecules are more challenging, but their crystal structures can also now be predicted with increasing rates of success. These advances are ushering in a new era where crystal structure prediction drives the experimental discovery of new solid forms. After briefly discussing the computational methods that enable successful crystal structure prediction, this perspective presents case studies from the literature that demonstrate how state-of-the-art crystal structure prediction can transform how scientists approach problems involving the organic solid state. Applications to pharmaceuticals, porous organic materials, photomechanical crystals, organic semi-conductors, and nuclear magnetic resonance crystallography are included. Finally, efforts to improve our understanding of which predicted crystal structures can actually be produced experimentally and other outstanding challenges are discussed.

Polymorphs of Substituted p-Toluenesulfonanilide: Synthesis, Single-Crystal Analysis, Hirshfeld Surface Exploration, and Theoretical Investigation

Polymorphism is an exciting feature of chemical systems where a compound can exist in different crystal forms. The present investigation is focused on the two polymorphic forms, triclinic (MSBT) and monoclinic (MSBM), of ethyl 3-iodo-4-((4-methylphenyl)sulfonamido)benzoate prepared from ethyl 4-amino-3-iodobenzoate. The prepared polymorphs were unambiguously confirmed by single-crystal X-ray diffraction (SC-XRD) analysis. According to the SC-XRD results, the molecular configurations of both structures are stabilized by intramolecular N-H···I and C-H···O bonding. The crystal packing of MSBT is different as compared to the crystal packing of MSBM because MSBT is crystallized in the triclinic crystal system with the space group P1̅, whereas MSBM is crystallized in the monoclinic crystal system with the space group P21/c. The molecules of MSBT are interlinked in the form of dimers through N-H···O bonding to form R22(8) loops, while the MSBM molecules are connected with each other in the form of an infinite chain through C-H···O bonding. The crystal packing of both compounds is further stabilized by off-set π···π stacking interactions between phenyl rings, which is found stronger in MSBM as compared to in MSBT. Moreover, Hirshfeld surface exploration of the polymorphs was carried out, and the results were compared with the closely related literature structure. Accordingly, the supramolecular assembly of these polymorphs is mainly stabilized by noncovalent interactions or intermolecular interactions. Furthermore, a density functional theory (DFT) study was also carried out, which provided good support for the SC-XRD and Hirshfeld studies, suggesting the formation of both intramolecular and intermolecular interactions for both compounds.

Behavior of Occupied and Void Space in Molecular Crystal Structures at High Pressure

We report a Monte Carlo algorithm for calculation of occupied ("network") and unoccupied ("void") space in crystal structures. The variation of the volumes of the voids and the network of intermolecular contacts with pressure sensitively reveals discontinuities associated with first- and second-order phase transitions, providing insights into the effect of compression (and, in principle, other external stimuli) at a level between those observed in individual contact distances and the overall unit cell dimensions. The method is shown to be especially useful for the correlation of high-pressure crystallographic and spectroscopic data, illustrated for naphthalene, where a phase transition previously detected by vibrational spectroscopy, and debated in the literature for over 80 years, has been revealed unambiguously in crystallographic data for the first time. Premonitory behavior before a phase transition and crystal collapse at the end of a compression series has also been detected. The network and void volumes for 129 high-pressure studies taken from the Cambridge Structural Database (CSD) were fitted to equation of state to show that networks typically have bulk moduli between 40 and 150 GPa, while those of voids fall into a much smaller range, 2-5 GPa. These figures are shown to reproduce the narrow range of overall bulk moduli of molecular solids (ca. 5-20 GPa). The program, called CellVol, has been written in Python using the CSD Python API and can be run through the command line or through the Cambridge Crystallographic Data Centre's Mercury interface.

A first-order phase transition in Blatter's radical at high pressure

The crystal structure of Blatter's radical (1,3-diphenyl-1,4-dihydrobenzo[e][1,2,4]triazin-4-yl) has been investigated between ambient pressure and 6.07 GPa. The sample remains in a compressed form of the ambient-pressure phase up to 5.34 GPa, the largest direction of strain being parallel to the direction of π-stacking interactions. The bulk modulus is 7.4 (6) GPa, with a pressure derivative equal to 9.33 (11). As pressure increases, the phenyl groups attached to the N1 and C3 positions of the triazinyl moieties of neighbouring pairs of molecules approach each other, causing the former to begin to rotate between 3.42 to 5.34 GPa. The onset of this phenyl rotation may be interpreted as a second-order phase transition which introduces a new mode for accommodating pressure. It is premonitory to a first-order isosymmetric phase transition which occurs on increasing pressure from 5.34 to 5.54 GPa. Although the phase transition is driven by volume minimization, rather than relief of unfavourable contacts, it is accompanied by a sharp jump in the orientation of the rotation angle of the phenyl group. DFT calculations suggest that the adoption of a more planar conformation by the triazinyl moiety at the phase transition can be attributed to relief of intramolecular H...H contacts at the transition. Although no dimerization of the radicals occurs, the π-stacking interactions are compressed by 0.341 (3) Å between ambient pressure and 6.07 GPa.

How many more polymorphs of ROY remain undiscovered

With 12 crystal forms, 5-methyl-2-[(2-nitrophenyl)amino]-3-thiophenecabonitrile (a.k.a. ROY) holds the current record for the largest number of fully characterized organic crystal polymorphs. Four of these polymorph structures have been reported since 2019, raising the question of how many more ROY polymorphs await future discovery. Employing crystal structure prediction and accurate energy rankings derived from conformational energy-corrected density functional theory, this study presents the first crystal energy landscape for ROY that agrees well with experiment. The lattice energies suggest that the seven most stable ROY polymorphs (and nine of the twelve lowest-energy forms) on the Z' = 1 landscape have already been discovered experimentally. Discovering any new polymorphs at ambient pressure will likely require specialized crystallization techniques capable of trapping metastable forms. At pressures above 10 GPa, however, a new crystal form is predicted to become enthalpically more stable than all known polymorphs, suggesting that further high-pressure experiments on ROY may be warranted. This work highlights the value of high-accuracy crystal structure prediction for solid-form screening and demonstrates how pragmatic conformational energy corrections can overcome the limitations of conventional density functionals for conformational polymorphs.

Use of a miniature diamond-anvil cell in a joint X-ray and neutron high-pressure study on copper sulfate pentahydrate

Single-crystal X-ray and neutron diffraction data are usually collected using separate samples. This is a disadvantage when the sample is studied at high pressure because it is very difficult to achieve exactly the same pressure in two separate experiments, especially if the neutron data are collected using Laue methods where precise absolute values of the unit-cell dimensions cannot be measured to check how close the pressures are. In this study, diffraction data have been collected under the same conditions on the same sample of copper(II) sulfate pentahydrate, using a conventional laboratory diffractometer and source for the X-ray measurements and the Koala single-crystal Laue diffractometer at the ANSTO facility for the neutron measurements. The sample, of dimensions 0.40 × 0.22 × 0.20 mm3 and held at a pressure of 0.71 GPa, was contained in a miniature Merrill-Bassett diamond-anvil cell. The highly penetrating diffracted neutron beams passing through the metal body of the miniature cell as well as through the diamonds yielded data suitable for structure refinement, and compensated for the low completeness of the X-ray measurements, which was only 24% on account of the triclinic symmetry of the sample and the shading of reciprocal space by the cell. The two data-sets were combined in a single 'XN' structure refinement in which all atoms, including H atoms, were refined with anisotropic displacement parameters. The precision of the structural parameters was improved by a factor of up to 50% in the XN refinement compared with refinements using the X-ray or neutron data separately.

Direct evidence for distinct colour origins in ROY polymorphs

ROY is one of the most well-studied families of crystal structures owing to it being the most polymorphic organic material on record. The various red, orange, and yellow colours of its crystal structures are widely-believed to originate from molecular conformation, though the orange needle (ON) polymorph is thought to be an exception. We report high-pressure, single-crystal X-ray measurements which provide direct experimental evidence that the colour origin in ON is intermolecular, revealing that the molecule undergoes minimal deformation but still exhibits a pronounced, reversible, pale orange → dark red colour change between ambient pressure and 4.18 GPa. Our experimental data are rationalised with band structures, calculated using an accurate hybrid DFT approach, where we are able to account for the variation in colour for five polymorphs of ROY. We highlight the outlier behaviour of ON which shows marked π⋯π stacking interactions that are directly modified through application of pressure. Band structure calculations confirm these intermolecular interactions as the origin of the colour change.

Prolific Polymorph Generator ROY in Its Liquid and Glass: Two Conformational Populations Mirroring the Crystalline-State Distribution

5-Methyl-2-[(2-nitrophenyl)amino]-3-thiophenecarbonitrile, dubbed ROY for its numerous crystal polymorphs of red, orange, and yellow colors, has been studied in its liquid and glassy state by infrared spectroscopy. Two populations of conformers are observed, whose equilibrium is characterized by ΔH = 2.4 kJ/mol and ΔS = 8.0 J/K/mol. The two populations correspond to the global and local minima of the torsional energy surface and to the conformational preference of the 13 crystal polymorphs. The local minimum features a more coplanar arrangement of the two aromatic rings, greater π conjugation, and lower CN stretch frequency. In the gas phase, the lowest-energy path between the two minima has an energy barrier 3.9 kJ/mol above the global minimum, consistent with the rapid equilibration between the two populations. The relevance of our result for understanding the prolific polymorphism of ROY is discussed.

Probing the Influence of Defects, Hydration, and Composition on Prussian Blue Analogues with Pressure

The vast compositional space of Prussian blue analogues (PBAs), formula A_xM[M′(CN)₆]_y·nH₂O, allows for a diverse range of functionality. Yet, the interplay between composition and physical properties—e.g., flexibility and propensity for phase transitions—is still largely unknown, despite its fundamental and industrial relevance. Here we use variable-pressure X-ray and neutron diffraction to explore how key structural features, i.e., defects, hydration, and composition, influence the compressibility and phase behavior of PBAs. Defects enhance the flexibility, manifesting as a remarkably low bulk modulus (B₀ ≈ 6 GPa) for defective PBAs. Interstitial water increases B₀ and enables a pressure-induced phase transition in defective systems. Conversely, hydration does not alter the compressibility of stoichiometric MnPt(CN)₆, but changes the high-pressure phase transitions, suggesting an interplay between low-energy distortions. AMnCo(CN)₆ (A^I = Rb, Cs) transition from F4̅3m to P4̅n2 upon compression due to octahedral tilting, and the critical pressure can be tuned by the A-site cation. At 1 GPa, the symmetry of Rb_0.87Mn[Co(CN)₆]_0.91 is further lowered to the polar space group Pn by an improper ferroelectric mechanism. These fundamental insights aim to facilitate the rational design of PBAs for applications within a wide range of fields.

MrPIXEL: automated execution of Pixel calculations via the Mercury interface

A program is described that enables the energetics of crystal packing to be analysed quickly and easily through a user interface, using Pixel, Python scripts and the widely known visualizer Mercury.

The interpretation of crystal structures in terms of intermolecular interaction energies enables phase stability and polymorphism to be rationalized in terms of quantitative thermodynamic models, while also providing insight into the origin of physical and chemical properties including solubility, compressibility and host–guest formation. The Pixel method is a semi-empirical procedure for the calculation of intermolecular interactions and lattice energies based only on crystal structure information. Molecules are represented as blocks of undistorted ab initio molecular electron and nuclear densities subdivided into small volume elements called pixels. Electrostatic, polarization, dispersion and Pauli repulsion terms are calculated between pairs of pixels and nuclei in different molecules, with the accumulated sum equating to the intermolecular interaction energy, which is broken down into physically meaningful component terms. The MrPIXEL procedure enables Pixel calculations to be carried out with minimal user intervention from the graphical interface of Mercury, which is part of the software distributed with the Cambridge Structural Database (CSD). Following initial setup of a crystallographic model, one module assigns atom types and writes necessary input files. A second module then submits the required electron-density calculation either locally or to a remote server, downloads the results, and submits the Pixel calculation itself. Full lattice energy calculations can be performed for structures with up to two molecules in the crystallographic asymmetric unit. For more complex cases, only molecule–molecule energies are calculated. The program makes use of the CSD Python API, which is also distributed with the CSD.

ROY Reclaims Its Crown: New Ways To Increase Polymorphic Diversity

Chemical compounds that exist in multiple crystalline forms are said to exhibit polymorphism. Polymorphs have the same composition, but their structures and properties can vary markedly. In many fields, conditions for crystallizing compounds of interest are screened exhaustively to generate as many polymorphs as possible, from which the most advantageous form can be selected. We report new ways to search for polymorphs and increase polymorphic diversity, based on crystallization induced by suitably designed mixed-crystal seeds. The potential of the strategy has been demonstrated by using it to produce new polymorphs of the benchmark compound ROY as single crystals structurally characterized by X-ray diffraction. This allows ROY to reclaim its crown as the most polymorphic compound in the Cambridge Structural Database. More generally, the methods promise to become valuable tools for polymorphic screening in all fields where crystalline solids are used.

High-pressure polymorphism in pyridine

Single crystals of the high-pressure phases II and III of pyridine have been obtained by in situ crystallization at 1.09 and 1.69 GPa, revealing the crystal structure of phase III for the first time using X-ray diffraction. Phase II crystallizes in P212121 with Z' = 1 and phase III in P41212 with Z' = ½. Neutron powder diffraction experiments using pyridine-d5 establish approximate equations of state of both phases. The space group and unit-cell dimensions of phase III are similar to the structures of other simple compounds with C 2v molecular symmetry, and the phase becomes stable at high pressure because it is topologically close-packed, resulting in a lower molar volume than the topologically body-centred cubic phase II. Phases II and III have been observed previously by Raman spectroscopy, but have been mis-identified or inconsistently named. Raman spectra collected on the same samples as used in the X-ray experiments establish the vibrational characteristics of both phases unambiguously. The pyridine molecules interact in both phases through CH⋯π and CH⋯N interactions. The nature of individual contacts is preserved through the phase transition between phases III and II, which occurs on decompression. A combination of rigid-body symmetry mode analysis and density functional theory calculations enables the soft vibrational lattice mode which governs the transformation to be identified.

Pressure-induced inclusion of neon in the crystal structure of a molecular Cu2(pacman) complex at 4.67 GPa

Crystals of Cu₂(pacman) inflate on taking up neon at 46 000 atm through a switch in the ligand conformation.

The Effect of Pressure on Halogen Bonding in 4-Iodobenzonitrile

The crystal structure of 4-iodobenzonitrile, which is monoclinic (space group I2/a) under ambient conditions, contains chains of molecules linked through C≡N···I halogen-bonds. The chains interact through CH···I, CH···N and π-stacking contacts. The crystal structure remains in the same phase up to 5.0 GPa, the b axis compressing by 3.3%, and the a and c axes by 12.3 and 10.9 %. Since the chains are exactly aligned with the crystallographic b axis these data characterise the compressibility of the I···N interaction relative to the inter-chain interactions, and indicate that the halogen bond is the most robust intermolecular interaction in the structure, shortening from 3.168(4) at ambient pressure to 2.840(1) Å at 5.0 GPa. The π∙∙∙π contacts are most sensitive to pressure, and in one case the perpendicular stacking distance shortens from 3.6420(8) to 3.139(4) Å. Packing energy calculations (PIXEL) indicate that the π∙∙∙π interactions have been distorted into a destabilising region of their potentials at 5.0 GPa. The structure undergoes a transition to a triclinic ( P 1 ¯ ) phase at 5.5 GPa. Over the course of the transition, the initially colourless and transparent crystal darkens on account of formation of microscopic cracks. The resistance drops by 10% and the optical transmittance drops by almost two orders of magnitude. The I···N bond increases in length to 2.928(10) Å and become less linear [<C-I∙∙∙N = 166.2(5)°]; the energy stabilises by 2.5 kJ mol^-1 and the mixed C-I/I..N stretching frequency observed by Raman spectroscopy increases from 249 to 252 cm^-1. The driving force of the transition is shown to be relief of strain built-up in the π∙∙∙π interactions rather than minimisation of the molar volume. The triclinic phase persists up to 8.1 GPa.