Electrostatically embedded many-body method for dipole moments, partial atomic charges, and charge transfer

2012 ◽  
Vol 14 (21) ◽  
pp. 7669 ◽  
Author(s):  
Hannah R. Leverentz ◽  
Katie A. Maerzke ◽  
Samuel J. Keasler ◽  
J. Ilja Siepmann ◽  
Donald G. Truhlar
2020 ◽  
Author(s):  
Bowen Han ◽  
Christine Isborn ◽  
Liang Shi

Partial atomic charges provide an intuitive and efficient way to describe the charge distribution and the resulting intermolecular electrostatic interactions in liquid water. Many charge models exist and it is unclear which model provides the best assignment of partial atomic charges in response to the local molecular environment. In this work, we systematically scrutinize various electronic structure methods and charge models (Mulliken, Natural Population Analysis, CHelpG, RESP, Hirshfeld, Iterative Hirshfeld, and Bader) by evaluating their performance in predicting the dipole moments of isolated water, water clusters, and liquid water as well as charge transfer in the water dimer and liquid water. Although none of the seven charge models is capable of fully capturing the dipole moment increase from isolated water (1.85 D) to liquid water (about 2.9 D), the Iterative Hirshfeld method performs best for liquid water, reproducing its experimental average molecular dipole moment, yielding a reasonable amount of intermolecular charge transfer, and showing modest sensitivity to the local water environment. The performance of the charge model is dependent on the choice of the density functional and the quantum treatment of the environment. The computed molecular dipole moment of water generally increases with the percentage of the exact Hartree-Fock exchange in the functional, whereas the amount of charge transfer between molecules decreases. For liquid water, including two full solvation shells of surrounding water molecules (within about 5.5 A of the central water) in the quantum-chemical calculation converges the charges of the central water molecule. Our final pragmatic quantum-chemical charge assigning protocol for liquid water is the Iterative Hirshfeld method with M06-HF/aug-cc-pVDZ and a quantum region cutoff radius of 5.5 A.<br>


2020 ◽  
Author(s):  
Bowen Han ◽  
Christine Isborn ◽  
Liang Shi

Partial atomic charges provide an intuitive and efficient way to describe the charge distribution and the resulting intermolecular electrostatic interactions in liquid water. Many charge models exist and it is unclear which model provides the best assignment of partial atomic charges in response to the local molecular environment. In this work, we systematically scrutinize various electronic structure methods and charge models (Mulliken, Natural Population Analysis, CHelpG, RESP, Hirshfeld, Iterative Hirshfeld, and Bader) by evaluating their performance in predicting the dipole moments of isolated water, water clusters, and liquid water as well as charge transfer in the water dimer and liquid water. Although none of the seven charge models is capable of fully capturing the dipole moment increase from isolated water (1.85 D) to liquid water (about 2.9 D), the Iterative Hirshfeld method performs best for liquid water, reproducing its experimental average molecular dipole moment, yielding a reasonable amount of intermolecular charge transfer, and showing modest sensitivity to the local water environment. The performance of the charge model is dependent on the choice of the density functional and the quantum treatment of the environment. The computed molecular dipole moment of water generally increases with the percentage of the exact Hartree-Fock exchange in the functional, whereas the amount of charge transfer between molecules decreases. For liquid water, including two full solvation shells of surrounding water molecules (within about 5.5 A of the central water) in the quantum-chemical calculation converges the charges of the central water molecule. Our final pragmatic quantum-chemical charge assigning protocol for liquid water is the Iterative Hirshfeld method with M06-HF/aug-cc-pVDZ and a quantum region cutoff radius of 5.5 A.<br>


2019 ◽  
Author(s):  
Alex Zhou ◽  
Michael Schauperl ◽  
Paul Nerenberg

<p>The accuracy of classical molecular mechanics (MM) force fields used for condensed phase molecular simulations depends strongly on the accuracy of modeling nonbonded interactions between atoms, such as electrostatic interactions. Some popular fixed-charge MM force fields use partial atomic charges derived from gas phase electronic structure calculations using the Hartree-Fock method with the relatively small 6-31G* basis set (HF/6-31G*). It is generally believed that HF/6-31G* generates fortuitously overpolarized electron distributions, as would be expected in the higher dielectric environment of the condensed phase. Using a benchmark set of 47 molecules we show that HF/6-31G* overpolarizes molecules by just under 10% on average with respect to experimental gas phase dipole moments. The overpolarization of this method/basis set combination varies significantly though and, in some cases, even leads to molecular dipole moments that are lower than experimental gas phase measurements. We further demonstrate that using computationally inexpensive density functional theory (DFT) methods, together with appropriate augmented basis sets and a continuum solvent model, can yield molecular dipole moments that are both more strongly and more uniformly overpolarized. These data suggest that these methods – or ones similar to them – should be adopted for the derivation of accurate partial atomic charges for next-generation MM force fields.<br></p>


2019 ◽  
Author(s):  
Alex Zhou ◽  
Michael Schauperl ◽  
Paul Nerenberg

<p>The accuracy of classical molecular mechanics (MM) force fields used for condensed phase molecular simulations depends strongly on the accuracy of modeling nonbonded interactions between atoms, such as electrostatic interactions. Some popular fixed-charge MM force fields use partial atomic charges derived from gas phase electronic structure calculations using the Hartree-Fock method with the relatively small 6-31G* basis set (HF/6-31G*). It is generally believed that HF/6-31G* generates fortuitously overpolarized electron distributions, as would be expected in the higher dielectric environment of the condensed phase. Using a benchmark set of 47 molecules we show that HF/6-31G* overpolarizes molecules by just under 10% on average with respect to experimental gas phase dipole moments. The overpolarization of this method/basis set combination varies significantly though and, in some cases, even leads to molecular dipole moments that are lower than experimental gas phase measurements. We further demonstrate that using computationally inexpensive density functional theory (DFT) methods, together with appropriate augmented basis sets and a continuum solvent model, can yield molecular dipole moments that are both more strongly and more uniformly overpolarized. These data suggest that these methods – or ones similar to them – should be adopted for the derivation of accurate partial atomic charges for next-generation MM force fields.<br></p>


2019 ◽  
Author(s):  
Alex Zhou ◽  
Michael Schauperl ◽  
Paul Nerenberg

<div><div><div><p>The accuracy of classical molecular mechanics (MM) force fields used for condensed phase molecular simulations depends strongly on the accuracy of modeling nonbonded interactions between atoms, such as electrostatic interactions. Some popular fixed-charge MM force fields use partial atomic charges derived from gas phase electronic structure calculations using the Hartree-Fock method with the relatively small 6-31G* basis set (HF/6-31G*). It is generally believed that HF/6-31G* generates fortuitously overpolarized electron distributions, as would be expected in the higher dielectric environment of the condensed phase. Using a benchmark set of 47 molecules, we show that HF/6-31G* does not uniformly overpolarize molecules and in some cases even leads to molecular dipole moments that are lower than experimental gas phase measurements. We further demonstrate that using computationally inexpensive density functional theory (DFT) methods, together with appropriate augmented basis sets and a continuum solvent model, can yield molecular dipole moments that are both more strongly and more uniformly overpolarized. These data suggest that these methods – or ones similar to them – should be adopted for the derivation of accurate partial atomic charges for next-generation MM force fields.</p></div></div></div>


2019 ◽  
Author(s):  
Brandon B. Bizzarro ◽  
Colin K. Egan ◽  
Francesco Paesani

<div> <div> <div> <p>Interaction energies of halide-water dimers, X<sup>-</sup>(H<sub>2</sub>O), and trimers, X<sup>-</sup>(H<sub>2</sub>O)<sub>2</sub>, with X = F, Cl, Br, and I, are investigated using various many-body models and exchange-correlation functionals selected across the hierarchy of density functional theory (DFT) approximations. Analysis of the results obtained with the many-body models demonstrates the need to capture important short-range interactions in the regime of large inter-molecular orbital overlap, such as charge transfer and charge penetration. Failure to reproduce these effects can lead to large deviations relative to reference data calculated at the coupled cluster level of theory. Decompositions of interaction energies carried out with the absolutely localized molecular orbital energy decomposition analysis (ALMO-EDA) method demonstrate that permanent and inductive electrostatic energies are accurately reproduced by all classes of XC functionals (from generalized gradient corrected (GGA) to hybrid and range-separated functionals), while significant variance is found for charge transfer energies predicted by different XC functionals. Since GGA and hybrid XC functionals predict the most and least attractive charge transfer energies, respectively, the large variance is likely due to the delocalization error. In this scenario, the hybrid XC functionals are then expected to provide the most accurate charge transfer energies. The sum of Pauli repulsion and dispersion energies are the most varied among the XC functionals, but it is found that a correspondence between the interaction energy and the ALMO EDA total frozen energy may be used to determine accurate estimates for these contributions. </p> </div> </div> </div>


1987 ◽  
Vol 105 ◽  
Author(s):  
Hisham Z. Massoud

AbstractThe magnitude of the dipole moment at the Si-SiO2 interface resulting from partial charge transfer that takes place upon the formation of interface bonds has been calculated. The charge transfer occurs because of the difference in electronegativity between silicon atoms and SiO2 molecules which are present across the interface. Results obtained for (100) and (111) silicon substrates indicate that the magnitude of the interface dipole moment is dependent on substrate orientation and the interface chemistry. Dipole moments at the Si-SiO2 and gate-SiO2 interfaces should be included in the definition of the flatband voltage VFB of MOS structures. CV-based measurements of the metal-semiconductor workfunction difference φms on (100) and (111) silicon oxidized in dry oxygen and metallized with Al agree with the predictions of this model. Other types of interface dipoles and their processing dependence are briefly discussed.


Sign in / Sign up

Export Citation Format

Share Document