We present an extensive study of spectroscopically-derived chemical abundances for M87 and its globular cluster (GC) system. Using observations from the Mitchell spectrograph at McDonald, LRIS at Keck, and Hectospec on the MMT, we derive new metallicity gradients from $\sim 2$ to $140$ kpc. We use a novel hierarchical statistical framework to simultaneously separate the GC system into subpopulations while measuring the metallicity gradients of those subpopulations. We create physically-motivated spectral stacks of the GC subpopulations by leveraging the output of this statistical framework to perform the first application of abundance tagging in a massive ETG to better constrain the origins of the GC subpopulations and, thus, the assembly history of M87. We find a metal-poor, $\alpha$-enhanced population of GCs in both in the inner and outer halo unanticipated by current cosmological simulations of galaxy evolution. We use the remarkably flat metallicity gradients we find for both the metal-rich and metal-poor GC subpopulations in the inner halo as tentative evidence that some amount of the metal-poor GCs formed directly in the halo of M87 at high redshift.

We present spatially-resolved stellar kinematics of the well-studied ultra diffuse galaxy (UDG) Dragonfly 44, as determined from 25.3 hrs of observations with the Keck Cosmic Web Imager. The luminosity-weighted dispersion within the half-light radius is $\sigma_{1/2}=33^{+3}_{-3}$ km/s. There is no evidence for rotation, with $V/\sigma<0.12$ (90% confidence) along the major axis, in apparent conflict with models where UDGs are the high-spin tail of the normal dwarf galaxy distribution. The spatially-averaged line profile is more peaked than a Gaussian, with Gauss-Hermite coefficient $h_4=0.13\pm 0.05$. The mass-to-light ratio within the effective radius is $M/L=26^{+7}_{-6}$, similar to other UDGs and higher by a factor of six than normal galaxies of the same luminosity. This difference between UDGs and other galaxies is, however, sensitive to the aperture that is used, and is much reduced when the $M/L$ ratios are measured within a fixed radius of 10 kpc. Dragonfly 44 has a rising velocity dispersion profile, from $\sigma=26^{+4}_{-4}$ km/s at R=0.2 kpc to $\sigma=41^{+8}_{-8}$ km/s at R=5.1 kpc. The profile can only be fit with a cuspy NFW profile if the orbital distribution has strong tangential anisotropy, with $\beta=-0.8^{+0.4}_{-0.5}$. An alternative explanation is that the dark matter profile has a core: a Di Cintio et al. (2014) density profile with a mass-dependent core provides a very good fit to the kinematics for a halo mass of $\log (M_{200}/{\rm M}_{\odot})=11.2^{+0.6}_{-0.6}$ and $\beta=-0.1^{+0.2}_{-0.3}$, i.e., isotropic orbits. This model predicts a slight positive kurtosis, in qualitative agreement with the measured $h_4$ parameter. UDGs such as Dragonfly 44 are dark matter dominated even in their centers, and can constrain the properties of dark matter in a regime where baryons usually dominate the kinematics: small spatial scales in massive halos.

The observed characteristics of globular cluster (GC) systems, such as metallicity distributions, are commonly used to place constraints on galaxy formation models. However, obtaining reliable metallicity values is particularly difficult because of our limited means to obtain high quality spectroscopy of extragalactic GCs. Often, `color--metallicity relations' are invoked to convert easier-to-obtain photometric measurements into metallicities, but there is no consensus on what form these relations should take. In this paper we make use of multiple photometric datasets and iron metallicity values derived from applying full-spectrum stellar population synthesis models to deep Keck/LRIS spectra of 177 GCs centrally located around M87 to obtain a new color--metallicity relation. Our new relation differs substantially from previous relations in the blue, and we present evidence that the M87 relation differs from that of the Milky Way GCs, suggesting environmental dependence of GC properties. We use our color--metallicity relation to derive a new GC metallicity-host galaxy luminosity relation for red and blue GCs and find a shallower relation for the blue GCs than what previous work has found and that the metal-poor GCs are more enriched than what was previously found. This could indicate that the progenitor satellite galaxies that now make up the stellar halos of early-type galaxies are more massive and formed later than previously thought, or that the properties of metal-poor GCs are less dependent on their present-day host, indicating a common origin.

Given the absence of directly detected dark matter (DM) as weakly interacting massive particles, there is strong interest in the possibility that DM is an ultra-light scalar field, here denoted as fuzzy DM. Ultra-diffuse galaxies, with the sizes of giant galaxies and the luminosities of dwarf galaxies, have a wide range of DM halo masses, thus providing new opportunities for exploring the connections between galaxies and their DM halos. Following up on new integral field unit spectroscopic observations and dynamics modeling of the DM-dominated ultra-diffuse galaxy Dragonfly 44 in the outskirts of the Coma Cluster, we present models of fuzzy DM constrained by the stellar dynamics of this galaxy. We infer a scalar field mass of $\sim 3 \times 10^{-22}$ eV, consistent with other constraints from galaxy dynamics but in tension with constraints from Ly$\alpha$ forest power spectrum modeling. While we are unable to statistically distinguish between fuzzy DM and normal cold DM models, we find that the inferred properties of the fuzzy DM halo satisfy a number of predictions for halos in a fuzzy DM cosmology. In particular, we find good agreement with the predicted core size--halo mass relation and the predicted transition radius between the quantum pressure-dominated inner region and the outer halo region.

Ultra diffuse galaxies are a population of extended galaxies but with relatively low luminosities. The origin of these objects remains unclear, largely due to the observational challenges of the low surface brightness Universe. We present here a detailed stellar population analysis of a relatively isolated UDG, DGSAT I, based on spectroscopic data from the Keck Cosmic Web Imager integral field unit. The star formation history of DGSAT I seems to be extended, with a mean luminosity-weighted age of ~3 Gyr, in agreement with previous photometric studies. However, we find a very high [Mg/Fe] abundance ratio, which is extreme even in the context of the highly alpha-enhanced massive ellipticals and ultra-faint dwarfs. The [Mg/Fe]-enhancement of DGSAT I appears to be 10 times higher than the most magnesium-enhanced stellar systems discovered to date, and suggests that the chemical enrichment of this object was dominated by core-collapse supernovae. Intriguingly, this breaks the canonical relation between [Mg/Fe] and star formation time-scale. With a measured velocity dispersion of 56 +/- 10 km/s, DGSAT I also shows a high mass-to-light ratio, which indicates that it is highly dark matter-dominated. The metal-poor conditions of DGSAT I may have enhanced the formation of massive stars, while at the same time, additional mechanisms are needed to prevent iron-rich yields from being recycled into stars. These results suggest that some ultra-diffuse galaxies could have experienced chemical enrichment episodes similar to the first building blocks of galaxies.

The stellar content of the intracluster light (ICL) provides unique insight into the hierarchical assembly process of galaxy clusters.However, the ICL is difficult to study due to its low surface brightness and large physical extent. We present optical spectra of three ICL regions in the Coma cluster, located between 100-180kpc from their nearest BCGs: NGC4889 and NGC4874. The mean surface brightness of the three ICL regions are {\mu}$_g$~25.3-26.2mag arcsec$^{-2}$. IFU spectroscopy with 13.5 hr on-source integration time were acquired as part of an ancillary program within the SDSS-IV MaNGA survey. We stacked the 127 individual fiber spectra in each IFU in order to achieve a 1{\sigma} limiting surface brightness of 27.9mag arcsec$^{-2}$, corresponding to a mean S/N in the optical of 21.6,9.6,and 11.6\AA$^{-1}$. We apply stellar population models to the stacked spectra, and measure the recession velocities, velocity dispersions ($\sigma$), stellar ages, and [Fe/H]. Our results show that the $\sigma$ of ICL regions are very high, indicating the stars are tracing the gravitational potential of the cluster, instead of any individual galaxy. The line-of-sight velocities of the three ICL regions are different from each other by ~700km/s, while the velocity of each region is similar to the closest BCG. This suggests that the ICL regions are associated with two distinct subclusters centered on NGC4889 and NGC4874.The stellar populations of these regions are old and metal poor, with ages of 7-12Gyr, and [Fe/H] of -0.8 to -0.6 dex. From the derived age and [Fe/H], the build-up of ICL in Coma is likely to be through the accretion of low mass galaxies or the tidal stripping of the outskirts of massive galaxies that have ended their star formation early on, instead of directly from major mergers of massive galaxies.

Using data from the WiFeS Atlas of Galactic Globular cluster Spectra we study the behaviour of the calcium triplet (CaT), a popular metallicity indicator in extragalactic stellar population studies. A major caveat of these studies is that the potential sensitivity to other stellar population parameters such as age, calcium abundance and the initial mass function has not yet been empirically evaluated. Here we present measurements of the strength of the CaT feature for 113 globular clusters in the Milky Way and its satellite galaxies. We derive empirical calibrations between the CaT index and both the iron abundance ([Fe/H]) and calcium abundance ([Ca/H]), finding a tighter relationship for [Ca/H] than for [Fe/H]. For stellar populations 3 Gyr and older the CaT can be used to reliably measure [Ca/H] at the 0.1 dex level but becomes less reliable for ages of $\sim 2$ Gyr and younger. We find that the CaT is relatively insensitive to the horizontal branch morphology. The stellar mass function however affects the CaT strengths significantly only at low metallicities. Using our newly derived empirical calibration, we convert our measured CaT indices into [Ca/H] values for the globular clusters in our sample.

We present observational constraints on the stellar populations of two ultra-diffuse galaxies (UDGs) using optical through near-infrared (NIR) spectral energy distribution (SED) fitting. Our analysis is enabled by new $Spitzer$-IRAC 3.6 $\mu$m and 4.5 $\mu$m imaging, archival optical imaging, and the prospector fully Bayesian SED fitting framework. Our sample contains one field UDG (DGSAT I), one Virgo cluster UDG (VCC 1287), and one Virgo cluster dwarf elliptical for comparison (VCC 1122). We find that the optical--NIR colors of the three galaxies are significantly different from each other. We infer that VCC 1287 has an old ($\gtrsim7.7$ Gyr) and surprisingly metal-poor ($[Z/Z_{\odot}]\lesssim-1.0$) stellar population, even after marginalizing over uncertainties on diffuse interstellar dust. In contrast, the field UDG DGSAT I shows evidence of being younger than the Virgo UDG, with an extended star formation history and an age posterior extending down to $\sim3$ Gyr. The stellar metallicity of DGSAT I is sub-solar but higher than that of the Virgo UDG, with $[Z/Z_{\odot}]=-0.63^{+0.35}_{-0.62}$; in the case of exactly zero diffuse interstellar dust, DGSAT I may even have solar metallicity. With VCC 1287 and several Coma UDGs, a general picture is emerging where cluster UDGs may be "failed" galaxies, but the field UDG DGSAT I seems more consistent with a stellar feedback-induced expansion scenario. In the future, our approach can be applied to a large and diverse sample of UDGs down to faint surface brightness limits, with the goal of constraining their stellar ages, stellar metallicities, and circumstellar and diffuse interstellar dust content.

We present updated stellar population models appropriate for old ages (>1 Gyr) and covering a wide range in metallicities (-1.5<[Fe/H]<0.3). These models predict the full spectral variation associated with individual element abundance variation as a function of metallicity and age. The models span the optical-NIR wavelength range (0.37-2.4um), include a range of initial mass functions (IMFs) and contain the flexibility to vary 18 individual elements including C, N, O, Mg, Si, Ca, Ti, and Fe. To test the fidelity of the models we fit them to integrated light optical spectra of 41 Galactic globular clusters (GCs). The value of testing models against GCs is that their ages, metallicities, and detailed abundance patterns have been derived from the HR diagram in combination with high resolution spectroscopy of individual stars. We determine stellar population parameters from fits to all wavelengths simultaneously ("full spectrum fitting"), and demonstrate explicitly with mock tests that this approach produces smaller uncertainties at fixed S/N ratio than fitting a standard set of 14 line indices. Comparison of our integrated-light results to literature values reveals good agreement in metallicity, [Fe/H]. When restricting to GCs without prominent blue horizontal branch populations we also find good agreement with literature values for ages, [Mg/Fe], [Si/Fe], and [Ti/Fe].

We investigate the dark matter density profile of the massive elliptical galaxy, NGC 1407, by constructing spherically symmetric Jeans models of its field star and globular cluster systems. Two major challenges in such models are the degeneracy between the stellar mass and the dark matter halo profiles, and the degeneracy between the orbital anisotropy of the tracer population and the total mass causing the observed motions. We address the first issue by using new measurements of the mass-to-light ratio profile from stellar population constraints that include a radially varying initial mass function. To deal with the mass--anisotropy degeneracy, we make use of multiple kinematic tracers, including two subpopulations of globular clusters in addition to the galaxy's field stars. We create a hierarchical Bayesian model that addresses several often neglected systematic uncertainties, such as the statistical weight given to various datasets and the adopted distance. After sampling the posterior probability distribution with a Markov Chain Monte Carlo method, we find evidence for a central cusp with a log-slope of $\gamma = 1.0^{+0.3}_{-0.4}$. This is lower than expected for dark matter halos that have undergone adiabatic contraction, supporting inferences from gravitational lensing that some process has suppressed the steepening of halos in massive galaxies. We also confirm radially-biased orbits for the metal-rich globular clusters and tangentially-biased orbits for the metal-poor globular clusters, which remains a puzzling finding for an accretion-dominated halo.

Analyses of strong gravitational lenses, galaxy-scale kinematics, and absorption line stellar population synthesis (SPS) have all concluded that the stellar initial mass function (IMF) varies within the massive early-type galaxy (ETG) population. However, the physical mechanism that drives variation in the IMF is an outstanding question. Here we use new SPS models to consider a diverse set of compact, low-velocity dispersion stellar systems: globular clusters (GCs), an ultra-compact dwarf (UCD), and the compact elliptical (cE) galaxy M32. We compare our results to massive ETGs and available dynamical measurements. We find that the GCs have stellar mass-to-light ratios (M/L) that are either consistent with a Kroupa IMF or are slightly bottom-light while the UCD and cE have mildly elevated M/L. The separation in derived IMFs for systems with similar metallicities and abundance patterns indicates that our SPS models can distinguish abundance and IMF effects. Variation among the sample in this paper is only $\sim 50\%$ in normalized M/L compared to the $\sim 4\times$ among the ETG sample. This suggests that metallicity is not the sole driver of IMF variability and additional parameters need to be considered.

We present an analysis of stellar populations in passive galaxies in seven massive X-ray clusters at z=0.19-0.89. Based on absorption line strengths measured from our high signal-to-noise spectra, the data support primarily passive evolution of the galaxies. We use the scaling relations between velocity dispersions and the absorption line strengths to determine representative mean line strengths for the clusters. From the age determinations based on the line strengths (and stellar population models), we find a formation redshift of z_form=1.96(-0.19,+0.24). Based on line strength measurements from high signal-to-noise composite spectra of our data, we establish the relations between velocity dispersion, ages, metallicities [M/H] and abundance ratios [alpha/Fe] as a function of redshift. The [M/H]-velocity dispersion and [alpha/Fe]-velocity dispersion relations are steep and tight. The age-velocity dispersion relation is flat, with zero point changes reflecting passive evolution. The scatter in all three parameters are within 0.08-0.15 dex at fixed velocity dispersions, indicating a large degree of synchronization in the evolution of the galaxies. We find indication of cluster-to-cluster differences in metallicities and abundance ratios. However, variations in stellar populations with the cluster environment can only account for a very small fraction of the intrinsic scatter in the scaling relations. Thus, within these very massive clusters the main driver of the properties of the stellar populations in passive galaxies appears to be the galaxy velocity dispersion.

We present new observations of the three nearest early-type galaxy (ETG) strong lenses discovered in the SINFONI Nearby Elliptical Lens Locator Survey (SNELLS). Based on their lensing masses, these ETGs were inferred to have a stellar initial mass function (IMF) consistent with that of the Milky Way, not the bottom-heavy IMF that has been reported as typical for high-$\sigma$ ETGs based on lensing, dynamical, and stellar population synthesis techniques. We use these unique systems to test the consistency of IMF estimates derived from different methods. We first estimate the stellar $M_*/L$ using lensing and stellar dynamics. We then fit high-quality optical spectra of the lenses using an updated version of the stellar population synthesis models developed by Conroy & van Dokkum. When examined individually, we find good agreement among these methods for one galaxy. The other two galaxies show 2-3$\sigma$ tension with lensing estimates, depending on the dark matter contribution, when considering IMFs that extend to 0.08 Msol. Allowing a variable low-mass cutoff or a nonparametric form of the IMF reduces the tension among the IMF estimates to $<2\sigma$. There is moderate evidence for a reduced number of low-mass stars in the SNELLS spectra, but no such evidence in a composite spectrum of matched-$\sigma$ ETGs drawn from the SDSS. Such variation in the form of the IMF at low stellar masses (m <~ 0.03 Msol), if present, could reconcile lensing/dynamical and spectroscopic IMF estimates for the SNELLS lenses and account for their lighter $M_*/L$ relative to the mean matched-$\sigma$ ETG. We provide the spectra used in this study to facilitate future comparisons.

We present a $0.7-2.5\mu m$ spectral library of 284 stars observed with the medium-resolution infrared spectrograph, SpeX, at the 3.0 meter NASA Infrared Telescope Facility (IRTF) on Maunakea, Hawaii. This library extends the metallicity range of the IRTF Cool Star library beyond solar metallicity to $-1.7 <$ [Fe/H] $< 0.6$. All of the observed stars are also in the MILES optical stellar library, providing continuous spectral coverage for each star from $0.35-2.5\mu m$. The spectra are absolute flux calibrated using Two Micron All Sky Survey photometry and the continuum shape of the spectra is preserved during the data reduction process. Synthesized $JHK_S$ colors agree with observed colors at the $1-2\%$ level, on average. We also present a spectral interpolator that uses the library to create a data-driven model of spectra as a function of $teff$, $logg$, and [Fe/H]. We use the library and interpolator to compare empirical trends with theoretical predictions of spectral feature behavior as a function of stellar parameters. These comparisons extend to the previously difficult to access low-metallicity and cool dwarf regimes, as well as the previously poorly sampled super-solar metallicity regime. The library and interpolator are publicly available.

There is good evidence that the centers of massive early-type galaxies have a bottom-heavy stellar initial mass function (IMF) compared to the IMF of the Milky Way. Here we study the radial variation of the IMF within such galaxies, using a combination of high quality Keck spectroscopy and a new suite of stellar population synthesis models that cover a wide range in metallicity. As in the previous studies in this series, the models are fitted directly to the spectra and treat all elemental abundance ratios as free parameters. Using newly obtained spectroscopy for six galaxies, including deep data extending to ~1Re for the galaxies NGC1407, NGC1600, and NGC2695, we find that the IMF varies strongly with galactocentric radius. For all six galaxies the IMF is bottom-heavy in the central regions, with average mass-to-light ratio "mismatch" parameter a~2.5 at R=0. The IMF rapidly becomes more bottom-light with increasing radius, flattening off near the Milky Way value (a~1.1) at R>0.4Re. A consequence is that the luminosity-weighted average IMF depends on the measurement aperture: within R=Re we find <a>=1.3-1.5, consistent with recent lensing and dynamical results from SLACS and ATLAS-3D. Our results are also consistent with several earlier studies that were based on analyses of radial gradients of line indices. The observed IMF gradients support galaxy formation models in which the central regions of massive galaxies had a different formation history than their outer parts. Finally, we make use of the high signal-to-noise central spectra of NGC1407 and NGC2695 to demonstrate how we can disentangle IMF effects and abundance effects.

We present the WiFeS Atlas of Galactic Globular cluster Spectra, a library of integrated spectra of Milky Way and Local Group globular clusters. We used the WiFeS integral field spectrograph on the Australian National University 2.3 m telescope to observe the central regions of 64 Milky Way globular clusters and 22 globular clusters hosted by the Milky Way's low mass satellite galaxies. The spectra have wider wavelength coverage (3300 {\AA} to 9050 {\AA}) and higher spectral resolution (R = 6800) than existing spectral libraries of Milky Way globular clusters. By including Large and Small Magellanic Cloud star clusters, we extend the coverage of parameter space of existing libraries towards young and intermediate ages. While testing stellar population synthesis models and analysis techniques is the main aim of this library, the observations may also further our understanding of the stellar populations of Local Group globular clusters and make possible the direct comparison of extragalactic globular cluster integrated light observations with well understood globular clusters in the Milky Way. The integrated spectra are publicly available via the project website.

It is now well-established that the stellar initial mass function (IMF) can be determined from the absorption line spectra of old stellar systems, and this has been used to measure the IMF and its variation across the early-type galaxy population. Previous work focused on measuring the slope of the IMF over one or more stellar mass intervals, implicitly assuming that this is a good description of the IMF and that the IMF has a universal low-mass cutoff. In this work we consider more flexible IMFs, including two-component power-laws with a variable low-mass cutoff and a general non-parametric model. We demonstrate with mock spectra that the detailed shape of the IMF can be accurately recovered as long as the data quality are high (S/N$\gtrsim300$) and cover a wide wavelength range (0.4um-1.0um). We apply these flexible IMF models to a high S/N spectrum of the center of the massive elliptical galaxy NGC 1407. Fitting the spectrum with non-parametric IMFs, we find that the IMF in the center shows a continuous rise extending toward the hydrogen-burning limit, with a behavior that is well-approximated by a power-law with an index of -2.7. These results provide strong evidence for the existence of extreme (super-Salpeter) IMFs in the cores of massive galaxies.

We utilise the DEIMOS instrument on the Keck telescope to measure the wide-field stellar kinematics of early-type galaxies as part of the SAGES Legacy Unifying Globulars and GalaxieS (SLUGGS) survey. In this paper, we focus on some of the lowest stellar mass lenticular galaxies within this survey, namely NGC 2549, NGC 4474, NGC 4459 and NGC 7457, performing detailed kinematic analyses out to large radial distances of $\sim 2-3$ effective radii. For NGC 2549, we present the first analysis of data taken with the SuperSKiMS (Stellar Kinematics from Multiple Slits) technique. To better probe kinematic variations in the outskirts of the SLUGGS galaxies, we have defined a local measure of stellar spin. We use this parameter and identify a clear separation in the radial behaviour of stellar spin between lenticular and elliptical galaxies, thereby reinforcing the physically meaningful nature of their morphological classifications. We compare the kinematic properties of our galaxies with those from various simulated galaxies to extract plausible formation scenarios. By doing this for multiple simulations, we assess the consistency of the theoretical results. Comparisons to binary merger simulations show that low-mass lenticular galaxies generally resemble the spiral progenitors more than the merger remnants themselves, an indication that these galaxies are not formed through merger events. We find, however, that recent mergers cannot be ruled out for some lenticular galaxies.

Here we present positions and radial velocities for over 4000 globular clusters (GCs) in 27 nearby early-type galaxies from the SLUGGS survey. The SLUGGS survey is designed to be representative of elliptical and lenticular galaxies in the stellar mass range 10 $<$ log M$_{\ast}$/M$_{\odot}$ $<$ 11.7. The data have been obtained over many years, mostly using the very stable multi-object spectrograph DEIMOS on the Keck II 10m telescope. Radial velocities are measured using the calcium triplet lines with a velocity accuracy of $\pm$ 10-15 km/s. We use phase space diagrams (i.e. velocity--position diagrams) to identify contaminants such as foreground stars and background galaxies, and to show that the contribution of GCs from neighboring galaxies is generally insignificant. Likely ultra-compact dwarfs are tabulated separately. We find that the mean velocity of the GC system is close to that of the host galaxy systemic velocity, indicating that the GC system is in overall dynamical equilibrium within the galaxy potential. We also find that the GC system velocity dispersion scales with host galaxy stellar mass in a similar manner to the Faber-Jackson relation for the stellar velocity dispersion. Publication of these GC radial velocity catalogs should enable further studies in many areas, such as GC system substructure, kinematics, and host galaxy mass measurements.

We study mass distributions within and beyond 5~effective radii ($R_{\rm e}$) in 23 early-type galaxies from the SLUGGS survey, using their globular cluster (GC) kinematic data. The data are obtained with Keck/DEIMOS spectrograph, and consist of line-of-sight velocities for ~$3500$ GCs, measured with a high precision of ~15 $\rm km\ s^{-1}$ per GC and extending out to $~13 R_{\rm e}$. We obtain the mass distribution in each galaxy using the tracer mass estimator of Watkins et al. and account for kinematic substructures, rotation of the GC systems and galaxy flattening in our mass estimates. The observed scatter between our mass estimates and results from the literature is less than 0.2 dex. The dark matter fraction within $5R_{\rm e}$ ($f_{\rm DM}$) increases from ~$0.6$ to ~$0.8$ for low- and high-mass galaxies, respectively, with some intermediate-mass galaxies ($M_*{\sim}10^{11}M_\odot$) having low $f_{\rm DM}\sim0.3$, which appears at odds with predictions from simple galaxy models. We show that these results are independent of the adopted orbital anisotropy, stellar mass-to-light ratio, and the assumed slope of the gravitational potential. However, the low $f_{\rm DM}$ in the ~$10^{11}M_\odot$ galaxies agrees with the cosmological simulations of Wu et al. where the pristine dark matter distribution has been modified by baryons during the galaxy assembly process. We find hints that these $M_*\sim10^{11}M_\odot$ galaxies with low $f_{\rm DM}$ have very diffuse dark matter haloes, implying that they assembled late. Beyond $5R_{\rm e}$, the $M/L$ gradients are steeper in the more massive galaxies and shallower in both low and intermediate mass galaxies.

We reconsider the origin and processing of dust in elliptical galaxies. We theoretically formulate the evolution of grain size distribution, taking into account dust supply from asymptotic giant branch (AGB) stars and dust destruction by sputtering in the hot interstellar medium (ISM), whose temperature evolution is treated by including two cooling paths: gas emission and dust emission (i.e. gas cooling and dust cooling). With our new full treatment of grain size distribution, we confirm that dust destruction by sputtering is too efficient to explain the observed dust abundance even if AGB stars continue to supply dust grains, and that, except for the case where the initial dust-to-gas ratio in the hot gas is as high as $\sim 0.01$, dust cooling is negligible compared with gas cooling. However, we show that, contrary to previous expectations, cooling does not help to protect the dust; rather, the sputtering efficiency is raised by the gas compression as a result of cooling. We additionally consider grain growth after the gas cools down. Dust growth by the accretion of gas-phase metals in the cold medium increase the dust-to-gas ratio up to $\sim 10^{-3}$ if this process lasts >10/(n_H/10^3 cm^{-3}) Myr, where $n_\mathrm{H}$ is the number density of hydrogen nuclei. We show that the accretion of gas-phase metals is a viable mechanism of increasing the dust abundance in elliptical galaxies to a level consistent with observations, and that the steepness of observed extinction curves is better explained with grain growth by accretion.

Motivated by the recent, serendipitous discovery of the densest known galaxy, M60-UCD1, we present two initial findings from a follow-up search, using the Sloan Digital Sky Survey, Subaru/Suprime-Cam and Hubble Space Telescope imaging, and SOuthern Astrophysical Research (SOAR)/Goodman spectroscopy. The first object discovered, M59-UCD3, has a similar size to M60-UCD1 (half-light radius of r_h ~ 20 pc) but is 40% more luminous (M_V ~ -14.6), making it the new densest-known galaxy. The second, M85-HCC1, has a size like a typical globular cluster (GC; r_h ~ 1.8 pc) but is much more luminous (M_V ~ -12.5). This hypercompact cluster is by far the densest confirmed free-floating stellar system, and is equivalent to the densest known nuclear star clusters. From spectroscopy, we find that both objects are relatively young (~9 Gyr and ~3 Gyr, respectively), with metal-abundances that resemble those of galaxy centers. Their host galaxies show clear signs of large-scale disturbances, and we conclude that these dense objects are the remnant nuclei of recently accreted galaxies. M59-UCD3 is an ideal target for follow-up with high-resolution imaging and spectroscopy to search for an overweight central supermassive black hole as was discovered in M60-UCD1. These findings also emphasize the potential value of ultra-compact dwarfs and massive GCs as tracers of the assembly histories of galaxies.

We recently identified a population of low surface brightness objects in the field of the z=0.023 Coma cluster, using the Dragonfly Telephoto Array. Here we present Keck spectroscopy of one of the largest of these "ultra-diffuse galaxies" (UDGs), confirming that it is a member of the cluster. The galaxy has prominent absorption features, including the Ca II H+K lines and the G-band, and no detected emission lines. Its radial velocity of cz=6280 +- 120 km/s is within the 1 sigma velocity dispersion of the Coma cluster. The galaxy has an effective radius of 4.3 +- 0.3 kpc and a Sersic index of 0.89 +- 0.06, as measured from Keck imaging. We find no indications of tidal tails or other distortions, at least out to a radius of ~2 r_e. We show that UDGs are located in a previously sparsely populated region of the size - magnitude plane of quiescent stellar systems, as they are ~6 magnitudes fainter than normal early-type galaxies of the same size. It appears that the luminosity distribution of large quiescent galaxies is not continuous, although this could largely be due to selection effects. Dynamical measurements are needed to determine whether the dark matter halos of UDGs are similar to those of galaxies with the same luminosity or to those of galaxies with the same size.

Stellar population synthesis (SPS) models are used to infer many galactic properties including star formation histories, metallicities, and stellar and dust masses. However, most SPS models neglect the effect of circumstellar dust shells around evolved stars and it is unclear to what extent they impact the analysis of SEDs. To overcome this shortcoming we have created a new set of circumstellar dust models, using the radiative transfer code DUSTY Ivezic et al. 1999, for asymptotic giant branch (AGB) stars and incorporated them into the Flexible Stellar Population Synthesis code. The circumstellar dust models provide a good fit to individual AGB stars as well as the IR color-magnitude diagrams of the Large and Small Magellanic Clouds. IR luminosity functions from the Large and Small Magellanic Clouds are not well-fit by the 2008 Padova isochrones when coupled to our circumstellar dust models, and so we adjusted the lifetimes of AGB stars in the models to provide a match to the data. We show, in agreement with previous work, that circumstellar dust from AGB stars can make a significant contribution to the IR ($\gtrsim4\mu m$) emission from galaxies that contain relatively little diffuse dust, including low-metallicity and/or non-star forming galaxies. Our models provide a good fit to the mid-IR spectra of early-type galaxies. Circumstellar dust around AGB stars appears to have a small effect on the IR SEDs of metal-rich star-forming galaxies (i.e., when A$_{\rm V}$ $\gtrsim$~0.1). Stellar population models that include circumstellar dust will be needed to accurately interpret data from the James Webb Space Telescope (JWST) and other IR facilities.

We present new wide-field photometry and spectroscopy of the globular clusters (GCs) around NGC 4649 (M60), the third brightest galaxy in the Virgo cluster. Imaging of NGC 4649 was assembled from a recently-obtained HST/ACS mosaic, and new Subaru/Suprime-Cam and archival CFHT/MegaCam data. About 1200 sources were followed up spectroscopically using combined observations from three multi-object spectrographs: Keck/DEIMOS, Gemini/GMOS and MMT/Hectospec. We confirm 431 unique GCs belonging to NGC 4649, a factor of 3.5 larger than previous datasets and with a factor of 3 improvement in velocity precision. We confirm significant GC colour bimodality and find that the red GCs are more centrally concentrated, while the blue GCs are more spatially extended. We infer negative GC colour gradients in the innermost 20 kpc and flat gradients out to large radii. Rotation is detected along the galaxy major axis for all tracers: blue GCs, red GCs, galaxy stars and planetary nebulae. We compare the observed properties of NGC 4649 with galaxy formation models. We find that formation via a major merger between two gas-poor galaxies, followed by satellite accretion, can consistently reproduce the observations of NGC 4649 at different radii. We find no strong evidence to support an interaction between NGC 4649 and the neighbouring spiral galaxy NGC 4647. We identify interesting GC kinematic features in our data, such as counter-rotating subgroups and bumpy kinematic profiles, which encode more clues about the formation history of NGC 4649.

We introduce and provide the scientific motivation for a wide-field photometric and spectroscopic chemodynamical survey of nearby early-type galaxies (ETGs) and their globular cluster (GC) systems. The SLUGGS (SAGES Legacy Unifying Globulars and GalaxieS) survey is being carried out primarily with Subaru/Suprime-Cam and Keck/DEIMOS. The former provides deep gri imaging over a 900 arcmin^2 field-of-view to characterize GC and host galaxy colors and spatial distributions, and to identify spectroscopic targets. The NIR CaII triplet provides GC line-of-sight velocities and metallicities out to typically ~8 R_e, and to ~15 R_e in some cases. New techniques to extract integrated stellar kinematics and metallicities to large radii (~2-3 R_e) are used in concert with GC data to create two-dimensional velocity and metallicity maps for comparison with simulations of galaxy formation. The advantages of SLUGGS compared with other, complementary, 2D-chemodynamical surveys are its superior velocity resolution, radial extent, and multiple halo tracers. We describe the sample of 25 nearby ETGs, the selection criteria for galaxies and GCs, the observing strategies, the data reduction techniques, and modeling methods. The survey observations are nearly complete and more than 30 papers have so far been published using SLUGGS data. Here we summarize some initial results, including signatures of two-phase galaxy assembly, evidence for GC metallicity bimodality, and a novel framework for the formation of extended star clusters and ultracompact dwarfs. An integrated overview of current chemodynamical constraints on GC systems points to separate, in-situ formation modes at high redshifts for metal-poor and metal-rich GCs.