Abstract:Molecular crystals occur widely in pharmaceuticals, foods, explosives, organic semiconductors, and many other applications. Thanks to substantial progress in electronic structure modeling of molecular crystals, attention is now shifting from basic crystal structure prediction and lattice energy modeling toward the accurate prediction of experimentally observable properties at finite temperatures and pressures. This Account discusses how fragment-based electronic structure methods can be used to model a variety of experimentally relevant molecular crystal properties. First, it describes the coupling of fragment electronic structure models with quasi-harmonic techniques for modeling the thermal expansion of molecular crystals, and what effects this expansion has on thermochemical and mechanical properties. Excellent agreement with experiment is demonstrated for the molar volume, sublimation enthalpy, entropy, and free energy, and the bulk modulus of phase I carbon dioxide when large basis second-order Møller-Plesset perturbation theory (MP2) or coupled cluster theories (CCSD(T)) are used. In addition, physical insight is offered into how neglect of thermal expansion affects these properties. Zero-point vibrational motion leads to an appreciable expansion in the molar volume; in carbon dioxide, it accounts for around 30% of the overall volume expansion between the electronic structure energy minimum and the molar volume at the sublimation point. In addition, because thermal expansion typically weakens the intermolecular interactions, neglecting thermal expansion artificially stabilizes the solid and causes the sublimation enthalpy to be too large at higher temperatures. Thermal expansion also frequently weakens the lower-frequency lattice phonon modes; neglecting thermal expansion causes the entropy of sublimation to be overestimated. Interestingly, the sublimation free energy is less significantly affected by neglecting thermal expansion because the systematic errors in the enthalpy and entropy cancel somewhat. Second, because solid state nuclear magnetic resonance (NMR) plays an increasingly important role in molecular crystal studies, this Account discusses how fragment methods can be used to achieve higher-accuracy chemical shifts in molecular crystals. Whereas widely used plane wave density functional theory models are largely restricted to generalized gradient approximation (GGA) functionals like PBE in practice, fragment methods allow the routine use of hybrid density functionals with only modest increases in computational cost. In extensive molecular crystal benchmarks, hybrid functionals like PBE0 predict chemical shifts with 20-30% higher accuracy than GGAs, particularly for 1H, 13C, and 15N nuclei. Due to their higher sensitivity to polarization effects, 17O chemical shifts prove slightly harder to predict with fragment methods. Nevertheless, the fragment model results are still competitive with those from GIPAW. The improved accuracy achievable with fragment approaches and hybrid density functionals increases discrimination between different potential assignments of individual shifts or crystal structures, which is critical in NMR crystallography applications. This higher accuracy and greater discrimination are highlighted in application to the solid state NMR of different acetaminophen and testosterone crystal forms.

Prediction of the Crystalline Densities of Aliphatic Nitrates by Quantum Chemistry Methods

Reliable prediction of crystal density of high nitrogen-containing organic compounds as powerful, less sensitive, eco-friendly energetic materials for dependable assessment of their performance

Benchmark Study of Density Functional Theory Methods in Geometry Optimization of Transition Metal-Dinitrogen Complexes

Theoretical studies on the thermodynamic properties, densities, detonation properties, and pyrolysis mechanisms of trinitromethyl-substituted aminotetrazole compounds

Theoretical Study of Heterocyclic Nitro Compounds for Use As Novel Explosives and Propellants

A novel all-nitrogen molecular crystal N16 as a promising high-energy density material

First-principle calculations of the electronic, vibrational, and thermodynamic properties of nitrogen-rich salt of 3,6-dinitramino-1,2,4,5-tetrazine [(NH4)2(DNAT)]

Density Prediction Models for Energetic Compounds Merely Using Molecular Topology

Accurate Predictions of the Nmr Parameters in Organic and Biological Crystallines

Density functional methods studies on geometry, vibrational frequencies and thermodynamic properties of propyl nitrate

Ab Initio Neural Network MD Simulation of Thermal Decomposition of a High Energy Material CL-20/TNT.

New approach for predicting crystal densities of energetic materials

A simple, fast and convenient new method for predicting the stability of nitro compounds

Theoretical simulation study on crystal property and hygroscopicity of ADN doping with nitramine explosives (RDX, HMX, and CL-20)

Theoretical Prediction on Properties of 3,4-Bisnitrofurazanfuroxan (DNTF) Crystal and its Polymer Bonded Explosives (PBXs) Through Molecular Dynamics (MD) Simulation

ReaxFF Molecular Dynamics Simulation of Hydrostatic and Uniaxial Compression of Nitrate Energetic Materials

A theoretical investigation on the chemical environment of pyrazine-2,3-dicarbonitrile and phthalonitrile: Density functional theory (DFT) calculation and experimental verification

Molecular Dynamic Simulation and Density Functional Theory Insight into the Nitrogen Rich Explosive 1,5-Diaminotetrazole(dat)

The interplay of density functional selection and crystal structure for accurate NMR chemical shift predictions

Theoretical Prediction of Structures and Properties of 2,4,6-Trinitro-1,3,5-Triazine (TNTA) Green Energetic Materials from DFT and ReaxFF Molecular Modeling

Predicting Molecular Crystal Properties from First Principles: Finite-Temperature Thermochemistry to NMR Crystallography