Repulsion-dispersion parameters for the modelling of organic molecular crystals containing N, O, S and Cl
File(s)c8fd00064f.pdf (1.19 MB)
Published version
Author(s)
Adjiman, CSJ
Pantelides, constantinos
Gatsiou, CA
Type
Journal Article
Abstract
In lattice energy models that combine ab initio and empirical components, it is important to ensure
consistency between these components so that meaningful quantitative results are obtained. A
method for deriving parameters of atom-atom repulsion dispersion potentials for crystals, tailored
to different ab initio models is presented. It is based on minimization of the sum of squared de-
viations between experimental and calculated structures and energies. The solution algorithm
is designed to avoid convergence to local minima in the parameter space by combining a deter-
ministic low-discrepancy sequence for the generation of multiple initial parameter guesses with
an efficient local minimization algorithm. The proposed approach is applied to derive transferable
exp-6 potential parameters suitable for use in conjunction with a distributed multipole electrostatics
model derived from isolated molecule charge densities calculated at the M06/6-31G(d,p) level of
theory. Data for hydrocarbons, azahydrocarbons, oxohydrocarbons, organosulphur compounds
and chlorohydrocarbons are used for the estimation. A good fit is achieved for the new set of
parameters with a mean absolute error in sublimation enthalpies of 4.1 kJ/mol and an average
rmsd
15
of 0.31 Å. The parameters are found to perform well on a separate cross-validation set of
39 compounds.
consistency between these components so that meaningful quantitative results are obtained. A
method for deriving parameters of atom-atom repulsion dispersion potentials for crystals, tailored
to different ab initio models is presented. It is based on minimization of the sum of squared de-
viations between experimental and calculated structures and energies. The solution algorithm
is designed to avoid convergence to local minima in the parameter space by combining a deter-
ministic low-discrepancy sequence for the generation of multiple initial parameter guesses with
an efficient local minimization algorithm. The proposed approach is applied to derive transferable
exp-6 potential parameters suitable for use in conjunction with a distributed multipole electrostatics
model derived from isolated molecule charge densities calculated at the M06/6-31G(d,p) level of
theory. Data for hydrocarbons, azahydrocarbons, oxohydrocarbons, organosulphur compounds
and chlorohydrocarbons are used for the estimation. A good fit is achieved for the new set of
parameters with a mean absolute error in sublimation enthalpies of 4.1 kJ/mol and an average
rmsd
15
of 0.31 Å. The parameters are found to perform well on a separate cross-validation set of
39 compounds.
Date Issued
2018-10-01
Date Acceptance
2018-05-01
Citation
Faraday Discussions, 2018, 211, pp.297-323
ISSN
1359-6640
Publisher
Royal Society of Chemistry
Start Page
297
End Page
323
Journal / Book Title
Faraday Discussions
Volume
211
Copyright Statement
© The Royal Society of Chemistry 2018. This article is licensed under a Creative Commons Attribution 3.0 Unported Licence (https://creativecommons.org/licenses/by/3.0/)
Sponsor
Engineering & Physical Science Research Council (EPSRC)
Engineering & Physical Science Research Council (EPSRC)
Grant Number
EP/J014958/1
EP/J003840/1
Subjects
Science & Technology
Physical Sciences
Chemistry, Physical
Chemistry
NONBONDED POTENTIAL PARAMETERS
INTERMOLECULAR FORCE-FIELD
DISTRIBUTED MULTIPOLE ANALYSIS
STRUCTURE PREDICTION
GLOBAL OPTIMIZATION
ENERGY FUNCTIONS
HYDROGEN-BOND
BLIND TEST
ATOM
PROTEINS
0306 Physical Chemistry (incl. Structural)
0904 Chemical Engineering
Chemical Physics
Publication Status
Published
Date Publish Online
2018-07-12