Improved methods for semiempirical solvation models

Daniel A. Liotard, Gregory D. Hawkins, Gillian C. Lynch, Christopher J. Cramer, Donald G. Truhlar

Research output: Contribution to journalArticle

111 Citations (Scopus)

Abstract

We present improved algorithms for the SMx (x = 1, 1a, 2, 3) solvation models presented previously [see the overview in C. J. Cramer and D. G. Truhlar, J. Comp.‐Aided Mol. Design, 6, 629 (1992)]. These models estimate the free energy of solvation by augmenting a semiempirical Hartree‐Fock calculation on the solute with the generalized Born (GB) model for electric polarization of the solvent and a surface tension term based on solvent‐accessible surface area. This article presents three improvements in the algorithms used to carry out such calculations, namely (1) an analytical accessible surface area algorithm, (2) a more efficient radial integration scheme for the dielectric screening computation in the GB model, and (3) a damping algorithm for updating the GB contribution to the Fock update during the iterations to achieve a self‐consistent field. Improvements (1) and (2) decrease the computer time, and improvement (3) leads to more stable convergence. Improvement (2) removes a small systematic numerical error that was explicitly absorbed into the parameterization in the SMx models. Therefore, we have adjusted the parameters for one of the previous models to yield essentially identical performance as was obtained originally while simultaneously taking advantage of improvement (2). The resulting model is called SM2.1. The fact that we obtain similar results after removing the systematic quadrature bias attests to the robustness of the original parameterization. © 1995 by John Wiley & Sons, Inc.

Original languageEnglish (US)
Pages (from-to)422-440
Number of pages19
JournalJournal of Computational Chemistry
Volume16
Issue number4
DOIs
StatePublished - 1995
Externally publishedYes

Fingerprint

Solvation
Surface area
Parameterization
Model
Stable Convergence
Surface Tension
Quadrature
Free energy
Updating
Screening
Surface tension
Free Energy
Damping
Polarization
Update
Robustness
Iteration
Decrease
Term
Estimate

ASJC Scopus subject areas

  • Chemistry(all)
  • Computational Mathematics

Cite this

Liotard, D. A., Hawkins, G. D., Lynch, G. C., Cramer, C. J., & Truhlar, D. G. (1995). Improved methods for semiempirical solvation models. Journal of Computational Chemistry, 16(4), 422-440. https://doi.org/10.1002/jcc.540160405

Improved methods for semiempirical solvation models. / Liotard, Daniel A.; Hawkins, Gregory D.; Lynch, Gillian C.; Cramer, Christopher J.; Truhlar, Donald G.

In: Journal of Computational Chemistry, Vol. 16, No. 4, 1995, p. 422-440.

Research output: Contribution to journalArticle

Liotard, DA, Hawkins, GD, Lynch, GC, Cramer, CJ & Truhlar, DG 1995, 'Improved methods for semiempirical solvation models', Journal of Computational Chemistry, vol. 16, no. 4, pp. 422-440. https://doi.org/10.1002/jcc.540160405
Liotard, Daniel A. ; Hawkins, Gregory D. ; Lynch, Gillian C. ; Cramer, Christopher J. ; Truhlar, Donald G. / Improved methods for semiempirical solvation models. In: Journal of Computational Chemistry. 1995 ; Vol. 16, No. 4. pp. 422-440.
@article{88f705357a0f4cf9950866f2a9620229,
title = "Improved methods for semiempirical solvation models",
abstract = "We present improved algorithms for the SMx (x = 1, 1a, 2, 3) solvation models presented previously [see the overview in C. J. Cramer and D. G. Truhlar, J. Comp.‐Aided Mol. Design, 6, 629 (1992)]. These models estimate the free energy of solvation by augmenting a semiempirical Hartree‐Fock calculation on the solute with the generalized Born (GB) model for electric polarization of the solvent and a surface tension term based on solvent‐accessible surface area. This article presents three improvements in the algorithms used to carry out such calculations, namely (1) an analytical accessible surface area algorithm, (2) a more efficient radial integration scheme for the dielectric screening computation in the GB model, and (3) a damping algorithm for updating the GB contribution to the Fock update during the iterations to achieve a self‐consistent field. Improvements (1) and (2) decrease the computer time, and improvement (3) leads to more stable convergence. Improvement (2) removes a small systematic numerical error that was explicitly absorbed into the parameterization in the SMx models. Therefore, we have adjusted the parameters for one of the previous models to yield essentially identical performance as was obtained originally while simultaneously taking advantage of improvement (2). The resulting model is called SM2.1. The fact that we obtain similar results after removing the systematic quadrature bias attests to the robustness of the original parameterization. {\circledC} 1995 by John Wiley & Sons, Inc.",
author = "Liotard, {Daniel A.} and Hawkins, {Gregory D.} and Lynch, {Gillian C.} and Cramer, {Christopher J.} and Truhlar, {Donald G.}",
year = "1995",
doi = "10.1002/jcc.540160405",
language = "English (US)",
volume = "16",
pages = "422--440",
journal = "Journal of Computational Chemistry",
issn = "0192-8651",
publisher = "John Wiley and Sons Inc.",
number = "4",

}

TY - JOUR

T1 - Improved methods for semiempirical solvation models

AU - Liotard, Daniel A.

AU - Hawkins, Gregory D.

AU - Lynch, Gillian C.

AU - Cramer, Christopher J.

AU - Truhlar, Donald G.

PY - 1995

Y1 - 1995

N2 - We present improved algorithms for the SMx (x = 1, 1a, 2, 3) solvation models presented previously [see the overview in C. J. Cramer and D. G. Truhlar, J. Comp.‐Aided Mol. Design, 6, 629 (1992)]. These models estimate the free energy of solvation by augmenting a semiempirical Hartree‐Fock calculation on the solute with the generalized Born (GB) model for electric polarization of the solvent and a surface tension term based on solvent‐accessible surface area. This article presents three improvements in the algorithms used to carry out such calculations, namely (1) an analytical accessible surface area algorithm, (2) a more efficient radial integration scheme for the dielectric screening computation in the GB model, and (3) a damping algorithm for updating the GB contribution to the Fock update during the iterations to achieve a self‐consistent field. Improvements (1) and (2) decrease the computer time, and improvement (3) leads to more stable convergence. Improvement (2) removes a small systematic numerical error that was explicitly absorbed into the parameterization in the SMx models. Therefore, we have adjusted the parameters for one of the previous models to yield essentially identical performance as was obtained originally while simultaneously taking advantage of improvement (2). The resulting model is called SM2.1. The fact that we obtain similar results after removing the systematic quadrature bias attests to the robustness of the original parameterization. © 1995 by John Wiley & Sons, Inc.

AB - We present improved algorithms for the SMx (x = 1, 1a, 2, 3) solvation models presented previously [see the overview in C. J. Cramer and D. G. Truhlar, J. Comp.‐Aided Mol. Design, 6, 629 (1992)]. These models estimate the free energy of solvation by augmenting a semiempirical Hartree‐Fock calculation on the solute with the generalized Born (GB) model for electric polarization of the solvent and a surface tension term based on solvent‐accessible surface area. This article presents three improvements in the algorithms used to carry out such calculations, namely (1) an analytical accessible surface area algorithm, (2) a more efficient radial integration scheme for the dielectric screening computation in the GB model, and (3) a damping algorithm for updating the GB contribution to the Fock update during the iterations to achieve a self‐consistent field. Improvements (1) and (2) decrease the computer time, and improvement (3) leads to more stable convergence. Improvement (2) removes a small systematic numerical error that was explicitly absorbed into the parameterization in the SMx models. Therefore, we have adjusted the parameters for one of the previous models to yield essentially identical performance as was obtained originally while simultaneously taking advantage of improvement (2). The resulting model is called SM2.1. The fact that we obtain similar results after removing the systematic quadrature bias attests to the robustness of the original parameterization. © 1995 by John Wiley & Sons, Inc.

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

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

U2 - 10.1002/jcc.540160405

DO - 10.1002/jcc.540160405

M3 - Article

VL - 16

SP - 422

EP - 440

JO - Journal of Computational Chemistry

JF - Journal of Computational Chemistry

SN - 0192-8651

IS - 4

ER -