Path Integral Monte Carlo Simulations of Flexible Water Clusters with Normal-Mode Sampling in Jacobi Coordinates