TY - JOUR
T1 - Potential parameterization from proxy systems
AU - Zhu, Wuming
AU - Runge, Keith
AU - Trickey, S. B.
N1 - Funding Information:
This work was supported in part by NSF ITR (medium) grant DMR-0325553 and, prior to that, by NSF KDI grant DMR-9980015. KR acknowledges informative correspondence with S.T. Bromley.
PY - 2006/10
Y1 - 2006/10
N2 - We give an analysis and critique of a procedure for sequential multi-scale parameterization of so-called ab initio potentials. Such potentials are based upon a small selection of reference data generated at relatively low levels refinement of quantum mechanical (QM) methodology. Two closely related potentials for silica, BKS and TTAM, are explored in some depth because both have found frequent use in the literature of molecular dynamics simulation. The use of a single small molecular cluster as a proxy for the behavior of larger silica structures is shown to be perilous. Choices of ionic charge, basis set size, inclusion of electron correlation at some prescribed level of approximation, and parameterization data are each shown to have a substantial impact on the final potential parameterization. The analysis and critique provide the logical foundation for the different approach to the problem of parameterization and use of classical potentials that is treated in detail elsewhere in this collection. The use of extended model systems, force data calculated from QM methods, and a focused (limited) range of application are shown to be essential for a consistent (and ultimately, predictive) embedding approach to concurrent multi-scale materials simulation.
AB - We give an analysis and critique of a procedure for sequential multi-scale parameterization of so-called ab initio potentials. Such potentials are based upon a small selection of reference data generated at relatively low levels refinement of quantum mechanical (QM) methodology. Two closely related potentials for silica, BKS and TTAM, are explored in some depth because both have found frequent use in the literature of molecular dynamics simulation. The use of a single small molecular cluster as a proxy for the behavior of larger silica structures is shown to be perilous. Choices of ionic charge, basis set size, inclusion of electron correlation at some prescribed level of approximation, and parameterization data are each shown to have a substantial impact on the final potential parameterization. The analysis and critique provide the logical foundation for the different approach to the problem of parameterization and use of classical potentials that is treated in detail elsewhere in this collection. The use of extended model systems, force data calculated from QM methods, and a focused (limited) range of application are shown to be essential for a consistent (and ultimately, predictive) embedding approach to concurrent multi-scale materials simulation.
KW - Force fields
KW - Interatomic potentials
KW - Potential parameters
KW - ab initio potentials
UR - http://www.scopus.com/inward/record.url?scp=33749003135&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=33749003135&partnerID=8YFLogxK
U2 - 10.1007/s10820-006-9016-y
DO - 10.1007/s10820-006-9016-y
M3 - Article
AN - SCOPUS:33749003135
SN - 0928-1045
VL - 13
SP - 75
EP - 84
JO - Journal of Computer-Aided Materials Design
JF - Journal of Computer-Aided Materials Design
IS - 1-3
ER -