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.

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

The dwarf satellites of `giant' Milky Way (MW)-mass galaxies are our primary probes of low-mass dark matter halos. The number and velocities of the satellite galaxies of the MW and M31 initially puzzled galaxy formation theorists, but are now reproduced well by many models. Yet, are the MW's and M31's satellites representative? Were galaxy formation models `overfit'? These questions motivate deep searches for satellite galaxies outside the Local Group. We present a deep survey of the `classical' satellites ($M_{\star}$$\geqslant$4$\times$10$^5 M_{\odot}$) of the MW-mass galaxy M94 out to 150 kpc projected distance. We find $only\ two$ satellites, each with $M_{\star}{\sim}10^6 M_{\odot}$, compared with 6-12 such satellites in the four other MW-mass systems with comparable data (MW, M31, M81, M101). Using a `standard' prescription for occupying dark matter halos (taken from the fully hydrodynamical EAGLE simulation) with galaxies, we find that such a sparse satellite population occurs in $<$0.2% of MW-mass systems $-$ a $<$1% probability among a sample of five (known systems + M94). In order to produce an M94-like system more frequently we make satellite galaxy formation much more stochastic than is currently predicted by dramatically increasing the slope and scatter of the stellar mass-halo mass (SMHM) relation. Surprisingly, the SMHM relation must be altered even for halos masses up to 10$^{11}M_{\odot} -$ significantly above the mass scales predicted to have increased scatter from current hydrodynamical models. The sparse satellite population of this `lonely giant' thus advocates for an important modification to ideas of how the satellites around MW-mass galaxies form.

Although the proximity of the Andromeda galaxy (M31) offers an opportunity to understand how mergers affect galaxies, uncertainty remains about M31's most important mergers. Previous studies focused individually on the giant stellar stream or the impact of M32 on M31's disk, thereby suggesting many substantial satellite interactions. Yet models of M31's disk heating and the similarity between the stellar populations of different tidal substructures in M31's outskirts both suggested a single large merger. M31's stellar halo (its outer low-surface-brightness regions) is built up from the tidal debris of satellites and provides information about its important mergers. Here we use cosmological models of galaxy formation to show that M31's massive and metal-rich stellar halo, containing intermediate-age stars, dramatically narrows the range of allowed interactions, requiring a single dominant merger with a large galaxy (with stellar mass about 2.5 x 10^10 solar masses, the third largest member of the Local Group) about 2 Gyr ago. This single event explains many observations that were previously considered separately: M31's compact and metal-rich satellite M32 is likely to be the stripped core of the disrupted galaxy, its rotating inner stellar halo contains most of the merger debris, and the giant stellar stream is likely to have been thrown out during the merger. This interaction may explain M31's global burst of star formation about 2 Gyr ago in which approximately a fifth of its stars were formed. Moreover, M31's disk and bulge were already in place, suggesting that mergers of this magnitude need not dramatically affect galaxy structure.

We investigate the origin of the type III (up-bending) discs based on a sample of 1808 galaxies from SDSS and a sub-sample of 286 galaxies with HI data from ALFALFA. We examine how the type III fraction f3, the fraction of disc galaxies which host up-bending disc breaks, depends on other galactic properties. We confirm that f3 strongly depends on the stellar concentration of galaxies. We find that HI-rich galaxies with low spins tend to have significantly more type III disc breaks than control galaxies, which are matched in concentration and stellar mass. This trend is independent of the existence of strong bars or environment of the galaxies. This result is broadly consistent with predictions from theoretical simulations, and indicates in-situ star formation fueled by gas accretion to be an important mechanism that builds the outer discs of type III galaxies.

There is increasing observational and theoretical evidence for a correlation between the metallicity and the mass of the stellar halo for galaxies with Milky Way-like stellar masses. Using the Illustris cosmological hydrodynamical simulations, we find that this relationship arises because a single massive massive progenitor contributes the bulk of the mass to the accreted stellar component as well as sets its metallicity. Moreover, in the Illustris simulations, this relationship extends over 3 orders of magnitude in accreted stellar mass for central galaxies. We show that for Milky Way-like mass galaxies, the scatter in accreted metallicity at a fixed accreted stellar mass encodes information about the stellar mass of the dominant accreted progenitor, while the radial density and metallicity gradients of the accreted stellar component encodes information about the time of accretion of the dominant progenitor. We demonstrate that for Milky Way-like mass galaxies, the Illustris simulations predict that the metallicity and the stellar mass of the total accreted stellar component can be reconstructed from aperture measurements of the stellar halo along the minor axis of edge-on disk galaxies. These correlations highlight the potential for observational studies of stellar halos to quantify our understanding of the most dominant events in the growth history of galaxies. We explore the implications of our model for our understanding of the accretion histories of the Milky Way, M31 and NGC 5128. In particular, a relatively late and massive accretion is favoured for M31; additionally, we provide a first estimate of the accreted stellar mass for NGC 5128.

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 conducted a study of extra-planar diffuse ionized gas using the first year data from the MaNGA IFU survey. We have stacked spectra from 49 edge-on, late-type galaxies as a function of distance from the midplane of the galaxy. With this technique we can detect the bright emission lines Halpha, Hbeta, [OII]3726, 3729, [OIII]5007, [NII]6549, 6584, and [SII]6717, 6731 out to about 4 kpc above the midplane. With 16 galaxies we can extend this analysis out to about 9 kpc, i.e. a distance of ~2R_e, vertically from the midplane. In the halo, the surface brightnesses of the [OII] and Halpha emission lines are comparable, unlike in the disk where Halpha dominates. When we split the sample by specific star formation rate, concentration index, and stellar mass, each subsample's emission line surface brightness profiles and ratios differ, indicating that extra-planar gas properties can vary. The emission line surface brightnesses of the gas around high specific star formation rate galaxies are higher at all distances, and the line ratios are closer to ratios characteristic of HII regions compared with low specific star formation rate galaxies. The less concentrated and lower stellar mass samples exhibit line ratios that are more like HII regions at larger distances than their more concentrated and higher stellar mass counterparts. The largest difference between different subsamples occurs when the galaxies are split by stellar mass. We additionally infer that gas far from the midplane in more massive galaxies has the highest temperatures and steepest radial temperature gradients based on their [NII]/Halpha and [OII]/Halpha ratios between the disk and the halo.

Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) is an optical fiber-bundle integral-field unit (IFU) spectroscopic survey that is one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV). With a spectral coverage of 3622 - 10,354 Angstroms and an average footprint of ~ 500 arcsec^2 per IFU the scientific data products derived from MaNGA will permit exploration of the internal structure of a statistically large sample of 10,000 low redshift galaxies in unprecedented detail. Comprising 174 individually pluggable science and calibration IFUs with a near-constant data stream, MaNGA is expected to obtain ~ 100 million raw-frame spectra and ~ 10 million reduced galaxy spectra over the six-year lifetime of the survey. In this contribution, we describe the MaNGA Data Reduction Pipeline (DRP) algorithms and centralized metadata framework that produces sky-subtracted, spectrophotometrically calibrated spectra and rectified 3-D data cubes that combine individual dithered observations. For the 1390 galaxy data cubes released in Summer 2016 as part of SDSS-IV Data Release 13 (DR13), we demonstrate that the MaNGA data have nearly Poisson-limited sky subtraction shortward of ~ 8500 Angstroms and reach a typical 10-sigma limiting continuum surface brightness mu = 23.5 AB/arcsec^2 in a five arcsec diameter aperture in the g band. The wavelength calibration of the MaNGA data is accurate to 5 km/s rms, with a median spatial resolution of 2.54 arcsec FWHM (1.8 kpc at the median redshift of 0.037) and a median spectral resolution of sigma = 72 km/s.

The MaNGA Survey (Mapping Nearby Galaxies at Apache Point Observatory) is one of three core programs in the Sloan Digital Sky Survey IV. It is obtaining integral field spectroscopy (IFS) for 10K nearby galaxies at a spectral resolution of R~2000 from 3,622-10,354A. The design of the survey is driven by a set of science requirements on the precision of estimates of the following properties: star formation rate surface density, gas metallicity, stellar population age, metallicity, and abundance ratio, and their gradients; stellar and gas kinematics; and enclosed gravitational mass as a function of radius. We describe how these science requirements set the depth of the observations and dictate sample selection. The majority of targeted galaxies are selected to ensure uniform spatial coverage in units of effective radius (Re) while maximizing spatial resolution. About 2/3 of the sample is covered out to 1.5Re (Primary sample), and 1/3 of the sample is covered to 2.5Re (Secondary sample). We describe the survey execution with details that would be useful in the design of similar future surveys. We also present statistics on the achieved data quality, specifically, the point spread function, sampling uniformity, spectral resolution, sky subtraction, and flux calibration. For our Primary sample, the median r-band signal-to-noise ratio is ~73 per 1.4A pixel for spectra stacked between 1-1.5 Re. Measurements of various galaxy properties from the first year data show that we are meeting or exceeding the defined requirements for the majority of our science goals.

We derive average flux corrections to the \texttt{Model} magnitudes of the Sloan Digital Sky Survey (SDSS) galaxies by stacking together mosaics of similar galaxies in bins of stellar mass and concentration. Extra flux is detected in the outer low surface brightness part of the galaxies, leading to corrections ranging from 0.05 to 0.32 mag for the highest stellar mass galaxies. We apply these corrections to the MPA-JHU (Max-Planck Institute for Astrophysics - John Hopkins University) stellar masses for a complete sample of half a million galaxies from the SDSS survey to derive a corrected galaxy stellar mass function at $z=0.1$ in the stellar mass range $9.5<\log(M_\ast/M_\odot)<12.0$. We find that the flux corrections and the use of the MPA-JHU stellar masses have a significant impact on the massive end of the stellar mass function, making the slope significantly shallower than that estimated by Li \& White (2009), but steeper than derived by Bernardi et al. (2013). This corresponds to a mean comoving stellar mass density of galaxies with stellar masses $\log(M_\ast/M_\odot) \ge 11.0$ that is a factor of 3.36 larger than the estimate by Li \& White (2009), but is 43\% smaller than reported by Bernardi et al. (2013).

A projected mass estimator and a maximum likelihood estimator are developed to best determine the mass of stellar clusters from proper motion data. Both methods effectively account and correct for errors in the velocities, and provide an unbiased and robust estimate of the cluster mass. Using an extensive proper motion study of $\omega$ Centauri by van Leeuwen et al. (2000), we estimate the mass of $\omega$ Centauri using these two different methods to be $2.8 \times 10^{6} M_{\odot} \left[d/5.1 kpc\right]^{3}$. Within this modeling context, the statistical error is $3\%$ and the systematic error is $4\%$.

We study the stellar haloes of galaxies out to 70-100 kpc as a function of stellar mass and galaxy type by stacking aligned $r$ and $g$ band images from a sample of 45508 galaxies from SDSS DR9 in the redshift range $0.06\,\le\,z\,\le\,0.1$ and in the mass range $10^{10.0} M_{\odot} < M_{*} < 10^{11.4} M_{\odot}$r. We derive surface brightness profiles to a depth of almost $\mu_r \sim 32 \,\mathrm{mag\,arcsec}^{-2}$. We find that the ellipticity of the stellar halo is a function of galaxy stellar mass and that the haloes of high concentration ($C > 2.6$) galaxies are more elliptical than those of low concentration ($C < 2.6$) galaxies. The $g$-$r$ colour profile of high concentration galaxies reveals that the $g$-$r$ colour of the stellar population in the stellar halo is bluer than in the main galaxy, and the colour of the stellar halo is redder for higher mass galaxies. We further demonstrate that the full two-dimensional surface intensity distribution of our galaxy stacks can only be fit through multi-component S\'{e}rsic models. Using the fraction of light in the outer component of the models as a proxy for the fraction of accreted stellar light, we show that this fraction is a function of stellar mass and galaxy type. For high concentration galaxies, the fraction of accreted stellar light rises from $30\%$ to $70\%$ for galaxies in the stellar mass range from $10^{10.0} M_{\odot}$ to $10^{11.4} M_{\odot}$. The fraction of accreted light is much smaller in low concentration systems, increasing from $2\%$ to $25\%$ over the same mass range. This work provides important constraints for the theoretical understanding of the formation of stellar haloes of galaxies.

We have combined the semi-analytic galaxy formation model of Guo et al. (2011) with the particle-tagging technique of Cooper et al. (2010) to predict galaxy surface brightness profiles in a representative sample of ~1900 massive dark matter haloes (10^12--10^14 M_sol) from the Millennium II Lambda-CDM N-body simulation. Here we present our method and basic results focusing on the outer regions of galaxies, consisting of stars accreted in mergers. These simulations cover scales from the stellar haloes of Milky Way-like galaxies to the 'cD envelopes' of groups and clusters, and resolve low surface brightness substructure such as tidal streams. We find that the surface density of accreted stellar mass around the central galaxies of dark matter haloes is well described by a Sersic profile, the radial scale and amplitude of which vary systematically with halo mass (M_200). The total stellar mass surface density profile breaks at the radius where accreted stars start to dominate over stars formed in the galaxy itself. This break disappears with increasing M_200 because accreted stars contribute more of the total mass of galaxies, and is less distinct when the same galaxies are averaged in bins of stellar mass, because of scatter in the relation between M_star and M_200. To test our model we have derived average stellar mass surface density profiles for massive galaxies at z~0.08 by stacking SDSS images. Our model agrees well with these stacked profiles and with other data from the literature, and makes predictions that can be more rigorously tested by future surveys that extend the analysis of the outer structure of galaxies to fainter isophotes. We conclude that it is likely that the outer structure of the spheroidal components of galaxies is largely determined by collisionless merging during their hierarchical assembly

We lay out and apply methods to use proper motions of individual kinematic tracers for estimating the dynamical mass of star clusters. We first describe a simple projected mass estimator and then develop an approach that evaluates directly the likelihood of the discrete kinematic data given the model predictions. Those predictions may come from any dynamical modelling approach, and we implement an analytic King model, a spherical isotropic Jeans equation model and an axisymmetric, anisotropic Jeans equation model.We apply these approaches to the enigmatic globular cluster omega Centauri, combining the proper motion from van Leeuwen et al (2000) with improved photometric cluster membership probabilities. We show that all mass estimates based on spherical isotropic models yield $(4.55\pm 0.1) \times 10^6 M_{\odot} [D/5.5 \pm 0.2 kpc]^3$, where our modelling allows us to show how the statistical precision of this estimate improves as more proper motion data of lower signal-to-noise are included. MLM predictions, based on an anisotropic axisymmetric Jeans model, indicate for $\omega$ Cen that the inclusion of anisotropies is not important for the mass estimates, but that accounting for the flattening is: flattened models imply $(4.05\pm 0.1) \times 10^6 M_{\odot} [D/5.5 \pm 0.2 kpc]^3$, 10% lower than when restricting the analysis to a spherical model. The best current distance estimates imply an additional uncertainty in the mass estimate of 12%.