Search

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., (open access.)

A seventh blind test of crystal structure prediction was organized by the Cambridge Crystallographic Data Centre featuring seven target systems of varying complexity: a silicon and iodine-containing molecule, a copper coordination complex, a near-rigid molecule, a cocrystal, a polymorphic small agrochemical, a highly flexible polymorphic drug candidate, and a polymorphic morpholine salt. In this first of two parts focusing on structure generation methods, many crystal structure prediction (CSP) methods performed well for the small but flexible agrochemical compound, successfully reproducing the experimentally observed crystal structures, while few groups were successful for the systems of higher complexity. A powder X-ray diffraction (PXRD) assisted exercise demonstrated the use of CSP in successfully determining a crystal structure from a low-quality PXRD pattern. The use of CSP in the prediction of likely cocrystal stoichiometry was also explored, demonstrating multiple possible approaches. Crystallographic disorder emerged as an important theme throughout the test as both a challenge for analysis and a major achievement where two groups blindly predicted the existence of disorder for the first time. Additionally, large-scale comparisons of the sets of predicted crystal structures also showed that some methods yield sets that largely contain the same crystal structures., (open access.)

The Hohenberg-Kohn theorem of density-functional theory establishes the existence of a bijection between the ground-state electron density and the external potential of a many-body system. This guarantees a one-to-one map from the electron density to all observables of interest including electronic excited-state energies. Time-Dependent Density-Functional Theory (TDDFT) provides one framework to resolve this map; however, the approximations inherent in practical TDDFT calculations, together with their computational expense, motivate finding a cheaper, more direct map for electronic excitations. Here, we show that determining density and energy functionals via machine learning allows the equations of TDDFT to be bypassed. The framework we introduce is used to perform the first excited-state molecular dynamics simulations with a machine-learned functional on malonaldehyde and correctly capture the kinetics of its excited-state intramolecular proton transfer, allowing insight into how mechanical constraints can be used to control the proton transfer reaction in this molecule. This development opens the door to using machine-learned functionals for highly efficient excited-state dynamics simulations., (© 2022. The Author(s).)

Imidacloprid, the world's leading insecticide, has been approved recently for controlling infectious disease vectors; yet, in agricultural settings, it has been implicated in the frightening decline of pollinators. This argues for strategies that sharply reduce the environmental impact of imidacloprid. When used as a contact insecticide, the effectiveness of imidacloprid relies on physical contact between its crystal surfaces and insect tarsi. Herein, seven new imidacloprid crystal polymorphs are reported, adding to two known forms. Anticipating that insect uptake of imidacloprid molecules would depend on the respective free energies of crystal polymorph surfaces, measurements of insect knockdown times for the metastable crystal forms were as much as nine times faster acting than the commercial form against Aedes, Anopheles , and Culex mosquitoes as well as Drosophila (fruit flies). These results suggest that replacement of commercially available imidacloprid crystals (a.k.a. Form I) in space-spraying with any one of three new polymorphs, Forms IV, VI, IX, would suppress vector-borne disease transmission while reducing environmental exposure and harm to nontarget organisms.

Predicting structures of organic molecular cocrystals is a challenging task when considering the immense number of possible intermolecular orientations. Use of the Shannon information entropy, constructed from an intermolecular orientational spatial distribution function, to drive a search for crystal structures via enhanced molecular dynamics can be an efficient way to map out a landscape of putative polymorphs. Here, the Shannon entropy is used to generate a set of collective variables for differentiating polymorphs of a 1:1 cocrystal of resorcinol and urea. We show that driven adiabatic free energy dynamics, a particular enhanced-sampling approach, combined with these entropy variables, can transform the stable phase into alternate polymorphs. Density functional theory calculations confirm that a structure obtained from the enhanced molecular dynamics is stable at pressures above 1 GPa. We thus show that enhanced sampling should be considered an integral component of crystal structure searching protocols for systems with multiple independent molecules.

Kohn-Sham density functional theory (DFT) is a standard tool in most branches of chemistry, but accuracies for many molecules are limited to 2-3 kcal ⋅ mol -1 with presently-available functionals. Ab initio methods, such as coupled-cluster, routinely produce much higher accuracy, but computational costs limit their application to small molecules. In this paper, we leverage machine learning to calculate coupled-cluster energies from DFT densities, reaching quantum chemical accuracy (errors below 1 kcal ⋅ mol -1 ) on test data. Moreover, density-based Δ-learning (learning only the correction to a standard DFT calculation, termed Δ-DFT ) significantly reduces the amount of training data required, particularly when molecular symmetries are included. The robustness of Δ-DFT  is highlighted by correcting "on the fly" DFT-based molecular dynamics (MD) simulations of resorcinol (C 6 H 4 (OH) 2 ) to obtain MD trajectories with coupled-cluster accuracy. We conclude, therefore, that Δ-DFT  facilitates running gas-phase MD simulations with quantum chemical accuracy, even for strained geometries and conformer changes where standard DFT fails.

A method for calculating the free energy difference between two structurally defined conformational states of a chemical system is developed. A path is defined using a previously reported collective variable that interpolates between two or more conformations, and a restraint is introduced in order to keep the system close to the path. The evolution of the system along the path, which typically presents a high free energy barrier, is generated using enhanced sampling schemes. Although the formulation of the method in terms of a path is quite general, an important advance in this work is the demonstration that prior knowledge of the path is, in fact, not needed and that the free energy difference can be obtained using a simplified definition of the path collective variable that only involves the endpoints. We first validate this method on cyclohexane isomerization. The method is then tested for an extensive conformational change in a realistic molecular system by calculating the free energy difference between the α -helix and β -hairpin conformations of deca-alanine in solution. Finally, the method is applied to a biologically relevant system to calculate the free energy difference of an observed and a hypothetical conformation of an antigenic peptide bound to a major histocompatibility complex.