Atomic level anisotropy in the electrostatic modeling of lone pairs for a polarizable force field based on the classical drude oscillator

Edward Harder, Victor M. Anisimov, Igor Vorobyov, Pedro E.M. Lopes, Sergei Y. Noskov, Alexander D. MacKerell, Benoît Roux

Research output: Contribution to journalArticle

108 Citations (Scopus)

Abstract

Electron pairs in the valence shell of an atom that do not participate in the bonding of a molecule ("lone pairs") give rise to a concentrated electron density away from the atom center. To account for the asymmetry in the electron charge density that arises from lone pairs, an electrostatic model is developed that is parametrically anisotropic at the atomic level. The model uses virtual interaction sites with partial charges that are associated but not coincident with the nuclei. In addition, the model incorporates anisotropic atomic polarizabilities. The protocol previously outlined in Anisimov et al. [J. Chem. Theory Comput. 2005, 1, 153] for parametrizing the electrostatic potential energy of a polarizable force field using classical Drude oscillators is extended to incorporate additional lone pair parameters. To probe the electrostatic environment around the lone pairs, the static (molecule alone) and perturbed (molecule in the presence of a test charge) electrostatic potential (ESP) are evaluated and compared to high level quantum mechanical (QM) electronic structure calculations. The parametrization of the virtual sites relies on data from the QM static ESP. The contribution to the perturbed ESP from the electronic polarization of the molecule is used to resolve the components of the atomic polarizability tensor. The model is tested in the case of four molecules: methanol, acetone, methylamine, and pyridine. Interaction energies with water and sodium are used to assess the accuracy of the model. The results are compared with simpler models placing all the charge on the nuclei as well as using only isotropic atomic polarizabilities. Analysis shows that the addition of virtual sites reduces the average error relative to the QM calculations. In contrast to models with atom centered charges, the virtual site models correctly predict the minimum energy conformation for acetone and methanol, with water, to be closely coordinated with the lone pair direction. Furthermore, addition of anisotropic atomic polarizabilities to the virtual site model allows for precise fitting to the local perturbed QM ESP.

Original languageEnglish (US)
Pages (from-to)1587-1597
Number of pages11
JournalJournal of Chemical Theory and Computation
Volume2
Issue number6
DOIs
StatePublished - Dec 1 2006
Externally publishedYes

Fingerprint

field theory (physics)
Electrostatics
Anisotropy
oscillators
electrostatics
anisotropy
Molecules
molecules
Acetone
Atoms
acetone
Methanol
methyl alcohol
atoms
electrostatic charge
nuclei
Electrons
Water
Potential energy
Charge density

ASJC Scopus subject areas

  • Computer Science Applications
  • Physical and Theoretical Chemistry

Cite this

Atomic level anisotropy in the electrostatic modeling of lone pairs for a polarizable force field based on the classical drude oscillator. / Harder, Edward; Anisimov, Victor M.; Vorobyov, Igor; Lopes, Pedro E.M.; Noskov, Sergei Y.; MacKerell, Alexander D.; Roux, Benoît.

In: Journal of Chemical Theory and Computation, Vol. 2, No. 6, 01.12.2006, p. 1587-1597.

Research output: Contribution to journalArticle

Harder, Edward ; Anisimov, Victor M. ; Vorobyov, Igor ; Lopes, Pedro E.M. ; Noskov, Sergei Y. ; MacKerell, Alexander D. ; Roux, Benoît. / Atomic level anisotropy in the electrostatic modeling of lone pairs for a polarizable force field based on the classical drude oscillator. In: Journal of Chemical Theory and Computation. 2006 ; Vol. 2, No. 6. pp. 1587-1597.
@article{4b04f85e85094a0e9100e2567e4a3a73,
title = "Atomic level anisotropy in the electrostatic modeling of lone pairs for a polarizable force field based on the classical drude oscillator",
abstract = "Electron pairs in the valence shell of an atom that do not participate in the bonding of a molecule ({"}lone pairs{"}) give rise to a concentrated electron density away from the atom center. To account for the asymmetry in the electron charge density that arises from lone pairs, an electrostatic model is developed that is parametrically anisotropic at the atomic level. The model uses virtual interaction sites with partial charges that are associated but not coincident with the nuclei. In addition, the model incorporates anisotropic atomic polarizabilities. The protocol previously outlined in Anisimov et al. [J. Chem. Theory Comput. 2005, 1, 153] for parametrizing the electrostatic potential energy of a polarizable force field using classical Drude oscillators is extended to incorporate additional lone pair parameters. To probe the electrostatic environment around the lone pairs, the static (molecule alone) and perturbed (molecule in the presence of a test charge) electrostatic potential (ESP) are evaluated and compared to high level quantum mechanical (QM) electronic structure calculations. The parametrization of the virtual sites relies on data from the QM static ESP. The contribution to the perturbed ESP from the electronic polarization of the molecule is used to resolve the components of the atomic polarizability tensor. The model is tested in the case of four molecules: methanol, acetone, methylamine, and pyridine. Interaction energies with water and sodium are used to assess the accuracy of the model. The results are compared with simpler models placing all the charge on the nuclei as well as using only isotropic atomic polarizabilities. Analysis shows that the addition of virtual sites reduces the average error relative to the QM calculations. In contrast to models with atom centered charges, the virtual site models correctly predict the minimum energy conformation for acetone and methanol, with water, to be closely coordinated with the lone pair direction. Furthermore, addition of anisotropic atomic polarizabilities to the virtual site model allows for precise fitting to the local perturbed QM ESP.",
author = "Edward Harder and Anisimov, {Victor M.} and Igor Vorobyov and Lopes, {Pedro E.M.} and Noskov, {Sergei Y.} and MacKerell, {Alexander D.} and Beno{\^i}t Roux",
year = "2006",
month = "12",
day = "1",
doi = "10.1021/ct600180x",
language = "English (US)",
volume = "2",
pages = "1587--1597",
journal = "Journal of Chemical Theory and Computation",
issn = "1549-9618",
publisher = "American Chemical Society",
number = "6",

}

TY - JOUR

T1 - Atomic level anisotropy in the electrostatic modeling of lone pairs for a polarizable force field based on the classical drude oscillator

AU - Harder, Edward

AU - Anisimov, Victor M.

AU - Vorobyov, Igor

AU - Lopes, Pedro E.M.

AU - Noskov, Sergei Y.

AU - MacKerell, Alexander D.

AU - Roux, Benoît

PY - 2006/12/1

Y1 - 2006/12/1

N2 - Electron pairs in the valence shell of an atom that do not participate in the bonding of a molecule ("lone pairs") give rise to a concentrated electron density away from the atom center. To account for the asymmetry in the electron charge density that arises from lone pairs, an electrostatic model is developed that is parametrically anisotropic at the atomic level. The model uses virtual interaction sites with partial charges that are associated but not coincident with the nuclei. In addition, the model incorporates anisotropic atomic polarizabilities. The protocol previously outlined in Anisimov et al. [J. Chem. Theory Comput. 2005, 1, 153] for parametrizing the electrostatic potential energy of a polarizable force field using classical Drude oscillators is extended to incorporate additional lone pair parameters. To probe the electrostatic environment around the lone pairs, the static (molecule alone) and perturbed (molecule in the presence of a test charge) electrostatic potential (ESP) are evaluated and compared to high level quantum mechanical (QM) electronic structure calculations. The parametrization of the virtual sites relies on data from the QM static ESP. The contribution to the perturbed ESP from the electronic polarization of the molecule is used to resolve the components of the atomic polarizability tensor. The model is tested in the case of four molecules: methanol, acetone, methylamine, and pyridine. Interaction energies with water and sodium are used to assess the accuracy of the model. The results are compared with simpler models placing all the charge on the nuclei as well as using only isotropic atomic polarizabilities. Analysis shows that the addition of virtual sites reduces the average error relative to the QM calculations. In contrast to models with atom centered charges, the virtual site models correctly predict the minimum energy conformation for acetone and methanol, with water, to be closely coordinated with the lone pair direction. Furthermore, addition of anisotropic atomic polarizabilities to the virtual site model allows for precise fitting to the local perturbed QM ESP.

AB - Electron pairs in the valence shell of an atom that do not participate in the bonding of a molecule ("lone pairs") give rise to a concentrated electron density away from the atom center. To account for the asymmetry in the electron charge density that arises from lone pairs, an electrostatic model is developed that is parametrically anisotropic at the atomic level. The model uses virtual interaction sites with partial charges that are associated but not coincident with the nuclei. In addition, the model incorporates anisotropic atomic polarizabilities. The protocol previously outlined in Anisimov et al. [J. Chem. Theory Comput. 2005, 1, 153] for parametrizing the electrostatic potential energy of a polarizable force field using classical Drude oscillators is extended to incorporate additional lone pair parameters. To probe the electrostatic environment around the lone pairs, the static (molecule alone) and perturbed (molecule in the presence of a test charge) electrostatic potential (ESP) are evaluated and compared to high level quantum mechanical (QM) electronic structure calculations. The parametrization of the virtual sites relies on data from the QM static ESP. The contribution to the perturbed ESP from the electronic polarization of the molecule is used to resolve the components of the atomic polarizability tensor. The model is tested in the case of four molecules: methanol, acetone, methylamine, and pyridine. Interaction energies with water and sodium are used to assess the accuracy of the model. The results are compared with simpler models placing all the charge on the nuclei as well as using only isotropic atomic polarizabilities. Analysis shows that the addition of virtual sites reduces the average error relative to the QM calculations. In contrast to models with atom centered charges, the virtual site models correctly predict the minimum energy conformation for acetone and methanol, with water, to be closely coordinated with the lone pair direction. Furthermore, addition of anisotropic atomic polarizabilities to the virtual site model allows for precise fitting to the local perturbed QM ESP.

UR - http://www.scopus.com/inward/record.url?scp=33846185169&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=33846185169&partnerID=8YFLogxK

U2 - 10.1021/ct600180x

DO - 10.1021/ct600180x

M3 - Article

AN - SCOPUS:33846185169

VL - 2

SP - 1587

EP - 1597

JO - Journal of Chemical Theory and Computation

JF - Journal of Chemical Theory and Computation

SN - 1549-9618

IS - 6

ER -