Skip to main content


  • Oral presentation
  • Open Access

Assessment of a variety of dispersion-corrected density functional theory calculations used in molecular crystal structure prediction

Journal of Cheminformatics20124 (Suppl 1) :O15

  • Published:


  • Density Functional Theory
  • Density Functional Theory Calculation
  • Cambridge Crystallographic Data
  • Neutral Molecule
  • Lattice Energy

For a long time, the prediction of molecular crystal structures from first principles has been one of the most challenging computational problems. Because of the large number of degrees of freedom, methods with low computational effort have to be used. The price for this efficiency is, however, a loss of accuracy. In recent years, important progress has been made in the development of dispersion corrections [13], which raise the accuracy of density functional theory (DFT) calculations to a much higher level.

Based on these developments, the software package GRACE [4] was the first program ever to correctly predict all four crystal structures at the 4th crystal structure prediction blind test, organized by G. M. Day and the Cambridge Crystallographic Data Centre in 2007 [5]. GRACE uses dispersion-corrected density functional theory calculations (DFT-D) to generate reference data to which a tailor-made force field (TMFF) is fitted for each respective molecule [6]. Crystal structures are generated with a Monte-Carlo parallel-tempering algorithm that uses the TMFF for the evaluation of lattice energies. Only the most stable crystal structures according to the TMFF are then re-optimized and re-ranked using DFT-D. The DFT calculations are performed using an interface to the two common ab-inito programs VASP [7] and QuantumESPRESSO [8].

In the most recent, 5th blind test, organized in 2010, the excellent agreement between experimental and predicted structures for small, neutral molecules was confirmed. However, it turned out that the crystal structure prediction of molecular salts and hydrates is still a challenging task. The inability of DFT without exact exchange to describe anions properly, as recently published by Jensen [9], suggests that the inclusion of the Hartree-Fock exchange should improve the accuracy for molecular salts and hydrates. Therefore, we will present an assessment of a variety of current density functionals, in which the dispersion-correction parameters were fitted with respect to a minimization of the deviation between experimental and computed crystal structures.

Authors’ Affiliations

Avant-garde Materials Simulation, Merzhauser Straße 177, 79100 Freiburg im Breisgau, Germany


  1. Grimme S: . J Comput Chem. 2004, 25: 1463-10.1002/jcc.20078.View ArticleGoogle Scholar
  2. Neumann MA, Perrin M-A: . J Phys Chem B. 2005, 109: 15531-10.1021/jp050121r.View ArticleGoogle Scholar
  3. Grimme S, et al: . J Chem Phys. 2010, 132: 154104-10.1063/1.3382344.View ArticleGoogle Scholar
  4. . []
  5. Day GM, et al: . Acta Crystallogr B. 2009, 65: 107-10.1107/S0108768109004066.View ArticleGoogle Scholar
  6. Neumann MA: . J Phys Chem B. 2008, 112: 9810-10.1021/jp710575h.View ArticleGoogle Scholar
  7. Kresse G, Hafner J: . Phys Rev B. 1993, 47: 558-10.1103/PhysRevB.47.558.View ArticleGoogle Scholar
  8. Giannozzi P, et al: . J Phys Condens Matter. 2009, 21: 395502-10.1088/0953-8984/21/39/395502.View ArticleGoogle Scholar
  9. Jensen F: . J Chem Theory Comput. 2010, 6: 2726-10.1021/ct1003324.View ArticleGoogle Scholar


© Doser and Neumann; licensee BioMed Central Ltd. 2012

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.