Energy ranking of molecular crystals using density functional theory calculations and an empirical van der Waals correction

被引:255
作者
Neumann, MA
Perrin, MA
机构
[1] SARL, Avant Garde Mat Simulat, F-78100 St Germain En Laye, France
[2] Aventis Pharm SA, Ctr Rech Paris, F-94400 Vitry Sur Seine, France
关键词
D O I
10.1021/jp050121r
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
By combination of high level density functional theory (DFT) calculations with an empirical van der Waals correction, a hybrid method has been designed and parametrized that provides unprecedented accuracy for the structure optimization and the energy ranking of molecular crystals. All DFT calculations are carried out using the VASP program. The van der Waals correction is expressed as the sum over atom-atom pair potentials with each pair potential for two atoms A and B being the product of an asymptotic C-6,C-A,C-B/r(6) term and a damping function d(A,B)(r). Empirical parameters are provided for the elements H, C, N, O, F, Cl, and S. Following Wit and Yang, the C-6 coefficients have been determined by least-squares fitting to molecular C6 coefficients derived by Meath and co-workers from dipole oscillator strength distributions. The damping functions d(A,B)(r) guarantee the crossover from the asymptotic C-6,C-A,C-B/r(6) behavior at large interatomic distances to a constant interaction energy at short distances. The careful parametrization of the damping functions is of crucial importance to obtain the correct balance between the DFT part of the lattice energy and the contribution from the empirical van der Waals correction. The damping functions have been adjusted to yield the best possible agreement between the unit cells of a set of experimental low temperature crystal structures and their counterparts obtained by lattice energy optimization using the hybrid method. On average, the experimental and the calculated unit cell lengths deviate by 1%. To assess the performance of the hybrid method with respect to the lattice energy ranking of molecular crystals, various crystal packings of ethane, ethylene, acetylene, methanol, acetic acid, and urea have been generated with Accelrys' Polymorph Predictor in a first step and optimized with the hybrid method in a second step. In five out of six cases, the experimentally observed low-temperature crystal structure corresponds to the most stable calculated structure.
引用
收藏
页码:15531 / 15541
页数:11
相关论文
共 43 条
[1]   van der Waals interactions in density-functional theory [J].
Andersson, Y ;
Langreth, DC ;
Lundqvist, BI .
PHYSICAL REVIEW LETTERS, 1996, 76 (01) :102-105
[2]   The generation and use of delocalized internal coordinates in geometry optimization [J].
Baker, J ;
Kessi, A ;
Delley, B .
JOURNAL OF CHEMICAL PHYSICS, 1996, 105 (01) :192-212
[3]   VAN DER WAALS VOLUMES + RADII [J].
BONDI, A .
JOURNAL OF PHYSICAL CHEMISTRY, 1964, 68 (03) :441-+
[4]   VALENCE SHELL ABSOLUTE PHOTOABSORPTION OSCILLATOR-STRENGTHS, CONSTRAINED DIPOLE OSCILLATOR STRENGTH DISTRIBUTIONS, AND DIPOLE PROPERTIES FOR CH3NH2, (CH3)(2)NH, AND (CH3)(3)N [J].
BURTON, GR ;
CHAN, WF ;
COOPER, G ;
BRION, CE ;
KUMAR, A ;
MEATH, WJ .
CANADIAN JOURNAL OF CHEMISTRY-REVUE CANADIENNE DE CHIMIE, 1994, 72 (03) :529-546
[5]   Constraint satisfaction in local and gradient susceptibility approximations: Application to a van der Waals density functional [J].
Dobson, JF ;
Dinte, BP .
PHYSICAL REVIEW LETTERS, 1996, 76 (11) :1780-1783
[6]  
DONOHUE J, 1961, ACTA CRYSTALLOGR, V14, P1001
[7]   An approximate DFT method for QM/MM simulations of biological structures and processes [J].
Elstner, M ;
Frauenheim, T ;
Suhai, S .
JOURNAL OF MOLECULAR STRUCTURE-THEOCHEM, 2003, 632 :29-41
[8]   Hydrogen bonding and stacking interactions of nucleic acid base pairs: A density-functional-theory based treatment [J].
Elstner, M ;
Hobza, P ;
Frauenheim, T ;
Suhai, S ;
Kaxiras, E .
JOURNAL OF CHEMICAL PHYSICS, 2001, 114 (12) :5149-5155
[9]   ENERGY FUNCTIONS FOR PEPTIDES AND PROTEINS .1. DERIVATION OF A CONSISTENT FORCE-FIELD INCLUDING HYDROGEN-BOND FROM AMIDE CRYSTALS [J].
HAGLER, AT ;
HULER, E ;
LIFSON, S .
JOURNAL OF THE AMERICAN CHEMICAL SOCIETY, 1974, 96 (17) :5319-5327
[10]   REPRESENTATION OF VANDERWAALS (VDW) INTERACTIONS IN MOLECULAR MECHANICS FORCE-FIELDS - POTENTIAL FORM, COMBINATION RULES, AND VDW PARAMETERS [J].
HALGREN, TA .
JOURNAL OF THE AMERICAN CHEMICAL SOCIETY, 1992, 114 (20) :7827-7843