Investigation of water-mediated intermolecular interactions with the adaptive resolution simulation technique

Nuclear quantum effects in fullerene-fullerene aggregation in water

We studied the effects of the quantum delocalization in space of the hydrogen atoms of water in the aggregation process of two fullerene molecules. We considered a case using a purely repulsive water-fullerene interaction, as such a situation has shown that water-mediated effects play a key role in the aggregation process. This study becomes feasible, at a reduced computational price, by combining the path integral (PI) molecular dynamics (MD) method with a recently developed open-system MD technique. Specifically, only the mandatory solvation shell of the two fullerene molecules was considered at full quantum resolution, while the rest of the system was represented as a mean-field macroscopic reservoir of particles and energy. Our results showed that the quantum nature of the hydrogen atoms leads to a sizable difference in the curve of the free energy of aggregation; that is, that nuclear quantum effects play a relevant role.

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.