Search for: All records

The ability to theoretically predict accurate NMR chemical shifts in solids is increasingly important due to the role such shifts play in selecting among proposed model structures. Herein, two theoretical methods are evaluated for their ability to assign¹⁵N shifts from guanosine dihydrate to one of the two independent molecules present in the lattice. The NMR data consist of¹⁵N shift tensors from 10 resonances. Analysis using periodic boundary or fragment methods consider a benchmark dataset to estimate errors and predict uncertainties of 5.6 and 6.2 ppm, respectively. Despite this high accuracy, only one of the five sites were confidently assigned to a specific molecule of the asymmetric unit. This limitation is not due to negligible differences in experimental data, as most sites exhibit differences of >6.0 ppm between pairs of resonances representing a given position. Instead, the theoretical methods are insufficiently accurate to make assignments at most positions.

 

Photomechanical molecular crystals have garnered attention for their ability to transform light into mechanical work, but difficulties in characterizing the structural changes and mechanical responses experimentally have hindered the development of practical organic crystal engines. This study proposes a new computational framework for predicting the solid-state crystal-to-crystal photochemical transformations entirely from first principles, and it establishes a photomechanical engine cycle that quantifies the anisotropic mechanical performance resulting from the transformation. The approach relies on crystal structure prediction, solid-state topochemical principles, and high-quality electronic structure methods. After validating the framework on the well-studied [4 + 4] cycloadditions in 9-methyl anthracene and 9- tert -butyl anthracene ester, the experimentally-unknown solid-state transformation of 9-carboxylic acid anthracene is predicted for the first time. The results illustrate how the mechanical work is done by relaxation of the crystal lattice to accommodate the photoproduct, rather than by the photochemistry itself. The large ∼10 7 J m −3 work densities computed for all three systems highlight the promise of photomechanical crystal engines. This study demonstrates the importance of crystal packing in determining molecular crystal engine performance and provides tools and insights to design improved materials in silico . 

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.

 

Second-order Møller–Plesset perturbation theory (MP2) provides a valuable alternative to density functional theory for modeling problems in organic and biological chemistry. However, MP2 suffers from known limitations in the description of van der Waals (London) dispersion interactions and reaction thermochemistry. Here, a spin-component-scaled, dispersion-corrected MP2 model (SCS-MP2D) is proposed that addresses these weaknesses. The dispersion correction, which is based on Grimme's D3 formalism, replaces the uncoupled Hartree–Fock dispersion inherent in MP2 with a more robust coupled Kohn–Sham treatment. The spin-component scaling of the residual MP2 correlation energy then reduces the remaining errors in the model. This two-part correction strategy solves the problem found in earlier spin-component-scaled MP2 models where completely different spin-scaling parameters were needed for describing reaction energies versus intermolecular interactions. Results on 18 benchmark data sets and two challenging potential energy curves demonstrate that SCS-MP2D considerably improves upon the accuracy of MP2 for intermolecular interactions, conformational energies, and reaction energies. Its accuracy and computational cost are competitive with state-of-the-art density functionals such as DSD-BLYP-D3(BJ), revDSD-PBEP86-D3(BJ), ωB97X-V, and ωB97M-V for systems with ∼100 atoms. 

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.