We present a new approach to study galaxy evolution in a cosmological context. We combine cosmological merger trees and semi-analytic models of galaxy formation to provide the initial conditions for multi-merger hydrodynamic simulations. In this way
we exploit the advantages of merger simulations (high resolution and inclusion of the gas physics) and semi-analytic models (cosmological background and low computational cost), and integrate them to create a novel tool. This approach allows us to study the evolution of various galaxy properties, including the treatment of the hot gaseous halo from which gas cools and accretes onto the central disc, which has been neglected in many previous studies. This method shows several advantages over other methods. As only the particles in the regions of interest are included, the run time is much shorter than in traditional cosmological simulations, leading to greater computational efficiency. Using cosmological simulations, we show that multiple mergers are expected to be more common than sequences of isolated mergers, and therefore studies of galaxy mergers should take this into account. In this pilot study, we present our method and illustrate the results of simulating ten Milky Way-like galaxies since z=1. We find good agreement with observations for the total stellar masses, star formation rates, cold gas fractions and disc scale length parameters. We expect that this novel numerical approach will be very useful for pursuing a number of questions pertaining to the transformation of galaxy internal structure through cosmic time.
We present a new statistical method to determine the relationship between the stellar masses of galaxies and the masses of their host dark matter haloes over the entire cosmic history from z~4 to the present. This multi-epoch abundance matching (MEAM
) model self-consistently takes into account that satellite galaxies first become satellites at times earlier than they are observed. We employ a redshift-dependent parameterization of the stellar-to-halo mass relation to populate haloes and subhaloes in the Millennium simulations with galaxies, requiring that the observed stellar mass functions at several redshifts be reproduced simultaneously. Using merger trees extracted from the dark matter simulations in combination with MEAM, we predict the average assembly histories of galaxies, separating into star formation within the galaxies (in-situ) and accretion of stars (ex-situ). The peak star formation efficiency decreases with redshift from 23% at z=0 to 9% at z=4 while the corresponding halo mass increases from 10^11.8Modot to 10^12.5Modot. The star formation rate of central galaxies peaks at a redshift which depends on halo mass; for massive haloes this peak is at early cosmic times while for low-mass galaxies the peak has not been reached yet. In haloes similar to that of the Milky-Way about half of the central stellar mass is assembled after z=0.7. In low-mass haloes, the accretion of satellites contributes little to the assembly of their central galaxies, while in massive haloes more than half of the central stellar mass is formed ex-situ with significant accretion of satellites at z<2. We find that our method implies a cosmic star formation history and an evolution of specific star formation rates which are consistent with those inferred directly. We present convenient fitting functions for stellar masses, star formation rates, and accretion rates as functions of halo mass and redshift.
Cosmological hydrodynamical simulations as well as observations indicate that spiral galaxies are comprised of five different components: dark matter halo, stellar disc, stellar bulge, gaseous disc and gaseous halo. While the first four components ha
ve been extensively considered in numerical simulations of binary galaxy mergers, the effect of a hot gaseous halo has usually been neglected even though it can contain up to 80% of the total gas within the galaxy virial radius. We present a series of hydrodynamic simulations of major mergers of disc galaxies, that for the first time include a diffuse, rotating, hot gaseous halo. Through cooling and accretion, the hot halo can dissipate and refuel the cold gas disc before and after a merger. This cold gas can subsequently form stars, thus impacting the morphology and kinematics of the remnant. Simulations of isolated systems with total mass M~10^12Msun show a nearly constant star formation rate of ~5Msun/yr if the hot gaseous halo is included, while the star formation rate declines exponentially if it is neglected. We conduct a detailed study of the star formation efficiency during mergers and find that the presence of a hot gaseous halo reduces the starburst efficiency (e=0.5) compared to simulations without a hot halo (e=0.68). Moreover we find cases where the stellar mass of the merger remnant is lower than the sum of the stellar mass of the two progenitor galaxies when evolved in isolation. This suggests a revision to semi-analytic galaxy formation models which assume that a merger always leads to enhanced star formation. We show that adding the hot gas component has a significant effect on the kinematics and internal structure of the merger remnants, like an increased abundance of fast rotators and an r^(1/4) surface brightness profile at small scales.
Deep pencil beam surveys (<1 deg^2) are of fundamental importance for studying the high-redshift universe. However, inferences about galaxy population properties are in practice limited by cosmic variance. This is the uncertainty in observational est
imates of the number density of galaxies arising from the underlying large-scale density fluctuations. This source of uncertainty can be significant, especially for surveys which cover only small areas and for massive high-redshift galaxies. Cosmic variance for a given galaxy population can be determined using predictions from cold dark matter theory and the galaxy bias. In this paper we provide tools for experiment design and interpretation. For a given survey geometry we present the cosmic variance of dark matter as a function of mean redshift z and redshift bin size Dz. Using a halo occupation model to predict galaxy clustering, we derive the galaxy bias as a function of mean redshift for galaxy samples of a given stellar mass range. In the linear regime, the cosmic variance of these galaxy samples is the product of the galaxy bias and the dark matter cosmic variance. We present a simple recipe using a fitting function to compute cosmic variance as a function of the angular dimensions of the field, z, Dz and stellar mass m*. We also provide tabulated values and a software tool. We find that for GOODS at z=2 and with Dz=0.5 the relative cosmic variance of galaxies with m*>10^11 Msun is ~38%, while it is ~27% for GEMS and ~12% for COSMOS. For galaxies of m*~10^10 Msun the relative cosmic variance is ~19% for GOODS, ~13% for GEMS and ~6% for COSMOS. This implies that cosmic variance is a significant source of uncertainty at z=2 for small fields and massive galaxies, while for larger fields and intermediate mass galaxies cosmic variance is less serious.
We study the effect of dissipational gas physics on the vertical heating and thickening of disc galaxies during minor mergers. We produce a suite of minor merger simulations for Milky Way-like galaxies. This suite consists of collisionless simulation
s as well as hydrodynamical runs including a gaseous component in the galactic disc. We find that in dissipationless simulations minor mergers cause the scale height of the disc to increase by up to a factor of ~2. When the presence of gas in the disc is taken into account this thickening is reduced by 25% (50%) for an initial disc gas fraction of 20% (40%), leading to a final scale height z0 between 0.6 and 0.7 kpc, regardless of the initial scale height. We argue that the presence of gas reduces disc heating via two mechanisms: absorption of kinetic impact energy by the gas and/or formation of a new thin stellar disc that can cause heated stars to recontract towards the disc plane. We show that in our simulations most of the gas is consumed during the merger and thus the regrowth of a new thin disc has a negligible impact on the z0 of the post merger galaxy. Final disc scale heights found in our simulations are in good agreement with studies of the vertical structure of spiral galaxies where the majority of the systems are found to have scale heights of 0.4 kpc < z0 < 0.8 kpc. We also found no tension between recent measurements of the scale height of the Milky Way thin disc and results coming from our hydrodynamical simulations. We conclude that the existence of a thin disc in the Milky Way and in external galaxies is not in obvious conflict with the predictions of the CDM model.
