Path Integral Molecular Dynamics of Liquid Water in a Mean-Field Particle Reservoir

We present a simulation scheme for path integral simulation of molecular liquids where a small open region is embedded in a large reservoir of non interacting point-particles. The scheme is based on the latest development of the adaptive resolution technique AdResS and allows for the space-dependent change of molecular resolution from a path integral representation with 120 degrees of freedom to a point particle that does not interact with other molecules and vice versa. The method is applied to liquid water and implies a sizable gain regarding the request of computational resources compared to full path integral simulations. Given the role of water as universal solvent with a specific hydrogen bonding network, the path integral treatment of water molecules is important to describe the quantum effects of hydrogen atoms’ delocalization in space on the hydrogen bonding network. The method presented here implies feasible computational efforts compared to full path integral simulations of liquid water which, on large scales, are often prohibitive.