An official website of the United States government
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
Rubrene untwisted: common density functional theory calculations overestimate its deviant tendencies
The exceptionally high carrier mobility of rubrene derives from the combination of its intrinsic electronic properties and favorable crystal packing that facilitates charge transport. Unlike the planar conformations adopted by rubrene single crystals, however, many rubrene derivatives crystallize with a twisted tetracene core and exhibit poor carrier mobility. Typical density functional theory (DFT) calculations suggest that the twisted conformation is preferred by ∼10–14 kJ mol −1 or more in the gas phase. However, the present work shows that those calculations overestimate the twisting energy by several kJ mol −1 due to density-driven delocalization error, and that the twisting energies are actually only ∼8–10 kJ mol −1 for typical rubrene derivatives when computed with higher-level correlated wave function models. This result has two significant implications for crystal engineering with rubrene derivatives: first, DFT calculations can erroneously predict polymorphs containing twisted rubrene conformations to be more stable, when in fact structures with planar conformations are preferred, as is demonstrated here for perfluororubrene. Second, the smaller twisting energies make it more likely that solid form screening could discover new planar-core polymorphs of rubrene derivatives that have previously been crystallized only in a twisted conformation. These in turn might exhibit better organic semiconducting properties.
more »
« less
- Award ID(s):
- 1955554
- PAR ID:
- 10252565
- Date Published:
- Journal Name:
- Journal of Materials Chemistry C
- Volume:
- 9
- Issue:
- 8
- ISSN:
- 2050-7526
- Page Range / eLocation ID:
- 2848 to 2857
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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.more » « less
-
Abstract The crystal structures of the charge‐transfer (CT) cocrystals formed by the π‐electron acceptor 1,3,4,5,7,8‐hexafluoro‐11,11,12,12‐tetracyanonaphtho‐2,6‐quinodimethane (F6TNAP) with the planar π‐electron‐donor molecules triphenylene (TP), benzo[b]benzo[4,5]thieno[2,3‐d]thiophene (BTBT), benzo[1,2‐b:4,5‐b′]dithiophene (BDT), pyrene (PY), anthracene (ANT), and carbazole (CBZ) have been determined using single‐crystal X‐ray diffraction (SCXRD), along with those of two polymorphs of F6TNAP. All six cocrystals exhibit 1:1 donor/acceptor stoichiometry and adopt mixed‐stacking motifs. Cocrystals based on BTBT and CBZ π‐electron donor molecules exhibit brickwork packing, while the other four CT cocrystals show herringbone‐type crystal packing. Infrared spectroscopy, molecular geometries determined by SCXRD, and electronic structure calculations indicate that the extent of ground‐state CT in each cocrystal is small. Density functional theory calculations predict large conduction bandwidths and, consequently, low effective masses for electrons for all six CT cocrystals, while the TP‐, BDT‐, and PY‐based cocrystals are also predicted to have large valence bandwidths and low effective masses for holes. Charge‐carrier mobility values are obtained from space‐charge limited current (SCLC) measurements and field‐effect transistor measurements, with values exceeding 1 cm2V−1s1being estimated from SCLC measurements for BTBT:F6TNAP and CBZ:F6TNAP cocrystals.more » « less
-
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−1at 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.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/D0TC05463A
