We use Monte Carlo techniques and analytical methods to study the phase diagram of multicomponent Widom-Rowlinson models on a square lattice:there are M species all with the same fugacity z and a nearest neighbor hard core exclusion between unlike particles. Simulations show that for M between two and six there is a direct transition from the gas phase at z < z d (M ) to a demixed phase consisting mostly of one species at z > z d (M ) while for M ≥ 7 there is an intermediate "crystal phase" for z lying between z c (M ) and z d (M ). In this phase, which is driven by entropy, particles, independent of species, preferentially occupy one of the sublattices, i.e. spatial symmetry