Conformational polymorphs of organic molecular crystals represent a challenging test for quantum chemistry because they require careful balancing of the intra- and intermolecular interactions. This study examines 54 molecular conformations from 20 sets of conformational polymorphs, along with the relative lattice energies and 173 dimer interactions taken from six of the polymorph sets. These systems are studied with a variety of van der Waals-inclusive density functionals theory models; dispersion-corrected spin-component-scaled second-order Møller–Plesset perturbation theory (SCS-MP2D); and domain local pair natural orbital coupled cluster singles, doubles, and perturbative triples [DLPNO-CCSD(T)]. We investigate how delocalization error in conventional density functionals impacts monomer conformational energies, systematic errors in the intermolecular interactions, and the nature of error cancellation that occurs in the overall crystal. The density functionals B86bPBE-XDM, PBE-D4, PBE-MBD, PBE0-D4, and PBE0-MBD are found to exhibit sizable one-body and two-body errors vs DLPNO-CCSD(T) benchmarks, and the level of success in predicting the relative polymorph energies relies heavily on error cancellation between different types of intermolecular interactions or between intra- and intermolecular interactions. The SCS-MP2D and, to a lesser extent, ωB97M-V models exhibit smaller errors and rely less on error cancellation. Implications for crystal structure prediction of flexible compounds are discussed. Finally, the one-body and two-body DLPNO-CCSD(T) energies taken from these conformational polymorphs establish the CP1b and CP2b benchmark datasets that could be useful for testing quantum chemistry models in challenging real-world systems with complex interplay between intra- and intermolecular interactions, a number of which are significantly impacted by delocalization error.
- Award ID(s):
- 1665212
- NSF-PAR ID:
- 10165368
- Date Published:
- Journal Name:
- Chemical Science
- Volume:
- 11
- Issue:
- 8
- ISSN:
- 2041-6520
- Page Range / eLocation ID:
- 2200 to 2214
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
Molecular crystal engineering seeks to tune the material properties by controlling the crystal packing. However, the range of achievable properties is constrained by the limited energy range of polymorphs which can be crystallized. Here, computational modeling highlights that a solid-state crystal-to-crystal chemical reaction in 9- tert -butyl anthracene ester (9TBAE) nanorods [Al-Kaysi et al. , J. Am. Chem. Soc. , 2006, 128 , 15938] imparts “synthetic memory” into the crystal structure that allows reproducible formation of a highly metastable, yet long-lived polymorph. Specifically, whereas the vast majority of known polymorphs exhibit lattice energy differences below 10 kJ mol −1 , the conformational polymorph formed via solid state reaction chemistry lies 14 kJ mol −1 higher in energy than the form grown from solution, according to calculations that combine a dispersion-corrected second-order Møller–Plesset perturbation theory (MP2D) treatment of the monomer and photodimer with a density functional theory treatment (B86bPBE-XDM) of the intermolecular interactions in the crystal. Moreover, the solid-state reaction environment traps a highly unstable intramolecular photodimer conformation which defies the conventional wisdom surrounding conformational polymorphs. These observations suggest that solid-state reaction chemistry represents an under-appreciated strategy for producing polymorphs that would likely be unobtainable otherwise.more » « less
-
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.more » « less
-
more » « less
Abstract We present a graph‐theoretic approach to adaptively compute many‐body approximations in an efficient manner to perform (a) accurate post‐Hartree–Fock (HF) ab initio molecular dynamics (AIMD) at density functional theory (DFT) cost for medium‐ to large‐sized molecular clusters, (b) hybrid DFT electronic structure calculations for condensed‐phase simulations at the cost of pure density functionals, (c) reduced‐cost on‐the‐fly basis extrapolation for gas‐phase AIMD and condensed phase studies, and (d) accurate post‐HF‐level potential energy surfaces at DFT cost for quantum nuclear effects. The salient features of our approach are ONIOM‐like in that (a) the full system (cluster or condensed phase) calculation is performed at a lower level of theory (pure DFT for condensed phase or hybrid DFT for molecular systems), and (b) this approximation is improved through a correction term that captures all many‐body interactions up to any given order within a higher level of theory (hybrid DFT for condensed phase; CCSD or MP2 for cluster), combined through graph‐theoretic methods. Specifically, a region of chemical interest is coarse‐grained into a set of nodes and these nodes are then connected to form edges based on a given definition of local envelope (or threshold) of interactions. The nodes and edges together define a graph, which forms the basis for developing the many‐body expansion. The methods are demonstrated through (a) ab initio dynamics studies on protonated water clusters and polypeptide fragments, (b) potential energy surface calculations on one‐dimensional water chains such as those found in ion channels, and (c) conformational stabilization and lattice energy studies on homogeneous and heterogeneous surfaces of water with organic adsorbates using two‐dimensional periodic boundary conditions.
-
In an attempt to grow 8-hydroxyquinoline–acetaminophen co-crystals from equimolar amounts of conformers in a chloroform–ethanol solvent mixture at room temperature, the title compound, C 9 H 7 NO, was obtained. The molecule is planar, with the hydroxy H atom forming an intramolecular O—H...N hydrogen bond. In the crystal, molecules form centrosymmetric dimers via two O—H...N hydrogen bonds. Thus, the hydroxy H atoms are involved in bifurcated O—H...N hydrogen bonds, leading to the formation of a central planar four-membered N 2 H 2 ring. The dimers are bound by intermolecular π–π stacking [the shortest C...C distance is 3.2997 (17) Å] and C—H...π interactions into a three-dimensional framework. The crystal grown represents a new monoclinic polymorph in the space group P 2 1 / n . The molecular structure of the present monoclinic polymorph is very similar to that of the orthorhombic polymorph (space group Fdd 2) studied previously [Roychowdhury et al. (1978). Acta Cryst. B 34 , 1047–1048; Banerjee & Saha (1986). Acta Cryst. C 42 , 1408–1411]. The structures of the two polymorphs are distinguished by the different geometries of the hydrogen-bonded dimers, which in the crystal of the orthorhombic polymorph possess twofold axis symmetry, with the central N 2 H 2 ring adopting a butterfly conformation.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/C9SC05689K
