A coupled density functional-molecular mechanics Monte Carlo simulation method: The water molecule in liquid water

A theoretical model to investigate chemical processes in solution is described. It is based on the use of a coupled density functional/molecular mechanics Hamiltonian. The most interesting feature of the method is that it allows a detailed study of the solute's electronic distribution and of its fluctuations. We present the results for isothermal-isobaric constant-NPT Monte Carlo simulation of a water molecule in liquid water. The quantum subsystem is described using a double-zeta quality basis set with polarization orbitals and nonlocal exchange-correlation corrections. The classical system is constituted by 128 classical TIP3P or Simple Point Charge (SPC) water molecules. The atom-atom radial distribution functions present a good agreement with the experimental curves. Differences with respect to the classical simulation are discussed. The instantaneous and the averaged polarization of the quantum molecule are also analyzed.