Charge Model 5: An Extension of Hirshfeld Population Analysis for the Accurate Description of Molecular Interactions in Gaseous and Condensed Phases

Citations Over TimeTop 1% of 2012 papers

Abstract

We propose a novel approach to deriving partial atomic charges from population analysis. The new model, called Charge Model 5 (CM5), yields class IV partial atomic charges by mapping from those obtained by Hirshfeld population analysis of density functional electronic charge distributions. The CM5 model utilizes a single set of parameters derived by fitting to reference values of the gas-phase dipole moments of 614 molecular structures. An additional test set (not included in the CM5 parametrization) contained 107 singly charged ions with nonzero dipole moments, calculated from the accurate electronic charge density, with respect to the center of nuclear charges. The CM5 model is applicable to any charged or uncharged molecule composed of any element of the periodic table in the gas phase or in solution. The CM5 model predicts dipole moments for the tested molecules that are more accurate on average than those from the original Hirshfeld method or from many other popular schemes including atomic polar tensor and Löwdin, Mulliken, and natural population analyses. In addition, the CM5 charge model is essentially independent of a basis set. It can be used with larger basis sets, and thereby this model significantly improves on our previous charge models CMx (x = 1-4 or 4M) and other methods that are prone to basis set sensitivity. CM5 partial atomic charges are less conformationally dependent than those derived from electrostatic potentials. The CM5 model does not suffer from ill conditioning for buried atoms in larger molecules, as electrostatic fitting schemes sometimes do. The CM5 model can be used with any level of electronic structure theory (Hartree-Fock, post-Hartree-Fock, and other wave function correlated methods or density functional theory) as long as an accurate electronic charge distribution and a Hirshfeld analysis can be computed for that level of theory.

Related Papers

• → Theoretical study on charge distribution in cetylpyridinium cationic surfactant(2021)15 cited
• → Molecular Dynamics Simulation for the Calibration of the OPLS Force Field Using DFT Derived Partial Charges for the Screening of Alkyl Phosphate Ligands by Studying Structure, Dynamics, and Thermodynamics(2017)23 cited
• → Rapid Estimation of Basis Set Error and Correlation Energy Based on Mulliken Charges and Mulliken Matrix with the Small 6-31g* Basis Set(2005)5 cited
• → Effect of charge distribution on electrostatic chromophore?protein interactions inBacteriorhodopsin(1997)5 cited
• → Synthesis, Molecular Structure, HOMO-LUMO and Spectroscopic Investigation of(E)-1-(2,4-Dichloro-5-fluorophenyl)-3-(2,6-dichlorophenyl)prop-2-en-1-one:A DFT Based Computational Exploration(2020)1 cited