Optimization of Crystal Structures of Archetypical Pharmaceutical Compounds: A Plane-Wave DFT‑D Study Using Quantum Espresso

• Published in: Crystal growth & design Vol. 13; no. 5; pp. 2181 - 2189
• Main Authors: Lund, Albert M, Orendt, Anita M, Pagola, Gabriel I, Ferraro, Marta B, Facelli, Julio C
• Format: Journal Article
• Language: English
• Published: Washington,DC American Chemical Society 01.05.2013
• Subjects: Condensed matter: electronic structure, electrical, magnetic, and optical properties; Condensed matter: structure, mechanical and thermal properties; Electron states; Exact sciences and technology; Methods of electronic structure calculations; Physics; Structure of solids and liquids; crystallography; Structure of specific crystalline solids; Electronic structure; Computer software; Force field method; Force field; Theoretical study; Density functional method; Plane waves; Dispersions; Optimization; Genetic algorithms; Crystal structure
• ISSN: 1528-7483; 1528-7505
• DOI: 10.1021/cg4002797

Previously, it was shown that crystal structure prediction based on genetic algorithms (MGAC program) coupled with force field methods could consistently find experimental structures of crystals. However, inaccuracies in the force field potentials often resulted in poor energetic ranking of the experimental structure, limiting the usefulness of the method. In this work, dispersion-corrected density functional theory is employed to improve the accuracy of the energy rankings, using the software package Quantum Espresso. The best choices of running parameters for this application were determined, followed by completion of crystal optimizations on a test set of archetypical pharmaceutical molecules. It is shown here that the variable cell optimization of experimental structures reproduces the experimental structure with high accuracy (RMS < 0.5 Å) for this test set. It is also shown that the use of electronic structure theory based methods greatly improves the energetic ranking of structures produced by MGAC when used with a force field method, such that the experimental match is found with a high degree of accuracy.