Nuclear quantum effects in water clusters: the role of the molecular flexibility

With the objective of establishing the importance of water flexibility in empirical models which explicitly include nuclear quantum effects, we have carried out path integral Monte Carlo simulations in water clusters with up to seven molecules. Two recently developed models have been used for comparison: the rigid TIP4PQ/2005 and the flexible q-TIP4P/F models, both inspired by the rigid TIP4P/2005 model. To obtain a starting configuration for our simulations, we have located the global minima for the rigid TIP4P/2005 and TIP4PQ/2005 models and for the flexible q-TIP4P/F model. All the structures are similar to those predicted by the rigid TIP4P potential showing that the charge distribution mainly determines the global minimum structure. For the flexible q-TIP4P/F model, we have studied the geometrical distortion upon isotopic substitution by studying tritiated water clusters. Our results show that tritiated water clusters exhibit an r(OT) distance shorter than the r(OH) distance in water clusters, not significant changes in the Phi(HOH) angle, and a lower average dipole moment than water clusters. We have also carried out classical simulations with the rigid TIP4PQ/2005 model showing that the rotational kinetic energy is greatly affected by quantum effects, but the translational kinetic energy is only slightly modified. The potential energy is also noticeably higher than in classical simulations. Finally, as a concluding remark, we have calculated the formation energies of water clusters using both models, finding that the formation energies predicted by the rigid TIP4PQ/2005 model are lower by roughly 0.6 kcal/mol than those of the flexible q-TIP4P/F model for clusters of moderate size, the origin of this difference coming mainly from the geometrical distortion of the water molecule in the clusters that causes an increase in the intramolecular potential energy.