I present our latest attempts at simulating the observed galaxy population at high and low redshifts using high-resolution cosmological hydrodynamic simulations together with population synthesis models, and compare them with observations. At high redshift, our simulations agree reasonably well with the observed "Lyman break" galaxy population provided a significant amount of dust is included. They indicate that Lyman break galaxies are the most massive galaxies at z~3, and have been forming stars quiescently for over a Gyr. Galaxy clustering from z=3->0 shows little evolution, while the dark matter clustering grows significantly, indicating that galaxies at high redshift are much more highly biased tracers than those at low redshift.