  We study the properties of galaxies at z=2 in a Lambda CDM universe, using
two different types of hydrodynamic simulation methods (Eulerian TVD and SPH)
and a spectrophotometric analysis in the Un, G, R filter set. The simulated
galaxies at z=2 satisfy the color-selection criteria proposed by Adelberger et
al. (2004) when we assume Calzetti extinction with E(B-V)=0.15. We find that
the number density of simulated galaxies brighter than R<25.5 at z=2 is about
2e-2 h^3/Mpc^3, roughly one order of magnitude larger than that of Lyman break
galaxies at z=3. The most massive galaxies at z=2 have stellar masses >~1e11
Msun, and their observed-frame G-R colors lie in the range 0.0<G-R<1.0. They
typically have been continuously forming stars with a rate exceeding 30 Msun/yr
over a few Gyrs from z=10 to z=2, although the TVD simulation indicates a more
sporadic star formation history than the SPH simulations. Of order half of
their stellar mass was already assembled by z~4. The reddest massive galaxies
at z=2 with G-R >= 1.0 and Mstar>1e10 Msun/h finished the build-up of their
stellar mass by z~3. Interestingly, our study suggests that the majority of the
most massive galaxies at z=2 should be detectable at rest-frame UV wavelengths,
contrary to some recent claims made on the basis of near-IR studies of galaxies
at the same epoch, provided the median extinction is less than E(B-V)<0.3.
However, our results also suggest that the fraction of stellar mass contained
in galaxies that pass the color-selection criteria could be as low as 50% of
the total stellar mass in the Universe at z=2. Our simulations suggest that the
missing stellar mass is contained in fainter (R>25.5) and intrinsically redder
galaxies. Our results do not suggest that hierarchical galaxy formation fails
to account for the massive galaxies at z>=1. (abridged)
