The multi-particle-move Monte Carlo methodology is implemented for Gibbs ensemble Monte Carlo (GEMC) simulations and applied to two recent Baranyai-Kiss polarizable models of water. In addition to obtaining and reporting the GEMC data for the considered models, attention is paid also to technical parts of simulations, particularly to the long-range corrections. A detailed analysis identifies the reasons for discrepancies found between the present results and those reported in literature. The problem lies in the application of incomplete long-range corrections in previous studies. It turns out that in fact both considered models are in better agreement with experiment than originally reported. (C) 2013 Elsevier B.V. All rights reserved.