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 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.

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.

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.

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 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.

We present detailed analysis of color-magnitude diagrams of NGC2403, obtained from a deep (m<28) Hubble Space Telescope (HST) Wide Field Planetary Camera 2 observation of the outer disk of NGC2403, supplemented by several shallow (m<26) HST Advanced Camera for Surveys fields. We derive the spatially resolved star formation history of NGC2403 out to 11 disk scale lengths. In the inner portions of the galaxy, we compare the recent star formation rates (SFRs) we derive from the resolved stars with those measured using GALEX FUV + Spitzer 24-micron fluxes, finding excellent agreement between the methods. Our measurements also show that the radial gradient in recent SFR mirror s the disk exponential profile to 11 scale lengths with no break, extending to SFR densities a factor of 100 lower than those that can be measured with GALEX and Spitzer (2x10^{-6} M_{\sun} yr^{-1} kpc^{-2}). Furthermore, we find that the cumulative stellar mass of the disk was formed at similar times at all radii. We compare these characteristics of NGC2403 to those of its "morphological twins," NGC300 and M33, showing that the structure and age distributions of the NGC2403 disk are more similar to those of the relatively isolated system NGC300 than to those of the Local Group analog M33. We also discuss the environments and HI morphologies of these three nearby galaxies, comparing them to integrated light studies of larger samples of more distant galaxy disks. Taken together, the physical properties and evolutionary history of NGC2403 suggest that the galaxy has had no close encounters with other M81 group members and may be falling into the group for the first time.

We analyzed the radial surface brightness profile of the spiral galaxy NGC 7793 using HST/ACS images from the GHOSTS survey and a new HST/WFC3 image across the disk break. We used the photometry of resolved stars to select distinct populations covering a wide range of stellar ages. We found breaks in the radial profiles of all stellar populations at 280" (~5.1 kpc). Beyond this disk break, the profiles become steeper for younger populations. This same trend is seen in numerical simulations where the outer disk is formed almost entirely by radial migration. We also found that the older stars of NGC 7793 extend significantly farther than the underlying HI disk. They are thus unlikely to have formed entirely at their current radii, unless the gas disk was substantially larger in the past. These observations thus provide evidence for substantial stellar radial migration in late-type disks.

We present Herschel/SPIRE images at 250, 350, and 500 {\mu}m of NGC 4244, a typical low-mass, disk-only and edge-on spiral galaxy. The dust disk is clumpy and shows signs of truncation at the break radius of the stellar disk. This disk coincides with the densest part of the Hi disk. We compare the Spectral Energy Distribution, including the new SPIRE fluxes, to 3D radiative transfer models; a smooth model disk and a clumpy model with embedded heating. Each model requires a very high value for the dust scale-length (h(dust) = 2 - 5 h(stars)), higher dust masses than previous models of NGC 4244 (Md = 0.47 - 1.39 \times 10e7 Msun) and a face-on optical depth of {\tau}(V) = 0.4 - 1.12, in agreement with previous disk opacity studies. The vertical scales of stars and dust are similar. The clumpy model much better mimics the general morphology in the submm images and the general SED. The inferred gas-to-dust mass ratio is compatible with those of similar low-mass disks. The relatively large radial scale-length of the dust disk points to radial mixing of the dusty ISM within the stellar disk. The large vertical dust scale and the clumpy dust distribution of our SED model are both consistent with a scenario in which the vertical structure of the ISM is dictated by the balance of turbulence and self-gravity.

The presence of a well-defined and narrow dust lane in an edge-on spiral galaxy is the observational signature of a thin and dense molecular disk, in which gravitational collapse has overcome turbulence. Using a sample of galaxies out to z~1 extracted from the COSMOS survey, we identify the fraction of massive disks that display a dust lane. Our goal is to explore the evolution in the stability of the molecular ISM disks in spiral galaxies over a cosmic timescale. We check the reliability of our morphological classifications against changes in restframe wavelength, resolution, and cosmic dimming with (artificially redshifted) images of local galaxies from SDSS. We find that the fraction of L* disks with dust lanes in COSMOS is consistent with the local fraction (~80%) out to z~0.7. At z=0.8, the dust lane fraction is only slightly lower. A somewhat lower dust lane fraction in starbursting galaxies tentatively supports the notion that a high specific star formation rate can efficiently destroy or inhibit a dense molecular disk. A small subsample of higher redshift COSMOS galaxies display low internal reddening (E[B-V]), as well as a low incidence of dust lanes. These may be disks in which the growth of the dusty ISM disk lags behind that of the stellar disk. We note that at z=0.8, the most massive galaxies display a lower dust lane fraction than lower mass galaxies. A small contribution of recent mergers or starbursts to this most massive population may be responsible. The fact that the fraction of galaxies with dust lanes in COSMOS is consistent with little or no evolution implies that models to explain the Spectral Energy Distribution or the host galaxy dust extinction of supernovae based on local galaxies are still applicable to higher redshift spirals. It also suggests that dust lanes are long lived phenomena or can be reformed over very short time-scales.

Edge-on spiral galaxies offer a unique perspective on the vertical structure of spiral disks, both stars and the iconic dark dustlanes. The thickness of these dustlanes can now be resolved for the first time with Herschel in far-infrared and sub-mm emission. We present NHEMESES, an ongoing project that targets 12 edge-on spiral galaxies with the PACS and SPIRE instruments on Herschel. These vertically resolved observations of edge-on spirals will impact on several current topics. First and foremost, these Herschel observations will settle whether or not there is a phase change in the vertical structure of the ISM with disk mass. Previously, a dramatic change in dustlane morphology was observed as in massive disks the dust collapses into a thin lane. If this is the case, the vertical balance between turbulence and gravity dictates the ISM structure and consequently star-formation and related phenomena (spiral arms, bars etc.). We specifically target lower mass nearby edge-ons to complement existing Herschel observations of high-mass edge-on spirals (the HEROES project). Secondly, the combined data-set, together with existing Spitzer observations, will drive a new generation of spiral disk Spectral Energy Distribution models. These model how dust reprocesses starlight to thermal emission but the dust geometry remains the critical unknown. And thirdly, the observations will provide an accurate and unbiased census of the cold dusty structures occasionally seen extending out of the plane of the disk, when backlit by the stellar disk. To illustrate the NHEMESES project, we present early results on NGC 4244 and NGC 891, two well studies examples of a low and high-mass edge-on spiral.

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.

Our GHOSTS survey measures the stellar envelope properties of 14 nearby disk galaxies by imaging their resolved stellar populations with HST/ACS&WFPC2. Most of the massive galaxies in the sample (Vrot >200 km/s) have very extended stellar envelopes with Sersic law profiles or mu(r)~r^-2.5 power law profiles in the outer regions. For these massive galaxies we can fit the central bulge light and the outer halo out to 30 kpc with one and the same Sersic profile and the stellar surface density of the profiles correlate with Hubble type and bulge- to-disk ratio. This suggests that the central bulges and inner halos are created in the same process. Smaller galaxies (Vrot ~100 km/s) have much smaller stellar envelopes, but depending on geometry, they could still be more luminous than expected from satellite remnants in hierarchical galaxy formation models. Alternatively, they could be created by disk heating through the bombardment of small dark matter sub-halos. The halos we fit are highly flattened, with minor- over-ma jor axis ratios c/a=0.6. The halos are somewhat more compact than hierarchical model predictions. The halos show small metallicity gradient out to 30 kpc and the massive galaxies have typical [Fe/H]~-0.8. We find indications of halo substructure in many galaxies, but some halos seem remarkable smooth.

In hierarchical galaxy formation the stellar halos of galaxies are formed by the accretion of minor satellites and therefore contain valuable information about the (early) assembly process of galaxies. Our GHOSTS survey measures the stellar envelope properties of 14 nearby disk galaxies by imaging their resolved stellar populations with HST/ACS&WFPC2. Most of the massive galaxies in the sample (Vrot>200 km/s) have very extended stellar envelopes with mu(r) ~ r^{-2.5} power law profiles in the outer regions. For these massive galaxies there is some evidence that the stellar surface density of the profiles correlates with Hubble type and bulge-to-disk ratio, begging the question whether these envelopes are more related to bulges than to a Milky Way-type stellar halo. Smaller galaxies (Vrot 100 km/s) have much smaller stellar envelopes, but depending on geometry, they could still be more luminous than expected from satellite remnants in hierarchical galaxy formation models. Alternatively, they could be created by disk heating through the bombardment of small dark matter sub-halos. We find that galaxies show varying amounts of halo substructure.

In a hierarchical merging scenario, the outer parts of a galaxy are a fossil record of the galaxy's early history. Observations of the outer disks and halos of galaxies thus provide a tool to study individual galaxy histories and test formation theories. Locally, an impressive effort has been made to understand the halo of the Milky Way, Andromeda, and M33. However, due to the stochastic nature of halo formation, a better understanding of this process requires a large sample of galaxies with known halo properties. The GHOSTS project (PI: R. de Jong) aims to characterize the halos and outer portions of 14 nearby (D=4-14 Mpc) spiral galaxies using the Hubble Space Telescope. Detection of individual stars in the outer parts of these galaxies enables us to study both the morphological properties of the galaxies, and determine the stars' metallicity and age.

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 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.

To further our understanding of the Great Attractor (GA), we have undertaken a redshift survey using the 2dF on the AAT. Clusters and filaments in the GA region were targeted with 25 separate pointings resulting in approximately 2600 new redshifts. Targets included poorly studied X-ray clusters from the CIZA catalogue as well as the Cen-Crux and PKS 1343-601 clusters, both of which lie close to the classic GA centre. For nine clusters in the region, we report velocity distributions as well as virial and projected mass estimates. The virial mass of CIZA J1324.7-5736, now identified as a separate structure from the Cen-Crux cluster, is found to be ~3x10^14 M_sun, in good agreement with the X-ray inferred mass. In the PKS 1343-601 field, five redshifts are measured of which four are new. An analysis of redshifts from this survey, in combination with those from the literature, reveals the dominant structure in the GA region to be a large filament, which appears to extend from Abell S0639 (l=281\deg, b=+11\deg) to (l~5\deg, b~-50\deg), encompassing the Cen-Crux, CIZA J1324.7-5736, Norma and Pavo II clusters. Behind the Norma Cluster at cz~15000 km/s, the masses of four rich clusters are calculated. These clusters (Triangulum-Australis, Ara, CIZA J1514.6-4558 and CIZA J1410.4-4246) may contribute to a continued large-scale flow beyond the GA. The results of these observations will be incorporated into a subsequent analysis of the GA flow.

We compare the measured peculiar velocities of 98 local (<150/h Mpc) type Ia supernovae (SNIa) with predictions derived from the PSCz. There is excellent agreement between the two datasets with a best fit beta_I (=Omega_m^0.6/b_I) of 0.55+-0.06. Subsets of the SNIa dataset are further analysed and the above result is found to be robust with respect to culls by distance, host-galaxy extinction and to the reference frame in which the analysis is carried out. We briefly review the peculiar motions in the direction of the Great Attractor. Most clusters in this part of the sky out to a distance of 14,000 km/s, i.e. those closer than the Shapley Concentration, have sizable positive peculiar velocities, i.e. (~ +400 km/s). There are nine local SNIa in the GA direction that are in the foreground of Shapley. All these SNIa have positive peculiar velocities. Hence both the cluster and local SNIa data strongly support the idea of a sizable flow into Shapley.

We compare the measured peculiar velocities of 98 local (<150/h Mpc) type Ia supernovae with predictions derived from the PSCz survey. There is excellent agreement between the two datasets with a best fit beta_I (=Omega_m^0.6/b_I) of 0.55+-0.06. Subsets of the supernovae dataset are further analysed and the above result is found to be robust with respect to culls by distance, host-galaxy extinction and to the choice of reference frame in which the analysis is carried out. Alternative methods of determining beta_I including density-density comparisons, dipole measurements and WMAP-based results are also discussed. We conclude that most recent determinations are consistent with a value of beta_I=0.5.