Quantitative assessment of electrostatic embedding in density functional theory calculations of biomolecular systems

J. L. Fattebert, R. J. Law, B. Bennion, Edmond Y Lau, E. Schwegler, Felice C Lightstone

Research output: Contribution to journalArticle

6 Citations (Scopus)

Abstract

We evaluate the accuracy of density functional theory quantum calculations of biomolecular subsystems using a simple electrostatic embedding scheme. Our scheme is based on dividing the system of interest into a primary and secondary subsystem. A finite difference discretization of the Kohn-Sham equations is used for the primary subsystem, while its electrostatic environment is modeled with a simple one-electron potential. Force-field atomic partial charges are used to generate smeared Gaussian charge densities and to model the secondary subsystem. We illustrate the utility of this approach with calculations of truncated dipeptide chains. We analyze quantitatively the accuracy of this approach by calculating atomic forces and comparing results with full QM calculations. The impact of the choice made in terminating dangling bonds at the frontier of the QM region is also investigated.

Original languageEnglish (US)
Pages (from-to)2257-2264
Number of pages8
JournalJournal of Chemical Theory and Computation
Volume5
Issue number9
DOIs
StatePublished - Sep 2009
Externally publishedYes

Fingerprint

embedding
Density functional theory
Electrostatics
electrostatics
density functional theory
stopping
field theory (physics)
Dangling bonds
Dipeptides
Charge density
Electrons
electrons

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry
  • Computer Science Applications

Cite this

Quantitative assessment of electrostatic embedding in density functional theory calculations of biomolecular systems. / Fattebert, J. L.; Law, R. J.; Bennion, B.; Lau, Edmond Y; Schwegler, E.; Lightstone, Felice C.

In: Journal of Chemical Theory and Computation, Vol. 5, No. 9, 09.2009, p. 2257-2264.

Research output: Contribution to journalArticle

@article{2928b3393b93409fb180e30687c26670,
title = "Quantitative assessment of electrostatic embedding in density functional theory calculations of biomolecular systems",
abstract = "We evaluate the accuracy of density functional theory quantum calculations of biomolecular subsystems using a simple electrostatic embedding scheme. Our scheme is based on dividing the system of interest into a primary and secondary subsystem. A finite difference discretization of the Kohn-Sham equations is used for the primary subsystem, while its electrostatic environment is modeled with a simple one-electron potential. Force-field atomic partial charges are used to generate smeared Gaussian charge densities and to model the secondary subsystem. We illustrate the utility of this approach with calculations of truncated dipeptide chains. We analyze quantitatively the accuracy of this approach by calculating atomic forces and comparing results with full QM calculations. The impact of the choice made in terminating dangling bonds at the frontier of the QM region is also investigated.",
author = "Fattebert, {J. L.} and Law, {R. J.} and B. Bennion and Lau, {Edmond Y} and E. Schwegler and Lightstone, {Felice C}",
year = "2009",
month = "9",
doi = "10.1021/ct900209y",
language = "English (US)",
volume = "5",
pages = "2257--2264",
journal = "Journal of Chemical Theory and Computation",
issn = "1549-9618",
publisher = "American Chemical Society",
number = "9",

}

TY - JOUR

T1 - Quantitative assessment of electrostatic embedding in density functional theory calculations of biomolecular systems

AU - Fattebert, J. L.

AU - Law, R. J.

AU - Bennion, B.

AU - Lau, Edmond Y

AU - Schwegler, E.

AU - Lightstone, Felice C

PY - 2009/9

Y1 - 2009/9

N2 - We evaluate the accuracy of density functional theory quantum calculations of biomolecular subsystems using a simple electrostatic embedding scheme. Our scheme is based on dividing the system of interest into a primary and secondary subsystem. A finite difference discretization of the Kohn-Sham equations is used for the primary subsystem, while its electrostatic environment is modeled with a simple one-electron potential. Force-field atomic partial charges are used to generate smeared Gaussian charge densities and to model the secondary subsystem. We illustrate the utility of this approach with calculations of truncated dipeptide chains. We analyze quantitatively the accuracy of this approach by calculating atomic forces and comparing results with full QM calculations. The impact of the choice made in terminating dangling bonds at the frontier of the QM region is also investigated.

AB - We evaluate the accuracy of density functional theory quantum calculations of biomolecular subsystems using a simple electrostatic embedding scheme. Our scheme is based on dividing the system of interest into a primary and secondary subsystem. A finite difference discretization of the Kohn-Sham equations is used for the primary subsystem, while its electrostatic environment is modeled with a simple one-electron potential. Force-field atomic partial charges are used to generate smeared Gaussian charge densities and to model the secondary subsystem. We illustrate the utility of this approach with calculations of truncated dipeptide chains. We analyze quantitatively the accuracy of this approach by calculating atomic forces and comparing results with full QM calculations. The impact of the choice made in terminating dangling bonds at the frontier of the QM region is also investigated.

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

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

U2 - 10.1021/ct900209y

DO - 10.1021/ct900209y

M3 - Article

AN - SCOPUS:73349091858

VL - 5

SP - 2257

EP - 2264

JO - Journal of Chemical Theory and Computation

JF - Journal of Chemical Theory and Computation

SN - 1549-9618

IS - 9

ER -