The seventh blind test of crystal structure prediction: structure ranking methods

A seventh blind test of crystal structure prediction has been organized by the Cambridge Crystallographic Data Centre. The results are presented in two parts, with this second part focusing on methods for ranking crystal structures in order of stability. The exercise involved standardized sets of structures seeded from a range of structure generation methods. Participants from 22 groups applied several periodic DFT-D methods, machine learned potentials, force fields derived from empirical data or quantum chemical calculations, and various combinations of the above. In addition, one non-energy-based scoring function was used. Results showed that periodic DFT-D methods overall agreed with experimental data within expected error margins, while one machine learned model, applying system-specific AIMnet potentials, agreed with experiment in many cases demonstrating promise as an efficient alternative to DFT-based methods. For target XXXII, a consensus was reached across periodic DFT methods, with consistently high predicted energies of experimental forms relative to the global minimum (above 4 kJ mol-1 at both low and ambient temperatures) suggesting a more stable polymorph is likely not yet observed. The calculation of free energies at ambient temperatures offered improvement of predictions only in some cases (for targets XXVII and XXXI). Several avenues for future research have been suggested, highlighting the need for greater efficiency considering the vast amounts of resources utilized in many cases.

Assessment of the exchange-hole dipole moment dispersion correction for the energy ranking stage of the seventh crystal structure prediction blind test

The seventh blind test of crystal structure prediction (CSP) methods substantially increased the level of complexity of the target compounds relative to the previous tests organized by the Cambridge Crystallographic Data Centre. In this work, the performance of density-functional methods is assessed using numerical atomic orbitals and the exchange-hole dipole moment dispersion correction (XDM) for the energy-ranking phase of the seventh blind test. Overall, excellent performance was seen for the two rigid molecules (XXVII, XXVIII) and for the organic salt (XXXIII). However, for the agrochemical (XXXI) and pharmaceutical (XXXII) targets, the experimental polymorphs were ranked fairly high in energy amongst the provided candidate structures and inclusion of thermal free-energy corrections from the lattice vibrations was found to be essential for compound XXXI. Based on these results, it is proposed that the importance of vibrational free-energy corrections increases with the number of rotatable bonds.

Unprecedented Packing Polymorphism of Oxindole: An Exploration Inspired by Crystal Structure Prediction

Crystal polymorphism, characterized by different packing arrangements of the same compound, strongly ties to the physical properties of a molecule. Determining the polymorphic landscape is complex and time-consuming, with the number of experimentally observed polymorphs varying widely from molecule to molecule. Furthermore, disappearing polymorphs, the phenomenon whereby experimentally observed forms cannot be reproduced, pose a significant challenge for the pharmaceutical industry. Herein, we focused on oxindole (OX), a small rigid molecule with four known polymorphs, including a reported disappearing form. Using crystal structure prediction (CSP), we assessed OX solid-state landscape and thermodynamic stability by comparing predicted structures with experimentally known forms. We then performed melt and solution crystallization in bulk and nanoconfinement to validate our predictions. These experiments successfully reproduced the known forms and led to the discovery of four novel polymorphs. Our approach provided insights into reconstructing disappearing polymorphs and building more comprehensive polymorph landscapes. These results also establish a new record of packing polymorphism for rigid molecules.

Contributions beyond direct random-phase approximation in the binding energy of solid ethane, ethylene, and acetylene

The random-phase approximation (RPA) includes a subset of higher than second-order correlation-energy contributions, but stays in the same complexity class as the second-order Møller-Plesset perturbation theory (MP2) in both Gaussian-orbital and plane-wave codes. This makes RPA a promising ab initio electronic structure approach for the binding energies of molecular crystals. Still, some issues stand out in practical applications of RPA. Notably, compact clusters of nonpolar molecules are poorly described, and the interaction energies strongly depend on the reference single-determinant state. Using the many-body expansion of the binding energy of a crystal, we investigate those issues and the effect of beyond-RPA corrections. We find the beneficial effect of quartic-scaling exchange and non-ring coupled-cluster doubles corrections. The nonadditive interactions in compact trimers of molecules are improved by using the self-consistent Hartree-Fock orbitals instead of the usual Kohn-Sham states, but this kind of orbital input also leads to underestimated dimer energies. Overall, a substantial improvement over the RPA with a renormalized singles approach is possible at a modest quartic-scaling cost, which encourages further research into additional RPA corrections.

Linking solid-state phenomena via energy differences in `archetype crystal structures'

Categorization underlies understanding. Conceptualizing solid-state structures of organic molecules with `archetype crystal structures' bridges established categories of disorder, polymorphism and solid solutions and is herein extended to special position and high-Z' structures. The concept was developed in the context of disorder modelling [Dittrich, B. (2021). IUCrJ, 8, 305-318] and relies on adding quantum chemical energy differences between disorder components to other criteria as an explanation as to why disorder - and disappearing disorder - occurs in an average structure. Part of the concept is that disorder, as probed by diffraction, affects entire molecules, rather than just the parts of a molecule with differing conformations, and the finding that an R·T energy difference between disorder archetypes is usually not exceeded. An illustrative example combining disorder and special positions is the crystal structure of oestradiol hemihydrate analysed here, where its space-group/subgroup relationship is required to explain its disorder of hydrogen-bonded hydrogen atoms. In addition, we show how high-Z' structures can also be analysed energetically and understood via archetypes: high-Z' structures occur when an energy gain from combining different rather than overall alike conformations in a crystal significantly exceeds R·T, and this finding is discussed in the context of earlier explanations in the literature. Twinning is not related to archetype structures since it involves macroscopic domains of the same crystal structure. Archetype crystal structures are distinguished from crystal structure prediction trial structures in that an experimental reference structure is required for them. Categorization into archetype structures also has practical relevance, leading to a new practice of disorder modelling in experimental least-squares refinement alluded to in the above-mentioned publication.

Organic crystal structure prediction via coupled generative adversarial networks and graph convolutional networks

Organic crystal structures exert a profound impact on the physicochemical properties and biological effects of organic compounds. Quantum mechanics (QM)-based crystal structure predictions (CSPs) have somewhat alleviated the dilemma that experimental crystal structure investigations struggle to conduct complete polymorphism studies, but the high computing cost poses a challenge to its widespread application. The present study aims to construct DeepCSP, a feasible pure machine learning framework for minute-scale rapid organic CSP. Initially, based on 177,746 data entries from the Cambridge Crystal Structure Database, a generative adversarial network was built to conditionally generate trial crystal structures under selected feature constraints for the given molecule. Simultaneously, a graph convolutional attention network was used to predict the density of stable crystal structures for the input molecule. Subsequently, the distances between the predicted density and the definition-based calculated density would be considered to be the crystal structure screening and ranking basis, and finally, the density-based crystal structure ranking would be output. Two such distinct algorithms, performing the generation and ranking functionalities, respectively, collectively constitute the DeepCSP, which has demonstrated compelling performance in marketed drug validations, achieving an accuracy rate exceeding 80% and a hit rate surpassing 85%. Inspiringly, the computing speed of the pure machine learning methodology demonstrates the potential of artificial intelligence in advancing CSP research.

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.

Predicting crystal form stability under real-world conditions

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.

Solid Form Screenings in Pharmaceutical Development: a Perspective on Current Practices

Solid form screening is a crucial step in new drug development because solid forms of a drug substance significantly affect stability, dissolution and manufacturing processes of its drug products. This perspective introduces solid-state science from a practical standpoint, aiming to reduce knowledge gaps and promote communications among scientists with diverse background. This perspective starts with a concise overview that followed by discussion on timeline and goals of solid form screening. Techniques for solid from identification and characterization are then discussed. Subsequently, the perspective presents commonly used methods in solid form screening and introduces criteria and strategies to effectively select a favorable solid form based on screening results. The last section summarizes current practices in pharmaceutical industries and suggests potential opportunities for future research and development.

An unexpected high-pressure stability domain for a lower density polymorph of benzophenone

For over a century, it was thought that the crystalline polymorph II of benzophenone does not possess a stable domain in the pressure-temperature phase diagram. With a combination of new experimental results and literature data, this case of crystalline dimorphism has finally been solved and it is shown that form II possesses a stable domain at high pressure and high temperature, even though its density is lower than that of form I, the stable form under ordinary pressure and temperature conditions. The phase diagram of benzophenone is a clear demonstration of the fact that to understand the phase behaviour of a chemical substance both the exchange of heat (due to the change in intermolecular interactions) and work (due to the change of volume at a given pressure) need to be taken into account.

Phase diagrams-Why they matter and how to predict them

Understanding the thermodynamic stability and metastability of materials can help us to, for example, gauge whether crystalline polymorphs in pharmaceutical formulations are likely to be durable. It can also help us to design experimental routes to novel phases with potentially interesting properties. In this Perspective, we provide an overview of how thermodynamic phase behavior can be quantified both in computer simulations and machine-learning approaches to determine phase diagrams, as well as combinations of the two. We review the basic workflow of free-energy computations for condensed phases, including some practical implementation advice, ranging from the Frenkel-Ladd approach to thermodynamic integration and to direct-coexistence simulations. We illustrate the applications of such methods on a range of systems from materials chemistry to biological phase separation. Finally, we outline some challenges, questions, and practical applications of phase-diagram determination which we believe are likely to be possible to address in the near future using such state-of-the-art free-energy calculations, which may provide fundamental insight into separation processes using multicomponent solvents.

Minimal Basis Set Hartree-Fock Corrected with Atom-Centered Potentials for Molecular Crystal Modeling and Crystal Structure Prediction

Crystal structure prediction (CSP), determining the experimentally observable structure of a molecular crystal from the molecular diagram, is an important challenge with technologically relevant applications in materials manufacturing and drug design. For the purpose of screening the randomly generated candidate crystal structures, CSP protocols require energy ranking methods that are fast and can accurately capture the small energy differences between molecular crystals. In addition, a good ranking method should also produce accurate equilibrium geometries, both intramolecular and intermolecular. In this article, we explore the combination of minimal-basis-set Hartree-Fock (HF) with atom-centered potentials (ACPs) as a method for modeling the structure and energetics of molecular crystals. The ACPs are developed for the H, C, N, and O atoms and fitted to a set of reference data at the B86bPBE-XDM level in order to mitigate basis-set incompleteness and missing correlation. In particular, ACPs are developed in combination with two methods: HF-D3/MINIs and HF-3c. The application of ACPs greatly improves the performance of HF-D3/MINIs for lattice energies, crystal energy differences, energy-volume and energy-strain relations, and crystal geometries. In the case of HF-3c, the improvement in the crystal energy differences is much smaller than in HF-D3/MINIs, but lattice energies and particularly crystal geometries are considerably better when ACPs are used. The resulting methods may be useful for CSP but also for quick calculation of molecular crystal lattice energies and geometries.

Efficient Crystal Structure Prediction for Structurally Related Molecules with Accurate and Transferable Tailor-Made Force Fields

Crystal structure prediction (CSP) his generally used to complement experimental solid form screening and applied to individual molecules in drug development. The fast development of algorithms and computing resources offers the opportunity to use CSP earlier and for a broader range of applications in the drug design cycle. This study presents a novel paradigm of CSP specifically designed for structurally related molecules, referred to as Quick-CSP. The approach prioritizes more accurate physics through robust and transferable tailor-made force fields (TMFFs), such that significant efficiency gains are achieved through the reduction of expensive ab initio calculations. The accuracy of the TMFF is increased by the introduction of electrostatic multipoles, and the fragment-based force field parameterization scheme is demonstrated to be transferable for a family of chemically related molecules. The protocol is benchmarked with structurally related compounds from the Bromodomain and Extraterminal (BET) domain inhibitors series. A new convergence criterion is introduced that aims at performing only as many ab initio optimizations of crystal structures as required to locate the bottom of the crystal energy landscape within a user-defined accuracy. The overall approach provides significant cost savings ranging from three- to eight-fold less than the full-CSP workflow. The reported advancements expand the scope and utility of the underlying CSP building blocks as well as their novel reassembly to other applications earlier in the drug design cycle to guide molecule design and selection.

From Powders to Single Crystals: A Crystallographer's Toolbox for Small-Molecule Structure Determination

Although the crystal structures of small-molecule compounds are often determined from single-crystal X-ray diffraction (scXRD), recent advances in three-dimensional electron diffraction (3DED) and crystal structure prediction (CSP) methods promise to expand the structure elucidation toolbox available to the crystallographer. Herein, a comparative assessment of scXRD, 3DED, and CSP in combination with powder X-ray diffraction is carried out on two former drug candidate compounds and a multicomponent crystal of a key building block in the synthesis of gefapixant citrate.

Structure determination from unindexed powder data from scratch by a global optimization approach using pattern comparison based on cross-correlation functions

A method of ab initio crystal structure determination from powder diffraction data for organic and metal-organic compounds, which does not require prior indexing of the powder pattern, has been developed. Only a reasonable molecular geometry is required, needing knowledge of neither unit-cell parameters nor space group. The structures are solved from scratch by a global fit to the powder data using the new program FIDEL-GO (`FIt with DEviating Lattice parameters - Global Optimization'). FIDEL-GO uses a similarity measure based on cross-correlation functions, which allows the comparison of simulated and experimental powder data even if the unit-cell parameters deviate strongly. The optimization starts from large sets of random structures in various space groups. The unit-cell parameters, molecular position and orientation, and selected internal degrees of freedom are fitted simultaneously to the powder pattern. The optimization proceeds in an elaborate multi-step procedure with built-in clustering of duplicate structures and iterative adaptation of parameter ranges. The best structures are selected for an automatic Rietveld refinement. Finally, a user-controlled Rietveld refinement is performed. The procedure aims for the analysis of a wide range of `problematic' powder patterns, in particular powders of low crystallinity. The method can also be used for the clustering and screening of a large number of possible structure candidates and other application scenarios. Examples are presented for structure determination from unindexed powder data of the previously unknown structures of the nanocrystalline phases of 4,11-difluoro-, 2,9-dichloro- and 2,9-dichloro-6,13-dihydro-quinacridone, which were solved from powder patterns with 14-20 peaks only, and of the coordination polymer dichloro-bis(pyridine-N)copper(II).

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.

Rationalising the difference in crystallisability of two sulflowers using efficient in silico methods

The molecular structures of the first and second generation sulflowers, sulflower and persulfurated coronene (PSC), are remarkably similar: carbon ring structures decorated with sulfur atoms, without any additional moiety.

Along the road to Crystal Structure Prediction (CSP) of pharmaceutical-like molecules

Computational methods used for predicting crystal structures of organic compounds are mature enough to be routinely used with many rigid and semi-rigid organic molecules. The usefulness of Crystal Structure Prediction...

Pre-Nucleation Clusters Predict Crystal Structures in Models of Chiral Molecules

Kinetics can play an important role in the crystallization of molecules and can give rise to polymorphism, the tendency of molecules to form more than one crystal structure. Current computational methods of crystal structure prediction, however, focus almost exclusively on identifying the thermodynamically stable polymorph. Kinetic factors of nucleation and growth are often neglected because the underlying microscopic processes can be complex and accurate rate calculations are numerically cumbersome. In this work, we use molecular dynamics computer simulations to study simple molecular models that reproduce the crystallization behavior of real chiral molecules, including the formation of enantiopure and racemic crystals, as well as polymorphism. A significant fraction of these molecules forms crystals that do not have the lowest free energy. We demonstrate that at high supersaturation crystal formation can be accurately predicted by considering the similarities between oligomeric species in solution and molecular motifs in the crystal structure. For the case of racemic mixtures, we even find that knowledge of crystal free energies is not necessary and kinetic considerations are sufficient to determine if the system will undergo spontaneous chiral separation. Our results suggest conceptually simple ways of improving current crystal structure prediction methods.

Implications of the Conformationally Flexible, Macrocyclic Structure of the First-Generation, Direct-Acting Anti-Viral Paritaprevir on Its Solid Form Complexity and Chameleonic Behavior

Direct-acting antiviral regimens have transformed therapeutic management of hepatitis C across all prevalent genotypes. Most of the chemical matter in these regimens comprises molecules well outside the traditional drug development chemical space and presents significant challenges. Herein, the implications of high conformational flexibility and the presence of a 15-membered macrocyclic ring in paritaprevir are studied through a combination of advanced computational and experimental methods with focus on molecular chameleonicity and crystal form complexity. The ability of the molecule to toggle between high and low 3D polar surface area (PSA) conformations is underpinned by intramolecular hydrogen bonding (IMHB) interactions and intramolecular steric effects. Computational studies consequently show a very significant difference of over 75 Ų in 3D PSA between polar and apolar environments and provide the structural basis for the perplexingly favorable passive permeability of the molecule. Crystal packing and protein binding resulting in strong intermolecular interactions disrupt these intramolecular interactions. Crystalline Form I benefits from strong intermolecular interactions, whereas the weaker intermolecular interactions in Form II are partially compensated by the energetic advantage of an IMHB. Like Form I, no IMHB is observed within the receptor-bound conformation; instead, an intermolecular H-bond contributes to the potency of the molecule. The choice of metastable Form II is derisked through strategies accounting for crystal surface and packing features to manage higher form specific solid-state chemical reactivity and specific processing requirements. Overall, the results show an unambiguous link between structural features and derived properties from crystallization to dissolution, permeation, and docking into the protein pocket.

Improved quantitative crystal-structure comparison using powder diffractograms via anisotropic volume correction

A new anisotropic volume correction improves quantitative crystal structure comparison. Benchmarking against the 6th crystal structure prediction blind test data results in identification of two previously uncredited matching structures.

Molecular conformational evolution mechanism during nucleation of crystals in solution

Nucleation of crystals from solution is fundamental to many natural and industrial processes. In this work, the molecular mechanism of conformational polymorphism nucleation and the links between the molecular conformation in solutions and in crystals were investigated in detail by using 5-nitro-furazone as the model compound. Different polymorphs were prepared, and the conformations in solutions obtained by dissolving different polymorphs were analysed and compared. The solutions of 5-nitro-furazone were proven to contain multiple conformers through quantum chemical computation, Raman spectra analysis, 2D nuclear Overhauser effect spectroscopy spectra analysis and molecular dynamics simulation. The conformational evolution and desolvation path was illustrated according to the 1H NMR spectra of solutions with different concentrations. Finally, based on all the above analysis, the molecular conformational evolution path during nucleation of 5-nitro-furazone was illustrated. The results presented in this work shed a new light on the molecular mechanism of conformational polymorphism nucleation in solution.

Overcoming the difficulties of predicting conformational polymorph energetics in molecular crystals via correlated wavefunction methods

Molecular crystal structure prediction is increasingly being applied to study the solid form landscapes of larger, more flexible pharmaceutical molecules. Despite many successes in crystal structure prediction, van der Waals-inclusive density functional theory (DFT) methods exhibit serious failures predicting the polymorph stabilities for a number of systems exhibiting conformational polymorphism, where changes in intramolecular conformation lead to different intermolecular crystal packings. Here, the stabilities of the conformational polymorphs of o-acetamidobenzamide, ROY, and oxalyl dihydrazide are examined in detail. DFT functionals that have previously been very successful in crystal structure prediction perform poorly in all three systems, due primarily to the poor intramolecular conformational energies, but also due to the intermolecular description in oxalyl dihydrazide. In all three cases, a fragment-based dispersion-corrected second-order Møller-Plesset perturbation theory (MP2D) treatment of the crystals overcomes these difficulties and predicts conformational polymorph stabilities in good agreement with experiment. These results highlight the need for methods which go beyond current-generation DFT functionals to make crystal polymorph stability predictions truly reliable.

The Polymorphism of Drugs: New Approaches to the Synthesis of Nanostructured Polymorphs

Among the significant problems of modern pharmacology are the low solubility and bioavailability of drugs. One way to resolve this problem is to obtain new polymorphic forms of drugs with improved physicochemical properties. Various approaches have been developed with this aim, including the preparation of co-crystals, the use of nanoparticles, or the use of compounds in the form of a salt. A promising direction in pharmacology concerns the production of new stable polymorphic structures. In this mini-review, we consider certain aspects of drug polymorphism, methods for the synthesis of polymorphs, and the stability, size, and transformation of crystalline polymorphs. Moreover, we summarize our results from several studies demonstrating the problems associated with the synthesis of new polymorphous modifications based on inert gases and cryotemperatures. The results indicate that the problems specific to drug polymorphisms have only been partly resolved, are of current interest, and require further development.

Experimental and computational approaches to rationalise multicomponent supramolecular assemblies: dapsone monosolvates

The monosolvate crystal energy landscapes of dapsone (DDS) including the solvents carbon tetrachloride, acetone, cyclohexanone, dimethyl formamide, tetrahydrofuran, methyl ethyl ketone, 1,2-dichloroethane, 1,4-dioxane, dichloromethane and chloroform were established using experimental and computational approaches. To rationalise and understand solvate formation, solvate stability and desolvation reactions a careful control of the experimental crystallisation and storage conditions, a range of thermoanalytical methods and crystal structure prediction were required. Six of the eight DDS monosolvates are reported and characterised for the first time. Structural similarity and diversity of the at ambient conditions unstable monosolvates were apparent from the computed crystal energy landscapes, which had the experimental packings as lowest energy structures. The computed structures were used as input for Rietveld refinements and isostructurality of four of the monosolvates was confirmed. Packing comparisons of the solvate structures and molecular properties of the solvent molecules indicated that both size/shape of the solvent molecule and the possible DDSsolvent interactions are the important factors for DDS solvate formation. Through the combination of experiment and theory solvate stability and structural features have been rationalised.

A jumping crystal predicted with molecular dynamics and analysed with TLS refinement against powder diffraction data

By running a temperature series of molecular dynamics (MD) simulations starting from the known low-temperature phase, the experimentally observed phase transition in a 'jumping crystal' was captured, thereby providing a prediction of the unknown crystal structure of the high-temperature phase and clarifying the phase-transition mechanism. The phase transition is accompanied by a discontinuity in two of the unit-cell parameters. The structure of the high-temperature phase is very similar to that of the low-temperature phase. The anisotropic displacement parameters calculated from the MD simulations readily identified libration as the driving force behind the phase transition. Both the predicted crystal structure and the phase-transition mechanism were verified experimentally using TLS (translation, libration, screw) refinement against X-ray powder diffraction data.

Accuracy and reproducibility in crystal structure prediction: the curious case of ROY

Because of excessive electron delocalization, the polymorphs of ROY constitute a surprisingly challenging system for crystal structure prediction.