Se presenta un estudio comparativo de cinco modelos moleculares diferentes para el agua evaluando el equilibrio de fases en la saturación mediante el método de simulación Monte Carlo en el ensamble de Gibbs. Los modelos moleculares comparados contienen tres o cuatro sitios con cargas. Se calculan las densidades de coexistencia líquido y vapor así como la presión de saturación a diferentes temperaturas. El punto crítico y las entalpías de vaporización se obtienen a partir de los resultados de la simulación. La comparación de los resultados de la simulación con datos experimentales obtenidos de la literatura muestra una buena concordancia.