On the re-engineered TIP4P water models for the prediction of vapor-liquid equilibrium

We perform extensive Gibbs Ensemble Monte Carlo simulations to study the capability of some recently re-parameterizations of the original TIP4P model intended to predict accurately the vapor-liquid coexistence envelope of water, its critical point, and its temperature dependence for the vapor pressure and second virial coefficient, and complement this analysis with the characterization of some specific crystalline faces of ice. We also disclose some trends between the resulting dipole moment of the models and the Lennard-Jones parameters, the location of the negative charge, as well as the estimated critical temperature. Finally, we discuss the inability of these models to predict accurately and simultaneously the melting temperature and the temperature of maximum density.