Atomic charges of the water molecule and the water dimer (original) (raw)
Related papers
Molecular Theories and Simulation of Ions and Polar Molecules in Water
The Journal of Physical Chemistry A, 1998
Recent developments in molecular theories and simulation of ions and polar molecules in water are reviewed. The hydration of imidazole and imidazolium solutes is used to exemplify the theoretical issues. The treatment of long-ranged electrostatic interactions in simulations is discussed extensively. It is argued that the Ewald approach is an easy way to get correct hydration free energies in the thermodynamic limit from molecular calculations; and that molecular simulations with Ewald interactions and periodic boundary conditions can also be more efficient than many common alternatives. The Ewald treatment permits a conclusive extrapolation to infinite system size. Accurate results for well-defined models have permitted careful testing of simple theories of electrostatic hydration free energies, such as dielectric continuum models. The picture that emerges from such testing is that the most prominent failings of the simplest theories are associated with solvent proton conformations that lead to non-gaussian fluctuations of electrostatic potentials. Thus, the most favorable cases for second-order perturbation theories are monoatomic positive ions. For polar and anionic solutes, continuum or gaussian theories are less accurate. The appreciation of the specific deficiencies of those simple models have led to new concepts, multistate gaussian and quasi-chemical theories, that address the cases for which the simpler theories fail. It is argued that, relative to direct dielectric continuum treatments, the quasi-chemical theories provide a better theoretical organization for the computational study of the electronic structure of solution species.
Polarizable Atomic Multipole Water Model for Molecular Mechanics Simulation
The Journal of Physical Chemistry B, 2003
A new classical empirical potential is proposed for water. The model uses a polarizable atomic multipole description of electrostatic interactions. Multipoles through the quadrupole are assigned to each atomic center based on a distributed multipole analysis (DMA) derived from large basis set molecular orbital calculations on the water monomer. Polarization is treated via self-consistent induced atomic dipoles. A modified version of Thole's interaction model is used to damp induction at short range. Repulsion-dispersion (vdW) effects are computed from a buffered 14-7 potential. In a departure from most current water potentials, we find that significant vdW parameters are necessary on hydrogen as well as oxygen. The new potential is fully flexible and has been tested versus a variety of experimental data and quantum calculations for small clusters, liquid water, and ice. Overall, excellent agreement with experimental and high level ab initio results is obtained for numerous properties, including cluster structures and energetics and bulk thermodynamic and structural measures. The parametrization scheme described here is easily extended to other molecular systems, and the resulting water potential should provide a useful explicit solvent model for organic solutes and biopolymer modeling.
Polarizable Water Potential Derived from a Model Electron Density
Journal of Chemical Theory and Computation, 2021
A new empirical potential for efficient, large scale molecular dynamics simulation of water is presented. The HIPPO (Hydrogen-like Intermolecular Polarizable POtential) force field is based upon the model electron density of a hydrogen-like atom. This framework is used to derive and parametrize individual terms describing charge penetration damped permanent electrostatics, damped polarization, charge transfer, anisotropic Pauli repulsion, and damped dispersion interactions. Initial parameter values were fit to Symmetry Adapted Perturbation Theory (SAPT) energy components for ten water dimer configurations, as well as the radial and angular dependence of the canonical dimer. The SAPTbased parameters were then systematically refined to extend the treatment to water bulk phases. The final HIPPO water model provides a balanced representation of a wide variety of properties of gas phase clusters, liquid water, and ice polymorphs, across a range of temperatures and pressures. This water potential yields a rationalization of water structure, dynamics, and thermodynamics explicitly correlated with an ab initio energy decomposition, while providing a level of accuracy comparable or superior to previous polarizable atomic multipole force fields. The HIPPO water model serves as a cornerstone around which similarly detailed physicsbased models can be developed for additional molecular species.
Chemical Physics, 2002
Electric field gradients at the oxygen and hydrogen nuclei of water have been calculated using high level ab initio methods. Systematic studies of basis set truncation errors have been carried out at the Hartree-Fock and coupled cluster singles and doubles ͑CCSD͒ levels using extended correlation consistent basis sets with up to 398 basis functions. Correlation effects are investigated using a hierarchy of correlation methods extending up to the approximate inclusion of triples excitations by means of the CCSD͑T͒ method. Rovibrational effects have been calculated combining accurate ab initio electric field gradient data and accurate experimental force fields. On the basis of the most accurate results for the electric field gradients, the nuclear quadrupole coupling constants for deuterium and oxygen-17 have been discussed including the temperature dependence. The final results are discussed in view of existing experimental data. Our best values for the nuclear quadrupole coupling constants are in excellent agreement ͑within 1%͒ of recent experimental results, while some earlier experimental values are shown to be less reliable.
The Journal of Physical Chemistry A, 2007
Adequate representation of the interactions that take place between water molecules has long been a goal of force field design. A full understanding of how the molecular charge distribution of water is altered by adjacent water molecules and by the hydrogen-bonding environment is a vital step toward achieving this task. For this purpose we generated ab initio electron densities of pure water clusters and hydrated serine and tyrosine. Quantum chemical topology enabled the study of a well-defined water molecule inside these clusters, by means of its volume, energy, and multipole moments. Intra-and intermolecular charge transfer was monitored and related to the polarization of water in hydrogen-bonded networks. Our analysis affords a way to define different types of water molecules in clusters.
A Molecular Density Functional Theory of Water
2016
Three dimensional implementations of liquid state theories offer an efficient alternative to computer simulations for the atomic-level description of aqueous solutions in complex environments. In this context, we present a (classical) molecular density functional theory (MDFT) of water that is derived from first principles and is based on two classical density fields, a scalar one, the particle density, and a vectorial one, the multipolar polarization density. Its implementa-tion requires as input the partial charge distribution of a water molecule and three measurable bulk properties, namely the structure factor and the k-dependent lon-gitudinal and transverse dielectric constants. It has to be complemented by a solute-solvent three-body term that reinforces tetrahedral order at short range. The approach is shown to provide the correct three-dimensional microscopic sol-vation profile around various molecular solutes, possibly possessing H-bonding sites, at a computer cost two-three...
Molecular Density Functional Theory of Water
The Journal of Physical Chemistry Letters, 2013
Three-dimensional implementations of liquid-state theories offer an efficient alternative to computer simulations for the atomic-level description of aqueous solutions in complex environments. In this context, we present a (classical) molecular density functional theory (MDFT) of water that is derived from first principles and is based on two classical density fields, a scalar one, the particle density, and a vectorial one, the multipolar polarization density. Its implementation requires as input the partial charge distribution of a water molecule and three measurable bulk properties, namely, the structure factor and the k-dependent longitudinal and transverse dielectric constants. It has to be complemented by a solute−solvent threebody term that reinforces tetrahedral order at short-range. The approach is shown to provide the correct 3-D microscopic solvation profile around various molecular solutes, possibly possessing H-bonding sites, at a computer cost two to three orders of magnitude lower than with explicit simulations.
Intermolecular Dynamics of Water: Suitability of Reactive Interatomic Potential
The Journal of Physical Chemistry B
Review of interatomic potentials TIP4P/2005f TIP4P exible model 1 (TIP4P/2005f) is extension of TIP4P model 2 (TIP4P/2005) with relaxation incorporated for intramolecular degree of freedom. TIP4P/2005f model consists of 4-point interaction site, H-O-H along with M site lying on bisector of H-O-H angle. Here two positive and one negative charges are placed at two hydrogens and at M-site respectively. Only oxygen is considered as Lennard-Jones interaction site. Hence LJ interaction between two water molecule is dened as
Systematic Improvement on the Classical Molecular Model of Water
Biophysical Journal, 2014
We report the iAMOEBA ("inexpensive AMOEBA") classical polarizable water model. The iAMOEBA model uses a direct approximation to describe electronic polarizability, in which the induced dipoles are determined directly from the permanent multipole electric fields and do not interact with one another. The direct approximation reduces the computational cost relative to a fully polarizable model such as AMOEBA. The model is parameterized using ForceBalance, a systematic optimization method that simultaneously utilizes training data from experimental measurements and high-level ab initio calculations. We show that iAMOEBA is a highly accurate model for water in the solid, liquid, and gas phases, with the ability to fully capture the effects of electronic polarization and predict a comprehensive set of water properties beyond the training data set including the phase diagram. The increased accuracy of iAMOEBA over the fully polarizable AMOEBA model demonstrates ForceBalance as a method that allows the researcher to systematically improve empirical models by efficiently utilizing the available data.
The Journal of Chemical Physics, 2010
We build on previous work ͓S. Y. Liem and P. L. A. Popelier, J. Chem. Theory Comput. 4, 353 ͑2008͔͒, where for the first time, a high-rank multipolar electrostatic potential was used in molecular dynamics simulations of liquid water at a wide range of pressures and temperatures, and using a multipolar Ewald summation. Water is represented as a rigid body, with atomic multipole moments defined by quantum chemical topology partitioning its gas phase electron density. The effect of the level of theory on the local structure of liquid water is systematically addressed. Values for Lennard-Jones ͑LJ͒ parameters are optimized, for both oxygen and hydrogen atoms, against bulk properties. The best LJ parameters were then used in a set of simulations at 30 different temperatures ͑1 atm͒ and another set at 11 different pressures ͑at 298 K͒. Inclusion of the hydrogen LJ parameters significantly increases the self-diffusion coefficient. The behavior of bulk properties was studied and the local water structure analyzed by both radial and spatial distribution functions. Comparisons with familiar point-charge potentials, such as TIP3P, TIP4P, TIP5P, and simple point charge, show the benefits of multipole moments.