Characterization of Complicated New Polymorphs of Chlorothalonil by X-ray Diffraction and Computer Crystal Structure Prediction

Isostructurality of quinoxaline crystal phases: the interplay of weak hydrogen bonds and halogen bonding

Creation of isostructural systems is a balance between thermodynamic and kinetic factors, shown by a set of substituted quinoxalines, where the lighter halogens form a set of metastable isostructural crystals, different to the heavy substitutions.

The use of variable temperature 13 C solid-state MAS NMR and GIPAW DFT calculations to explore the dynamics of diethylcarbamazine citrate

Experimental ¹³ C solid-state magic-angle spinning (MAS) Nuclear Magnetic Resonance (NMR) as well as Density-Functional Theory (DFT) gauge-including projector augmented wave (GIPAW) calculations were used to probe disorder and local mobility in diethylcarbamazine citrate, (DEC)⁺ (citrate)^- . This compound has been used as the first option drug for the treatment of filariasis, a disease endemic in tropical countries and caused by adult worms of Wuchereria bancrofti, which is transmitted by mosquitoes. We firstly present 2D ¹³ C─¹ H dipolar-coupling-mediated heteronuclear correlation spectra recorded at moderate spinning frequency, to explore the intermolecular interaction between DEC and citrate molecules. Secondly, we investigate the dynamic behavior of (DEC)⁺ (citrate)^- by varying the temperature and correlating the experimental MAS NMR results with DFT GIPAW calculations that consider two (DEC)⁺ conformers (in a 70:30 ratio) for crystal structures determined at 293 and 235 K. Solid-state NMR provides insights on slow exchange dynamics revealing conformational changes involving particularly the DEC ethyl groups.

Modeling Polymorphic Molecular Crystals with Electronic Structure Theory

Interest in molecular crystals has grown thanks to their relevance to pharmaceuticals, organic semiconductor materials, foods, and many other applications. Electronic structure methods have become an increasingly important tool for modeling molecular crystals and polymorphism. This article reviews electronic structure techniques used to model molecular crystals, including periodic density functional theory, periodic second-order Møller-Plesset perturbation theory, fragment-based electronic structure methods, and diffusion Monte Carlo. It also discusses the use of these models for predicting a variety of crystal properties that are relevant to the study of polymorphism, including lattice energies, structures, crystal structure prediction, polymorphism, phase diagrams, vibrational spectroscopies, and nuclear magnetic resonance spectroscopy. Finally, tools for analyzing crystal structures and intermolecular interactions are briefly discussed.

Why don't we find more polymorphs?

Crystal structure prediction (CSP) studies are not limited to being a search for the most thermodynamically stable crystal structure, but play a valuable role in understanding polymorphism, as shown by interdisciplinary studies where the crystal energy landscape has been explored experimentally and computationally. CSP usually produces more thermodynamically plausible crystal structures than known polymorphs. This article illustrates some reasons why: because (i) of approximations in the calculations, particularly the neglect of thermal effects (see §1.1); (ii) of the molecular rearrangement during nucleation and growth (see §1.2); (iii) the solid-state structures observed show dynamic or static disorder, stacking faults, other defects or are not crystalline and so represent more than one calculated structure (see §1.3); (iv) the structures are metastable relative to other molecular compositions (see §1.4); (v) the right crystallization experiment has not yet been performed (see §1.5) or (vi) cannot be performed (see §1.6) and the possibility (vii) that the polymorphs are not detected or structurally characterized (see §1.7). Thus, we can only aspire to a general predictive theory for polymorphism, as this appears to require a quantitative understanding of the kinetic factors involved in all possible multi-component crystallizations. For a specific molecule, analysis of the crystal energy landscape shows the potential complexity of its crystallization behaviour.

The complexity of hydration of phloroglucinol: a comprehensive structural and thermodynamic characterization

Hydrate formation is of great importance as the inclusion of water molecules affects many solid state properties and hence determines the required chemical processing, handling, and storage. Phloroglucinol is industrially important, and the observed differences in the morphology and diffuse scattering effects with growth conditions have been scientifically controversial. We have studied the anhydrate and dihydrate of phloroglucinol and their transformations by a unique combination of complementary experimental and computational techniques, namely, moisture sorption analysis, hot-stage microscopy, differential scanning calorimetry, thermogravimetry, isothermal calorimetry, single crystal and powder X-ray diffractometry, and crystal energy landscape calculations. The enthalpically stable dihydrate phase is unstable below 16% relative humidity (25 °C) and above 50 °C (ambient humidity), and the kinetics of hydration/dehydration are relatively rapid with a small hysteresis. A consistent atomistic picture of the thermodynamics of the hydrate/anhydrate transition was derived, consistent with the disordered single X-ray crystal structure and crystal energy landscape showing closely related low energy hydrate structures. These structures provide models for proton disorder and show stacking faults as intergrowth of different layers are possible. This indicates that the consequent variability in crystal surface features and diffuse scattering with growth conditions is not a practical concern.

Substitutional and orientational disorder in organic crystals: a symmetry-adapted ensemble model

Modelling of disorder in organic crystals is highly desirable since it would allow thermodynamic stabilities and other disorder-sensitive properties to be estimated for such systems. Two disordered organic molecular systems are modeled using a symmetry-adapted ensemble approach, in which the disordered system is treated as an ensemble of the configurations of a supercell with respect to substitution of one disorder component for another. Computation time is kept manageable by performing calculations only on the symmetrically inequivalent configurations. Calculations are presented on a substitutionally disordered system, the dichloro/dibromobenzene solid solution, and on an orientationally disordered system, eniluracil, and the resultant free energies, disorder patterns, and system properties are discussed. The results are found to be in agreement with experiment, when some physically implausible configurations are removed from the ensemble average for eniluracil, highlighting the dangers of a completely automated approach to organic crystal thermodynamics which ignores the barriers to equilibration once the crystal has been formed.

Successful prediction of a model pharmaceutical in the fifth blind test of crystal structure prediction

The range of target structures in the fifth international blind test of crystal structure prediction was extended to include a highly flexible molecule, (benzyl-(4-(4-methyl-5-(p-tolylsulfonyl)-1,3-thiazol-2-yl)phenyl)carbamate, as a challenge representative of modern pharmaceuticals. Two of the groups participating in the blind test independently predicted the correct structure. The methods they used are described and contrasted, and the implications of the capability to tackle molecules of this complexity are discussed.

DFT computations and spectroscopic analysis of a pesticide: chlorothalonil

NIR FT-Raman and IR spectra of the biologically active molecule, chlorothalonil have been recorded and analyzed. The molecular geometry and vibrational wavenumbers of the title compound in the ground state have been calculated by density functional theory (DFT) with 6-31G(d) basis set. In order to obtain the information about the influence of intramolecular interaction on the molecule, the calculated geometries of chlorothalonil molecule was compared with the experimental data. The results of the optimized molecular structure gave clear evidence for the intramolecular charge transfer (ICT). Time-dependent density functional theory (TD-DFT) calculation of the electronic spectra has been performed and compared with the experimental UV-visible spectrum. Mulliken's net charges have been calculated and compared with the atomic natural charges. The effects of chlorine and cyanide group substituent in benzene ring in the vibrational wavenumbers have been analyzed. NBO analysis is useful to understand the intramolecular hyperconjugative interaction between lone pair Cl and sigma*(C-C) bond orbital.

Modelling organic crystal structures using distributed multipole and polarizability-based model intermolecular potentials

Crystal structure prediction for organic molecules requires both the fast assessment of thousands to millions of crystal structures and the greatest possible accuracy in their relative energies. We describe a crystal lattice simulation program, DMACRYS, emphasizing the features that make it suitable for use in crystal structure prediction for pharmaceutical molecules using accurate anisotropic atom-atom model intermolecular potentials based on the theory of intermolecular forces. DMACRYS can optimize the lattice energy of a crystal, calculate the second derivative properties, and reduce the symmetry of the spacegroup to move away from a transition state. The calculated terahertz frequency k = 0 rigid-body lattice modes and elastic tensor can be used to estimate free energies. The program uses a distributed multipole electrostatic model (Q, t = 00,...,44s) for the electrostatic fields, and can use anisotropic atom-atom repulsion models, damped isotropic dispersion up to R(-10), as well as a range of empirically fitted isotropic exp-6 atom-atom models with different definitions of atomic types. A new feature is that an accurate model for the induction energy contribution to the lattice energy has been implemented that uses atomic anisotropic dipole polarizability models (alpha, t = (10,10)...(11c,11s)) to evaluate the changes in the molecular charge density induced by the electrostatic field within the crystal. It is demonstrated, using the four polymorphs of the pharmaceutical carbamazepine C(15)H(12)N(2)O, that whilst reproducing crystal structures is relatively easy, calculating the polymorphic energy differences to the accuracy of a few kJ mol(-1) required for applications is very demanding of assumptions made in the modelling. Thus DMACRYS enables the comparison of both known and hypothetical crystal structures as an aid to the development of pharmaceuticals and other speciality organic materials, and provides a tool to develop the modelling of the intermolecular forces involved in molecular recognition processes.

Computed crystal energy landscapes for understanding and predicting organic crystal structures and polymorphism

The phenomenon of polymorphism, the ability of a molecule to adopt more than one crystal structure, is a well-established property of crystalline solids. The possible variations in physical properties between polymorphs make the reliable reproduction of a crystalline form essential for all research using organic materials, as well as quality control in manufacture. Thus, the last two decades have seen both an increase in interest in polymorphism and the availability of the computer power needed to make the computational prediction of organic crystal structures a practical possibility. In the past decade, researchers have made considerable improvements in the theoretical basis for calculating the sets of structures that are within the energy range of possible polymorphism, called crystal energy landscapes. It is common to find that a molecule has a wide variety of ways of packing with lattice energy within a few kilojoules per mole of the most stable structure. However, as we develop methods to search for and characterize "all" solid forms, it is also now usual for polymorphs and solvates to be found. Thus, the computed crystal energy landscape reflects and to an increasing extent "predicts" the emerging complexity of the solid state observed for many organic molecules. This Account will discuss the ways in which the calculation of the crystal energy landscape of a molecule can be used as a complementary technique to solid form screening for polymorphs. Current methods can predict the known crystal structure, even under "blind test" conditions, but such successes are generally restricted to those structures that are the most stable over a wide range of thermodynamic conditions. The other low-energy structures can be alternative polymorphs, which have sometimes been found in later experimental studies. Examining the computed structures reveals the various compromises between close packing, hydrogen bonding, and pi-pi stacking that can result in energetically feasible structures. Indeed, we have observed that systems with many almost equi-energetic structures that contain a common interchangeable motif correlate with a tendency to disorder and problems with control of the crystallization product. Thus, contrasting the computed crystal energy landscape with the known crystal structures of a given molecule provides a valuable complement to solid form screening, and the examination of the low-energy structures often leads to a rationalization of the forms found.