Multi-cell Monte Carlo method for phase prediction