Simulation of liquid imidazole using a high-rank quantum topological electrostatic potential

Majeed S. Shaik, Steven Y. Liem, Yongna Yuan, Paul L A Popelier

    Research output: Contribution to journalArticlepeer-review

    174 Downloads (Pure)


    Rigid body molecular dynamics simulations were carried out on pure liquid imidazole at four different temperatures and at 1 atm. Imidazole, which is important both in life science and materials science, is one of the simplest molecules to possess both a lone pair and a π system. These two features are known to benefit from multipolar electrostatics. Here the electrostatic interaction is governed by atomic multipole moments obtained from topologically partitioned ab initio electron densities. The non-electrostatic terms are modeled with Lennard-Jones parameters adjusted to fit the experimental liquid density. All σ values are incrementally increased by one single scaling factor. We report on how the presence of multipolar electrostatics influences the local structure, dynamics and thermodynamics of the liquid compared to electrostatics by atomic point charges. The point charge force field exaggerates the number of π-stacked dimers in the liquid, and underestimates the number of hydrogen-bonded dimers. The effect of the temperature on the local structure of liquid imidazole was analysed using radial and spatial distribution functions. © the Owner Societies 2010.
    Original languageEnglish
    Pages (from-to)15040-15055
    Number of pages15
    JournalPhysical Chemistry Chemical Physics
    Issue number45
    Publication statusPublished - 7 Dec 2010


    Dive into the research topics of 'Simulation of liquid imidazole using a high-rank quantum topological electrostatic potential'. Together they form a unique fingerprint.

    Cite this