Alyson M. Brooks

The Milky Way Tomography with SDSS. I. Stellar Number Density Distribution

Using the photometric parallax method we estimate the distances to ~48 million stars detected by the Sloan Digital Sky Survey (SDSS) and map their three-dimensional number density distribution in the Galaxy. The currently available data sample the distance range from 100 pc to 20 kpc and cover 6500 deg2 of sky, mostly at high Galactic latitudes (|b| > 25). These stellar number density maps allow an investigation of the Galactic structure with no a priori assumptions about the functional form of its components. The data show strong evidence for a Galaxy consisting of an oblate halo, a disk component, and a number of localized overdensities. The number density distribution of stars as traced by M dwarfs in the solar neighborhood (D < 2 kpc) is well fit by two exponential disks (the thin and thick disk) with scale heights and lengths, bias corrected for an assumed 35% binary fraction, of H1 = 300 pc and L1 = 2600 pc, and H2 = 900 pc and L2 = 3600 pc, and local thick-to-thin disk density normalization ρthick(R☉)/ρthin(R☉) = 12% . We use the stars near main-sequence turnoff to measure the shape of the Galactic halo. We find a strong preference for oblate halo models, with best-fit axis ratio c/a = 0.64, ρH ∝ r−2.8 power-law profile, and the local halo-to-thin disk normalization of 0.5%. Based on a series of Monte Carlo simulations, we estimate the errors of derived model parameters not to be larger than ~20% for the disk scales and ~10% for the density normalization, with largest contributions to error coming from the uncertainty in calibration of the photometric parallax relation and poorly constrained binary fraction. While generally consistent with the above model, the measured density distribution shows a number of statistically significant localized deviations. In addition to known features, such as the Monoceros stream, we detect two overdensities in the thick disk region at cylindrical galactocentric radii and heights (R,Z) ~ (6.5,1.5) kpc and (R,Z) ~ (9.5,0.8) kpc and a remarkable density enhancement in the halo covering over 1000 deg2 of sky toward the constellation of Virgo, at distances of ~6-20 kpc. Compared to counts in a region symmetric with respect to the l = 0° line and with the same Galactic latitude, the Virgo overdensity is responsible for a factor of 2 number density excess and may be a nearby tidal stream or a low-surface brightness dwarf galaxy merging with the Milky Way. The u − g color distribution of stars associated with it implies metallicity lower than that of thick disk stars and consistent with the halo metallicity distribution. After removal of the resolved overdensities, the remaining data are consistent with a smooth density distribution; we detect no evidence of further unresolved clumpy substructure at scales ranging from ~50 pc in the disk to ~1-2 kpc in the halo.

Forming disc galaxies in ΛCDM simulations

We used fully cosmological, high resolution N-body + SPH simulations to follow the formation of disk galaxies with rotational velocities between 135 and 270 km/sec in a ΛCDM universe. The simulations include gas cooling, star formation, the effects of a uniform UV background and a physically motivated description of feedback from supernovae. The host dark matter halos have a spin and last major merger redshift typical of galaxy sized halos as measured in recent large scale N–Body simulations. The simulated galaxies form rotationally supported disks with realistic exponential scale lengths and fall on both the I-band and baryonic Tully Fisher relations. An extended stellar disk forms inside the Milky Way sized halo immediately after the last major merger. The combination of UV background and SN feedback drastically reduces the number of visible satellites orbiting inside a Milky Way sized halo, bringing it in fair agreement with observations. Our simulations predict that the average age of a primary galaxy’s stellar population decreases with mass, because feedback delays star formation in less massive galaxies. Galaxies have stellar masses and current star formation rates as a function of total mass that are in good agreement with observational data. We discuss how both high mass and force resolution and a realistic description of star formation and feedback are important ingredients to match the observed properties of galaxies.

THE ROLE OF COLD FLOWS IN THE ASSEMBLY OF GALAXY DISKS

We use high-resolution cosmological hydrodynamical simulations to demonstrate that cold flow gas accretion, particularly along filaments, modifies the standard picture of gas accretion and cooling onto galaxy disks. In the standard picture, all gas is initially heated to the virial temperature of the galaxy as it enters the virial radius. Low-mass galaxies are instead dominated by accretion of gas that stays well below the virial temperature, and even when a hot halo is able to develop in more massive galaxies there exist dense filaments that penetrate inside of the virial radius and deliver cold gas to the central galaxy. For galaxies up to ~L*, this cold accretion gas is responsible for the star formation (SF) in the disk at all times to the present. Even for galaxies at higher masses, cold flows dominate the growth of the disk at early times. Within this modified picture, galaxies are able to accrete a large mass of cold gas, with lower initial gas temperatures leading to shorter cooling times to reach the disk. Although SF in the disk is mitigated by supernovae feedback, the short cooling times allow for the growth of stellar disks at higher redshifts than predicted by the standard model.

BARYONS MATTER: WHY LUMINOUS SATELLITE GALAXIES HAVE REDUCED CENTRAL MASSES

Using high-resolution cosmological hydrodynamical simulations of Milky Way-massed disk galaxies, we demonstrate that supernovae feedback and tidal stripping lower the central masses of bright (–15 < MV < –8) satellite galaxies. These simulations resolve high-density regions, comparable to giant molecular clouds, where stars form. This resolution allows us to adopt a prescription for H2 formation and destruction that ties star formation to the presence of shielded, molecular gas. Before infall, supernova feedback from the clumpy, bursty star formation captured by this physically motivated model leads to reduced dark matter (DM) densities and shallower inner density profiles in the massive satellite progenitors (M vir ≥ 109 M ☉, M * ≥ 107 M ☉) compared with DM-only simulations. The progenitors of the lower mass satellites are unable to maintain bursty star formation histories, due to both heating at reionization and gas loss from initial star-forming events, preserving the steep inner density profile predicted by DM-only simulations. After infall, gas stripping from satellites reduces the total central masses of satellites simulated with DM+baryons relative to DM-only satellites. Additionally, enhanced tidal stripping after infall due to the baryonic disk acts to further reduce the central DM densities of the luminous satellites. Satellites that enter with cored DM halos are particularly vulnerable to the tidal effects of the disk, exacerbating the discrepancy in the central masses predicted by baryon+DM and DM-only simulations. We show that DM-only simulations, which neglect the highly non-adiabatic evolution of baryons described in this work, produce denser satellites with larger central velocities. We provide a simple correction to the central DM mass predicted for satellites by DM-only simulations. We conclude that DM-only simulations should be used with great caution when interpreting kinematic observations of the Milky Ways dwarf satellites.

Damped Lyman α systems in galaxy formation simulations

We investigate the population of z = 3 damped Lyman alpha systems (DLAs) in a recent series of high resolution galaxy formation simulations. The simulations are of interest because they form at z = 0 some of the most realistic disk galaxies to date. No free parameters are available in our study: the simulation parameters have been fixed by physical and z = 0 observational constraints, and thus our work provides a genuine consistency test. The precise role of DLAs in galaxy formation remains in debate, but they provide a number of strong constraints on the nature of our simulated bound systems at z = 3 because of their coupled information on neutral H I densities, kinematics, metallicity and estimates of star f ormation activity. Our results, without any parameter-tuning, closely match the observed incidence rate and column density distributions of DLAs. Our simulations are the first to reproduce the distribution of metallicities (with a median of ZDLA ≃ Z⊙/20) without invoking observationally unsupported mechanisms such as significant dust biasin g. This is especially encouraging given that these simulations have previously been shown to have a realistic 0 < z < 2 stellar mass-metallicity relation. Additionally, we see a strong p ositive correlation between sightline metallicity and low-ion velocity width, the normalization and slope of which comes close to matching recent observational results. However, we somewhat underestimate the number of observed high velocity width systems; the severity of this disagreement is comparable to other recent DLA-focused studies. DLAs in our simulations are predominantly associated with dark matter haloes with virial masses in the range 10 9 < Mvir/M⊙ < 10 11 . We are able to probe DLAs at high resolution, irrespective of their masses, by using a range of simulation s of differing volumes. The fully constrained feedback prescription in use causes the majority of DLA haloes to form stars at a very low rate, accounting for the low metallicities. It i s also responsible for the massmetallicity relation which appears essential for reproduc ing the velocity-metallicity correlation. By z = 0 the majority of the z = 3 neutral gas forming the DLAs has been converted into stars, in agreement with rough physical expectations.

Hierarchical formation of bulgeless galaxies: why outflows have low angular momentum

Using high resolution, fully cosmological smoothed particle hydrodynamical simulations of dwarf galaxies in a Lambda cold dark matter Universe, we show how high redshift gas outflows can modify the baryon angular momentum distribution and allow pure disc galaxies to form. We outline how galactic outflows preferentially remove low angular momentum material due a combination of (a) star formation peaking at high redshift in shallow dark matter potentials, an epoch when accreted gas has relatively low angular momentum, (b) the existence of an extended reservoir of high angular momentum gas in the outer disc to provide material for prolonged SF at later times and (c) the tendency for outflows to follow the path of least resistance which is perpendicular to the disc. We also show that outflows are enhanced during mergers, thus expelling much of the gas which has lost its angular momentum during these events, and preventing the formation of ‘classical’, merger driven bulges in low-mass systems. Stars formed prior to such mergers form a diffuse, extended stellar halo component similar to those detected in nearby dwarfs.

The Origin and Evolution of the Mass-Metallicity Relationship for Galaxies: Results from Cosmological N-Body Simulations

We examine the origin and evolution of the mass-metallicity relationship (MZR, M*-Z) for galaxies using high-resolution cosmological smoothed particle hydrodynamics (SPH) + N-body simulations that include a physically motivated description of supernova feedback and subsequent metal enrichment. We discriminate between two sources that may contribute to the origin of the MZR: (1) metal and baryon loss due to gas outflow or (2) inefficient star formation at the lowest galaxy masses. Our simulated galaxies reproduce the observed MZR in shape and normalization at both z = 0 and 2. We find that baryon loss occurs due to UV heating before star formation turns on in galaxies with Mbar < 108 M☉, but that some gas loss due to supernova-induced winds is required to subsequently reproduce the low effective chemical yield observed in low-mass galaxies. Despite this, we show that low star formation efficiencies, regulated by supernova feedback, are primarily responsible for the lower metallicities of low-mass galaxies and the overall M*-Z trend. We find that the shape of the MZR is relatively constant with redshift but that its normalization increases with time. Simulations with no energy feedback from supernovae overproduce metals at low galaxy masses by rapidly transforming a large fraction of their gas into stars. Despite the fact that our low-mass galaxies have lost a majority of their baryons, they are still the most gas-rich objects in our simulations due to their low star formation efficiencies.

The SAMI Galaxy Survey: instrument specification and target selection

The SAMI Galaxy Survey will observe 3400 galaxies with the Sydney-AAO Multi- object Integral-field spectrograph (SAMI) on the Anglo-Australian Telescope (AAT) in a 3-year survey which began in 2013. We present the throughput of the SAMI system, the science basis and specifications for the target selection, the survey observation plan and the combined properties of the selected galaxies. The survey includes four volume-limited galaxy samples based on cuts in a proxy for stellar mass, along with low-stellar-mass dwarf galaxies all selected from the Galaxy And Mass Assembly (GAMA) survey. The GAMA regions were selected because of the vast array of ancillary data available, including ultraviolet through to radio bands. These fields are on the celestial equator at 9, 12, and 14.5 hours, and cover a total of 144 square degrees (in GAMA-I). Higher density environments are also included with the addition of eight clusters. The clusters have spectroscopy from 2dFGRS and SDSS and photometry in regions covered by the Sloan Digital Sky Survey (SDSS) and/or VLT Survey Telescope/ATLAS. The aim is to cover a broad range in stellar mass and environment, and therefore the primary survey targets cover redshifts 0.004 < z < 0.095, magnitudes rpet < 19.4, stellar masses 107– 1012M⊙, and environments from isolated field galaxies through groups to clusters of _ 1015M⊙.

Why Baryons Matter: The Kinematics of Dwarf Spheroidal Satellites

We use some of the highest resolution cosmological simulations ever produced of Milky Way-mass galaxies that include both baryons and dark matter to show that baryonic physics (energetic feedback from supernovae and subsequent tidal stripping) significantly reduces the dark matter mass in the central regions of luminous satellite galaxies. The reduced central masses of the simulated satellites reproduce the observed internal dynamics of Milky Way and M31 satellites as a function of luminosity. Including baryonic physics in Cold Dark Matter models naturally explains the observed low dark matter densities in the Milky Way’s dwarf spheroidal population. Our simulations therefore resolve the tension between kinematics predicted in Cold Dark Mater theory and observations of satellites, without invoking alternative forms of dark matter.

A Baryonic Solution to the Missing Satellites Problem

It has been demonstrated that the inclusion of baryonic physics can alter the dark matter densities in the centers of low-mass galaxies, making the central dark matter slope more shallow than predicted in pure cold dark matter simulations. This flattening of the dark matter profile can occur in the most luminous subhalos around Milky Way mass galaxies. Zolotov et al. have suggested a correction to be applied to the central masses of dark matter-only satellites in order to mimic the affect of (1) the flattening of the dark matter cusp due to supernova feedback in luminous satellites and (2) enhanced tidal stripping due to the presence of a baryonic disk. In this paper, we apply this correction to the z = 0 subhalo masses from the high resolution, dark matter-only Via Lactea II (VL2) simulation, and find that the number of massive subhalos is dramatically reduced. After adopting a stellar mass to halo mass relationship for the VL2 halos, and identifying subhalos that are (1) likely to be destroyed by stripping and (2) likely to have star formation suppressed by photo-heating, we find that the number of massive, luminous satellites around a Milky Way mass galaxy is in agreement with the number of observed satellites around the Milky Way or M31. We conclude that baryonic processes have the potential to solve the missing satellites problem.