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.)

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

While density-functional theory (DFT) remains one of the most widely used tools in computational chemistry, most functionals fail to properly account for the effects of London dispersion. Hence, there are many popular post-self-consistent methods to add a dispersion correction to the DFT energy. Until now, the most popular methods have never been compared on equal footing due to not being implemented in the same electronic structure packages. In this work, we performed a large-scale benchmarking study, directly comparing the accuracy of the exchange-hole dipole moment (XDM), D3BJ, D4, TS, MBD, and MBD-NL dispersion models when applied to the recent DES15K database of nearly 15,000 molecular complexes at both expanded and compressed geometries. Our study showed similarly good performance for all dispersion methods (except TS) when applied to neutral complexes. However, they all performed worse for ionic complexes, particularly those involving dications of alkaline earth metals, due to systematic overbinding by the base PBE0 density functional. Investigation of the largest outliers also revealed that only the MBD and MBD-NL methods demonstrate surprising errors for complexes involving alkali metal cations at compressed geometries where they tended to significantly overbind. As we would expect minimal dispersion binding for such complexes, we further investigated the origins of these errors for the potential energy curve of a model cation-π complex. Overall, there is little choice between the XDM, D3BJ, D4, MBD, and MBD-NL dispersion methods for most systems. However, the MBD-based methods are not recommended for complexes involving organic species and alkali or alkaline earth metal cations, for example when modeling Li + intercalation into graphite.

Molecular crystals are important for many applications, including energetic materials, organic semiconductors, and the development and commercialization of pharmaceuticals. The exchange-hole dipole moment (XDM) dispersion model has shown good performance in the calculation of relative and absolute lattice energies of molecular crystals, although it has traditionally been applied in combination with plane-wave/pseudopotential approaches. This has limited XDM to use with semilocal functional approximations, which suffer from delocalization error and poor quality conformational energies, and to systems with a few hundreds of atoms at most due to unfavorable scaling. In this work, we combine XDM with numerical atomic orbitals, which enable the efficient use of XDM-corrected hybrid functionals for molecular crystals. We test the new XDM-corrected functionals for their ability to predict the lattice energies of molecular crystals for the X23 set and 13 ice phases, the latter being a particularly stringent test. A composite approach using a XDM-corrected, 25% hybrid functional based on B86bPBE achieves a mean absolute error of 0.48 kcal mol -1 per molecule for the X23 set and 0.19 kcal mol -1 for the total lattice energies of the ice phases, compared to recent diffusion Monte-Carlo data. These results make the new XDM-corrected hybrids not only far more computationally efficient than previous XDM implementations, but also the most accurate density-functional methods for molecular crystal lattice energies to date., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

Despite dominating industrial processes, heterogeneous catalysts remain challenging to characterize and control. This is largely attributable to the diversity of potentially active sites at the catalyst-reactant interface and the complex behaviour that can arise from interactions between active sites. Surface-supported, single-site molecular catalysts aim to bring together benefits of both heterogeneous and homogeneous catalysts, offering easy separability while exploiting molecular design of reactivity, though the presence of a surface is likely to influence reaction mechanisms. Here, we use metal-organic coordination to build reactive Fe-terpyridine sites on the Ag(111) surface and study their activity towards CO and C 2 H 4 gaseous reactants using low-temperature ultrahigh-vacuum scanning tunnelling microscopy, scanning tunnelling spectroscopy, and atomic force microscopy supported by density-functional theory models. Using a site-by-site approach at low temperature to visualize the reaction pathway, we find that reactants bond to the Fe-tpy active sites via surface-bound intermediates, and investigate the role of the substrate in understanding and designing single-site catalysts on metallic supports., (© 2022. The Author(s).)

Many crystal structure prediction protocols only concern themselves with the electronic energy of molecular crystals. However, vibrational contributions to the free energy (F vib ) can be significant in determining accurate stability rankings for crystal candidates. While force-field studies have been conducted to gauge the magnitude of these free-energy corrections, highly accurate results from quantum mechanical methods, such as density-functional theory (DFT), are desirable. Here, we introduce the PV17 set of 17 polymorphic pairs of organic molecular crystals, for which plane wave DFT is used to calculate the vibrational free energies and free-energy differences (ΔF vib ) between each pair. Our DFT results confirm that the vibrational free-energy corrections are small, having a mean value of 1.0 kJ/mol and a maximum value of 2.3 kJ/mol for the PV17 set. Furthermore, we assess the accuracy of a series of lower-cost DFT, semi-empirical, and force-field models for computing ΔF vib that have been proposed in the literature. It is found that calculating F vib using the Γ-point frequencies does not provide ΔF vib values of sufficiently high quality. In addition, ΔF vib values calculated using various approximate methods have mean absolute errors relative to our converged DFT results of equivalent or larger magnitude than the vibrational free-energy corrections themselves. Thus, we conclude that, in a crystal structure prediction protocol, it is preferable to forego the inclusion of vibrational free-energy corrections than to estimate them with any of the approximate methods considered here.

Post-self-consistent dispersion corrections are now the norm when applying density-functional theory to systems where non-covalent interactions play an important role. However, there is a wide range of base functionals and dispersion corrections available from which to choose. In this work, we opine on the most desirable requirements to ensure that both the base functional and dispersion correction, individually, are as accurate as possible for non-bonded repulsion and dispersion attraction. The base functional should be dispersionless, numerically stable, and involve minimal delocalization error. Simultaneously, the dispersion correction should include finite damping, higher-order pairwise dispersion terms, and electronic many-body effects. These criteria are essential for avoiding reliance on error cancellation and obtaining correct results from correct physics.