Rapid estimation of electronic degrees of freedom in Monte Carlo calculations for polarizable models of liquid water

An improved simple polarisable water model for use in biomolecular simulation

The accuracy of biomolecular simulations depends to some degree on the accuracy of the water model used to solvate the biomolecules. Because many biomolecules such as proteins are electrostatically rather inhomogeneous, containing apolar, polar, and charged moieties or side chains, a water model should be able to represent the polarisation response to a local electrostatic field, while being compatible with the force field used to model the biomolecules or protein. The two polarisable water models, COS/G2 and COS/D, that are compatible with the GROMOS biomolecular force fields leave room for improvement. The COS/G2 model has a slightly too large dielectric permittivity and the COS/D model displays a much too slow dynamics. The proposed COS/D2 model has four interaction sites: only one Lennard-Jones interaction site, the oxygen atom, and three permanent charge sites, the two hydrogens, and one massless off-atom site that also serves as charge-on-spring (COS) polarisable site with a damped or sub-linear dependence of the induced dipole on the electric field strength for large values of the latter. These properties make it a cheap and yet realistic water model for biomolecular solvation.

Proton momentum distribution in water: an open path integral molecular dynamics study

Recent neutron Compton scattering experiments have detected the proton momentum distribution in water. The theoretical calculation of this property can be carried out via "open" path integral expressions. In this work, present an extension of the staging path integral molecular dynamics method, which is then employed to calculate the proton momentum distributions of water in the solid, liquid, and supercritical phases. We utilize a flexible, single point charge empirical force field to model the system's interactions. The calculated momentum distributions depict both agreement and discrepancies with experiment. The differences may be explained by the deviation of the force field from the true interactions. These distributions provide an abundance of information about the environment and interactions surrounding the proton.

Polarization and charge-transfer effects in aqueous solution via ab initio QM/MM simulations

Combined ab initio quantum mechanical and molecular mechanical (QM/MM) simulations coupled with the block-localized wave function energy decomposition (BLW-ED) method have been conducted to study the solvation of two prototypical ionic systems, acetate and methylammonium ions in aqueous solution. Calculations reveal that the electronic polarization between the targeted solutes and water is the primary many-body effect, whereas the charge-transfer term only makes a small fraction of the total solute-solvent interaction energy. In particular, the polarization effect is dominated by the solvent (water) polarization.

A computationally inexpensive modification of the point dipole electrostatic polarization model for molecular simulations

We present an approximation, which allows reduction of computational resources needed to explicitly incorporate electrostatic polarization into molecular simulations utilizing empirical force fields. The proposed method is employed to compute three-body energies of molecular complexes with dipolar electrostatic probes, gas-phase dimerization energies, and pure liquid properties for five systems that are important in biophysical and organic simulations-water, methanol, methylamine, methanethiol, and acetamide. In all the cases, the three-body energies agreed with high level ab initio data within 0.07 kcal/mol, dimerization energies-within 0.43 kcal/mol (except for the special case of the CH(3)SH), and computed heats of vaporization and densities differed from the experimental results by less than 2%. Moreover, because the presented method allows a significant reduction in computational cost, we were able to carry out the liquid-state calculations with Monte Carlo technique. Comparison with the full-scale point dipole method showed that the computational time was reduced by 3.5 to more than 20 times, depending on the system in hand and on the desired level of the full-scale model accuracy, while the difference in energetic results between the full-scale and the presented approximate model was not great in the most cases. Comparison with the nonpolarizable OPLS-AA force field for all the substances involved and with the polarizable POL3 and q90 models for water and methanol, respectively, demonstrates that the presented technique allows reduction of computational cost with no sacrifice of accuracy. We hope that the proposed method will be of benefit to research employing molecular modeling technique in the biophysical and physical organic chemistry areas.