Affiliation:
1. State Key Laboratory of Precision Spectroscopy, Department of Physics, East China Normal University, Shanghai 200062, China
2. Department of Chemistry, New York University, New York, NY 10003, USA
Abstract
It is known that electrostatic potential (ESP) derived atomic charges for molecular systems suffer from rank deficiency in least-square fitting. In this paper, we studied numerical properties of atomic charges by fitting to ESP and electric field (EF) for a series of molecular systems, ranging from small molecules to peptides. Our study shows that although atomic charges of certain atoms in a molecule may differ a lot resulting from different fitting schemes, the effect on some observed physical properties of the molecules, such as solvation free energy and radial distributions of the solvent, are negligible. For all the molecules studied in this paper, it is shown that by incorporating EF terms in the least-square fitting, the rank of the least-squares matrix is increased, and the full rank is achieved when only EF terms are included. The current study demonstrates that by using EF, instead of the widely used ESP, in least-square fitting of atomic charges, one can obtain numerically more stable and better defined atomic charges in molecules. Such derived atomic charges may have more desirable properties and better numerical stabilities in studying detailed molecular processes such as in molecular dynamics simulation of macromolecules.
Publisher
World Scientific Pub Co Pte Lt
Subject
Computational Theory and Mathematics,Physical and Theoretical Chemistry,Computer Science Applications
Cited by
2 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献