We study the transitions from band insulator to metal to Mott insulator in the ionic Hubbard model on a two dimensional square lattice using determinant Quantum Monte Carlo. Evaluation of the temperature dependence of the conductivity demonstrates that the metallic region extends for a finite range of interaction values. The Mott phase at strong coupling is accompanied by antiferromagnetic (AF) order. Inclusion of these intersite correlations changes the phase diagram qualitatively compared to dynamical mean field theory.