Modeling the interplay of inter- and intramolecular hydrogen bonding in conformational polymorphs
Modeling the interplay of inter- and intramolecular hydrogen bonding in conformational polymorphs
The predicted stability differences of the conformational polymorphs of oxalyl dihydrazide and ortho-acetamidobenzamide are unrealistically large when the modeling of intermolecular energies is solely based on the isolated-molecule charge density, neglecting charge density polarization. Ab initio calculated crystal electron densities showed qualitative differences depending on the spatial arrangement of molecules in the lattice with the greatest variations observed for polymorphs that differ in the extent of inter- and intramolecular hydrogen bonding. We show that accounting for induction dramatically alters the calculated stability order of the polymorphs and reduces their predicted stability differences to be in better agreement with experiment. Given the challenges in modeling conformational polymorphs with marked differences in hydrogen bonding geometries, we performed an extensive periodic density functional study with a range of exchange-correlation functionals using both atomic and plane wave basis sets. Although such electronic structure methods model the electrostatic and polarization contributions well, the underestimation of dispersion interactions by current exchange-correlation functionals limits their applicability. The use of an empirical dispersion-corrected density functional method consistently reduces the structural deviations between the experimental and energy minimized crystal structures and achieves plausible stability differences. Thus, we have established which types of models may give worthwhile relative energies for crystal structures and other condensed phases of flexible molecules with intra- and intermolecular hydrogen bonding capabilities, advancing the possibility of simulation studies on polymorphic pharmaceuticals.
ab initio calculations, band structure, crystal structure, density functional theory, electron density, hydrogen bonds, polymorphism
244708-[17pp]
Karamertzanis, Panagiotis G.
1b376223-33ec-4e35-87bc-1c4e39959580
Day, Graeme M.
e3be79ba-ad12-4461-b735-74d5c4355636
Welch, Gareth W.A.
9b0e95de-06e7-4afc-b366-b62a6540f4a1
Kendrick, John
94d1a2a2-6c7e-41f8-8ef2-036d5bdfa930
Leusen, Frank J.J.
d41462dd-3804-4380-ad6f-c45e9faeec45
Neumann, Marcus A.
f004825f-59bd-489d-8337-55148a4da3ca
Price, Sarah L.
ab33d469-c548-4a15-918f-b0614ce6129a
24 June 2008
Karamertzanis, Panagiotis G.
1b376223-33ec-4e35-87bc-1c4e39959580
Day, Graeme M.
e3be79ba-ad12-4461-b735-74d5c4355636
Welch, Gareth W.A.
9b0e95de-06e7-4afc-b366-b62a6540f4a1
Kendrick, John
94d1a2a2-6c7e-41f8-8ef2-036d5bdfa930
Leusen, Frank J.J.
d41462dd-3804-4380-ad6f-c45e9faeec45
Neumann, Marcus A.
f004825f-59bd-489d-8337-55148a4da3ca
Price, Sarah L.
ab33d469-c548-4a15-918f-b0614ce6129a
Karamertzanis, Panagiotis G., Day, Graeme M., Welch, Gareth W.A., Kendrick, John, Leusen, Frank J.J., Neumann, Marcus A. and Price, Sarah L.
(2008)
Modeling the interplay of inter- and intramolecular hydrogen bonding in conformational polymorphs.
Journal of Chemical Physics, 128 (24), .
(doi:10.1063/1.2937446).
Abstract
The predicted stability differences of the conformational polymorphs of oxalyl dihydrazide and ortho-acetamidobenzamide are unrealistically large when the modeling of intermolecular energies is solely based on the isolated-molecule charge density, neglecting charge density polarization. Ab initio calculated crystal electron densities showed qualitative differences depending on the spatial arrangement of molecules in the lattice with the greatest variations observed for polymorphs that differ in the extent of inter- and intramolecular hydrogen bonding. We show that accounting for induction dramatically alters the calculated stability order of the polymorphs and reduces their predicted stability differences to be in better agreement with experiment. Given the challenges in modeling conformational polymorphs with marked differences in hydrogen bonding geometries, we performed an extensive periodic density functional study with a range of exchange-correlation functionals using both atomic and plane wave basis sets. Although such electronic structure methods model the electrostatic and polarization contributions well, the underestimation of dispersion interactions by current exchange-correlation functionals limits their applicability. The use of an empirical dispersion-corrected density functional method consistently reduces the structural deviations between the experimental and energy minimized crystal structures and achieves plausible stability differences. Thus, we have established which types of models may give worthwhile relative energies for crystal structures and other condensed phases of flexible molecules with intra- and intermolecular hydrogen bonding capabilities, advancing the possibility of simulation studies on polymorphic pharmaceuticals.
This record has no associated files available for download.
More information
Published date: 24 June 2008
Keywords:
ab initio calculations, band structure, crystal structure, density functional theory, electron density, hydrogen bonds, polymorphism
Organisations:
Organic Chemistry: Synthesis, Catalysis and Flow, Computational Systems Chemistry
Identifiers
Local EPrints ID: 343440
URI: http://eprints.soton.ac.uk/id/eprint/343440
ISSN: 0021-9606
PURE UUID: e1f9ff9a-5040-45ac-b1de-1962c9bc4fac
Catalogue record
Date deposited: 08 Oct 2012 10:12
Last modified: 15 Mar 2024 03:44
Export record
Altmetrics
Contributors
Author:
Panagiotis G. Karamertzanis
Author:
Gareth W.A. Welch
Author:
John Kendrick
Author:
Frank J.J. Leusen
Author:
Marcus A. Neumann
Author:
Sarah L. Price
Download statistics
Downloads from ePrints over the past year. Other digital versions may also be available to download e.g. from the publisher's website.
View more statistics