We analyze the cold dark matter density profiles of 54 galaxy halos simulated with FIRE-2 galaxy formation physics, each resolved within $0.5\%$ of the halo virial radius. These halos contain galaxies with masses that range from ultra-faint dwarfs ($M_\star \simeq 10^{4.5} M_{\odot}$) to the largest spirals ($M_\star \simeq 10^{11} M_{\odot}$) and have density profiles that are both cored and cuspy. We characterize our results using a new analytic density profile that extends the standard Einasto form to allow for a pronounced constant-density core in the resolved innermost radius. With one additional core-radius parameter, $r_{c}$, this "core-Einasto" profile is able to characterize the shape and normalization of our feedback-impacted dark matter halos. In order to enable comparisons with observations, we provide fitting functions for $r_{c}$ and other profile parameters as a function of both $M_\star$ and $M_{\star}/M_{\rm halo}$. In agreement with similar studies done in the literature, we find that dark matter core formation is most efficient at the characteristic stellar-mass to halo-mass ratio $M_\star/M_{\rm halo} \simeq 5 \times 10^{-3}$, or $M_{\star} \sim 10^9 \, M_{\odot}$, with cores that are roughly the size of the galaxy half-light radius, $r_{c} \simeq 1-5$ kpc. Furthermore, we find no evidence for core formation at radii $\gtrsim 100\ \rm pc$ in galaxies with $M_{\star}/M_{\rm halo} < 5\times 10^{-4}$ or $M_\star \lesssim 10^6 \, M_{\odot}$. For Milky Way-size galaxies, baryonic contraction often makes halos significantly more concentrated and dense at the stellar half-light radius than dark matter only runs. However, even at the Milky Way scale, FIRE-2 galaxy formation still produces small dark matter cores of $\simeq 0.5-2$ kpc in size. Recent evidence for a ${\sim} 2$ kpc core in the Milky Way's dark matter halo is consistent with this expectation.

While many tensions between Local Group (LG) satellite galaxies and LCDM cosmology have been alleviated through recent cosmological simulations, the spatial distribution of satellites remains an important test of physical models and physical versus numerical disruption in simulations. Using the FIRE-2 cosmological zoom-in baryonic simulations, we examine the radial distributions of satellites with Mstar > 10^5 Msun around 8 isolated Milky Way- (MW) mass host galaxies and 4 hosts in LG-like pairs. We demonstrate that these simulations resolve the survival and physical destruction of satellites with Mstar >~ 10^5 Msun. The simulations broadly agree with LG observations, spanning the radial profiles around the MW and M31. This agreement does not depend strongly on satellite mass, even at distances <~ 100 kpc. Host-to-host variation dominates the scatter in satellite counts within 300 kpc of the hosts, while time variation dominates scatter within 50 kpc. More massive host galaxies within our sample have fewer satellites at small distances, likely because of enhanced tidal destruction of satellites via the baryonic disks of host galaxies. Furthermore, we quantify and provide fits to the tidal depletion of subhalos in baryonic relative to dark matter-only simulations as a function of distance. Our simulated profiles imply observational incompleteness in the LG even at Mstar >~ 10^5 Msun: we predict 2-10 such satellites to be discovered around the MW and possibly 6-9 around M31. To provide cosmological context, we compare our results with the radial profiles of satellites around MW analogs in the SAGA survey, finding that our simulations are broadly consistent with most SAGA systems.

Starting with the spherical Jeans equation, we show that there exists a radius where the mass enclosed depends only on the projected tangential velocity dispersion, assuming that the anisotropy profile slowly varies. This is well-approximated at the radius where the log-slope of the stellar tracer profile is $-2$: $r_{-2}$. The associated mass is $M(r_{-2}) = 2 G^{-1} \langle \sigma_{\mathcal{T}}^{2}\rangle^{*} r_{-2}$ and the circular velocity is $V^{2}({r_{-2}}) = 2\langle \sigma_{\mathcal{T}}^{2}\rangle^{*}$. For a Plummer profile $r_{-2} \simeq 4R_e/5$. Importantly, $r_{-2}$ is smaller than the characteristic radius for line-of-sight velocities derived by Wolf et al. 2010. Together, the two estimators can constrain the mass profiles of dispersion-supported galaxies. We illustrate its applicability using published proper motion measurements of dwarf galaxies Draco and Sculptor, and find that they are consistent with inhabiting cuspy NFW subhalos of the kind predicted in CDM but we cannot rule out a core. We test our combined mass estimators against previously-published, non-spherical cosmological dwarf galaxy simulations done in both CDM and SIDM. For CDM, the estimates for the dynamic rotation curves are found to be accurate to $10\%$ while SIDM are accurate to $15\%$. Unfortunately, this level of accuracy is not good enough to measure slopes at the level required to distinguish between cusps and cores of the type predicted in viable SIDM models without stronger priors. However, we find that this provides good enough accuracy to distinguish between the normalization differences predicted at small radii ($r \simeq r_{-2} < r_{\rm core}$) for interesting SIDM models. As the number of galaxies with internal proper motions increases, mass estimators of this kind will enable valuable constraints on SIDM and CDM models.

The past two decades have seen a tremendous investment in observational facilities that promise to reveal new and unprecedented discoveries about the universe. In comparison, the investment in theoretical work is completely dwarfed, even though theory plays a crucial role in the interpretation of these observations, predicting new types of phenomena, and informing observing strategies. In this white paper, we argue that in order to reach the promised critical breakthroughs in astrophysics over the next decade and well beyond, the national agencies must take a serious approach to investment in theoretical astrophysics research. We discuss the role of theory in shaping our understanding of the universe, and then we provide a multi-level strategy, from the grassroots to the national, to address the current underinvestment in theory relative to observational work.

We study stellar-halo formation using six Milky Way-mass galaxies in FIRE-2 cosmological zoom simulations. We find that $5-40\%$ of the outer ($50-300$ kpc) stellar halo in each system consists of $\textit{in-situ}$ stars that were born in outflows from the main galaxy. Outflow stars originate from gas accelerated by super-bubble winds, which can be compressed, cool, and form co-moving stars. The majority of these stars remain bound to the halo and fall back with orbital properties similar to the rest of the stellar halo at $z=0$.In the outer halo, outflow stars are more spatially homogeneous, metal rich, and alpha-element-enhanced than the accreted stellar halo. At the solar location, up to $\sim 10 \%$ of our kinematically-identified halo stars were born in outflows; the fraction rises to as high as $\sim 40\%$ for the most metal-rich local halo stars ([Fe/H] $> -0.5$). We conclude that the Milky Way stellar halo could contain local counterparts to stars that are observed to form in molecular outflows in distant galaxies. Searches for such a population may provide a new, near-field approach to constraining feedback and outflow physics. A stellar halo contribution from outflows is a phase-reversal of the classic halo formation scenario of Eggen, Lynden-Bell $\&$ Sandange, who suggested that halo stars formed in rapidly $\textit{infalling}$ gas clouds. Stellar outflows may be observable in direct imaging of external galaxies and could provide a source for metal-rich, extreme velocity stars in the Milky Way.

We present a suite of FIRE-2 cosmological zoom-in simulations of isolated field dwarf galaxies, all with masses of $M_\mathrm{halo} \approx 10^{10}\,$M$_\odot$ at $z=0$, across a range of dark matter models. For the first time, we compare how both self-interacting dark matter (SIDM) and/or warm dark matter (WDM) models affect the assembly histories as well as the central density structure in fully hydrodynamical simulations of dwarfs. Dwarfs with smaller stellar half-mass radii (r$_{1/2}<500$ pc) have lower $\sigma_\star/V_\mathrm{max}$ ratios, reinforcing the idea that smaller dwarfs may reside in halos that are more massive than is naively expected. The majority of dwarfs simulated with self-interactions actually experience contraction of their inner density profiles with the addition of baryons relative to the cores produced in dark-matter-only runs, though the simulated dwarfs are always less centrally dense than in $\Lambda$CDM. The V$_{1/2}-$r$_{1/2}$ relation across all simulations is generally consistent with observations of Local Field dwarfs, though compact objects such as Tucana provide a unique challenge. Spatially-resolved rotation curves in the central regions ($<400$ pc) of small dwarfs could provide a way to distinguish between CDM, WDM, and SIDM, however: at the masses probed in this simulation suite, cored density profiles in dwarfs with small r$_{1/2}$ values can only originate from dark matter self-interactions.

We use a dedicated 0.7-m telescope to image the halos of 119 galaxies in the Local Volume to $\mu_r \sim 28-30$ mag/arcsec$^2$. The sample is primarily from the 2MASS Large Galaxy Atlas and extended to include nearby dwarf galaxies and more distant giant ellipticals, and spans fully the galaxy colour-magnitude diagram including the blue cloud and red sequence. We present an initial overview, including deep images of our galaxies. Our observations reproduce previously reported low surface brightness structures, including extended plumes in M51, and a newly discovered tidally extended dwarf galaxy in NGC7331. Low surface brightness structures, or "envelopes", exceeding 50 kpc in diameter are found mostly in galaxies with $M_V<-20.5$, and classic interaction signatures are infrequent. Defining a halo diameter at the surface brightness 28 mag/arcsec$^2$, we find that halo diameter is correlated with total galaxy luminosity. Extended signatures of interaction are found throughout the galaxy colour-magnitude diagram without preference for the red or blue sequences, or the green valley. Large envelopes may be found throughout the colour-magnitude diagram with some preference for the bright end of the red sequence. Spiral and S0 galaxies have broadly similar sizes, but ellipticals extend to notably greater diameters, reaching 150 kpc. We propose that the extended envelopes of disk galaxies are dominated by an extension of the disk population rather than by a classical population II halo.

We present the star formation history (SFH) of the isolated (D~970 kpc) Local Group dwarf galaxy WLM measured from color-magnitude diagrams constructed from deep Hubble Space Telescope imaging. Our observations include a central (0.5 $r_h$) and outer field (0.7 $ r_h$) that reach below the oldest main sequence turnoff. WLM has no early dominant episode of star formation: 20% of its stellar mass formed by ~12.5 Gyr ago (z~5). It also has an SFR that rises to the present with 50% of the stellar mass within the most recent 5 Gyr (z<0.7). There is evidence of a strong age gradient: the mean age of the outer field is 5 Gyr older than the inner field despite being only 0.4 kpc apart. Some models suggest such steep gradients are associated with strong stellar feedback and dark matter core creation. The SFHs of real isolated dwarf galaxies and those from the the Feedback In Realistic Environment suite are in good agreement for $M_{\star}(z=0) \sim 10^7-10^9 M_{\odot}$, but in worse agreement at lower masses ($M_{\star}(z=0) \sim 10^5-10^7 M_{\odot}$). These differences may be explainable by systematics in the models (e.g., reionization model) and/or observations (HST field placement). We suggest that a coordinated effort to get deep CMDs between HST/JWST (crowded central fields) and WFIRST (wide-area halo coverage) is the optimal path for measuring global SFHs of isolated dwarf galaxies.

The dark energy plus cold dark matter ($\Lambda$CDM) cosmological model has been a demonstrably successful framework for predicting and explaining the large-scale structure of Universe and its evolution with time. Yet on length scales smaller than $\sim 1$ Mpc and mass scales smaller than $\sim 10^{11} M_{\odot}$, the theory faces a number of challenges. For example, the observed cores of many dark-matter dominated galaxies are both less dense and less cuspy than naively predicted in $\Lambda$CDM. The number of small galaxies and dwarf satellites in the Local Group is also far below the predicted count of low-mass dark matter halos and subhalos within similar volumes. These issues underlie the most well-documented problems with $\Lambda$CDM: Cusp/Core, Missing Satellites, and Too-Big-to-Fail. The key question is whether a better understanding of baryon physics, dark matter physics, or both will be required to meet these challenges. Other anomalies, including the observed planar and orbital configurations of Local Group satellites and the tight baryonic/dark matter scaling relations obeyed by the galaxy population, have been less thoroughly explored in the context of $\Lambda$CDM theory. Future surveys to discover faint, distant dwarf galaxies and to precisely measure their masses and density structure hold promising avenues for testing possible solutions to the small-scale challenges going forward. Observational programs to constrain or discover and characterize the number of truly dark low-mass halos are among the most important, and achievable, goals in this field over then next decade. These efforts will either further verify the $\Lambda$CDM paradigm or demand a substantial revision in our understanding of the nature of dark matter.

Within $\Lambda$CDM, dwarf galaxies like the Large Magellanic Cloud (LMC) are expected to host numerous dark matter subhalos, several of which should host faint dwarf companions. Recent Gaia proper motions confirm new members of the LMC-system in addition to the previously known SMC, including two classical dwarf galaxies ($M_\ast$ > $10^5$ M$_{\odot}$; Carina and Fornax) as well as several ultra-faint dwarfs (Car2, Car3, Hor1, and Hyd1). We use the Feedback In Realistic Environments (FIRE) simulations to study the dark and luminous (down to ultrafaint masses, $M_\ast$ ~ $6\times10^3$ M$_{\odot}$) substructure population of isolated LMC-mass hosts ($M_\text{200}$ = $1-3\times10^{11}$ M$_{\odot}$) and place the Gaia + DES results in a cosmological context. By comparing number counts of subhalos in simulations with and without baryons, we find that, within 0.2 $r_\text{200}$, LMC-mass hosts deplete ~30% of their substructure, significantly lower than the ~70% of substructure depleted by Milky Way (MW) mass hosts. For our highest resolution runs ($m_\text{bary}$ = 880 M$_{\odot}$), ~5-10 subhalos form galaxies with $M_\ast$ > $10^4$ M$_{\odot}$, in agreement with the 7 observationally inferred pre-infall LMC companions. However, we find steeper simulated luminosity functions than observed, hinting at observation incompleteness at the faint end. The predicted DM content for classical satellites in FIRE agrees with observed estimates for Carina and Fornax, supporting the case for an LMC association. We predict that tidal stripping within the LMC potential lowers the inner dark matter density of ultra faint companions of the LMC. Thus, in addition to their orbital consistency, the low densities of dwarfs Car2, Hyd1, and Hyd2 reinforce their likelihood of Magellanic association.

The Maunakea Spectroscopic Explorer is a next-generation massively multiplexed spectroscopic facility currently under development in Hawaii. It is completely dedicated to large-scale spectroscopic surveys and will enable transformative science. In this white paper we summarize the science case and describe the current state of the project.

Observations of low-mass satellite galaxies in the nearby Universe point towards a strong dichotomy in their star-forming properties relative to systems with similar mass in the field. Specifically, satellite galaxies are preferentially gas poor and no longer forming stars, while their field counterparts are largely gas rich and actively forming stars. Much of the recent work to understand this dichotomy has been statistical in nature, determining not just that environmental processes are most likely responsible for quenching these low-mass systems but also that they must operate very quickly after infall onto the host system, with quenching timescales $\lesssim 2~ {\rm Gyr}$ at ${M}_{\star} \lesssim 10^{8}~{\rm M}_{\odot}$. This work utilizes the newly-available $Gaia$ DR2 proper motion measurements along with the Phat ELVIS suite of high-resolution, cosmological, zoom-in simulations to study low-mass satellite quenching around the Milky Way on an object-by-object basis. We derive constraints on the infall times for $37$ of the known low-mass satellite galaxies of the Milky Way, finding that $\gtrsim~70\%$ of the `classical' satellites of the Milky Way are consistent with the very short quenching timescales inferred from the total population in previous works. The remaining classical Milky Way satellites have quenching timescales noticeably longer, with $\tau_{\rm quench} \sim 6 - 8~{\rm Gyr}$, highlighting how detailed orbital modeling is likely necessary to understand the specifics of environmental quenching for individual satellite galaxies. Additionally, we find that the $6$ ultra-faint dwarf galaxies with publicly available $HST$-based star-formation histories are all consistent with having their star formation shut down prior to infall onto the Milky Way -- which, combined with their very early quenching times, strongly favors quenching driven by reionization.

We introduce an extension of the ELVIS project to account for the effects of the Milky Way galaxy on its subhalo population. Our simulation suite, Phat ELVIS, consists of twelve high-resolution cosmological dark matter-only (DMO) zoom simulations of Milky Way-size $\Lambda$CDM~ haloes ($M_{\rm v} = 0.7-2 \times 10^{12} \,\mathrm{M}_\odot$) along with twelve re-runs with embedded galaxy potentials grown to match the observed Milky Way disk and bulge today. The central galaxy potential destroys subhalos on orbits with small pericenters in every halo, regardless of the ratio of galaxy mass to halo mass. This has several important implications. 1) Most of the $\mathtt{Disk}$ runs have no subhaloes larger than $V_{\rm max} = 4.5$ km s$^{-1}$ within $20$ kpc and a significant lack of substructure going back $\sim 8$ Gyr, suggesting that local stream-heating signals from dark substructure will be rare. 2) The pericenter distributions of Milky Way satellites derived from $\mathit{Gaia}$ data are remarkably similar to the pericenter distributions of subhaloes in the $\mathtt{Disk}$ runs, while the DMO runs drastically over-predict galaxies with pericenters smaller than 20 kpc. 3) The enhanced destruction produces a tension opposite to that of the classic `missing satellites' problem: in order to account for ultra-faint galaxies known within $30$ kpc of the Galaxy, we must populate haloes with $V_\mathrm{peak} \simeq 7$ km s$^{-1}$ ($M \simeq 3 \times 10^{7} \,\mathrm{M}_\odot$ at infall), well below the atomic cooling limit of $V_\mathrm{peak} \simeq 16$ km s$^{-1}$ ($M \simeq 5 \times 10^{8} \,\mathrm{M}_\odot$ at infall). 4) If such tiny haloes do host ultra-faint dwarfs, this implies the existence of $\sim 1000$ satellite galaxies within 300 kpc of the Milky Way.

The Great Andromeda Galaxy (M31) is the nexus of the near-far galaxy evolution connection and a principal data point for near-field cosmology. Due to its proximity (780 kpc), M31 can be resolved into individual stars like the Milky Way (MW). Unlike the MW, we have the advantage of a global view of M31, enabling M31 to be observed with techniques that also apply to more distant galaxies. Moreover, recent evidence suggests that M31 may have survived a major merger within the last several Gyr, shaping the morphology of its stellar halo and triggering a starburst, while leaving the stellar disk largely intact. The MW and M31 thus provide complementary opportunities for in-depth studies of the disks, halos, and satellites of L* galaxies. Our understanding of the M31 system will be transformed in the 2020s if they include wide field facilities for both photometry (HST-like sensitivity and resolution) and spectroscopy (10-m class telescope, >1 sq. deg. field, highly multiplexed, R~ 3000 to 6000). We focus here on the power of these facilities to constrain the past, present, and future merger history of M31, via chemo-dynamical analyses and star formation histories of phase-mixed stars accreted at early times, as well as stars in surviving tidal debris features, M31's extended disk, and intact satellite galaxies that will eventually be tidally incorporated into the halo. This will yield an unprecedented view of the hierarchical formation of the M31 system and the subhalos that built it into the L* galaxy we observe today.

The Milky Way, the Andromeda galaxy, and Centaurus A host flattened distributions of satellite galaxies which exhibits coherent velocity trends indicative of rotation. Comparably extreme satellite structures are very rare in cosmological LCDM simulations, giving rise to the `satellite plane problem'. As a possible explanation it has been suggested that earlier-forming, higher concentration host halos contain more flattened and kinematically coherent satellite planes. We have tested for such a proposed correlation between the satellite plane and host halo properties in the ELVIS suite of simulations. We find evidence neither for a correlation of plane flattening with halo concentration or formation time, nor for a correlation of kinematic coherence with concentration. The height of the thinnest sub-halo planes does correlate with the host virial radius and with the radial extent of the sub-halo system. This can be understood as an effect of not accounting for differences in the radial distribution of sub-halos, and selecting them from different volumes than covered by the actual observations. Being part of a halo pair like the Local Group does not result in more narrow or more correlated satellite planes either. Additionally, using the PhatELVIS simulations we show that the presence of a central galaxy potential does not favor more narrow or more correlated satellite planes, it rather leads to slightly wider planes. Such a central potential is a good approximation of the dominant effect baryonic physics in cosmological simulations has on a sub-halo population. This suggests that, in contrast to other small-scale problems, the planes of satellite galaxies issue is made worse by accounting for baryonic effects.

We study star formation histories (SFHs) of $\simeq500$ dwarf galaxies (stellar mass $M_\ast = 10^5 - 10^9\,M_\odot$) from FIRE-2 cosmological zoom-in simulations. We compare dwarfs around individual Milky Way (MW)-mass galaxies, dwarfs in Local Group (LG)-like environments, and true field (i.e. isolated) dwarf galaxies. We reproduce observed trends wherein higher-mass dwarfs quench later (if at all), regardless of environment. We also identify differences between the environments, both in terms of "satellite vs. central" and "LG vs. individual MWvs. isolated dwarf central." Around the individual MW-mass hosts, we recover the result expected from environmental quenching: central galaxies in the "near field" have more extended SFHs than their satellite counterparts, with the former more closely resemble isolated ("true field") dwarfs (though near-field centrals are still somewhat earlier forming). However, this difference is muted in the LG-like environments, where both near-field centrals and satellites have similar SFHs, which resemble satellites of single MW-mass hosts. This distinction is strongest for $M_\ast = 10^6 - 10^7\,M_\odot$ but exists at other masses. Our results suggest that the paired halo nature of the LG may regulate star formation in dwarf galaxies even beyond the virial radii of the MW and Andromeda. Caution is needed when comparing zoom-in simulations targeting isolated dwarf galaxies against observed dwarf galaxies in the LG.

Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution.

We perform high-resolution simulations of a MW-like galaxy in a self-interacting cold dark matter model with elastic cross section over mass of $1~\rm cm^2/g$ (SIDM) and compare to a model without self-interactions (CDM). We run our simulations with and without a time-dependent embedded potential to capture effects of the baryonic disk and bulge contributions. The CDM and SIDM simulations with the embedded baryonic potential exhibit remarkably similar host halo profiles, subhalo abundances and radial distributions within the virial radius. The SIDM host halo is denser in the center than the CDM host and has no discernible core, in sharp contrast to the case without the baryonic potential (core size $\sim 7 \, \rm kpc$). The most massive subhalos (with $V_{\mathrm{peak}}> 20 \, \rm km/s$) in our SIDM simulations, expected to host the classical satellite galaxies, have density profiles that are less dense than their CDM analogs at radii less than 500 pc but the deviation diminishes for less massive subhalos. With the baryonic potential included in the CDM and SIDM simulations, the most massive subhalos do not display the too-big-to-fail problem. However, the least dense among the massive subhalos in both these simulations tend to have the smallest pericenter values, a trend that is not apparent among the bright MW satellite galaxies.

We explore the radial variation of star formation histories in dwarf galaxies simulated with Feedback In Realistic Environments (FIRE) physics. The sample contains 9 low-mass field dwarfs with M_ star = 10^5 - 10^7 M_sun from previous FIRE results, and a new suite of 17 higher mass field dwarfs with M_star = 10^7 - 10^9 M_sun introduced here. We find that age gradients are common in our dwarfs, with older stars dominant at large radii. The strength of the gradient correlates with overall galaxy age such that earlier star formation produces a more pronounced gradient. The relation between formation time and strength of the gradient is driven by both mergers and star-formation feedback. Mergers can both steepen and flatten the age gradient depending on the timing of the merger and star formation history of the merging galaxy. In galaxies without significant mergers, early feedback pushes stars to the outskirts at early times. Interestingly, among galaxies without mergers, those with large dark matter cores have flatter age gradients because these galaxies have more late-time feedback. If real galaxies have age gradients as we predict, stellar population studies that rely on sampling a limited fraction of a galaxy can give a biased view of its global star formation history. We show that central fields can be biased young by a few Gyrs while outer fields are biased old. Fields positioned near the 2D half-light radius will provide the least biased measure of a dwarf galaxy's global star formation history.

The predominantly ancient stellar populations observed in the lowest-mass galaxies (i.e. ultra-faint dwarfs) suggest that their star formation was suppressed by reionization. Most of the well-studied ultra-faint dwarfs, however, are within the central half of the Milky Way dark matter halo, such that they are consistent with a population that was accreted at early times and thus potentially quenched via environmental processes. To study the potential role of environment in suppressing star formation on the smallest scales, we utilize the Exploring the Local Volume in Simulations (ELVIS) suite of $N$-body simulations to constrain the distribution of infall times for low-mass subhalos likely to host the ultra-faint population. For the ultra-faint satellites of the Milky Way with star-formation histories inferred from $Hubble~Space~Telescope$ imaging, we find that environment is highly unlikely to play a dominant role in quenching their star formation. Even when including the potential effects of pre-processing, there is a $\lesssim 0.1\%$ probability that environmental processes quenched all of the known ultra-faint dwarfs early enough to explain their observed star-formation histories. Instead, we argue for a mass floor in the effectiveness of satellite quenching at roughly $M_{\star} \sim 10^{5}~M_{\odot}$, below which star formation in surviving galaxies is globally suppressed by reionization. We predict a large population of quenched ultra-faint dwarfs in the Local Field ($1 < R/R_{\rm{vir}} < 2$), with as many as $\sim250$ to be discovered by future wide-field imaging surveys.

We study a suite of extremely high-resolution cosmological FIRE simulations of dwarf galaxies ($M_{\rm halo} \lesssim 10^{10}$$M_{\odot}$), run to $z=0$ with $30 M_{\odot}$ resolution, sufficient (for the first time) to resolve the internal structure of individual supernovae remnants within the cooling radius. Every halo with $M_{\rm halo} \gtrsim 10^{8.6} M_{\odot}$ is populated by a resolved {\em stellar} galaxy, suggesting very low-mass dwarfs may be ubiquitous in the field. Our ultra-faint dwarfs (UFDs; $M_{\ast}<10^{5}\,M_{\odot}$) have their star formation truncated early ($z\gtrsim2$), likely by reionization, while classical dwarfs ($M_{\ast}>10^{5} M_{\odot}$) continue forming stars to $z<0.5$. The systems have bursty star formation (SF) histories, forming most of their stars in periods of elevated SF strongly clustered in both space and time. This allows our dwarf with $M_{\ast}/M_{\rm halo} > 10^{-4}$ to form a dark matter core $>200$pc, while lower-mass UFDs exhibit cusps down to $\lesssim100$pc, as expected from energetic arguments. Our dwarfs with $M_{\ast}>10^{4}\,M_{\odot}$ have half-mass radii ($R_{\rm 1/2}$) in agreement with Local Group (LG) dwarfs; dynamical mass vs. $R_{1/2}$ and the degree of rotational support also resemble observations. The lowest-mass UFDs are below surface brightness limits of current surveys but are potentially visible in next-generation surveys (e.g. LSST). The stellar metallicities are lower than in LG dwarfs; this may reflect pre-enrichment of the LG by the massive hosts or Pop-III stars. Consistency with lower resolution studies implies that our simulations are numerically robust (for a given physical model).

The Feedback In Realistic Environments (FIRE) project explores feedback in cosmological galaxy formation simulations. Previous FIRE simulations used an identical source code (FIRE-1) for consistency. Motivated by the development of more accurate numerics - including hydrodynamic solvers, gravitational softening, and supernova coupling algorithms - and exploration of new physics (e.g. magnetic fields), we introduce FIRE-2, an updated numerical implementation of FIRE physics for the GIZMO code. We run a suite of simulations and compare against FIRE-1: overall, FIRE-2 improvements do not qualitatively change galaxy-scale properties. We pursue an extensive study of numerics versus physics. Details of the star-formation algorithm, cooling physics, and chemistry have weak effects, provided that we include metal-line cooling and star formation occurs at higher-than-mean densities. We present new resolution criteria for high-resolution galaxy simulations. Most galaxy-scale properties are robust to numerics we test, provided: (1) Toomre masses are resolved; (2) feedback coupling ensures conservation, and (3) individual supernovae are time-resolved. Stellar masses and profiles are most robust to resolution, followed by metal abundances and morphologies, followed by properties of winds and circum-galactic media (CGM). Central (~kpc) mass concentrations in massive (L*) galaxies are sensitive to numerics (via trapping/recycling of winds in hot halos). Multiple feedback mechanisms play key roles: supernovae regulate stellar masses/winds; stellar mass-loss fuels late star formation; radiative feedback suppresses accretion onto dwarfs and instantaneous star formation in disks. We provide all initial conditions and numerical algorithms used.

A standard prediction of galaxy formation theory is that the ionizing background suppresses galaxy formation in haloes with peak circular velocities smaller than Vpeak ~ 20 km/s, rendering the majority of haloes below this scale completely dark. We use a suite of cosmological zoom simulations of Milky Way-like haloes that include central Milky Way disk galaxy potentials to investigate the relationship between subhaloes and ultrafaint galaxies. We find that there are far too few subhaloes within 50 kpc of the Milky Way that had Vpeak > 20 km/s to account for the number of ultrafaint galaxies already known within that volume today. In order to match the observed count, we must populate subhaloes down to Vpeak ~ 6 km/s with ultrafaint dwarfs. The required haloes have peak virial temperatures as low as 1,500 K, well below the atomic hydrogen cooling limit of 10^4 K. Allowing for the possibility that the Large Magellanic Cloud contributes several of the satellites within 50 kpc could potentially raise this threshold to 10 km/s (4,000 K), still below the atomic cooling limit and far below the nominal reionization threshold.

Binary black holes are the primary endpoint of massive stellar evolution. Their properties provide a unique opportunity to constrain binary evolution, which is still poorly understood. In this paper, we predict the inventory of binary black holes and their merger products in/around the Milky Way, and detail their main properties. We present the first combination of a high-resolution cosmological simulation of a Milky Way-mass galaxy with a binary population synthesis model. The hydrodynamic simulation, taken from the FIRE project, provides a cosmologically realistic star formation history for the galaxy and its stellar halo and satellites. We apply a metallicity-dependent evolutionary model to the star particles to produce individual binary black holes. We find that a million binary black holes have merged in the model Milky Way, and 3 million binaries are still present, with an average mass of 28 Msun per binary. Because the black hole progenitors are biased towards low metallicity stars, half reside in the stellar halo and satellites and 40 per cent of the binaries were formed outside the main galaxy. This trend increases with the masses of the black holes. The numbers and mass distribution of the merged systems is compatible with the LIGO/Virgo detections. Observations of these black holes will be challenging, both with electromagnetic methods and LISA. We find that a cosmologically realistic star formation history, with self-consistent metal enrichment and Galactic accretion history, are key ingredients for determining binary black hole rates that can be compared with observations to constrain massive binary evolution.

Building upon results of cosmological simulations of ultra-light scalar field dark matter (SFDM), we present a comprehensive model for the density profiles of SFDM haloes as a function of halo virial mass $M_{\rm h}$ and scalar field mass $m$. The central regions of SFDM haloes are dominated by solitons with characteristic densities that increase with increasing halo mass and asymptote to CDM-like profiles at large radii. For scalar field masses $m \sim 10^{-22}$ eV, consistent with large-scale structure observations, $M_{\rm h} \sim 10^{10} \,M_\odot$ haloes have lower core densities than their Cold Dark Matter (CDM) counterparts and this alleviates the Too Big to Fail problem (TBTF) in a regime where feedback is less effective. However, higher-mass SFDM haloes with $M_{\rm h} \sim 10^{11} \,M_\odot$ are denser than their CDM counterparts at small, observationally relevant radii. We use rotation curves of $V \sim 100$ km s$^{-1}$ galaxies from the SPARC database to show that SFDM exacerbates the cusp/core and central density problems seen in CDM at this scale. We conclude that if the conventional cosmological SFDM scaling relations are correct, then baryonic feedback is required to lower densities in SFDM haloes even more so than in CDM. This motivates cosmological and self-consistent hydrodynamic simulations of SFDM to determine whether central soliton structure can be altered by realistic feedback implementations.

We present a new set of high-resolution hydrodynamic cosmological zoom-in simulations that apply the Feedback In Realistic Environments (FIRE) physics to both Local Group (LG)-like and isolated Milky Way (MW)-like volumes (ten host systems in total with baryonic particle mass $\simeq 3,500-7,000\,M_\odot$). We study the stellar mass functions, circular velocity or mass profiles, and velocity dispersions of the dwarf galaxy populations. The simulations reproduce the stellar mass function and central densities of MW satellite dwarfs for $M_\ast \geq10^{5.5}\,M_\odot$ and predict the existence of $\sim3$ unidentified galaxies with $M_\ast\sim10^5\,M_\odot$ within $300$ kpc of the MW. Overall, we find no evidence for the classical missing satellites or too-big-to-fail (TBTF) problems for satellite galaxies in our sample. Among the satellites, TBTF is resolved primarily by subhalo disruption and overall mass loss; central density profiles of subhalos are of secondary importance. For non-satellite galaxies, our LG-like simulations predict as many as $\sim10$ as-of-yet unseen galaxies at distances $0.3-1$ Mpc from both hosts, with $M_\ast\simeq10^{5-6}\,M_\odot$ (in halos with $V_\mathrm{max}\sim 20~\mathrm{km\,s}^{-1}$), albeit with large halo-to-halo variance. None of our simulations produces a compact, baryon-dominated, high-density dwarf elliptical-type galaxy (with $V_\mathrm{circ} \gtrsim 35~\mathrm{km\,s}^{-1}$ at $r<1\,$kpc), of which six may appear in the LG (but none in the MW). It may therefore remain a challenge to reproduce the full diversity of the dwarf population, including both the highest and lowest density systems.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

In the local Universe, there is a strong division in the star-forming properties of low-mass galaxies, with star formation largely ubiquitous amongst the field population while satellite systems are predominantly quenched. This dichotomy implies that environmental processes play the dominant role in suppressing star formation within this low-mass regime (${M}_{\star} \sim 10^{5.5-8}~{\rm M}_{\odot}$). As shown by observations of the Local Volume, however, there is a non-negligible population of passive systems in the field, which challenges our understanding of quenching at low masses. By applying the satellite quenching models of Fillingham et al. (2015) to subhalo populations in the Exploring the Local Volume In Simulations (ELVIS) suite, we investigate the role of environmental processes in quenching star formation within the nearby field. Using model parameters that reproduce the satellite quenched fraction in the Local Group, we predict a quenched fraction -- due solely to environmental effects -- of $\sim 0.52 \pm 0.26$ within $1< R/R_{\rm vir} < 2$ of the Milky Way and M31. This is in good agreement with current observations of the Local Volume and suggests that the majority of the passive field systems observed at these distances are quenched via environmental mechanisms. Beyond $2~R_{\rm vir}$, however, dwarf galaxy quenching becomes difficult to explain through an interaction with either the Milky Way or M31, such that more isolated, field dwarfs may be self-quenched as a result of star-formation feedback.

We study the impact of a warm dark matter (WDM) cosmology on dwarf galaxy formation through a suite of cosmological hydrodynamical zoom-in simulations of $M_{\rm halo} \approx10^{10}\,M_{\odot}$ dark matter halos as part of the Feedback in Realistic Environments (FIRE) project. A main focus of this paper is to evaluate the combined effects of dark matter physics and stellar feedback on the well-known small-scale issues found in cold dark matter (CDM) models. We find that the $z=0$ stellar mass of a galaxy is strongly correlated with the central density of its host dark matter halo at the time of formation, $z_{\rm f}$, in both CDM and WDM models. WDM halos follow the same $M_{\star}(z=0)-V_{\rm max}(z_{\rm f})$ relation as in CDM, but they form later, are less centrally dense, and therefore contain galaxies that are less massive than their CDM counterparts. As a result, the impact of baryonic effects on the central gravitational potential is typically diminished relative to CDM. However, the combination of delayed formation in WDM and energy input from stellar feedback results in dark matter profiles with lower overall densities. The WDM galaxies studied here have a wider diversity of star formation histories (SFHs) than the same systems simulated in CDM, and the two lowest $M_{\star}$ WDM galaxies form all of their stars at late times. The discovery of young ultra-faint dwarf galaxies with no ancient star formation -- which do not exist in our CDM simulations -- would therefore provide evidence in support of WDM.

The existence of stellar halos around galaxies is a natural prediction of the hierarchical nature of the LambdaCDM model. Recent observations of Milky Way-like galaxies have revealed a wide range in stellar halo mass, including cases with no statistically significant detection of a stellar halo, as in the case of M101, NGC3351 and NGC1042. We use the Illustris simulation to investigate the scatter in stellar halo content and, in particular, to study the formation of galaxies with the smallest fraction of this diffuse component. Stellar halos are far from spherical, which diminishes the surface brightness of the stellar halo for face-on disks. Once accounting for projection effects, we find that the stellar halo fraction f_SH correlates strongly with galaxy morphology and star formation rate, but not with environment, in agreement with observations. Galaxies with the lowest stellar halo fractions are disk-dominated, star-forming and assemble their dark matter halos on average earlier than galaxies with similar stellar masses. Accreted satellites are also lower in stellar mass and have earlier infall times than centrals with high f_SH. In situ rather than accreted stars dominate the stellar halos of galaxies with the lowest stellar halo fractions, with a transition radius from in situ to accretion-dominated r~45 kpc. Our results extrapolated to real galaxies such as M101 may indicate that these galaxies inhabit old halos which endured mergers only at higher redshifts and evolved relatively unperturbed in the last ~10 Gyrs.

We investigate the merger histories of isolated dwarf galaxies based on a suite of 15 high-resolution cosmological zoom-in simulations, all with masses of $M_{\rm halo} \approx 10^{10}\,{\rm M}_{\odot}$ (and M$_\star\sim10^5-10^7\,{\rm M}_{\odot}$) at $z=0$, from the Feedback in Realistic Environments (FIRE) project. The stellar populations of these dwarf galaxies at $z=0$ are formed essentially entirely "in situ": over 90$\%$ of the stellar mass is formed in the main progenitor in all but two cases, and all 15 of the galaxies have >70$\%$ of their stellar mass formed in situ. Virtually all galaxy mergers occur prior to $z\sim3$, meaning that accreted stellar populations are ancient. On average, our simulated dwarfs undergo 5 galaxy mergers in their lifetimes, with typical pre-merger galaxy mass ratios that are less than 1:10. This merger frequency is generally comparable to what has been found in dissipationless simulations when coupled with abundance matching. Two of the simulated dwarfs have a luminous satellite companion at $z=0$. These ultra-faint dwarfs lie at or below current detectability thresholds but are intriguing targets for next-generation facilities. The small contribution of accreted stars make it extremely difficult to discern the effects of mergers in the vast majority of dwarfs either photometrically or using resolved-star color-magnitude diagrams (CMDs). The important implication for near-field cosmology is that star formation histories of comparably massive galaxies derived from resolved CMDs should trace the build-up of stellar mass in one main system across cosmic time as opposed to reflecting the contributions of many individual star formation histories of merged dwarfs.

We use hydrodynamic cosmological zoom-in simulations from the FIRE project to explore the morphologies and kinematics of fifteen Milky Way (MW)-mass galaxies. Our sample ranges from compact, bulge-dominated systems with 90% of their stellar mass within 2.5 kpc to well-ordered disks that reach $\gtrsim15$ kpc. The gas in our galaxies always forms a thin, rotation-supported disk at $z=0$, with sizes primarily determined by the gas mass. For stars, we quantify kinematics and morphology both via the fraction of stars on disk-like orbits and with the radial extent of the stellar disk. In this mass range, stellar morphology and kinematics are poorly correlated with the properties of the halo available from dark matter-only simulations (halo merger history, spin, or formation time). They more strongly correlate with the gaseous histories of the galaxies: those that maintain a high gas mass in the disk after $z\sim1$ develop well-ordered stellar disks. The best predictor of morphology we identify is the spin of the gas in the halo at the time the galaxy formed 1/2 of its stars (i.e. the gas that builds the galaxy). High-$z$ mergers, before a hot halo emerges, produce some of the most massive bulges in the sample (from compact disks in gas-rich mergers), while later-forming bulges typically originate from internal processes, as satellites are stripped of gas before the galaxies merge. Moreover, most stars in $z=0$ MW-mass galaxies (even $z=0$ bulge stars) form in a disk: $\gtrsim$60-90% of stars begin their lives rotationally supported.

We present the velocity dispersion of red giant branch (RGB) stars in M31's halo, derived by modeling the line of sight velocity distribution of over 5000 stars in 50 fields spread throughout M31's stellar halo. The dataset was obtained as part of the SPLASH (Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo) Survey, and covers projected radii of 9 to 175 kpc from M31's center. All major structural components along the line of sight in both the Milky Way (MW) and M31 are incorporated in a Gaussian Mixture Model, including all previously identified M31 tidal debris features in the observed fields. The probability an individual star is a constituent of M31 or the MW, based on a set of empirical photometric and spectroscopic diagnostics, is included as a prior probability in the mixture model. The velocity dispersion of stars in M31's halo is found to decrease only mildly with projected radius, from 108 km/s in the innermost radial bin (8.2 to 14.1 kpc) to $\sim 80$ to 90 km/s at projected radii of $\sim 40$ to 130 kpc, and can be parameterized with a power-law of slope $-0.12\pm 0.05$. The quoted uncertainty on the power-law slope reflects only the precision of the method, although other sources of uncertainty we consider contribute negligibly to the overall error budget.

A fundamental prediction of the cold dark matter cosmology is the existence of a large number of dark subhalos around galaxies, most of which should be entirely devoid of stars. Confirming the existence of dark substructures stands among the most important empirical challenges in modern cosmology: if they are found and quantified with the mass spectrum expected, then this would close the door on a vast array of competing theories. But in order for observational programs of this kind to reach fruition, we need robust predictions. Here we explore substructure predictions for lensing using galaxy lens-like hosts at z=0.2 from the Illustris simulations both in full hydrodynamics and dark matter only. We quantify substructures more massive than ~ 10^9 M_sun, comparable to current lensing detections derived from HST, Keck, and ALMA. The addition of full hydrodynamics reduces the overall subhalo mass function by about a factor of two. Even for the dark matter only runs, most (~ 85%) lines of sight through projected cylinders of size close to an Einstein radius contain no substructures larger than 10^9 M_sun. The fraction of empty sight lines rises to ~ 95% in full physics simulations. This suggests we will likely need hundreds of strong lensing systems suitable for substructure studies, as well as predictions that include the effects of baryon physics on substructure, to properly constrain cosmological models. Fortunately, the field is poised to fulfill these requirements.

The spatial distribution of satellite galaxies around pairs of galaxies in the Sloan Digital Sky Survey (SDSS) have been found to bulge significantly towards the respective partner. Highly anisotropic, planar distributions of satellite galaxies are in conflict with expectations derived from cosmological simulations. Does the lopsided distribution of satellite systems around host galaxy pairs constitute a similar challenge to the standard model of cosmology? We investigate whether such satellite distributions are present around stacked pairs of hosts extracted from the $\Lambda$CDM simulations Millennium-I, Millennium-II, ELVIS, and Illustris-1. By utilizing this set of simulations covering different volumes, resolutions, and physics, we implicitly test whether a lopsided signal exists for different ranges of satellite galaxy masses, and whether the inclusion of hydrodynamical effects produces significantly different results. All simulations display a lopsidedness similar to the observed situation. The signal is highly significant for simulations containing a sufficient number of hosts and resolved satellite galaxies (up to $5\,\sigma$ for Millennium-II). We find a projected signal that is up to twice as strong as that reported for the SDSS systems for certain opening angles ($\sim16\%$ more satellites in the direction between the pair than expected for uniform distributions). Considering that the SDSS signal is a lower limit owing to likely back- and foreground contamination, the $\Lambda$CDM simulations appear to be consistent with this particular empirical property of galaxy pairs.

We present an empirical approach for interpreting gravitational wave signals of binary black hole mergers under the assumption that the underlying black hole population is sourced by remnants of stellar evolution. Using the observed relationship between galaxy mass and stellar metallicity, we predict the black hole count as a function of galaxy stellar mass. We show, for example, that a galaxy like the Milky Way should host millions of $\sim 30~M_\odot$ black holes and dwarf satellite galaxies like Draco should host $\sim 100$ such remnants, with weak dependence on the assumed IMF and stellar evolution model. Most low-mass black holes ($\sim10 M_\odot$) typically reside within massive galaxies ($M_\star \simeq 10^{11} M_\odot$) while massive black holes ($\sim 50~M_\odot$) typically reside within dwarf galaxies ($M_\odot \simeq 10^9 M_\odot$) today. If roughly $1\%$ of black holes are involved in a binary black hole merger, then the reported merger rate densities from Advanced LIGO can be accommodated for a range of merger timescales, and the detection of mergers with $> 50~M_\odot$ black holes should be expected within the next decade. Identifying the host galaxy population of the mergers provides a way to constrain both the binary neutron star or black hole formation efficiencies and the merger timescale distributions; these events would be primarily localized in dwarf galaxies if the merger timescale is short compared to the age of the universe and in massive galaxies otherwise. As more mergers are detected, the prospect of identifying the host galaxy population, either directly through the detection of electromagnetic counterparts of binary neutron star mergers or indirectly through the anisotropy of the events, will become a realistic possibility.

We present new proper motion (PM) measurements of the dwarf spheroidal galaxies (dSphs) Draco and Sculptor using multi-epoch images obtained with the Hubble Space Telescope ACS/WFC. Our PM results have uncertainties far lower than previous measurements, even made with the same instrument. The PM results for Draco and Sculptor are (mu_W,mu_N)_Dra = (-0.0562+/-0.0099,-0.1765+/-0.0100) mas/yr and (mu_W,mu_N)_Scl = (-0.0296+/-0.0209,-0.1358 +/-0.0214) mas/yr. The implied Galactocentric velocity vectors for Draco and Sculptor have radial and tangential components: (V_rad,V_tan)_Dra = (-88.6,161.4) +/- (4.4,5.6) km/s; and (V_rad,V_tan)_Scl = (72.6,200.2) +/- (1.3,10.8) km/s. We study the detailed orbital history of both Draco and Sculptor via numerical orbit integrations. Orbital periods of Draco and Sculptor are found to be 1-2 and 2-5 Gyrs, respectively, accounting for uncertainties in the MW mass. We also study the influence of the Large Magellanic Cloud (LMC) on the orbits of Draco and Sculptor. Overall, the inclusion of the LMC increases the scatter in the orbital results. Based on our calculations, Draco shows a rather wide range of orbital parameters depending on the MW mass and inclusion/exclusion of the LMC, but Sculptor's orbit is very well constrained with its most recent pericentric approach to the MW being 0.3-0.4 Gyr ago. Our new PMs imply that the orbital trajectories of both Draco and Sculptor are confined within the Disk of Satellites (DoS), better so than implied by earlier PM measurements, and likely rule out the possibility that these two galaxies were accreted together as part of a tightly bound group.

We present a suite of 15 cosmological zoom-in simulations of isolated dark matter halos, all with masses of $M_{\rm halo} \approx 10^{10}\,{\rm M}_\odot$ at $z=0$, in order to understand the relationship between halo assembly, galaxy formation, and feedback's effects on the central density structure in dwarf galaxies. These simulations are part of the Feedback in Realistic Environments (FIRE) project and are performed at extremely high resolution. The resultant galaxies have stellar masses that are consistent with rough abundance matching estimates, coinciding with the faintest galaxies that can be seen beyond the virial radius of the Milky Way ($M_\star/{\rm M}_\odot\approx 10^5-10^7$). This non-negligible spread in stellar mass at $z=0$ in halos within a narrow range of virial masses is strongly correlated with central halo density or maximum circular velocity $V_{\rm max}$. Much of this dependence of $M_\star$ on a second parameter (beyond $M_{\rm halo}$) is a direct consequence of the $M_{\rm halo}\sim10^{10}\,{\rm M}_\odot$ mass scale coinciding with the threshold for strong reionization suppression: the densest, earliest-forming halos remain above the UV-suppression scale throughout their histories while late-forming systems fall below the UV-suppression scale over longer periods and form fewer stars as a result. In fact, the latest-forming, lowest-concentration halo in our suite fails to form any stars. Halos that form galaxies with $M_\star\gtrsim2\times10^{6}\,{\rm M}_\odot$ have reduced central densities relative to dark-matter-only simulations, and the radial extent of the density modifications is well-approximated by the galaxy half-mass radius $r_{1/2}$. This apparent stellar mass threshold of $M_\star \approx 2\times 10^{6} \approx 2\times 10^{-4} \,M_{\rm halo}$ is broadly consistent with previous work and provides a testable prediction of FIRE feedback models in LCDM.

We investigate angular momentum acquisition in Milky Way-sized galaxies by comparing five high resolution zoom-in simulations, each implementing identical cosmological initial conditions, but utilizing different hydrodynamic codes: Enzo, Art, Ramses, Arepo, and Gizmo-PSPH. Each code implements a distinct set of feedback and star formation prescriptions. We find that while many galaxy and halo properties vary between the different codes (and feedback prescriptions), there is qualitative agreement on the process of angular momentum acquisition in the galaxy's halo. In all simulations, cold filamentary gas accretion to the halo results in ~4 times more specific angular momentum in cold halo gas ($\lambda_{cold} \gtrsim 0.1$) than in the dark matter halo. At z>1, this inflow takes the form of inspiraling cold streams that are co-directional in the halo of the galaxy and are fueled, aligned, and kinematically connected to filamentary gas infall along the cosmic web. Due to the qualitative agreement among disparate simulations, we conclude that the buildup of high angular momentum halo gas and the presence of these inspiraling cold streams are robust predictions of Lambda Cold Dark Matter galaxy formation, though the detailed morphology of these streams is significantly less certain. A growing body of observational evidence suggests that this process is borne out in the real universe.

Among the most important goals in cosmology is detecting and quantifying small ($M_{\rm halo}\simeq10^{6-9}~\mathrm{M}_\odot$) dark matter (DM) subhalos. Current probes around the Milky Way (MW) are most sensitive to such substructure within $\sim20$ kpc of the halo center, where the galaxy contributes significantly to the potential. We explore the effects of baryons on subhalo populations in $\Lambda$CDM using cosmological zoom-in baryonic simulations of MW-mass halos from the Latte simulation suite, part of the Feedback In Realistic Environments (FIRE) project. Specifically, we compare simulations of the same two halos run using (1) DM-only (DMO), (2) full baryonic physics, and (3) DM with an embedded disk potential grown to match the FIRE simulation. Relative to baryonic simulations, DMO simulations contain $\sim2\times$ as many subhalos within 100 kpc of the halo center; this excess is $\gtrsim5\times$ within 25 kpc. At $z=0$, the baryonic simulations are completely devoid of subhalos down to $3\times10^6~\mathrm{M}_\odot$ within $15$ kpc of the MW-mass galaxy, and fewer than 20 surviving subhalos have orbital pericenters <20 kpc. Despite the complexities of baryonic physics, the simple addition of an embedded central disk potential to DMO simulations reproduces this subhalo depletion, including trends with radius, remarkably well. Thus, the additional tidal field from the central galaxy is the primary cause of subhalo depletion. Subhalos on radial orbits that pass close to the central galaxy are preferentially destroyed, causing the surviving subhalo population to have tangentially biased orbits compared to DMO predictions. Our method of embedding a disk potential in DMO simulations provides a fast and accurate alternative to full baryonic simulations, thus enabling suites of cosmological simulations that can provide accurate and statistical predictions of substructure populations.

We compare a suite of four simulated dwarf galaxies formed in 10$^{10} M_{\odot}$ haloes of collisionless Cold Dark Matter (CDM) with galaxies simulated in the same haloes with an identical galaxy formation model but a non-zero cross-section for dark matter self-interactions. These cosmological zoom-in simulations are part of the Feedback In Realistic Environments (FIRE) project and utilize the FIRE-2 model for hydrodynamics and galaxy formation physics. We find the stellar masses of the galaxies formed in Self-Interacting Dark Matter (SIDM) with $\sigma/m= 1\, cm^2/g$ are very similar to those in CDM (spanning $M_{\star} \approx 10^{5.7 - 7.0} M_{\odot}$) and all runs lie on a similar stellar mass -- size relation. The logarithmic dark matter density slope ($\alpha=d\log \rho / d\log r$) in the central $250-500$ pc remains steeper than $\alpha= -0.8$ for the CDM-Hydro simulations with stellar mass $M_{\star} \sim 10^{6.6} M_{\odot}$ and core-like in the most massive galaxy. In contrast, every SIDM hydrodynamic simulation yields a flatter profile, with $\alpha >-0.4$. Moreover, the central density profiles predicted in SIDM runs without baryons are similar to the SIDM runs that include FIRE-2 baryonic physics. Thus, SIDM appears to be much more robust to the inclusion of (potentially uncertain) baryonic physics than CDM on this mass scale, suggesting SIDM will be easier to falsify than CDM using low-mass galaxies. Our FIRE simulations predict that galaxies less massive than $M_{\star} < 3 \times 10^6 M_{\odot}$ provide potentially ideal targets for discriminating models, with SIDM producing substantial cores in such tiny galaxies and CDM producing cusps.

We use a suite of high-resolution cosmological dwarf galaxy simulations to test the accuracy of commonly-used mass estimators from Walker et al.(2009) and Wolf et al.(2010), both of which depend on the observed line-of-sight velocity dispersion and the 2D half-light radius of the galaxy, $Re$. The simulations are part of the the Feedback in Realistic Environments (FIRE) project and include twelve systems with stellar masses spanning $10^{5} - 10^{7} M_{\odot}$ that have structural and kinematic properties similar to those of observed dispersion-supported dwarfs. Both estimators are found to be quite accurate: $M_{Wolf}/M_{true} = 0.98^{+0.19}_{-0.12}$ and $M_{Walker}/M_{true} =1.07^{+0.21}_{-0.15}$, with errors reflecting the 68% range over all simulations. The excellent performance of these estimators is remarkable given that they each assume spherical symmetry, a supposition that is broken in our simulated galaxies. Though our dwarfs have negligible rotation support, their 3D stellar distributions are flattened, with short-to-long axis ratios $ c/a \simeq 0.4-0.7$. The accuracy of the estimators shows no trend with asphericity. Our simulated galaxies have sphericalized stellar profiles in 3D that follow a nearly universal form, one that transitions from a core at small radius to a steep fall-off $\propto r^{-4.2}$ at large $r$, they are well fit by S\'ersic profiles in projection. We find that the most important empirical quantity affecting mass estimator accuracy is $Re$ . Determining $Re$ by an analytic fit to the surface density profile produces a better estimated mass than if the half-light radius is determined via direct summation.

We use Local Group galaxy counts together with the ELVIS N-body simulations to explore the relationship between the scatter and slope in the stellar mass vs. halo mass relation at low masses, $M_\star \simeq 10^5 - 10^8 M_\odot$. Assuming models with log-normal scatter about a median relation of the form $M_\star \propto M_\mathrm{halo}^\alpha$, the preferred log-slope steepens from $\alpha \simeq 1.8$ in the limit of zero scatter to $\alpha \simeq 2.6$ in the case of $2$ dex of scatter in $M_\star$ at fixed halo mass. We provide fitting functions for the best-fit relations as a function of scatter, including cases where the relation becomes increasingly stochastic with decreasing mass. We show that if the scatter at fixed halo mass is large enough ($\gtrsim 1$ dex) and if the median relation is steep enough ($\alpha \gtrsim 2$), then the "too-big-to-fail" problem seen in the Local Group can be self-consistently eliminated in about $\sim 5-10\%$ of realizations. This scenario requires that the most massive subhalos host unobservable ultra-faint dwarfs fairly often; we discuss potentially observable signatures of these systems. Finally, we compare our derived constraints to recent high-resolution simulations of dwarf galaxy formation in the literature. Though simulation-to-simulation scatter in $M_\star$ at fixed $M_\mathrm{halo}$ is large among separate authors ($\sim 2$ dex), individual codes produce relations with much less scatter and usually give relations that would over-produce local galaxy counts.

We perform a systematic Bayesian analysis of rotation vs. dispersion support ($v_{\rm rot} / \sigma$) in $40$ dwarf galaxies throughout the Local Volume (LV) over a stellar mass range $10^{3.5} M_{\rm \odot} < M_{\star} < 10^8 M_{\rm \odot}$. We find that the stars in $\sim 80\%$ of the LV dwarf galaxies studied -- both satellites and isolated systems -- are dispersion-supported. In particular, we show that $6/10$ *isolated* dwarfs in our sample have $v_{\rm rot} / \sigma < 1.0$. All have $v_{\rm rot} / \sigma \lesssim 2.0$. These results challenge the traditional view that the stars in gas-rich dwarf irregulars (dIrrs) are distributed in cold, rotationally-supported stellar disks, while gas-poor dwarf spheroidals (dSphs) are kinematically distinct in having dispersion-supported stars. We see no clear trend between $v_{\rm rot} / \sigma$ and distance to the closest $\rm L_{\star}$ galaxy, nor between $v_{\rm rot} / \sigma$ and $M_{\star}$ within our mass range. We apply the same Bayesian analysis to four FIRE hydrodynamic zoom-in simulations of isolated dwarf galaxies ($10^9 M_{\odot} < M_{\rm vir} < 10^{10} M_{\rm \odot}$) and show that the simulated *isolated* dIrr galaxies have stellar ellipticities and stellar $v_{\rm rot} / \sigma$ ratios that are consistent with the observed population of dIrrs *and* dSphs without the need to subject these dwarfs to any external perturbations or tidal forces. We posit that most dwarf galaxies form as puffy, dispersion-dominated systems, rather than cold, angular momentum-supported disks. If this is the case, then transforming a dIrr into a dSph may require little more than removing its gas.

We have used dedicated 0.7m telescopes in California and Israel to image the halos of ~200 galaxies in the Local Volume to 29 mag/sq arcsec, the sample mainly drawn from the 2MASS Large Galaxy Atlas (LGA). We supplement the LGA sample with dwarf galaxies and more distant giant ellipticals. Low surface brightness halos exceeding 50 kpc in diameter are found only in galaxies more luminous than L* and classic interaction signatures are relatively infrequent. Halo diameter is correlated with total galaxy luminosity. Extended low surface brightness halos are present even in galaxies as faint as M_V=-18. Edge-on galaxies with boxy bulges tend to lack extended spheroidal halos, while those with large classical bulges exhibit extended round halos, supporting the notions that boxy or barlike bulges originate from disks. Most face-on spiral galaxies present features that appear to be irregular extensions of spiral arms, although rare cases show smooth boundaries with no sign of star formation. Although we serendipitously discovered a dwarf galaxy undergoing tidal disruption in the halo of NGC 4449, we found no comparable examples in our general survey. A search for similar examples in the Local Volume identified hcc087, a tidally disrupting dwarf galaxy in the Hercules Cluster, but we do not confirm an anomalously large half-light radius reported for the dwarf VCC 1661.

We investigate the response of self-interacting dark matter (SIDM) halos to the growth of galaxy potentials using idealized simulations, each run in tandem with standard collisionless Cold Dark Matter (CDM). We find a greater diversity in the SIDM halo profiles compared to the CDM halo profiles. If the stellar gravitational potential strongly dominates in the central parts of a galaxy, then SIDM halos can be as dense as CDM halos on observable scales. For extreme cases with highly compact disks core collapse can occur, leading to SIDM halos that are denser and cuspier than their CDM counterparts. If the stellar potential is not dominant, then SIDM halos retain constant density cores with densities far below CDM predictions. When a disk potential is present, the inner SIDM halo becomes \em{more flattened} in the disk plane than the CDM halo. These results are in excellent quantitative agreement with the predictions of Kaplinghat et al. (2014). We also simulated a galaxy cluster halo with a central stellar distribution similar to the brightest central galaxy of the cluster A2667. A SIDM halo simulated with cross section over mass $\sigma/m = 0.1\ \mathrm{cm^2 g^{-1}}$ provides a good match to the measured dark matter density profile of A2667, while an adiabatically-contracted CDM halo is denser and cuspier. The cored profile of the same halo simulated with $\sigma/m = 0.5\ \mathrm{cm^2 g^{-1}}$ is not dense enough to match A2667. Our findings are in agreement with previous results that $\sigma/m \gtrsim 0.1\ \mathrm{cm^2 g^{-1}}$ is disfavored for dark matter collision velocities in excess of about 1500 km/s. More generally, the predictive cross-talk between baryonic potentials and SIDM density distributions offers new directions for constraining SIDM cross sections in massive galaxies where baryons are dynamically important.

Recent studies of galaxies in the local Universe, including those in the Local Group, find that the efficiency of environmental (or satellite) quenching increases dramatically at satellite stellar masses below ~ $10^8\ {\rm M}_{\odot}$. This suggests a physical scale where quenching transitions from a slow "starvation" mode to a rapid "stripping" mode at low masses. We investigate the plausibility of this scenario using observed HI surface density profiles for a sample of 66 nearby galaxies as inputs to analytic calculations of ram-pressure and viscous stripping. Across a broad range of host properties, we find that stripping becomes increasingly effective at $M_{*} < 10^{8-9}\ {\rm M}_{\odot}$, reproducing the critical mass scale observed. However, for canonical values of the circumgalactic medium density ($n_{\rm halo} < 10^{-3.5}$ ${\rm cm}^{-3}$), we find that stripping is not fully effective; infalling satellites are, on average, stripped of < 40 - 70% of their cold gas reservoir, which is insufficient to match observations. By including a host halo gas distribution that is clumpy and therefore contains regions of higher density, we are able to reproduce the observed HI gas fractions (and thus the high quenched fraction and short quenching timescale) of Local Group satellites, suggesting that a host halo with clumpy gas may be crucial for quenching low-mass systems in Local Group-like (and more massive) host halos.

Near-field cosmology -- using detailed observations of the Local Group and its environs to study wide-ranging questions in galaxy formation and dark matter physics -- has become a mature and rich field over the past decade. There are lingering concerns, however, that the relatively small size of the present-day Local Group ($\sim 2$ Mpc diameter) imposes insurmountable sample-variance uncertainties, limiting its broader utility. We consider the region spanned by the Local Group's progenitors at earlier times and show that it reaches $3' \approx 7$ co-moving Mpc in linear size (a volume of $\approx 350\,{\rm Mpc}^3$) at $z=7$. This size at early cosmic epochs is large enough to be representative in terms of the matter density and counts of dark matter halos with $M_{\rm vir}(z=7) \lesssim 2\times 10^{9}\,M_{\odot}$. The Local Group's stellar fossil record traces the cosmic evolution of galaxies with $10^{3} \lesssim M_{\star}(z=0) / M_{\odot} \lesssim 10^{9}$ (reaching $M_{1500} > -9$ at $z\sim7$) over a region that is comparable to or larger than the Hubble Ultra-Deep Field (HUDF) for the entire history of the Universe. It is highly complementary to the HUDF, as it probes much fainter galaxies but does not contain the intrinsically rarer, brighter sources that are detectable in the HUDF. Archaeological studies in the Local Group also provide the ability to trace the evolution of individual galaxies across time as opposed to evaluating statistical connections between temporally distinct populations. In the JWST era, resolved stellar populations will probe regions larger than the HUDF and any deep JWST fields, further enhancing the value of near-field cosmology.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

The Frontier Fields are a director's discretionary time campaign with HST and the Spitzer Space Telescope to see deeper into the universe than ever before. The Frontier Fields combine the power of HST and Spitzer with the natural gravitational telescopes of massive high-magnification clusters of galaxies to produce the deepest observations of clusters and their lensed galaxies ever obtained. Six clusters - Abell 2744, MACSJ0416.1-2403, MACSJ0717.5+3745, MACSJ1149.5+2223, Abell S1063, and Abell 370 - were selected based on their lensing strength, sky darkness, Galactic extinction, parallel field suitability, accessibility to ground-based facilities, HST, Spitzer and JWST observability, and pre-existing ancillary data. These clusters have been targeted by the HST ACS/WFC and WFC3/IR with coordinated parallels of adjacent blank fields for over 840 HST orbits. The Spitzer Space Telescope has dedicated > 1000 hours of director's discretionary time to obtain IRAC 3.6 and 4.5 micron imaging to ~26.5, 26.0 ABmag 5-sigma point-source depths in the six cluster and six parallel Frontier Fields. The Frontier Field parallel fields are the second-deepest observations thus far by HST with ~29th ABmag 5-sigma point source depths in seven optical - near-infrared bandpasses. Galaxies behind the Frontier Field cluster lenses experience typical magnification factors of a few, with small regions near the critical curves magnified by factors 10-100. Therefore, the Frontier Field cluster HST images achieve intrinsic depths of ~30-33 magnitudes over very small volumes. Early studies of the Frontier Fields have probed galaxies fainter than any seen before during the epoch of reionization 6 < z < 10, mapped out the cluster dark matter to unprecedented resolution, and followed lensed transient events.

We present a combination of new and archival neutral hydrogen (HI) observations and new ionized gas spectroscopic observations for sixteen galaxies in the statistically representative EDGES kinematic sample. HI rotation curves are derived from new and archival radio synthesis observations from the Very Large Array (VLA) as well as processed data products from the Westerbork Radio Synthesis Telescope (WSRT). The HI rotation curves are supplemented with optical spectroscopic integral field unit (IFU) observations using SparsePak on the WIYN 3.5 m telescope to constrain the central ionized gas kinematics in twelve galaxies. The full rotation curves of each galaxy are decomposed into baryonic and dark matter halo components using 3.6$\mu$m images from the Spitzer Space Telescope for the stellar content, the neutral hydrogen data for the atomic gas component, and, when available, CO data from the literature for the molecular gas component. Differences in the inferred distribution of mass are illustrated under fixed stellar mass-to-light ratio (M/L) and maximum disc/bulge assumptions in the rotation curve decomposition.

Sterile neutrinos comprise an entire class of dark matter models that, depending on their production mechanism, can be hot, warm, or cold dark matter. We simulate the Local Group and representative volumes of the Universe in a variety of sterile neutrino models, all of which are consistent with the possible existence of a radiative decay line at ~3.5 keV. We compare models of production via resonances in the presence of a lepton asymmetry (suggested by Shi & Fuller 1999) to "thermal" models. We find that properties in the highly nonlinear regime - e.g., counts of satellites and internal properties of halos and subhalos - are insensitive to the precise fall-off in power with wavenumber, indicating that nonlinear evolution essentially washes away differences in the initial (linear) matter power spectrum. In the quasi-linear regime at higher redshifts, however, quantitative differences in the 3D matter power spectra remain, raising the possibility that such models can be tested with future observations of the Lyman-alpha forest. While many of the sterile neutrino models largely eliminate multiple small-scale issues within the Cold Dark Matter (CDM) paradigm, we show that these models may be ruled out in the near future via discoveries of additional dwarf satellites in the Local Group.

The anomalous 3.55 keV X-ray line recently detected towards a number of massive dark matter objects may be interpreted as the radiative decays of 7.1 keV mass sterile neutrino dark matter. Depending on its parameters, the sterile neutrino can range from cold to warm dark matter with small-scale suppression that differs in form from commonly-adopted thermal warm dark matter. Here, we numerically investigate the subhalo properties for 7.1 keV sterile neutrino dark matter produced via the resonant Shi-Fuller mechanism. Using accurate matter power spectra, we run cosmological zoom-in simulations of a Milky Way-sized halo and explore the abundance of massive subhalos, their radial distributions, and their internal structure. We also simulate the halo with thermal 2.0 keV warm dark matter for comparison and discuss quantitative differences. We find that the resonantly produced sterile neutrino model for the 3.55 keV line provides a good description of structures in the Local Group, including the number of satellite dwarf galaxies and their radial distribution, and largely mitigates the too-big-to-fail problem. Future searches for satellite galaxies by deep surveys, such as the Dark Energy Survey, Large Synoptic Survey Telescope, and Wide Field Infrared Survey Telescope, will be a strong direct test of warm dark matter scenarios.

We present ultra-high resolution cosmological hydrodynamic simulations of $M_*\simeq10^{4-6.3}M_{\odot}$ dwarf galaxies that form within $M_{v}=10^{9.5-10}M_{\odot}$ dark matter halos. Our simulations rely on the FIRE implementation of star formation feedback and were run with high enough force and mass resolution to directly resolve stellar and dark matter structure on the ~200 pc scales of interest for classical and ultra-faint dwarfs in the Local Group. The resultant galaxies sit on the $M_*$ vs. $M_{v}$ relation required to match the Local Group stellar mass function. They have bursty star formation histories and also form with half-light radii and metallicities that broadly match those observed for local dwarfs at the same stellar mass. We demonstrate that it is possible to create a large (~1 kpc) dark matter core in a cosmological simulation of an $M_*\simeq10^{6.5}M_{\odot}$ dwarf galaxy that resides within an $M_{v}=10^{10}M_{\odot}$ halo -- precisely the scale of interest for resolving the Too Big to Fail problem. However, these large cores are not ubiquitous and appear to correlate closely with the star formation histories of the dwarfs: dark matter cores are largest in systems that form their stars late ($z\lesssim2$), after the early epoch of cusp building mergers has ended. Our $M_*\simeq10^4M_{\odot}$ dwarf retains a cuspy dark matter halo density profile that matches almost identically that of a dark-matter only run of the same system. Despite forming in a field environment, this very low mass dwarf has observable properties that match closely to those of ultra-faint satellite galaxies of the Milky Way, including a uniformly old stellar population (>10 Gyr). Though ancient, most of the stars in our ultra-faint form after reionization; the UV field acts mainly to suppress fresh gas accretion, not to boil away gas that is already present in the proto-dwarf.

We constrain the evolution of the galaxy stellar mass function from 2 < z < 5 for galaxies with stellar masses as low as 10^5 Msun by combining star formation histories of Milky Way satellite galaxies derived from deep Hubble Space Telescope observations with merger trees from the ELVIS suite of N-body simulations. This approach extends our understanding more than two orders of magnitude lower in stellar mass than is currently possible by direct imaging. We find the faint end slopes of the mass functions to be alpha= -1.42(+0.07/-0.05) at z = 2 and alpha = -1.57^(+0.06/-0.06) at z = 5, and show the slope only weakly evolves from z = 5 to z = 0. Our findings are in stark contrast to a number of direct detection studies that suggest slopes as steep as alpha = -1.9 at these epochs. Such a steep slope would result in an order of magnitude too many luminous Milky Way satellites in a mass regime that is observationally complete (Mstar > 2*10^5 Msun at z = 0). The most recent studies from ZFOURGE and CANDELS also suggest flatter faint end slopes that are consistent with our results, but with a lower degree of precision. This work illustrates the strong connections between low and high-z observations when viewed through the lens of LCDM numerical simulations.

We infer the UV luminosities of Local Group galaxies at early cosmic times ($z \sim 2$ and $z \sim 7$) by combining stellar population synthesis modeling with star formation histories derived from deep color-magnitude diagrams constructed from Hubble Space Telescope (HST) observations. Our analysis provides a basis for understanding high-$z$ galaxies - including those that may be unobservable even with the James Webb Space Telescope (JWST) - in the context of familiar, well-studied objects in the very low-$z$ Universe. We find that, at the epoch of reionization, all Local Group dwarfs were less luminous than the faintest galaxies detectable in deep HST observations of blank fields. We predict that JWST will observe $z \sim 7$ progenitors of galaxies similar to the Large Magellanic Cloud today; however, the HST Frontier Fields initiative may already be observing such galaxies, highlighting the power of gravitational lensing. Consensus reionization models require an extrapolation of the observed blank-field luminosity function at $z \approx 7$ by at least two orders of magnitude in order to maintain reionization. This scenario requires the progenitors of the Fornax and Sagittarius dwarf spheroidal galaxies to be contributors to the ionizing background at $z \sim 7$. Combined with numerical simulations, our results argue for a break in the UV luminosity function from a faint-end slope of $\alpha \sim -2$ at $M_{\rm UV} < -13$ to $\alpha \sim -1.2$ at lower luminosities. Applied to photometric samples at lower redshifts, our analysis suggests that HST observations in lensing fields at $z \sim 2$ are capable of probing galaxies with luminosities comparable to the expected progenitor of Fornax.

The vast majority of dwarf satellites orbiting the Milky Way and M31 are quenched, while comparable galaxies in the field are gas-rich and star-forming. Assuming that this dichotomy is driven by environmental quenching, we use the ELVIS suite of N-body simulations to constrain the characteristic timescale upon which satellites must quench following infall into the virial volumes of their hosts. The high satellite quenched fraction observed in the Local Group demands an extremely short quenching timescale (~ 2 Gyr) for dwarf satellites in the mass range Mstar ~ 10^6-10^8 Msun. This quenching timescale is significantly shorter than that required to explain the quenched fraction of more massive satellites (~ 8 Gyr), both in the Local Group and in more massive host halos, suggesting a dramatic change in the dominant satellite quenching mechanism at Mstar < 10^8 Msun. Combining our work with the results of complementary analyses in the literature, we conclude that the suppression of star formation in massive satellites (Mstar ~ 10^8 - 10^11 Msun) is broadly consistent with being driven by starvation, such that the satellite quenching timescale corresponds to the cold gas depletion time. Below a critical stellar mass scale of ~ 10^8 Msun, however, the required quenching times are much shorter than the expected cold gas depletion times. Instead, quenching must act on a timescale comparable to the dynamical time of the host halo. We posit that ram-pressure stripping can naturally explain this behavior, with the critical mass (of Mstar ~ 10^8 Msun) corresponding to halos with gravitational restoring forces that are too weak to overcome the drag force encountered when moving through an extended, hot circumgalactic medium.

We compare the dynamics of satellite galaxies in the Sloan Digital Sky Survey to simple models in order to test the hypothesis that a large fraction of satellites co-rotate in coherent planes. We confirm the previously-reported excess of co-rotating satellite pairs located near diametric opposition with respect to the host, but show that this signal is unlikely to be due to rotating discs (or planes) of satellites. In particular, no overabundance of co-rotating satellites pairs is observed within $\sim 20^{\circ}-50^{\circ}$ of direct opposition, as would be expected for planar distributions inclined relative to the line-of-sight. Instead, the excess co-rotation for satellite pairs within $\sim 10^{\circ}$ of opposition is consistent with random noise associated with undersampling of an underlying isotropic velocity distribution. We conclude that at most $10\%$ of the hosts in our sample harbor co-rotating satellite planes (as traced by the luminous satellite population).

The TMT Detailed Science Case describes the transformational science that the Thirty Meter Telescope will enable. Planned to begin science operations in 2024, TMT will open up opportunities for revolutionary discoveries in essentially every field of astronomy, astrophysics and cosmology, seeing much fainter objects much more clearly than existing telescopes. Per this capability, TMT's science agenda fills all of space and time, from nearby comets and asteroids, to exoplanets, to the most distant galaxies, and all the way back to the very first sources of light in the Universe. More than 150 astronomers from within the TMT partnership and beyond offered input in compiling the new 2015 Detailed Science Case. The contributing astronomers represent the entire TMT partnership, including the California Institute of Technology (Caltech), the Indian Institute of Astrophysics (IIA), the National Astronomical Observatories of the Chinese Academy of Sciences (NAOC), the National Astronomical Observatory of Japan (NAOJ), the University of California, the Association of Canadian Universities for Research in Astronomy (ACURA) and US associate partner, the Association of Universities for Research in Astronomy (AURA).

We present FIRE/Gizmo hydrodynamic zoom-in simulations of isolated dark matter halos, two each at the mass of classical dwarf galaxies ($M_{\rm vir} \simeq 10^{10} M_{\odot}$) and ultra-faint galaxies ($M_{\rm vir} \simeq 10^9 M_{\odot}$), and with two feedback implementations. The resultant central galaxies lie on an extrapolated abundance matching relation from $M_{\star} \simeq 10^6$ to $10^4 M_{\odot}$ without a break. Every host is filled with subhalos, many of which form stars. Our dwarfs with $M_{\star} \simeq 10^6 M_{\odot}$ each have 1-2 well-resolved satellites with $M_{\star} = 3-200 \times 10^3 M_{\odot}$. Even our isolated ultra-faint galaxies have star-forming subhalos. If this is representative, dwarf galaxies throughout the universe should commonly host tiny satellite galaxies of their own. We combine our results with the ELVIS simulations to show that targeting $\sim 50~ \rm kpc$ regions around nearby isolated dwarfs could increase the chances of discovering ultra-faint galaxies by $\sim 35\%$ compared to random halo pointings, and specifically identify the region around the Phoenix dwarf galaxy as a good potential target. The well-resolved ultra-faint galaxies in our simulations ($M_{\star} \simeq 3 - 30 \times 10^3 M_{\odot}$) form within $M_{\rm peak} \simeq 0.5 - 3 \times 10^9 M_{\odot}$ halos. Each has a uniformly ancient stellar population ($ > 10~ \rm Gyr$) owing to reionization-related quenching. More massive systems, in contrast, all have late-time star formation. Our results suggest that $M_{\rm halo} \simeq 5 \times 10^9 M_{\odot}$ is a probable dividing line between halos hosting reionization "fossils" and those hosting dwarfs that can continue to form stars in isolation after reionization.

We propose that luminous transients, including novae and supernovae, can be used to detect the faintest galaxies in the universe. Beyond a few Mpc, dwarf galaxies with stellar masses $<10^6 M_{\odot}$ will likely be too faint and/or too low in surface brightness to be directly detected in upcoming large area ground-based photometric surveys. However, single epoch LSST photometry will be able to detect novae to distances of $\sim30$ Mpc and SNe to Gpc-scale distances. Depending on the form of the stellar mass-halo mass relation and the underlying star formation histories of low mass dwarfs, the expected nova rates will be a few to $\sim100$ yr$^{-1}$ and the expected SN rates (including both type Ia and core-collapse) will be $\sim10^2-10^4$ within the observable ($4\pi$ sr) volume. The transient rate associated with intrahalo stars will be comparably large, but these transients will be located close to bright galaxies, in contrast to the dwarfs, which should trace the underlying large scale structure of the cosmic web. Aggressive follow-up of hostless transients has the potential to uncover the predicted enormous population of low mass field dwarf galaxies.

In this paper we present the MOSFIRE Deep Evolution Field (MOSDEF) survey. The MOSDEF survey aims to obtain moderate-resolution (R=3000-3650) rest-frame optical spectra (~3700-7000 Angstrom) for ~1500 galaxies at 1.37<z<3.80 in three well-studied CANDELS fields: AEGIS, COSMOS, and GOODS-N. Targets are selected in three redshift intervals: 1.37<z<1.70, 2.09<z<2.61, and 2.95<z<3.80, down to fixed H_AB (F160W) magnitudes of 24.0, 24.5 and 25.0, respectively, using the photometric and spectroscopic catalogs from the 3D-HST survey. We target both strong nebular emission lines (e.g., [OII], Hbeta, [OIII], 5008, Halpha, [NII], and [SII]) and stellar continuum and absorption features (e.g., Balmer lines, Ca-II H and K, Mgb, 4000 Angstrom break). Here we present an overview of our survey, the observational strategy, the data reduction and analysis, and the sample characteristics based on spectra obtained during the first 24 nights. To date, we have completed 21 masks, obtaining spectra for 591 galaxies. For ~80% of the targets we derive a robust redshift from either emission or absorption lines. In addition, we confirm 55 additional galaxies, which were serendipitously detected. The MOSDEF galaxy sample includes unobscured star-forming, dusty star-forming, and quiescent galaxies and spans a wide range in stellar mass (~10^9-10^11.5 Msol) and star formation rate (~10^0-10^3 Msol/yr). The spectroscopically confirmed sample is roughly representative of an H-band limited galaxy sample at these redshifts. With its large sample size, broad diversity in galaxy properties, and wealth of available ancillary data, MOSDEF will transform our understanding of the stellar, gaseous, metal, dust, and black hole content of galaxies during the time when the universe was most active.

We present results from multiwavelength observations of the galaxy NGC5005. We use new neutral hydrogen (HI) observations from the Very Large Array to examine the neutral gas morphology and kinematics. We find an HI disk with a well-behaved flat rotation curve in the radial range 20\arcsec-140\arcsec. Ionized gas observations from the SparsePak integral field unit on the WIYN 3.5m telescope provide kinematics for the central 70\arcsec. We use both the SparsePak and HI velocity fields to derive a rotation curve for NGC5005. Deep 3.6{\mu}m observations from the Spitzer Space Telescope probe the faint extended stellar population of NGC5005. The images reveal a large stellar disk with a high surface brightness component that transitions to a low surface brightness component at a radius nearly 1.6 times farther than the extent of the gas disk detected in HI. The 3.6{\mu}m image is also decomposed into bulge and disk components to account for the stellar light distribution. Optical broadband B and R and narrowband H{\alpha} from the WIYN 0.9m telescope complement the 3.6{\mu}m data by providing information about the dominant stellar population and current star formation activity. The neutral and ionized gas rotation curve is used along with the stellar bulge and disk light profiles to decompose the mass distributions in NGC5005 and determine a dark matter halo model. The maximum stellar disk contribution to the total rotation curve is only about 70\%, suggesting that dark matter makes a significant contribution to the dynamics at all radii.

Highly ionized, z=0 metal absorption lines detected in the X-ray spectra of background active galactic nuclei (AGNs) provide an effective method to probe the hot ($T\sim10^6$ K) gas and its metal content in and around the Milky Way. We present an all-sky survey of the $K_{\alpha}$ transition of the local O VII absorption lines obtained by Voigt-profile fitting archival XMM-Newton observations. A total of 43 AGNs were selected, among which 12 are BL Lac-type AGNs, and the rest are Seyfert 1 galaxies. At above the $3\sigma$ level the local O VII absorption lines were detected in 21 AGNs, among which 7 were newly discovered in this work. The sky covering fraction, defined as the ratio between the number of detections and the sample size, increases from at about 40% for all targets to 100% for the brightest targets, suggesting a uniform distribution of the O VII absorbers. We correlate the line equivalent width with the Galactic coordinates and do not find any strong correlations between these quantities. Some AGNs have warm absorbers that may complicate the analysis of the local X-ray absorber since the recession velocity can be compensated by the outflow velocity, especially for the nearby targets. We discuss the potential impact of the warm absorbers on our analysis. A comprehensive theoretical modelling of the X-ray absorbers will be presented in a later paper.

We investigate the effect of self-interacting dark matter (SIDM) on the density profiles of $V_{\rm max} \simeq 40~km~s^{-1}$ isolated dwarf dark matter halos -- the scale of relevance for the too big to fail problem (TBTF) -- using very high-resolution cosmological zoom simulations. Each halo has millions of particles within its virial radius. We find that SIDM models with cross sections per unit mass spanning the range \sigma/m = $0.5 - 50$ $cm^2~g^{-1}$ alleviate TBTF and produce constant density cores of size 300-1000 pc, comparable to the half-light radii of $M_\star$ ~ $10^{5-7}$ $M_\odot$ dwarfs. The largest, lowest density cores develop for cross sections in the middle of this range, \sigma/m ~ $5-10~cm^2~g^{-1}$. Our largest SIDM cross section run (\sigma/m = $50~cm^2~g^{-1}$) develops a slightly denser core owing to mild core-collapse behavior, but it remains less dense than the CDM case and retains a constant density core profile. Our work suggests that SIDM cross sections as large or larger than $50~cm^2~g^{-1}$ remain viable on velocity scales of dwarf galaxies ($v_{\rm rms}$ ~ $40~km~s^{-1}$). The range of SIDM cross sections that alleviate TBTF and the cusp/core problem spans at least two orders of magnitude and therefore need not be particularly fine-tuned.

The BICEP2 results, when interpreted as a gravitational wave signal and combined with other CMB data, suggest a roll-off in power towards small scales in the primordial matter power spectrum. Among the simplest possibilities is a running of the spectral index. Here we show that the preferred level of running alleviates small-scale issues within the $\Lambda$CDM model, more so even than viable WDM models. We use cosmological zoom-in simulations of a Milky Way-size halo along with full-box simulations to compare predictions among four separate cosmologies: a BICEP2-inspired running index model ($\alpha_s$ = -0.024), two fixed-tilt $\Lambda$CDM models motivated by Planck, and a 2.6 keV thermal WDM model. We find that the running BICEP2 model reduces the central densities of large dwarf-size halos ($V_\mathrm{max}$ ~ 30 - 80 km s$^{-1}$) and alleviates the too-big-to-fail problem significantly compared to our adopted Planck and WDM cases. Further, the BICEP2 model suppresses the count of small subhalos by ~50% relative to Planck models, and yields a significantly lower "boost" factor for dark matter annihilation signals. Our findings highlight the need to understand the shape of the primordial power spectrum in order to correctly interpret small-scale data.

We present a series of high-resolution cosmological simulations of galaxy formation to z=0, spanning halo masses ~10^8-10^13 M_sun, and stellar masses ~10^4-10^11. Our simulations include fully explicit treatment of both the multi-phase ISM (molecular through hot) and stellar feedback. The stellar feedback inputs (energy, momentum, mass, and metal fluxes) are taken directly from stellar population models. These sources of stellar feedback, with zero adjusted parameters, reproduce the observed relation between stellar and halo mass up to M_halo~10^12 M_sun (including dwarfs, satellites, MW-mass disks, and small groups). By extension, this leads to reasonable agreement with the stellar mass function for M_star<10^11 M_sun. We predict weak redshift evolution in the M_star-M_halo relation, consistent with current constraints to z>6. We find that the M_star-M_halo relation is insensitive to numerical details, but is sensitive to the feedback physics. Simulations with only supernova feedback fail to reproduce the observed stellar masses, particularly in dwarf and high-redshift galaxies: radiative feedback (photo-heating and radiation pressure) is necessary to disrupt GMCs and enable efficient coupling of later supernovae to the gas. Star formation rates agree well with the observed Kennicutt relation at all redshifts. The galaxy-averaged Kennicutt relation is very different from the numerically imposed law for converting gas into stars in the simulation, and is instead determined by self-regulation via stellar feedback. Feedback reduces star formation rates considerably and produces a reservoir of gas that leads to rising late-time star formation histories significantly different from the halo accretion history. Feedback also produces large short-timescale variability in galactic SFRs, especially in dwarfs. Many of these properties are not captured by common 'sub-grid' galactic wind models.

We present the metallicity distribution of red giant branch (RGB) stars in M31's stellar halo, derived from photometric metallicity estimates for over 1500 spectroscopically confirmed RGB halo stars. The stellar sample comes from 38 halo fields observed with the Keck/DEIMOS spectrograph, ranging from 9 to 175 kpc in projected distance from M31's center, and includes 52 confirmed M31 halo stars beyond 100 kpc. While a wide range of metallicities is seen throughout the halo, the metal-rich peak of the metallicity distribution function becomes significantly less prominent with increasing radius. The metallicity profile of M31's stellar halo shows a continuous gradient from 9 to ~100 kpc, with a magnitude of -0.01 dex/kpc. The stellar velocity distributions in each field are used to identify stars that are likely associated with tidal debris features. The removal of tidal debris features does not significantly alter the metallicity gradient in M31's halo: a gradient is maintained in fields spanning 10 to 90 kpc. We analyze the halo metallicity profile, as well as the relative metallicities of stars associated with tidal debris features and the underlying halo population, in the context of current simulations of stellar halo formation. We argue that the large scale gradient in M31's halo implies M31 accreted at least one relatively massive progenitor in the past, while the field to field variation seen in the metallicity profile indicates that multiple smaller progenitors are likely to have contributed substantially to M31's outer halo.

Using the Sloan Digital Sky Survey, we examine the quenching of satellite galaxies around isolated Milky Way-like hosts in the local Universe. We find that the efficiency of satellite quenching around isolated galaxies is low and roughly constant over two orders of magnitude in satellite stellar mass ($M_{*}$ = $10^{8.5}-10^{10.5} \, M_{\odot}$), with only $\sim~20\%$ of systems quenched as a result of environmental processes. While largely independent of satellite stellar mass, satellite quenching does exhibit clear dependence on the properties of the host. We show that satellites of passive hosts are substantially more likely to be quenched than those of star-forming hosts, and we present evidence that more massive halos quench their satellites more efficiently. These results extend trends seen previously in more massive host halos and for higher satellite masses. Taken together, it appears that galaxies with stellar masses larger than about $10^{8}~M_{\odot}$ are uniformly resistant to environmental quenching, with the relative harshness of the host environment likely serving as the primary driver of satellite quenching. At lower stellar masses ($< 10^{8}~M_{\odot}$), however, observations of the Local Group suggest that the vast majority of satellite galaxies are quenched, potentially pointing towards a characteristic satellite mass scale below which quenching efficiency increases dramatically.

Reionizing the Universe with galaxies appears to require significant star formation in low-mass halos at early times, while local dwarf galaxy counts tell us that star formation has been minimal in small halos around us today. Using simple models and the ELVIS simulation suite, we show that reionization scenarios requiring appreciable star formation in halos with $M_{\rm vir} \approx 10^{8}\,M_{\odot}$ at $z=8$ are in serious tension with galaxy counts in the Local Group. This tension originates from the seemingly inescapable conclusion that 30 - 60 halos with $M_{\rm vir} > 10^{8}\,M_{\odot}$ at $z=8$ will survive to be distinct bound satellites of the Milky Way at $z = 0$. Reionization models requiring star formation in such halos will produce dozens of bound galaxies in the Milky Way's virial volume today (and 100 - 200 throughout the Local Group), each with $\gtrsim 10^{5}\,M_{\odot}$ of old stars ($\gtrsim 13$ Gyr). This exceeds the stellar mass function of classical Milky Way satellites today, even without allowing for the (significant) post-reionization star formation observed in these galaxies. One possible implication of these findings is that star formation became sharply inefficient in halos smaller than $\sim 10^9 \,M_{\odot}$ at early times, implying that the high-$z$ luminosity function must break at magnitudes brighter than is often assumed (at ${\rm M_{UV}} \approx -14$). Our results suggest that JWST (and possibly even HST with the Frontier Fields) may realistically detect the faintest galaxies that drive reionization. It remains to be seen how these results can be reconciled with the most sophisticated simulations of early galaxy formation at present, which predict substantial star formation in $M_{\rm vir} \sim 10^8 \, M_{\odot}$ halos during the epoch of reionization.

We use $N$-body simulations to show that high-redshift galaxy counts provide an interesting constraint on the nature of dark matter, specifically Warm Dark Matter (WDM), owing to the lack of early structure formation these models. Our simulations include three WDM models with thermal-production masses of 0.8 keV, 1.3 keV, and 2.6 keV, as well as CDM. Assuming a relationship between dark halo mass and galaxy luminosity that is set by the observed luminosity function at bright magnitudes, we find that 0.8 keV WDM is disfavored by direct galaxy counts in the Hubble Ultra Deep Field at $>\!\!10\sigma$. Similarly, 1.3 keV WDM is statistically inconsistent at $2.2\sigma$. Future observations with JWST (and possibly HST via the Frontier Fields) could rule out $1.3$ keV WDM at high significance, and may be sensitive to WDM masses greater than 2.6 keV. We also examine the ability of galaxies in these WDM models to reionize the universe, and find that 0.8 keV and 1.3 keV WDM produce optical depths to the Cosmic Microwave Background (CMB) that are inconsistent at 68% C.L. with current Planck results, even with extremely high ionizing radiation escape fractions, and 2.6 keV WDM requires an optimistic escape fraction to yield an optical depth consistent with Planck data. Although CMB optical depth calculations are model dependent, we find a strong challenge for stellar processes alone to reionize the universe in a 0.8 keV and 1.3 keV WDM cosmology.

We study dwarf satellite galaxy quenching using observations from the Geha et al. (2012) NSA/SDSS catalog together with LCDM cosmological simulations to facilitate selection and interpretation. We show that fewer than 30% of dwarfs (M* ~ 10^8.5-10^9.5 Msun) identified as satellites within massive host halos (Mhost ~ 10^12.5-10^14 Msun) are quenched, in spite of the expectation from simulations that half of them should have been accreted more than 6 Gyr ago. We conclude that whatever the action triggering environmental quenching of dwarf satellites, the process must be highly inefficient. We investigate a series of simple, one-parameter quenching models in order to understand what is required to explain the low quenched fraction and conclude that either the quenching timescale is very long (> 9.5 Gyr, a "slow starvation" scenario) or that the environmental trigger is not well matched to accretion within the virial volume. We discuss these results in light of the fact that most of the low mass dwarf satellites in the Local Group are quenched, a seeming contradiction that could point to a characteristic mass scale for satellite quenching.

We compare the dynamical masses of dwarf galaxies in the Local Group (LG) to the predicted masses of halos in the ELVIS suite of $\Lambda$CDM simulations, a sample of 48 Galaxy-size hosts, 24 of which are in paired configuration similar to the LG. We enumerate unaccounted-for dense halos ($V_\mathrm{max} \gtrsim 25$ km s$^{-1}$) in these volumes that at some point in their histories were massive enough to have formed stars in the presence of an ionizing background ($V_\mathrm{peak} > 30$ km s$^{-1}$). Within 300 kpc of the Milky Way, the number of unaccounted-for massive halos ranges from 2 - 25 over our full sample. Moreover, this "too big to fail" count grows as we extend our comparison to the outer regions of the Local Group: within 1.2 Mpc of either giant we find that there are 12-40 unaccounted-for massive halos. This count excludes volumes within 300 kpc of both the MW and M31, and thus should be largely unaffected by any baryonically-induced environmental processes. According to abundance matching -- specifically abundance matching that reproduces the Local Group stellar mass function -- all of these missing massive systems should have been quite bright, with $M_\star > 10^6M_\odot$. Finally, we use the predicted density structure of outer LG dark matter halos together with observed dwarf galaxy masses to derive an $M_\star-V_\mathrm{max}$ relation for LG galaxies that are outside the virial regions of either giant. We find that there is no obvious trend in the relation over three orders of magnitude in stellar mass (a "common mass" relation), from $M_\star \sim 10^8 - 10^5 M_\odot$, in drastic conflict with the tight relation expected for halos that are unaffected by reionization. Solutions to the too big to fail problem that rely on ram pressure stripping, tidal effects, or statistical flukes appear less likely in the face of these results.

We compare the kinematics of M31's satellite galaxies to the mass profiles of the subhaloes they are expected to inhabit in LCDM. We consider the most massive subhaloes of an approximately M31-sized halo, following the assumption of a monotonic galaxy luminosity-to-subhalo mass mapping. While this abundance matching relation is consistent with the kinematic data for galaxies down to the luminosity of the bright satellites of the Milky Way and M31, it is not consistent with kinematic data for fainter dwarf galaxies (those with L <~ 10^8 Lsun). Comparing the kinematics of M31's dwarf Spheroidal (dSph) satellites to the subhaloes reveals that M31's dSph satellites are too low density to be consistent with the subhaloes' mass profiles. A similar discrepancy has been reported between Milky Way dSphs and their predicted subhaloes, the "too big to fail" problem. By contrast, total mass profiles of the dwarf Elliptical (and similarly bright) satellites are consistent with the subhaloes. However, they suffer from large systematic uncertainties in their dark matter content because of substantial (and potentially dominant) contributions from baryons within their half-light radii.

We measured velocities of 862 individual red giant stars in seven isolated dwarf galaxies in the Local Group: NGC 6822, IC 1613, VV 124 (UGC 4879), the Pegasus dwarf irregular galaxy (DDO 216), Leo A, Cetus, and Aquarius (DDO 210). We also computed velocity dispersions, taking into account the measurement uncertainties on individual stars. None of the isolated galaxies is denser than the densest Local Group satellite galaxy. Furthermore, the isolated dwarf galaxies have no obvious distinction in the velocity dispersion--half-light radius plane from the satellite galaxies of the Milky Way and M31. The similarity of the isolated and satellite galaxies' dynamics and structural parameters imposes limitations on environmental solutions to the too-big-to-fail problem, wherein there are fewer dense dwarf satellite galaxies than would be expected from cold dark matter simulations. This data set also has many other applications for dwarf galaxy evolution, including the transformation of dwarf irregular into dwarf spheroidal galaxies. We intend to explore these issues in future work.

We introduce a set of high-resolution dissipationless simulations that model the Local Group (LG) in a cosmological context: Exploring the Local Volume in Simulations (ELVIS). The suite contains 48 Galaxy-size halos, each within high-resolution volumes that span 2-5 Mpc in size, and each resolving thousands of systems with masses below the atomic cooling limit. Half of the ELVIS galaxy halos are in paired configurations similar to the Milky Way (MW) and M31; the other half are isolated, mass-matched analogs. We find no difference in the abundance or kinematics of substructure within the virial radii of isolated versus paired hosts. On Mpc scales, however, LG-like pairs average almost twice as many companions and the velocity field is kinematically hotter and more complex. We present a refined abundance matching relation between stellar mass and halo mass that reproduces the observed satellite stellar mass functions of the MW and M31 down to the regime where incompleteness is an issue, $M_\star \sim 5\times 10^5 \, M_\odot$. Within a larger region spanning approximately 3 Mpc, the same relation predicts that there should be $\sim$ 1000 galaxies with $M_\star > 10^{3}\,M_\odot$ awaiting discovery. We show that up to 50% of halos within 1 Mpc of the MW or M31 could be systems that have previously been within the virial radius of either giant. By associating never-accreted halos with gas-rich dwarfs, we show that there are plausibly 50 undiscovered dwarf galaxies with HI masses $> 10^5\,M_\odot$ within the Local Volume. The radial velocity distribution of these predicted gas-rich dwarfs can be used to inform follow-up searches based on ultra-compact high-velocity clouds found in the ALFALFA survey.

We present 3D kinematic observations of stars within the central 0.5 pc of the Milky Way nuclear star cluster using adaptive optics imaging and spectroscopy from the Keck telescopes. Recent observations have shown that the cluster has a shallower surface density profile than expected for a dynamically relaxed cusp, leading to important implications for its formation and evolution. However, the true three dimensional profile of the cluster is unknown due to the difficulty in de-projecting the stellar number counts. Here, we use spherical Jeans modeling of individual proper motions and radial velocities to constrain for the first time, the de-projected spatial density profile, cluster velocity anisotropy, black hole mass ($M_\mathrm{BH}$), and distance to the Galactic center ($R_0$) simultaneously. We find that the inner stellar density profile of the late-type stars, $\rho(r)\propto r^{-\gamma}$ to have a power law slope $\gamma=0.05_{-0.60}^{+0.29}$, much more shallow than the frequently assumed Bahcall $\&$ Wolf slope of $\gamma=7/4$. The measured slope will significantly affect dynamical predictions involving the cluster, such as the dynamical friction time scale. The cluster core must be larger than 0.5 pc, which disfavors some scenarios for its origin. Our measurement of $M_\mathrm{BH}=5.76_{-1.26}^{+1.76}\times10^6$ $M_\odot$ and $R_0=8.92_{-0.55}^{+0.58}$ kpc is consistent with that derived from stellar orbits within 1$^{\prime\prime}$ of Sgr A*. When combined with the orbit of S0-2, the uncertainty on $R_0$ is reduced by 30% ($8.46_{-0.38}^{+0.42}$ kpc). We suggest that the MW NSC can be used in the future in combination with stellar orbits to significantly improve constraints on $R_0$.

We examine the star formation properties of bright (~0.1 L*) satellites around isolated ~L* hosts in the local Universe using spectroscopically confirmed systems in the Sloan Digital Sky Survey DR7. Our selection method is carefully designed with the aid of N-body simulations to avoid groups and clusters. We find that satellites are significantly more likely to be quenched than a stellar mass-matched sample of isolated galaxies. Remarkably, this quenching occurs only for satellites of hosts that are themselves quenched: while star formation is unaffected in the satellites of star-forming hosts, satellites around quiescent hosts are more than twice as likely to be quenched than stellar-mass matched field samples. One implication of this is that whatever shuts down star formation in isolated, passive L* galaxies also plays at least an indirect role in quenching star formation in their bright satellites. The previously-reported tendency for "galactic conformity" in color/morphology may be a by-product of this host-specific quenching dichotomy. The S\'ersic indices of quenched satellites are statistically identical to those of field galaxies with the same specific star formation rates, suggesting that environmental and secular quenching give rise to the same morphological structure. By studying the distribution of pairwise velocities between the hosts and satellites, we find dynamical evidence that passive host galaxies reside in dark matter halos that are ~45% more massive than those of star-forming host galaxies of the same stellar mass. We emphasize that even around passive hosts, the mere fact that galaxies become satellites does not typically result in star formation quenching: we find that only ~30% of ~0.1 L* galaxies that fall in from the field are quenched around passive hosts, compared with ~0% around star forming hosts.

We perform a suite of multimass cosmological zoom simulations of individual dark matter halos and explore how to best select Lagrangian regions for resimulation without contaminating the halo of interest with low-resolution particles. Such contamination can lead to significant errors in the gas distribution of hydrodynamical simulations, as we show. For a fixed Lagrange volume, we find that the chance of contamination increases systematically with the level of zoom. In order to avoid contamination, the Lagrangian volume selected for resimulation must increase monotonically with the resolution difference between parent box and the zoom region. We provide a simple formula for selecting Lagrangian regions (in units of the halo virial volume) as a function of the level of zoom required. We also explore the degree to which a halo's Lagrangian volume correlates with other halo properties (concentration, spin, formation time, shape, etc.) and find no significant correlation. There is a mild correlation between Lagrange volume and environment, such that halos living in the most clustered regions have larger Lagrangian volumes. Nevertheless, selecting halos to be isolated is not the best way to ensure inexpensive zoom simulations. We explain how one can safely choose halos with the smallest Lagrangian volumes, which are the least expensive to resimulate, without biasing one's sample.

We present spectroscopic metallicities of individual stars in seven gas-rich dwarf irregular galaxies (dIrrs), and we show that dIrrs obey the same mass-metallicity relation as the dwarf spheroidal (dSph) satellites of both the Milky Way and M31: Z_* ~ M_*^(0.30 +/- 0.02). The uniformity of the relation is in contradiction to previous estimates of metallicity based on photometry. This relationship is roughly continuous with the stellar mass-stellar metallicity relation for galaxies as massive as M_* = 10^12 M_sun. Although the average metallicities of dwarf galaxies depend only on stellar mass, the shapes of their metallicity distributions depend on galaxy type. The metallicity distributions of dIrrs resemble simple, leaky box chemical evolution models, whereas dSphs require an additional parameter, such as gas accretion, to explain the shapes of their metallicity distributions. Furthermore, the metallicity distributions of the more luminous dSphs have sharp, metal-rich cut-offs that are consistent with the sudden truncation of star formation due to ram pressure stripping.

We have devised a method to select galaxies that are isolated in their dark matter halo (N=1 systems) and galaxies that reside in a group of exactly two (N=2 systems). Our N=2 systems are widely-separated (up to $\sim$\,200\,$h^{-1}$\,kpc), where close galaxy-galaxy interactions are not dominant. We apply our selection criteria to two volume-limited samples of galaxies from SDSS DR6 with $M_{r}-5 \log_{10} h \leq$ -19 and -20 to study the effects of the environment of very sparse groups on galaxy colour. For satellite galaxies in a group of two, we find a red excess attributed to star formation quenching of 0.15\,$\pm$\,0.01 and 0.14\,$\pm$\,0.01 for the -19 and -20 samples, respectively, relative to isolated galaxies of the same stellar mass. Assuming N=1 systems are the progenitors of N=2 systems, an immediate-rapid star formation quenching scenario is inconsistent with these observations. A delayed-then-rapid star formation quenching scenario with a delay time of 3.3 and 3.7\,Gyr for the -19 and -20 samples, respectively, yields a red excess prediction in agreement with the observations. The observations also reveal that central galaxies in a group of two have a slight blue excess of 0.06\,$\pm$\,0.02 and 0.02\,$\pm$\,0.01 for the -19 and -20 samples, respectively, relative to N=1 populations of the same stellar mass. Our results demonstrate that even the environment of very sparse groups of luminous galaxies influence galaxy evolution and in-depth studies of these simple systems are an essential step towards understanding galaxy evolution in general.

We use high-resolution cosmological hydrodynamic simulations to study the angular momentum acquisition of gaseous halos around Milky Way sized galaxies. We find that cold mode accreted gas enters a galaxy halo with ~70% more specific angular momentum than dark matter averaged over cosmic time (though with a very large dispersion). In fact, we find that all matter has a higher spin parameter when measured at accretion than when averaged over the entire halo lifetime, and is well characterized by \lambda~0.1, at accretion. Combined with the fact that cold flow gas spends a relatively short time (1-2 dynamical times) in the halo before sinking to the center, this naturally explains why cold flow halo gas has a specific angular momentum much higher than that of the halo and often forms "cold flow disks". We demonstrate that the higher angular momentum of cold flow gas is related to the fact that it tends to be accreted along filaments.

The observed central densities of Milky Way dwarf spheroidal galaxies (dSphs) are significantly lower than the densities of the largest (Vmax about 35 km/s) subhalos found in dissipationless simulations of Galaxy-size dark matter hosts. One possible explanation is that gas removal from feedback can lower core densities enough to match observations. We model the dynamical effects of supernova feedback through the use of a time-varying central potential in high resolution, idealized numerical simulations and explore the resulting impact on the mass distributions of dwarf dark matter halos. We find that in order to match the observed central masses of M_star about 10^6 M_sun dSphs, the energy equivalent of more than 40,000 supernovae must be delivered with 100% efficiency directly to the dark matter. This energy requirement exceeds the number of supernovae that have ever exploded in most dSphs for typical initial mass functions. We also find that, per unit energy delivered and per cumulative mass removed from the galaxy, single blow-out events are more effective than repeated small bursts in reducing central dark matter densities. We conclude that it is unlikely that supernova feedback alone can solve the "Too Big to Fail" problem for Milky Way subhalos.

The cold dark matter (CDM) cosmological model has been remarkably successful in explaining cosmic structure over an enormous span of redshift, but it has faced persistent challenges from observations that probe the innermost regions of dark matter halos and the properties of the Milky Way's dwarf galaxy satellites. We review the current observational and theoretical status of these "small scale controversies." Cosmological simulations that incorporate only gravity and collisionless CDM predict halos with abundant substructure and central densities that are too high to match constraints from galaxy dynamics. The solution could lie in baryonic physics: recent numerical simulations and analytic models suggest that gravitational potential fluctuations tied to efficient supernova feedback can flatten the central cusps of halos in massive galaxies, and a combination of feedback and low star-formation efficiency could explain why most of the dark matter subhalos orbiting the Milky Way do not host visible galaxies. However, it is not clear that this solution can work in the lowest mass galaxies where discrepancies are observed. Alternatively, the small-scale conflicts could be evidence of more complex physics in the dark sector itself. For example, elastic scattering from strong dark matter self-interactions can alter predicted halo mass profiles, leading to good agreement with observations across a wide range of galaxy mass. Gravitational lensing and dynamical perturbations of tidal streams in the stellar halo provide evidence for an abundant population of low mass subhalos in accord with CDM predictions. These observational approaches will get more powerful over the next few years.

Segue 2, discovered by Belokurov et al. (2009), is a galaxy with a luminosity of only 900 L_sun. We present Keck/DEIMOS spectroscopy of 25 members of Segue 2--a threefold increase in spectroscopic sample size. The velocity dispersion is too small to be measured with our data. The upper limit with 90% (95%) confidence is sigma_v < 2.2 (2.6) km/s, the most stringent limit for any galaxy. The corresponding limit on the mass within the 3-D half-light radius (46 pc) is M_1/2 < 1.5 (2.1) x 10^5 M_sun. Segue 2 is the least massive galaxy known. We identify Segue 2 as a galaxy rather than a star cluster based the wide dispersion in [Fe/H] (from -2.85 to -1.33) among the member stars. The stars' [alpha/Fe] ratios decline with increasing [Fe/H], indicating that Segue 2 retained Type Ia supernova ejecta despite its presently small mass and that star formation lasted for at least 100 Myr. The mean metallicity, <[Fe/H]> = -2.22 +/- 0.13 (about the same as the Ursa Minor galaxy, 330 times more luminous than Segue 2), is higher than expected from the luminosity-metallicity relation defined by more luminous dwarf galaxy satellites of the Milky Way. Segue 2 may be the barest remnant of a tidally stripped, Ursa Minor-sized galaxy. If so, it is the best example of an ultra-faint dwarf galaxy that came to be ultra-faint through tidal stripping. Alternatively, Segue 2 could have been born in a very low-mass dark matter subhalo (v_max < 10 km/s), below the atomic hydrogen cooling limit.

We compare spherically-averaged radial number counts of bright (> 10^5 Lsun) dwarf satellite galaxies within 400 kpc of the Milky Way (MW) and M31 and find that the MW satellites are much more centrally concentrated. Remarkably, the two satellite systems are almost identical within the central 100 kpc, while M31 satellites outnumber MW satellites by about a factor of four at deprojected distances spanning 100 - 400 kpc. We compare the observed distributions to those predicted for LCDM suhbalos using a suite of 44 high-resolution ~10^12 halo zoom simulations, 22 of which are in pairs like the MW and M31. We find that the radial distribution of satellites around M31 is fairly typical of those predicted for subhalos, while the Milky Way's distribution is more centrally concentrated that any of our simulated LCDM halos. One possible explanation is that our census is bright (> 10^5 Lsun) MW dwarf galaxies is significantly incomplete beyond ~ 100 kpc of the Sun. If there were ~8 - 20 more bright dwarfs orbiting undetected at 100 - 400 kpc, then the Milky Way's radial distribution would fall within the range expected from subhalo distributions and alos look very much like the known M31 system. We use our simulations to demonstrate that there is enough area left unexplored by the Sloan Digital Sky Survey and its extensions that the discovery of ~10 new bright dwarfs is not implausible given the expected range of angular anisotropy of subhalos in the sky.

We present the first absolute proper motion measurement of Leo I, based on two epochs of HST ACS/WFC images separated by ~5 years. The average shift of Leo I stars with respect to ~100 background galaxies implies a proper motion of (mu_W, mu_N) = (0.1140 +/- 0.0295, -0.1256 +/- 0.0293) mas/yr. The implied Galactocentric velocity vector, corrected for the reflex motion of the Sun, has radial and tangential components V_rad = 167.9 +/- 2.8 km/s and V_tan = 101.0 +/- 34.4 km/s, respectively. We study the detailed orbital history of Leo I by solving its equations of motion backward in time for a range of plausible mass models for the Milky Way and its surrounding galaxies. Leo I entered the Milky Way virial radius 2.33 +/- 0.21 Gyr ago, most likely on its first infall. It had a pericentric approach 1.05 +/- 0.09 Gyr ago at a Galactocentric distance of 91 +/- 36 kpc. We associate these time scales with characteristic time scales in Leo I's star formation history, which shows an enhanced star formation activity ~2 Gyr ago and quenching ~1 Gyr ago. There is no indication from our calculations that other galaxies have significantly influenced Leo I's orbit, although there is a small probability that it may have interacted with either Ursa Minor or Leo II within the last ~1 Gyr. For most plausible Milky Way masses, the observed velocity implies that Leo I is bound to the Milky Way. However, it may not be appropriate to include it in models of the Milky Way satellite population that assume dynamical equilibrium, given its recent infall. Solution of the complete (non-radial) timing equations for the Leo I orbit implies a Milky Way mass M_MW,vir = 3.15 (-1.36, +1.58) x 10^12 Msun, with the large uncertainty dominated by cosmic scatter. In a companion paper, we compare the new observations to the properties of Leo I subhalo analogs extracted from cosmological simulations.

We combine our Hubble Space Telescope measurement of the proper motion of the Leo I dwarf spheroidal galaxy (presented in a companion paper) with the highest resolution numerical simulations of Galaxy-size dark matter halos in existence to constrain the mass of the Milky Way's dark matter halo (M_MW). Despite Leo I's large Galacto-centric space velocity (200 km/s) and distance (261 kpc), we show that it is extremely unlikely to be unbound if Galactic satellites are associated with dark matter substructure, as 99.9% of subhalos in the simulations are bound to their host. The observed position and velocity of Leo I strongly disfavor a low mass Milky Way: if we assume that Leo I is the least bound of the Milky Way's classical satellites, then we find that M_MW > 10^{12} M_sun at 95% confidence for a variety of Bayesian priors on M_MW. In lower mass halos, it is vanishingly rare to find subhalos at 261 kpc moving as fast as Leo I. Should an additional classical satellite be found to be less bound than Leo I, this lower limit on M_MW would increase by 30%. Imposing a mass weighted LCDM prior, we find a median Milky Way virial mass of M_MW=1.6 x 10^{12} M_sun, with a 90% confidence interval of [1.0-2.4] x 10^{12} M_sun. We also confirm a strong correlation between subhalo infall time and orbital energy in the simulations and show that proper motions can aid significantly in interpreting the infall times and orbital histories of satellites.

We present Keck/DEIMOS spectroscopy of resolved stars in the M31 satellites And XXVIII & And XXIX. We show that these are likely self-bound galaxies based on 18 and 24 members in And XXVIII & And XXIX, respectively. And XXVIII has a systemic velocity of -331.1 +/- 1.8 km/s and velocity dispersion of 4.9+/-1.6 km/s, implying a mass-to-light ratio (within r_1/2) of ~ 44 +/- 41. And XXIX has a systemic velocity of -194.4 +/- 1.5 km/s and velocity dispersion of 5.7+/-1.2 km/s, implying a mass-to-light ratio (within r_1/2) of ~ 124 +/- 72. The internal kinematics and implied masses of And XXVIII & And XXIX are similar to dwarf spheroidals (dSphs) of comparable luminosities, implying that these objects are dark matter-dominated dwarf galaxies. Despite the large projected distances from their host (380 and 188 kpc), the kinematics of these dSph suggest that they are bound M31 satellites.

Recently, Widrow and collaborators announced the discovery of vertical density waves in the Milky Way disk. Here we investigate a scenario where these waves were induced by the Sagittarius dwarf galaxy as it plunged through the Galaxy. Using numerical simulations, we find that the Sagittarius impact produces North-South asymmetries and vertical wave-like behavior that qualitatively agrees with what is observed. The extent to which vertical modes can radially penetrate into the disc, as well as their amplitudes, depend on the mass of the perturbing satellite. We show that the mean height of the disc is expected to vary more rapidly in the radial than in the azimuthal direction. If the observed vertical density asymmetry is indeed caused by vertical oscillations, we predict radial and azimuthal variations of the mean vertical velocity, correlating with the spatial structure. These variations can have amplitudes as large as 8 km/s.

We present the surface brightness profile of M31's stellar halo out to a projected radius of 175 kpc. The surface brightness estimates are based on confirmed samples of M31 red giant branch stars derived from Keck/DEIMOS spectroscopic observations. A set of empirical spectroscopic and photometric M31 membership diagnostics is used to identify and reject foreground and background contaminants. This enables us to trace the stellar halo of M31 to larger projected distances and fainter surface brightnesses than previous photometric studies. The surface brightness profile of M31's halo follows a power-law with index -2.2 +/- 0.2 and extends to a projected distance of at least ~175 kpc (~ 2/3 of M31's virial radius), with no evidence of a downward break at large radii. The best-fit elliptical isophotes have b/a=0.94 with the major axis of the halo aligned along the minor axis of M31's disk, consistent with a prolate halo, although the data are also consistent with M31's halo having spherical symmetry. The fact that tidal debris features are kinematically cold is used to identify substructure in the spectroscopic fields out to projected radii of 90 kpc, and investigate the effect of this substructure on the surface brightness profile. The scatter in the surface brightness profile is reduced when kinematically identified tidal debris features in M31 are statistically subtracted; the remaining profile indicates a comparatively diffuse stellar component to M31's stellar halo exists to large distances. Beyond 90 kpc, kinematically cold tidal debris features can not be identified due to small number statistics; nevertheless, the significant field-to-field variation in surface brightness beyond 90 kpc suggests that the outermost region of M31's halo is also comprised to a significant degree of stars stripped from accreted objects.

The Milky Way appears to be missing baryons, as the observed mass in stars and gas is well below the cosmic mean. One possibility is that a substantial fraction of the Galaxy's baryons are embedded within an extended, million-degree hot halo, an idea supported indirectly by observations of warm gas clouds in the halo and gas-free dwarf spheroidal satellites. X-ray observations have established that hot gas does exist in our Galaxy beyond the local hot bubble; however, it may be distributed in a hot disk configuration. Moreover, recent investigations into the X-ray constraints have suggested that any Galactic corona must be insignificant. Here we re-examine the observational data, particularly in the X-ray and radio bands, in order to determine whether it is possible for a substantial fraction of the Galaxy's baryons to exist in ~ 10^6 K gas. In agreement with past studies, we find that a baryonically closed halo is clearly ruled out if one assumes that the hot corona is distributed with a cuspy NFW profile. However, if the hot corona of the galaxy is in an extended, low-density distribution with a large central core, as expected for an adiabatic gas in hydrostatic equilibrium, then it may contain up to 10^11 M_sun of material, possibly accounting for all of the missing Galactic baryons. We briefly discuss some potential avenues for discriminating between a massive, extended hot halo and a local hot disk.

This document summarizes the results of a community-based discussion of the potential science impact of the Mayall+BigBOSS highly multiplexed multi-object spectroscopic capability. The KPNO Mayall 4m telescope equipped with the DOE- and internationally-funded BigBOSS spectrograph offers one of the most cost-efficient ways of accomplishing many of the pressing scientific goals identified for this decade by the "New Worlds, New Horizons" report. The BigBOSS Key Project will place unprecedented constraints on cosmological parameters related to the expansion history of the universe. With the addition of an open (publicly funded) community access component, the scientific impact of BigBOSS can be extended to many important astrophysical questions related to the origin and evolution of galaxies, stars, and the IGM. Massive spectroscopy is the critical missing ingredient in numerous ongoing and planned ground- and space-based surveys, and BigBOSS is unique in its ability to provide this to the US community. BigBOSS data from community-led projects will play a vital role in the education and training of students and in maintaining US leadership in these fields of astrophysics. We urge the NSF-AST division to support community science with the BigBOSS multi-object spectrograph through the period of the BigBOSS survey in order to ensure public access to the extraordinary spectroscopic capability.

We present kinematical profiles and metallicity for the M31 dwarf spheroidal (dSph) satellite galaxy Andromeda II (And II) based on Keck DEIMOS spectroscopy of 531 red giant branch stars. Our kinematical sample is among the largest for any M31 satellite and extends out to two effective radii (r_eff = 5.3' = 1.1 kpc). We find a mean systemic velocity of -192.4+-0.5 km/s and an average velocity dispersion of sigma_v = 7.8+-1.1 km/s. While the rotation velocity along the major axis of And II is nearly zero (<1 km/s), the rotation along the minor axis is significant with a maximum rotational velocity of v_max=8.6+-1.8 km/s. We find a kinematical major axis, with a maximum rotational velocity of v_max=10.9+-2.4 km/s, misaligned by 67 degrees to the isophotal major axis. And II is thus the first dwarf galaxy with evidence for nearly prolate rotation with a v_max/sigma_v = 1.1, although given its ellipticity of epsilon = 0.10, this object may be triaxial. We measured metallicities for a subsample of our data, finding a mean metallicity of [Fe/H] = -1.39+- 0.03 dex and an internal metallicity dispersion of 0.72+-0.03 dex. We find a radial metallicity gradient with metal-rich stars more centrally concentrated, but do not observe a significant difference in the dynamics of two metallicity populations. And II is the only known dwarf galaxy to show minor axis rotation making it a unique system whose existence offers important clues on the processes responsible for the formation of dSphs.

We use cosmological simulations to study the effects of self-interacting dark matter (SIDM) on the density profiles and substructure counts of dark matter halos from the scales of spiral galaxies to galaxy clusters, focusing explicitly on models with cross sections over dark matter particle mass \sigma/m = 1 and 0.1 cm^2/g. Our simulations rely on a new SIDM N-body algorithm that is derived self-consistently from the Boltzmann equation and that reproduces analytic expectations in controlled numerical experiments. We find that well-resolved SIDM halos have constant-density cores, with significantly lower central densities than their CDM counterparts. In contrast, the subhalo content of SIDM halos is only modestly reduced compared to CDM, with the suppression greatest for large hosts and small halo-centric distances. Moreover, the large-scale clustering and halo circular velocity functions in SIDM are effectively identical to CDM, meaning that all of the large-scale successes of CDM are equally well matched by SIDM. From our largest cross section runs we are able to extract scaling relations for core sizes and central densities over a range of halo sizes and find a strong correlation between the core radius of an SIDM halo and the NFW scale radius of its CDM counterpart. We construct a simple analytic model, based on CDM scaling relations, that captures all aspects of the scaling relations for SIDM halos. Our results show that halo core densities in \sigma/m = 1 cm^2/g models are too low to match observations of galaxy clusters, low surface brightness spirals (LSBs), and dwarf spheroidal galaxies. However, SIDM with \sigma/m ~ 0.1 cm^2/g appears capable of reproducing reported core sizes and central densities of dwarfs, LSBs, and galaxy clusters without the need for velocity dependence. (abridged)

If dark matter has a large self-interaction scattering cross section, then interactions among dark-matter particles will drive galaxy and cluster halos to become spherical in their centers. Work in the past has used this effect to rule out velocity-independent, elastic cross sections larger than sigma/m ~ 0.02 cm^2/g based on comparisons to the shapes of galaxy cluster lensing potentials and X-ray isophotes. In this paper, we use cosmological simulations to show that these constraints were off by more than an order of magnitude because (a) they did not properly account for the fact that the observed ellipticity gets contributions from the triaxial mass distribution outside the core set by scatterings, (b) the scatter in axis ratios is large and (c) the core region retains more of its triaxial nature than estimated before. Including these effects properly shows that the same observations now allow dark matter self-interaction cross sections at least as large as sigma/m = 0.1 cm^2/g. We show that constraints on self-interacting dark matter from strong-lensing clusters are likely to improve significantly in the near future, but possibly more via central densities and core sizes than halo shapes.

We present a resolved-star spectroscopic survey of 15 dwarf spheroidal (dSph) satellites of the Andromeda Galaxy (M31) as part of the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo (SPLASH) project. We filter foreground contamination from Milky Way (MW) stars, noting that MW substructure is evident in this contaminant sample. We also filter M31 halo field giant stars, and identify the remainder as probable dSph members. We then use these members to determine the kinematical properties of the dSphs. For the first time, we confirm that And XVIII, XXI, and XXII show kinematics consistent with bound, dark matter-dominated galaxies. From the velocity dispersions for the full sample of dSphs we determine masses, which we combine with the size and luminosity of the galaxies to produce mass-size-luminosity scaling relations. With these scalings we determine that the M31 dSphs are fully consistent with the MW dSphs, suggesting that the well-studied MW satellite population provides a fair sample for broader conclusions. We also estimate dark matter halo masses of the satellites, and find that there is no sign that the luminosity of these galaxies depends on their dark halo mass, a result consistent with what is seen for MW dwarfs. Two of the M31 dSphs (And XV, XVI) have estimated maximum circular velocities smaller than 12 km/s (to 1sigma), which likely places them within the lowest mass dark matter halos known to host stars (along with Bootes I of the MW). Finally, we use the systemic velocities of the M31 satellites to estimate the total mass of the M31 halo: within 139 kpc, M31's mass is 8 +4.1 - 3.7 x 10^11 corresponding to a virial mass for M31's dark matter halo of 1.2 +0.9 -0.7 x 10^12, consistent with previous results.

It is now known that minor mergers are capable of creating structure in the phase-space distribution of their host galaxy's disc. In order to search for such imprints in the Milky Way, we analyse the SEGUE F/G-dwarf and the Schuster et al. (2006) stellar samples. We find similar features in these two completely independent stellar samples, consistent with the predictions of a Milky Way minor-merger event. We next apply the same analyses to high-resolution, idealised N-body simulations of the interaction between the Sagittarius dwarf galaxy and the Milky Way. The energy distributions of stellar particle samples in small spatial regions in the host disc reveal strong variations of structure with position. We find good matches to the observations for models with a mass of Sagittarius' dark matter halo progenitor $\lessapprox 10^{11}$ M$_{\odot}$. Thus, we show that this kind of analysis could be used to provide unprecedentedly tight constraints on Sagittarius' orbital parameters, as well as place a lower limit on its mass.

The combination of large size, high stellar density, high metallicity, and Sersic surface brightness profile of the spheroidal component of the Andromeda galaxy (M31) within R_proj ~ 20 kpc suggest that it is unlike any subcomponent of the Milky Way. In this work we capitalize on our proximity to and external view of M31 to probe the kinematical properties of this "inner spheroid." We employ a Markov chain Monte Carlo (MCMC) analysis of resolved stellar kinematics from Keck/DEIMOS spectra of 5651 red giant branch stars to disentangle M31's inner spheroid from its stellar disk. We measure the mean velocity and dispersion of the spheroid in each of five spatial bins after accounting for a locally cold stellar disk as well as the Giant Southern Stream and associated tidal debris. For the first time, we detect significant spheroid rotation (v_rot ~ 50 km/s) beyond R_proj ~ 5 kpc. The velocity dispersion decreases from about 140 km/s at R_proj = 7 kpc to 120 km/s at R_proj = 14 kpc, consistent to 2 sigma with existing measurements and models. We calculate the probability that a given star is a member of the spheroid and find that the spheroid has a significant presence throughout the spatial extent of our sample. Lastly, we show that the flattening of the spheroid is due to velocity anisotropy in addition to rotation. Though this suggests that the inner spheroid of M31 more closely resembles an elliptical galaxy than a typical spiral galaxy bulge, it should be cautioned that our measurements are much farther out (2 - 14 r_eff) than for the comparison samples.

We use the Aquarius simulations to show that the most massive subhalos in galaxy-mass dark matter halos in LCDM are grossly inconsistent with the dynamics of the brightest Milky Way dwarf spheroidal galaxies. While the best-fitting hosts of the dwarf spheroidals all have 12 < Vmax < 25 km/s, LCDM simulations predict at least ten subhalos with Vmax > 25 km/s. These subhalos are also among the most massive at earlier times, and significantly exceed the UV suppression mass back to z ~ 10. No LCDM-based model of the satellite population of the Milky Way explains this result. The problem lies in the satellites' densities: it is straightforward to match the observed Milky Way luminosity function, but doing so requires the dwarf spheroidals to have dark matter halos that are a factor of ~5 more massive than is observed. Independent of the difficulty in explaining the absence of these dense, massive subhalos, there is a basic tension between the derived properties of the bright Milky Way dwarf spheroidals and LCDM expectations. The inferred infall masses of these galaxies are all approximately equal and are much lower than standard LCDM predictions for systems with their luminosities. Consequently, their implied star formation efficiencies span over two orders of magnitude, from 0.2% to 20% of baryons converted into stars, in stark contrast with expectations gleaned from more massive galaxies. We explore possible solutions to these problems within the context of LCDM and find them to be unconvincing. In particular, we use controlled simulations to demonstrate that the small stellar masses of the bright dwarf spheroidals make supernova feedback an unlikely explanation for their low inferred densities.

We investigate the mass distributions within eight classical Milky Way dwarf spheroidal galaxies (MW dSphs) using an equilibrium Jeans analysis and we compare our results to the mass distributions predicted for subhalos in dissipationless \Lambda CDM simulations. In order to match the dark matter density concentrations predicted, the stars in these galaxies must have a fairly significant tangential velocity dispersion anisotropy (\beta ~-1.5). For the limiting case of an isotropic velocity dispersion (\beta =0), the classical MW dSphs predominantly prefer to live in halos that are less concentrated than \Lambda CDM predictions. We also investigate whether the dSphs prefer to live in halos with constant density cores in the limit of isotropic velocity dispersion. Interestingly, even in this limit, not all of the dSphs prefer large constant-density cores: the Sculptor dSph prefers a cusp while Carina, Draco and Leo I prefer cores. The other four dSphs do not show a statistically significant preference for either cuspy or cored profiles. Finally, we re-examine the hypothesis that the density profiles of these eight MW dSphs can be quantified by a common dark matter halo.

Previous studies have shown that galaxies with minor companions exhibit an elevated star formation rate. We reverse this inquiry, constructing a volume-limited sample of \simL\star (Mr \leq -19.5 + 5 log h) galaxies from the Sloan Digital Sky Survey that are isolated with respect to other luminous galaxies. Cosmological simulations suggest that 99.8% of these galaxies are alone in their dark matter haloes with respect to other luminous galaxies. We search the area around these galaxies for photometric companions. Matching strongly star forming (EW(H{\alpha})\geq 35 \AA) and quiescent (EW(H{\alpha})< 35 \AA) samples for stellar mass and redshift using a Monte Carlo resampling technique, we demonstrate that rapidly star-forming galaxies are more likely to have photometric companions than other galaxies. The effect is relatively small; about 11% of quiescent, isolated galaxies have minor photometric companions at radii \leq 60 kpc h$^{-1}$ kpc while about 16% of strongly star-forming ones do. Though small, the cumulative difference in satellite counts between strongly star-forming and quiescent galaxies is highly statistically significant (PKS = 1.350 \times10$^{-3}$) out to to radii of \sim 100 h$^{-1}$ kpc. We discuss explanations for this excess, including the possibility that \sim 5% of strongly star-forming galaxies have star formation that is causally related to the presence of a minor companion.

We analyze subhalos in the Via Lactea II (VL2) cosmological simulation to look for correlations among their infall times and z = 0 dynamical properties. We find that the present day orbital energy is tightly correlated with the time at which subhalos last crossed into the virial radius. This energy-infall correlation provides a means to infer infall times for Milky Way satellite galaxies. Assuming that the Milky Way's assembly can be modeled by VL2, we show that the infall times of some satellites are well constrained given only their Galactocentric positions and line-of-sight velocities. The constraints sharpen for satellites with proper motion measurements. We find that Carina, Ursa Minor, and Sculptor were all accreted early, more than 8 Gyr ago. Five other dwarfs, including Sextans and Segue 1, are also probable early accreters, though with larger uncertainties. On the other extreme, Leo T is just falling into the Milky Way for the first time while Leo I fell in \sim 2 Gyr ago and is now climbing out of the Milky Way's potential after its first perigalacticon. The energies of several other dwarfs, including Fornax and Hercules, point to intermediate infall times, 2 - 8 Gyr ago. We compare our infall time estimates to published star formation histories and find hints of a dichotomy between ultrafaint and classical dwarfs. The classical dwarfs appear to have quenched star formation after infall but the ultrafaint dwarfs tend to be quenched long before infall, at least for the cases in which our uncertainties allow us to discern differences. Our analysis suggests that the Large Magellanic Cloud crossed inside the Milky Way virial radius recently, within the last \sim 4 billion years.

Like many galaxies of its size, the Milky Way is a disk with prominent spiral arms rooted in a central bar, although our knowledge of its structure and origin is incomplete. Traditional attempts to understand the Galaxy's morphology assume that it has been unperturbed by major external forces. Here we report simulations of the response of the Milky Way to the infall of the Sagittarius dwarf galaxy (Sgr), which results in the formation of spiral arms, influences the central bar and produces a flared outer disk. Two ring-like wrappings emerge towards the Galactic anti-Center in our model that are reminiscent of the low- latitude arcs observed in the same area of the Milky Way. Previous models have focused on Sgr itself to reproduce the dwarf's orbital history and place associated constraints on the shape of the Milky Way gravitational potential, treating the Sgr impact event as a trivial influence on the Galactic disk. Our results show that the Milky Way's morphology is not purely secular in origin and that low-mass minor mergers predicted to be common throughout the Universe probably have a similarly important role in shaping galactic structure.

We use cosmological SPH simulations to study the kinematic signatures of cool gas accretion onto a pair of well-resolved galaxy halos. Cold-flow streams and gas-rich mergers produce a circum-galactic component of cool gas that generally orbits with high angular momentum about the galaxy halo before falling in to build the disk. This signature of cosmological accretion should be observable using background-object absorption line studies as features that are offset from the galaxy's systemic velocity by ~100 km/s. Accreted gas typically co-rotates with the central disk in the form of a warped, extended cold flow disk, such that the observed velocity offset is in the same direction as galaxy rotation, appearing in sight lines that avoid the galactic poles. This prediction provides a means to observationally distinguish accreted gas from outflow gas: the accreted gas will show large one-sided velocity offsets in absorption line studies while radial/bi-conical outflows will not (except possibly in special polar projections). This rotation signature has already been seen in studies of intermediate redshift galaxy-absorber pairs; we suggest that these observations may be among the first to provide indirect observational evidence for cold accretion onto galactic halos. Cold mode halo gas typically has ~3-5 times more specific angular momentum than the dark matter. The associated cold mode disk configurations are likely related to extended HI/XUV disks seen around galaxies in the local universe. The fraction of galaxies with extended cold flow disks and associated offset absorption-line gas should decrease around bright galaxies at low redshift, as cold mode accretion dies out.

We introduce a comprehensive analysis of multi-epoch stellar line-of-sight velocities to determine the intrinsic velocity dispersion of the ultrafaint satellites of the Milky Way. Our method includes a simultaneous Bayesian analysis of both membership probabilities and the contribution of binary orbital motion to the observed velocity dispersion within a 14-parameter likelihood. We apply our method to the Segue 1 dwarf galaxy and conclude that Segue 1 is a dark-matter-dominated galaxy at high probability with an intrinsic velocity dispersion of 3.7^{+1.4}_{-1.1} km/sec. The dark matter halo required to produce this dispersion must have an average density of 2.5^{+4.1}_{-1.9} solar mass/pc^3 within a sphere that encloses half the galaxy's stellar luminosity. This is the highest measured density of dark matter in the Local Group. Our results show that a significant fraction of the stars in Segue 1 may be binaries with the most probable mean period close to 10 years, but also consistent with the 180 year mean period seen in the solar vicinity at about 1 sigma. Despite this binary population, the possibility that Segue 1 is a bound star cluster with the observed velocity dispersion arising from the orbital motion of binary stars is disfavored by the multi-epoch stellar velocity data at greater than 99% C.L. Finally, our treatment yields a projected (two-dimensional) half-light radius for the stellar profile of Segue 1 of 28^{+5}_{-4} pc, in excellent agreement with photometric measurements.

We present and analyze the positions, distances, and radial velocities for over 4000 blue horizontal-branch (BHB) stars in the Milky Way's halo, drawn from SDSS DR8. We search for position-velocity substructure in these data, a signature of the hierarchical assembly of the stellar halo. Using a cumulative "close pair distribution" (CPD) as a statistic in the 4-dimensional space of sky position, distance, and velocity, we quantify the presence of position-velocity substructure at high statistical significance among the BHB stars: pairs of BHB stars that are close in position on the sky tend to have more similar distances and radial velocities compared to a random sampling of these overall distributions. We make analogous mock-observations of 11 numerical halo formation simulations, in which the stellar halo is entirely composed of disrupted satellite debris, and find a level of substructure comparable to that seen in the actually observed BHB star sample. This result quantitatively confirms the hierarchical build-up of the stellar halo through a signature in phase (position-velocity) space. In detail, the structure present in the BHB stars is somewhat less prominent than that seen in most simulated halos, quite possibly because BHB stars represent an older sub-population. BHB stars located beyond 20 kpc from the Galactic center exhibit stronger substructure than at $\rm r_{gc} < 20$ kpc.

We use cosmological SPH simulations to study the cool, accreted gas in two Milky Way-size galaxies through cosmic time to z=0. We find that gas from mergers and cold flow accretion results in significant amounts of cool gas in galaxy halos. This cool circum-galactic component drops precipitously once the galaxies cross the critical mass to form stable shocks, Mvir = Msh ~ 10^12 Msun. Before reaching Msh, the galaxies experience cold mode accretion (T<10^5 K) and show moderately high covering fractions in accreted gas: f_c ~ 30-50% for R<50 co-moving kpc and N_HI>10^16 cm^-2. These values are considerably lower than observed covering fractions, suggesting that outflowing gas (not included here) is important in simulating galaxies with realistic gaseous halos. Within ~500 Myr of crossing the Msh threshold, each galaxy transitions to hot mode gas accretion, and f_c drops to ~5%. The sharp transition in covering fraction is primarily a function of halo mass, not redshift. This signature should be detectable in absorption system studies that target galaxies of varying host mass, and may provide a direct observational tracer of the transition from cold flow accretion to hot mode accretion in galaxies.

We use a volume-limited spectroscopic sample of isolated galaxies in the Sloan Digital Sky Survey (SDSS) to investigate the frequency and radial distribution of luminous (M_r <~ -18.3) satellites like the Large Magellanic Cloud (LMC) around ~L* Milky Way analogs and compare our results object-by-object to LCDM predictions based on abundance matching in simulations. We show that 12% of Milky Way-like galaxies host an LMC-like satellite within 75 kpc (projected), and 42 % within 250 kpc (projected). This implies ~10% have a satellite within the distance of the LMC, and ~40% of L* galaxies host a bright satellite within the virialized extent of their dark matter halos. Remarkably, the simulation reproduces the observed frequency, radial dependence, velocity distribution, and luminosity function of observed secondaries exceptionally well, suggesting that LCDM provides an accurate reproduction of the observed Universe to galaxies as faint as L~10^9 Lsun on ~50 kpc scales. When stacked, the observed projected pairwise velocity dispersion of these satellites is sigma~160 km/s, in agreement with abundance-matching expectations for their host halo masses. Finally, bright satellites around L* primaries are significantly redder than typical galaxies in their luminosity range, indicating that environmental quenching is operating within galaxy-size dark matter halos that typically contain only a single bright satellite. This redness trend is in stark contrast to the Milky Way's LMC, which is unusually blue even for a field galaxy. We suggest that the LMC's discrepant color might be further evidence that it is undergoing a triggered star-formation event upon first infall.

We show that dissipationless LCDM simulations predict that the majority of the most massive subhaloes of the Milky Way are too dense to host any of its bright satellites (L_V > 10^5 L_sun). These dark subhaloes have circular velocities at infall of 30-70 km/s and infall masses of [0.2-4] x 10^10 M_sun. Unless the Milky Way is a statistical anomaly, this implies that galaxy formation becomes effectively stochastic at these masses. This is in marked contrast to the well-established monotonic relation between galaxy luminosity and halo circular velocity (or halo mass) for more massive haloes. We show that at least two (and typically four) of these massive dark subhaloes are expected to produce a larger dark matter annihilation flux than Draco. It may be possible to circumvent these conclusions if baryonic feedback in dwarf satellites or different dark matter physics can reduce the central densities of massive subhaloes by order unity on a scale of 0.3 - 1 kpc.

We present the results of a comprehensive Keck/DEIMOS spectroscopic survey of the ultra-faint Milky Way satellite galaxy Segue 1. We have obtained velocity measurements for 98.2% of the stars within 67 pc (10 arcmin, or 2.3 half-light radii) of the center of Segue 1 that have colors and magnitudes consistent with membership, down to a magnitude limit of r=21.7. Based on photometric, kinematic, and metallicity information, we identify 71 stars as probable Segue 1 members, including some as far out as 87 pc. After correcting for the influence of binary stars using repeated velocity measurements, we determine a velocity dispersion of 3.7^{+1.4}_{-1.1} km/s, with a corresponding mass within the half-light radius of 5.8^{+8.2}_{-3.1} x 10^5 Msun. The stellar kinematics of Segue 1 require very high mass-to-light ratios unless the system is far from dynamical equilibrium, even if the period distribution of unresolved binary stars is skewed toward implausibly short periods. With a total luminosity less than that of a single bright red giant and a V-band mass-to-light ratio of 3400 Msun/Lsun, Segue 1 is the darkest galaxy currently known. We critically re-examine recent claims that Segue 1 is a tidally disrupting star cluster and that kinematic samples are contaminated by the Sagittarius stream. The extremely low metallicities ([Fe/H] < -3) of two Segue 1 stars and the large metallicity spread among the members demonstrate conclusively that Segue 1 is a dwarf galaxy, and we find no evidence in favor of tidal effects. We also show that contamination by the Sagittarius stream has been overestimated. Segue 1 has the highest measured dark matter density of any known galaxy and will therefore be a prime testing ground for dark matter physics and galaxy formation on small scales.

The ubiquity of substructure in the stellar halo has already been demonstrated by the SDSS and 2MASS and future surveys promise to explore the halo in ever more detail. This paper examines what can be learnt from current and future photometric-databases using group-finding techniques. We compare groups recovered from a sample of M-giants from 2MASS with those found in synthetic surveys of simulated $\Lambda$CDM stellar halos and demonstrate broad consistency. We also find that these recovered groups are likely to represent the majority of high-luminosity ($\log(L/L_{Sun})>6.67$) satellites accreted within the last 10 Gyr and on orbits with apocenters within 100 kpc. However the sensitivity of the M-giant survey to accretion events that were either ancient, from low-luminosity objects or those on radial orbits is limited because of the low number of stars, bias towards high-metallicity stars and the shallow depth. We examine the extent to which these limitations are addressed by current and future surveys, in particular catalogues of main-sequence turn-off (MSTO) stars from SDSS and LSST, and of RR-Lyrae stars from LSST or PanSTARRS. The MSTO surveys are more sensitive to low-luminosity events ($\log(L/L_{Sun}~5$ or less) than the 2MASS M-giant sample, while the RR-Lyrae surveys, with superior depth, are good at detecting events on highly eccentric orbits. When combined we expect these photometric surveys to provide a comprehensive picture of the last 10 Gyr of Galactic accretion. Events older than this are too phase mixed to be discovered and would require additional information such as velocity and chemical abundance.

We examine scaling relations of dispersion-supported galaxies over more than eight orders of magnitude in luminosity by transforming standard fundamental plane parameters into a space of mass (M1/2), radius (r1/2), and luminosity (L1/2). We find that from ultra-faint dwarf spheroidals to giant cluster spheroids, dispersion-supported galaxies scatter about a one-dimensional "fundamental curve" through this MRL space. The weakness of the M1/2-L1/2 slope on the faint end may imply that potential well depth limits galaxy formation in small galaxies, while the stronger dependence on L1/2 on the bright end suggests that baryonic physics limits galaxy formation in massive galaxies. The mass-radius projection of this curve can be compared to median dark matter halo mass profiles of LCDM halos in order to construct a virial mass-luminosity relationship (Mvir-L) for galaxies that spans seven orders of magnitude in Mvir. Independent of any global abundance or clustering information, we find that (spheroidal) galaxy formation needs to be most efficient in halos of Mvir ~ 10^12 Msun and to become inefficient above and below this scale. Moreover, this profile matching technique is most accurate at the high and low luminosity extremes (where dark matter fractions are highest) and is therefore quite complementary to statistical approaches that rely on having a well-sampled luminosity function. We also consider the significance and utility of the scatter about this relation, and find that in the dSph regime observational errors are almost at the point where we can explore the intrinsic scatter in the luminosity-virial mass relation. Finally, we note that purely stellar systems like Globular Clusters and Ultra Compact Dwarfs do not follow the fundamental curve relation. This allows them to be easily distinguished from dark-matter dominated dSph galaxies in MRL space. (abridged)

Environmental statistics provide a necessary means of comparing the properties of galaxies in different environments and a vital test of models of galaxy formation within the prevailing, hierarchical cosmological model. We explore counts-in-cylinders, a common statistic defined as the number of companions of a particular galaxy found within a given projected radius and redshift interval. Galaxy distributions with the same two-point correlation functions do not necessarily have the same companion count distributions. We use this statistic to examine the environments of galaxies in the Sloan Digital Sky Survey, Data Release 4. We also make preliminary comparisons to four models for the spatial distributions of galaxies, based on N-body simulations, and data from SDSS DR4 to study the utility of the counts-in-cylinders statistic. There is a very large scatter between the number of companions a galaxy has and the mass of its parent dark matter halo and the halo occupation, limiting the utility of this statistic for certain kinds of environmental studies. We also show that prevalent, empirical models of galaxy clustering that match observed two- and three-point clustering statistics well fail to reproduce some aspects of the observed distribution of counts-in-cylinders on 1, 3 and 6-Mpc/h scales. All models that we explore underpredict the fraction of galaxies with few or no companions in 3 and 6-Mpc/h cylinders. Roughly 7% of galaxies in the real universe are significantly more isolated within a 6 Mpc/h cylinder than the galaxies in any of the models we use. Simple, phenomenological models that map galaxies to dark matter halos fail to reproduce high-order clustering statistics in low-density environments.

The Missing Satellites Problem (MSP) broadly refers to the overabundance of predicted Cold Dark Matter (CDM) subhalos compared to satellite galaxies known to exist in the Local Group. The most popular interpretation of the MSP is that the smallest dark matter halos in the universe are extremely inefficient at forming stars. The question from that standpoint is to identify the feedback source that makes small halos dark and to identify any obvious mass scale where the truncation in the efficiency of galaxy formation occurs. Among the most exciting developments in near-field cosmology in recent years is the discovery of a new population satellite galaxies orbiting the Milky Way and M31. Wide field, resolved star surveys have more than doubled the dwarf satellite count in less than a decade, revealing a population of ultrafaint galaxies that are less luminous that some star clusters. For the first time, there are empirical reasons to believe that there really are missing satellite galaxies in the Local Group, lurking just beyond our ability to detect them, or simply inhabiting a region of the sky that has yet to have been surveyed. Both kinematic studies and completeness-correction studies seem to point to a characteristic potential well depth for satellite subhalos that is quite close to the mass scale where photoionization and atomic cooling should limit galaxy formation. Among the more pressing problems associated with this interpretation is to understand the selection biases that limit our ability to detect the lowest mass galaxies. The least massive satellite halos are likely to host stealth galaxies with very-low surface brightness and this may be an important limitation in the hunt for low-mass fossils from the epoch of reionization.

We derive an accurate mass estimator for dispersion-supported stellar systems and demonstrate its validity by analyzing resolved line-of-sight velocity data for globular clusters, dwarf galaxies, and elliptical galaxies. Specifically, by manipulating the spherical Jeans equation we show that the dynamical mass enclosed within the 3D deprojected half-light radius r_1/2 can be determined with only mild assumptions about the spatial variation of the stellar velocity dispersion anisotropy. We find M_1/2 = 3 \sigma_los^2 r_1/2 / G ~ 4 \sigma_los^2 R_eff / G, where \sigma_los^2 is the luminosity-weighted square of the line-of-sight velocity dispersion and R_eff is the 2D projected half-light radius. While deceptively familiar in form, this formula is not the virial theorem, which cannot be used to determine accurate masses unless the radial profile of the total mass is known a priori. We utilize this finding to show that all of the Milky Way dwarf spheroidal galaxies (MW dSphs) are consistent with having formed within a halo of mass approximately 3 x 10^9 M_sun in Lambda CDM cosmology. The faintest MW dSphs seem to have formed in dark matter halos that are at least as massive as those of the brightest MW dSphs, despite the almost five orders of magnitude spread in luminosity. We expand our analysis to the full range of observed dispersion-supported stellar systems and examine their I-band mass-to-light ratios (M/L). The M/L vs. M_1/2 relation for dispersion-supported galaxies follows a U-shape, with a broad minimum near M/L ~ 3 that spans dwarf elliptical galaxies to normal ellipticals, a steep rise to M/L ~ 3,200 for ultra-faint dSphs, and a more shallow rise to M/L ~ 800 for galaxy cluster spheroids.

A density based hierarchical group-finding algorithm is used to identify stellar halo structures in a catalog of M-giants from the Two Micron All Sky Survey (2MASS). The intrinsic brightness of M-giant stars means that this catalog probes deep into the halo where substructures are expected to be abundant and easy to detect. Our analysis reveals 16 structures at high Galactic latitude (greater than 15 degree), of which 10 have been previously identified. Among the six new structures two could plausibly be due to masks applied to the data, one is associated with a strong extinction region and one is probably a part of the Monoceros ring. Another one originates at low latitudes, suggesting some contamination from disk stars, but also shows protrusions extending to high latitudes, implying that it could be a real feature in the stellar halo. The last remaining structure is free from the defects discussed above and hence is very likely a satellite remnant. Although the extinction in the direction of the structure is very low, the structure does match a low temperature feature in the dust maps. While this casts some doubt on its origin, the low temperature feature could plausibly be due to real dust in the structure itself. The angular position and distance of this structure encompass the Pisces overdensity traced by RR Lyraes in Stripe 82 of the Sloan Digital Sky Survey (SDSS). However, the 2MASS M-giants indicate that the structure is much more extended than what is visible with the SDSS, with the point of peak density lying just outside Stripe 82. The morphology of the structure is more like a cloud than a stream and reminiscent of that seen in simulations of satellites disrupting along highly eccentric orbits.

We predict that there is a population of low-luminosity dwarf galaxies orbiting within the halo of the Milky Way that have surface brightnesses low enough to have escaped detection in star-count surveys. The overall count of stealth galaxies is sensitive to the presence (or lack) of a low-mass threshold in galaxy formation. These systems have luminosities and stellar velocity dispersions that are similar to those of known ultrafaint dwarf galaxies but they have more extended stellar distributions (half light radii greater than about 100 pc) because they inhabit dark subhalos that are slightly less massive than their higher surface brightness counterparts. As a result, the typical peak surface brightness is fainter than 30 mag per square arcsec. One implication is that the inferred common mass scale for Milky Way dwarfs may be an artifact of selection bias. If there is no sharp threshold in galaxy formation at low halo mass, then ultrafaint galaxies like Segue 1 represent the high-mass, early forming tail of a much larger population of objects that could number in the hundreds and have typical peak circular velocities of about 8 km/s and masses within 300 pc of about 5 million solar masses. Alternatively, if we impose a low-mass threshold in galaxy formation in order to explain the unexpectedly high densities of the ultrafaint dwarfs, then we expect only a handful of stealth galaxies in the halo of the Milky Way. A complete census of these objects will require deeper sky surveys, 30m-class follow-up telescopes, and more refined methods to identify extended, self-bound groupings of stars in the halo.

Warm dark matter (WDM) and self-interacting dark matter (SIDM) are often motivated by the inferred cores in the dark matter halos of low surface brightness (LSB) galaxies. We test thermal WDM, non-thermal WDM, and SIDM using high-resolution rotation curves of nine LSB galaxies. We fit these dark matter models to the data and determine the halo core radii and central densities. While the minimum core size in WDM models is predicted to decrease with halo mass, we find that the inferred core radii increase with halo mass and also cannot be explained with a single value of the primordial phase space density. Moreover, if the core size is set by WDM particle properties, then even the smallest cores we infer would require primordial phase space density values that are orders of magnitude smaller than lower limits obtained from the Lyman alpha forest power spectra. We also find that the dark matter halo core densities vary by a factor of about 30 from system to system while showing no systematic trend with the maximum rotation velocity of the galaxy. This strongly argues against the core size being directly set by large self-interactions (scattering or annihilation) of dark matter. We therefore conclude that the inferred cores do not provide motivation to prefer WDM or SIDM over other dark matter models.

Minor accretion events with mass ratio M_sat : M_host ~ 1:10 are common in the context of LCDM cosmology. We use high-resolution simulations of Galaxy-analogue systems to show that these mergers can dynamically eject disk stars into a diffuse light component that resembles a stellar halo both spatially and kinematically. For a variety of orbital configurations, we find that ~3-5e8 M_sun of primary stellar disk material is ejected to a distance larger than 5 kpc above the galactic plane. This ejected contribution is similar to the mass contributed by the tidal disruption of the satellite galaxy itself, though it is less extended. If we restrict our analysis to the approximate solar neighborhood in the disk plane, we find that ~1% of the initial disk stars in that region would be classified kinematically as halo stars. Our results suggest that the inner parts of galactic stellar halos contain ancient disk stars and that these stars may have been liberated in the very same events that delivered material to the outer stellar halo.

We report on the observations of 14 dwarf spheroidal galaxies with the Fermi Gamma-Ray Space Telescope taken during the first 11 months of survey mode operations. The Fermi telescope provides a new opportunity to test particle dark matter models through the expected gamma-ray emission produced by pair annihilation of weakly interacting massive particles (WIMPs). Local Group dwarf spheroidal galaxies, the largest galactic substructures predicted by the cold dark matter scenario, are attractive targets for such indirect searches for dark matter because they are nearby and among the most extreme dark matter dominated environments. No significant gamma-ray emission was detected above 100 MeV from the candidate dwarf galaxies. We determine upper limits to the gamma-ray flux assuming both power-law spectra and representative spectra from WIMP annihilation. The resulting integral flux above 100 MeV is constrained to be at a level below around 10^-9 photons cm^-2 s^-1. Using recent stellar kinematic data, the gamma-ray flux limits are combined with improved determinations of the dark matter density profile in 8 of the 14 candidate dwarfs to place limits on the pair annihilation cross-section of WIMPs in several widely studied extensions of the standard model. With the present data, we are able to rule out large parts of the parameter space where the thermal relic density is below the observed cosmological dark matter density and WIMPs (neutralinos here) are dominantly produced non-thermally, e.g. in models where supersymmetry breaking occurs via anomaly mediation. The gamma-ray limits presented here also constrain some WIMP models proposed to explain the Fermi and PAMELA e^+e^- data, including low-mass wino-like neutralinos and models with TeV masses pair-annihilating into muon-antimuon pairs. (Abridged)

We show that measured velocity dispersions of dwarf spheroidal galaxies from about 4 to 10 km/s are unlikely to be inflated by more than 20% due to the orbital motion of binary stars, and demonstrate that the intrinsic velocity dispersions can be determined to within a few percent accuracy using two-epoch observations with 1-2 years as the optimal time interval. The crucial observable is the threshold fraction--the fraction of stars that show velocity changes larger than a given threshold between measurements. The threshold fraction is tightly correlated with the dispersion introduced by binaries, independent of the underlying binary fraction and distribution of orbital parameters. We outline a simple procedure to correct the velocity dispersion to within a few percent accuracy by using the threshold fraction and provide fitting functions for this method. We also develop a methodology for constraining properties of binary populations from both single- and two-epoch velocity measurements by including the binary velocity distribution in a Bayesian analysis.

To date, the identification of interactions at z ~ 3 and above has relied on morphological analysis. Here, we present five serendipitous spectroscopic z ~ 3 Lyman break galaxy (LBG) pairs with projected proper separations < 15 h^-1 kpc in our survey of nine separate Keck fields. The data consist of 140 of our highest signal-to-noise ratio LBG spectra and ~500 of our most confident colour-selected LBGs. We show that the pairs are composed of two distinct close and/or interacting LBGs from a detailed analysis of the rest-frame ultraviolet spectra and images. In addition, we show that the pair number and separation distribution is expected from (1) the angular correlation function when applied to our survey and ~2500 colour-selected LBGs from the literature and (2) an analysis of a carefully matched high-resolution hybrid numerical and analytical cosmological simulation. Because the spectroscopic slitlets have random orientations with respect to the close pairs on the sky, the serendipitous pairs provide an unbiased sampling of the underlying close pair fraction. Finally, we discover two Ly-a emitters (LAEs) in our slitlets and find that they reside within 50 projected h^-1 kpc of LBGs. In this work, we uncover a strong relationship between Ly-a emission and pair separation. All confirmed and all candidate LBG pairs with separations of < 15 projected h^-1 kpc exhibit Ly-a in emission and we find an overabundance of Ly-a emission in pairs with < 50 projected h^-1 kpc separations. This relationship suggests a picture in which a measurable fraction of the Ly-a emission of LBGs, and potentially LAEs, is generated via interaction mechanisms. As a result, serendipitous spectroscopic close pairs provide a unique means to help identify and study high-redshift galaxy interactions using conventional ground-based optical data.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

We employ a high-resolution LCDM N-body simulation to present merger rate predictions for dark matter halos and investigate how common merger-related observables for galaxies--such as close pair counts, starburst counts, and the morphologically disturbed fraction--likely scale with luminosity, stellar mass, merger mass ratio, and redshift from z=0 to z=4. We provide a simple 'universal' fitting formula that describes our derived merger rates for dark matter halos a function of dark halo mass, merger mass ratio, and redshift, and go on to predict galaxy merger rates using number density-matching to associate halos with galaxies. For example, we find that the instantaneous merger rate of m/M>0.3 mass ratio events into typical L > f L* galaxies follows the simple relation dN/dt=0.03(1+f)(1+z)^2.1 Gyr^-1. Despite the rapid increase in merger rate with redshift, only a small fraction of >0.4 L* high-redshift galaxies (~3% at z=2) should have experienced a major merger (m/M >0.3) in the very recent past (t< 100 Myr). This suggests that short-lived, merger-induced bursts of star formation should not contribute significantly to the global star formation rate at early times, in agreement with observational indications. We emphasize that great care must be made in comparisons to observations because the predicted observables depend very sensitively on galaxy luminosity, redshift, overall mass ratio, and uncertain relaxation timescales for merger remnants. We show that the majority of bright galaxies at z=3 should have undergone a major merger (>0.3) in the last 700 Myr and conclude that mergers almost certainly play an important role in delivering baryons and influencing the kinematic properties of Lyman Break Galaxies (LBGs). (abridged)

We use N-body simulations and observationally-normalized relations between dark matter halo mass, stellar mass, and cold gas mass to derive robust, arguably inevitable expectations about the baryonic content of major mergers out to redshift z~2. First, we find that the majority of major mergers (m/M > 0.3) experienced by Milky Way size dark matter halos should have been gas-rich, and that gas-rich mergers are increasingly common at high redshift. Though the frequency of major mergers into galaxy halos in our simulations greatly exceeds the observed late-type galaxy fraction, the frequency of gas-poor major mergers is consistent with the observed fraction of spheroid-dominated galaxies across the halo mass range M_DM ~ 10^11-10^13 Msun. These results lend support to the conjecture that mergers with high baryonic gas fractions play an important role in building and/or preserving disk galaxies in the universe. Also, we find that the overall fraction of a galaxy's cold baryons deposited directly via major mergers is substantial. Approximately ~30% of the cold baryonic material in M_star ~ 10^10 Msun$ (M_DM ~ 10^11.5 Msun) galaxies is accreted as cold gas in major mergers. For more massive galaxies with M_star ~ 10^11 Msun (M_DM ~ 10^13 Msun) the fraction of baryons amassed in mergers is even higher, ~50%, but most of these accreted baryons are delivered directly in the form of stars. This baryonic mass deposition is almost unavoidable, and provides a limit on the fraction of a galaxy's cold baryons that can originate in cold flows or from hot halo cooling. (Abridged)

Modern redshift surveys enable the identification of large samples of galaxies in pairs, taken from many different environments. Meanwhile, cosmological simulations allow a detailed understanding of the statistical properties of the selected pair samples. Using these tools in tandem leads to a quantitative understanding of the effects of galaxy-galaxy interactions and, separately, the effects quenching processes in the environments of even very small groups. In the era of the next generation of large telescopes, detailed studies of interactions will be enabled to much higher redshifts.

Supernovae have been confirmed to redshift z ~ 1.7 for type Ia (thermonuclear detonation of a white dwarf) and to z ~ 0.7 for type II (collapse of the core of the star). The subclass type IIn supernovae are luminous core-collapse explosions of massive stars and, unlike other types, are very bright in the ultraviolet, which should enable them to be found optically at redshifts z ~ 2 and higher. In addition, the interaction of the ejecta with circumstellar material creates strong, long-lived emission lines that allow spectroscopic confirmation of many events of this type at z ~ 2 for 3 - 5 years after explosion. Here we report three spectroscopically confirmed type IIn supernovae, at redshifts z = 0.808, 2.013 and 2.357, detected in archival data using a method designed to exploit these properties at z ~ 2. Type IIn supernovae directly probe the formation of massive stars at high redshift. The number found to date is consistent with the expectations of a locally measured stellar initial mass function, but not with an evolving initial mass function proposed to explain independent observations at low and high redshift.

Massive satellite accretions onto early galactic disks can lead to the deposition of dark matter in disk-like configurations that co-rotate with the galaxy. This phenomenon has potentially dramatic consequences for dark matter detection experiments. We utilize focused, high-resolution simulations of accretion events onto disks designed to be Galaxy analogues, and compare the resultant disks to the morphological and kinematic properties of the Milky Way's thick disk in order to bracket the range of co-rotating accreted dark matter. We find that the Milky Way's merger history must have been unusually quiescent compared to median LCDM expectations and therefore its dark disk must be relatively small: the fraction of accreted dark disk material near the Sun is about 20% of the host halo density or smaller and the co-rotating dark matter fraction near the Sun, defined as particles moving with a rotational velocity lag less than 50 km/s, is enhanced by about 30% or less compared to a standard halo model. Such a dark disk could contribute dominantly to the low energy (of order keV for a dark matter particle with mass 100 GeV) nuclear recoil event rate of direct dectection experiments, but it will not change the likelihood of detection significantly. These dark disks provide testable predictions of weakly-interacting massive particle dark matter models and should be considered in detailed comparisons to experimental data. Our findings suggest that the dark disk of the Milky Way may provide a detectable signal for indirect detection experiments, contributing up to about 25% of the dark matter self-annihilation signal in the direction of the center of the Galaxy, lending the signal a noticeably oblate morphology.

(Abridged) We perform dissipationless N-body simulations to elucidate the dynamical response of thin disks to bombardment by cold dark matter (CDM) substructure. Our method combines (1) cosmological simulations of the formation of Milky Way (MW)-sized CDM halos to derive the properties of substructure and (2) controlled numerical experiments of consecutive subhalo impacts onto an initially-thin, fully-formed MW type disk galaxy. The present study is the first to account for the evolution of satellite populations over cosmic time in such an investigation of disk structure. We find that accretions of massive subhalos onto the central regions of host halos, where the galactic disks reside, since z~1 should be common. One host halo accretion history is used to initialize the controlled simulations of satellite-disk encounters. We show that these accretion events severely perturb the thin galactic disk and produce a wealth of distinctive dynamical signatures on its structure and kinematics. These include (1) considerable thickening and heating at all radii, with the disk thickness and velocity ellipsoid nearly doubling at the solar radius; (2) prominent flaring associated with an increase in disk thickness greater than a factor of 4 in the disk outskirts; (3) surface density excesses at large radii, beyond ~5 disk scale lengths, resembling those of observed antitruncated disks; (4) lopsidedness at levels similar to those measured in observational samples of disk galaxies; and (5) substantial tilting. The interaction with the most massive subhalo drives the disk response while subsequent bombardment is much less efficient at disturbing the disk. We conclude that substructure-disk encounters of the kind expected in the LCDM paradigm play a significant role in setting the structure of disk galaxies and driving galaxy evolution.

We present a general methodology for determining the gamma-ray flux from annihilation of dark matter particles in Milky Way satellite galaxies, focusing on two promising satellites as examples: Segue 1 and Draco. We use the SuperBayeS code to explore the best-fitting regions of the Constrained Minimal Supersymmetric Standard Model (CMSSM) parameter space, and an independent MCMC analysis of the dark matter halo properties of the satellites using published radial velocities. We present a formalism for determining the boost from halo substructure in these galaxies and show that its value depends strongly on the extrapolation of the concentration-mass (c(M)) relation for CDM subhalos down to the minimum possible mass. We show that the preferred region for this minimum halo mass within the CMSSM with neutralino dark matter is ~10^-9-10^-6 solar masses. For the boost model where the observed power-law c(M) relation is extrapolated down to the minimum halo mass we find average boosts of about 20, while the Bullock et al (2001) c(M) model results in boosts of order unity. We estimate that for the power-law c(M) boost model and photon energies greater than a GeV, the Fermi space-telescope has about 20% chance of detecting a dark matter annihilation signal from Draco with signal-to-noise greater than 3 after about 5 years of observation.

We present a pair of high-resolution smoothed particle hydrodynamics (SPH) simulations that explore the evolution and cooling behavior of hot gas around Milky-Way size galaxies. The simulations contain the same total baryonic mass and are identical other than their initial gas density distributions. The first is initialised with a low entropy hot gas halo that traces the cuspy profile of the dark matter, and the second is initialised with a high-entropy hot halo with a cored density profile as might be expected in models with pre-heating feedback. Galaxy formation proceeds in dramatically different fashion depending on the initial setup. While the low-entropy halo cools rapidly, primarily from the central region, the high-entropy halo is quasi-stable for ~4 Gyr and eventually cools via the fragmentation and infall of clouds from ~100 kpc distances. The low-entropy halo's X-ray surface brightness is ~100 times brighter than current limits and the resultant disc galaxy contains more than half of the system's baryons. The high-entropy halo has an X-ray brightness that is in line with observations, an extended distribution of pressure-confined clouds reminiscent of observed populations, and a final disc galaxy that has half the mass and ~50% more specific angular momentum than the disc formed in the low-entropy simulation. The final high-entropy system retains the majority of its baryons in a low-density hot halo. The hot halo harbours a trace population of cool, mostly ionised, pressure-confined clouds that contain ~10% of the halo's baryons after 10 Gyr of cooling. The covering fraction for HI and MgII absorption clouds in the high-entropy halo is ~0.4 and ~0.6, respectively, although most of the mass that fuels disc growth is ionised, and hence would be under counted in HI surveys.

Most Galaxy-sized systems (M_host ~ 10^12 M_sun) in the LCDM cosmology are expected to have accreted at least one satellite with a total mass M_sat ~ 10^11 M_sun = 3M_disk in the past 8 Gyr. Analytic and numerical investigations suggest that this is the most precarious type of merger for the survival of thin galactic disks because more massive accretion events are relatively rare and less massive ones preserve thin disk components. We use high-resolution, dissipationless N-body simulations to study the response of an initially-thin, fully-formed Milky-Way type stellar disk to these cosmologically common events and show that the thin disk does not survive. Regardless of orbital configuration, the impacts transform the disks into structures that are roughly three times as thick and more than twice as kinematically hot as the observed dominant thin disk component of the Milky Way. We conclude that if the Galactic thin disk is a representative case, then the presence of a stabilizing gas component is the only recourse for explaining the preponderance of disk galaxies in an LCDM universe; otherwise, the disk of the Milky Way must be uncommonly cold and thin for its luminosity, perhaps as a consequence of an unusually quiescent accretion history.

Over the past five years, searches in Sloan Digital Sky Survey data have more than doubled the number of known dwarf satellite galaxies of the Milky Way, and have revealed a population of ultra-faint galaxies with luminosities smaller than typical globular clusters, L ~ 1000 Lsun. These systems are the faintest, most dark matter dominated, and most metal poor galaxies in the universe. Completeness corrections suggest that we are poised on the edge of a vast discovery space in galaxy phenomenology, with hundreds more of these extreme galaxies to be discovered as future instruments hunt for the low-luminosity threshold of galaxy formation. Dark matter dominated dwarfs of this kind probe the small-scale power-spectrum, provide the most stringent limits on the phase-space packing of dark matter, and offer a particularly useful target for dark matter indirect detection experiments. Full use of dwarfs as dark matter laboratories will require synergy between deep, large-area photometric searches; spectroscopic and astrometric follow-up with next-generation optical telescopes; and subsequent observations with gamma-ray telescopes for dark matter indirect detection.

There is a vast menagerie of plausible candidates for the constituents of dark matter, both within and beyond extensions of the Standard Model of particle physics. Each of these candidates may have scattering (and other) cross section properties that are consistent with the dark matter abundance, BBN, and the most scales in the matter power spectrum; but which may have vastly different behavior at sub-galactic "cutoff" scales, below which dark matter density fluctuations are smoothed out. The only way to quantitatively measure the power spectrum behavior at sub-galactic scales at distances beyond the local universe, and indeed over cosmic time, is through probes available in multiply imaged strong gravitational lenses. Gravitational potential perturbations by dark matter substructure encode information in the observed relative magnifications, positions, and time delays in a strong lens. Each of these is sensitive to a different moment of the substructure mass function and to different effective mass ranges of the substructure. The time delay perturbations, in particular, are proving to be largely immune to the degeneracies and systematic uncertainties that have impacted exploitation of strong lenses for such studies. There is great potential for a coordinated theoretical and observational effort to enable a sophisticated exploitation of strong gravitational lenses as direct probes of dark matter properties. This opportunity motivates this white paper, and drives the need for: a) strong support of the theoretical work necessary to understand all astrophysical consequences for different dark matter candidates; and b) tailored observational campaigns, and even a fully dedicated mission, to obtain the requisite data.

The hierarchical theory of galaxy formation rests on the idea that smaller galactic structures merge to form the galaxies that we see today. The past decade has provided remarkable observational support for this scenario, driven in part by advances in spectroscopic instrumentation. Multi-object spectroscopy enabled the discovery of kinematically cold substructures around the Milky Way and M31 that are likely the debris of disrupting satellites. Improvements in high-resolution spectroscopy have produced key evidence that the abundance patterns of the Milky Way halo and its dwarf satellites can be explained by Galactic chemical evolution models based on hierarchical assembly. These breakthroughs have depended almost entirely on observations of nearby stars in the Milky Way and luminous red giant stars in M31 and Local Group dwarf satellites. In the next decade, extremely large telescopes will allow observations far down the luminosity function in the known dwarf galaxies, and they will enable observations of individual stars far out in the Galactic halo. The chemical abundance census now available for the Milky Way will become possible for our nearest neighbor, M31. Velocity dispersion measurements now available in M31 will become possible for systems beyond the Local Group such as Sculptor and M81 Group galaxies. Detailed studies of a greater number of individual stars in a greater number of spiral galaxies and their satellites will test hierarchical assembly in new ways because dynamical and chemical evolution models predict different outcomes for halos of different masses in different environments.

The fact that dark matter (DM), thus far, has revealed itself only on scales of galaxies and larger, again thrusts onto astrophysics the opportunity and the responsibility to confront the age old mystery "What is the nature of matter?" By deriving basic data on the nature of DM - e.g., mass of its particle(s), present mean temperature, distribution in galaxies and other structures in the universe, and capacity for dissipational collapse - we will be uncovering the properties of the dominant species of matter in the universe and significantly extending the standard models of particle physics. Determining the mass of the DM particle to an order of magnitude would help to sort out the particle family to which it (or they) belongs. Beyond mass, there are issues of stability. The DM particle may be unstable with a measurable half-life, or it may become unstable after absorbing a certain amount of energy from collisions. In both cases it would contribute to the present hot dark matter component. Some key parameters of DM can most accurately be measured in the very nearby universe because DM dominates the mass in the outer Milky Way (MW), in other galaxies in the Local Group, and in the Local Group in its entirety. The presence and amount of DM can be quantified by study of dynamical processes observable in fine detail within these entities. Precise measurements of 3-D velocities for stars, coherent star streams, and stars in satellite stellar systems out to the edge of the Galaxy can reveal "what is the shape, orientation, density law, and lumpiness of the dark matter halo" as well as "what is the total mass of the Galaxy?"

The full spatio-chemo-dynamical structure of galaxies of all types and environments at low redshift provides a critical accompaniment to observations of galaxy formation at high redshift. The next decade brings the observational opportunity to strongly constrain nearby galaxies' histories of star formation and assembly, especially in the spheroids that comprise the large majority of the stellar mass in the Universe but have until now been difficult to study. In order to constrain the pathways to building up the spheroidal "red-sequence", various standard techniques in photometry and spectroscopy, particularly with resolved tracer populations like globular clusters and planetary nebulae, can be scaled up to comprehensive surveys as improved wide-field instrumentation is increasingly available. At the same time, progress in adaptive optics on giant telescopes could for the first time permit deep, resolved photometric and spectroscopic analysis of large samples of individual stars in these systems, thereby revolutionizing galaxy studies. Strong theoretical support is needed in order to understand the new observational constraints via detailed modeling and self-consistent simulations of star and galaxy formation throughout cosmic time.

The concordance Lambda Cold Dark Matter (Lambda-CDM) model for the formation of structure in the Universe, while remarkably successful at describing observations of structure on large scales, continues to be challenged by observations on galactic scales. Fortunately, CDM models and their various proposed alternatives make a rich variety of testable predictions that make the Local Group and our own Milky Way Galaxy key laboratories for exploring dark matter (DM) in this regime. However, some of the most definitive tests of local DM require microarcsecond astrometry of faint sources, an astrometric regime that is a unique niche of SIM Lite. This chapter explores the important and distinct contributions that can be made by SIM Lite in the exploration of galaxy dynamics and DM on galaxy scales and that have cosmological consequences. Key areas of potential SIM Lite exploration include (1) measuring the shape, orientation, density law, and lumpiness of the dark halo of the Milky Way and other nearby galaxies, (2) determining the orbits of Galactic satellites, which may be representatives of late infall from the hierarchical formation of the Milky Way, (3) ascertaining the distribution of angular momentum and orbital anisotropy of stars and globular clusters to the outer reaches of the Galactic halo, dynamical properties that hold clues to the early hierarchical formation of the Galaxy, (4) measuring the physical nature of DM by placing strong constraints on the phase space density in the cores of nearby dSph galaxies, and (5) reconstructing the dynamical history of the Local Group through the determination of orbits and masses of its constituent galaxies.

Disk galaxies are common in our universe and this is a source of concern for hierarchical formation models like LCDM. Here we investigate this issue as motivated by raw merger statistics derived for galaxy-size dark matter halos from LCDM simulations. Our analysis shows that a majority (~70%) of galaxy halos with M = 10^12 M_sun at z=0 should have accreted at least one object with mass m > 10^11 M_sun ~ 3 M_disk over the last 10 Gyr. Mergers involving larger objects should have been rare, and this pinpoints m/M = 0.1 mass-ratio mergers as the most worrying ones for the survival of thin galactic disks. Motivated by these results, we use use high-resolution, dissipationless N-body simulations to study the response of stellar Milky-Way type disks to these common mergers and show that thin disks do not survive the bombardment. The remnant galaxies are roughly three times as thick and twice as kinematically hot as the observed thin disk of the Milky Way. Finally, we evaluate the suggestion that disks may be preserved if the mergers involve gas-rich progenitors. Using empirical measures to assign stellar masses and gas masses to dark matter halos as a function of redshift, we show that the vast majority of large mergers experienced by 10^12 M_sun halos should be gas-rich (f_gas > 0.5), suggesting that this is a potentially viable solution to the disk formation conundrum. Moreover, gas-rich mergers should become increasingly rare in more massive halos > 10^12.5 M_sun, and this suggest that merger gas fractions may play an important role in establishing morphological trends with galaxy luminosity.

We correct the observed Milky Way satellite luminosity function for luminosity bias using published completeness limits for the Sloan Digital Sky Survey DR5. Assuming that the spatial distribution of Milky Way satellites tracks the subhalos found in the Via Lactea LCDM N-body simulation, we show that there should be between ~300 and ~600 satellites within 400 kpc of the Sun that are brighter than the faintest known dwarf galaxies, and that there may be as many as ~ 1000, depending on assumptions. By taking into account completeness limits, we show that the radial distribution of known Milky Way dwarfs is consistent with our assumption that the full satellite population tracks that of subhalos. These results alleviate the primary worries associated with the so-called missing satellites problem in CDM. We show that future, deep wide-field surveys such as SkyMapper, the Dark Energy Survey (DES), PanSTARRS, and the Large Synoptic Survey Telescope (LSST) will deliver a complete census of ultra-faint dwarf satellites out to the Milky Way virial radius, offer new limits on the free-streaming scale of dark matter, and provide unprecedented constraints on the low-luminosity threshold of galaxy formation.

We make predictions for the metallicity of diffuse stellar components in systems ranging from small spiral galaxies to rich galaxy clusters. We extend the formalism of Purcell et al. (2007), in which diffuse stellar mass is produced via galaxy disruption, and we convolve this result with the observed mass-metallicity relation for galaxies in order to analyze the chemical abundance of intrahalo light (IHL) in host halos with virial mass 10^10.5 M_sun < M_host < 10^15 M_sun. We predict a steep rise of roughly two dex in IHL metallicity from the scales of small to large spiral galaxies. In terms of the total dynamical mass M_host of the host systems under consideration, we predict diffuse light metallicities ranging from Z_IHL < -2.5 for M_host ~ 10^11 M_sun, to Z_IHL ~ -1.0 for M_host ~ 10^12 M_sun. In larger systems, we predict a gradual flattening of this trend with Z_IHL ~ -0.4 for M_host ~ 10^13 M_sun, increasing to Z_IHL ~ 0.1 for M_host ~ 10^15 M_sun. This behavior is coincident with a narrowing of the intrahalo metallicity distribution as host mass increases. The observable distinction in surface brightness between old, metal-poor IHL stars and more metal-rich, dynamically-younger tidal streams is of crucial importance when estimating the chemical abundance of an intrahalo population with multiple origins.

The Milky Way has at least twenty-three known satellite galaxies that shine with luminosities ranging from about a thousand to a billion times that of the Sun. Half of these galaxies were discovered in the past few years in the Sloan Digital Sky Survey, and they are among the least luminous galaxies in the known Universe. A determination of the mass of these galaxies provides a test of galaxy formation at the smallest scales and probes the nature of the dark matter that dominates the mass density of the Universe. Here we use new measurements of the velocities of the stars in these galaxies to show that they are consistent with them having a common mass of about 10^7 M_\odot within their central 300 parsecs. This result demonstrates that the faintest of the Milky Way satellites are the most dark matter-dominated galaxies known, and could be a hint of a new scale in galaxy formation or a characteristic scale for the clustering of dark matter.

Integral field spectroscopy of galaxies at redshift z~2 has revealed a population of early-forming, rotationally-supported disks. These high-redshift systems provide a potentially important clue to the formation processes that build disk galaxies in the universe. A particularly well-studied example is the z=2.38 galaxy BzK-15504, which was shown by Genzel et al. (2006) to be a rotationally supported disk despite the fact that its high star formation rate and short gas consumption timescale require a very rapid acquisition of mass. Previous kinematical analyses have suggested that z~2 disk galaxies like BzK-15504 did not form through mergers because their line-of-sight velocity fields display low levels of asymmetry. We perform the same kinematical analysis on a set of simulated disk galaxies formed in gas-rich mergers of the type that may be common at high redshift, and show that the remnant disks display low velocity field asymmetry and satisfy the criteria that have been used to classify high-redshift galaxies as disks observationally. Further, we compare one of our remnants to the bulk properties of BzK-15504 and show that it has a star formation rate, gas surface density, and a circular velocity-to-velocity dispersion ratio that matches BzK-15504 remarkably well. We suggest that observations of high-redshift disk galaxies like BzK-15504 are consistent with the hypothesis that gas-rich mergers play an important role in disk formation at high redshift.

This paper explores the mapping between the observable properties of a stellar halo in phase- and abundance-space and the parent galaxy's accretion history in terms of the characteristic epoch of accretion and mass and orbits of progenitor objects. The study utilizes a suite of eleven stellar halo models constructed within the context of a standard LCDM cosmology. The results demonstrate that coordinate-space studies are sensitive to the recent (0-8 Gyears ago) merger histories of galaxies (this timescale corresponds to the last few to tens of percent of mass accretion for a Milky-Way-type galaxy). Specifically, the {\it frequency, sky coverage} and {\it fraction of stars} in substructures in the stellar halo as a function of surface brightness are indicators of the importance of recent merging and of the luminosity function of infalling dwarfs. The {\it morphology} of features serves as a guide to the orbital distribution of those dwarfs. Constraints on the earlier merger history (> 8 Gyears ago) can be gleaned from the abundance patterns in halo stars: within our models, dramatic differences in the dominant epoch of accretion or luminosity function of progenitor objects leave clear signatures in the [alpha/Fe] and [Fe/H] distributions of the stellar halo - halos dominated by very early accretion have higher average [alpha/Fe], while those dominated by high luminosity satellites have higher [Fe/H]. This intuition can be applied to reconstruct much about the merger histories of nearby galaxies from current and future data sets.

We present an updated mass model for M31 that makes use of a Spitzer 3.6 $\mu$m image, a mass-to-light ratio gradient based on the galaxy's B-R colour profile, and observed rotation curve data from a variety of sources. We examine cases where the dark matter follows a pure NFW profile and where an initial NFW halo contracts adiabatically in response to the formation of the galaxy. We find that both of these scenarios can produce a reasonable fit to the observed rotation curve data. However, a pure NFW model requires a concentration c_{vir}=51 that is well outside the range predicted in LCDM cosmology and is therefore disfavoured. An adiabatically contracted NFW halo favors an initial concentration c_{vir}=20 and virial mass 8.2x10^{11} M_{\odot}, and this is in line with the cosmological expectations for a galaxy of the size of M31. The best-fit mass is consistent with published estimates from Andromeda Stream kinematics, satellite galaxy radial velocities, and planetary nebulae studies. Finally, using the known linear correlation between rotation curve shear and spiral arm pitch angle, we show that the stellar spiral arm pitch angle of M31 (which cannot be deduced from imaging data due to the galaxy's inclination) is P=24.7\pm4.4 degrees.

We perform a set of high-resolution, dissipationless N-body simulations to investigate the influence of cold dark matter (CDM) substructure on the dynamical evolution of thin galactic disks. Our method combines cosmological simulations of galaxy-sized CDM halos to derive the properties of substructure populations and controlled numerical experiments of consecutive subhalo impacts onto initially-thin, fully-formed disk galaxies. We demonstrate that close encounters between massive subhalos and galactic disks since z~1 should be common occurrences in LCDM models. In contrast, extremely few satellites in present-day CDM halos are likely to have a significant impact on the disk structure. One typical host halo merger history is used to seed controlled N-body experiments of subhalo-disk encounters. As a result of these accretion events, the disk thickens considerably at all radii with the disk scale height increasing in excess of a factor of 2 in the solar neighborhood. We show that interactions with the subhalo population produce a wealth of distinctive morphological signatures in the disk stars including: conspicuous flares; bars; low-lived, ring-like features in the outskirts; and low-density, filamentary structures above the disk plane. We compare a resulting dynamically-cold, ring-like feature in our simulations to the Monoceros ring stellar structure in the MW. The comparison shows quantitative agreement in both spatial distribution and kinematics, suggesting that such observed complex stellar components may arise naturally as disk stars are excited by encounters with subhalos. These findings highlight the significant role of CDM substructure in setting the structure of disk galaxies and driving galaxy evolution.

(Abridged) We conduct a series of high-resolution, dissipationless N-body simulations to investigate the cumulative effect of substructure mergers onto thin disk galaxies in the context of the LCDM paradigm of structure formation. Our simulation campaign is based on a hybrid approach. Substructure properties are culled directly from cosmological simulations of galaxy-sized cold dark matter (CDM) halos. In contrast to what can be inferred from statistics of the present-day substructure populations, accretions of massive subhalos onto the central regions of host halos, where the galactic disk resides, since z~1 should be common occurrences. One host halo merger history is subsequently used to seed controlled numerical experiments of repeated satellite impacts on an initially-thin Milky Way-type disk galaxy. We show that these accretion events produce several distinctive observational signatures in the stellar disk including: a ring-like feature in the outskirts; a significant flare; a central bar; and faint filamentary structures that (spuriously) resemble tidal streams. The final distribution of disk stars exhibits a complex vertical structure that is well-described by a standard ``thin-thick'' disk decomposition. We conclude that satellite-disk encounters of the kind expected in LCDM models can induce morphological features in galactic disks that are similar to those being discovered in the Milky Way, M31, and in other disk galaxies. These results highlight the significant role of CDM substructure in setting the structure of disk galaxies and driving galaxy evolution. Upcoming galactic structure surveys and astrometric satellites may be able to distinguish between competing cosmological models by testing whether the detailed structure of galactic disks is as excited as predicted by the CDM paradigm.

(abridged) The Observatory for Multi-Epoch Gravitational Lens Astrophysics (OMEGA) is a mission concept for a 1.5-m near-UV through near-IR space observatory that will be dedicated to frequent imaging and spectroscopic monitoring of ~100 multiply-imaged active galactic nuclei over the whole sky. Using wavelength-tailored dichroics with extremely high transmittance, efficient imaging in six channels will be done simultaneously during each visit to each target. The separate spectroscopic mode, engaged through a flip-in mirror, uses an image slicer spectrograph. After a period of many visits to all targets, the resulting multidimensional movies can then be analyzed to a) measure the mass function of dark matter substructure; b) measure precise masses of the accreting black holes as well as the structure of their accretion disks and their environments over several decades of physical scale; and c) measure a combination of the local Hubble expansion and cosmological distances to unprecedented precision.

We study the formation of fifty-three galaxy cluster-size dark matter halos formed within a pair of cosmological LCDM N-body simulations, and track the accretion histories of cluster subhalos with masses large enough to host 0.1L* galaxies. By associating subhalos with cluster galaxies, we find the majority of galaxies in clusters experience no pre-processing in the group environment prior to their accretion into the cluster. On average, ~70% of cluster galaxies fall into the cluster potential directly from the field, with no luminous companions in their host halos at the time of accretion; and less than ~12% are accreted as members of groups with five or more galaxies. Moreover, we find that cluster galaxies are significantly less likely to have experienced a merger in the recent past (~6 Gyr) than a field halo of the same mass. These results suggest that local, cluster processes like ram-pressure stripping, galaxy harassment, or strangulation play the dominant role in explaining the difference between cluster and field populations at a fixed stellar mass; and that pre-evolution or past merging in the group environment is of secondary importance for setting cluster galaxy properties for most clusters. The accretion times for z = 0 cluster members are quite extended, with ~20% incorporated into the cluster halo more than 7 Gyr ago and ~20% within the last 2 Gyr. By comparing the observed morphological fractions in cluster and field populations, we estimate an approximate timescale for late-type to early-type transformation within the cluster environment to be ~6 Gyr.

We study the merger histories of galaxy dark matter halos using a high resolution LCDM N-body simulation. Our merger trees follow ~17,000 halos with masses M_0 = (10^11--10^13) Msun at z=0 and track accretion events involving objects as small as m = 10^10 Msun. We find that mass assembly is remarkably self-similar in m/M_0, and dominated by mergers that are ~10% of the final halo mass. While very large mergers, m > 0.4 M_0, are quite rare, sizeable accretion events, m ~ 0.1 M_0, are common. Over the last 10 Gyr, an overwhelming majority (~95%) of Milky Way-sized halos with M_0 = 10^12 Msun have accreted at least one object with greater total mass than the Milky Way disk (m > 5x10^10 Msun), and approximately 70% have accreted an object with more than twice that mass (m > 10^11 Msun). Our results raise serious concerns about the survival of thin-disk dominated galaxies within the current paradigm for galaxy formation in a CDM universe. In order to achieve a ~70% disk-dominated fraction in Milky Way-sized CDM halos, mergers involving m ~ 2x10^11 Msun objects must not destroy disks. Considering that most thick disks and bulges contain old stellar populations, the situation is even more restrictive: these mergers must not heat disks or drive gas into their centers to create young bulges.

We report the discovery of a very bright (m_R = 22.2) Lyman break galaxy at z = 3.03 that appears to be a massive system in a late stage of merging. Deep imaging reveals multiple peaks in the brightness profile with angular separations of ~0.''8 (~25 h^-1 kpc comoving). In addition, high signal-to-noise ratio rest-frame UV spectroscopy shows evidence for ~5 components based on stellar photospheric and ISM absorption lines with a velocity dispersion of sigma ~460 km s^-1 for the three strongest components. Both the dynamics and high luminosity, as well as our analysis of a LCDM numerical simulation, suggest a very massive system with halo mass M ~ 10^13 M_solar. The simulation finds that all halos at z = 3 of this mass contain sub-halos in agreement with the properties of these observed components and that such systems typically evolve into M ~ 10^14 M_solar halos in groups and clusters by z = 0. This discovery provides a rare opportunity to study the properties and individual components of z ~ 3 systems that are likely to be the progenitors to brightest cluster galaxies.

We have used data from the Sloan Digital Sky Survey (SDSS) Data Release 5 to explore the overall structure and substructure of the stellar halo of the Milky Way using about 4 million color-selected main sequence turn-off stars. We fit oblate and triaxial broken power-law models to the data, and found a `best-fit' oblateness of the stellar halo 0.5<c/a<0.8, and halo stellar masses between Galactocentric radii of 1 and 40kpc of (3.7+/-1.2)x10^8 M_sun. The density profile of the stellar halo is approximately r^{-3}; it is possible that the power law slope is shallower inside 20kpc and steeper outside that radius. Yet, we found that all smooth and symmetric models were very poor fits to the distribution of stellar halo stars because the data exhibit a great deal of spatial substructure. We quantified deviations from a smooth oblate/triaxial model using the RMS of the data around the model profile on scales >~100pc, after accounting for the (known) contribution of Poisson uncertainties. The fractional RMS deviation of the actual stellar distribution from any smooth, parameterized halo model is >~40%: hence, the stellar halo is highly structured. We compared the observations with simulations of galactic stellar halos formed entirely from the accretion of satellites in a cosmological context by analysing the simulations in the same way as the data. While the masses, overall profiles, and degree of substructure in the simulated stellar halos show considerable scatter, the properties and degree of substructure in the Milky Way's halo match well the properties of a `typical' stellar halo built exclusively out of the debris from disrupted satellite galaxies. Our results therefore point towards a picture in which an important fraction of the Milky Way's stellar halo has been accreted from satellite galaxies.

Observations of local X-ray absorbers, high-velocity clouds, and distant quasar absorption line systems suggest that a significant fraction of baryons may reside in multi-phase, low-density, extended, ~100 kpc, gaseous halos around normal galaxies. We present a pair of high-resolution SPH (smoothed particle hydrodynamics) simulations that explore the nature of cool gas infall into galaxies, and the physical conditions necessary to support the type of gaseous halos that seem to be required by observations. The two simulations are identical other than their initial gas density distributions: one is initialized with a standard hot gas halo that traces the cuspy profile of the dark matter, and the other is initialized with a cored hot halo with a high central entropy, as might be expected in models with early pre-heating feedback. Galaxy formation proceeds in dramatically different fashions in these two cases. While the standard cuspy halo cools rapidly, primarily from the central region, the cored halo is quasi-stable for ~4 Gyr and eventually cools via the fragmentation and infall of clouds from ~100 kpc distances. After 10 Gyr of cooling, the standard halo's X-ray luminosity is ~100 times current limits and the resultant disk galaxy is twice as massive as the Milky Way. In contrast, the cored halo has an X-ray luminosity that is in line with observations, an extended cloud population reminiscent of the high-velocity cloud population of the Milky Way, and a disk galaxy with half the mass and ~50% more specific angular momentum than the disk formed in the low-entropy simulation. These results suggest that the distribution and character of halo gas provides an important testing ground for galaxy formation models and may be used to constrain the physics of galaxy formation.

Numerical simulations of Milky-Way size Cold Dark Matter (CDM) halos predict a steeply rising mass function of small dark matter subhalos and a substructure count that greatly outnumbers the observed satellites of the Milky Way. Several proposed explanations exist, but detailed comparison between theory and observation in terms of the maximum circular velocity (Vmax) of the subhalos is hampered by the fact that Vmax for satellite halos is poorly constrained. We present comprehensive mass models for the well-known Milky Way dwarf satellites, and derive likelihood functions to show that their masses within 0.6 kpc (M_0.6) are strongly constrained by the present data. We show that the M_0.6 mass function of luminous satellite halos is flat between ~ 10^7 and 10^8 M_\odot. We use the ``Via Lactea'' N-body simulation to show that the M_0.6 mass function of CDM subhalos is steeply rising over this range. We rule out the hypothesis that the 11 well-known satellites of the Milky Way are hosted by the 11 most massive subhalos. We show that models where the brightest satellites correspond to the earliest forming subhalos or the most massive accreted objects both reproduce the observed mass function. A similar analysis with the newly-discovered dwarf satellites will further test these scenarios and provide powerful constraints on the CDM small-scale power spectrum and warm dark matter models.

We use a series of N-body/smoothed particle hydrodynamics simulations and analytic arguments to show that the presence of an effective temperature floor in the interstellar medium at T_F ~ 10^4 K naturally explains the tendency for low-mass galaxies to be more spheroidal, more gas rich, and less efficient in converting baryons into stars than larger galaxies. The trend arises because gas pressure support becomes important compared to angular momentum support in small dark matter haloes. We suggest that dwarf galaxies with rotational velocities ~ 40 km/s do not originate as thin discs, but rather are born as thick, puffy systems. If accreted on to larger haloes, tenuous dwarfs of this kind will be more susceptible to gas loss or tidal transformation than scaled-down versions of larger spirals. For a constant temperature floor, pressure support becomes less important in large haloes, and this produces a tendency for massive isolated galaxies to have thinner discs and more efficient star formation than their less massive counterparts, as observed.

In the course of our survey of the outer halo of the Andromeda Galaxy we have discovered a remote, possible satellite of that system at a projected 162 kpc (11.7 degrees) radius. The fairly elongated (0.31 +/- 0.09 ellipticity) dwarf can be fit with a King profile of 1.07 kpc (d/784 kpc) limiting radius, where the satellite distance, d, is estimated at ~630-850 kpc from the tip of the red giant branch. The newfound galaxy, ``Andromeda XIV'' (``AndXIV''), distinguishes itself from other Local Group galaxies by its extreme dynamics: Keck/DEIMOS spectroscopy reveals it to have a large heliocentric radial velocity (-481 km/s), or -206 km/s velocity relative to M31. Even at its projected radius AndXIV already is at the M31 escape velocity based on the latest M31 mass models. If AndXIV is bound to M31, then recent models with reduced M31 virial masses need revision upward. If not bound to M31, then AndXIV is just now falling into the Local Group for the first time and represents a dwarf galaxy that formed and spent almost its entire life in isolation.

Recent observations have revealed that high surface-brightness, metal-rich debris is present over large regions of the Andromeda (M31) stellar halo. We use a set of numerical models to determine whether extended metal-rich debris is expected to exist in galaxy halos formed in a hierarchical LambdaCDM universe. We identify tidal debris in the simulations according to the current best surface brightness detection limits in M31 and demonstrate that bright features in the outer halo regions usually originate from single satellites, with very little contamination from other sources due to the low probability of tidal streams from two overlapping accretion events. In contrast, high-surface brightness features in the inner halo often originate from multiple progenitors. We also compare the age and metallicity distribution of the debris with the well-mixed stellar halos in which they reside. We find that high-surface brightness tidal debris is produced almost exclusively by relatively high mass progenitors (Mstar ~ 10^7- 10^9 Msun) and thus is expected to be of moderate- to high-metallicity. Similarly, in our models the smooth inner halo is expected to be metal-rich as this region has been built up mainly from massive satellites. Our results imply that the stellar populations of substructure observed around external galaxies with current techniques should not resemble old and metal-poor dwarf spheroidal satellites, nor the underlying component of the stellar halo.

We use kinematic data from three new, nearby, extremely low-luminosity Milky Way dwarf galaxies (Ursa Major II, Willman 1, and Coma Berenices) to constrain the properties of their dark matter halos, and from these make predictions for the gamma-ray flux from annihilation of dark matter particles in these halos. We show that these 10^3 solar luminosity dwarfs are the most dark matter dominated galaxies in the Universe, with total masses within 100 pc in excess of 10^6 solar masses. Coupled with their relative proximity, their large masses imply that they should have mean gamma-ray fluxes comparable to or greater than any other known satellite galaxy of the Milky Way. Our results are robust to both variations of the inner slope of the density profile and the effect of tidal interactions. The fluxes could be boosted by up to two orders of magnitude if we include the density enhancements caused by surviving dark matter substructure.

Galaxy pairs provide a potentially powerful means of studying triggered star formation from galaxy interactions. We use a large cosmological N-body simulation coupled with a well-tested semi-analytic substructure model to demonstrate that the majority of galaxies in close pairs reside within cluster or group-size halos and therefore represent a biased population, poorly suited for direct comparison to ``field'' galaxies. Thus, the frequent observation that some types of galaxies in pairs have redder colors than ``field'' galaxies is primarily a selection effect. We select galaxy pairs that are isolated in their dark matter halos with respect to other massive subhalos (N=2 halos) and a control sample of isolated galaxies (N=1 halos) for comparison. We then apply these selection criteria to a volume-limited subset of the 2dF Galaxy Redshift Survey with M_Bj <= -19 and obtain the first clean measure of the typical fraction of galaxies affected by triggered star formation and the average elevation in the star formation rate. We find that 24% (30.5%) of these L^\star and sub-L^{\star} galaxies in isolated 50 (30) kpc/h pairs exhibit star formation that is boosted by a factor of >~ 5 above their average past value, while only 10% of isolated galaxies in the control sample show this level of enhancement. Thus, 14% (20 %) of the galaxies in these close pairs show clear triggered star formation. The isolation criteria we develop provide a means to constrain star formation and feedback prescriptions in hydrodynamic simulations and a very general method of understanding the importance of triggered star formation in a cosmological context. (Abridged.)

We use HST/ACS to study the resolved stellar populations of the nearby, nearly edge-on galaxy NGC4244 across its outer disk surface density break. The stellar photometry allows us to study the distribution of different stellar populations and reach very low equivalent surface brightnesses. We find that the break occurs at the same radius for young, intermediate age, and old stars. The stellar density beyond the break drops sharply by a factor of at least 600 in 5 kpc. The break occurs at the same radius independent of height above the disk, but is sharpest in the midplane and nearly disappears at large heights. These results make it unlikely that truncations are caused by a star formation threshold alone: the threshold would have to keep the same radial position from less than 100 Myr to 10 Gyr ago, in spite of potential disturbances such as infall and redistribution of gas by internal processes. A dynamical interpretation of truncation formation is more likely such as due to angular momentum redistribution by bars or density waves, or heating and stripping of stars caused by the bombardment of dark matter sub-halos. The latter explanation is also in quantitative agreement with the small diffuse component we see around the galaxy.

We present radial mass profiles within 0.3 r_vir for 16 relaxed galaxy groups-poor clusters (kT range 1-3 keV) selected for optimal mass constraints from the Chandra and XMM data archives. After accounting for the mass of hot gas, the resulting mass profiles are described well by a two-component model consisting of dark matter (DM), represented by an NFW model, and stars from the central galaxy. The stellar component is required only for 8 systems, for which reasonable stellar mass-to-light ratios (M/L_K) are obtained, assuming a Kroupa IMF. Modifying the NFW dark matter halo by adiabatic contraction does not improve the fit and yields systematically lower M/L_K. In contrast to previous results for massive clusters, we find that the NFW concentration parameter (c_vir) for groups decreases with increasing M_vir and is inconsistent with no variation at the 3 sigma level. The normalization and slope of the c_vir-M_vir relation are consistent with the standard LambdaCDM cosmological model with sigma_8 = 0.9. The small intrinsic scatter measured about the c_vir-M_vir relation implies the groups represent preferentially relaxed, early forming systems. The mean gas fraction (f =0.05 +/- 0.01) of the groups measured within an overdensity Delta=2500 is lower than for hot, massive clusters, but the fractional scatter (sigma_f/f=0.2) for groups is larger, implying a greater impact of feedback processes on groups, as expected.

We make predictions for diffuse stellar mass fractions in dark matter halos from the scales of small spiral galaxies to those of large galaxy clusters. We use an extensively-tested analytic model for subhalo infall and evolution and empirical constraints from galaxy survey data to set the stellar mass in each accreted subhalo to model diffuse light. We add stellar mass to the diffuse light as subhalos become disrupted due to interactions within their host halos. We predict that the stellar mass fraction in diffuse, intrahalo light should rise on average from ~0.5% to approximately 20% from small galaxy halos to poor groups. The trend with mass flattens considerably beyond the group scale, increasing weakly from a fraction of ~20% in poor galaxy clusters (~10^14 M_sun) to roughly ~30% in massive clusters (~10^15 M_sun). The mass-dependent diffuse light fraction is governed primarily by the empirical fact that the mass-to-light ratio in galaxy halos must vary as a function of halo mass. Galaxy halos have little diffuse light because they accrete most of their mass in small subhalos that themselves have high mass-to-light ratios; stellar halos around galaxies are built primarily from disrupted dwarf-irregular-type galaxies with M*~10^8.5 M_sun. The diffuse light in group and cluster halos is built from satellite galaxies that form stars efficiently and have correspondingly low mass-to-light ratios; intracluster light is dominated by material liberated from massive galaxies with M*~10^11 M_sun. Our results are consistent with existing observations spanning the galaxy, group, and cluster scale; however, they can be tested more rigorously in future deep surveys for faint diffuse light.

We present the concentration (c)-virial mass (M) relation of 39 galaxy systems ranging in mass from individual early-type galaxies up to the most massive galaxy clusters, (0.06-20) x 10^{14} M_sun. We selected for analysis the most relaxed systems possessing the highest quality data currently available in the Chandra and XMM public data archives. A power-law model fitted to the X-ray c-M relation requires at high significance (6.6 sigma) that c decreases with increasing M, which is a general feature of CDM models. The median and scatter of the c-M relation produced by the flat, concordance LCDM model (Omega_m=0.3, sigma_8=0.9) agrees with the X-ray data provided the sample is comprised of the most relaxed, early forming systems, which is consistent with our selection criteria. Holding the rest of the cosmological parameters fixed to those in the concordance model the c-M relation requires 0.76< sigma_8 <1.07 (99% conf.), assuming a 10% upward bias in the concentrations for early forming systems. The tilted, low-sigma_8 model suggested by a new WMAP analysis is rejected at 99.99% confidence, but a model with the same tilt and normalization can be reconciled with the X-ray data by increasing the dark energy equation of state parameter to w ~ -0.8. When imposing the additional constraint of the tight relation between sigma_8 and Omega_m from studies of cluster abundances, the X-ray c-M relation excludes (>99% conf.) both open CDM models and flat CDM models with Omega_m ~1. This result provides novel evidence for a flat, low-Omega_m universe with dark energy using observations only in the local (z << 1) universe. Possible systematic errors in the X-ray mass measurements of a magnitude ~10% suggested by CDM simulations do not change our conclusions.

We examine the prospects for detecting gamma-rays from dark matter annihilation in the six most promising dwarf spheroidal (dSph) satellite galaxies of the Milky Way. We use recently-measured velocity dispersion profiles to provide a systematic investigation of the dark matter mass distribution of each galaxy, and show that the uncertainty in the gamma-ray flux from mass modeling is less than a factor of ~ 5 for each dSph if we assume a smooth NFW profile. We show that Ursa Minor and Draco are the most promising dSphs for gamma-ray detection with GLAST and other planned observatories. For each dSph, we investigate the flux enhancement resulting from halo substructure, and show that the enhancement factor relative to a smooth halo flux cannot be greater than about 100. This enhancement depends very weakly on the lower mass cut-off scale of the substructure mass function. While the amplitude of the expected flux from each dSph depends sensitively on the dark matter model, we show that the flux ratios between the six Sphs are known to within a factor of about 10. The flux ratios are also relatively insensitive to the current theoretical range of cold dark matter halo central slopes and substructure fractions.

(Abridged) We present a semi-analytic model for Cold Dark Matter halo substructure that can be used as a framework for studying galaxy formation and an ingredient in halo models of galaxy clustering. We perform a comprehensive comparison of the model to the results of a suite of high-resolution cosmological simulations. The comparisons reveal that subhalo statistics, such as velocity and mass functions, radial distributions, and the halo occupation distributions, agree well over three orders of magnitude in host halo mass at various redshifts. Both in the simulations and in our model, the radial distributions of subhalos are significantly shallower than that of the dark matter density. The abundance of subhalos in a host is set by competition between destruction and new accretion. High mass halos and halos at high redshift tend to host more subhalos because the subhalos have, on average, been accreted more recently. Similarly, at fixed mass, halos that formed more recently host more subhalos. Observed "fossil groups" may represent an extreme tail of this correlation. We find a related correlation between host halo concentration and subhalo abundance at fixed host mass, Nsat ~ c^-a, where a changes with redshift and host-to-subhalo mass ratio. Lastly, we use our model to populate host halos in one of our high-resolution cosmological simulations, replacing the subhalos resolved in the simulation with subhalos computed according to our model. We show that the resulting 2pt correlation function of such a "hybrid" halo ensemble is indistinguishable from that measured directly in the simulation. This supports one of the key tenets of the standard halo model -- the assumption that the halo occupation distribution is statistically independent of host halo environment.

We show that dark matter emerging from late decays (z < 1000) produces a linear power spectrum identical to that of Cold Dark Matter (CDM) on all observationally relavant scales (> 0.1 Mpc), and simultaneously generates observable constant-density cores in small dark matter halos. We refer to this class of models as meta-Cold Dark Matter (mCDM), because it is born with non-relativistic velocities from the decays of cold thermal relics. The constant-density cores are a result of the low phase-space density of mCDM at birth. Warm dark matter cannot produce similar size phase-space limited cores without saturating the Ly-alpha power spectrum bounds. Dark matter dominated galaxy rotation curves and stellar velocity dispersion profiles may provide the best means to discriminate between mCDM and CDM. mCDM candidates are motivated by the particle spectrum of supersymmetric and extra dimensional extensions to the standard model of particle physics.

We show initial results from our ongoing HST GHOSTS survey of the resolved stellar envelopes of 14 nearby, massive disk galaxies. In hierarchical galaxy formation the stellar halos and thick disks of galaxies are formed by accretion of minor satellites and therefore contain valuable information about the (early) assembly process of galaxies. We detect for the first time the very small halo of NGC4244, a low mass edge-on galaxy. We find that massive galaxies have very extended halos, with equivalent surface brightnesses of 28-29 V-mag/arcsec^2 at 20-30 kpc from the disk. The old RGB stars of the thick disk in the NGC891 and NGC4244 edge-on galaxies truncate at the same radius as the young thin disk stars, providing insights into the formation of both disk truncations and thick disks. We furthermore present the stellar populations of a very low surface brightness stream around M83, the first such a stream resolved into stars beyond those of the Milky Way and M31.

We show that measurements of stellar proper motions in dwarf spheroidal galaxies provide a powerful probe of the nature of dark matter. Allowing for general dark matter density profiles and stellar velocity anisotropy profiles, we show that the log-slope of the dark matter profile at about twice the stellar core (King) radius can be measured to within \pm 0.2 when the proper motions of 200 stars are added to standard line-of-sight velocity dispersion data. This measurement of the log-slope provides a test of Cold and Warm Dark Matter theories at a sensitivity not possible with line-of-sight velocity dispersion measurements alone. The upcoming SIM PlanetQuest will have the sensitivity to obtain the required number of proper motions in Milky Way dwarf spheroidal galaxies.

Standard LCDM predicts that the major merger rate of galaxy-size dark matter halos rises rapidly with redshift. The average number of close companions per galaxy, Nc, is often used to infer the galaxy merger rate, however, recent observational studies suggest that Nc evolves very little with redshift. Here we use a "hybrid" N- body simulation plus analytic substructure model to predict Nc directly. We identify dark matter subhalos with galaxies and show that the observed lack of close pair count evolution arises because the high merger rate per halo at early times is counteracted by a decrease in the number of halos massive enough to host a galaxy pair. We compare our results to data compiled from the DEEP2, SSRS2, and the UZC redshift surveys. Observed pair counts match our predictions if we assume a monotonic mapping between galaxy luminosity and the maximum circular velocity that each subhalo had when it was first accreted onto its host halo. This suggests that satellite galaxies are significantly more resilient to destruction than are dissipationless dark matter subhalos. We argue that while Nc does not provide a direct measure of the halo merger rate, it offers a powerful means to constrain the Halo Occupation Distribution and the spatial distribution of galaxies within halos. Interpreted in this way, close pair counts provide a useful test of galaxy formation processes on < 100 kpc scales.

We investigate the dependence of dark matter halo clustering on halo formation time, density profile concentration, and subhalo occupation number, using high-resolution numerical simulations of a LCDM cosmology. We confirm results that halo clustering is a function of halo formation time, and that this trend depends on halo mass. For the first time, we show unequivocally that halo clustering is a function of halo concentration and show that the dependence of halo bias on concentration, mass, and redshift can be accurately parameterized in a simple way: b(c,M|z) = b(M|z) b(c|M/M*). The scaling between bias and concentration changes sign with the value of M/M*: high concentration (early forming) objects cluster more strongly for M <~ M* while low concentration (late forming) objects cluster more strongly for rare high-mass halos, M >~ M*. We show the first explicit demonstration that host dark halo clustering depends on the halo occupation number (of dark matter subhalos) and discuss implications for halo model calculations of dark matter power spectra and galaxy clustering statistics. The effect of these halo properties on clustering is strongest for early-forming dwarf-mass halos, which are significantly more clustered than typical halos of their mass. Our results suggest that isolated low-mass galaxies (e.g. low surface-brightness dwarfs) should have more slowly-rising rotation curves than their clustered counterparts, and may have consequences for the dearth of dwarf galaxies in voids. They also imply that self calibrating richness-selected cluster samples with their clustering properties might overestimate cluster masses and bias cosmological parameter estimation.

We present an X-ray analysis of the radial mass profile of the radio-quiet galaxy cluster A2589 between 0.015-0.25 r_vir using an XMM-Newton observation. Except for a ~16 kpc shift of the X-ray center of the R=45-60 kpc annulus, A2589 possesses a remarkably symmetrical X-ray image and is therefore an exceptional candidate for precision studies of its mass profile by applying hydrostatic equilibrium. The total gravitating matter profile is well described by the NFW model (fractional residuals <~10%) with c_vir=6.1 +/- 0.3 and M_vir = 3.3 +/- 0.3 x 10^{14} M_sun (r_vir = 1.74 +/- 0.05 Mpc) in excellent agreement with LCDM. When the mass of the hot ICM is subtracted from the gravitating matter profile, the NFW model fitted to the resulting dark matter (DM) profile produces essentially the same result. However, if a component accounting for the stellar mass (M_*) of the cD galaxy is included, then the NFW fit to the DM profile is substantially degraded in the central r ~50 kpc for reasonable M_*/L_V. Modifying the NFW DM halo by adiabatic contraction arising from the early condensation of stellar baryons in the cD galaxy further degrades the fit. The fit is improved substantially with a Sersic-like model recently suggested by high resolution N-body simulations but with an inverse Sersic index, alpha ~0.5, a factor of ~3 higher than predicted. We argue that neither random turbulent motions nor magnetic fields can provide sufficient non-thermal pressure support to reconcile the XMM mass profile with adiabatic contraction of a CDM halo assuming reasonable M_*/L_V. Our results support the scenario where, at least for galaxy clusters, processes during halo formation counteract adiabatic contraction so that the total gravitating mass in the core approximately follows the NFW profile.

We present a Chandra study of mass profiles in 7 elliptical galaxies, of which 3 have galaxy-scale and 4 group-scale halos demarcated at 1E13Msun. These represent the best available data for nearby objects with comparable X-ray luminosities. We measure ~flat mass-to-light (M/L) profiles within an optical half-light radius (Reff), rising by an order of magnitude at ~10Reff, which confirms the presence of dark matter (DM). The data indicate hydrostatic equilibrium, which is also supported by agreement with studies of stellar kinematics in elliptical galaxies. The data are well-fitted by a model comprising an NFW DM profile and a baryonic component following the optical light. The distribution of DM halo concentration parameters (c) versus Mvir agrees with LCDM predictions and our observations of bright groups. Concentrations are slightly higher than expected, which is most likely a selection effect. Omitting the stellar mass drastically increases c, possibly explaining large concentrations found by some past observers. The stellar M/LK agree with population synthesis models, assuming a Kroupa IMF. Allowing adiabatic compression (AC) of the DM halo by baryons made M/L more discrepant, casting some doubt on AC. Our best-fitting models imply total baryon fractions \~0.04--0.09, consistent with models of galaxy formation incorporating strong feedback. The groups exhibit positive temperature gradients, consistent with the "Universal" profiles found in other groups and clusters, whereas the galaxies have negative gradients, suggesting a change in the evolutionary history of the systems around Mvir=1E13 Msun.

Large galaxies are thought to form hierarchically, from the accretion and disruption of many smaller galaxies. Such a scenario should naturally lead to galactic phase-space distributions containing some degree of substructure. We examine the errors in mass estimates of galaxies and their dark halos made using the projected phase-space distribution of a tracer population (such as a globular cluster system or planetary nebulae) due to falsely assuming that the tracers are distributed randomly. The level of this uncertainty is assessed by applying a standard mass estimator to samples drawn from 11 random realizations of galaxy halos containing levels of substructure consistent with current models of structure formation. We find that substructure will distort our mass estimates by up to ~20% - a negligible error compared to statistical and measurement errors in current derivations of masses for our own and other galaxies. However, this represents a fundamental limit to the accuracy of any future mass estimates made under the assumption that the tracer population is distributed randomly, regardless of the size of the sample or the accuracy of the measurements.

We use measurements of the stellar velocity dispersion profile of the Fornax dwarf spheroidal galaxy to derive constraints on its dark matter distribution. Though the data are unable to distinguish between models with small cores and those with cusps, we show that a large > 1 kpc dark matter core in Fornax is highly implausible. Irrespective of the origin of the core, reasonable dynamical limits on the mass of the Fornax halo constrain its core radius to be no larger than 700 pc. We derive an upper limit core radius of 300 pc by demanding that the central phase space density of Fornax not exceed that directly inferred from the rotation curves of low-mass spiral galaxies. Further, if the halo is composed of warm dark matter then phase-space constraints force the core to be quite small in order to avoid conservative limits from the Ly alpha forest power spectrum, implying a core radius < 85 pc. We discuss our results in the context of the idea that the extended globular cluster distribution in Fornax can be explained by the presence of a large 1.5 kpc core. A self-consistent core of this size would be drastically inconsistent with the expectations of standard warm or cold dark matter models, and would also require an unreasonably massive dark matter halo, with a maximum circular velocity of 200 km/s.

(Abridged) The violent hierarchical nature of the LCDM cosmology poses serious difficulties for the formation of disk galaxies. To help resolve these issues, we describe a new, merger-driven scenario for the cosmological formation of disk galaxies at high redshifts that supplements the standard model based on dissipational collapse.In this picture, large gaseous disks may be produced from high-angular momentum mergers of systems that are gas-dominated, i.e. M_gas/(M_gas +M_star > 0.5 at the height of the merger. Pressurization from the multiphase structure of the interstellar medium prevents the complete conversion of gas into stars during the merger, and if enough gas remains to form a disk, the remnant eventually resembles a disk galaxy. We perform numerical simulations of galaxy mergers to study how supernovae feedback strength, supermassive black hole growth and feedback, progenitor gas fraction, merger mass-ratio, and orbital geometry impact the formation of remnant disks. We find that disks can build angular momentum through mergers and the degree of rotational support of the baryons in the merger remnant is primarily related to feedback processes associated with star formation. Disk-dominated remnants are restricted to form in mergers that are gas-dominated at the time of final coalescence and gas-dominated mergers typically require extreme progenitor gas fractions (>80%). We also show that the formation of rotationally-supported stellar systems in mergers is not restricted to idealized orbits, or major or minor mergers. We suggest that the hierarchical nature of the LCDM cosmology and the physics of the interstellar gas may act together to form spiral galaxies by building the angular momentum of disks through early, gas-dominated mergers.

We investigate the use of spiral arm pitch angles as a probe of disk galaxy mass profiles. We confirm our previous result that spiral arm pitch angles (P) are well-correlated with the rate of shear (S) in disk galaxy rotation curves, by using a much larger sample (51 galaxies) than used previously (17 galaxies). We use this correlation to argue that imaging data alone can provide a powerful probe of galactic mass distributions out to large lookback times. In contrast to previous work, we show that observed spiral arm pitch angles are similar when measured in the optical (at 0.4 um) and the near-infrared (at 2.1 um) with a mean difference of 2.3+/-2.7 degrees. This is then used to strengthen the known correlation between P and S using B-band images. We then use two example galaxies to demonstrate how an inferred shear rate coupled with a bulge-disk decomposition model and a Tully-Fisher derived velocity normalization can be used to place constraints on a galaxy's baryon fraction and dark matter halo profile. We show that ESO 582-G12, a galaxy with a high shear rate (slightly declining rotation curve) at ~10 kpc, favors an adiabatically contracted halo, with high initial NFW concentration (c_vir > 16) and a high fraction of halo baryons in the form of stars (~15-40%). In contrast, IC 2522 has a low shear rate (rising rotation curve) at ~10 kpc and favors non-adiabatically contracted models with low NFW concentrations (c_vir ~ 2-8) and a low stellar baryon fraction <10%.

[Abridged] Motivated by upcoming data from astrometric and spectroscopic surveys of the Galaxy, we explore the chemical abundance properties and phase-space distributions in hierarchically-formed stellar halo simulations set in a LambdaCDM Universe. Our sample of Milky-Way type stellar halo simulations result in average metallicities that range from [Fe/H] = -1.3 to -0.9, with the most metal poor halos resulting from accretion histories that lack destructive mergers with massive (metal rich) satellites. Our stellar halo metallicities increase with stellar halo mass. The slope of the [Fe/H]-stellar mass trend mimics that of the satellite galaxies that were destroyed to build the halos, implying that the relation propagates hierarchically. All simulated halos contain a significant fraction of old stellar populations accreted more than 10 Gyr ago and in a few cases, some intermediate age populations exist. In contrast with the Milky Way, many of our simulated stellar halos contain old stellar populations which are metal rich, originating in the early accretion of massive satellites. We suggest that the (metal rich) stellar halo of M31 falls into this category, while the more metal poor halo of the Milky Way is lacking in early massive accretion events. Interestingly, our hierarchically-formed stellar halos often have non-negligible metallicity gradients in both [Fe/H] and [alpha/Fe]. These gradients extend a few tens of kpc, and can be as large as 0.5 dex in [Fe/H] and 0.2 dex in [alpha/Fe]. Finally, we find that chemical abundances can act as a rough substitute for time of accretion of satellite galaxies. We propose a criterion for identifying tidal streams spatially by selecting stars with [alpha/Fe] ratios below solar.

The search for the first illuminated astronomical sources in the universe is at the edge of the cosmic frontier. Promising techniques for discovering the first objects and their effects span the electromagnetic spectrum and include gravitational waves. We summarize a workshop on discovering and understanding these sources which was held in May 2005 through the Center for Cosmology at the University of California, Irvine.

(Abridged) The hierarchical formation scenario for the stellar halo requires the accretion and disruption of dwarf galaxies, yet low-metallicity halo stars are enriched in alpha-elements compared to similar, low-metallicity stars in dwarf spheroidal (dSph) galaxies. We address this primary challenge for the hierarchical formation scenario for the stellar halo by combining chemical evolution modelling with cosmologically-motivated mass accretion histories for the Milky Way dark halo and its satellites. We demonstrate that stellar halo and dwarf galaxy abundance patterns can be explained naturally within the LCDM framework. Our solution relies fundamentally on the LCDM model prediction that the majority of the stars in the stellar halo were formed within a few relatively massive, ~5 x 10^10 Msun, dwarf irregular (dIrr)-size dark matter halos, which were accreted and destroyed ~10 Gyr in the past. These systems necessarily have short-lived, rapid star formation histories, are enriched primarily by Type II supernovae, and host stars with enhanced [a/Fe] abundances. In contrast, dwarf spheroidal galaxies exist within low-mass dark matter hosts of ~10^9 Msun, where supernovae winds are important in setting the intermediate [a/Fe] ratios observed. Our model includes enrichment from Type Ia and Type II supernovae as well as stellar winds, and includes a physically-motivated supernovae feedback prescription calibrated to reproduce the local dwarf galaxy stellar mass - metallicity relation. We use representative examples of the type of dark matter halos we expect to host a destroyed ``stellar halo progenitor'' dwarf, a surviving dIrr, and a surviving dSph galaxy, and show that their derived abundance patterns, stellar masses, and gas masses are consistent with those observed for each type of system.

Cosmological simulations of galaxy formation predict a universal form for the mass profile of dark matter (DM) halos from cluster to galaxy scales. Remarkably few interesting constraints exist, however, on DM halos in early-type galaxies. Using Chandra we present the temperature, density and mass profiles of a small sample of early-type galaxies, revealing significant DM in each case. When a component is included to account for stellar mass and the DM halo is allowed to respond adiabatically to the baryonic condensation into stars, the mass profiles are well-fitted by the universal profile, with Virial masses and concentrations in agreement with simulations. However, only ~half, or less, of the mass within R_e seems attributable to the stars, implying stellar M_*/L_B ~1-5.

The Chandra and XMM data for a sample of 19 relaxed groups/poor clusters, covering the temperature range 1-3 keV and selected to have the best available data, reveal a remarkable similarity in their temperature profile: cool cores outside of which the temperatures reach a peak for radii less than 0.1 of the virial radius and then decline. We fitted the derived mass profiles using an NFW model, which provides a good fit to the data when accounting for the central galaxy in the inner region. The concentration parameters and virial masses are in the range c = 5-22 and M = 2 x 10^13 - 4 x 10^14 solar masses, in general agreement with the concentrations found in numerical simulations.

We have analyzed an XMM-Newton observation of the cluster Abell 2589. Apart from a low-level asymmetry in the central region, the cluster appears very relaxed and does not show presence of a central AGN. We derived constraints for the radial temperature, density and, assuming hydrostatic equilibrium, mass profiles. We find that the best fit to the dark matter profile is given by the Sersic-like profile proposed by Navarro et al. (2004). The NFW model does not provide a good fit. We also tested whether the central stellar component could affect the profile through the adiabatic contraction model but were unable to distinguish it from a simple "dark matter + stars" modeling.

(Abridged) We study predictions for galaxy cluster observables that can test the statistics of dark matter halo shapes expected in a flat LCDM universe. We present a simple analytical model for the prediction of cluster-scale X-ray observations, approximating clusters as isothermal systems in hydrostatic equilibrium, and dark matter haloes as ellipsoids with uniform axial ratios. We test the model against high-resolution, hydrodynamic cluster simulations to gauge its reliability. We find that this simple prescription does a good job of predicting the distribution of cluster X-ray ellipticities compared to the simulations as long as one focuses on cluster regions that are less sensitive to recent mergers. Based on this simple model, the distribution of cluster-size halo shapes expected in the concordance LCDM cosmology implies an X-ray ellipticity distribution with a mean of 0.32 +- 0.01 and a scatter of 0.14 +- 0.01 for the mass range (1-4)x10^{14} Msun/h. We find it important to include the mass dependence of halo shape to make comparisons to observational samples that contain many, very massive clusters. We analyse the systematics of four observational samples of cluster ellipticities and find that our results are statistically compatible with observations. In particular, we find remarkably good agreement between two recent ROSAT samples and LCDM predictions that DO NOT include gas cooling. We also test how well our analytical model can predict Sunyaev-Zel'dovich decrement maps and find that it is less successful although still useful; the model does not perform as well as a function of flux level in this case because of the changing triaxiality of dark matter haloes as a function of radial distance. Both this effect and the changing alignment of isodensity shells of dark matter haloes leave an imprint on cluster gas...

Using six high resolution dissipationless simulations with a varying box size in a flat LCDM universe, we study the mass and redshift dependence of dark matter halo shapes for M_vir = 9.0e11 - 2.0e14, over the redshift range z=0-3, and for two values of sigma_8=0.75 and 0.9. Remarkably, we find that the redshift, mass, and sigma_8 dependence of the mean smallest-to-largest axis ratio of halos is well described by the simple power-law relation <s> = (0.54 +- 0.02)(M_vir/M_*)^(-0.050 +- 0.003), where s is measured at 0.3 R_vir and the z and sigma_8 dependences are governed by the characteristic nonlinear mass, M_*=M_*(z,sigma_8). We find that the scatter about the mean s is well described by a Gaussian with sigma ~ 0.1, for all masses and redshifts. We compare our results to a variety of previous works on halo shapes and find that reported differences between studies are primarily explained by differences in their methodologies. We address the evolutionary aspects of individual halo shapes by following the shapes of the halos through ~100 snapshots in time. We determine the formation scalefactor a_c as defined by Wechsler et al. (2002) and find that it can be related to the halo shape at z = 0 and its evolution over time.

Recently, it has been shown that a correlation exists between the rate of shear and the spiral arm pitch angle in disk galaxies. The rate of shear depends upon the shape of the rotation curve, which is dependent upon the mass distribution in spiral galaxies. Here, we present an imporoved correlation between shear rate and spiral arm pitch angle, by increasing the sample size. We also use an adiabatic infall code to show that the rate of shear is most strongly correlated with the central mass concentration, c_m. The spin parameter, $\lambda$, and the fraction of baryons that cool, F, cause scatter in this correlation. Limiting this scatter, such that it is equal to that in the correlation between shear rate and pitch angle, and using a value of F=0.1 to 0.2, the spin parameter must be in the range 0.03<lambda<0.09 for spiral galaxies. We also derive an equation which links spiral arm pitch angle directly to c_m.

We present a cosmologically motivated model for the hierachical formation of the stellar halo which includes a semi-analytic treatment of galactic chemical enrichment coupled to numerical simulations that track the orbital evolution and tidal disruption of satellites. A major motivating factor in this investigation is the observed systematic difference between the chemical abundances of stars in satellite galaxies and those in the Milky Way halo. Specifically, for the same [Fe/H] values, stars in neighboring satellite galaxies display significantly lower [alpha/Fe] ratios than stars in the halo. We find that the observed chemical abundance patterns are a natural outcome of the process of hierarchical assembly of the Galaxy. This result follows because the stellar halo in this context is built up from satellite galaxies accreted early on (more than 8-9 Gyr ago) and enriched in alpha-elements produced in Type II supernovae. In contrast, satellites which survive today are typically accreted late (within the last 4-5 Gyr) with nearly solar [alpha/Fe] values as a result of contributions from both Type II and Type Ia supernovae. We use our model to investigate the abundance distribution functions (using both [Fe/H] and [alpha/Fe] ratios) for stars in the halo and in surviving satellites. Our results suggest that the shapes and peaks in the abundance distribution functions provide a direct probe of the accretion histories of galaxies.

If the favored hierarchical cosmological model is correct, then the Milky Way system should have accreted ~100-200 luminous satellite galaxies in the past \~12 Gyr. We model this process using a hybrid semi-analytic plus N-body approach which distinguishes explicitly between the evolution of light and dark matter in accreted satellites. This distinction is essential to our ability to produce a realistic stellar halo, with mass and density profile much like that of our own Galaxy, and a surviving satellite population that matches the observed number counts and structural parameter distributions of the satellite galaxies of the Milky Way. Our model stellar halos have density profiles which typically drop off with radius faster than those of the dark matter. They are assembled from the inside out, with the majority of mass (~80%) coming from the \~15 most massive accretion events. The satellites that contribute to the stellar halo have median accretion times of ~9 Gyr in the past, while surviving satellite systems have median accretion times of ~5 Gyr in the past. This implies that stars associated with the inner halo should be quite different chemically from stars in surviving satellites and also from stars in the outer halo or those liberated in recent disruption events. We briefly discuss the expected spatial structure and phase space structure for halos formed in this manner. Searches for this type of structure offer a direct test of whether cosmology is indeed hierarchical on small scales.

We discuss how measurements of the absorption of gamma-rays from GeV to TeV energies via pair production on the extragalactic background light (EBL) can probe important issues in galaxy formation. Semi-analytic models (SAMs) of galaxy formation, based on the flat LCDM hierarchical structure formation scenario, are used to make predictions of the EBL from 0.1 to 1000 microns. SAMs incorporate simplified physical treatments of the key processes of galaxy formation -- including gravitational collapse and merging of dark matter halos, gas cooling and dissipation, star formation, supernova feedback and metal production. We will summarize SAM successes and failures in accounting for observations at low and high redshift. New ground- and space-based gamma ray telescopes will help to determine the EBL, and also help to explain its origin by constraining some of the most uncertain features of galaxy formation theory, including the stellar initial mass function, the history of star formation, and the reprocessing of light by dust. On a separate topic concerning gamma ray cosmology, we discuss a new theoretical insight into the distribution of dark matter at the center of the Milky Way, and its implications concerning the high energy gamma rays observed from the Galactic center.

We investigate the effect of dark energy on the density profiles of dark matter haloes with a suite of cosmological N-body simulations and use our results to test analytic models. We consider constant equation of state models, and allow both w>-1 and w<-1. Using five simulations with w ranging from -1.5 to -0.5, and with more than ~1600 well-resolved haloes each, we show that the halo concentration model of Bullock et al. (2001) accurately predicts the median concentrations of haloes over the range of w, halo masses, and redshifts that we are capable of probing. We find that the Bullock et al. (2001) model works best when halo masses and concentrations are defined relative to an outer radius set by a cosmology-dependent virial overdensity. For a fixed power spectrum normalization and fixed-mass haloes, larger values of w lead to higher concentrations and higher halo central densities, both because collapse occurs earlier and because haloes have higher virial densities. While precise predictions of halo densities are quite sensitive to various uncertainties, we make broad comparisons to galaxy rotation curve data. At fixed power spectrum normalization (fixed sigma_8), w>-1 quintessence models seem to exacerbate the central density problem relative to the standard w=-1 model. Meanwhile w<-1 models help to reduce the apparent discrepancy. We confirm that the Jenkins et al. (2001) halo mass function provides an excellent approximation to the abundance of haloes in our simulations and extend its region of validity to include models with w<-1.

The standard treatment of cooling in Cold Dark Matter halos assumes that all of the gas within a ``cooling radius'' cools and contracts monolithically to fuel galaxy formation. Here we take into account the expectation that the hot gas in galactic halos is thermally unstable and prone to fragmentation during cooling and show that the implications are more far-reaching than previously expected: allowing multi-phase cooling fundamentally alters expectations about gas infall in halos and naturally explains the bright-end cutoff in the galaxy luminosity function. We argue that cooling should proceed via the formation of high-density, 10^4 K clouds, pressure-confined within a hot gas background. The background medium has a low density, and can survive as a stable corona with a long cooling time. The fraction of baryons contained in the residual hot core grows with halo mass because the cooling density increases, and this leads to an upper-mass limit in quiescent, non-merged galaxies of ~10^11 Msun. In this scenario, galaxy formation is fueled by the infall of pressure-supported clouds. For Milky-Way-size systems, clouds of mass ~ 5x10^6 Msun that formed or merged within the last several Gyrs should still exist as a residual population in the halo, with a total mass in clouds of ~ 2 x 10^10 Msun. The mass of the Milky Way galaxy is explained naturally in this model, and is a factor of two smaller than would result in the standard treatment without feedback. We expect clouds in galactic halos to be ~ 1 kpc in size and to extend ~150 kpc from galactic centers. The predicted properties of clouds match well the observed radial velocities, angular sizes, column densities, and velocity widths of High Velocity Clouds around our Galaxy. The clouds also explain high-ion absorption systems at z<1.

(ABRIDGED) We present a model in which the seeds of supermassive black holes form from the lowest angular momentum gas in proto-galaxies at high redshift. We show that this leads to a correlation between black hole masses and spheroid properties, as observed today. We assume that gas in early-forming, rare halos has a distribution of specific angular momentum similar to that expected for dark matter halos. This distribution has a significant low angular momentum tail which implies that every proto-galaxy should contain gas that ends up in a high-density disc. In halos more massive than a critical threshold of \~10^8Msun, the discs are gravitationally unstable, and experience an efficient Lin-Pringle viscosity that allows mass inflow. This process continues until the first massive stars disrupt the disc. Black holes are created with a characteristic mass of ~10^5 Msun, independent of the redshift of formation. This serves as a lower bound for black-hole masses in galaxies today. The comoving mass density of black hole seeds grows with time, and saturates when cosmic reionization stops gas cooling in these low-mass systems. By z~15, the comoving mass density becomes comparable to that inferred from observations, with room for appropriate additional luminous growth during a later quasar accretion phase. Hierarchical merging after z~15 naturally leads to a linear correlation between black-hole mass and stellar spheroid mass, with negligible black hole masses in disc-dominated galaxies.

If the favored hierarchical cosmological model is correct, then the Milky Way system should have accreted and subsequently tidally destroyed ~100 low-mass galaxies in the past ~12 Gyr. We model this process using a hybrid semi-analytic plus N-body approach and show that the disrupted systems lead naturally to stellar halos with masses and density profiles much like the stellar halo of our own Galaxy. We present predictions for the properties of stellar halos and show that ours is likely dominated by substructure beyond \~50kpc and more spatially smooth within that radius. The average stellar halo density profile is expected to drop off with radius more quickly than that of the dark matter because the stellar halo is formed from the most tightly bound material in infalling systems, while the majority of the accreted dark matter is stripped and deposited at larger radii. We argue that stars associated with the inner halo should be quite different chemically from stars in surviving satellites and also different from stars in the outer halo or those liberated in recent disruption events. We discuss how deep halo surveys and chemical probes may be useful tools for uncovering evidence of accretion. Searches of this kind offer a direct test of whether cosmology is indeed hierarchical on small scales.

(ABRIDGED) We present a semi-analytic model to explore merger histories, destruction rates, and survival probabilities of substructure in dark matter halos and use it to study the substructure populations of galaxy-sized halos as a function of the power spectrum. We successfully reproduce the subhalo velocity function and radial distribution seen in N-body simulations for standard LCDM. We explore the implications of spectra with normalizations and tilts spanning sigma_8 = 0.65-1 and n = 0.8-1. We also study a running index (RI) model with dn/dlnk=-0.03, as discussed in the first year WMAP report, and several WDM models with masses m_W = 0.75, 1.5, 3.0 keV. The substructure mass fraction is relatively insensitive to the tilt and overall normalization of the power spectrum. All CDM-type models yield projected substructure mass fractions that are consistent with, but on the low side of, estimates from strong lens systems: f = 0.4-1.5% (64 percentile) in systems M_sub < 10^9 Msun. Truncated models produce significantly smaller fractions and are disfavored by lensing results. We compare our predicted subhalo velocity functions to the dwarf satellite population of the Milky Way. Assuming isotropic velocity dispersions, we find the standard n=1 model overpredicts the number of MW satellites as expected. Models with less small-scale power are more successful because there are fewer subhalos of a given circular velocity and the mapping between observed velocity dispersion and halo circular velocity is markedly altered. The RI model, or a fixed tilt with sigma_8=0.75, can account for the MW dwarfs without the need for differential feedback; however, these comparisons depend sensitively on the assumption of isotropic velocities in satellite galaxies.

Reducing the power on small scales relative to the `standard' LCDM model alleviates a number of possible discrepancies with observations, and is favored by the recent analysis of WMAP plus galaxy and Lyman-alpha forest data. Here, we investigate the epoch of reionization in several models normalized to WMAP on large scales, and with sufficiently reduced power on small scales to solve the halo concentration and substructure problems. These include a tilted model, the WMAP running-index model, and a Warm Dark Matter model. We assume that the Universe was reionized by stellar sources comprised of a combination of supermassive (~200 M_sun) Pop III stars and Pop II stars with a `normal' IMF. We find that in all of these models, structure formation and hence reionization occurs late, certainly at redshifts below ten, and more probably at z < 6. This is inconsistent (at 2 sigma) with the determination of z_reion ~ 17 from the WMAP temperature-polarization data and is only marginally consistent with SDSS quasar observations. The tension between the galactic-scale observations, which favor low-power models, and the early reionization favored by WMAP can only be resolved if the efficiency of Pop III star formation is dramatically higher than any current estimate, or if there is an exotic population of ionizing sources such as mini-quasars. Otherwise, we may have to live with the standard LCDM power spectrum, and solve the small-scale problems in some other way.

In this proceeding, we present the results of a semi-analytic study of CDM substructure as a function of the primordial power spectrum. We apply our method to several tilted models in the LCDM framework with n=0.85-1.1, sigma_8=0.65-1.2 when COBE normalized. We also study a more extreme, warm dark matter-like spectrum that is sharply truncated below a scale of 10^10 h^-1 Msun. We show that the mass fraction of halo substructure is not a strong function of spectral slope, so it likely will be difficult to constrain tilt using flux ratios of gravitationally lensed quasars. On the positive side, all of our CDM-type models yield projected mass fractions in good agreement with strong lensing estimates: f \sim 1.5% at M \sim 10^8 Msun. The truncated model produces a significantly smaller fraction, f \lsim 0.3%, suggesting that warm dark matter-like spectra may be distinguished from CDM spectra using lensing. We also discuss the issue of dwarf satellite abundances, with emphasis on the cosmological dependence of the map between the observed central velocity dispersion of Milky Way satellites and the maximum circular velocities of their host halos. In agreement with earlier work, we find that standard LCDM over-predicts the estimated count of Milky Way satellites at fixed Vmax by an order of magnitude, but tilted models do better because subhalos are less concentrated. Interestingly, under the assumption that dwarfs have isotropic velocity dispersion tensors, models with significantly tilted spectra (n \lsim 0.85, sigma_8 \lsim 0.7) may under-predict the number of large Milky Way satellites with Vmax \gsim 40 km/s.

We determine the evolution of the faint, high-redshift, optical luminosity function (LF) of AGN implied by several observationally-motivated models of the ionizing background. Our results depend crucially on whether we use the total ionizing rate measured by the proximity effect technique or the lower determination from the flux decrement distribution of Ly alpha forest lines. Assuming a faint-end LF slope of 1.58 and the SDSS estimates of the bright-end slope and normalization, we find that the LF must break at M_B*=-24.2,-22.3, -20.8 at z=3,4, 5 if we adopt the lower ionization rate and assume no stellar contribution to the background. The break must occur at M_B*=-20.6,-18.7, -18.7 for the proximity effect background estimate. These values brighten by as much as ~2 mag if high-z galaxies contribute to the background with an escape fraction of ionizing photons consistent with recent estimates: f_e=0.16. By comparing to faint AGN searches, we find that the typically-quoted proximity effect estimates of the background imply an over-abundance of faint AGN (even with f_e=1). Even adopting the lower bound on proximity effect measurements, the stellar escape fraction must be high: f_e>0.2. Conversely, the lower flux- decrement-derived background requires a limited stellar contribution: f_e<0.05. Our derived LFs together with the locally-estimated black hole density suggest that the efficiency of converting mass to light in optically-unobscured AGN is somewhat lower than expected, <0.05. Comparison with similar estimates based on X-ray counts suggests that more than half of all AGN are obscured in the UV/optical. We also derive lower limits on typical AGN lifetimes and obtain >10^7 yrs for favored cases.

We propose a new explanation for the origin of angular momentum in galaxies and their dark halos, in which the halos obtain their spin through the cumulative acquisition of angular momentum from satellite accretion. In our model, the build-up of angular momentum is a random walk process associated with the mass assembly history of the halo's major progenitor. We assume no correlation between the angular momenta of accreted objects. Using the extended Press-Schechter approximation, we calculate the growth of mass, angular momentum, and spin parameter $\lambda$ for many halos. Our random walk model reproduces the key features of the angular momentum of halos found in N-body simulations: a lognormal distribution in $\lambda$ with an average of $<\lambda> \approx 0.04$, independent of mass and redshift. The evolution of the spin parameter in individual halos in this model is quite different from the steady increase with time of angular momentum in the tidal torque picture. We find both in N-body simulations and in our random walk model that the value of $\lambda$ changes significantly with time for a halo's major progenitor. It typically has a sharp increase due to major mergers, and a steady decline during periods of gradual accretion of small satellites. The model predicts that on average the $\lambda$ of halos which had major mergers after redshift $z=2$ should be substantially larger than the $\lambda$ of those which did not. Perhaps surprisingly, this suggests that halos that host late-forming elliptical galaxies should rotate faster than halos of spiral galaxies.

A problem with high central densities in dark halos has arisen in the context of LCDM cosmologies with scale-invariant initial power spectra. Although n=1 is often justified by appealing to the inflation scenario, inflationary models with mild deviations from scale-invariance are not uncommon and models with significant running of the spectral index are plausible. Even mild deviations from scale-invariance can be important because halo collapse times and densities depend on the relative amount of small-scale power. We choose several popular models of inflation and work out the ramifications for galaxy central densities. For each model, we calculate its COBE-normalized power spectrum and deduce the implied halo densities using a semi-analytic method calibrated against N-body simulations. We compare our predictions to a sample of dark matter-dominated galaxies using a non-parametric measure of the density. While standard n=1, LCDM halos are overdense by a factor of 6, several of our example inflation+CDM models predict halo densities well within the range preferred by observations. We also show how the presence of massive (0.5 eV) neutrinos may help to alleviate the central density problem even with n=1. We conclude that galaxy central densities may not be as problematic for the CDM paradigm as is sometimes assumed: rather than telling us something about the nature of the dark matter, galaxy rotation curves may be telling us something about inflation and/or neutrinos. An important test of this idea will be an eventual consensus on the value of sigma_8, the rms overdensity on the scale 8 h^-1 Mpc. Our successful models have values of sigma_8 approximately 0.75, which is within the range of recent determinations. Finally, models with n>1 (or sigma_8 > 1) are highly disfavored.

Although the currently favored cold dark matter plus cosmological constant model for structure formation assumes an n=1 scale-invariant initial power spectrum, most inflation models produce at least mild deviations from n=1. Because the lever arm from the CMB normalization to galaxy scales is long, even a small ``tilt'' can have important implications for galactic observations. Here we calculate the COBE-normalized power spectra for several well-motivated models of inflation and compute implications for the substructure content and central densities of galaxy halos. Using an analytic model, normalized against N-body simulations, we show that while halos in the standard (n=1) model are overdense by a factor of ~6 compared to observations, several of our example inflation+LCDM models predict halo densities well within the range of observations, which prefer models with n ~ 0.85. We go on to use a semi-analytic model (also normalized against N-body simulations) to follow the merger histories of galaxy-sized halos and track the orbital decay, disruption, and evolution of the merging substructure. Models with n ~0.85 predict a factor of ~3 fewer subhalos at a fixed circular velocity than the standard $n = 1$ case. Although this level of reduction does not resolve the ``dwarf satellite problem'', it does imply that the level of feedback required to match the observed number of dwarfs is sensitive to the initial power spectrum. Finally, the fraction of galaxy-halo mass that is bound up in substructure is consistent with limits imposed by multiply imaged quasars for all models considered: f_sub > 0.01 even for an effective tilt of n ~0.8.We conclude that, at their current level, lensing constraints of this kind do not provide strong limits on the primordial power spectrum.

(abridged) We study the relation between the density profiles of dark matter halos and their mass assembly histories, using a statistical sample of halos in a high-resolution N-body simulation of the LCDM cosmology. For each halo at z=0, we identify its merger-history tree, and determine concentration parameters c_vir for all progenitors, thus providing a structural merger tree for each halo. We fit the mass accretion histories by a universal function with one parameter, the formation epoch a_c, defined when the log mass accretion rate dlogM/dloga falls below a critical value S. We find that late forming galaxies tend to be less concentrated, such that c_vir ``observed'' at any epoch a_o is strongly correlated with a_c via c_vir=c_1*a_o/a_c. Scatter about this relation is mostly due to measurement errors in c_v and a_c, implying that the actual spread in c_vir for halos of a given mass can be mostly attributed to scatter in a_c. We demonstrate that this relation can also be used to predict the mass and redshift dependence of c_v, and the scatter about the median c_vir(M,z), using accretion histories derived from the Extended Press-Schechter (EPS) formalism, after adjusting for a constant offset between the formation times as predicted by EPS and as measured in the simulations;this new ingredient can thus be easily incorporated into semi-analytic models of galaxy formation. The correlation found between halo concentration and mass accretion rate suggests a physical interpretation: for high mass infall rates the central density is related to the background density; when the mass infall rate slows, the central density stays approximately constant and the halo concentration just grows as R_vir. The tight correlation demonstrated here provides an essential new ingredient for galaxy formation modeling.

We study the relation between the density profiles of dark matter halos and their mass assembly histories, using a statistical sample of halos in a high-resolution N-body simulation of the LCDM cosmology. For each halo at z=0, we identify its merger-history tree, and determine concentration parameters c_vir for all progenitors, thus providing a structural merger tree for each halo. We fit the mass accretion histories by a universal function with one parameter, the formation epoch a_c, defined when the log mass accretion rate dlog M/doga falls below a critical value. We find that late forming galaxies tend to be less concentrated, such that c_vir ``observed'' at any epoch a_obs is strongly correlated with a_c via c_vir=c_1 a_obs/a_c. Scatter about this relation is mostly due to measurement errors in c_vir and a_c, implying that the actual spread in c_vir for halos of a given mass can be mostly attributed to scatter in a_c. Because of the direct connection between halo concentration and velocity rotation curves, and because of probable connections between halo mass assembly history and star formation history, the tight correlation between these properties provides an essential new ingredient for galaxy formation modeling.

Although the currently favored cold dark matter plus cosmological constant model (LCDM) has proven to be remarkably successful on large scales, on subgalactic scales it faces some potentially fatal difficulties; these include over-producing dwarf satellite galaxies and predicting excessive central densities in dark halos. Among the most natural cosmological solutions to these problems is to replace cold dark matter with a warm species (LWDM). The warm component acts to reduce the small-scale power, resulting in fewer galactic subhalos and lower halo central densities. An alternative model with a mild ``tilt'' in the inflationary power spectrum (TLCDM; n =0.9) similarly reduces the central densities of dark halos, although the substructure abundance remains relatively high. Here I argue that because dwarf galaxy formation should be suppressed in the presence of a strong ionizing background, favored LWDM models will generally under-predict the observed abundance of dwarf galaxies. The satellite count for TLCDM fairs much better, as long as the photoionization effect is taken into account. TLCDM provides a more successful alternative to LWDM on subgalactic scales with the added attraction that it relies on only a minor, natural adjustment to the standard framework of CDM.

We compare the specific angular momentum profiles of virialized dark halos in cold dark matter (CDM) and warm dark matter (WDM) models using high-resolution dissipationless simulations. The simulations were initialized using the same set of modes, except on small scales, where the power was suppressed in WDM below the filtering length. Remarkably, WDM as well as CDM halos are well-described by the two-parameter angular momentum profile of Bullock et al. (2001), even though the halo masses are below the filtering scale of the WDM. Although the best-fit shape parameters change quantitatively for individual halos in the two simulations, we find no systematic variation in profile shapes as a function of the dark matter type. The scatter in shape parameters is significantly smaller for the WDM halos, suggesting that substructure and/or merging history plays a role producing scatter about the mean angular momentum distribution, but that the average angular momentum profiles of halos originate from larger-scale phenomena or a mechanism associated with the virialization process. The known mismatch between the angular momentum distributions of dark halos and disk galaxies is therefore present in WDM as well as CDM models. Our WDM halos tend to have a less coherent (more misaligned) angular momentum structure and smaller spin parameters than do their CDM counterparts, although we caution that this result is based on a small number of halos.

Using an analytic model calibrated against numerical simulations, we calculate the central densities of dark matter halos in a ``conventional'' cold dark matter model with a cosmological constant (LCDM) and in a ``tilted'' model (TLCDM) with slightly modified parameters motivated by recent analyses of Ly-alpha forest data. We also calculate how warm dark matter (WDM) would modify these predicted densities by delaying halo formation and imposing phase space constraints. As a measure of central density, we adopt the quantity D_{V/2}, the density within the radius R_{V/2} at which the halo rotation curve falls to half of its maximum value, in units of the critical density. We compare the theoretical predictions to values of D_{V/2} estimated from the rotation curves of dark matter dominated disk galaxies. Assuming that dark halos are described by NFW profiles, our results suggest that the conventional LCDM model predicts excessively high dark matter densities, unless there is some selection bias in the data toward the low-concentration tail of the halo distribution. A WDM model with particle mass 0.5-1 keV provides a better match to the observational data. However, the modified cold dark matter model, TLCDM, fits the data equally well, suggesting that the solution to the ``halo cores'' problem might lie in moderate changes to cosmological parameters rather than radical changes to the properties of dark matter. If CDM halos have the steeper density profiles found by Moore et al., then neither conventional LCDM nor TLCDM can reproduce the observed central densities.

I present an analysis of the density shapes of dark matter halos in LCDM and LWDM cosmologies. The main results are derived from a statistical sample of galaxy-mass halos drawn from a high resolution LCDM N-body simulation. Halo shapes show significant trends with mass and redshift: low-mass halos are rounder than high mass halos, and, for a fixed mass, halos are rounder at low z. Contrary to previous expectations, which were based on cluster-mass halos and non-COBE normalized simulations, LCDM galaxy-mass halos at z=0 are not strongly flattened, with short to long axis ratios of s = 0.70 +/- 0.17. I go on to study how the shapes of individual halos change when going from a LCDM simulation to a simulation with a warm dark matter power spectrum (LWDM). Four halos were compared, and, on average, the WDM halos are more spherical than their CDM counterparts (s =0.77 compared to s = 0.71). A larger sample of objects will be needed to test whether the trend is significant.

We discuss how current and future data on the clustering and number density of z~3 Lyman-break galaxies (LBGs) can be used to constrain their relationship to dark matter haloes. We explore a three-parameter model in which the number of LBGs per dark halo scales like a power-law in the halo mass: N(M) = (M/M_1)^S for M>M_m. Here, M_m is the minimum mass halo that can host an LBG, M_1 is a normalization parameter, associated with the mass above which haloes host more than one observed LBG, and S determines the strength of the mass dependence. We show how these three parameters are constrained by three observable properties of LBGs: the number density, the large-scale bias, and the fraction of objects in close pairs. Given these three quantities, the three unknown model parameters may be estimated analytically, allowing a full exploration of parameter space. As an example, we assume an LCDM cosmology and consider the observed properties of a recent sample of spectroscopically confirmed LBGs. We find that the favored range for our model parameters is M_m = (0.4-8)x10^10 Msun/h, M_1 = (6-10)x10^12 Msun/h, and 0.9 < S < 1.1. The preferred region in M_m expands by an order of magnitude and slightly shallower slopes are acceptable if the allowed range of bias is permitted to span all recent observational estimates. We also discuss how the observed clustering of LBGs as a function of luminosity can be used to constrain halo occupation, although due to current observational uncertainties we are unable to reach any strong conclusions. Our methods and results can be used to constrain more realistic models that aim to derive the occupation function N(M) from first principles, and offer insight into how basic physical properties affect the observed properties of LBGs.

We examine the detectability and interpretation of debris trails caused by satellite disruption in external galaxies using semi-analytic approximations for the dependence of streamer length, width and surface brightness on satellite and primary galaxy characteristics. The semi-analytic method is tested successfully against N-body simulations and then applied to three representative astronomical applications. First, we show how streamer properties can be used to estimate mass-to-light ratios $\Upsilon$ and streamer ages of totally disrupted satellites, and apply the method to the stellar arc in NGC 5907. Second, we discuss how the lack of observed tidal debris around a satellite can provide an upper limit on its mass-loss rate, and, as an example, derive the implied limits on mass-loss rates for M32 and NGC 205 around Andromeda. Finally, we point out that a statistical analysis of streamer properties might be applied to test and refine cosmological models of hierarchical galaxy formation, and use the predicted debris from a standard $\Lambda$CDM realization to test the feasibility of such a study. Using the Local Group satellites and the few known examples of debris trails in the Galaxy and in external systems, we estimate that the best current techniques could characterize the brightest ($R < 29$ mag/ arcsec$^{2}$) portions of the youngest (3 dynamical periods) debris streamers. If systematics can be controlled, planned large-aperture telescopes such as CELT and OWL may allow fainter trails to be detected routinely and thus used for statistical studies such as those required for tests of galaxy formation.

We study the effects of mass and force resolution on the density profiles of galaxy-size Cold Dark Matter (CDM) halos in a flat, low-density cosmological model with vacuum energy. We show thatalthough increasing the mass and force resolution allows us to probe deeper into the inner halo regions, it does not lead to steeper inner density profiles. Instead, the halo profiles converge at scales larger than four times the formal resolution or the radius containing more than 200 particles, whichever is larger. In the simulations presented in this paper, we are able to probe density profile of a relaxed isolated galaxy-size halo at scales r=(0.005-1)r_vir. We find that the density distribution can be well approximated by the profile suggested by Moore etal (1998): rho= x^{-1.5}(1+x^{1.5})^-1, where x=r/r_s and r_s is the characteristic radius. The analytical profile proposed by Navarro et al. (1996) rho= x^{-1}(1+x)^-2, also provides a good fit, with the same relative errors of about 10% for radii larger than 1% of the virial radius. Both analytical profiles fit equally well because for high-concentration galaxy-size halos the differences between these profiles become significant only at scales well below 0.01r_vir. We also find that halos of similar mass may have somewhat different parameters (characteristic radius, maximum rotation velocity, etc.) and shapes of their density profiles. We associate this scatter in properties with differences in halo merger histories and the amount of substructure present in the analyzed halos.

We use N-body simulations combined with semi-analytic models to compute the clustering properties of modeled galaxies at z~3, and confront these predictions with the clustering properties of the observed population of Lyman-break galaxies (LBGs). Several scenarios for the nature of LBGs are explored, which may be broadly categorized into models in which high-redshift star formation is driven by collisional starbursts and those in which quiescent star formation dominates. For each model, we make predictions for the LBG overdensity distribution, the variance of counts-in-cells, the correlation length, and close pair statistics. Models which assume a one-to-one relationship between massive dark-matter halos and galaxies are disfavored by close pair statistics, as are models in which colliding halos are associated with galaxies in a simplified way. However, when modeling of gas consumption and star formation is included using a semi-analytic treatment, the quiescent and collisional starburst models predict similar clustering properties and none of these models can be ruled out based on the available clustering data. None of the ``realistic'' models predict a strong dependence of clustering amplitude on the luminosity threshold of the sample, in apparent conflict with some observational results.

We discuss how measurements of the absorption of $\gamma$-rays from GeV to TeV energies via pair production on the extragalactic background light (EBL) can probe important issues in galaxy formation. We use semi-analytic models (SAMs) of galaxy formation, set within the hierarchical structure formation scenario, to obtain predictions of the EBL for 0.1-1000$\mu$m. SAMs incorporate simplified physical treatments of the key processes of galaxy formation --- including gravitational collapse and merging of dark matter halos, gas cooling and dissipation, star formation, supernova feedback and metal production --- and have been shown to reproduce key observations at low and high redshift. Here we also introduce improved modelling of the spectral energy distributions in the mid-to-far-IR arising from emission by dust grains. Assuming a flat \lcdm cosmology with $\Omega_m=0.3$ and Hubble parameter $h=0.65$, we investigate the consequences of variations in input assumptions such as the stellar initial mass function (IMF) and the efficiency of converting cold gas into stars. We conclude that observational studies of the absorption of $\gamma$-rays with energies from 10s of Gev to 10s of TeV will help to determine the EBL, and also help to explain its origin by constraining some of the most uncertain features of galaxy formation theory, including the IMF, the history of star formation, and the reprocessing of light by dust.

[Abridged] We study the angular-momentum profiles of a statistical sample of halos drawn from a high-resolution N-body simulation of the LCDM cosmology. We find that the cumulative mass distribution of specific angular momentum, j, in a halo of mass Mv is well fit by a universal function, M(<j) = Mv \mu j/(j_0+j). This profile is defined by one shape parameter (\mu or j_0) in addition to the global spin parameter \lambda. It follows a power-law over most of the mass, and flattens at large j, with the flattening more pronounced for small values of \mu. Compared to a uniform sphere in solid-body rotation, most halos have a higher fraction of their mass in the low- and high-j tails of the distribution. The spatial distribution of angular momentum in halos tends to be cylindrical and is well-aligned within each halo for ~80% of the halos. We investigate two ideas for the origin of this profile. The first is based on a revised version of linear tidal-torque theory combined with extended Press-Schechter mass accretion, and the second focuses on j transport in minor mergers. Finally, we briefly explore implications of the M(<j) profile on the formation of galactic disks assuming that j is conserved during an adiabatic baryonic infall. The implied gas density profile deviates from an exponential disk, with a higher density at small radii and a tail extending to large radii. The steep central density profiles may imply disk scale lengths that are smaller than observed. This is reminiscent of the "angular-momentum problem" seen in hydrodynamic simulations, even though we have assumed perfect j conservation. A possible solution is to associate the central excesses with bulge components and the outer regions with extended gaseous disks.

We summarize recent developments in the study of the origin of halo spin profiles and preliminary implications on disk formation. The specific angular-momentum distributions within halos in N-body simulations match a universal shape, M(<j) \propto j/(j_0+j). It is characterized by a power law over most of the mass, and one shape parameter in addition to the spin parameter \lambda. The angular momentum tends to be aligned throughout the halo and of cylindrical symmetry. Even if angular momentum is conserved during baryonic infall, the resultant disk density profile is predicted to deviate from exponential, with a denser core and an extended tail. A slightly corrected version of the scaling relation due to linear tidal-torque theory is used to explain the origin of a typical power-law profile in shells, j(M) \propto M^s with s \gsim 1. While linear theory crudely predicts the amplitudes of halo spins, it is not a good predictor of their directions. Independently, mergers of halos are found to produce a similar profile due to j transfer from the orbit to the product halo via dynamical friction and tidal stripping. The halo spin is correlated with having a recent major merger, though this correlation is weakened by mass loss. These two effects are due to a correlation between the spins of neighboring halos and their orbit, leading to prograde mergers.

[Abridged] The interaction rates of dark-matter halos and subhalos, are computed using high-resolution cosmological N-body simulations of the Lambda-CDM model. Although the number fraction of subhalos of mass $>2x10^{11}\hsolmass$ is only ~10%, we find that the interaction rate of such subhalos is relatively high because they reside in high density environments. At low redshift, the subhalo collisions dominate the total collision rate, and even at z=3 they are involved in more than 30% of all collisions. About 40% of the "major" mergers are between subhalos. Therefore subhalo interactions must be incorporated in models of structure formation. We find that the collision rate between halos in physical density units, is $\propto (1+z)^\delta $, with $\delta = 3-4$, in agreement with earlier simulations and most observational data. We test previous analytic estimates of the interaction rates of subhalos based on statistical models, which could be very inaccurate because of the small number of subhalos and the variation of conditions within small host halos. We find that, while such statistical estimates may severely overestimate the rate within hosts $< 10^{13} \hsolmass$, typical of high redshifts, they are valid for larger hosts regardless of the number of subhalos in them. We find the Makino & Hut (1997) estimate of the subhalo merger rate to be valid for hosts $ > 6x10^{11}\hsolmass$ at all redshifts. The collision rate between subhalos and the central object of their host halo is approximated relatively well using the timescale for dynamical friction in circular orbits. This approximation fails in ~40% of the cases, partly because of deviations from the assumption of circular orbits and partly because of the invalidity of the assumption that the subhalo mass is negligible.

[Abridged] We analyze the exchange of dark matter between halos, subhalos, and their environments in a high-resolution cosmological N-body simulation of a Lambda CDM cosmology. At each analyzed redshift z we divide the dark matter particles into 4 components: (i) isolated galactic halos, (ii) subhalos, (iii) the diffuse medium of group and cluster halos, and (iv) the background outside of virialized halos. We follow the time evolution of the mass distribution and flows between these components and provide fitting functions for the exchange rates. We use our derived exchange rates to gauge the importance of metal redistribution in the universe due solely to gravity-induced interactions. The diffuse metallicity in clusters is predicted to be ~40% that in isolated galaxies (~55% of groups) at z=0, and should be lower only slightly by z=1, consistent with observations. The metallicity of the diffuse media in poor groups is expected to be lower by a factor of 5 by z~2, in agreement with the observed metallicity of damped Ly$\alpha$ systems. The metallicity of the background IGM is predicted to be (1-3)x10^{-4} that of z=0 clusters, also consistent with observations. The agreement of predicted and observed trends indicates that gravitational interaction alone may play an important role in metal enrichment of the intra-cluster and intergalactic media.

We study dark-matter halo density profiles in a high-resolution N-body simulation of an LCDM cosmology. Our statistical sample contains ~5000 haloes in the range 10^{11}-10^{14} M_sun. The profiles are parameterized by an NFW form with two parameters, an inner radius r_s and a virial radius r_v, and we define the halo concentration c_v = r_v/r_s. First, we find that, for a given halo mass, the redshift dependence of the median concentration is c_v ~ 1/(1+z), corresponding to a roughly constant r_s with redshift. We present an improved analytic treatment of halo formation that fits the measured relations between halo parameters and their redshift dependence. The implications are that high-redshift galaxies are predicted to be more extended and dimmer than expected before. Second, we find that the scatter in log(c_v) is ~0.18, corresponding to a scatter in maximum rotation velocities of dV/V ~ 0.12. We discuss implications for modelling the Tully-Fisher relation, which has a smaller reported intrinsic scatter. Third, haloes in dense environments tend to be more concentrated than isolated haloes. These results suggest that c_v is an essential parameter for the theory of galaxy modelling, and we briefly discuss implications for the universality of the Tully-Fisher relation, the formation of low surface brightness galaxies, and the origin of the Hubble sequence.

We develop an explicit model for the formation of the stellar halo from tidally disrupted, accreted dwarf satellites in the cold dark matter (CDM) framework, focusing on predictions testable with the Sloan Digital Sky Survey (SDSS) and other wide-field surveys. Subhalo accretion and orbital evolution are calculated using a semi-analytic approach within the Press-Schechter formalism. Motivated by our previous work, we assume that low-mass subhalos (v < 30 km/s) can form significant populations of stars only if they accreted a substantial fraction of their mass before the epoch of reionization. With this assumption, the model reproduces the observed velocity function of galactic satellites in the Local Group, solving the ``dwarf satellite problem'' without modifying the popular LCDM cosmology. The disrupted satellites yield a stellar distribution with a total mass and radial density profile consistent with those observed for the Milky Way stellar halo. Most significantly, the model predicts the presence of many large-scale, coherent substructures in the outer halo. These substructures are remnants of individual, tidally disrupted dwarf satellite galaxies. Substructure is more pronounced at large galactocentric radii because of the smaller number density of tidal streams and the longer orbital times. This model provides a natural explanation for the coherent structures in the outer stellar halo found in the SDSS commissioning data, and it predicts that many more such structures should be found as the survey covers more of the sky. The detection (or non-detection) and characterization of such structures could eventually test variants of the CDM scenario, especially those that aim to solve the dwarf satellite problem by enhancing satellite disruption.

We present a method for constraining the evolution of the galaxy luminosity-velocity (LV) relation in hierarchical scenarios of structure formation. The comoving number density of dark-matter halos with circular velocity of 200 km/s is predicted in favored CDM cosmologies to be nearly constant over the redshift range 0<z<5. Any observed evolution in the density of bright galaxies implies in turn a corresponding evolution in the LV relation. We consider several possible forms of evolution for the zero-point of the LV relation and predict the corresponding evolution in galaxy number density. The Hubble Deep Field suggests a large deficit of bright (M_V < -19) galaxies at 1.4 < z < 2. If taken at face value, this implies a dimming of the LV zero-point by roughly 2 magnitudes. Deep, wide-field, near-IR selected surveys will provide more secure measurements to compare with our predictions.

We study the distribution functions of mass and circular velocity for dark matter halos in N-body simulations of the $\Lambda$CDM cosmology, addressing redshift and environmental dependence. The dynamical range enables us to resolve subhalos and distinguish them from "distinct" halos. The mass function is compared to analytic models, and is used to derive the more observationally relevant circular velocity function. The distribution functions in the velocity range 100--500 km/s are well fit by a power-law with two parameters, slope and amplitude. We present the parameter dependence on redshift and provide useful fitting formulae. The amplitudes of the mass functions decrease with z, but, contrary to naive expectation, the comoving density of halos of a fixed velocity ~200 km/s actually increases out to z=5. This is because high-z halos are denser, so a fixed velocity corresponds to a smaller mass. The slope of the velocity function at z=0 is as steep as ~ -4, and the mass and velocity functions of distinct halos steepen with increasing z, while the functions of subhalos do not steepen with $z$, and become even flatter at z>2. A simple observable prediction is that the slope of the velocity function of isolated galaxies is steeper than that of galaxies in groups by as much as unity, reflecting the density biasing of high-velocity halos. We confirm that the Press-Schechter approximation typically overestimates the halo mass function by a factor of ~2, while modified approximations provide improved predictions.

One of the main challenges facing standard hierarchical structure formation models is that the predicted abundance of galactic subhalos with circular velocities of 10-30 km/s is an order of magnitude higher than the number of satellites actually observed within the Local Group. Using a simple model for the formation and evolution of dark halos, based on the extended Press-Schechter formalism and tested against N-body results, we show that the theoretical predictions can be reconciled with observations if gas accretion in low-mass halos is suppressed after the epoch of reionization. In this picture, the observed dwarf satellites correspond to the small fraction of halos that accreted substantial amounts of gas before reionization. The photoionization mechanism naturally explains why the discrepancy between predicted halos and observed satellites sets in at about 30 km/s, and for reasonable choices of the reionization redshift (z_re = 5-12) the model can reproduce both the amplitude and shape of the observed velocity function of galactic satellites. If this explanation is correct, then typical bright galaxy halos contain many low-mass dark matter subhalos. These might be detectable through their gravitational lensing effects, through their influence on stellar disks, or as dwarf satellites with very high mass-to-light ratios. This model also predicts a diffuse stellar component produced by large numbers of tidally disrupted dwarfs, perhaps sufficient to account for most of the Milky Way's stellar halo.

We investigate the clustering properties of high-redshift galaxies within three competing scenarios for assigning luminous galaxies to dark matter halos from N-body simulations: a one galaxy per massive halo model, a quiescent star formation model, and a collisional starburst model. We compare these models to observations of Lyman-Break galaxies at z~3$ With current data and the simple statistic used here, one cannot rule out any of these models, but we see potential for finding distinguishing features using statistics that are sensitive to the tails of the distribution, and statistics based on the number of multiple galaxies per halo, which we explore in an ongoing study.

This paper represents a step toward constraining galaxy formation models via TeV gamm a ray observations. We use semi-analytic models of galaxy formation to predict a spectral distribution for the intergalactic infrared photon field, which in turn yields information about the absorption of TeV gamma rays from extra-galactic sources. By making predictions for integral flux observations at >200 GeV for several known EGRE T sources, we directly compare our models with current observational upper limits obtained by Whipple. In addition, our predictions may offer a guide to the observing programs for the current population of TeV gamma ray observatories.

We present a galaxy circular velocity function, Psi(log v), derived from existing luminosity functions and luminosity-velocity relations. Such a velocity function is desirable for several reasons. First, it enables an objective comparison of luminosity functions obtained in different bands and for different galaxy morphologies, with a statistical correction for dust extinction. In addition, the velocity function simplifies comparison of observations with predictions from high-resolution cosmological N-body simulations. We derive velocity functions from five different data sets and find rough agreement among them, but about a factor of 2 variation in amplitude. These velocity functions are then compared with N-body simulations of a LCDM model (corrected for baryonic infall) in order to demonstrate both the utility and current limitations of this approach. The number density of dark matter halos and the slope of the velocity function near v_*, the circular velocity corresponding to an ~L_* spiral galaxy, are found to be comparable to that of observed galaxies. The primary sources of uncertainty in construction of Psi(log v) from observations and N-body simulations are discussed and explanations are suggected to account for these discrepancies.

Lyman break galaxies (LBGs) at $z \sim 3$ exhibit number densities and clustering similar to local $L_*$ galaxies; however, their star formation rates (SFRs) are much higher. We explore the scenario in which LBGs are starburst galaxies triggered by collisions, and thus provide an explanation for these key properties. The relative paucity of starburst galaxies at low redshift can be attributed to a much higher collision rate in the past. We use high-resolution cosmological $N$-body simulations and a hierarchical halo finder to estimate the galaxy collision rate as a function of time in the popular $\Lambda$CDM cosmological model. We find that bright collisional starbursts are frequent enough to account for most of the high-$z$ ($2.5-4.5$) LBGs. Moreover, many of the objects are of relatively small mass, but they cluster about large-mass halos. They therefore exhibit strong clustering, similar to that observed and stronger than that of the relevant massive halos.

We present here the extragalactic background light (EBL) predicted by semi-analytic models of galaxy formation, and show how measurements of the absorption of gamma rays of $\sim$ TeV energies via pair production on the EBL can probe cosmology and the formation of galaxies. Semi-analytic models permit a physical treatment of the key processes of galaxy formation -- including gravitational collapse and merging of dark matter halos, gas cooling and dissipation, star formation, supernova feedback and metal production -- and have been shown to reproduce key observations at low and high redshift. Using this approach, we investigate the consequences of variations in input assumptions such as the stellar initial mass function and the underlying cosmology. We conclude that observational studies of the absorption of $\sim 10^{-2}-10^{2}$ TeV gamma rays will help to constrain the star formation history of the universe, and the nature and extent of the extinction of starlight due to dust and reradiation of the absorbed energy at infrared wavelengths.

Gamma rays of ~TeV energies from distant sources suffer attenuation due to pair production off of ~1 micron EBL photons. We may exploit this process in order to indirectly measure the EBL and constrain models of galaxy formation. Here, using semi-analytic models of galaxy formation, we examine how gamma ray absorption may be used as an indirect probe of the stellar initial mass function (IMF), although there is a degeneracy with dust modeling. We point out that with the new generation of gamma ray telescopes including STACEE, MAGIC, HESS, VERITAS, and Milagro, we should soon possess a wealth of new data and a new method for probing the nature of the IMF.

This article concerns the formation and structure of dark matter halos, including (1) their radial density profiles, (2) their abundance, and (3) their merger rates. The last topic may be relevant to the nature of the small, bright, high-redshift galaxies discovered by the Lyman break technique. (1) Study of a statistical sample of galaxy-mass dark halos in high-resolution Adaptive Refinement Tree simulations shows that they have a central density profile p(r) \propto 1/r^g with g \approx 0.2, in agreement with data on dark-matter-dominated disk galaxies. We present recent, higher resolution results on this. (2) Another important new result is that the Press-Schechter approximation predicts about twice as many galaxy-mass halos at z=0 as are present in large dissipationless N-body simulations; more generally, PS overpredicts the abundance of M ~< 0.1M_* halos at all redshifts. (3) Finally, we discuss the assembly of these halos, in particular the merger rate of (sub-)halos at high redshift and the distribution of the starbursts that these mergers are likely to trigger. If most of the Lyman-break galaxies are such starbursts, this perhaps resolves the apparent paradox that these galaxies appear to cluster like massive halos (~10^{12} M_\odot), while their relatively low linewidths and their spectral energy distributions suggest that they have relatively low mass (few X 10^{10} M_\odot) and young ages (few X 10^8 yr). It also predicts much more star formation at high redshift in CDM-type hierarchical models for structure formation than if only quiescent star formation is included.

We review the basic arguments for the likelihood of non-Gaussian density perturbations in inflation models with primordial black hole (PBH) production. We discuss our derived distributions of field fluctuations and their implications, specifically commenting on the fine-tuning problem. We also discuss how the derived distributions may be affected when linked to metric perturbations. While linking the metric perturbations to field fluctuations in a nonlinear way may be important for determining exact probability distributions, the correct mapping is not self-evident. The calculation of P. Ivanov, which yields skew positive distribution, is based on an ansatz for the behavior of the nonlinear metric perturbation. We note that the ``natural'' generalization of the gauge-invariant formalism favored by Bond and Salopek yields an effective linear link between the distribution of field fluctuations and metric perturbations during inflation.

We use the rotation curves of a sample of dark matter dominated dwarf and low-surface brightness (LSB) late-type galaxies to study their radial mass distributions. We find that the shape of the rotation curves is remarkably similar for all (both dwarf and LSB) galaxies in the sample, suggesting a self-similar density distribution of their dark matter (DM) halos. This shape can be reproduced well by a density profile with a shallow central cusp (rho(r)\propto 1/r^gamma, gamma\approx 0.2-0.4), corresponding to a steeply rising velocity curve (v(r)\propto r^g, g\approx 0.9-0.8). We further show that the observed shape of the rotation curves is well matched by the density profiles of DM halos formed in very high resolution simulations of the CDM, LCDM, and CHDM models of structure formation. This is surprising in light of several previous studies which suggested that the structure of simulated DM halos is inconsistent with the dynamics of dwarf galaxies. We discuss possible explanations for this discrepancy. We show that although the mass distribution in the hierarchically formed halos is on average consistent with the shape of rotation curves of dark matter dominated galaxies, the scatter of the individual profiles around the average is substantial and should not be neglected in comparisons with the data. Finally, we show that the DM halos in our hierarchical simulations and the real galaxies in our sample exhibit a similar decrease in their characteristic densities with increasing characteristic radial scales, and increase in their maximum rotation velocities with increasing radii at which their maximum velocity occurs. (Abridged)

We explore the role of non-Gaussian fluctuations in primordial black hole (PBH) formation and show that the standard Gaussian assumption, used in all PBH formation papers to date, is not justified. Since large spikes in power are usually associated with flat regions of the inflaton potential, quantum fluctuations become more important in the field dynamics, leading to mode-mode coupling and non-Gaussian statistics. Moreover, PBH production requires several sigma (rare) fluctuations in order to prevent premature matter dominance of the universe, so we are necessarily concerned with distribution tails, where any intrinsic skewness will be especially important. We quantify this argument by using the stochastic slow-roll equation and a relatively simple analytic method to obtain the final distribution of fluctuations. We work out several examples with toy models that produce PBH's, and show that the naive Gaussian assumption can result in errors of many orders of magnitude. For models with spikes in power, our calculations give sharp cut-offs in the probability of large positive fluctuations, meaning that Gaussian distributions would vastly over-produce PBH's. The standard results that link inflation-produced power spectra and PBH number densities must then be reconsidered, since they rely quite heavily on the Gaussian assumption. We point out that since the probability distributions depend on the nature of the potential, it is impossible to obtain results for general models. However, calculating the distribution of fluctuations for any specific model seems to be relatively straightforward, at least in the single inflaton case.