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.

Models of galaxy formation in a cosmological context predict that massive disk galaxies should have structured extended stellar halos. Recent studies in integrated light, however, report a few galaxies, including the nearby disk galaxy M101, that have no measurable stellar halos to the detection limit. We aim to quantify the stellar content and structure of M101's outskirts by resolving its stars. We present the photometry of its stars based on deep F606W and F814W images taken with Hubble Space Telescope as part of the GHOSTS survey. The constructed CMDs of stars reach down to two magnitudes below the tip of the red giant branch. We derived radial number density profiles of the bright red giant branch (RGB) stars. The mean color of the RGB stars at $R \sim$ 40 -- 60 kpc is similar to those of metal-poor globular clusters in the Milky Way. We also derived radial surface brightness profiles using the public image data provided by the Dragonfly team. Both the radial number density and surface brightness profiles were converted to radial mass density profiles and combined. We find that the mass density profiles show a weak upturn at the very outer region, where surface brightness is as faint as $\mu_g\approx 34$ mag arcsec$^{-1}$. An exponential disk + power-law halo model on the mass density profiles finds the total stellar halo mass of $M_{halo}=8.2_{-2.2}^{+3.5}\times 10^7M_\odot$. The total stellar halo mass does not exceed $M_{halo} = 3.2 \times 10^8$ $M_{\odot}$ when strongly truncated disk models are considered. Combining the halo mass with the total stellar mass of M101, we obtain the stellar halo mass fraction of $M_{halo}/M_{gal} = 0.20_{-0.08}^{+0.10}\%$ with an upper limit of 0.78\%. We compare the halo properties of M101 with those of six GHOSTS survey galaxies as well as the Milky Way and M31 and find that M101 has an anemic stellar halo.

Surveys of the Milky Way (MW) and M31 enable detailed studies of stellar populations across ages and metallicities, with the goal of reconstructing formation histories across cosmic time. These surveys motivate key questions for galactic archaeology in a cosmological context: when did the main progenitor of a MW/M31-mass galaxy form, and what were the galactic building blocks that formed it? We investigate the formation times and progenitor galaxies of MW/M31-mass galaxies using the FIRE-2 cosmological simulations, including 6 isolated MW/M31-mass galaxies and 6 galaxies in Local Group (LG)-like pairs at z = 0. We examine main progenitor "formation" based on two metrics: (1) transition from primarily ex-situ to in-situ stellar mass growth and (2) mass dominance compared to other progenitors. We find that the main progenitor of a MW/M31-mass galaxy emerged typically at z ~ 3-4 (11.6-12.2 Gyr ago), while stars in the bulge region (inner 2 kpc) at z = 0 formed primarily in a single main progenitor at z < 5 (< 12.6 Gyr ago). Compared with isolated hosts, the main progenitors of LG-like paired hosts emerged significantly earlier (\Delta z ~ 2, \Delta t ~ 1.6 Gyr), with ~ 4x higher stellar mass at all z > 4 (> 12.2 Gyr ago). This highlights the importance of environment in MW/M31-mass galaxy formation, especially at early times. Overall, about 100 galaxies with M_star > 10^5 M_sun formed a typical MW/M31-mass system. Thus, surviving satellites represent a highly incomplete census (by ~ 5x) of the progenitor population.

Recent work has shown that Milky Way-mass galaxies display an incredible range of stellar halo properties. Yet, the origin of this diversity is unclear. The nearby galaxy M81 $-$ currently interacting with M82 and NGC 3077 $-$ sheds unique light on this problem. We present a Subaru Hyper Suprime-Cam survey of the resolved stellar populations in M81's halo. Using a unique three-filter strategy, and numerous HST calibration fields, we reveal M81's halo in never-before-seen detail. We resolve the halo to unprecedented $V$-band equivalent surface brightnesses of $>$34 mag arcsec$^{-2}$, and produce the first-ever global stellar mass density map for a Milky Way-mass stellar halo outside of the Local Group (LG). Using the minor axis, we confirm the previous assessment of M81 as one of the lowest mass and metal-poorest stellar halos known ($M_{\star} \sim 10^9 M_{\odot}$, [Fe/H] $\sim {-}1.2$) $-$ indicating a relatively quiet prior accretion history. Yet, in our global stellar mass census we find that tidally unbound material from M82 and NGC 3077 provides a substantial infusion of metal-rich material ($M_{\star} \simeq 6{\times}10^8$ $M_{\odot}$, [Fe/H] $\simeq {-}$0.9). We further show that, following the accretion of its massive satellite M82 (and the LMC-like NGC 3077), M81 will host one of the most massive and metal-rich stellar halos in the nearby universe. Thus, the saga of M81: following a relatively passive history, M81's merger with M82 will completely transform its halo from a low-mass, anemic halo rivaling the MW, to a metal-rich behemoth rivaled only by systems such as M31. This dramatic transformation indicates that the observed diversity in stellar halo properties is primarily driven by diversity in the largest mergers these galaxies have experienced.

We present an up-to-date catalog of intrinsic iron abundance spreads in the 55 Milky Way globular clusters for which sufficiently precise spectroscopic measurements are available. Our method combines multiple datasets when possible to improve the statistics, taking into account the fact that different methods and instruments can lead to systematically offset metallicities. Only high spectral resolution (R>14,000) studies that measure the equivalent widths of individual iron lines are found to have uncertainties on the metallicities of the individual stars that can be calibrated sufficiently well for the intrinsic dispersion to be separated cleanly from random measurement error. The median intrinsic iron spread is found to be 0.045 dex, which is small but unambiguously measured to be non-zero in most cases. There is large variation between clusters, but more luminous globular clusters, above 10^5 L_sun, have increasingly large iron spreads on average; no trend between iron spread and metallicity is found.

We report on the spread of [Fe/H] values in the massive Large Magellanic Cloud cluster NGC1978, recently confirmed to harbor multiple populations of nearly the same age. We used accurate Stromgren photometry of carefully selected cluster red giant branch stars along with a high-dispersion spectroscopy-based calibration of the metallicity-sensitive index m1. Once accounting for the photometry quality, assessed from extensive artificial stars tests to trace the photometric uncertainties as a function of the position to the cluster's center as well as the stellar brightness, and those from the metallicity calibration, we found that NGC1978 exhibits a small metallicity spread of 0.035 dex (+- 0.019-0.023), depending on whether stars with individual sigma[Fe/H] < 0.15 dex or those located in the cluster's outer areas are considered. Such a spread in [Fe/H] is consistent with a cluster formation model with self-enrichment, if mass loss higher than 90 per cent due to stellar evolutionary and galactic tidal effects is assumed. Nevertheless, scenarios where the apparent [Fe/H] variation reflects CN abundance anomalies or less extreme mass loss models with environmentally-dependent self enrichment should not be ruled out.

Detailed observations of globular clusters (GCs) have revealed evidence of self-enrichment: some of the heavy elements that we see in stars today were produced by cluster stars themselves. Moreover, GCs have internal subpopulations with different elemental abundances, including, in some cases, in elements such as iron that are produced by supernovae. This paper presents a theoretical model for GC formation motivated by observations of Milky Way star forming regions and simulations of star formation, where giant molecular clouds fragment into multiple clumps which undergo star formation at slightly different times. Core collapse supernovae from earlier-forming clumps can enrich later-forming clumps to the degree that the ejecta can be retained within the gravitational potential well, resulting in subpopulations with different total metallicities once the clumps merge to form the final cluster. The model matches the mass-metallicity relation seen in GC populations around massive elliptical galaxies, and predicts metallicity spreads within clusters in excellent agreement with those seen in Milky Way GCs, even for those whose internal abundance spreads are so large that their entire identity as a GC is in question. The internal metallicity spread serves as an excellent measurement of how much self-enrichment has occurred in a cluster, a result that is very robust to variation in the model parameters.

We have explored the outskirts of dark matter haloes out to 2.5 times the virial radius using a large sample of halos drawn from Illustris, along with a set of zoom simulations (MUGS). Using these, we make a systematic exploration of the shape profile beyond R$_{vir}$. In the mean sphericity profile of Illustris halos we identify a dip close to the virial radius, which is robust across a broad range of masses and infall rates. The inner edge of this feature may be related to the virial radius and the outer edge with the splashback radius. Due to the high halo-to-halo variation this result is visible only on average. However, in four individual halos in the MUGS sample, a decrease in the sphericity and a subsequent recovery is evident close to the splashback radius. We find that this feature persists for several Gyr, growing with the halo. This feature appears at the interface between the spherical halo density distribution and the filamentary structure in the environment. The shape feature is strongest when there is a high rate of infall, implying that the effect is due to the mixing of accreting and virializing material. The filamentary velocity field becomes rapidly mixed in the halo region inside the virial radius, with the area between this and the splashback radius serving as the transition region. We also identify a long-lasting and smoothly evolving splashback region in the radial density gradient in many of the MUGS halos.

We present the discovery of d1005+68, a new faint dwarf galaxy in the M81 Group, using observations taken with the Subaru Hyper Suprime-Cam. d1005+68's color-magnitude diagram is consistent with a distance of $3.98_{-0.43}^{+0.39}$ Mpc, establishing group membership. We derive an absolute $V$-band magnitude, from stellar isochrone fitting, of $M_{V} = -7.94_{-0.50}^{+0.38}$, with a half-light radius of $r_{h} = 188_{-41}^{+39}$ pc. These place d1005+68 within the radius-luminosity locus of Local Group and M81 satellites and among the faintest confirmed satellites outside the Local Group. Assuming an age of 12 Gyr, d1005+68's red giant branch is best fit by an isochrone of [Fe/H] $= -1.90 \pm 0.24$. It has a projected separation from nearby M81 satellite BK5N of only 5 kpc. As this is well within BK5N's virial radius, we speculate that d1005+68 may be a satellite of BK5N. If confirmed, this would make d1005+68 one of the first detected satellites-of-a-satellite.

Galaxies with Milky Way-like stellar masses have a wide range of bulge and black hole masses; in turn, these correlate with other properties such as star formation history. While many processes may drive bulge formation, major and minor mergers are expected to play a crucial role. Stellar halos offer a novel and robust measurement of galactic merger history; cosmologically-motivated models predict that mergers with larger satellites produce more massive, higher metallicity stellar halos, reproducing the recently-observed stellar halo metallicity-mass relation. We quantify the relationship between stellar halo mass and bulge or black hole prominence using a sample of eighteen Milky Way-mass galaxies with newly-available measurements of (or limits on) stellar halo properties. There is an order of magnitude range in bulge mass, and two orders of magnitude in black hole mass, at a given stellar halo mass (or, equivalently, merger history). Galaxies with low mass bulges show a wide range of quiet merger histories, implying formation mechanisms that do not require intense merging activity. Galaxies with massive 'classical' bulges and central black holes also show a wide range of merger histories. While three of these galaxies have massive stellar halos consistent with a merger origin, two do not - merging appears to have had little impact in making these two massive 'classical' bulges. Such galaxies may be ideal laboratories to study massive bulge formation through pathways such as early gas-rich accretion, violent disk instabilities or misaligned infall of gas throughout cosmic time.

We present new deep photometry of the rich globular cluster (GC) systems around the Brightest Cluster Galaxies UGC 9799 (Abell 2052) and UGC 10143 (Abell 2147), obtained with the HST ACS and WFC3 cameras. For comparison, we also present new reductions of similar HST/ACS data for the Coma supergiants NGC 4874 and 4889. All four of these galaxies have huge cluster populations (to the radial limits of our data, comprising from 12000 to 23000 clusters per galaxy). The metallicity distribution functions (MDFs) of the GCs can still be matched by a bimodal-Gaussian form where the metal-rich and metal-poor modes are separated by ~0.8 dex, but the internal dispersions of each mode are so large that the total MDF becomes very broad and nearly continuous from [Fe/H] = -2.4 to Solar. There are, however, significant differences between galaxies in the relative numbers of \emph{metal-rich} clusters, suggesting that they underwent significantly different histories of mergers with massive, gas-rich halos. Lastly, the proportion of metal-poor GCs rises especially rapidly outside projected radii R > 4 R_eff, suggesting the importance of accreted dwarf satellites in the outer halo. Comprehensive models for the formation of GCs as part of the hierarchical formation of their parent galaxies will be needed to trace the systematic change in structure of the MDF with galaxy mass, from the distinctly bimodal form in smaller galaxies up to the broad continuum that we see in the very largest systems.

The number, distribution, and properties of dwarf satellites are crucial probes of the physics of galaxy formation at low masses and the response of satellite galaxies to the tidal and gas dynamical effects of their more massive parent.To make progress, it is necessary to augment and solidify the census of dwarf satellites of galaxies outside the Local Group. M\"uller et al. (2015) presented 16 dwarf galaxy candidates near M83, but lacking reliable distances, it is unclear which candidates are M83 satellites. Using red giant branch stars from the HST/GHOSTS survey in conjunction with ground-based images from VLT/VIMOS, we confirm that one of the candidates, dw1335-29-- with a projected distance of 26 kpc from M83 and a distance modulus of $(m - M)_0 = 28.5^{+0.3}_{-0.1}$ -- is a satellite of M83. We estimate an absolute magnitude $M_V = -10.1 \pm{0.4}$, an ellipticity of $0.40^{+0.14}_{-0.22}$, a half light radius of $656^{+121}_{-170 }$ pc, and [Fe/H] = $-1.3^{+0.3}_{-0.4}$. Owing to dw1335-29's somewhat irregular shape and possible young stars, we classify this galaxy as a dwarf irregular or transition dwarf. This is curious, as with a projected distance of 26 kpc from M83, dw1335-29 is expected to lack recent star formation. Further study of M83's dwarf population will reveal if star formation in its satellites is commonplace (suggesting a lack of a hot gas envelope for M83 that would quench star formation) or rare (suggesting that dw1335-29 has a larger M83-centric distance, and is fortuitously projected to small radii).

We have examined the resolved stellar populations at large galactocentric distances along the minor axis (from 10 kpc up to between 40 and 75 kpc), with limited major axis coverage, of six nearby highly-inclined Milky Way-mass disc galaxies using HST data from the GHOSTS survey. We select red giant branch stars to derive stellar halo density profiles. The projected minor axis density profiles can be approximated by power laws with projected slopes of between $-2$ and $-3.7$ and a diversity of stellar halo masses of $1-6\times 10^{9}M_{\odot}$, or $2-14\%$ of the total galaxy stellar masses. The typical intrinsic scatter around a smooth power law fit is $0.05-0.1$ dex owing to substructure. By comparing the minor and major axis profiles, we infer projected axis ratios $c/a$ at $\sim 25$ kpc between $0.4-0.75$. The GHOSTS stellar haloes are diverse, lying between the extremes charted out by the (rather atypical) haloes of the Milky Way and M31. We find a strong correlation between the stellar halo metallicities and the stellar halo masses. We compare our results with cosmological models, finding good agreement between our observations and accretion-only models where the stellar haloes are formed by the disruption of dwarf satellites. In particular, the strong observed correlation between stellar halo metallicity and mass is naturally reproduced. Low-resolution hydrodynamical models have unrealistically high stellar halo masses. Current high-resolution hydrodynamical models appear to predict stellar halo masses somewhat higher than observed but with reasonable metallicities, metallicity gradients and density profiles.

The spatial and velocity distributions of dark matter particles in the Milky Way Halo affect the signals expected to be observed in searches for dark matter. Results from direct detection experiments are often analyzed assuming a simple isothermal distribution of dark matter, the Standard Halo Model (SHM). Yet there has been skepticism regarding the validity of this simple model due to the complicated gravitational collapse and merger history of actual galaxies. In this paper we compare the SHM to the results of cosmological hydrodynamical simulations of galaxy formation to investigate whether or not the SHM is a good representation of the true WIMP distribution in the analysis of direct detection data. We examine two Milky Way-like galaxies from the MaGICC cosmological simulations (a) with dark matter only and (b) with baryonic physics included. The inclusion of baryons drives the shape of the DM halo to become more spherical and makes the velocity distribution of dark matter particles less anisotropic especially at large heliocentric velocities, thereby making the SHM a better fit. We also note that we do not find a significant disk-like rotating dark matter component in either of the two galaxy halos with baryons that we examine, suggesting that dark disks are not a generic prediction of cosmological hydrodynamical simulations. We conclude that in the Solar neighborhood, the SHM is in fact a good approximation to the true dark matter distribution in these cosmological simulations (with baryons) which are reasonable representations of the Milky Way, and hence can also be used for the purpose of dark matter direct detection calculations.

Measurement of the Galactic neutral atomic hydrogen (HI) column density, NHI, and brightness temperatures, Tb, is of high scientific value for a broad range of astrophysical disciplines. In the past two decades, one of the most-used legacy HI datasets has been the Leiden/Argentine/Bonn Survey (LAB). We release the HI 4$\pi$ survey (HI4PI), an all-sky database of Galactic HI, which supersedes the LAB survey. The HI4PI survey is based on data from the recently completed first coverage of the Effelsberg-Bonn HI Survey (EBHIS) and from the third revision of the Galactic All-Sky Survey (GASS). EBHIS and GASS share similar angular resolution and match well in sensitivity. Combined, they are ideally suited to be a successor to LAB. The new HI4PI survey outperforms the LAB in angular resolution (16.2', FWHM) and sensitivity (RMS: 43 mK). Moreover, it has full spatial sampling and thus overcomes a major drawback of LAB, which severely undersamples the sky. We publish all-sky column density maps of the neutral atomic hydrogen in the Milky Way, along with full spectroscopic data, in several map projections including HEALPix.

We study the stellar halo color properties of six nearby massive highly inclined disk galaxies using Hubble Space Telescope Advanced Camera for Surveys and Wide Field Camera 3 observations in both F606W and F814W filters from the GHOSTS survey. The observed fields, placed both along the minor and major axis of each galaxy, probe the stellar outskirts out to projected distances of ~ 50-70 kpc from their galactic centre along the minor axis. The 50% completeness levels of the color magnitude diagrams are typically at two mag below the tip of the red giant branch. We find that all galaxies have extended stellar halos out to ~ 50 kpc and two out to ~ 70 kpc. We determined the halo color distribution and color profile for each galaxy using the median colors of stars in the RGB. Within each galaxy we find variations in the median colors as a function of radius which likely indicates population variations, reflecting that their outskirts were built from several small accreted objects. We find that half of the galaxies (NGC 0891, NGC 4565, and NGC 7814) present a clear negative color gradient, reflecting a declining metallicity in their halos; the other have no significant color or population gradient. In addition, notwithstanding the modest sample size of galaxies, there is no strong correlation between their halo color/metallicity or gradient with galaxy's properties such as rotational velocity or stellar mass. The diversity in halo color profiles observed in the GHOSTS galaxies qualitatively supports the predicted galaxy-to-galaxy scatter in halo stellar properties; a consequence of the stochasticity inherent in the assembling history of galaxies.

Aims. Using a suite of cosmological chemodynamical disc galaxy simulations, we assess how (a) radial metallicity gradients evolve with scaleheight; (b) the vertical metallicity gradients change through the thick disc; and (c) the vertical gradient of the stellar rotation velocity varies through the disc. We compare with the Milky Way to search for analogous trends. Methods. We analyse five simulated spiral galaxies with masses comparable to the Milky Way. The simulations span a range of star formation and energy feedback strengths and prescriptions, particle- and grid-based hydrodynamical implementations, as well as initial conditions/assembly history. Results. Consistently, we find that the steeper, negative, radial metallicity gradients seen in the mid-plane flatten with increasing height away from the plane. In simulations with stronger (and/or more spatially-extended) feedback, the negative radial gradients invert, becoming positive for heights in excess of 1 kpc. Such behaviour is consistent with that inferred from recent observations. Our measurements of the vertical metallicity gradients show no clear correlation with galactocentric radius, and are in good agreement with those observed in the Milky Way's thick disc (locally). Conclusions. Simulations employing stronger/more extended feedback prescriptions possess radial and vertical metallicity and kinematic gradients more in line with recent observations. The inverted, positive, radial metallicity gradients seen in the simulated thick stellar discs originate from a population of younger, more metal-rich, stars formed in-situ, superimposed upon a background population of older migrators from the inner disc; the contrast provided by the former increases radially, due to the inside-out growth of the disc. A similar behaviour may be responsible for the same flattening seen in the radial gradients with scaleheight in the Milky Way.

We explore the chemical distribution of stars in a simulated galaxy. Using simulations of the same initial conditions but with two different feedback schemes (MUGS and MaGICC), we examine the features of the age-metallicity relation (AMR), and the three-dimensional age-metallicity-[O/Fe] distribution, both for the galaxy as a whole and decomposed into disc, bulge, halo, and satellites. The MUGS simulation, which uses traditional supernova feedback, is replete with chemical substructure. This sub- structure is absent from the MaGICC simulation, which includes early feedback from stellar winds, a modified IMF and more efficient feedback. The reduced amount of substructure is due to the almost complete lack of satellites in MaGICC. We identify a significant separation between the bulge and disc AMRs, where the bulge is considerably more metal-rich with a smaller spread in metallicity at any given time than the disc. Our results suggest, however, that identifying the substructure in observations will require exquisite age resolution, on the order of 0.25 Gyr. Certain satellites show exotic features in the AMR, even forming a 'sawtooth' shape of increasing metallicity followed by sharp declines which correspond to pericentric passages. This fact, along with the large spread in stellar age at a given metallicity, compromises the use of metallicity as an age indicator, although alpha abundance provides a more robust clock at early times. This may also impact algorithms that are used to reconstruct star formation histories from resolved stellar populations, which frequently assume a monotonically-increasing AMR.

We examine metallicities, ages and orbital properties of halo stars in a Milky-Way like disk galaxy formed in the cosmological hydrodynamical MaGICC simulations. Halo stars were either accreted from satellites or they formed in situ in the disk or bulge of the galaxy and were then kicked up into the halo ("in situ/ kicked-up" stars). Regardless of where they formed both types show surprisingly similar orbital properties: the majority of stars of both types are on short-axis tubes with the same sense of rotation as the disk -- implying that a large fraction of satellites are accreted onto the halo with the same sense of angular momentum as the disk.

We study the individual evolution histories of three nearby low-mass edge-on galaxies (IC 5052, NGC4244, and NGC5023). Using resolved stellar populations, we constructed star count density maps for populations of different ages and analyzed the change of structural parameters with stellar age within each galaxy. We do not detect a separate thick disk in any of the three galaxies, even though our observations cover a wider range in equivalent surface brightness than any integrated light study. While scale heights increase with age, each population can be well described by a single disk. Two of the galaxies contain a very weak additional component, which we identify as the faint halo. The mass of these faint halos is lower than 1% of the mass of the disk. The three galaxies show low vertical heating rates, which are much lower than the heating rate of the Milky Way. This indicates that heating agents, such as giant molecular clouds and spiral structure, are weak in low-mass galaxies. All populations in the three galaxies exhibit no or only little flaring. While this finding is consistent with previous integrated light studies, it poses strong constraints on galaxy simulations, where strong flaring is often found as a result of interactions or radial migration.

We present the first results from our HST Brightest Cluster Galaxy (BCG) survey of seven central supergiant cluster galaxies and their globular cluster (GC) systems. We measure a total of 48000 GCs in all seven galaxies, representing the largest single GC database. We find that a log-normal shape accurately matches the observed luminosity function (LF) of the GCs down to the GCLF turnover point, which is near our photometric limit. In addition, the LF has a virtually identical shape in all seven galaxies. Our data underscore the similarity in the formation mechanism of massive star clusters in diverse galactic environments. At the highest luminosities (log L > 10^7 L_Sun) we find small numbers of "superluminous" objects in five of the galaxies; their luminosity and color ranges are at least partly consistent with those of UCDs (Ultra-Compact Dwarfs). Lastly, we find preliminary evidence that in the outer halo (R > 20 kpc), the LF turnover point shows a weak dependence on projected distance, scaling as L_0 ~ R^-0.2, while the LF dispersion remains nearly constant.

We present HI spectral-line imaging of the extremely metal-poor galaxy DDO 68. This system has a nebular oxygen abundance of only 3% Z$_{\odot}$, making it one of the most metal-deficient galaxies known in the local volume. Surprisingly, DDO 68 is a relatively massive and luminous galaxy for its metal content, making it a significant outlier in the mass-metallicity and luminosity-metallicity relationships. The origin of such a low oxygen abundance in DDO 68 presents a challenge for models of the chemical evolution of galaxies. One possible solution to this problem is the infall of pristine neutral gas, potentially initiated during a gravitational interaction. Using archival HI spectral-line imaging obtained with the Karl G. Jansky Very Large Array, we have discovered a previously unknown companion of DDO 68. This low-mass (M$_{\rm HI}$ $=$ 2.8$\times$10$^{7}$ M$_{\odot}$), recently star-forming (SFR$_{\rm FUV}$ $=$ 1.4$\times$10$^{-3}$ M$_{\odot}$ yr$^{-1}$, SFR$_{\rm H\alpha}$ $<$ 7$\times$10$^{-5}$ M$_{\odot}$ yr$^{-1}$) companion has the same systemic velocity as DDO 68 (V$_{\rm sys}$ $=$ 506 km s$^{-1}$; D $=$ 12.74$\pm$0.27 Mpc) and is located at a projected distance of 42 kpc. New HI maps obtained with the 100m Robert C. Byrd Green Bank Telescope provide evidence that DDO 68 and this companion are gravitationally interacting at the present time. Low surface brightness HI gas forms a bridge between these objects.

We present an analysis of the role of feedback in shaping the neutral hydrogen (HI) content of simulated disc galaxies. For our analysis, we have used two realisations of two separate Milky Way-like (~L*) discs - one employing a conservative feedback scheme (MUGS), the other significantly more energetic (MaGICC). To quantify the impact of these schemes, we generate zeroth moment (surface density) maps of the inferred HI distribution; construct power spectra associated with the underlying structure of the simulated cold ISM, in addition to their radial surface density and velocity dispersion profiles. Our results are compared with a parallel, self-consistent, analysis of empirical data from THINGS (The HI Nearby Galaxy Survey). Single power-law fits (P~k^gamma) to the power spectra of the stronger-feedback (MaGICC) runs (over spatial scales corresponding to 0.5 kpc to 20 kpc) result in slopes consistent with those seen in the THINGS sample (gamma = -2.5). The weaker-feedback (MUGS) runs exhibit shallower power law slopes (gamma = -1.2). The power spectra of the MaGICC simulations are more consistent though with a two-component fit, with a flatter distribution of power on larger scales (i.e., gamma = -1.4 for scales in excess of 2 kpc) and a steeper slope on scales below 1 kpc (gamma = -5), qualitatively consistent with empirical claims, as well as our earlier work on dwarf discs. The radial HI surface density profiles of the MaGICC discs show a clear exponential behaviour, while those of the MUGS suite are essentially flat; both behaviours are encountered in nature, although the THINGS sample is more consistent with our stronger (MaGICC) feedback runs.

The limits of available computing power have forced models for the structure of stellar halos to adopt one or both of the following simplifying assumptions: (1) stellar mass can be "painted" onto dark matter particles in progenitor satellites; (2) pure dark matter simulations that do not form a luminous galaxy can be used. We estimate the magnitude of the systematic errors introduced by these assumptions using a controlled set of stellar halo models where we independently vary whether we look at star particles or painted dark matter particles, and whether we use a simulation in which a baryonic disk galaxy forms or a matching pure dark matter simulation that does not form a baryonic disk. We find that the "painting" simplification reduces the halo concentration and internal structure, predominantly because painted dark matter particles have different kinematics than star particles even when both are buried deep in the potential well of the satellite. The simplification of using pure dark matter simulations reduces the concentration further, but increases the internal structure, and results in a more prolate stellar halo. These differences can be a factor of 1.5-7 in concentration (as measured by the half-mass radius) and 2-7 in internal density structure. Given this level of systematic uncertainty, one should be wary of overinterpreting differences between observations and the current generation of stellar halo models based on dark matter only simulations when such differences are less than an order of magnitude.

The determination of stellar metallicity and its gradient in external galaxies is a difficult task, but crucial for the understanding of galaxy formation and evolution. The color of the Red Giant Branch (RGB) can be used to determine metallicities of stellar populations that have only shallow photometry. We will quantify the relation between metallicity and color in the widely used HST ACS filters F606W and F814W. We use a sample of globular clusters from the ACS Globular Cluster Survey and measure their RGB color at given absolute magnitudes to derive the color-metallicity relation. We especially investigate the scatter and the uncertainties in this relation and show its limitations. There is a clear relation between metallicity and RGB color. A comparison with isochrones shows reasonably good agreement with BaSTI models, a small offset to Dartmouth models, and a larger offset to Padua models. Even for the best globular cluster data available, the metallicity of a simple stellar population can be determined from the RGB alone only with an accuracy of 0.3dex for [M/H]< -1, and 0.15dex for [M/H]> -1. For mixed populations, as they are observed in external galaxies, the uncertainties will be even larger due to uncertainties in extinction, age, etc. Therefore caution is necessary when interpreting photometric metallicities.

We have mapped the distribution of young and old stars in the gaseous HI warp of NGC 4565. We find a clear correlation of young stars (<600 Myr) with the warp, but no coincident old stars (>1 Gyr), which places an upper limit on the age of the structure. The formation rate of the young stars, which increased ~300 Myr ago relative to the surrounding regions, is (6.3 +2.5/-1.5) x 10^-5 M_sol/yr/kpc^2. This implies a ~60+/-20 Gyr depletion time of the HI warp, similar to the timescales calculated for the outer HI disks of nearby spiral galaxies. While some stars associated with the warp fall into the asymptotic giant branch (AGB) region of the color magnitude diagram, where stars could be as old as 1 Gyr, further investigation suggests that they may be interlopers rather than real AGB stars. We discuss the implications of these age constraints for the formation of HI warps, and the gas fueling of disk galaxies.

We report the discovery of a new faint dwarf galaxy, GHOSTS I, using HST/ACS data from one of our GHOSTS (Galaxy Halos, Outer disks, Substructure, Thick disk, and Star clusters) fields. Its detected individual stars populate an approximately one magnitude range of its luminosity function (LF). Using synthetic color-magnitude diagrams (CMDs) to compare with the galaxy's CMD, we find that the colors and magnitudes of GHOSTS I's individual stars are most consistent with being young helium-burning and asymptotic giant branch stars at a distance of 12 +/- 2 Mpc. Morphologically, GHOSTS I appears to be actively forming stars, so we tentatively classify it as a dwarf irregular (dIrr) galaxy, although future HST observations deep enough to resolve a larger magnitude range in its LF are required to make a more secure classification. GHOSTS I's absolute magnitude is $M_V = -9.85^{+ 0.40}_{- 0.33}$, making it one of the least luminous dIrr galaxies known, and its metallicity is lower than [Fe/H] =-1.5 dex. The half-light radius of GHOSTS I is 226 +/- 38 pc and its ellipticity is 0.47 +/- 0.07, similar to Milky Way and M31 dwarf satellites at comparable luminosity. There are no luminous massive galaxies or galaxy clusters within ~ 4 Mpc from GHOSTS I that could be considered as its host, making it a very isolated dwarf galaxy in the Local Universe.

We examine the role of energy feedback in shaping the distribution of metals within cosmological hydrodynamical simulations of L* disc galaxies. While negative abundance gradients today provide a boundary condition for galaxy evolution models, in support of inside-out disc growth, empirical evidence as to whether abundance gradients steepen or flatten with time remains highly contradictory. We made use of a suite of L* discs, realised with and without `enhanced' feedback. All the simulations were produced using the smoothed particle hydrodynamics code Gasoline, and their in situ gas-phase metallicity gradients traced from redshift z~2 to the present-day. Present-day age-metallicity relations and metallicity distribution functions were derived for each system. The `enhanced' feedback models, which have been shown to be in agreement with a broad range of empirical scaling relations, distribute energy and re-cycled ISM material over large scales and predict the existence of relatively `flat' and temporally invariant abundance gradients. Enhanced feedback schemes reduce significantly the scatter in the local stellar age-metallicity relation and, especially, the [O/Fe]-[Fe/H] relation. The local [O/Fe] distribution functions for our L* discs show clear bimodality, with peaks at [O/Fe]=-0.05 and +0.05 (for stars with [Fe/H]>-1), consistent with our earlier work on dwarf discs. Our results with `enhanced' feedback are inconsistent with our earlier generation of simulations realised with `conservative' feedback. We conclude that spatially-resolved metallicity distributions, particularly at high-redshift, offer a unique and under-utilised constraint on the uncertain nature of stellar feedback processes.

We analyze the orbits of stars and dark matter particles in the halo of a disk galaxy formed in a cosmological hydrodynamical simulation. The halo is oblate within the inner ~20 kpc and triaxial beyond this radius. About 43% of orbits are short axis tubes - the rest belong to orbit families that characterize triaxial potentials (boxes, long-axis tubes and chaotic orbits), but their shapes are close to axisymmetric. We find no evidence that the self-consistent distribution function of the nearly oblate inner halo is comprised primarily of axisymmetric short-axis tube orbits. Orbits of all families, and both types of particles are highly eccentric with mean eccentricity >0.6. We find that randomly selected samples of halo stars show no substructure in "integrals of motion" space. However individual accretion events can be clearly identified in plots of metallicity versus formation time. Dynamically young tidal debris is found primarily on a single type of orbit. However, stars associated with older satellites become chaotically mixed during the formation process (possibly due to scattering by the central bulge and disk, and baryonic processes), and appear on all four types of orbits. We find that the tidal debris in cosmological hydrodynamical simulations experiences significantly more chaotic evolution than in collisionless simulations, making it much harder to identify individual progenitors using phase space coordinates alone. However by combining information on stellar ages and chemical abundances with the orbital properties of halo stars in the underlying self-consistent potential, the identification of progenitors is likely to be possible.

We study the properties of the stellar populations in M81's outermost part, which hereafter we will term the stellar halo, using HST ACS/WFC observations of 19 fields from the GHOSTS survey. The observed fields probe the stellar halo out to a projected distance of ~ 50 kpc from the galactic center. Each field was observed in both F606W and F814W filters. The 50% completeness levels of the color magnitude diagrams (CMDs) are typically at 2 mag below the tip of the red giant branch (TRGB). Fields at distances closer than 15 kpc show evidence of disk-dominated populations whereas fields at larger distances are mostly populated by halo stars. The RGB of the M81's halo CMDs is well matched with isochrones of ~ 10 Gyr and metallicities [Fe/H] ~ -1.2 dex, suggesting that the dominant stellar population of M81's halo has a similar age and metallicity. The halo of M81 is characterized by a color distribution of width ~ 0.4 mag and an approximately constant median value of (F606W - F814W) ~ 1 mag measured using stars within the magnitude range 23.7 < F814W < 25.5. When considering only fields located at galactocentric radius R > 15 kpc, we detect no color gradient in the stellar halo of M81. We place a limit of 0.03+/-0.11 mag difference between the median color of RGB M81 halo stars at ~ 15 and at 50 kpc, corresponding to a metallicity difference of 0.08+/-0.35 dex over that radial range for an assumed constant age of 10 Gyr. We compare these results with model predictions for the colors of stellar halos formed purely via accretion of satellite galaxies. When we analyze the cosmologically motivated models in the same way as the HST data, we find that they predict no color gradient for the stellar halos, in good agreement with the observations.

The high-resolution, SPH galaxies of the McMaster Unbiased Galaxy Survey (MUGS) are used to examine the satellite systems of sixteen model host galaxies. Each galaxy has a different mass, angular momentum and merger history that yield a rich set of satellite luminosity functions. With new observations of distant satellite systems, we can compare these luminosity functions to satellite systems beyond the Local Group. We find that the luminosity functions of our simulations compare well to observations when the luminosity functions are scaled according to host mass. We use the recently-found relationship between dwarf satellites and host mass in distant satellite systems (Trentham & Tully 2009) to normalize a theoretical, complete luminosity function for the Milky Way (Koposov et al. 2008). The luminosity function of satellites, expressed as a function of the host mass, is given by dN/dM_V= 3.5 x M_host^(0.91) x 10^(0.1 x M_V-10.2), where host mass is given in M_sun. The mass of a host galaxy can be used to predict the number of dwarf satellites and even when considering spiral and elliptical hosts separately this relation holds.

By means of high-resolution cosmological hydrodynamical simulations of Milky Way-like disc galaxies, we conduct an analysis of the associated stellar metallicity distribution functions (MDFs). After undertaking a kinematic decomposition of each simulation into spheroid and disc sub-components, we compare the predicted MDFs to those observed in the solar neighbourhood and the Galactic bulge. The effects of the star formation density threshold are visible in the star formation histories, which show a modulation in their behaviour driven by the threshold. The derived MDFs show median metallicities lower by 0.2-0.3 dex than the MDF observed locally in the disc and in the Galactic bulge. Possible reasons for this apparent discrepancy include the use of low stellar yields and/or centrally-concentrated star formation. The dispersions are larger than the one of the observed MDF; this could be due to simulated discs being kinematically hotter relative to the Milky Way. The fraction of low metallicity stars is largely overestimated, visible from the more negatively skewed MDF with respect to the observational sample. For our fiducial Milky Way analog, we study the metallicity distribution of the stars born "in situ" relative to those formed via accretion (from disrupted satellites), and demonstrate that this low-metallicity tail to the MDF is populated primarily by accreted stars. Enhanced supernova and stellar radiation energy feedback to the surrounding interstellar media of these pre-disrupted satellites is suggested as an important regulator of the MDF skewness.

We examine the chemical properties of 5 cosmological hydrodynamical simulations of an M33-like disc galaxy which have been shown to be consistent with the morphological characteristics and bulk scaling relations expected of late-type spirals. These simulations are part of the Making Galaxies In a Cosmological Context (MaGICC) Project, in which stellar feedback is tuned to match the stellar mass -- halo mass relationship. Each realisation employed identical initial conditions and assembly histories, but differed from one another in their underlying baryonic physics prescriptions, including (a) the efficiency with which each supernova energy couples to the ISM, (b) the impact of feedback associated with massive star radiation pressure, (c) the role of the minimum shut-off time for radiative cooling of Type II SNe remnants, (d) the treatment of metal diffusion, and (e) varying the IMF. Our analysis focusses on the resulting stellar metallicity distribution functions (MDFs) in each simulated (analogous) `solar neighbourhood' and central `bulge' region. We compare the simulated MDFs' skewness, kurtosis, and dispersion (inter-quartile, inter-decile, inter-centile, and inter-tenth-percentile regions) with that of the empirical solar neighbourhood MDF and Local Group dwarfs. We find that the MDFs of the simulated discs are more negatively skewed, with higher kurtosis, than those observed locally. We can trace this difference to the simulations' tight and correlated age-metallicity relations (compared with that of the Milky Way), suggesting that these relations within `dwarf' discs might be steeper than in L* discs and/or the degree of stellar orbital re-distribution and migration inferred locally has not been captured in their entirety, at the resolution of our simulations. The important role of metal diffusion in ameliorating the over-production of extremely metal-poor stars is highlighted.

We examine radial and vertical metallicity gradients using a suite of disk galaxy simulations, supplemented with two classic chemical evolution approaches. We determine the rate of change of gradient and reconcile differences between extant models and observations within the `inside-out' disk growth paradigm. A sample of 25 disks is used, consisting of 19 from our RaDES (Ramses Disk Environment Study) sample, realised with the adaptive mesh refinement code RAMSES. Four disks are selected from the MUGS (McMaster Unbiased Galaxy Simulations) sample, generated with the smoothed particle hydrodynamics (SPH) code GASOLINE, alongside disks from Rahimi et al. (GCD+) and Kobayashi & Nakasato (GRAPE-SPH). Two chemical evolution models of inside-out disk growth were employed to contrast the temporal evolution of their radial gradients with those of the simulations. We find that systematic differences exist between the predicted evolution of radial abundance gradients in the RaDES and chemical evolution models, compared with the MUGS sample; specifically, the MUGS simulations are systematically steeper at high-redshift, and present much more rapid evolution in their gradients. We find that the majority of the models predict radial gradients today which are consistent with those observed in late-type disks, but they evolve to this self-similarity in different fashions, despite each adhering to classical `inside-out' growth. We find that radial dependence of the efficiency with which stars form as a function of time drives the differences seen in the gradients; systematic differences in the sub-grid physics between the various codes are responsible for setting these gradients. Recent, albeit limited, data at redshift z=1.5 are consistent with the steeper gradients seen in our SPH sample, suggesting a modest revision of the classical chemical evolution models may be required.

Encoded within the morphological structure of galaxies are clues related to their formation and evolutionary history. Recent advances pertaining to the statistics of galaxy morphology include sophisticated measures of concentration (C), asymmetry (A), and clumpiness (S). In this study, these three parameters (CAS) have been applied to a suite of simulated galaxies and compared with observational results inferred from a sample of nearby galaxies. The simulations span a range of late-type systems, with masses between ~1e10 Msun and ~1e12 Msun, and employ star formation density thresholds between 0.1 cm^-3 and 100 cm^-3. We have found that the simulated galaxies possess comparable concentrations to their real counterparts. However, the results of the CAS analysis revealed that the simulated galaxies are generally more asymmetric, and that the range of clumpiness values extends beyond the range of those observed. Strong correlations were obtained between the three CAS parameters and colour (B-V), consistent with observed galaxies. Furthermore, the simulated galaxies possess strong links between their CAS parameters and Hubble type, mostly in-line with their real counterparts.

We have obtained Magellan/IMACS and HST/ACS imaging data that resolve red giant branch stars in the stellar halo of the starburst galaxy NGC 253. The HST data cover a small area, and allow us to accurately interpret the ground-based data, which cover 30% of the halo to a distance of 30 kpc, allowing us to make detailed quantitative measurements of the global properties and structure of a stellar halo outside of the Local Group. The geometry of the halo is significantly flattened in the same sense as the disk, with a projected axis ratio of b/a ~ 0.35 +/- 0.1. The total stellar mass of the halo is estimated to be M_halo ~ 2.5 +/- 1.5 x 10^9 M_sun, or 6% of the total stellar mass of the galaxy, and has a projected radial dependence that follows a power law of index -2.8 +/- 0.6, corresponding to a three-dimensional power law index of ~ -4. The total luminosity and profile shape that we measure for NGC 253 are somewhat larger and steeper than the equivalent values for the Milky Way and M31, but are well within the scatter of model predictions for the properties of stellar halos built up in a cosmological context. Structure within the halo is seen at a variety of scales: there is small kpc-scale density variation and a large shelf-like feature near the middle of the field. The techniques that have been developed will be essential for quantitatively comparing our upcoming larger sample of observed stellar halos to models of halo formation.

We present a study of satellites in orbit around a high-resolution, smoothed particle hydrodynamics (SPH) galaxy simulated in a cosmological context. The simulated galaxy is approximately the same mass as the Milky Way. The cumulative number of luminous satellites at z = 0 is similar to the observed system of satellites orbiting the Milky Way although an analysis of the satellite mass function reveals an order of magnitude more dark satellites than luminous. Some of the dark subhalos are more massive than some of the luminous subhalos at z = 0. What separates luminous and dark subhalos is not their mass at z = 0, but the maximum mass the subhalos ever achieve. We study the effect of four mass loss mechanisms on the subhalos: ultraviolet (UV) ionising radiation, ram pressure stripping, tidal stripping, and stellar feedback, and compare the impact of each of these four mechanisms on the satellites. In the lowest mass subhalos, UV is responsible for the majority of the baryonic mass loss. Ram pressure stripping removes whatever mass remains from the low mass satellites. More massive subhalos have deeper potential wells and retain more mass during reionisation. However, as satellites pass near the centre of the main halo, tidal forces cause significant mass loss from satellites of all masses. Satellites that are tidally stripped from the outside can account for the luminous satellites that are lower mass than some of the dark satellites. Stellar feedback has the greatest impact on medium mass satellites that had formed stars, but lost all their gas by z = 0. Our results demonstrate that the missing satellite problem is not an intractable issue with the cold dark matter cosmology, but is rather a manifestation of baryonic processes.

We present the McMaster Unbiased Galaxy Simulations (MUGS), the first 9 galaxies of an unbiased selection ranging in total mass from 5$\times10^{11}$ M$_\odot$ to 2$\times10^{12}$ M$_\odot$ simulated using n-body smoothed particle hydrodynamics (SPH) at high resolution. The simulations include a treatment of low temperature metal cooling, UV background radiation, star formation, and physically motivated stellar feedback. Mock images of the simulations show that the simulations lie within the observed range of relations such as that between color and magnitude and that between brightness and circular velocity (Tully-Fisher). The greatest discrepancy between the simulated galaxies and observed galaxies is the high concentration of material at the center of the galaxies as represented by the centrally peaked rotation curves and the high bulge-to-total ratios of the simulations determined both kinematically and photometrically. This central concentration represents the excess of low angular momentum material that long has plagued morphological studies of simulated galaxies and suggests that higher resolutions and a more accurate description of feedback will be required to simulate more realistic galaxies. Even with the excess central mass concentrations, the simulations suggest the important role merger history and halo spin play in the formation of disks.

It has been hypothesized that the cosmic microwave background (CMB) provides a temperature floor for collapsing protostars that can regulate the process of star formation and result in a top-heavy initial mass function at high metallicity and high redshift. We examine whether this hypothesis has any testable observational consequences. First we determine, using a set of hydrodynamic galaxy formation simulations, that the CMB temperature floor would have influenced the majority of stars formed at redshifts between z=3 and 6, and probably even to higher redshift. Five signatures of CMB-regulated star formation are: (1) a higher supernova rate than currently predicted at high redshift; (2) a systematic discrepancy between direct and indirect measurements of the high redshift star formation rate; (3) a lack of surviving globular clusters that formed at high metallicity and high redshift; (4) a more rapid rise in the metallicity of cosmic gas than is predicted by current simulations; and (5) an enhancement in the abundances of alpha elements such as O and Mg at metallicities -2 < [Fe/H] < -0.5. Observations are not presently able to either confirm or rule out the presence of these signatures. However, if correct, the top-heavy IMF of high-redshift high-metallicity globular clusters could provide an explanation for the observed bimodality of their metallicity distribution.

Here we present HST/ACS imaging, in the B and I bands, of the edge-on Sb/Sc galaxy NGC 5170. Excluding the central disk region region, we detect a 142 objects with colours and sizes typical of globular clusters (GCs). Our main result is the discovery of a `blue tilt' (a mass-metallicity relation), at the 3sigma level, in the metal-poor GC subpopulation of this Milky Way like galaxy. The tilt is consistent with that seen in massive elliptical galaxies and with the self enrichment model of Bailin & Harris. For a linear mass-metallicity relation, the tilt has the form Z ~ L^{0.42 +/- 0.13}. We derive a total GC system population of 600 +/- 100, making it much richer than the Milky Way. However when this number is normalised by the host galaxy luminosity or stellar mass it is similar to that of M31. Finally, we report the presence of a potential Ultra Compact Dwarf of size ~ 6 pc and luminosity M_I ~ -12.5, assuming it is physically associated with NGC 5170.

Images from the HST ACS are used to carry out a new photometric study of the globular clusters (GCs) in M104, the Sombrero galaxy. The primary focus of our study is the characteristic distribution function of linear sizes (SDF) of the GCs. We measure the effective radii for 652 clusters with PSF-convolved King and Wilson dynamical model fits. The SDF is remarkably similar to those measured for other large galaxies of all types, adding strong support to the view that it is a "universal" feature of globular cluster systems. We develop a more general interpretation of the size distribution function for globular clusters, proposing that the shape of the SDF that we see today for GCs is strongly influenced by the early rapid mass loss during their star forming stage, coupled with stochastic differences from cluster to cluster in the star formation efficiency (SFE) and their initial sizes. We find that the observed SDF shape can be accurately predicted by a simple model in which the protocluster clouds had characteristic sizes of $0.9 \pm 0.1$ pc and SFEs of $0.3 \pm 0.07$. The colors and luminosities of the M104 clusters show the clearly defined classic bimodal form. The blue sequence exhibits a mass/metallicity relation (MMR), following a scaling of heavy-element abundance with luminosity of $Z \sim L^{0.3}$ very similar to what has been found in most giant elliptical galaxies. A quantitative self-enrichment model provides a good first-order match to the data for the same initial SFE and protocluster size that were required to explain the SDF. We also discuss various forms of the globular cluster Fundamental Plane (FP) of structural parameters, and show that useful tests of it can be extended to galaxies beyond the Local Group.

We develop a model for stochastic pre-enrichment and self-enrichment in globular clusters (GCs) during their formation process. GCs beginning their formation have an initial metallicity determined by the pre-enrichment of their surrounding protocloud, but can also undergo internal self-enrichment during formation. Stochastic variations in metallicity arise because of the finite numbers of supernova. We construct an analytic formulation of the combined effects of pre-enrichment and self-enrichment and use Monte Carlo models to verify that the model accurately encapsulates the mean metallicity and metallicity spread among real GCs. The predicted metallicity spread due to self-enrichment alone, a robust prediction of the model, is much smaller than the observed spread among real GCs. This result rules out self-enrichment as a significant contributor to the metal content in most GCs, leaving pre-enrichment as the viable alternative. Self-enrichment can, however, be important for clusters with masses well above 10^6 Msun, which are massive enough to hold in a significant fraction of their SN ejecta even without any external pressure confinement. This transition point corresponds well to the mass at which a mass-metallicity relationship ("blue tilt") appears in the metal-poor cluster sequence in many large galaxies. We therefore suggest that self-enrichment is the primary driver for the mass-metallicity relation. Other predictions from our model are that the cluster-to-cluster metallicity spread decreases amongst the highest mass clusters; and that the red GC sequence should also display a more modest mass-metallicity trend if it can be traced to similarly high mass.

The Parkes Galactic All-Sky Survey (GASS) is a survey of Galactic atomic hydrogen (HI) emission in the Southern sky covering declinations $\delta \leq 1^{\circ}$ using the Parkes Radio Telescope. The survey covers $2\pi$ steradians with an effective angular resolution of ~16', at a velocity resolution of 1.0 km/s, and with an rms brightness temperature noise of 57 mK. GASS is the most sensitive, highest angular resolution survey of Galactic HI emission ever made in the Southern sky. In this paper we outline the survey goals, describe the observations and data analysis, and present the first-stage data release. The data product is a single cube at full resolution, not corrected for stray radiation. Spectra from the survey and other data products are publicly available online.

We identify satellites of isolated galaxies in SDSS and examine their angular distribution. Using mock catalogues generated from cosmological N-body simulations, we demonstrate that the selection criteria used to select isolated galaxies and their satellites must be very strict in order to correctly identify systems in which the primary galaxy dominates its environment. The criteria used in many previous studies instead select predominantly group members. We refine a set of selection criteria for which the group contamination is estimated to be less than 7% and present a catalogue of the resulting sample. The angular distribution of satellites about their host is biased towards the major axes for spheroidal galaxies and probably also for red disc galaxies, but is isotropic for blue disc galaxies, i.e. it is the colour of the host that determines the distribution of its satellites rather than its morphology. The similar anisotropy measured in this study as in studies that were dominated by groups implies that group-specific processes are not responsible for the angular distribution. Satellites that are most likely to have been recently accreted show a tendancy to lie along the same axis as the surrounding large scale structure. The orientations of isolated early and intermediate-type galaxies also align with the surrounding large scale structures. We discuss the origin of the anisotropic satellite distribution and consider the implications of our results, critically assessing the respective roles played by the orientation of the visible galaxy within its dark matter halo; anisotropic accretion of satellites from the larger scale environment; and the biased nature of satellites as tracers of the underlying dark matter subhalo population. (Abridged)

We have detected over 400 HI clouds in the lower halo of the Galaxy within the pilot region of the Galactic All-Sky Survey (GASS), a region of the fourth quadrant that spans 18 degrees in longitude, 40 degrees in latitude and is centered on the Galactic equator. These clouds have a median peak brightness temperature of 0.6 K, a median velocity width of 12.8 km/s, and angular sizes <1 degree. The motion of these clouds is dominated by Galactic rotation with a random cloud-to-cloud velocity dispersion of 18 km/s. A sample of clouds likely to be near tangent points was analyzed in detail. These clouds have radii on the order of 30 pc and a median HI mass of 630 Msun. The population has a vertical scale height of 400 pc and is concentrated in Galactocentric radius, peaking at R=3.8 kpc. This confined structure suggests that the clouds are linked to spiral features, while morphological evidence that many clouds are aligned with loops and filaments is suggestive of a relationship with star formation. The clouds might result from supernovae and stellar winds in the form of fragmenting shells and gas that has been pushed into the halo rather than from a galactic fountain.

We present the second data release of the Radial Velocity Experiment (RAVE), an ambitious spectroscopic survey to measure radial velocities (RVs) and stellar atmosphere parameters of up to one million stars using the 6dF multi-object spectrograph on the 1.2-m UK Schmidt Telescope of the Anglo-Australian Observatory (AAO). It is obtaining medium resolution spectra (median R=7,500) in the Ca-triplet region (8,410--8,795 \AA) for southern hemisphere stars in the magnitude range 9<I<12. Following the first data release (Steinmetz et al. 2006) the current release doubles the sample of published RVs, now containing 51,829 RVs for 49,327 individual stars observed on 141 nights between April 11 2003 and March 31 2005. Comparison with external data sets shows that the new data collected since April 3 2004 show a standard deviation of 1.3 km/s, about twice better than for the first data release. For the first time this data release contains values of stellar parameters from 22,407 spectra of 21,121 individual stars. They were derived by a penalized \chi^2 method using an extensive grid of synthetic spectra calculated from the latest version of Kurucz models. From comparison with external data sets, our conservative estimates of errors of the stellar parameters (for a spectrum with S/N=40) are 400 K in temperature, 0.5 dex in gravity, and 0.2 dex in metallicity. We note however that the internal errors estimated from repeat RAVE observations of 822 stars are at least a factor 2 smaller. We demonstrate that the results show no systematic offsets if compared to values derived from photometry or complementary spectroscopic analyses. The data release includes proper motion and photometric measurements. It can be accessed via the RAVE webpage: http://www.rave-survey.org and through CDS.

We present a new method to classify galaxies from large surveys like the Sloan Digital Sky Survey using inclination-corrected concentration, inclination-corrected location on the color-magnitude diagram, and apparent axis ratio. Explicitly accounting for inclination tightens the distribution of each of these parameters and enables simple boundaries to be drawn that delineate three different galaxy populations: Early-type galaxies, which are red, highly concentrated, and round; Late-type galaxies, which are blue, have low concentrations, and are disk dominated; and Intermediate-type galaxies, which are red, have intermediate concentrations, and have disks. We have validated our method by comparing to visual classifications of high-quality imaging data from the Millennium Galaxy Catalogue. The inclination correction is crucial to unveiling the previously unrecognized Intermediate class. Intermediate-type galaxies, roughly corresponding to lenticulars and early spirals, lie on the red sequence. The red sequence is therefore composed of two distinct morphological types, suggesting that there are two distinct mechanisms for transiting to the red sequence. We propose that Intermediate-type galaxies are those that have lost their cold gas via strangulation, while Early-type galaxies are those that have experienced a major merger that either consumed their cold gas, or whose merger progenitors were already devoid of cold gas (the ``dry merger'' scenario).

We have analysed the distribution of inclination-corrected galaxy concentrations in the Sloan Digital Sky Survey. We find that unlike most galaxy properties, which are distributed bimodally, the distribution of concentrations is trimodal: it exhibits three distinct peaks. The newly-discovered intermediate peak, which consists of early-type spirals and lenticulars, may contain ~60% of the number density and ~50% of the luminosity density of M_r < -17 galaxies in the local universe. These galaxies are generally red and quiescent, although the distribution contains a tail of blue star-forming galaxies and also shows evidence of dust. The intermediate-type galaxies have higher apparent ellipticities than either disc or elliptical galaxies, most likely because some of the face-on intermediate types are misidentified as ellipticals. Their physical half-light radii are smaller than the radii of either the disc or elliptical galaxies, which may be evidence that they form from disc fading. The existence of a distinct peak in parameter space associated with early-type spiral galaxies and lenticulars implies that they have a distinct formation mechanism and are not simply the smooth transition between disc-dominated and spheroid-dominated galaxies.

We measure the distribution of velocities for prograde and retrograde satellite galaxies using a combination of published data and new observations for 78 satellites of 63 extremely isolated disc galaxies (169 satellites total). We find that the velocity distribution is non-Gaussian (>99.9% confidence), but that it can be described as the sum of two Gaussians, one of which is broad (sigma = 176 \pm 15 km/s), has a mean prograde velocity of 86 \pm 30 km/s, and contains ~55% of the satellites, while the other is slightly retrograde with a mean velocity of -21 \pm 22 km/s and sigma = 74 \pm 18 km/s and contains ~45% of the satellites. Both of these components are present over all projected radii and found in the sample regardless of cuts on primary inclination or satellite disc angle. The double-Gaussian shape, however, becomes more pronounced among satellites of more luminous primaries. We remove the potential dependence of satellite velocity on primary luminosity using the Tully-Fisher relation and still find the velocity distribution to be asymmetric and even more significantly non-Gaussian. The asymmetric velocity distribution demonstrates a connection between the inner, visible disc galaxy and the kinematics of the outer, dark halo. The reach of this connection, extending even beyond the virial radii, suggests that it is imprinted by the satellite infall pattern and large-scale effects, rather than by higher-level dynamical processes in the formation of the central galaxy or late-term evolution of the satellites.

Galactic disks in triaxial dark matter halos become deformed by the elliptical potential in the plane of the disk in such a way as to counteract the halo ellipticity. We develop a technique to calculate the equilibrium configuration of such a disk in the combined disk-halo potential, which is based on the method of Jog (2000) but accounts for the radial variation in both the halo potential and the disk ellipticity. This crucial ingredient results in qualitatively different behavior of the disk: the disk circularizes the potential at small radii, even for a reasonably low disk mass. This effect has important implications for proposals to reconcile cuspy halo density profiles with low surface brightness galaxy rotation curves using halo triaxiality. The disk ellipticities in our models are consistent with observational estimates based on two-dimensional velocity fields and isophotal axis ratios.

Neutral hydrogen observations towards the Bootes dwarf spheroidal galaxy, a very low luminosity metal-poor Galactic satellite, were obtained using the Parkes Radio Telescope. We do not detect any HI in or around Bootes to a 3sigma upper limit of 180 Msun within the optical half light radius and 8000 Msun within 1.6 kpc. Its HI mass-to-light ratio is less than 0.002 Msun/Lsun, making Bootes one of the most gas-poor galaxies known. Either reionisation severely inhibited gas infall onto the proto-Bootes, or large amounts of gas have been removed by ram pressure and/or tidal stripping. Since Bootes lies on the mass-metallicity fundamental line, this relation and the inefficiency of star formation at the faintest end of the galaxy luminosity function must be partly driven, or at least not disrupted, by extreme gas loss in such low luminosity galaxies. We also do not detect any HI associated with the leading tidal tail of the Sagittarius dSph galaxy, which fortuitously passes through the observed field, to a 3sigma column density limit of 2 x 10^17 cm^-2. This suggests that either the leading gaseous tail is ionised, or the gas in the trailing tail was removed before the current tidal disruption of the parent dSph began.

We report that neutral hydrogen (HI) gas clouds, resembling High Velocity Clouds (HVCs) observed in the Milky Way (MW), appear in MW-sized disk galaxies formed in high-resolution Lambda Cold Dark Matter (LCDM) cosmological simulations which include gas-dynamics, radiative cooling, star formation, supernova feedback, and metal enrichment. Two such disk galaxies are analyzed, and HI column density and velocity distributions in all-sky Aitoff projections are constructed. The simulations demonstrate that LCDM is able to create galaxies with sufficient numbers of anomalous velocity gas clouds consistent with the HVCs observed within the MW, and that they are found within a galactocentric radius of 150 kpc. We also find that one of the galaxies has a polar gas ring, with radius 30 kpc, which appears as a large structure of HVCs in the Aitoff projection. Such large structures may share an origin similar to extended HVCs observed in the MW, such as Complex C.

We present the first data release of the Radial Velocity Experiment (RAVE), an ambitious spectroscopic survey to measure radial velocities and stellar atmosphere parameters (temperature, metallicity, surface gravity) of up to one million stars using the 6dF multi-object spectrograph on the 1.2-m UK Schmidt Telescope of the Anglo-Australian Observatory (AAO). The RAVE program started in 2003, obtaining medium resolution spectra (median R=7,500) in the Ca-triplet region ($\lambda\lambda$ 8,410--8,795 \AA) for southern hemisphere stars drawn from the Tycho-2 and SuperCOSMOS catalogs, in the magnitude range 9<I<12. The first data release is described in this paper and contains radial velocities for 24,748 individual stars (25,274 measurements when including re-observations). Those data were obtained on 67 nights between 11 April 2003 to 03 April 2004. The total sky coverage within this data release is $\sim$4,760 square degrees. The average signal to noise ratio of the observed spectra is 29.5, and 80% of the radial velocities have uncertainties better than 3.4 km/s. Combining internal errors and zero-point errors, the mode is found to be 2 km/s. Repeat observations are used to assess the stability of our radial velocity solution, resulting in a variance of 2.8 km/s. We demonstrate that the radial velocities derived for the first data set do not show any systematic trend with color or signal to noise. The RAVE radial velocities are complemented in the data release with proper motions from Starnet 2.0, Tycho-2 and SuperCOSMOS, in addition to photometric data from the major optical and infrared catalogs (Tycho-2, USNO-B, DENIS and 2MASS). The data release can be accessed via the RAVE webpage: http://www.rave-survey.org.

New images of M31 at 24, 70, and 160 micron taken with the Multiband Imaging Photometer for Spitzer (MIPS) reveal the morphology of the dust in this galaxy. This morphology is well represented by a composite of two logarithmic spiral arms and a circular ring (radius ~10 kpc) of star formation offset from the nucleus. The two spiral arms appear to start at the ends of a bar in the nuclear region and extend beyond the star forming ring. As has been found in previous work, the spiral arms are not continuous but composed of spiral segments. The star forming ring is very circular except for a region near M32 where it splits. The lack of well defined spiral arms and the prominence of the nearly circular ring suggests that M31 has been distorted by interactions with its satellite galaxies. Using new dynamical simulations of M31 interacting with M32 and NGC 205 we find that, qualitatively, such interactions can produce an offset, split ring like that seen in the MIPS images.

Increasingly large populations of disk galaxies are now being observed at increasingly high redshifts, providing new constraints on our knowledge of how such galaxies evolve. Are these observations consistent with a cosmology in which structures form hierarchically? To probe this question, we employ SPH/N-body galaxy scale simulations of late-type galaxies. We examine the evolution of these simulated disk galaxies from redshift 1 to 0, looking at the mass-size and luminosity-size relations, and the thickness parameter, defined as the ratio of scale-height to scale-length. The structural parameters of our simulated disks settle down quickly, and after redshift z=1 the galaxies evolve to become only slightly flatter. Our present day simulated galaxies are larger, more massive, less bright, and redder than at z=1. The inside-out nature of the growth of our simulated galaxies reduces, and perhaps eliminates, expectations of evolution in the size-mass relation.

Seven cosmological hydrodynamic simulations of disk galaxy formation are analyzed to determine the alignment of the disk within the dark matter halo and the internal structure of the halo. We find that the orientation of the outer halo, beyond ~0.1 r_vir, is unaffected by the presence of the disk. In contrast, the inner halo is aligned such that the halo minor axis aligns with the disk axis. The relative orientation of these two regions of the halo are uncorrelated. The alignment of the disk and inner halo appears to take place simultaneously through their joint evolution. The disconnect between these two regions of the halo should be taken into account when modelling tidal streams in the halos of disk galaxies and when calculating intrinsic alignments of disk galaxies based on the properties of dark matter halos.

(Abridged) We investigate how the shapes and angular momenta of galaxy and group mass dark matter halos in a LCDM N-body simulation are correlated internally, and how they are aligned with respect to the location and properties of surrounding halos. We explore these relationships down to halos of much lower mass (10^11/h Msun) than previous studies. The halos are triaxial, with c/a ratios of 0.6+-0.1. More massive halos are more flattened. The principal axes are very well aligned within 0.6 r_vir. The angular momentum vectors are also reasonably well aligned except between the very outermost and very innermost regions of the halo. The angular momentum vectors tend to align with the minor axes, with a mean misalignment of ~25 degrees, and lie perpendicular to the major and intermediate axes. The properties of a halo at 0.4 r_vir are quite characteristic of the properties at most other radii within the halo. There is a very strong tendency for the minor axes of halos to lie perpendicular to large scale filaments. This alignment extends to much larger separations for group and cluster mass halos than for galaxy mass halos. As a consequence, the intrinsic alignments of galaxies are likely weaker than previous predictions, which were based on the shapes of cluster mass halos. The angular momenta of the highest concentration halos tend to point toward other halos. The angular momenta of galaxy mass halos point parallel to filaments, while those of group and cluster mass halos show a very strong tendency to point perpendicular to the filaments. This suggests that group and cluster mass halos acquire most of their angular momentum from major mergers along filaments, while the accretion history of mass and angular momentum onto galaxy mass halos has been smoother.

We investigate whether a correlation exists between the concentration and spin parameters of dark matter haloes forming in LCDM N-body simulations. In particular, we focus on haloes with virial masses characteristic of the hosts of Low Surface Brightness (LSB) and High Surface Brightness (HSB) galaxies. Our analysis demonstrates that concentration and spin are anti-correlated. If we assume that a galaxy disk forms in one of these haloes from baryons that approximately conserved angular momentum during their dissipation into the rotationally supported disk, then it is possible to estimate the disk's central surface density. For an appropriate choice of mass-to-light ratio, we can place constraints on the central surface brightness of the disk and hence identify the analogues of the host haloes of LSB and HSB galaxies. We find that our LSB galaxy analogues occupy haloes that have lower concentrations than might be expected based on consideration of the mass-concentration relation for the LCDM cosmology. The distribution of concentrations peaks at c~6, in good agreement with observational data, although there are important differences between the shapes of the simulated and observationally inferred distributions. This suggests that LSB galaxies inhabit a biased subsample of the halo population, and that this bias may be an important ingredient for resolving the current debate regarding the structure of LSB galaxy dark matter haloes.

We have analyzed galaxy and group-sized dark matter halos formed in a high resolution LCDM numerical N-body simulation in order to study the rotation of the triaxial figure, a property in principle independent of the angular momentum of the particles themselves. Such figure rotation may have observational consequences, such as triggering spiral structure in extended gas disks. The orientation of the major axis is compared at 5 late snapshots of the simulation. Halos with significant substructure or that appear otherwise disturbed are excluded from the sample. We detect smooth figure rotation in 278 of the 317 halos in the sample. The pattern speeds follow a log normal distribution centred at 0.148 h km/s/kpc with a width of 0.83. These speeds are an order of magnitude smaller than required to explain the spiral structure of galaxies such as NGC 2915. The axis about which the figure rotates aligns very well with the halo minor axis, and also reasonably well with its angular momentum vector. The pattern speed is correlated with the halo spin parameter, but shows no correlation with the halo mass. The halos with the highest pattern speeds show particularly strong alignment between their angular momentum vectors and their figure rotation axes. The figure rotation is coherent outside 0.12 r_vir. The measured pattern speed and degree of internal alignment of the figure rotation axis drops in the innermost region of the halo, which may be an artifact of the numerical force softening. The axis ratios show a weak tendency to become more spherical with time.

The perturbations of satellite galaxies, in particular the Large Magellanic Cloud (LMC), have been repeatedly proposed and discounted as the cause of the Milky Way warp. While the LMC may excite a wake in the Galactic dark matter halo that could provide sufficient torque to excite a warp of the observed magnitude, its orbit may be incompatible with the orientation of the warp's line of nodes. The Sgr dSph galaxy has an appropriately-oriented orbit, and due to its closer orbit may produce a stronger tidal effect than the LMC. Evidence that Sgr may be responsible for the warp comes from its orbital angular momentum, which has the same magnitude as the angular momentum of the warped component of the Galactic disk and is anti-aligned with it. We have run high resolution numerical N-body simulations of Sgr-sized satellites around a Milky Way-sized disk to test this idea. Preliminary results suggest that Sgr-sized satellites can indeed excite warps with properties similar to that observed for the Milky Way.

Using recent determinations of the mass and orbit of Sagittarius, I calculate its orbital angular momentum. From the latest observational data, I also calculate the angular momentum of the Milky Way's warp. I find that both angular momenta are directed toward l=270, b=0, and have magnitude 2-8x10^12 M_Sun kpc km s^-1, where the range in both cases reflects uncertainty in the mass. The coincidence of the angular momenta is suggestive of a coupling between these systems. Direct gravitational torque of Sgr on the disk is ruled out as the coupling mechanism. Gravitational torque due to a wake in the halo and the impulsive deposition of momentum by a passage of Sgr through the disk are still both viable mechanisms pending better simulations to test their predictions on the observed Sgr-MW system.

(Abridged) We investigate how galactic disks react to external tidal torques. We calculate the strength and radial dependence of torques on disks that arise from a misalignment between the disk and the main axis system of a flattened dark matter halo. We find that except for in the very inner regions, the torques are well-described by a power law of the form $\tau \propto r^{-2.5}$. The magnitude of the torque is large enough for the entire disk to react to the torque in less than the Hubble time. We demonstrate analytically that disks which are originally located in the $xy$-plane and which are subjected to a torque around the x-axis tilt around the y-axis, as also found in fully non-linear N-body simulations. We further demonstrate that that the torque causes the radius of a chosen particle to increase with time. Investigations of tilting disks which treat the disk as a set of solid rings thus may systematically overestimate the effects of the torque by a factor of two. For torques of the form we investigate, the inner regions of the disk react to the torque faster than the outer regions, resulting in a trailing warp. We then study the effect of the self-gravity of the disk in such a scenario using numerical N-body models. Self-gravity flattens out the inner regions of the disk, but these regions are tilted with respect to their initial plane followed by a non-flat outer region whose tilt decreases with radius. The ``warp radius,'' which marks the end of the inner flat disk, grows throughout the disk at a rate that depends only on the strength of the torque and the local surface density of the disk.

We investigate the tilting and warping of galactic disks in response to tidal torquing. The strength of the torque is determined from cosmological N-body simulations. We find the tidal torques to be dominated by substructure in the galactic halo, such as dwarf satellites, and by a misalignment between the disk angular momentum and the figure axes of the dark matter halo. The radial dependence of the torque can approximated by a power law of index -2.5. A massless disk subjected to such a torque develops a trailing warp similar to those seen in a large number of disk galaxies, i.e. the inner regions of the disk tilt faster than the outer disk. In the case of massive disks, self gravity causes the inner disk to stay locally flat and only the outer regions show the signatures of warps. The radius outside of which a massive disk is warped depends on the local surface density of the disk and on the product of the strength of the torque and the exposure time to this torque.

We show that, in the context of Moore's Law, overall productivity can be increased for large enough computations by `slacking' or waiting for some period of time before purchasing a computer and beginning the calculation.