We discuss cosmological hydrodynamic simulations of galaxy formation performed with the new moving-mesh code AREPO, which promises higher accuracy compared with the traditional SPH technique that has been widely employed for this problem. We use an i
dentical set of physics in corresponding simulations carried out with the well-tested SPH code GADGET, adopting also the same high-resolution gravity solver. We are thus able to compare both simulation sets on an object-by-object basis, allowing us to cleanly isolate the impact of different hydrodynamical methods on galaxy and halo properties. In accompanying papers, we focus on an analysis of the global baryonic statistics predicted by the simulation codes, (Vogelsberger et al. 2011) and complementary idealized simulations that highlight the differences between the hydrodynamical schemes (Sijacki et al. 2011). Here we investigate their influence on the baryonic properties of simulated galaxies and their surrounding haloes. We find that AREPO leads to significantly higher star formation rates for galaxies in massive haloes and to more extended gaseous disks in galaxies, which also feature a thinner and smoother morphology than their GADGET counterparts. Consequently, galaxies formed in AREPO have larger sizes and higher specific angular momentum than their SPH correspondents. The more efficient cooling flows in AREPO yield higher densities and lower entropies in halo centers (and the opposite trend in halo outskirts) leading to higher star formation rates of massive galaxies. While both codes agree to acceptable accuracy on a number of baryonic properties of cosmic structures, our results clearly demonstrate that galaxy formation simulations greatly benefit from the use of more accurate hydrodynamical techniques such as AREPO.
We study the formation of galaxies in a (50 Mpc/h)^3 cosmological simulation (2x288^3 particles), evolved using the entropy conserving SPH code Gadget-2. Most of the baryonic mass in galaxies of all masses is originally acquired through filamentary c
old mode accretion of gas that was never shock heated to its halo virial temperature, confirming the key feature of our earlier results obtained with a different SPH code (Keres et al. 2005). Atmospheres of hot, virialized gas develop in halos above ~2.5e11 Msun, a transition mass that is nearly constant from z=3 to z=0. Cold accretion persists in halos above the transition mass, especially at z>=2. It dominates the growth of galaxies in low mass halos at all times, and it is the main driver of the cosmic star formation history. Satellite galaxies have accretion rates similar to central galaxies of the same baryonic mass at high redshifts, but they have less accretion than comparable central galaxies at low redshift. Relative to our earlier results, the Gadget-2 simulations predict much lower rates of hot mode accretion from the virialized gas component of massive halos. At z<=1, typical hot accretion rates in halos above 5e12 Msun are below 1 Msun/yr, even though our simulation does not include AGN heating or other forms of preventive feedback. The inner density profiles of hot gas in these halos are shallow, with long associated cooling times. The cooling recipes typically used in semi-analytic models can overestimate the accretion rates in these halos by orders of magnitude, so such models may overemphasize the role of preventive feedback in producing observed galaxy masses and colors. A fraction of the massive halos develop cuspy profiles and significant cooling rates between z=1 and z=0, a redshift trend similar to the observed trend in the frequency of cooling flow clusters.
We compare the properties of galaxies that form in a cosmological simulation without strong feedback to observations at z=0. We confirm previous findings that models without strong feedback overproduce the observed galaxy baryonic mass function, espe
cially at the low and high mass extremes. Through post-processing we investigate what kinds of feedback would be required to reproduce observed galaxy masses and star formation rates. To mimic an extreme form of preventive feedback (e.g., AGN radio mode) we remove all baryonic mass that was originally accreted via hot mode from shock-heated gas. This does not bring the high mass end of the galaxy mass function into agreement with observations because much of the stellar mass in these systems formed at high redshift from baryons that originally accreted via cold mode onto lower mass progenitors. An efficient ejective feedback mechanism, such as supernova driven winds, must reduce the masses of these progenitors. Feedback must also reduce the masses of lower mass z=0 galaxies, which assemble at lower redshifts and have much lower star formation rates. If we monotonically re-map galaxy masses to reproduce the observed mass function, but retain the simulations predicted star formation rates, we obtain fairly good agreement with the observed sequence of star-forming galaxies but fail to recover the observed population of passive, low star formation rate galaxies. Suppressing all hot mode accretion improves agreement for high mass galaxies but worsens the agreement at intermediate masses. Reproducing these z=0 observations requires a feedback mechanism that dramatically suppresses star formation in a fraction of galaxies, increasing with mass, while leaving star formation rates of other galaxies essentially unchanged.
