We show that the cool gas masses of galactic discs reach a steady state that lasts many Gyr after their last major merger in cosmological hydrodynamic simulations. The mass of disc gas, M$_{rm gas}$, depends upon a galaxy halos spin and virial mass,
but not upon stellar feedback. Halos with low spin have high star formation efficiency and lower disc gas mass. Similarly, lower stellar feedback leads to more star formation so the gas mass ends up nearly the same irregardless of stellar feedback strength. Even considering spin, the M$_{rm gas}$ relation with halo mass, M$_{200}$ only shows a factor of 3 scatter. The M$_{rm gas}$--M$_{200}$ relation show a break at M$_{200}$=$2times10^{11}$ M$_odot$ that corresponds to an observed break in the M$_{rm gas}$--M$_star$ relation. The constant disc mass stems from a shared halo gas density profile in all the simulated galaxies. In their outer regions, the profiles are isothermal. Where the profile rises above $n=10^{-3}$ cm$^{-3}$, the gas readily cools and the profile steepens. Inside the disc, rotation supports gas with a flatter density profile except where supernova explosions disrupt the disc. Energy injection from stellar feedback also provides pressure support to the halo gas to prevent runaway cooling flows. The resulting constant gas mass makes simpler models for galaxy formation possible, either using a bathtub model for star formation rates or when modeling chemical evolution.
We analyse the structure and chemical enrichment of a Milky Way-like galaxy with a stellar mass of 2 10^{10} M_sun, formed in a cosmological hydrodynamical simulation. It is disk-dominated with a flat rotation curve, and has a disk scale length simil
ar to the Milky Ways, but a velocity dispersion that is ~50% higher. Examining stars in narrow [Fe/H] and [alpha/Fe] abundance ranges, we find remarkable qualitative agreement between this simulation and observations: a) The old stars lie in a thickened distribution with a short scale length, while the young stars form a thinner disk, with scale lengths decreasing, as [Fe/H] increases. b) Consequently, there is a distinct outward metallicity gradient. c) Mono-abundance populations exist with a continuous distribution of scale heights (from thin to thick). However, the simulated galaxy has a distinct and substantive very thick disk (h_z~1.5 kpc), not seen in the Milky Way. The broad agreement between simulations and observations allows us to test the validity of observational proxies used in the literature: we find in the simulation that mono-abundance populations are good proxies for single age populations (<1 Gyr) for most abundances.
We introduce the Making Galaxies in a Cosmological Context (MaGICC) program of smoothed particle hydrodynamics (SPH) simulations. We describe a parameter study of galaxy formation simulations of an L* galaxy that uses early stellar feedback combined
with supernova feedback to match the stellar mass--halo mass relationship. While supernova feedback alone can reduce star formation enough to match the stellar mass--halo mass relationship, the galaxy forms too many stars before z=2 to match the evolution seen using abundance matching. Our early stellar feedback is purely thermal and thus operates like a UV ionization source as well as providing some additional pressure from the radiation of massive, young stars. The early feedback heats gas to >10^6 K before cooling to 10^4 K. The pressure from this hot gas creates a more extended disk and prevents more star formation prior to z=1 than supernovae feedback alone. The resulting disk galaxy has a flat rotation curve, an exponential surface brightness profile, and matches a wide range of disk scaling relationships. The disk forms from the inside-out with an increasing exponential scale length as the galaxy evolves. Overall, early stellar feedback helps to simulate galaxies that match observational results at low and high redshifts.
We explore the circumgalactic medium (CGM) of two simulated star-forming galaxies with luminosities L ~ 0.1 and 1 L* generated using the smooth particle hydrodynamic code GASOLINE. These simulations are part of the Making Galaxies In a Cosmological C
ontext (MAGICC) program in which the stellar feedback is tuned to match the stellar mass-halo mass relationship. For comparison, each galaxy was also simulated using a lower feedback (LF) model which has strength comparable to other implementations in the literature. The MAGICC feedback (MF) model has a higher incidence of massive stars and an approximately two times higher energy input per supernova. Apart from the low-mass halo using LF, each galaxy exhibits a metal-enriched CGM that extends to approximately the virial radius. A significant fraction of this gas has been heated in supernova explosions in the disc and subsequently ejected into the CGM where it is predicted to give rise to substantial O VI absorption. The simulations do not yet address the question of what happens to the O VI when the galaxies stop forming stars. Our models also predict a reservoir of cool H I clouds that show strong Lyalpha absorption to several hundred kpc. Comparing these models to recent surveys with the Hubble Space Telescope, we find that only the MF models have sufficient O VI and H I gas in the CGM to reproduce the observed distributions. In separate analyses, these same MF models also show better agreement with other galaxy observables (e.g. rotation curves, surface brightness profiles and H I gas distribution). We infer that the CGM is the dominant reservoir of baryons for galaxy haloes.
We present the McMaster Unbiased Galaxy Simulations (MUGS), the first 9 galaxies of an unbiased selection ranging in total mass from 5$times10^{11}$ M$_odot$ to 2$times10^{12}$ M$_odot$ simulated using n-body smoothed particle hydrodynamics (SPH) at
high resolution. The simulations include a treatment of low temperature metal cooling, UV background radiation, star formation, and physically motivated stellar feedback. Mock images of the simulations show that the simulations lie within the observed range of relations such as that between color and magnitude and that between brightness and circular velocity (Tully-Fisher). The greatest discrepancy between the simulated galaxies and observed galaxies is the high concentration of material at the center of the galaxies as represented by the centrally peaked rotation curves and the high bulge-to-total ratios of the simulations determined both kinematically and photometrically. This central concentration represents the excess of low angular momentum material that long has plagued morphological studies of simulated galaxies and suggests that higher resolutions and a more accurate description of feedback will be required to simulate more realistic galaxies. Even with the excess central mass concentrations, the simulations suggest the important role merger history and halo spin play in the formation of disks.
Stellar population studies show that low mass galaxies in all environments exhibit stellar halos that are older and more spherically distributed than the main body of the galaxy. In some cases, there is a significant intermediate age component that e
xtends beyond the young disk. We examine a suite of Smoothed Particle Hydrodynamic (SPH) simulations and find that elevated early star formation activity combined with supernova feedback can produce an extended stellar distribution that resembles these halos for model galaxies ranging from $v_{200}$ = 15 km s$^{-1}$ to 35 km s$^{-1}$, without the need for accretion of subhalos.
We simulate the collapse of isolated dwarf galaxies using SPH + N-Body simulations including a physically motivated description of the effects of supernova feedback. As the gas collapses and stars form, the supernova feedback disrupts enough gas to t
emporarily quench star formation. The gas flows outward into a hot halo, where it cools until star formation can continue once more and the cycle repeats. The star formation histories of isolated Local Group dwarf galaxies exhibit similar episodic bursts of star formation. We examine the mass dependence of the stellar velocity dispersions and find that they are no less than half the velocity of the halos measured at the virial radius.
